Analysis of 60 reported glioma risk SNPs replicates published GWAS findings but fails to replicate associations from published candidate-gene studies.

Genomewide association studies (GWAS) and candidate-gene studies have implicated single-nucleotide polymorphisms (SNPs) in at least 45 different genes as putative glioma risk factors. Attempts to validate these associations have yielded variable results and few genetic risk factors have been consistently replicated. We conducted a case-control study of Caucasian glioma cases and controls from the University of California San Francisco (810 cases, 512 controls) and the Mayo Clinic (852 cases, 789 controls) in an attempt to replicate previously reported genetic risk factors for glioma. Sixty SNPs selected from the literature (eight from GWAS and 52 from candidate-gene studies) were successfully genotyped on an Illumina custom genotyping panel. Eight SNPs in/near seven different genes (TERT, EGFR, CCDC26, CDKN2A, PHLDB1, RTEL1, TP53) were significantly associated with glioma risk in the combined dataset (P < 0.05), with all associations in the same direction as in previous reports. Several SNP associations showed considerable differences across histologic subtype. All eight successfully replicated associations were first identified by GWAS, although none of the putative risk SNPs from candidate-gene studies was associated in the full case-control sample (all P values > 0.05). Although several confirmed associations are located near genes long known to be involved in gliomagenesis (e.g., EGFR, CDKN2A, TP53), these associations were first discovered by the GWAS approach and are in noncoding regions. These results highlight that the deficiencies of the candidate-gene approach lay in selecting both appropriate genes and relevant SNPs within these genes.

Inherited variation in immune genes and pathways and glioblastoma risk.

To determine whether inherited variations in immune function single-nucleotide polymorphisms (SNPs), genes or pathways affect glioblastoma risk, we analyzed data from recent genome-wide association studies in conjunction with predefined immune function genes and pathways. Gene and pathway analyses were conducted on two independent data sets using 6629 SNPs in 911 genes on 17 immune pathways from 525 glioblastoma cases and 602 controls from the University of California, San Francisco (UCSF) and a subset of 6029 SNPs in 893 genes from 531 cases and 1782 controls from MD Anderson (MDA). To further assess consistency of SNP-level associations, we also compared data from the UK (266 cases and 2482 controls) and the Mayo Clinic (114 cases and 111 controls). Although three correlated epidermal growth factor receptor (EGFR) SNPs were consistently associated with glioblastoma in all four data sets (Mantel-Haenzel P values = 1 × 10⁻5 to 4 × 10⁻³), independent replication is required as genome-wide significance was not attained. In gene-level analyses, eight immune function genes were significantly (minP < 0.05) associated with glioblastoma; the IL-2RA (CD25) cytokine gene had the smallest minP values in both UCSF (minP = 0.01) and MDA (minP = 0.001) data sets. The IL-2RA receptor is found on the surface of regulatory T cells potentially contributing to immunosuppression characteristic of the glioblastoma microenvironment. In pathway correlation analyses, cytokine signaling and adhesion-extravasation-migration pathways showed similar associations with glioblastoma risk in both MDA and UCSF data sets. Our findings represent the first systematic description of immune genes and pathways that characterize glioblastoma risk.

Inherited variant on chromosome 11q23 increases susceptibility to IDH-mutated but not IDH-normal gliomas regardless of grade or histology.

INTRODUCTION: Recent discoveries of inherited glioma risk loci and acquired IDH mutations are providing new insights into glioma etiology. IDH mutations are common in lower grade gliomas and secondary glioblastomas and uncommon in primary glioblastomas. Because the inherited variant in 11q23 has been associated with risk of lower grade glioma and not with glioblastomas, we hypothesized that this variant increases susceptibility to IDH-mutated gliomas, but not to IDH-wild-type gliomas. METHODS: We tested this hypothesis in patients with glioma and controls from the San Francisco Adult Glioma Study, the Mayo Clinic, and Illumina controls (1102 total patients, 5299 total controls). Case-control additive associations of 11q23 risk alleles (rs498872, T allele) were calculated using logistic regression, stratified by tumor IDH status (mutated or wild-type) and by histology and grade. We also adjusted for the recently discovered 8q24 glioma risk locus rs55705857 G allele. RESULTS: The 11q23 glioma risk locus was associated with increased risk of IDH-mutated gliomas of all histologies and grades (odds ratio [OR] = 1.50; 95% confidence interval [CI] = 1.29-1.74; P = 1.3X10(-7)) but not with IDH-wild-type gliomas of any histology or grade (OR = 0.91; 95% CI = 0.81-1.03; P = 0.14). The associations were independent of the rs55705857 G allele. CONCLUSION: A variant at the 11q23 locus increases risk for IDH-mutated but not IDH-wild-type gliomas, regardless of grade or histology.

Genetic variants in telomerase-related genes are associated with an older age at diagnosis in glioma patients: evidence for distinct pathways of gliomagenesis.

