Q-space truncation and sampling in diffusion spectrum imaging

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.

Comprehensive diffusion MRI dataset for in vivo human brain microstructure mapping using 300 mT/m gradients

Strong gradient systems can improve the signal-to-noise ratio of diffusion MRI measurements and enable a wider range of acquisition parameters that are beneficial for microstructural imaging. We present a comprehensive diffusion MRI dataset of 26 healthy participants acquired on the MGH-USC 3 T Connectome scanner equipped with 300 mT/m maximum gradient strength and a custom-built 64-channel head coil. For each participant, the one-hour long acquisition systematically sampled the accessible diffusion measurement space, including two diffusion times (19 and 49 ms), eight gradient strengths linearly spaced between 30 mT/m and 290 mT/m for each diffusion time, and 32 or 64 uniformly distributed directions. The diffusion MRI data were preprocessed to correct for gradient nonlinearity, eddy currents, and susceptibility induced distortions. In addition, scan/rescan data from a subset of seven individuals were also acquired and provided. The MGH Connectome Diffusion Microstructure Dataset (CDMD) may serve as a test bed for the development of new data analysis methods, such as fiber orientation estimation, tractography and microstructural modelling.

High-fidelity approximation of grid- and shell-based sampling schemes from undersampled DSI using compressed sensing: Post mortem validation

While many useful microstructural indices, as well as orientation distribution functions, can be obtained from multi-shell dMRI data, there is growing interest in exploring the richer set of microstructural features that can be extracted from the full ensemble average propagator (EAP). The EAP can be readily computed from diffusion spectrum imaging (DSI) data, at the cost of a very lengthy acquisition. Compressed sensing (CS) has been used to make DSI more practical by reducing its acquisition time. CS applied to DSI (CS-DSI) attempts to reconstruct the EAP from significantly undersampled q-space data. We present a post mortem validation study where we evaluate the ability of CS-DSI to approximate not only fully sampled DSI but also multi-shell acquisitions with high fidelity. Human brain samples are imaged with high-resolution DSI at 9.4T and with polarization-sensitive optical coherence tomography (PSOCT). The latter provides direct measurements of axonal orientations at microscopic resolutions, allowing us to evaluate the mesoscopic orientation estimates obtained from diffusion MRI, in terms of their angular error and the presence of spurious peaks. We test two fast, dictionary-based, L2-regularized algorithms for CS-DSI reconstruction. We find that, for a CS acceleration factor of R=3, i.e., an acquisition with 171 gradient directions, one of these methods is able to achieve both low angular error and low number of spurious peaks. With a scan length similar to that of high angular resolution multi-shell acquisition schemes, this CS-DSI approach is able to approximate both fully sampled DSI and multi-shell data with high accuracy. Thus it is suitable for orientation reconstruction and microstructural modeling techniques that require either grid- or shell-based acquisitions. We find that the signal-to-noise ratio (SNR) of the training data used to construct the dictionary can have an impact on the accuracy of CS-DSI, but that there is substantial robustness to loss of SNR in the test data. Finally, we show that, as the CS acceleration factor increases beyond R=3, the accuracy of these reconstruction methods degrade, either in terms of the angular error, or in terms of the number of spurious peaks. Our results provide useful benchmarks for the future development of even more efficient q-space acceleration techniques.

Comparison of q-Space Reconstruction Methods for Undersampled Diffusion Spectrum Imaging Data

PURPOSE

To compare different q-space reconstruction methods for undersampled diffusion spectrum imaging data.

MATERIALS AND METHODS

We compared the quality of three methods: Mean Apparent Propagator (MAP); Compressed Sensing using Identity (CSI) and Compressed Sensing using Dictionary (CSD) with simulated data and in vivo acquisitions. We used retrospective undersampling so that the fully sampled reconstruction could be used as ground truth. We used the normalized mean squared error (NMSE) and the Pearson's correlation coefficient as reconstruction quality indices. Additionally, we evaluated two propagator-based diffusion indices: mean squared displacement and return to zero probability. We also did a visual analysis around the centrum semiovale.

