Dietary protein to starch metabolizable energy ratios alter growth performance and gastrointestinal microbiota of calves

The diet structure is very important for the growth and development of calves. This study aimed to investigate the effects of dietary protein-to-starch metabolizable energy ratios (DPSRs) on growth performance, blood index, and gastrointestinal microbiota of calves. Forty-eight Holstein bull calves were fed six dietary DPSRs including A20-35 (20% CP and 35% starch), B20-30, C20-25, D22-35, E22-30, and F22-25 at d 4 to d 60, and then changed to another six dietary DPSRs at d 61 to d 180 (A18-30, B18-27, C18-24, D20-30, E20-27, and F20-24). Twelve calves (d 60) from groups A20-35, C20-25, D22-35, and F22-25 (n = 3) and another twelve calves (d 180) from groups A18-30, C18-24, D20-30, and F20-24 (n = 3) were euthanized. The growth performance parameters were measured. Blood, ruminal fluid, and cecum digesta were collected for further analysis. Results showed heart girth gain of B18-27 was significantly higher than A18-30, C18-24, and heart girth gain (d 180) was significantly affected by protein × starch (DPSRs; p < 0.05). Blood urea nitrogen (BUN; d 60) in C20-25 was significantly higher than A20-35 and B20-30 (p < 0.05). The BUN (d 180) in D20-30 was significantly higher than A18-30 (p < 0.05). The BUN was significantly affected by protein × starch (p < 0.05) on d 60. The albumin (ALB) levels in C20-25 and C18-24 were significantly higher than that in A20-35 on d 60 and A18-30 on d 180, respectively (p < 0.05). The ALB level in D22-35 on d 60 and E20-27 on d 180 was significantly higher than that in other groups (p < 0.05). The ALB level was significantly affected by protein and starch, respectively, on d 60 (p < 0.05). In the rumen, the genera Roseburia (C20-25) and Dialister (D22-35), Prevotellaceae UCG-001 (C18-24), Erysipelotrichaceae UCG-002, and Anaerovorax (F20-24) were found in significant higher relative abundances than those in other groups (p < 0.05). In the cecum, the genera Bacteroides and Eisenbergiella (F22-25), Ruminiclostridium_1 and Candidatus Stoquefichus (A18-30), Erysipelotrichaceae UCG-004 and Tyzzerella 4 (D20-30), and Prevotellaceae UCG-003 and Klebsiella (F20-24) were found in significant higher abundances than those in other groups (p < 0.05). Collectively, these results indicated that the heart girth, BUN, ALB, and gastrointestinal microbiota responded distinctly to differing DPSRs.


Introduction
Protein is an essential nutrient for the growth and development of calves. Hence, it is important to accurately determine the protein requirements of dairy calves to be able to provide the appropriate amino acids for rapid structural growth and lean tissue deposition (Blome et al., 2003). However, the optimum protein content of a starter diet for dairy calves remains controversial. The NRC (2001) recommended having 18% crude protein in calf starters for optimum growth. Many studies have been conducted to determine the optimum protein level in calf starter for proper growth and performance. Akayezu et al. (1994) reported that maximum growth was supported by 19.6% protein in the starter and that calf growth did not increase with higher (22.4%) protein content. Hill et al. (2008) suggested that calves should be provided 18% crude protein (CP) in starter up to 8 weeks of age and 15%-16% CP from 8 to 12 weeks. Moreover, starch is considered as the most important energy source required for the growth of calves (Dennis et al., 2018;Aragona et al., 2020). Hence, the starch was used as the main energy source to provide metabolizable energy in the current study. Rauba et al. (2018) indicated that calves that had greater average daily gain (ADG, >0.80 kg/d) consumed more combined protein and metabolizable energy (ME) from milk replacer and starter than calves with lower ADG. Brown et al. (2005) reported that increased energy and protein intake can increase the rate of body growth of heifer calves and potentially reduce rearing costs. Gabler and Heinrichs (2003) reported varying the CP:ME range from 48.3 to 76.5 did not affect the ADG prepubertal Holstein heifers. According to Lammers and Heinrichs (2000), with the ratio of dietary protein to energy increasing, the average daily gain, and feed e fficiency of heifers were increased linearly. So understanding the influence of protein: energy (starch metabolizable energy) on a calf is valuable.
The ruminal microbiota is important for nutrient absorption and energy metabolism of the host (Belanche et al., 2012) and different microbes play different roles in the rumen. For example, bacteria and protozoa break down starch and fiber by releasing cellulolytic enzymes, which are then anaerobically fermented to metabolic intermediates such as volatile fatty acids, ammonia, carbon dioxide, and hydrogen (Kamra, 2005;Kebreab et al., 2009). Ruminants are born without a functional rumen, however, microbial communities exist in a newborn's rumen, and rumen bacterial communities are thought to be influenced by the diet and age of the animal (Jami et al., 2013;Jewell et al., 2015). Malmuthuge et al. (2013) found that different calf starters affected the bacterial community among gut regions and sample types. In addition, Jami et al. (2013) observed that rumen bacterial diversity and species richness increased with increasing age. Malmuthuge et al. (2014) investigated the composition of the bacteria in the gastrointestinal tracts of pre-weaned bull calves and indicated that the bacterial community was the most abundant in the rumen, followed by the large intestine and then the small intestine.
The relationship between dietary protein and energy to calf growth is well-documented above; however, very limited information exists on how the rumen and the cecum microbiomes are affected by different diet structures during pre-weaning and post-weaning periods. We hypothesized that diet with appropriate crude protein-to-starch metabolizable energy ratios (DPSRs) could improve the growth performance and gastrointestinal microbiota of calves. Therefore, the purpose of this study was to investigate the effects of starters with different DPSRs on the growth performance and rumen and cecum microbiomes of calves.

