Skip to main content

Comparison of gut microbiota diversity between wild and captive bharals (Pseudois nayaur)



Gastrointestinal microbiota play an important role in animal host immunity, nutrient metabolism, and energy acquisition, and have therefore drawn increasing attentions. This study compared the diversity of the gut microbiota of both wild and captive bharals, which is an ungulate herbivore of caprid from the Qinghai-Tibet plateau.


The sequencing of the V4-V5 region of the 16S rRNA gene via high-throughput sequencing technology showed that the dominant bacterial phyla are Firmicutes and Bacteroides both in wild and captive bharals. However, their abundance differed significantly between groups. Firmicutes were significantly higher in wild bharals, while Bacteroides were significantly higher in captive bharals. Different diets are likely a key influencing factor in the diversity and abundance of gut microbiota in bharals.


Changes in diets affect the diversity of gut microbiota and the relative abundance of pathogenic bacteria, increasing the risk of diseases outbreak in captive bharals. The results of this study suggest that the structure and function of the gut microbiota should be regulated via dietary intervention, accurate provision of an individualized diet, and optimization of the functional network of gut microbiota and its interaction with the host. This will improve the ex situ protection of wild animals.


Due to the use of modern molecular technology, the gut microbiota can be utilized as a signal hub that combines environmental inputs (e.g., diet) with genetic and immune signals. All of this affect the host metabolism, immunity, and infection responses [1], and play an important role in the development of the immune system and in animal health [2,3,4,5,6]. The diversity and abundance of the host intestinal flora are influenced by factors such as species [7], food [8,9,10], genotype, and age of the host [11]. Further research indicated that changes in the dietary patterns of the host can lead to rapid changes in the structure of the microbial community [12], which in turn exerts a profound impact on the health of the host [13]. The role of the diet in the regulation of the composition and metabolic activities of the gut microbiota has been increasingly recognized. Furthermore, the relationship between gut microbiota and animal health has been extensively studied.

The bharal (Pseudois nayaur) is a member of the Artiodactyla, Bovidae, Caprinae, Pseudois, and a national second-class protected animal on the list of key wildlife protection in China [14]. The bharal is mainly distributed throughout the Qinghai-Tibet Plateau and its surrounding areas, including Tibet, Yunnan, Sichuan, Xinjiang, Qinghai, Gansu, Inner Mongolia, Ningxia and Shaanxi [14]. It is one of the large hoofed animals in China and presents the most widely distributed cloven hoofed animal with the largest number on the Qinghai-Tibet Plateau [15, 16]. The major predator of wild bharals is the snow leopard (Uncia uncia) [17, 18], and this predation is important to maintain both the stability of ecosystems and species diversity [15]. At present, research on bharals focuses on population ecology [14, 19], behavioral ecology [20], and system evolution [21, 22]. During winter, bharals often face food shortage, which even lead to death [14]. Therefore, it is very important to study the digestibility and utilization rate of food in bharals. Its gut microbiota exerts a very important influence on this function. Artificial captive protection is an effective means of protection; however, it remains unclear whether the health status of captive bharals will be affected when animals are faced with novel living conditions. Therefore, it is necessary to study the relationship between changes in the gut microbiota and animal health.

This study investigated the composition and structure of the gut microbiota of bharals. The results lead to a better understanding of the digestive mechanism, provide a theoretical basis for the monitoring of abnormal physiological status, control the occurrence of diseases, and optimize the energy conversion rate of the food. The findings provide an important theoretical basis for the research digestive physiology of captive bharals. The results provide information to improve the diet of animals as well as for diagnosis and treatment of intestinal diseases. Furthermore, new research directions are initiated for the development of intestinal microecological agents.


Sequencing data

A total of 8,217,442 high-quality reads were obtained after data quality control, and fecal samples were classified into 3,878 OTUs (operational taxonomic units), 2,443 of which were in the captive group and 3,166 in the wild group.

