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−3), 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

Epigenetic biomarkers of T-cells in human glioma

Immune factors are thought to influence glioma risk and outcomes, but immune profiling studies to further our understanding of the immune response are limited by current immunodiagnostic methods. We developed a new assay to capture glioma immune biology based on quantitative methylation specific PCR (qMSP) of two T-cell genes (CD3Z: T-cells, and FOXP3: Tregs). Flow cytometry of T-cells correlated well with the CD3Z demethylation assay (r = 0.93; p < 2.2 × 10(−16)), demonstrating the validity of the assay. Furthermore, there was a high correlation between qMSP and immunohistochemistry (IHC) in quantifying tumor infiltrating T-cells (r = 0.85; p = 3.4 × 10(−11)). Applying our qMSP methods to archival whole blood from 65 glioblastoma multiforme (GBM) cases and 94 non-diseased controls, GBM cases had highly statistically significantly lower T-cells (p = 1.7 × 10(−9)) as well as Tregs (p = 5.2 × 10(−11)) and a modestly lower ratio of Tregs/T-cells (p = 0.024). Applying the methods to 120 excised glioma tumors, we observed that tumor infiltrating CD3+ T-cells were positively correlated with glioma tumor grade (p = 5.7 × 10(−7)), and that Tregs were enriched in tumors compared with peripheral blood indicating active chemoattraction of suppressive Tregs into the tumor compartment. Poorer patient survival was correlated with higher levels of tumor infiltrating T-cells (p = 0.01) and Tregs (p = 0.04). DNA methylation based immunodiagnostics represent a new generation of powerful laboratory tools offering many advantages over conventional methods that will facilitate large clinical epidemiologic studies and capitalize on stored archival blood and tissue banks

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⁻⁵ 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 &lt; 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

SSBP2 Variants Are Associated with Survival in Glioblastoma Patients

PurposeGlioblastoma is a devastating, incurable disease with few known prognostic factors. Here, we present the first genome-wide survival and validation study for glioblastoma.Experimental designCox 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.ResultsFrom 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)).ConclusionThe 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