A Composite Bioinformatic Analysis to Explore Endoplasmic Reticulum Stress-Related Prognostic Marker and Potential Pathogenic Mechanisms in Glioma by Integrating Multiomics Data

Fan, Xin; Nie, Xiyi; Huang, Junwen; Zhang, Lingling; Wang, Xifu; Lu, Min

doi:https://doi.org/10.1155/2022/9886044

Journal of Oncology

On this page

Abstract Introduction Methods Results Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

Prognostic Models Based on Machine Learning for Clinical Cancer Research 2022

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 9886044 | https://doi.org/10.1155/2022/9886044

A Composite Bioinformatic Analysis to Explore Endoplasmic Reticulum Stress-Related Prognostic Marker and Potential Pathogenic Mechanisms in Glioma by Integrating Multiomics Data

Xin Fan,^1,2Xiyi Nie,³Junwen Huang,⁴Lingling Zhang,⁵Xifu Wang,¹and Min Lu¹

Academic Editor: Quan Cheng

Received07 Apr 2022

Revised18 Jul 2022

Accepted14 Aug 2022

Published04 Oct 2022

Abstract

In recent years, abnormal endoplasmic reticulum stress (ERS) response, as an important regulator of immunity, may play a vital role in the occurrence, development, and treatment of glioma. Weighted correlation network analysis (WGCNA) based on six glioma datasets was used to screen eight prognostic-related differentially expressed ERS-related genes (PR-DE-ERSGs) and to construct a prognostic model. BMP2 and HEY2 were identified as protective factors (HR < 1), and NUP107, DRAM1, F2R, PXDN, RNF19A, and SCG5 were identified as risk factors for glioma (HR > 1). QRT-PCR further supported significantly higher DRAM1 and lower SCG5 relative mRNA expression in gliomas. Our model has demonstrated excellent performance in predicting the prognosis of glioma patients from numerous datasets. In addition, the model shows good stability in multiple tests. Our model also shows broad clinical promise in predicting drug treatment effects. More immune cells/processes in the high-risk population with poor prognosis illustrate the importance of the tumor immunosuppressive environment in glioma. The potential role of the HEY2-based competitive endogenous RNA (ceRNA) regulatory network in glioma was validated and revealed the possible important role of glycolysis in glioma ERS. IDH1 and TP53 mutations with better prognosis were strongly associated with the risk score and PR-DE-ERSGs expression in the model. mDNAsi was also closely related to the risk score and clinical characteristics.

1. Introduction

Gliomas are the most common primary intracranial tumor, accounting for 81% of malignant brain tumor [1]. It is a very invasive brain tumor, which causes significant mortality and morbidity [2, 3]. With the development of molecular biology techniques, much evidence shows that gliomas’ genetic and epigenetic features have changed significantly, facilitating many new diagnostic and therapeutic approaches, including targeted therapy and immunotherapy. Conventional therapies, including surgery, chemotherapy, and radiotherapy, have achieved limited improvements in the prognosis of glioma patients [4, 5]. Compared with other solid tumors, gliomas pose significant challenges to the development of novel tumor treatments due to biological factors such as the blood-brain barrier and unique tumor and immune microenvironments [6]. A remarkable level of genetic, epigenetic, and environmental heterogeneity exists within each glioma, providing multiple mechanisms of therapeutic resistance and forming this highly adaptable and resilient disease [7]. The current rapid development of in-silico technology provides new computational strategies for potential drug targets for many diseases [8–10]. Therefore, in order to improve the level of diagnosis and treatment of glioma and prognosis prediction, it is particularly necessary to identify effective biomarkers combined with clinical features.

The endoplasmic reticulum (ER) is the primary site for protein synthesis, processing, and trafficking, and many genetic and environmental impairments impede the cell’s ability to fold correctly and post-translationally modify secreted and transmembrane proteins in the ER, resulting in the accumulation of misfolded proteins in the ER [11]. The above situation is called ER stress (ERS). Beyond a tolerable threshold, the accumulation of misfolded proteins triggers the unfolded protein response (UPR) to improve the folding ability of ER proteins during transcription and translation [12]. Robust ERS responses have been documented in most major types of human cancer, including breast, pancreatic, lung, skin, prostate, brain, and liquid malignancies [11]. Hostile microenvironmental conditions within tumor masses, such as nutrient deprivation, oxygen limitation, high metabolic demand, and oxidative stress, disturb the protein folding capacity of the ER, thereby provoking a cellular state of ERS [13]. Sustained activation of ERS sensors endows malignant cells with greater tumorigenic, metastatic, and drug resistant capacity [14].

Many studies have also reported the pathogenic role of ERS in the gliomas’ initiation and progression and the potential therapeutic guiding value. ERS has been reported to be associated with self-renewal, differentiation, and drug resistance of glioma stem cells, and plays a vital role in the resistance of glioblastoma (GBM) to temozolomide toxicity (TMZ) [15–17]. General UPR activators or selective GRP78, ATF6, and PERK inducers have been detected to modulate cell proliferation and induce apoptosis in glioma cells [18]. Additionally, increased UPR activity strongly affects many intracellular metabolic pathways, which in turn shape the tumor microenvironment [18]. UPR favors increased 13C-glucose uptake and results in higher levels of lactate, alanine, and uridine diphosphate glucose, which are often associated with tumor aggressiveness [19]. Recent studies have shown that abnormal ERS response emerges as critical regulators of immune cell function in the tumor microenvironment, and it can subvert the protective function of innate immune cells in the tumor microenvironment to cripple the development of antitumor immunity [14]. Therefore, understanding the relationship between ERS and tumor microenvironment/immune cells and exploring the related mechanisms may become a new way to supplement and improve standard chemotherapy and immunotherapy in the clinic.

The competing endogenous RNA (ceRNA) hypothesis postulates that any RNA transcript with microRNA (miRNA)-response elements (MREs) can sequester miRNAs from other targets sharing the same MREs, thereby regulating their expression [20]. Therefore, cross-talk between RNAs, both coding and noncoding, through MREs forms a large-scale regulatory network across the transcriptome [21]. Mounting evidence has shown that various types of RNAs, including pseudogenes, long noncoding RNAs (lncRNAs), circular RNAs, and messenger RNAs, can function as ceRNAs in distinct physiological and pathophysiological states [22]. Recent studies on solid tumors and hematopoietic malignancies have shown that ceRNAs play significant roles in cancer pathogenesis by altering the expression of key tumorigenic or tumor-suppressive genes [23], including gastric cancer [24], lung adenocarcinoma [25], pancreatic cancer [26], and gallbladder [27]. While many studies investigate tumors, including gliomas, from the perspective of ERS, few studies have investigated gliomas from the perspective of ERS gene-related ceRNA. Therefore, exploring ERS-related ceRNA regulatory networks may contribute to the understanding of ERS-related biological roles in gliomas.

In this study, with the primary purpose of constructing a biomarker with high predictive performance and clinical value, we deeply analyzed ERS-related genes that are highly correlated with glioma progression and prognosis. To further elucidate the important role of ERS in glioma, a multi-omics analysis from the perspective of ERS-related immune microenvironment, cell stemness, and ceRNA regulatory network was also included in this study.

2. Material and Methods

2.1. Recruitment of Glioma Samples for Obtaining DE-ERSGs

The schematic overview of our research process is shown in Figure 1. Our research enrolled the RNA sequencing data and clinical information of glioma samples from 3 public databases, including The Cancer Genome Atlas database (TCGA), Gene Expression Omnibus (GEO), and Chinese Glioma Genome Atlas (CGGA). The TCGA and CGGA databases contain 703 samples (698 glioma and 5 adjacent normal brain tissues) and 1018 glioma samples, respectively. The following 4 datasets from the GEO database were also included: GSE108474 (376 astrocytoma, 67 oligodendroglioma, and 28 normal brain tissues), GSE4290 (26 astrocytoma, 81 glioblastoma, 50 oligodendroglioma, and 23 normal brain tissues), GSE4412 (85 glioma), and GSE43378 (5 astrocytoma, 32 glioblastoma, and 13 anaplastic oligoastrocytoma).

In the meantime, we obtained 5909 endoplasmic reticulum stress-related genes (ERSGs) with a relevance score ≥1 from the gene card database (GeneCards), which contains genomic, transcriptomic, proteomic, genetic, clinical, and functional information integrated from 150 web sources [28]. Subsequently, we obtained the RNA sequencing level of ERSGs in those 6 datasets (TCGA: 5828, GSE108474: 5353, GSE4290: 5353, GSE4412: 4748, GSE43378: 5351, CGGA: 5254).

To obtain the differentially expressed endoplasmic reticulum stress-related genes (DE-ERSGs) between gliomas and normal tissues, we set filter conditions for datasets from different databases (TCGA: |log 2 fold change| (|log 2FC|) > 1 and false discovery rate (FDR) < 0.05, GEO: FDR < 0.05), respectively. We intersected DE-ERSGs from TCGA, GSE108474, and GSE4290 datasets with ERSGs from GSE4412, GSE43378, and CGGA datasets to obtain common genes. This process was shown via the R “Venn” package.

2.2. Screening for DE-ERSGs Strongly Related to Glioma

Weighted gene coexpression network analysis (WGCNA) was used to explore the modules of highly correlated genes among samples for relating modules to external sample traits [29]. We identified 5, 3, and 2 gene modules by WGCNA for the TCGA, GSE108474, and GSE4290 datasets, respectively. The values of best soft power β used in these three WGCNA processes are 11, 3, and 1, respectively. Among all modules, the turquoise module of GSE4290 has the strongest correlation with glioma. Therefore, we extracted DE-ERSGs strongly related to glioma from the modules with the strongest positive correlation with glioma from the TCGA, GSE108474, and GSE4290 datasets, respectively. The R “Venn” package was again used to obtain the common genes of the three gene modules.

2.3. Screening PR-DE-ERSGs for Prognostic Model’s Construction Based on LASSO Regression

