Single cell RNA sequencing reveals differentiation related genes with drawing implications in predicting prognosis and immunotherapy response in gliomas

Differentiation states of glioma cells correlated with prognosis and tumor-immune microenvironment (TIME) in patients with gliomas. We aimed to identify differentiation related genes (DRGs) for predicting the prognosis and immunotherapy response in patients with gliomas. We identified three differentiation states and the corresponding DRGs in glioma cells through single-cell transcriptomics analysis. Based on the DRGs, we separated glioma patients into three clusters with distinct clinicopathological features in combination with bulk RNA-seq data. Weighted correlation network analysis, univariate cox regression analysis and least absolute shrinkage and selection operator analysis were involved in the construction of the prognostic model based on DRGs. Distinct clinicopathological characteristics, TIME, immunogenomic patterns and immunotherapy responses were identified across three clusters. A DRG signature composing of 12 genes were identified for predicting the survival of glioma patients and nomogram model integrating the risk score and multi-clinicopathological factors were constructed for clinical practice. Patients in high-risk group tended to get shorter overall survival and better response to immune checkpoint blockage therapy. We obtained 9 candidate drugs through comprehensive analysis of the differentially expressed genes between the low and high-risk groups in the model. Our findings indicated that the risk score may not only contribute to the determination of prognosis but also facilitate in the prediction of immunotherapy response in glioma patients.