Skip to main content

Fecal microbiome of horses transitioning between warm-season and cool-season grass pasture within integrated rotational grazing systems

Abstract

Background

Diet is a key driver of equine hindgut microbial community structure and composition. The aim of this study was to characterize shifts in the fecal microbiota of grazing horses during transitions between forage types within integrated warm- (WSG) and cool-season grass (CSG) rotational grazing systems (IRS). Eight mares were randomly assigned to two IRS containing mixed cool-season grass and one of two warm-season grasses: bermudagrass [Cynodon dactylon (L.) Pers.] or crabgrass [Digitaria sanguinalis (L.) Scop.]. Fecal samples were collected during transitions from CSG to WSG pasture sections (C–W) and WSG to CSG (W–C) on days 0, 2, 4, and 6 following pasture rotation and compared using 16S rRNA gene sequencing.

Results

Regardless of IRS or transition (C–W vs. W–C), species richness was greater on day 4 and 6 in comparison to day 0 (P < 0.05). Evenness, however, did not differ by day. Weighted UniFrac also did not differ by day, and the most influential factor impacting β-diversity was the individual horse (R2 ≥ 0.24; P = 0.0001). Random forest modeling was unable to accurately predict days within C–W and W–C, but could predict the individual horse based on microbial composition (accuracy: 0.92 ± 0.05). Only three differentially abundant bacterial co-abundance groups (BCG) were identified across days within all C–W and W–C for both IRS (W ≥ 126). The BCG differing by day for all transitions included amplicon sequence variants (ASV) assigned to bacterial groups with known fibrolytic and butyrate-producing functions including members of Lachnospiraceae, Clostridium sensu stricto 1, Anaerovorax the NK4A214 group of Oscillospiraceae, and Sarcina maxima. In comparison, 38 BCG were identified as differentially abundant by horse (W ≥ 704). The ASV in these groups were most commonly assigned to genera associated with degradation of structural carbohydrates included Rikenellaceae RC9 gut group, Treponema, Christensenellaceae R-7 group, and the NK4A214 group of Oscillospiraceae. Fecal pH also did not differ by day.

Conclusions

Overall, these results demonstrated a strong influence of individual horse on the fecal microbial community, particularly on the specific composition of fiber-degraders. The equine fecal microbiota were largely stable across transitions between forages within IRS suggesting that the equine gut microbiota adjusted at the individual level to the subtle dietary changes imposed by these transitions. This adaptive capacity indicates that horses can be managed in IRS without inducing gastrointestinal dysfunction.

Background

Traditional cool-season grass rotational equine pastures in temperate regions of the United States are most productive during the spring, early summer, and fall months, but undergo a period of low forage productivity during hot, dry summer months [1]. Conversely, due to differences in photosynthetic mechanisms, warm-season grasses are most vigorous during this same “summer slump” period, when cool-season grasses are semi-dormant [2,3,4]. Thus, an integrated rotational grazing system that incorporates both warm- and cool-season grasses can increase pasture yield throughout this period of the growing season [5,6,7].

Production of prussic acid associated with many warm-season annual species and cold-sensitivity of perennials limit options for integrating warm-season grasses into equine grazing systems. However, some annuals such as teff and improved forage varieties of crabgrass may be grazed by horses with a lesser concern for forage-related disorders (i.e. prussic acid or nitrate toxicity, etc.) [8, 9]. Additionally, cultivation of cold-tolerant bermudagrass varieties has allowed this perennial to be grown in the transition zone corresponding to United States Department of Agriculture (USDA) plant hardiness zones 5–7 [10, 11]. These warm-season annuals and perennials, therefore, could be utilized in equine pastures to increase forage availability during summer months.

While an integrated warm- and cool-season grass rotational grazing management strategy has demonstrated benefits for pasture yield, potential impacts on equine gastrointestinal health have not been evaluated. The equine hindgut microbial ecosystem is particularly sensitive to sudden dietary change [12, 13], as evidenced by the well-documented role of the gut microbiota in carbohydrate induction models of laminitis [14,15,16]. Additionally, in epidemiological studies, recent dietary change, including changes in hay or forage, has been consistently associated with increased risk of colic [13, 17, 18]. While the role of the hindgut microbiome in the etiology of colic has not been fully elucidated, fluctuations in the gut microbiota preceding and following episodes of colic have been documented [19, 20], and Stewart et al. [21] found differences in the hindgut microbial communities of horses admitted to a veterinary hospital for treatment of colic versus horses admitted for other elective procedures unrelated to gastrointestinal disease. Laminitis is the leading cause of foot lameness in horses [22], and colic is the leading cause of veterinary emergencies and mortality in adult horses [23, 24]. Characterizing shifts in gut microbial community structure that could occur when horses are transitioned between warm- and cool-season grasses and how the equine hindgut microbiota adapts to differing forage types is necessary to fully understand the impact of integrated warm- and cool-season grass rotational grazing management on horse health.

However, a decreased risk of colic has been reported for horses with pasture access [13], potentially suggesting greater stability of the equine hindgut microbiome in grazing horses. Furthermore, the microbial communities of grazing horses have been found to be more diverse than in horses fed mixed diets of conserved forage and concentrate [25]. This research seeks to understand the extent to which this diversity confers enhanced resilience through the transition between forage types for horses managed in integrated rotational grazing systems.

Dietary nutrient composition is a key driver of gut microbial community structure in humans and across animal species including the horse [26,27,28]. In addition to differences in photosynthetic processes, warm- and cool-season grasses also differ in mechanisms for storage of soluble carbohydrates such as sugars, starches, and fructans [2, 29]. Differences in soluble carbohydrate storage contribute to varying nutrient composition between these two forage types, with non-structural carbohydrate (NSC; NSC = fructans + sugars + starch) concentrations typically greater in cool-season grasses in comparison to warm-season grasses [5, 29, 30]. Prior studies have reported higher fiber and lower protein and digestibility in warm- vs. cool-season grasses [3, 4, 31]; however, the hindgut microbiomes of horses grazing cool- vs. warm-season pasture grasses have not been previously characterized. Recommendations for integrating warm- and cool-season grasses in horse pasture systems should be informed by an understanding of potential impacts of this management practice on the hindgut microbiome. Therefore, the objectives of this study were to characterize shifts in the fecal microbiota of grazing horses during transitions between forage types within integrated warm- and cool-season grass rotational grazing systems.

Results

Samples and initial 16s rRNA gene sequence analysis

Fecal samples were collected from eight horses assigned to one of two integrated warm-and cool-season rotational grazing systems containing either the warm-season perennial Wrangler bermudagrass [(BER); Cynodon dactylon (L.) Pers.; Johnston Seed Company, Enid, OK] or the warm-season annual Quick-N-Big crabgrass [(CRB); Digitaria sanguinalis (L.) Scop.; Dalrymple Farms, Thomas, OK], with systems denoted as BRS and CRS, respectively. Samples were collected during transitions from cool-season grass to warm-season grass (C–W) and warm-season grass to cool-season grass (W–C). Samples were collected on the day of rotation (D0) and subsequently on D2, D4, and D6 for both C–W and W–C. Prior to each transition, horses had been grazing the respective forage (cool-season grass for C–W and warm-season grass for W–C) for 17–21 days preceding the rotation. The duration of grazing prior to transition was selected to represent the longest duration that horses were likely to continuously consume an individual forage type within the context of integrated warm- and cool-season grass rotational grazing management.

In the 64 fecal samples analyzed for this study, there were a total of 1,289,921 reads prior to initial quality and chimera filtering of paired reads in Qiime 2 [32]. A total of 496,679 reads passed the initial filtering steps. The initial minimum frequency count per sample was 4617 and the maximum was 12,340, with a mean frequency of 7761 and a median frequency of 7750. Following additional filtering to remove low abundance features, 395,520 reads remained, with a minimum frequency of 3711 and a maximum of 10,541. The mean frequency of the final filtered sequence data was 6180, and the median frequency was 6104. The total number of distinct amplicon sequence variants (ASV) generated in this final dataset was 1030. Taxonomy of all ASV included in the final dataset is detailed in Additional File 1.

Diversity analyses

For each transition, α-diversity metrics by day, including the Shannon Diversity Index, Faith’s Phylogenetic Diversity, Pielou’s Evenness, and Observed ASVs are shown in Fig. 1a–d. The microbial species richness (Observed ASVs) differed by day and transition (C–W vs. W–C; mixed model AOVA with Tukey’s post hoc adjustment; P < 0.02) but did not differ by grazing system. Regardless of system or transition, species richness was greater on D4 and D6 in comparison to D0 (P < 0.05). Overall, species richness was greater in W–C than in C–W transitions (P = 0.0008). Evenness, however, did not differ by day or transition, but there was a trend for difference by system, with evenness greater in CRS vs. BRS (P < 0.10). The Shannon Diversity Index did not vary by day, transition, or grazing system. Faith’s Phylogenetic Diversity differed by transition, with greater diversity in W–C vs. C–W (P = 0.0015). There was also a trend for a difference in Faith’s Diversity by day, with greater diversity on D6 in comparison to D0 (P = 0.07). There were no interactions between main effects terms for any of the α-diversity metrics evaluated.

Fig. 1
figure 1

Fecal microbiota α-diversity following transitions between cool-season and warm-season grass. Metrics including a Observed ASVs (richness) b Pielou’s Evenness, c Shannon Diversity Index, and d Faith’s Phylogenetic Diversity were analyzed by mixed model ANOVA with Tukey’s post hoc adjustment for transitions from cool-season grass to warm-season grass (C–W) and warm-season grass to cool-season grass (W–C). There was an effect of day and transition for Observed ASVs (P < 0.02). There was an effect of transition (P = 0.0015), but not for Faith’s Phylogenetic Diversity. Shannon Diversity and Pielou’s Evenness did not differ by transition or day, and there were no significant interactions between any of the main effects for any of the metrics analyzed