For subsequent survival analysis, we curated data from 6 datasets to extract glioma samples with full overall survival (OS) data. TCGA dataset’s 668, GSE4412 dataset’s 85, GSE43378 dataset’s 50, and CGGA dataset’s 983 glioma samples satisfy this requirement, as counted in Tables 1 and 2. Univariate Cox analyses identified differentially expressed endoplasmic reticulum stress-related genes associated with prognosis (PR-DE-ERSGs) based on TCGA and CGGA datasets data and filter threshold (), respectively. Common PR-DE-ERSGs from two datasets were extracted. Expression heatmap, forest plot reflecting univariate Cox regression results, and Kaplan–Meier survival curves based on TCGA data were all used to display the 14 common PR-DE-ERSGs.

The 668 TCGA samples were randomly equally divided into 2 subsets (training and test sets). Based on the RNA sequencing level of common PR-DE-ERSGs from the training set, the least absolute shrinkage and selection operator (LASSO) regression analysis was utilized to screen out highly related genes and minimize the risk of overfitting for screening signatures [30]. Eventually, 8 PR-DE-ERSGs and corresponding coefficients were obtained for constructing the prognostic model by calculating the optimal penalty parameter (λ). The risk score for each sample = sum of corresponding coefficients × matching gene’s RNA sequencing level.

2.4. Validation of Differential Expression of 8 Genes by External Data

To maintain the stability of the model, we tried to validate the differential expression of eight genes in the model between glioma and normal brain tissues using data from an external database. GEPIA is a web tool server for cancer and normal gene expression profiling and interactive analyses [31]. Finally, we used this database to compare the eight gene expression differences between glioma (glioblastoma multiforme + lower grade glioma) and normal brain tissues.

2.5. Model’s Performance and Stability Testing

After assigning a risk score for each sample, we got the median risk score. This median risk score directly split all samples into high-risk and low-risk groups. Risk curve and survival point plot were matched to show each sample’s risk score and survival status. Kaplan–Meier survival curve and receiver operator characteristic (ROC) curve were used to test the performance of the model to discriminate and predict prognosis, respectively. Univariate and multivariate Cox regression analyses were used to determine the stability of the model as a prognostic factor before and after adjustment for confounders. To maintain the stability of the test results, we repeated the above analysis based on 3 sets from the TCGA database, GSE4412 and GSE43378 sets from the GEO database, and CGGA set, respectively.

As one of the most popular machine learning algorithms for regression and classification tasks [32], the random forest (RF) algorithm is often used to determine the importance of eigengenes [33]. To further assess the relative importance of the eight PR-DE-ERSGs used to build the model in glioma, we ran the RF algorithm. The RF algorithm was run based on the minimum point of cross-validation error to further rank the importance of the 8 PR-DE-ERSGs. This step was performed synchronously in the TCGA, GSE108474, and GSE4290 cohorts. In addition, we also plotted ROC curves based on 8 PR-DE-ERSGs expression of glioma and normal samples from the TCGA dataset to evaluate their diagnostic value in glioma.

2.6. Deep Validation of Model’s Performance and Stability

The heatmap visually shows the clinical features matched to the risk score of each sample from TCGA. We observed differences in risk scores across subgroups of different clinical characteristics. Kaplan–Meier survival curves further tested the model’s ability to distinguish prognosis within each clinical subgroup.

We attempted to obtain glioma and normal brain tissues’ immunohistochemical (IHC) staining images of modeled 8 PR-DE-ERSGs from the Human Protein Atlas database. These images were used to compare further the expression differences of the 8 PR-DE-ERSGs at the protein expression level to verify the stability of the model.

Unfortunately, IHC staining images still failed to confirm the differential expression of BMP2, DRAM1, NUP107, and SCG5 between gliomas and normal tissues. Therefore, we further carried out Quantitative Real-time PCR (QRT-PCR) experiments to detect the mRNA expression levels of these four genes.

After obtaining the approval from the Ethics Committee of Shangrao People’s Hospital (2022-04-041) and patients’ informed consent, we obtained 6 matched pairs of gliomas and adjacent paracancerous tissues from the hospital. All primer sequences were listed in the Supplementary Table 2. Total RNA was isolated from tissues using the TransZol Up Plus RNA Kit (TRANS, Beijing, China). According to the manufacturer’s instructions, cDNA was synthesized by using EasyScript First-Strand cDNA Synthesis SuperMix (TRANS, Beijing, China). QRT-PCR was performed by using PerfectStart® Green qPCR SuperMix (TRANS, Beijing, China) on the Roche Light Cycler 96 Real-time Fluorescent Quantitative PCR System. The detection values of 4 genes were normalized to the β-actin gene with the 2^−ΔΔCt method. After calculating the relative mRNA expression values of the four genes, we compared their differences between glioma and normal paracancerous tissues using paired t-test.

2.7. Enrichment Analysis for Biological Functions and Pathways Related to Different Risk Groups

Gene Set Enrichment Analysis (GSEA) was utilized to perform Kyoto encyclopedia of genes and genomes (KEGG) and Gene ontology (GO) analyses to screen the signal pathways and biological functions involved in the different risk groups based on expression matrix data [34]. The R package “clusterProfiler” and gene sets “c2.cp.kegg.v7.4.symbols.gmt” and “c5.go.v7.4.symbols.gmt” were applied to run enrichment analysis [34].

2.8. Exploring the Potential Role of Tumor Immunosuppressive Environment/Immune Infiltrates Types

After calculating the overall immune and stromal cell scores for each sample using the R package “estimate,” we further quantified 16 immune cells and 13 immune function scores by single-sample gene set enrichment analysis (ssGSEA) based on R packages “GSEABase” and “gsva” [35, 36]. Correlation analysis was run between overall immune cell/stromal cell/16 immune cell/13 immune function scores and risk scores to explore potential relationships between risk scores and immune microenvironment. As reported previously, six types of immune infiltration were identified in human tumors, which corresponded from tumor-promoting to tumor-suppressing, respectively, namely, C1 (wound healing), C2 (INF-g dominant), C3 (inflammatory), C4 (lymphocyte depleted), C5 (immunologically quiet), and C6 (TGF-b dominant) [37]. Heatmaps showed differences in the distribution of different immune infiltrating subtypes between high- and low-risk groups. Boxplots were used to further visualize differences in risk scores between different immune infiltrating subtypes.

2.9. Exploring the Potential Role of Cell Stemness

The acquisition of progenitor-like, stem-cell-like features is strongly connected with cancer progression. Stemness, the attribute of self-renewal and differentiation from the cell of origin, can be described by two stemness indexes: a gene expression-based stemness index (mRNAsi) and a DNA methylation-based stemness index (mDNAsi), which are calculated by one-class logistic regression machine learning algorithm (OCLR) training on stem cell (ESC, embryonic stemcell; iPSC, induced pluripotent stem cell) classes and their differentiated ecto-, meso-, and endoderm progenitors [38]. Malta et al. calculated the mDNAsi for each TCGA sample based on OCLR-based transcriptomic and epigenetic signatures [38]. Correlation and difference analyses were again run to explore the potential association of mDNAsi with risk score/various clinical characteristics.

2.10. Exploring the Potential Role of HEY2-Related ceRNA Regulatory Network

Gene names of miRNAs for RNA sequencing data from the TCGA database (530 gliomas and 5 adjacent normal tissues) were matched with annotation files of mature miRNAs obtained from miRbase. Considering the critical role of hairy and enhancer of split-related with YRPW motif protein 2 (HEY2) in glioma, we predicted miRNAs and lncRNAs that may regulate HEY2 expression.

To improve the accuracy of the prediction results, we selected miRNAs that appeared more than 3 times in several target gene prediction programs of starBase. The coexpression network between all predicted miRNAs and HEY2 was plotted by Cytoscape (v3.7.2). TCGA data was used to further screen miRNAs with significant roles in glioma. The correlation analysis between miRNAs and HEY2, miRNAs difference analysis between glioma and normal brain tissue, and the Kaplan–Meier survival analysis of miRNAs were performed under the corresponding filtering conditions (correlation analysis: correlation coefficient < −0.2, ; difference analysis: |log 2FC| > 1 and ; Kaplan–Meier survival analysis: ). Eventually, hsa-miR-369-3p and hsa-miR-181a-5p were identified as candidate miRNAs upstream of HEY2.

Afterward, we again predicted potential lncRNAs upstream of hsa-miR-369-3p and hsa-miR-181a-5p, respectively. Similar methods and filtering conditions were used to screen candidate lncRNAs, respectively. Only hsa-miR-369-3p obtained candidate lncRNAs. Cytoscape was again used to map the regulatory network composed of 5 candidate lncRNAs, hsa-miR-369-3p, and HEY2. We again ran correlation analysis between lncRNAs and hsa-miR-369-3p (correlation coefficient < −0.35 and ), correlation analysis between lncRNAs and HEY2 (correlation coefficient < −0.15 and ) to filter lncRNAs as the finalized lncRNAs deeply. Two finalized lncRNAs were highlighted in yellow in regulatory networks.

2.11. Exploring the Potential Role of Mutation

After detecting and counting somatic gene mutation data in glioma samples in MAF files from the TCGA database, waterfall plots were used to visualize the corresponding results. The correlation analysis between risk score and tumor mutation burden (TMB), the difference analysis of TMB between different risk groups, and the Kaplan–Meier survival analysis of TMB were all used to explore the potential relationship between the model and TMB.

The statistical results show that IDH1/TP53 mutation, the most frequently mutated gene in glioma, may play an indispensable role in glioma. Therefore, after obtaining the effect of IDH1/TP53 mutation on prognosis, we further compared the difference in the risk score/8 PR-DE-ERSGs expression between IDH1/TP53 mutation and wild groups. Considering the predictive value of CD274 in immunotherapy effect, we also compared the CD274 expression difference between the two groups to explore the guiding value of IDH1/TP53 mutation in clinical treatment.

2.12. Evaluation of the Guiding Value of the Model in Therapy

