Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 30 July 2019
Sec. Systems Microbiology

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

\r\nGuolei SunGuolei Sun1Honghai Zhang*Honghai Zhang1*Qinguo WeiQinguo Wei1Chao ZhaoChao Zhao1Xiufeng YangXiufeng Yang1Xiaoyang WuXiaoyang Wu1Tian XiaTian Xia1Guangshuai LiuGuangshuai Liu1Lei ZhangLei Zhang1Ying GaoYing Gao2Weilai ShaWeilai Sha1Ying LiYing Li2
  • 1College of Life Science, Qufu Normal University, Qufu, China
  • 2Wild World Jinan, Jinan, China

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.

TABLE 1
www.frontiersin.org

Table 1. Information for fecal samples from European mouflon and blue sheep.

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 ddH2O. 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.

TABLE 2
www.frontiersin.org

Table 2. Alpha-diversity of gut microbiota in fecal samples from European mouflon and blue sheep.

FIGURE 1
www.frontiersin.org

Figure 1. Rarefaction curves (A) and rank abundance curves (B) of each sample. The saturated rarefaction curves and species richness indices indicate that the sampling was comprehensive. In rank abundance curves, a wider span of curves reflects a higher relative species abundance, and the smoother the curve on the y-axis reveals a higher evenness of bacterial species in the fecal samples of European mouflon and blue sheep.

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
www.frontiersin.org

Figure 2. Fecal microbial composition of European mouflon and blue sheep at the phylum (A) and genus (B) levels. K, low-altitude mouflon; L, high-altitude mouflon; M, low-altitude blue sheep; N, high-altitude blue sheep.

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
www.frontiersin.org

Figure 3. A heatmap of species abundance clustering. The heatmap shows the hierarchical clustering of samples based on the relative abundance of the top ranked 35 genera of fecal microbiota in of European mouflon and blue sheep. The relative values in the heatmap (after normalization), depicted by colors, indicate the aggregation degree or content of bacterial species among samples at the genus level. The relative values were normalized in four individual groups and displayed in figure on the right-hand side. K, low-altitude mouflon; L, high-altitude mouflon; M, low-altitude blue sheep; N, high-altitude blue sheep.

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
www.frontiersin.org

Figure 4. Cluster analysis with the Unifracdistance matrix. The top figure is based on the unweighted Unifrac distance and the clustering in groups. The bottom figure is based on the weighted Unifrac distance. K, low-altitude mouflon; L, high-altitude mouflon; M, low-altitude blue sheep; N, high-altitude blue sheep.

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.

TABLE 3
www.frontiersin.org

Table 3. ANOSIM analysis between the four groups of European mouflon and blue sheep.

TABLE 4
www.frontiersin.org

Table 4. Multi-Response Permutation Procedures difference analysis between the four groups of European mouflon and blue sheep.

FIGURE 5
www.frontiersin.org

Figure 5. Comparisons for alpha-diversity (Shannon and observed species index) and beta-diversity (with weighted and unweighted Unifrac distance matrix) of fecal the microbiota of between the low- and high-altitude groups of sheep. K, low-altitude mouflon; L, high-altitude mouflon; M, low-altitude blue sheep; N, high-altitude blue sheep.

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
www.frontiersin.org

Figure 6. Heatmap of the beta-diversity of microbiota. The numbers in the grid blocks represent the dissimilarity coefficients between fecal samples. The top number in each block represents the weighted Unifrac distance and the bottom number represents the unweighted Unifrac distance. K, low-altitude mouflon; L, high-altitude mouflon; M, low-altitude blue sheep; N, high-altitude blue sheep.

FIGURE 7
www.frontiersin.org

Figure 7. Non-metric multidimensional scaling (NMDS) and principal component analysis (PCA) of bacterial community structures. Each symbol represents the gut microbiota of an individual and the different colors represent the different groups. K, low-altitude mouflon; L, high-altitude mouflon; M, low-altitude blue sheep; N, high-altitude blue sheep.

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).

TABLE 5
www.frontiersin.org

Table 5. Relative abundance of the top 10 most abundant bacterial phyla in the European mouflon sheep.

TABLE 6
www.frontiersin.org

Table 6. Relative abundance of the top 10 most abundant bacterial phyla in blue sheep.

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
www.frontiersin.org