Rarefaction curves and rank abundance curves are commonly used to describe the diversity of samples within a group. The rarefaction curves directly reflect the rationality of the sequencing data volume, and indirectly reflect the species richness of samples, as shown in Fig. 1a. Since the curves are smooth, a higher data volume would only yield a low number of OTUs, indicating that the volume of sequencing data is sufficiently reasonable. Rank abundance curves provide a visual representation of species richness and sample uniformity. In the horizontal direction, a greater span of the curve indicates a higher species richness, while in the vertical direction, a smoother curve indicates a more homogeneous species distribution [23] (Fig. 1b).

Fig. 1

Rarefaction curves (a) and rank abundance curves (b). CW for captive group, WW for wild group

Bacteria composition and relative abundance

In order to show the relative abundance of bacteria communities more intuitively, we selected the top 10 taxa and generated the relative abundance superposition histogram at phylum and genus level respectively in Fig. 2.

Fig. 2

Relative abundance histogram. Fecal microbial composition of wild (WW) and captive (CW) bharals at the phylum (a) and genus (b) level. Each bar represents the top ten bacterial species ranked by the relative abundance in each group

At the phylum level, Firmicutes and Bacteroidetes are the dominant phyla of both wild and captive groups; however, the abundance of Firmicutes (60.35%) in wild bharals was significantly higher (P < 0.01) than that of captive bharals (49.68%). The abundance of Bacteroides (40.98%) was significantly higher (P < 0.01) in captive bharals than in wild bharals (31.35%). In wild and captive groups, the main genera are Bacteroides and Alistipes, and the abundances of Bacteroides and Alistipes in captive bharals were significantly higher (P < 0.05) than that in wild bharals.

Analysis of discrepancies between groups

The goods coverage index exceeded 99%, indicating a high level of diversity coverage in the samples. The observed species, Shannon index, and Simpson index in the wild group were significantly higher than in the captive group (PObserved species < 0.01, PShannon < 0.01, and PSimpson < 0.01).

To further analyze whether significant differences affected the microbiome structure between the wild and captive group, the multi response permutation procedure (MRPP) significance analysis method was used. The result was A = 0.1153 > 0, indicating that differences between groups were higher than within groups; consequently, the applied study grouping was reasonable and a significant difference (P = 0.001) was found between wild and captive groups.

Furthermore, to assess the differences between wild and captive groups, principal component analysis (PCA) (Fig. 3a) and principal coordinates analysis (PCoA) (Fig. 3b) were used. The distance between the points in all graphs reflects the degree of similarity of their respective microbial flora structures. All graphs reflect the obvious difference of samples in each group.

Fig. 3

PCA (a) and PCoA (b) plot of the bacterial population structures. The green and red dots represented wild (WW) and captive (CW) bharals samples respectively

Using the LDA (local-density approximation) Effect Size (LEfSe) analysis method, species with significant differences between groups were selected. The results included three parts: a LDA value distribution histogram, an evolutionary branch diagram (phylogenetic distribution), and an abundance comparison diagram of biomarkers with statistical differences (LDA SCORE > 4) between groups (Fig. 4). According to the LDA scores, biomarkers with statistically significant differences between both groups are listed. The microbial communities that play an important role in each group were identified with the evolutionary branch diagram. In the wild group, Firmicutes, Clostridia, Clostridiales, Ruminococcaceae, Christensenellaceae and Lachnospiraceae are very important. In the captive group, important species are Bacteroidetes, Bacteroidia, Bacteroidales, Rikenellaceae and Alistipes.

Fig. 4

The results of LEfSe (LDA Effect Size) analysis. The histogram of LDA score(a) showed the biomarkers with significant differences between groups. The length of the column (LDA Score) represents the influencing degree of biomarkers, In the cladogram(b), the circle radiated inside-out represented the classification of phylum to genus level. Each small circle at different classification levels represented a taxon and the diameter of small circle is proportional to the relative abundance. The species not with significant differences were colored by yellow and biomarkers were colored by different groups. CW for captive group, WW for wild group


At the phylum level, the core microflora of the gut microbiota consisted of Firmicutes and Bacteroides in both captive and wild bharals, accounting for more than 90% of the total gut microbiota. This result coincided with previous studies of the gut microbiota in other ruminants [24,25,26,27]. However, their abundances showed significant differences (P < 0.01) between groups. In ruminants, Firmicutes play an important role in fiber and cellulose degradation, and can degrade cellulose into volatile fatty acids, which can be used by the host [28]. Firmicutes and Bacteroidetes increase with increasing hay content in the food, and this characteristic of Firmicutes is particularly prominent [29]. Although wild bharals have access to a wide variety of food sources, due to the harshness of the environment in winter, they mainly feed on grasses with high fibre content. Compared to captive bharals, wild bharals showed a higher abundance of Firmicutes, which leads to improved digestion and absorption of nutrients. Bacteroides can promote digestion, decompose polysaccharides and proteins, improve the utilization rate of assimilated nutrients [30, 31], and maintain the balance of intestinal microecological system [32]. The diet of captive bharals consists of semi-dried oat grass (Arrhenatherum elatius), carrot (Daucus carota), and artificial fodder. The utilized recipes are simple in structure; however, the protein, polysaccharide, and fat contents of the artificial fodder are relatively high. This might be the reason for the higher abundance of Bacteroides in the intestinal flora of captive bharals compared to wild bharals.

LEfSe analysis showed a number of bacteria that play an important role in metabolism of nutrients and the health of the host. Previous studies have reported that Alistipes shows a correlation with metabolites such as short-chain fatty acids, oligosaccharides, and amino acids. Lachnospiraceae and Ruminococcaceae belong to butyrate-producing bacteria. These bacteria and metabolites exert an important impact on host health [33, 34].

In addition, the Alpha-diversity of the captive group was significantly lower than that of the wild group. The food sources accessible to wild animals are more diverse, thus providing bharals with many different types of nutrition, which may require a more diversified gut microbiota to help bharals complete digestion and utilization of these nutrients [35, 36].

Further research has shown that many chronic diseases are caused by a decrease in the diversity of gut microbiota, including inflammatory gut disease and diabetes [37]. Comparison of the inter-group differences of the intestinal flora showed that the relative abundances of Spirochaetes (P < 0.01), Acidobacteria (P < 0.01) and Gemmatimonadetes (P < 0.01) in the captive group were significantly higher than in the wild group; these bacteria showed a certain pathogenicity [38, 39]. However, Ruminococcaceae, Lachnospiraceae, and Christensenellaceae in the captive group were significantly lower than in the wild group; these bacteria are beneficial microbiota in the host intestines [40]. These results suggest that captive bharals may face a higher risk to catch diseases than wild bharals.

This relationship between gut microbiota and host health, and in particular, the impact of microbial dissonance on host health, indicates the importance of feasible strategies to optimize the gut microbiota via diet. This means that diseases caused by changes in the gut microbiome can be regulated through the diet of the host. Food types can be reasonably chosen, a number of nutrients can be appropriately supplemented, and the healthy development of gut microbiota can be promoted. This will be conducive to the protection of wildlife ex situ.


This study described and compared the gut microbiota of wild and captive bharals and significant differences were found between both groups. The difference in diet may have caused a decrease of gut microbiotic diversity and an increase of relative abundance of pathogenic bacteria. The result would be an increased susceptibility to diseases in captive bharals, which is extremely unfavorable for their protection. These findings provide further research directions for the study of the effect of gut microbiota on the growth and development of bharals. Moreover, this study provides theoretical guidance for the diagnosis and detection of intestinal diseases in bharals, and can also be used as a reasonable and balanced reference for the diet of captive bharals. The findings provide a theoretical basis for the screening of probiotics and development of intestinal microecological agents.


Feces samples collection

The difficulties of wildlife sampling and the religious beliefs of ethnic minorities are considered, feces samples were selected as test material in this study. At present, most of the research on wildlife in this aspect takes fecal samples as research objects, the results describe the structure of the gut microbiota of bharals. Therefore, during January 2018, a total of 240 fresh feces samples of wild bharals were collected from different regions near the Donggeicuona Lake in Maduo County, Qing Hai, China. The average temperature of Maduo in January is − 16.8 °C, which retained the freshness of the feces of wild bharals as much as possible. The collected samples were temporarily stored in an in-vehicle refrigerator (− 20 °C).

