An Integrated Immune-Related Bioinformatics Analysis in Glioma: Prognostic Signature’s Identification and Multi-Omics Mechanisms’ Exploration

As the traditional treatment for glioma, the most common central nervous system malignancy with poor prognosis, the efficacy of high-intensity surgery combined with radiotherapy and chemotherapy is not satisfactory. The development of individualized scientific treatment strategy urgently requires the guidance of signature with clinical predictive value. In this study, five prognosis-related differentially expressed immune-related genes (PR-DE-IRGs) (CCNA2, HMGB2, CASP3, APOBEC3C, and BMP2) highly associated with glioma were identified for a prognostic model through weighted gene co-expression network analysis, univariate Cox and lasso regression. Kaplan-Meier survival curves, receiver operating characteristic curves and other methods have shown that the model has good performance in predicting the glioma patients’ prognosis. Further combined nomogram provided better predictive performance. The signature’s guiding value in clinical treatment has also been verified by multiple analysis results. We also constructed a comprehensive competing endogenous RNA (ceRNA) regulatory network based on the protective factor BMP2 to further explore its potential role in glioma progression. Numerous immune-related biological functions and pathways were enriched in a high-risk population. Further multi-omics integrative analysis revealed a strong correlation between tumor immunosuppressive environment/IDH1 mutation and signature, suggesting that their cooperation plays an important role in glioma progression.