Abstract
OBJECTIVE:
This study aims to assess the value of ultrasound real-time shear wave elastography (US-SWE) for evaluation of nonalcoholic fatty liver disease (NAFLD) in a rabbit model compared with multislice computed tomography (MSCT).
MATERIAL AND METHODS:
Twenty-six rabbits were fed with high-fat, high-cholesterol diet and six rabbits were fed with a standard diet. All rabbits were performed with MSCT and US-SWE at various time points to measure changes in liver parenchyma. The diagnostic efficiency of US-SWE was analyzed using receiver operating characteristics (ROC) curves compared with MSCT based on the liver pathology.
RESULTS:
The statistically significant differences in the areas under the ROC curves between using MSCT and US-SWE modalities were detected to discriminate between normal vs. NAFLD or higher severity pathology. Similarly, for normal or NAFLD vs. borderline or NASH livers, statistically significant differences between using US-SWE and MSCT modalities were also detected for nonalcoholic steatohepatitis (NASH) vs. lower severity pathology.
CONCLUSIONS:
MSCT, but not US-SWE, had a better ability to differentiate normal or NAFLD livers from higher severity NAFLD livers. However, the diagnostic efficiency of US-SWE was superior to that of MSCT for differentiating NASH from normal or lower severity NAFLD.
Keywords
Introduction
NAFLD is defined as intracellular deposition of fat droplets exceeding 5% of triglycerides in the patients without clinically significant alcohol intake, infection or any other liver disease induced by specific cause [1–4].The earliest and most significant histological characteristic of NAFLD is the presence of intracellular triglyceride accumulation within hepatocytes. Oxidative stress is caused by hepatic intracellular triglycerides deposition, leading to simple steatosis, liver injury, inflammation, NASH, and eventually progresses to hepatocellular carcinoma and cirrhosis [5, 6].NAFLD has been confirmed as a risk factor for cardiovascular morbidity and metabolic disease (including type-2 diabetes) [7–9]. Specifically, NAFLD severity groups were classified into normal, NAFLD, borderline NASH, and NASH according to NAFLD activity score as following [18]: (1) Normal: activity score = 0; (2) NAFLD: activity score = 1, 2; (3) Borderline NASH: activity score = 3, 4; and (4) NASH: activity score ≥5. Therefore, early diagnosis of severity and degree of NAFLD is of great importance for clinical practice.
Currently, there are no specific serum biomarkers for hepatic fat content quantification. Aminotransferases are not sensitive and specific enough for detection of liver fat content. Liver biopsy is considered the gold standard for evaluation of NAFLD. However, the procedure is not widely available because of its invasiveness. It is crucial for clinical practice to accurately quantify the liver with NAFLD or higher severity livers noninvasively.
Non-invasive imaging modalities for diagnosis of NAFLD include ultrasound (US), computed tomography (CT) and magnetic resonance imaging (MRI) [10, 11]. However, conventional imaging modalities do not detect specific signs of NAFLD owing to many factors. The non-imaging technique of controlled attenuation parameters (CAP) is an objective assessment for ultrasound attenuation, which can be performed to measure the hepatic steatosis by using the appropriately selected FibroScan probe [12, 13]. This parameter can be used as a validated marker of steatosis. MR spectroscopy (MRS) is a technique that the relative strength signals of water and fat were determined by frequency-domain analysis of MRS data. Liver fat quantification can be evaluated by MRS, which requires physicist or technologist to set up the sequence and acquire data, the MRS analysis software is also needed to analyze the data. Liver fat content was assessed with MRS, however, MRS is expensive and not broadly available, the accuracy of magnetic resonance-derived methods may decrease with significant fibrosis [14]. In comparison to MRI, the advantages of US or MSCT are its relatively low cost and ease of access [11]. A recent study showed that a novel technique, quantitative US-SWE, may be useful to measure the severity and degree of liver cirrhosis [15–17].
As for the quantitative assessment of liver density by MSCT, which is based on X-ray attenuation value or density that corresponds to the “brightness” of the pixels. MSCT is less prone to anatomical constraints than US allowing complete visualization of the entire liver within a short breath-hold. With increasing fat accumulation, the CT value of liver parenchymal decreases. Therefore, the severity of hepatocellular steatosis may be represented by means of a continuous range of CT attenuation values, the CT attenuation values of fatty liver range from that of normal liver parenchyma to that of fat tissue [11]. According to this theory, so we studied the diagnostic efficiency of MSCT for assessment of NAFLD in the rabbit model in comparison to US-SWE in a rabbit NAFLD model.
Materials and methods
This study was approved by the animal care and use committee of Guizhou Medical University.
The study design for rabbit NAFLD model
Thirty-six New Zealand white rabbits (male = 18, female = 18), 8-weeks-old, were randomly classified into control group (n = 9) or experimental groups (9 rabbits/group). Randomization procedure: the randomization procedure was based on the random number and block randomization by using SPSS software (version 23.0). The control group was fed a standard diet. The experimental group was fed a high-fat and high-cholesterol diet (standard diet with an additional 2% cholesterol, 6% yolk powder, 10% lard, and 2% maltose dextrin). All rabbits were performed US-SWE and MSCT before experiment begun, and then underwent US and MSCT at various time points: 4 weeks (group 4W), 8 weeks (group 8W), or 12 weeks (group 12W) for different groups, respectively. The body weight also was also recorded. This process was conducted by the Laboratory Animal Research Center of Guizhou Medical University (Guiyang, China). The study design was as following (Fig. 1):

