Comparative Analyses of Fecal Microbiota in European Mouflon (Ovis orientalis musimon) and Blue Sheep (Pseudois nayaur) Living at Low or High Altitudes

Abstract

The gut microbiota is a complex and essential system organ that plays an integrative role in balancing key vital functions in the host. Knowledge of the impact of altitude on the gut microbiota of European mouflon (Ovis orientalis musimon) and blue sheep (Pseudois nayaur) is currently limited. In this study, we compared the characteristics of gut microbiota in 5 mouflon at low altitude (K group), 4 mouflon at high altitude (L group), 4 blue sheep at low altitude (M group), and 4 blue sheep at high altitude (N group). The V3–V4 region of the 16S rRNA gene was analyzed using high-throughput sequencing. Analyses based on the operational taxonomic units showed significant changes in the gut microbial communities between groups at different altitudes. At the phylum level, groups at the high altitudes had a higher relative abundance of Firmicutes and a lower relative abundance of Bacteroidetes than those at the low altitudes. A higher Firmicutes:Bacteroidetes ratio is beneficial to animals in terms of the gut microbiota-mediated energy harvest. The relative abundance of Proteobacteria was significantly higher in the gut microbiota of mouflon sheep at high altitudes. At the genus level, the Bacteroides:Prevotella ratio was significantly higher in the low-altitude group (than the high-altitude group) of mouflon sheep and the ratio was significantly higher in the high-altitude group (than the low-altitude group) in blue sheep. In addition, the Ruminococcaceae_UCG-005 related to cellulose and starch digestion was the predominant genus in blue sheep and the relative abundance of the genus was significant higher in the high-altitude group than the low-altitude group of blue sheep (P < 0.01). In conclusion, our results suggested that the gut microbiota of high-altitude groups of sheep had stronger abilities related to energy metabolism and the decomposition of substances, e.g., fiber and cellulose, and that such abilities are associated with high-altitude adaptation.

Introduction

Natural environmental factors can have significant effects on all biological aspects of an animal. Animals are often able to adapt to different environments through physiological or behavioral changes. In many animals, a huge number of microbes reside in the gastrointestinal tract and form a complex gut microbiome. Hosts provide a living environment for the microbial communities of the gut, while gut microbiota contributes to nutrient absorption of the host. The gut microbiota plays an integrative role in balancing key vital functions in the host, such as metabolism, inflammation, and immune function (Round and Mazmanian, 2009; Sekirov et al., 2010; Qin et al., 2012; Tremaroli and Bäckhed, 2012). A number of studies have assessed the bacterial community structure and function of the gut microbiota in herbivores (Castillo et al., 2007; Bergmann et al., 2015; Miller et al., 2016; Avila-Jaime et al., 2018). Coevolutionary modes tend to occur when interactions develop between the gut microbiota and hosts (Maynard et al., 2012). Considering the importance of gut microbiota, the identification of factors that affect animal gut microbial flora may have broad implications in the fields of biology, ecology, veterinary science, and more.

The blue sheep (Pseudois nayaur) is a medium-sized herbivore and a species of Bovidae that includes two subspecies: P. nayaur and P. szechuanensis (Marschall, 1998). Blue sheep are distributed throughout Bhutan, China (Tibetan Plateau and the surrounding mountain regions), northern India and Myanmar, Nepal, and northern Pakistan (Shackleton, 1997; Rabinowitz and Khaing, 1998). In addition, blue sheep inhabit high mountainous areas [>4000 m above sea level (a.s.l.)] with woodland steppe or bare rocks as the main habitat types, that are far from water resources (Liu et al., 2018).

Although blue sheep are listed as a Least Concern (LC, Version 3.1) species by the International Union for Conservation of Nature (Harris, 2014) due to its wide distribution, the population structure of the species is still unclear due to the difficulty in conducting estimates in the steep habitats. The total population size of this species is roughly estimated to be 47000–414000 (Fox et al., 1991; Schaller, 1998; Wang, 1998; Liu et al., 2005). In China, the number of blue sheep in the Tibetan Plateau, the main habitat for blue sheep, were estimated and found to have been dramatically declining from 1991 to 2002 (Harris and Loggers, 2004). As a consequence, blue sheep were listed as a second-grade state protection animal and an Endangered class.

The European mouflon (Ovis orientalis musimon) is a subspecies of O. orientalis, a member of the Bovidae, and is considered to be the ancestor of the modern-day sheep (Valdez and Batten, 1982; Groves, 1993; Rezaei et al., 2010). Due to the drastic decline in the population as a result of hunting, environmental changes, and habitat deterioration, the European mouflon is listed as a Vulnerable species under criterion A2cde by the Valdez (2008). Notably, O. orientalis is listed on CITES Appendix I. The mouflon is a herbivorous animal that feeds on grasses and shrubs and inhabits moderate or very arid grasslands that generally occur below 2500 m a.s.l. (Hoefs, 1985).

Blue sheep and mouflon can adapt to complex and hostile environments. Although the two species have strong adaptive abilities to different conditions, their population size were still small. And previous studies have mainly concentrated on their taxonomy, population, distribution, habitats, and ecology (Yuqun, 1990; Hermans, 1996; Shackleton et al., 1997; Namgail et al., 2009; Tan et al., 2012). However, few studies have been conducted on their molecular biology to determine the taxonomic and biological mechanisms (Zeng et al., 2008; Tan et al., 2012; Barbato et al., 2017; Yang et al., 2019).

“Altitude” is an important factor that can affect and shape gut microbial communities. Specifically, altitude can affect the composition and structure of gut microbiota in mammals (Li K. et al., 2016; Karl et al., 2018; Zhao et al., 2018). Low barometric pressure, low temperature, and high radiation conditions are physiologic stressors that can affect both the survival and reproduction of organisms (Hochachka et al., 1996; Cheviron et al., 2012). In addition, intense evolutionary selection pressures can cause a variety of physiological changes in animals (Storz et al., 2010; Ivy and Scott, 2015; Simonson, 2015). Environmental change can also alter the gut microbiota, studies have shown that short exposure to high-altitude can change the composition of gut microbiota.

