TE dependent Diffusion Imaging (TEdDI) distinguishes between compartmental T2 relaxation times

Filter-exchange spectroscopy is sensitive to gradual cell membrane degradation

Transmembrane water permeability changes occur after initialization of necrosis and are a mechanism for early detection of cell death. Filter-exchange spectroscopy (FEXSY) is sensitive to transmembrane water permeability and enables its quantification by magnetic resonance via the apparent exchange rate (AXR). In this study, we investigate AXR changes during necrotic cell death. FEXSY measurements of yeast cells in different necrotic stages were performed and compared with established fluorescence cell death markers and pulsed gradient spin echo measurements. Furthermore, the influence of T2 relaxation on AXR was examined in a two-compartment system. The AXR of yeast cells increased slightly after incubation with 20% isopropanol, whereas it peaked sharply after incubation with 25% isopropanol. At this point, almost all the yeast cells were vital but showed compromised membranes. After incubation with 30% isopropanol, AXR measurements showed high variability, at a point corresponding to a majority of the yeast cells being in late-stage necrosis with disrupted cell membranes. Simulations revealed that, for FEXSY measurements in a two-compartment system, a long filter echo time (TEf), compared with the T2 of the slow-diffusing compartment, filters out a fraction of the slow-diffusing compartment signal and leads to overestimation of apparent diffusion coefficient (ADC) and underestimation of AXR. Our results demonstrate that AXR is sensitive to gradual permeabilization of the cell membrane of living cells in different permeabilization stages without exogenous contrast agents. AXR measurements were sensitive to permeability changes induced by relatively low concentrations of isopropanol, at levels for which no measurable effect was detectable by ADC measurements. TEf may act as a signal filter that affects the estimated AXR value of a system consisting of a variety of local diffusivities and a range of T2 that includes T2 values shorter or comparable with the TEf.

Constrained alternating minimization for parameter mapping (CAMP)

OBJECTIVE

To improve image quality in highly accelerated parameter mapping by incorporating a linear constraint that relates consecutive images.

APPROACH

In multi-echo T1 or T2 mapping, scan time is often shortened by acquiring undersampled but complementary measures of k-space at each TE or TI. However, residual undersampling artifacts from the individual images can then degrade the quality of the final parameter maps. In this work, a new reconstruction method, dubbed Constrained Alternating Minimization for Parameter mapping (CAMP), is introduced. This method simultaneously extracts T2 or T1* maps in addition to an image for each TE or TI from accelerated datasets, leveraging the constraints of the decay to improve the reconstructed image quality. The model enforces exponential decay through a linear constraint, resulting in a biconvex objective function that lends itself to alternating minimization. The method was tested in four in vivo volunteer experiments and validated in phantom studies and healthy subjects, using T2 and T1 mapping, with accelerations of up to 12.

MAIN RESULTS

CAMP is demonstrated for accelerated radial and Cartesian acquisitions in T2 and T1 mapping. The method is even applied to generate an entire T2 weighted image series from a single TSE dataset, despite the blockwise k-space sampling at each echo time. Experimental undersampled phantom and in vivo results processed with CAMP exhibit reduced artifacts without introducing bias.

SIGNIFICANCE

For a wide array of applications, CAMP linearizes the model cost function without sacrificing model accuracy so that the well-conditioned and highly efficient reconstruction algorithm improves the image quality of accelerated parameter maps.

Time-division multiplexing (TDM) sequence removes bias in T2 estimation and relaxation-diffusion measurements

Purpose

To compare the performance of multi-echo (ME) and time-division multiplexing (TDM) sequences for accelerated relaxation-diffusion MRI (rdMRI) acquisition and to examine their reliability in estimating accurate rdMRI microstructure measures.

Method

The ME, TDM, and the reference single-echo (SE) sequences with six echo times (TE) were implemented using Pulseq with single-band (SB-) and multi-band 2 (MB2-) acceleration factors. On a diffusion phantom, the image intensities of the three sequences were compared, and the differences were quantified using the normalized root mean squared error (NRMSE). For the in-vivo brain scan, besides the image intensity comparison and T2-estimates, different methods were used to assess sequence-related effects on microstructure estimation, including the relaxation diffusion imaging moment (REDIM) and the maximum-entropy relaxation diffusion distribution (MaxEnt-RDD).

Results

TDM performance was similar to the gold standard SE acquisition, whereas ME showed greater biases (3-4× larger NRMSEs for phantom, 2× for in-vivo). T2 values obtained from TDM closely matched SE, whereas ME sequences underestimated the T2 relaxation time. TDM provided similar diffusion and relaxation parameters as SE using REDIM, whereas SB-ME exhibited a 60% larger bias in the map and on average 3.5× larger bias in the covariance between relaxation-diffusion coefficients.

Conclusion

Our analysis demonstrates that TDM provides a more accurate estimation of relaxation-diffusion measurements while accelerating the acquisitions by a factor of 2 to 3.

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.

Normal Brain and Brain Tumor ADC: Changes Resulting From Variation of Diffusion Time and/or Echo Time in Pulsed-Gradient Spin Echo Diffusion Imaging

OBJECTIVES

Increasing gradient performance on modern magnetic resonance imaging scanners has profoundly reduced the attainable diffusion and echo times for clinically available pulsed-gradient spin echo (PGSE) sequences. This study investigated how this may impact the measured apparent diffusion coefficient (ADC), which is considered an important diagnostic marker for differentiation between normal and abnormal brain tissue and for therapeutic follow-up.

MATERIALS AND METHODS

Diffusion time and echo time dependence of the ADC were evaluated on a high-performance 3 T magnetic resonance imaging scanner. Diffusion PGSE brain scans were performed in 10 healthy volunteers and in 10 brain tumor patients using diffusion times of 16, 40, and 70 ms, echo times of 60, 75, and 104 ms at 3 b-values (0, 100, and 1000 s/mm 2 ), and a maximum gradient amplitude of 68 mT/m. A low gradient performance system was also emulated by reducing the diffusion encoding gradient amplitude to 19 mT/m. In healthy subjects, the ADC was measured in 6 deep gray matter regions and in 6 white matter regions. In patients, the ADC was measured in the solid part of the tumor.

RESULTS

With increasing diffusion time, a small but significant ADC increase of up to 2.5% was observed for 6 aggregate deep gray matter structures. With increasing echo time or reduced gradient performance, a small but significant ADC decrease of up to 2.6% was observed for 6 aggregate white matter structures. In tumors, diffusion time-related ADC changes were inconsistent without clear trend. For tumors with diffusivity above 1.0 μm 2 /ms, with prolonged echo time, there was a pronounced ADC increase of up to 12%. Meanwhile, for tumors with diffusivity at or below 1.0 μm 2 /ms, no change or a reduction was observed. Similar results were observed for gradient performance reduction, with an increase of up to 21%. The coefficient of variation determined in repeat experiments was 2.4%.

CONCLUSIONS

For PGSE and the explored parameter range, normal tissue ADC changes seem negligible. Meanwhile, observed tumor ADC changes can be relevant if ADC is used as a quantitative biomarker and not merely assessed by visual inspection. This highlights the importance of reporting all pertinent timing parameters in ADC studies and of considering these effects when building scan protocols for use in multicenter investigations.

Romer-EPTI: rotating-view motion-robust super-resolution EPTI for SNR-efficient distortion-free in-vivo mesoscale dMRI and microstructure imaging

Purpose

To overcome the major challenges in dMRI acquisition, including low SNR, distortion/blurring, and motion vulnerability.

Methods

A novel Romer-EPTI technique is developed to provide distortion-free dMRI with significant SNR gain, high motion-robustness, sharp spatial resolution, and simultaneous multi-TE imaging. It introduces a ROtating-view Motion-robust supEr-Resolution technique (Romer) combined with a distortion/blurring-free EPTI encoding. Romer enhances SNR by a simultaneous multi-thick-slice acquisition with rotating-view encoding, while providing high motion-robustness through a motion-aware super-resolution reconstruction, which also incorporates slice-profile and real-value diffusion, to resolve high-isotropic-resolution volumes. The in-plane encoding is performed using distortion/blurring-free EPTI, which further improves effective spatial resolution and motion robustness by preventing not only T2/T2*-blurring but also additional blurring resulting from combining encoded volumes with inconsistent geometries caused by dynamic distortions. Self-navigation was incorporated to enable efficient phase correction. Additional developments include strategies to address slab-boundary artifacts, achieve minimal TE for SNR gain at 7T, and achieve high robustness to strong phase variations at high b-values.