At present, chemotherapy and immunotherapy are still the main drug treatment methods for glioma patients [39, 40]. Immunophenoscore (IPS) from The Cancer Immunome Atlas (TCIA) database and Tumor Immune Dysfunction and Exclusion (TIDE) indexes from the TIDE online website have previously been used as predictors of the immunotherapy efficacy [41, 42]. A lower TIDE score and higher IPS represent a better patient response to immunotherapy [41–43]. Therefore, we obtained 3 kinds of IPSs and TIDE/MSI/T cell dysfunction/T cell exclusion scores from these two websites as indicators of immunotherapy’s effect. Correlations between 3 IPSs and 8 PR-DE-ERSGs expressions, and correlations between TIDE indexes and 8 PR-DE-ERSGs expressions were analyzed to explore the model’s guiding value in immunotherapy efficacy.

The NCI-60 database containing 60 cancer cell lines from 9 different types of tumors was accessed through the CellMiner interface [35]. Then, the association between 8 PR-DE-ERSGs and 263 drugs approved by the FDA or in clinical trials was examined by correlation analyses, respectively.

2.13. Extending the Model’s Clinical Value

To construct a more accurate comprehensive nomogram to predict the survival probability of glioma patients, we used independent prognostic factors obtained by Cox regression analysis. The calibration curve and ROC curve at 1, 2, and 3 years were used for performance testing of this tool.

Recent studies have shown that immune checkpoint inhibitors (ICIs), N6-methyladenosine (m6A), and multidrug resistance-related genes play an important role in the diagnosis and therapy of gliomas [44–46]. The correlation between the expression levels of 46 ICIs/12 m6A/2 multidrug resistance-related genes and risk scores, and the expression differences of these genes between different risk groups were analyzed to further explore the model’s extended clinical value.

2.14. Statistical Analysis

This analysis was performed primarily in the R programming language (version 4.0.3) based on numerous packages including “limma,” “WGCNA,” “ComplexHeatmap,” “GSEABase,” “ggExtra,” “maftools,” “regplot,” and “ggplot2.” In data comparison between/among groups, student’s t-test/chi-square test or nonparametric test/fisher’s exact test was used when appropriate. Unless otherwise stated, statistical significance was set at .

3. Results

3.1. Recruitment of Glioma Samples for Obtaining DE-ERSGs