However, studies using high-throughput sequencing technology to assess the gut microbiota of Bovidae are relatively insufficient, and studies on the gut microbiota of blue sheep and mouflon sheep are greatly lacking. Therefore, we conducted a gut microbiota analysis using samples of blue sheep and mouflon sheep collected from low- (<100 m a.s.l.) and high- (>2300 m a.s.l.) altitude habitats. In this study, we aimed to identify the effects of altitudinal variation on the gut microbial flora of mouflon and blue sheep to provide greater insights into the gut microbiota community composition and structure of animals that normally inhabit different altitudes. And the research will be benefit to adaptive management and the evaluation of conservation efficacy of the two species.

Materials and Methods

Ethics Statement

This study was carried out and executed in accordance with the recommendations of the Guide to Animal Experiments of the Ministry of Science and Technology (Beijing, China). All experiments involving animals were approved by the Qufu Normal University Institutional Animal Care and Use Committee (Permit Number: QFNU2018-031).

Sample Collection

Samples were collected from the Ji’nan Wild Zoo (Shandong Province, China) and Xining Wildlife Park (Qinghai Province, China), which are positioned at low altitudes (<100 m a.s.l.) and high altitudes (>2000 m a.s.l.), respectively. Detail information of the two study species and the groupings are listed in Table 1. The age of all study animals ranged from 3 to 8 years, and the subjects were healthy, free of diseases, and had not been on any medications for at least 3 months prior to the sampling. Both the species from the different locations were fed roughly the same diets, which were primarily comprised of cotton grass, alfalfa, carrots, and pellet feed. Fresh fecal samples were obtained in the early mornings after the barns were cleaned. Samples were collected within half an hour of the aseptic cleaning process. The fecal samples were frozen to cryogenic temperatures for immediate storage and then frozen to −80°C before further processing and analysis.

DNA Extraction

Genomic DNA extraction was performed using QIAamp^® Stool Mini Kit (Qiagen, Germany) according to manufacturer’s recommendations. A Nanodrop UV-Vis Spectrophotometer (Thermo Fisher scientific) was used to determine the amount of genomic DNA.

16S rRNA Gene PCR and Sequencing

To characterize the resident microbial community, we amplified and analyzed the 16S rRNA gene of the V3–V4 region with the forward primer (CCTAYGGGRBGCASCAG) and the reverse primer (GGACTACNNGGGTATCTAAT). The final volume of the amplification reaction mixture was 50 μl and contained 5 μl microbial genomic DNA, 5 μl of each primer, 25 μl KAPA HiFi HotStart Ready Mix (KAPA Biosystems), and 10 μl ddH₂O. The following PCR conditions were used to acquire the products: 95°C for 10 min, followed by 25 cycles at 94°C for 30 s, annealing at 55°C for 30 s, elongation at 72°C for 30 s, and a final elongation step at 72°C for 5 min. The DNA fragments of the products were ∼410 bp in size and were analyzed using agarose gel electrophoresis to confirm the quantity and quality of the products. The PCR products were then purified using AMPure XP beads according to the manufacturer’s instructions. Before sequencing, we used TruSeq DNA PCR-free Sample Preparation Kit (Illumina, United States) to generate the sequencing library. Samples were sequenced using the Illumina HiSeq 2500 platform (Illumina) with the 250 bp. A pair-ended running mode.

Sequence Processing and Analysis

Raw reads from the Illumina MiSeqPE250 were cleaned up and assembled in the following steps. Barcodes and primer sequences were cut off and raw tags were obtained and performed using FLASH (Magoč and Salzberg, 2011) and QIIME (Version 1.7.0) (Caporaso et al., 2010). The UCHIME algorithm (Edgar et al., 2011) was used to remove the chimeric sequences from the raw tags and effective tags were obtained. We assessed the community composition by classifying the sequences (Uparsev7.0.1001) into operational taxonomic units (OTUs) defined by whether the 16S rRNA sequence similarity was greater than or equal to 97% (Edgar, 2013). Finally, species annotations were performed using the GreenGene Database (DeSantis et al., 2006) with the RDP classifier algorithm (Version 2.2) (Wang et al., 2007).

The alpha-diversity of microbial communities was determined using different indices (Chao1, ACE, Shannon, Simpson, Coverage) and calculated with QIIME (V1.7.0). Furthermore, the bacterial community diversity was analyzed using rarefaction and rank-abundance plots, and displayed using R software. Tukey and Wilcoxon’s tests were used for statistical analysis. Differences among samples were assessed using a Principal Component Analysis (PCA), Principal Coordinates Analysis (PCoA), and non-metric multidimensional scaling (NMDS), and displayed using R software. Distance-based methods, such as the Unweighted Pair Group Method with Arithmetic Mean (UPGMA) have been used to conduct cluster analyses based on the similarity and dissimilarity of bacterial communities among samples. Linear discriminant analysis (LDA) effect size (LEfSe) analysis (Segata et al., 2011) was used to detect statistically significant differences in the species between the study groups.

Next-generation sequencing datasets have been deposited in the NCBI sequence read archive (SRA) with accession number: PRJNA511517.

Results

Overview of the Sequencing Data

A total of 1,377,737 high quality reads were produced from 17 fecal samples and were classified into 2774 OTUs after conducting quality control with 97% similarity. The alpha diversity of the microbial communities was assessed and calculated (Table 2). Rarefaction curves (Figure 1A) of OTU richness were calculated using the vegan library in the R statistical computing language (Oksanen et al., 2013). When the rarefaction curves approached a plateau, this suggested that the number of OTUs was sufficient to reveal the authentic bacterial communities within each sample. Furthermore, ranked abundance curves (Figure 1B) were used to analyze the community diversity, which indicated both the evenness and abundance of species in the samples.

FIGURE 1

Bacteria Composition and Relative Abundance

A total of 33 phyla, 64 classes, 121orders, 209 families, 378 genera, and 108 species were detected in the prokaryotic microbiota communities from all fecal samples from both blue sheep and the mouflon.