Results

Using Romer-EPTI, we demonstrate distortion-free whole-brain mesoscale in-vivo dMRI at both 3T (500-μm-iso) and 7T (485-μm-iso) for the first time, with high SNR efficiency (e.g., 25 × ), and high image quality free from distortion and slab-boundary artifacts with minimal blurring. Motion experiments demonstrate Romer-EPTI's high motion-robustness and ability to recover sharp images in the presence of motion. Romer-EPTI also demonstrates significant SNR gain and robustness in high b-value (b=5000s/mm2) and time-dependent dMRI.

Conclusion

Romer-EPTI significantly improves SNR, motion-robustness, and image quality, providing a highly efficient acquisition for high-resolution dMRI and microstructure imaging.

Exploiting gradient-echo frequency evolution: Probing white matter microstructure and extracting bulk susceptibility-induced frequency for quantitative susceptibility mapping

PURPOSE

This work is to investigate the microstructure-induced frequency shift in white matter (WM) with crossing fibers and to separate the microstructure-related frequency shift from the bulk susceptibility-induced frequency shift by model fitting the gradient-echo (GRE) frequency evolution for potentially more accurate quantitative susceptibility mapping (QSM).

METHODS

A hollow-cylinder fiber model (HCFM) with two fiber populations was developed to investigate GRE frequency evolutions in WM voxels with microstructural orientation dispersion. The simulated and experimentally measured TE-dependent local frequency shift was then fitted to a simplified frequency evolution model to obtain a microstructure-related frequency difference parameter (∆ f $$ \Delta f $$ ) and a TE-independent bulk susceptibility-induced frequency shift (C f $$ {C}_f $$ ). The obtainedC f $$ {C}_f $$ was then used for QSM reconstruction. Reconstruction performances were evaluated using a numerical head phantom and in vivo data and then compared to other multi-echo combination methods.

RESULTS

GRE frequency evolutions and∆ f $$ \Delta f $$ -based tissue parameters in both parallel and crossing fibers determined from our simulations were comparable to those observed in vivo. The TE-dependent frequency fitting method outperformed other multi-echo combination methods in estimatingC f $$ {C}_f $$ in simulations. The fitted∆ f $$ \Delta f $$ ,C f $$ {C}_f $$ , and QSM could be improved further by navigator-based B0 fluctuation correction.

CONCLUSION

A HCFM with two fiber populations can be used to characterize microstructure-induced frequency shifts in WM regions with crossing fibers. HCFM-based TE-dependent frequency fitting provides tissue contrast related to microstructure (∆ f $$ \Delta f $$ ) and in addition may help improve the quantification accuracy ofC f $$ {C}_f $$ and the corresponding QSM.

Maximizing SNR per unit time in diffusion MRI with multiband T-Hex spirals

PURPOSE

The characterization of tissue microstructure using diffusion MRI (dMRI) signals is rapidly evolving, with increasing sophistication of signal representations and microstructure models. However, this progress often requires signals to be acquired with very high b-values (e.g., b > 30 ms/μm2 ), along many directions, and using multiple b-values, leading to long scan times and extremely low SNR in dMRI images. The purpose of this work is to boost the SNR efficiency of dMRI by combining three particularly efficient spatial encoding techniques and utilizing a high-performance gradient system (Gmax  ≤ 300 mT/m) for efficient diffusion encoding.

METHODS

Spiral readouts, multiband imaging, and sampling on tilted hexagonal grids (T-Hex) are combined and implemented on a 3T MRI system with ultra-strong gradients. Image reconstruction is performed through an iterative cg-SENSE algorithm incorporating static off-resonance distributions and field dynamics as measured with an NMR field camera. Additionally, T-Hex multiband is combined with a more conventional EPI-readout and compared with state-of-the-art blipped-CAIPIRINHA sampling. The advantage of the proposed approach is furthermore investigated for clinically available gradient performance and diffusion kurtosis imaging.

RESULTS

High fidelity in vivo images with b-values up to 40 ms/μm2 are obtained. The approach provides superior SNR efficiency over other state-of-the-art multiband diffusion readout schemes.

CONCLUSION

The demonstrated gains hold promise for the widespread dissemination of advanced microstructural scans, especially in clinical populations.

Insights and improvements in correspondence between axonal volume fraction measured with diffusion-weighted MRI and electron microscopy

Biophysical diffusion-weighted imaging (DWI) models are increasingly used in neuroscience to estimate the axonal water fraction (f AW ), which in turn is key for noninvasive estimation of the axonal volume fraction (f A ). These models require thorough validation by comparison with a reference method, for example, electron microscopy (EM). While EM studies often neglect the unmyelinated axons and solely report the fraction of myelinated axons, in DWI both myelinated and unmyelinated axons contribute to the DWI signal. However, DWI models often include simplifications, for example, the neglect of differences in the compartmental relaxation times or fixed diffusivities, which in turn might affect the estimation off AW . We investigate whether linear calibration parameters (scaling and offset) can improve the comparability between EM- and DWI-based metrics off A . To this end, we (a) used six DWI models based on the so-called standard model of white matter (WM), including two models with fixed compartmental diffusivities (e.g., neurite orientation dispersion and density imaging, NODDI) and four models that fitted the compartmental diffusivities (e.g., white matter tract integrity, WMTI), and (b) used a multimodal data set including ex vivo diffusion DWI and EM data in mice with a broad dynamic range of fibre volume metrics. We demonstrated that the offset is associated with the volume fraction of unmyelinated axons and the scaling factor is associated with different compartmentalT 2 and can substantially enhance the comparability between EM- and DWI-based metrics off A . We found that DWI models that fitted compartmental diffusivities provided the most accurate estimates of the EM-basedf A . Finally, we introduced a more efficient hybrid calibration approach, where only the offset is estimated but the scaling is fixed to a theoretically predicted value. Using this approach, a similar one-to-one correspondence to EM was achieved for WMTI. The method presented can pave the way for use of validated DWI-based models in clinical research and neuroscience.

Assessment of Precision and Accuracy of Brain White Matter Microstructure using Combined Diffusion MRI and Relaxometry

Joint modeling of diffusion and relaxation has seen growing interest due to its potential to provide complementary information about tissue microstructure. For brain white matter, we designed an optimal diffusion-relaxometry MRI protocol that samples multiple b-values, B-tensor shapes, and echo times (TE). This variable-TE protocol (27 min) has as subsets a fixed-TE protocol (15 min) and a 2-shell dMRI protocol (7 min), both characterizing diffusion only. We assessed the sensitivity, specificity and reproducibility of these protocols with synthetic experiments and in six healthy volunteers. Compared with the fixed-TE protocol, the variable-TE protocol enables estimation of free water fractions while also capturing compartmental T 2 relaxation times. Jointly measuring diffusion and relaxation offers increased sensitivity and specificity to microstructure parameters in brain white matter with voxelwise coefficients of variation below 10%.

An open relaxation-diffusion MRI dataset in neurosurgical studies

Diffusion MRI (dMRI) is a safe and noninvasive technique that provides insight into the microarchitecture of brain tissue. Relaxation-diffusion MRI (rdMRI) is an extension of traditional dMRI that captures diffusion imaging data at multiple TEs to detect tissue heterogeneity between relaxation and diffusivity. rdMRI has great potential in neurosurgical research including brain tumor grading and treatment response evaluation. However, the lack of available data has limited the exploration of rdMRI in clinical settings. To address this, we are sharing a high-quality rdMRI dataset from 18 neurosurgical patients with different types of lesions, as well as two healthy individuals as controls. The rdMRI data was acquired using 7 TEs, where at each TE multi-shell dMRI with high spatial and angular resolutions is obtained at each TE. Each rdMRI scan underwent thorough artifact and distortion corrections using a specially designed processing pipeline. The dataset's quality was assessed using standard practices, including quality control and assurance. This resource is a valuable addition to neurosurgical studies, and all data are openly accessible.