Materials and methods
The animals were managed according to the guidelines for the care and use of experimental animals of Jilin Agricultural University Care and Use Committee and commercial dairy farms (JLAU-ACUC2019-018, Changchun, China).

Animal housing and feeding
The calves during the pre-weaning period were housed in individual calf hutches (2.20 m × 1.32 m × 1.39 m) with pens (1.75 m × 1.30 m × 1.10 m) until d 60. The calf hutches were bedded with sand and pens were equipped with a trough (diameter: 25 cm, depth: 16 cm) for the calf starter and a bucket for milk and water. The calves were given colostrum (IgG > 50 mg/ mL) at 5 L/d for 4 days and the calf starter was offered from d 4. The calves were fed up to 6 L/d of commercial milk replacer (Sprayfo Red, Sloten BV, Antwerpenweg 7, 7418 CR Deventer, Netherlands) from d 5 to d 7. The amount of milk replacer fed to calves was gradually increased to 7 L/d and 8 L/d from d 8 to d 14 and d 15 to d 21, respectively. From d 22 to d 37 and d 38 to d 51, calves were fed milk replacer at 9 and 10 L/d, respectively. On day 52, weaning was initiated by gradually reducing the milk replacer to 1 L/d. So, the calves were completely weaned off milk replacer at d 60. Each calf drank milk twice a day during the pre-weaning period. At d 61, the calves were not fed milk replacer and moved to 6 pens (6.2 m × 5.0 m) according to their treatments. The d 61 to d 66 constituted the transition period to using the two-stage diets according to the management system; the calf starter was gradually decreased and the concentrate (post-weaning period) was gradually increased in different ratios (proportion of starter to concentrate: 6:1, 5:2, 4:3, 3:4, 2:5, and 1:6 on d 61, d 62, d 63, d 64, d 65, and d 66, respectively). Each post-weaning calf began to feed on oat grass (0.2 kg/d) at d 90 and the calves were fed twice daily. Water, calf starter, and concentrate were always freely available.