At the phylum level, the top five ranked abundance-based phyla in the mouflon sheep were Firmicutes (57.25%), Bacteroidetes (33.66%), Proteobacteria (2.14%), Spirochaetes (1.63%), and Verrucomicrobia (1.57%) at low altitudes (group K); and Firmicutes (62.13%), Bacteroidetes (21.07%), Proteobacteria (9.21%), Fibrobacteres (2.90%), and Spirochaetes (1.37%) at high altitudes (group L).

The top five ranked abundance-based phyla in the mouflon were Firmicutes (57.31%), Bacteroidetes (33.48%), Proteobacteria (2.29%), Spirochaetes (2.14%), and Fibrobacteres (1.30%) at low altitudes (group M); and Firmicutes (61.09%), Bacteroidetes (28.30%), Proteobacteria (2.94%), Fibrobacteres (2.62%), and Spirochaetes (1.89%) at high altitudes (group N).

In summary, the most predominant phyla in all groups (K, L, M, and N) were Firmicutes, Bacteroidetes, and Proteobacteria; the proportion of which comprised more than 92% of the total composition of the samples. It is note worthy that Verrucomicrobia was a unique phylum that occurred in the top five ranked phyla in group K, while the others contained the phylum Fibrobacteres. The rank of the community structure and composition of the top five phyla was completely consistent within mouflon and blue sheep.

At the genus level, the most common genera in all groups mainly contained Ruminococcaceae_UCG-005, Ruminococcaceae_UCG-010, Fibrobacter,Christensenellaceae_R-7_group, Rikenellaceae_RC9_gut_group, and Alistipes; of which, the first four belonged to Firmicutes and the last two genera belonged to Bacteroidetes.

We chose the top ten phyla and genera based on species abundance to generate a histogram, which showed the percentages of relative abundance in each sample or group (Figures 2A,B).

FIGURE 2

According to the relative abundance of genera in the bacterial community compositions, we selected the top ranked 35 genera to generate a clustering heatmap (Figure 3) based on relative abundance values. The values are shown in the heatmap grid with increasing abundance from blue to red. The fecal samples from mouflon and blue sheep at low altitudes (groups K and M) were grouped together while those from high altitudes (groups L and N) were clustered into another clade.

FIGURE 3

Cluster analysis using the unweighted Unifrac distance matrix was conducted following the UPGMA clustering method (Figure 4). The results were basically consistent with those from the clustering heatmap analysis.

FIGURE 4

Analysis of Discrepancies for Between Groups

Turkey and Wilcoxon’s tests were used to calculate the alpha diversity (Shannon and observed species) and beta diversity in the two paired groups (Figure 5). Analysis of similarity (ANOSIM) (Table 3) and Multi-Response Permutation Procedures (MRPP) (Table 4) indicated significant differences in the bacterial communities between fecal samples collected from low and high altitudes. Notably, the differences between the samples collected from the different altitudes were considered highly significant (P < 0.01) in mouflon sheep and statistically significant (P < 0.05) in blue sheep. However, the most interesting finding was that the gut microbiota did not differ significantly (P > 0.05) between the mouflon and blue sheep at the same altitudes.

FIGURE 5

To compare the similarity between samples and groups, a heatmap of the beta-diversity index (Figure 6) was constructed to demonstrate the correlativity based on the dissimilarity coefficient. An NMDS ordination plot (Figure 7) of bacterial taxonomy data that demonstrated the microbial differences in samples and the distance between samples reflected the degree of the discrepancy. In addition, PCAs were used to clarify the similarity of the bacterial communities in the samples collected from the same altitude (Figure 7). The gut microbiota in the mouflon and blue sheep at different altitudes were significantly different, and the bacterial population structures in the samples from the different species at the same altitudes were similar.

FIGURE 6

FIGURE 7

Additionally, we employed a T-test and LDA LEfSe to determine significantly different taxa between groups from different altitudes. We calculated the statistical significance using the relative abundance of the main (top 10 ranked) bacterial phyla in the mouflon (Table 5) and blue sheep (Table 6).

Significant differences were observed at the phylum level between groups K and L of mouflon sheep. The relative abundances of Bacteroidetes and Verrucomicrobia were significantly higher in the low-altitude samples (P = 0.021 and 0.033), and those of Acidobacteria and Chloroflexi were significantly higher in the high-altitude samples (P = 0.042 and 0.041). At the genus level, there were 61 genera that differed significantly between groups K and L, of which 22 of were significantly higher in the low-altitude samples while the others were significantly higher in the high-altitude samples.

Regarding differences in blue sheep at the phylum level, the relative abundance of WCHB1-60 was significantly (P = 0.018) higher in group N. In addition, 34 genera showed significant differences between groups M and N, of which 22 were significantly higher in group M and 12 were significantly higher in group N.

It is remarkable that the relative abundances of Firmicutes, Proteobacteria, and Fibrobacteres were slightly higher in the gut microbiota of samples collected from high altitudes in both studied animals. However, the relative abundances of Bacteroidetes and Spirochaetes were slightly higher in samples from low altitudes.

Furthermore, LEfSe used LDA to estimate the magnitude of the effect of altitude and provide a list of the species that showed significant differences. Seven (between group K and group L) and three (between group M and group N) biomarkers were significantly different in terms of their relative abundance in the mouflon (Figure 8A) and blue sheep (Figure 9A). Circular cladograms show the main taxa that were significantly different in the mouflon (Figure 8B) and blue sheep (Figure 9B).

FIGURE 8

FIGURE 9

Discussion

