# Weaning age influences the severity of gastrointestinal microbiome shifts in dairy calves

Sarah Jade Meale, S. C. Li, P. Azevedo, H. Derakhshani, T. J. Devries, J. C.

Plaizier, M. A. Steele, E. Khafipour

► To cite this version:

Sarah Jade Meale, S. C. Li, P. Azevedo, H. Derakhshani, T. J. Devries, et al. Weaning age influences the severity of gastrointestinal microbiome shifts in dairy calves. Scientific Reports, 2017, 7, 10.1038/s41598-017-00223-7. hal-01595259

HAL Id: hal-01595259

https://hal.science/hal-01595259

Submitted on 26 Sep 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Weaning age influences the severity of gastrointestinal microbiome shifts in dairy calves

Received: 1 November 2016

Accepted: 14 February 2017

Published online: 15 March 2017

S. J. Meale<fnref n="1" />, S. C. Li², P. Azevedo², H. Derakhshani©², T. J. DeVries³, J. C. Plaizier², M. A. Steele⁴ & E. Khafipour²,⁵

Ruminants microbial consortium is responsible for ruminal fermentation, a process which converts fibrous feeds unsuitable for human consumption into desirable dairy and meat products, begins to establish soon after birth. However, it undergoes a significant transition when digestion shifts from the lower intestine to ruminal fermentation. We hypothesised that delaying the transition from a high milk diet to an exclusively solid food diet (weaning) would lessen the severity of changes in the gastrointestinal microbiome during this transition. β-diversity of ruminal and faecal microbiota shifted rapidly in early-weaned calves (6 weeks), whereas, a more gradual shift was observed in late-weaned calves (8 weeks) up to weaning. Bacteroidetes and Firmicutes were the most abundant ruminal phyla in pre- and post-weaned calves, respectively. Yet, the relative abundance of these phyla remained stable in faeces (P ≥ 0.391). Inferred gene families assigned to KEGG pathways revealed an increase in ruminal carbohydrate metabolism (P ≤ 0.009) at 9, compared to 5 weeks. Conversely, carbohydrate metabolism in faeces declined (P ≤ 0.002) following a change in weaning status (i.e., the shift from pre- to post-weaning). Our results indicate weaning later facilitates a more gradual shift in microbiota and could potentially explain the negative effects of early-weaning associated with feeding a high-plane of pre-weaning nutrition.

The progressive microbial colonization of the ruminal intestinal tract begins at birth with the random acquisition of microbes from the animals' surroundings. A number of factors may influence gastrointestinal microbial colonization including the cow's vaginal microbiome, colostrum microbiota and co-inhabitance with other calves. Colostrum and milk harbor a rich microbial consortium imparting the gastrointestinal tract with beneficial microbes for metabolism and immune development, much like that in the human<fnref n="1" />. However, in young calves milk bypasses the rumen following activation of the esophageal groove²,³, resulting in the initial acquisition of nutrients from the intestinal tract.

Activation of ruminal fermentative processes commences with the introduction of solid feed into the diet, with a dramatic shift occurring when milk is completely removed from the diet (weaning), as the source from which an animal attains its nutrition shifts, greatly altering the composition of the ruminal and intestinal microbiomes⁴. The resulting ruminal microbial fermentation products are essential for the development of the rumen wall papillae, providing niche environments⁵,⁶, which aid the successive microbial colonization of the rumen with age. Early research using culture-based microbiology showed that aerobic and facultative anaerobic microbiota are gradually replaced by exclusively anaerobic taxa up to 6 and 8 weeks of age⁷,⁸. This was confirmed by recent studies characterising the natural progression of the ruminal microbiome towards that of a mature ruminant⁹, even when fed a solely milk replacer diet¹⁰, suggesting the rumen and, thus, the intestine develop temporarily to prepare for solid feed ingestion, but also indicates stimulation of ruminal development extends beyond just solid feed intake.

The stressful weaning transition has previously been linked with depressions in growth and intake. However, we recently showed that the age at which an animal is weaned can lessen these effects¹¹, possibly a result of greater

<footnote n="1">¹UMR Herbivores, INRA, Vetagro Sup, 63122, Saint-Genès-Chamanelle, France. ²Department of Animal Science, University of Manitoba, Winnipeg, MB, Canada. ³Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada. ⁴Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada. ⁵Department of Medical Microbiology, University of Manitoba, Winnipeg, MB, Canada. S. J. Meale and S. C. Li contributed equally to this work. Correspondence and requests for materials should be addressed to M.A.S. (email: masteele@ualberta.ca) or E.K. (email: Ehsan.Khafipour@umanitoba.ca)</footnote>
<table><caption>**Table 1.** Alpha diversity indices.</caption><thead><tr><th rowspan="3"></th><th colspan="9">Calf age (weeks)</th><th colspan="3">P-value</th></tr><tr><th colspan="3">Early-weaned</th><th colspan="3">Late-weaned</th><th rowspan="2">SEM</th><th rowspan="2">Weaning age</th><th rowspan="2">Calf age</th><th rowspan="2">Weaning age * Calf age</th></tr><tr><th>5</th><th>7</th><th>9</th><th>5</th><th>7</th><th>9</th></tr></thead><tbody><tr><td>Rumen</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Chao1</td><td>1065</td><td>914</td><td>883</td><td>981</td><td>1009</td><td>938</td><td>55.1</td><td>0.70</td><td>0.18</td><td>0.24</td></tr><tr><td>Shannon</td><td>4.1</td><td>3.5</td><td>3.2</td><td>3.8</td><td>3.9</td><td>3.4</td><td>0.15</td><td>0.22</td><td>&lt;0.001</td><td>0.05</td></tr><tr><td>Inverse Simpson</td><td>23.3</td><td>13.0</td><td>9.8</td><td>16.1</td><td>19.2</td><td>13.0</td><td>2.81</td><td>0.72</td><td>0.02</td><td>0.03</td></tr><tr><td>Phylogenetic diversity</td><td>48.2</td><td>35.7</td><td>34.1</td><td>45.7</td><td>45.4</td><td>38.5</td><td>2.03</td><td>0.01</td><td>&lt;0.001</td><td>0.01</td></tr><tr><td>Faeces</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Chao1</td><td>2054</td><td>2320</td><td>2601</td><td>2058</td><td>2945</td><td>2236</td><td>311</td><td>0.86</td><td>0.04</td><td>0.30</td></tr><tr><td>Shannon</td><td>4.2</td><td>5.0</td><td>5.1</td><td>3.9</td><td>4.7</td><td>4.9</td><td>0.18</td><td>0.11</td><td>&lt;0.001</td><td>0.80</td></tr><tr><td>Inverse Simpson</td><td>18.8</td><td>44.8</td><td>51.9</td><td>18.7</td><td>25.9</td><td>40.9</td><td>5.42</td><td>0.10</td><td>&lt;0.001</td><td>0.32</td></tr><tr><td>Phylogenetic diversity</td><td>52.1</td><td>66.3</td><td>69.9</td><td>46.1</td><td>62.6</td><td>64.1</td><td>4.12</td><td>0.17</td><td>&lt;0.001</td><td>0.94</td></tr></tbody></table>

ruminal and intestinal maturation at 8, compared to 6 weeks of age. Considering this, in the current study, we investigated whether weaning age is a determining factor in the extent of microbial shifts in the rumen and faeces during weaning to better understand why calf growth slows during the weaning transition. We hypothesised that earlier weaned calves (6 week of age) would show a more rapid transition of ruminal and faecal microbiota, compared to calves weaned later in life (8 weeks of age), increasing the severity of this transition, potentially explaining the greater growth depressions observed throughout weaning of calves fed a high plane of pre-weaning nutrition.

