Diffusion MRI head motion correction methods are highly accurate but impacted by denoising and sampling scheme

Head motion correction is particularly challenging in diffusion-weighted MRI (dMRI) scans due to the dramatic changes in image contrast at different gradient strengths and directions. Head motion correction is typically performed using a Gaussian Process model implemented in FSL's Eddy. Recently, the 3dSHORE-based SHORELine method was introduced that does not require shell-based acquisitions, but it has not been previously benchmarked. Here we perform a comprehensive evaluation of both methods on realistic simulations of a software fiber phantom that provides known ground-truth head motion. We demonstrate that both methods perform remarkably well, but that performance can be impacted by sampling scheme and the extent of head motion and the denoising strategy applied before head motion correction. Furthermore, we find Eddy benefits from denoising the data first with MP-PCA. In sum, we provide the most extensive known benchmarking of dMRI head motion correction, together with extensive simulation data and a reproducible workflow. PRACTITIONER POINTS: Both Eddy and SHORELine head motion correction methods performed quite well on a large variety of simulated data. Denoising with MP-PCA can improve head motion correction performance when Eddy is used. SHORELine effectively corrects motion in non-shelled diffusion spectrum imaging data.

Applications of advanced diffusion MRI in early brain development: a comprehensive review

Brain development follows a protracted developmental timeline with foundational processes of neurodevelopment occurring from the third trimester of gestation into the first decade of life. Defining structural maturational patterns of early brain development is a critical step in detecting divergent developmental trajectories associated with neurodevelopmental and psychiatric disorders that arise later in life. While considerable advancements have already been made in diffusion magnetic resonance imaging (dMRI) for pediatric research over the past three decades, the field of neurodevelopment is still in its infancy with remarkable scientific and clinical potential. This comprehensive review evaluates the application, findings, and limitations of advanced dMRI methods beyond diffusion tensor imaging, including diffusion kurtosis imaging (DKI), constrained spherical deconvolution (CSD), neurite orientation dispersion and density imaging (NODDI) and composite hindered and restricted model of diffusion (CHARMED) to quantify the rapid and dynamic changes supporting the underlying microstructural architectural foundations of the brain in early life.

Test-retest reliability of diffusion tensor imaging scalars in 5-year-olds

Diffusion tensor imaging (DTI) has provided great insights into the microstructural features of the developing brain. However, DTI images are prone to several artifacts and the reliability of DTI scalars is of paramount importance for interpreting and generalizing the findings of DTI studies, especially in the younger population. In this study, we investigated the intrascan test–retest repeatability of four DTI scalars: fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) in 5‐year‐old children (N = 67) with two different data preprocessing approaches: a volume censoring pipeline and an outlier replacement pipeline. We applied a region of interest (ROI) and a voxelwise analysis after careful quality control, tensor fitting and tract‐based spatial statistics. The data had three subsets and each subset included 31, 32, or 33 directions thus a total of 96 unique uniformly distributed diffusion encoding directions per subject. The repeatability of DTI scalars was evaluated with intraclass correlation coefficient (ICC(3,1)) and the variability between test and retest subsets. The results of both pipelines yielded good to excellent (ICC(3,1) > 0.75) reliability for most of the ROIs and an overall low variability (<10%). In the voxelwise analysis, FA and RD had higher ICC(3,1) values compared to AD and MD and the variability remained low (<12%) across all scalars. Our results suggest high intrascan repeatability in pediatric DTI and lend confidence to the use of the data in future cross‐sectional and longitudinal studies.

A Comparative Multimodal Meta-analysis of Anisotropy and Volume Abnormalities in White Matter in People Suffering From Bipolar Disorder or Schizophrenia

