Reducing susceptibility distortion related image blurring in diffusion MRI EPI data

Diffusion magnetic resonance imaging (MRI) is an increasingly popular technique in basic and clinical neuroscience. One promising application is to combine diffusion MRI with myelin maps from complementary MRI techniques such as multi-parameter mapping (MPM) to produce g-ratio maps that represent the relative myelination of axons and predict their conduction velocity. Statistical Parametric Mapping (SPM) can process both diffusion data and MPMs, making SPM the only widely accessible software that contains all the processing steps required to perform group analyses of g-ratio data in a common space. However, limitations have been identified in its method for reducing susceptibility-related distortion in diffusion data. More generally, susceptibility-related image distortion is often corrected by combining reverse phase-encoded images (blip-up and blip-down) using the arithmetic mean (AM), however, this can lead to blurred images. In this study we sought to (1) improve the susceptibility-related distortion correction for diffusion MRI data in SPM; (2) deploy an alternative approach to the AM to reduce image blurring in diffusion MRI data when combining blip-up and blip-down EPI data after susceptibility-related distortion correction; and (3) assess the benefits of these changes for g-ratio mapping. We found that the new processing pipeline, called consecutive Hyperelastic Susceptibility Artefact Correction (HySCO) improved distortion correction when compared to the standard approach in the ACID toolbox for SPM. Moreover, using a weighted average (WA) method to combine the distortion corrected data from each phase-encoding polarity achieved greater overlap of diffusion and more anatomically faithful structural white matter probability maps derived from minimally distorted multi-parameter maps as compared to the AM. Third, we showed that the consecutive HySCO WA performed better than the AM method when combined with multi-parameter maps to perform g-ratio mapping. These improvements mean that researchers can conveniently access a wide range of diffusion-related analysis methods within one framework because they are now available within the open-source ACID toolbox as part of SPM, which can be easily combined with other SPM toolboxes, such as the hMRI toolbox, to facilitate computation of myelin biomarkers that are necessary for g-ratio mapping

A protocol for the analysis of DTI data collected from young children

Analysis of scalar maps obtained by diffusion tensor imaging (DTI) produce valuable information about the microstructure of the brain white matter. The DTI scanning of child populations, compared with adult groups, requires specifically designed data acquisition protocols that take into consideration the trade-off between the scanning time, diffusion strength, number of diffusion directions, and the applied analysis techniques. Furthermore, inadequate normalization of DTI images and non-robust tensor reconstruction have profound effects on data analyses and may produce biased statistical results. Here, we present an acquisition sequence that was specifically designed for pediatric populations, and describe the analysis steps of the DTI data collected from extremely preterm-born young school-aged children and their age- and gender-matched controls. The protocol utilizes multiple software packages to address the effects of artifacts and to produce robust tensor estimation. The computation of a population-specific template and the nonlinear registration of tensorial images with this template were implemented to improve alignment of brain images from the children.Peer reviewe

Sampling strategies and integrated reconstruction for reducing distortion and boundary slice aliasing in high-resolution 3D diffusion MRI

Purpose: To develop a new method for high-fidelity, high-resolution 3D multi-slab diffusion MRI with minimal distortion and boundary slice aliasing. Methods: Our method modifies 3D multi-slab imaging to integrate blip-reversed acquisitions for distortion correction and oversampling in the slice direction (kz) for reducing boundary slice aliasing. Our aim is to achieve robust acceleration to keep the scan time the same as conventional 3D multi-slab acquisitions, in which data are acquired with a single direction of blip traversal and without kz-oversampling. We employ a two-stage reconstruction. In the first stage, the blip-up/down images are respectively reconstructed and analyzed to produce a field map for each diffusion direction. In the second stage, the blip-reversed data and the field map are incorporated into a joint reconstruction to produce images that are corrected for distortion and boundary slice aliasing. Results: We conducted experiments at 7T in six healthy subjects. Stage 1 reconstruction produces images from highly under-sampled data (R = 7.2) with sufficient quality to provide accurate field map estimation. Stage 2 joint reconstruction substantially reduces distortion artifacts with comparable quality to fully-sampled blip-reversed results (2.4× scan time). Whole-brain in-vivo results acquired at 1.22 mm and 1.05 mm isotropic resolutions demonstrate improved anatomical fidelity compared to conventional 3D multi-slab imaging. Data demonstrate good reliability and reproducibility of the proposed method over multiple subjects. Conclusion: The proposed acquisition and reconstruction framework provide major reductions in distortion and boundary slice aliasing for 3D multi-slab diffusion MRI without increasing the scan time, which can potentially produce high-quality, high-resolution diffusion MRI

