Abstract
BACKGROUND: Glioblastoma (GBM) is a molecularly and phenotypically heterogeneous tumor, with regional invasion patterns driving treatment resistance and universal recurrence. In infiltrative regions, GBM cells adopt at least two predominant invasive phenotypes: Neuronal-enriched (NEU), with synaptic and neurodevelopment programs, and Glycolytic/Plurimetabolic (GPM), with heightened immune and metabolic activity. These subtypes differ in therapeutic vulnerability and microenvironmental interactions; features not captured by conventional contrast-enhanced (T1+C) MRI due to minimal blood-brain barrier (BBB) disruption. Our group sought to develop predictive models to distinguish variability of NEU/GPM subtypes within infiltrative GBM, utilizing graph convolutional networks (GCN), advanced multiparametric MRI (mpMRI), and spatially-matched image-localized biopsies with transcriptomic profiling. METHODS: We analyzed 18 mpMRI features —including conventional, diffusion-weighted, and Dynamic Susceptibility Contrast (DSC) perfusion MRI— of 101 image-localized biopsies from MRI non-enhancing infiltrative margins of 64 GBM patient tumors. Features were hierarchically clustered per modality; each cluster was represented as a node in a Neo4j graph. Edges connected nodes sharing overlapping samples. Nodes were annotated with NEU or GPM subtype proportion. GCNs were trained to classify subtype-enriched nodes using graph topology and node features. RESULTS: GCNs achieved high predictive performance on held-out test data (AUC = 0.95 for NEU, 0.88 for GPM; 70:30 train-test split) with similar 5-fold cross-validation performance and generalized well to an independent cohort (AUC > 0.85, n = 45). Importantly, performance was robust without relying on T1+C imaging, enabling prediction beyond BBB-disrupted regions. GCNs outperformed traditional classifiers (logistic regression, random forest). CONCLUSION: Our results demonstrate that GCN-modeling of mpMRI features enables non-invasive prediction of invasive GBM subtypes at the tumor margin. Unlike convolutional neural networks (CNN), which rely on spatial continuity, GCNs leverage relational structure between heterogeneous imaging features, making them ideal for modeling sparse, non-grid data. This may support non-invasive regional vulnerability detection, surgical targeting, and improved response assessment.