McKay, James

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 Lung cancer (LC) is the primary cause of cancer-related deaths in the United States. Whereas smoking and other environmental factors strongly increase LC risk, multiple genetic variants also contribute to risk in smokers. Furthermore, among smokers, some families have been identified with an abnormally high prevalence of LC, suggesting that unknown genetic factors can greatly increase LC risk in smokers. The typically short survival after LC diagnosis impedes collection of detailed genotypic information on any single large family pedigree, impairing the identification of putative high-risk factors. Therefore, the Genetic Epidemiology of LC Consortium has collected epidemiological and genetic data from a number of families with high numbers of LC cases from eight different sites across the US. In this study, we have obtained whole exome sequences (WES) from 290 members of 28 families, including 66 LC cases. We used a gene-based approach to allow for the possibility that different families may contain different variants of the same gene. Variants were filtered for i) allele frequency, ii) functional effect using combined annotation-dependent depletion (CADD), and, iii) affecting a gene with either a known or suspected role in cancer. We further selected variants based on their segregation in family pedigrees in a pattern consistent with a large effect on LC risk. Candidate LC risk genes were then identified as those represented in at least two families by the same or different variants. We further culled the list of genes by requiring the presence of at least one rare, functional variant enriched in the WES of 1060 cases relative to 899 controls from the Transdisciplinary Research on Cancer of the Lung consortium. This analysis narrowed our results to two genes, one being E2A, a member of the E family of bHLH transcription factors. Whereas loss-of-function mutations in E2A drive lymphoid cancers, the E2A protein also participates in an oncogenic heterodimer with TWIST1 that promotes the epithelial-mesenchymal transition and is implicated in multiple cancer types. Furthermore, the E2A/TWIST1 heterodimer is the primary TWIST1-containing complex implicated in oncogenesis, and silencing of E2A in KRAS-mutant non-small cell lung cancer (NSCLC) cell lines results in oncogene-stimulated senescence and apoptosis. Our data identified three distinct E2A variants present in all 10 sequenced LC cases in the 5 families in which those variants are found. Finally, two of these E2A variants are located only 57 nucleotides from each other, immediately adjacent to sequences encoding a transcription activation domain, suggesting that both variants alter the same specific protein function. These data identify E2A variants as likely high risk factors for LC in smokers and validate our general approach for identifying genetic factors with a large impact on LC risk. Citation Format: Claudio W. Pikielny, Anthony M. Musolf, Mariza de Andrade, Diptasri Mandal, Colette Gaba, Ping Yang, yafang Li, Ming You, Richard Wilson, Elena Y. Kupert, Marshall W. Anderson, Ann G. Schwartz, Susan M. Pinney, Ambrose I. Granizo-Mackenzie, Yanhong Liu, Ramaswamy Govindan, James McKay, Rayjean Hung, John K. Field, David C. Christiani, Joan E. Bailey-Wilson, Christopher I. Amos. Familial studies identify variants in the E2A transcription factor as putative risk factors for lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr LB-053.

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.

Key Points PRS, based on the known CLL loci, predicts CLL risk with high discrimination. This PRS predicts risk of monoclonal B-cell lymphocytosis, a precursor to CLL and a condition that has clinical impact beyond risk for CLL.

Abstract Given the frequent and largely incurable occurrence of multiple myeloma, identification of germline genetic mutations that predispose cells to multiple myeloma may provide insight into disease etiology and the developmental mechanisms of its cell of origin, the plasma cell (PC). Here, we identified familial and early-onset multiple myeloma kindreds with truncating mutations in lysine-specific demethylase 1 (LSD1/KDM1A), an epigenetic transcriptional repressor that primarily demethylates histone H3 on lysine 4 and regulates hematopoietic stem cell self-renewal. In addition, we found higher rates of germline truncating and predicted deleterious missense KDM1A mutations in patients with multiple myeloma unselected for family history compared with controls. Both monoclonal gammopathy of undetermined significance (MGUS) and multiple myeloma cells have significantly lower KDM1A transcript levels compared with normal PCs. Transcriptome analysis of multiple myeloma cells from KDM1A mutation carriers shows enrichment of pathways and MYC target genes previously associated with myeloma pathogenesis. In mice, antigen challenge followed by pharmacologic inhibition of KDM1A promoted PC expansion, enhanced secondary immune response, elicited appearance of serum paraprotein, and mediated upregulation of MYC transcriptional targets. These changes are consistent with the development of MGUS. Collectively, our findings show that KDM1A is the first autosomal-dominant multiple myeloma germline predisposition gene providing new insights into its mechanistic roles as a tumor suppressor during post-germinal center B-cell differentiation. Significance: KDM1A is the first germline autosomal dominant predisposition gene identified in multiple myeloma and provides new insights into multiple myeloma etiology and the mechanistic role of KDM1A as a tumor suppressor during post-germinal center B-cell differentiation. Cancer Res; 78(10); 2747–59. ©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 Waldenström macroglobulinemia (WM) is a unique subset of lymphoplasmacytic lymphoma (LPL) that is defined by the presence of an LPL infiltrate in the bone marrow together with a monoclonal IgM protein in the serum. A somatic activating mutation, MYD88 L265P, occurs in 85+% of WM and in 25%-50% of patients with the precursor condition, IgM monoclonal gammopathy of undetermined significance (MGUS); however, germline MYD88 mutations have not been observed in WM patients, and the genetic basis for WM predisposition remains undefined. To identify novel WM susceptibility loci we conducted a two-stage genome-wide association study (GWAS) in over 450 WM cases and 4300 controls of European ancestry. Discovery (stage 1) included 217 WM cases (40% familial) and 3798 controls genotyped on the Illumina Omni Express or Illumina Omni2.5 platforms following standard quality control procedures. The genotyped data were imputed using the Haplotype Reference Consortium panel as a reference and analyzed using logistic regression. In stage 1, we identified three loci on chromosomes 6, 14 and 3 significantly associated (P<5.0x10-8) with risk of WM. Eleven promising SNPs in these and other suggestive loci (P<5.0x10-7) were selected for replication (stage 2) in 269 WM or LPL cases (4% familial) and 571 controls, and genotyping was conducted using standard methods on Taqman and Sequenom platforms or Sanger sequencing (1 SNP). Preliminary results confirm replication of the chromosome 6 locus. Stratification on familial status will illuminate the contribution of familial disease. These results will provide insight into the underlying genetic basis of WM susceptibility. Citation Format: Mary L. McMaster, Sonja I. Berndt, Shengchao A. Li, Susan Slager, Joseph Vijai, Charles C. Chung, Bin Zhu, Laurie Burdette, Brenda Birmann, Elizabeth E. Brown, James R. Cerhan, Karin Ekstrom-Smedby, Henrik Hjalgrim, Geffen Kleinstern, Brian K. Link, James McKay, Alain Monnereau, Lindsay M. Morton, Alexandra Nieters, Nathaniel Rothman, Christine F. Skibola, Alex Smith, Lauren R. Teras, Claire M. Vajdic, Roel Vermeulen, Belynda Hicks, Lynn R. Goldin, Neil E. Caporaso. A genome-wide association study of Waldenström macroglobulinemia/lymphoplasmacytic lymphoma demonstrates association with chromosome 6 [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1318. doi:10.1158/1538-7445.AM2017-1318