Oscillating gradient diffusion MRI reveals unique microstructural information in normal and hypoxia-ischemia injured mouse brains

Low b-value Diffusion Tensor Imaging to Analyze the Dynamics of Cerebrospinal Fluid: Resolving Intravoxel Pseudorandom Motion into Ordered and Disordered Motions

PURPOSE

Analysis of cerebrospinal fluid (CSF) dynamics may be beneficial for understanding the mechanisms and diagnosis of several neurological diseases. Low b-value diffusion tensor imaging (low-b DTI) is useful for observing the slow and complex motion of the CSF. Theoretically, a mathematical framework suggests that low-b DTI provides the variance of the pseudorandom motion of the CSF. Furthermore, low-b DTI could provide comprehensive information on fluid dynamics. Accordingly, we proposed an analysis technique that resolves intravoxel pseudorandom motion into ordered (linear) and disordered (random) motions based on the mathematical framework.

METHODS

The proposed analysis technique helps measure low-b DTI with multiple diffusion times and linearly fits its mean diffusivity (MD) with the diffusion time to obtain two parameters, double-slope Vv and y-intersect Dr, which represent the variance of the velocity distribution of linear motion and the diffusion coefficient of random motion, respectively. Seven healthy subjects were scanned to evaluate the proposed technique and investigate fluid dynamics in several representative ROIs.

RESULTS

The obtained data showed the validity of the technique, repeatability, and consistency across the subjects in ROIs, such as the lateral ventricle (LV), third ventricle (3V), fourth ventricle (4V), and Sylvian fissure (SF). The obtained parameters Vv and Dr highlighted different characteristics of fluid dynamics in the representative ROIs: low Vv and low Dr in the LV, high Vv and moderate Dr in the 3V, and moderate Vv and moderate Dr in the 4V and SF.

CONCLUSION

The proposed analysis technique will facilitate a comprehensive investigation of the complex dynamics of the CSF using resolved parameters representing ordered and disordered motions.

In vivo rat brain mapping of multiple gray matter water populations using nonparametric D(ω)-R1-R2 distributions MRI

Massively multidimensional diffusion magnetic resonance imaging combines tensor-valued encoding, oscillating gradients, and diffusion-relaxation correlation to provide multicomponent subvoxel parameters depicting some tissue microstructural features. This method was successfully implemented ex vivo in microimaging systems and clinical conditions with tensor-valued gradient waveform of variable duration giving access to a narrow diffusion frequency (ω) range. We demonstrate here its preclinical in vivo implementation with a protocol of 389 contrast images probing a wide diffusion frequency range of 18 to 92 Hz at b-values up to 2.1 ms/μm2 enabled by the use of modulated gradient waveforms and combined with multislice high-resolution and low-distortion echo planar imaging acquisition with segmented and full reversed phase-encode acquisition. This framework allows the identification of diffusion ω-dependence in the rat cerebellum and olfactory bulb gray matter (GM), and the parameter distributions are shown to resolve two water pools in the cerebellum GM with different diffusion coefficients, shapes, ω-dependence, relaxation rates, and spatial repartition whose attribution to specific microstructure could modify the current understanding of the origin of restriction in GM.

In vivo mapping of cellular resolution neuropathology in brain ischemia with diffusion MRI

Noninvasive mapping of cellular pathology can provide critical diagnostic and prognostic information. Recent advances in diffusion magnetic resonance imaging enabled in vivo examination of tissue microstructures well beyond the imaging resolution. Here, we proposed to use diffusion time-dependent diffusion kurtosis imaging (tDKI) to simultaneously assess cellular morphology and transmembrane permeability in hypoxic-ischemic (HI) brain injury. Through numerical simulations and organoid imaging, we demonstrated the feasibility of capturing effective size and permeability changes using tDKI. In vivo MRI of HI-injured mouse brains detected a shift of the tDKI peak to longer diffusion times, suggesting swelling of the cellular processes. Furthermore, we observed a faster decrease of the tDKI tail, reflecting increased transmembrane permeability associated with up-regulated water exchange or necrosis. Such information, unavailable from a single diffusion time, can predict salvageable tissues. Preliminary applications of tDKI in patients with ischemic stroke suggested increased transmembrane permeability in stroke regions, illustrating tDKI's potential for detecting pathological changes in the clinics.

Tensor-valued diffusion MRI of human acute stroke

PURPOSE

Tensor-valued diffusion encoding can disentangle orientation dispersion and subvoxel anisotropy, potentially offering insight into microstructural changes after cerebral ischemia. The purpose was to evaluate tensor-valued diffusion MRI in human acute ischemic stroke, assess potential confounders from diffusion time dependencies, and compare to Monte Carlo diffusion simulations of axon beading.

METHODS

Linear (LTE) and spherical (STE) b-tensor encoding with inherently different effective diffusion times were acquired in 21 acute ischemic stroke patients between 3 and 57 h post-onset at 3 T in 2.5 min. In an additional 10 patients, STE with 2 LTE yielding different effective diffusion times were acquired for comparison. Diffusional variance decomposition (DIVIDE) was used to estimate microscopic anisotropy (μFA), as well as anisotropic, isotropic, and total diffusional variance (MKA , MKI , MKT ). DIVIDE parameters, and diffusion tensor imaging (DTI)-derived mean diffusivity and fractional anisotropy (FA) were compared in lesion versus contralateral white matter. Monte Carlo diffusion simulations of various cylindrical geometries for all b-tensor protocols were used to interpret parameter measurements.

RESULTS

MD was ˜40% lower in lesions for all LTE/STE protocols. The DIVIDE parameters varied with effective diffusion time: higher μFA and MKA in lesion versus contralateral white matter for STE with longer effective diffusion time LTE, whereas the shorter effective diffusion time LTE protocol yielded lower μFA and MKA in lesions. Both protocols, regardless of diffusion time, were consistent with simulations of greater beading amplitude and intracellular volume fraction.

CONCLUSION

DIVIDE parameters depend on diffusion time in acute stroke but consistently indicate neurite beading and larger intracellular volume fraction.

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

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

Differentiating primary central nervous system lymphoma from glioblastoma by time-dependent diffusion using oscillating gradient

BACKGROUND

This study aimed to elucidate the impact of effective diffusion time setting on apparent diffusion coefficient (ADC)-based differentiation between primary central nervous system lymphomas (PCNSLs) and glioblastomas (GBMs) and to investigate the usage of time-dependent diffusion magnetic resonance imaging (MRI) parameters.

METHODS

A retrospective study was conducted involving 21 patients with PCNSLs and 66 patients with GBMs using diffusion weighted imaging (DWI) sequences with oscillating gradient spin-echo (Δeff = 7.1 ms) and conventional pulsed gradient (Δeff = 44.5 ms). In addition to ADC maps at the two diffusion times (ADC7.1 ms and ADC44.5 ms), we generated maps of the ADC changes (cADC) and the relative ADC changes (rcADC) between the two diffusion times. Regions of interest were placed on enhancing regions and non-enhancing peritumoral regions. The mean and the fifth and 95th percentile values of each parameter were compared between PCNSLs and GBMs. The area under the receiver operating characteristic curve (AUC) values were used to compare the discriminating performances among the indices.

RESULTS

In enhancing regions, the mean and fifth and 95th percentile values of ADC44.5 ms and ADC7.1 ms in PCNSLs were significantly lower than those in GBMs (p = 0.02 for 95th percentile of ADC44.5 ms, p = 0.04 for ADC7.1 ms, and p < 0.01 for others). Furthermore, the mean and fifth and 95th percentile values of cADC and rcADC were significantly higher in PCNSLs than in GBMs (each p < 0.01). The AUC of the best-performing index for ADC7.1 ms was significantly lower than that for ADC44.5 ms (p < 0.001). The mean rcADC showed the highest discriminating performance (AUC = 0.920) among all indices. In peritumoral regions, no significant difference in any of the three indices of ADC44.5 ms, ADC7.1 ms, cADC, and rcADC was observed between PCNSLs and GBMs.

CONCLUSIONS

Effective diffusion time setting can have a crucial impact on the performance of ADC in differentiating between PCNSLs and GBMs. The time-dependent diffusion MRI parameters may be useful in the differentiation of these lesions.

Frequency-dependent diffusion kurtosis imaging in the human brain using an oscillating gradient spin echo sequence and a high-performance head-only gradient

Measuring the time/frequency dependence of diffusion MRI is a promising approach to distinguish between the effects of different tissue microenvironments, such as membrane restriction, tissue heterogeneity, and compartmental water exchange. In this study, we measure the frequency dependence of diffusivity (D) and kurtosis (K) with oscillating gradient diffusion encoding waveforms and a diffusion kurtosis imaging (DKI) model in human brains using a high-performance, head-only MAGNUS gradient system, with a combination of b-values, oscillating frequencies (f), and echo time that has not been achieved in human studies before. Frequency dependence of diffusivity and kurtosis are observed in both global and local white matter (WM) and gray matter (GM) regions and characterized with a power-law model ∼Λ*fθ. The frequency dependences of diffusivity and kurtosis (including changes between fmin and fmax, Λ, and θ) vary over different WM and GM regions, indicating potential microstructural differences between regions. A trend of decreasing kurtosis over frequency in the short-time limit is successfully captured for in vivo human brains. The effects of gradient nonlinearity (GNL) on frequency-dependent diffusivity and kurtosis measurements are investigated and corrected. Our results show that the GNL has prominent scaling effects on the measured diffusivity values (3.5∼5.5% difference in the global WM and 6∼8% difference in the global cortex) and subsequently affects the corresponding power-law parameters (Λ, θ) while having a marginal influence on the measured kurtosis values (<0.05% difference) and power-law parameters (Λ, θ). This study expands previous OGSE studies and further demonstrates the translatability of frequency-dependent diffusivity and kurtosis measurements to human brains, which may provide new opportunities to probe human brain microstructure in health and disease.

In vivo Mapping of Cellular Resolution Neuropathology in Brain Ischemia by Diffusion MRI

Non-invasive mapping of cellular pathology can provide critical diagnostic and prognostic information. Recent developments in diffusion MRI have produced new tools for examining tissue microstructure at a level well below the imaging resolution. Here, we report the use of diffusion time ( t )-dependent diffusion kurtosis imaging ( t DKI) to simultaneously assess the morphology and transmembrane permeability of cells and their processes in the context of pathological changes in hypoxic-ischemic brain (HI) injury. Through Monte Carlo simulations and cell culture organoid imaging, we demonstrate feasibility in measuring effective size and permeability changes based on the peak and tail of t DKI curves. In a mouse model of HI, in vivo imaging at 11.7T detects a marked shift of the t DKI peak to longer t in brain edema, suggesting swelling and beading associated with the astrocytic processes and neuronal neurites. Furthermore, we observed a faster decrease of the t DKI tail in injured brain regions, reflecting increased membrane permeability that was associated with upregulated water exchange upon astrocyte activation at acute stage as well as necrosis with disrupted membrane integrity at subacute stage. Such information, unavailable with conventional diffusion MRI at a single t, can predict salvageable tissues. For a proof-of-concept, t DKI at 3T on an ischemic stroke patient suggested increased membrane permeability in the stroke region. This work therefore demonstrates the potential of t DKI for in vivo detection of the pathological changes in microstructural morphology and transmembrane permeability after ischemic injury using a clinically translatable protocol.

Differentiating brain metastasis from glioblastoma by time-dependent diffusion MRI

BACKGROUND

This study was designed to investigate the use of time-dependent diffusion magnetic resonance imaging (MRI) parameters in distinguishing between glioblastomas and brain metastases.

METHODS

A retrospective study was conducted involving 65 patients with glioblastomas and 27 patients with metastases using a diffusion-weighted imaging sequence with oscillating gradient spin-echo (OGSE, 50 Hz) and a conventional pulsed gradient spin-echo (PGSE, 0 Hz) sequence. In addition to apparent diffusion coefficient (ADC) maps from two sequences (ADC50Hz and ADC0Hz), we generated maps of the ADC change (cADC): ADC50Hz - ADC0Hz and the relative ADC change (rcADC): (ADC50Hz - ADC0Hz)/ ADC0Hz × 100 (%).

RESULTS

The mean and the fifth and 95th percentile values of each parameter in enhancing and peritumoral regions were compared between glioblastomas and metastases. The area under the receiver operating characteristic curve (AUC) values of the best discriminating indices were compared. In enhancing regions, none of the indices of ADC0Hz and ADC50Hz showed significant differences between metastases and glioblastomas. The mean cADC and rcADC values of metastases were significantly higher than those of glioblastomas (0.24 ± 0.12 × 10-3mm2/s vs. 0.14 ± 0.03 × 10-3mm2/s and 23.3 ± 9.4% vs. 14.0 ± 4.7%; all p < 0.01). In peritumoral regions, no significant difference in all ADC indices was observed between metastases and glioblastomas. The AUC values for the mean cADC (0.877) and rcADC (0.819) values in enhancing regions were significantly higher than those for ADC0Hz5th (0.595; all p < 0.001).

CONCLUSIONS

The time-dependent diffusion MRI parameters may be useful for differentiating brain metastases from glioblastomas.

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.

Quantification of Diffusion Magnetic Resonance Imaging for Prognostic Prediction of Neonatal Hypoxic-Ischemic Encephalopathy

Neonatal hypoxic-ischemic encephalopathy (HIE) is the leading cause of acquired neonatal brain injury with the risk of developing serious neurological sequelae and death. An accurate and robust prediction of short- and long-term outcomes may provide clinicians and families with fundamental evidence for their decision-making, the design of treatment strategies, and the discussion of developmental intervention plans after discharge. Diffusion tensor imaging (DTI) is one of the most powerful neuroimaging tools with which to predict the prognosis of neonatal HIE by providing microscopic features that cannot be assessed by conventional magnetic resonance imaging (MRI). DTI provides various scalar measures that represent the properties of the tissue, such as fractional anisotropy (FA) and mean diffusivity (MD). Since the characteristics of the diffusion of water molecules represented by these measures are affected by the microscopic cellular and extracellular environment, such as the orientation of structural components and cell density, they are often used to study the normal developmental trajectory of the brain and as indicators of various tissue damage, including HIE-related pathologies, such as cytotoxic edema, vascular edema, inflammation, cell death, and Wallerian degeneration. Previous studies have demonstrated widespread alteration in DTI measurements in severe cases of HIE and more localized changes in neonates with mild-to-moderate HIE. In an attempt to establish cutoff values to predict the occurrence of neurological sequelae, MD and FA measurements in the corpus callosum, thalamus, basal ganglia, corticospinal tract, and frontal white matter have proven to have an excellent ability to predict severe neurological outcomes. In addition, a recent study has suggested that a data-driven, unbiased approach using machine learning techniques on features obtained from whole-brain image quantification may accurately predict the prognosis of HIE, including for mild-to-moderate cases. Further efforts are needed to overcome current challenges, such as MRI infrastructure, diffusion modeling methods, and data harmonization for clinical application. In addition, external validation of predictive models is essential for clinical application of DTI to prognostication.

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

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

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

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

