Introduction. Glioblastoma exhibits profound tumor heterogeneity, which causes inconsistent treatment response. The aim of this study was to propose an interpretation method of diffusion tensor imaging (DTI) using joint histogram analysis of DTI-p and -q. With this method we explored the patterns of tumor infiltration which causes disruption of brain microstructure, and examined the prognostic value of tumor infiltrative patterns for patient survival. Materials and methods. A total of 115 primary glioblastoma patients (mean age 59.3 years, 87 males) were prospectively recruited from July 2010 to August 2015. Patients underwent preoperative MRI scans and maximal safe resection. DTI was processed and decomposed into p and q components. The univariate and joint histograms of DTI-p and -q were constructed using the voxels of contrast-enhancing and non-enhancing regions respectively. Eight joint histogram features were obtained and correlated with tumor progression and patient survival. Their prognostic values were compared with clinical factors using receiver operating characteristic curves. Results. The subregion of increased DTI-p and decreased DTI-q accounted for the largest proportion. Additional diffusion patterns can be identified via joint histogram analysis. Particularly, higher proportion of decreased DTI-p and increased DTI-q in non-enhancing region contributed to worse progression-free survival (hazard ratio = 1.08, p < 0.001) and overall survival (hazard ratio = 1.11, p < 0.001). Conclusions. Joint histogram analysis of DTI can provide a comprehensive measure of tumor infiltration and microstructure change, which showed prognostic values. The subregion of decreased DTI-p and increased DTI-q in non-enhancing regions may indicate a more invasive habitat.

Introduction: Glioblastoma is characterized by its remarkable heterogeneity and dismal prognosis. Histogram analysis of quantitative magnetic resonance imaging (MRI) is an important in vivo method to study intratumoral heterogeneity. With large amounts of histogram features generated, integrating these modalities effectively for clinical decision remains a challenge. Methods: A total of 80 patients with supratentorial primary glioblastoma were recruited. All patients received surgery and standard regimen of temozolomide chemoradiotherapy. Diagnosis was confirmed by pathology. Anatomical T2-weighted, T1-weighted post-contrast and FLAIR images, as well as dynamic susceptibility contrast (DSC), diffusion tensor imaging (DTI) and chemical shift imaging were acquired preoperatively using a 3T MRI scanner. DTI-p, DTI-q, relative cerebral blood volume (rCBV), mean transit time (MTT) and relative cerebral blood flow (rCBF) maps were generated. Contrast-enhancing (CE) and non-enhancing (NE) regions of interest were manually delineated. Voxel intensity histograms were constructed from the CE and NE regions independently. Patient clustering was performed by the Multi-View Biological Data Analysis (MVDA) approach. Kaplan-Meier and Cox proportional hazards regression analyses were performed to evaluate the relevance of the patient clustering to survival. The histogram features selected from MVDA approach were evaluated using receiver operator characteristics (ROC) curve analysis. The metabolic signatures of the patient clusters were analyzed by multivoxel MR spectroscopy (MRS). Results: The MVDA approach yielded two final patient clusters, consisting of 53 and 27 patients respectively. The two patient subgroups showed significance for overall survival (p = 0.007, HR = 0.32) and progression-free survival (p < 0.001, HR = 0.33) in multivariate Cox regression analysis. Among the features selected by MVDA, higher mean value of DTI-q in the non-enhancing region contributed to a worse OS (HR = 1.40, p = 0.020) and worse PFS (HR = 1.36, p = 0.031). Multivoxel MRS showed N-acetylaspartate/creatine (NAA/Cr) ratio between the two clusters, both in the CE region (p < 0.001) and NE region (p = 0.013). Glutamate/Cr (Glu/Cr) ratio and glutamate + glutamine/Cr (Glx/Cr) of the cluster 1 was significantly lower than cluster 2 (p = 0.037, and 0.027 respectively) In the NE region. Discussion: This study demonstrated that integrating multi-parametric and multi-regional MRI histogram features may help to stratify patients. The histogram features selected from the proposed approach may be used as potential imaging markers in personalized treatment strategy and response determination.

Background: Glioblastoma exhibits profound intratumoral heterogeneity in blood perfusion, which may cause inconsistent therapy response. Particularly, low perfusion may create hypoxic microenvironment and induce resistant clones. Thus, developing validated imaging approaches that define low perfusion compartments is crucial for clinical management. Methods: A total of 112 newly-diagnosed supratentorial glioblastoma patients were prospectively recruited for maximal safe resection. Preoperative MRI included anatomical, dynamic susceptibility contrast (DSC), diffusion tensor imaging (DTI) and chemical shift imaging (CSI). The apparent diffusion coefficient (ADC) and relative cerebral blood volume (rCBV) were calculated from DTI and DSC respectively. Using thresholding methods, two low perfusion compartments (ADCH-rCBVL and ADCL-rCBVL) were identified. Volumetric analysis was performed. Lactate and macromolecule/lipid levels were determined from multivoxel spectroscopy. Progression-free survival (PFS) and overall survival (OS) were analysed using Kaplan-Meier and multivariate Cox regression analyses. Results: Two compartments displayed higher lactate and macromolecule/lipid levels than normal controls (each P < 0.001), suggesting hypoxic and pro-inflammatory microenvironment. The proportional volume of ADCL-rCBVL compartment was associated with a larger infiltration area (P < 0.001, rho = 0.42). Lower lactate in this compartment was associated with a less invasive phenotype visualized on DTI. Multivariate Cox regression showed higher lactate level in the ADCL-rCBVL compartment was associated with a worse survival (PFS: HR 2.995, P = 0.047; OS: HR 4.974, P = 0.005). Conclusions: The ADCL-rCBVL compartment represent a treatment resistant sub-region associated with glioblastoma invasiveness. This approach was based on clinically available imaging modalities and could thus provide crucial pretreatment information for clinical decision making.