Hansen, Helen M.

Abstract Background Temozolomide (TMZ) treatment has demonstrated a variable impact on glioma prognosis. This study examines associations of survival with DNA repair gene germline polymorphisms among glioma patients who did and did not have TMZ treatment. Identifying genetic markers which sensitize tumor cells to TMZ could personalize therapy and improve outcomes. Methods We evaluated TMZ-related survival associations of pathogenic germline SNPs and genetically predicted transcript levels within 34 DNA repair genes among 1504 glioma patients from the University of California San Francisco Adult Glioma Study (UCSF AGS) and Mayo Clinic, whose diagnoses spanned pre- and post-TMZ eras within the major known glioma prognostic molecular subtypes. Results Among those who received TMZ, 5 SNPs were associated with overall survival, but not in those who did not receive TMZ. Only rs2308321-G, in MGMT, was associated with decreased survival (hazard ratio = 1.21, P = .019) for all glioma subtypes. Rs73191162-T (near UNG), rs13076508-C (near PARP3), rs7840433-A (near NEIL2), and rs3130618-A (near MSH5) were only associated with survival and TMZ treatment for certain subtypes, suggesting subtype-specific germline chemosensitization. Genetically predicted elevated expression of PNKP, compared to normal brain expression, was associated with markedly poor survival in TMZ-treated patients with isocitrate dehydrogenase (IDH)-mutant and 1p/19q non-codeleted gliomas (P = .015), with a median difference of over 70 months in overall survival times. Similarly, NEIL2 and TDG expressions were associated with altered TMZ-related survival only among certain subtypes. Conclusions Functional germline alterations within DNA repair genes were associated with TMZ sensitivity, measured by overall survival among adults with glioma. These variants should be evaluated in prospective analyses and functional studies.

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

Glioblastoma (GBM) is the most common malignant brain tumor in the United States. Incidence of GBM increases with age, and younger age‐at‐diagnosis is significantly associated with improved prognosis. While the relationship between candidate GBM risk SNPs and age‐at‐diagnosis has been explored, genome‐wide association studies (GWAS) have not previously been stratified by age. Potential age‐specific genetic effects were assessed in autosomal SNPs for GBM patients using data from four previous GWAS. Using age distribution tertiles (18–53, 54–64, 65+) datasets were analyzed using age‐stratified logistic regression to generate p values, odds ratios (OR), and 95% confidence intervals (95%CI), and then combined using meta‐analysis. There were 4,512 total GBM cases, and 10,582 controls used for analysis. Significant associations were detected at two previously identified SNPs in 7p11.2 (rs723527 [p54–63 = 1.50x10−9, OR54–63 = 1.28, 95%CI54–63 = 1.18–1.39; p64+ = 2.14x10−11, OR64+ = 1.32, 95%CI64+ = 1.21–1.43] and rs11979158 [p54–63 = 6.13x10−8, OR54–63 = 1.35, 95%CI54–63 = 1.21–1.50; p64+ = 2.18x10−10, OR64+ = 1.42, 95%CI64+ = 1.27–1.58]) but only in persons >54. There was also a significant association at the previously identified lower grade glioma (LGG) risk locus at 8q24.21 (rs55705857) in persons ages 18–53 (p18–53 = 9.30 × 10−11, OR18–53 = 1.76, 95%CI18–53 = 1.49–2.10). Within The Cancer Genome Atlas (TCGA) there was higher prevalence of ‘LGG’‐like tumor characteristics in GBM samples in those 18–53, with IDH1/2 mutation frequency of 15%, as compared to 2.1% [54–63] and 0.8% [64+] (p = 0.0005). Age‐specific differences in cancer susceptibility can provide important clues to etiology. The association of a SNP known to confer risk for IDH1/2 mutant glioma and higher prevalence of IDH1/2 mutation within younger individuals 18–53 suggests that more younger individuals may present initially with ‘secondary glioblastoma.’