- 1Department of Animal Science, University of Manitoba, Winnipeg, MB, Canada
- 2Department of Agriculture, Food and Nutrition, University of Alberta, Edmonton, AB, Canada
- 3Diamond V, Cedar Rapids, IA, United States
Background: High-yielding dairy cows are commonly fed high-grain rations. However, this can cause subacute ruminal acidosis (SARA), a metabolic disorder in dairy cows that is usually accompanied by dysbiosis of the rumen microbiome. Postbiotics that contain functional metabolites provide a competitive niche for influential members of the rumen microbiome, may stabilize and promote their populations, and, therefore, may attenuate the adverse effects of SARA.
Methods: This study used a total of 32 rumen-cannulated lactating dairy cows, which were randomly assigned into four treatments: no SCFP (control), 14 g/d Original XPC (SCFPa), 19 g/d NutriTek (SCFPb-1X), and 38 g/d NutriTek (SCFPb-2X) (Diamond V, Cedar Rapids, IA) from 4 weeks before until 12 weeks after parturition. Grain-based SARA challenges were conducted during week 5 (SARA1) and week 8 (SARA2) after parturition by replacing 20% dry matter of the base total mixed ration (TMR) with pellets containing 50% ground barley and 50% ground wheat. The DNA of rumen solids digesta was extracted and subjected to V3-V4 16S rRNA gene sequencing. The characteristics of rumen solids microbiota were compared between non-SARA (Pre-SARA1, week 4; Post-SARA1, week 7; and Post-SARA2, weeks 10 and 12) and SARA stages (SARA1/1, SARA1/2, SARA2/1, SARA2/2), as well as among treatments.
Results: Both SARA challenges reduced the richness and diversity of the microbiota and the relative abundances of the phylum Fibrobacteres. Supplementation with SCFP promoted the growth of several fibrolytic bacteria, including Lachnospiraceae UCG-009, Treponema, unclassified Lachnospiraceae, and unclassified Ruminococcaceae during the SARA challenges. These challenges also reduced the positive interactions and the numbers of hub taxa in the microbiota. The SCFPb treatment increased positive interactions among microbial members of the solids digesta and the number of hub taxa during the SARA and non-SARA stages. The SCFPb-2X treatment prevented changes in the network characteristics, including the number of components, clustering coefficient, modularity, positive edge percentage, and edge density of the microbiota during SARA challenges. These challenges reduced predicted carbohydrate and nitrogen metabolism in microbiota, whereas SCFP supplementation attenuated those reductions.
Conclusions: Supplementation with SCFP, especially the SCFPb-2X attenuated the adverse effects of grain-based SARA on the diversity and predicted functionality of rumen solids microbiota.
1 Introduction
High-grain diets are widely fed to high milk-producing dairy cows to meet their high energy requirements. These diets can result in the accumulation of volatile fatty acids (VFAs) such as acetate, propionate, and butyrate from rumen fermentation and reduced rumen buffering, resulting in rumen pH depressions (Plaizier et al., 2012). When the rumen pH remains for extended periods, e.g., more than 180 min/d, below 5.6, it can result in subacute ruminal acidosis (SARA) (Cooper et al., 1997; Gozho et al., 2005; Plaizier et al., 2008). Further decline of rumen pH below 5.0 results in acute ruminal acidosis, which involves the accumulation of lactate acid in the rumen (Zebeli and Metzler-Zebeli, 2012; Plaizier et al., 2018, 2022). SARA is especially prevalent during early and mid-lactation due to switches from high-forage to high-grain diets (Kleen et al., 2003; Plaizier et al., 2008; Valente et al., 2017). A depressed rumen pH can interrupt the barrier function of the rumen epithelium. This disruption causes the release of immunogenic compounds, such as lipopolysaccharides (LPS), from gram-negative bacteria and their translocation from the rumen into the systemic circulation (Khafipour et al., 2009a; Plaizier et al., 2022). LPS binds to toll-like receptor 4 on cytomembranes and triggers the release of the pro-inflammatory cytokines and the production of acute phase proteins in the liver, causing a systemic immune response in dairy cows (Eckel and Ametaj, 2016; Li et al., 2016).
Rumen bacteria produce a wide array of enzymes to break down structural and non-structural carbohydrates, protein and peptides, and triglycerides into absorbable substrates that can then be utilized by dairy cows (Russell and Rychlik, 2001; Firkins and Yu, 2015; Plaizier et al., 2018; Gruninger et al., 2019). The composition and functionality of the rumen microbiota can be changed because of changes in the diet and/or the availability of substrates for microbiota since some microbes can take advantage of the newly available substrates while others cannot (Khafipour et al., 2009b; Mao et al., 2013; Agler et al., 2016). The colonization of microbes in the rumen relies on other members of the community, such that their competitive advantage can increase when others provide them with metabolites. Consequently, direct interactions among microbes play an important role in determining the rumen microbiome structure (Faust and Raes, 2012, 2016; Agler et al., 2016; Manirajan et al., 2018).
To predict the microbe-microbe interactions from repeated measurements of their presence or abundance, co-occurrence analysis is commonly used as microbial network inference (Faust and Raes, 2016). From microbial interactions, “hub taxa,” which are observed as the most interactive OTUs/ASVs in a microbial network (e.g., having the most positive or negative connections with other members), are key in shaping the microbiome structures. These hub taxa control the abundance of many other microorganisms as they can promote or suppress the growth and diversity of those microbes (Faust and Raes, 2012; Agler et al., 2016; Manirajan et al., 2018; Derakhshani et al., 2020). Therefore, assessing the shifts in microbial connections and hub taxa aids in better understanding the effect of the diet or host phenotype on the entire microbial community. The abundance of several hub taxa, such as fibrolytic bacteria that are pH sensitive, decreases when the rumen pH declines (Nagaraja and Titgemeyer, 2007; Zhou et al., 2015). Accordingly, SARA affects the composition and functionality of the rumen microbiota by reducing its richness, evenness, and the relative abundances of several hub taxa, including those of fibrolytic bacteria (Mao et al., 2013; Brede et al., 2020). Although SARA negatively affects milk production, animal health, and welfare, its clinical symptoms can be easily unnoticed. Hence, it is critical to prevent SARA and attenuate its adverse effects.
Several strategies have been developed to attenuate the negative effect of SARA, such as supplementing diets with buffers, ionophores, prebiotics, probiotics, and postbiotics. Among these, postbiotics have received more attention in recent years. Postbiotics are “a preparation of inanimate microorganisms and/or their components that confers a health benefit on the host” (Salminen et al., 2021). As such, a postbiotic formulation may contain microbial cell components, microbial metabolites produced through the anabolic activity of microbes, and intermediate or end-products of microbial fermentation produced through the catabolic activity of microbes (Duysburgh et al., 2024). Some of these metabolites and compounds affect the cross-feeding patterns among microbes, while others act as signaling molecules or neurotransmitters affecting the microbe-microbe or microbe-host interactions (Salminen et al., 2021). As a result, postbiotics potentially have a wider variety of modes of action compared to other biotic strategies (Salminen et al., 2021). Several studies have shown that supplementation with postbiotics reduces inflammation, promotes milk production, improves the efficiency of feed utilization in lactating dairy cows and calves, and stabilizes rumen pH and rumen fermentation during high-grain feeding (Acharya et al., 2017; Mahmoud et al., 2020; Vailati-Riboni et al., 2021; Guo et al., 2022).
Rumen bacteria comprises a large part (>95%) of the rumen microbiota, which can be divided into different subpopulations, such as rumen liquid microbiota and rumen solids microbiota, based on their adhering and colonization sites (Zhou et al., 2015). Previous studies have mostly assessed the microbiota changes only in rumen liquid fraction due to easier sampling methods, and few have reported shifts in the rumen solids microbial community, especially when cows were exposed to metabolic stressors such as SARA (Martin et al., 2001; Michalet-Doreau et al., 2001; Petri et al., 2013; McCann et al., 2016; Ji et al., 2017; Brede et al., 2020). Additionally, the dynamic changes of the rumen bacteria during early- and mid-lactation when supplemented with postbiotics were rarely examined in previous studies.
In this study, we tested two postbiotics from Saccharomyces cerevisiae fermentation (SCFP) (Original XPC and NutriTek, Diamond V, Cedar Rapids, IA), where NutriTek contained higher concentrations of antioxidants and polyphenol compounds among other metabolites compared to XPC. In our companion papers, we reported that NutriTek was more efficient in stabilizing rumen pH (Khalouei et al., 2020) and rumen liquid microbiota (Guo et al., 2024) and further reduced inflammation in dairy cows during SARA compared to XPC (Guo et al., 2022). Therefore, this study aimed to investigate: (1) the effect of grain-based SARA challenges on microbial diversity and predicted functionalities of the rumen solids microbiota, and (2) the effects of two SCFPs on composition, functionality, and interaction among rumen solids microbes in lactating dairy cows subjected to repeated grain-based SARA challenges.
2 Materials and methods
The protocol used in this study was approved by the University of Manitoba Animal Care Committee (Protocol # F14-038) and followed the guidelines of the Canadian Council for Animal Care (CCAC, 1993).
2.1 Animals, diet, and experimental design
Complete details of the animal management and experimental design are described by Guo et al. (2022). Briefly, a total of 32 rumen-cannulated lactating dairy cows were used in a randomized complete block design with 8 blocks. Cows were placed in blocks based on their parity, previous milk yield, and expected calving date. Cows were fitted with cannulas approximately 12 weeks before calving and had fully recovered before this experiment started, as confirmed by veterinarians. Within each block, cows were randomly assigned to one of four treatments: a basal diet without supplementation (control) or the basal diet supplemented with 14 g/d Diamond V Original XPC (SCFPa, Diamond V), 19 g/d NutriTek (SCFPb-1X, Diamond V), or 38 g/d NutriTek (SCFPb-2X, Diamond V) mixed with 140, 126, 121, and 102 g/d of ground corn, respectively.
The SCFP was supplemented once daily as a top-dress after morning delivery of the diet. A detailed description of the feed sample collection and analyses has been reported by Khalouei et al. (2020). The chemical and nutrient composition of diets are shown in companion papers (Khalouei et al., 2020; Guo et al., 2022). Briefly, cows were fed individually with a prepartum diet containing 38.7% dry matter (DM) of neutral detergent fiber (NDF), 15.5% crude protein (CP), and 17.6% starch (DM basis) from 4 weeks before parturition. They were then switched to a lactation diet containing 34.9% DM NDF, 7.9% DM CP, and 18.6% DM starch until 12 weeks after parturition, except for SARA challenge weeks. On weeks 5 and 8, SARA challenges were conducted (SARA1, SARA2) by gradually replacing 20% DM of the base TMR with pellets containing 50% ground barley and 50% ground wheat over 3 days. The SARA induction diet contained 28.2% DM NDF, 17.2% DM CP, and 27.9% DM starch. The replacement was completed gradually over 3 days before SARA1 and SARA2 weeks. During the 2 weeks between the SARA challenges, cows returned to the base TMR and the treatment supplementations to wash out the first SARA effect. Cows were fed the TMR ad libitum once daily at 0900 h, with feed refusal allowances of 5%−10%. Cows were housed in individual stalls and had free access to fresh water during the whole experiment.
2.2 Sample collection and processing
2.2.1 Sample collection
Rumen digesta samples were taken 6 h after feed delivery once weekly on weeks −4, −1, 1, 3, 4 (Pre-SARA1), 7 (Post-SARA1), 10 (Post-SARA2), and 12 (Post-SARA2) from five sites including cranial, caudal, dorsal, caudal ventral, and caudal dorsal of the rumen. On weeks 5 (SARA 1) and 8 (SARA2), rumen samples were taken on the second and fifth days of these weeks (timepoints were indicated as SARA1/1, SARA1/2, SARA2/1, SARA2/2). After that, rumen solids and liquid digesta were separated using a Bodum coffee filter plunger (Bodum Inc. Triengen, Switzerland) as described in our companion paper (Guo et al., 2022). An approximately 8 g of rumen solids sample was aliquoted into plastic bags, snapped frozen immediately in liquid nitrogen, and stored at −80°C for further microbial analysis.
2.2.2 DNA extraction
Genomic DNA of rumen solids was extracted using Quick-DNA ZR Fecal/Soil DNA kits (D6010; Zymo Research Corp., Orange, CA) following the manufacturer's procedures as described previously (Guo et al., 2024). The procedure included a 2-min bead-beating step at 1750 strokes per minute using a Geno/Grinder equipment (2010, SPEX SamplePrep, Metuchen, NJ, USA). Genomic DNA was quantified using a NanoDrop 2000 spectrophotometer (Thermo Scientific, Waltham, MA, USA) and was normalized to 20 ng/μl followed by MiSeq Illumina sequencing. The DNA samples were quality checked by gel electrophoresis and PCR amplification of the 16S rRNA gene using universal primers 27F (5′-GAAGAGTTTGATCATGGCTCAG-3′) and 342R (5′-CTGCTGCCTCCCGTAG-3′) as described previously (Khafipour et al., 2009b). Amplicons were verified by agarose gel electrophoresis.
2.2.3 PCR amplification and construction of sequencing libraries
These procedures were described by Guo et al. (2024) for rumen liquid digesta samples. The PCR was targeted to amplify the V3-V4 regions of the bacterial 16S rRNA genes using modified F338/R306 primers. The PCR reactions included an initial denaturing step at 94°C for 3 min, followed by 32 amplification cycles at 94°C for 30 s, 55°C for 20 s, and 72°C for 20 s, with a final extension step at 72°C for 5 min in an Eppendorf Mastercycler pro-S (Eppendorf, Hamburg, Germany). Subsequently, the sequencing library was generated and sequenced using a MiSeq Reagent Kit V3 (600-cycle; Illumina, San Diego, CA, USA) at the Gut Microbiome and Large Animal Biosecurity Laboratories, Department of Animal Science, University of Manitoba, Winnipeg, MB, Canada.
2.3 Bioinformatics and statistical analyses
2.3.1 Bioinformatics analyses
Bioinformatics and statistical analysis were performed as described by Guo et al. (2024). Briefly, the DADA2 algorithm was used for quality control, and the feature table was constructed. Sequences were assigned into amplicon sequence variants (ASVs), and taxonomy was classified with QIIME 2 2023.2 (Bolyen et al., 2019). Only samples with sequencing depth >6,072 were kept. Community α-diversity (Shannon's diversity index, Observed Features index, Faith's Phylogenetic Diversity (PD) index, and Pielou's Evenness index) and β-diversity (Jaccard distance, Bray-Curtis distance, unweighted UniFrac distance and weighted UniFrac distance indices) metrics were computed using QIIME 2 default scripts, at an even depth per sample. Non-metric multidimensional scaling (nMDS) was applied on the resulting Bray-Curtis distance matrices to generate two-dimensional plots using default settings of the PRIMER-E software ver. 7.0.17 (Clarke and Gorley, 2015).
2.3.2 Statistical analyses of rumen solids microbiota composition and diversity
Similar to that described by Guo et al. (2024), the univariate procedure of SAS (v 9.4, SAS Institute Inc., Cary, NC, USA) was used for testing the normality of residuals for α-diversity and the relative abundances of bacteria analyses. If the residuals were not normal, the data was log or Box-Cox transformed for normalization before being subjected to the mixed procedure of SAS as described in the results tables. The block was considered a random effect, the stage was repeated measure, and parity and treatment were fixed effects. The effect of parity was removed from the model when its p > 0.10. All pairwise comparisons between the groups were tested using the Tukey studentized range adjustment. Contrast comparisons were made between SARA 1 and SARA 2 stages, as well as control and SCFP, SCFPa and SCFPb, and SCFPb-1X and SCFPb-2X. Significant effects were considered at p < 0.05, and tendencies at 0.05 ≤ p < 0.1 were discussed. Data were presented as means from the original data in this study.
Permutational multivariate analysis of variance (PERMANOVA; implemented in PRIMER-E software v.7.0.17) (Clarke and Gorley, 2015) was used to detect significant differences between β-diversity metrics of bacterial communities as described by Guo et al. (2024). Non-metric multidimensional scaling (nMDS) was applied to the resulting distance matrices to generate two-dimensional plots using the default settings of the PRIMERE. Further to PERMANOVA, permutational multivariate analysis of dispersion (PERMDISP) was performed in PRIMER-E to detect the homogeneity of the dispersions among treatments and stages in rumen solids microbiota (Clarke and Gorley, 2015).
Compositional dynamics of solids digesta bacterial communities of the rumen were assessed using the Metagenomic Longitudinal Differential Abundance (MetaLonDA) method with edgeR package in R (Metwally et al., 2018), similarly as for rumen liquid microbiota (Guo et al., 2024). Taxonomic profiles were normalized with cumulative sum scaling (CSS), and pairwise comparisons were conducted between control vs. SCFPb-2X, control vs. SCFPb-1X, and control vs. SCFPa treatments. Within each comparison, the longitudinal profiles were fitted with a negative binomial smoothing spline. Significant time intervals were identified when p < 0.05 after multiple testing corrections by Benjamini and Hochberg (1995). Data were presented both at the phylum and genus levels.
2.3.3 Functional prediction of rumen solids microbiota
Predicted functions of the rumen solids microbiome, including amino acid, lipid, and carbohydrate metabolic pathways, were assessed by the Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt; CowPI) (Wilkinson et al., 2018). The differences between the control and SCFP supplementation treatment groups during SARA challenges were analyzed by STAMP v2.1.3 (Parks et al., 2014) using an ANOVA test. Significant differences were adjusted using the Tukey-Kramer test and identified when p < 0.05 after multiple testing corrections by Benjamini-Hochberg's false discovery rate (FDR) correction.
2.3.4 Co-occurrence analysis
Correlation network analysis (CoNet) (Derakhshani et al., 2018; Zhang et al., 2018) was used to determine microbial co-occurrence/mutual-exclusion relationships among rumen solids microbial members at the genus level under SARA (SARA1 and SARA2) and non-SARA (pre-SARA1, post-SARA1, post-SARA2) conditions. To determine the influential capacity of bacterial taxa (Trosvik and de Muinck, 2015), their degree of connectedness was calculated by dividing the total number of edges (connections) observed for each phylum by its relative abundance in the community. Hub ASVs in each niche were identified when they had more than 15 edges with other members in the community. The correlation between hub ASVs and biodiversity metrics, rumen VFAs, and ammonia concentrations was conducted by Spearman's rank correlation coefficient (rho) and Corrplot package of R (Wei et al., 2017).
The NetCoMi (network construction and comparison for microbiome data) was used to determine if associations between paired microbes and the overall network structure differed among treatments (Peschel et al., 2021). The same data filtering and normalization were used as in CoNet. Subsequently, group comparisons were conducted by calculating network summary characteristics, including the number of components, clustering coefficients, modularity, positive edge percentage, edge density, and natural connectivity. The network comparison module assessed the overall network structure differences between the treatment groups and experimental periods, using permutation tests with 1,000 replicates, and the differential associations with an alpha of 0.05 were selected (Peschel et al., 2021).
3 Results
Three companion papers reported the effects of SCFPs and SARA on feed intake, milk production, rumen fermentation, and nutrient digestibilities (Khalouei et al., 2020), bacterial endotoxin concentrations and markers of inflammation (Guo et al., 2022), and composition and functionality of microbiota in rumen liquid digesta in lactating dairy cows (Guo et al., 2024). In this study, we detected the effects of SCFPs and SARA on the characteristics of rumen solids microbiota.
3.1 Alpha- and beta-diversity dynamics of rumen solids microbiota
A total of 22,622 ASVs were found in samples with an average of 19,738 ± 6,072 reads per sample through 16S rRNA gene sequencing. The SARA challenges reduced the richness and diversity of microbiota in rumen solids (p < 0.001, Supplementary Table 1, Figure 1). There was no interaction effect of treatment and stage on alpha-diversity indices (p > 0.1). Treatments tended to affect evenness indices of the rumen solids microbiota (p = 0.07). The SCFPa treatment had lower evenness than the control treatment (p = 0.05), whereas no differences were observed between the control and other SCFP treatments across experimental stages. No treatment effect was observed on Shannon's diversity, Faith's PD, or the observed features of the rumen solids microbiota (p > 0.1). The PERMANOVA analysis showed that SCFP (p = 0.0001, Figure 2A) and SARA challenges (p = 0.0001, Figure 2B) affected the β-diversity of rumen solids microbiota. However, the PERMDISP analysis showed a significant treatment effect (p = 0.004, Figure 2A) but not a significant SARA challenge effect (p = 0.07, Figure 2B), confirming homogeneity of dispersions among SARA stages but not treatments. The composition of the rumen solids microbiota differed among treatment groups across the SARA stages.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409659/fmicb-15-1409659-HTML/image_m/fmicb-15-1409659-g001.jpg)
Figure 1. Alpha-diversity dynamics of rumen solids microbiota. Dynamics of (A) Shannon diversity index, (B) Pielou's evenness index, (C) Faith's phylogenetic diversity, and (D) Observed features within treatments (control, SCFPa, SCFPb-1X, and SCFPb-2X) from 4 weeks before until 12 weeks after calving. SARA challenges were conducted on week 5 and week 8 after parturition. Rumen samples were taken weekly but two times during SARA weeks (SARA1/1, SARA1/2, SARA2/1, SARA2/2). Week 4 was considered as pre-SARA1, week 7 as post-SARA1, and weeks 10 and 12 as post-SARA2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409659/fmicb-15-1409659-HTML/image_m/fmicb-15-1409659-g002.jpg)
Figure 2. nMDS of Bray-Curtis distances of rumen solids microbial communities among treatments (A) and stages (B). The ASV table was normalized using cumulative sum scaling (CSS) transformation. Permutational multivariate analysis of variance (PERMANOVA) was used to detect the distinction of clustering patterns between treatments and stages. PERMDISP analysis in Primer 7 was used to detect the effect of the dispersions on the β-diversity of the rumen solids microbiota. p < 0.05 was considered a significant difference. The effect of the block was considered as a random factor in all comparisons. The experimental stage was started 4 weeks before and until 12 weeks after calving. SARA challenges were conducted on week 5 and week 8 after parturition. Rumen samples were taken weekly but two times during SARA weeks (SARA1/1, SARA1/2, SARA2/1, SARA2/2). Week 4 was considered as Pre-SARA1, week 7 as Post-SARA1, and weeks 10 and 12 as Post-SARA2.
3.2 Compositional dynamics of rumen solids microbiota
In total, 14 phyla and 445 genera were identified in the microbiota of rumen solids. The most abundant bacterial phyla in rumen solids were Firmicutes (77–80%), Bacteroidetes (16–19%), Actinobacteria (0.8–1.1%), Proteobacteria (0.5–0.8%), and Fibrobacteres (0.1–0.3%). The average relative abundances of these phyla during the non-SARA and SARA stages in each treatment group are summarized in Table 1. There was a treatment effect on the relative abundance of Fibrobacteres across experimental stages (p < 0.01), and the relative abundance of this phylum was lower in the SCFPb-1X treatment compared with the control and SCFPb-2X treatments (p < 0.01). No treatment effect was observed on the relative abundances of other phyla (p > 0.05). The effect of the SARA stage was significant for all phyla (p < 0.05). There were no interaction effects of treatment and SARA stage with the relative abundance of any identified phyla (p > 0.1). The relative abundances of Actinobacteria were higher in the SARA2 stage than in the post-SARA2 stage (1.15 vs. 0.84%, p < 0.05). During the SARA2 stage, the relative abundance of Proteobacteria was higher than the non-SARA stages (p < 0.05), while SARA1 tended to increase more than post-SARA2 (0.71 vs. 0.46%, p = 0.07). Meanwhile, both grain-based SARA challenges reduced the relative abundance of Fibrobacteres (p < 0.001). Neither the SARA stages nor the SCFP treatment affected the relative abundances of Firmicutes and Bacteroidetes. There was also no treatment or SARA stage effect on the ratio of Firmicutes-to-Bacteroidetes (F:B, p > 0.1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409659/fmicb-15-1409659-HTML/image_m/fmicb-15-1409659-t001.jpg)
Table 1. Effects of treatment (Control, SCFPa, SCFPb-1X, and SCFPb-2X), stage of subacute ruminal acidosis (SARA) induction (Pre-SARA1, SARA1, Post-SARA1, SARA2, and Post-SARA2) on the relative abundances (%) of the main phyla within the rumen solids microbiota in lactating dairy cows.
Comparisons of compositional shifts in proportions of Firmicutes and Bacteroidetes between the control and SCFP treatment groups are shown in Supplementary Figure 1. No significant differences in the fitting splines of Firmicutes and Bacteroidetes were observed between the control and SCFP treatments during the SARA challenges.
At the genus level, taxa that were significantly more abundant within each stage are shown in Figures 3–5. The relative abundances of 33 taxa were increased by either control or SCFPb-2X treatments across the 16-week experimental period (Figure 3). Briefly, the SCFPb-2X treatment increased the relative abundances of several genera, including members of Ruminiclostridium, Lachanospiraceae UCG-009, Candidatus Hepatincola, and Lachnoclostridium during the first SARA challenge (SARA1). The SCFPb-2X treatment also increased members of the family Bacteroidales S24-7 group, Treponema, Lactobacillus, and family Peptococcaceae during SARA2 (p < 0.05). When comparing the control and SCFPb-1X treatments, the relative abundances of 43 bacterial taxa were increased by one of these two treatment groups within each SARA stage (Figure 4). Briefly, the SCFPb-1X treatment increased the relative abundance of Sharpea during the first SARA challenge (SARA1), as well as members of Erysipelotrichaceae UCG-009, Megasphaera, and Lachospiraceae (Syntrophococcus) compared to the control during SARA2 (p < 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409659/fmicb-15-1409659-HTML/image_m/fmicb-15-1409659-g003.jpg)
Figure 3. Differentially abundant taxa in rumen solids microbiota in control vs. SCFPb-2X. The graph summarizes data that were generated using the Metagenomic Longitudinal Differential Abundance (MetaLonDA) method. The significant time intervals were identified when p < 0.05 after multiple testing corrections using the Benjamini-Hochberg False Discovery Rate. X-axes represent the stage relative to parturition. Y-axes represent taxa that were promoted by control (blue lines) or SCFPb-2X groups (red lines) during corresponding stages.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409659/fmicb-15-1409659-HTML/image_m/fmicb-15-1409659-g004.jpg)
Figure 4. Differentially abundant taxa in rumen solids microbiota in control vs. SCFPb-1X. The graph summarizes data that were generated using the Metagenomic Longitudinal Differential Abundance (MetaLonDA) method. The significant time intervals were identified when p < 0.05 after multiple testing corrections using the Benjamini-Hochberg False Discovery Rate. X-axes represent the stage relative to parturition. Y-axes represent taxa that were promoted by control (blue lines) or SCFPb-1X groups (red lines) during corresponding stages.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409659/fmicb-15-1409659-HTML/image_m/fmicb-15-1409659-g005.jpg)
Figure 5. Differentially abundant taxa in rumen solids microbiota in control vs. SCFPa. The graph summarizes data that were generated using the Metagenomic Longitudinal Differential Abundance (MetaLonDA) method. The significant time intervals were identified when p < 0.05 after multiple testing corrections using the Benjamini-Hochberg False Discovery Rate. X-axes represent the stage relative to parturition. Y-axes represent taxa that were promoted by control (blue lines) or SCFPa groups (red lines) during corresponding stages.
A comparison of the control and SCFPa treatments revealed that the relative abundances of 39 bacterial taxa were increased by one of these two treatments (Figure 5). Briefly, the SCFPa treatment increased the relative abundances of Prevotella, Butyrivibrio, and Anaerovibrio during the SARA1 stage. It increased these abundances of members of Ruminococcaceae, Succinivibrionaceae, Syntrophococcus, Olsenella, and Lactobacillus during the SARA2 stage, as well as Corynebacterium, [Eubacterium] ruminantium and Streptococcus during the first week after calving (p < 0.05).
3.3 Predicted functionality of rumen solids microbiota
The predicted functionality of rumen solids microbiota revealed 254 endogenous third-level KEGG pathways using CowPi. Of these, 150 were considered rumen microbial metabolic pathways. The comparison between non-SARA and SARA stages in each treatment group is shown in Supplementary Figures 2–5. Five pathways related to carbohydrate metabolism, including butanoate metabolism, propanoate metabolism, pyruvate metabolism, carbohydrate digestion and absorption, and starch and sucrose metabolism, were inhibited by the SARA challenges in all treatment groups. These challenges inhibited four amino acid metabolic pathways, including alanine, aspartate, and glutamate metabolism; lysine biosynthesis, cysteine, and methionine metabolism; and valine, leucine, and isoleucine biosynthesis in all treatment groups. However, these challenges only inhibited histidine metabolism in the control and SCFPa treatment groups. Additionally, SARA promoted the synthesis and degradation of ketone bodies in the control and SCFPb-1X treatment groups but not in the SCFPa and SCFPb-2X groups.
The comparison between the control and SCFP groups during the SARA challenges is shown in Figure 6. Compared with control, SCFP supplementation promoted three carbohydrate metabolic pathways, including ascorbate and aldarate metabolism, fructose and mannose metabolism, and inositol phosphate metabolism. It inhibited two carbohydrate metabolic pathways, including starch and sucrose metabolism and citrate cycle (TCA cycle) during SARA challenges. Furthermore, SCFP supplementation inhibited three amino acid metabolic pathways, including valine, leucine, and isoleucine biosynthesis; histidine metabolism; and phenylalanine metabolism. It promoted one amino acid metabolic pathway (amino acid metabolism) compared with the control treatment during SARA challenges. In addition, SCFP supplementation promoted three lipid metabolic pathways, including primary bile acid biosynthesis, secondary bile acid biosynthesis, and synthesis and degradation of ketone bodies, but inhibited one lipid metabolic pathway (arachidonic acid metabolism) compared with the control treatment during the SARA challenges. A comparison of predicted microbial functionalities between the control treatment and each SCFP treatment group during SARA showed that three carbohydrate pathways (including ascorbate and aldarate metabolism, fructose and mannose metabolism, and inositol phosphate metabolism), one amino acid metabolism, and three lipid metabolic pathways (including primary bile acid biosynthesis, secondary bile acid biosynthesis, and synthesis and degradation of ketone bodies) were promoted by the SCFPa treatment (p < 0.05, Supplementary Figure 6). However, no differences in these predicted functionalities were observed between the control and SCFPb treatment groups.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409659/fmicb-15-1409659-HTML/image_m/fmicb-15-1409659-g006.jpg)
Figure 6. Differences in predicted functions of rumen solids microbiota between control and SCFP groups during SARA. Functionalities of rumen solids microbiota were predicted using CowPi. Output was analyzed using STAMP following log_transformation and False Discovery Rate correction. Significant differences were considered as p < 0.05.
3.4 Effect of SCFP supplementation and SARA inductions on co-occurrence patterns of rumen solids microbiota
The SARA challenges reduced the total number of significant associations among bacterial taxa in rumen solids microbiota (Figure 7A). No negative connections were detected in the SCFPb treatment group during the SARA challenges (Figure 7B). The relative degree of connectedness of each phylum (total number of positive and negative edges observed for each phylum divided by its relative abundance in the community) varied among treatment groups during the SARA stages.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409659/fmicb-15-1409659-HTML/image_m/fmicb-15-1409659-g007.jpg)
Figure 7. Microbial interaction networks. The degree of connections for each phylum was normalized by dividing the total number of positive and negative edges observed for each phylum by their relative abundance in the community. (A) Normalized positive connections for the dominant bacteria phyla within each treatment group during non-SARA (Pre-SARA1, Post-SARA1, and Post-SARA2) and SARA (SARA1/1, SARA1/2, SARA2/1, SARA2/2) stages. (B) Normalized negative connections for the dominant bacteria phyla with each treatment group during non-SARA and SARA stages.
During non-SARA stages, the degrees of positive connections were higher in SCFP treatment groups compared to the control treatment group, and these connections were higher in the SCFPb treatment group than in the SCFPa treatment group (Figure 7A). Bacteroidetes had the highest degree of positive connections in the SCFP treatment groups. In descending order, Bacteroidetes, Actinobacteria, and Proteobacteria showed high degrees of positive connections in the SCFPb treatment groups, whereas Firmicutes and Fibrobacteres showed high degrees of positive connections in the SCFPa treatment group. Meanwhile, Firmicutes had the highest degree of connections in the control treatment group.
Across SARA challenges, the SCFP treatment groups had higher positive connections than the control treatment group, and these connections were higher in the SCFPb treatment group than in the SCFPa treatment group (Figure 7A). Actinobacteria contributed the most to the positive connections in SCFPb treatment groups, followed by Bacteroidetes and Firmicutes. In the SCFPa treatment group, Bacteroidetes had the highest number of positive connections. Bacteroidetes and Firmicutes contributed equally to positive connections in the control treatment group (Figure 7A). Actinobacteria, Bacteroidetes, and Firmicutes had equal numbers of negative connections in the control group, whereas in the SCFPa group, Firmicutes was the major contributor to the negative connections, followed by Bacteroidetes and Actinobacteria (Figure 7B).
Hub taxa were identified as taxa with a high number (>15) of positive or negative connections with other members of the community (Figure 8). The SARA challenges reduced the number of hub taxa in the control and SCFPb treatment groups but not in the SCFPa treatment group. In the control treatment group, all hub taxa were from the Firmicutes phylum during both SARA and non-SARA stages. During non-SARA stages, cows in the control group had nine negatively connected hub taxa (three from the Lachnospiraceae NK3A20 group, one from the family Lachnospiraceae, four from the Ruminococcaceae NK4A214 group, and one from Ruminococcus) and one positively connected hub taxon (from genera Oribacterium) (Figure 8A). During the SARA stages, control cows had two negatively connected hub taxa (one from Oribacterium and one from Mitsuokella) and one positively connected taxon from the family Lachnospiraceae (Figure 8B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409659/fmicb-15-1409659-HTML/image_m/fmicb-15-1409659-g008.jpg)
Figure 8. Distribution of hub taxa within each treatment group during non-SARA and SARA stages. Hub taxa were identified as OTUs with more than 15 connections with other members of rumen solids microbiota within each treatment group during non-SARA (Pre-SARA1, Post-SARA1, and Post-SARA2) and SARA (SARA1/1, SARA1/2, SARA2/1, SARA2/2) stages. (A, B) control; (C, D) SCFPb-2X; (E, F) SCFPb-1X; (G, H) SCFPa. Orange color represents negative connections while blue color represents positive connections.
In the SCFPb-2X treatment group, all positively connected hub taxa were from the Bacteroidetes and Firmicutes phylum during non-SARA stages and from the Bacteroidetes, Firmicutes, and Actinobacteria phyla during SARA stages. During non-SARA stages, SCFPb-2X cows had positively connected hub taxa from Prevotellaceae (one from Prevotellaceae UCG-003, five from Prevotella), Christensenellaceae (four from genera Christensenellaceae R-7 group), Lachnospiraceae (one from Lachnospiraceae NK3A20 group, five from family Lachnospiraceae, one from Lachnospiraceae UCG-006), Ruminococcaceae (two from Ruminococcaceae NK4A214 group, two from family Ruminococcaceae and two from Saccharofermentans) and one from genera Rikenellaceae RC9 gut group (Figure 8C). During SARA stages, SCFPb-2X cows had positively connected hub taxa from Lachnospiraceae (one from genera Lachnospiraceae NK3A20 group and two from family Lachnospiraceae), Ruminococcaceae (one from genera Anaerotruncus), one from family Prevotellaceae, family Bacteroidales Incertae Sedis (one from genara Phocaeicola) and family Coriobacteriaceae (one from genera Senegalimassilia) (Figure 8D).
Similarly, all positively connected hub taxa in the SCFPb-1X treatment group were from the Bacteroidetes and Firmicutes phyla during non-SARA stags and from the Firmicutes phylum during SARA stages. During the SARA stages, hub taxa were from Ruminococcaceae (one from the genera Anaerotruncus, three from the genera Saccharofermentans, one from family Ruminococcaceae, and one from Ruminococcaceae UCG-014), Christensenellaceae (five from Christensenellaceae R-7 group), Lachnospiraceae (one from [Eubacterium] xylanophilum group, one from Lachnospiraceae XPB1014 group, and three from family Lachnospiraceae), Erysipelotrichaceae (one from Erysipelotrichaceae UCG-004, one from Erysipelotrichaceae UCG-007, one from Solobaterium, one from [Anaerorhabdus] furcosa group), Bacillaceae (one from Bacillus), Eubacteriaceae (one from Anaerofustis), Rikenellaceae (three from Rikenellaceae RC9 gut group), Prevotellaceae (one from Prevotella, one from Prevotellaceae UCG-003, one from family Prevotellaceae) and one from Family XIII AD3011 group (Figure 8E). During SARA stages, hub taxa were from Ruminococcaceae (one from Caproiciproducens, one from Saccharofermentans, two from family Ruminococcaceae) and two from family Lachnospiraceae (Figure 8F).
In the SCFPa treatment group, hub taxa that were negatively connected were from Bacteroidetes and Firmicutes during non-SARA stages and had a combination of positively and negatively connected hub taxa from Firmicutes during SARA stages. During non-SARA stages, SCFPa cows had hub taxa from Lachnospiraceae (one from Lachnospiraceae XPB1014 group, one from [Ruminococcus] gauvreauii group, one from Lachnospiraceae NK3A20 group and two from family Lachnospiraceae), Ruminococcaceae (one from Ruminococcus, one from Ruminococcaceae NK4A214 group, one from Saccharofermentans and one from family Ruminococcaceae), Erysipelotrichaceae (one from Solobacterium and one from Sharpea), Christensenellaceae (one from Christensenellaceae R-7 group), and Prevotellaceae (one from Prevotella) (Figure 8G). During SARA challenges, cows had six positively connected hub taxa (one from Lachnospiraceae NK3A20 group, four from genera Succiniclasticum, and one from Ruminococcaceae NK4A214 group) and seven negatively connected taxa (three from [Ruminococcus] gauvreauii group, one from family Lachnospiraceae, one from family Bacteroidales S24-7 group, one from Erysipelotrichaceae UCG-007, and one from Ruminococcus) (Figure 8H).
As shown in Table 2, results from the NetCoMi analysis indicated that network characteristics, including modularity (p = 0.04), clustering coefficient (p = 0.09), positive edge percentage (p = 0.07), and edge density were increased during SARA stages compared to non-SARA stages in the control group. In the SCFPb-1X treatment group, the clustering coefficient (p = 0.02), positive edge percentage (p < 0.001), edge density (p < 0.001), and natural connectivity (p < 0.001) were higher during SARA stages compared to non-SARA stages. In the SCFPa treatment group, the clustering coefficient (p = 0.02), positive edge percentage (p < 0.001), edge density (p < 0.001), and natural connectivity (p < 0.001) were also increased, while modularity decreased (p = 0.04) during SARA stages compared to non-SARA stages. However, there was no difference in network characteristics between non-SARA and SARA stages in the SCFPb-2X treatment group.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409659/fmicb-15-1409659-HTML/image_m/fmicb-15-1409659-t002.jpg)
Table 2. Comparison of network properties in each treatment group (control, SCFPa, SCFPb-1X, and SCFPb-2X) between non-SARA and SARA stages.
3.5 Correlation between hub taxa, biodiversity metrics of rumen solids microbiota, and rumen fermentation characteristics
The association between hub taxa and biodiversity metrics of rumen solids microbiota, rumen VFA, and ammonia concentrations that were published in a companion paper (Khalouei et al., 2020) are shown in Figure 9. A group of taxa from Bacteroidetes, including Rikenellaceae RC9 gut group, Firmicutes, including Christensenellaceae R-7 group, Family XIII AD3011 group, Lachnospiraceae XPB1014 group, unclassified family Lachnospiraceae, Anaerotruncus, Ruminococcaceae NK4A214 group, unclassified family Ruminococcaceae, [Anaerorhabdus] furcosa group, and genera Saccharofermentans were negatively correlated with rumen concentrations of propionate, butyric, valerate, and lactate (p < 0.05). However, they were positively correlated with rumen concentrations of acetate and ammonia and with α-, β-diversity metrics of the microbiota in solid rumen digesta (p < 0.05). A group of taxa from Firmicutes that included members of the Ruminococcus gauvreauii and Sharpea genera were positively correlated with rumen propionate and butyrate concentrations but were negatively correlated with rumen acetate and ammonia concentrations and with α-, β-diversity metrics of the rumen solids microbiota (p < 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409659/fmicb-15-1409659-HTML/image_m/fmicb-15-1409659-g009.jpg)
Figure 9. (A, B) Correlation among hub taxa and rumen fermentation characteristics and biodiversity metrics of rumen solids microbiota. Spearman's correlation coefficient was used to explore the relationships between the relative abundances of rumen solids hub taxa and community α-diversity (Shannon, Evenness, Faith_PD, and Observed_Features), β-diversity (Bray-Curtis dissimilarities) and rumen fermentation characteristics (VFAs and ammonia concentrations). *indicates p < 0.05. The color ramp and the size of the squares indicate the type and strength of Spearman's correlation coefficient (rho): rho = 1 shows a strong positive correlation, and rho = −1 shows a strong negative correlation between the two parameters.
4 Discussion
Postbiotics contain bioactive compounds and functional metabolites with a dual mode of action that on the one hand can prime the immune response (Vailati-Riboni et al., 2021; Guo et al., 2024) and, on the other hand can improve rumen health and fermentation by supporting the populations of hub taxa, such as fibrolytic bacteria, therefore increasing the stability and robustness of microbial community during dietary, metabolic or infectious stressors (Zhu et al., 2017; Tun et al., 2020; Ganda et al., 2023). Bacteria are the most abundant kingdom in the rumen microbial community, and rumen solids-associated bacteria have been estimated to comprise more than 70% of bacterial biomass in the rumen (Craig et al., 1987; Mullins et al., 2013). These bacteria play key roles in the degradation and fermentation of the digesta (Bickhart and Weimer, 2018). The relationship between the changes in the diet and the composition of the microbial community in the rumen liquid has been widely studied. Our study revealed the effects of two commercial postbiotics from Saccharomyces cerevisiae fermentation on diversity, composition, predicted functionality, and network structure of the bacterial community in rumen solids in dairy cows subjected to repeated grain-based SARA challenges.
4.1 Effects of grain-based SARA challenges on diversity, composition, and predicted functionality of rumen solids microbiota
Companion papers to this work showed that SARA was induced successfully, as the SARA challenges increased the duration of rumen pH below 5.6 from 8 to 186 min/d (Khalouei et al., 2020) and the concentration of rumen free-LPS from 5,012 to 63,596 endotoxin unit (EU)/mL (Guo et al., 2022), both of which are indicative of SARA (Plaizier et al., 2008, 2012). Previous studies have shown that grain-based SARA challenges reduce the rumen pH and the relative abundances of pH-sensitive microorganisms, such as fibrolytic bacteria, but increase the relative abundance of low pH tolerant and amylolytic bacteria (Cherdthong et al., 2010; Mao et al., 2013; Petri et al., 2013; Plaizier et al., 2017). These changes in fiber- and starch-degrading bacteria during SARA may further influence the rumen pH because of the lower fiber and higher starch fermentation. In our study, SARA was induced by increasing the starch content of the diet from 17.6% to 27.9% DM while decreasing the dietary NDF from 55.4% to 48.1%. Hence, we expected similar changes in the populations of rumen bacteria as in previous work.
Earlier studies reported that SARA reduces the richness and diversity of the rumen liquid microbiota (Fernando et al., 2010; Mao et al., 2013; McCann et al., 2016; Guo et al., 2024). The impact of SARA on the diversity of rumen solids microbiota has been less studied. McCann et al. (2016) reported a surprise increase in the diversity of rumen solids following 6 days of SARA induction, whereas Brede et al. (2020) observed a decrease in the diversity indices when SARA was induced in an in vitro RUSITEC system. In the current study, we observed a reduction in the richness and diversity of the microbiota in the solid fraction during grain-based SARA challenges. This implies that induction of SARA was able to disturb the rumen solids microbial community regardless of treatments. The PERMANOVA and PERMDISP analyses of Bray-Curtis distances also demonstrated that SARA challenges had major effects on the overall composition and beta diversity of the rumen solids microbiota.
In total, 14 phyla and 445 genera were identified in the rumen solids microbiota. The predominant phyla were Firmicutes and Bacteroidetes, which comprised more than 90% of the community. The relative abundance of Firmicutes was higher than that of Bacteroidetes in the solid fraction. It has been recognized that microbiota of rumen solids plays a major role in fiber degradation, as cellulolytic bacteria are more abundant, resulting in a greater prevalence of carbohydrate active enzymes (CAZy) in this fraction compared to the liquid fraction (Williams et al., 1989; Michalet-Doreau et al., 2001). Earlier studies have also suggested that perhaps the abundance of bacteria in the rumen solids fraction is greater than in the liquid fraction (>60% vs. < 20% of total bacteria; Yang et al., 2001); however, this claim needs to be re-examined using more accurate methods. Martin et al. (2001) reported that supplementation with barley grains decreased the rumen degradation rate of hay and the fibrolytic activities of bacteria in the rumen solids (Martin et al., 2001). Increased grain feeding is expected to increase the proportion of Firmicutes and decrease that of Bacteroidetes, at least in the rumen liquid environment (Khafipour et al., 2009b; Huo et al., 2014). McCann et al. (2016) reported that the relative abundance of Firmicutes in the rumen solids fraction decreased, and that of Bacteroidetes increased following 6 days of SARA induction. Similarly, Fernando et al. (2010) observed a gradual increase in the proportion of Bacteroidetes in the rumen when beef steers were adapted to a high-concentrate diet. In contrast with previous studies, we did not observe significant changes in relative abundances of Firmicutes and Bacteroidetes in rumen solids during grain-based SARA. This discrepancy may be due to different sampling strategies, experimental designs, and microbiome evaluation methods. We separated rumen solids and assessed the microbiota of this fraction, whereas the study of Fernando et al. (2010) evaluated the whole rumen content. Further, we used a complete randomized block design to avoid any crossover effects, whereas the above-mentioned studies incorporated a Latin square design.
The Lachnospiraceae family belongs to the phylum Firmicutes and appears to be more abundant in rumen solids fraction than in liquid fraction in cows fed in pastures (De Menezes et al., 2011). Several members of the Lachnospiraceae family have been reported to be major butyrate producers (Louis et al., 2010; Meehan and Beiko, 2014). Studies have been conducted on the correlation between the abundance of this family and feed efficiency, but the results are inconsistent (Guan et al., 2008; Myer et al., 2015; Li and Guan, 2017). These discrepancies among studies could be due to limited butyrate production by some members of this family (De Menezes et al., 2011). In our companion paper, SARA challenges increased the concentration of butyrate in the rumen liquid (Khalouei et al., 2020). Also, we observed that the abundance of members of the Lachnospiraceae family either positively or negatively correlated with the butyrate concentration in the rumen liquid environment (Guo et al., 2024). However, we found that SARA challenges did not affect the Firmicutes-to-Bacteroidetes ratio in the rumen solids or the liquid fraction.
Fibrobacteres, which are an important phylum of the cellulolytic bacteria, were reduced during SARA in this study. Consistently, Lourenco et al. (2020) reported that supplementation with grain decreases the abundance of Fibrobacteres in the rumen liquid digesta of beef calves. Fernando et al. (2010) detected a higher proportion of Fibrobacteres in rumen solids and liquid digesta in hay-fed animals than in grain-fed animals. The reduction of the relative abundance of this phylum during SARA stages was expected as the dietary starch content was higher, and that of fiber was lower in the SARA diet compared with those during non-SARA stages. Actinobacteria, a gram-positive phylum, also increased in rumen solids during SARA in our study. In this phylum, the relative abundance of the most abundant genus Bifidobacterium from the family Bifidobacteriaceae, a lactic acid producer, has been shown to increase during SARA challenges (Nagaraja and Titgemeyer, 2007; Plaizier et al., 2017; Monteiro and Faciola, 2020). Previous studies have also reported that the phylum Proteobacteria is more abundant in the rumen when animals are fed a high-grain diet (Fernando et al., 2010; Petri et al., 2013; Auffret et al., 2017). The relative abundance of Proteobacteria increased during SARA2 in our study, which is consistent with previous observations. The Proteobacteria phylum contains many opportunistic and pathogenic bacteria, such as Escherichia coli. The relative abundances of several of these pathogens can be enriched when feeding high-grain diets (Diez-Gonzalez et al., 1998; Khafipour et al., 2011; Bäumler and Sperandio, 2016; Auffret et al., 2017). Consequently, the increase in the abundance of Proteobacteria could be indicative of the dysbiosis of the rumen microbiota and that greater prevalence of opportunistic and pathogenic microorganisms may increase the risk of inflammation and metabolic dysfunction in the host (Khafipour et al., 2011; Petri et al., 2013; Shin et al., 2015).
Due to the changes in the composition of rumen solids microbiota during SARA challenges, differences in microbiome functionalities were expected. It has been reported that feeding highly fermentable carbohydrates promotes the production of propionate, butyrate, and lactate and reduces the ratio of acetate-to-propionate in the rumen (Gozho et al., 2007; Li et al., 2016; Khalouei et al., 2020). Liu et al. (2022) found that cattle fed with a grain diet had higher relative abundances of the starch-fermenting bacteria, such as Succinivibrionaceae and Succinimonas, and higher relative abundances of the lactate-utilizing bacteria, such as Megasphaera and Acetobacter. Our study found that the relative abundance of genera Succiniclasticum was positively correlated with the propionate concentration in the rumen solids digesta. As a major succinate producer, the family Succinivibrionaceae competes with methanogens for hydrogen (Pope et al., 2011). Succinate is the precursor for propionate, which is a substrate for gluconeogenesis (Yost et al., 1977; Li and Guan, 2017). Hence, a higher ratio of propionate-to-acetate may reduce methane emissions from the rumen (Russell, 1998). In our companion paper, a higher propionate-to-acetate ratio during grain-based SARA challenges was observed (Khalouei et al., 2020). Gagen et al. (2015) found acetogens in both families of Lachnospiraceae and Ruminococcaceae, which can provide hydrogen, and their abundance may increase while methane production declines. It has been reported that increases in rumen acetate can increase milk fat production as acetate is the major source of energy and substrate for milk fat synthesis in dairy cows (Urrutia and Harvatine, 2017; Urrutia et al., 2019). We determined that the relative abundances of taxa from the Lachnospiraceae NK3A20 group, Ruminococcaceae NK4A214 group, and Ruminococcus and Ruminococcaceae groups were positively related to the acetate concentration in rumen liquid and with the diversity and richness of the microbial community in the rumen solids. Therefore, the growth of those bacteria was inhibited during grain-based SARA.
Several studies have observed a positive relationship between less diversity in the rumen microbial community and feed efficiency (McCann et al., 2014; Myer et al., 2015; Li and Guan, 2017). Li and Guan (2017) showed that nitrogen metabolism activities were inhibited in efficient cows. Similarly, we found predicted nitrogen metabolism pathways such as alanine, aspartate, and glutamate metabolism; lysine biosynthesis; phenylalanine metabolism; glycine, serine, and threonine metabolism; and valine, leucine, and isoleucine biosynthesis were inhibited during SARA. Valine, leucine, and isoleucine are important contributors to microbial protein synthesis, and microbial proteins are precursors for the synthesis of milk protein in the mammary gland (Allison et al., 1966; Xue et al., 2020). Thus, the challenges of grain-based SARA in our study may reduce milk protein synthesis. That being said, our predicted functionalities of the rumen microbiome can be biased as they are based on amplicon sequencing, which can only distinguish bacteria at the genus level. Different strains from the same species can perform different functions, while several species can contribute to similar or different functions (Benson et al., 2010; Qi et al., 2024). To properly assess the functional shifts in the rumen solids microbiota during SARA induction, more comprehensive omics studies such as shotgun metagenomics and metabolomics need to be performed.
4.2 Effects of SCFP supplementation on diversity, composition, and predicted functionality of rumen solids microbiota
The SCFP are commonly used as rumen fermentation modifiers. It has been reported that SCFP influences rumen fermentation by promoting the growth and functionality of fibrolytic bacteria and lactate-utilizing bacteria, as well as increasing microbial protein synthesis and reducing lactate accumulation (Callaway and Martin, 1997; Zhu et al., 2017). Different from our companion study, where we showed SCFP supplementation attenuates the fluctuations in the abundances of main phyla in the rumen liquid microbiota (Guo et al., 2024), such effects were not observed in the rumen solids microbiota during SARA challenges. We speculate this could be due to differences in relative abundances of phyla between rumen liquid and solid fractions that made them respond differently to SCFP supplementation (Petri et al., 2013; Brede et al., 2020).
It has been reported that SCFP supplementation increases rumen VFA production, resulting in greater milk yield and feed efficiency in lactating dairy cows (Callaway and Martin, 1997; Hristov et al., 2010; Poppy et al., 2012). Xiao et al. (2016, 2018) found that the SCFP supplementation (XPC) increased ruminal butyrate concentrations and Butyrivibrio abundance, a genera within the Lachnospriraceae family. The authors also reported increased richness in rumen liquid microbiota but did not observe any effect on the richness and diversity of rumen solids microbiota in calves. In the current study, we also observed that SCFPa supplementation promoted the growth of Butyrivibrio in rumen solids during SARA challenges; however, as reported before, SCFP supplementation did not increase butyrate concentration (Khalouei et al., 2020). Butyrate and propionate provide energy for the rumen epithelial cells; thus, higher concentrations of these VFAs may promote the development of rumen papillae and the absorption area for VFA (Meehan and Beiko, 2014). However, excessive accumulation of these VFA in the rumen will reduce the rumen pH (Gozho et al., 2005). It has been demonstrated that SCFP promotes the length and weight of rumen papillae (Lesmeister et al., 2004; Brewer et al., 2014; Xiao et al., 2016) and stabilizes the rumen pH during high-starch diet feeding or under grain-based SARA conductions (Allen and Ying, 2012; Li et al., 2016; Dias et al., 2018). In agreement with the above, our companion papers (Khalouei et al., 2020; Guo et al., 2022) reported that SCFPb-2X supplementation increased the rumen pH, reduced ruminal propionate and free-LPS concentrations during grain-based SARA challenges, indicating the benefit of SCFP on improving rumen health and reducing the adverse effects of SARA.
Different members of the rumen microbiota utilize substrates different from those in the diet. Thus, maintaining a diverse microbial community during periods of dietary or metabolic stress can improve the efficiency of nutrient utilization (Henderson et al., 2015; Khafipour et al., 2016). Tun et al. (2020) found that SCFP supplementation attenuated the reduction in richness and diversity of the rumen liquid microbiota during the SARA challenge. Similarly, our parallel study showed that the reduction in the evenness of the rumen liquid microbiota following exposure to SARA was attenuated by SCFP supplementation (Guo et al., 2024). However, we did not observe a treatment effect on the alpha diversity of rumen solid microbiota in the current study. We did, however, observe a treatment effect on the beta diversity as determined by PERMANOVA. That being said, the significant treatment effect obtained by PERMDISP analysis, which takes into account dispersion in data, suggests that the observed effect may have been influenced by the differences in the taxonomic composition within each group.
The effect of SCFP supplementation on the functionality of the rumen microbiota was also studied. Nitrogen metabolism in the rumen includes protein degradation that provides amino acids and nitrogen for bacteria and microbial protein synthesis that is the major source of protein for the cow when it reaches the small intestine (Storm and Ørskov, 1983; Clark et al., 1992). Protein degradation in the rumen requires the combination of proteolytic and nonproteolytic enzymes and amino acid uptake by rumen microbes (Bach et al., 2005). Positive effects of amylases on protein degradation have been reported in earlier studies (Assoumani et al., 1992; Tománková and Kopečný, 1995). Bach et al. (2005) reported that the abundance of cellulolytic bacteria decreases while proteolytic bacteria are less affected as rumen pH declines. When fed high-concentrate diets, starch-degrading bacteria predominate in the rumen while the number of cellulolytic bacteria decreases, and as a consequence, fiber and protein degradation reduces (Mould and Ørskov, 1983; Cherdthong et al., 2010). Carbohydrates are the main energy resources for rumen bacteria, and they can also be used as carbon skeletons for microbial protein synthesis. The nitrogen uptake by the ruminal microbes can be improved when the readily fermentable carbohydrates increase (Stern and Hoover, 1979; Casper and Schingoethe, 1989; Cameron et al., 1991). When the rate of carbohydrate fermentation exceeds the protein degradation rate, microbial protein synthesis decreases (Nocek and Russell, 1988). The reduction of ruminal nitrogen metabolism when cows are fed a high-grain diet may be related to the increase in digesta passage rate in the rumen that reduces the proteolytic activity of rumen microbiota (Sullivan and Martin, 1999; Bach et al., 2005).
Several studies have been conducted to determine the effect of SCFP on ruminal nitrogen metabolism. It has been reported that SCFP increases nitrogen utilization by rumen bacteria and reduces ammonia and methane emissions by dairy cows (Hristov et al., 2010). Tun et al. (2020) found that SCFP attenuated the reduction of nitrogen metabolism during the SARA challenge. Our results partly agreed with these previous studies, as we found that SCFP promoted predicted amino acid metabolism and carbohydrate metabolism, including ascorbate and aldarate metabolism, fructose and mannose metabolism, and inositol phosphate metabolism during SARA challenges. Additionally, we also observed that SCFP promoted the colonization of several fibrolytic bacteria, including Lachnospiraceae UCG-009, Treponema, unclassified Lachnospiraceae, and unclassified Ruminococcaceae during the SARA challenges. As fibrolytic bacteria have a high preference for ammonia (Bryant, 1973), the SCFP supplementation may have increased the nitrogen metabolism in our study.
Microbe-to-microbe association networks offer high-level insights into the changes in the global structure of microbial communities (Peschel et al., 2021). We found that only the SCFPb-2X treatment maintained the characteristics of the network inside the rumen solids microbial community under grain-based SARA challenges, indicating a more resilient and stable microbial community during this treatment. Besides, members from the Lachnospiraceae NK3A20 group, Phocaeicola, Anaerotruncus, Anaerofustis, Senegalimassilis, and unclassified Lachnospiraceae played important hub roles in the network derived from the SCFPb-2X treatment during SARA and had strong positive associations with richness and diversity of the microbiota in rumen solids microbiota. It is conceivable that an increased abundance of the above taxa will assist in maintaining the diversity and richness of the microbial community during dietary or metabolic stress periods.
5 Conclusions
Repeated grain-based SARA challenges reduced the richness and diversity and changed the β-diversity of the bacterial community in rumen solids. Grain-based SARA challenges also reduced the relative abundance of cellulolytic bacteria, such as members of Fibrobacteres phyla, reduced the interactions among rumen bacteria, and the diversity of hub taxa in solids microbiota. Supplementation with SCFPb increased the positive connectedness among microbial members in the solid fraction and the diversity of hub taxa. In particular, SCFPb-2X maintained the pre-SARA properties of the rumen solids microbial network during the SARA challenges. Grain-based SARA challenges inhibited the growth of several fibrolytic bacteria, while SCFP supplementation also promoted the growth of several fibrolytic taxa during SARA challenges. SARA challenges also reduced the predicted nitrogen metabolism in the microbiota of rumen solids, while SCFP supplementation attenuated this decline, which may prevent the reduction in microbial protein synthesis. Due to the limitations of the 16S rRNA gene sequencing to classify taxa at the species level and limitations of current databases when using CowPi to predict the rumen microbiome functionality, further shotgun metagenomic and metabolomic approaches are needed to better understand the consequences of SARA and benefits of SCFPs on the rumen microbiome.
Data availability statement
The rumen solids sequencing data were deposited into NCBI Sequence Read Archive (SRA) and are available under the Bioproject ID: PRJNA1139131.
Ethics statement
All procedures used in this study was approved by the University of Manitoba Animal Care Committee (Protocol # F14-038) and followed the guidelines of the Canadian Council for Animal Care (CCAC, 1993).
Author contributions
JG: Investigation, Methodology, Formal analysis, Visualization, Writing – original draft. ZZ: Formal analysis, Investigation, Methodology, Writing – review & editing. LG: Formal analysis, Methodology, Writing – review & editing, Resources, Supervision. MZ: Methodology, Writing – review & editing. IY: Writing – review & editing, Conceptualization. EK: Conceptualization, Writing – review & editing, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation. JP: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing, Data curation.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by grants from the Natural Science and Engineering Research Council (NSERC) of Canada Collaborative Research and Development (CRD) program, NSERC Discovery program, Dairy Farmers of Manitoba, Canada, and Diamond V, USA, to EK and JP.
Acknowledgments
We thank the staff of the Dairy Research Unit at the Glenlea Research Station, University of Manitoba, for the maintenance and care of the animals. We also thank the members of the Microbiome Laboratory at the University of Manitoba, including Drs. Paula Azevedo, Hein Min Tun, Shucong Li, and Hooman Derakhshani, as well as Dr. Anirikh Chakrabarti from Cargill Inc. for their assistance with laboratory and experimental setup and data analyses.
Conflict of interest
IY was employed by Diamond V.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1409659/full#supplementary-material
Abbreviations
ASV, amplicon sequence variants; CAZy, carbohydrate active enzymes; CoNet, correlation network; CP, crude protein; CSS, cumulative sum scaling; FDR, false discovery rate; MetaLonDA, Metagenomic Longitudinal Differential Abundance; NDF, detergent fiber; NetCoMi, Network Construction and comparison for Microbiome data; nMDS, Non-metric multidimensional scaling; PERMANOVA, Permutational multivariate analysis of variance; PERMDISP, permutational multivariate analysis of dispersion; PICRUSt, Phylogenetic Investigation of Communities by Reconstruction of Unobserved States; SARA, subacute ruminal acidosis; SCFP, Saccharomyces cerevisiae fermentation product; TMR, total mixed ration; VFA, volatile fatty acid.
References
Acharya, S., Pretz, J., Yoon, I., Scott, M., and Casper, D. (2017). Effects of Saccharomyces cerevisiae fermentation products on the lactational performance of mid-lactation dairy cows. Transl. Anim Sci. 1, 221–228. doi: 10.2527/tas2017.0028
Agler, M. T., Ruhe, J., Kroll, S., Morhenn, C., Kim, S.-T., Weigel, D., et al. (2016). Microbial hub taxa link host and abiotic factors to plant microbiome variation. PLoS Biol. 14:e1002352. doi: 10.1371/journal.pbio.1002352
Allen, M., and Ying, Y. (2012). Effects of Saccharomyces cerevisiae fermentation product on ruminal starch digestion are dependent upon dry matter intake for lactating cows. J. Dairy Sci. 95, 6591–6605. doi: 10.3168/jds.2012-5377
Allison, M. J., Bucklin, J. A., and Robinson, I. (1966). Importance of the isovalerate carboxylation pathway of leucine biosynthesis in the rumen. Appl. Microbiol. 14, 807–814. doi: 10.1128/am.14.5.807-814.1966
Assoumani, M., Vedeau, F., Jacquot, L., and Sniffen, C. (1992). Refinement of an enzymatic method for estimating the theoretical degradability of proteins in feedstuffs for ruminants. Anim. Feed Sci. Tech. 39, 357–368. doi: 10.1016/0377-8401(92)90054-A
Auffret, M. D., Dewhurst, R. J., Duthie, C.-A., Rooke, J. A., John Wallace, R., Freeman, T. C., et al. (2017). The rumen microbiome as a reservoir of antimicrobial resistance and pathogenicity genes is directly affected by diet in beef cattle. Microbiome 5, 1–11. doi: 10.1186/s40168-017-0378-z
Bach, A., Calsamiglia, S., and Stern, M. (2005). Nitrogen metabolism in the rumen. J. Dairy Sci. 88, E9–E21. doi: 10.3168/jds.S0022-0302(05)73133-7
Bäumler, A. J., and Sperandio, V. (2016). Interactions between the microbiota and pathogenic bacteria in the gut. Nature 535, 85–93. doi: 10.1038/nature18849
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R Stat Soc B. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Benson, A. K., Kelly, S. A., Legge, R., Ma, F., Low, S. J., Kim, J., et al. (2010). Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors. PNAS USA. 107, 18933–18938. doi: 10.1073/pnas.1007028107
Bickhart, D., and Weimer, P. (2018). Symposium review: host–rumen microbe interactions may be leveraged to improve the productivity of dairy cows. J. Dairy Sci. 101, 7680–7689. doi: 10.3168/jds.2017-13328
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nature Biotechnol. 37, 852–857. doi: 10.1038/s41587-019-0209-9
Brede, M., Orton, T., Pinior, B., Roch, F.-F., Dzieciol, M., Zwirzitz, B., et al. (2020). PacBio and Illumina MiSeq amplicon sequencing confirm full recovery of the bacterial community after subacute ruminal acidosis challenge in the RUSITEC system. Front. Microbiol. 11:1813. doi: 10.3389/fmicb.2020.01813
Brewer, M. T., Anderson, K. L., Yoon, I., Scott, M. F., and Carlson, S. A. (2014). Amelioration of salmonellosis in pre-weaned dairy calves fed Saccharomyces cerevisiae fermentation products in feed and milk replacer. Vet. Microbiol. 172, 248–255. doi: 10.1016/j.vetmic.2014.05.026
Bryant, M. P. (1973). Nutritional requirements of the predominant rumen cellulolytic bacteri. Federat. Proc. 32, 1809–1813.
Callaway, E., and Martin, S. (1997). Effects of a Saccharomyces cerevisiae culture on ruminal bacteria that utilize lactate and digest cellulose. J. Dairy Sci. 80, 2035–2044. doi: 10.3168/jds.S0022-0302(97)76148-4
Cameron, M., Klusmeyer, T., Lynch, G., Clark, J., and Nelson, D. (1991). Effects of urea and starch on rumen fermentation, nutrient passage to the duodenum, and performance of cows. J. Dairy Sci. 74, 1321–1336. doi: 10.3168/jds.S0022-0302(91)78288-X
Casper, D. P., and Schingoethe, D. J. (1989). Lactational response of dairy cows to diets varying in ruminal solubilities of carbohydrate and crude protein. J. Dairy Sci. 72, 928–941. doi: 10.3168/jds.S0022-0302(89)79186-4
CCAC (1993). Guide to the Care and Use of Experimental Animals. Vol. 1, 2nd Edn. Ottawa, ON: Canadian Council on Animal Care.
Cherdthong, A., Wanapat, M., Kongmun, P., Pilajun, R., and Khejornsart, P. (2010). Rumen fermentation, microbial protein synthesis and cellulolytic bacterial population of swamp buffaloes as affected by roughage to concentrate ratio. J. Anim. Vet. Adv. 9, 1667–1675. doi: 10.3923/javaa.2010.1667.1675
Clark, J., Klusmeyer, T., and Cameron, M. (1992). Symposium: Nitrogen metabolism and amino acid nutrition in dairy cattle: microbial protein synthesis and flows of nitrogen fractions to the duodenum of dairy cows. J. Dairy Sci. 75, 2304–2323. doi: 10.3168/jds.S0022-0302(92)77992-2
Clarke, K., and Gorley, R. (2015). Getting started with PRIMER v7. Plymouth: Plymouth Marine Laboratory, 20.
Cooper, R., Klopfenstein, T. J., Stock, R., Parrott, C., and Herold, D. (1997). “Effect of Rumensin and feed intake variation on ruminal pH,” in Nebraska Beef Cattle Reports (Lincoln: University of Nebraska), 430.
Craig, W. M., Broderick, G. A., and Ricker, D. B. (1987). Quantitation of microorganisms associated with the particulate phase of ruminal ingesta. J. Nutri. 117, 56–62. doi: 10.1093/jn/117.1.56
De Menezes, A. B., Lewis, E., 'Donovan, M., 'Neill, B. F., Clipson, N., and Doyle, E. M. (2011). Microbiome analysis of dairy cows fed pasture or total mixed ration diets. FEMS. Microb Ecol. 78, 256–265. doi: 10.1111/j.1574-6941.2011.01151.x
Derakhshani, H., Plaizier, J. C., De Buck, J., Barkema, H. W., and Khafipour, E. (2018). Association of bovine major histocompatibility complex (BoLA) gene polymorphism with colostrum and milk microbiota of dairy cows during the first week of lactation. Microbiome 6, 1–18. doi: 10.1186/s40168-018-0586-1
Derakhshani, H., Plaizier, J. C., De Buck, J., Barkema, H. W., and Khafipour, E. (2020). Composition and co-occurrence patterns of the microbiota of different niches of the bovine mammary gland: potential associations with mastitis susceptibility, udder inflammation, and teat-end hyperkeratosis. Anim. Microbiome 2, 1–17. doi: 10.1186/s42523-020-00028-6
Dias, A., Freitas, J., Micai, B., Azevedo, R., Greco, L., and Santos, J. (2018). Effects of supplementing yeast culture to diets differing in starch content on performance and feeding behavior of dairy cows. J. Dairy Sci. 101, 186–200. doi: 10.3168/jds.2017-13240
Diez-Gonzalez, F., Callaway, T. R., Kizoulis, M. G., and Russell, J. B. (1998). Grain feeding and the dissemination of acid-resistant Escherichia coli from cattle. Science 281, 1666–1668. doi: 10.1126/science.281.5383.1666
Duysburgh, C., Miclotte, L., Green, J. B., Watts, K. T., Sardi, M. I., Chakrabarti, A., et al. (2024). Saccharomyces cerevisiae derived postbiotic alters gut microbiome metabolism in the human distal colon resulting in immunomodulatory potential in vitro. Front. Microbiol. 15:1358456. doi: 10.3389/fmicb.2024.1358456
Eckel, E. F., and Ametaj, B. N. (2016). Invited review: Role of bacterial endotoxins in the etiopathogenesis of periparturient diseases of transition dairy cows. J. Dairy Sci. 99, 5967–5990. doi: 10.3168/jds.2015-10727
Faust, K., and Raes, J. (2012). Microbial interactions: from networks to models. Nat. Rev Microb. 10, 538–550. doi: 10.1038/nrmicro2832
Faust, K., and Raes, J. (2016). CoNet app: inference of biological association networks using Cytoscape. F1000Res 5, 1519. doi: 10.12688/f1000research.9050.1
Fernando, S. C., Purvis, H., Najar, F., Sukharnikov, L., Krehbiel, C., Nagaraja, T., et al. (2010). Rumen microbial population dynamics during adaptation to a high-grain diet. Appl. Environ Microb. 76, 7482–7490. doi: 10.1128/AEM.00388-10
Firkins, J., and Yu, Z. (2015). Ruminant nutrition symposium: how to use data on the rumen microbiome to improve our understanding of ruminant nutrition. J. Anim Sci. 93, 1450–1470. doi: 10.2527/jas.2014-8754
Gagen, E. J., Padmanabha, J., Denman, S. E., and McSweeney, C. S. (2015). Hydrogenotrophic culture enrichment reveals rumen Lachnospiraceae and Ruminococcaceae acetogens and hydrogen-responsive Bacteroidetes from pasture-fed cattle. FEMS. Microbiol Lett. 362:fnv104. doi: 10.1093/femsle/fnv104
Ganda, E., Chakrabarti, A., Sardi, M. I., Tench, M., Kozlowicz, B. K., Norton, S. A., et al. (2023). Saccharomyces cerevisiae fermentation product improves robustness of equine gut microbiome upon stress. Front. Vet Sci. 10:1134092. doi: 10.3389/fvets.2023.1134092
Gozho, G., Krause, D., and Plaizier, J. (2007). Ruminal lipopolysaccharide concentration and inflammatory response during grain-induced subacute ruminal acidosis in dairy cows. J. Dairy Sci. 90, 856–866. doi: 10.3168/jds.S0022-0302(07)71569-2
Gozho, G. N., Plaizier, J. C., Krause, D. O., Kennedy, A. D., and Wittenberg, K. M. (2005). Subacute ruminal acidosis induces ruminal lipopolysaccharide endotoxin release and triggers an inflammatory response. J. Dairy Sci. 88, 1399–1403. doi: 10.3168/jds.S0022-0302(05)72807-1
Gruninger, R., Ribeiro, G., Cameron, A., and McAllister, T. (2019). Invited review: application of meta-omics to understand the dynamic nature of the rumen microbiome and how it responds to diet in ruminants. Animal 13, 1843–1854. doi: 10.1017/S1751731119000752
Guan, L. L., Nkrumah, J. D., Basarab, J. A., and Moore, S. S. (2008). Linkage of microbial ecology to phenotype: correlation of rumen microbial ecology to cattl's feed efficiency. FEMS. Microbiol Lett. 288, 85–91. doi: 10.1111/j.1574-6968.2008.01343.x
Guo, J., Xu, L., Khalouei, H., Fehr, K., Senaratne, V., Ghia, J., et al. (2022). Saccharomyces cerevisiae fermentation products reduce bacterial endotoxin concentrations and inflammation during grain-based subacute ruminal acidosis in lactating dairy cows. J. Dairy Sci. 105, 2354–2368. doi: 10.3168/jds.2021-20572
Guo, J., Zhang, Z., Guan, L. L., Yoon, I., Plaizier, J. C., and Khafipour, E. (2024). Postbiotics from Saccharomyces cerevisiae fermentation stabilize microbiota in rumen liquid digesta during grain-based subacute ruminal acidosis (SARA) in lactating dairy cows. J. Anim. Sci. Biotechnol. 15:101. doi: 10.1186/s40104-024-01056-x
Henderson, G., Cox, F., Ganesh, S., Jonker, A., Young, W., and Janssen, P. H. (2015). Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci. Rep. 5, 1–15. doi: 10.1038/srep14567
Hristov, A., Varga, G., Cassidy, T., Long, M., Heyler, K., Karnati, S., et al. (2010). Effect of Saccharomyces cerevisiae fermentation product on ruminal fermentation and nutrient utilization in dairy cows. J. Dairy Sci. 93, 682–692. doi: 10.3168/jds.2009-2379
Huo, W., Zhu, W., and Mao, S. (2014). Impact of subacute ruminal acidosis on the diversity of liquid and solid-associated bacteria in the rumen of goats. World J. Microb Biot. 30, 669–680. doi: 10.1007/s11274-013-1489-8
Ji, S., Zhang, H., Yan, H., Azarfar, A., Shi, H., Alugongo, G., et al. (2017). Comparison of rumen bacteria distribution in original rumen digesta, rumen liquid and solid fractions in lactating Holstein cows. J. Anim Sci Biot. 8, 1–7. doi: 10.1186/s40104-017-0142-z
Khafipour, E., Krause, D. O., and Plaizier, J. C. (2009a). A grain-based subacute ruminal acidosis challenge causes translocation of lipopolysaccharide and triggers inflammation. J. Dairy Sci. 92, 1060–1070. doi: 10.3168/jds.2008-1389
Khafipour, E., Li, S., Plaizier, J. C., and Krause, D. O. (2009b). Rumen microbiome composition determined using two nutritional models of subacute ruminal acidosis. Appl. Environ Microb. 75, 7115–7124. doi: 10.1128/AEM.00739-09
Khafipour, E., Li, S., Tun, H. M., Derakhshani, H., Moossavi, S., and Plaizier, J. C. (2016). Effects of grain feeding on microbiota in the digestive tract of cattle. Anim Frontiers. 6, 13–19 doi: 10.2527/af.2016-0018
Khafipour, E., Plaizier, J. C., Aikman, P. C., and Krause, D. O. (2011). Population structure of rumen Escherichia coli associated with subacute ruminal acidosis (SARA) in dairy cattle. J. Dairy Sci. 94, 351–360. doi: 10.3168/jds.2010-3435
Khalouei, H., Seranatne, V., Fehr, K., Guo, J., Yoon, I., Khafipour, E., et al. (2020). Effects of Saccharomyces cerevisiae fermentation products and subacute ruminal acidosis on feed intake, fermentation, and nutrient digestibilities in lactating dairy cows. Can J. Anim Sci. 101, 143–157. doi: 10.1139/cjas-2020-0018
Kleen, J., Hooijer, G., Rehage, J., and Noordhuizen, J. (2003). Subacute ruminal acidosis (SARA): a review. J. Vet Med A. 50, 406–414. doi: 10.1046/j.1439-0442.2003.00569.x
Lesmeister, K., Heinrichs, A. J., and Gabler, M. (2004). Effects of supplemental yeast (Saccharomyces cerevisiae) culture on rumen development, growth characteristics, and blood parameters in neonatal dairy calves. J. Dairy Sci. 87, 1832–1839. doi: 10.3168/jds.S0022-0302(04)73340-8
Li, F., and Guan, L. L. (2017). Metatranscriptomic profiling reveals linkages between the active rumen microbiome and feed efficiency in beef cattle. Appl. Environ Microb. 83, e00061–e00017. doi: 10.1128/AEM.00061-17
Li, S., Yoon, I., Scott, M., Khafipour, E., and Plaizier, J. (2016). Impact of Saccharomyces cerevisiae fermentation product and subacute ruminal acidosis on production, inflammation, and fermentation in the rumen and hindgut of dairy cows. Anim. Feed Sci Tech. 211, 50–60. doi: 10.1016/j.anifeedsci.2015.10.010
Liu, J., Bai, Y., Liu, F., Kohn, R. A., Tadesse, D. A., Sarria, S., et al. (2022). Rumen microbial predictors for short-chain fatty acid levels and the grass-fed regimen in angus cattle. Animals 12:2995. doi: 10.3390/ani12212995
Louis, P., Young, P., Holtrop, G., and Flint, H. J. (2010). Diversity of human colonic butyrate-producing bacteria revealed by analysis of the butyryl-CoA: acetate CoA-transferase gene. Environ Microb. 12, 304–314. doi: 10.1111/j.1462-2920.2009.02066.x
Lourenco, J. M., Kieran, T. J., Seidel, D. S., Glenn, T. C., Silveira, M. F. d., et al. (2020). Comparison of the ruminal and fecal microbiotas in beef calves supplemented or not with concentrate. PLoS ONE 15:e0231533. doi: 10.1371/journal.pone.0231533
Mahmoud, A. H., Slate, J. R., Hong, S., Yoon, I., and McGill, J. L. (2020). Supplementing a Saccharomyces cerevisiae fermentation product modulates innate immune function and ameliorates bovine respiratory syncytial virus infection in neonatal calves. J. Anim Sci. 98:skaa252. doi: 10.1093/jas/skaa252
Manirajan, B. A., Maisinger, C., Ratering, S., Rusch, V., Schwiertz, A., Cardinale, M., et al. (2018). Diversity, specificity, co-occurrence and hub taxa of the bacterial–fungal pollen microbiome. FEMS. Microb Ecol. 94:fiy112. doi; 10.1093/femsec/fiy112
Mao, S., Zhang, R., Wang, D., and Zhu, W. (2013). Impact of subacute ruminal acidosis (SARA) adaptation on rumen microbiota in dairy cattle using pyrosequencing. Anaerobe 24, 12–19. doi: 10.1016/j.anaerobe.2013.08.003
Martin, C., Millet, L., Fonty, G., and Michalet-Doreau, B. (2001). Cereal supplementation modified the fibrolytic activity but not the structure of the cellulolytic bacterial community associated with rumen solid digesta. Reprod. Nutr Develop. 41, 413–424. doi: 10.1051/rnd:2001141
McCann, J. C., Luan, S., Cardoso, F. C., Derakhshani, H., Khafipour, E., and Loor, J. J. (2016). Induction of subacute ruminal acidosis affects the ruminal microbiome and epithelium. Front. Microb. 7:701. doi: 10.3389/fmicb.2016.00701
McCann, J. C., Wiley, L. M., Forbes, T. D., Rouquette Jr, F. M., and Tedeschi, L. O. (2014). Relationship between the rumen microbiome and residual feed intake-efficiency of Brahman bulls stocked on bermudagrass pastures. PLoS ONE 9:e91864. doi: 10.1371/journal.pone.0091864
Meehan, C. J., and Beiko, R. G. (2014). A phylogenomic view of ecological specialization in the Lachnospiraceae, a family of digestive tract-associated bacteria. Genome. Biol Evol. 6, 703–713. doi: 10.1093/gbe/evu050
Metwally, A. A., Yang, J., Ascoli, C., Dai, Y., Finn, P. W., and Perkins, D. L. (2018). MetaLonDA: a flexible R package for identifying time intervals of differentially abundant features in metagenomic longitudinal studies. Microbiome 6, 1–12. doi: 10.1186/s40168-018-0402-y
Michalet-Doreau, B., Fernandez, I., Peyron, C., Millet, L., and Fonty, G. (2001). Fibrolytic activities and cellulolytic bacterial community structure in the solid and liquid phases of rumen contents. Reprod. Nutr Develop. 41, 187–194. doi: 10.1051/rnd:2001122
Monteiro, H. F., and Faciola, A. P. (2020). Ruminal acidosis, bacterial changes, and lipopolysaccharides. J. Anim Sci. 98:skaa248. doi: 10.1093/jas/skaa248
Mould, F., and Ørskov, E. (1983). Manipulation of rumen fluid pH and its influence on cellulolysis in sacco, dry matter degradation and the rumen microflora of sheep offered either hay or concentrate. Anim. Feed Sci Tech. 10, 1–14. doi: 10.1016/0377-8401(83)90002-0
Mullins, C., Mamedova, L., Carpenter, A., Ying, Y., Allen, M., Yoon, I., et al. (2013). Analysis of rumen microbial populations in lactating dairy cattle fed diets varying in carbohydrate profiles and Saccharomyces cerevisiae fermentation product. J. Dairy Sci. 96, 5872–5881. doi: 10.3168/jds.2013-6775
Myer, P. R., Smith, T. P., Wells, J. E., Kuehn, L. A., and Freetly, H. C. (2015). Rumen microbiome from steers differing in feed efficiency. PLoS ONE 10:e0129174. doi: 10.1371/journal.pone.0129174
Nagaraja, T., and Titgemeyer, E. (2007). Ruminal acidosis in beef cattle: the current microbiological and nutritional outlook. J. Dairy Sci. 90, E17–E38. doi: 10.3168/jds.2006-478
Nocek, J., and Russell, J. (1988). Protein and energy as an integrated system. Relationship of ruminal protein and carbohydrate availability to microbial synthesis and milk production. J. Dairy Sci. 71, 2070–2107. doi: 10.3168/jds.S0022-0302(88)79782-9
Parks, D. H., Tyson, G. W., Hugenholtz, P., and Beiko, R. G. (2014). STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics 30, 3123–3124. doi: 10.1093/bioinformatics/btu494
Peschel, S., Müller, C. L., von Mutius, E., Boulesteix, A.-L., and Depner, M. (2021). NetCoMi: network construction and comparison for microbiome data in R. Brief. Bioinform. 22:bbaa290. doi: 10.1093/bib/bbaa290
Petri, R. M., Schwaiger, T., Penner, G. B., Beauchemin, K. A., Forster, R. J., McKinnon, J. J., et al. (2013). Characterization of the core rumen microbiome in cattle during transition from forage to concentrate as well as during and after an acidotic challenge. PLoS ONE 8:e83424. doi: 10.1371/journal.pone.0083424
Plaizier, J., Khafipour, E., Li, S., Gozho, G., and Krause, D. (2012). Subacute ruminal acidosis (SARA), endotoxins and health consequences. Anim. Feed Sci Tech. 172, 9–21. doi: 10.1016/j.anifeedsci.2011.12.004
Plaizier, J., Mesgaran, M. D., Derakhshani, H., Golder, H., Khafipour, E., Kleen, J., et al. (2018). Enhancing gastrointestinal health in dairy cows. Animal 12, s399–s418. doi: 10.1017/S1751731118001921
Plaizier, J., Mulligan, F., Neville, E., Guan, L., Steele, M., and Penner, G. (2022). Invited review: Effect of subacute ruminal acidosis on gut health of dairy cows. J. Dairy Sci. 105, 7141–7160. doi: 10.3168/jds.2022-21960
Plaizier, J. C., Krause, D. O., Gozho, G. N., and McBride, B. W. (2008). Subacute ruminal acidosis in dairy cows: the physiological causes, incidence and consequences. Vet J. 176, 21–31. doi: 10.1016/j.tvjl.2007.12.016
Plaizier, J. C., Li, S., Danscher, A. M., Derakshani, H., Andersen, P. H., and Khafipour, E. (2017). Changes in microbiota in rumen digesta and feces due to a grain-based subacute ruminal acidosis (SARA) challenge. Microb. Ecol. 74, 485–495. doi: 10.1007/s00248-017-0940-z
Pope, P., Smith, W., Denman, S., Tringe, S., Barry, K., Hugenholtz, P., et al. (2011). Isolation of Succinivibrionaceae implicated in low methane emissions from Tammar wallabies. Science 333, 646–648. doi: 10.1126/science.1205760
Poppy, G., Rabiee, A., Lean, I., Sanchez, W., Dorton, K., and Morley, P. (2012). A meta-analysis of the effects of feeding yeast culture produced by anaerobic fermentation of Saccharomyces cerevisiae on milk production of lactating dairy cows. J. Dairy Sci. 95, 6027–6041. doi: 10.3168/jds.2012-5577
Qi, W., Xue, M. Y., Jia, M. H., Zhang, S., Yan, Q., and Sun, H. Z. (2024). Understanding the functionality of the rumen microbiota: searching for better opportunities for rumen microbial manipulation. Anim Biosci 37, 370. doi: 10.5713/ab.23.0308
Russell, J. (1998). The importance of pH in the regulation of ruminal acetate to propionate ratio and methane production in vitro. J. Dairy Sci. 81, 3222–3230. doi: 10.3168/jds.S0022-0302(98)75886-2
Russell, J. B., and Rychlik, J. L. (2001). Factors that alter rumen microbial ecology. Science 292, 1119–1122. doi: 10.1126/science.1058830
Salminen, S., Collado, M. C., Endo, A., Hill, C., Lebeer, S., Quigley, E. M. M., et al. (2021). The International Scientific Association of Probiotics and Prebiotics (ISAPP) consensus statement on the definition and scope of postbiotics. Nat. Rev. Gastroenterol. Hepatol. 18, 649–667. doi: 10.1038/s41575-021-00440-6
Shin, N.-R., Whon, T. W., and Bae, J.-W. (2015). Proteobacteria: microbial signature of dysbiosis in gut microbiota. Trends Biotechnol. 33, 496–503. doi: 10.1016/j.tibtech.2015.06.011
Stern, M. D., and Hoover, W. H. (1979). Methods for determining and factors affecting rumen microbial protein synthesis: a review. J. Anim Sci. 49, 1590–1603. doi: 10.2527/jas1979.4961590x
Storm, E., and Ørskov, E. (1983). The nutritive value of rumen micro-organisms in ruminants: 1. Large-scale isolation and chemical composition of rumen micro-organisms. Br. J. Nutri. 50, 463–470. doi: 10.1079/BJN19830114
Sullivan, H., and Martin, S. (1999). Effects of a Saccharomyces cerevisiae culture on in vitro mixed ruminal microorganism fermentation. J. Dairy Sci. 82, 2011–2016. doi: 10.3168/jds.S0022-0302(99)75438-X
Tománková, O., and Kopečný, J. (1995). Prediction of feed protein degradation in the rumen with bromelain. Anim. Feed Sci Tech. 53, 71–80. doi: 10.1016/0377-8401(94)00735-R
Trosvik, P., and de Muinck, E. J. (2015). Ecology of bacteria in the human gastrointestinal tract—identification of keystone and foundation taxa. Microbiome 3, 1–12. doi: 10.1186/s40168-015-0107-4
Tun, H. M., Li, S., Yoon, I., Meale, S. J., Azevedo, P. A., Khafipour, E., et al. (2020). Saccharomyces cerevisiae fermentation products (SCFP) stabilize the ruminal microbiota of lactating dairy cows during periods of a depressed rumen pH. BMC. Vet Res. 16, 1–17. doi: 10.1186/s12917-020-02437-w
Urrutia, N., Bomberger, R., Matamoros, C., and Harvatine, K. (2019). Effect of dietary supplementation of sodium acetate and calcium butyrate on milk fat synthesis in lactating dairy cows. J. Dairy Sci. 102, 5172–5181. doi: 10.3168/jds.2018-16024
Urrutia, N. L., and Harvatine, K. J. (2017). Acetate dose-dependently stimulates milk fat synthesis in lactating dairy cows. J. Nutri. 147, 763–769. doi: 10.3945/jn.116.245001
Vailati-Riboni, M., Coleman, D., Lopreiato, V., Alharthi, A., Bucktrout, R., Abdel-Hamied, E., et al. (2021). Feeding a Saccharomyces cerevisiae fermentation product improves udder health and immune response to a Streptococcus uberis mastitis challenge in mid-lactation dairy cows. J. Anim Sci Biotechnol. 12, 1–19. doi: 10.1186/s40104-021-00560-8
Valente, T. N. P., Sampaio, C. B., da Silva Lima, E., Deminicis, B. B., Cezário, A. S., and dos Santos, W. B. R. (2017). Aspects of acidosis in ruminants with a focus on nutrition: a review. J. Agric Sci. 9, 90. doi: 10.5539/jas.v9n3p90
Wei, T., Simko, V., Levy, M., Xie, Y., Jin, Y., and Zemla, J. (2017). Package ‘corrplot'. Statistician. 56:e24.
Wilkinson, T. J., Huws, S. A., Edwards, J. E., Kingston-Smith, A. H., Siu-Ting, K., Hughes, M., et al. (2018). CowPI: a rumen microbiome focussed version of the PICRUSt functional inference software. Front Microb. 9:1095. doi: 10.3389/fmicb.2018.01095
Williams, A., Withers, S. E., and Strachan, N. (1989). Postprandial variations in the activity of polysaccharide-degrading enzymes in microbial populations from the digesta solids and liquor fractions of rumen contents. J. Appl Bacteriol. 66, 15–26. doi: 10.1111/j.1365-2672.1989.tb02449.x
Xiao, J., Alugongo, G., Chung, R., Dong, S., Li, S., Yoon, I., et al. (2016). Effects of Saccharomyces cerevisiae fermentation products on dairy calves: ruminal fermentation, gastrointestinal morphology, and microbial community. J. Dairy Sci. 99, 5401–5412. doi: 10.3168/jds.2015-10563
Xiao, J., Alugongo, G. M., Ji, S., Wu, Z., Dong, S., Li, S., et al. (2018). Effects of Saccharomyces cerevisiae fermentation products on the microbial community throughout the gastrointestinal tract of calves. Animals 9:4. doi: 10.3390/ani9010004
Xue, M.-Y., Sun, H.-Z., Wu, X.-H., Liu, J.-X., and Guan, L. L. (2020). Multi-omics reveals that the rumen microbiome and its metabolome together with the host metabolome contribute to individualized dairy cow performance. Microbiome 8, 1–19. doi: 10.1186/s40168-020-00819-8
Yang, W. Z., Beauchemin, K. A., and Rode, L. M. (2001). Effects of grain processing, forage to concentrate ratio, and forage particle size on rumen pH and digestion by dairy cows. J. Dairy Sci. 84, 2203–2216. doi: 10.3168/jds.S0022-0302(01)74667-X
Yost, W. M., Young, J. W., Schmidt, S. P., and McGilliard, A. D. (1977). Gluconeogenesis in ruminants: propionic acid production from a high-grain diet fed to cattle. J. Nutri. 107, 2036–2043. doi: 10.1093/jn/107.11.2036
Zebeli, Q., and Metzler-Zebeli, B. (2012). Interplay between rumen digestive disorders and diet-induced inflammation in dairy cattle. Res. Vet Sci. 93, 1099–1108. doi: 10.1016/j.rvsc.2012.02.004
Zhang, Z., Tun, H. M., Li, R., Gonzalez, B. J., Keenes, H. C., Nyachoti, C. M., et al. (2018). Impact of xylanases on gut microbiota of growing pigs fed corn-or wheat-based diets. Anim Nutri. 4, 339–350. doi: 10.1016/j.aninu.2018.06.007
Zhou, M., Chen, Y., and Guan, L. (2015). “Rumen bacteria,” in Rumen Microbiology: From Evolution to Revolution (Cham: Springer), 79–95. doi: 10.1007/978-81-322-2401-3_6
Zhu, W., Wei, Z., Xu, N., Yang, F., Yoon, I., Chung, Y., et al. (2017). Effects of Saccharomyces cerevisiae fermentation products on performance and rumen fermentation and microbiota in dairy cows fed a diet containing low quality forage. J. Anim Sci Biotechnol. 8, 1–9. doi: 10.1186/s40104-017-0167-3
Keywords: dairy cows, postbiotics, Saccharomyces cerevisiae fermentation product, SARA, rumen solids microbiota
Citation: Guo J, Zhang Z, Guan LL, Zhou M, Yoon I, Khafipour E and Plaizier JC (2024) Postbiotics from Saccharomyces cerevisiae fermentation stabilize rumen solids microbiota and promote microbial network interactions and diversity of hub taxa during grain-based subacute ruminal acidosis (SARA) challenges in lactating dairy cows. Front. Microbiol. 15:1409659. doi: 10.3389/fmicb.2024.1409659
Received: 30 March 2024; Accepted: 03 July 2024;
Published: 16 August 2024.
Edited by:
Zheng Ruan, Nanchang University, ChinaReviewed by:
Todd Riley Callaway, University of Georgia, United StatesMichael D. Flythe, United States Department of Agriculture (USDA), United States
Hao Lizhuang, Qinghai University, China
Copyright © 2024 Guo, Zhang, Guan, Zhou, Yoon, Khafipour and Plaizier. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ehsan Khafipour, Ehsan.Khafipour@gmail.com; Jan C. Plaizier, kees.plaizier@umanitoba.ca
†Present addresses: Zhengxiao Zhang, College of Food and Biological Engineering, Jimei University, Xiamen, Fujian, China
Le Luo Guan, Faculty of Land and Food Systems, University of British Columbia, Vancouver, BC, Canada
Ehsan Khafipour, Cargill Animal Nutrition, Cargill Inc., Wayzata, MN, United States