Wrensch, Margaret

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 BACKGROUND The prognostic importance of maximal resection of contrast enhancing and non-contrast enhancing disease has been established. Nonetheless, glioblastomas exist within the framework of complex neural circuitry serving cognition, movement, and behavior consequential leading to neurological impairments. The prognostic importance of neurological impairments on survival remains poorly understood. METHODS This is a retrospective, single cohort study from UCSF including 316 eligible patients diagnosed over 20 years with 9.6 years of follow-up. All patients underwent surgical resection for newly diagnosed glioblastoma for whom survival, molecular, preoperative and postoperative MRI images, and clinical data were available. All patients had chemoradiation treated IDH-wild-type glioblastoma with available preoperative and 1-month post-surgical resection neurological outcomes. We employed survival models and recursive partitioning (RPA) to investigate multivariate relationships of overall survival (OS). RESULTS Preoperative neurological impairments were present in 75.6% (n= 239) and new post resection impairments were identified in 37.3% (n=117). Univariate analysis confirmed that new postoperative cognitive impairment [HR 7.91, 95% CI 2.47-25.33] and hemiplegia [HR 3.38, 95% CI 0.83-13.67] (not hemiparesis) impact OS. Risk stratified grouping by RPA demonstrated that gross total resection of contrast enhancing tumor in patients with no new postoperative neurological impairments confers the longest OS (median OS 27.1 months 95%CI 21.5-33.7). Patients with any residual tumor volume after surgery but no new neurological deficits experience a similar survival to younger patients (under 65) with 1 or more new postoperative neurological deficits (median OS 16.6 months 95%CI 15.2-19.2). Shortest OS is identified in patients with any volume of residual tumor plus 1 or more new postoperative neurological deficits and age over 65 (median OS 11.4 months, 95%CI 9.3-13.5). CONCLUSIONS This study confirms that new postoperative neurological impairments impact overall survival in patients with chemoradiation treated IDH-wild-type glioblastoma.

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.

Abstract Background The etiological basis of glioma is poorly understood. We have used genetic markers in a Mendelian Randomisation (MR) framework to examine if lifestyle, cardiometabolic and inflammatory factors influence the risk of glioma. This methodology reduces bias from confounding and is not affected by reverse causation. Methods We identified genetic instruments for 37 potentially modifiable risk factors and evaluated their association with glioma risk using data from a genome-wide association study of 12,488 glioma patients and 18,169 controls. We used the estimated odds ratio of glioma associated with each of the genetically defined traits to infer evidence for a causal relationship with the following exposures: lifestyle and dietary factors (height, plasma IGF-1, blood carnitine, blood methionine, blood selenium, blood zinc, circulating adiponectin, circulating carotenoids, iron status, serum calcium, vitamin [A1, B12, B6, E and 25-hydroxyvitamin D], fatty acids levels [mono-unsaturated, omega-3 and omega-6] and circulating fetuin-A); cardiometabolic factors (birth weight, HDL cholesterol, LDL cholesterol, total cholesterol, total triglycerides, basal metabolic rate, body fat percentage, body mass index, fasting glucose, fasting proinsulin, HbA1C levels, diastolic and systolic blood pressure, waist circumference, waist-to-hip ratio) were included; inflammatory factors (C-reactive protein (CRP), plasma IL-6 sRa and serum IgE). Results After correction for the testing of multiple potential risk factors and excluding associations driven by one single nucleotide polymorphism (SNP) no significant association with glioma risk was observed (i.e. PCorrected > 0.05). Conclusions This study did not provide evidence supporting any of the 37 factors examined as having a significant influence on glioma risk.

Abstract Little is known about the interaction between constitutive and somatic genetic variation in carcinogenesis. Genome-wide association studies (GWAS) of childhood acute lymphoblastic leukemia (ALL) have identified 6 genetic associations, including the single nucleotide polymorphism (SNP) rs3731217, located between the CDKN2A and CDKN2B genes. Somatic deletions of this chromosome 9p21.3 region are common in cancer, and occur in ~30% of childhood ALL. We replicated the rs3731217 association with childhood ALL in Hispanics (P = 0.021) via a GWAS of 323 childhood ALL cases and 454 controls, using the Illumina OmniExpress SNP array. SNP imputation analysis across chromosome 9p21.3 revealed a stronger association at SNP rs3731249 (P = 6.4x10-4; OR = 2.77), a low-frequency missense variant in exon 2 (p.Ala148Thr) of CDKN2A (MAFcontrols = 2.1%). These associations were replicated in an additional 378 ALL cases and 536 controls using Taqman assays, where the missense variant again had a greater magnitude of effect (Prs3731249 = 6.5x10-3; OR rs3731249 = 2.00 versus Prs3731217 = 0.029; OR rs3731217= 1.38). We assessed the effect of rs3731249 on risk of other cancers with known chr9p21.3 GWAS associations, but observed no association between this variant and either adult glioblastoma or melanoma risk, suggesting an ALL-specific role for rs3731249. We hypothesized that cases heterozygous for rs3731249 may show allelic imbalance, with preferential loss of the wildtype (WT) allele when hemizygous chr9p21.3 deletions occur in leukemic cells. To test this, we developed a novel method termed “SMART-ddPCR” (Somatic Mutation Allelic Ratio Test using Droplet Digital PCR). Concentration of risk and WT alleles was measured by ddPCR using the rs3731249 Taqman assay in diagnostic bone marrow (i.e. tumor) DNA from 35 known heterozygotes. Allelic imbalance due to somatic CDKN2A alterations was identified in 17 cases, of which 14 had a higher risk:WT allele ratio (P = 0.006, 1-sided binomial significance test), thus demonstrating preferential loss of WT rs3731249 and retention of the missense allele. Subclonal heterogeneity in CDKN2A copy number was also evident, suggesting that CDKN2A alterations are secondary events in leukemogenesis. Preferential allelic imbalance towards the rs3731249 risk allele in ALL tumor DNA provides strong evidence of an important role for this missense variant in leukemogenesis. We show, for the first time, a direct relationship between constitutive and somatic genetic variants underlying ALL development, suggesting that such constitutive variants can provide an additional substrate for selection during tumor evolution. Disclosures No relevant conflicts of interest to declare.