Figure 8. LEfSe (LDA Effect Size) analysis (A) and the cladogram (B) in European mouflon. The different colors represent the different groups. The histogram (A) shows the biomarkers with statistical differences between groups and the lengths of the bars indicate the influential degree of the species. The cladogram (B) demonstrates the classification of taxa at five levels, and the red and green circles indicate the differences in relative abundance. Non-significant differences are expressed by the yellow circles. K, low-altitude mouflon; L, high-altitude mouflon.

FIGURE 9
www.frontiersin.org

Figure 9. LEfSe (LDA Effect Size) analysis (A) and the cladogram (B) in blue sheep. The different colors represent different groups. The histogram (A) shows the biomarkers with statistical difference between groups and the lengths of the bars indicate the influential degree of species. The cladogram (B) demonstrates the classification of taxa at five levels, and the red and green circles indicate the differences in relative abundance. Non-significant differences are expressed by the yellow circles. M, low-altitude blue sheep; N, high-altitude blue sheep.

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.

Author Contributions

HZ, QW, GL, CZ, and WS designed the study. GS, XY, TX, and LZ conducted the research. GS, XW, YG, and YL analyzed the data. GS and QW prepared the manuscript.

Funding

This work was supported by the Special Fund for Forest Scientific Research in the Public Welfare (201404420) and the National Natural Science Foundation of China (31872242, 31672313, and 31372220).

Conflict of Interest Statement

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.

Acknowledgments

GS would like to show his deepest gratitude to his supervisor, HZ, a respectable, responsible, and resourceful scholar, who has provided him with the valuable guidance in every stage of the writing of this thesis.

References

Avila-Jaime, B., Kawas, J., and Garcia-Mazcorro, J. (2018). Prediction of functional metagenomic composition using archived 16S rDNA sequence data from the gut microbiota of livestock. Livest. Sci. 213, 28–34. doi: 10.1016/j.livsci.2018.04.017

CrossRef Full Text | Google Scholar

Barbato, M., Hailer, F., Orozco-terWengel, P., Kijas, J., Mereu, P., Cabras, P., et al. (2017). Genomic signatures of adaptive introgression from European mouflon into domestic sheep. Sci. Rep. 7:7623. doi: 10.1038/s41598-017-07382-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Bergmann, G. T., Craine, J. M., Robeson, M. S. II, and Fierer, N. (2015). Seasonal shifts in diet and gut microbiota of the American bison (Bison bison). PLoS One 10:e0142409. doi: 10.1371/journal.pone.0142409

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

Castillo, M., Martín-Orúe, S., Anguita, M., Pérez, J., and Gasa, J. (2007). Adaptation of gut microbiota to corn physical structure and different types of dietary fibre. Livest. Sci. 109, 149–152. doi: 10.1016/j.livsci.2007.01.129

CrossRef Full Text | Google Scholar

Cheviron, Z. A., Bachman, G. C., Connaty, A. D., McClelland, G. B., and Storz, J. F. (2012). Regulatory changes contribute to the adaptive enhancement of thermogenic capacity in high-altitude deer mice. Proc. Natl. Acad. Sci. U.S.A. 109, 8635–8640. doi: 10.1073/pnas.1120523109

PubMed Abstract | CrossRef Full Text | Google Scholar

Costea, P. I., Hildebrand, F., Arumugam, M., Bäckhed, F., Blaser, M. J., Bushman, F. D., et al. (2018). Enterotypes in the landscape of gut microbial community composition. Nat. Microbiol. 3, 8–16. doi: 10.1038/s41564-017-0072-8

PubMed Abstract | CrossRef Full Text | Google Scholar

De Filippo, C., Cavalieri, D., Di Paola, M., Ramazzotti, M., Poullet, J. B., Massart, S., 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. 107, 14691–14696. doi: 10.1073/pnas.1005963107

PubMed Abstract | CrossRef Full Text | Google Scholar

den Besten, G., van Eunen, K., Groen, A. K., Venema, K., Reijngoud, D.-J., and Bakker, B. M. (2013). The role of short-chain fatty acids in the interplay between diet, gut microbiota, and host energy metabolism. J. Lipid Res. 54, 2325–2340. doi: 10.1194/jlr.R036012