Considering the importance of gut microbiota, the identification of factors that affected the gut microbial flora of animals may have broad implications for various fields of science, including microbiology, ecology, veterinary science, and others. Mouflon and blue sheep are considered to be primarily grazers, and in the present study we compared the gut microbial community of these two species in two independent altitudinal transects. Organisms face physiological challenges at high altitudes, such as decreases in barometric pressure, oxygen content, and temperatures. These challenges and the variations in responses to them provided unique opportunities to investigate their effects, e.g., different altitudes, on microbial composition and structure. Our results indicated that altitudinal variation was a factor that shaped the composition and structure of gut microbiota in both studied herbivores. Our findings were consistent with those of previous studies that demonstrated that altitudinal variation can affect gut microbial communities in mammals (Li H. et al., 2016; Sun et al., 2016; Zhang et al., 2016; Lan et al., 2017).

The cold and hypobaric hypoxic conditions in high-altitude habitats serve as physiologic stressors that can affect both the survival and reproduction of animals (Hochachka et al., 1996; Cheviron et al., 2012). Such stressors can cause an array of physiological changes in animals (Storz et al., 2010; Ivy and Scott, 2015; Simonson, 2015). For instance, exposure to hypobaric hypoxia can lead to increased inflammation and incidences of illness (Hartmann et al., 2000; Khanna et al., 2018). Convergent patterns have also been observed in multiple of species at high altitudes, which were consistent with the hypoxia treatment conducted in previous laboratory experiments (Suzuki et al., 2018). Moreover, weight loss was common in animals based on long sojourn times at high altitudes (Hamad and Travis, 2006; Westerterp and Kayser, 2006). The overall alpha-diversity and beta-diversity measurements, based on the rarefied OTUs in this study, were higher in the low-altitude groups than in the high-altitude groups, and a higher diversity is generally considered to be a mark of a healthy and resilient microbiota (Petersen and Round, 2014).

The Firmicutes and Bacteroidetes were the prominent phyla in both the mouflon and blue sheep, which is consistent with the findings of other studies regarding mammals (Duncan et al., 2008; Ley et al., 2008; Kong et al., 2014). Our results indicated that the groups at higher elevations had a higher relative abundance of Firmicutes and a lower relative abundance of Bacteroidetes than those at low elevations. However, the relative abundance of Firmicutes and Bacteroidetes of the gut microbiota between the mouflon and blue sheep were similar at the same altitude. Based on the understanding that microbes belonging to phylum Firmicutes contain genes related to energy metabolism and the decomposition of substances, such as fiber and cellulose (Kaakoush, 2015), and because the main functions of microbes that belong to Bacteroidetes are to degrade proteins and carbohydrates (Fernando et al., 2010; Waite and Taylor, 2014), our results suggested that the gut microbiota of high-altitude groups of sheep had stronger abilities related to energy metabolism and the decomposition of substances, e.g., fiber and cellulose, and that such abilities are associated with high-altitude adaptation. In addition, we also found that the groups at high altitudes had an increased Firmicutes:Bacteroidetes ratio for both of the two species. The high ratio of these two phyla is beneficial to animals via a gut microbiota-mediated energy harvest (Li H. et al., 2016; Li K. et al., 2016; Zhang et al., 2016; Lan et al., 2017) that can assist the host to maintain a metabolic balance and core body temperature during environmental temperature challenges (Ley et al., 2006; Turnbaugh et al., 2006; Fernando et al., 2010; Murphy et al., 2010).

At the genus level, Christensenellaceae_R-7_group, Bacteroides, Ruminococcaceae_UCG-013, Prevotellaceae_UCG-003, and Akkermansia were identified as the genera that differed significantly between the low- and high-altitude groups of mouflon sheep, while Ruminococcaceae_UCG-005, Bacteroides, Ruminococcaceae_UCG-013, Prevotellaceae_UCG-003, and Prevotella_1 were the genera that showed predominant significant differences in the gut microbial flora of blue sheep. Bacteroides, Ruminococcaceae_UCG-013, and Prevotellaceae_UCG-003 were identified in the guts of both the study species and their relative abundances were found to be higher in animals from low altitudes.

Regarding the mouflon sheep, the relative abundance of the genus Christensenellaceae_R-7_group (phyla: Firmicutes) was greater in the high-altitude group (P ≤ 0.05). Previous studies in mice, humans, and Tibetan macaques showed that the high abundance of this genus was related to a lean body shape, and some species belonging to family Christensenellaceae can help reduce weight gain (Goodrich et al., 2014; Oki et al., 2016; Zhao et al., 2018). Genus Akkermansia was the most negatively correlated bacterial genus with altitude and was an aerotolerant mucin-degrader that colonizes the mucus layer (Reunanen et al., 2015; Ouwerkerk et al., 2016). In addition, Ruminococcaceae_UCG-005, which is related to cellulose and starch digestion, was the predominant genus in blue sheep and the relative abundance the genus was significant higher in animals at high altitudes than in low altitudes (P < 0.01) (Jindou et al., 2006; Ze et al., 2012). The result suggested that Ruminococcaceae_UCG-005 played an important role in cellulose-degradation in blue sheep at high altitudes. Interestingly, the relative abundance of Ruminococcaceae_UCG-013 was significant higher in both studied animals at low altitudes.

Bacteroides are the most frequently isolated microbes from human infections, and they can cause endogenous infections in the host (Veeranagouda et al., 2014). Prevotella, a genus within the family Prevotellaceae, is a saccharolytic anaerobe that is associated with the consumption of a diet rich in carbohydrates and with the production of short-chain fatty acids (SCFAs) (De Filippo et al., 2010; Wu et al., 2011; Ramakrishna, 2013). Prevotella are known SCFA-producers, and SCFAs, such as acetic acid, can provide energy for the host and mediate blood pressure by acting on SCFA-receptors (Tremaroli and Bäckhed, 2012; den Besten et al., 2013). Moreover, a higher relative abundance of Prevotella can aid in the degradation of plant fibers and decrease the lipolytic and proteolytic fermentation potential (Vieira-Silva et al., 2016). The relative abundance of Prevotella increased dramatically in the mouflon at high altitudes, and a similar pattern was observed in pika (Li H. et al., 2016), some ruminants (Zhang et al., 2016), and humans (Li and Zhao, 2015; Lan et al., 2017). Interestingly, the relative abundance of Prevotella between groups of blue sheep at different altitudes was similar. However, previous studies have demonstrated negative correlations between Prevotella abundance and altitude (Li and Zhao, 2015; Zhao et al., 2018). We speculated that the low relative abundance of Prevotella in the high-altitude group may have been correlated with the short migration duration from the low to high altitudes.

