Skip to main content

Identification of 3 subpopulations of tumor-infiltrating immune cells for malignant transformation of low-grade glioma



Tumor-infiltrating immune cells (TIICs) are highly relevant to clinical outcome of glioma. However, previous studies cannot account for the diverse functions that make up the immune response in malignant transformation (MT) from low-grade glioma (LGG) to high-grade glioma (HGG).


Transcriptome level, genomic profiles and its relationship with clinical practice were obtained from TCGA and CGGA database. The “Cell type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT)” algorithm was used to estimate the fraction of 22 immune cell types. We divided the TCGA and CGGA set into an experiment set (n = 174) and a validation set (n = 74) by random number table method. Univariate and multivariate analyses were performed to evaluate the 22 TIICs’ value for MT in LGG. ROC curve was plotted to calculate area under curve (AUC) and cut-off value.


Heterogeneity between TIICs exists in both intra- and inter-groups. Several TIICs are notably associated with tumor grade, molecular subtypes and survival. T follicular helper (TFH) cells, activated NK Cells and M0 macrophages were screened out to be independent predictors for MT in LGG and formed an immune risk score (IRS) (AUC = 0.732, p < 0.001, 95% CI 0.657–0.808 cut-off value = 0.191). In addition, the IRS model was validated by validation group, Immunohistochemistry (IHC) and functional enrichment analyses.


The proposed IRS model provides promising novel signatures for predicting MT from LGG to HGG and may bring a better design of glioma immunotherapy studies in years to come.


Gliomas account for 60% of all primary and other central nervous system (CNS) tumor diagnoses, and make up ~ 80% of all malignant brain tumors [1]. The World Health Organization (WHO) classifies gliomas according to histology and molecular subtype, and grades them by the scale of I, II, III, IV. low-grade gliomas (LGG) typically range from grades I–II, while high-grade gliomas (HGG) are categorized as grades III–IV. Glioblastoma multiforme (GBM) is grade IV glioma subtype which often spontaneously appears in the CNS, but can also progress from LGG. GBM takes up half of CNS tumors, and is a fatal disease with no curable therapy [2]. Even with a comprehensive therapy, such as surgical resection, adjuvant radiotherapy, and alkylating agent temozolomide chemotherapy, patients who suffer from gliomas still have short median survival time, due to the aggressiveness of tumors, resistance to treatments, and recurrence over time [3]. In particular, patients with GBM approximately has a median survival of only 14–16 months [4]. In the past decade, studies on the anticancer immune responses for other tumors have promoted clinical advances in the limited success of conventional therapies. Meanwhile, the discovery of CNS lymphatic system has provided a new theoretical basis and opportunity for brain tumor immunotherapy [5].

Tumor-infiltrating immune cells (TIICs), whose function and composition subtly altered with the immune status of the host have been reported to be effectively targeted by drugs correlate with clinical outcome [6]. Melanoma and non-small-cell lung cancer are the two solid tumors in which immunotherapy has proved to be effective [7]. However, compared with these two tumors, glioma harbors a lower burden of somatic mutations and a more immunosuppressive tumor microenvironment [8]. Unique challenges should be overcome before immunotherapy applied to CNS. First, anatomically, the blood–brain barrier (BBB) restricts the entry of immune cells to the brain parenchyma. Also, the tumor cells themselves secrete a variety of immunosuppressive factors that influence macrophage polarization, dendritic cell (DC) maturation, regulatory T cell recruitment, inhibition of neutrophil and natural killer (NK) cell function. Previous studies have revealed that glioblastomas are heavily infiltrated with monocytes/microglia, although TIICs are relatively rare. Reports suggest that these cells account for 10–30% of viable cells within the tumor mass. They appear to be affected by tumors and have positive immunosuppressive effects. For example, Rodrigues et al. demonstrated that normal monocytes that come into contact with glioblastoma cells secrete multiple immunosuppressive factors (IF-10, TGF-β, B7-H1), have reduced phagocytic ability and induce apoptosis in activated T cells [9]. While preclinical data shows the success of immunotherapy for gliomas, the profiles of TIICs in glioma and their clinical value still remain to be explained.

