Intravoxel incoherent motion magnetic resonance imaging : basic principles and clinical applications

The purpose of this article was to show basic principles, acquisition, advantages, disadvantages, and clinical applications of intravoxel incoherent motion (IVIM) magnetic resonance imaging (MRI). IVIM MRI as a method was introduced in the late 1980s, but recently it started attracting more interest thanks to its applications in many fields, particularly in oncology and neuroradiology. This imaging technique has been developed with the objective of obtaining not only a functional analysis of different organs but also different types of lesions. Among many accessible tools in diagnostic imaging, IVIM MRI aroused the interest of many researchers in terms of studying its applicability in the evaluation of abdominal organs and diseases. The major conclusion of this article is that IVIM MRI seems to be a very auspicious method to investigate the human body, and that nowadays the most promising clinical application for IVIM perfusion MRI is oncology. However, due to lack of standardisation of image acquisition and analysis, further studies are needed to validate this method in clinical practice

Signal to Noise and b-value Analysis for Optimal Intra-Voxel Incoherent Motion Imaging in the Brain

Intravoxel incoherent motion (IVIM) is a method that can provide quantitative information about perfusion in the human body, in vivo, and without contrast agent. Unfortunately, the IVIM perfusion parameter maps are known to be relatively noisy in the brain, in particular for the pseudo-diffusion coefficient, which might hinder its potential broader use in clinical applications. Therefore, we studied the conditions to produce optimal IVIM perfusion images in the brain. IVIM imaging was performed on a 3-Tesla clinical system in four healthy volunteers, with 16 b values 0, 10, 20, 40, 80, 110, 140, 170, 200, 300, 400, 500, 600, 700, 800, 900 s/mm2, repeated 20 times. We analyzed the noise characteristics of the trace images as a function of b-value, and the homogeneity of the IVIM parameter maps across number of averages and sub-sets of the acquired b values. We found two peaks of noise of the trace images as function of b value, one due to thermal noise at high b-value, and one due to physiological noise at low b-value. The selection of b value distribution was found to have higher impact on the homogeneity of the IVIM parameter maps than the number of averages. Based on evaluations, we suggest an optimal b value acquisition scheme for a 12 min scan as 0 (7), 20 (4), 140 (19), 300 (9), 500 (19), 700 (1), 800 (4), 900 (1) s/mm2.Comment: 26 pages, 5 Figure

Perfusion MRI for brain tumors

Purpose : To compare data on brain tumors derived from intravoxel incoherent motion (IVIM) and arterial spin labeling (ASL) imaging with multiple parameters obtained on dynamic susceptibility contrast (DSC) perfusion MRI and to clarify the characteristics of IVIM and ASL perfusion data from the viewpoint of cerebral blood flow (CBF) analysis. Methods : ASL-CBF and IVIM techniques as well as DSC examination were performed in 24 patients with brain tumors. The IVIM data were analyzed with the two models. The relative blood flow (rBF), relative blood volume (rBV) corrected relative blood volume (crBV), mean transit time (MTT), and leakage coefficient (K2) were obtained from the DSC MRI data. Results : The ASL-CBF had the same tendency as the perfusion parameters derived from the DSC data, but the permeability from the vessels had less of an effect on the ASL-CBF. The diffusion coefficient of the fast component on IVIM contained more information on permeability than the f value. Conclusion : ASL-CBF is more suitable for the evaluation of perfusion in brain tumors than IVIM parameters. ASL-CBF and IVIM techniques should be carefully selected and the biological significance of each parameter should be understood for the correct comprehension of the pathological status of brain tumors

Rapid measurement of intravoxel incoherent motion (IVIM) derived perfusion fraction for clinical magnetic resonance imaging

Objective This study aimed to investigate the reliability of intravoxel incoherent motion (IVIM) model derived parameters D and f and their dependence on b value distributions with a rapid three b value acquisition protocol. Materials and methods Diffusion models for brain, kidney, and liver were assessed for bias, error, and reproducibility for the estimated IVIM parameters using b values 0 and 1000, and a b value between 200 and 900, at signal-to-noise ratios (SNR) 40, 55, and 80. Relative errors were used to estimate optimal b value distributions for each tissue scenario. Sixteen volunteers underwent brain DW-MRI, for which bias and coefficient of variation were determined in the grey matter. Results Bias had a large influence in the estimation of D and f for the low-perfused brain model, particularly at lower b values, with the same trends being confirmed by in vivo imaging. Significant differences were demonstrated in vivo for estimation of D (P = 0.029) and f (P < 0.001) with [300,1000] and [500,1000] distributions. The effect of bias was considerably lower for the high-perfused models. The optimal b value distributions were estimated to be brain500,1000, kidney300,1000, and liver200,1000. Conclusion IVIM parameters can be estimated using a rapid DW-MRI protocol, where the optimal b value distribution depends on tissue characteristics and compromise between bias and variability

Intravoxel Incoherent Motion Metrics as Potential Biomarkers for Survival in Glioblastoma.

Intravoxel incoherent motion (IVIM) is an MRI technique with potential applications in measuring brain tumor perfusion, but its clinical impact remains to be determined. We assessed the usefulness of IVIM-metrics in predicting survival in newly diagnosed glioblastoma. Fifteen patients with glioblastoma underwent MRI including spin-echo echo-planar DWI using 13 b-values ranging from 0 to 1000 s/mm2. Parametric maps for diffusion coefficient (D), pseudodiffusion coefficient (D*), and perfusion fraction (f) were generated for contrast-enhancing regions (CER) and non-enhancing regions (NCER). Regions of interest were manually drawn in regions of maximum f and on the corresponding dynamic susceptibility contrast images. Prognostic factors were evaluated by Kaplan-Meier survival and Cox proportional hazards analyses. We found that fCER and D*CER correlated with rCBFCER. The best cutoffs for 6-month survival were fCER&gt;9.86% and D*CER&gt;21.712 x10-3mm2/s (100% sensitivity, 71.4% specificity, 100% and 80% positive predictive values, and 80% and 100% negative predictive values; AUC:0.893 and 0.857, respectively). Treatment yielded the highest hazard ratio (5.484; 95% CI: 1.162-25.88; AUC: 0.723; P = 0.031); fCER combined with treatment predicted survival with 100% accuracy. The IVIM-metrics fCER and D*CER are promising biomarkers of 6-month survival in newly diagnosed glioblastoma

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

Correction for fast pseudo-diffusive fluid motion contaminations in diffusion tensor imaging

In this prospective study, we quantified the fast pseudo-diffusion contamination by blood perfusion or cerebrospinal fluid (CSF) intravoxel incoherent movements on the measurement of the diffusion tensor metrics in healthy brain tissue. Diffusion-weighted imaging (TR/TE = 4100 ms/90 ms; b-values: 0, 5, 10, 20, 35, 55, 80, 110, 150, 200, 300, 500, 750, 1000, 1300 s/mm2, 20 diffusion-encoding directions) was performed on a cohort of five healthy volunteers at 3 Tesla. The projections of the diffusion tensor along each diffusion-encoding direction were computed using a two b-value approach (2b), by fitting the signal to a monoexponential curve (mono), and by correcting for fast pseudo-diffusion compartments using the biexponential intravoxel incoherent motion model (IVIM) (bi). Fractional Anisotropy (FA) and Mean Diffusivity (MD) of the diffusion tensor were quantified in regions of interest drawn over white matter areas, gray matter areas, and the ventricles. A significant dependence of the MD from the evaluation method was found in all selected regions. A lower MD was computed when accounting for the fast-diffusion compartments. A larger dependence was found in the nucleus caudatus (bi: median 0.86 10-3 mm2/s, Δ2b: -11.2%, Δmono: -14.4%; p = 0.007), in the anterior horn (bi: median 2.04 10-3 mm2/s, Δ2b: -9.4%, Δmono: -11.5%, p = 0.007) and in the posterior horn of the lateral ventricles (bi: median 2.47 10-3 mm2/s, Δ2b: -5.5%, Δmono: -11.7%; p = 0.007). Also for the FA, the signal modeling affected the computation of the anisotropy metrics. The deviation depended on the evaluated region with significant differences mainly in the nucleus caudatus (bi: median 0.15, Δ2b: +39.3%, Δmono: +14.7%; p = 0.022) and putamen (bi: median 0.19, Δ2b: +3.1%, Δmono: +17.3%; p = 0.015). Fast pseudo-diffusive regimes locally affect diffusion tensor imaging (DTI) metrics in the brain. Here, we propose the use of an IVIM-based method for correction of signal contaminations through CSF or perfusion

Application of intravoxel incoherent motion (IVIM) model for differentiation between metastatic and non-metastatic head and neck lymph nodes

Background: Application of intravoxel incoherent motion (IVIM) model parameters, including: true diffusion (D), pseudodiffusion (D*), and perfusion fraction (Fp), for differentiation between metastatic and non-metastatic head and neck lymph nodes. Material/Methods: Diffusion-weighted images/apparent diffusion coefficient (DWI/ADC) images of 86 lymph nodes from 31 cancer patients were analyzed. DWI images were obtained with a 1.5T MRI scanner (Magnetom Avanto); b=0,50, 150, 300, 500, 750, 1000, 1200 s/mm2. Results: In the study group, there were 32 (37%) and 54 (67%) metastatic and non-metastatic lymph nodes, respectively. The mean values of D, D*, and Fp did not differ significantly between metastatic and non-metastatic lymph nodes. Conclusions: IVIM parameters are not useful for differentiation between metastatic and non-metastatic head and neck lymph nodes

Fitting IVIM with Variable Projection and Simplicial Optimization

Fitting multi-exponential models to Diffusion MRI (dMRI) data has always been challenging due to various underlying complexities. In this work, we introduce a novel and robust fitting framework for the standard two-compartment IVIM microstructural model. This framework provides a significant improvement over the existing methods and helps estimate the associated diffusion and perfusion parameters of IVIM in an automatic manner. As a part of this work we provide capabilities to switch between more advanced global optimization methods such as simplicial homology (SH) and differential evolution (DE). Our experiments show that the results obtained from this simultaneous fitting procedure disentangle the model parameters in a reduced subspace. The proposed framework extends the seminal work originated in the MIX framework, with improved procedures for multi-stage fitting. This framework has been made available as an open-source Python implementation and disseminated to the community through the DIPY project