Xinyue Li1,2†
Xiaochen Wang2,3†Ziming Wang2,3†Mingyi Zhang1Song Wang4
Zuofu Xiang5
Huijuan Pan1*
Ming Li2,6- 1School of Ecology and Nature Conservation, Beijing Forestry University, Beijing, China
- 2CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Beijing, China
- 3College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China
- 4Nanning Zoo, Nanning, China
- 5College of Life Sciences and Technology, Central South University of Forestry and Technology, Changsha, China
- 6Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming, China
Primates have evolved a variety of feeding habits and intestinal physiological structure. Gut microbiome act as metabolic organs in many biological processes and play a vital role in adaptation to dietary niches. Gut microbiome also convert primary bile acids (BAs) to secondary. BAs profile and gut microbiome are together influenced by diets and play a significant role in nutrient absorption. The regulation between gut microbiome and BAs metabolism is bidirectional although the relationship in primates consuming diverse diets is still unclear. Here, we investigated gut microbiome structures, fecal BAs profile, and their relationship in primates preferring three distinct diets. We found that gut microbiome communities are well differentiated among dietary groups. Folivorous primates had higher Firmicutes abundance and lower Prevotella to Bacaeroides ratios, possibly related to fiber consumption. Frugivorous primates are colonized predominantly by Prevotella and Bacteroides, pointing to an increased adaptation to high-sugar and simple carbohydrate diets. Likewise, BA profiles differ according to diet in a manner predictable from the known effects of BAs on metabolism. Folivorous primates have high conjugated bile acid levels and low unconjugated to conjugated BA ratios, consistent with their fiber-rich leaf-eating diet. Much of the differentiation in secondary and unconjugated BAs is associated with microbiome composition shifts and individual bile acid concentrations are correlated with the abundance of distinct bacterial taxonomic groups. Omnivores have higher concentrations of secondary BAs, mainly lithocholic acid (LCA). These levels are significantly positively correlated with the presence of Clostrida species, showing that the digestion requirements of omnivores are different from plant-eating primates. In conclusion, gut microbiome and BAs can respond to changes in diet and are associated with nutrient component consumption in each diet primate group. Our study is the first to demonstrate BA profile differentiation among primates preferring diverse diets. BAs thus appear to work with gut microbiome to help primates adapt to their diet.
Introduction
Over the long course of evolution primates have adapted to diverse diets and display a diverse array of anatomical specializations in their stomachs, caeca, and colons (Lambert, 1998). In addition to these structural organ changes, the set of microorganisms (commonly known as the microbiome) inhabiting gastrointestinal tracts can contribute to metabolic versatility and play an important role in host dietary adaptation, for example, digesting substrates otherwise unavailable to hosts. This can help expand host’s dietary options (Greene et al., 2020). There is mounting evidence that microbial metabolites are the major mediators of the microbiome’s influence on host physiology (Thursby and Juge, 2017; Muller et al., 2021). Bile acids (BAs), a class of metabolites produced in the liver from cholesterol and metabolized in the intestine by the gut microbiome (Wahlström et al., 2016), play an important role in lipid metabolism (De Aguiar Vallim et al., 2013). Studies have shown that the structure and function of the gut microbiome, and BAs composition all affected by diet, they play an important role in the maintenance of intestinal homeostasis and health (Clayton et al., 2018; Katsidzira et al., 2019; Trefflich et al., 2019; Ocvirk and O'Keefe, 2021). Moreover, the gut microbiome/BA profile interplay complicated by microbial biotransformation of BA types on the one hand, and the BA-mediated gut microbiome composition change on the other (Monte et al., 2009; Wahlström et al., 2016; Ocvirk and O'Keefe, 2017; Kriaa et al., 2019).
Host hepatocytes synthesize primary BAs from cholesterol (Winston and Theriot, 2020). These primary products as then conjugated with either glycine or taurine to produce conjugated BAs (Wahlström et al., 2016; Singh et al., 2019). Once host-derived primary BAs enter the gastrointestinal tract, they expedite the digestion and absorption of dietary lipids and lipophilic vitamins by forming micelles (Li and Chiang, 2014; Kiriyama and Nochi, 2019) and regulate glycolipid metabolism by activate receptors expressed in the liver and intestine, also influence immune cell function (De Aguiar Vallim et al., 2013; Fiorucci and Distrutti, 2015; Jia et al., 2019; Kiriyama and Nochi, 2019; Sato et al., 2021; Paik et al., 2022). Gut microbiome influence the composition of host BA pools by chemically modifying them into secondary BAs, mainly by deconjugation, dehydroxylation, dehydrogenation, and epimerization (Spence et al., 2006; Lucas et al., 2021). In turn, BAs also regulate the composition of gut microbiome (Winston and Theriot, 2020). The effects of BAs on gut microbiome include bacterial cell membranes damage (Begley et al., 2005; Kurdi et al., 2006), macromolecule stability disturbance (Begley et al., 2005; Merritt and Donaldson, 2009), and intracellular acidification (Kurdi et al., 2006). These important features of BAs lead to inhibition of the growth of some bacteria in the intestinal tract (Begley et al., 2005; Kurdi et al., 2006; Di Ciaula et al., 2017). Some bile acids have a variety of antibacterial activities (Begley et al., 2005; Kurdi et al., 2006). Thus, the BAs pool is a function of a collaborative metabolism of the host and the gut microbiome and the relationship between the gut microbiome and BAs is bidirectional. For example, changes in gut microbiome composition can alter BA profiles (Chen et al., 2020; Connors et al., 2020; Tanaka et al., 2020; Kamp et al., 2021). A comparative study showed that BA pools of conventional mice are more chemically diverse than gnotobiotic mouse pools (Sayin et al., 2013). BA profiles also vary according to dominant microbiome composition in the first 3 years in Japanese infants (Tanaka et al., 2020). Conversely, insufficient bile secretion, such as biliary obstruction or liver cirrhosis, leads to dramatic gut microbiome proliferation and increased bacterial translocation, highlighting bile acids’ antibacterial activity (Ding et al., 1993; Inagaki et al., 2006).
Given what we know from studies within humans, we would expect that divergence in diets among closely related species is accompanied by changes in gut microbiome and BA composition (Di Ciaula et al., 2017; Singh et al., 2019; Winston and Theriot, 2020). Primate diets include three major groups of nutrients: proteins, lipids (fats), and carbohydrates. These come from a diverse set of foods, including fruit pulp, seeds, leaves (and other structural components of plants), plant exudates, and animal matter (Lambert, 1998). These nutrition sources can be classified into three main diets: folivorous (e.g., Colobidae), omnivorous (e.g., Cercopithecinae), and frugivorous (e.g., Hylobatidae; Struhsaker, 1995). Colobines are a unique group of old world monkeys from southeast Asia and Africa who consume a folivorous diet rich in structural carbohydrates (crude fiber content up to 52%) that is difficult to digest (Lambert, 1998; Nijboer and Clauss, 2006). These primates evolved a foregut fermentation system characterized by a large, sacculated forestomach with a diverse array of microflora to digest the crude fiber (Milton, 1993; Lambert, 1998; Nijboer and Clauss, 2006), chemical secondary compounds (such as tannin, alkaloids, and terpenoids), and toxic secondary metabolites (Milton, 1993; Lambert, 1998). Sichuan snub-nosed monkeys (Rhinopithecus roxellana) are representative colobidae species (Jablonski and Peng, 1993). They have an unusual foregut fermentation system different from other common primates (Milton, 1993; Lambert, 1998). Macaques (Macaca mulatta) are omnivores and representative species of the Cercopithecinae subfamily (Fooden, 2000; Brandon-Jones et al., 2004; Sengupta and Radhakrishna, 2016). They generally feed on twigs, spears, plant seeds, herbs, fruits, leaves, and some bird eggs (Fooden, 2000). The wide food preference of macaques is one of the reasons they are widely distributed and adapted to complex habitats (Hill, 1997). Gibbons (Hylobatidae, except for Hylobates syndactylus), commonly known as lesser apes, are representative primate frugivores that inhabit tropical and subtropical rainforests (Chatterjee, 2009; Hon et al., 2018). Fruit is the primary component of the gibbon diet (40 ~ 78%), with high concentrations of total non-structural carbohydrates (Lambert, 1998; Kim et al., 2011; Ni et al., 2014; Dooley and Judge, 2015; Hon et al., 2018), especially fructose, as well as water and lipids (Lambert, 1998). Gibbons supplement their diets with young leaves, flowers, mature leaves, and insects (Elder, 2009). Such diversification of diet preferences by closely related species makes primates a good model for studying diet adaptation (Milton, 1993).
To elucidate the role that BAs and gut microbiome composition play in evolutionary adaptation to diverse diets, we analyzed fecal BAs and gut microbiome community profiles of representative primate species preferring one of the major diet groups. We aim to uncover correlations between BAs and the gut microbiome that accompany adaptation to specific diets.
Materials and Methods
Sample Collection
We collected 30 fecal samples from three primate species with different dietary preferences: the folivorous Rhinopithecus roxellanae (n = 10, Fol), the omnivorous Macaca mulatta (n = 10, Omn), and the frugivorous Hylobates pileatus, Nomascus annamensis, and Nomascus leucogenys (n = 10, Fru). Additional sample information is in the Supplementary Table S1. Fecal samples of folivorous and omnivorous primates were collected from wild groups, while frugivorous primate fecal samples were collected from captive groups because of difficulties in collecting samples in the field. All fresh fecal samples from primates were collected immediately after defecation and divided into two tubes. One contained RNA-later (Tiangen biotech, Bejing, China) for shotgun sequencing and the other was used to quantify BAs using ultra-high performance liquid chromatography coupled with tandem mass spectrometry (UHPLC–MS/MS). All samples were kept on dry ice and transported to the laboratory while frozen, where they were preserved at −80°C until processing.
Metagenomic Community and Functional Analysis
DNA was extracted using the PowerSoil® DNA Isolation kit (MO BIO Laboratories, Carlsbad, CA, United States). DNA quality and quantity DNA were determined using a nanodrop (ND-1000) spectrophotometer (Nanodrop Technologies, Wilmington, DE, United States) and agarose gel electrophoresis. Shotgun sequencing was performed using Illumina NovaSeq 6000, with at least 10 Gb per sample (Supplementary Table S2). Adaptor and low-quality reads were trimmed using Trimmomatic (version 0.36; Bolger et al., 2014): 3′ tailing sequences were removed when the average quality over a 4-b sliding window was less than 20 and reads shorter than 70 bp were discarded. Host and human contaminations were filtered by aligning to reference genomes (R. roxellana: assembly ASM756505v1; Homo sapiens: assembly GRCh38.p13; M. mulatta: assembly Mmul_10; N. leucogenys assembly Asia_NLE_v1; and Hylobates moloch: assembly HMol_V2) using bowtie2 (v2.3.5; Langmead and Salzberg, 2012) with the option “–very-sensitivelocal.”
Taxonomic profiles were generated with the mOTUs profiler (v2.0.0) 52 (Milanese et al., 2019) using the following parameters: -l 75; -g 2; and -c. Single-sample metagenomic assemblies were produced with MEGAHIT (v1.2.6) (Li et al., 2015) and gene finding was performed on contigs longer than 300 bp using Prodigal. Protein sequences longer than 100 bp were functionally annotated using diamond (v0.9.24; Buchfink et al., 2015) against KEGG (v50) database (Kanehisa et al., 2016) with the parameters: -d -q -e 1e-5 -k 1.
MOTUs profiles and K-numbers were first converted to relative abundance to account for library size, then were analyzed using a nonparametric test and correlation (Supplementary Tables S4 S5). Species abundance was used to calculate alpha and beta diversities. Alpha diversity analysis was performed using observed species (Ace and Shannon index), implemented in R (vegan). Venn and principal component analyses (PCA) were performed using the OmicStudio tools.1 Linear discriminant analysis (LDA) effect size (LEfSe) was estimated using Galaxy.2
Metabolomic Profiling
A UHPLC–MS/MS system was used to quantify BAs at Novogene Co., Ltd (Beijing, China; ExionLC™ AD UHPLC-QTRAP 6500+, AB SCIEX Corp., Boston, MA, United States). We analyzed 33 kinds of BAs along with six stable isotope-labeled internal standards obtained from ZZ Standards Co., LTD (Shanghai, China). Ammonium acetate was of analytical grade and obtained from Sigma-Aldrich (St. Louis, MO, United States). Methanol (Optima LC–MS), acetonitrile (Optima LC–MS), and formic acid (Optima LC–MS) were purchased from Thermo-Fisher Scientific (FairLawn, NJ, United States). Standard operating procedures were followed to prepare standard solutions and extract BAs from fecal samples.
Separation was performed on an Agela Venusil MP C18 column (2.1 × 100 mm, 2.5 μm) maintained at 50°C. The mobile phase, consisting of 0.1% formic acid in water (solvent A) and acetonitrile (solvent B), was delivered at a flow rate of 0.50 ml/min. The solvent gradient was set as follows: initial 20% B, 0.5 min; 20–35% B, 1 min; 35–37% B, 2.5 min; 37–38% B, 4.1 min; 38–39% B, 6 min; 39–40% B, 6.5 min; 40–44% B, 8.5 min; 44–45% B, 9 min; 45–52% B, 9.5 min; 52–55% B, 11.5 min; 55–100% B, 12.5 min; 100–20% B, 15.1 min; and 20% B, 17.1 min. The mass spectrometer was operated in negative multiple reaction mode (MRM). Parameters were as follows: IonSpray Voltage (−4,500 V), Curtain Gas (30 psi), Ion Source Temp (550°C), and Ion Source Gas of 1 and 2 (65 psi; Supplementary Table S3).
Quality control (QC) was used to evaluate the stability of methods and machines (Supplementary Table S6–1). Data files generated by HPLC-MS/MS were processed using the SCIEX OS Version 1.4 to integrate and correct peaks (Supplementary Table S6–2). Data pre-processing was done according to the “80% rule” (Bijlsma et al., 2006) and the remaining missing values were imputed using mean values (Supplementary Table S6–3). The processed data set was then entered into the SIMCA-P software package (v13.0) and was used to perform orthogonal to partial least squares-discriminate analysis (OPLS-DA).
Statistical Analysis
Statistical analyses were performed using GraphPad Prism version 5.0 (GraphPad Software). Differences in characteristics between groups were analyzed using the Kruskal–Wallis test with Dunn post hoc tests (pFO, Folivorous-Omnivorous; pFF, Folivorous-Frugivorous; and pOF, Omnivorous-Frugivorous) used the adjusted p-value. Spearman correlations were used to explore the relationship between fecal concentrations of individual BAs metabolites and the relative abundance of bacterial taxonomic groups obtained by the LEfSe analysis, performed using OmicStudio tools. Value of p < 0.05 was considered statistically significant.
Results
Among Dietary Group Gut Microbiome Differentiation
We identified 473 species from fecal samples of three dietary groups of primates. These species belong to 14 phyla, 24 classes, 38 orders, 59 families, and 110 genera (Supplementary Table S4). We found 192 species of bacteria in the folivorous (Fol) group, 237 in the omnivorous (Omn) group, and 336 in frugivorous (Fru) group (Figure 1A). Interestingly, only 91 species of bacteria are shared among all three groups, indicating appreciable species composition differentiation accompanying shifts in dietary preference (Figure 1A). Moreover, alpha diversity estimates identify significant differences among the three groups (Shannon index: pFF = 0.0034, pOF = 0.0412; ACE index: pFF < 0.0001, pFO = 0.0045; Figure 1B). Beta diversity analysis likewise shows tight, but well-separated, within-group clusters of individuals (R = 0.8333, p = 0.001; Figure 1C), reinforcing the conclusion that microbiome compositions have significantly diverged according to primate dietary preference.
Figure 1. Gut microbiome structure in three primate groups. (A) Venn diagram at the species level. (B) Box-plots of alpha diversity (ACE and Shannon index) at the species level. (C) Microbiome beta diversity (Principal component analysis, PCA) plot among primates at the species level. (D) Cladogram of differential microbial prevalence identified by Linear discriminant analysis (LDA), effect size (LEfSe) analysis (LDA score > 4, one against all), from the phylum (innermost ring) to genus (outermost ring). Differentiated bacterial taxa are marked by lowercase letters. Each small circle at different taxonomic levels represents a taxon at that level, and the diameter of the circle is proportional to relative abundance. Taxa with no significant difference are yellow. Differentiated taxa are colored according to the group with the highest species abundance. Colors represent groups and nodes are colored according to the communities that play an important role in the group represented by the color. A value of p < 0.05 was considered statistically significant. *p < 0.05; **p < 0.01.
We further used LEfSe to draw a cladogram that reflects differential occurrence of taxonomic terms in individual samples (Supplementary Figure S1; Figure 1D). Nine taxa are significantly more overrepresented among folivores, which belong to Firmicutes, Verrucomicrobia, Proteobacteria, and Euryarchaeota (Figure 1D; Supplementary Figure S1). Clostridia and Clostridiales are among the nine taxa preferentially found in omnivorous primates (Figure 1D; Supplementary Figure S1), while 21 taxonomic terms overabundant in frugivores belong to Bacteroidetes, Chloroflexi, Firmicutes, and Proteobacteria (Figure 1D; Supplementary Figure S1). Folivores have the highest log F/B (pFF < 0.0001, pOF = 0.0049) and the lowest log P/B (pFO = 0.0008, pFF = 0.0177; Supplementary Figures S2, S3). In contrast, omnivores and frugivores have substantially higher log P/B (pFO = 0.0008 and pFF = 0.0177; Supplementary Figure S3). These results show that gut microbiome communities are influenced by host diets and may play an important role in host adaptation.
Among Dietary Group BAs Composition Differentiation
We next compared BAs profiles among groups of non-human primates consuming different diets. We observed obvious differences in BAs composition among diet groups, when viewed directly (Figure 2A) or through an OPLS-DA analysis (Figure 2B). While total BAs amounts do not appear to differ among groups (Figure 2C), we noticed differences in BAs composition (Figures 2D–H). There are more conjugated BAs (pFO = 0.0034; Figure 2F) and a lower unconjugated to conjugated BAs ratio (pFO = 0.0034; Figure 2H) in folivores than in other groups, while omnivores produce more total secondary BAs (pFO = 0.0169).
Figure 2. BAs metabolome. (A) BAs composition column. (CA, cholic acid; GCA, glycocholic acid hydrate; TCA, taurocholic acid; CDCA, chenodeoxycholic acid; GCDCA, glycochenodeoxycholic acid sodium salt; CDCA_3Gln, chenodeoxycholic acid-3-β-D-glucuronide; TCDCA, taurochenodeoxycholic acid sodium salt; DCA, deoxycholic acid; GDCA, glycodeoxycholic acid; TDCA, taurodeoxycholic acid sodium salt; UDCA, ursodeoxycholic acid; GUDCA, glycoursodeoxycholic acid; TUDCA, tauroursodeoxycholic acid dihydrate; βUDCA, 3β-ursodeoxycholic acid; LCA, lithocholic acid; TLCA, taurolithocholic acid sodium salt; GLCA, glycolithocholic acid; LCA_S, lithocholic acid 3-sulfate sodium salt; HCA, hyocholic acid; HDCA, hyodeoxycholic acid; 7_ketoLCA, 7-ketolithocholic acid; NorCA, nor cholic acid; 23norDCA, 23-nordeoxycholic acid; isoLCA, isolithocholic acid; 12_ketoLCA, 12-ketolithocholic acid; AlloLCA, allolithocholic acid; 3_DHCA, 3-dehydrocholic acid; α-MCA, alpha-muricholic acid; and β-MCA, beta-muricholic acid). (B) Orthogonal to partial least squares-discriminate analysis (OPLS-DA) score plots of fecal samples from the three primate dietary groups. (C) Total bile acids (BAs) levels. (D) Total primary BAs (%). (E) Total secondary BAs (%). (F) Total conjugated BAs (%). (G) Total unconjugated BAs (%). (H) The ratio of conjugated/unconjugated BAs. A value of p < 0.05 was considered statistically significant. *p < 0.05; **p < 0.01. All data are shown in box-and-whisker plots according to Tukey.
We used the Kruskal–Wallis test with Dunn post-hoc tests to identify individual BAs whose levels differ among groups. We found that 15 such compounds (Figure 3). GCA (pFO = 0.0026, pFF = 0.0105), TCA (pFF = 0.0168), GDCA (pFO = 0.0006, pFF = 0.0014), TDCA (pFO = 0.0287, pFF = 0.0010), HCA (pFF = 0.0228), and NorCA (pFO = 0.0003, pFF = 0.0068) are highest in folivores, while individuals in this group are relatively deficient in alloLCA (pFO = 0.0133, pFF = 0.0015; Figure 3). Omnivores produce an abundance of LCA (pFO = 0.0287, pOF = 0.0041), isoLCA (pOF = 0.0063), and 3_DHCA (pFO = 0.0002), but relatively little CDCA (pOF = 0.0002), β-MCA (pFO = 0.0248), and DCA (pFO = 0.0475, pOF = 0.0412; Figure 3). Finally, frugivores produce a relatively high levels of CDCA (pOF = 0.0002) and TUDCA (pOF < 0.0003), but relatively little TLCA (pFF = 0.0111; Figure 3). It thus appears that BAs composition does vary with dietary preference.
Figure 3. Individual BAs amounts. BAs in feces were analyzed by the Kruskal–Wallis test with Dunn post hoc tests (for nonparametric unpaired data) using GraphPad Prism version 5.0. A value of p < 0.05 was considered statistically significant. *p < 0.05; **p < 0.01. All data are shown in box-and-whisker plots according to Tukey.
Individual BA Concentrations Correlate With the Gut Microbiome
To investigate whether microbiome variation is associated with BAs composition, we estimated Spearman correlations within diet groups. We found 306 significant (p < 0.05) associations. HCA, NorCA, and most conjugated BAs (primary BAs: GCA and TCA; secondary BAs: GDCA and TDCA) are positively correlated with Enterobacterales (Proteobacteria) and Coraliomargarita (Verrucomicrobia) in folivores (Figure 4A). In contrast, alloLCA and 3_DHCA abundances are negatively correlated with the prevalence of these bacteria (Figure 4A). LCA is positively correlated with Clostridia and Helicobacter (Campylobacterales) in omnivores, while CDCA and TUDCA are negatively correlated with Intestinibacter (Clostridia) and Helicobacter (Campylobacterales; Figure 4A). Most conjugated BAs are negatively correlated with Chloroflexi and Prevotella in frugivores, whereas CDCA and TUDCA are positively correlated with Phascolarctobacterium abundance (Figure 4A).
Figure 4. Association between BAs composition and gut microbial diversity. (A) Heat map summarizing the correlations of fecal bile acid concentrations and differentially abundant taxa relative abundance between groups (*p < 0.05; **p < 0.01). (B) Secondary BAs metabolic pathways modified by the microbiome. (a) BSH deconjugation of T/GCA and T/GCDCA. (b) 7α/β-dehydroxylation of CA and CDCA. (c) 7α/β-dehydrogenation of CA and DCA. (d) 3α/β-dehydrogenation of CA and CDCA. The hydrolysis of taurine or glycine conjugated BAs to free BAs is performed by the bacterial bile salt hydrolase (BSH; K01442; Di Ciaula et al., 2017). In addition, following deconjugation, gut microbiota can perform 7-dehydroxylation, involving a multistep biochemical pathway found only in anaerobic gut bacteria (7-dehydroxylation: K15868, K15869, K15870, K15871, K15872, K15873, K15874, and K07007; Monte et al., 2009), bacterial dehydratases of the anaerobic flora from this region attack and remove the hydroxyl group to form 7-deoxy BAs (Monte et al., 2009), including 7α-dehydroxylation of CA and CDCA yielding DCA and LCA, respectively; BAs 7β-dehydroxylation of UDCA yielding LCA (Di Ciaula et al., 2017). Another well-recognized transformation is carried out by hydroxysteroid dehydrogenases (HSD) of intestinal bacteria (Monte et al., 2009; Ridlon et al., 2016; Lucas et al., 2021). HSDs catalyze reversible oxidation of hydroxyl groups on the C-3 (3α-HSD: K22604, K22605; 3β-HSD: K22606, K22607) and C-7 (7α-HSD: K00076; 7β-HSD: K23231) carbon positions of the bile acid steroid core (Monte et al., 2009; Ridlon et al., 2016; Lucas et al., 2021).
To see if there is any evidence that correlations between BAs and bacterial taxon abundance may be driven by microbiome metabolism, we focused on genes in our sequence pool that may be involved in BAs metabolism. We found many genes whose products are predicted to be involved in BSH deconjugation (K01442) and 7-dehydroxylation (K15868, K15869, K15870, K15871, K15872, K15873, K15874, and K07007). We also identified genes coding for the HSD class of enzymes that catalyze the reversible oxidation of hydroxyl groups on the C-3 (3α-HSD: K22604 and K22605; 3β-HSD: K22606 and K22607) and C-7 (7α-HSD: K00076; 7β-HSD: K23231) carbon positions of the BAs steroid core. Furthermore, we found genes whose products belong to 13 KEGG ontologies involved in secondary BAs synthesis. Products of these genes cover all the major metabolic steps, except for K15869 and K23231. These were not found possibly due to detection and processing problems. Of the genes, we did identify, only three (K01442, K22605, and K07007) were detected across all samples, and one (K22607) was absent in omnivores (Supplementary Table S5). The BSH deconjugation process is both a prerequisite and a limitation for other transformation processes. There was no difference in the prevalence of the gene (K01442) that codes for BSH among dietary groups. However, we found more genes coding for 7-dehydroxylation enzymes among folivores (Figure 4B). This increase in enzyme-coding gene prevalence is broad, including K15868 (pFO = 0.0039), K15870 (pFO = 0.0183), K15872 (pFO = 0.0065), and K07007 (pFF = 0.0002, pOF = 0.0309). K15871, K15873, and K15874 are no significantly different among dietary groups (Figure 4B). Furthermore, genes coding for 7α/β-HSD and 3α/β-HSD processing enzymes, K00076 (pFF = 0.0021), K22604 (pFO = 0.0436), and K22607 (pFO = 0.0203) are also overrepresented in folivores (Figure 4B). In contrast, 7α/β-HSD coding gene K00076 (pFF = 0.0021) is underrepresented among frugivores (Figure 4B).
Discussion
Using microbiome shotgun sequencing and metabolite mass spectrometry, we examined gut microbiome and BAs metabolite profiles of non-human primates consuming a variety of diets. We found significant differences in both gut microbiome composition and BAs profiles among three major dietary groups (folivores, omnivores, and frugivores). Furthermore, correlations between bacterial taxon abundance and metabolite concentration vary among three groups. We further focused on the abundance of individual bacterial genes coding for BA metabolizing enzymes and found that some were associated bile acid profile differences among dietary groups. This suggests that the microbiome plays a direct role in diet-driven differentiation in BAs prevalence among primate species.
Gut Microbiome Differentiation Among Primate Species Grouped by Dietary Preference
We found clear separation in microbial profiles among the three major dietary groups of primates, with very little overlap. This observation is consistent with previous studies that found correlations of structure and function of intestinal microbiome phylogenies and diets (Clayton et al., 2018). Firmicutes are overrepresented in folivores, along with a low Prevotella to Bacaeroides ratio. This combination is thought to promote dietary fiber digestion (Figure 1D; Supplementary Figures S1, S3; Chen et al., 2017; Nobel et al., 2018), which is important for this group since a folivorous diet includes up to 52% crude fiber (Lambert, 1998; Nijboer and Clauss, 2006). Omnivores harbor relatively more Clostrida, bacteria involved in lipid metabolism (Figure 1D; Supplementary Figure S1; Petersen et al., 2019). Both omni- and frugivores have high Prevotella to Bacaeroides ratios, suggesting a decrease in fiber digestion capacity (Supplementary Figure S3; Chen et al., 2017). High Firmicutes to Bacteroidetes ratio promote efficient energy harvesting (Supplementary Figure S2; Turnbaugh et al., 2006). It is thus unsurprising to find that omnivores and frugivores exhibit this trend. In particular, frugivores consume a diet high in non-structural carbohydrates (Lambert, 1998; Kim et al., 2011; Ni et al., 2014; Dooley and Judge, 2015; Hon et al., 2018). Consistent with this, we found that Prevotella (Bacteroidetes) and Bacteroides (Bacteroidetes) are enriched in frugivores. These bacteria help to degrade simple sugars and carbohydrates (Figure 1D; Supplementary Figure S1; Hale et al., 2019).
Primate Diets Are Associated With Distinct Bile Microbial Transformation Profiles
We found that folivorous primates have high levels of conjugated BAs (Figure 2H). Furthermore, abundance of these metabolites (primary BAs: GCA and TCA; secondary BAs: GDCA and TDCA) were positively correlated with the prevalence of Enterobacterales (Proteobacteria) and Puniceicoccaceae (Verrucomicrobia; Figure 4A). Enterobacterales (Proteobacteria) and Coraliomargarita (Verrucomicrobia) are probably crucial intestinal flora influencing levels on the metabolic health biomarker HCA (Zheng et al., 2021). Narrow diet preference often leads to low intestinal diversity in folivorous animals. Indeed, we found that folivores in our sample exhibit lower α-diversity, possibly resulting in fewer taxa contributing bsh genes (Connors et al., 2020). However, we observed no difference in bsh gene (K01442) abundance among dietary groups (Figure 4B), showing that BAs deconjugation is not deficient in folivores despite low taxonomic diversity.
We found a positive association between LCA and Clostridia in omnivores (Figure 4A). LCA is transformed from CDCA by 7α-dehydroxylation, a multistep pathway that is highly conserved and has primarily been observed in Clostridia species (Vital et al., 2019). In addition, although increased fat intake leads to increases in LCA (Zeng et al., 2019), Clostridia can block lipid absorption and thus regulate lipid metabolism (Petersen et al., 2019). Recent studies have observed the loss of Clostridia colonization and function in individuals with metabolic syndrome and obesity. Amounts of most conjugated BAs are negatively correlated with Chloroflexi and Prevotella prevalence in frugivores (Figure 4A). While high-sugar diets cause metabolic disorders (Qiu et al., 2021), aggravate chronic colitis, and increase the incidence of colitis related tumors (51), Prevotella (Bacteroidetes) can help degrade simple sugars and carbohydrates, alleviating the adverse effects of fruit-derived sugars (Lambert, 1998; Kim et al., 2011; Ni et al., 2014; Dooley and Judge, 2015; Hon et al., 2018; Hale et al., 2019). CDCA, a metabolite associated with glucose and lipid metabolism regulation (Yu et al., 2015), is positively correlated with Phascolarctobacterium (Negativicutes) in frugivores. CDCA is beneficial to the host because it can prevent the growth of Clostridioides difficile (Nagao-Kitamoto et al., 2020), reduce inflammation, and improve gastrointestinal function (Nagao-Kitamoto et al., 2020).
Gut microbiome can carry out numerous biotransformations of BAs during their enterohepatic circulation (Ridlon et al., 2016; Singh et al., 2019). BAs metabolism by the gut microbiome mainly involve several categories of reactions, including deconjugation, dehydroxylation, dehydrogenation, epimerization, and oxidation (Monte et al., 2009; Nobel et al., 2018; Lucas et al., 2021). On the other hand, BAs have antimicrobial properties and can influence the species composition of gut microbiome (Di Ciaula et al., 2017). We found that individual bile acid and microbiome taxa are correlated, likely together playing a role in promoting nutrient absorption and metabolism, as well as maintenance of intestinal health under diverse diets.
BA Profiles Differ According to Primate Feeding Habits
Primates have adapted to a wide variety of diets, with concomitant physiological, metabolic, and behavioral strategies to cope with difficulties presented by each feeding habit (Milton, 1993). Morphological differentiation of the digestive tract is one of the prominent adaptations driven by dietary shifts. Cercopithecine monkeys are distinguished from other primates by the presence of cheek pouches where they temporarily store and begin to digest food (Murray, 2011). Amylase found in the saliva of cercopithecine cheek pouches appears to be more active in starch digestion than the human homolog (Lambert, 1998). Colobines evolved a large sacculated forestomach with a diverse array of microbiome to digest the crude fiber in their folivorous diet (Milton, 1993; Struhsaker, 1995; Lambert, 1998). Other common primates are usually hindgut fermenters with an enlarged colon or cecum (Milton, 1993; Lambert, 1998). These two basic evolutionary strategies, morphological and behavioral, to cope with foraging difficulties are not mutually exclusive. Species vary in the extent to which they employ one or the other (Milton, 1993). For example, frugivorous and folivorous primates also supplement their diet with insects and other animal foods (Garber, 1987). This supplemental diet provides them with some of the nutrients they lack and is part of their dietary strategy (Milton, 1993).
Most previous studies have directly studied the relationship between dietary habits and gut microbiome, but rarely considered the existence of other mediating factors. We found significant differences in gut microbiome composition accompanied by shifts in BAs metabolite profiles. Folivores have more conjugated BAs (Figure 2F). These may be preferentially excreted in feces because of binding to fiber molecules before being deconjugated (Singh et al., 2019). This binding may prevent BAs re-entry into enterohepatic circulation and lead to excretion in feces (Naumann et al., 2019), thereby reducing secondary BAs production. The higher concentration of HCA in folivores indicates that these individuals are at a lower risk of metabolic problems. HCA species play critical roles in glucose homeostasis and are good predictive biomarkers of metabolic health (Zheng et al., 2021). LCA, a secondary BAs metabolite, is elevated in omnivore primates (Figure 3). LCA is one of the most hydrophobic BAs and is one of the most toxic BAs (Zheng et al., 2021). LCA accumulation is associated with increased intake of fat (Zheng et al., 2021). CDCA, a primary BAs, is prevalent in frugivorous primates (Figure 3). CDCA may be closely associated with glucose homeostasis and the metabolic status (Yu et al., 2015). Increase in CDCA levels may lead to stronger activity in regulating glucolipids and energy metabolism. Thus, CDCA elevation may be a mechanism of self-protection by frugivores in the face of a high-sugar diet (Yu et al., 2015). Gut microbiome and BAs can respond to changes in diet and are associated with fiber, lipid, sugar, and carbohydrate consumption in each diet primate group. Similar patterns have been observed in other systems. Remission induced by a therapeutic hydrolyzed protein diet is linked to improved microbiome structure in canine chronic inflammatory enteropathy, marked by decreased relative abundance of pathobionts (e.g., Escherichia coli and Clostridium perfringens), increased abundance of a secondary BAs producer (Clostridium hiranonis), and increased levels of secondary BAs, LCA, and DCA (Wang et al., 2019). Physiological levels of these BAs inhibit in vitro growth of disease-associated taxa (Wang et al., 2019). High-fat diet results in BAs profile alterations, particularly in elevated unconjugated and secondary BAs. In association with gut microbiome taxa, this elevation likely confers unfavorable impacts on colonic and host cardiometabolic health in healthy young adults (Wan et al., 2020). In addition, long-term diet influences intestinal health and correlation between BAs levels and microbiome composition may be a maintenance and regulatory mechanism of intestinal health under diverse diets (Ocvirk and O'Keefe, 2017; Oluwagbemigun et al., 2021).
The current study has several potential limitations. Stool samples were collected from different primate families, while fecal samples from frugivores were obtained from captive gibbon individuals. This forced difference in sampling may have confounding effects on our conclusions. Further experimental validation of the results in other primate populations is necessary. Our study is the first to demonstrate differences in BA profiles in primates with diverse diets. We also examine correlations between BAs and gut microbiome, highlighting a promising avenue for primate diet adaption.
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 at: NCBI, BioProject ID PRJNA799478 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA799478), and Supplementary Material.
Ethics Statement
The animal study was reviewed and approved by the Committee for Animal Experiments of the Institute of Zoology (CAS).
Author Contributions
ML and HP designed the study. XL, MZ, and SW collected samples in the field. XL and XW performed the experiments. XL, XW, and ZW analyzed the data. XL, XW, ZW, HP, and ML wrote the manuscript. ZX provided constructive comments. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Natural Science Foundation of China and the Strategic Priority Research Program of Chinese Academy of Sciences (31821001, XDB31000000, and 32070404).
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.
Acknowledgments
The authors thank Yao Hui for help in sample collecting in the field. We also thank Li Guanghao and Yuan Shenglei for their constructive suggestion of this paper.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.899102/full#supplementary-material
Supplementary Table S1 | Samples information.
Supplementary Table S2 | Sample sequencing data evaluation.
Supplementary Table S3 | Parameters of ultra-high performance liquid chromatography coupled to tandem mass spectrometry (UHPLC–MS/MS).
Supplementary Table S4 | Taxonomy annotation.
Supplementary Table S5 | KEGG annotation.
Supplementary Table S6 | Data of ultra-high performance liquid chromatography coupled to tandem mass spectrometry (UHPLC-MS/MS). S6-1: Quality control data. S6-2: BAs quantitative results of samples. S6-3: Results of BAs data processing.
Supplementary Figure S1 | Microbiome taxa with linear discriminant analysis (LDA) score greater than four of different primate groups.
Supplementary Figures S2 and S3 | S2: Log Firmicutes/Bacteroidetes (Log F/B) and S3: Log Prevotella/Bacteroides (Log P/B) were analyzed by analyzed by the Kruskal–Wallis test with Dunn post hoc tests (for nonparametric unpaired data) using GraphPad Prism version 5.0.
Footnotes
References
Begley, M., Gahan, C. G., and Hill, C. (2005). The interaction between bacteria and bile. FEMS Microbiol. Rev. 29, 625–651. doi: 10.1016/j.femsre.2004.09.003
Bijlsma, S., Bobeldijk, I., Verheij, E. R., Ramaker, R., Kochhar, S., Macdonald, I. A., et al. (2006). Large-scale human metabolomics studies: a strategy for data (pre-) processing and validation. Anal. Chem. 78, 567–574. doi: 10.1021/ac051495j
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170
Brandon-Jones, D., Eudey, A. A., Geissmann, T., Groves, C. P., Melnick, D. J., Morales, J. C., et al. (2004). Asian primate classification. Int. J. Primatol. 25, 97–164. doi: 10.1023/B:IJOP.0000014647.18720.32
Buchfink, B., Xie, C., and Huson, D. H. (2015). Fast and sensitive protein alignment using diamond. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176
Chatterjee, H. J. (2009). “Evolutionary relationships among the gibbons: a biogeographic perspective,” in The Gibbons. Developments in Primatology: Progress and Prospects. eds. D. Whittaker and S. Lappan (New York, NY: Springer), 13–36.
Chen, T., Long, W., Zhang, C., Liu, S., Zhao, L., and Hamaker, B. R. (2017). Fiber-utilizing capacity varies in Prevotella-versus Bacteroides-dominated gut microbiota. Sci. Rep. 7:2594. doi: 10.1038/s41598-017-02995-4
Chen, W., Wei, Y., Xiong, A., Li, Y., Guan, H., Wang, Q., et al. (2020). Comprehensive analysis of serum and fecal BAs profiles and interaction with gut microbiota in primary biliary cholangitis. Clin. Rev. Allergy Immunol. 58, 25–38. doi: 10.1007/s12016-019-08731-2
Clayton, J. B., Gomez, A., Amato, K., Knights, D., Travis, D. A., Blekhman, R., et al. (2018). The gut microbiome of nonhuman primates: lessons in ecology and evolution. Am. J. Primatol. 80:e22867. doi: 10.1002/ajp.22867
Connors, J., Dunn, K. A., Allott, J., Bandsma, R., Rashid, M., Otley, A. R., et al. (2020). The relationship between fecal BAs and microbiome community structure in pediatric Crohn's disease. ISME J. 14, 702–713. doi: 10.1038/s41396-019-0560-3
De Aguiar Vallim, T. Q., Tarling, E. J., and Edwards, P. A. (2013). Pleiotropic roles of BAs in metabolism. Cell Metab. 17, 657–669. doi: 10.1016/j.cmet.2013.03.013
Di Ciaula, A., Garruti, G., Baccetto, R. L., Molina-Molina, E., Bonfrate, L., Wang, D. Q., et al. (2017). BAs physiology. Ann. Hepatol. 16, S4–S14. doi: 10.5604/01.3001.0010.5493
Ding, J. W., Andersson, R., Soltesz, V., Willén, R., and Bengmark, S. (1993). The role of bile and BAs in bacterial translocation in obstructive jaundice in rats. Eur. Surg. Res. 25, 11–19. doi: 10.1159/000129252
Dooley, H. M., and Judge, D. S. (2015). Kloss gibbon (Hylobates klossii) behavior facilitates the avoidance of human predation in the Peleonan forest, Siberut Island, Indonesia. Am. J. Primatol. 77, 296–308. doi: 10.1002/ajp.22345
Elder, A. A. (2009). “Hylobatid diets revisited: the importance of body mass, fruit availability, and interspecific competition,” in The Gibbons. Developments in Primatology: Progress and Prospects. eds. D. Whittaker and S. Lappan (New York, NY: Springer Press), 133–159.
Fiorucci, S., and Distrutti, E. (2015). BAs-activated receptors, intestinal microbiota, and the treatment of metabolic disorders. Trends Mol. Med. 21, 702–714. doi: 10.1016/j.molmed.2015.09.001
Fooden, J. (2000). Systematic review of the rhesus macaque, Macaca mulatta (Zimmermann, 1780). Field. Zool. 96, 1–180.
Garber, P. A. (1987). Foraging strategies among living primates. Annu. Rev. Anthropol. 16, 339–364. doi: 10.2307/2155875
Greene, L. K., Williams, C. V., Junge, R. E., Mahefarisoa, K. L., Rajaonarivelo, T., Rakotondrainibe, H., et al. (2020). A role for gut microbiota in host niche differentiation. ISME 14, 1675–1687. doi: 10.1038/s41396-020-0640-4
Hale, V. L., Tan, C. L., Niu, K., Yang, Y., Zhang, Q., Knight, R., et al. (2019). Gut microbiota in wild and captive Guizhou snub-nosed monkeys, Rhinopithecus brelichi. Am. J. Primatol. 81:e22989. doi: 10.1002/ajp.22989
Hill, D. A. (1997). Seasonal variation in the feeding behavior and diet of Japanese macaques (Macaca fuscata yakui) in lowland forest of Yakushima. Am. J. Primatol. 43, 305–320. doi: 10.1002/(SICI)1098-2345(1997)43:4<305::AID-AJP2>3.0.CO;2-0
Hon, N., Behie, A. M., Rothman, J. M., and Ryan, K. G. (2018). Nutritional composition of the diet of the northern yellow-cheeked crested gibbon (Nomascus annamensis) in northeastern Cambodia. Primates 59, 339–346. doi: 10.1007/s10329-018-0663-x
Inagaki, T., Moschetta, A., Lee, Y. K., Peng, L., Zhao, G., Downes, M., et al. (2006). Regulation of antibacterial defense in the small intestine by the nuclear BAs receptor. Proc. Natl. Acad. Sci. U. S. A. 103, 3920–3925. doi: 10.1073/pnas.0509592103
Jablonski, N. G., and Peng, Y. Z. (1993). The phylogenetic relationships and classification of the doucs and snub-nosed langurs of China and Vietnam. Folia Primatol. 60, 36–55. doi: 10.1159/000156674
Jia, E. T., Liu, Z. Y., Pan, M., Lu, J. F., and Ge, Q. Y. (2019). Regulation of BAs metabolism-related signaling pathways by gut microbiota in diseases. J Zhejiang Univ Sci B 20, 781–792. doi: 10.1631/jzus.B1900073
Kamp, K. J., Cain, K. C., Utleg, A., Burr, R. L., Raftery, D., Luna, R. A., et al. (2021). BAs and microbiome among individuals with irritable bowel syndrome and healthy volunteers. Biol. Res. Nurs. 23, 65–74. doi: 10.1177/1099800420941255
Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. (2016). KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 44, D457–D462. doi: 10.1093/nar/gkv1070
Katsidzira, L., Ocvirk, S., Wilson, A., Li, J., Mahachi, C. B., Soni, D., et al. (2019). Differences in fecal gut microbiota, short-chain fatty acids and BAs link colorectal cancer risk to dietary changes associated with urbanization among Zimbabweans. Nutr. Cancer 71, 1313–1324. doi: 10.1080/01635581.2019.1602659
Kim, S., Lappan, S., and Choe, J. C. (2011). Diet and ranging behavior of the endangered Javan gibbon (Hylobates moloch) in a submontane tropical rainforest. Am. J. Primatol. 73, 270–280. doi: 10.1002/ajp.20893
Kiriyama, Y., and Nochi, H. (2019). The biosynthesis, signaling, and neurological functions of BAs. Biomolecules 9:232. doi: 10.3390/biom9060232
Kriaa, A., Bourgin, M., Potiron, A., Mkaouar, H., Jablaoui, A., Gérard, P., et al. (2019). Microbial impact on cholesterol and BAs metabolism: current status and future prospects. J. Lipid Res. 60, 323–332. doi: 10.1194/jlr.R088989
Kurdi, P., Kawanishi, K., Mizutani, K., and Yokota, A. (2006). Mechanism of growth inhibition by free BAs in lactobacilli and bifidobacteria. J. Bacteriol. 188, 1979–1986. doi: 10.1128/JB.188.5.1979-1986.2006
Lambert, J. E. (1998). Primate digestion: interactions among anatomy, physiology, and feeding ecology. Evol. Anthropol. Issues News Rev. 7, 8–20. doi: 10.1002/(SICI)1520-6505(1998)7:13.0.CO;2-C
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923
Li, T., and Chiang, J. Y. (2014). BAs signaling in metabolic disease and drug therapy. Pharmacol. Rev. 66, 948–983. doi: 10.1124/pr.113.008201
Li, D., Liu, C. M., Luo, R., Sadakane, K., and Lam, T. W. (2015). MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31, 1674–1676. doi: 10.1093/bioinformatics/btv033
Lucas, L. N., Barrett, K., Kerby, R. L., Zhang, Q., Cattaneo, L. E., Stevenson, D., et al. (2021). Dominant bacterial phyla from the human gut show widespread ability to transform and conjugate BAs. mSystems 31:e0080521. doi: 10.1128/mSystems.00805-21
Merritt, M. E., and Donaldson, J. R. (2009). Effect of bile salts on the DNA and membrane integrity of enteric bacteria. J. Med. Microbiol. 58, 1533–1541. doi: 10.1099/jmm.0.014092-0
Milanese, A., Mende, D. R., Paoli, L., Salazar, G., Ruscheweyh, H. J., Cuenca, M., et al. (2019). Microbial abundance, activity and population genomic profiling with mOTUs2. Nat. Commun. 10:1014. doi: 10.1038/s41467-019-08844-4
Milton, K. (1993). Diet and primate evolution. Sci. Am. 269, 86–93. doi: 10.1038/scientificamerican0893-86
Monte, M. J., Marin, J. J., Antelo, A., and Vazquez-Tato, J. (2009). BAs: chemistry, physiology, and pathophysiology. World J. Gastroenterol. 15, 804–816. doi: 10.3748/wjg.15.804
Muller, E., Algavi, Y. M., and Borenstein, E. (2021). A meta-analysis study of the robustness and universality of gut microbiome-metabolome associations. Microbiome 9:203. doi: 10.1186/s40168-021-01149-z
Murray, P. (2011). “The role of cheek pouches in cercopithecine monkey adaptive strategy,” in Primate Functional Morphology and Evolution (Berlin, New York: De Gruyter Mouton), 151–194.
Nagao-Kitamoto, H., Leslie, J. L., Kitamoto, S., Jin, C., Thomsson, K. A., Gillilland, M. G., et al. (2020). Interleukin-22-mediated host glycosylation prevents Clostridioides difficile infection by modulating the metabolic activity of the gut microbiota. Nat. Med. 26, 608–617. doi: 10.1038/s41591-020-0764-0
Naumann, S., Schweiggert-Weisz, U., Eglmeier, J., Haller, D., and Eisner, P. (2019). In vitro interactions of dietary fibre enriched food ingredients with primary and secondary BAs. Nutrients 11:1424. doi: 10.3390/nu11061424
Ni, Q. Y., Huang, B., Liang, Z. L., Wang, X. W., and Jiang, X. L. (2014). Dietary variability in the western black crested gibbon (Nomascus concolor) inhabiting an isolated and disturbed forest fragment in southern Yunnan, China. Am. J. Primatol. 76, 217–229. doi: 10.1002/ajp.22224
Nijboer, J., and Clauss, M. (2006). “The digestive physiology of colobine primates,” in Fibre Intake and Faeces Quality in Leaf-Eating Monkeys. ed. J. Nijboer (Ridder: Utrecht University Repository), 9–28.
Nobel, Y. R., Snider, E. J., Compres, G., Freedberg, D. E., Khiabanian, H., Lightdale, C. J., et al. (2018). Increasing dietary fiber intake is associated with a distinct esophageal microbiome. Clin. Transl. Gastroenterol. 9:199. doi: 10.1038/s41424-018-0067-7
Ocvirk, S., and O'Keefe, S. J. (2017). Influence of BAs on colorectal cancer risk: potential mechanisms mediated by diet-gut microbiota interactions. Curr. Nutr. Rep. 6, 315–322. doi: 10.1007/s13668-017-0219-5
Ocvirk, S., and O'Keefe, S. (2021). Dietary fat, BAs metabolism and colorectal cancer. Semin. Cancer Biol. 73, 347–355. doi: 10.1016/j.semcancer.2020.10.003
Oluwagbemigun, K., O'Donovan, A. N., Berding, K., Lyons, K., Alexy, U., Schmid, M., et al. (2021). Long-term dietary intake from infancy to late adolescence is associated with gut microbiota composition in young adulthood. Am. J. Clin. Nutr. 113, 647–656. doi: 10.1093/ajcn/nqaa340
Paik, D., Yao, L., Zhang, Y., Bae, S., D'Agostino, G. D., Zhang, M., et al. (2022). Human gut bacteria produce ΤΗ17-modulating BAs metabolites. Nature 603, 907–912. doi: 10.1038/s41586-022-04480-z
Petersen, C., Bell, R., Klag, K. A., Lee, S. H., Soto, R., Ghazaryan, A., et al. (2019). T cell-mediated regulation of the microbiota protects against obesity. Science 365:eaat9351. doi: 10.1126/science.aat9351
Qiu, Y., Xu, G. E., and Xiao, J. (2021). High-fat/high-sugar diet causes metabolic disorders and cardiac dysfunction: a new model in primates. J. Cardiovasc. Transl. Res. 14, 1019–1020. doi: 10.1007/s12265-021-10119-7
Ridlon, J. M., Harris, S. C., Bhowmik, S., Kang, D. J., and Hylemon, P. B. (2016). Consequences of bile salt biotransformations by intestinal bacteria. Gut Microbes 7, 22–39. doi: 10.1080/19490976.2015.1127483
Sato, Y., Atarashi, K., Plichta, D. R., Arai, Y., Sasajima, S., Kearney, S. M., et al. (2021). Novel BAs biosynthetic pathways are enriched in the microbiome of centenarians. Nature 599, 458–464. doi: 10.1038/s41586-021-03832-5
Sayin, S. I., Wahlström, A., Felin, J., Jäntti, S., Marschall, H. U., Bamberg, K., et al. (2013). Gut microbiota regulates BAs metabolism by reducing the levels of tauro-beta-muricholic acid, a naturally occurring FXR antagonist. Cell Metab. 17, 225–235. doi: 10.1016/j.cmet.2013.01.003
Sengupta, A., and Radhakrishna, S. (2016). Influence of fruit availability on fruit consumption in a generalist primate, the rhesus macaque macaca mulatta. Int. J. Primatol. 37, 703–717. doi: 10.1007/s10764-016-9933-x
Singh, J., Metrani, R., Shivanagoudra, S. R., Jayaprakasha, G. K., and Patil, B. S. (2019). Review on BAs: effects of the gut microbiome, interactions with dietary fiber, and alterations in the bioaccessibility of bioactive compounds. J. Agric. Food Chem. 67, 9124–9138. doi: 10.1021/acs.jafc.8b07306
Spence, C., Wells, W. G., and Smith, C. J. (2006). Characterization of the primary starch utilization operon in the obligate anaerobe Bacteroides fragilis: regulation by carbon source and oxygen. J. Bacteriol. 188, 4663–4672. doi: 10.1128/JB.00125-06
Struhsaker, T. T. (1995). Colobine monkeys. Their ecology, behaviour and evolution. Int. J. Primatol. 16, 1035–1037. doi: 10.1007/BF02696118
Tanaka, M., Sanefuji, M., Morokuma, S., Yoden, M., Momoda, R., Sonomoto, K., et al. (2020). The association between gut microbiota development and maturation of intestinal BAs metabolism in the first 3 y of healthy Japanese infants. Gut Microbes 11, 205–216. doi: 10.1080/19490976.2019.1650997
Thursby, E., and Juge, N. (2017). Introduction to the human gut microbiota. Biochem. J. 474, 1823–1836. doi: 10.1042/BCJ20160510
Trefflich, I., Marschall, H. U., Giuseppe, R. D., Ståhlman, M., Michalsen, A., Lampen, A., et al. (2019). Associations between dietary patterns and BAs-results from a cross-sectional study in vegans and omnivores. Nutrients 12:47. doi: 10.3390/nu12010047
Turnbaugh, P. J., Ley, R. E., Mahowald, M. A., Magrini, V., Mardis, E. R., and Gordon, J. I. (2006). An obesity-associated gut microbiome with increased capacity for energy harvest. Nature 444, 1027–1031. doi: 10.1038/nature05414
Vital, M., Rud, T., Rath, S., Pieper, D. H., and Schlüter, D. (2019). Diversity of bacteria exhibiting BAs-inducible 7α-dehydroxylation genes in the human gut. Comput. Struct. Biotechnol. J. 17, 1016–1019. doi: 10.1016/j.csbj.2019.07.012
Wahlström, A., Sayin, S. I., Marschall, H. U., and Bäckhed, F. (2016). Intestinal crosstalk between BAs and microbiota and its impact on host metabolism. Cell Metab. 24, 41–50. doi: 10.1016/j.cmet.2016.05.005
Wan, Y., Yuan, J., Li, J., Li, H., Zhang, J., Tang, J., et al. (2020). Unconjugated and secondary BAs profiles in response to higher-fat, lower-carbohydrate diet and associated with related gut microbiota: a 6-month randomized controlled-feeding trial. Clin. Nutr. 39, 395–404. doi: 10.1016/j.clnu.2019.02.037
Wang, S., Martins, R., Sullivan, M. C., Friedman, E. S., Misic, A. M., et al. (2019). Diet-induced remission in chronic enteropathy is associated with altered microbial community structure and synthesis of secondary BAs. Microbiome 7:126. doi: 10.1186/s40168-019-0740-4
Winston, J. A., and Theriot, C. M. (2020). Diversification of host BAs by members of the gut microbiota. Gut Microbes 11, 158–171. doi: 10.1080/19490976.2019.1674124
Yu, H., Ni, Y., Bao, Y., Zhang, P., Zhao, A., Chen, T., et al. (2015). Chenodeoxycholic acid as a potential prognostic marker for roux-en-y gastric bypass in Chinese obese patients. J. Clin. Endocrinol. Metab. 100, 4222–4230. doi: 10.1210/jc.2015-2884
Zeng, H., Umar, S., Rust, B., Lazarova, D., and Bordonaro, M. (2019). Secondary BAs and short chain fatty acids in the colon: a focus on colonic microbiome, cell proliferation, inflammation, and cancer. Int. J. Mol. Sci. 20:1214. doi: 10.3390/ijms20051214
Keywords: bile acids, diet, gut microbiome, primates, metabolomics
Citation: Li X, Wang X, Wang Z, Zhang M, Wang S, Xiang Z, Pan H and Li M (2022) The Relationship Between Gut Microbiome and Bile Acids in Primates With Diverse Diets. Front. Microbiol. 13:899102. doi: 10.3389/fmicb.2022.899102
Edited by:
Anna Gałązka, Institute of Soil Science and Plant Cultivation, PolandReviewed by:
Artur Banach, The John Paul II Catholic University of Lublin, PolandIram Liaqat, Government College University, Lahore, Pakistan
Copyright © 2022 Li, Wang, Wang, Zhang, Wang, Xiang, Pan and Li. 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: Huijuan Pan, cGhqamFuaW5lQGJqZnUuZWR1LmNu
†These authors have contributed equally to this work