|Home | About | Journals | Submit | Contact Us | Français|
Conceived and designed the experiments: JRE AM JJP. Performed the experiments: JRE AR IS AM JJP. Analyzed the data: JRE AR IS AM JJP. Contributed reagents/materials/analysis tools: JGH MSB NG CDJ. Wrote the paper: JRE IS CDJ AM JJP.
Glioblastoma (GBM) is a highly malignant brain tumor with a dismal prognosis. Gene expression profiling of GBM has revealed clinically relevant tumor subtypes, and this provides exciting opportunities to better understand disease pathogenesis. Results from an increasing number of studies demonstrate a role for the immune response in cancer progression, yet it is unclear how the immune response differs across tumor subtypes and how it affects outcome. Utilizing gene expression data from The Cancer Genome Atlas Project and the Gene Expression Omnibus database, we demonstrate an enrichment of immune response-related gene expression in the mesenchymal subtype of adult GBM (n=173) and pediatric high-grade gliomas (n=53). In an independent cohort of pediatric astrocytomas (n=24) from UCSF, we stratified tumors into subtypes and confirmed these findings. Using novel immune cell-specific gene signatures we demonstrate selective enrichment of microglia/macrophage-related genes in adult and pediatric GBM tumors of the mesenchymal subtype. Furthermore, immunostaining of adult GBM tumors showed significantly higher cell numbers of microglia/macrophages in mesenchymal versus non-mesenchymal tumors (p=0.04). Interestingly, adult GBM tumors with the shortest survival had significant enrichment of microglia/macrophage-related genes but this was not true for pediatric GBMs. Consistent with an association with poor outcome, immune response-related genes were highly represented in an adult poor prognosis gene signature, with the expression of genes related to macrophage recruitment and activation being most strongly associated with survival (p<0.05) using CoxBoost multivariate modeling. Using a microglia/macrophage high gene signature derived from quantification of tumor-infiltrating cells in adult GBM, we identified enrichment of genes characteristic of CD4 T cells, granulocytes, and microglia/macrophages (n=573). These studies support a role for the immune response, particularly the microglia/macrophage response, in the biology of an important subset of GBM. Identification of this subset may be important for future therapeutic stratification.
Glioblastoma (GBM), WHO grade IV astrocytoma, is a highly malignant disease with a poor prognosis despite aggressive treatment strategies . Results from molecular profiling studies of GBM suggest that stratification of these tumors into clinically relevant subtypes will lead to improved therapy outcomes as a result of individualizing treatment based on tumor molecular signatures. In adults, GBM is the most common primary malignant brain tumor. In children, high-grade astrocytomas, including GBM, are less common, comprising 8–12% of central nervous system (CNS) malignancies, but are equally heterogeneous and present with a similarly dismal outcome .
Gene expression profiling of adult GBM has greatly aided our understanding of underlying tumor pathogenesis –, and has promoted GBM stratification into tumor subtypes , . While many different subtype classifications have been performed, the mesenchymal subtype as defined by both Phillips et al.  and Verhaak et al.  shows the greatest concordance . The mesenchymal GBM subtype is associated with reduced survival , and it is characterized by upregulation of genes involved in tumor microenvironment interactions and processes , , , .
Although adult and pediatric high-grade astrocytomas are widely considered as having distinct clinical and molecular characteristics, a recent gene expression profiling study identified three major subtypes of pediatric high-grade glioma that recapitulated subtypes identified by Phillips et al. in adult high grade astrocytoma : HC1/proliferative, HC2/proneural and HC3/mesenchymal . Similar to adult tumors, genes upregulated in the HC3/mesenchymal subtype were associated with extracellular matrix-receptor interactions and cell adhesion . Given the subtype-specific increase in expression of genes involved in tumor microenvironment interactions, we reasoned that a subset of adult GBM and pediatric high-grade astrocytoma may have a pronounced and potentially important immune cell component.
Several lines of evidence suggest that the immune response is important to glioma biology. First, there is a significant inverse correlation between atopic disease and serologic markers of allergy, including IgE, and the sporadic incidence of glioma , . Secondly, single nucleotide polymorphisms of immune system genes are significantly associated with glioma risk . Lastly, manipulation of the immune response in experimental models for glioma alter disease outcome –. As a result of such observations, a number of novel immunomodulatory strategies are currently in clinical trials for treating GBM patients .
The predominant immune cell infiltrate in human GBM consists of CD45+CD11b+ microglia/macrophages ,  and in experimental models of glioma, tumor-associated microglia/macrophages directly influence tumor growth ,  and invasion . GBM also elicits an adaptive immune response and there is an increased infiltrate of CD4+FoxP3+ regulatory T cells in GBM relative to lower-grade gliomas and normal brain , , . These regulatory T cells are thought to contribute to the profound lack of activated effector CD8+ T cells in GBM  and their depletion in model systems of glioma have been shown to promote host survival .
Given the clinical and experimental evidence supporting a role for inflammation in GBM, we combined transcriptional expression profiling and immunoprofiling data on human high-grade astrocytoma to investigate the tumor-associated immune response across tumor subtypes and in relation to disease outcome in adult and pediatric astrocytoma.
Collection of samples used in this study was approved by the Human Research Protection Program Committee on Human Research of the University of California San Francisco. Each patient provided written consent for tissue collection, banking and research use by the UCSF Brain Tumor Research Center (BTRC).
Total RNA was extracted from 24 flash-frozen pediatric astrocytoma samples described previously , and included 20 high-grade astrocytomas (WHO grade III and IV astrocytoma) and 4 WHO grade II diffuse astrocytomas. Total cellular RNA was extracted using Qiagen’s RNA isolation kit according to the manufacturer’s protocol, and genomic DNA contamination was removed by an on-column DNase digestion step. Expression profiling was conducted by the UCSF Sandler Center Functional Genomics Core Facility using Agilent 4×44 K arrays (G42514F). Normalization of data was performed as follows: Raw Cy3 median signal intensities were recovered from Agilent Feature Extraction software output, converted to Log2 scale and quantile normalized between arrays. Duplicated probe signals were merged by taking their median values. The microarray data have been deposited in the GEO database (accession number GSE38330) and described in accordance with MIAME guidelines. Unsupervised hierarchical clustering was performed using Multiple Experiment Viewer Software (MeV; http://www.tm4.org/mev/). The gene signatures identified by Paugh et al.  for the subtype classification of pediatric high-grade astrocytomas were applied to the UCSF pediatric astrocytoma dataset. Spearman correlation was used as the distance metric and weighted average as the linkage method.
Gene expression data on 173 subtyped adult GBM patients  (Level 3, Agilent G4502A) was obtained from the TCGA Data Portal (http://cancergenome.nih.gov) accessed on April 15, 2010. Gene expression data from the entire TCGA cohort of adult GBM patients (n=506) as of June 16, 2011 was downloaded from the TCGA Data Portal (Level 3, Agilent G4502A). Raw gene expression data from the entire TCGA cohort of adult GBM patients (573 tumors and 10 normal brain controls as of December, 2011) was downloaded from the TCGA Data Portal (Level 1, Affymetrix HT Human Genome U133 Array). The data was pre-processed by RMA procedure  using custom CDF file, mapping good quality probes to unique NCBI Entrez genes [http://http://brainarray.mbni.med.umich.edu/Brainarray/Database/CustomCDF/CDF_download.asp, Version 14]. Raw gene expression data on 53 pediatric high-grade glioma patients (Affymetrix Human U133 Plus 2.0)  was accessed through the Gene Expression Omnibus (GEO) Website (http://www.ncbi.nlm.nih.gov/geo/, accession No. GSE19578) and was Log2 transformed and normalized using RMA.
Gene set enrichment analysis (GSEA)  was implemented using the Broad Institute GSEA v2.07 software (http://www.broadinstitute.org/gsea), the molecular signatures database (http://www.broadinstitute.org/gsea/msigdb), and the C5: GO gene sets database, comprised of 1454 gene sets named by GO terms and contains genes annotated by that term (www.geneontology.org). For all GSEA analyses 1000 phenotype permutations were performed and a Signal2Noise ranking metric was used to create the ranked list of genes. A false discovery rate (FDR) q-value of less than 0.25 (25%) was considered statistically significant. There was minimal overlap of genes between the subtype-specific signatures and the immune response-related gene sets: 28 of 300 genes in the adult mesenchymal signature and 29 of 236 genes in the pediatric HC3/mesenchymal signature. Gene signatures suggestive of specific immune cell subsets (Table S1) were also used to analyze the data. Cell-specific gene signatures of M1 and M2 polarized macrophages were compiled from analyzing multiple studies characterizing M1 and M2 macrophages ,  with the majority from a profiling study of in vitro polarized human monocytes by Martinez et al. . Hematopoietic cell lineage gene signatures (including monocytes, granulocytes, B cells, NK cells, naïve and activated T cells) were obtained from recent murine profiling , in which cells from peripheral blood, bone marrow or spleen were isolated and expression “fingerprints” specific to each cell type were determined. The glioma-infiltrating microglia/macrophage (GIM) gene signature was derived from differentially expressed genes (2-fold) in tumor-associated cells enriched for microglia/macrophages versus bulk tumor from a single GBM tumor .
CoxBoost was designed to develop proportional hazards models from microarray data and clinical covariates using a boosting approach . CoxBoost modeling to identify specific immune genes associated with survival was performed in R (http://www.rprojects.org). Microarray data used to develop the model included gene expression data of 58 genes derived from the GO terms Inflammatory Response (GO:0006954) and Response to Wounding (GO:0009611). Clinical variables included age at diagnosis, Karnofsky Performance Score (KPS) and gender. Overall survival was the outcome of interest and patients alive at the time of analysis were censored. The model was created using a training set of n=205 GBM patients from Carro et al. . Variables found to be associated with overall survival at p<0.05 were verified using a permutation test and the final model was subsequently assessed via the Brier score [defined as the squared difference between an event (i.e. death) occurrence and its predicted probability] in an independent cohort of n=301 GBM patients (TCGA Data Portal, 2011). A low Brier Score (i.e. close to 0) indicates increased predictability of the model .
Formalin-fixed paraffin embedded tissue microarray’s (TMAs) were generated from 34 adult GBM patients and 13 pediatric astrocytoma patients who had been transcriptionally profiled. Each tumor was represented at least twice on the array. Immunohistochemistry was performed according to standard methods and immunostaining for Iba1 (019-19741 at 12000, Wako Chemicals USA, Inc.), CA9 (NB100-417 at 11500, Novus USA) and CD34 (CBL496 at 1800, Chemicon) was performed on the Ventana Medical Systems Benchmark XT. Iba1 cell counts were performed using ImageJ Software, where a minimum of three fields at 400× magnification for each tumor were counted for Iba1 positive cells and then averaged. CA9 and CD34 staining was scored in a blinded fashion in a semi-quantitative manner as follows. CA9: 0, zero positivity; 1, <10% positivity; 2, 10–25% positivity; 3, >25% positive staining. CD34: 0, density of positive blood vessels less than normal brain; 1, density equal to normal brain; 2, density greater than normal brain. Mann-Whitney Rank Sum analysis, implemented in GraphPad Prism v5.0, was used to determine differences in immunostaining between patient groups. Correlations were performed using Spearman’s rho. A test-for-trend via a linear regression model was used to evaluate a trend in Iba1 staining between subtypes in the direction of mesenchymal, classical, neural, proneural. A p-value of <0.05 was considered statistically significant.
Iba1 scored tumors were divided into quartiles based on their microglia/macrophage cell counts and the highest quartile (high microglia/macrophages) was compared with the lowest quartile (low microglia/macrophages) for differential gene expression analysis. Significance Analysis of Microarrays [SAM ] was used to determine differentially expressed genes. The six genes identified as significantly differentially expressed were then applied to the whole TCGA GBM cohort (n=573), median centered, and hierarchical clustering was performed based on these genes. Clustering analysis was performed with Spearman correlation based distance and Ward linkage.
Adult GBM of the mesenchymal subtype have upregulation of genes involved in interactions with the tumor microenvironment , , . To address potential relationships with immune system response in this subtype we performed gene set enrichment analysis (GSEA) on gene expression data from adult GBM previously subtyped by Verhaak et al.  (n=173, TCGA Data Portal, 2010) (c5: GO gene sets). The top 50 enriched gene sets (NES: 1.74–1.89, FDR: 2.8–14%) in the mesenchymal subtype (n=56) versus the classical, neural and proneural subtypes (n=117) were grouped into broad GO classifications (Figure 1A and Table S2A). Strikingly, 48% of the enriched gene sets in the mesenchymal subtype were related to immune response processes, and included gene sets associated with development of the immune system, lymphocyte activation, response to infection or injury, and the adaptive immune response. Gene sets associated with signal transduction ranked second (8%) in the mesenchymal subtype and, of these, 3 of 4 were related to NF-κB signaling (Table S2A).
Recent gene expression profiling of pediatric high-grade gliomas (WHO grades III and IV/GBM, n=53)  identified a tumor subtype (defined as HC3/mesenchymal) with a gene expression signature similar to that of the adult mesenchymal subtype proposed by Phillips et al. . Similar to adult tumors, GSEA analysis of grade III and IV pediatric gliomas demonstrated enrichment of immune response-related gene sets in the HC3/mesenchymal subtype (Figure 1B and Table S2B). As the adult mesenchymal subtype proposed by Phillips et al.  and Verhaak et al.  are not identical, we used GSEA analysis to demonstrate enrichment of subtype-specific signature genes. Tumors defined as HC3/mesenchymal by Paugh et al.  were significantly enriched for the top 50 mesenchymal signature genes as defined by Verhaak et al.  (NES: 1.49, FDR: 5.4%, Figure S1A). Conversely, tumors defined as mesenchymal by Verhaak et al.  were significantly enriched for HC3/mesenchymal signature genes as defined by Paugh et al. (NES: 1.65, FDR: <0.001%, Figure S1B). Thereby, demonstrating similarities between the adult Verhaak mesenchymal subtype and the pediatric Paugh HC3/mesenchymal subtype.
In an independent cohort of pediatric astrocytomas from UCSF, we performed gene expression profiling (grades II, III and IV, n=24), applied 1035 signature genes proposed by Paugh et al.  to define tumor subtypes, and performed unsupervised hierarchical clustering. Three major clusters representative of the three previously defined subtypes were identified: HC1/proliferative, HC2/proneural, HC3/mesenchymal (Figure 2A–B). In the UCSF pediatric HC3/mesenchymal cohort, immune response GO gene sets were also significantly enriched (Table S3). Figure 2C, D, E demonstrates selective enrichment of the Immune Response gene set (GO:0006955) in the HC3/mesenchymal group (NES: 1.22, FDR: 23.7%), further supporting an association between the mesenchymal subtype and increased expression of immune response-related genes.
Having demonstrated an upregulation of immune response-related genes in the mesenchymal subtype we aimed to characterize the cell types potentially driving this pattern of gene expression. Immune cell-specific gene signatures (Table S1) were identified from recent profiling studies – and utilized for GSEA analysis. In adult GBM, six of 14 cell-specific gene signatures, were enriched in the mesenchymal subtype (NES: 1.53–1.64, FDR: 1.5–3.2%) (Figure 3A and Table S4A). The six enriched gene sets included gene signatures associated with microglia/macrophages, M1 and M2 macrophages and glioma infiltrating microglia/macrophages (GIM), along with the monocyte, granulocyte and naïve CD4 T cell gene sets. Similar to adults, the HC3/mesenchymal subtype of pediatric GBM (grade IV, n=38) demonstrated significant enrichment of the gene signatures associated with microglia/macrophages and monocytes (NES: 1.25–1.67, FDR: 1.6–24%) (Figure 3B and ). Unlike adults, the gene sets representing granulocytes and naïve CD4 T cells were not significantly enriched and two gene sets were negatively enriched; NK cells and nucleated erythrocytes (NES: −1.21 and −1.74, FDR: 22% and 0.1%, respectively). A similar analysis was performed with the independent UCSF cohort of pediatric astrocytomas. Due to the small number of GBMs (n=2) in the HC3/mesenchymal subtype (Figure S2), we included all astrocytoma grades (II, III and IV) in the analysis. Similar to the larger pediatric GBM data set, there was significant enrichment of the M2 macrophage gene signature (NES: 1.5, FDR: 18%) and significant negative enrichment of nucleated erythrocytes (NES: −1.49, FDR: 9.7%), in the HC3/mesenchymal subtype (Figure S3A and Table S4C).
Increased expression of microglia/macrophage-related genes could reflect increased numbers of immune cells or increased transcript expression levels per cell. To determine whether the mesenchymal subtype contained a greater number of microglia/macrophages, we generated tissue microarrays from adult subtyped GBM (n=24) and from UCSF pediatric subtyped astrocytomas (n=13), and analyzed microglia/macrophage number, as determined by Iba1 immunohistochemistry. While there was inter-tumoral heterogeneity, the mean microglia/macrophage number was significantly greater in the mesenchymal subtype of adult GBMs than in the non-mesenchymal subtype (Mann Whitney p=0.04, Figure 3C). Comparing across the four tumor subtypes there was a linear trend of decreasing numbers of microglia/macrophages in the order of mesenchymal, classical, neural and proneural (test-for-trend p=0.02, Figure S4A). Microglia/macrophage recruitment can be regulated via multiple mechanisms including levels of HIF1α . Interestingly, hypoxia, as determined by CA9, was significantly higher in tumors of the mesenchymal versus the non-mesenchymal subtype, (Mann-Whitney p=0.027, Figure 3D). While both immune cell recruitment and HIF1α can promote angiogenesis , we did not identify a significant difference in vascularity, as determined by CD34 immunostaining, between tumors of the mesenchymal and non-mesenchymal subtype (Mann-Whitney p=0.173, Figure 3E). In the UCSF pediatric cohort the tumors represented a diversity of grades (grades II, III and IV) and microglia/macrophage number did not differ between mesenchymal and non-mesenchymal tumors (Mann Whitney p=0.414, Figure S3B). Subtype specific comparisons are shown in Figure S4.
Previous studies have identified a gene signature associated with worse prognosis in adult GBM . To examine whether immune response-related genes are associated with survival, independent of tumor subtype, we performed Gene Ontology (GO) analysis of this signature gene set and 2 of the top 3 GO terms were immune response-related (“Inflammatory Response” and “Response to Wounding”, including a total of 58 independent genes, Table S5). To identify the immune response-related genes in this list of 58 that are most closely associated with survival in a robust and multivariate manner, we used CoxBoost modeling . The training set was comprised of gene expression data from 205 adult GBM previously used to define the worst prognosis signature  and the independent validation set was comprised of data from 301 independent GBM tumors (TCGA Data Portal, 2011). Clinical variables were included in the model, and as expected, age and KPS score were found to be significantly associated with survival (p<0.005, Table 1). In addition, the expression of 12 potential immune response genes were also significantly associated with overall survival: SPP1, PDPN, IL1RN, IL6, CCL8, MDK, IL1RAP, IRF7, MASP2, CASP3, TFRC and HRH4 (Table 1). Validation of the model in the independent data set resulted in a low Brier Score of 0.0793, demonstrating the strength of the model and the association between immune response gene expression and poor prognosis.
To assess if these survival-related immune response genes may be expressed by the microglia/macrophage population within the tumor mass we examined gene expression profiling data from tumor-associated cells enriched for microglia/macrophages compared with bulk tumor from a single GBM sample . As shown in Figure S5, eight of these immune response genes (SPP1, PDPN, IL1RN, IL6, CCL8, MDK, IRF7, and HRH4) were expressed between 1.25–1.59 fold higher in the microglia/macrophage enriched population compared to bulk tumor. This fold change was comparable to the microglia/macrophage specific markers CD45, CD11b and Iba1 (1.28–1.69 fold).
Rather than stratify tumors by subtype, adult GBM patients (n=506) were stratified based on survival (short survival, <6 months, n=98 and long survival, >24 months, n=78), and GSEA analysis demonstrated an enrichment of all macrophage, monocyte, and granulocyte gene signatures in short survival patients (NES: 1.39–1.53, FDR: 15–21%) (Figure 4A and Table S6A). In contrast to adult GBM, a similar stratification of pediatric GBM into short (<6 months, n=9) and long (>24 months, n=12) survival groups did not reveal a significant enrichment of immune cell-specific gene signatures in the short survival group (Figure 4B and Table S6B).
In the subset of adult GBM for which tissue was available for analysis, there was not a significant increase in mean number of microglia/macrophages in tumors from patients with shorter survival (p=0.252, Figure 4C). Interestingly, there was also no significant increase in hypoxia (p=0.209, Figure 4D) or vascularity (p=0.303, Figure 4E).
To identify factors associated with increased microglia/macrophage infiltrates, we stratified adult tumors (n=34) by number of tumor-associated microglia/macrophages (Iba1) and analyzed them by immunohistochemistry and differential gene expression analysis. Ignoring tumor subtype, there was no significant association between microglia/macrophage number and levels of hypoxia or vascularity in our cohort of adult GBM (p=0.674 and p=0.987 respectively, Figure S6). Comparing tumors with the highest and the lowest number of microglia/macrophages (top 25%, n=9 vs bottom 25%, n=9, Figure 5A) we identified 3 genes significantly upregulated and 3 genes significantly downregulated (Figure 5B and Table 2). Using these six genes as a surrogate marker for microglia/macrophage infiltrate, we identified four distinct tumor clusters based on unsupervised hierarchical clustering of gene expression data (n=573 tumors and n=10 normal samples, TCGA; Figure 5C). Eight of nine high microglia/macrophage tumors were in Cluster 3 and six of nine low microglia/macrophage tumors were in Cluster 2.
GSEA analysis of Cluster 3 (surrogate high microglia/macrophage, n=150) versus Cluster 2 (surrogate low microglia/macrophages, n=129) demonstrated significant enrichment of the M2 macrophage and monocyte gene signatures (NES: 1.47 and 1.62, FDR: 18% and 9.9%, respectively) (Figure 5D and Table S7), validating our immune cell type specific gene signatures. Furthermore, Cluster 3 had significant enrichment of the granulocyte and CD4 naïve T cell signatures (NES: 1.32 and 1.62, FDR: 24% and 19%, respectively), suggesting a potential association between high microglia/macrophage infiltrate and increased expression of characteristic granulocyte and CD4 T cell genes.
Protein analysis demonstrated tumors of the mesenchymal subtype exhibited high microglia/macrophage numbers in contrast to non-mesenchymal tumors that had prominent inter-tumoral heterogeneity of microglia/macrophage number (Figure 3). Using our high microglia/macrophage gene signature, tumors in Cluster 3 were enriched for tumors of the mesenchymal subtype, with 62% (34/55) of mesenchymal tumors falling into Cluster 3 (Figure 5E). Consistent with heterogeneity within tumor subtypes, Cluster 3 also contained non-mesenchymal tumors, including classical, neural and proneural tumors (Figure 5E). There was no significant difference in survival between patients in Cluster 3 versus Cluster 2 as a group or when analyzed within individual subtypes (data not shown).
In this study we analyzed primary human astrocytomas to investigate immune response indicators across tumor subtypes, and in relation to disease outcome. We demonstrate a striking enrichment of immune response-related genes and genes characteristic of myeloid cells including microglia, macrophages, and monocytes in the mesenchymal subtype of adult and pediatric GBM. Importantly, in adult GBMs, microglia/macrophage, monocyte and granulocyte cell signatures were significantly enriched in tumors from patients with the shortest survival. Interestingly, this was not the case for pediatric GBM: i.e., immune cell-specific signature genes were not enriched in children with the shortest survival. In adults, tumors having a gene expression signature associated with high microglia/macrophage numbers exhibited enhanced expression of characteristic myeloid, granulocyte and CD4 T cell genes and were enriched for tumors of the mesenchymal subtype. These data support marked inter-tumoral differences in the microglia/macrophage response.
Molecular stratification of adult GBM into subtypes has provided a framework for investigations into differential disease pathogenesis and tumor response to therapy –. However, to date, tumor profiling of GBM has concentrated on analysis of bulk tumor with little attention directed to the non-neoplastic elements, including immune cells, present in the tumor. In the current study we examined the pattern of gene expression signatures characteristic of specific immune cell populations. In adult GBM, immune response genes, including genes representative of microglia/macrophage populations, were enriched in tumors of the mesenchymal subtype. This pattern of gene expression was consistent with the increased numbers of microglia/macrophages identified in mesenchymal as compared to non-mesenchymal GBMs. Furthermore, we identified enrichment of microglia/macrophage-related signature genes in patients with the shortest survival. We suggest microglia/macrophage function may be particularly pronounced in a subset of GBM with aggressive clinical behavior. Further characterization of this subset may be particularly important for identifying patients for clinical trials involving immunomodulatory therapies. Indeed, a recent study of dendritic cell vaccination demonstrated improved responses in patients of the mesenchymal subtype compared to other subtypes . PLX3397, an inhibitor of macrophage colony stimulating factor 1, and other such therapies aimed at inhibiting the microglia/macrophage response, may also prove particularly effective in this molecularly defined subset of GBM.
Our data indicate an enrichment of microglia/macrophage-related signature genes in tumors from patients with short survival. As increased gene expression may reflect both differences in cell number and differences in microglia/macrophage phenotype or activation state, we examined microglia/macrophage number in available UCSF tumors. Increased absolute numbers of microglia/macrophages was not associated with short survival in this GBM cohort. This discrepancy with the gene expression data may be a result of intra-tumoral heterogeneity, alternatively it may suggest possible differences in gene expression levels per cell. To identify immune response-related genes most strongly associated with survival we performed CoxBoost modeling on immune response genes identified in a worse prognosis gene signature of adult GBM . In addition to established clinical variables, 12 potential immune response-related genes were found to be significantly associated with poor survival. These included genes known to be expressed by microglia (CCL8 , SPP1 , IRF7 , IL1RAP ), macrophages (IL1RN , SPP1 , PDPN , IL1RAP , MDK , TFRC , HRH4 ), and tumor-associated macrophages (IL6 , SPP1 ). Many of the corresponding gene products are also known to influence monocyte/microglia/macrophage recruitment and activation (IL6 , SPP1 , CCL8 , MDK ).
Some of the genes identified by CoxBoost modeling have been previously associated with GBM disease outcome. Notably, PDPN was included in a 9 gene prognostic signature defined by Colman et al. , the inflammatory cytokine IL6 has been linked to poor prognosis in GBM  and high serum levels of the secreted glycoprotein SPP1 (osteopontin) correlate with poor survival in GBM . Our studies, similar to all studies that analyze expression data from bulk tumor, are limited by the inability to definitively identify the cellular source of a particular transcript. However, analysis of gene expression data from purified human tumor-associated microglia/macrophages demonstrated increased expression of some of these genes (8/12) relative to bulk tumor. Thus total expression levels of many of these survival-related immune genes is likely derived in part from the microglia/macrophage infiltrate.
To elucidate other factors associated with a high microglia/macrophage infiltrate, adult tumors were stratified by microglia/macrophage number into high and low subsets and based on the differentially expressed genes, hierarchical clustering was performed on a large set of profiled GBMs. Tumors with a high microglia/macrophage gene expression signature demonstrated specific enrichment of macrophage, granulocyte and CD4 T cell gene sets. In mammary adenocarcinoma, a significant interplay between CD4+ T cells and macrophages in promoting invasion and metastasis has been demonstrated , and it is possible that such interactions may also occur in a subset of GBM. While tumors with a high microglia/macrophage gene expression signature were most commonly tumors of the mesenchymal subtype they also included other tumor subtypes suggesting potential biologic differences within the tumor subtypes.
The gene expression signature associated with high microglia/macrophage number was most strongly driven by two genes: KCNN2 and C14orf139. KCNN2 or SK2, a calcium activated potassium channel, is induced in activated microglia and helps to regulate the microglia respiratory burst function . KCNN2 expression has been previously reported in human glioma , and in melanoma KCNN2 is thought to help regulate hypoxia-induced cell proliferation . The function of C14orf139 is currently unknown.
While there exist some similarities between pediatric and adult high-grade astrocytomas, a substantial literature attests to anatomic, histologic, and molecular differences in adult versus pediatric disease , . Our data suggest the immune response may represent one of these differences. Recently, pediatric high-grade gliomas have been divided into three subtypes based on patterns of gene expression. While the pediatric HC3/mesenchymal subtype-specific signature genes are not identical to those used for adult tumors , , they share common elements with the adult mesenchymal subtype of Phillips et al. ,  and we demonstrate similarities with the mesenchymal subtype of Verhaak et al. . Using these signature genes  we were able to stratify an independent cohort of pediatric grade II, III and IV astrocytomas into the three defined subtypes. To our knowledge, this is the first study validating the stratification of pediatric astrocytomas into transcriptionally defined classes. These data emphasize the robustness of the gene signatures and subtypes defined by Paugh et al. .
Similar to adult tumors, the pediatric HC3/mesenchymal subtype was enriched in immune response-related genes in both the UCSF and Paugh et al.  pediatric cohorts. Furthermore, pediatric HC3/mesenchymal GBM were also enriched for gene signatures of microglia/macrophages and monocytes, emphasizing a potential subtype-specific role for these immune cells in adult and pediatric tumors. In contrast to adult GBM, the microglia/macrophage gene signatures were not significantly enriched in tumors from short survival patients among the pediatric GBM cohort. While these results must be interpreted with caution given the relatively small number of pediatric tumors analyzed, it emphasizes potentially important differences in tumor biology between adult and pediatric patients.
Taken together, our study demonstrates increased expression of immune response related genes, including microglia/macrophage signature genes, in a subset of adult and pediatric GBM. Understanding the factors that drive this differential immune response and its implications for therapeutic decision-making is critical. Future studies will be designed to elucidate these differences.
Comparison of gene expression across different mesenchymal tumor subtypes.
Distribution of tumor grade between tumor subtypes in two pediatric astrocytoma cohorts.
GSEA analysis of cell-specific gene signatures and number of microglia/macrophages in mesenchymal versus non-mesenchymal subtypes in the UCSF pediatric astrocytoma cohort.
Subtype specific comparisons of microglia/macrophage cell number, hypoxia and vascularity in adult and pediatric astrocytomas.
Expression of survival-associated immune response-related genes in glioma infiltrating microglia/macrophages relative to bulk tumor.
Correlation between microglia/macrophage cell number and hypoxia or vascularity in adult GBM.
Cell-specific gene signatures used for GSEA analysis.
Enrichment scores and statistics of the top 50 enriched GO gene sets from GSEA analysis of the mesenchymal versus non-mesenchymal adult and pediatric astrocytoma subtypes.
Enrichment scores and statistics of defense/immune response GO gene sets from GSEA analysis of the mesenchymal versus non-mesenchymal subtypes in the UCSF pediatric astrocytoma cohort.
Enrichment scores and statistics of cell-specific gene signatures from GSEA analysis of the mesenchymal versus non-mesenchymal subtypes in adult and pediatric astrocytomas.
Gene Ontology (Biological Process) analysis of the TCGA Worst Prognosis Signature.
Enrichment scores and statistics of cell-specific gene signatures from GSEA analysis of short versus long survival patients in adult and pediatric astrocytomas.
Enrichment scores and statistics of cell-specific gene signatures from GSEA analysis of Cluster 3 (surrogate high microglia/macrophage) versus Cluster 2 (surrogate low microglia/macrophage) in adult GBM.
We thank the BTRC Tissue Core and Cynthia Cowdrey for assistance with obtaining human tissue, preparing tissue microarrays, and immunohistochemistry.