BACKGROUND: Genome-wide association studies have implicated single nucleotide polymorphisms (SNPs) in 7 genes as glioma risk factors, including 2 (TERT, RTEL1) involved in telomerase structure/function. We examined associations of these 7 established glioma risk loci with age at diagnosis among patients with glioma. METHODS: SNP genotype data were available for 2286 Caucasian glioma patients from the University of California, San Francisco (n = 1434) and the Mayo Clinic (n = 852). Regression analyses were performed to test for associations between "number of risk alleles" and "age at diagnosis," adjusted for sex and study site and stratified by tumor grade/histology where appropriate. RESULTS: Four SNPs were significantly associated with age at diagnosis. Carrying a greater number of risk alleles at rs55705857 (CCDC26) and at rs498872 (PHLDB1) was associated with younger age at diagnosis (P = 1.4 × 10(-22) and P = 9.5 × 10(-7), respectively). These SNPs are stronger risk factors for oligodendroglial tumors, which tend to occur in younger patients, and their association with age at diagnosis varied across tumor subtypes. In contrast, carrying more risk alleles at rs2736100 (TERT) and at rs6010620 (RTEL1) was associated with older age at diagnosis (P = 6.2 × 10(-4) and P = 2.5 × 10(-4), respectively). These SNPs are risk factors for all glioma grades/histologies, and their association with age at diagnosis was consistent across tumor subgroups. CONCLUSIONS: Carrying a greater number of risk alleles might be expected to decrease age at diagnosis. However, glioma susceptibility conferred by variation in telomerase-related genes did not follow this pattern. This supports the hypothesis that telomerase-related mechanisms of telomere maintenance are more associated with gliomas that develop later in life than those utilizing telomerase-independent mechanisms (ie, alternative lengthening of telomeres).

SSBP2 variants are associated with survival in glioblastoma patients.

PURPOSE: Glioblastoma is a devastating, incurable disease with few known prognostic factors. Here, we present the first genome-wide survival and validation study for glioblastoma. EXPERIMENTAL DESIGN: Cox regressions for survival with 314,635 inherited autosomal single-nucleotide polymorphisms (SNP) among 315 San Francisco Adult Glioma Study patients for discovery and three independent validation data sets [87 Mayo Clinic, 232 glioma patients recruited from several medical centers in Southeastern United States (GliomaSE), and 115 The Cancer Genome Atlas patients] were used to identify SNPs associated with overall survival for Caucasian glioblastoma patients treated with the current standard of care, resection, radiation, and temozolomide (total n = 749). Tumor expression of the gene that contained the identified prognostic SNP was examined in three separate data sets (total n = 619). Genotype imputation was used to estimate hazard ratios (HR) for SNPs that had not been directly genotyped. RESULTS: From the discovery and validation analyses, we identified a variant in single-stranded DNA-binding protein 2 (SSBP2) on 5q14.1 associated with overall survival in combined analyses (HR, 1.64; P = 1.3 × 10(-6)). Expression of SSBP2 in tumors from three independent data sets also was significantly related to patient survival (P = 5.3 × 10(-4)). Using genotype imputation, the SSBP2 SNP rs17296479 had the strongest statistically significant genome-wide association with poorer overall patient survival (HR, 1.79; 95% CI, 1.45-2.22; P = 1.0 × 10(-7)). CONCLUSION: The minor allele of SSBP2 SNP rs17296479 and the increased tumor expression of SSBP2 were statistically significantly associated with poorer overall survival among glioblastoma patients. With further confirmation, previously unrecognized inherited variations influencing survival may warrant inclusion in clinical trials to improve randomization. Unaccounted for genetic influence on survival could produce unwanted bias in such studies.

A low-frequency variant at 8q24.21 is strongly associated with risk of oligodendroglial tumors and astrocytomas with IDH1 or IDH2 mutation.

Variants at 8q24.21 have been shown to be associated with glioma development. By means of tag SNP genotyping and imputation, pooled next-generation sequencing using long-range PCR and subsequent validation SNP genotyping, we identified seven low-frequency SNPs at 8q24.21 that were strongly associated with glioma risk (P=1×10(-25) to 1×10(-14)). The most strongly associated SNP, rs55705857, remained highly significant after individual adjustment for the other top six SNPs and two previously published SNPs. After stratifying by histological and tumor genetic subtype, the most significant associations of rs55705857 were with oligodendroglial tumors and gliomas with mutant IDH1 or IDH2 (odds ratio (OR)=5.1, P=1.1×10(-31) and OR=4.8, P=6.6×10(-22), respectively). Strong associations were observed for astrocytomas with mutated IDH1 or IDH2 (grades 2-4) (OR=5.16-6.66, P=4.7×10(-12) to 2.2×10(-8)) but not for astrocytomas with wild-type IDH1 and IDH2 (smallest P=0.26). The conserved sequence block that includes rs55705857 is consistently modeled as a microRNA.

Epigenetic profiles distinguish pleural mesothelioma from normal pleura and predict lung asbestos burden and clinical outcome.

Mechanisms of action of nonmutagenic carcinogens such as asbestos remain poorly characterized. As pleural mesothelioma is known to have limited numbers of genetic mutations, we aimed to characterize the relationships among gene-locus-specific methylation alterations, disease status, asbestos burden, and survival in this rapidly fatal asbestos-associated tumor. Methylation of 1505 CpG loci associated with 803 cancer-related genes were studied in 158 pleural mesotheliomas and 18 normal pleura. After false-discovery rate correction, 969 CpG loci were independently associated with disease status (Q < 0.05). Classifying samples based on CpG methylation profile with a mixture model approach, methylation classes discriminated tumor from normal pleura (permutation P < 0.0001). In a random forests classification, the overall misclassification error rate was 3.4%, with <1% (n = 1) of tumors misclassified as normal (P < 0.0001). Among tumors, methylation class membership was significantly associated with lung tissue asbestos body burden (P < 0.03), and significantly predicted survival (likelihood ratio P < 0.01). Consistent with prior work, asbestos burden was associated with an increased risk of death (hazard ratio, 1.4; 95% confidence interval, 1.1-1.8). Our results have shown that methylation profiles powerfully differentiate diseased pleura from nontumor pleura and that asbestos burden and methylation profiles are independent predictors of mesothelioma patient survival. We have added to the growing body of evidence that cellular epigenetic dysregulation is a critical mode of action for asbestos in the induction of pleural mesothelioma. Importantly, these findings hold great promise for using epigenetic profiling in the diagnosis and prognosis of human cancers.