Schizophrenia (SZ) and bipolar disorder (BD) share some similarities in terms of genetic-risk genes and abnormalities of gray-matter structure in the brain, but white matter (WM) abnormalities have not been studied in depth. We undertook a comparative multimodal meta-analysis to identify common and disorder-specific abnormalities in WM structure between SZ and BD. Anisotropic effect size-signed differential mapping software was used to conduct a comparative meta-analysis of 68 diffusion tensor imaging (DTI) and 34 voxel-based morphometry (VBM) studies comparing fractional anisotropy (FA) and white matter volume (WMV), respectively, between patients with SZ (DTI: N = 1543; VBM: N = 1068) and BD (DTI: N = 983; VBM: N = 518) and healthy controls (HCs). The bilateral corpus callosum (extending to the anterior and superior corona radiata) showed shared decreased WMV and FA in SZ and BD. Compared with BD patients, SZ patients showed remarkable disorder-specific WM abnormalities: decreased FA and increased WMV in the left cingulum, and increased FA plus decreased WMV in the right anterior limb of the internal capsule. SZ patients showed more extensive alterations in WM than BD cases, which may be the pathophysiological basis for the clinical continuity of both disorders. The disorder-specific regions in the left cingulum and right anterior limb of the internal capsule provided novel insights into both disorders. Our study adds value to further understanding of the pathophysiology, classification, and differential diagnosis of SZ and BD.

Feasibility of diffusion and probabilistic white matter analysis in patients implanted with a deep brain stimulator

Deep brain stimulation (DBS) for Parkinson's disease (PD) is an established advanced therapy that produces therapeutic effects through high frequency stimulation. Although this therapeutic option leads to improved clinical outcomes, the mechanisms of the underlying efficacy of this treatment are not well understood. Therefore, investigation of DBS and its postoperative effects on brain architecture is of great interest. Diffusion weighted imaging (DWI) is an advanced imaging technique, which has the ability to estimate the structure of white matter fibers; however, clinical application of DWI after DBS implantation is challenging due to the strong susceptibility artifacts caused by implanted devices. This study aims to evaluate the feasibility of generating meaningful white matter reconstructions after DBS implantation; and to subsequently quantify the degree to which these tracts are affected by post-operative device-related artifacts. DWI was safely performed before and after implanting electrodes for DBS in 9 PD patients. Differences within each subject between pre- and post-implantation FA, MD, and RD values for 123 regions of interest (ROIs) were calculated. While differences were noted globally, they were larger in regions directly affected by the artifact. White matter tracts were generated from each ROI with probabilistic tractography, revealing significant differences in the reconstruction of several white matter structures after DBS. Tracts pertinent to PD, such as regions of the substantia nigra and nigrostriatal tracts, were largely unaffected. The aim of this study was to demonstrate the feasibility and clinical applicability of acquiring and processing DWI post-operatively in PD patients after DBS implantation. The presence of global differences provides an impetus for acquiring DWI shortly after implantation to establish a new baseline against which longitudinal changes in brain connectivity in DBS patients can be compared. Understanding that post-operative fiber tracking in patients is feasible on a clinically-relevant scale has significant implications for increasing our current understanding of the pathophysiology of movement disorders, and may provide insights into better defining the pathophysiology and therapeutic effects of DBS.

Comparison of manual and automated fiber quantification tractography in patients with temporal lobe epilepsy

•

Tractography approaches showed moderate to good agreement for tract morphology.

•

Along- and whole-tract diffusivity was significantly correlated across approaches.

•

Whole-tract AFQ but not manual tract diffusivity correlated with clinical variables.

•

Absence of excellent agreement between approaches warrants caution.

Objective

To investigate the agreement between manually and automatically generated tracts from diffusion tensor imaging (DTI) in patients with temporal lobe epilepsy (TLE). Whole and along-the-tract diffusivity metrics and correlations with patient clinical characteristics were analyzed with respect to tractography approach.

Methods

We recruited 40 healthy controls and 24 patients with TLE who underwent conventional T1-weighted imaging and 60-direction DTI. An automated (Automated Fiber Quantification, AFQ) and manual (TrackVis) deterministic tractography approach was used to identify the uncinate fasciculus (UF) and parahippocampal white matter bundle (PHWM). Tract diffusion scalar metrics were analyzed with respect to agreement across automated and manual approaches (Dice Coefficient and Spearman correlations), to side of onset of epilepsy and patient clinical characteristics, including duration of epilepsy, age of onset and presence of hippocampal sclerosis.

Results