## Results

**Alpha-diversity Measures.** Richness of rumen microbiota, as indicated by the Chao1 index, was not affected by the interaction of weaning age x calf age (P = 0.24; Table 1). Additionally, the microbiota of all calves displayed similar evenness (Shannon and Inverse Simpson) across all ages, with the exception of early week 5 and late week 7 calves, which both had a greater (P ≤ 0.047) evenness than ruminal microbiota of early week 9 calves. Phylogenetic diversity (PD) of ruminal microbiota was greater (P = 0.01) in pre-weaned calves, compared to post weaned calves, with the exception of late-weaned week 5 calves which only showed a tendency (P = 0.06) to differ from late-weaned week 9 calves. Conversely, α-diversity of faecal microbiota was not affected by the interaction of weaning age x calf age or weaning age (P ≥ 0.10; Table 1). As expected, all indices were affected by calf age, where Chao1 was greater (P = 0.04) in calves at 7 weeks, compared to 5 weeks of age. Shannon, Inverse Simpson and PD were lower (P < 0.001) in faecal microbiota of calves at 5 weeks, compared to both 7 and 9 week old calves, which were similar (P > 0.05; Table 1).

**OTU Diversity and Similarity Analysis.** Community OTU comparisons were visualised by PCoA analysis (OTU ≥ 97% identity, species level similarity) using weighted (Fig. 1) and unweighted (Fig. S1) UniFrac distance metrics revealing a rapid shift in both ruminal and faecal microbiota of early-weaned calves at the time of weaning, compared to a more gradual shift in late-weaned calves (P ≤ 0.005; PERMANOVA), as evidenced by the proximity of early-weaned week 7 calves to that of post-weaned calves (early and late week 9 calves; Fig. 1), compared to the intermediate cluster of late week 7 calves between that of other pre-weaned calves (early and late week 5) and all post-weaned calves. In contrast to late-weaned calves, where faecal communities shifted (P = 0.004) between each age, faecal microbiota of early-weaned calves shifted (P < 0.002) at weaning and then remained stable (P = 0.257) in post-weaned calves regardless of their age. Interestingly, ruminal and faecal microbiota were similar across early- and late-weaned calves at 9 weeks of age (P ≥ 0.280).

**Bacterial Composition of Rumen and Faecal Microbiota.** The number of input 16S rRNA sequence reads generated from rumen and faecal microbiota were 26 772 ± 2 651 and 26 605 ± 4 503 (mean ± SD), respectively. While the majority of OTUs were identified at the genus (g.) level, some were only classified at the phylum (p.), class (c.), order (o.), or family (f.) level. Phylogenetic analysis of these sequences identified 24 phyla from the rumen, 12 of which had a relative abundance >0.1% of the community. The most abundant phyla in the rumen of pre-weaned calves was Bacteroidetes, which decreased (P ≤ 0.039) from week 5 to 7 in all calves, regardless of treatment, and also showed differences (P < 0.001) between week 5 and 9 calves, but remained stable (P ≥ 0.112) thereafter in all calves, despite their weaning status (Fig. 2; Table S1). In contrast, the relative abundance of p. Firmicutes increased (P = 0.017) in early-weaned calves from week 5 to 9 and increased (P < 0.014) in late-weaned calves after weaning, making it the most abundant ruminal phyla in all post-weaned calves (Fig. 2; Table S1). Despite different weaning ages, the magnitude of change at weaning in these two ruminal phyla was similar across all calves (P = 0.969). Additionally, the ratio of Bacteroidetes to Firmicutes in the rumen decreased (P < 0.001) by 61 vs. 67% at the respective times of weaning in early- and late-weaned calves, respectively. However, no difference (P > 0.05) between the two weaning groups was observed in this ratio at the time of weaning. Relative abundances of p. Proteobacteria, p. Actinobacteria and p. Planctomycetes were similar across all calves, regardless of weaning age or calf age, averaging 6.03%, 4.96% and 0.01% of total sequences, respectively (Fig. 2; Table S1). Several other phyla were affected (P ≤ 0.035) by weaning in early-weaned calves including
<figure><img src="image_3.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_4.png" /></figure>

**Figure 1.** Similarities of ruminal and faecal microbiota compared across weaning age [early-weaned (Wk 6) vs. late-weaned (Wk 8)] and calf age (Wk 5, Wk 7 and Wk 9). Distance between the samples, based on similarity in OTU composition (OTU similarity ≥97%) calculated using weighted UniFrac distance in the a) rumen, and b) faeces, were visualised in principal coordinates analysis (PCoA) plots. A greater distance between two points infers a lower similarity. The impact of weaning age × calf age on the clustering pattern of microbial communities was tested using PERMANOVA (implemented in PRIMER-6 software; significance was declared when P < 0.05). Within the rumen, weaning age × calf age differed (P = 0.035) across all ages and groups, with the exception (P > 0.05) of early- versus late-weaned calves at 9 weeks of age, which were similar. Weaning age × calf age also differed (P = 0.004) in the faeces, with the exception (P > 0.05) of early- versus late-weaned calves at 5 and 9 weeks of age, and in early-weaned calves aged 7, compared to 9 weeks of age.

p. Cyanobacteria, p. Spirochaetes, p. Elusimicrobia, p. Synergistetes, p. Verrucomicrobia and p. Fibrobacteres, where all except the later decreased in calves at 7, compared to 5 weeks of age. Similarly, p. Elusimicrobia, and p. Fibrobacteres decreased (P ≤ 0.02) in late-weaned calves at the time of weaning, whereas, p. Tenericutes increased (P = 0.014) post-weaning in late-weaned calves.

Weaning resulted in faecal phyla shifting from the dominant p. Firmicutes in pre-weaned calves to p. Bacteroidetes in post-weaned calves (Fig. 2; Table S2). In early-weaned calves, dominance shifted again in 9 week old calves, returning to Firmicutes as the dominant phyla. Despite these shifts in dominance, weaning age did not affect (P ≥ 0.54) faecal abundance of p. Firmicutes or p. Bacteroidetes. The ratio of Bacteroidetes to Firmicutes in the faeces showed differences to that in the rumen, where a 25% reduction was observed between week 5 and 7 for all calves (1.16 to 0.87 and 1.07 to 0.81 for early- and late-weaned calves, respectively), followed by an increase of 19 and 6% from week 7 to 9. Relative abundances of p. Tenericutes and p. Cyanobacteria were greater (P ≤ 0.002) in late week 7 and 9, compared to late week 5 calves, but showed no differences (P ≥ 0.059) in early-weaned
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Phylum level composition. Bar plots showing average relative abundance of bacterial phyla (%) in (a) rumen and (b) faeces of early- (week 6) and late-weaned (week 8) calves for each weaning age.</figcaption></figure>

calves. Comparatively, p. Spirochaetes was more abundant (P < 0.001) in late week 9 calves, compared to late week 5 and 7 calves, and was greater (P = 0.021) in early week 9 calves, compared to early week 5 calves. Similarly, p. Fusobacteria was greater (P = 0.046) in early week 9 calves, compared to early week 5 calves, yet showed an increase (P = 0.001) from week 7 to week 9 in late-weaned calves. The relative abundances of p. Proteobacteria, p. Actinobacteria, p. Elusimicrobia and p. Verrucomicrobia were similar (P ≥ 0.08) across all calves and ages.

