Abstract PURPOSE To comprehensively unravel the heterogeneity of glioblastoma, we developed a subtyping method using integrated imaging and genomic data. METHODS We assembled a retrospective cohort of 571 IDH-wildtype glioblastoma patients, obtaining pre-operative multi-parametric MRI (T1, T1-GD, T2, T2-FLAIR, DSC, DTI) scans (available for 462 patients) and targeted next-generation sequencing (NGS) data (available for 355 patients). We extracted 971 radiomic features from these MRI scans, selecting 12 significant dimensions through L21-norm minimization guided by 13 key driver genes from the five most frequently altered pathways in glioblastoma compared to healthy controls (RB1, P53, MAPK, PI3K, and RTK). Subtypes were identified using a joint clustering approach that integrated radiomics and genomics. RESULTS Our method identified three glioblastoma subtypes with distinct risk levels: high-risk (subtype 1), medium-risk (subtype 2), and low-risk (subtype 3), each characterized by their differential overall survival outcomes (Kaplan-Meier analysis, p < 0.05). The intensities of axial diffusivity in enhancing tumor regions and radial diffusivity in non-enhancing cores increased across three subtypes. All subtypes shared a common tendency for more frequent mutation co-occurrence pattern in [TP53,RB1] but exhibited unique molecular characteristics. For example, subtype 1 displayed a relatively high frequency of co-occurring mutations in [TP53,KDR] and [NOTCH2,MDM4]; subtype 2 had six co-occurring pairs [TP53,KDR], [PTPN11,NF1], [RB1,FUBP1], [PTEN,CIC], [PTEN,KRAS], and [NOTCH2,CDKN2A]; subtype 3 showed six co-occurring pairs [PTPN11,NF1], [PDGFRA,NF1], [NTRK1,NOTCH2], [MET,ATRX], [KRAS,KDR], and [FGFR2,DNMT3A] plus four significant patterns of mutual exclusivity in [TP53,PTPN11], [TP53,EGFR], [EGFR,RB1], and [EGFR,NF1]. CONCLUSION Our study discovered distinct glioblastoma subtypes, revealing complex patterns and structures across multiple data modalities without relying on prior clinical assumptions. This is promising to help enhance the precision of diagnosis and treatment of glioblastoma.
Read full abstract