Abstract BACKGROUND: Approximately 5% of gliomas occur in individuals with a family history of glioma, and first-degree relatives of brain tumor cases have a two-fold increase in risk of brain tumor. Recent somatic characterization has shown that tumors from familial cases are indistinguishable from sporadic cases, suggesting that familial cases may arise through similar mechanisms of gliomagenesis, and therefore may be associated with common variants as well as rare mutations. In this analysis, we assessed whether previously identified common risk variants are associated with familial glioma. METHODS: Data were obtained from the Glioma International Case Control (GICC) and Gliogene studies for 448 cases with reported family history, 4,405 cases without reported family history, and 3,288 controls. We assessed 25 risk loci previously identified by glioma GWAS, and odds ratios (OR) and 95% confidence intervals (95%CI) were calculated using an additive genetic logistic regression model adjusted for age, sex, and the first two principal components for familial cases versus unaffected controls, and non-familial cases versus controls. Results were considered significant at p&lt;0.002 (Bonferroni correction for 25 tests). RESULTS: Significant associations were detected at 5/25 loci, including TERT, EGFR, CCDC26, CDKN2B, and RTEL1. The strongest association was at rs55705857 (CCDC26, OR=2.7, p=7.49x10-17). For GBM (222 familial cases), significant associations were detected at 6/26 loci (TERT, EGFR, CDKN2B, TP53 and RTEL1), while in non-GBM (205 familial cases) significant associations were detected at 3/25 loci (LRIG1, CCDC26, PHLDB1). These SNPs were further examined using a case-only approach comparing familial to non-familial cases, and there was no significant difference in allele frequencies by family history status. There was a strong correlation between log(OR) for familial cases only versus non-familial cases (adjusted R2=0.88). CONCLUSIONS: In this analysis we identified a significant association between familial glioma and five common risk loci previously identified by glioma GWAS. This provides further evidence of shared pathways of genetic risk and gliomagenesis between familial and non-familial glioma. Further exploration is necessary to determine the overall contribution of common genetic variation to risk of familial glioma. Citation Format: Quinn T. Ostrom, Georgina Armstrong, Christopher I. Amos, Jonine L. Bernstein, Elizabeth B. Claus, Jeanette E. Eckel-Passow, Dora Il'yasova, Christoffer Johansen, Daniel H. Lachance, Rose K. Lai, Ryan T. Merrell, Sara H. Olson, Joellen H. Schildkraut, Sanjay S. Shete, Richard S. Houlston, Robert B. Jenkins, Margaret R. Wrensch, Beatrice Melin, Jill S. Barnholtz-Sloan, Melissa L. Bondy. Previously identified common glioma risk SNPs are associated with familial glioma [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 4173.

Glioma incidence is highest in non‐Hispanic Whites, and to date, glioma genome‐wide association studies (GWAS) to date have only included European ancestry (EA) populations. African Americans and Hispanics in the US have varying proportions of EA, African (AA) and Native American ancestries (NAA). It is unknown if identified GWAS loci or increased EA is associated with increased glioma risk. We assessed whether EA was associated with glioma in African Americans and Hispanics. Data were obtained for 832 cases and 675 controls from the Glioma International Case–Control Study and GliomaSE Case–Control Study previously estimated to have <80% EA, or self‐identify as non‐White. We estimated global and local ancestry using fastStructure and RFMix, respectively, using 1,000 genomes project reference populations. Within groups with ≥40% AA (AFR ≥0.4 ), and ≥15% NAA (AMR ≥0.15 ), genome‐wide association between local EA and glioma was evaluated using logistic regression conditioned on global EA for all gliomas. We identified two regions (7q21.11, p = 6.36 × 10 −4 ; 11p11.12, p = 7.0 × 10 −4 ) associated with increased EA, and one associated with decreased EA (20p12.13, p = 0.0026) in AFR ≥0.4 . In addition, we identified a peak at rs1620291 ( p = 4.36 × 10 −6 ) in 7q21.3. Among AMR ≥0.15 , we found an association between increased EA in one region (12q24.21, p = 8.38 × 10 −4 ), and decreased EA in two regions (8q24.21, p = 0. 0010; 20q13.33, p = 6.36 × 10 −4 ). No other significant associations were identified. This analysis identified an association between glioma and two regions previously identified in EA populations (8q24.21, 20q13.33) and four novel regions (7q21.11, 11p11.12, 12q24.21 and 20p12.13). The identifications of novel association with EA suggest regions to target for future genetic association studies.