A data-driven intravoxel mean diffusivities distribution approach for molecular classifications and MIB-1 prediction of gliomas.

ABSTRACT

BACKGROUND: Measuring non-parametric intravoxel mean diffusivity distributions (MDDs) using magnetic resonance imaging (MRI) is a sensitive method for detecting intracellular diffusivity changes during physiological alterations. Histological and molecular glioma classifications are essential for prognosis and treatment, with distinct water diffusion dynamics among subtypes. PURPOSE: We developed a data-driven approach using a fully connected network (FCN) to enhance the speed and stability of calculating MDDs across varying SNRs, enable tumor microstructural mapping, and test its reliability in identifying MIB-1 labeling index (LI) levels and molecular status of gliomas. METHODS: An FCN was trained to learn the mapping between the simulated diffusion decay curves and the ground truth MDDs. We performed 5 000 000 simulation curves with various diffusivity components and random SNR ∈ [ 30 , 300 ] $ \in [ {30,\ 300} ]$ . Eighty percent of simulation curves were used for the FCN training, 10% for validation, and the others were external tests for the FCN performance evaluation. In vivo data were collected to evaluate its clinical reliability. One hundred one patients (44 years ± $ \pm $ 14, 67 men) with gliomas and six healthy controls underwent a 3.0 T MRI examination with a spin echo-echo planar imaging (SE-EPI) diffusion-weighted imaging (DWI) sequence. The trained FCN was employed to calculate MDDs of each brain voxel by voxel. We used the Fuzzy C-means algorithm to cluster the MDDs of tumor voxels, facilitating the characterization of distinct glioma tissues. Quantitative assessments were conducted through sectional integrals of the MDDs, demarcated by six bands to derive signal fractions ( f n , n = 1 - 6 ${{f}_n},\ n = 1 -6$ ) and diffusivities of the maximum peaks ( D p e a k ${{D}_{peak}}$ ). Cosine similarity scores (CSS) were used for MDD similarity. ANOVA and Mann-Whitney U test were used for difference analysis. Logistic regression and area under the receiver operator characteristic curve (AUC) were used for classification evaluation. RESULTS: The simulation results showed that the FCN-based MDD approach (FCN-MDD) achieved higher CSS than non-negative least squares-based MDD (NNLS-MDD). For in vivo data, the spectra of ET and NET obtained by FCN-MDD are more distinguishable than NNLS-MDD. Fraction maps delineate the characteristics of different tumor tissues (enhancing and non-enhancing tumor, edema, and necrosis). f 3 , f 4 , D p e a k ${{f}_3},\ {{f}_4},{{D}_{peak}}$ showed a positive and negative correlation with MIB-1 respectively ( r = 0.568 , r = - 0.521 , r = - 0.654 $r = 0.568,\ r = - 0.521,\ r = - 0.654$ , all p < 0.001 $p < 0.001$ ). The AUC of D p e a k ${{D}_{peak}}$ for predicting MIB-1 LI levels was 0.900 (95% CI, 0.826-0.974), versus 0.781 (0.677-0.886) of ADC. The highest AUC of isocitrate dehydrogenase (IDH) mutation status, assessed by a logistic regression model ( f 1 + f 3 ${{f}_1} + {{f}_3}$ ) was 0.873 (95% CI, 0.802-0.944). CONCLUSION: The proposed FCN-MDD method was more robust to variations in SNR and less reliant on empirically set regularization values than the NNLS-MDD method. FCN-MDD also enabled qualitative and quantitative evaluation of the composition of gliomas.