High-gradient diffusion MRI reveals distinct estimates of axon diameter index within different white matter tracts in the in vivo human brain

Improving MR axon radius estimation in human white matter using spiral acquisition and field monitoring

PURPOSE

To compare MR axon radius estimation in human white matter using a multiband spiral sequence combined with field monitoring to the current state-of-the-art echo-planar imaging (EPI)-based approach.

METHODS

A custom multiband spiral sequence was used for diffusion-weighted imaging at ultra-highb $$ b $$ -values. Field monitoring and higher order image reconstruction were employed to greatly reduce artifacts in spiral images. Diffusion weighting parameters were chosen to match a state-of-the art EPI-based axon radius mapping protocol. The spiral approach was compared to the EPI approach by comparing the image signal-to-noise ratio (SNR) and performing a test-retest study to assess the respective variability and repeatability of axon radius mapping. Effective axon radius estimates were compared over white matter voxels and along the left corticospinal tract.

RESULTS

Increased SNR and reduced artifacts in spiral images led to reduced variability in resulting axon radius maps, especially in low-SNR regions. Test-retest variability was reduced by a factor of approximately 1.5 using the spiral approach. Reduced repeatability due to significant bias was found for some subjects in both spiral and EPI approaches, and attributed to scanner instability, pointing to a previously unknown limitation of the state-of-the-art approach.

CONCLUSION

Combining spiral readouts with field monitoring improved mapping of the effective axon radius compared to the conventional EPI approach.

Comparisons of MR and EM inferred tissue microstructure properties using a human autopsy corpus callosum sample

Degeneration of white matter (WM) microstructure in the central nervous system is characteristic of many neurodegenerative conditions. Previous research indicates that axonal degeneration visible in ex vivo electron microscopy (EM) photomicrographs precede the onset of clinical symptoms. Measuring WM microstructural features, such as axon diameter and packing fraction, currently require these highly invasive methods of analysis and it is therefore of great importance to develop methods for in vivo measurements. Diffusion weighted Magnetic Resonance Imaging (MRI) is a non-invasive method which can be used in conjunction with temporal diffusion spectroscopy (TDS) and an oscillating gradient spin echo (OGSE) pulse sequence to probe micron-scale structures within neural tissue. The current experiment aims to compare axon diameter measurements, mean effective axon diameter (AxD¯), and packing fractions calculated from EM histopathological analysis and inferred values from MR images. Mathematical models of axon diameters used for analysis include the ActiveAx Frequency-Dependent Extra-Axonal Diffusion (AAD) model and the AxCaliber Frequency-Dependent Extra-Axonal Diffusion (ACD) model using ROI (Region of Interest) based analysis (RBA) and voxel-based analysis (VBA), respectively. Overall, it was observed that MRI inferred WM microstructural parameters overestimate those calculated from EM. This may be attributable to tissue shrinkage during EM dehydration, the sensitivity of MR pulse sequences to larger diameter axons, and/or inaccurate model assumptions. The results of the current study provide a means to characterize the precision and accuracy of RBA-ACD and VBA-AAD OGSE-TDS and highlight the need for further research investigating the relationship between ex vivo MRI and EM, with the goal of reaching in vivo MRI.

The use of 7T MRI in multiple sclerosis: review and consensus statement from the North American Imaging in Multiple Sclerosis Cooperative

The use of ultra-high-field 7-Tesla (7T) MRI in multiple sclerosis (MS) research has grown significantly over the past two decades. With recent regulatory approvals of 7T scanners for clinical use in 2017 and 2020, the use of this technology for routine care is poised to continue to increase in the coming years. In this context, the North American Imaging in MS Cooperative (NAIMS) convened a workshop in February 2023 to review the previous and current use of 7T technology for MS research and potential future research and clinical applications. In this workshop, experts were tasked with reviewing the current literature and proposing a series of consensus statements, which were reviewed and approved by the NAIMS. In this review and consensus paper, we provide background on the use of 7T MRI in MS research, highlighting this technology's promise for identification and quantification of aspects of MS pathology that are more difficult to visualize with lower-field MRI, such as grey matter lesions, paramagnetic rim lesions, leptomeningeal enhancement and the central vein sign. We also review the promise of 7T MRI to study metabolic and functional changes to the brain in MS. The NAIMS provides a series of consensus statements regarding what is currently known about the use of 7T MRI in MS, and additional statements intended to provide guidance as to what work is necessary going forward to accelerate 7T MRI research in MS and translate this technology for use in clinical practice and clinical trials. This includes guidance on technical development, proposals for a universal acquisition protocol and suggestions for research geared towards assessing the utility of 7T MRI to improve MS diagnostics, prognostics and therapeutic efficacy monitoring. The NAIMS expects that this article will provide a roadmap for future use of 7T MRI in MS.

Testing retrogenesis and physiological explanations for tract-wise white matter aging: links to developmental order, fibre calibre, and vascularization

To understand the consistently observed spatial distribution of white-matter (WM) aging, developmentally driven theories termed "retrogenesis" have gained traction, positing that the order of WM tract development predicts the order of declines. Regions that develop first are expected to deteriorate the last, i.e. "last-in-first-out". Alternatively, regions which develop most rapidly may also decline most rapidly in aging, or "gains-predict-loss". The validity of such theories remains uncertain, in part due to lack of clarity on the definition of developmental order. Importantly, our recent findings suggest that WM aging is also associated with physiological parameters such as perfusion, which may be linked to fibre metabolic need, which in turn varies with fibre size. Here we address the extent to which the degree of WM aging is determined by development trajectory (i.e. retrogenesis) and/or by physiological state. We obtained microstructural and perfusion measures using data from the Human Connectome Project in Aging (HCP-A), complemented by a meta-analysis involving maps of fibre calibre and macrovascular volume. Our results suggest that (1) while tracts that appear last or finish myelinating first in development display the slowest aging, the pattern of aging is not fully explained by retrogenesis; in fact, time courses of tract emergence and myelination give rise to opposite associations with WM decline; (2) tracts that appear earlier also have higher mean axon calibre and are also associated with lower degrees of WM microstructural aging; (3) such tracts also tend to exhibit relatively sustained CBF with a higher rate of lengthening of the arterial transit times (ATT), suggestive of collateral blood supply. These findings were also sex dependent in a tract-specific manner. Future work will investigate whether these are ultimately influenced by each tract's metabolic demand and the role of macrovascular collateral flow.

