Linear multi-scale modeling of diffusion MRI data: A framework for characterization of oriented structures across length scales.

Diffusion-weighted magnetic resonance imaging (DW-MRI) has evolved to provide increasingly sophisticated investigations of the human brain's structural connectome in vivo. Restriction spectrum imaging (RSI) is a method that reconstructs the orientation distribution of diffusion within tissues over a range of length scales. In its original formulation, RSI represented the signal as consisting of a spectrum of Gaussian diffusion response functions. Recent technological advances have enabled the use of ultra-high b-values on human MRI scanners, providing higher sensitivity to intracellular water diffusion in the living human brain. To capture the complex diffusion time dependence of the signal within restricted water compartments, we expand upon the RSI approach to represent restricted water compartments with non-Gaussian response functions, in an extended analysis framework called linear multi-scale modeling (LMM). The LMM approach is designed to resolve length scale and orientation-specific information with greater specificity to tissue microstructure in the restricted and hindered compartments, while retaining the advantages of the RSI approach in its implementation as a linear inverse problem. Using multi-shell, multi-diffusion time DW-MRI data acquired with a state-of-the-art 3 T MRI scanner equipped with 300 mT/m gradients, we demonstrate the ability of the LMM approach to distinguish different anatomical structures in the human brain and the potential to advance mapping of the human connectome through joint estimation of the fiber orientation distributions and compartment size characteristics.

SDnDTI: Self-supervised deep learning-based denoising for diffusion tensor MRI.

Diffusion tensor magnetic resonance imaging (DTI) is a widely adopted neuroimaging method for the in vivo mapping of brain tissue microstructure and white matter tracts. Nonetheless, the noise in the diffusion-weighted images (DWIs) decreases the accuracy and precision of DTI derived microstructural parameters and leads to prolonged acquisition time for achieving improved signal-to-noise ratio (SNR). Deep learning-based image denoising using convolutional neural networks (CNNs) has superior performance but often requires additional high-SNR data for supervising the training of CNNs, which reduces the feasibility of supervised learning-based denoising in practice. In this work, we develop a self-supervised deep learning-based method entitled "SDnDTI" for denoising DTI data, which does not require additional high-SNR data for training. Specifically, SDnDTI divides multi-directional DTI data into many subsets of six DWI volumes and transforms DWIs from each subset to along the same diffusion-encoding directions through the diffusion tensor model, generating multiple repetitions of DWIs with identical image contrasts but different noise observations. SDnDTI removes noise by first denoising each repetition of DWIs using a deep 3-dimensional CNN with the average of all repetitions with higher SNR as the training target, following the same approach as normal supervised learning based denoising methods, and then averaging CNN-denoised images for achieving higher SNR. The denoising efficacy of SDnDTI is demonstrated in terms of the similarity of output images and resultant DTI metrics compared to the ground truth generated using substantially more DWI volumes on two datasets with different spatial resolutions, b-values and numbers of input DWI volumes provided by the Human Connectome Project (HCP) and the Lifespan HCP in Aging. The SDnDTI results preserve image sharpness and textural details and substantially improve upon those from the raw data. The results of SDnDTI are comparable to those from supervised learning-based denoising and outperform those from state-of-the-art conventional denoising algorithms including BM4D, AONLM and MPPCA. By leveraging domain knowledge of diffusion MRI physics, SDnDTI makes it easier to use CNN-based denoising methods in practice and has the potential to benefit a wider range of research and clinical applications that require accelerated DTI acquisition and high-quality DTI data for mapping of tissue microstructure, fiber tracts and structural connectivity in the living human brain.

Estimating axial diffusivity in the NODDI model.