Principle coordinate analysis of β-diversity metrics did not reveal distinct clustering by day within each transition (Fig. 2a–d). This was confirmed with statistical analysis by Adonis PERMANOVA, which did not find a significant effect of day for Weighted UniFrac (R2 = 0.013; P = 0.13). There was, however, as significant effect of day for Unweighted UniFrac (R2 = 0.013; P = 0.039), as well as for both Bray–Curtis Dissimilarity (R2 = 0.014; P = 0.014) and Jaccard Index (R2 = 0.015; P = 0.011). There was also a significant effect of transition (C–W vs. W–C) for Weighted UniFrac (R2 = 0.032), Unweighted UniFrac (R2 = 0.052), Bray–Curtis (R2 = 0.037), and Jaccard (R2 = 0.042; P = 0.001). The influence of grazing system (Weighted UniFrac: R2 = 0.16; Unweighted UniFrac: R2 = 0.10; Bray–Curtis: R2 = 0.11; Jaccard: R2 = 0.081; P = 0.001) was stronger than the influence of either day or transition. However, the most influential factor was the individual horse, with strong R2 values for all β-diversity metrics (Weighted UniFrac: R2 = 0.24; Unweighted UniFrac: R2 = 0.31; Bray–Curtis: R2 = 0.34; Jaccard: R2 = 0.29; P = 0.001), which was also reflected in visual clustering in the PCoA plots. There were also significant interactions between main terms, which are detailed in full in Table 1. Subsequent application of PERMDISP (permutational analysis of dispersion) confirmed that differences in the main effects of day, transition, grazing system, and horse were not due to differences of variance or dispersion within groups.

Fig. 2
figure 2

Fecal microbiota β-diversity following transitions between cool-season and warm-season grass. Metrics including a Weighted UniFrac, b Unweighted UniFrac, c Bray Curtis Dissimilarity, and d Jaccard Index were analyzed for transitions from cool-season to warm-season grass (C–W) and warm-season to cool-season grass (W–C) in both the bermudagrass integrated rotational grazing system (BRS) and the crabgrass integrated system (CRS). Individual horses are designated with two-letter abbreviations and arranged such that horses on the darker end of the color spectrum grazed in BRS (n = 4) and horses on the lighter end of the color spectrum grazed in CRS (n = 4). Analysis by PERMANOVA with the Adonis action in Qiime 2 (v.2020.8) found that the influence of individual horse was the most influential factor shaping β-diversity across all metrics (R2 ≥ 24; P = 0.001)

Table 1 Results of Adonis1 multivariate PERMANOVA analysis of β-diversity showing influence of main effects and interactions

Differential abundance

Application of Sparce Cooccurrence Network Investigation for Compositional Data (SCNIC) identified 149 distinct bacterial co-abundance groups (BCG), with 621 individual ASV remaining ungrouped. Random forest classification models were able to accurately predict grazing system (0.95 ± 0.06), transition (0.97 ± 0.06), and horse (0.92 ± 0.05) based on microbial composition. However, day could not be accurately predicted in either the C–W (0.25 ± 0.13) or the W–C transitions (0.22 ± 0.18). Accordingly, Analysis of Composition of Microbes (ANCOM) found three BCG/ungrouped ASV that differed by transition (W ≥ 264), seven BCG that differed by grazing system (W ≥ 701), and 42 BCG and ungrouped ASV that differed across horses (W ≥ 704). Conversely, there was only one differentially abundant BCG across days within each of the transitions (C–W and W–C) for horses in CRS (W ≥ 126). There was also only one BCG that differed across days within the W–C transition for horses in BRS (≥ 299), while no BCG or ungrouped ASV differed across days within the C–W transition for BRS based on analysis by ANCOM.

The BCG and ungrouped ASV identified as differing by grazing system, transition, or day as well as taxonomic classifications are presented in Tables 2, 3 and 4. Features differing across horses, including taxonomy, are detailed in Additional File 2. The ungrouped ASV and ASV within BCG that differed by transition were assigned to taxa including the genera Christensenellaceae R-7 group and Rikenellaceae RC9 gut group as well as ASV assigned only to the family level for Lachnospiraceae, and the class level for Bacilli and Clostridia. These ungrouped ASV and BCG were all more abundant in the C–W transition, but their summed relative abundance comprised < 1% of total microbial community on average across all samples. The BCG that differed between grazing systems also represented a small portion of the total microbial community at < 7%. Of these BCG, five were more abundant in BRS. The BCG that were more abundant in BRS included ASV members assigned to the genera Fibrobacter, Treponema, Christensellaceae R-7 group, Prevotellaceae 6a6A1 group, Methanocorpusculum, the NK4A214 group within Oscillospiraceae, and the UCG-010 genus and family within Oscillospirales. Additional ASV within these groups were mapped only to the order level for Bacteroidales and the class level for Bacteroidia. Two BCG were more abundant in CRS, with ASV members from Denitrobacterium detoxificans and the genera Desulfovibrio and Phascolarctobacterium.

Table 2 Differentially abundant1 bacterial co-abundance groups (BCG)2 and ungrouped amplicon sequence variants (ASV) across transitions3
Table 3 Differentially abundant1 bacterial co-abundance groups (BCG)2 and ungrouped amplicon sequence variants (ASV) across grazing systems3
Table 4 Differentially abundant1 bacterial co-abundance groups (BCG)2 and ungrouped amplicon sequence variants (ASV) across days3

Similarly, the three BCG that differed by day within transitions and systems comprised a small percentage (< 2%) of the total microbial community. The BCG that differed by day for the C–W transition within CRS (BCG_114) included ASV mapped to the Lachnospiraceae family and the NK4A214 group within Oscillospiraceae, with this BCG increasing across days within C–W (Fig. 3a). The BCG that differed by day (increasing) for the W–C transition within CRS (BCG_124; Fig. 3b) contained two ASV assigned to Sarcina maxima, while the BCG that differed by day (decreasing) for the W–C transition within BRS (BCG_10; Fig. 3c) included member ASV assigned to Clostridium sensu stricto 1, Anaerovorax, and Lachnospiraceae.

Fig. 3
figure 3

Differentially abundant bacterial co-abundance groups (BCG) following transitions between cool-season and warm-season grass. Relative abundances are shown across days within transitions between a cool-season grass to warm-season grass pasture and b warm-season to cool-season grass within the crabgrass integrated rotational grazing system and the c warm-season to cool-season grass transition within the bermudagrass integrated rotational grazing system. The BCG were identified as differentially abundant using Analysis of Composition of Microbes (ANCOM) in Qiime 2 (v.2020.8) (W ≥ 126)

In comparison to the relatively small number of BCG for which abundance changed in response to transition, system, and day, 38 BCG and 4 ungrouped ASV were identified as differentially abundant by horse. When averaged across all samples, the BCG and ungrouped ASV that differed by horse comprised > 32% of the total microbial community. Abundance of BCG and ungrouped ASV identified as differing by horse are shown in Fig. 4. The genera to which these ungrouped ASV and ASV members of the BCG were most commonly assigned included Rikenellaceae RC9 gut group (18 ASV), Treponema (10 ASV), Christensenellaceae R-7 group (10 ASV), and the NK4A214 group of Oscillospiraceae (8 ASV). Shifts in relative abundance of the twelve most abundant BCG that differed by horse are shown in Additional File 3.

Fig. 4
figure 4

Differentially abundant bacterial co-abundance groups (BCG) and ungrouped amplicon sequenc variants (ASV) by horse. Features (BCG and ASV) were dentified as differentially abundant (W ≥ 704) across horses using Analysis of Composition of Microbes (ANCOM) in Qiime 2 (v.2020.8). Individual horses are designated with two-letter abbreviations and arranged such that each horse is represented by an individual facet grid, with horses on the left-hand side assigned to the bermudagrass integrated rotational grazing system (n = 4) and horses on the right-hand side assigned to the crabgrass integrated system (n = 4). Sample IDs on the x-axis are arranged such that the first four samples within each horse represent transitions from cool-season grass to warm-season grass and the final four samples within each horse represent transitions from warm-season grass to cool-season grass. Within transitions, samples are arranged from D0 through D6

Analysis of microbial composition at the genus-level confirmed results presented for BCG analysis. Grazing system, transition, and horse could be predicted based on microbial composition through random forest classification with accuracies ≥ 0.80, while day could not be accurately predicted for either transition (accuracy ≤ 0.22). Additionally, ANCOM identified a limited number of genera groups and/or ungrouped genera as differentially abundant by grazing system, transition, and day (within transition for each system) in comparison to genera groups that differed by horse (see Additional File 4). However, predictive accuracy of the random forest classifiers was over 10% lower than for the BCG analysis and the percentage of the total microbial communities represented by the differential abundant genera was also lower (see analytical comparison in Additional File 5).

Fecal pH

Fecal pH across the days of each transition is shown in Fig. 5. Consistent with above results for microbiome data, fecal pH differed by transition and grazing system (mixed model AOVA with Tukey’s post hoc adjustment; P ≤ 0.02), but did not differ by day, nor were there any significant interactions between the main effects. Fecal pH was lower in the C–W transition (7.09 ± 0.10) than in the W–C transition (7.42 ± 0.10; P = 0.02). Fecal pH was also lower in horses within BRS (6.97 ± 0.13) in comparison to CRS (7.55 ± 0.10; P = 0.02). Overall, the lowest recorded fecal pH value was 6.0 for two horses in BRS, one on D2 and the other on D4 of the W–C transition. The highest recorded value for fecal pH was 8.0 for three horses within CRS on various days of both the C–W and the W–C transitions.

Fig. 5
figure 5

Fecal pH following transitions between cool-season and warm-season grass. Fecal pH of horses across six days following transitions from cool-season to warm-season grass (C–W) and warm-season to cool-season grass (W–C) in both the bermudagrass integrated rotational grazing system and the crabgrass integrated system. There was an effect of transition and grazing system (mixed model ANOVA with Tukey’s post hoc adjustment; P ≤ 0.02), but pH did not differ by day. There were no significant interactions between main effects

Discussion

Results of this study suggest that the equine microbiome is largely stable during transition between forages within an integrated warm- and cool-season grass rotational grazing system. While species richness did increase by D4 of the transitions and there was a trend for increased phylogenetic diversity (Faith’s Diversity) by D6, evenness and the Shannon Index did not differ over the course of the transitions. Multivariate β-diversity analysis also indicated minimal influence of day in comparison to other variables such as the specific transition between forages (C–W vs. W–C), grazing system (BRS vs. CRS), and individual horse. Additionally, random forest classification was unable to accurately predict days within the C–W and W–C transitions. This was reinforced by the identification of only 3 BCG or ungrouped ASV as differentially abundant across days within specific transitions, with these bacteria representing less than 2% of the total fecal microbial community. Furthermore, fecal pH, as a marker of microbial community function, did not differ by day within transitions between pasture forages. Stability of the microbiota during transitions between pasture forage types indicates that the equine hindgut microbiota is capable of adapting to subtle shifts in nutrient composition between warm- and cool-season grasses.