A practical guide for combining functional regions of interest and white matter bundles

Diffusion-weighted imaging (DWI) is the primary method to investigate macro- and microstructure of neural white matter in vivo. DWI can be used to identify and characterize individual-specific white matter bundles, enabling precise analyses on hypothesis-driven connections in the brain and bridging the relationships between brain structure, function, and behavior. However, cortical endpoints of bundles may span larger areas than what a researcher is interested in, challenging presumptions that bundles are specifically tied to certain brain functions. Functional MRI (fMRI) can be integrated to further refine bundles such that they are restricted to functionally-defined cortical regions. Analyzing properties of these Functional Sub-Bundles (FSuB) increases precision and interpretability of results when studying neural connections supporting specific tasks. Several parameters of DWI and fMRI analyses, ranging from data acquisition to processing, can impact the efficacy of integrating functional and diffusion MRI. Here, we discuss the applications of the FSuB approach, suggest best practices for acquiring and processing neuroimaging data towards this end, and introduce the FSuB-Extractor, a flexible open-source software for creating FSuBs. We demonstrate our processing code and the FSuB-Extractor on an openly-available dataset, the Natural Scenes Dataset.

Age-related alterations in human cortical microstructure across the lifespan: Insights from high-gradient diffusion MRI

The human brain undergoes age-related microstructural alterations across the lifespan. Soma and Neurite Density Imaging (SANDI), a novel biophysical model of diffusion MRI, provides estimates of cell body (soma) radius and density, and neurite density in gray matter. The goal of this cross-sectional study was to assess the sensitivity of high-gradient diffusion MRI toward age-related alterations in cortical microstructure across the adult lifespan using SANDI. Seventy-two cognitively unimpaired healthy subjects (ages 19-85 years; 40 females) were scanned on the 3T Connectome MRI scanner with a maximum gradient strength of 300mT/m using a multi-shell diffusion MRI protocol incorporating 8 b-values and diffusion time of 19 ms. Intra-soma signal fraction obtained from SANDI model-fitting to the data was strongly correlated with age in all major cortical lobes (r = -0.69 to -0.60, FDR-p < 0.001). Intra-soma signal fraction (r = 0.48-0.63, FDR-p < 0.001) and soma radius (r = 0.28-0.40, FDR-p < 0.04) were significantly correlated with cortical volume in the prefrontal cortex, frontal, parietal, and temporal lobes. The strength of the relationship between SANDI metrics and age was greater than or comparable to the relationship between cortical volume and age across the cortical regions, particularly in the occipital lobe and anterior cingulate gyrus. In contrast to the SANDI metrics, all associations between diffusion tensor imaging (DTI) and diffusion kurtosis imaging metrics and age were low to moderate. These results suggest that high-gradient diffusion MRI may be more sensitive to underlying substrates of neurodegeneration in the aging brain than DTI and traditional macroscopic measures of neurodegeneration such as cortical volume and thickness.

DIMOND: DIffusion Model OptimizatioN with Deep Learning

Diffusion magnetic resonance imaging is an important tool for mapping tissue microstructure and structural connectivity non-invasively in the in vivo human brain. Numerous diffusion signal models are proposed to quantify microstructural properties. Nonetheless, accurate estimation of model parameters is computationally expensive and impeded by image noise. Supervised deep learning-based estimation approaches exhibit efficiency and superior performance but require additional training data and may be not generalizable. A new DIffusion Model OptimizatioN framework using physics-informed and self-supervised Deep learning entitled "DIMOND" is proposed to address this problem. DIMOND employs a neural network to map input image data to model parameters and optimizes the network by minimizing the difference between the input acquired data and synthetic data generated via the diffusion model parametrized by network outputs. DIMOND produces accurate diffusion tensor imaging results and is generalizable across subjects and datasets. Moreover, DIMOND outperforms conventional methods for fitting sophisticated microstructural models including the kurtosis and NODDI model. Importantly, DIMOND reduces NODDI model fitting time from hours to minutes, or seconds by leveraging transfer learning. In summary, the self-supervised manner, high efficacy, and efficiency of DIMOND increase the practical feasibility and adoption of microstructure and connectivity mapping in clinical and neuroscientific applications.

Biophysical Principles and Computational Modeling of Deep Brain Stimulation

BACKGROUND

Deep brain stimulation (DBS) has revolutionized the treatment of neurological disorders, yet the mechanisms of DBS are still under investigation. Computational models are important in silico tools for elucidating these underlying principles and potentially for personalizing DBS therapy to individual patients. The basic principles underlying neurostimulation computational models, however, are not well known in the clinical neuromodulation community.

OBJECTIVE

In this study, we present a tutorial on the derivation of computational models of DBS and outline the biophysical contributions of electrodes, stimulation parameters, and tissue substrates to the effects of DBS.

RESULTS

Given that many aspects of DBS are difficult to characterize experimentally, computational models have played an important role in understanding how material, size, shape, and contact segmentation influence device biocompatibility, energy efficiency, the spatial spread of the electric field, and the specificity of neural activation. Neural activation is dictated by stimulation parameters including frequency, current vs voltage control, amplitude, pulse width, polarity configurations, and waveform. These parameters also affect the potential for tissue damage, energy efficiency, the spatial spread of the electric field, and the specificity of neural activation. Activation of the neural substrate also is influenced by the encapsulation layer surrounding the electrode, the conductivity of the surrounding tissue, and the size and orientation of white matter fibers. These properties modulate the effects of the electric field and determine the ultimate therapeutic response.

CONCLUSION

This article describes biophysical principles that are useful for understanding the mechanisms of neurostimulation.

The effects of axonal beading and undulation on axonal diameter estimation from diffusion MRI: Insights from simulations in human axons segmented from three-dimensional electron microscopy

