- Open Access
Comparison of rumen bacteria distribution in original rumen digesta, rumen liquid and solid fractions in lactating Holstein cows
Journal of Animal Science and Biotechnology volume 8, Article number: 16 (2017)
Original rumen digesta, rumen liquid and solid fractions have been frequently used to assess the rumen bacterial community. However, bacterial profiles in rumen original digesta, liquid and solid fractions vary from each other and need to be better established.
To compare bacterial profiles in each fraction, samples of rumen digesta from six cows fed either a high fiber diet (HFD) or a high energy diet (HED) were collected via rumen fistulas. Rumen digesta was then squeezed through four layers of cheesecloth to separate liquid and solid fractions. The bacterial profiles of rumen original digesta, liquid and solid fractions were analyzed with High-throughput sequencing technique.
Rumen bacterial diversity was mainly affected by diet and individual cow (P > 0.05) rather than rumen fraction. Bias distributed bacteria were observed in solid and liquid fractions of rumen content using Venn diagram and LEfSe analysis. Fifteen out of 16 detected biomarkers (using LEfSe analysis) were found in liquid fraction, and these 15 biomarkers contributed the most to the bacterial differences among rumen content fractions.
Similar results were found when using samples of original rumen digesta, rumen liquid or solid fractions to assess diversity of rumen bacteria; however, more attention should be draw onto bias distributed bacteria in different ruminal fractions, especially when liquid fraction has been used as a representative sample for rumen bacterial study.
The bovine rumen harbors a diverse population of microorganisms that convert ingested plant biomass into microbial protein and volatile fatty acids, and their fermentation end-products provide the host with essential nutrients for metabolism. Rumen microbes, therefore, play a key role in the productivity and health of ruminants .
Collection and sampling of ruminal content are important in both scientific research and diagnosis of diseases in ruminants . Ruminal microbial diversity has been investigated in numbers of studies using different ruminal fractions including original rumen digesta, rumen liquid or solid fractions [3–8]. It has been demonstrated that rumen sampling methods and/or sampling pre-treatments could affect on the results of rumen microbial community [9, 10]. The relationships among microbial communities in different fractions of rumen content have been studied previously [11, 12]; however, it is still a controversial topic and deserves further investigations. In late studies, development of high-throughput sequencing techniques have allowed subtle effects on microbial community components to be detected as changes in relative numbers of bacterial community . The objective of this study was to assess the differences and similarities of bacterial community in original rumen digesta, rumen liquid and solid fractions using a high-throughput sequencing technique.
Animals and sampling
Six ruminally fistulated lactating Holstein cows were housed in a free stall pen at the Zhongdi Dairy Research Center (Beijing, China) and were cared for according to the practices outlined in the Guide for the Care and Use of Agriculture Animals in Agriculture Research and Teaching (FASS, 2010).
Cows were randomly assigned to two groups with three cows in each group and individually fed by Roughage Intake Control System (RIC, Insentec B.V, Netherland). One group of cows were fed with a high fiber containing diet (HFD group) and the other group was fed with a high energy containing diet (HED group) (Table 1). After a 14-days adaptation period to the experimental diets, approximately 500 g of original rumen digesta of each cow was collected 5 h after morning feeding via rumen fistulas from the middle part of the ventral sac. After the original rumen digesta sampled, the solid and liquid fractions were obtained by squeezing the original digesta through four layers of sterile cheesecloth. All samples were snap-frozen in liquid nitrogen and were stored at −80 °C until DNA extraction.
Genomic DNA (gDNA) was extracted from 1 g of original or solid fraction, and from 1 mL of ruminal liquid with Qiagen DNA Extraction Kit™ (Qiagen, Hilden, Germany) using a repeated bead beating method followed by phenol-chloroform extraction according to the manufacturer’s protocol. The DNA was re-suspended after being precipitated with ethanol. The quality of extracted DNA was assessed based on the absorbance ratios of 260/280 nm and 260/230 nm using a NanoDrop ND-1000 Spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). The values for A260/A280 ratio in the present study were ranged from 1.8 to 2.0.
PCR amplification and purification
For illumina MiSeq sequencing, bacterial 16S rRNA gene were amplified using primers covering the V3 region (343 F, 5′-GATCCTACGGGAGGCAGCA-3′ and 534R, 5′-GCTTACCGCGGCTGCTGGC-3′) with barcodes. All PCR reactions were carried out in 30 μL reaction mixtures with 15 μL of Phusion High-Fidelity PCR Master Mix (New England Biolabs), 0.2 μL of both forward and reverse primers and 10 ng of template DNA. PCR amplification was carried out according to the following protocol: initial denaturation for 5 min at 95 °C, followed by 25 cycles of denaturation at 95 °C for 1 min, annealing at 50 °C for 1 min, and elongation at 72 °C for 1 min, with a final elongation step at 72 °C for 7 min. To qualify and quantify PCR products, the same volume of 1 × loading buffer (containing SYBR green) mixed with PCR products and electrophoresis on 2% agarose gel. Samples with bright main strip between 200 and 210 bp were then chosen for further analysis. PCR products from samples for sequencing in the same MiSeq run were pooled at equal molality. The pooled mixture was purified with a QIAquick PCR Purification Kit (Qiagen, Hilden, Germany) and re-quantified with Agilent DNA 1000 Kit (Agilent Technologies Inc.).
Sequencing with high-throughput sequencing technique
Sequencing libraries were generated using NEBNext ultra DNA sample preparation kit (NEB, USA), following standard Illumina sample-preparation protocol. The quality of library was assessed on the Qubit 2.0 Fluorometer (Life technologies, Grand Island, NY, USA) and Agilent Bioanalyzer 2100 system (Agilent Technologies, Palo Alto, Calif.). The library was sequenced on illumina MiSeq platform and ~250 bp to ~300 bp paired-end reads were generated.
Quality control of raw data and data processing
Controlling the quality of raw data was done by FastQC (version 0.11.3). Reads with quality score higher than 30 were retained for further analysis. Paired-end reads from the original DNA fragments were merged using FLASH (version 1.2.7) . Paired-end reads was assigned to each sample according to the unique barcodes. Concatenated sequences were detected using USEARCH (v6.1), and subsequently filtered out. Sequences analyses were performed using QIIME pipeline (version 1.5.0) . Generated sequences were distributed into different samples based on barcodes, and the OTUs were defined by clustering sequences together with a 97% identity cut-off at UCLUST software  after removing the barcode. The RDP classifier  was used for taxonomic classification of generated OTUs . To ensure the comparability of the species diversity between the samples, standardized OTU documents were used to analyze the species and diversity indexes. The threshold for the number of standardized sequences was set at 150,000 sequences.
Alpha diversity indices were calculated using QIIME pipeline (version 1.5.0) . Beta diversity indices between samples were determined based on Bray-Curtis metric, relationship network of each fraction was calculated using Pearson correlation. LEfSe (LDA Effect Size) analysis was performed online (https://huttenhower.sph.harvard.edu/galaxy) to find differentially abundant taxa (biomarkers) with P- value higher than 0.05 and LDA score higher than 2. ANOSIM analysis was performed with R software (version 3.1.2). Comparisons of bacterial abundance in different experimental groups were performed using Wilcoxon or Kruskal-Wallis test with R software (version 3.1.2).
Bacterial diversities and community composition in different rumen fractions
Six ruminally fistulated Holstein dairy cows were grouped by two different diets containing either high level of fiber or high level of energy. Rumen contents were collected and separated to liquid and solid fractions. After sequencing on the Illumina MiSeq platform, 4,650,173 quality reads were generated, and 23,896 OTUs in total in all samples were detected with an average of 258,343 ± 26,890 reads and 5,003 ± 678 OTUs in each sample.
Within each dietary treatment, the number of OTUs and Chao1 index between different fractions of rumen samples were similar (P > 0.05); however, the number of OTUs and Chao1 index in rumen liquid fraction was higher in HFD group than that of HED group (P < 0.05; Additional file 1: Figure S1a and b). Simpson index has no difference among fractions within and between dietary treatments (P > 0.05; Additional file 1: Figure S1c).
The microbial diversity difference was displayed by heatmap, the similarity index with Pearson correlation, and ANOSIM analysis with Bray-Curtis similarity (Fig. 1; Additional file 1: Figure S2). We found that bacteria communities were clustered mainly by diets (Fig. 1a, b; Additional file 1: Figure S2a) and individual cows (Fig. 1c; Additional file 1: Figure S2b).
Venn plot was used to illustrate the distribution of bacteria in different ruminal fractions. After filtering out the rare OTUs (defined as OUTs that only appear in one sample in each group), we identified 2,022 OTUs appeared in ruminal fractions. More unique OTUs (appeared in one fraction but not appeared in others) were detected in HFD group compared with that in HED group (Additional file 1: Figure S3a). Rumen original digesta, solid and liquid fractions shared 53% of detected OTUs, but unique OTUs was also found in different fractions although they were evenly distributed among the rumen fractions (Additional file 1: Figure S3b).
Rumen bacterial taxa change in original and fractional rumen digesta
In total, 22 phyla, 38 classes, 62 orders, 96 families, 127 genera were detected regardless of fractional and dietary treatments (Fig. 2). To identify the taxon distributions in different fractions, LEfSe analysis was performed and biomarkers of liquid and/or solid fractions were found in both HFD and HED groups (Fig. 2). Fifteen taxa were found as biomarker in HFD group (Fig. 2a; Additional file 1: Figure S4a), while only one taxon was detected as biomarker in HED group (Fig. 2b; Additional file 1: Figure S4b).
Within HFD group, 14 taxa increased in liquid fraction and one taxon increased in solid fraction (Fig. 2a; Additional file 1: Figure S5). Among the changed taxa, genus Coprococcus and Oscillospira were found to be predominant bacteria (appeared in all samples and relative abundance ≥ 1% in at least one sample; Additional file 1: Figure S6) in both rumen solid and liquid fractions in HFD group, respectively. Within HFD group, relative abundance of genus Succinivibirio among ruminal fractions were found different (Fig. 2b) with the highest abundance found in solid fraction (P < 0.05; Additional file 1: Figure S5).
Influence of rumen content fractions on bacterial diversity
Rumen, harboring large number of inhabiting microbes (approximately 1011 bacterial cells per g of rumen content) which play important roles in providing necessary nutrients (such as proteins and energy yielding substrates) to the host animal [18, 19]. Many factors such as age , diet [11, 19], and animal individual  affect rumen microbial community. Original rumen digesta, liquid fraction and solid fraction samples have been frequently used to assess rumen microbes [1, 4]; however, the fractions used may also cause biased observation in rumen microbial studies [1, 20]. In present study, we examined the differences and similarities of bacterial communities in different rumen fractions from lactating dairy cows fed with either HFD or HED using high throughput sequencing technology.
Many researchers proved that diet is one of the main factors that affect rumen microbial diversity, and that the nutritional plane and/or feed ingredients have great impacts on rumen microbial communities [3, 21–25], which may due to bacteria preferring particular metabolic substrates and rumen environment . This was supported by our findings that HFD- and HED-fed cows harbored different bacterial communities (Fig. 1a; Additional file 1: Figure S3a), which was illustrated by different α-diversity indices (Additional file 1: Figure S1) and different clusters (Fig. 1b; Additional file 1: Figure S2a) between dietary treatment groups.
Regardless the dietary effect, bacterial community of each individual cow was quite unique (Fig. 1c). A clear grouping resulted by individual cows (Additional file 1: Figure S2b) indicating that animal individual had its own distinct microbiota . This finding was also supported by a previous study in which cows were switched from a high forage diet to a high concentrate diet (during acidosis and after recovery) and the bacterial populations exhibited a low taxonomic variability as less than 5% of total identified OTUs differed .
Between different ruminal content fractions, the diversity difference was not observed (Fig. 1b and c; Additional file 1: Figure S1; Figure S2c), while different bacteria distribution was detected (Fig. 1a; Additional file 1: Figure S3b) indicating that each fraction had some unique bacteria.
Influence of ruminal fractions on bacterial taxon distribution
Different bacteria prefer particular metabolic substrates and rumen environment and thus might be distributed differently in disparate phases . In the current study, results of LEfSe illustrated that distribution difference of rumen bacteria in different fractions can be observed in both HFD group (Fig. 2a) and HED group (Fig. 2b). Biomarkers of liquid and solid fractions were detected (Fig. 2), and bacterial community of original digesta were largely likely displaying an intermediate state of liquid and solid fractions (Additional file 1: Figure S5). This verifies the common theory that bacterial community in original digesta represents the real rumen bacterial community the best.
Crucially, we found that liquid fraction contributed most to the bacteria difference in different fractions of HFD group (Fig. 2; Additional file 1: Figure S5). Usually, for studying rumen microbial diversity, original rumen digesta sampling through a fistula or from slaughtered animals was described the best method to have a representative sample . Alternatively, rumen liquid collected via a stomach tube has become a routinely used method for rumen sample collection because of its easy achievement [1, 28, 29]. A previous study showed significant differences in bacterial communities between rumen solid and liquid content . In current study, liquid fraction had similar bacterial diversity with original rumen digesta and solid fraction, but particular bacteria abundance (including predominant bacteria) in liquid fraction differed with original content and solid fraction (Additional file 1: Figure S5; Figure S6). The biased distribution of bacteria should be taken into consideration when samples of liquid fraction have been used to assess rumen bacterial community.
This study investigated and compared rumen bacteria community of original rumen digesta, liquid and solid fractions from cows fed with HFD and HED. Rumen bacteria diversity was mainly affected by diet and cow individuals rather than by rumen fractions. Bias distributed bacteria were observed in different fractions of rumen content. Liquid fraction contributed most to bacterial differences among rumen content fractions of HFD group. Results indicated that using different fractions to assess rumen bacterial diversity will generate similar results; however, more attention should be paid to bias distributed bacteria in different fractions, especially when liquid fraction is used as a representative sample.
Acid detergent fiber
Dried distillers grains with soluble
High energy diet
High fiber diet
Neutral detergent fiber
- NEL :
Net energy requirement for lactation
Operational Taxonomic Units
Henderson G, Cox F, Kittelmann S, Miri VH, Zethof M, Noel SJ, et al. Effect of DNA extraction methods and sampling techniques on the apparent structure of cow and sheep rumen microbial communities. PLoS One. 2013;8:e74787.
Shen JS, Chai Z, Song LJ, Liu JX, Wu YM. Insertion depth of oral stomach tubes may affect the fermentation parameters of ruminal fluid collected in dairy cows1. J Dairy Sci. 2012;95:5978–84.
Golder HM, Denman SE, McSweeney C, Celi P, Lean IJ. Ruminal bacterial community shifts in grain-, sugar-, and histidine-challenged dairy heifers. J Dairy Sci. 2014;97:5131–50.
Jami E, Israel A, Kotser A, Mizrahi I. Exploring the bovine rumen bacterial community from birth to adulthood. ISME J. 2013;7:1069–79.
Li RW, Connor EE, Li C, Baldwin RL, Sparks ME. Characterization of the rumen microbiota of pre-ruminant calves using metagenomic tools. Environ Microbiol. 2012;14:129–39.
Lima FS, Oikonomou G, Lima SF, Bicalho MLS, Ganda EK, de Oliveira Filho JC, et al. Prepartum and postpartum rumen fluid microbiomes: characterization and correlation with production traits in dairy cows. Appl Environ Microbiol. 2015;81:1327–37.
Nieman CC, Steensma KM, Rowntree JE, Beede DK, Utsumi SA. Differential response to stocking rates and feeding by two genotypes of Holstein-Friesian cows in a pasture-based automatic milking system. Animal. 2015;9:2039–49.
Sandri M, Manfrin C, Pallavicini A, Stefanon B. Microbial biodiversity of the liquid fraction of rumen content from lactating cows. Animal. 2014;8:572–9.
Fliegerova K, Tapio I, Bonin A, Mrazek J, Callegari ML, Bani P, et al. Effect of DNA extraction and sample preservation method on rumen bacterial population. Anaerobe. 2014;29:80–4.
McKain N, Genc B, Snelling TJ, Wallace RJ. Differential recovery of bacterial and archaeal 16S rRNA genes from ruminal digesta in response to glycerol as cryoprotectant. J Microbiol Methods. 2013;95:381–3.
Petri RM, Schwaiger T, Penner GB, Beauchemin KA, Forster RJ, McKinnon JJ, et al. Changes in the rumen epimural bacterial diversity of beef cattle as affected by diet and induced ruminal acidosis. Appl Environ Microbiol. 2013;79:3744–55.
Singh KM, Jisha TK, Reddy B, Parmar N, Patel A, Patel AK, et al. Microbial profiles of liquid and solid fraction associated biomaterial in buffalo rumen fed green and dry roughage diets by tagged 16S rRNA gene pyrosequencing. Mol Biol Rep. 2015;42:95–103.
Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci. 2011;108:4516–22.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.
DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with arb. Appl Environ Microbiol. 2006;72:5069–72.
Russell JB. Factors that alter rumen microbial ecology. Science. 2001;292:1119–22.
Henderson G, Cox F, Ganesh S, Jonker A, Young W, Abecia L, et al. Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci Rep. 2015;5:14567.
Brulc JM, Antonopoulos DA, Miller ME, Wilson MK, Yannarell AC, Dinsdale EA, et al. Gene-centric metagenomics of the fiber-adherent bovine rumen microbiome reveals forage specific glycoside hydrolases. Proc Natl Acad Sci U S A. 2009;106:1948–53.
de Menezes AB, Lewis E, O’Donovan M, O’Neill BF, Clipson N, Doyle EM. Microbiome analysis of dairy cows fed pasture or total mixed ration diets. FEMS Microbiol Ecol. 2011;78:256–65.
Thoetkiattikul H, Mhuantong W, Laothanachareon T, Tangphatsornruang S, Pattarajinda V, Eurwilaichitr L, et al. Comparative analysis of microbial profiles in cow rumen fed with different dietary fiber by tagged 16S rRNA gene pyrosequencing. Curr Microbiol. 2013;67:130–7.
Tajima K, Aminov RI, Nagamine T, Matsui H, Nakamura M, Benno Y. Diet-dependent shifts in the bacterial population of the rumen revealed with real-time PCR. Appl Environ Microbiol. 2001;67:2766–74.
Hook SE, Steele MA, Northwood KS, Dijkstra J, France J, Wright AG, et al. Impact of subacute ruminal acidosis (SARA) adaptation and recovery on the density and diversity of bacteria in the rumen of dairy cows. FEMS Microbiol Ecol. 2011;78:275–84.
Fernando SC, Purvis HT, Najar FZ, Sukharnikov LO, Krehbiel CR, Nagaraja TG, et al. Rumen microbial population dynamics during adaptation to a high-grain diet. Appl Environ Microbiol. 2010;76:7482–90.
Zebeli Q, Tafaj M, Weber I, Steingass H, Drochner W. Effects of dietary forage particle size and concentrate level on fermentation profile, in vitro degradation characteristics and concentration of liquid- or solid-associated bacterial mass in the rumen of dairy cows. Anim Feed Sci Technol. 2008;140:307–25.
Castro-Carrera T, Toral PG, Frutos P, McEwan NR, Hervás G, Abecia L, et al. Rumen bacterial community evaluated by 454 pyrosequencing and terminal restriction fragment length polymorphism analyses in dairy sheep fed marine algae. J Dairy Sci. 2014;97:1661–9.
Ramos-Morales E, Arco-Pérez A, Martín-García AI, Yáñez-Ruiz DR, Frutos P, Hervás G. Use of stomach tubing as an alternative to rumen cannulation to study ruminal fermentation and microbiota in sheep and goats. Anim Feed Sci Technol. 2014;198:57–66.
Lodge-Ivey SL, Browne-Silva J, Horvath MB. Technical note: Bacterial diversity and fermentation end products in rumen fluid samples collected via oral lavage or rumen cannula. 2009. 87: 2333-7.
The authors thank L.F. Zhang for analyses of diet nutrition. The receipt of help in samples collection from Y.Q. Wu, J. Mao, and Y. Du, is gratefully acknowledged.
This research was supported by National Dairy Industry and Technology System (CARS-37) and National Natural Science Foundation of China (31402099).
Availability of data and materials
SKJ, SLL, YJW, and ZJC designed the research; SKJ and HTS collected the samples; SKJ and HTZ performed the research, analyzed data and prepared the manuscript; HY, AA and GMA helped improving the data analysis and manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
All authors read and approved the final manuscript for publication.
Experimental protocol was reviewed and approved by Ethical Committee of College of Animal Science and Technology of China Agricultural University.
Microbial diversity in different fractions of rumen content. a, the OTU numbers in original, solid or liquid fraction samples. b, Chao1 index in original, solid or liquid fraction samples. c, Simpson index based on OTUs in original, solid, and liquid fraction samples. HFD: High fiber diet; HED: High energy diet. Data are presented as Mean ± SD. Figure S2. Analysis of similarity (ANOSIM) in different groups. ANOSIM results are presented with box plot when bacteria communities are grouped by diet (a), cows (b), and ruminal content fractions (c) using Bray-Curtis metric based on OTUs. Figure S3. Venn plot for shared OTUs. a, OTUs in HFD and HED. b, OTUs in original, solid and liquid fractions. Figure S4. Ruminal bacteria change in different fractions of rumen content at genera level. LEfSe histogram demonstrating taxonomic differences among different fractions in HFD group (a) and HED group (b) respectively, LDA scores above 2 and P value smaller than 0.05 were shown. LEfSe: linear discriminant analysis (LDA) effect size. Figure S5. Influence of rumen fractions on biomarker taxa abundance. p_: phylum; c_: class; o_: order; f_: family; g_: genus. Data was presented as Mean ± SD. Figure S6. Predominant rumen bacteria at genera level. a, predominant genera higher than 1% in proportion in all samples. b, distribution of predominant genera in each fractions. (DOC 1371 kb)
OTUs distribution in each sample. (XLSX 2056 kb)