Abstract
The current study used a culture-independent methodology to investigate the rumen microbiome composition in two genetically divergent groups of spring-calving dairy cows, high (€218) and low (€157) Economic Breeding Index (EBI), grazing two sward treatments—perennial ryegrass (grass-only) or perennial ryegrass and white clover (grass-clover)—at three time points across spring, summer, and autumn of a single grazing season. The analysis indicated that the EBI status had no significant effect on the rumen microbial community within the statistical power of this study. Beta diversity between the microbiomes was different (p < 0.001) between the two sward treatments only in autumn, when the clover proportion was highest (50.2%). Season had a significant effect on microbiome beta diversity across sward treatments (p < 0.001). There were only minor differences in the composition of the rumen microbiomes between the two sward treatments. Many bacterial genera were differentially abundant between spring and the two later time points. Bacterial genera that were more abundant in spring were positively correlated with rumen propionate levels, while those more abundant in summer and autumn were negatively correlated with propionate and positively correlated with acetate and butyrate. Methanogenic archaeal abundance was greater in summer and autumn compared to spring, and they were negatively correlated with propionate and positively correlated with methane (CH4) production. The results of this study demonstrate that genetic selection using the EBI does not affect the rumen microbial community and the core rumen microbial community is similar in cows grazing either grass-only or grass-clover swards. The results also demonstrate that the rumen bacterial community shifts across the grazing season, providing more favorable conditions for methanogenesis in summer and autumn compared to spring.
1 Introduction
The rumen of dairy cows has a complex and diverse microbial ecosystem that plays an important role in the digestion of plant material from which energy is provided to its host (Van Soest, 1994). This symbiotic relationship is essential for the production of nutrient-dense milk for human consumption (Mottet et al., 2018). However, hydrogen (H2) is a by-product of rumen fermentation that is removed to maintain redox balance and efficient fermentation in the rumen (). Hydrogen removal is achieved primarily by methanogenic archaea, which utilize excess H2 in the reduction of carbon dioxide (CO2) to methane (CH4), the majority of which is ultimately eructated by the animal into the atmosphere (). Enteric CH4 is a potent greenhouse gas (GHG) and a major contributor to agricultural emissions (IPCC 2019).
Due to the fact that microbial populations associate with different substrates, diet has a major influence on the rumen microbiota of their host (). In temperate regions with consistent annual grass growth, the diet of dairy cows largely consists of grazed pasture (; ). In recent years, the inclusion of legumes in pasture-based systems, such as white clover (Trifolium repens L.; clover), has gained increasing importance, as their nitrogen (N) fixation ability offers the potential to offset the amount of artificial N fertilizer required (). In addition, grass-clover swards have improved feed quality compared to grass-only swards, increasing dry matter intake (DMI), which often leads to increased animal performance (; ). However, the increased DMI associated with clover has been reported to increase enteric CH4 emissions (; ), although some studies have reported that clover reduces CH4 yield (CH4 per unit of DMI; ; ). Previous studies have reported subtle differences in the rumen microbial communities of animals grazing grass-only and grass-clover swards (; Woodmartin et al., 2024). Furthermore, grass-only and grass-clover swards both vary in quality across the grazing season (). Previous studies have also noted variation in CH4 emissions from dairy cows over the grazing season, the cause of which has been postulated to be partly due to shifts in pasture quality (; ). The rumen microbial community of cows grazing grass-clover swards has been previously shown to undergo minor alterations over the grazing season with shifting sward quality and changes in lactation stage (Noel et al., 2017). However, these alterations have not been characterized in grass-only swards or with a focus on methanogenesis. To better understand seasonal variation in the CH4 yield of pasture-based dairy systems, the dynamics of the rumen microbial community structure across the grazing season in both grass-only and grass-clover swards must be further investigated.
In addition to diet, there is some evidence for host influences on the rumen microbial community, whereby differences have been observed between dairy cow breeds (; ; ) and between production, efficiency, and CH4 phenotypes within breeds (Tong et al., 2018; ; Zhu et al., 2021; ). However, differences in microbiome composition between phenotypes are not always consistent, nor are they necessarily genetically driven (). Previous research has shown that dairy cows ranking higher on the Irish total merit index, the Economic Breeding Index (EBI; www.icbf.com), have greater milk solids (milk fat + protein) production and/or greater feed efficiency (milk solids per kg of DMI) (; ; ). The EBI is a total merit index that combines various genetic traits and subindexes into a single monetary value that represents an animal’s potential profitability within the context of Irish pasture-based dairy systems (Veerkamp et al., 2002). Research has demonstrated that improvements in performance associated with the EBI are not associated with greater CH4 emissions. However, there is currently limited information on whether selected genetic differences can influence the rumen microbial community, particularly in the context of enteric CH4 emissions.
Therefore, this study aimed to compare the rumen bacterial and methanogenic archaeal communities of high EBI and low EBI spring-calving Holstein Friesian dairy cows grazing grass-only or grass-clover swards at different stages of the grazing season.
2 Materials and methods
2.1 Experiment design
All experimental procedures conducted on the animals in the study were approved by the Health Products and Regulatory Authority (Dublin, Ireland) under project authorization AE19132/P138. The reported experiment was undertaken in 2022 at the Teagasc, Animal and Grassland Research and Innovation Centre, Moorepark, Fermoy, Co. Cork, Ireland. This study was part of a larger farm systems study described by , which involved two Holstein Friesian dairy cow genotypes of high and low genetic merit (€218 and €157 EBI, respectively). The cows in each EBI group were distributed evenly across two sward treatments: perennial ryegrass (Lolium perenne L.; PRG) receiving 225 kg artificial N/ha/yr (grass-only) and PRG with white clover (Trifolium repens L.) receiving 150 kg artificial N/ha/yr (grass-clover). The experiment consisted of three 14-day measurement periods: spring (27th of April to the 10th of May), summer (6th of July to the 19th of July), and autumn (21 September to 4 October). The study parameters comprised 28 spring-calving Holstein-Friesian dairy cows with average “days in milk” values of 70 in spring, 140 in summer, and 217 in autumn (SD = 15.6). A total of 14 cows were randomly selected from each sward treatment, balanced for EBI group, parity, milk production, milk composition, body weight, and body condition score. The cows were managed in a rotational grazing system similar to that described by . The target pre-grazing herbage mass (>4 cm) was 1,300–1,600 kg DM/ha, and daily herbage allocations (>4 cm) were approximately 17 kg DM/cow/day. Residency time in each allocation was determined by a target post-grazing sward height of 4.0 to 4.5 cm measured using a rising plate meter (Jenquip Ltd.; Fielding, New Zealand).
2.2 Animal measurements
Animal measurements have been previously outlined by . In brief, CH4 emissions were measured using two Greenfeed units (C-lock Inc.; South Dakota, USA), which were swapped between the treatments twice during each measurement period. The cows received a target concentrate supplementation of approximately 1 kg/cow from the Greenfeed units per day. Milk production was measured daily, and milk fat and protein content were measured from subsequent AM and PM milk samples collected twice during each measurement period. Dry matter intake was estimated using the n-alkane method described by from days 6 to 10 during each measurement period. All measurements were then averaged per day across the measurement period.
2.3 Rumen fluid
Rumen fluid samples were collected from each animal on the final day of each measurement period using a transesophageal sampling device (FLORA rumen scoop; Guelph, ON, Canada). Rumen fluid was passed through cheesecloth, and 8 mL was added to 2 mL of 50% trichloroacetic acid. The samples were then stored at −18 °C until analysis. The remaining non-acidified sample was stored separately at −80 °C. The acidified samples were analyzed for volatile fatty acids (VFAs; acetate, propionate, butyrate, valerate, isobutyrate, isovalerate) and ammonia concentration, as previously described by . Insufficient fluid was collected from one cow; therefore, it was excluded from the analysis. Three cows were also excluded due to health reasons not related to the experiment, and an additional two were excluded as they had insufficient methane data. A breakdown of the number of subjects used in the final analysis is presented in Table 1.
Table 1
| GO | GC | Total | |||
|---|---|---|---|---|---|
| High EBI | Low EBI | High EBI | Low EBI | ||
| Spring | 6 | 7 | 5 | 4 | 22 |
| Summer | 6 | 7 | 5 | 4 | 22 |
| Autumn | 6 | 7 | 5 | 4 | 22 |
| Total | 18 | 21 | 15 | 12 | 66 |
The number of subjects (n) analyzed within each breed group and sward treatment across the three measurement periods.
GO, grass-only; GC, grass-clover; EBI, economic breeding index.
2.4 DNA extraction, library preparation, and sequencing
DNA was extracted from 250 mg of the non-acidified rumen fluid samples using the repeated bead beating and column purification method described by Yu and Morrison (2004). DNA quality was assessed on agarose gels, and DNA concentration was quantified using a Nanodrop 1,000 spectrophotometer. The DNA samples were then standardized to a volume of 25 μL at a concentration of 20 ng/μl of DNA per sample. They were then shipped on ice to Novogene (UK) Co., Ltd. for PCR and sequencing.
Amplicons targeting the V4 region of the 16S rRNA gene were generated using the 515F/806R primers with Illumina technology. Sequencing libraries were generated, and indexes were added. The libraries were checked using Qubit and real-time PCR for quantification and a bioanalyzer for size distribution detection. The quantified libraries were pooled and sequenced on Illumina platforms, according to effective library concentration and the required data amount. Paired-end reads were merged using FLASH (). Quality filtering on the raw tags was performed using fastp (). The SILVA database1 was used as a reference to detect chimeric sequences. The chimeric sequences were removed using the vsearch package (). Effective tags were cleaned using DADA2 to obtain initial amplicon sequence variants (ASVs). Taxonomic annotation was performed using QIIME 2 software, primarily with the SILVA database. For sequences that could not be annotated using SILVA, NCBI BLAST analysis was used to supplement the taxonomic information. The ASVs were separated into their individual kingdoms (Bacteria and Archaea). The archaeal ASVs were further classified against the RefSeq database2 using DADA2 in R to determine the proportion belonging to the Methanobrevibacter SGMT (smithii, gottschalkii, millerae, and thaueri) and RO (ruminantium and olleyae) clades, as defined by . Separate ASV tables were then constructed for the bacteria and archaea.
2.5 Data and statistical analysis
To investigate alpha diversity, the bacterial and archaeal ASV tables were each converted into phyloseq objects in R using phyloseq (). Each was then rarefied to its respective minimum sequencing depth (43,009 and 604 reads, respectively), following which observed ASVs and Shannon diversity were calculated for each sample’s respective bacterial and archaeal populations. Due to poor species-level classification, bacterial and archaeal population composition was aggregated to the genus level and the clade level in the case of Methanobrevibacter. Principal components and the Procrustes test were used to check if the aggregated data were similar in structure with the original data (p < 0.001) using the vegan R package (). The aggregated tables were converted into phyloseq objects and rarefied to their respective minimum read depths after aggregation (35,502 and 604 reads for bacteria and archaea, respectively).
2.5.1 Bacterial alpha and beta diversity
Bacterial alpha diversity was compared between the genotypes, the sward treatments and seasons using the Kruskal–Wallis test, followed by pairwise comparisons with the Dunn test adjusted using the Benjamini–Hochberg (BH) method. Only bacterial genera present in >5% of the samples were included in the beta diversity analysis. Principal coordinates analysis (PCoA) based on the Bray–Curtis dissimilarity was performed to visualize differences in bacterial beta diversity at the genus level. Co-inertia plots were generated, displaying BH-adjusted significant Spearman correlations between bacterial genera and the principal axes. A plot showing significant correlations (p < 0.05) of animal performance, gaseous emissions, rumen fermentation characteristics, and experimental factors with the PCoA ordination was generated using the envfit function from the vegan package in R. The experimental factors tested were genotype, sward treatment, season, and parity. Statistical differences in beta diversity were identified using PERMANOVA with the adonis2 function from the vegan package in R.
2.5.2 Differential abundance testing
The analysis of composition of microbiomes (ANCOM) test was carried out on the non-rarefied bacterial ASVs present in >5% of the samples. Differences were identified at the genus level between the genotypes, sward treatments, and seasons using the ANCOM-BC package in R with a BH adjustment ().
2.5.3 Archaeal population
The ratio of archaeal to bacterial reads, archaeal alpha diversity, and the relative abundance of the archaeal genera and clades were analyzed using the MIXED procedure in SAS, with genotype, sward treatment, measurement period, and their interactions as fixed effects (Woodmartin et al., 2024). A linear mixed model was used, as all archaeal data were normally distributed and taxonomic diversity was low both at the genus (two genera) and clade levels (two clades). Individual animal was included as a random effect, with measurement period as a repeated measure and a compound symmetry correlation structure. Other effects tested were parity and calving day of the year, as well as all possible interactions, but none of which were included in the final model as they were not significant for any variable tested (p > 0.05). The results were expressed as least square means, and the multiple comparisons were adjusted using the Tukey–Kramer method.
2.5.4 Partial correlations
Spearman partial correlations, adjusted for genotype, sward treatment, measurement period, and parity, were used to investigate the relationship between microbial taxa and CH4 emissions, animal performance, and VFA concentrations using R.
3 Results
Grazing results, herbage chemical composition, enteric emissions, animal performance, and rumen VFA results have been reported previously by . A total of 6,386,864 reads were obtained from the 16S amplicon sequencing analysis, with a mean (standard deviation) number of reads per sample of 96,770 (26,160.3). Following quality filtering, merging, and removal of chimeric sequences, there was a total of 6,103,700 reads, with a mean of 92,480 (24,979.7) reads per sample. These reads were mapped to 8,510 ASVs, with a mean of 1,292 (340.8) ASVs per sample.
3.1 Rumen bacterial microbiome analysis
Across both sward treatments and all measurement periods, Prevotella (32.2%) was the most prevalent bacterial genus, followed by Olsenella (11.0%), the Christensenellaceae R7 group (9.49%), and Kandleria (6.62%; Supplementary Figure 1). On average, Erysipelotrichaceae UCG 002 was also a dominant genus (5.52%); however, this genus was more prevalent in grass-only swards during spring (20.3%) and was skewed by two samples with particularly high abundances (78.5 and 58.7%, respectively). There was no effect of the EBI group on the bacterial alpha diversity. The effects of sward treatment and measurement period on bacterial alpha-diversity are presented in Figure 1. Grass-clover swards had a greater (p < 0.05) number of observed ASVs than grass-only swards, but there was no difference in Shannon diversity. Within the grass-only treatment, season tended (p = 0.08) to have an effect on observed ASVs, while there was no effect within the grass-clover treatment. Within the grass-only treatment, spring had lower Shannon diversity than summer and autumn (p < 0.05), while summer and autumn did not differ from each other. Within the grass-clover treatment, spring tended to have lower Shannon diversity than summer (p = 0.08) and autumn (p = 0.09), with no significant difference between summer and autumn.
Figure 1
Beta-diversity of the EBI groups and sward treatments across the measurement periods was investigated at the genus level using PCoA (Figure 2A). PERMANOVA showed that the EBI group had no effect on beta diversity. Grass-only and grass-clover treatments tended to separate diagonally across both the first and second principal axes in spring (p = 0.052) and also separated diagonally across both axes in the opposite direction in autumn (p < 0.01). The separation in spring was primarily associated with two divergent values within the grass-only treatment. Figure 2B shows genera associated with these separations based on their Spearman correlations with both the first and second principal axes (p < 0.05). A total of two genera of the Erysipelotrichaceae family were negatively correlated with the first axis and positively correlated with the second axis, while 11 genera were negatively correlated with both axes, five of which belonged to the Lachnospiraceae family. The Prevotellaceae YAB2003 group showed the strongest negative correlation with both axes (ρ = −0.80 and −0.52 for the first and second axes, respectively). In total, nine genera were positively correlated with both axes, two of which were members of the Lachnospiraceae family which showed the strongest correlations. PERMANOVA also showed separation within each sward treatment across the seasons, primarily along the first principal axis (p < 0.001). Figure 2C shows the genera that were Spearman correlated (ρ > 0.5, p < 0.05) with the first axis only. A total of 16 genera were negatively correlated with the first principal axis, with Prevotella 7, Erysipelotrichaceae UCG 002, and the Prevotellaceae YAB2003 group showing the strongest correlations (ρ = −0.95, −0.85, and −0.80, respectively). A total of 33 genera were positively correlated with the first principal axis, with the Christensenellaceae R7 group, the Oscillospiraceae NK4A214 group, Mogibacterium, Flexilinea, and Eggerthellaceae DNF00809 showing the strongest correlations (ρ = 0.95, 0.91, 0.86, 0.86, and 0.80, respectively). Figure 2D shows the animal measurements and experimental factors that correlated (p < 0.05) with the rumen bacterial community. Methane, CH4 per kg of milk solids, CH4 yield, rumen ammonia, acetate, and isobutyrate were positively correlated with the first principal axis toward autumn, while feed efficiency (milk solids per kg DMI) and propionate were negatively correlated with the first principal axis toward spring. Valerate was associated with two divergent values within the grass-only treatment in spring. Isobutyrate was associated with the grass-clover treatment in autumn. Of the experimental factors tested, only season and treatment were significant (p < 0.05), and their positions matched those presented in Figure 2A.
Figure 2
The analysis of the microbiome composition found that no genera were differentially abundant between the two EBI groups. A total of 21 genera were differentially abundant between the sward treatments, although only four genera were abundant in the rumen (>1%; Table 2). Erysipelotrichaceae UCG 002 and Prevotella 7 were reduced, while the Christensenellaceae R7 group and Oscillospiraceae NK4A214 group were increased in the grass-clover treatment compared to the grass-only treatment (p < 0.05). Within individual measurement periods, there were no differences between the sward treatments. Compared with spring, 54 and 55 genera were differentially abundant in summer and autumn, respectively (p < 0.05; Supplementary Table S1). In agreement with the PCoA, Prevotella 7, Erysipelotrichaceae UCG 002, and the Prevotellaceae YAB2003 group were reduced, while the Christensenellaceae R7 group, Oscillospiraceae NK4A214 group, Olsenella, and Butyrivibrio were increased in both summer and autumn compared to spring. In total, there were 12 and 14 differentially abundant genera with relative abundances >1% in summer and autumn, respectively, compared with spring (Table 2). Autumn had 24 differentially abundant genera compared to summer. Among the most abundant (>1%), Kandleria and the Eubacterium cellulosolvens group decreased, while Succiniclasticum increased.
Table 2
| Genus | T | G | S | |||||||
|---|---|---|---|---|---|---|---|---|---|---|
| GC – GO | Low EBI – High EBI | Summer – Spring | Autumn – Spring | Autumn – Summer | ||||||
| W | p-value | W | p-value | W | p-value | W | p-value | W | p-value | |
| Prevotella | 0.187 | 0.99 | −1.901 | 0.57 | −0.067 | 1.00 | −1.717 | 1.00 | −1.616 | 1.00 |
| Olsenella | 0.464 | 0.83 | −0.093 | 1.00 | 4.943 | <0.001 | 4.763 | <0.001 | 0.329 | 1.00 |
| Christensenellaceae R 7 group | 2.641 | <0.05 | 0.518 | 0.98 | 3.784 | <0.001 | 5.883 | <0.001 | 1.990 | 0.09 |
| Kandleria | −2.105 | 0.11 | 0.635 | 0.96 | 0.189 | 1.00 | −5.921 | <0.001 | −5.194 | <0.001 |
| NK4A214 group | 3.217 | <0.05 | −0.259 | 1.00 | 5.033 | <0.001 | 6.340 | <0.001 | 1.486 | 0.25 |
| Erysipelotrichaceae UCG 002 | −5.232 | <0.001 | 1.128 | 0.78 | −4.770 | <0.001 | −4.647 | <0.001 | 0.646 | 0.92 |
| Prevotella 7 | −4.315 | <0.001 | 0.458 | 1.00 | −5.366 | <0.001 | −7.287 | <0.001 | −1.662 | 0.18 |
| Prevotellaceae UCG 001 | −1.087 | 0.46 | −0.594 | 0.96 | 0.437 | 1.00 | −1.225 | 1.00 | −1.501 | 1.00 |
| Succiniclasticum | 2.131 | 0.11 | −1.074 | 0.78 | −2.704 | <0.05 | 1.101 | 0.49 | 3.689 | <0.01 |
| Butyrivibrio | 0.760 | 0.62 | 0.055 | 1.00 | 3.695 | <0.01 | 0.785 | 0.77 | −2.167 | 0.09 |
| Lachnospiraceae NK3A20 group | 1.138 | 0.44 | 0.600 | 0.96 | 0.021 | 1.00 | −0.876 | 1.00 | −0.752 | 1.00 |
| Ruminococcus | 0.746 | 0.63 | 0.608 | 0.96 | −2.975 | <0.05 | −4.468 | 0.00 | −1.716 | 0.16 |
| Acetitomaculum | 1.974 | 0.14 | 1.009 | 0.81 | 1.414 | 0.37 | 2.739 | <0.05 | 1.275 | 0.37 |
| Eubacterium cellulosolvens group | −2.377 | 0.07 | 0.967 | 0.81 | −2.080 | 0.07 | −5.342 | <0.001 | −3.059 | <0.01 |
| Rikenellaceae RC9 gut group | 0.794 | 0.62 | −1.650 | 0.57 | 4.009 | <0.001 | 2.685 | <0.05 | −0.505 | 1.00 |
| Ruminococcus gauvreauii group | 0.157 | 1.00 | 1.276 | 0.71 | 0.864 | 0.69 | 3.133 | <0.05 | 2.157 | 0.09 |
The effect of sward type (T), dairy cow genotype (G), and season (S) on the differential abundance of the most abundant bacterial genera (>1%) in the rumen fluid of the dairy cows.
The results are expressed as the effect size (W) for the grass-clover (GC) treatment compared to the grass-only (GO) treatment, low Economic Breeding Index (EBI) compared to high EBI, and pairwise comparisons between spring, summer, and autumn.
3.2 Rumen archaeal community
The effects of EBI group, sward treatment, and measurement period on the rumen archaeal population are presented in Table 3. There was no effect of the EBI group on the archaeal community. The ratio of archaea to bacteria was greatest (p < 0.01) in the grass-clover treatment, although this was only observed in autumn (p < 0.001; Figure 3). The ratio of archaea to bacteria was lowest (p < 0.001) in spring, intermediate in summer, and greatest in autumn. Observed archaeal ASVs were greater (p < 0.01) in the grass-clover treatment than in the grass-only treatment. There was also an interaction with the measurement period, where grass-only and grass-clover treatments were only different in autumn (p < 0.001; Figure 3). Observed archaeal ASVs also increased (p < 0.001) from spring to summer to autumn across the sward treatments. Archaeal Shannon diversity tended to be greater (p < 0.1) in the grass-clover treatment than in the grass-only treatment. Shannon diversity was reduced in spring (p < 0.001) compared to both summer and autumn, but summer and autumn did not differ from each other. Neither the sward treatment nor the measurement period had any effect on the relative abundance of archaeal genera. Across all samples, Methanobrevibacter and Methanosphaera accounted for 91.48 and 8.51% of archaea, respectively. The relative abundance of the Methanobrevibacter SGMT clade was not different between the swards, but the RO clade was slightly greater (p < 0.05) in the grass-only treatment compared to the grass-clover treatment due to an increase in spring. The SGMT clade was reduced (p < 0.001) in spring compared to summer and autumn, which did not differ from each other. The RO clade was increased (p < 0.001) in spring compared to summer and autumn, while autumn tended to be lower than summer (p < 0.1).
Table 3
| Item | T | SE | G | SE | S | SE | p-value | |||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| GO | GC | High | Low | Spring | Summer | Autumn | T1 | G | S | T × S | ||||
| Archaea:Bacteria | 0.036 | 0.055 | 0.0547 | 0.044 | 0.046 | 0.0054 | 0.021a | 0.047b | 0.067c | 0.0673 | <0.05 | 0.81 | <0.001 | <0.001 |
| Alpha diversity | ||||||||||||||
| Observed ASVs | 10.01 | 12.43 | 0.607 | 10.99 | 11.45 | 0.558 | 8.51a | 11.38b | 13.77c | 0.500 | <0.01 | 0.56 | <0.001 | <0.001 |
| Shannon | 1.83 | 1.92 | 0.039 | 1.88 | 1.87 | 0.036 | 1.74a | 1.91b | 1.98b | 0.034 | 0.11 | 0.89 | <0.001 | 0.52 |
| Genus (%) | ||||||||||||||
| Methanobrevibacter | 92.09 | 90.87 | 0.787 | 91.25 | 91.71 | 0.724 | 91.09 | 91.54 | 91.83 | 0.793 | 0.24 | 0.65 | 0.78 | 0.27 |
| Methanosphaera | 7.90 | 9.11 | 0.789 | 8.74 | 8.27 | 0.725 | 8.91 | 8.44 | 8.17 | 0.792 | 0.25 | 0.65 | 0.78 | 0.28 |
| Clade (%) | ||||||||||||||
| SGMT | 62.64 | 68.52 | 2.859 | 66.80 | 64.36 | 2.629 | 52.17a | 69.68b | 74.89b | 2.403 | 0.12 | 0.51 | <0.001 | 0.09 |
| RO | 27.55 | 20.39 | 2.599 | 22.39 | 25.55 | 2.390 | 37.89a | 19.70b | 14.32b | 2.232 | <0.05 | 0.35 | <0.001 | <0.05 |
The effect of sward type (T), dairy cow genotype (G), and season (S) on the rumen archaeal community.
1GO = grass-only, GC = grass-clover. a,b,cValues in the same row of S not sharing a common superscript are significantly different (p < 0.05; Tukey–Kramer adjusted). SE, standard error of the mean.
Figure 3
3.3 Relationship between rumen bacteria and enteric CH4 emissions, animal performance, and rumen VFAs
Partial correlations between the most relatively abundant rumen bacterial genera (>1%) with enteric emissions, animal performance, and rumen VFA proportions are shown in Tables 4, 5. The remaining genera are shown in Supplementary Table S2. Methane production (218–446 g/d) was positively correlated (p < 0.05) with Acetitomaculum, the Christensenellaceae R7 group, and Olsenella and negatively correlated with Prevotella 7. Methane yield (11.8–24.5 g/kg) was positively correlated (p < 0.05) with Acetitomaculum, and the ratio of milk solids to DMI was positively correlated with Acetitomaculum and the Rikenellaceae RC9 gut group. Total VFA concentration (116–233 mmoL/L) was positively correlated (p < 0.05) with Prevotella and negatively correlated with Acetitomaculum, Erysipelotrichaceae UCG 002, and the Lachnospiraceae NK3A20 group. Acetate (64.3–77.1%) was positively correlated (p < 0.05) with Butyrivibrio, the Eubacterium cellulosolvens group, the Lachnospiraceae NK3A20 group, Prevotellaceae UCG 001, the Rikenellaceae RC9 gut group, Ruminococcus, and Succiniclasticum and negatively correlated with Erysipelotrichaceae UCG 002 and Prevotella 7. Propionate (12.6–23.8%) had positive relationships (p < 0.05) with Erysipelotrichaceae UCG 002 and Prevotella 7 and negative relationships with Acetitomaculum, Butyrivibrio, the Christensenellaceae R7 group, the Lachnospiraceae NK3A20 group, the Oscillospiraceae NK4A214 group, Olsenella, the Rikenellaceae RC9 gut group, and the Ruminococcus gauvreauii group. Butyrate (5.6–10.1%) had a positive correlation (p < 0.05) with Prevotella and negative relationships with Butyrivibrio, the Eubacterium cellulosolvens group, the Lachnospiraceae NK3A20 group, and Ruminococcus. The correlations between the ratio of acetate to propionate (2.7–6.0 mmoL/L) and bacterial genera were the inverse of those of propionate (p < 0.05).
Table 4
| CH4 | CH4/MS | CH4/DMI | MS | DMI | MS/DMI | |
|---|---|---|---|---|---|---|
| Acetitomaculum | 0.32** | 0.09 | 0.28* | 0.2 | 0.01 | 0.24* |
| Butyrivibrio | 0.19 | −0.09 | 0.03 | 0.22† | 0.13 | 0.14 |
| Christensenellaceae R7 group | 0.29* | 0.2 | 0.22† | 0.13 | 0.07 | 0.09 |
| Erysipelotrichaceae UCG 002 | −0.23† | 0.09 | −0.12 | −0.2 | −0.16 | −0.11 |
| Eubacterium cellulosolvens group | −0.07 | −0.03 | −0.13 | −0.13 | −0.01 | −0.09 |
| Kandleria | −0.09 | 0.05 | −0.05 | −0.22† | −0.11 | −0.14 |
| Lachnospiraceae NK3A20 group | 0.16 | −0.02 | 0.13 | 0.07 | −0.09 | 0.18 |
| Oscillospiraceae NK4A214 group | 0.13 | −0.04 | 0.09 | 0.15 | 0.01 | 0.15 |
| Olsenella | 0.25* | 0.08 | 0.19 | 0.16 | 0.02 | 0.08 |
| Prevotella | 0.21† | 0.04 | 0.14 | 0.16 | 0.14 | 0.06 |
| Prevotella 7 | −0.26* | 0.01 | −0.18 | −0.22† | −0.11 | −0.09 |
| Prevotellaceae UCG 001 | 0.21† | 0.08 | 0.19 | 0.08 | 0.01 | 0.05 |
| Rikenellaceae RC9 gut group | 0.04 | −0.15 | 0.08 | 0.16 | −0.11 | 0.25* |
| Ruminococcus | −0.08 | −0.19 | −0.15 | 0 | 0.03 | 0.04 |
| Ruminococcus gauvreauii group | 0.14 | 0.17 | 0.13 | 0.04 | 0.08 | −0.05 |
| Succiniclasticum | −0.01 | 0.06 | 0.09 | −0.13 | −0.11 | −0.04 |
Partial correlations between rumen bacterial genera and enteric CH4 emissions and animal performance1.
Significant relationships are indicated by ‘*’ (p < 0.05), ‘**’ (p < 0.01), and ‘***’ (p < 0.001), and tendencies are indicated by ‘†’ (P < 0.1).
MS, milk solids (fat + protein); DMI, dry matter intake.
Table 5
| Total VFA | Acetate | Propionate | Butyrate | A:P | |
|---|---|---|---|---|---|
| Acetitomaculum | −0.28* | 0.13 | −0.47*** | 0.09 | 0.44*** |
| Butyrivibrio | −0.22† | 0.62*** | −0.35** | −0.47*** | 0.37** |
| Christensenellaceae R7 group | −0.19 | 0.17 | −0.63*** | 0.19 | 0.64*** |
| Erysipelotrichaceae UCG 002 | −0.26* | −0.46*** | 0.34** | −0.09 | −0.27* |
| Eubacterium cellulosolvens group | −0.12 | 0.25* | 0.15 | −0.47*** | −0.11 |
| Kandleria | −0.12 | 0.01 | 0.20 | −0.21† | −0.18 |
| Lachnospiraceae NK3A20 group | −0.48*** | 0.28* | −0.32** | −0.34** | 0.39*** |
| Oscillospiraceae NK4A214 group | −0.17 | 0.02 | −0.38** | 0.14 | 0.41*** |
| Olsenella | −0.11 | 0.02 | −0.26* | 0.16 | 0.22† |
| Prevotella | 0.66*** | 0.13 | 0.09 | 0.36** | −0.18 |
| Prevotella 7 | −0.14 | −0.41*** | 0.52*** | −0.11 | −0.42*** |
| Prevotellaceae UCG 001 | 0.19 | 0.34** | −0.23† | 0.07 | 0.23† |
| Rikenellaceae RC9 gut group | −0.11 | 0.28* | −0.33** | −0.09 | 0.36** |
| Ruminococcus | −0.05 | 0.25* | 0.07 | −0.30* | −0.04 |
| Ruminococcus gauvreauii group | −0.13 | 0.19 | −0.44*** | 0.08 | 0.44*** |
| Succiniclasticum | −0.01 | 0.25* | −0.02 | −0.18 | 0.00 |
Partial correlations between rumen bacterial genera and rumen volatile fatty acids (VFA)1.
Significant relationships are indicated by ‘*’ (p < 0.05), ‘**’ (p < 0.01), and ‘***’ (p < 0.001), and tendencies are indicated by ‘†’ (p < 0.1).
A: P, acetate to propionate ratio.
3.4 Relationship between rumen archaea and enteric CH4 emissions, animal performance, and rumen VFAs
The partial correlations of the ratio of archaea to bacteria, archaeal genera, and Methanobrevibacter clades with enteric CH4 emissions, animal performance traits, and rumen VFA proportions are shown in Tables 6, 7. The ratio of archaea to bacteria had a positive relationship with CH4 yield and CH4 per kg of MS (p < 0.05), as well as a positive tendency with CH4 production (p < 0.1). At the genus level, Methanobrevibacter abundance had a positive relationship with CH4 production and the ratio of acetate to propionate and a negative relationship with propionate (p < 0.05). Methanosphaera abundance had opposite relationships with Methanobrevibacter (p < 0.05). The SGMT clade had a positive relationship with CH4 production, DMI, total VFA concentration, and butyrate (p < 0.05). The RO clade tended to have a negative relationship with CH4 production and butyrate (p < 0.1) and had a negative relationship with DMI and total VFA concentration (p < 0.05).
Table 6
| CH4 | CH4/MS | CH4/DMI | MS | DMI | MS/DMI | |
|---|---|---|---|---|---|---|
| Archaea:Bacteria | 0.24† | 0.26* | 0.31* | −0.03 | −0.11 | 0.08 |
| Genus | ||||||
| Methanobrevibacter | 0.29* | 0.04 | 0.19 | 0.18 | 0.09 | 0.14 |
| Methanosphaera | −0.29* | −0.04 | −0.19 | −0.18 | −0.09 | −0.14 |
| Clade | ||||||
| SGMT | 0.28* | −0.04 | 0.03 | 0.24† | 0.33** | 0.02 |
| RO | −0.22† | 0.04 | 0.01 | −0.20 | −0.32** | 0.01 |
Partial correlations between rumen archaea and enteric CH4 emissions and animal performance1.
Significant relationships are indicated by ‘*’ (p < 0.05), ‘**’ (p < 0.01), and ‘***’ (p < 0.001), and tendencies are indicated by ‘†’ (p < 0.1).
MS, milk solids (fat + protein); DMI, dry matter intake; SGMT, Methanobrevibacter smithii, gottschalkii, millerae, and thaueri; RO, Methanobrevibacter ruminantium and olleyae.
Table 7
| Total VFA | Acetate | Propionate | Butyrate | A:P | |
|---|---|---|---|---|---|
| Archaea:Bacteria | −0.23† | −0.01 | −0.46*** | 0.19 | 0.45*** |
| Genus | |||||
| Methanobrevibacter | −0.16 | −0.16 | −0.34** | 0.21† | 0.29* |
| Methanosphaera | 0.16 | 0.16 | 0.34** | −0.21† | −0.29* |
| Clade | |||||
| SGMT | 0.46*** | 0.09 | −0.22† | 0.25* | 0.09 |
| RO | −0.55*** | −0.05 | 0.16 | −0.24† | −0.03 |
Partial correlations between rumen archaea and rumen volatile fatty acids (VFA)1.
Significant relationships are indicated by ‘*’ (p < 0.05), ‘**’ (p < 0.01), and ‘***’ (p < 0.001), and tendencies are indicated by ‘†’ (p < 0.1).
A: P, acetate to propionate ratio; SGMT, Methanobrevibacter smithii, gottschalkii, millerae, and thaueri; RO, Methanobrevibacter ruminantium and olleyae.
4 Discussion
The current study examined ruminal bacteria and archaea populations of genetically divergent dairy cows grazing grass-only or grass-clover swards at different stages of the grazing season. Similar to previous findings in autumn, there was no difference in bacterial Shannon diversity between the sward treatments across all three measurement periods (). In addition, the abundances of the main rumen bacterial genera were similar between the cows grazing the two sward treatments in the current study. Prevotella dominated the microbiome of the cows grazing both sward types, which is consistently reported across a wide variety of diets (; ; Zhu et al., 2021). As such, Prevotella is often labeled a generalist, owing to its ability to degrade a wide variety of plant cell wall and storage polysaccharides (). However, there was a difference in Erysipelotrichaceae UCG002 between the two swards, although this genus was primarily associated with two divergent values in the grass-only treatment in spring. The co-inertia plots also showed that these values were associated with Sharpea and Erysipelotrichaceae UCG007. Sharpea and other members of the Erysipelotrichaceae family are associated with lactic acid production (; ; Zhou et al., 2024); therefore, these animals may have been experiencing rumen dysbiosis due to low rumen pH at the time of sampling (). Low rumen pH is common among spring-calving dairy cows grazing highly digestible spring pastures, although proliferation of lactic acid bacteria is not (; ). However, rumen pH was not measured in the current study; therefore, no conclusions can be drawn.
Similar to , we observed an overall difference in microbial community composition between the two sward treatments in autumn. This was likely due to the increased clover proportions (50.2%) during this measurement period, providing a greater diversity of substrates available in the rumen, such as pectin and its subsequent metabolites (Ulyatt, 1971; Van Soest, 1994). The digestion of pectin in animals grazing grass-clover swards has been previously associated with methanol production, a substrate for methylotrophic methanogenesis (). Previous studies have linked pectin-derived methanol with a greater abundance of the methylotrophic methanogen, Methanosphaera, in animals consuming grass-clover swards (; Woodmartin et al., 2024). No such relationship was observed in the current study. Generally, feeds with lower neutral detergent fiber (NDF) levels, such as the grass-clover sward in the current study, are associated with reduced CH4 yield, as they are associated with lower rumen H2 levels (Van Soest, 1994; ). However, this was not observed in our previously reported data (). We previously postulated that pectin digestion may also contribute to hydrogenotrophic methanogenesis in the rumen (), as it has been reported to produce acetate and formate when digested (). Formate and H2 associated with acetate are utilized by Methanobrevibacter in hydrogenotrophic methanogenesis (). This may explain why the grass-clover sward does not have the lower CH4 yield that would be expected of a feed with lower NDF levels. In general, archaeal abundance was greater in the grass-clover sward during autumn in the current study. It is unclear why this occurred and why it was not associated with an increase in CH4 yield compared to the grass-only sward (). It could theoretically be a result of differences in digestive kinetics between the two swards (Wang et al., 2016). Although the digestive products of pectin are similar to those of NDF, its digestion rate is much faster and may contribute to an increased digestion rate in cows grazing grass-clover swards (; Van Soest, 1994). A limitation of the current study is that the rumen samples were collected from all cows at one time point only (01:30 p.m.) during each measurement period. The spike in methanogen abundance observed in the grass-clover treatment during autumn may be a result of more rapid digestion of their morning grazing compared to the grass-only treatment and may not be representative of methanogen abundance across the entire day (Wang et al., 2016; ). This is important considering that cows have a partial preference for clover over PRG, which is reported to be greater in the morning compared to later in the day (). Future investigation is required into the rumen digestive kinetics of cows grazing grass-clover swards and their associations with the rumen microbial community.
The composition of the rumen microbiome has a major influence on the end products of ruminal fermentation (Russell and Hespell, 1981). Roehe et al. (2016) noted a greater abundance of methanogenic genes in the rumen metatranscriptome of cattle phenotypically high for CH4 production. These authors also suggested that archaeal abundance is under host genetic control and is an accurate predictor of CH4 output. Therefore, the fact that we observed no differences in the rumen microbial communities of the two dairy cow genotypes in the current study could be expected based on their similar rumen fermentation characteristics and enteric CH4 emissions (). These results suggest that the improved feed efficiency and greater milk constituents of the high EBI group are largely a result of physiological, as opposed to digestive, mechanisms, although a greater sample size could potentially reveal some subtle differences in the microbiome between the two genotypes. Nonetheless, our observation is similar to the findings of studies with comparable sample sizes investigating cows divergent in residual feed intake, where no differences in the rumen microbial community were observed (; ). Despite the similar microbial populations in the current study, future research should investigate the rumen metagenome of genetically divergent dairy cows to identify potential associations between host genetics, phenotypic performance, and microbial gene content, which can differ significantly among strains of the same microbial species ().
A novel aspect of the current study is its investigation of the rumen microbial community of spring-calving dairy cows at different stages of the grazing season. Noel et al. (2017) reported minor seasonal fluctuations in the rumen microbiome composition of grazing dairy cows, albeit that study focused on fiber-adherent microbes and included only four subjects. The current study is important in the context of previous research that reported variation in CH4 yield across the grazing season, with particularly lower levels observed in spring (; ). We acknowledge that the effect of season cannot be separated from lactation stage within such systems; however, previous studies on consistent total mixed ration (TMR) diets have reported only a minor effect of lactation stage on CH4 yield and the rumen microbial community (; ; Zhu et al., 2021; ). Across various analyses in the current study, we observed significant differences in both rumen bacterial and archaeal communities across seasons. With the exception of grass-clover swards in autumn, the ratios of archaea to bacteria across the three seasons in the current study align with our CH4 yield data () and those reported by previous studies (; ). Our results suggest that differences in archaeal abundance across seasons may be driven by shifts in rumen VFA proportions.
Changes in sward quality across the grazing season are associated with a shift in rumen VFA proportions, from greater propionate levels in spring to greater acetate levels in summer and autumn (). Previous data indicate that methanogens are negatively associated with propionate in the rumen (), which is in agreement with the results obtained from the current study. There are a number of probable reasons for this association: firstly, propionate production at the expense of acetate and butyrate production reduces H2 production and, subsequently, availability for methanogen growth (; Wang et al., 2014). Secondly, propionate production is generally associated with more rapidly digestible feeds, which reduce rumen pH and increase rumen passage rates, both of which are detrimental for methanogen growth (Van Soest, 1994; Van Kessel and Russell, 1996; ). Furthermore, the H2−utilizing bacterial genus Acetitomaculum () showed a similar trend in abundance to the archaea in the current study, which further supports the idea that their abundances are related to rumen H2 concentrations across seasons. Elevated propionate levels are driven by higher sward water-soluble carbohydrate levels as a result of lower NDF content, which is indicative of high-quality spring pasture (; Taweel et al., 2005). This shift in sward quality was associated with a shift in the rumen bacterial community, particularly between spring and the other two seasons, in the current study. The two bacterial genera with the greatest reduction in abundance from spring to the later seasons were Erysipelotrichaceae UCG002 and Prevotella 7. The Erysipelotrichaceae family has been shown to be highly abundant in grazing dairy cows in spring, a pattern associated with lower CH4 yield compared to a TMR diet (; ). Species in the Prevotella 7 genus are predicted to digest non-structural carbohydrates, such as starch, and produce propionate (), which is in line with the partial correlations observed in the current study. Furthermore, the bacterial genera with the greatest increase in abundance in summer and autumn compared to spring in the current study—Olsnella, the Christensenellaceae R7 group, the Oscillospiraceae NK4A214 group, the Rikenellaceae RC9 group, and Butyrivibrio—have been associated with forage digestion and/or acetate production (; Kraatz et al., 2011; ; ; ). This was further substantiated by the partial correlations observed in the current study, where each genus was negatively correlated with propionate and positively correlated with either acetate or butyrate, both of which would theoretically enhance H2 availability for methanogenesis ().
Our results show that the core rumen microbiome of grazing dairy cows is robust to differences in sward composition, facilitating similar levels of animal performance across grass-only and grass-clover swards (; ). In addition, based on our data, genetic selection using the EBI has no effect on rumen bacterial and archaeal populations. The results from the current study also indicate that archaeal abundance has the strongest association with the CH4 yield of any microbiota parameter tested and that there is a significant change in their abundance across the grazing season, irrespective of sward treatment, aligning with seasonal variations in CH4 yield. These changes seem to be largely driven by changes in sward quality across the grazing season, shifting the composition of the rumen bacterial community, which alters rumen fermentation characteristics. Future investigations into the rumen metagenome and the expression of genes related to hydrogen production and methanogenesis would further improve our understanding of this relationship.
Statements
Data availability statement
The datasets presented in this study are publicly available. This data can be found at: https://www.ncbi.nlm.nih.gov/, accession number PRJNA1273030.
Ethics statement
The animal study was approved by The Teagasc Animal Ethics Committee. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
CD: Investigation, Writing – review & editing, Formal analysis, Visualization, Writing – original draft, Data curation. AD: Methodology, Writing – review & editing, Formal analysis. PO’T: Writing – review & editing, Methodology. TO’C: Supervision, Writing – review & editing, Funding acquisition, Conceptualization. DM: Writing – review & editing, Methodology, Formal analysis. DH: Writing – review & editing, Conceptualization, Funding acquisition. HI: Writing – review & editing, Data curation. FB: Funding acquisition, Conceptualization, Writing – review & editing, Supervision. BL: Supervision, Project administration, Conceptualization, Writing – review & editing, Funding acquisition, Writing – original draft.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded through the Irish Dairy Levy administered by Dairy Research Ireland as well as a grant (16/RC/3835; VistaMilk) from Research Ireland funded by the Department of Agriculture, Food and Marine on behalf of the Government of Ireland. The financial contributions of Carberry, Ornua, Dairygold, and Tirlán are also gratefully acknowledged. This research forms part of the Teagasc Climate Centre.
Acknowledgments
The authors would like to acknowledge the Moorepark farm staff for their care and management of the animals used in this experiment.
Conflict of interest
The 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.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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.2025.1642486/full#supplementary-material
References
1
AsmaZ.SylvieC.LaurentC.JérômeM.ChristopheK.OlivierB.et al. (2013). Microbial ecology of the rumen evaluated by 454 GS FLX pyrosequencing is affected by starch and oil supplementation of diets. FEMS Microbiol. Ecol.83, 504–514. doi: 10.1111/1574-6941.12011
2
BainbridgeM. L.CersosimoL. M.WrightA.-D. G.KraftJ. (2016). Rumen bacterial communities shift across a lactation in Holstein, Jersey and Holstein × Jersey dairy cows and correlate to rumen function, bacterial fatty acid composition and production parameters. FEMS Microbiol. Ecol.92:fiw059. doi: 10.1093/femsec/fiw059
3
Betancur-MurilloC. L.Aguilar-MarínS. B.JovelJ. (2023). Prevotella: a key player in ruminal metabolism. Microorganisms11:1. doi: 10.3390/microorganisms11010001
4
BokulichN.MillsD. A. (2012). Next-generation approaches to the microbial ecology of food fermentations. BMB Rep.45:389. doi: 10.5483/bmbrep.2012.45.7.148
5
ChenX.SuX.LiJ.YangY.WangP.YanF.et al. (2021). Real-time monitoring of ruminal microbiota reveals their roles in dairy goats during subacute ruminal acidosis. npj Biofilms Microb7:45. doi: 10.1038/s41522-021-00215-6
6
De FilippisF.PasolliE.TettA.TaralloS.NaccaratiA.De AngelisM.et al. (2019). Distinct genetic and functional traits of human intestinal Prevotella copri strains are associated with different habitual diets. Cell Host Microbe25, 444–453.e3. doi: 10.1016/j.chom.2019.01.004
7
de MenezesA. B.LewisE.O’DonovanM.O’NeillB. F.ClipsonN.DoyleE. M. (2011). Microbiome analysis of dairy cows fed pasture or total mixed ration diets. FEMS Microbiol. Ecol.78, 256–265. doi: 10.1111/j.1574-6941.2011.01151.x
8
DillonP.StakelumG. (1989). Herbage and dosed alkanes as a grass measurement technique for dairy cows. Ir. J. Agric. Res.28:104.
9
DineenM.DelabyL.GillilandT.McCarthyB. (2018). Meta-analysis of the effect of white clover inclusion in perennial ryegrass swards on milk production. J. Dairy Sci.101, 1804–1816. doi: 10.3168/jds.2017-12586
10
DoyleN.MbandlwaP.KellyW. J.AttwoodG.LiY.RossR. P.et al. (2019). Use of lactic acid Bacteria to reduce methane production in ruminants, a critical review. Front. Microbiol.10:2207. doi: 10.3389/fmicb.2019.02207
11
DuškováD.MarounekM. (2001). Fermentation of pectin and glucose, and activity of pectin-degrading enzymes in the rumen bacterium Lachnospira multiparus. Lett. Appl. Microbiol.33, 159–163. doi: 10.1046/j.1472-765x.2001.00970.x
12
DwanC.HennessyD.IrishH.ShallooL.BuckleyF.DelabyL.et al. (2025). Evaluating the effect of white clover inclusion and dairy cow genetic merit on enteric methane emissions within rotational grazing systems. J Dairy Sci.108, 9844–9860. doi: 10.3168/jds.2024-26106
13
EdgarR. C.HaasB. J.ClementeJ. C.QuinceC.KnightR. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics27, 2194–2200. doi: 10.1093/bioinformatics/btr381
14
EganM.GalvinN.HennessyD. (2018). Incorporating white clover (Trifolium repens L.) into perennial ryegrass (Lolium perenne L.) swards receiving varying levels of nitrogen fertilizer: effects on milk and herbage production. J. Dairy Sci.101, 3412–3427. doi: 10.3168/jds.2017-13233
15
Enriquez-HidalgoD.GillilandT.DeightonM. H.O’DonovanM.HennessyD. (2014). Milk production and enteric methane emissions by dairy cows grazing fertilized perennial ryegrass pasture with or without inclusion of white clover. J. Dairy Sci.97, 1400–1412. doi: 10.3168/jds.2013-7034
16
FreguliaP.CamposM. M.DhakalR.DiasR. J. P.NevesA. L. A. (2024). Feed efficiency and enteric methane emissions indices are inconsistent with the outcomes of the rumen microbiome composition. Sci. Total Environ.949:175263. doi: 10.1016/j.scitotenv.2024.175263
17
FrescoS.BoichardD.FritzS.LefebvreR.BarbeyS.GaboritM.et al. (2023). Comparison of methane production, intensity, and yield throughout lactation in Holstein cows. J. Dairy Sci.106, 4147–4157. doi: 10.3168/jds.2022-22855
18
GreeningR. C.LeedleJ. A. (1989). Enrichment and isolation of Acetitomaculum ruminis, gen. Nov., sp. nov.: acetogenic bacteria from the bovine rumen. Arch. Microbiol.151, 399–406. doi: 10.1007/bf00416597
19
HaoY.XiaJ.WangW.WangY.CaoZ.YangH.et al. (2024). Diurnal shifts of rumen fermentation and microbial profiles revealed circadian rhythms of rumen bacteria, methanogens, and protozoa under high-grain and high-forage diets. JDS Commun5, 700–706. doi: 10.3168/jdsc.2023-0526
20
HearnC.EganM.LynchM. B.FlemingC.O’DonovanM. (2022). Seasonal variations in nutritive and botanical composition properties of multispecies grazing swards over an entire dairy grazing season. Grassl. Res.1, 221–233. doi: 10.1002/glr2.12037
21
HendersonG.CoxF.GaneshS.JonkerA.YoungW.AbeciaL.et al. (2015). Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci. Rep.5:14567. doi: 10.1038/srep14567
22
HitchT. C. A.BisdorfK.AfrizalA.RiedelT.OvermannJ.StrowigT.et al. (2022). A taxonomic note on the genus Prevotella: description of four novel genera and emended description of the genera Hallella and Xylanibacter. Syst. Appl. Microbiol.45:126354. doi: 10.1016/j.syapm.2022.126354
23
JanssenP. H. (2010). Influence of hydrogen on rumen methane formation and fermentation balances through microbial growth kinetics and fermentation thermodynamics. Anim. Feed Sci. Technol.160, 1–22. doi: 10.1016/j.anifeedsci.2010.07.002
24
KellyW. J.LeahyS. C.KamkeJ.SoniP.KoikeS.MackieR.et al. (2019). Occurrence and expression of genes encoding methyl-compound production in rumen bacteria. Animal Microb1:15. doi: 10.1186/s42523-019-0016-0
25
KingE. E.SmithR. P.St-PierreB.WrightA.-D. G. (2011). Differences in the rumen methanogen populations of lactating Jersey and Holstein dairy cows under the same diet regimen. Appl. Environ. Microbiol.77, 5682–5687. doi: 10.1128/AEM.05130-11
26
KolverE. S.de VethM. J. (2002). Prediction of ruminal pH from pasture-based diets. J. Dairy Sci.85, 1255–1266. doi: 10.3168/jds.S0022-0302(02)74190-8
27
KraatzM.WallaceR. J.SvenssonL. (2011). Olsenella umbonata sp. nov., a microaerotolerant anaerobic lactic acid bacterium from the sheep rumen and pig jejunum, and emended descriptions of Olsenella, Olsenella uli and Olsenella profusa. Int. J. Syst. Evol. Microbiol.61, 795–803. doi: 10.1099/ijs.0.022954-0
28
LahartB.BuckleyF.HerronJ.FitzgeraldR.FitzpatrickE.GalvinN.et al. (2024). Evaluating enteric methane emissions within a herd of genetically divergent grazing dairy cows. J. Dairy Sci.107, 383–397. doi: 10.3168/jds.2022-22646
29
LeeJ.WoodwardS.WaghornG., andClarkD. (2004). Methane emissions by dairy cows fed increasing proportions of white clover (Trifolium repens) in pasture, in: Proceedings of the New Zealand Grassland Association, 151–155.
30
LinH.PeddadaS. D. (2024). Multigroup analysis of compositions of microbiomes with covariate adjustments and repeated measures. Nat. Methods21, 83–91. doi: 10.1038/s41592-023-02092-7
31
LopesD. R. G.de Souza DuarteM.La ReauA. J.ChavesI. Z.de Oliveira MendesT. A.DetmannE.et al. (2021). Assessing the relationship between the rumen microbiota and feed efficiency in Nellore steers. J Anim Sci Biotechnol12:79. doi: 10.1186/s40104-021-00599-7
32
LyonsT.BielakA.DoyleE.KuhlaB. (2018). Variations in methane yield and microbial community profiles in the rumen of dairy cows as they pass through stages of first lactation. J. Dairy Sci.101, 5102–5114. doi: 10.3168/jds.2017-14200
33
MacdonaldK.PennoJ.NicholasP.LileJ.CoulterM.LancasterJ. (2001). Farm systems-impact of stocking rate on dairy farm efficiency, in: Proceedings of the conference-New Zealand GRASSLAND association, 223–228.
34
Mackie RoderickI.Aminov RustamI.HuW.Klieve AtholV.OuwerkerkD.Sundset MonicaA.et al. (2003). Ecology of uncultivated Oscillospira species in the rumen of cattle, sheep, and reindeer as assessed by microscopy and molecular approaches. Appl. Environ. Microbiol.69, 6808–6815. doi: 10.1128/AEM.69.11.6808-6815.2003
35
MagočT.SalzbergS. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics27, 2957–2963. doi: 10.1093/bioinformatics/btr507
36
McAuliffeS.MeeJ. F.LewisE.GalvinN.HennessyD. (2022). Feeding system effects on dairy cow rumen function and Milk production. Animals12:523. doi: 10.3390/ani12040523
37
McMurdieP. J.HolmesS. (2013). Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One8:e61217. doi: 10.1371/journal.pone.0061217
38
MonteiroH. F.FaciolaA. P. (2020). Ruminal acidosis, bacterial changes, and lipopolysaccharides. J. Anim. Sci.98:248. doi: 10.1093/jas/skaa248
39
MorotomiM.NagaiF.WatanabeY. (2012). Description of Christensenella minuta gen. Nov., sp. nov., isolated from human faeces, which forms a distinct branch in the order Clostridiales, and proposal of Christensenellaceae fam. Nov. Int. J. Syst. Evol. Microbiol.62, 144–149. doi: 10.1099/ijs.0.026989-0
40
MottetA.TeillardF.BoettcherP.De’ BesiG.BesbesB. (2018). Review: Domestic herbivores and food security: current contribution, trends and challenges for a sustainable development. Animal12, s188–s198. doi: 10.1017/S1751731118002215
41
MurrayÁ.DelabyL.GillilandT. J.McCarthyB. (2023). Nitrogen offset potential in a multi-year farmlet-scale study; milk and herbage production from grazed perennial ryegrass/white clover swards. J. Dairy Sci.107, 2129–2142. doi: 10.3168/jds.2023-23822
42
NoelS. J.AttwoodG. T.RakonjacJ.MoonC. D.WaghornG. C.JanssenP. H. (2017). Seasonal changes in the digesta-adherent rumen bacterial communities of dairy cattle grazing pasture. PLOS ONE12, e0173819. doi: 10.1371/journal.pone.0173819
43
NoelS. J.OlijhoekD. W.McleanF.LøvendahlP.LundP.HøjbergO. (2019). Rumen and Fecal microbial community structure of Holstein and Jersey dairy cows as affected by breed, diet, and residual feed intake. Animals9:498. doi: 10.3390/ani9080498
44
O’GradyL.DohertyM. L.MulliganF. J. (2008). Subacute ruminal acidosis (SARA) in grazing Irish dairy cows. Vet. J.176, 44–49. doi: 10.1016/j.tvjl.2007.12.017
45
O’NeillB. F.DeightonM. H.O’LoughlinB. M.MulliganF. J.BolandT. M.O’DonovanM.et al. (2011). Effects of a perennial ryegrass diet or total mixed ration diet offered to spring-calving Holstein-Friesian dairy cows on methane emissions, dry matter intake, and milk production. J. Dairy Sci.94, 1941–1951. doi: 10.3168/jds.2010-3361
46
O’BrienD.MoranB.ShallooL. (2018). A national methodology to quantify the diet of grazing dairy cows. J. Dairy Sci.101, 8595–8604. doi: 10.3168/jds.2017-13604
47
OksanenJ.KindtR.LegendreP.O’HaraB.StevensM. H. H.OksanenM. J.et al. (2007). The vegan package. Commun Ecol Package10:719. Available at: https://www.researchgate.net/publication/228975085_The_Vegan_Package
48
OlijhoekD. W.HellwingA. L. F.NoelS. J.LundP.LarsenM.WeisbjergM. R.et al. (2022). Feeding up to 91% concentrate to Holstein and Jersey dairy cows: effects on enteric methane emission, rumen fermentation and bacterial community, digestibility, production, and feeding behavior. J. Dairy Sci.105, 9523–9541. doi: 10.3168/jds.2021-21676
49
O’SullivanM.DillonP.O’SullivanK.PierceK. M.GalvinN.EganM.et al. (2019). Intake, efficiency, and feeding behavior characteristics of Holstein-Friesian cows of divergent economic breeding index evaluated under contrasting pasture-based feeding treatments. J. Dairy Sci.102, 8234–8246. doi: 10.3168/jds.2019-16371
50
PalevichN.Kelly WilliamJ.Leahy SineadC.DenmanS.AltermannE.RakonjacJ.et al. (2019). Comparative genomics of rumen Butyrivibrio spp. uncovers a continuum of polysaccharide-degrading capabilities. Appl. Environ. Microbiol.86:e01993-01919. doi: 10.1128/AEM.01993-19
51
RiusA. G.KittelmannS.MacdonaldK. A.WaghornG. C.JanssenP. H.SikkemaE. (2012). Nitrogen metabolism and rumen microbial enumeration in lactating cows with divergent residual feed intake fed high-digestibility pasture. J. Dairy Sci.95, 5024–5034. doi: 10.3168/jds.2012-5392
52
RobertsonL. J.WaghornG. C. (2002). Dairy industry perspectives of methane emissions and production from cattle fed pasture or total mixed rations in New Zealand, in: Proceedings of the New Zealand Society of Animal Production. (Palmerston North: New Zealand Society of Animal Production).
53
RocheJ. R.BerryD. P.BryantA. M.BurkeC. R.ButlerS. T.DillonP. G.et al. (2017). A 100-year review: a century of change in temperate grazing dairy systems. J. Dairy Sci.100, 10189–10233. doi: 10.3168/jds.2017-13182
54
RoeheR.DewhurstR. J.DuthieC.-A.RookeJ. A.McKainN.RossD. W.et al. (2016). Bovine Host Genetic Variation Influences Rumen Microbial Methane Production with Best Selection Criterion for Low Methane Emitting and Efficiently Feed Converting Hosts Based on Metagenomic Gene Abundance. PLOS Genetics12:e1005846. doi: 10.1371/journal.pgen.1005846
55
RussellJ. B.HespellR. B. (1981). Microbial Rumen Fermentation. J. dairy sci. 64, 1153–1169. doi: 10.3168/jds.S0022-0302(81)82694-X
56
RutterS. M.OrrR. J.YarrowN. H.ChampionR. A. (2004). Dietary preference of dairy cows grazing ryegrass and white clover. J. Dairy Sci.87, 1317–1324. doi: 10.3168/jds.S0022-0302(04)73281-6
57
SmithP. E.Enriquez-HidalgoD.HennessyD.McCabeM. S.KennyD. A.KellyA. K.et al. (2020). Sward type alters the relative abundance of members of the rumen microbial ecosystem in dairy cows. Sci. Rep.10:9317. doi: 10.1038/s41598-020-66028-3
58
StegA.van StraalenW. M.HindleV. A.WensinkW. A.DooperF. M. H.SchilsR. L. M. (1994). Rumen degradation and intestinal digestion of grass and clover at two maturity levels during the season in dairy cows. Grass Forage Sci.49, 378–390. doi: 10.1111/j.1365-2494.1994.tb02014.x
59
StepanchenkoN.StefenoniH.HennessyM.NagarajuI.WassonD. E.CuevaS. F.et al. (2023). Microbial composition, rumen fermentation parameters, enteric methane emissions, and lactational performance of phenotypically high and low methane-emitting dairy cows. J. Dairy Sci.106, 6146–6170. doi: 10.3168/jds.2022-23190
60
St-PierreB.CersosimoL. M.IshaqS. L.WrightA.-D. G. (2015). Toward the identification of methanogenic archaeal groups as targets of methane mitigation in livestock animalsr. Front. Microbiol.6:776. doi: 10.3389/fmicb.2015.00776
61
TapioI.SnellingT. J.StrozziF.WallaceR. J. (2017). The ruminal microbiome associated with methane emissions from ruminant livestock. J Anim Sci Biotechnol8:7. doi: 10.1186/s40104-017-0141-0
62
TaweelH. Z.TasB. M.SmitH. J.ElgersmaA.DijkstraJ.TammingaS. (2005). Effects of feeding perennial ryegrass with an elevated concentration of water-soluble carbohydrates on intake, rumen function and performance of dairy cows. Anim. Feed Sci. Technol.121, 243–256. doi: 10.1016/j.anifeedsci.2005.02.024
63
TongJ.ZhangH.YangD.ZhangY.XiongB.JiangL. (2018). Illumina sequencing analysis of the ruminal microbiota in high-yield and low-yield lactating dairy cows. PLoS One13:e0198225. doi: 10.1371/journal.pone.0198225
64
UlyattM. J. (1971). Studies on the causes of the differences in pasture quality between perennial ryegrass, short-rotation ryegrass, and white clover. N. Z. J. Agric. Res.14, 352–367. doi: 10.1080/00288233.1971.10427100
65
Van KesselJ. A. S.RussellJ. B. (1996). The effect of pH on ruminal methanogenesis. FEMS Microbiol. Ecol.20, 205–210. doi: 10.1111/j.1574-6941.1996.tb00319.x
66
Van SoestP. J. (1994). Nutritional ecology of the ruminant. Ithaca, New York, USA: Cornell University Press.
67
VeerkampR. F.DillonP.KellyE.CromieA. R.GroenA. F. (2002). Dairy cattle breeding objectives combining yield, survival and calving interval for pasture-based systems in Ireland under different milk quota scenarios. Livest. Prod. Sci.76, 137–151. doi: 10.1016/S0301-6226(02)00006-4
68
WangM.SunX. Z.JanssenP. H.TangS. X.TanZ. L. (2014). Responses of methane production and fermentation pathways to the increased dissolved hydrogen concentration generated by eight substrates in in vitro ruminal cultures. Anim. Feed Sci. Technol.194, 1–11. doi: 10.1016/j.anifeedsci.2014.04.012
69
WangY.JanssenP. H.LynchT. A.BruntB. V.PachecoD. (2016). A mechanistic model of hydrogen–methanogen dynamics in the rumen. J. Theor. Biol.393, 75–81. doi: 10.1016/j.jtbi.2015.12.026
70
WoodmartinS.SmithP. E.CreightonP.BolandT. M.DunneE.McGovernF. M. (2024). Sward type alters enteric methane emissions, nitrogen output and the relative abundance of the rumen microbial ecosystem in sheep. J. Anim. Sci.102:256. doi: 10.1093/jas/skae256
71
YuZ.MorrisonM. (2004). Improved extraction of PCR-quality community DNA from digesta and fecal samples. BioTechniques36, 808–812. doi: 10.2144/04365st04
72
ZhouM.WuJ.WuL.SunX.ChenC.HuangL. (2024). The utilization of N-acetylgalactosamine and its effect on the metabolism of amino acids in Erysipelotrichaceae strain. BMC Microbiol.24:397. doi: 10.1186/s12866-024-03505-z
73
ZhuZ.DiffordG. F.NoelS. J.LassenJ.LøvendahlP.HøjbergO. (2021). Stability assessment of the rumen bacterial and archaeal communities in dairy cows within a single lactation and its association with host phenotype. Front. Microbiol.12:636223. doi: 10.3389/fmicb.2021.636223
Summary
Keywords
dairy, rumen microbiome, pasture, season, enteric methane, clover, breeding, rumen fermentation characteristics
Citation
Dwan C, Das A, O’Toole PW, O’Callaghan TF, Meehan D, Hennessy D, Irish H, Buckley F and Lahart B (2025) Evaluating the rumen microbial community of genetically divergent spring-calving dairy cows grazing grass-only or grass-clover swards at different stages of the grazing season. Front. Microbiol. 16:1642486. doi: 10.3389/fmicb.2025.1642486
Received
06 June 2025
Accepted
01 October 2025
Published
06 November 2025
Volume
16 - 2025
Edited by
Anusorn Cherdthong, Khon Kaen University, Thailand
Reviewed by
Biswarup Sen, Independent Researcher, Bengaluru, India; Chanadol Supapong, Rajamangala University of Technology Isan, Thailand
Updates
Copyright
© 2025 Dwan, Das, O’Toole, O’Callaghan, Meehan, Hennessy, Irish, Buckley and Lahart.
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: Charles Dwan, charles.dwan@teagasc.ie
Disclaimer
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.