The increasing availability of high-performance gradient systems in human MRI scanners has generated great interest in diffusion microstructural imaging applications such as axonal diameter mapping. Practically, sensitivity to axon diameter in diffusion MRI is attained at strong diffusion weightings b , where the deviation from the expected 1 / b scaling in white matter yields a finite transverse diffusivity, which is then translated into an axon diameter estimate. While axons are usually modeled as perfectly straight, impermeable cylinders, local variations in diameter (caliber variation or beading) and direction (undulation) are known to influence axonal diameter estimates and have been observed in microscopy data of human axons. In this study, we performed Monte Carlo simulations of diffusion in axons reconstructed from three-dimensional electron microscopy of a human temporal lobe specimen using simulated sequence parameters matched to the maximal gradient strength of the next-generation Connectome 2.0 human MRI scanner ( ≲ 500 mT/m). We show that axon diameter estimation is accurate for nonbeaded, nonundulating fibers; however, in fibers with caliber variations and undulations, the axon diameter is heavily underestimated due to caliber variations, and this effect overshadows the known overestimation of the axon diameter due to undulations. This unexpected underestimation may originate from variations in the coarse-grained axial diffusivity due to caliber variations. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.

Response to "Validity of Quantitative, Multi-Parametric MRI in the Diagnosis of Polyneuropathies"

Level of Evidence

5

Technical Efficacy Stage

1

Diffusion MRI with free gradient waveforms on a high-performance gradient system: Probing restriction and exchange in the human brain

The dependence of the diffusion MRI signal on the diffusion time carries signatures of restricted diffusion and exchange. Here we seek to highlight these signatures in the human brain by performing experiments using free gradient waveforms designed to be selectively sensitive to the two effects. We examine six healthy volunteers using both strong and ultra-strong gradients (80, 200 and 300 mT/m). In an experiment featuring a large set of 150 gradient waveforms with different sensitivities to restricted diffusion and exchange, our results reveal unique and different time-dependence signatures in grey and white matter. Grey matter was characterised by both restricted diffusion and exchange and white matter predominantly by restricted diffusion. Exchange in grey matter was at least twice as fast as in white matter, across all subjects and all gradient strengths. The cerebellar cortex featured relatively short exchange times (115 ms). Furthermore, we show that gradient waveforms with tailored designs can be used to map exchange in the human brain. We also assessed the feasibility of clinical applications of the method used in this work and found that the exchange-related contrast obtained with a 25-minute protocol at 300 mT/m was preserved in a 4-minute protocol at 300 mT/m and a 10-minute protocol at 80 mT/m. Our work underlines the utility of free waveforms for detecting time dependence signatures due to restricted diffusion and exchange in vivo, which may potentially serve as a tool for studying diseased tissue.

In vivo evidence of microstructural hypo-connectivity of brain white matter in 22q11.2 deletion syndrome

22q11.2 deletion syndrome, or 22q11.2DS, is a genetic syndrome associated with high rates of schizophrenia and autism spectrum disorders, in addition to widespread structural and functional abnormalities throughout the brain. Experimental animal models have identified neuronal connectivity deficits, e.g., decreased axonal length and complexity of axonal branching, as a primary mechanism underlying atypical brain development in 22q11.2DS. However, it is still unclear whether deficits in axonal morphology can also be observed in people with 22q11.2DS. Here, we provide an unparalleled in vivo characterization of white matter microstructure in participants with 22q11.2DS (12-15 years) and those undergoing typical development (8-18 years) using a customized magnetic resonance imaging scanner which is sensitive to axonal morphology. A rich array of diffusion MRI metrics are extracted to present microstructural profiles of typical and atypical white matter development, and provide new evidence of connectivity differences in individuals with 22q11.2DS. A recent, large-scale consortium study of 22q11.2DS identified higher diffusion anisotropy and reduced overall diffusion mobility of water as hallmark microstructural alterations of white matter in individuals across a wide age range (6-52 years). We observed similar findings across the white matter tracts included in this study, in addition to identifying deficits in axonal morphology. This, in combination with reduced tract volume measurements, supports the hypothesis that abnormal microstructural connectivity in 22q11.2DS may be mediated by densely packed axons with disproportionately small diameters. Our findings provide insight into the in vivo white matter phenotype of 22q11.2DS, and promote the continued investigation of shared features in neurodevelopmental and psychiatric disorders.

Single-subject electroencephalography measurement of interhemispheric transfer time for the in-vivo estimation of axonal morphology

Assessing axonal morphology in vivo opens new avenues for the combined study of brain structure and function. A novel approach has recently been introduced to estimate the morphology of axonal fibers from the combination of magnetic resonance imaging (MRI) data and electroencephalography (EEG) measures of the interhemispheric transfer time (IHTT). In the original study, the IHTT measures were computed from EEG data averaged across a group, leading to bias of the axonal morphology estimates. Here, we seek to estimate axonal morphology from individual measures of IHTT, obtained from EEG data acquired in a visual evoked potential experiment. Subject-specific IHTTs are computed in a data-driven framework with minimal a priori constraints, based on the maximal peak of neural responses to visual stimuli within periods of statistically significant evoked activity in the inverse solution space. The subject-specific IHTT estimates ranged from 8 to 29 ms except for one participant and the between-session variability was comparable to between-subject variability. The mean radius of the axonal radius distribution, computed from the IHTT estimates and the MRI data, ranged from 0 to 1.09 μm across subjects. The change in axonal g-ratio with axonal radius ranged from 0.62 to 0.81 μm-α . The single-subject measurement of the IHTT yields estimates of axonal morphology that are consistent with histological values. However, improvement of the repeatability of the IHTT estimates is required to improve the specificity of the single-subject axonal morphology estimates.

A Method for In-Vivo Mapping of Axonal Diameter Distributions in the Human Brain Using Diffusion-Based Axonal Spectrum Imaging (AxSI)

In this paper we demonstrate a generalized and simplified pipeline called axonal spectrum imaging (AxSI) for in-vivo estimation of axonal characteristics in the human brain. Whole-brain estimation of the axon diameter, in-vivo and non-invasively, across all fiber systems will allow exploring uncharted aspects of brain structure and function relations with emphasis on connectivity and connectome analysis. While axon diameter mapping is important in and of itself, its correlation with conduction velocity will allow, for the first time, the explorations of information transfer mechanisms within the brain. We demonstrate various well-known aspects of axonal morphometry (e.g., the corpus callosum axon diameter variation) as well as other aspects that are less explored (e.g., axon diameter-based separation of the superior longitudinal fasciculus into segments). Moreover, we have created an MNI based mean axon diameter map over the entire brain for a large cohort of subjects providing the reference basis for future studies exploring relation between axon properties, its connectome representation, and other functional and behavioral aspects of the brain.

