A Standardized Parameter-Free Algorithm for Combined Intravoxel Incoherent Motion and Diffusion Kurtosis Analysis of Diffusion Imaging Data. Invest Radiol. 2016;51:203-210. [PMID: 26561050 DOI: 10.1097/rli.0000000000000223]

Influence of Magnetic Field Strength on Intravoxel Incoherent Motion Parameters in Diffusion MRI of the Calf

Background: The purpose of this study was to investigate the dependence of Intravoxel Incoherent Motion (IVIM) parameters measured in the human calf on B0. Methods: Diffusion-weighted image data of eight healthy volunteers were acquired using five b-values (0-600 s/mm2) at rest and after muscle activation at 0.55 and 7 T. The musculus gastrocnemius mediale (GM, activated) was assessed. The perfusion fraction f and diffusion coefficient D were determined using segmented fits. The dependence on field strength was assessed using Student's t-test for paired samples and the Wilcoxon signed-rank test. A biophysical model built on the three non-exchanging compartments of muscle, venous blood, and arterial blood was used to interpret the data using literature relaxation times. Results: The measured perfusion fraction of the GM was significantly lower at 7 T, both for the baseline measurement and after muscle activation. For 0.55 and 7 T, the mean f values were 7.59% and 3.63% at rest, and 14.03% and 6.92% after activation, respectively. The biophysical model estimations for the mean proton-density-weighted perfusion fraction were 3.37% and 6.50% for the non-activated and activated states, respectively. Conclusions: B0 may have a significant effect on the measured IVIM parameters. The blood relaxation times suggest that 7 T IVIM may be arterial-weighted whereas 0.55 T IVIM may exhibit an approximately equal weighting of arterial and venous blood.

Influence of saturation effects on biexponential liver intravoxel incoherent motion

PURPOSE

Studies on intravoxel incoherent motion (IVIM) imaging in the liver have been carried out with different acquisition protocols. The number of acquired slices and the distances between slices can influence IVIM measurements due to saturation effects, but these effects have often been disregarded. This study investigated differences in biexponential IVIM parameters between two slice settings.

METHODS

Fifteen healthy volunteers (21-30 years) were examined at a field strength of 3 T. Diffusion-weighted images of the abdomen were acquired with 16 b values (0-800 s/mm² ), with four slices for the few slices setting and 24-27 slices for the many slices setting. Regions of interest were manually drawn in the liver. The data were fitted with a monoexponential signal curve and a biexponential IVIM curve, and biexponential IVIM parameters were determined. The dependence on the slice setting was assessed with Student's t test for paired samples (normally distributed IVIM parameters) and the Wilcoxon signed-rank test (non-normally distributed parameters).

RESULTS

None of the parameters were significantly different between the settings. For few slices and many slices, respectively, the mean values (SDs) for D $$ D $$ were 1.21 μm 2 / ms $$ 1.21{\upmu \mathrm{m}}^2/\mathrm{ms} $$ ( 0.19 μm 2 / ms $$ 0.19\kern0.3em {\upmu \mathrm{m}}^2/\mathrm{ms} $$ ) and 1.20 μm 2 / ms $$ 1.20{\upmu \mathrm{m}}^2/\mathrm{ms} $$ ( 0.11 μm 2 / ms $$ 0.11\kern0.3em {\upmu \mathrm{m}}^2/\mathrm{ms} $$ ); for f $$ f $$ they were 29.7% (6.2%) and 27.7% (3.6%); and for D * $$ {D}^{\ast } $$ they were 8.76 ⋅ 10 - 2 mm 2 / s $$ 8.76\cdot {10}^{-2}{\mathrm{mm}}^2/\mathrm{s} $$ ( 4.54 ⋅ 10 - 2 mm 2 / s $$ 4.54\cdot {10}^{-2}\kern0.3em {\mathrm{mm}}^2/\mathrm{s} $$ ) and 8.71 ⋅ 10 - 2 mm 2 / s $$ 8.71\cdot {10}^{-2}{\mathrm{mm}}^2/\mathrm{s} $$ ( 4.06 ⋅ 10 - 2 mm 2 / s $$ 4.06\cdot {10}^{-2}\kern0.3em {\mathrm{mm}}^2/\mathrm{s} $$ ).

CONCLUSION

Biexponential IVIM parameters in the liver are comparable among IVIM studies that use different slice settings, with mostly negligible saturation effects. However, this may not hold for studies that use much shorter TR.

Echo time dependence of biexponential and triexponential intravoxel incoherent motion parameters in the liver

PURPOSE

Intravoxel incoherent motion (IVIM) studies are performed with different acquisition protocols. Comparing them requires knowledge of echo time (TE) dependencies. The TE-dependence of the biexponential perfusion fraction f is well-documented, unlike that of its triexponential counterparts f₁ and f₂ and the biexponential and triexponential pseudodiffusion coefficients D^* , D 1 ∗ , and D 2 ∗ . The purpose was to investigate the TE-dependence of these parameters and to check whether the triexponential pseudodiffusion compartments are associated with arterial and venous blood.

METHODS

Fifteen healthy volunteers (19-58 y; mean: 24.7 y) underwent diffusion-weighted imaging of the abdomen with 24 b-values (0.2-800 s/mm² ) at TEs of 45, 60, 75, and 90 ms. Regions of interest (ROIs) were manually drawn in the liver. One set of bi- and triexponential IVIM parameters per volunteer and TE was determined. The TE-dependence was assessed with the Kruskal-Wallis test.

RESULTS

TE-dependence was observed for f (P < .001), f₁ (P = .001), and f₂ (P < .001). Their median values at the four measured TEs were: f: 0.198/0.240/0.274/0.359, f₁ : 0.113/0.139/0.146/0.205, f₂ : 0.115/0.155/0.182/0.194. D, D^* , D 1 ∗ , and D 2 ∗ showed no significant TE-dependence. Their values were: diffusion coefficient D (10^-4 mm² /s): 9.45/9.63/9.75/9.41, biexponential D^* (10^-2 mm² /s): 5.26/5.52/6.13/5.82, triexponential D 1 ∗ (10^-2 mm² /s): 1.73/2.91/2.25/2.51, triexponential D 2 ∗ (mm² /s): 0.478/1.385/0.616/0.846.

CONCLUSION

f₁ and f₂ show similar TE-dependence as f, ie, increase with rising TE; an effect that must be accounted for when comparing different studies. The diffusion and pseudodiffusion coefficients might be compared without TE correction. Because of the similar TE-dependence of f₁ and f₂ , the triexponential pseudodiffusion compartments are most probably not associated to venous and arterial blood.

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.

Understanding diffusion-weighted MRI analysis: Repeatability and performance of diffusion models in a benign breast lesion cohort

Diffusion-weighted MRI (DWI) is an important tool for oncology research, with great clinical potential for the classification and monitoring of breast lesions. The utility of parameters derived from DWI, however, is influenced by specific analysis choices. The purpose of this study was to critically evaluate repeatability and curve-fitting performance of common DWI signal representations, for a prospective cohort of patients with benign breast lesions. Twenty informed, consented patients with confirmed benign breast lesions underwent repeated DWI (3 T) using: sagittal single-shot spin-echo echo planar imaging, bipolar encoding, TR/TE: 11,600/86 ms, FOV: 180 x 180 mm, matrix: 90 x 90, slices: 60 x 2.5 mm, iPAT: GRAPPA 2, fat suppression, and 13 b-values: 0-700 s/mm² . A phase-reversed scan (b = 0 s/mm² ) was acquired for distortion correction. Voxel-wise repeat-measures coefficients of variation (CoVs) were derived for monoexponential (apparent diffusion coefficient [ADC]), biexponential (intravoxel incoherent motion: f, D, D*) and stretched exponential (α, DDC) across the parameter histograms for lesion regions of interest (ROIs). Goodness-of-fit for each representation was assessed by Bayesian information criterion. The volume of interest (VOI) definition was repeatable (CoV 13.9%). Within lesions, and across both visits and the cohort, there was no dominant best-fit model, with all representations giving the best fit for a fraction of the voxels. Diffusivity measures from the signal representations (ADC, D, DDC) all showed good repeatability (CoV < 10%), whereas parameters associated with pseudodiffusion (f, D*) performed poorly (CoV > 50%). The stretching exponent α was repeatable (CoV < 12%). This pattern of repeatability was consistent over the central part of the parameter percentiles. Assumptions often made in diffusion studies about analysis choices will influence the detectability of changes, potentially obscuring useful information. No single signal representation prevails within or across lesions, or across repeated visits; parameter robustness is therefore a critical consideration. Our results suggest that stretched exponential representation is more repeatable than biexponential, with pseudodiffusion parameters unlikely to provide clinically useful biomarkers.

MRI phantom for tissue simulation with respect to diffusion coefficient and kurtosis - Validation with injection of liposomal theranostics

This study presents gelatine-based and agar-based phantoms with an addition of glycerol, safflower oil, silicone oil and cellulose microcrystalline with a potential to cover the entire range of tissue diffusion coefficients and kurtosis values. Forty types of phantoms were prepared and examined for NMR relaxation times T₁ and T₂ and diffusional metrics D, K and ADC. Wide ranges of values of D (0.0003-0.0031 mm²s^-1), K (0.00-7.24) and ADC (0.0002-0.0031 mm²s^-1) were observed. Two of the phantoms closely mimic muscle and cortical gray matter with respect to water diffusion parameters. Although many of the presented phantoms display both D and K values within the range of human tissues, they match different tissues with respect to D and K. The imaging results for the gray matter simulating phantom injected with the liposomal solution, bear a resemblance to the particle size effect described in the literature. The phantoms presented in this work are simple in preparation and affordable tissue-simulating materials to be used primarily in development of diffusion kurtosis-based MRI methods and possibly in a preliminary assessment of MRI contrast agents. Further adjustments of the chemical compositions could potentially lead to development of new types of phantoms mimicking diffusional properties of more kinds of soft tissues.

IVIM-DKI for differentiation between prostate cancer and benign prostatic hyperplasia: comparison of 1.5 T vs. 3 T MRI

OBJECTIVE

To implement an advanced spatial penalty-based reconstruction to constrain the intravoxel incoherent motion (IVIM)-diffusion kurtosis imaging (DKI) model and investigate whether it provides a suitable alternative at 1.5 T to the traditional IVIM-DKI model at 3 T for clinical characterization of prostate cancer (PCa) and benign prostatic hyperplasia (BPH).

MATERIALS AND METHODS

Thirty-two patients with biopsy-proven PCa were recruited for MRI examination (n = 16 scanned at 1.5 T, n = 16 scanned at 3 T). Diffusion-weighted imaging (DWI) with 13 b values (b = 0 to 2000 s/mm² up to 3 averages, 1.5 T: TR = 5.774 s, TE = 81 ms and 3 T: TR = 4.899 s, TE = 100 ms), T2-weighted, and T1-weighted imaging were used on the 1.5 T and 3 T MRI scanner, respectively. The IVIM-DKI signal was modeled using the traditional IVIM-DKI model and a novel model in which the total variation (TV) penalty function was combined with the traditional model to optimize non-physiological variations. Paired and unpaired t-tests were used to compare intra-scanner and scanner group differences in IVIM-DKI parameters obtained using the novel and the traditional models. Analysis of variance with post hoc test and receiver operating characteristic (ROC) curve analysis were used to assess the ability of parameters obtained using the novel model (at 1.5 T) and the traditional model (at 3 T) to characterize prostate lesions.

RESULTS

IVIM-DKI modeled using novel model with TV spatial penalty function at 1.5 T, produced parameter maps with 50-78% lower coefficient of variation (CV) than traditional model at 3 T. Novel model estimated higher D with lower D*, f and k values at both field strengths compared to traditional model. For scanner differences, the novel model at 1.5 T estimated lower D* and f values as compared to traditional model at 3 T. At 1.5 T, D and f values were significantly lower with k values significantly higher in tumor than BPH and healthy tissue. D (AUC: 0.98), f (AUC: 0.82), and k (AUC: 0.91) parameters estimated using novel model showed high diagnostic performance in cancer lesion detection at 1.5 T.

DISCUSSION

In comparison with the IVIM-DKI model at 3 T, IVIM-DKI signal modeled with the TV penalty function at 1.5 T showed lower estimation errors. The proposed novel model can be utilized for improved detection of prostate lesions.

The value of intravoxel incoherent motion imaging in predicting the survival of patients with astrocytoma

BACKGROUND

The evaluation of the prognosis of gliomas may have great value in individualized treatment.

PURPOSE

To evaluate the value of intravoxel incoherent motion (IVIM) in predicting the survival of patients with astrocytoma and comparing it to apparent diffusion coefficients (ADC).

MATERIAL AND METHODS

Sixty patients with pathologically confirmed cerebral astrocytomas underwent IVIM scans before any treatment was performed. Patients were divided into death group and survival group according to a two-year follow-up. ADC and quantitative parameters of IVIM including D, D*, and f were measured. Independent sample t test was used to compare the two groups of parameters. The accuracy of each parameter for two-year survival rate was analyzed by receiver operating characteristic (ROC) curve and Kaplan-Meier survival curves. The correlation between quantitative parameters and survival days was analyzed by Pearson correlation analysis.

RESULTS

The ADC, D*, and f values were statistically significant different between the death and the survival groups (P < 0.05). The AUC of the ADC, D*, and f were 0.811, 0.858, and 0.892, respectively. The ADC cut-off value of 0.668 × 10^-3 mm²/s corresponded to 82.6% sensitivity and 73% specificity. The D* cut-off value of 3.913 × 10^-3 mm²/s corresponded to 78.4% sensitivity and 87% specificity. The f cut-off value of 0.487 corresponded to 83.8% sensitivity and 87% specificity. Significant log rank test was performed for each parameter to predict overall survival (P < 0.05). There was a correlation between ADC (r = 0.625, P = 0.023), D* (r = -0.655, P = 0.012), f (r = -0.725, P = 0.000) and survival days.

CONCLUSION

The D* and f values demonstrated great potential in predicting the two-year survival rate for patients with astrocytoma.

Analysis of Renal Diffusion-Weighted Imaging (DWI) Using Apparent Diffusion Coefficient (ADC) and Intravoxel Incoherent Motion (IVIM) Models

Analysis of renal diffusion-weighted imaging (DWI) data to derive markers of tissue properties requires careful consideration of the type, extent, and limitations of the acquired data. Alongside data quality and general suitability for quantitative analysis, choice of diffusion model, fitting algorithm, and processing steps can have consequences for the precision, accuracy, and reliability of derived diffusion parameters. Here we introduce and discuss important steps for diffusion-weighted image processing, and in particular give example analysis protocols and pseudo-code for analysis using the apparent diffusion coefficient (ADC) and intravoxel incoherent motion (IVIM) models. Following an overview of general principles, we provide details of optional steps, and steps for validation of results. Illustrative examples are provided, together with extensive notes discussing wider context of individual steps, and notes on potential pitfalls.This publication is based upon work from the COST Action PARENCHIMA, a community-driven network funded by the European Cooperation in Science and Technology (COST) program of the European Union, which aims to improve the reproducibility and standardization of renal MRI biomarkers. This analysis protocol chapter is complemented by two separate chapters describing the basic concepts and experimental procedure.

Advanced diffusion imaging of abdominal organs in different hydration states of the human body: stability of biomarkers

BACKGROUND

MR diffusion weighted imaging (DWI) may provide important information regarding the pathophysiology of parenchymal abdominal organs. The purpose of our study was to investigate the stability of imaging biomarkers of diffusion weighted imaging (DWI), intravoxel incoherent motion (IVIM) and diffusion kurtosis imaging (DKI) in abdominal parenchymal organs regarding two body hydration states.

METHODS

Ten healthy volunteers twice underwent DWI of abdominal organs using a double-refocused spin-echo echo-planar imaging sequences with 11 different b-values (ranging from 0 to 1,500 s/mm2): after 4 h of fluid deprivation; 45 min following 1000 ml of water intake. Four different diffusion models were evaluated and compared: standard DWI, DKI with mono-exponential fitting, multistep algorithm with variable b-value threshold for IVIM, combined IVIM-Kurtosis; in four abdominal organs: kidneys, liver, spleen and psoas muscle.

RESULTS

Diffusion parameters from all four models remained similar for the renal parenchyma before and after the water challenge. Significant differences were found for the liver, spleen, and psoas muscle. The largest effects were seen for: the liver parenchyma after the water challenge by means of IVIM model's true diffusion (p < 0.02); the spleen, for IVIM's perfusion fraction (p < 0.03), the psoas muscle for the ADC value (p < 0.02).

CONCLUSIONS

Herein, we showed that diffusion parameters of the kidney remain remarkably stable regarding the hydration status. This may be attributed to the kidney-specific compensatory mechanisms. For the liver, spleen and psoas muscle the diffusion parameters were sensitive to changes of the hydration. This phenomenon needs to be considered when evaluating diffusion data of these organs.

Role of intravoxel incoherent motion parameters in gastroesophageal cancer: relationship with 18F-FDG-positron emission tomography, computed tomography perfusion and magnetic resonance perfusion imaging parameters

BACKGROUND

Identification of pretherapeutic predictive markers in gastro-esophageal cancer is essential for individual-oriented treatment. This study evaluated the relationship of multimodality parameters derived from intravoxel incoherent motion method (IVIM), 18F-FDG-positron emission tomography (PET), computed tomography (CT) perfusion and dynamic contrast enhanced magnetic resonance imaging (MRI) in patients with gastro-esophageal cancer and investigated their histopathological correlation.

METHODS

Thirty-one consecutive patients (28 males; median age 63.9 years; range 37-84 years) with gastro-esophageal adenocarcinoma (N.=22) and esophageal squamous cell carcinoma (N.=9) were analyzed. IVIM parameters: pseudodiffusion (D*), perfusion fraction (fp), true diffusion (D) and the threshold b-value (bval); PET-parameters: SUV<inf>max</inf>, metabolic tumor volume (MTV) and total lesion glycolysis (TLG); CT perfusion parameters: blood flow (BF), blood volume (BV) and mean transit time (MTT); and MR perfusion parameters: time to enhance, positive enhancement integral, time-to-peak (TTP), maximum-slope-of-increase, and maximum-slope-of-decrease were determined, and correlated to each other and to histopathology.

RESULTS

IVIM and PET parameters showed significant negative correlations: MTV and bval (r<inf>s</inf> =-0.643, P=0.002), TLG and bval (r<inf>s</inf>=-0.699, P<0.01) and TLG and fp (r<inf>s</inf>=-0.577, P=0.006). Positive correlation was found for TLG and D (r<inf>s</inf>=0.705, P=0.000). Negative correlation was found for bval and staging (r<inf>s</inf>=0.590, P=0.005). Positive correlation was found for positive enhancement interval and BV (r<inf>s</inf>=0.547, P=0.007), BF and regression index (r<inf>s</inf>=0.753, P=0.005) and for time-to-peak and staging (r<inf>s</inf>=0.557, P=0.005).

CONCLUSIONS

IVIM parameters (bval, fp, D) provide quantitative information and correlate with PET parameters (MTV, TLG) and staging. IVIM might be a useful tool for additional characterization of gastro-esophageal cancer.

Pancreatic imaging: Current status of clinical practices and small animal studies

Different causative factors acting on the pancreas can result in diseases such as pancreatitis, diabetes and pancreatic tumors. The high incidence and mortality of pancreatic diseases have placed diagnostic imaging in a crucial position in daily clinical practice. In this mini-review article different pancreatic imaging techniques are discussed, from the standard clinical imaging modalities and state of the art clinical magnetic resonance imaging techniques to current situations in pre-clinical pancreatic imaging studies. In particular, the challenges of pre-clinical rodent pancreatic imaging are addressed, with both the image acquisition techniques and the post-processing methods for rodent pancreatic imaging elaborated.

Liver intravoxel incoherent motion (IVIM) magnetic resonance imaging: a comprehensive review of published data on normal values and applications for fibrosis and tumor evaluation

A comprehensive literature review was performed on liver intravoxel incoherent motion (IVIM) magnetic resonance imaging (MRI) technique and its applications. Heterogeneous data have been reported. IVIM parameters are magnetic field strength dependent to a mild extent. A lower Dslow (D) value at 3 T than at 1.5 T and higher perfusion fraction (PF) value at 3 T than at 1.5 T were noted. An increased number of b values are associated with increased IVIM parameter measurement accuracy. With the current status of art, IVIM technique is not yet capable of detecting early stage liver fibrosis and diagnosing liver fibrosis grades, nor can it differentiate liver tumors. Though IVIM parameters show promise for tumor treatment monitoring, till now how PF and Dfast (D*) add diagnostic value to Dslow or apparent diffusion coefficient (ADC) remains unclear. This paper shows the state-of-art IVIM MR technique is still not able to offer reliable measurement for liver. More works on the measurement robustness are warranted as they are essential to justify follow-up clinical studies on patients.