- Short report
- Open Access
Effects of live yeast on differential genetic and functional attributes of rumen microbiota in beef cattle
Journal of Animal Science and Biotechnologyvolume 10, Article number: 68 (2019)
Several studies have evaluated the effects of live yeast supplementation on rumen microbial population; however, its effect on differential microbial genes and their functional potential has not been described. Thus, this study applied shotgun metagenomic sequencing to evaluate the effects of live yeast supplementation on genetic and functional potential of the rumen microbiota in beef cattle. Eight rumen-cannulated Holstein steers were randomly assigned to two treatments in a cross-over design with two 25-day experimental periods and a 10-day wash-out between the two periods. The steers were housed in individual pens and fed 50% concentrate-mix and 50% red clover/orchard hay ad libitum. Treatments were (1) control (CON; basal diet without additive) and (2) yeast (YEA; basal diet plus 15 g/d of live yeast product). Rumen fluid samples were collected at 3, 6, and 9 h after feeding on the last d of each period. Sequencing was done on an Illumina HiSeq 2500 platform. Dietary yeast supplementation increased the relative abundance of carbohydrate-fermenting bacteria (such as Ruminococcus albus, R. champanellensis, R. bromii, and R. obeum) and lactate-utilizing bacteria (such as Megasphaera elsdenii, Desulfovibrio desulfuricans, and D. vulgaris). A total of 154 differentially abundant genes (DEGs) were obtained (false discovery rate < 0.01). Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation analysis of the DEGs revealed that 10 pathways, including amino sugar and nucleotide sugar metabolism, oxidative phosphorylation, lipopolysaccharide biosynthesis, pantothenate and coenzyme A biosynthesis, glutathione metabolism, beta-alanine metabolism, polyketide sugar unit biosynthesis, protein export, ribosome, and bacterial secretory system, were enriched in steers fed YEA. Annotation analysis of the DEGs in the carbohydrate-active enzymes (CAZy) database revealed that the abundance of genes coding for enzymes belonging to glycoside hydrolases, glycosyltransferases, and carbohydrate binding modules were enriched in steers fed YEA. These results confirm the effectiveness of a live S. cerevisiae product for improving rumen function in beef steers by increasing the abundance of cellulolytic bacteria, lactic acid-utilizing bacteria, and carbohydrate-active enzymes in the rumen.
Live yeast products containing Saccharomyces cerevisiae are feed additives used to alter rumen fermentation for improved performance, health, and feed efficiency of ruminants . Yeast products favor the growth of fiber-degrading bacteria, stabilize rumen pH by stimulating the population of lactate-consuming bacteria, and increase abundance of rumen microbes by reducing the rumen’s redox potential .
Given the role of live yeast product on rumen function, earlier studies utilized culture-based methods to assess the effects of live yeast in the rumen ; however, this technique is unable to cultivate most naturally occurring species in the rumen . In recent years, several 16S ribosomal ribonucleic acid (rRNA) gene-based sequencing studies have advanced our knowledge of ruminal bacterial diversity and have reported the effects of S. cerevisiae products in the rumen . However, 16S rRNA gene sequencing reveals the bacterial diversity without any insight into the functional potential of microbial communities . In our recent companion study, belonging to the same feeding experiment as the current study , we integrated 16S rRNA gene sequencing and liquid chromatography-mass spectrometry-based metabolomics. Our results revealed that some bacterial genera that responded to treatment with a live S. cerevisiae-based additive showed positive correlations with metabolites involved in amino acid metabolism and biosynthesis as well as metabolism of energy substrates. Therefore, this study applied a shotgun metagenomic sequencing approach to offer taxonomic classification of sequences at the species level and reveal, for the first time, the effects of live yeast supplementation on differentially abundant genes and their functional potential.
All experimental animals were managed according to guidelines approved by the Institutional Animal Care and Use Committee of Kentucky State University (protocol number 18–001). Detailed information about cows and feeding have been reported in our companion paper . Briefly, eight rumen-cannulated Holstein steers were assigned randomly to one of two treatments in a cross-over design with two 25-day experimental periods and a 10-day wash-out between periods. The steers were housed in individual pens and fed 50% concentrate-mix and 50% red clover/orchard hay ad libitum. Mineral mix (Hubbard Feeds; Mankato, MN) was fed free-choice. Dietary treatments were (1) control (CON; basal diet without additive) and (2) yeast (YEA; basal diet plus 15 g/d of Peloton live yeast product; PMI, Arden Hills, MN, USA). The yeast additive was top-dressed on the concentrate mix from d 1 to 25 of each period.
Rumen fluid collection
On d 25 of each experimental period, representative samples (100 mL) of the ruminal contents were collected via the cannula at approximately 3, 6, and 9 h after feeding. At the time of collection, the samples were hand-strained through 4 layers of sterile cheesecloth to separate the liquid and solid samples. Daily composited samples (solid and liquid samples) were mixed 1:1 (w/w) and stored at − 80 °C for subsequent shotgun metagenomic sequencing.
Shotgun metagenomic sequencing
Deoxyribonucleic acid (DNA) extraction and sequencing were performed as reported by Ogunade et al. . Briefly, rumen fluid samples were thawed at room temperature and centrifuged at 15,000×g. DNA was extracted and purified from the pellets using a PowerSoil DNA Isolation Kit (MO BIO Laboratories Inc., Carlsbad, CA). The integrity of the DNA was verified by agarose (0.7%) gel electrophoresis and the DNA was stored at − 20 °C until further use. The sequencing library was constructed by DNA fragmentation . The sheared DNA fragments were then sequenced on an Illumina HiSeq 2500 platform using an Illumina HiSeq – paired-end 150 base pair (bp) strategy. Quality control was performed at each step of the procedure to guarantee the reliability of the sequence data.
Bioinformatics and statistical analysis
The raw sequence data files were demultiplexed and stored as fastq format. Quality filtering, such as removal of low-quality (Q-score ≤ 10) and contaminative reads, was performed. Thereafter, the quality-filtered sequence data were assembled using Iterative De Bruijn Graph De Novo Assembler . The resulting contig sequences (≥ 500 bp) were evaluated using Quality Assessment Tool for Genome Assemblies (QUAST) software . Thereafter, the contig sequences were aligned to National Center for Biotechnology Information’s reference sequence database using Kraken software  to obtain taxonomic composition and relative abundance information. Differences in the relative abundance of taxa at the species level (> 0.1% of total sequences) were analyzed using the generalized linear mixed model procedure in SAS version 9.4 (SAS Institute Inc., Cary, NC, USA); the model included the effects of treatment, period, and their interaction. Differences between means were determined using the Fisher’s test and significant differences were declared at P ≤ 0.05 and tendency was declared at 0.05 < P ≤ 0.10. Normality was tested by examining the distribution of residuals.
The open reading frames in the assembled contigs were predicted using GeneMark software . Redundancy was removed and all predicted genes were clustered using Cluster Database at High Identity with Tolerance (CD-HIT) software (http://www.bioinformatics.org/cd-hit/) with the following parameters: 90% similarity, 90% coverage threshold . Differential gene abundance analysis between CON and YEA was performed using DESeq . Fold change (the ratio of abundance between CON and YEA treatments) ≥ 2 and false discovery rate (FDR) < 0.01, obtained using the Benjamini-Hochberg calibration method, were set as screening criteria to obtain differential gene sets. The differential gene sets were annotated in Kyoto Encyclopedia of Genes and Genomes database (KEGG) pathways to reveal information on differential gene functional potential and the carbohydrate-active enzymes (CAZy) database (http://www.cazy.org/) to obtain information on carbohydrate-active enzymes.
Results and discussion
The modes of action of several yeast products vary due to differences in yeast strain, composition, processing effect, and yeast cell wall constituents. Therefore, comparisons of responses to yeast supplementation among studies require caution and clarification . The product used in this study is a thermal-stable live S. cerevisiae (5.7 × 106 colony forming units/g) with high levels of mannan oligosaccharides and beta-glucans, which were sourced from human-grade food plant (PMI, Arden Hills, MN, USA).
After sequencing of the 16 samples and assembly of the quality-controlled sequence data, 4,333,577 contigs, with an average read length of approximately 1,000 bp, were generated (see Additional file 1). Dietary yeast supplementation increased the relative abundance of 37 microbial species (see Additional file 2). The microbial species (> 1% of the total community) that responded to dietary treatment are shown in Table 1. Ruminococcus albus and R. flavefaciens are cellulose-degrading bacteria , R. bromii degrades resistant starch and xylan, while R. obeum degrades complex deoxy sugars, such as rhamnose and fucose . Several studies have reported the effectiveness of certain strains of S. cerevisiae at increasing the abundance and activities of dominant fiber-degrading rumen microorganisms, such as Fibrobacter and Ruminococcus spp. . This study also revealed that the abundance of relatively unknown and less abundant ruminal cellulolytic microbial species, such as Rhodopseudomonas palustris and Sorangium cellulosum, responded to YEA treatment (see Additional file 2). Some of the microbial species such as R. obeum, Faecalibacterium prausnitzii, and Enterococcus casseliflavus, that responded to YEA treatment are common human gut microbes and are not often detected in the rumen at high abundance (> 1%), therefore, further studies are needed to quantify these bacteria via quantitative polymerase chain reaction.
Yeast additives maintain efficient rumen function necessary for the growth and activities of cellulolytic bacteria by increasing ruminal pH . Yeast additives stimulate the growth of lactate-fermenting bacteria via supply of vitamins and amino acids essential for their growth  and/or outcompeting lactate-producers (such as Streptococcus bovis) for sugar utilization, thereby preventing ruminal lactate accumulation . In this study, the relative abundance of Megasphaera elsdenii, a lactate-fermenting bacterium, and two species of Desulfovibrio (D. desulfuricans and D. vulgaris; see Additional file 2) known to oxidize lactate to acetate , were increased by YEA. It is important to note that lactate concentration was not influenced, probably because the diet fed was non-acidotic (as reported in our companion paper ).
The increased relative abundance of Methanobrevibacter ruminantium, a prominent member of methanogenic archaea, observed in this study is probably due to the increased abundance of fiber-degrading bacteria, especially Ruminococcus albus, a carbohydrate-fermenting, H2-producing organism . This might be linked to H2 transfer between H2-producing cellulolytic bacteria and methanogenic archaea . Fiber is fermented to volatile fatty acids, ammonia, H2, and CO2 in the rumen . Rumen methanogens, such as M. ruminantium, use H2 as a substrate to reduce CO2 to methane . Results of studies on the effects of strains of S. cerevisiae on ruminal methane production have been conflicting. An in vitro study by Chaucheyras-Durand et al.  revealed that a certain live yeast strain increased growth of acetogenic bacteria that can outcompete methanogens by utilizing H2 and CO2 to produce acetate. Newbold and Rode  evaluated the effects of several yeast strains on methane production and reported either no effect or decreased methane production. Our results suggest that YEA treatment might increase the amount of ruminal CH4 produced due to the increased abundance of M. ruminantium; however, we believe that CH4 production per unit of milk or meat will probably be reduced through improved ruminal fiber degradation and feed conversion efficiency.
To our knowledge, this study is the first to evaluate the influence of dietary treatment with yeast product on differentially abundant gene (DEG). A total of 154 DEGs were obtained (fold change > 2; FDR < 0.01). Exactly 139 genes were upregulated in the YEA treatment whereas 15 genes were upregulated in CON. The details of the DEGs, including gene IDs and statistical information, are shown in Additional file 3. None of the genes enriched in CON were annotated in KEGG. KEGG annotation analysis of the upregulated genes in YEA revealed that 10 pathways were enriched (Table 2). These 10 pathways were related to folding, sorting and degradation, translation, membrane transport, and metabolisms of vitamins and co-factors, energy, terpenoids and polyketides, amino acids, and carbohydrate pathways at KEGG orthology level 2 (Table 2). These confirm the results of our companion paper  in which we reported positive correlations between some bacterial genera that responded to yeast treatment and metabolites involved in metabolism of carbohydrates and amino acids. The DEGs related to carbohydrate metabolism were UDP-N-acetyl-2-amino-2-deoxyglucuronate dehydrogenase and β-N-acetylhexosaminidase, both of which are involved in metabolism of amino and nucleotide sugars, which are components of structural polysaccharides. The DEGs related to oxidative phosphorylation were ubiquinol-cytochrome c reductase cytochrome b subunit and cytochrome c oxidase subunit 2, both of which are components of the mitochondrial respiratory chain of S. cerevisiae that catalyzes the reduction of oxygen to water . This indicates an oxygen-scavenging activity, one of the modes of action by which S. cerevisiae improves rumen function. Oxygen is toxic to rumen microbes; therefore, lowering the redox potential in the rumen improves the growing conditions for anaerobic rumen microbes . This probably explains increased relative abundance of strictly anaerobic bacteria, such as Oscillibacter valericigenes and F. prausnitzii, observed in this study.
Lipopolysaccharides are components of gram-negative bacterial cell walls and are required for growth, virulence, and drug resistance of these bacteria . The population of gram-negative bacteria accounts for up to 80%–90% of the total ruminal bacterial population . In this study, increased lipopolysaccharide biosynthesis probably indicates improved ruminal bacterial growth. Feeding live yeast creates favorable conditions for growth and activities of ruminal bacteria through increased pH, lower redox potential, and supply of growth factors .
Rumen microbes are able to synthesize pantothenic acid and other B vitamins ; however, the net vitamin synthesis in the rumen may not be adequate to satisfy the nutritional requirements of the animals, especially in high-producing or growing animals . Saccharomyces cerevisiae is capable of de novo biosynthesis of pantothenic acid and co-enzyme A via the activities of pantoate-beta-alanine ligase (using pantoate and alanine as substrates) and dephospho-CoA kinase (using dephospho-CoA as a substrate), respectively . This explains the enriched pantothenic acid and CoA biosynthesis pathway and increased abundance of Megasphaera elsdenii and Desulfovibrio spp., in YEA treatment because lactate-utilizers require vitamins for their growth and activities . Aminopeptidase N, also called alanyl aminopeptidase, is a hydrolase enzyme that catalyzes the cleavage of alanine from N-terminus of a dipeptide such as alanylglycine . Enrichment of aminopeptidase N is probably due to increased demand of alanine as a substrate for pantothenic acid synthesis by S. cerevisiae.
Enrichment of pathways related to ribosome and protein export probably indicates increased microbial protein synthesis. Ribosomes are organelles that catalyze mRNA-directed protein synthesis in all organisms . In our companion study, we observed lower ammonia-N in steers receiving YEA, which was partly due to its increased uptake by rumen microbes for protein synthesis. In this study, yeast supplementation improved rumen microbial growth, specifically the cellulolytic bacteria. Since cellulolytic bacteria have a high preference for ammonia-N , it was expected that utilization of ammonia-N by the rumen microbes for protein synthesis would increase.
The rumen contains several carbohydrate-active enzymes that are produced by rumen microorganisms and required for hydrolysis of plant cell-wall polysaccharides. A total of 148,586 predicted genes were annotated to the CAZy database, a database containing 6 classes of CAZy: glycoside hydrolases (GHs), glycosyltransferases (GTs), polysaccharide lyases (PLs), carbohydrate binding modules (CBM), carbohydrate esterases (CEs), and auxiliary activities (AAs). The distribution of the annotated genes is shown in Fig. 1. Annotation of the DEGs revealed that 3 genes belonged to GH3, GH16, and GH23 families; 4 genes belonged to GT2, GT4, GT9, and GT83 families; and 1 gene belonged to the CBM37 family (Table 3). Families of enzymes belonging to GHs and CBMs are polysaccharides-degrading enzymes produced by ruminal cellulolytic bacteria . Increased abundance of GHs and CBMs was due to increased growth of cellulolytic bacteria in beef steer fed YEA. Glycosyltransferases represent a class of enzymes involved in synthesizing glycosidic bonds for the biosynthesis of polysaccharides such as lipopolysaccharides . This explains the enriched lipopolysaccharides biosynthesis pathway observed in steers fed YEA.
This study provides new insights into the mode of action of live yeast, such as evidence of oxygen consumption and increased abundance of carbohydrate-active enzymes. Taken together, these results confirm the effectiveness of a live S. cerevisiae product at improving rumen function by increasing the abundance of cellulolytic bacteria, lactic acid-utilizing bacteria, and carbohydrate-active enzymes in the rumen.
Availability of data and materials
All data generated or analyzed are available from the corresponding author on request.
Cluster database at high identity with tolerance
Differentially abundant genes
False discovery rate
Kyoto Encyclopedia of Genes and Genomes
Quality Assessment Tool for Genome Assemblies
Ribosomal ribonucleic acid
Live yeast treatment
Chaucheyras-Durand F, Chevaux E, Martin C, Forano E. Use of yeast probiotics in ruminants: Effects and mechanisms of action on rumen pH, fiber degradation, and microbiota according to the diet. In: Rigobelo E, editor. Probiotic in Animals. Rijeka, Croatia: IntechOpen; 2012. p. 119–52. https://doi.org/10.5772/50192.
Nisbet DJ, Martin SA. The effect of Saccharomyces cerevisiae culture on lactate utilization by the ruminal bacterium Selenomonas ruminantium. J Anim Sci. 1991;69:4628–33.
Ercolini D. High-throughput sequencing and metagenomics: moving forward in the culture-independent analysis of food microbial ecology. Appl Environ Microbiol. 2013;79:3148–55.
Jiang Y, Ogunade IM, Qi S, Hackmann TJ, Staples CR, Adesogan AT. Effects of the dose and viability of Saccharomyces cerevisiae. 1. Diversity of ruminal microbes as analyzed by Illumina MiSeq sequencing and quantitative PCR. J Dairy Sci. 2017;100:325–42.
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.
Ogunade I, Schweickart H, McCoun M, Cannon K, McManus C. Integrating 16S rRNA sequencing and LC–MS-based metabolomics to evaluate the effects of live yeast on rumen function in beef cattle. Animals. 2019;9(1):28.
Ogunade IM, Schweickart H, Andries K, Lay J. Monensin alters the functional and metabolomics profile of rumen microbiota in beef cattle. Animals. 2018;8:211.
Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010;464:59–65.
Peng Y, Leung HCM, Yiu SM, Chin FYL. IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics. 2012;28:1420–8.
Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29(8):1072–5.
Wood DE, Salzberg SL. Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 2014;15:46.
Zhu W, Lomsadze A, Borodovsky M. Ab initio gene identification in metagenomic sequences. Nucleic Acids Res. 2010;38:132.
Li W, Godzik A. CD-HIT: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22:1658–9.
Anders S, Huber W. Differential expression of RNA-Seq data at the gene level–the DESeq package. Heidelberg: European Molecular Biology Laboratory (EMBL); 2012.
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:121–31.
Reichardt N, Duncan SH, Young P, Belenguer A, Leitch CM, Scott KP, et al. Phylogenetic distribution of three pathways for propionate production within the human gut microbiota. ISME J. 2014;8:1323.
Vita N, Valette O, Brasseur G, Lignon S, Denis Y, Ansaldi M, et al. The primary pathway for lactate oxidation in Desulfovibrio vulgaris. Front Microbiol. 2015;6:606.
Iannotti EL, Kafkewitz D, Wolin MJ, Bryant MP. Glucose fermentation products of Ruminococcus albus grown in continuous culture with Vibrio succinogenes: changes caused by interspecies transfer of H2. J Bacteriol. 1973;114:1231–40.
Chassard C, Delmas E, Robert C, Bernalier-Donadille A. The cellulose-degrading microbial community of the human gut varies according to the presence or absence of methanogens. FEMS Microbiol Ecol. 2010;74:205–13.
Loor JJ, Elolimy AA, McCann JC. Dietary impacts on rumen microbiota in beef and dairy production. Anim Front. 2016;6:22–9.
Hook SE, Wright ADG, McBride BW. Methanogens: methane producers of the rumen and mitigation strategies. Archaea. 2010;2010:945785.
Chaucheyras-Durand F, Fonty G, Bertin G, Gouet P. In vitro H2 utilization by a ruminal acetogenic bacterium cultivated alone or in association with an archaea methanogen is stimulated by a probiotic strain of Saccharomyces cerevisiae. Appl Environ Microbiol. 1995;61:3466–7.
Newbold CJ, Rode LM. Dietary additives to control methanogenesis in the rumen. Int Congr Ser. 2006;1293:138–47.
Beattie DS, Jenkins HC, Howton MM. Biochemical evidence for the orientation of cytochrome b in the yeast mitochondrial membrane in the eight-helix model. Arch Biochem Biophys. 1994;312:292–300.
Chaucheyras-Durand F, Walker ND, Bach A. Effects of active dry yeasts on the rumen microbial ecosystem: past, present and future. Anim Feed Sci Technol. 2008;145:5–26.
Kondo KI, Doi H, Adachi H, Nishimura Y. Synergistic effect of CMP/KDO synthase inhibitors with antimicrobial agents on inhibition of production and release of vero toxin by enterohaemorrhagic Escherichia coli O157: H7. Bioorg Med Chem Lett. 2004;14:467–70.
Krause DO, Nagaraja TG, Wright AD, Callaway TR. Board invited review: rumen microbiology: leading the way in microbial ecology. J Anim Sci. 2013;91:331–41.
Volker D, Hüther L, Daş G, Abel H. Pantothenic acid supplementation to support rumen microbes? Arch Anim Nutr. 2011;65:163–73.
NRC. Nutrient Requirements of Dairy Cattle. 7th rev. ed. Washington, DC: National Academies Press; 2001.
White WH, Gunyuzlu PL, Toyn JH. Saccharomyces cerevisiae is capable of de novo pantothenic acid biosynthesis involving a novel pathway of beta-alanine production from spermine. J Biol Chem. 2001;276:10794–800.
Wolin MJ, Miller TL. Microbe-microbe interactions. In: Hobson PN, editor. The rumen microbial ecosystem. London: Elsevier Applied Science; 1988. p. 343–59.
Luan Y, Xu W. The structure and main functions of aminopeptidase N. Curr Med Chem. 2007;14:639–47.
Alberts B, Johnson A, Lewis J, Raff M, Roberts K, Walter P. The Endoplasmic Reticulum. In Molecular Biology of the Cell. 4th ed. New York: Garland Science; 2002.
Bryant MP, Robinson IM. Studies on the nitrogen requirements of some ruminal cellulolytic bacteria. Appl Microbiol. 1961;9:96–103.
Flint HJ, Scott KP, Duncan SH, Louis P, Forano E. Microbial degradation of complex carbohydrates in the gut. Gut Microbes. 2012;3(4):1–18.
Lairson LL, Henrissat B, Davies GJ, Withers SG. Glycosyltransferases: structures, functions, and mechanisms. Annu Rev Biochem. 2008;77:521–55.
The authors are grateful to Hank Schweickart for the daily care and treatment of the animals.
This study (KYSU-000058) was partially supported by funds from Land O′ Lakes Inc. and the US Department of Agriculture’s National Institute of Food and Agriculture Evans-Allen project 1008985.
The research protocol (protocol number 18–001) was reviewed and approved by the Institutional Animal Care and Use Committee of Kentucky State University.
Consent for publication
The authors declare that they have no competing interests.