Purdue, Mark P.

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 The shared inherited genetic contribution to risk of different cancers is not fully known. In this study, we leverage results from 12 cancer genome-wide association studies (GWAS) to quantify pairwise genome-wide genetic correlations across cancers and identify novel cancer susceptibility loci. Methods We collected GWAS summary statistics for 12 solid cancers based on 376 759 participants with cancer and 532 864 participants without cancer of European ancestry. The included cancer types were breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancers. We conducted cross-cancer GWAS and transcriptome-wide association studies to discover novel cancer susceptibility loci. Finally, we assessed the extent of variant-specific pleiotropy among cancers at known and newly identified cancer susceptibility loci. Results We observed widespread but modest genome-wide genetic correlations across cancers. In cross-cancer GWAS and transcriptome-wide association studies, we identified 15 novel cancer susceptibility loci. Additionally, we identified multiple variants at 77 distinct loci with strong evidence of being associated with at least 2 cancer types by testing for pleiotropy at known cancer susceptibility loci. Conclusions Overall, these results suggest that some genetic risk variants are shared among cancers, though much of cancer heritability is cancer-specific and thus tissue-specific. The increase in statistical power associated with larger sample sizes in cross-disease analysis allows for the identification of novel susceptibility regions. Future studies incorporating data on multiple cancer types are likely to identify additional regions associated with the risk of multiple cancer types.

