White matter microstructural disruption in minimal hepatic encephalopathy: a neurite orientation dispersion and density imaging (NODDI) study.

PURPOSE: To evaluate the ability of neurite orientation dispersion and density imaging (NODDI) for detecting white matter (WM) microstructural abnormalities in minimal hepatic encephalopathy (MHE). METHODS: Diffusion-weighted images, enabling the estimation of NODDI and diffusion tensor imaging (DTI) parameters, were acquired from 20 healthy controls (HC), 22 cirrhotic patients without MHE (NHE), and 15 cirrhotic patients with MHE. Tract-based spatial statistics were used to determine differences in DTI (including fractional anisotropy [FA] and mean/axial/radial diffusivity [MD/AD/RD]) and NODDI parameters (including neurite density index [NDI], orientation dispersion index [ODI], and isotropic volume fraction [ISO]). Voxel-wise analyses of correlations between diffusion parameters and neurocognitive performance determined by Psychometric Hepatic Encephalopathy Score (PHES) were completed. RESULTS: MHE patients had extensive NDI reduction and rare ODI reduction, primarily involving the genu and body of corpus callosum and the bilateral frontal lobe, corona radiata, external capsule, anterior limb of internal capsule, temporal lobe, posterior thalamic radiation, and brainstem. The extent of NDI and ODI reduction expanded from NHE to MHE. In both MHE and NHE groups, the extent of NDI change was quite larger than that of FA change. No significant intergroup difference in ISO/MD/AD/RD was observed. Tissue specificity afforded by NODDI revealed the underpinning of FA reduction in MHE. The NDI in left frontal lobe was significantly correlated with PHES. CONCLUSION: MHE is characterized by diffuse WM microstructural impairment (especially neurite density reduction). NODDI can improve the detection of WM microstructural impairments in MHE and provides more precise information about MHE-related pathology than DTI.

Characterizing Sensorimotor-Related Area Abnormalities in Amyotrophic Lateral Sclerosis: An Intravoxel Incoherent Motion Magnetic Resonance Imaging Study.

RATIONALE AND OBJECTIVES: To investigate the microperfusion and water molecule diffusion alterations in sensorimotor-related areas in amyotrophic lateral sclerosis (ALS) using intravoxel incoherent motion (IVIM) magnetic resonance imaging. MATERIALS AND METHODS: IVIM data were obtained from 43 ALS patients and 31 controls. This study employed the revised ALS Functional Rating Scale (ALSFRS-R) in evaluating disease severity. IVIM-derived metrics were calculated, including diffusion coefficient (D), pseudo-diffusion coefficient, and perfusion fraction. Conventional apparent diffusion coefficient was also computed. Atlas-based analysis was conducted to detect between-group difference in these metrics in sensorimotor-related gray/white matter areas. Spearman correlation analysis was employed to establish correlation between various metrics and ALSFRS-R. RESULTS: ALS patients had perfusion fraction (× 10-3) reduction in the left presupplementary motor area (60.72 ± 16.15 vs. 71.15 ± 12.98, p = 0.016), right presupplementary motor area (61.35 ± 17.02 vs. 72.18 ± 14.22, p = 0.016), left supplementary motor area (55.73 ± 12.29 vs. 64.12 ± 9.17, p = 0.015), and right supplementary motor area (56.53 ± 11.93 vs. 63.67 ± 10.03, p = 0.020). Patients showed D (× 10-6 mm2/s) increase in a white matter tract projecting to the right ventral premotor cortex (714.20 ± 39.75 vs. 691.01 ± 24.53, p = 0.034). A negative correlation between D of right ventral premotor cortex tract and ALSFRS-R score was observed (r = -0.316, p = 0.039). CONCLUSION: These findings suggest aberrant microperfusion and water molecule diffusion in the sensorimotor-related areas in ALS patients, which are associated with motor impairment in ALS.