Single-shot Echo Planar Time-resolved Imaging for multi-echo functional MRI and distortion-free diffusion imaging

Purpose

To develop EPTI, a multi-shot distortion-free multi-echo imaging technique, into a single-shot acquisition to achieve improved robustness to motion and physiological noise, increased temporal resolution, and high SNR efficiency for dynamic imaging applications.

Methods

A new spatiotemporal encoding was developed to achieve single-shot EPTI by enhancing spatiotemporal correlation in k-t space. The proposed single-shot encoding improves reconstruction conditioning and sampling efficiency, with additional optimization under various accelerations to achieve optimized performance. To achieve high SNR efficiency, continuous readout with minimized deadtime was employed that begins immediately after excitation and extends for an SNR-optimized length. Moreover, k-t partial Fourier and simultaneous multi-slice acquisition were integrated to further accelerate the acquisition and achieve high spatial and temporal resolution.

Results

We demonstrated that ss-EPTI achieves higher tSNR efficiency than multi-shot EPTI, and provides distortion-free imaging with densely-sampled multi-echo images at resolutions ~1.25-3 mm at 3T and 7T-with high SNR efficiency and with comparable temporal resolutions to ss-EPI. The ability of ss-EPTI to eliminate dynamic distortions common in EPI also further improves temporal stability. For fMRI, ss-EPTI also provides early-TE images (e.g., 2.9ms) to recover signal-intensity and functional-sensitivity dropout in challenging regions. The multi-echo images provide TE-dependent information about functional fluctuations, successfully distinguishing noise-components from BOLD signals and further improving tSNR. For diffusion MRI, ss-EPTI provides high-quality distortion-free diffusion images and multi-echo diffusion metrics.

Conclusion

ss-EPTI provides distortion-free imaging with high image quality, rich multi-echo information, and enhanced efficiency within comparable temporal resolution to ss-EPI, offering a robust and efficient acquisition for dynamic imaging.

Quantitative MRI reveals widespread, network-specific myelination change during generalized epilepsy progression

Activity-dependent myelination is a fundamental mode of brain plasticity which significantly influences network function. We recently discovered that absence seizures, which occur in multiple forms of generalized epilepsy, can induce activity-dependent myelination, which in turn promotes further progression of epilepsy. Structural alterations of myelin are likely to be widespread, given that absence seizures arise from an extensive thalamocortical network involving frontoparietal regions of the bilateral hemispheres. However, the temporal course and spatial extent of myelin plasticity is unknown, due to limitations of gold-standard histological methods such as electron microscopy (EM). In this study, we leveraged magnetization transfer and diffusion MRI for estimation of g-ratios across major white matter tracts in a mouse model of generalized epilepsy with progressive absence seizures. EM was performed on the same brains after MRI. After seizure progression, we found increased myelination (decreased g-ratios) throughout the anterior portion (genu-to-body) of the corpus callosum but not in the posterior portion (body-splenium) nor in the fornix or the internal capsule. Curves obtained from averaging g-ratio values at every longitudinal point of the corpus callosum were statistically different with p<0.001. Seizure-associated myelin differences found in the corpus callosum body with MRI were statistically significant (p = 0.0027) and were concordant with EM in the same region (p = 0.01). Notably, these differences were not detected by diffusion tensor imaging. This study reveals widespread myelin structural change that is specific to the absence seizure network. Furthermore, our findings demonstrate the potential utility and importance of MRI-based g-ratio estimation to non-invasively detect myelin plasticity.

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.

The impact of head orientation with respect to B0 on diffusion tensor MRI measures

Diffusion tensor MRI (DT-MRI) remains the most commonly used approach to characterise white matter (WM) anisotropy. However, DT estimates may be affected by tissue orientation w.r.t. B → 0 due to local gradients and intrinsic T 2 orientation dependence induced by the microstructure. This work aimed to investigate whether and how diffusion tensor MRI-derived measures depend on the orientation of the head with respect to the static magnetic field, B → 0 . By simulating WM as two compartments, we demonstrated that compartmental T 2 anisotropy can induce the dependence of diffusion tensor measures on the angle between WM fibres and the magnetic field. In in vivo experiments, reduced radial diffusivity and increased axial diffusivity were observed in white matter fibres perpendicular to B → 0 compared to those parallel to B → 0 . Fractional anisotropy varied by up to 20 % as a function of the angle between WM fibres and the orientation of the main magnetic field. To conclude, fibre orientation w.r.t. B → 0 is responsible for up to 7 % variance in diffusion tensor measures across the whole brain white matter from all subjects and head orientations. Fibre orientation w.r.t. B → 0 may introduce additional variance in clinical research studies using diffusion tensor imaging, particularly when it is difficult to control for (e.g., fetal or neonatal imaging, or when the trajectories of fibres change due to, e.g., space occupying lesions).

Fiber-orientation independent component of R2* obtained from single-orientation MRI measurements in simulations and a post-mortem human optic chiasm

The effective transverse relaxation rate (R2*) is sensitive to the microstructure of the human brain like the g-ratio which characterises the relative myelination of axons. However, the fibre-orientation dependence of R2* degrades its reproducibility and any microstructural derivative measure. To estimate its orientation-independent part (R2,iso*) from single multi-echo gradient-recalled-echo (meGRE) measurements at arbitrary orientations, a second-order polynomial in time model (hereafter M2) can be used. Its linear time-dependent parameter, β1, can be biophysically related to R2,iso* when neglecting the myelin water (MW) signal in the hollow cylinder fibre model (HCFM). Here, we examined the performance of M2 using experimental and simulated data with variable g-ratio and fibre dispersion. We found that the fitted β1 can estimate R2,iso* using meGRE with long maximum-echo time (TEmax ≈ 54 ms), but not accurately captures its microscopic dependence on the g-ratio (error 84%). We proposed a new heuristic expression for β1 that reduced the error to 12% for ex vivo compartmental R2 values. Using the new expression, we could estimate an MW fraction of 0.14 for fibres with negligible dispersion in a fixed human optic chiasm for the ex vivo compartmental R2 values but not for the in vivo values. M2 and the HCFM-based simulations failed to explain the measured R2*-orientation-dependence around the magic angle for a typical in vivo meGRE protocol (with TEmax ≈ 18 ms). In conclusion, further validation and the development of movement-robust in vivo meGRE protocols with TEmax ≈ 54 ms are required before M2 can be used to estimate R2,iso* in subjects.

Estimating axon radius using diffusion-relaxation MRI: calibrating a surface-based relaxation model with histology

Axon radius is a potential biomarker for brain diseases and a crucial tissue microstructure parameter that determines the speed of action potentials. Diffusion MRI (dMRI) allows non-invasive estimation of axon radius, but accurately estimating the radius of axons in the human brain is challenging. Most axons in the brain have a radius below one micrometer, which falls below the sensitivity limit of dMRI signals even when using the most advanced human MRI scanners. Therefore, new MRI methods that are sensitive to small axon radii are needed. In this proof-of-concept investigation, we examine whether a surface-based axonal relaxation process could mediate a relationship between intra-axonal T2 and T1 times and inner axon radius, as measured using postmortem histology. A unique in vivo human diffusion-T1-T2 relaxation dataset was acquired on a 3T MRI scanner with ultra-strong diffusion gradients, using a strong diffusion-weighting (i.e., b = 6,000 s/mm2) and multiple inversion and echo times. A second reduced diffusion-T2 dataset was collected at various echo times to evaluate the model further. The intra-axonal relaxation times were estimated by fitting a diffusion-relaxation model to the orientation-averaged spherical mean signals. Our analysis revealed that the proposed surface-based relaxation model effectively explains the relationship between the estimated relaxation times and the histological axon radius measured in various corpus callosum regions. Using these histological values, we developed a novel calibration approach to predict axon radius in other areas of the corpus callosum. Notably, the predicted radii and those determined from histological measurements were in close agreement.

