Multi-sampling allows intra-tumoral heterogeneity querying and vulnerability profiling in glioblastoma

Glioblastoma (GBM) remains a devastating cancer with limited treatment options, largely due to its heterogeneity. While supramaximal resection has recently provided survival benefits, therapeutic profiling of different tumor compartments, particularly its infiltrative edge remains largely unexplored.

Here, we leveraged magnetic resonance imaging (MRI)-guided multi-sampling, collecting 2 cores and 2 margins per case, to query GBM heterogeneity. Whole-exome and RNA-seq with drug testing in two patient-derived 3D models were used to reveal similarities and differences in genomic and transcriptomic makeups, cellular compositions, and drug responses across cores and margins. Bioinformatics interrogations further identified response biomarkers.

Mutation analysis showed that oncogenes exhibited a higher degree of spatial heterogeneity than tumor suppressor genes, regardless of MRI status. While the mesenchymal transcriptional subtype with extracellular matrix remodeling, stress response, and immune programs were preferentially enriched in enhancing cores, proneural tumors with neurological processes favored non-enhancing margins. Using a 15-drug GBM-targeted panel, ERK (ulixertinib) and PI3K pathway (paxalisib, CC-115) inhibitors showed preferential efficacy in enhancing cores and non-enhancing margins, respectively. The anti-apoptosis, pan-Bcl2 agent navitoclax and the epigenetic drug trotabresib represented the most effective, tumor-wide monotherapies. Importantly, drug combinations generally outperformed single agents across all regions.

This work demonstrates the regional heterogeneity of therapeutic vulnerabilities in GBM ex vivo, showing various drugs with tumor-wide or MRI-enhancement informed activity. These findings offer preclinical bases of numerous monotherapies and drug combinations for future clinical trial design.