Effect of multiple perfusion components on pseudo-diffusion coefficient in intravoxel incoherent motion imaging

Reproducibility of intravoxel incoherent motion quantification in the liver across field strengths and gradient hardware

PURPOSE

To evaluate reproducibility and interlobar agreement of intravoxel incoherent motion (IVIM) quantification in the liver across field strengths and MR scanners with different gradient hardware.

METHODS

Cramer-Rao lower bound optimization was performed to determine optimized monopolar and motion-robust 2D (b-value and first-order motion moment [M1]) IVIM-DWI acquisitions. Eleven healthy volunteers underwent diffusion MRI of the liver, where each optimized acquisition was obtained five times across three MRI scanners. For each data set, IVIM estimates (diffusion coefficient (D), pseudo-diffusion coefficients (d 1 * $$ {d}_1^{\ast } $$ andd 2 * $$ {d}_2^{\ast } $$ ), blood velocity SDs (Vb1 and Vb2), and perfusion fractions [f1 and f2]) were obtained in the right and left liver lobes using two signal models (pseudo-diffusion and M1-dependent physical) with and without T2 correction (fc1 and fc2) and three fitting techniques (tri-exponential region of interest-based full and segmented fitting and blood velocity SD distribution fitting). Reproducibility and interlobar agreement were compared across methods using within-subject and pairwise coefficients of variation (CVw and CVp), paired sample t-tests, and Bland-Altman analysis.

RESULTS

Using a combination of motion-robust 2D (b-M1) data acquisition, M1-dependent physical signal modeling with T2 correction, and blood velocity SD distribution fitting, multiscanner reproducibility with median CVw = 5.09%, 11.3%, 9.20%, 14.2%, and 12.6% for D, Vb1, Vb2, fc1, and fc2, respectively, and interlobar agreement with CVp = 8.14%, 11.9%, 8.50%, 49.9%, and 42.0%, respectively, was achieved.

CONCLUSION

Recently proposed advanced IVIM acquisition, signal modeling, and fitting techniques may facilitate reproducible IVIM quantification in the liver, as needed for establishment of IVIM-based quantitative biomarkers for detection, staging, and treatment monitoring of diseases.

Image denoising and model-independent parameterization for IVIM MRI

Objective. To improve intravoxel incoherent motion imaging (IVIM) magnetic resonance Imaging quality using a new image denoising technique and model-independent parameterization of the signal versusb-value curve.Approach. IVIM images were acquired for 13 head-and-neck patients prior to radiotherapy. Post-radiotherapy scans were also acquired for five of these patients. Images were denoised prior to parameter fitting using neural blind deconvolution, a method of solving the ill-posed mathematical problem of blind deconvolution using neural networks. The signal decay curve was then quantified in terms of several area under the curve (AUC) parameters. Improvements in image quality were assessed using blind image quality metrics, total variation (TV), and the correlations between parameter changes in parotid glands with radiotherapy dose levels. The validity of blur kernel predictions was assessed by the testing the method's ability to recover artificial 'pseudokernels'. AUC parameters were compared with monoexponential, biexponential, and triexponential model parameters in terms of their correlations with dose, contrast-to-noise (CNR) around parotid glands, and relative importance via principal component analysis.Main results. Image denoising improved blind image quality metrics, smoothed the signal versusb-value curve, and strengthened correlations between IVIM parameters and dose levels. Image TV was reduced and parameter CNRs generally increased following denoising.AUCparameters were more correlated with dose and had higher relative importance than exponential model parameters.Significance. IVIM parameters have high variability in the literature and perfusion-related parameters are difficult to interpret. Describing the signal versusb-value curve with model-independent parameters like theAUCand preprocessing images with denoising techniques could potentially benefit IVIM image parameterization in terms of reproducibility and functional utility.

The Feasibility of Using Tri-Exponential Intra-Voxel Incoherent Motion DWI for Identifying the Microvascular Invasion in Hepatocellular Carcinoma

Purpose

To assess the effectiveness of tri-exponential Intra-Voxel Incoherent Motion (tri-IVIM) MRI in preoperatively identifying microvascular invasion (MVI) in hepatocellular carcinoma (HCC).

Patients and Methods