White matter microstructural impairments in amyotrophic lateral sclerosis: A mean apparent propagator MRI study.

BACKGROUND: White matter (WM) impairment is a hallmark of amyotrophic lateral sclerosis (ALS). This study evaluated the capacity of mean apparent propagator magnetic resonance imaging (MAP-MRI) for detecting ALS-related WM alterations. METHODS: Diffusion images were obtained from 52 ALS patients and 51 controls. MAP-derived indices [return-to-origin/-axis/-plane probability (RTOP/RTAP/RTPP) and non-Gaussianity (NG)/perpendicular/parallel NG (NG⊥/NG||)] were computed. Measures from diffusion tensor/kurtosis imaging (DTI/DKI) and neurite orientation dispersion and density imaging (NODDI) were also obtained. Voxel-wise analysis (VBA) was performed to determine differences in these parameters. Relationship between MAP parameters and disease severity (assessed by the revised ALS Functional Rating Scale (ALSFRS-R)) was evaluated by Pearson's correlation analysis in a voxel-wise way. ALS patients were further divided into two subgroups: 29 with limb-only involvement and 23 with both bulbar and limb involvement. Subgroup analysis was then conducted to investigate diffusion parameter differences related to bulbar impairment. RESULTS: The VBA (with threshold of P < 0.05 after family-wise error correction (FWE)) showed that ALS patients had significantly decreased RTOP/RTAP/RTPP and NG/ NG⊥/NG|| in a set of WM areas, including the bilateral precentral gyrus, corona radiata, posterior limb of internal capsule, midbrain, middle corpus callosum, anterior corpus callosum, parahippocampal gyrus, and medulla. MAP-MRI had the capacity to capture WM damage in ALS, which was higher than DTI and similar to DKI/NODDI. RTOP/RTAP/NG/NG⊥/NG|| parameters, especially in the bilateral posterior limb of internal capsule and middle corpus callosum, were significantly correlated with ALSFRS-R (with threshold of FWE-corrected P < 0.05). The VBA (with FWE-corrected P < 0.05) revealed the significant RTAP reduction in subgroup with both bulbar and limb involvement, compared with those with limb-only involvement. CONCLUSIONS: Microstructural impairments in corticospinal tract and corpus callosum represent the consistent characteristic of ALS. MAP-MRI could provide alternative measures depicting ALS-related WM alterations, complementary to the common diffusion imaging methods.

Mean Apparent Propagator MRI Is Better Than Conventional Diffusion Tensor Imaging for the Evaluation of Parkinson's Disease: A Prospective Pilot Study.

BACKGROUND AND PURPOSE: Mean apparent propagator (MAP) MRI is a novel diffusion imaging method to map tissue microstructure. The purpose of this study was to evaluate the diagnostic value of the MAP MRI in Parkinson's disease (PD) in comparison with conventional diffusion tensor imaging (DTI). METHODS: 23 PD patients and 22 age- and gender-matched healthy controls were included. MAP MRI and DTI were performed on a 3T MR scanner with a 20-channel head coil. The MAP metrics including mean square displacement (MSD), return to the origin probability (RTOP), return to the axis probability (RTAP), and return to the plane probability (RTPP), and DTI metrics including fractional anisotropy (FA), and mean diffusivity (MD), were measured in subcortical gray matter and compared between the two groups. The receiver operating characteristic (ROC) curve was used to analyze the diagnostic performance of all the metrics. The association between the diffusion metrics and disease severity was assessed by Pearson correlation analysis. RESULTS: For MAP MRI, the mean values of MSD in the bilateral caudate, pallidum, putamen, thalamus and substantia nigra (SN) were higher in PD patients than in healthy controls (p FDR ≤ 0.001); the mean values of the zero displacement probabilities (RTOP, RTAP, and RTPP) in the bilateral caudate, pallidum, putamen and thalamus were lower in PD patients (p FDR < 0.001). For DTI, only FA in the bilateral SN was significantly higher in PD patients than those in the controls (p FDR < 0.001). ROC analysis showed that the areas under the curves of MAP MRI metrics (MSD, RTOP, RTAP, and RTPP) in the bilateral caudate, pallidum, putamen and thalamus (range, 0.85-0.94) were greater than those of FA and MD of DTI (range, 0.55-0.69) in discriminating between PD patients and healthy controls. RTAP in the ipsilateral pallidum (r = -0.56, p FDR = 0.027), RTOP in the bilateral and contralateral putamen (r = -0.58, p FDR = 0.019; r = -0.57, p FDR = 0.024) were negatively correlated with UPDRS III motor scores. CONCLUSION: MAP MRI outperformed the conventional DTI in the diagnosis of PD and evaluation of the disease severity.