A flow chart showing our study design of NAFLD rabbit models. 36 rabbits were randomly classified into a control group (n = 9) or an experimental group (n = 27). Control group was fed with a standard diet, and all experimental groups were fed with a high-fat and high-cholesterol diet; All the rabbits underwent US-SWE and MSCT prior to being fed with standard diet or high-fat and high-cholesterol diet. Rabbits in the experimental groups were sacrificed either at 4 weeks (n = 9), 8 weeks (n = 9), or 12 weeks (n = 9). Control rabbits were sacrificed at 8 weeks. All liver tissue specimens were stained with H&E and Masson trichrome staining. Histological features were classified as follows: normal, NAFLD, borderline NASH, and NASH in our study.
After the US-SWE examination, all rabbits underwent 128-slice MSCT scanning (SOMATOM definition AS+, Siemens, Germany). The scan parameters were as follows: Tube voltage 120 kVp, tube current 250 mA, collimation 128×0.6 mm, slice thickness1 mm, FOV 500 mm.
Ultrasound examination
All rabbits underwent ultrasound imaging by the Supersonic Imaging Aixplorer with a linear array probe (L15-4) and a frequency of 4–15 Hz. Skin preparation was performed by the ultrasonographer. Firstly, we obtained a conventional B-mode gray-scale image for assessment of the entire liver parenchyma. Secondly, SWE imaging was performed on the liver parenchyma for assessment of liver stiffness.
Image data analysis
All the rabbits underwent MSCT and US-SWE on the same day. All CT images datasets were reconstructed at a workstation (Syngo CT Workplace VE40B). Three ROIs (diameter 0.7 cm–1.0 cm) were measured on the liver parenchyma of each rabbit by an experienced abdominal radiologist who was blinded to the rabbit’s diet. Care was taken to avoid large bile ducts and hepatic vessels on the CT image. One ROI was placed in the left lateral lobe, one in the left medial lobe, and the third in the right lobe, using the left hepatic vein and the gallbladder as anatomic landmarks. CT attenuation values of each ROI were recorded, average CT value of three ROIs was calculated for analysis. As for the ultrasound images of SWE, sample frames were full-color (≥95%) for stabilization and color homogeneity, three ROIs (diameter 0.7 cm–1.0 cm) were measured in the liver parenchyma by Q-Box software, avoiding intrahepatic vascular and bile ducts, the mean elastic modulus values were recorded and calculated for statistical analysis.
Histological analysis
Following all US-SWE and MSCT examinations, all rabbits were scarified under deep anesthesia with 5 ml of potassium chloride injected into the auricular veins. All liver tissue specimens were stained with Masson staining and hematoxylin and eosin staining, respectively. Pathological sections were reviewed by an experienced pathologist who was blinded to rabbit’s diet. The reference standards for pathological specimen were as follows [18]: NAFLD activity score was given for each pathological section: the unweighted sum of steatosis (score 0–3), lobular inflammation (score 0–3), and hepatocellular ballooning (score 0–2).
After the NAFLD activity score was recorded, each study subject was classified into various NAFLD severity groups [18]: Normal: activity score = 0; NAFLD: activity score = 1, 2; Borderline: activity score = 3, 4; and NASH: activity score ≥5.
The METAVIR score system was employed for evaluation of liver fibrosis [19]: F0 = no fibrosis; F1 = portal fibrosis without septa; F2 = portal fibrosis with few septa; F3 = a number of septa without cirrhosis; F4 = cirrhosis.
Data analysis
Statistical analysis was performed using SPSS software (version 23.0, SPSS Inc., Chicago, IL) and MedCalc software (MedCalc 15.2.2, Mariakerke, Belgium). Continuous variables were expressed as mean±SD. The normality assumption of variables was conducted according to the Kolmogorov-Smirnov Z test. The body weight, CT value, and elastic modulus value were tested by two-tailed paired-samples T Test before and after experiment for different groups. The CT value and elastic modulus value of liver parenchyma among the grades of different histological features were analyzed by one-way ANOVA or two-tailed independent-sample.
Diagnostic efficiency of MSCT was tested and compared with US-SWE; Corresponding sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and the area under ROC curve (AUC) were also analyzed. 95% confidence interval (CI) was recorded. P < 0.05 was considered statistically significant difference.
Results
Characterization of NAFLD severity according to histologic features
Three rabbits died in control group before sacrifice date, and one rabbit died in group 12 W. At last, 32 rabbits (control group n = 6, experiment groups n = 26) were enrolled in our study. For the six animals in the control group, histological features were designated as normal liver.
For the nine rabbits in group 4W, all rabbits were fed 4 weeks on a high-fat, high-cholesterol diet, mean NAFLD activity score (NAS) was 2.56 (range, 2–4). Six rabbits were classified as NAFLD, three rabbits were borderline, but no rabbits were NASH or fibrotic.
Group 8W (n = 9) was maintained for 8 weeks on a high-fat, high-cholesterol diet. Mean NAS was 4.00 (range, 2–7). Two rabbits were classified as NAFLD, five as borderline, two as NASH. Two of five borderline and two NASH rabbits had the fibrosis of F1.
For group 12 W (n = 8), eight rabbits were fed 12 weeks of a high-fat, high-cholesterol diet. Mean NAS was 5.63 (range, 2–7). One rabbit was classified as NAFLD, one as borderline, six as NASH. One borderline and six NASH rabbits had pericellular fibrosis of F1.
In summary, according to histological feature scores, six rabbits were normal, nine had NAFLD, nine were borderline, eight had NASH, and eleven showed fibrosis of F1 (Table 1).
The NAFLD severity livers at different time points
The NAFLD severity livers at different time points
In addition, the body weight was increased in the control group, group 4 W, and group 12 W (all P < 0.05). As for the CT value of liver parenchymal, the CT value was decreased in all experiment groups (all P < 0.05). However, the elastic modulus value was increased in group 8 W and group 12 W (all P < 0.05). (Table 2 and Fig. 2).

