Background

Solid organ transplant recipients have heightened risk for diffuse large B cell lymphoma (DLBCL). The role of donor-recipient HLA mismatch and recipient HLA type on DLBCL risk are not well established.

Methods

We examined 172 231 kidney, heart, pancreas, and lung recipients transplanted in the United States between 1987 and 2010, including 902 with DLBCL. Incidence rate ratios (IRRs) were calculated using Poisson regression for DLBCL risk in relation to HLA mismatch, types, and zygosity, adjusting for sex, age, race/ethnicity, year, organ, and transplant number.

Results

Compared with recipients who had 2 HLA-DR mismatches, those with zero or 1 mismatch had reduced DLBCL risk, (zero: IRR, 0.76, 95% confidence interval [95% CI], 0.61-0.95; one: IRR, 0.83; 95% CI, 0.69-1.00). In stratified analyses, recipients matched at either HLA-A, -B, or -DR had a significantly reduced risk of late-onset (>2 years after transplantation), but not early-onset DLBCL, and there was a trend for decreasing risk with decreasing mismatch across all 3 loci (P = 0.0003). Several individual recipient HLA-A, -B, -C, -DR, and -DQ antigens were also associated with DLBCL risk, including DR13 (IRR, 0.74; 95% CI, 0.57-0.93) and B38 (IRR, 1.48; 95% CI, 1.10-1.93), confirming prior findings that these 2 antigens are associated with risk of infection-associated cancers.

Conclusions

In conclusion, variation in HLA is related to susceptibility to DLBCL, perhaps reflecting intensity of immunosuppression, control of Epstein-Barr virus infection among transplant recipients or chronic immune stimulation.

Background

Epidemiological studies have demonstrated a positive association between chronic lymphocytic leukaemia (CLL) and non-melanoma skin cancer (NMSC). We hypothesized that shared genetic risk factors between CLL and NMSC could contribute to the association observed between these diseases.

Methods

We examined the association between (i) established NMSC susceptibility loci and CLL risk in a meta-analysis including 3100 CLL cases and 7667 controls and (ii) established CLL loci and NMSC risk in a study of 4242 basal cell carcinoma (BCC) cases, 825 squamous cell carcinoma (SCC) cases and 12802 controls. Polygenic risk scores (PRS) for CLL, BCC and SCC were constructed using established loci. Logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs).

Results

Higher CLL-PRS was associated with increased BCC risk (OR4th-quartile-vs-1st-quartile = 1.13, 95% CI: 1.02-1.24, Ptrend = 0.009), even after removing the shared 6p25.3 locus. No association was observed with BCC-PRS and CLL risk (Ptrend = 0.68). These findings support a contributory role for CLL in BCC risk, but not for BCC in CLL risk. Increased CLL risk was observed with higher SCC-PRS (OR4th-quartile-vs-1st-quartile = 1.22, 95% CI: 1.08-1.38, Ptrend = 1.36 × 10-5), which was driven by shared genetic susceptibility at the 6p25.3 locus.

Conclusion

These findings highlight the role of pleiotropy regarding the pathogenesis of CLL and NMSC and shows that a single pleiotropic locus, 6p25.3, drives the observed association between genetic susceptibility to SCC and increased CLL risk. The study also provides evidence that genetic susceptibility for CLL increases BCC risk.

The use of publicly available sequencing datasets as controls (hereafter, "public controls") in studies of rare variant disease associations has great promise but can increase the risk of false-positive discovery. The specific factors that could contribute to inflated distribution of test statistics have not been systematically examined. Here, we leveraged both public controls, gnomAD v2.1 and several datasets sequenced in our laboratory to systematically investigate factors that could contribute to the false-positive discovery, as measured by λΔ95, a measure to quantify the degree of inflation in statistical significance. Analyses of datasets in this investigation found that 1) the significantly inflated distribution of test statistics decreased substantially when the same variant caller and filtering pipelines were employed, 2) differences in library prep kits and sequencers did not affect the false-positive discovery rate and, 3) joint vs. separate variant-calling of cases and controls did not contribute to the inflation of test statistics. Currently available methods do not adequately adjust for the high false-positive discovery. These results, especially if replicated, emphasize the risks of using public controls for rare-variant association tests in which individual-level data and the computational pipeline are not readily accessible, which prevents the use of the same variant-calling and filtering pipelines on both cases and controls. A plausible solution exists with the emergence of cloud-based computing, which can make it possible to bring containerized analytical pipelines to the data (rather than the data to the pipeline) and could avert or minimize these issues. It is suggested that future reports account for this issue and provide this as a limitation in reporting new findings based on studies that cannot practically analyze all data on a single pipeline.