Comparison of distortion correction preprocessing pipelines for DTI in the upper limb

PURPOSE: DTI characterizes tissue microstructure and provides proxy measures of nerve health. Echo-planar imaging is a popular method of acquiring DTI but is susceptible to various artifacts (e.g., susceptibility, motion, and eddy currents), which may be ameliorated via preprocessing. There are many pipelines available but limited data comparing their performance, which provides the rationale for this study. METHODS: DTI was acquired from the upper limb of heathy volunteers at 3T in blip-up and blip-down directions. Data were independently corrected using (i) FSL's TOPUP & eddy, (ii) FSL's TOPUP, (iii) DSI Studio, and (iv) TORTOISE. DTI metrics were extracted from the median, radial, and ulnar nerves and compared (between pipelines) using mixed-effects linear regression. The geometric similarity of corrected b = 0 images and the slice matched T1-weighted (T1w) images were computed using the Sörenson-Dice coefficient. RESULTS: Without preprocessing, the similarity coefficient of the blip-up and blip-down datasets to the T1w was 0·80 and 0·79, respectively. Preprocessing improved the geometric similarity by 1% with no difference between pipelines. Compared to TOPUP & eddy, DSI Studio and TORTOISE generated 2% and 6% lower estimates of fractional anisotropy, and 6% and 13% higher estimates of radial diffusivity, respectively. Estimates of anisotropy from TOPUP & eddy versus TOPUP were not different but TOPUP reduced radial diffusivity by 3%. The agreement of DTI metrics between pipelines was poor. CONCLUSIONS: Preprocessing DTI from the upper limb improves geometric similarity but the choice of the pipeline introduces clinically important variability in diffusion parameter estimates from peripheral nerves

PyHySCO: GPU-enabled susceptibility artifact distortion correction in seconds

Over the past decade, reversed gradient polarity (RGP) methods have become a popular approach for correcting susceptibility artifacts in echo-planar imaging (EPI). Although several post-processing tools for RGP are available, their implementations do not fully leverage recent hardware, algorithmic, and computational advances, leading to correction times of several minutes per image volume. To enable 3D RGP correction in seconds, we introduce PyTorch Hyperelastic Susceptibility Correction (PyHySCO), a user-friendly EPI distortion correction tool implemented in PyTorch that enables multi-threading and efficient use of graphics processing units (GPUs). PyHySCO uses a time-tested physical distortion model and mathematical formulation and is, therefore, reliable without training. An algorithmic improvement in PyHySCO is its use of the one-dimensional distortion correction method by Chang and Fitzpatrick to initialize the non-linear optimization. PyHySCO is published under the GNU public license and can be used from the command line or its Python interface. Our extensive numerical validation using 3T and 7T data from the Human Connectome Project suggests that PyHySCO can achieve accuracy comparable to that of leading RGP tools at a fraction of the cost. We also validate the new initialization scheme, compare different optimization algorithms, and test the algorithm on different hardware and arithmetic precisions

What's new and what's next in diffusion MRI preprocessing