Changes were observed in (A) body weight, (B) liver stiffness, and (C) liver density at various time points.
Dynamic changes of body weight, US-SWE, and MSCT for evaluation of liver parenchyma at different time points for various groups
“*” All the data were tested with two-tailed paired-samples T Test.
MSCT and US-SWE measurements for evaluation of histological features were shown in Fig. 3. There was statistically significant difference in steatosis grades for liver density and stiffness (P < 0.001 and P = 0.007, respectively). Statistically significant difference was observed in fibrosis grades (P = 0.002 and P = 0.026, respectively). Statistically significant difference was also existed in inflammation grades for MSCT and US-SWE (P = 0.001 and P < 0.001, respectively).

Box-plots of MSCT and SWE measurements for evaluation of histological features, respectively, including steatosis grade (A and B), Fibrosis grade (C and D), and inflammation grade (E and F). Boundaries of the boxes indicate the lower and upper percentiles; lines within the boxes indicate medians; and error bars indicate ranges for each grade of histological feature.
The statistically significant difference was observed to discriminate normal liver parenchyma from NAFLD livers between MSCT and US-SWE (Z = 4.599, P < 0.001). The difference of AUC between MSCT (AUC 0.949) and US-SWE (AUC 0.551) was 0.398 (95% CI 0.228, 0.567). For discrimination between normal liver or NAFLD versus borderline or NASH, there were also significant differences between US-SWE and MSCT (Z = 2.631, P = 0.009). The difference of AUC between MSCT (AUC 0.710) and US-SWE (AUC 0.941) was 0.231 (95% CI 0.059, 0.404); For discrimination of livers with NASH versus those with borderline or lower severity liver pathology, we also saw a statistically significant difference (Z = 2.003, P < 0.045). The difference of AUC between US-SWE (AUC 0.953) and MSCT (AUC 0.773) was 0.180 (95% CI 0.004, 0.355). The AUC, cutoff values for MSCT and US-SWE parameters, sensitivities, specificities, PPV, and NPV are displayed in Table 3, Figs. 4, and Figs. 5.