Abstract

Compared with the general population, hematopoietic cell transplantation (HCT) survivors are at elevated risk for developing solid subsequent neoplasms (SNs). The Center for International Blood and Marrow Transplant Research (CIBMTR) is a key resource for quantifying solid SN incidence following HCT, but the completeness of SN ascertainment is uncertain. Within a cohort of 18 450 CIBMTR patients linked to the California Cancer Registry (CCR), we evaluated the completeness of solid SN data reported to the CIBMTR from 1991 to 2018 to understand the implications of using CIBMTR data alone or combined with CCR data to quantify the burden of solid SNs after HCT. We estimated the cumulative incidence of developing a solid SN, accounting for the competing risk of death. Within the cohort, solid SNs were reported among 724 patients; 15.6% of these patients had an SN reported by CIBMTR only, 36.9% by CCR only, and 47.5% by both. The corresponding cumulative incidence of developing a solid SN at 10 years following a first HCT was 4.0% (95% confidence interval [CI], 3.5-4.4) according to CIBMTR data only, 5.3% (95% CI, 4.9-5.9) according to CCR data only, and 6.3% (95% CI, 5.7-6.8) according to both sources combined. The patterns were similar for allogeneic and autologous HCT recipients. Linking detailed HCT information from CIBMTR with comprehensive SN data from cancer registries provides an opportunity to optimize SN ascertainment for informing follow-up care practices and evaluating risk factors in the growing population of HCT survivors.

Non-Hodgkin lymphoma (NHL) represents a diverse group of hematological malignancies, of which follicular lymphoma (FL) is a prevalent subtype. A previous genome-wide association study has established a marker, rs10484561 in the human leukocyte antigen (HLA) class II region on 6p21.32 associated with increased FL risk. Here, in a three-stage genome-wide association study, starting with a genome-wide scan of 379 FL cases and 791 controls followed by validation in 1,049 cases and 5,790 controls, we identified a second independent FL-associated locus on 6p21.32, rs2647012 (OR(combined)  = 0.64, P(combined)  = 2 × 10(-21)) located 962 bp away from rs10484561 (r(2)<0.1 in controls). After mutual adjustment, the associations at the two SNPs remained genome-wide significant (rs2647012:OR(adjusted)  = 0.70, P(adjusted)  =  4 × 10(-12); rs10484561:OR(adjusted)  = 1.64, P(adjusted)  = 5 × 10(-15)). Haplotype and coalescence analyses indicated that rs2647012 arose on an evolutionarily distinct haplotype from that of rs10484561 and tags a novel allele with an opposite (protective) effect on FL risk. Moreover, in a follow-up analysis of the top 6 FL-associated SNPs in 4,449 cases of other NHL subtypes, rs10484561 was associated with risk of diffuse large B-cell lymphoma (OR(combined)  = 1.36, P(combined)  =  1.4 × 10(-7)). Our results reveal the presence of allelic heterogeneity within the HLA class II region influencing FL susceptibility and indicate a possible shared genetic etiology with diffuse large B-cell lymphoma. These findings suggest that the HLA class II region plays a complex yet important role in NHL.

Autoimmune conditions and immune system-related genetic variations are associated with risk of non-Hodgkin lymphoma (NHL). In a pooled analysis of 8,692 NHL cases and 9,260 controls from 14 studies (1988-2007) within the International Lymphoma Epidemiology Consortium, we evaluated the interaction between immune system genetic variants and autoimmune conditions in NHL risk. We evaluated the immunity-related single nucleotide polymorphisms rs1800629 (tumor necrosis factor gene (TNF) G308A), rs1800890 (interleukin-10 gene (IL10) T3575A), rs6457327 (human leukocyte antigen gene (HLA) class I), rs10484561 (HLA class II), and rs2647012 (HLA class II)) and categorized autoimmune conditions as primarily mediated by B-cell or T-cell responses. We constructed unconditional logistic regression models to measure associations between autoimmune conditions and NHL with stratification by genotype. Autoimmune conditions mediated by B-cell responses were associated with increased NHL risk, specifically diffuse large B-cell lymphoma (odds ratio (OR) = 3.11, 95% confidence interval (CI): 2.25, 4.30) and marginal zone lymphoma (OR = 5.80, 95% CI: 3.82, 8.80); those mediated by T-cell responses were associated with peripheral T-cell lymphoma (OR = 2.14, 95% CI: 1.35, 3.38). In the presence of the rs1800629 AG/AA genotype, B-cell-mediated autoimmune conditions increased NHL risk (OR = 3.27, 95% CI: 2.07, 5.16; P-interaction = 0.03) in comparison with the GG genotype (OR = 1.82, 95% CI: 1.31, 2.53). This interaction was consistent across major B-cell NHL subtypes, including marginal zone lymphoma (P-interaction = 0.02) and follicular lymphoma (P-interaction = 0.04).

