Patterns of gray and white matter functional networks involvement in glioblastoma patients: indirect mapping from clinical MRI scans

Background: Resting-state functional-MRI studies identified several cortical gray matter functional networks (GMNs) and white matter functional networks (WMNs) with precise anatomical localization. Here, we aimed at describing the relationships between brain’s functional topological organization and glioblastoma (GBM) location. Furthermore, we assessed whether GBM distribution across these networks was associated with overall survival (OS). Materials and methods: We included patients with histopathological diagnosis of IDH-wildtype GBM, presurgical MRI and survival data. For each patient, we recorded clinical-prognostic variables. GBM core and edema were segmented and normalized to a standard space. Pre-existing functional connectivity-based atlases were used to define network parcellations: 17 GMNs and 12 WMNs were considered in particular. We computed the percentage of lesion overlap with GMNs and WMNs, both for core and edema. Differences between overlap percentages were assessed through descriptive statistics, ANOVA, post-hoc tests, Pearson’s correlation tests and canonical correlations. Multiple linear and non-linear regression tests were employed to explore relationships with OS. Results: 99 patients were included (70 males, mean age 62 years). The most involved GMNs included ventral somatomotor, salient ventral attention and default-mode networks; the most involved WMNs were ventral frontoparietal tracts, deep frontal white matter, and superior longitudinal fasciculus system. Superior longitudinal fasciculus system and dorsal frontoparietal tracts were significantly more included in the edema (p < 0.001). 5 main patterns of GBM core distribution across functional networks were found, while edema localization was less classifiable. ANOVA showed significant differences between mean overlap percentages, separately for GMNs and WMNs (p-values<0.0001). Core-N12 overlap predicts higher OS, although its inclusion does not increase the explained OS variance. Discussion and conclusion: Both GBM core and edema preferentially overlap with specific GMNs and WMNs, especially associative networks, and GBM core follows five main distribution patterns. Some inter-related GMNs and WMNs were co-lesioned by GBM, suggesting that GBM distribution is not independent of the brain’s structural and functional organization. Although the involvement of ventral frontoparietal tracts (N12) seems to have some role in predicting survival, network-topology information is overall scarcely informative about OS. fMRI-based approaches may more effectively demonstrate the effects of GBM on brain networks and survival.