Cox, David G.

Abstract Obesity has been associated with non‐Hodgkin lymphoma (NHL), but the evidence is inconclusive. We examined the association between genetically determined adiposity and four common NHL subtypes: diffuse large B‐cell lymphoma (DLBCL), follicular lymphoma, chronic lymphocytic leukemia, and marginal zone lymphoma, using eight genome‐wide association studies of European ancestry ( N = 10,629 cases, 9505 controls) and constructing polygenic scores for body mass index (BMI), waist‐to‐hip ratio (WHR), and waist‐to‐hip ratio adjusted for BMI (WHRadjBMI). Higher genetically determined BMI was associated with an increased risk of DLBCL [odds ratio (OR) per standard deviation (SD) = 1.18, 95% confidence interval (95% CI): 1.05–1.33, p = .005]. This finding was consistent with Mendelian randomization analyses, which demonstrated a similar increased risk of DLBCL with higher genetically determined BMI (OR per SD = 1.12, 95% CI: 1.02–1.23, p = .03). No significant associations were observed with other NHL subtypes. Our study demonstrates a positive link between a genetically determined BMI and an increased risk of DLBCL, providing additional support for increased adiposity as a risk factor for DLBCL.

Abstract Background Relatively few studies have evaluated the role of host (germline) genetics in predicting outcome in diffuse large B-cell lymphoma (DLBCL), and of these, most were small candidate gene studies without replication and/or from the pre-R-CHOP era. We report the first genome-wide association study (GWAS) of outcome in DLBCL patients treated with rituximab and anthracycline-based chemotherapy. We used two independent studies for discovery, and validated the results in two additional independent studies. Methods In the discovery phase, we conducted a meta-analysis of two GWAS studies of DLBCL patients treated with rituximab and anthracycline-based chemotherapy, the first study as part of the prospective clinical trials of the LNH2003B program (N=540) of the Lymphoma Study Association (LYSA) and the second as part of the Molecular Epidemiology Resource (N=316), a prospective observational study from the University of Iowa/Mayo Clinic Lymphoma Specialized Program of Research Excellence (SPORE). Genotyping for the GWAS was performed using the Illumina 610 (LYSA) and 660 Quad (SPORE) BeadArrays (Illumina, San Diego, USA). HapMap2 was used to impute additional SNPs. Association testing was conducted assuming a log-additive genetic model and using Cox regression analyses. We calculated the hazard ratios (HRs) and 95% CIs for event-free survival (EFS), with adjustment for age, sex and age-adjusted IPI. The two studies were then combined using the fixed effects inverse variance method based on the beta estimates and standard errors from each study. From the meta-analysis, 89 SNPs were selected for validation, based on a minor allele frequency (MAF) >0.05 and LD pattern for multiple SNPs from the same region. Validation cohorts were an independent set of DLBCL patients from the SPORE (SPORE-2, N=366) and DLBCL patients from the prospective GOELAMS 075 phase III study (N=294), and validation genotyping was conducted using the Illumina VeraCode platform. We then conducted a meta-analysis of the 4 studies. We report SNPs with both a statistical significance of p<5x10-6 and effect sizes in the same direction across the 4 studies in the meta-analysis. For the latter SNPs, we also report the results for overall survival (OS), adjusted for age, sex and age-adjusted IPI. Results In the final meta-analysis of all four studies for the prediction of EFS, the top SNP was rs7765004 (MAF=0.30; HR=1.45; p=9.7 x 10-8) in the 6q21 region, which is near the MARCKS (myristoylated alanine-rich protein kinase C substrate) and HDAC2 genes. MARCKS encodes for a protein involved in cell mobility, while HDAC2 belongs to the histone deacetylase family. Additional SNPs from this region included rs6918103 (MAF=0.25; HR=1.38; p=7.6 x 10-6), which marks an EZH2 histone binding site, and rs11969684 (MAF=0.12; HR=1.43; p=3.3 x 10-5). The next hit was rs7712513 (MAF=0.31; HR=1.40; p=3.2 x 10-7), which is near arginine-fifty homeobox pseudogene 1 (ARGFXP1) on chromosome 5q23.2. This was followed by a region on chromosome 8q13.3, with the SNP rs9298183 (MAF=0.29; HR=1.41; p=1.4 x 10-6) from LOC100132891 showing the strongest association, along with 3 other SNPs from this region (rs2035376, rs6472637, and rs1600797; all p<5 x 10-6). These SNPs are near the musculin gene (MSC) a transcriptional repressor of E2A frequently silenced by promoter methylation in DLBCL. Finally, rs2253986 (MAF=0.27; HR=1.39; p=3.4 x 10-6) at chromosome 16q12.2, along with rs2253971 and rs6499810, were associated with EFS. These SNPs are located in CES5A (carboxylesterase 5A), which encodes a member of the carboxylesterase family. The family members are responsible for the hydrolysis or transesterification of various xenobiotics. The top SNPs associated with EFS also predicted OS, including rs7765004 (HR=1.41; p=1.2 x 10-5), rs7712513 (HR=1.39; p=8.4 x 10-6), rs9298183 (HR=1.36; p=3.5 x 10-5), and rs2253986 (HR=1.35; p=5.2 x 10-5). Conclusions Although our top SNPs on did not reach genome-wide significance (p<5x10-8), they showed clear consistency across the four large prospective studies and predicted both EFS and OS independently of IPI. The strongest findings were from loci on 6q21, 8q13.3, 5q23.2, and 16q12.2, and several genes from these regions are involved in methylation and xenobiotic metabolism, supporting a role for host genetic variation in these pathways in DLBCL prognosis. Disclosures: No relevant conflicts of interest to declare.