C-NODDI: a constrained NODDI model for axonal density and orientation determinations in cerebral white matter

Purpose

Neurite orientation dispersion and density imaging (NODDI) provides measures of neurite density and dispersion through computation of the neurite density index (NDI) and the orientation dispersion index (ODI). However, NODDI overestimates the cerebrospinal fluid water fraction in white matter (WM) and provides physiologically unrealistic high NDI values. Furthermore, derived NDI values are echo-time (TE)-dependent. In this work, we propose a modification of NODDI, named constrained NODDI (C-NODDI), for NDI and ODI mapping in WM.

Methods

Using NODDI and C-NODDI, we investigated age-related alterations in WM in a cohort of 58 cognitively unimpaired adults. Further, NDI values derived using NODDI or C-NODDI were correlated with the neurofilament light chain (NfL) concentration levels, a plasma biomarker of axonal degeneration. Finally, we investigated the TE dependence of NODDI or C-NODDI derived NDI and ODI.

Results

ODI derived values using both approaches were virtually identical, exhibiting constant trends with age. Further, our results indicated a quadratic relationship between NDI and age suggesting that axonal maturation continues until middle age followed by a decrease. This quadratic association was notably significant in several WM regions using C-NODDI, while limited to a few regions using NODDI. Further, C-NODDI-NDI values exhibited a stronger correlation with NfL concentration levels as compared to NODDI-NDI, with lower NDI values corresponding to higher levels of NfL. Finally, we confirmed the previous finding that NDI estimation using NODDI was dependent on TE, while NDI derived values using C-NODDI exhibited lower sensitivity to TE in WM.

Conclusion

C-NODDI provides a complementary method to NODDI for determination of NDI in white matter.

High-resolution diffusion-weighted imaging at 7 Tesla: single-shot readout trajectories and their impact on signal-to-noise ratio, spatial resolution and accuracy

Diffusion MRI (dMRI) is a valuable imaging technique to study the connectivity and microstructure of the brain in vivo. However, the resolution of dMRI is limited by the low signal-to-noise ratio (SNR) of this technique. Various multi-shot acquisition strategies have been developed to achieve sub-millimeter resolution, but they require long scan times which can be restricting for patient scans. Alternatively, the SNR of single-shot acquisitions can be increased by using a spiral readout trajectory to minimize the sequence echo time. Imaging at ultra-high fields (UHF) could further increase the SNR of single-shot dMRI; however, the shorter T2* of brain tissue and the greater field non-uniformities at UHFs will degrade image quality, causing image blurring, distortions, and signal loss. In this study, we investigated the trade-off between the SNR and resolution of different k-space trajectories, including echo planar imaging (EPI), partial Fourier EPI, and spiral trajectories, over a range of dMRI resolutions at 7T. The effective resolution, spatial specificity and sharpening effect were measured from the point spread function (PSF) of the simulated diffusion sequences for a nominal resolution range of 0.6-1.8 mm. In-vivo partial brain scans at a nominal resolution of 1.5 mm isotropic were acquired using the three readout trajectories to validate the simulation results. Field probes were used to measure dynamic magnetic fields offline up to the 3^rd order of spherical harmonics. Image reconstruction was performed using static ΔB₀ field maps and the measured trajectories to correct image distortions and artifacts, leaving T2* effects as the primary source of blurring. The effective resolution was examined in fractional anisotropy (FA) maps calculated from a multi-shell dataset with b-values of 300, 1000, and 2000 s/mm² in 5, 16, and 48 directions, respectively. In-vivo scans at nominal resolutions of 1, 1.2, and 1.5 mm were acquired and the SNR of the different trajectories calculated using the multiple replica method to investigate the SNR. Finally, in-vivo whole brain scans with an effective resolution of 1.5 mm isotropic were acquired to explore the SNR and efficiency of different trajectories at a matching effective resolution. FA and intra-cellular volume fraction (ICVF) maps calculated using neurite orientation dispersion and density imaging (NODDI) were used for the comparison. The simulations and in vivo imaging results showed that for matching nominal resolutions, EPI trajectories had the highest specificity and effective resolution with maximum image sharpening effect. However, spirals have a significantly higher SNR, in particular at higher resolutions and even when the effective image resolutions are matched. Overall, this work shows that the higher SNR of single-shot spiral trajectories at 7T allows us to achieve higher effective resolutions compared to EPI and PF-EPI to map the microstructure and connectivity of small brain structures.

Integrating multimodality magnetic resonance imaging to the Allen Mouse Brain Common Coordinate Framework

High-resolution magnetic resonance imaging (MRI) affords unique image contrasts to nondestructively probe the tissue microstructure; validation of MRI findings with conventional histology is essential to better understand the MRI contrasts. However, the dramatic difference in the spatial resolution and image contrast of these two techniques impedes accurate comparison between MRI metrics and traditional histology. To better validate various MRI metrics, we acquired whole mouse brain multigradient recalled-echo and multishell diffusion MRI datasets at 25-μm isotropic resolution. The recently developed Allen Mouse Brain Common Coordinate Framework (CCFv3) provides opportunities to integrate multimodal and multiscale datasets of the whole mouse brain in a common three-dimensional (3D) space. The T2*, quantitative susceptibility mapping, diffusion tensor imaging, and neurite orientation dispersion and density imaging parameters were compared with both serial two-photon tomography images and 3D Nissl staining images in the CCFv3 at the same spatial resolution. The correlation between MRI and Nissl staining strongly depends on different metrics and different regions of the brain. Integrating different imaging modalities to the same space may substantially improve our understanding of the complexity of the brain at different scales.

Orientation dependent proton transverse relaxation in the human brain white matter: The magic angle effect on a cylindrical helix

PURPOSE

To overcome some limitations of previous proton orientation-dependent transverse relaxation formalisms in human brain white matter (WM) by a generalized magic angle effect function.

METHODS

A cylindrical helix model was developed embracing anisotropic rotational and translational diffusion of restricted molecules in WM, with the former characterized by an axially symmetric system. Transverse relaxation rates R₂ and R₂^∗ were divided into isotropic R₂ⁱ and anisotropic parts, R₂^a ∗ f(α,Φ - ε₀), with α denoting an open angle and ε₀ an orientation (Φ) offset from DTI-derived primary diffusivity direction. The proposed framework (Fit A) was compared to prior models without ε₀ on previously published water and methylene proton transverse relaxation rates from developing, healthy, and pathological WM at 3 T. Goodness of fit was represented by root-mean-square error (RMSE). F-test and linear correlation were used with statistical significance set to P ≤ 0.05.

RESULTS

Fit A significantly (P < 0.01) outperformed prior models as demonstrated by reduced RMSEs, e.g., 0.349 vs. 0.724 in myelin water. Fitted ε₀ was in good agreement with calculated ε₀ from directional diffusivities. Compared with those from healthy adult, the fitted R₂ⁱ, R₂^a, and α from neonates were substantially reduced but ε₀ increased, consistent with developing myelination. Significant positive (R₂ⁱ) and negative (α and R₂^a) correlations were found with aging (demyelination) in elderly.

CONCLUSION

The developed framework can better characterize orientation dependences from a wide range of proton transverse relaxation measurements in the human brain WM, thus shedding new light on myelin microstructural alterations at the molecular level.

Tensor denoising of multidimensional MRI data

PURPOSE

To develop a denoising strategy leveraging redundancy in high-dimensional data.

THEORY AND METHODS

The SNR fundamentally limits the information accessible by MRI. This limitation has been addressed by a host of denoising techniques, recently including the so-called MPPCA: principal component analysis of the signal followed by automated rank estimation, exploiting the Marchenko-Pastur distribution of noise singular values. Operating on matrices comprised of data patches, this popular approach objectively identifies noise components and, ideally, allows noise to be removed without introducing artifacts such as image blurring, or nonlocal averaging. The MPPCA rank estimation, however, relies on a large number of noise singular values relative to the number of signal components to avoid such ill effects. This condition is unlikely to be met when data patches and therefore matrices are small, for example due to spatially varying noise. Here, we introduce tensor MPPCA (tMPPCA) for the purpose of denoising multidimensional data, such as from multicontrast acquisitions. Rather than combining dimensions in matrices, tMPPCA uses each dimension of the multidimensional data's inherent tensor-structure to better characterize noise, and to recursively estimate signal components.

RESULTS

Relative to matrix-based MPPCA, tMPPCA requires no additional assumptions, and comparing the two in a numerical phantom and a multi-TE diffusion MRI data set, tMPPCA dramatically improves denoising performance. This is particularly true for small data patches, suggesting that tMPPCA can be especially beneficial in such cases.

CONCLUSIONS

The MPPCA denoising technique can be extended to high-dimensional data with improved performance for smaller patch sizes.

The distortions of the free water model for diffusion MRI data when assuming single compartment relaxometry and proton density

Objective.To document the bias of thesimplifiedfree water model of diffusion MRI (dMRI) signal vis-à-vis aspecificmodel which, in addition to diffusion, incorporates compartment-specific proton density (PD), T1 recovery during repetition time (TR), and T2 decay during echo time (TE).Approach.Both models assume that volume fractionfof the total signal in any voxel arises from the free water compartment (fw) such as cerebrospinal fluid or edema, and the remainder (1-f) from hindered water (hw) which is constrained by cellular structures such as white matter (WM). Thespecificandsimplifiedmodels are compared on a synthetic dataset, using a range of PD, T1 and T2 values. We then fit the models to anin vivohealthy brain dMRI dataset. For bothsyntheticandin vivodata we use experimentally feasible TR, TE, signal-to-noise ratio (SNR) and physiologically plausible diffusion profiles.Main results.From the simulations we see that the difference between the estimatedsimplified fandspecific fis largest for mid-range ground-truthf, and it increases as SNR increases. The estimation of volume fractionfis sensitive to the choice of model,simplifiedorspecific, but the estimated diffusion parameters are robust to small perturbations in the simulation.Specific fis more accurate and precise thansimplified f. In the white matter (WM) regions of thein vivoimages,specific fis lower thansimplified f.Significance.In dMRI models for free water, accounting for compartment specific PD, T1 and T2, in addition to diffusion, improves the estimation of model parameters. This extra model specification attenuates the estimation bias of compartmental volume fraction without affecting the estimation of other diffusion parameters.

Joint estimation of relaxation and diffusion tissue parameters for prostate cancer with relaxation-VERDICT MRI

This work presents a biophysical model of diffusion and relaxation MRI for prostate called relaxation vascular, extracellular and restricted diffusion for cytometry in tumours (rVERDICT). The model includes compartment-specific relaxation effects providing T1/T2 estimates and microstructural parameters unbiased by relaxation properties of the tissue. 44 men with suspected prostate cancer (PCa) underwent multiparametric MRI (mp-MRI) and VERDICT-MRI followed by targeted biopsy. We estimate joint diffusion and relaxation prostate tissue parameters with rVERDICT using deep neural networks for fast fitting. We tested the feasibility of rVERDICT estimates for Gleason grade discrimination and compared with classic VERDICT and the apparent diffusion coefficient (ADC) from mp-MRI. The rVERDICT intracellular volume fraction f_ic discriminated between Gleason 3 + 3 and 3 + 4 (p = 0.003) and Gleason 3 + 4 and ≥ 4 + 3 (p = 0.040), outperforming classic VERDICT and the ADC from mp-MRI. To evaluate the relaxation estimates we compare against independent multi-TE acquisitions, showing that the rVERDICT T2 values are not significantly different from those estimated with the independent multi-TE acquisition (p > 0.05). Also, rVERDICT parameters exhibited high repeatability when rescanning five patients (R² = 0.79-0.98; CV = 1-7%; ICC = 92-98%). The rVERDICT model allows for accurate, fast and repeatable estimation of diffusion and relaxation properties of PCa sensitive enough to discriminate Gleason grades 3 + 3, 3 + 4 and ≥ 4 + 3.

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.

Estimation of free water-corrected microscopic fractional anisotropy

Water diffusion anisotropy MRI is sensitive to microstructural changes in the brain that are hallmarks of various neurological conditions. However, conventional metrics like fractional anisotropy are confounded by neuron fiber orientation dispersion, and the relatively low resolution of diffusion-weighted MRI gives rise to significant free water partial volume effects in many brain regions that are adjacent to cerebrospinal fluid. Microscopic fractional anisotropy is a recent metric that can report water diffusion anisotropy independent of neuron fiber orientation dispersion but is still susceptible to free water contamination. In this paper, we present a free water elimination (FWE) technique to estimate microscopic fractional anisotropy and other related diffusion indices by implementing a signal representation in which the MRI signal within a voxel is assumed to come from two distinct sources: a tissue compartment and a free water compartment. A two-part algorithm is proposed to rapidly fit a set of diffusion-weighted MRI volumes containing both linear- and spherical-tensor encoding acquisitions to the representation. Simulations and in vivo acquisitions with four healthy volunteers indicated that the FWE method may be a feasible technique for measuring microscopic fractional anisotropy and other indices with greater specificity to neural tissue characteristics than conventional methods.

Minimal number of sampling directions for robust measures of the spherical mean diffusion weighted signal: Effects of sampling directions, b-value, signal-to-noise ratio, hardware, and fitting strategy

Several recent multi-compartment diffusion MRI investigations and modeling strategies have utilized the orientationally-averaged, or spherical mean, diffusion-weighted signal to study tissue microstructure of the central nervous system. Most experimental designs sample a large number of diffusion weighted directions in order to calculate the spherical mean signal, however, sampling a subset of these directions may increase scanning efficiency and enable either a decrease in scan time or the ability to sample more diffusion weightings. Here, we aim to determine the minimum number of gradient directions needed for a robust measurement of the spherical mean signal. We used computer simulations to characterize the variation of the measured spherical mean signal as a function of the number of gradient directions, while also investigating the effects of diffusion weighting (b-value), signal-to-noise ratio (SNR), available hardware, and spherical mean fitting strategy. We then utilize empirically acquired data in the brain and spinal cord to validate simulations, showing experimental results are in good agreement with simulations. We summarize these results by providing an intuitive lookup table to facilitate the determination of the minimal number of sampling directions needed for robust spherical mean measurements, and give recommendations based on SNR and experimental conditions.

Estimating axial diffusivity in the NODDI model

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

MR Fingerprinting with b-Tensor Encoding for Simultaneous Quantification of Relaxation and Diffusion in a Single Scan

PURPOSE

Although both relaxation and diffusion imaging are sensitive to tissue microstructure, studies have reported limited sensitivity and robustness of using relaxation or conventional diffusion alone to characterize tissue microstructure. Recently, it has been shown that tensor-valued diffusion encoding and joint relaxation-diffusion quantification enable more reliable quantification of compartment-specific microstructural properties. However, scan times to acquire such data can be prohibitive. Here, we aim to simultaneously quantify relaxation and diffusion using MR fingerprinting (MRF) and b-tensor encoding in a clinically feasible time.

METHODS

We developed multidimensional MRF scans (mdMRF) with linear and spherical b-tensor encoding (LTE and STE) to simultaneously quantify T1, T2, and ADC maps from a single scan. The image quality, accuracy, and scan efficiency were compared between the mdMRF using LTE and STE. Moreover, we investigated the robustness of different sequence designs to signal errors and their impact on the maps.

RESULTS

T1 and T2 maps derived from the mdMRF scans have consistently high image quality, while ADC maps are sensitive to different sequence designs. Notably, the fast imaging steady state precession (FISP)-based mdMRF scan with peripheral pulse gating provides the best ADC maps that are free of image distortion and shading artifacts.

CONCLUSION

We demonstrated the feasibility of quantifying T1, T2, and ADC maps simultaneously from a single mdMRF scan in around 24 s/slice. The map quality and quantitative values are consistent with the reference scans.

Reproducibility of the Standard Model of diffusion in white matter on clinical MRI systems

Estimating intra- and extra-axonal microstructure parameters, such as volume fractions and diffusivities, has been one of the major efforts in brain microstructure imaging with MRI. The Standard Model (SM) of diffusion in white matter has unified various modeling approaches based on impermeable narrow cylinders embedded in locally anisotropic extra-axonal space. However, estimating the SM parameters from a set of conventional diffusion MRI (dMRI) measurements is ill-conditioned. Multidimensional dMRI helps resolve the estimation degeneracies, but there remains a need for clinically feasible acquisitions that yield robust parameter maps. Here we find optimal multidimensional protocols by minimizing the mean-squared error of machine learning-based SM parameter estimates for two 3T scanners with corresponding gradient strengths of 40and80mT/m. We assess intra-scanner and inter-scanner repeatability for 15-minute optimal protocols by scanning 20 healthy volunteers twice on both scanners. The coefficients of variation all SM parameters except free water fraction are ≲10% voxelwise and 1-4% for their region-averaged values. As the achieved SM reproducibility outcomes are similar to those of conventional diffusion tensor imaging, our results enable robust in vivo mapping of white matter microstructure in neuroscience research and in the clinic.

Neurite Exchange Imaging (NEXI): A minimal model of diffusion in gray matter with inter-compartment water exchange

Biophysical models of diffusion in white matter have been center-stage over the past two decades and are essentially based on what is now commonly referred to as the "Standard Model" (SM) of non-exchanging anisotropic compartments with Gaussian diffusion. In this work, we focus on diffusion MRI in gray matter, which requires rethinking basic microstructure modeling blocks. In particular, at least three contributions beyond the SM need to be considered for gray matter: water exchange across the cell membrane - between neurites and the extracellular space; non-Gaussian diffusion along neuronal and glial processes - resulting from structural disorder; and signal contribution from soma. For the first contribution, we propose Neurite Exchange Imaging (NEXI) as an extension of the SM of diffusion, which builds on the anisotropic Kärger model of two exchanging compartments. Using datasets acquired at multiple diffusion weightings (b) and diffusion times (t) in the rat brain in vivo, we investigate the suitability of NEXI to describe the diffusion signal in the gray matter, compared to the other two possible contributions. Our results for the diffusion time window 20-45 ms show minimal diffusivity time-dependence and more pronounced kurtosis decay with time, which is well fit by the exchange model. Moreover, we observe lower signal for longer diffusion times at high b. In light of these observations, we identify exchange as the mechanism that best explains these signal signatures in both low-b and high-b regime, and thereby propose NEXI as the minimal model for gray matter microstructure mapping. We finally highlight multi-b multi-t acquisition protocols as being best suited to estimate NEXI model parameters reliably. Using this approach, we estimate the inter-compartment water exchange time to be 15 - 60 ms in the rat cortex and hippocampus in vivo, which is of the same order or shorter than the diffusion time in typical diffusion MRI acquisitions. This suggests water exchange as an essential component for interpreting diffusion MRI measurements in gray matter.

Assessing potential correlation between T2 relaxation and diffusion of lactate in the mouse brain

Purpose

While diffusion and T₂ relaxation are intertwined, little or no correlation exists between diffusion and T₂ relaxation of intracellular metabolites in the rodent brain, as measured by diffusion‐weighted MRS at different TEs. However, situation might be different for lactate, since it is present in both extracellular and intracellular spaces, which exhibit different diffusion properties and may also exhibit different T₂. Such a TE dependence would be crucial to account for when interpreting or modeling lactate diffusion. Here we propose to take advantage of a new diffusion sequence, where J‐modulation of lactate is canceled even at long TE, thus retaining excellent signal, to assess potential T₂ dependence on diffusion of lactate in the mouse brain.

Methods

Using a frequency‐selective diffusion‐weighted spin‐echo sequence that removes J‐modulation at 1.3 ppm, thus preserving lactate signal even at long TE, we investigate the effect of TE between 50.9 and 110.9 ms (while keeping diffusion time constant) on apparent diffusivity and kurtosis in the mouse brain.

Results

Regardless of the metabolites, no difference appears for the diffusion‐weighted signal attenuation with increasing TE. For lactate, apparent diffusivity and kurtosis remain unchanged as TE increases.

Conclusion

No significant TE dependence of diffusivity and kurtosis is measured for lactate in the 50–110 ms TE range, confirming that potential T₂ effects can be ignored when interpreting or modeling lactate diffusion.

Click here for author‐reader discussions

Performance of orientation distribution function-fingerprinting with a biophysical multicompartment diffusion model

PURPOSE

Orientation Distribution Function (ODF) peak finding methods typically fail to reconstruct fibers crossing at shallow angles below 40°, leading to errors in tractography. ODF-Fingerprinting (ODF-FP) with the biophysical multicompartment diffusion model allows for breaking this barrier.

METHODS

A randomized mechanism to generate a multidimensional ODF-dictionary that covers biologically plausible ranges of intra- and extra-axonal diffusivities and fraction volumes is introduced. This enables ODF-FP to address the high variability of brain tissue. The performance of the proposed approach is evaluated on both numerical simulations and a reconstruction of major fascicles from high- and low-resolution in vivo diffusion images.

RESULTS

ODF-FP with the suggested modifications correctly identifies fibers crossing at angles as shallow as 10 degrees in the simulated data. In vivo, our approach reaches 56% of true positives in determining fiber directions, resulting in visibly more accurate reconstruction of pyramidal tracts, arcuate fasciculus, and optic radiations than the state-of-the-art techniques. Moreover, the estimated diffusivity values and fraction volumes in corpus callosum conform with the values reported in the literature.

CONCLUSION

The modified ODF-FP outperforms commonly used fiber reconstruction methods at shallow angles, which improves deterministic tractography outcomes of major fascicles. In addition, the proposed approach allows for linearization of the microstructure parameters fitting problem.

Improved diffusion parameter estimation by incorporating T2 relaxation properties into the DKI-FWE model

The free water elimination (FWE) model and its kurtosis variant (DKI-FWE) can separate tissue and free water signal contributions, thus providing tissue-specific diffusional information. However, a downside of these models is that the associated parameter estimation problem is ill-conditioned, necessitating the use of advanced estimation techniques that can potentially bias the parameter estimates. In this work, we propose the T₂-DKI-FWE model that exploits the T₂ relaxation properties of both compartments, thereby better conditioning the parameter estimation problem and providing, at the same time, an additional potential biomarker (the T₂ of tissue). In our approach, the T₂ of tissue is estimated as an unknown parameter, whereas the T₂ of free water is assumed known a priori and fixed to a literature value (1573 ms). First, the error propagation of an erroneous assumption on the T₂ of free water is studied. Next, the improved conditioning of T₂-DKI-FWE compared to DKI-FWE is illustrated using the Cramér-Rao lower bound matrix. Finally, the performance of the T₂-DKI-FWE model is compared to that of the DKI-FWE and T₂-DKI models on both simulated and real datasets. The error due to a biased approximation of the T₂ of free water was found to be relatively small in various diffusion metrics and for a broad range of erroneous assumptions on its underlying ground truth value. Compared to DKI-FWE, using the T₂-DKI-FWE model is beneficial for the identifiability of the model parameters. Our results suggest that the T₂-DKI-FWE model can achieve precise and accurate diffusion parameter estimates, through effective reduction of free water partial volume effects and by using a standard nonlinear least squares approach. In conclusion, incorporating T₂ relaxation properties into the DKI-FWE model improves the conditioning of the model fitting, while only requiring an acquisition scheme with at least two different echo times.

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.

Accelerating joint relaxation-diffusion MRI by integrating time division multiplexing and simultaneous multi-slice (TDM-SMS) strategies

PURPOSE

To accelerate the acquisition of relaxation-diffusion imaging by integrating time-division multiplexing (TDM) with simultaneous multi-slice (SMS) for EPI and evaluate imaging quality and diffusion measures.

METHODS

The time-division multiplexing (TDM) technique and SMS method were integrated to achieve a high slice-acceleration (e.g., 6×) factor for acquiring relaxation-diffusion MRI. Two variants of the sequence, referred to as TDM3e-SMS and TDM2s-SMS, were developed to simultaneously acquire slice groups with three distinct TEs and two slice groups with the same TE, respectively. Both sequences were evaluated on a 3T scanner with in vivo human brains and compared with standard single-band (SB) -EPI and SMS-EPI using diffusion measures and tractography results.

RESULTS

Experimental results showed that the TDM3e-SMS sequence with total slice acceleration of 6 (multiplexing factor (MP) = 3 × multi-band factor (MB) = 2) provided similar image intensity and microstructure measures compared to standard SMS-EPI with MB = 2, and yielded less bias in intensity compared to standard SMS-EPI with MB = 4. The three sequences showed a similar positive correlation between TE and mean kurtosis (MK) and a negative correlation between TE and mean diffusivity (MD) in white matter. Multi-fiber tractography also shows consistency of results in TE-dependent measures between different sequences. The TDM2s-SMS sequence (MP = 2, MB = 2) also provided imaging measures similar to standard SMS-EPI sequences (MB = 2) for single-TE diffusion imaging.

CONCLUSIONS

The TDM-SMS sequence can provide additional 2× to 3× acceleration to SMS without degrading imaging quality. With the significant reduction in scan time, TDM-SMS makes joint relaxation-diffusion MRI a feasible technique in neuroimaging research to investigate new markers of brain disorders.

Distortion-Free Diffusion Imaging Using Self-Navigated Cartesian Echo-Planar Time Resolved Acquisition and Joint Magnitude and Phase Constrained Reconstruction

Echo-planar time resolved imaging (EPTI) is an effective approach for acquiring high-quality distortion-free images with a multi-shot EPI (ms-EPI) readout. As with traditional ms-EPI acquisitions, inter-shot phase variations present a main challenge when incorporating EPTI into a diffusion-prepared pulse sequence. The aim of this study is to develop a self-navigated Cartesian EPTI-based (scEPTI) acquisition together with a magnitude and phase constrained reconstruction for distortion-free diffusion imaging. A self-navigated Cartesian EPTI-based diffusion-prepared pulse sequence is designed. The different phase components in EPTI diffusion signal are analyzed and an approach to synthesize a fully phase-matched navigator for the inter-shot phase correction is demonstrated. Lastly, EPTI contains richer magnitude and phase information than conventional ms-EPI, such as the magnitude and phase correlations along the temporal dimension. The potential of these magnitude and phase correlations to enhance the reconstruction is explored. The reconstruction results with and without phase matching and with and without phase or magnitude constraints are compared. Compared with reconstruction without phase matching, the proposed phase matching method can improve the accuracy of inter-shot phase correction and reduce signal corruption in the final diffusion images. Magnitude constraints further improve image quality by suppressing the background noise and thereby increasing SNR, while phase constraints can mitigate possible image blurring from adding magnitude constraints. The high-quality distortion-free diffusion images and simultaneous diffusion-relaxometry imaging capacity provided by the proposed EPTI design represent a highly valuable tool for both clinical and neuroscientific assessments of tissue microstructure.

What's New and What's Next in Diffusion MRI Preprocessing

Diffusion MRI (dMRI) provides invaluable information for the study of tissue microstructure and brain connectivity, but suffers from a range of imaging artifacts that greatly challenge the analysis of results and their interpretability if not appropriately accounted for. This review will cover dMRI artifacts and preprocessing steps, some of which have not typically been considered in existing pipelines or reviews, or have only gained attention in recent years: brain/skull extraction, B-matrix incompatibilities w.r.t the imaging data, signal drift, Gibbs ringing, noise distribution bias, denoising, between- and within-volumes motion, eddy currents, outliers, susceptibility distortions, EPI Nyquist ghosts, gradient deviations, B₁ bias fields, and spatial normalization. The focus will be on “what’s new” since the notable advances prior to and brought by the Human Connectome Project (HCP), as presented in the predecessing issue on “Mapping the Connectome” in 2013. In addition to the development of novel strategies for dMRI preprocessing, exciting progress has been made in the availability of open source tools and reproducible pipelines, databases and simulation tools for the evaluation of preprocessing steps, and automated quality control frameworks, amongst others. Finally, this review will consider practical considerations and our view on “what’s next” in dMRI preprocessing.

The direction-dependence of apparent water exchange rate in human white matter

Transmembrane water exchange is a potential biomarker in the diagnosis and understanding of cancers, brain disorders, and other diseases. Filter-exchange imaging (FEXI), a special case of diffusion exchange spectroscopy adapted for clinical applications, has the potential to reveal different physiological water exchange processes. However, it is still controversial whether modulating the diffusion encoding gradient direction can affect the apparent exchange rate (AXR) measurements of FEXI in white matter (WM) where water diffusion shows strong anisotropy. In this study, we explored the diffusion-encoding direction dependence of FEXI in human brain white matter by performing FEXI with 20 diffusion-encoding directions on a clinical 3T scanner in-vivo. The results show that the AXR values measured when the gradients are perpendicular to the fiber orientation (0.77 ± 0.13 s ^- 1, mean ± standard deviation of all the subjects) are significantly larger than the AXR estimates when the gradients are parallel to the fiber orientation (0.33 ± 0.14 s ^- 1, p < 0.001) in WM voxels with coherently-orientated fibers. In addition, no significant correlation is found between AXRs measured along these two directions, indicating that they are measuring different water exchange processes. What's more, only the perpendicular AXR rather than the parallel AXR shows dependence on axonal diameter, indicating that the perpendicular AXR might reflect transmembrane water exchange between intra-axonal and extra-cellular spaces. Further finite difference (FD) simulations having three water compartments (intra-axonal, intra-glial, and extra-cellular spaces) to mimic WM micro-environments also suggest that the perpendicular AXR is more sensitive to the axonal water transmembrane exchange than parallel AXR. Taken together, our results show that AXR measured along different directions could be utilized to probe different water exchange processes in WM.

Nonparametric D-R1-R2 distribution MRI of the living human brain

Diffusion-relaxation correlation NMR can simultaneously characterize both the microstructure and the local chemical composition of complex samples that contain multiple populations of water. Recent developments on tensor-valued diffusion encoding and Monte Carlo inversion algorithms have made it possible to transfer diffusion-relaxation correlation NMR from small-bore scanners to clinical MRI systems. Initial studies on clinical MRI systems employed 5D D-R₁ and D-R₂ correlation to characterize healthy brain in vivo. However, these methods are subject to an inherent bias that originates from not including R₂ or R₁ in the analysis, respectively. This drawback can be remedied by extending the concept to 6D D-R₁-R₂ correlation. In this work, we present a sparse acquisition protocol that records all data necessary for in vivo 6D D-R₁-R₂ correlation MRI across 633 individual measurements within 25 min-a time frame comparable to previous lower-dimensional acquisition protocols. The data were processed with a Monte Carlo inversion algorithm to obtain nonparametric 6D D-R₁-R₂ distributions. We validated the reproducibility of the method in repeated measurements of healthy volunteers. For a post-therapy glioblastoma case featuring cysts, edema, and partially necrotic remains of tumor, we present representative single-voxel 6D distributions, parameter maps, and artificial contrasts over a wide range of diffusion-, R₁-, and R₂-weightings based on the rich information contained in the D-R₁-R₂ distributions.

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.

Neural networks for parameter estimation in microstructural MRI: Application to a diffusion-relaxation model of white matter