Detection of grey matter microstructural substrates of neurodegeneration in multiple sclerosis

Multiple sclerosis features complex pathological changes in grey matter that begin early and eventually lead to diffuse atrophy. Novel approaches to image grey-matter microstructural alterations in vivo are highly sought after and would enable more sensitive monitoring of disease activity and progression. This cross-sectional study aimed to assess the sensitivity of high-gradient diffusion MRI for microstructural tissue damage in cortical and deep grey matter in people with multiple sclerosis and test the hypothesis that reduced cortical cell body density is associated with cortical and deep grey-matter volume loss. Forty-one people with multiple sclerosis (age 24-72, 14 females) and 37 age- and sex-matched healthy controls were scanned on a 3 T Connectom MRI scanner equipped with 300 mT/m gradients using a multi-shell diffusion MRI protocol. The soma and neurite density imaging model was fitted to high-gradient diffusion MRI data to obtain estimates of intra-neurite, intra-cellular and extra-cellular signal fractions and apparent soma radius. Cortical and deep grey-matter microstructural imaging metrics were compared between multiple sclerosis and healthy controls and correlated with grey-matter volume, clinical disability and cognitive outcomes. People with multiple sclerosis showed significant cortical and deep grey-matter volume loss compared with healthy controls. People with multiple sclerosis showed trends towards lower cortical intra-cellular signal fraction and significantly lower intra-cellular and higher extra-cellular signal fractions in deep grey matter, especially the thalamus and caudate, compared with healthy controls. Changes were most pronounced in progressive disease and correlated with the Expanded Disability Status Scale, but not the Symbol Digit Modalities Test. In multiple sclerosis, normalized thalamic volume was associated with thalamic microstructural imaging metrics. Whereas thalamic volume loss did not correlate with cortical volume loss, cortical microstructural imaging metrics were significantly associated with thalamic volume, and not with cortical volume. Compared with the short diffusion time (Δ = 19 ms) achievable on the Connectom scanner, at the longer diffusion time of Δ = 49 ms attainable on clinical scanners, multiple sclerosis-related changes in imaging metrics were generally less apparent with lower effect sizes in cortical and deep grey matter. Soma and neurite density imaging metrics obtained from high-gradient diffusion MRI data provide detailed grey-matter characterization beyond cortical and thalamic volumes and distinguish multiple sclerosis-related microstructural pathology from healthy controls. Cortical cell body density correlates with thalamic volume, appears sensitive to the microstructural substrate of neurodegeneration and reflects disability status in people with multiple sclerosis, becoming more pronounced as disability worsens.

Novel insights into axon diameter and myelin content in late childhood and adolescence

White matter microstructural development in late childhood and adolescence is driven predominantly by increasing axon density and myelin thickness. Ex vivo studies suggest that the increase in axon diameter drives developmental increases in axon density observed with pubertal onset. In this cross-sectional study, 50 typically developing participants aged 8-18 years were scanned using an ultra-strong gradient magnetic resonance imaging scanner. Microstructural properties, including apparent axon diameter $({d}_a)$, myelin content, and g-ratio, were estimated in regions of the corpus callosum. We observed age-related differences in ${d}_a$, myelin content, and g-ratio. In early puberty, males had larger ${d}_a$ in the splenium and lower myelin content in the genu and body of the corpus callosum, compared with females. Overall, this work provides novel insights into developmental, pubertal, and cognitive correlates of individual differences in apparent axon diameter and myelin content in the developing human brain.

The influence of axonal beading and undulation on axonal diameter mapping

We consider the effect of non-cylindrical axonal shape on axonal diameter mapping with diffusion MRI. Practical sensitivity to axon diameter is attained at strong diffusion weightings b , where the deviation from the 1 / b scaling yields the finite transverse diffusivity, which is then translated into axon diameter. While axons are usually modeled as perfectly straight, impermeable cylinders, the local variations in diameter (caliber variation or beading) and direction (undulation) have been observed in microscopy data of human axons. Here we quantify the influence of cellular-level features such as caliber variation and undulation on axon diameter estimation. For that, we simulate the diffusion MRI signal in realistic axons segmented from 3-dimensional electron microscopy of a human brain sample. We then create artificial fibers with the same features and tune the amplitude of their caliber variations and undulations. Numerical simulations of diffusion in fibers with such tunable features show that caliber variations and undulations result in under- and over-estimation of axon diameters, correspondingly; this bias can be as large as 100%. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.

Diffusion MRI with free gradient waveforms on a high-performance gradient system: Probing restriction and exchange in the human brain

The dependence of the diffusion MRI signal on the diffusion time carries signatures of restricted diffusion and exchange. Here we seek to highlight these signatures in the human brain by performing experiments using free gradient waveforms that are selectively sensitive to the two effects. We examine six healthy volunteers using both strong and ultra-strong gradients (80, 200 and 300 mT/m). In an experiment featuring a large set of gradient waveforms with different sensitivities to restricted diffusion and exchange (150 samples), our results reveal unique time-dependence signatures in grey and white matter, where the former is characterised by both restricted diffusion and exchange and the latter predominantly exhibits restricted diffusion. Furthermore, we show that gradient waveforms with independently varying sensitivities to restricted diffusion and exchange can be used to map exchange in the human brain. We consistently find that exchange in grey matter is at least twice as fast as in white matter, across all subjects and all gradient strengths. The shortest exchange times observed in this study were in the cerebellar cortex (115 ms). We also assess the feasibility of future clinical applications of the method used in this work, where we find that the grey-white matter exchange contrast obtained with a 25-minute 300 mT/m protocol is preserved by a 4-minute 300 mT/m and a 10-minute 80 mT/m protocol. Our work underlines the utility of free waveforms for detecting time-dependence signatures due to restricted diffusion and exchange in vivo, which may potentially serve as a tool for studying diseased tissue.

