Front. Endocrinol. Frontiers in Endocrinology Front. Endocrinol. 1664-2392 Frontiers Media S.A. 10.3389/fendo.2023.1250865 Endocrinology Original Research Species variations in the gut microbiota of captive snub-nosed monkeys Xi Li 1 2 Han Jincheng 1 2 Wen Xiaohui 3 * Zhao Longfei 1 * Qin Xinxi 1 4 * Luo Shengjun 3 Lv Dianhong 3 Song Shuai 3 1 College of Biology and Food, Shangqiu Normal University, Shangqiu, China 2 Henan Engineering Research Center of Development and Application of Green Feed Additives, College of Biology and Food, Shangqiu Normal University, Shangqiu, China 3 Institute of Animal Health, Guangdong Academy of Agricultural Sciences, Guangzhou, China 4 College of Veterinary Medicine, Northwest A&F University, Yangling, China

Edited by: Fazul Nabi, Southwest University, China

Reviewed by: Xin-He Lai, University of Colorado Anschutz Medical Campus, United States; Nurgulsim Kaster, S. Seifullin Kazakh Agro Technical University, Kazakhstan; Reetesh Kumar, GLA University, India

*Correspondence: Xiaohui Wen, wenxiaohui@gdaas.cn; Longfei Zhao, hnzhaolongfei@163.com; Xinxi Qin, qinxinxi@163.com

13 09 2023 2023 14 1250865 30 06 2023 29 08 2023 Copyright © 2023 Xi, Han, Wen, Zhao, Qin, Luo, Lv and Song 2023 Xi, Han, Wen, Zhao, Qin, Luo, Lv and Song

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.

 Introduction

Snub-nosed monkeys are species in danger of extinction due to habitat fragmentation and human activities. Captivity has been suggested as an Auxiliary Conservation Area (ASA) strategy. However, little is known about the adaptation of different species of snub-nosed monkeys to captive environments.

 Methods

This study compared the gut microbiota between Rhinopithecus bieti, R. brelichi, and R. roxellana under identical captive conditions to provide insights for improving captive conservation strategies.

 Results

The results showed that these three Rhinopithecus species shared 80.94% of their Operational Taxonomic Unit (OTU), indicating high similarity in gut microbiota composition. The predominant phyla were Firmicutes and Bacteroidetes for all three Rhinopithecus species, but differences were observed in diversity, characteristic bacterial communities, and predicted function. Significant enrichment of cellulolytic families, including Ruminococcaceae, Clostridiales vadinBB60 group, Christensenellaceae, and Erysipelotrichaceae, and pathways involved in propionate and butyrate metabolism in the gut of R. bieti suggested that it may have a superior dietary fiber utilization capacity. In contrast, Bacteroidetes, Ruminoccaceae, and Trichospiraceae were more abundant in R. brelichi and R. roxellana, and were associated with saccharide and glycan metabolic pathways. Moreover, R. brelichi and R. roxellana also had higher similarity in microbiota composition and predicted function.

 Discussion

In conclusion, the results demonstrate that host species are associated with the composition and function of the gut microbiota in snub-nosed monkeys. Thus, host species should be considered when formulating nutritional strategies and disease surveillance in captive snub-nosed monkeys.

 snub-nosed monkey gut microbiota captivity species conservation section-in-acceptance Gut Endocrinology Introduction

Hosts and gut microbiota rely on each other and maintain a dynamic balance through long-term mutual adaptation, forming a mutualistic relationship. This dynamic equilibrium is important for maintaining host health and metabolic stability, particularly the intestinal environment’s homeostasis (1). Gut microbes can provide substrates and enzymes for host metabolism, help digest nutrients in food, synthesize vitamins, and regulate energy metabolism (2, 3). The host’s physiological structure and immune system largely regulate and restrain the colonization of gut microbes, and the gut microbiota composition changes based on the host’s diet, physiological state, and external environment (4).

Snub-nosed monkeys belong to Rhinopithecus in the Primate, Cercopithecidae and Colobinae, which includes five endangered species: R. bieti, R. brelichi, R. roxellana, R. strykeri, and R. avunculu. Zhou et al. (2014) found that the northern species (R. brelichi and R. roxellana) and the Himalayan species (R. bieti and R. strykeri) diverged about 1.6 million years ago (5). Snub-nosed monkeys are typical phytophagous animals that mainly feed on shoots, young buds, fresh leaves, mosses and fruits. Snub-nosed monkeys have a multicapsular stomach, and the microbes that live in the stomach can promote the digestion of plant fibers (5). In addition, the hindgut of snub-nosed monkeys is also dominated by bacteria capable of producing complex carbohydrate-degrading enzymes (6), which helps snub-nosed monkeys to maintain physiological functions and consistently gain energy from high-fiber foods during the long winters and low temperatures (7).