Fecal microbial communities were clustered into enterotypes, determined primarily by levels of Bacteroides and Prevotella (Wu et al., 2011). The Bacteroides:Prevotella ratio, which is known as a dietary and lifestyle biomarker, is correlated with long-term habitual diets and provides a useful intervention target for improving host responses to environmental challenges (Gorvitovskaia et al., 2016; Costea et al., 2018). A lower Bacteroides:Prevotella ratio indicates a healthy, high-fiber, plant-rich diet (Gorvitovskaia et al., 2016). Previous studies revealed that a higher Bacteroides:Prevotella ratio is likely to indicate a high fat, high protein, and low fiber diet in animals (Wu et al., 2011). Our results indicated that the Bacteroides:Prevotella ratio in mouflon sheep was significantly higher in the low-altitude group (12.76) than in the high-altitude group (5.28), but in startling contrast, the Bacteroides:Prevotella ratio was significant higher in the high-altitude group (16.03) than in the low-altitude group (2.07) of blue sheep. The results suggested that the mouflon in the high-altitude areas were more effective in digesting high-fiber, plant-rich diets than the mouflon living at low altitudes. The differences in the Bacteroides:Prevotella ratio between the mouflon and blue sheep suggest that these animals may have slight differences in their nutrient digestion and metabolic mechanisms.

Moreover, there were significant differences in the relative abundance of Proteobacteria between the two altitudinal groups of mouflon sheep (P < 0.05). In addition, the relative abundance of Proteobacteria in blue sheep was slightly higher in the high-altitude group than in the low-altitude group, although the difference was not significant (P = 0.568). Proteobacteria was found to be the most predominant phylum in the giant panda and was related to lignin digestion as well as the catabolizing of various components (Evans et al., 2011; Fang et al., 2012). Due to their similar dietary composition to that of sheep, i.e., includes fiber and cellulose and other substances that are not easily digested, we inferred that the higher relative abundance of Proteobacteria observed in the present study would assist the host animals with effectively degrading their food, to obtain higher amounts of nutrients and energy.

Because altitude is an important factor that can affect and shape the gut microbial community, we attempted to understand how altitudinal variations affect the gut microbial flora in mouflon and blue sheep in the present study. Animals in high-altitude areas are exposed to hypobaric, hypoxic, and hypothermic conditions; they show increased physical activity, and even face severe dietary challenges. In our study, we collected fecal samples from study animals (under controlled dietary protocols) using non-invasive methods. Structure and cluster analyses of gut microbiota revealed that groups in the higher altitudes had a higher relative abundance of Firmicutes and a lower relative abundance of Bacteroidetes than those at low altitudes. The high Firmicutes:Bacteroidetes ratio observed in the present study is presumably beneficial to the sheep by producing gut microbiota-mediated energy. The abundance of Proteobacteria was significant higher in the gut microbiota from mouflon sheep living at high altitudes. Compared to the high-altitude groups, the Bacteroides:Prevotella ratio in mouflon sheep was significant higher in the low-altitude group and the ratio in blue sheep was slightly higher in the low-altitude group. The differences between the high and low-altitude groups indicated that high-altitude groups had stronger abilities related to energy metabolism and the decomposition of substances, such as fiber and cellulose. However, our understanding of the importance of altitude for the gut microbial communities of animals was limited because samples were only collected from captive animals. Furthermore, there may potentially be bacterial strains that are not identifiable using the routine procedures of next-generation sequencing technology. Thus, further investigations are needed to comprehensively understand the composition, structure, and function of gut microbial flora of European mouflon and blue sheep.

References

• 1

  Avila-JaimeB.KawasJ.Garcia-MazcorroJ. (2018). Prediction of functional metagenomic composition using archived 16S rDNA sequence data from the gut microbiota of livestock.Livest. Sci.21328–34. 10.1016/j.livsci.2018.04.017

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  BarbatoM.HailerF.Orozco-terWengelP.KijasJ.MereuP.CabrasP.et al (2017). Genomic signatures of adaptive introgression from European mouflon into domestic sheep.Sci. Rep.7:7623. 10.1038/s41598-017-07382-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BergmannG. T.CraineJ. M.RobesonM. S.IIFiererN. (2015). Seasonal shifts in diet and gut microbiota of the American bison (Bison bison).PLoS One10:e0142409. 10.1371/journal.pone.0142409

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  CaporasoJ. G.KuczynskiJ.StombaughJ.BittingerK.BushmanF. D.CostelloE. K.et al (2010). QIIME allows analysis of high-throughput community sequencing data.Nat. Methods7335–336.

  • Google Scholar
  • View reference in article

