Guerra, Geno A

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