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.

Natural killer (NK) cells are now recognized to exhibit characteristics akin to cells of the adaptive immune system. The generation of adaptive memory is linked to epigenetic reprogramming including alterations in DNA methylation. The study herein found reproducible genome wide DNA methylation changes associated with human NK cell activation. Activation led predominately to CpG hypomethylation (81% of significant loci). Bioinformatics analysis confirmed that non-coding and gene-associated differentially methylated sites (DMS) are enriched for immune related functions (i.e., immune cell activation). Known DNA methylation-regulated immune loci were also identified in activated NK cells (e.g., TNFA, LTA, IL13, CSF2). Twenty-one loci were designated high priority and further investigated as potential markers of NK activation. BHLHE40 was identified as a viable candidate for which a droplet digital PCR assay for demethylation was developed. The assay revealed high demethylation in activated NK cells and low demethylation in naïve NK, T- and B-cells. We conclude the NK cell methylome is plastic with potential for remodeling. The differentially methylated region signature of activated NKs revealed similarities with T cell activation, but also provided unique biomarker candidates of NK activation, which could be useful in epigenome-wide association studies to interrogate the role of NK subtypes in global methylation changes associated with exposures and/or disease states.

Quantitating the copy number of demethylated CpG promoter sites of the CD3Z gene can be used to estimate the numbers and proportions of T cells in human blood and tissue. Quantitative methylation specific PCR (qPCR) is useful for studying T cells but requires extensive calibration and is imprecise at low copy numbers. Here we compared the performance of a new digital PCR platform (droplet digital PCR or ddPCR) to qPCR using bisulfite converted DNA from 157 blood specimens obtained from ambulatory care controls and patients with primary glioma. We compared both ddPCR and qPCR with conventional flow cytometry (FACS) evaluation of CD3 positive T cells. Repeated measures on the same blood sample revealed ddPCR to be less variable than qPCR. Both qPCR and ddPCR correlated significantly with FACS evaluation of peripheral blood CD3 counts and CD3/total leukocyte values. However, statistical measures of agreement showed that linear concordance was stronger for ddPCR than for qPCR and the absolute values were closer to FACS for ddPCR. Both qPCR and ddPCR could distinguish clinically significant differences in T cell proportions and performed similarly to FACS. Given the higher precision, greater accuracy, and technical simplicity of ddPCR, this approach appears to be a superior DNA methylation based method than conventional qPCR for the assessment of T cells.

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.

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

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.

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.