Water Diffusion in the Live Human Brain is Gaussian at the Mesoscale.

Imaging the live human brain at the mesoscopic scale is a desideratum in basic and clinical neurosciences. Despite the promise of diffusion MRI, the lack of an accurate model relating the measured signal and the associated microstructure has hampered its success. The widely used diffusion tensor MRI (DTI) model assumes an anisotropic Gaussian diffusion process in each voxel, but lacks the ability to capture intravoxel heterogeneity. This study explores the extension of the DTI model to mesoscopic length scales by use of the diffusion tensor distribution (DTD) model, which assumes a Gaussian diffusion process in each subvoxel. DTD MRI has shown promise in addressing some limitations of DTI, particularly in distinguishing among different types of brain cancers and elucidating multiple fiber populations within a voxel. However, its validity in live brain tissue has never been established. Here, multiple diffusion-encoded (MDE) data were acquired in the living human brain using a 3 Tesla MRI scanner with large diffusion weighting factors. Two different diffusion times (Δ = 37, 74 ms) were employed, with other scanning parameters fixed to assess signal decay differences. In vivo diffusion-weighted signals in gray and white matter were nearly identical at the two diffusion times. Fitting the signals to the DTD model yielded indistinguishable results, except in the cerebrospinal fluid (CSF)-filled voxels likely due to pulsatile flow. Overall, the study supports the time invariance of water diffusion at the mesoscopic scale in live brain parenchyma, extending the validity of the anisotropic Gaussian diffusion model in clinical brain imaging.

Eddy current-induced artifact correction in high b-value ex vivo human brain diffusion MRI with dynamic field monitoring.

PURPOSE: To investigate whether spatiotemporal magnetic field monitoring can correct pronounced eddy current-induced artifacts incurred by strong diffusion-sensitizing gradients up to 300 mT/m used in high b-value diffusion-weighted (DW) EPI. METHODS: A dynamic field camera equipped with 16 1 H NMR field probes was first used to characterize field perturbations caused by residual eddy currents from diffusion gradients waveforms in a 3D multi-shot EPI sequence on a 3T Connectom scanner for different gradient strengths (up to 300 mT/m), diffusion directions, and shots. The efficacy of dynamic field monitoring-based image reconstruction was demonstrated on high-gradient strength, submillimeter resolution whole-brain ex vivo diffusion MRI. A 3D multi-shot image reconstruction framework was developed that incorporated the nonlinear phase evolution measured with the dynamic field camera. RESULTS: Phase perturbations in the readout induced by residual eddy currents from strong diffusion gradients are highly nonlinear in space and time, vary among diffusion directions, and interfere significantly with the image encoding gradients, changing the k-space trajectory. During the readout, phase modulations between odd and even EPI echoes become non-static and diffusion encoding direction-dependent. Superior reduction of ghosting and geometric distortion was achieved with dynamic field monitoring compared to ghosting reduction approaches such as navigator- and structured low-rank-based methods or MUSE followed by image-based distortion correction with the FSL tool "eddy." CONCLUSION: Strong eddy current artifacts characteristic of high-gradient strength DW-EPI can be well corrected with dynamic field monitoring-based image reconstruction.

Eddy current-induced artifacts correction in high gradient strength diffusion MRI with dynamic field monitoring: demonstration in ex vivo human brain imaging.

