- State Key Laboratory of Food Science and Resources, School of Food Science, Nanchang University, Nanchang, China
Dysbiosis of gut microbiota and metabolic pathway disorders are closely related to the ulcerative colitis. Through network pharmacology, we found that puerarin is a potential ingredient that can improve the crypt deformation and inflammatory infiltration in mice, and decrease the levels of IL-1β, IL-6 and TNF-α significantly. Listeria, Alistipes and P. copri gradually became dominant bacteria in UC mice, which were positively correlated with inflammatory factors. Puerarin effectively improved dysbiosis by reducing the abundance of Alistipes, P. copri and Veillonella, and increasing the level of Desulfovibrionacea. Correlation network and metabolic function prediction analysis of the microbiota showed that they formed a tightly connected network and were widely involved in carbohydrate metabolism and amino acid metabolism. Specifically, we observed significant changes in the tryptophan metabolism pathway in DSS mice, with an increase in the abundance of Bacteroidetes and Enterobacteriaceae involved in tryptophan metabolism. However, this metabolic disorder was alleviated after puerarin treatment, including the reversal of 3-HAA levels and an increase in the abundance of Rhodobacteraceae and Halomonadaceae involved in kynurenine metabolism, as well as a significant increase in the purine metabolite guanosine. In conclusion, our study suggests that puerarin has a good therapeutic effect on UC, which is partially achieved by restoring the composition and abundance of gut microbiota and their metabolism.
1. Introduction
Ulcerative colitis (UC) is a relapsing disease characterized by chronic inflammation. Clinical symptoms include diarrhea, abdominal pain, and hematochezia. According to the latest data, the prevalence of UC was estimated to be 5 million cases around the world, and there is an increasing trend in some emerging industrialized countries (Le Berre et al., 2023; Liu J. et al., 2023).
The pathogenesis of UC is multifactorial, of which the epithelial barrier defects and disorder of immune response are important factors. Persistent inflammation in the gut can lead to dysbiosis of gut microbiota and abnormal accumulation of microbiota metabolites. These abnormal metabolites can further stimulate the immune system, ultimately resulting in impaired mucosal barrier function and exacerbated inflammation (Neurath, 2014).
The current mainstream treatment for UC is the use of salicylic acid drugs to inhibit the synthesis of prostaglandins and leukotrienes. However, they have significant side effects, such as nausea, rash and liver dysfunction, which make them unsuitable for long-term use (Kucharzik et al., 2020). Given the long course and high recurrence rate of UC, it is necessary to search other natural active substances for intervention to achieve a relieving effect.
Radix Puerariae is the dried root of the leguminous plant Pueraria lobata, which contains abundant isoflavones and has been widely used in the food and health product industries (Wang et al., 2022). Previous studies have found that water extracts containing Radix Puerariae can increase tight junction protein expression (Wang et al., 2023), lower NF-κB activation, reduce inflammatory factor expression in the colon and relieve oxidative stress (Li et al., 2016). These findings suggest that Radix Puerariae has a beneficial regulatory effect on intestinal barrier function and inflammation. However, its specific active ingredients and their effects on the composition of the intestinal microbiome and metabolite levels in UC have not been fully elucidated.
In recent years, with the proposal of disease network and systems pharmacology theory (Nogales et al., 2022), it has become a convenient and efficient method to systematically analyze and predict the “substance-gene-target-disease” action network and find potential effective active ingredients based on technologies such as genomics, high-throughput screening, network visualization, and network analysis (Wang et al., 2021).
Based on this, this article first uses network pharmacology to mine the single active ingredients in Radix Puerariae, and then establishes a DSS-induced colitis mouse model to verify its improvement effect. In order to further explore possible mechanisms of action, untargeted metabolomics and 16S rRNA technology are used to investigate the improvement of gut microbiota and metabolites by active ingredients.
2. Materials and methods
2.1. Virtual compound screening
We utilized the TCMSP database1 to retrieve compounds by searching for ‘Radix Puerariae’ as the primary keywords. To screen the active ingredients and corresponding targets, we applied OB ≥ 30% and DL ≥ 0.18 as screening conditions (Jiao et al., 2022). Additionally, we included substances with active ingredients that have been widely reported in Radix Puerariae for subsequent analysis. Then combined these active targets with compound active targets obtained from the Drugbank database2 to obtain the total active targets for pueraria. Subsequently, we used uniprot database3 to convert the corresponding gene names (Zhang H. et al., 2021).
For searching the disease targets, we used DisGeNET4 and GeneCards database.5 Then used Cytoscape software to map the cross-gene information between UC and Radix Puerariae. Additionally, we conducted PPI analysis using the STRING database,6 as well as GO and KEGG analyses using the DAVID database7 (Li et al., 2022).
2.2. Molecular docking
We searched the PubChem database8 to obtain the 3D structure files of puerarin, genistein, formononetin, and daidzein. And downloaded the crystal structures of TNF-α, IL-1B, and IL-10 from the PDB database.9 We then docked the small molecules with these proteins and visualized the results of small molecule-protein binding pockets and interactions between the small molecules and amino acid residues by using PyMol 2.5.5 software.
2.3. Animals and experimental treatment
C57BL/6J mice (male, 6–7 weeks old, 20–23 g) were obtained from Spafu Biotechnology (Beijing, China) and housed in a SPF animal facility with a 12 h light/dark cycle. The mice were provided with standard rodent feed and purified water, and the ambient relative humidity was maintained between 40 and 60% at a temperature of 25°C ± 2°C. The experiment was approved by the Animal Experimentation Ethical Committee of Nanchang University (License number: SYXK 2021-0001), and all animal feeding and procedures were conducted in accordance with the experimental animal welfare ethical code of Nanchang University.
The animal experiment procedures are depicted in Figure 1C. After 1 week of adaptive feeding, mice were randomly assigned to four groups: Control group (CON), Model group (UC), Positive group (Sulfasalazine, SASP), and Puerarin group (PUE). All mice were fed the same diet. During the first week, the Model group, Positive group, and Puerarin group received intragastric administration of DSS aqueous solution three times a day, while the Control group received an equivalent amount of normal saline. For the next 2 weeks, the Puerarin group received 200 mg/kg·BW puerarin via gavage, while the Positive group received 200 mg/kg·BW SASP via gavage (Huang et al., 2017). The Model and Control groups were administered an equivalent volume of normal saline via gavage. The weight and diet of all mice were monitored on a weekly basis.
Figure 1. Active ingredient puerarin was screened out by network pharmacology, and it has a good effect on the growth performance of mice. (A) Diagram of overlapping genes’ network between UC and active components in Radix Puerariae. (B) Protein–protein interaction analysis of overlapping genes. (C) Schematic diagram of mouse modeling. (D) Changes in body weight of mice during the experiment. (E) Changes of food intake of mice compared with CON group during the experiment. **p < 0.01, compared to the CON group.
2.4. Colitis induction and evaluation
To ensure that each mouse receives an equal dose of DSS during modeling, we induced colitis by gavage with DSS in this study. During the first week, mice in the model group were gavaged with 0.1 mL/g (0.23 mL) of 15% DSS water solution at 10:00, 15:00 and 20:00 each day until the end of modeling, with an average cumulative intake of DSS of 31.5 mg/g per mouse. After modeling, the DAI index was used for preliminary evaluation of the model (Peng et al., 2020). The detailed scoring criteria can be found in Supplementary Table S1.
2.5. Biochemical assays of the serum
Whole blood was obtained by puncturing the submandibular vein. To be specific, we inserted a blood collection needle into the posterior aspect of the top of the mandible, where the vein was located, and collected blood. After that, we applied a dry cotton ball to the puncture site to stop bleeding. The whole blood of mice was placed in the tube and centrifuged at 3000 rpm for 15 min at 4°C to obtain serum. The concentrations of TNF-α, IL-6, and IL-1β in the serum were measured using an ELISA kit in accordance with the manufacturer’s instructions (Zhao et al., 2021).
2.6. Hematoxylin and eosin staining
The colon tissue was fixed in 4% paraformaldehyde, embedded in paraffin to make 4 mm thick sections, stained with hematoxylin and eosin (H&E) staining solution, and observed under the orthostatic optical microscope (Zhang et al., 2022).
2.7. 16S rRNA and microbiota data analysis
We used QIAamp DNA Stool Kit (Qiagen, Hilden, Germany) and extracted the DNA according to the instructions. The V3-V4 region of the gene was amplified using primers with barcodes, and the resulting amplicons were sequenced on the Illumina NovaSeq platform (Yan et al., 2020).
The initial sequencing data were processed using QIIME2 (2021.8) to generate a sample metadata file (Ding et al., 2021). DADA2 was used for denoising analysis to obtain corresponding feature tables and ASVs data. A Venn diagram was constructed to display the diversity information of ASVs data, showing the abundance of different species and groups. Beta diversity analysis was performed to reveal differences in community structure among different sample groups using principal coordinates analysis (PCoA) and orthogonal partial least squares discriminant analysis (OPLS-DA). Linear discriminant analysis effect size (LEfSe) was used to identify significant differences in microbial composition between different groups (Xu et al., 2022).
2.8. Untargeted metabolomics
In untargeted metabolomics experiments, an X500R HPLC-Q-TOF/MS (SCIEX) was used with an XSelect Premier HSS T3 column (100 mm × 2.1 mm, Waters Corporation, Milford, Massachusetts, USA) for chromatographic separation (Jing et al., 2021). The mobile phase consisted of 0.1% formic acid in water (solvent A) and acetonitrile (solvent B), with a gradient duration of 18 min and a flow rate of 0.3 mL/min (Mi et al., 2023). ESI was used for ionization in both positive and negative ion modes, with ion source gasses 1 and 2 set at 55 psi and curtain gas at 35 psi. The ion source temperature was set to 550 C, while the positive and negative ion source voltages were set to 5,000 V and − 4,500 V, respectively. TOF/MS scans were performed over a m/z range of 50–1,200 Da, while TOF/MS/MS scans were performed over a m/z range of 30–1,200 Da. MS/MS data was detected using information-dependent acquisition (IDA).
Data acquisition was controlled using SCIEX OS 2.0 software (AB SCIEX). Metabolite annotation was performed using MetDNA and SCIEX OS 2.0 software. Statistical analysis was conducted using unsupervised principal component analysis (PCA) and supervised orthogonal partial least squares discriminant analysis (OPLS-DA) in MetaboAnalyst 5.010 to identify differences in the studied metabolites and screen for potential biomarkers (Zhang Q. et al., 2021). After verifying the differential metabolites in the public database HMDB, we performed metabolic pathway analysis based on the KEGG database (Zhou et al., 2021).
2.9. Statistical analysis
The data collected in this study were analyzed by using SPSS statistical software (24.0). One-way analysis of variance (ANOVA) was used to compare mean differences among the groups, with significance set at p < 0.05. The results were expressed as mean ± SD.
3. Results
3.1. Puerarin is a potential ingredient that restored the growth performance of colitis mice
A total of 8 potential active components have been screened from Radix Puerariae (detailed information can be found in Supplementary Table S2) and identified 46 overlapping genes between core disease targets (1090) and active ingredient targets (113). Using Cytoscape software, we constructed a network diagram of the ingredient-disease-target interactions (Figure 1A), which revealed that genistein, puerarin, and daidzein had the highest degree of overlap with disease targets. We also performed protein interaction analysis on the overlapping genes and identified key genes such as IL1β, TNF, TP53, and INS (Figure 1B), and are involved in TNF, PI3K-Akt, and VEGF signaling pathways (Supplementary Figures S1–S3).
Subsequent molecular docking revealed that puerarin had the most stable docking with the core target (Figure 2; Supplementary Table S3). And considering the low OB and DL values of genistein (OB = 17.93%, DL = 0.21), we prioritized the use of puerarin, which showed similar good effect, for subsequent animal experiments.
Figure 2. Docking of active ingredients in Radix Puerariae with target molecules. (A) Schematic diagram of docking of puerarin with TNF-α. (B) Schematic diagram of docking of Puerarin with IL1B. (C) Schematic diagram of docking of Puerarin with IL10.
As shown in Figures 1D,E, mice that were orally administered DSS in the first week exhibited a significant decrease in body weight and food intake compared to the control group (p < 0.01). However, after stopping the modeling and supplementing with puerarin or sulfasalazine, the body weight of mice gradually recovered, and their food intake also increased. After 2 weeks of oral administration, the food intake of mice in the modeling group gradually returned to normal levels and approached those of mice in the control group.
3.2. Puerarin alleviates the symptoms of DSS induced colitis in mice
We performed DAI scoring on mice at the end of both modeling and treatment (Figure 3A). Representative images of fecal samples from mice at the end of modeling/treatment are shown in Supplementary Figure S4A. The control group showed no abnormalities throughout the experiment, while severe bloody stools were observed in the other three groups after modeling. Some mice in the UC group still had mild red stool, while some mice in the PUE group had mild soft stools without blood and had returned to normal status.
Figure 3. Puerarin improved disease score, colonic tissue morphology and reduced inflammatory factor expression in mice. (A) DAI score of each group at the end of modeling and administration. (B) Colon length of mice in each group. *p < 0.05, compared to the CON group; #p < 0.05, compared to the UC group. (C) HE section of colon of mice in each group. (D) The expression levels of IL-1β, IL-6 and TNF-α in each group (n = 8). **p < 0.01, compared to the CON group; ##p < 0.01, compared to the UC group; &&p < 0.01, compared to the SASP group.
Observation of colon tissue morphology revealed that DSS-induced colitis significantly shortened colon length (p < 0.05), while supplementation with puerarin significantly increased colon length compared to the UC group (p < 0.05), and showed no difference compared to the control and SASP groups (Figure 3B; Supplementary Figure S4B).
3.3. Puerarin improved colonic tissue morphology and reduced inflammatory factor expression
Figure 3C shows the HE staining of intestinal sections from mice in each group. In the control group, the colon tissue appeared normal, with intact crypt epithelial cells, goblet cells, and submucosa. However, in the UC group induced by DSS, ulcers appeared in the colon and spread to the mucosa and submucosa. The mucosal layer showed infiltration of inflammatory cells, and there was significant damage to the crypts with structural deformation. In contrast, mice treated with PUE did not show obvious colonic ulcers, and only a small amount of inflammatory cell infiltration was observed in the mucosal layer. The deformation of crypts also improved.
As shown in Figure 3D, compared to the CON group, the concentrations of IL-1β, IL-6 and TNF-α were significantly increased (p < 0.01) in the UC group. This suggests that there was significant inflammation in the intestines of UC mice. After treatment with SASP or PUE, the concentrations of these three cytokines significantly decreased (p < 0.01), but did not fully recover to normal levels (vs CON). Notably, PUE showed a more significant decrease in IL-1β and IL-6 compared to SASP, indicating that puerarin has a better effect on improving inflammation caused by DSS-induced cytokine surge.
3.4. Effect of puerarin on intestinal microorganisms in mice with colitis
By examining the species accumulation curve in microbial analysis (Supplementary Figure S5) and conducting α-diversity, β-diversity analysis on different microbial groups, the sample size and community richness were evaluated. The PCoA (Figure 4A) and NMDS plot (Figure 4B) showed significant separation between the UC group and the control group, indicating that DSS significantly altered the composition of the microbiota in mice, which was distinct from that of the control group. Moreover, there was no overlap between the PUE group, SASP group and UC group, suggesting that administration induced significant changes in gut microbiota composition with a trend toward normalizing to that of the control group.
Figure 4. Changes of microbiota composition in mice after puerarin intervention. (A) PcoA diagram. (B) NMDS chart. (C) ASV Venn diagram of mice in each group. (D) Comparison of α diversity between UC and PUE groups. (E) Species composition at phylum level. (F) Species composition at genus level.
A Venn diagram was constructed to characterize the overlap of ASVs among the microbiota based on the ASV information table obtained from QIIME2 denoising analysis. As shown in Figure 4C, the four groups contained 264 (CON), 97 (UC), 147 (SASP), and 237 (PUE) unique ASVs respectively, with a total of 105 ASVs shared as their core microbiota. The results indicated that induced colitis significantly reduced gut species composition, while the number of ASVs increased after SASP intervention and was significantly higher after PUE intervention.
Figure 4D displays the changes in α-diversity of gut microbiota after modeling and intervention. It was observed that the Shannon, ACE, and Chao indices of mice in the UC group were significantly lower than those in the PUE group, indicating that puerarin intervention partially restored the species composition and abundance of gut microbiota in mice.
Furthermore, the top 5 species at the phylum level and top 10 species at the genus level were analyzed (Figures 4E,F). The major microbiota included Firmicutes, Bacteroidetes, Proteobacteria, Deferribacteres etc. DSS modeling increased the level of Bacteroidetes and decreased the level of Firmicutes. However, compared to the UC group, the PUE group reversed the change in F/B ratio and showed a trend toward normalization to that of the control group.
The major species at the genus level included Allobaculum, S24-7, Lactobacillus, Prevotella, Bacteroides, and Turicibacter. Comparing the changes in gut microbiota at the genus level revealed that DSS modeling significantly increased the levels of Prevotella, S24-7, Bacteroides and Mucispirillum, while decreased the levels of Allobaculum and Lactobacillus. Compared to the UC group, PUE intervention reversed the changes in Allobaculum, Prevotella, Bacteroides, and Mucispirillum levels while decreased S24-7 levels and moving toward a microbiota structure similar to that of the control group.
In addition, we used linear discriminant analysis effect size (LEfSe) to examine the differences in fecal species abundance between two groups of mice at different levels (Figures 5A,B). A total of 26 differentially abundant species were identified, including 11 in the UC group: Prevotella, S24-7, Bacteroides, Veillonella, Veillonellaceae, Listeria, Eysipelotrichaceae, Ovatus, Alistipes, Copri and Enterobacteriaceae. The PUE group had 6 differentially abundant species: Rhodothermi, Alteromonadaceae, Halomonas, Halomonadaceae, Marinobacter and Bilophila. The SASP intervention had 6 significant species: Turicibacter, Turicibacteraceae, Turicibacterales, Reuteri, Parabacteroides, and Johnsonii.
Figure 5. (A) Evolutionary branching diagram. (B) LEfSe multi-level species discrimination. (C) Correlation network of gut microbiota in PUE Group. (D) Correlation network of gut microbiota in UC Group. The red lines indicate negative correlation, while the blue lines indicate positive correlation. The red circles represent the bacteria with a larger number of associations. (E,F) The heat map of correlation analysis between inflammatory factors and significantly different microorganisms. * indicates p < 0.05; ** indicates p < 0.01.
In order to further elucidate the potential relationships among gut microbiota, we conducted a correlation analysis of the gut microbiota in UC and PUE mice. The resulting visualization network diagram is shown in Figures 5C,D. Our analysis revealed that in the UC state, the abundance of Allobaculum, which is notably reduced in the UC state, exhibits a significant negative correlation with the dominant bacterium Prevotella. While various pathogenic bacteria such as Alistipes, Listeria, and Rikenellaceaeg formed a tightly connected network with multiple bacteria. After puerarin intervention, we observed a significant enrichment of Halomonas, while Bilophila and Halomonas were negatively correlated with Prevotella. Additionally, Allobaculum was negatively correlated with Marinobacter and Marinimicrobium. These correlation analysis results indicate that there are extensive connections among gut microbiota, and the dominance of different bacteria in different states may be due to their interactions with each other.
Similarly, to understand the relationship between differentially abundant species and inflammatory factors, we conducted a correlation analysis between inflammatory factors and significantly different species (Figures 5E,F). We found that in the UC group, Enterobacteriaceae, Alistipes and Listeria (p < 0.01), Copri, Veillonellaceae and S24_7 (p < 0.05) were positively correlated with the concentration of pro-inflammatory factors, while Eysipelotrichaceae (p < 0.01) was negatively correlated with them. In the PUE group, Alteromonadaceae and Halomonadaceae (p < 0.05) were negatively correlated with the concentration of pro-inflammatory factors, while Balneola and Marinobacter were negatively correlated with TNF-α (p < 0.05).
3.5. Effect of puerarin on intestinal metabolism in mice with colitis
In unsupervised mode (Figure 6A), the UC group was clearly separated from the control group, indicating significant differences in fecal metabolites between the two groups; PUE and SASP interventions caused their metabolite characteristics to move away from the UC group and toward the control group. At the same time, in the supervised OPLS-DA mode, each group could also be well separated (Figures 6B–D): UC vs. control group (R2 = 0.618, Q2 = 0.342), SASP vs. UC group (R2 = 0.793, Q2 = 0.364), PUE vs. UC group (R2 = 0.602, Q2 = 0.35). The results indicate that there are differences in metabolites in feces of each group.
Figure 6. Changes of metabolite levels in each group. (A) PCA diagram between four groups. (B–D) OPLS-DA diagram between two groups. CON vs. UC; UC vs. SASP; UC vs. PUE. (E) Volcano plot of untargeted differential metabolites in UC mice compared to control group. (F) Volcano plot of differential metabolites in SASP mice compared to UC group. (G) Volcano plot of differential metabolites in PUE mice compared to UC group. (H) Bubble plot of differential metabolic pathways analysis in UC mice compared to control group.
By setting the conditions of FC >1 or FC <1/2 and p < 0.1, differentially abundant metabolites were screened out, and volcano plots were used to represent the differentially abundant metabolites obtained from different groups (Figures 6E–G). From the figure, it can be seen that compared to the CON group, there were 13 metabolites (such as 3-hydroxybenzoic acid, docosahexaenoic acid, indolelactic acid, etc.) increased and 18 metabolites decreased (such as azelaic acid, sebacic acid, uridine, etc.) in the UC group. Compared to the UC group, there were 7 metabolites increased and 1 metabolite (L-methionine) decreased in the SASP group; while there were 17 metabolites increased and 1 metabolite (docosahexaenoic acid) decreased in the PUE group. Specific information on differentially abundant metabolites can be found in Supplementary Tables S4–S6.
Further setting VIP value >1, 26 differentially abundant metabolites obtained from the secondary screening of the UC group were subjected to metabolic pathway analysis (Figure 6H), and it was found that these metabolites had a major impact on 7 metabolic pathways, including the biosynthesis of phenylalanine, tyrosine and tryptophan, D-glutamine and D-glutamate metabolism, taurine and hypotaurine metabolism, phenylalanine metabolism, alanine, aspartate and glutamate metabolism, arginine biosynthesis and cysteine and methionine metabolism.
To investigate the relationship between differentially abundant metabolites and differentially abundant species in gut microbiota, we performed correlation clustering analysis, as shown in Figure 7A. Enterobacteriaceae and Listeria were significantly negatively correlated with undecanedioic, sebacic acid, azelacic acid, sebaric acid, guanosine, and uridine. While bilophila was positively correlated with taurine, 4-Hydroxy-3-methylbenzoic acid and negtively correlated with various fatty acids levels. On the other hand, the enrichment of several bacterial taxa may affect the levels of 3-Hydroxybenzoic acid, azelaic acid and sebacic acid.
Figure 7. (A) Correlation analysis chart of gut microbiota and fecal metabolites. * indicates p < 0.05; ** indicates p < 0.01. (B) Prediction and analysis of gut microbiota function in mice.
To gain a deeper understanding of the roles and functions of gut microbiota in mice, we conducted a comprehensive functional prediction analysis (Figure 7B). Our findings revealed that these gut microbiota are widely involved in various metabolic activities, including carbohydrate metabolism, amino acid metabolism, metabolism of cofactors and vitamins etc. These results reveal that gut microbiota play a significant role in maintaining the host’s homeostasis and health, and this impact is closely related to their metabolic activities.
4. Discussion
Nowadays, there is a widespread and profound understanding of the role of natural products in disease intervention. Due to the large number and variety of small molecules, it is still necessary to seek more efficient and convenient methods to explore potential effective small molecule substances. Network pharmacology, which utilizes multiple databases and integrates various bioinformatics analyses, can effectively aid in addressing these issues (Wang et al., 2021).
In this study, we first utilized network pharmacology techniques to analyze the active ingredients in Radix Puerariae, and based on GO and PPI analysis, core protein targets for ulcerative colitis were identified, and molecular docking was performed between active ingredients and targets, resulting in the identification of puerarin as the most promising potential small molecule substance. We confirmed that it can effectively improve the pathological symptoms of DSS-induced colitis in mice, including restoring mouse weight and colon length, inhibiting intestinal tissue damage, and reducing DAI scores. In our experiments, the decrease in concentrations of IL-1β, IL-6 and TNF-α indicated that puerarin has an inhibitory effect on intestinal inflammation. Previous studies have also shown that puerarin can promote goblet cell secretion of mucin to rebuild the mucus layer, forming a barrier between bacteria and epithelial cells to alleviate colitis (Wu et al., 2021).
The microbial community inhabiting the gut plays an essential role in maintaining intestinal health. This study revealed that UC mice exhibited a reduction in Firmicutes levels and an increase in Bacteroidetes levels within their gut microbiota. The change in the F/B ratio is considered a manifestation of gut dysbiosis. Firmicutes are known to produce more energy from equivalent amounts of food compared to other bacterial groups, thus their reduction may lead to weight loss observed in colitis mice. Conversely, Bacteroidetes exhibit pro-inflammatory properties due to their lipopolysaccharides (Shen et al., 2018). The observed decrease in F/B ratio is closely associated with inflammatory bowel disease (Stojanov et al., 2020).
In DSS-induced colitis mice, Alistipes, P. copri, Enterobacteriaceae, Veillonella and other bacterial groups are significantly enriched. Studies have shown that they can enhance susceptibility to colorectal cancer and arthritis by interfering with IL-6/STAT3 signaling pathway and bile acid metabolism (David et al., 2014; Hofer, 2014; Parker et al., 2020). The differential bacteria Listeria has a higher infectivity for immunodeficient patients and also has a higher pathogenicity in IBD (Azimi et al., 2018). Recent studies have also suggested that Enterobacteriaceae can use mucosal inflammation to gain a competitive advantage in obtaining energy among symbiotic bacteria, thereby enriching and exacerbating the development of colitis (Garrett et al., 2010). In addition, our correlation analysis between inflammatory factors and differential species also shows a highly positive correlation between these inflammatory pathogenic bacteria and the concentrations of IL-1β, IL-6 and TNF-α. All of these indicate that the enrichment of these pathogenic and pro-inflammatory microorganisms as significant differential species in the UC group is a key factor leading to dysbiosis in the gut microbiota of colitis mice.
Compared with the UC group, after treatment with puerarin, the F/B ratio was restored, inflammation subsided, and the dominant pro-inflammatory bacterial groups were no longer present. Furthermore, it further reversed the significant changes in abundance of bacterial groups in the model group, such as increasing the level of Desulfovibrionaceae, which is also a type of bacteria that performs well in anti-inflammatory and anti-obesity aspects (Popli et al., 2022).
For better understanding the role of these bacteria in the gut and their interrelationships, we first conducted Spearman correlation analysis on the microbiota of mice in the model group and the puerarin intervention group. We found that the bacteria formed a tightly connected network. Interestingly, some bacteria that significantly decreased in disease, such as Allobaculum, showed significant negative correlations with pathogenic bacteria such as Alistipes and Listeria. This suggests that in a diseased state, bacteria may resist the growth of other bacterial species by taking advantage of their ecological niche, ultimately resulting in increased abundance, while bacteria negatively correlated with them decrease in abundance. In addition, using PICRUSt to predict the function of the mice microbiota, we found that the top-ranking abundance were involved in various metabolic activities, such as carbohydrate metabolism, amino acid metabolism, and metabolism of cofactors and vitamins. These results all indicate that the gut microbiota, as a special type of organism within the body, has a crucial impact on the host due to its close interrelationships and rich metabolic activities.
Recently, there have been continuous discussions about how metabolites from gut microbiota affect immune homeostasis, host energy metabolism and intestinal mucosal integrity in inflammatory bowel disease (Lavelle and Sokol, 2020). Our previous study found that indole-3-propionic acid derived from the tryptophan pathway can enhance mucous barrier by increasing secretion of mucin (MUC2 and MUC4) and goblet cell products (Li J. et al., 2021). Song et al. also reported that 3-oxo-deoxycholic acid, a bile acid metabolite, can regulate Foxp3 Treg cell homeostasis in the intestinal lamina propria to exert anti-inflammatory properties (Song et al., 2020). Although the oral bioavailability of screened puerarin did not reach ideal standards, recent reports suggest that these polyphenolic small molecules can remain in the intestinal lumen and participate in gut microbiota metabolism to exert their functions (Luca et al., 2020).
In this experiment, it was found that levels of fecal metabolites such as glutamate and 4-pyridoxic acid decreased while taurine levels significantly increased in UC mice. Taurine is known to activate NLRP6 inflammasomes leading to inflammatory reactions in the intestinal mucosal layer and activation of vitamin B6 degradation pathways during cellular immune responses (Levy et al., 2015). However, no changes were observed in these metabolic pathways after treatment with puerarin. Additionally, tryptophan was not detected but its secondary metabolite indole derivative significantly increased in DSS-induced colitis mice.
Further metabolic pathway analysis of UC mice showed significant changes in the tryptophan metabolism pathway. In recent studies, it has similarly been found that DSS-induced colitis disrupts tryptophan/kynurenine metabolism through the regulation of the rate-limiting enzyme IDO1, and this phenomenon is closely associated with alterations in the gut microbiota (Zhao et al., 2023). Differential species analysis of fecal microbiota also revealed a significant increase in the abundance of the Bacteroides and Enterobacteriaceae in the UC group. Previous reports have suggested that Bacteroides is involved in many important metabolic activities in the colon, including carbohydrate fermentation, nitrogenous substance utilization, and biotransformation of bile acids and other steroids. Some species of Bacteroides can produce indole-3-lactic acid, indole-3-acetic acid, and tryptamine (Liu Y. et al., 2023; Nie et al., 2023). Enterobacteriaceae have also been reported to have the ability to metabolize tryptophan in vivo (Zhang J. et al., 2021). These findings suggest that there may be an imbalance and excessive utilization of tryptophan in the UC mice. Therefore, metabolic disorders such as the indole pathway, urea cycle, glutamate, vitamins, and taurine may exacerbate the development of intestinal inflammation in colitis mice.
After administration of puerarin, the abundance of the Bacteroides in the gut was found to be decreased, and the metabolic levels of 3-HAA, picolinic acid, and niacin metabolites were reversed. Among them, 3-HAA is a type of anti-inflammatory tryptophan metabolite, an increase in its concentration can inhibit the development of inflammation in collagen-induced arthritis models. Other studies have shown that 3-HAA can induce apoptosis in liver cancer cells as a ligand-activated transcription factor YY1 and has anti-cancer effects (Kolodziej, 2013; Shi et al., 2021). Puerarin may have an inhibitory effect on intestinal inflammation by regulating the tryptophan metabolism pathway. Yang et al. found that the active ingredient, theabrownin, in Fu Brick Tea can effectively correct intestinal damage in UC mice and promote the secretion of indole derivatives, which are downstream products of the tryptophan metabolic pathway (Yang et al., 2023). This effective improvement of the disease by alleviating the dysregulation of tryptophan metabolism in the microbiota was also observed after turmeric polysaccharides intervention (Yang et al., 2021).
In our results, the differential bacterial groups Rhodobacteraceae and Halomonadaceae after puerarin treatment contain enzymes that break down tryptophan into kynurenine and downstream metabolites (Vujkovic-Cvijin et al., 2013). Therefore, we speculate that the improvement effect of puerarin on UC may be partially achieved by restoring the function of gut microbiota and thus restoring tryptophan metabolism.
In addition, we correlated differential species and metabolites and found that Enterobacteriaceae and Listeria were significantly and negatively correlated with undecanedioic, sebacic acid, azelacic acid, sebaric acid, guanosine, and uridine, which are pathogenic bacteria that dominate under disease conditions. After treatment, we also observed a significant increase in purine metabolites such as uridine and guanosine, which have been reported in previous studies to activate the PPARγ signaling pathway in human colonic epithelial cells and restore mucosal barrier function (Li D. et al., 2021). These results well demonstrate that puerarin alleviates DSS-induced colitis by restoring the abundance and metabolic advantage of specific bacterial communities in mice.
Overall, our research indicates that puerarin intervention effectively reduced pathological damage and inflammatory status in DSS-induced colitis mice, and partially restored gut microbiota and metabolic disorders. Our study expands on the potential effects of puerarin in improving intestinal inflammation based on previous research, providing data support for the exploration and development of such natural small molecule substances in the future.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Ethics statement
The animal study was approved by the Animal Experimentation Ethical Committee of Nanchang University (License number: SYXK 2021-0001). The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
YZ: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft. WD: Conceptualization, Data curation, Methodology, Visualization, Writing – original draft. YW: Methodology, Software, Supervision, Writing – original draft. TC: Data curation, Formal analysis, Investigation, Writing – original draft. ZR: Funding acquisition, Resources, Supervision, Validation, Writing – review & editing, Project administration.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the National Natural Science Foundation of China (32060539).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1279029/full#supplementary-material
Footnotes
1. ^ https://old.tcmsp-e.com/tcmsp.php
4. ^ https://www.disgenet.org/home/
5. ^ https://www.genecards.org/
7. ^ https://DAVID.ncifcrf.gov/
References
Azimi, T., Nasiri, M. J., Chirani, A. S., Pouriran, R., and Dabiri, H. (2018). The role of bacteria in the inflammatory bowel disease development: a narrative review. APMIS 126, 275–283. doi: 10.1111/apm.12814
David, L. A., Maurice, C. F., Carmody, R. N., Gootenberg, D. B., Button, J. E., Wolfe, B. E., et al. (2014). Diet rapidly and reproducibly alters the human gut microbiome. Nature 505, 559–563. doi: 10.1038/nature12820
Ding, Y., Li, X., Liu, Y., Wang, S., and Cheng, D. (2021). Protection mechanisms underlying oral administration of chlorogenic acid against cadmium-induced hepatorenal injury related to regulating intestinal flora balance. J. Agric. Food Chem. 69, 1675–1683. doi: 10.1021/acs.jafc.0c06698
Garrett, W. S., Gallini, C. A., Yatsunenko, T., Michaud, M., DuBois, A., Delaney, M. L., et al. (2010). Enterobacteriaceae act in concert with the gut microbiota to induce spontaneous and maternally transmitted colitis. Cell Host Microbe 8, 292–300. doi: 10.1016/j.chom.2010.08.004
Hofer, U. (2014). Microbiome: pro-inflammatory Prevotella? Nat. Rev. Microbiol. 12:5. doi: 10.1038/nrmicro3180
Huang, Y. F., Zhou, J. T., Qu, C., Dou, Y. X., Huang, Q. H., Lin, Z. X., et al. (2017). Anti-inflammatory effects of Brucea javanica oil emulsion by suppressing NF-κB activation on dextran sulfate sodium-induced ulcerative colitis in mice. J. Ethnopharmacol. 198, 389–398. doi: 10.1016/j.jep.2017.01.042
Jiao, W., Mi, S., Sang, Y., Jin, Q., Chitrakar, B., Wang, X., et al. (2022). Integrated network pharmacology and cellular assay for the investigation of an anti-obesity effect of 6-shogaol. Food Chem. 374:131755. doi: 10.1016/j.foodchem.2021.131755
Jing, W., Dong, S., Luo, X., Liu, J., Wei, B., Du, W., et al. (2021). Berberine improves colitis by triggering AhR activation by microbial tryptophan catabolites. Pharmacol. Res. 164:105358. doi: 10.1016/j.phrs.2020.105358
Kolodziej, L. (2013). Systemic metabolism of tryptophan and its catabolites, kynurenine and 3-HAA, in mice with inflammatory arthritis. Gene 512, 23–27. doi: 10.1016/j.gene.2012.09.122
Kucharzik, T., Koletzko, S., Kannengiesser, K., and Dignass, A. (2020). Ulcerative colitis-diagnostic and therapeutic algorithms. Dtsch. Arztebl. Int. 117, 564–574. doi: 10.3238/arztebl.2020.0564
Lavelle, A., and Sokol, H. (2020). Gut microbiota-derived metabolites as key actors in inflammatory bowel disease. Nat. Rev. Gastroenterol. Hepatol. 17, 223–237. doi: 10.1038/s41575-019-0258-z
Le Berre, C., Honap, S., and Peyrin-Biroulet, L. (2023). Ulcerative colitis. Lancet 402, 571–584. doi: 10.1016/s0140-6736(23)00966-2
Levy, M., Thaiss, C. A., Zeevi, D., Dohnalová, L., Zilberman-Schapira, G., Mahdi, J. A., et al. (2015). Microbiota-modulated metabolites shape the intestinal microenvironment by regulating NLRP6 inflammasome signaling. Cells 163, 1428–1443. doi: 10.1016/j.cell.2015.10.048
Li, R., Chen, Y., Shi, M., Xu, X., Zhao, Y., Wu, X., et al. (2016). Gegen Qinlian decoction alleviates experimental colitis via suppressing TLR4/NF-κB signaling and enhancing antioxidant effect. Phytomedicine 23, 1012–1020. doi: 10.1016/j.phymed.2016.06.010
Li, D., Feng, Y., Tian, M., Ji, J., Hu, X., and Chen, F. (2021). Gut microbiota-derived inosine from dietary barley leaf supplementation attenuates colitis through PPARγ signaling activation. Microbiome 9:83. doi: 10.1186/s40168-021-01028-7
Li, X., Wei, S., Niu, S., Ma, X., Li, H., Jing, M., et al. (2022). Network pharmacology prediction and molecular docking-based strategy to explore the potential mechanism of Huanglian Jiedu decoction against sepsis. Comput. Biol. Med. 144:105389. doi: 10.1016/j.compbiomed.2022.105389
Li, J., Zhang, L., Wu, T., Li, Y., Zhou, X., and Ruan, Z. (2021). Indole-3-propionic acid improved the intestinal barrier by enhancing epithelial barrier and mucus barrier. J. Agric. Food Chem. 69, 1487–1495. doi: 10.1021/acs.jafc.0c05205
Liu, J., Di, B., and Xu, L. L. (2023). Recent advances in the treatment of IBD: targets, mechanisms and related therapies. Cytokine Growth Factor Rev. 71-72, 1–12. doi: 10.1016/j.cytogfr.2023.07.001
Liu, Y., Pei, Z., Pan, T., Wang, H., Chen, W., and Lu, W. (2023). Indole metabolites and colorectal cancer: gut microbial tryptophan metabolism, host gut microbiome biomarkers, and potential intervention mechanisms. Microbiol. Res. 272:127392. doi: 10.1016/j.micres.2023.127392
Luca, S. V., Macovei, I., Bujor, A., Miron, A., Skalicka-Woźniak, K., Aprotosoaie, A. C., et al. (2020). Bioactivity of dietary polyphenols: the role of metabolites. Crit. Rev. Food Sci. Nutr. 60, 626–659. doi: 10.1080/10408398.2018.1546669
Mi, Y., Xu, X., Zeng, F., Li, N., Tan, X., and Gong, Z. (2023). Prebiotics alleviates cartilage degradation and inflammation in post-traumatic osteoarthritis mice by modulating the gut barrier and fecal metabolomics. Food Funct. 14, 4065–4077. doi: 10.1039/D3FO00775h
Neurath, M. F. (2014). Cytokines in inflammatory bowel disease. Nat. Rev. Immunol. 14, 329–342. doi: 10.1038/nri3661
Nie, Q., Sun, Y., Li, M., Zuo, S., Chen, C., Lin, Q., et al. (2023). Targeted modification of gut microbiota and related metabolites via dietary fiber. Carbohydr. Polym. 316:120986. doi: 10.1016/j.carbpol.2023.120986
Nogales, C., Mamdouh, Z. M., List, M., Kiel, C., Casas, A. I., and Schmidt, H. (2022). Network pharmacology: curing causal mechanisms instead of treating symptoms. Trends Pharmacol. Sci. 43, 136–150. doi: 10.1016/j.tips.2021.11.004
Parker, B. J., Wearsch, P. A., Veloo, A. C. M., and Rodriguez-Palacios, A. (2020). The genus Alistipes: gut Bacteria with emerging implications to inflammation, Cancer, and mental health. Front. Immunol. 11:906. doi: 10.3389/fimmu.2020.00906
Peng, L., Gao, X., Nie, L., Xie, J., Dai, T., Shi, C., et al. (2020). Astragalin attenuates dextran sulfate sodium (DSS)-induced acute experimental colitis by alleviating gut microbiota Dysbiosis and inhibiting NF-κB activation in mice. Front. Immunol. 11:2058. doi: 10.3389/fimmu.2020.02058
Popli, S., Badgujar, P. C., Agarwal, T., Bhushan, B., and Mishra, V. (2022). Persistent organic pollutants in foods, their interplay with gut microbiota and resultant toxicity. Sci. Total Environ. 832:155084. doi: 10.1016/j.scitotenv.2022.155084
Shen, Z. H., Zhu, C. X., Quan, Y. S., Yang, Z. Y., Wu, S., Luo, W. W., et al. (2018). Relationship between intestinal microbiota and ulcerative colitis: mechanisms and clinical application of probiotics and fecal microbiota transplantation. World J. Gastroenterol. 24, 5–14. doi: 10.3748/wjg.v24.i1.5
Shi, Z., Gan, G., Xu, X., Zhang, J., Yuan, Y., Bi, B., et al. (2021). Kynurenine derivative 3-HAA is an agonist ligand for transcription factor YY1. J. Hematol. Oncol. 14:153. doi: 10.1186/s13045-021-01165-4
Song, X., Sun, X., Oh, S. F., Wu, M., Zhang, Y., Zheng, W., et al. (2020). Microbial bile acid metabolites modulate gut RORγ(+) regulatory T cell homeostasis. Nature 577, 410–415. doi: 10.1038/s41586-019-1865-0
Stojanov, S., Berlec, A., and Štrukelj, B. (2020). The influence of probiotics on the Firmicutes/Bacteroidetes ratio in the treatment of obesity and inflammatory bowel disease. Microorganisms 8:1715. doi: 10.3390/microorganisms8111715
Vujkovic-Cvijin, I., Dunham, R. M., Iwai, S., Maher, M. C., Albright, R. G., Broadhurst, M. J., et al. (2013). Dysbiosis of the gut microbiota is associated with HIV disease progression and tryptophan catabolism. Sci. Transl. Med. 5:193ra91. doi: 10.1126/scitranslmed.3006438
Wang, D., Bu, T., Li, Y., He, Y., Yang, F., and Zou, L. (2022). Pharmacological activity, pharmacokinetics, and clinical research progress of puerarin. Antioxidants 11:2121. doi: 10.3390/antiox11112121
Wang, X., Gao, Y., Wang, L., Yang, D., Bu, W., Gou, L., et al. (2021). Troxerutin improves dextran sulfate sodium-induced ulcerative colitis in mice. J. Agric. Food Chem. 69, 2729–2744. doi: 10.1021/acs.jafc.0c06755
Wang, Y., Zhang, J., Zhang, B., Lu, M., Ma, J., Liu, Z., et al. (2023). Modified Gegen Qinlian decoction ameliorated ulcerative colitis by attenuating inflammation and oxidative stress and enhancing intestinal barrier function in vivo and in vitro. J. Ethnopharmacol. 313:116538. doi: 10.1016/j.jep.2023.116538
Wu, Y., Li, J., Ding, W., Ruan, Z., and Zhang, L. (2021). Enhanced intestinal barriers by puerarin in combination with tryptophan. J. Agric. Food Chem. 69, 15575–15584. doi: 10.1021/acs.jafc.1c05830
Xu, D., Zhuang, L., Gao, S., Ma, H., Cheng, J., Liu, J., et al. (2022). Orally administered Ginkgolide C attenuates DSS-induced colitis by maintaining gut barrier integrity, inhibiting inflammatory responses, and regulating intestinal Flora. J. Agric. Food Chem. 70, 14718–14731. doi: 10.1021/acs.jafc.2c06177
Yan, X., Jin, J., Su, X., Yin, X., Gao, J., Wang, X., et al. (2020). Intestinal Flora modulates blood pressure by regulating the synthesis of intestinal-derived corticosterone in high salt-induced hypertension. Circ. Res. 126, 839–853. doi: 10.1161/circresaha.119.316394
Yang, C., Du, Y., Ren, D., Yang, X., and Zhao, Y. (2021). Gut microbiota-dependent catabolites of tryptophan play a predominant role in the protective effects of turmeric polysaccharides against DSS-induced ulcerative colitis. Food Funct. 12, 9793–9807. doi: 10.1039/d1fo01468d
Yang, W., Ren, D., Shao, H., Zhang, X., Li, T., Zhang, L., et al. (2023). Theabrownin from Fu brick tea improves ulcerative colitis by shaping the gut microbiota and modulating the tryptophan metabolism. J. Agric. Food Chem. 71, 2898–2913. doi: 10.1021/acs.jafc.2c06821
Zhang, R., Zhang, Q., Chen, Y., Zhao, Q., Zhang, B., Wang, L., et al. (2022). Combined treatment with Rg1 and adipose-derived stem cells alleviates DSS-induced colitis in a mouse model. Stem Cell Res Ther 13:272. doi: 10.1186/s13287-022-02940-x
Zhang, H., Zhang, Y., Li, Y., Wang, Y., Yan, S., Xu, S., et al. (2021). Bioinformatics and network pharmacology identify the therapeutic role and potential mechanism of melatonin in AD and Rosacea. Front. Immunol. 12:756550. doi: 10.3389/fimmu.2021.756550
Zhang, Q., Zhong, D., Ren, Y. Y., Meng, Z. K., Pegg, R. B., and Zhong, G. (2021). Effect of konjac glucomannan on metabolites in the stomach, small intestine and large intestine of constipated mice and prediction of the KEGG pathway. Food Funct. 12, 3044–3056. doi: 10.1039/d0fo02682d
Zhang, J., Zhu, S., Ma, N., Johnston, L. J., Wu, C., and Ma, X. (2021). Metabolites of microbiota response to tryptophan and intestinal mucosal immunity: a therapeutic target to control intestinal inflammation. Med. Res. Rev. 41, 1061–1088. doi: 10.1002/med.21752
Zhao, Y., Luan, H., Jiang, H., Xu, Y., Wu, X., Zhang, Y., et al. (2021). Gegen Qinlian decoction relieved DSS-induced ulcerative colitis in mice by modulating Th17/Treg cell homeostasis via suppressing IL-6/JAK2/STAT3 signaling. Phytomedicine 84:153519. doi: 10.1016/j.phymed.2021.153519
Zhao, L. P., Wu, J., Quan, W., Zhou, Y., Hong, H., Niu, G. Y., et al. (2023). DSS-induced acute colitis causes dysregulated tryptophan metabolism in brain: an involvement of gut microbiota. J. Nutr. Biochem. 115:109282. doi: 10.1016/j.jnutbio.2023.109282
Zhou, R., He, D., Xie, J., Zhou, Q., Zeng, H., Li, H., et al. (2021). The synergistic effects of polysaccharides and Ginsenosides from American ginseng (Panax quinquefolius L.) ameliorating cyclophosphamide-induced intestinal immune disorders and gut barrier dysfunctions based on microbiome-metabolomics analysis. Front. Immunol. 12:665901. doi: 10.3389/fimmu.2021.665901
Keywords: puerarin, ulcerative colitis, network pharmacology, gut microbiota, untargeted metabonomics
Citation: Zou Y, Ding W, Wu Y, Chen T and Ruan Z (2023) Puerarin alleviates inflammation and pathological damage in colitis mice by regulating metabolism and gut microbiota. Front. Microbiol. 14:1279029. doi: 10.3389/fmicb.2023.1279029
Edited by:
Huaxi Yi, Ocean University of China, ChinaCopyright © 2023 Zou, Ding, Wu, Chen and Ruan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zheng Ruan, ruanzheng@ncu.edu.cn