Skip to main content

Molecular diversity of rumen bacterial communities from tannin-rich and fiber-rich forage fed domestic Sika deer (Cervus nippon) in China



Sika deer (Cervus nippon) have different dietary preferences to other ruminants and are tolerant to tannin-rich plants. Because the rumen bacteria in domestic Sika deer have not been comprehensively studied, it is important to investigate its rumen bacterial population in order to understand its gut health and to improve the productivity of domestic Sika deer.


The rumen bacterial diversity in domestic Sika deer (Cervus nippon) fed oak leaves- (OL group) and corn stalks-based diets (CS group) were elucidated using 16S rRNA gene libraries and denaturing gradient gel electrophoresis (DGGE). Overall, 239 sequences were examined from the two groups, 139 clones from the OL group were assigned to 57 operational taxonomic units (OTUs) and 100 sequences from the CS group were divided into 50 OTUs. Prevotella-like sequences belonging to the phylum Bacteroidetes were the dominant bacteria in both groups (97.2% OL and 77% CS), and sequences related to Prevotella brevis were present in both groups. However, Prevotella shahii-like, Prevotella veroralis-like, Prevotella albensis-like, and Prevotella salivae-like sequences were abundant in the OL group compared to those in the CS group, while Succinivibrio dextrinosolvens-like and Prevotella ruminicola-like sequences were prevalent in the CS group. PCR-DGGE showed that bacterial communities clustered with respect to diets and the genus Prevotella was the dominant bacteria in the rumen of domestic Sika deer. However, the distribution of genus Prevotella from two groups was apparent. In addition, other fibrolytic bacteria, such as Clostridium populeti and Eubacterium cellulosolvens were found in the rumen of domestic Sika deer.


The rumen of domestic Sika deer harbored unique bacteria which may represent novel species. The bacterial composition appeared to be affected by diet, and sequences related to Prevotella spp. may represent new species that may be related to the degradation of fiber biomass or tannins. Moreover, the mechanism and biological functions of Prevotella spp. in the rumen ecosystem, and synergistic interactions with other microorganisms should be noticed.


Sika deer (Cervus nippon) represent the most ancient and primitive members of the genus Cervus because of the simple structure of their antlers, which is very distinct from those of reindeer. Velvet antlers are one of the main products from Sika deer, and are used in traditional Chinese medicine. In addition, Sika deer yield high quality meat and skin. Domestication of Sika deer began much later than for other ruminants. At present, the number of domesticated Sika deer in China is approximately 550,000 head, most of which are distributed in northwestern China.

In nature, Sika deer graze a wide range of forage types, such as Amur grape, elm, maple, bamboo and some toxic species including Chinese Stellera roots and large flowered larkspurs. Moreover, grazing Sika deer have been observed to prefer tannin-rich plants, such as oak leaves. Similar behavior has also been observed in wild Sika deer (Cervus nippon yesoensis) inhabiting the Shiretoko Peninsula of Hokkaido Island in Japan, and in the roe deer (Capreolus capreolus) [1, 2]. However, domesticated Sika deer held in captivity are commonly fed corn stalks containing a much higher fibrous content. Like other ruminants, Sika deer depend on the rumen for fermentation that involves the conversion of plant fiber to volatile fatty acids. This involves a diverse and dense array of microorganisms, including bacteria, fungi, archaea and protozoa [3]. Among these microorganisms bacterial populations have been extensively studied for many years since rumen bacteria have important roles in the efficient degradation of plant biomass and detoxification of secondary compounds in plants [1, 47]. This has led to a variety of studies investigating rumen bacterial structure have been conducted on domestic cows, sheep, yak, Reindeer in Norway and wild Sika deer in Japan [4, 5, 810]. Moreover, rumen bacterial communities are affected by the host and diet [11, 12]. To our knowledge, very little is known about the rumen bacterial community of domesticated Sika deer in China. A comprehensive understanding of bacterial ecology in the rumen of domesticated Sika deer is necessary to increase the efficiency of fiber digestion and to improve the productivity of velvet antlers.

Thus, we hypotheses the bacterial communities in the rumen of domesticated Sika deer may be unique. And the objectives of the present study were: (1) to describe the bacterial diversity in the rumen from domesticated Sika deer ingesting different diets based on 16S rRNA gene sequence libraries and PCR-DGGE; and (2) to compare the unique rumen bacterial populations of domesticated Sika deer ingesting tannin-rich and fiber-rich materials.


Comparative analysis of 16S rRNA gene libraries from two groups

A total of 239 non-chimeric sequences were analyzed, 139 sequences from the OL 16S rRNA clone library and 100 sequences from the CS clone library. The two rumen bacterial populations were distinct according to the RDP classifier tool at a confidence threshold of 80% (Figure 1). Within the two groups, members of the phylum Bacteroidetes were the predominant bacteria (99.3% and 85% of clones in the OL and CS groups, respectively). Domesticated Sika deer consuming corn stalks has Firmicutes present whereas they were not found in oak leaves fed domesticated Sika deer. Similarly, Proteobacteria were more expressed in corn stalks than oak leaves diets. The Chao1 (114.2 vs 143.5) and Shannon-Wiener (3.5 vs 3.7) indices of domesticated Sika deer consuming oak leaves were decreased compared to those feeding on corn stalks (Table 1). Moreover, the Libshuff analysis also showed that the bacterial communities between two diets were significantly differed (P<0.0001). Rarefaction curves at 3% distance levels revealed 74% and 66% coverage for the OL and CS groups, respectively (Figure 2).

Figure 1
figure 1

Composition of 16S rRNA gene libraries at the phylum level. Clones obtained from the OL and CS groups representing by black and grey bars, respectively.

Table 1 Number of OTUs, diversity and coverage at 3% distance level using the MOTHUR platform
Figure 2
figure 2

