The effect of local perturbation fields on human DTI: characterisation, measurement and correction

Efficient approximate signal reconstruction for correction of gradient nonlinearities in diffusion-weighted imaging

In diffusion weighted MRI (DW-MRI), hardware nonlinearities lead to spatial variations in the orientation and magnitude of diffusion weighting. While the correction of these spatial distortions has been well established for analyses of DW-MRI, the existing voxel-wise empirical correction for gradient nonlinearities requires reimplementation of existing models, as the resultant gradients vary by voxel. Herein, we propose a two-step signal approximation after voxel-wise correction of gradient nonlinearity effects in DW-MRI. The proposed technique (1) scales the diffusion signal and (2) resamples the gradient orientations. This results in uniform gradients across the corrected image and provides the key advantage of seamless integration into current diffusion workflows. We investigated the validity of our technique by fitting a multi-compartment neurite orientation dispersion and density imaging (NODDI) model to the empirical correction and proposed approximation in five subjects from the MASiVar pediatric dataset. We evaluated intra-cellular volume fraction (iVF), CSF volume fraction (cVF), and orientation dispersion index (ODI) from NODDI. The Cohen's d of iVF, cVF and ODI between the techniques was <0.2 indicating the proposed technique does not exhibit significant differences from the voxel-wise correction technique. Our two-step signal approximation is an efficient representation of the voxel-wise gradient table correction. Using this approximation, correction of gradient nonlinearities can be easily incorporated into existing diffusion preprocessing pipelines and is implemented in "PreQual: An automated pipeline for integrated preprocessing and quality assurance of diffusion weighted MRI images".

Hypertensive Adults Exhibit Lower Myelin Content: A Multicomponent Relaxometry and Diffusion Magnetic Resonance Imaging Study

BACKGROUND

It is unknown whether hypertension plays any role in cerebral myelination. To fill this knowledge gap, we studied 90 cognitively unimpaired adults, age range 40 to 94 years, who are participants in the Baltimore Longitudinal Study of Aging and the Genetic and Epigenetic Signatures of Translational Aging Laboratory Testing to look for potential associations between hypertension and cerebral myelin content across 14 white matter brain regions.

METHODS

Myelin content was probed using our advanced multicomponent magnetic resonance relaxometry method of myelin water fraction, a direct and specific magnetic resonance imaging measure of myelin content, and longitudinal and transverse relaxation rates (R1 and R2), 2 highly sensitive magnetic resonance imaging metrics of myelin content. We also applied diffusion tensor imaging magnetic resonance imaging to measure fractional anisotropy, mean diffusivity, radial diffusivity, and axial diffusivity values, which are metrics of cerebral microstructural tissue integrity, to provide context with previous magnetic resonance imaging findings.

RESULTS

After adjustment of age, sex, systolic blood pressure, smoking status, diabetes status, and cholesterol level, our results indicated that participants with hypertension exhibited lower myelin water fraction, fractional anisotropy, R1 and R2 values and higher mean diffusivity, radial diffusivity, and axial diffusivity values, indicating lower myelin content and higher impairment to the brain microstructure. These associations were significant across several white matter regions, particularly in the corpus callosum, fronto-occipital fasciculus, temporal lobes, internal capsules, and corona radiata.

CONCLUSIONS

These original findings suggest a direct association between myelin content and hypertension and form the basis for further investigations including longitudinal assessments of this relationship.

Mapping the impact of nonlinear gradient fields with noise on diffusion MRI

In diffusion MRI, gradient nonlinearities cause spatial variations in the magnitude and direction of diffusion gradients. Studies have shown artifacts from these distortions can results in biased diffusion tensor information and tractography. Here, we investigate the impact of gradient nonlinearity correction in the presence of noise. We introduced empirically derived gradient nonlinear fields at different signal-to-noise ratio (SNR) levels in two experiments: tensor field simulation and simulation of the brain. For each experiment, this work compares two techniques empirically: voxel-wise gradient table correction and approximate correction by scaling the signal directly. The impact was assessed through diffusion metrics including mean diffusivity (MD), fractional anisotropy (FA), axial diffusivity (AD), radial diffusivity (RD), and principal eigen vector (V1). The study shows (1) the correction of gradient nonlinearities will not lead to substantively incorrect estimation of diffusion metrics in a linear system, (2) gradient nonlinearity correction does not interact adversely with noise, (3) nonlinearity correction suppresses the impact of nonlinearities in typical SNR data, (4) for SNR below 30, the performance of both the gradient nonlinearity correction techniques were similar, and (5) larger impacts are seen in regions where the gradient nonlinearities are distinct. Thus, this study suggests that there were greater beneficial effects than adverse effects due to the correction of nonlinearities. Additionally, correction of nonlinearities is recommended when region of interests are in areas with pronounced nonlinearities.