To estimate microstructure-related parameters from diffusion MRI data, biophysical models make strong, simplifying assumptions about the underlying tissue. The extent to which many of these assumptions are valid remains an open research question. This study was inspired by the disparity between the estimated intra-axonal axial diffusivity from literature and that typically assumed by the Neurite Orientation Dispersion and Density Imaging (NODDI) model (d∥=1.7µm2/ms). We first demonstrate how changing the assumed axial diffusivity results in considerably different NODDI parameter estimates. Second, we illustrate the ability to estimate axial diffusivity as a free parameter of the model using high b-value data and an adapted NODDI framework. Using both simulated and in vivo data we investigate the impact of fitting to either real-valued or magnitude data, with Gaussian and Rician noise characteristics respectively, and what happens if we get the noise assumptions wrong in this high b-value and thus low SNR regime. Our results from real-valued human data estimate intra-axonal axial diffusivities of ∼2-2.5µm2/ms, in line with current literature. Crucially, our results demonstrate the importance of accounting for both a rectified noise floor and/or a signal offset to avoid biased parameter estimates when dealing with low SNR data.

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.

Improving in vivo human cerebral cortical surface reconstruction using data-driven super-resolution.

Accurate and automated reconstruction of the in vivo human cerebral cortical surface from anatomical magnetic resonance (MR) images facilitates the quantitative analysis of cortical structure. Anatomical MR images with sub-millimeter isotropic spatial resolution improve the accuracy of cortical surface and thickness estimation compared to the standard 1-millimeter isotropic resolution. Nonetheless, sub-millimeter resolution acquisitions require averaging multiple repetitions to achieve sufficient signal-to-noise ratio and are therefore long and potentially vulnerable to subject motion. We address this challenge by synthesizing sub-millimeter resolution images from standard 1-millimeter isotropic resolution images using a data-driven supervised machine learning-based super-resolution approach achieved via a deep convolutional neural network. We systematically characterize our approach using a large-scale simulated dataset and demonstrate its efficacy in empirical data. The super-resolution data provide improved cortical surfaces similar to those obtained from native sub-millimeter resolution data. The whole-brain mean absolute discrepancy in cortical surface positioning and thickness estimation is below 100 µm at the single-subject level and below 50 µm at the group level for the simulated data, and below 200 µm at the single-subject level and below 100 µm at the group level for the empirical data, making the accuracy of cortical surfaces derived from super-resolution sufficient for most applications.

Improved cortical surface reconstruction using sub-millimeter resolution MPRAGE by image denoising.

Automatic cerebral cortical surface reconstruction is a useful tool for cortical anatomy quantification, analysis and visualization. Recently, the Human Connectome Project and several studies have shown the advantages of using T1-weighted magnetic resonance (MR) images with sub-millimeter isotropic spatial resolution instead of the standard 1-mm isotropic resolution for improved accuracy of cortical surface positioning and thickness estimation. Nonetheless, sub-millimeter resolution images are noisy by nature and require averaging multiple repetitions to increase the signal-to-noise ratio for precisely delineating the cortical boundary. The prolonged acquisition time and potential motion artifacts pose significant barriers to the wide adoption of cortical surface reconstruction at sub-millimeter resolution for a broad range of neuroscientific and clinical applications. We address this challenge by evaluating the cortical surface reconstruction resulting from denoised single-repetition sub-millimeter T1-weighted images. We systematically characterized the effects of image denoising on empirical data acquired at 0.6 mm isotropic resolution using three classical denoising methods, including denoising convolutional neural network (DnCNN), block-matching and 4-dimensional filtering (BM4D) and adaptive optimized non-local means (AONLM). The denoised single-repetition images were found to be highly similar to 6-repetition averaged images, with a low whole-brain averaged mean absolute difference of ~0.016, high whole-brain averaged peak signal-to-noise ratio of ~33.5 dB and structural similarity index of ~0.92, and minimal gray matter-white matter contrast loss (2% to 9%). The whole-brain mean absolute discrepancies in gray matter-white matter surface placement, gray matter-cerebrospinal fluid surface placement and cortical thickness estimation were lower than 165 µm, 155 µm and 145 µm-sufficiently accurate for most applications. These discrepancies were approximately one third to half of those from 1-mm isotropic resolution data. The denoising performance was equivalent to averaging ~2.5 repetitions of the data in terms of image similarity, and 1.6-2.2 repetitions in terms of the cortical surface placement accuracy. The scan-rescan variability of the cortical surface positioning and thickness estimation was lower than 170 µm. Our unique dataset and systematic characterization support the use of denoising methods for improved cortical surface reconstruction at sub-millimeter resolution.