Rarefaction curves for bacterial 16S rRNA gene libraries. Dark and gray represent Sika deer feeding on oak leaves-based (OL group) and corn stalks-based (CS group) diets, respectively. Rarefaction curves were generated from the platform MOTHUR using the furthest neighbor method.

Using the software program MOTHUR and a sequence identity criterion cut off of 97%, the 139 OL clone sequences were assigned to 57 OTUs and the 100 CS clone sequences were assigned to 50 OTUs (Table 1).To determine the nearest valid related species, the 16S rRNA gene sequences were compared using GenBank’s Basic Local Alignment Search Tool (BLAST). Within the OL library, 53 of the 57 OTUs (i.e. 97.2% of clones) had 85% or greater sequence identities to genus Prevotella (Table 2). Within these OTUs, 23 OTUs (38.1% of clones) showed 87-92% sequence identities to P. brevis, 11 OTUs (16.5% of clones) had 86-90% sequence identities to P. shahii, 3 OTUs (23.8% of clones) had 91-92% sequence identities to P. veroralis, 6 OTUs (12.3% of clones) had distant sequence identities to P. salivae, and the remaining 9 OTUs (6.5% of clones) showed sequence identities to several Prevotella species including P. albensis, P. dentalis, P. ruminicola, P. multiformis, P. stercorea, P. bryantii and P. copri (Table 2). Of the remaining 4 OTUs (of the 57 total OTUs), 2 OTUs (1.4% of clones) were distantly related (85%) to Alistipes shahii, 1 OTU (0.7% of clones) had 84% identity to Barnesiella intestinihominis, and 1 OTU (0.7% of clones) had 97% sequence identity to S. dextrinosolvens.

Table 2 Comparison of 16S rRNA gene libraries between the OL and CS groups

Within the CS clone library, 36 of the 50 OTUs were 85-98% related to species belonging to genus Prevotella. Within these 36 OTUs, only one OTU (2% of clones) had >97% sequence identity to P. brevis, 14 OTUs (36% of clones) had 90-93% identity to P. brevis and 11 OTUs (27% of clones) had 91-95% identity to P. ruminicola making them the dominant bacterial species, whereas the remaining 10 OTUs (12% of clones) exhibited distant sequence identity to P. shahii, P. veroralis, P. albensis, P. salivae and P. dentalis. Of the remaining 14 OTUs (of the 50 total), 3 OTUs (3% of clones) were distantly related (89%) to Paraprevotella clara, 1 OTU (9% of clones) showed 97% identity to S. dextrinosolvens, 3 OTUs (3% of clones) had 90-95% identity to Ruminococcus bromii, 2 OTUs (2% of clones) had 84% identity to Parabacteroides merdae, 1 OTU (1% of clones) was 86% related to Clostridium aldrichii, and 1 OTU (1% of clones) was 91% related to Clostridium bolteae, 4 other OTUs (4% of clones) showed distant sequence identities to Roseburia hominis, Proteiniphilum acetatigenes, A. shahii and Sporanaerobacter acetigenes, respectively.

Overall, phylogenetic analysis revealed that the 107 OTUs were divided into six distinct phylogenetic groups (Figure 3).In addition, the comparison between Norwegian reindeer, Svalbard reindeer and domesticated Sika deer at community level with Fast Unifrac [13], which analyze phylogenetic lineages, showed that the bacterial composition in the rumen of domesticated Sika deer fed oak leaves based diets was more similar to that of domesti-cated Sika deer fed corn stalks based diets, and differed from Svalbard reindeer and Norwegian reindeer (Figure 4). However, there were also shared bacterial communities between domestic Sika deer and Reindeer.

Figure 3
figure 3

Phylogenetic tree of bacterial 16S rRNA sequences from two groups using the Neighbor-Joining method and Kimura two-parameter model in MEGA. Clones from Sika deer fed oak leaves beginning with SDMOL, followed by clone number, and from corn stalks beginning with SDCS, followed by clone number. Aquifex pyrophilus was used as the outgroup. Statistical significance was verified by bootstrapping 1000 replicates.

Figure 4
figure 4

PCoA analysis generating from the UniFrac software coloured by host animals and diets. Norwegian reindeer fed concentrates were represented by dark-blue square; Norwegian reindeer consuming summer pasture were represented by yellow triangle; Svalbard reindeer consuming summer pasture were represented by purple triangle; domesicated Sika deer fed oak leaves were represented by green triangle and domesicated Sika deer fed corn stalks were represented by red round spot.

Rumen bacterial diversity based on the PCR-DGGE profile

PCR-DGGE banding profiles showed that the bacterial communities clustered with respect to diets (Figure 5). However, considerable animal-to-animal variation was also observed. A distinct difference in the bacterial structure was observed between two diets. By comparing the PCR-DGGE profiles between the two diets, the number of DGGE bands from CS group was considerably abundant compared to those from OL group (Figure 5). There were also several bands that were common for all domestic Sika deer.

Figure 5
figure 5

PCR-DGGE profiles of the rumen bacterial 16S rNA gene (V3 region) from domestic Sika deer fed oak leaves (Sika deer A and B) and corn stalks (Sika deer C and D). OL and CS represented Sika deer fed oak leaves and corn stalks, respectively. Three replicates (1, 2 and 3) were taken from each Sika deer. Bionumerics software generated the clustering dendrogram using the UPGMA method.

