Pre-surgery immune profiles of adult glioma patients

INTRODUCTION

Although immunosuppression is a known characteristic of glioma, no previous large studies have reported peripheral blood immune cell profiles prior to patient surgery and chemoradiation. This report describes blood immune cell characteristics and associated variables prior to surgery among typical glioma patients seen at a large University practice.

METHODS

We analyzed pre-surgery blood samples from 139 glioma patients diagnosed with a new or recurrent grade II/III glioma (LrGG, n = 64) or new glioblastoma (GBM, n = 75) and 454 control participants without glioma. Relative cell fractions of CD4, CD8, B-cells, Natural Killer cells, monocytes, and neutrophils, were estimated via a validated deconvolution algorithm from blood DNA methylation measures from Illumina EPIC arrays.

RESULTS

Dexamethasone use at time of blood draw varied by glioma type being highest among patients with IDH wild-type (wt) GBM (75%) and lowest for those with oligodendroglioma (14%). Compared to controls, glioma patients showed statistically significant lower cell fractions for all immune cell subsets except for neutrophils which were higher (all p-values < 0.001), in part because of the higher prevalence of dexamethasone use at time of blood draw for IDHwt GBM. Patients who were taking dexamethasone were more likely to have a low CD4 count (< 200, < 500), increased neutrophils, low absolute lymphocyte counts, higher total cell count and higher NLR.

CONCLUSION

We show that pre-surgery blood immune profiles vary by glioma subtype, age, and more critically, by use of dexamethasone. Our results highlight the importance of considering dexamethasone exposures in all studies of immune profiles and of obtaining immune measures prior to use of dexamethasone, if possible.

Cytokine Levels at Birth in Children Who Developed Acute Lymphoblastic Leukemia

BACKGROUND

Prenatal immune development may play an important role in the etiology of childhood acute lymphoblastic leukemia (ALL).

METHODS

Seven cytokines, IL1β, IL4, IL6, IL8, GM-CSF, TNFα, and VEGF, were analyzed in blood spots collected at birth from 1,020 ALL cases and 1,003 controls participating in the California Childhood Leukemia Study. ORs and 95% confidence intervals (95% CI) associated with an interquartile range increment in cytokine levels were calculated using logistic regression, adjusting for sociodemographic and birth characteristics.

RESULTS

We found that patients with ALL were born with higher levels of a group of correlated cytokines than controls [IL1β: OR of 1.18 (95% confidence interval [CI], 1.03-1.35); IL8: 1.19 (1.03-1.38); TNFα: 1.15 (1.01-1.30); VEGF: 1.16 (1.01-1.33)], especially among children of Latina mothers (ORs from 1.31 to 1.40) and for ALL with high hyperdiploidy (ORs as high as 1.27). We found that neonatal cytokine levels were correlated with neonatal levels of endogenous metabolites which had been previously associated with ALL risk; however, there was no evidence that the cytokines were mediating the relationship between these metabolites and ALL risk.

CONCLUSIONS

We posit that children born with altered cytokine levels are set on a trajectory towards an increased risk for subsequent aberrant immune reactions that can initiate ALL.

IMPACT

This is the first study to evaluate the interplay between levels of immunomodulatory cytokines at birth, prenatal exposures, and the risk of childhood ALL.

Adult diffuse glioma GWAS by molecular subtype identifies variants in D2HGDH and FAM20C

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 1193: Adult diffuse glioma GWAS by molecular subtype identifies variants in D2HGDH, FAM20C and GMEB2

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.

Using germline variants to estimate glioma and subtype risks

BACKGROUND

Twenty-five single nucleotide polymorphisms (SNPs) are associated with adult diffuse glioma risk. We hypothesized that the inclusion of these 25 SNPs with age at diagnosis and sex could estimate risk of glioma as well as identify glioma subtypes.

METHODS

Case-control design and multinomial logistic regression were used to develop models to estimate the risk of glioma development while accounting for histologic and molecular subtypes. Case-case design and logistic regression were used to develop models to predict isocitrate dehydrogenase (IDH) mutation status. A total of 1273 glioma cases and 443 controls from Mayo Clinic were used in the discovery set, and 852 glioma cases and 231 controls from UCSF were used in the validation set. All samples were genotyped using a custom Illumina OncoArray.

RESULTS

Patients in the highest 5% of the risk score had more than a 14-fold increase in relative risk of developing an IDH mutant glioma. Large differences in lifetime absolute risk were observed at the extremes of the risk score percentile. For both IDH mutant 1p/19q non-codeleted glioma and IDH mutant 1p/19q codeleted glioma, the lifetime risk increased from almost null to 2.3% and almost null to 1.7%, respectively. The SNP-based model that predicted IDH mutation status had a validation concordance index of 0.85.

CONCLUSIONS

These results suggest that germline genotyping can provide new tools for the initial management of newly discovered brain lesions. Given the low lifetime risk of glioma, risk scores will not be useful for population screening; however, they may be useful in certain clinically defined high-risk groups.

Sex-specific gene and pathway modeling of inherited glioma risk

Background

To date, genome-wide association studies (GWAS) have identified 25 risk variants for glioma, explaining 30% of heritable risk. Most histologies occur with significantly higher incidence in males, and this difference is not explained by currently known risk factors. A previous GWAS identified sex-specific glioma risk variants, and this analysis aims to further elucidate risk variation by sex using gene- and pathway-based approaches.

Methods

Results from the Glioma International Case-Control Study were used as a testing set, and results from 3 GWAS were combined via meta-analysis and used as a validation set. Using summary statistics for nominally significant autosomal SNPs (P < 0.01 in a previous meta-analysis) and nominally significant X-chromosome SNPs (P < 0.01), 3 algorithms (Pascal, BimBam, and GATES) were used to generate gene scores, and Pascal was used to generate pathway scores. Results were considered statistically significant in the discovery set when P < 3.3 × 10-6 and in the validation set when P < 0.001 in 2 of 3 algorithms.

Results

Twenty-five genes within 5 regions and 19 genes within 6 regions reached statistical significance in at least 2 of 3 algorithms in males and females, respectively. EGFR was significantly associated with all glioma and glioblastoma in males only and a female-specific association in TERT, all of which remained nominally significant after conditioning on known risk loci. There were nominal associations with the BioCarta telomeres pathway in both males and females.

Conclusions

These results provide additional evidence that there may be differences by sex in genetic risk for glioma. Additional analyses may further elucidate the biological processes through which this risk is conferred.

Age-specific genome-wide association study in glioblastoma identifies increased proportion of 'lower grade glioma'-like features associated with younger age

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

Sex-specific glioma genome-wide association study identifies new risk locus at 3p21.31 in females, and finds sex-differences in risk at 8q24.21

Incidence of glioma is approximately 50% higher in males. Previous analyses have examined exposures related to sex hormones in women as potential protective factors for these tumors, with inconsistent results. Previous glioma genome-wide association studies (GWAS) have not stratified by sex. Potential sex-specific genetic effects were assessed in autosomal SNPs and sex chromosome variants for all glioma, GBM and non-GBM patients using data from four previous glioma GWAS. Datasets were analyzed using sex-stratified logistic regression models and combined using meta-analysis. There were 4,831 male cases, 5,216 male controls, 3,206 female cases and 5,470 female controls. A significant association was detected at rs11979158 (7p11.2) in males only. Association at rs55705857 (8q24.21) was stronger in females than in males. A large region on 3p21.31 was identified with significant association in females only. The identified differences in effect of risk variants do not fully explain the observed incidence difference in glioma by sex.

Adult infiltrating gliomas with WHO 2016 integrated diagnosis: additional prognostic roles of ATRX and TERT

