Next Article in Journal
Comparative Analysis of the Codon Usage Pattern in the Chloroplast Genomes of Gnetales Species
Previous Article in Journal
MicroRNA Inhibiting Atheroprotective Proteins in Patients with Unstable Angina Comparing to Chronic Coronary Syndrome
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ursolic Acid Regulates Immune Balance, Modulates Gut Microbial Metabolism, and Improves Liver Health in Mice

1
College of Animal Science and Technology, Hebei Agricultural University, Baoding 071051, China
2
College of Traditional Chinese Veterinary Medicine, Hebei Agricultural University, Baoding 071051, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2024, 25(19), 10623; https://doi.org/10.3390/ijms251910623
Submission received: 18 August 2024 / Revised: 20 September 2024 / Accepted: 23 September 2024 / Published: 2 October 2024

Abstract

:
Ursolic acid (UA) has demonstrated significant immunomodulatory and hepatoprotective effects; however, the underlying mechanisms remain unclear. This study aims to analyze the impact of UA on the gut microbiome, metabolome, and liver transcriptome, investigate UA’s role in maintaining gut immune homeostasis and liver health, and evaluate the potential contributions of gut microbes and their metabolites to these beneficial effects. Our findings indicate that UA enhances immune balance in the jejunum, fortifies intestinal barrier function, and promotes overall gut health. UA modulates the intestinal microbiota and its metabolic processes, notably increasing the abundance of beneficial bacteria such as Odoribacter and Parabacteroides, along with their metabolites, including ornithine and lactucin. Additionally, UA inhibits the expression of interleukin-1 receptor 1 (IL1R1) and calcium (Ca2+) voltage-gated channel auxiliary subunit beta 2 (CACNB2) while enhancing the synthesis pathways of retinol and ascorbic acid, thereby exerting a protective influence on liver function. In summary, UA enhances intestinal immune homeostasis and promotes liver health, with these advantageous effects potentially mediated by beneficial bacteria (Odoribacter and Parabacteroides) and their metabolites (ornithine and lactucin).

Graphical Abstract

1. Introduction

Body health is closely related to gut health, relying on moderate immunity, an intact gut barrier, and a stable gut microbe. The gut harbors a diverse array of immune cells, including T helper (Th) cells and T regulatory (Treg) cells, which work together to maintain the balance of the body’s immune response [1]. An overproduction of immune-related cytokines, such as interleukin (IL)-6, IL-1β, and IL-17A, can lead to intestinal cell damage and compromise the integrity of the gut barrier, consequently triggering systemic inflammation. In contrast, cytokines such as IL-10 and transforming growth factor-β (TGF-β) play a role in inhibiting inflammation [1]. Ursolic acid (UA) is a naturally occurring pentacyclic triterpene carboxylic acid found in various plants, either in free acid form or as an aglycone of triterpene saponins [2]. In recent years, there has been a surge of interest in UA, primarily concerning its pharmacological properties, which include antioxidant [3], antitumor [4], anti-inflammatory [5], and antimicrobial [2] effects. Our previous meta-analysis confirmed that UA can reduce levels of pro-inflammatory cytokines like IL-1β and IL-6 [6]. Furthermore, UA has been shown to activate the immune system and enhance the proliferation of immune cells [7]. Nevertheless, the regulatory mechanisms by which UA influences intestinal immune homeostasis and gut health remain unclear.
The gut microbiota plays a critical role in maintaining intestinal immune homeostasis and barrier function [8,9,10]. Research has demonstrated that gut microbes can modulate immune equilibrium through the regulation of Th17/Treg cell dynamics. For example, Lactobacillus rhamnosus GG is capable of modulating Th17/Treg homeostasis, thus aiding in the preservation of host intestinal health [11]. Similarly, Bacteroides fragilis enhances Treg cell activation by increasing Foxp3 expression [12]. UA has exhibited significant regulatory effects on the gut microbiota [13,14,15]. Consequently, understanding how UA influences immune balance and gut health via the modulation of gut microbes presents a compelling area of investigation. Additionally, UA has been shown to mitigate liver damage by modulating enterohepatic axis homeostasis [13] and is recognized as a contributing factor to the hepatoprotective properties of various foods and medicinal plants [13,16]. The gut microbiota, a key component of the gut–liver axis (GLA), is also closely associated with liver health [8]. Studies indicate that Lactobacillus and Bifidobacterium can enhance metabolic processes in rats through the GLA [17]. However, the specific mechanisms through which UA exerts protective effects on the liver necessitate further exploration.
To address these questions, this study utilizes a multi-omics approach to analyze the effects of UA on the intestinal microbiota, metabolome, and liver transcriptome in mice. We evaluate the role of UA in immune balance and liver protection, elucidating the potential contributions of gut microbes and their metabolites to these processes.

2. Results

2.1. Effect of UA on Immune-Related Cytokines in Mice

As shown in Figure 1A, the body weight (BW) of mice in the UA 5 and UA 25 groups was significantly higher than that of the CON group (p < 0.05). Compared to the CON group, the UA 5 and UA 25 groups had significantly lower levels of interleukin (IL)-1β, IL-6, and tumor necrosis factor (TNF)-α (p < 0.05) and higher levels of transforming growth factor (TGF)-β and IL-10 (p < 0.05) (Figure 1B–G).

2.2. UA Promotes Intestinal Development, Intestinal Epithelial Barrier Integrity, and Immune Balance in Mice

Mouse intestinal tissues (duodenum, jejunum, and ileum) were stained with H&E to explore the effect of UA on intestinal development in mice (Figure 2A–F). The results indicated that UA treatment (UA 5 and UA 100) significantly increased duodenal villus length (V), while (UA 5, UA 25, and UA 100) significantly increased duodenal crypt depth (C), (UA 5, UA 25, UA 100, and UA 250) significantly increased duodenal wall thickness (p < 0.05), and (UA 5, UA 25, UA 100, and UA 250) significantly increased duodenal crypt depth. However, there was no significant effect on the V/C ratio (p > 0.05). Compared to the CON group, UA treatment (UA 5, UA 25, UA 100, and UA 250) significantly increased jejunal villus length, (UA 25) crypt depth, and (UA 25) V/C and (UA 5 and UA 250) intestinal wall thickness (p < 0.05). UA treatment had no significant effect on the intestinal morphology (villus length, crypt depth, V/C, and intestinal wall thickness) of the ileum (p > 0.05). It is concluded that the effect of UA on the small intestine is mainly concentrated in the jejunum, with the optimal dose being the UA 25 group. Based on this, we further investigated the effects of optimal UA treatment (UA 25) on the intestinal barrier and immune balance in mice. The results showed that UA treatment significantly increased the mRNA expression of zonula occludens-1 (ZO-1) and Occludin in the jejunum of mice (p < 0.05) but had no significant effect on the mRNA expression of Claudin1 (p > 0.05) (Figure 2G–I). Compared with the CON group, the UA-treated group had significantly higher forkhead box P3 (Foxp3) mRNA expression levels and significantly lower signal transducer and activator of transcription 3 (Stat3) mRNA expression levels (p < 0.05). UA had no significant effect on retinoid-related orphan receptor gamma t (RORγt) mRNA expression (p > 0.05) (Figure 2J–L) based on the optimal effects of the UA 25 group on body weight, serum immune factor levels, and intestinal tissue morphology, as well as beneficial effects on genes related to the intestinal barrier and immune balance. We selected the UA 25 group (UA group) for the omics study.

2.3. UA Changes the Composition and Structure of the Gut Microbiota in Mice