Purpose: To demonstrate the advantages of spatiotemporal magnetic field monitoring to correct eddy current-induced artifacts (ghosting and geometric distortions) in high gradient strength diffusion MRI (dMRI). Methods: A dynamic field camera with 16 NMR field probes was used to characterize eddy current fields induced from diffusion gradients for different gradients strengths (up to 300 mT/m), diffusion directions, and shots in a 3D multi-shot EPI sequence on a 3T Connectom scanner. The efficacy of dynamic field monitoring-based image reconstruction was demonstrated on high-resolution whole brain ex vivo dMRI. A 3D multi-shot image reconstruction framework was informed with the actual nonlinear phase evolution measured with the dynamic field camera, thereby accounting for high-order eddy currents fields on top of the image encoding gradients in the image formation model. Results: Eddy current fields from diffusion gradients at high gradient strength in a 3T Connectom scanner are highly nonlinear in space and time, inducing high-order spatial phase modulations between odd/even echoes and shots that are not static during the readout. Superior reduction of ghosting and geometric distortion was achieved with dynamic field monitoring compared to ghosting approaches such as navigator- and structured low-rank-based methods or MUSE, followed by image-based distortion correction with eddy. Improved dMRI analysis is demonstrated with diffusion tensor imaging and high-angular resolution diffusion imaging. Conclusion: Strong eddy current artifacts characteristic of high gradient strength dMRI can be well corrected with dynamic field monitoring-based image reconstruction, unlike the two-step approach consisting of ghosting correction followed by geometric distortion reduction with eddy.

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.

High-fidelity, high-spatial-resolution diffusion magnetic resonance imaging of ex vivo whole human brain at ultra-high gradient strength with structured low-rank echo-planar imaging ghost correction.

Diffusion magnetic resonance imaging (dMRI) of whole ex vivo human brain specimens enables three-dimensional (3D) mapping of structural connectivity at the mesoscopic scale, providing detailed evaluation of fiber architecture and tissue microstructure at a spatial resolution that is difficult to access in vivo. To account for the short T2 and low diffusivity of fixed tissue, ex vivo dMRI is often acquired using strong diffusion-sensitizing gradients and multishot/segmented 3D echo-planar imaging (EPI) sequences to achieve high spatial resolution. However, the combination of strong diffusion-sensitizing gradients and multishot/segmented EPI readout can result in pronounced ghosting artifacts incurred by nonlinear spatiotemporal variations in the magnetic field produced by eddy currents. Such ghosting artifacts cannot be corrected with conventional correction solutions and pose a significant roadblock to leveraging human MRI scanners with ultrahigh gradients for ex vivo whole-brain dMRI. Here, we show that ghosting-correction approaches that correct for either polarity-related ghosting or shot-to-shot variations in a separate manner are suboptimal for 3D multishot diffusion-weighted EPI experiments in fixed human brain specimens using strong diffusion-sensitizing gradients on the 3-T Connectom MRI scanner, resulting in orientationally biased dMRI estimates. We apply a recently developed advanced k-space reconstruction method based on structured low-rank matrix (SLM) modeling that handles both polarity-related ghosting and shot-to-shot variation simultaneously, to mitigate artifacts in high-angular resolution multishot dMRI data acquired in several fixed human brain specimens at 0.7-0.8-mm isotropic spatial resolution using b-values up to 10,000 s/mm2 and gradient strengths up to 280 mT/m. We demonstrate the improved mapping of diffusion tensor imaging and fiber orientation distribution functions in key neuroanatomical areas distributed across the whole brain using SLM-based EPI ghost correction compared with alternative techniques.

Slice-direction geometric distortion evaluation and correction with reversed slice-select gradient acquisitions.

Accurate spatial alignment of MRI data acquired across multiple contrasts in the same subject is often crucial for data analysis and interpretation, but can be challenging in the presence of geometric distortions that differ between acquisitions. It is well known that single-shot echo-planar imaging (EPI) acquisitions suffer from distortion in the phase-encoding direction due to B0 field inhomogeneities arising from tissue magnetic susceptibility differences and other sources, however there can be distortion in other encoding directions as well in the presence of strong field inhomogeneities. High-resolution ultrahigh-field MRI typically uses low bandwidth in the slice-encoding direction to acquire thin slices and, when combined with the pronounced B0 inhomogeneities, is prone to an additional geometric distortion in the slice direction as well. Here we demonstrate the presence of this slice distortion in high-resolution 7T EPI acquired with a novel pulse sequence allowing for the reversal of the slice-encoding gradient polarity that enables the acquisition of pairs of images with equal magnitudes of distortion in the slice direction but with opposing polarities. We also show that the slice-direction distortion can be corrected using gradient reversal-based method applying the same software used for conventional corrections of phase-encoding direction distortion.

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.