The "integrated diagnosis" for infiltrating gliomas in the 2016 revised World Health Organization (WHO) classification of tumors of the central nervous system requires assessment of the tumor for IDH mutations and 1p/19q codeletion. Since TERT promoter mutations and ATRX alterations have been shown to be associated with prognosis, we analyzed whether these tumor markers provide additional prognostic information within each of the five WHO 2016 categories. We used data for 1206 patients from the UCSF Adult Glioma Study, the Mayo Clinic and The Cancer Genome Atlas (TCGA) with infiltrative glioma, grades II-IV for whom tumor status for IDH, 1p/19q codeletion, ATRX, and TERT had been determined. All cases were assigned to one of 5 groups following the WHO 2016 diagnostic criteria based on their morphologic features, and IDH and 1p/19q codeletion status. These groups are: (1) Oligodendroglioma, IDH-mutant and 1p/19q-codeleted; (2) Astrocytoma, IDH-mutant; (3) Glioblastoma, IDH-mutant; (4) Glioblastoma, IDH-wildtype; and (5) Astrocytoma, IDH-wildtype. Within each group, we used univariate and multivariate Cox proportional hazards models to assess associations of overall survival with patient age at diagnosis, grade, and ATRX alteration status and/or TERT promoter mutation status. Among Group 1 IDH-mutant 1p/19q-codeleted oligodendrogliomas, the TERT-WT group had significantly worse overall survival than the TERT-MUT group (HR: 2.72, 95% CI 1.05-7.04, p = 0.04). In both Group 2, IDH-mutant astrocytomas and Group 3, IDH-mutant glioblastomas, neither TERT mutations nor ATRX alterations were significantly associated with survival. Among Group 4, IDH-wildtype glioblastomas, ATRX alterations were associated with favorable outcomes (HR: 0.36, 95% CI 0.17-0.81, p = 0.01). Among Group 5, IDH-wildtype astrocytomas, the TERT-WT group had significantly better overall survival than the TERT-MUT group (HR: 0.48, 95% CI 0.27-0.87), p = 0.02). Thus, we present evidence that in certain WHO 2016 diagnostic groups, testing for TERT promoter mutations or ATRX alterations may provide additional useful prognostic information.

Immunomethylomic approach to explore the blood neutrophil lymphocyte ratio (NLR) in glioma survival

Background

Differentially methylated regions (DMRs) within DNA isolated from whole blood can be used to estimate the proportions of circulating leukocyte subtypes. We use the term “immunomethylomics” to describe the application of these immune lineage DMRs to studying leukocyte profiles. Here, we applied this approach to peripheral blood DNA from 72 glioma patients with molecularly defined brain tumors, representing common patient groups with defined characteristic survival times and risk factors. We first estimated the proportions of leukocyte subtypes in samples using deconvolution algorithms with reference DMR libraries from isolated leukocyte populations and Illumina 450K DNA methylation data. Then, we calculated the neutrophil to lymphocyte ratio (NLR) using methylation-derived cell composition estimates (mdNLR). The NLR is considered an indicator of immunosuppressive cells in cancer patients.

Results

Elevated mdNLR scores were observed in glioma patients compared to mdNLR values of published controls. Significantly decreased survival times were associated with mdNLR ≥ 4.0 in Cox proportional hazards models adjusted for age, gender, tumor grade, and molecular subtype (HR 2.02, 95% CI, 1.11–3.69). We also identified five myeloid-related CpGs that were highly correlated with the mdNLR (adjusted R² ≥ 0.80). Each of the five myeloid CpG loci was associated with survival when adjusted for the above covariates and offer a simplified approach for utilizing fresh or archived peripheral blood samples for interrogating a very small number of methylation markers to estimate myeloid immune influences in glioma survival.

Conclusions

The mdNLR (based on DNA methylation) is a novel candidate methylation biomarker that represents immunosuppressive myeloid cells within the blood of glioma patients with potential application in clinical trials and future epidemiologic studies of glioma risk and survival.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-017-0316-8) contains supplementary material, which is available to authorized users.

Approaching a Scientific Consensus on the Association between Allergies and Glioma Risk: A Report from the Glioma International Case-Control Study

BACKGROUND

Several previous studies have found inverse associations between glioma susceptibility and a history of allergies or other atopic conditions. Some evidence indicates that respiratory allergies are likely to be particularly relevant with regard to glioma risk. Using data from the Glioma International Case-Control Study (GICC), we examined the effects of respiratory allergies and other atopic conditions on glioma risk.

METHODS

The GICC contains detailed information on history of atopic conditions for 4,533 cases and 4,171 controls, recruited from 14 study sites across five countries. Using two-stage random-effects restricted maximum likelihood modeling to calculate meta-analysis ORs, we examined the associations between glioma and allergy status, respiratory allergy status, asthma, and eczema.

RESULTS