DNA extraction and sequencing
Total genomic DNA samples were extracted using Fast DNA SPIN extraction kits (MP Biotechnology, CA, United States) according to the manufacturer's specifications. The DNA concentration was measured using a NanoDrop ND-1000 (Thermo Scientific, Massachusetts, United States). The V3-V4 regions of the 16S rRNA genes were amplified using the universal primers 338F (5′-ACTCCTACGGGAGGCAGCA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). The PCR was conducted according to a previous method (Dill-McFarland et al., 2017). The PCR products were purified and quantified using VAHTSTM DNA Clean Beads (Vazyme, Nanjing, China) and PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, United States), respectively. The purified PCR products were pooled and sequenced using the Illumina MiSeq platform with a MiSeq Reagent Kit v3 at Shanghai Personal Biotechnology Co., Ltd. (Shanghai, China). The sequencing data were processed using the Quantitative Insights into Microbial Ecology (QIIME 2) pipeline (Bolyen et al., 2018). Briefly, the primer sequence was trimmed off and the remaining sequence was then quality filtered, denoised, and merged, and chimeras were removed using the DADA2 plugin (Callahan et al., 2016). The remaining high-quality sequences were re-clustered at 97% to generate amplicon sequence variants (ASVs) using the clustering program Vsearch (v2.13.4;Rognes et al., 2016). The sample sequence was flattened based on ASV, and the average score of the maximum flattening depth was selected to calculate the Alpha diversity indices such as Chao1, Shannon, and Pielou_e. The amplicon sequence data were deposited with the National Center for Biotechnology Information (Accession No. SRP348482).

Statistical analyses
Data were analyzed by the SPSS software version 23.0, the data were analyzed as a linear mixed-effects model, which included crude protein, Frontiers in Microbiology 04 frontiersin.org starch (starch metabolizable energy), and their interaction (CP × starch; DPSRs) as fixed effects. Data were analyzed using the model: where μ is the overall mean, F i is the effect of crude protein (i = 1-2), V j is the effect of starch (j = 1-3), F i × V j is the dietary crude protein × starch level interaction, and e ijk is the residual effect. Then, the treatment groups of pre-weaning and post-weaning period were analyzed by one-way ANOVA, respectively. The results were expressed as the mean ± SE (standard error), and the differences among means were tested for significance by Duncan's multiple range test. p < 0.05 was considered to indicate a statistically significant difference and tendency toward significance was declared at 0.05 < p ≤ 0.10.

Growth performance
No significant difference was observed in overall growth performance during the experimental period (p > 0.05, Table 3). For lower protein level (pre-weaning period: 20%, post-weaning: 18%) treatments, the average body weight on days 60 and 180, and ADG from d 4 to d 60 of group C20-25 were higher than other groups, whereas, ADG of group B18-27 was highest from d 61 to d 180 (p > 0.05). For higher protein level (pre-weaning period: 22%, post-weaning period: 20%) treatments, average body weight on d 60 and ADG of group D22-35 were higher than other groups before weaning; however, average body weight on d 180 and ADG of group E20-27 were higher than other groups after weaning (p > 0.05). We observed a significant protein × starch interaction (DPSRs) for heart girth (p = 0.002) on d 180. Heart girth gain of group B18-27 was significantly higher than groups A18-30 and C18-24 on d 180 (p < 0.05). When the starch level in the diet was higher, the body length gain and cannon bone circumference gain of group A18-30 were significantly higher than for group D20-30 (p < 0.05). Cannon bone circumference gain was significantly affected independently by protein (p = 0.02) and starch (p = 0.01) on d 180. Therefore, according to these data, the calves fed the 20% CP in the starter with a ratio of 155.81 g (group A20-35) and 218.28 g (group C20-25) of CP/Mcal of starch and 22% CP in the starter with a ratio of 150.54 g (group D22-35) and 240.22 g (group F22-25) of CP/Mcal of starch were euthanized on d 61 for further research.

Blood index
The total protein (TP) levels in serum were not significantly different among the groups before weaning when calves were fed different diets (Table 4). When the protein amount in the diet was lower (18% or 20%), the ALB level increased with DPSRs increasing during the experimental period, and the ALB level in group C20-25 and C18-24 was significantly higher than that in group A20-35 on d 60 and group A18-30 on d 180, respectively (p < 0.05). ALB level was significantly affected independently by protein and starch, respectively, on d 60 (p < 0.05). For groups with a higher CP level, the ALB level in the group D22-35 on d 60 and group E20-27 on d 180 was significantly higher than that in other groups (p < 0.05). For the groups with the same starch level, the urea amount was significantly affected by starch and protein × starch on d 60 (p < 0.05)

