Identification of an Immune signature assisted prognosis, and immunotherapy prediction for IDH wildtype glioblastoma

IDH-wildtype glioblastoma (GBM) is the most common and malignant primary brain tumor. The purpose of this study is to establish a prognostic gene signature for IDH-wildtype GBM. RNA sequencing data of normal brain tissue and GBM patients were obtained from TCGA, CGGA, GEO and the GTEx databases. Identification of prognostic differentially expressed genes (DEGs) with | log2 fold change | > 0.5 and adjust p < 0.05 in TCGA and CGGA databases by "limma" method. By LASSO regression analysis and multivariate Cox analysis, a 3-gene prognostic signature composed of FMOD, MXRA5 and RAB36 was established. The 3-gene prognostic risk model is validated by TCGA and GSE43378 datasets. The expression of FMOD, MXRA5 and RAB36 in GBM patients was significantly higher than that in normal brain tissues in CCGA, TCGA and GSE29796 data sets. In order to further verify this result, total RNA was extracted from tumors and paracancerous tissues of 9 GBM patients. RT-PCR results showed that the expression of FMOD, MXRA5 and RAB36 in tumor tissues of most patients was higher than that in paracancerous tissues. The results of GSEA showed that the pathway enrichment of the 3-gene signature was mainly related to tumor immunity. Immune cell infiltration analyzed by ssGSEA showed that there were significant differences in macrophages between high- and low-risk groups. Immune checkpoint genes correlation analysis showed that PD-L1 gene expression is closely related to risk score. Our study identifies a prognostic-associated risk model and provides a potential effective immunotherapy target for IDH-wildtype GBM patients.