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

Maternal immunoglobulin E and childhood leukemia.

Childhood leukemia, particularly acute lymphoblastic leukemia (ALL), has long been hypothesized to be affected by abnormal immune responses to microbial challenges stemming from a lack of immune modulation in early childhood. Studies of allergies suggest that a child's immune development may be modulated by maternal immune status. We conducted a study to explore the relationship between maternal immunoglobulin E (IgE) and childhood leukemia and to investigate whether maternal immune status can influence childhood leukemia risk. Serum total and specific IgE (respiratory and food) were measured in biological mothers of 352 children (193 healthy controls and 159 leukemia cases, including 139 ALL cases) ages <8 years who were enrolled in the Northern California Childhood Leukemia Study. Odds ratios associated with maternal IgE were calculated using unconditional logistic regression adjusted for child's age, sex, race/ethnicity, and annual household income. A positive association between childhood leukemia or ALL and elevated levels of maternal serum total IgE was observed, especially among Hispanics. In addition, a positive association was observed between childhood leukemia or ALL and maternal respiratory or food IgE status. These results suggest that maternal immune function may play a crucial role in the etiology of childhood leukemia, although additional studies need to be conducted to confirm the results of this study and provide a perspective on mechanisms.

Variants in the CDKN2B and RTEL1 regions are associated with high-grade glioma susceptibility.

The causes of glioblastoma and other gliomas remain obscure. To discover new candidate genes influencing glioma susceptibility, we conducted a principal component-adjusted genome-wide association study (GWAS) of 275,895 autosomal variants among 692 adult high-grade glioma cases (622 from the San Francisco Adult Glioma Study (AGS) and 70 from the Cancer Genome Atlas (TCGA)) and 3,992 controls (602 from AGS and 3,390 from Illumina iControlDB (iControls)). For replication, we analyzed the 13 SNPs with P < 10(-6) using independent data from 176 high-grade glioma cases and 174 controls from the Mayo Clinic. On 9p21, rs1412829 near CDKN2B had discovery P = 3.4 x 10(-8), replication P = 0.0038 and combined P = 1.85 x 10(-10). On 20q13.3, rs6010620 intronic to RTEL1 had discovery P = 1.5 x 10(-7), replication P = 0.00035 and combined P = 3.40 x 10(-9). For both SNPs, the direction of association was the same in discovery and replication phases.

Associations between polymorphisms in DNA repair genes and glioblastoma.

A pooled analysis was conducted to examine the association between select variants in DNA repair genes and glioblastoma multiforme, the most common and deadliest form of adult brain tumors. Genetic data for approximately 1,000 glioblastoma multiforme cases and 2,000 controls were combined from four centers in the United States that have conducted case-control studies on adult glioblastoma multiforme, including the National Cancer Institute, the National Institute for Occupational Safety and Health, the University of Texas M. D. Anderson Cancer Center, and the University of California at San Francisco. Twelve DNA repair single-nucleotide polymorphisms were selected for investigation in the pilot collaborative project. The C allele of the PARP1 rs1136410 variant was associated with a 20% reduction in risk for glioblastoma multiforme (odds ratio(CT or CC), 0.80; 95% confidence interval, 0.67-0.95). A 44% increase in risk for glioblastoma multiforme was found for individuals homozygous for the G allele of the PRKDC rs7003908 variant (odds ratio(GG), 1.44; 95% confidence interval, 1.13-1.84); there was a statistically significant trend (P = 0.009) with increasing number of G alleles. A significant, protective effect was found when three single-nucleotide polymorphisms (ERCC2 rs13181, ERCC1 rs3212986, and GLTSCR1 rs1035938) located near each other on chromosome 19 were modeled as a haplotype. The most common haplotype (AGC) was associated with a 23% reduction in risk (P = 0.03) compared with all other haplotypes combined. Few studies have reported on the associations between variants in DNA repair genes and brain tumors, and few specifically have examined their impact on glioblastoma multiforme. Our results suggest that common variation in DNA repair genes may be associated with risk for glioblastoma multiforme.