Scan-rescan repeatability of axonal imaging metrics using high-gradient diffusion MRI and statistical implications for study design.

Axon diameter mapping using diffusion MRI in the living human brain has attracted growing interests with the increasing availability of high gradient strength MRI systems. A systematic assessment of the consistency of axon diameter estimates within and between individuals is needed to gain a comprehensive understanding of how such methods extend to quantifying differences in axon diameter index between groups and facilitate the design of neurobiological studies using such measures. We examined the scan-rescan repeatability of axon diameter index estimation based on the spherical mean technique (SMT) approach using diffusion MRI data acquired with gradient strengths up to 300 mT/m on a 3T Connectom system in 7 healthy volunteers. We performed statistical power analyses using data acquired with the same protocol in a larger cohort consisting of 15 healthy adults to investigate the implications for study design. Results revealed a high degree of repeatability in voxel-wise restricted volume fraction estimates and tract-wise estimates of axon diameter index derived from high-gradient diffusion MRI data. On the region of interest (ROI) level, across white matter tracts in the whole brain, the Pearson's correlation coefficient of the axon diameter index estimated between scan and rescan experiments was r = 0.72 with an absolute deviation of 0.18 µm. For an anticipated 10% effect size in studies of axon diameter index, most white matter regions required a sample size of less than 15 people to observe a measurable difference between groups using an ROI-based approach. To facilitate the use of high-gradient strength diffusion MRI data for neuroscientific studies of axonal microstructure, the comprehensive multi-gradient strength, multi-diffusion time data used in this work will be made publicly available, in support of open science and increasing the accessibility of such data to the greater scientific community.

Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome.

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

Distortion-free, high-isotropic-resolution diffusion MRI with gSlider BUDA-EPI and multicoil dynamic B0 shimming.

PURPOSE: We combine SNR-efficient acquisition and model-based reconstruction strategies with newly available hardware instrumentation to achieve distortion-free in vivo diffusion MRI of the brain at submillimeter-isotropic resolution with high fidelity and sensitivity on a clinical 3T scanner. METHODS: We propose blip-up/down acquisition (BUDA) for multishot EPI using interleaved blip-up/blip-down phase encoding and incorporate B0 forward-modeling into structured low-rank reconstruction to enable distortion-free and navigator-free diffusion MRI. We further combine BUDA-EPI with an SNR-efficient simultaneous multislab acquisition (generalized slice-dithered enhanced resolution ["gSlider"]), to achieve high-isotropic-resolution diffusion MRI. To validate gSlider BUDA-EPI, whole-brain diffusion data at 860-µm and 780-µm data sets were acquired. Finally, to improve the conditioning and minimize noise penalty in BUDA reconstruction at very high resolutions where B0 inhomogeneity can have a detrimental effect, the level of B0 inhomogeneity was reduced by incorporating slab-by-slab dynamic shimming with a 32-channel AC/DC coil into the acquisition. Whole-brain 600-µm diffusion data were then acquired with this combined approach of gSlider BUDA-EPI with dynamic shimming. RESULTS: The results of 860-µm and 780-µm datasets show high geometry fidelity with gSlider BUDA-EPI. With dynamic shimming, the BUDA reconstruction's noise penalty was further alleviated. This enables whole-brain 600-µm isotropic resolution diffusion imaging with high image quality. CONCLUSIONS: The gSlider BUDA-EPI method enables high-quality, distortion-free diffusion imaging across the whole brain at submillimeter resolution, where the use of multicoil dynamic B0 shimming further improves reconstruction performance, which can be particularly useful at very high resolutions.

