Generalization of intravoxel incoherent motion model by introducing the notion of continuous pseudodiffusion variable

Probing Renal Microstructure and Function with Advanced Diffusion MRI: Concepts, Applications, Challenges, and Future Directions

Diffusion measurements in the kidney are affected not only by renal microstructure but also by physiological processes (i.e., glomerular filtration, water reabsorption, and urine formation). Because of the superposition of passive tissue diffusion, blood perfusion, and tubular pre-urine flow, the limitations of the monoexponential apparent diffusion coefficient (ADC) model in assessing pathophysiological changes in renal tissue are becoming apparent and motivate the development of more advanced diffusion-weighted imaging (DWI) variants. These approaches take advantage of the fact that the length scale probed in DWI measurements can be adjusted by experimental parameters, including diffusion-weighting, diffusion gradient directions and diffusion time. This forms the basis by which advanced DWI models can be used to capture not only passive diffusion effects, but also microcirculation, compartmentalization, tissue anisotropy. In this review, we provide a comprehensive overview of the recent advancements in the field of renal DWI. Following a short introduction on renal structure and physiology, we present the key methodological approaches for the acquisition and analysis of renal DWI data, including intravoxel incoherent motion (IVIM), diffusion tensor imaging (DTI), non-Gaussian diffusion, and hybrid IVIM-DTI. We then briefly summarize the applications of these methods in chronic kidney disease and renal allograft dysfunction. Finally, we discuss the challenges and potential avenues for further development of renal DWI. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 2.

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/mm2 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  mm2 /sec, D* was 29.7 ± 38.1 × 10-3  mm2 /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.

Short-Term Repeatability of in Vivo Cardiac Intravoxel Incoherent Motion Tensor Imaging in Healthy Human Volunteers

BACKGROUND

Intravoxel incoherent motion (IVIM) tensor imaging is a promising technique for diagnosis and monitoring of cardiovascular diseases. Knowledge about measurement repeatability, however, remains limited.

PURPOSE

To evaluate short-term repeatability of IVIM tensor imaging in normal in vivo human hearts.

STUDY TYPE

Prospective.

POPULATION

Ten healthy subjects without history of heart diseases.

FIELD STRENGTH/SEQUENCE

Balanced steady-state free-precession cine sequence and single-shot spin-echo echo planar IVIM tensor imaging sequence (9 b-values, 0-400 seconds/mm² and six diffusion-encoding directions) at 3.0 T.

ASSESSMENT

Subjects were scanned twice with an interval of 15 minutes, leaving the scanner between studies. The signal-to-noise ratio (SNR) was evaluated in anterior, lateral, septal, and inferior segments of the left ventricle wall. Fractional anisotropy (FA), mean diffusivity (MD), mean fraction (MF), and helix angle (HA) in the four segments were independently measured by five radiologists.

STATISTICAL TESTS

IVIM tensor indexes were compared between observers using a one-way analysis of variance or between scans using a paired t-test (normal data) or a Wilcoxon rank-sum test (non-normal data). Interobserver agreement and test-retest repeatability were assessed using the intraclass correlation coefficient (ICC), within-subject coefficient of variation (WCV), and Bland-Altman limits of agreements.

RESULTS

SNR of inferior segment was significantly lower than the other three segments, and inferior segment was therefore excluded from repeatability analysis. Interobserver repeatability was excellent for all IVIM tensor indexes (ICC: 0.886-0.972; WCV: 0.62%-4.22%). Test-retest repeatability was excellent for MD of the self-diffusion tensor (D) and MF of the perfusion fraction tensor (f_p ) (ICC: 0.803-0.888; WCV: 1.42%-9.51%) and moderate for FA and MD of the pseudo-diffusion tensor (D^* ) (ICC: 0.487-0.532; WCV: 6.98%-10.89%). FA of D and f_p and HA of D presented good test-retest repeatability (ICC: 0.732-0.788; WCV: 3.28%-8.71%).

DATA CONCLUSION

The D and f_p indexes exhibited satisfactory repeatability, but further efforts were needed to improve repeatability of D^* indexes.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 1.

Evidence of Tri-Exponential Decay for Liver Intravoxel Incoherent Motion MRI: A Review of Published Results and Limitations

Diffusion weighted imaging (DWI) and intravoxel incoherent motion (IVIM) have been explored to assess liver tumors and diffused liver diseases. IVIM reflects the microscopic translational motions that occur in voxels in magnetic resonance (MR) DWI. In biologic tissues, molecular diffusion of water and microcirculation of blood in the capillary network can be assessed using IVIM DWI. The most commonly applied model to describe the DWI signal is a bi-exponential model, with a slow compartment of diffusion linked to pure molecular diffusion (represented by the coefficient D_slow), and a fast compartment of diffusion, related to microperfusion (represented by the coefficient D_fast). However, high variance in D_fast estimates has been consistently shown in literature for liver IVIM, restricting its application in clinical practice. This variation could be explained by the presence of another very fast compartment of diffusion in the liver. Therefore, a tri-exponential model would be more suitable to describe the DWI signal. This article reviews the published evidence of the existence of this additional very fast diffusion compartment and discusses the performance and limitations of the tri-exponential model for liver IVIM in current clinical settings.