Mapping the human connectome using diffusion MRI at 300 mT/m gradient strength: Methodological advances and scientific impact

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength diffusion MRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for diffusion MRI and where the field is headed in the coming years.

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, B1 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.

Cerebral Oxygen Metabolic Stress, Microstructural Injury, and Infarction in Adults With Sickle Cell Disease

OBJECTIVE

To determine the patient- and tissue-based relationships between cerebral hemodynamic and oxygen metabolic stress, microstructural injury, and infarct location in adults with sickle cell disease (SCD).

METHODS

Control and SCD participants underwent brain MRI to quantify cerebral blood flow (CBF), oxygen extraction fraction (OEF), mean diffusivity (MD), and fractional anisotropy (FA) within normal-appearing white matter (NAWM), and infarcts on FLAIR. Multivariable linear regression examined the patient- and voxel-based associations between hemodynamic and metabolic stress (defined as elevated CBF and OEF, respectively), white matter microstructure, and infarct location.

RESULTS

Of 83 control and SCD participants, adults with SCD demonstrated increased CBF (50.9 vs 38.8 mL/min/100g, p<0.001), increased OEF (0.35 vs 0.25, p<0.001), increased MD (0.76 vs 0.72 x 10^-3mm² s^-1, p=0.005), and decreased FA (0.40 vs 0.42, p=0.021) within NAWM compared to controls. In multivariable analysis, increased OEF (β=0.19, p=0.035), but not CBF (β=0.00, p=0.340), independently predicted increased MD in the SCD cohort, while neither were predictors in controls. On voxel-wise regression, the SCD cohort demonstrated widespread OEF elevation, encompassing deep white matter regions of elevated MD and reduced FA, which spatially extended beyond high density infarct locations from the SCD cohort.

CONCLUSION

Elevated OEF, a putative index of cerebral oxygen metabolic stress, may provide a metric of ischemic vulnerability which could enable individualization of therapeutic strategies in SCD. The patient- and tissue-based relationships between elevated OEF, elevated MD, and cerebral infarcts suggest that oxygen metabolic stress may underlie microstructural injury prior to the development of cerebral infarcts in SCD.

In vivo human whole-brain Connectom diffusion MRI dataset at 760 µm isotropic resolution

We present a whole-brain in vivo diffusion MRI (dMRI) dataset acquired at 760 μm isotropic resolution and sampled at 1260 q-space points across 9 two-hour sessions on a single healthy participant. The creation of this benchmark dataset is possible through the synergistic use of advanced acquisition hardware and software including the high-gradient-strength Connectom scanner, a custom-built 64-channel phased-array coil, a personalized motion-robust head stabilizer, a recently developed SNR-efficient dMRI acquisition method, and parallel imaging reconstruction with advanced ghost reduction algorithm. With its unprecedented resolution, SNR and image quality, we envision that this dataset will have a broad range of investigational, educational, and clinical applications that will advance the understanding of human brain structures and connectivity. This comprehensive dataset can also be used as a test bed for new modeling, sub-sampling strategies, denoising and processing algorithms, potentially providing a common testing platform for further development of in vivo high resolution dMRI techniques. Whole brain anatomical T1-weighted and T2-weighted images at submillimeter scale along with field maps are also made available.

The effect of gradient nonlinearities on fiber orientation estimates from spherical deconvolution of diffusion magnetic resonance imaging data

Gradient nonlinearities in magnetic resonance imaging (MRI) cause spatially varying mismatches between the imposed and the effective gradients and can cause significant biases in rotationally invariant diffusion MRI measures derived from, for example, diffusion tensor imaging. The estimation of the orientational organization of fibrous tissue, which is nowadays frequently performed with spherical deconvolution techniques ideally using higher diffusion weightings, can likewise be biased by gradient nonlinearities. We explore the sensitivity of two established spherical deconvolution approaches to gradient nonlinearities, namely constrained spherical deconvolution (CSD) and damped Richardson-Lucy (dRL). Additionally, we propose an extension of dRL to take into account gradient imperfections, without the need of data interpolation. Simulations show that using the effective b-matrix can improve dRL fiber orientation estimation and reduces angular deviations, while CSD can be more robust to gradient nonlinearity depending on the implementation. Angular errors depend on a complex interplay of many factors, including the direction and magnitude of gradient deviations, underlying microstructure, SNR, anisotropy of the effective response function, and diffusion weighting. Notably, angular deviations can also be observed at lower b-values in contrast to the perhaps common assumption that only high b-value data are affected. In in vivo Human Connectome Project data and acquisitions from an ultrastrong gradient (300 mT/m) scanner, angular differences are observed between applying and not applying the effective gradients in dRL estimation. As even small angular differences can lead to error propagation during tractography and as such impact connectivity analyses, incorporating gradient deviations into the estimation of fiber orientations should make such analyses more reliable.