Abstract Background Mendelian randomization (MR) studies are susceptible to metadata errors (e.g. incorrect specification of the effect allele column) and other analytical issues that can introduce substantial bias into analyses. We developed a quality control (QC) pipeline for the Fatty Acids in Cancer Mendelian Randomization Collaboration (FAMRC) that can be used to identify and correct for such errors. Methods We collated summary association statistics from fatty acid and cancer genome-wide association studies (GWAS) and subjected the collated data to a comprehensive QC pipeline. We identified metadata errors through comparison of study-specific statistics to external reference data sets (the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue and 1000 genome super populations) and other analytical issues through comparison of reported to expected genetic effect sizes. Comparisons were based on three sets of genetic variants: (i) GWAS hits for fatty acids, (ii) GWAS hits for cancer and (iii) a 1000 genomes reference set. Results We collated summary data from 6 fatty acid and 54 cancer GWAS. Metadata errors and analytical issues with the potential to introduce substantial bias were identified in seven studies (11.6%). After resolving metadata errors and analytical issues, we created a data set of 219 842 genetic associations with 90 cancer types, generated in analyses of 566 665 cancer cases and 1 622 374 controls. Conclusions In this large MR collaboration, 11.6% of included studies were affected by a substantial metadata error or analytical issue. By increasing the integrity of collated summary data prior to their analysis, our protocol can be used to increase the reliability of downstream MR analyses. Our pipeline is available to other researchers via the CheckSumStats package (https://github.com/MRCIEU/CheckSumStats).

Abstract Background: Lipid traits have been inconsistently linked to risk of non-Hodgkin lymphoma (NHL). We examined the association of genetically predicted lipid traits with risk of diffuse large B-cell lymphoma (DLBCL), chronic lymphocytic leukemia (CLL), follicular lymphoma (FL), and marginal zone lymphoma (MZL) using Mendelian randomization (MR) analysis. Methods: Genome-wide association study data from the InterLymph Consortium were available for 2,661 DLBCLs, 2,179 CLLs, 2,142 FLs, 824 MZLs, and 6,221 controls. SNPs associated (P < 5 × 10−8) with high-density lipoprotein (HDL, n = 164), low-density lipoprotein (LDL, n = 137), total cholesterol (TC, n = 161), and triglycerides (TG, n = 123) were used as instrumental variables (IV), explaining 14.6%, 27.7%, 16.8%, and 12.8% of phenotypic variation, respectively. Associations between each lipid trait and NHL subtype were calculated using the MR inverse variance–weighted method, estimating odds ratios (OR) per standard deviation and 95% confidence intervals (CI). Results: HDL was positively associated with DLBCL (OR = 1.14; 95% CI, 1.00–1.30) and MZL (OR = 1.09; 95% CI, 1.01–1.18), while TG was inversely associated with MZL risk (OR = 0.90; 95% CI, 0.83–0.99), all at nominal significance (P < 0.05). A positive trend was observed for HDL with FL risk (OR = 1.08; 95% CI, 0.99–1.19; P = 0.087). No associations were noteworthy after adjusting for multiple testing. Conclusions: We did not find evidence of a clear or strong association of these lipid traits with the most common NHL subtypes. While these IVs have been previously linked to other cancers, our findings do not support any causal associations with these NHL subtypes. Impact: Our results suggest that prior reported inverse associations of lipid traits are not likely to be causal and could represent reverse causality or confounding.

Abstract We previously identified five single nucleotide polymorphisms (SNPs) at four susceptibility loci for diffuse large B-cell lymphoma (DLBCL) in individuals of European ancestry through a large genome-wide association study (GWAS). To further elucidate genetic susceptibility to DLBCL, we sought to validate two loci at 3q13.33 and 3p24.1 that were suggestive in the original GWAS with additional genotyping. In the meta-analysis (5662 cases and 9237 controls) of the four original GWAS discovery scans and three replication studies, the 3q13.33 locus (rs9831894; minor allele frequency [MAF] = 0.40) was associated with DLBCL risk [odds ratio (OR) = 0.83, P = 3.62 × 10−13]. rs9831894 is in linkage disequilibrium (LD) with additional variants that are part of a super-enhancer that physically interacts with promoters of CD86 and ILDR1. In the meta-analysis (5510 cases and 12 817 controls) of the four GWAS discovery scans and four replication studies, the 3p24.1 locus (rs6773363; MAF = 0.45) was also associated with DLBCL risk (OR = 1.20, P = 2.31 × 10−12). This SNP is 29 426-bp upstream of the nearest gene EOMES and in LD with additional SNPs that are part of a highly lineage-specific and tumor-acquired super-enhancer that shows long-range interaction with AZI2 promoter. These loci provide additional evidence for the role of immune function in the etiology of DLBCL, the most common lymphoma subtype.

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 Background: We previously identified 5 SNPs at 4 susceptibility loci for diffuse large B-cell lymphoma (DLBCL) in individuals of European ancestry through a large genome-wide association study (Cerhan et al, Nat Gen 2014;46:1233-8); however, much of the heritability remains unexplained. To further elucidate genetic susceptibility to DLBCL, we sought to validate 2 loci at 3q13.33 and 3p24.1 that were suggestive in the original report. Methods: We selected two SNPs rs9831894 (3q13.33) and rs6773363 (3p24.1) for de novo replication from studies at the Mayo Clinic, MD Anderson, and Memorial Sloan Kettering. Logistic regression was used to estimate odds ratios (ORs), using a log-additive model, and adjusting for age, gender and significant eigenvectors for each genotyping center separately. For both SNPs, we then conducted a meta-analysis of all replication studies (n=3) and our original GWAS studies (n=4), encompassing in total 5662 cases and 9237 controls for rs9831894, and 5510 cases and 12,817 controls for rs6773363. Meta-analysis was conducted using the fixed-effects inverse variance method based on the β estimates and standard errors from each study. Results: In a meta-analysis of the 4 original GWAS scans (3856 cases and 7666 controls), rs9831894 (MAF=0.40) was associated with DLBCL risk (OR=0.84, P=4.52 x 10-9) and this replicated in a meta-analysis of the 3 studies with de novo genotyping (OR=0.80, P=4.17 x 10-5); the overall meta-analysis showed a strong association with DLBCL risk (OR=0.83, P=3.62 x 10-13). This locus maps near a plausible candidate gene, CD86, a protein coding gene and a member of the immunoglobulin superfamily that encodes a type I membrane protein. Binding of CD86 with cytotoxic T-lymphocyte-associated protein 4 negatively regulates T-cell activation and diminishes the immune response, while binding of CD86 with CD28 antigen is a costimulatory signal for activation of the T-cell. In a meta-analysis of the same 4 GWAS scans, rs6773363 (MAF=0.45) was tentatively associated with DLBCL risk (OR=1.17, P= 3.68 x 10-7) and this replicated in a meta-analysis of the studies with de novo genotyping (OR=1.27, P=3.78 x 10-7); the overall meta-analysis showed a strong association with DLBCL risk (OR=1.20, P=2.31 x 10-12). This locus also maps near a plausible candidate gene, eomesodermin (EOMES), a transcription factor crucial for embryonic development of the central nervous system and also putatively involved in T-cell differentiation in viral infection defense. Conclusion: In this follow-up analysis of our initial GWAS, we have identified two additional loci associated with risk of DLBCL, the most common lymphoma subtype. These loci provide additional evidence for the role of immune function in the etiology of DLBCL. Citation Format: Geffen Kleinstern, Michelle Hildebrand, Vijai Joseph, Sonja I. Berndt, Hervé Ghesquières, James McKay, Sophia S. Wang, Alexandra Nieters, David Cox, Alain Monnereau, Angela R. Brooks-Wilson, Qing Lan, Mads Melbye, Rebecca D. Jackson, Lauren R. Teras, Mark P. Purdue, Claire M. Vajdic, Demetrius Albanes, Anne Zeleniuch-Jacquotte, Simon Crouch, Yawei Zhang, Susan L. Slager, Xifeng Wu, Karin E. Smedby, Gilles Salles, Christine F. Skibola, Nathaniel Rothman, Stephen J. Chanock, James R. Cerhan. Inherited variants at 3q13.33 and 3p24.1 influences risk of diffuse large B-cell lymphoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 221.

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 < 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.