Having a history of respiratory allergies was associated with an approximately 30% lower glioma risk, compared with not having respiratory allergies (mOR, 0.72; 95% confidence interval, 0.58-0.90). This association was similar when restricting to high-grade glioma cases. Asthma and eczema were also significantly protective against glioma.

CONCLUSION

A substantial amount of data on the inverse association between atopic conditions and glioma has accumulated, and findings from the GICC study further strengthen the existing evidence that the relationship between atopy and glioma is unlikely to be coincidental.

IMPACT

As the literature approaches a consensus on the impact of allergies in glioma risk, future research can begin to shift focus to what the underlying biologic mechanism behind this association may be, which could, in turn, yield new opportunities for immunotherapy or cancer prevention.

Longer genotypically-estimated leukocyte telomere length is associated with increased adult glioma risk

Telomere maintenance has emerged as an important molecular feature with impacts on adult glioma susceptibility and prognosis. Whether longer or shorter leukocyte telomere length (LTL) is associated with glioma risk remains elusive and is often confounded by the effects of age and patient treatment. We sought to determine if genotypically-estimated LTL is associated with glioma risk and if inherited single nucleotide polymorphisms (SNPs) that are associated with LTL are glioma risk factors. Using a Mendelian randomization approach, we assessed differences in genotypically-estimated relative LTL in two independent glioma case-control datasets from the UCSF Adult Glioma Study (652 patients and 3735 controls) and The Cancer Genome Atlas (478 non-overlapping patients and 2559 controls). LTL estimates were based on a weighted linear combination of subject genotype at eight SNPs, previously associated with LTL in the ENGAGE Consortium Telomere Project. Mean estimated LTL was 31bp (5.7%) longer in glioma patients than controls in discovery analyses (P = 7.82x10-8) and 27bp (5.0%) longer in glioma patients than controls in replication analyses (1.48x10-3). Glioma risk increased monotonically with each increasing septile of LTL (O.R.=1.12; P = 3.83x10-12). Four LTL-associated SNPs were significantly associated with glioma risk in pooled analyses, including those in the telomerase component genes TERC (O.R.=1.14; 95% C.I.=1.03-1.28) and TERT (O.R.=1.39; 95% C.I.=1.27-1.52), and those in the CST complex genes OBFC1 (O.R.=1.18; 95% C.I.=1.05-1.33) and CTC1 (O.R.=1.14; 95% C.I.=1.02-1.28). Future work is needed to characterize the role of the CST complex in gliomagenesis and further elucidate the complex balance between ageing, telomere length, and molecular carcinogenesis.

A Heritable Missense Polymorphism in CDKN2A Confers Strong Risk of Childhood Acute Lymphoblastic Leukemia and Is Preferentially Selected during Clonal Evolution

Genome-wide association studies (GWAS) have identified SNPs in six genes that are associated with childhood acute lymphoblastic leukemia (ALL). A lead SNP was found to occur on chromosome 9p21.3, a region that is deleted in 30% of childhood ALLs, suggesting the presence of causal polymorphisms linked to ALL risk. We used SNP genotyping and imputation-based fine-mapping of a multiethnic ALL case-control population (Ncases = 1,464, Ncontrols = 3,279) to identify variants of large effect within 9p21.3. We identified a CDKN2A missense variant (rs3731249) with 2% allele frequency in controls that confers three-fold increased risk of ALL in children of European ancestry (OR, 2.99; P = 1.51 × 10(-9)) and Hispanic children (OR, 2.77; P = 3.78 × 10(-4)). Moreover, of 17 patients whose tumors displayed allelic imbalance at CDKN2A, 14 preferentially retained the risk allele and lost the protective allele (PBinomial = 0.006), suggesting that the risk allele provides a selective advantage during tumor growth. Notably, the CDKN2A variant was not significantly associated with melanoma, glioblastoma, or pancreatic cancer risk, implying that this polymorphism specifically confers ALL risk but not general cancer risk. Taken together, our findings demonstrate that coding polymorphisms of large effect can underlie GWAS "hits" and that inherited polymorphisms may undergo directional selection during clonal expansion of tumors.

Serum macrophage-derived chemokine/CCL22 levels are associated with glioma risk, CD4 T cell lymphopenia and survival time

