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

Intravoxel incoherent motion as a surrogate marker of perfused vascular density in rat brain

Intravoxel incoherent motion (IVIM) MRI has emerged as a valuable technique for the assessment of tissue characteristics and perfusion. However, there is limited knowledge about the relationship between IVIM-derived measures and changes at the level of the vascular network. In this study, we investigated the potential use of IVIM MRI as a noninvasive tool for measuring changes in cerebral vascular density. Variations in quantitative immunohistochemical measurements of the vascular density across different regions in the rat brain (cortex, corpus callosum, hippocampus, thalamus, and hypothalamus) were related to the pseudo-diffusion coefficient D* and the flowing blood fraction f in healthy Wistar rats. We assessed whether region-wise differences in the vascular density are reflected by variations in the IVIM measurements and found a significant positive relationship with the pseudo-diffusion coefficient (p < 0.05, β = 0.24). The effect of cerebrovascular alterations, such as blood-brain barrier (BBB) disruption on the perfusion-related IVIM parameters, is not well understood. Therefore, we investigated the effect of BBB disruption on the IVIM measures in a rat model of metabolic and vascular comorbidities (ZSF1 obese rat) and assessed whether this affects the relationship between the cerebral vascular density and the noninvasive IVIM measurements. We observed increased vascular permeability without detecting any differences in diffusivity, suggesting that BBB leakage is present before changes in the tissue integrity. We observed no significant difference in the relationship between cerebral vascular density and the IVIM measurements in our model of comorbidities compared with healthy normotensive rats.

Reproducibility of intravoxel incoherent motion quantification in the liver across field strengths and gradient hardware

PURPOSE

To evaluate reproducibility and interlobar agreement of intravoxel incoherent motion (IVIM) quantification in the liver across field strengths and MR scanners with different gradient hardware.

METHODS

Cramer-Rao lower bound optimization was performed to determine optimized monopolar and motion-robust 2D (b-value and first-order motion moment [M1]) IVIM-DWI acquisitions. Eleven healthy volunteers underwent diffusion MRI of the liver, where each optimized acquisition was obtained five times across three MRI scanners. For each data set, IVIM estimates (diffusion coefficient (D), pseudo-diffusion coefficients (d 1 * $$ {d}_1^{\ast } $$ andd 2 * $$ {d}_2^{\ast } $$ ), blood velocity SDs (Vb1 and Vb2), and perfusion fractions [f1 and f2]) were obtained in the right and left liver lobes using two signal models (pseudo-diffusion and M1-dependent physical) with and without T2 correction (fc1 and fc2) and three fitting techniques (tri-exponential region of interest-based full and segmented fitting and blood velocity SD distribution fitting). Reproducibility and interlobar agreement were compared across methods using within-subject and pairwise coefficients of variation (CVw and CVp), paired sample t-tests, and Bland-Altman analysis.

RESULTS

Using a combination of motion-robust 2D (b-M1) data acquisition, M1-dependent physical signal modeling with T2 correction, and blood velocity SD distribution fitting, multiscanner reproducibility with median CVw = 5.09%, 11.3%, 9.20%, 14.2%, and 12.6% for D, Vb1, Vb2, fc1, and fc2, respectively, and interlobar agreement with CVp = 8.14%, 11.9%, 8.50%, 49.9%, and 42.0%, respectively, was achieved.

CONCLUSION

Recently proposed advanced IVIM acquisition, signal modeling, and fitting techniques may facilitate reproducible IVIM quantification in the liver, as needed for establishment of IVIM-based quantitative biomarkers for detection, staging, and treatment monitoring of diseases.

Comparison of methods for intravoxel incoherent motion parameter estimation in the brain from flow-compensated and non-flow-compensated diffusion-encoded data

PURPOSE

Joint analysis of flow-compensated (FC) and non-flow-compensated (NC) diffusion MRI (dMRI) data has been suggested for increased robustness of intravoxel incoherent motion (IVIM) parameter estimation. For this purpose, a set of methods commonly used or previously found useful for IVIM analysis of dMRI data obtained with conventional diffusion encoding were evaluated in healthy human brain.

METHODS

