The dot-compartment revealed? Diffusion MRI with ultra-strong gradients and spherical tensor encoding in the living human brain

A multimodal submillimeter MRI atlas of the human cerebellum

The human cerebellum is engaged in a broad array of tasks related to motor coordination, cognition, language, attention, memory, and emotional regulation. A detailed cerebellar atlas can facilitate the investigation of the structural and functional organization of the cerebellum. However, existing cerebellar atlases are typically limited to a single imaging modality with insufficient characterization of tissue properties. Here, we introduce a multifaceted cerebellar atlas based on high-resolution multimodal MRI, facilitating the understanding of the neurodevelopment and neurodegeneration of the cerebellum based on cortical morphology, tissue microstructure, and intra-cerebellar and cerebello-cerebral connectivity.

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.

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.

Implications of fitting a two-compartment model in single-shell diffusion MRI

It is becoming increasingly common for studies to fit single-shell diffusion MRI data to a two-compartment model, which comprises a hindered cellular compartment and a freely diffusing isotropic compartment. These studies consistently find that the fraction of the isotropic compartment (f) is sensitive to white matter (WM) conditions and pathologies, although the actual biological source of changes infhas not been validated. In this work we put aside the biological interpretation offand study the sensitivity implications of fitting single-shell data to a two-compartment model. We identify a nonlinear transformation between the one-compartment model (diffusion tensor imaging, DTI) and a two-compartment model in which the mean diffusivities of both compartments are effectively fixed. While the analytic relationship implies that fitting this two-compartment model does not offer any more information than DTI, it explains why metrics derived from a two-compartment model can exhibit enhanced sensitivity over DTI to certain types of WM processes, such as age-related WM differences. The sensitivity enhancement should not be viewed as a substitute for acquiring multi-shell data. Rather, the results of this study provide insight into the consequences of choosing a two-compartment model when only single-shell data is available.

Imaging local diffusion in microstructures using NV-based pulsed field gradient NMR

Understanding diffusion in microstructures plays a crucial role in many scientific fields, including neuroscience, medicine, or energy research. While magnetic resonance (MR) methods are the gold standard for diffusion measurements, spatial encoding in MR imaging has limitations. Here, we introduce nitrogen-vacancy (NV) center-based nuclear MR (NMR) spectroscopy as a powerful tool to probe diffusion within microscopic sample volumes. We have developed an experimental scheme that combines pulsed gradient spin echo (PGSE) with optically detected NV-NMR spectroscopy, allowing local quantification of molecular diffusion and flow. We demonstrate correlated optical imaging with spatially resolved PGSE NV-NMR experiments probing anisotropic water diffusion within an individual model microstructure. Our optically detected PGSE NV-NMR technique opens up prospects for extending the current capabilities of investigating diffusion processes with the future potential of probing single cells, tissue microstructures, or ion mobility in thin film materials for battery applications.

A novel framework for in-vivo diffusion tensor distribution MRI of the human brain

Neural tissue microstructure plays an important role in developmental, physiological and pathophysiological processes. Diffusion tensor distribution (DTD) MRI helps probe subvoxel heterogeneity by describing water diffusion within a voxel using an ensemble of non-exchanging compartments characterized by a probability density function of diffusion tensors. In this study, we provide a new framework for acquiring multiple diffusion encoding (MDE) images and estimating DTD from them in the human brain in vivo. We interfused pulsed field gradients (iPFG) in a single spin echo to generate arbitrary b-tensors of rank one, two, or three without introducing concomitant gradient artifacts. Employing well-defined diffusion encoding parameters we show that iPFG retains salient features of a traditional multiple-PFG (mPFG/MDE) sequence while reducing the echo time and coherence pathway artifacts thereby extending its applications beyond DTD MRI. Our DTD is a maximum entropy tensor-variate normal distribution whose tensor random variables are constrained to be positive definite to ensure their physicality. In each voxel, the second-order mean and fourth-order covariance tensors of the DTD are estimated using a Monte Carlo method that synthesizes micro-diffusion tensors with corresponding size, shape, and orientation distributions to best fit the measured MDE images. From these tensors we obtain the spectrum of diffusion tensor ellipsoid sizes and shapes, and the microscopic orientation distribution function (μODF) and microscopic fractional anisotropy (μFA) that disentangle the underlying heterogeneity within a voxel. Using the DTD-derived μODF, we introduce a new method to perform fiber tractography capable of resolving complex fiber configurations. The results revealed microscopic anisotropy in various gray and white matter regions and skewed MD distributions in cerebellar gray matter not observed previously. DTD MRI tractography captured complex white matter fiber organization consistent with known anatomy. DTD MRI also resolved some degeneracies associated with diffusion tensor imaging (DTI) and elucidated the source of diffusion heterogeneity which may help improve the diagnosis of various neurological diseases and disorders.

Ultra-strong diffusion-weighted MRI reveals cerebellar grey matter abnormalities in movement disorders

