Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 11332.
Published online 2017 September 12. doi:  10.1038/s41598-017-11770-4
PMCID: PMC5596057

Analysis and comparison of the wolf microbiome under different environmental factors using three different data of Next Generation Sequencing


Next Generation Sequencing has been widely used to characterize the prevalence of fecal bacteria in many different species. In this study, we attempted to employ a low-cost and high-throughput sequencing model to discern information pertaining to the wolf microbiota. It is hoped that this model will allow researchers to elucidate potential protective factors in relation to endangered wolf species. We propose three high-throughput sequencing models to reveal information pertaining to the micro-ecology of the wolf. Our analyses advised that, among the three models, more than 100,000 sequences are more appropriate to retrieve the communities’ richness and diversity of micro-ecology. In addition, the top five wolf microbiome OTUs (99%) were members of the following five phyla: Bacteroidetes, Fusobacteria, Firmicutes, Proteobacteria, and Actinobacteria. While Alloprevotella, Clostridium_sensu_stricto_1, Anaerobiospirillum, Faecalibactreium and Streptococcus were shared by all samples, their relative abundances were differentially represented between domestic dogs and other wolves. Our findings suggest that altitude, human interference, age, and climate all contribute towards the micro-ecology of the wolf. Specifically, we observed that genera Succinivibrio and Turicibacter are significantly related to altitude and human interference (including hunting practices).


A large number of bacterial species are known to colonize various anatomical sites within the body13. For instance, the colon is estimated to contain between 1010 and 1014 bacterial cells4. These extremely large microbial populations display significant diversity and have evolved novel mechanisms that facilitate their proliferation and maintenance5. The associated mechanisms have been the focus of a large number of studies worldwide6. Many studies have demonstrated that microbiota play important roles in the biogeochemical cycling of carbon, nitrogen, and phosphorus. Furthermore, these microbiological communities are capable of facilitating the decomposition of organic material and the extraction of nutrients from the resultant matter.

The focus of this research involved the microbiota of the native wolf with a particular emphasis on microbial diversity. It has been reported that microbial diversity plays an important role in the maintenance of host health7. It has also been suggested that factors such as age and gender affect the composition of microbial populations8. It is likely that mutualism between -specific characteristics, microbial richness, and diversity is an important phenomenon in host health9.

The traditional method that has been used to study microbial diversity in animals involves morphological and biochemical examination of colonizing microorganisms following in vitro isolation and purification techniques10. However, this method has many inherent drawbacks and limitations. For instance the vast majority of fecal microbes cannot be cultured. Furthermore, only 0.1%–10% of microbial species can be differentiated using these methodologies11. Another method that can be used to analyze microbial diversity depends on traditional molecular biology techniques, including denaturing gradient gel electrophoresis (DGGE)12, temperature gradient gel electrophoresis (TGGE)13, restriction fragment length polymorphism (RFLP) analysis14, and terminal restriction fragment length polymorphism (T-RFLP) analysis15. These methods do not exhibit many of the drawbacks associated with morphological and biochemical methods. Each of these approaches has contributed to our understanding of the importance of microbial diversity in hosts; however, all of these approaches have limitations in relation to both the quantity and quality of the data generated.

In recent years, a cross-disciplinary approach that has utilized both bioinformatics and molecular biological techniques is being adopted to study microbial diversity with dedicated academic programs featuring utilization of associated techniques. DNA sequencing technologies and applications have evolved extremely rapidly and the associated platforms currently facilitate the assimilation of big genome data, which are crucial for many research areas and applications16. From Sanger sequencing technologies in the mid-twentieth century to the utilization of High-Throughput sequencing platforms at the beginning of the twenty-first century, sequencing analyses have helped to develop our understanding of fecal bacterial profiles. Indeed, Sanger sequencing technologies have helped scientists to overcome several of the limitations pertaining to more traditional approaches that permit bacterial identification17. Furthermore, high-throughput sequencing technologies have provided quicker turnaround times and increased sequence data assembly rates. The latter technologies also facilitate low-cost utilization, higher coverage rates, higher accuracy, increasing read-lengths, and paired-end sequencing18. Thus, Next Generation Sequencing (NGS) platforms are now more widely applied in the analysis of fecal bacteria1922. These analyses favor low-cost, high-throughput methodologies; however, the associated cost and throughput is somewhat dependent on the application requirements. In this article, we analyzed a substantial body of High-Throughput sequencing data to determine an appropriate balance that meets high-throughput and low-cost requirements.

Previous studies have demonstrated the existence of microbiota diversity and richness in animals including the panda and the monkey2328. These studies investigated the latter parameters in both species and subspecies29, 30. Several studies pertaining to both the wolf and the dog have also been conducted in this field. Zhang et al. reported data relating to microbiota diversity and richness in the wolf following an analysis that facilitated the cloning of bacterial 16S rRNA gene amplicons31. Suchodolski et al. showed microbiota diversity in the intestinal segments of dogs32, 33. Additional studies reported data pertaining to microbiomes of ill-conditioned fecal32, 34, 35. However, the methods that were used in these analyses and the associated data sizes that were generated are outdated. Furthermore, there is no reference to the dog microbiome in these wolf-specific studies. The dog is a subspecies of the wolf; thus, to attain a more generic understanding of microbiota diversity and richness in related subspecies, samples from both the dog and the wolf should be studied together.