A comparative study of gradient nonlinearity correction strategies for processing diffusion data obtained with ultra-strong gradient MRI scanners

Purpose

The analysis of diffusion data obtained under large gradient nonlinearities necessitates corrections during data reconstruction and analysis. While two such preprocessing pipelines have been proposed, no comparative studies assessing their performance exist. Furthermore, both pipelines neglect the impact of subject motion during acquisition, which, in the presence of gradient nonlinearities, induces spatio‐temporal B‐matrix variations. Here, spatio‐temporal B‐matrix tracking (STB) is proposed and its performance compared to established pipelines.

Methods

Diffusion tensor MRI (DT‐MRI) was performed using a 300 mT/m gradient system. Data were acquired with volunteers positioned in regions with pronounced gradient nonlinearities, and used to compare the performance of six different processing pipelines, including STB.

Results

Up to 30% errors were observed in DT‐MRI parameter estimates when neglecting gradient nonlinearities. Moreover, the order in which B0 inhomogeneity, eddy current and gradient nonlinearity corrections were performed was found to impact the consistency of parameter estimates significantly. Although, no pipeline emerged as a clear winner, the STB approach seemed to yield the most consistent parameter estimates under large gradient nonlinearities.

Conclusions

Under large gradient nonlinearities, the choice of preprocessing pipeline significantly impacts the estimated diffusion parameters. Motion‐induced spatio‐temporal B‐matrix variations can lead to systematic bias in the parameter estimates, that can be ameliorated using the proposed STB framework.

Click here for author‐reader discussions

A comprehensive approach for correcting voxel-wise b-value errors in diffusion MRI

PURPOSE

In diffusion MRI, the actual b-value played out on the scanner may deviate from the nominal value due to magnetic field imperfections. A simple image-based correction method for this problem is presented.

METHODS

The apparent diffusion constant (ADC) of a water phantom was measured voxel-wise along 64 diffusion directions at b = 1000 s/mm² . The true diffusion constant of water was estimated, considering the phantom temperature. A voxel-wise correction factor, providing an effective b-value including any magnetic field deviations, was determined for each diffusion direction by relating the measured ADC to the true diffusion constant. To test the method, the measured b-value map was used to calculate the corrected voxel-wise ADC for additionally acquired diffusion data sets on the same water phantom and data sets acquired on a small water phantom at three different positions. Diffusion tensor was estimated by applying the measured b-value map to phantom and in vivo data sets.

RESULTS

The b-value-corrected ADC maps of the phantom showed the expected spatial uniformity as well as a marked improvement in consistency across diffusion directions. The b-value correction for the brain data resulted in a 5.8% and 5.5% decrease in mean diffusivity and angular differences of the primary diffusion direction of 2.71° and 0.73° inside gray and white matter, respectively.

CONCLUSION

The actual b-value deviates significantly from its nominal setting, leading to a spatially variable error in the common diffusion outcome measures. The suggested method measures and corrects these artifacts.

Microstructural imaging of the human brain with a 'super-scanner': 10 key advantages of ultra-strong gradients for diffusion MRI

The key component of a microstructural diffusion MRI 'super-scanner' is a dedicated high-strength gradient system that enables stronger diffusion weightings per unit time compared to conventional gradient designs. This can, in turn, drastically shorten the time needed for diffusion encoding, increase the signal-to-noise ratio, and facilitate measurements at shorter diffusion times. This review, written from the perspective of the UK National Facility for In Vivo MR Imaging of Human Tissue Microstructure, an initiative to establish a shared 300 mT/m-gradient facility amongst the microstructural imaging community, describes ten advantages of ultra-strong gradients for microstructural imaging. Specifically, we will discuss how the increase of the accessible measurement space compared to a lower-gradient systems (in terms of Δ, b-value, and TE) can accelerate developments in the areas of 1) axon diameter distribution mapping; 2) microstructural parameter estimation; 3) mapping micro-vs macroscopic anisotropy features with gradient waveforms beyond a single pair of pulsed-gradients; 4) multi-contrast experiments, e.g. diffusion-relaxometry; 5) tractography and high-resolution imaging in vivo and 6) post mortem; 7) diffusion-weighted spectroscopy of metabolites other than water; 8) tumour characterisation; 9) functional diffusion MRI; and 10) quality enhancement of images acquired on lower-gradient systems. We finally discuss practical barriers in the use of ultra-strong gradients, and provide an outlook on the next generation of 'super-scanners'.