• 5

  CastilloM.Martín-OrúeS.AnguitaM.PérezJ.GasaJ. (2007). Adaptation of gut microbiota to corn physical structure and different types of dietary fibre.Livest. Sci.109149–152. 10.1016/j.livsci.2007.01.129

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  ChevironZ. A.BachmanG. C.ConnatyA. D.McClellandG. B.StorzJ. F. (2012). Regulatory changes contribute to the adaptive enhancement of thermogenic capacity in high-altitude deer mice.Proc. Natl. Acad. Sci. U.S.A.1098635–8640. 10.1073/pnas.1120523109

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  CosteaP. I.HildebrandF.ArumugamM.BäckhedF.BlaserM. J.BushmanF. D.et al (2018). Enterotypes in the landscape of gut microbial community composition.Nat. Microbiol.38–16. 10.1038/s41564-017-0072-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  De FilippoC.CavalieriD.Di PaolaM.RamazzottiM.PoulletJ. B.MassartS.et al (2010). Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa.Proc. Natl. Acad. Sci. U.S.A.10714691–14696. 10.1073/pnas.1005963107

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  den BestenG.van EunenK.GroenA. K.VenemaK.ReijngoudD.-J.BakkerB. M. (2013). The role of short-chain fatty acids in the interplay between diet, gut microbiota, and host energy metabolism.J. Lipid Res.542325–2340. 10.1194/jlr.R036012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  DeSantisT. Z.HugenholtzP.LarsenN.RojasM.BrodieE. L.KellerK.et al (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB.Appl. Environ. Microbiol.725069–5072. 10.1128/aem.03006-05

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  DuncanS. H.LobleyG.HoltropG.InceJ.JohnstoneA.LouisP.et al (2008). Human colonic microbiota associated with diet, obesity and weight loss.Int. J. Obes.321720–1724. 10.1038/ijo.2008.155

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  EdgarR. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads.Nat. Methods10996–998. 10.1038/nmeth.2604

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  EdgarR. C.HaasB. J.ClementeJ. C.QuinceC.KnightR. (2011). UCHIME improves sensitivity and speed of chimera detection.Bioinformatics272194–2200. 10.1093/bioinformatics/btr381

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  EvansN. J.BrownJ. M.MurrayR. D.GettyB.BirtlesR. J.HartC. A.et al (2011). Characterization of novel bovine gastrointestinal tract Treponema isolates and comparison with bovine digital dermatitis treponemes.Appl. Environ. Microbiol.77138–147. 10.1128/AEM.00993-10

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  FangW.FangZ.ZhouP.ChangF.HongY.ZhangX.et al (2012). Evidence for lignin oxidation by the giant panda fecal microbiome.PLoS One7:e50312. 10.1371/journal.pone.0050312

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  FernandoS. C.PurvisH.NajarF.SukharnikovL.KrehbielC.NagarajaT.et al (2010). Rumen microbial population dynamics during adaptation to a high-grain diet.Appl. Environ. Microbiol.767482–7490. 10.1128/AEM.00388-10

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  FoxJ. L.NurbuC.ChundawatR. S. (1991). The mountain ungulates of Ladakh, India.Biol. Conserv.58167–190. 10.1016/0006-3207(91)90118-s

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  GoodrichJ. K.WatersJ. L.PooleA. C.SutterJ. L.KorenO.BlekhmanR.et al (2014). Human genetics shape the gut microbiome.Cell159789–799. 10.1016/j.cell.2014.09.053

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  GorvitovskaiaA.HolmesS. P.HuseS. M. (2016). Interpreting Prevotella and Bacteroides as biomarkers of diet and lifestyle.Microbiome4:15. 10.1186/s40168-016-0160-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  GrovesC. (1993). “Order Diprotodontia,” in Mammal Species of the World: a Taxonomic and Geographic Reference, edsWilsonD. E.ReederD. A. (Washington, DC: Smithsonian Institution Press).

  • Google Scholar
  • View reference in article

• 21

  HamadN.TravisS. P. (2006). Weight loss at high altitude: pathophysiology and practical implications.Eur. J. Gastroenterol. Hepatol.185–10. 10.1097/00042737-200601000-00002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  HarrisR. B. (2014). Pseudois nayaur. The IUCN red list of threatened species. 2014: e.T61513537A64313015.

  • Google Scholar
  • View reference in article

• 23

  HarrisR. B.LoggersC. O. (2004). Status of Tibetan plateau mammals in Yeniugou, China.Wildlife Biol.1091–100.

  • Google Scholar
  • View reference in article

• 24

  HartmannG.TschöpM.FischerR.BidlingmaierC.RieplR.TschöpK.et al (2000). High altitude increases circulating interleukin-6, interleukin-1 receptor antagonist and C-reactive protein.Cytokine12246–252. 10.1006/cyto.1999.0533

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  HermansW. A. (1996). [The European mouflon, Ovis musimon].Tijdschr Diergeneeskd121515–517.

  • Google Scholar
  • View reference in article

• 26

  HochachkaP.BuckL.DollC.LandS. (1996). Unifying theory of hypoxia tolerance: molecular/metabolic defense and rescue mechanisms for surviving oxygen lack.Proc. Natl. Acad. Sci. U.S.A.939493–9498. 10.1073/pnas.93.18.9493

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  HoefsM. (1985). Wild Sheep: Distribution, Abundance, Management and Conservation of the Sheep of the World and Closely Related Mountain Ungulates.Whitehorse: Northern Wild Sheep and Goat Council.

  • Google Scholar
  • View reference in article

• 28

  IvyC. M.ScottG. R. (2015). Control of breathing and the circulation in high-altitude mammals and birds.Comp. Biochem. Physiol. A Mol. Integr. Physiol.18666–74. 10.1016/j.cbpa.2014.10.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  JindouS.BorovokI.RinconM. T.FlintH. J.AntonopoulosD. A.BergM. E.et al (2006). Conservation and divergence in cellulosome architecture between two strains of Ruminococcus flavefaciens.J. Bacteriol.1887971–7976. 10.1128/jb.00973-06

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  KaakoushN. O. (2015). Insights into the role of Erysipelotrichaceae in the human host.Front. Cell. Infect. Microbiol.5:84. 10.3389/fcimb.2015.00084

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  KarlJ. P.BerrymanC. E.YoungA. J.RadcliffeP. N.BranckT. A.Pantoja-FelicianoI. G.et al (2018). Associations between the gut microbiota and host responses to high altitude.Am. J. Physiol. Gastrointest. Liver Physiol.315G1003–G1015. 10.1152/ajpgi.00253.2018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  KhannaK.MishraK.GanjuL.KumarB.SinghS. B. (2018). High-altitude-induced alterations in gut-immune axis: a review.Int. Rev. Immunol.37119–126. 10.1080/08830185.2017.1407763

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  KongF.ZhaoJ.HanS.ZengB.YangJ.SiX.et al (2014). Characterization of the gut microbiota in the red panda (Ailurus fulgens).PLoS One9:e87885. 10.1371/journal.pone.0087885

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  LanD.JiW.LinB.ChenY.HuangC.XiongX.et al (2017). Correlations between gut microbiota community structures of Tibetans and geography.Sci. Rep.7:16982.

  • Google Scholar
  • View reference in article

• 35

  LeyR. E.HamadyM.LozuponeC.TurnbaughP. J.RameyR. R.BircherJ. S.et al (2008). Evolution of mammals and their gut microbes.Science3201647–1651. 10.1126/science.1155725

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  LeyR. E.TurnbaughP. J.KleinS.GordonJ. I. (2006). Microbial ecology: human gut microbes associated with obesity.Nature4441022–1023.

  • Pubmed Abstract
  • Google Scholar
  • View reference in article

• 37

  LiH.LiT.BeasleyD. E.HedìnecP.XiaoZ.ZhangS.et al (2016). Diet diversity is associated with beta but not alpha diversity of pika gut microbiota.Front. Microbiol.7:1169. 10.3389/fmicb.2016.01169

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  LiK.DanZ.GesangL.WangH.ZhouY.DuY.et al (2016). Comparative analysis of gut microbiota of native Tibetan and Han populations living at different altitudes.PLoS One11:e0155863. 10.1371/journal.pone.0155863

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  LiL.ZhaoX. (2015). Comparative analyses of fecal microbiota in Tibetan and Chinese Han living at low or high altitude by barcoded 454 pyrosequencing.Sci. Rep.5:14682. 10.1038/srep14682

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  LiuZ.WangX.LiZ.CuiD.LiX. (2005). Comparison of seasonal feeding habitats by blue sheep (Pseudois nayaur) during winter and spring in Helan mountain, China.Zool. Res.26580–589.

  • Google Scholar
  • View reference in article

• 41

  LiuZ.ZhuZ.GaoH.ZhaoC.SunY.TengL. (2018). Comparative analysis of winter diets and habitat use by the sympatric blue sheep (Pseudois nayaur) and Alashan red deer (Cervus alashanicus) in the Helan Mountains, China.Folia Zool.6743–54.

  • Google Scholar
  • View reference in article

• 42

  MagočT.SalzbergS. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies.Bioinformatics272957–2963. 10.1093/bioinformatics/btr507

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  MarschallL. A. (1998). Wildlife of the Tibetan Steppe.Sciences3853–55.

  • Google Scholar
  • View reference in article

• 44

  MaynardC. L.ElsonC. O.HattonR. D.WeaverC. T. (2012). Reciprocal interactions of the intestinal microbiota and immune system.Nature489231–241. 10.1038/nature11551

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  MillerA. W.OakesonK. F.DaleC.DearingM. D. (2016). Effect of dietary oxalate on the gut microbiota of the mammalian herbivore Neotoma albigula.Appl. Environ. Microbiol.822669–2675. 10.1128/AEM.00216-16

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  MurphyE.CotterP.HealyS.MarquesT. M.O’sullivanO.FouhyF.et al (2010). Composition and energy harvesting capacity of the gut microbiota: relationship to diet, obesity and time in mouse models.Gut591635–1642. 10.1136/gut.2010.215665

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  NamgailT.MishraC.De JongC. B.Van WierenS. E.PrinsH. H. (2009). Effects of herbivore species richness on the niche dynamics and distribution of blue sheep in the Trans-Himalaya.Divers. Distrib.15940–947. 10.1111/j.1472-4642.2009.00611.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  OkiK.ToyamaM.BannoT.ChonanO.BennoY.WatanabeK. (2016). Comprehensive analysis of the fecal microbiota of healthy Japanese adults reveals a new bacterial lineage associated with a phenotype characterized by a high frequency of bowel movements and a lean body type.BMC Microbiol.16:284. 10.1186/s12866-016-0898-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  OksanenJ.BlanchetF. G.KindtR.LegendreP.MinchinP. R.O’haraR.et al (2013). Package ‘Vegan’. Community Ecology Package, Version 2(9).

  • Google Scholar
  • View reference in article

• 50

  OuwerkerkJ. P.van der ArkK. C.DavidsM.ClaassensN. J.FinestraT. R.de VosW. M.et al (2016). Adaptation of Akkermansia muciniphila to the oxic-anoxic interface of the mucus layer.Appl. Environ. Microbiol.826983–6993. 10.1128/AEM.01641-16

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  PetersenC.RoundJ. L. (2014). Defining dysbiosis and its influence on host immunity and disease.Cell. Microbiol.161024–1033. 10.1111/cmi.12308

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  QinJ.LiY.CaiZ.LiS.ZhuJ.ZhangF.et al (2012). A metagenome-wide association study of gut microbiota in type 2 diabetes.Nature49055–60. 10.1038/nature11450

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  RabinowitzA.KhaingS. T. (1998). Status of selected mammal species in North Myanmar.Oryx32201–208. 10.1046/j.1365-3008.1998.d01-37.x

  • CrossRef
  • Google Scholar

• 54

  RamakrishnaB. S. (2013). Role of the gut microbiota in human nutrition and metabolism.J. Gastroenterol. Hepatol.289–17. 10.1111/jgh.12294

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  ReunanenJ.KainulainenV.HuuskonenL.OttmanN.BelzerC.HuhtinenH.et al (2015). Akkermansia muciniphila adheres to enterocytes and strengthens the integrity of the epithelial cell layer.Appl. Environ. Microbiol.813655–3662. 10.1128/AEM.04050-14

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  RezaeiH. R.NaderiS.Chintauan-MarquierI. C.TaberletP.VirkA. T.NaghashH. R.et al (2010). Evolution and taxonomy of the wild species of the genus Ovis (Mammalia, Artiodactyla, Bovidae).Mol. Phylogenet. Evol.54315–326. 10.1016/j.ympev.2009.10.037

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  RoundJ. L.MazmanianS. K. (2009). The gut microbiota shapes intestinal immune responses during health and disease.Nat. Rev. Immunol.9313–323. 10.1038/nri2515

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  SchallerG. (1998). Wildlife of the Tibetan Steppe.Chicago: University of Chicago.

  • Google Scholar

• 59

  SegataN.IzardJ.WaldronL.GeversD.MiropolskyL.GarrettW. S.et al (2011). Metagenomic biomarker discovery and explanation.Genome Biol.12:R60. 10.1186/gb-2011-12-6-r60

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  SekirovI.RussellS. L.AntunesL. C. M.FinlayB. B. (2010). Gut microbiota in health and disease.Physiol. Rev.90859–904. 10.1152/physrev.00045.2009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  ShackletonD. M. (1997). Wild Sheep and Goats and their Relatives.Gland: IUCN.

  • Google Scholar

• 62

  ShackletonD. M.International Union for Conservation of Nature, Natural Resources. Species Survival Commission, and Caprinae Specialist Group (1997). Wild Sheep and Goats and their Relatives : Status Survey and Conservation Action Plan for Caprinae.Gland: IUCN.

  • Google Scholar

• 63

  SimonsonT. S. (2015). Altitude adaptation: a glimpse through various lenses.High Altit. Med. Biol.16125–137. 10.1089/ham.2015.0033

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  StorzJ. F.ScottG. R.ChevironZ. A. (2010). Phenotypic plasticity and genetic adaptation to high-altitude hypoxia in vertebrates.J. Exp. Biol.2134125–4136. 10.1242/jeb.048181

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65

  SunB.WangX.BernsteinS.HuffmanM. A.XiaD.-P.GuZ.et al (2016). Marked variation between winter and spring gut microbiota in free-ranging Tibetan Macaques (Macaca thibetana).Sci. Rep.6:26035. 10.1038/srep26035

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

  SuzukiT. A.MartinsF. M.NachmanM. W. (2018). Altitudinal variation of the gut microbiota in wild house mice.Mol. Ecol.282378–2390. 10.1111/mec.14905

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 67

  TanS.ZouD.TangL.WangG.PengQ.ZengB.et al (2012). Molecular evidence for the subspecific differentiation of blue sheep (Pseudois nayaur) and polyphyletic origin of dwarf blue sheep (Pseudois schaeferi).Genetica140159–167. 10.1007/s10709-012-9667-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68

  TremaroliV.BäckhedF. (2012). Functional interactions between the gut microbiota and host metabolism.Nature489242–249. 10.1038/nature11552

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69

  TurnbaughP. J.LeyR. E.MahowaldM. A.MagriniV.MardisE. R.GordonJ. I. (2006). An obesity-associated gut microbiome with increased capacity for energy harvest.Nature4441027–1031.

  • Pubmed Abstract
  • Google Scholar

• 70

  ValdezR. (2008). Ovis orientalis. The IUCN red list of threatened species. 2008: e.T15739A5076068.

  • Google Scholar

• 71

  ValdezR.BattenJ. (1982). The Wild Sheep of the World.Mesilla, NM: Wild Sheep & Goat International.

  • Google Scholar

• 72

  VeeranagoudaY.HusainF.TenorioE. L.WexlerH. M. (2014). Identification of genes required for the survival of B. fragilis using massive parallel sequencing of a saturated transposon mutant library.BMC Genomics15:429. 10.1186/1471-2164-15-429

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  Vieira-SilvaS.FalonyG.DarziY.Lima-MendezG.YuntaR. G.OkudaS.et al (2016). Species–function relationships shape ecological properties of the human gut microbiome.Nat. Microbiol.1:16088. 10.1038/nmicrobiol.2016.88

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74

  WaiteD. W.TaylorM. W. (2014). Characterizing the avian gut microbiota: membership, driving influences, and potential function.Front. Microbiol.5:223. 10.3389/fmicb.2014.00223

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75

  WangQ.GarrityG. M.TiedjeJ. M.ColeJ. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy.Appl. Environ. Microbiol.735261–5267. 10.1128/aem.00062-07

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76

  WangS. (1998). China Red Data Book of Endangered Animals: Mammalia.Beijing: Science Press.

  • Google Scholar

• 77

  WesterterpK. R.KayserB. (2006). Body mass regulation at altitude.Eur. J. Gastroenterol. Hepatol.181–3. 10.1097/00042737-200601000-00001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 78

  WuG. D.ChenJ.HoffmannC.BittingerK.ChenY.-Y.KeilbaughS. A.et al (2011). Linking long-term dietary patterns with gut microbial enterotypes.Science334105–108. 10.1126/science.1208344

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 79

  YangJ.LiX.CaoY.-H.PokharelK.HuX.-J.ChenZ.-H.et al (2019). Comparative mRNA and miRNA expression in European mouflon (Ovis musimon) and sheep (Ovis aries) provides novel insights into the genetic mechanisms for female reproductive success.Heredity122172–186. 10.1038/s41437-018-0090-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 80

  YuqunR. J. Y. (1990). A study on the population structure and life table of blue sheep in yushu and golog, qinghai province.Acta Theriol. Sin.3189–193. 10.16829/j.slxb.1990.03.006

  • CrossRef
  • Google Scholar

• 81

  ZeX.DuncanS. H.LouisP.FlintH. J. (2012). Ruminococcus bromii is a keystone species for the degradation of resistant starch in the human colon.ISME J.61535–1543. 10.1038/ismej.2012.4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 82

  ZengB.XuL.YueB.LiZ.ZouF. (2008). Molecular phylogeography and genetic differentiation of blue sheep Pseudois nayaurszechuanensis and Pseudois schaeferi in China.Mol. Phylogenet. Evol.48387–395. 10.1016/j.ympev.2008.05.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 83

  ZhangZ.XuD.WangL.HaoJ.WangJ.ZhouX.et al (2016). Convergent evolution of rumen microbiomes in high-altitude mammals.Curr. Biol.261873–1879. 10.1016/j.cub.2016.05.012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 84

  ZhaoJ.YaoY.LiD.XuH.WuJ.WenA.et al (2018). Characterization of the Gut microbiota in six geographical populations of chinese Rhesus macaques (Macaca mulatta), implying an adaptation to high-altitude environment.Microb. Ecol.76565–577. 10.1007/s00248-018-1146-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar