key: cord-0062483-3brh04gq
authors: Liu, Liling; Yang, Yuanzhong; Duan, Hao; He, Jiahua; Sun, Lu; Hu, Wanming; Zeng, Jing
title: CHI3L2 Is a Novel Prognostic Biomarker and Correlated With Immune Infiltrates in Gliomas
date: 2021-04-15
journal: Front Oncol
DOI: 10.3389/fonc.2021.611038
sha: 7d26bde4ed7b59bbb73a1c9b88aeb1f9fc557311
doc_id: 62483
cord_uid: 3brh04gq

CHI3L2 (Chitinase-3-Like Protein 2) is a member of chitinase-like proteins (CLPs), which belong to the glycoside hydrolase 18 family. Its homologous gene, CHI3L1, has been extensively studied in various tumors and has been shown to be related to immune infiltration in breast cancer and glioblastoma. High CHI3L2 expression was reported to be associated with poor prognosis in breast cancer and renal cell carcinoma. However, the prognostic significance of CHI3L2 in glioma and its correlation between immune infiltration remains unclear. In this study, we examined 288 glioma samples by immunohistochemistry to find that CHI3L2 is expressed in tumor cells and macrophages in glioma tissues and highly expressed in glioblastoma and IDH wild-type gliomas. Relationships between CHI3L2 expression and clinical features (grade, age, Ki67 index, P53, PHH3 (mitotic figures), ATRX, TERTp, MGMTp, IDH, and 1p/19q co-deleted status) were evaluated. Kaplan-Meier survival was conducted to show high CHI3L2 expression in tumor cells (TC) and macrophage cells (MC) indicated poor prognosis in diffusely infiltrating glioma (DIG), lower-grade glioma (LGG), and IDH wild-type gliomas (IDH-wt). The overall survival time was higher in patients with dual-low CHI3L2 expression in TC and MC compared to those in patients with non-dual CHI3L2 expression and dual high expression in DIG and IDH wild-type gliomas. By univariate and multivariate analysis, we found that high CHI3L2 expression in tumor cells was an independent unfavorable prognostic factor in glioma patients. Moreover, we used two datasets (TCGA and CGGA) to verify the results of our study and explore the potential functional role of CHI3L2 by GO and KEGG analyses in gliomas. TIMER platform analysis indicated CHI3L2 expression was closely related to diverse marker genes of tumor immune infiltrating cells, including monocytes, TAMs, M1 macrophages, M2 macrophages, TGFβ1+ Treg and T cell exhaustion in GBM and LGG. Western Blot validated CHI3L2 is expressed in glioma cells and microglia cells. The results of flow cytometry showed that CHI3L2 induces the apoptosis of CD8+ T cells. In conclusion, these results demonstrate CHI3L2 is related to poor prognosis and immune infiltrates in gliomas, suggesting it may serve as a promising prognostic biomarker and represent a new target for glioma patients.

Gliomas comprise the bulk of primary brain tumors in adults (1) . Diffuse glioma is histopathologically classified into grade II-IV according to morphological criteria, including mitotic count, nuclear atypia, microvascular proliferation, and necrosis. Glioblastoma multiform (GBM) is categorized as one of the most malignant subtypes (2, 3) . The 2016 World Health Organization (WHO) classification of adult diffuse glioma combines tumor histological morphology and molecular features, including the isocitrate dehydrogenase (IDH) mutation and the chromosomal arms 1p and 19q complete deletion (1p/19q co-deletion) (4) . Even combining maximal surgical resection and radiotherapy with adjuvant temozolomide, tumor recurrence is inevitable and the prognosis of gliomas remains very poor (5) . Consequently, it is an urgent demand to discover the potential molecular characteristics of gliomas and look for more effective treatment strategies.