To explore whether the physiological effects of UA were related to the gut microbiota, we further investigated the effects of UA on the gut microbiota. Figure S1A–C illustrates the rank abundance curves, rarefaction curves, and Shannon curves derived from this gene sequencing. There was no significant difference in the Shannon index between the two groups (p > 0.05) (Figure 3A). The results of principal coordinate analysis (PCoA) indicated that the first principal component explained 42.32% of the variance between groups, and the second principal component explained 25.31% of the variance between groups. The results of the Adonis analysis demonstrated that UA treatment significantly regulated the structure of intestinal flora in mice (p < 0.05) (Figure 3B). Figure S1D,E presents the top 10 bacteria with relative abundance at the phylum level and the top 15 bacteria with relative abundance at the genus level. UA significantly reduced the relative abundance of Romboutsia and significantly increased the relative abundance of Odoribacter, Coriobacteriaceae_UCG-002, Parabacteroides, Anaeroplasma, and Norank_f__Eggerthellaceae (p < 0.05) (Figure 3C). Linear discriminant analysis (LDA) effect size (LEfSe) was employed to identify the core taxa (LDA ≥ 3, p < 0.05) that were most likely to explain the differences between groups. The results revealed that there were six taxa-marker bacteria in the CON group and nine taxa-marker bacteria in the UA group (Figure 3D and Figure S1F). Notably, Romboutsia was further identified as the marker bacteria in the CON group, while Odoribacter, Coriobacteriaceae_UCG-002, and Parabacteroides were identified as the marker bacteria in the UA group.

2.4. UA Also Altered Metabolism of Gut Microbiota

A total of 355 metabolites were identified in the mouse fecal metabolome. Partial least squares discriminant analysis (PLS-DA) indicated that the mice in the CON group and the UA group exhibited different metabolic characteristics (Figure S2A,B). Sixteen metabolites were significantly different between the two groups, among which indole carboxylic acid sulfate and cholic acid were significantly enriched in the CON group. Dihydroconiferyl alcohol, 5′-phosphoribosyl-N-formylglycinamide, 5-hydroperoxypent-1-enylbenzene, Asn Asn Leu Asn Val, 3-hydroxy-N6,N6,N6-trimethyl-L-lysine, ornithine, lactucin, D-erythro-sphingosine C-17, dodemorph, wyerone, neoporrigenin B, Xi-4-hydroxy-4-methyl-2-cyclohexen-1-one, prolylglycine, and PA(24:0/22:6(5Z,8E,10Z,13Z,15E,19Z)-2OH(7S, 17S)) were significantly enriched in the UA group (variable importance in projection (VIP) > 1, p < 0.05) (Figure 4A). Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis based on the 16 differential metabolites showed that three pathways were significantly enriched. It includes D-amino acid metabolism, arginine biosynthesis, and biosynthesis of various other secondary metabolites (Benjamini–Hochberg false discovery rate (FDR) < 0.05) (Figure 4B).

2.5. UA Treatment Alters the Gut Microbiota and Its Metabolic Characteristics and Affects the Host Phenotype

