Fast and simplified mapping of mean axon diameter using temporal diffusion spectroscopy

Assessing tumor microstructure with time-dependent diffusion imaging: Considerations and feasibility on clinical MRI and MRI-Linac

BACKGROUND

Quantitative imaging biomarkers (QIBs) can characterize tumor heterogeneity and provide information for biological guidance in radiotherapy (RT). Time-dependent diffusion MRI (TDD-MRI) derived parameters are promising QIBs, as they describe tissue microstructure with more specificity than traditional diffusion-weighted MRI (DW-MRI). Specifically, TDD-MRI can provide information about both restricted diffusion and diffusional exchange, which are the two time-dependent effects affecting diffusion in tissue, and relevant in tumors. However, exhaustive modeling of both effects can require long acquisitions and complex model fitting. Furthermore, several introduced TDD-MRI measurements can require high gradient strengths and/or complex gradient waveforms that are possibly not available in RT settings.

PURPOSE

In this study, we investigated the feasibility of a simple analysis framework for the detection of restricted diffusion and diffusional exchange effects in the TDD-MRI signal. To promote the clinical applicability, we use standard gradient waveforms on a conventional 1.5 T MRI system with moderate gradient strength (Gmax = 45 mT/m), and on a hybrid 1.5 T MRI-Linac system with low gradient strength (Gmax = 15 mT/m).

METHODS

Restricted diffusion and diffusional exchange were simulated in geometries mimicking tumor microstructure to investigate the DW-MRI signal behavior and to determine optimal experimental parameters. TDD-MRI was implemented using pulsed field gradient spin echo with the optimized parameters on a conventional MRI system and a MRI-Linac. Experiments in green asparagus and 10 patients with brain lesions were performed to evaluate the time-dependent diffusion (TDD) contrast in the source DW-images.

RESULTS

Simulations demonstrated how the TDD contrast was able to differentiate only dominating diffusional exchange in smaller cells from dominating restricted diffusion in larger cells. The maximal TDD contrast in simulations with typical cancer cell sizes and in asparagus measurements exceeded 5% on the conventional MRI but remained below 5% on the MRI-Linac. In particular, the simulated TDD contrast in typical cancer cell sizes (r = 5-10 µm) remained below or around 2% with the MRI-Linac gradient strength. In patients measured with the conventional MRI, we found sub-regions reflecting either dominating restricted diffusion or dominating diffusional exchange in and around brain lesions compared to the noisy appearing white matter.

CONCLUSIONS

On the conventional MRI system, the TDD contrast maps showed consistent tumor sub-regions indicating different dominating TDD effects, potentially providing information on the spatial tumor heterogeneity. On the MRI-Linac, the available TDD contrast measured in asparagus showed the same trends as with the conventional MRI but remained close to typical measurement noise levels when simulated in common cancer cell sizes. On conventional MRI systems with moderate gradient strengths, the TDD contrast could potentially be used as a tool to identify which time-dependent effects to include when choosing a biophysical model for more specific tumor characterization.

Nonparametric distributions of tensor-valued Lorentzian diffusion spectra for model-free data inversion in multidimensional diffusion MRI

Magnetic resonance imaging (MRI) is the method of choice for noninvasive studies of micrometer-scale structures in biological tissues via their effects on the time- and frequency-dependent (restricted) and anisotropic self-diffusion of water. While new designs of time-dependent magnetic field gradient waveforms have enabled disambiguation between different aspects of translational motion that are convolved in traditional MRI methods relying on single pairs of field gradient pulses, data analysis for complex heterogeneous materials remains a challenge. Here, we propose and demonstrate nonparametric distributions of tensor-valued Lorentzian diffusion spectra, or "D(ω) distributions," as a general representation with sufficient flexibility to describe the MRI signal response from a wide range of model systems and biological tissues investigated with modulated gradient waveforms separating and correlating the effects of restricted and anisotropic diffusion.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency,ω $$ \omega $$ , in addition to the diffusion tensor,D $$ \mathbf{D} $$ , and relaxation,R 1 $$ {R}_1 $$ ,R 2 $$ {R}_2 $$ , correlations. AD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on theirD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency, ω , in addition to the diffusion tensor, D , and relaxation, R 1 , R 2 , correlations. A D ( ω ) - R 1 - R 2 clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on their D ( ω ) - R 1 - R 2 distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