CHI3L2 (Chitinase-3-Like Protein 2), also known as YKL39, is a kind of secretory protein. It is a member of chitinase-like proteins (CLPs) which include CHI3L1, CHI3L2, SI-CLP, YM1 and YM2. CHI3L2 was originally isolated from the culture medium of primary human articular cartilage cells (6) . It has two physiological activities, one is to induce autoimmune response (7) , the other is to participate in tissue remodeling, both of which may lead to disease progression. Previous studies showed CHI3L2 mRNA is significantly up-regulated in osteoarthritis, Alzheimer's disease, multiple sclerosis, and amyotrophic lateral sclerosis patients (8) (9) (10) (11) . CHI3L2 was secreted by microglia/astrocytes and could increase monocyte/ macrophage infiltration, angiogenesis, and neuronal death in amyotrophic lateral sclerosis (11) . It is not yet clear what type of cells CHI3L2 is expressed in gliomas. However, previous studies on CHI3L2 have shown that macrophages are a possible source of CHI3L2 in tumors (12) (13) (14) (15) . CHI3L2 has a high degree of sequence identity with CHI3L1, but no cross-reactivity has been observed (16) (17) (18) . There have been many studies on the correlations between CHI3L1 and the progression of a number of cancers (19) (20) (21) (22) . In recent years, the relationship between tumor immune microenvironment and immunotherapy has received more and more attention. It was also reported CHI3L1 was related to immune infiltration in breast cancer and glioblastoma (23, 24) . However, the data about the role of CHI3L2 in cancers and its association with immune infiltrates are fragmentary. Previous studies reported that CHI3L2 was overexpressed in tumor-associated macrophages and related to poor outcomes in breast cancer and renal cell carcinoma (13) (14) (15) . Studies have been shown that CHI3L2 mRNA expression was increased in gliomas (18, 25, 26) . However, the prognostic significance of CHI3L2 and its correlation with immune infiltrates in glioma remain unclear.

To systematically explore the CHI3L2 protein expression in diffusely infiltrating glioma, we first evaluated the CHI3L2 expression levels of 288 glioma tissues by immunohistochemistry (IHC) and analyzed the association between CHI3L2 levels and clinicopathological parameters. Moreover, we took advantage of CHI3L2 transcriptional data of gliomas in The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA) datasets to validate our findings. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were used to explore the potential biological process and pathways of CHI3L2 in glioma. The Tumor IMmune Estimation Resource (TIMER) platform was used to explore the correlations between CHI3L2 and diverse marker genes of tumor immune infiltrates. Finally, we further verified the results through Western Blot and flow cytometry. In gliomas, this is the first comprehensive study to elaborate on the clinical significance of CHI3L2, its influence on prognosis and its correlation with immune infiltrates.

We enrolled 288 glioma patients (WHO grade II-IV) operated at the Sun Yat-sen University Cancer Center (Guangzhou, China) from January 2009 to January 2016. The median follow-up time was 54 months. Follow-up was last done in June 2019. The detailed clinical data are listed in Table S1 . There were 167 males and 121 females. The median age of all patients at initial diagnosis was 43 years (range 7-78years). According to MRI imaging, 278 cases of glioma were located on the supratentorial and 10 cases of glioma were located infratentorial. 264 out of 288 patients received postoperative adjuvant treatment (radiotherapy or chemotherapy). The median overall survival time of all patients was 27 months (range 0-110 months). This cohort included 112 cases of astrocytoma, 45 cases of oligodendroglial gliomas and 131 cases of glioblastoma (WHO IV). All samples were ethically approved for use based on informed consent.

The human glioma cell lines and a human microglia cell line were purchased from the American Type Culture Collection resource center. The human glioma cells were maintained in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% heat-inactivated fetal bovine serum (FBS) and the human microglia cells were cultured in Minimum Essential Medium (MEM) with 10% FBS at 37°C in a humidified incubator containing 5% CO2. Peripheral blood mononuclear cells (PBMC) were isolated by Ficoll-Hypaque density gradient centrifugation (Solarbio, Beijing, China). CD8+ T cell were separated by positive selection from PBMCs with CD8 magnetic beads and cultured in RPMI-1640 medium supplemented with 10% human serum, 5% L-glutamine-penicillin-streptomycin solution (Sigma-Aldrich, USA), CD3/ CD28 antibody (Biolegend, USA) (25ul/ml) and IL-2 (100IU/ ml) in 24-well plates. After culturing for 24 hours, add the corresponding concentration of human CHI3L2 protein (Sino Biological, Beijing, China) to T cells and culture for 72 hours at 37°C in a humidified incubator containing 5% CO2.