As part of this study, eighteen dogs and wolves were selected to characterize microbiota composition using multi-group sequencing analysis. The wolf is now listed on the International Union for Conservation of Nature and Natural Resources endangered species list of threatened species36. Over the last number of centuries, the wolf has been one of the most widely distributed animal species in the world31, 37. However, the population of this species has rapidly declined and is currently threatened by habitat loss and hunting practices. For all this, scientists will continue to work hard to explore the protection mechanism38. Wolves are independent and are able to adapt to the environment. However, its independent and adaptation mechanisms are unknown. At present only semi-artificial, semi-wild, and bowel disease is a significant aspect of health. All of this needs to be solved. Intestinal microbes are part of a complex ecosystem. They have a mutual relationship with the host and play an essential role in maintaining the host’s health. Following a comprehensive analysis of multi-group sequencing data, this report reveals the composition of the wolf fecal microbiota while also describing factors that influence this composition. This study will provide valuable basic data that should help in efforts related to the future conservation of wolves.


OTUs and Taxonomic composition of the fecal microbiomes

From the three data sets, we obtained >50000, >100000 and >150000 raw sequences, across a total of 18 samples, respectively (Table 1). To avoid analytical variation, identical protocols, including the mothur MiSeq-SOP, were utilized. Following denoising steps, effective and unique sequences were obtained. The average numbers of unique sequences were 1916, 3339, and 4977 (Table 1). Clustering was performed using Uparse, and a total of 118, 151 and 170 OTUs were generated for the three data sets, respectively (Table 1). To explore the dominant bacterial species, the results were annotated using GraPhlAn (Fig. 1A)39.

Table 1
Comparison of phylotype coverage and diversity estimation of the 16S rRNA gene libraries at 3% dissimilarity from the analysis of the basic statistics information.
Figure 1
(A) The circle from inside to outside in turn was on behalf of different classification levels, and species abundance is proportional to the size of the circle, the different colors represent different phyla, and abundance of the top 40 species was showed ...

Associated histograms were generated with species annotation to six classification levels. Twenty phyla were observed following this analysis. The five most prominent OTUs (99%) were from the following five phyla; Bacteroidetes, Fusobacteria, Firmicutes, Proteobacteria, and Actinobacteria. This observation was consistent among the three data sets (Fig. 1E). Upon comparison of the associated phyla, the relative abundance of Bacteroidetes was as follows: WJ > WN > WH > DN; the relative abundance of Fusobacteria was as follows: WN > DN > WH > WJ; the relative abundance of Firmicutes was as follows: WH > WN > DN > WJ; the relative abundance of Proteobacteria was as follows: DN > WJ > WH > WN; and the relative abundance of Actinobacteria was as follows: WJ > WH > DN > WN. Upon analysis of climate-mediated effects in relation to microbiota composition, we observed positive and negative correlations with respect to climate for Bacteroidetes and Firmicutes, respectively. Proteobacteria and Actinobacteria were observed to positively correlate with wolves exposed to human interference practices, while Fusobacteria were observed to be negatively correlated with these interference effects. Interestingly, Cyanobacteria were only observed in the WN, and Verrucomicrobia were observed exclusively in the DN. At the genus level, the most frequently detected genera were Bacteroides, Alloprevotella, Sutterella, Clostridium_sensu_stricto_1, Anaerobiospirillum, Prevotellaceae_Ga6A1_group, Helicobacter, Faecalibacterium, Phascolarctobacterium, and Lachnoclostridium. In addition, Clostridium_sensu_stricto_1 was most prevalent in the WN group and Phascolarctobacterium and Lachnoclostridium were the most prevalent genera in the WH group.

Comparison of microbial group diversity

Following alpha diversity analysis, the indices for bacterial richness and the diversity of OTUs at a 3% sequence dissimilarity level are summarized in Table 1. Increased community richness was observed following the analysis of increased amounts of data (Table 1). The richness index values for SET 2 and SET 3 were almost equal and were greater than the associated value for SET 1 (Table 1). In SET 3, there was no significant difference in community richness between the different samples; however, the Shannon and Simpson indices for wolves in Inner Mongolia were reduced compared with other areas (Table 1). Following analysis of the rarefaction curve (Fig. 1B) and species accumulation boxplot (Fig. 1C) diagram, it is noticeable that the curves are on a trajectory towards a constant equilibrium.

We are told from two model OTU diagrams (linear model PCA and nonlinear model NMDS) (Fig. 1D) that the samples from one groups are collected together, so the difference and consistency were showed phenomenally. In order to elucidate further the cluster structures for wolves, we attempted to map UPGMA hierarchical clustering for SET 3 (Fig. 1E). However, the cluster analysis generated following Weighted Unifrac and Unweighted Unifrac analyses was more complex (Fig. 1E). The heatmap from the Spearman analysis suggested that Succinvibrio richness negatively correlated with altitude and positively correlated with pressures (Fig. 2). Conversely, Turicibacter richness positively correlated with altitude and negatively correlated with human interference (Fig. 2). We also observed that the prevalence of some bacterial phyla correlated both negatively and positively with age. Additional details relating to observed correlations are presented in Fig. 2.

Figure 2
The ordinate is the environment information and the abscissa is species. The corresponding values of heatmap are spearman rank relational coefficient. The * means significant (P < 0.05) and the ** means remarkable significant (P < 0.01). ...

Optimization selection of the three models

Utilizing identical protocols, we get more diversity index (ace, chao1, simpson, shonnon) from three models. We realize of visualizing processing to the basic statistics data, such as minimum, maximum, median, averages and so on (Fig. 3A). From Fig. 3A, the community richness (ace, chao1) increased with the data increases. However, the community diversity is different from the patterns of community richness. The community diversity of SET 1 has significantly different from the SET2 and SET3. Among them, the median, average and numeric range is similar in SET2 and SET3 and is higher than the SET1. In other words, the development of SET2 and SET3 has a similar stable trend. In particular, the SET3 can reflect the outliers. In the Goodness of fit on three data set of analysis, we added the scatter plot and regression curve. Our normalization shows the probability distributions. The data of SET 1 is decentralized but the data of SET 2 and SET 3 were more centrally. Moreover, the data of SET 2 and SET 3 were almost unanimously and tend to be stable. On the other hand, SET 2 and SET 3 reflect the community diversity is largely consistent (Shannon: R2 = 99.5%, Simpson: R2 = 99.6%) (Fig. 3B).

Figure 3
(A): Boxplot with whiskers from minimum to maximum; (B): scatter diagram and Goodness of Fit.

Comparison of the gut microbiota of domestic dog and wolf from the SET 3

We characterized the gut microbiota of 14 wolves and 4 domestic dogs. The wolves’ number of OTUs is 167 and the domestic dogs’ number of OTUs is 170. Five phyla including Bacteroidetes, Fusobacteria, Firmicutes, Proteobacteria, Actinobacteria were predominant bacterial. In addition, Cyanobacteria (0.1391%) was only detected in wolf and Verrucomicrobia (0.2312%) only detected in domestic dog. Particular phyla were at low relative abundances. We also use GraPhlAn (Graphical Phylogenetic Analysis), a computational tool that produces high-quality, compact visualizations of microbial metagenomes(Fig. 4)39. Exploiting the shared and unique bacterial taxa between the gut microbiota of the domestic dog and the wolf is also our aim. It was unexpected that Cyanobacteria exist only in the wolf and Verrucomicrobia exists only in the domestic dog. We used linear discriminant analysis effect size (LEfSe) to identify genus differentially represented between the wolf and the domestic dog. While Alloprevotella, Clostridium_sensu_stricto_1, Anaerobiospirillum, Faecalibactreium and Streptococcus were shared by all samples, their relative abundances were differentially represented between the two (Fig. 5). The genus Alloprevotella and Clostridium_sensu_stricto_1 is significantly higher in the wolf than in the domestic dog. In contrast, Anaerobiospirillum, Faecalibactreium and Streptococcus is more abundant in the domestic dog (Fig. 5). Streptococcus only exists in the domestic dog and was absent from the wolf. Whether or not the gut microbiota is involved within digestion needs further investigation.

Figure 4
GraPhlAn visualization of annotated phylogenies and taxonomies. We comprise of microbial community abundances between dog and wolf using the phylogenetic tree on all available microbiota. Colors and background annotation highlight bacterial phyla.
Figure 5
Bacteria differentially represented between the domestic dogs and the wolves identified by linear discriminant analysis coupled with effect size (LEfSe). Histogram showing Bacteria that are more abundant in the domestic dogs (red color) or the wolves ...


The objective of this study was to determine an appropriate balance between sequence data generation and cost in relation to an analytical model that could help to protect the valuable wildlife-wolf. The study explored the composition of microbiota of wolves following high-throughput sequencing analysis.

With the rapid development of Next Generation Sequencing (NGS) technology, there is an onus on micro-ecologists to use high-throughput sequencing platforms to help solve complex biological problems. However, high-throughput sequencing can be relatively costly. Nevertheless, high-throughput sequencing generates accurate results that facilitate complex sequence analysis. Upon analysis of these results, we advised that, >100,000 sequences can restore the information pertaining to community richness and micro-ecological diversity.

Traditional paradigms suggest that the generation of greater amounts of sequencing data will facilitate more accurate results. Thus, we chose the SET3 data to explore the microbiota of wolves. Previous work in this field performed using phylogenetic analysis of 16S rRNA gene sequences revealed distal fecal bacterial diversity in wild wolves and dogs4042. The latter studies on bacterial diversity were predominantly based on cloning bacterial 16S rRNA gene amplicons. This method was deemed the most advanced for these analysis types; however, there are limitations associated with this form of analysis including limitations relating to the amount of data generated. Moreover, the gray wolf (Canis lupus) encompasses both wolf (Canis lupus linnaeus) and dog (Canis lupus familiaris) species. Unfortunately, these studies only revealed information regarding singular dog or wolf species. Therefore, we attempted to revisit this area of research using Next Generation Sequencing (NGS) technology and the Illumina MiSeq (Illumina MiSeq, USA) platform. Our results demonstrated that Bacteroidetes (46.48%), Fusobacteria (30.54%), Firmicutes (13.46%), Proteobacteria (8.83%), Actinobacteria (0.53%), Cyanobacteria (0.10%), Verrucomicrobia (0.01%), Tenericutes (<0.01%), Saccharibacteria (<0.01%) and Lentisphaerae (<0.01%) were the most prevalent phyla in wolf microbiota. However, Chen et al. identified the presence of five phyla (Firmicutes (60%), Bacteroidetes (16.9%), Proteobacteria (9.2%), Fusobacteria (9.2%) and Actinobacteria (4.6%)) and other reports only detected seven bacterial phyla including the additional phyla Spirochaetes and Tenericutes33, 34, 41, 43, 44.

Upon analysis of different variables, we observed that microbiota composition in wolves is dependent upon a number of factors including age, altitude, pressure, and climate. Another interesting observation included the fact that one of the groups, DN, exhibited several noticeable differences from the other groups. This occurrence is most likely because the associated group consisted of dogs, which are a subspecies of wolf. Following on from our analysis and contributions from other studies pertaining to affecting factors2, 3, 7, 4557, it is apparent that the micro-ecological environment of wolves is responsible for defense against unfavorable environmental factors. For instance, Succinivibrio, Turicibacter, and Prevotellaceae_Ga6A1_group appear to be involved in protection against pressures associated with human interference of wolf species (reference Spearman analysis diagram, Fig. 2).

