Drucker, Kristen

Establishing causal links between inherited polymorphisms and cancer risk is challenging. Here, we focus on the single-nucleotide polymorphism rs55705857, which confers a sixfold greater risk of isocitrate dehydrogenase ( IDH) –mutant low-grade glioma (LGG). We reveal that rs55705857 itself is the causal variant and is associated with molecular pathways that drive LGG. Mechanistically, we show that rs55705857 resides within a brain-specific enhancer, where the risk allele disrupts OCT2/4 binding, allowing increased interaction with the Myc promoter and increased Myc expression. Mutating the orthologous mouse rs55705857 locus accelerated tumor development in an Idh1 R132H -driven LGG mouse model from 472 to 172 days and increased penetrance from 30% to 75%. Our work reveals mechanisms of the heritable predisposition to lethal glioma in ~40% of LGG patients.

Abstract Knowledge about inherited and acquired genetics of adult diffuse glioma has expanded significantly over the past decade. Genomewide association studies (GWAS) stratified by histologic subtype identified six germline variants that were associated specifically with glioblastoma (GBM) and 12 that were associated with lower grade glioma. A GWAS performed using the 2016 WHO criteria, stratifying patients by IDH mutation and 1p/19q codeletion (as well as TERT promoter mutation), discovered that many of the known variants are associated with specific WHO glioma subtypes. In addition, the GWAS stratified by molecular group identified two additional novel regions: variants in D2HGDH that were associated with tumors that had an IDH mutation and a variant near FAM20C that was associated with tumors that had both IDH mutation and 1p/19q codeletion. The results of these germline associations have been used to calculate polygenic risk scores, from which to estimate relative and absolute risk of overall glioma and risk of specific glioma subtypes. We will review the concept of polygenic risk models and their potential clinical utility, as well as discuss the published adult diffuse glioma polygenic risk models. To date, these prior genetic studies have been done on European populations. Using the published glioma polygenic risk model, we show that the genetic associations published to date do not generalize across genetic ancestries, demonstrating that genetic studies need to be done on more diverse populations.

Abstract Low-grade glioma (LGG) are generally slowly growing brain cancers, that frequently undergo malignant progression to aggressive, secondary glioblastoma with a dismal prognosis. By combining genetically engineered Idh1-mutant mice with in vivo CRISPR gene editing we generated a mouse model faithfully recapitulating the founder mutations of LGG. Clonal activation of the neomorphic Idh1 R132H mutation cooperates with Trp53 and Atrx mutations to trigger development of brain tumors but only with ~30% penetrance and very long latency. To elucidate the molecular mechanisms underlying the malignant progression of IDH1-mutant LGG, we devised and deployed a direct in vivo CRISPR screen targeting genes commonly mutated in human IDH-mutant secondary glioblastoma. Stereotaxic delivery of a lentiviral sgRNA library targeting the mouse orthologs of these genes into the brain of Idh1 R132H ;Trp53;Atrx;Cas9 and control Idh1 wt ;Trp53;Atrx;Cas9 compound mutant mice resulted in rapid formation of tumors that recapitulate human Idh1-mutant glioblastoma. Deconvoluting the screen showed that PI3K pathway members Pten and Pik3ca as well as Notch1, Smarca4 and Fat1 are preferentially enriched in Idh1 R132H-tumors, while Rb1 and NF2 were enriched in Idh1 wt tumors. Co-mutation analysis further identified additional co-occurring driver combinations such as Bcor-Met, Olig2-Met, Olig2-Med12 or Bcor-Olig2. We validated the tumor suppressive function of Notch1 and Pten using conventional floxed knock-out alleles and found that Notch1 functions in a haploinsufficient manner. Interestingly, Idh1 R132H did not alter tumor latency or pathology in a high grade p53;Pten;Rb1 mutant background, indicating that the neomorphic IDH-mutations can drive low but not high grade glioma development. Our study provides a functional landscape of gliomagenesis suppressors in vivo.

