Integrative Diffusion-Weighted Imaging and Radiogenomic Network Analysis of Glioblastoma multiforme

In the past, changes of the Apparent Diffusion Coefficient in glioblastoma multiforme have been shown to be related to specific genes and described as being associated with survival. The purpose of this study was to investigate diffusion imaging parameters in combination with genome-wide expression data in order to obtain a comprehensive characterisation of the transcriptomic changes indicated by diffusion imaging parameters. Diffusion-weighted imaging, molecular and clinical data were collected prospectively in 21 patients. Before surgery, MRI diffusion metrics such as axial (AD), radial (RD), mean diffusivity (MD) and fractional anisotropy (FA) were assessed from the contrast enhancing tumour regions. Intraoperatively, tissue was sampled from the same areas using neuronavigation. Transcriptional data of the tissue samples was analysed by Weighted Gene Co-Expression Network Analysis (WGCNA) thus classifying genes into modules based on their network-based affiliations. Subsequent Gene Set Enrichment Analysis (GSEA) identified biological functions or pathways of the expression modules. Network analysis showed a strong association between FA and epithelial-to-mesenchymal-transition (EMT) pathway activation. Also, patients with high FA had a worse clinical outcome. MD correlated with neural function related genes and patients with high MD values had longer overall survival. In conclusion, FA and MD are associated with distinct molecular patterns and opposed clinical outcomes.