Pathway analysis of single-nucleotide polymorphisms potentially associated with glioblastoma multiforme susceptibility using random forests.

Glioma is a complex disease that is unlikely to result from the effect of a single gene. Genetic analysis at the pathway level involving multiple genes may be more likely to capture gene-disease associations than analyzing genes one at a time. The current pilot study included 112 Caucasians with glioblastoma multiforme and 112 Caucasian healthy controls frequency matched to cases by age and gender. Subjects were genotyped using a commercially available (ParAllele/Affymetrix) assay panel of 10,177 nonsynonymous coding single-nucleotide polymorphisms (SNP) spanning the genome known at the time the panel was constructed. For this analysis, we selected 10 pathways potentially involved in gliomagenesis that had SNPs represented on the panel. We performed random forests (RF) analyses of SNPs within each pathway group and logistic regression to assess interaction among genes in the one pathway for which the RF prediction error was better than chance and the permutation P < 0.10. Only the DNA repair pathway had a better than chance classification of case-control status with a prediction error of 45.5% and P = 0.09. Three SNPs (rs1047840 of EXO1, rs12450550 of EME1, and rs799917 of BRCA1) of the DNA repair pathway were identified as promising candidates for further replication. In addition, statistically significant interactions (P < 0.05) between rs1047840 of EXO1 and rs799917 or rs1799966 of BRCA1 were observed. Despite less than complete inclusion of genes and SNPs relevant to glioma and a small sample size, RF analysis identified one important biological pathway and several SNPs potentially associated with the development of glioblastoma.

Nonsynonymous coding single-nucleotide polymorphisms spanning the genome in relation to glioblastoma survival and age at diagnosis.

PURPOSE: Our aim was to discover possible inherited factors associated with glioblastoma age at diagnosis and survival. Although new genotyping technologies allow greatly expanded exploration of such factors, they pose many challenges. EXPERIMENTAL DESIGN: In this pilot study, we (a) genotyped 112 newly diagnosed glioblastoma patients ascertained through a population-based study (group 1) with the ParAllele assay panel of approximately 10,000 nonsynonymous coding single-nucleotide polymorphisms (SNP), (b) used several statistical and bioinformatic techniques to identify 17 SNPs potentially related to either glioblastoma age at diagnosis or survival, and (c) genotyped 16 of these SNPs using conventional PCR methods in an independent group of 195 glioblastoma patients (group 2). RESULTS: In group 2, only one of the 16 SNPs, rs8057643 (located on 16p13.2), was significantly associated with glioblastoma age at diagnosis (nominal P = 0.0017; Bonferroni corrected P = 0.054). Median ages at diagnosis for those with 0, 1, or 2 T alleles were 66, 57, and 59 years in group 1 and 64, 57, and 55 years in group 2 (combined P = 0.001). Furthermore, Cox regression analyses of time to death with number of T alleles adjusted for gender and patient group yielded a hazard ratio of 0.82 (95% confidence interval, 0.68-0.98; P = 0.03). CONCLUSIONS: Although limited by a relatively small sample size, this pilot study, using well-characterized, unambiguous disease characteristics, illustrates the necessity of independent replication owing to the likelihood of false positives. Several other challenges are discussed, including attempts to incorporate information on the potential functional importance of SNPs in genome-disease association studies.

Serum IgE, tumor epidermal growth factor receptor expression, and inherited polymorphisms associated with glioma survival.

