Intravoxel incoherent motion diffusion-weighted MRI of the abdomen: The effect of fitting algorithms on the accuracy and reliability of the parameters

IVIM parameters mapping with artificial neural network based on mean deviation prior

BACKGROUND

The diffusion and perfusion parameters derived from intravoxel incoherent motion (IVIM) imaging provide promising biomarkers for noninvasively quantifying and managing various diseases. Nevertheless, due to the distribution gap between simulated and real datasets, the out-of-distribution (OOD) problem occurred in supervised learning-based methods degrades their performance and hinders their real applications.

PURPOSE

To address the OOD problem in supervised methods and to further improve the accuracy and stability of IVIM parameter estimation, this work proposes a novel learning framework called IterANN, based on mean deviation prior (MDP) between training and estimated IVIM parameters on the test set.

METHODS

Specifically, MDP indicates that the mean of the estimated IVIM parameters always locates between the mean of IVIM parameters in the test and train sets. In IterANN, we adopt a very simple artificial neural network (ANN) architecture of two hidden layers with 12 neurons per hidden layer, an input layer containing the signals acquired at multiple b-values and an output layer composed of three IVIM parameters ( D $D$ , F $F$ andD S t a r $DStar$ ). Inspired by MDP, the distribution of IVIM parameters in the training set (simulated data) is iteratively updated so that their mean gradually approaches the predicted values of the real data. This aims to achieve a strong correlation between the simulated data and the real data. To validate the effectiveness of IterANN, we compare it with several methods on both simulation and real acquisition datasets, including 21 healthy and 3 tumor subjects, in terms of residual errors of IVIM parameters or DW signals, the coefficients of variation (CV) of IVIM parameters, and the parameter contrast-to-noise ratio (PCNR) between normal and tumor tissues.

RESULTS