The lack of differences in these variables across days within transitions may be reflective of the fact that there were only minimal differences in nutritional composition between warm- and cool-season grasses in either system (mean nutrient composition of pasture forages is shown in Table 5). The most pronounced difference in nutrient concentration was the 10% greater crude protein in CRB compared to cool-season grass for the C–W transition in CRS. Water-soluble carbohydrate (WSC; WSC = fructans + sugars) in cool-season grasses were twice that of BER (in BRS) and CRB (in CRS) during the C–W transition. However, the magnitude of these differences was minimal, with WSC concentrations of 5.5 – 7% for cool-season grasses and 3 – 3.5% for BER and CRB. It should be noted that forage processing protocols following sample collection (i.e. forced-air or radiant heat vs. freeze drying) can impact measured plant soluble carbohydrates and that processing samples by forced-air drying as was conducted in the current study would result in lower NSC concentrations due to the continuation of plant respiration through the drying process [33]. However, all samples of all forages in the present study were processed similarly, so this impact of processing method would have been similar in all samples regardless of forage type.

Table 5 Nutrient composition1 of pasture forages

While minimal changes in microbial species composition were observed across days in the current study, differences did occur in BCG containing ASV from the Lachnospiraceae family as well as the NK4A214 group within Oscillospiraceae, Clostridium sensu stricto 1, Anaerovorax, and Sarcina. Differences in Sarcina across the W–C transition in CRS were also reflected in the genus-level analysis in addition to fluctuation in Bacillus for the W–C in BRS. Garber et al. [34] also found Sarcina enriched in horses transitioning to cool-season grass from a hay diet. These results indicate that taxa represented by these ASV, and the Sarcina and Bacillus genera more broadly, may be more sensitive to subtle changes in pasture forage type and nutrient composition.

Few studies have investigated longitudinal shifts in the hindgut or fecal microbiome of horses during transitions between all-forage diets [34,35,36,37], with prior studies primarily focused on abrupt dietary changes between concentrates [38] or abrupt inclusion/elimination of supplemental concentrates in horses maintained on hay or pasture [12, 25, 39, 40]. Garber et al., [34] reported no significant changes in α-diversity (either species richness or evenness) over a 14-d transition between hay and pasture but did identify broad phylum-level changes in microbial composition as well as lower taxa enriched on specific days of the transition periods. In contrast, the current study found changes in species richness, but limited changes in microbial species composition. These differences could be attributable to the fact that Garber et al., [34] evaluated hay and pasture diets, while the current study assessed transition between two pasture forages. To the authors' knowledge, this is the first study to report on transitions between pasture forages or cool- vs. warm-season grasses. The minimal change in microbial species composition in this study is similar to results of previous studies in horses transitioned either between hays [36] or from hay to either haylage or silage of similar botanical composition [41] or between silages varying in crude protein [35].

However, differences in experimental design (i.e. total duration of treatment or adaptation periods, sampling timepoints, etc.) and analytical approaches (culture-based vs. culture-independent, methods of statistical analysis, etc.) make inter-study comparisons difficult. The current study utilized a guild-based analytical approach, grouping individual ASV into BCG prior to random forest classification or analysis of differential abundance, rather than more conventional taxon-based analyses. Given the genetic variation between bacteria within a given taxa (even strains within a single species), members of taxonomic groups may not be functionally homologous and may respond differentially to experimental treatments such as diet or environmental variables [42,43,44,45]. Results of this study support this suggestion, as distinct ASV assigned to the same taxonomic lineage were identified in BCG enriched in different transitions, grazing systems, and/or horses. Similarly, Mach et al. [46] reported differential responses of ASV with identical taxonomic assignments in horses transitioned from pasture to a stabled environment for performance training. A number of other recent equine microbiome studies have also evaluated microbial composition at the level of individual ASV [47,48,49]. Strategies similar to the analyses employed in the current study (grouping bacteria by co-abundance instead of taxonomic rank for purposes of dimensionality reduction) have been previously implemented in studies of the human [50] and mouse gut microbiomes [43, 51, 52]. The use of this BCG-based approach in the current study is further supported by the greater predictive accuracy of random forest classifiers and the greater percentage of the total microbial community captured by differential abundance analysis in comparison to the concurrent taxon-based analysis at the genus level.

Another primary finding of this study was the strong influence of the individual horse on microbial diversity and composition. Analyses of β-diversity metrics by PERMANOVA with Adonis yielded R2 values from 0.24 to 0.34, which were 1.5–3.6 times the R2 values for the next greatest variable, depending upon the metric analyzed. Application of a random forest classification model indicated that individual horses could be predicted based on microbial community composition with high accuracy, and BCG (and ungrouped ASV) comprising > 32% of the total microbial community composition was identified as differing by horse. Substantial inter-horse variation of the hindgut microbiota has been previously documented in ponies [34, 53] and horses [54,55,56]. Gomez et al. [48] demonstrated a differential response of individual horses to dietary treatments consisting of a conventional or reduced-lignin alfalfa hay. This individualized response to dietary intervention has also been reported in humans [57] and mice [58].

In the current study, differential abundance analysis found that the most common genera represented in BCG (as well as ungrouped ASV) that differed across horses were Rikenellaceae RC9 gut group, Treponema, Christensenellaceae R-7 group, and the NK4A214 group of Oscillospiraceae. These genera have all been previously identified in the equine fecal microbiota of forage-fed horses [48, 59, 60] and have been associated with degradation of a diverse array of structural carbohydrates [61,62,63,64,65] as well as production of butyrate and other short-chain fatty acids [66,67,68]. Interestingly, Christensenellaceae has been identified as a highly heritable taxa in the gut of humans and other animal species [61, 69,70,71], which may explain the inter-individual variability in ASV assigned to this group. Inter-individual variation has also been found in Treponema, with Gomez et al. [48] reporting that ASV within this genus displayed differing responses to forage lignin concentration across individual horses. Rikenellaceae RC9 gut group was also identified as differentially abundant across horses in the genus-level analysis (in addition to another unclassified genus within Rikenellaceae). Rikenellaceae RC9 gut group has been identified across a range of herbivorous animal species [72,73,74] including the horse [75, 76] and has been found to increase in horses transitioned from pasture to a hay-based diet [59]. Members of the Rikenellaceae family are specialized to the animal gut [77], and while the role and function of the RC9 gut group within the gastrointestinal system has not been fully elucidated [76], this genus has been associated with degradation of structural carbohydrates. Asma et al. [78] found that abundance of Rikenellaceae RC9 group decreased as cattle maintained on silage-based diets adapted to supplemental concentrates higher in starch. Members of the Rikenellaceae family also degrade mucin, with a documented preference for mucin over simple sugars [78, 79]. It is possible that inter-horse variation in these fiber-degrading genera is indicative of functional redundancy within the hindgut microbial ecosystem and may serve as a mechanism for resilience, allowing horses to withstand subtle changes within the pasture environment.

It should be noted that a limitation of the current study is the relatively small sample size (8 horses). The smaller sample size in this study is similar to sample sizes utilized in prior equine hindgut microbiome studies that have emphasized a strong influence of individual horse [34, 48, 54]. However, it is possible that these small sample sizes are contributing to the strength of the influence of individual horse in comparison to other study variables, and studies in a greater number of horses would provide additional insight regarding any potential impact of transitions between forage types within an integrated grazing system. Additionally, Gomez et al. [35] described individualized responses to dietary intervention, the analysis of which was made possible by twice-daily serial sampling over a 5-d fecal collection period, yielding 10 samples per horse per treatment. In the current study, only one sample was collected on each of 4 days within each transition. While variations in BCG and ungrouped ASV relative abundance were evident in graphical representations of individual horses across days within transitions, additional replications would be required to statistically analyze shifts in bacterial abundances within individual horses. Further research is thus necessary to evaluate whether the microbiota of individual horses are indeed functionally stable across transitions between warm- and cool-season grasses, or to more definitively identify specific bacteria that are shifting in an individualized manner.

Results presented for the current study may have also been impacted by seasonality, as a seasonal control, such as horses maintained on an all-hay diet in which nutrient composition was constant, was not utilized. Seasonal shifts in the gut microbiota have been reported for other non-foraging species including mice [80], monkeys [81], and humans [82], but the effect of seasonality on the equine hindgut microbiome has not been extensively investigated. Theelen et al., [49] reported differences in β-diversity in horses sampled in the summer and winter, regardless of other animal and management variables. In pastured horses assessed monthly over a 12-mon period, Salem et al. [56] found fluctuations in microbial diversity and composition, which were associated with weather variables. This previous study did not evaluate the impact of nutritional composition of the pasture forage. The nutritional profile of pasture forages varies across the grazing season [1, 30, 83] and even over the course of a given day [84,85,86]. Given the established influence of diet and dietary nutrients on the gut microbiota [26,27,28] it is likely that changes in pasture forage nutrients also influenced microbial composition of grazing horses evaluated by Salem et al. [56]. Because horses in the current study were rotated through the grazing systems based on forage availability according to rotational grazing best management practices, transitions in each field could not be conducted on the same day. Furthermore, the establishment of BER was delayed in comparison to that of CRB, and thus the timing of productive grazing for these warm-season grasses did not align. However, analysis of individual transitions within each field also did not reveal substantive changes by day.

Additionally, many environmental and management factors can potentially influence nutrient content of pasture forages. Evaluation of the microbiome of horses grazing other cool-season and warm-season grass species and or varieties would be necessary to determine if the stability of the hindgut microbiome found in horses the current study extend to horses grazing in integrated rotational systems containing other forages. Collecting data from multiple years in multiple regional sites would also provide greater insight into the impacts of integrated warm- and cool-season grass rotational grazing systems on the equine hindgut microbial community.