Composition of rumen and cecum microbiomes
At the phylum level, the rumens of different groups had similar dominant microbial communities, which included Firmicutes, Bacteroidetes, and Actinobacteria ( Figure 1A). However, the abundance of the same phyla in the pre-weaning and post-weaning periods showed differences. For instance, Firmicutes was the predominant phylum, followed by Bacteroidetes, in the rumen during the pre-weaning period, but the opposite result was detected in groups during the post-weaning period, during which Bacteroidetes was the predominant phylum. The abundance of Bacteroidetes in the rumen of groups C20-25 and D20-30 was higher than that in other groups in pre-weaning and post-weaning periods. However, Bacteroidetes was more abundant in the cecum of groups F22-25 (F20-24) than that in other groups in pre-weaning and post-weaning periods ( Figures 1A,B). Visible differences in the microbial relative abundance were detected between the rumen and the cecum. Firmicutes, Bacteroidetes, Tenericutes, and Verrucomicrobia were the four most dominant phyla in the cecum ( Figure 1B). In both the pre-weaning and post-weaning periods, the abundance of Firmicutes in the cecum was higher than that of Bacteroidetes, and the abundance of Firmicutes groups C20-25 and G20-30 was more abundant than that in other groups in pre-weaning and post-weaning periods.
At the genus level, 284 genera in the rumen and 381 genera in the cecum were found in the weaning period. During the post-weaning period, 388 and 454 genera were detected in the rumen and cecum, respectively. Additionally, most genera in both the rumen and cecum belonged to Firmicutes and Bacteroidetes. Figures 2A,B show the genera whose relative abundance was more than 1%. In the rumen, Prevotella 7 was the most abundant genus during the pre-weaning period (average relative abundance more than 35% in different groups), whereas Prevotella 1 in groups C18-24 (average relative abundance was 33%) and D20-30 (average relative abundance was 46%), the genus Bifidobacterium in group A18-30 (average relative abundance was 22%), and the genus The postfix of the groups mean the levels of protein and starch in diet. For example, A20-35 means 20% CP and 35% starch in diet was supplied in group A.
Frontiers in Microbiology 07 frontiersin.org Prevotella 7 in group F20-24 (average relative abundance was 20%) were prominent during the post-weaning period, respectively. The abundance of Succiniclasticum was higher during the pre-weaning period than in the post-weaning period. In the cecum, the relative abundance of Ruminococcaceae UCG-005 was the highest (average relative abundance was from 18 to 38%) during the pre-weaning and post-weaning periods. Genus Ruminococcaceae UCG-005 was more abundant in high DPSRs groups (groups C20-25 and F22-25) than that in low DPSRs groups (groups A20-35 and D22-35) on d 60; however, opposite results were observed on d 180, which was Ruminococcaceae UCG-005 was more abundant in low DPSRs groups (groups A18-30 and D20-30) than that in high DPSRs groups (groups C18-24 and F20-24) on d 180. Bacteroides and [Eubacterium]coprostanoligenes were the second most abundant genera in the pre-weaning and post-weaning periods, respectively. The result showed that the ASV flower petal map of rumen microbes was increased during the post-weaning period when compared with the pre-weaning period; however, the opposite results were found in the cecum (Figures 3A,B). The number of ASV in the rumen of group A18-30 was higher than that of other groups during the post-weaning period. The number of ASV in the cecum of group A20-35 was the highest and that of group C20-25 was the lowest on d 60; however, the opposite results of ASV were shown on d 180, which group C18-24 was the highest and group A18-30 was the lowest. No significant difference in Alpha diversity indexes about rumen and cecum was obtained in the Frontiers in Microbiology 08 frontiersin.org present research, however, the variation trend of microbes could be observed (Figures 4A,B). For the lower protein groups (pre-weaning period: 20% CP; post-weaning period: 18% CP), Chao1, Shannon, and Pielou_e in the rumen of group A20-35 were higher than those of other groups (p > 0.05) during the pre-weaning period, and Chao1, Shannon, and Pielou_e of group D20-30 were higher than those of other groups (p > 0.05) during the post-weaning period. Chao1 in the cecum of group A20-35 and Shannon and Pielou of group C20-25were higher than those of other groups (p > 0.05) during the pre-weaning period. Chao1 in the cecum of group C18-24, and Shannon and Pielou of group A18-30 were higher than those of other groups (p > 0.05) during the post-weaning period. For the higher protein groups (pre-weaning period: 22% CP; post-weaning period: 20% CP), Chao1, Shannon, and Pielou_e in the rumen of group C20-25 were higher than those of other groups (p > 0.05) during the pre-weaning period, and Chao1, Shannon, and Pielou_e of group F20-24 were higher than those of other groups (p > 0.05) during the post-weaning period. Chao1, Shannon, and Pielou_e in the cecum of group F22-25 were higher than those of other groups (p > 0.05) during the pre-weaning period, and Chao1of group F20-24 and Shannon and Pielou_e of group D20-30 were higher than those of other groups (p > 0.05).
To identify specific microbial candidates present in the different groups during the pre-weaning and post-weaning periods, LEFSE (LDA Score > 2) and p < 0.05 were used to determine biomarkers belonging to different taxa. In the rumen, the families Lactobacillaceae, Streptococcaceae, and Enterobacteriaceae were higher in the group A20-35 during the pre-weaning period. The genera Roseburia and Prevotellaceae UCG-001 were higher in groups C20-25 and C18-24 during the pre-weaning and post-weaning periods, respectively. The genera Dialister and U29-B03 were higher in groups D22-35 and D20-30 during the pre-weaning and post-weaning periods, respectively. The order Bacillales was the biomarker in group F22-25 during the pre-weaning period, and the genera Erysipelotrichaceae UCG-002 and Anaerovorax in group F20-24 were found in higher abundances during the post-weaning period ( Figures 5A,B). In the cecum, the family Bacteroidaceae and the genus Bacteroides and Eisenbergiella were the biomarkers in group F22-25 during the pre-weaning period. The genera Ruminiclostridium_1 and Candidatus Stoquefichus in group A18-30; the family Eubacteriaceae in group C18-24; the phylum Tenericutes and the genera Erysipelotrichaceae UCG-004 and Tyzzerella 4 in group D20-30; and the family Clostridiales vadin BB60 and the genera Prevotellaceae UCG-003 and Klebsiella in group F20-24 during the post-weaning were found in higher abundances ( Figures 5C,D).