Of the 312 genera observed in the rumen, 56 had a relative abundance over 1%, and were affected either by weaning age or calf age (Fig. 3; Table S1). The relative abundance of the dominant g. *Prevotella* remained stable (P ≥ 0.28) in the rumen of early- and late-weaned calves. The most abundant genera became *Sharpea* in both early- and late-weaned calves at 9 weeks of age following an increase (P < 0.001) as a result of weaning (week 5 vs. 7) in early-weaned calves, and an additional increase (P < 0.001) in the relative abundance from week 7 to 9 of life. Conversely, no change (P = 0.90) occurred in the relative abundance of g. *Sharpea* until weaning in late-weaned calves, where an increase was observed (P < 0.001; Fig. 3 and Table S1). The relative abundance of g. *Ruminococcus* decreased (P < 0.001) at weaning in early-weaned calves and remained stable (P ≥ 0.64) thereafter, and only showed a reduction (P < 0.001) between week 5 and week 9 late-weaned calves. The time of weaning did not affect the size of this shift (P = 1.00; Table S1). Additionally, weaning at 6 weeks (early-weaning) caused
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Correlation between ruminal microbiota and calf age and weaning strategy. The heatmap shows bacterial abundance in rumen varied with calf age (5, 7 and 9 weeks) and weaning age (6 vs. 8 weeks for early- and late-weaning, respectively). Data represent taxa present at greater than 1% of community that were significantly different (P < 0.05) between at least two ages or between treatment groups (weaning age). See Table S1 for additional data. The OTU count data was subjected to variance stabilizing transformation and summarised to calculate genus abundant data. Hierarchically clustering was then performed based on complete linkage of Euclidean distances matrix and visualized over heatmap using the heatmap.2 function of the "gplots" R package.</figcaption></figure>

a sharp increase (P < 0.001) in the relative abundance of g. *Shuttleworthia* at weaning, remaining stable there-
after, whereas late-weaned calves exhibited only a tendency (P = 0.06) to increase from week 5 to 9 (Fig. 3 and
Table S1).

A total of 190 genera were identified in faecal microbiota, 46 of which were present at >1% of the total
sequences and showed a difference relative to either weaning age of calf age (Fig. 4; Table S2). The relative abun-
dance of g. *Bacteroides* was dominant and stable (P = 1.00; Table S2) across all ages regardless of weaning age,
averaging 25.1% of total sequences. *Blautia* was the second most abundant genera in pre-weaned calves, regard-
less of weaning age. Relative abundance of this genera declined (P < 0.001) at weaning; however, the magnitude
of this decrease was similar despite different weaning ages (av. 84.8% decrease; P = 1.00; Fig. 4 and Table S2).
Conversely, an unclassified genus in f. *Ruminococcaceae* was the second most abundant taxa in post-weaned
calves following similar increases (P ≤ 0.036) at weaning, regardless of weaning age.

**Correlation between Rumen and Faecal Microbiota and Calves Physiological Variables.** Of the genera detected in ruminal and faecal microbiota, the relative abundance of those present in at least 50% of calves and represented at least 0.1% of the bacterial community in at least one calf, were compared to calf physiological variables and rumen fermentation characteristics (Fig. 5). The relative abundance of g. *Prevotella* in the rumen was negatively correlated (P ≤ 0.021) with calf body weight (Spearman’ ρ = −0.53), starter intake (Spearman’ ρ = −0.77) and serum βHBA (Spearman’ ρ = −0.81) and positively correlated with ruminal pH (Spearman’ ρ = 0.49; Fig. 5a). Conversely, in the faeces, g. *Prevotella* was positively correlated (Spearman’ ρ = 0.56; P = 0.001) with starter intake, but showed no correlation with body weight (Fig. 5b). Similarly, g. *Dialister* in both rumen and faeces were positively correlated (P ≤ 0.002) with calf physiological variables including body weight (Spearman’ ρ = 0.55; Spearman’ ρ = 0.55) and starter intake (Spearman’ ρ = 0.67; Spearman’ ρ = 0.66), respectively, and serum βHBA (Spearman’ ρ = 0.63) with ruminal g. *Dialister*.

Similarly, g. *Sharpea* in the rumen was positively correlated (P ≤ 0.001) with body weight (Spearman’ ρ = 0.59),
starter intake (Spearman’ ρ = 0.66) and serum βHBA (Spearman’ ρ = 0.64) and, in addition, showed a positive
correlation (P = 0.03) with forage intake (Spearman’ ρ = 0.48). Conversely, no correlation (P > 0.05) between the
relative abundance of g. *Sharpea* in faeces and calf physiological variables was observed.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Correlation between faecal microbiota and calf age and weaning strategy. The heatmap shows bacterial abundance in faeces varied with calf age (5, 7 and 9 weeks) and weaning age (6 vs. 8 weeks for early- and late-weaning, respectively). Data presented represent taxa present at greater than 1% of community that were significantly different (*P* < 0.05) between at least two ages or between treatment groups (weaning age). See Table S2 for additional data. The OTU count data was subjected to variance stabilising transformation and summarised to calculate genus abundant data. Hierarchically clustering was then performed based on complete linkage of Euclidean distances matrix and visualised over heatmap using the heatmap.2 function of the "gplots" R package.</figcaption></figure>

**Metabolic Functions and Capacity of Ruminal and Faecal Microbiomes.** To assess the metabolic potentials of ruminal and faecal microbiomes, OTUs were entered into PICRUSt and the inferred gene families were annotated against KOs then collapsed into KEGG pathways. Gene families annotated to bacterial carbohydrate metabolism increased (*P* ≤ 0.009) in the rumen of all calves at 9 weeks of age, compared to 5 weeks of age (Fig. 6; Table S3). Conversely, the number of gene families annotated to bacterial carbohydrate metabolism in faeces only declined (*P* ≤ 0.002) following a change in weaning status (i.e., pre- vs. post-weaning) in both early- and late-weaned calves (Fig. 6; Table S4). The relative abundance of KEGG L2 pathways associated with metabolism of amino acids, energy or lipids remained stable (*P* ≥ 0.114) in ruminal bacteria across all calf ages, regardless of weaning age. Similarly, the abundance of genes associated with lipid metabolism was stable (*P* ≥ 0.201) in the faeces of all calves. However, the abundance of genes associated with both energy and amino acid metabolism tended to be lower (*P* ≤ 0.059) in faeces of early-weaned week 7, compared to early-weaned week 5 calves. Additionally, genes annotated to amino acid metabolism showed a tendency to decrease (*P* = 0.094), and decreased (*P* ≤ 0.036) in faeces of calves at 9 weeks, compared to those at 5 weeks of age, across early- and late-weaned calves, respectively.

A total of 328 KEGG level 3 pathways were identified in the rumen and 266 from faeces (Tables S4 and S5). Of the OTU inferred gene families assigned to KEGG pathways in the rumen, the three most abundant related to carbohydrate metabolism that changed in both early- and late-weaned, included amino sugar and nucleotide sugar metabolism (hsa00520), glycolysis/gluconeogenesis (hsa00010), fructose and mannose metabolism (hsa00051; Fig. 6). Fructose and mannose metabolism (*P* ≤ 0.002) increased following weaning in early-weaned calves. However, only a difference was observed between calves at week 5 and week 9 in late-weaned calves. Glycolysis/ gluconeogenesis and amino sugar and nucleotide sugar metabolism were greater (*P* ≤ 0.01) in calves at 9 weeks of age, compared to 5 or 7 weeks for late-weaned calves. Yet, early-weaned calves at 9 week of age showed a greater (*P* ≤ 0.014) abundance of OTUs assigned to glycolysis/gluconeogenesis and amino sugar and nucleotide sugar metabolism, compared to those aged 5 weeks, with those aged 7 weeks intermediary between week 5 and week 9.