Is Apparent Diffusion Coefficient Established as an Imaging Biomarker for Stereotactic Body Radiation Therapy Assessment in Hepatocellular Carcinoma?

ABSTRACT

In this article, as part of this special issue on biomarkers of early response, we review currently available reports regarding magnetic resonance imaging apparent diffusion coefficient (ADC) changes in hepatocellular carcinoma (HCC) in response to stereotactic body radiation therapy. We compare diffusion image acquisition, ADC analysis, methods for HCC response assessment, and statistical methods for prediction of local tumor progression by ADC metrics. We discuss the pros and cons of these studies. Following detailed analyses of existing investigations, we cannot conclude that ADC is established as an imaging biomarker for stereotactic body radiation therapy assessment in HCC.

Diffusion time dependency of extracellular diffusion

PURPOSE

To quantify the variations of the power-law dependences on diffusion time t or gradient frequencyf $$ f $$ of extracellular water diffusion measured by diffusion MRI (dMRI).

METHODS

Model cellular systems containing only extracellular water were used to investigate thet / f $$ t/f $$ dependence ofD ex $$ {D}_{ex} $$ , the extracellular diffusion coefficient. Computer simulations used a randomly packed tissue model with realistic intracellular volume fractions and cell sizes. DMRI measurements were performed on samples consisting of liposomes containing heavy water(D2 O, deuterium oxide) dispersed in regular water (H2 O).D ex $$ {D}_{ex} $$ was obtained over a broadt $$ t $$ range (∼1-1000 ms) and then fit power-law equationsD ex ( t ) = D const + const · t - ϑ t $$ {D}_{ex}(t)={D}_{\mathrm{const}}+\mathrm{const}\cdotp {t}^{-{\vartheta}_t} $$ andD ex ( f ) = D const + const · f ϑ f $$ {D}_{ex}(f)={D}_{\mathrm{const}}+\mathrm{const}\cdotp {f}^{\vartheta_f} $$ .

RESULTS

Both simulated and experimental results suggest that no single power-law adequately describes the behavior ofD ex $$ {D}_{ex} $$ over the range of diffusion times of most interest in practical dMRI. Previous theoretical predictions are accurate over only limitedt $$ t $$ ranges; for example,θ t = θ f = - 1 2 $$ {\theta}_t={\theta}_f=-\frac{1}{2} $$ is valid only for short times, whereasθ t = 1 $$ {\theta}_t=1 $$ orθ f = 3 2 $$ {\theta}_f=\frac{3}{2} $$ is valid only for long times but cannot describe other ranges simultaneously. For the specifict $$ t $$ range of 5-70 ms used in typical human dMRI measurements,θ t = θ f = 1 $$ {\theta}_t={\theta}_f=1 $$ matches the data well empirically.

CONCLUSION

The optimal power-law fit of extracellular diffusion varies with diffusion time. The dependency obtained at short or longt $$ t $$ limits cannot be applied to typical dMRI measurements in human cancer or liver. It is essential to determine the appropriate diffusion time range when modeling extracellular diffusion in dMRI-based quantitative microstructural imaging.

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.

Fourier transform temporal diffusion spectroscopy

Temporal diffusion spectroscopy (TDS) currently uses the oscillating gradient spin echo (OGSE) experiment to measure the spectral density of translational velocity autocorrelation at single frequencies. Due to timing restrictions imposed by the transverse relaxation, the frequency selectivity and the sampling density of OGSE are limited, especially at low frequencies. We propose to overcome this problem by adopting the principles of Fourier transform spectroscopy. The new method of Fourier transform TDS (FTDS) uses two broadband gradient waveforms with different relative delays to make the spin echo attenuation sensitive to a broad range of diffusion frequencies with different harmonic modulations and calculates the spectrum by discrete Fourier transform. The method was validated by a measurement of diffusion spectra in highly restrictive tissues of a celery stalk and provided results consistent with OGSE, however, on a denser frequency grid.

Inversion-Recovery-Prepared Oscillating Gradient Sequence Improves Diffusion-Time Dependency Measurements in the Human Brain