Immunohistochemistry was essentially performed as previously reported (27) . These tissue specimens were incubated with CHI3L2 rabbit polyclonal antibody (#22164, SAB, Maryland, USA). Immunohistochemical evaluation was independently conducted by two pathologists blinded for patient characteristics and outcome, and CHI3L2 expression by tumor cells and macrophage cells was scored separately. The discrepancies were resolved by consensus under a microscope for multi-viewing. A semi-quantitative IHC scoring criterion was used to determine the CHI3L2 protein expression levels in tumor cells. The percentage of positive cells and staining intensity were assessed to improve accuracy. The percent positivity of staining cells range from 0 to 4: 0, none; 1, 1%-25%; 2, 26%-50%; 3, 51%-75%; 4, 76%-100%. The intensity of staining was graded from 0 to 3 (0, no staining; 1, weak; 2, moderate and 3, strong). Then, we obtained the final IHC score by multiplying the proportion score by the intensity score of staining. We chose 4.5, which was determined by the Youden index as an optimal cutoff point to separate low CHI3L2 expression (score of 0-4.5) from high CHI3L2 expression (score>4.5) in tumor cells. For the macrophage cells, we only count the number of CHI3L2 positive staining macrophages, regardless of the staining intensity. We designed 7.5 determined by Youden index as an optimal cutoff point to differentiate low expression (number≤ 7.5) from high expression (number>7.5) in macrophages.

The other antibody markers including PHH3, P53, Ki67, ATRX, CD163, CD4, CD8, and CD20 were also used by immunohistochemistry tests. We detected MGMT promoter methylated status, TERT promoters and IDH mutation status by Sanger sequencing. 1p and 19q deletion status was detected using fluorescent in situ hybridization (FISH). The detailed protocol and assessment standard were described as a previous study (28) .

The CHI3L2 RNA-seq data were downloaded from http://www. cgga.org.cn/. We totally analyzed 601 TCGA RNA-seq cohort and 608 CGGA RNA-seq cohort of gliomas, ranging from WHO grade II to grade IV.

To identify the CHI3L2-related genes, the limma package of R software was used to screen out the differentially expressed genes (DEGs). The top 26 hub genes of the overlapping DEGs were built via the plug-in molecular complex detection and cytoHubba of Cytoscape. To explore the functions and pathways of CHI3L2-related genes, we performed GO and KEGG analyses on ClueGo and Metascape websites. The TIMER platform (https://cistrome.shinyapps.io/timer/) (29, 30) was performed to explore the association between CHI3L2 and marker sets of tumor immune infiltrates in GBM and LGG (31) .

Total protein was extracted from seven human glioma cell lines and one human microglia cell line (HMC3). 30 ug of protein was loaded onto 10% SDSPAGE and electrophoretically transferred to PVDF membranes. After blocking, the membranes were incubated with primary antibody against CHI3L2 (1:1000 dilutions, rabbit polyclonal anti-CHI3L2, #22164, SAB, Maryland, USA). The membranes were then incubated with horseradish peroxidase-linked anti-rabbit antibody (at a 1:3000 dilution, Santa Cruz Biotechnology, Inc., Santa Cruz, Calif., USA). B-Actin was served as a loading control. 

GraphPad Prism 8 and SPSS 22 software were performed for statistical analyses. The measurement data are represented as mean ± SD. The Chi-square test was conducted to explore the correlations between CHI3L2 levels and clinicopathological features. Kaplan-Meier analysis was conducted for the overall survival of glioma patients with the log-rank test. The Cox proportional hazards regression model was used for univariate and multivariate analyses to evaluate the independence of CHI3L2 in predicting prognosis. The association between CHI3L2 and marker genes of immune infiltrating cells was assessed by Spearman's correlation coefficients. P < 0.05 was regarded as statistically significant.

We detected CHI3L2 protein expression levels in histological sections from patients with different glioma grades by immunohistochemistry (IHC). Among the 288 glioma specimens inspected, we found CHI3L2 was mainly stained in tumor cells, as well as macrophage cells. cells were upregulated with increasing WHO grade of gliomas ( Figure 2A ), but there was no significant difference in CHI3L2+ macrophage density between WHO II and WHO III gliomas ( Figure 2D ). The CHI3L2 expression levels of GBM were significantly increased compared with LGG (WHO II-III) (P<0.001) in tumor cells ( Figure 2B ) and macrophages ( Figure 2E) . A significant increase of CHI3L2 expression levels was found in IDH-wildtype gliomas compared with IDH-mutant gliomas (P<0.001) in tumor cells ( Figure 2C ) and macrophages ( Figure 2F ). We further analyze the CHI3L2 IHC score and density of CHI3L2+ macrophages in diffusely infiltrating glioma of new molecular classification, including IDH mutant without 1p/19q codeleted gliomas, IDH mutant with 1p/19q codeleted gliomas, and IDH wild-type gliomas ( Figures S1A, B) . We found the expression of CHI3L2 is not related to the status of 1p/19q codeleted in IDH mutant gliomas. Based on the expression levels of CHI3L2 in tumor cells and macrophage cells, we evaluated the association between CHI3L2 staining and clinicopathological factors, as listed in Table 1 . In tumor cells, we found significant correlations between CHI3L2 expression and WHO grade (P<0.001), age (P=0.001), Ki67 (P<0.001), P53 (P=0.034), PHH3 (mitotic figures) (P<0.001), ATRX protein expression (P=0.026), IDH (P<0.001) and 1p/19q codeleted (P=0.002). In macrophage cells, CHI3L2 is strongly correlated with WHO grade (P<0.001), gender (P=0.008), age (P=0.006), Ki67 (P<0.001), PHH3 (mitotic figures) (P<0.001), and IDH status (P<0.001). However, CHI3L2 expression levels of glioma cells were not significantly related to gender, location, TERT promoter mutation status, and MGMT promoter methylated status. In macrophage cells, CHI3L2 expression has no correlation with location, P53, ATRX protein expression, MGMT promoter methylated status, TERT promoter mutation, and 1p/19q codeleted status.

To explore the prognostic significance of CHI3L2 in gliomas, we performed the Kaplan-Meier method and log-rank test. We found high CHI3L2 expression levels of tumor cells and macrophages significantly predicted worse overall survival in diffusely infiltrating glioma (DIG) (Figures 3A, D) and lowergrade glioma (LGG) patients ( Figures 3B, E) . However, there was no statistical significance difference in GBM in our cohort ( Figures 3C, F) . When considering the CHI3L2 expression of tumor cells and macrophages together, we found the overall survival time was higher in patients with dual-low CHI3L2 expression in TC and MC compared to those in patients with non-dual CHI3L2 expression and dual high expression in DIG ( Figure 3G ), but this difference is not statistically significant in LGG and GBM ( Figures 3H, I) . Similarly, we also analyze the effect of CHI3L2 on the prognosis in the new molecular classification of glioma. We found CHI3L2 expression in tumor cells is closely related to the prognosis of all new molecular classification of glioma, and high CHI3L2 expression in tumor cells, macrophages and TC + MC predicted poor outcome for IDH wild-type gliomas ( Figure S2 ). Furthermore, we found regardless of whether patients with glioma have methylation of the MGMT promoter or have received adjuvant therapy, high CHI3L2 indicates a poor prognosis for glioma ( Figure S3) . Additionally, to evaluate the independent risk factors for prognosis of glioma, we conducted the univariate analysis ( Table 2 ) and multivariate analysis ( Table 3) . In univariate analysis, CHI3L2 expression in tumor cells, CHI3L2+ macrophage cells density, CHI3L2 expression in both tumor cells and macrophage cells (TC + MC), grade, age, location, adjuvant therapy, Ki67 index, PHH3 (mitotic figures), IDH, and 1p/19q codeleted status were shown to be prognostic variables for the prognosis of overall survival in glioma patients ( Table 2) . Then we included the prognostic variables in the univariate analysis (P<0.05) into the multivariate analysis. We found CHI3L2 expression in tumor cells, location, Ki67, IDH, 1p/19q codeleted were independent prognostic factors in gliomas ( Table 3) .

To further verify the results of our study, we collected a total of 601 glioma samples from the TCGA dataset and 608 glioma samples from the CGGA dataset to analyze the CHI3L2 mRNA expression. In the TCGA dataset, CHI3L2 mRNA levels were significantly increased in GBM (WHO IV) compared with WHO II, WHO III, and LGG patients (Figures 4A, B) . CHI3L2 mRNA expression levels were significantly higher in IDH wild-type gliomas compared with IDH mutant gliomas ( Figure 4C) . Similar results were also obtained in the CGGA dataset ( Figures 4D-F) . We further analyze the CHI3L2 mRNA expression levels in new molecular classification of diffusely infiltrating glioma, including IDH mutant without 1p/19q codeleted gliomas, IDH mutant with 1p/19q codeleted gliomas, and IDH wild-type gliomas, in TCGA and CGGA database ( Figure S4A and Figure S4B ). We found CHI3L2 mRNA expression levels in IDH wild-type gliomas are higher than IDH mutant gliomas. Gliomas with IDH mutant and 1p/19q codeleted have higher CHI3L2 mRNA levels than gliomas with IDH mutant and non-1p/19q codeleted. Moreover, we performed Kaplan-Meier analysis to confirm whether CHI3L2 mRNA levels could predict poor prognosis of gliomas in datasets. As shown in Figure 5 , patients with high CHI3L2 mRNA levels correspond to shorter survival time in all glioma subgroups both in the TCGA (Figures 5A-C) and CGGA ( Figures 5D-F) datasets. Similarly, we also verified the effect of CHI3L2 on the prognosis in the new molecular classification of glioma in database ( Figure S5) . Except for gliomas with IDH mutant and non-1p/19q codeleted in the TCGA dataset, high levels of CHI3L2 mRNA in any other subgroup indicate a poor prognosis, whether in the TCGA or CGGA dataset.

The GBM and LGG RNA-seq data were from TCGA and CGGA datasets. Limma package in R software was conducted to screen out the differentially expressed genes (DEGs) with the cut-off criterion of adjusted P< 0.05 and |log2FC| > 1. We identified 1356 overlapping DEGs which were aberrantly expressed in TCGA and CGGA datasets ( Figure 6A ). The top 26 hub genes were screened via the plug-in molecular complex detection and cytoHubba of Cytoscape ( Figure 6B ). GO analysis was performed to show the overlapping DEGs were involved in several biological processes, including angiogenesis, immune, and inflammatory response ( Figure 6C ). The KEGG pathways enriched in several classic signaling pathways, such as cell adhesion molecules (CAMs) and PI3K-Akt signaling pathways ( Figure 6D ).

Infiltrating immune cells are important components of the tumor microenvironment and are frequently associated with tumor behavior and patient outcomes. Since GO analysis revealed that CHI3L2 was related to the immune response, we further explored the infiltration of immune cells in gliomas. To estimate the relevance of CHI3L2 and diverse immune cell markers, we used the TIMER platform to investigate correlations between CHI3L2 levels and markers of diverse immune cells, included monocytes, TAMs (tumor-associated macrophages), M1 and M2 macrophages,

LGG ( Table 4) . We found CHI3L2 was significantly associated with marker sets of monocytes, TAMs, and M2 macrophages in GBM and LGG. Particularly, we showed the scatter plots of association between CHI3L2 and the marker sets of monocytes, TAMs, M1 phenotype, and M2 phenotype in GBM and LGG ( Figures 7A-H) . We also found significant correlations between CHI3L2 and some markers of Treg and T cell exhaustion, such as TGFb1, CTLA4, TIM-3, and GZMB. Since Treg and T cell exhaustion play an important role in tumor immune escape. We believe CHI3L2 may also play an immunomodulation role in gliomas. In addition, we used several clinical commonly immune cell markers, including CD163, CD4, CD8, CD20, to perform immunohistochemical test on glioma samples, and found that CHI3L2+ macrophages have a certain correlation with CD163+ M2 macrophages (r=0.547, p<0.001), CD4+ T cells (r=0.330, p<0.001), CD8+ T cells (r=0.389, p<0.001), and CD20+ B cells (r=0.237, p<0.001) in gliomas ( Figure S6 ).

The expression of CHI3L2 in glioblastoma cells (U251, U87, T98G, DBTRG, A172, LN229) and microglia cell (HMC3) has been verified by Western Blot ( Figure 8A ). Figure 8A shows that CHI3L2 is expressed in glioblastoma cell lines and a microglia cell line. It is strongly expressed in the glioblastoma cell lines U87, U251, LN229, A172 and the microglia cell line HMC3, while the expression in the glioblastoma cell lines T98G and DBTRG is weak. Figure 8B is the result of flow cytometry analysis. The left image is a representative sorting that lists the percentage of cells in each quadrant: bottom left-live cells; top left-mechanically damaged cells; bottom right-early apoptosis; top right-late apoptosis. The cellular apoptotic rate was a sum of early and late apoptotic rates. The proportion of apoptotic cells in the control group was 28.1%, the proportion of apoptotic cells in the 0.5ug/ml CHI3L2 group was 33.6%, and the proportion of apoptotic cells in the 2.5ug/ml CHI3L2 

At present, the outcome of most glioma is very poor, even with the use of comprehensive treatment strategies. It has been widely reported that the therapeutic resistance of glioma is closely related to its unique metabolic mechanism and the surrounding complex immunosuppressive microenvironment (32) (33) (34) . Therefore, exploring reliable prognostic biomarkers and personalized treatment strategies for this disease are urgently needed. In the present study, CHI3L2 has been identified as a novel prognostic biomarker and associated with tumor immune infiltration markers in gliomas, which indicate CHI3L2 may serve as a target for glioma treatment in the future. CHI3L2, as a member of the glycoside hydrolases 18 family, can act as a cytokine and growth factor but lacks chitinase activity (17) . It was found produced by tumor-associated macrophages in breast cancer (13, 14) . In renal cell carcinoma, MAP kinase signaling cascade in 293 and U87 MG cells (18, 26) . However, the correlations between CHI3L2 expression and clinicopathological features, the association with tumorinfiltrating immune cells, the prognostic value of CHI3L2, and its other functions in gliomas are still unknown.

Our study showed CHI3L2 expressed in tumor cells and macrophage cells in glioma tissues and particularly up-regulated in GBM and IDH wild-type gliomas. The Kaplan-Meier curves reveal higher CHI3L2 expression levels correlated with short overall survival in diffusely infiltrating glioma, lower-grade glioma, and IDH wild-type gliomas. High CHI3L2 expression indicates a poor prognosis for glioma patients, regardless of whether the MGMT promoter is methylated or has received adjuvant therapy. Cox proportional hazards regression model indicates CHI3L2 expression in tumor cells is an independent prognostic indicator of glioma. TCGA and CGGA datasets further confirm our findings. However, it should be pointed out that no significant differences between high CHI3L2 and poor prognosis in GBM were achieved in our cohort, which is different from the results (The relationship between high CHI3L2 mRNA expression and short overall survival time were statistically significant in all subgroups) in TCGA and CGGA datasets. We believe there are two reasons for the inconsistent results. One probable reason is the difference in detection levels. The relative expression levels of CHI3L2 mRNA were detected using high-throughput sequencing in TCGA and CGGA datasets, but CHI3L2 expression levels in our samples were assessed at protein levels by immunohistochemistry. Another possible reason is the difference in sample size. Accordingly, we intend to enlarge our sample size in the following study.

Based on these results, we further performed GO and KEGG pathway analyses to conclude the CHI3L2-related genes were involved in several biological processes, including angiogenesis, immune, and inflammatory response, and enriched in several classic signaling pathways, including cell adhesion molecules and PI3K-Akt signaling pathway. It has been reported that CHI3L2 acts as a powerful monocyte chemotactic factor and angiogenesis stimulating factor in breast cancer (14) . A recent review also reported CHI3L2 acts as a new target for anti-angiogenic therapy in breast cancer patients (35) . The angiogenesis function of CHI3L2 may be responsible for the poor prognosis of glioma, which needs further confirmation by follow-up studies. A recent review reported the role of cell adhesion molecules (CAM) in immune responses and tumor microenvironment--Cell adhesion molecules affect the antigen-presenting function, and inhibit the development of regulated cells and the leaching of regulatory cells into tumors, thus promoting tumor immune escape (36) . Our study showed CHI3L2 expressed in tumor cells and macrophages in glioma tissues. A previous study suggests human glioma-infiltrating macrophages have similar functions to CAM in mediated immune responses (37) . It was also reported that CAMs are potential prognostic biomarkers and attractive therapeutic targets for glioblastoma (38) . A previous study suggested PI3K-Akt signaling pathway activation, to some extent, affects the activity of most immune cell types. PI3K-AKT-mTOR pathway plays a certain role in regulating immunosuppression in tumor microenvironment (39) . We speculate CHI3L2 may be able to act as immunomodulation through this pathway. Additionally, we have learned in previous studies that CHI3L1 (the homologous gene of CHI3L2) may be used as an immunomodulatory factor to affect the therapeutic efficacy of PI3K/AKT-based pathway inhibitors in glioblastoma (40) . CHI3L1 also plays a key role in inducing immunosuppression and metastasis in breast cancer. CHI3L1 up-regulates pro-inflammatory mediators, CCL2, CXCL2 and MMP-9, all of which contribute to tumor growth and metastasis, and treatment with chitin can significantly reduce these effects (23) . These studies provide a reference for us to further explore the internal mechanism of CHI3L2 in immune infiltration.

Based on the analysis of the TIMER platform, the correlations between CHI3L2 and markers of immune cells imply CHI3L2 may play a part in immunomodulation in GBM and LGG. Our results suggest CHI3L2 expression has strong correlations with marker sets, include CD86 and CSF1R of monocytes, CCL2, CD68, and IL 10 of TAMs and CD163, VSIG4, and MS4A4A of M2 macrophages in GBM and LGG. A study has shown that purified CHI3L2 strongly induces the migration of freshly isolated human CD14+ monocytes (14) . It was reported that CD163 could act as a regulator of immune response and the potential to be a target to suppress immune escape and recover the function of T-cell populations in gliomas (41) . In our experiments, we also found a strong correlation between CHI3L2 and CD163 ( Figure S6) . Additionally, there was a Treg  and T cell exhaustion, including TGFb1, CTLA4, TIM-3, and GZMB. In the tumor microenvironment, TGFb can serve as an anti-tumor immunosuppressive factor and play an essential part in Treg cells (42) . It was reported that only TGFb, the key regulatory factor of tumor progression, was able to stimulate CHI3L2 mRNA levels in human macrophages in vitro (43) . TGF-b, which can be secreted by both microglial cells and glioma cells, participates in the functional transformation of macrophages into immunosuppressive and pro-invasive phenotypes, which supports tumor growth (44, 45) . Macrophages are believed to be activated microglia within the central nervous system. Our data show that CHI3L2 is expressed in tumor cells and macrophages in glioma tissues, and has a certain correlation with TGF-b, further suggesting that CHI3L2 may also have the function of inhibiting tumor immune regulation and promoting tumor growth. Similarly, CTLA4 and TIM-3 can induce T cell exhaustion via direct mechanisms through the interactions with their ligands, leading to impaired T cell activation, inhibition of T cell proliferation, and impaired cytokine release (46) . The correlation between CHI3L2 and T cell exhaustion markers indicates CHI3L2 may play a part in mediating T cell depletion. Our flow cytometry results further confirmed that CHI3L2 can induce CD8+ T cell apoptosis, indicating that CHI3L2 has the potential to promote tumor immune escape. However, how CHI3L2 protein promotes the apoptosis of CD8+ T cells needs to be further explored in future research. The limitations of this study are as follows: First, the sample size of gliomas for IHC is limited. Second, we did not accurately define the type of CHI3L2 + macrophages. In addition, further experimental investigation and analysis are needed to gain insights into the underlying mechanisms.

In conclusion, our study suggests CHI3L2 may be a promising prognostic biomarker that contributes to poor prognosis for gliomas. CHI3L2 may also play an important role in immunomodulation, suggesting CHI3L2 may serve as a novel therapeutic target for glioma patients.

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

The studies involving human participants were reviewed and approved by Sun Yat-Sen University Cancer Center. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

LL and WH designed the study and wrote the manuscript. LL, WH, JH, HD, and LS acquired the data. YY and HD provided help for the IHC test and Flow cytometry analysis. LL, YY, and WH performed the data analysis. WH and JZ reviewed the manuscript. All authors contributed to the article and approved the submitted version.

CBTRUS Statistical Report: Primary Brain and Other Central Nervous System Tumors Diagnosed in the United States in 2011-2015

WHO Classification of Tumours of the Central Nervous System

Genetic and Molecular Epidemiology of Adult Diffuse Glioma

The 2016 World Health Organization Classification of Tumors of the Central Nervous System: A Summary

Exciting New Advances in Neuro-Oncology The Avenue to a Cure for Malignant Glioma

Enhanced Expression of the Human Chitinase 3-like 2 Gene (YKL-39) but not Chitinase 3-like 1 Gene (YKL-40) in Osteoarthritic Cartilage

The Prevalence of Autoantibodies Against Cartilage Intermediate Layer Protein, YKL-39, Osteopontin, and Cyclic Citrullinated Peptide in Patients with Early-Stage Knee Osteoarthritis: Evidence of a Variety of Autoimmune Processes

Human YKL39 (Chitinase 3-like Protein 2), an Osteoarthritis-Associated Gene, Enhances Proliferation and Type II Collagen Expression in ATDC5 Cells

Chitinase Expression in Alzheimer's Disease and Non-Demented Brains Regions

Chitinase 3-like Proteins as Diagnostic and Prognostic Biomarkers of Multiple Sclerosis

CHI3L1 and CHI3L2 Overexpression in Motor Cortex and Spinal Cord of sALS Patients

Expression Profiles for Macrophage Alternative Activation Genes in AD and in Mouse Models of AD

Expression of M2 Macrophage Markers YKL-39 and CCL18 in Breast Cancer is Associated With the Effect of Neoadjuvant Chemotherapy

Tumor-Associated Macrophages in Human Breast Cancer Produce New Monocyte Attracting and Pro-Angiogenic Factor YKL-39 Indicative for Increased Metastasis After Neoadjuvant Chemotherapy

M2 Macrophage Marker Chitinase 3-Like 2 (CHI3L2) Associates With Progression of Conventional Renal Cell Carcinoma

Isolation and Sequence of a Novel Human Chondrocyte Protein Related to Mammalian Members of the Chitinase Protein Family

Role of Chitinase-Like Proteins in Cancer

Overexpression of YKL-39 Gene in Glial Brain Tumors

High Serum YKL-40 Level in Patients With Small Cell Lung Cancer is Related to Early Death. Lung Cancer

YKL-40 and Matrix Metalloproteinase-9 as Potential Serum Biomarkers for Patients with High-Grade Gliomas

Serum YKL-40, A New Prognostic Biomarker in Cancer Patients?

Overexpression of YKL-40 is an Independent Prognostic Marker in Gastric Cancer

Induction of Proinflammatory Mediators by CHI3L1 is Reduced by Chitin Treatment: Decreased Tumor Metastasis in a Breast Cancer Model

Prognostic Value and Immune Cell Infiltration of Hypoxic Phenotype-Related Gene Signatures in Glioblastoma Microenvironment

Experimental Anti-Angiogenesis Causes Upregulation of Genes Associated With Poor Survival in Glioblastoma

Chitinase 3-like Protein 2 (CHI3L2, YKL-39) Activates Phosphorylation of Extracellular Signal-Regulated Kinases ERK1/ ERK2 in Human Embryonic Kidney (HEK293) and Human Glioblastoma (U87 MG) Cells

Expression of CPEB4 in Human Glioma and Its Correlations With Prognosis

Practice of the New Integrated Molecular Diagnostics in Gliomas: Experiences and New Findings in a Single Chinese Center

TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells

0 for Analysis of Tumor-Infiltrating Immune Cells

LAYN Is a Prognostic Biomarker and Correlated With Immune Infiltrates in Gastric and Colon Cancers

Inhibition of Colony Stimulating Factor-1 Receptor Abrogates Microenvironment-Mediated Therapeutic Resistance in Gliomas

Cancer Stem Cells and Immunosuppressive Microenvironment in Glioma

YKL-39 as a Potential New Target for Anti-Angiogenic Therapy in Cancer

Cell Adhesion Molecules and Their Roles and Regulation in the Immune and Tumor Microenvironment

The Role of Human Glioma-Infiltrating Microglia/Macrophages in Mediating Antitumor Immune Responses

Adhering Towards Tumorigenicity: Altered Adhesion Mechanisms in Glioblastoma Cancer Stem Cells

PI3K-AKT-mTOR Inhibition in Cancer Immunotherapy

Differential Regulation of the Pro-Inflammatory Biomarker, YKL-40/CHI3L1, by PTEN/ Phosphoinositide 3-kinase and JAK2/STAT3 Pathways in Glioblastoma

Molecular and Clinical Characterization of CD163 Expression via Large-Scale Analysis in Glioma

TGFb in T Cell Biology and Tumor Immunity: Angel or Devil?

Expression of Osteoarthritis Marker YKL-39 is Stimulated by Transforming Growth Factor Beta (TGF-beta) and IL-4 in Differentiating Macrophages

Microglia-derived TGF-beta as an Important Regulator of Glioblastoma Invasion-an Inhibition of TGF-Beta-Dependent Effects by shRNA Against Human TGF-Beta Type II Receptor

Immune microenvironment of gliomas

Expression of Immune Checkpoints and T Cell Exhaustion Markers in Early and Advanced Stages of Colorectal Cancer

Cor, R-value of Spearman's correlation. (*P < 0.05; **P < 0.01; ***P < 0.001).Liu et al.

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.611038/ full#supplementary-material