Abstract Abnormalities on brain MRI may be identified in ~15% of individuals aged 50-66. It can be difficult to discriminate glioma, CNS lymphoma, Inflammatory Demyelinating Disease (CNSIDD) and solitary metastasis, and misdiagnosis may expose patients to unnecessary anxiety, surgery, or radiotherapy. CNS lymphoma requires only biopsy, solitary metastasis may be resected and radiated or radiated empirically, and high-grade glioma requires maximal safe resection followed by chemoradiation. CNSIDD should only be biopsied when diagnostic uncertainty requires it, and resection and radiotherapy are unnecessary, introducing unwarranted morbidity. Polygenic risk models can identify patients at the highest risk of developing glioma; we hypothesized that these polygenic models could help with differential diagnosis of indeterminate brain lesions. We also hypothesized that race would be an important contributing factor in the models. In the initial discovery and validation European (EUR) cohorts the mean probability of glioma for IDHmut non-codeleted glioma was 0.55 and 0.52, respectively, and in healthy controls was 0.19 and 0.21, respectively. To further evaluate sensitivity, we analyzed additional genotype data from 867 gliomas (764 EUR, 54 AFR, 24 AMR, 15 EAS, 10 SAS) from The Cancer Genome Atlas (TCGA). Across 764 EUR IDHmut non-codeleted glioma, the mean probability 0.53, whereas across 54 AFR and 24 AMR the mean was 0.22 and 0.32, respectively. To evaluate specificity, we analyzed 3200 TCGA patients with primary tumor types that commonly metastasize to the brain (2676 EUR, 365 AFR, 46 AMR, 95 EAS, 18 SAS), and 236 AFR healthy controls. For patients with non-CNS primary tumors, the mean probability ranged from 0.09-0.18 for EUR and 0.04-0.07 for AFR. For AFR healthy controls, the mean was 0.05. Overall, race is a significant factor for polygenic risk models. Further work entails evaluating polygenic risk models in IDHwt glioma, CNS lymphoma and CNSIDD cohorts as well as in additional races.

Abstract BACKGROUND Determination of the causation of germline single nucleotide polymorphisms (SNPs) located in non-coding regions of the genome is challenging. The genomic region of 8q24 has been identified as important in many kinds of cancer, linked to a topologically associated domain (TAD) encompassing MYC; this TAD contains a GWAS SNP (rs55705857) associated with IDH-mutant glioma. METHODS Germline genotyping data from 622 IDH-mutant glioma and 668 controls were used to fine map the rs55705857 locus by detailed haplotype analysis. Chromatin immunoprecipitation sequencing (ChIP-seq) of histone markers H3K4me1, H3K4me3, H3K27ac and H3K36me3 was performed on normal brain samples (n=8) and human glioma samples (n=11 IDH-wt and 52 IDH-mut). RNAseq from 9 normal and 83 brain tumors (n=26 IDH-wt and 55 IDH-mut) were used to assess differential gene expression. RESULTS Fine-mapping identified rs55705857 SNP as the most likely causative allele (OR=8.69; p<0.001) within 8q24 for the development of IDH-mutant glioma. At rs55705857, both H3K27ac and H3K4me1 in IDH-mutant vs IDH-wt tumors were increased 3.05- and 1.58-fold, respectively (DiffBind; p=5.81×10-7 and p=2.31×10-3). ChromHMM analysis of the marks indicated that promoter and enhancer functions were significantly increased, and the activity broadened at rs55705857 in IDH-mut gliomas compared to IDH-wt tumors and normal brain samples. This enhancement correlated with significant increased MYC expression in IDH-mut gliomas (p=3.1×10-13), as well as alterations of Myc signaling targets. Publicly available ATACseq, ChIPseq and long-range DNA interaction data demonstrated that the rs55705857 locus is open and interacts with the MYC promoter. CONCLUSIONS Fine-mapping of the 8q24 locus provided strong evidence that rs55705857 is the causative 8q24 locus associated with IDH-mut glioma. Functional experiments suggest that IDH mutation facilitates rs55705857 interaction with MYC to alter downstream MYC targets.

Abstract Establishing causal links between genetic polymorphisms and increased heritable risk of developing brain cancer is a major challenge. The non-coding single nucleotide polymorphism rs55705857 (A >G) is associated with a ~6-fold increased risk to develop IDH-mutant low-grade glioma (LGG). The rs55705857 G allele has a minor allele frequency of only ~5% in the general population but is found in ~40% of patients with IDH-mutant LGG and patients carrying risk-alleles are diagnosed on average 7-12 years earlier than those carrying non-risk A alleles. This makes rs55705857 one of the highest reported genetic associations with cancer, comparable with inherited BRCA1 gene mutations and the risk of developing breast cancer or other familial glioma genes such as NF1/2, CDKN2A or p53. To generate a LGG mouse model, we combined clonal activation of IDH1R132H with mutations of Trp53 and Atrx, which resulted in the development of LGG-like brain tumors in 25% of mice. Mutating the highly conserved, orthologous mouse rs55705857 locus to mimic the human risk allele dramatically accelerated tumor development from 463 to 172 days and increased penetrance to 75%. The resulting tumors exhibit elevated R-2-hydroxyglutarate levels, well-differentiated fibrillary neoplastic histology and metabolic rewiring, recapitulating histopathological and molecular hallmarks of human LGG. Mechanistically, we show that the rs55705857 locus resides within a brain-specific enhancer, which shows enhanced activity in IDH-mutant tumors. In addition, we found that the risk allele disrupts OCT2/4 binding, allowing increased interaction with the Myc promoter and increased Myc expression. The hyperactive chromatin status combined with the tissue specificity of this enhancer explains the cooperativity between mutant IDH and rs55705857 and why rs55705857 is associated specifically with IDH-mutant glioma, but not other cancers. Overall, we generated new LGG mouse models, which provide insights into the pathophysiology of this deadly disease and shed light into the heritable predisposition to LGG development.