Subsequently, we evaluated the correlation among gut microbiota, serum metabolites, and host phenotype to further explore the possible regulatory mechanisms between gut microbiota characteristics and host phenotype following UA treatment. As indicated by spearman’s correlation coefficient, the four bacteria identified as marker bacteria were significantly correlated with 13 differential metabolites. These metabolites were further associated with BW, TNF-α, IL-6, IL-1β, IL-10, and TGF-β (representing the physiological features of mice) (|R| > 0.6, p < 0.05) (Figure 4C). Odoribacter, Coriobacteriaceae_UCG-002, and Parabacteroides, which were identified as the marker bacteria of the UA group, were positively correlated with metabolites (lactucin, ornithine, PA(24:0/22:6(5Z,8E,10Z,13Z,15E,19Z)-2OH(7S, 17S), prolylglycine, and wyerone) enriched in the UA group and were negatively correlated with Indole carboxylic acid sulfate enriched in the CON group. Romboutsia, which was identified as the marker bacteria in the CON group, was negatively correlated with 5-hydroperoxypent-1-enylbenzene, Asn Asn Leu Asn Val, dihydroconiferyl alcohol, neoporrigenin B, ornithine, and PA(24:0/22:6(5Z,8E,10Z,13Z,15E,19Z)-2OH(7S, 17S). Prolylglycin and Xi-4-hydroxy-4-methyl-2-cyclohexen-1-one) enriched in the UA group were positively correlated with indole carboxylic acid sulfate enriched in the CON group. Lactucin was negatively correlated with IL-6 and IL-1β and positively correlated with BW. Ornithine was negatively correlated with TNF-α, IL-6, and IL-1β, and positively correlated with BW and IL-10. In addition, there was a direct correlation between the differential bacteria and the host phenotype. Odoribacter and Parabacteroides were significantly negatively correlated with TNF-α, IL-6, and IL-1β and significantly positively correlated with BW, TGF-β, and IL-10 (Figure S3). Coriobacteriaceae_UCG-002 was negatively correlated with IL-6 and IL-1β (p < 0.05) but positively correlated with IL-10 (p < 0.05). Romboutsia was significantly positively correlated with IL-6 and negatively correlated with IL-10.

2.6. UA Alters the Liver Transcriptome

To investigate the effects of UA on the mouse liver, we conducted RNA-seq analysis. A total of 29,094 genes were identified in the mouse liver. The results of the PCA analysis indicated a clear separation between the two groups of samples. The first principal component explained 48.12% of the variance between the two groups, and the second principal component explained 15.22% of the variance (Figure 5A). According to the filter threshold |log2FC| ≥ 1 and FDR < 0.05, we identified 636 differentially expressed genes (DEGs) in the UA vs. CON group. Compared to the CON group, there were 167 up-regulated DEGs and 469 down-regulated DEGs in the UA group (Figure 5B). In addition, the expression results of eight random mRNAs detected by qRT-PCR are presented (Figure S4A). The results of RNA-Seq (log2 fold change) and qRT-PCR (log2 fold change) were significantly positively correlated (correlation value = 0.84, p < 0.01) (Figure S4B). This indicates the high accuracy of both the sequencing and differential expression analysis results. Gene Ontology (GO) enrichment analysis was performed separately for down-regulated and up-regulated DEGs, and the results showed that the down-regulated genes were significantly enriched to 47 GO terms, such as cell projection organization, chemical synaptic transmission, anterograde trans-synaptic signaling, neuromuscular junction development, centrosome, and triplet codon-amino acid adaptor activity (Figure S5A). The up-regulated genes were significantly enriched in 191 GO terms, including myofilament, skeletal muscle thin filament assembly, and FATZ binding (Figure S5B).
The up-regulated and down-regulated genes were further subjected to KEGG pathway analysis to determine the functional categories. The results showed that there were no significantly enriched pathways among the down-regulated genes, but the pathway with the most enriched genes was the MAPK signaling pathway (DEGs = 11, FDR > 0.05) (Figure 5C). The top five significantly enriched pathways of up-regulated genes were retinol metabolism (DEGs = 6), ascorbate and aldarate metabolism (DEGs = 4), oxidative phosphorylation (DEGs = 7), chemical carcinogenesis-reactive oxygen species (DEGs = 9), and arginine and proline metabolism (DEGs = 5) (FDR < 0.01) (Figure 5D). Based on all DEGs, we further constructed the regulatory mechanism of UA on the ko04010, ko04668, ko00053, and ko00190 pathways (Figure 5E). Interestingly, all DEGs involved in the ko04010 and ko04668 pathways were down-regulated genes (Figure S6A), and all DEGs involved in the ko00053 (Figure S6B) and ko00190 (Figure S6C) pathways were up-regulated genes.

2.7. Multi-Omics Integrated Analysis Reveals the Possible Regulatory Mechanism of UA

Connections between bacteria, metabolites, liver genes, gut genes, and host phenotypes were established through an association study (Spearman correlation coefficient, |R| > 0.6, p < 0.05) (Figure 6). BW was positively correlated with intestinal barrier-related genes (Occludin and ZO-1) and negatively correlated with inflammation-related genes (Mitogen-activated protein kinase (MAPK)13, MAPK11, Interleukin-1 receptor 1 (IL1R1), Neurotrophic factor-3 (NTF3), and Calcium (Ca2+) voltage-gated channel auxiliary subunit beta 2 (CACNB2)). The immune-related cytokines (IL-1β, IL-6, and TNF-α) were negatively correlated with the intestinal barrier-related genes (Occludin and ZO-1), and all of them were positively correlated with the inflammation-related gene (CACNB2). In addition, IL-1β and TNF-α exhibited a strong positive correlation with IL1R1, of which IL-1β had the strongest correlation with IL1R1 (R = −0.799). Stat3, an immune balance-related gene, inhibited the expression of TGF-β and promoted the expression of IL-17A, while Foxp3 facilitated the expression of TGF-β. Odoribacter and Parabacteroides were positively correlated with Foxp3 and Parabacteroides was negatively correlated with Stat3. Coriobacteriaceae_UCG-002, Odoribacter, Parabacteroides, wyerone, ornithine, and lactucin could enhance intestinal barrier function by promoting the expression of Occludin and ZO-1. Odoribacter, wyerone, ornithine, and dodemorph showed significant negative correlations with IL1R1 and CACNB2. Notably, Odoribacter (16 correlations: 3 negative and 13 positive), Parabacteroides (15 correlations: 3 negative and 12 positive), ornithine (10 correlations: 4 negative and 6 positive), wyerone (8 correlations: 4 negative and 4 positive), and lactucin (7 correlations: 4 negative and 3 positive) had the highest number of correlations, suggesting that they may serve as crucial connections in the anti-inflammatory pathway of UA.

3. Discussion

UA has garnered considerable attention for its array of pharmacological benefits [6,15,18,19]. In this study, the significant effects of UA in the jejunum were observed, including a notable increase in the V/C ratio and enhancement of intestinal barrier function. These findings highlight UA’s role in promoting digestion and absorption as well as in improving gut health in mice. Th17 and Treg cells, which differentiate from CD4+ T cells, are pivotal in the immune response and synergistically contribute to intestinal immune function [20]. The balance between Th17 and Treg cells is crucial for maintaining intestinal immune homeostasis [21]. Disruption of this balance can adversely affect nutrient digestion and absorption, thus compromising overall health [1,21,22].
IL-6-mediated Stat3 signaling stimulates the differentiation of Th17 cells by inducing the nuclear receptor RORγt, leading to the production of the inflammatory cytokine IL-17A [23]. Subsequently, IL-17A can enhance the secretion of other pro-inflammatory cytokines, including TNF-α, IL-6, and IL-1β [24]. Conversely, TGF-β promotes the transcription of Foxp3, facilitating the differentiation of Treg cells and the production of anti-inflammatory cytokines such as IL-10 and TGF-β [25]. Moreover, Foxp3 can inhibit the differentiation of Th17 cells by restricting their transcriptional potential [23]. Our results indicate that UA treatment resulted in the down-regulation of Stat3 mRNA expression, up-regulation of Foxp3 mRNA expression, and increased levels of TGF-β and IL-10 in mice. These findings suggest that UA can improve intestinal immune homeostasis.
The gut microbiota serves as a critical link between environmental factors and host health [26,27]. Investigating whether UA’s functional role is connected to intestinal microbiota composition has become a topic of interest [13,15,28]. In our study, UA treatment significantly altered the composition and metabolites of the gut microbiota. Notably, the abundance of Odoribacter and Parabacteroides was significantly increased in the UA-treated mice. Odoribacter, a potential probiotic, is known to enhance mucosal immunity and maintain intestinal homeostasis [29]. Its colonization can augment Foxp3/RORγt regulatory T cells (Treg cells) and induce IL-10 production, thereby mitigating colitis in mice [30,31]. Numerous studies have highlighted the immunomodulatory effects of Parabacteroides [32,33,34], which have also been shown to promote CD4+ T cell differentiation toward the CD4/Foxp3/IL-10 regulatory pathway [35,36]. Our correlation analysis demonstrated associations between Odoribacter, Parabacteroides, and key immune factors, including Foxp3, Occludin, ZO-1, IL-1β, IL-6, TNF-α, and IL-10. These observations suggest that UA’s effects on immune balance and intestinal barrier function may be mediated by Odoribacter and Parabacteroides.
Furthermore, gut microbiota-derived ornithine has been shown to stimulate RORγt (+) IL-22 (+) cells and promote gut mucus formation, playing a vital role in maintaining intestinal homeostasis [37]. Previous research has demonstrated that lactucin can modulate gut microbiota and restore intestinal barrier function in inflammatory rat models [38]. Our findings revealed that both ornithine and lactucin were significantly enriched in UA-treated mice, indicating that these metabolites may be key factors in UA-mediated immune balance and intestinal barrier enhancement.
Additionally, our study investigates the effects of UA on liver health in mice. The MAPK signaling pathway is integral to regulating immune responses and plays a critical role in innate immunity [39]. Genes such as IL1R1, MAPK11, and MAPK13, expressed by immune or inflammatory cells, are crucial for inflammatory signaling [40,41]. Ras-guanine nucleotide-releasing protein 1 (RasGRP1) serves as a key regulator of the inflammatory response and is implicated in various inflammatory diseases [42]. Calcium (Ca2+) voltage-gated channel auxiliary subunit beta 2 (CACNB2) may contribute to aberrant RAS-MAPK activation [43]. Our study found that IL1R1, MAPK11, MAPK13, CACNB2, and RasGRP1 expression levels were significantly downregulated in the livers of mice in the UA group, suggesting that UA may lower the risk of liver inflammation and improve liver homeostasis. Moreover, the upregulated genes in the UA group were significantly enriched in pathways related to retinol metabolism, ascorbic acid and aldarate metabolism, and oxidative phosphorylation. Retinol, a form of vitamin A, plays an unexpected but crucial role in regulating immune responses [44]. Ascorbate, or vitamin C, acts as a potent antioxidant that supports various immune functions and enhances immune defenses [45]. Mitochondrial oxidative phosphorylation is responsible for generating the majority of ATP, which is essential for sustaining life and maintaining metabolic homeostasis [46]. These results imply that UA may enhance liver health by promoting the synthesis of retinol and ascorbic acid, thereby regulating metabolic homeostasis in the liver.
Interestingly, we observed that the IL1R1 and CACNB2 genes, which are associated with the MAPK signaling pathway, are broadly linked to immune-related factors, bacteria (Odoribacter and Parabacteroides), and metabolites (ornithine and lactucin). Parabacteroides distasonis has been shown to reduce liver damage and inflammation by inhibiting the NF-κB/MAPK pathway and activating the Nrf2 pathway, which may be related to the improvement in amino acid and bile acid metabolism in intestinal microorganisms [47]. Parabacteroides distasonis has a protective effect against non-alcoholic steatohepatitis (NASH) in mice by producing beneficial metabolites from dietary fiber [48]. In addition, Odoribacter has also been shown to be significantly negatively associated with measures of NASH progression [49]. Previous studies suggest that the liver-protective mechanism of ornithine may involve its conversion into the antioxidant glutathione and the nitric oxide precursor L-arginine, both of which are crucial for improving liver microcirculation [50,51]. Additionally, lactucin has been shown to downregulate the MAPK pathway, further enhancing liver health [52,53]. Lactucin may also enhance GLA circulation and reduce liver inflammation [38]. These observations suggest that the effects of UA on liver health may be mediated by regulating gut microbiome composition and its metabolism. However, our correlation analysis offers only indicative and predictive insights and does not definitively establish causal relationships between metabolites and hepatoprotective effects. Further studies are required to validate the relationship between these metabolites and liver health.
Our data can help us understand the function of UA in regulating immune homeostasis and improving liver health and identify the key role that gut microbiota and its metabolites may play in this process. However, we acknowledge that some limitations should be noted. We did not assess changes in gut microbial function and its association with inflammation in UA-treated mice. Further systematic studies, such as metagenomic, metatranscriptomic, or other multi-omics investigations, are needed to identify the function of gut microbiota and further explain the impact of changes in microbial metabolic function on inflammation.
In conclusion, UA can improve immune balance and enhance intestinal barrier function in mice. UA is able to increase beneficial bacteria (Odoribacter and Parabacteroides) and their metabolites (ornithin and lactucin). UA can inhibit the expression of IL1R1 and CACNB2 mRNA in the liver, upregulate the synthesis pathway of retinol and ascorbic acid in the liver, and promote liver health.

4. Materials and Methods

4.1. Ethics Statement

All experiments were performed in accordance with the ARRIVE guidelines and were carried out in accordance with the U.K. Animals (Scientific Procedures) Act, 1986 and associated guidelines, EU Directive 2010/63/EU for animal experiments, and approved by the Ethics Committee of Animal Experimentation of Hebei Agricultural University (Protocol 2023058).

4.2. Animals and Experimental Design

Eight-week-old female C57BL/6J mice were purchased from SPF Biotechnology Co., Ltd. (Beijing, China) and housed in pathogen-free facilities with a 12-h light and 12-h dark cycle at 22 °C.
After acclimatization to the laboratory conditions for 1 week, mice were randomly divided into 5 groups (n = 6 per group), including the control (CON) group, the UA 5, UA 25, UA 100, and UA 250 groups. Mice in the CON group were provided with a daily oral administration of 100 μL of PBS. Mice from the UA 5, UA 25, UA 100, and UA 250 groups received 5, 25, 100, and 250 mg/kg body weight of UA (purity ≥ 98%, Shenniu Pharmaceutical Co., Ltd., Dezhou, China; Figure S7) added to 100 μL PBS by daily gavage for 21 days. All group mice were administered menstruum, which was isovolumetric with 0.5% sodium carboxymethyl cellulose and 1% Tween 80. UA was dissolved in that menstruum (the menstruum and UA were administered directly into the stomach using an oral gavage needle).
Body weight was measured every three days throughout the study. On the 21st day, the mice were anesthetized with ether, the eyeballs were removed, blood samples were taken in 1.5 mL centrifuge tubes, and serum was separated by centrifugation for 15 min (3500× g, 4 °C), after which the mice were sacrificed. Fecal samples from all mice were collected aseptically and stored at −80 °C for microbiome and metabolomic analyses. The middle segments (0.5 cm) of the duodenum, jejunum, and ileum were flushed with PBS and fixed in a 10% formaldehyde solution for histological observation. Jejunum and liver tissue were washed with PBS, cut into slices, snap-frozen in liquid nitrogen, and stored at −80 °C for future analysis.

4.3. Enzyme-Linked Immunosorbent Assay (ELISA)

Serum samples were measured for the levels of IL-1β (MM-0040M2), IL-6 (MM-0163M2), TNF-α (MM-0132M2), TGF-β (MM-0689M2), IL-10 (MM-0176M2), and IL-17A (MM-0759M2) according to the instructions of specific ELISA kits (Jiangsu Meimian Industrial Co., Ltd., Yancheng, China).

4.4. Histological Analysis

The fixed duodenum, jejunum, and ileum samples were embedded in paraffin and stained with hematoxylin–eosin to obtain cross-sectional sections. In each section, intact and representative villi were selected, and intestinal morphology was measured. The images were evaluated using the ImageJ software 1.53j (US National Institutes of Health, Bethesda, MD, USA). Villus height (VH), crypt depth (CD), and intestinal wall thickness were measured. The ratio of villus height to crypt depth (V/C) was calculated.

4.5. Quantitative Reverse Transcription-Polymerase Chain Reaction (qRT-PCR)

Total RNA extraction and quantitative real-time PCR were performed as described previously. The extracted RNA was dissolved in RNase-free water and quantified using a NanoDrop-2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). RNA purity was evaluated by measuring the absorbance ratio at 260:280 nm, while RNA integrity was confirmed by the detection of the 18S and 28S bands after electrophoresis in 1% agarose gels. The primer sequences are shown in Table S1. The standard curve was drawn for each gene, and the linear relationship and amplification efficiency were calculated. Primers were confirmed as follows: correlation coefficient (R2) > 0.98; the slope of the standard curve: −3–−3.5; and the amplification efficiency: 0.9–1.2. The mRNA levels of target genes were normalized to β-actin mRNA. Based on the cycle threshold (CT) value, the β-actin mRNA level was stable across the treatments in this study (p > 0.1) (Figure S8). The results of the relative mRNA expression of genes were calculated using the 2−ΔΔCT method [54].