Across approaches the analysis of tract morphology similarity revealed Dice coefficients at moderate to good agreement (0.54 - 0.6) and significant correlations between diffusion values (Spearman's Rho=0.4–0.9). However, within bilateral PHWM, AFQ yielded significantly lower FA (left: Z = 4.4, p<0.001; right: Z = 5.1, p<0.001) and higher MD values (left: Z=-4.7, p<0.001; right: Z=-3.7, p<0.001) compared to the manual approach. Whole tract DTI metrics determined using AFQ were significantly correlated with patient characteristics, including age of epilepsy onset in FA (R = 0.6, p = 0.02) and MD of the ipsilateral PHWM (R=-0.6, p = 0.02), while duration of epilepsy corrected for age correlated with MD in ipsilateral PHWM (R = 0.7, p<0.01). Correlations between clinical metrics and diffusion values extracted using the manual whole tract technique did not survive correction for multiple comparisons. Both manual and automated along-the-tract analyses demonstrated significant correlations with patient clinical characteristics such as age of onset and epilepsy duration. The strongest and most widespread localized ipsi- and contralateral diffusivity alterations were observed in patients with left TLE and patients with HS compared to controls, while patients with right TLE and patients without HS did not show these strong effects.

Conclusions

Manual and AFQ tractography approaches revealed significant correlations in the reconstruction of tract morphology and extracted whole and along-tract diffusivity values. However, as non-identical methods they differed in the respective yield of significant results across clinical correlations and group-wise statistics. Given the absence of excellent agreement between manual and AFQ techniques as demonstrated in the present study, caution should be considered when using AFQ particularly when used without reference to benchmark manual measures.

Minimisation of Signal Intensity Differences in Distortion Correction Approaches of Brain Magnetic Resonance Diffusion Tensor Imaging

OBJECTIVES

To evaluate the effects of signal intensity differences between the b0 image and diffusion tensor imaging (DTI) in the image registration process.

METHODS

To correct signal intensity differences between the b0 image and DTI data, a simple image intensity compensation (SIMIC) method, which is a b0 image re-calculation process from DTI data, was applied before the image registration. The re-calculated b0 image (b0_ext) from each diffusion direction was registered to the b0 image acquired through the MR scanning (b0_nd) with two types of cost functions and their transformation matrices were acquired. These transformation matrices were then used to register the DTI data. For quantifications, the dice similarity coefficient (DSC) values, diffusion scalar matrix, and quantified fibre numbers and lengths were calculated.

RESULTS

The combined SIMIC method with two cost functions showed the highest DSC value (0.802 ± 0.007). Regarding diffusion scalar values and numbers and lengths of fibres from the corpus callosum, superior longitudinal fasciculus, and cortico-spinal tract, only using normalised cross correlation (NCC) showed a specific tendency toward lower values in the brain regions.

CONCLUSION

Image-based distortion correction with SIMIC for DTI data would help in image analysis by accounting for signal intensity differences as one additional option for DTI analysis.

KEY POINTS

• We evaluated the effects of signal intensity differences at DTI registration. • The non-diffusion-weighted image re-calculation process from DTI data was applied. • SIMIC can minimise the signal intensity differences at DTI registration.

Temporal slice registration and robust diffusion-tensor reconstruction for improved fetal brain structural connectivity analysis

Diffusion weighted magnetic resonance imaging, or DWI, is one of the most promising tools for the analysis of neural microstructure and the structural connectome of the human brain. The application of DWI to map early development of the human connectome in-utero, however, is challenged by intermittent fetal and maternal motion that disrupts the spatial correspondence of data acquired in the relatively long DWI acquisitions. Fetuses move continuously during DWI scans. Reliable and accurate analysis of the fetal brain structural connectome requires careful compensation of motion effects and robust reconstruction to avoid introducing bias based on the degree of fetal motion. In this paper we introduce a novel robust algorithm to reconstruct in-vivo diffusion-tensor MRI (DTI) of the moving fetal brain and show its effect on structural connectivity analysis. The proposed algorithm involves multiple steps of image registration incorporating a dynamic registration-based motion tracking algorithm to restore the spatial correspondence of DWI data at the slice level and reconstruct DTI of the fetal brain in the standard (atlas) coordinate space. A weighted linear least squares approach is adapted to remove the effect of intra-slice motion and reconstruct DTI from motion-corrected data. The proposed algorithm was tested on data obtained from 21 healthy fetuses scanned in-utero at 22-38 weeks gestation. Significantly higher fractional anisotropy values in fiber-rich regions, and the analysis of whole-brain tractography and group structural connectivity, showed the efficacy of the proposed method compared to the analyses based on original data and previously proposed methods. The results of this study show that slice-level motion correction and robust reconstruction is necessary for reliable in-vivo structural connectivity analysis of the fetal brain. Connectivity analysis based on graph theoretic measures show high degree of modularity and clustering, and short average characteristic path lengths indicative of small-worldness property of the fetal brain network. These findings comply with previous findings in newborns and a recent study on fetuses. The proposed algorithm can provide valuable information from DWI of the fetal brain not available in the assessment of the original 2D slices and may be used to more reliably study the developing fetal brain connectome.

Assessing the performance of different DTI motion correction strategies in the presence of EPI distortion correction

Diffusion tensor imaging (DTI) is susceptible to several artifacts due to eddy currents, echo planar imaging (EPI) distortion and subject motion. While several techniques correct for individual distortion effects, no optimal combination of DTI acquisition and processing has been determined. Here, the effects of several motion correction techniques are investigated while also correcting for EPI distortion: prospective correction, using navigation; retrospective correction, using two different popular packages (FSL and TORTOISE); and the combination of both methods. Data from a pediatric group that exhibited incidental motion in varying degrees are analyzed. Comparisons are carried while implementing eddy current and EPI distortion correction. DTI parameter distributions, white matter (WM) maps and probabilistic tractography are examined. The importance of prospective correction during data acquisition is demonstrated. In contrast to some previous studies, results also show that the inclusion of retrospective processing also improved ellipsoid fits and both the sensitivity and specificity of group tractographic results, even for navigated data. Matches with anatomical WM maps are highest throughout the brain for data that have been both navigated and processed using TORTOISE. The inclusion of both prospective and retrospective motion correction with EPI distortion correction is important for DTI analysis, particularly when studying subject populations that are prone to motion. Hum Brain Mapp 37:4405-4424, 2016. © 2016 Wiley Periodicals, Inc.

Harmonizing DTI measurements across scanners to examine the development of white matter microstructure in 803 adolescents of the NCANDA study

Neurodevelopment continues through adolescence, with notable maturation of white matter tracts comprising regional fiber systems progressing at different rates. To identify factors that could contribute to regional differences in white matter microstructure development, large samples of youth spanning adolescence to young adulthood are essential to parse these factors. Recruitment of adequate samples generally relies on multi-site consortia but comes with the challenge of merging data acquired on different platforms. In the current study, diffusion tensor imaging (DTI) data were acquired on GE and Siemens systems through the National Consortium on Alcohol and NeuroDevelopment in Adolescence (NCANDA), a multi-site study designed to track the trajectories of regional brain development during a time of high risk for initiating alcohol consumption. This cross-sectional analysis reports baseline Tract-Based Spatial Statistic (TBSS) of regional fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (L1), and radial diffusivity (LT) from the five consortium sites on 671 adolescents who met no/low alcohol or drug consumption criteria and 132 adolescents with a history of exceeding consumption criteria. Harmonization of DTI metrics across manufacturers entailed the use of human-phantom data, acquired multiple times on each of three non-NCANDA participants at each site's MR system, to determine a manufacturer-specific correction factor. Application of the correction factor derived from human phantom data measured on MR systems from different manufacturers reduced the standard deviation of the DTI metrics for FA by almost a half, enabling harmonization of data that would have otherwise carried systematic error. Permutation testing supported the hypothesis of higher FA and lower diffusivity measures in older adolescents and indicated that, overall, the FA, MD, and L1 of the boys were higher than those of the girls, suggesting continued microstructural development notable in the boys. The contribution of demographic and clinical differences to DTI metrics was assessed with General Additive Models (GAM) testing for age, sex, and ethnicity differences in regional skeleton mean values. The results supported the primary study hypothesis that FA skeleton mean values in the no/low-drinking group were highest at different ages. When differences in intracranial volume were covaried, FA skeleton mean reached a maximum at younger ages in girls than boys and varied in magnitude with ethnicity. Our results, however, did not support the hypothesis that youth who exceeded exposure criteria would have lower FA or higher diffusivity measures than the no/low-drinking group; detecting the effects of excessive alcohol consumption during adolescence on DTI metrics may require longitudinal study.