DeepDTI: High-fidelity six-direction diffusion tensor imaging using deep learning.

Diffusion tensor magnetic resonance imaging (DTI) is unsurpassed in its ability to map tissue microstructure and structural connectivity in the living human brain. Nonetheless, the angular sampling requirement for DTI leads to long scan times and poses a critical barrier to performing high-quality DTI in routine clinical practice and large-scale research studies. In this work we present a new processing framework for DTI entitled DeepDTI that minimizes the data requirement of DTI to six diffusion-weighted images (DWIs) required by conventional voxel-wise fitting methods for deriving the six unique unknowns in a diffusion tensor using data-driven supervised deep learning. DeepDTI maps the input non-diffusion-weighted (b â€‹= â€‹0) image and six DWI volumes sampled along optimized diffusion-encoding directions, along with T1-weighted and T2-weighted image volumes, to the residuals between the input and high-quality output b = 0 image and DWI volumes using a 10-layer three-dimensional convolutional neural network (CNN). The inputs and outputs of DeepDTI are uniquely formulated, which not only enables residual learning to boost CNN performance but also enables tensor fitting of resultant high-quality DWIs to generate orientational DTI metrics for tractography. The very deep CNN used by DeepDTI leverages the redundancy in local and non-local spatial information and across diffusion-encoding directions and image contrasts in the data. The performance of DeepDTI was systematically quantified in terms of the quality of the output images, DTI metrics, DTI-based tractography and tract-specific analysis results. We demonstrate rotationally-invariant and robust estimation of DTI metrics from DeepDTI that are comparable to those obtained with two b â€‹= â€‹0 images and 21 DWIs for the primary eigenvector derived from DTI and two b â€‹= â€‹0 images and 26-30 DWIs for various scalar metrics derived from DTI, achieving 3.3-4.6 נ​acceleration, and twice as good as those of a state-of-the-art denoising algorithm at the group level. The twenty major white-matter tracts can be accurately identified from the tractography of DeepDTI results. The mean distance between the core of the major white-matter tracts identified from DeepDTI results and those from the ground-truth results using 18 â€‹b â€‹= â€‹0 images and 90 DWIs measures around 1-1.5 â€‹mm. DeepDTI leverages domain knowledge of diffusion MRI physics and power of deep learning to render DTI, DTI-based tractography, major white-matter tracts identification and tract-specific analysis more feasible for a wider range of neuroscientific and clinical studies.

Axon diameter index estimation independent of fiber orientation distribution using high-gradient diffusion MRI.

Axon diameter mapping using high-gradient diffusion MRI has generated great interest as a noninvasive tool for studying trends in axonal size in the human brain. One of the main barriers to mapping axon diameter across the whole brain is accounting for complex white matter fiber configurations (e.g., crossings and fanning), which are prevalent throughout the brain. Here, we present a framework for generalizing axon diameter index estimation to the whole brain independent of the underlying fiber orientation distribution using the spherical mean technique (SMT). This approach is shown to significantly benefit from the use of real-valued diffusion data with Gaussian noise, which reduces the systematic bias in the estimated parameters resulting from the elevation of the noise floor when using magnitude data with Rician noise. We demonstrate the feasibility of obtaining whole-brain orientationally invariant estimates of axon diameter index and relative volume fractions in six healthy human volunteers using real-valued diffusion data acquired on a dedicated high-gradient 3-Tesla human MRI scanner with 300 mT/m maximum gradient strength. The trends in axon diameter index are consistent with known variations in axon diameter from histology and demonstrate the potential of this generalized framework for revealing coherent patterns in axonal structure throughout the living human brain. The use of real-valued diffusion data provides a viable solution for eliminating the Rician noise floor and should be considered for all spherical mean approaches to microstructural parameter estimation.