Specific features of white matter microstructure can be investigated by using biophysical models to interpret relaxation-diffusion MRI brain data. Although more intricate models have the potential to reveal more details of the tissue, they also incur time-consuming parameter estimation that may converge to inaccurate solutions due to a prevalence of local minima in a degenerate fitting landscape. Machine-learning fitting algorithms have been proposed to accelerate the parameter estimation and increase the robustness of the attained estimates. So far, learning-based fitting approaches have been restricted to microstructural models with a reduced number of independent model parameters where dense sets of training data are easy to generate. Moreover, the degree to which machine learning can alleviate the degeneracy problem is poorly understood. For conventional least-squares solvers, it has been shown that degeneracy can be avoided by acquisition with optimized relaxation-diffusion-correlation protocols that include tensor-valued diffusion encoding. Whether machine-learning techniques can offset these acquisition requirements remains to be tested. In this work, we employ artificial neural networks to vastly accelerate the parameter estimation for a recently introduced relaxation-diffusion model of white matter microstructure. We also develop strategies for assessing the accuracy and sensitivity of function fitting networks and use those strategies to explore the impact of the acquisition protocol. The developed learning-based fitting pipelines were tested on relaxation-diffusion data acquired with optimal and sub-optimal acquisition protocols. Networks trained with an optimized protocol were observed to provide accurate parameter estimates within short computational times. Comparing neural networks and least-squares solvers, we found the performance of the former to be less affected by sub-optimal protocols; however, model fitting networks were still susceptible to degeneracy issues and their use could not fully replace a careful design of the acquisition protocol.

Combined diffusion-relaxometry microstructure imaging: Current status and future prospects

Microstructure imaging seeks to noninvasively measure and map microscopic tissue features by pairing mathematical modeling with tailored MRI protocols. This article reviews an emerging paradigm that has the potential to provide a more detailed assessment of tissue microstructure-combined diffusion-relaxometry imaging. Combined diffusion-relaxometry acquisitions vary multiple MR contrast encodings-such as b-value, gradient direction, inversion time, and echo time-in a multidimensional acquisition space. When paired with suitable analysis techniques, this enables quantification of correlations and coupling between multiple MR parameters-such as diffusivity, T 1 , T 2 , and T 2 ∗ . This opens the possibility of disentangling multiple tissue compartments (within voxels) that are indistinguishable with single-contrast scans, enabling a new generation of microstructural maps with improved biological sensitivity and specificity.

Differentiation of symptomatic and asymptomatic carotid intraplaque hemorrhage using 3D high-resolution diffusion-weighted stack of stars imaging

Ischemic events related to carotid disease are far more strongly associated with plaque instability than stenosis. 3D high-resolution diffusion-weighted (DW) imaging can provide quantitative diffusion measurements on carotid atherosclerosis and may improve detection of vulnerable intraplaque hemorrhage (IPH). The 3D DW-stack of stars (SOS) sequence was implemented with 3D SOS acquisition combined with DW preparation. After simulation of signals created from 3D DW-SOS, phantom studies were performed. Three healthy subjects and 20 patients with carotid disease were recruited. Apparent diffusion coefficient (ADC) values were statistically analyzed on three subgroups by using a two-group comparison Wilcoxon-Mann-Whitney U test with p values less than 0.05: symptomatic versus asymptomatic; IPH-positive versus IPH-negative; and IPH-positive symptomatic versus asymptomatic plaques to determine the relationship with plaque vulnerability. ADC values calculated by 3D DW-SOS provided values similar to those calculated from other techniques. Mean ADC of symptomatic plaque was significantly lower than asymptomatic plaque (0.68 ± 0.18 vs. 0.98 ± 0.16 x 10-3  mm2 /s, p < 0.001). ADC was also significantly lower in IPH-positive versus IPH-negative plaque (0.68 ± 0.13 vs. 1.04 ± 0.11 x 10-3  mm2 /s, p < 0.001). Additionally, ADC was significantly lower in symptomatic versus asymptomatic IPH-positive plaque (0.57 ± 0.09 vs. 0.75 ± 0.11 x 10-3  mm2 /s, p < 0.001). Our results provide strong evidence that ADC measurements from 3D DW-SOS correlate with the symptomatic status of extracranial internal carotid artery plaque. Further, ADC improved discrimination of symptomatic plaque in IPH. These data suggest that diffusion characteristics may improve detection of destabilized plaque leading to elevated stroke risk.

Quantifying myelin in crossing fibers using diffusion-prepared phase imaging: Theory and simulations

PURPOSE

Myelin has long been the target of neuroimaging research. However, most available techniques can only provide a voxel-averaged estimate of myelin content. In the human brain, white matter fiber pathways connecting different brain areas and carrying different functions often cross each other in the same voxel. A measure that can differentiate the degree of myelination of crossing fibers would provide a more specific marker of myelination.

THEORY AND METHODS

One MRI signal property that is sensitive to myelin is the phase accumulation. This sensitivity is used by measuring the phase accumulation of the signal remaining after diffusion-weighting, which is called diffusion-prepared phase imaging (DIPPI). Including diffusion-weighting before estimating the phase accumulation has two distinct advantages for estimating the degree of myelination: (1) It increases the relative contribution of intra-axonal water, whose phase is related linearly to the thickness of the surrounding myelin (in particular the log g-ratio); and (2) it gives directional information, which can be used to distinguish between crossing fibers. Here the DIPPI sequence is described, an approach is proposed to estimate the log g-ratio, and simulations are used and DIPPI data acquired in an isotropic phantom to quantify other sources of phase accumulation.

RESULTS

The expected bias is estimated in the log g-ratio for reasonable in vivo acquisition parameters caused by eddy currents (~4%-10%), remaining extra-axonal signal (~15%), and gradients in the bulk off-resonance field (<10% for most of the brain).

CONCLUSION

This new sequence may provide a g-ratio estimate per fiber population crossing within a voxel.

Assessment of the effects of mimicking tissue microstructural properties on apparent diffusion coefficient and apparent exchange rate in diffusion MRI via a series of specially designed phantoms

PURPOSE

Diffusion MRI provides a valuable tool for imaging tissue microstructure. However, due to the lack of related experimental methods and specially designed phantoms, no experimental study has been conducted yet to quantitatively assess the effects of membrane permeability, intracellular volume fraction (IVF), and intracellular diffusivity on the apparent diffusion coefficient (ADC) obtained from diffusion weighted imaging (DWI), and the effects of membrane permeability on the apparent exchange rate (AXR) obtained from filter exchange imaging (FEXI).

METHODS

A series of phantoms with three adjustable parameters was designed to mimic tissue microstructural properties including membrane permeability, IVF, and intracellular diffusivity. Quantitative experiments were conducted to assess the effects of these properties on ADC and AXR. DWI scans were performed to obtain axial and radial ADC values. FEXI scans were performed to obtain AXR values.

RESULTS

Axial ADC values range from 1.148 μm² /ms to 2.157 μm² /ms, and radial ADC values range from 0.904 μm² /ms to 2.067 μm² /ms. Radial ADC decreased with a decrease in fiber permeability. Decreased axial and radial ADC values with increased intra-fiber volume fraction, and increased polyvinylpyrrolidone (PVP) concentration of the intra-fiber space were observed. AXR values range from 2.1 s^-1 to 4.9 s^-1 . AXR increases with fiber permeability.

CONCLUSION

The proposed phantoms can quantitatively evaluate the effects of mimicking tissue microstructural properties on ADC and AXR. This new phantom design provides a potential method for further understanding the biophysical mechanisms underlying the change in ADC and diffusion exchange.

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.

Diffusion in Sephadex Gel Structures: Time Dependency Revealed by Multi-Sequence Acquisition over a Broad Diffusion Time Range

It has been increasingly reported that in biological tissues diffusion-weighted MRI signal attenuation deviates from mono-exponential decay, especially at high b-values. A number of diffusion models have been proposed to characterize this non-Gaussian diffusion behavior. One of these models is the continuous-time random-walk (CTRW) model, which introduces two new parameters: a fractional order time derivative α and a fractional order spatial derivative β. These new parameters have been linked to intravoxel diffusion heterogeneities in time and space, respectively, and are believed to depend on diffusion times. Studies on this time dependency are limited, largely because the diffusion time cannot vary over a board range in a conventional spin-echo echo-planar imaging sequence due to the accompanying T2 decays. In this study, we investigated the time-dependency of the CTRW model in Sephadex gel phantoms across a broad diffusion time range by employing oscillating-gradient spin-echo, pulsed-gradient spin-echo, and pulsed-gradient stimulated echo sequences. We also performed Monte Carlo simulations to help understand our experimental results. It was observed that the diffusion process fell into the Gaussian regime at extremely short diffusion times whereas it exhibited a strong time dependency in the CTRW parameters at longer diffusion times.