Abstract
PURPOSE: Glioblastoma is one of the most aggressive nervous system neoplasms. Immunotherapy represents a hot spot and has not been included in standard treatments of glioblastoma. So in this study, we aim to filtrate an immune-related gene pairs (IRGPs) signature for predicting survival and immune heterogeneity. METHODS: We used gene expression profiles and clinical information of glioblastoma patients in the TCGA and CGGA datasets, dividing into discovery and validation cohorts. IRGPs significantly correlative with prognosis were selected to conduct an IRGPs signature. Low and high risk groups were separated by this IRGPs signature. Univariate and multivariate cox analysis were adopted to check whether risk can be a independent prognostic factor. Immune heterogeneity between different risk groups was analyzed via immune infiltration and gene set enrichment analysis (GSEA). Some different expressed genes between groups were selected to determine their relationship with immune cells and immune checkpoints. RESULTS: We found an IRGPs signature consisting of 5 IRGPs. Different risk based on IRGPs signature is a independent prognostic factor both in the discovery and validation cohorts. High risk group has some immune positive cells and more immune repressive cells than low risk group by means of immune infiltration. We discovered some pathways are more active in the high risk group, leading to immune suppression, drug resistance and tumor evasion. In two specific signaling, some genes are over expressed in high risk group and positive related to immune repressive cells and immune checkpoints, which indicate aggression and immunotherapy resistance. CONCLUSION: We identified a robust IRGPs signature to predict prognosis and immune heterogeneity in glioblastoma patients. Some potential targets and pathways need to be further researched to make different patients benefit from personalized immunotherapy.