Specnuezhenide Ameliorates Age-Related Hepatic Lipid Accumulation via Modulating Bile Acid Homeostasis and Gut Microbiota in D-Galactose-Induced Mice

Age-related hepatic lipid accumulation has become a major health problem in the elderly population. Specnuezhenide (SPN) is a major active iridoid glycoside from an edible herb Fructus Ligustri Lucidi, which is commonly used for preventing age-related diseases. However, the beneficial effects of SPN on age-related liver injury remain unknown. This study aimed to reveal the effect of SPN on age-related hepatic lipid accumulation and the underlying mechanism. D-galactose (D-gal)-induced aging mice were treated with vehicle or SPN for 12 weeks. Treatment of SPN decreased lipid accumulation and inflammation in the liver of D-gal–induced mice. Untargeted and targeted metabolomics showed that the SPN could regulate the bile acid (BA) synthesis pathway and restore the BA compositions in serum, livers, and feces of the D-gal–induced mice. Furthermore, SPN enhanced the protein and mRNA levels of hepatic BAs synthesis enzymes cytochrome P45027A1, cytochrome P4507A1, cytochrome P4507B1, and cytochrome P4508B1. Meanwhile, SPN alleviated D-gal-induced gut dysbiosis and reversed the proportions of microbes associated with bile salt hydrolase activity, including Lactobacillus, Ruminiclostridium, and Butyrivibrio. Our study revealed that SPN attenuated age-related hepatic lipid accumulation by improving BA profiles via modulating hepatic BA synthesis enzymes and gut microbiota.


Introduction
As human life expectancy keeps increasing, aging-related health deterioration has been a growing challenge for clinical practices.Nonalcoholic fatty liver disease (NAFLD) is a major cause of chronic liver disease and has been prevalent in the elderly population [1].NLFLD is a liver metabolic disorder that results from increased total cholesterol (TC) and total triglyceride (TG) accumulation [2].Bile acids (BAs) are hepatic metabolites and have a close relationship with the pathophysiology of NAFLD.BAs have a well-established role in the regulation of cholesterol and lipid metabolism [3].Disordered BA homeostasis commonly accompanies aging [4], which increases the risk of aging-related hepatic lipid accumulation [5].Clinical studies proved that altered BA profiles underlie the development of NAFLD.Dietary BA supplementation can ameliorate high dietary lipid-induced liver injury by reducing lipid accumulation and improving liver injury [6].Thus, targeting BA metabolism has the potential to prevent hepatic lipid accumulation and NAFLD.Specnuezhenide (SPN), an iridoid glycoside, is the main component of Fructus Ligustri Lucidi, which is an edible herb commonly used to prevent aging-related diseases [7].The extract of Fructus Ligustri Lucidi could inhibit hepatic injury by carbon tetrachloride-induced mice [8], regulated fat balance in ovariectomized mice [9], and preserved bone mass in D-galactose (D-gal)/sodium nitrite-induced mice [10].Iridoid glycosides and SPN could reduce the development of liver injury in concanavalin A-induced mice [11] and carbon tetrachloride-induced mice [12], respectively.However, it is not clear whether SPN can attenuate age-related hepatic lipid accumulation and regulate BA homeostasis.
The current study aimed to describe the effect of SPN on age-related hepatic lipid accumulation and determine the mechanism of SPN with BA homeostasis in D-gal-induced mice.We established that SPN could reduce D-gal-induced hepatic lipid accumulation and regulate BA homeostasis.Moreover, SPN improved BA profiles by increasing the expressions of BA synthesis enzymes in the liver and modulating the gut microbiota community structure of D-gal-induced mice.
MS analysis was performed with a heated electrospray ionization source in both positive and negative ionization modes.The mass parameters were set as follows: auxiliary gas flow, 10 arb; sheath gas flow, 40 arb; spray voltage, 3.5 kV; m/z, 60-900; auxiliary gas heater temperature, 350 • C; and capillary temperature, 320 • C. For multivariate and univariable analysis, the metabolites were searched in the Human Metabolome Database (http://www.hmdb.ca/,accessed on 16 September 2022).MetaboAnalyst (http://www.metaboanalyst.ca/,accessed on 27 December 2022) was used to explore the key metabolic pathways represented by the differential metabolites.The multivariate data matrix was imported into SIMCA 14.0 (Umetrics, Umeå, Sweden) for visualization.

RNA Isolation and Real-Time Reverse Transcription Polymerase Chain Reaction (RT-PCR)
Total RNA was extracted using TRIzol (Invitrogen, Carlsbad, CA, USA) and reverse transcribed using a TOBOBlue qRT Premix with a gDNA Eraser 2.0 kit (#RTQ202; Toroivd Technology, Shanghai, China).RT-PCR (7500 RT-PCR; Applied Biosystems, Waltham, MA, USA) was used to analyze the levels of CYP7A1, CYP27A1, CYP7B1, and CYP8B1 using SYBR Green reagents (#QPS201; Toyobo, Osaka, Japan).The primer sequences are listed in Table S2.The relative mRNA expression was quantified by comparing cycle threshold values.GAPDH was used as a housekeeping gene.Data are shown as the fold change relative to controls.

16S rDNA Gene High-Throughput Sequencing
The total bacterial genomic DNA was extracted using a MagaBio Feces Genomic DNA Purification kit (#BSC48; Hangzhou Bioer Technology, Zhejiang, China) according to the manufacturer's protocol.The quality and quantity of the DNA were analyzed using a microvolume spectrophotometer (NanoDrop One; Thermo Fisher).DNA sequences were amplified using a PCR (S1000; Bio-Rad, Hercules, CA, USA) with primers F515/R806 with a dual indexing approach.A total of 50 µL of PCR reactions was prepared in duplicate and contained 50 ng DNA template, 200 nM primers (Invitrogen, Carlsbad, CA, USA), 2 × Premix Taq (#RR902A; TaKaRa, Dalian, China), and nuclease-free water.Reactions were carried out with an initial denaturation at 94 • C for 5 min, followed by 30 cycles of 30 s at 94 • C, 30 s at 52 • C, 30 s at 72 • C, and then at 72 • C for 10 min.The quality of the pooled sample was evaluated using a high-throughput nucleic acid protein analysis system (Houze Biological Technology, Hangzhou, China) and a Qubit 4.0 Fluorometer (Life Technologies, Grand Island, NY, USA).

Statistical Analysis
Statistical analysis was carried out using GraphPad Prism 7.0 (GraphPad, San Diego, CA, USA).Statistical analysis was conducted using one-way or two-way analysis of variance, followed by Dunnett corrections to compare multiple groups.Data were presented as means ± SD. p < 0.05 was considered statistically significant.

SPN Reduces Hepatic Lipid Accumulation in D-Gal-Induced Mice
After 12 weeks of treatment, there were no significant differences in body weight or liver index among the groups (Figure 1A,B).However, lobular inflammation and focal necrosis were observed in the histopathology of the liver in the MOD group (Figure 1C).SPN treatment could attenuate these injuries in the D-gal-induced liver.Furthermore, the hepatic lipid accumulation was examined by measuring TC, TG, LDL-C, and HDL-C in the liver (Figure 1D).TC, TG, and LDL-C contents in the MOD group increased compared with the control group (p < 0.01), and the content of HDL-C in the MOD group decreased (p < 0.01).SPN treatments decreased the levels of TC (p < 0.01), TG (p < 0.01), and LDL-C (p < 0.01) but increased the HDL-C contents in the D-gal-induced mice (p < 0.01), indicating that SPN could decrease hepatic lipid accumulation in D-gal-induced mice.

SPN Reverses BA Profile in Serum, Liver, and Feces in D-Gal-Induced Mice
Metabolomics analysis was carried out to investigate the mechanism underlying the effect of SPN on hepatic lipid accumulation.Untargeted metabolomics data showed that SPN groups were located closer to the CON group in both principal component analysis (Figure 2A, PCA) and orthogonal partial least-squares discriminant analysis (OPLS-DA).These data indicated that the serum metabolic profiles were reversed by SPN treatment.The metabolic pathway enrichment results indicated that primary BA biosynthesis was closely related to the effect of SPN (Figure 2B).

SPN Reverses BA Profile in Serum, Liver, and Feces in D-Gal-Induced Mice
Metabolomics analysis was carried out to investigate the mechanism underlying the effect of SPN on hepatic lipid accumulation.Untargeted metabolomics data showed that SPN groups were located closer to the CON group in both principal component analysis (Figure 2A, PCA) and orthogonal partial least-squares discriminant analysis (OPLS-DA).These data indicated that the serum metabolic profiles were reversed by SPN treatment.The metabolic pathway enrichment results indicated that primary BA biosynthesis was closely related to the effect of SPN (Figure 2B).BAs were quantified on the basis of the respective standard curves (Table S3).The total BAs decreased in the serum of the MOD group compared with the CON group (Figure 2C, p < 0.05).SPN treatment increased the total BAs in the serum of the D-gal mice (p < 0.05).Levels of total BAs in the liver and feces were also reduced in the MOD group (p < 0.01).SPN treatment increased the total BAs in the liver of the D-gal mice (p < 0.01).Compared with the MOD group, more BAs were excreted through the feces in the SPN-treated mice (p < 0.01).BAs were quantified on the basis of the respective standard curves (Table S3).T total BAs decreased in the serum of the MOD group compared with the CON group (Fi ure 2C, p < 0.05).SPN treatment increased the total BAs in the serum of the D-gal mice < 0.05).Levels of total BAs in the liver and feces were also reduced in the MOD group < 0.01).SPN treatment increased the total BAs in the liver of the D-gal mice (p < 0.0 Compared with the MOD group, more BAs were excreted through the feces in the SPN treated mice (p < 0.01).

SPN Enhanced Bile Acid Synthesis by Stimulating Hepatic Enzymes
Since primary BA synthesis occurs in the liver, we examined the effect of SPN on the hepatic enzymes for bile acid synthesis, including CYP7A1, CYP7B1, CYP8B1, and CYP27A1.IHC results showed that the protein levels of CYP7A1 (Figure 4A, p < 0.01), CYP7B1 (p < 0.01), CYP8B1 (p < 0.01), and CYP27A1 (p < 0.01) were downregulated in the MOD group compared to the CON group, while SPN treatment restored these changes (p < 0.01).The expressions of CYP7A1, CYP7B1, CYP8B1, and CYP27A1 in the liver were also analyzed using western blo ing (Figure 4B) and RT-PCR (Figure 5).The levels of these

SPN Enhanced Bile Acid Synthesis by Stimulating Hepatic Enzymes
Since primary BA synthesis occurs in the liver, we examined the effect of SPN on the hepatic enzymes for bile acid synthesis, including CYP7A1, CYP7B1, CYP8B1, and CYP27A1.IHC results showed that the protein levels of CYP7A1 (Figure 4A, p < 0.01), CYP7B1 (p < 0.01), CYP8B1 (p < 0.01), and CYP27A1 (p < 0.01) were downregulated in the MOD group compared to the CON group, while SPN treatment restored these changes (p < 0.01).The expressions of CYP7A1, CYP7B1, CYP8B1, and CYP27A1 in the liver were also analyzed using western blotting (Figure 4B) and RT-PCR (Figure 5).The levels of these enzymes decreased in the MOD group compared with the CON group (p < 0.01), which was consistent with a previous report [19].Importantly, the expression levels of CYP7A1 (p < 0.01), CYP7B1 (p < 0.01), CYP8B1 (p < 0.01), and CYP27A1 (p < 0.01) in the D-gal-induced livers were reversed in response to the SPN treatment.The above results suggested that SPN activated both the classical and alternative pathways of BA synthesis.

13, x FOR PEER REVIEW 8 of 14
was consistent with a previous report [19].Importantly, the expression levels of CYP7A1 (p < 0.01), CYP7B1 (p < 0.01), CYP8B1 (p < 0.01), and CYP27A1 (p < 0.01) in the D-galinduced livers were reversed in response to the SPN treatment.The above results suggested that SPN activated both the classical and alternative pathways of BA synthesis.

SPN Remodels the Gut Microbiota Genera Associated with Bile Salt Hydrolases
Given the critical crosstalk between the gut microbiota and BA metabolism, 16S rDNA high-throughput sequencing of the feces was used to evaluate the effect of SPN on gut microbiota.The Bray-Curtis PCoA revealed a total separation between the bacterial profiles of the CON, MOD, and SPN-H groups (Figure 6A).These groups shared 794 common OTUs (Figure 6B).Among the predicted functional metabolic pathways (Figure 6C), secondary BA biosynthesis and primary BA biosynthesis were affected by gut microbiota.
Bile salt hydrolase (BSH) was found in three phyla in D-gal-induced mice, including Firmicutes, Bacteroidetes, and Actinobacteria.D-gal injection significantly reduced Firmicutes (Figure 6D).There was no significant difference in the Firmicutes-to-Bacteroides ratio or αdiversity indices between groups (Figure S1).The relative abundance of the family Lactobacillaceae from Firmicutes displayed different profiles in the MOD and CON groups (Figure S2).Specifically, SPN increased the abundance of the BSH-related genera Lactobacillus from Lactobacillaceae (Figure S3, p < 0.05).Lactobacillus might contain BSH-T0, BSH-T2, and BSH-T3 [20].Additionally, SPN upregulated the abundance of Ruminiclostridium (p < 0.01) and Butyrivibrio (p < 0.01), which might contain BSH-T.These data indicated that SPN could affect BA compositions by regulating gut microbiota in D-gal mice.

SPN Remodels the Gut Microbiota Genera Associated with Bile Salt Hydrolases
Given the critical crosstalk between the gut microbiota and BA metabolism, 16S rDNA high-throughput sequencing of the feces was used to evaluate the effect of SPN on gut microbiota.The Bray-Curtis PCoA revealed a total separation between the bacterial profiles of the CON, MOD, and SPN-H groups (Figure 6A).These groups shared 794 common OTUs (Figure 6B).Among the predicted functional metabolic pathways (Figure 6C), secondary BA biosynthesis and primary BA biosynthesis were affected by gut microbiota.Bile salt hydrolase (BSH) was found in three phyla in D-gal-induced mice, including Firmicutes, Bacteroidetes, and Actinobacteria.D-gal injection significantly reduced Firmicutes (Figure 6D).There was no significant difference in the Firmicutes-to-Bacteroides ratio or α-diversity indices between groups (Figure S1).The relative abundance of the family Lactobacillaceae from Firmicutes displayed different profiles in the MOD and CON (Figure S2).Specifically, SPN increased the abundance of the BSH-related genera Lactobacillus from Lactobacillaceae (Figure S3, p < 0.05).Lactobacillus might contain BSH-T0, BSH-T2, and BSH-T3 [20].Additionally, SPN upregulated the abundance of Ruminiclostridium (p < 0.01) and Butyrivibrio (p < 0.01), which might contain BSH-T.These data indicated that SPN could affect BA compositions by regulating gut microbiota in D-gal mice.

Discussion
Hepatic lipid accumulation implies injury to the liver and most likely causes NAFLD in elderly populations.Recently, natural hepatoprotective compounds have attracted increasing attention.SPN is a major active component from an edible herb Fructus Ligustri Lucidi, which can tonify the liver.However, the protective effects of SPN against agerelated liver injury remain unclear.In this study, SPN treatment could attenuate the hepatic lipid accumulation and inflammation in D-gal-induced mice.Untargeted and targeted metabolomics analysis proved that the BA synthesis pathway was crucial for the protection effects of SPN.SPN treatment restored BA compositions in the serum, liver, and feces of D-gal-induced mice.Furthermore, we found that SPN could upregulate the expressions of the hepatic BAs synthesis enzymes and modulate the gut microbiota.
Mass spectrometry is a powerful technique for investigating metabolic features associated with aging [21].A non-targeted analysis of human blood metabolome disclosed a close link between metabolites and aging [22].Consistent with this, our untargeted metabolomics revealed that shifts in metabolomics profiles were observed in D-gal-induced mice.Treatment with SPN moderately corrected this metabolic change.An altered BA profile is tightly related to the development of age-related liver disease [23].The enrichment analysis showed that primary BA biosynthesis was the most significant metabolic pathway in response to the SPN treatment.BA metabolism is a promising target for developing antisenescencebased therapies.We established a method to simultaneously quantify 21 BAs in the serum, liver, and feces.Targeted-BA metabolomics demonstrated that changes in serum BA contents appeared in response to D-gal injection.SPN restored six BAs in the D-gal-induced mice.Our data showed that SPN treatment increased CDCA in the D-gal-induced mice.Additionally, NAFLD is characterized by increased DCA and LCA [24].Excessive DCA produces reactive oxygen species that cause DNA damage and senescence-associated secretory phenotypes, inducing a pro-inflammatory response in steatosis hepatocytes [25], which facilitates the development of aging-related liver disease [26].Additionally, LCA shows a deleterious effect on the liver [27], but these effects can be neutralized by UDCA [28].In our study, SPN reduced the contents of DCA and LCA in the serum, liver, and feces of the D-gal-induced mice.These results indicated that SPN may attenuate D-gal-induced hepatic lipid accumulation and injury via two pathways: (1) increasing the CDCA and UDCA content, which is positively related to cholesterol absorption [29]; (2) reducing the DCA level, which contributes to hepatotoxicity in the mice [30].
BAs are synthesized as primary BAs in the liver through enzymatic pathways [31].Cholesterol is hydroxylated in the 7α-position by CYP7A1 or converted to an oxysterol before being in the 7α-position by CYP7B1 and arises in the liver by CYP27A.However, liver function deteriorates with age, leading to downregulated BA synthesis enzymes and reduced BA amount [32].In our study, D-gal injection reduced the expression of CYP7A1, CYP27A1, CYP7B1, and CYP8B1 in mouse liver.The pathophysiological changes of these BA synthesis enzymes are risk factors for the development of lipid accumulation and cholesterol metabolism [33].Iridoid glycosides are an important source of BA synthesis modulators.For example, swertiamarin [34] and gentiopicroside [35] showed hepatoprotective effects in alpha-naphthylisothiocyanate-induced mice via upregulating CYP8B1 and CYP27A1.Our data showed that SPN enhanced the expression of these four key BAs enzymes (CYP7A1, CYP27A1, CYP7B1, and CYP8B1) in the liver.These findings suggest that SPN promoted BA synthesis by stimulating the hepatic BA synthesis enzymes in D-gal-induced mice.
Furthermore, primary BAs are modified in the intestine through the actions of bacteria, such as deconjugation via BSH activity [36].Since the abundance and diversity of the gut microbiota in mice are changed upon senescence [37], the BA composition also changes during aging.In our study, the age-related model is commonly established by D-gal injection.The action of D-gal-induced models is largely attributed to changes in gut microbial populations [38].In the present study, D-gal injection changed gut microbiota compositions with BSHs, which has essential effects on secondary BA production.A very intriguing finding was that SPN could reverse three BSH-related genera, including Lactobacillus, Ruminiclostridium, and Butyrivibrio.The decreased level of Lactobacillus results in increased levels of Tauro-β-MCA and a substrate of BSHs [39].Ruminiclostridium and Butyrivibrio might contain BSH-T1, and BSH-T1 exerted vigorous enzyme activity when the substrates were TCDCA and GCDCA, resulting in a significant change in CDCA content [40].Our data indicated that SPN regulated the relative abundance of BAs at least partially by increasing the gut microbes that possess BSHs.

Conclusions
In this work, SPN, one of the iridoid glucoside representatives, the most potential active component in Fructus Ligustri Lucidi, could attenuate D-gal-induced hepatic lipid accumulation and injury through regulating the expression of BA synthesis enzymes and the structure of the gut microbiota.Our study provides a theoretical basis that SPN may be a functional factor supplement for treating age-related hepatic lipid accumulation via modulation of BA homeostasis.Future work may identify the toxicity and the clinical application of this compound.

Figure 2 .
Figure 2. Specnuezhenide (SPN) regulated the metabolic profile and bile acid (BA) composition (A) The plot depicts separation of the principal component analysis (PCA) and orthogonal partial lea squares discriminant analysis (OPLS−DA) of serum from different groups at positive and negati modes.(B) Metabolic pathway analysis of identified potential marker.(C) Untargeted metabolom analyses of serum, liver, and feces were performed to measure the concentration of total BAs (n = Data are expressed as mean ± SD. * p < 0.05, ** p < 0.01 compared with the model (MOD) group.

Figure 2 .
Figure 2. Specnuezhenide (SPN) regulated the metabolic profile and bile acid (BA) compositions.(A) The plot depicts separation of the principal component analysis (PCA) and orthogonal partial least-squares discriminant analysis (OPLS−DA) of serum from different groups at positive and negative modes.(B) Metabolic pathway analysis of identified potential marker.(C) Untargeted metabolomic analyses of serum, liver, and feces were performed to measure the concentration of total BAs (n = 5).Data are expressed as mean ± SD. * p < 0.05, ** p < 0.01 compared with the model (MOD) group.

Figure 6 .
Figure 6.Specnuezhenide (SPN) treatment improved gut microbiota disbalance.(A) Principal correlation coefficient analysis (PCoA) analysis of gut microbiota at operational taxonomic units (OTUSA) level based on Bray−Curtis.(B) Shared and unique OTUs among the different groups.(C) The predicted functional metabolic pathways.(D) The relative abundance of bile salt hydrolase (BSH) related phylum (n = 5).Data are expressed as mean ± SD. * p < 0.05 compared with the MOD group.