In population-based glioma patients, we examined survival in relation to potentially pertinent constitutive polymorphisms, serologic factors, and tumor genetic and protein alterations in epidermal growth factor receptor (EGFR), MDM2, and TP53. Subjects were newly diagnosed adults residing in the San Francisco Bay Surveillance Epidemiology and End Results Area during 1991 to 1994 and 1997 to 1999 with central neuropathology review (n = 873). Subjects provided blood for serologic studies of IgE and IgG to four herpes viruses and constitutive specimens for genotyping 22 polymorphisms in 13 genes (n = 471). We obtained 595 of 697 astrocytic tumors for marker studies. We determined treatments, vital status, and other factors using registry, interview, medical record, and active follow-up data. Cox regressions for survival were adjusted for age, gender, ethnicity, study series, resection versus biopsy only, radiation, and chemotherapy. Using a stringent P < 0.001, glioma survival was associated with ERCC1 C8092A [hazard ratio (HR), 0.72; 95% confidence limits (95% CL), 0.60-0.86; P = 0.0004] and GSTT1 deletion (HR, 1.64; 95% CL, 1.25-2.16; P = 0.0004); glioblastoma patients with elevated IgE had 9 months longer survival than those with normal or borderline IgE levels (HR, 0.62; 95% CL, 0.47-0.82; P = 0.0007), and EGFR expression in anaplastic astrocytoma was associated with nearly 3-fold poorer survival (HR, 2.97; 95% CL, 1.70-5.19; P = 0.0001). Based on our and others' findings, we recommend further studies to (a) understand relationships of elevated IgE levels and other immunologic factors with improved glioblastoma survival potentially relevant to immunologic therapies and (b) determine which inherited ERCC1 variants or other variants in the 19q13.3 region influence survival. We also suggest that tumor EGFR expression be incorporated into clinical evaluation of anaplastic astrocytoma patients.

Molecular features of adult glioma associated with patient race/ethnicity, age, and a polymorphism in O6-methylguanine-DNA-methyltransferase.

BACKGROUND: Risk factors for adult glioma in the San Francisco Bay Area include well-known demographic features such as age and race/ethnicity, and our previous studies indicated that these characteristics are associated with the TP53 mutation status of patients' tumors. We enlarged our study to assess the relationships of risk factors with TP53 as well as epidermal growth factor receptor (EGFR) and murine double minute-2 (MDM2) gene amplification and expression and the germ line Leu84Phe polymorphism in the DNA repair protein O6-methylguanine-DNA-methyltransferase (MGMT). MGMT expression may depend on the TP53 status of cells. METHODS: Molecular analyses were carried out on 556 incident astrocytic tumors. MGMT genotype data were collected on germ line DNA from 260 of these cases. RESULTS: The tumor data confirm the inverse relationships between TP53 mutation and MDM2 (P = 0.04) or EGFR (P = 0.004) amplification and that patients whose tumors contain TP53 mutations are younger than those without (P < 0.001). Although there was little difference in age of patient by EGFR amplification or expression among glioblastoma multiforme cases, EGFR gene amplification was associated with much older age of onset of anaplastic astrocytoma; for example, EGFR-amplified anaplastic astrocytoma cases were on average 63 years old compared with 48 years for nonamplified cases (P = 0.005). An increased prevalence of TP53 mutation positive glioblastoma multiforme was noted among nonwhites (African American and Asian) compared with whites (Latino and non-Latino; P = 0.004). Carriers of the MGMT variant 84Phe allele were significantly less likely to have tumors with TP53 overexpression (odds ratio, 0.30; 95% confidence interval, 0.13-0.71) and somewhat less likely to have tumors with any TP53 mutation (odds ratio, 0.47; 95% confidence interval, 0.13-1.69) after adjusting for age, gender, and ethnicity. Interestingly, EGFR gene amplification and EGFR protein overexpression were also inversely associated with the MGMT 84Phe allele. CONCLUSIONS: Our results are consistent with ethnic variation in glioma pathogenesis. The data on MGMT show that an inherited factor involving the repair of methylation and other alkylation damage, specifically to the O6 position of guanine, may be associated with the development of tumors that proceed in their development without TP53 mutations or accumulation of TP53 protein and possibly also those that do not involve amplification of the EGFR locus.