Nevertheless, Immunohistochemistry and flow cytometry are the two most commonly used techniques that depend on a single marker for detecting TIICs in previous studies. Obviously, these approaches can be misleading and are not comprehensive as many marker proteins are expressed in different cell types. The “Cell type Identification By Estimating Relative Subsets Of RNA Transcripts” (CIBERSORT) employs deconvolution of bulk gene expression data and a sophisticate algorithm for in silico quantification of many immune cell types in heterogeneous samples as tumor stroma. Here, we used CIBERSORT, for the first time, to quantify the 22 TIICs subpopulations of immune response in glioma based on the patients’ gene expression profiling from TCGA and CGGA public databases in order to investigate its relationship between clinicals factors, with the final goal of developing new immunotherapeutic strategies.

Materials and methods


We examined expression data and clinical variables from the following main sources: The Cancer Genome Atlas (TCGA) dataset (,Chinese Glioma Genome Atlas (CGGA) dataset ( (up to April 10, 2019). We searched the supplements and contacted the investigators to get the missing information for samples. In TCGA dataset, we obtained mRNAseq data of 160 GBM samples and 528 LGG samples. 325 samples generated by Illumina HiSeq platform were collected form CGGA dataset, ranging from WHO grade II to grade IV. Then, expression profiles of each samples and corresponding clinical data were manually organized. Besides, 5 non-GBM patient data from the TCGA-GBM dataset, patients with any missing or insufficient data on age or survival data were excluded from subsequent processing. RNA sequencing data were firstly transformed using “voom” (variance modeling at the observational level) for the two datasets. Details of the study design are illustrated in Fig. 1 as a flowchart.

Fig. 1

Details of the study design. TCGA The Cancer Genome Atlas, CGGA Chinese Glioma Genome Atlas, CIBERSORT Cell type Identification By Estimating Relative Subsets Of RNA Transcripts

CIBERSORT estimation

The gene expression with standard annotation were uploaded to the CIBERSORT web portal (, and the algorithm was running the LM22 signature and 1000 permutations. Cases with a CIBERSORT output of p < 0.05, indicating that the inferred fractions of TIICs populations produced by CIBERSORT are accurate, were considered to be eligible for further analysis. For each sample, the final CIBERSORT output estimates were normalized to sum up to one and thus can be interpreted directly as cell fractions for comparison across different immune cell types and dataset.

Immunohistochemical detection of immune cell types

5 LGG and 5GBM tissue from 10 patients who received surgery in the First Affiliated Hospital of Nanjing Medical University (Nanjing, Jiangsu province, China) were constructed for immunohistochemistry. Specimens were all confirmed by pathological analysis as glioma. IHC was performed as described earlier, using monoclonal antibodies against CXCR5, CD4, CD68, CD11b, CD57 and CD56 (H-132; Santa Cruz Biotechnology, Santa Cruz, CA). Isotype-matched mouse monoclonal antibodies were used as negative controls. Slides were analyzed using an image analysis workstation (Spot Browser, ALPHELYS). Polychromatic high-resolution spot-images (740 × 540 pixel, 1.181 μm/pixel resolution) were obtained (200× fold magnification). The density was recorded as the number of positive cells per unit tissue surface area. For each duplicate, the mean density was used for statistical analysis.

Gene oncology (GO) and Kyoto encyclopedia of genes and genomes (KEGG)

GO was applied to determine the function of differentially expressed genes and pathway enrichment was analyzed by KEGG (

Statistical analyses

Statistical analyses were conducted using R software version 3.5.3 ( and SPSS 19.0 for windows (IBM, NY, USA). All statistical tests were two-sided and a p value < 0.05 is considered as significant.

Hierarchical clustering of immune cell proportions was conducted to compare distinct immune cell infiltration in different samples. The proportions of various TIICs were defined as a change from 0 to 1 based on our observation. R packages “Corrplot” (, “Pheatmap” ( and “Vioplot” ( also used to investigate differences in the composition of immune cells within and between groups. Wilcoxon test was used to evaluate the relationship between tumor grades, tumor molecular subtypes and TIICs. The association between TIICs and survival were analyzed using log rank test and Kaplan–Meier (K–M) curve visualized the results. Multivariable analyses were further operated to screen independently predictors. AUC and cut-off value were obtained by conducting ROC curve. “Limma” package was used to analysis the differential expressed gene, |log2FC| > 1.3219 and FDR < 0.05 were set as filters.


Composition of immune cells in LGG and GBM

“Limma” package [10] run firstly to normalize the gene expression data and to accommodate the operational requirements of CIBERSORT. Then, CIBERSORT algorithm was used to analysis the difference of immune infiltration between LGG and HGG samples in 22 subpopulations of immune cells. 269 out of the total 1008 samples from TCGA and CGGA datasets with p-value < 0.05 were included for subsequent processing, of which 81 samples were grouped into LGG cohort and 188 samples into the HGG cohort. The total value of all immune cells in each sample was set at one, Fig. 2a showed the proportion of all 22 subpopulations of immune cells in these samples (Fig. 2a). Obviously, the proportions of immune cells in glioma varied significantly between both intra- and inter-group. Resting NK cells and T cells regulatory (Tregs) exhibited a significant positive correlation, while there was a distinctive negative correlation between M0 macrophages and monocytes by average linkage clustering (Fig. 2b). Through hierarchical clustering based on the above data, we can find that TIICs, such as monocytes, M0 macrophages showed striking distribution differences in LGG and HGG (Fig. 2c). The violin plot (Fig. 2d) showed that there were marked differences in the distribution of 10 out of 22 immune cells, such as monocytes (p < 0.001), M0 macrophages (p < 0.001), activated NK cells (p < 0.01), between LGG and HGG cohorts. Taken together, these results suggest that the heterogeneity of TIICs in gliomas is evident and may play a role in the malignant progression of LGG to HGG.

Fig. 2

The landscape of immune infiltration in glioma. a The difference of immune infiltration between low- and high-grade glioma. b Correlation matrix of all 22 immune cell proportions in TCGA and CGGA datasets. c Heatmap of 22 immune cell proportions. The horizontal axis shows the clustering information of samples which were divided into two major clusters. d The distribution of same immune cells between low- and high-grade glioma. p-values show the significance of the distribution

Immune cells associated with tumor grade and molecular subtypes

Wilcoxon tests were used to evaluated the relationship between tumor grades, tumor molecular subtypes and TIICs. p-values are shown in Table 1. TIICs such as M0 Macrophages and TFH cells are positively correlated with elevated levels (Fig. 3a, b), while cell subpopulations like monocytes and activated NK cells were negatively correlated (Fig. 3c, d). Subsequently, due to the lack of molecular subtypes information in TCGA dataset, we analyzed the relationship between glioma molecular subtypes and TIICs in samples from CGGA dataset. p-values are also shown in Table 1. Differences exist in the distribution of tumor-infiltrating cells of various glioma molecular subtypes (Fig. 3e–h). These results further demonstrate that TIICs may affect the progression of glioma to some extent.

Table 1 Comparison of CIBERSORT immune cells fractions between grade and molecular subtypes of glioma
Fig. 3

Immune cells associated with tumor grade and molecular subtypes. ad Relationship between the fraction of M0 macrophages (p < 0.001), TFH cells (p = 0.001), monocytes (p < 0.001) and activated NK cells (p < 0.001) and glioma WHO grade. eh Relationship between the fraction of M0 macrophages (p < 0.001), activated Mast cells (p < 0.001), activated NK cells (p < 0.001) and monocytes (p = 0.001) and molecular subtypes of glioma

Immune cell comparison responding to the prognosis of gliomas

We then obtained clinical data from TCGA and CGGA databases, samples with a survival time less than 30 days were excluded. Then expression profiles of each samples and corresponding clinical data were manually organized. To further analyze the factors associated with patient prognosis and perform subsequent validation, we randomly divided the total sample into 70% of the experiment group and 30% of the validation group. For immune cell infiltration and corresponding survival time, we performed univariate analyses (Table 2) and nine immune cells showed significant prognostic value in both databases. Highly expressed T follicular helper cells, resting NK cells, M0 macrophages, M1 macrophages and resting Dendritic cells predicted poor overall survival, while high expression of the other 4 cell subpopulations, including plasma cells, activated NK cells, monocytes and activated dendritic cells predicted a better prognostic value. Kaplan–Meier curves visualizes the above results (Fig. 4a–i). The consequences of the univariate analyses further sought out that these nine TIICs specifically influence patients’ prognosis and paved the way for further screening independent predictors.

Table 2 Univariate and multivariate analysis between 22 subpopulations of TIICs and survival in experiment group
Fig. 4

Nine subpopulations of TIICs significantly associated with the prognosis of patients with glioma in the experiment group. ae Kaplan–Meier plots visualized high level of immune cells, including TFH cells (p < 0.001), resting NK cells (p = 0.036), M0 macrophages (p < 0.001), M1 macrophages (p = 0.019) and resting dendritic cells (p = 0.004), associated with poor OS. fi Kaplan–Meier plots visualized high level of immune cells, including plasma cells (p = 0.028, activated NK cells (p = 0.003), monocytes (p < 0.001) and activated dendritic cells (p = 0.04), associated with good OS

Identify several independent prognostic factors

Multivariate analysis was conducted to identify independently prognostic factors and the p-values were also shown in Table 2. The result suggested that these three TIICs, including TFH cells, activated NK cells and M0 macrophages, may serve as independent predictors of the progression of glioma, especially in the MT of LGG to GBM. Based on the correlation coefficients obtained by multivariate analysis, we constructed an immune risk score (IRS) model based on these three immune cells: IRS = 14.339*TFH cells + 2.354* M0 macrophages − 22.69* activated NK cells. ROC curve was than drawn and we therefore got the AUC = 0.732 (Fig. 5a) and figured out cut-off value = − 0.43124. Next, we divided the experiment group into high- and low-risk groups by cut-off value after calculating IRS. K–M curve indicated a significant difference in survival between high- and low-risk groups (Fig. 5b). What’s more, IRS scores showed statistical differences with molecular subtypes of glioma (Table 3).

Fig. 5

IRS construction and validation. a ROC curve of the IRS model in experiment group. b Kaplan–Meier curve visualized the overall survival of experiment group based of the level of IRS. c Kaplan–Meier curve visualized the overall survival of validation group based of the level of IRS

Table 3 Relationships between IRS scores and molecular subtypes of glioma

Verify the prognostic model in the validation group

For the validation group, we calculated IRSs and also divided it into high- and low-risk groups based on the cut-off value. The K–M curve showed the survival rate of patients in the high-risk group was significantly worse than that in the low-risk group (Fig. 5c). This was a good proof of the validity of the IRS model we constructed.

IHC confirmed the CIBERSORT result

In order to verify the explorative data obtained for TFH cells, activated NK cells and M0 macrophages, we evaluated these cells density by immunohistochemistry in 5 human LGG tumor tissues and 5 human GBM tumor tissues Examples of these cells tryptase staining and quantification summary are shown in Fig. 6a–r. In agreement with CIBERSORT results, activated NK cells was reduced in GBM while TFH cells and M0 macrophages were increased in GBM.

Fig. 6

IHC of three significant TIICs. ac, gi Level of TFH cells, M0 macrophages and activated NK cells in LGG group. df, jl Level of TFH cells, M0 macrophages and activated NK cells in HGG group. mr The quantification of TFH cells (p < 0.001), M0 macrophages (p < 0.001) and activated NK cells (p < 0.001) are shown

Differential expressed genes and enrichment analysis based on prognostic model

We calculated IRSs for all samples and divided them into high- and low-risk group for TCGA and CGGA databases. After obtaining the differential genes of the two databases separately, we found that we got 118 common differential expressed genes for the intersection of the results (Fig. 7). For these genes, online tools “STRING” was used for GO/KEGG enrichment analysis and found that they enriched in the following biological processes (Table 4) which were mainly related to immune response.

Fig. 7

Venn diagram of the differential expressed genes. 249 DEGs from CGGA datasets and 166 DEGs from TCGA were taken to obtain the 39 common DEGs

Table 4 GO enrichment and KEGG pathway analysis results with grouped all samples into high- and low-IRS group


For a long time, although it is well known that immune cells play an important role in tumor initiation and development, these insights have few influence on clinical practice [11, 12]. In addition, the role of genes that are abnormally expressed in tumor tissues in diagnosis and prognosis has also attracted widespread attention; however, few studies have focused on the differential distribution of immune cells between different components. In this paper, we firstly established an immune risk score model based on the fractions of three subpopulations of TIICs. Compared with the high-IRS group based on our model, the low-IRS group has a significantly better survival rate (p < 0.001). This finding suggests that our IRS model can better predict progression of glioma, especially in the MT from LGG to GBM. Validation group, IHC and functional enrichment analyses further illustrate the validity of the model. This study opens a door for a better understanding of new diagnosis strategy from the perspective of TIICs. We acknowledge that there exist limitations in this research, particularly no precise analysis of the effect of single TIICs. Besides, studies on TIICs in initiation of glioma weren’t carried on due to lack of sequencing samples from normal people in these public databases. Therefore, further studies are urgently needed to analysis single TIISs and whether it is possible to detect the real-time progression of tumor through the state of immune cells in the circulatory system.

Gliomas are tumors of CNS, originating from transformed neural stem or progenitor glial cells [13]. On the basis of histopathological characteristics WHO classified gliomas into groups: low-grade gliomas (LGG, grades I and II) are well differentiated, slow-growing tumors, whereas high-grade gliomas (HGG, grades III and IV) are less differentiated or anaplastic, and strongly infiltrate brain parenchyma [14]. Glioblastoma (GBM) is categorized as the most malignant type (grade IV). It accounts for 50% of CNS tumors and is a deadly disease without curable therapy. Despite aggressive treatments, such as extensive resection combined with radiation and/or chemotherapy, patients with GBM eventually die of their disease [4]. In another aspect, patients with LGG may survive for many years, but after transformed to GBM, survival rates rapidly decline [15, 16]. A population-based study showed that the mean period of malignant transformation from LGG to GBM was 5.3 years and for anaplastic astrocytoma to GBM was 1.4 years [17]. Most of the predictive models established in previous studies on glioma development and malignant transformation were based on differential expressed genes, but they neglected that immune cells may also play an important role in tumorigenesis. Due to technical limitations, previous researches were limiting to a narrow insight of tumor-infiltrating cells. Immunohistochemistry and flow cytometry which depend on a single surface marker were used to evaluated TIICs. Apparently, these techniques may have misidentified other cell with the same surface markers as TIICs and are subjectively affected by observers. Thus, in the current study, we employed a silicon analysis, known as CIBERSORT, to infer the proportions of 22 immune cell subpopulations from glioma transcriptomes. CIBERSORT is a deconvolution algorithm for charactering TIICs composition of complex tissues by analyzing 547 gene expression, introduced by Newman etc. in 2015. They firstly employed a novel application of liner support vector regression to deconvolve the tissue composition. In order to assess the feasibility of TIICs deconvolution from bulk tumors, they then designed and validated a TIICs gene signature matrix, termed LM22. By using LM22 to deconvolve 3061 human transcriptomes, they therefore proved CIBERSORT has great specificity and sensitivity [11]. As an emerging technology, CIBERSORT have already conducted in breast cancer [18], lung cancer [19], colon cancer [6] and so on, all these studies demonstrated the effectiveness and accuracy of this tool when analyzing TIICs.

Univariate and multivariate analyses indicated TFH cells, activated NK cells and M0 macrophages as independent predictors. Then, based on their correlation coefficients, we firstly constructed such an IRS model in glioma. Among these correlation coefficients, or the degree to which the cell distribution correlated with tumor progression, the coefficient of activated NK cells is negative, while the coefficients of the other two TIICs are positive. This is consistent with our previous analyses between TIICs and tumor grade. Hence, we have adequate reason to believe that this model can predict MT between LGG and GBM well.

Immune system can be functionally divided into innate immunity and adaptive immunity, where adaptive immunity is antigen-specific. It mainly consists of B cell-mediated humoral immunity and cytotoxic T cell-mediated cellular immune responses, and both these two adaptive immunity processes require signals from CD4 T cells [9]. In one aspect, some CD4 T cell subpopulations such as Th1 cells can exert anti-tumor immunity by overcoming the tolerance of autoantibodies expressed by tumors, and these effectors T cells are advantageous for tumor immunotherapy [20]. However, other subsets of CD4 T cells, particularly regulatory T cells and TFH cells, inhibit tumor immunity, thereby promoting cancer growth [21,22,23]. In our study, although there was no significant difference in the composition of TFH cells between LGG group and HGG group, but in Fig. 2b we can see that its level in the GBM group is higher than that of the lower grade gliomas. At present, there is no research on the role of TFH in the immune microenvironment of glioma, which is the problem we need to think about and solve next.

Unlike T cells, NK cells play a unique role in innate and adaptive immune responses without the involvement of major histocompatibility complex (MHC) antigens or antibodies [24], and monitor status of intracellular bacteria, viruses-infected cells and transformed cells. Activated NK cells are one of two types of lymphokine-activated killer (LAK) cells. When stimulated by IL-2, they become activated against tumor cells. Although no randomized controlled trails of immunotherapy with HGG by LAK has been performed to date, one study showed patients treated with LAK cells had longer survivals than control groups [25]. Due to the difficulties in producing sufficient LAK cells, researches on activated NK cells for immunotherapy of glioma have been restricted. We pointed out a significant difference in the distribution of activated NK cells between low- and high-grade gliomas (p < 0.001) (Fig. 2c, d), and the lower the level of activated NK cells in the higher grade of gliomas (p < 0.001) (Fig. 3d). In studies of the association with glioma molecular subtype, the level of activated NK cells was the lowest in the mesenchymal subtype, which has the worst prognosis, while the other three subtypes harbor relatively higher level of it (p < 0.001) (Fig. 3g). In addition, as previously stated, the correlation coefficient of activated NK cells is also negative. These results indicate that activated NK cells may induce favorable clinical outcome of glioma, in another word, it may also be a vital suppressor for MT in LGG.

TAMs are macrophages infiltrating in tumor tissues which are the main composition in tumor microenvironment (TME). They differentiate through alternative pathways, among which the most common one is Notch pathway [26, 27]. What’s more, they facilitate tumor progression [28]. Once activated, monocytes continue to differentiate, first differentiated into M0 macrophages and then M1 and M2 arisen from M0. Others have shown that increased level of M0 is associated with poor clinical outcomes of lung adenocarcinoma [29]. So far, no clear experiments have been conducted to demonstrate the relationship between TAMs and glioma prognosis. Some people believe that TAMs in gliomas may be affected by tumor tissues and show immunosuppressive effects [19]. According to our work, the contents of M0 (p < 0.001) in GBM is higher than that of LGG. Wilcoxon test result also exhibits a gradual increase in the level of M0 from LGG to GBM. Moreover, coefficient of our IRS model also indicate that M0 comes under the influence of tumor development and promote malignant progression.

To again insight into the immune-related biological processes during glioma progression, we performed GO biological process (GOBP), GO cellular component (GOCC), GO molecular function (GOMF) and KEGG pathway analysis. Not only do the top results are immune-related, it is particularly worth mentioning that the first of each analysis are all immunologically relevant. This proves the validity of our consequences to some extent, on the other hand, it also finds some hub pathways in the MT of glioma, which indicates a path for future researches.


In summary, our study expounded the distinct composition of tumor-infiltrating immune cells in different grades and molecular subtypes of glioma. The complex intersection between TIICs and MT was quantified by our IRS model. Finally, we pointed out some relevant pathways related to progression and MT of glioma. These findings deepen the understanding of immune responses in CNS tumors and may enable to develop more effective immunotherapeutic strategies.



tumor-infiltrating immune cells


malignant transformation


low-grade glioma


high-grade glioma


Cell type Identification By Estimating Relative Subsets Of RNA Transcripts


area under curve


T follicular helper




central nervous system


World Health Organization


glioblastoma multiforme


dendritic cell


natural killer


The Cancer Genome Atlas


Chinese Glioma Genome Atlas


gene oncology


Kyoto encyclopedia of genes and genomes


immune risk score




major histocompatibility complex


lymphokine-activated killer


tumor microenvironment


  1. 1.

    Ostrom QT, Gittleman H, Truitt G, Boscia A, Kruchko C, Barnholtz-Sloan JS. CBTRUS Statistical Report: primary brain and other central nervous system tumors diagnosed in the United States in 2011–2015. Neuro-oncology. 2018;20(suppl4):1–86.

    Article  Google Scholar 

  2. 2.

    Louis DN, Perry A, Reifenberger G, von Deimling A, Figarella-Branger D, Cavenee WK, et al. The 2016 World Health Organization Classification of tumors of the central nervous system: a summary. Acta Neuropathol. 2016;131(6):803–20.

    Article  Google Scholar 

  3. 3.

    Jiang T, Mao Y, Ma W, Mao Q, You Y, Yang X, et al. CGCG clinical practice guidelines for the management of adult diffuse gliomas. Cancer Lett. 2016;375(2):263–73.

    CAS  Article  Google Scholar 

  4. 4.

    Wong ET, Lok E, Swanson KD. Clinical benefit in recurrent glioblastoma from adjuvant NovoTTF-100A and TCCC after temozolomide and bevacizumab failure: a preliminary observation. Cancer Med. 2015;4(3):383–91.

    CAS  Article  Google Scholar 

  5. 5.

    Louveau A, Smirnov I, Keyes TJ, Eccles JD, Rouhani SJ, Peske JD, et al. Structural and functional features of central nervous system lymphatic vessels. Nature. 2015;523(7560):337–41.

    CAS  Article  Google Scholar 

  6. 6.

    Xiong Y, Wang K, Zhou H, Peng L, You W, Fu Z. Profiles of immune infiltration in colorectal cancer and their clinical significant: a gene expression-based study. Cancer Med. 2018;7(9):4496–508.

    CAS  Article  Google Scholar 

  7. 7.

    Young JS, Dayani F, Morshed RA, Okada H, Aghi MK. Immunotherapy for high grade gliomas: a clinical update and practical considerations for neurosurgeons. World Neurosurg. 2019.

    Article  PubMed  Google Scholar 

  8. 8.

    Zhao J, Chen AX, Gartrell RD, Silverman AM, Aparicio L, Chu T, et al. Immune and genomic correlates of response to anti-PD-1 immunotherapy in glioblastoma. Nat Med. 2019;25(3):462–9.

    Article  Google Scholar 

  9. 9.

    Parney IF. Basic concepts in glioma immunology. Adv Exp Med Biol. 2012;746:42–52.

    CAS  Article  Google Scholar 

  10. 10.

    Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

    Article  Google Scholar 

  11. 11.

    Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.

    CAS  Article  Google Scholar 

  12. 12.

    Fridman WH, Pages F, Sautes-Fridman C, Galon J. The immune contexture in human tumours: impact on clinical outcome. Nat Rev Cancer. 2012;12(4):298–306.

    CAS  Article  Google Scholar 

  13. 13.

    Galon J, Mlecnik B, Bindea G, Angell HK, Berger A, Lagorce C, et al. Towards the introduction of the ‘Immunoscore’ in the classification of malignant tumours. J Pathol. 2014;232(2):199–209.

    CAS  Article  Google Scholar 

  14. 14.

    Ilkhanizadeh S, Lau J, Huang M, Foster DJ, Wong R, Frantz A, et al. Glial progenitors as targets for transformation in glioma. Adv Cancer Res. 2014;121:1–65.

    CAS  Article  Google Scholar 

  15. 15.

    Wilcox JA, Ramakrishna R, Magge R. Immunotherapy in glioblastoma. World Neurosurg. 2018;116:518–28.

    Article  Google Scholar 

  16. 16.

    Huse JT, Phillips HS, Brennan CW. Molecular subclassification of diffuse gliomas: seeing order in the chaos. Glia. 2011;59(8):1190–9.

    Article  Google Scholar 

  17. 17.

    Ohgaki H, Kleihues P. Population-based studies on incidence, survival rates, and genetic alterations in astrocytic and oligodendroglial gliomas. J Neuropathol Exp Neurol. 2005;64(6):479–89.

    CAS  Article  Google Scholar 

  18. 18.

    Ali HR, Chlon L, Pharoah PD, Markowetz F, Caldas C. Patterns of immune infiltration in breast cancer and their clinical implications: a gene-expression-based retrospective study. PLoS Med. 2016;13(12):e1002194.

    Article  Google Scholar 

  19. 19.

    Liu X, Wu S, Yang Y, Zhao M, Zhu G, Hou Z. The prognostic landscape of tumor-infiltrating immune cell and immunomodulators in lung cancer. Biomed Pharmacother. 2017;95:55–61.

    CAS  Article  Google Scholar 

  20. 20.

    Tran E, Turcotte S, Gros A, Robbins PF, Lu YC, Dudley ME, et al. Cancer immunotherapy based on mutation-specific CD4+ T cells in a patient with epithelial cancer. Science. 2014;344(6184):641–5.

    CAS  Article  Google Scholar 

  21. 21.

    Shirota H, Klinman DM, Ito SE, Ito H, Kubo M, Ishioka C. IL4 from T follicular helper cells downregulates antitumor immunity. Cancer Immunol Res. 2017;5(1):61–71.

    CAS  Article  Google Scholar 

  22. 22.

    Munn DH, Mellor AL. The tumor-draining lymph node as an immune-privileged site. Immunol Rev. 2006;213:146–58.

    Article  Google Scholar 

  23. 23.

    Willimsky G, Blankenstein T. Sporadic immunogenic tumours avoid destruction by inducing T-cell tolerance. Nature. 2005;437(7055):141–6.

    CAS  Article  Google Scholar 

  24. 24.

    Kim HJ, Cantor H. CD4 T-cell subsets and tumor immunity: the helpful and the not-so-helpful. Cancer Immunol Res. 2014;2(2):91–8.

    CAS  Article  Google Scholar 

  25. 25.

    Moretta A, Bottino C, Vitale M, Pende D, Biassoni R, Mingari MC, et al. Receptors for HLA class-I molecules in human natural killer cells. Annu Rev Immunol. 1996;14:619–48.

    CAS  Article  Google Scholar 

  26. 26.

    Everson RG, Antonios JP, Liau LM. Cell-based immunotherapy of gliomas. Prog Neurol Surg. 2018;32:90–100.

    Article  Google Scholar 

  27. 27.

    Kim J, Bae JS. Tumor-associated macrophages and neutrophils in tumor microenvironment. Mediators Inflamm. 2016;2016:6058147.

    PubMed  PubMed Central  Google Scholar 

  28. 28.

    Sica A, Larghi P, Mancino A, Rubino L, Porta C, Totaro MG, et al. Macrophage polarization in tumour progression. Semin Cancer Biol. 2008;18(5):349–55.

    CAS  Article  Google Scholar 

  29. 29.

    Kurahara H, Takao S, Maemura K, Mataki Y, Kuwahata T, Maeda K, et al. M2-polarized tumor-associated macrophage infiltration of regional lymph nodes is associated with nodal lymphangiogenesis and occult nodal involvement in pN0 pancreatic cancer. Pancreas. 2013;42(1):155–9.

    CAS  Article  Google Scholar 

Download references


Not applicable.

Web links and URLs

The Cancer Genome Atlas (TCGA) dataset ( Chinese Glioma Genome Atlas (CGGA) dataset ( CIBERSORT web portal ( R software version 3.5.3 ( “Corrplot” ( “Pheatmap” ( “Vioplot” (


This work was supported by the National Natural Science Foundation of China (81772681 and 81670153), the Natural Science Foundation of Jiangsu Province (BK20160098), the Program for Development of Innovative Research Team in the First Affiliated Hospital of NJMU, and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Author information




JL, HL, ZC and HW: Conception and design. JL, HL, ZC and HW: Development of the methodology. JL, HL, ZC and LF: Acquisition of data. JL, SF and XC: Analysis and interpretation of data (e.g., statistical, biostatistics, computational analysis). JL and HW: Writing, review, and/or revision of the manuscript. JL, HL, ZC and HW: Study supervision. All authors read and approved the final manscript.

Corresponding author

Correspondence to Huibo Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Lu, J., Li, H., Chen, Z. et al. Identification of 3 subpopulations of tumor-infiltrating immune cells for malignant transformation of low-grade glioma. Cancer Cell Int 19, 265 (2019).

Download citation


  • Glioma
  • Tumor immune
  • Malignant transformation
  • Prognosis