Comparatively, in the faecal microbiome, of the inferred gene families assigned to KEGG pathways, the three most abundant related to carbohydrate metabolism were the same as that of the rumen, with amino sugar and nucleotide sugar metabolism, and fructose and mannose metabolism decreasing (*P* < 0.006) following the respective weaning of both early- and late-weaned calves (Fig. 6). Glycolysis/gluconeogenesis showed a tendency to
a

<figure><img src="image_8.png" /></figure>

b

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Correlations between production variables and relative taxa abundance. Spearman non-parametric rank correlation matrix of the dominant bacterial genera across (a) rumen and (b) faecal samples. The genera were included in the matrix if they were in at least 50% of calves and represented at least 0.1% of the bacterial community in at least one calf. All correlations presented were statistically significant (P < 0.05), with strong correlations indicated by large circles and weaker correlations indicated by small circles. The scale colours denote whether the correlation is positive (closer to 1, blue circles) or negative (closer to -1, red circles) between the taxa and production variables. The OTU count data was subjected to variance stabilising transformation then pairwise Spearman correlations between bacterial genus and biological parameters at corresponding ages were analysed. P-values were adjusted for multiple comparisons.</figcaption></figure>

decrease following weaning in early-weaned calves (P ≤ 0.082), but remained stable (P ≥ 0.482) in late-weaned calves (Fig. 6).

Discussion

Recently, we showed that weaning age (6 vs. 8 weeks of age) differentially affects physiological development of a calf fed a high plane of pre-weaning nutrition (>8 L/d), including growth during the transitional period, where later weaning promoted beneficial effects on growth and intake during weaning and thereafter¹¹. Early studies suggested that the natural transition to strictly anaerobic taxa in the rumen ceases at around 6 to 8 weeks of age indicating that maturity of the ruminal and intestinal microbiomes at weaning may influence the impact this transition has on calf development and growth. Subsequently, in this study we sought to characterise the shift in ruminal and faecal microbiota, representative of the intestinal microbiome, which occurs across weaning in calves weaned at 6 and 8 weeks of age to elucidate if changes in the microbial structure could be responsible for the observed production effects (reported previously¹¹).

A comparison of all calves showed that species richness (chao1 index) within the rumen and faeces were con-
stant across all ages, regardless of calf weaning age (Table. 1). However, the microbial community in the rumen
of calves weaned at an earlier age became less even, as determined by the Shannon index, between week 5 and 9
of age, a change which was not observed in late-weaned calves. Similarly, weaning calves at 6 weeks of age caused
a rapid shift in β-diversity of rumen and faeces from a pre- to post-weaned state, evidenced by the similar diver-
sity in those calves at week 7 and 9 (Fig. 1). Such a shift reflects the abrupt change in the source from which the
calves obtained nutrients and was perhaps related to the pre-mature development of the rumen, paralleled by a
<figure><img src="image_10.png" /><figcaption>**Figure 6.** PICRUSt predicted summary of COG categories from 16S amplicon sequencing of (a) ruminal and (b) faecal microbiomes. Relative abundances of Level 2 KEGG pathways are depicted by calf age within treatment group (weaning age) in circular barplots. Relative abundances of the three most abundant L3 KEGG pathways involved in Carbohydrate Metabolism (Fructose and Mannose Metabolism, Glycose/gluconeogenesis and, Amino sugar and nucleotide sugar Metabolism) at different ages (5, 7, or 9) subjected to weaning at either 6 weeks (Early) or 8 weeks (Late) of age are also presented within the boxes. \*P < 0.05 for the shown comparisons of L3 pathways (See Supplementary Tables 3–6 for complete KEGG data).</figcaption></figure>

reduction in calf growth (44 vs. 12% reduction for early- and late-weaned calves, respectively) observed by these calves during weaning¹¹. Conversely, late-weaned calves displayed a gradual transition in β-diversity (Fig. 1) with calves at week 7 displaying an intermediate diversity between calves aged 5 and 9 weeks, reflecting the gradual increase in consumption of starter (Fig. S2) as calves aged and thus, a progressive development of the rumen. Interestingly, all calves at 9 weeks of age displayed similar diversity when using weighted UniFrac measurement; however, a difference was observed in the rumen of early- vs. late-weaned calves at 9 weeks of age when considering unweighted UniFrac (Fig. S1), which accounts for differences in unique evolutionary history of either group. This may indicate that although the microbial profile was similar, its functional potential may have evolved differently as a result of different weaning ages.

Weaning of calves shifted the dominant ruminal phyla from Bacteroidetes in pre-weaned calves, to Firmicutes in post-weaned calves, regardless of their age at weaning (Fig. 2). The abundance of p. Bacteroidetes; however, only declined between week 5 and 7, and remained stable thereafter. This corresponds to 2 month old pre-weaned calves fed milk replacer with access to calf starter where Bacteroidetes was also the dominant phyla⁹. Conversely, the increase in p. Firmicutes at the respective weaning times in our study, agrees with the increase observed between the abundance of p. Firmicutes in weaned-calves at 6 months of age, compared to pre-weaned calves at 2 months of age⁹. Similarly, Li et al.¹⁰ observed a decline in the abundance of Bacteroidetes in pre-weaned calves to that of 1 year old post-weaned bull calves fed a hay diet and a corresponding increase in p. Firmicutes. These results, in conjunction with the abundance of p. Proteobacteria, concur with previous studies in cattle⁴, ⁹ which suggest that the developing rumen (pre-weaning) comprises the same three dominant phylum (Bacteroidetes, Firmicutes and Proteobacteria) as the more mature gut of post-weaned calves, although the relative abundances of each phyla vary with the stage of development¹².

The dominance of g. *Prevotella* in the rumen of all calves was not unexpected, as this genus comprises species involved in degradation of non-structural carbohydrates and xylan¹³. Additionally, the stability of this genus across weaning is in agreement with Rey *et al.*¹⁴, where the Bacteroidetes phylum was dominated by *Prevotella* in all calves older than 15 d of age when consumption of starter increased above 100 g/d, resulting in the acidification of the rumen¹⁵. Similarly, in our study, calves began consuming starter at 2 weeks of age with intakes averaging 83 g/d, and reaching 1307 and 504 g/d at week 7 for early- and late-weaned calves, respectively. Similarly, our results agree with those obtained in steers¹⁰, ¹⁶ and in adult dairy cows¹⁶, ¹⁷ suggesting that this genus is associated with diets containing solid food. Conversely, Li *et al.*¹⁸ and Wu *et al.*¹⁷ noted in calves fed solely milk
replacer at 42 d of age that g. *Bacteroides* was dominant within the rumen, yet g. *Prevotella* was only present at low abundance rates. In accordance, within our study we observed a negative correlation between g. *Bacteroidetes* and starter intake, where g. *Bacteroides* decreased from week 5 to 7 in all calves, coinciding with the timing of starter consumption reaching above 100 g/d. However, calves weaned earlier displayed a greater reduction in g. *Bacteroides* at the time of weaning (6 weeks of age), compared to later-weaned calves (8 weeks of age), likely reflecting the greater shift in dietary composition from high milk replacer and low starter intakes to solely starter in early-weaned calves, compared to a more moderate starter and milk replacer diet consumed by late-weaned calves. Two additional ruminal genera which showed contrasting changes depending on the age at weaning included *Dialister* and *Shuttleworthia*, both of which increased dramatically in early-weaned calves, but showed no differences in late-weaned calves across weaning (Fig. 3). Collectively, these results indicate that, in agreement with earlier studies¹⁹, ²⁰, the sooner a calf consumes solid feed, the sooner a ruminal bacterial community more representative of a mature rumen develops. However, the rumen may not be as capable of handling such dramatic changes in calves at 6 weeks of age, and could possibly account for the negative production effects observed in these calves.