4.6. 16S rRNA Gene Sequencing

Microbial genomic DNA was isolated from fecal specimens utilizing the E.Z.N.A.® soil DNA kit (Omega Bio-tek, Norcross, GA, USA). The V3-V4 hypervariable region of the bacterial 16S rRNA gene underwent amplification with the primer combinations 338F (5’-ACTCCTACGGGAGGCAGCAG-3’) and 806R (5’-GGACTACHVGGGTWTCTAAT-3’) via the T100 Thermal Cycler PCR thermocycler (BIO-RAD, Hercules, CA, USA). Subsequently, the purified amplicons were combined in equimolar proportions and subjected to paired-end sequencing on an Illumina PE300 platform (Illumina, San Diego, CA, USA) in accordance with standard operating procedures established by Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China).

4.7. Microbiota Data Analysis

Initial sequencing data underwent quality filtering using Fastp version 0.19.6, followed by read merging using FLASH version 1.2.11. Subsequently, the sequences were optimized and grouped into operational taxonomic units (OTUs) at a 97% sequence similarity threshold utilizing UPARSE 11. The taxonomic classification of each OTU representative sequence was conducted, employing RDP Classifier version 2.13 against the Silva v138 16S rRNA gene database, with a confidence threshold set at 0.7.
Analysis of fecal microbiota data was executed on the Majorbio Cloud platform (https://cloud.majorbio.com (accessed on 15 April 2024)). Rarefaction curves and alpha diversity indices were computed using Mothur v1.30.1. Principal coordinate analysis (PCoA) and Adonis analysis, based on Bray–Curtis dissimilarity, were employed to assess the similarity among microbial communities across different samples, facilitated by the “Vegan” package. The statistical significance of non-parametric relative abundance profiles was evaluated utilizing the Mann–Whitney test. Additionally, linear discriminant analysis (LDA) effect size (LEfSe) analysis was conducted to identify significantly enriched bacterial taxa among distinct groups (LDA score > 3, p < 0.05).

4.8. Liquid Chromatograph–Mass Spectrometer Analysis (LC–MS Analysis)

A solid sample weighing 50 mg underwent metabolite extraction utilizing 400 μL of a 4:1 (v:v) methanol–water extract solution. LC–MS/MS analysis was performed on a SCIEX UPLC-Triple TOF 6600 system featuring an ACQUITY HSS T3 column (100 mm × 2.1 mm i.d., 1.8 μm; Waters, Milford, MA, USA) at Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China). Mobile phase A is 95% water + 5% acetonitrile (containing 0.1% formic acid); mobile phase B is 47.5% acetonitrile + 47.5% isopropyl alcohol + 5% water (containing 0.1% formic acid). The flow rate was 0.40 mL/min, the sample size was 10 μL, and the column temperature was 45 °C. Positive ion mode separation gradient 0–3 min and mobile phase B increased from 0% to 20%; in 3–4.5 min, mobile phase B increased from 20% to 35%; in 4.5–5 min, mobile phase B increased from 35% to 100%; at 5–6.3 min, mobile phase B remained 100%; and in 6.3–6.4 min, mobile phase B decreased from 100% to 0%. At 6.4–8 min, mobile phase B maintained 0%. Negative ion mode separation gradient: 0–1.5 min, mobile phase B from 0 to 5%; 1.5–2 min, mobile phase B increased from 5% to 10%; 2–4.5 min, mobile phase B increased from 10% to 30%; 4.5–5 min, mobile phase B increased from 30% to 100%; and at 5–6.3 min, the linearity of mobile phase B remained 100%. In 6.3–6.4 min, mobile phase B decreases from 100% to 0%. The linearity of mobile phase B remained at 0% from 6.4–8 min. Positive and negative ion scanning modes are used to collect the sample quality spectrum signal, and the quality scanning range is 50–1200 m/z. The spray gas flow rate is 50 psi, the auxiliary gas flow rate is 50 psi, the curtain gas flow rate is 35 psi, the ion source temperature is 500 °C, the positive mode ion spray voltage is set to 5500 V, the negative mode ion spray voltage is set to −4500 V, and the ion transfer tube temperature is 325 °C. Normalized collision energy is 20-40-60 eV cyclic collision energy. DDA mode was used to collect data.
Raw UHPLC–MS data were processed using Progenesis QI 2.3 software (Nonlinear Dynamics, Waters, USA) for peak detection and alignment. Metabolic features detected in at least 80% of samples were retained, and minimum imputation was applied to samples falling below the quantitation limit. The normalization of each metabolic feature was conducted by sum and quality control. Variables with relative standard deviation (RSD) exceeding 30% of QC samples were excluded, and log10 transformation was performed to obtain the final data matrix.
Conduct variance analysis on the processed data matrix file. The “ropls” R package (Version 1.6.2) was utilized for partial least squares discriminant analysis (PLS-DA). Metabolites with a variable importance in projection (VIP) value exceeding 1 and a significance level of p < 0.05 were deemed significantly different. Differential metabolites between the two groups were subjected to metabolic enrichment and pathway analysis based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genome.jp/kegg/ (accessed on 15 April 2024)) to map their biochemical pathways.