Temporal Diffusion Ratio (TDR) for imaging restricted diffusion: Optimisation and pre-clinical demonstration

Temporal Diffusion Ratio (TDR) is a recently proposed dMRI technique (Dell'Acqua et al., proc. ISMRM 2019) which provides contrast between areas with restricted diffusion and areas either without restricted diffusion or with length scales too small for characterisation. Hence, it has a potential for informing on pore sizes, in particular the presence of large axon diameters or other cellular structures. TDR employs the signal from two dMRI acquisitions obtained with the same, large, b-value but with different diffusion gradient waveforms. TDR is advantageous as it employs standard acquisition sequences, does not make any assumptions on the underlying tissue structure and does not require any model fitting, avoiding issues related to model degeneracy. This work for the first time introduces and optimises the TDR method in simulation for a range of different tissues and scanner constraints and validates it in a pre-clinical demonstration. We consider both substrates containing cylinders and spherical structures, representing cell soma in tissue. Our results show that contrasting an acquisition with short gradient duration, short diffusion time and high gradient strength with an acquisition with long gradient duration, long diffusion time and low gradient strength, maximises the TDR contrast for a wide range of pore configurations. Additionally, in the presence of Rician noise, computing TDR from a subset (50% or fewer) of the acquired diffusion gradients rather than the entire shell as proposed originally further improves the contrast. In the last part of the work the results are demonstrated experimentally on rat spinal cord. In line with simulations, the experimental data shows that optimised TDR improves the contrast compared to non-optimised TDR. Furthermore, we find a strong correlation between TDR and histology measurements of axon diameter. In conclusion, we find that TDR has great potential and is a very promising alternative (or potentially complement) to model-based approaches for informing on pore sizes and restricted diffusion in general.

Mapping Human Brain Pathways: Challenges and Opportunities in the Integration of Scales

The human brain is composed of a complex web of pathways. Diffusion magnetic resonance (MR) tractography is a neuroimaging technique that relies on the principle of diffusion to reconstruct brain pathways. Its tractography is broadly applicable to a range of problems as it is amenable for study in individuals of any age and from any species. However, it is well known that this technique can generate biologically implausible pathways, especially in regions of the brain where multiple fibers cross. This review highlights potential misconnections in two cortico-cortical association pathways with a focus on the aslant tract and inferior frontal occipital fasciculus. The lack of alternative methods to validate observations from diffusion MR tractography means there is a need to develop new integrative approaches to trace human brain pathways. This review discusses integrative approaches in neuroimaging, anatomical, and transcriptional variation as having much potential to trace the evolution of human brain pathways.

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.

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.

Short superficial white matter and aging: a longitudinal multi-site study of 1293 subjects and 2711 sessions

It is estimated that short association fibers running immediately beneath the cortex may make up as much as 60% of the total white matter volume. However, these have been understudied relative to the long-range association, projection, and commissural fibers of the brain. This is largely because of limitations of diffusion MRI fiber tractography, which is the primary methodology used to non-invasively study the white matter connections. Inspired by recent anatomical considerations and methodological improvements in superficial white matter (SWM) tractography, we aim to characterize changes in these fiber systems in cognitively normal aging, which provide insight into the biological foundation of age-related cognitive changes, and a better understanding of how age-related pathology differs from healthy aging. To do this, we used three large, longitudinal and cross-sectional datasets (N = 1293 subjects, 2711 sessions) to quantify microstructural features and length/volume features of several SWM systems. We find that axial, radial, and mean diffusivities show positive associations with age, while fractional anisotropy has negative associations with age in SWM throughout the entire brain. These associations were most pronounced in the frontal, temporal, and temporoparietal regions. Moreover, measures of SWM volume and length decrease with age in a heterogenous manner across the brain, with different rates of change in inter-gyri and intra-gyri SWM, and at slower rates than well-studied long-range white matter pathways. These features, and their variations with age, provide the background for characterizing normal aging, and, in combination with larger association pathways and gray matter microstructural features, may provide insight into fundamental mechanisms associated with aging and cognition.

Advancements in Gradient System Performance for Clinical and Research MRI

In magnetic resonance imaging (MRI), spatial field gradients are applied along each axis to encode the location of the nuclear spin in the frequency domain. During recent years, the development of new gradient technologies has been focused on the generation of stronger and faster gradient fields for imaging with higher spatial and temporal resolution. This benefits imaging methods, such as brain diffusion and functional MRI, and enables human imaging at ultra-high field MRI. In addition to improving gradient performance, new technologies have been presented to minimize peripheral nerve stimulation and gradient-related acoustic noise, both generated by the rapid switching of strong gradient fields. This review will provide a general background on the gradient system and update on the state-of-the-art gradient technology. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 1.

Diffusion MRI with pulsed and free gradient waveforms: Effects of restricted diffusion and exchange

Monitoring time dependence with diffusion MRI yields observables sensitive to compartment sizes (restricted diffusion) and membrane permeability (water exchange). However, restricted diffusion and exchange have opposite effects on the diffusion-weighted signal, which can lead to errors in parameter estimates. In this work, we propose a signal representation that incorporates the effects of both restricted diffusion and exchange up to second order in b-value and is compatible with gradient waveforms of arbitrary shape. The representation features mappings from a gradient waveform to two scalars that separately control the sensitivity to restriction and exchange. We demonstrate that these scalars span a two-dimensional space that can be used to choose waveforms that selectively probe restricted diffusion or exchange, eliminating the correlation between the two phenomena. We found that waveforms with specific but unconventional shapes provide an advantage over conventional pulsed and oscillating gradient acquisitions. We also show that parametrization of waveforms into a two-dimensional space can be used to understand protocols from other approaches that probe restricted diffusion and exchange. For example, we found that the variation of mixing time in filter-exchange imaging corresponds to variation of our exchange-weighting scalar at a fixed value of the restriction-weighting scalar. The proposed signal representation was evaluated using Monte Carlo simulations in identical parallel cylinders with hexagonal and random packing as well as parallel cylinders with gamma-distributed radii. Results showed that the approach is sensitive to sizes in the interval 4-12 μm and exchange rates in the simulated range of 0 to 20 s - 1 , but also that there is a sensitivity to the extracellular geometry. The presented theory constitutes a simple and intuitive description of how restricted diffusion and exchange influence the signal as well as a guide to protocol design capable of separating the two effects.

Axonal and myelin changes and their inter-relationship in the optic radiations in people with multiple sclerosis

Background

The imaging g-ratio, estimated from axonal volume fraction (AVF) and myelin volume fraction (MVF), is a novel biomarker of microstructural tissue integrity in multiple sclerosis (MS).

Objective

To assess axonal and myelin changes and their inter-relationship as measured by g-ratio in the optic radiations (OR) in people with MS (pwMS) with and without previous optic neuritis (ON) compared to healthy controls (HC).

Methods

Thirty pwMS and 17 HCs were scanned on a 3Tesla Connectom scanner. AVF and MVF, derived from a multi-shell diffusion protocol and macromolecular tissue volume, respectively, were measured in normal-appearing white matter (NAWM) and lesions within the OR and used to calculate imaging g-ratio.

Results

OR AVF and MVF were decreased in pwMS compared to HC, and in OR lesions compared to NAWM, whereas the g-ratio was not different. Compared to pwMS with previous ON, AVF and g-ratio tended to be higher in pwMS without prior ON. AVF and MVF, particularly in NAWM, were positively correlated with retinal thickness, which was more pronounced in pwMS with prior ON.

Conclusion

Axonal measures reflect microstructural tissue damage in the OR, particularly in the setting of remote ON, and correlate with established metrics of visual health in MS.

Associations between corpus callosum damage, clinical disability, and surface-based homologous inter-hemispheric connectivity in multiple sclerosis

Axonal damage in the corpus callosum is prevalent in multiple sclerosis (MS). Although callosal damage is associated with disrupted functional connectivity between hemispheres, it is unclear how this relates to cognitive and physical disability. We investigated this phenomenon using advanced measures of microstructural integrity in the corpus callosum and surface-based homologous inter-hemispheric connectivity (sHIC) in the cortex. We found that sHIC was significantly decreased in primary motor, somatosensory, visual, and temporal cortical areas in a group of 36 participants with MS (29 relapsing-remitting, 4 secondary progressive MS, and 3 primary-progressive MS) compared with 42 healthy controls (cluster level, p < 0.05). In participants with MS, global sHIC correlated with fractional anisotropy and restricted volume fraction in the posterior segment of the corpus callosum (r = 0.426, p = 0.013; r = 0.399, p = 0.020, respectively). Lower sHIC, particularly in somatomotor and posterior cortical areas, was associated with cognitive impairment and higher disability scores on the Expanded Disability Status Scale (EDSS). We demonstrated that higher levels of sHIC attenuated the effects of posterior callosal damage on physical disability and cognitive dysfunction, as measured by the EDSS and Brief Visuospatial Memory Test-Revised (interaction effect, p < 0.05). We also observed a positive association between global sHIC and years of education (r = 0.402, p = 0.018), supporting the phenomenon of "brain reserve" in MS. Our data suggest that preserved sHIC helps prevent cognitive and physical decline in MS.

Post mortem mapping of connectional anatomy for the validation of diffusion MRI

Diffusion MRI (dMRI) is a unique tool for the study of brain circuitry, as it allows us to image both the macroscopic trajectories and the microstructural properties of axon bundles in vivo. The Human Connectome Project ushered in an era of impressive advances in dMRI acquisition and analysis. As a result of these efforts, the quality of dMRI data that could be acquired in vivo improved substantially, and large collections of such data became widely available. Despite this progress, the main limitation of dMRI remains: it does not image axons directly, but only provides indirect measurements based on the diffusion of water molecules. Thus, it must be validated by methods that allow direct visualization of axons but that can only be performed in post mortem brain tissue. In this review, we discuss methods for validating the various features of connectional anatomy that are extracted from dMRI, both at the macro-scale (trajectories of axon bundles), and at micro-scale (axonal orientations and other microstructural properties). We present a range of validation tools, including anatomic tracer studies, Klingler's dissection, myelin stains, label-free optical imaging techniques, and others. We provide an overview of the basic principles of each technique, its limitations, and what it has taught us so far about the accuracy of different dMRI acquisition and analysis approaches.

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.

Soma and Neurite Density MRI (SANDI) of the in-vivo mouse brain and comparison with the Allen Brain Atlas

Diffusion MRI (dMRI) provides unique insights into the neural tissue milieu by probing interactions between diffusing molecules and tissue microstructure. Most dMRI techniques focus on white matter (WM) tissues, nevertheless, interest in gray matter characterizations is growing. The Soma and Neurite Density MRI (SANDI) methodology harnesses a model incorporating water diffusion in spherical objects (assumed to be associated with cell bodies) and in impermeable "sticks" (assumed to represent neurites), which potentially enables the characterization of cellular and neurite densities. Recognising the importance of rodents in animal models of development, aging, plasticity, and disease, we here employ SANDI for in-vivo preclinical imaging and provide a first validation of the methodology by comparing SANDI metrics with cellular density reflected by the Allen mouse brain atlas. SANDI was implemented on a 9.4T scanner equipped with a cryogenic coil, and in-vivo experiments were carried out on N = 6 mice. Pixelwise, ROI-based, and atlas comparisons were performed, magnitude vs. real-valued analyses were compared, and shorter acquisitions with reduced the number of b-value shells were investigated. Our findings reveal good reproducibility of the SANDI parameters, including the sphere and stick fractions, as well as sphere size (CoV < 7%, 12% and 3%, respectively). Additionally, we find a very good rank correlation between SANDI-driven sphere fraction and Allen mouse brain atlas contrast that represents cellular density. We conclude that SANDI is a viable preclinical MRI technique that can greatly contribute to research on brain tissue microstructure.

Long-Term Effects of Repeated Blast Exposure in United States Special Operations Forces Personnel: A Pilot Study Protocol

Emerging evidence suggests that repeated blast exposure (RBE) is associated with brain injury in military personnel. United States (U.S.) Special Operations Forces (SOF) personnel experience high rates of blast exposure during training and combat, but the effects of low-level RBE on brain structure and function in SOF have not been comprehensively characterized. Further, the pathophysiological link between RBE-related brain injuries and cognitive, behavioral, and physical symptoms has not been fully elucidated. We present a protocol for an observational pilot study, Long-Term Effects of Repeated Blast Exposure in U.S. SOF Personnel (ReBlast). In this exploratory study, 30 active-duty SOF personnel with RBE will participate in a comprehensive evaluation of: 1) brain network structure and function using Connectome magnetic resonance imaging (MRI) and 7 Tesla MRI; 2) neuroinflammation and tau deposition using positron emission tomography; 3) blood proteomics and metabolomics; 4) behavioral and physical symptoms using self-report measures; and 5) cognition using a battery of conventional and digitized assessments designed to detect subtle deficits in otherwise high-performing individuals. We will identify clinical, neuroimaging, and blood-based phenotypes that are associated with level of RBE, as measured by the Generalized Blast Exposure Value. Candidate biomarkers of RBE-related brain injury will inform the design of a subsequent study that will test a diagnostic assessment battery for detecting RBE-related brain injury. Ultimately, we anticipate that the ReBlast study will facilitate the development of interventions to optimize the brain health, quality of life, and battle readiness of U.S. SOF personnel.

Resolution and b value dependent Structural Connectome in ex vivo Mouse Brain

Diffusion magnetic resonance imaging has been widely used in both clinical and preclinical studies to characterize tissue microstructure and structural connectivity. The diffusion MRI protocol for the Human Connectome Project (HCP) has been developed and optimized to obtain high-quality, high-resolution diffusion MRI (dMRI) datasets. However, such efforts have not been fully explored in preclinical studies, especially for rodents. In this study, high quality dMRI datasets of mouse brains were acquired at 9.4T system from two vendors. In particular, we acquired a high-spatial resolution dMRI dataset (25 μm isotropic with 126 diffusion encoding directions), which we believe to be the highest spatial resolution yet obtained; and a high-angular resolution dMRI dataset (50 μm isotropic with 384 diffusion encoding directions), which we believe to be the highest angular resolution compared to the dMRI datasets at the microscopic resolution. We systematically investigated the effects of three important parameters that affect the final outcome of the connectome: b value (1000s/mm² to 8000 s/mm²), angular resolution (10 to 126), and spatial resolution (25 μm to 200 μm). The stability of tractography and connectome increase with the angular resolution, where more than 50 angles is necessary to achieve consistent results. The connectome and quantitative parameters derived from graph theory exhibit a linear relationship to the b value (R² > 0.99); a single-shell acquisition with b value of 3000 s/mm² shows comparable results to the multi-shell high angular resolution dataset. The dice coefficient decreases and both false positive rate and false negative rate gradually increase with coarser spatial resolution. Our study provides guidelines and foundations for exploration of tradeoffs among acquisition parameters for the structural connectome in ex vivo mouse brain.