The abundance of the two dominant faecal phyla Firmicutes and Bacteroidetes remained stable across weaning in all calves. This corresponds to a recent study in humans²¹, which explored the faecal microbiome as solid feed was introduced into the diet of babies fed breast milk or formula. The relative abundance of p. Firmicutes remained stable whereas, p. Bacteroidetes increased with solid feed consumption, as also reported by Koenig *et al.*²². As the calves in the current study were already consuming solid feeds at the first sampling (64 vs. 103 g/d for early- and late-weaned calves, respectively) it is reasonable to conclude the shift in p. Bacteroidetes may have already occurred. Additionally, it has been suggested that the dietary balance between carbohydrates and protein or animal fat is the primary factor shifting faecal microbiota, where greater g. *Bacteroides* abundance is associated with high protein and animal fat diets, compared to g. *Prevotella* which predominates with greater carbohydrate intake¹⁹, ²¹, ²³. In our study, the dominant genera within the faeces was *Bacteroides* which also remained stable across all calves. Additionally, g. *Blautia* from the f. Lachnospiraceae, which constitutes a major taxonomic group capable of degrading complex polysaccharides to short chain fatty acids including acetate, butyrate, and propionate in the human gut²⁴, ²⁵ was observed to decrease post-weaning in our study, compared to an increase observed following solid feed consumption in human babies²¹, likely reflecting the different digestive systems and the shift to ruminal fermentation in post-weaned calves.

The observed shifts in bacterial community have previously been noted to correspond to the nature of the fermentation substrate more so than the amount consumed¹⁴, with increases in solid feed correlating to increased production of short chain fatty acids (SCFA)²¹, ²². However, here we observed that despite being a producer of SCFA, ruminal abundance of g. *Bacteroides* was negatively correlated with acetate and butyrate concentration as well as calf body weight, starter intake and serum βHBA; however, no correlations were observed between *Bacteroides* and propionate concentration in the rumen. We observed a co-linearity between starter intake and water intake (R² = 0.90), and starter intake and serum βHBA (R² = 0.85; data not shown) in calves across all ages, providing a possible explanation for the strong association between bacterial taxa abundance and serum βHBA levels, and further suggests that in this study, the abundance of g. *Bacteroidetes* strongly relates to starter consumption. Furthermore, increasing serum βHBA levels in relation to increasing starter consumption is not surprising, as the end products of ruminal fermentation, SCFA, are absorbed and metabolised by ketogenesis in ruminal epithelium²⁶.

Within the Firmicutes phyla, g. *Sharpea* increased in agreement with our previous work⁴ and showed a positive correlation to starter intake and serum βHBA (Fig. 5). A strong association between the ruminal abundance of *Sharpea* and methane was found recently, where lower emitting sheep had a higher abundance of g. *Sharpea* and an increased lactate production, which in turn, resulted in an increased conversion of lactate to butyrate by *Megasphaera* spp.²⁷. Whist an association was not observed between ruminal *Sharpea* and fermentation parameters here, it is possible that such an increase in butyrate contributed to the observed overall increase in serum βHBA levels.

Similarly, ruminal and faecal abundances of g. *Dialister* were positively correlated with starter intake and calf body weight in the current study, corresponding to a study by Meyer *et al.*²⁸ which suggested that ruminal abundance of this genera were positively associated with average daily gain in beef steers, suggesting that perhaps this genera is responsive to body weight. Additionally, g. *Dialister* has been positively correlated with total VFAs and propionate, with the authors indicating this genera may be involved in VFA metabolism and starch degradation due to observed increases in amylase and carboxymethylcellulase activities²⁹.

In addition to the observed changes in bacterial abundances, diet induced differences may also contribute to host metabolism and immune function by regulating genes involved in lipid and carbohydrate metabolism, altering endocrine functions, increasing inflammatory responses, and influencing energy balance and body weight³⁰. Using predicted metagenomes (PICRUSt) of rumen and faeces microbiomes we inferred that there were several functional pathways associated with weaning (Fig. 6). Most notably, a significant reduction in genes associated with carbohydrate metabolism in the faecal microbiome following weaning mirrors the shift in nutrient metabolism away from the lower GI tract to ruminal fermentation, as previously seen in Meale *et al.*⁴. However, increases in the abundance of genes associated with carbohydrate metabolism in the rumen only occurred in 9-week old calves, compared to those at 5 weeks, rather than with changes in weaning status. This suggests, in agreement with our previous study⁴, that consumption of starter initiated a shift in microbial taxa towards that of the mature rumen and potentially lessened any effects of a shift at weaning. Li *et al.*¹⁰ noted that in the absence of solid feed, carbohydrate transport and metabolism, as indicated by the Clusters of Orthologous Groups (COG) classification, did not change with increasing age in pre-weaned calves suggesting that ingestion of solid feed is essential in developing the functionality of the rumen. The increase observed in calves at 9 weeks of age, indicates that perhaps the functional role of the rumen was still developing. It is possible that at this age the rumen is better
developed to handle increasing consumption of starter and can regulate carbohydrate metabolism accordingly, providing some insights as to why calves weaned at 6 weeks of age exhibited a greater decline in growth through the weaning transition. Importantly, however, we must note that PICRUSt predictions are based only on known functions of the microbial communities present in the GI tract of humans and animals obtained from the whole genome shotgun sequencing of such samples. As there are limited numbers of shotgun sequencing studies in ruminants, the functionality of the members of ruminant GI microbiota may be over- or underestimated in PICRUSt.

In conclusion, our results suggest that weaning later results in a more gradual shift in microbial diversity in both the rumen and faeces of dairy calves and this shift, in conjunction with that observed in functional data, indicates an increased state of readiness for the weaning transition in calves weaned at 8 weeks of age; this may account for the lower declines in growth observed in these calves reported earlier¹¹. Despite no changes in the relative abundance of dominant intestinal phyla, as represented by the faeces, across weaning at either weaning age, the gradual increase in solid feed consumption of late-weaned calves resulted in the rumen harbouring a microbial consortium sufficient to perform the major fermentative and metabolic functions prior to weaning, and as such, the level of physiological stress associated with the weaning transition was not as severe as observed in early-weaned calves.

## Methods

