Glioblastoma multiforme (GBM) is the most aggressive form of brain tumor characterized by inter and intra-tumor heterogeneity and complex tumor microenvironment. To uncover the molecular targets in this milieu, we systematically identified immune and stromal interactions at the glial cell type level that leverages on RNA-sequencing data of GBM patients from The Cancer Genome Atlas. The perturbed genes between the high vs low immune and stromal scored patients were subjected to weighted gene co-expression network analysis to identify the glial cell type specific networks in immune and stromal infiltrated patients. The intramodular connectivity analysis identified the highly connected genes in each module. Combining it with univariable and multivariable prognostic analysis revealed common vital gene ITGB2, between the immune and stromal infiltrated patients enriched in microglia and newly formed oligodendrocytes. We found following unique hub genes in immune infiltrated patients; COL6A3 (microglia), ITGAM (oligodendrocyte precursor cells), TNFSF9 (microglia), and in stromal infiltrated patients, SERPINE1 (microglia) and THBS1 (newly formed oligodendrocytes, oligodendrocyte precursor cells). To validate these hub genes, we used external GBM patient single cell RNA-sequencing dataset and this identified ITGB2 to be significantly enriched in microglia, newly formed oligodendrocytes, T-cells, macrophages and adipocyte cell types in both immune and stromal datasets. The tumor infiltration analysis of ITGB2 showed that it is correlated with myeloid dendritic cells, macrophages, monocytes, neutrophils, B-cells, fibroblasts and adipocytes. Overall, the systematic screening of tumor microenvironment components at glial cell types uncovered ITGB2 as a potential target in primary GBM.
Keywords: Co-expression networks; Glial cell-types; Glioblastoma multiforme; Prognostic marker; Tumor microenvironment.
Copyright © 2023 Elsevier Ltd. All rights reserved.