Abstract
Glioblastoma multiforme (GBM) is the most aggressive primary tumor of the central nervous system. As biomedicine advances, the researcher has found the development of GBM is closely related to immunity. In this study, we evaluated the GBM tumor immunoreactivity and defined the Immune-High (IH) and Immune-Low (IL) immunophenotypes using transcriptome data from 144 tumors profiled by The Cancer Genome Atlas (TCGA) project based on the single-sample gene set enrichment analysis (ssGSEA) of five immune expression signatures (IFN-γ response, macrophages, lymphocyte infiltration, TGF-β response, and wound healing). Next, we identified six immunophenotype-related long non-coding RNA biomarkers (im-lncRNAs, USP30-AS1, HCP5, PSMB8-AS1, AL133264.2, LINC01684, and LINC01506) by employing a machine learning computational framework combining minimum redundancy maximum relevance algorithm (mRMR) and random forest model. Moreover, the expression level of identified im-lncRNAs was converted into an im-lncScore using the normalized principal component analysis. The im-lncScore showed a promising performance for distinguishing the GBM immunophenotypes with an area under the curve (AUC) of 0.928. Furthermore, the im-lncRNAs were also closely associated with the levels of tumor immune cell infiltration in GBM. In summary, the im-lncRNA signature had important clinical implications for tumor immunophenotyping and guiding immunotherapy in glioblastoma patients in future.
Introduction
Glioblastoma multiforme (GBM) is the most aggressive type of primary brain tumor in adults, with a median survival of 14.6 months (). The emergence of tumor immunotherapy has revolutionized GBM treatment and its success is highly dependent on the development and activation of immune cells in the host microenvironment (; ). In the GBM microenvironment, the non-neoplastic cells are mainly from the innate immune system, which can interact with neoplastic tumor cells and play an important role in tumor growth and progression (; ; ). Therefore, evaluation of GBM tumor immunoreactivity is critical in determining personalized treatment.
Long non-coding RNAs (lncRNAs) are defined as non-coding RNAs of more than 200 nt in length (). The discovery of lncRNAs has revealed a new dimension to the pathological processes of many diseases, including the occurrence and development of cancer (; ). Moreover, recent studies showed that lncRNAs play an important role in tumor immune escape (; ; ). For example, UCA1 is able to promote proliferation, migration, immune escape, and inhibit apoptosis in gastric cancer (); SNHG1 can regulate the differentiation of Treg cells and affect the immune escape of breast cancer (). Besides, immune-associated lncRNAs can also serve as improving prognosis and immunotherapy response biomarkers (; ). Therefore, identification of lncRNA biomarkers for tumor immunoreactivity may provide new insights into the treatment of GBM patients.
In this study, we systemically characterized the GBM tumor immune microenvironment in the TCGA GBM cohort. Moreover, we defined the GBM Immune-High (IH) and Immune-Low (IL) subtype based on five immune expression signatures including macrophages, lymphocyte infiltration, TGF-β response, IFN-γ response, wound healing. Furthermore, we identified six immunophenotype-related lncRNA signatures (im-lncRNAs, including USP30-AS1, HCP5, PSMB8-AS1, AL133264.2, LINC01684, and LINC01506) using the minimum redundancy maximum relevance (mRMR) feature selection method and the random forest model. The im-lncRNAs showed good performance in distinguishing tumor immunophenotypes and were closely associated with the levels of tumor immune cell infiltration. These results suggested the im-lncRNAs had the promising potential for clinical diagnosis of GBM immunophenotypes.
Materials and Methods
Data Acquisition and Pre-processing
All glioblastoma multiforme tissue samples were obtained from the surgical resection tissue of GBM patients (n = 10), non-tumor brain tissue was used as the negative control group (n = 5). The tissue samples were stored in liquid nitrogen. All patients have signed informed consent, and the study was supervised and approved by the Ethics Committee of The First Affiliated Hospital of Harbin Medical University.
The Cancer Genome Atlas (TCGA) level 3 gene/lncRNA expression data, and clinical data of GBM (n = 149, 144 cancer samples, 5 normal samples) were obtained from the Genomic Data Commons (GDC, available at https://www.cancer.gov/tcga). Two independent datasets GSE79671 () and GSE121810 () were used for the validation of im-lncRNAs. For the gene/lncRNA expression data, we removed the genes/lncRNAs that were not expressed over 70% of the samples. The remaining 18,094 genes and 18,567 lncRNAs were used for subsequent analysis.
Total RNA Extraction and Quantitative Real-Time PCR
According to the manufacturer’s instructions, total RNA was extracted from the GBM tissues and non-tumor brain tissues using TRIzol Reagent (Invitrogen, Carlsbad, CA, United States). The concentration of the total RNA was detected by spectrophotometer (Thermo ScientificTM NanoDrop 2000c). Total RNA (1000 ng) was reverse transcribed into cDNA using qPCR RT Kit (TOYOBO, Japan). The relative level of lncRNAs to the housekeeping gene GAPDH was determined by qRT-PCR using FastStart Universal SYBR Green Master (ROX) (Roche, Germany). All primers used in this study is showed in Supplementary Table 1. Analysis between the two groups was performed by an unpaired t-test, P < 0.05 was considered statistically significant.
Identification of Tumor Immune Subtypes of GBM
Based on five immune expression signatures reorganized by Vesteinn et al. () including IFN-γ response (), macrophages/monocytes (), overall lymphocyte infiltration (dominated by T and B cells) (), TGF-β response (), wound healing (), we evaluated the enrichment scores (ESs) of GBM samples using the single-sample gene set enrichment analysis (ssGSEA) (). The ssGSEA was based on the R package “GSVA.” Furthermore, we used the ESs of immune expression signatures to perform a consensus clustering on 149 cancer samples using the R package “ConsensusClusterPlus” ().
Evaluation of Tumor Purity, Tumor-Infiltrating Immune Cells, and Cytolytic Activity
The tumor purity of corresponding TCGA samples was evaluated using the ESTIMATE score calculated by the R package “ESTIMATE” (). The higher ESTIMATE score, the lower tumor purity. The tumor immune cell infiltration levels were estimated based on the gene expression profile by Tumor Immune Estimation Resource (TIMER) (). Here, six tumor-infiltrating immune cells (B cells, CD4 T cells, CD8 T cells, macrophages, neutrophils, and myeloid dendritic cells) were considered. Cytolytic activity (CYT) was calculated as the geometric mean of the GZMA gene and PRF1 gene (as expressed in FPKM) ().
Differential Expression Analysis of lncRNAs
We first calculated the log2(fold change) (log2(FC)) of each lncRNA between the IL and normal samples, and between the IH and normal samples, respectively. Then we scaled the expression level (log2FPKM) of each lncRNA and into a Z-score. Next, we compared lncRNA expression differences between the IL and normal samples, and between the IH and normal samples, using the Student’s t-test, respectively. The P-values were corrected using the Benjamini-Hochberg adjustment. The lncRNAs with FDR < 0.01 and | log2FC| > 2 were considered as the differentially expressed lncRNAs.
Identification of im-lncRNAs
We first divided the GBM cancer samples into three parts (two “training” sets and one “test” set) to apply three-fold cross-validation. Next, we screened the lncRNA features with minimal redundancy using the mRMR feature selection method in the training set (). Further, we trained a random forest model based on the top 5% mRMR lncRNA features. The performance of the random forest model was assessed through prediction making in the test set and the computation of the balanced error rate (BER). For a more robust estimation of the BER, three-fold cross-validation was applied 1,000 times and for each run, randomized data were used as the negative control. The signature size, for which no more improvement of the BER was observed (6 features signature size), was selected as the final size. This process generated 3 × 1000 output signatures. The distance (D) between these signatures was defined as ():
where cor represents the Spearman’s correlation coefficient (Rho); F1i to the ith feature from signature 1 and F2i to the ith feature from signature 2 after sorting the features to maximize the sum of the Rho based on the changes in the Gini index. For each signature, the sum of its pairwise distance to all other output signatures was computed, and the signature with the smallest sum was assumed to be the most representative and chosen as the final lncRNA signature (im-lncRNA).
Construction of im-lncScore
To conveniently evaluate the GBM tumor immunophenotypes, we constructed the im-lncScore. Firstly, we applied principal component analysis (PCA) on the Z-scores of im-lncRNAs. Then the first component was used as the final im-lncScore for the cancer samples.
Analysis of Association Between im-lncRNAs and Tumor Immune Cell Infiltration
Firstly, we calculated the median infiltration levels of each immune cell; if the sample infiltration level was higher than the median level, the sample was defined as a high immune infiltration sample; if the sample infiltration level was lower than the median level, the sample was defined as a low immune infiltration sample. Then, the univariate logistic regression was performed to assess the association between each im-lncRNA and the infiltration levels of each immune cell. The OR, 95% confidence level of the OR, and P-values were calculated for each immune cell. The logistic regression was based on the R package “epiDisplay.”
Identification of Co-expressed Genes With im-lncRNAs
Based on the expression profiles of im-lncRNAs and genes, we calculated the Spearman’s correlation coefficient (Rho) between im-lncRNAs and genes. The raw P-values (Pr) were adjusted by multiple hypotheses using a permutation method. For each gene, the expression value (FPKM) was held consistent, and 1,000 random im-lncRNAs were used to perform the same Spearman’s correlation test, generating a set of 1,000 permutation P-values (Pp). Finally, an empirical P-value (Pe) was corrected using the following formula (which introduces a pseudo-count of 1). The gene with the Rho>0.6 and Pe <0.01 were treated as the co-expressed genes of im-lncRNAs.
Functional Enrichment Analysis
To annotated the biological functions of im-lncRNAs, we performed functional enrichment analysis on the co-expressed genes of im-lncRNAs using Metascape (). For each co-expressed gene list, pathway and process enrichment analysis have been carried out with the following ontology sources: KEGG Pathway, GO Biological Processes, Reactome Gene Sets, Canonical Pathways, and Hallmark Gene Sets.
Results
Characterizing the Immune Microenvironment of GBM
We analyzed 149 GBM RNA-seq expression profiles from TCGA. To evaluate the tumor immune activity, we used a previously described technique employing ssGSEA () based on the five immune expression signatures reorganized by Vesteinn et al. () including IFN-γ response (), lymphocyte infiltration (), macrophages/monocytes (), TGF-β response (), wound healing (). The result showed that there were higher ESs of all immune expression signatures in cancer than in normal samples (Figure 1A, IFN-γ P = 1.01e-03, leukocyte infiltration P = 1.83e-04, macrophages P = 5.80e-08, TGF-β P = 1.55e-05, and wound healing P = 3.33e-08). Moreover, based on the ESs of immune expression signatures, we subclassified the cancer samples using the consensus clustering method. The analysis resulted in 2 robust clusters: C1 and C2. Notably, the ESs of IFN-γ, leukocyte infiltration, macrophages in C1 were significantly higher than in C2 (Figure 1A, IFN-γ P = 1.46e-29, leukocyte infiltration P = 3.43e-18, and macrophages P = 2.60e-16). And, there was no significant difference in the ESs of TGF-β (P = 3.46e-01) and wound healing (P = 1.09e-01) between the two clusters. Furthermore, we evaluated the levels of tumor purity and tumor-infiltrating immune cells between the two clusters. There were lower tumor purity (P = 4.75e-16) and higher percent of tumor-infiltrating immune cells (B cell P = 1.48e-05, T cell CD4 P = 8.26e-4, Neutrophil P = 3.58e-06, and macrophage P = 7.14e-10) in C1 than C2 (Figure 1A). Therefore, we annotated the C1 sample was the Immune-High (IH) subtype, and the C2 sample was the Immune-Low (IL) subtype.
FIGURE 1
To further verify the levels of immune activation in different immune subtypes, we examined the expression levels of common effectors of immunity, such as granzyme A (GZMA), granzyme B (GZMB), and perforin (PRF1) (Figure 1B; ) and the immune cytolytic activity (CYT, an indicator of tumor local immunity, Figure 1C; ). Remarkably, these effectors of immunity and CYT were much higher in the IH subtype compared with the IL subtype.
Glioblastoma multiforme can be subclassified into distinct molecular subtypes based on their expression profiles: classical (CL), mesenchymal (ME), neural (NE), and proneural (PN) (; ). Here, we also enriched the tumor immune subtypes into the GBM molecular subtypes using Fisher’s exact test. The previous study indicated ME GBM was the most immunogenic among the four subclasses while the PN subtype was the least immunogenic (). Our result also showed that ME GBM was significantly enriched in the IH subtype, while CL and PN GBM tumors were significantly enriched in the IL subtype (Figure 1D). Besides, by analyzing the survival of ME GBM patients between IL and IH subtypes, we found the survival of ME with IH patients was significantly better than ME with IL patients (Figure 1E).
Identification of Immunophenotype-Related lncRNA Biomarkers in GBM
lncRNA, an emerging biomarker, plays an important role in tumor immune regulation (). However, few studies focus on the ability of lncRNA in tumor immunophenotyping. To identify the immunophenotype-related lncRNA biomarkers (im-lncRNAs), we first characterized the differentially expressed lncRNAs between the IH/IL and normal samples, respectively (FDR < 0.01 and | log2FC| > 2, see section “Materials and Methods,” Figure 2A). We identified 261 “Constitutive” lncRNAs differentially expressed in both immune subtypes (142 upregulated and 119 downregulated), 145 “IH-specific” lncRNAs only differentially expressed in IH subtype (72 upregulated and 73 downregulated), and 70 “IL-specific” lncRNAs only differentially expressed in IL subtype (55 upregulated and 7 downregulated, Figure 2B).
FIGURE 2
Next, we applied a machine learning method in differentially expressed lncRNAs to identify the im-lncRNAs (Figure 2C). Firstly, under three-fold cross-validation (dividing 144 cancer samples into three parts, two “training” sets [96 samples], and one “test” set [48 samples]), the mRMR feature selection method was used to establish a small signature with minimal redundancy and selected the top 5% lncRNA features to train the random forest models. Next, in the test set, the balanced-error rate (BER) was calculated to evaluate the model performance. For a more robust estimation of the BER, three-fold cross-validation was applied 1,000 times. In each run, randomized data were used as the negative control. The signature size, for which no more improvement of the BER was observed (6 features signature size), was selected as the final size (Figure 2D). This pipeline generated 3 × 1000 output signatures and the signature with the minimum distance summed was assumed to be the most representative (see section “Materials and Methods”). Based on the approach, we identified 6 im-lncRNAs (USP30-AS1, HCP5, PSMB8-AS1, AL133264.2, LINC01684, LINC01506). Notably, USP30-AS1, HCP5, PSMB8-AS1, and LINC01506 were “IH-specific” lncRNAs, and AL133264.2, LINC01684 were “Constitutive” lncRNAs. The expression levels of all im-lncRNAs were significantly higher in IH than IL samples (USP30-AS1 P = 1.12e-18, HCP5 P = 8.07e-15, PSMB8-AS1 P = 1.15e-15, AL133264.2 P = 4.07e-10, LINC01684 P = 1.50e-11, and LINC01506 P = 3.00e-10, Figure 2E). Besides, all of 6 im-lncRNAs were also upregulated in GBM cancer than normal samples, which had been validated by RT-qPCR in five non-tumor brain tissues and ten GBM tissues (Figure 2F). To ensure that the recognized im-lncRNAs were robust, we also employed the same way on an independent dataset [GSE79671 ()]. The results showed that six im-lncRNAs closely associated with the GBM immunophenotypes and four of the six (USP30-AS1, HCP5, AL133264.2, and LINC01506) were consistent with the im-lncRNAs identified in the TCGA GBM cohort (Supplementary Figures 1A–C).
Evaluation of GBM Tumor Immunophenotyping Efficacy of im-lncRNAs
To further evaluate the relationship between im-lncRNAs and GBM immunophenotyping, we transformed the individual expression values of the im-lncRNAs into a score (im-lncScore) by applying a principal component analysis (PCA). We assessed the potential of the im-lncScore to predict GBM immunophenotypes in 144 cancer samples. Compared with the immune expression signatures, the im-lncScore also showed a promising performance. An AUC of 0.928 (95% CI, 0.87–0.97) suggested a predictive value for the im-lncScore (Figure 3A). Moreover, the optimal cutoff point determined by the ROC curve analysis was 0.0 (95% CI, 0.83–0.90). We also found that the im-lncScores of IH samples were usually greater than the optimal cutoff, while the opposite was observed for the IL samples (Figure 3B).
FIGURE 3
Besides, we also validated the immunophenotyping ability of im-lncScore in an independent dataset [GSE121810 ()]. The dataset included 29 GBM samples. We first calculated the im-lncScore to subclassify the GBM samples into IH/IL subtype. 13 IH and 16 IL samples were identified in the dataset (Figure 3C). Next, we also evaluated the ESs of five immune expression signatures and infiltration levels of tumor immune cells. As described above, there were higher ESs of IFN-γ, leukocyte infiltration and macrophages signatures, and higher levels of tumor-infiltrating immune cells in IH than IL samples (Figure 2C). These results suggested that the im-lncScore does not require a complex algorithm to effectively subclassify the GBM tumor immunophenotypes, which also further indicated the important role of im-lncRNAs in GBM tumor immunity.
Im-lncRNAs Are Associated With the Tumor Immune Cell Infiltration
To evaluate whether the im-lncRNAs associated with the levels of tumor immune cell infiltration, we first subclassified the cancer samples into high and low immune infiltration groups by comparing the sample immune infiltration levels to the median immune infiltration level of each immune cell. And then, the univariate logistic regression was performed based on the six im-lncRNAs expression value and im-lncScore. We found that the im-lncRNAs significantly correlated with the infiltration level of multiple immune cells (Figures 4A–F). Notably, the im-lncScore also showed the significantly correlation with multiple immune cell infiltration levels (except for T cell CD8). Besides, HCP5 and PSMB8-AS1 have been demonstrated could be the tumor-infiltrating immune-related lncRNA signature of non-small cell lung cancer and closely associated with outcome and immune cell infiltrates (). These results suggested that the im-lncRNAs played crucial roles in the tumor immune infiltration.
FIGURE 4
The Functional Enrichment Analysis of im-lncRNAs
To further explore the biological functions of im-lncRNAs, we identified the co-expressed genes with the im-lncRNAs using the Spearman’s correlation test. The P-values were adjusted by multiple hypotheses. A total of 459 co-expressed im-lncRNA-gene pairs were identified (Figure 5A). Furthermore, we performed functional enrichment analysis on the co-expressed genes using the Metascape (). The result showed that the functions of co-expressed genes of each im-lncRNAs were all significantly enriched in the immune-related terms, such as hallmark interferon-gamma response (M5913), myeloid leukocyte activation (GO:0002274), tumor necrosis factor superfamily cytokine production (GO:0071706), ER-Phagosome pathway (R-HSA-1236974), etc. (Figure 5B). Moreover, we also found the im-lncRNAs were closely correlated with the GBM-related immune pathways (). For instance, the HCP5 and PSMB8-AS1 were related to the Antigen Processing and Presentation, Antimicrobials, and Natural Killer Cell Cytotoxicity; AL133264.2 was related to Interleukins; the LINC01684 and USP30-AS1 were related to Antimicrobials. These results further validated the important role of im-lncRNAs in the GBM immune regulation.
FIGURE 5
Discussion
Accumulating evidence suggests that lncRNA serves as a specific molecular marker for tumor immunoreactivity (; ). In this study, we analyze the role of lncRNAs in IH and IL tumor immunophenotypes. Moreover, we identify im-lncRNAs based on the machine learning method. Furthermore, we construct an im-lncScore using the expression value of im-lncRNAs. The im-lncScore shows a good performance for distinguishing the GBM tumor immunophenotypes (AUC = 0.928, 95%CI: 0.885–0.970). The im-lncScore does not need a complex algorithm to effectively reflect the patient tumor immunoreactivity. Furthermore, these im-lncRNAs are also closely associated with the levels of tumor immune cell infiltration. This evidence indicates that the im-lncRNAs have the potential to be an important indicator for future clinical diagnosis of GBM immunophenotypes. However, these results are still at the level of the initial calculation, so to ensure accuracy the biology experiments are required to further validate the role of im-lncRNAs. Besides, due to the limited scale, we only use the TCGA data to train our models. Therefore, as the scale of data increases, we will continue to validate the efficiency of im-lncRNAs in GBM.
In summary, im-lncRNAs play an important role in tumor immunophenotyping. Identification of GBM immunophenotypes will provide us a novel insight to improve the therapeutic strategy of GBM. Therefore, the im-lncRNAs has the promising potential for clinical diagnosis of GBM immunophenotypes in the future.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: The Cancer Genome Atlas (TCGA) project (https://www.cancer.gov/tcga), GEO (GSE121810, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE121810).
Author contributions
MG, XC, and SZ conceived and designed the experiments. MG, XW, DH, and EL analyzed the data. JZ, LW, and QY collected the data. QJ and JW validated the method and data. CZ checked the writing standard of the manuscript. SZ, XC, and MG wrote this manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (81972363) and Research Fund for the Postdoctoral Science Foundation of China (2017M620119).
Acknowledgments
The authors gratefully thank the TCGA Research Network for providing data for this work.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.604655/full#supplementary-material
References
1
BarbieD. A.TamayoP.BoehmJ. S.KimS. Y.MoodyS. E.DunnI. F.et al (2009). Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1.Nature462108–112. 10.1038/nature08460
2
BeckA. H.EspinosaI.EdrisB.LiR.MontgomeryK.ZhuS.et al (2009). The macrophage colony-stimulating factor 1 response signature in breast carcinoma.Clin. Cancer Res.15778–787. 10.1158/1078-0432.CCR-08-1283
3
CalabroA.BeissbarthT.KunerR.StojanovM.BennerA.AsslaberM.et al (2009). Effects of infiltrating lymphocytes and estrogen receptor on gene expression and prognosis in breast cancer.Breast Cancer Res. Treat.11669–77. 10.1007/s10549-008-0105-3
4
CeccarelliM.BarthelF. P.MaltaT. M.SabedotT. S.SalamaS. R.MurrayB. A.et al (2016). Molecular profiling reveals biologically discrete subsets and pathways of progression in diffuse glioma.Cell164550–563. 10.1016/j.cell.2015.12.028
5
ChangH. Y.SneddonJ. B.AlizadehA. A.SoodR.WestR. B.MontgomeryK.et al (2004). Gene expression signature of fibroblast serum response predicts human cancer progression: similarities between tumors and wounds.PLoS Biol.2:E7. 10.1371/journal.pbio.0020007
6
CloughesyT. F.MochizukiA. Y.OrpillaJ. R.HugoW.LeeA. H.DavidsonT. B.et al (2019). Neoadjuvant anti-PD-1 immunotherapy promotes a survival benefit with intratumoral and systemic immune responses in recurrent glioblastoma.Nat. Med.25477–486. 10.1038/s41591-018-0337-7
7
DaudA. I.WolchokJ. D.RobertC.HwuW. J.WeberJ. S.RibasA.et al (2016). Programmed death-ligand 1 expression and response to the anti-programmed death 1 antibody pembrolizumab in melanoma.J. Clin. Oncol.344102–4109. 10.1200/JCO.2016.67.2477
8
DoucetteT.RaoG.RaoA.ShenL.AldapeK.WeiJ.et al (2013). Immune heterogeneity of glioblastoma subtypes: extrapolation from the cancer genome atlas.Cancer Immunol. Res.1112–122. 10.1158/2326-6066.CIR-13-0028
9
EnglerJ. R.RobinsonA. E.SmirnovI.HodgsonJ. G.BergerM. S.GuptaN.et al (2012). Increased microglia/macrophage gene expression in a subset of adult and pediatric astrocytomas.PLoS One7:e43339. 10.1371/journal.pone.0043339
10
FengX.SzulzewskyF.YerevanianA.ChenZ.HeinzmannD.RasmussenR. D.et al (2015). Loss of CX3CR1 increases accumulation of inflammatory monocytes and promotes gliomagenesis.Oncotarget615077–15094. 10.18632/oncotarget.3730
11
HanchuanP.FuhuiL.DingC. (2005). Feature selection based on mutual information criteria of max-dependency, max-relevance, and min-redundancy.IEEE Trans. Pattern Anal. Mach. Intellig.271226–1238. 10.1109/TPAMI.2005.159
12
HuF.DzayeO.HahnA.YuY.ScavettaR. J.DittmarG.et al (2015). Glioma-derived versican promotes tumor expansion via glioma-associated microglial/macrophages Toll-like receptor 2 signaling.Neuro Oncol.17200–210. 10.1093/neuonc/nou324
13
JeschkeJ.BizetM.DesmedtC.CalonneE.DedeurwaerderS.GaraudS.et al (2017). DNA methylation-based immune response signature improves patient diagnosis in multiple cancers.J. Clin. Invest.1273090–3102. 10.1172/JCI91095
14
JinK. T.YaoJ. Y.FangX. L.DiH.MaY. Y. (2020). Roles of lncRNAs in cancer: focusing on angiogenesis.Life Sci.252:117647. 10.1016/j.lfs.2020.117647
15
KaffesI.SzulzewskyF.ChenZ.HertingC. J.GabanicB.Velazquez VegaJ. E.et al (2019). Human Mesenchymal glioblastomas are characterized by an increased immune cell presence compared to proneural and classical tumors.Oncoimmunology8:e1655360. 10.1080/2162402X.2019.1655360
16
LekM.KarczewskiK. J.MinikelE. V.SamochaK. E.BanksE.FennellT.et al (2016). Analysis of protein-coding genetic variation in 60,706 humans.Nature536285–291. 10.1038/nature19057
17
LiT.FanJ.WangB.TraughN.ChenQ.LiuJ. S.et al (2017). TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells.Cancer Res.77e108–e110. 10.1158/0008-5472.CAN-17-0307
18
LiY.JiangT.ZhouW.LiJ.LiX.WangQ.et al (2020). Pan-cancer characterization of immune-related lncRNAs identifies potential oncogenic biomarkers.Nat. Commun.11:1000. 10.1038/s41467-020-14802-2
19
MandalR.SenbabaogluY.DesrichardA.HavelJ. J.DalinM. G.RiazN.et al (2016). The head and neck cancer immune landscape and its immunotherapeutic implications.JCI Insight1:e89829. 10.1172/jci.insight.89829
20
Martens-UzunovaE. S.BottcherR.CroceC. M.JensterG.VisakorpiT.CalinG. A. (2014). Long noncoding RNA in prostate, bladder, and kidney cancer.Eur. Urol.651140–1151. 10.1016/j.eururo.2013.12.003
21
MontiS.TamayoP.MesirovJ.GolubT. (2003). Consensus clustering: a resampling-based method for class discovery and visualization of gene expression microarray data.Mach. Learn.5291–118. 10.1023/A:1023949509487
22
PardollD. M. (2012). The blockade of immune checkpoints in cancer immunotherapy.Nat. Rev. Cancer12252–264. 10.1038/nrc3239
23
PeiX.WangX.LiH. (2018). LncRNA SNHG1 regulates the differentiation of Treg cells and affects the immune escape of breast cancer via regulating miR-448/IDO.Int. J. Biol. Macromol.118(Pt A), 24–30. 10.1016/j.ijbiomac.2018.06.033
24
RooneyM. S.ShuklaS. A.WuC. J.GetzG.HacohenN. (2015). Molecular and genetic properties of tumors associated with local immune cytolytic activity.Cell16048–61. 10.1016/j.cell.2014.12.033
25
SunJ.ZhangZ.BaoS.YanC.HouP.WuN.et al (2020). Identification of tumor immune infiltration-associated lncRNAs for improving prognosis and immunotherapy response of patients with non-small cell lung cancer.J. Immunother. Cancer8:e000110. 10.1136/jitc-2019-000110
26
TeschendorffA. E.GomezS.ArenasA.El-AshryD.SchmidtM.GehrmannM.et al (2010). Improved prognostic classification of breast cancer defined by antagonistic activation patterns of immune response pathway modules.BMC Cancer10:604. 10.1186/1471-2407-10-604
27
UrupT.StaunstrupL. M.MichaelsenS. R.Vitting-SeerupK.BennedbaekM.ToftA.et al (2017). Transcriptional changes induced by bevacizumab combination therapy in responding and non-responding recurrent glioblastoma patients.BMC Cancer17:278. 10.1186/s12885-017-3251-3
28
VerhaakR. G.HoadleyK. A.PurdomE.WangV.QiY.WilkersonM. D.et al (2010). Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1.Cancer Cell1798–110. 10.1016/j.ccr.2009.12.020
29
WangC. J.ZhuC. C.XuJ.WangM.ZhaoW. Y.LiuQ.et al (2019). The lncRNA UCA1 promotes proliferation, migration, immune escape and inhibits apoptosis in gastric cancer by sponging anti-tumor miRNAs.Mol. Cancer18:115. 10.1186/s12943-019-1032-0
30
WolfD. M.LenburgM. E.YauC.BoudreauA.van ’t VeerL. J. (2014). Gene co-expression modules as clinically relevant hallmarks of breast cancer diversity.PLoS One9:e88309. 10.1371/journal.pone.0088309
31
YoshiharaK.ShahmoradgoliM.MartinezE.VegesnaR.KimH.Torres-GarciaW.et al (2013). Inferring tumour purity and stromal and immune cell admixture from expression data.Nat. Commun.4:2612. 10.1038/ncomms3612
32
ZhaiX.ZhaoJ.WangY.WeiX.LiG.YangY.et al (2018). Bibliometric analysis of global scientific research on lncRNA: a swiftly expanding trend.Biomed. Res. Int.2018:7625078. 10.1155/2018/7625078
33
ZhouM.ZhangZ.ZhaoH.BaoS.ChengL.SunJ. (2018). An immune-related six-lncRNA signature to improve prognosis prediction of glioblastoma multiforme.Mol. Neurobiol.553684–3697. 10.1007/s12035-017-0572-9
34
ZhouQ.ChungA. C.HuangX. R.DongY.YuX.LanH. Y. (2014). Identification of novel long noncoding RNAs associated with TGF-beta/Smad3-mediated renal inflammation and fibrosis by RNA sequencing.Am. J. Pathol.184409–417. 10.1016/j.ajpath.2013.10.007
35
ZhouY.ZhouB.PacheL.ChangM.KhodabakhshiA. H.TanaseichukO.et al (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets.Nat. Commun.10:1523. 10.1038/s41467-019-09234-6
Summary
Keywords
long non-coding RNA, biomarker, immunophenotype, machine learning, glioblastoma multiforme
Citation
Gao M, Wang X, Han D, Lu E, Zhang J, Zhang C, Wang L, Yang Q, Jiang Q, Wu J, Chen X and Zhao S (2021) A Six-lncRNA Signature for Immunophenotype Prediction of Glioblastoma Multiforme. Front. Genet. 11:604655. doi: 10.3389/fgene.2020.604655
Received
10 September 2020
Accepted
26 October 2020
Published
13 January 2021
Volume
11 - 2020
Edited by
Suman Ghosal, National Institutes of Health, United States
Reviewed by
Provas Das, Baylor University, United States; Jayshree Advani, National Institutes of Health, United States
Updates
Copyright
© 2021 Gao, Wang, Han, Lu, Zhang, Zhang, Wang, Yang, Jiang, Wu, Chen and Zhao.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xin Chen, chenxin_tracy@yeah.netShiguang Zhao, guangsz@hotmail.com
This article was submitted to Computational Genomics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.