Marginal zone lymphoma (MZL) is the third most common subtype of B-cell non-Hodgkin lymphoma. Here we perform a two-stage GWAS of 1,281 MZL cases and 7,127 controls of European ancestry and identify two independent loci near BTNL2 (rs9461741, P=3.95 × 10(-15)) and HLA-B (rs2922994, P=2.43 × 10(-9)) in the HLA region significantly associated with MZL risk. This is the first evidence that genetic variation in the major histocompatibility complex influences MZL susceptibility.

To identify susceptibility loci for non-Hodgkin lymphoma subtypes, we conducted a three-stage genome-wide association study. We identified two variants associated with follicular lymphoma at 6p21.32 (rs10484561, combined P = 1.12 x 10(-29) and rs7755224, combined P = 2.00 x 10(-19); r(2) = 1.0), supporting the idea that major histocompatibility complex genetic variation influences follicular lymphoma susceptibility. We also found confirmatory evidence of a previously reported association between chronic lymphocytic leukemia/small lymphocytic lymphoma and rs735665 (combined P = 4.24 x 10(-9)).

Genome-wide association studies (GWASs) of follicular lymphoma (FL) have previously identified human leukocyte antigen (HLA) gene variants. To identify additional FL susceptibility loci, we conducted a large-scale two-stage GWAS in 4,523 case subjects and 13,344 control subjects of European ancestry. Five non-HLA loci were associated with FL risk: 11q23.3 (rs4938573, p = 5.79 × 10(-20)) near CXCR5; 11q24.3 (rs4937362, p = 6.76 × 10(-11)) near ETS1; 3q28 (rs6444305, p = 1.10 × 10(-10)) in LPP; 18q21.33 (rs17749561, p = 8.28 × 10(-10)) near BCL2; and 8q24.21 (rs13254990, p = 1.06 × 10(-8)) near PVT1. In an analysis of the HLA region, we identified four linked HLA-DRβ1 multiallelic amino acids at positions 11, 13, 28, and 30 that were associated with FL risk (pomnibus = 4.20 × 10(-67) to 2.67 × 10(-70)). Additional independent signals included rs17203612 in HLA class II (odds ratio [OR(per-allele)] = 1.44; p = 4.59 × 10(-16)) and rs3130437 in HLA class I (OR(per-allele) = 1.23; p = 8.23 × 10(-9)). Our findings further expand the number of loci associated with FL and provide evidence that multiple common variants outside the HLA region make a significant contribution to FL risk.

Genome-wide association studies (GWAS) have previously identified 13 loci associated with risk of chronic lymphocytic leukemia or small lymphocytic lymphoma (CLL). To identify additional CLL susceptibility loci, we conducted the largest meta-analysis for CLL thus far, including four GWAS with a total of 3,100 individuals with CLL (cases) and 7,667 controls. In the meta-analysis, we identified ten independent associated SNPs in nine new loci at 10q23.31 (ACTA2 or FAS (ACTA2/FAS), P=1.22×10(-14)), 18q21.33 (BCL2, P=7.76×10(-11)), 11p15.5 (C11orf21, P=2.15×10(-10)), 4q25 (LEF1, P=4.24×10(-10)), 2q33.1 (CASP10 or CASP8 (CASP10/CASP8), P=2.50×10(-9)), 9p21.3 (CDKN2B-AS1, P=1.27×10(-8)), 18q21.32 (PMAIP1, P=2.51×10(-8)), 15q15.1 (BMF, P=2.71×10(-10)) and 2p22.2 (QPCT, P=1.68×10(-8)), as well as an independent signal at an established locus (2q13, ACOXL, P=2.08×10(-18)). We also found evidence for two additional promising loci below genome-wide significance at 8q22.3 (ODF1, P=5.40×10(-8)) and 5p15.33 (TERT, P=1.92×10(-7)). Although further studies are required, the proximity of several of these loci to genes involved in apoptosis suggests a plausible underlying biological mechanism.