Differentiation between high-grade gliomas and solitary brain metastases: a comparison of five diffusion-weighted MRI models.

BACKGROUND: To compare the diagnostic performance of neurite orientation dispersion and density imaging (NODDI), mean apparent propagator magnetic resonance imaging (MAP-MRI), diffusion kurtosis imaging (DKI), diffusion tensor imaging (DTI) and diffusion-weighted imaging (DWI) in distinguishing high-grade gliomas (HGGs) from solitary brain metastases (SBMs). METHODS: Patients with previously untreated, histopathologically confirmed HGGs (n = 20) or SBMs (n = 21) appearing as a solitary and contrast-enhancing lesion on structural MRI were prospectively recruited to undergo diffusion-weighted MRI. DWI data were obtained using a q-space Cartesian grid sampling procedure and were processed to generate parametric maps by fitting the NODDI, MAP-MRI, DKI, DTI and DWI models. The diffusion metrics of the contrast-enhancing tumor and peritumoral edema were measured. Differences in the diffusion metrics were compared between HGGs and SBMs, followed by receiver operating characteristic (ROC) analysis and the Hanley and McNeill test to determine their diagnostic performances. RESULTS: NODDI-based isotropic volume fraction (Viso) and orientation dispersion index (ODI); MAP-MRI-based mean-squared displacement (MSD) and q-space inverse variance (QIV); DKI-generated radial, mean diffusivity and fractional anisotropy (RDk, MDk and FAk); and DTI-generated radial, mean diffusivity and fractional anisotropy (RD, MD and FA) of the contrast-enhancing tumor were significantly different between HGGs and SBMs (p < 0.05). The best single discriminative parameters of each model were Viso, MSD, RDk and RD for NODDI, MAP-MRI, DKI and DTI, respectively. The AUC of Viso (0.871) was significantly higher than that of MSD (0.736), RDk (0.760) and RD (0.733) (p < 0.05). CONCLUSION: NODDI outperforms MAP-MRI, DKI, DTI and DWI in differentiating between HGGs and SBMs. NODDI-based Viso has the highest performance.

FeAture Explorer (FAE): A tool for developing and comparing radiomics models.

In radiomics studies, researchers usually need to develop a supervised machine learning model to map image features onto the clinical conclusion. A classical machine learning pipeline consists of several steps, including normalization, feature selection, and classification. It is often tedious to find an optimal pipeline with appropriate combinations. We designed an open-source software package named FeAture Explorer (FAE). It was programmed with Python and used NumPy, pandas, and scikit-learning modules. FAE can be used to extract image features, preprocess the feature matrix, develop different models automatically, and evaluate them with common clinical statistics. FAE features a user-friendly graphical user interface that can be used by radiologists and researchers to build many different pipelines, and to compare their results visually. To prove the effectiveness of FAE, we developed a candidate model to classify the clinical-significant prostate cancer (CS PCa) and non-CS PCa using the PROSTATEx dataset. We used FAE to try out different combinations of feature selectors and classifiers, compare the area under the receiver operating characteristic curve of different models on the validation dataset, and evaluate the model using independent test data. The final model with the analysis of variance as the feature selector and linear discriminate analysis as the classifier was selected and evaluated conveniently by FAE. The area under the receiver operating characteristic curve on the training, validation, and test dataset achieved results of 0.838, 0.814, and 0.824, respectively. FAE allows researchers to build radiomics models and evaluate them using an independent testing dataset. It also provides easy model comparison and result visualization. We believe FAE can be a convenient tool for radiomics studies and other medical studies involving supervised machine learning.

Computer-aided diagnosis of prostate cancer using a deep convolutional neural network from multiparametric MRI.

BACKGROUND: Deep learning is the most promising methodology for automatic computer-aided diagnosis of prostate cancer (PCa) with multiparametric MRI (mp-MRI). PURPOSE: To develop an automatic approach based on deep convolutional neural network (DCNN) to classify PCa and noncancerous tissues (NC) with mp-MRI. STUDY TYPE: Retrospective. SUBJECTS: In all, 195 patients with localized PCa were collected from a PROSTATEx database. In total, 159/17/19 patients with 444/48/55 observations (215/23/23 PCas and 229/25/32 NCs) were randomly selected for training/validation/testing, respectively. SEQUENCE: T2 -weighted, diffusion-weighted, and apparent diffusion coefficient images. ASSESSMENT: A radiologist manually labeled the regions of interest of PCas and NCs and estimated the Prostate Imaging Reporting and Data System (PI-RADS) scores for each region. Inspired by VGG-Net, we designed a patch-based DCNN model to distinguish between PCa and NCs based on a combination of mp-MRI data. Additionally, an enhanced prediction method was used to improve the prediction accuracy. The performance of DCNN prediction was tested using a receiver operating characteristic (ROC) curve, and the area under the ROC curve (AUC), sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) were calculated. Moreover, the predicted result was compared with the PI-RADS score to evaluate its clinical value using decision curve analysis. STATISTICAL TEST: Two-sided Wilcoxon signed-rank test with statistical significance set at 0.05. RESULTS: The DCNN produced excellent diagnostic performance in distinguishing between PCa and NC for testing datasets with an AUC of 0.944 (95% confidence interval: 0.876-0.994), sensitivity of 87.0%, specificity of 90.6%, PPV of 87.0%, and NPV of 90.6%. The decision curve analysis revealed that the joint model of PI-RADS and DCNN provided additional net benefits compared with the DCNN model and the PI-RADS scheme. DATA CONCLUSION: The proposed DCNN-based model with enhanced prediction yielded high performance in statistical analysis, suggesting that DCNN could be used in computer-aided diagnosis (CAD) for PCa classification. LEVEL OF EVIDENCE: 3 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2018;48:1570-1577.

A diffusional kurtosis imaging study of idiopathic generalized epilepsy with unilateral interictal epileptiform discharges in children.

BACKGROUND AND PURPOSE: To investigate brain abnormalities in children with a clinical diagnosis of idiopathic generalized epilepsy (IGE) and unilateral interictal epileptiform discharges (IED) demonstrated on electroencephalography (EEG) by diffusional kurtosis imaging (DKI). MATERIALS AND METHODS: DKI images were obtained from 18 patients (n=9 each in the left and right hemispheres). Fractional anisotropy (FA), mean diffusivity (MD), and mean kurtosis (MK) maps were estimated through voxel-based analyses, and compared with 18 normal controls matched for age and sex. RESULTS: In the left side group, the significant differences of FA were in the left fusiform gyrus and occipital lobe of the white matter (WM). The significant differences of MD were in the left pons. The significant differences of MK were in the anterior cingulate gyrus, limbic lobe, gray matter (GM) and WM of the right cerebrum. In the right side group, the significant differences of FA were in the WM of the left cerebrum. MD identified differences in the frontal, temporal, occipital, and parietal lobes of both hemispheres, especially in the limbic system, fusiform gyrus, uncus, and parahippocampal gyrus. The significant differences of MK were in the GM of the right cerebrum, particularly in the rolandic operculum and frontal lobe. CONCLUSIONS: DKI is sensitive for the detection of diffusion abnormalities in both WM and GM of IGE in children. Secondary brain abnormalities may exist in regions outside the unilateral epileptogenic zone through the limbic epileptic network, and can be detected by DKI indices FA, MD and MK.