Habitat variations often trigger adaptive changes in the gut microbiota of snub-nosed monkeys. Hale et al. (2019) found that the richness of gut microbiota and the population of cellulolytic bacteria (Lachnospiraceae and Ruminococcaceae) decreased significantly in captive R. brelichi compared to their wild counterparts, while Prevotella and Bacteroides, which can break down simple sugars and carbohydrates, increased significantly to adapt to the captive diet (8). Wang et al. (2021) similarly observed lower diversity and richness of gut microbiota in captive R. roxellana compared to the wild ones, accompanied by a significant increase in Prevotella/Bacteroides ratio, suggesting an enhanced ability to digest simple sugars (9). Even when their habitat remains unchanged, artificial feeding still affects the gut microbiota of snub-nosed monkeys. Feeding peanuts, apples, and carrots to wild R. roxellana and R. bieti significantly reduced Firmicutes population and increased Bacteroidetes population (10, 11). Therefore, while captivity or artificial assisted feeding can be important for the conservation of snub-nosed monkeys, changes in habitat and food structure can affect their physiology and gut health. R. bieti, R. brelichi, and R. roxellana are only distributed in the mountainous areas of southwest China. Investigating the adaptation of these three Rhinopithecus species to captive environments can facilitate the development of conservation strategies for them. In this study, Illumina NovaSeq high-throughput sequencing was used to identify the characteristics of gut microbiota in R. bieti, R. brelichi, and R. roxellana that were captivated in the Beijing Zoo. By comparing similarities and differences in their microbial composition and predicted functions, this research provides a reference for disease surveillance and dietary adjustment during the conservation of captive snub-nosed monkeys.

 Materials and methods Sample collection

Fresh feces were collected from male R. bieti, R. brelichi, and R. roxellana, aged between 2-9 years old, with a sample size of 5 for each species, in the Beijing Zoo in February 2023. The snub-nosed monkeys were kept in three adjacent enclosures and provided with a uniform diet including fresh leaves, fruits, vegetables, cooked corn cakes, eggs, and beef strips, etc. Three months prior to stool collection, no snub-nosed monkeys had a history of disease or medication.

 High-throughput sequencing of the fecal microbiota

The total DNA of fecal microbiota was extracted using the GenElute™ Stool DNA Isolation Kit (Sigma-Aldrich, USA) according to the instructions. DNA quality was detected using a Nanodrop 2000 ultra-trace spectrophotometer (Thermo Scientific, USA) and 1% agarose gel electrophoresis. The V3-V4 region of the bacterial 16S rRNA gene was amplified using universal primers (338F: 5’-ACTCCTACGGGAGGCAGCAG-3’ and 806R: 5’-GGACTACHVGGGTWTCTAAT-3’) and three repeats for each sample. The PCR was performed in a total reaction volume of 10 μL: template DNA 10 ng, 2.0 mmol/L dNTPs 2 μL, KOD FX Neo Buffer 5 μL, KOD FX Neo 0.2 μL, 10 μmol/L primer 338F 0.3 μL, 10 μmol/L primer 806R 0.3 μL, and ddH2O up to 10 μL.The amplification procedure consisted of pre-denaturation at 95°C for 3 min, followed by 25 cycles of denaturation at 95°C for 30 s, annealing at 50°C for 30 s, and extension at 72°C for 40 s. Lastly, a final step was performed at 72°C for 7 min. The triplicate products were mixed and electrophoresed in 2% agarose gel at 110 V for 20 min for quality detection ( Supplementary Figure 1 ). The PCR recovered products were mixed equally according to their concentration. After clone libraries were constructed using the TruSeq® Nano DNA Kit (Illumina, USA), 2 × 250 bp paired-end high-throughput sequencing was performed using the Illumina NovaSeq platform (Illumina, USA).

 Bioinformatics analysis