Box-plots show elastic modulus and density of liver parenchyma for normal, NAFLD, borderline, and NASH. (A) US-SWE for normal, NAFLD or higher severity pathology. The lowest line indicates the minimum value, the highest line indicates the maximum value, lines in the box indicate median, boundary of box indicate lower and upper quartile. (B) MSCT for normal, NAFLD or higher severity pathology.

ROC curve analysis for US-SWE and MSCT. (A) ROC curves show that there was significant statistical difference in the area under the ROC curve between MSCT and US-SWE to evaluate for normal vs. ≥NAFLD severity livers. (B) there was significant statistical difference between MSCT and US-SWE with respect to discrimination of normal or NAFLD vs. borderline or NASH livers. (C) Significant statistical difference was also existed between US-SWE and MSCT for NASH vs. lower severity livers.
Receive operating characteristic curve analysis for NAFLD severity livers by US-SWE and MSCT
“*” In comparison to US-SWE, there was significant statistical difference in the area under ROC curve, P < 0.05.
The results we got were that we detected a significant difference between US-SWE and MSCT to evaluate the NAFLD severity in a rabbit model. For the ability to discriminate 1) normal liver versus NAFLD, and 2) normal liver or NAFLD versus borderline liver or NASH, the diagnostic efficiency of MSCT was superior to that of US-SWE. However, for discrimination of livers with NASH versus borderline or lower pathological severity livers, the diagnostic performance of US-SWE was superior to that of MSCT. This suggests that MSCT and US-SWE provide different diagnostic efficiencies in diagnosis of the severity of NAFLD.
We chose the rabbit as NAFLD model in our study due to its size. The body weight of rabbit was observed in our study because we needed to know the volume of pentobarbital sodium for anesthesia before we performed US-SWE and MSCT for assessing the liver stiffness and density dynamically. In our study, the body weight was decreased at 8 weeks compared with control group. Histological results shown that the rabbit NAFLD livers had pericellular fibrosis by Masson staining after 8 weeks. It is hepatic fibrosis and damage that previous studies suggested may be associated with body weight; In addition, liver fibrosis presented with anorexia and emaciation were observed in animal study [20, 21]. In our study, the liver had fibrosis in the borderline NASH or NASH rabbit model after 8 weeks (Fig. 6).

