Abstract
This study investigated a combined use of intravoxel incoherent motion (IVIM) parameters, Dslow (
Keywords
Introduction
Chronic liver disease is a major public health problem worldwide. The epidemic trend of chronic liver disease is expected to increase due to an aging population, the growing epidemic of obesity, and nonalcoholic steatohepatitis (NASH). Viral hepatitis is the most common blood-borne infection worldwide.1,2 Chronic viral hepatitis can lead to hepatic fibrosis, cirrhosis, and hepatocellular carcinoma. 3 Liver fibrosis, a common feature of almost all chronic liver diseases, involves the accumulation of collagen, proteoglycans, and other macromolecules in the extracellular matrix. 4 Clinically, liver fibrosis usually has an insidious onset and progresses slowly over decades. Originally considered to be irreversible, hepatic fibrosis is now regarded as a dynamic process with the potential for regression. 4 Treatment with combined therapies on underlying etiology and fibrosis simultaneously might expedite the regression of liver fibrosis and promote liver regeneration.5–7 Earlier-stage liver fibrosis is more amenable to therapeutic intervention. Even when the underlying etiology of liver fibrosis could not be eradicated, therapies on liver fibrosis might help delay the progression of the disease to cirrhosis.
To date, noninvasive diagnostic tests available from clinical practice are not sensitive or specific enough to detect occult liver injury at an early stage. 8 Liver biopsy is currently the standard of reference for the diagnosis and staging of liver fibrosis. However, liver biopsy is an invasive procedure with several contraindications and with a risk of complications, such as pain, hemorrhage, bile peritonitis, penetration of abdominal viscera, pneumothorax, and even death.9,10 The mortality rate associated with needle biopsy was estimated to be between 0.009% and 0.12%. 10 A noninvasive and quantitative technique for detecting liver fibrosis is highly desirable.
In diffusion-weighted (DW) MRI, the intensity of the acquired MR signal depends on the self-diffusion of the excited spins, that is, on the microscopic stochastic Brownian molecular motion, and the extent and orientation of molecular motion is influenced by the microscopic structure and organization of biological tissues.11–14 Perfusion can contribute to the diffusion measurements significantly because of the incoherent motion of blood in the pseudorandom capillary network at the macroscopic level.15–18 Intravoxel incoherent motion (IVIM) reflects the random microscopic motion that occurs in voxels on MR images of water molecules (either intracellular or extracellular) and the microcirculation of blood. Le Bihan et al.15,16 proposed the principle of IVIM, which enables the quantitative parameters that separately reflect tissue diffusivity and tissue microcapillary perfusion to be estimated. IVIM signal attenuation is modeled according to the equation
where
Molecular water diffusion in fibrotic liver would be restricted by the presence of collagen fibers in the distorted lobular structure. Given the relatively high blood volume fraction of <25–30 mL of blood per 100 g in liver,
19
perfusion can contribute to the diffusion measurements significantly because of the incoherent motion of blood in the pseudorandom capillary network at the macroscopic level. It is well accepted that liver fibrosis is associated with reduced liver perfusion.20–23 Recently, there has been greater interest in using the IVIM technique to study diffused liver diseases, such as liver fibrosis. However, so far the literature has shown that IVIM is unable to detect early liver fibrosis reliably.
24
We noticed that the potential optimal combination of three IVIM parameters, Dslow, PF, and Dfast, for the detection of liver fibrosis has not been explored in sufficient detail. In this study, we set out to explore whether a combination of Dslow, PF, and Dfast can be used to separate fibrotic livers from healthy livers. We reanalyze our previously reported cohort, the
Materials and Methods
The characteristics of the Shenzhen 2012/2013 dataset have been previously reported.
25
The MRI data were acquired during the period from August 1, 2012, to August 15, 2013. The study was approved by the ethical committee of Shenzhen No. 3 Hospital, and informed consent was obtained, which included secondary analysis for the data acquired. The IVIM images of one volunteer and one patient were adjudged to contain substantial motion artifacts and were therefore excluded for analysis. Sixteen healthy volunteers (10 males, 6 females; mean age 36.4 years old; range 21–79 years old) and 33 consecutively viral hepatitis-b patients were included in the current study. The patient cohort had 15 stage F1 subjects (mean age 31.8 years, range 22–53 years) and 18 stage F2–4 subjects (mean age 42 years, range 22–53 years). The histology diagnosis for liver fibrosis was based on the
MRI was performed with a 1.5 T magnet (Achieva, Philips Healthcare, Best, Netherlands). The IVIM DW imaging sequence was based on a single-shot DW spin-echo-type echo-planar imaging sequence, with 10
All curve-fitting algorithms were implemented in an accustom program developed on MatLab (Mathworks, Natick, MA). Regions of interest (ROIs) were positioned to cover a large portion of liver parenchyma while avoiding large vessels ( Fig. 1 ). For ROI analysis, the IVIM parameters were calculated based on the mean signal intensity of the whole ROI, which has been shown to offer better estimation than pixel-wise fitting when the signal-to-noise ratio (SNR) of the DW images is low.30,31

(
The following measurement modifications were made compared with our previous report.
25
As the left lobe of the liver is more likely to suffer from artifacts associated with cardiac motion and B0 inhomogeneity susceptibility due to its proximity to the stomach and the air inside, in the current study, only the right lobe of the liver was measured (
Fig. 1
).
Figure 1
demonstrates that the ROI was carefully drawn to cover only the right liver parenchyma while avoiding vasculature and artifacts. All six slices per subject were evaluated, while the slices with notable motion artifacts and those demonstrating notable outliers with signal
A three-dimensional tool was programed using IBM SPSS 23 for Windows (SPSS Inc., Chicago, IL), and the measures of Dslow, PF, and Dfast were placed along the
where x(i) is the original data and z(i) is the normalized data, and xmax and xmin are the maximum and the minimum value of x(i), respectively. Note the range of z(i) after normalization is from 0 to 1 for each dimension.
The classification and regression tree (CART) model was used to find the cutoff values for PF, Dslow, and Dfast to differentiate F0 livers versus F1–4 livers and F0 livers versus F2–4 livers. 33
Results
Livers of healthy volunteers had a PF of 16.6% ± 3.6% (mean ± standard deviation [SD]), Dslow of 1.14 ± 0.22 (×10−3 mm2/s), and Dfast of 12.3 ± 3.1 (×10−3 mm2/s), respectively. The coefficient of variance (CoV, SD/mean) for PF, Dslow, and Dfast in healthy volunteers was 0.19, 0.22, and 0.25, respectively.
When the study participants were grouped into three groups, (1) healthy volunteers (F0), (2) insignificant liver fibrosis (F1), and (3) significant liver fibrosis (F2–4), it was seen that PF offered the best differentiation of the three groups, followed by Dslow (
Table 1
,
Fig. 2
). By adjusting the viewing angel, the three-dimensional visual tool demonstrated healthy volunteers (F0,
Mean, SD, and CoV of PF, Dslow, and Dfast of Healthy Volunteers, F1 Liver Fibrosis Patients, and F2–4 Liver Fibrosis Patients.

Scattered plots and mean of PF, Dslow, and Dfast of healthy volunteers, F1 liver fibrosis patients, and F2–4 liver fibrosis patients.

Three-dimensional display of healthy volunteer group (blue balls), F1 patient group (pink balls), and F2–3 patient group (red balls). Each ball represents one participant. The differentiation of the volunteer group and patient group can be better visualized by rotation in a three-dimensional space (dotted yellow line or purple plane).

Three-dimensional display of healthy volunteer group (blue balls) and F2–3 patient group (red balls). Each ball represents one participant. The differentiation of the volunteer group and the patient group can be better visualized by rotation in a three-dimensional space (dotted yellow line or purple plane).
Quantitative analysis with SVM showed healthy volunteers, and all patients with liver fibrosis (F1–4) were differentiated with a plane defined by (166.58*PF) + (8.90*Dslow) – (0.98*Dfast) –19.71 = 0 ( Fig. 3 ); healthy volunteers and patients with significant liver fibrosis (F2–4) were more clearly differentiated with a plane defined by (29.56*PF) + (4.33*Dslow) – (0.12*Dfast) – 6.67 = 0 ( Fig. 4 ). The mean distance of the data points for F0 versus F1–4 to the central plane was 0.0021F0 + 0.0026F1–2 = 0.0047, and for F0 and F2–4 to the central plane, 0.0149F0 + 0.0138F2–4 = 0.0287.
The three-dimensional visual tool demonstrated better differentiation of healthy livers and fibrotic livers than the two-dimensional plot using PF and Dslow values ( Figs. 3 – 5 ), indicating that Dfast contributed to differentiating healthy volunteers and patients with liver fibrosis.

Two-dimensional demonstration of healthy volunteer group (blue balls), F1 patient group (pink balls), and F2–3 patient group (red balls) using the PF axis and Dslow axis. Each ball represents one participant. A combination of the PF axis and Dslow axis is insufficient to differentiate F0 subjects vs. F1–4 patients, in contrast to demonstrations in
Figs. 3
The CART analysis result is shown in Table 2 ; a combination of PF (cutoff value: PF < 12.55%), Dslow (Dslow < 1.152 × 10−3 mm2/s), and Dfast (Dfast < 13.36 × 10−3 mm2/s) can differentiate healthy subjects (F0) and fibrotic livers (F1–4) with an area under the curve of logistic regression (AUC) of 0.986. The AUC for differentiation of healthy subjects (F0) from significantly fibrotic ones (F2–4) was 1.
CART Analysis of AUC of F1–4 or F2–4 (compared with F0) on PF, Dslow, and Dfast.
Unit of Dslow and Dfast: ×10−3 mm2/s.
Discussion
Currently, the most clinically used imaging technique for evaluation of liver fibrosis is ultrasound elastography, while the investigational technique of MR elastography has undergone many clinical trials. 34 The ultrasound elastography method is inexpensive, fast to acquire, and does not require postprocessing. In particular, one-dimensional transient elastography has been adopted clinically. The diagnostic accuracy of ultrasound elastography has been assessed in numerous studies and pooled in meta-analyses. The reported diagnostic accuracy for various ultrasound elastography techniques has an AUC in the range of 0.84–0.87 for fibrosis stage ≥F2, 0.89–0.91 for fibrosis stage ≥F3, and 0.93–0.96 for fibrosis stage F4.35–37 MR elastography provides higher overall diagnostic accuracy than ultrasound-based elastography. 34 Meta-analyses report an AUC for MR elastography in the range of 0.84–0.95 for diagnosing fibrosis stage ≥F1, 0.88–0.98 for fibrosis stage ≥F2, 0.93–0.98 for fibrosis stage ≥F3, and 0.92–0.99 for fibrosis stage F4.38–40 The limitations of ultrasound elastography include that it is an operator-dependent technique, while MR elastography requires additional hardware (an elastography driver) to be added to the MR scanner. As the IVIM imaging sequence is widely available in clinical MR scanners, it represents a promising alternative to existing techniques for liver fibrosis evaluation.
The PF and Dslow measurements obtained in this study broadly agreed with previous reports.24,41–46 For the 27 studies that reported measurement for healthy livers, the median value for Dslow was 1.11 × 10−3 mm2/s and 1.02 × 10−3 mm2/s at 1.5 and 3 T, respectively; the median value for PF was 22.00% and 22.65% at 1.5 and 3 T, respectively.
24
Using an acquisition of 25
Although much research has been published on the evaluation of liver fibrosis using IVIM, how to optimally combine Dslow, PF, and Dfast to get diagnostic information is not yet explored. The most important result of the current study is that despite the Shezhen 2012/2013 ivim dataset not being acquired with an optimized protocol, we were still able to demonstrate that healthy volunteers and patients with liver fibrosis could be differentiated. Another important finding of this study is that among the IVIM parameters of PF, Dslow, and Dfast, PF offers the best diagnostic value and Dfast can provide additional differentiation value, although it is a less stable measurement. Overall, the cluster of F1 subjects was located between F0 and F2–4; however, it was not possible to differentiate patients of different stages of F2–4 in this study. We expect this is at least partially due to the fact that the histological diagnosis is also not clear-cut; a high-end F1 liver will be similar to a low-end F2 liver. 47 The findings of this study are important, as until now it has been considered that there is no reliable noninvasive method, whether imaging or serum biomarkers, that can reliably detect early-stage liver fibrosis. In the meantime, we are looking into further validating our approach with other datasets or new prospective studies. Additionally, a two-dimensional flat plane was used in this study to separate healthy livers versus fibrotic livers; theoretically curved planes can also be used if curved planes offer better separation of healthy livers versus fibrotic livers, or better staging of fibrotic livers.
IVIM parameters strongly depend on the choice of the
There are a few ways to improve the measurement accuracy in this study. This study used only 10
In conclusion, a combination of PF, Dslow, and Dfast shows the potential of IVIM to detect early-stage liver fibrosis. Among the three parameters, PF offers the best diagnostic value, followed by Dslow; however, all three parameters contribute to liver fibrosis evaluation. Further research shall improve image data postprocessing, denoise or remove poorly fitted regions in the liver, and validate our approach with additional datasets.
Footnotes
Acknowledgements
The authors thank Dr. Jing Yuan, Medical Physics and Research Department, Hong Kong Sanatorium and Hospital, for setting up the data acquisition protocol at the Shenzhen No. 3 People’s Hospital.
Supplementary material is available online with this article.
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by a direct grant for research from the Chinese University of Hong Kong (No. 4054167).
Notes
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