Supplementary Figures 1A and 1B show 1438 DE-ERSGs (up-regulated: 785; down-regulated: 653) are identified from the TCGA dataset based on the 5828 ERSGs’ RNA sequencing data. Similarly, there are 2828 DE-ERSGs (up-regulated: 1507; down-regulated: 1321; Supplementary Figures 2A and 2B) extracted from the GSE108474 dataset and 3375 DE-ERSGs (up-regulated: 1886; down-regulated: 1489; Supplementary Figures 3A and 3B) extracted from the GSE4290 dataset. Five hundred eighty-four common DE-ERSGs were obtained by intersecting DE-ERSGs from these 3 datasets and ERSGs from the remaining 3 datasets (Figure 2(a)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

(n)

(o)

(p)

(q)

(r)

(s)

3.2. Screening for DE-ERSGs Strongly Related to Glioma

After setting the best soft power β, the dynamic pruning trees merged the similar modules and matched each gene with the corresponding gene module based on TCGA, GSE108474, and GSE4290 datasets, respectively (Figures 3(a), 3(c), and 3(e)). The heatmaps in Figures 3(b), 3(d), and 3(f) show the correlation between gene modules identified by WGCNA and glioma status based on these three datasets, respectively. We extracted 24 common DE-ERSGs strongly related to glioma from TCGA grey module, GSE108474 blue module, and GSE4290 turquoise module using Venn diagram (Figure 2(b)).

(a)

(b)

(c)

(d)

(e)

(f)

3.3. Screening PR-DE-ERSGs for Prognostic Model’s Construction Based on LASSO Regression

Fourteen common PR-DE-ERSGs were obtained from 20 PR-DE-ERSGs in the TCGA dataset and 17 PR-DE-ERSGs in the CGGA dataset (Figure 2(c)). The expression levels of these 14 genes in glioma and normal brain tissues are shown in Figure 2(d). Except for CD22, MTX1, BMP2, and HEY2 (HR < 1), the remaining 10 genes (HR > 1) were identified as risk factors in the TCGA dataset, as shown in the forest plot (Figure 2(e)) and Kaplan–Meier survival curves (Figures 2(f)–2(s)). It can be seen that random grouping did not cause differences in the distribution of data between groups (Table 1, ). By setting the optimal value of λ by the LASSO regression analysis (Supplementary Figures 4A and 4B), we acquired 8 PR-DE-ERSGs to construct the model, eventually. The corresponding coefficients of the 8 PR-DE-ERSGs are detailed in Supplementary Table 1.

3.4. Validation of Differential Expression of 8 Genes by External Data

Supplementary Figures 5A–5H show the differential analysis results of 8 genes from the GEPIA database. Consistent with our previous findings, the expressions of NUP107, DRAM1, RNF19A, PXDN, HEY2, F2R, and BMP2 were up-regulated in gliomas (glioblastoma multiforme + lower grade glioma). The stable expression trends of these genes all support the good stability of our model.

3.5. Model’s Performance and Stability Testing

Samples from all three sets of TCGA exhibited worse prognosis as the risk score increased (Figures 4(a)–4(c)). This conclusion was also supported in three external datasets samples (Supplementary Figures 6A–6C). Figures 4(d)–4(f) show that our model achieved an area under the curve (AUC) value of >0.85 in three sets of TCGA data. Surprisingly, such excellent performance can still stand in the test of external datasets (all AUC values > 0.65, Supplementary Figures 6D–6F). Kaplan–Meier survival curves show significant differences in survival probability between different risk groups. This conclusion is applicable to the corresponding tests based on all internal and external datasets (Figures 4(g)–4(i) and Supplementary Figures 6G–6I). Both univariate and multivariate analyses based on 6 sets all confirmed that risk score was a prognostic factor for patients with glioma (Figures 4(j)–4(o) and Supplementary Figures 6J–6O). Together, the above results support the excellent performance of our model in predicting the prognosis of glioma patients.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

(n)

(o)

Figure 4

Model performance testing based on training dataset (a), (d), (g), (j), and (m), test dataset (b), (e), (h), (k), and (n), whole set (c), (f), (i), (l), and (o). (a)–(c) Changes in sample’s survival with increasing risk score. (d)–(f) ROC curves reflecting model’s performance. (g)–(i) Kaplan–Meier survival curves reflecting model’s prognostic discrimination performance. (j)–(l) Results of univariate Cox regression analysis to determine influencing factors of prognosis. (m)–(o) Results of multivariate Cox regression analysis to determine independent influencing factors of prognosis.

The effect of the number of decision trees on the error rate from three datasets is shown in Supplementary Figures 7A–7C. From the 3 datasets, we observed that the importance scores of the 8 ERSGs were relatively high (all importance scores > 0.5), which suggested that these genes played relatively important roles in the model (Supplementary Figures 7D–7F). The diagnostic ROC curves of the genes used for modeling all verified their high diagnostic value in glioma (all AUC > 0.85, Supplementary Figures 7G–7M). These results support the stability of the model constructed by these eight genes.

3.6. Deep Validation of Model Performance and Stability

The heatmap visually shows 3 clinical features matched to the risk score of each sample from TCGA (Figure 5(a)). Additionally, the following Wilcoxon signed-rank test further provided many exciting results with clinical value. The risk scores between survival status, age, and grade subtypes show dramatical differences (Figures 5(b), 5(c), and 5(e)), demonstrating the close association between survival status/age/grade and risk score. It is worth pointing out that patients in the dead group, age ༞60 group, and G3 group were observed to have a higher risk score. However, gender characteristics did not show a similar conclusion (Figure 5(d)). More importantly, there were also significant differences in survival probability between samples from different risk groups within each clinical feature subgroup (Figures 5(f)–5(k)), indicating the guiding significance of risk score for prognosis.

Figure 5

Deep validation of model’s performance and stability based on TCGA dataset. (a) Heatmap showing clinical features matched to each sample’s risk score. . (b)–(e) Risk score differences between subgroups with different clinical features. (f)–(k) Kaplan–Meier survival analysis reflecting model’s prognostic discrimination performance in each clinical feature subgroup. (l)–(p) IHC staining images of F2R, HEY2, PXDN, RNF19A, and SCG5 in normal brain and glioma tissues, respectively.

We only obtained IHC staining images reflecting the protein expression of five PR-DE-ERSGs (F2R, HEY2, PXDN, RNF19A and SCG5) (Figures 5(l)–5(p)). As shown in the stained images, the four PR-DE-ERSGs (F2R, HEY2, PXDN, and RNF19A) exhibited differences in protein expression levels between glioma and normal brain tissues, consistent with their differential expression results based on the TCGA dataset (Figures 5(l)–5(o)). Unfortunately, the image of SCG5 fails to support the analysis results based on the TCGA dataset (Figure 5(p)).

Figures 6(a) and 6(b)showed the significantly higher DRAM1 and lower SCG5 relative mRNA expression in gliomas. But we still failed to verify the significant differential relative mRNA expression of BMP2 and NUP107 between gliomas and adjacent normal tissues (Figures 6(c) and 6(d)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

3.7. Enrichment Analysis for Biological Functions and Pathways Related to Different Risk Groups

To explore the relationship between different risk groups and biological pathways, GSEA was utilized to catalogue biological pathways modulated by risk groups. In the high-risk group, most of the GO analysis results and KEGG analysis results were immune-related functions (positive regulation of T cell proliferation, regulation of T cell mediated immunity, response to interferon-beta, T cell activation, T helper 2 cell differentiation, toll-like receptor signaling pathway) and pathways (antigen processing and presentation, JAK-STAT signaling pathway, natural killer cell mediated cytotoxicity, p53 signaling pathway) (Figures 6(e) and 6(g)). Likewise, the biological functions in the low-risk group are shown in Figure 6(f) (including cell differentiation in hindbrain, regulation of postsynaptic membrane potential, neurotransmitter receptor complex, and postsynaptic density membrane).

3.8. Exploring the Potential Role of Immune Infiltrates Types and Cell Stemness

Both correlation and difference analysis results supported a positive association between the overall immune/stromal cell score and the risk score (Figures 7(a)–7(d)). The distribution of 16 immune cells and 13 immune function scores as risk scores increased preliminarily showed their positive association with risk score (Figure 7(e)). Figures 7(f) and 7(g) further show the specific association of each immune cell/function score with the risk score. Notably, regulatory T (Treg) cells, CD8+ T cells, immature dendritic cells (iDCs), T helper 2 (Th2) cells, and macrophages were all positively associated with risk scores, and their higher scores in the high-risk group were also observed.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 7

Exploratory results related to tumor immunosuppressive environment. (a), (c), and (f) Correlation analysis results between risk score and overall immune cell/overall stromal cell/16 immune cells and 13 immune functions scores. (b), (d), and (g) Overall immune cell/overall stromal cell/16 immune cells and 13 immune functions scores difference between different risk groups. (e) The distribution of 16 immune cells and 13 immune functions scores as risk scores increased. ns: no significance; ; ; .

3.9. Exploring the Potential Role of Tumor Immunosuppressive Environment

The heatmap showed significant differences in the number of samples of C3, C4, and C5 immune infiltrating subtypes between different risk groups (Figure 8(a)). Higher risk scores in the C3/C4 subgroup and lower risk scores in the C5 subgroup further confirmed this distribution (Figure 8(b)). The above results indicated that the C5 subtype was dominant in the low-risk group, and the C3 and C4 subtypes were dominant in the high-risk group.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figures 8(c) and 8(d) show that mDNAsi is positively correlated with the risk score. We also explored the relationship between mDNAsi and four clinical features (survival state, age, gender, and grade). We found that mDNAsi had significantly high scores in the dead (Figure 8(e)), age > 60 (Figure 8(f)) and G3 grade (Figure 8(h)) populations. However, no significant association between gender and mDNAsi was found (Figure 8(g)). The above results imply the role of cell stemness in promoting glioma’s progression.

3.10. Exploring the Potential Role of HEY2-Related ceRNA Regulatory Network

Of all predicted miRNAs, only 25 miRNAs potentially regulating HEY2 exceeded three times in Starbase database prediction procedures. The corresponding coexpression network is shown in Figure 9(a). The mechanism whereby miRNAs regulate the expression of target genes elucidates that the potential miRNAs should be inversely correlated with HEY2 expression [47]. The negative correlation between miR-181b-5p/miR-369-3p and HEY2 expression level and the significantly lower expression of these two genes in glioma both support their negative regulation of HEY2 expression (Figures 9(b)–9(e)). The lower survival probability observed in the high expression group of these two miRNAs further confirms their essential roles in glioma progression (Figures 9(f) and 9(g)). The above analysis results well supported the roles of miR-181a-5p and miR-369-3p as upstream candidate miRNAs in regulating HEY2 expression.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

(n)

(o)

(p)

Figure 9

Exploratory results of HEY2-related ceRNA regulatory network. (a) The potential network of 25 miRNAs regulating HEY2 expression. (b) and (c) Correlation analysis result of HEY2 with miR-181a-5p/miR-369-3p, respectively. (d) and (e) Difference analysis result of the miR-181a-5p/miR-369-3p expression between normal and tumor tissue, respectively. (f) and (g) The effects of miR-181a-5p/miR-369-3p on prognosis of glioma patients, respectively. (h) and (i) Correlation analysis results of miR-369-3p with GAS5/LINC00689, respectively. (j) and (k) Correlation analysis results of HEY2 with GAS5/LINC00689, respectively. (l) and (m) Difference analysis result of the GAS5/LINC00689 expression between normal and tumor tissue, respectively. (n) and (o) The effects of GAS5/LINC00689 on prognosis of glioma patients, respectively. (p) A ceRNA regulatory network composed of 5 candidate lncRNAs, miR-369-3p and HEY2 (yellow box: correlation coefficient < −0.35, green box: correlation coefficient < −0.2).

The prediction procedures predicted 122/54 lncRNAs might regulate the expression of miR-181a-5p/miR-369-3p, respectively. Stringent filtering conditions ultimately provided only two final candidate lncRNAs for miR-369-3p. Eventually, the entire coexpression network is visualized in Figure 9(p). Corresponding correlation and differential analysis results further confirmed the inverse regulatory relationship between these two lncRNAs (LINC00689 and GAS5) and miR-369-3p in glioma (Figures 9(h)–9(m)). These two lncRNAs also showed a positive effect on the prognosis of glioma patients (Figures 9(n) and 9(o)).

3.11. Exploring the Potential Role of Mutation

The waterfall charts provide an overview of the mutation distribution of the top 20 genes with mutation frequency in the samples of different risk groups (Figures 10(a) and 10(b)). TP53 and IDH1 are the two genes with the highest mutation frequency in those two charts. In addition, 17 of 24 DE-ERSGs strongly related to glioma were mutated (Figure 10(c)). Correlation and difference analysis both supported the close association between TMB and the model (Figures 10(d) and 10(e)). Clearly, TMB has a negative prognostic effect on gliomas (Figure 10(f)). Figures 10(i)–10(k) and 10(m) reveal a strong association between TP53/IDH1 mutation and 8 PR-DE-ERSGs/risk score, suggesting cooperative roles of TP53/IDH1 mutation and ERS in glioma. Figures 10(l) and 10(n) show the positive effect of TP53/IDH1 mutation on the prognosis of glioma patients. Lower CD274 expression in the TP53/IDH1 mutation group revealed the predictive value of TP53/IDH1 mutation in response to immunotherapy (Figures 10(g) and 10(h)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

(n)

Figure 10

Exploratory results of mutation. (a) and (b) An overview of the mutation distribution in the samples of different risk groups of the top 20 genes with mutation frequency, respectively. (c) A mutation distribution overview of 17 mutated DE-ERSGs strongly related to glioma. (d) Correlation analysis result between risk score and TMB. (e) Difference analysis result of TMB between different risk groups. (f) Kaplan–Meier survival analysis showed the effect of TMB on prognosis. (g) and (i) CD274/8 PR-DE-ERSGs expression difference analysis between IDH1 mutant and wild groups. (h) and (j) CD274/8 PR-DE-ERSGs expression difference analysis between TP53 mutant and TP53 wild groups. (k) and (l) Risk score and survival difference between IDH1 mutant and wild groups. (m) and (n) Risk score and survival difference between TP53 mutant and TP53 wild groups.

3.12. Evaluation of the Guiding Value of the Model in Therapy

Figures 11(a)–11(h) show the corresponding results for the top six most strongly correlated with 8 PR-DE-ERSGs expression among 263 FDA approved or clinically tested drugs. It is worth mentioning that we observed an inverse correlation between DRAM1 expression and carmustine (a chemotherapeutic drug recommended by the latest National Comprehensive Cancer Network (NCCN) guidelines for the treatment of glioma) (Figure 11(a)). This suggests that the 8 PR-DE-ERSGs in the model may serve as guide markers for the rational use of these drugs. Combining the correlation results, we found that the expressions of PXDN, DRAM1, and BMP2 were all closely correlated with IPSs and TIDE scores, implying the guiding value of these genes in immunotherapy strategies (Figures 12(a) and 12(b)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(a)

(b)

3.13. Extending the Model’s Clinical Value

To expand the clinical utility scope of the model, we developed a novel nomogram to facilitate the prediction of the 1-, 2-, and 3-year survival probability of glioma patients (Figure 13(a)). Our nomogram considers the combined effects of composite factors (risk group, age, and grade) on prognosis. The calibration curves demonstrate the excellent internal consistency between the predictions made by our nomogram and the actual outcomes (Figures 13(b)–13(d)). Also, ROC curves show high AUC values in three sets of TCGA (Figures 13(e)–13(g)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 13

Construction of a nomogram that effectively extends the clinical value of the model. (a) A novel nomogram to accurately predict the 1-, 2- and 3-year survival probability of glioma patients. (b)–(d) Internal calibration curves of 1, 2 and 3 years based on three sets, respectively. Higher degree of fitting with the grey line represents higher accuracy of the Nomogram prediction. (e)–(g) ROC curves of 1, 2 and 3 years based on three sets, respectively. AUC > 0.5 is considered to have a predictive value.

Interestingly, with the exception of VTCN1, TMIGD2, CD200, HHLA2, and ADORA2A, most ICIs related gene expression showed positively correlated with the risk score (Supplementary Figure 8A). In addition, half of the m6A-related genes demonstrated a positive correlation with the risk score (Supplementary Figures 8A and 8C). The following difference analysis results shown in Supplementary Figures 8B and 8C further support the above conclusions.

After analyzing the correlation between multidrug resistance-related genes and the prognosis model, we found that ABCC1 and ABCC3 gene expressions were significantly positively correlated with the risk score (Supplementary Figures 8D and 8F). And both ABCC1 and ABCC3 were highly expressed in the high-risk group (Supplementary Figures 8E and 8G).

4. Discussion

WGCNA, a systems biology method for describing the correlation patterns among genes across microarray samples, can be used for finding clusters (modules) of highly correlated genes and relating modules to one another and to external sample traits [48]. WGCNA is widely used to identify tumor candidate biomarkers and explore the mechanism, such as gastric cancer [49], bladder cancer [50], gastric adenocarcinoma [51], and breast cancer [52]. This study utilized WGCNA and univariate Cox regression to screen out 8 PR-DE-ERSGs highly correlated with glioma to establish a prognostic model. We also established a nomogram by combining the prognosis model with clinicopathological factors to predict survival in glioma patients. In repeated validation tests, both the prognostic model and nomogram demonstrated excellent prognostic performance and substantial clinical value. The strong association between risk score and immune cells/functions, cell stemness, and IDH1/TP53 mutation further elucidates the important role of ERS in glioma progression and treatment. The TCGA data further confirmed the important regulatory value of the predicted HEY2/miR-369-3p/lncRNAs networks in glioma.

To better understand the potential molecular mechanism of our prognostic model, we have expanded on these genes. As an essential component of the nuclear pore complex, Nucleoporin 107 (NUP107) may contribute to the regulation of cell fate in aged and transformed cells by modulating nuclear trafficking of signal based on undifferentiated oligodendroglioma cells in vitro [53]. Damage Regulated Autophagy Modulator 1 (DRAM1), a p53 target gene encoding a lysosomal protein that induces macroautophagy, plays a central role in p53/TP53-mediated apoptosis [54]. Galavotti et al. found that high levels of DRAM1 were associated with shorter overall survival in GBM patients, which is consistent with our findings [55]. Ring Finger Protein 19A (RNF19A) encodes a protein that belongs to the ring fingers protein family, an E3 ubiquitin ligase [56]. RNF19A decreases p53 expression and its downstream signaling, binds to p53, and promotes its ubiquitination, thereby promoting nonsmall cell lung carcinoma (NSCLC) growth and progression [57]. Peroxidasin (PXDN) encodes a heme-containing peroxidase that is involved in extracellular matrix formation associated with prostate cancer [58]. Overexpression of PXDN is found in ovarian cancer [59], melanoma [60], oral squamous cell carcinoma [61], and prostate cancer [58], associated with poor prognosis. All aforementioned studies support our findings that PXDN is a prognostic risk factor in glioma. Coagulation factor 2 thrombin receptor (F2R), a key component in the thrombosis process, has been demonstrated up-regulated in various cancers and was associated with several protumoral responses, including primary growth, invasion, metastasis, and angiogenesis [62–66]. In glioma, the expression of F2R was up-regulated and predicted poor prognosis. F2R promotes glioma cell proliferation and metastasis under SOX2 and actives the WNT/β-catenin signaling pathway [62]. The epigenetic inactivation of secretogranin V (SCG5/SGNE1/7B2) was a frequent early event in glioma formation, resulting in significant downregulation of SCG5/SGNE1/7B2 expression [67]. Bone morphogenetic protein 2 (BMP2), a member of the transforming growth factor-β (TGF-β) super-family, is one of the main chondrogenic growth factors involved in cartilage regeneration [68]. In the new glioma grading model constructed by Zhou et al., BMP2 shows reliable prognostic value as a protective factor in GBM, consistent with our results [69]. In our study, the results of human tissue-based IHC staining and QRT-PCR experiments well confirmed the biological value of F2R, HEY2, PXDN, RNF19A, SCG5, and DRAM1 in glioma.

Recently, the close functional relationship between ERS and noncoding RNAs (ncRNAs), including miRNAs, lncRNAs, and circRNAs, has been widely reported in cancer development [70]. MiRNAs and lncRNAs directly or indirectly act on UPR pathway molecules to regulate intracellular homeostasis and affect carcinogenic processes, including survival, apoptosis, invasion, metastasis, cancer stem cell characteristics, and the tumor microenvironment [70]. However, few studies focus on ERS-related ceRNA regulatory network in gliomas. After prediction and verification of TCGA data, a ceRNA regulatory network composed of HEY2, miR-369-3p, and GAS5/LINC00689 was screened out.

HEY2, one of the most prominent Notch pathway target genes, encodes a transcription factor [71]. Halani et al. proposed that the loss of the Notch pathway activity and particularly of Hey2 levels were correlated with oligodendroglioma [72]. Moreover, Giachino et al. found that high HEY2 expression correlates with a better prognosis for patients with grades II–III astrocytoma and GBM, further supporting our observation [73]. Previous studies have shown that overexpression of growth arrest-specific transcript 5 (GAS5) inhibits the malignant phenotype of glioma cells, including proliferation, migration, and invasion, acting as a tumor suppressor of human gliomas [74]. GAS5 may blunt the resistance of glioma cells to cisplatin by suppressing excessive autophagy through the activation of mTOR signaling, implying a promising therapeutic strategy against chemoresistance in glioma [75]. However, increased LINC00689 level was associated with poor clinical features and decreased overall survival in glioma patients, and the overexpression of miR-369-3p inhibited the proliferation and migration in glioblastoma cells [76, 77]. Those appear to be different from ours. However, the predominant roles of HEY2, miR-181a-5p, and GAS5 and the cross-talk between multiple ceRNAs may explain our results. Interestingly, miR-369-3p, miR-181a-5p, and LINC00689 were found to be related to glycolysis [77–79]. And the high rate of aerobic glycolysis is a key characteristic of cancer cells, which promotes uncontrolled proliferation [80]. Recently, accumulating studies have illustrated that activated aerobic glycolysis participated in various cellular and clinical activities of glioma, thus influencing the efficacy of radiotherapy and chemotherapy [80]. Taken together, miRNA dysregulation in glioma disables proper tumor suppression, increases glycolytic metabolism, and augments tumor malignancy through multiple effectors and signaling pathways [81].

The majority of immune cell and function scores in the high-risk group further explain the enrichment of numerous immune-related biological processes in the high-risk group. Infiltration of Treg cells is often enriched in tumor tissue, and a high proportion of Treg cells and effector T cells often predicts poorer patient outcomes [82]. It has been suggested that Treg cells in tumor-associated tertiary lymphoid structures can suppress endogenous immune responses against tumors, resulting in poor prognosis [83]. More and more studies have shown that tumor-associated myeloid cells, such as iDCs, can be recruited in the tumor environment and form an organic whole to maintain an immunosuppressive state [84]. The Siglec-9 expressed on this cell can bind to the Mucins1-related STn antigen released by cancer cells and lead to an increase in interleukin-10 and a decrease in the production of IL-12 [85]. This prevents the DC’s ability to activate T-helper 1 responses [86]. According to the tissue environment, macrophages can be divided into M1 macrophages with proinflammatory and tumor-suppressive effects and M2 macrophages with immunosuppressive and tumor-promoting effects [87]. While tumor-associated macrophages generally acquire the same phenotype as M2 macrophages, their roles in cancer are often associated with poor prognosis [88]. Kakizaki et al. found that the cells, together with Tregs, play a role in maintaining the tumor microenvironment [89]. In conclusion, the reason for the poor prognosis of the high-risk group but the enrichment of more immune-related processes may be related to the fact that these immune cells can maintain the tumor immunosuppressive environment.

The waterfall plot showed that the two most frequently mutated genes were isocitrate dehydrogenases 1 (IDH1) and Tumor Protein p53 (TP53) in the high-risk and low-risk groups. IDH1 mutation is an early event in glioma development and positively correlates with TP53 mutation [90, 91]. And the co-occurrence of IDH1 mutations and TP53 alterations is widespread in gliomas [92]. An increasing number of reports have shown that gliomas patients with IDH1 mutations have a better prognosis than the one with wild-type IDH1 [93, 94]. IDH1 is increasingly recognized as an independent prognostic marker for gliomas, consistent with our findings [95]. Importantly, Amankulor et al. suggested that differences in immune cell content may contribute to this prognostic difference [96]. By creating mouse models, they found reduced levels of many types of immune cells in IDH1-mutated tumors, such as macrophages, microglia, monocytes, and neutrophils [96]. Also, previous reports suggested that significant infiltration of these immune cells leads to poor prognosis in a variety of cancers [88, 97–99]. Surprisingly, immune cells such as macrophages were also less in the patients of the low-risk group with a better prognosis. Based on the mutual support between these studies and our findings, we speculate that the favorable prognosis of the IDH1 mutant group may also be related to the reduction of immune cells capable of maintaining the tumor immunosuppressive environment.

As the results show, our model exhibited excellent capabilities in predicting the prognosis of patients with glioma and facilitating the optimal selection of clinical treatment strategies. However, the research still has numerous limitations. For example, the difference in protein expression of SCG5 is still not supported by IHC staining images obtained based on experiments. In addition, our research failed to combine basic experiments to verify the study results. Considering the complexity of the ceRNA network, we only selected HEY2 to explore the ceRNA axis to avoid involving a large number of RNA transcription samples that make the detailed study laborious. Despite many shortcomings, our study still provides a rich exploration and in-depth explanation of the underlying mechanism. These may provide new ideas for follow-up research.

5. Conclusions

In this study, 8 PR-DE-ERSGs strongly related to glioma were screened out by WGCNA to establish a prognostic model with excellent prognostic performance and important clinical value. The strong association between risk score and immune cells/functions, cell stemness, and IDH1/TP53 mutation further elucidates the important role of ERS in glioma progression and treatment. The TCGA data further confirmed the important regulatory value of the predicted HEY2/miR-369-3p/lncRNAs networks in glioma. All in all, the above results all supported that biological processes related to HEY2-based ceRNA regulatory network, tumor immunosuppressive environment, cell stemness, and IDH1 mutation may be involved in the progression and treatment of glioma.

Data Availability

The datasets analyzed in this study came from databases shared publicly. Data can be obtained from CellMiner https://discover.nci.nih.gov/cellminer, CGGA http://www.cgga.org.cn/, GeneCards http://www.genecards.org/, GEO https://www.ncbi.nlm.nih.gov/geo/, HPA http://www.proteinatlas.org/, mirbase https://www.mirbase.org/, StarBase http://starbase.sysu.edu.cn/, TCGA https://cancergenome.nih.gov/, TCIA https://tcia.at/home, TIDE http://tide.dfci.harvard.edu/.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

The authors sincerely thank all the participants and the supporting from the Research Fund Project of Jiangxi Provincial Department of Education (180017).

Supplementary Materials

Supplementary Figure 1. Differential expression analysis of ERSGs in the TCGA dataset. (A) Heap map of the differential expression of ERSGs in the TCGA dataset. (B) Volcano plot of DE-ERSGs in the TCGA dataset. (red dots: up-regulated genes, black dots: genes with no significant differences, green dots: down-regulated genes). Supplementary Figure 2. Differential expression analysis of ERSGs in the GSE108474 dataset. (A) Heap map of the differential expression of ERSGs in the GSE108474 set. (B) Volcano plot of DE-ERSGs in the GSE108474 set. (red dots: up-regulated genes, black dots: genes with no significant differences, green dots: down-regulated genes). Supplementary Figure 3. Differential expression analysis of ERSGs in the GSE4290 dataset. (A) Heap map of the differential expression of ERSGs in the GSE4290 set. (B) Volcano plot of DE-ERSGs in the GSE4290 dataset. (red dots: up-regulated genes, black dots: genes with no significant differences, green dots: down-regulated genes). Supplementary Figure 4. LASSO analysis based on the training set. (A) The minimum standard. (B) Coefficients. Supplementary Figure 5. Validation of differential expression of 8 genes by GEPIA database’s data. Supplementary Figure 6. The testing results of model’s performance based on GSE4412, GSE43378, and CGGA datasets. (A–C) The changes in sample’s survival state with increasing risk score. (D–F) ROC curves reflecting model’s prognostic predictive performance. (G–I) Kaplan–Meier survival curves reflecting model’s prognostic discrimination performance. (J–L) The results of univariate Cox regression analysis to determine influencing factors of prognosis. (M–O) The results of multivariate Cox regression analysis to determine independent influencing factors of prognosis. Supplementary Figure 7. Calculate the relative importance of 8 PR-DE-ERSGs in glioma based on RF algorithm. (A–C) The error rates are influenced by the number of decision trees in three datasets, respectively. The x-axis represents the number of decision trees and the y-axis is the error rate. (D–F) The importance score of the 8 PR-DE-ERSGs in three datasets, respectively. The x-axis represents the importance index, and the y-axis represents the genes. (G–N) The ROC curves of 8 PR-DE-ERSGs from TCGA dataset. Supplementary Figure 8. Extending the model’s clinical value. (A) Correlation analysis of risk scores with ICIs and m6A-related genes expression. (B–C) Differential expression analysis of ICIs and m6A-related genes in different risk-score subtypes, respectively. (D, F) Correlation analysis result between risk-score and ABCC1/ABCC3, respectively. (E, G) Differences in ABCC1/ABCC3 expression between different risk groups, respectively. Supplementary Table 1: the DE-FRGs used to construct the model and the corresponding coefficients. Supplementary Table 2. All primer sequences for qRT-PCR. (Supplementary Materials)

References

Q. T. Ostrom, P. Liao, K. T. Vecchione, Y. Wolinsky, and C. Kruchko, “CBTRUS Statistical Report: primary brain and other central nervous system tumors diagnosed in the United States in 2010-2014,” Neuro-Oncology, vol. 19, no. 5, 2017.
View at: Publisher Site | Google Scholar
Q. T. Ostrom, L. Bauchet, F. G. Davis et al., “The epidemiology of glioma in adults: a “state of the science” review,” Neuro-Oncology, vol. 16, no. 7, pp. 896–913, 2014.
View at: Publisher Site | Google Scholar
J. Chen, R. M. McKay, and L. F. Parada, “Malignant glioma: lessons from genomics, mouse models, and stem cells,” Cell, vol. 149, no. 1, pp. 36–47, 2012.
View at: Publisher Site | Google Scholar
S. Xu, L. Tang, X. Li, F. Fan, and Z. Liu, “Immunotherapy for glioma: current management and future application,” Cancer Letters, vol. 476, pp. 1–12, 2020.
View at: Publisher Site | Google Scholar
L. Lin, J. Cai, and C. Jiang, “Recent Advances in targeted therapy for glioma,” Current Medicinal Chemistry, vol. 24, no. 13, pp. 1365–1381, 2017.
View at: Publisher Site | Google Scholar
A. C. Tan, D. M. Ashley, G. Y. López, M. Malinzak, H. S. Friedman, and M. Khasraw, “Management of glioblastoma: state of the art and future directions,” CA: A Cancer Journal for Clinicians, vol. 70, no. 4, pp. 299–312, 2020.
View at: Publisher Site | Google Scholar
J. G. Nicholson and H. A. Fine, “Diffuse glioma heterogeneity and its therapeutic implications,” Cancer Discovery, vol. 11, no. 3, pp. 575–590, 2021.
View at: Publisher Site | Google Scholar
V. Rajendran, R. Purohit, and R. Sethumadhavan, “In silico investigation of molecular mechanism of laminopathy caused by a point mutation (R482W) in lamin A/C protein,” Amino Acids, vol. 43, no. 2, pp. 603–615, 2012.
View at: Publisher Site | Google Scholar
V. Rajendran and R. Sethumadhavan, “Drug resistance mechanism of PncA in Mycobacterium tuberculosis,” Journal of Biomolecular Structure and Dynamics, vol. 32, no. 2, pp. 209–221, 2014.
View at: Publisher Site | Google Scholar
S. Kumar, V. K. Bhardwaj, R. Singh, P. Das, and R. Purohit, “Identification of acridinedione scaffolds as potential inhibitor of DENV-2 C protein: an in silico strategy to combat dengue,” Journal of Cellular Biochemistry, vol. 123, no. 5, pp. 935–946, May 2022.
View at: Publisher Site | Google Scholar
S. A. Oakes and F. R. Papa, “The role of endoplasmic reticulum stress in human pathology,” Annual Review of Pathology: Mechanisms of Disease, vol. 10, no. 1, pp. 173–194, 2015.
View at: Publisher Site | Google Scholar
M. Wang and R. J. Kaufman, “The impact of the endoplasmic reticulum protein-folding environment on cancer development,” Nature Reviews Cancer, vol. 14, no. 9, pp. 581–597, 2014.
View at: Publisher Site | Google Scholar
X. Chen and J. R. Cubillos-Ruiz, “Endoplasmic reticulum stress signals in the tumour and its microenvironment,” Nature Reviews Cancer, vol. 21, no. 2, pp. 71–88, 2021.
View at: Publisher Site | Google Scholar
J. R. Cubillos-Ruiz, S. E. Bettigole, and L. H. Glimcher, “Tumorigenic and immunosuppressive effects of endoplasmic reticulum stress in cancer,” Cell, vol. 168, no. 4, pp. 692–706, Feb. 2017.
View at: Publisher Site | Google Scholar
Y. He, J. Su, B. Lan, Y. Gao, and J. Zhao, “Targeting off-target effects: endoplasmic reticulum stress and autophagy as effective strategies to enhance temozolomide treatment,” OncoTargets and Therapy, vol. 12, pp. 1857–1865, 2019.
View at: Publisher Site | Google Scholar
N. M. Peñaranda-Fajardo, C. Meijer, Y. Liang et al., “ER stress and UPR activation in glioblastoma: identification of a noncanonical PERK mechanism regulating GBM stem cells through SOX2 modulation,” Cell Death & Disease, vol. 10, no. 10, p. 690, 2019.
View at: Publisher Site | Google Scholar
M. Cattaneo, S. Baronchelli, D. Schiffer et al., “Down-modulation of SEL1L, an unfolded protein response and endoplasmic reticulum-associated degradation protein, sensitizes glioma stem cells to the cytotoxic effect of valproic acid,” Journal of Biological Chemistry, vol. 289, no. 5, pp. 2826–2838, 2014.
View at: Publisher Site | Google Scholar
C. Giampietri, S. Petrungaro, S. Conti, A. Facchiano, A. Filippini, and E. Ziparo, “Cancer microenvironment and endoplasmic reticulum stress response,” Mediators of Inflammation, vol. 2015, Article ID 417281, 11 pages, 2015.
View at: Publisher Site | Google Scholar
M. Markouli, D. Strepkos, A. G. Papavassiliou, and C. Piperi, “Targeting of endoplasmic reticulum (ER) stress in gliomas,” Pharmacological Research, vol. 157, Article ID 104823, 2020.
View at: Publisher Site | Google Scholar
F. A. Karreth and P. P. Pandolfi, “ceRNA cross-talk in cancer: when ce-bling rivalries go awry,” Cancer Discovery, vol. 3, no. 10, pp. 1113–1121, 2013.
View at: Publisher Site | Google Scholar
L. Salmena, L. Poliseno, Y. Tay, L. Kats, and P. P. Pandolfi, “A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language?” Cell, vol. 146, no. 3, pp. 353–358, Aug. 2011.
View at: Publisher Site | Google Scholar
L. J. Li, W. Zhao, S. S. Tao et al., “Competitive endogenous RNA network: potential implication for systemic lupus erythematosus,” Expert Opinion on Therapeutic Targets, vol. 21, no. 6, pp. 639–648, 2017.
View at: Publisher Site | Google Scholar
Y. Wang, J. Hou, D. He et al., “The emerging function and mechanism of ceRNAs in cancer,” Trends in Genetics, vol. 32, no. 4, pp. 211–224, Apr. 2016.
View at: Publisher Site | Google Scholar
X.-Z. Yang, T. T. Cheng, Q. J. He et al., “LINC01133 as ceRNA inhibits gastric cancer progression by sponging miR-106a-3p to regulate APC expression and the Wnt/β-catenin pathway,” Molecular Cancer, vol. 17, no. 1, p. 126, 2018.
View at: Publisher Site | Google Scholar
J. Wang, X. Zhao, Y. Wang et al., “circRNA-002178 act as a ceRNA to promote PDL1/PD1 expression in lung adenocarcinoma,” Cell Death & Disease, vol. 11, no. 1, p. 32, 2020.
View at: Publisher Site | Google Scholar
Q. Han, J. Li, J. Xiong, and Z. Song, “Long noncoding RNA LINC00514 accelerates pancreatic cancer progression by acting as a ceRNA of miR-28-5p to upregulate Rap1b expression,” Journal of Experimental & Clinical Cancer Research, vol. 39, no. 1, p. 151, 2020.
View at: Publisher Site | Google Scholar
X.-S. Wu, F. Wang, H. Li et al., “LncRNA-PAGBC acts as a microRNA sponge and promotes gallbladder tumorigenesis,” EMBO Reports, vol. 18, no. 10, pp. 1837–1853, 2017.
View at: Publisher Site | Google Scholar
G. Stelzer, N. Rosen, I. Plaschkes et al., “The GeneCards suite: from gene data mining to disease genome sequence analyses,” Current Protocols in Bioinformatics, vol. 54, p. 1, 2016.
View at: Publisher Site | Google Scholar
C. Y. Li, J. H. Cai, J. J. P. Tsai, and C. C. N. Wang, “Identification of Hub genes associated with development of Head and Neck squamous cell carcinoma by integrated Bioinformatics analysis,” Frontiers in Oncology, vol. 10, p. 681, 2020.
View at: Publisher Site | Google Scholar
R. Tibshirani, “The lasso method for variable selection in the Cox model,” Statistics in Medicine, vol. 16, no. 4, pp. 385–395, 1997.
View at: Publisher Site | Google Scholar
Z. Tang, C. Li, B. Kang, G. Gao, C. Li, and Z. Zhang, “GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses,” Nucleic Acids Research, vol. 45, no. W1, pp. W98–W102, 2017.
View at: Publisher Site | Google Scholar
U. Sunarya, Y. Sun Hariyani, T. Cho et al., “Feature analysis of Smart Shoe sensors for classification of Gait patterns,” Sensors, vol. 20, pp. 6253–21, 2020.
View at: Publisher Site | Google Scholar
E. A. Clayton, T. A. Pujol, J. F. McDonald, and P. Qiu, “Leveraging TCGA gene expression data to build predictive models for cancer drug response,” BMC Bioinformatics, vol. 21, no. S14, p. 364, 2020.
View at: Publisher Site | Google Scholar
M. Lu, J. Li, X. Fan, F. Xie, J. Fan, and Y. Xiong, “Novel immune-related Ferroptosis signature in esophageal cancer: an Informatics exploration of biological processes related to the TMEM161B-AS1/hsa-miR-27a-3p/GCH1 regulatory network,” Frontiers in Genetics, vol. 13, Article ID 829384, 2022.
View at: Publisher Site | Google Scholar
Z. Lin, Q. Xu, D. Miao, and F. Yu, “An inflammatory response-related gene signature can impact the immune status and predict the prognosis of Hepatocellular carcinoma,” Frontiers in Oncology, vol. 11, Article ID 644416, 2021.
View at: Publisher Site | Google Scholar
K. Yoshihara, M. Shahmoradgoli, E. Martinez et al., “Inferring tumour purity and stromal and immune cell admixture from expression data,” Nature Communications, vol. 4, no. 1, p. 2612, 2013.
View at: Publisher Site | Google Scholar
D. Tamborero, C. Rubio-Perez, F. Muinos et al., “A Pan-cancer Landscape of interactions between solid tumors and infiltrating immune cell populations,” Clinical Cancer Research, vol. 24, no. 15, pp. 3717–3728, 2018.
View at: Publisher Site | Google Scholar
T. M. Malta, A. Sokolov, A. J. Gentles et al., “Machine learning Identifies stemness features associated with Oncogenic Dedifferentiation,” Cell, vol. 173, no. 2, pp. 338–354.e15, 2018.
View at: Publisher Site | Google Scholar
M. Lim, Y. Xia, C. Bettegowda, and M. Weller, “Current state of immunotherapy for glioblastoma,” Nature Reviews Clinical Oncology, vol. 15, no. 7, pp. 422–442, 2018.
View at: Publisher Site | Google Scholar
W. Roth and M. Weller, “Chemotherapy and immunotherapy of malignant glioma: molecular mechanisms and clinical perspectives,” Cellular and Molecular Life Sciences, vol. 56, no. 5–6, pp. 481–506, Oct. 1999.
View at: Publisher Site | Google Scholar
P. Charoentong, F. Finotello, M. Angelova et al., “Pan-cancer Immunogenomic analyses reveal Genotype-Immunophenotype relationships and predictors of response to checkpoint Blockade,” Cell Reports, vol. 18, no. 1, pp. 248–262, Jan. 2017.
View at: Publisher Site | Google Scholar
P. Jiang, S. Gu, D. Pan et al., “Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response,” Nature Medicine, vol. 24, no. 10, pp. 1550–1558, 2018.
View at: Publisher Site | Google Scholar
Y. Chen, Z. Y. Li, G. Q. Zhou, and Y. Sun, “An immune-related gene prognostic index for Head and Neck squamous cell carcinoma,” Clinical Cancer Research, vol. 27, no. 1, pp. 330–341, Jan. 2021.
View at: Publisher Site | Google Scholar
Y. Qi, B. Liu, Q. Sun, X. Xiong, and Q. Chen, “Immune checkpoint targeted therapy in glioma: status and Hopes,” Frontiers in Immunology, vol. 11, Article ID 578877, 2020.
View at: Publisher Site | Google Scholar
Y. Zhang, X. Geng, Q. Li et al., “m6A modification in RNA: biogenesis, functions and roles in gliomas,” Journal of Experimental & Clinical Cancer Research, vol. 39, no. 1, p. 192, 2020.
View at: Publisher Site | Google Scholar
S. Hasegawa, K. Taniguchi, A. Yokomizo et al., “Possible involvement of multidrug-resistance-associated protein (MRP) gene expression in spontaneous drug resistance to vincristine, etoposide and adriamycin in human glioma cells,” International Journal of Cancer, vol. 58, no. 6, pp. 860–864, Sep. 1994.
View at: Publisher Site | Google Scholar
A. R. Gawronski, M. Uhl, Y. Zhang et al., “MechRNA: prediction of lncRNA mechanisms from RNA-RNA and RNA-protein interactions,” Bioinformatics, vol. 34, no. 18, pp. 3101–3110, 2018.
View at: Publisher Site | Google Scholar
P. Langfelder and S. Horvath, “WGCNA: an R package for weighted correlation network analysis,” BMC Bioinformatics, vol. 9, no. 1, p. 559, 2008.
View at: Publisher Site | Google Scholar
Z. H. Ren, G. P. Shang, K. Wu, C. Y. Hu, and T. Ji, “WGCNA Co-expression network analysis reveals ILF3-AS1 functions as a CeRNA to regulate PTBP1 expression by sponging miR-29a in gastric cancer,” Frontiers in Genetics, vol. 11, p. 39, 2020.
View at: Publisher Site | Google Scholar
Y. Wang, L. Chen, L. Ju et al., “Novel biomarkers associated with progression and prognosis of bladder cancer identified by Co-expression analysis,” Frontiers in Oncology, vol. 9, p. 1030, 2019.
View at: Publisher Site | Google Scholar
R. Tan, G. Zhang, R. Liu et al., “Identification of early diagnostic and prognostic biomarkers via WGCNA in Stomach adenocarcinoma,” Frontiers in Oncology, vol. 11, Article ID 636461, 2021.
View at: Publisher Site | Google Scholar
J. Tang, D. Kong, Q. Cui et al., “Prognostic genes of breast cancer identified by gene Co-expression network analysis,” Frontiers in Oncology, vol. 8, p. 374, 2018.
View at: Publisher Site | Google Scholar
S. Y. Kim, H. T. Kang, H. R. Choi, and S. C. Park, “Reduction of Nup107 attenuates the growth factor signaling in the senescent cells,” Biochemical and Biophysical Research Communications, vol. 401, no. 1, pp. 131–136, Oct. 2010.
View at: Publisher Site | Google Scholar
D. Crighton, S. Wilkinson, J. O’Prey et al., “DRAM, a p53-induced modulator of autophagy, is critical for apoptosis,” Cell, vol. 126, no. 1, pp. 121–134, 2006.
View at: Publisher Site | Google Scholar
S. Galavotti, S. Bartesaghi, D. Faccenda et al., “The autophagy-associated factors DRAM1 and p62 regulate cell migration and invasion in glioblastoma stem cells,” Oncogene, vol. 32, no. 6, pp. 699–712, 2013.
View at: Publisher Site | Google Scholar
J. i. Niwa, S. Ishigaki, M. Doyu, T. Suzuki, K. Tanaka, and G. Sobue, “A novel centrosomal ring-finger protein, dorfin, mediates ubiquitin ligase activity,” Biochemical and Biophysical Research Communications, vol. 281, no. 3, pp. 706–713, Mar. 2001.
View at: Publisher Site | Google Scholar
Y. Cheng, Y. Hu, H. Wang et al., “Ring finger protein 19A is overexpressed in non-small cell lung cancer and mediates p53 ubiquitin-degradation to promote cancer growth,” Journal of Cellular and Molecular Medicine, vol. 25, no. 16, pp. 7796–7808, 2021.
View at: Publisher Site | Google Scholar
J. Dougan, O. Hawsawi, L. J. Burton et al., “Proteomics-metabolomics combined Approach Identifies Peroxidasin as a protector against metabolic and oxidative stress in prostate cancer,” International Journal of Molecular Sciences, vol. 20, no. 12, p. 3046, 2019.
View at: Publisher Site | Google Scholar
Y.-Z. Zheng and L. Liang, “High expression of PXDN is associated with poor prognosis and promotes proliferation, invasion as well as migration in ovarian cancer,” Annals of Diagnostic Pathology, vol. 34, pp. 161–165, 2018.
View at: Publisher Site | Google Scholar
A. Jayachandran, P. Prithviraj, P. H. Lo et al., “Identifying and targeting determinants of melanoma cellular invasion,” Oncotarget, vol. 7, no. 27, pp. 41186–41202, 2016.
View at: Publisher Site | Google Scholar
M. Kurihara-Shimomura, T. Sasahira, H. Shimomura, and T. Kirita, “Peroxidan plays a tumor-promoting role in oral squamous cell carcinoma,” International Journal of Molecular Sciences, vol. 21, no. 15, p. 5416, 2020.
View at: Publisher Site | Google Scholar
G. Gao, M. Yang, F. Wang et al., “Coagulation factor 2 thrombin receptor promotes malignancy in glioma under SOX2 regulation,” Aging (Albany NY), vol. 12, no. 11, pp. 10594–10613, 2020.
View at: Publisher Site | Google Scholar
F. Schaffner and W. Ruf, “Tissue factor and PAR2 signaling in the tumor microenvironment,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 29, no. 12, p. 2004, 2009.
View at: Publisher Site | Google Scholar
M. Jaber, M. Maoz, A. Kancharla et al., “Protease-activated-receptor-2 affects protease-activated-receptor-1-driven breast cancer,” Cellular and Molecular Life Sciences, vol. 71, no. 13, pp. 2517–2533, 2014.
View at: Publisher Site | Google Scholar
A. Boire, L. Covic, A. Agarwal, S. Jacques, S. Sherifi, and A. Kuliopulos, “PAR1 is a matrix metalloprotease-1 receptor that promotes invasion and tumorigenesis of breast cancer cells,” Cell, vol. 120, no. 3, pp. 303–313, 2005.
View at: Publisher Site | Google Scholar
F. S. Ribeiro, T. A. Simao, N. D. Amoedo et al., “Evidence for increased expression of tissue factor and protease-activatedreceptor-1 in human esophageal cancer,” Oncology Reports, vol. 21, no. 6, pp. 1599–1604, 2009.
View at: Publisher Site | Google Scholar
A. Waha, J. Felsberg, W. Hartmann et al., “Frequent epigenetic inactivation of the chaperone SGNE1/7B2 in human gliomas,” International Journal of Cancer, vol. 131, no. 3, pp. 612–622, Aug. 2012.
View at: Publisher Site | Google Scholar
N. Zhou, Q. Li, X. Lin et al., “BMP2 induces chondrogenic differentiation, osteogenic differentiation and endochondral ossification in stem cells,” Cell and Tissue Research, vol. 366, no. 1, pp. 101–111, Oct. 2016.
View at: Publisher Site | Google Scholar
K. Zhou, Z. Zhao, S. Li, Y. Liu, G. Li, and T. Jiang, “A new glioma grading model based on histopathology and Bone Morphogenetic Protein 2 mRNA expression,” Scientific Reports, vol. 10, no. 1, Article ID 18420, 2020.
View at: Publisher Site | Google Scholar
T. Zhao, J. Du, and H. Zeng, “Interplay between endoplasmic reticulum stress and non-coding RNAs in cancer,” Journal of Hematology & Oncology, vol. 13, no. 1, p. 163, 2020.
View at: Publisher Site | Google Scholar
C. Barone, M. Buccarelli, F. Alessandrini et al., “Sox2-dependent maintenance of mouse oligodendroglioma involves the Sox2-mediated downregulation of Cdkn2b, Ebf1, Zfp423, and Hey2,” Glia, vol. 69, no. 3, pp. 579–593, Mar. 2021.
View at: Publisher Site | Google Scholar
S. H. Halani, S. Yousefi, J. Velazquez Vega et al., “Multi-faceted computational assessment of risk and progression in oligodendroglioma implicates NOTCH and PI3K pathways,” Npj Precision Oncology, vol. 2, no. 1, p. 24, 2018.
View at: Publisher Site | Google Scholar
C. Giachino, J. L. Boulay, R. Ivanek et al., “A tumor suppressor function for Notch signaling in Forebrain tumor subtypes,” Cancer Cell, vol. 28, no. 6, pp. 730–742, 2015.
View at: Publisher Site | Google Scholar
Q. Liu, W. Yu, S. Zhu et al., “Long noncoding RNA GAS5 regulates the proliferation, migration, and invasion of glioma cells by negatively regulating miR-18a-5p,” Journal of Cellular Physiology, vol. 234, no. 1, pp. 757–768, 2018.
View at: Publisher Site | Google Scholar
J.-F. Huo and X.-B. Chen, “Long noncoding RNA growth arrest-specific 5 facilitates glioma cell sensitivity to cisplatin by suppressing excessive autophagy in an mTOR-dependent manner,” Journal of Cellular Biochemistry, vol. 120, no. 4, pp. 6127–6136, 2019.
View at: Publisher Site | Google Scholar
T. Shahar, A. Granit, D. Zrihan et al., “Expression level of miRNAs on chromosome 14q32.31 region correlates with tumor aggressiveness and survival of glioblastoma patients,” Journal of Neuro-Oncology, vol. 130, no. 3, pp. 413–422, Dec. 2016.
View at: Publisher Site | Google Scholar
X. Liu, Q. Zhu, Y. Guo, Z. Xiao, L. Hu, and Q. Xu, “LncRNA LINC00689 promotes the growth, metastasis and glycolysis of glioma cells by targeting miR-338-3p/PKM2 axis,” Biomedicine & Pharmacotherapy, vol. 117, Article ID 109069, 2019.
View at: Publisher Site | Google Scholar
H. Shen, X. D. Weng, X. H. Liu et al., “miR-181a-5p is downregulated and inhibits proliferation and the cell cycle in prostate cancer,” International Journal of Clinical and Experimental Pathology, vol. 11, no. 8, pp. 3969–3976, 2018.
View at: Google Scholar
J. Wang, H. Wang, A. Liu, C. Fang, J. Hao, and Z. Wang, “Lactate dehydrogenase A negatively regulated by miRNAs promotes aerobic glycolysis and is increased in colorectal cancer,” Oncotarget, vol. 6, no. 23, Article ID 19468, 2015.
View at: Publisher Site | Google Scholar
W. Han, J. Shi, J. Cao, B. Dong, and W. Guan, “Emerging roles and therapeutic Interventions of aerobic glycolysis in glioma,” OncoTargets and Therapy, vol. 13, pp. 6937–6955, 2020.
View at: Publisher Site | Google Scholar
A. M. Reuss, D. Groos, M. Buchfelder, and N. Savaskan, “The Acidic brain-glycolytic Switch in the microenvironment of malignant glioma,” International Journal of Molecular Sciences, vol. 22, no. 11, p. 5518, May 2021.
View at: Publisher Site | Google Scholar
W. H. Fridman, F. Pagès, C. Sautès-Fridman, and J. Galon, “The immune contexture in human tumours: impact on clinical outcome,” Nature Reviews Cancer, England, vol. 12, no. 4, pp. 298–306, 2012.
View at: Publisher Site | Google Scholar
N. S. Joshi, E. Akama-Garren, Y. Lu et al., “Regulatory T cells in tumor-associated tertiary lymphoid structures suppress anti-tumor T cell responses,” Immunity, vol. 43, no. 3, pp. 579–590, Sep. 2015.
View at: Publisher Site | Google Scholar
M. Jinushi, S. Chiba, H. Yoshiyama et al., “Tumor-associated macrophages regulate tumorigenicity and anticancer drug responses of cancer stem/initiating cells,” Proceedings of the National Academy of Sciences of the United States of America, vol. 108, no. 30, pp. 12425–12430, 2011.
View at: Publisher Site | Google Scholar
R. Bhatia, S. K. Gautam, A. Cannon et al., “Cancer-associated mucins: role in immune modulation and metastasis,” Cancer & Metastasis Reviews, vol. 38, no. 1–2, pp. 223–236, 2019.
View at: Publisher Site | Google Scholar
M. Ohta, A. Ishida, M. Toda et al., “Immunomodulation of monocyte-derived dendritic cells through ligation of tumor-produced mucins to Siglec-9,” Biochemical and Biophysical Research Communications, vol. 402, no. 4, pp. 663–669, Nov. 2010.
View at: Publisher Site | Google Scholar
Y. Sunakawa, S. Stremitzer, S. Cao et al., “Association of variants in genes encoding for macrophage-related functions with clinical outcome in patients with locoregional gastric cancer,” Annals of Oncology, vol. 26, no. 2, pp. 332–339, Feb. 2015.
View at: Publisher Site | Google Scholar
B.-Z. Qian and J. W. Pollard, “Macrophage diversity enhances tumor progression and metastasis,” Cell, vol. 141, no. 1, pp. 39–51, 2010.
View at: Publisher Site | Google Scholar
A. Kakizaki, T. Fujimura, S. Furudate et al., “Immunomodulatory effect of peritumorally administered interferon-beta on melanoma through tumor-associated macrophages,” OncoImmunology, vol. 4, no. 11, p. e1047584, 2015.
View at: Publisher Site | Google Scholar
R. Gupta, R. Webb-Myers, S. Flanagan, and M. E. Buckland, “Isocitrate dehydrogenase mutations in diffuse gliomas: clinical and aetiological implications,” Journal of Clinical Pathology, vol. 64, no. 10, pp. 835–844, Oct. 2011.
View at: Publisher Site | Google Scholar
Z. J. Reitman and H. Yan, “Isocitrate dehydrogenase 1 and 2 mutations in cancer: alterations at a crossroads of cellular metabolism,” JNCI Journal of the National Cancer Institute, vol. 102, no. 13, pp. 932–941, Jul. 2010.
View at: Publisher Site | Google Scholar
B. Murnyak and L. E. Huang, “Association of TP53 alteration with tissue Specificity and patient outcome of IDH1-mutant glioma,” Cells, vol. 10, no. 8, p. 2116, 2021.
View at: Publisher Site | Google Scholar
F. E. Bleeker, N. A. Atai, S. Lamba et al., “The prognostic IDH1(R132) mutation is associated with reduced NADP+-dependent IDH activity in glioblastoma,” Acta Neuropathologica, vol. 119, no. 4, pp. 487–494, Apr. 2010.
View at: Publisher Site | Google Scholar
J. Balss, J. Meyer, W. Mueller, A. Korshunov, C. Hartmann, and A. von Deimling, “Analysis of the IDH1 codon 132 mutation in brain tumors,” Acta Neuropathologica, vol. 116, no. 6, pp. 597–602, Dec. 2008.
View at: Publisher Site | Google Scholar
H. Yan, D. W. Parsons, G. Jin et al., “IDH1 and IDH2 mutations in gliomas,” New England Journal of Medicine, vol. 360, no. 8, pp. 765–773, 2009.
View at: Publisher Site | Google Scholar
N. M. Amankulor, Y. Kim, S. Arora et al., “Mutant IDH1 regulates the tumor-associated immune system in gliomas,” Genes & Development, vol. 31, no. 8, pp. 774–786, Apr. 2017.
View at: Publisher Site | Google Scholar
L. d. S. Chagas, P. C. Sandre, N. C. A. Ribeiro e Ribeiro et al., “Environmental signals on microglial function during brain development, Neuroplasticity, and disease,” International Journal of Molecular Sciences, vol. 21, no. 6, p. 2111, 2020.
View at: Publisher Site | Google Scholar
J. Xiang, L. Zhou, X. Li et al., “Preoperative monocyte-to-lymphocyte ratio in Peripheral blood predicts Stages, metastasis, and Histological Grades in patients with ovarian cancer,” Translational Oncology, vol. 10, no. 1, pp. 33–39, 2017.
View at: Publisher Site | Google Scholar
R. M. Bambury, M. Y. Teo, D. G. Power et al., “The association of pre-treatment neutrophil to lymphocyte ratio with overall survival in patients with glioblastoma multiforme,” Journal of Neuro-Oncology, vol. 114, no. 1, pp. 149–154, 2013.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Xin Fan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

360

Downloads

544

Citations