3.1. Effect of HNK on the hepatic biomarkers in the tested mice. Hepatic biomarkers of the mice were detected by commercial test kits. The results showed that HNK had no effect on weight loss and liver index of mice fed by MCD diet ((Fig. 1(a) and 1(b)). Compared with MCS group, serum ALT, AST and liver TG, MDA levels in MCD group were significantly increased. After HNK treatment, serum ALT activity, liver TG and MDA content decreased remarkably, AST level was reduced similarly, but there was no significant difference compared with MCD group. In addition, the levels of liver TG, MDA and serum AST in OCA group were obviously lower than those in MCD group (Fig. 1(c)-1(f)).
3.2. HNK attenuated hepatic steatosis, inflammation and fibrosis of NASH mice. Hepatic lipid accumulation, steatosis, inflammation and fibrosis of the mice were analyzed by Oil red O staining, HE staining and Masson staining (Fig. 2). Compared with MCS group, the liver cells of MCD mice were disordered and showed ballooning changes, with a large number of lipid droplets and fat vacuoles, accompanied by inflammatory cell infiltration and fibrosis. After treatment with HNK and OCA, the number of lipid droplets and fat vacuoles in liver cells were reduced significantly, the degree of inflammatory cell infiltration and fibrosis decreased, and the liver tissue morphology was closer to MCS group.
3.3. Effects of HNK on serum lipids of NASH mice. The principal component analysis (PCA) showed the separation trend of serum lipid composition in each group. As shown in Fig. 3(a), there was an obvious separation trend between the MCS group and the MCD group, indicating that MCD diet caused significant changes of serum lipid composition in mice. Compared with MCD group, the sample distribution of HNK group showed no change.
To further clarify the differential lipids among each group, we analyzed the serum lipid composition of MCS, MCD and HNK groups using orthogonal partial least squares discriminant analysis (OPLS-DA). As shown in Fig. 3(b), MCS and MCD group were significantly separated. R2Y = 0.997, Q2 = 0.967, and both R2Y and Q2 were greater than 0.4, indicating that the OPLS-DA model was effective. Figure 3(c) showed that HNK group and MCD group were significantly differentiated, with R2Y = 0.998 and Q2 = 0.719, indicating that the OPLS-DA model was effective.
Variable importance in projection (VIP) (VIP ≥ 1) in OPLS-DA analysis and fold change (fold change ≥ 2 or ≤ 0.5) of univariate analysis were used to select differential lipid metabolites. As shown in Fig. 3(d), MCD group showed 448 different lipid metabolites compared with MCS group, including 24 fatty acyls (FA), 297 glycerophospholipids (GP), 49 sphingolipids (SP), 15 sterol lipid (ST), 63 glycerides (GL). Compared with MCD group, HNK group showed 34 different lipid metabolites, 18 of them (including 17 FA and 1 GP) were down-regulated, and 16 of them (including 2 ST, 10 GP, 1 FA, 1 prenol lipids, 1 SP and 1 GL) were up-regulated (Fig. 3(e)).
In order to intuitively observe the change of the relative content of different lipid metabolites, we standardized 34 screened significantly different lipid metabolites between MCD and HNK groups, and showed in the heat map (Fig. 3(f)). Compared with MCS group, 3 eicosanoids (12,13-EPOME, 9,10-DIHOME and 14(S)-HDHA) and 3 FFAs (16:2, 18:3 and 22:6) were significantly up-regulated, while TxB3, PA(18:0_22:6), SE(28:1_22:6) and 8 PEs (20:1_18:0, O-18:0_16:0, O-20:0_16:0, O-16:1_16:0, O-20:0_22:6, O-15:1_22:6, P-18:0_22:6 and P-16:0_16:0) were down-regulated in MCD group. HNK treatment significantly back-regulated the levels of these lipid metabolites. Furthermore, HNK treatment also increased the levels of taurolithocholic acid-3-sulfate (TLCA-3-S), PE(20:2_18:0), Coenzyme Q10, TG(O-20:0_16:0_18:1) and Cer(t18:0/24:0), and decreased the levels of PC(20:4_22:6) and 11 FFAs (14:0, 16:1, 17:1, 18:1, 19:1, 20:2, 20:3, 22:4, 22:5, 24:6 and 16:3), as compared to MCD mice.
KEGG pathway enrichment analysis was used to determine the metabolic pathways related to the 34 different lipids mentioned above. Figure 3(g) showed that fatty acid biosynthesis and unsaturated fatty acid biosynthesis pathways were related to HNK alleviating lipid metabolism disorders in MCD mice.
3.4. Effects of HNK on serum BA levels and hepatic mRNA expression of BA metabolism-related genes in NASH mice. The serum levels of 50 BAs were shown in Table S2 and the significant changes of serum BA levels in mice were shown in Table 1. Compared with MCS group, the levels of 23-nor-deoxycholic acid (23-DCA), hyocholic acid (HCA), ursocholic acid (UCA), hyodeoxycholic acid (HDCA), deoxycholic acid (DCA), 7-ketodeoxycholic acid (7-KDCA), 3β-ursodeoxycholic acid (3β-UDCA) and 3-oxodeoxycholic acid (3-oxo-DCA) were increased in MCD group. HNK treatment markedly decreased the levels of 23-DCA, HDCA, glycocholic acid (GCA) and taurodeoxycholic acid (TDCA), as compared with MCD group.
Table 1
Levels of differential bile acids in serum of the tested mice
BAs(ng/ml)
|
MCS
|
MCD
|
HNK
|
7-KDCA
|
717.72 ± 609.34
|
2875.19 ± 1749.38*
|
1924.01 ± 1137.45
|
DCA
|
223.43 ± 110.20
|
540.44 ± 301.67*
|
271.33 ± 96.50
|
TDCA
|
197.29 ± 229.70
|
308.34 ± 70.54
|
147.28 ± 29.11##
|
23-DCA
|
27.22 ± 17.66
|
277.30 ± 55.15**
|
97.73 ± 33.74##
|
HDCA
|
22.94 ± 10.95
|
44.70 ± 11.27**
|
24.53 ± 3.13##
|
UCA
|
19.27 ± 9.60
|
91.55 ± 58.73*
|
35.31 ± 23.46
|
3-oxo-DCA
|
13.85 ± 8.02
|
44.12 ± 29.71*
|
23.08 ± 17.75
|
HCA
|
13.50 ± 10.40
|
57.77 ± 41.79*
|
47.82 ± 26.85
|
GCA
|
6.43 ± 10.42
|
11.91 ± 1.42
|
6.18 ± 1.54##
|
3β-UDCA
|
3.30 ± 1.19
|
6.33 ± 2.97*
|
5.71 ± 2.72
|
(n = 7). **p < 0.01, *p < 0.05 vs MCS group; ##p < 0.01, #p < 0.05 vs MCD group. MCD, methionine- and choline-deficient diet; MCS, methionine- and choline-sufficient diet; HNK, honokiol.
The mRNA expression of BAs metabolism-related genes in mice liver were shown in Fig. 4. Compared with MCS group, the mRNA expression of CYP7A1, CYP27A1, Bsep, Mrp2, Ntcp and Oatp1b2 in MCD group were significantly decreased. HNK treatment markedly increased Oatp1b2 mRNA expression of MCD mice (Fig. 4f), whereas there were no significant effect on the mRNA expression of CYP7A1, CYP27A1, Bsep, Mrp2 and Ntcp (Fig. 4(a)-(e)).
3.5. Effects of HNK on the composition of GM in MCD mice. The Exponential dilution curve showed that the amount of 16Sr DNA sequencing data in this experiment was progressive and reasonable, and the sequencing depth was reliable (Fig. 5(a)). The statistical analysis of the Alpha diversity index of different samples under the 97% consistency threshold were shown in Table S3. There were no significant differences in Chao1, ACE, Simpson and Shannon indexes of GM among the three groups, indicating that HNK had no significant effect on the diversity of GM. Beta diversity analysis was applied to further study the composition of GM in each group, and non-metric multi-dimensional scaling (NMDS) showed significant differences in the composition of GM among the three groups (Fig. 5(b)).
The composition of GM of the tested mice at phylum, family and genus levels were shown in Fig. 5(c)-(e). Compared with MCS group, the abundance of Firmicutes and Turicibacter were markedly increased, and the abundance of Eggerthellaceae was significantly decreased in MCD group. MCD mice treated by HNK significantly increased the relative abundance of Ruminococcaceae, Caulobacteraceae, Micrococcaceae and Brevundimonas, and the abundance of Firmicutes and Dubosiella significantly decreased. Based on LEFSe analysis, species with LDA>4 were considered as the set value, and the microbiota with statistical differences among three groups were screened out and showed in Fig. 5(f).
3.6. Spearman correlation analysis. Spearman correlation analysis was used to explore the relationship between the abundance of GM and the contents of serum lipid or BAs significantly changed by HNK treatment in MCD mice. At phylum level, Firmicutes was negatively correlated with TLCA-3-S and TxB3 (Fig. 6(a)). At family level, Ruminococcaceae was negatively correlated with the levels of 14 FFAs (16:1、16:2、14:0、18:1、17:1、22:4、20:3、18:3、20:2、19:1、16:3、22:6、22:5、24:6) (Fig. 6b), Caulobacteraceae was negatively correlated with 9 FFAs (18:1、17:1、22:4、20:3、18:3、16:3、22:6、22:5 and 24:6) and positively related with 6 PEs (O-16:1_16:0、P-16:0_16:0、O-18:0_16:0、O-20:0_16:0、O-20:0_22:6 and P-18:0_22:6) (Fig. 6b), Micrococcaceae was positively related with Cer(t18:0/24:0), TLCA-3-S and TxB3 (Fig. 6(b)). At genus level, Dubosiella was positively related with 4 FFAs (18:1、17:1、22:4、18:3 and 24:6) and negatively correlated with multiple 4 PEs (20:1_18:0、20:2_18:0、P-16:0_16:0 and O-16:1_16:0),Brevundimonas was negatively correlated with 2 oxidized lipids (12, 13-EPOME and 9,10-DiHOME) and 9 FFAs (18:1、17:1、16:3、20:3、22:4、18:3、22:5、22:6 and 24:6) (Fig. 6(c)). At phylum level, Firmicutes was positively correlated with 23-DCA, 7-KDCA, DCA, UCA, 3-oxo-DCA, HCA, 3β-UDCA(Fig. 6(d)). At family level, Caulobacteraceae was negatively correlated with 23-DCA and TDCA, Ruminococcaceae was negatively correlated with 23-DCA(Fig. 6(e)). At genus level, Brevundimonas were negatively correlated with 23-DCA and TDCA, Dubosiella was positively related with 23-DCA, HCA, 7-KDCA, 3β-UDCA (Fig. 6(f)).