Search

Background: Achieving inter-site / inter-scanner reproducibility of diffusion weighted magnetic resonance imaging (DW-MRI) metrics has been challenging given differences in acquisition protocols, analysis models, and hardware factors., Purpose: Magnetic field gradients impart scanner-dependent spatial variations in the applied diffusion weighting that can be corrected if the gradient nonlinearities are known. However, retrieving manufacturer nonlinearity specifications is not well supported and may introduce errors in interpretation of units or coordinate systems. We propose an empirical approach to mapping the gradient nonlinearities with sequences that are supported across the major scanner vendors., Study Type: Prospective observational study., Subjects: A spherical isotropic diffusion phantom, and a single human control volunteer., Field Strength/sequence: 3 T (two scanners). Stejskal-Tanner spin echo sequence with b-values of 1000, 2000 s/mm 2 with 12, 32, and 384 diffusion gradient directions per shell., Assessment: We compare the proposed correction with the prior approach using manufacturer specifications against typical diffusion pre-processing pipelines (i.e., ignoring spatial gradient nonlinearities). In phantom data, we evaluate metrics against the ground truth. In human and phantom data, we evaluate reproducibility across scans, sessions, and hardware., Statistical Tests: Wilcoxon rank-sum test between uncorrected and corrected data., Results: In phantom data, our correction method reduces variation in mean diffusivity across sessions over uncorrected data (p < 0.05). In human data, we show that this method can also reduce variation in mean diffusivity across scanners (p < 0.05)., Conclusion: Our method is relatively simple, fast, and can be applied retroactively. We advocate incorporating voxel-specific b-value and b-vector maps should be incorporated in DW-MRI harmonization preprocessing pipelines to improve quantitative accuracy of measured diffusion parameters., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Purpose: Diffusion-weighted magnetic resonance imaging (DW-MRI) is of critical importance for characterizing in-vivo white matter. Models relating microarchitecture to observed DW-MRI signals as a function of diffusion sensitization are the lens through which DW-MRI data are interpreted. Numerous modern approaches offer opportunities to assess more complex intra-voxel structures. Nevertheless, there remains a substantial gap between intra-voxel estimated structures and ground truth captured by 3-D histology., Methods: Herein, we propose a novel data-driven approach to model the non-linear mapping between observed DW-MRI signals and ground truth structures using a sequential deep neural network regression using residual block deep neural network (ResDNN). Training was performed on two 3-D histology datasets of squirrel monkey brains and validated on a third. A second validation was performed using scan-rescan datasets of 12 subjects from Human Connectome Project. The ResDNN was compared with multiple micro-structure reconstruction methods and super resolved-constrained spherical deconvolution (sCSD) in particular as baseline for both the validations., Results: Angular correlation coefficient (ACC) is a correlation/similarity measure and can be interpreted as accuracy when compared with a ground truth. The median ACC of ResDNN is 0.82 and median ACC's of different variants of CSD are 0.75, 0.77, 0.79. The mean, median and std. of ResDNN & sCSD ACC across 12 subjects from HCP are 0.74, 0.88, 0.31 and 0.61, 0.71, 0.31 respectively., Conclusion: This work highlights the ability of deep learning to capture linkages between ex-vivo ground truth data with feasible MRI sequences. The data-driven approach is applicable to human in-vivo data and results in intriguingly high reproducibility of orientation structure., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Neuroimaging often involves acquiring high-resolution anatomical images along with other low-resolution image modalities, like diffusion and functional magnetic resonance imaging. Performing gray matter statistics with low-resolution image modalities is a challenge due to registration artifacts and partial volume effects. Gray matter surface based spatial statistics (GS-BSS) has been shown to provide higher sensitivity using gray matter surfaces compared to that of skeletonization approach of gray matter based spatial statistics which is adapted from tract based spatial statistics in diffusion studies. In this study, we improve upon GS-BSS incorporating neurite orientation dispersion and density imaging (NODDI) based search (denoted N-GSBSS) by 1) enhancing metrics mapping from native space, 2) incorporating maximum orientation dispersion index (ODI) search along surface normal, and 3) proposing applicability to other modalities, such as functional MRI (fMRI). We evaluated the performance of N-GSBSS against three baseline pipelines: volume-based registration, FreeSurfer's surface registration and ciftify pipeline for fMRI and simulation studies. First, qualitative mean ODI results are shown for N-GSBSS with and without NODDI based search in comparison with ciftify pipeline. Second, we conducted one-sample t-tests on working memory activations in fMRI to show that the proposed method can aid in the analysis of low resolution fMRI data. Finally we performed a sensitivity test in a simulation study by varying percentage change of intensity values within a region of interest in gray matter probability maps. N-GSBSS showed higher sensitivity in the simulation test compared to the other methods capturing difference between the groups starting at 10% change in the intensity values. The computational time of N-GSBSS is 68 times faster than that of traditional surface-based or 86 times faster than that of ciftify pipeline analysis., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The ability to evaluate empirical diffusion MRI acquisitions for quality and to correct the resulting imaging metrics allows for improved inference and increased replicability. Previous work has shown promise for estimation of bias and variance of generalized fractional anisotropy (GFA) but comes at the price of computational complexity. This paper aims to provide methods for estimating GFA, bias of GFA and standard deviation of GFA quickly and accurately. In order to provide a method for bias and variance estimation that can return results faster than the previously studied statistical techniques, three deep, fully-connected neural networks are developed for GFA, bias of GFA, and standard deviation of GFA. The results of these networks are compared to the observed values of the metrics as well as those fit from the statistical techniques (i.e. Simulation Extrapolation (SIMEX) for bias estimation and wild bootstrap for variance estimation). Our GFA network provides predictions that are closer to the true GFA values than a Q-ball fit of the observed data (root-mean-square error (RMSE) 0.0077 vs 0.0082, p < .001). The bias network also shows statistically significant improvement in comparison to the SIMEX-estimated error of GFA (RMSE 0.0071 vs. 0.01, p < .001)., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Diffusion weighted MRI (DWMRI) and the myriad of analysis approaches (from tensors to spherical harmonics and brain tractography to body multi-compartment models) depend on accurate quantification of the apparent diffusion coefficient (ADC). Signal drift during imaging (e.g., due to b0 drift associated with heating) can cause systematic non-linearities that manifest as ADC changes if not corrected. Herein, we present a case study on two phantoms on one scanner. Different scan protocols exhibit different degrees of drift during similar scans and may be sensitive to the order of scans within an exam. Vos et al. recently reviewed the effects of signal drift in DWMRI acquisitions and proposed a temporal model for correction. We propose a novel spatial-temporal model to correct for higher order aspects of the signal drift and derive a statistically robust variant. We evaluate the Vos model and propose a method using two phantoms that mimic the ADC of the relevant brain tissue (0.36-2.2 × 10-3 mm 2 /s) on a single 3 T scanner. The phantoms are (1) a spherical isotropic sphere consisting of a single concentration of polyvinylpyrrolidone (PVP) and (2) an ice-water phantom with 13 vials of varying PVP concentrations. To characterize the impact of interspersed minimally weighted volumes ("b0's"), image volumes with b-value equal to 0.1 s/mm 2 are interspersed every 8, 16, 32, 48, and 96 diffusion weighted volumes in different trials. Signal drift is found to have spatially varying effects that are not accounted for with temporal-only models. The novel model captures drift more accurately (i.e., reduces the overall change per-voxel over the course of a scan) and results in more consistent ADC metrics., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Neurite Orientation Dispersion and Density Imaging (NODDI) has been gaining prominence for estimating multiple diffusion compartments from MRI data acquired in a clinically feasible time. To establish a pathway for adoption of NODDI in clinical studies, it is important to understand the sensitivity and reproducibility of NODDI metrics on empirical data in the context of acquisition protocol and brain anatomy. Previous studies addressed reproducibility across the 3 T scanners and within session and between subject reproducibility at 1.5 T and 3 T. However, empirical reproducibility on the performance of NODDI metrics based on b-value and diffusion-sensitized directions has not yet been addressed. In this study, we investigate a high angular resolution dataset with 11 repeats of a study with five b-values shells (1000, 1500, 2000, 2500 and 3000 s/mm 2 ) and 96 directions per shell on a single subject. We validated the findings with a dataset from second subject with 10 repeats and 3 b-value shells (1000, 2000, 3000 s/mm 2 ). The NODDI model was estimated using Accelerated Microstructure Imaging via Convex Optimization (AMICO) for different b-values and gradient directions on two-shell High Angular Resolution Density Imaging (HARDI) data fixing the lower shell at b = 1000 s/mm 2 . NODDI model applied to all acquired imaging data was used as a baseline gold standard for comparison. Additionally, we characterize orientation dispersion index (ODI) reproducibility using single-shell data. The experimental findings confirmed the sensitivity of intracellular volume fraction (V ic ) with the choice of outer shell b-value more than with the choice of gradient directions. On the other hand, ODI is more sensitive to the number of gradient directions compared to b-value selection. Single-shell results for ODI are more comparable to 2-shell data at lower b-values than higher b-values. Recommended settings by region of interest and acquisition time are reported for the researchers considering using NODDI in human studies and/or comparing results across acquisition protocols., (Copyright © 2018 Elsevier Inc. All rights reserved.)