4.9. RNA Extraction, Library Construction, and Sequencing

Total RNA was isolated from the tissue using TRIzol® Reagent as per the manufacturer’s protocol. Subsequently, RNA integrity was assessed using a 5300 Bioanalyser (Agilent, Santa Clara, CA, USA) and quantification was carried out using the ND-2000 (NanoDrop Technologies, Wilmington, DE, USA). RNA purification, reverse transcription, library preparation, and sequencing procedures were conducted at Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China) following the guidelines provided by Illumina (San Diego, CA, USA). The RNA-seq transcriptome library was constructed using Illumina® Stranded mRNA Prep and Ligation Kit from Illumina (San Diego, CA, USA), utilizing 1 μg of total RNA. Libraries underwent size selection for cDNA fragments of 300 bp on 2% low-range ultra-agarose, followed by PCR amplification with Phusion DNA polymerase (NEB) for 15 PCR cycles. Following quantification with Qubit 4.0, paired-end RNA-seq libraries were sequenced using the NovaSeq 6000 sequencer (Illumina, San Diego, CA, USA) (2 × 150 bp read length).
The raw paired-end reads underwent trimming and quality control using Fastp with default parameters. Subsequently, clean reads were individually aligned to the reference genome in orientation mode employing HISAT2 software 2.2.1. The mapped reads of each sample were assembled by StringTie v2.1.4 utilizing a reference-based approach.
To identify differential expression genes (DEGs) between distinct samples, the expression level of each transcript was quantified utilizing the transcripts per million reads (TPM) method via RSEM. Differential expression analysis was conducted using DESeq2 (version 1.12.3), with DEGs meeting the criteria of |log2FC| ≥ 1 and false discovery rate (FDR) < 0.05 considered significantly differentially expressed genes. Additionally, functional enrichment analysis, encompassing Gene Ontology (GO) and KEGG pathways, was performed to discern which DEGs were enriched in GO terms and metabolic pathways.

4.10. Correlation Analysis

Given the extensive diversity, dimensional complexity, and discontinuous distribution of omics data, we applied the Spearman correlation examination in this investigation to investigate the associations among the bacteria, metabolites, liver genes, gut genes, and host phenotypes. A threshold of |R| > 0.6, with p < 0.05, was utilized to discern noteworthy correlations. Cytoscape version 3.9.0 was employed for network construction. Visualization of the interrelations between microbiome, metabolome, and physiological attributes was accomplished utilizing the R package “ggplot2”.

4.11. Statistical Analysis

The body weight, ELISA, and histological observation data were analyzed and one-way analysis of variance (ANOVA) by Statistical Package for the Social Sciences (SPSS) version 26; the differences between means were compared by Duncan’s multiple range test. Gene relative expression was compared using the t-test. The results were expressed as mean ± standard deviation (SD). The significance was set at p < 0.05.

Supplementary Materials

The supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms251910623/s1.

Author Contributions

M.Z.: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Visualization, Writing—original draft, and Writing—review and editing. Y.C.: Data curation, Formal analysis, Methodology, and Writing—original draft. F.W. (Fengxia Wang): Data curation, Formal analysis, Methodology, and Writing—original draft. F.W. (Fengyang Wu): Data curation, Validation, and Methodology. C.L.: Data curation, Validation, and Methodology. S.L.: Data curation, Validation, and Methodology. B.C.: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, and Writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Modern Agriculture Industry Technology System of Rabbit, China (CARS-43-B-2) and the Postgraduate Innovation Funding project of Hebei Province (CXZZBS2024075).

Institutional Review Board Statement