In this prospective study, 67 patients with HCC were included. Metrics from bi-exponential IVIM (bi-IVIM) and tri-IVIM were calculated. Subgroup comparisons were analyzed using the independent Student's t-test or Mann-Whitney U-test. Logistic regression was performed to explore clinical risk factors. Diagnostic performance was assessed using receiver operating characteristic (ROC) curves, calibration curves and decision curve analysis.

Results

MVI-positive HCCs exhibited significantly lower true diffusion coefficient (Dt) from bi-IVIM, as well as fast-diffusion coefficients (Df) and slow-diffusion coefficients (Ds) from tri-IVIM, compared to MVI-negative HCCs (p < 0.05). Tumor size and alpha-fetoprotein (AFP) were identified as risk factors. The combination of tri-IVIM-derived metrics (Ds and Df) yielded higher diagnostic accuracy (AUC = 0.808) compared to bi-IVIM (AUC = 0.741). A predictive model based on a nomogram was constructed using Ds, Df, tumor size, and AFP, resulting in the highest diagnostic accuracy (AUC = 0.859). Decision curve analysis indicated that the constructed model, provided the highest net benefit by accurately stratifying the risk of MVI, followed by tri-IVIM and bi-IVIM.

Conclusion

Tri-IVIM can provide information on perfusion and diffusion for evaluating MVI in HCC. Additionally, tri-IVIM outperformed bi-IVIM in identifying MVI-positive HCC. By integrating clinical risk factors and metrics from tri-IVIM, a predictive nomogram exhibited the highest diagnostic accuracy, potentially aiding in the noninvasive and preoperative assessment of MVI.

Precision of region of interest-based tri-exponential intravoxel incoherent motion quantification and the role of the Intervoxel spatial distribution of flow velocities

PURPOSE

The purpose of this work was to obtain precise tri-exponential intravoxel incoherent motion (IVIM) quantification in the liver using 2D (b-value and first-order motion moment [M1 ]) IVIM-DWI acquisitions and region of interest (ROI)-based fitting techniques.

METHODS

Diffusion MRI of the liver was performed in 10 healthy volunteers using three IVIM-DWI acquisitions: conventional monopolar, optimized monopolar, and optimized 2D (b-M1 ). For each acquisition, bi-exponential and tri-exponential full, segmented, and over-segmented ROI-based fitting and a newly proposed blood velocity SDdistribution (BVD) fitting technique were performed to obtain IVIM estimates in the right and left liver lobes. Fitting quality was evaluated using corrected Akaike information criterion. Precision metrics (test-retest repeatability, inter-reader reproducibility, and inter-lobar agreement) were evaluated using Bland-Altman analysis, repeatability/reproducibility coefficients (RPCs), and paired sample t-tests. Precision was compared across acquisitions and fitting methods.

RESULTS

High repeatability and reproducibility was observed in the estimations of the diffusion coefficient (Dtri  = [1.03 ± 0.11] × 10-3  mm2 /s; RPCs ≤ 1.34 × 10-4  mm2 /s), perfusion fractions (F1  = 3.19 ± 1.89% and F2  = 16.4 ± 2.07%; RPCs ≤ 2.51%), and blood velocity SDs (Vb,1  = 1.44 ± 0.14 mm/s and Vb,2  = 3.62 ± 0.13 mm/s; RPCs ≤ 0.41 mm/s) in the right liver lobe using the 2D (b-M1 ) acquisition in conjunction with BVD fitting. Using these methods, significantly larger (p < 0.01) estimates of Dtri and F1 were observed in the left lobe in comparison to the right lobe, while estimates of Vb,1 and Vb,2 demonstrated high interlobar agreement (RPCs ≤ 0.45 mm/s).

CONCLUSIONS

The 2D (b-M1 ) IVIM-DWI data acquisition in conjunction with BVD fitting enables highly precise tri-exponential IVIM quantification in the right liver lobe.

Machine Learning Models and Multiparametric Magnetic Resonance Imaging for the Prediction of Pathologic Response to Neoadjuvant Chemotherapy in Breast Cancer

Simple Summary