Diffusion MRI (dMRI) provides invaluable information for the study of tissue microstructure and brain connectivity, but suffers from a range of imaging artifacts that greatly challenge the analysis of results and their interpretability if not appropriately accounted for. This review will cover dMRI artifacts and preprocessing steps, some of which have not typically been considered in existing pipelines or reviews, or have only gained attention in recent years: brain/skull extraction, B-matrix incompatibilities w.r.t the imaging data, signal drift, Gibbs ringing, noise distribution bias, denoising, between- and within-volumes motion, eddy currents, outliers, susceptibility distortions, EPI Nyquist ghosts, gradient deviations, bias fields, and spatial normalization. The focus will be on “what’s new” since the notable advances prior to and brought by the Human Connectome Project (HCP), as presented in the predecessing issue on “Mapping the Connectome” in 2013. In addition to the development of novel strategies for dMRI preprocessing, exciting progress has been made in the availability of open source tools and reproducible pipelines, databases and simulation tools for the evaluation of preprocessing steps, and automated quality control frameworks, amongst others. Finally, this review will consider practical considerations and our view on “what’s next” in dMRI preprocessing

Central white matter integrity alterations in 2-3-year-old children following prenatal alcohol exposure.

BACKGROUND: Prenatal alcohol exposure (PAE) remains a potentially preventable, but pervasive risk factor to neurodevelopment. Yet, evidence is lacking on the impact of alcohol on brain development in toddlers. This study aimed to investigate the impact of PAE on brain white matter integrity in 2-3-year-old children. METHODS: Children (n = 83, 30-37 months old) of the Drakenstein Child Health Study birth cohort, underwent diffusion MRI on a 3 T Siemens scanner during natural sleep. Parameters were extracted in children with PAE (n = 25, 56 % boys) and unexposed controls (n = 58, 62 % boys) using Tract-based Spatial Statistics, and compared by group. The contribution of maternal tobacco smoking to white matter differences was also explored. RESULTS: Children with PAE had altered fractional anisotropy, radial diffusivity and axial diffusivity in brain stem, limbic and association tracts compared to unexposed controls. Notably lower fractional anisotropy was found in the uncinate fasciculus, and lower mean and radial diffusivity were found in the fornix stria terminalis and corticospinal tract (FDR corrected p < 0.05). There was a significant interaction effect of PAE and prenatal tobacco exposure which lowered mean, radial and axial diffusivity in the corticospinal tract significantly in the PAE group but not controls. CONCLUSION: Widespread altered white matter microstructural integrity at 2-3 years of age is consistent with findings in neonates in the same and other cohorts, indicating persistence of effects of PAE through early life. Findings also highlight that prenatal tobacco exposure impacts the association of PAE on white matter alterations, amplifying effects in tracts underlying motor function

Single- and Multiple-Shell Uniform Sampling Schemes for Diffusion MRI Using Spherical Codes

In diffusion MRI (dMRI), a good sampling scheme is important for efficient acquisition and robust reconstruction. Diffusion weighted signal is normally acquired on single or multiple shells in q-space. Signal samples are typically distributed uniformly on different shells to make them invariant to the orientation of structures within tissue, or the laboratory coordinate frame. The Electrostatic Energy Minimization (EEM) method, originally proposed for single shell sampling scheme in dMRI, was recently generalized to multi-shell schemes, called Generalized EEM (GEEM). GEEM has been successfully used in the Human Connectome Project (HCP). However, EEM does not directly address the goal of optimal sampling, i.e., achieving large angular separation between sampling points. In this paper, we propose a more natural formulation, called Spherical Code (SC), to directly maximize the minimal angle between different samples in single or multiple shells. We consider not only continuous problems to design single or multiple shell sampling schemes, but also discrete problems to uniformly extract sub-sampled schemes from an existing single or multiple shell scheme, and to order samples in an existing scheme. We propose five algorithms to solve the above problems, including an incremental SC (ISC), a sophisticated greedy algorithm called Iterative Maximum Overlap Construction (IMOC), an 1-Opt greedy method, a Mixed Integer Linear Programming (MILP) method, and a Constrained Non-Linear Optimization (CNLO) method. To our knowledge, this is the first work to use the SC formulation for single or multiple shell sampling schemes in dMRI. Experimental results indicate that SC methods obtain larger angular separation and better rotational invariance than the state-of-the-art EEM and GEEM. The related codes and a tutorial have been released in DMRITool.Comment: Accepted by IEEE transactions on Medical Imaging. Codes have been released in dmritool https://diffusionmritool.github.io/tutorial_qspacesampling.htm