Ultrasound SWE image for assessment of liver parenchyma. The normal liver of (A) MSCT, (B) US-SWE, (C) hematoxylin and eosin staining, and (D) Masson staining; the NAFLD liver of (E) MSCT, (F) US-SWE, (G) hematoxylin and eosin staining, and (H) Masson staining.
Histological features for NAFLD livers were classified as normal, NAFLD, borderline NASH, and NASH according to the unweighted sum of steatosis, lobular inflammation, and hepatocellular ballooning. In our study, the CT value of liver parenchyma was decreased when increasing the steatosis grades, but the stiffness of liver parenchyma was increased. Previous human or animal studies have shown that the liver stiffness increased when increasing the fibrosis stages or inflammation grades [22–24]. In our study, the elastic modulus value of F1 was higher than that of F0, and the elastic modulus value was increased when increasing the inflammation grades. However, the histological features of steatosis, fibrosis, and inflammation coexisted in the NAFLD liver specimen in our study (11/26). Therefore, it is very difficult to evaluate the difference of liver density or stiffness induced only by steatosis, fibrosis or inflammation, respectively.
The main advantages of conventional ultrasound are ease of access and low cost. However, many factors affect the quality and accuracy of ultrasound, such as type of equipment, the experience of operator, and anatomical constraints [25]. In addition, sonographic images are confounded in advanced NAFLD because of coexisting inflammation and fibrosis, therefore, conventional ultrasound is not recommended for diagnosis of mild steatosis or for assessing the severity of NAFLD [26–28]. For mild steatosis, the diagnostic sensitivity and specificity of sonographic features were 62.2–82.1% and 76.2% –90.1%, respectively [10]. To overcome the subjectivity of conventional ultrasound, several quantitative ultrasound techniques for assessment of liver lesions have been utilized [29–31]. Real-time SWE is a quantitative technique with excellent severity grading performance and reproducibility [32]. Several studies have shown that SWE may be used in quantitative ultrasound elastography to assess liver stiffness associated with cirrhotic features [33–35].
In general, fatty liver disease develops from a simple fatty liver to steatohepatitis, and eventually to fatty cirrhosis [36, 37]. However, there are exceptions, for example, fatty liver disease may not exhibit clear steatohepatitis, but directly evolve to cirrhosis [36]. In group fed with high-fat and high-cholesterol diet for 12 weeks, we suspect that higher elastic modulus value may be associated with the histological feature of fibrosis.
MSCT is based on the attenuation values which are standardized by means of Housfield Units (HU). On the unenhanced CT image, the attenuation value of liver parenchyma is approximately 60±10 HU. Subcutaneous adipose tissue is approximately –90 HU. Therefore, the severity of hepatocellular steatosis may be represented by means of a continuous range of CT attenuation values from normal liver parenchyma to fat tissue [11]. However, CT attenuation value can be confounded by many factors for evaluation of hepatic steatosis [38]. In addition, accumulation of water, iron, collagen, and glycogen often occurs in the diffuse liver disease [39]. All these factors may affect the CT attenuation value of liver parenchyma and may confound the evaluation of hepatic steatosis. In our study, the CT attenuation value of rabbit liver parenchyma decreased gradually in all experimental groups.
The mean elastic modulus value increased as the severity of NAFLD increased (P < 0.05). There was a statistically significant difference in discrimination capability for NASH from borderline or lower severity pathology, with an AUC value of 0.953 (0.814, 0.997). In comparison to MSCT, the diagnostic efficiency of SWE was superior to that of MSCT (P < 0.05). Using MSCT, the CT number for liver was decreased compared with normal liver as the severity of NAFLD increased (P < 0.05). There was a statistically significant difference between modalities for discrimination between normal liver and liver with various severities of pathology. However, the diagnostic efficiency of MSCT was superior to that of US-SWE (P < 0.05) for differentiating 1) normal liver from NAFLD, and 2) normal or NAFLD from borderline or NASH; AUC value 0.949 (0.808, 0.996) and 0.940 (0.797, 0.993), respectively.
There are some limitations in our study. First, the sample size was small, including the number of subjects in the control group and experimental group. Second, we used a semiquantitative grading system for evaluation of histological features, including steatosis, hepatocellular ballooning, lobular inflammation, and fibrosis. Our study showed that the histological feature of rabbit liver specimens with NAFLD was micro-vesicular steatosis predominately (24/26). Previous studies demonstrated a correlation between histological features of micro-vesicular steatosis and advanced histological results of NAFLD [40]. Third, the information about fat content of liver tissue was not measured directly in our study, which may better reveal the steatosis of NAFLD livers. Fourth, ROIs were not a consistent one-to-one match for measuring the density and stiffness of the liver parenchyma by MSCT and US-SWE, respectively. Last, the repeatability and reproducibility of US-SWE and MSCT were not performed in this study.
In conclusion, the results we got are promising. As for these two quantitative techniques of US-SWE and MSCT, different modalities with different diagnostic efficiencies for assessing the different stages of NAFLD livers. MSCT, but not US-SWE, had a better ability to differentiate normal or early NAFLD livers from higher severity NAFLD livers, which suggests that MSCT could be recommended to evaluate or monitor the early stage of NAFLD livers. However, the diagnostic efficiency of US-SWE was superior to that of MSCT for differentiating NASH from normal or lower severity NAFLD, which suggests that US-SWE could be recommended to assess or follow up the dynamic changes for the advanced stage of NAFLD livers. It means that the diagnostic efficiencies of these two techniques are mutually supplementary for assessment of NAFLD.
Conflict of interest
The authors do not have any disclosures to report.
Footnotes
Acknowledgments
We are thankful to Dr. Lujun Dai and Dr. Yining Xiang who supported the techniques of H&E staining and Masson trichrome staining.
This study was supported by Joint Fund of Guizhou Province Department of Science and Technology (No. LH[2017]7208); National Natural Science Foundation of China (No. 81760312); Guizhou Province Clinical Discipline Subject Fund (No. QWBF2011[52]-03); Doctoral Research Initiation Fund of the Affiliated Hospital of Guizhou Medical University.