In total, 47 dominant bands were excised from the PCR-DGGE profile and sequenced, of which 20 and 27 bands obtained from the OL and CS groups, respectively (see Additional file 1). Sequences from the excised bands from the OL group belonged to the phyla Firmicutes, Bacteroidetes and Proteobacteria, whereas DGGE sequences from the CS group belonged to the phyla Firmicutes, Bacteroidetes, Proteobacteria and Synergistetes. Among the 47 bands, 13 bands in two groups were identified as known species based on ≥ 97% sequence similarity (Table 3). Bands O-1, C-3 and C-5 showed ≥ 98% similarity with known species of C. populeti 743A. Bands O-3 and O-18 were identified as Streptococcus pasteurianus CIP 107122, while bands O-9 and C-14 showed 98% similarity with of Eubacterium cellulosolvens 6. Band O-12 displayed 97% similarity with known species of Moryella indoligenes AIP 220.04, and band O-13 showed species-level sequence similarity to Pseudobutyrivibrio ruminis DSM9787. Bands O-10 and C-10 displayed 98% similarity to Succinivibrio dextrinosolvens 0554, while bands C-18 and C-1 had 98% sequence similarity to Coprococcus eutactus ATCC 27759 and Prevotella ruminicola ATCC 19189, respectively. Moreover, band C-21 had the 93% similarity with known species of Eubacterium ruminan-tium GA 195. Bands C-13 and C-22 were distantly related to Galbibacter mesophilus Mok-17 with 88% and 91% similarity, respectively. Band C-24 displayed 88% similarity with Capnocytophaga cynodegmi CIP 103937, and band C-27 showed 94% similarity with known species of Bacteroides uniformis JCM 5828. Bands C-19 and C-20 had 92% similarity with known species of Dethiosulfovibrio acidaminovorans sr15. The remaining 30 bands from two groups had 92-96% sequence similarities with several species belonging to genus Prevotella including P. loescheii, P. pleuritidis, P. corporis, P. buccalis, P. dentalis, P. melani-nogenica, P. salivae, P. copri, P. denticola, P. oulorum and P. histicola.

Table 3 Sequences analysis of V3 region of 16S rRNA gene from PCR-DGGE


In the present study, two 16S rRNA gene libraries and PCR-DGGE were used to study the rumen bacteria in the rumen of domesticated Sika deer feeding on oak leaves-based (OL) and corn stalks-based (CS) diets. Sequences from the two clone libraries and PCR-DGGE bands indicated that the majority of sequences belonged to phylum Bacteroidetes. The findings from the current study are similar to previous findings for other ruminants, such as Reindeer, yaks, cattle and goats [1418]. The predominance of sequences belonging to the phylum Bacteroidetes highlights their important role in the rumen fermentation of domesticated Sika deer. While, the phylum Firmicutes being prevalent in other ruminants were not found in the OL library [19], which could be caused by the tannins contained in oak leaves, because some studies reported that the tannins in oak leaves may have a negative effect on some of the bacterial species [20, 21], and the growth of proteolytic bacteria, such as Butyrivibrio fibrisolvens, Ruminococcus albus and Streptococcus bovis, were inhibited by tannins [22, 23]. This may also indicate that some species belonging to phylum Firmicutes in the rumen of domestic Sika deer may be sensitive to tannins.

Within the phylum Bacteroidetes, Prevotella-like clones accounted for 97.2% of the clones in the OL group and 77% in the CS group. Moreover, the PCR-DGGE results also showed the genus Prevotella represented the predominant bacteria in rumen of domesticated Sika deer (Table 3), which is in agreement with other studies [19, 2428] . The prevalence of Prevotella spp. in rumen fermentation of domesticated Sika deer was likely because they utilize a wide variety of polysaccharides, and are thought to be important contributors to xylan degradation in the rumen [2932]. Although other studies found that concentrate diets increased the numbers of clones related to Prevotella spp. [33, 34], however, in comparison with other ruminants, there was an apparent difference in the proportion of Prevotella spp. [6, 25, 27, 28]. Prevotella spp. belonged to the hydrogen-consuming bacteria, which could produce propionate via succinate or acrylate pathways though fermentation of sugars and lactate, respectively [3537]. Therefore, the dominant genus Prevotella in the rumen of domesticated Sika deer suggested that the propionate pathway may be relatively vital in the rumen fermentation of domestic Sika deer, which, in turn, may lead to the decreased production of methane, since the succinate-propionate pathway could compete with methanogens for hydrogen [38]. The relationship between Prevotella spp. and methanogens in the rumen of domesticated Sika deer was worth of further investigating. In addition, the bacterial communities in the rumen between domesticated Sika deer, Svalbard reindeer and Norwegian reindeer, all cervids, were compared using Fast UniFrac, which can be used to determine whether communities are significantly different [13]. The results of Principal coordinate analysis (PCoA) between domesticated Sika deer and Reindeer using the Fast Unifrac platform clearly showed that the rumen bacterial communities were distinct, which can be attributed to the host-species (Figure 5) [13, 26, 39].

It is important to note, that fibrolytic bacteria, such as C. populeti, E. cellulosolvens and Ps. ruminis were discovered in our analysis based on PCR-DGGE, rather than the predominant fibrolytic bacteria, B. fibrisolvens, Fibrobacter succinogenes, Ruminococcus flavefaciens and R. albus. This may suggest that the rumen of domesticated Sika deer depend on unique bacterial communities in rumen fermentation. In contrast, the absence of R. flavefaciens, B. fibrisolvens, F. succinogenes and R. albus in the present work may be attributed to the small number of clones may have missed some other members of the bacterial community, and the weak or unidentifiable bands in DGGE. Future work will employ next generation sequencing to effectively elucidate the bacterial diversity present in the rumen of domesticated Sika deer and other livestock. Collectively, these data indicated that the rumen of domesticated Sika deer harbored unique bacterial populations for the fermentation of plant biomass and concentrate diet.

Interestingly, in both clone libraries, none of the sequences were 100% identical. Rather, most clones were in the range of 83-98% identify to known species in both libraries. These results suggested that the rumen bacteria of domesticated Sika deer were not previously characterized and that these clones related to Prevotella spp. in the rumen represented new species. This agrees with previous findings suggesting that most of the bacterial species in rumen of other cervids (96% for Hokkaido Sika deer and 100% for Svalbard reindeer) are unknown [26, 40]. Despite the diets and geographic location are important factors affecting bacterial diversity in the rumen, however, the presence of these unknown or unidentified species may be the result of co-evolution between microbial communities and the host.