Forty-four fresh feces samples were collected from 11 bharals in captivity at the Qinghai-Tibet Plateau Wildlife Zoo during December 2017. To prevent sample contamination, the bharals enclosure was cleaned in advance and none of the sample donors received antibiotic or probiotic therapy for the past three months. The collected samples were temporarily stored in dry ice (− 50 °C). All feces samples were eventually frozen and stored at − 80 °C for further analyses.

DNA extraction

Total genome DNA of samples was extracted by CTAB method using QIAamp® Stool Mini Kit (Qiagen, Germany). DNA concentration and purity were detected with 1% agarose gel. Diluted the DNA to 1 ng/μL with sterile water.

PCR amplification, purification and sequencing

16S rRNA genes of distinct regions (V4-V5) were amplified used specific primer (515F:5′-GTGCCAGCMGCCGCGGTAA-3′; 907R:5′-CCGTCAATTCCTTTGAGTTT-3′) with the barcode [41]. All PCR reactions were 30 μL systems, including 15 μL Phusion® High-Fidelity PCR Master Mix (New England Biolabs), 3 μL (2 μM) forward and reverse primers, 10 μL (10 ng) template DNA and 2 μL ddH2O. The thermal cycling consisted of pre-denatured at 98 °C for 1 min, followed by 30 cycles, including denaturation at 98 °C for 10 s, annealing at 50 °C for 30 s, and extension at 72 °C for 30 s. Finally, 72 °C for 5 min.

The 1 × loading buffer (contained SYB green) of the same volume was mixed with the PCR products and detected by 2% agarose gel electrophoresis. PCR products was mixed in equidensity ratios. Then, mixture PCR products was purified with GeneJET™ Gel Extraction Kit (Thermo Scientific).

Sequencing libraries were generated using Ion Plus Fragment Library Kit 48 rxns (Thermo Scientific). The library was sequenced on an Ion S5™ XL platform to produce 400 bp/600 bp single-end reads.

Sequence processing and statistical analysis

Get single-end reads and filtered the quality of the raw reads to obtain high-quality clean reads [42] (V1.9.1, By using UCHIME algorithm ( [43], the reads were compared with the reference database to detect chimera sequences, and then removed [44]. Then the Clean Reads finally obtained. Sequences analysis were performed by Uparse software (Uparse v7.0.1001, [45]. Sequences with ≥97% similarity were assigned to the same operational taxonomic units (OTUs). Representative sequence for each OTU was screened for further annotation [41, 46].

Qiime software (Version 1.9.1) was used to analyze the differences of alpha diversity indexes between groups. R software (Version 2.15.3) was used to plot rarefaction curves, rank abundance curves, and stacked histograms of the relative abundance. Ade4 and ggplot2 packages in R were used for PCA and PCoA analyses, respectively. LEfSe software (LEfSe 1.0) was used for LEfSe analysis, which defaulted to a filter value of 4 for the LDA Score. For MRPP analysis, the MRPP function of the vegan package for R was used.

Availability of data and materials

The datasets generated and analysed during the current study are not publicly available due to further analysis at later stage, but are available from the corresponding author on reasonable request.



Linear discriminant analysis coupled with effect size


Multi Response Permutation Procedure


Operational taxonomic unit


Principal component analysis


Principal co-ordinate analysis


Polymerase chain reaction


Ribosomal ribonucleic acid


  1. 1.

    Thaiss CA, Zmora N, Levy M, Elinav E. The microbiome and innate immunity. Nature. 2016;535(7610):65–74.

    CAS  Article  Google Scholar 

  2. 2.

    Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Human gut microbes associated with obesity. Nature. 2006;444(7122):1022–3.

    CAS  Article  Google Scholar 

  3. 3.

    Mølbak L, Johnsen K, Boye M, Jensen TK, Johansen M, Møller K, Leser TD. The microbiota of pigs influenced by diet texture and severity of Lawsonia intracellularis infection, Veterinary Microbiology. 2008:128(1–2):96–107.

    Article  Google Scholar 

  4. 4.

    Buddington RK, Sangild PT. COMPANION ANIMALS SYMPOSIUM: development of the mammalian gastrointestinal tract, the resident microbiota, and the role of diet in early life. J Anim Sci. 2011;89(5):1506–19.

    CAS  Article  Google Scholar 

  5. 5.

    Gillor O, Etzion A, Riley MA. The dual role of bacteriocins as anti- and probiotics. Appl Microbiol Biotechnol. 2008;81(4):591–606.

    CAS  Article  Google Scholar 

  6. 6.

    Schuijt TJ, van der Poll T, de Vos WM, Wiersinga WJ. The intestinal microbiota and host immune interactions in the critically ill. Trends Microbiol. 2013;21(5):221–9.

    CAS  Article  Google Scholar 

  7. 7.

    Koboziev I, Reinoso Webb C, Furr KL, Grisham MB. Role of the enteric microbiota in intestinal homeostasis and inflammation. Free Radic Biol Med. 2014;68:122–33.

    CAS  Article  Google Scholar 

  8. 8.

    Ley RE, Lozupone CA, Hamady M, Knight R, Gordon JI. Worlds within worlds: evolution of the vertebrate gut microbiota. Nat Rev Microbiol. 2008;6(10):776–88.

    CAS  Article  Google Scholar 

  9. 9.

    Schwab C, Cristescu B, Northrup JM, Stenhouse GB, Ganzle M. Diet and environment shape fecal bacterial microbiota composition and enteric pathogen load of grizzly bears. PLoS One. 2011;6(12):e27905.

    CAS  Article  Google Scholar 

  10. 10.

    Navarrete P, Magne F, Araneda C, Fuentes P, Barros L, Opazo R, Espejo R, Romero J. PCR-TTGE analysis of 16S rRNA from rainbow trout (Oncorhynchus mykiss) gut microbiota reveals host-specific communities of active bacteria. PLoS One. 2012;7(2):e31335.

    CAS  Article  Google Scholar 

  11. 11.

    Kovács F, Nagy B, Sinkovics G. The gut bacterial flora of healthy early weaned piglets, with special regard to factors influencing its composition. Acta Vet Acad Sci Hung. 1972;22(4):327.

    PubMed  Google Scholar 

  12. 12.

    McKenney PT, Pamer EG. From hype to Hope: the gut microbiota in enteric infectious disease. Cell. 2015;163(6):1326–32.

    CAS  Article  Google Scholar 

  13. 13.

    Cox LA, Olivier M, Spradlingreeves K, Karere GM, Comuzzie AG, Vandeberg JL. Nonhuman Primates and translational research-cardiovascular disease. ILAR J. 2017;58(2):235–50.

    CAS  Article  Google Scholar 

  14. 14.

    Wang XM, Li M, Tang S, Liu Z. A preliminary study of some characters of blue sheep population ecology in spring. Acta Theriologica Sinica. 1998;18(1):27–33.

    Google Scholar 

  15. 15.

    Schaller GB. Mountain monarchs. wild sheep and goats of the Himalaya. Chicago: University of Chicago Press; 1977. p. 425.

    Google Scholar 

  16. 16.

    Wang XM, Schaller GB. Status of Large Mammals in Western Inner Mongolia, China. Journal of East China Normal University (Special Issue of Mammals). 1996;6:93–104.

    Google Scholar 

  17. 17.

    Schaller GB. Wildlife of the Tibetan steppe, vol. 373. Chicago: University of Chicago Press; 1998.

    Google Scholar 

  18. 18.

    Liu CG, Zhang SW, Ren JR. Research foods and foods source about snow leopard(Panthera uncia). Journal of Shaanxi No rmal University (Natural Science Edition). 2003;31(S2):154–9.

    Google Scholar 

  19. 19.

    Cheng X, Yin X, Xia G, Yu Y, Hou X. Analysis of age structure of blue sheep population. Shaanxi Journal of Agricultural Sciences. 2014;60(08):40–2.

    Google Scholar 

  20. 20.

    Liu ZS, Wang XM, Li ZG, Cui DY, Li XQ. Seasonal variation of diurnal activity budgets by blue sheep (Pseudois nayaur) with different age-sex classes in Helan Mountain. Zool Res. 2005;26(4):350–7.

    Google Scholar 

  21. 21.

    Zhu M, Zhou C, He Y, Huang Y, Lu D, Zeng X. Exploring the phylogenetic relationships among species of tribe caprini based on mitochondrial Cyt b and nuclear gene ZFY. Acta Theriologica Sinica. 2014;34(4):366–73.

    Google Scholar 

  22. 22.

    Cao LR, Wang XM, Fang SG. A molecular phylogeny of Bharal and dwarf blue sheep based on mitochondrial cytochrome b gene sequences. Acta Zool Sin. 2003;49(2):198–204.

    CAS  Google Scholar 

  23. 23.

    Lundberg DS, Yourstone S, Mieczkowski P, Jones CD, Dangl JL. Practical innovations for high-throughput amplicon sequencing. Nat Methods. 2013;10(10):999–1002.

    CAS  Article  Google Scholar 

  24. 24.

    Sundset MA, Praesteng KE, Cann IK, Mathiesen SD, Mackie RI. Novel rumen bacterial diversity in two geographically separated sub-species of reindeer. Microb Ecol. 2007;54(3):424–38.

    Article  Google Scholar 

  25. 25.

    Gruninger RJ, Sensen CW, McAllister TA, Forster RJ. Diversity of rumen bacteria in Canadian cervids. PLoS One. 2014;9(2):e89682.

    Article  Google Scholar 

  26. 26.

    Ishaq SL, Wright AD. High-throughput DNA sequencing of the ruminal bacteria from moose (Alces alces) in Vermont, Alaska, and Norway. Microb Ecol. 2014;68(2):185–95.

    CAS  Article  Google Scholar 

  27. 27.

    Koike S, Yoshitani S, Kobayashi Y, Tanaka K. Phylogenetic analysis of fiber-associated rumen bacterial community and PCR detection of uncultured bacteria. FEMS Microbiol Lett. 2003;229(1):23–30.

    CAS  Article  Google Scholar 

  28. 28.

    Thoetkiattikul H, Mhuantong W, Laothanachareon T, Tangphatsornruang S, Pattarajinda V, Eurwilaichitr L, Champreda V. Comparative analysis of microbial profiles in cow rumen fed with different dietary fiber by tagged 16S rRNA gene pyrosequencing. Curr Microbiol. 2013;67(2):130–7.

    CAS  Article  Google Scholar 

  29. 29.

    Chen Y, Penner GB, Li M, Oba M, Guan LL. Changes in bacterial diversity associated with epithelial tissue in the beef cow rumen during the transition to a high-grain diet. Appl Environ Microbiol. 2011;77(16):5770–81.

    CAS  Article  Google Scholar 

  30. 30.

    Spence C, Wells WG, Smith CJ. Characterization of the primary starch utilization operon in the obligate anaerobe Bacteroides fragilis: regulation by carbon source and oxygen. J Bacteriol. 2006;188(13):4663–72.

    CAS  Article  Google Scholar 

  31. 31.

    Backhed F, Ding H, Wang T, Hooper LV, Koh GY, Nagy A, Semenkovich CF, Gordon JI. The gut microbiota as an environmental factor that regulates fat storage. Proc Natl Acad Sci U S A. 2004;101(44):15718–23.

    Article  Google Scholar 

  32. 32.

    Sears CL. A dynamic partnership: celebrating our gut flora. Anaerobe. 2005;11(5):247–51.

    Article  Google Scholar 

  33. 33.

    Zhao Y. Study on the correlation between antibiotic induced changes in mouse flora and host metabolism: Doctoral dissertation. Huazhong university of science and technology; 2013.

  34. 34.

    Kang C, Hui SC, Zhu XH, Mi MT: Study on the mechanism of intestinal flora in capsaicin improving the signs of obesity induced by high-fat diet in mice. In: The 13th western China nutrition and health summit BBS: 2018; Chengdu, sichuan, China; 2018: 1.

  35. 35.

    Bourquin LD, Titgemeyer EC, Fahey GC Jr. Vegetable fiber fermentation by human fecal bacteria: cell wall polysaccharide disappearance and short-chain fatty acid production during in vitro fermentation and water-holding capacity of unfermented residues. J Nutr. 1993;123(5):860–9.

    CAS  Article  Google Scholar 

  36. 36.

    Bourquin LD, Titgemeyer EC, Fahey GC. Fermentation of various dietary fiber sources by human fecal bacteria. Nutr Res. 1996;16(7):1119–31.

    Article  Google Scholar 

  37. 37.

    Bolnick DI, Snowberg LK, Hirsch PE, Lauber CL, Org E, Parks B, Lusis AJ, Knight R, Caporaso JG, Svanbäck R. Individual diet has sex-dependent effects on vertebrate gut microbiota. Nat Commun. 2015;5:4500.

    Article  Google Scholar 

  38. 38.

    Konturek PC, Haziri D, Brzozowski T, Hess T, Heyman S, Kwiecien S, Konturek SJ, Koziel J. Emerging role of fecal microbiota therapy in the treatment of gastrointestinal and extra-gastrointestinal diseases. J Physiol Pharmacol. 2015;66(4):483–91.

    CAS  PubMed  Google Scholar 

  39. 39.

    Rieger MA, Parlesak A, Pool-Zobel BL, Rechkemmer G, Bode C. A diet high in fat and meat but low in dietary fibre increases the genotoxic potential of 'faecal water'. Carcinogenesis. 1999;20(12):2311–6.

    CAS  Article  Google Scholar 

  40. 40.

    Zhang Y: Antibiotics-induced mice microbial alteration and their relationship with host metabolic profiles. Master's thesis. Shandong Normal University; 2018.

  41. 41.

    Luo M, Liu Y, Wu P, Luo DX, Sun Q, Zheng H, Hu R, Pandol SJ, Li QF, Han YP. Alternation of gut microbiota in patients with pulmonary tuberculosis. Front Physiol. 2017;8(822).

  42. 42.

    Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. Embnet Journal. 2011;17(1).

    Article  Google Scholar 

  43. 43.

    Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200.

    CAS  Article  Google Scholar 

  44. 44.

    Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, Ciulla D, Tabbaa D, Highlander SK, Sodergren E, et al. Chimeric 16S rRNA sequence formation and detection in sanger and 454-pyrosequenced PCR amplicons. Genome Res. 2011;21(3):494–504.

    CAS  Article  Google Scholar 

  45. 45.

    Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8.

    CAS  Article  Google Scholar 

  46. 46.

    Bokulich NA, Subramanian S, Faith JJ, Gevers D, Gordon JI, Knight R, Mills DA, Caporaso JG. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Methods. 2013;10(1):57–9.

    CAS  Article  Google Scholar 

Download references


We would like to express our heartfelt thanks to the director of the Yellow-River-Sources Park Management Station in Three-River-Sources National Park of Mardo county and all bharals breeders at the Qinghai-Tibet Plateau Wildlife Zoo in Xining for their active cooperation and their valuable suggestions on the collection of bharals feces samples.


This work was supported by following grants: National Key R&D Program of China (2017YFC0506405) research design; The Strategic Priority Research Program of the Chinese Academy of Sciences (XDA23060602; XDA2002030302) data collection; Construction Fund for Qinghai Key Laboratories (2017-ZJ-Y23) data analysis and interpretation; Qinghai Key R&D and Transformation Program (2019-SF-150) manuscript writing.

Author information




XWC, HMG, WQ, JRC, ZYC and TZZ contributed conception and design of the study; XWC, HMG, WQ, LW and JRC carried out the DNA extraction and data analysis; XWC, HMG, WQ, PFS, LW, and GSW participated in samples collection; HMG, GSW, WQ, PFS, LW and TZZ advised on manuscript content; XWC wrote the first draft of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Tongzuo Zhang.

Ethics declarations

Ethics approval and consent to participate

All applicable international, national and institutional guidelines for animal care and use were observed. All procedures related to animal’s research were approved by the Ethics and Welfare of Experiment Animals Committee affiliated to Northwest Institute of Plateau Biology.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chi, X., Gao, H., Wu, G. et al. Comparison of gut microbiota diversity between wild and captive bharals (Pseudois nayaur). BMC Vet Res 15, 243 (2019).

Download citation


  • Bharal (Pseudois nayaur)
  • 16S rRNA gene
  • Gut microbiota
  • Health assessment
  • High-throughput sequence analysis