Achieving pathological complete response (pCR) to neoadjuvant chemotherapy (NAC) in breast cancer (BC) is crucial, as pCR is a surrogate marker for survival. However, only 10–30% of patients achieve it. It is therefore essential to develop tools that enable to accurately predict response. Recently, different studies have demonstrated the feasibility of applying machine learning approaches to non-invasively predict pCR before NAC administration from magnetic resonance imaging (MRI) data. Some of those models are based on radiomics, an emerging field that allows the automated extraction of clinically relevant information from radiologic images. However, the research is still at an early stage and further data are needed. Here, we determine whether the combination of imaging data (radiomic features and perfusion/diffusion imaging biomarkers) extracted from multiparametric MRIs and clinical variables can improve pCR prediction to NAC compared to models only including imaging or clinical data, potentially avoiding unnecessary treatment and delays to surgery.

Abstract

Background: Most breast cancer (BC) patients fail to achieve pathological complete response (pCR) after neoadjuvant chemotherapy (NAC). The aim of this study was to evaluate whether imaging features (perfusion/diffusion imaging biomarkers + radiomic features) extracted from pre-treatment multiparametric (mp)MRIs were able to predict, alone or in combination with clinical data, pCR to NAC. Methods: Patients with stage II-III BC receiving NAC and undergoing breast mpMRI were retrospectively evaluated. Imaging features were extracted from mpMRIs performed before NAC. Three different machine learning models based on imaging features, clinical data or imaging features + clinical data were trained to predict pCR. Confusion matrices and performance metrics were obtained to assess model performance. Statistical analyses were conducted to evaluate differences between responders and non-responders. Results: Fifty-eight patients (median [range] age, 52 [45–58] years) were included, of whom 12 showed pCR. The combined model improved pCR prediction compared to clinical and imaging models, yielding 91.5% of accuracy with no false positive cases and only 17% false negative results. Changes in different parameters between responders and non-responders suggested a possible increase in vascularity and reduced tumour heterogeneity in patients with pCR, with the percentile 25th of time-to-peak (TTP), a classical perfusion parameter, being able to discriminate both groups in a 75% of the cases. Conclusions: A combination of mpMRI-derived imaging features and clinical variables was able to successfully predict pCR to NAC. Specific patient profiles according to tumour vascularity and heterogeneity might explain pCR differences, where TTP could emerge as a putative surrogate marker for pCR.

b value and first-order motion moment optimized data acquisition for repeatable quantitative intravoxel incoherent motion DWI

PURPOSE

To design a b value and first-order motion moment (M1 ) optimized data acquisition for repeatable intravoxel incoherent motion (IVIM) quantification in the liver.

METHODS

Cramer-Rao lower bound optimization was performed to determine optimal monopolar and optimal 2D samplings of the b-M1 space based on noise performance. Monte Carlo simulations were used to evaluate the bias and variability in estimates obtained using the proposed optimal samplings and conventional monopolar sampling. Diffusion MRI of the liver was performed in 10 volunteers using 3 IVIM acquisitions: conventional monopolar, optimized monopolar, and b-M1 -optimized gradient waveforms (designed based on the optimal 2D sampling). IVIM parameter maps of diffusion coefficient, perfusion fraction, and blood velocity SD were obtained using nonlinear least squares fitting. Noise performance (SDs), stability (outlier percentage), and test-retest or scan-rescan repeatability (intraclass correlation coefficients) were evaluated and compared across acquisitions.

RESULTS

Cramer-Rao lower bound and Monte Carlo simulations demonstrated improved noise performance of the optimal 2D sampling in comparison to monopolar samplings. Evaluating the designed b-M1 -optimized waveforms in healthy volunteers, significant decreases (p < 0.05) in the SDs and outlier percentages were observed for measurements of diffusion coefficient, perfusion fraction, and blood velocity SD in comparison to measurements obtained using monopolar samplings. Good-to-excellent repeatability (intraclass correlation coefficients ≥ 0.77) was observed for all 3 parameters in both the right and left liver lobes using the b-M1 -optimized waveforms.

CONCLUSIONS

2D b-M1 -optimized data acquisition enables repeatable IVIM quantification with improved noise performance. 2D acquisitions may advance the establishment of IVIM quantitative biomarkers for liver diseases.

Diffusion-Weighted Magnetic Resonance Imaging in Ovarian Cancer: Exploiting Strengths and Understanding Limitations

Detection, characterization, staging, and response assessment are key steps in the imaging pathway of ovarian cancer. The most common type, high grade serous ovarian cancer, often presents late, so that accurate disease staging and response assessment are required through imaging in order to improve patient management. Currently, computerized tomography (CT) is the most common method for these tasks, but due to its poor soft-tissue contrast, it is unable to quantify early response within lesions before shrinkage is observed by size criteria. Therefore, quantifiable techniques, such as diffusion-weighted magnetic resonance imaging (DW-MRI), which generates high contrast between tumor and healthy tissue, are increasingly being explored. This article discusses the basis of diffusion-weighted contrast and the technical issues that must be addressed in order to achieve optimal implementation and robust quantifiable diffusion-weighted metrics in the abdomen and pelvis. The role of DW-MRI in characterizing adnexal masses in order to distinguish benign from malignant disease, and to differentiate borderline from frankly invasive malignancy is discussed, emphasizing the importance of morphological imaging over diffusion-weighted metrics in this regard. Its key role in disease staging and predicting resectability in comparison to CT is addressed, including its valuable use as a biomarker for following response within individual lesions, where early changes in the apparent diffusion coefficient in peritoneal metastases may be detected. Finally, the task of implementing DW-MRI into clinical trials in order to validate this biomarker for clinical use are discussed, along with the trials that include it within their protocols.

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.

Application of Multiparametric Magnetic Resonance Imaging to Monitor the Early Antitumor Effect of CuS@GOD Nanoparticles in a 4 T1 Breast Cancer Xenograft Model

BACKGROUND

We have developed hybrid nanoparticles (NPs) by co-loading copper sulfide (CuS) NPs and glucose oxidase (GOD) (CuS@GOD NPs) to explore their antitumor properties.

PURPOSE

To investigate the feasibility of using multiparametric magnetic resonance imaging (MRI) including intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) and R₂ * mapping to quantitatively assess the early antitumor effect of CuS@GOD NPs.

STUDY TYPE

Prospective.

ANIMAL MODEL

The orthotopic BALB/c mice 4 T1 breast cancer model. The 4 T1 xenografts in group 1 mice received normal saline, group 2 received CuS@GOD NPs, group 3 received CuS NPs plus laser, and group 4 received CuS@GOD NPs plus laser (n = 28 for each group).

FIELD STRENGTH/SEQUENCE

A 3.0 T/IVIM-DWI MRI single-shot echo-planar imaging, R₂ * mapping spoiled gradient recalled echo (SPGR) sequence, T2-weighted images (T2WI) and T1-weighted images (T1WI) fast spin echo (FSE) sequence.

ASSESSMENT

The IVIM-DWI and R₂ * mapping were performed before and after treatment at 0 hour, 0.5 hour, 1 hour, 2 hours, 4 hours, and 24 hours in four groups and the MRI parameters were obtained. Correlation analysis between the MRI parameters and histological analyses was conducted.

STATISTICAL TESTS

One-way ANOVA, Pearson's correlation analysis, two independent samples t test, intraclass correlation coefficient. P < 0.05 was considered to be statistically significant.

RESULTS

In group 4, the tumoral D value was significantly higher than that of group 2 at 24 hours (0.541 ± 0.065 vs. 0.492 ± 0.051). The f value of group 4 was significantly lower than that of groups 1 and 2 at 2 hours (10.83 ± 2.16 vs. 14.28 ± 1.86, 16.67 ± 3.53, respectively). The R₂ * value was significantly increased at 0 hour in group 4 compared to that of groups 1 and 2 (64.552 ± 4.663 vs. 42.441 ± 1.516, 43.165 ± 1.709, respectively). D, f, and R₂ * were correlated with the histological staining results (r = 0.695-0.970).

DATA CONCLUSION

The IVIM-DWI-derived D and f and R₂ * mapping-derived R₂ * could monitor early response to CuS@GOD NPs treatment in vivo.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 2.

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

BACKGROUND

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

PURPOSE

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

STUDY TYPE

Prospective.

POPULATION

Ten healthy subjects without history of heart diseases.

FIELD STRENGTH/SEQUENCE

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

ASSESSMENT

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

STATISTICAL TESTS

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

RESULTS

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

DATA CONCLUSION

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

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 1.

Crossed cerebellar diaschisis after acute ischemic stroke detected by intravoxel incoherent motion magnetic resonance imaging

PURPOSE

To study the value of 3.0 T magnetic resonance imaging with intravoxel incoherent motion (IVIM) in the diagnosis of the crossed cerebellar diaschisis (CCD) after the unilateral supratentorial acute ischemic stroke.

METHODS

Seventy-four patients with acute ischemic stroke who underwent intravoxel incoherent motion (IVIM), arterial spin labeling (ASL), and conventional magnetic resonance imaging (MRI) scanning were enrolled. Intravoxel incoherent motion-derived perfusion-related parameters including fast diffusion coefficient (D*), slow diffusion coefficient (D), vascular volume fraction (f), and arterial spin-labeling-derived cerebral blood flow (CBF) of bilateral cerebellum were measured.

RESULTS

In the CCD-positive group, D*, D, and CBF values of the contralateral cerebellum decreased compared with those of the ipsilesional cerebellum (P < 0.05), whereas f significantly increased (P < 0.05). A positive correlation was detected between the slow diffusion coefficient-based asymmetry index (AI-D) and the cerebral blood flow-based asymmetry index (AI-CBF) (r = 0.515, P < 0.01), whereas the vascular volume fraction-based asymmetry index (AI-f) had a negative correlation with the cerebral blood flow-based asymmetry index (AI-CBF) (r =  - 0.485, P < 0.01). Furthermore, the area under the receiver operating characteristic (ROC) curve value of AI-D and AI-f was 0.81 and 0.76, respectively.

CONCLUSIONS

The IVIM is feasible for the detection of CCD. This technique might provide opportunities to further investigate the pathophysiology of CCD.

Continuous diffusion spectrum computation for diffusion-weighted magnetic resonance imaging of the kidney tubule system

Background

The use of rigid multi-exponential models (with a priori predefined numbers of components) is common practice for diffusion-weighted MRI (DWI) analysis of the kidney. This approach may not accurately reflect renal microstructure, as the data are forced to conform to the a priori assumptions of simplified models. This work examines the feasibility of less constrained, data-driven non-negative least squares (NNLS) continuum modelling for DWI of the kidney tubule system in simulations that include emulations of pathophysiological conditions.

Methods

Non-linear least squares (LS) fitting was used as reference for the simulations. For performance assessment, a threshold of 5% or 10% for the mean absolute percentage error (MAPE) of NNLS and LS results was used. As ground truth, a tri-exponential model using defined volume fractions and diffusion coefficients for each renal compartment (tubule system: D_tubules , f_tubules ; renal tissue: D_tissue , f_tissue ; renal blood: D_blood , f_blood ;) was applied. The impact of: (I) signal-to-noise ratio (SNR) =40-1,000, (II) number of b-values (n=10-50), (III) diffusion weighting (b-range_small =0-800 up to b-range_large =0-2,180 s/mm²), and (IV) fixation of the diffusion coefficients D_tissue and D_blood was examined. NNLS was evaluated for baseline and pathophysiological conditions, namely increased tubular volume fraction (ITV) and renal fibrosis (10%: grade I, mild) and 30% (grade II, moderate).

Results

NNLS showed the same high degree of reliability as the non-linear LS. MAPE of the tubular volume fraction (f_tubules ) decreased with increasing SNR. Increasing the number of b-values was beneficial for f_tubules precision. Using the b-range_large led to a decrease in MAPE _ftubules compared to b-range_small. The use of a medium b-value range of b=0-1,380 s/mm² improved f_tubules precision, and further b_max increases beyond this range yielded diminishing improvements. Fixing D_blood and D_tissue significantly reduced MAPE _ftubules and provided near perfect distinction between baseline and ITV conditions. Without constraining the number of renal compartments in advance, NNLS was able to detect the (fourth) fibrotic compartment, to differentiate it from the other three diffusion components, and to distinguish between 10% vs. 30% fibrosis.

Conclusions

This work demonstrates the feasibility of NNLS modelling for DWI of the kidney tubule system and shows its potential for examining diffusion compartments associated with renal pathophysiology including ITV fraction and different degrees of fibrosis.

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.

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.

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.

An optimized b-value distribution for triexponential intravoxel incoherent motion (IVIM) in the liver

PURPOSE

To find an optimized b-value distribution for reproducible triexponential intravoxel incoherent motion (IVIM) exams in the liver.

METHODS

A numeric optimization of b-value distributions was performed using the triexponential IVIM equation and 27 different IVIM parameter sets. Starting with an initially optimized distribution of 6 b-values, the number of b-values was increased stepwise. Each new b-value was chosen from a set of 64 predefined b-values based on the computed summed relative mean error of the fitted triexponential IVIM parameters. This process was repeated for up to 100 b-values. In simulations and in vivo measurements, optimized b-value distributions were compared to 4 representative distributions found in literature.

RESULTS

The first 16 optimized b-values were 0, 0.3, 0.3, 70, 200, 800, 70, 1, 3.5, 5, 70, 1.2, 6, 45, 1.5, and 60 in units of s/mm² . Low b-values were much more frequent than high b-values. The optimized b-value distribution resulted in a higher fit stability compared to distributions used in literature in both, simulation and in vivo measurements. Using more than 6 b-values, ideally 16 or more, increased the fit stability considerably.

CONCLUSION

Using optimized b-values, the fit uncertainty in triexponential IVIM can be largely reduced. Ideally, 16 or more b-values should be acquired.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Multiparametric Functional Magnetic Resonance Imaging for Evaluating Renal Allograft Injury

Kidney transplantation is the treatment of choice for patients with end-stage renal disease, as it extends survival and increases quality of life in these patients. However, chronic allograft injury continues to be a major problem, and leads to eventual graft loss. Early detection of allograft injury is essential for guiding appropriate intervention to delay or prevent irreversible damage. Several advanced MRI techniques can offer some important information regarding functional changes such as perfusion, diffusion, structural complexity, as well as oxygenation and fibrosis. This review highlights the potential of multiparametric MRI for noninvasive and comprehensive assessment of renal allograft injury.

On the Field Strength Dependence of Bi- and Triexponential Intravoxel Incoherent Motion (IVIM) Parameters in the Liver

BACKGROUND

Studies on intravoxel incoherent motion (IVIM) imaging are carried out with different acquisition protocols.

PURPOSE

To investigate the dependence of IVIM parameters on the B₀ field strength when using a bi- or triexponential model.

STUDY TYPE

Prospective.

STUDY POPULATION

20 healthy volunteers (age: 19-28 years).

FIELD STRENGTH/SEQUENCE

Volunteers were examined at two field strengths (1.5 and 3T). Diffusion-weighted images of the abdomen were acquired at 24 b-values ranging from 0.2 to 500 s/mm² .

ASSESSMENT

ROIs were manually drawn in the liver. Data were fitted with a bi- and a triexponential IVIM model. The resulting parameters were compared between both field strengths.

STATISTICAL TESTS

One-way analysis of variance (ANOVA) and Kruskal-Wallis test were used to test the obtained IVIM parameters for a significant field strength dependency.

RESULTS

At b-values below 6 s/mm² , the triexponential model provided better agreement with the data than the biexponential model. The average tissue diffusivity was D = 1.22/1.00 μm² /msec at 1.5/3T. The average pseudodiffusion coefficients for the biexponential model were D^* = 308/260 μm² /msec at 1.5/3T; and for the triexponential model D 1 * = 81.3/65.9 μm² /msec, D 2 * = 2453/2333 μm² /msec at 1.5/3T. The average perfusion fractions for the biexponential model were f = 0.286/0.303 at 1.5/3T; and for the triexponential model f₁ = 0.161/0.174 and f₂ = 0.152/0.159 at 1.5/3T. A significant B₀ dependence was only found for the biexponential pseudodiffusion coefficient (ANOVA/KW P = 0.037/0.0453) and tissue diffusivity (ANOVA/KW: P < 0.001).

DATA CONCLUSION

Our experimental results suggest that triexponential pseudodiffusion coefficients and perfusion fractions obtained at different field strengths could be compared across different studies using different B₀ . However, it is recommended to take the field strength into account when comparing tissue diffusivities or using the biexponential IVIM model. Considering published values for oxygenation-dependent transversal relaxation times of blood, it is unlikely that the two blood compartments of the triexponential model represent venous and arterial blood.

LEVEL OF EVIDENCE

1 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2019;50:1883-1892.

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

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