Screening seven hub genes associated with prognosis and immune infiltration in glioblastoma

Abstract

Glioblastoma (GBM) is the most common and deadly primary brain tumor in adults. Diagnostic and therapeutic challenges have been raised because of poor prognosis. Gene expression profiles of GBM and normal brain tissue samples from GSE68848, GSE16011, GSE7696, and The Cancer Genome Atlas (TCGA) were downloaded. We identified differentially expressed genes (DEGs) by differential expression analysis and obtained 3,800 intersected DEGs from all datasets. Enrichment analysis revealed that the intersected DEGs were involved in the MAPK and cAMP signaling pathways. We identified seven different modules and 2,856 module genes based on the co-expression analysis. Module genes were used to perform Cox and Kaplan-Meier analysis in TCGA to obtain 91 prognosis-related genes. Subsequently, we constructed a random survival forest model and a multivariate Cox model to identify seven hub genes (KDELR2, DLEU1, PTPRN, SRBD1, CRNDE, HPCAL1, and POLR1E). The seven hub genes were subjected to the risk score and survival analyses. Among these, CRNDE may be a key gene in GBM. A network of prognosis-related genes and the top three differentially expressed microRNAs with the largest fold-change was constructed. Moreover, we found a high infiltration of plasmacytoid dendritic cells and T helper 17 cells in GBM. In conclusion, the seven hub genes were speculated to be potential prognostic biomarkers for guiding immunotherapy and may have significant implications for the diagnosis and treatment of GBM.