The efficiency of retrospective artifact correction methods in improving the statistical power of between-group differences in spinal cord DTI

Diffusion tensor imaging (DTI) is a promising approach for investigating the white matter microstructure of the spinal cord. However, it suffers from severe susceptibility, physiological, and instrumental artifacts present in the cord. Retrospective correction techniques are popular approaches to reduce these artifacts, because they are widely applicable and do not increase scan time.

In this paper, we present a novel outlier rejection approach (reliability masking) which is designed to supplement existing correction approaches by excluding irreversibly corrupted and thus unreliable data points from the DTI index maps. Then, we investigate how chains of retrospective correction techniques including (i) registration, (ii) registration and robust fitting, and (iii) registration, robust fitting, and reliability masking affect the statistical power of a previously reported finding of lower fractional anisotropy values in the posterior column and lateral corticospinal tracts in cervical spondylotic myelopathy (CSM) patients.

While established post-processing steps had small effect on the statistical power of the clinical finding (slice-wise registration: −0.5%, robust fitting: +0.6%), adding reliability masking to the post-processing chain increased it by 4.7%. Interestingly, reliability masking and registration affected the t-score metric differently: while the gain in statistical power due to reliability masking was mainly driven by decreased variability in both groups, registration slightly increased variability. In conclusion, reliability masking is particularly attractive for neuroscience and clinical research studies, as it increases statistical power by reducing group variability and thus provides a cost-efficient alternative to increasing the group size.

•

A novel outlier rejection technique (reliability masking) is introduced.

•

Standard artifact correction has little effect on the statistical power of between-group differences.

•

Reliability masking improves the statistical power of between-group differences.

•

This improvement is driven by decreased group-level variability.

Whole-Brain In-vivo Measurements of the Axonal G-Ratio in a Group of 37 Healthy Volunteers

The g-ratio, quantifying the ratio between the inner and outer diameters of a fiber, is an important microstructural characteristic of fiber pathways and is functionally related to conduction velocity. We introduce a novel method for estimating the MR g-ratio non-invasively across the whole brain using high-fidelity magnetization transfer (MT) imaging and single-shell diffusion MRI. These methods enabled us to map the MR g-ratio in vivo across the brain's prominent fiber pathways in a group of 37 healthy volunteers and to estimate the inter-subject variability. Effective correction of susceptibility-related distortion artifacts was essential before combining the MT and diffusion data, in order to reduce partial volume and edge artifacts. The MR g-ratio is in good qualitative agreement with histological findings despite the different resolution and spatial coverage of MRI and histology. The MR g-ratio holds promise as an important non-invasive biomarker due to its microstructural and functional relevance in neurodegeneration.

Experimental considerations for fast kurtosis imaging

PURPOSE

The clinical use of kurtosis imaging is impeded by long acquisitions and postprocessing. Recently, estimation of mean kurtosis tensor W¯ and mean diffusivity ( D¯) was made possible from 13 distinct diffusion weighted MRI acquisitions (the 1-3-9 protocol) with simple postprocessing. Here, we analyze the effects of noise and nonideal diffusion encoding, and propose a new correction strategy. We also present a 1-9-9 protocol with increased robustness to experimental imperfections and minimal additional scan time. This refinement does not affect computation time and also provides a fast estimate of fractional anisotropy (FA).

THEORY AND METHODS

1-3-9/1-9-9 data are acquired in rat and human brains, and estimates of D¯, FA, W¯ from human brains are compared with traditional estimates from an extensive diffusion kurtosis imaging data set. Simulations are used to evaluate the influence of noise and diffusion encodings deviating from the scheme, and the performance of the correction strategy. Optimal b-values are determined from simulations and data.

RESULTS

Accuracy and precision in D¯ and W¯ are comparable to nonlinear least squares estimation, and is improved with the 1-9-9 protocol. The compensation strategy vastly improves parameter estimation in nonideal data.

CONCLUSION

The framework offers a robust and compact method for estimating several diffusion metrics. The protocol is easily implemented. Magn Reson Med 76:1455-1468, 2016. © 2015 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

POAS4SPM: a toolbox for SPM to denoise diffusion MRI data

We present an implementation of a recently developed noise reduction algorithm for dMRI data, called multi-shell position orientation adaptive smoothing (msPOAS), as a toolbox for SPM. The method intrinsically adapts to the structures of different size and shape in dMRI and hence avoids blurring typically observed in non-adaptive smoothing. We give examples for the usage of the toolbox and explain the determination of experiment-dependent parameters for an optimal performance of msPOAS.

Investigation of vibration-induced artifact in clinical diffusion-weighted imaging of pediatric subjects