PCR-DGGE analysis showed that the bacterial diversity in domesticated Sika deer fed corn stalks differed from the domesticated Sika deer consuming oak leaves (Figure 5), indicating forage affected the relative abundance and composition of the bacteria. Moreover, the difference in the Prevotella species between the two groups was very apparent (Table 3). For instance, the results of clone library showed that the proportion of P. ruminicola-like clones (27%) was abundant in the CS group comparing with those in the OL group, and sequences analysis of PCR-DGGE also indicated that P. ruminicola was only presented in CS group. Interestingly, Prevotella species in the rumen could contribute to cell wall degradation through synergistic interactions with species of cellulolytic bacteria [41]. Therefore, considering the relatively high fiber content (about 36%) in corn stalks, these P. ruminicola-like clones in the CS group may play a role in the degradation of cellulose. This explanation is partly supported by recent metagenomics data from the Svalbard reindeer rumen microbiome, where the presence of polysaccharide utilizing glycoside hydrolase and other carbohydrate-active enzyme families target various polysaccharides including cellulose, xylan and pectin [18].

In the OL group, the distribution of P. shahii-like clones (16.5%), P. veroralis-like clones (23.8%) and P. salivae-like clones (12.3%) were several times higher in the OL library than in the CS library, and several bands in the PCR-DGGE analysis showed sequence similarities to P. salivae (Table 3). Previous study reported that P. ruminicola may tolerate condensed tannins [22]. Considering the genetic diversity of Prevotella spp. [27, 42], it is assumed that the tolerance to tannins of domestic Sika deer may be related to the abundance of Prevotella spp. in the OL group. In addition, we found two bands (O-3 and O-18) were identified as St. pasteurianus using PCR-DGGE. Thus this species may also be important in the process of degrading tannins in diets, because tannin-degrading capability of Streptococcus sp. have been demonstrated in other studies [4346]. However, these assumptions need to be investigated in future studies.

Phylogenetic analysis indicated the presence of diet-specific subpopulations of Prevotella. Prevotella clusters 1 and 2 not only demonstrated the genetic diversity of Prevotella spp., but also confirmed the above assumption that clones grouped within clusters 1 or 2 may be related to the degradation of fiber (cluster 1) or tannins (cluster 2), whereas, the clones in cluster 3 may have common features of degrading starch and proteins contained in concentrate diets (Figure 3). However, clones related to the bacterial genera Sporanaerobacter, Parabacteroides and Proteiniphilum were found in the rumen of domesticated Sika deer fed corn stalks that were not previously reported in the rumen from other ruminants. Sporanaerobacter acetigenes is an acetogenic and a sulfur-reducing bacterium that was isolated from an anaerobic sludge blanket reactor in Mexico [47, 48]. The rumen has considerable capacity to convert sulfate into sulfur-containing amino acids. Similarly, little is known about Proteiniphilum acetatigenes, which was originally isolated from a UASB reactor treating brewery wastewater in China [49]. These bacteria in rumen of domesticated Sika deer may have other biological functions and is worthy of further investigation.


In conclusion, this study is the first to report the rumen bacteria in Chinese domesticated Sika deer, consuming either oak leaves-based or corn stalks-based diets. Sequences analysis from 16S rRNA clone libraries and PCR-DGGE revealed that the domesticated Sika deer harbored unique rumen bacterial populations, most of which may present novel species, and that the bacterial compositions were affected by forage. It is speculated that the possible new species of Prevotella may be related to the degradation of tannins or fiber biomass. Moreover, the species diversity of Prevotella sp. in the rumen combined with their synergistic interactions with other microorganisms requires further in depth investigation.


Animals and sampling

Four male rumen-cannulated domestic Sika deer (Cervus nippon) maintained at the research farm (44.04° N, 129.09° E) of the Institute of Special Animal and Plant Sciences, Chinese Academy of Agricultural Sciences, in Jilin Province, were used in this study. From September to October, four domestic Sika deer were offered the same concentrated diets (64.5% corn, 19.7% soybean meal, 12.8% distiller dried grains with solubles and a 3% mixture of vitamins and mineral salts) and mixed with either oak leaves (OL) or corn stalks (CS). All domestic Sika deer were fed twice each day at 8:00 AM and 4:00 PM and had free access to water. The whole rumen contents, which included solid and fluid fractions, were collected via rumen cannula before the morning feeding, and stored at -20°C for analysis. All domestic Sika deer used in present experiment must be performed according to the animal health and well-being regulations, all animal procedures were approved and authorized by the Chinese Academy of Agricultural Sciences Animal Care and Use Committee, and by the Wild Animal and Plant Subcommittee, Institute of Special Animal and Plant Sciences.

DNA extraction

Total DNA was directly extracted from rumen contents containing solid and liquid fraction according to methods described by LaMontagne [50] with few modifications. In brief, 800 μl lysis buffer (0.15 M NaCl, 0.2 M EDTA, 10 lysozyme, pH8.0), 20 μl of 20 proteinase K (Sigma, Germany), and 0.3 g glass beads (0.1 mm, Sigma, Germany) were added to 0.5 g of whole rumen contents. After shaking at 37°C for 1 h, 300 μl heated lysis buffer (10% SDS, 0.1 M NaCl, 0.5 M Tris–HCl, pH8.0) at 65°C, 300 μl phosphate buffer (pH8.0) and 600 μl chloroform-isoamyl alcohol (24:1, V/V) were added, and the mixture was incubated at 65°C in a water bath for 30 min with intense shaking 30 s at 10 min intervals. After centrifugation at 5,000 rpm for 6 min, the supernatant was transferred to a clean tube. DNA was then precipitated with a 0.6 volume of isopropanol at -80°C for 15 min, and the pellet was washed several times with 75% ethanol. The DNA was dried and dissolved in TE buffer (pH 8.0). The DNA quality was assessed by 0.8% agarose gel electrophoresis, and the purity was determined by spectrophotometry (SPECORD 50, analytikjena, Germany), after which it was purified using a QIAEX II Gel Extraction Kit (QIAGEN, Germany).