Defects in antigen presenting cell function have been implicated in glioma immunosuppression. We measured peripheral CCL22, a dendritic cell/macrophage derived T cell trafficking chemokine, in sera from 1,208 glioma cases and 976 controls to assess whether it might provide a biomarker of glioma risk, survival and immune dysfunction. Cluster models were used to examine the relationship between CCL22 and glioma risk. Patient survival was assessed using Cox regression models. We also examined the relationship between CCL22 levels and CD4 cell counts, as well as allergy history and IgE levels. CCL22 levels were significantly lower among glioma cases compared with controls (Mean ± SEM: 1.23 ± 0.03 ng/mL in cases vs. 1.60 ± 0.03 ng/mL in controls, p < 0.0001) and this difference remained significant even after controlling for other covariates in the cluster models (highest quartile versus lowest Odds Ratio = 0.21, p < 0.0001). CD4 cell counts were positively correlated with CCL22 in glioma cases (Spearman r(2) = 0.51, p < 0.01) and were significantly lower in cases compared with controls. Higher CCL22 levels were associated with longer survival in all cases combined and in GBM cases (hazard ratio(allcases) = 0.81; 95% CI: 0.72-0.91, p = 0.0003). CCL22 levels were not associated with IgE level or self-reported allergies. Circulating CCL22 levels are related to both glioma risk and survival duration independent of age, histology, grade and IDH mutation status. CCL22 should be considered a marker of immune status with potential prognostic value.

Glioma Groups Based on 1p/19q, IDH, and TERT Promoter Mutations in Tumors

BACKGROUND

The prediction of clinical behavior, response to therapy, and outcome of infiltrative glioma is challenging. On the basis of previous studies of tumor biology, we defined five glioma molecular groups with the use of three alterations: mutations in the TERT promoter, mutations in IDH, and codeletion of chromosome arms 1p and 19q (1p/19q codeletion). We tested the hypothesis that within groups based on these features, tumors would have similar clinical variables, acquired somatic alterations, and germline variants.

METHODS

We scored tumors as negative or positive for each of these markers in 1087 gliomas and compared acquired alterations and patient characteristics among the five primary molecular groups. Using 11,590 controls, we assessed associations between these groups and known glioma germline variants.

RESULTS

Among 615 grade II or III gliomas, 29% had all three alterations (i.e., were triple-positive), 5% had TERT and IDH mutations, 45% had only IDH mutations, 7% were triple-negative, and 10% had only TERT mutations; 5% had other combinations. Among 472 grade IV gliomas, less than 1% were triple-positive, 2% had TERT and IDH mutations, 7% had only IDH mutations, 17% were triple-negative, and 74% had only TERT mutations. The mean age at diagnosis was lowest (37 years) among patients who had gliomas with only IDH mutations and was highest (59 years) among patients who had gliomas with only TERT mutations. The molecular groups were independently associated with overall survival among patients with grade II or III gliomas but not among patients with grade IV gliomas. The molecular groups were associated with specific germline variants.

CONCLUSIONS

Gliomas were classified into five principal groups on the basis of three tumor markers. The groups had different ages at onset, overall survival, and associations with germline variants, which implies that they are characterized by distinct mechanisms of pathogenesis. (Funded by the National Institutes of Health and others.).

Variants near TERT and TERC influencing telomere length are associated with high-grade glioma risk

Glioma, the most common central nervous system cancer in adults, has poor prognosis. Here we identify a new SNP associated with glioma risk, rs1920116 (near TERC), that reached genome-wide significance (Pcombined = 8.3 × 10(-9)) in a meta-analysis of genome-wide association studies (GWAS) of high-grade glioma and replication data (1,644 cases and 7,736 controls). This region has previously been associated with mean leukocyte telomere length (LTL). We therefore examined the relationship between LTL and both this new risk locus and other previously established risk loci for glioma using data from a recent GWAS of LTL (n = 37,684 individuals). Alleles associated with glioma risk near TERC and TERT were strongly associated with longer LTL (P = 5.5 × 10(-20) and 4.4 × 10(-19), respectively). In contrast, risk-associated alleles near RTEL1 were inconsistently associated with LTL, suggesting the presence of distinct causal alleles. No other risk loci for glioma were associated with LTL. The identification of risk alleles for glioma near TERC and TERT that also associate with telomere length implicates telomerase in gliomagenesis.

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

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.

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.

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.