An estimation of the absolute number of axons indicates that human cortical areas are sparsely connected

The tracts between cortical areas are conceived as playing a central role in cortical information processing, but their actual numbers have never been determined in humans. Here, we estimate the absolute number of axons linking cortical areas from a whole-cortex diffusion MRI (dMRI) connectome, calibrated using the histologically measured callosal fiber density. Median connectivity is estimated as approximately 6,200 axons between cortical areas within hemisphere and approximately 1,300 axons interhemispherically, with axons connecting functionally related areas surprisingly sparse. For example, we estimate that <5% of the axons in the trunk of the arcuate and superior longitudinal fasciculi connect Wernicke's and Broca's areas. These results suggest that detailed information is transmitted between cortical areas either via linkage of the dense local connections or via rare, extraordinarily privileged long-range connections.

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.

Using diffusion MRI data acquired with ultra-high gradient strength to improve tractography in routine-quality data

The development of scanners with ultra-high gradient strength, spearheaded by the Human Connectome Project, has led to dramatic improvements in the spatial, angular, and diffusion resolution that is feasible for in vivo diffusion MRI acquisitions. The improved quality of the data can be exploited to achieve higher accuracy in the inference of both microstructural and macrostructural anatomy. However, such high-quality data can only be acquired on a handful of Connectom MRI scanners worldwide, while remaining prohibitive in clinical settings because of the constraints imposed by hardware and scanning time. In this study, we first update the classical protocols for tractography-based, manual annotation of major white-matter pathways, to adapt them to the much greater volume and variability of the streamlines that can be produced from today’s state-of-the-art diffusion MRI data. We then use these protocols to annotate 42 major pathways manually in data from a Connectom scanner. Finally, we show that, when we use these manually annotated pathways as training data for global probabilistic tractography with anatomical neighborhood priors, we can perform highly accurate, automated reconstruction of the same pathways in much lower-quality, more widely available diffusion MRI data. The outcomes of this work include both a new, comprehensive atlas of WM pathways from Connectom data, and an updated version of our tractography toolbox, TRActs Constrained by UnderLying Anatomy (TRACULA), which is trained on data from this atlas. Both the atlas and TRACULA are distributed publicly as part of FreeSurfer. We present the first comprehensive comparison of TRACULA to the more conventional, multi-region-of-interest approach to automated tractography, and the first demonstration of training TRACULA on high-quality, Connectom data to benefit studies that use more modest acquisition protocols.

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.

MRI with ultrahigh field strength and high-performance gradients: challenges and opportunities for clinical neuroimaging at 7 T and beyond

Research in ultrahigh magnetic field strength combined with ultrahigh and ultrafast gradient technology has provided enormous gains in sensitivity, resolution, and contrast for neuroimaging. This article provides an overview of the technical advantages and challenges of performing clinical neuroimaging studies at ultrahigh magnetic field strength combined with ultrahigh and ultrafast gradient technology. Emerging clinical applications of 7-T MRI and state-of-the-art gradient systems equipped with up to 300 mT/m gradient strength are reviewed, and the impact and benefits of such advances to anatomical, structural and functional MRI are discussed in a variety of neurological conditions. Finally, an outlook and future directions for ultrahigh field MRI combined with ultrahigh and ultrafast gradient technology in neuroimaging are examined.

Dissecting whole-brain conduction delays through MRI microstructural measures

Network models based on structural connectivity have been increasingly used as the blueprint for large-scale simulations of the human brain. As the nodes of this network are distributed through the cortex and interconnected by white matter pathways with different characteristics, modeling the associated conduction delays becomes important. The goal of this study is to estimate and characterize these delays directly from the brain structure. To achieve this, we leveraged microstructural measures from a combination of advanced magnetic resonance imaging acquisitions and computed the main determinants of conduction velocity, namely axonal diameter and myelin content. Using the model proposed by Rushton, we used these measures to calculate the conduction velocity and estimated the associated delays using tractography. We observed that both the axonal diameter and conduction velocity distributions presented a rather constant trend across different connection lengths, with resulting delays that scale linearly with the connection length. Relying on insights from graph theory and Kuramoto simulations, our results support the approximation of constant conduction velocity but also show path- and region-specific differences.

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.

Bundle-Specific Axon Diameter Index as a New Contrast to Differentiate White Matter Tracts

In the central nervous system of primates, several pathways are characterized by different spectra of axon diameters. In vivo methods, based on diffusion-weighted magnetic resonance imaging, can provide axon diameter index estimates non-invasively. However, such methods report voxel-wise estimates, which vary from voxel-to-voxel for the same white matter bundle due to partial volume contributions from other pathways having different microstructure properties. Here, we propose a novel microstructure-informed tractography approach, COMMITAxSize, to resolve axon diameter index estimates at the streamline level, thus making the estimates invariant along trajectories. Compared to previously proposed voxel-wise methods, our formulation allows the estimation of a distinct axon diameter index value for each streamline, directly, furnishing a complementary measure to the existing calculation of the mean value along the bundle. We demonstrate the favourable performance of our approach comparing our estimates with existing histologically-derived measurements performed in the corpus callosum and the posterior limb of the internal capsule. Overall, our method provides a more robust estimation of the axon diameter index of pathways by jointly estimating the microstructure properties of the tissue and the macroscopic organisation of the white matter connectivity.

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.

The variability of MR axon radii estimates in the human white matter

The noninvasive quantification of axonal morphology is an exciting avenue for gaining understanding of the function and structure of the central nervous system. Accurate non-invasive mapping of micron-sized axon radii using commonly applied neuroimaging techniques, that is, diffusion-weighted MRI, has been bolstered by recent hardware developments, specifically MR gradient design. Here the whole brain characterization of the effective MR axon radius is presented and the inter- and intra-scanner test-retest repeatability and reproducibility are evaluated to promote the further development of the effective MR axon radius as a neuroimaging biomarker. A coefficient-of-variability of approximately 10% in the voxelwise estimation of the effective MR radius is observed in the test-retest analysis, but it is shown that the performance can be improved fourfold using a customized along-tract analysis.

RUN-UP: Accelerated multishot diffusion-weighted MRI reconstruction using an unrolled network with U-Net as priors

PURPOSE

To accelerate and improve multishot diffusion-weighted MRI reconstruction using deep learning.

METHODS

An unrolled pipeline containing recurrences of model-based gradient updates and neural networks was introduced for accelerating multishot DWI reconstruction with shot-to-shot phase correction. The network was trained to predict results of jointly reconstructed multidirection data using single-direction data as input. In vivo brain and breast experiments were performed for evaluation.

RESULTS

The proposed method achieves a reconstruction time of 0.1 second per image, over 100-fold faster than a shot locally low-rank reconstruction. The resultant image quality is comparable to the target from the joint reconstruction with a peak signal-to-noise ratio of 35.3 dB, a normalized root-mean-square error of 0.0177, and a structural similarity index of 0.944. The proposed method also improves upon the locally low-rank reconstruction (2.9 dB higher peak signal-to-noise ratio, 29% lower normalized root-mean-square error, and 0.037 higher structural similarity index). With training data from the brain, this method also generalizes well to breast diffusion-weighted imaging, and fine-tuning further reduces aliasing artifacts.

CONCLUSION

A proposed data-driven approach enables almost real-time reconstruction with improved image quality, which improves the feasibility of multishot DWI in a wide range of clinical and neuroscientific studies.

A simple estimate of axon size with diffusion MRI

Noninvasive estimation of mean axon diameter presents a new opportunity to explore white matter plasticity, development, and pathology. Several diffusion-weighted MRI (DW-MRI) methods have been proposed to measure the average axon diameter in white matter, but they typically require many diffusion encoding measurements and complicated mathematical models to fit the signal to multiple tissue compartments, including intra- and extra-axonal spaces. Here, Monte Carlo simulations uncovered a straightforward DW-MRI metric of axon diameter: the change in radial apparent diffusion coefficient estimated at different effective diffusion times, ΔD_⊥. Simulations indicated that this metric increases monotonically within a relevant range of effective mean axon diameter while being insensitive to changes in extra-axonal volume fraction, axon diameter distribution, g-ratio, and influence of myelin water. Also, a monotonic relationship was found to exist for signals coming from both intra- and extra-axonal compartments. The slope in ΔD_⊥ with effective axon diameter increased with the difference in diffusion time of both oscillating and pulsed gradient diffusion sequences.

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.

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.

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/mm²). 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.

Advanced Multicompartment Diffusion MRI Models and Their Application in Multiple Sclerosis

Conventional MR imaging techniques are sensitive to pathologic changes of the brain and spinal cord seen in MS, but they lack specificity for underlying axonal and myelin integrity. By isolating the signal contribution from different tissue compartments, newly developed advanced multicompartment diffusion MR imaging models have the potential to detect specific tissue subtypes and associated injuries with increased pathologic specificity. These models include neurite orientation dispersion and density imaging, diffusion basis spectrum imaging, multicompartment microscopic diffusion MR imaging with the spherical mean technique, and models enabled through high-gradient diffusion MR imaging. In this review, we provide an appraisal of the current literature on the physics principles, histopathologic validation, and clinical applications of each of these techniques in both brains and spinal cords of patients with MS. We discuss limitations of each of the methods and directions that future research could take to provide additional validation of their roles as biomarkers of axonal and myelin injury in MS.