In order to investigate the relationship between bacterial populations and environmental factors, we performed a literature search focusing on specific bacterial genera. The Spearman’s rank correlation for the researched bacterial genera was greater than 0.7 and was extremely significant (P < 0.01). Succinivibrio and Turicibacter were observed to correlate with altitude and human interference. Members of the genus Succinivibrio require carbon dioxide for growth. These bacteria are anaerobic and ferment organic matter produced by the Krebs cycle to generate acetic acid and succinic acid. Turicibacter species are facultative. In addition, they produce acid but not gas. Thus, it is likely that a greater abundance of Succinivibrio and Turicibacter can facilitate Krebs Cycle progression, thereby resulting in greater absorption of carbon dioxide. Therefore, as environmental pressure increases, fecal bacteria can help to reduce the amount of energy produced, thereby maintaining homeostasis in the body. These findings relating to the prevalence of Succinivibrio and Turicibacter in the microbiota of wolves may help to explain how this energy conservation manifests itself under adverse conditions.

Gut microbiota of wolf is negatively associated with the pressure from humans while domestic dogs are just the opposite. The wolves are not interested in human social cues and will be a threat. Conversely, the domestic dogs care about the human social cues. After thirty-thousand years of domesticating, the dog has gradually understood and adapted to humans while the wolf continued in a state of enmity. In order to satisfy the different relations, the organisms may make some changes, which include the gut microbial.

In summary, the three data set determined an appropriate balance between sequence data generation and cost for the primary scientific research workers. The findings obtained in our study also provided a special insight into the ecology and biodiversity of the wolf gut microbiome. We observed that the microbiota composition was dependent on habitat and four different analyzed factors were important in determining the prevalence of microbiological genera in the wolf fecal. These microbial communities co-exist with their host and play important roles in the long-term evolution of the host3. These findings provided a powerful tool for characterization of the micro-ecological environment of threatened wildlife species, thereby allowing us to identify factors that might be important in population maintenance and protection.

Our samples in this study were taken from several representative geographical areas in China. We all know that wolves are globally distributed and inhabited in different ecological environment. Further research is necessary to test our findings in more wolves, as well as to better understand the trade-off between nutrition and health via shifts in gut microbiota composition. The ultimate goal is that all wild animals around us can survive and reproduce better.


Sample collection

Fecal samples were collected from wolves in China. They were all raised semi-freely in the Care Centre with raw meat and water before the fecal samples were collected. None of the wolves received any treatment (e.g., antibiotic therapy) that would be expected to have an impact on the composition of the intestinal microbial community. All procedures that were performed on animals were conducted in accordance with the ethical standards of the Qufu Normal University Animal Care and Use Committee (Permit Number: QFNU2015-002). None of the animals were harmed during the collection of fecal samples. Most scientists divide the wolves in China into five subspecies and forms: Canis lupus desertorum Bogdanow, C. I. filchneri Matschie, C. I. chanco Gray, Inner Mongolia (eastern part), C. I. Nei-Mongol form (western and mid part) and C. I. South China form. In the economically developed eastern China, the human disturbance is very powerful. The result is that the wolves’ population has declined drastically. Especially in the southern provinces, it is not clear whether wolves recorded in the southern provinces represent permanent populations, or a steady stream of individuals migrating from the northern provinces5860. However, the extreme conditions would have a serious impact on body function and gut microbiome. So, our study collected three subspecies and forms: C. I. chanco Gray (WH), C. I. Nei-Mongol form (WN and DN) and C. I. South China form (WJ). These environments are less extreme. Details pertaining to the animals chosen are presented in Table 2. Different species-specific territories exhibit different climates and the three regions chosen as part of this analysis represented the different climates in China. In China, rainfall is one of the most important climate-mediated effects and we used rainfall as an indicator of climate. From previous studies, we hypothesized that human interference affects nervous system and sample microbiome composition as well6164. Thus, pressures exerted by human interference were investigated for their effects on the micro-ecology of the dog and the wolf. The fecal matter of one sample is collected in triplicate in three days. Fecal collections were immediately made after defecation during the early morning. And fecal samples were immediately placed into sterile plastic tubes frozen in sample containers and carried to the lab and stored at −80 °C until DNA extraction.

Table 2
Details of the animal information.

DNA extraction, PCR amplification, and 16S rRNA sequencing

Total community genomic DNA was extracted from all fecal samples using a QIAamp DNA Stool Mini Kit (Qiagen, Germany) as per the manufacturer’s instructions. The DNA was amplified using specific primers that targeted the V3-V4 region of the 16S rRNA bacterial gene. The primers also carried the Illumian MiSeq sequencing adapter (16S Amplicon PCR Forward Primer: CTACGGGNGGCWGCAG and 16S Amplicon PCR Reverse Primer: GACTACHVGGGTATCTAATCC)7. The PCR mix was prepared using the KAPA HiFi Hot Start Ready Mix (2×) (TaKaRa Bio Inc., Japan). PCR conditions for MiSeq are described in Wu et al. (2016). Amplicons of 16S rDNA were purified using AMPure XP beads (Beckman, USA). The final DNA samples that were extracted from the fecal samples were pooled in equal concentrations prior to sequencing on the Illumina MiSeq platform (Illumina MiSeq sequencing system, USA) in our laboratory.

Data selection and arrangement

Greater than 150,000 reads per sample were generated using the Illumina MiSeq platform. In order to assess an optimal balance between throughput and cost, we varied the numbers of reads per sample. The variations, which were randomly selected, included 50,000 reads per sample, 100,000 reads per sample, and 150,000 reads per sample, respectively. During sequencing, single DNA molecules are randomly bound to the surface of the flow cell and bridge-amplified to form clusters. Reads in the fastq file are subsequently randomly generated. The head command line subsequently permits selection of the number of sequences using the Linux operation system. The three data sets represent three different data sizes (30%, 60%, and 100%). The different sequencing data sizes represent differing sequencing depths. Researchers in the area of intestinal micro-ecology predominantly use between 20,000 and 100,000 sequencing reads3, 7, 16, 31, 51, 65. Thus, our three data sets expand the upper range and encompass much of the conventional range of use.

OTUs and fecal bacteria

First, operational taxonomical units (OTUs) were analyzed for each sample with a 97% sequence similarity cutoff value. Secondly, a summary of all taxonomic information was generated using RDP Classifier version 2.266. The phylogenetic relationship was elucidated using GraPhlAn. Finally, to standardize results, the lowest number of sequences from each sample was randomly selected and different data sets were observed to contain different homogeneous sequences.

Diversity analysis

Alpha diversity analysis facilitated the construction of a rarefaction curve and species accumulation boxplot. These were used to describe the number of OTUs and species as a function of sampling effort7. Next, richness of the associated communities was compared based on the ACE estimator and the Chao1 estimator. Community diversity was subsequently analyzed using both the Simpson index and the Shannon index.

Beta diversity analysis was used to determine microbiota composition diversity between the individuals using the linear Principal Component Analysis (PCA) model and the nonlinear Non-Metric Multi-Dimensional Scaling (NMDS) model. Wolf microbiota phylogenetic analysis was performed using an Unweighted Pair-group Method with Arithmetic Mean (UPGMA) and the associated phylogenetic trees were based on the Weighted Unifrac and Unweighted Unifrac values.

To understand the correlation between parameters including attitude, pressure, climate and age, the Spearman correlation and the Mantel test correlation were calculated. LDA Effect Size (LEfSe) can search for a Metagenomic biomarker between the two groups and the biomarker is statistically significant (P < 0.05).

We analyzed the three data sets using the Mothur (Version 1.36.1)67, Qiime (Version 1.7.0)68, Uparse (Uparse v7.0.1001)69, PyNAST (Version 1.2)70 and the following disgrams were made by R (R version 3.3.1)71.

Statistical analysis of three models

For statistical analysis of the three models, we use identical protocols to calculate four kinds of diversity indices. Using the boxplot, we performed the preliminary observation on some basic statistics. The goodness of fit describes how well it fits a set of observations. A scatter plot can suggest various kinds of correlations between variables and analyzing join level for regression model. At last, we make an optimization selection. All analysis was carried out using R (R version 3.3.1).

Availability of data and materials

We upload our raw sequences about this research on the Sequence Read Archive (accession number SRP089855).


The authors would like to thank the two anonymous reviewers and the associate editor for their valuable comments and suggestions that have helped improve the manuscript greatly. We also are thankful to Professor Katie Amo, Skidmore College, for polishing the manuscript. This work was supported in part by grants from the Special Fund for Forest Scientific Research in the Public Welfare (NO. 201404420), the National Natural Science Fund of China (NO. 31672313, NO. 31372220), the Natural Science Fund of Shandong (2013GSF11707). The authors thank all of supports.

Author Contributions

Author Contributions

H.Z., X.T. and X.W. designed and conducted the research. X.W. and J.C. performed the statistical analysis, analyzed the data and wrote the paper. X.W. and H.Z. was responsible for the DNA extraction. X.W. S.S. Y.C. and J. Y. were involved in sequencing and processing the microbial information. All authors read and approved the final manuscript.


Competing Interests

The authors declare that they have no competing interests.


Xiaoyang Wu and Huanxin Zhang contributed equally to this work.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contributor Information

Xuexi Tang, ude.cuo@xxgnat.

Honghai Zhang, moc.621@76iahgnohgnahz.


1. Dethlefsen L, McFall-Ngai M, Relman DA. An ecological and evolutionary perspective on human-microbe mutualism and disease. Nature. 2007;449:811–818. doi: 10.1038/nature06245. [PubMed] [Cross Ref]
2. Proctor LM. The National Institutes of Health Human Microbiome Project. Seminars in fetal & neonatal medicine. 2016
3. Ling Z, et al. Pyrosequencing analysis of the human microbiota of healthy Chinese undergraduates. BMC Genomics. 2013;14:63–70. doi: 10.1186/1471-2164-14-63. [PMC free article] [PubMed] [Cross Ref]
4. Pace NR. A molecular view of microbial diversity and the biosphere. Science (New York, N.Y.) 1997;276:734–740. doi: 10.1126/science.276.5313.734. [PubMed] [Cross Ref]
5. Chan. YK, et al. High fat diet induced atherosclerosis is accompanied with low colonic bacterial diversity and altered abundances that correlates with plaque size, plasma A-FABP and cholesterol: a pilot study of high fat diet and its intervention with Lactobacillus rhamnosus GG (LGG) or telmisartan in ApoE mice. BMC microbiology. 2016;16:264. doi: 10.1186/s12866-016-0883-4. [PMC free article] [PubMed] [Cross Ref]
6. Yang X, Xie L, Li Y, Wei C. More than 9,000,000 unique genes in human gut bacterial community: estimating gene numbers inside a human body. PloS one. 2009;4:e6074. doi: 10.1371/journal.pone.0006074. [PMC free article] [PubMed] [Cross Ref]
7. Wu X, et al. Comparison of the fecal microbiota of dholes high-throughput Illumina sequencing of the V3–V4 region of the 16S rRNA gene. Applied Microbiology & Biotechnology. 2016;100:1–10. doi: 10.1007/s00253-015-7019-x. [PubMed] [Cross Ref]
8. Bolnick DI, et al. Individual diet has sex-dependent effects on vertebrate gut microbiota. Nature communications. 2014;5:4500–4500. doi: 10.1038/ncomms5500. [PMC free article] [PubMed] [Cross Ref]
9. Ley RE, et al. Evolution of mammals and their gut microbes. Science. 2008;320:1647. doi: 10.1126/science.1155725. [PMC free article] [PubMed] [Cross Ref]
10. Martín V, et al. Sharing of bacterial strains between breast milk and infant feces. Journal of Human Lactation Official Journal of International Lactation Consultant Association. 2012;28:36–44. doi: 10.1177/0890334411424729. [PubMed] [Cross Ref]
11. Ringø E, Sperstad S, Kraugerud OF, Krogdahl Å. Use of 16S rRNA gene sequencing analysis to characterize culturable intestinal bacteria in Atlantic salmon (Salmo salar) fed diets with cellulose or non-starch polysaccharides from soy. Aquaculture Research. 2008;39:1087–1100. doi: 10.1111/j.1365-2109.2008.01972.x. [Cross Ref]
12. Wintzingerode FV, Göbel UB, Stackebrandt E. Determination of microbial diversity in environmental samples: pitfalls of PCR-based rRNA analysis. Fems Microbiology Reviews. 1997;21:213–229. doi: 10.1111/j.1574-6976.1997.tb00351.x. [PubMed] [Cross Ref]
13. Teske A, Wawer C, Muyzer G, Ramsing NB. Distribution of sulfate-reducing bacteria in a stratified fjord (Mariager Fjord, Denmark) as evaluated by most-probable-number counts and denaturing gradient gel electrophoresis of PCR-amplified ribosomal DNA fragments. Applied & Environmental Microbiology. 1996;62:1405–1415. [PMC free article] [PubMed]
14. Lukow T, Dunfield PF, Liesack W. Use of the T-RFLP technique to assess spatial and temporal changes in the bacterial community structure within an agricultural soil planted with transgenic and non-transgenic potato plants. FEMS microbiology ecology. 2000;32:241–247. doi: 10.1111/j.1574-6941.2000.tb00717.x. [PubMed] [Cross Ref]
15. Leckie SE. Methods of microbial community profiling and their application to forest soils. Forest Ecology & Management. 2005;220:88–106. doi: 10.1016/j.foreco.2005.08.007. [Cross Ref]
16. Liu L, et al. Comparison of next-generation sequencing systems. Journal of Biomedicine & Biotechnology. 2012;2012:355–355.
17. Eckburg PB, et al. Diversity of the human intestinal microbial flora. Science. 2005;308:1635–1638. doi: 10.1126/science.1110591. [PMC free article] [PubMed] [Cross Ref]
18. Mardis ER. Next-generation DNA sequencing methods. Annual Review of Genomics & Human Genetics. 2008;9:387–402. doi: 10.1146/annurev.genom.9.081307.164359. [PubMed] [Cross Ref]
19. Chen J, et al. Characterization of the gut microbiota in the golden takin (Budorcas taxicolor bedfordi) Amb Express. 2017;7:81. doi: 10.1186/s13568-017-0374-5. [PMC free article] [PubMed] [Cross Ref]
20. Chen J, Shang S, Wu X, Zhong H, Zhang H. The Microbial Community in the Feces of Cape Oryx (Oryx gazella) as Determined by Highthroughput Illumina Sequencing Technology. Pakistan J. Zool. 2017;49:0–9.
21. Schuster SC. Next-generation sequencing transforms today's biology. Nature Methods. 2008;5:16–18. doi: 10.1038/nmeth1156. [PubMed] [Cross Ref]
22. Tun HM, et al. Gene-centric metagenomics analysis of feline intestinal microbiome using 454 junior pyrosequencing. Journal of Microbiological Methods. 2012;88:369–376. doi: 10.1016/j.mimet.2012.01.001. [PubMed] [Cross Ref]
23. Wu CF, et al. Study of fecal bacterial diversity in Yunnan snub-nosed monkey (Rhinopithecus bieti) using phylogenetic analysis of cloned 16S rRNA gene sequences. African Journal of Biotechnology. 2010;9:6278–6289.
24. Mckenna P, et al. The Macaque Gut Microbiome in Health, Lentiviral Infection, and Chronic Enterocolitis. Plos Pathogens. 2008;4:20. doi: 10.1371/journal.ppat.0040020. [PMC free article] [PubMed] [Cross Ref]
25. Wei F, Wang X, Wu Q. The giant panda gut microbiome. Trends in Microbiology. 2015;23:450–452. doi: 10.1016/j.tim.2015.06.004. [PubMed] [Cross Ref]
26. Zhu L, Wu Q, Dai J, Zhang S, Wei F. Evidence of cellulose metabolism by the giant panda gut microbiome. Proceedings of the National Academy of Sciences of the United States of America. 2011;108:17714–17719. doi: 10.1073/pnas.1017956108. [PubMed] [Cross Ref]
27. Tun HM, et al. Microbial diversity and evidence of novel homoacetogens in the gut of both geriatric and adult giant pandas (Ailuropoda melanoleuca) Plos One. 2014;9:e79902. doi: 10.1371/journal.pone.0079902. [PMC free article] [PubMed] [Cross Ref]
28. Chan YK, et al. High fat diet induced atherosclerosis is accompanied with low colonic bacterial diversity and altered abundances that correlates with plaque size, plasma A-FABP and cholesterol: a pilot study of high fat diet and its intervention with Lactobacillus rhamnosus GG (LGG) or telmisartan in ApoE−/− mice. BMC microbiology. 2016;16:264. doi: 10.1186/s12866-016-0883-4. [PMC free article] [PubMed] [Cross Ref]
29. Bian G, Ma L, Su Y, Zhu W. The Microbial Community in the Feces of the White Rhinoceros (Ceratotherium simum) as Determined by Barcoded Pyrosequencing Analysis. PloS one. 2013;8:e70103–e70103. doi: 10.1371/journal.pone.0070103. [PMC free article] [PubMed] [Cross Ref]
30. Yu L, et al. Genomic analysis of snub-nosed monkeys (Rhinopithecus) identifies genes and processes related to high-altitude adaptation. Nature Genetics. 2016;48:947. doi: 10.1038/ng.3615. [PubMed] [Cross Ref]
31. Zhang H, Lei C. Phylogenetic analysis of 16S rRNA gene sequences reveals distal gut bacterial diversity in wild wolves (Canis lupus) Molecular Biology Reports. 2010;37:4013–4022. doi: 10.1007/s11033-010-0060-z. [PubMed] [Cross Ref]
32. Suchodolski JS, Dowd SE, Wilke V, Steiner JM, Jergens AE. 16S rRNA Gene Pyrosequencing Reveals Bacterial Dysbiosis in the Duodenum of Dogs with Idiopathic Inflammatory Bowel Disease. PloS one. 2012;7:e39333. doi: 10.1371/journal.pone.0039333. [PMC free article] [PubMed] [Cross Ref]
33. Suchodolski JS, Camacho J, Steiner JM. Analysis of bacterial diversity in the canine duodenum, jejunum, ileum, and colon by comparative 16S rRNA gene analysis. FEMS Microbiol Ecol. 2008;66:567–578. doi: 10.1111/j.1574-6941.2008.00521.x. [PubMed] [Cross Ref]
34. Suchodolski JS, et al. The effect of the macrolide antibiotic tylosin on microbial diversity in the canine small intestine as demonstrated by massive parallel 16S rRNA gene sequencing. BMC microbiology. 2009;9:210. doi: 10.1186/1471-2180-9-210. [PMC free article] [PubMed] [Cross Ref]
35. Minamoto Y, Dhanani N, Markel ME, Steiner JM, Suchodolski JS. Prevalence of Clostridium perfringens, Clostridium perfringens enterotoxin and dysbiosis in fecal samples of dogs with diarrhea. Veterinary Microbiology. 2014;174:463–473. doi: 10.1016/j.vetmic.2014.10.005. [PubMed] [Cross Ref]
36. Mech, L. D. & Boitani, L. The IUCN Red List of Threatened Species 2010: e.T3746A10049204. IUCN SSC Wolf Specialist Group (2017).
37. Kaminski, J. D (Canis familiaris) are Adapted to Receive Human Communication. (Springer Berlin Heidelberg, 2009).
38. Shang, S. et al. The repertoire of bitter taste receptor genes in canids. Amino Acids, 1–9 (2017). [PubMed]
39. Asnicar F, Weingart G, Tickle TL, Huttenhower C, Segata N. Compact graphical representation of phylogenetic data and metadata with GraPhlAn. Peerj. 2015;3:77–78. doi: 10.7717/peerj.1029. [PMC free article] [PubMed] [Cross Ref]
40. Zhang H, Chen L. Phylogenetic analysis of 16S rRNA gene sequences reveals distal gut bacterial diversity in wild wolves (Canis lupus) Molecular biology reports. 2010;37:4013–4022. doi: 10.1007/s11033-010-0060-z. [PubMed] [Cross Ref]
41. Suchodolski JS. Companion animal’s symposium: microbes and gastrointestinal health of dogs and cats. Journal of animal science. 2011;89:1520–1530. doi: 10.2527/jas.2010-3377. [PubMed] [Cross Ref]
42. Middelbos IS, et al. Phylogenetic characterization of fecal microbial communities of dogs fed diets with or without supplemental dietary fiber using 454 pyrosequencing. PloS one. 2010;5:e9768. doi: 10.1371/journal.pone.0009768. [PMC free article] [PubMed] [Cross Ref]
43. Nutrition, N. C. C. O. A. Nutrient requirements of dogs and cats. (National Academies Press, 2003).
44. Smith J, Liles C. Assessment of the canine intestinal microflora using molecular methods and serum markers. Texas A & M University. 2007;16:662–671.
45. Andersson AF, et al. Comparative analysis of human gut microbiota by barcoded pyrosequencing. PloS one. 2008;3:e2836. doi: 10.1371/journal.pone.0002836. [PMC free article] [PubMed] [Cross Ref]
46. Wei-Lin, et al. Application of metagenomics in the human gut microbiome. World journal of gastroenterology. 2015;21:803–814. doi: 10.3748/wjg.v21.i3.803. [PMC free article] [PubMed] [Cross Ref]
47. Kinumaki A, et al. Characterization of the gut microbiota of Kawasaki disease patients by metagenomic analysis. Frontiers in microbiology. 2015;6:824. doi: 10.3389/fmicb.2015.00824. [PMC free article] [PubMed] [Cross Ref]
48. Xu Z, Knight R. Dietary effects on human gut microbiome diversity. The British journal of nutrition. 2015;113(Suppl):S1–5. doi: 10.1017/S0007114514004127. [PMC free article] [PubMed] [Cross Ref]
49. Davidovics ZH, et al. The Fecal Microbiome in Pediatric Patients with Short Bowel Syndrome. JPEN. Journal of parenteral and enteral nutrition. 2015 [PMC free article] [PubMed]
50. Llorente C, Schnabl B. The gut microbiota and liver disease. Cellular and molecular gastroenterology and hepatology. 2015;1:275–284. doi: 10.1016/j.jcmgh.2015.04.003. [PMC free article] [PubMed] [Cross Ref]
51. Indugu, N., Bittinger, K., Kumar, S., Vecchiarelli, B. & Pitta, D. A comparison of rumen microbial profiles in dairy cows as retrieved by 454 Roche and Ion Torrent (PGM) sequencing platforms. Peerj4 (2016). [PMC free article] [PubMed]
52. Shakya M, et al. Comparative metagenomic and rRNA microbial diversity characterization using archaeal and bacterial synthetic communities. Environmental microbiology. 2013;15:1882–1899. doi: 10.1111/1462-2920.12086. [PMC free article] [PubMed] [Cross Ref]
53. M. M. OD, et al. The core faecal bacterial microbiome of Irish Thoroughbred racehorses. Letters in applied microbiology. 2013;57:492–501. doi: 10.1111/lam.12137. [PubMed] [Cross Ref]
54. Suchodolski JS, Dowd SE, Wilke V, Steiner JM, Jergens AE. 16S rRNA gene pyrosequencing reveals bacterial dysbiosis in the duodenum of dogs with idiopathic inflammatory bowel disease. PloS one. 2012;7:e39333. doi: 10.1371/journal.pone.0039333. [PMC free article] [PubMed] [Cross Ref]
55. Lagier JC, Million M, Hugon P, Armougom F, Raoult D. Human gut microbiota: repertoire and variations. Frontiers in cellular and infection microbiology. 2012;2:136. doi: 10.3389/fcimb.2012.00136. [PMC free article] [PubMed] [Cross Ref]
56. Sturgeon A, Stull JW, Costa MC, Weese JS. Metagenomic analysis of the canine oral cavity as revealed by high-throughput pyrosequencing of the 16S rRNA gene. Veterinary microbiology. 2013;162:891–898. doi: 10.1016/j.vetmic.2012.11.018. [PubMed] [Cross Ref]
57. Handl S, Dowd SE, Garcia-Mazcorro JF, Steiner JM, Suchodolski JS. Massive parallel 16S rRNA gene pyrosequencing reveals highly diverse fecal bacterial and fungal communities in healthy dogs and cats. FEMS microbiology ecology. 2011;76:301–310. doi: 10.1111/j.1574-6941.2011.01058.x. [PubMed] [Cross Ref]
58. Ginsberg, J. R., Macdonald, D. W., Group, I. S. C. S. & Group, I. S. W. S. Foxes, wolves, jackals, and dogs: an action plan for the conservation of canids. (International Union for Conservation of Nature and Natural Resources, 1990).
59. Wang L, et al. The geographical distribution of grey wolves (Canis lupus) in China: a systematic review. Zoological. Research. 2016;36:315–326. [PMC free article] [PubMed]
60. Lau WN, Fellowes JR, Chan BPL. Carnivores (Mammalia: Carnivora) in South China: a status review with notes on the commercial trade. Mammal Review. 2010;40:247–292. doi: 10.1111/j.1365-2907.2010.00163.x. [Cross Ref]
61. Wang Y, Kasper LH. The role of microbiome in central nervous system disorders. Brain Behavior & Immunity. 2014;38:1. doi: 10.1016/j.bbi.2013.12.015. [PMC free article] [PubMed] [Cross Ref]
62. Fukuda S, Ohno H. Gut microbiome and metabolic diseases. Seminars in Immunopathology. 2014;36:103–114. doi: 10.1007/s00281-013-0399-z. [PubMed] [Cross Ref]
63. Santocchi E, et al. Gut to brain interaction in Autism Spectrum Disorders: a randomized controlled trial on the role of probiotics on clinical, biochemical and neurophysiological parameters. BMC Psychiatry. 2016;16:1–16. doi: 10.1186/s12888-016-0887-5. [PMC free article] [PubMed] [Cross Ref]
64. Matsumoto M, et al. Impact of Intestinal Microbiota on Intestinal Luminal Metabolome. Scientific Reports. 2012;2:233. doi: 10.1038/srep00233. [PMC free article] [PubMed] [Cross Ref]
65. Yatsunenko T, et al. Human gut microbiome viewed across age and geography. Nature. 2012;486:222–227. [PMC free article] [PubMed]
66. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naïve Bayesian Classifier for Rapid Assignment of rRNA Sequences into the New Bacterial Taxonomy. Applied & Environmental Microbiology. 2007;73:5261. doi: 10.1128/AEM.00062-07. [PMC free article] [PubMed] [Cross Ref]
67. Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a Dual-Index Sequencing Strategy and Curation Pipeline for Analyzing Amplicon Sequence Data on the MiSeq Illumina Sequencing Platform. Appl. Environ. Microbiol. 2013;79:5112–5120. doi: 10.1128/AEM.01043-13. [PMC free article] [PubMed] [Cross Ref]
68. Caporaso JG, et al. QIIME allows analysis of high-throughput community sequencing data. Nature methods. 2010;7:335–336. doi: 10.1038/nmeth.f.303. [PMC free article] [PubMed] [Cross Ref]
69. Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nature methods. 2013;10:996. doi: 10.1038/nmeth.2604. [PubMed] [Cross Ref]
70. P Y, et al. The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucleic Acids Research. 2014;42:643–648. doi: 10.1093/nar/gkt888. [PMC free article] [PubMed] [Cross Ref]
71. Team CR, Team RDC. R. A Language And Environment For Statistical Computing. R Foundation for Statistical Computing: Vienna, Austria. Computing. 2011;14:12–21. [PubMed]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group