Imaging High-Risk Atherothrombosis Using a Novel Fibrin-Binding Positron Emission Tomography Probe.

BACKGROUND AND PURPOSE: High-risk atherosclerosis is an underlying cause of cardiovascular events, yet identifying the specific patient population at immediate risk is still challenging. Here, we used a rabbit model of atherosclerotic plaque rupture and human carotid endarterectomy specimens to describe the potential of molecular fibrin imaging as a tool to identify thrombotic plaques. METHODS: Atherosclerotic plaques in rabbits were induced using a high-cholesterol diet and aortic balloon injury (N=13). Pharmacological triggering was used in a group of rabbits (n=9) to induce plaque disruption. Animals were grouped into thrombotic and nonthrombotic plaque groups based on gross pathology (gold standard). All animals were injected with a novel fibrin-specific probe 68Ga-CM246 followed by positron emission tomography (PET)/magnetic resonance imaging 90 minutes later. 68Ga-CM246 was quantified on the PET images using tissue-to-background (back muscle) ratios and standardized uptake value. RESULTS: Both tissue-to-background (back muscle) ratios and standardized uptake value were significantly higher in the thrombotic versus nonthrombotic group (P<0.05). Ex vivo PET and autoradiography of the abdominal aorta correlated positively with in vivo PET measurements. Plaque disruption identified by 68Ga-CM246 PET agreed with gross pathology assessment (85%). In ex vivo surgical specimens obtained from patients undergoing elective carotid endarterectomy (N=12), 68Ga-CM246 showed significantly higher binding to carotid plaques compared to a D-cysteine nonbinding control probe. CONCLUSIONS: We demonstrated that molecular fibrin PET imaging using 68Ga-CM246 could be a useful tool to diagnose experimental and clinical atherothrombosis. Based on our initial results using human carotid plaque specimens, in vivo molecular imaging studies are warranted to test 68Ga-CM246 PET as a tool to stratify risk in atherosclerotic patients.

A 31-channel integrated "AC/DC" B0 shim and radiofrequency receive array coil for improved 7T MRI.

PURPOSE: To test an integrated "AC/DC" array approach at 7T, where B0 inhomogeneity poses an obstacle for functional imaging, diffusion-weighted MRI, MR spectroscopy, and other applications. METHODS: A close-fitting 7T 31-channel (31-ch) brain array was constructed and tested using combined Rx and ΔB0 shim channels driven by a set of rapidly switchable current amplifiers. The coil was compared to a shape-matched 31-ch reference receive-only array for RF safety, signal-to-noise ratio (SNR), and inter-element noise correlation. We characterize the coil array's ability to provide global and dynamic (slice-optimized) shimming using ΔB0 field maps and echo planar imaging (EPI) acquisitions. RESULTS: The SNR and average noise correlation were similar to the 31-ch reference array. Global and slice-optimized shimming provide 11% and 40% improvements respectively compared to baseline second-order spherical harmonic shimming. Birdcage transmit coil efficiency was similar for the reference and AC/DC array setups. CONCLUSION: Adding ΔB0 shim capability to a 31-ch 7T receive array can significantly boost 7T brain B0 homogeneity without sacrificing the array's rdiofrequency performance, potentially improving ultra-high field neuroimaging applications that are vulnerable to off-resonance effects.

Rapid simultaneous acquisition of macromolecular tissue volume, susceptibility, and relaxometry maps.