Five methods for joint IVIM analysis of FC and NC dMRI data were compared: (1) direct non-linear least squares fitting, (2) a segmented fitting algorithm with estimation of the diffusion coefficient from higher b-values of NC data, (3) a Bayesian algorithm with uniform prior distributions, (4) a Bayesian algorithm with spatial prior distributions, and (5) a deep learning-based algorithm. Methods were evaluated on brain dMRI data from healthy subjects and simulated data at multiple noise levels. Bipolar diffusion encoding gradients were used with b-values 0-200 s/mm2 and corresponding flow weighting factors 0-2.35 s/mm for NC data and by design 0 for FC data. Data were acquired twice for repeatability analysis.

RESULTS

Measurement repeatability as well as estimation bias and variability were at similar levels or better with the Bayesian algorithm with spatial prior distributions and the deep learning-based algorithm for IVIM parametersD $$ D $$ andf $$ f $$ , and for the Bayesian algorithm only forv d $$ {v}_d $$ , relative to the other methods.

CONCLUSION

A Bayesian algorithm with spatial prior distributions is preferable for joint IVIM analysis of FC and NC dMRI data in the healthy human brain, but deep learning-based algorithms appear promising.

In Vivo and Post-mortem Comparisons of IVIM/Time-dependent Diffusion MR Imaging Parameters in Melanoma and Breast Cancer Xenograft Models

PURPOSE

We aimed to investigate the changes in intravoxel incoherent motion (IVIM) and diffusion parameters between in vivo and post-mortem conditions and the time dependency of these parameters using two different mouse tumor models with different vessel lumen sizes.

METHODS

Six B16 and six MDA-MB-231 xenograft mice were scanned using 7 Tesla MRI under both in vivo/post-mortem conditions. Diffusion weighted imaging with 17 b-values (0-3000 s/mm2) were obtained at two diffusion times (9 and 27.6 ms). The shifted apparent diffusion coefficient (sADC) using 2 b-values (200 and 1500 s/mm2), non-Gaussian diffusion and IVIM parameters (ADC0, K, fIVIM) were estimated at each of the diffusion times. The results were evaluated by repeated measures two-way analysis of variance and post hoc Bonferroni test.

RESULTS

In B16 tumors, fIVIM significantly decreased with post-mortem conditions (from 12.6 ± 6.5% to 5.2 ± 1.9%, P < 0.05 at long diffusion time; from 11.0 ± 2.4% to 4.6 ± 2.7%, P < 0.05 at short diffusion time). In MDA-MB-231 tumors, fIVIM also significantly decreased (from 8.8 ± 3.8% to 2.6 ± 1.1%, P < 0.05 at long; from 7.9 ± 5.4% to 2.9 ± 1.1%, P < 0.05 at short). No diffusion time dependency was observed (P = 0.59 in B16 and P = 0.77 in MDA-MB-231). The sADC and ADC0 values tended to decrease and the K value tended to increase after sacrificing and when increasing the diffusion time.

CONCLUSION

The fIVIM values dropped after sacrificing, confirming that IVIM MRI is a promising quantitative parameter to evaluate blood microcirculation. The presence of residual post-mortem fIVIM values suggested that the influence of water molecule diffusion in the blood lumen may contribute to the IVIM effect. Diffusion MRI parameter's time dependency and those changes after sacrificing could possibly provide additional insights into diffusion hindrance mechanisms.

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.

Characterization of breast lesions using multi-parametric diffusion MRI and machine learning

Objective. To investigate quantitative imaging markers based on parameters from two diffusion-weighted imaging (DWI) models, continuous-time random-walk (CTRW) and intravoxel incoherent motion (IVIM) models, for characterizing malignant and benign breast lesions by using a machine learning algorithm.Approach. With IRB approval, 40 women with histologically confirmed breast lesions (16 benign, 24 malignant) underwent DWI with 11b-values (50 to 3000 s/mm2) at 3T. Three CTRW parameters,Dm,α, andβand three IVIM parametersDdiff,Dperf, andfwere estimated from the lesions. A histogram was generated and histogram features of skewness, variance, mean, median, interquartile range; and the value of the 10%, 25% and 75% quantiles were extracted for each parameter from the regions-of-interest. Iterative feature selection was performed using the Boruta algorithm that uses the Benjamin Hochberg False Discover Rate to first determine significant features and then to apply the Bonferroni correction to further control for false positives across multiple comparisons during the iterative procedure. Predictive performance of the significant features was evaluated using Support Vector Machine, Random Forest, Naïve Bayes, Gradient Boosted Classifier (GB), Decision Trees, AdaBoost and Gaussian Process machine learning classifiers.Main Results. The 75% quantile, and median ofDm; 75% quantile off;mean, median, and skewness ofβ;kurtosis ofDperf; and 75% quantile ofDdiffwere the most significant features. The GB differentiated malignant and benign lesions with an accuracy of 0.833, an area-under-the-curve of 0.942, and an F1 score of 0.87 providing the best statistical performance (p-value < 0.05) compared to the other classifiers.Significance. Our study has demonstrated that GB with a set of histogram features from the CTRW and IVIM model parameters can effectively differentiate malignant and benign breast lesions.

Loss of aquaporin-4 results in glymphatic system dysfunction via brain-wide interstitial fluid stagnation

The glymphatic system is a fluid transport network of cerebrospinal fluid (CSF) entering the brain along arterial perivascular spaces, exchanging with interstitial fluid (ISF), ultimately establishing directional clearance of interstitial solutes. CSF transport is facilitated by the expression of aquaporin-4 (AQP4) water channels on the perivascular endfeet of astrocytes. Mice with genetic deletion of AQP4 (AQP4 KO) exhibit abnormalities in the brain structure and molecular water transport. Yet, no studies have systematically examined how these abnormalities in structure and water transport correlate with glymphatic function. Here, we used high-resolution 3D magnetic resonance (MR) non-contrast cisternography, diffusion-weighted MR imaging (MR-DWI) along with intravoxel-incoherent motion (IVIM) DWI, while evaluating glymphatic function using a standard dynamic contrast-enhanced MR imaging to better understand how water transport and glymphatic function is disrupted after genetic deletion of AQP4. AQP4 KO mice had larger interstitial spaces and total brain volumes resulting in higher water content and reduced CSF space volumes, despite similar CSF production rates and vascular density compared to wildtype mice. The larger interstitial fluid volume likely resulted in increased slow but not fast MR diffusion measures and coincided with reduced glymphatic influx. This markedly altered brain fluid transport in AQP4 KO mice may result from a reduction in glymphatic clearance, leading to enlargement and stagnation of fluid in the interstitial space. Overall, diffusion MR is a useful tool to evaluate glymphatic function and may serve as valuable translational biomarker to study glymphatics in human disease.

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.

Intravoxel Incoherent Motion Magnetic Resonance Imaging in Skeletal Muscle: Review and Future Directions

Throughout the body, muscle structure and function can be interrogated using a variety of noninvasive magnetic resonance imaging (MRI) methods. Recently, intravoxel incoherent motion (IVIM) MRI has gained momentum as a method to evaluate components of blood flow and tissue diffusion simultaneously. Much of the prior research has focused on highly vascularized organs, including the brain, kidney, and liver. Unique aspects of skeletal muscle, including the relatively low perfusion at rest and large dynamic range of perfusion between resting and maximal hyperemic states, may influence the acquisition, postprocessing, and interpretation of IVIM data. Here, we introduce several of those unique features of skeletal muscle; review existing studies of IVIM in skeletal muscle at rest, in response to exercise, and in disease states; and consider possible confounds that should be addressed for muscle-specific evaluations. Most studies used segmented nonlinear least squares fitting with a b-value threshold of 200 sec/mm² to obtain IVIM parameters of perfusion fraction (f), pseudo-diffusion coefficient (D*), and diffusion coefficient (D). In healthy individuals, across all muscles, the average ± standard deviation of D was 1.46 ± 0.30 × 10^-3  mm² /sec, D* was 29.7 ± 38.1 × 10^-3  mm² /sec, and f was 11.1 ± 6.7%. Comparisons of reported IVIM parameters in muscles of the back, thigh, and leg of healthy individuals showed no significant difference between anatomic locations. Throughout the body, exercise elicited a positive change of all IVIM parameters. Future directions including advanced postprocessing models and potential sequence modifications are discussed. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY: Stage 2.

Characterisation of microvessel blood velocity and segment length in the brain using multi-diffusion-time diffusion-weighted MRI

Multi-diffusion-time diffusion-weighted MRI can probe tissue microstructure, but the method has not been widely applied to the microvasculature. At long diffusion-times, blood flow in capillaries is in the diffusive regime, and signal attenuation is dependent on blood velocity (v) and capillary segment length (l). It is described by the pseudo-diffusion coefficient (D*=vl/6) of intravoxel incoherent motion (IVIM). At shorter diffusion-times, blood flow is in the ballistic regime, and signal attenuation depends on v, and not l. In theory, l could be estimated using D* and v. In this study, we compare the accuracy and repeatability of three approaches to estimating v, and therefore l: the IVIM ballistic model, the velocity autocorrelation model, and the ballistic approximation to the velocity autocorrelation model. Twenty-nine rat datasets from two strains were acquired at 7 T, with b-values between 0 and 1000 smm^-2 and diffusion times between 11.6 and 50 ms. Five rats were scanned twice to assess scan-rescan repeatability. Measurements of l were validated using corrosion casting and micro-CT imaging. The ballistic approximation of the velocity autocorrelation model had lowest bias relative to corrosion cast estimates of l, and had highest repeatability.

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.

MR cell size imaging with temporal diffusion spectroscopy

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

Improved gradient waveforms for oscillating gradient spin-echo (OGSE) diffusion tensor imaging

The dependence of the diffusion tensor on frequency is of great interest in studies of tissue microstructure because it reveals restrictions to the Brownian motion of water molecules caused by cell membranes. Oscillating gradient spin-echo (OGSE) sequences can sample this dependence with gradient shapes for which the power spectrum of the zeroth moment is focused at a target frequency. In order to maintain the total spectral power (ie the b-value), oscillating gradient amplitudes must grow with the frequency squared. For this reason, OGSE applications on clinical MRI scanners are limited to low frequencies, for which it is difficult to obtain a narrow frequency bandwidth of the gradient moment in a useful echo time. In particular, the commonly used pair of single-period trapezoidal-cosine pulses separated by a half-period produces significant side lobes away from the target frequency. To mitigate this effect, improved OGSE waveforms are proposed, which reduce the gap between the two gradient pulses to the minimum duration required for the refocusing RF pulse. Additionally, a slight deviation from the periodicity of the waveforms is proposed in order to permit using the maximum slew rate of the gradient system for all lobes of trapezoidal waveforms while maintaining advantageous spectral properties, which is not the case for the currently used OGSE sequences. Numerical calculations validate these changes, showing that both modifications significantly narrow the gradient moment power spectrum and increase the contribution of its main lobe to the b-value, thus improving the specificity of the measurement. The utility of the new shapes is demonstrated by diffusion tensor measurements of human white matter in vivo over the range of 30-75 Hz with a b-value of nearly 1000 s/mm² , using a high-performance gradient insert. However, the improvement should increase the sampling precision of OGSE experiments for all gradient systems.

Probing the ballistic microcirculation in placenta using flow-compensated and non-compensated intravoxel incoherent motion imaging

PURPOSE

Intravoxel incoherent motion (IVIM) imaging is widely used to evaluate microcirculatory flow, which consists of diffusive and ballistic flow components. We proposed a joint use of flow-compensated (FC) and non-compensated (NC) diffusion gradients to probe the fraction and velocity of ballistic flow in the placenta.

METHODS

Forty pregnant women were included in this study and scanned on a 1.5T clinical scanner. FC and NC diffusion MRI (dMRI) sequences were achieved using a pair of identical or mirrored bipolar gradients. A joint FC-NC model was established to estimate the fraction (f_b ) and velocity (v_b ) of the ballistic flow. Conventional IVIM parameters (f, D, and D*) were obtained from the FC and NC data, separately. The v_b and f·D*, as placental flow velocity measurements, were correlated with the umbilical-artery Doppler ultrasound indices and gestational ages.

RESULTS

The ballistic flow component can be observed from the difference between the FC and NC dMRI signal decay curves. v_b fitted from the FC-NC model showed strong correlations with umbilical-artery impedance indices, the systolic-to-diastolic (SD) ratio and pulsatility index (PI), with correlation coefficients of 0.65 and 0.62. The f·D* estimated from the NC data positively correlated with SD and PI, while the FC-based f·D* values showed weak negative correlations. Significant gestational-age dependence was also found in the flow velocity measurements.

CONCLUSION

Our results demonstrated the feasibility of using FC and NC dMRI to noninvasively measure ballistic flow velocity in the placenta, which may be used as a new marker to evaluate placenta microcirculation.