BACKGROUND

Oscillating gradient diffusion MRI (dMRI) enables measurements at a short diffusion-time (t_d ), but it is challenging for clinical systems. Particularly, the low b-value and low resolution may give rise to cerebrospinal fluid (CSF) contamination.

PURPOSE

To assess the effect of CSF partial volume on t_d -dMRI measurements and efficacy of inversion-recovery (IR) prepared oscillating and pulsed gradient dMRI sequence to improve t_d -dMRI measurements in the human brain.

STUDY TYPE

Prospective.

SUBJECTS

Ten normal volunteers and six glioma patients.

FIELD STRENGTH/SEQUENCE

A 3 T; three-dimensional (3D) IR-prepared oscillating gradient-prepared gradient spin-echo (GRASE) and two-dimensional (2D) IR-prepared oscillating gradient echo-planar imaging (EPI) sequences.

ASSESSMENT

We assessed the t_d -dependent patterns of apparent diffusion coefficient (ADC) in several gray and white matter structures, including the hippocampal subfields (head, body, and tail), cortical gray matter, thalamus, and posterior white matter in normal volunteers. Pulsed gradient (0 Hz) and oscillating gradients at frequencies of 20 Hz, 40 Hz, and 60 Hz dMRI were acquired with GRASE and EPI sequences with or without the IR module. We also tested the t_d -dependency patterns in glioma patients using the EPI sequence with or without the IR module.

STATISTICAL TESTS

The differences in ADC across the different t_d s were compared by one-way ANOVA followed by post hoc pairwise t-tests with Bonferroni correction.

RESULTS

In the healthy subjects, brain regions that were possibly contaminated by CSF signals, such as the hippocampus (head, body, and tail) and cortical gray matter, t_d -dependent ADC changes were only significant with the IR-prepared 2D and 3D sequences but not with the non-IR sequences. In brain glioblastomas patients, significantly higher t_d -dependence was observed in the tumor region with the IR module than that without IR (slope = 0.0196 μm² /msec² vs. 0.0034 μm² /msec² ).

CONCLUSION

The IR-prepared sequence effectively suppressed the CSF partial volume effect and significantly improved the t_d -dependent measurements in the human brain.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 1.

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.

MR cell size imaging with temporal diffusion spectroscopy

Cytological features such as cell size and intracellular morphology provide fundamental information on cell status and hence may provide specific information on changes that arise within biological tissues. Such information is usually obtained by invasive biopsy in current clinical practice, which suffers several well-known disadvantages. Recently, novel MRI methods such as IMPULSED (imaging microstructural parameters using limited spectrally edited diffusion) have been developed for direct measurements of mean cell size non-invasively. The IMPULSED protocol is based on using temporal diffusion spectroscopy (TDS) to combine measurements of water diffusion over a wide range of diffusion times to probe cellular microstructure over varying length scales. IMPULSED has been shown to provide rapid, robust, and reliable mapping of mean cell size and is suitable for clinical imaging. More recently, cell size distributions have also been derived by appropriate analyses of data acquired with IMPULSED or similar sequences, which thus provides MRI-cytometry. This review summarizes the basic principles, practical implementations, validations, and example applications of MR cell size imaging based on TDS and demonstrates how cytometric information can be used in various applications. In addition, the limitations and potential future directions of MR cytometry are identified including the diagnosis of nonalcoholic steatohepatitis of the liver and the assessment of treatment response of cancers.

Probing neural tissues at small scales: Recent progress of oscillating gradient spin echo (OGSE) neuroimaging in humans

The detection sensitivity of diffusion MRI (dMRI) is dependent on diffusion times. A shorter diffusion time can increase the sensitivity to smaller length scales. However, the conventional dMRI uses the pulse gradient spin echo (PGSE) sequence that probes relatively long diffusion times only. To overcome this, the oscillating gradient spin echo (OGSE) sequence has been developed to probe much shorter diffusion times with hardware limitations on preclinical and clinical MRI systems. The OGSE sequence has been previously used on preclinical animal MRI systems. Recently, several studies have translated the OGSE sequence to humans on clinical MRI systems and achieved new information that is invisible using conventional PGSE sequence. This paper provides an overview of the recent progress of the OGSE neuroimaging in humans, including the technical improvements in the translation of the OGSE sequence to human imaging and various applications in different neurological disorders and stroke. Some possible future directions of the OGSE sequence are also discussed.

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.

Oscillating diffusion-encoding with a high gradient-amplitude and high slew-rate head-only gradient for human brain imaging

PURPOSE

We investigate the importance of high gradient-amplitude and high slew-rate on oscillating gradient spin echo (OGSE) diffusion imaging for human brain imaging and evaluate human brain imaging with OGSE on the MAGNUS head-gradient insert (200 mT/m amplitude and 500 T/m/s slew rate).

METHODS

Simulations with cosine-modulated and trapezoidal-cosine OGSE at various gradient amplitudes and slew rates were performed. Six healthy subjects were imaged with the MAGNUS gradient at 3T with OGSE at frequencies up to 100 Hz and b = 450 s/mm² . Comparisons were made against standard pulsed gradient spin echo (PGSE) diffusion in vivo and in an isotropic diffusion phantom.

RESULTS

Simulations show that to achieve high frequency and b-value simultaneously for OGSE, high gradient amplitude, high slew rates, and high peripheral nerve stimulation limits are required. A strong linear trend for increased diffusivity (mean: 8-19%, radial: 9-27%, parallel: 8-15%) was observed in normal white matter with OGSE (20 Hz to 100 Hz) as compared to PGSE. Linear fitting to frequency provided excellent correlation, and using a short-range disorder model provided radial long-term diffusivities of D_∞,MD = 911 ± 72 µm² /s, D_∞,PD = 1519 ± 164 µm² /s, and D_∞,RD = 640 ± 111 µm² /s and correlation lengths of l_{c ,MD} = 0.802 ± 0.156 µm, l_{c ,PD} = 0.837 ± 0.172 µm, and l_{c ,RD} = 0.780 ± 0.174 µm. Diffusivity changes with OGSE frequency were negligible in the phantom, as expected.

CONCLUSION

The high gradient amplitude, high slew rate, and high peripheral nerve stimulation thresholds of the MAGNUS head-gradient enables OGSE acquisition for in vivo human brain imaging.

Diffusion time dependency along the human corpus callosum and exploration of age and sex differences as assessed by oscillating gradient spin-echo diffusion tensor imaging

Conventional diffusion imaging uses pulsed gradient spin echo (PGSE) waveforms with diffusion times of tens of milliseconds (ms) to infer differences of white matter microstructure. The combined use of these long diffusion times with short diffusion times (<10 ​ms) enabled by oscillating gradient spin echo (OGSE) waveforms can enable more sensitivity to changes of restrictive boundaries on the scale of white matter microstructure (e.g. membranes reflecting the axon diameters). Here, PGSE and OGSE images were acquired at 4.7 ​T from 20 healthy volunteers aged 20-73 years (10 males). Mean, radial, and axial diffusivity, as well as fractional anisotropy were calculated in the genu, body and splenium of the corpus callosum (CC). Monte Carlo simulations were also conducted to examine the relationship of intra- and extra-axonal radial diffusivity with diffusion time over a range of axon diameters and distributions. The results showed elevated diffusivities with OGSE relative to PGSE in the genu and splenium (but not the body) in both males and females, but the OGSE-PGSE difference was greater in the genu for males. Females showed positive correlations of OGSE-PGSE diffusivity difference with age across the CC, whereas there were no such age correlations in males. Simulations of radial diffusion demonstrated that for axon sizes in human brain both OGSE and PGSE diffusivities were dominated by extra-axonal water, but the OGSE-PGSE difference nonetheless increased with area-weighted outer-axon diameter. Therefore, the lack of OGSE-PGSE difference in the body is not entirely consistent with literature that suggests it is composed predominantly of axons with large diameter. The greater OGSE-PGSE difference in the genu of males could reflect larger axon diameters than females. The OGSE-PGSE difference correlation with age in females could reflect loss of smaller axons at older ages. The use of OGSE with short diffusion times to sample the microstructural scale of restriction implies regional differences of axon diameters along the corpus callosum with preliminary results suggesting a dependence on age and sex.

Single- and double-Diffusion encoding MRI for studying ex vivo apparent axon diameter distribution in spinal cord white matter

Mapping average axon diameter (AAD) and axon diameter distribution (ADD) in neuronal tissues non-invasively is a challenging task that may have a tremendous effect on our understanding of the normal and diseased central nervous system (CNS). Water diffusion is used to probe microstructure in neuronal tissues, however, the different water populations and barriers that are present in these tissues turn this into a complex task. Therefore, it is not surprising that recently we have witnessed a burst in the development of new approaches and models that attempt to obtain, non-invasively, detailed microstructural information in the CNS. In this work, we aim at challenging and comparing the microstructural information obtained from single diffusion encoding (SDE) with double diffusion encoding (DDE) MRI. We first applied SDE and DDE MR spectroscopy (MRS) on microcapillary phantoms and then applied SDE and DDE MRI on an ex vivo porcine spinal cord (SC), using similar experimental conditions. The obtained diffusion MRI data were fitted by the same theoretical model, assuming that the signal in every voxel can be approximated as the superposition of a Gaussian-diffusing component and a series of restricted components having infinite cylindrical geometries. The diffusion MRI results were then compared with histological findings. We found a good agreement between the fittings and the experimental data in white matter (WM) voxels of the SC in both diffusion MRI methods. The microstructural information and apparent AADs extracted from SDE MRI were found to be similar or somewhat larger than those extracted from DDE MRI especially when the diffusion time was set to 40 ms. The apparent ADDs extracted from SDE and DDE MRI show reasonable agreement but somewhat weaker correspondence was observed between the diffusion MRI results and histology. The apparent subtle differences between the microstructural information obtained from SDE and DDE MRI are briefly discussed.

Diffusion dispersion imaging: Mapping oscillating gradient spin-echo frequency dependence in the human brain

PURPOSE

Oscillating gradient spin-echo (OGSE) diffusion MRI provides information about the microstructure of biological tissues by means of the frequency dependence of the apparent diffusion coefficient (ADC). ADC dependence on OGSE frequency has been explored in numerous rodent studies, but applications in the human brain have been limited and have suffered from low contrast between different frequencies, long scan times, and a limited exploration of the nature of the ADC dependence on frequency.

THEORY AND METHODS

Multiple frequency OGSE acquisitions were acquired in healthy subjects at 7T to explore the power-law frequency dependence of ADC, the "diffusion dispersion." Furthermore, a method for optimizing the estimation of the ADC difference between different OGSE frequencies was developed, which enabled the design of a highly efficient protocol for mapping diffusion dispersion.

RESULTS

For the first time, evidence of a linear dependence of ADC on the square root of frequency in healthy human white matter was obtained. Using the optimized protocol, high-quality, full-brain maps of apparent diffusion dispersion rate were also demonstrated at an isotropic resolution of 2 mm in a scan time of 6 min.

CONCLUSIONS

This work sheds light on the nature of diffusion dispersion in the healthy human brain and introduces full-brain diffusion dispersion mapping at clinically relevant scan times. These advances may lead to new biomarkers of pathology or improved microstructural modeling.

Diffusion MRI microstructural models in the cervical spinal cord - Application, normative values, and correlations with histological analysis

Multi-compartment tissue modeling using diffusion magnetic resonance imaging has proven valuable in the brain, offering novel indices sensitive to the tissue microstructural environment in vivo on clinical MRI scanners. However, application, characterization, and validation of these models in the spinal cord remain relatively under-studied. In this study, we apply a diffusion "signal" model (diffusion tensor imaging, DTI) and two commonly implemented "microstructural" models (neurite orientation dispersion and density imaging, NODDI; spherical mean technique, SMT) in the human cervical spinal cord of twenty-one healthy controls. We first provide normative values of DTI, SMT, and NODDI indices in a number of white matter ascending and descending pathways, as well as various gray matter regions. We then aim to validate the sensitivity and specificity of these diffusion-derived contrasts by relating these measures to indices of the tissue microenvironment provided by a histological template. We find that DTI indices are sensitive to a number of microstructural features, but lack specificity. The microstructural models also show sensitivity to a number of microstructure features; however, they do not capture the specific microstructural features explicitly modelled. Although often regarded as a simple extension of the brain in the central nervous system, it may be necessary to re-envision, or specifically adapt, diffusion microstructural models for application to the human spinal cord with clinically feasible acquisitions - specifically, adjusting, adapting, and re-validating the modeling as it relates to both theory (i.e. relevant biology, assumptions, and signal regimes) and parameter estimation (for example challenges of acquisition, artifacts, and processing).

Challenges in diffusion MRI tractography - Lessons learned from international benchmark competitions

Diffusion MRI (dMRI) fiber tractography has become a pillar of the neuroimaging community due to its ability to noninvasively map the structural connectivity of the brain. Despite widespread use in clinical and research domains, these methods suffer from several potential drawbacks or limitations. Thus, validating the accuracy and reproducibility of techniques is critical for sound scientific conclusions and effective clinical outcomes. Towards this end, a number of international benchmark competitions, or "challenges", has been organized by the diffusion MRI community in order to investigate the reliability of the tractography process by providing a platform to compare algorithms and results in a fair manner, and evaluate common and emerging algorithms in an effort to advance the state of the field. In this paper, we summarize the lessons from a decade of challenges in tractography, and give perspective on the past, present, and future "challenges" that the field of diffusion tractography faces.

Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation

We review, systematize and discuss models of diffusion in neuronal tissue, by putting them into an overarching physical context of coarse-graining over an increasing diffusion length scale. From this perspective, we view research on quantifying brain microstructure as occurring along three major avenues. The first avenue focusses on transient, or time-dependent, effects in diffusion. These effects signify the gradual coarse-graining of tissue structure, which occurs qualitatively differently in different brain tissue compartments. We show that transient effects contain information about the relevant length scales for neuronal tissue, such as the packing correlation length for neuronal fibers, as well as the degree of structural disorder along the neurites. The second avenue corresponds to the long-time limit, when the observed signal can be approximated as a sum of multiple nonexchanging anisotropic Gaussian components. Here, the challenge lies in parameter estimation and in resolving its hidden degeneracies. The third avenue employs multiple diffusion encoding techniques, able to access information not contained in the conventional diffusion propagator. We conclude with our outlook on future directions that could open exciting possibilities for designing quantitative markers of tissue physiology and pathology, based on methods of studying mesoscopic transport in disordered systems.

Microstructural imaging of the human brain with a 'super-scanner': 10 key advantages of ultra-strong gradients for diffusion MRI

The key component of a microstructural diffusion MRI 'super-scanner' is a dedicated high-strength gradient system that enables stronger diffusion weightings per unit time compared to conventional gradient designs. This can, in turn, drastically shorten the time needed for diffusion encoding, increase the signal-to-noise ratio, and facilitate measurements at shorter diffusion times. This review, written from the perspective of the UK National Facility for In Vivo MR Imaging of Human Tissue Microstructure, an initiative to establish a shared 300 mT/m-gradient facility amongst the microstructural imaging community, describes ten advantages of ultra-strong gradients for microstructural imaging. Specifically, we will discuss how the increase of the accessible measurement space compared to a lower-gradient systems (in terms of Δ, b-value, and TE) can accelerate developments in the areas of 1) axon diameter distribution mapping; 2) microstructural parameter estimation; 3) mapping micro-vs macroscopic anisotropy features with gradient waveforms beyond a single pair of pulsed-gradients; 4) multi-contrast experiments, e.g. diffusion-relaxometry; 5) tractography and high-resolution imaging in vivo and 6) post mortem; 7) diffusion-weighted spectroscopy of metabolites other than water; 8) tumour characterisation; 9) functional diffusion MRI; and 10) quality enhancement of images acquired on lower-gradient systems. We finally discuss practical barriers in the use of ultra-strong gradients, and provide an outlook on the next generation of 'super-scanners'.

Studying brain microstructure with magnetic susceptibility contrast at high-field