PURPOSE: A major obstacle to the clinical implementation of quantitative MR is the lengthy acquisition time required to derive multi-contrast parametric maps. We sought to reduce the acquisition time for QSM and macromolecular tissue volume by acquiring both contrasts simultaneously by leveraging their redundancies. The joint virtual coil concept with GRAPPA (JVC-GRAPPA) was applied to reduce acquisition time further. METHODS: Three adult volunteers were imaged on a 3 Tesla scanner using a multi-echo 3D GRE sequence acquired at 3 head orientations. Macromolecular tissue volume, QSM, R2∗ , T1 , and proton density maps were reconstructed. The same sequence (GRAPPA R = 4) was performed in subject 1 with a single head orientation for comparison. Fully sampled data was acquired in subject 2, from which retrospective undersampling was performed (R = 6 GRAPPA and R = 9 JVC-GRAPPA). Prospective undersampling was performed in subject 3 (R = 6 GRAPPA and R = 9 JVC-GRAPPA) using gradient blips to shift k-space sampling in later echoes. RESULTS: Subject 1's multi-orientation and single-orientation macromolecular tissue volume maps were not significantly different based on RMSE. For subject 2, the retrospectively undersampled JVC-GRAPPA and GRAPPA generated similar results as fully sampled data. This approach was validated with the prospectively undersampled images in subject 3. Using QSM, R2∗ , and macromolecular tissue volume, the contributions of myelin and iron content to susceptibility were estimated. CONCLUSION: We have developed a novel strategy to simultaneously acquire data for the reconstruction of 5 intrinsically coregistered 1-mm isotropic resolution multi-parametric maps, with a scan time of 6 min using JVC-GRAPPA.

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.

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.

A 48-channel receive array coil for mesoscopic diffusion-weighted MRI of ex vivo human brain on the 3 T connectome scanner.

In vivo diffusion-weighted magnetic resonance imaging is limited in signal-to-noise-ratio (SNR) and acquisition time, which constrains spatial resolution to the macroscale regime. Ex vivo imaging, which allows for arbitrarily long scan times, is critical for exploring human brain structure in the mesoscale regime without loss of SNR. Standard head array coils designed for patients are sub-optimal for imaging ex vivo whole brain specimens. The goal of this work was to design and construct a 48-channel ex vivo whole brain array coil for high-resolution and high b-value diffusion-weighted imaging on a 3T Connectome scanner. The coil was validated with bench measurements and characterized by imaging metrics on an agar brain phantom and an ex vivo human brain sample. The two-segment coil former was constructed for a close fit to a whole human brain, with small receive elements distributed over the entire brain. Imaging tests including SNR and G-factor maps were compared to a 64-channel head coil designed for in vivo use. There was a 2.9-fold increase in SNR in the peripheral cortex and a 1.3-fold gain in the center when compared to the 64-channel head coil. The 48-channel ex vivo whole brain coil also decreases noise amplification in highly parallel imaging, allowing acceleration factors of approximately one unit higher for a given noise amplification level. The acquired diffusion-weighted images in a whole ex vivo brain specimen demonstrate the applicability and advantage of the developed coil for high-resolution and high b-value diffusion-weighted ex vivo brain MRI studies.

Ultra-high spatial resolution BOLD fMRI in humans using combined segmented-accelerated VFA-FLEET with a recursive RF pulse design.

PURPOSE: To alleviate the spatial encoding limitations of single-shot echo-planar imaging (EPI) by developing multi-shot segmented EPI for ultra-high-resolution functional MRI (fMRI) with reduced ghosting artifacts from subject motion and respiration. THEORY AND METHODS: Segmented EPI can reduce readout duration and reduce acceleration factors, however, the time elapsed between segment acquisitions (on the order of seconds) can result in intermittent ghosting, limiting its use for fMRI. Here, "FLEET" segment ordering, where segments are looped over before slices, was combined with a variable flip angle progression (VFA-FLEET) to improve inter-segment fidelity and maximize signal for fMRI. Scaling a sinc pulse's flip angle for each segment (VFA-FLEET-Sinc) produced inconsistent slice profiles and ghosting, therefore, a recursive Shinnar-Le Roux (SLR) radiofrequency (RF) pulse design was developed (VFA-FLEET-SLR) to generate unique pulses for every segment that together produce consistent slice profiles and signals. RESULTS: The temporal stability of VFA-FLEET-SLR was compared against conventional-segmented EPI and VFA-FLEET-Sinc at 3T and 7T. VFA-FLEET-SLR showed reductions in both intermittent and stable ghosting compared to conventional-segmented and VFA-FLEET-Sinc, resulting in improved image quality with a minor trade-off in temporal SNR. Combining VFA-FLEET-SLR with acceleration, we achieved a 0.6-mm isotropic acquisition at 7T, without zoomed imaging or partial Fourier, demonstrating reliable detection of blood oxygenation level-dependent (BOLD) responses to a visual stimulus. To counteract the increased repetition time from segmentation, simultaneous multi-slice VFA-FLEET-SLR was demonstrated using RF-encoded controlled aliasing. CONCLUSIONS: VFA-FLEET with a recursive RF pulse design supports acquisitions with low levels of artifact and spatial blur, enabling fMRI at previously inaccessible spatial resolutions with a "full-brain" field of view.

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.

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.

Axonal damage in the optic radiation assessed by white matter tract integrity metrics is associated with retinal thinning in multiple sclerosis.

INTRODUCTION: White matter damage in the visual pathway is common in multiple sclerosis (MS) and is associated with retinal thinning, although the underlying mechanism of association remains unclear. The goal of this work was to evaluate the presence and extent of white matter tract integrity (WMTI) alterations in the optic radiation (OR) in people with MS and to investigate the association between WMTI metrics and retinal thinning in the eyes of MS patients without a history of optic neuritis (ON) as measured by optical coherence tomography (OCT). We hypothesized that WMTI metrics would reflect axonal damage that occurs in the OR in MS, and that axonal alterations revealed by WMTI would be associated with retinal thinning. METHODS: Twenty-nine MS patients without previous ON in at least one eye and twenty-nine age-matched healthy controls (HC) were scanned on a dedicated high-gradient 3-Tesla MRI scanner with 300 mT/m maximum gradient strength using a multi-shell diffusion MRI protocol (b = 800, 1500, 2400 s/mm2). The patients were divided into two subgroups according to history without ON (N = 18) or with ON in one eye (N = 11). Diffusion tensor imaging (DTI) metrics and WMTI metrics derived from diffusion kurtosis imaging were assessed in normal-appearing white matter (NAWM) of the OR and in focal lesions. Retinal thickness in the eyes of MS patients was measured by OCT. Student's t-test was used to assess group differences between MRI metrics. Linear regression was used to study the relationship between OCT metrics, including retinal nerve fiber layer (RNFL) and combined ganglion cell and inner plexiform layer thickness (GCL/IPL), visual acuity measures and DTI and WMTI metrics. RESULTS: OR NAWM in MS showed significantly decreased axonal water fraction (AWF) compared to HC (0.36 vs 0.39, p < 0.001), with similar trends observed in AWF of lesions compared to NAWM (0.27 vs 0.36, p < 0.001). Fractional anisotropy (FA) was lower in OR NAWM of MS patients compared to HC (0.49 vs 0.52, p < 0.001). In patients without ON, AWF was the only diffusion MRI metric that was significantly associated with average RNFL (r = 0.68, p = 0.005), adjusting for age, sex and disease duration and correcting for multiple comparisons. Of all the DTI and WMTI metrics, AWF was the strongest and most significant predictor of average RNFL thickness in MS patients without ON. There was no significant correlation between visual acuity scores and DTI or WMTI metrics after correction for multiple comparisons. CONCLUSION: Axonal damage may be the substrate of previously observed DTI alterations in the OR, as supported by the significant reduction in AWF within both NAWM and lesions of the OR in MS. Our results support the concept that axonal damage is widespread throughout the visual pathway in MS and may be mediated through trans-synaptic degeneration.

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.