AbstractBackgroundLarge-scale genome-wide association studies (GWAS) have implicated thousands of germline genetic variants in modulating individuals’ risk to various diseases, including cancer. At least 25 risk loci have been identified for low-grade gliomas (LGGs), but their molecular functions remain largely unknown.MethodsWe hypothesized that GWAS loci contain causal single nucleotide polymorphisms (SNPs) that reside in accessible open chromatin regions and modulate the expression of target genes by perturbing the binding affinity of transcription factors (TFs). We performed an integrative analysis of genomic and epigenomic data from The Cancer Genome Atlas and other public repositories to identify candidate causal SNPs within linkage disequilibrium blocks of LGG GWAS loci. We assessed their potential regulatory role via in silico TF binding sequence perturbations, convolutional neural network trained on TF binding data, and simulated annealing–based interpretation methods.ResultsWe built an interactive website (http://education.knoweng.org/alg3/) summarizing the functional footprinting of 280 variants in 25 LGG GWAS regions, providing rich information for further computational and experimental scrutiny. We identified as case studies PHLDB1 and SLC25A26 as candidate target genes of rs12803321 and rs11706832, respectively, and predicted the GWAS variant rs648044 to be the causal SNP modulating ZBTB16, a known tumor suppressor in multiple cancers. We showed that rs648044 likely perturbed the binding affinity of the TF MAFF, as supported by RNA interference and in vitro MAFF binding experiments.ConclusionsThe identified candidate (causal SNP, target gene, TF) triplets and the accompanying resource will help accelerate our understanding of the molecular mechanisms underlying genetic risk factors for gliomas.

Abstract Background: Genome-wide association studies (GWAS) revealed that 25 regions in 24 genes are associated with adult diffuse glioma development. These regions were identified by performing GWAS of glioma overall and by pathology (GBM and nonGBM). The 2016 WHO Classification of Tumors of the Central Nervous System utilizes two somatic alterations to molecularly-classify adult diffuse glioma: IDH mutation and 1p/19q codeletion. We hypothesized that germline variants may increase susceptibility to, or interact with, these somatic alterations to accelerate the development of specific molecular subtypes of glioma. We further hypothesize that germline variants associated with IDH-mutated glioma might be associated with other IDH-mutated tumors, namely, cholangiocarcinoma, acute myeloid leukemia (AML) and melanoma. Methods: We performed a GWAS by glioma molecular subtype - as defined by presence or absence of IDH somatic mutation and 1p/19q codeletion. A total of 1320 glioma cases and 1889 controls were used in the discovery set, and 799 glioma cases and 808 controls in the validation set. A meta-analysis was performed with a genome-wide p-value threshold of 5 × 10−8. GTEx data were used to perform an expression quantitative trait loci (eQTL) analysis. For germline variants that were significantly associated with IDH-mutated glioma, we evaluated pleiotropy with cholangiocarcinoma, AML and melanoma using TCGA and Mayo Biobank controls. Results: Variants in or near D2HGDH on chromosome 2 were genome-wide significant in IDH-mutated glioma (meta p-value = 2.82 × 10−10). TCGA reported that the D2HGDH region was commonly deleted in IDH-mutated gliomas that do not have 1p/19q codeletion. In TCGA data for IDH-mutated, non-codeleted glioma, we observed that the D2HGDH variant was inversely associated with tumor deletions of D2HGDH (odds ratio=0.57, p-value=0.015). The eQTL analyses demonstrated significant associations between D2HGDH germline variant and expression of D2HGDH (p=2.2 × 10−11). Further stratifying IDH-mutated glioma by 1p/19q codeletion status, one variant near FAM20C on chromosome 7 was genome-wide significant in gliomas that have IDH mutation and 1p/19q codeletion (meta p-value=9.56 × 10−9). Analyses are currently underway to evaluate pleiotropy of these IDH-mutated glioma germline variants with other IDH-mutated tumors including cholangiocarcinoma, AML and melanoma. Variants in or near GMEB2 on chromosome 20 were genome-wide significant in IDH wild-type glioma (meta p-value=2.60 × 10−10). The most significant variant in the GMEB2 region remained significant after adjustment for the known RTEL1 glioma risk variant nearby on chromosome 20 (p=0.029). Conclusions: We identified and validated novel germline variants in two genes that are associated with etiology of IDH-mutated and one gene that is associated with IDH wild-type adult diffuse glioma. Citation Format: Jeanette E. Eckel-Passow, Kristen L. Drucker, Thomas M. Kollmeyer, Matthew L. Kosel, Paul A. Decker, Annette M. Molinaro, Terri Rice, Corrine E. Praska, Lauren E. Clark, Alissa A. Caron, Alexej Abyzov, Anthony Batzler, Jun S. Song, Melike Pekmezci, Helen M. Hansen, Lucie S. McCoy, Paige M. Bracci, Joseph Wiemels, John K. Wiencke, Stephen Francis, Terence C. Burns, Caterina Giannini, Daniel H. Lachance, Margaret Wrensch, Robert B. Jenkins. Adult diffuse glioma GWAS by molecular subtype identifies variants in D2HGDH, FAM20C and GMEB2 [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1193.

Abstract Background Twenty-five germline variants are associated with adult diffuse glioma, and some of these variants have been shown to be associated with particular subtypes of glioma. We hypothesized that additional germline variants could be identified if a genome-wide association study (GWAS) were performed by molecular subtype. Methods A total of 1320 glioma cases and 1889 controls were used in the discovery set and 799 glioma cases and 808 controls in the validation set. Glioma cases were classified into molecular subtypes based on combinations of isocitrate dehydrogenase (IDH) mutation, telomerase reverse transcriptase (TERT) promoter mutation, and 1p/19q codeletion. Logistic regression was applied to the discovery and validation sets to test for associations of variants with each of the subtypes. A meta-analysis was subsequently performed using a genome-wide P-value threshold of 5 × 10−8. Results Nine variants in or near D-2-hydroxyglutarate dehydrogenase (D2HGDH) on chromosome 2 were genome-wide significant in IDH-mutated glioma (most significant was rs5839764, meta P = 2.82 × 10−10). Further stratifying by 1p/19q codeletion status, one variant in D2HGDH was genome-wide significant in IDH-mutated non-codeleted glioma (rs1106639, meta P = 4.96 × 10−8). Further stratifying by TERT mutation, one variant near FAM20C (family with sequence similarity 20, member C) on chromosome 7 was genome-wide significant in gliomas that have IDH mutation, TERT mutation, and 1p/19q codeletion (rs111976262, meta P = 9.56 × 10−9). Thirty-six variants in or near GMEB2 on chromosome 20 near regulator of telomere elongation helicase 1 (RTEL1) were genome-wide significant in IDH wild-type glioma (most significant was rs4809313, meta P = 2.60 × 10−10). Conclusions Performing a GWAS by molecular subtype identified 2 new regions and a candidate independent region near RTEL1, which were associated with specific glioma molecular subtypes.

Abstract Genome-wide association studies (GWAS) have revealed that 25 regions in 24 genes are associated with adult diffuse glioma development. These regions were identified by performing GWAS of glioma overall and GWAS by pathology (GBM and nonGBM). Subsequently, these regions have been evaluated for associations with specific molecular subtypes. The 2016 WHO Classification of Tumors of the Central Nervous System utilizes two somatic alterations to molecularly-classify adult diffuse glioma: IDH mutation and 1p/19q codeletion. TERT promoter mutation has also been shown to be associated with age at diagnosis and patient outcome. We hypothesized that germline variants may increase susceptibility to, or interact with, these somatic alterations to accelerate the development of specific molecular subtypes of glioma. To test our hypothesis, we performed a GWAS by glioma molecular subtype – as defined by presence or absence of IDH and TERT somatic mutation and 1p/19q codeletion – utilizing a two-stage design and subsequent meta analysis that included 3001 total glioma cases and 2697 total controls. Data were imputed using the Michigan Server and logistic regression was used, adjusting for age and sex. The Cancer Genome Atlas (TCGA) data were used to perform an expression quantitative trait loci (eQTL) analysis on candidate germline variants. Variants in 2q37 and 7p22 were associated with IDH-mutated glioma (meta analysis p< 5x10-8). The eQTL analyses demonstrated significant associations between 2q37 variants and expression of nearby genes as well as associations between 7p22 variants and nearby genes (p< 0.0001). In conclusion, we identified and validated novel germline variants in two genes that are associated with etiology of IDH-mutated adult diffuse glioma.