Construction of 16S rRNA gene clone libraries and sequences analyses

Universal primers 27F (5-AGAGTTTGATCMTGGCTCAG-3) and 1492R (5-TACGGYTACCTTGTTACGACTT-3) were used to amplify the 16S rRNA gene (approximately 1.5 kb) [51]. Each 50 ul reaction contained 50 ng template DNA, 0.25 mM of each primer, 250 mM dNTPs, 1.25 U of Ex Taq and 5 μl Ex Taq buffer (TaKaRa, Dalian). PCR was performed on a 2720 Thermal Cycler (Applied Biosystems, USA) with hot start at 94°C for 5 min, followed by 20 cycles of 30 s at 94°C, 1 min at 55°C and 2 min at 72°C; and a final extension at 72°C for 10 min. The PCR product was assessed using 2% agarose gel electrophoresis (approximately 1.5 kb), and were purified using a TaKaRa MiniBEST DNA Fragment Purification Kit (TaKaRa, Dalian) and then pooled within each group. Two 16S rRNA gene clone libraries were constructed from the pooled PCR products using the TOPO® TA Cloning® Kit (Invitrogen, USA). Positive (white) clones were screened by colony PCR with the M13 Forward and M13 Reverse primers, and sequenced using an ABI 3730XL DNA Analyzer.

The chimera check program Bellerophon was used to identify chimeric sequences [52]. The remaining sequences were assigned using the Classifier tool available at Ribosomal Database Project (RDP) Release 10 at a confidence threshold of 80% [53]. OTUs based on 97% sequence identity, and the Shannon-Wiener index-based diversity estimator and the Chao1 based index of richness were calculated using MOTHUR platform to determine the diversity and richness of bacterial communities in each group based on the 16S rRNA gene libraries [54]. Libshuff analysis was performed to estimate the similarity between libraries from two diets based on evolutionary distance of all sequences. Coverage and rarefaction curves were also determined using the MOTHUR platform [54]. The 16S rRNA gene sequences were screened using GenBank’s BLAST program [55]. The closest related sequences were retrieved and aligned with sequences from the present study using the CLUSTALW 1.83 program in MEGA 5.05 software [56]. A phylogenetic tree was constructed using the Kimura two-parameter model and the Neighbor-Joining method as part of the MEGA 5.05 software. The statistical significance was verified by 1000 bootstrapped replicates. The sequences obtained from this study were submitted to GenBank under the accession numbers JX889268 to JX889378. Furthermore, an unweighted UniFrac distance matrix was constructed from the phylogenetic tree of clone libraries of Norwegian reindeer, Svalbard reindeer and Sika deer, and was visualized using PCoA [13, 26, 39].

PCR-DGGE banding profiles and statistical analysis

The variable region (V3) of the bacterial 16S rRNA gene was amplified using the primers of F341GC and R534, and PCR condition was described previously [57]. A 40 bp GC-clamp (5′-CGCCCGGGGCGCGCCCCGGGCGGGGCGGGGGCACGGGGGG-3′) was on the 5′ end of the F341 primer. The PCR products were loaded onto 8% polyacrylamide gels (37.5:1) with a denaturing gradient of 40–60% at 80V over 16 h at 60°C. Electrophoresis was performed using Bio-Rad’s DCode detection system. The gels were stained with SYBR Green I (Invitrogen, USA) for 25 min and gel images were captured using the Gel Doc™ XR+ system (BIO-RAD, CA). Cluster analysis was performed using a Dice similarity coefficient at 0.5% optimization and 1% tolerance following the unweighted pair-group method using arithmetic averages (UPGMA) on BioNumerics 6.0 software (Applied-Maths, Kortrijk, Belgium).

Dominant bands were excised from DGGE gel and eluted overnight in 500 μl of sterilized ddH2O at 4°C. Extracted DNA was re-amplified using PCR primers F341 and R534 without GC-clamp. The size of PCR products were determined using agarose gel and were purified using QIAquick® PCR Purification Kit (Qiagen, USA). The PCR products were cloned into TOPO® TA Cloning® Kit with TOP 10 according to the manufacturer’s instruction (Invitrogen, San Diego, CA, USA). Recombinant plasmids of positive clones (white) were sequenced using ABI 3730XL DNA Analyzer. The sequences were compared with those sequences deposited in NCBI web site using BLAST program [55].



Denaturing gradient gel electrophoresis


Operational taxonomic units


Volatile fatty acids


Basic local alignment search tool


Ribosomal database project


Principal coordinates analysis.


  1. Hiura T, Hashidoko Y, Kobayashi Y, Tahara S: Effective degradation of tannic acid by immobilized rumen microbes of a sika deer (Cervus nippon yesoensis) in winter. Anim Feed Sci Technol. 2010, 155 (1): 1-8. 10.1016/j.anifeedsci.2009.09.015.

    Article  CAS  Google Scholar 

  2. Clauss M, Lason K, Gehrke J, Lechner-Doll M, Fickel J, Grune T, Jurgen Streich W: Captive roe deer (Capreolus capreolus) select for low amounts of tannic acid but not quebracho: fluctuation of preferences and potential benefits. Comp Biochem Physiol B Biochem Mol Biol. 2003, 136 (2): 369-382. 10.1016/S1096-4959(03)00244-6.

    Article  PubMed  Google Scholar 

  3. Wright A-DG, Klieve AK: Does the complexity of the rumen microbial ecology preclude methane mitigation?. Animal Feed Sci. Technol. 2011, 166-167: 248-253.

    Article  CAS  Google Scholar 

  4. Tajima K, Arai S, Ogata K, Nagamine T, Matsui H, Nakamura M, Aminov RI, Benno Y: Rumen bacterial community transition during adaptation to high-grain diet. Anaerobe. 2000, 6 (5): 273-284. 10.1006/anae.2000.0353.

    Article  CAS  Google Scholar 

  5. An DD, Dong XZ, Dong ZY: Prokaryote diversity in the rumen of yak (Bos grunniens) and Jinnan cattle (Bos taurus) estimated by 16S rDNA homology analyses. Anaerobe. 2005, 11 (4): 207-215. 10.1016/j.anaerobe.2005.02.001.

    Article  PubMed  CAS  Google Scholar 

  6. Pei CX, Liu QA, Dong CS, Li HQ, Jiang JB, Gao WJ: Diversity and abundance of the bacterial 16S rRNA gene sequences in forestomach of alpacas (Lama pacos) and sheep (Ovis aries). Anaerobe. 2010, 16 (4): 426-432. 10.1016/j.anaerobe.2010.06.004.

    Article  PubMed  CAS  Google Scholar 

  7. Yang LY, Chen J, Cheng XL, Xi DM, Yang SL, Deng WD, Mao HM: Phylogenetic analysis of 16S rRNA gene sequences reveals rumen bacterial diversity in Yaks (Bos grunniens). Mol Biol Rep. 2010, 37 (1): 553-562. 10.1007/s11033-009-9794-x.

    Article  PubMed  CAS  Google Scholar 

  8. Aagnes TH, Sormo W, Mathiesen SD: Ruminal microbial digestion in free-living, in captive lichen-ded, and in Starved Reindeer (Rangifer tarandus tarandus) in winter. Appl Environ Microbiol. 1995, 61 (2): 583-591.

    PubMed  CAS  PubMed Central  Google Scholar 

  9. Edwards JE, McEwan NR, Travis AJ, Wallace RJ: 16S rDNA library-based analysis of ruminal bacterial diversity. Antonie Leeuwenhoek Int J Gen Mol Microbiol. 2004, 86 (3): 263-281.

    Article  CAS  Google Scholar 

  10. Ichimura Y, Yamano H, Takano T, Koike S, Kobayashi Y, Tanaka K, Ozaki N, Suzuki M, Okada H, Yamanaka M: Rumen microbes and fermentation of wild sika deer on the Shiretoko peninsula of Hokkaido Island, Japan. Ecol Res. 2004, 19 (4): 389-395. 10.1111/j.1440-1703.2004.00649.x.

    Article  Google Scholar 

  11. Kocherginskaya SA, Aminov RI, White BA: Analysis of the rumen bacterial diversity under two different diet conditions using denaturing gradient gel electrophoresis, random sequencing, and statistical ecology approaches. Anaerobe. 2001, 7 (3): 119-134. 10.1006/anae.2001.0378.

    Article  CAS  Google Scholar 

  12. Shi PJ, Meng K, Zhou ZG, Wang YR, Diao QY, Yao B: The host species affects the microbial community in the goat rumen. Lett Appl Microbiol. 2008, 46 (1): 132-135.

    PubMed  CAS  Google Scholar 

  13. Lozupone C, Knight R: UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005, 71 (12): 8228-8235. 10.1128/AEM.71.12.8228-8235.2005.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  14. Cho SJ, Cho KM, Shin EC, Lim WJ, Hong SY, Choi BR, Kang JM, Lee SM, Kim YH, Kim H, et al.: 16S rDNA analysis of bacterial diversity in three fractions of cow rumen. J Microbiol Biotechnol. 2006, 16 (1): 92-101.

    CAS  Google Scholar 

  15. Yang SL, Ma SC, Chen J, Mao HM, He YD, Xi DM, Yang LY, He TB, Deng WD: Bacterial diversity in the rumen of Gayals (Bos frontalis), Swamp buffaloes (Bubalus bubalis) and Holstein cow as revealed by cloned 16S rRNA gene sequences. Mol Biol Rep. 2010, 37 (4): 2063-2073. 10.1007/s11033-009-9664-6.

    Article  PubMed  CAS  Google Scholar 

  16. Cunha IS, Barreto CC, Costa OYA, Bomfim MA, Castro AP, Kruger RH, Quirino BF: Bacteria and archaea community structure in the rumen microbiome of goats (Capra hircus) from the semiarid region of Brazil. Anaerobe. 2011, 17 (3): 118-124. 10.1016/j.anaerobe.2011.04.018.

    Article  PubMed  Google Scholar 

  17. Li MJ, Zhou M, Adamowicz E, Basarab JA, Guan LL: Characterization of bovine ruminal epithelial bacterial communities using 16S rRNA sequencing, PCR-DGGE, and qRT-PCR analysis. Vet Microbiol. 2012, 155 (1): 72-80. 10.1016/j.vetmic.2011.08.007.

    Article  PubMed  CAS  Google Scholar 

  18. Pope PB, Mackenzie AK, Gregor I, Smith W, Sundset MA, McHardy AC, Morrison M, Eijsink VG: Metagenomics of the Svalbard reindeer rumen microbiome reveals abundance of polysaccharide utilization loci. PLoS One. 2012, 7 (6): e38571-10.1371/journal.pone.0038571.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Kim M, Morrison M, Yu Z: Status of the phylogenetic diversity census of ruminal microbiomes. FEMS Microbiol Ecol. 2011, 76 (1): 49-63. 10.1111/j.1574-6941.2010.01029.x.

    Article  PubMed  CAS  Google Scholar 

  20. Bae HD, McAllister TA, Yanke J, Cheng KJ, Muir AD: Effects of condensed tannins on endoglucanase activity and filter paper digestion by Fibrobacter succinogenes S85. Appl Environ Microbiol. 1993, 59 (7): 2132-2138.

    PubMed  CAS  PubMed Central  Google Scholar 

  21. McSweeney CS, Palmer B, McNeill DM, Krause DO: Microbial interactions with tannins: nutritional consequences for ruminants. Anim Feed Sci Technol. 2001, 91 (1–2): 83-93.

    Article  CAS  Google Scholar 

  22. Jones GA, McAllister TA, Muir AD, Cheng KJ: Effects of sainfoin (Onobrychis viciifolia Scop.) condensed tannins on growth and proteolysis by four strains of ruminal bacteria. Appl Environ Microbiol. 1994, 60 (4): 1374-1378.

    PubMed  CAS  PubMed Central  Google Scholar 

  23. Min BR, Attwood GT, McNabb WC, Molan AL, Barry TN: The effect of condensed tannins from Lotus corniculatus on the proteolytic activities and growth of rumen bacteria. Anim Feed Sci Technol. 2005, 121 (1–2): 45-58.

    Article  CAS  Google Scholar 

  24. 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. 10.1016/S0378-1097(03)00760-2.

    Article  PubMed  CAS  Google Scholar 

  25. Stevenson DM, Weimer PJ: Dominance of Prevotella and low abundance of classical ruminal bacterial species in the bovine rumen revealed by relative quantification real-time PCR. Appl Microbiol Biotechnol. 2007, 75 (1): 165-174. 10.1007/s00253-006-0802-y.

    Article  PubMed  CAS  Google Scholar 

  26. 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-438. 10.1007/s00248-007-9254-x.

    Article  PubMed  Google Scholar 

  27. Bekele AZ, Koike S, Kobayashi Y: Genetic diversity and diet specificity of ruminal Prevotella revealed by 16S rRNA gene-based analysis. FEMS Microbiol Lett. 2010, 305 (1): 49-57. 10.1111/j.1574-6968.2010.01911.x.

    Article  PubMed  CAS  Google Scholar 

  28. Wu S, Baldwin RL, Li W, Li C, Connor EE, Li RW: The bacterial community composition of the bovine rumen detected using pyrosequencing of 16S rRNA genes. Metagenomics. 2012, 1: 1-11.

    Article  Google Scholar 

  29. Cotta MA: Interaction of ruminal bacteria in the production and utilization of maltooligosaccharides from starch. Appl Environ Microbiol. 1992, 58 (1): 48-54.

    PubMed  CAS  PubMed Central  Google Scholar 

  30. Gardner RG, Wells JE, Russell JB, Wilson DB: The cellular location of Prevotella ruminicola beta-1,4-D-endoglucanase and its occurrence in other strains of ruminal bacteria. Appl Environ Microbiol. 1995, 61 (9): 3288-3292.

    PubMed  CAS  PubMed Central  Google Scholar 

  31. Matsui H, Ogata K, Tajima K, Nakamura M, Nagamine T, Aminov RI, Benno Y: Phenotypic characterization of polysaccharidases produced by four Prevotella type strains. Curr Microbiol. 2000, 41 (1): 45-49. 10.1007/s002840010089.

    Article  PubMed  CAS  Google Scholar 

  32. Krause DO, Denman SE, Mackie RI, Morrison M, Rae AL, Attwood GT, McSweeney CS: Opportunities to improve fiber degradation in the rumen: microbiology, ecology, and genomics. FEMS Microbiol Rev. 2003, 27 (5): 663-693. 10.1016/S0168-6445(03)00072-X.

    Article  PubMed  CAS  Google Scholar 

  33. Fernando SC, Purvis HT, Najar FZ, Sukharnikov LO, Krehbiel CR, Nagaraja TG, Roe BA, DeSilva U: Rumen microbial population dynamics during adaptation to a high-grain diet. Appl Environ Microbiol. 2010, 76 (22): 7482-7490. 10.1128/AEM.00388-10.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Sadet-Bourgeteau S, Martin C, Morgavi DP: Bacterial diversity dynamics in rumen epithelium of wethers fed forage and mixed concentrate forage diets. Vet Microbiol. 2010, 146 (1–2): 98-104.

    Article  PubMed  CAS  Google Scholar 

  35. Strobel HJ: Vitamin B12-dependent propionate production by the ruminal bacterium Prevotella ruminicola 23. Appl Environ Microbiol. 1992, 58 (7): 2331-2333.

    PubMed  CAS  PubMed Central  Google Scholar 

  36. Purushe J, Fouts DE, Morrison M, White BA, Mackie RI, North American Consortium for Rumen B, Coutinho PM, Henrissat B, Nelson KE: Comparative genome analysis of Prevotella ruminicola and Prevotella bryantii: insights into their environmental niche. Microb Ecol. 2010, 60 (4): 721-729. 10.1007/s00248-010-9692-8.

    Article  PubMed  Google Scholar 

  37. Flint HJ, Bayer EA, Rincon MT, Lamed R, White BA: Polysaccharide utilization by gut bacteria: potential for new insights from genomic analysis. Nat Rev Microbiol. 2008, 6 (2): 121-131. 10.1038/nrmicro1817.

    Article  PubMed  CAS  Google Scholar 

  38. Newbold CJ, Lopez S, Nelson N, Ouda JO, Wallace RJ, Moss AR: Propionate precursors and other metabolic intermediates as possible alternative electron acceptors to methanogenesis in ruminal fermentation in vitro. Br J Nutr. 2005, 94 (1): 27-35. 10.1079/BJN20051445.

    Article  PubMed  CAS  Google Scholar 

  39. Sundset M, Edwards J, Cheng Y, Senosiain R, Fraile M, Northwood K, Praesteng KE, Glad T, Mathiesen S, Wright AD: Molecular diversity of the rumen microbiome of Norwegian Reindeer on natural summer pasture. Microb Ecol. 2009, 57 (2): 335-348. 10.1007/s00248-008-9414-7.

    Article  PubMed  CAS  Google Scholar 

  40. Kobayashi Y: Inclusion of novel bacteria in rumen microbiology: need for basic and applied science. Anim Sci J. 2006, 77 (4): 375-385. 10.1111/j.1740-0929.2006.00362.x.

    Article  CAS  Google Scholar 

  41. Whitehead TR: Analyses of the gene and amino acid sequence of the Prevotella (Bacteroides) ruminicola 23 xylanase reveals unexpected homology with endoglucanases from other genera of bacteria. Curr Microbiol. 1993, 27 (1): 27-33. 10.1007/BF01576830.

    Article  PubMed  CAS  Google Scholar 

  42. Ramsak A, Peterka M, Tajima K, Martin JC, Wood J, Johnston MEA, Aminov RI, Flint HJ, Avgustin G: Unravelling the genetic diversity of ruminal bacteria belonging to the CFB phylum. FEMS Microbiol Ecol. 2000, 33 (1): 69-79. 10.1111/j.1574-6941.2000.tb00728.x.

    Article  PubMed  CAS  Google Scholar 

  43. Brooker JD, O’Donovan LA, Skene I, Clarke K, Blackall L, Muslera P: Streptococcus caprinus sp. nov, a tannin-resistant ruminal bacterium from feral goats. Lett Appl Microbiol. 1994, 18 (6): 313-318. 10.1111/j.1472-765X.1994.tb00877.x.

    Article  CAS  Google Scholar 

  44. Sly LI, Cahill MM, Osawa R, Fujisawa T: The tannin-degrading species Streptococcus gallolyticus and Streptococcus caprinus are subjective synonyms. Int J Syst Bacteriol. 1997, 47 (3): 893-894. 10.1099/00207713-47-3-893.

    Article  PubMed  CAS  Google Scholar 

  45. Chamkha M, Patel BK, Traore A, Garcia JL, Labat M: Isolation from a shea cake digester of a tannin-degrading Streptococcus gallolyticus strain that decarboxylates protocatechuic and hydroxycinnamic acids, and emendation of the species. Int J Syst Evol Microbiol. 2002, 52 (Pt 3): 939-944.

    PubMed  CAS  Google Scholar 

  46. Goel G, Puniya AK, Singh K: Tannic acid resistance in ruminal streptococcal isolates. J Basic Microbiol. 2005, 45 (3): 243-245. 10.1002/jobm.200410517.

    Article  PubMed  CAS  Google Scholar 

  47. Eaton HL, De Lorme M, Chaney RL, Craig AM: Ovine ruminal microbes are capable of biotransforming hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX). Microb Ecol. 2011, 62 (2): 274-286. 10.1007/s00248-011-9809-8.

    Article  PubMed  CAS  Google Scholar 

  48. Hernandez-Eugenio G, Fardeau ML, Cayol JL, Patel BK, Thomas P, Macarie H, Garcia JL, Ollivierx P: Sporanaerobacter acetigenes gen. nov., sp. nov., a novel acetogenic, facultatively sulfur-reducing bacterium. Int J Syst Evol Microbiol. 2002, 52 (Pt 4): 1217-1223.

    PubMed  CAS  Google Scholar 

  49. Chen S, Dong X: Proteiniphilum acetatigenes gen. nov., sp. nov., from a UASB reactor treating brewery wastewater. Int J Syst Evol Microbiol. 2005, 55 (Pt 6): 2257-2261.

    Article  PubMed  CAS  Google Scholar 

  50. LaMontagne MG, Michel FC, Holden PA, Reddy CA: Evaluation of extraction and purification methods for obtaining PCR-amplifiable DNA from compost for microbial community analysis. J Microbiol Methods. 2002, 49 (3): 255-264. 10.1016/S0167-7012(01)00377-3.

    Article  PubMed  CAS  Google Scholar 

  51. Lane DJ: 16S/23S rRNA sequencing. Nucleic acid techniques in bacteria systematics. Edited by: Stackebrandt EGM. 1991, New York: Wiley, 115-175.

    Google Scholar 

  52. Huber T, Faulkner G, Hugenholtz P: Bellerophon: a program to detect chimeric sequences in multiple sequence alignments. Bioinformatics. 2004, 20 (14): 2317-2319. 10.1093/bioinformatics/bth226.

    Article  PubMed  CAS  Google Scholar 

  53. Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ, Kulam-Syed-Mohideen AS, McGarrell DM, Marsh T, Garrity GM, et al.: The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 2009, 37: D141-D145. 10.1093/nar/gkn879.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  54. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al.: Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009, 75 (23): 7537-7541. 10.1128/AEM.01541-09.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  55. Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  56. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  57. Muyzer G, de Waal EC, Uitterlinden AG: Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl Environ Microbiol. 1993, 59 (3): 695-700.

    PubMed  CAS  PubMed Central  Google Scholar 

Download references


Special thanks to Dr. Yanfeng Cheng in the analysis of 16S rRNA gene sequences and Dr. Zhi Gang Cui in the analysis of PCR-DGGE. This work is supported by a Young Researcher funded project (201101086), Science and Technology Development project (20090238) and a Leading Talent and Creative Team project (20121810), all from Jilin province, the Ministry of Agriculture Public Sector (Agriculture) Special Research Project (200903014) and Key Projects in the National Science & Technology Pillar Program (2011BAI03B02).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Guang Yu Li.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ZPL sampled rumen contents, extracted DNA, constructed the clone library, data analysis and drafted the manuscript. ADGW was involved with interpretation of data and with preparing the manuscript. HLL designed the study and drafted the paper. KB, YFY, CX and KYW contributed to sample rumen contents and all of lab works. GYL and FHYconceived the study. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Dominant bands of PCR-DGGE banding patterns of bacteria 16SrRNA gene (V3 region). In the text, bands from OL group were defined as O and followed by bands number, bands from CS group begin with C and followed by bands numbers. (PDF 86 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Li, Z.P., Liu, H.L., Li, G.Y. et al. Molecular diversity of rumen bacterial communities from tannin-rich and fiber-rich forage fed domestic Sika deer (Cervus nippon) in China. BMC Microbiol 13, 151 (2013).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: