Abstract
A previous International Lymphoma Epidemiology (InterLymph) Consortium evaluation of joint associations between five immune gene variants and autoimmune conditions reported interactions between B-cell response-mediated autoimmune conditions and the rs1800629 genotype on risk of B-cell non–Hodgkin lymphoma (NHL) subtypes. Here, we extend that evaluation using NHL subtype-specific polygenic risk scores (PRS) constructed from loci identified in genome-wide association studies of three common B-cell NHL subtypes.
In a pooled analysis of NHL cases and controls of Caucasian descent from 14 participating InterLymph studies, we evaluated joint associations between B-cell–mediated autoimmune conditions and tertile (T) of PRS for risk of diffuse large B-cell lymphoma (DLBCL; n = 1,914), follicular lymphoma (n = 1,733), and marginal zone lymphoma (MZL; n = 407), using unconditional logistic regression.
We demonstrated a positive association of DLBCL PRS with DLBCL risk [T2 vs. T1: OR = 1.24; 95% confidence interval (CI), 1.08–1.43; T3 vs. T1: OR = 1.81; 95% CI, 1.59–2.07; P-trend (Ptrend) < 0.0001]. DLBCL risk also increased with increasing PRS tertile among those with an autoimmune condition, being highest for those with a B-cell–mediated autoimmune condition and a T3 PRS [OR = 6.46 vs. no autoimmune condition and a T1 PRS, Ptrend < 0.0001, P-interaction (Pinteraction) = 0.49]. Follicular lymphoma and MZL risk demonstrated no evidence of joint associations or significant Pinteraction.
Our results suggest that PRS constructed from currently known subtype-specific loci may not necessarily capture biological pathways shared with autoimmune conditions.
Targeted genetic (PRS) screening among population subsets with autoimmune conditions may offer opportunities for identifying those at highest risk for (and early detection from) DLBCL.
Introduction
Since its inception in 2001, the International Lymphoma Epidemiology (InterLymph) Consortium has sought to identify genetic and nongenetic risk factors for non–Hodgkin lymphoma (NHL) and its heterogenous subtypes by conducting pooled analyses comprising thousands of patients and controls from case–control epidemiologic studies across multiple countries (https://epi.grants.cancer.gov/interlymph). Consortium members have demonstrated statistical associations between autoimmune conditions and NHL subtypes (1–3), most notably between B-cell–mediated autoimmune conditions and diffuse large B-cell lymphoma (DLBCL) and marginal zone lymphoma (MZL; refs. 4, 5). Although autoimmune conditions are generally rare, their individual associations with NHL subtypes are among the most robust and strongest (e.g., over 2-fold) in NHL etiology (1). In contrast, genome-wide association studies (GWAS) have identified statistically significant associations with NHL subtypes for multiple common genetic susceptibility loci, but their associations have largely exhibited modest magnitudes of risk (e.g., less than 2-fold; refs. 6–11), as is typical for genetic associations studies.
In a previous InterLymph Consortium evaluation of joint associations between autoimmune conditions and five putative susceptibility loci among genes potentially linked with immune function, we demonstrated a statistically significant interaction between TNF -308G-A (rs1800629) and B-cell–mediated autoimmune conditions for DLBCL risk (3). Since that publication, several other NHL GWASs have yielded multiple highly significant susceptibility loci; notably, the growing list of implicated putative loci differs by NHL subtype. Here, we construct polygenic risk scores (PRS) derived from NHL subtype-specific GWAS results for three common NHL subtypes, including DLBCL, follicular lymphoma, and MZL, and evaluate the joint associations between these subtype-specific PRS and autoimmune conditions to determine whether autoimmune conditions and genetic susceptibility to these NHL subtypes share common biological underpinnings.
Materials and Methods
Study population
The present analysis included data from 14 case–control studies that previously participated in an InterLymph Consortium (https://epi.grants.cancer.gov/interlymph/) pooled analysis of autoimmune conditions in relation to NHL risk (Supplementary Table S1; refs. 1–3; 12–25). As previously described, NHL cases and controls from the 14 studies were eligible if they were: age 17 years or older, not known to be human immunodeficiency virus (HIV)-positive, and had no history of organ transplantation (3). Among the eligible participants, inclusion further required genotype data from NHL GWAS studies published from 2014 to 15 (6–8). To minimize confounding by race, the GWAS studies were conducted among Caucasians. Specifically, ancestry was assessed using the Genotyping Library and Utilities (GLU- http://code.google.com/p/glu-genetics/) struct.admix module based on the method by Pritchard and colleagues (26) and participants with less than 80% European ancestry were excluded (6). The final study sample for the present analysis thus reflects a subset of those from previously published studies who also had DNA and passed quality control for GWAS analyses. A comparison of select demographic characteristics of participants included in this analysis and those in our previous publication is summarized in Supplementary Table S2 (3). The final analytic dataset comprised 5,886 cases and 5,687 controls.
Exposure assessment
As previously described, harmonization of autoimmune conditions was conducted previously as part of a consortium-wide effort and detailed in the resulting publications (1, 3). Briefly, self-reported history of autoimmune conditions was collected in each participating study using structured questionnaires during in-person or telephone interviews. In most studies (70%), respondents were asked whether any autoimmune condition had been diagnosed by a physician. Consistent with the original InterLymph Consortium study on autoimmune conditions (1), we included the following: primary Sjögren syndrome, systemic lupus erythematosus (SLE), rheumatoid arthritis, systemic sclerosis or scleroderma, poly- or dermatomyositis, immune thrombocytopenic purpura, type-1 diabetes (defined as diabetes diagnosed at age ≤30 years), pernicious anemia, multiple sclerosis, myasthenia gravis, celiac disease, psoriasis, sarcoidosis, Crohn disease, ulcerative colitis, autoimmune hemolytic anemia, and Hashimoto thyroiditis. No imputation was conducted for missing data; studies and individuals with missing data were excluded from each of the respective analyses. Based on this harmonization, autoimmune conditions were then categorized on the basis of the type of primary immune response involved in mediating autoimmunity: specifically, predominance of B-cell activation versus predominance of T-cell activation, based on a consensus panel comprised of rheumatologists, immunologists, and hematologist-oncologists (3). Autoimmune conditions were also categorized by organ involvement as multiple-organ–targeted versus primarily single-organ–targeted, with further organ-specific evaluations for pancreatic, gastrointestinal/hepatobiliary, dermatologic, hematologic, neurologic, and endocrine organs. B-cell–mediated autoimmune conditions included autoimmune hemolytic anemia, Hashimoto thyroiditis, primary Sjögren syndrome, SLE, rheumatoid arthritis, pernicious anemia, and myasthenia gravis. T-cell–mediated autoimmune conditions included systemic sclerosis or scleroderma, poly- or dermatomyositis, immune thrombocytopenic purpura, type-1 diabetes, multiple sclerosis, celiac disease, psoriasis, sarcoidosis, Crohn disease, and ulcerative colitis. These categorizations were based on the type of primary immune response believed to be involved in mediating autoimmunity. Categorizations by organ involvement (multiple-organ– versus primarily single-organ–targeted) were also evaluated (Supplementary Table S3; ref. 3). Because the focus of the present manuscript is on the three major B-cell lymphoma subtypes, our analyses and results are focused on the B-cell–mediated autoimmune conditions for which prior associations were reported (3).
PRS
GWAS data were ascertained through a number of efforts and platforms, including Illumina 317K [Scandinavian Lymphoma Etiology (SCALE) study], Illumina Human CNV370-Duo BeadChip [University of California, San Francisco (UCSF), San Francisco, CA), OmniExpress [NCI–Surveillance Epidemiology and End Results (SEER); New South Wales, Australia; Yale, New Haven, CN; British Columbia, Canada; UCSF, United Kingdom], Illumina660 (Mayo Clinic), Sequenom MassARRAY iPLEX (SF1B), Illumina GoldenGate, and Pyrosequencing (EpiLymph). Collectively, we included loci for which call rates were ≥95% and sample completion rates were ≥90% (6–9). To evaluate the genome across studies, we thus imputed all GWAS data using the 1000 Genomes Project (1kGP) v.3 (March 2012 release, http://www.1000genomes.org/) reference panel and IMPUTE2 (http://mathgen.stats.ox.ac.uk/impute/impute_v2.html). The imputation analysis was restricted to common SNPs (cut-off minor allele frequency > 0.01 with imputation accuracy INFO score >0.3). Wherever possible, either the best directly genotyped SNP or the most significant imputed SNP (with information score > 0.8) for the locus was used. We note that genotyping of SNPs by TaqMan (Applied Biosystems) in a subset of subjects yielded more than 88.9% concordance with the imputed dosages (median concordance = 99.6%), indicating that imputation accuracy was high (4, 7, 8).
SNPs included for each PRS comprised those statistically significantly associated with the respective NHL subtype in recently published GWAS that included InterLymph studies. Specifically, we included the following SNPs and corresponding reported ORs to derive subtype-specific PRS: DLBCL [rs116446171 (OR = 2.2), rs2523607 (OR = 1.32), rs79480871 (OR = 1.34); rs13255292 (OR = 1.22); rs4733601 (OR = 1.18)] (6); follicular lymphoma [rs4938573 (OR = 1.34), rs4937362 (OR = 1.19), rs6444305 (OR = 1.21), rs17749561 (OR = 1.34), rs13254990 (OR = 1.18), rs12195582/rs115374828 (OR = 1.44)] (7); and MZL [rs9461741 (OR = 2.66), rs2922994/rs116778584 (OR = 1.64)] (8). For each SNP, there was no statistically significant heterogeneity by study (6–8).
To construct the PRS, homozygous risk alleles were assigned a value of 2, heterozygotes a value of 1, and homozygous nonrisk alleles a value of 0 (27, 28). This value was subsequently multiplied by the log of the published OR, and the resulting values were summed across all SNPs included in each subtype-specific PRS. Subtype-specific PRSs were then categorised into tertiles using the PRS distribution among controls for each subtype. We note that numbers of controls vary for each NHL subtype due to the different numbers of participants with complete data for each of the SNPs required to construct the given PRS.
NHL classification
NHL subtypes were grouped as per InterLymph Pathology Working Group guidelines (29–31). Results are presented for two NHL subtypes for which associations with several autoimmune conditions have been consistently reported: DLBCL and MZL. We also included follicular lymphoma because we previously observed an excess risk of follicular lymphoma among those with a history of Sjogren syndrome (32) and an association with autoimmune conditions among those with a variant rs1800629 allele (3). Our analytic dataset included 1914 DLBCLs, 1733 follicular lymphomas, and 407 MZLs.
Statistical methods
Independent associations
We first confirmed associations between autoimmune conditions and NHL subtypes in our subset of eligible cases and controls. Unconditional logistic regression, adjusted for age as a continuous variable, sex, and geographical region/study center (Supplementary Table S2), were used to calculate ORs and 95% confidence intervals (CI) as estimates of NHL subtype risk. ORs and 95% CIs for each NHL subtype were similarly calculated for each subtype-specific PRS, comparing tertiles (T) 2 and T3 to T1. We calculated P-trend (Ptrend) using the Wilcoxon rank-sum test across the ordered categories to measure whether magnitudes of risk across PRS tertiles T2 and T3 yielded a significant trend. For both independent models (e.g., models examining only autoimmune conditions or PRS), we adjusted for age. Potential confounders, such as socioeconomic status, smoking status, and family history of hematologic malignancies did not change risk estimates more than 10% and were thus not retained in any of the models. Heterogeneity was evaluated using χ2 tests of interaction between the studies and variables of interest; as there was no strong evidence of heterogeneity by study or any other variable of interest, all data were examined in a logistic regression model (by subtype).
Stratified associations
We tested for associations between PRS and autoimmune conditions on an additive scale for each subtype (DLBCL, follicular lymphoma, MZL) by calculating ORs and 95% CIs relative to those without autoimmune conditions and a PRS in T1 as the common referent group (Tables 1,2–3). We calculated Ptrends (Wilcoxon rank-sum test) across ordered categories to test for trends. Potential interaction on the multiplicative scale was assessed by stratifying on PRS tertile within each category of autoimmune condition. The P value for interaction was estimated using the Wald test for homogeneity of the associations of autoimmune conditions with NHL subtype risk by PRS tertile. Parallel analyses by individual SNPs were also conducted (Supplementary Tables S4–S6).
Autoimmune conditions . | DLBCL PRS . | Controls . | DLBCL . | OR (95% CI) . | Pinteractiona . |
---|---|---|---|---|---|
None | 5,064 | 1,808 | 1.00 (Ref) | ||
Any | 193 | 106 | 1.42 (1.10–1.83) | – | |
No B-cell conditions | 5,046 | 1,796 | 1.00 (Ref) | ||
B-cell conditions | 53 | 50 | 2.90 (1.94–4.35) | – | |
No multiple organs | 5,053 | 1,797 | 1.00 (Ref) | ||
Multiple organs | 64 | 46 | 2.24 (1.50–3.33) | – | |
None/any | T1 | 1,965 | 549 | 1.00 (Ref) | |
T2 | 1,717 | 578 | 1.24 (1.08–1.43) | ||
T3 | 1,575 | 787 | 1.81 (1.59–2.07) | – | |
None | T1 | 1,888 | 513 | 1.00 (Ref) | |
T2 | 1,656 | 550 | 1.25 (1.09–1.44) | ||
T3 | 1,520 | 745 | 1.82 (1.59–2.08) | ||
Any | T1 | 77 | 36 | 1.52 (0.99–2.33) | |
T2 | 61 | 28 | 1.73 (1.07–2.79) | Pinteraction | |
T3 | 55 | 42 | 2.65 (1.71–4.10) | = 0.82 | |
No B-cell conditions | T1 | 1,881 | 512 | 1.00 (Ref) | |
T2 | 1,647 | 545 | 1.24 (1.08–1.43) | ||
T3 | 1,518 | 739 | 1.80 (1.57–2.07) | ||
B-cell conditions | T1 | 22 | 18 | 2.94 (1.52–5.70) | |
T2 | 18 | 14 | 3.27 (1.56–6.84) | Pinteraction | |
T3 | 13 | 18 | 6.46 (3.04–13.7) | = 0.49 | |
No multiple organs | T1 | 1,885 | 512 | 1.00 (ref) | |
T2 | 1,650 | 544 | 1.24 (1.08–1.43) | ||
T3 | 1,518 | 741 | 1.81 (1.58–2.07) | ||
Multiple organs | T1 | 31 | 17 | 2.18 (1.16–4.11) | |
T2 | 19 | 12 | 2.54 (1.17–5.48) | Pinteraction | |
T3 | 14 | 17 | 5.92 (2.80–12.5) | = 0.34 |
Autoimmune conditions . | DLBCL PRS . | Controls . | DLBCL . | OR (95% CI) . | Pinteractiona . |
---|---|---|---|---|---|
None | 5,064 | 1,808 | 1.00 (Ref) | ||
Any | 193 | 106 | 1.42 (1.10–1.83) | – | |
No B-cell conditions | 5,046 | 1,796 | 1.00 (Ref) | ||
B-cell conditions | 53 | 50 | 2.90 (1.94–4.35) | – | |
No multiple organs | 5,053 | 1,797 | 1.00 (Ref) | ||
Multiple organs | 64 | 46 | 2.24 (1.50–3.33) | – | |
None/any | T1 | 1,965 | 549 | 1.00 (Ref) | |
T2 | 1,717 | 578 | 1.24 (1.08–1.43) | ||
T3 | 1,575 | 787 | 1.81 (1.59–2.07) | – | |
None | T1 | 1,888 | 513 | 1.00 (Ref) | |
T2 | 1,656 | 550 | 1.25 (1.09–1.44) | ||
T3 | 1,520 | 745 | 1.82 (1.59–2.08) | ||
Any | T1 | 77 | 36 | 1.52 (0.99–2.33) | |
T2 | 61 | 28 | 1.73 (1.07–2.79) | Pinteraction | |
T3 | 55 | 42 | 2.65 (1.71–4.10) | = 0.82 | |
No B-cell conditions | T1 | 1,881 | 512 | 1.00 (Ref) | |
T2 | 1,647 | 545 | 1.24 (1.08–1.43) | ||
T3 | 1,518 | 739 | 1.80 (1.57–2.07) | ||
B-cell conditions | T1 | 22 | 18 | 2.94 (1.52–5.70) | |
T2 | 18 | 14 | 3.27 (1.56–6.84) | Pinteraction | |
T3 | 13 | 18 | 6.46 (3.04–13.7) | = 0.49 | |
No multiple organs | T1 | 1,885 | 512 | 1.00 (ref) | |
T2 | 1,650 | 544 | 1.24 (1.08–1.43) | ||
T3 | 1,518 | 741 | 1.81 (1.58–2.07) | ||
Multiple organs | T1 | 31 | 17 | 2.18 (1.16–4.11) | |
T2 | 19 | 12 | 2.54 (1.17–5.48) | Pinteraction | |
T3 | 14 | 17 | 5.92 (2.80–12.5) | = 0.34 |
Note: See Supplementary Table S3 for definition of autoimmune condition categories.
aWald test for homogeneity.
Autoimmune conditions . | FL PRS . | Controls . | FL . | OR (95% CI) . | Pinteractiona . |
---|---|---|---|---|---|
None | 5,142 | 1,666 | 1.00 (Ref) | ||
Any | 196 | 67 | 0.91 (0.67–1.23) | – | |
No B-cell conditions | 5,123 | 1,660 | 1.00 (Ref) | ||
B-cell conditions | 51 | 19 | 1.12 (0.65–1.93) | – | |
No multiple organs | 5,131 | 1,661 | 1.00 (Ref) | ||
Multiple organs | 65 | 18 | 1.10 (0.66–1.82) | – | |
None/any | T1 | 2,367 | 298 | 1.00 (Ref) | |
T2 | 1,680 | 549 | 2.38 (2.03–2.80) | ||
T3 | 1,127 | 832 | 4.89 (4.18–5.73) | – | |
None | T1 | 2,352 | 292 | 1.00 (Ref) | |
T2 | 1,675 | 543 | 2.37 (2.01–2.79) | ||
T3 | 1,115 | 829 | 4.99 (4.24–5.86) | ||
Any | T1 | 87 | 12 | 1.01 (1.68–4.84) | |
T2 | 55 | 24 | 2.85 (1.97–5.22) | Pinteraction | |
T3 | 54 | 31 | 3.21 (1.07–1.37) | = 0.30 | |
No B-cell conditions | T1 | 2,345 | 292 | 1.00 (Ref) | |
T2 | 1,666 | 543 | 2.37 (2.01–2.79) | ||
T3 | 1,112 | 825 | 4.98 (4.23–5.85) | ||
B-cell conditions | T1 | 22 | 6 | 1.96 (0.77–5.02) | |
T2 | 14 | 6 | 3.61 (1.33–9.76) | Pinteraction | |
T3 | 15 | 7 | 2.92 (1.15–7.39) | = 0.10 | |
No multiple organs | T1 | 2,346 | 292 | 1.00 (Ref) | |
T2 | 1,672 | 543 | 2.36 (2.01–2.78) | ||
T3 | 1,113 | 826 | 4.98 (4.24–5.86) | ||
Multiple organs | T1 | 30 | 8 | 2.02 (0.89–4.58) | |
T2 | 18 | 10 | 4.15 (1.79–9.61) | Pinteraction | |
T3 | 17 | 5 | 2.11 (0.75–5.94) | = 0.20 | |
No organ-specific conditions | T1 | 2,349 | 291 | 1.00 (Ref) | |
T2 | 1,668 | 545 | 2.39 (2.03–2.81) | ||
T3 | 1,113 | 825 | 4.99 (4.24–5.97) | ||
Organ-specific conditions | T1 | 59 | 6 | 0.73 (0.31–1.72) | |
T2 | 38 | 15 | 2.44 (1.28–4.66) | Pinteraction | |
T3 | 38 | 26 | 3.57 (2.07–6.17) | = 0.78 |
Autoimmune conditions . | FL PRS . | Controls . | FL . | OR (95% CI) . | Pinteractiona . |
---|---|---|---|---|---|
None | 5,142 | 1,666 | 1.00 (Ref) | ||
Any | 196 | 67 | 0.91 (0.67–1.23) | – | |
No B-cell conditions | 5,123 | 1,660 | 1.00 (Ref) | ||
B-cell conditions | 51 | 19 | 1.12 (0.65–1.93) | – | |
No multiple organs | 5,131 | 1,661 | 1.00 (Ref) | ||
Multiple organs | 65 | 18 | 1.10 (0.66–1.82) | – | |
None/any | T1 | 2,367 | 298 | 1.00 (Ref) | |
T2 | 1,680 | 549 | 2.38 (2.03–2.80) | ||
T3 | 1,127 | 832 | 4.89 (4.18–5.73) | – | |
None | T1 | 2,352 | 292 | 1.00 (Ref) | |
T2 | 1,675 | 543 | 2.37 (2.01–2.79) | ||
T3 | 1,115 | 829 | 4.99 (4.24–5.86) | ||
Any | T1 | 87 | 12 | 1.01 (1.68–4.84) | |
T2 | 55 | 24 | 2.85 (1.97–5.22) | Pinteraction | |
T3 | 54 | 31 | 3.21 (1.07–1.37) | = 0.30 | |
No B-cell conditions | T1 | 2,345 | 292 | 1.00 (Ref) | |
T2 | 1,666 | 543 | 2.37 (2.01–2.79) | ||
T3 | 1,112 | 825 | 4.98 (4.23–5.85) | ||
B-cell conditions | T1 | 22 | 6 | 1.96 (0.77–5.02) | |
T2 | 14 | 6 | 3.61 (1.33–9.76) | Pinteraction | |
T3 | 15 | 7 | 2.92 (1.15–7.39) | = 0.10 | |
No multiple organs | T1 | 2,346 | 292 | 1.00 (Ref) | |
T2 | 1,672 | 543 | 2.36 (2.01–2.78) | ||
T3 | 1,113 | 826 | 4.98 (4.24–5.86) | ||
Multiple organs | T1 | 30 | 8 | 2.02 (0.89–4.58) | |
T2 | 18 | 10 | 4.15 (1.79–9.61) | Pinteraction | |
T3 | 17 | 5 | 2.11 (0.75–5.94) | = 0.20 | |
No organ-specific conditions | T1 | 2,349 | 291 | 1.00 (Ref) | |
T2 | 1,668 | 545 | 2.39 (2.03–2.81) | ||
T3 | 1,113 | 825 | 4.99 (4.24–5.97) | ||
Organ-specific conditions | T1 | 59 | 6 | 0.73 (0.31–1.72) | |
T2 | 38 | 15 | 2.44 (1.28–4.66) | Pinteraction | |
T3 | 38 | 26 | 3.57 (2.07–6.17) | = 0.78 |
aWald test for homogeneity.
Autoimmune conditions . | MZL PRS . | Controls . | MZL . | OR (95% CI) . | Pinteractiona . |
---|---|---|---|---|---|
None | 2,774 | 376 | 1.00 (ref) | ||
Any | 109 | 31 | 1.97 (1.33–2.94) | – | |
No B-cell conditions | 2,764 | 373 | 1.00 (ref) | ||
B-cell conditions | 28 | 24 | 5.88 (3.55–9.74) | – | |
No multiple organs | 2,770 | 374 | 1.00 (ref) | ||
Multiple organs | 33 | 24 | 4.96 (3.04–8.12) | – | |
None/any | T1 | 2,123 | 294 | 1.00 (ref) | |
T2 | 586 | 86 | 1.10 (0.84–1.43) | ||
T3 | 174 | 27 | 1.09 (0.71–1.68) | – | |
None | T1 | 2,046 | 272 | 1.00 (ref) | |
T2 | 567 | 80 | 1.10 (0.83–1.44) | ||
T3 | 161 | 24 | 1.10 (0.70–1.74) | ||
Any | T1 | 77 | 22 | 2.13 (1.28–3.55) | |
T2 | 19 | 6 | 2.64 (1.02–6.83) | Pinteraction | |
T3 | 13 | 3 | 1.42 (0.39–5.10) | = 0.58 | |
No B-cell conditions | T1 | 2,037 | 272 | 1.00 (ref) | |
T2 | 566 | 78 | 1.06 (0.81–1.40) | ||
T3 | 161 | 23 | 1.48 (0.66–1.67) | ||
B-cell conditions | T1 | 16 | 17 | 7.95 (3.79–16.7) | |
T2 | 7 | 5 | 6.07 (1.80–20.5) | Pinteraction | |
T3 | 5 | 2 | 2.59 (0.48–13.9) | = 0.21 | |
No multiple organs | T1 | 2,042 | 272 | 1.00 (ref) | |
T2 | 567 | 79 | 1.08 (0.82–1.42) | ||
T3 | 161 | 23 | 1.05 (0.66–1.67) | ||
Multiple organs | T1 | 19 | 18 | 7.43 (3.74–14.8) | |
T2 | 9 | 4 | 4.90 (1.44–16.7) | Pinteraction | |
T3 | 5 | 2 | 3.24 (0.60–17.5) | =0.27 | |
No organ-specific conditions | T1 | 2,040 | 268 | 1.00 (ref) | |
T2 | 566 | 77 | 1.06 (0.80–1.40) | ||
T3 | 161 | 24 | 1.11 (0.70–1.75) | ||
Organ-specific conditions | T1 | 60 | 5 | 0.59 (0.23–1.51) | |
T2 | 11 | 2 | 1.26 (0.27–5.88) | Pinteraction | |
T3 | 8 | 1 | 0.65 (0.08–5.27) | =0.80 |
Autoimmune conditions . | MZL PRS . | Controls . | MZL . | OR (95% CI) . | Pinteractiona . |
---|---|---|---|---|---|
None | 2,774 | 376 | 1.00 (ref) | ||
Any | 109 | 31 | 1.97 (1.33–2.94) | – | |
No B-cell conditions | 2,764 | 373 | 1.00 (ref) | ||
B-cell conditions | 28 | 24 | 5.88 (3.55–9.74) | – | |
No multiple organs | 2,770 | 374 | 1.00 (ref) | ||
Multiple organs | 33 | 24 | 4.96 (3.04–8.12) | – | |
None/any | T1 | 2,123 | 294 | 1.00 (ref) | |
T2 | 586 | 86 | 1.10 (0.84–1.43) | ||
T3 | 174 | 27 | 1.09 (0.71–1.68) | – | |
None | T1 | 2,046 | 272 | 1.00 (ref) | |
T2 | 567 | 80 | 1.10 (0.83–1.44) | ||
T3 | 161 | 24 | 1.10 (0.70–1.74) | ||
Any | T1 | 77 | 22 | 2.13 (1.28–3.55) | |
T2 | 19 | 6 | 2.64 (1.02–6.83) | Pinteraction | |
T3 | 13 | 3 | 1.42 (0.39–5.10) | = 0.58 | |
No B-cell conditions | T1 | 2,037 | 272 | 1.00 (ref) | |
T2 | 566 | 78 | 1.06 (0.81–1.40) | ||
T3 | 161 | 23 | 1.48 (0.66–1.67) | ||
B-cell conditions | T1 | 16 | 17 | 7.95 (3.79–16.7) | |
T2 | 7 | 5 | 6.07 (1.80–20.5) | Pinteraction | |
T3 | 5 | 2 | 2.59 (0.48–13.9) | = 0.21 | |
No multiple organs | T1 | 2,042 | 272 | 1.00 (ref) | |
T2 | 567 | 79 | 1.08 (0.82–1.42) | ||
T3 | 161 | 23 | 1.05 (0.66–1.67) | ||
Multiple organs | T1 | 19 | 18 | 7.43 (3.74–14.8) | |
T2 | 9 | 4 | 4.90 (1.44–16.7) | Pinteraction | |
T3 | 5 | 2 | 3.24 (0.60–17.5) | =0.27 | |
No organ-specific conditions | T1 | 2,040 | 268 | 1.00 (ref) | |
T2 | 566 | 77 | 1.06 (0.80–1.40) | ||
T3 | 161 | 24 | 1.11 (0.70–1.75) | ||
Organ-specific conditions | T1 | 60 | 5 | 0.59 (0.23–1.51) | |
T2 | 11 | 2 | 1.26 (0.27–5.88) | Pinteraction | |
T3 | 8 | 1 | 0.65 (0.08–5.27) | =0.80 |
aWald test for homogeneity.
All analyses were conducted using SAS 9.3 (SAS Institute, Inc.). All tests were two-sided and P < 0.05 were considered statistically significant. For P-interactions (Pinteractions) < 0.05, Bonferoni adjustment was applied to account for multiple comparisons.
Results
NHL cases and controls representing 14 participating InterLymph Consortium studies in North America, Europe, and Australia were included, comprising 1,914 DLBCL (5,257 controls), 1,733 follicular lymphoma (5,338 controls), and 407 MZL (2,883 controls; Supplementary Table S1). As noted, controls differed by subtype due to differences in the available SNP data to construct the corresponding PRS across participants. Comparisons with our previous publication showed similar distributions of sex, socioeconomic status, and smoking; in addition to being restricted to Caucasian race, the present population was slightly older (Supplementary Table S2).
DLBCL
The highest magnitude of (nongenetic) DLBCL risk was observed for individuals with B-cell–mediated autoimmune conditions (OR = 2.90; 95% CI, 1.94–4.35) or autoimmune conditions that affected multiple organs (OR = 2.24; 95% CI, 1.50–3.33). We also observed an increased DLBCL risk with increasing tertile of DLBCL PRS (vs. T1, T2: OR = 1.24; 95% CI, 1.08–1.43; T3: OR = 1.81; 95% CI, 1.59–2.07; Table 1). DLBCL risk among those with no B-cell–mediated autoimmune conditions increased in a dose-dependent manner for each subsequent level of PRS (Ptrend < 0.0001). Among those with B-cell–mediated autoimmune conditions, ORs were further elevated, with the highest magnitude of DLBCL risk among those with a PRS in T3 (OR = 6.46; 95% CI, 3.04–13.7), compared with the referent group. A similar pattern of elevated DLBCL risk was observed for autoimmune conditions affecting multiple organs and increasing PRS, though we note that B-cell–mediated autoimmune conditions comprise approximately 80% of autoimmune conditions affecting multiple organs (Table 1). There was no statistically significant evidence for interaction with the PRS (Table 1) or in any analyses stratified by the individual SNPs (Supplementary Table S4).
Follicular lymphoma
No association with any category of autoimmune conditions was observed. However, the PRS for follicular lymphoma demonstrated an increasing trend in risk with each PRS tertile (vs. T1, T2: OR = 2.38; 95% CI, 2.03–2.80; T3 OR = 4.89; 95% CI, 4.18–5.73; Ptrend < 0.0001; Table 2). Evaluation of stratified associations showed that the increasing follicular lymphoma risk by PRS did not differ by autoimmune condition category. Evaluation of autoimmune associations by individual loci yielded a significant Pinteraction for rs115374828, whereby associations between autoimmune conditions and follicular lymphoma risk were observed and significant only among those with the CC genotype (Supplementary Table S5). However, we note the Pinteraction was not statistically significant after adjustment for multiple comparisons.
MZL
MZL risk was significantly elevated specifically among those reporting B-cell–mediated autoimmune conditions (OR = 5.88; 95% CI, 3.55–9.74) and conditions that affect multiple organ systems (OR = 4.96; 95% CI, 3.04–8.12; Table 3). We observed no association with increasing tertile of the MZL PRS. Evaluation of autoimmune conditions by PRS yielded statistically significant risk increases for B-cell–mediated conditions and conditions affecting multiple organ systems, regardless of PRS tertile. Further evaluation of associations for autoimmune conditions by individual SNPs did not yield evidence of individual SNP-level interaction (Supplementary Table S6).
Discussion
In this analysis, we observed DLBCL risk to be positively associated with increasing PRS tertiles among those with a history of B-cell–related autoimmune conditions. We did not observe significant joint associations between individual SNPs and autoimmune conditions for increased DLBCL, follicular lymphoma, or MZL risk. Although increased follicular lymphoma risk related to autoimmune conditions appeared restricted to one locus (rs115374828), the Pinteraction was not statistically significant. For MZL, the lack of joint association between the PRS and history of autoimmune conditions was not entirely surprising given the lack of association between the PRS and MZL risk in our population.
As a reflection of the currently established genetic risk, the PRS provides improved statistical power for evaluating joint associations. This power is evident with the significant increase in the magnitudes of risk for each PRS tertile for DLBCL and follicular lymphoma. The suggestive joint association between DLBCL PRS and autoimmune conditions supports the hypothesis that susceptibility loci identified to date independently contribute to DLBCL susceptibility in an additive manner to having autoimmune conditions. It remains to be determined whether their associations are within shared biological pathways, but the additive effects at a minimum suggest that their associations are not overlapping. This is supported by the curious observation that while HLA associations have been reported for both DLBCL and various autoimmune conditions, these associated loci remain distinct for both outcomes. Even each of the purported loci comprising the DLBCL PRS appear to function independently from one another at the moment. Briefly, the susceptibility locus at 6p25.3 (rs116446171) maps near exocyst complex component 2 (EXOC2) which functions at the interface between host defense and cell death regulation and interacts with Ral proteins, which play a crucial role in the maintenance of epithelial cell polarity, cell motility, and cytokinesis, and in proliferation and metastasis (4). The two 8q24.21 variants, rs13255292 and rs4736601, are in close proximity to PVT1, an ncRNA implicated in MYC (oncogene) activation. The susceptibility locus at 2p23.3 (rs79480871) maps near nuclear receptor coactivator 1 (NCOA1) and intersectin 2 (ITSN2). NCOA1 acts as a transcriptional coactivator for steroid and nuclear hormone receptors and ITSN2 encodes a protein involved in clathrin-mediated endocytosis (4). The HLA-B SNP rs2523607 (HLA-B*08:01) plays a central role in presenting intracellularly processed self or foreign antigens to CD8+ cytotoxic T lymphocytes and is carried by the ancestral 8.1 haplotype which has been associated with other autoimmune conditions, such as type-I diabetes. The additive effects between autoimmune conditions and the DLBCL PRS thus appears statistical in nature; larger sample sizes would be required to determine whether biological or multiplicative interaction with specific loci could pinpoint biological pathways of interest and high susceptibility.
Though confirmation is required, the association between rs115374828 and autoimmune conditions and follicular lymphoma is worth noting; rs115374828 reflects HLA DRB1 loci, which has also been implicated in various autoimmune conditions (33, 34), providing plausible biological rationale for a joint association. Although we broadly categorized autoimmune conditions by their purported biological pathways to enhance power, a similar task cannot yet be taken with susceptibility loci, as the culpable loci and their functions remain largely speculative. This is consistent with a recent analysis of InterLymph Conosrtium data by Din and colleagues (35) which showed that the PRS of each of three autoimmune conditions that most consistently demonstrate an association with NHL risk (RA, SLE, MS) was not associated with risk of any major NHL subtype, and that PRS of NHL subtypes were, likewise, not associated with these three autoimmune conditions. With little overall overlap evident among the respective lists of GWAS-identified susceptibility loci, identifying novel ways to pinpoint common but specific biologic processes (e.g., chronic inflammation) may reveal critical interactions and show a causative cascade effect (35). This notion is consistent with the recent report by Ben Eli (2019) that found family history of lymphoma to be associated with NHL, family history of autoimmune conditions to be associated with autoimmune conditions, but family history of autoimmune conditions not to be associated with NHL and vice versa (36)
Study limitations include the use of self-reported autoimmune conditions (although most studies queried for personal history of physician-diagnosed conditions; ref. 3), the low prevalence of autoimmune conditions, and lack of data on severity or treatment; it is possible that true interactions with individual conditions might exist (e.g., for Sjogren syndrome and MZL), for which we lack the power to investigate. Another limitation is the target population in which our study was conducted, as the eligible studies were generally conducted in geographical areas with small minority populations. Our evaluation of PRS for MZL was limited, with only two confirmed SNPs to date and a small sample size. In addition, there is also clinical and biological heterogeneity within MZL that has not yet been evaluated in the context of GWAS. This is also applicable to DLBCL where tumor molecular subtypes such as cell of origin (37) or double hit status (38) used to subclassify DLBCL cases into clinically meaningful categories was not available.
Ongoing efforts include the identification of additional genetic susceptibility loci for NHL subtypes through an expanded GWAS with additional studies. Despite the large sample size compiled from our present consortial collaboration and the strength of our biologically-based groupings of autoimmune conditions to enhance power, even larger sample sizes are needed to confirm these results and to power further exploration of specific autoimmune conditions with specific genetic loci. Moreover, additional loci from these analyses will refine the PRS and may prove instrumental for identifying common pathways of interest that contribute to disease risk. Clinically, leveraging a refined PRS coupled with established risk factors (e.g., autoimmune conditions) may also prove relevant for identifying those with autoimmune conditions who are at highest risk for developing DLBCL based on their genetic profile. It is not implausible to envision enhanced screening among targeted populations, such as those with specific autoimmune conditions known to be at highest risk for DLBCL so that early diagnosis may improve subsequent prognosis. Future efforts to expand the evaluation of NHL risk factors to include potential similar immune responses (e.g., infections) may also be warranted.
Authors' Disclosures
C.M. Vajdic reports grants from Cancer Council NSW during the conduct of the study. A. Nieters reports grants from Jose Carreras Foundation and European Union during the conduct of the study. J.R. Cerhan reports grants from NIH during the conduct of the study as well as grants from Genentech, NanoString, Bristol-Myers Squibb and other support from Bristol-Myers Squibb and Regeneron Genetics Center outside the submitted work. O. Martinez-Maza reports grants from NIH during the conduct of the study. P.M. Bracci reports grants from NIH during the conduct of the study. H. Hjalgrim reports grants from NIH, Danish Cancer Society, and Danish National Research Council during the conduct of the study. C.R. Flowers reports personal fees from AbbVie, Bayer HealthCare, BeiGene, Celgene, Denovo Biopharma, Epizyme/Incyte, Genentech/Roche, Genmab, Gilead, Karyopharm, MEI Pharmaceuticals, MorphoSys AG Pharmacyclics/Janssen, SeaGen, and Spectrum and grants from 4D, AbbVie, Acerta, Adaptimmune, Allogene, Amgen, Bayer HealthCare, Celgene, Cellectis, EMD, Gilead, Genentech/Roche, Guardant, Iovance, Janssen Pharmaceutical, Kite, Morphosys, Nektar, Novartis, Pfizer, Pharmacyclics, Sanofi-Aventis, Takeda, TG Therapeutics, Xencor, Ziopharm, Burroughs Wellcome Fund, Eastern Cooperative Oncology Group, NCI, V Foundation, and Cancer Prevention and Research Institute of Texas: CPRIT Scholar in Cancer Research outside the submitted work. E. Kane reports grants from Blood Cancer UK during the conduct of the study. M. Maynadié reports expertise for Janssen. E.C. Breen reports grants from NIH during the conduct of the study. M.T. Smith reports personal fees from law firms outside the submitted work. E. Roman reports grants from Blood Cancer UK during the conduct of the study. No disclosures were reported by the other authors.
Authors' Contributions
S.S. Wang: Conceptualization, resources, formal analysis, supervision, funding acquisition, methodology, writing–original draft, project administration, writing–review and editing. C.M. Vajdic: Data curation, writing–review and editing. M.S. Linet: Conceptualization, data curation, writing–review and editing. S.L. Slager: Data curation, writing–review and editing. J. Voutsinas: Formal analysis, writing–review and editing. A. Nieters: Conceptualization, data curation, writing–review and editing. D. Casabonne: Data curation, writing–review and editing. J.R. Cerhan: Conceptualization, data curation, writing–review and editing. W. Cozen: Conceptualization, data curation, writing–review and editing. G. Alarcón: Conceptualization, formal analysis, methodology, writing–review and editing. O. Martínez-Maza: Methodology, writing–review and editing. E.E. Brown: Methodology, writing–review and editing. P.M. Bracci: Data curation, writing–review and editing. J. Turner: Data curation, writing–review and editing. H. Hjalgrim: Conceptualization, data curation, methodology, writing–review and editing. P. Bhatti: Data curation, writing–review and editing. Y. Zhang: Data curation, writing–review and editing. B.M. Birmann: Data curation, methodology, writing–review and editing. C.R. Flowers: Conceptualization, methodology, writing–review and editing. O. Paltiel: Methodology, writing–review and editing. E.A. Holly: Data curation, writing–review and editing. E. Kane: Data curation, writing–review and editing. D.D. Weisenburger: Methodology, writing–review and editing. M. Maynadié: Data curation, writing–review and editing. P. Cocco: Data curation, writing–review and editing. L. Foretova: Data curation, writing–review and editing. E.C. Breen: Conceptualization, methodology, writing–review and editing. Q. Lan: Writing-review and editing. A. Brooks-Wilson: Data curation, writing–review and editing. A.J. De Roos: Methodology, writing–review and editing. M.T. Smith: Data curation. E. Roman: Data curation, writing–review and editing. P. Boffetta: Data curation, writing–review and editing. A. Kricker: Data curation, writing–review and editing. T. Zheng: Data curation, writing–review and editing. C.F. Skibola: Writing-review and editing. J. Clavel: Data curation, writing–review and editing. A. Monnereau: Data curation, writing–review and editing. S.J. Chanock: Data curation, funding acquisition, writing–review and editing. N. Rothman: Data curation. Y. Benavente: Data curation, writing–review and editing. P. Hartge: Conceptualization, data curation, methodology, writing–review and editing. K.E. Smedby: Conceptualization, data curation, methodology, writing–original draft, writing–review and editing.
Acknowledgments
This work was supported by NIH grants R03 CA179558 [to S. S. Wang, principal investigator (PI), City of Hope, Duarte, CA]. This work was also supported by:
British Columbia Study: National Cancer Institute of Canada, Canadian Cancer Society; the Canadian Institutes of Health Research (CIHR), and the Michael Smith Foundation for Health Research (to P. Bhatti, A. Brooks-Wilson, PI: John Spinelli).
ENGELA: Association pour la Recherche contre le Cancer (ARC), Institut National du Cancer (INCa), Fondation de France, Fondation contre la Leucémie, Agence nationale de sécurité sanitaire de l'alimentation, de l'environnement et du travail (ANSES; to A. Monnereau, M. Maynadié, J. Clavel).
EpiLymph: European Commission (grant references QLK4-CT-2000-00422 and FOOD-CT-2006-023103); this work was partially supported by the Institut de Salud Carlos III-ISCIII (Spanish Government) cofunded by FEDER funds/European Regional Development Fund (ERDF) – a way to build Europe (CIBERESP CB06/02/0073, PI17/01280, PI20/00288), AGENCIA DE GESTIO D'AJUTS UNIVERSITARIS I DE RECERCA (2017SGR1085), European Commission (grant references 051210) and had no role in the data collection, analysis, or interpretation of the results; the NIH (contract NO1-CO-12400); the Compagnia di San Paolo—Programma Oncologia; the Federal Office for Radiation Protection grants StSch4261 and StSch4420, the José Carreras Leukemia Foundation grant DJCLS-R12/23, the German Federal Ministry for Education and Research (BMBF-01-EO-1303); the Health Research Board, Ireland and Cancer Research Ireland; Czech Republic supported by MH CZ – DRO (MMCI, 00209805) and MEYS – NPS I – LO1413; Fondation de France and Association de Recherche Contre le Cancer (France: M. Maynadié, A. Monnereau; Germany: A. Nieters; Spain: Y. Benavente, D. Casabonne; Czech Republic: L. Foretova; Italy: P. Boffetta).
Mayo Clinic Case-Control Study: NIH (R01 CA92153; R01 CA200703); NCI (P30 CA015083, to J.R. Cerhan, S.L. Slager).
NCI-SEER study: Intramural Research Program of the NCI, NIH, and Public Health Service (N01-PC-65064, N01-PC-67008, N01-PC-67009, N01-PC-67010, N02-PC-71105; to P. Hartge, N. Rothman, S.J. Chanock).
New South Wales study: The Australian National Health and Medical Research Council (ID990920), the Cancer Council NSW, and the University of Sydney Foundation Program (to C.M. Vadjic, J. Turner, A. Kricker).
SCALE: Swedish Cancer Society (2009/659), Stockholm County Council (20110209), and the Strategic Research Program in Epidemiology at Karolinska Institutet; Swedish Cancer Society grant (02 6661); NIH (5R01 CA69669-02); Plan Denmar. (to K.E. Smedby, H. Hjalgrim).
UCSF1 and UCSF2: The UCSF studies were supported by the NCI, NIH (grant nos. CA45614, CA89745, CA87014, CA1046282, and CA154643). The collection of cancer incidence data used in this study was supported by the California Department of Health Services as part of the statewide cancer reporting program mandated by California Health and Safety Code Section 103885; the NCI's SEER Program under contract HHSN261201000140C awarded to the Cancer Prevention Institute of California (to E.A. Holly, P.M. Bracci, C.F Skibola, M.T Smith).
United Kingdom study: Blood Cancer UK (to E. Kane, E. Roman)
Yale University Study: NCI (CA62006 and CA165923; to T. Zheng, Y. Zhang)
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Note: Supplementary data for this article are available at Cancer Epidemiology, Biomarkers & Prevention Online (http://cebp.aacrjournals.org/).