# Integrated pathway analysis identifies prognostically relevant subtypes of glioblastoma characterized by abnormalities in multi‐omics

**Authors:** Pei Zhang, Dan Liu, Tonghui Yu, Yanlin Zhang, Lu Zhong, Xiao Ouyang, Qin Xia, Lei Dong

PMC · DOI: 10.1002/ctm2.70517 · 2025-11-25

## TL;DR

This study identifies two new glioblastoma subtypes based on tumor-driving mutations and immune cell presence, using multi-omics data to better understand the disease and its treatment.

## Contribution

The study introduces two newly identified glioblastoma subtypes through integrated multi-omics analysis, revealing their distinct biological and therapeutic features.

## Key findings

- Two new glioblastoma subtypes were identified: tumor-driving and immune-blockage.
- The immune-inflamed subtype shows upregulated immune-related pathways like PD-1 and IFN-γ signaling.
- IDH1 mutations co-occur with TP53, suggesting a protective role in TP53 mutant patients.

## Abstract

Gene expression‐based molecular subtypes in glioblastoma from The Cancer Genome Atlas Network (TCGA‐GBM) unraveled the pathological origins by identifying tumour cell driver genes. However, the causal inference between molecular subtype origins and their therapeutic efficacy remains obscure.

We integrated TCGA‐GBM multi‐omics (DNA, mRNA, and protein profiles) using correlation analysis to identify cis‐regulation. We analyzed the exposure‐mediated base substitution‐level mutations and their potential triggers. Importantly, we performed Consensus Clustering based on the MSigDB database with Silhouette‐correction to identify prognostically relevant pathway‐based MSig subtypes. The tumour driver mutations (co‐occurrence mutation pattern), aberrant pathways (tumour hallmarks), immune microenvironment (xCell), and pseudo‐time analysis (dyno) were used to characterize the MSig subtype landscape. Furthermore, we evaluated potential drug sensitivities across MSig subtypes using the Genomics of Drug Sensitivity in Cancer database.

We classified five MSig subtypes, characterized by neural‐like, tumour‐driving, low tumour evolution, immune‐inflamed, and classical tumour features. We observed several key features in ‘tumour‐driving’ GBM patients: (1) mutual exclusivity between prognostic factors TP53 and EGFR; and (2) IDH1 mutations co‐occurring with TP53, which account for the protective role of IDH1 in TP53 mutant patients. The immune‐inflamed GBM, characterized as a ‘hot’ tumour, exhibited upregulation of immune‐related pathways, including PD‐1 and IFN‐γ signalling responses. DNA methylation landscape revealed 14 MGMT CpG‐rich regions regulating expression. Evolutionary trajectories revealed progression from a primary tumour state (close to normal tissue) to two distinct endpoints (tumour‐driving and immune‐inflamed subtypes).

Our findings reveal interactions between tumour cells and their surrounding immune environment, classifying GBM into two newly identified subtypes: (1) the tumour‐driving subtype is characterized by multiple oncogenic mutations, while (2) the immune‐blockade subtype is marked by a high presence of immune cells. We highlight the importance of integrating multi‐type data (somatic mutations, DNA methylation, and RNA transcripts, etc.) to decipher GBM biology and potential therapeutic implications.

We report the interaction between tumor cells and environmental immune cells, classifying GBM into two main subtypes: 1) The tumor‐driving subtype is characterized by multiple oncogenic mutations, while 2) the immune‐blockage subtype is marked by a high presence of immune cells. We used integrated multidimensional analyses of somatic mutations, DNA methylation, and RNA transcripts to gain a deeper understanding of GBM biology and potential therapeutic implications.

We report the interaction between tumour cells and environmental immune cells, classifying GBM into two main subtypes: (1) The tumour‐driving subtype is characterized by multiple oncogenic mutations, while (2) the immune‐blockage subtype is marked by a high presence of immune cells. We used integrated multi‐dimensional analyses of somatic mutations, DNA methylation, and RNA transcripts to gain a deeper understanding of GBM biology and potential therapeutic implications.

## Linked entities

- **Genes:** TP53 (tumor protein p53) [NCBI Gene 7157], EGFR (epidermal growth factor receptor) [NCBI Gene 1956], IDH1 (isocitrate dehydrogenase (NADP(+)) 1) [NCBI Gene 3417], MGMT (O-6-methylguanine-DNA methyltransferase) [NCBI Gene 4255]
- **Diseases:** glioblastoma (MONDO:0018177)

## Full-text entities

- **Genes:** IFNG (interferon gamma) [NCBI Gene 3458] {aka IFG, IFI, IMD69}, TP53 (tumor protein p53) [NCBI Gene 7157] {aka BCC7, BMFS5, LFS1, P53, TRP53}, PDCD1 (programmed cell death 1) [NCBI Gene 5133] {aka ADMIO4, AIMTBS, CD279, PD-1, PD1, SLEB2}, MGMT (O-6-methylguanine-DNA methyltransferase) [NCBI Gene 4255], IDH1 (isocitrate dehydrogenase (NADP(+)) 1) [NCBI Gene 3417] {aka HEL-216, HEL-S-26, IDCD, IDH, IDP, IDPC}, EGFR (epidermal growth factor receptor) [NCBI Gene 1956] {aka ERBB, ERBB1, ERRP, HER1, NISBD2, NNCIS}
- **Diseases:** GBM (MESH:D005910), glioblastoma (MESH:D005909), Cancer (MESH:D009369)
- **Species:** Homo sapiens (human, species) [taxon 9606]

## Figures

11 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12647366/full.md

---
Source: https://tomesphere.com/paper/PMC12647366