PubMed Abstract | CrossRef Full Text | Google Scholar

DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/aem.03006-05

PubMed Abstract | CrossRef Full Text | Google Scholar

Duncan, S. H., Lobley, G., Holtrop, G., Ince, J., Johnstone, A., Louis, P., et al. (2008). Human colonic microbiota associated with diet, obesity and weight loss. Int. J. Obes. 32, 1720–1724. doi: 10.1038/ijo.2008.155

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381

PubMed Abstract | CrossRef Full Text | Google Scholar

Evans, N. J., Brown, J. M., Murray, R. D., Getty, B., Birtles, R. J., Hart, C. A., et al. (2011). Characterization of novel bovine gastrointestinal tract Treponema isolates and comparison with bovine digital dermatitis treponemes. Appl. Environ. Microbiol. 77, 138–147. doi: 10.1128/AEM.00993-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, W., Fang, Z., Zhou, P., Chang, F., Hong, Y., Zhang, X., et al. (2012). Evidence for lignin oxidation by the giant panda fecal microbiome. PLoS One 7:e50312. doi: 10.1371/journal.pone.0050312

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernando, S. C., Purvis, H., Najar, F., Sukharnikov, L., Krehbiel, C., Nagaraja, T., et al. (2010). Rumen microbial population dynamics during adaptation to a high-grain diet. Appl. Environ. Microbiol. 76, 7482–7490. doi: 10.1128/AEM.00388-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Fox, J. L., Nurbu, C., and Chundawat, R. S. (1991). The mountain ungulates of Ladakh, India. Biol. Conserv. 58, 167–190. doi: 10.1016/0006-3207(91)90118-s

CrossRef Full Text | Google Scholar

Goodrich, J. K., Waters, J. L., Poole, A. C., Sutter, J. L., Koren, O., Blekhman, R., et al. (2014). Human genetics shape the gut microbiome. Cell 159, 789–799. doi: 10.1016/j.cell.2014.09.053

PubMed Abstract | CrossRef Full Text | Google Scholar

Gorvitovskaia, A., Holmes, S. P., and Huse, S. M. (2016). Interpreting Prevotella and Bacteroides as biomarkers of diet and lifestyle. Microbiome 4:15. doi: 10.1186/s40168-016-0160-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Groves, C. (1993). “Order Diprotodontia,” in Mammal Species of the World: a Taxonomic and Geographic Reference, eds D. E. Wilson and D. A. Reeder (Washington, DC: Smithsonian Institution Press).

Google Scholar

Hamad, N., and Travis, S. P. (2006). Weight loss at high altitude: pathophysiology and practical implications. Eur. J. Gastroenterol. Hepatol. 18, 5–10. doi: 10.1097/00042737-200601000-00002

PubMed Abstract | CrossRef Full Text | Google Scholar

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

Google Scholar

Harris, R. B., and Loggers, C. O. (2004). Status of Tibetan plateau mammals in Yeniugou, China. Wildlife Biol. 10, 91–100.

Google Scholar

Hartmann, G., Tschöp, M., Fischer, R., Bidlingmaier, C., Riepl, R., Tschöp, K., et al. (2000). High altitude increases circulating interleukin-6, interleukin-1 receptor antagonist and C-reactive protein. Cytokine 12, 246–252. doi: 10.1006/cyto.1999.0533

PubMed Abstract | CrossRef Full Text | Google Scholar

Hermans, W. A. (1996). [The European mouflon, Ovis musimon]. Tijdschr Diergeneeskd 121, 515–517.

Google Scholar