Network analysis of microbial metabolic pathways
Microbiomes of different groups were imported into the MetaCyc database 1 for pathway analysis, and the significant metabolic pathways (p < 0.05) with microbes are shown in Figure 6 when compared among different levels of DPSRs. A total of 6, 18, 10, and 6 pathways were related to microbiomes in the rumen (d 60, d 180) and the cecum (d 60, d 180), respectively. In the lower protein groups (pre-weaning period: 20% CP; post-weaning period, 18% CP; Figure 6A, A vs. C, Figure 6B, A vs. D), microbiomes in the rumen (d 60) and cecum (d 60, d 180) were involved in the common pathways of aromatic compound degradation, amine and polyamine degradation, and amino acid degradation, and the rumen (d 180) and cecum (d 60, d 180) were involved in the common pathway of carbohydrate degradation. For lower CP groups, the pathways of the cofactor, prosthetic group, electron carrier, vitamin biosynthesis(group C20-25), amine and polyamine degradation (group C20-25), and amino acid degradation (group C18-24) were promoted in the cecum in high DPSRs groups compared with low DPSRs group; however, other related pathways in the rumen and cecum were restrained. For the higher CP groups (pre-weaning period: 22% CP; post-weaning period: 20% CP; Figure 6D, D vs. F), the pathways of aromatic compound degradation (group F22-25), ethylmalonyl-CoA pathway (group F22-25), and nucleoside and nucleotide biosynthesis (group F22-25) were promoted in the cecum in higher DPSRs group when compared with the lower DPSRs group, whereas other related pathways in the rumen and the cecum were restrained. For the 35% starch level groups ( Figure 6B, A vs. D), only the pathway of carbohydrate degradation, related to the genus Bacteroides in the rumen on d 180, was promoted with increasing CP. For the 25% starch level groups ( Figure 6C, C vs. F), the pathways of amine and polyamine 1 https://metacyc.org/
Frontiers in Microbiology 09 frontiersin.org degradation (groups F22-25 and F20-24) in the cecum, and aromatic compound degradation (groups F20-24) in cecum and rumen were promoted in the 22% CP groups when compared with 20% CP groups.

