Abstract
Purpose: The molecular heterogeneity of glioblastoma has been well recognized and has resulted in the generation of molecularly defined subtypes. These subtypes (classical, neural, mesenchymal, and proneural) are associated with particular signaling pathways and differential patient survival. Less understood is the correlation between these glioblastoma subtypes with immune system effector responses, immunosuppression, and tumor-associated and tumor-specific antigens. The role of the immune system is becoming increasingly relevant to treatment as new agents are being developed to target mediators of tumor-induced immunosuppression, which is well documented in glioblastoma.
Experimental Design: To ascertain the association of antigen expression, immunosuppression, and effector response genes within glioblastoma subtypes, we analyzed the Cancer Genome Atlas (TCGA) glioblastoma database.
Results: We found an enrichment of genes within the mesenchymal subtype that are reflective of antitumor proinflammatory responses, including both adaptive and innate immunity and immunosuppression.
Conclusions: These results indicate that distinct glioma antigens and immune genes show differential expression between glioblastoma subtypes and this may influence responses to immunotherapeutic strategies in patients depending on the subtype of glioblastoma they harbor. Cancer Immunol Res; 1(2); 112–22. ©2013 AACR.
Introduction
Glioblastoma remains the most common malignant primary brain tumor in adults. Despite aggressive treatment with surgical resection, radio-, and chemotherapy, the tumor ultimately recurs. A major obstacle to treatment is the subversion of the immune system by the tumor to facilitate proliferation and malignant degeneration of tumor cells. Immunosuppression is thought to play a major role in the aggressive nature of gliomas and their resistance to current therapies. The presence of immunosuppressive infiltrates such as FoxP3+ regulatory T cells (Treg) and M2 macrophages has been documented in gliomas and in some cases, correlates with prognosis (1–3). In addition, immune responses to glioma can be impaired by the tumor itself through expression of immunosuppressive cytokines such as TGF-β and the upregulation of the STAT3 (4). Genomic profiling of glioblastomas has shown that up to 4 genotypic subtypes exist, with 2 showing marked differences in gene expression and patient survival (proneural and mesenchymal; ref. 5). The STAT3 pathway has been shown to be the key molecular driver of the mesenchymal transformation within glioblastoma (6, 7). In addition, STAT3 has been implicated in many mechanisms of tumor-mediated immunosuppression (8–13) and is a negative prognosticator for survival in mice (14) and in human patients with anaplastic astrocytoma (15). These data would suggest that the mesenchymal glioblastoma subtype may be more immunosuppressive compared with other glioblastoma subtypes and possibly more refractory to immunotherapy.
Immunotherapy is an appealing treatment of gliomas because it allows for tumor specificity while minimizing collateral damage to normal brain tissue. Clinical trials using dendritic cell or peptide vaccines to target glioma cells have shown promising results (16–19). However, as with other treatments, only subsets of patients respond. This may be due to molecular and genomic factors that affect interactions between the immune system and the tumor. A recent study using glioma lysate-pulsed dendritic cell vaccination showed that glioblastomas of the mesenchymal phenotype had higher levels of CD3+ and CD8+ tumor-infiltrating lymphocytes than glioblastomas of other subtypes (20). Furthermore, patients whose tumors had the mesenchymal gene signature, survived longer after dendritic cell vaccination than controls of the same genetic subtype (20). Although this finding tends to contradict what is thought about the immunosuppressive nature of mesenchymal gliomas with respect to aggressiveness and poor patient survival, it suggests that mesenchymal gliomas may be more immunogenic and more responsive to immunotherapy. Interestingly, robust immunosuppressive Treg infiltration into the glioblastoma microenvironment is almost always observed with a corresponding influx of T effector cells (2). Thus, one may speculate that the tumor becomes selectively more immunosuppressive as a reaction to antitumor effector responses. To investigate this paradox, we analyzed mRNA expression levels of immune system genes among the various glioblastoma subtypes using The Cancer Genome Atlas (TCGA) database and found not only a preferential enrichment of immunosuppressive genes but also an enrichment of immune effector genes within the mesenchymal subset.
Materials and Methods
TCGA data acquisition
It should be noted that the composition of the analyzed glioblastoma tissue from the TCGA database contains genetic information from both glioma cells and tumor-supportive stroma cells including infiltrating immune cells because these components were not selectively eliminated from the specimen, which was used in the construction of the database. The glioblastoma cancer study set from the TCGA database consists of defined glioblastoma subtypes (21). The TCGA database was analyzed in 2 distinct ways: (i) using online knowledge bases [Ingenuity Pathway Analysis (IPA) and Uniprot Protein Knowledgebase (www.uniprot.org)] to define genes associated with immune responses; and (ii) using a collated list of immune response genes from the literature with an emphasis on those previously documented to have a role in glioblastoma. For the first analysis, mRNA expression levels (available from the TCGA Research Network website: cancergenome.nih.gov), of the proneural (n = 107), mesenchymal (n = 119), classical (n = 115), and neural (n = 58) glioblastoma subsets were used as the source data available as of February 2013. This first analysis compares mRNA expression levels using the mesenchymal subset as the reference relative to the other subsets. The second analysis was through the open access cBio Cancer Genomics Portal at www.cbioportal.org (22), in which proneural (n = 141), mesenchymal (n = 160), classical (n = 147), and neural (n = 96) were used as the source data. This second analysis compares specific mRNA expression levels between glioblastoma subsets. To analyze mRNA expression based on Agilent microarray of selected genes in the 4 subtypes of glioma, the z score threshold for all genes (the number of SDs above the mean expression level of the selected gene) was set to more than 1 for each of the subtypes [an example of the syntax used for the search is: interleukin (IL)-10: EXP>1].
Compilation of immune genes and tumor antigens
The immune effector gene sets were defined by:
The IPA tool to identify immune activators with the following characteristics: cytokine, enzyme, G-protein–coupled receptor, growth factor, kinase, peptidase, phosphatase, transcription factor, transcription regulator, transmembrane receptor, transporter, and unknown and secondarily the Uniprot knowledge base (www.uniprot.org) for the terms “activation of immune response” and “immune effector process.” Of 814 genes identified by IPA and Uniprot, 734 could be found in the TCGA database and were used for the analysis.
By collating a list from the literature of documented proinflammatory effector cytokines (IFN-γ, IL-1, IL-2, IL-4, IL-7, IL-12, IL-15, and TNF-α); surface markers reflective of the presence of immune effector cell responses (CD3 and CD8) and their associated immune-activating markers (CD80, CD86, CD40, and HLA); signaling pathways reflective of adaptive immune activation (NF-κB, STAT1, IRAK, STAT4, T-bet, DNAM1, and IRF7); and innate immunity and markers reflective of activation (NKp46, NKG2D, NKp30, NKp44, KLRD1, TREM1, TREM2, MIF, TLR2, TLR3, TLR4, and TLR9) that participate in antitumor effector responses.
The immunosuppressive genes were defined by:
The IPA tool to identify mediators of immunosuppression with the following characteristics: cytokine, enzyme, G-protein–coupled receptor, growth factor, kinase, peptidase, phosphatase, transcription factor, transcription regulator, transmembrane receptor, transporter, and unknown and secondarily the Uniprot knowledge base for the search terms related to tumor-mediated immunosuppression including: “inhibition of immunity,” “inhibition of immune system process,” “inhibition of effector immunity,” and “inhibition of immune response to tumor cell.” Of 235 genes identified by IPA and Uniprot, 218 could be found in the TCGA database and were used for the analysis.
By collating a list from the literature of documented immunosuppressive cytokines and mechanisms (galectin-3, VEGF, IL-10, IL-23, TGF-β, PD-1, PD-L1, and CTLA-4); chemokines (CSF-1, CCL2, and CCL22); tumor-supportive and immunosuppressive myeloid- and monocyte-related genes (CD163, CD204, MIC-1, arginase, and CD47); immunosuppressive signaling pathways (IL-6, gp130, Jak2, STAT3, Pim-1, SOCS3, STAT5A, and STAT5B); and markers related to Tregs (CD4, ICOS, IDO, and FoxP3) that participate in tumor-mediated immunosuppression.
Statistical analysis
A χ2 test was used to compare the number of patients between glioblastoma subtypes in which mRNA expression levels of a selected gene were at least one SD above mean expression in a particular subtype. To identify differentially expressed genes between the subtypes, a modified two-sample t test using the limma package was applied. Genes that were differentially expressed (DE) on the basis of the false discovery rate (FDR) of at least <0.01 in the mesenchymal subset relative to the proneural (∼7,000), neural (∼9,000), and classical (∼6,000) subsets were then identified. The β-uniform mixture (BUM) model, described by Pounds and Morris (26), was used to control FDR. This list of differentially expressed genes was then dichotomized into 2 subgroups: genes (mRNAs) overexpressed in the mesenchymal subtype [designated DE(m>x); where m = mesenchymal and x = subtype] and vice versa [DE(x>m)]. Pearson χ2 with Yates' continuity correction and hypergeometric tests were used to determine the enrichment of the immune gene sets within these 2 groups.
Results
Immune genes have differential expression in glioblastoma subtypes
Analysis of the IPA-selected immune activators and suppressors revealed that the greatest immunologic diversity exists between the proneural and mesenchymal subsets (Fig. 1). Only 17% (n = 123) of the immune activators and 19% (n = 41) of the immunosuppressive genes shared differential expression among all glioblastoma subtypes (proneural, classical, and neuronal) relative to the mesenchymal subset. This diversity of the immune activating and immunosuppressive genes can distinguish the various glioblastoma subtypes on heatmaps (Fig. 2). To ascertain which immune genes were differentially enriched in the mesenchymal subset in comparison with the proneural subset, the immune genes were ranked and then compared with the overall 17,000 genes from the TCGA dataset (Table 1). Using both χ2 and hypergeometric tests to ascertain if the immune genes were selected by chance, a preferential enrichment of both immune activators (Pearson χ2 test with Yates' continuity correction, P < 2.2 × 10−16; hypergeometric test, P < 2.2 × 10−16) and immunosuppressors (Pearson χ2 test with Yates' continuity correction, P = 3.5 × 10−5; hypergeometric test, P = 2.04 × 10−5) was found in the mesenchymal subset relative to the proneural subset. Similarly, in the comparison between the classical and mesenchymal subset, immune activators (Pearson χ2 test with Yates' continuity correction, P < 2.2 × 10−16; hypergeometric test, P <2.2 × 10−16) and immunosuppressors (Person χ2 test with Yates' continuity correction, P = 2.2 × 10−08; hypergeometric test, P = 4.88 × 10−8) were found preferentially enriched in the mesenchymal subset. However, this preferential enrichment was not as evident in the comparison between the mesenchymal and neural subsets [immune activators (Pearson χ2 test with Yates' continuity correction, P = 0.2501; hypergeometric test, P = 0.8915) and immunosuppressors (Pearson χ2 test with Yates' continuity correction, P = 0.5736; hypergeometric test, P = 0.2847)].
Venn diagrams showing immune activators (left) and immunosuppressors (right) that are differentially expressed at a FDR of 0.001.
Venn diagrams showing immune activators (left) and immunosuppressors (right) that are differentially expressed at a FDR of 0.001.
Heatmaps of the 734 immune activator genes (left) and 218 immunosuppressor genes (right) compared between glioblastoma subtypes. The expression values shown on the heatmap have been standardized and at ±2 SDs for display purposes. The scale of the values is indicated in the color key.
Heatmaps of the 734 immune activator genes (left) and 218 immunosuppressor genes (right) compared between glioblastoma subtypes. The expression values shown on the heatmap have been standardized and at ±2 SDs for display purposes. The scale of the values is indicated in the color key.
Ranking of differentially expressed immune genes between the mesenchymal (M) and proneural (P) subsets
Gene . | Role . | Fold change M relative to P . | P value . | Rank among all genes . |
---|---|---|---|---|
Immune activators | ||||
TIMP1 | Regulates resistance to infection | 4.3 | 4.2 × 10−45 | 12 |
SERPING1 | Complement activation | 2.3 | 1.1 × 10−34 | 71 |
TNFRSF14 | TNF family member | 2.3 | 3.1 × 10−34 | 81 |
LGALS3 | Galectin-3 that induces immunosuppression | 4.8 | 3.0 × 10−33 | 93 |
TNFAIP3 | Terminates TNF-induced NF-κB responses | 2.0 | 2.6 × 10−30 | 139 |
CCR2 | Monocyte chemoattractant protein | 2.5 | 7.0 × 10−28 | 213 |
IL-15 | Enhances immunity of CD8+ T cells | 2.3 | 2.2 × 10−27 | 232 |
LPXN | Substrate for tyrosine kinase | 1.7 | 4.0 × 10−27 | 238 |
FAS | T-cell apoptosis | 2.7 | 8.3 × 10−26 | 294 |
TNFSF4 | Encodes OX40 ligand | 2.3 | 1.8 × 1 0−25 | 305 |
PTGER4 | Prostaglandin E receptor 4 | 2.2 | 2.2 × 10−24 | 346 |
CCL2 | Monocyte chemotactic protein | 3.7 | 9.9 × 10−24 | 386 |
SOX11 | Signal transducer molecule | −2.5 | 1.3 × 10−23 | 401 |
LYN | Inhibitory role in myeloid proliferation | 1.5 | 3.2 × 10−22 | 492 |
THYBS1 | Thrombospondin 1; assists in tumor death | 1.9 | 8.1 × 10−22 | 532 |
IRAK3 | Negative regulator of TLR signaling | 2.0 | 2.3 × 10−21 | 568 |
GDF15 | TGF-β superfamily member | 2.8 | 7.1 × 10−21 | 611 |
PRDM1 | Represses IFN-β gene expression | 1.6 | 7.8 × 10−21 | 617 |
MICB | Activates cytolytic response of NK and CD8 T cells | 1.5 | 5.6 × 10−20 | 696 |
IL-6R | Key immunosuppressive pathway | 1.7 | 6.4 × 10−20 | 707 |
Immunosuppressors | ||||
TIMP1 | Regulates resistance to infection | 4.3 | 4.2 × 10−45 | 12 |
ANXA1 | Enhances/inhibits adaptive immunity | 3.3 | 3.4 × 10−41 | 23 |
TNFRSF1A | Tumor necrosis superfamily protein | 2.2 | 5.4 × 10−40 | 31 |
MR1 | Antigen presentation function | 2.1 | 2.2 × 10−37 | 45 |
KLRC2 | Encodes NKG2-C (expressed on NK cells) | −10.9 | 1.3 × 10−36 | 50 |
SERPING1 | Complement activation | 2.3 | 1.1 × 10−34 | 71 |
LTBP1 | Activates TGF-β | 3.2 | 1.0 × 10−32 | 100 |
SWAP70 | Mediates IgE responses/antibody switch | 2.1 | 2.1 × 10−32 | 109 |
SATB1 | Silencing mimics the effects of IFN-γ treatment | −2.4 | 4.5 × 10−32 | 114 |
TGFBR2 | TGF-β2 receptor | 1.7 | 6.0 × 10−32 | 120 |
SHC1 | Inhibits adaptive immune responses | 1.9 | 1.0 × 10−31 | 125 |
LGALS1 | Galectin 1—deactivates M1 macrophages | 2.3 | 1.1 × 10−31 | 127 |
RELB | Interacts with NF-κβ2 | 2.0 | 4.1 × 10−31 | 135 |
DLL1 | T-cell activation | −2.0 | 4.4 × 10−30 | 142 |
BCL3 | Regulates NF-κβ | 2.7 | 6.3 × 10−30 | 146 |
RAB27A | Granulocyte exocytosis | 2.2 | 1.3 × 10−29 | 157 |
NOD1 | Innate immunity | 1.8 | 1.4 × 10−29 | 160 |
KLRC1 | Recognition of MHC class I molecules by NK cells | −3.0 | 1.5 × 10−28 | 161 |
CHI3L1 | YKL-40 immunosuppressor in gliomas | 3.0 | 2.2 × 10−29 | 170 |
ICAM1 | Transmigration receptor for immune cells | 2.7 | 4.9 × 10−29 | 182 |
Gene . | Role . | Fold change M relative to P . | P value . | Rank among all genes . |
---|---|---|---|---|
Immune activators | ||||
TIMP1 | Regulates resistance to infection | 4.3 | 4.2 × 10−45 | 12 |
SERPING1 | Complement activation | 2.3 | 1.1 × 10−34 | 71 |
TNFRSF14 | TNF family member | 2.3 | 3.1 × 10−34 | 81 |
LGALS3 | Galectin-3 that induces immunosuppression | 4.8 | 3.0 × 10−33 | 93 |
TNFAIP3 | Terminates TNF-induced NF-κB responses | 2.0 | 2.6 × 10−30 | 139 |
CCR2 | Monocyte chemoattractant protein | 2.5 | 7.0 × 10−28 | 213 |
IL-15 | Enhances immunity of CD8+ T cells | 2.3 | 2.2 × 10−27 | 232 |
LPXN | Substrate for tyrosine kinase | 1.7 | 4.0 × 10−27 | 238 |
FAS | T-cell apoptosis | 2.7 | 8.3 × 10−26 | 294 |
TNFSF4 | Encodes OX40 ligand | 2.3 | 1.8 × 1 0−25 | 305 |
PTGER4 | Prostaglandin E receptor 4 | 2.2 | 2.2 × 10−24 | 346 |
CCL2 | Monocyte chemotactic protein | 3.7 | 9.9 × 10−24 | 386 |
SOX11 | Signal transducer molecule | −2.5 | 1.3 × 10−23 | 401 |
LYN | Inhibitory role in myeloid proliferation | 1.5 | 3.2 × 10−22 | 492 |
THYBS1 | Thrombospondin 1; assists in tumor death | 1.9 | 8.1 × 10−22 | 532 |
IRAK3 | Negative regulator of TLR signaling | 2.0 | 2.3 × 10−21 | 568 |
GDF15 | TGF-β superfamily member | 2.8 | 7.1 × 10−21 | 611 |
PRDM1 | Represses IFN-β gene expression | 1.6 | 7.8 × 10−21 | 617 |
MICB | Activates cytolytic response of NK and CD8 T cells | 1.5 | 5.6 × 10−20 | 696 |
IL-6R | Key immunosuppressive pathway | 1.7 | 6.4 × 10−20 | 707 |
Immunosuppressors | ||||
TIMP1 | Regulates resistance to infection | 4.3 | 4.2 × 10−45 | 12 |
ANXA1 | Enhances/inhibits adaptive immunity | 3.3 | 3.4 × 10−41 | 23 |
TNFRSF1A | Tumor necrosis superfamily protein | 2.2 | 5.4 × 10−40 | 31 |
MR1 | Antigen presentation function | 2.1 | 2.2 × 10−37 | 45 |
KLRC2 | Encodes NKG2-C (expressed on NK cells) | −10.9 | 1.3 × 10−36 | 50 |
SERPING1 | Complement activation | 2.3 | 1.1 × 10−34 | 71 |
LTBP1 | Activates TGF-β | 3.2 | 1.0 × 10−32 | 100 |
SWAP70 | Mediates IgE responses/antibody switch | 2.1 | 2.1 × 10−32 | 109 |
SATB1 | Silencing mimics the effects of IFN-γ treatment | −2.4 | 4.5 × 10−32 | 114 |
TGFBR2 | TGF-β2 receptor | 1.7 | 6.0 × 10−32 | 120 |
SHC1 | Inhibits adaptive immune responses | 1.9 | 1.0 × 10−31 | 125 |
LGALS1 | Galectin 1—deactivates M1 macrophages | 2.3 | 1.1 × 10−31 | 127 |
RELB | Interacts with NF-κβ2 | 2.0 | 4.1 × 10−31 | 135 |
DLL1 | T-cell activation | −2.0 | 4.4 × 10−30 | 142 |
BCL3 | Regulates NF-κβ | 2.7 | 6.3 × 10−30 | 146 |
RAB27A | Granulocyte exocytosis | 2.2 | 1.3 × 10−29 | 157 |
NOD1 | Innate immunity | 1.8 | 1.4 × 10−29 | 160 |
KLRC1 | Recognition of MHC class I molecules by NK cells | −3.0 | 1.5 × 10−28 | 161 |
CHI3L1 | YKL-40 immunosuppressor in gliomas | 3.0 | 2.2 × 10−29 | 170 |
ICAM1 | Transmigration receptor for immune cells | 2.7 | 4.9 × 10−29 | 182 |
Abbreviation: IgE, immunoglobulin E.
Although 4 subtypes of glioblastoma were identified by Verhaak and colleagues (21) and 3 subsets by Phillips and colleagues (5), the proneural and mesenchymal subsets identified using distinct methodologies and sample sets seem to be the most robust and concordant (27, 28). Because the greatest immunologic diversity exists between the proneural and mesenchymal subset, we compared the immune genes of interest between these subsets. Many of the top tiers of preferentially expressed immune genes in the mesenchymal glioblastoma subtype relative to the proneural subtype have been shown previously to have a biologic role in glioblastoma (Table 1).
Immunosuppression predominates in mesenchymal glioblastomas
To determine if the list of genes generated from IPA and Uniprot was artificial or if it contained candidates that might correlate to the growth and maintenance of mesenchymal glioblastoma, we generated a list of known immunosuppressive genes associated with glioma biology and then surveyed the TCGA database for overexpression within subsets. To ascertain if there was a selective enrichment of the immunosuppressive genes among the overexpressed genes in the mesenchymal subset or if this was merely by chance, we compared the immunosuppressive gene set with the overall set of differentially expressed genes from the TCGA data set using the χ2 and hypergeometric tests and found selective enrichment of immunosuppressive genes within those overexpressed in the mesenchymal subset relative to the other subsets (P < 0.05), similar to the findings described earlier.
We parsed the TCGA dataset using the cBio cancer genomics portal analysis tool to show that the proportion of patients with glioblastoma in the mesenchymal subset had significantly higher (relative to the proneural subtype) mRNA expression of previously well-characterized glioma-mediated immunosuppressive cytokines such as galectin-3 (χ2 test; P < 0.0001), IL-10 (P < 0.0001), IL-23 (P = 0.0012), TGF-β (P < 0.0001), and the immune activation inhibitor PD-L1 (P < 0.0001) relative to expression in the proneural subtype (Table 2). Many of these immunosuppressive genes are associated with the monocyte and macrophage family including: colony-stimulating factor 1 (CSF-1; P < 0.0001), a cytokine that controls the production, differentiation, and function of macrophages; CCL2 (P < 0.0001) and CCL-22 (P = 0.0008), chemokines that attract monocytes; CD163 (P < 0.0001), a marker of immunosuppressive M2 macrophages; CD204 (P < 0.0001), a macrophage scavenger receptor; and macrophage inhibitory cytokine-1 (MIC-1; P < 0.0001). CD47, a block used by solid cancers to prevent phagocytosis (29), was not preferentially expressed in the glioblastoma subtype. Arginase, produced by myeloid-derived suppressor cells and an inhibitor of T-cell responses in the glioma microenvironment (30), is also enriched in the mesenchymal subset (P = 0.0248). Genes in the IL-6/STAT3 immunosuppressive signaling axis was preferentially expressed in the mesenchymal subset; they include: gp130 (P < 0.0001), part of the IL-6 receptor family; IL-6 (P < 0.0001); STAT3 (P = 0.0004); Pim1 (P < 0.0001), a STAT3 downstream-regulated proto-oncogene; and SOCS (P < 0.0001), a negative regulator of cytokine signaling.
Enrichment of immunosuppression in the mesenchymal glioblastoma subset
. | . | Number of cases; % of cases with mRNA overexpression . | ||||
---|---|---|---|---|---|---|
. | . | . | Proneural . | Mesenchymal . | Classical . | Neural . |
Immunosuppressor/gene . | PMID reference . | n = 141 . | n = 160 . | n = 147 . | n = 96 . | |
Immunosuppressive cytokines and checkpoints | Galectin-3/LGALS3 | 22672152 | 2; 1 | 28; 18 | 13; 9 | 6; 6 |
VEGF/VEGFA | 20549821 | 16; 11 | 26; 16 | 32; 22 | 3; 3 | |
IL-10/IL-10 | 22981868 | 4; 3 | 39; 24 | 5; 3 | 13; 14 | |
IL-23/IL-23A | 20404142 | 4; 3 | 21; 13 | 12; 8 | 5; 5 | |
TGF-β/TGFB1 | 9597127 | 5; 4 | 50; 31 | 14; 10 | 2; 2 | |
PD-1/SPATA2 | 11209085 | 28; 20 | 14; 9 | 58; 39 | 27; 28 | |
PD-L1/PD-L1 | 22180678 | 0; 0 | 25; 16 | 14; 10 | 5; 5 | |
CTLA-4/CTLA-4 | 20578982 | 12; 9 | 30; 19 | 8; 5 | 11; 11 | |
Tumor-supportive macrophage chemotactic and skewing molecules | CSF-1/CSF | 14709771 | 3; 2 | 30; 19 | 4, 3 | 1, 1 |
CCL2/CCL2 | 22162712 | 5; 4 | 53; 33 | 9; 6 | 7; 7 | |
CCL-22/CCL-22 | 20518016 | 10; 7 | 33; 21 | 17; 12 | 12; 13 | |
CD163/CD163 | 15478309 | 8; 6 | 60; 38 | 2; 1 | 11; 11 | |
CD204/MSR1 | 22083206 | 5; 4 | 53; 33 | 3; 2 | 8; 8 | |
MIC-1/GDF15 | 20534737 | 7; 5 | 43; 27 | 25; 17 | 14; 15 | |
Arginase/ARG1 | 20643302 | 9; 6 | 23; 14 | 16; 11 | 22; 23 | |
CD47/CD47 | 19666525 | 15; 11 | 30; 19 | 10; 7 | 19; 20 | |
Immunosuppressive signaling pathways | IL-6/IL-6 | 23248265 | 32; 23 | 83; 52 | 16; 11 | 15; 16 |
gp130/IL-6ST | 20610800 | 0; 0 | 25; 16 | 17; 12 | 8; 8 | |
Jak2 | 22684105 | 6; 4 | 22; 14 | 9; 6 | 11; 11 | |
STAT3/STAT3 | 20053772 | 8; 6 | 31; 19 | 26; 18 | 0; 0 | |
Pim-1/PIM1 | 22384197 | 4; 3 | 44; 28 | 13; 9 | 6; 6 | |
SOCS3/SOCS3 | 10837055 | 5; 4 | 36; 23 | 10; 7 | 3; 3 | |
STAT5A/STAT5A | 12835478 | 4; 3 | 48; 30 | 10; 7 | 2; 2 | |
Markers of Tregs | CD4/CD4 | 20605226 | 5; 4 | 57; 36 | 0; 0 | 9; 9 |
CD278/ICOS | 23026134 | 8; 6 | 23; 14 | 9; 6 | 9; 9 | |
IDO/IDO1 | 22932670 | 6; 4 | 25; 16 | 14; 10 | 4; 4 | |
FoxP3/FOXP3 | 20068105 | 0; 0 | 0; 0 | 0; 0 | 0; 0 |
. | . | Number of cases; % of cases with mRNA overexpression . | ||||
---|---|---|---|---|---|---|
. | . | . | Proneural . | Mesenchymal . | Classical . | Neural . |
Immunosuppressor/gene . | PMID reference . | n = 141 . | n = 160 . | n = 147 . | n = 96 . | |
Immunosuppressive cytokines and checkpoints | Galectin-3/LGALS3 | 22672152 | 2; 1 | 28; 18 | 13; 9 | 6; 6 |
VEGF/VEGFA | 20549821 | 16; 11 | 26; 16 | 32; 22 | 3; 3 | |
IL-10/IL-10 | 22981868 | 4; 3 | 39; 24 | 5; 3 | 13; 14 | |
IL-23/IL-23A | 20404142 | 4; 3 | 21; 13 | 12; 8 | 5; 5 | |
TGF-β/TGFB1 | 9597127 | 5; 4 | 50; 31 | 14; 10 | 2; 2 | |
PD-1/SPATA2 | 11209085 | 28; 20 | 14; 9 | 58; 39 | 27; 28 | |
PD-L1/PD-L1 | 22180678 | 0; 0 | 25; 16 | 14; 10 | 5; 5 | |
CTLA-4/CTLA-4 | 20578982 | 12; 9 | 30; 19 | 8; 5 | 11; 11 | |
Tumor-supportive macrophage chemotactic and skewing molecules | CSF-1/CSF | 14709771 | 3; 2 | 30; 19 | 4, 3 | 1, 1 |
CCL2/CCL2 | 22162712 | 5; 4 | 53; 33 | 9; 6 | 7; 7 | |
CCL-22/CCL-22 | 20518016 | 10; 7 | 33; 21 | 17; 12 | 12; 13 | |
CD163/CD163 | 15478309 | 8; 6 | 60; 38 | 2; 1 | 11; 11 | |
CD204/MSR1 | 22083206 | 5; 4 | 53; 33 | 3; 2 | 8; 8 | |
MIC-1/GDF15 | 20534737 | 7; 5 | 43; 27 | 25; 17 | 14; 15 | |
Arginase/ARG1 | 20643302 | 9; 6 | 23; 14 | 16; 11 | 22; 23 | |
CD47/CD47 | 19666525 | 15; 11 | 30; 19 | 10; 7 | 19; 20 | |
Immunosuppressive signaling pathways | IL-6/IL-6 | 23248265 | 32; 23 | 83; 52 | 16; 11 | 15; 16 |
gp130/IL-6ST | 20610800 | 0; 0 | 25; 16 | 17; 12 | 8; 8 | |
Jak2 | 22684105 | 6; 4 | 22; 14 | 9; 6 | 11; 11 | |
STAT3/STAT3 | 20053772 | 8; 6 | 31; 19 | 26; 18 | 0; 0 | |
Pim-1/PIM1 | 22384197 | 4; 3 | 44; 28 | 13; 9 | 6; 6 | |
SOCS3/SOCS3 | 10837055 | 5; 4 | 36; 23 | 10; 7 | 3; 3 | |
STAT5A/STAT5A | 12835478 | 4; 3 | 48; 30 | 10; 7 | 2; 2 | |
Markers of Tregs | CD4/CD4 | 20605226 | 5; 4 | 57; 36 | 0; 0 | 9; 9 |
CD278/ICOS | 23026134 | 8; 6 | 23; 14 | 9; 6 | 9; 9 | |
IDO/IDO1 | 22932670 | 6; 4 | 25; 16 | 14; 10 | 4; 4 | |
FoxP3/FOXP3 | 20068105 | 0; 0 | 0; 0 | 0; 0 | 0; 0 |
NOTE: Red denotes preferential statistically significant enrichment (P < 0.05) of mRNA expression in the comparison of the mesenchymal vs. proneural subset.
Proinflammatory responses are more frequent in the mesenchymal glioblastoma subtype
We generated a list of key immune genes that reflect immune stimulation/effector responses of both adaptive and innate immunity. Again, we found a selective enrichment of the overexpressed immune effector genes in the mesenchymal subset relative to all other subsets (P < 0.05). Focusing on the differences between immune effector genes expressed between the mesenchymal and proneural subset, genes encoding key receptor interactions for T-cell activation, including CD3 (P < 0.0001), CD40 (P < 0.0001), CD80 (P ≤ 0.0001), CD86 (P < 0.0001), MHC HLA-B (P = 0.0001), MHC HLA-DRA (P < 0.0001), MHC HLA-DQA1 (P < 0.0001), and MHC HLA-DPB1 (P < 0.0001) and downstream signaling pathways reflective of T-cell activation, effector function, and immune-activating transcription factors, such as STAT4 (P = 0.0034), DNAM-1 (P < 0.0001), and IRF7 (P = 0.0088), were enhanced in the proportion of patients with glioblastoma with the mesenchymal phenotype (Table 3). Although there was a trend toward preferential expression of the adaptive proinflammatory cytokines in the proportion of patients with glioblastoma with the mesenchymal subtype, only IL-1 (P = 0.0007), IL-15 (P < 0.0001), and IL-7 (P < 0.0001) were significantly expressed at the mRNA level.
Enrichment of proinflammatory responses in the mesenchymal glioblastoma subset
. | . | Number of cases; % of cases with mRNA overexpression . | ||||
---|---|---|---|---|---|---|
. | . | . | Proneural . | Mesenchymal . | Classical . | Neural . |
Immune effector/gene . | PMID reference . | n = 141 . | n = 160 . | n = 147 . | n = 96 . | |
Proinflammatory cytokines | IFN-γ/IFNG | 21705455 | 17; 12 | 27; 17 | 17; 12 | 12; 13 |
IL-1/IL-1A | 17596594 | 11; 8 | 35; 22 | 6; 4 | 17; 18 | |
IL-2/IL-2 | 20227890 | 14; 10 | 14; 9 | 12; 8 | 11; 11 | |
IL-4/IL-4 | 17389011 | 16; 11 | 29; 18 | 19; 13 | 16: 17 | |
IL-7/IL-7 | 8473051 | 5; 4 | 35; 22 | 20; 14 | 15; 16 | |
IL-12/IL-12A | 22515181 | 10; 7 | 13; 8 | 35; 24 | 12; 13 | |
IL-15/IL-15 | 22032984 | 7; 5 | 42; 26 | 4; 3 | 17; 18 | |
TNF-α/TNF | 22075379 | 19; 13 | 21; 13 | 5; 3 | 13; 14 | |
T-cell effector markers | CD3/CD3D | 23009144 | 8; 6 | 40; 25 | 3; 2 | 9; 9 |
CD8/CD8A | 18698034 | 12; 9 | 29; 18 | 12; 8 | 9; 9 | |
CD8/CD8B | 18698034 | 3; 3 | 22; 14 | 9; 6 | 6; 6 | |
CD80/CD80 | 22798683 | 9; 6 | 36; 23 | 8; 5 | 13; 14 | |
CD86/CD86 | 20380573 | 4; 3 | 53; 33 | 2; 1 | 19; 20 | |
CD40/CD40 | 22903346 | 7; 3 | 38; 24 | 6; 4 | 10; 10 | |
MHC/HLA-A | 23108078 | 12; 9 | 21; 13 | 25; 17 | 10; 10 | |
MHC/HLA-B | 22426959 | 5; 4 | 28; 18 | 30; 20 | 9; 9 | |
MHC/HLA-C | 17054674 | 10; 7 | 27; 17 | 25; 17 | 12; 13 | |
MHC/HLA-DRA | 15688398 | 5; 4 | 42; 26 | 9; 6 | 11; 11 | |
MHC/HLA-DQA1 | 17016821 | 13; 9 | 46; 29 | 17; 12 | 15; 16 | |
MHC/HLA-DPB1 | 21716314 | 5; 4 | 45; 28 | 10; 17 | 8; 8 | |
β2-Microglobulin/B2M | 22353804 | 8; 6 | 23; 14 | 21; 14 | 8; 8 | |
Immune effector signaling pathways | NF-κB/RELA | 16724054 | 10; 7 | 17; 11 | 17; 12 | 3; 3 |
STAT1/STAT1 | 22805310 | 11; 8 | 20; 13 | 30; 20 | 10; 10 | |
IRAK/IRAK1 | 17265049 | 0 | 0 | 0 | 0 | |
STAT4/STAT4 | 22121102 | 9; 6 | 28; 18 | 5; 3 | 15; 16 | |
T-bet/TBX21 | 22186896 | 12; 9 | 20; 13 | 17; 12 | 6; 6 | |
DNAM-1/CD226 | 16015041 | 8; 6 | 42; 26 | 4; 3 | 11; 11 | |
IRF7/IRF7 | 22295238 | 8; 6 | 24; 15 | 16; 11 | 17; 18 | |
Markers of innate immunity | NKp46/NCR1 | 22308311 | 14; 10 | 22; 14 | 22; 15 | 7; 7 |
NKG2D/KLRK1 | 22530569 | 42; 30 | 12; 8 | 6; 4 | 2; 2 | |
NKp30/NCR3 | 21877119 | 18; 13 | 24; 15 | 15; 10 | 11; 11 | |
NKp44/NCR2 | 15728472 | 12; 9 | 16; 10 | 18; 12 | 8; 8 | |
KLRD1/KLRD1 | 15805295 | 16; 11 | 27; 17 | 11; 7 | 16; 17 | |
NK/CD244 | 19638467 | 8; 6 | 42; 26 | 4; 3 | 11; 11 | |
TLR2/TLR2 | 23324344 | 1; 1 | 56; 35 | 3; 2 | 13; 14 | |
TLR3/TLR3 | 23197495 | 7; 5 | 34; 21 | 16; 11 | 22; 23 | |
TLR4/TLR4 | 21129170 | 17; 12 | 26; 16 | 8; 5 | 17; 18 | |
TLR9/TLR9 | 22169598 | 13; 9 | 23; 14 | 23; 16 | 6; 6 | |
Myeloid markers | TREM1/TREM1 | 22719066 | 12; 9 | 63; 39 | 9; 6 | 2; 2 |
TREM2/TREM2 | 4; 3 | 27; 17 | 6; 4 | 31; 32 | ||
MIF/MIF | 21773885 | 13; 9 | 25; 16 | 8; 5 | 14; 15 |
. | . | Number of cases; % of cases with mRNA overexpression . | ||||
---|---|---|---|---|---|---|
. | . | . | Proneural . | Mesenchymal . | Classical . | Neural . |
Immune effector/gene . | PMID reference . | n = 141 . | n = 160 . | n = 147 . | n = 96 . | |
Proinflammatory cytokines | IFN-γ/IFNG | 21705455 | 17; 12 | 27; 17 | 17; 12 | 12; 13 |
IL-1/IL-1A | 17596594 | 11; 8 | 35; 22 | 6; 4 | 17; 18 | |
IL-2/IL-2 | 20227890 | 14; 10 | 14; 9 | 12; 8 | 11; 11 | |
IL-4/IL-4 | 17389011 | 16; 11 | 29; 18 | 19; 13 | 16: 17 | |
IL-7/IL-7 | 8473051 | 5; 4 | 35; 22 | 20; 14 | 15; 16 | |
IL-12/IL-12A | 22515181 | 10; 7 | 13; 8 | 35; 24 | 12; 13 | |
IL-15/IL-15 | 22032984 | 7; 5 | 42; 26 | 4; 3 | 17; 18 | |
TNF-α/TNF | 22075379 | 19; 13 | 21; 13 | 5; 3 | 13; 14 | |
T-cell effector markers | CD3/CD3D | 23009144 | 8; 6 | 40; 25 | 3; 2 | 9; 9 |
CD8/CD8A | 18698034 | 12; 9 | 29; 18 | 12; 8 | 9; 9 | |
CD8/CD8B | 18698034 | 3; 3 | 22; 14 | 9; 6 | 6; 6 | |
CD80/CD80 | 22798683 | 9; 6 | 36; 23 | 8; 5 | 13; 14 | |
CD86/CD86 | 20380573 | 4; 3 | 53; 33 | 2; 1 | 19; 20 | |
CD40/CD40 | 22903346 | 7; 3 | 38; 24 | 6; 4 | 10; 10 | |
MHC/HLA-A | 23108078 | 12; 9 | 21; 13 | 25; 17 | 10; 10 | |
MHC/HLA-B | 22426959 | 5; 4 | 28; 18 | 30; 20 | 9; 9 | |
MHC/HLA-C | 17054674 | 10; 7 | 27; 17 | 25; 17 | 12; 13 | |
MHC/HLA-DRA | 15688398 | 5; 4 | 42; 26 | 9; 6 | 11; 11 | |
MHC/HLA-DQA1 | 17016821 | 13; 9 | 46; 29 | 17; 12 | 15; 16 | |
MHC/HLA-DPB1 | 21716314 | 5; 4 | 45; 28 | 10; 17 | 8; 8 | |
β2-Microglobulin/B2M | 22353804 | 8; 6 | 23; 14 | 21; 14 | 8; 8 | |
Immune effector signaling pathways | NF-κB/RELA | 16724054 | 10; 7 | 17; 11 | 17; 12 | 3; 3 |
STAT1/STAT1 | 22805310 | 11; 8 | 20; 13 | 30; 20 | 10; 10 | |
IRAK/IRAK1 | 17265049 | 0 | 0 | 0 | 0 | |
STAT4/STAT4 | 22121102 | 9; 6 | 28; 18 | 5; 3 | 15; 16 | |
T-bet/TBX21 | 22186896 | 12; 9 | 20; 13 | 17; 12 | 6; 6 | |
DNAM-1/CD226 | 16015041 | 8; 6 | 42; 26 | 4; 3 | 11; 11 | |
IRF7/IRF7 | 22295238 | 8; 6 | 24; 15 | 16; 11 | 17; 18 | |
Markers of innate immunity | NKp46/NCR1 | 22308311 | 14; 10 | 22; 14 | 22; 15 | 7; 7 |
NKG2D/KLRK1 | 22530569 | 42; 30 | 12; 8 | 6; 4 | 2; 2 | |
NKp30/NCR3 | 21877119 | 18; 13 | 24; 15 | 15; 10 | 11; 11 | |
NKp44/NCR2 | 15728472 | 12; 9 | 16; 10 | 18; 12 | 8; 8 | |
KLRD1/KLRD1 | 15805295 | 16; 11 | 27; 17 | 11; 7 | 16; 17 | |
NK/CD244 | 19638467 | 8; 6 | 42; 26 | 4; 3 | 11; 11 | |
TLR2/TLR2 | 23324344 | 1; 1 | 56; 35 | 3; 2 | 13; 14 | |
TLR3/TLR3 | 23197495 | 7; 5 | 34; 21 | 16; 11 | 22; 23 | |
TLR4/TLR4 | 21129170 | 17; 12 | 26; 16 | 8; 5 | 17; 18 | |
TLR9/TLR9 | 22169598 | 13; 9 | 23; 14 | 23; 16 | 6; 6 | |
Myeloid markers | TREM1/TREM1 | 22719066 | 12; 9 | 63; 39 | 9; 6 | 2; 2 |
TREM2/TREM2 | 4; 3 | 27; 17 | 6; 4 | 31; 32 | ||
MIF/MIF | 21773885 | 13; 9 | 25; 16 | 8; 5 | 14; 15 |
NOTE: Red denotes preferential statistically significant enrichment (P < 0.05) of mRNA expression in the comparison of the mesenchymal vs. proneural subset.
In contrast, genes associated with innate immunity, specifically natural killer (NK) cell markers such as NKp46, NKG2D, NKp30, NKp44, and KLRD1 were more uniformly distributed across the subtypes, with the exception of the NK-activating receptor CD244 (P < 0.0001) that was preferentially expressed in the proneural subtype. Other genes involved in innate immune response, such as Toll-like receptor (TLR) 2 (P < 0.0001) and 3 (P < 0.0001) were preferentially overexpressed in patients with glioblastoma with the mesenchymal subset. The triggering receptors expressed on myeloid cells (TREM) 1 (P < 0.0001) and 2 (P < 0.0001) that enhance monocyte/macrophage–inflammatory responses were also enhanced in the mesenchymal subset. Cumulatively, these data suggest that there is a preferential distribution of both proinflammatory and immunosuppressive genes within the mesenchymal subset.
Glioma antigens segregate differentially within glioblastoma subtypes
We postulated that the mesenchymal subset may have greater incidence of tumor-associated and -specific antigens that could contribute to the increased propensity of the immune effector genes present within the mesenchymal subset. Furthermore, there would be a predicted enhancement in immunosuppression to counteract the antitumor effector responses. Thus, to evaluate if tumor antigens are selectively enriched in the mesenchymal subset, a list of known glioma antigens was compiled, and the frequency of each mRNA overexpression was determined within the various glioblastoma subtypes (Table 4). Although there were predilections of specific tumor antigens for glioblastoma subtypes, there was not a preferential enrichment of antigens found within the mesenchymal subset. Specifically, the EGFR family of antigens including EGFR and ERBB2 were frequently found to be overexpressed at 80% and 24%, respectively, within the classical subtype, consistent with previous reports (21). Within the neural subtype, EGFR (54%) was frequently expressed, whereas within the proneural subtype survivin (32%) and Sart-1 (25%) were more commonly expressed. SART-2 (36%) was the most commonly expressed antigen within the mesenchymal subset. Thus, preferential antigen enrichment in the mesenchymal subset does not seem to be the underlying etiology for enrichment of the immune activator and suppressor genes within this subset.
Antigenic diversity among glioblastoma subset
. | Number of cases; % of cases with mRNA overexpression . | |||
---|---|---|---|---|
. | Proneural . | Mesenchymal . | Classical . | Neural . |
Glioma antigen/gene . | n = 141 . | n = 160 . | n = 147 . | n = 96 . |
EGFR/EGFR | 32; 23 | 54; 34 | 118; 80 | 52; 54 |
Her2/ERBB2 | 4; 3 | 25; 16 | 35; 24 | 7; 7 |
Survivin/BIRC5 | 43; 32 | 22; 14 | 6; 4 | 19; 20 |
Nucleolin/NCL | 30; 21 | 22; 14 | 23; 16 | 2; 2 |
Epha2/EPHA2 | 2; 1 | 27; 17 | 28; 19 | 4; 4 |
Telomerase/TERT | 6; 4 | 14; 9 | 18; 12 | 5; 4 |
B-cyclin/CCNB1 | 34; 24 | 20; 13 | 7; 5 | 12; 13 |
Sart-1/SART1 | 35; 25 | 9; 6 | 18; 12 | 2; 2 |
Sart-2/DSE | 2; 1 | 57; 36 | 2; 1 | 3; 3 |
Sart-3/SART3 | 18; 13 | 19; 12 | 20; 14 | 3; 3 |
Aim-2/AIM2 | 32; 23 | 38; 24 | 5; 3 | 10; 10 |
Trp-1/TYRP1 | 13; 9 | 23; 14 | 8; 5 | 15; 16 |
Tyrosinase/TYR | 11; 8 | 15; 9 | 7; 5 | 11; 11 |
GnT-V/MGAT5 | 5; 4 | 8; 5 | 5; 3 | 1; 1 |
GP100/PMEL | 7; 5 | 18; 11 | 7; 5 | 6; 6 |
Mart-1/MLANA | 11; 8 | 24; 15 | 11; 7 | 10; 10 |
Mage-1/MAGEA1 | 0 | 0 | 0 | 0 |
Gage-1/GAGE1 | 0 | 0 | 0 | 0 |
. | Number of cases; % of cases with mRNA overexpression . | |||
---|---|---|---|---|
. | Proneural . | Mesenchymal . | Classical . | Neural . |
Glioma antigen/gene . | n = 141 . | n = 160 . | n = 147 . | n = 96 . |
EGFR/EGFR | 32; 23 | 54; 34 | 118; 80 | 52; 54 |
Her2/ERBB2 | 4; 3 | 25; 16 | 35; 24 | 7; 7 |
Survivin/BIRC5 | 43; 32 | 22; 14 | 6; 4 | 19; 20 |
Nucleolin/NCL | 30; 21 | 22; 14 | 23; 16 | 2; 2 |
Epha2/EPHA2 | 2; 1 | 27; 17 | 28; 19 | 4; 4 |
Telomerase/TERT | 6; 4 | 14; 9 | 18; 12 | 5; 4 |
B-cyclin/CCNB1 | 34; 24 | 20; 13 | 7; 5 | 12; 13 |
Sart-1/SART1 | 35; 25 | 9; 6 | 18; 12 | 2; 2 |
Sart-2/DSE | 2; 1 | 57; 36 | 2; 1 | 3; 3 |
Sart-3/SART3 | 18; 13 | 19; 12 | 20; 14 | 3; 3 |
Aim-2/AIM2 | 32; 23 | 38; 24 | 5; 3 | 10; 10 |
Trp-1/TYRP1 | 13; 9 | 23; 14 | 8; 5 | 15; 16 |
Tyrosinase/TYR | 11; 8 | 15; 9 | 7; 5 | 11; 11 |
GnT-V/MGAT5 | 5; 4 | 8; 5 | 5; 3 | 1; 1 |
GP100/PMEL | 7; 5 | 18; 11 | 7; 5 | 6; 6 |
Mart-1/MLANA | 11; 8 | 24; 15 | 11; 7 | 10; 10 |
Mage-1/MAGEA1 | 0 | 0 | 0 | 0 |
Gage-1/GAGE1 | 0 | 0 | 0 | 0 |
NOTE: Red denotes preferential statistically significant enrichment (P < 0.05) of mRNA expression in the comparison of the mesenchymal vs. proneural subset.
Discussion
The TCGA database consists of an unselected cellular population that includes glioma-infiltrating immune cells within the genetic composition. The analysis of the mRNA overexpression of immune genes may reflect either gene amplification (increased expression on the immune population) or a relative increase in a designated infiltrating immune population. The interaction between the immune system and molecular subtypes of glioma remains largely unstudied. Although a previous report identified an immune signature that is prognostic for survival in patients with glioblastoma, especially within the proneural subtype (31), this is the first report to show the concordant association of proinflammatory and immunosuppression in the mesenchymal subset. Specifically, patients with glioblastoma that have a preexisting induced anti-glioma effector response and an actionable immunosuppressive target (i.e., the mesenchymal subset) may be “immune reactive” and therefore particularly amenable to immunotherapeutic approaches including those targeting immune checkpoints. This is further supported by a retrospective analysis of patients with glioblastoma receiving dendritic cell immunotherapy who were more likely to have a mesenchymal subtype glioblastoma (20). These data indicate that the glioblastoma subtype may be a confounding variable in therapeutic response analysis and should be considered in stratification. In addition, our data indicate that glioblastoma subtype may influence the interpretation of posttreatment analysis of immune infiltration. Specifically, because glioblastomas are more likely to transition to the mesenchymal subtype upon recurrence (5), the intratumoral immune analysis after immunotherapy may be more reflective of the biology of the underlying subtype rather than a direct immunotherapy-induced response.
The proinflammatory and immunosuppressive gene sets derived from online-curated databases identified a set of common genes (e.g., TIMP1). In some cases, the role of these genes is contextual based on the disease state (i.e., autoimmunity vs. malignancy) or model system studied. For example, ANXA1 has been shown to both inhibit (32, 33) and enhance inflammation (34, 35). In other instances, that attribution to a particular category simply seemed erroneous (i.e., LTBP1 that activates the immunosuppressive cytokine TGF-β or TGFBR2 that encodes the TGF-β2 receptor). Thus, we selected a proinflammatory and immunosuppressive list of genes based on the documented roles of these genes in the setting of malignancy with a special emphasis on those operational in glioblastoma. This defined gene list also showed a preferential enrichment of both proinflammatory and immunosuppressive genes within the mesenchymal subtype. Although we had concerns about the selection of genes by the online-curated databases, they did provide novel insights into immune genes that had, thus far, not been appreciated as playing a role in glioblastoma immune biology. Furthermore, these databases revealed the marked diversity of the preferential enrichment of proinflammatory or immunosuppressive genes in the mesenchymal subtype relative to the others.
In many instances the findings of the immune genes associating with a particular glioblastoma subtype were consistent with and validated previous observations. For example, as STAT3 has been previously shown to be a key molecular hub driving the mesenchymal transformation (6), STAT3 and related downstream targets, such as Pim-1 and VEGF, were found to be enriched in the mesenchymal subtype relative to the other subsets. Furthermore, multiple monocyte genes were enriched within the mesenchymal subtype as reflected in the expression of CCL2 (a known monocyte chemotactic protein), CD163 (a marker of cells of the monocyte/macrophage lineage), and CD204 (a marker of macrophage scavenger receptors). These data are consistent with our findings in genetically engineered murine models in which the mesenchymal transition was shown to correlate with increased macrophage infiltration (7, 14). In addition, the proinflammatory T-cell cytokine IL-2 was minimally enriched in any subtype, consistent with our previous reports showing that while T cells may be active in the periphery, upon encountering the local tumor microenvironment, this immune response is markedly downregulated (3). Interestingly, markers of Tregs such as CD4, ICOS, IDO1, and CTLA-4 seemed to be preferentially enriched in the mesenchymal subset. Because robust immune effector response is also present in the mesenchymal subset, these patients may be specifically predisposed to respond to immunotherapeutics targeting the Treg population. However, while these markers have been shown to be important for Treg function and are associated with immune inhibition, they do not independently define the Treg population and mechanistic conclusions cannot be drawn from this type of analysis.
Although CCL22 (a chemokine that attracts Tregs; ref. 36), IDO and TGF-β (inducible Treg factors; ref. 37), and ICOS (critical for the functional stability of Tregs; ref. 38) were present in gliomas, no differential expression of FoxP3 was observed in glioblastoma subtypes. This seems inconsistent with the immunohistochemical data showing a great deal of heterogeneity in the presence of Tregs within glioblastomas (2). However, because the TCGA database is generated by cDNA microarray analysis, which represents gene regulation at the transcriptional level, immune genes that are regulated at the posttranscriptional, translational, or posttranslational level may not appear through this screening. A specific example of this is FoxP3, which is regulated either at the mRNA transcription level or at its protein stability level, which is degraded through the proteasome rapidly (39). Thus, the discrepancy between the immunohistochemical data and the TCGA data can not only be explained on the basis of FoxP3 mRNA and protein stability, but also shows a limitation of solely relying on analysis of mRNA data.
Another limitation of the current study is that the analysis does not directly reflect the systemic immune status or functional status of these immune responses. Some of the immunosuppressive targets may be markedly enhanced within the tumor microenvironment, with marginal elevations systemically, and agents without significant glioma penetration may fail to show a correlation of treatment response with tumor expression levels. Furthermore, the analysis of the frequency of overexpression is relative to the mean expression levels of a marker across all gliomas, and thus more clinical responders may be identified within any given phenotype because the threshold of minimal expression necessary to result in a therapeutic response has not yet been characterized. In addition, the mRNA expression levels analyzed in our study may not correspond to the protein expression levels especially in the case of B7-H1, which has been shown to be posttranscriptionally regulated by cytokines (40). Finally, although there is an association of immune effector and suppressor genes within the mesenchymal subtype, this has not yet been immunologically functionally validated. Not previously reported, but as part of a previous analysis, we found that 26% of newly diagnosed glioblastoma tumors have very little immune infiltration, whereas approximately 50% have both effector (CD8+ T cells) and suppressor (FoxP3+) immune responses in the glioblastoma (2). The correlation of these immunohistochemical findings with glioblastoma subtype is an area of future investigation along with the association of observed immune gene signatures and patient outcome.
Although the mesenchymal glioblastoma subtype may be “immunologically reactive,” as defined by the preexistence of immune effector responses and possessing the therapeutic targets of immunosuppressive modulators, it is unlikely that this is the sole subtype that can potentially benefit from immunotherapeutic strategies. We have previously shown that a peptide vaccine targeting the EGFR variant III (EGFRvIII) in newly diagnosed patients with glioblastoma significantly increased median survival time to more than 26 months (17, 18). On the basis of the TCGA data, it is likely that most of these patients have glioblastomas within the classical subtype because EGFRvIII expression almost always coincides with EGFR amplification (24, 41). This indicates that with specific immunotherapeutic strategies, patients who have distinct glioblastoma subtypes may preferentially benefit.
In conclusion, the analysis of these selected immune effector and suppressor genes may provide a way to programmatically prioritize different types of potentially competing immunotherapeutic strategies, and to identify subsets of patients that may respond to a particular strategy and thus selectively enrich for potential responders during early or small-scale clinical trials.
Disclosure of Potential Conflicts of Interest
M. Gilbert is a consultant/advisory board member of Merck, Genentech, EMD Serono, and Novartis. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: G. Rao, A.B. Heimberger
Development of methodology: T. Doucette, G. Rao, A. Rao, L. Shen, K. Aldape, J. Wei, A.B. Heimberger
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): G. Rao, L. Shen, K. Aldape, K. Dziurzynski, A.B. Heimberger
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): T. Doucette, G. Rao, A. Rao, L. Shen, K. Aldape, M. Gilbert, A.B. Heimberger
Writing, review, and/or revision of the manuscript: T. Doucette, G. Rao, A. Rao, L. Shen, K. Aldape, K. Dziurzynski, M. Gilbert, A.B. Heimberger
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): A.B. Heimberger
Study supervision: A.B. Heimberger
Acknowledgments
We thank Audria Patrick and David M. Wildrick for editorial assistance.
Grant Support
Grant support was provided by the NIH RO1-CA1208113 (to A.B. Heimberger), 5P50 CA127001 (to T. Doucette, G. Rao, K. Aldape, and A.B. Heimberger), Dr. Marnie Rose Foundation (to A.B. Heimberger), the Cynthia and George Mitchell Foundation (to A.B. Heimberger), and the Vaughn Foundation (to A.B. Heimberger).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.