RESULTS

All methods had reconstruction errors below 5% with low undersampling factors and with a wide range of noise levels. However, the CSD method had at least 1-2% lower NMSE than the other reconstruction methods at higher noise levels. MAP was the second-best method when using a sufficiently high number of q-space samples. MAP reconstruction showed better propagator-based diffusion indices for in vivo acquisitions. With undersampling factors greater than 4, MAP and CSI have noticeably more reconstruction error than CSD.

CONCLUSION

Undersampled data were best reconstructed by means of CSD in simulations and in vivo. MAP was more accurate in the extraction of propagator-based indices, particularly for in vivo data.

Modelling white matter with spherical deconvolution: How and why?

Since the realization that diffusion MRI can probe the microstructural organization and orientation of biological tissue in vivo and non-invasively, a multitude of diffusion imaging methods have been developed and applied to study the living human brain. Diffusion tensor imaging was the first model to be widely adopted in clinical and neuroscience research, but it was also clear from the beginning that it suffered from limitations when mapping complex configurations, such as crossing fibres. In this review, we highlight the main steps that have led the field of diffusion imaging to move from the tensor model to the adoption of diffusion and fibre orientation density functions as a more effective way to describe the complexity of white matter organization within each brain voxel. Among several techniques, spherical deconvolution has emerged today as one of the main approaches to model multiple fibre orientations and for tractography applications. Here we illustrate the main concepts and the reasoning behind this technique, as well as the latest developments in the field. The final part of this review provides practical guidelines and recommendations on how to set up processing and acquisition protocols suitable for spherical deconvolution.

Generalized diffusion spectrum magnetic resonance imaging (GDSI) for model-free reconstruction of the ensemble average propagator

Diffusion spectrum MRI (DSI) provides model-free estimation of the diffusion ensemble average propagator (EAP) and orientation distribution function (ODF) but requires the diffusion data to be acquired on a Cartesian q-space grid. Multi-shell diffusion acquisitions are more flexible and more commonly acquired but have, thus far, only been compatible with model-based analysis methods. Here, we propose a generalized DSI (GDSI) framework to recover the EAP from multi-shell diffusion MRI data. The proposed GDSI approach corrects for q-space sampling density non-uniformity using a fast geometrical approach. The EAP is directly calculated in a preferable coordinate system by multiplying the sampling density corrected q-space signals by a discrete Fourier transform matrix, without any need for gridding. The EAP is demonstrated as a way to map diffusion patterns in brain regions such as the thalamus, cortex and brainstem where the tissue microstructure is not as well characterized as in white matter. Scalar metrics such as the zero displacement probability and displacement distances at different fractions of the zero displacement probability were computed from the recovered EAP to characterize the diffusion pattern within each voxel. The probability averaged across directions at a specific displacement distance provides a diffusion property based image contrast that clearly differentiates tissue types. The displacement distance at the first zero crossing of the EAP averaged across directions orthogonal to the primary fiber orientation in the corpus callosum is found to be larger in the body (5.65 ± 0.09 μm) than in the genu (5.55 ± 0.15 μm) and splenium (5.4 ± 0.15 μm) of the corpus callosum, which corresponds well to prior histological studies. The EAP also provides model-free representations of angular structure such as the diffusion ODF, which allows estimation and comparison of fiber orientations from both the model-free and model-based methods on the same multi-shell data. For the model-free methods, detection of crossing fibers is found to be strongly dependent on the maximum b-value and less sensitive compared to the model-based methods. In conclusion, our study provides a generalized DSI approach that allows flexible reconstruction of the diffusion EAP and ODF from multi-shell diffusion data and data acquired with other sampling patterns.

Simultaneous NODDI and GFA parameter map generation from subsampled q-space imaging using deep learning

PURPOSE

To develop a robust multidimensional deep-learning based method to simultaneously generate accurate neurite orientation dispersion and density imaging (NODDI) and generalized fractional anisotropy (GFA) parameter maps from undersampled q-space datasets for use in stroke imaging.

METHODS

Traditional diffusion spectrum imaging (DSI) capable of producing accurate NODDI and GFA parameter maps requires hundreds of q-space samples which renders the scan time clinically untenable. A convolutional neural network (CNN) was trained to generated NODDI and GFA parameter maps simultaneously from 10× undersampled q-space data. A total of 48 DSI scans from 15 stroke patients and 14 normal subjects were acquired for training, validating, and testing this method. The proposed network was compared to previously proposed voxel-wise machine learning based approaches for q-space imaging. Network-generated images were used to predict stroke functional outcome measures.

RESULTS

The proposed network achieves significant performance advantages compared to previously proposed machine learning approaches, showing significant improvements across image quality metrics. Generating these parameter maps using CNNs also comes with the computational benefits of only needing to generate and train a single network instead of multiple networks for each parameter type. Post-stroke outcome prediction metrics do not appreciably change when using images generated from this proposed technique. Over three test participants, the predicted stroke functional outcome scores were within 1-6% of the clinical evaluations.

CONCLUSIONS

Estimates of NODDI and GFA parameters estimated simultaneously with a deep learning network from highly undersampled q-space data were improved compared to other state-of-the-art methods providing a 10-fold reduction scan time compared to conventional methods.

Diffusion MRI tractography for improved transcranial MRI-guided focused ultrasound thalamotomy targeting for essential tremor

Purpose

To evaluate the use of diffusion magnetic resonance imaging (MRI) tractography for neurosurgical guidance of transcranial MRI-guided focused ultrasound (tcMRgFUS) thalamotomy for essential tremor (ET).

Materials and methods

Eight patients with medication-refractory ET were treated with tcMRgFUS targeting the ventral intermediate nucleus (Vim) of the thalamus contralateral to their dominant hand. Diffusion and structural MRI data and clinical evaluations were acquired pre-treatment and post-treatment. To identify the optimal target location, tractography was performed on pre-treatment diffusion MRI data between the treated thalamus and the hand-knob region of the ipsilateral motor cortex, the entire ipsilateral motor cortex and the contralateral dentate nucleus. The tractography-identified locations were compared to the lesion location delineated on 1 year post-treatment T₂-weighted MR image. Their overlap was correlated with the clinical outcomes measured by the percentage change of the Clinical Rating Scale for Tremor scores acquired pre-treatment, as well as 1 month, 3 months, 6 months and 1 year post-treatment.

Results

The probabilistic tractography was consistent from subject-to-subject and followed the expected anatomy of the thalamocortical radiation and the dentatothalamic tract. Higher overlap between the tractography-identified location and the tcMRgFUS treatment-induced lesion highly correlated with better treatment outcome (r = −0.929, −0.75, −0.643, p = 0.00675, 0.0663, 0.139 for the tractography between the treated thalamus and the hand-knob region of the ipsilateral motor cortex, the entire ipsilateral motor cortex and the contralateral dentate nucleus, respectively, at 1 year post-treatment). The correlation for the tractography between the treated thalamus and the hand-knob region of the ipsilateral motor cortex is the highest for all time points (r = −0.719, −0.976, −0.707, −0.929, p = 0.0519, 0.000397, 0.0595, 0.00675 at 1 month, 3 months, 6 months and 1 year post-treatment, respectively).

Conclusion

Our data support the use of diffusion tractography as a complementary approach to current targeting methods for tcMRgFUS thalamotomy.

•

Retrospectively used tractography to define a target for MRgFUS thalamotomy for ET.

•

Larger overlap between tractography and lesion correlates with better outcomes.

•

Strongest correlations for tract between the thalamus and motor hand-knob region

•

Diffusion tractography is a complementary approach to current targeting methods.

Double diffusion encoding MRI for the clinic

PURPOSE