Finally, a complete analysis of microbial functionality also was not conducted in the present study. While limited differences in the structure and composition of the microbiota were observed across days within transitions, it is possible that function could have differed even as bacterial composition remained relatively stable [87]. A number of ASV within the limited BCG varying over days within transitions between cool- and warm-season grasses were assigned to taxa known to degrade fiber as well as other carbohydrates and produce butyrate in addition to other short-chain fatty acids [88,89,90]. Groups of fibrolytic and butyrate-producing bacteria were differentially abundant across C–W transitions while separate and distinct groups containing ASV assigned to similar taxa and/or taxa with similar ascribed functions varied over the W–C transitions, suggesting functional redundancy and further reinforcing the overall stability of microbial communities during changes in pasture forages. While no differences were found in fecal pH, which is often cited as an indicator of bacterial activity in the hindgut [7, 91, 92], analyses of metabolites such as lactate and short-chain fatty acids or culture-based assays of functional communities (i.e. cellulytics, amylolytics, lactate-utilizers, etc.) would provide a more comprehensive evaluation and deeper understanding of hindgut microbiota function during transitions between warm- and cool-season pasture grasses.

Conclusions

In conclusion, equine fecal microbial community structure and composition as well as fecal pH are largely stable across transitions between warm-season grass and cool-season grass pasture sections within integrated warm- and cool-season rotational grazing systems. The capacity of the equine hindgut microbiota to adapt to these forages suggests that it is possible to manage healthy adult horses integrated rotational grazing systems of warm- and cool-season grasses without inducing dysfunction, although additional research is necessary to determine if these findings extend to other species and/or varieties of grasses under various environmental and management conditions. The individual horse was the strongest factor influencing the structure and composition of the gut microbiota during transitions, with bacteria representing over 30% of the microbial community differing between horses pointing to possible mechanisms of resilience in response to modest nutrient changes. Fiber-degraders were heavily represented in differentially abundant bacterial groups, regardless of factor (transition, grazing system, or individual horse), but further research is needed to determine if these findings indicate differences in fibrolytic capacity or are reflective of functional redundancy occurring in equine hindgut microbial communities.

Methods

Research was conducted in 2018 at the Ryders Lane Environmental Best Management Practices Demonstration Horse Farm (Rutgers, The State University of New Jersey; New Brunswick, New Jersey). Weather data for the New Brunswick station nearest to the site was obtained from the Historical Monthly Station Data portal of the Office of the New Jersey State Climatologist website [93]. Monthly average temperatures and precipitation totals across the study period and historical averages can be found in Additional File 6.

Grazing systems

Two separate 1.5 ha integrated warm- and cool-season rotational grazing systems were utilized: CRS and BRS, as described above. In CRS, 3 sections contained an established cool-season grass mix and the remaining 3 sections were planted with CRB. In BRS, 3 sections contained the established cool-season grass mix and the remaining 3 sections were planted with BER. The cool-season grass mix included Inavale Orchardgrass [Dactylis glomerata (L.)], Tower Tall Fescue (endophyte-free) [Lolium arundinaceum (Schreb.) Darbysh.], and Argyle Kentucky Bluegrass [Poa pratensis (L.)] (DLF Pickseed, Halsey, OR).

Animal and grazing management

Use of animals in this study was approved by the Rutgers University Institutional Animal Care and Use Committee protocol #PROTO201800013. Eight adult Standardbred mares with a body condition score (BCS; [94]) of 5–7 out of 9 were used for the study. Prior to grazing, horses were weighed and then grouped by age and BW. Groups of four horses were randomly assigned to each of the grazing systems. Initial mean age, BCS, and body weight (BW) for each group is as follows: BRS – 18.0 ± 0.85 yr; 5.6 ± 0.31 (BCS), 549 ± 24 kg; CRS – 17.7 ± 0.44 yr; 5.6 ± 0.13 (BCS), 533 ± 28 kg. Horses received regular veterinary and dental health care administered through Rutgers University Animal Care Program, with the most recent dental examinations conducted within six months prior to the start of the study.

In the spring of 2018 prior to first availability of pasture forage for grazing, each group of horses was housed in a separate dry lot and fed mixed cool-season grass hay at 2.5% BW per day. Horses began grazing on June 7, 2018, when cool-season grass pasture sections reached adequate height for grazing (at least 15.2 cm). Sequential grazing of warm-season grass sections began once the planted forage reached the same 15.2 cm minimum height. In CRS, the CRB was available for first grazing on August 5th. Despite being planted on the same date as CRB, the BER in BRS was slower to establish, with first grazing on September 18th. Weather data across Pastures were managed according to established best management practices for equine rotational grazing [1]. Forage availability dictated rotation frequency, with horses allowed to graze a given pasture section until forage was reduced to approximately 7.6 cm sward height, at which time horses were moved to a new section. After horses were removed from a grazed section, any remaining tall weeds were mowed to a height of 7.6 cm. Pasture sections were then dragged to evenly spread-out manure from defecation areas. Horses grazed on pastures until cool-season grasses became dormant at the end of the growing season in November. When on pasture, horses were allowed 24-h ad libitum access to pasture forage. If adequate pasture forage was not available at any point during the grazing season, horses were confined to a stress lot and supplemental grass hay was provided at 2.5% of body weight (BW) on a dry-matter (DM) basis. Horses had unlimited access to shelters, automatic waterers, and salt blocks throughout the grazing season.

Forage sampling

Representative hand-clipped forage samples were collected on the day of rotation. Pasture samples were collected between 0800 and 1000 by walking in a random zig-zag pattern throughout each pasture section, stopping to clip forage (at 7.6 cm height) every 30 paces [1]. Pasture samples were composited and then dried at 60 °C for at least 36 h in a Thelco oven (Precision Scientific, Chicago, IL) and ground to 1 mm using a Wiley Mill. Hay and pasture samples were then submitted to Equi-Analytical Laboratories (Ithaca, NY) to be analyzed by near-infrared spectroscopy.

Fecal and sample collection

Manual grab fecal samples were collected rectally from horses during C–W and W–C transitions on the D0, D2, D4, and D6 (see Additional file 7 for a diagram of experimental design and sampling protocol). Prior to the C–W rotation, horses had been grazing within the cool-season pasture sections for a minimum of 21 uninterrupted days (no confinement to stress lots or hay feeding). Similarly, horses grazed the CRB pasture sections in CRS for a minimum of 21 d prior to the W–C transitions. In BRS, only 17 days of grazing BER were possible prior to the W–C transition due to delayed establishment of BER. The duration of grazing prior to transition was selected to represent the longest duration that horses were likely to continuously consume an individual forage type within the context of integrated warm- and cool-season grass rotational grazing management. An adaptation period of 2–3 weeks has been utilized in prior studies evaluating relationships between the diet and the equine microbiome, including in grazing horses, with this period sufficient for stabilization of microbial communities [25, 37, 38]. All samples were collected at 0800 h. Samples were immediately placed on ice and transported from the field to the laboratory (5 min drive). Samples were then stored in a − 80 °C freezer until subsequent analysis.

Fecal sample analyses

Fecal pH was measured with a handheld Accumet pH meter (Fisher Scientific; Waltham, MA) calibrated with standard solutions at a pH of 4 and 7. Each sample was analyzed in duplicate, with fecal slurries of 10 mL of dH2O and 10 g of fecal matter thoroughly mixed by vortex in two separate 50 mL conical tubes. Two measurements were then collected from each of the duplicate tubes. The pH probe was inserted into the slurry and a reading was collected. The probe was then rinsed with dH2O to clean away any fecal material and reinserted to take a second pH reading. Readings were averaged for each horse at each sample point.

DNA was extracted from fecal samples in triplicate using Quick-DNA Fecal/Soil Microbe Kits (Zymo Research; Irvine, CA). Extracted DNA was quantified using a Qubit 2.0 Flourometer (Invitrogen; Carlsbad, CA). For each sample, the highest yielding replicate was selected and submitted to a commercial laboratory for amplification and sequencing of the V4-V5 region of the 16S rRNA gene (RTL Genomics; Lubbock, TX). Region specific primers (515F/926R) were used [95], and sequencing was conducted by Illumina MiSeq. Sequencing blanks were also prepared beginning with the extraction protocol, with sequencing producing zero reads per run for each blank.

Sequence and statistical analysis

Sequence and statistical analyses were performed in QIIME 2 (Quantitative Insights Into Microbial Ecology, v. 2020.8) [32] and R (v. 4.0.2) [96]. For all analyses, animal was used as the experimental unit. A record of code utilized in these analyses is included in Additional File 8.

Quality and chimera filtering of paired-end reads was conducted using DADA2 in Qiime 2, with read length set at 260 bp for forward reads and 200 bp for reverse reads [97, 98]. Trees for phylogenetic diversity analyses were generated using the mafft and FastTree pipeline from the q2-phylogeny Qiime 2 plugin [99,100,101]. Feature tables were filtered so that the lower quartile of ASV based on absolute abundance were removed, with filtering criteria set at a minimum frequency of 9 and presence in a minimum of 4 samples. Bacterial α-diversity and β-diversity analyses were also conducted in Qiime 2, with the feature table rarefied to an even minimum sampling depth of 3700. Differences in α-diversity metrics (Shannon Diversity Index, Faith’s Phylogenetic Diversity, Pielou’s Evenness, and Observed ASVs) were analyzed by mixed model ANOVA in R. System, transition, day, and their interactions were set as fixed factors and transition within horse as the random effect, with means separated by Tukey’s method [102,103,104,105,106]. The β-diversity metrics (Bray–Curtis, Jaccard, Unweighted UniFrac and Weighted Unifrac distance matrices) were analyzed by permutational ANOVA (PERMANOVA) with the Adonis action for multivariate analysis comparisons in Qiime 2 [107,108,109,110,111,112]. Influence of within-group variance on β-diversity analysis was evaluated by testing homogeneity of dispersion using PERMDISP for individual factors [113].

To further explore differential abundances, ASV were then grouped into BCG based on abundance profiles using SCNIC, with the correlation method set as Spearman and rs ≥ 0.50 [114, 115] with the q2-SCNIC plugin in QIIME 2. Random forest classifiers with nested cross validation (q2-sample-classifier plugin) were applied to determine if metadata variables (grazing system, field, day, horse) could be predicted based on BCG composition [116, 117]. Relative abundances of BCG and ASV that remained ungrouped following application of SCNIC were evaluated by ANCOM [118]. For ANCOM, H0(ij): mean(log[xi/xj) = mean(log[yi/yj). Strength of the ANCOM statistical test is denoted by W values, which indicate the number of times H0(ij) is rejected for the ith species. For example, if W = 770, the feature was significantly different relative to 770 other features [118]. Taxonomy was assigned to individual ASV using the latest SILVA database (SSU 138) in Qiime 2 [117, 119,120,121,122].

Shifts in microbial composition occurring at the genus level were also evaluated by (1.) collapsing ASV by genus, (2.) utilizing SCNIC to group genera by co-abundance, (3.) determining if the variables transition, grazing system, day (for each transition within each grazing system), and horse could be predicted based on microbial composition (random forest classification with nested cross-validation) and (4.) applying ANCOM to identify genera groups that were differentially abundant.

Fecal pH was also analyzed by mixed model ANOVA in R. System, transition, day, and their interactions were set as fixed factors and transition within horse as the random effect, with means separated by Tukey’s method. For α-diversity metrics and pH analyzed by mixed model, normality of model residuals was assessed using the Shapiro–Wilk test, and means were separated using Tukey’s method. Data for variables analyzed by mixed model are presented as means ± SEM. For all analyses which generated P-values, results were considered significant at P ≤ 0.05, with trends considered at P ≤ 0.10.

Availability of data and materials

The datasets generated and analyzed during the current study are available in the NCBI Sequence Read Archive at https://www.ncbi.nlm.nih.gov/sra, Bioproject: PRJNA804247, Accession numbers: SAMN25718984- SAMN25719047.

Abbreviations

ANCOM:

Analysis of Composition of Microbes

ASV:

Amplicon sequence variants

BER:

Wrangler Bermudagrass

BCG:

Bacterial co-abundance groups

BCS:

Body condition score

BRS:

Bermudagrass and cool-season grass integrated rotational grazing system

BW:

Body weight

CRB:

Quick-N-Big Crabgrass

CRS:

Crabgrass and cool-season grass integrated rotational grazing system

C–W:

Cool-season grass to warm-season grass transition

D0:

Day zero of rotation

D2:

Day two of rotation

D4:

Day four of rotation

D6:

Day six of rotation

DM:

Dry-matter

NSC:

Non-structural carbohydrates

PERMANOVA:

Permutational analysis of variance

PERMDISP:

Permutational analysis of dispersion

SCNIC:

Sparce Cooccurence Network Investigation for Compositional Data

WSC:

Water-soluble carbohydrates

W–C:

Warm-season grass to cool-season grass transition

References

  1. Williams CA, Kenny LB, Weinert JR, Sullivan K, Meyer W, Robson MG. Effects of 27 mo of rotational vs. continuous grazing on horse and pasture condition. Transl Anim Sci. 2020. https://doi.org/10.1093/tas/txaa084.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Taiz L, Zeiger E. Photosynthesis: carbon reactions. In: Taiz L, Zeiger E, editors. Plant physiology. 3rd ed. Sunderland: Sinauer Associates, Inc; 2002. p. 111–43.

    Google Scholar 

  3. Moore KJ, White TA, Hintz RL, Patrick PK, Brummer EC. Sequential grazing of cool-and warm-season pastures. Agron J. 2004;96(4):1103–11. https://doi.org/10.2134/agronj1991.00021962008300060027x.

    Article  Google Scholar 

  4. Tracy BF, Maughan M, Post N, Faulkner DB. Integrating annual and perennial warm-season grasses in a temperate grazing system. Crop Sci. 2010;50(5):2171–7. https://doi.org/10.2135/cropsci2010.02.0110.

    Article  Google Scholar 

  5. DeBoer ML, Sheaffer CC, Grev AM, Catalano DN, Wells MS, Hathaway MR, Martinson KL. Yield, nutritive value, and preference of annual warm-season grasses grazed by horses. Agron J. 2017;109(5):2136–48. https://doi.org/10.2134/agronj2017.02.0099.

    Article  CAS  Google Scholar 

  6. Ritz KE, Heins BJ, Moon R, Sheaffer C, Weyers SL. Forage yield and nutritive value of cool-season and warm-season forages for grazing organic dairy cattle. Agronomy. 2020;10(12):1963. https://doi.org/10.3390/agronomy10121963.

    Article  CAS  Google Scholar 

  7. Weinert-Nelson JR, Meyer WA, Williams CA. Yield, nutritive value, and horse condition in integrated crabgrass and cool-season grass rotational grazing pasture systems. Transl Anim Sci. 2021. https://doi.org/10.1093/tas/txab208 (in press).

    Article  PubMed  PubMed Central  Google Scholar 

  8. Teutsch C. Warm-season annual grasses for summer forage. Publication 418-004. Communication and marketing, College of Agriculture and Life Sciences, Virginia Polytechnic Inst. and State Univ.: Blacksburg; 2006.

  9. Teff KS (Eragrostis teff (Zucc.)). Trotter. Promoting the Conservation and use of the under utilized crops. In: vol. 12. Institute of Plant Genetics and Crop Plant Research, Garersleben/International Plant Genetic Resource Institute. Rome, Italy; 1997.

  10. Taliaferro CM. Breeding forage bermudagrass for the US Transition zone. In: Proceedings 59th southern pasture and forage crop improvement conference, Philadelphia, MS; 2005. p. 11–13.

  11. Ditsch DC, Smith SR, Lacefield GD. Bermudagrass: a summer forage in Kentucky. Publication #AGR-48. University of Kentucky College of Agriculture, Lexington, KY; 2011.

  12. Goodson J, Tyznik WJ, Cline JH, Dehority BA. Effects of an abrupt diet change from hay to concentrate on microbial numbers and physical environment in the cecum of the pony. Appl Environ Microbiol. 1988;54:1946–50. https://doi.org/10.1128/aem.54.8.1946-1950.1988.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Hudson JM, Cohen ND, Gibbs PG, Thompson JA. Feeding practices associated with colic in horses. J Am Vet Med Assoc. 2001;219(10):1419–25. https://doi.org/10.2460/javma.2001.219.1419.

    Article  CAS  PubMed  Google Scholar 

  14. Garner HE, Moore JN, Johnson JH, ClarkL AJF, Tritschler LG, Coffmann JR, Sprouse RF, Hutcheson DP, Salem CA. Changes in the caecal flora associated with the onset of laminitis. Equine Vet J. 1978;10:249–52. https://doi.org/10.1111/j.2042-3306.1978.tb02273.x.

    Article  CAS  PubMed  Google Scholar 

  15. Millinovich GJ, Burrell PC, Pollitt CC, Klieve AV, Blackall LL, Ouwerkerk D, Woodland E, Trott DJ. Microbial ecology of the equine hindgut during oliofructose-induced laminitis. ISME J. 2008;2:1089–100. https://doi.org/10.1038/ismej.2008.67.

    Article  CAS  Google Scholar 

  16. Tuniyazi M, He J, Guo J, Li S, Zhang N, Hu X, Fu Y. Changes of microbial and metabolome of the equine hindgut during oligofructose-induced laminitis. BMC Vet Res. 2021;17(1):1–13. https://doi.org/10.1186/s12917-020-02686-9.

    Article  CAS  Google Scholar 

  17. Cohen ND, Matejka PL, Honnas CM, Hooper RN. Case-control study of the association between various management factors and development of colic in horses. Texas equine colic study group. J Am Vet Med Assoc. 1995;206(5):667–73.

    CAS  PubMed  Google Scholar 

  18. Tinker MK, White NA, Lessard P, Thatcher CD, Pelzer KD, Davis B, Carmel DK. Prospective study of equine colic risk factors. Equine Vet J. 1997;29(6):454–8. https://doi.org/10.1111/j.2042-3306.1997.tb03158.x.

    Article  CAS  PubMed  Google Scholar 

  19. Venable E, Kerley MS, Raub R. Assessment of equine fecal microbial profiles during and after a colic episode using pyrosequencing. J Equine Vet Sci. 2013;33:347. https://doi.org/10.1016/j.jevs.2013.03.066.

    Article  Google Scholar 

  20. Weese JS, Holcombe SJ, Embertson RM, Kurtz KA, Roessner HA, Jalali M, Wismer SE. Changes in the faecal microbiota ofmares precede the development of post partum colic. Equine Vet J. 2015;47:641–9. https://doi.org/10.1111/evj.12361.

    Article  CAS  PubMed  Google Scholar 

  21. Stewart HL, Southwood LL, Indugu N, Vecchiarelli B, Engiles JB, Pitta D. Differences in the equine faecal microbiota between horses presenting to a tertiary referral hospital for colic compared with an elective surgical procedure. Equine Vet J. 2019;51(3):336–42. https://doi.org/10.1111/evj.13010.

    Article  CAS  PubMed  Google Scholar 

  22. United States Department of Agriculture. Lameness and laminitis in US horses. USDA: APHIS: US, CEAH, National Animal Health Monitoring System. United States Department of Agriculture, Washington DC; 2000.

  23. United States Department of Agriculture. Baseline reference of equine health and management in the United States, 2015. USDA: APHIS: US, CEAH, National Animal Health Monitoring System. US Department of Agriculture, Washington DC. 2016.

  24. Troya L, Blanco J, Romero I, Re M. Comparison of the colic incidence in a horse population with or without inclusion of germinated barley in the diet. Equine Vet Educ. 2020;32:28–32. https://doi.org/10.1111/eve.13274.

    Article  Google Scholar 

  25. Fernandes KA, Kittelmann S, Rogers CW, Gee EK, Bolwell CF, Thomas BEN, DG,. Faecal microbiota of forage-fed horses in New Zealand and the population dynamics of microbial communities following dietary change. PLoS ONE. 2014;9(11):e112846. https://doi.org/10.1371/journal.pone.0112846.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Zhang C, Zhang M, Wang S, Han R, Cao Y, Hua W, Mao Y, Zhang X, Pang X, Wei C. Interactions between gut microbiota, host genetics and diet relevant to development of metabolic syndromes in mice. ISME J. 2010;4(2):232. https://doi.org/10.1038/ismej.2009.112.

    Article  CAS  PubMed  Google Scholar 

  27. Zhang C, Li S, Yang L, Huang P, Li W, Wang S, Zhao G, Zhang M, Pang X, Yan Z. Structural modulation of gut microbiota in life-long calorie-restricted mice. Nat Commun. 2013;4:2163. https://doi.org/10.1038/ncomms3163.

    Article  CAS  PubMed  Google Scholar 

  28. Dougal K, de la Fuente G, Harris PA, Girdwood SE, Pinloche E, Geor RJ, Nielsen BD, Schott HC II, Elzinga S, Newbold CJ. Characterisation of the faecal bacterial community in adult and elderly horses fed a high fibre, high oil or high starch diet using 454 pyrosequencing. PLoS ONE. 2014;9(2):e87424. https://doi.org/10.1371/journal.pone.0087424.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Chatterton NJ, Harrison PA, Bennett JH, Asay KH. Carbohydrate partitioning in 185 accessions of gramineae grown under warm and cool temperatures. J Plant Physiol. 1989;134(2):169–79. https://doi.org/10.1016/S0176-1617(89)80051-3.

    Article  CAS  Google Scholar 

  30. Jensen KB, Harrison P, Chatterton NJ, Bushman BS, Creech JE. Seasonal trends in nonstructural carbohydrates in cool-and warm-season grasses. Crop Sci. 2014;54(5):2328–40. https://doi.org/10.2135/cropsci2013.07.0465.

    Article  CAS  Google Scholar 

  31. Hudson DJ, Leep RH, Dietz TS, Ragavendran A, Kravchenko A. Integrated warm-and cool-season grass and legume pastures: I. seasonal forage dynamics. Agron J. 2010;102(1):303–9. https://doi.org/10.2134/agronj2009.0204.

    Article  Google Scholar 

  32. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, Alexander H, Alm EJ, Arumugam M, Asnicar F, Bai Y. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7. https://doi.org/10.1038/s41587-019-0209-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Pelletier S, Tremblay GF, Bertrand A, Belanger G, Castonguay Y, Michaud R. Drying procedures affect non-structural carbohydrates and other nutritive value attributes in forage samples. Anim Feed Sci Technol. 2010;157:139–50. https://doi.org/10.1016/j.anifeedsci.2010.02.010.

    Article  CAS  Google Scholar 

  34. Garber A, Hastie P, McGuinness D, Malarange P, Murray JA. Abrupt dietary changes between grass and hay alter faecal microbiota of ponies. PLoS ONE. 2020;15(8):e0237869. https://doi.org/10.1371/journal.pone.0237869.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Muhonen S, Connysson M, Lindberg JE, Julliand V, Bertilsson J, Jansson A. Effects of crude protein intake from grass silage-only diets on the equine colon ecosystem after an abrupt feed change. J Anim Sci. 2008;86(12):3465–72. https://doi.org/10.2527/jas.2007-0374.

    Article  CAS  PubMed  Google Scholar 

  36. Grimm P, Philippeau C, Julliand V. Faecal parameters as biomarkers of the equine hindgut microbial ecosystem under dietary change. Animal. 2017;11(7):1136–45. https://doi.org/10.1017/S1751731116002779.

    Article  CAS  PubMed  Google Scholar 

  37. Fitzgerald DM, Spence RJ, Stewart ZK, Prentis PJ, Sillence MN, De Laat MA. The effect of diet change and insulin dysregulation on the faecal microbiome of ponies. J Exper Biol. 2020;223(7):jeb219154. https://doi.org/10.1242/jeb.219154.

    Article  Google Scholar 

  38. Respondek F, Goachet A, Julliand RFV. Effects of short-chain fructo-oligosaccharides on the microbial and biochemical profile of different segments of the gastro-intestinal tract in horses. Pferdeheilkunde. 2008;23(2):146. https://doi.org/10.21836/PEM20070206.

    Article  Google Scholar 

  39. De Fombelle A, Julliand V, Drogoul C, Jacotot E. Feeding and microbial disorders in horses: 1-effects of an abrupt incorporation of two levels of barley in a hay diet on microbial profile and activities. J Equine Vet Sci. 2001;21:439–45. https://doi.org/10.1016/S0737-0806(01)70018-4.

    Article  Google Scholar 

  40. Warzecha CM, Coverdale JA, Janecka JE, Leatherwood JL, Pinchak WE, Wickersham TA, McCann JC. Influence of short-term dietary starch inclusion on the equine cecal microbiome. J Anim Sci. 2017;95(11):5077–90. https://doi.org/10.2527/jas2017.1754.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Muhonen S, Julliand V, Lindberg JE, Bertilsson J, Jansson A. Effects on the equine colon ecosystem of grass silage and haylage diets after an abrupt change from hay. J Anim Sci. 2009;87(7):2291–8. https://doi.org/10.2527/jas.2008-1461.

    Article  CAS  PubMed  Google Scholar 

  42. Zhang C, Zhao L. Strain-level dissection of the contribution of the gut microbiome to human metabolic disease. Genome Med. 2016;8(1):1–10. https://doi.org/10.1186/s13073-016-0304-1.

    Article  Google Scholar 

  43. Pan F, Zhang L, Li M, Hu Y, Zeng B, Yuan H, Zhao L, Zhang C. Predominant gut Lactobacillus murinus strain mediates anti-inflammaging effects in calorie-restricted mice. Microbiome. 2018;6(1):1–17. https://doi.org/10.1186/s40168-018-0440-5.

    Article  CAS  Google Scholar 

  44. Zhai R, Xue X, Zhang L, Yang X, Zhao L, Zhang C. Strain-specific anti-inflammatory properties of two Akkermansia muciniphila strains on chronic colitis in mice. Front Cell Infect Microbiol. 2019;9:239. https://doi.org/10.3389/fcimb.2019.00239.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Wu G, Zhao N, Zhang C, Lam YY, Zhao L. Guild-based analysis for understanding gut microbiome in human health and diseases. Genome Med. 2021;13(1):1–12. https://doi.org/10.1186/s13073-021-00840-y.

    Article  Google Scholar 

  46. Mach N, Ruet A, Clark A, Bars-Cortina D, Ramayo-Caldas Y, Crisci E, Pennarun S, Dhorne-Pollet S, Foury A, Moisan MP, Lansade L. Priming for welfare: gut microbiota is associated with equitation conditions and behavior in horse athletes. Sci Rep. 2020;10(1):1–19. https://doi.org/10.1038/s41598-020-65444-9.

    Article  CAS  Google Scholar 

  47. Husso A, Jalanka J, Alipour MJ, Huhti P, Kareskoski M, Pessa-Morikawa T, Iivanainen A, Niku M. The composition of the perinatal intestinal microbiota in horse. Sci Rep. 2020;10(1):1–12. https://doi.org/10.1038/s41598-019-57003-8.

    Article  CAS  Google Scholar 

  48. Gomez A, Sharma AK, Grev A, Sheaffer C, Martinson K. The horse gut microbiome responds in a highly individualized manner to forage lignification. J Equine Vet Sci. 2021;96: 103306. https://doi.org/10.1016/j.jevs.2020.103306.

    Article  PubMed  Google Scholar 

  49. Theelen MJ, Luiken RE, Wagenaar JA, Sloet van Oldruitenborgh-Oosterbaan MM, Rossen JW, Zomer AL. The equine faecal microbiota of healthy horses and ponies in The Netherlands: impact of host and environmental factors. Animals. 2021;11(6):1762. https://doi.org/10.3390/ani11061762.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Zhang C, Yin A, Li H, Wang R, Wu G, Shen J, Zhang M, Wang L, Hou Y, Ouyang H, Zhang Y. Dietary modulation of gut microbiota contributes to alleviation of both genetic and simple obesity in children. EBioMedicine. 2015;2(8):968–84. https://doi.org/10.1016/j.ebiom.2015.07.007.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Zhao L, Zhang F, Ding X, Wu G, Lam YY, Wang X, Fu H, Xue X, Lu C, Ma J, Yu L. Gut bacteria selectively promoted by dietary fibers alleviate type 2 diabetes. Science. 2018;359(6380):1151–6. https://doi.org/10.1126/science.aao5774.

    Article  CAS  PubMed  Google Scholar 

  52. Chen T, Liu AB, Sun S, Ajami NJ, Ross MC, Wang H, Zhang L, Reuhl K, Kobayashi K, Onishi JC, Zhao L, Yang CS. Green tea polyphenols modify the gut microbiome in db/db mice as co-abundance grouips correlating with the blood glucose lowering effect. Mol Nutr Food Res. 2019;63(8): 180164. https://doi.org/10.1002/mnfr.201801064.

    Article  CAS  Google Scholar 

  53. Blackmore TM, Dugdale A, Argo CM, Curtis G, Pinloche E, Harris PA, Worgan HJ, Girwood SE, Dougal K, Newbold CJ, McEwan NR. Strong stability and host specific bacterial community in faeces of ponies. PLoS ONE. 2013;8(9): e75079. https://doi.org/10.1371/journal.pone.0075079.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Costa MC, Weese JS. The equine intestinal microbiome. Anim Health Res Rev. 2012;13(1):121–8. https://doi.org/10.1017/S1466252312000035.

    Article  PubMed  Google Scholar 

  55. Proudman A, Darby C, Escalona E. Faecal microbiome of the Thoroughbred racehorse and its response to dietary amylase supplementation. Equine Vet J. 2014;46(S46):35. https://doi.org/10.1111/evj.12267_107.

    Article  Google Scholar 

  56. Salem SE, Maddox TW, Berg A, Antczak P, Ketley JM, Williams NJ, Archer DC. Variation in faecal microbiota in a group of horses managed at pasture over a 12-month period. Sci Rep. 2018;8(1):8510. https://doi.org/10.1038/s41598-018-26930-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Johnson AJ, Vangay P, Al-Ghalith GA, Hillman BM, Ward TL, Shields-Cutler RR, Kim AD, Shmagel AK, Syed AN, Personalized Microbiome Students, Walter J. Daily sampling reveals personalized diet-microbiome associations in humans. Cell Host Microbe. 2019;25(6):789–802. https://doi.org/10.1016/j.chom.2019.05.005.

    Article  CAS  PubMed  Google Scholar 

  58. Smits SA, Marcobal A, Higginbottom S, Sonnenburg JL, Kashyap PC. Individualized responses of gut microbiota to dietary intervention modeled in humanized mice. mSystems. 2016;1(5):e00098. https://doi.org/10.1128/mSystems.00098-16.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Ericsson AC, Johnson PJ, Gieche LM, Zobrist C, Bucy K, Townsend KS, Martin LM, LaCarrubba AM. The influence of diet change and oral metformin on blood glucose regulation and the fecal microbiota of healthy horses. Animals. 2021;11(4):976. https://doi.org/10.3390/ani11040976.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Zhu Y, Wang X, Deng L, Chen S, Zhu C, Li J. Effects of pasture grass, silage, and hay diet on equine fecal microbiota. Animals. 2021;11(5):1330. https://doi.org/10.3390/ani11051330.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Goodrich JK, Waters JL, Poole AC, Sutter JL, Koren O, Blekhman R, Beaumont M, Van Treuren W, Knight R, Bell JT, Spector TD, Clark AG, Ley RE. Human genetics shape the gut microbiome. Cell. 2014;159(4):789–99.

    Article  CAS  Google Scholar 

  62. Svartström O, Alneberg J, Terrapon N, Lombard V, de Bruijn I, Malmsten J, Dalin A, Muller EEL, Shah P, Wilmes P, Henrissat B, Aspeborg H, Andersson AF. Ninety-nine de novo assembled genomes from the moose (Alces alces) rumen microbiome provide new insights into microbial plant biomass degradation. ISME J. 2017;11:2538–51. https://doi.org/10.1038/ismej.2017.108.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. La Reau AJ, Suen G. The Ruminocci: key symbionts of the gut ecosystem. J Microbiol. 2018;56(3):199–208. https://doi.org/10.1007/s12275-018-8024-4.

    Article  CAS  PubMed  Google Scholar 

  64. Tokuda G, Mikaelyan A, Fukui C, Watanabe H, Funishima M, Brune A. Fiber-associated spirochetes are major agents of hemicellulose degradation in the hindgut of wood-feeding higher termites. PNAS. 2018;115(51):E11996–2004. https://doi.org/10.1073/pnas.1810550115.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Ren Q, Si H, Yan X, Liu C, Ding L, Long R, Li Z, Qiu Q. Bacterial communities in the solid, liquid, dorsal, and ventral epithelium fractions of yak (Bos grunniens) rumen. Microbiologyopen. 2020;9(2):e963. https://doi.org/10.1002/mbo3.963.

    Article  CAS  PubMed  Google Scholar 

  66. Vital M, Jairong G, Rizzo R, Harrison T, Tiedje JM. Diet is a major factor governing the fecal butyrate-producing community structure across Mammalia, Aves and Reptilia. ISME J. 2015;9:832–43. https://doi.org/10.1038/ismej.2014.179.

    Article  CAS  PubMed  Google Scholar 

  67. Perea K, Perz K, Olivo SK, Williams A, Lachman M, Ishaq SL, Thomson J, Yeoman CJ. J Anim Sci. 2017;95(6):2585–92. https://doi.org/10.2527/jas.2016.1222.

    Article  CAS  PubMed  Google Scholar 

  68. Gharechahi J, Vahidi MF, Ding X-Z, Han J-L, Salekdeh GH. Temporal changes in microbial communities attached to forages with different lignocellulosic compositions in cattle rumen. FEMS Microbiol Ecol. 2020. https://doi.org/10.1093/femsec/fiaa069.

    Article  PubMed  Google Scholar 

  69. Goodrich JK, Davenport ER, Waters JL, Clark AG, Ley RE. Cross-species comparisons of host genetic associations with the microbiome. Science. 2016;352:532–5. https://doi.org/10.1126/science.aad9379.DOI:10.1016/j.cell.2014.09.053.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Lim MY, You HJ, Yoon HS, Kwon B, Lee JY, Lee S, Song Y, Lee K, Sung J, Ko G. The effect of heritability and host genetics on the gut microbiota and metabolic syndrome. Gut. 2017;66:1031–8. https://doi.org/10.1136/gutjnl-2015-311326.

    Article  CAS  PubMed  Google Scholar 

  71. Waters JL, Ley RE. The human gut bacteria Christensenellaceae are widespread, heritable, and associated with health. BMC Biol. 2019;17:83. https://doi.org/10.1186/s12915-019-0699-4.

    Article  PubMed  PubMed Central  Google Scholar 

  72. Ilmberger N, Güllert S, Dannenberg J, Rabausch U, Torres J, Wemheuer B, Alawi M, Poehlein A, Chow J, Turaev D, Rattei T. A comparative metagenome survey of the fecal microbiota of a breast- an a plant-fed Asian elephant reveals an unexpectedly high diversity of glycoside hydrolase family enzymes. PLoS ONE. 2014;9(9): e106707. https://doi.org/10.1371/journal.pone.0106707.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Li Y, Hu X, Yang S, Zhou J, Zhang T, Qi L, Sun X, Fan M, Xu S, Cha M, Zhang M. Comparative analysis of the gut microbiota composition between captive and wild forest musk deer. Front Microbiol. 2017;8:1705. https://doi.org/10.3389/fmicb.2017.01705.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Huang Q, Holman BD, Alexander T, Hu T, Jin L, Xu Z, McAllister TA, Acharya S, Zhao G, Wang Y. Fecal microbiota of lambs fed purple prairie clover (Dalea purpurea Vent) and alfalfa (Medicago sativa). Arch Microbiol. 2018;200(1):137–45. https://doi.org/10.1007/s00203-017-1427-5.

    Article  CAS  PubMed  Google Scholar 

  75. Rodriquez C, Taminiau B, Brévers B, Avesani V, Van Broeck J, Leroux A, Gallot M, Bruwier A, Amory H, Delmée M, Daube G. Faecal microbiota characterisation of horses using 16 rdna barcoded pyrosequencing, and carriage rate of clostridium difficile at hospital admission. BMC Microbiol. 2015;15(1):1–14. https://doi.org/10.1186/s12866-015-0514-5.

    Article  CAS  Google Scholar 

  76. Li Y, Zhang K, Yang L, Kai L, Defu H, Wronski T. Community composition and diversity of intestinal microbiota in captive and re-introduced Prezwalski’s Horse (Equus ferus prezwalskii). Front Microbiol. 2019;10:1821. https://doi.org/10.3389/fmicb.2019.01821.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Graf J. The family Rikenellaceae. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, editors. The prokaryotes. Berlin: Springer Berlin Heidelberg; 2014. p. 857–9. https://doi.org/10.1007/978-3-642-38954-2_134.

    Chapter  Google Scholar 

  78. Asma Z, Sylvie C, Laurent C, Jérôme M, Christophe K, Oliver B, Annabelle TM, Francis E. Microbial ecology of the rumen evaluated by 454 GS FLX pyrosequencing is affected by starch and oil supplementation of diets. FEMS Microbio Ecol. 2013;83(2):504–14. https://doi.org/10.1111/1574-6941.12011.

    Article  CAS  Google Scholar 

  79. Bomar I, Malz M, Colston S, Graf J. Directed culturing of microorganisms using metatranscriptomics. MBio. 2011;2(2):e00012-e111. https://doi.org/10.1128/mBio.00012-11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Maurice CF, Knowles SC, Ladau J, Pollard KS, Fenton A, Pedersen AB, Turnbaugh PJ. Marked seasonal variation in the wild mouse gut microbiota. ISME J. 2015;9(11):2423–34. https://doi.org/10.1038/ismej.2015.53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Amato KR, Leigh SR, Kent A, Mackie RI, Yeoman CJ, Stumpf RM, Wilson BA, Nelson KE, White BA, Garber PA. The gut microbiota appears to compensate for seasonal diet variation in the wild black howler monkey (Alouatta pigra). Microb Ecol. 2015;69(2):434–43. https://doi.org/10.1007/s00248-014-0554-7.

    Article  CAS  PubMed  Google Scholar 

  82. Parfrey LW, Knight R. Spatial and temporal variability of the human microbiota. Clin Microbiol Infect. 2012;18(S4):5–7. https://doi.org/10.1111/j.1469-0691.2012.03861.x.

    Article  Google Scholar 

  83. Williams CA, Kenny LB, Burk AO. Effects of grazing system, season, and forage carbohydrates on glucose and insulin dynamics of the grazing horse. J Anim Sci. 2019;97(6):2541–54. https://doi.org/10.1093/jas/skz103.

    Article  PubMed  PubMed Central  Google Scholar 

  84. Kagan IA, Kirch BH, Thatcher CD, Strickland JR, Teutsch CD, Elvinger F, Pleasant RS. Seasonal and diurnal variation in simple sugar and fructan composition of orchardgrass pasture and hay in the Piedmont region of the United States. J Equine Vet Sci. 2011;31(8):488–97. https://doi.org/10.1016/j.jevs.2011.03.004.

    Article  Google Scholar 

  85. Kagan IA, Kirch BH, Thatcher CD, Teutsch CD, Elvinger F, Shepherd DM, Pleasant S. Seasonal and diurnal changes in starch content and sugar profiles of Bermudagrass in the Piedmont region of the United States. J Equine Veterinary Sci. 2011;31(9):521–9. https://doi.org/10.1016/j.jevs.2011.08.010.

    Article  Google Scholar 

  86. Weinert-Nelson JR, Meyer WA, Williams CA. Diurnal variation in forage nutrient composition of mixed cool-season grass, crabgrass, and bermudagrass pastures. J Equine Vet Sci. 2022;110: 103836. https://doi.org/10.1016/j.jevs.2021.103836.

    Article  PubMed  Google Scholar 

  87. Berg EL, Fu CJ, Porter JH, Kerley MS. Fructooligosaccharide supplementation in the yearling horse: effects on fecal pH, microbial content, and volatile fatty acid concentrations. J Anim Sci. 2005;83(7):1549–53. https://doi.org/10.2527/2005.8371549x.

    Article  CAS  PubMed  Google Scholar 

  88. Biddle AS, Stewart L, Blanchard J, Leschine S. Untangling the genetic basis of fibrolytic specialization by Lachnospiraceae and Ruminococcaceae in Diverse Gut Communities. Diversity. 2013;5(3):627–40. https://doi.org/10.3390/d5030627.

    Article  Google Scholar 

  89. Lawson PA, Rainey FA. Proposal to restrict the genus Clostridium Prazmowski to Clostridium butyricum and related species. Int J Syst Evol. 2016;66(2):1009–16. https://doi.org/10.1099/ijsem.0.000824.

    Article  CAS  Google Scholar 

  90. La Reau AJ, Suen G. The Ruminococci: key symbionts of the gut ecosystem. J Microbiol. 2018;56(3):199–208. https://doi.org/10.1007/s12275-018-8024-4.

    Article  CAS  PubMed  Google Scholar 

  91. Willing B, Vörös A, Roos S, Jones C, Jansson A, Lindberg J. Changes in faecal bacteria associated with concentrate and forage-only diets fed to horses in training. Equine Vet J. 2009;41:908–14. https://doi.org/10.2746/042516409X447806.

    Article  CAS  PubMed  Google Scholar 

  92. Sorensen RJ, Drouillard JS, Douthit TL, Ran Q, Marthaler DG, Kang Q, Vahl CI, Lattimer JM. Effect of hay type on cecal and fecal microbiome and fermentation parameters in horses. J Anim Sci. 2021. https://doi.org/10.1093/jas/skaa407.

    Article  PubMed  Google Scholar 

  93. Office of the New Jersey State climatologist at Rutgers University: Rutgers New Jersey weather network. https://www.njweather.org/data (2021). Accessed 12 Jul 2021.

  94. Henneke DR, Potter GD, Kreider JL, Yeates BF. Relationship between condition score, physical measurements and body fat percentage in mares. Equine Vet J. 1983;15(4):371–2. https://doi.org/10.1111/j.2042-3306.1983.tb01826.x.

    Article  CAS  PubMed  Google Scholar 

  95. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, Owens M, Betley J, Fraser L, Bauer M, Gormley N. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6:1621–4. https://doi.org/10.1038/ismej.2012.8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  96. R Development Core Team. R: A language and environment for statistical computing. 2010. http://cran.r-project.org.

  97. McDonald D, Clemente JC, Kuczynski J, Rideout JR, Stombaugh J, Wendel D, Wilke A, Huse S, Hufnagle J, Meyer F, Knight R. The Biological observation matrix (BIOM) format or: how I learned to stop worrying and love the ome-ome. Gigascience. 2012;1(1):2047–217. https://doi.org/10.1186/2047-217X-1-7.

    Article  Google Scholar 

  98. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581–3. https://doi.org/10.1038/nmeth.3869.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  99. Lane DJ. 16S/23S rRNA Sequencing. In: Stakebrandt E, Goodfellow M, editors. Nucleic acid techniques in bacterial systematics. New York City: John Wiley and Sons; 1991. p. 115–75.

    Google Scholar 

  100. Price MN, Dehal PS, Arkin AP. FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS ONE. 2010;5(3):e9490. https://doi.org/10.1371/journal.pone.0009490.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  101. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80. https://doi.org/10.1093/molbev/mst010.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  102. Shannon CE. A mathematical theory of communication. Bell Sys Tech J. 1948;27(3):379–423. https://doi.org/10.1002/j.1538-7305.1948.tb01338.x.

    Article  Google Scholar 

  103. Pielou EC. The measurement of diversity in different types of biological collections. J Theor Biol. 1966;13:131–44. https://doi.org/10.1016/0022-5193(66)90013-0.

    Article  Google Scholar 

  104. Faith D. Conservation evaluation and phylogenetic diversity. Biol Conserv. 1992;61(1):1–10. https://doi.org/10.1016/0006-3207(92)91201-3.

    Article  Google Scholar 

  105. McKinney W. Data structures for statistical computing in python. In: van der Walt S, Millman J, editors. Proceedings of the 9th python in science conference; 2010. p. 51–6.

  106. Weiss S, Xu ZZ, Peddada S, Amir A, Bittinger K, Gonzalez A, Lozupone C, Zaneveld JR, Vázquez-Baeza Y, Birmingham A, Hyde ER. Normalization and microbial differential abundance strategies depend upon data characteristics. Microbiome. 2017;5(1):1–18. https://doi.org/10.1186/s40168-017-0237-y.

    Article  Google Scholar 

  107. Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Envir Microbiol. 2005;71(12):8228–35. https://doi.org/10.1128/AEM.71.12.8228-8235.2005.

    Article  CAS  Google Scholar 

  108. Lozupone CA, Hamady M, Kelley ST, Knight R. Quantitative and uqalitative β diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol. 2007;73(5):1576–85. https://doi.org/10.1128/AEM.01996-06.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  109. Hamady M, Lozupone C, Knight R. Fast unifrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequening and PhyloChip data. ISME J. 2010;4(1):17–27. https://doi.org/10.1038/ismej.2009.97.

    Article  CAS  PubMed  Google Scholar 

  110. Chang Q, Luan Y, Sun F. Variance adjusted weighted UniFrac: a powerful beta diversity measure for comparing communities based on phylogeny. BMC Bioinform. 2011. https://doi.org/10.1186/1471-2105-12-118.

    Article  Google Scholar 

  111. Chen J, Bittinger K, Charlson ES, Hofmann C, Lewis J, Wu GD, Collman G, Bushman FD, Li H. Associating microbiome composition with environmental covariates using generalized UniFrac distances. Bioinformatics. 2012;28(16):2106–13. https://doi.org/10.1093/bioinformatics/bts342.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  112. McDonald D, Vázquez-Baeza Y, Koslicki D, McClelland J, Reeve N, Zhenjiang X, Gonzalez A, Knight R. Striped UniFrac: enabling microbiome analysis at unprecedented scale. Nat Methods. 2018;15(11):847–8. https://doi.org/10.1038/s41592-018-0187-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  113. Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26(1):32–46. https://doi.org/10.1111/j.1442-9993.2001.01070.pp.x.

    Article  Google Scholar 

  114. Hagberg AA, Shult DA, Swart PJ. Exploring network structure, dynamics, and function using NetworkX. In: Varoquaux G, Vaught T, Millman J, editors. Proceedings of the 7th Python in Science Conference; 2008. p. 11–15.

  115. Shaffer M, Thurimella K, Lozupone CA. SCNIC: Sparse correlation network investigation for compositional data. bioRxiv. 2020. https://doi.org/10.1101/2020.11.13.380733.

    Article  PubMed  PubMed Central  Google Scholar 

  116. Bokulich N, Dillon M, Bolyen E, Kaehler BD, Huttley GA, Caporaso JG. q2-sample-classifier: machine-learning tools for microbiome classification and regression. J Open Source Softw. 2018;3(30):934. https://doi.org/10.21105/joss.00934.

    Article  Google Scholar 

  117. Pedregosa F, Varoquaux G, Gramfort A, Michel B, Thirion B, Grisel O, Blondel M, Prettenhofer P, Weiss R, Dubourg V, Vanderplas J, Passos A, Cournapeau D, Brucher M, Perrot M, Duchesnay E. Scikit-learn: machine learning in Python. J Mach Learn Res. 2011;12:2825–30.

    Google Scholar 

  118. Mandal S, Van Treuren W, White RA, Eggesbø M, Knight R, Peddada SD. Microb Ecol Health Dis. 2015;26(1):27663. https://doi.org/10.3402/mehd.v26.27663.

    Article  PubMed  Google Scholar 

  119. Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig K, Peplies J, Glockner FO. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucl Acids Res. 2007;35:7188–96.

    Article  CAS  Google Scholar 

  120. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Pablo J, Glockner FO. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucl Acids Res. 2013;41:D590–6.

    Article  CAS  Google Scholar 

  121. Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4: e2584. https://doi.org/10.7717/peerj.2584.

    Article  PubMed  PubMed Central  Google Scholar 

  122. Bokulich NA, Kaehler BD, Rideout JR, Dillon M, Bolyen E, Knight R, Huttley GA, Caparaso JG. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-eature-classifier plugin. Microbiome. 2018;1(6):90. https://doi.org/10.1186/s40168-018-0470-z.

    Article  Google Scholar 

Download references

Acknowledgements

We thank DLF Pickseed (Halsey, OR) and Dalrymple Farms (Thomas, OK) for donations of grass seed used to establish the pastures grazed in this study. We also thank the Rutgers Equine Science Center for providing office space and computers utilized for data analyses.

Funding

This work was supported by the USDA National Institute of Food and Agriculture Hatch project 1003557 through the New Jersey Agricultural Experiment Station, Hatch project NJ06170 and the USDA Northeast Sustainable Agriculture, Research, and Education project gne17-162. Additional funding was provided by the Rutgers University Equine Science Center (57 US 1, New Brunswick, NJ 08901).

Author information

Authors and Affiliations

Authors

Contributions

JRW designed and conducted the study as well as bioinformatics analysis and wrote the majority of the manuscript. ASB assisted with study design, the workflow for bioinformatics analysis of data, and editing of the manuscript. CAW secured funding for the project in addition to assisting with experimental design, subject sampling, and manuscript editing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jennifer R. Weinert-Nelson.

Ethics declarations

Ethics approval and consent to participate

All animal procedures were conducted following Rutgers University Institutional Animal Care and Use Committee guidelines (protocol #PROTO201800013).

Consent for publication

Not applicable.

Competing interests

No conflict of interest declared by any of the authors.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Taxonomic classification of amplicon sequence variants (ASV) retained after quality filtering, chimera checking, and filtering low-abundance ASV.

Additional file 2:

Differentially abundant bacterial co-abundance groups (BCG) and ungrouped amplicon sequence variants (ASV) across horses.

Additional file 3:

Relative abundance of the twelve most abundant BCG that differed by horse.

Additional file 4:

Genera groups and/or ungrouped genera identified as differentially abundant by grazing system, transition, day (within each grazing system for each transition) and horse.

Additional file 5:

Analysis comparison of grouping amplicon sequence variants into bacterial co-abundance groups (BCG) versus genus-level groupings.

Additional file 6:

Monthly average temperatures and precipitation totals across the study period and historical averages.

Additional file 7:

Diagram of experimental design and sampling protocol.

Additional file 8:

Code used for sequence and statistical analysis in Qiime 2 (v.2020.8) and R (v. 4.0.2).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Weinert-Nelson, J.R., Biddle, A.S. & Williams, C.A. Fecal microbiome of horses transitioning between warm-season and cool-season grass pasture within integrated rotational grazing systems. anim microbiome 4, 41 (2022). https://doi.org/10.1186/s42523-022-00192-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s42523-022-00192-x

Keywords