It has been reported that mechanical vibrations of the magnetic resonance imaging scanner could produce spurious signal dropouts in diffusion-weighted images resulting in artifactual anisotropy in certain regions of the brain with red appearance in the Directionally Encoded Color maps. We performed a review of the frequency of this artifact across pediatric studies, noting differences by scanner manufacturer, acquisition protocol, as well as weight and position of the subject. We also evaluated the ability of automated and quantitative methods to detect this artifact. We found that the artifact may be present in over 50% of data in certain protocols and is not limited to one scanner manufacturer. While a specific scanner had the highest incidence, low body weight and positioning were also associated with appearance of the artifact for both scanner types evaluated, making children potentially more susceptible than adults. Visual inspection remains the best method for artifact identification. Software for automated detection showed very low sensitivity (10%). The artifact may present inconsistently in longitudinal studies. We discuss a published case report that has been widely cited and used as evidence to set policy about diagnostic criteria for determining vegetative state. That report attributed longitudinal changes in anisotropy to white matter plasticity without considering the possibility that the changes were caused by this artifact. Our study underscores the need to check for the presence of this artifact in clinical studies, analyzes circumstances for when it may be more likely to occur, and suggests simple strategies to identify and potentially avoid its effects.

Diffusion-MRI in neurodegenerative disorders

The ability to image the whole brain through ever more subtle and specific methods/contrasts has come to play a key role in understanding the basis of brain abnormalities in several diseases. In magnetic resonance imaging (MRI), "diffusion" (i.e. the random, thermally-induced displacements of water molecules over time) represents an extraordinarily sensitive contrast mechanism, and the exquisite structural detail it affords has proven useful in a vast number of clinical as well as research applications. Since diffusion-MRI is a truly quantitative imaging technique, the indices it provides can serve as potential imaging biomarkers which could allow early detection of pathological alterations as well as tracking and possibly predicting subtle changes in follow-up examinations and clinical trials. Accordingly, diffusion-MRI has proven useful in obtaining information to better understand the microstructural changes and neurophysiological mechanisms underlying various neurodegenerative disorders. In this review article, we summarize and explore the main applications, findings, perspectives as well as challenges and future research of diffusion-MRI in various neurodegenerative disorders including Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis, Huntington's disease and degenerative ataxias.

Prospective motion correction of 3D echo-planar imaging data for functional MRI using optical tracking

We evaluated the performance of an optical camera based prospective motion correction (PMC) system in improving the quality of 3D echo-planar imaging functional MRI data. An optical camera and external marker were used to dynamically track the head movement of subjects during fMRI scanning. PMC was performed by using the motion information to dynamically update the sequence's RF excitation and gradient waveforms such that the field-of-view was realigned to match the subject's head movement. Task-free fMRI experiments on five healthy volunteers followed a 2 × 2 × 3 factorial design with the following factors: PMC on or off; 3.0 mm or 1.5 mm isotropic resolution; and no, slow, or fast head movements. Visual and motor fMRI experiments were additionally performed on one of the volunteers at 1.5 mm resolution comparing PMC on vs PMC off for no and slow head movements. Metrics were developed to quantify the amount of motion as it occurred relative to k-space data acquisition. The motion quantification metric collapsed the very rich camera tracking data into one scalar value for each image volume that was strongly predictive of motion-induced artifacts. The PMC system did not introduce extraneous artifacts for the no motion conditions and improved the time series temporal signal-to-noise by 30% to 40% for all combinations of low/high resolution and slow/fast head movement relative to the standard acquisition with no prospective correction. The numbers of activated voxels (p < 0.001, uncorrected) in both task-based experiments were comparable for the no motion cases and increased by 78% and 330%, respectively, for PMC on versus PMC off in the slow motion cases. The PMC system is a robust solution to decrease the motion sensitivity of multi-shot 3D EPI sequences and thereby overcome one of the main roadblocks to their widespread use in fMRI studies.

•

Data quality was assessed in terms of tSNR, RMSE, and task fMRI t-statistics.

•

Correction improved tSNR by 30% to 40%.

•

Percent reduction in RMSE as a function of motion level was characterized.

•

Number of active voxels for motor and visual tasks increased significantly.

High-resolution diffusion kurtosis imaging at 3T enabled by advanced post-processing

Diffusion Kurtosis Imaging (DKI) is more sensitive to microstructural differences and can be related to more specific micro-scale metrics (e.g., intra-axonal volume fraction) than diffusion tensor imaging (DTI), offering exceptional potential for clinical diagnosis and research into the white and gray matter. Currently DKI is acquired only at low spatial resolution (2–3 mm isotropic), because of the lower signal-to-noise ratio (SNR) and higher artifact level associated with the technically more demanding DKI. Higher spatial resolution of about 1 mm is required for the characterization of fine white matter pathways or cortical microstructure. We used restricted-field-of-view (rFoV) imaging in combination with advanced post-processing methods to enable unprecedented high-quality, high-resolution DKI (1.2 mm isotropic) on a clinical 3T scanner. Post-processing was advanced by developing a novel method for Retrospective Eddy current and Motion ArtifacT Correction in High-resolution, multi-shell diffusion data (REMATCH). Furthermore, we applied a powerful edge preserving denoising method, denoted as multi-shell orientation-position-adaptive smoothing (msPOAS). We demonstrated the feasibility of high-quality, high-resolution DKI and its potential for delineating highly myelinated fiber pathways in the motor cortex. REMATCH performs robustly even at the low SNR level of high-resolution DKI, where standard EC and motion correction failed (i.e., produced incorrectly aligned images) and thus biased the diffusion model fit. We showed that the combination of REMATCH and msPOAS increased the contrast between gray and white matter in mean kurtosis (MK) maps by about 35% and at the same time preserves the original distribution of MK values, whereas standard Gaussian smoothing strongly biases the distribution.

Adaptive smoothing of multi-shell diffusion weighted magnetic resonance data by msPOAS

We present a novel multi-shell position-orientation adaptive smoothing (msPOAS) method for diffusion weighted magnetic resonance data. Smoothing in voxel and diffusion gradient space is embedded in an iterative adaptive multiscale approach. The adaptive character avoids blurring of the inherent structures and preserves discontinuities. The simultaneous treatment of all q-shells improves the stability compared to single-shell approaches such as the original POAS method. The msPOAS implementation simplifies and speeds up calculations, compared to POAS, facilitating its practical application. Simulations and heuristics support the face validity of the technique and its rigorousness. The characteristics of msPOAS were evaluated on single and multi-shell diffusion data of the human brain. Significant reduction in noise while preserving the fine structure was demonstrated for diffusion weighted images, standard DTI analysis and advanced diffusion models such as NODDI. MsPOAS effectively improves the poor signal-to-noise ratio in highly diffusion weighted multi-shell diffusion data, which is required by recent advanced diffusion micro-structure models. We demonstrate the superiority of the new method compared to other advanced denoising methods.

•

Method for structure preserving smoothing multi-shell dMRI data

•

Does not rely on any dMRI diffusion model

•

Outperforms naive single-shell POAS and other approaches

•

Feasible for real data application

•

Implemented within a freely available package dti for the R Language

MR scanner systems should be adequately characterized in diffusion-MRI of the breast

Breast imaging represents a relatively recent and promising field of application of quantitative diffusion-MRI techniques. In view of the importance of guaranteeing and assessing its reliability in clinical as well as research settings, the aim of this study was to specifically characterize how the main MR scanner system-related factors affect quantitative measurements in diffusion-MRI of the breast. In particular, phantom acquisitions were performed on three 1.5 T MR scanner systems by different manufacturers, all equipped with a dedicated multi-channel breast coil as well as acquisition sequences for diffusion-MRI of the breast. We assessed the accuracy, inter-scan and inter-scanner reproducibility of the mean apparent diffusion coefficient measured along the main orthogonal directions (<ADC>) as well as of diffusion-tensor imaging (DTI)-derived mean diffusivity (MD) measurements. Additionally, we estimated spatial non-uniformity of <ADC> (NU_<ADC>) and MD (NU_MD) maps. We showed that the signal-to-noise ratio as well as overall calibration of high strength diffusion gradients system in typical acquisition sequences for diffusion-MRI of the breast varied across MR scanner systems, introducing systematic bias in the measurements of diffusion indices. While <ADC> and MD values were not appreciably different from each other, they substantially varied across MR scanner systems. The mean of the accuracies of measured <ADC> and MD was in the range [−2.3%,11.9%], and the mean of the coefficients of variation for <ADC> and MD measurements across MR scanner systems was 6.8%. The coefficient of variation for repeated measurements of both <ADC> and MD was < 1%, while NU_<ADC> and NU_MD values were <4%. Our results highlight that MR scanner system-related factors can substantially affect quantitative diffusion-MRI of the breast. Therefore, a specific quality control program for assessing and monitoring the performance of MR scanner systems for diffusion-MRI of the breast is highly recommended at every site, especially in multicenter and longitudinal studies.

The role of diffusion tensor imaging in detecting microstructural changes in prodromal Alzheimer's disease