Associations of high-grade glioma with glioma risk alleles and histories of allergy and smoking

Glioma risk has consistently been inversely associated with allergy history but not with smoking history despite putative biologic plausibility. Data from 855 high-grade glioma cases and 1,160 controls from 4 geographic regions of the United States during 1997-2008 were analyzed for interactions between allergy and smoking histories and inherited variants in 5 established glioma risk regions: 5p15.3 (TERT), 8q24.21 (CCDC26/MLZE), 9p21.3 (CDKN2B), 11q23.3 (PHLDB1/DDX6), and 20q13.3 (RTEL1). The inverse relation between allergy and glioma was stronger among those who did not (odds ratio(allergy-glioma) = 0.40, 95% confidence interval: 0.28, 0.58) versus those who did (odds ratio(allergy-glioma) = 0.76, 95% confidence interval: 0.59, 0.97; P(interaction) = 0.02) carry the 9p21.3 risk allele. However, the inverse association with allergy was stronger among those who carried (odds ratio(allergy-glioma) = 0.44, 95% confidence interval: 0.29, 0.68) versus those who did not carry (odds ratio(allergy-glioma) = 0.68, 95% confidence interval: 0.54, 0.86) the 20q13.3 glioma risk allele, but this interaction was not statistically significant (P = 0.14). No relation was observed between glioma risk and smoking (odds ratio = 0.92, 95% confidence interval: 0.77, 1.10; P = 0.37), and there were no interactions for glioma risk of smoking history with any of the risk alleles. The authors' observations are consistent with a recent report that the inherited glioma risk variants in chromosome regions 9p21.3 and 20q13.3 may modify the inverse association of allergy and glioma.

Distinct germ line polymorphisms underlie glioma morphologic heterogeneity

Two recent genome-wide association studies reported that single nucleotide polymorphisms (SNPs) in (or near) TERT (5p15), CCDC26 (8q24), CDKN2A/B (9p21), PHLDB1 (11q23), and RTEL1 (20q13) are associated with infiltrating glioma. From these reports, it was not clear whether the single nucleotide polymorphism associations predispose to glioma in general or whether they are specific to certain glioma grades or morphologic subtypes. To identify hypothesized associations between susceptibility loci and tumor subtype, we genotyped two case-control groups composed of the spectrum of infiltrating glioma subtypes and stratified the analyses by type. We report that specific germ line polymorphisms are associated with different glioma subtypes. CCDC26 (8q24) region polymorphisms are strongly associated with oligodendroglial tumor risk (rs4295627, odds ratio [OR] = 2.05, P = 8.3 × 10(-11)) but not glioblastoma risk. The opposite is true of RTEL (20q13) region polymorphisms, which are significantly associated with glioblastoma (rs2297440, OR = 0.56, P = 4.6 × 10(-10)) but not oligodendroglial tumor. The SNPs in or near CCDC26 (8q24) are associated with oligodendroglial tumors regardless of combined 1p and 19q deletion status; however, the association is greatest for those with combined deletion (rs4295627, OR = 2.77, P = 2.6 × 10(-9)). These observations generate hypotheses concerning the possible mechanisms by which specific SNPs (or alterations in linkage disequilibrium with such SNPs) are associated with glioma development.

Circulating levels of the innate and humoral immune regulators CD14 and CD23 are associated with adult glioma

Allergy history has been consistently inversely associated with glioma risk. Two serologic markers, soluble CD23 (sCD23) and soluble CD14 (sCD14), are part of the innate and adaptive humoral immune systems and modulate allergic responses in opposite directions, with sCD23 enhancing and sCD14 blunting inflammatory responses. We measured sCD23 and sCD14 in serum from blood that was drawn at a single time point from 1,079 glioma patients postdiagnosis and 736 healthy controls. Glioma was strongly associated with high sCD14 [highest versus lowest quartile odds ratio (OR), 3.94; 95% confidence interval (95% CI), 2.98-5.21] and low sCD23 (lowest versus highest quartile OR, 2.5; 95% CI, 1.89-3.23). Results were consistent across glioma histologic types and grades, but were strongest for glioblastoma. Whereas temozolomide treatment was not associated with either sCD14 or sCD23 levels among cases, those taking dexamethasone had somewhat lower sCD23 levels than those not taking dexamethasone. However, sCD23 was associated with case status regardless of dexamethasone treatment. These results augment the long-observed association between allergies and glioma and support a role for the innate and adaptive humoral functions of the immune system, in particular immunoregulatory proteins, in gliomagenesis.

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