Technical Advancements in Abdominal Diffusion-weighted Imaging

Since its first observation in the 18th century, the diffusion phenomenon has been actively studied by many researchers. Diffusion-weighted imaging (DWI) is a technique to probe the diffusion of water molecules and create a MR image with contrast based on the local diffusion properties. The DWI pixel intensity is modulated by the hindrance the diffusing water molecules experience. This hindrance is caused by structures in the tissue and reflects the state of the tissue. This characteristic makes DWI a unique and effective tool to gain more insight into the tissue's pathophysiological condition. In the past decades, DWI has made dramatic technical progress, leading to greater acceptance in clinical practice. In the abdominal region, however, acquiring DWI with good quality is challenging because of several reasons, such as large imaging volume, respiratory and other types of motion, and difficulty in achieving homogeneous fat suppression. In this review, we discuss technical advancements from the past decades that help mitigate these problems common in abdominal imaging. We describe the use of scan acceleration techniques such as parallel imaging and compressed sensing to reduce image distortion in echo planar imaging. Then we compare techniques developed to mitigate issues due to respiratory motion, such as free-breathing, respiratory-triggering, and navigator-based approaches. Commonly used fat suppression techniques are also introduced, and their effectiveness is discussed. Additionally, the influence of the abovementioned techniques on image quality is demonstrated. Finally, we discuss the current and future clinical applications of abdominal DWI, such as whole-body DWI, simultaneous multiple-slice excitation, intravoxel incoherent motion, and the use of artificial intelligence. Abdominal DWI has the potential to develop further in the future, thanks to scan acceleration and image quality improvement driven by technological advancements. The accumulation of clinical proof will further drive clinical acceptance.

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.

Investigation of time-dependent diffusion in extra-axial brain tumors using oscillating-gradient spin-echo

Oscillating gradient spin-echo (OGSE) sequences provide access to short diffusion times and may provide insight into micro-scale internal structures of pathologic lesions based on an analysis of changes in diffusivity with differing diffusion times. We hypothesized that changes in diffusivity acquired with a shorter diffusion time may permit elucidation of properties related to the internal structure of extra-axial brain tumors. This study aimed to investigate the utility of changes in diffusivity between short and long diffusion times for characterizing extra-axial brain tumors. In total, 12 patients with meningothelial meningiomas, 13 patients with acoustic neuromas, and 11 patients with pituitary adenomas were scanned with a 3 T magnetic resonance imaging (MRI) scanner with diffusion-weighted imaging (DWI) using OGSE and pulsed gradient spin-echo (PGSE) (effective diffusion times [Δ_eff]: 6.5 ms and 35.2 ms) with b-values of 0 and 1000 s/mm². Relative percentage changes between shorter and longer diffusion times were calculated using region-of-interest (ROI) analysis of brain tumors on λ₁, λ₂, λ₃, and mean diffusivity (MD) maps. The diffusivities of PGSE, OGSE, and relative percentage changes were compared among each tumor type using a multiple comparisons Steel-Dwass test. The mean (standard deviation) MD at Δ_eff of 6.5 ms was 1.07 ± 0.23 10^-3 mm²/s, 1.19 ± 0.18 10^-3 mm²/s, 1.19 ± 0.21 10^-3 mm²/s for meningothelial meningiomas, acoustic neuromas, and pituitary adenomas, respectively. The mean (standard deviation) MD at Δ_eff of 35.2 ms was 0.93 ± 0.22 10^-3 mm²/s, 1.07 ± 0.19 10^-3 mm²/s, 0.82 ± 0.21 10^-3 mm²/s for meningothelial meningiomas, acoustic neuromas, and pituitary adenomas, respectively. The mean (standard deviation) of the relative percentage change was 15.7 ± 4.4%, 12.4 ± 8.2%, 46.8 ± 11.3% for meningothelial meningiomas, acoustic neuromas, and pituitary adenomas, respectively. Compared to meningiomas and acoustic neuromas, pituitary adenoma exhibited stronger diffusion time-dependence with diffusion times between 6.5 ms and 35.2 ms (P < 0.05). In conclusion, differences in diffusion time-dependence may be attributed to differences in the internal structures of brain tumors. DWI with a short diffusion time may provide additional information on the microstructure of each tumor and contribute to tumor diagnosis.

Development of a composite diffusion tensor imaging score correlating with short-term neurological status in neonatal hypoxic-ischemic encephalopathy

Hypoxic-ischemic encephalopathy (HIE) is the most common cause of neonatal acquired brain injury. Although conventional MRI may predict neurodevelopmental outcomes, accurate prognostication remains difficult. As diffusion tensor imaging (DTI) may provide an additional diagnostic and prognostic value over conventional MRI, we aimed to develop a composite DTI (cDTI) score to relate to short-term neurological function. Sixty prospective neonates treated with therapeutic hypothermia (TH) for HIE were evaluated with DTI, with a voxel size of 1 × 1 × 2 mm. Fractional anisotropy (FA) and mean diffusivity (MD) from 100 neuroanatomical regions (FA/MD *100 = 200 DTI parameters in total) were quantified using an atlas-based image parcellation technique. A least absolute shrinkage and selection operator (LASSO) regression was applied to the DTI parameters to generate the cDTI score. Time to full oral nutrition [short-term oral feeding (STO) score] was used as a measure of short-term neurological function and was correlated with extracted DTI features. Seventeen DTI parameters were selected with LASSO and built into the final unbiased regression model. The selected factors included FA or MD values of the limbic structures, the corticospinal tract, and the frontotemporal cortices. While the cDTI score strongly correlated with the STO score (rho = 0.83, p = 2.8 × 10-16), it only weakly correlated with the Sarnat score (rho = 0.27, p = 0.035) and moderately with the NICHD-NRN neuroimaging score (rho = 0.43, p = 6.6 × 10-04). In contrast to the cDTI score, the NICHD-NRN score only moderately correlated with the STO score (rho = 0.37, p = 0.0037). Using a mixed-model analysis, interleukin-10 at admission to the NICU (p = 1.5 × 10-13) and tau protein at the end of TH/rewarming (p = 0.036) and after rewarming (p = 0.0015) were significantly associated with higher cDTI scores, suggesting that high cDTI scores were related to the intensity of the early inflammatory response and the severity of neuronal impairment after TH. In conclusion, a data-driven unbiased approach was applied to identify anatomical structures associated with some aspects of neurological function of HIE neonates after cooling and to build a cDTI score, which was correlated with the severity of short-term neurological functions.

Nonparametric 5D D-R2 distribution imaging with single-shot EPI at 21.1 T: Initial results for in vivo rat brain

In vivo human diffusion MRI is by default performed using single-shot EPI with greater than 50-ms echo times and associated signal loss from transverse relaxation. The individual benefits of the current trends of increasing B0 to boost SNR and employing more advanced signal preparation schemes to improve the specificity for selected microstructural properties eventually may be cancelled by increased relaxation rates at high B0 and echo times with advanced encoding. Here, initial attempts to translate state-of-the-art diffusion-relaxation correlation methods from 3 T to 21.1 T are made to identify hurdles that need to be overcome to fulfill the promises of both high SNR and readily interpretable microstructural information.

Impact of tissue properties on time-dependent alterations in apparent diffusion coefficient: a phantom study using oscillating-gradient spin-echo and pulsed-gradient spin-echo sequences

Purpose

The purpose of this study was to investigate whether the changes in apparent diffusion coefficients (ADCs) due to differences in diffusion time reflect tissue properties in actual measurements of phantoms.

Materials and methods

Various n-alkane phantoms and sucrose/collagen phantoms with various collagen densities were set up with and without polyvinyl alcohol (PVA) foam with an average pore diameter of 300 μm. Thus, n-alkanes or sucrose/collagen represented substrate viscosity and the presence of PVA foam represented tissue structure with septum. Diffusion-weighted images with various diffusion times (7.71–60 ms) were acquired using pulsed-gradient spin-echo (PGSE) and oscillating-gradient spin-echo (OGSE) sequences. The ADCs of the phantoms with and without PVA foam were calculated.

Results

The ADCs of some of the phantoms without PVA decreased with diffusion times decreased. In the n-alkane phantoms, only C₈H₁₈ showed significantly different ADCs depending on the use of PVA foam in the OGSE sequence. On the other hand, sucrose/collagen phantoms showed significant differences according to diffusion time. The ADCs of the phantoms decreased as the molecular size of the n-alkanes or collagen density of the sucrose/collagen phantom increased. Compared to phantoms without PVA foam, the ADC of the phantoms with PVA foam decreased as the diffusion time increased.

Conclusion

Changes in ADCs due to differences in diffusion time reflect tissue properties in actual measurements of phantoms. These changes in ADCs can be used for tissue characterization in vivo.

Investigation of breast cancer microstructure and microvasculature from time-dependent DWI and CEST in correlation with histological biomarkers

We investigated the associations of time-dependent DWI, non-Gaussian DWI, and CEST parameters with histological biomarkers in a breast cancer xenograft model. 22 xenograft mice (7 MCF-7 and 15 MDA-MB-231) were scanned at 4 diffusion times [T_d = 2.5/5 ms with 11 b-values (0–600 s/mm²) and T_d = 9/27.6 ms with 17 b-values (0–3000 s/mm²), respectively]. The apparent diffusion coefficient (ADC) was estimated using 2 b-values in different combinations (ADC^0–600 using b = 0 and 600 s/mm² and shifted ADC [sADC^200–1500] using b = 200 and 1500 s/mm²) at each of those diffusion times. Then the change (Δ) in ADC/sADC between diffusion times was evaluated. Non-Gaussian diffusion and intravoxel incoherent motion (IVIM) parameters (ADC₀, the virtual ADC at b = 0; K, Kurtosis from non-Gaussian diffusion; f, the IVIM perfusion fraction) were estimated. CEST images were acquired and the amide proton transfer signal intensity (APT SI) were measured. The ΔsADC_9–27.6 (between \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{sADC}}_{{9\,{\text{ms}}}}^{200{-}1500}$$\end{document}sADC9ms200-1500 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{sADC}}_{{27.6\,{\text{ms}}}}^{200{-}1500}$$\end{document}sADC27.6ms200-1500 and ΔADC_2.5_sADC_27.6 (between \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{ADC}}_{{2.5\, {\text{ms}}}}^{0{-}600}$$\end{document}ADC2.5ms0-600 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{sADC}}_{{27.6\,{\text{ms}}}}^{200{-}1500}$$\end{document}sADC27.6ms200-1500) was significantly larger for MCF-7 groups, and ΔADC_2.5_sADC_27.6 was positively correlated with Ki67_max and APT SI. ADC₀ decreased significantly in MDA-MB-231 group and K increased significantly with T_d in MCF-7 group. APT SI and cellular area had a moderately strong positive correlation in MDA-MB-231 and MCF-7 tumors combined, and there was a positive correlation in MDA-MB-231 tumors. There was a significant negative correlation between APT SI and the Ki-67-positive ratio in MDA-MB-231 tumors and when combined with MCF-7 tumors. The associations of ΔADC_2.5_sADC_27.6 and API SI with Ki-67 parameters indicate that the T_d-dependent DW and CEST parameters are useful to predict the histological markers of breast cancers.

Advanced Diffusion MR Imaging for Multiple Sclerosis in the Brain and Spinal Cord

Diffusion tensor imaging (DTI) has been established its usefulness in evaluating normal-appearing white matter (NAWM) and other lesions that are difficult to evaluate with routine clinical MRI in the evaluation of the brain and spinal cord lesions in multiple sclerosis (MS), a demyelinating disease. With the recent advances in the software and hardware of MRI systems, increasingly complex and sophisticated MRI and analysis methods, such as q-space imaging, diffusional kurtosis imaging, neurite orientation dispersion and density imaging, white matter tract integrity, and multiple diffusion encoding, referred to as advanced diffusion MRI, have been proposed. These are capable of capturing in vivo microstructural changes in the brain and spinal cord in normal and pathological states in greater detail than DTI.

This paper reviews the current status of recent advanced diffusion MRI for assessing MS in vivo as part of an issue celebrating two decades of magnetic resonance in medical sciences (MRMS), an official journal of the Japanese Society of Magnetic Resonance in Medicine.

Correlation Tensor MRI deciphers underlying kurtosis sources in stroke

Noninvasively detecting and characterizing modulations in cellular scale micro-architecture remains a desideratum for contemporary neuroimaging. Diffusion MRI (dMRI) has become the mainstay methodology for probing microstructure, and, in ischemia, its contrasts have revolutionized stroke management. Diffusion kurtosis imaging (DKI) has been shown to significantly enhance the sensitivity of stroke detection compared to its diffusion tensor imaging (DTI) counterparts. However, the interpretation of DKI remains ambiguous as its contrast may arise from competing kurtosis sources related to the anisotropy of tissue components, diffusivity variance across components, and microscopic kurtosis (e.g., arising from cross-sectional variance, structural disorder, and restriction). Resolving these sources may be fundamental for developing more specific imaging techniques for stroke management, prognosis, and understanding its pathophysiology. In this study, we apply Correlation Tensor MRI (CTI) - a double diffusion encoding (DDE) methodology recently introduced for deciphering kurtosis sources based on the unique information captured in DDE's diffusion correlation tensors - to investigate the underpinnings of kurtosis measurements in acute ischemic lesions. Simulations for the different kurtosis sources revealed specific signatures for cross-sectional variance (representing neurite beading), edema, and cell swelling. Ex vivo CTI experiments at 16.4 T were then performed in an experimental photothrombotic stroke model 3 h post-stroke (N = 10), and successfully separated anisotropic, isotropic, and microscopic non-Gaussian diffusion sources in the ischemic lesions. Each of these kurtosis sources provided unique contrasts in the stroked area. Particularly, microscopic kurtosis was shown to be a primary "driver" of total kurtosis upon ischemia; its large increases, coupled with decreases in anisotropic kurtosis, are consistent with the expected elevation in cross-sectional variance, likely linked to beading effects in small objects such as neurites. In vivo experiments at 9.4 T at the same time point (3 h post ischemia, N = 5) demonstrated the stability and relevance of the findings and showed that fixation is not a dominant confounder in our findings. In future studies, the different CTI contrasts may be useful to address current limitations of stroke imaging, e.g., penumbra characterization, distinguishing lesion progression form tissue recovery, and elucidating pathophysiological correlates.

Test-retest reproducibility of in vivo oscillating gradient and microscopic anisotropy diffusion MRI in mice at 9.4 Tesla

BACKGROUND AND PURPOSE

Microstructure imaging with advanced diffusion MRI (dMRI) techniques have shown increased sensitivity and specificity to microstructural changes in various disease and injury models. Oscillating gradient spin echo (OGSE) dMRI, implemented by varying the oscillating gradient frequency, and microscopic anisotropy (μA) dMRI, implemented via tensor valued diffusion encoding, may provide additional insight by increasing sensitivity to smaller spatial scales and disentangling fiber orientation dispersion from true microstructural changes, respectively. The aims of this study were to characterize the test-retest reproducibility of in vivo OGSE and μA dMRI metrics in the mouse brain at 9.4 Tesla and provide estimates of required sample sizes for future investigations.

METHODS

Twelve adult C57Bl/6 mice were scanned twice (5 days apart). Each imaging session consisted of multifrequency OGSE and μA dMRI protocols. Metrics investigated included μA, linear diffusion kurtosis, isotropic diffusion kurtosis, and the diffusion dispersion rate (Λ), which explores the power-law frequency dependence of mean diffusivity. The dMRI metric maps were analyzed with mean region-of-interest (ROI) and whole brain voxel-wise analysis. Bland-Altman plots and coefficients of variation (CV) were used to assess the reproducibility of OGSE and μA metrics. Furthermore, we estimated sample sizes required to detect a variety of effect sizes.

RESULTS

Bland-Altman plots showed negligible biases between test and retest sessions. ROI-based CVs revealed high reproducibility for most metrics (CVs < 15%). Voxel-wise CV maps revealed high reproducibility for μA (CVs ~ 10%), but low reproducibility for OGSE metrics (CVs ~ 50%).

CONCLUSION

Most of the μA dMRI metrics are reproducible in both ROI-based and voxel-wise analysis, while the OGSE dMRI metrics are only reproducible in ROI-based analysis. Given feasible sample sizes (10-15), μA metrics and OGSE metrics may provide sensitivity to subtle microstructural changes (4-8%) and moderate changes (> 6%), respectively.

Time-dependent Diffusion in Brain Abscesses Investigated with Oscillating-gradient Spin-echo

Oscillating-gradient spin-echo sequences enable the measurement of diffusion weighting with a short diffusion time and can provide indications of internal structures. We report two cases of brain abscess in which the apparent diffusion coefficient (ADC) values appear higher at short diffusion times in comparison with those at long diffusion times. Diffusion time dependence of the ADC in brain abscesses suggests not only substrate viscosity but also restricted diffusion due to the structure within the lesions.

Feasibility of oscillating and pulsed gradient diffusion MRI to assess neonatal hypoxia-ischemia on clinical systems

Diffusion-time- (t_d) dependent diffusion MRI (dMRI) extends our ability to characterize brain microstructure by measuring dMRI signals at varying t_d. The use of oscillating gradient (OG) is essential for accessing short t_d but is technically challenging on clinical MRI systems. This study aims to investigate the clinical feasibility and value of t_d-dependent dMRI in neonatal hypoxic-ischemic encephalopathy (HIE). Eighteen HIE neonates and six normal term-born neonates were scanned on a 3 T scanner, with OG-dMRI at an oscillating frequency of 33 Hz (equivalent t_d ≈ 7.5 ms) and pulsed gradient (PG)-dMRI at a t_d of 82.8 ms and b-value of 700 s/mm². The t_d-dependence, as quantified by the difference in apparent diffusivity coefficients between OG- and PG-dMRI (ΔADC), was observed in the normal neonatal brains, and the ΔADC was higher in the subcortical white matter than the deep grey matter. In HIE neonates with severe and moderate injury, ΔADC significantly increased in the basal ganglia (BG) compared to the controls (23.7% and 10.6%, respectively). In contrast, the conventional PG-ADC showed a 12.6% reduction only in the severe HIE group. White matter edema regions also demonstrated increased ΔADC, where PG-ADC did not show apparent changes. Our result demonstrated that t_d-dependent dMRI provided high sensitivity in detecting moderate-to-severe HIE.

Time-dependent diffusion MRI probes cerebellar microstructural alterations in a mouse model of Down syndrome

The cerebellum is a complex system with distinct cortical laminar organization. Alterations in cerebellar microstructure are common and associated with many factors such as genetics, cancer and ageing. Diffusion MRI (dMRI) provides a non-invasive tool to map the brain structural organization, and the recently proposed diffusion-time (td )-dependent dMRI further improves its capability to probe the cellular and axonal/dendritic microstructures by measuring water diffusion at multiple spatial scales. The td -dependent diffusion profile in the cerebellum and its utility in detecting cerebellar disorders, however, are not yet elucidated. Here, we first deciphered the spatial correspondence between dMRI contrast and cerebellar layers, based on which the cerebellar layer-specific td -dependent dMRI patterns were characterized in both euploid and Ts65Dn mice, a mouse model of Down syndrome. Using oscillating gradient dMRI, which accesses diffusion at short td 's by modulating the oscillating frequency, we detected subtle changes in the apparent diffusivity coefficient of the cerebellar internal granular layer and Purkinje cell layer of Ts65Dn mice that were not detectable by conventional pulsed gradient dMRI. The detection sensitivity of oscillating gradient dMRI increased with the oscillating frequency at both the neonatal and adult stages. The td -dependence, quantified by ΔADC map, was reduced in Ts65Dn mice, likely associated with the reduced granule cell density and abnormal dendritic arborization of Purkinje cells as revealed from histological evidence. Our study demonstrates superior sensitivity of short-td diffusion using oscillating gradient dMRI to detect cerebellar microstructural changes in Down syndrome, suggesting the potential application of this technique in cerebellar disorders.

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

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

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.

In vivo observation and biophysical interpretation of time-dependent diffusion in human cortical gray matter

The dependence of the diffusion MRI signal on the diffusion time t is a hallmark of tissue microstructure at the scale of the diffusion length. Here we measure the time-dependence of the mean diffusivity D(t) and mean kurtosis K(t) in cortical gray matter and in 25 ​gray matter sub-regions, in 10 healthy subjects. Significant diffusivity and kurtosis time-dependence is observed for t=21.2-100 ​ms, and is characterized by a power-law tail ∼t-ϑ with dynamical exponent ϑ. To interpret our measurements, we systematize the relevant scenarios and mechanisms for diffusion time-dependence in the brain. Using the effective medium theory formalism, we derive an exact relation between the power-law tails in D(t) and K(t). The estimated dynamical exponent ϑ≃1/2 in both D(t) and K(t) is consistent with one-dimensional diffusion in the presence of randomly positioned restrictions along neurites. We analyze the short-range disordered statistics of synapses on axon collaterals in the cortex, and perform one-dimensional Monte Carlo simulations of diffusion restricted by permeable barriers with a similar randomness in their placement, to confirm the ϑ=1/2 exponent. In contrast, the Kärger model of exchange is less consistent with the data since it does not capture the diffusivity time-dependence, and the estimated exchange time from K(t) falls below our measured t-range. Although we cannot exclude exchange as a contributing factor, we argue that structural disorder along neurites is mainly responsible for the observed time-dependence of diffusivity and kurtosis. Our observation and theoretical interpretation of the t-1/2 tail in D(t) and K(t) altogether establish the sensitivity of a macroscopic MRI signal to micrometer-scale structural heterogeneities along neurites in human gray matter in vivo.

Challenges for biophysical modeling of microstructure

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

Diffusion-prepared 3D gradient spin-echo sequence for improved oscillating gradient diffusion MRI

PURPOSE

Oscillating gradient (OG) enables the access of short diffusion times for time-dependent diffusion MRI (dMRI); however, it poses several technical challenges for clinical use. This study proposes a 3D oscillating gradient-prepared gradient spin-echo (OGprep-GRASE) sequence to improve SNR and shorten acquisition time for OG dMRI on clinical scanners.

METHODS

The 3D OGprep-GRASE sequence consisted of global saturation, diffusion encoding, fat saturation, and GRASE readout modules. Multiplexed sensitivity-encoding reconstruction was used to correct the phase errors between multiple shots. We compared the scan time and SNR of the proposed sequence and the conventional 2D-EPI sequence for OG dMRI at 30-90-mm slice coverage. We also examined the time-dependent diffusivity changes with OG dMRI acquired at frequencies of 50 Hz and 25 Hz and pulsed-gradient dMRI at diffusion times of 30 ms and 60 ms.

RESULTS

The OGprep-GRASE sequence reduced the scan time by a factor of 1.38, and increased the SNR by 1.74-2.27 times compared with 2D EPI for relatively thick slice coverage (60-90 mm). The SNR gain led to improved diffusion-tensor reconstruction in the multishot protocols. Image distortion in 2D-EPI images was also reduced in GRASE images. Diffusivity measurements from the pulsed-gradient dMRI and OG dMRI showed clear diffusion-time dependency in the white matter and gray matter of the human brain, using both the GRASE and EPI sequences.

CONCLUSION

The 3D OGprep-GRASE sequence improved scan time and SNR and reduced image distortion compared with the 2D multislice acquisition for OG dMRI on a 3T clinical system, which may facilitate the clinical translation of time-dependent dMRI.

Differentiation of high-grade and low-grade intra-axial brain tumors by time-dependent diffusion MRI

INTRODUCTION

Oscillating gradient spin-echo (OGSE) sequences enable acquisitions with shorter diffusion times. There is growing interest in the effect of diffusion time on apparent diffusion coefficient (ADC) values in patients with cancer. However, little evidence exists regarding its usefulness for differentiating between high-grade and low-grade brain tumors. The purpose of this study is to investigate the utility of changes in the ADC value between short and long diffusion times in distinguishing low-grade and high-grade brain tumors.

MATERIAL AND METHODS

Eleven patients with high-grade brain tumors and ten patients with low-grade brain tumors were scanned using a 3 T magnetic resonance imaging with diffusion-weighted imaging (DWI) using OGSE and PGSE (effective diffusion time [Δ_eff]: 6.5 ms and 35.2 ms) and b-values of 0 and 1000 s/mm². Using a region of interest (ROI) analysis of the brain tumors, we measured the ADC for two Δeff (ADC_Δeff) values and computed the subtraction ADC (ΔADC = ADC_{6.5 ms} - ADC_{35.2 ms}) and the relative ADC (ΔADC = (ADC_{6.5 ms} - ADC_{35.2 ms}) / ADC_{35.2 ms} × 100). The maximum values for the subtraction ADC (ΔADC_max) and the relative ADC (rADC_max) on the ROI were compared between low-grade and high-grade tumors using the Wilcoxon rank-sum test. A P-value <.05 was considered significant. The ROIs were also placed in the normal white matter of patients with high- and low-grade brain tumors, and ΔADC_max values were determined.

RESULTS

High-grade tumors had significantly higher ΔADC_max and rADC_max than low-grade tumors. The ΔADC_max values of the normal white matter were lower than the ΔADC_max of high- and low-grade brain tumors.

CONCLUSION

The dependence of ADC values on diffusion time between 6.5 ms and 35.2 ms was stronger in high-grade tumors than in low-grade tumors, suggesting differences in internal tissue structure. This finding highlights the importance of reporting diffusion times in ADC evaluations and might contribute to the grading of brain tumors using DWI.

Continuous Video Electroencephalogram during Hypoxia-Ischemia in Neonatal Mice

Hypoxia ischemia is the most common cause of neonatal seizures. Animal models are crucial for understanding the mechanisms and physiology underlying neonatal seizures and hypoxia ischemia. This manuscript describes a method for continuous video electroencephalogram (EEG) monitoring in neonatal mice to detect seizures and analyze EEG background during hypoxia ischemia. Use of video and EEG in conjunction allows description of seizure semiology and confirmation of seizures. This method also allows analysis of power spectrograms and EEG background pattern trends over the experimental time period. In this hypoxia ischemia model, the method allows EEG recording prior to injury to obtain a normative baseline and during injury and recovery. Total monitoring time is limited by the inability to separate pups from the mother for longer than four hours. Although, we have used a model of hypoxic-ischemic seizures in this manuscript, this method for neonatal video EEG monitoring could be applied to diverse disease and seizure models in rodents.

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.

Evidence of the diffusion time dependence of intravoxel incoherent motion in the brain

PURPOSE

To investigate the diffusion time (T_D ) dependence of intravoxel incoherent motion (IVIM) signals in the brain.

METHODS

A 3-compartment IVIM model was proposed to characterize 2 types of microcirculatory flows in addition to tissue water in the brain: flows that cross multiple vascular segments (pseudo-diffusive) and flows that stay in 1 segment (ballistic) within T_D . The model was first evaluated using simulated flow signals. Experimentally, flow-compensated (FC) pulsed-gradient spin-echo (PGSE) and oscillating-gradient spin-echo (OGSE) sequences were tested using a flow phantom and then used to examine IVIM signals in the mouse brain with T_D ranging from ~2.5 ms to 40 ms on an 11.7T scanner.

RESULTS

By fitting the model to simulated flow signals, we demonstrated the T_D dependency of the estimated fraction of pseudo-diffusive flow and the pseudo-diffusion coefficient (D*), which were dictated by the characteristic timescale of microcirculatory flow (τ). Flow phantom experiments validated that the OGSE and FC-PGSE sequences were not susceptible to the change in flow velocity. In vivo mouse brain data showed that both the estimated fraction of pseudo-diffusive flow and D* increased significantly as T_D increased.

CONCLUSION

We demonstrated that IVIM signals measured in the brain are T_D -dependent, potentially because more microcirculatory flows approach the pseudo-diffusive limit as T_D increases with respect to τ. Measuring the T_D dependency of IVIM signals may provide additional information on microvascular flows in the brain.

Oscillating-gradient diffusion magnetic resonance imaging detects acute subcellular structural changes in the mouse forebrain after neonatal hypoxia-ischemia

The recently developed oscillating-gradient diffusion MRI (OG-dMRI) technique extends our ability to examine brain structures at different spatial scales. In this study, we investigated the sensitivity of OG-dMRI in detecting cellular and subcellular structural changes in a mouse model of neonatal hypoxia ischemia (HI). Neonatal mice received unilateral HI injury or sham injury at postnatal day 10, followed by in vivo T₂-weighted and diffusion MRI of the brains at 3-6 h and 24 h after HI. Apparent diffusion coefficient (ADC) maps were acquired using conventional pulsed-gradient dMRI (PG-dMRI) and OG-dMRI with oscillating frequencies from 50 to 200 Hz. Pathology at cellular and subcellular levels was evaluated using neuronal, glial, and mitochondrial markers. We found significantly higher rates of ADC increase with oscillating frequencies (Δ_fADC) in the ipsilateral edema region, compared to the contralateral side, starting as early as 3 h after HI. Even in injured regions that showed no apparent change in PG-ADC or pseudo-normalized PG-ADC measurements, Δ_fADC remained significantly elevated. Histopathology showed swelling of sub-cellular structures in these regions with no apparent whole-cell level change. These results suggest that OG-dMRI is sensitive to subcellular structural changes in the brain after HI and is less susceptible to pseudo-normalization than PG-dMRI.

Reproducibility of Neurite Orientation Dispersion and Density Imaging (NODDI) in rats at 9.4 Tesla

PURPOSE

Neurite Orientation Dispersion and Density Imaging (NODDI) is a diffusion MRI (dMRI) technique used to characterize tissue microstructure by compartmental modelling of neural water fractions. Intra-neurite, extra-neurite, and cerebral spinal fluid volume fractions are measured. The purpose of this study was to determine the reproducibility of NODDI in the rat brain at 9.4 Tesla.

METHODS

Eight data sets were successfully acquired on adult male Sprague Dawley rats. Each rat was scanned twice on a 9.4T Agilent MRI with a 7 ± 1 day separation between scans. A multi-shell diffusion protocol was implemented consisting of 108 total directions varied over two shells (b-values of 1000 s/mm2 and 2000 s/mm2). Three techniques were used to analyze the NODDI scalar maps: mean region of interest (ROI) analysis, whole brain voxel-wise analysis, and targeted ROI analyses (voxel-wise within a given ROI). The coefficient of variation (CV) was used to assess the reproducibility of NODDI and provide insight into necessary sample sizes and minimum detectable effect size.

RESULTS

CV maps for orientation dispersion index (ODI) and neurite density index (NDI) showed high reproducibility both between and within subjects. Furthermore, it was found that small biological changes (<5%) may be detected with feasible sample sizes (n < 6-10). In contrast, isotropic volume fraction (IsoVF) was found to have low reproducibility, requiring very large sample sizes (n > 50) for biological changes to be detected.

CONCLUSIONS

The ODI and NDI measured by NODDI in the rat brain at 9.4T are highly reproducible and may be sensitive to subtle changes in tissue microstructure.

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.

Effects of nongaussian diffusion on "isotropic diffusion" measurements: An ex-vivo microimaging and simulation study

Designing novel diffusion-weighted pulse sequences to probe tissue microstructure beyond the conventional Stejskal-Tanner family is currently of broad interest. One such technique, multidimensional diffusion MRI, has been recently proposed to afford model-free decomposition of diffusion signal kurtosis into terms originating from either ensemble variance of isotropic diffusivity or microscopic diffusion anisotropy. This ability rests on the assumption that diffusion can be described as a sum of multiple Gaussian compartments, but this is often not strictly fulfilled. The effects of nongaussian diffusion on single shot isotropic diffusion sequences were first considered in detail by de Swiet and Mitra in 1996. They showed theoretically that anisotropic compartments lead to anisotropic time dependence of the diffusion tensors, which causes the measured isotropic diffusivity to depend on gradient frame orientation. Here we show how such deviations from the multiple Gaussian compartments assumption conflates orientation dispersion with ensemble variance in isotropic diffusivity. Second, we consider additional contributions to the apparent variance in isotropic diffusivity arising due to intracompartmental kurtosis. These will likewise depend on gradient frame orientation. We illustrate the potential importance of these confounds with analytical expressions, numerical simulations in simple model geometries, and microimaging experiments in fixed spinal cord using isotropic diffusion encoding waveforms with 7.5 ms duration and 3000 mT/m maximum amplitude.

Diffusion time dependence of microstructural parameters in fixed spinal cord

Biophysical modelling of diffusion MRI is necessary to provide specific microstructural tissue properties. However, estimating model parameters from data with limited diffusion gradient strength, such as clinical scanners, has proven unreliable due to a shallow optimization landscape. On the other hand, estimation of diffusion kurtosis (DKI) parameters is more robust, and its parameters may be connected to microstructural parameters, given an appropriate biophysical model. However, it was previously shown that this procedure still does not provide sufficient information to uniquely determine all model parameters. In particular, a parameter degeneracy related to the relative magnitude of intra-axonal and extra-axonal diffusivities remains. Here we develop a model of diffusion in white matter including axonal dispersion and demonstrate stable estimation of all model parameters from DKI in fixed pig spinal cord. By employing the recently developed fast axisymmetric DKI, we use stimulated echo acquisition mode to collect data over a two orders of magnitude diffusion time range with very narrow diffusion gradient pulses, enabling finely resolved measurements of diffusion time dependence of both net diffusion and kurtosis metrics, as well as model intra- and extra-axonal diffusivities, and axonal dispersion. Our results demonstrate substantial time dependence of all parameters except volume fractions, and the additional time dimension provides support for intra-axonal diffusivity to be larger than extra-axonal diffusivity in spinal cord white matter, although not unambiguously. We compare our findings for the time-dependent compartmental diffusivities to predictions from effective medium theory with reasonable agreement.

When tractography meets tracer injections: a systematic study of trends and variation sources of diffusion-based connectivity