The MRI technique diffusion tensor imaging (DTI) is reviewed along with microstructural changes associated with prodromal Alzheimer's disease (AD) as a potential biomarker for clinical applications. The prodromal stage of AD is characterized by mild cognitive impairment (MCI), representing a transitional state between normal aging and AD. Microstructural abnormalities on DTI are promising in vivo biomarkers of gray and white matter changes associated with the progression of AD pathology. Elevated mean diffusivity and decreased fractional anisotropy are consistently found in prodromal AD, and even in cognitively normal elderly who progress to MCI. However, quality of parameter maps may be affected by artifacts of motion, susceptibility, and eddy current-induced distortions. The DTI maps are typically analyzed by region-of-interest or voxel-based analytic techniques such as tract-based spatial statistics. DTI-based index of diffusivity is complementary to macrostructural gray matter changes in the hippocampus in detecting prodromal AD. Breakdown of structural connectivity measured with DTI may impact cognitive performance during early AD. Furthermore, assessment of hippocampal connections may help in understanding the cerebral organization and remodeling associated with treatment response.

Quantitative multi-parameter mapping of R1, PD(*), MT, and R2(*) at 3T: a multi-center validation

Multi-center studies using magnetic resonance imaging facilitate studying small effect sizes, global population variance and rare diseases. The reliability and sensitivity of these multi-center studies crucially depend on the comparability of the data generated at different sites and time points. The level of inter-site comparability is still controversial for conventional anatomical T1-weighted MRI data. Quantitative multi-parameter mapping (MPM) was designed to provide MR parameter measures that are comparable across sites and time points, i.e., 1 mm high-resolution maps of the longitudinal relaxation rate (R1 = 1/T1), effective proton density (PD(*)), magnetization transfer saturation (MT) and effective transverse relaxation rate (R2(*) = 1/T2(*)). MPM was validated at 3T for use in multi-center studies by scanning five volunteers at three different sites. We determined the inter-site bias, inter-site and intra-site coefficient of variation (CoV) for typical morphometric measures [i.e., gray matter (GM) probability maps used in voxel-based morphometry] and the four quantitative parameters. The inter-site bias and CoV were smaller than 3.1 and 8%, respectively, except for the inter-site CoV of R2(*) (<20%). The GM probability maps based on the MT parameter maps had a 14% higher inter-site reproducibility than maps based on conventional T1-weighted images. The low inter-site bias and variance in the parameters and derived GM probability maps confirm the high comparability of the quantitative maps across sites and time points. The reliability, short acquisition time, high resolution and the detailed insights into the brain microstructure provided by MPM makes it an efficient tool for multi-center imaging studies.

The impact of post-processing on spinal cord diffusion tensor imaging

Diffusion tensor imaging (DTI) provides information about the microstructure in the brain and spinal cord. While new neuroimaging techniques have significantly advanced the accuracy and sensitivity of DTI of the brain, the quality of spinal cord DTI data has improved less. This is in part due to the small size of the spinal cord (ca. 1cm diameter) and more severe instrumental (e.g. eddy current) and physiological (e.g. cardiac pulsation) artefacts present in spinal cord DTI. So far, the improvements in image quality and resolution have resulted from cardiac gating and new acquisition approaches (e.g. reduced field-of-view techniques). The use of retrospective correction methods is not well established for spinal cord DTI. The aim of this paper is to develop an improved post-processing pipeline tailored for DTI data of the spinal cord with increased quality. For this purpose, we compared two eddy current and motion correction approaches using three-dimensional affine (3D-affine) and slice-wise registrations. We also introduced a new robust-tensor-fitting method that controls for whole-volume outliers. Although in general 3D-affine registration improves data quality, occasionally it can lead to misregistrations and biassed tensor estimates. The proposed robust tensor fitting reduced misregistration-related bias and yielded more reliable tensor estimates. Overall, the combination of slice-wise motion correction, eddy current correction, and robust tensor fitting yielded the best results. It increased the contrast-to-noise ratio (CNR) in FA maps by about 30% and reduced intra-subject variation in fractional anisotropy (FA) maps by 18%. The higher quality of FA maps allows for a better distinction between grey and white matter without increasing scan time and is compatible with any multi-directional DTI acquisition scheme.

Role of magnetic resonance imaging, cerebrospinal fluid, and electroencephalogram in diagnosis of sporadic Creutzfeldt-Jakob disease

