Phenotypic characteristics of mouse models of diet-induced NASH
Mice fed a normal chow (NC) diet were used as the control group. Figure 1A presents the nutritional composition of each diet and the percentage of total energy contributed by each dietary component. The feeding schedule is depicted in Fig. 1B. The changes in the body weights (BWs) of the mice in the three experimental diet groups were recorded. A rapid increase was noted in the BW of all three groups; however, measurements performed when the mice were aged 24 wks revealed that BW gain was faster in the HFD and WD groups than in the HFC group (Fig. 1C). After 16 and 32 wks, the livers were larger and lighter (in color) in the experimental groups than in the control group (Fig. 1D).
The ratio of liver weight (LW) to BW is a robust and vital indicator of liver homeostasis. Normally, the LW/BW ratio ranges from 4.5–5%. Changes in the LW/BW ratio indicate disease progression. After 16 wks of feeding, we noted no significant change in the LW/BW ratio in the four (HFD, WD, HFC, and NC) groups. By contrast, after 32 wks of feeding, the ratio increased significantly in the three experimental groups compared with the ratio in the control group (Fig. 1E). Thus, prolonged intake of the experimental diets led to pathological changes in the mouse liver that accelerated diet-induced NAFLD/NASH. Notably, BW gain was slower and the LW/BW ratio was higher in the HFC group than in the other groups; this indicates a more rapid progression of NASH in the HFC group than in the HFD and WD groups. In summary, all three experimental diets induced NAFLD/NASH with varying degrees of severity.
Importance of clinically relevant diagnostic methods for the severity of diet-induced NASH in mice
Clinically relevant quantitative diagnosis was performed to evaluate the severity of diet-induced NASH. The results of histopathological staining performed using the liver tissue samples of each group are presented in Fig. 2A. The results of H&E staining revealed no abnormalities in the control group but gradual pathological changes (e.g., hepatic steatosis and hepatocyte ballooning; after 16 and 32 wks) in the experimental groups (Fig. 2A). Moreover, Masson’s trichrome staining revealed prominent collagen fibers (liver fibrosis) in the WD and HFC groups (after 16 and 32 wks; Fig. 2A). The METAVIR scoring system was used to evaluate the degree of liver fibrosis.[45] The METAVIR grading criteria were as follows: F0, no fibrosis; F1, portal fibrosis without septa; F2, portal fibrosis with a few septa; F3, portal fibrosis with numerous septa without cirrhosis; F4, cirrhosis and single-blinded diagnosis by a certified pathologist. No significant difference was found among the four groups with respect to fibrosis scores after 16 wks of feeding (Fig. 2B). However, after 32 wks of feeding, the fibrosis scores of the WD and HFC groups were higher than those of the NC and HFD groups; this finding suggests that liver fibrosis was more severe in the WD and HFC groups than in the NC and HFD groups (Fig. 2B). After 16 wks of feeding, no significant changes were noted in the level of collagen deposition in the four groups. However, after 32 wks of feeding, the level of collagen deposition increased significantly in the WD and HFC groups compared with the levels in the NC and HFD groups; this finding indicates severe fibrosis in the WD and HFC groups (Fig. 2C).
In addition to histopathological analysis, this study quantified a clinically used serological biomarker of NAFLD/NASH—hyaluronan.[45] The serum level of hyaluronan was measured through sandwich ELISA (Fig. 2D). No significant changes were noted in the level of hyaluronan in the NC group. Nevertheless, the level of hyaluronan increased significantly in the HFD and WD groups; the highest degree of increase was noted in the HFC group. However, the gradual changes observed in the level of hyaluronan in the HFD group were nonsignificant. These findings indicate prolonged intakes of the WD and HFC increased the level of hyaluronan in the blood of the experimental mice (Fig. 2D). Furthermore, the WD and HFC accelerated NASH progression, leading to advanced, conditions such as liver fibrosis and cirrhosis. In summary, the severity of NAFLD/NASH induced by the three experimental diets could be quantified and scored using clinically relevant diagnostic methods.
Lipoprotein metabolites are key biomarkers of diet-induced NASH in mice
The NAFLD/NASH phenotype induced by experimental diets is quantifiable and scalable; thus, the corresponding serological metabolomic characteristics can be aligned. Mouse serum samples were subjected to NMR spectroscopy–based metabolomic profiling followed by bioinformatics analyses. The analytical logic and diagram are presented in Fig. 3. The data corresponding to serological metabolome and collagen deposition (Masson’s trichrome staining) were discovered to be correlated. A total of 41 biomarkers and 121 lipoprotein metabolites were assessed (Fig. 3). PCA and heatmap analysis were performed to identify the overall pattern and trends of changes in metabolites. Multiple comparisons were performed to determine the significant differences between the groups (significance was observed after 32 wks of feeding; Kruskal–Walls test) and between the feeding durations (16 and 32 wks; time effect; Wilcoxon rank-sum test). Additional multiple comparisons were performed for the collagen-related quantitative data (Fig. 3) to determine the differences between the feeding durations; various metabolites were identified in this analysis (Fig. 5A,B). In the case of significant between-group differences in collagen deposition level, correlation and comparative analyses of metabolites and collagen deposition levels were performed. A trimmed list of significant metabolites was obtained (Fig. 6A). Finally, the most significant metabolites (n = 17) were selected from the WD and HFC groups (Fig. 6B).
PCA was performed for the unbiased clustering of the obtained metabolomic data. All metabolites (162 biomarkers), the lipoprotein metabolites alone (121 biomarkers), and the small metabolites alone (41 biomarkers) were assessed. For all metabolites, the first principle component (PC1) accounted for 18.9% of the overall variability, whereas the second principle component (PC2) accounted for 35.3%. The PCA plot based on the top two principal components indicated that all samples could be categorized into NC, 16-wk-experimental-diet, and 32-wk-experimental-diet groups (Supplementary Fig. 1A). For the lipoprotein metabolites, PC1 accounted for 44.1% of the overall variability, whereas PC2 accounted for 22.5%. The PCA plot indicated that all samples could be categorized into NC, 16-wk-experimental-diet, and 32-wk-experimental-diet groups (Fig. 4A). For the small metabolites, PC1 accounted for 18% of the overall variability, whereas PC2 accounted for 12%. The PCA plot indicated that the samples could not be categorized (Fig. 4B). Heatmaps were constructed to visualize the distribution of the significant metabolites across groups. The heatmap corresponding to all metabolites revealed significant differences among the NC, 16-wk-experimental-diet, and 32-wk-experimental-diet groups (Supplementary Fig. 1B). The heatmap corresponding to the lipoprotein metabolites revealed significant differences among the NC, 16-wk-experimental-diet, and 32-wk-experimental-diet groups (Fig. 4C). The heatmaps revealed no significant difference between the groups in terms of serological biomarkers (Fig. 4D). These findings are consistent with the trends of the corresponding PCA results. In summary, the PCA and heatmap analysis indicated that lipoproteins undergo gradual changes and play major roles in the progression of NAFLD/NASH.
VLDL–VLDLR axis is associated with lipid deposition and steatosis symptoms in mouse models of diet-induced NASH
VLDL and LDL are the predominant metabolites in mouse models of diet-induced NASH
The association between phenotypic and metabolomic characteristics indicated that lipoproteins can play a pivotal role in the development of NAFLD/NASH. Differential expression analysis was performed to understand the gradual changes in lipoproteins during the progression of diet-induced NASH. The nonnominal method (Wilcoxon rank-sum test) was used for statistical analysis.[39] Significant biomarkers were selected on the basis of a p value of < 0.05 and an log2FC value of > 1. Pie charts depicting the results of differential gene expression analysis are presented in Fig. 5A. In the HFD group, significant changes were observed in a total of 36 biomarkers, accounting for 14% of all primary biomarkers. VLDL, LDL, and high-density lipoprotein (HDL) accounted for 8%, 58%, and 20%, respectively, of all biomarkers (Fig. 5A). In the WD group, the primary biomarkers, VLDL, LDL, and HDL accounted for 13%, 15%, 61%, and 11%, respectively, of all biomarkers (Fig. 5A). In the HFC group, the primary biomarkers, VLDL, LDL, and HDL accounted for 14%, 20%, 42%, and 24%, respectively, of all biomarkers (Fig. 5A). The changes in lipoprotein subfractions were ranked and are shown in Fig. 5B. The names, abbreviations, and sizes of all lipoprotein biomarkers are presented in Supplementary Table 3. In the HFD group, the top five metabolites (small LDLs) were L4TG, L4PN, L4AB, L4CH, and L4PL; the expression of these metabolites was considerably upregulated with time. However, the expression of small HDLs, such as H4PL, H4FC, and H4CH, was downregulated. In the WD group, the expression of the following three predominant metabolites was upregulated: the large VLDLs V2TG and V2CH and the large LDL L4TG. The size and number of various lipoproteins exhibited an average distribution; significant increases were particularly observed for LDLs, including for L1TG, L4TG, L2T, L3AB, L3PN, L3CH, L1FC, and L4FC. Regarding HDL, low levels of increases were observed in large HDLs, such as H1CH, H1FC, and H1A1. In the HFC group, considerable changes were noted with time in five VLDLs and two LDLs: the large VLDLs V2TG, V2CH, V1FC, and V4CH and the medium to large LDLs L1TG and L4FC. A significant reduction was noted in one HDL: H4CH. In summary, we observed gradual changes in lipoproteins with the changing severity of NAFLD/NASH. In the HFD, WD, and HFC groups, major changes were observed in small LDLs, large VLDLs and medium LDLs, and large VLDLs and medium to large LDLs, respectively.
To identify the correlation between metabolites and NAFLD/NASH severity, we performed correlation analysis of significant metabolites with collagen deposition (Fig. 6A). A total of 25 significant biomarkers of NAFLD were identified in all groups fed for 32 wks. In total, 34 significant biomarkers of NASH were identified in the WD group; moreover, 60 significant biomarkers of NASH/cirrhosis were identified in the HFC group (Fig. 6A). A comparison of the significant metabolites identified from the experimental groups with the metabolites that were significantly correlated with collagen deposition revealed a high degree of overlapping (Fig. 6B). Pathological analysis (Fig. 2) revealed that although both WD and HFC induced NASH, HFC induced a more severe condition—liver fibrosis (e.g., cirrhosis-like phenotype). Therefore, we further compared the WD group (mild fibrosis) with the HFC group (severe fibrosis) to explore fibrosis-specific metabolites. A total of 17 metabolites were found to be associated with severe liver fibrosis (Fig. 6C). VLDLs and LDLs accounted for 35% and 47%, respectively, of the aforementioned metabolites; both VLDLs and LDLs were found to be predominant in the liver of mice with NASH/cirrhosis (Fig. 6C). The 17 markers included the large VLDLs V1CH, V1PL, V2CH, V2TG, and V4CH; the small LDLs L5PN, L5CH, L5PL, and L5AB; and the HDL HDTG (Fig. 6C). Taken together, the results indicate VLDL and LDL are involved in the development of NAFLD/NASH and may induce cirrhosis. These findings elucidate both the roles of VLDL and LDL as biomarkers of NASH/cirrhosis and the pathophysiological changes that occur during the progression of NASH/cirrhosis.
Upregulation of VLDLR expression in mouse models of diet-induced NASH
The liver is the most prominent contributor of lipoproteins because this organ is responsible for both the production and recycling of lipoproteins.[46] Lipoprotein receptors are crucial for systemic lipid metabolism. The expression of lipoprotein receptors, such as VLDLR, LDL receptor (LDLR), and HDL receptor (SR-B1), in normal organs has been studied in humans and mice.[46, 47] VLDL is believed to be produced only by the liver; VLDLR is expressed in the periphery of but not within the liver.[48] In our study, the expression of both LDL and VLDL was upregulated in mouse models of NASH/cirrhosis, which prompted us to investigate the receptors of these lipoproteins in diseased liver tissues. The expression and distribution of VLDLR, LDLR, and SR-B1 were evaluated through immunohistochemical analysis (Fig. 7A) and was quantified using ImageJ (Fig. 7B). The expression of VLDLR was not similar between the diet- or feeding time–based groups, with the exception of the HCF group, which was fed for 32 wks (Fig. 7B). After 16 wks of feeding, the expression of LDLR was markedly downregulated in the experimental groups compared with that in the control group; nonetheless, the expression was gradually restored in the WD and HFC groups after 32 wks of feeding (Fig. 7B). Notably, the expression of SR-B1 remained high and did not change with diet (Fig. 7B). In summary, the expression of VLDLR is considerably upregulated in severe liver fibrosis. The findings of increases in the levels of serological VLDL and LDL and the upregulation of VLDLR expression in the severe NASH/cirrhosis of this study indicate a feedforward mechanism for lipid deposition.
VLDLs serve as the biomarkers of liver fibrosis/cirrhosis in humans
The expression of large VLDLs and VLDLR are upregulated in mouse models of diet-induced NASH with fibrosis/cirrhosis (Figs. 6C and 7). We analyzed the clinical specimens of a retrospective cohort of patients with liver fibrosis/cirrhosis to identify the correlation between serological metabolites and clinical features. FIB-4 is the most widely used noninvasive tool for evaluating the degree of liver fibrosis. This tool was developed to noninvasively predict severe fibrosis in patients with decompensated liver disease.[49, 50] The demographic characteristics of our cohort are summarized in Supplementary Table 4. On the basis of their METAVIR scores, the patients were stratified into mild and severe disease groups. A comparison of the metabolome and differential expression of relevant genes were performed. The results revealed significant increases in the levels of the following metabolites (very large VLDLs) in patients with severe fibrosis: XXL_VLDL_CE, XXL_VLDL_C, L_VLDL_CE, and L_VLDL_C (Fig. 8A). Human and experimental (mouse) NASH disease were similar in terms of the upregulation of the expression of very large VLDLs. In summary, the severity of diet-induced NASH in mouse models can be evaluated to align with clinical diagnostic methods. Metabolomic profiling revealed a likely mechanism of VLDL recycling through VLDLR, which may be involved the pathogenesis of liver fibrosis/cirrhosis.