Hochachka, P., Buck, L., Doll, C., and Land, S. (1996). Unifying theory of hypoxia tolerance: molecular/metabolic defense and rescue mechanisms for surviving oxygen lack. Proc. Natl. Acad. Sci. U.S.A. 93, 9493–9498. doi: 10.1073/pnas.93.18.9493

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoefs, M. (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

Ivy, C. M., and Scott, G. R. (2015). Control of breathing and the circulation in high-altitude mammals and birds. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 186, 66–74. doi: 10.1016/j.cbpa.2014.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Jindou, S., Borovok, I., Rincon, M. T., Flint, H. J., Antonopoulos, D. A., Berg, M. E., et al. (2006). Conservation and divergence in cellulosome architecture between two strains of Ruminococcus flavefaciens. J. Bacteriol. 188, 7971–7976. doi: 10.1128/jb.00973-06

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Karl, J. P., Berryman, C. E., Young, A. J., Radcliffe, P. N., Branck, T. A., Pantoja-Feliciano, I. G., et al. (2018). Associations between the gut microbiota and host responses to high altitude. Am. J. Physiol. Gastrointest. Liver Physiol. 315, G1003–G1015. doi: 10.1152/ajpgi.00253.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Khanna, K., Mishra, K., Ganju, L., Kumar, B., and Singh, S. B. (2018). High-altitude-induced alterations in gut-immune axis: a review. Int. Rev. Immunol. 37, 119–126. doi: 10.1080/08830185.2017.1407763

PubMed Abstract | CrossRef Full Text | Google Scholar

Kong, F., Zhao, J., Han, S., Zeng, B., Yang, J., Si, X., et al. (2014). Characterization of the gut microbiota in the red panda (Ailurus fulgens). PLoS One 9:e87885. doi: 10.1371/journal.pone.0087885

PubMed Abstract | CrossRef Full Text | Google Scholar

Lan, D., Ji, W., Lin, B., Chen, Y., Huang, C., Xiong, X., et al. (2017). Correlations between gut microbiota community structures of Tibetans and geography. Sci. Rep. 7:16982.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Ley, R. E., Turnbaugh, P. J., Klein, S., and Gordon, J. I. (2006). Microbial ecology: human gut microbes associated with obesity. Nature 444, 1022–1023.

PubMed Abstract | Google Scholar

Li, H., Li, T., Beasley, D. E., Hedìnec, P., Xiao, Z., Zhang, S., et al. (2016). Diet diversity is associated with beta but not alpha diversity of pika gut microbiota. Front. Microbiol. 7:1169. doi: 10.3389/fmicb.2016.01169

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, K., Dan, Z., Gesang, L., Wang, H., Zhou, Y., Du, Y., et al. (2016). Comparative analysis of gut microbiota of native Tibetan and Han populations living at different altitudes. PLoS One 11:e0155863. doi: 10.1371/journal.pone.0155863

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., and Zhao, X. (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. doi: 10.1038/srep14682

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z., Wang, X., Li, Z., Cui, D., and Li, X. (2005). Comparison of seasonal feeding habitats by blue sheep (Pseudois nayaur) during winter and spring in Helan mountain, China. Zool. Res. 26, 580–589.

Google Scholar

Liu, Z., Zhu, Z., Gao, H., Zhao, C., Sun, Y., and Teng, L. (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. 67, 43–54.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Marschall, L. A. (1998). Wildlife of the Tibetan Steppe. Sciences 38, 53–55.

Google Scholar

Maynard, C. L., Elson, C. O., Hatton, R. D., and Weaver, C. T. (2012). Reciprocal interactions of the intestinal microbiota and immune system. Nature 489, 231–241. doi: 10.1038/nature11551

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, A. W., Oakeson, K. F., Dale, C., and Dearing, M. D. (2016). Effect of dietary oxalate on the gut microbiota of the mammalian herbivore Neotoma albigula. Appl. Environ. Microbiol. 82, 2669–2675. doi: 10.1128/AEM.00216-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Murphy, E., Cotter, P., Healy, S., Marques, T. M., O’sullivan, O., Fouhy, F., et al. (2010). Composition and energy harvesting capacity of the gut microbiota: relationship to diet, obesity and time in mouse models. Gut 59, 1635–1642. doi: 10.1136/gut.2010.215665

PubMed Abstract | CrossRef Full Text | Google Scholar

Namgail, T., Mishra, C., De Jong, C. B., Van Wieren, S. E., and Prins, H. H. (2009). Effects of herbivore species richness on the niche dynamics and distribution of blue sheep in the Trans-Himalaya. Divers. Distrib. 15, 940–947. doi: 10.1111/j.1472-4642.2009.00611.x

CrossRef Full Text | Google Scholar

Oki, K., Toyama, M., Banno, T., Chonan, O., Benno, Y., and Watanabe, K. (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. doi: 10.1186/s12866-016-0898-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O’hara, R., et al. (2013). Package ‘Vegan’. Community Ecology Package, Version 2(9).

Google Scholar

Ouwerkerk, J. P., van der Ark, K. C., Davids, M., Claassens, N. J., Finestra, T. R., de Vos, W. M., et al. (2016). Adaptation of Akkermansia muciniphila to the oxic-anoxic interface of the mucus layer. Appl. Environ. Microbiol. 82, 6983–6993. doi: 10.1128/AEM.01641-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Petersen, C., and Round, J. L. (2014). Defining dysbiosis and its influence on host immunity and disease. Cell. Microbiol. 16, 1024–1033. doi: 10.1111/cmi.12308

PubMed Abstract | CrossRef Full Text | Google Scholar

Qin, J., Li, Y., Cai, Z., Li, S., Zhu, J., Zhang, F., et al. (2012). A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature 490, 55–60. doi: 10.1038/nature11450

PubMed Abstract | CrossRef Full Text | Google Scholar

Rabinowitz, A., and Khaing, S. T. (1998). Status of selected mammal species in North Myanmar. Oryx 32, 201–208. doi: 10.1046/j.1365-3008.1998.d01-37.x

CrossRef Full Text | Google Scholar

Ramakrishna, B. S. (2013). Role of the gut microbiota in human nutrition and metabolism. J. Gastroenterol. Hepatol. 28, 9–17. doi: 10.1111/jgh.12294

PubMed Abstract | CrossRef Full Text | Google Scholar

Reunanen, J., Kainulainen, V., Huuskonen, L., Ottman, N., Belzer, C., Huhtinen, H., et al. (2015). Akkermansia muciniphila adheres to enterocytes and strengthens the integrity of the epithelial cell layer. Appl. Environ. Microbiol. 81, 3655–3662. doi: 10.1128/AEM.04050-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Rezaei, H. R., Naderi, S., Chintauan-Marquier, I. C., Taberlet, P., Virk, A. T., Naghash, H. R., et al. (2010). Evolution and taxonomy of the wild species of the genus Ovis (Mammalia, Artiodactyla, Bovidae). Mol. Phylogenet. Evol. 54, 315–326. doi: 10.1016/j.ympev.2009.10.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Round, J. L., and Mazmanian, S. K. (2009). The gut microbiota shapes intestinal immune responses during health and disease. Nat. Rev. Immunol. 9, 313–323. doi: 10.1038/nri2515

PubMed Abstract | CrossRef Full Text | Google Scholar

Schaller, G. (1998). Wildlife of the Tibetan Steppe. Chicago: University of Chicago.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Sekirov, I., Russell, S. L., Antunes, L. C. M., and Finlay, B. B. (2010). Gut microbiota in health and disease. Physiol. Rev. 90, 859–904. doi: 10.1152/physrev.00045.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Shackleton, D. M. (1997). Wild Sheep and Goats and their Relatives. Gland: IUCN.

Google Scholar

Shackleton, D. 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

Simonson, T. S. (2015). Altitude adaptation: a glimpse through various lenses. High Altit. Med. Biol. 16, 125–137. doi: 10.1089/ham.2015.0033

PubMed Abstract | CrossRef Full Text | Google Scholar

Storz, J. F., Scott, G. R., and Cheviron, Z. A. (2010). Phenotypic plasticity and genetic adaptation to high-altitude hypoxia in vertebrates. J. Exp. Biol. 213, 4125–4136. doi: 10.1242/jeb.048181

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, B., Wang, X., Bernstein, S., Huffman, M. A., Xia, D.-P., Gu, Z., et al. (2016). Marked variation between winter and spring gut microbiota in free-ranging Tibetan Macaques (Macaca thibetana). Sci. Rep. 6:26035. doi: 10.1038/srep26035

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, T. A., Martins, F. M., and Nachman, M. W. (2018). Altitudinal variation of the gut microbiota in wild house mice. Mol. Ecol. 28, 2378–2390. doi: 10.1111/mec.14905

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, S., Zou, D., Tang, L., Wang, G., Peng, Q., Zeng, B., et al. (2012). Molecular evidence for the subspecific differentiation of blue sheep (Pseudois nayaur) and polyphyletic origin of dwarf blue sheep (Pseudois schaeferi). Genetica 140, 159–167. doi: 10.1007/s10709-012-9667-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Tremaroli, V., and Bäckhed, F. (2012). Functional interactions between the gut microbiota and host metabolism. Nature 489, 242–249. doi: 10.1038/nature11552

PubMed Abstract | CrossRef Full Text | Google Scholar

Turnbaugh, P. J., Ley, R. E., Mahowald, M. A., Magrini, V., Mardis, E. R., and Gordon, J. I. (2006). An obesity-associated gut microbiome with increased capacity for energy harvest. Nature 444, 1027–1031.

PubMed Abstract | Google Scholar

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

Google Scholar

Valdez, R., and Batten, J. (1982). The Wild Sheep of the World. Mesilla, NM: Wild Sheep & Goat International.

Google Scholar

Veeranagouda, Y., Husain, F., Tenorio, E. L., and Wexler, H. M. (2014). Identification of genes required for the survival of B. fragilis using massive parallel sequencing of a saturated transposon mutant library. BMC Genomics 15:429. doi: 10.1186/1471-2164-15-429

PubMed Abstract | CrossRef Full Text | Google Scholar

Vieira-Silva, S., Falony, G., Darzi, Y., Lima-Mendez, G., Yunta, R. G., Okuda, S., et al. (2016). Species–function relationships shape ecological properties of the human gut microbiome. Nat. Microbiol. 1:16088. doi: 10.1038/nmicrobiol.2016.88

PubMed Abstract | CrossRef Full Text | Google Scholar

Waite, D. W., and Taylor, M. W. (2014). Characterizing the avian gut microbiota: membership, driving influences, and potential function. Front. Microbiol. 5:223. doi: 10.3389/fmicb.2014.00223

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive 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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

Google Scholar

Westerterp, K. R., and Kayser, B. (2006). Body mass regulation at altitude. Eur. J. Gastroenterol. Hepatol. 18, 1–3. doi: 10.1097/00042737-200601000-00001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, G. D., Chen, J., Hoffmann, C., Bittinger, K., Chen, Y.-Y., Keilbaugh, S. A., et al. (2011). Linking long-term dietary patterns with gut microbial enterotypes. Science 334, 105–108. doi: 10.1126/science.1208344

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Li, X., Cao, Y.-H., Pokharel, K., Hu, X.-J., Chen, Z.-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. Heredity 122, 172–186. doi: 10.1038/s41437-018-0090-1

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Ze, X., Duncan, S. H., Louis, P., and Flint, H. J. (2012). Ruminococcus bromii is a keystone species for the degradation of resistant starch in the human colon. ISME J. 6, 1535–1543. doi: 10.1038/ismej.2012.4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, B., Xu, L., Yue, B., Li, Z., and Zou, F. (2008). Molecular phylogeography and genetic differentiation of blue sheep Pseudois nayaurszechuanensis and Pseudois schaeferi in China. Mol. Phylogenet. Evol. 48, 387–395. doi: 10.1016/j.ympev.2008.05.010

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, J., Yao, Y., Li, D., Xu, H., Wu, J., Wen, A., 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. 76, 565–577. doi: 10.1007/s00248-018-1146-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: European mouflon, blue sheep, gut microbiota, 16S rRNA gene, altitude

Citation: Sun G, Zhang H, Wei Q, Zhao C, Yang X, Wu X, Xia T, Liu G, Zhang L, Gao Y, Sha W and Li Y (2019) Comparative Analyses of Fecal Microbiota in European Mouflon (Ovis orientalis musimon) and Blue Sheep (Pseudois nayaur) Living at Low or High Altitudes. Front. Microbiol. 10:1735. doi: 10.3389/fmicb.2019.01735

Received: 25 December 2018; Accepted: 15 July 2019;
Published: 30 July 2019.

Edited by:

George Tsiamis, University of Patras, Greece

Reviewed by:

Aggeliki Saridaki, Technical University of Crete, Greece
Edward Pajarillo, Florida Agricultural and Mechanical University, United States

Copyright © 2019 Sun, Zhang, Wei, Zhao, Yang, Wu, Xia, Liu, Zhang, Gao, Sha and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Honghai Zhang, zhanghonghai67@126.com

Disclaimer: 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.