On two simulation datasets, the proposed IterANN achieves the lowest residual error in IVIM parameters, especially in the case of low signal-to-noise ratio (SNR = 10), the residual error of D $D$ , F $F$ andD S t a r $DStar$ is decreased by15.82 % / 14.92 % , 81.19 % / 74.04 % , 50.77 % / 1.549 % $15.82\%/14.92\%, 81.19\%/74.04\%, 50.77\%/1.549\%$ (Gaussian distribution /realistic distribution) respectively comparing to the suboptimal method. On real dataset, the IterANN achieves the highest PCNR when comparing the normal and tumor regions. Additionally, the proposed IterANN demonstrated better stability, with its CV being significantly lower than that of other methods in the vast majority of cases (p < 0.01 $p<0.01$ , paired-sample Student's t-test).

CONCLUSIONS

The superior performance of IterANN demonstrates that updating the distribution of the train set based on MDP can effectively solve the OOD problem, which allows us not only to improve the accuracy and stability of the estimated IVIM parameters, but also to increase the potential of IVIM in disease diagnosis.

Effect of simultaneous multislice imaging, slice properties, and repetition time on the measured magnetic resonance biexponential intravoxel incoherent motion in the liver

OBJECTIVES

This study aims to investigate the previously reported dependency of intravoxel incoherent motion (IVIM) parameters on simultaneous multislice (SMS) acquisition and repetition time (TR). This includes the influence of slice thickness, slice gaps, and slice order on measured IVIM parameters.

MATERIALS AND METHODS

Diffusion-weighted imaging (DWI) of the liver was performed on 10 healthy volunteers (aged 20-30 years) at 3T with a slice thickness of 5 mm, a slice gap of 5 mm, and a linear slice order. Diffusion-weighted images were acquired with 19 b-values (0-800 s/mm2) using both conventional slice excitation with an acceleration factor of one (AF1) and SMS excitation with an acceleration factor of three (AF3). Each of these measurements were carried out with two repetition times (TRs)- 1,300 ms (prefix s) and 4,500 ms (prefix l)-resulting in four different combinations: sAF1, sAF3, lAF1, and lAF3. Five volunteers underwent additional measurements using a 10 mm slice thickness and with AF1. Median signal values in the liver were used to determine the biexponential IVIM parameters. Statistical significances were assessed using the Kruskal-Wallis test, Wilcoxon signed-rank test, and Student's t-test. In-silico investigations were also used to interpret the data.

RESULTS

There were no significant differences between the biexponential IVIM parameters acquired from sAF1, sAF3, lAF1, and lAF3. Median values of the perfusion fraction f were as follows: 29.9% (sAF1), 26.9% (sAF3), 28.1% (lAF1), and 27.5% (lAF3). In the 10 mm-thick slices, f decreased from 31.3% (lAF1) to 27.4% (sAF1) (p = 0.141).

CONCLUSION

The slice excitation mode did not appear to have any significant influence on the biexponential IVIM parameters. However, our simulations, as well as values reported from previous publications, show that slice thickness, slice gaps, and slice order are relevant and should thus be reported in IVIM studies.

Velocity-compensated intravoxel incoherent motion of the human calf muscle

PURPOSE

To determine whether intravoxel incoherent motion (IVIM) describes the blood perfusion in muscles better, assuming pseudo diffusion (Bihan Model 1) or ballistic motion (Bihan Model 2).

METHODS

IVIM parameters were measured in 18 healthy subjects with three different diffusion gradient time profiles (bipolar with two diffusion times and one with velocity compensation) and 17 b-values (0-600 s/mm2) at rest and after muscle activation. The diffusion coefficient, perfusion fraction, and pseudo-diffusion coefficient were estimated with a segmented fit in the gastrocnemius medialis (GM) and tibialis anterior (TA) muscles.

RESULTS

Velocity-compensated gradients resulted in a decreased perfusion fraction (6.9% ± 1.4% vs. 4.4% ± 1.3% in the GM after activation) and pseudo-diffusion coefficient (0.069 ± 0.046 mm2/s vs. 0.014 ± 0.006 in the GM after activation) compared to the bipolar gradients with the longer diffusion encoding time. Increased diffusion coefficients, perfusion fractions, and pseudo-diffusion coefficients were observed in the GM after activation for all gradient profiles. However, the increase was significantly smaller for the velocity-compensated gradients. A diffusion time dependence was found for the pseudo-diffusion coefficient in the activated muscle.

CONCLUSION

Velocity-compensated diffusion gradients significantly suppress the IVIM effect in the calf muscle, indicating that the ballistic limit is mostly reached, which is supported by the time dependence of the pseudo-diffusion coefficient.

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.

Repeatability and reproducibility comparisons of liver IVIM imaging with free-breathing or respiratory-triggered sequences

For liver intravoxel incoherent motion (IVIM) data acquisition, respiratory-triggering (RT) MRI is commonly used, and there are strong motivations to shorten the scan duration. For the same scan duration, more b values or higher numbers of excitations can be allowed for free-breathing (FB) imaging than for RT. We studied whether FB can be used to replace RT when careful IVIM image acquisition and image processing are conducted. MRI data of 22 healthy participants were acquired using a 3.0 T scanner. Diffusion imaging was based on a single-shot spin-echo-type echo-planar sequence and 16 b values of 0, 2, 4, 7, 10, 15, 20, 30, 46, 60, 72, 100, 150, 200, 400, and 600 s/mm2 . Each subject attended two scan sessions with an interval of 10-20 days. For each scan session, a subject was scanned twice, first with RT and then with FB. The mean image acquisition time was 5.4 min for FB and 10.8 min for RT. IVIM parameters were calculated with bi-exponential model segmented fitting with a threshold b value of 60 s/mm2 , and fitting started from b = 2 s/mm2 . There was no statistically significant difference between IVIM parameters measured with FB imaging or RT imaging. Perfusion fraction ICC (intraclass correlation coefficient) for FB imaging and RT imaging in the same scan session was 0.824. For perfusion fraction, wSD (within-subject standard deviation), BA (Bland-Altman) difference, BA 95% limit, and ICC were 0.022, 0.0001, -0.0635~0.0637, and 0.687 for FB and 0.031, 0.0122, -0.0723~0.0967, and 0.611 for RT. For Dslow (×10-3  s/mm2 ), wSD, BA difference, BA 95% limit, and ICC were 0.057, 0.0268, -0.1258~0.1793, and 0.471 for FB and 0.073, -0.0078, -0.2170-0.2014, and <0.4 for RT. The Dfast coefficient of variation was 0.20 for FB imaging and 0.28 for RT imaging. All reproducibility indicators slightly favored FB imaging.

Quantitative tumor burden imaging parameters of the spleen at MRI for predicting treatment response in patients with acute leukemia

Objectives

To study the value of standardized volume and intravoxel incoherent motion (IVIM) parameters of the spleen based on tumor burden for predicting treatment response in newly diagnosed acute leukemia (AL).

Methods

Patients with newly diagnosed AL were recruited and underwent abdominal IVIM diffusion-weighted imaging within one week before the first induction chemotherapy. Quantitative parameters of magnetic resonance imaging (MRI) included the standardized volume (representing volumetric tumor burden) and IVIM parameters (standard apparent diffusion coefficient [sADC]; pure diffusion coefficient [D]; pseudo-diffusion coefficient [D∗]; and pseudo-perfusion fraction [f], representing functional tumor burden) of the spleen. Clinical biomarkers of tumor burden were collected. Patients were divided into complete remission (CR) and non-CR groups according to the treatment response after the first standardized induction chemotherapy, and the MRI and clinical parameters were compared between the two groups. The correlations of MRI parameters with clinical biomarkers were analyzed. Multivariate logistic regression was performed to determine the independent predictors for treatment response. Receiver operating characteristic curves were used to analyze the predicted performance.

Results

76 AL patients (CR: n = 43; non-CR: n = 33) were evaluated. Standardized spleen volume, sADC, D, f, white blood cell counts, and lactate dehydrogenase were significantly different between CR and non-CR groups (all p < 0.05). Standardized spleen volume, sADC, and D were correlated with white blood cell and lactate dehydrogenase, and f was correlated with lactate dehydrogenase (all p < 0.05). Standardized spleen volume (hazard ratio = 4.055, p = 0.042), D (hazard ratio = 0.991, p = 0.027), and f (hazard ratio = 1.142, p = 0.008) were independent predictors for treatment response, and the combination of standardized spleen volume, D, and f showed more favorable discrimination (area under the curve = 0.856) than individual predictors.

Conclusion

Standardized volume, D, and f of the spleen could be used to predict treatment response in newly diagnosed AL, and the combination of morphological and functional parameters would further improve the predicted performance. IVIM parameters of the spleen may be viable indicators for evaluating functional tumor burden in AL.

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.

Three-dimensional amide proton transfer-weighted and intravoxel incoherent motion imaging for predicting bone metastasis in patients with prostate cancer: A pilot study

PURPOSE

To explore the value of 3-dimensional amide proton transfer-weighted (APTw) and intravoxel incoherent motion (IVIM) imaging in predicting bone metastasis (BM) of prostate cancer (PCa) in addition to routine diffusion-weighted imaging (DWI).

METHODS

The clinical and imaging data of 39 PCa patients who were pathologically confirmed in our hospital from March 2019 to February 2022 were retrospectively analyzed, and they were divided into BM-negative (27 patients) and BM-positive (12 patients) groups. MR examination included APTw, DWI and IVIM imaging. The IVIM data was fitted by single-exponential IVIM model (IVIM_mono) and double-exponential IVIM model (IVIM_bi), respectively. The APTw, ADC, IVIM_mono (D_mono, D*_mono, and f_mono), and IVIM_bi (D_bi, D*_bi, and f_bi) parameters were independently measured by two radiologists. The synthetic minority oversampling technique (SMOTE) was conducted to balance the minority group. Mann-Whitney U test or Student's t-test was used to compare above values between the BM-negative and BM-positive groups. The diagnostic performance was evaluated with receiver operating characteristic (ROC) analysis of each parameter and their combination. The Delong test was used for ROC curve comparison.The relationship between APTw and IVIM was explored through Spearman's rank correlation analysis.

RESULTS

The APTw and D*_mono values were higher, and the ADC, f_mono, and f_bi values were lower in the BM-positive group than in the BM-negative group (all P < 0.05). Among the individual parameters, the AUC of f_mono was the highest (AUC = 0.865), and AUC (f_mono) was significantly higher than AUC (f_bi), AUC (D*_mono), and AUC (ADC) (all P < 0.05). The AUC (IVIM_mono) was higher than the AUC (IVIM_bi) (P = 0.0068). The combination of APTw and IVIM_mono further improved diagnostic capability, and the AUC of APTw+IVIM_mono was significantly higher than those of APTw and DWI (all P < 0.05). No correlation was found between IVIM-derived parameters and APTw value.

CONCLUSION

Both 3D APTw and IVIM imaging could predict BM of PCa. IVIM showed better performance than APTw and DWI, and the single-exponential IVIM model was superior to the double-exponential IVIM model. The combination of APTw and IVIM could further improve diagnostic performance.

Differential value of diffusion kurtosis imaging and intravoxel incoherent motion in benign and malignant solitary pulmonary lesions

Objective

To investigate the diagnostic value of diffusion kurtosis imaging (DKI) and intravoxel incoherent motion (IVIM) whole-lesion histogram parameters in differentiating benign and malignant solitary pulmonary lesions (SPLs).

Materials and Methods

Patients with SPLs detected by chest CT examination and with further routine MRI, DKI and IVIM-DWI functional sequence scanning data were recruited. According to the pathological results, SPLs were divided into a benign group and a malignant group. Independent samples t tests (normal distribution) or Mann‒Whitney U tests (nonnormal distribution) were used to compare the differences in DKI (Dk, K), IVIM (D, D*, f) and ADC whole-lesion histogram parameters between the benign and malignant SPL groups. The receiver operating characteristic (ROC) curve was used to evaluate the diagnostic efficiency of the histogram parameters and determine the optimal threshold. The area under the curve (AUC) of each histogram parameter was compared by the DeLong method. Spearman rank correlation was used to analyze the correlation between histogram parameters and malignant SPLs.

Results

Most of the histogram parameters for diffusion-related values (Dk, D, ADC) of malignant SPLs were significantly lower than those of benign SPLs, while most of the histogram parameters for the K value of malignant SPLs were significantly higher than those of benign SPLs. DKI (Dk, K), IVIM (D) and ADC were effective in differentiating benign and malignant SPLs and combined with multiple parameters of the whole-lesion histogram for the D value, had the highest diagnostic efficiency, with an AUC of 0.967, a sensitivity of 90.00% and a specificity of 94.03%. Most of the histogram parameters for the Dk, D and ADC values were negatively correlated with malignant SPLs, while most of the histogram parameters for the K value were positively correlated with malignant SPLs.

Conclusions

DKI (Dk, K) and IVIM (D) whole-lesion histogram parameters can noninvasively distinguish benign and malignant SPLs, and the diagnostic performance is better than that of DWI. Moreover, they can provide additional information on SPL microstructure, which has important significance for guiding clinical individualized precision diagnosis and treatment and has potential clinical application value.

A supervised deep neural network approach with standardized targets for enhanced accuracy of IVIM parameter estimation from multi-SNR images

Extraction of intravoxel incoherent motion (IVIM) parameters from noisy diffusion-weighted (DW) images using a biexponential fitting model is computationally challenging, and the reliability of the estimated perfusion-related quantities represents a limitation of this technique. Artificial intelligence can overcome the current limitations and be a suitable solution to advance use of this technique in both preclinical and clinical settings. The purpose of this work was to develop a deep neural network (DNN) approach, trained on numerical simulated phantoms with different signal to noise ratios (SNRs), to improve IVIM parameter estimation. The proposed approach is based on a supervised fully connected DNN having 3 hidden layers, 18 inputs and 3 targets with standardized values. 14 × 10³ simulated DW images, based on a Shepp-Logan phantom, were randomly generated with varying SNRs (ranging from 10 to 100). 7 × 10³ images (1000 for each SNR) were used for training. Performance accuracy was assessed in simulated images and the proposed approach was compared with the state-of-the-art Bayesian approach and other DNN algorithms. The DNN approach was also evaluated in vivo on a high-field MRI preclinical scanner. Our DNN approach showed an overall improvement in accuracy when compared with the Bayesian approach and other DNN methods in most of the simulated conditions. The in vivo results demonstrated the feasibility of the proposed approach in real settings and generated quantitative results comparable to those obtained using the Bayesian and unsupervised approaches, especially for D and f, and with lower variability in homogeneous regions. The DNN architecture proposed in this work outlines two innovative features as compared with other studies: (1) the use of standardized targets to improve the estimation of parameters, and (2) the implementation of a single DNN to enhance the IVIM fitting at different SNRs, providing a valuable alternative tool to compute IVIM parameters in conditions of high background noise.

Bi-exponential fitting excluding b=0 data improves the scan-rescan stability of liver IVIM parameter measures and particularly so for the perfusion fraction

BACKGROUND

A prerequisite to translating intravoxel incoherent motion (IVIM) imaging into meaningful clinical applications is sufficient scan-rescan reproducibility. This study aims to confirm the hypothesis that IVIM data fitting by not using b=0 images will improve the stability of liver IVIM measurement.

METHODS

Healthy volunteers' liver IVIM images were prospectively acquired using a 1.5-T magnet or a 3.0 T with 16 b-values. Repeatability study subjects were scanned twice during the same session, resulted in 35 paired scans for 35 subjects (11 men, mean age: 41.82 years, range: 32-60 years; 24 women, mean age: 42.67 years, range: 20-71 years). IVIM analysis was performed with full-fitting and segmented-fitting with a threshold b-value of 60 s/mm2, and fitting started from b=0 s/mm2 or from b=2 s/mm2. Reproducibility study subjects were scanned and then rescanned with an interval of 5-18 days, resulted in 20 paired scans for 11 subjects (4 men, mean age: 26.25 years, range: 25-27 years; 7 women, mean age: 25.57 years, range: 24-27 years). IVIM analysis was performed with segmented-fitting with a threshold b-value of 50 s/mm2, and fitting started from b=0 s/mm2 or from b=3 s/mm2.

RESULTS

Fitting without b=0 data generally improved the repeatability and reproducibility for both PF and Dslow, and particularly so for PF. For with b=0 data segmented fitting repeatability, PF had within-subject standard deviation of 0.019, bland-Atman 75% agreement limit of -31.52% to 28.35%, and ICC of 0.647, while these values were 0.009, -20.78% to 16.86%, and 0.837 for without b=0 analysis. Though the repeatability and reproducibility for Dfast generally also improved, they remained suboptimal. Measurement stability was better for repeatability than for reproducibility.

CONCLUSIONS

Scan-rescan repeatability and reproducibility of liver IVIM parameters can be improved by fitting without b=0 data, which is particularly so for PF.

Advances in Imaging in Evaluating the Efficacy of Neoadjuvant Chemotherapy for Breast Cancer

Neoadjuvant chemotherapy (NAC) is increasingly widely used in breast cancer treatment, and accurate evaluation of its response provides essential information for treatment and prognosis. Thus, the imaging tools used to quantify the disease response are critical in evaluating and managing patients treated with NAC. We discussed the recent progress, advantages, and disadvantages of common imaging methods in assessing the efficacy of NAC for breast cancer.

Perfusion-Diffusion Ratio: A New IVIM Approach in Differentiating Solid Benign and Malignant Primary Lesions of the Liver

Introduction

In order to improve the efficacy of intravoxel incoherent motion (IVIM) parameters in characterising specific tissues, a new concept is introduced: the perfusion–diffusion ratio (PDR), which expresses the relationship between the signal S(b) decline rate as a result of IVIM and the rate of signal S(b) decline due to diffusion. The aim of this study was to investigate this novel approach in the differentiation of solid primary liver lesions. Material and Methods. Eighty-three patients referred for liver MRI between August 2017 and January 2020 with a suspected liver tumour were prospectively examined with the standard liver MRI protocol extended by DWI-IVIM sequence. Patients with no liver lesions, haemangiomas, or metastases were excluded. The final study population consisted of 34 patients with primary solid liver masses, 9 with FNH, 4 with regenerative nodules, 10 with HCC, and 11 with CCC. The PDR coefficient was introduced, defined as the ratio of the rate of signal S(b) decrease due to the IVIM effect to the rate of signal S(b) decrease due to the diffusion process, for b = 0.

Results

No significant differences were found between benign and malignant lesions in the case of IVIM parameters (f, D, or D^∗) and ADC. Significant differences were observed only for PDR, with lower values for malignant lesions (p = 0.03). The ROC analysis yielded an AUC value for PDR equal to 0.74, with a cut-off value of 5.06, sensitivity of 81%, specificity of 77%, and accuracy of 79%.

Conclusion

PDR proved to be more effective than IVIM parameters and ADC in the differentiation of solid benign and malignant primary liver lesions.

Associations between angiogenic factors and intravoxel incoherent motion-derived parameters in diffusion-weighted magnetic resonance imaging of breast cancer

Intravoxel incoherent motion (IVIM) diffusion-weighted magnetic resonance imaging (MRI) can be used to estimate perfusion-related parameters, but these parameters may differ, based on the curve-fitting algorithm used for IVIM. Microvessel density (MVD) and vascular endothelial growth factor (VEGF) status are used as angiogenic factors in breast cancer. We aimed to investigate the relationship between MVD, VEGF, and intravoxel incoherent motion (IVIM)-derived parameters, obtained by 4 curve-fitting algorithms, in patients with invasive breast cancers.This retrospective study investigated IVIM-derived parameters, D (ie, tissue diffusivity), D∗ (ie, pseudodiffusivity), and f (ie, perfusion fraction), of 55 breast cancers, using 10 b values (range, 0-800 s/mm2) and 4 curve-fitting algorithms: algorithm 1, linear fitting of D and f first, followed by D∗; algorithm 2, linear fitting of D and f and nonlinear fitting of D∗; algorithm 3, linear fitting of D and f, linear fitting of D∗, and ignoring D contribution for low b values; and algorithm 4, full nonlinear fitting of D, f, and D∗. We evaluated whole-tumor histograms of D, f, and D∗ for their association with MVD and VEGF.D∗10, D∗25, D∗50, D∗mean, D∗75, D∗90, f10, and f25, derived using algorithm 3, were associated with VEGF expression (P = .043, P = 0.012, P = .019, P = .024, P = .044, P = .041, P = .010, and P = .005, respectively). However, no correlation existed between MVD and IVIM-derived parameters.Perfusion-related IVIM parameters obtained by curve-fitting algorithm 3 may reflect VEGF expression.

Intravoxel Incoherent Motion Model in Differentiating the Pathological Grades of Esophageal Carcinoma: Comparison of Mono-Exponential and Bi-Exponential Fit Model

PURPOSE

To compare the diagnostic efficiency of the mono-exponential model and bi-exponential model deriving from intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) in differentiating the pathological grade of esophageal squamous cell carcinoma (ESCC).

METHODS

Fifty-four patients with ESCC were divided into three groups of poorly-differentiated (PD), moderately-differentiated (MD), and well-differentiated (WD), and underwent the IVIM-DWI scan. Mono-exponential (Dmono, D*mono, and fmono) and bi-exponential fit parameters (Dbi, D*bi, and fbi) were calculated using the IVIM data for the tumors. Mean parameter values of three groups were compared using a one-way ANOVA followed by post hoc tests. The receiver operating characteristic curve was drawn for differentiating pathological grade of ESCC. Correlations between pathological grades and IVIM parameters were analyzed.

RESULTS

There were significant differences in fmono and fbi among the PD, MD and WD ESCC groups (all p<0.05). The fmono were 0.32 ± 0.07, 0.23 ± 0.08, and 0.16 ± 0.05, respectively, and the fbi were 0.35 ± 0.08, 0.26 ± 0.10, and 0.18 ± 0.07, respectively. There was a significant difference in the Dmono between the WD and the PD group (1.48 ± 0.51* 10-3 mm2/s versus 1.05 ± 0.44*10-3 mm2/s, p<0.05), but there was no significant difference between the WD and MD groups, MD and PD groups (all p>0.05). The D*mono, Dbi, and D*bi showed no significant difference among the three groups (all p>0.05). The area under the curve (AUC) of Dmono, fmono and fbi in differentiating WD from PD ESCC were 0.764, 0.961 and 0.932, and the sensitivity and specificity were 92.9% and 60%, 92.9% and 90%, 85.7% and 100%, respectively. The AUC of fmono and fbi in differentiating MD from PD ESCC were 0.839 and 0.757, and the sensitivity and specificity were 78.6% and 80%, 85.7% and 70%, respectively. The AUC of fmono and fbi in differentiating MD from WD ESCC were 0.746 and 0.740, and the sensitivity and specificity were 65% and 85%, 80% and 60%, respectively. The pathologically differentiated grade was correlated with all IVIM parameters (all p<0.05).

CONCLUSIONS

The mono-exponential IVIM model is superior to the bi-exponential IVIM model in differentiating pathological grades of ESCC, which may be a promising imaging method to predict pathological grades of ESCC.

Age and gender dependence of liver diffusion parameters and the possibility that intravoxel incoherent motion modeling of the perfusion component is constrained by the diffusion component

The aim of this study was to establish reference values for middle-aged subjects and to investigate the age and gender dependence of liver diffusion MRI parameters. The IVIM type of liver diffusion scan was based on a single-shot spin-echo-type echo-planar sequence using a 1.5 T magnet with 16 b-values. Diffusion-derived vessel density (DDVD)(b0b2) or DDVD(b0b10) was the signal difference between b = 0 and b = 2 (or b = 10) s/mm² images after removing visible vessels. IVIM analysis was performed with full fitting and segmented fitting, and with a threshold b-value of 60 or 200 s/mm² , and fitting started from b = 2 s/mm² . Thirty-one men (age range: 25-71 years) and 26 men (age: 22-69 years) had DDVD and IVIM analysis, respectively, while 37 women (age: 20-71 years) and 36 women (age: 20-71 years) had DDVD and IVIM analysis, respectively. DDVD results showed a significant age-related reduction for women. IVIM results for full fitting showed excellent agreement with those for segmented fitting using a threshold b of 60 s/mm² , but this was less good for results with a threshold b of 200 s/mm² . As age increased, female subjects' D_slow measure showed a significant reduction, while their PF and D_fast measures showed a significant increase. For the age group of 40-55 years, DDVD(b0b2), DDVD(b0b10), D_slow , PF and D_fast were 12.26 ± 3.90 au/pixel, 16.95 ± 5.45 au/pixel, 1.072 ± 0.067 (10^-3 mm² /s), 0.141 ± 0.025 and 61.0 ± 14.0 (10^-3 mm² /s) for men, and 13.35 ± 3.6 au/pixel, 17.20 ± 3.62 au/pixel, 1.069 ± 0.074 (10^-3 mm² /s), 0.119 ± 0.014 and 57.1 ± 13.2 (10^-3 mm² /s) for women, respectively. DDVD measure of this study suggest that aging is associated with a reduction in liver perfusion. There is a possibility that a lower D_slow measure is associated with artificially higher PF and D_fast measures, and that IVIM modeling of the perfusion component is constrained by the diffusion component.

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.

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.

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.

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.

A novel bayesian approach with conditional autoregressive specification for intravoxel incoherent motion diffusion-weighted MRI

The Intra-Voxel Incoherent Motion (IVIM) model is largely adopted to estimate slow and fast diffusion coefficients of water molecules in biological tissues, which are used in cancer applications. The most reported fitting approach is a voxel-wise segmented non-linear least square, whereas Bayesian approaches with a direct fit, also considering spatial regularization, were proposed too. In this work a novel segmented Bayesian method was proposed, also in combination with a spatial regularization through a Conditional Autoregressive (CAR) prior specification. The two segmented Bayesian approaches, with and without CAR specification, were compared with two standard least-square and a direct Bayesian fitting methods. All approaches were tested on simulated images and real data of patients with head-and-neck and rectal cancer. Estimation accuracy and maps noisiness were quantified on simulated images, whereas the coefficient of variation and the goodness of fit were evaluated for real data. Both versions of the segmented Bayesian approach outperformed the standard methods on simulated images for pseudo-diffusion (D^∗ ) and perfusion fraction (f), whilst the segmented least-square fitting remained the less biased for the diffusion coefficient (D). On real data, Bayesian approaches provided the less noisy maps, and the two Bayesian methods without CAR generally estimated lower values for f and D^∗ coefficients with respect to the other approaches. The proposed segmented Bayesian approaches were superior, in terms of estimation accuracy and maps quality, to the direct Bayesian model and the least-square fittings. The CAR method improved the estimation accuracy, especially for D^∗ .

Optimal acquisition scheme for flow-compensated intravoxel incoherent motion diffusion-weighted imaging in the abdomen: An accurate and precise clinically feasible protocol

PURPOSE

Flow-compensated (FC) diffusion-weighted MRI (DWI) for intravoxel-incoherent motion (IVIM) modeling allows for a more detailed description of tissue microvasculature than conventional IVIM. The long acquisition time of current FC-IVIM protocols, however, has prohibited clinical application. Therefore, we developed an optimized abdominal FC-IVIM acquisition with a clinically feasible scan time.

METHODS

Precision and accuracy of the FC-IVIM parameters were assessed by fitting the FC-IVIM model to signal decay curves, simulated for different acquisition schemes. Diffusion-weighted acquisitions were added subsequently to the protocol, where we chose the combination of b-value, diffusion time and gradient profile (FC or bipolar) that resulted in the largest improvement to its accuracy and precision. The resulting two optimized FC-IVIM protocols with 25 and 50 acquisitions (FC-IVIMopt25 and FC-IVIMopt50 ), together with a complementary acquisition consisting of 50 diffusion-weighting (FC-IVIMcomp ), were acquired in repeated abdominal free-breathing FC-IVIM imaging of seven healthy volunteers. Intersession and intrasession within-subject coefficient of variation of the FC-IVIM parameters were compared for the liver, spleen, and kidneys.

RESULTS

Simulations showed that the performance of FC-IVIM improved in tissue with larger perfusion fraction and signal-to-noise ratio. The scan time of the FC-IVIMopt25 and FC-IVIMopt50 protocols were 8 and 16 min. The best in vivo performance was seen in FC-IVIMopt50 . The intersession within-subject coefficients of variation of FC-IVIMopt50 were 11.6%, 16.3%, 65.5%, and 36.0% for FC-IVIM model parameters diffusivity, perfusion fraction, characteristic time and blood flow velocity, respectively.

CONCLUSIONS

We have optimized the FC-IVIM protocol, allowing for clinically feasible scan times (8-16 min).

Intravoxel incoherent motion imaging for diagnosing and staging the liver fibrosis and inflammation

PURPOSE

To evaluate the diagnostic accuracy of intravoxel incoherent motion (IVIM) model parameters for the diagnosis and staging of liver fibrosis and inflammation in patients with chronic hepatitis B.

METHODS

Fifty-four patients with chronic hepatitis B and 42 healthy volunteers were included in the study. All subjects were examined by 3 T magnetic resonance imaging. Diffusion-weighted imaging was undertaken with sixteen b values. IVIM parameters [D (true diffusion coefficient), D* (pseudo-diffusion coefficient), f (perfusion fraction)] were calculated. Histological evaluation of biopsy samples was considered the reference standard for the staging of liver fibrosis and inflammation. Differences in IVIM parameters between patient and control groups were analyzed. In the patient group, fibrosis stage and inflammation grade groups were analyzed with respect to IVIM parameters. The correlation was assessed between IVIM parameters and Ishak-modified scale of fibrosis stages and inflammation grades.

RESULTS

The D was significantly lower in the patient group than the control group, p = 0.038 with Cohen's d effect size of 0.452. D was significantly different between fibrosis stage levels. D values decreased in fibrosis stages from the minimal to moderate to marked fibrosis. Fibrosis grades significantly negatively correlated with D and D* values, p = 0.001, and 0.021, respectively. In addition, inflammation grades negatively correlated with f values, p = 0.047.

CONCLUSION

D values measured with IVIM imaging may help to diagnose liver fibrosis. IVIM imaging could be an alternative to liver biopsy for the staging of liver fibrosis.

Comparison of tri-exponential decay versus bi-exponential decay and full fitting versus segmented fitting for modeling liver intravoxel incoherent motion diffusion MRI

OBJECTIVES

To determine whether bi- or tri-exponential models, and full or segmented fittings, better fit the intravoxel incoherent motion (IVIM) imaging signal of healthy livers.

METHODS

Diffusion-weighted images were acquired with a 3 T scanner using a respiratory-triggered echo-planar sequence and 16 b-values (0-800 s/mm² ). Eighteen healthy volunteers had their livers scanned twice in the same session, and then once in another session. Liver parenchyma region-of-interest-based measurements were processed with bi-exponential and tri-exponential models, with both full fitting and segmented fitting (threshold b-value = 200 s/mm² ).

RESULTS

With the signal of all scans averaged, bi-exponential model full fitting showed D_slow  = 1.14 × 10^-3  mm² /s, D_fast  = 193.6 × 10^-3  mm² /s, and perfusion fraction (PF) = 16.9%, and segmented fitting showed D_slow  = 0.98 × 10^-3  mm² /s, D_fast  = 42.2 × 10^-3  mm² /s, and PF = 23.3%. IVIM parameters derived from the tri-exponential model were similar for full fitting and segmented fitting, with slow (D'_slow  = 0.98 × 10^-3  mm² /s; F'_slow  = 76.4 or 76.6%), fast (D'_fast  = 15.1 or 15.4 × 10^-3  mm² /s; F'_fast  = 11.8 or 11.7%) and very fast (D'_Vfast  = 445.0 or 448.8 × 10^-3  mm² /s; F'_Vfast  = 11.8 or 11.7%) diffusion compartments. The tri-exponential model provided an overall better fit than the bi-exponential model. For the bi-exponential model, full fitting provided a better fit at very low and low b-values compared with segmented fitting, with the latter tending to underestimate D_fast ; however, the segmented method demonstrated lower error in signal prediction for high b-values. Compared with full fitting, tri-exponential segmented fitting offered better scan-rescan reproducibility.

CONCLUSION

For healthy liver, tri-exponential modeling is preferred to bi-exponential modeling. For the bi-exponential model, segmented fitting underestimates D_fast , but offers a more accurate estimation of D_slow .

Spontaneous rupture of solid pseudopapillary tumor of pancreas: A case report and review of literature

INTRODUCTION

Solid pseudopapillary tumors (SPT) account for 1% to 3% of all pancreatic tumors. They have low malignant potential with a favorable prognosis, and predominantly occur in young women. The pathogenesis and clinical behavior of SPT are still uncertain. In addition, most ruptures of SPT were associated with blunt abdominal trauma, while spontaneous ruptures seemed to be quite rare. Up to now, there have been only 3 spontaneous ruptured SPT cases reported worldwide.

PATIENT CONCERNS

Here, we reported a 22-year-old female patient with left lower abdominal pain. Computed tomography (CT) showed that a hemorrhagic complex solid cystic mass located in the lesser omentum sac.

DIAGNOSIS

According to pathological findings of tumor specimen, the diagnosis of solid pseudopapillary tumor (SPT) of the pancreas was made.

INTERVENTIONS

Distal pancreatectomy and splenectomy was carried out.

OUTCOMES

The patient recovered to normal status within 10 days after surgery.

CONCLUSION

Besides, we reviewed about 50 cases in literatures to find out the clinical characteristics and differential diagnostic strategies of SPT.

Topics on quantitative liver magnetic resonance imaging

Liver magnetic resonance imaging (MRI) is subject to continuous technical innovations through advances in hardware, sequence and novel contrast agent development. In order to utilize the abilities of liver MR to its full extent and perform high-quality efficient exams, it is mandatory to use the best imaging protocol, to minimize artifacts and to select the most adequate type of contrast agent. In this article, we review the routine clinical MR techniques applied currently and some latest developments of liver imaging techniques to help radiologists and technologists to better understand how to choose and optimize liver MRI protocols that can be used in clinical practice. This article covers topics on (I) fat signal suppression; (II) diffusion weighted imaging (DWI) and intravoxel incoherent motion (IVIM) analysis; (III) dynamic contrast-enhanced (DCE) MR imaging; (IV) liver fat quantification; (V) liver iron quantification; and (VI) scan speed acceleration.

Impact of b-Value Sampling Scheme on Brain IVIM Parameter Estimation in Healthy Subjects

Purpose:

Intravoxel incoherent motion (IVIM) analysis has attracted the interest of the clinical community due to its close relationship with microperfusion. Nevertheless, there is no clear reference protocol for its implementation; one of the questions being which b-value distribution to use. This study aimed to stress the importance of the sampling scheme and to show that an optimized b-value distribution decreases the variance associated with IVIM parameters in the brain with respect to a regular distribution in healthy volunteers.

Methods:

Ten volunteers were included in this study; images were acquired on a 1.5T MR scanner. Two distributions of 16 b-values were used: one considered ‘regular’ due to its close association with that used in other studies, and the other considered ‘optimized’ according to previous studies. IVIM parameters were adjusted according to the bi-exponential model, using two-step method. Analysis was undertaken in ROI defined using in the Automated Anatomical Labeling atlas, and parameters distributions were compared in a total of 832 ROI.

Results:

Maps with fewer speckles were obtained with the ‘optimized’ distribution. Coefficients of variation did not change significantly for the estimation of the diffusion coefficient D but decreased by approximately 39% for the pseudo-diffusion coefficient estimation and by 21% for the perfusion fraction. Distributions of adjusted parameters were found significantly different in 50% of the cases for the perfusion fraction, in 80% of the cases for the pseudo-diffusion coefficient and 17% of the cases for the diffusion coefficient. Observations across brain areas show that the range of average values for IVIM parameters is smaller in the ‘optimized’ case.

Conclusion:

Using an optimized distribution, data are sampled in a way that the IVIM signal decay is better described and less variance is obtained in the fitted parameters. The increased precision gained could help to detect small variations in IVIM parameters.

The Dmipy Toolbox: Diffusion MRI Multi-Compartment Modeling and Microstructure Recovery Made Easy

Non-invasive estimation of brain microstructure features using diffusion MRI (dMRI)—known as Microstructure Imaging—has become an increasingly diverse and complicated field over the last decades. Multi-compartment (MC)-models, representing the measured diffusion signal as a linear combination of signal models of distinct tissue types, have been developed in many forms to estimate these features. However, a generalized implementation of MC-modeling as a whole, providing deeper insights in its capabilities, remains missing. To address this fact, we present Diffusion Microstructure Imaging in Python (Dmipy), an open-source toolbox implementing PGSE-based MC-modeling in its most general form. Dmipy allows on-the-fly implementation, signal modeling, and optimization of any user-defined MC-model, for any PGSE acquisition scheme. Dmipy follows a “building block”-based philosophy to Microstructure Imaging, meaning MC-models are modularly constructed to include any number and type of tissue models, allowing simultaneous representation of a tissue's diffusivity, orientation, volume fractions, axon orientation dispersion, and axon diameter distribution. In particular, Dmipy is geared toward facilitating reproducible, reliable MC-modeling pipelines, often allowing the whole process from model construction to parameter map recovery in fewer than 10 lines of code. To demonstrate Dmipy's ease of use and potential, we implement a wide range of well-known MC-models, including IVIM, AxCaliber, NODDI(x), Bingham-NODDI, the spherical mean-based SMT and MC-MDI, and spherical convolution-based single- and multi-tissue CSD. By allowing parameter cascading between MC-models, Dmipy also facilitates implementation of advanced approaches like CSD with voxel-varying kernels and single-shell 3-tissue CSD. By providing a well-tested, user-friendly toolbox that simplifies the interaction with the otherwise complicated field of dMRI-based Microstructure Imaging, Dmipy contributes to more reproducible, high-quality research.

Quantitative pancreatic MRI: a pathology-based review

MRI plays an important role in the clinical management of pancreatic disorders and interpretation is reliant on qualitative assessment of anatomy. Conventional sequences capturing pancreatic structure can however be adapted to yield quantitative measures which provide more diagnostic information, with a view to increasing diagnostic accuracy, improving patient stratification, providing robust non-invasive outcome measures for therapeutic trials and ultimately personalizing patient care. In this review, we evaluate the use of established techniques such as secretin-enhanced MR cholangiopancreatography, diffusion-weighted imaging, T ₁, T ₂* and fat fraction mapping, but also more experimental methods such as MR elastography and arterial spin labelling, and their application to the assessment of diffuse pancreatic disease (including chronic, acute and autoimmune pancreatitis/IgG4 disease, metabolic disease and iron deposition disorders) and cystic/solid focal pancreatic masses. Finally, we explore some of the broader challenges to their implementation and future directions in this promising area.

Optimization of b-value schemes for estimation of the diffusion coefficient and the perfusion fraction with segmented intravoxel incoherent motion model fitting

PURPOSE

Intravoxel incoherent motion (IVIM) modeling for estimation of the diffusion coefficient (D) and perfusion fraction (f) is increasingly popular, but no consensus on standard protocols exists. This study provides a framework for optimization of b-value schemes for reduced estimation uncertainty of D and f from segmented model fitting.

THEORY

Analytical expressions for uncertainties of D and f from segmented model fitting were derived as Cramer-Rao lower bounds (CRLBs).

METHODS

Optimized b-value schemes were obtained for 3 to 12 acquisitions and in the limit of infinitely many acquisitions through constrained minimization of the CRLBs, with b-values constrained to be 0 or 200 to 800 s/mm² . The optimized b-value scheme with eight acquisitions was compared with b-values linearly distributed in the allowed range using simulations and in vivo liver data from seven healthy volunteers.

RESULTS

All optimized b-value schemes contained exactly three unique b-values regardless of the total number of acquisitions (0, 200, and 800 s/mm² ) with repeated acquisitions distributed approximately as 1:2:2. Compared with linearly distributed b-values, the variability of estimates of D and f was reduced by approximately 30% as seen both in simulations and in repeated in vivo measurements.

CONCLUSION

The uncertainty of IVIM D and f estimates can be reduced by the use of optimized b-value schemes.

Whole-lesion ADC histogram analysis is not able to reflect microvessel density in HNSCC

Diffusion-weighted imaging (DWI) is a functional imaging technique sensitive to microstructure in tissues. It is widely acknowledged to reflect cellularity in tumors. A small part of DWI is also sensitive to perfusion-related information and might therefore be also be able to reflect microvessel density in tumor tissues. Aim of the present study was to elucidate possible correlations between microvessel density and apparent diffusion coefficient (ADC) values in head and neck squamous cell carcinoma (HNSCC).Thirty-four patients with histologically proven primary HNSCC were included in the study. DWI was performed with a 3 T magnetic resonance imaging (MRI) (b-values 0 and 800 s/mm) and histogram analysis was calculated with a whole lesion measurement. In every case, microvessel density was estimated with CD105-stained specimens.There were no statistically significant correlations between ADC histogram parameters and microvessel density. The calculated correlation coefficients ranged from r = -0.27, P = .13 for entropy and vessel area to r = 0.16, P = .40 for ADCmin and vessel count.Whole-lesion histogram analysis of ADC values cannot reflect microvessel density in HNSCC.

Intravoxel incoherent motion derived liver perfusion/diffusion readouts can be reliable biomarker for the detection of viral hepatitis B induced liver fibrosis

BACKGROUND

Recent two studies reported that intravoxel incoherent motion (IVIM) analysis can separate healthy livers and viral hepatitis B (VHB) induced liver fibrosis. However, in these two studies the starting b value for bi-exponential decay analysis was b =10 and 15 s/mm2 respectively. The current study has two primary aims. The first is to further confirm the diagnostic value of IVIM in detecting liver fibrosis. The second is to test whether by sampling very low b value densely, then b =0 s/mm2 image could be included to improve IVIM's diagnostic performance.

METHODS

This was a prospective study with data acquired at the Third Xiangya Hospital of Central South University, Changsha, China. Healthy volunteers and patients suspected of VHB induced liver fibrosis with liver biopsy performed, as well as hepatocellular carcinoma patients scheduled for surgery, were recruited. All the hepatocellular carcinoma patients had liver fibrosis. After exclusions based on pre-defined criteria for image data quality, for IVIM analysis this study included 20 healthy volunteers; 4 chronic VHB patients with biopsy showing no liver fibrosis; 11 stage-1 liver fibrosis patients, 10 stage-2 liver fibrosis patients, 2 stage-3 liver fibrosis patients, and 5 stage-4 liver fibrosis patients. In the liver fibrosis patients, 1, 19, and 8 cases had inflammation grade-0, grade-1, and grade-2 respectively. The reference IVIM bi-exponential decay curve fitting analysis was segmented fitting performed with b =2 s/mm2 image as the starting point and a threshold-b of 60 s/mm2. This reference fitting method was compared with threshold-b of 40 s/mm2, full fitting, fitting starting from b =0, 5, and 10 s/mm2 respectively. The potential correlation between IVIM readouts and liver function was assessed for the liver fibrosis patients.

RESULTS

Based on the smaller coefficient of variation (CoV) for the volunteer group and the smaller patient/volunteer ratios [= (mean measurement for patient groups)/(mean measurement for healthy volunteers)], the comparison of fitting methods favored the reference approach starting from b =2 s/mm2 with a threshold-b of 60 s/mm2. The IVIM measures of four patients without liver fibrosis resembled those of healthy subjects. PF offered the best diagnostic value for separating healthy livers and fibrotic livers, and a threshold of PF =0.1406 separated all fibrotic livers and healthy livers with an exception of one hepatocellular carcinoma patient (fibrosis grade-2/inflammation grade-2). The correlation between fibrosis grading and inflammation grading was weakly positive; while compared with fibrotic livers with inflammation grade-1, fibrotic livers with inflammation grade-2 showed a trend of higher Dfast. A weak correlation is shown with lower PF and lower Dfast associated with lower total protein, lower albumin; higher alanine transaminase, higher aspartate transaminase; higher total bilirubin, and higher direct bilirubin.

CONCLUSIONS

Segmented-fitting with threshold-b =60 s/mm2 and starting from non-zero very low b value outperforms other methods. IVIM has high sensitivity in detecting liver fibrosis, and PF and Dfast have potential correlation with serum liver function biomarkers. IVIM measures and liver fibrosis grading are not in a linear relationship.

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.

Dependence of intravoxel incoherent motion diffusion MR threshold b-value selection for separating perfusion and diffusion compartments and liver fibrosis diagnostic performance

BACKGROUND

Intravoxel incoherent motion (IVIM) tissue parameters depend on the threshold b-value.

PURPOSE

To explore how threshold b-value impacts PF ( f), D_slow ( D), and D_fast ( D*) values and their performance for liver fibrosis detection.

MATERIAL AND METHODS

Fifteen healthy volunteers and 33 hepatitis B patients were included. With a 1.5-T magnetic resonance (MR) scanner and respiration gating, IVIM data were acquired with ten b-values of 10, 20, 40, 60, 80, 100, 150, 200, 400, and 800 s/mm². Signal measurement was performed on the right liver. Segmented-unconstrained analysis was used to compute IVIM parameters and six threshold b-values in the range of 40-200 s/mm² were compared. PF, D_slow, and D_fast values were placed along the x-axis, y-axis, and z-axis, and a plane was defined to separate volunteers from patients.

RESULTS

Higher threshold b-values were associated with higher PF measurement; while lower threshold b-values led to higher D_slow and D_fast measurements. The dependence of PF, D_slow, and D_fast on threshold b-value differed between healthy livers and fibrotic livers; with the healthy livers showing a higher dependence. Threshold b-value = 60 s/mm² showed the largest mean distance between healthy liver datapoints vs. fibrotic liver datapoints, and a classification and regression tree showed that a combination of PF (PF < 9.5%), D_slow (D_slow < 1.239 × 10^-3 mm²/s), and D_fast (D_fast < 20.85 × 10^-3 mm²/s) differentiated healthy individuals and all individual fibrotic livers with an area under the curve of logistic regression (AUC) of 1.

CONCLUSION

For segmented-unconstrained analysis, the selection of threshold b-value = 60 s/mm² improves IVIM differentiation between healthy livers and fibrotic livers.

Removal of evidential motion-contaminated and poorly fitted image data improves IVIM diffusion MRI parameter scan-rescan reproducibility

Background It has been reported that intravoxel incoherent motion (IVIM) diffusion magnetic resonance imaging (MRI) scan-rescan reproducibility is unsatisfactory. Purpose To study IVIM MRI parameter reproducibility for liver parenchyma after the removal of motion-contaminated and/or poorly fitted image data. Material and Methods Eighteen healthy volunteers had liver scans twice in the same session to assess scan-rescan repeatability, and again in another session after an average interval of 13 days to assess reproducibility. Diffusion-weighted images were acquired with a 3-T scanner using respiratory-triggered echo-planar sequence and 16 b-values (0-800 s/mm²). Measurement was performed on the right liver with segment-unconstrained least square fitting. Image series with evidential anatomical mismatch, apparent artifacts, and poorly fitted signal intensity vs. b-value curve were excluded. A minimum of three slices was deemed necessary for IVIM parameter estimation. Results With a total 54 examinations, six did not satisfy inclusion criteria, leading to a success rate of 89%, and 14 volunteers were finally included for the repeatability/reproducibility study. A total of 3-10 slices per examination (mean = 5.3 slices, median = 5 slices) were utilized for analysis. Using threshold b-value = 80 s/mm², the coefficient of variation and within-subject coefficient of variation for repeatability were 2.86% and 3.36% for D_slow, 3.81% and 4.24% for perfusion fraction (PF), 18.16% and 24.88% for D_fast; and those for reproducibility were 2.48% and 3.24% for D_slow, 4.91% and 5.38% for PF, and 21.18% and 30.89% for D_fast. Conclusion Removal of motion-contaminated and/or poorly fitted image data improves IVIM parameter reproducibility.

Evaluation of Six Diffusion-weighted MRI Models for Assessing Effects of Neoadjuvant Chemoradiation in Pancreatic Cancer Patients

PURPOSE

To compare 6 diffusion-weighted imaging (DWI) MRI models for response evaluation in patients with pancreatic ductal adenocarcinoma (PDAC).

MATERIALS AND METHODS

DWI images were acquired at 3T for b = 0-600 s/mm² in fourteen patients with advanced PDAC during 2 separate pretreatment sessions and 9 patients with (borderline) resectable PDAC pre and post neoadjuvant chemoradiation. Data was fitted with a mono-exponential (ADC), double mono-exponential to b = 0 and 100 s/mm² (ADC_fast), and b = 100 and 600 s/mm² (ADC_slow), IVIM model with D* free (D, f, D*) and fixed (D, f), tri-exponent (D, f₁, f₂), and stretched exponent model (DDC, α). Goodness of fit (adjusted R²), tumor to normal tissue contrast, repeatability (coefficient of variation), and parameter correlations (Spearman's rho) were assessed for the repeated measures. Treatment induced changes were assessed and compared to the repeatability.

RESULTS

The mono-exponential model had the lowest goodness of fit in both tumor (R² = 0.94) and normal-appearing pancreas (R² = 0.88). Tumour to normal tissue contrast was higher for the 'non-diffusion' parameters (ADC_fast, f, D*, f₁, f₂, α), with better repeatability for the diffusion parameters (ADC, ADC_slow, D, DDC). Diffusion parameters were strongly correlated between the models (rho ≥0.81) and showed a general treatment associated increase. All models were able to identify individual treatment effects, showing a change greater than the repeatability in 5 out of 9 patients for at least one of the parameters.

CONCLUSIONS

Individual treatment evaluation is possible with all investigated DWI models, with treatment associated changes exceeding the repeatability. The double monoexponential fit with ADC_fast and ADC_slow is able to discriminate between non-diffusion and diffusion related effects, is measured fast and can be performed on most commercial scanners, making it an attractive alternative for the more advanced multiparametric models in radiotherapy treatment evaluation.

Quantification of hepatic perfusion and hepatocyte function with dynamic gadoxetic acid-enhanced MRI in patients with chronic liver disease

The purpose of the present study was to develop and perform initial validation of dynamic MRI enhanced with gadoxetic acid as hepatobiliary contrast agent to quantify hepatic perfusion and hepatocyte function in patients with chronic liver disease. Free-breathing, dynamic gadoxetic acid-enhanced MRI was performed at 3.0 T using a 3D time-resolved angiography sequence with stochastic trajectories during 38 min. A dual-input three-compartment model was developed to derive hepatic perfusion and hepatocyte function parameters. Method feasibility was assessed in 23 patients with biopsy-proven chronic liver disease. Parameter analysis could be performed in 21 patients (91%). The hepatocyte function parameters were more discriminant than the perfusion parameters to differentiate between patients with minimal fibrosis (METAVIR F0-F1), intermediate fibrosis (F2-F3) and cirrhosis (F4). The areas under the receiver operating characteristic curves (ROCs) to diagnose significant fibrosis (METAVIR F ≥ 2) were: 0.95 (95% CI: 0.87-1; P<0.001) for biliary efflux, 0.88 (95% CI: 0.73-1; P<0.01) for sinusoidal backflux, 0.81 (95% CI: 0.61-1; P<0.05) for hepatocyte uptake fraction and 0.75 (95% CI: 0.54-1; P<0.05) for hepatic perfusion index (HPI), respectively. These initial results in patients with chronic liver diseases show that simultaneous quantification of hepatic perfusion and hepatocyte function is feasible with free breathing dynamic gadoxetic acid-enhanced MRI. Hepatocyte function parameters may be relevant to assess liver fibrosis severity.

Intravoxel incoherent motion diffusion-weighted MR imaging parameters predict pathological classification in thymic epithelial tumors

We evaluated the performance of intravoxel incoherent motion (IVIM) parameters for preoperatively predicting the subtype and Masaoka stage of thymic epithelial tumors (TETs). Seventy-seven patients with pathologically confirmed TETs underwent a diffusion weighted imaging (DWI) sequence with 9 b values. Differences in the slow diffusion coefficient (D), fast perfusion coefficient (D), and perfusion fraction (f) IVIM parameters, as well as the multi b-value fitted apparent diffusion coefficient (ADC_mb), were compared among patients with low-risk (LRT) and high-risk thymomas (HRT) and thymic carcinomas (TC), and between early stage (stages I and II) and advanced stage (stages III and IV) TET patients. ADC_mb, D, and D values were higher in the LRT group than in the HRT or TC group, but did not differ between the HRT and TC groups. The mean ADC_mb, D, and D values were higher in the early stage TETs group than the advanced stage TETs group. The f values did not differ among the groups. These results suggest that IVIM DWI could be used to preoperatively predict subtype and Masaoka stage in TET patients.

Comparison of six fit algorithms for the intra-voxel incoherent motion model of diffusion-weighted magnetic resonance imaging data of pancreatic cancer patients

The intravoxel incoherent motion (IVIM) model for diffusion-weighted imaging (DWI) MRI data bears much promise as a tool for visualizing tumours and monitoring treatment response. To improve the currently poor precision of IVIM, several fit algorithms have been suggested. In this work, we compared the performance of two Bayesian IVIM fit algorithms and four other IVIM fit algorithms for pancreatic cancer imaging. DWI data were acquired in 14 pancreatic cancer patients during two MRI examinations. Three different measures of performance of the fitting algorithms were assessed: (i) uniqueness of fit parameters (Spearman’s rho); (ii) precision (within-subject coefficient of variation, wCV); and (iii) contrast between tumour and normal-appearing pancreatic tissue. For the diffusivity D and perfusion fraction f, a Bayesian fit (IVIM-Bayesian-lin) offered the best trade-off between tumour contrast and precision. With the exception for IVIM-Bayesian-lin, all algorithms resulted in a very poor precision of the pseudo-diffusion coefficient D* with a wCV of more than 50%. The pseudo-diffusion coefficient D* of the Bayesian approaches were, however, significantly correlated with D and f. Therefore, the added value of fitting D* was considered limited in pancreatic cancer patients. The easier implemented least squares fit with fixed D* (IVIM-fixed) performed similar to IVIM-Bayesian-lin for f and D. In conclusion, the best performing IVIM fit algorithm was IVM-Bayesian-lin, but an easier to implement least squares fit with fixed D* performs similarly in pancreatic cancer patients.

Intra-voxel incoherent motion MRI of the living human foetus: technique and test-retest repeatability

Background

Our purpose was to test the within-subject (test–retest) reproducibility of the perfusion fraction, diffusion coefficient, and pseudo-diffusion coefficient measurements in various foetus organs and in the placenta based on the intra-voxel incoherent motion (IVIM) principle.

Methods

In utero diffusion-weighted IVIM magnetic resonance imaging (MRI) was performed in 15 pregnant women (pregnancy age 21–36 weeks) on 1.5-T and 3.0-T clinical scanners with b-factors in the range of 0–900 s/mm² in 16 steps. A bi-exponential model was fitted on the volume-averaged diffusion values. Perfusion fraction (f), diffusion coefficient (d), and pseudo-diffusion coefficient (D*) were calculated. Within-subject reproducibility was evaluated as test–retest variability (VAR %) of the IVIM parameters in the foetal frontal cortex, frontal white matter, cerebellum, lungs, kidneys, liver, and in the placenta.

Results

For the foetal lungs, liver and the placenta, test–retest variability was in the range of 14–20% for f, 12–14% for d, and 17–25% for D*. The diffusion coefficients of the investigated brain regions were moderately to highly reproducible (VAR 5–15%). However, f and D* showed inferior reproducibility compared to corresponding measures for the lungs, liver, and placenta. The IVIM parameters of the foetal kidney were revealed to be highly variable across scans.

Conclusions

IVIM MRI potentially provides a novel method for examining microvascular perfusion and diffusion in the developing human foetus. However, reproducibility of perfusion and diffusion parameters depends greatly upon data quality, foetal and maternal movements, and foetal-specific image post-processing.

Liver fibrosis: in vivo evaluation using intravoxel incoherent motion-derived histogram metrics with histopathologic findings at 3.0 T

PURPOSE

To retrospectively determine the feasibility of intravoxel incoherent motion (IVIM) imaging based on histogram analysis for the staging of liver fibrosis (LF) using histopathologic findings as the reference standard.

METHODS

56 consecutive patients (14 men, 42 women; age range, 15-76, years) with chronic liver diseases (CLDs) were studied using IVIM-DWI with 9 b-values (0, 25, 50, 75, 100, 150, 200, 500, 800 s/mm²) at 3.0 T. Fibrosis stage was evaluated using the METAVIR scoring system. Histogram metrics including mean, standard deviation (Std), skewness, kurtosis, minimum (Min), maximum (Max), range, interquartile (Iq) range, and percentiles (10, 25, 50, 75, 90th) were extracted from apparent diffusion coefficient (ADC), true diffusion coefficient (D), pseudo-diffusion coefficient (D*), and perfusion fraction (f) maps. All histogram metrics among different fibrosis groups were compared using one-way analysis of variance or nonparametric Kruskal-Wallis test. For significant parameters, receivers operating characteristic curve (ROC) analyses were further performed for the staging of LF.

RESULTS

Based on their METAVIR stage, the 56 patients were reclassified into three groups as follows: F0-1 group (n = 25), F2-3 group (n = 21), and F4 group (n = 10). The mean, Iq range, percentiles (50, 75, and 90th) of D* maps between the groups were significant differences (all P < 0.05). Area under the ROC curve (AUC) of the mean, Iq range, 50, 75, and 90th percentile of D* maps for identifying significant LF (≥F2 stage) was 0.901, 0.859, 0.876, 0.943, and 0.886 (all P < 0.0001), respectively; for diagnosing severe fibrosis or cirrhosis (F4), AUC was 0.917, 0.922, 0.943, 0.985, and 0.939 (all P < 0.0001), respectively. The histogram metrics of ADC, D, and f maps demonstrated no significant difference among the groups (all P > 0.05).

CONCLUSION

Histogram analysis of D* map derived from IVIM can be used to stage liver fibrosis in patients with CLDs and provide more quantitative information beyond the mean value.

Diffusion-weighted MRI and intravoxel incoherent motion model for diagnosis of pediatric solid abdominal tumors

Background

Pediatric retroperitoneal tumors in the renal bed are often large and heterogeneous, and their diagnosis based on conventional imaging alone is not possible. More advanced imaging methods, such as diffusion‐weighted (DW) MRI and the use of intravoxel incoherent motion (IVIM), have the potential to provide additional biomarkers that could facilitate their noninvasive diagnosis.

Purpose

To assess the use of an IVIM model for diagnosis of childhood malignant abdominal tumors and discrimination of benign from malignant lesions.

Study Type

Retrospective.

Population

Forty‐two pediatric patients with abdominal lesions (n = 32 malignant, n = 10 benign), verified by histopathology.

Field Strength/Sequence

1.5T MRI system and a DW‐MRI sequence with six b‐values (0, 50, 100, 150, 600, 1000 s/mm²).

Assessment

Parameter maps of apparent diffusion coefficient (ADC), and IVIM maps of slow diffusion coefficient (D), fast diffusion coefficient (D*), and perfusion fraction (f) were computed using a segmented fitting model. Histograms were constructed for whole‐tumor regions of each parameter.

Statistical Tests

Comparison of histogram parameters of and their diagnostic performance was determined using Kruskal–Wallis, Mann–Whitney U, and receiver‐operating characteristic (ROC) analysis.

Results

IVIM parameters D* and f were significantly higher in neuroblastoma compared to Wilms' tumors (P < 0.05). The ROC analysis showed that the best diagnostic performance was achieved with D* 90^th percentile (area under the curve [AUC] = 0.935; P = 0.002; cutoff value = 32,376 × 10⁻⁶ mm²/s) and f mean values (AUC = 1.00; P < 0.001; cutoff value = 14.7) in discriminating between neuroblastoma (n = 11) and Wilms' tumors (n = 8). Discrimination between tumor types was not possible with IVIM D or ADC parameters. Malignant tumors revealed significantly lower ADC, D, and higher D* values than in benign lesions (all P < 0.05).

Data Conclusion

IVIM perfusion parameters could distinguish between malignant childhood tumor types, providing potential imaging biomarkers for their diagnosis.

Level of Evidence: 4

Technical Efficacy: Stage 2

J. Magn. Reson. Imaging 2018;47:1475–1486.

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 brain_500,1000, kidney_300,1000, and liver_200,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.

Impact of prior distributions and central tendency measures on Bayesian intravoxel incoherent motion model fitting

PURPOSE

Bayesian model fitting has been proposed as a robust alternative for intravoxel incoherent motion (IVIM) model-fitting parameter estimation. However, consensus regarding choice of prior distribution and posterior distribution central tendency measure is needed. The aim of this study was to compare the quality of IVIM parameter estimates produced by different prior distributions and central tendency measures, and to gain knowledge about the effect of these choices.

METHODS

Three prior distributions (uniform, reciprocal, and lognormal) and two measures of central tendency (mean and mode) found in the literature were studied using simulations and in vivo data from a tumor mouse model.

RESULTS

Simulations showed that the uniform and lognormal priors were superior to the reciprocal prior, especially for the parameters D and f and clinically relevant SNR levels. The choice of central tendency measure had less effect on the results, but had some effects on estimation bias. Results based on simulations and in vivo data agreed well, indicating high validity of the simulations.

CONCLUSIONS

Choice of prior distribution and central tendency measure affects the results of Bayesian IVIM parameter estimates. This must be considered when comparing results from different studies. The best overall quality of IVIM parameter estimates was obtained using the lognormal prior. Magn Reson Med 79:1674-1683, 2018. © 2017 International Society for Magnetic Resonance in Medicine.