Growth performance
The proper nutrition level and diet structure are critical for the growth and development of calves, especially the level of dietary protein and energy, and the ratio of protein to energy. Zothanpuii et al. (2015) found that crossbred calves fed lower protein calf starter with higher nitrogen-free extract (22% CP) showed a significantly increased body weight compared to those fed calf starters with higher protein (24% and 26% CP) and lower nitrogen-free extract. Our results showed the similar growth tendency; however, no significant difference among different groups, which we hypothesized the BW appears to be affected by many factors, including the individual difference and interaction of the level of protein and starch in the diet. Stamey et al. (2012) found that a high-CP starter (25.5% CP; ME: 3.21 Mcal/kg) could significantly increase the ADG of calves compared to those fed conventional starter (19.6% CP, ME: 3.16 Mcal/kg) before weaning. Yousefinejad et al. (2021) reported that female Holstein dairy calves supplied with 22% CP in starter showed greater overall ADG, pre-weaning BW, and feed efficiency than that supplied with 18% CP. In the present study, all the calves were fed the same amount of milk before weaning, and it should be noted that, when the diet with lower CP (20% CP in the pre-weaning period and 18% CP in the post-weaning periods) was provided to the calves during the experimental period, their growth was not constrained during the developmental stage, i.e., when their requirements for protein were high. This is probably because the digestible energy intake that was provided by starch was sufficient. For the groups with the same level of protein, lower CP (20%) with lower starch (25%) of 218.48 g of CP/Mcal of starch ME and higher CP (22%) with higher starch (35%; 170.54 g of CP/Mcal of Starch ME) showed the higher BW during the pre-weaning period. Similar results were found by Lammers and Heinrichs (2000), which found that calf starters with 35.3%, 33.4%, and 31.4% starch and 27.1%, 25%, and 24.7% CP, respectively, did not significantly affect ADG. However, in the current study, higher CP (20%) with higher starch (35%) did not affect the body weight of calves when compared with the lower levels of starch (30 and 25%) during the post-weaning period. Accordingly, different levels of CP could influence calf performance, but a proper combination of dietary nutrient levels and CP: ME ratio balance for calves at different stages of development may reduce or eliminate this effect and may even benefit calf performance. Lower CP (20%) with lower starch (25%; 218.48 g of CP/ Mcal of Starch ME) or higher CP with higher starch (170.54 of CP/Mcal of Starch ME) before weaning, and lower CP with medium-level starch (180.05 of CP/Mcal of Starch ME) or higher CP with medium-level starch (205.05 of CP/Mcal of Starch ME) showed greater benefits for calf performance.

Blood index
The serum TP did not change considerably with different DPSRs, but the ALB concentrations in the groups C20-25 and C18-24 (higher DPSRs) was higher than that in groups A20-35 (A18-30) and B20-30 (B18-27), and the group D22-35 (lower DPSRs) were higher than that in groups E22-30 and F22-25 in the current study. In accordance with the present results, Daneshvar et al. (2016) and Makizadeh et al. (2020) reported no significant change in TP after feeding calves with different levels of CP; furthermore, both results showed no significant effect on ALB concentration. In contrast to our results, Sharma et al. (2019) found that TP concentration increased with higher levels of protein in the starter, but no significant change was observed in ALB. Blood urea N (BUN) concentration is generally considered an indicator of protein status (Kohn et al., 2005). In the present study, the BUN concentration was higher in the 20% CP groups than in the 18% CP groups during the post-weaning period, which may be due to excess ruminal nitrogen concentrations and lower protein requirements for the growth than in the pre-weaning period (Kazemi-Bonchenari et al., 2018). For the 20% CP groups, the concentration of BUN was decreasing with the starch decreasing; however, seldom research was about the effect of starch on the blood index. We hypothesized that the changes in blood might be related to growth performance and the ADG was lowest when the level of starch was 27% in the current study.