**Animal Experiment and Sample Collection.** Experimental procedures used in this study were approved by the Nutreco Canada Agresearch Animal Care Committee in accordance with the Canadian Council on Animal Care guidelines³¹. A detailed description of the experimental treatments and growth measurements was published previously¹¹. Briefly, 20 female Holstein dairy calves were randomly assigned at birth to be weaned at 6 (early) or 8 (late) weeks of age. Calves were fed Shur-Gain Optivia High Performance (26% CP, 16% fat; ME = 4.58 Mcal/kg) Milk Replacer (Nutreco Canada Inc., Guelph, ON, Canada), mixed at a rate of 150 g/L of 40 °C water twice daily (0730 and 1630 h) by nipple bottle (Super Calf Nipple). At d 7 ± 3 calves were transitioned to a gate-mounted artificial teat (Peach Teats; Skellerup Industries Ltd., Woolston, New Zealand), through which they were fed until weaning. Calves were offered 6 L from d 1 to 3, 7 L from d 4 to 6, and 8 L from d 7 until their respected time of step-down weaning. The early-weaning group was stepped down to 4 L/d on d 36 and weaned on d 43. The late-weaning group was stepped down to 4 L/d on d 50 and weaned on d 57 [see Eckert et al.¹¹ for figure]. Calves had *ad libitum* access to water, starter (Optivia Advantage Calf Starter; 22% CP; ME = 2.63 Mcal/kg), chopped straw (8% CP) and oat straw (8% CP; ME = 1.55 Mcal/kg) chopped to a length of 3 cm, from d 0–70. Orts were weighed and discarded daily.

Rumen and faecal samples were collected on d 35, 49 and 63 (Week 5, 7 and 9) of life. Rumen fluid was collected, on those same days, at 1100 h using a modified Geishauser oral probe (1.3 cm diameter)¹¹, ³². The probe was inserted 50 cm inside the calf and pH measurements of all samples were above 5.0, indicative of the rumen not the abomasum. Samples were collected in 15 mL tubes, flash frozen in liquid nitrogen and stored at -80 °C. Calves were rectally finger-stimulated with sterile-gloved hand to facilitate the collection of a 50 g faecal sample, which was immediately frozen at -80 °C.

Calf body weight (BW), starter intake, straw intake, faecal starch concentrations, serum β-hydroxybutyrate (βHBA), behavioural observations, rumen acetate, butyrate and propionate, and total rumen volatile fatty acids (VFA) concentrations, reported previously¹¹, were reinterpreted in the context of their correlations with microbial changes evaluated in this research.

**DNA Extraction.** Rumen fluid and faecal samples were thawed at room temperature and kept on ice during the extraction process. Two hundred milligram faecal samples or the sediment collected from 1 mL rumen fluid by centrifuging at 15,000 × g for 5 min was used for DNA extraction using a ZR faecal DNA kit (D6010; Zymo Research Corp., Irvine, CA, USA) that included a bead-beating step for mechanical disruption of microbial cells. DNA was eluted from the column with elution buffer, and DNA concentration was quantified using a NanoDrop 2000 spectrophotometer (Thermo Scientific, Waltham, MA, USA). DNA samples were normalized to 50 ng/µL, and quality verified by PCR amplification of the 16S rRNA gene using universal primers 27 F (5'-GAAGAGTTTGATCATGGCTCAG-3') and 342 R (5'-CTGCTGCCTCCCGTAG-3'), as described before³³. Amplicons were verified by agarose gel electrophoresis. All DNA samples were stored at -80 °C.