Brain fibronectin expression in prenatally irradiated mice

Activation of gene transcription by radiation has been recently demonstrated in vitro. However, little is known on the specificity of these alterations on gene transcription. Prenatal irradiation is a known teratogen that affects the developing mammalian central nervous system (CNS). Altered neuronal migration has been suggested as a mechanism for abnormal development of prenatally irradiated brains. Fibronectin (FN), an extracellular glycoprotein, is essential for neural crest cell migration and neural cell growth. In addition, elevated levels of FN have been found in the extracellular matrix of irradiated lung. To test whether brain FN is affected by radiation, either FN level in insoluble matrix fraction or expression of FN mRNA was examined pre- and postnatally after irradiation. Mice (CD1), at 13 d of gestation (DG), served either as controls or were irradiated with gamma rays at 0.5 or 1 Gy. Control and irradiated animals were killed either at 13 DG, 14 DG, 17 DG, or 5, 6, or 14 d postnatal. Brain and liver were collected from offspring and analyzed for either total FN protein levels or relative mRNAs for FN and tubulin. Results of prenatal irradiation on reduction of postnatal brain weight relative to whole body weight and morphological reduction in cerebral cortex regions of postnatal brains are comparable to that reported by others. Insoluble matrix fraction (IMF) per gram of brain, liver, lung, and heart weight was not significantly different either between control and irradiated groups or between postnatal stages, suggesting that radiation did not affect the IMF. However, total amounts of FN in brain IMF at 17 DG were significantly different (p < .02) between normal (1.66 +/- 0.80 micrograms) and irradiated brains (0.58 +/- 0.22 microgram). FN mRNA was detectable at 13, 14, and 17 DG, but was not detectable at 6 and 14 d postnatal, indicating that FN mRNA is developmentally regulated. After 0.5 Gy of irradiation, expression of FN mRNA was reduced to 36% +/- 22% (1 h), 52% +/- 10% (1 d), and 76% +/- 10% (4 d) of the control level. After 1 Gy of irradiation, relative FN mRNA was 62% +/- 28% (1 h) and 75% +/- 3% (4 days) to the control level, respectively. This reduction was comparable to that reported by others for the cytoskeletal protein beta-actin. In contrast, mRNA for tubulin, another cytoskeletal protein, increased at 1 h after irradiation but then approached normal postnatally. The longer lasting alteration of FN may be more directly related to neural development, particularly if the reduction in FN is nonuniform.

The use of alkaline elution procedures to measure DNA damage in spermiogenic stages of mice exposed to methyl methanesulfonate

The ability of methyl methanesulfonate (MMS) to induce DNA breakage in spermiogenic stages of the mouse was studied using an alkaline elution technique. At daily intervals over a 3-week period following i.p. injection of 50 mg MMS/kg, mature spermatozoa were recovered from treated (3H-labeled) and control (14C-labeled) animals, lysed together on polycarbonate filters, and eluted with a high pH (12.2) buffer. Elution of germ-cell DNA from MMS-treated animals was found to increase in stages in which genetic damage from MMS is greatest. In general, the pattern of DNA elution from treated, spermiogenic stages paralleled the pattern of sensitivity to dominant lethals, specific-locus mutations and heritable translocations found by other investigators. It also paralleled the pattern of sperm-head methylation and protamine methylation measured in an earlier study (Sega and Owens, 1983). At 9 days post treatment (sperm sampled were in mid-to late-spermatid stages at the time of MMS exposure) the elution of sperm DNA did not change significantly over a pH range of 11.6-12.8, suggesting that, at the time of assay, DNA breaks were already present in the sperm. Because of the parallelism found between increased sperm DNA elution and increased genetic damage after mutagen treatment, alkaline elution may prove useful in monitoring potential genetic damage in human sperm.