Wrensch, Margaret

ABSTRACT Background Most genomic studies of adult-type diffuse gliomas have focused on predominantly European ancestry populations, limiting the generalizability of molecular classifications and precision medicine approaches. We assembled a multi-institutional glioma cohort of diverse patients to investigate how germline ancestry, molecular subtypes, and mutational processes shape tumor biology and clinical outcomes. Methods We analyzed 1,102 adults with WHO 2021–classified diffuse gliomas (IDH-mutant, 1p/19q-codeleted oligodendroglioma; IDH-mutant astrocytoma; IDH-wildtype glioma) from seven U.S. institutions. Whole-exome sequencing (WES) of FFPE tumors identified somatic alterations and COSMIC SBS v3.2 mutational signatures. Genetic ancestry was estimated from WES using 1000 Genomes reference populations. Overall survival was assessed using Kaplan–Meier and multivariable models. Results The cohort included 66.9% European (EUR), 21.1% Admixed American/Hispanic (AMR), 10.3% Admixed African (AFR), and 1.6% Asian (AS) ancestry. Survival followed expected molecular hierarchy (median overall survival: oligodendroglioma 15.7 years, astrocytoma 10.6 years, IDH-wildtype glioma 1.9 years). Within oligodendroglioma, AMR patients showed improved survival versus EUR (HR 0.67, 95% CI 0.48–0.94; p=0.011), with similar trends across subtypes. Somatic profiling confirmed canonical subtype-defining alterations and revealed higher ATRX alterations in AFR and AMR IDH-wildtype tumors compared with EUR. ATRX alterations were associated with improved survival only in AFR (p=0.003). Mutational signature analysis identified subtype-specific signatures, including therapy-associated signatures. Chemotherapy-related signatures were more frequent in EUR and AMR than in AFR. Conclusions This ancestrally diverse glioma cohort confirms established molecular classifications and identifies ancestry-associated differences in survival, somatic alterations, and mutational processes, indicating the critical need for broad representation to inform precision neuro-oncology. Key Points A multi-institutional glioma cohort validates subtype and survival patterns across ancestries. Therapy-associated mutational signatures differ by ancestry, suggesting distinct treatment-related mutational processes. Admixed American patients show improved survival, particularly in oligodendroglioma. Importance of the Study Most genomic studies of adult-type diffuse glioma have focused on populations of predominantly European ancestry which limits the ability to examine variation in tumor biology and clinical outcomes across populations. In this study, we assembled one of the largest ancestrally diverse cohorts of molecularly characterized adult diffuse gliomas, integrating germline ancestry inference with tumor whole-exome sequencing and mutational signature analysis. We confirm that established molecular classifications and survival hierarchies remain robust across ancestry groups. However, we also identified ancestry-associated differences in survival within specific tumor subtypes, higher ATRX alteration frequencies in African American and admixed American patients with IDH-wildtype tumors, and variation in therapy-associated mutational signatures across ancestry groups. These findings highlight the importance of incorporating population differences into genomic studies of glioma and provide a resource for future multi-ancestry investigations of glioma risk, tumor evolution, and treatment response, ultimately supporting more inclusive precision neuro-oncology.

Abstract Background The CCDC26 germline variant rs55705857 is causal for development of IDH mutant (IDHmut) adult glioma. However, ~60% of IDHmut patients do not carry the rs55705857 risk allele. We aimed to identify variants associated with developing IDHmut glioma among patients that do not have the rs55705857 risk allele and to further understand development of IDHwt glioma. Methods We used three datasets that included 1216 IDHmut and 1442 IDHwt glioma patients and a case–case design to perform genome-wide association (GWAS) analyses comparing IDHmut versus IDHwt glioma. Analyses were performed overall and stratified by rs55705857 genotype and sex. Multivariable logistic regression and regression trees were used to develop models to predict IDH status using germline variants, age, and contrast enhancement on MRI. Results Three regions were identified comparing IDHmut versus IDHwt: rs55705857 (meta P-value [P] = 1.35 × 10−43), PHLDB1 (rs7125115, P = 3.46 × 10−17), and D2HGDH (rs71430382, P = 2.43 × 10−12). When analyzing only patients that do not have the rs55705857 risk allele, PHLDB1 (rs7125115, P = 1.73 × 10−13) and D2HGDH (rs71430382, P = 8.86 × 10−10) were identified. Among patients who have the rs55705857 risk allele, four variants between ROBO1 and ROBO2 (rs4680975, P = 4.65 × 10−7) increased the likelihood of having an IDHwt tumor. Tumor expression of ROBO1 was associated with rs4680975 genotype in IDHwt patients that have the rs55705857 risk allele (P = 0.034). Multivariable logistic analysis demonstrated that rs55705857, rs71430382 (D2HGDH), and age predicted IDH mutation status. Conclusions To understand the development of adult glioma, we demonstrate that D2HGDH and PHLDB1 should be prioritized for functional studies in IDHmut tumors. The ROBO1 region warrants further investigation in IDHwt tumors.

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 Acquired tumor alterations are part of the WHO diagnosis of glioma. However, germline predisposition to these alterations is unknown. GWAS in European and Asian populations identified germline variants in 30 regions, many in or near genes that also acquire alterations, suggesting that interactions between inherited and acquired genetics are involved in glioma development. To evaluate such interactions, we used 3686 glioma cases and 1889 controls. A case-control design and multinomial logistic regression were used to identify germline variants associated with development of TERT-promoter mutant gliomas, generating two interesting observations. First, TERT germline variants were associated with predisposition to IDHwt glioma that have a TERT promoter mutation, but not with IDH-mutant codeleted gliomas that have a TERT promoter mutation. Second, two novel regions were observed in this more homogeneous subgroup of TERT mutated tumors: variants in TXNDC11 (OR=2.66, p=5.7x10-8) and RUVBL2 (OR=3.1, p=5.8x10-8). High tumor expression of TXNDC11 is associated with poor prognosis in glioma and RUVBL2 interacts with TERT to affect telomerase activity. Similar analyses are being performed on additional clinically relevant acquired alterations. In parallel, case-case GWAS analyses were performed to evaluate germline variants associated with other tumor alterations. Comparing IDH-mutant versus IDHwt stratified by rs55705857 genotype, the most significant variants were PHLDB1 (p=1.2x10-13) and D2HGDH (p=2.1x10-10), highlighting the importance of these two variants independent of CCDC26 for the development of IDH-mutant glioma. As a parallel approach, we utilized recursive partitioning and regression trees (rpart) to classify IDH mutation status using the 30 known variants. The top four branches of the regression tree included CCDC26, D2HGDH, PHLDB1, and MAML2, as well as two Asian variants (RAB27A, CYP4F12). This highlights the importance of utilizing variants discovered in non-European populations in analyses. Overall, understanding these interactions are necessary to develop therapeutics, stratify clinical trials, and to develop relevant tumor models.