**Library Construction and Illumina Sequencing.** The V4 region of 16S rRNA gene was targeted for PCR amplification using modified F515/R806 primers³⁴ as described before³⁵. The reverse PCR primer was indexed with 12-base Golay barcodes allowing for multiplexing of samples. PCR reaction for each sample was performed in duplicate and contained 1.0 µL of pre-normalized DNA, 1.0 µL of each forward and reverse primers (10 µM), 12 µL HPLC grade water (Fisher Scientific, Ottawa, ON, Canada) and 10 µL 5 Prime Hot MasterMix (5 Prime, Inc., Gaithersburg, USA). Reactions consisted of an initial denaturing step at 94 °C for 3 min followed by 35 amplification cycles at 94 °C for 45 sec, 50 °C for 60 sec, and 72 °C for 90 sec; finalized by an extension step at 72 °C for 10 min in an Eppendorf Mastercycler pro (Eppendorf, Hamburg, Germany). PCR products were then purified using ZR-96 DNA Clean-up Kit (ZYMO Research, Irvine, CA, USA) to remove primers, dNTPs and reaction components. The V4 library was generated by pooling 200 ng of each sample, quantified by Picogreen dsDNA (Invitrogen, Grand Island, NY, USA). This was followed by multiple dilution steps using pre-chilled hybridization buffer (HT1; Illumina, San Diego, CA, USA) to bring the pooled amplicons to a final concentration of 5 pM, measured by Qubit 2.0 Fluorometer (Life Technologies, Burlington, ON, Canada). Finally, 15% of PhiX control library was spiked into the amplicon pool to improve the unbalanced and biased base composition, a known characteristic of low diversity 16S rRNA libraries. Customized sequencing primers for read1 (5'-TATGGTAATTGTGTGCCAGCMGCCGCGGTAA
-3'), read2 (5'-AGTCAGTCAGCCGGACTACHVGGGTWTCTAAT-3') and index read (5'-ATTAGAWACCCBDGTAGTCCGGCTGACTGACT-3') were synthesized and purified by polyacrylamide gel electrophoresis (Integrated DNA Technologies, Coralville, IA, USA) and added to the MiSeq Reagent Kit V2 (300-cycle; Illumina, San Diego, CA, USA). The 150 paired-end sequencing reaction was performed on a MiSeq platform (Illumina, San Diego, CA, USA) at the Gut Microbiome Laboratory (Department of Animal Science, University of Manitoba, Winnipeg, Canada). The sequencing data were deposited into the Sequence Read Archive (SRA) of NCBI (http://www.ncbi.nlm.nih.gov/sra) and can be accessed via accession number SRR4733877.

**Bioinformatic Analysis.** The PANDAseq assembler³⁶ was used to merge and fix the overlapping paired-end Illumina fastq files. All the sequences with low quality base calling scores, as well as those containing uncalled bases (N) in the overlapping region were discarded. The output fastq file was then analysed by downstream computational pipelines of the open source software package QIIME³⁷. The default minimum quality threshold of 25 was used. Chimeric sequences were detected using the UCHIME algorithm (USEARCH 6.1) to run both *de novo* and reference based chimera detection. The number of chimeric sequences identified and consequently removed by both detection methods was 1.7% of total high quality sequences. Sequences were clustered at the 97% sequence similarity level using the Greengenes database (Version 13.5)³⁸ using an open reference-based OTU picking approach with the QIIME algorithm and usearch61 method with default parameters³⁸, ³⁹. Those sequences that failed to cluster were subsampled for *de novo* OTU picking. All picked OTUs were subsequently aligned by PyNAST⁴⁰ and a phylogenetic tree was built using FastTree method⁴¹ to calculate UniFrac distances⁴² within QIIME. Representative OTUs were assigned to bacterial taxonomies using RDP classifier via QIIME with a confidence threshold of 0.8⁴³. Finally, open source software PICRUSt (phylogenetic investigation of communities by reconstruction of unobserved states)⁴⁴ was used on 16S rRNA gene sequencing data to predict functional genes of the classified members of the rumen and faecal microbiota resulting from reference-based OTU picking against Greengenes database. Predicted genes were then hierarchically clustered and categorized under Kyoto Encyclopaedia of Genes and Genomes (KEGG)⁴⁵ orthologs (KOs) and pathways (level 3; Tables S4 and S6).

**Alpha- and Beta-Diversity Analyses.** The phylogenetic tree was rooted using *Methanococcus jannaschii* (L77117) as an outgroup. Subsequently an OTU table was generated by QIIME, which along with the mapping file was assembled into a Phyloseq object⁴⁶. For within community diversity (α-diversity) calculations, the same number of sequences for each sample were randomly selected (corresponding to the sample with lowest number of sequences being 14,148 and 20,031 for rumen and faecal samples, respectively), to eliminate the bias caused by different sample sizes⁴⁷. Alpha-diversity analysis was conducted with standard diversity metrics accessed via Phyloseq, including observed richness, Chao1, Shannon index, Simpson index, Inverse Simpson index and Phylogenetic Diversity (PD). The dataset was also subsampled to the minimum⁴⁸ to compare microbial compositions between samples (β-diversity). Beta-diversity was measured by calculating the weighted and unweighted UniFrac distances⁴⁹ using Phyloseq default scripts. Principal coordinate analysis (PCoA) was applied on the resulting distance matrices to generate two-dimensional plots using PRIMER v6 software⁵⁰. Permutational multivariate analysis of variance (PERMANOVA)⁵¹ was used to calculate P-values and to test differences of β-diversity among treatment groups for significance. Both weighted and unweighted UniFrac distances were used to compute the test statistic and P-values⁵².

**Bacterial Community Composition and Metagenome prediction.** Bacterial community composition at the phylum and genus levels, and predictive metagenome profiles were compared among treatments (early vs. late weaning) and weeks (5, 7, and 9) using the nbinomWaldTest method of DESeq²⁵³ according to a randomized design with pre-weaning as covariate and containing dietary treatment, day and calf; these were considered significant at P < 0.05 (DESeq2, R package version 1.8.1, 2014). All P-values were corrected for false discovery rates using R software package (version 3.3.1, Vienna, Austria).

**Correlations between Calf Physiological Variables and Bacterial Abundance.** Non-parametric Spearman rank correlation coefficient analysis implemented in PAST software⁵⁴ was used to test the relationship between BW, average daily gain, starter intake, forage intake, water intake, metabolizable energy intake, ruminal pH, faecal starch, rumen acetate, rumen butyrate, rumen propionate, total rumen VFA and the bacterial communities in rumen fluid and faeces. Physiological data was published previously (Eckert *et al*.¹¹). Data from the three corresponding ages (5, 7 and 9 weeks of age) of microbial sampling were used to produce the correlation matrix. The resulting correlation matrix was visualized in a heatmap format generated by the corrplot package of R [Corrplot: visualization of a correlation matrix; R package version 02-0. 2010 (http://CRAN)].

**References**

1. Jeurink, P. V. et al. Human milk: a source of more life than we imagine. *Benef Microbes* 4, 17–30, doi:10.3920/BM2012.0040 (2013).

2. Black, J. L. & Sharkey, M. J. Reticular groove (sulcis reticuli): an obligatory adaptation in ruminant-like herbivores. *Mammalia* 34, 294–302, doi:10.1515/mamm.1970.34.2.294 (1970).

3. Church, D. C. *Digestive Physiology and Nutrition of Ruminants*. 2 edn, Vol. 1 (Corvallis OR: O&B Books, 1976).

4. Meale, S. J. et al. Development of ruminal and fecal microbiomes are affected by weaning but not weaning strategy in dairy calves. *Front Microbiol* 7, doi:10.3389/fmicb.2016.00582 (2016).

5. Klein, R. D. et al. Dietary fiber and early weaning on growth and rumen development of calves. *J Dairy Sci* 70, 2095–2104 (1987).

6. Beharka, A. A., Nagaraja, T. G., Morrill, J. L., Kennedy, G. A. & Klemm, R. D. Effects of form of the diet on anatomical, microbial, and fermentative development of the rumen of neonatal calves. *J Dairy Sci* 81, 1946–1955, doi:10.3168/jds.S0022-0302(98)75768-6 (1998).

7. Fonty, G., Gouet, P., Jouany, J. P. & Senaud, J. Establishment of the Microflora and Anaerobic Fungi in the Rumen of Lambs. *J Gen Microbiol* 133, 1835–1843 (1987).

8. Minato, H., Otsuka, M., Shirasaka, S., Itabashi, H. & Mitsumori, M. Colonization of Microorganisms in the Rumen of Young Calves. *J Gen Appl Microbiol* 38, 447–456, doi:10.2323/Jgam.38.447 (1992).
9. Jami, E., Israel, A., Kotser, A. & Mizrahi, I. Exploring the bovine rumen bacterial community from birth to adulthood. *ISME J* 7, 1069–1079, doi:10.1038/ismej.2013.2 (2013).

10. Li, R. W., Connor, E. E., Li, C., Baldwin, V. I. R. L. & Sparks, M. E. Characterization of the rumen microbiota of pre-ruminant calves using metagenomic tools. *Environ Microbiol* 14, 129–139, doi:10.1111/j.1462-2920.2011.02543.x (2012).

11. Eckert, E., Brown, H. E., Leslie, K. E., DeVries, T. J. & Steele, M. A. Weaning age affects growth, feed intake, gastrointestinal development, and behavior in Holstein calves fed an elevated plane of nutrition during the preweaning stage. *J Dairy Sci* 98, 6315–6326 (2015).

12. Jami, E. & Mizrahi, I. Composition and Similarity of Bovine Rumen Microbiota across Individual Animals. *PLoS One* 7, e33306, doi:10.1371/journal.pone.0033306 (2012).

13. Belanche, A. *et al.* Shifts in the rumen microbiota due to the type of carbohydrate and level of protein ingested by dairy cattle are associated with changes in rumen fermentation. *J Nutr* 142, 1684–1692, doi:10.3945/jn.112.159574 (2012).

14. Rey, M. *et al.* Establishment of ruminal bacterial community in dairy calves from birth to weaning is sequential. *J Appl Microbiol* 116, 245–257, doi:10.1111/jam.12405 (2014).

15. Rey, M., Enjalbert, F. & Monteils, V. Establishment of ruminal enzyme activities and fermentation capacity in dairy calves from birth through weaning. *J Dairy Sci* 95, 1500–1512 http://dx.doi.org/10.3168/jds.2011-4902 (2012).

16. Pitta, D. W. *et al.* Rumen bacterial diversity dynamics associated with changing from bermudagrass hay to grazed winter wheat diets. *Microb Ecol* 59, doi:10.1007/s00248-009-9609-6 (2010).

17. Wu, S. *et al.* The Bacterial Community Composition of the Bovine Rumen Detected Using Pyrosequencing of 16S rRNA Genes. *Metagenomics* 1, 1–11 (2012).

18. Li, R. W. In *Metagenomics and Its Applications in Agriculture, Biomedicine, and Environmental Studies* (ed. R. W., Li) (Nova Science Publisher's, 2011).

19. Hungate, R. E. (Academic Press, New York, 1966).

20. Biesheuvel, M. H., Bijker, P. G. H. & Urlings, H. A. P. Some Aspects of the Gastrointestinal Microflora of Veal Calves Fed Different Rations - a Pilot-Study. *Vet Quart* 13, 97–104 (1991).

21. Thompson, A. L., Monteagudo-Mera, A., Cadenas, M. B., Lampl, M. L. & Azcarate-Peril, M. A. Milk- and solid-feeding practices and daycare attendance are associated with differences in bacterial diversity, predominant communities, and metabolic and immune function of the infant gut microbiome. *Front Cell Infect Microbiol* 5, doi:10.3389/fcimb.2015.00003 (2015).

22. Koenig, K. M. *et al.* Effect of dietary *Enterolobium cyclocarpum* on microbial protein flow and nutrient digestibility in sheep maintained fauna-free, with total mixed fauna or with *Entodinium caudatum* monofauna. *Brit J Nutr* 98, 504–516 (2007).

23. De Filippo, C. *et al.* Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. *PNAS* 107, 14691–14696, doi:10.1073/pnas.1005963107 (2010).

24. Biddle, A., Stewart, L., Blanchard, J. & Leschine, S. Untangling the genetic basis of fibrolytic specialization by *Lachnospiraceae* and *Ruminococcaceae* in diverse gut communities. *Diversity* 5, doi:10.3390/d5030627 (2013).

25. Eren, A. M. *et al.* A single genus in the gut microbiome reflects host preference and specificity. *ISME J* 9, doi:10.1038/ismej.2014.97 (2015).

26. Baldwin, R. L., McLeod, K. R., Klotz, J. L. & Heitmann, R. N. Rumen development, intestinal growth and hepatic metabolism in the pre- and postweaning ruminant. *Journal of Dairy Science* 87, E55–E65 (2004).

27. Kamke, J. *et al.* Rumen metagenome and metatranscriptome analyses of low methane yield sheep reveals a Sharpea-enriched microbiome characterised by lactic acid formation and utilisation. *Microbiome* 4, 56, doi:10.1186/s40168-016-0201-2 (2016).

28. Myer, P. R., Smith, T. P. L., Wells, J. E., Kuehn, L. A. & Freetly, H. C. Rumen Microbiome from Steers Differing in Feed Efficiency. *PLoS One* 10, e0129174, doi:10.1371/journal.pone.0129174 (2015).

29. Wang, W. *et al.* Effects of early feeding on the host rumen transcriptome and bacterial diversity in lambs. *Scientific Reports* 6, 32479, doi:10.1038/srep32479 http://www.nature.com/articles/srep32479#supplementary-information (2016).

30. Gill, S. R. *et al.* Metagenomic analysis of the human distal gut microbiome. *Science* 312, doi:10.1126/science.1124234 (2006).

31. Canadian Grain Commission. Official Grain Grading Guide (ISSN 1704-5118, 2013).

32. Geishauser, T. An instrument for the collection and transfer of ruminal fluid and for the administration of water soluble drugs in adult cattle. *Bovine Practice* 27, 38–42 (1993).

33. Khafipour, E., Li, S., Plaizier, J. C. & Krause, D. O. Rumen microbiome composition determined using two nutritional models of subacute ruminal acidosis. *Appl Environ Microbiol* 75, 7115–7124 (2009).

34. Caporaso, J. G. *et al.* Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. *ISME J* 6, 1621–1624, doi:10.1038/ismej.2012.8 (2012).

35. Derakhshani, H., Tun, H. M. & Khafipour, E. An extended single-index multiplexed 16S rRNA sequencing for microbial community analysis on MiSeq illumina platforms. *J Basic Microbiol* 56, 321–326, doi:10.1002/jobm.201500420 (2016).

36. Masella, A. P., Bartram, A. K., Truszkowski, J. M., Brown, D. G. & Neufeld, J. D. PANDAseq: paired-end assembler for illumina sequences. *BMC Bioinformatics* 13, 1–7, doi:10.1186/1471-2105-13-31 (2012).

37. Caporaso, J. G. *et al.* QIIME allows analysis of high-throughput community sequencing data. *Nat methods* 7, 335–336, doi:10.1038/nnmeth.f303 (2010).

38. Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C. & Knight, R. UCHIME improves sensitivity and speed of chimera detection. *Bioinformatics* 27, 2194–2200, doi:10.1093/bioinformatics/btr381 (2011).

39. McDonald, D. *et al.* An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. *ISME J* 6, 610–618, http://www.nature.com/ismej/journal/v6/n3/suppinfo/ismej2011139s1.html (2012).

40. Caporaso, J. G. *et al.* PyNAST: a flexible tool for aligning sequences to a template alignment. *Bioinformatics* 26, 266–267, doi:10.1093/bioinformatics/btp636 (2010).

41. Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree: Computing Large Minimum Evolution Trees with Profiles instead of a Distance Matrix. *Mol Biol Evol* 26, 1641–1650, doi:10.1093/molbev/msp077 (2009).

42. Lozupone, C. A., Lladser, M. E., Knights, D., Stombaugh, J. & Knight, R. UniFrac: an effective distance metric for microbial community comparison. *ISME J* 5, doi:10.1038/ismej.2010.133 (2011).

43. Wang, Q., Garrity, G. M., Tiedje, J. M. & Cole, J. R. Naïve Bayesian Classifier for Rapid Assignment of rRNA Sequences into the New Bacterial Taxonomy. *Appl Environ Microbiol* 73, 5261–5267, doi:10.1128/aem.00062-07 (2007).

44. Langille, M. G. I. *et al.* Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. *Nat Biotech* 31, 814–821, doi:10.1038/nbt.2676 (2013).

45. Kanehisa, M. & Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. *Nucleic Acids Res* 28, 27–30, doi:10.1093/nar/28.1.27 (2000).

46. McMurdie, P. J. & Holmes, S. phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. *PLoS One* 8, e61217, doi:10.1371/journal.pone.0061217 (2013).

47. Roesch, L. F. W. *et al.* Pyrosequencing enumerates and contrasts soil microbial diversity. *ISME J* 1, 283–290, doi:10.1038/ismej/journal/v1/n4/suppinfo/ismej200753s1.html (2007).

48. Aguirre de Cárcer, D., Denman, S. E., McSweeney, C. & Morrison, M. Evaluation of subsampling-based normalization strategies for tagged high-throughput sequencing datasets from gut microbiomes. *Appl Environ Microbiol*, doi:10.1128/aem.05491-11 (2011).

49. Lozupone, C. A. & Knight, R. Unifrac: a new phylogenetic method for comparing microbial communities. *Appl Environ Microbiol* 71, doi:10.1128/aem.71.12.8228-8235.2005 (2005).
50. Warwick, R. & Clarke, K. PRIMER 6 (PRIMER-E Ltd, 2006).

51. Anderson, M. PERMANOVA: a FORTRAN Computer Program for Permutational Multivariate Analysis of Variance (Department of Statistics, University of Auckland, 2005).

52. Kuczynski, J. et al. Direct sequencing of the human microbiome readily reveals community differences. *Genome Biol* 11, 1-9, doi:10.1186/gb-2010-11-5-210 (2010).

53. Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biol* 15, 1-21, doi:10.1186/s13059-014-0550-8 (2014).

54. Hammer, O., Harper, D. & Ryan, P. Past: paleontological statistics software package for education and data analysis. *Palaeontol Electron* 4, 9 (2012).

## Acknowledgements

The authors gratefully acknowledge the technical support and guidance of Dr. John Metcalf, Dr. John Doelman, Michelle Carson and Carolyn Innes of Nutreco Canada Inc. (Guelph, ON, Canada). The present research was supported by grants from Natural Science and Engineering Council of Canada and Internal Funding from University of Manitoba.

## Author Contributions

M.S., T.D., and E.K. made substantial contributions to the conception or design of the work. H.D. performed sequencing. S.M. and S.L. conducted data analysis and all authors were involved in the interpretation of data. S.M. and S.L. drafted the manuscript and prepared all figures and tables, and P.A., H.D., J.P., E.K., T.D., and M.S. revised it critically for important intellectual content. All authors approved the submitted versions and agree to be accountable for all aspects of the work ensuring that questions related to the accuracy or integrity of any part of the work were appropriately investigated and resolved.

## Additional Information

Supplementary information accompanies this paper at doi:10.1038/s41598-017-00223-7

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2017