Investigation of intravoxel incoherent motion tensor imaging for the characterization of the in vivo human heart

PURPOSE

To investigate intravoxel incoherent motion (IVIM) tensor imaging of the in vivo human heart and elucidate whether the estimation of IVIM tensors is affected by the complexity of pseudo-diffusion components in myocardium.

METHODS

The cardiac IVIM data of 10 healthy subjects were acquired using a diffusion weighted spin-echo echo-planar imaging sequence along 6 gradient directions with 10 b values (0~400 s/mm² ). The IVIM data of left ventricle myocardium were fitted to the IVIM tensor model. The complexity of myocardial pseudo-diffusion components was reduced through exclusion of low b values (0 and 5 s/mm² ) from the IVIM curve-fitting analysis. The fractional anisotropy, mean fraction/mean diffusivity, and Westin measurements of pseudo-diffusion tensors (f_p and D*) and self-diffusion tensor (D), as well as the angle between the main eigenvector of f_p (or D*) and that of D, were computed and compared before and after excluding low b values.

RESULTS

The fractional anisotropy values of f_p and D* without low b value participation were significantly higher (P < .001) than those with low b value participation, but an opposite trend was found for the mean fraction/diffusivity values. Besides, after removing low b values, the angle between the main eigenvector of f_p (or D*) and that of D became small, and both f_p and D* tensors presented significant decrease of spherical components and significant increase of linear components.

CONCLUSION

The presence of multiple pseudo-diffusion components in myocardium indeed influences the estimation of IVIM tensors. The IVIM tensor model needs to be further improved to account for the complexity of myocardial microcirculatory network and blood flow.

Accurate intravoxel incoherent motion parameter estimation using Bayesian fitting and reduced number of low b-values

PURPOSE

Intravoxel incoherent motion (IVIM) magnetic resonance imaging is a potential noninvasive technique for the diagnosis of brain tumors. However, perfusion-related parameter mapping is a persistent problem. The purpose of this paper is to investigate the IVIM parameter mapping of brain tumors using Bayesian fitting and low b-values.

METHODS

Bayesian shrinkage prior (BSP) fitting method and different low b-value distributions were used to estimate IVIM parameters (diffusion D, pseudo-diffusion D*, and perfusion fraction F). The results were compared to those obtained by least squares (LSQ) on both simulated and in vivo brain data. Relative error (RE) and reproducibility were used to evaluate the results. The differences of IVIM parameters between brain tumor and normal regions were compared and used to assess the performance of Bayesian fitting in the IVIM application of brain tumor.

RESULTS

In tumor regions, the value of D* tended to be decreased when the number of low b-values was insufficient, especially with LSQ. BSP required less low b-values than LSQ for the correct estimation of perfusion parameters of brain tumors. The IVIM parameter maps of brain tumors yielded by BSP had smaller variability, lower RE, and higher reproducibility with respect to those obtained by LSQ. Obvious differences were observed between tumor and normal regions in parameters D (P < 0.05) and F (P < 0.001), especially F. BSP generated fewer outliers than LSQ, and distinguished better tumors from normal regions in parameter F.

CONCLUSIONS

Intravoxel incoherent motion parameters clearly allow brain tumors to be differentiated from normal regions. Bayesian fitting yields robust IVIM parameter mapping with fewer outliers and requires less low b-values than LSQ for the parameter estimation.

Comparing mono-exponential, bi-exponential, and stretched-exponential diffusion-weighted MR imaging for stratifying non-alcoholic fatty liver disease in a rabbit model

OBJECTIVES

To compare diffusion parameters obtained from mono-exponential, bi-exponential, and stretched-exponential diffusion-weighted imaging (DWI) in stratifying non-alcoholic fatty liver disease (NAFLD).

METHODS

Thirty-two New Zealand rabbits were fed a high-fat/cholesterol or standard diet to obtain different stages of NAFLD before 12 b-values (0-800 s/mm²) DWI. The apparent diffusion coefficient (ADC) from the mono-exponential model; pure water diffusion (D), pseudo-diffusion (D*), and perfusion fraction (f) from bi-exponential DWI; and distributed diffusion coefficient (DDC) and water molecular diffusion heterogeneity index (α) from stretched-exponential DWI were calculated for hepatic parenchyma. The goodness of fit of the three models was compared. NAFLD severity was pathologically graded as normal, simple steatosis, borderline, and non-alcoholic steatohepatitis (NASH). Spearman rank correlation analysis and receiver operating characteristic curves were used to assess NAFLD severity.

RESULTS

Upon comparison, the goodness of fit chi-square from stretched-exponential fitting (0.077 ± 0.012) was significantly lower than that for the bi-exponential (0.110 ± 0.090) and mono-exponential (0.181 ± 0.131) models (p < 0.05). Seven normal, 8 simple steatosis, 6 borderline, and 11 NASH livers were pathologically confirmed from 32 rabbits. Both α and D increased with increasing NAFLD severity (r = 0.811 and 0.373, respectively; p < 0.05). ADC, f, and DDC decreased as NAFLD severity increased (r = - 0.529, - 0.717, and - 0.541, respectively; p < 0.05). Both α (area under the curve [AUC] = 0.952) and f (AUC = 0.931) had significantly greater AUCs than ADC (AUC = 0.727) in the differentiation of NASH from borderline or less severe groups (p < 0.05).

CONCLUSIONS

Stretched-exponential DWI with higher fitting efficiency performed, as well as bi-exponential DWI, better than mono-exponential DWI in the stratification of NAFLD severity.

KEY POINTS

• Stretched-exponential diffusion model fitting was more reliable than the bi-exponential and mono-exponential diffusion models (p = 0.039 and p < 0.001, respectively). • As NAFLD severity increased, the diffusion heterogeneity index (α) increased, while the perfusion fraction (f) decreased (r = 0.811, - 0.717, p < 0.05). • Both α and f showed superior NASH diagnostic performance (AUC = 0.952, 0.931) compared with ADC (AUC = 0.727, p < 0.05).

Estimation of intravoxel incoherent motion parameters using low b-values

Intravoxel incoherent motion (IVIM) imaging is a magnetic resonance imaging (MRI) technique widely used in clinical applications for various organs. However, IVIM imaging at low b-values is a persistent problem. This paper aims to investigate in a systematic and detailed manner how the number of low b-values influences the estimation of IVIM parameters. To this end, diffusion-weighted (DW) data with different low b-values were simulated to get insight into the distributions of subsequent IVIM parameters. Then, in vivo DW data with different numbers of low b-values and different number of excitations (NEX) were acquired. Finally, least-squares (LSQ) and Bayesian shrinkage prior (BSP) fitting methods were implemented to estimate IVIM parameters. The influence of the number of low b-values on IVIM parameters was analyzed in terms of relative error (RE) and structural similarity (SSIM). The results showed that the influence of the number of low b-values on IVIM parameters is variable. LSQ is more dependent on the number of low b-values than BSP, but the latter is more sensitive to noise.

Determination of Hepatocellular Carcinoma and Characterization of Hepatic Focal Lesions with Adaptive Multi-Exponential Intravoxel Incoherent Motion Model

PURPOSE: To distinguish hepatocellular carcinoma (HCC) from other types of hepatic lesions with the adaptive multi-exponential IVIM model. METHODS: 94 hepatic focal lesions, including 38 HCC, 16 metastasis, 12 focal nodular hyperplasia, 13 cholangiocarcinoma, and 15 hemangioma, were examined in this study. Diffusion-weighted images were acquired with 13 b values (b = 0, 3, …, 500 s/mm²) to measure the adaptive multi-exponential IVIM parameters, namely, pure diffusion coefficient (D), diffusion fraction (f_d), pseudo-diffusion coefficient (D_i*) and perfusion-related diffusion fraction (f_i) of the ith perfusion component. Comparison of the parameters of and their diagnostic performance was determined using Mann-Whitney U test, independent-sample t test, one-way analysis of variance, Z test and receiver-operating characteristic analysis. RESULTS: D, D₁* and D₂* presented significantly difference between HCCs and other hepatic lesions, whereas f_d, f₁ and f₂ did not show statistical differences. In the differential diagnosis of HCCs from other hepatic lesions, D₂* (AUC, 0.927) provided best diagnostic performance among all parameters. Additionally, the number of exponential terms in the model was also an important indicator for distinguishing HCCs from other hepatic lesions. In the benign and malignant analysis, D gave the greatest AUC values, 0.895 or 0.853, for differentiation between malignant and benign lesions with three or two exponential terms. Most parameters were not significantly different between hypovascular and hypervascular lesions. For multiple comparisons, significant differences of D, D₁* or D₂* were found between certain lesion types. CONCLUSION: The adaptive multi-exponential IVIM model was useful and reliable to distinguish HCC from other hepatic lesions.

Intravoxel incoherent motion MRI as a means to measure in vivo perfusion: A review of the evidence

The idea that in vivo intravoxel incoherent motion magnetic resonance signal is influenced by blood motion in the microvasculature is exciting, because it suggests that local and quantitative perfusion information can be obtained in a simple and elegant way from a few diffusion-weighted images, without contrast injection. When the method was proposed in the late 1980s some doubts appeared as to its feasibility, and, probably because the signal to noise and image quality at the time was not sufficient, no obvious experimental evidence could be produced to alleviate them. Helped by the tremendous improvements seen in the last three decades in MR hardware, pulse design, and post-processing capabilities, an increasing number of encouraging reports on the value of intravoxel incoherent motion perfusion imaging have emerged. The aim of this article is to review the current published evidence on the feasibility of in vivo perfusion imaging with intravoxel incoherent motion MRI.