Evaluation of non-local means based denoising filters for diffusion kurtosis imaging using a new phantom.

Image denoising has a profound impact on the precision of estimated parameters in diffusion kurtosis imaging (DKI). This work first proposes an approach to constructing a DKI phantom that can be used to evaluate the performance of denoising algorithms in regard to their abilities of improving the reliability of DKI parameter estimation. The phantom was constructed from a real DKI dataset of a human brain, and the pipeline used to construct the phantom consists of diffusion-weighted (DW) image filtering, diffusion and kurtosis tensor regularization, and DW image reconstruction. The phantom preserves the image structure while minimizing image noise, and thus can be used as ground truth in the evaluation. Second, we used the phantom to evaluate three representative algorithms of non-local means (NLM). Results showed that one scheme of vector-based NLM, which uses DWI data with redundant information acquired at different b-values, produced the most reliable estimation of DKI parameters in terms of Mean Square Error (MSE), Bias and standard deviation (Std). The result of the comparison based on the phantom was consistent with those based on real datasets.

A fast schema for parameter estimation in diffusion kurtosis imaging.

Diffusion kurtosis imaging (DKI) is a new model in magnetic resonance imaging (MRI) characterizing restricted diffusion of water molecules in living tissues. We propose a method for fast estimation of the DKI parameters. These parameters - apparent diffusion coefficient (ADC) and apparent kurtosis coefficient (AKC) - are evaluated using an alternative iteration schema (AIS). This schema first roughly estimates a pair of ADC and AKC values from a subset of the DKI data acquired at 3 b-values. It then iteratively and alternately updates the ADC and AKC until they are converged. This approach employs the technique of linear least square fitting to minimize estimation error in each iteration. In addition to the common physical and biological constrains that set the upper and lower boundaries of the ADC and AKC values, we use a smoothing procedure to ensure that estimation is robust. Quantitative comparisons between our AIS methods and the conventional methods of unconstrained nonlinear least square (UNLS) using both synthetic and real data showed that our unconstrained AIS method can significantly accelerate the estimation procedure without compromising its accuracy, with the computational time for a DKI dataset successfully reduced to only 1 or 2min. Moreover, the incorporation of the smoothing procedure using one of our AIS methods can significantly enhance the contrast of AKC maps and greatly improve the visibility of details in fine structures.

Evaluation of optimized b-value sampling schemas for diffusion kurtosis imaging with an application to stroke patient data.

Diffusion kurtosis imaging (DKI) is a new method of magnetic resonance imaging (MRI) that provides non-Gaussian information that is not available in conventional diffusion tensor imaging (DTI). DKI requires data acquisition at multiple b-values for parameter estimation; this process is usually time-consuming. Therefore, fewer b-values are preferable to expedite acquisition. In this study, we carefully evaluated various acquisition schemas using different numbers and combinations of b-values. Acquisition schemas that sampled b-values that were distributed to two ends were optimized. Compared to conventional schemas using equally spaced b-values (ESB), optimized schemas require fewer b-values to minimize fitting errors in parameter estimation and may thus significantly reduce scanning time. Following a ranked list of optimized schemas resulted from the evaluation, we recommend the 3b schema based on its estimation accuracy and time efficiency, which needs data from only 3 b-values at 0, around 800 and around 2600 s/mm2, respectively. Analyses using voxel-based analysis (VBA) and region-of-interest (ROI) analysis with human DKI datasets support the use of the optimized 3b (0, 1000, 2500 s/mm2) DKI schema in practical clinical applications.