Quality control of the sequences was performed through various steps. Firstly, low-quality sequences were removed using Trimmomatic v0.33 software (window size: 50 bp) (12). Primer sequences were then eliminated using Cutadapt v1.9.1 software (Maximum mismatch accepted: 20%; Minimum coverage: 80%) (13). Paired-end reads were assembled using Usearch v10 software (Minimum length of overlap: 10 bp; Minimum similarity within overlapping region: 90%; Maximum mismatch accepted: 5 bp) (14). For OTU clustering and taxonomic annotation of quality control sequences, QIIME 2 software (15) was used to analyze sequences with a similarity of 97% based on the SILVA database (Release132, http://www.arb-silva.de) (16). Alpha diversity indices were calculated using the QIIME 2 software. To assess the species richness of bacterial communities, the ACE index was adopted, with higher values indicating a richer community (17). Furthermore, the Shannon index, which considers both species richness and evenness of bacterial communities, was employed to assess alpha diversity, with higher values indicating greater diversity (18). One-way ANOVA with Tukey’s post-hoc test was subsequently employed to analyze the differences in alpha diversity among the three Rhinopithecus species. The Principal Coordinates Analysis (PCoA) (19) and Unweighted Pair-group Method with Arithmetic Mean (UPGMA) (20) based on the weighted UniFrac distance were used to analyze the beta diversity. The microbiota differences and similarities between the three Rhinopithecus species were analyzed using two statistical methods, specifically Permutational Multivariate Analysis of Variance (PERMANOVA) (21) and Analysis of Similarities (ANOSIM) (22), based on the weighted UniFrac distances in the vegan package of R software. Linear discriminant analysis (LDA) effect size (LEfSe) analysis (23) was used to look for biomarkers that have significant statistical differences among these three species. The Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2 (PICRUSt 2) software (24) was used to predict the functional gene composition of gut microbiota in the three Rhinopithecus species based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (25).

 Results Sequencing information

Of the 1,182,216 effective reads obtained from the samples, 745 OTUs were identified. According to the Venn diagram ( Figure 1 ), 603 OTUs were shared among the three Rhinopithecus species, largely belonging to Firmicutes (433 OTUs) and Bacteroidetes (61 OTUs). The core bacterial families identified were Ruminococcaceae, Lachnospiraceae, Christensenellaceae, Muribaculaceae, Rikenellaceae, Clostridiales vadinBB60 group, and Bacteroidaceae. Furthermore, R. bieti, R. brelichi, and R. roxellana had 27, 7, and 5 unique OTUs, respectively.

Shared and unique OTUs among the three Rhinopithecus species were visualized by Venn diagram.

 Analysis of the gut microbiota diversity

As shown in Figure 2 , the gut microbiota of R. roxellana displayed significantly higher ACE and Shannon index values than R. bieti and R. brelichi, indicating that the gut microbiota of R. roxellana was more diverse and richer in microbial taxa. R. brelichi, however, displayed the lowest alpha diversity. The weighted UniFrac distance-based PCoA plot in Figure 3A showed that the gut microbiota of the Rhinopithecus species varied significantly. PERMANOVA analysis (R. bieti vs. R. brelichi: R2 = 0.487, P = 0.007; R. bieti vs. R. roxellana: R2 = 0.609, P = 0.007; R. brelichi vs. R. roxellana: R2 = 0.320, P = 0.016) indicated that the most significant difference in gut microbiota existed between R. bieti and R. roxellana. In additiona, the UPGMA clustering tree, based on the same distance measure, illustrated a higher similarity in species composition between R. brelichi and R. roxellana ( Figure 3B ). This similarity was also verified by the ANOSIM analysis (R. bieti vs. R. brelichi: R = 0.788, P = 0.007; R. bieti vs. R. roxellana: R = 1.000, P = 0.007; R. brelichi vs. R. roxellana: R = 0.480, P = 0.027). Comparative genomics assigned R. brelichi and R. roxellana to northern species and R. bieti to Himalayan species (5), suggesting a mutual evolution between the monkeys and their gut microbiota.

Differences in Alpha diversity of gut microbiota among the three Rhinopithecus species. (A) Pairwise comparisons of the ACE index between the groups. (B) Pairwise comparisons of the Shannon index between the groups. One-way ANOVA with Tukey’s post-hoc test: *P < 0.05, **P < 0.01.

Beta diversity of the gut microbiota in captive snub-nosed monkeys. (A) PCoA plot. (B) UPGMA clustering tree. Weighted UniFrac algorithm.

 Comparative analysis of the gut microbiota composition

The top 10 phyla in the gut microbiota of the three Rhinopithecus species were shown in Figure 4A . The predominant phyla were Firmicutes and Bacteroidetes, accounting for 78.44% in R. bieti, 82.56% in R. brelichi, and 86.38% in R. roxellana. The relative abundance of Firmicutes in the gut of R. brelichi was significantly lower than that of R. bieti and R. roxellana (P < 0.01 for R. brelichi vs. R. bieti; P < 0.05 for R. brelichi vs. R. roxellana), but it had the highest proportion of Bacteroidetes (34.49% for R. brelichi; 29.94% for R. roxellana; 17.57% for R. bieti, P < 0.01 for R. brelichi vs. R. bieti, P > 0.05 for R. brelichi vs. R. roxellana). Among the other nondominant phyla, Kiritimatiellaeota was highest in R. bieti (P < 0.05 for R. bieti vs. R. brelichi and R. bieti vs. R. roxellana), Proteobacteria was lowest in R. brelichi (P < 0.05 for R. brelichi vs. R. bieti and R. brelichi vs. R. roxellana), and Fibrobacteres, Elusimicrobia, and Acidobacteria were highest in R. roxellana (P < 0.05 for R. roxellana vs. R. bieti and R. roxellana vs. R. brelichi) ( Supplementary Table 1 ). At the genus level, the gut microbiota of R. bieti, R. brelichi, and R. roxellana were dominated by Ruminococcaceae UCG-005 (11.06%), uncultured f Muribaculacea (16.04%), and uncultured f Lachnospiraceae (10.52%), respectively ( Figure 4B ). Out of the top 50 genera, 25 showed significant differences between the three Rhinopithecus species ( Supplementary Table 2 ). LEfSe analysis showed that the biomarkers of R. bieti were Ruminococcaceae (including Ruminococcaceae UCG-005 and Ruminococcaceae UCG-002), Clostridiales vadinBB60 group (uncultured f Clostridiales vadinBB60 group), Christensenellaceae (Christensenellaceae R-7 group), Erysipelotrichaceae (Erysipelotrichaceae UCG-004) belonging to Firmicutes; Kiritimatiellae belonging to Kiritimatiellaeota; uncultured family and genus of WCHB1-41 belonging to Verrucomicrobia; and Rikenellaceae RC9 gut group belonging to Bacteroidetes. Meanwhile, Muribaculaceae (uncultured f Muribaculacea) and Bacteroidaceae (Bacteroides) belonging to Bacteroidetes; Spirochaetaceae (Treponema 2) belonging to Spirochaetes; Eubacterium coprostanoligenes group belonging to Firmicutes were significantly enriched in R. brelichi. The uncultured f Lachnospiraceae belonging to Firmicutes and Prevotella 9 belonging to Bacteroidetes were significantly enriched in R. roxellana ( Figure 5 ).

Species distribution of the gut microbiota in captive snub-nosed monkeys. (A) The distribution histogram of the top ten phyla in the three groups. (B) The distribution histogram of the top ten genera in the three groups.

LDA value distribution histogram of the gut microbiota in the captive snub-nosed monkeys. The default LDA threshold is 4.0. Larger LDA values indicate a greater influence of the abundance of bacterial species on the differences in the community between groups.

 Functional prediction of the gut microbiota

High-throughput sequencing data of the 16S rRNA V3-V4 gene were analyzed using PICRUSt 2 software based on the KEGG database. The comparison between different Rhinopithecus species indicated that the proportion of 15 functional genes varied significantly between R. bieti and R. brelichi ( Figure 6A ). Specifically, The relative abundance of amino acid metabolism, terpenoids and polyketides metabolism, xenobiotics biodegradation and metabolism, and the global and overview maps pathways of R. bieti were significantly higher than that of R. brelichi. However, R. brelichi showed a significantly higher proportion of the biosynthesis of other secondary metabolites and the glycan biosynthesis and metabolism pathways when compared to R. bieti. In addition, within human disease pathways, R. bieti showed a higher proportion of parasitic and viral infectious disease, and cardiovascular disease pathways, while R. brelichi exhibited a higher proportion of antimicrobial drug resistance, and immune disease pathways. There were 11 functional genes that differed significantly between R. bieti and R. roxellana ( Figure 6B ). Similar to R. brelichi, the glycine biosynthesis and metabolism, antimicrobial drug resistance, and immune diseases were significantly higher in R. roxellana than in R. bieti. Furthermore, other pathways of amino acid metabolism and carbohydrate metabolism were also significantly higher in R. roxellana than in R. bieti. No significant difference was observed in carbohydrate metabolism between R. brelichi and R. roxellana. Overall, R. bieti may have different carbohydrate metabolic pathway strategies from the other two Rhinopithecus species ( Supplementary Table 3 ). There were no significant differences in other secondary and tertiary metabolic pathways between the gut microbiota of R. brelichi and R. roxellana.

Differences in the KEGG secondary metabolic pathways between the three Rhinopithecus species. (A) The proportion of different functions between R. bieti and R. brelichi. (B) The abundance proportion of different functions between R. bieti and R. roxellana. One-way ANOVA (Tukey’s post-hoc test), the P-value threshold was 0.05.

 Discussion

The gut microbiota of the three Rhinopithecus species in captive conditions have different characteristics and differences that can be informative for developing captive diet structures, and disease surveillance measures necessary to better protect these endangered species. The study showed that the predominant phyla in snub-nosed monkeys were the Firmicutes and Bacteroidetes, which matches earlier studies (26, 27). Most of the bacterial communities in the Firmicutes can use cellulose, hemicellulose, and xylan as energy sources (28). The relative abundance of Firmicutes in R. bieti was significantly higher than in R. brelichi and R. roxellana, and it was enriched with cellulose or hemicellulose-degrading bacterial families, including Ruminococcaceae (29), Clostridiales vadinBB60 group (30), Christensenellaceae (31) and Erysipelotrichaceae (32). The proportion of the Propanoate and Butanoate metabolism pathway was also significantly higher in R. bieti. It is speculated that R. bieti may have a stronger ability to decompose plant fibers to produce energy materials. The Bacteroidetes has rich polysaccharide lyase and glycoside hydrolyase genes necessary to break down polysaccharides (33, 34). The Bacteroidetes had the highest abundance in R. brelichi, followed by R. roxellana, and R. bieti showed the lowest percentage. Muribaculaceae and Bacteroidaceae were enriched in R. brelichi, and Prevotella 9 was enriched in R. roxellana. Muribaculaceae bacteria are the main users of mucosal sugars (35). Members of Bacteroidaceae provide nutrients to the host by breaking down different glycans (36, 37). The glycan biosynthesis and metabolism pathway was found to be significantly more pronounced in R. brelichi and R. roxellana, as compared to R. bieti. Moreover, these two species demonstrated greater enrichment of important carbohydrate tertiary metabolism pathways such as pentose, glucuronate, starch, and sucrose metabolism, indicating species-specific strategies for carbohydrate metabolism among snub-nosed monkeys.

Apart from the divergence of the predominant phyla, 57.14% of the nondominant phyla showed host differences. Kiritimatiellaeota was found to be the third largest phylum in the gut of R. bieti with the highest proportion of genus belonging to order WCHB 1-41. The arginine and fatty acid biosynthesis pathways encoded by the WCHB 1-41 genome were found to play a role in host energy metabolism and nitrogen utilization (38). The comparison of lipid metabolism pathways among the three species of snub-nosed monkeys revealed that the proportion of the fatty acid metabolism in the gut microbiota of R. bieti was significantly higher than that in R. brelichi and R. roxellana. The proportion of Spirochaetes in R. brelichi was significantly higher than that in R. roxellana (P < 0.05). Treponema 2 was found to be the main genus responsible for this difference (9.69% for R. brelichi vs. 1.22% for R. roxellana, P < 0.05). Treponema 2 has enzymes that mediate pyruvate oxidation and decarboxylation to enter the citrate cycle, promoting the biosynthesis of arginine and fatty acids (39). The proportion of Proteobacteria was significantly higher in R. roxellana, and the relative abundance of the genus Parasutterella was significantly higher than in R. bieti and R. brelichi (P < 0.01). The Parasutterell occupied a specific gut niche and could cause changes in bacterial metabolites, such as aromatic amino acids and bile acids in the mouse gut (40). These low-abundance bacterial communities in this study not only caused differences in the composition of the consortium but might also have played a relevant role in interspecies differences in metabolic pathways.

Studies have revealed that northern species (R. brelichi and R. roxellana) and the Himalayan species (R. bieti and R. strykeri) diverged about 1.6 million years ago (5). R. bieti, which is the largest monkey in the genus Rhinopithecus, inhabits virgin alpine forests between 2,500-5,000 meters in southeastern Tibet and northwestern Yunnan (41). R. strykeri can only be found in the virgin forests of 1,700-3,100 meters in northwestern Yunnan and northeastern Myanmar (42). R. brelichi lives in forests on Fanjing Mountain at an altitude of 500-800 meters in Guizhou Province (43). R. roxellana is distributed in forests at an altitude of 1,500-3,300 meters in Sichuan, Gansu, Shaanxi and Hubei provinces of China (44). R. bieti, R. strykeri, and R. brelichi primarily feed on plant such as leaves, shoots, buds, fruits, bark, and lichen (45–47). R. roxellana has a more diverse diet and occasionally preys on small animals, birds, eggs, or insects (48). Long-term dietary and habitat differences encourage different hosts to adopt distinct strategies to acquire microbes from nature and attain a healthy and balanced symbiotic relationship (49). This symbiotic relationship makes the microbial communities species-specific and resistant to disturbance over time (50). An analysis of adaptive variations in the gut microbiota of 18 nonhuman primates revealed that the physiological evolution of the host has a greater effect on the construction of the microbiota than the dietary niche (51). However, differences were found when comparing the gut microbiota of wild R. bieti, R. strykeri, and R. roxellana by Wang et al. (2023). Although the gut microbiota of the three Rhinopithecus species had similar core bacterial communities involved in degrading cellulose, hemicellulose and lignin, R. strykeri had significantly higher levels of enzymes related to pectin and glucose metabolism than R. bieti and R. roxellana. This difference may be attributed to R. strykeri’s greater consumption of fruits, seeds, buds, and flowers, leading to adaptive changes in its gut microbiota. The gut microbiota of R. bieti and R. roxellana, which share similar lifestyles, were more alike than that of R. strykeri (52). Therefore, it can be seen that host species, diet, and habitat have significant effects on the gut microbiota of nonhuman primates. In this study, all snub-nosed monkeys were kept at the Beijing Zoo under the same living environment and dietary structure, which minimized the influence of environmental and dietary factors on the gut microbiota. All samples shared 80.94% of OTUs, mainly annotated to Firmicutes and Bacteroidetes, indicating that snub-nosed monkeys have the same predominant phyla. PCoA and UPGMA clustering analysis showed that R. brelichi and R. roxellana had a more similar bacterial community composition, and the differential analysis of the KEGG metabolic pathways between these two species did not show significant differences in the metabolic function. However, R. bieti (Himalayan species) showed greater differences in both microbial composition and function. Therefore, the profound effects of host species on gut microbiota stability and adaptability must be considered when studying nutritional strategies or the effects of niche on the gut microbiota in snub-nosed monkeys.

 Conclusions

This study investigated the gut microbiota characteristics of three endangered snub-nosed monkey species (R. bieti, R. brelichi, and R. roxellana) in China, under identical captive conditions. Although the predominant phyla of the gut microbiota of the three Rhinopithecus species were Firmicutes and Bacteroidetes, the proportions and bacterial communities were found to differ among them. Since most of the bacteria within Firmicutes and Bacteroidetes are involved in carbohydrate metabolism, differences in carbohydrate metabolism pathways were identified between R. bieti and the other two Rhinopithecus species. Furthermore, 57.14% of the non-dominant phyla were found to differ significantly between the three species. This affected not only microbiota composition but also played a role in the formation of interspecific differences in metabolic function. Finally, R. brelichi and R. roxellana exhibited greater similarity in gut microbiota composition and metabolic function than R. bieti, which is evolutionarily more distant. In conclusion, we obtained valuable insights into the gut microbiota characteristics of endangered snub-nosed monkeys. This study demonstrates the association of host species with gut microbiota, which can be applied to disease surveillance and development of dietary protocols in captive endangered snub-nosed monkeys.

 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 Ethics Committee of the Shangqiu Normal University. The study was conducted in accordance with the local legislation and institutional requirements.

 Author contributions

LX designed the experiments and wrote the manuscript. SL and DL collected the samples. JH, SS, and XQ analyzed the data. XW and LZ revised the manuscript. All authors contributed to the article and approved the submitted version.

 Funding

This study was funded by the key research and development projects of Henan Province (222102320024); Guangdong provincial special fund for modern agriculture industry technology innovation teams (Department of Agriculture and Rural Affairs of Guangdong Province: 2023KJ119); and National natural science foundation of China (32072753).

 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/fendo.2023.1250865/full#supplementary-material

References Duncan SH Conti E Ricci L Walker AW . Links between diet, intestinal anaerobes, microbial metabolites and health. Biomedicines (2023) 11:1338. doi: 10.3390/Biomedicines11051338 Meyers RG Samouda H Bohn T . Short chain fatty acid metabolism in relation to gut microbiota and genetic variability. Nutrients (2022) 14:5361. doi: 10.3390/Nu14245361 Zhan Q Wang R Thakur K Feng JY Zhu YY Zhang JG . Unveiling of dietary and gut-microbiota derived B vitamins: metabolism patterns and their synergistic functions in gut-brain homeostasis. Crit Rev Food Sci Nutr (2022), 1–13. doi: 10.1080/10408398.2022.2138263 Flynn JK Ortiz AM Herbert R Brenchley JM . Host genetics and environment shape the composition of the gastrointestinal microbiome in nonhuman primates. Microbiol Spectr (2022) 11:E0213922. doi: 10.1128/Spectrum.02139-22 Zhou X Wang B Pan Q Zhang J Kumar S Sun X . Whole-genome sequencing of the snub-nosed monkey provides insights into folivory and evolutionary history. Nat Genet (2014) 46:1303–10. doi: 10.1038/Ng.3137 Liu R Amato K Hou R Gomez A Dunn DW Zhang J . Specialized digestive adaptations within the hindgut of A colobine monkey. Innovation (Camb) (2022) 3:100207. doi: 10.1016/J.Xinn.2022.100207 Li Y Yan Y Fu H Jin S He S Wang Z . Does diet or macronutrients intake drive the structure and function of gut microbiota? Front Microbiol (2023) 14:1126189. doi: 10.3389/Fmicb.2023.1126189 Hale VL Tan CL Niu K Yang Y Zhang Q Knight R . Gut microbiota in wild and captive guizhou snub-nosed monkeys, Rhinopithecus Brelichi . Am J Primatol (2019) 81:E22989. doi: 10.1002/Ajp.22989 Wang X Wang Z Pan H Qi J Li D Zhang L . Captivity influences the gut microbiome of Rhinopithecus Roxellana . Front Microbiol (2021) 12:763022. doi: 10.3389/Fmicb.2021.763022 Li H Xia W Liu X Wang X Liu G Chen H . Food provisioning results in functional, but not compositional, convergence of the gut microbiomes of two wild Rhinopithecus species: evidence of functional redundancy in the gut microbiome. Sci Total Environ (2023) 858:159957. doi: 10.1016/J.Scitotenv.2022.159957 Xia W Liu G Wang D Chen H Zhu L Li D . Functional convergence of yunnan snub-nosed monkey and bamboo-eating panda gut microbiomes revealing the driving by dietary flexibility on mammal gut microbiome. Comput Struct Biotechnol J (2022) 20:685–99. doi: 10.1016/J.Csbj.2022.01.011 Bolger AM Lohse M Usadel B . Trimmomatic: A flexible trimmer for illumina sequence data. Bioinformatics (2014) 30:2114–20. doi: 10.1093/Bioinformatics/Btu170 Martin M . Cutadapt removes adapter sequences from high-throughput sequencing reads. Emb Net J (2011) 17:10–2. doi: 10.14806/Ej.17.1.200 Edgar RC . Uparse: highly accurate otu sequences from microbial amplicon reads. Nat Methods (2013) 10:996–8. doi: 10.1038/Nmeth.2604 Bolyen E Rideout Jr Dillon MR Bokulich NA Abnet CC Al-Ghalith GA . Reproducible, interactive, scalable and extensible microbiome data science using qiime 2. Nat Biotechnol (2019) 37:852–7. doi: 10.1038/S41587-019-0209-9 Quast C Pruesse E Yilmaz P Gerken J Schweer T Yarza P . The silva ribosomal rna gene database project: improved data processing and web-based tools. Nucleic Acids Res (2013) 41:590–6. doi: 10.1093/Nar/Gks1219 Chao A Yang MCK . Stopping rules and estimation for recapture debugging with unequal failure rates. Biometrika (1993) 80:193–201. doi: 10.2307/2336768 Shannon CE . A mathematical theory of communication. Bell Sys Tech J (1948) 27:379–423. doi: 10.1002/J.1538-7305.1948.Tb01338.X Anderson MJ Willis TJ . Canonical analysis of principal coordinates: A useful method of constrained ordination for ecology. Ecology (2003) 84:511–25. doi: 10.1890/0012-9658(2003)084[0511:Caopca]2.0.Co;2 Sourdis J Krimbas C . Accuracy of phylogenetic trees estimated from dna sequence data. Mol Biol Evol (1987) 4:159–66. doi: 10.1016/0047-2484(87)90009-1 Anderson MJ . A new method for non-parametric multivariate analysis of variance. Aust Ecol (2001) 26:32–46. doi: 10.1111/J.1442-9993.2001.01070.Pp.X Clarke KR . Non-parametric multivariate analyses of changes in community structure. Aust J Ecol (1993) 18:117–43. doi: 10.1111/J.1442-9993.1993.Tb00438.X Segata N Izard J Waldron L Gevers D Miropolsky L Garrett WS . Metagenomic biomarker discovery and explanation. Genome Biol (2011) 12:R60. doi: 10.1186/Gb-2011-12-6-R60 Douglas GM Maffei VJ Zaneveld J Yurgel SN Brown Jr Taylor CM . Picrust2: an improved and extensible approach for metagenome inference. Biorxiv (2019). doi: 10.1101/672295 Kanehisa M Goto S . Kegg: kyoto encyclopedia of genes and genomes. Nucleic Acids Res (2000) 28:27–30. doi: 10.1093/Nar/28.1.27 Liu X Fan P Che R Li H Yi L Zhao N . Fecal bacterial diversity of wild sichuan snub-nosed monkeys (Rhinopithecus roxellana). Am J Primatol (2018) 80:E22753. doi: 10.1002/Ajp.22753 Xu B Xu W Li J Dai L Xiong C Tang X . Metagenomic analysis of the Rhinopithecus Bieti fecal microbiome reveals A broad diversity of bacterial and glycoside hydrolase profiles related to lignocellulose degradation. BMC Genomics (2015) 16:174. doi: 10.1186/S12864-015-1378-7 Mangi MH Hussain T Shahid MS Sabir N Kalhoro MS Zhou X . Effects of flaxseed and multi-carbohydrase enzymes on the cecal microbiota and liver inflammation of laying hens. Anim (Basel) (2021) 11:600. doi: 10.3390/Ani11030600 Biddle A Stewart L Blanchard J Leschine S . Untangling the genetic basis of fibrolytic specialization by lachnospiraceae and ruminococcaceae in diverse gut communities. Diversity (2013) 5:627–40. doi: 10.3390/D5030627 Baniel A Amato KR Beehner JC Bergman TJ Mercer A Perlman RF . Seasonal shifts in the gut microbiome indicate plastic responses to diet in wild Geladas . Microbiome (2021) 9:26. doi: 10.1186/S40168-020-00977-9 Waters JL Ley RE . The human gut bacteria christensenellaceae are widespread, heritable, and associated with health. BMC Biol (2019) 17:83. doi: 10.1186/S12915-019-0699-4 Ziemer CJ . Newly cultured bacteria with broad diversity isolated from eight-week continuous culture enrichments of cow feces on complex polysaccharides. Appl Environ Microbiol (2014) 80:574–85. doi: 10.1128/Aem.03016-13 Kibegwa FM Bett RC Gachuiri CK Machuka E Stomeo F Mujibi FD . Diversity and functional analysis of rumen and fecal microbial communities associated with dietary changes in crossbreed dairy cattle. PloS One (2023) 18:E0274371. doi: 10.1371/Journal.Pone.0274371 El Kaoutari A Armougom F Leroy Q Vialettes B Million M Raoult D . Development and validation of A microarray for the investigation of the cazymes encoded by the human gut microbiome. PloS One (2013) 8:E84033. doi: 10.1371/Journal.Pone.0084033 Pereira FC Wasmund K Cobankovic I Jehmlich N Herbold CW Lee KS . Rational design of A microbial consortium of mucosal sugar utilizers reduces Clostridiodes difficile colonization. Nat Commun (2020) 11:5104. doi: 10.1038/S41467-020-18928-1 Morita H Kano C Ishii C Kagata N Ishikawa T Hirayama A . Bacteroides uniformis and its preferred substrate, α-cyclodextrin, enhance endurance exercise performance in mice and human males. Sci Adv (2023) 9:Eadd2120. doi: 10.1126/Sciadv.Add2120 Zhang M Wang X Wang Z Mao S Zhang J Li M . Metatranscriptomic analyses reveal important roles of the gut microbiome in primate dietary adaptation. Genes (Basel) (2023) 14:228. doi: 10.3390/Genes14010228 Guo N Wu Q Shi F Niu J Zhang T Degen AA . Seasonal dynamics of diet-gut microbiota interaction in adaptation of yaks to life at high altitude. NPJ Biofilms Microbiomes (2021) 7:38. doi: 10.1038/S41522-021-00207-6 Liu H Han X Zhao N Hu L Wang X Luo C . The gut microbiota determines the high-altitude adaptability of tibetan wild asses (Equus Kiang) in qinghai-tibet plateau. Front Microbiol (2022) 13:949002. doi: 10.3389/Fmicb.2022.949002 Ju T Jy K Stothard P Willing BP . Defining the role of Parasutterella, A previously uncharacterized member of the core gut microbiota. Isme J (2019) 13:1520–34. doi: 10.1038/S41396-019-0364-5 Liu Z Ren B Wei F Long Y Hao Y Li M . Phylogeography and population structure of the yunnan snub-nosed monkey (Rhinopithecus Bieti) inferred from mitochondrial control region dna sequence analysis. Mol Ecol (2007) 16:3334–49. doi: 10.1111/J.1365-294x.2007.03383.X Chen YX Yu Y Li C Xiao ZS Zhou GW Zhang ZJ . Population and conservation status of A transboundary group of black snub-nosed monkeys (Rhinopithecus Strykeri) between China and Myanmar. Zool Res (2022) 43:523–7. doi: 10.24272/J.Issn.2095-8137.2021.424 Niu K Cl T Yang Y . Altitudinal movements of guizhou snub-nosed monkeys (Rhinopithecus Brelichi) in fanjingshan national nature reserve, China: implications for conservation management of A flagship species. Folia Primatol (Basel) (2010) 81:233–44. doi: 10.1159/000322115 Hou R Chapman CA Rothman JM Zhang H Huang K Guo S . The geometry of resource constraint: an empirical study of the golden snub-nosed monkey. J Anim Ecol (2022) 90:751–65. doi: 10.1111/1365-2656.13408 Huang ZP Scott MB Li YP Ren GP Xiang ZF Cui LW . Black-and-white snub-nosed monkey (Rhinopithecus Bieti) feeding behavior in A degraded forest fragment: clues to A stressed population. Primates (2017) 58:517–24. doi: 10.1007/S10329-017-0618-7 Guo YQ Zhou J Xie JH Garber PA Bruford M Ren BP . Altitudinal ranging of the guizhou golden monkey (Rhinopithecus Brelichi): patterns of habitat selection and habitat use. Global Ecol And Conserv (2018) 16:E00473. doi: 10.1016/J.Gecco.2018.E00473 Yang Y Groves C Garber P Wang X Li H Long Y . First insights into the feeding habits of the critically endangered black snub-nosed monkey, Rhinopithecus Strykeri (Colobinae, primates). Primates (2019) 60:143–53. doi: 10.1007/S10329-019-00717-0 Hou R He S Wu F Chapman CA Pan R Garber PA . Seasonal variation in diet and nutrition of the northern-most population of Rhinopithecus Roxellana . Am J Primatol (2018) 80:E22755. doi: 10.1002/Ajp.22755 Que T Pang X Huang H Chen P Wei Y Hua Y . Comparative gut microbiome in Trachypithecus Leucocephalus and other primates in Guangxi, China, based on metagenome sequencing. Front Cell Infect Microbiol (2022) 12:872841. doi: 10.3389/Fcimb.2022.872841 Mccord AI Chapman CA Weny G Tumukunde A Hyeroba D Klotz K . Fecal microbiomes of non-human primates in western Uganda reveal species-specific communities largely resistant to habitat perturbation. Am J Primatol (2014) 76:347–54. doi: 10.1002/Ajp.22238 Amato KR Sanders GJ Song SJ Nute M Metcalf JL Thompson LR . Evolutionary trends in host physiology outweigh dietary niche in structuring primate gut microbiomes. Isme J (2019) 13:576–87. doi: 10.1038/S41396-018-0175-0 Wang XC Zhang JL Pan HJ Chen YX Mao SX Qi JW . Unique characteristics of gut microbiota in black snub-nosed monkeys (Rhinopithecus Strykeri) reveal an enzymatic mechanism of adaptation to dietary vegetation. Zool Res (2023) 44:357–60. doi: 10.24272/J.Issn.2095-8137