- Open Access
Effects of dietary replacement of soybean meal with dried distillers grains with solubles on the microbiota occupying different ecological niches in the rumen of growing Hu lambs
Journal of Animal Science and Biotechnology volume 11, Article number: 93 (2020)
Diet has a profound impact on the rumen microbiota, and the impact can vary among the different rumen ecological niches (REN). This study investigated the effects of dietary replacement of soybean meal (SBM) with dried distillers grains with solubles (DDGS) on the rumen microbiota occupying different REN of growing Hu lambs. After a 9-week feeding trial, 6 lambs from each dietary treatment (SBM vs. DDGS-based diets) were slaughtered for sample collection. The microbiota of the rumen solid, liquid, and epithelium fractions was examined using amplicon sequencing analysis of bacterial 16S rRNA gene, functional prediction, and qPCR.
No interaction of dietary protein source (PS) and REN were detected for virtually all the measurements made in this study. The DDGS substitution resulted in very limited influence on bacterial community structure. However, the metabolic pathways predicted from 16S rRNA gene sequences varied greatly between SBM- and DDGS-based diets. The populations of rumen total bacteria, fungi, sulfate-reducing bacteria (SRB), and methanogens were not influenced by DDGS substitution, but the population of protozoa was reduced. The bacterial communities in rumen solid (RS) and liquid (RL) were similar in taxonomic composition but were different in relative abundance of some taxa. In contrast, the bacterial composition and relative abundance of rumen epithelium (RE) were greatly distinct from those of the RS and the RL. In alignment with the bacterial relative abundance, the metabolic pathways predicted from 16S rRNA genes also varied greatly among the different REN. The populations of total bacteria, protozoa, and methanogens attached to the RE were smaller than those in the RS and RL, and the fungal population on the rumen epithelium was smaller than that in the RS but similar to that in the RL. On the contrary, the SRB population on the RE was greater than that in the RS and RL.
Substitution of SBM with DDGS had greater impact to the protozoa than to the other microbes, and the microbial community structure and functions at different REN are distinct and niche-adapted.
Rumen microbiota determines animal health and production performance, and a better understanding of this complex and diverse microbiota, especially the relationship between the structure and function under different dietary conditions can lead to innovative interventions to improve animal productivity . Dietary changes can shift rumen microbiota with respect to its composition and structure and fermentation patterns [2, 3]. Dried distillers grains with solubles (DDGS) is an abundant byproduct of ethanol production from grains (primarily corn), which are high in crude protein (CP), ether extract (EE) and energy content [4, 5], and low in cost . It has been widely used as a substitute for corn and soybean meal (SBM) in ruminant production [7, 8]. Several studies have evaluated the impact of feeding DDGS on rumen bacterial community in steers [9,10,11] and dairy cattle [12, 13]. However, the influences of DDGS addition on rumen microbial community structure and function in growing lambs remain poorly understood. In our previous study using growing lambs, substituting DDGS for SBM in an isonitrogenous diet increased dietary EE, neutral detergent fiber (NDF) and acid detergent fiber (ADF) contents, and significantly altered rumen fermentation parameters such as volatile fatty acid (VFA) and ammonia concentrations .
Within the rumen ecosystem, the environments create three different ruminal ecological niches (REN), namely the solid, liquid, and epithelium [14, 15]. The three REN have different community structure [15,16,17,18] and rumen function [17, 19]. Previous studies have compared the bacterial community occupying the three REN using polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE) [17, 20,21,22], cloning [16, 21] or high-throughput sequencing [15, 18, 19, 23,24,25]. To our knowledge, most of these studies using dairy cattle or steers as experimental animals, and only one study, which used PCR-DGGR , investigated the bacterial community occupying the different REN of growing lambs. Recently, we also evaluated the changes of the bacterial communities at the three different REN of growing Hu lambs in response to dietary urea supplementation . However, most of the studies focused only on the community structure and function of rumen bacteria, without considering the eukaryotes . This study, by integrating real-time quantitative polymerase chain reaction (qPCR), high-throughput sequencing, and functional prediction, investigated the effect of dietary replacement of SBM with DDGS on community structure and function of the microbiota at the three different REN of growing Hu lambs.
Animals, diets and experimental design
The experimental design, diets, and management have been reported previously . Briefly, The DDGS replaced all the SBM and a portion of the ground corn in the diets. The feeding trial was conducted for 10 weeks, with the first 1 week for adaptation followed by 9 weeks of dietary treatment. At the end of week 9, six lambs each were randomly selected only from the SBM-Control and the DDGS-Control groups and slaughtered for sample collection.
On days 6 and 7 of week 9, the six lambs each selected from the SBM-Control and DDGS-Control were slaughtered at 4–6 h after morning feeding. The whole rumen contents of each Hu lamb were first homogenized, mixed, and then strained through four layers of cheesecloth to separate the rumen liquid (RL) and the rumen solid (RS) fractions. Approximately 30 mL each of liquid and the solid fraction was collected into a sterilized tube and immediately stored in liquid nitrogen. To collect the rumen epithelial (RE) samples, the rumen walls were rinsed with cold sterile saline solution three times after the removal of the rumen contents. Epithelial tissue samples from an approximately 8 cm2 area of the rumen epithelium were scraped using a sterilized glass slide and stored in liquid nitrogen immediately after collection. The RS, RL, and RE samples were stored at − 80 °C until further analysis.
Metagenomic DNA of the rumen solid, liquid, and epithelium samples was extracted using the bead-beating and phenol-chloroform extraction method as previously described . The quality of the DNA extracts was visually checked using electrophoresis on 1.2% agarose gel (w/v) containing Goldview™ (SaiBaiSheng, Shanghai, China), and the DNA concentration of each sample was determined using a Nanodrop 2000 (Thermo Fisher Scientific, Inc., Madison, USA). The DNA samples were stored at − 20 °C until analyses.
Illumina sequencing of 16S rRNA gene amplicons and data analysis
The V3-V4 hypervariable regions of the 16S rRNA gene were amplified using primers 338F (5′-ACTCCTACGGGAGGCAGCA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). Unique barcodes were added to the 5′ end of both primers for multiplexing. PCR products were verified on agarose gel (2%, w/v), and the expected bands were each extracted and purified using the AxyPrepDNA Gel Extraction Kit (Axygen Biosciences, CA, USA). The concentrations of the purified DNA amplicons were each quantified using a QuantiFluor® dsDNA kit (Promega, Madison, WI, USA). Amplicons from different samples were mixed in equal ratio and sequenced using the 2 × 300 paired-end kit on an Illumina MiSeq platform. The raw sequence reads were deposited into the NCBI Sequence Read Archive (SRA) database under the accession number PRJNA565493.
Raw fastq files were de-multiplexed, quality-filtered, and analyzed using QIIME 1.9.1 . Operational taxonomic units (OTUs) were de novo clustered with a 97% sequence similarity cutoff using UPARSE (version 7.1 http://drive5.com/uparse/), and possible chimeric sequences were identified and removed using UCHIME . The most abundant sequence within each OTU was selected as the ‘representative sequence’ and was taxonomically classified based on the SILVA database (version 128) . A PH Lane mask supplied by QIIME was used to remove the hypervariable regions from the aligned sequences. FASTTREE  was used to create a phylogenetic tree of the representative sequences for each sample. Sequences identified as of chloroplasts or mitochondria were removed before further analysis. Alpha diversity measurements including observed OTUs, Chao1 richness estimate, and Shannon diversity index, as well as Good’s coverage, were calculated using QIIME 1.9.1 . Principal coordinates analysis (PCoA) was performed based on weighted UniFrac distance and Bray-Curtis dissimilarity to reveal overall differences in prokaryotic communities among the rumen solid, liquid and epithelial fractions from SBM and DDGS groups (SBM-RS, SBM-RL, SBM-RE, DDGS-RS, DDGS-RL, and DDGS-RE). The functional profiles of the rumen microbiota from different REN of the growing Hu lambs were predicted from the 16S rRNA gene data using Tax4Fun . The functional profiles were summarized at hierarchy level 2 of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways.
Quantitative real-time PCR analysis
The PCR primers used for real-time qPCR of total bacteria , fungi , protozoa , methanogens , and SRB  are listed in Table S1. Real-time PCR was performed on a StepOnePlus system (Applied Biosystems, California, USA) using the SYBR Premix Ex Taq dye (Takara Bio Inc.). Copies of 16S rRNA gene (total bacteria), 18S rRNA gene (fungi and protozoa), methyl coenzyme-M reductase alpha subunit gene (mcrA, for methanogens), and dissimilatory sulfite reductase alpha subunit gene (dsrA, for SRB) in each sample was quantified in three technical replicate against respective standards, which were purified PCR products of known length and concentration. The absolute abundance of each microbial population was expressed as copies of the target gene/g of samples.
Analysis of similarity (ANOSIM) was used to assess the statistical significance of the PCoA analysis of overall microbiota across the treatments. The real-time PCR data were log transformed to improve normality. Residual analysis was used to determine if a transformation of variables was needed. If needed, cubic root transformations were performed. All data (absolute abundance quantified by qPCR, alpha diversity measurements, relative abundances of microbial populations at the phylum and genus levels, and the relative abundance of level 2 KEGG pathways) were analyzed using the MIXED procedure of SAS version 9.4 (SAS Institute Inc., Cary, NC). The model included PS and REN as fixed effects, with block as random effects. Degrees of freedom were calculated using the Kenward-Roger option. The mean separation test was performed using the Tukey multiple range test. Differences were considered statistically significant at P ≤ 0.05.
Effect on bacterial alpha diversity
Across all 36 samples, a total of 1,275,639 quality-checked 16S rRNA gene sequences were obtained. On average, each sample had at least 35,434 sequences. Greater than 99.3% depth coverage was achieved for all the treatments (Table 1). There were no interactions (P ≥ 0.06) between PS and REN for OTU numbers, Chao 1, or Shannon index. None of the alpha diversity measurements was affected (P ≥ 0.36) by dietary PS. The number of OTUs and Chao 1 richness estimate in the RS and the RL fractions were similar (P > 0.05), which were higher (P < 0.05) than in the RE fraction. In contrast, The Shannon index was higher (P > 0.05) in the RS than in the RL and the RE fraction.
Effect on overall bacteria
The PCoA analysis, based on either Weighted UniFrac or Bray-Curtis dissimilarity, indicated a clear separation between the epithelium and the liquid or solid fraction, but no dietary PS effects were observed (Fig. 1). Moreover, the samples from the RS and the RL fractions clustered closely together, indicating very similar community structures. These findings were further supported by the ANOSIM based on the Bray-Curtis distance (P < 0.001; Figure S1).
Effect on bacteria at phylum and genus levels
In total 19 bacterial phyla were identified among all the treatments, with Bacteroidetes, Firmicutes, Proteobacteria, Spirochaetes, Fibrobacteres, and Actinobacteria being the six most predominant phyla, each of which was represented by more than 0.5% of the total sequences in at least one treatment (Table 2). No interaction (P ≥ 0.59) of PS with REN was detected with respect to any of the bacterial phyla. None of the bacterial phyla was affected (P ≥ 0.11) by dietary PS, but all were greatly influenced by REN. The relative abundance of Bacteroidetes (52.15–55.01%) in the RS and the RL was higher (P < 0.01) than that in the RE fraction (27.41%). On the contrary, the relative abundance of Proteobacteria (0.64–0.93%) and Spirochaetes (0.76–2.45%) in the RS and the RL was lower (P < 0.01) than that in the RE fraction (18.06% and 6.23%). The relative abundance of Fibrobacteres was higher (1.16%) in the RS than in the RL and the RE fractions (0.12–0.30%). The relative abundance of Firmicutes (41.41–45.63%) and Actinobacteria (0.65–0.85%) were not influenced (P ≥ 0.55) by REN.
A total of 238 bacterial genera were identified among all the treatments, but only 49 of them were each represented by more than 0.5% of the total sequences in at least one treatment (Table 3) and they were regarded as the “major genera”. No interaction (P ≥ 0.13) of PS with REN was detected on any of the bacterial abundance at the genus level except for the genera Pseudobutyrivibrio and Suttonella (P < 0.01). The relative abundance of Ruminococcaceae UCG-005 was higher (P < 0.01) for the lambs receiving SBM than those fed DDGS, while the relative abundance of other major bacterial genera was not influenced (P ≥ 0.08) by PS. In contrast, 39 out of the 49 major genera were significantly different (P < 0.05) among the three REN. Compared with the RS and the RL fraction, the relative abundance of 17 genera (including Butyrivibrio 2, Prevotellaceae UCG-001, Rikenellaceae U29-B03, Rikenellaceae Blvii28, Lachnospiraceae UCG-008, Defluviitaleaceae UCG-011, Eubacterium nodatum group, Erysipelotrichaceae uncultured, Lachnospiraceae UCG-010, Howardella, Halomonas, Desulfobulbus, Comamonas, Suttonella, Neisseriaceae uncultured, Campylobacter, and Treponema 2) was higher (P < 0.01), while that of five genera (including Prevotella 1, Prevotellaceae UCG-003, Christensenellaceae R-7 group, Ruminococcaceae NK4A214 group, and Lachnospiraceae XPB1014 group) was lower (P < 0.01) in the RE fraction than in the RS and the RL fractions. Seven genera, including Bacteroidales S24-7 group_norank, Succiniclasticum, Ruminococcus 1, Saccharofermentans, Eubacterium ruminantium group, Lachnospiraceae probable genus 10, and Fibrobacter, had a greater relative abundance (P < 0.01) in the RS than in the RL and RE fractions. In contrast, the RL fraction had a higher (P < 0.01) relative abundance of Bacteroidales BS11 gut group_norank, Bacteroidales RF16 group_norank, Ruminococcus 2, Eubacterium coprostanoligenes group, and Pseudobutyrivibrio, but a lower (P < 0.01) relative abundance of Ruminococcaceae UCG-014, than the RS and the RE fractions.
Effect on inferred functional pathways
The top 20 predominant level 2 metabolic pathways were shown in Table 4. No interaction (P ≥ 0.10) of PS with REN was detected for any of the pathways. However, the relative abundance of the predominant pathways was altered differently by dietary PS or REN. Among the pathways, 9 out of 20 were significantly influenced by PS. Compared with the SBM-based diets, the DDGS-based diets resulted in higher (P ≤ 0.02) relative abundance of the functions related to amino acid metabolism, lipid metabolism, and metabolism of terpenoids and polyketides, but lower (P ≤ 0.04) relative abundance of the functions related to nucleotide metabolism, replication and repair, translation, cell motility, folding, sorting and degradation, infectious disease (bacterial), and cell growth and death. In contrast, 16 out of the predominant 20 pathways were significantly different (P < 0.05) among the three REN. Compared with the RS and RL, the RE fraction had greater (P < 0.01) predominance of the pathways related to amino acid metabolism, lipid metabolism, xenobiotics biodegradation and metabolism, and metabolism of terpenoids and polyketides, but lower (P < 0.01) predominance of the function related to membrane transport, signal transduction, and cell motility. The pathways related to replication and repair, biosynthesis of other secondary metabolites, and endocrine system were higher (P ≤ 0.02) in the RS than in the RL and RE fractions. In contrast, a higher (P < 0.01) relative abundance of the pathways related to the metabolism of cofactors and vitamins, and a lower (P < 0.01) relative abundance of the pathway related to the metabolism of other amino acids were observed in the RL than in the RS or the RE fractions.
Effect on total bacteria, sulfate-reducing bacteria, methanogens, fungi, and protozoa
No interaction (P ≥ 0.16) between PS and REN was detected with respect to the absolute abundance of total bacteria, SRB, methanogens, fungi, or protozoa (Table 5). The DDGS-fed lambs had a smaller (P = 0.04) population of protozoa than those fed SBM, but the population of total bacteria, SRB, fungi, and methanogens were similar (P ≥ 0.41). The populations of total bacteria and methanogens were lower (P < 0.05) in the RL fraction than in the RS fraction, but higher (P < 0.05) than in the RE fraction. The RS had a larger (P < 0.05) fungal population than the RL and the RE fractions, but no difference (P > 0.05) in the fungal population was noted between the RL and the RE fractions. Compared with the RS and the RL fractions, the RE fraction had a larger (P < 0.05) SRB population but a smaller (P < 0.05) protozoal population. However, no difference (P > 0.05) was observed in the populations SRB or protozoa between the RS and the RL fractions.
Effects of dietary protein source on rumen microbiota structure and function
Diet is one of the main factors influencing the structure and function of rumen microbiota [2, 3]. In our previous feeding trial, replacing SBM in an isonitrogenous lamb diet with 20% DDGS increased dietary EE, NDF and ADF contents . Therefore, we expected shifts of rumen bacterial community and structure. However, high throughput sequencing results (Tables 1, 2 and 3) revealed very limited influence of dietary PS on the bacterial community structure. These results are similar to a previous study that reported no or little difference in microbiota diversity and relative abundance of most bacteria in the rumen of crossbred steers fed 19.5% DDGS replacing corn bran  or dairy cows fed 20% DDGS replacing SBM . In contrast, Callaway et al.  and Ramirez-Ramirez et al.  found drastic changes in rumen bacterial community structure when DDGS were added at 30% and 50% of the diet [dry matter (DM)]. The disparity in bacterial community structure among different studies may be due to differences in DDGS feeding level and the chemical composition of the diets . Additionally, the difference of ruminant species, duration of experimental period, sequencing and data analysis methods may also contribute to the discrepancy in the bacterial community among studies.
Although the bacterial community structure in the rumen did not differ much between the two different protein sources, the metabolic pathways predicted from the 16S rRNA gene sequences varied greatly (Table 4). These results demonstrated that the substitution of SBM with DDGS had a greater impact to the rumen microbial function than to the bacterial composition. It is also suggested that small changes in the rumen bacterial community may lead to greater variations in metabolic pathways. Among the predicted metabolic pathways, the increase of lipid metabolism pathway in the DDGS treatment group might be related to the higher EE content in the diet. The changes in other metabolic pathways, such as amino acid metabolism, metabolism of terpenoids and polyketides, nucleotide metabolism, replication and repair, translation, cell motility, might also be related to changes in dietary composition, but the specific mechanism is still unclear. The application of multiple meta-omics, such as metagenomics, metatranscriptomics, metaproteomics, and metabolomics  in future studies may help further explain the effect of different protein sources on rumen metabolism.
Feed digestion and fermentation are the concerted functions of a variety of microorganisms including bacteria, fungi, protozoa, archaea and phages . Therefore, the change of rumen fermentation parameters found in our previous feeding trial  might be closely associated with shifts of members of the rumen microbiota besides bacteria. This study showed that the populations of rumen total bacteria, fungi, SRB, and methanogens were not influenced by the dietary protein sources, but the population of protozoa was reduced by the DDGS substitution for SBM (Table 5). Protozoa can account for up to 50% of the rumen biomass and play an important role in the degradation of dietary fiber and protein . Moreover, protozoa are also known to increase ammonia production . In a previous meta-analysis, it has been reported that the elimination of protozoa from the rumen significantly decreased VFA and ammonia concentration . Therefore, the decreased protozoal population might have primarily contributed to the reduced VFA and ammonia concentration in the DDGS group . These findings indicate that, in addition to bacteria, the structure and function of other microbes such as protozoa, fungi, methanogens, and phages should also be studied in nutritional research.
Effects of ruminal ecological niche on microbiota structure and function
Consistent with previous studies using PCR-DGGE  or high-throughput sequencing [15, 18, 26], RS and RL shared more bacteria when compared to RE. In addition, the predicted functions and the qPCR results of the microbial groups further demonstrated that there were great differences in the structure and function of the microbiota occupying different REN. This corroborates the niche partitioning of the rumen microbiota.
Differences in microbiota between rumen solid and liquid fractions
The Shannon diversity index in the RS is higher than in the RL (Table 1), which agrees with the findings of previous studies [15, 18]. The bacterial richness and taxonomic composition in the RS and RL were similar, suggesting continuous exchange between these two fractions . However, some taxa had different relative abundance in these two niches, a finding consistent with other reports [15, 18, 26] and likely reflecting specialized niches related to digestion of soluble components vs. dietary fiber .
Rumen bacteria found in the RS fraction are mainly responsible for the initial and secondary degradation of feed and play an important role in fiber digestion, while the bacteria in the RL fraction are mainly involved in the fermentation of soluble nutrients and metabolic end products of feed digestion [15, 41]. Ruminococcus and Fibrobacter are the two main known fibrolytic bacterial genera in the rumen . We indeed found a higher predominance of Ruminococcus 1 and Fibrobacter spp. in the RS than in the RL fractions of the growing lambs (Table 3). Some other genera, such as Bacteroidales S24-7 group, Succiniclasticum, Saccharofermentans, Eubacterium ruminantium group, and Lachnospiraceae probable genus 10, were also more predominant in the RS than in the RL, which indicates that these bacteria may play an important role in the initial and secondary degradation of feed. In contrast, some other bacterial genera, such as Bacteroidales BS11 gut group, Bacteroidales RF16 group, Ruminococcus 2, Eubacterium coprostanoligenes group, and Pseudobytyrivibrio showed the opposite trend, which indicates that these bacteria probably are mainly involved in the catabolism of soluble nutrients. Moreover, the results of functional analysis further demonstrated that there were great differences in multiple bacterial metabolic pathways between the RS and the RL.
It should be noted that bacteria with similar relative abundance in the RS and the RL fraction may differ in their absolute abundance. In the present study, the population of total bacteria in the RS is much higher than that in the RL (Table 5), indicating that the population of bacteria might be greater in the RS than in the RL. In addition to bacteria, other rumen microorganisms such as protozoa, fungi and methanogens also play an important role in feed digestion . However, recent studies mainly focused on the diversity and function of rumen bacteria, ignoring the structure and function of other microbial microbes . The rumen fungi can account for up to 20% of the microbial biomass and play an important role in ruminal fiber degradation . Methanogens are the main hydrogen utilizing microorganisms in the rumen, and hydrogen produced by hydrogen-producing bacteria such as cellulolytic bacteria and fungi are used by methanogens to reduce CO2 to methane . Therefore, it is not surprising that the population of fungi and methanogens is higher in the RS than in the RL. In contrast, the similar population of protozoa between the RS and the RL in the rumen is probably due to the ability of protozoa to freely attach to and dissociate from feed particles . Generally, the population of protozoa is more abundant than fungi in rumen . In the present study, similar results was also found in the RL fraction. However, the population of fungi in RS fraction is much higher than protozoa. Therefore, the disparity with previous reports in the population of rumen fungi and protozoa may be due to differences in REN.
Differences in microbiota between epithelial and the solid or the liquid fractions
In line with previous studies [15, 18, 26], the RE bacterial diversity was lower than that in the RS or the RL fraction (Table 2), and the bacterial composition was greatly distinct from that of the RS or the RL (Table 3). At the phylum level, contrary to the finding in the RS or the RL, Firmicutes was the first while Bacteroidetes the second largest phyla in the RE microbiota. It is worth noting that Proteobacteria was the third largest phylum (with a relative abundance of 18.06%) of the RE microbiota, which is consistent with previous studies on dairy cattle [19, 24]. Members of Proteobacteria are mostly facultative anaerobes . Therefore, the higher predominance of Proteobacteria on RE could be explained by the trace amounts of oxygen diffused through the rumen tissue . Interestingly, all the dominant bacterial genera classified to Proteobacteria in the epithelial microbiota, such as Campylobacter, Desulfobulbus, Neisseriaceae uncultured, Comamonas, and Halomonas were barely detected in the RS or the RL. In addition to bacteria of Proteobacteria, some bacterial genera assigned to other phyla were predominant on the RE (e.g., Prevotellaceae UCG-001, Butyrivibrio 2, and Treponema 2), and again, they are virtually undetected or only detected at much lower relative abundance in the RS and or the RL. These results support the notion that the distinctive epithelial bacteria assigned to different phyla may have additional functions other than feed digestion .
Corroborating the finding of previous studies [26, 48, 49], Desulfobulbus, one genus of SRB , was the predominant SRB genus (5.22%) in RE microbiota (Table 3). In line with the high-throughput sequencing result, the qPCR results also demonstrated that SRB were more predominant on RE than in the RS or the RL fraction. Previous studies have shown that sulfides produced by SRB can disrupt the gut epithelial tissues and induce DNA damage, and adversely affect host health [51, 52]. However, it was found that the relative abundance of Butyrivibrio 2 (15.51%) on RE was about 3 times higher than that of Desulfobulbus. Members of Butyrivibrio (e.g., Butyrivibrio fibrisolvens) are important butyrate producing bacteria ubiquitous in the rumen . In vitro cell culture experiments have shown that butyrate could regulate colonic proliferation and treat ulcerative colitis caused by sulfide . Therefore, it is also reasonable to believe that butyrate, a metabolite of Butyrivibrio spp., may help repair the damage to the RE caused by the sulfide produced by SRB (e.g., Desulfobulbus spp.). Based on this, we speculate that there is a self-regulating mechanism between the epithelial microbes and the host that helps the homeostasis of RE function.
Rumen epithelial microbiota play an important role in digesting and recycling the keratinized distal cells of the epithelium . In the present study, the enhanced pathway of amino acid metabolism and lipid metabolism in the RE microbiota compared with the RS and the RL microbiota (Table 4) are probably related with the special nutrient composition of the keratinized distal cells of the epithelium. One recent metatranscriptomic analysis of RE microbiota showed that many metabolic genes encoding enzymes involved in N metabolism such as glutamate dehydrogenase, glutamine synthase and glutamate synthase were highly expressed, which demonstrated the importance of RE bacteria in N metabolism . Similarly, Mao et al.  also found enhanced amino acid metabolism in the RE microbiota compared with the microbiota of rumen content. In the present study, the qPCR results confirmed the existence of fungi, protozoa and methanogens in the epithelial microbiota. However, Mann et al.  detected fungi and methanogens but no protozoa on the RE in their metatranscriptomic analysis. The inconsistent results might be explained by the more sensitive detection by qPCR than by RNA-seq.
This study, combing high-throughput sequencing, functional prediction, and qPCR, for the first time, revealed the difference of structure and function of rumen microbiota occupying different REN of growing Hu lambs in response to alterations in dietary PS. The results of the present study indicate that substitution of SBM with DDGS had greater impact to protozoa than to other microbes, and the microbiota structure and function in different REN are specific and niche-adapted. In order to have a better understanding of the complex rumen ecosystem, in addition to bacteria, the structure and function of other microbes such as protozoa, fungi, and methanogens should also be studied.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its supplementary information files].
Acid detergent fiber
Analysis of similarity for multivariate
Dried distillers grains with solubles
Dissimilatory sulfite reductase alpha subunit gene
Kyoto Encyclopedia of Genes and Genomes
Methyl coenzyme-M reductase alpha subunit gene
Neutral detergent fiber
Operational taxonomic units
Principal coordinates analysis
Polymerase chain reaction-Denaturing gradient gel electrophoresis
Real-time quantitative polymerase chain reaction
Ruminal ecological niche
Total mixed ration
Volatile fatty acid
Cammack KM, Austin KJ, Lamberson WR, Conant GC, Cunningham HC. RUMINANT NUTRITION SYMPOSIUM: Tiny but mighty: the role of the rumen microbes in livestock production. J Anim Sci 2018;96(10):4481–4481.
Loor JJ, Elolimy AA, McCann JC. Dietary impacts on rumen microbiota in beef and dairy production. Anim Front. 2016;6(3):22–9.
Wolff SM, Ellison MJ, Hao Y, Cockrum RR, Austin KJ, Baraboo M, et al. Diet shifts provoke complex and variable changes in the metabolic networks of the ruminal microbiome. Microbiome. 2017;5:60.
Kleinschmit DH, Anderson JL, Schingoethe DJ, Kalscheur KF, Hippen AR. Ruminal and intestinal degradability of distillers grains plus solubles varies by source. J Dairy Sci. 2007;90(6):2909–18.
Shen JS, Chen YY, Moraes LE, Yu ZT, Zhu WY. Effects of dietary protein sources and nisin on rumen fermentation, nutrient digestion, plasma metabolites, nitrogen utilization, and growth performance in growing lambs. J Anim Sci. 2018;96(5):1929–38.
Held J: Feeding soy hulls and dried distillers grain with solubles to sheep. South Dakota State University, Extension Extra Archives. 2006. http://openprairie.sdstate.edu/extension_extra/75. Accessed 1 Jan 2006.
Klopfenstein TJ, Erickson GE, Bremer VR. BOARD-INVITED REVIEW: Use of distillers by-products in the beef cattle feeding industry. J Anim Sci 2008;86(5):1223–1231.
Schingoethe DJ, Kalscheur KF, Hippen AR, Garcia AD. Invited review: the use of distillers products in dairy cattle diets. J Dairy Sci. 2009;92(12):5802–13.
Callaway TR, Dowd SE, Edrington TS, Anderson RC, Krueger N, Bauer N, et al. Evaluation of bacterial diversity in the rumen and feces of cattle fed different levels of dried distillers grains plus solubles using bacterial tag-encoded FLX amplicon pyrosequencing. J Anim Sci. 2010;88(12):3977–83.
Aldai N, Klieve AV, Dugan MER, Kramer JKG, Ouwerkerk D, Aalhus JL, et al. Evaluation of rumen fatty acid hydrogenation intermediates and differences in bacterial communities after feeding wheat- or corn-based dried distillers grains to feedlot cattle. J Anim Sci. 2012;90(8):2699–709.
Castillo-Lopez E, Ramirez Ramirez HA, Klopfenstein TJ, Anderson CL, Aluthge ND, Fernando SC, et al. Effect of feeding dried distillers grains with solubles on ruminal biohydrogenation, intestinal fatty acid profile, and gut microbial diversity evaluated through DNA pyro-sequencing. J Anim Sci. 2014;92(2):733–43.
Ramirez HAR, Nestor K, Tedeschi LO, Callaway TR, Dowd SE, Fernando SC, et al. The effect of brown midrib corn silage and dried distillers’ grains with solubles on milk production, nitrogen utilization and microbial community structure in dairy cows. Can J Anim Sci. 2012;92(3):365–80.
Castillo-Lopez E, Jenkins CJR, Aluthge ND, Tom W, Kononoff PJ, Fernando SC. The effect of regular or reduced-fat distillers grains with solubles on rumen methanogenesis and the rumen bacterial community. J Appl Microbiol. 2017;123(6):1381–95.
Cheng KJ, McCowan RP, Costerton JW. Adherent epithelial bacteria in ruminants and their roles in digestive tract. Am J Clin Nutr. 1979;32(1):139–48.
Mulder TD, Goossens K, Peiren N, Vandaele L, Haegeman A, Tender CD, et al. Exploring the methanogen and bacterial communities of rumen environments: solid adherent, fluid and epimural. FEMS Microbiol Ecol. 2017;93(3):fiw251. https://doi.org/10.1093/femsec/fiw251.
Cho SJ, Cho KM, Shin EC, Lim WJ, Hong SY, Choi BR, et al. 16S rDNA analysis of bacterial diversity in three fractions of cow rumen. J Microbiol Biotechnol. 2006;16(1):92–101.
Sadet S, Martin C, Meunier B, Morgavi DP. PCR-DGGE analysis reveals a distinct diversity in the bacterial population attached to the rumen epithelium. Animal. 2007;1(7):939–44.
Scharen M, Kiri K, Riede S, Gardener M, Meyer U, Hummel J, et al. Alterations in the rumen liquid-, particle- and epithelium-associated microbiota of dairy cows during the transition from a silage- and concentrate-based ration to pasture in spring. Front Microbiol. 2017;8:744.
Mao SY, Zhang ML, Liu JH, Zhu WY. Characterising the bacterial microbiota across the gastrointestinal tracts of dairy cattle: membership and potential function. Sci Rep. 2015;5:16116.
Chen YH, Oba M, Guan LL. Variation of bacterial communities and expression of toll-like receptor genes in the rumen of steers differing in susceptibility to subacute ruminal acidosis. Vet Microbiol. 2012;159(3–4):451–9.
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.
Malmuthuge N, Li MJ, Chen YH, Fries P, Griebel PJ, Baurhoo B, et al. Distinct commensal bacteria associated with ingesta and mucosal epithelium in the gastrointestinal tracts of calves and chickens. FEMS Microbiol Ecol. 2012;79(2):337–47.
Pitta DW, Pinchak WE, Dowd SE, Osterstock J, Gontcharova V, Youn E, et al. Rumen bacterial diversity dynamics associated with changing from Bermudagrass hay to grazed winter wheat diets. Microb Ecol. 2010;59(3):511–22.
Liu JH, Zhang ML, Zhang RY, Zhu WY, Mao SY. Comparative studies of the composition of bacterial microbiota associated with the ruminal content, ruminal epithelium and in the faeces of lactating dairy cows. Microb Biotechnol. 2016;9(2):257–68.
Ren QM, Si HZ, Yan XT, Liu C, Ding LM, Long RJ, et al. Bacterial communities in the solid, liquid, dorsal, and ventral epithelium fractions of yak (Bos grunniens) rumen. Microbiologyopen. 2020;9(2):e963.
Li ZP, Mu CL, Xu YX, Shen JS, Zhu WY. Changes in the solid-, liquid-, and epithelium-associated bacterial communities in the rumen of hu lambs in response to dietary urea supplementation. Front Microbiol. 2020;11:244.
Huws SA, Creevey CJ, Oyama LB, Mizrahi I, Denman SE, Popova M, et al. Addressing global ruminant agricultural challenges through understanding the rumen microbiome: past, present, and future. Front Microbiol. 2018;9:2161.
Dai ZL, Zhang J, Wu GY, Zhu WY. Utilization of amino acids by bacteria from the pig small intestine. Amino Acids. 2010;39(5):1201–15.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.
Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26(7):1641–50.
Asshauer KP, Wemheuer B, Daniel R, Meinicke P. Tax4Fun: predicting functional profiles from metagenomic 16S rRNA data. Bioinformatics. 2015;31(17):2882–4.
Denman SE, McSweeney CS. Development of a real-time PCR assay for monitoring anaerobic fungal and cellulolytic bacterial populations within the rumen. FEMS Microbiol Ecol. 2006;58(3):572–82.
Sylvester JT, Karnati SKR, Yu ZT, Morrison M, Firkins JL. Development of an assay to quantify rumen ciliate protozoal biomass in cows using real-time PCR. J Nutr. 2004;134(12):3378–84.
Denman SE, Tomkins NW, McSweeney CS. Quantitation and diversity analysis of ruminal methanogenic populations in response to the antimethanogenic compound bromochloromethane. FEMS Microbiol Ecol. 2007;62(3):313–22.
Devkota S, Wang YW, Musch MW, Leone V, Fehlner-Peach H, Nadimpalli A, et al. Dietary-fat-induced taurocholic acid promotes pathobiont expansion and colitis in Il10(−/−) mice. Nature. 2012;487(7405):104–8.
Williams AG, Coleman GS. The rumen protozoa. In: Hobson PN, Stewart CS, editors. The rumen microbial ecosystem. NY: Springer; 1997. p. 73–139.
Firkins JL, Yu Z, Morrison M. Ruminal nitrogen metabolism: perspectives for integration of microbiology and nutrition for dairy. J Dairy Sci. 2007;90(Suppl 1):E1–16.
Newbold CJ, de la Fuente G, Belanche A, Ramos-Morales E, McEwan NR. The role of ciliate protozoa in the rumen. Front Microbiol. 2015;6:1313.
McAllister TA, Bae HD, Jones GA, Cheng KJ. Microbial attachment and feed digestion in the rumen. J Anim Sci. 1994;72(11):3004–18.
Krause DO, Denman SE, Mackie RI, Morrison M, Rae AL, Attwood GT, et al. Opportunities to improve fiber degradation in the rumen: microbiology, ecology, and genomics. FEMS Microbiol Rev. 2003;27(5):663–93.
Edwards JE, Forster RJ, Callaghan TM, Dollhofer V, Dagar SS, Cheng Y, et al. PCR and omics based techniques to study the diversity, ecology and biology of anaerobic fungi: insights, challenges and opportunities. Front Microbiol. 2017;8:1657.
Jeyanathan J, Martin C, Morgavi DP. The use of direct-fed microbials for mitigation of ruminant methane emissions: a review. Animal. 2014;8(2):250–61.
Choudhury PK, Salem AZM, Jena R, Kumar S, Sigh R, and Puniya AK. Rumen microbiology: an overview In: Puniya AK, Singh R, Kamra DN, editors. The Rumen microbiology: from evolution to revolution. New Delhi: Springer. 2015. p. 3–16.
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.
Liu JH, Bian GR, Zhu WY, Mao SY. High-grain feeding causes strong shifts in ruminal epithelial bacterial community and expression of toll-like receptor genes in goats. Front Microbiol. 2015;6:167.
Jiao JZ, Huang JY, Zhou CS, Tan ZL. Taxonomic identification of ruminal epithelial bacterial diversity during rumen development in goats. Appl Environ Microb. 2015;81(10):3502–9.
Liu J, Bian G, Sun D, Zhu W, Mao S. Starter feeding altered ruminal epithelial bacterial communities and some key immune-related genes’ expression before weaning in lambs. J Anim Sci. 2017;95(2):910–21.
Vos P, Garrity G, Jones D, Krieg NR, Ludwig W, Rainey FA, et al. Bergey’s manual of systematic bacteriology. Vol 3, The Firmicutes New York, NY: Springer Science & Business Media; 2009.
Attene-Ramos MS, Wagner ED, Plewa MJ, Gaskins HR. Evidence that hydrogen sulfide is a genotoxic agent. Mol Cancer Res. 2006;4(1):9–14.
Yang G. Hydrogen sulfide in cell survival: a double-edged sword. Expert Rev Clin Pharmacol. 2011;4(1):33–47.
Stewart CS, Flint HJ, Bryant MP. The rumen bacteria. In: Hobson PN, Stewart CS, editors. The rumen microbial ecosystem. NY: Springer; 1997. p. 10–72.
Christl SU, Eisner HD, Dusel G, Kasper H, Scheppach W. Antagonistic effects of sulfide and butyrate on proliferation of colonic mucosa: a potential role for these agents in the pathogenesis of ulcerative colitis. Digest Dis Sci. 1996;41(12):2477–81.
Mann E, Wetzels SU, Wagner M, Zebeli Q, Schmitz-Esser S. Metatranscriptome sequencing reveals insights into the gene expression and functional potential of rumen wall bacteria. Front Microbiol. 2018;9:43.
The authors thank all of the staff of the Shanghai Yonghui Sheep Industry Co., Ltd. for their assistance in animal care and sample collection.
This work was partially supported by the National Natural Science Foundation of China (award No: 31402101), the “One Belt and One Road” Technical Cooperation Project of Jiangsu Province (award No: BZ2018055), and the Jiangsu Agriculture Science and Technology Innovation Fund (award No: CX(19)3023).
Ethics approval and consent to participate
The experimental procedures used in this study were approved by the Animal Care and Use Committee of Nanjing Agricultural University (protocol number: SYXK2017–0007) .
Consent for publication
The authors declare no competing interests.
. Primers used for real-time PCR quantification of rumen target organisms. Figure S1. Box plots showing within-group similarity and between-group dissimilarity of rumen microbiota based on Bray-Curtis dissimilarity in different ruminal ecological niches of growing Hu lambs in response to alterations in dietary protein sources. The different letters denote significant differences (Kruskal-Wallis tests, FDR-adjusted q < 0.05). DDGS: dried distillers grains with solubles, SBM: soybean meal, RE: rumen epithelium, RS: rumen solid, RL: rumen liquid.
About this article
Cite this article
Shen, J., Li, Z., Yu, Z. et al. Effects of dietary replacement of soybean meal with dried distillers grains with solubles on the microbiota occupying different ecological niches in the rumen of growing Hu lambs. J Animal Sci Biotechnol 11, 93 (2020). https://doi.org/10.1186/s40104-020-00499-2
- Distillers dried grains with solubles
- Growing lamb
- Ruminal ecological niche
- Soybean meal