All experiments were performed in accordance with the ARRIVE guidelines and were carried out in accordance with the U.K. Animals (Scientific Procedures) Act, 1986, associated guidelines, EU Directive 2010/63/EU for animal experiments and approved by the Ethics Committee of Animal Experimentation of Hebei Agricultural University (Protocol 2023058).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in the NCBI Sequence Read Archive (SRA), reference number PRJNA1115940.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ma, Z.; Akhtar, M.; Pan, H.; Liu, Q.; Chen, Y.; Zhou, X.; You, Y.; Shi, D.; Liu, H. Fecal microbiota transplantation improves chicken growth performance by balancing jejunal Th17/Treg cells. Microbiome 2023, 11, 137. [Google Scholar] [CrossRef] [PubMed]
  2. Cargnin, S.T.; Gnoatto, S.B. Ursolic acid from apple pomace and traditional plants: A valuable triterpenoid with functional properties. Food Chem. 2017, 220, 477–489. [Google Scholar] [CrossRef] [PubMed]
  3. Fu, Y.; Liu, T.; He, S.; Zhang, Y.; Tan, Y.; Bai, Y.; Shi, J.; Deng, W.; Qiu, J.; Wang, Z.; et al. Ursolic acid reduces oxidative stress injury to ameliorate experimental autoimmune myocarditis by activating Nrf2/HO-1 signaling pathway. Front. Pharmacol. 2023, 14, 1189372. [Google Scholar] [CrossRef] [PubMed]
  4. Gou, W.; Luo, N.; Yu, B.; Wu, H.; Wu, S.; Tian, C.; Guo, J.; Ning, H.; Bi, C.; Wei, H.; et al. Ursolic Acid Derivative UA232 Promotes Tumor Cell Apoptosis by Inducing Endoplasmic Reticulum Stress and Lysosomal Dysfunction. Int. J. Biol. Sci. 2022, 18, 2639–2651. [Google Scholar] [CrossRef] [PubMed]
  5. Wan, Y.; Zhang, W.; Huang, C.; Jian, J.; Zhang, Y.; Liu, Q.; Chen, P.; Zhu, X. Ursolic acid alleviates Kupffer cells pyroptosis in liver fibrosis by the NOX2/NLRP3 inflammasome signaling pathway. Int. Immunopharmacol. 2022, 113, 109321. [Google Scholar] [CrossRef]
  6. Zhao, M.; Wu, F.; Tang, Z.; Yang, X.; Liu, Y.; Wang, F.; Chen, B. Anti-inflammatory and antioxidant activity of ursolic acid: A systematic review and meta-analysis. Front. Pharmacol. 2023, 14, 1256946. [Google Scholar] [CrossRef]
  7. Zhang, B.C.; Lai, C.M.; Luo, B.Y.; Shao, J.W. Triterpenoids-templated self-assembly nanosystem for biomimetic delivery of CRISPR/Cas9 based on the synergy of TLR-2 and ICB to enhance HCC immunotherapy. Acta Pharm. Sin. B 2024, 14, 3205–3217. [Google Scholar] [CrossRef]
  8. Larabi, A.; Barnich, N.; Nguyen, H.T.T. New insights into the interplay between autophagy, gut microbiota and inflammatory responses in IBD. Autophagy 2020, 16, 38–51. [Google Scholar] [CrossRef]
  9. Singh, V.; Yeoh, B.S.; Walker, R.E.; Xiao, X.; Saha, P.; Golonka, R.M.; Cai, J.; Bretin, A.C.A.; Cheng, X.; Liu, Q.; et al. Microbiota fermentation-NLRP3 axis shapes the impact of dietary fibres on intestinal inflammation. Gut 2019, 68, 1801–1812. [Google Scholar] [CrossRef]
  10. Cheng, H.; Guan, X.; Chen, D.; Ma, W. The Th17/Treg Cell Balance: A Gut Microbiota-Modulated Story. Microorganisms 2019, 7, 583. [Google Scholar] [CrossRef]
  11. Zhang, L.L.; Xu, J.Y.; Xing, Y.; Wu, P.; Jin, Y.W.; Wei, W.; Zhao, L.; Yang, J.; Chen, G.C.; Qin, L.Q. Lactobacillus rhamnosus GG alleviates radiation-induced intestinal injury by modulating intestinal immunity and remodeling gut microbiota. Microbiol. Res. 2024, 286, 127821. [Google Scholar] [CrossRef] [PubMed]
  12. Yan, J.-B.; Luo, M.-M.; Chen, Z.-Y.; He, B.-H. The Function and Role of the Th17/Treg Cell Balance in Inflammatory Bowel Disease. J. Immunol. Res. 2020, 2020, 8813558. [Google Scholar] [CrossRef] [PubMed]
  13. Yan, X.; Ren, X.; Liu, X.; Wang, Y.; Ma, J.; Song, R.; Wang, X.; Dong, Y.; Fan, Q.; Wei, J.; et al. Dietary Ursolic Acid Prevents Alcohol-Induced Liver Injury via Gut–Liver Axis Homeostasis Modulation: The Key Role of Microbiome Manipulation. J. Agric. Food Chem. 2021, 69, 7074–7083. [Google Scholar] [CrossRef]
  14. Sheng, Q.; Li, F.; Chen, G.; Li, J.; Li, J.; Wang, Y.; Lu, Y.; Li, Q.; Li, M.; Chai, K. Ursolic Acid Regulates Intestinal Microbiota and Inflammatory Cell Infiltration to Prevent Ulcerative Colitis. J. Immunol. Res. 2021, 2021, 6679316. [Google Scholar] [CrossRef]
  15. Nie, Y.; Liu, Q.; Zhang, W.; Wan, Y.; Huang, C.; Zhu, X. Ursolic acid reverses liver fibrosis by inhibiting NOX4/NLRP3 inflammasome pathways and bacterial dysbiosis. Gut Microbes 2021, 13, 1972746. [Google Scholar] [CrossRef]
  16. Kashyap, D.; Tuli, H.S.; Sharma, A.K. Ursolic acid (UA): A metabolite with promising therapeutic potential. Life Sci. 2016, 146, 201–213. [Google Scholar] [CrossRef]
  17. Chen, C.; Gao, K.; Chen, Z.; Zhang, Q.; Ke, X.; Mao, B.; Fan, Q.; Li, Y.; Chen, S. The supplementation of the multi-strain probiotics WHHPRO™ alleviates high-fat diet-induced metabolic symptoms in rats via gut-liver axis. Front. Nutr. 2023, 10, 1324691. [Google Scholar] [CrossRef]
  18. Li, J.; Li, N.; Yan, S.; Liu, M.; Sun, B.; Lu, Y.; Shao, Y. Ursolic acid alleviates inflammation and against diabetes-induced nephropathy through TLR4-mediated inflammatory pathway. Mol. Med. Rep. 2018, 18, 4675–4681. [Google Scholar] [CrossRef] [PubMed]
  19. Li, H.; Yu, Y.; Liu, Y.; Luo, Z.; Law, B.Y.K.; Zheng, Y.; Huang, X.; Li, W. Ursolic acid enhances the antitumor effects of sorafenib associated with Mcl-1-related apoptosis and SLC7A11-dependent ferroptosis in human cancer. Pharmacol. Res. 2022, 182, 106306. [Google Scholar] [CrossRef]
  20. Knochelmann, H.M.; Dwyer, C.J.; Bailey, S.R.; Amaya, S.M.; Elston, D.M.; Mazza-McCrann, J.M.; Paulos, C.M. When worlds collide: Th17 and Treg cells in cancer and autoimmunity. Cell. Mol. Immunol. 2018, 15, 458–469. [Google Scholar] [CrossRef]
  21. Liu, Y.J.; Tang, B.; Wang, F.C.; Tang, L.; Lei, Y.Y.; Luo, Y.; Huang, S.J.; Yang, M.; Wu, L.Y.; Wang, W.; et al. Parthenolide ameliorates colon inflammation through regulating Treg/Th17 balance in a gut microbiota-dependent manner. Theranostics 2020, 10, 5225–5241. [Google Scholar] [CrossRef] [PubMed]
  22. Zhong, Y.; Liu, W.; Xiong, Y.; Li, Y.; Wan, Q.; Zhou, W.; Zhao, H.; Xiao, Q.; Liu, D. Astragaloside IV alleviates ulcerative colitis by regulating the balance of Th17/Treg cells. Phytomedicine 2022, 104, 154287. [Google Scholar] [CrossRef] [PubMed]
  23. Chung, Y.; Chang, S.H.; Martinez, G.J.; Yang, X.O.; Nurieva, R.; Kang, H.S.; Ma, L.; Watowich, S.S.; Jetten, A.M.; Tian, Q.; et al. Critical regulation of early Th17 cell differentiation by interleukin-1 signaling. Immunity 2009, 30, 576–587. [Google Scholar] [CrossRef] [PubMed]
  24. Noack, M.; Miossec, P. Th17 and regulatory T cell balance in autoimmune and inflammatory diseases. Autoimmun. Rev. 2014, 13, 668–677. [Google Scholar] [CrossRef]
  25. Wang, J.; Zhao, X.; Wan, Y.Y. Intricacies of TGF-β signaling in Treg and Th17 cell biology. Cell. Mol. Immunol. 2023, 20, 1002–1022. [Google Scholar] [CrossRef]
  26. Kamada, N.; Seo, S.-U.; Chen, G.Y.; Núñez, G. Role of the gut microbiota in immunity and inflammatory disease. Nat. Rev. Immunol. 2013, 13, 321–335. [Google Scholar] [CrossRef]
  27. Wu, Z.; Huang, S.; Li, T.; Li, N.; Han, D.; Zhang, B.; Xu, Z.Z.; Zhang, S.; Pang, J.; Wang, S.; et al. Gut microbiota from green tea polyphenol-dosed mice improves intestinal epithelial homeostasis and ameliorates experimental colitis. Microbiome 2021, 9, 184. [Google Scholar] [CrossRef]
  28. Hao, W.; Kwek, E.; He, Z.; Zhu, H.; Liu, J.; Zhao, Y.; Ma, K.Y.; He, W.S.; Chen, Z.Y. Ursolic acid alleviates hypercholesterolemia and modulates the gut microbiota in hamsters. Food Funct. 2020, 11, 6091–6103. [Google Scholar] [CrossRef]
  29. Zhou, J.; Li, M.; Chen, Q.; Li, X.; Chen, L.; Dong, Z.; Zhu, W.; Yang, Y.; Liu, Z.; Chen, Q. Programmable probiotics modulate inflammation and gut microbiota for inflammatory bowel disease treatment after effective oral delivery. Nat. Commun. 2022, 13, 3432. [Google Scholar] [CrossRef]
  30. Lima, S.F.; Gogokhia, L.; Viladomiu, M.; Chou, L.; Putzel, G.; Jin, W.B.; Pires, S.; Guo, C.J.; Gerardin, Y.; Crawford, C.V.; et al. Transferable Immunoglobulin A–Coated Odoribacter splanchnicus in Responders to Fecal Microbiota Transplantation for Ulcerative Colitis Limits Colonic Inflammation. Gastroenterology 2021, 162, 166–178. [Google Scholar] [CrossRef]
  31. Xing, C.; Wang, M.; Ajibade, A.A.; Tan, P.; Fu, C.; Chen, L.; Zhu, M.; Hao, Z.Z.; Chu, J.; Yu, X.; et al. Microbiota regulate innate immune signaling and protective immunity against cancer. Cell Host Microbe 2021, 29, 959–974.E7. [Google Scholar] [CrossRef] [PubMed]
  32. Kverka, M.; Zakostelska, Z.; Klimesova, K.; Sokol, D.; Hudcovic, T.; Hrncir, T.; Rossmann, P.; Mrazek, J.; Kopecny, J.; Verdu, E.F.; et al. Oral administration of Parabacteroides distasonis antigens attenuates experimental murine colitis through modulation of immunity and microbiota composition. Clin. Exp. Immunol. 2011, 163, 250–259. [Google Scholar] [CrossRef]
  33. Koh, G.Y.; Kane, A.V.; Wu, X.; Crott, J.W. Parabacteroides distasonis attenuates tumorigenesis, modulates inflammatory markers and promotes intestinal barrier integrity in azoxymethane-treated A/J mice. Carcinogenesis 2020, 41, 909–917. [Google Scholar] [CrossRef]
  34. Sun, H.; Guo, Y.; Wang, H.; Yin, A.; Hu, J.; Yuan, T.; Zhou, S.; Xu, W.; Wei, P.; Yin, S.; et al. Gut commensal Parabacteroides distasonis alleviates inflammatory arthritis. Gut 2023, 72, 1664–1677. [Google Scholar] [CrossRef]
  35. Cuffaro, B.; Assohoun, A.L.W.; Boutillier, D.; Súkeníková, L.; Desramaut, J.; Boudebbouze, S.; Salomé-Desnoulez, S.; Hrdý, J.; Waligora-Dupriet, A.J.; Maguin, E.; et al. In Vitro Characterization of Gut Microbiota-Derived Commensal Strains: Selection of Parabacteroides distasonis Strains Alleviating TNBS-Induced Colitis in Mice. Cells 2020, 9, 2104. [Google Scholar] [CrossRef] [PubMed]
  36. Ezeji, J.C.; Sarikonda, D.K.; Hopperton, A.; Erkkila, H.L.; Cohen, D.E.; Martinez, S.P.; Cominelli, F.; Kuwahara, T.; Dichosa, A.E.K.; Good, C.E.; et al. Parabacteroides distasonis: Intriguing aerotolerant gut anaerobe with emerging antimicrobial resistance and pathogenic and probiotic roles in human health. Gut Microbes 2021, 13, 1922241. [Google Scholar] [CrossRef]
  37. Qi, H.; Li, Y.; Yun, H.; Zhang, T.; Huang, Y.; Zhou, J.; Yan, H.; Wei, J.; Liu, Y.; Zhang, Z.; et al. Lactobacillus maintains healthy gut mucosa by producing L-Ornithine. Commun. Biol. 2019, 2, 171. [Google Scholar] [CrossRef] [PubMed]
  38. Han, C.; Wu, X.; Zou, N.; Zhang, Y.; Yuan, J.; Gao, Y.; Chen, W.; Yao, J.; Li, C.; Hou, J.; et al. Cichorium pumilum Jacq Extract Inhibits LPS-Induced Inflammation via MAPK Signaling Pathway and Protects Rats from Hepatic Fibrosis Caused by Abnormalities in the Gut-Liver Axis. Front. Pharmacol. 2021, 12, 683613. [Google Scholar] [CrossRef]
  39. Krishna, M.; Narang, H. The complexity of mitogen-activated protein kinases (MAPKs) made simple. Cell. Mol. Life Sci. 2008, 65, 3525–3544. [Google Scholar] [CrossRef]
  40. Kyriakis, J.M.; Avruch, J. Mammalian MAPK signal transduction pathways activated by stress and inflammation: A 10-year update. Physiol. Rev. 2012, 92, 689–737. [Google Scholar] [CrossRef]
  41. Petrasek, J.; Bala, S.; Csak, T.; Lippai, D.; Kodys, K.; Menashy, V.; Barrieau, M.; Min, S.Y.; Kurt-Jones, E.A.; Szabo, G. IL-1 receptor antagonist ameliorates inflammasome-dependent alcoholic steatohepatitis in mice. J. Clin. Investig. 2012, 122, 3476–3489. [Google Scholar] [CrossRef] [PubMed]
  42. Wang, C.; Li, X.; Xue, B.; Yu, C.; Wang, L.; Deng, R.; Liu, H.; Chen, Z.; Zhang, Y.; Fan, S.; et al. RasGRP1 promotes the acute inflammatory response and restricts inflammation-associated cancer cell growth. Nat. Commun. 2022, 13, 7001. [Google Scholar] [CrossRef] [PubMed]
  43. Pandian, V.D.; Giovannucci, D.R.; Vazquez, G.; Kumarasamy, S. CACNB2 is associated with aberrant RAS-MAPK signaling in hypertensive Dahl Salt-Sensitive rats. Biochem. Biophys. Res. Commun. 2019, 513, 760–765. [Google Scholar] [CrossRef]
  44. Mora, J.R.; Iwata, M.; von Andrian, U.H. Vitamin effects on the immune system: Vitamins A and D take centre stage. Nat. Rev. Immunol. 2008, 8, 685–698. [Google Scholar] [CrossRef]
  45. Carr, A.C.; Maggini, S. Vitamin C and Immune Function. Nutrients 2017, 9, 1211. [Google Scholar] [CrossRef]
  46. Wilson, D.F. Oxidative phosphorylation: Regulation and role in cellular and tissue metabolism. J. Physiol. 2017, 595, 7023–7038. [Google Scholar] [CrossRef] [PubMed]
  47. Wang, L.; Zheng, W.; Sun, Y.; Ren, X.; Yan, C.; Song, S.; Ai, C. Fucoidan ameliorates alcohol-induced liver injury in mice through Parabacteroides distasonis-mediated regulation of the gut-liver axis. Int. J. Biol. Macromol. 2024, 279, 135309. [Google Scholar] [CrossRef]
  48. Zhao, Q.; Dai, M.Y.; Huang, R.Y.; Duan, J.Y.; Zhang, T.; Bao, W.M.; Zhang, J.Y.; Gui, S.Q.; Xia, S.M.; Dai, C.T.; et al. Parabacteroides distasonis ameliorates hepatic fibrosis potentially via modulating intestinal bile acid metabolism and hepatocyte pyroptosis in male mice. Nat. Commun. 2023, 14, 1829. [Google Scholar] [CrossRef]
  49. Shu, Y.; Huang, Y.; Dong, W.; Fan, X.; Sun, Y.; Chen, G.; Zeng, X.; Ye, H. The polysaccharides from Auricularia auricula alleviate non-alcoholic fatty liver disease via modulating gut microbiota and bile acids metabolism. Int. J. Biol. Macromol. 2023, 246, 125662. [Google Scholar] [CrossRef]
  50. Canbay, A.; Sowa, J.P. L-Ornithine L-Aspartate (LOLA) as a Novel Approach for Therapy of Non-alcoholic Fatty Liver Disease. Drugs 2019, 79 (Suppl. 1), 39–44. [Google Scholar] [CrossRef]
  51. Ding, Y.; Koda, Y.; Shashni, B.; Takeda, N.; Zhang, X.; Tanaka, N.; Nishikawa, Y.; Nagasaki, Y. An orally deliverable ornithine-based self-assembling polymer nanomedicine ameliorates hyperammonemia in acetaminophen-induced acute liver injury. Acta Biomater. 2023, 168, 515–528. [Google Scholar] [CrossRef] [PubMed]
  52. Imam, K.M.S.U.; Tian, Y.; Xin, F.; Xie, Y.; Wen, B. Lactucin, a Bitter Sesquiterpene from Cichorium intybus, Inhibits Cancer Cell Proliferation by Downregulating the MAPK and Central Carbon Metabolism Pathway. Molecules 2022, 27, 7358. [Google Scholar] [CrossRef] [PubMed]
  53. Dang, T.; Zheng, G.; Zhang, Q.; Jin, P.; Zhang, H.; Su, L.; Qin, D.; Yao, G. Sesquiterpenoids with diverse carbon skeletons from the roots of Cichorium glandulosum and their anti-inflammatory activities. Fitoterapia 2019, 136, 104170. [Google Scholar] [CrossRef] [PubMed]
  54. Zhou, W.; Zhang, H.; Huang, L.; Sun, C.; Yue, Y.; Cao, X.; Jia, H.; Wang, C.; Gao, Y. Disulfiram with Cu2+ alleviates dextran sulfate sodium-induced ulcerative colitis in mice. Theranostics 2023, 13, 2879–2895. [Google Scholar] [CrossRef]
Figure 1. Effects of ursolic acid on body weight and immune-related cytokines in mice. (A) Body weight (*, p < 0.05, and ***, p < 0.001). (B) IL-1β, (C) IL-6, (D) TNF-α, (E) TGF-β, (F) IL-10, and (G) IL-17A cytokine levels in serum. Values with different letters (a, b, and c) differed significantly within a mice trial (p < 0.05). Data are shown as mean ± SD. n = 6.
Figure 1. Effects of ursolic acid on body weight and immune-related cytokines in mice. (A) Body weight (*, p < 0.05, and ***, p < 0.001). (B) IL-1β, (C) IL-6, (D) TNF-α, (E) TGF-β, (F) IL-10, and (G) IL-17A cytokine levels in serum. Values with different letters (a, b, and c) differed significantly within a mice trial (p < 0.05). Data are shown as mean ± SD. n = 6.
Ijms 25 10623 g001
Figure 2. Effect of ursolic acid on the intestinal health of mice. Effect of ursolic acid on the intestinal morphology of the duodenum (A,B), jejunum (C,D), and ileum (E,F) in mice. Effects of ursolic acid on the expression of genes related to the intestinal barrier (GI) and Th17/Treg cell-related factor (JL) in the jejunum of mice. Values with different letters (a, b, c, and d) differed significantly within a mice trial (p < 0.05). Data are shown as mean ± SD. n = 6. Scale bars in (A,C,E) = 100 μm.
Figure 2. Effect of ursolic acid on the intestinal health of mice. Effect of ursolic acid on the intestinal morphology of the duodenum (A,B), jejunum (C,D), and ileum (E,F) in mice. Effects of ursolic acid on the expression of genes related to the intestinal barrier (GI) and Th17/Treg cell-related factor (JL) in the jejunum of mice. Values with different letters (a, b, c, and d) differed significantly within a mice trial (p < 0.05). Data are shown as mean ± SD. n = 6. Scale bars in (A,C,E) = 100 μm.
Ijms 25 10623 g002
Figure 3. Effect of ursolic acid on gut microbiota in mice. (A) The comparison of microbiota diversity (Shannon index) and the Adonis analysis. Boxes denote the interquartile (IQR) between the first and third quartiles (25th and 75th percentiles, respectively) and the line inside denotes the median. (B) Principal coordinate analysis (PCoA) based on Bray–Curtis distance. (C) Comparison of differences in the relative abundance of bacteria at the genus level (*, p < 0.05, and **, p < 0.01). (D) Differentiating bacteria identified by LEfSe with LDA effect size ≥ 3. n = 6.
Figure 3. Effect of ursolic acid on gut microbiota in mice. (A) The comparison of microbiota diversity (Shannon index) and the Adonis analysis. Boxes denote the interquartile (IQR) between the first and third quartiles (25th and 75th percentiles, respectively) and the line inside denotes the median. (B) Principal coordinate analysis (PCoA) based on Bray–Curtis distance. (C) Comparison of differences in the relative abundance of bacteria at the genus level (*, p < 0.05, and **, p < 0.01). (D) Differentiating bacteria identified by LEfSe with LDA effect size ≥ 3. n = 6.
Ijms 25 10623 g003
Figure 4. Effect of ursolic acid on the fecal metabolism of mice. (A) Heatmap of 16 different metabolites. The right side shows the variable importance in the projection (VIP) bar graph of the metabolite. The color of the bars indicates the significant difference in metabolites between the two groups of samples, and the smaller the p-value, the darker the color (*, p < 0.05; **, p < 0.01). (B) KEGG pathway enrichment analysis of differential metabolites. (C) Interrelationship between microbiota, metabolomics, and physiological features. Red connections indicate a positive correlation (Spearman correlation test, |R| > 0.6, p < 0.05,), while blue connections show correlations that were negative (Spearman correlation test, |R| > 0.6, p < 0.05). In the microbiota column, the orange stratum represents bacteria that were highly enriched in the UA group, and the stratum-colored blue was increased in the control group. In the metabolomics column, the purple stratum represents metabolites that were highly enriched in the UA group, and the stratum-colored ice blue was increased in the control group. In the physiological feature column, the green stratum represents bacteria that were highly enriched in the UA group, and the pink stratum was increased in the control group. n = 6.
Figure 4. Effect of ursolic acid on the fecal metabolism of mice. (A) Heatmap of 16 different metabolites. The right side shows the variable importance in the projection (VIP) bar graph of the metabolite. The color of the bars indicates the significant difference in metabolites between the two groups of samples, and the smaller the p-value, the darker the color (*, p < 0.05; **, p < 0.01). (B) KEGG pathway enrichment analysis of differential metabolites. (C) Interrelationship between microbiota, metabolomics, and physiological features. Red connections indicate a positive correlation (Spearman correlation test, |R| > 0.6, p < 0.05,), while blue connections show correlations that were negative (Spearman correlation test, |R| > 0.6, p < 0.05). In the microbiota column, the orange stratum represents bacteria that were highly enriched in the UA group, and the stratum-colored blue was increased in the control group. In the metabolomics column, the purple stratum represents metabolites that were highly enriched in the UA group, and the stratum-colored ice blue was increased in the control group. In the physiological feature column, the green stratum represents bacteria that were highly enriched in the UA group, and the pink stratum was increased in the control group. n = 6.
Ijms 25 10623 g004
Figure 5. Effect of ursolic acid on the liver transcriptome of mice. (A) PCA analysis. (B) The number of differentially expressed genes (DEGs). KEGG enrichment pathways of down-regulated (C) and up-regulated (D) DEGs. (E) DEGs related to enrichment pathways (ko04010, ko04668, ko00053, and ko00190). ko04010, MAPK signaling pathway; ko04668, TNF signaling pathway; ko00053, Ascorbate and aldarate metabolism; and ko00190, Oxidative phosphorylation; I, Nicotinamide adenine dinucleotide (NADH) dehydrogenase; II, Succinate dehydrogenase; III, Cytochrome C Reductase; IV, Cytochrome c oxidase; V, ATP synthase. n = 6.
Figure 5. Effect of ursolic acid on the liver transcriptome of mice. (A) PCA analysis. (B) The number of differentially expressed genes (DEGs). KEGG enrichment pathways of down-regulated (C) and up-regulated (D) DEGs. (E) DEGs related to enrichment pathways (ko04010, ko04668, ko00053, and ko00190). ko04010, MAPK signaling pathway; ko04668, TNF signaling pathway; ko00053, Ascorbate and aldarate metabolism; and ko00190, Oxidative phosphorylation; I, Nicotinamide adenine dinucleotide (NADH) dehydrogenase; II, Succinate dehydrogenase; III, Cytochrome C Reductase; IV, Cytochrome c oxidase; V, ATP synthase. n = 6.
Ijms 25 10623 g005
Figure 6. Spearman correlation network among the bacteria, metabolites, liver genes, gut genes, and host phenotypes (|R| > 0.6, p < 0.05). The width of the edges is proportional to the correlation strength. The color of the edges: red, positive; blue, and negative, n = 6.
Figure 6. Spearman correlation network among the bacteria, metabolites, liver genes, gut genes, and host phenotypes (|R| > 0.6, p < 0.05). The width of the edges is proportional to the correlation strength. The color of the edges: red, positive; blue, and negative, n = 6.
Ijms 25 10623 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, M.; Cui, Y.; Wang, F.; Wu, F.; Li, C.; Liu, S.; Chen, B. Ursolic Acid Regulates Immune Balance, Modulates Gut Microbial Metabolism, and Improves Liver Health in Mice. Int. J. Mol. Sci. 2024, 25, 10623. https://doi.org/10.3390/ijms251910623

AMA Style

Zhao M, Cui Y, Wang F, Wu F, Li C, Liu S, Chen B. Ursolic Acid Regulates Immune Balance, Modulates Gut Microbial Metabolism, and Improves Liver Health in Mice. International Journal of Molecular Sciences. 2024; 25(19):10623. https://doi.org/10.3390/ijms251910623

Chicago/Turabian Style

Zhao, Man, Yali Cui, Fengxia Wang, Fengyang Wu, Chong Li, Shudong Liu, and Baojiang Chen. 2024. "Ursolic Acid Regulates Immune Balance, Modulates Gut Microbial Metabolism, and Improves Liver Health in Mice" International Journal of Molecular Sciences 25, no. 19: 10623. https://doi.org/10.3390/ijms251910623

APA Style

Zhao, M., Cui, Y., Wang, F., Wu, F., Li, C., Liu, S., & Chen, B. (2024). Ursolic Acid Regulates Immune Balance, Modulates Gut Microbial Metabolism, and Improves Liver Health in Mice. International Journal of Molecular Sciences, 25(19), 10623. https://doi.org/10.3390/ijms251910623

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop