Abstract
Glioblastoma (GBM) is the most malignant and prevalent primary brain tumor. In this study, weighted gene coexpression network analysis (WGCNA) was performed to analyze RNA binding protein (RBP) expression data from The Cancer Genome Atlas (TCGA) for the IDH-wild type GBM cohort. The CIBERSORT algorithm quantified the cellular composition of immune cells and was used to identify key modules associated with CD8+ T cell infiltration. Coexpression networks analysis and protein-protein interaction (PPI) network analysis was used to filter out central RBP genes. Eleven RBP genes, including MYEF2, MAPT, NOVA1, MAP2, TUBB2B, CDH10, TTYH1, PTPRZ1, SOX2, NOVA2 and SCG3, were identified as candidate CD8+ T cell infiltration-associated central genes. A Cox proportional hazards regression model and Kaplan-Meier analysis were applied to identify candidate biomarkers. MYEF2 was selected as a prognostic biomarker based on the results of prognostic analysis. Flow Cytometric Analysis indicated that MYEF2 expression was negatively correlated with dysfunctional CD8+ T cell markers. Kaplan-Meier survival analysis (based on IHC staining) revealed that GBM patients with elevated MYEF2 expression have a better prognosis. Knockdown of MYEF2 in GBM cells via in vitro assays was observed to promote cell proliferation and migration. Our study suggests that MYEF2 expression negatively correlates with T cell exhaustion and tumor progression, rendering it a potentially valuable prognostic biomarker for GBM.