AbstractEpidemiologic studies show an increased risk of non‐Hodgkin lymphoma (NHL) in patients with autoimmune disease (AD), due to a combination of shared environmental factors and/or genetic factors, or a causative cascade: chronic inflammation/antigen‐stimulation in one disease leads to another. Here we assess shared genetic risk in genome‐wide‐association‐studies (GWAS).Secondary analysis of GWAS of NHL subtypes (chronic lymphocytic leukemia, diffuse large B‐cell lymphoma, follicular lymphoma, and marginal zone lymphoma) and ADs (rheumatoid arthritis, systemic lupus erythematosus, and multiple sclerosis). Shared genetic risk was assessed by (a) description of regional genetic of overlap, (b) polygenic risk score (PRS), (c)"diseasome", (d)meta‐analysis.Descriptive analysis revealed few shared genetic factors between each AD and each NHL subtype. The PRS of ADs were not increased in NHL patients (nor vice versa). In the diseasome, NHLs shared more genetic etiology with ADs than solid cancers (p = .0041). A meta‐analysis (combing AD with NHL) implicated genes of apoptosis and telomere length.This GWAS‐based analysis four NHL subtypes and three ADs revealed few weakly‐associated shared loci, explaining little total risk. This suggests common genetic variation, as assessed by GWAS in these sample sizes, may not be the primary explanation for the link between these ADs and NHLs.

Abstract A growing number of loci within the human leukocyte antigen (HLA) region have been implicated in non-Hodgkin lymphoma (NHL) etiology. Here, we test a complementary hypothesis of “heterozygote advantage” regarding the role of HLA and NHL, whereby HLA diversity is beneficial and homozygous HLA loci are associated with increased disease risk. HLA alleles at class I and II loci were imputed from genome-wide association studies (GWAS) using SNP2HLA for 3,617 diffuse large B-cell lymphomas (DLBCL), 2,686 follicular lymphomas (FL), 2,878 chronic lymphocytic leukemia/small lymphocytic lymphomas (CLL/SLL), 741 marginal zone lymphomas (MZL), and 8,753 controls of European descent. Both DLBCL and MZL risk were elevated with homozygosity at class I HLA-B and -C loci (OR DLBCL = 1.31, 95% CI = 1.06–1.60; OR MZL = 1.45, 95% CI = 1.12–1.89) and class II HLA-DRB1 locus (OR DLBCL = 2.10, 95% CI = 1.24–3.55; OR MZL = 2.10, 95% CI = 0.99–4.45). Increased FL risk was observed with the overall increase in number of homozygous HLA class II loci (P trend &lt; 0.0001, FDR = 0.0005). These results support a role for HLA zygosity in NHL etiology and suggests that distinct immune pathways may underly the etiology of the different NHL subtypes. Significance: HLA gene diversity reduces risk for non-Hodgkin lymphoma. Cancer Res; 78(14); 4086–96. ©2018 AACR.

Objective Determinants of the increased risk of diffuse large B-cell lymphoma (DLBCL) in SLE are unclear. Using data from a recent lymphoma genome-wide association study (GWAS), we assessed whether certain lupus-related single nucleotide polymorphisms (SNPs) were also associated with DLBCL. Methods GWAS data on European Caucasians from the International Lymphoma Epidemiology Consortium (InterLymph) provided a total of 3857 DLBCL cases and 7666 general-population controls. Data were pooled in a random-effects meta-analysis. Results Among the 28 SLE-related SNPs investigated, the two most convincingly associated with risk of DLBCL included the CD40 SLE risk allele rs4810485 on chromosome 20q13 (OR per risk allele=1.09, 95% CI 1.02 to 1.16, p=0.0134), and the HLA SLE risk allele rs1270942 on chromosome 6p21.33 (OR per risk allele=1.17, 95% CI 1.01 to 1.36, p=0.0362). Of additional possible interest were rs2205960 and rs12537284. The rs2205960 SNP, related to a cytokine of the tumour necrosis factor superfamily TNFSF4, was associated with an OR per risk allele of 1.07, 95% CI 1.00 to 1.16, p=0.0549. The OR for the rs12537284 (chromosome 7q32, IRF5 gene) risk allele was 1.08, 95% CI 0.99 to 1.18, p=0.0765. Conclusions These data suggest several plausible genetic links between DLBCL and SLE.

Abstract Background: Unnecessary intervention and overtreatment of indolent disease are common challenges in clinical management of prostate cancer. Improved tools to distinguish lethal from indolent disease are critical. Methods: We performed a genome-wide survival analysis of cause-specific death in 24,023 prostate cancer patients (3,513 disease-specific deaths) from the PRACTICAL and BPC3 consortia. Top findings were assessed for replication in a Norwegian cohort (CONOR). Results: We observed no significant association between genetic variants and prostate cancer survival. Conclusions: Common genetic variants with large impact on prostate cancer survival were not observed in this study. Impact: Future studies should be designed for identification of rare variants with large effect sizes or common variants with small effect sizes. Cancer Epidemiol Biomarkers Prev; 24(11); 1796–800. ©2015 AACR.

Abstract Background: Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma subtype and is clinically aggressive. Relatives of DLBCL patients are at elevated risk for DLBCL. In small genome-wide association studies (GWAS) in which all NHL subtypes were combined, no conclusive loci for DLBCL have been previously identified in individuals of European background, whereas a recent study conducted in East Asia identified a locus at 3q27. Methods: To identify genetic susceptibility loci for DLBCL, we conducted a meta-analysis of three new GWAS and one prior scan, totaling 3,857 DLBCL cases and 7,666 controls of European ancestry. Imputation was used to combine all the data using the 1000 Genomes Project release v3 and IMPUTE2. Based on these findings, we selected 6 promising SNPs for genotyping in an additional 202 cases and 3,431 controls, which were then included in a second meta-analysis. SNPTEST was used to estimate odds ratios (ORs), using the additive model and adjusting for age, gender and significant eigenvectors. Results: In our multi-stage analysis, five independent SNPs in four loci achieved genome-wide significance marked by rs116446171 at 6p25.3 (EXOC2; OR=2.33, P=1.60x10-20), rs13255292 and rs4733601 at 8q24.21 (PVT1; OR=1.23, P=4.86x10-12 and OR=1.18, P=1.63x10-8, respectively), rs2523607 at 6p21.33 (HLA-B; OR=1.34, 3.35x10-9) and rs2681416 at 3q13.33 (CD86; OR=1.20, P=6.04x10-9). The two 8q24.21 SNPs displayed minimal linkage disequilibrium (r2=0.03) and conditional analysis supported their independence. We did not observe a notable signal for a locus previously reported for DLBCL on 3q27 in East Asia, rs6773854 (OR=1.06, P=0.81). We estimated that common SNPs, including but not limited to the loci discovered in this study, explain approximately 16% of the variance for DLBCL risk overall. Discussion: The susceptibility locus at 6p25.3 maps near a plausible DLBCL candidate gene, EXOC2 (exocyst complex component 2), which interacts with Ral proteins, and disruption of the Ral-exocyst regulatory node impacts cellular and developmental processes associated with malignant transformation and progression. The two 8q24.21 variants are approximately 1Mb telomeric to the 8q24 region linked with multiple cancers including CLL. Both variants are positioned in close proximity to PVT1, which is a non-coding RNA implicated in the MYC activation, and MYC is known to be deregulated in some DLBCLs. The locus at 3q13.33 maps to CD86 (B7-2), a B-7 family member that along with CD80 (B7-1) is essential in cognate interactions between T- and B-lymphocytes; both play key roles in immunodeficiency, autoimmune diseases, and in anti-tumor immunity. Finally, we had multiple signals in the HLA region, the strongest of which was at 6p21.33 near HLA-B. HLA-B encodes the HLA class I heavy chain paralogue, which heterodimerizes with a light chain (β2 microglobulin) and plays a central role in presenting intracellularly processed self or foreign antigens to CD8+ cytotoxic T lymphocytes. Class I molecules have been extensively linked to the risk of a variety of immune-mediated diseases and cancers including DLBCL. Conclusions: These data provide substantial new evidence for genetic susceptibility to this B-cell malignancy, and point towards pathways involved in immune recognition and immune function in the pathogenesis of DLBCL. Citation Format: James R. Cerhan, Sonja I. Berndt, Joseph Vijai, Hervé Ghesquières, James McKay, Sophia S. Wang, Zhaoming Wang, Meredith Yeager, Alexandra Nieters, David Cox, Alain Monnereau, Christopher Flowers, Anneclaire J. De Roos, Angela R. Brooks-Wilson, Qing Lan, Gianluca Severi, Mads Melbye, Rebecca D. Jackson, Lauren R. Teras, Mark Purdue, Claire Vajdic, Demetrius Albanes, Kimberly A. Bertrand, Anne Zeleniuch-Jacquotte, Simon Crouch, Yawei Zhang, Paolo Vineis, Susan L. Slager, Karin E. Smedby, Gilles Salles, Christine F. Skibola, Nathaniel Rothman, Stephen J. Chanock, on behalf of the NHL GWAS Project. Genome-wide association study identifies multiple susceptibility loci for diffuse large B-cell lymphoma. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr LB-272. doi:10.1158/1538-7445.AM2014-LB-272

Abstract Previous studies have suggested that minor alleles for ERCC4 rs744154, TNF rs361525, CASP10 rs13010627, PGR rs1042838, and BID rs8190315 may influence breast cancer risk, but the evidence is inconclusive due to their small sample size. These polymorphisms were genotyped in more than 30,000 breast cancer cases and 30,000 controls, primarily of European descent, from 30 studies in the Breast Cancer Association Consortium. We calculated odds ratios (OR) and 95% confidence intervals (95% CI) as a measure of association. We found that the minor alleles for these polymorphisms were not related to invasive breast cancer risk overall in women of European descent: ECCR4 per-allele OR (95% CI) = 0.99 (0.97-1.02), minor allele frequency = 27.5%; TNF 1.00 (0.95-1.06), 5.0%; CASP10 1.02 (0.98-1.07), 6.5%; PGR 1.02 (0.99-1.06), 15.3%; and BID 0.98 (0.86-1.12), 1.7%. However, we observed significant between-study heterogeneity for associations with risk for single-nucleotide polymorphisms (SNP) in CASP10, PGR, and BID. Estimates were imprecise for women of Asian and African descent due to small numbers and lower minor allele frequencies (with the exception of BID SNP). The ORs for each copy of the minor allele were not significantly different by estrogen or progesterone receptor status, nor were any significant interactions found between the polymorphisms and age or family history of breast cancer. In conclusion, our data provide persuasive evidence against an overall association between invasive breast cancer risk and ERCC4 rs744154, TNF rs361525, CASP10 rs13010627, PGR rs1042838, and BID rs8190315 genotypes among women of European descent. (Cancer Epidemiol Biomarkers Prev 2009;18(5):1610–6)