Zhichao Zhou1
Liping Tang1
Liping Yan1
Huiping Jia1
Yu Xiong1Jin Shang1Changliang Shao2Qiangwei Zhang3Hongjun Wang3Lun He4
Defu Hu1
Dong Zhang1*- 1School of Ecology and Nature Conservation, Beijing Forestry University, Beijing, China
- 2Mt. Kalamaili Ungulate Nature Reserve, Changji, China
- 3Gansu Endangered Animals Protection Center, Wuwei, China
- 4China Wildlife Conservation Association, Beijing, China
To explore how the living environment influences the establishment of gut microbiota in different species, as well as the extent to which changes in the living environment caused by captive breeding affect wildlife’s gut microbiota and health, we used 16S rRNA gene amplicon sequencing and shotgun metagenomic sequencing to compare the gut microbiome of two species of threatened equids, the Przewalski’s Horse and the Asian wild ass, in the wild and captivity. The results revealed that different species of Equidae living in the same environment showed remarkable convergence of gut microflora. At the same time, captive populations exhibited significantly “unhealthy” microbiota, such as low Alpha diversity, high levels of potentially pathogenic bacteria and biomarkers of physical or psychological disease, and enrichment of microbial functions associated with exogenous exposure and susceptibility, implying that the artificial environment created by captivity may adversely impact the health of wildlife to some extent. Our findings demonstrate the importance of the environmental factors for the establishment of gut microbiota and host health and provide new insights into the conservation of wildlife in captivity from the perspective of the microbiome.
Introduction
The mammalian gut is home to a large number of symbiotic microorganisms that play a vital role in maintaining the normal life processes of the host animal, such as nutrient digestion and absorption, pathogen prevention, and toxic compound degradation (Roggenbuck et al., 2014; Hanning and Diaz-Sanchez, 2015; Blyton et al., 2019; Li et al., 2021; Peixoto et al., 2021), and even the ability of gut microbiota to produce neuroactive molecules that affect the host animal’s brain, thus causing alterations in mood and behavior (Vuong et al., 2017; Sarkar et al., 2018). The establishment of commensal microbiota is usually affected by the phylogeny of the host which means that taxonomically different hosts tend to form their unique microflora structure during the long-term evolutionary process (Ley et al., 2008; Groussin et al., 2017). However, there has been accumulating evidence in recent years that environmental factors such as diet, habitat, and lifestyle play a dominant role in the composition of host gut microbiota (Zarrinpar et al., 2014; Ren et al., 2017; Kartzinel et al., 2019). Changes in the living environment can have a significant impact on the host’s gut microbiota. Jin et al. (2021) demonstrated that as carnivores, giant pandas emerged with numerous lignin-digesting microbial taxa in their gut microbiota to adapt to bamboo as a special diet. Some studies even suggest that the same living environment might drive the gut microflora of distinct host species to convergence (Muegge et al., 2011; Angelakis et al., 2016; Clayton et al., 2016; Greene et al., 2019). For instance, a study on herbivores showed that the gut microbiota of yaks and Tibetan sheep, which share a high altitude, converged to adapt to the harsh plateau environment (Zhang et al., 2016). In non-human primates, sympatric gorillas and chimpanzees also have similar gut microflora (Moeller et al., 2013).
As an essential way to conserve threatened wildlife, captive breeding can effectively prevent species extinction (Bowkett, 2009). At the same time, the dramatic change in the living environment undergone by wildlife as they transfer from the wild to captivity becomes a quintessential model for exploring how environmental factors affect the gut microbiota of wildlife (McKenzie et al., 2017). Furthermore, numerous studies have confirmed that host microbiota disorders due to environmental alterations may cause various metabolic diseases in hosts, such as obesity, irritable bowel syndrome (IBS), and inflammatory bowel disease (IBD; Cox et al., 2015; Shreiner et al., 2015). Therefore, we can further understand the effects of captive breeding on wildlife health by examining differences in gut microbiota composition caused to environmental changes (West et al., 2019; Zhu et al., 2021). However, current comparative microbiological studies of wildlife in the wild and captivity have focused on single species (Gibson et al., 2019; Schmidt et al., 2019; Sun et al., 2019). Rarely have multiple species been compared simultaneously to examine whether uniform conditions of captivity have similar impacts on the gut microbiota of different species.
Therefore, we utilized 16S rRNA gene amplicon sequencing and shotgun metagenomic sequencing in this study to compare the diversity, composition, and function of the gut microbiome of two threatened Equidae species, the Przewalski’s Horse (Equus ferus przewalskii; PH) and the Asian wild ass (Equus hemionus; AWA), living in the wild and captive environments. To explore the effects of field sympatry and uniform captive environments on the gut microbiota of different species of Equidae and reveal the potential impact of anthropogenic alterations in the living environment on the health of threatened equids from the perspective of the microbiome, allowing for a more comprehensive understanding of the captive breeding program for wildlife conservation management.
Materials and Methods
Study Area and Sample Collection
All fecal samples were collected from the adult in September of the same year to avoid the effect of ages and different seasons on the microbiota of equids (Supplementary Table 1). Fresh fecal samples of PHs (n = 12) and AWAs (n = 10) in the wild environment were collected on 21–22 September 2020 at Kalamaili Nature Reserve (KNR), located in the Altay Prefecture of Xinjiang Uygur Autonomous Region, China (44°36′-46°00′ N, 88°30′-90°03′ E). Fresh fecal samples of PHs (n = 13) and AWAs (n = 6) in the captive environment were collected on 26 September 2020 at Gansu Endangered Animals Protection Center (GEAPC) in Wuwei, Gansu Province, China (37°52′50″ N, 102°52′48″ E). The KNR is geographically located in the arid desert region of the Junggar Basin, covering an area of 18,000 km2, where PHs and AWAs roam freely, mainly feeding on bunchgrasses like needlegrass (Stipa capillata) and subshrubs such as Pamirian winterfat (Krascheninnikovia ceratoides), wormwood (Artemisia spp.), and Anabasis brevifolia that distributed in reserve (Meng, 2007). The GEAPC lies 1,400 km from KNR, at the southern edge of the Tengger Desert, where PHs and AWAs are kept in enclosures and fed daily with alfalfa (Medicago sativa) hay. Fresh feces samples were collected into sterile centrifuge tubes, labeled, and immediately placed in a portable fridge before being transported to the laboratory and stored at −20°C until DNA extraction.
The animal study was reviewed and approved by the Ethics and Animal Welfare Committee of Beijing Forestry University (EAWC_BJFU_2021012). Fecal sample collection was approved by the WHBRC and GEAPC.
DNA Extraction
Total DNA extraction from the microbial community of feces samples was performed according to the manufacturer’s instructions of the E.Z.N.A.® soil DNA kit (Omega Bio-tek, Norcross, GA, United States). The TissueLyser (Onebio tech, China) was used to add a mechanical disruption stage of bacterial cells by bead-beating using disruptor tubes (Omega Bio-tech, United States) at a frequency of 45 Hz for 250 s. The quality of extracted DNA was examined using 1% agarose gel electrophoresis, and DNA concentration and purity were determined with NanoDrop 2000.
16S rRNA Gene Amplicon Sequencing and Data Processing
PCR amplification of the V3-V4 hyper-variable region of the 16S rRNA gene of all DNA samples was carried out using 338F and 806R (Mori et al., 2014) with the following amplification procedure: Initial denaturation at 95°C for 3 min, 27 cycles (denaturation at 95°C for 30 s, annealing at 55°C for 30 s, elongation at 72°C for 45 s), followed by extension at 72°C for 10 min, and finally storage at 4°C. The PCR mixtures consist of template DNA 10 ng, 0.8 μl of 5 μM of each forward and reverse primer, 2 μl of 2.5 mM dNTPs, 0.4 μl of TransStart FastPfu DNA polymerase, 4 μl of 5 × TransStart FastPfu buffer, and finally add ddH2O supplement to 20 μl. PCR reaction of each sample was performed in triplicate. The PCR products of each sample were mixed and extracted from 2% agarose gel, then purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, United States) according to manufacturer’s instructions, detected by 2% agarose gel electrophoresis, and quantified using the QuantiFluor® dsDNA System (Promega, United States). Purified amplicons were pooled in equimolar and paired-end sequenced on an Illumina Miseq PE300 platform. 16S rRNA gene amplicon sequencing data were deposited into the NCBI Sequence Read Archive (SRA) database under accession numbers PRJNA784453 (AWAs and captive PHs) and PRJNA701711 (wild PHs).
The raw 16S rRNA gene sequencing reads were quality-filtered by fastp (version 0.20.0; Chen et al., 2018) and merged by FLASH (version 1.2.11; Magoč and Salzberg, 2011) with a minimum overlap length of 10 bp. Operational taxonomic units (OTUs) were clustered with a 97% similarity criterion using UPARSE (version 7.1; Edgar, 2013), and chimeric sequences were identified and eliminated. The taxonomy of each OTU representative sequence was annotated by RDP Classifier (version 2.11; Wang et al., 2007) against the Silva 16S rRNA database (v138) with a confidence threshold of 70%.
Shotgun Metagenomic Sequencing and Bioinformatics Analysis
Shotgun metagenomic sequencing was carried out on the same DNA extracts as the 16S rRNA analysis. Twelve samples (three each for wild PHs, wild AWAs, captive PHs, and captive AWAs) with high DNA quality were selected. Metagenomic DNA was broken into about 400 bp fragments with Covaris M220 (Gene Company Limited, China). A paired-end library was then constructed using NEXTFLEX Rapid DNA-Seq (Bioo Scientific, Austin, TX, United States). Paired-end sequencing was performed using the NovaSeq Reagent Kits on the NovaSeq 6,000 platform (Illumina Inc., San Diego, CA, United States). Metagenomic sequence data for this project have been deposited in the NCBI Sequence Read Archive (Accession Number: PRJNA784453).
The raw metagenomic data were quality-filtered using the Fastp (version 0.20.0), whereby the adapters were eliminated and reads with low quality (Q < 20) and shorter than 20 bp were removed. Reads that mapped to the PH and AWA genome were also filtered out using BWA (version 0.7.17; Li and Durbin, 2009). Metagenomic data were assembled using MEGAHIT (version 1.1.2; Li et al., 2015) with the minimum contig length set to 300 bp. A non-redundant gene catalog was constructed with CD-HIT (version 4.6.1; Fu et al., 2012) with a 90% sequence identity (90% coverage). Bacterial and archaeal genes were filtered from the non-redundant gene catalog to create a prokaryotic gene catalog. High-quality reads were aligned against the prokaryotic gene catalog with 95% identity using SOAPaligner (version 2.21; Li et al., 2009) to determine the gene abundance in each sample. Representative sequences were mapped to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database using Diamond (version 0.8.35; Buchfink et al., 2014) for functional annotations with an e-value of 1e-5.
Statistical Analysis
16S rRNA gene amplicon sequences were rarefied at a minimum sequencing depth before further analysis and obtained the practical OTU set. Alpha diversity of each sample was calculated using mothur (version 1.30.2; Schloss et al., 2009), and R presented Shannon rarefaction curves. The Wilcoxon rank-sum test was utilized to detect a statistically significant difference in Alpha diversity between groups. R was used to count the bacterial abundance of each sample at each taxonomic level and to visualize the bacterial community composition. The QIIME pipeline (Caporaso et al., 2010) was applied to calculate the beta diversity metrics, and NMDS analysis was performed using the R Vegan package. The differences between groups were tested via ANOSIM analysis with a Monte Carlo permutation test (999 permutations). Linear discriminant analysis (LDA) effect size (LEfSe) analyses were performed using LEfSe software1 (Segata et al., 2011) to compare the microbial composition differences between wild and captive equids and to compare the metagenomic functional differences based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database.
Results
Sequence Statistics
A total of 2,064,221 raw reads were obtained after 16S rRNA gene sequencing of 41 fecal samples from twelve wild PHs, ten captive PHs, thirteen wild AWAs, and six captive AWAs, with an average of 50,347 ± 8,304 reads per sample, and there were 1,007,730 reads (24,580 reads per sample) after quality control and sequences rarefied. A total of 3,660 OTUs were obtained after a 97% clustering threshold and classified into 25 phyla, 57 classes, 142 orders, 258 families, and 561 genera of bacteria. Shotgun metagenomic sequencing produced a total of 892,958,226 reads from 12 samples, with an average of 74,413,186 ± 4,493,192 reads per sample. In total, 5,093,229 catalog genes were constructed with an average sequence length of 539 bp in the non-redundant gene catalog.
Microbial Diversity Analysis
There were varying levels of Alpha diversity exhibited across PHs and AWAs in different living environments. The Shannon rarefaction curves revealed that each sample’s Shannon diversity increased and eventually leveled off with increasing sequencing depth, indicating the amount of sequencing data is adequate; nevertheless, the Shannon diversity of captive PHs and captive AWAs presented lower levels (Figure 1). The Wilcoxon rank-sum test was utilized to determine the richness (Chao index), diversity (Shannon index), and evenness (Shannon even index) of the gut microbial communities of wild and captive equids. There were significant differences in the Alpha diversity of equid gut microbiota between the two living environments. Except for microbial richness between AWAs, practically all diversity indexes exhibited that the wild population was significantly higher than the captive population. However, there was no significant difference in the microflora Alpha diversity between the two species of Equidae under the same living environment (i.e., wild or captive; Figure 2). These results indicate that the same survival environment confers convergent gut microbial Alpha diversity to PHs and AWAs, although the wild environment manifests a higher contribution in maintaining high Alpha diversity of gut microbiota.
Figure 1. Shannon index’s Alpha diversity rarefaction curves on the OTU level. The abscissa represents the number of reads selected randomly; the ordinate represents the Alpha diversity measured.
Figure 2. Wilcoxon rank-sum test for differences in Alpha diversity (Shannon index) between groups on the OTU level. **p < 0.01 and ***p < 0.001.
The Beta diversity of the two species of Equidae was calculated at the OTU level by the NMDS analysis based on weighted UniFrac distance (stress: 0.085). The results displayed that the microbiota composition of cohabiting PHs and AWAs was clustered, whereas the microflora structure of equids living in different environments had apparent separation (Figure 3).
Figure 3. Non-metric multidimensional scaling (NMDS) plot of the gut microbiota in equids on the OTU level. The closer the points, the more similar of microbial community composition.
Microbial Composition Analysis
The analysis of microbial community composition demonstrated that PHs and AWAs inhabiting the same environment had convergent microbiota structures, whereas the microflora of equids in different environments had significant differences. At the Phylum level, the gut microbiota of PHs and AWAs in the wild was primarily composed of Firmicutes (PH: 61.14%; AWA: 57.62%), Verrucomicrobiota (PH: 17.73%; AWA: 16.01%), and Bacteroidota (PH: 10.26%; AWA: 7.80%), while the microflora of equids in the captivity was principally constituted of Firmicutes (PH: 66.26%; AWA: 64.74%) and Actinobacteriota (PH: 24.03%; AWA: 27.66%). Meanwhile, equids in the captive environment had a higher ratio of Firmicutes to Bacteroidota (F/B; Figure 4A). At the Family level, the microbial composition of equids in the wild exhibited higher community evenness, with Lachnospiraceae, Christensenellaceae, norank_o__WCHB1-41, Oscillospiraceae, and Akkermansiaceae all accounting for considerable percentages. The symbiotic microbiota of equids in captivity, on the other hand, presented a fair preponderance of Planococcaceae and Micrococcaceae (Figure 4B).
Figure 4. Fecal microbial community composition on the Phylum (A) and family (B) level in equids in the wild and captivity and values are averaged within groups.
Linear discriminant analysis (LDA) effect size (LEfSe) analyses with the same LDA score threshold were performed on the microbial communities of PHs and AWAs in the wild and captivity, respectively, to explore the effects of environmental change on the relative abundance of gut microbiota in equids (Figure 5). The results exhibited that in comparison with equids living in different environments, different species of Equidae living in the same environment have fairly similar gut microbiota structures. In detail, only 7 (1 phylum, 1 class, 1 order, 2 families, 2 genera) and 12 taxa (1 phylum, 0 class, 2 orders, 4 families, 5 genera) displayed significant differences in relative abundance in the same living environment (wild or captive, respectively) of the gut microbiota between PHs and AWAs. However, the relative abundance of 51 (3 phyla, 7 classes, 12 orders, 14 families, 15 genera) and 62 taxa (5 phyla, 9 classes, 13 orders, 17 families, 18 genera) of the gut microbiota of PHs and AWAs, respectively, exhibited significant differences between different environments. Specifically, at the family level, Lachnospiraceae, Christensenellaceae, norank_o__WCHB1-41, Oscillospiraceae, Akkermansiaceae, Rikenellaceae, Eubacterium_coprostanoligenes_group, and Anaerovoracaceae exhibited higher relative abundance in wild PHs and AWAs. Meanwhile, Planococcaceae, Micrococcaceae, Nocardiaceae, and Corynebacteriaceae showed higher relative abundance in captive populations. On the side, Carnobacteriaceae and Enterococcaceae were abundant in captive PHs, whereas Dietziaceae and Moraxellaceae had a significantly higher trend in captive AWAs. These results suggested that the effect of environmental variation from wild to captive on the composition and abundance of microflora in PHs and AWAs were remarkably comparable. In other words, cohabitation converged the establishment of the gut microbiota in different Equidae species.
Figure 5. Linear discriminant analysis (LDA) effect size (LEfSe) analyses (LDA score threshold >4.0) on fecal microbial communities between wild PHs and AWAs (A), captive PHs and AWAs (B), wild and captive PHs (C), and wild and captive AWAs (D).
Microbial Functional Analysis
Although PHs and AWAs living in the same environment have remarkably similar gut microbiota structures, equids in different environments have significantly different microflora. Therefore, based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, the gut microbiome function (KEGG pathway level 2) of equids living in the wild and captivity was analyzed by the linear discriminant analysis (LDA) effect size (LEfSe) algorithm to investigate how this variation in microbiota structure would affect the functional pathway (Figure 6). The results indicated that the abundance of pathways for digestion and utilization of nutrients, including “Digestive system,” “Amino acid metabolism,” “Metabolism of cofactors and vitamins,” and “Lipid metabolism,” and pathways related to sub-health and disease, including “Xenobiotics biodegradation and metabolism,” “Aging,” “Substance dependence,” “Infectious disease: parasitic,” and “Immune disease,” were higher in captive PHs and AWAs than in wild populations. Wild equids, however, had a higher abundance of pathways for metabolism and development, including “Glycan biosynthesis and metabolism,” “Nucleotide metabolism,” “Cell growth and death,” “Translation,” and “Transcription,” as well as pathways about handling various adverse factors such as “Replication and repair,” “Immune system,” and “Environmental adaptation,” than captive populations. This analysis revealed that equids in different environments exhibited apparent divergence in microbiome functions in addition to the diversity and composition of the gut microbiota.
Figure 6. Linear discriminant analysis (LDA) effect size (LEfSe) analyses (LDA score threshold >2.0) of KEGG functional pathways for the microbial metagenome of the wild and captive equids.
Discussion
In addition to host phylogeny, environmental factors such as diet, habitat, and lifestyle play a significant role in establishing gut microflora in wildlife (Trevelline et al., 2019; Yu et al., 2022). Meanwhile, the considerable environmental changes that captivity imposes on wildlife remarkably impact their symbiotic gut microbiota (Chong et al., 2019). In this study, we compared the diversity, composition, and function of the gut microbiota of two Equidae species (the Przewalski’s Horse; PH, and the Asian wild ass; AWA) living in different environmental conditions (wild or captive) by 16S rRNA gene sequencing, to explore how the same living environment shapes the gut microflora of different species of Equidae and to reveal the implications of changes in survival environments on equine health from the perspective of the symbiotic microbiome.
The Alpha and Beta diversity of the microflora of PHs and AWAs in the same environment exhibited remarkable similarity, indicating that the cohabitation had a convergence effect on the gut microbiota diversity in different species of Equidae. Equids in captivity had significantly lower Alpha diversity, which is typical with most captive wildlife research — captive animals’ symbiotic microbiome diversity has a considerably decreased tendency than their wild counterparts (Kohl et al., 2014; Wang et al., 2017; Quiroga-González et al., 2021; San Juan et al., 2021). The decline in diversity might be due to the artificial captivity of wildlife greatly hinders their exposure to microorganisms in the natural environment and increases the risk of exposure to antibiotics and anthelmintics, which have detrimental effects on developing the diversity of the microflora (McKenzie et al., 2017; Kunz et al., 2019). At the same time, the limited food variety in the captive lifestyle may also lead to the loss of Alpha diversity in the gut microbiota (McKenzie et al., 2017; Reese and Dunn, 2018). However, low Alpha diversity of the symbiotic microbiome tends to affect the host’s regular metabolism and reduce the host’s ability to resist adverse environments, thus negatively affecting the health of wildlife (Li et al., 2018; Xing et al., 2019). A previous comparative study of captive and wild Tibetan wild asses showed that captivity reduced the Alpha diversity of their gut microbiota and increased their risk of disease (Gao et al., 2019).
The dramatic differences in the bacterial Phylum level between wild and captive equids suggested that various environments had distinctive shaping effects on the host symbiotic microbiota. In the meantime, the captive PHs and AWAs had a considerably higher ratio of Firmicutes/Bacteroidota (F/B) in their gut microbiota and apparent enrichment of Actinobacteriota compared to wild populations. However, many studies have found that as a critical “obesity marker,” a high level of F/B ratio is frequently associated with obesity in humans or animals (Ley et al., 2005; Stojanov et al., 2020). Meanwhile, the relative abundance of Actinobacteriota in the intestinal tract is also positively correlated with obesity (Turnbaugh et al., 2009). With the constant supply of “free” food and the severe drop in exercise, obesity has become a prevalent and intractable problem in captive wildlife, which considerably impacts animals’ regular metabolism and health (D’Eath et al., 2009). Therefore, a larger ratio of F/B and a higher relative abundance of Actinobacteriota may indicate a greater risk of obesity in captive equids, as well as a threat to the health and conservation of captive threatened equids. It suggests that the F/B ratio of gut microbiota can be utilized as a selective indicator for identifying candidate individuals in the reintroduction of endangered equids and that gut microflora monitoring is required as part of the preliminary preparation for the wildlife reintroduction.
The result of LEfSe analysis of the gut microbial taxa of PHs and AWAs in the two environments showed that the composition and relative abundance of the gut microbiota of distinct species of sympatric equids were essentially the same; however, the gut microflora has dramatically changed while living in different environments, even in the same species. Meanwhile, it is noteworthy that most taxa with significant differences in relative abundance under different environmental conditions play essential roles in host health and disease. Some probiotics that play a critical role in the degradation and fermentation of plant materials were more abundant in wild populations. Lachnospiraceae and Rikenellaceae, for example, can transform polysaccharides such as plant fibers that animals cannot digest into short-chain fatty acids (SCFAs) that can be absorbed and utilized by the host, mainly acetic, propionic, and butyric acids, which serve as the host’s principal energy source (Boutard et al., 2014; Graf, 2014; Koh et al., 2016). Meanwhile, wild equids showed a higher relative abundance of Akkermansiaceae and Norank_o__WCHB1-41, and the enterotype represented by these two bacteria can effectively improve the energy and nitrogen utilization efficiency of wildlife (Guo et al., 2021). Due to animals in the field being subject to poorer food conditions, both in quantity and quality, more bacteria associated with energy harvesting and utilization may provide them with a greater capacity to adapt to the harsh wild environment (Yan et al., 2021). In addition, Christensenellaceae and Oscillospiraceae, more abundant in wild equids, are strongly negatively correlated with metabolic diseases such as inflammation and are essential anti-inflammatory beneficial bacteria (Gophna et al., 2017; Waters and Ley, 2019). Notably, wild equids had substantially more Eubacterium_coprostanoligenes_group, which was negatively correlated with anxiety index and positively correlated with anxiety reduction. (Chen et al., 2019).
On the contrary, Enterococcaceae is a potential pathogenic pathogen, of which the genus Enterococcus is a biomarker for psychological diseases such as depression (Nikolova et al., 2021). However, the relative abundance of Enterococcus showed a significant increasing trend in captive PHs (Figure 5C). Previous research has shown that captivity implies more confined living spaces, more anthropogenic stimuli, and less social interaction, thus predisposing wildlife to psychological disorders such as depression and anxiety (Ferdowsian et al., 2011; Lecorps et al., 2021). Thus, variations in the relative abundance of relevant bacterial taxa suggest that the captive environment may impair the mental health of threatened equids to some extent. Planococcaceae were significantly increased and had the highest relative abundance in captive equids, and an animal model study suggests that elevated levels of this bacterium may be associated with a high-calorie diet and may exacerbate local inflammatory responses (Bai et al., 2020). Forage for captive Equidae is mainly alfalfa hay, which, as a high-quality forage grass worldwide used, has a high carbohydrate and protein content (Scholtz et al., 2009; Bouton, 2012), as well as a constant supply of feed, which may cause a large caloric intake of captive Equidae, resulting in a dramatic increase in the level of Planococcaceae. Potentially pathogenic bacteria such as Micrococcaceae (Popovici et al., 1978), Nocardiaceae (Goodfellow, 2014), Carnobacteriaceae (Aranaz et al., 2021), and Dietziaceae (Koerner et al., 2009) were more abundant in PHs, AWAs, or both in the captive environment than in wild populations. These results indicated that PHs and AWAs in the wild possessed more “healthy” gut microbiota, whereas captive populations trended more pro-inflammatory and pathogenic bacteria, as well as bacteria associated with psychological disorders.
Microbial functional analysis revealed that equids in captivity exhibited significant enrichment in nutrient metabolisms such as amino acid, lipid, cofactors, and vitamins metabolism, implying again that the food of captive equids is of higher nutritional value and more abundant food resources. Gut microbes play an essential role in addiction formation (Qin et al., 2021). The increase in the “substance dependence” pathway in captive populations may be related to gut microbiota disturbance caused by excess antibiotics and anthelmintics interventions. Meanwhile, the “Xenobiotics biodegradation and metabolism” pathway displayed an unusually high trend in captive populations, which further proved the exposure risk of captive wildlife to exogenous substances such as drugs, thus increasing the probability of gut microbial dysbiosis and even disease. Higher “Infectious disease: Parasitic” and “Immune disease” pathways in captive equids demonstrated significant subhealth status. Captive populations enriched in the “Aging” pathway, indicating that PHs and AWAs in captivity may have shorter lifespans than wild populations, which may be due to more significant captive stress than captive welfare (Mason, 2010). Equids in the wild presented enrichment in pathways such as “Glycan biosynthesis and metabolism,” “Nucleotide metabolism,” and “Cell growth and death,” indicating higher levels of growth and development (Viinikangas et al., 2021), which might be attributed to the higher levels of short-chain fatty acid (SCFA) generating bacteria which are involved in energy-yielding. In addition, wild equids were more abundant in “Replication and repair,” “Environmental adaptation,” and “Immune system” pathways suggesting that wild populations may have more outstanding capabilities to adapt to complex environments and resist pathogens.
Conclusion
Overall, PHs and AWAs had significantly different microbial diversities, compositions, and predicted functions in the wild and captivity, and all three aspects exhibited considerable convergence in sympatric equids. It suggests that the same living environment could have a convergence effect on the shaping process of the gut microbiota of different species of Equidae to some extent. Furthermore, the gut microflora of wild equids contained more beneficial bacteria associated with energy generation and anti-inflammatory activity, indicating that the evolutionary adaptation of gut microbiota endowed the host with greater adaptability and a healthier physique. However, captive equids showed significantly “unhealthy” microbiota, with lower Alpha diversity, higher levels of pro-inflammatory bacteria, and potentially pathogenic bacteria, including more biomarkers that may represent psychological disorders such as anxiety or depression, as well as enrichment of microbial functions related to exposure to exogenous substances and susceptibility. These results indicate that the captive environment has apparent adverse effects on threatened equids’ physiological and mental health, which may be due to high-intensity drug exposure, artificial stimulation, limited living space, etc. It hints that in the captive breeding of endangered equids, we should pay attention to the suitable combination of coarse and refined grains, avoid excessive human exposure, appropriate use of antibiotics and anthelmintics, and provide ample active space. Our findings demonstrate that environmental factors play a dominant role in establishing the host’s gut microbiota and warrant further consideration and improvement in the future conservation of captive wildlife.
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 reviewed and approved by the Ethics and Animal Welfare Committee of Beijing Forestry University (EAWC_BJFU_2021012).
Author Contributions
DZ, ZZ, DH, and LY conceived and designed the experiments. HJ, YX, and CS performed the sample collection and DNA extraction. ZZ, LT, HW, and LH analyzed the data. ZZ, LT, QZ, JS, and DZ wrote and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was funded by the Beijing Forestry University Outstanding Young Talent Cultivation Project (no. 2019JQ0318).
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 are grateful to rangers and staff members at Kalamaili Nature Reserve and Gansu Endangered Animals Protection Center for their support in sample collection.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb. 2022.832410/full#supplementary-material
Footnotes
References
Angelakis, E., Yasir, M., Bachar, D., Azhar, E. I., Lagier, J. C., Bibi, F., et al. (2016). Gut microbiome and dietary patterns in different Saudi populations and monkeys. Sci. Rep. 6:32191. doi: 10.1038/srep32191
Aranaz, P., Ramos-Lopez, O., Cuevas-Sierra, A., Martinez, J. A., Milagro, F. I., and Riezu-Boj, J. I. (2021). A predictive regression model of the obesity-related inflammatory status based on gut microbiota composition. Int. J. Obes. 45, 2261–2268. doi: 10.1038/s41366-021-00904-4
Bai, C., Liu, T., Xu, J., Ma, X., Huang, L., Liu, S., et al. (2020). Effect of High Calorie Diet on Intestinal Flora in LPS-Induced Pneumonia Rats. Sci. Rep. 10:1701. doi: 10.1038/s41598-020-58632-0
Blyton, M. D. J., Soo, R. M., Whisson, D., Marsh, K. J., Pascoe, J., Le Pla, M., et al. (2019). Faecal inoculations alter the gastrointestinal microbiome and allow dietary expansion in a wild specialist herbivore, the koala. Anim. Microbiome 1:6. doi: 10.1186/s42523-019-0008-0
Boutard, M., Cerisy, T., Nogue, P. Y., Alberti, A., Weissenbach, J., Salanoubat, M., et al. (2014). Functional Diversity of Carbohydrate-Active Enzymes Enabling a Bacterium to Ferment Plant Biomass. PLoS Genet. 10:e1004773. doi: 10.1371/journal.pgen.1004773
Bouton, J. H. (2012). An overview of the role of lucerne (Medicago sativa L.) in pastoral agriculture. Crop Pasture Sci. 63, 734–738. doi: 10.1071/CP12127
Bowkett, A. E. (2009). Recent captive-breeding proposals and the return of the ark concept to global species conservation. Conserv. Biol. 23, 773–776. doi: 10.1111/j.1523-1739.2008.01157.x
Buchfink, B., Xie, C., and Huson, D. H. (2014). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Chen, Y. H., Bai, J., Wu, D., Yu, S. F., Qiang, X. L., Bai, H., et al. (2019). Association between fecal microbiota and generalized anxiety disorder: Severity and early treatment response. J. Affect. Disord. 259, 56–66. doi: 10.1016/j.jad.2019.08.014
Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). Fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890. doi: 10.1093/bioinformatics/bty560
Chong, R., Grueber, C. E., Fox, S., Wise, P., Barrs, V. R., Hogg, C. J., et al. (2019). Looking like the locals - gut microbiome changes post-release in an endangered species. Anim. Microbiome 1:8. doi: 10.1186/s42523-019-0012-4
Clayton, J. B., Vangay, P., Huang, H., Ward, T., Hillmann, B. M., Al-Ghalith, G. A., et al. (2016). Captivity humanizes the primate microbiome. Proc. Natl. Acad. Sci. U. S. A. 113, 10376–10381. doi: 10.1073/pnas.1521835113
Cox, A. J., West, N. P., and Cripps, A. W. (2015). Obesity, inflammation, and the gut microbiota. Lancet Diabetes Endocrinol. 3, 207–215. doi: 10.1016/S2213-8587(14)70134-2
D’Eath, R. B., Tolkamp, B. J., Kyriazakis, I., and Lawrence, A. B. (2009). “Freedom from hunger” and preventing obesity: the animal welfare implications of reducing food quantity or quality. Anim. Behav. 77, 275–288. doi: 10.1016/j.anbehav.2008.10.028
Edgar, R. C. (2013). UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Ferdowsian, H. R., Durham, D. L., Kimwele, C., Kranendonk, G., Otali, E., Akugizibwe, T., et al. (2011). Signs of mood and anxiety disorders in chimpanzees. PLoS One 6:e19855. doi: 10.1371/journal.pone.0019855
Fu, L., Niu, B., Zhu, Z., Wu, S., and Li, W. (2012). CD-HIT: Accelerated for clustering the next-generation sequencing data. Bioinformatics 28, 3150–3152. doi: 10.1093/bioinformatics/bts565
Gao, H., Chi, X., Qin, W., Wang, L., Song, P., Cai, Z., et al. (2019). Comparison of the gut microbiota composition between the wild and captive Tibetan wild ass (Equus kiang). J. Appl. Microbiol. 126, 1869–1878. doi: 10.1111/jam.14240
Gibson, K. M., Nguyen, B. N., Neumann, L. M., Miller, M., Buss, P., Daniels, S., et al. (2019). Gut microbiome differences between wild and captive black rhinoceros – implications for rhino health. Sci. Rep. 9:7570. doi: 10.1038/s41598-019-43875-3
Goodfellow, M. (2014). “The family Nocardiaceae,” in The Prokaryotes: Actinobacteria. eds. E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin: Springer-Verlag Berlin Heidelberg), 595–650.
Gophna, U., Konikoff, T., and Nielsen, H. B. (2017). Oscillospira and related bacteria – From metagenomic species to metabolic features. Environ. Microbiol. 19, 835–841. doi: 10.1111/1462-2920.13658
Graf, J. (2014). “The family rikenellaceae,” in The Prokaryotes: Other Major Lineages of Bacteria and the Archaea. eds. E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin, Heidelberg: Springer), 857–859.
Greene, L. K., Clayton, J. B., Rothman, R. S., Semel, B. P., Semel, M. A., Gillespie, T. R., et al. (2019). Local habitat, not phylogenetic relatedness, predicts gut microbiota better within folivorous than frugivorous lemur lineages. Biol. Lett. 15, 5–11. doi: 10.1098/rsbl.2019.0028
Groussin, M., Mazel, F., Sanders, J. G., Smillie, C. S., Lavergne, S., Thuiller, W., et al. (2017). Unraveling the processes shaping mammalian gut microbiomes over evolutionary time. Nat. Commun. 8:14319. doi: 10.1038/ncomms14319
Guo, N., Wu, Q., Shi, F., Niu, J., Zhang, T., Degen, A. A., et al. (2021). Seasonal dynamics of diet–gut microbiota interaction in adaptation of yaks to life at high altitude. npj Biofilms Microbiomes 7:38. doi: 10.1038/s41522-021-00207-6
Hanning, I., and Diaz-Sanchez, S. (2015). The functionality of the gastrointestinal microbiome in non-human animals. Microbiome 3:51. doi: 10.1186/s40168-015-0113-6
Jin, L., Huang, Y., Yang, S., Wu, D., Li, C., Deng, W., et al. (2021). Diet, habitat environment and lifestyle conversion affect the gut microbiomes of giant pandas. Sci. Total Environ. 770:145316. doi: 10.1016/j.scitotenv.2021.145316
Kartzinel, T. R., Hsing, J. C., Musili, P. M., Brown, B. R. P., and Pringle, R. M. (2019). Covariation of diet and gut microbiome in African megafauna. Proc. Natl. Acad. Sci. U. S. A. 116, 23588–23593. doi: 10.1073/pnas.1905666116
Koerner, R. J., Goodfellow, M., and Jones, A. L. (2009). The genus Dietzia: a new home for some known and emerging opportunist pathogens. FEMS Immunol. Med. Microbiol. 55, 296–305. doi: 10.1111/j.1574-695X.2008.00513.x
Koh, A., De Vadder, F., Kovatcheva-Datchary, P., and Bäckhed, F. (2016). From dietary fiber to host physiology: Short-chain fatty acids as key bacterial metabolites. Cell 165, 1332–1345. doi: 10.1016/j.cell.2016.05.041
Kohl, K. D., Skopec, M. M., and Dearing, M. D. (2014). Captivity results in disparate loss of gut microbial diversity in closely related hosts. Conserv. Physiol. 2:cou009. doi: 10.1093/conphys/cou009
Kunz, I. G. Z., Reed, K. J., Metcalf, J. L., Hassel, D. M., Coleman, R. J., Hess, T. M., et al. (2019). Equine Fecal Microbiota Changes Associated With Anthelmintic Administration. J. Equine Vet. Sci. 77, 98–106. doi: 10.1016/j.jevs.2019.01.018
Lecorps, B., Weary, D. M., and von Keyserlingk, M. A. G. (2021). Captivity-Induced Depression in Animals. Trends Cogn. Sci. 25, 539–541. doi: 10.1016/j.tics.2021.03.010
Ley, R. E., Bäckhed, F., Turnbaugh, P., Lozupone, C. A., Knight, R. D., and Gordon, J. I. (2005). Obesity alters gut microbial ecology. Proc. Natl. Acad. Sci. U. S. A. 102, 11070–11075. doi: 10.1073/pnas.0504978102
Ley, R. E., Hamady, M., Lozupone, C., Turnbaugh, P. J., Ramey, R. R., Bircher, J. S., et al. (2008). Evolution of mammals and their gut microbes. Science 320, 1647–1651. doi: 10.1126/science.1155725
Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760. doi: 10.1093/bioinformatics/btp324
Li, Y., Hu, X., Yang, S., Zhou, J., Qi, L., Sun, X., et al. (2018). Comparison between the fecal bacterial microbiota of healthy and diarrheic captive musk deer. Front. Microbiol. 9:300. doi: 10.3389/fmicb.2018.00300
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
Li, R., Yu, C., Li, Y., Lam, T. W., Yiu, S. M., Kristiansen, K., et al. (2009). SOAP2: An improved ultrafast tool for short read alignment. Bioinformatics 25, 1966–1967. doi: 10.1093/bioinformatics/btp336
Li, Y., Zhang, T., Shi, M., Zhang, B., Hu, X., Xu, S., et al. (2021). Characterization of intestinal microbiota and fecal cortisol, T3, and IgA in forest musk deer (Moschus berezovskii) from birth to weaning. Integr. Zool. 16, 300–312. doi: 10.1111/1749-4877.12522
Magoč, T., and Salzberg, S. L. (2011). FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Mason, G. J. (2010). Species differences in responses to captivity: Stress, welfare and the comparative method. Trends Ecol. Evol. 25, 713–721. doi: 10.1016/j.tree.2010.08.011
McKenzie, V. J., Song, S. J., Delsuc, F., Prest, T. L., Oliverio, A. M., Korpita, T. M., et al. (2017). The effects of captivity on the mammalian gut microbiome. Integr. Comp. Biol. 57, 690–704. doi: 10.1093/icb/icx090
Meng, Y. (2007). Research on Feeding Plants, Feeding Habits and Feeding Strategies of Przewalski Horse. Beijing: Beijing Forestry University.
Moeller, A. H., Peeters, M., Ndjango, J. B., Li, Y., Hahn, B. H., and Ochman, H. (2013). Sympatric chimpanzees and gorillas harbor convergent gut microbial communities. Genome Res. 23, 1715–1720. doi: 10.1101/gr.154773.113
Mori, H., Maruyama, F., Kato, H., Toyoda, A., Dozono, A., Ohtsubo, Y., et al. (2014). Design and experimental application of a novel non-degenerate universal primer set that amplifies prokaryotic 16S rRNA genes with a low possibility to amplify eukaryotic rRNA genes. DNA Res. 21, 217–227. doi: 10.1093/dnares/dst052
Muegge, B. D., Kuczynski, J., Knights, D., Clemente, J. C., González, A., Fontana, L., et al. (2011). Diet drives convergence in gut microbiome functions across mammalian phylogeny and within humans. Science 332, 970–974. doi: 10.1126/science.1198719
Nikolova, V. L., Smith, M. R. B., Hall, L. J., Cleare, A. J., Stone, J. M., and Young, A. H. (2021). Perturbations in Gut Microbiota Composition in Psychiatric Disorders: A Review and Meta-analysis. JAMA Psychiat. 78, 1343–1354. doi: 10.1001/jamapsychiatry.2021.2573
Peixoto, R. S., Harkins, D. M., and Nelson, K. E. (2021). Advances in Microbiome Research for Animal Health. Annu. Rev. Anim. Biosci. 9, 289–311. doi: 10.1146/annurev-animal-091020-075907
Popovici, M., Greceanu, V., and Vianu, I. (1978). Value of different diagnostic tests for Micrococcaceae and their implication in pathology. Rev. Ig. Bacteriol. Virusol. Parazitol. Epidemiol. Pneumoftiziol. Bacteriol. Virusol. Parazitol. Epidemiol. 23, 231–236.
Qin, C., Hu, J., Wan, Y., Cai, M., Wang, Z., Peng, Z., et al. (2021). Narrative review on potential role of gut microbiota in certain substance addiction. Prog. Neuro-Psychopharmacology Biol. Psychiatry 106:110093. doi: 10.1016/j.pnpbp.2020.110093
Quiroga-González, C., Cardenas, L. A. C., Ramírez, M., Reyes, A., González, C., and Stevenson, P. R. (2021). Monitoring the variation in the gut microbiota of captive woolly monkeys related to changes in diet during a reintroduction process. Sci. Rep. 11:6522. doi: 10.1038/s41598-021-85990-0
Reese, A. T., and Dunn, R. R. (2018). Drivers of microbiome biodiversity: A review of general rules, feces, and ignorance. MBio 9, e01294–e01218. doi: 10.1128/mBio.01294-18
Ren, T., Boutin, S., Humphries, M. M., Dantzer, B., Gorrell, J. C., Coltman, D. W., et al. (2017). Seasonal, spatial, and maternal effects on gut microbiome in wild red squirrels. Microbiome 5:163. doi: 10.1186/s40168-017-0382-3
Roggenbuck, M., Bærholm Schnell, I., Blom, N., Bælum, J., Bertelsen, M. F., Pontén, T. S., et al. (2014). The microbiome of New World vultures. Nat. Commun. 5:5498. doi: 10.1038/ncomms6498
San Juan, P. A., Castro, I., and Dhami, M. K. (2021). Captivity reduces diversity and shifts composition of the Brown Kiwi microbiome. Anim. Microbiome 3:48. doi: 10.1186/s42523-021-00109-0
Sarkar, A., Harty, S., Lehto, S. M., Moeller, A. H., Dinan, T. G., Dunbar, R. I. M., et al. (2018). The Microbiome in Psychology and Cognitive Neuroscience. Trends Cogn. Sci. 22, 611–636. doi: 10.1016/j.tics.2018.04.006
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Schmidt, E., Mykytczuk, N., and Schulte-Hostedde, A. I. (2019). Effects of the captive and wild environment on diversity of the gut microbiome of deer mice (Peromyscus maniculatus). ISME J. 13, 1293–1305. doi: 10.1038/s41396-019-0345-8
Scholtz, G. D. J., van der Merwe, H. J., and Tylutki, T. P. (2009). The nutritive value of South African Medicago sativa L. hay. South African J. Anim. Sci. 39, 179–182. doi: 10.4314/sajas.v39i1.61269
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60
Shreiner, A. B., Kao, J. Y., and Young, V. B. (2015). The gut microbiome in health and in disease. Curr. Opin. Gastroenterol. 31, 69–75. doi: 10.1097/MOG.0000000000000139
Stojanov, S., Berlec, A., and Štrukelj, B. (2020). The influence of probiotics on the firmicutes/bacteroidetes ratio in the treatment of obesity and inflammatory bowel disease. Microorganisms 8, 1–16. doi: 10.3390/microorganisms8111715
Sun, C. H., Liu, H. Y., Liu, B., Yuan, B. D., and Lu, C. H. (2019). Analysis of the Gut Microbiome of Wild and Captive Père David’s Deer. Front. Microbiol. 10:2331. doi: 10.3389/fmicb.2019.02331
Trevelline, B. K., Fontaine, S. S., Hartup, B. K., and Kohl, K. D. (2019). Conservation biology needs a microbial renaissance: A call for the consideration of host-associated microbiota in wildlife management practices. Proc. R. Soc. B Biol. Sci. 286:20182448. doi: 10.1098/rspb.2018.2448
Turnbaugh, P. J., Hamady, M., Yatsunenko, T., Cantarel, B. L., Duncan, A., Ley, R. E., et al. (2009). A core gut microbiome in obese and lean twins. Nature 457, 480–484. doi: 10.1038/nature07540
Viinikangas, T., Khosrowabadi, E., and Kellokumpu, S. (2021). “N-Glycan Biosynthesis: Basic Principles and Factors Affecting Its Outcome,” in Experientia Supplementum (2012). ed. M. Pezer (Cham: Springer International Publishing), 237–257.
Vuong, H. E., Yano, J. M., Fung, T. C., and Hsiao, E. Y. (2017). The Microbiome and Host Behavior. Annu. Rev. Neurosci. 40, 21–49. doi: 10.1146/annurev-neuro-072116-031347
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wang, W., Zheng, S., Sharshov, K., Sun, H., Yang, F., Wang, X., et al. (2017). Metagenomic profiling of gut microbial communities in both wild and artificially reared Bar-headed goose (Anser indicus). Microbiology 6:e00429. doi: 10.1002/mbo3.429
Waters, J. L., and Ley, R. E. (2019). The human gut bacteria Christensenellaceae are widespread, heritable, and associated with health. BMC Biol. 17:83. doi: 10.1186/s12915-019-0699-4
West, A. G., Waite, D. W., Deines, P., Bourne, D. G., Digby, A., McKenzie, V. J., et al. (2019). The microbiome in threatened species conservation. Biol. Conserv. 229, 85–98. doi: 10.1016/j.biocon.2018.11.016
Xing, Y., Liu, J., Lu, F., Wang, L., Li, Y., and Ouyang, C. (2019). Dynamic distribution of gallbladder microbiota in rabbit at different ages and health states. PLoS One 14:2509. doi: 10.1371/journal.pone.0211828
Yan, L., Tang, L., Zhou, Z., Wei, L. U., Wang, B., Sun, Z., et al. (2021). Metagenomics reveals contrasting energy utilization efficiencies of captive and wild camels (Camelus ferus). Integr. Zool. 17, 333–345. doi: 10.1111/1749-4877.12585
Yu, Q., Li, G., and Li, H. (2022). Two community types occur in gut microbiota of large-sample wild plateau pikas (Ochotona curzoniae). Integr. Zool. 17, 366–378. doi: 10.1111/1749-4877.12575
Zarrinpar, A., Chaix, A., Yooseph, S., and Panda, S. (2014). Diet and Feeding Pattern Affect the Diurnal Dynamics of the Gut Microbiome. Cell Metab. 20, 1006–1017. doi: 10.1016/j.cmet.2014.11.008
Zhang, Z., Xu, D., Wang, L., Hao, J., Wang, J., Zhou, X., et al. (2016). Convergent Evolution of Rumen Microbiomes in High-Altitude Mammals. Curr. Biol. 26, 1873–1879. doi: 10.1016/j.cub.2016.05.012
Keywords: gut microbiota, sympatry, microbiome convergence, wild and captive, equid, health
Citation: Zhou Z, Tang L, Yan L, Jia H, Xiong Y, Shang J, Shao C, Zhang Q, Wang H, He L, Hu D and Zhang D (2022) Wild and Captive Environments Drive the Convergence of Gut Microbiota and Impact Health in Threatened Equids. Front. Microbiol. 13:832410. doi: 10.3389/fmicb.2022.832410
Edited by:
Lifeng Zhu, Nanjing Normal University, ChinaCopyright © 2022 Zhou, Tang, Yan, Jia, Xiong, Shang, Shao, Zhang, Wang, He, Hu and Zhang. 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: Dong Zhang, zhangdong_bjfu@bjfu.edu.cn; ernest8445@163.com