Registration-free analysis of diffusion MRI tractography data across subjects through the human lifespan.

Diffusion MRI tractography produces massive sets of streamlines that need to be clustered into anatomically meaningful white-matter bundles. Conventional clustering techniques group streamlines based on their proximity in Euclidean space. We have developed AnatomiCuts, an unsupervised method for clustering tractography streamlines based on their neighboring anatomical structures, rather than their coordinates in Euclidean space. In this work, we show that the anatomical similarity metric used in AnatomiCuts can be extended to find corresponding clusters across subjects and across hemispheres, without inter-subject or inter-hemispheric registration. Our proposed approach enables group-wise tract cluster analysis, as well as studies of hemispheric asymmetry. We evaluate our approach on data from the pilot MGH-Harvard-USC Lifespan Human Connectome project, showing improved correspondence in tract clusters across 184 subjects aged 8-90. Our method shows up to 38% improvement in the overlap of corresponding clusters when comparing subjects with large age differences. The techniques presented here do not require registration to a template and can thus be applied to populations with large inter-subject variability, e.g., due to brain development, aging, or neurological disorders.

Fast submillimeter diffusion MRI using gSlider-SMS and SNR-enhancing joint reconstruction.

PURPOSE: We evaluate a new approach for achieving diffusion MRI data with high spatial resolution, large volume coverage, and fast acquisition speed. THEORY AND METHODS: A recent method called gSlider-SMS enables whole-brain submillimeter diffusion MRI with high signal-to-noise ratio (SNR) efficiency. However, despite the efficient acquisition, the resulting images can still suffer from low SNR due to the small size of the imaging voxels. This work proposes to mitigate the SNR problem by combining gSlider-SMS with a regularized SNR-enhancing reconstruction approach. RESULTS: Illustrative results show that, from gSlider-SMS data acquired over a span of only 25 minutes on a 3T scanner, the proposed method is able to produce 71 MRI images (64 diffusion encoding orientations with b = 1500 s/ mm2 , and 7 images without diffusion weighting) of the entire in vivo human brain with nominal 0.66 mm spatial resolution. Using data acquired from 75 minutes of acquisition as a gold standard reference, we demonstrate that the proposed SNR-enhancement procedure leads to substantial improvements in estimated diffusion parameters compared to conventional gSlider reconstruction. Results also demonstrate that the proposed method has advantages relative to denoising methods based on low-rank matrix modeling. A theoretical analysis of the trade-off between spatial resolution and SNR suggests that the proposed approach has high efficiency. CONCLUSIONS: The combination of gSlider-SMS with advanced regularized reconstruction enables high-resolution quantitative diffusion MRI from a relatively fast acquisition.

Phase-matched virtual coil reconstruction for highly accelerated diffusion echo-planar imaging.

PURPOSE: To propose a virtual coil (VC) acquisition/reconstruction framework to improve highly accelerated single-shot EPI (SS-EPI) and generalized slice dithered enhanced resolution (gSlider) acquisition in high-resolution diffusion imaging (DI). METHODS: For robust VC-GRAPPA reconstruction, a background phase correction scheme was developed to match the image phase of the reference data with the corrupted phase of the accelerated diffusion-weighted data, where the corrupted phase of the diffusion data varies from shot to shot. A Gy prewinding-blip was also added to the EPI acquisition, to create a shifted-ky sampling strategy that allows for better exploitation of VC concept in the reconstruction. To evaluate the performance of the proposed methods, 1.5 mm isotropic whole-brain SS-EPI and 860 µm isotropic whole-brain gSlider-EPI diffusion data were acquired at an acceleration of 8-9 fold. Conventional and VC-GRAPPA reconstructions were performed and compared, and corresponding g-factors were calculated. RESULTS: The proposed VC reconstruction substantially improves the image quality of both SS-EPI and gSlider-EPI, with reduced g-factor noise and reconstruction artifacts when compared to the conventional method. This has enabled high-quality low-noise diffusion imaging to be performed at 8-9 fold acceleration. CONCLUSIONS: The proposed VC acquisition/reconstruction framework improves the reconstruction of DI at high accelerations. The ability to now employ such high accelerations will allow DI with EPI at reduced distortion and faster scan time, which should be beneficial for many clinical and neuroscience applications.

Age-related alterations in axonal microstructure in the corpus callosum measured by high-gradient diffusion MRI.

Cerebral white matter exhibits age-related degenerative changes during the course of normal aging, including decreases in axon density and alterations in axonal structure. Noninvasive approaches to measure these microstructural alterations throughout the lifespan would be invaluable for understanding the substrate and regional variability of age-related white matter degeneration. Recent advances in diffusion magnetic resonance imaging (MRI) have leveraged high gradient strengths to increase sensitivity toward axonal size and density in the living human brain. Here, we examined the relationship between age and indices of axon diameter and packing density using high-gradient strength diffusion MRI in 36 healthy adults (aged 22-72) in well-defined central white matter tracts in the brain. A recently validated method for inferring the effective axonal compartment size and packing density from diffusion MRI measurements acquired with 300 mT/m maximum gradient strength was applied to the in vivo human brain to obtain indices of axon diameter and density in the corpus callosum, its sub-regions, and adjacent anterior and posterior fibers in the forceps minor and forceps major. The relationships between the axonal metrics, corpus callosum area and regional gray matter volume were also explored. Results revealed a significant increase in axon diameter index with advancing age in the whole corpus callosum. Similar analyses in sub-regions of the corpus callosum showed that age-related alterations in axon diameter index and axon density were most pronounced in the genu of the corpus callosum and relatively absent in the splenium, in keeping with findings from previous histological studies. The significance of these correlations was mirrored in the forceps minor and forceps major, consistent with previously reported decreases in FA in the forceps minor but not in the forceps major with age. Alterations in the axonal imaging metrics paralleled decreases in corpus callosum area and regional gray matter volume with age. Among older adults, results from cognitive testing suggested an association between larger effective compartment size in the corpus callosum, particularly within the genu of the corpus callosum, and lower scores on the Montreal Cognitive Assessment, largely driven by deficits in short-term memory. The current study suggests that high-gradient diffusion MRI may be sensitive to the axonal substrate of age-related white matter degeneration reflected in traditional DTI metrics and provides further evidence for regionally selective alterations in white matter microstructure with advancing age.

Validation of diffusion MRI estimates of compartment size and volume fraction in a biomimetic brain phantom using a human MRI scanner with 300 mT/m maximum gradient strength.

Diffusion microstructural imaging techniques have attracted great interest in the last decade due to their ability to quantify axon diameter and volume fraction in healthy and diseased human white matter. The estimates of compartment size and volume fraction continue to be debated, in part due to the lack of a gold standard for validation and quality control. In this work, we validate diffusion MRI estimates of compartment size and volume fraction using a novel textile axon ("taxon") phantom constructed from hollow polypropylene yarns with distinct intra- and extra-taxonal compartments to mimic white matter in the brain. We acquired a comprehensive set of diffusion MRI measurements in the phantom using multiple gradient directions, diffusion times and gradient strengths on a human MRI scanner equipped with maximum gradient strength (Gmax) of 300 mT/m. We obtained estimates of compartment size and restricted volume fraction through a straightforward extension of the AxCaliber/ActiveAx frameworks that enables estimation of mean compartment size in fiber bundles of arbitrary orientation. The voxel-wise taxon diameter estimates of 12.2 ±â€¯0.9 µm were close to the manufactured inner diameter of 11.8 ±â€¯1.2 µm with Gmax = 300 mT/m. The estimated restricted volume fraction demonstrated an expected decrease along the length of the fiber bundles in accordance with the known construction of the phantom. When Gmax was restricted to 80 mT/m, the taxon diameter was overestimated, and the estimates for taxon diameter and packing density showed greater uncertainty compared to data with Gmax = 300 mT/m. In conclusion, the compartment size and volume fraction estimates resulting from diffusion measurements on a human scanner were validated against ground truth in a phantom mimicking human white matter, providing confidence that this method can yield accurate estimates of parameters in simplified but realistic microstructural environments. Our work also demonstrates the importance of a biologically analogous phantom that can be applied to validate a variety of diffusion microstructural imaging methods in human scanners and be used for standardization of diffusion MRI protocols for neuroimaging research.

High-resolution in vivo diffusion imaging of the human brain with generalized slice dithered enhanced resolution: Simultaneous multislice (gSlider-SMS).

PURPOSE: To develop an efficient acquisition for high-resolution diffusion imaging and allow in vivo whole-brain acquisitions at 600- to 700-µm isotropic resolution. METHODS: We combine blipped-controlled aliasing in parallel imaging simultaneous multislice (SMS) with a novel slab radiofrequency (RF) encoding gSlider (generalized slice-dithered enhanced resolution) to form a signal-to-noise ratio-efficient volumetric simultaneous multislab acquisition. Here, multiple thin slabs are acquired simultaneously with controlled aliasing, and unaliased with parallel imaging. To achieve high resolution in the slice direction, the slab is volumetrically encoded using RF encoding with a scheme similar to Hadamard encoding. However, with gSlider, the RF-encoding bases are specifically designed to be highly independent and provide high image signal-to-noise ratio in each slab acquisition to enable self-navigation of the diffusion's phase corruption. Finally, the method is combined with zoomed imaging (while retaining whole-brain coverage) to facilitate low-distortion single-shot in-plane encoding with echo-planar imaging at high resolution. RESULTS: A 10-slices-per-shot gSlider-SMS acquisition was used to acquire whole-brain data at 660 and 760 µm isotropic resolution with b-values of 1500 and 1800 s/mm2 , respectively. Data were acquired on the Connectome 3 Tesla scanner with 64-channel head coil. High-quality data with excellent contrast were achieved at these resolutions, which enable the visualization of fine-scale structures. CONCLUSIONS: The gSlider-SMS approach provides a new, efficient way to acquire high-resolution diffusion data. Magn Reson Med 79:141-151, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Motion-robust sub-millimeter isotropic diffusion imaging through motion corrected generalized slice dithered enhanced resolution (MC-gSlider) acquisition.

PURPOSE: To develop an efficient MR technique for ultra-high resolution diffusion MRI (dMRI) in the presence of motion. METHODS: gSlider is an SNR-efficient high-resolution dMRI acquisition technique. However, subject motion is inevitable during a prolonged scan for high spatial resolution, leading to potential image artifacts and blurring. In this study, an integrated technique termed Motion Corrected gSlider (MC-gSlider) is proposed to obtain high-quality, high-resolution dMRI in the presence of large in-plane and through-plane motion. A motion-aware reconstruction with spatially adaptive regularization is developed to optimize the conditioning of the image reconstruction under difficult through-plane motion cases. In addition, an approach for intra-volume motion estimation and correction is proposed to achieve motion correction at high temporal resolution. RESULTS: Theoretical SNR and resolution analysis validated the efficiency of MC-gSlider with regularization, and aided in selection of reconstruction parameters. Simulations and in vivo experiments further demonstrated the ability of MC-gSlider to mitigate motion artifacts and recover detailed brain structures for dMRI at 860 µm isotropic resolution in the presence of motion with various ranges. CONCLUSION: MC-gSlider provides motion-robust, high-resolution dMRI with a temporal motion correction sensitivity of 2 s, allowing for the recovery of fine detailed brain structures in the presence of large subject movements.

In vivo characterization of cortical and white matter neuroaxonal pathology in early multiple sclerosis.

Neuroaxonal pathology is a main determinant of disease progression in multiple sclerosis; however, its underlying pathophysiological mechanisms, including its link to inflammatory demyelination and temporal occurrence in the disease course are still unknown. We used ultra-high field (7 T), ultra-high gradient strength diffusion and T1/T2-weighted myelin-sensitive magnetic resonance imaging to characterize microstructural changes in myelin and neuroaxonal integrity in the cortex and white matter in early stage multiple sclerosis, their distribution in lesional and normal-appearing tissue, and their correlations with neurological disability. Twenty-six early stage multiple sclerosis subjects (disease duration ≤5 years) and 24 age-matched healthy controls underwent 7 T T2*-weighted imaging for cortical lesion segmentation and 3 T T1/T2-weighted myelin-sensitive imaging and neurite orientation dispersion and density imaging for assessing microstructural myelin, axonal and dendrite integrity in lesional and normal-appearing tissue of the cortex and the white matter. Conventional mean diffusivity and fractional anisotropy metrics were also assessed for comparison. Cortical lesions were identified in 92% of early multiple sclerosis subjects and they were characterized by lower intracellular volume fraction (P = 0.015 by paired t-test), lower myelin-sensitive contrast (P = 0.030 by related-samples Wilcoxon signed-rank test) and higher mean diffusivity (P = 0.022 by related-samples Wilcoxon signed-rank test) relative to the contralateral normal-appearing cortex. Similar findings were observed in white matter lesions relative to normal-appearing white matter (all P < 0.001), accompanied by an increased orientation dispersion (P < 0.001 by paired t-test) and lower fractional anisotropy (P < 0.001 by related-samples Wilcoxon signed-rank test) suggestive of less coherent underlying fibre orientation. Additionally, the normal-appearing white matter in multiple sclerosis subjects had diffusely lower intracellular volume fractions than the white matter in controls (P = 0.029 by unpaired t-test). Cortical thickness did not differ significantly between multiple sclerosis subjects and controls. Higher orientation dispersion in the left primary motor-somatosensory cortex was associated with increased Expanded Disability Status Scale scores in surface-based general linear modelling (P < 0.05). Microstructural pathology was frequent in early multiple sclerosis, and present mainly focally in cortical lesions, whereas more diffusely in white matter. These results suggest early demyelination with loss of cells and/or cell volumes in cortical and white matter lesions, with additional axonal dispersion in white matter lesions. In the cortex, focal lesion changes might precede diffuse atrophy with cortical thinning. Findings in the normal-appearing white matter reveal early axonal pathology outside inflammatory demyelinating lesions.

HIgh b-value and high Resolution Integrated Diffusion (HIBRID) imaging.

The parameter selection for diffusion MRI experiments is dominated by the "k-q tradeoff" whereby the Signal to Noise Ratio (SNR) of the images is traded for either high spatial resolution (determined by the maximum k-value collected) or high diffusion sensitivity (effected by b-value or the q vector) but usually not both. Furthermore, different brain regions (such as gray matter and white matter) likely require different tradeoffs between these parameters due to the size of the structures to be visualized or the length-scale of the microstructure being probed. In this case, it might be advantageous to combine information from two scans - a scan with high q but low k (high angular resolution in diffusion but low spatial resolution in the image domain) to provide maximal information about white matter fiber crossing, and one low q but high k (low angular resolution but high spatial resolution) for probing the cortex. In this study, we propose a method, termed HIgh b-value and high Resolution Integrated Diffusion (HIBRID) imaging, for acquiring and combining the information from these two complementary types of scan with the goal of studying diffusion in the cortex without compromising white matter fiber information. The white-gray boundary and pial surface obtained from anatomical scans are incorporated as prior information to guide the fusion. We study the complementary advantages of the fused datasets, and assess the quality of the HIBRID data compared to either alone.