Formation of parametric images using mixed-effects models: a feasibility study

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.

A robust deconvolution method to disentangle multiple water pools in diffusion MRI

The diffusion-weighted magnetic resonance imaging (dMRI) signal measured in vivo arises from multiple diffusion domains, including hindered and restricted water pools, free water and blood pseudo-diffusion. Not accounting for the correct number of components can bias metrics obtained from model fitting because of partial volume effects that are present in, for instance, diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI). Approaches that aim to overcome this shortcoming generally make assumptions about the number of considered components, which are not likely to hold for all voxels. The spectral analysis of the dMRI signal has been proposed to relax assumptions on the number of components. However, it currently requires a clinically challenging signal-to-noise ratio (SNR) and accounts only for two diffusion processes defined by hard thresholds. In this work, we developed a method to automatically identify the number of components in the spectral analysis, and enforced its robustness to noise, including outlier rejection and a data-driven regularization term. Furthermore, we showed how this method can be used to take into account partial volume effects in DTI and DKI fitting. The proof of concept and performance of the method were evaluated through numerical simulations and in vivo MRI data acquired at 3 T. With simulations our method reliably decomposed three diffusion components from SNR = 30. Biases in metrics derived from DTI and DKI were considerably reduced when components beyond hindered diffusion were taken into account. With the in vivo data our method determined three macro-compartments, which were consistent with hindered diffusion, free water and pseudo-diffusion. Taking free water and pseudo-diffusion into account in DKI resulted in lower mean diffusivity and higher fractional anisotropy values in both gray and white matter. In conclusion, the proposed method allows one to determine co-existing diffusion compartments without prior assumptions on their number, and to account for undesired signal contaminations within clinically achievable SNR levels.

Bayesian intravoxel incoherent motion parameter mapping in the human heart

BACKGROUND

Intravoxel incoherent motion (IVIM) imaging of diffusion and perfusion in the heart suffers from high parameter estimation error. The purpose of this work is to improve cardiac IVIM parameter mapping using Bayesian inference.

METHODS

A second-order motion-compensated diffusion weighted spin-echo sequence with navigator-based slice tracking was implemented to collect cardiac IVIM data in early systole in eight healthy subjects on a clinical 1.5 T CMR system. IVIM data were encoded along six gradient optimized directions with b-values of 0-300 s/mm2. Subjects were scanned twice in two scan sessions one week apart to assess intra-subject reproducibility. Bayesian shrinkage prior (BSP) inference was implemented to determine IVIM parameters (diffusion D, perfusion fraction F and pseudo-diffusion D*). Results were compared to least-squares (LSQ) parameter estimation. Signal-to-noise ratio (SNR) requirements for a given fitting error were assessed for the two methods using simulated data. Reproducibility analysis of parameter estimation in-vivo using BSP and LSQ was performed.

RESULTS

BSP resulted in reduced SNR requirements when compared to LSQ in simulations. In-vivo, BSP analysis yielded IVIM parameter maps with smaller intra-myocardial variability and higher estimation certainty relative to LSQ. Mean IVIM parameter estimates in eight healthy subjects were (LSQ/BSP): 1.63 ± 0.28/1.51 ± 0.14·10-3 mm2/s for D, 13.13 ± 19.81/13.11 ± 5.95% for F and 201.45 ± 313.23/13.11 ± 14.53·10-3 mm2/s for D ∗. Parameter variation across all volunteers and measurements was lower with BSP compared to LSQ (coefficient of variation BSP vs. LSQ: 9% vs. 17% for D, 45% vs. 151% for F and 111% vs. 155% for D ∗). In addition, reproducibility of the IVIM parameter estimates was higher with BSP compared to LSQ (Bland-Altman coefficients of repeatability BSP vs. LSQ: 0.21 vs. 0.26·10-3 mm2/s for D, 5.55 vs. 6.91% for F and 15.06 vs. 422.80·10-3 mm2/s for D*).

CONCLUSION

Robust free-breathing cardiac IVIM data acquisition in early systole is possible with the proposed method. BSP analysis yields improved IVIM parameter maps relative to conventional LSQ fitting with fewer outliers, improved estimation certainty and higher reproducibility. IVIM parameter mapping holds promise for myocardial perfusion measurements without the need for contrast agents.