Composition of rumen and cecum microbiomes
The current study characterized the ruminal and cecal bacterial composition in calves with different DPSRs in the feeding trial; we also identified variations in bacterial diversity and specific ruminal and cecal bacteria in different groups. Jami et al. (2013) found that the ruminal bacterial diversity of calves increased from 2 to 6 months of age, which is in agreement with the current results. The ruminal bacterial community during the post-weaning period was more abundant than in the pre-weaning period; however, the opposite trend was found in the cecal bacterial community, probably because the calf rumen is not fully matured by d 60 its function was similar as monogastric animals. Thus, the cecal bacterial diversity was greater than that of the rumen, which was consistent with the results of Malmuthuge et al. (2014). During the post-weaning (d 180), calves developed a more mature rumen, suggesting a more diverse but specific and homogeneous bacterial community (Jami et al., 2013). In the present study, the same dominant phyla were present Frontiers in Microbiology 11 frontiersin.org in the rumen and cecum, although the relative abundances varied depending on the calf starter and the stage of development. Moreover, the relative abundance of Bacteroidetes increased, and that of Firmicutes decreased in the rumen from pre-weaning to post-weaning, which is inconsistent with the findings of Jami et al. (2013) and Meale et al. (2016). Conversely, the relative abundance of Bacteroidetes decreased and that of Firmicutes increased in the cecum from d 60 to d 180. Prevotella belongs to the phylum Bacteroidetes, and Prevotella can produce a range of xylanases (Krause et al., 2003) and may contribute to the efficient utilization of hemicellulose pectin and protein (Osborne and Dehority, 1989;Satoshi and Yasuo, 2001). The increase in Bacteroidetes in the rumen is indicative of the gradual rumen development to maturity and is promoted by processes utilizing oat grass and calf starter during the post-weaning period.
Previous studies have reported that diet is one of the main factors that can change gut microbial diversity (Ley et al., 2008;Maslowski and Mackay, 2011;Dill-Mcfarland et al., 2019). The abundance of each phylum or genus was influenced when calves were supplied with different DPSRs in the present study. In agreement with a previous study in dairy calves by Meale et al. (2016), the genus Succiniclasticum was more abundant in pre-weaning than that in post-weaning. Succiniclasticum specializes in fermenting succinate, quantitatively converting it to propionate (Gylswyk, 1995); furthermore, a higher proportion of propionate will be produced in the rumen compared with the period of supplementation with oat grass after weaning (Kim et al., 2016). In this case, the currently observed decrease in abundance of Succiniclasticum may have been because no forage was supplemented for calves before weaning. Petri et al. (2013) reported that Succiniclasticum in heifers fed a high-grain diet was more abundant than in those fed forages or a mixed forage diet. Moreover, the current study showed that the phylum Bacteroidetes was dominated by the genus Prevotella in the rumen at all stages of development; however, in the cecum, the genus Bacteroides was the most abundant, although the abundance decreased after weaning, which was also found by Meale et al. (2016). Furthermore, for groups with the same level of CP, most of the pathways in which microbes in the rumen and cecum are involved became restrained when the starch concentration in the diet was increased; similar changes were obtained in the groups with the same starch level, and most of the pathways were restrained in higher CP groups when compared with lower CP groups.

Conclusion
This study found that heart girth gain (d 180) and BUN (d 60) was significantly affected by protein × starch (DPSRs). The level of ALB was significantly affected independently by protein and starch, respectively, on d 60 during the whole experimental period, the gastrointestinal microbiota were different among the dietary treatments. Compared with pre-weaning period, the rumen microbial abundance increased in postweaning period, whereas the cecum microbial abundance was decreased. In a word, these results indicated that the BUN, ALB, and gastrointestinal microbiota responded distinctly to differing DPSRs or period, and a reasonable DPSRs for calves in different stages of development is important.

Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, SRP348482.

Ethics statement
The animals were managed according to the guidelines for the care and use of experimental animals of Jilin Agricultural University Care and Use Committee and commercial dairy farms (JLAU-ACUC2019-018, Changchun, China).

Author contributions
YZ, ZS, TW, GQ, and XZ: conceptualization. YZ and XC: methodology. XC, ZH, JD, JX, and JR: formal analysis. XC, YZ, JX, WZ, NA, and ZS: data curation. GQ, YZ, TW, and ZS: supervision. XC and ZH: writing-original draft. XC, YZ, ZS, and TW: writing-review and editing. All authors contributed to the article and approved the submitted version.

Funding
This study was supported by the National Natural Science Foundation of China (grant number 32102574) and Scientific and Technological Developing Scheme of Jilin Province (grant numbers 20210202037NC and 20220202049NC).

Conflict of interest
YZ and ZS were employed by Changchun Borui Science and Technology Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.