Structural brain MRI has proven invaluable in understanding movement disorder pathophysiology. However, most work has focused on grey/white matter volumetric (macrostructural) and white matter microstructural effects, limiting understanding of frequently implicated grey matter microstructural differences. Using ultra-strong spherical tensor encoding diffusion-weighted MRI, a persistent MRI signal was seen in healthy cerebellar grey matter even at high diffusion-weightings (b ​≥ 10,000 s/mm2). Quantifying the proportion of this signal (denoted fs), previously ascertained to originate from inside small spherical spaces, provides a potential proxy for cell body density. In this work, this approach was applied for the first time to a clinical cohort, including patients with diagnosed movement disorders in which the cerebellum has been implicated in symptom pathophysiology. Five control participants (control group 1, median age 24.5 years (20-39 years), imaged at two timepoints, demonstrated consistency in measurement of all three measures - MD (Mean Diffusivity) fs, and Ds (dot diffusivity)- with intraclass correlation coefficients (ICC) of 0.98, 0.86 and 0.76, respectively. Comparison with an older control group (control group 2 (n = 5), median age 51 years (43-58 years)) found no significant differences, neither with morphometric nor microstructural (MD (p = 0.36), fs (p = 0.17) and Ds (p = 0.22)) measures. The movement disorder cohort (Parkinson's Disease, n = 5, dystonia, n = 5. Spinocerebellar Ataxia 6, n = 5) when compared to the age-matched control cohort (Control Group 2) identified significantly lower MD (p < 0.0001 and p < 0.0001) and higher fs values (p < 0.0001 and p < 0.0001) in SCA6 and dystonia cohorts respectively. Lobar division of the cerebellum found these same differences in the superior and inferior posterior lobes, while no differences were seen in either the anterior lobes or with Ds measurements. In contrast to more conventional measures from diffusion tensor imaging, this framework provides enhanced specificity to differences in restricted spherical spaces in grey matter (including small cells) by eliminating signals from cerebrospinal fluid and axons. In the context of human and animal histopathology studies, these findings potentially implicate the cerebellar Purkinje and granule cells as contributors to the observed signal differences, with both cell types having been implicated in several neurological disorders through both postmortem and animal model studies. This novel microstructural imaging approach shows promise for improving movement disorder diagnosis, prognosis, and treatment.

Low-field magnetic resonance image enhancement via stochastic image quality transfer

Low-field (<1T) magnetic resonance imaging (MRI) scanners remain in widespread use in low- and middle-income countries (LMICs) and are commonly used for some applications in higher income countries e.g. for small child patients with obesity, claustrophobia, implants, or tattoos. However, low-field MR images commonly have lower resolution and poorer contrast than images from high field (1.5T, 3T, and above). Here, we present Image Quality Transfer (IQT) to enhance low-field structural MRI by estimating from a low-field image the image we would have obtained from the same subject at high field. Our approach uses (i) a stochastic low-field image simulator as the forward model to capture uncertainty and variation in the contrast of low-field images corresponding to a particular high-field image, and (ii) an anisotropic U-Net variant specifically designed for the IQT inverse problem. We evaluate the proposed algorithm both in simulation and using multi-contrast (T1-weighted, T2-weighted, and fluid attenuated inversion recovery (FLAIR)) clinical low-field MRI data from an LMIC hospital. We show the efficacy of IQT in improving contrast and resolution of low-field MR images. We demonstrate that IQT-enhanced images have potential for enhancing visualisation of anatomical structures and pathological lesions of clinical relevance from the perspective of radiologists. IQT is proved to have capability of boosting the diagnostic value of low-field MRI, especially in low-resource settings.

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.

Axial and radial axonal diffusivities and radii from single encoding strongly diffusion-weighted MRI

We enable the estimation of the per-axon axial diffusivity from single encoding, strongly diffusion-weighted, pulsed gradient spin echo data. Additionally, we improve the estimation of the per-axon radial diffusivity compared to estimates based on spherical averaging. The use of strong diffusion weightings in magnetic resonance imaging (MRI) allows to approximate the signal in white matter as the sum of the contributions from only axons. At the same time, spherical averaging leads to a major simplification of the modeling by removing the need to explicitly account for the unknown distribution of axonal orientations. However, the spherically averaged signal acquired at strong diffusion weightings is not sensitive to the axial diffusivity, which cannot therefore be estimated although needed for modeling axons - especially in the context of multi-compartmental modeling. We introduce a new general method for the estimation of both the axial and radial axonal diffusivities at strong diffusion weightings based on kernel zonal modeling. The method could lead to estimates that are free from partial volume bias with gray matter or other isotropic compartments. The method is tested on publicly available data from the MGH Adult Diffusion Human Connectome project. We report reference values of axonal diffusivities based on 34 subjects, and derive estimates of axonal radii from only two shells. The estimation problem is also addressed from the angle of the required data preprocessing, the presence of biases related to modeling assumptions, current limitations, and future possibilities.

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.

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.

Physiological effects of human body imaging with 300 mT/m gradients

PURPOSE

The use of high-performance gradient systems (i.e., high gradient strength and/or high slew rate) for human MRI is limited by physiological effects (including the elicitation of magnetophosphenes and peripheral nerve stimulation (PNS)). These effects, in turn, depend on the interaction between time-varying magnetic fields and the body, and thus on the participant's position with respect to the scanner's isocenter. This study investigated the occurrence of magnetophosphenes and PNS when scanning participants on a high-gradient (300 mT/m) system, for different gradient amplitudes, ramp times, and participant positions.

METHODS

Using a whole-body 300 mT/m gradient MRI system, a cohort of participants was scanned with the head, heart, and prostate at magnet isocenter and a train of trapezoidal bipolar gradient pulses, with ramp times from 0.88 to 4.20 ms and gradient amplitudes from 60 to 300 mT/m. Reports of magnetophosphenes and incidental reports of PNS were obtained. A questionnaire was used to record any additional subjective effects.

RESULTS

Magnetophosphenes were strongly dependent on participant position in the scanner. 87% of participants reported the effect with the heart at isocenter, 33% with the head at isocenter, and only 7% with the prostate at isocenter. PNS was most widely reported by participants for the vertical gradient axis (67% of participants), and was the dominant physiological effect for ramp times below 2 ms.

CONCLUSION

This study evaluates the probability of eliciting magnetophosphenes during whole-body imaging using an ultra-strong gradient MRI system. It provides empirical guidance on the use of high-performance gradient systems for whole-body human MRI.

Separating Glioma Hyperintensities From White Matter by Diffusion-Weighted Imaging With Spherical Tensor Encoding

Background

Tumor-related hyperintensities in high b-value diffusion-weighted imaging (DWI) are radiologically important in the workup of gliomas. However, the white matter may also appear as hyperintense, which may conflate interpretation.

Purpose

To investigate whether DWI with spherical b-tensor encoding (STE) can be used to suppress white matter and enhance the conspicuity of glioma hyperintensities unrelated to white matter.

Materials and Methods

Twenty-five patients with a glioma tumor and at least one pathology-related hyperintensity on DWI underwent conventional MRI at 3 T. The DWI was performed both with linear and spherical tensor encoding (LTE-DWI and STE-DWI). The LTE-DWI here refers to the DWI obtained with conventional diffusion encoding and averaged across diffusion-encoding directions. Retrospectively, the differences in contrast between LTE-DWI and STE-DWI, obtained at a b-value of 2,000 s/mm², were evaluated by comparing hyperintensities and contralateral normal-appearing white matter (NAWM) both visually and quantitatively in terms of the signal intensity ratio (SIR) and contrast-to-noise ratio efficiency (CNR_eff).

Results

The spherical tensor encoding DWI was more effective than LTE-DWI at suppressing signals from white matter and improved conspicuity of pathology-related hyperintensities. The median SIR improved in all cases and on average by 28%. The median (interquartile range) SIR was 1.9 (1.6 – 2.1) for STE and 1.4 (1.3 – 1.7) for LTE, with a significant difference of 0.4 (0.3 –0.5) (p < 10^–4, paired U-test). In 40% of the patients, the SIR was above 2 for STE-DWI, but with LTE-DWI, the SIR was below 2 for all patients. The CNR_eff of STE-DWI was significantly higher than of LTE-DWI: 2.5 (2 – 3.5) vs. 2.3 (1.7 – 3.1), with a significant difference of 0.4 (−0.1 –0.6) (p < 10^–3, paired U-test). The STE improved CNR_eff in 70% of the cases. We illustrate the benefits of STE-DWI in three patients, where STE-DWI may facilitate an improved radiological description of tumor-related hyperintensity, including one case that could have been missed out if only LTE-DWI was inspected.

Conclusion

The contrast mechanism of high b-value STE-DWI results in a stronger suppression of white matter than conventional LTE-DWI, and may, therefore, be more sensitive and specific for assessment of glioma tumors and DWI-hyperintensities.

In vivo Correlation Tensor MRI reveals microscopic kurtosis in the human brain on a clinical 3T scanner

Diffusion MRI (dMRI) has become one of the most important imaging modalities for noninvasively probing tissue microstructure. Diffusional Kurtosis MRI (DKI) quantifies the degree of non-gaussian diffusion, which in turn has been shown to increase sensitivity towards, e.g., disease and orientation mapping in neural tissue. However, the specificity of DKI is limited as different sources can contribute to the total intravoxel diffusional kurtosis, including: variance in diffusion tensor magnitudes (K_iso), variance due to diffusion anisotropy (K_aniso), and microscopic kurtosis (μK) related to restricted diffusion, microstructural disorder, and/or exchange. Interestingly, μK is typically ignored in diffusion MRI signal modeling as it is assumed to be negligible in neural tissues. However, recently, Correlation Tensor MRI (CTI) based on Double-Diffusion-Encoding (DDE) was introduced for kurtosis source separation, revealing non negligible μK in preclinical imaging. Here, we implemented CTI for the first time on a clinical 3T scanner and investigated the sources of total kurtosis in healthy subjects. A robust framework for kurtosis source separation in humans is introduced, followed by estimation of μK (and the other kurtosis sources) in the healthy brain. Using this clinical CTI approach, we find that μK significantly contributes to total diffusional kurtosis both in gray and white matter tissue but, as expected, not in the ventricles. The first μK maps of the human brain are presented, revealing that the spatial distribution of μK provides a unique source of contrast, appearing different from isotropic and anisotropic kurtosis counterparts. Moreover, group average templates of these kurtosis sources have been generated for the first time, which corroborated our findings at the underlying individual-level maps. We further show that the common practice of ignoring μK and assuming the multiple gaussian component approximation for kurtosis source estimation introduces significant bias in the estimation of other kurtosis sources and, perhaps even worse, compromises their interpretation. Finally, a twofold acceleration of CTI is discussed in the context of potential future clinical applications. We conclude that CTI has much potential for future in vivo microstructural characterizations in healthy and pathological tissue.

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 dMRI 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 dMRI and where the field is headed in the coming years.

Diffusion time dependence, power-law scaling, and exchange in gray matter

Characterizing neural tissue microstructure is a critical goal for future neuroimaging. Diffusion MRI (dMRI) provides contrasts that reflect diffusing spins’ interactions with myriad microstructural features of biological systems. However, the specificity of dMRI remains limited due to the ambiguity of its signals vis-à-vis the underlying microstructure. To improve specificity, biophysical models of white matter (WM) typically express dMRI signals according to the Standard Model (SM) and have more recently in gray matter (GM) taken spherical compartments into account (the SANDI model) in attempts to represent cell soma. The validity of the assumptions underlying these models, however, remains largely undetermined, especially in GM. To validate these assumptions experimentally, observing their unique, functional properties, such as the b^−1/2 power-law associated with one-dimensional diffusion, has emerged as a fruitful strategy. The absence of this signature in GM, in turn, has been explained by neurite water exchange, non-linear morphology, and/or by obscuring soma signal contributions. Here, we present diffusion simulations in realistic neurons demonstrating that curvature and branching does not destroy the stick power-law behavior in impermeable neurites, but also that their signal is drowned by the soma signal under typical experimental conditions. Nevertheless, by studying the GM dMRI signal's behavior as a function of diffusion weighting as well as time, we identify an attainable experimental regime in which the neurite signal dominates. Furthermore, we find that exchange-driven time dependence produces a signal behavior opposite to that which would be expected from restricted diffusion, thereby providing a functional signature that disambiguates the two effects. We present data from dMRI experiments in ex vivo rat brain at ultrahigh field of 16.4T and observe a time dependence that is consistent with substantial exchange but also with a GM stick power-law. The first finding suggests significant water exchange between neurites and the extracellular space while the second suggests a small sub-population of impermeable neurites. To quantify these observations, we harness the Kärger exchange model and incorporate the corresponding signal time dependence in the SM and SANDI models.

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.

Disentangling the effects of restriction and exchange with diffusion exchange spectroscopy

Diffusion exchange spectroscopy (DEXSY) is a multidimensional NMR technique that can reveal how water molecules exchange between compartments within heterogeneous media, such as biological tissue. Data from DEXSY experiments is typically processed using numerical inverse Laplace transforms (ILTs) to produce a diffusion-diffusion spectrum. A tacit assumption of this ILT approach is that the signal behavior is Gaussian - i.e., the spin echo intensity decays exponentially with the degree of diffusion weighting. The assumptions that underlie Gaussian signal behavior may be violated, however, depending on the gradient strength applied and the sample under study. We argue that non-Gaussian signal behavior due to restrictions is to be expected in the study of biological tissue using diffusion NMR. Further, we argue that this signal behavior can produce confounding features in the diffusion-diffusion spectra obtained from numerical ILTs of DEXSY data - entangling the effects of restriction and exchange. Specifically, restricted signal behavior can result in broadening of peaks and in the appearance of illusory exchanging compartments with distributed diffusivities, which pearl into multiple peaks if not highly regularized. We demonstrate these effects on simulated data. That said, we suggest the use of features in the signal acquisition domain that can be used to rapidly probe exchange without employing an ILT. We also propose a means to characterize the non-Gaussian signal behavior due to restrictions within a sample using DEXSY measurements with a near zero mixing time or storage interval. We propose a combined acquisition scheme to independently characterize restriction and exchange with various DEXSY measurements, which we term Restriction and Exchange from Equally-weighted Double and Single Diffusion Encodings (REEDS-DE). We test this method on ex vivo neonatal mouse spinal cord - a sample consisting primarily of gray matter - using a low-field, static gradient NMR system. In sum, we highlight critical shortcomings of prevailing DEXSY analysis methods that conflate the effects of restriction and exchange, and suggest a viable experimental approach to disentangle them.

Multi-tissue spherical deconvolution of tensor-valued diffusion MRI

Multi-tissue constrained spherical deconvolution (MT-CSD) leverages the characteristic b-value dependency of each tissue type to estimate both the apparent tissue densities and the white matter fiber orientation distribution function from diffusion MRI data. In this work, we generalize MT-CSD to tensor-valued diffusion encoding with arbitrary b-tensor shapes. This enables the use of data encoded with mixed b-tensors, rather than being limited to the subset of linear (conventional) b-tensors. Using the complete set of data, including all b-tensor shapes, provides a categorical improvement in the estimation of apparent tissue densities, fiber ODF, and resulting tractography. Furthermore, we demonstrate that including multiple b-tensor shapes in the analysis provides improved contrast between tissue types, in particular between gray matter and white matter. We also show that our approach provides high-quality apparent tissue density maps and high-quality fiber tracking from data, even with sparse sampling across b-tensors that yield whole-brain coverage at 2 mm isotropic resolution in approximately 5:15 min.

Axonal T2 estimation using the spherical variance of the strongly diffusion-weighted MRI signal

In magnetic resonance imaging, the application of a strong diffusion weighting suppresses the signal contributions from the less diffusion-restricted constituents of the brain's white matter, thus enabling the estimation of the transverse relaxation time T₂ that arises from the more diffusion-restricted constituents such as the axons. However, the presence of cell nuclei and vacuoles can confound the estimation of the axonal T₂, as diffusion within those structures is also restricted, causing the corresponding signal to survive the strong diffusion weighting. We devise an estimator of the axonal T₂ based on the directional spherical variance of the strongly diffusion-weighted signal. The spherical variance T₂ estimates are insensitive to the presence of isotropic contributions to the signal like those provided by cell nuclei and vacuoles. We show that with a strong diffusion weighting these estimates differ from those obtained using the directional spherical mean of the signal which contains both axonal and isotropically-restricted contributions. Our findings hint at the presence of an MRI-visible isotropically-restricted contribution to the signal in the white matter ex vivo fixed tissue (monkey) at 7T, and do not allow us to discard such a possibility also for in vivo human data collected with a clinical 3T system.

A Pilot Study of Multidimensional Diffusion MRI for Assessment of Tissue Heterogeneity in Prostate Cancer

OBJECTIVES

The objectives of this exploratory study were to investigate the feasibility of multidimensional diffusion magnetic resonance imaging (MddMRI) in assessing diffusion heterogeneity at both a macroscopic and microscopic level in prostate cancer (PCa).

MATERIALS AND METHODS

Informed consent was obtained from 46 subjects who underwent 3.0-T prostate multiparametric MRI, complemented with a prototype spin echo-based MddMRI sequence in this institutional review board-approved study. Prostate cancer tumors and comparative normal tissue from each patient were contoured on both apparent diffusion coefficient and MddMRI-derived mean diffusivity (MD) maps (from which microscopic diffusion heterogeneity [MKi] and microscopic diffusion anisotropy were derived) using 3D Slicer. The discriminative ability of MddMRI-derived parameters to differentiate PCa from normal tissue was determined using the Friedman test. To determine if tumor diffusion heterogeneity is similar on macroscopic and microscopic scales, the linear association between SD of MD and mean MKi was estimated using robust regression (bisquare weighting). Hypothesis testing was 2 tailed; P values less than 0.05 were considered statistically significant.

RESULTS

All MddMRI-derived parameters could distinguish tumor from normal tissue in the fixed-effects analysis (P < 0.0001). Tumor MKi was higher (P < 0.05) compared with normal tissue (median, 0.40; interquartile range, 0.29-0.52 vs 0.20-0.18; 0.25), as was tumor microscopic diffusion anisotropy (0.55; 0.36-0.81 vs 0.20-0.15; 0.28). The MKi could not be predicted (no significant association) by SD of MD. There was a significant correlation between tumor volume and SD of MD (R2 = 0.50, slope = 0.008 μm2/ms per millimeter, P < 0.001) but not between tumor volume and MKi.

CONCLUSIONS

This explorative study demonstrates that MddMRI provides novel information on MKi and microscopic anisotropy, which differ from measures at the macroscopic level. MddMRI has the potential to characterize tumor tissue heterogeneity at different spatial scales.

Fiber tractography bundle segmentation depends on scanner effects, vendor effects, acquisition resolution, diffusion sampling scheme, diffusion sensitization, and bundle segmentation workflow

When investigating connectivity and microstructure of white matter pathways of the brain using diffusion tractography bundle segmentation, it is important to understand potential confounds and sources of variation in the process. While cross-scanner and cross-protocol effects on diffusion microstructure measures are well described (in particular fractional anisotropy and mean diffusivity), it is unknown how potential sources of variation effect bundle segmentation results, which features of the bundle are most affected, where variability occurs, nor how these sources of variation depend upon the method used to reconstruct and segment bundles. In this study, we investigate six potential sources of variation, or confounds, for bundle segmentation: variation (1) across scan repeats, (2) across scanners, (3) across vendors (4) across acquisition resolution, (5) across diffusion schemes, and (6) across diffusion sensitization. We employ four different bundle segmentation workflows on two benchmark multi-subject cross-scanner and cross-protocol databases, and investigate reproducibility and biases in volume overlap, shape geometry features of fiber pathways, and microstructure features within the pathways. We find that the effects of acquisition protocol, in particular acquisition resolution, result in the lowest reproducibility of tractography and largest variation of features, followed by vendor-effects, scanner-effects, and finally diffusion scheme and b-value effects which had similar reproducibility as scan-rescan variation. However, confounds varied both across pathways and across segmentation workflows, with some bundle segmentation workflows more (or less) robust to sources of variation. Despite variability, bundle dissection is consistently able to recover the same location of pathways in the deep white matter, with variation at the gray matter/ white matter interface. Next, we show that differences due to the choice of bundle segmentation workflows are larger than any other studied confound, with low-to-moderate overlap of the same intended pathway when segmented using different methods. Finally, quantifying microstructure features within a pathway, we show that tractography adds variability over-and-above that which exists due to noise, scanner effects, and acquisition effects. Overall, these confounds need to be considered when harmonizing diffusion datasets, interpreting or combining data across sites, and when attempting to understand the successes and limitations of different methodologies in the design and development of new tractography or bundle segmentation methods.

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.

Accuracy and precision in super-resolution MRI: Enabling spherical tensor diffusion encoding at ultra-high b-values and high resolution

Diffusion MRI (dMRI) can probe the tissue microstructure but suffers from low signal-to-noise ratio (SNR) whenever high resolution is combined with high diffusion encoding strengths. Low SNR leads to poor precision as well as poor accuracy of the diffusion-weighted signal; the latter is caused by the rectified noise floor and can be observed as a positive bias in magnitude signal. Super-resolution techniques may facilitate a beneficial tradeoff between bias and resolution by allowing acquisition at low spatial resolution and high SNR, whereafter high spatial resolution is recovered by image reconstruction. In this work, we describe a super-resolution reconstruction framework for dMRI and investigate its performance with respect to signal accuracy and precision. Using phantom experiments and numerical simulations, we show that the super-resolution approach improves accuracy by facilitating a more beneficial trade-off between spatial resolution and diffusion encoding strength before the noise floor affects the signal. By contrast, precision is shown to have a less straightforward dependency on acquisition, reconstruction, and intrinsic tissue parameters. Indeed, we find a gain in precision from super-resolution reconstruction is substantial only when some spatial resolution is sacrificed. Finally, we deployed super-resolution reconstruction in a healthy brain for the challenging combination of spherical b-tensor encoding at ultra-high b-values and high spatial resolution—a configuration that produces a unique contrast that emphasizes tissue in which diffusion is restricted in all directions. This demonstration showcased that super-resolution reconstruction enables a vastly superior image contrast compared to conventional imaging, facilitating investigations that would otherwise have prohibitively low SNR, resolution or require non-conventional MRI hardware.

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.

SPHERIOUSLY? The challenges of estimating sphere radius non-invasively in the human brain from diffusion MRI

The Soma and Neurite Density Imaging (SANDI) three-compartment model was recently proposed to disentangle cylindrical and spherical geometries, attributed to neurite and soma compartments, respectively, in brain tissue. There are some recent advances in diffusion-weighted MRI signal encoding and analysis (including the use of multiple so-called 'b-tensor' encodings and analysing the signal in the frequency-domain) that have not yet been applied in the context of SANDI. In this work, using: (i) ultra-strong gradients; (ii) a combination of linear, planar, and spherical b-tensor encodings; and (iii) analysing the signal in the frequency domain, three main challenges to robust estimation of sphere size were identified: First, the Rician noise floor in magnitude-reconstructed data biases estimates of sphere properties in a non-uniform fashion. It may cause overestimation or underestimation of the spherical compartment size and density. This can be partly ameliorated by accounting for the noise floor in the estimation routine. Second, even when using the strongest diffusion-encoding gradient strengths available for human MRI, there is an empirical lower bound on the spherical signal fraction and radius that can be detected and estimated robustly. For the experimental setup used here, the lower bound on the sphere signal fraction was approximately 10%. We employed two different ways of establishing the lower bound for spherical radius estimates in white matter. The first, examining power-law relationships between the DW-signal and diffusion weighting in empirical data, yielded a lower bound of 7μm, while the second, pure Monte Carlo simulations, yielded a lower limit of 3μm and in this low radii domain, there is little differentiation in signal attenuation. Third, if there is sensitivity to the transverse intra-cellular diffusivity in cylindrical structures, e.g., axons and cellular projections, then trying to disentangle two diffusion-time-dependencies using one experimental parameter (i.e., change in frequency-content of the encoding waveform) makes spherical radii estimates particularly challenging. We conclude that due to the aforementioned challenges spherical radii estimates may be biased when the corresponding sphere signal fraction is low, which must be considered.

Measuring compartmental T2-orientational dependence in human brain white matter using a tiltable RF coil and diffusion-T2 correlation MRI

The anisotropy of brain white matter microstructure manifests itself in orientational-dependence of various MRI contrasts, and can result in significant quantification biases if ignored. Understanding the origins of this orientation-dependence could enhance the interpretation of MRI signal changes in development, ageing and disease and ultimately improve clinical diagnosis. Using a novel experimental setup, this work studies the contributions of the intra- and extra-axonal water to the orientation-dependence of one of the most clinically-studied parameters, apparent transverse relaxation T2. Specifically, a tiltable receive coil is interfaced with an ultra-strong gradient MRI scanner to acquire multidimensional MRI data with an unprecedented range of acquisition parameters. Using this setup, compartmental T2 can be disentangled based on differences in diffusional-anisotropy, and its orientation-dependence further elucidated by re-orienting the head with respect to the main magnetic field B→0. A dependence of (compartmental) T2 on the fibre orientation w.r.t. B→0 was observed, and further quantified using characteristic representations for susceptibility- and magic angle effects. Across white matter, anisotropy effects were dominated by the extra-axonal water signal, while the intra-axonal water signal decay varied less with fibre-orientation. Moreover, the results suggest that the stronger extra-axonal T2 orientation-dependence is dominated by magnetic susceptibility effects (presumably from the myelin sheath) while the weaker intra-axonal T2 orientation-dependence may be driven by a combination of microstructural effects. Even though the current design of the tiltable coil only offers a modest range of angles, the results demonstrate an overall effect of tilt and serve as a proof-of-concept motivating further hardware development to facilitate experiments that explore orientational anisotropy. These observations have the potential to lead to white matter microstructural models with increased compartmental sensitivity to disease, and can have direct consequences for longitudinal and group-wise T2- and diffusion-MRI data analysis, where the effect of head-orientation in the scanner is commonly ignored.

Accelerated diffusion and relaxation-diffusion MRI using time-division multiplexing EPI

PURPOSE

To develop a time-division multiplexing echo-planar imaging (TDM-EPI) sequence for approximately two- to threefold acceleration when acquiring joint relaxation-diffusion MRI data with multiple TEs.

METHODS

The proposed TDM-EPI sequence interleaves excitation and data collection for up to 3 separate slices at different TEs and uses echo-shifting gradients to disentangle the overlapping echo signals during the readout period. By properly arranging the sequence event blocks for each slice and adjusting the echo-shifting gradients, diffusion-weighted images from separate slices can be acquired. Therefore, we present 2 variants of the sequence. A single-TE TDM-EPI is presented to demonstrate the concept. Next, a multi-TE TDM-EPI is presented to highlight the advantages of the TDM approach for relaxation-diffusion imaging. These sequences were evaluated on a 3 Tesla scanner with a water phantom and in vivo human brain data.

RESULTS

The single-TE TDM-EPI sequence can simultaneously acquire 2 slices with a maximum b value of 3000 s/mm² and 2.5 mm isotropic resolution using interleaved readout windows with TE ≈ 78 ms. With the same b value and resolution, the multi-TE TDM-EPI sequence can simultaneously acquire 2 or 3 separate slices using interleaved readout sections with shortest TE ≈ 70 ms and ΔTE ≈ 30 ms. Phantom and in vivo experiments have shown that the proposed TDM-EPI sequences can provide similar image quality and diffusion measures as conventional EPI readouts with multiple echoes but can reduce the overall relaxation-diffusion protocol scan time by approximately two- to threefold.

CONCLUSION

TDM-EPI is a novel approach to acquire diffusion imaging data at multiple TEs. This enables a significant reduction in acquisition time for relaxation-diffusion MRI experiments but without compromising image quality and diffusion measurements, thus removing a significant barrier to the adoption of relaxation-diffusion MRI in clinical research studies of neurological and mental disorders.

On the potential for mapping apparent neural soma density via a clinically viable diffusion MRI protocol

•

B-tensor encoding enables estimation of spherical cellular structures in the brain.

•

Spherical compartments may provide markers for apparent neural soma density.

•

Model parameters can be estimated in a fast and robust way using deep learning.

•

Practical acquisition times are achievable on widely available clinical scanners.

Diffusion MRI is a valuable tool for probing tissue microstructure in the brain noninvasively. Today, model-based techniques are widely available and used for white matter characterisation where their development is relatively mature. Conversely, tissue modelling in grey matter is more challenging, and no generally accepted models exist. With advances in measurement technology and modelling efforts, a clinically viable technique that reveals salient features of grey matter microstructure, such as the density of quasi-spherical cell bodies and quasi-cylindrical cell projections, is an exciting prospect. As a step towards capturing the microscopic architecture of grey matter in clinically feasible settings, this work uses a biophysical model that is designed to disentangle the diffusion signatures of spherical and cylindrical structures in the presence of orientation heterogeneity, and takes advantage of B-tensor encoding measurements, which provide additional sensitivity compared to standard single diffusion encoding sequences. For the fast and robust estimation of microstructural parameters, we leverage recent advances in machine learning and replace conventional fitting techniques with an artificial neural network that fits complex biophysical models within seconds. Our results demonstrate apparent markers of spherical and cylindrical geometries in healthy human subjects, and in particular an increased volume fraction of spherical compartments in grey matter compared to white matter. We evaluate the extent to which spherical and cylindrical geometries may be interpreted as correlates of neural soma and neural projections, respectively, and quantify parameter estimation errors in the presence of various departures from the modelling assumptions. While further work is necessary to translate the ideas presented in this work to the clinic, we suggest that biomarkers focussing on quasi-spherical cellular geometries may be valuable for the enhanced assessment of neurodevelopmental disorders and neurodegenerative diseases.

Whole-Brain Imaging of Subvoxel T1-Diffusion Correlation Spectra in Human Subjects

T1 relaxation and water mobility generate eloquent MRI tissue contrasts with great diagnostic value in many neuroradiological applications. However, conventional methods do not adequately quantify the microscopic heterogeneity of these important biophysical properties within a voxel, and therefore have limited biological specificity. We describe a new correlation spectroscopic (CS) MRI method for measuring how T1 and mean diffusivity (MD) co-vary in microscopic tissue environments. We develop a clinical pulse sequence that combines inversion recovery (IR) with single-shot isotropic diffusion encoding (IDE) to efficiently acquire whole-brain MRIs with a wide range of joint T1-MD weightings. Unlike conventional diffusion encoding, the IDE preparation ensures that all subvoxel water pools are weighted by their MDs regardless of the sizes, shapes, and orientations of their corresponding microscopic diffusion tensors. Accordingly, IR-IDE measurements are well-suited for model-free, quantitative spectroscopic analysis of microscopic water pools. Using numerical simulations, phantom experiments, and data from healthy volunteers we demonstrate how IR-IDE MRIs can be processed to reconstruct maps of two-dimensional joint probability density functions, i.e., correlation spectra, of subvoxel T1-MD values. In vivo T1-MD spectra show distinct cerebrospinal fluid and parenchymal tissue components specific to white matter, cortical gray matter, basal ganglia, and myelinated fiber pathways, suggesting the potential for improved biological specificity. The one-dimensional marginal distributions derived from the T1-MD correlation spectra agree well with results from other relaxation spectroscopic and quantitative MRI studies, validating the T1-MD contrast encoding and the spectral reconstruction. Mapping subvoxel T1-diffusion correlations in patient populations may provide a more nuanced, comprehensive, sensitive, and specific neuroradiological assessment of the non-specific changes seen on fluid-attenuated inversion recovery (FLAIR) and diffusion-weighted MRIs (DWIs) in cancer, ischemic stroke, or brain injury.

Mapping prostatic microscopic anisotropy using linear and spherical b-tensor encoding: A preliminary study

Purpose:

Tensor-valued diffusion encoding provides more specific information than conventional diffusion-weighted imaging (DWI), but has mainly been applied in neuroimaging studies. This study aimed to assess its potential for the imaging of prostate cancer (PCa).

Methods:

Seventeen patients with histologically proven PCa were enrolled. DWI of the prostate was performed with linear and spherical tensor encoding using a maximal b-value of 1.5 ms/μm² and a voxel size of 3 × 3 × 4 mm³. The gamma-distribution model was used to estimate the mean diffusivity (MD), the isotropic kurtosis (MK_I), and the anisotropic kurtosis (MK_A). Regions of interest were placed in MR-defined cancerous tissues, as well as in apparently healthy tissues in the peripheral and transitional zones (PZs and TZs).

Results:

DWI with linear and spherical encoding yielded different image contrasts at high b-values, which enabled the estimation of MK_A and MK_I. Compared with healthy tissue (PZs and TZs combined) the cancers displayed a significantly lower MD (P < .05), higher MK_I (P < 10⁻⁵), and lower MK_A (P < .05). Compared with the TZ, tissue in the PZ showed lower MD (P < 10⁻³) and higher MK_A (P < 10⁻³). No significant differences were found between cancers of different Gleason scores, possibly because of the limited sample size.

Conclusion:

Tensor-valued diffusion encoding enabled mapping of MK_A and MK_I in the prostate. The elevated MK_I in PCa compared with normal tissues suggests an elevated heterogeneity in the cancers. Increased in-plane resolution could improve tumor delineation in future studies.

Glioma grading, molecular feature classification, and microstructural characterization using MR diffusional variance decomposition (DIVIDE) imaging

OBJECTIVE

To evaluate the potential of diffusional variance decomposition (DIVIDE) for grading, molecular feature classification, and microstructural characterization of gliomas.

MATERIALS AND METHODS

Participants with suspected gliomas underwent DIVIDE imaging, yielding parameter maps of fractional anisotropy (FA), mean diffusivity (MD), anisotropic mean kurtosis (MK_A), isotropic mean kurtosis (MK_I), total mean kurtosis (MK_T), MK_A/MK_T, and microscopic fractional anisotropy (μFA). Tumor type and grade, isocitrate dehydrogenase (IDH) 1/2 mutant status, and the Ki-67 labeling index (Ki-67 LI) were determined after surgery. Statistical analysis included 33 high-grade gliomas (HGG) and 17 low-grade gliomas (LGG). Tumor diffusion metrics were compared between HGG and LGG, among grades, and between wild and mutated IDH types using appropriate tests according to normality assessment results. Receiver operating characteristic and Spearman correlation analysis were also used for statistical evaluations.

RESULTS

FA, MD, MK_A, MK_I, MK_T, μFA, and MK_A/MK_T differed between HGG and LGG (FA: p = 0.047; MD: p = 0.037, others p < 0.001), and among glioma grade II, III, and IV (FA: p = 0.048; MD: p = 0.038, others p < 0.001). All diffusion metrics differed between wild-type and mutated IDH tumors (MK_I: p = 0.003; others: p < 0.001). The metrics that best discriminated between HGG and LGGs and between wild-type and mutated IDH tumors were MK_T and FA respectively (area under the curve 0.866 and 0.881). All diffusion metrics except FA showed significant correlation with Ki-67 LI, and MK_I had the highest correlation coefficient (r_s = 0.618).

CONCLUSION

DIVIDE is a promising technique for glioma characterization and diagnosis.

KEY POINTS

• DIVIDE metrics MK_I is related to cell density heterogeneity while MK_A and μFA are related to cell eccentricity. • DIVIDE metrics can effectively differentiate LGG from HGG and IDH mutation from wild-type tumor, and showed significant correlation with the Ki-67 labeling index. • MK_I was larger than MK_A which indicates predominant cell density heterogeneity in gliomas. • MK_A and MK_I increased with grade or degree of malignancy, however with a relatively larger increase in the cell eccentricity metric MK_A in relation to the cell density heterogeneity metric MK_I.

Compartmental diffusion and microstructural properties of human brain gray and white matter studied with double diffusion encoding magnetic resonance spectroscopy of metabolites and water

Double diffusion encoding (DDE) of the water signal offers a unique ability to separate the effect of microscopic anisotropic diffusion in structural units of tissue from the overall macroscopic orientational distribution of cells. However, the specificity in detected microscopic anisotropy is limited as the signal is averaged over different cell types and across tissue compartments. Performing side-by-side water and metabolite DDE spectroscopic (DDES) experiments provides complementary measures from which intracellular and extracellular microscopic fractional anisotropies (μFA) and diffusivities can be estimated. Metabolites are largely confined to the intracellular space and therefore provide a benchmark for intracellular μFA and diffusivities of specific cell types. By contrast, water DDES measurements allow examination of the separate contributions to water μFA and diffusivity from the intra- and extracellular spaces, by using a wide range of b values to gradually eliminate the extracellular contribution. Here, we aimed to estimate tissue and compartment specific human brain microstructure by combining water and metabolites DDES experiments. We performed our DDES measurements in two brain regions that contain widely different amounts of white matter (WM) and gray matter (GM): parietal white matter (PWM) and occipital gray matter (OGM) in a total of 20 healthy volunteers at 7 Tesla. Metabolite DDES measurements were performed at b = 7199 s/mm², while water DDES measurements were performed with a range of b values from 918 to 7199 s/mm². The experimental framework we employed here resulted in a set of insights pertaining to the morphology of the intracellular and extracellular spaces in both gray and white matter. Results of the metabolite DDES experiments in both PWM and OGM suggest a highly anisotropic intracellular space within neurons and glia, with the possible exception of gray matter glia. The water μFA obtained from the DDES results at high b values in both regions converged with that of the metabolite DDES, suggesting that the signal from the extracellular space is indeed effectively suppressed at the highest b value. The μFA measured in the OGM significantly decreased at lower b values, suggesting a considerably lower anisotropy of the extracellular space in GM compared to WM. In PWM, the water μFA remained high even at the lowest b value, indicating a high degree of organization in the interstitial space in WM. Tortuosity values in the cytoplasm for water and tNAA, obtained with correlation analysis of microscopic parallel diffusivity with respect to GM/WM tissue fraction in the volume of interest, are remarkably similar for both molecules, while exhibiting a clear difference between gray and white matter, suggesting a more crowded cytoplasm and more complex cytomorphology of neuronal cell bodies and dendrites in GM than those found in long-range axons in WM.

Beyond the diffusion standard model in fixed rat spinal cord with combined linear and planar encoding

Information about tissue on the microscopic and mesoscopic scales can be accessed by modelling diffusion MRI signals, with the aim of extracting microstructure-specific biomarkers. The standard model (SM) of diffusion, currently the most broadly adopted microstructural model, describes diffusion in white matter (WM) tissues by two Gaussian components, one of which has zero radial diffusivity, to represent diffusion in intra- and extra-axonal water, respectively. Here, we reappraise these SM assumptions by collecting comprehensive double diffusion encoded (DDE) MRI data with both linear and planar encodings, which was recently shown to substantially enhance the ability to estimate SM parameters. We find however, that the SM is unable to account for data recorded in fixed rat spinal cord at an ultrahigh field of 16.4 T, suggesting that its underlying assumptions are violated in our experimental data. We offer three model extensions to mitigate this problem: first, we generalize the SM to accommodate finite radii (axons) by releasing the constraint of zero radial diffusivity in the intra-axonal compartment. Second, we include intracompartmental kurtosis to account for non-Gaussian behaviour. Third, we introduce an additional (third) compartment. The ability of these models to account for our experimental data are compared based on parameter feasibility and Bayesian information criterion. Our analysis identifies the three-compartment description as the optimal model. The third compartment exhibits slow diffusion with a minor but non-negligible signal fraction (∼12%). We demonstrate how failure to take the presence of such a compartment into account severely misguides inferences about WM microstructure. Our findings bear significance for microstructural modelling at large and can impact the interpretation of biomarkers extracted from the standard model of diffusion.

Gradient waveform design for tensor-valued encoding in diffusion MRI

Diffusion encoding along multiple spatial directions per signal acquisition can be described in terms of a b-tensor. The benefit of tensor-valued diffusion encoding is that it unlocks the 'shape of the b-tensor' as a new encoding dimension. By modulating the b-tensor shape, we can control the sensitivity to microscopic diffusion anisotropy which can be used as a contrast mechanism; a feature that is inaccessible by conventional diffusion encoding. Since imaging methods based on tensor-valued diffusion encoding are finding an increasing number of applications we are prompted to highlight the challenge of designing the optimal gradient waveforms for any given application. In this review, we first establish the basic design objectives in creating field gradient waveforms for tensor-valued diffusion MRI. We also survey additional design considerations related to limitations imposed by hardware and physiology, potential confounding effects that cannot be captured by the b-tensor, and artifacts related to the diffusion encoding waveform. Throughout, we discuss the expected compromises and tradeoffs with an aim to establish a more complete understanding of gradient waveform design and its impact on accurate measurements and interpretations of data.

Toward Quantification: Microstructure and Magnetic Resonance Fingerprinting

Quantitative magnetic resonance imaging (MRI) is a long-standing challenge. We advocate that the origin of the problem is the simplification applied in commonly used models of the MRI signal relation to the target parameters of biological tissues. Two research fields are briefly reviewed as ways to respond to the challenge of quantitative MRI, both experiencing an exponential growth right now. Microstructure MRI strives to build physiology-based models from cells to signal and, given the signal, back to the cells again. Magnetic resonance fingerprinting aims at efficient simultaneous determination of multiple signal parameters. The synergy of these yet disjoined approaches promises truly quantitative MRI with specific target-oriented diagnostic tools rather than universal imaging methods.

Double diffusion encoding and applications for biomedical imaging

Diffusion Magnetic Resonance Imaging (dMRI) is one of the most important contemporary non-invasive modalities for probing tissue structure at the microscopic scale. The majority of dMRI techniques employ standard single diffusion encoding (SDE) measurements, covering different sequence parameter ranges depending on the complexity of the method. Although many signal representations and biophysical models have been proposed for SDE data, they are intrinsically limited by a lack of specificity. Advanced dMRI methods have been proposed to provide additional microstructural information beyond what can be inferred from SDE. These enhanced contrasts can play important roles in characterizing biological tissues, for instance upon diseases (e.g. neurodegenerative, cancer, stroke), aging, learning, and development. In this review we focus on double diffusion encoding (DDE), which stands out among other advanced acquisitions for its versatility, ability to probe more specific diffusion correlations, and feasibility for preclinical and clinical applications. Various DDE methodologies have been employed to probe compartment sizes (Section 3), decouple the effects of microscopic diffusion anisotropy from orientation dispersion (Section 4), probe displacement correlations, study exchange, or suppress fast diffusing compartments (Section 6). DDE measurements can also be used to improve the robustness of biophysical models (Section 5) and study intra-cellular diffusion via magnetic resonance spectroscopy of metabolites (Section 7). This review discusses all these topics as well as important practical aspects related to the implementation and contrast in preclinical and clinical settings (Section 9) and aims to provide the readers a guide for deciding on the right DDE acquisition for their specific application.

Characterizing Intra-soma Diffusion with Spherical Mean Spectrum Imaging

Most brain microstructure models are dedicated to the quantification of white matter microstructure, using for example sticks, cylinders, and zeppelins to model intra- and extra-axonal environments. Gray matter presents unique micro-architecture with cell bodies (somas) exhibiting diffusion characteristics that differ from axons in white matter. In this paper, we introduce a method to quantify soma microstructure, giving measures such as volume fraction, diffusivity, and kurtosis. Our method captures a spectrum of diffusion patterns and scales and does not rely on restrictive model assumptions. We show that our method yields unique and meaningful contrasts that are in agreement with histological data. We demonstrate its application in the mapping of the distinct spatial patterns of soma density in the cortex.

Challenges for biophysical modeling of microstructure

The biophysical modeling efforts in diffusion MRI have grown considerably over the past 25 years. In this review, we dwell on the various challenges along the journey of bringing a biophysical model from initial design to clinical implementation, identifying both hurdles that have been already overcome and outstanding issues. First, we describe the critical initial task of selecting which features of tissue microstructure can be estimated using a model and which acquisition protocol needs to be implemented to make the estimation possible. The model performance should necessarily be tested in realistic numerical simulations and in experimental data - adapting the fitting strategy accordingly, and parameter estimates should be validated against complementary techniques, when/if available. Secondly, the model performance and validity should be explored in pathological conditions, and, if appropriate, dedicated models for pathology should be developed. We build on examples from tumors, ischemia and demyelinating diseases. We then discuss the challenges associated with clinical translation and added value. Finally, we single out four major unresolved challenges that are related to: the availability of a microstructural ground truth, the validation of model parameters which cannot be accessed with complementary techniques, the development of a generalized standard model for any brain region and pathology, and the seamless communication between different parties involved in the development and application of biophysical models of diffusion.

NODDI in clinical research

Diffusion MRI (dMRI) has proven to be a useful imaging approach for both clinical diagnosis and research investigating the microstructures of nervous tissues, and it has helped us to better understand the neurophysiological mechanisms of many diseases. Though diffusion tensor imaging (DTI) has long been the default tool to analyze dMRI data in clinical research, acquisition with stronger diffusion weightings beyond the DTI regimen is now possible with modern clinical scanners, potentially enabling even more detailed characterization of tissue microstructures. To take advantage of such data, neurite orientation dispersion and density imaging (NODDI) has been proposed as a way to relate the dMRI signal to tissue features via biophysically inspired modeling. The number of reports demonstrating the potential clinical utility of NODDI is rapidly increasing. At the same time, the pitfalls and limitations of NODDI, and general challenges in microstructure modeling, are becoming increasingly recognized by clinicians. dMRI microstructure modeling is a rapidly evolving field with great promise, where people from different scientific backgrounds, such as physics, medicine, biology, neuroscience, and statistics, are collaborating to build novel tools that contribute to improving human healthcare. Here, we review the applications of NODDI in clinical research and discuss future perspectives for investigations toward the implementation of dMRI microstructure imaging in clinical practice.

Microstructure with diffusion MRI: what scale we are sensitive to?

Diffusion-weighted MRI is the forerunner of the rapidly developed microstructure MRI (μMRI) aimed at in vivo evaluation of the cellular tissue architecture. This brief review focuses on the spatiotemporal scales of the microstructure that are accessible using different diffusion MRI techniques and the need to weight the measurability against the interpretability of results. Diffusion phenomena and models are first classified in two-dimensional space (the q-t-plane) of the measurement with narrow gradient pulses. Three-dimensional parameter space of the Stejskal-Tanner diffusion weighting adds more phenomena to this collection. Modern measurement techniques with larger number of parameters are briefly discussed under the overarching idea of diffusion weighting matching the geometry of the targeted cell species.

Direct and specific assessment of axonal injury and spinal cord microenvironments using diffusion correlation imaging

We describe a practical two-dimensional (2D) diffusion MRI framework to deliver specificity and improve sensitivity to axonal injury in the spinal cord. This approach provides intravoxel distributions of correlations of water mobilities in orthogonal directions, revealing sub-voxel diffusion components. Here we use it to investigate water diffusivities along axial and radial orientations within spinal cord specimens with confirmed, tract-specific axonal injury. First, we show using transmission electron microscopy and immunohistochemistry that tract-specific axonal beading occurs following Wallerian degeneration in the cortico-spinal tract as direct sequelae to closed head injury. We demonstrate that although some voxel-averaged diffusion tensor imaging (DTI) metrics are sensitive to this axonal injury, they are non-specific, i.e., they do not reveal an underlying biophysical mechanism of injury. Then we employ 2D diffusion correlation imaging (DCI) to improve discrimination of different water microenvironments by measuring and mapping the joint water mobility distributions perpendicular and parallel to the spinal cord axis. We determine six distinct diffusion spectral components that differ according to their microscopic anisotropy and mobility. We show that at the injury site a highly anisotropic diffusion component completely disappears and instead becomes more isotropic. Based on these findings, an injury-specific MR image of the spinal cord was generated, and a radiological-pathological correlation with histological silver staining % area was performed. The resulting strong and significant correlation (r = 0.70, p < 0.0001) indicates the high specificity with which DCI detects injury-induced tissue alterations. We predict that the ability to selectively image microstructural changes following axonal injury in the spinal cord can be useful in clinical and research applications by enabling specific detection and increased sensitivity to injury-induced microstructural alterations. These results also encourage us to translate DCI to higher spatial dimensions to enable assessment of traumatic axonal injury, and possibly other diseases and disorders in the brain.

Towards unconstrained compartment modeling in white matter using diffusion-relaxation MRI with tensor-valued diffusion encoding

PURPOSE

To optimize diffusion-relaxation MRI with tensor-valued diffusion encoding for precise estimation of compartment-specific fractions, diffusivities, and T₂ values within a two-compartment model of white matter, and to explore the approach in vivo.

METHODS

Sampling protocols featuring different b-values (b), b-tensor shapes (b_Δ ), and echo times (TE) were optimized using Cramér-Rao lower bounds (CRLB). Whole-brain data were acquired in children, adults, and elderly with white matter lesions. Compartment fractions, diffusivities, and T₂ values were estimated in a model featuring two microstructural compartments represented by a "stick" and a "zeppelin."

RESULTS

Precise parameter estimates were enabled by sampling protocols featuring seven or more "shells" with unique b/b_Δ /TE-combinations. Acquisition times were approximately 15 minutes. In white matter of adults, the "stick" compartment had a fraction of approximately 0.5 and, compared with the "zeppelin" compartment, featured lower isotropic diffusivities (0.6 vs. 1.3 μm² /ms) but higher T₂ values (85 vs. 65 ms). Children featured lower "stick" fractions (0.4). White matter lesions exhibited high "zeppelin" isotropic diffusivities (1.7 μm² /ms) and T₂ values (150 ms).

CONCLUSIONS

Diffusion-relaxation MRI with tensor-valued diffusion encoding expands the set of microstructure parameters that can be precisely estimated and therefore increases their specificity to biological quantities.

Nonivasive quantification of axon radii using diffusion MRI

Axon caliber plays a crucial role in determining conduction velocity and, consequently, in the timing and synchronization of neural activation. Noninvasive measurement of axon radii could have significant impact on the understanding of healthy and diseased neural processes. Until now, accurate axon radius mapping has eluded in vivo neuroimaging, mainly due to a lack of sensitivity of the MRI signal to micron-sized axons. Here, we show how - when confounding factors such as extra-axonal water and axonal orientation dispersion are eliminated - heavily diffusion-weighted MRI signals become sensitive to axon radii. However, diffusion MRI is only capable of estimating a single metric, the effective radius, representing the entire axon radius distribution within a voxel that emphasizes the larger axons. Our findings, both in rodents and humans, enable noninvasive mapping of critical information on axon radii, as well as resolve the long-standing debate on whether axon radii can be quantified.

On the need for bundle-specific microstructure kernels in diffusion MRI

Probing microstructure with diffusion magnetic resonance imaging (dMRI) on a scale orders of magnitude below the imaging resolution relies on biophysical modelling of the signal response in the tissue. The vast majority of these biophysical models of diffusion in white matter assume that the measured dMRI signal is the sum of the signals emanating from each of the constituent compartments, each of which exhibits a distinct behaviour in the b-value and/or orientation domain. Many of these models further assume that the dMRI behaviour of the oriented compartments (e.g. the intra-axonal space) is identical between distinct fibre populations, at least at the level of a single voxel. This implicitly assumes that any potential biological differences between fibre populations are negligible, at least as far as is measurable using dMRI. Here, we validate this assumption by means of a voxel-wise, model-free signal decomposition that, under the assumption above and in the absence of noise, is shown to be rank-1. We evaluate the effect size of signal components beyond this rank-1 representation and use permutation testing to assess their significance. We conclude that in the healthy adult brain, the dMRI signal is adequately represented by a rank-1 model, implying that biologically more realistic, but mathematically more complex fascicle-specific microstructure models do not capture statistically significant or anatomically meaningful structure, even in extended high-b diffusion MRI scans.