Intravoxel incoherent motion


Intravoxel incoherent motion imaging is a concept and a method initially introduced and developed by Le Bihan et al. to quantitatively assess all the microscopic translational motions that could contribute to the signal acquired with diffusion MRI. In this model, biological tissue contains two distinct environments: molecular diffusion of water in the tissue, and microcirculation of blood in the capillary network. The concept introduced by D. Le Bihan is that water flowing in capillaries mimics a random walk , as long as the assumption that all directions are represented in the capillaries is satisfied.
It is responsible for a signal attenuation in diffusion MRI, which depends on the velocity of the flowing blood and the vascular architecture. Similarly to molecular diffusion, the effect of pseudodiffusion on the signal attenuation depends on the b value. However, the rate of signal attenuation resulting from pseudodiffusion is typically an order of magnitude greater than molecular diffusion in tissues, so its relative contribution to the diffusion-weighted MRI signal becomes significant only at very low b values, allowing diffusion and perfusion effects to be separated.

Model

In the presence of the magnetic field gradient pulses of a diffusion MRI sequence, the MRI signal gets attenuated due to diffusion and perfusion effects. In a simple model, this signal attenuation, S/So, can be written as:
where is the volume fraction of incoherently flowing blood in the tissue, the signal attenuation from the IVIM effect and is the signal attenuation from molecular diffusion in the tissue.
Assuming blood water flowing in the randomly oriented vasculature changes several times direction during the measurement time, one has for :
where is the diffusion-sensitization of the MRI sequence, is the sum of the pseudo-diffusion coefficient associated to the IVIM effect and, the diffusion coefficient of water in blood:
where is the mean capillary segment length and is the blood velocity.
If blood water flows without changing direction while capillary segments are randomly and isotropically oriented, becomes:
where is a parameter linked to the gradient pulse amplitude and time course.
In both cases, the perfusion effect results in a curvature of the diffusion attenuation plot towards b=0. In a simple approach and under some approximations, the ADC calculated from 2 diffusion-weighted images acquired with b0=0 and b1, as ADC = ln/S ), is:
where is the tissue diffusion coefficient. The ADC thus only depends on the flowing vascular volume and not on the blood velocity and capillary geometry, which is a strong advantage. The contribution of perfusion to the ADC is larger when using small b values.
On the other hand, set of data obtained from images acquired with a multiple b values can be fitted with Eq. using either model 1 or model 2 to estimate and/or blood velocity.
The late part of the curve also presents some degree of curvature. This is because diffusion in biological tissues is not free, but can be hindered by many obstacles or even restricted. Several models have been proposed to describe this curvature at higher b-values, mainly the “biexponential” model which assumes the presence of 2 water compartments with fast and slow diffusion , the relative 'fast' and 'slow' labels referring to restricted and hindered diffusion, rather than pseudodiffusion/perfusion and true diffusion. Another alternative is the “kurtosis” model which quantifies the deviation from free diffusion in the parameter .
Biexponential model:
Where and are the relative fractions and diffusion coefficients of the fast and slow compartments. This general formulation of a biexponential decay of diffusion-weighted imaging signal with b-value can be used for IVIM, which requires sampling of low b-values to capture pseudodiffusion decay, or for restriction imaging, which requires higher b-value acquisitions to capture restricted diffusion.
Kurtosis model:
where is the tissue intrinsic diffusion coefficient and the Kurtosis parameter.
Both models can be related assuming some hypotheses about the tissue structure and the measurement conditions.
Separation of perfusion from diffusion requires good signal-to-noise ratios and there are some technical challenges to overcome. Also the “perfusion” parameters accessible with the IVIM method somewhat differs from the “classical” perfusion parameters obtained with tracer methods: “Perfusion” can be seen with the physiologist eyes or the radiologist eyes. Indeed, there is room to improve the IVIM model and understand better its relationship with the functional vascular architecture and its biological relevance.

Applications

IVIM MRI was initially introduced to evaluate perfusion and produce maps of brain perfusion, for brain activation studies and clinical applications. Recent work has proven the validity of the IVIM concept from fMRI, with an increase in the IVIM perfusion parameters in brain activated regions, and the potential of the approach to aid in our understanding of the different vascular contributions to the fMRI signal. IVIM MRI has also been used in the context of fMRI in a negative way.
A limitation of BOLD fMRI is its spatial resolution, as flow increase in somewhat large arteries or veins feed or drain large neuronal territories. By inserting “diffusion” gradient pulses in the MRI sequence, one may crush the contribution of the largest vessels in the BOLD signal and improve the spatial resolution of the activation maps. Several groups have relied on this trick, though not always considering referring to the IVIM concept. This IVIM concept has also been borrowed to improve other applications, for instance, arterial spin labeling or to suppress signal from extracellular flowing fluid in perfused cell systems.
However, IVIM MRI has recently undergone a striking revival for applications not in the brain, but throughout the body as well. Following earlier encouraging results in the kidneys, or even the heart, IVIM MRI really took off for liver applications. For instance, Luciani et al. found that D* was significantly reduced in cirrhotic patients, which, according to the IVIM model, points out to reduce blood velocity.. The perfusion fraction, f, which is linked to blood volume in the IVIM model, remained normal, confirming earlier results by Yamada et al. Though, blood volume is expected to be reduced in liver cirrhosis.
One has to keep in mind that IVIM imaging has a differential sensitivity to vessel types, according to the range of motion sensitization which are used. Signal from large vessels with rapid flow disappears quickly with very low b values, while smaller vessels with slower flow might still contribute to the IVIM signal acquired with b values larger than 200 s/mm². It has also been shown that the parameter f, often related to perfusion fraction, is sensitive to differential spin-spin relaxation rates in the two model compartments and can thus be overestimated in highly-perfused tissue. Correction of this effect is achieved by additional images at a different echo time.
Many more applications are now under investigation, especially for imaging of patients suspected of cancer in the body and human placenta. A key feature of IVIM diffusion MRI is that it does not involve contrast agents, and it may appear as an interesting alternative for perfusion MRI in some patients at risk for Nephrogenic Systemic Fibrosis.