A rapidly developing application of high field MRI is the study of brain anatomy and function with contrast based on the magnetic susceptibility of tissues. To study the subtle variations in susceptibility contrast between tissues and with changes in brain activity, dedicated scan techniques such as susceptibility-weighted MRI and blood-oxygen level dependent functional MRI have been developed. Particularly strong susceptibility contrast has been observed with systems that operate at 7T and above, and their recent widespread use has led to an improved understanding of contributing sources and mechanisms. To interpret magnetic susceptibility contrast, analysis approaches have been developed with the goal of extracting measures that report on local tissue magnetic susceptibility, a physical quantity that, under certain conditions, allows estimation of blood oxygenation, local tissue iron content, and quantification of its changes with disease. Interestingly, high field studies have also brought to light that not only the makeup of tissues affects MRI susceptibility contrast, but that also a tissue's sub-voxel structure at scales all the way down to the molecular level plays an important role as well. In this review, various ways will be discussed by which sub-voxel structure can affect the MRI signal in general, and magnetic susceptibility in particular, sometimes in a complex fashion. In the light of this complexity, it appears likely that accurate, brain-wide quantification of iron will require the combination of multiple contrasts that may include diffusion and magnetization transfer information with susceptibility-weighted contrast. On the other hand, this complexity also offers opportunities to use magnetic susceptibility contrast to inform about specific microstructural aspects of brain tissue. Details and several examples will be presented in this review.

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

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

Can increased spatial resolution solve the crossing fiber problem for diffusion MRI?

It is now widely recognized that voxels with crossing fibers or complex geometrical configurations present a challenge for diffusion MRI (dMRI) reconstruction and fiber tracking, as well as microstructural modeling of brain tissues. This "crossing fiber" problem has been estimated to affect anywhere from 30% to as many as 90% of white matter voxels, and it is often assumed that increasing spatial resolution will decrease the prevalence of voxels containing multiple fiber populations. The aim of this study is to estimate the extent of the crossing fiber problem as we progressively increase the spatial resolution, with the goal of determining whether it is possible to mitigate this problem with higher resolution spatial sampling. This is accomplished using ex vivo MRI data of the macaque brain, followed by histological analysis of the same specimen to validate these measurements, as well as to extend this analysis to resolutions not yet achievable in practice with MRI. In both dMRI and histology, we find unexpected results: the prevalence of crossing fibers increases as we increase spatial resolution. The problem of crossing fibers appears to be a fundamental limitation of dMRI associated with the complexity of brain tissue, rather than a technical problem that can be overcome with advances such as higher fields and stronger gradients.

Validation of surface-to-volume ratio measurements derived from oscillating gradient spin echo on a clinical scanner using anisotropic fiber phantoms

A diffusion measurement in the short-time surface-to-volume ratio (S/V) limit (Mitra et al., Phys Rev Lett. 1992;68:3555) can disentangle the free diffusion coefficient from geometric restrictions to diffusion. Biophysical parameters, such as the S/V of tissue membranes, can be used to estimate microscopic length scales non-invasively. However, due to gradient strength limitations on clinical MRI scanners, pulsed gradient spin echo (PGSE) measurements are impractical for probing the S/V limit. To achieve this limit on clinical systems, an oscillating gradient spin echo (OGSE) sequence was developed. Two phantoms containing 10 fiber bundles, each consisting of impermeable aligned fibers with different packing densities, were constructed to achieve a range of S/V values. The frequency-dependent diffusion coefficient, D(ω), was measured in each fiber bundle using OGSE with different gradient waveforms (cosine, stretched cosine, and trapezoidal), while D(t) was measured from PGSE and stimulated-echo measurements. The S/V values derived from the universal high-frequency behavior of D(ω) were compared against those derived from quantitative proton density measurements using single spin echo (SE) with varying echo times, and from magnetic resonance fingerprinting (MRF). S/V estimates derived from different OGSE waveforms were similar and demonstrated excellent correlation with both SE- and MRF-derived S/V measures (ρ  ≥  0.99). Furthermore, there was a smoother transition between OGSE frequency f and PGSE diffusion time when using teffS/V=9/64f, rather than the commonly used t_eff  = 1/(4f), validating the specific frequency/diffusion time conversion for this regime. Our well-characterized fiber phantom can be used for the calibration of OGSE and diffusion modeling techniques, as the S/V ratio can be measured independently using other MR modalities. Moreover, our calibration experiment offers an exciting perspective of mapping tissue S/V on clinical systems.