Abstract Gliomas are heterogeneous brain tumors with prognostically-significant subtypes and few risk factors. IDH-wildtype glioblastoma, the most common subtype, carries a dismal prognosis. Using the largest GWAS dataset with histological (11,304 cases, 304,523 controls) and molecular (3,418 cases, 8,156 controls) subtypes, we analyzed the plasma and brain proteome to uncover proteins associated with glioma susceptibility. Using protein quantitative trait loci (pQTL) for 1,776 proteins in brain tissue, we identified 21 candidate proteins, including 9 associations with high probability (>70%) of colocalization (shared causal genetic variants for glioma risk and protein abundance). We confirmed known risk genes such as EGFR (p=1.9×10-35) and D2HGDH (p=1.1×10-5). We also identified novel risk proteins, including ENPP6, a choline phosphodiesterase overexpressed in glioma, for IDH-wildtype glioblastoma (OR=0.62, p=7.6×10-5), and galectin-3, a lectin encoded by a tumor fitness gene in CRISPR screens of glioma proliferation and a clinical drug target, for glioma overall (OR=1.37, p=1.3×10-5). Analyses of the plasma proteome using pQTLs for 1,362 proteins from 3 cohorts (N=88,838) identified 9 colocalized associations in whole blood. Most risk signals were subtype-specific, such as CRELD1 (p=2.5×10-5) for triple-negative glioblastoma. We found three candidate risk proteins for IDH-mutant 1p19q-intact glioma: immunoglobulin glycoprotein BCAM (p=9.6×10-8), neurotrophin receptor SorCS2 (p=3.6×10-5) and immune checkpoint PD-1 (p=1.1×10-5). Interestingly, genetically-predicted PD-1 levels had a larger effect on IDH-mutant 1p19q-intact (OR=2.75, 95% CI=1.61-4.69) than IDH-wildtype (OR=0.73, 95% CI=0.51-1.05) tumors, suggesting germline differences in immune checkpoint mechanisms among subtypes. For all colocalized candidates, genetically-predicted protein abundance was correlated with cis-regulated transcriptional activity in brain tissue. For example, genetically-predicted SORCS2 gene expression in brain tissue was associated with plasma SorCS2 levels (Z=16.1, p=1.1×10-68). Our analysis identified candidate proteins involved in immune response and neuronal proliferation, which may enable prioritization of drug targets and provide insight into disease mechanisms for glioma.

Abstract Background Polygenic risk scores (PRS) aggregate the contribution of many risk variants to provide a personalized genetic susceptibility profile. Since sample sizes of glioma genome-wide association studies (GWAS) remain modest, there is a need to efficiently capture genetic risk using available data. Methods We applied a method based on continuous shrinkage priors (PRS-CS) to model the joint effects of over 1 million common variants on disease risk and compared this to an approach (PRS-CT) that only selects a limited set of independent variants that reach genome-wide significance (P < 5 × 10–8). PRS models were trained using GWAS stratified by histological (10 346 cases and 14 687 controls) and molecular subtype (2632 cases and 2445 controls), and validated in 2 independent cohorts. Results PRS-CS was generally more predictive than PRS-CT with a median increase in explained variance (R2) of 24% (interquartile range = 11–30%) across glioma subtypes. Improvements were pronounced for glioblastoma (GBM), with PRS-CS yielding larger odds ratios (OR) per standard deviation (SD) (OR = 1.93, P = 2.0 × 10–54 vs. OR = 1.83, P = 9.4 × 10–50) and higher explained variance (R2 = 2.82% vs. R2 = 2.56%). Individuals in the 80th percentile of the PRS-CS distribution had a significantly higher risk of GBM (0.107%) at age 60 compared to those with average PRS (0.046%, P = 2.4 × 10–12). Lifetime absolute risk reached 1.18% for glioma and 0.76% for IDH wildtype tumors for individuals in the 95th PRS percentile. PRS-CS augmented the classification of IDH mutation status in cases when added to demographic factors (AUC = 0.839 vs. AUC = 0.895, PΔAUC = 6.8 × 10–9). Conclusions Genome-wide PRS has the potential to enhance the detection of high-risk individuals and help distinguish between prognostic glioma subtypes.

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

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.