Sporadic Creutzfeldt-Jakob disease (sCJD) is a rapidly progressive dementia (RPD) that can be difficult to identify antemortem, with definitive diagnosis requiring tissue confirmation. We describe the clinical, magnetic resonance imaging (MRI), cerebrospinal fluid (CSF), and electroencephalogram (EEG) measures of a small cohort of 30 patients evaluated for RPD. Clinical and diagnostic measures were cross-sectionally obtained from 17 sCJD patients (15 definite, two probable), 13 non-prion rapidly progressive dementia patients (npRPD), and 18 unimpaired controls. In a subset of patients (nine sCJD and nine npRPD) diffusion tensor imaging (DTI) measures [fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD)] were also obtained for the caudate, corpus callosum, posterior limb of the internal capsule, pulvinar, precuneus, and frontal lobe. Differences among groups were assessed by an analysis of variance. Compared to npRPD individuals, sCJD patients had cerebellar dysfunction, significantly higher CSF tau, "positive" CSF 14-3-3, and hyperintensities on diffusion-weighted imaging (DWI) that met previously established imaging criteria for sCJD. EEG changes were similar for the two groups. In addition, sCJD patients had significant decreases in DTI measures (MD, AD, RD but not FA) within the caudate and pulvinar compared to either npRPD patients or unimpaired controls. Our results confirm that CSF abnormalities and MRI (especially DWI) can assist in distinguishing sCJD patients from npRPD patients. Future longitudinal studies using multiple measures (including CSF and MRI) are needed for evaluating pathological changes seen in sCJD patients.

Retrospective correction of physiological noise in DTI using an extended tensor model and peripheral measurements

Diffusion tensor imaging is widely used in research and clinical applications, but this modality is highly sensitive to artefacts. We developed an easy-to-implement extension of the original diffusion tensor model to account for physiological noise in diffusion tensor imaging using measures of peripheral physiology (pulse and respiration), the so-called extended tensor model. Within the framework of the extended tensor model two types of regressors, which respectively modeled small (linear) and strong (nonlinear) variations in the diffusion signal, were derived from peripheral measures. We tested the performance of four extended tensor models with different physiological noise regressors on nongated and gated diffusion tensor imaging data, and compared it to an established data-driven robust fitting method. In the brainstem and cerebellum the extended tensor models reduced the noise in the tensor-fit by up to 23% in accordance with previous studies on physiological noise. The extended tensor model addresses both large-amplitude outliers and small-amplitude signal-changes. The framework of the extended tensor model also facilitates further investigation into physiological noise in diffusion tensor imaging. The proposed extended tensor model can be readily combined with other artefact correction methods such as robust fitting and eddy current correction.

Tracking changes following spinal cord injury: insights from neuroimaging

Traumatic spinal cord injury is often disabling and recovery of function is limited. As a consequence of damage, both spinal cord and brain undergo anatomical and functional changes. Besides clinical measures of recovery, biomarkers that can detect early anatomical and functional changes might be useful in determining clinical outcome—during the course of rehabilitation and recovery—as well as furnishing a tool to evaluate novel treatment interventions and their mechanisms of action. Recent evidence suggests an interesting three-way relationship between neurological deficit and changes in the spinal cord and of the brain and that, importantly, noninvasive magnetic resonance imaging techniques, both structural and functional, provide a sensitive tool to lay out these interactions. This review describes recent findings from multimodal imaging studies of remote anatomical changes (i.e., beyond the lesion site), cortical reorganization, and their relationship to clinical disability. These developments in this field may improve our understanding of effects on the nervous system that are attributable to the injury itself and will allow their distinction from changes that result from rehabilitation (i.e., functional retraining) and from interventions affecting the nervous system directly (i.e., neuroprotection or regeneration).

The influence of spatial registration on detection of cerebral asymmetries using voxel-based statistics of fractional anisotropy images and TBSS

The sensitivity of diffusion tensor imaging (DTI) for detecting microstructural white matter alterations has motivated the application of voxel-based statistics (VBS) to fractional anisotropy (FA) images (FA-VBS). However, detected group differences may depend on the spatial registration method used. The objective of this study was to investigate the influence of spatial registration on detecting cerebral asymmetries in FA-VBS analyses with reference to data obtained using Tract-Based Spatial Statistics (TBSS). In the first part of this study we performed FA-VBS analyses using three single-contrast and one multi-contrast registration: (i) whole-brain registration based on T2 contrast, (ii) whole-brain registration based on FA contrast, (iii) individual-hemisphere registration based on FA contrast, and (iv) a combination of (i) and (iii). We then compared the FA-VBS results with those obtained from TBSS. We found that the FA-VBS results depended strongly on the employed registration approach, with the best correspondence between FA-VBS and TBSS results when approach (iv), the “multi-contrast individual-hemisphere” method was employed. In the second part of the study, we investigated the spatial distribution of residual misregistration for each registration approach and the effect on FA-VBS results. For the FA-VBS analyses using the three single-contrast registration methods, we identified FA asymmetries that were (a) located in regions prone to misregistrations, (b) not detected by TBSS, and (c) specific to the applied registration approach. These asymmetries were considered candidates for apparent FA asymmetries due to systematic misregistrations associated with the FA-VBS approach. Finally, we demonstrated that the “multi-contrast individual-hemisphere” approach showed the least residual spatial misregistrations and thus might be most appropriate for cerebral FA-VBS analyses.