Abstract
Glioblastoma multiforme (GBM) is a highly malignant primary brain tumor with a poor prognosis. Reactive oxygen species that accumulate during tumorigenesis can cause oxidative stress (OS), which plays a crucial role in cancer cell survival. Clinical and transcriptome data of TCGA-GBM dataset from UCSC Xena database were analyzed. Consensus clustering analysis was conducted to identify OS-related molecular subtypes for GBM. The immune infiltrate level between subtypes were characterized by ESTIMATE algorithm. Differentially expressed genes (DEGs) between subtypes were screened by DESeq2 package. Two OS-related molecular subtypes of GBM were identified, and cluster 2 had poorer overall survival and higher immune infiltration levels than cluster 1. Enrichment analysis showed that 54 DEGs in cluster 2 were significantly enriched in cytokine/chemokine-related functions or pathways. Ten hub genes (CSF2, CSF3, CCL7, LCN2, CXCL6, MMP8, CCR8, TNFSF11, IL22RA2, and ORM1) were identified in GBM subtype 2 through protein-protein interaction network, most of which were positively correlated with immune factors and immune checkpoints. A total of 55 small molecule drugs obtained from drug gene interaction database (DGIdb) may have potential therapeutic effects in GBM subtype 2 patients. Our study identified 10 hub genes as potential therapeutic targets in GBM subtype 2 patients, who have poorer overall survival and higher immune infiltration levels. These findings could pave the way for new treatments for this aggressive form of brain cancer.