Tractography is a powerful technique capable of non-invasively reconstructing the structural connections in the brain using diffusion MRI images, but the validation of tractograms is challenging due to lack of ground truth. Owing to recent developments in mapping the mouse brain connectome, high-resolution tracer injection-based axonal projection maps have been created and quickly adopted for the validation of tractography. Previous studies using tracer injections mainly focused on investigating the match in projections and optimal tractography protocols. Being a complicated technique, however, tractography relies on multiple stages of operations and parameters. These factors introduce large variabilities in tractograms, hindering the optimization of protocols and making the interpretation of results difficult. Based on this observation, in contrast to previous studies, in this work we focused on quantifying and ranking the amount of performance variation introduced by these factors. For this purpose, we performed over a million tractography experiments and studied the variability across different subjects, injections, anatomical constraints and tractography parameters. By using N-way ANOVA analysis, we show that all tractography parameters are significant and importantly performance variations with respect to the differences in subjects are comparable to the variations due to tractography parameters, which strongly underlines the importance of fully documenting the tractography protocols in scientific experiments. We also quantitatively show that inclusion of anatomical constraints is the most significant factor for improving tractography performance. Although this critical factor helps reduce false positives, our analysis indicates that anatomy-informed tractography still fails to capture a large portion of axonal projections.

Oscillating gradient diffusion kurtosis imaging of normal and injured mouse brains

Recent advances in diffusion MRI employ multiple diffusion encoding schemes with varying diffusion direction, weighting, and diffusion time to investigate specific microstructural properties in biological tissues. In this study, we examined time-dependent diffusion kurtosis contrast in adult mouse brains and in neonatal mouse brains after hypoxic-ischemic (HI) injury. In vivo diffusion kurtosis maps were acquired with a short diffusion time using an oscillating gradient spin echo (OGSE) sequence at 100 Hz and with a relatively long diffusion time (20 ms) using a pulsed gradient spin echo (PGSE) sequence. In the adult mouse brain, we found that the cortex and hippocampus showed larger differences between OGSE kurtosis and PGSE kurtosis than major white matter tracts. In neonatal mouse brains with unilateral HI injury, the OGSE kurtosis map overall provided stronger edema contrast than the PGSE kurtosis map, and the differences between OGSE and PGSE kurtosis measurements in the edema region reflected heterogeneity of injury. This is the first in vivo study that has demonstrated multi-direction OGSE kurtosis contrasts in the mouse brain. Comparing PGSE and OGSE kurtosis measures may provide additional information on microstructural changes after ischemic stroke.

Spatial Restriction within Intracranial Epidermoid Cysts Observed Using Short Diffusion-time Diffusion-weighted Imaging

We report two cases of pathologically proven intracranial epidermoid cysts. Both cases were scanned with diffusion-weighted imaging using pulsed gradient spin-echo (PGSE) and oscillating gradient spin-echo (OGSE; 50 Hz) prototype sequences with diffusion times of 47.3 ms and 8.5 ms, respectively. The apparent diffusion coefficient measured by OGSE was higher than that measured by PGSE, indicating the spatial restriction of water diffusion in the laminated keratin layers within the cyst as demonstrated by histopathology.

The Effect of Microcirculatory Flow on Oscillating Gradient Diffusion MRI and Diffusion Encoding with Dual-Frequency Orthogonal Gradients (DEFOG)

PURPOSE

We investigated the effect of microcirculatory flow on oscillating gradient spin echo (OGSE) diffusion MRI at low b-values and developed a diffusion preparation method called diffusion encoding with dual-frequency orthogonal gradients (DEFOG) to suppress the effect.

METHODS

Compared to conventional OGSE sequences, DEFOG adds a pulsed gradient that is orthogonal to the oscillating gradient and has a moderate diffusion weighting (e.g., 300 s/mm² ). In vivo MRI data were acquired from adult mouse brains (n = 5) on an 11.7 Tesla scanner, with diffusion times from 23.2 to 0.83 ms and b-values from 50 to 700 s/mm² .

RESULTS

Apparent diffusion coefficients (ADCs) measured using a conventional OGSE sequence at low b-values (< 200 mm² /s) were significantly higher than those measured at moderate b-values (> 300 mm² /s), potentially due to contributions from microcirculatory flow. In comparison, OGSE ADCs measured using the DEFOG method at low b-values were comparable to those measured at moderate b-values. The effect of microcirculatory flow on diffusion signals was diffusion time-dependent, and this dependency may reflect the capillary geometry and blood flow velocity in the mouse cortex.

CONCLUSION

Microcirculatory flow affects OGSE diffusion MRI measurements at low b-values, and this effect can be suppressed using the DEFOG method. Magn Reson Med 77:1583-1592, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Detection of acute nervous system injury with advanced diffusion-weighted MRI: a simulation and sensitivity analysis

Diffusion-weighted imaging (DWI) is a powerful tool to investigate the microscopic structure of the central nervous system (CNS). Diffusion tensor imaging (DTI), a common model of the DWI signal, has a demonstrated sensitivity to detect microscopic changes as a result of injury or disease. However, DTI and other similar models have inherent limitations that reduce their specificity for certain pathological features, particularly in tissues with complex fiber arrangements. Methods such as double pulsed field gradient (dPFG) and q-vector magic angle spinning (qMAS) have been proposed to specifically probe the underlying microscopic anisotropy without interference from the macroscopic tissue organization. This is particularly important for the study of acute injury, where abrupt changes in the microscopic morphology of axons and dendrites manifest as focal enlargements known as beading. The purpose of this work was to assess the relative sensitivity of DWI measures to beading in the context of macroscopic fiber organization and edema. Computational simulations of DWI experiments in normal and beaded axons demonstrated that, although DWI models can be highly specific for the simulated pathologies of beading and volume fraction changes in coherent fiber pathways, their sensitivity to a single idealized pathology is considerably reduced in crossing and dispersed fibers. However, dPFG and qMAS have a high sensitivity for beading, even in complex fiber tracts. Moreover, in tissues with coherent arrangements, such as the spinal cord or nerve fibers in which tract orientation is known a priori, a specific dPFG sequence variant decreases the effects of edema and improves specificity for beading. Collectively, the simulation results demonstrate that advanced DWI methods, particularly those which sample diffusion along multiple directions within a single acquisition, have improved sensitivity to acute axonal injury over conventional DTI metrics and hold promise for more informative clinical diagnostic use in CNS injury evaluation.

Reduction of Diffusion-Weighted Imaging Contrast of Acute Ischemic Stroke at Short Diffusion Times

Background and Purpose—

Diffusion-weighted imaging (DWI) of tissue water is a sensitive and specific indicator of acute brain ischemia, where reductions of the diffusion of tissue water are observed acutely in the stroke lesion core. Although these diffusion changes have been long attributed to cell swelling, the precise nature of the biophysical mechanisms remains uncertain.

Methods—

The potential cause of diffusion reductions after stroke was investigated using an advanced DWI technique, oscillating gradient spin-echo DWI, that enables much shorter diffusion times and can improve specificity for alterations of structure at the micron level.

Results—

Diffusion measurements in the white matter lesions of patients with acute ischemic stroke were reduced by only 8% using oscillating gradient spin-echo DWI, in contrast to a 37% decrease using standard DWI. Neurite beading has recently been proposed as a mechanism for the diffusion changes after ischemic stroke with some ex vivo evidence. To explore whether beading could cause such differential results, simulations of beaded cylinders and axonal swelling were performed, yielding good agreement with experiment.

Conclusions—

Short diffusion times result in dramatically reduced diffusion contrast of human stroke. Simulations implicate a combination of neuronal beading and axonal swelling as the key structural changes leading to the reduced apparent diffusion coefficient after stroke.