Tumor-Associated Tractography Derived from High-Angular-Resolution Q-Space MRI May Predict Patterns of Cellular Invasion in Glioblastoma

Background: The invasion of glioblastoma cells beyond the visible tumor margin depicted by conventional neuroimaging is believed to mediate recurrence and predict poor survival. Radiomic biomarkers that are associated with the direction and extent of tumor infiltration are, however, non-existent.

Methods: Patients from a single center with newly diagnosed glioblastoma (n = 7) underwent preoperative Q-space magnetic resonance imaging (QSI; 3T, 64 gradient directions, b = 1000 s/mm²) between 2018 and 2019. Tumors were manually segmented, and patterns of inter-voxel coherence spatially intersecting each segmentation were generated to represent tumor-associated tractography. One patient additionally underwent regional biopsy of diffusion tract- versus non-tract-associated tissue during tumor resection for RNA sequencing. Imaging data from this cohort were compared with a historical cohort of n = 66 glioblastoma patients who underwent similar QSI scans. Associations of tractography-derived metrics with survival were assessed using t-tests, linear regression, and Kaplan-Meier statistics. Patient-derived glioblastoma xenograft (PDX) mice generated with the sub-hippocampal injection of human-derived glioblastoma stem cells (GSCs) were scanned under high-field conditions (QSI, 7T, 512 gradient directions), and tumor-associated tractography was compared with the 3D microscopic reconstruction of immunostained GSCs.

Results: In the principal enrollment cohort of patients with glioblastoma, all cases displayed tractography patterns with tumor-intersecting tract bundles extending into brain parenchyma, a phenotype which was reproduced in PDX mice as well as in a larger comparison cohort of glioblastoma patients (n = 66), when applying similar methods. Reconstructed spatial patterns of GSCs in PDX mice closely mirrored tumor-associated tractography. On a Kaplan-Meier survival analysis of n = 66 patients, the calculated intra-tumoral mean diffusivity predicted the overall survival (p = 0.037), as did tractography-associated features including mean tract length (p = 0.039) and mean projecting tract length (p = 0.022). The RNA sequencing of human tissue samples (n = 13 tumor samples from a single patient) revealed the overexpression of transcripts which regulate cell motility in tract-associated samples.

Conclusions: QSI discriminates tumor-specific patterns of inter-voxel coherence believed to represent white matter pathways which may be susceptible to glioblastoma invasion. These findings may lay the groundwork for future work on therapeutic targeting, patient stratification, and prognosis in glioblastoma.