The purpose of this study is to develop double diffusion encoding (DDE) MRI methods for clinical use. Microscopic diffusion anisotropy measurements from DDE promise greater specificity to changes in tissue microstructure compared with conventional diffusion tensor imaging, but implementation of DDE sequences on whole-body MRI scanners is challenging because of the limited gradient strengths and lengthy acquisition times.

METHODS

A custom single-refocused DDE sequence was implemented on a 3T whole-body scanner. The DDE gradient orientation scheme and sequence parameters were optimized based on a Gaussian diffusion assumption. Using an optimized 5-min DDE acquisition, microscopic fractional anisotropy (μFA) maps were acquired for the first time in multiple sclerosis patients.

RESULTS

Based on simulations and in vivo human measurements, six parallel and six orthogonal diffusion gradient pairs were found to be the minimum number of diffusion gradient pairs necessary to produce a rotationally invariant measurement of μFA. Simulations showed that optimal precision and accuracy of μFA measurements were obtained using b-values between 1500 and 3000 s/mm² . The μFA maps showed improved delineation of multiple sclerosis lesions compared with conventional fractional anisotropy and distinct contrast from T₂ -weighted fluid attenuated inversion recovery and T₁ -weighted imaging.

CONCLUSION

The μFA maps can be measured using DDE in a clinical setting and may provide new opportunities for characterizing multiple sclerosis lesions and other types of tissue degeneration. Magn Reson Med 80:507-520, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

A theoretical signal processing framework for linear diffusion MRI: Implications for parameter estimation and experiment design

The data measured in diffusion MRI can be modeled as the Fourier transform of the Ensemble Average Propagator (EAP), a probability distribution that summarizes the molecular diffusion behavior of the spins within each voxel. This Fourier relationship is potentially advantageous because of the extensive theory that has been developed to characterize the sampling requirements, accuracy, and stability of linear Fourier reconstruction methods. However, existing diffusion MRI data sampling and signal estimation methods have largely been developed and tuned without the benefit of such theory, instead relying on approximations, intuition, and extensive empirical evaluation. This paper aims to address this discrepancy by introducing a novel theoretical signal processing framework for diffusion MRI. The new framework can be used to characterize arbitrary linear diffusion estimation methods with arbitrary q-space sampling, and can be used to theoretically evaluate and compare the accuracy, resolution, and noise-resilience of different data acquisition and parameter estimation techniques. The framework is based on the EAP, and makes very limited modeling assumptions. As a result, the approach can even provide new insight into the behavior of model-based linear diffusion estimation methods in contexts where the modeling assumptions are inaccurate. The practical usefulness of the proposed framework is illustrated using both simulated and real diffusion MRI data in applications such as choosing between different parameter estimation methods and choosing between different q-space sampling schemes.

Accelerated radial diffusion spectrum imaging using a multi-echo stimulated echo diffusion sequence

PURPOSE

Diffusion spectrum imaging (DSI) provides us non-invasively and robustly with anatomical details of brain microstructure. To achieve sufficient angular resolution, DSI requires a large number of q-space samples, leading to long acquisition times. This need is mitigated here by combining the beneficial properties of Radial q-space sampling for DSI with a Multi-Echo Stimulated Echo Sequence (MESTIM).

METHODS

Full 2D k-spaces for each of several q-space samples, along the same radially outward line in q-space, are acquired in one readout train with one spin and three stimulated echoes. RF flip angles are carefully chosen to distribute spin magnetization over the echoes and the DSI reconstruction is adapted to account for differences in diffusion time among echoes.

RESULTS

Individual datasets and bootstrapped reproducibility analysis demonstrate image quality and SNR of the more-than-twofold-accelerated RDSI MESTIM sequence. Orientation distribution functions (ODF) and tractography results benefit from the longer diffusion times of the latter echoes in the echo train.

CONCLUSION

A MESTIM sequence can be used to shorten RDSI acquisition times significantly without loss of image or ODF quality. Further acceleration is possible by combination with simultaneous multi-slice techniques. Magn Reson Med 79:306-316, 2018. © 2017 International Society for Magnetic Resonance in Medicine.