Prognostic Features of the Tumor Immune Microenvironment in Glioma and Their Clinical Applications: Analysis of Multiple Cohorts

Abstract

Glioma is the most common malignant tumor of the central nervous system. Tumor purity is a source of important prognostic factor for glioma patients, showing the key roles of the microenvironment in glioma prognosis. In this study, we systematically screened functional characterization related to the tumor immune microenvironment and constructed a risk model named Glioma MicroEnvironment Functional Signature (GMEFS) based on eight cohorts. The prognostic value of the GMEFS model was also verified in another two glioma cohorts, glioblastoma (GBM) and low-grade glioma (LGG) cohorts, from The Cancer Genome Atlas (TCGA). Nomograms were established in the training and testing cohorts to validate the clinical use of this model. Furthermore, the relationships between the risk score, intrinsic molecular subtypes, tumor purity, and tumor-infiltrating immune cell abundance were also evaluated. Meanwhile, the performance of the GMEFS model in glioma formation and glioma recurrence was systematically analyzed based on 16 glioma cohorts from the Gene Expression Omnibus (GEO) database. Based on multiple-cohort integrated analysis, risk subpathway signatures were identified, and a drug–subpathway association network was further constructed to explore candidate therapy target regions. Three subpathways derived from Focal adhesion (path: 04510) were identified and contained known targets including platelet derived growth factor receptor alpha (PDGFRA), epidermal growth factor receptor (EGFR), and erb-b2 receptor tyrosine kinase 2 (ERBB2). In conclusion, the novel functional signatures identified in this study could serve as a robust prognostic biomarker, and this study provided a framework to identify candidate therapeutic target regions, which further guide glioma patients’ clinical decision.