Fecal microbiota dynamics and its relationship to diarrhea and health in dairy calves

Hongwei Chen, Yalu Liu, Kailang Huang, Bin Yang, Yuanyuan Zhang, Zhongtang Yu and Jiakun Wang*

Abstract Background: Diarrhea is a major cause of morbidity and mortality in young calves, resulting in considerable economic loss for dairy farms. To determine if some gut microbes might have resistance to dysbiotic process with calf diarrhea by dictating the microbial co-occurrence patterns from birth to post-weaning, we examined the dynamic development of the gut microbiota and diarrhea status using two animal trials, with the first trial having 14 Holstein dairy calves whose fecal samples were collected 18 times over 78 d from birth to 15 d post-weaning and the second trial having 43 Holstein dairy calves whose fecal samples were collected daily from 8 to 18 days of age corresponding to the first diarrhea peak of trial 1.Results: Metataxonomic analysis of the fecal microbiota showed that the development of gut microbiota had three age periods with birth and weaning as the separatrices. Two diarrhea peaks were observed during the transition of the three age periods. Fusobacteriaceae was identified as a diarrhea-associated taxon both in the early stage and during weaning, and Clostridium_sensu_stricto_1 was another increased genus among diarrheic calves in the early stage. In the neonatal calves, Prevotella_2 (ASV4 and ASV26), Prevotella_9 (ASV43), and Alloprevotella (ASV14) were negatively associated with Clostridium_sensu_stricto_1 (ASV48), the keystone taxa of the diarrhea-phase module. During weaning, unclassified Muribaculaceae (ASV28 and ASV44), UBA1819 (ASV151), Barnesiella (ASV497), and Ruminococcaceae_UCG-005 (ASV254) were identified being associated with non-diarrheic status, and they aggregated in the non-diarrhea module of co-occurrence patterns wherein unclassified Muribaculaceae (ASV28) and Barnesiella (ASV497) had a direct negative relationship with the members of the diarrhea module.Conclusions: Taken together, our results suggest that the dynamic successions of calf gut microbiota and the interactions among some bacteria could influence calf diarrhea, and some species of Prevotella might be the core microbiota in both neonatal and weaning calves, while species of Muribaculaceae might be the core microbiota in weaning calves for preventing calf diarrhea. Some ASVs affiliated with Prevotella_2 (ASV4 and ASV26), Prevotella_9 (ASV43), Alloprevotella (AVS14), unclassified Muribaculaceae (ASV28 and ASV44), UBA1819 (ASV151), Ruminococcaceae_UCG-005 (ASV254), and Barnesiella (ASV497) might be proper probiotics for preventing calf diarrhea whereas Clostridium_sensu_stricto_1 (ASV48) might be the biomarker for diarrhea risk in specific commercial farms.

Keywords: Calf diarrhea, Co-occurrence pattern, Dynamic development, Fecal microbiota

Background

Diarrhea is one of the most common causes of morbidity and mortality in young calves, especially dairy calves younger than one month [1, 2]. In the US, the National Animal Health Monitoring System estimated that more than half of the calf deaths were due to diarrhea and related diseases [3]. A similar mortality rate among dairy calves attributable to diarrhea was also reported in Korea [1]. Calf diarrhea not only incurs treatment cost but also lowers daily weight gain in calves [4], delays first conceptions [5], decreases milk production in the first lactation [6], consequently resulting in a great economic loss for dairy farms. Many pathogens have been implicated in calf diarrhea, includingEscherichia coliK99+(E. coliK99+),Salmonellaspp.,Clostridium perfringens,Cryptosporidium parvum, bovine rotavirus (BRV), and bovine coronavirus (BCoV) [1]. Non-infectious factors such as improper colostrum management [7], weaning stress [8], and poor feeding environments can also increase the risk of calf diarrhea [9]. Moreover, oftentimes uncertain causality of diarrhea occurrence and cross infection together make diarrhea difficult to prevent and treat.

Diarrhea is much more common among young animals, especially after birth to shortly after weaning than among adult animals because of higher risk of enteric pathogens colonization in young animals than in adult animals [10, 11]. The high risk of colonization by pathogens in young animals stems from their underdeveloped gut and gut microbiota. Indeed, Kim et al. found that the absence of some species of Clostridia in the gut microbiota of neonatal mice made them unable to resist the colonization bySalmonella entericaTyphimurium orCitrobacter rodentium[12]. In addition, dysbiosis of gut microbiota induced by certain factors such as antibiotics and diseases can also lower resistance and increase enteric pathogen infection. Wu et al. demonstrated that dysbiosis of gut microbiota in rat induced by antibiotics lowered the resistance to colonization bySalmonella[13]. Compared to uninfected calves, calves infected with rotavirus have a higher abundance ofEscherichia,Clostridium_g21,Streptococcus, andClostridium, but a lower abundance ofLactobacillus,Subdoligranulum,Blautia, andCoprococcus_g2 in the fecal microbiota [14].Weaning is a stressful process for young animals and can lower resistance to pathogenic colonization. Haag et al. found that infant mice after weaning were deficient in preventing their gut from colonization byCampylobacter jejuni[15]. During weaning, dairy calves are stressed nutritionally, and their gut microbiota undergoes compositional changes [16]. Typically, during weaning, the gut microbiota increases its diversity, and some important taxa, such asBacteroides,Blautia,Ruminococcus, andSuccinivibrio, may change their abundance [17]. The nutritional stresses coupled with changes in gut microbiota during and immediately after weaning increase the risk of enteric diseases, particularly diarrhea [17, 18].

Interestingly, some calves are prone to diarrhea, whereas others are resistant to diarrhea. In the 56-day experiment of Ma et. al [19], 13 out 42 milk replacer fed calves never exhibited diarrhea, 18 calves exhibited diarrhea and recovered after treated with electrolyte, but 11 calves exhibited diarrhea and needed to be treated with therapeutic antimicrobials. Ma et al. defined the latter two diarrhea status as resistant to diarrhea-induced dysbiosis, and susceptible to diarrhea-induced dysbiosis, respectively. The age of first diarrhea of these calves varied from d 8 to 19. Using the random forest machine learning algorithm with the microbiota collected from health calves at 7, 14, and 21 days of age and diarrhea calves prior to the onset of diarrhea, Ma et al. suggest that diarrhea could be predicted by the microbiota shift in early life. Kim et al. [20] described the similar age (5–50 days of age) fecal microbiota transplantation could ameliorate calf diarrhea with increasing the family Porphyromonadaceae. Therefore, we hypothesized that the gut microbiota might strongly correlate to calf diarrhea for some gut microbes might have resistance to the dysbiotic process with calf diarrhea from birth to post-weaning. To test our hypothesis, the stages with high incidence of diarrhea were identified in cohort from birth to shortly after weaning, and the fecal microbiotas between diarrheic calves and non-diarrheic calves in stages of high incidence of diarrhea, and between diarrhea phases (prediarrhea, diarrhea and post-diarrhea) were compared, and the diarrheic status-and diarrhea phase-associated amplicon sequence variants (ASVs) and modules were identified.

Methods

Animal experiment and sample collection

Two animal trials were conducted on a commercial farm located in Shaoxing (more than 3800 cows in stock), Zhejiang Province, China from September to November 2017 (trial 1) and November to December 2017 (trial 2). In trial 1, 14 female Holstein calves (initial bodyweight = 38.2 ± 2.0 kg, mean ± SD) were enrolled at birth. Immediately after birth, the calves were separated from their dams and moved to individual pens (1.0 m × 1.2 m) bedded with dry wheat straw and each offered 4 L of thawed colostrum (pooled the colostrum from other mothers before the trial, and stored at – 80 oC) via an esophageal tube within 2 h after birthing. All neonatal calves in this trial were fed the same colostrum. From 1 to 5 days of age, each calf was bottle-fed 2 L of whole milk each at 0700, 1330, and 1800 h. During these 5 d, 4 g each of tylosin tartrate and sulfadimidine soluble (both powder) was mixed with the morning milk and fed to each calf for prophylaxis of pneumonia. At d 6, all the calves were moved to one group pen (8.0 m × 10.0 m) bedded with a rubber mat and dry wheat straw. The straw was replaced every other day. The calves had free access to a preset volume of whole milk (Additional file 1: Fig. S1) with an automated milk feeder system (F?rster-Technik, Engen, Germany). Briefly, the calves get 6 L/d of whole milk from d 6 and this amount increased evenly up to 8 L/d on d 22. From d 23 to 41, each calf received 8 L of whole milk daily, and from d 42 to 62, the milk allowance was decreased evenly from 8 to 2 L/d. Milk supply was denied using the automated milk feeder system once a calf consumed 2 L of milk within 2 h during each meal to avoid overconsumption of milk. All the calves were weaned off milk at d 63. Calf starter pellets and a hay mixture of oat and alfalfa were offered ad libitum in the group pen. All the calves had free access to clean drinking water throughout the trial. Fecal samples were collected directly via rectal stimulation from every calf of the trial at d 1, 3, 5, 7, 9, 12, 15, 18, 28, 38, 48, 58, 62, 63, 65, 68, 73, and 78 (18 times in total). The fecal samples were placed in 2 mL nucleasefree tubes, immediately frozen in liquid nitrogen, then each month the 15 L liquid nitrogen tank with samples was transported to the lab and samples were stored at – 80 °C. The DNA of all the stored fecal samples was extracted in one month after the trial. When fecal samples were collected, fecal scores were recorded based on fecal fluidity [21]: 1 = normal, 2 = soft, 3 = runny, or 4 = watery. Calves that had a fecal score of 3 or 4 were considered diarrheic. No additional antibiotics were offered to calves after the diarrhea episode.

In trial 2, 43 female Holstein calves (initial bodyweight = 36.6 ± 1.9 kg, mean ± SD) were enrolled at birth. The calves were fed, nursed, and sampled exactly the same as for those in trial 1, but fecal samples were collected daily from d 8 to 18, which corresponded to the age when the first diarrhea peak was observed in trial 1.

Metataxonomic analysis of the fecal microbiota

DNA was extracted from individual fecal samples according to the method of Zoetendal et al. with minimal modifications [22]. In brief, approximately 0.2 g of frozen fecal sample each was homogenized together with 1 mL of TE buffer using bead-beating (Biospec Products; Bartlesville, OK, United States) followed by phenol:chloroform-based DNA extraction. Agarose gel (1%) electrophoresis was performed to evaluate the DNA quality, and the DNA concentrations were determined using a NanoDrop 2000 spectrophotometer (Thermo Scientific, Waltham, MA, United States). The V3-V4 region of the 16S rRNA gene was amplified using primers 341F (5’-CCT AYG GGRBGCASCAG-3’) and 806R (5’-GGA CTA CNNGGG TAT CTAAT-3’) to prepare amplicon libraries with each being labeled with a unique barcode sequence [23]. The amplicon libraries of all the samples were pooled in equal molar ratio and paired-end (2 × 250) sequenced on the Illumina HiSeq platform by Novogene Bioinformatics Technology Co., Ltd. (Tianjin, China). After demultiplexing, the paired-end sequencing reads were processed using the DADA2 package (version 1.8.0) in R and its pipeline [24]. Briefly, barcodes and primers were removed from the reads. Reads with more than 2 expected errors (maxN = 0, maxEE = c(2, 2), truncQ = 2) were filtered out. Dereplication and inference were performed using the DADA2 pipeline. After merging the paired reads and chimera filtering, an ASV table was constructed (to resolve bacteria at the species level [25]). The ASVs were taxonomically assigned based on the SILVA 16S rRNA gene database (version 132) [26] using a naive Bayesian classifier method [27] implemented in DADA2. The sequences of each sample were rarefied to same size with the minimum number of sequences in samples (35,137 sequences/sample in trial 1, and 19,986 sequences/sample in trial 2) using the ‘rarefy_even_depth’ function in Phyloseq package (version 1.24.2) [28], and ASVs that appeared only once among all the samples were removed from the dataset. Alpha diversity metrics including observed species and Shannon diversity index were calculated using the Phyloseq package. Faith’s phylogenetic diversity (Faith’s PD) was calculated using the Picante package (version 1.8.2) in R, and evenness was calculated using the Microbiome package (version 1.12.0) in R. To minimize individual variance of calves, only the ASVs that were observed in at least 20% of calves at every single day (3 out of 14 calves in trial 1) or every single diarrheic status-associated phase (trial 2) were subjected to the downstream analysis.

Pathogen detection

Major pathogens that can lead to calf diarrhea were tested in the fecal samples of trial 1 when diarrhea (fecal score ≥ 3) was first noted in a calf.Salmonellaspp.was detected by PCR with specific primers (F:5′-TCG TCA TTC CAT TAC CTA CC-3′ and R:5′-AAA CGT TGA AAA ACT GAG GA-3′ [29] using fecal DNA samples. A commercial ELISA kit (BIO K 315, Bio-X Diagnostics, Rochefort, Belgique) each was used to detect BRV, BCoV, andE. coli(through its F5 attachment factor) in the fecal samples.

Evaluation of the effects of age, diarrheic status, and diarrheic phases on fecal microbiota

A generalized linear mixed-effects model implemented in the nlme package (version 3.1.137) [30] in R was used to evaluate the effect of age and diarrhea on the alpha diversity metrics of the fecal microbiota, and Tukey’s allpair comparison test using the ‘glht’ function in the multcomp package (version 1.4.8) [31] in R was used to do the multiple comparisons. In trial 1, the model included calf age and diarrheic status (diarrheic vs. non-diarrheic calves) among the calves as fixed effects and individual calves as random effect:

In trial 2, the model included diarrhea phases (prediarrhea, diarrhea, and post-diarrhea, see the results for the delineation of diarrhea phases) as fixed effect and individual calves as random effect:

whereYijkand Ymkrepresent the variable of interest;Aiis the fixed effect of calf age;Hjis the fixed effect of calf diarrheic status, defined as non-diarrheic or diarrheic status based on fecal score;Smis the fixed effect of calf diarrhea phases, defined as pre-diarrhea phase, diarrhea phase, or post-diarrhea phase based on the temporal changes of fecal scores of the same calves;Ikis the random effect of individual calves; andεijkandεmkare the residual error. Non-parametric Kruskal–Wallis test was used to assess the effects of age and diarrhea on bacterial relative abundance, and Dunn’s all-pairs rank comparison test withPadjusted by false discovery rate was used to conduct multiple comparisons. A significant change was declared withP< 0.05.

Fecal microbiota comparison among ages and identification of age?associated genera of bacteria

In trial 1, the overall fecal microbiotas between two ages were pairwise compared using analysis of similarity (ANOSIM) implemented in the Vegan package (version 2.5.3) [32] in R. WhenP< 0.05, the fecal microbiotas between two ages were considered completely different (R-value > 0.75), different (0.5

Random forest regression was used to identify the fecal bacterial genera that were associated with the age of the calves using the randomForest package (version 4.6.14) [33, 34] in R. The genus table of all the samples was the input data. The random forest algorithm was executed with the default parameters (ntree = 1000, default mtry of p/3, where p is the number of input genera (‘features’)). The importance of a genus was ranked in the order of its ‘feature importance’, with feature importance being the decrease in prediction accuracy (in percent) of the model when that genus was removed. To explore the ageassociated microbiota development, cross-validation was performed to estimate the optimal number of top-ranking age-associated genera required for prediction using the rfcv function implemented in the randomForest package in R. The identified genera were shown in a heatmap with their relative abundance. These genera were considered microbial markers of the respective ages.

Fecal microbiota comparison between diarrheic statuses, between diarrhea phases, and identification of their associated ASVs and modules

Fecal microbiotas between the diarrheic and non-diarrheic calves in trial 1 and between diarrhea phases (prediarrhea, diarrhea, and post-diarrhea) in trial 2 were compared using principal coordinates analysis (PCoA) and permutational multivariate analysis of variance (PERMANOVA/adonis) based on Bray–Curtis dissimilarity. When the fecal microbiotas of diarrheic calves were significantly different from those of age-matched nondiarrheic calves, in trial 1, and when the fecal microbiotas of the calves with diarrhea differed significantly from those at their pre-or post-diarrhea phase in trial 2, Linear discriminant analysis Effect Size (LEfSe) [35] and significance test with DESeq2 [36] were used to identify the ASVs that might be associated with one of the diarrheic statuses or one of the diarrhea phases. Of the ASVs with an LDA score > 2 in LEfSe or an adjustedP-value < 0.05 in DESeq2, those with a log2fold change > 1 (diarrheic/nondiarrheic calves, or diarrhea/pre-or post-diarrhea phase) were considered to be associated with diarrheic status or phase, whereas those with a log2fold change < – 1 (defined the same as above) were considered associated with non-diarrheic status or phase.

Co-occurrence patterns of the fecal microbiotas of the diarrheic and non-diarrheic calves (see the results for the delineation of the stages) in trial 1, or the fecal microbiota of the calves at different diarrhea phases were examined using the SparCC algorithm [37] with ASV count table as the input data. The pattern was visualized using the igraph package (version 1.2.5) [38] in R, and correlations with aP< 0.05 and a co-efficient R ≥ 0.5 or ≤ – 0.5 being considered positive and negative correlations, respectively. Modules of the co-occurrence patterns were generated using the walktrap algorithm [39] implemented in igraph. Modules with less than 3 nodes were deleted from the co-occurrence patterns. The identified ASVs associated with a diarrheic status or diarrhea phase were highlighted in the patterns. The modules aggregated with ASVs that were associated with diarrhea in trial 1 or with the diarrhea phase in trial 2 were considered as diarrhea modules and diarrhea-phase modules, respectively, whereas the modules formed with ASVs that were not with diarrhea in trial 1 or associated with pre-or post-diarrhea phases in trial 2 were considered as nondiarrhea modules or pre-or post-diarrhea modules.

Results

Trial 1

Development of fecal bacterial microbiota of calves and temporal microbial successions

In total, 15,283,464 quality-filtered amplicon sequences were obtained from 251 fecal samples (the fecal sample of calf Y03 on d 5 was not analyzed due to contamination with wheat straw) with an average of 60,890 ± 7193 (mean ± SD) sequences per sample. The sequencing depth coverage reached > 99.96% on average (99.89% to 100.00%). In the fecal samples collected from 1 day of age, 345 ASVs (referred to as species hereafter) on average per sample were identified with a Shannon diversity index of 4.03 (Fig. 1A). The number of observed species, Faith’s PD, Shannon diversity index, and evenness decreased from d 1 to 7 but recovered at d 9. Then, observed species, Faith’s PD, Shannon index, and evenness gradually increased, though with fluctuation, to about 509, 28.77, 5.03 and 0.81, respectively, at d 78 (Fig. 1A). Over this period, age significantly (P< 0.05) affected all the diversity metrics, whereas diarrheic status did not affect (P> 0.05) any of these four metrics.

Collectively, the ASVs were classified into 200 genera within 12 phyla. Bacteroidota, Firmicutes, Proteobacteria, and Fusobacteria each had a relative abundance > 1% in more than 60% of the fecal samples at any day (Fig. 1B). Of the 200 bacterial genera identified across all the fecal samples, 15 genera each had a relative abundance > 1% in at least 60% of the fecal samples at a single day. These genera includedAlloprevotella,Bacteroides,Escherichia/Shigella,Faecalibacterium,Fusobacterium,Acinetobacter,Prevotellaceae_UCG-003,Rikenellaceae_RC9_gut_group,Ruminococcaceae_UCG-005,Ruminococcaceae_UGC-010,Butyricicoccus,Lachnospiraceae_FCS020_group,Parabacteroides,Prevotella_9andSutterella.These predominant genera displayed temporal changes in relative abundance over the course of the trial (Fig. 1C).Alloprevotella,Faecalibacterium,Parabacteroides, andSutterellaincreased their relative abundance (P< 0.05) and maintained a higher abundance around d 15 to 58 compared to d 1 and then decreasing towards the end of the trial. Compared to d 1,Bacteroides,Escherichia/Shigella,Fusobacterium, andButyricicoccusincreased and then decreased their relative abundance sharply (P< 0.05) at around d 10.On the contrary,Prevotellaceae_UCG-003,Rikenellaceae_RC9_gut_group,Ruminococcaceae_UCG-005,Ruminococcaceae_UGC-010, andLachnospiraceae_FCS020_groupincreased their relative abundance (P< 0.05) during the later days of the trial and maintained or decreased their relative abundance until the end of the trial.Acinetobacterdiffered from all the other genera as it lost its initial high relative abundance dramatically by d 3 (P< 0.05) and never recovered.Prevotella_9only increased (P< 0.05) on d 9.

Composition and distribution of age?associated bacterial genera from birth to post?weaning

Pair-wise comparison of the fecal microbiotas at the ASV level among the 18 time points using ANOSIM revealed three age periods, with the first, second, and third age periods being from d 1 to 12, 15 to 63, and 58 to 78, respectively. The fecal microbiotas were similar (R-value < 0.5) within each age period but different between most of the two age periods (Table 1).

Thirty-five age-associated bacterial genera were identified by random forest regression (Fig. 2A and B), and they were distributed in 4 clusters (Fig. 2C) each corresponding to one of the age periods (Table 1). Of these age-associated genera,Klebsiella,Escherichia/Shigella,Enterococcus,Bacteroides,ButyricicoccusandMegamonasin the first cluster were predominant at the early age (d 1 to 12);Alloprevotella,Faecalibacterium,Intestinimonas,Paraprevotella,UBA1819,Lachnospiraceae_UCG-010,Pygmaiobacter, andSubdoligranulumin the second cluster were predominant at d 15 to 58;Blautia,Breznakia,Agathobacter,Anaeroplasma,Romboutsia,Erysipelotrichaceae_UCG-004,Prevotella_9, andSuccinivibrioin the third cluster andCandidatus_Stoquefichus,Parasutterella,Lachnospiraceae_NK4A136_group,Oscillibacter,Prevotellaceae_UCG-003,Family_XIII_AD3011_group,Lachnospiraceae_FCS020_group,Rikenellaceae_RC9_gut_group,Prevotellaceae_UCG-001,Ruminococcaceae_UCG-005,Ruminococcaceae_UCG-010,Negativibacillus, andTyzzerellain the fourth cluster were predominant at d 48 to 68 and d 58 to 78, respectively.

Table 1 A matrix of R-values of pair-wise comparison of the fecal microbiota of trial one at the ASV level using ANOSIM

(See figure on next page.)

Fig. 1 Dynamic changes of fecal bacterial microbiota of calves from birth to post-weaning (trial 1). Dynamic changes of alpha diversity metrics with the red boxes indicating the diarrhea peaks (A), bacterial phyla (B), and major genera (C) across ages. Only the phyla and genera each with a relative abundance > 1% in at least 60% of samples at any single age were shown. The relative abundance significantly differing from that of d 1 is indicated with a * (P< 0.05)

Fig. 1 (See legend on previous page.)

Diarrhea characteristics of the study cohort

Over the course of the trial 1, all the calves had fecal score ≥ 3 at least one sampling day (Fig. 3A). Based on the fecal scores of all the calves, the 78 d of trial was divided into five stages: stage 1: d 1 to 7, before the first diarrhea peak; stage 2: d 9 to 15, the first diarrhea peak; stage 3: d 18 to 38, the stage after the first peak but before the second diarrhea peak; stage 4: d 48 to 68, the second diarrhea peak; and stage 5: d 75 to 78, the stage after the second diarrhea peak. All the 14 calves were non-diarrheic in stages 1, 3, and 5. Pathogen detection showed that all the 20 diarrheic fecal samples (9 in the first and 11 in the second diarrhea peaks) wereSalmonellaspp. and BCoV negative, but eight wereE. coliK99+positive (2 in the first and 6 in the second diarrhea peaks). One of the diarrheic samples (in the second diarrhea peak) was both BRV andE. coliK99+positive (Fig. 3A).

ASVs and modules associated with diarrheic status

The fecal microbiotas differed among the five stages (P< 0.001, Fig. 3B). The fecal microbiotas of diarrheic and non-diarrheic calves did not differ (P= 0.450) in stage 2 (Fig. 3C) but did differ (P= 0.004) in stage 4 (Fig. 3D). Based on fold change and LEfSe or DESeq2 analysis, 147 diarrheic status-associated ASVs were identified in stage 4 including 91 diarrhea-associated ASVs and 56 non-diarrhea-associated ASVs (Additional file 2: Fig. S2). The diarrhea-associated ASVs mainly consisted ofBacteroides(15 ASVs),Ruminococcaceae_UCG-005(8 ASVs),Ruminococcaceae_UCG-010(6 ASVs),Ruminococcaceae_UCG-013(4 ASVs) andLachnospiraceae_FCS020_group(4 ASVs). The non-diarrhea-associated ASVs mainly consisted ofRuminococcaceae_UCG-005(8 ASVs),Flavonifractor(4 ASVs) andBacteroides(3 ASVs).

The co-occurrence pattern for stage 4 had 195 nodes, 603 edges, and 17 modules (Fig. 4A, Additional file 3: Table S1). Nineteen diarrhea-associated ASVs and 5 nondiarrhea-associated ASVs showed up in the co-occurrence pattern. The diarrhea-associated ASVs includingPrevotellaceae_UCG-003(ASV11, ASV43),Ruminococcaceae_UCG-010(ASV427, ASV574),Butyricicoccus(ASV130),Lachnospiraceae_FCS020_group(ASV234) and unclassified Ruminococcaceae (ASV557) were scattered without aggregation in any of the modules.Bacteroides(ASV170, ASV206),Rikenellaceae_RC9_gut_group(ASV33),Ruminococcaceae_UCG-010(ASV196),Lachnospiraceae_FCS020_group(ASV174), unclassified Bacteroidales (ASV177), unclassified Barnesiellaceae (ASV22) and unclassified Ruminococcaceae (ASV340) were aggregated in diarrhea module 1 (D-M1), whileBacteroides(ASV39, ASV41, ASV145, and ASV259) were aggregated in diarrhea module 2 (D-M2). The nondiarrhea-associated Muribaculaceae (ASV28 and ASV44) andUBA1819(ASV151) were aggregated in non-diarrhea module (ND-M). The non-diarrhea-associatedBarnesiella(ASV497) andRuminococcaceae_UCG-005(ASV254) were aggregated in D-M1, and they formed negative correlations with the ASVs in D-M1. The ASVs in D-M1 and D-M2 had a higher (P< 0.05) total relative abundance in diarrheic calves than in non-diarrheic calves (Fig. 4B). The D-M1 was mainly occupied byBacteroides, unclassified Barnesiellaceae,andRikenellaceae_RC9_gut_group, with 1.49%, 1.53% and 0.75% relative abundance, respectively; D-M2 was occupied byFusobacteriumandBacteroides, with 1.67% and 1.58% relative abundance, respectively, but the ND-M was mainly occupied by unclassified Muribaculaceae,Prevotellaceae_Ga6A1_group,Prevotellaceae_UCG-003, andPrevotella_9, with 1.90%, 0.69%, 0.57%, and 0.42% relative abundance, respectively (Fig. 4C).

Fig. 2 Age-associated bacterial genera from birth to post-weaning (trial 1). The top 35 age-associated bacterial genera ranked by importance to the accuracy of the random forest regression model (A). Ten-fold cross-validation error as a function of the number of input genera was used to regress against the chronologic age of calves. The dotted line indicates the 35 genera used in the model (B). Heatmap of the top 35 age-associated genera and the clusters they formed based on their relative abundance across ages (C)

Trial 2

Diarrhea characteristics of the study cohort

Over the course of the trial 2, all the calves had fecal score ≥ 3 at least 1 d. (Fig. 5A). In the 43 calves, 4, 20, 33, 31, 30, 18, 18, 16, 15, 4 and 5 calves had diarrhea on d 8 to 18 respectively. In which, 32 calves had one or more episodes of diarrhea 1 to 5 d after recovering from a previous episode (Fig. 5A). Except of d 9, 10 and 11, the microbial composition had no significant difference between diarrheic and non-diarrheic calves within the same age (Additional file 4: Fig. S3). At the same age, calves were in different days of diarrhea, so microbial changes were also out of synchronization (Fig. S3). So based on the temporal changes of the fecal scores, the fecal samples were divided into four phases: pre-diarrhea (when fecal score < 3), diarrhea (fecal scores ≥ 3 consecutively), post-diarrhea (fecal score falling below 3), and volatility (fecal score rising to ≥ 3 after 1 to 5 d below 3) (Fig. 5B). To track the microbial changes from the pre-diarrhea to post-diarrhea, samples of calves of different ages with or without diarrhea were combined into phases for analysis.

Changes in fecal bacterial communities of calves from pre? to post?diarrhea

In total, 18,706,378 quality-filtered amplicon sequences were obtained from 340 fecal samples (the fecal samples of pre-diarrhea, diarrhea, and post-diarrhea phases) with an average of 55,019 ± 10,329 (mean ± SD) sequences per sample. The sequencing depth coverage reached > 99.97% on average (99.91% to 99.99%). The fecal microbiotas differed among pre-diarrhea, diarrhea, and post-diarrhea phases (P< 0.01, Fig. 5C).

The lowest number of observed species (on average 117 ASVs per sample) and Shannon index (3.2) were found in the fecal samples of the diarrhea phase (Fig. 6A). Both metrics increased (P< 0.05) in the post-diarrhea phase. Of the major phyla and genera (each with a relative abundance > 1% in > 60% of the fecal samples at any phase), the phylum Bacteroidota decreased (P< 0.05), while the phylum Fusobacteria increased (P< 0.05) from the pre-diarrhea to diarrhea phase (Fig. 6B), and the generaBacteroides,Butyricicoccus, andLachnoclostridiumdecreased (P< 0.05), whileClostridium_sensu_stricto_1andFusobacteriumincreased (P< 0.05) during the phase transition (Fig. 6C). In the post-diarrhea phase, the phyla Bacteroidota and Fusobacteria and the genusFusobacteriumtended to return their relative abundance to the pre-diarrhea level (P< 0.05), while the genusClostridium_sensu_stricto_1lost the relative abundance gained in the diarrhea phase (P> 0.05), and the generaBacteroides,Butyricicoccus, andLachnoclostridiummaintained their relative abundance of the diarrhea phase. Although not having changed their relative abundance from the prediarrhea phase to the diarrhea phase, the generaFaecalibacterium,Prevotella_2, andAlloprevotellaincreased their relative abundance (P< 0.05), while the phylum Epsilonbacteraeota and the genusCampylobacterdecreased their relative abundance (P< 0.05) in the postdiarrhea phase.

ASVs and modules associated with the diarrhea phase

Comparison of the fecal microbiotas between pre-diarrhea and diarrhea phases revealed (based on fold change and analysis using LEfSe or DESeq2) 32 pre-diarrheaassociated ASVs and 29 diarrhea phase-associated ASVs (Additional file 5: Fig. S4). The pre-diarrhea phase-associated ASVs mainly consisted ofBacteroides(12 ASVs),

Parabacteroides(4 ASVs),Butyricicoccus(3 ASVs), andFlavonifractor(3 ASVs). The diarrhea phase-associated ASVs mainly consisted ofClostridium_sensu_stricto_1(8 ASVs), unclassified Fusobacteriaceae (7 ASVs),Fusobacterium(6 ASVs), and unclassified Enterobacteriaceae (3 ASVs). The co-occurrence pattern of the pre-diarrhea and diarrhea phases had 44 nodes, 78 edges, and 4 modules (Fig. 7A). The pre-diarrhea and diarrhea phase-associated modules (pre-D-M and D-P-M, respectively) were aggregated with six and 16 pre-diarrhea and diarrhea phase-associated ASVs, respectively, and both modules were independent. Module 3 and module 4 were aggregated with 14 and 5 ASVs, respectively. Both module 3 and module 4 had no diarrhea or pre-diarrhea phase-associated ASVs, but module 4 formed a negative correlation with the D-P-M throughPrevotella_2(ASV26) andClostridium_sensu_stricto_1(ASV48). The ASVs in the D-P-M had a higher (P< 0.05) total relative abundance in the diarrhea phase than in the pre-diarrhea phase, whereas the opposite was true for the ASVs in pre-D-M (Fig. 7B). The D-P-M was occupied by ASVs assigned toClostridium_sensu_stricto_1andFusobacterium, with a relative abundance of 29.11% and 5.42%, respectively. The pre-D-M was only occupied by ASVs ofBacteroides, with a relative abundance about 11.77%. The ASVs of module 3 mainly consisted ofAllprevotella,BacteroidesandPrevotella_9, with 2.94%, 1.65% and 0.55% relative abundance, respectively. The ASVs of module 4 were classified toPrevotella_2andLachnospiraceae_UCG-004, with 4.48% and 1.05% relative abundance, respectively (Fig. 7C).

Fig. 4 Co-occurrent pattern of the fecal microbiota at stage 4 in trial 1. Pattern showing the diarrhea status-associated ASVs and modules (A). Mean relative abundance of diarrhea status-associated modules in diarrheic and non-diarrheic calves (B). ASVs relative abundance at genus level in different modules (C). ASV254 and ASV497 in module 2 were not included in plots (B) or (C)

Fig. 5 Transition phase division of the 43 calves based on fecal score in trial 2. Fecal score of the 43 calves from d 8 to 18 (A). Phase division based on samples fecal score (B). Principal coordinates analysis (PCoA) plot of fecal microbiotas among pre-diarrhea, diarrhea, and post-diarrhea phases (C)

When compared the ASVs between diarrhea and post-diarrhea phases, 44 post-diarrhea phase-associated ASVs and 23 diarrhea phase-associated ASVs were identified based on fold change and analysis using LEfSe or DESeq2 (Additional file 4: Fig. S4). The post-diarrhea phase-associated ASVs mainly consisted ofPrevotella_9(13 ASVs),Alloprevotella(8 ASVs),Bacteroides(4 ASVs), andCollinsella(3 ASVs). The diarrhea phase-associated ASVs mainly consisted ofClostridium_sensu_stricto_1(6 ASVs), unclassified Fusobacteriaceae (7 ASVs), and unclassified Enterobacteriaceae (4 ASVs). The co-occurrence pattern constructed of the fecal microbiotas at post-diarrhea and diarrhea phases had 57 nodes, 116 edges, and 5 modules (Fig. 7D). The D-P-M was aggregated with 11 diarrhea phase-associated ASVs and two post-diarrhea phase-associated ASVs (ASV4 and ASV26, both belonging toPrevotella_2), both of which had a negative relationship withClostridium_sensu_stricto_1(ASV48), a diarrhea phase-associated ASV and the keystone of D-P-M. Two post-diarrhea phase-associated modules (post-D-M1 and post-D-M2) were identified in the co-occurrence network, and post-D-M1 and post-D-M2 were aggregated with 11 and 5 post-diarrhea phase-associated ASVs, respectively. The composition of post-D-M1 was similar to that of module 3 in Fig. 7A, but post-D-M1 had more betweenness and close centrality than module 3 and formed a negative correlation with D-P-M throughAlloprevotella(ASV14) andPrevotella_9(ASV43), both of which negatively associated withClostridium_sensu_stricto_1(ASV48) in the D-P-M. The ASV composition of module 2 was the same as that of pre-D-M in Fig. 7A, and module 2 was standalone from the other four modules (Fig. 7D). Module 5 consisted of six ASVs. Although having no diarrhea and post-diarrhea phase-associated ASVs, module 5 positively connected with D-P-M throughTyzzerella_4(ASV23) andClostridium_sensu_stricto_1(ASV48). The D-P-M was occupied byFusobacterium,Prevotella_2andClostridium_sensu_stricto_1with a relative abundance of 30.94%, 5.34% and 5.19%, respectively (Fig. 7F). The post-D-M1 was mainly occupied byAlloprevotella, BacteroidesandPrevotella_9, with a combined relative abundance of 4.80%, 2.32%, and 0.88%, respectively. The post-D-M2 was only occupied byPrevotella_9, with a relative abundance of 0.21%.

Fig. 6 Fecal microbiota changes from pre-diarrhea, diarrhea, and post-diarrhea phases in trial 2. Two alpha diversity metrics of the fecal microbiota (A), Relative abundance of predominant bacterial phyla (B) and genera (C). Only the phyla and genera each with a relative abundance > 1% in at least 60% of samples in a single phase were shown. Significant differences between two phases were indicated with a * (P < 0.05)

Fig. 7 Co-occurrent pattern of the fecal microbiota among three phases in trial 2. Pattern among pre-diarrhea and diarrhea phases (A). Mean relative abundance of phase-associated modules in pre-diarrhea and diarrhea (B). ASVs relative abundance at genus level in pre-diarrhea and diarrhea modules (C). Pattern among diarrhea and post-diarrhea phases (D). Mean relative abundance of phase-associated modules in diarrhea and post-diarrhea (E). ASVs relative abundance at genus level in diarrhea and post-diarrhea modules (F). ASV4 and ASV26 in diarrhea-M were divided into post-D-M1 in plot (E) and (F)

Discussion

A better understanding of the fecal microbiota in diarrheic and non-diarrheic calves can inform improved treatment and prevention strategies. Fecal microbiotas between diarrheic and non-diarrheic calves have been compared after interventions, such as trehalose supplementation [40], feeding waste milk containing antibiotic residues [41] or supplemented with sodium humate and glutamine combination [42], or single species [43] or multispecies probiotics [44]. These interventions affected the developing process of gut microbiota and increased the abundance ofBifidobacteriumandLactobacillus, which might conceal the natural development of resistance to pathogenic colonization in preweaned calves. Without these types of interventions, Kim et al. [20] described the ability of a fecal microbiota transplantation (inclusion was collected from the health calves of the similar age to diarrheic calves) could ameliorate diarrhea and restore gut microbial composition in pre-weaning calves. Most of these comparative studies focused on pre-weaning calves. In the present study, we first examined the dynamic development of the gut microbiota (represented by fecal microbiota) and diarrhea occurrence in dairy calves from birth to 15 d postweaning with frequently fecal sampling in trial 1, which allowed us to identify two diarrhea peaks. Then we analyzed the fecal microbiota and diarrhea occurrence of another group of dairy calves over 11 d corresponding to the first diarrhea peak with fecal samples collected daily. Trial 1 helped test our hypothesis that some gut microbes might have resistance to dysbiotic process with calf diarrhea by dictating the microbial co-occurrence patterns during weaning, while trial 2 allowed us to examine in details as the calves transitioned from pre-diarrhea to diarrhea and then to post-diarrhea phases shortly after birth.

The development of the gut microbiota appeared to have three age periods with birth and weaning as the separatrices [45, 46]. As antibiotics had been used from d 1 to 5, the drastically decreasing of bacterial richness after birth (Fig. 1A) might be the effects of tylosin tartrate and sulfadimidine. However, in the recent study of neonatal dairy calves, without antibiotics, a similar decrease of bacterial richness was reported by Kim et al. [47] and Klein-J?bstl et al. [48]. In the first two weeks after birth, the underdeveloped gut was colonized primarily by facultatively anaerobic microbes, especiallyEscherichia/Shigella, to render the intestinal environment suitable for anaerobic intestinal microbes to colonize, so the decrease of bacterial richness might be a result of bacterial adaptation.Escherichia/Shigellahad a low relative abundance in the fecal samples at any day after birth to 15 d after weaning, but its relative abundance increased up to 45.36% and 31.51% on d 3 and 5, respectively, and then decreased to around 2% since d 15. The increased colonization byEscherichia/Shigellain young calves has been described previously [49], and it might explain the susceptibility of calves to diarrhea caused byE. coli.Klebsiella, which belongs to the same family, Enterobacteriaceae,asEscherichia/Shigella, also appeared to be an age-associated bacterial genus for this period.Klebsiella pneumoniae[50] andKlebsiella oxytoca[51] are opportunistic pathogenic in humans and are associated with increased infection mortality rate, particularly in immunocompromised individuals, neonates, and the elderly. However, infections in calves caused byKlebsiellahave not been reported frequently. Glantz and Jacks reported thatKlebsiellaspp. occurred naturally in calves, and they might be responsible for some mortality [52]. Komatsu et al. reported fatal suppurative meningoencephalitis caused byK. pneumoniaein two calves [53]. Aslan et al. reported thatK. pneumoniaecould be isolated in calves suffering from respiratory tract infection, which was not cured by florfenicol [54]. In the present study, we observed changes of relative abundance ofKlebsiella, decreasing from 1.07% and 2.09% on d 1 and 3, respectively, to less than 0.1% at d 5. Tylosin tartrate and sulfadimine are both broad spectrum antibiotics exerting their antimicrobial action by inhibiting the bacterial protein synthesis [55], and competing with para-aminobenzoic acid for dihydropteroate synthase [56], respectively. The decrease inKlebsiellamight be a result of tylosin tartrate [57]. The peak ofKlebsiellaabundance did not correspond to a diarrhea peak, but it should be given special attention in calf industry because of its pathogenic significance.

Bacteroideswas the most abundant genus through the whole period, but its relative abundance increased sharply to 41.96%, 49.66%, and 40.38% at d 5, 7, and 9, respectively, so it could be an age-associated bacterial genus for the first two weeks after birth. Its high abundance was attributed to its stronger saccharolytic ability thanPrevotellain the gut of young calves [58, 59]. The acetate produced byBacteroidescould be consumed by other bacteria, such asButyricicoccusandMegamonas, to produce butyrate and propionate [60], both of which are the main source of energy for intestinal epithelial cells, and butyrate can also inhibit the signaling pathways of pro-inflammatory cytokines [61], enhance intestinal barrier function by increasing mucin secretion and enhancing the tight-junctions [62]. Thus, the colonization by these microbes in the gut of young calves might facilitate the establishment of a functional gut. From the third week onward, weaning separated the microbiota characteristic of the two periods. Prior to weaning for more than 40 d, the microbiotas remained similar, suggesting that the gut might have established a stable functional microbiota over that period, with limited recruitment. However, with the transition from liquid to total solid feed, some of the gut bacteria displayed considerable changes, as exemplified by the replacement ofAlloprevotella,Faecalibacterium, andParabacteroidesbyPrevotellaceae_UCG-003,Rikenellaceae_RC9_gut_group,Ruminococcaceae_UCG-005,Ruminococcaceae_UGC-010, andLachnospiraceae_FCS020_group. These five genera are within families that contain species capable of utilizing structural polysaccharides, and this replacement might facilitate the degradation of structural polysaccharides and the stability of the gut microbiota. Future research can help identify the species of these five genera and characterize their functions. Two diarrhea peaks were observed during the three developing periods of the gut microbiota, and less than half of the tested diarrheic fecal samples were pathogen positive, which suggests that microbiota homeostasis may be more important in preventing diarrhea than directly killing pathogens. Comparing fecal microbiota transplantation and antibiotic treatment for ameliorating calf diarrhea, Kim et al. [20] confirmed that gut microbial manipulation could offer another therapeutic paradigm, beyond antibiotic based therapies. Our data also suggest that prophylaxis/preventions with probiotics should be better administered in d 1 to 7 (stage 1, before the first diarrhea peak) and d 18 to 38 (stage 3, the stage after the first peak but before the second diarrhea peak). Supplementation of newborn calves withLactobacillusandBifidobacterium[63] orFaecalibacterium prausnitzii[64] within the first 7 d of life decreased diarrhea, but no study was found in the literature that had tested probiotic supplementation in stage 3. To promote the natural development of the gut microbiotas during the transition, the potential probiotics in the gut of calves should be identified.

Comparison of the fecal microbiota using LEfSe and the subsequent identification of the correlations between the differential bacterial genera and the core bacterial genera in the gut is a well-trodden path to reveal the effects of diarrhea on the microbiota and find potential probiotics in neonatal dairy calves for preventing diarrhea [65]. But many standard correlation analyses may lead to misleading results because 16S rRNA gene profiling data are sparse and compositional [37]. SparCC, which is tailored to the compositional and sparse features of genomic survey data and allows for inference of correlations between genes or species, has been used to elucidate the networks of interaction among microbial species living in or on the human body [37, 66]. Therefore, to reduce the incidence of false positive results, three levels of evaluation criteria including co-occurrence patterns examined using the SparCC algorithm (correlations with aP< 0.05 and a co-efficient R ≥ 0.5 or ≤ – 0.5 being considered positive and negative correlations, respectively), |log2fold change|> 1, LDA score > 2 in LEfSe or adjustedP-value < 0.05 in DESeq2 were used in the present study. During the weaning, some ASVs and modules were associated with diarrhea, while some were associated with non-diarrhea (Fig. 4 and Fig. S2).

With the three levels of evaluation criteria, the identified diarrhea-associated ASVs were aggregated inBacteroides(ASV39, ASV41, ASV145, ASV170, ASV206, and ASV259).These members might be biomarkers of diarrhea risk. Muribaculaceae (ASV28 and ASV44), UBA1819(ASV151),Barnesiella(ASV497),andRuminococcaceae_UCG-005(ASV254) were non-diarrhea-associated ASVs in the co-occurrence pattern, and ASV28 and ASV497 had a direct inhibitory relationship with the members of D-M1 (Fig. 4A). Muribaculaceae, which was previously assigned as family S24-7 or Homeothermaceae, is a common and abundant family of symbiotic bacteria in the gut and specialized in fermenting complex carbohydrates [67]. It responded most positively to acarbose treatment for diabetes [68] and was linked to longevity [69].Barnesiellabelongs to the family Porphyromonadaceae within the phylum Bacteroidota. It was found to suppress the growth of intestinal vancomycinresistantEnterococcus[70]. Members of Ruminococcaceae are mostly butyrate-producing bacteria. Weese et al. suggested that Firmicutes (particularly Lachnospiraceae and Ruminococcaceae)/Proteobacteria ratio might be used to potentially predict and prevent colic [71]. Although some members of Ruminococcacea were associated with diarrhea, both ASV151 and ASV254 (both assigned to Ruminococcaceae) were associated with non-diarrhea. So Muribaculaceae (ASV28 and ASV44),UBA1819(ASV151),Barnesiella(ASV497),andRuminococcaceae_UCG-005(ASV254) might be used as potential probiotics for this specific commercial farm, with this specific microbial colonization patterns. The sequences of these ASVs were included in Table 2, and these ASVs can be verified in future studies. Fusobacteriaceae dominated D-M2 (Fig. 4C). Fusobacteriaceae was reported to have high relative abundances in dairy calves suffering from diarrhea, either infected [72] or uninfected [65], which indicates that module analysis can help identify bacteria associated with diarrhea or otherwise. With this module analysis, Muribaculaceae andPrevotellawere identified as the core microbiota resisting diarrhea in weaning calves. Kim et al. reported that substituting fermented soybean meal (FSBM) for soybean meal (SBM) at 5% level in calf starter reduced the incidence of diarrhea and improved immunocompetence in neonatal calves after microbial infection [73]. The reason for the positive effects of FSBM on immunocompetence was not reported, but in a recent study, Feizi et al. found that FSBM increased the abundance ofPrevotella ruminicolain the rumen of dairy calves [74]. Essential oils showed similar effects on dairy calves, increasing the Prevotellaceae abundance in the rumen [75] and decreasing the morbidity of neonatal diarrhea among pre-weaning calves [76]. Furthermore, a recent study reported that sodium humate and glutamine in combination also elevated the abundance ofP. ruminicolain the rectum while reducing diarrhea incidence among dairy calves during the weaning period [42]. Tap et al. reported that Prevotellaceae enterotype was less susceptible to irritable bowel syndrome (IBS) compared with Bacteroidaceae enterotype [77]. Therefore, future research is warranted to investigate the relationship between calf diarrhea andPrevotellaas a genus and its species.Prevotellamay also be explored for its preventative ability to reduce calf diarrhea.

Consistent changes in relative abundance ofBacteroides,Butyricicoccus,Faecalibacterium,Alloprevotella, andFusobacteriumwere observed in both trial 1 and trial 2 (Fig. 1, 2 and 6). When the fecal microbiota was examined in detail as the calves transitioned from prediarrhea to diarrhea and then to post-diarrhea phases in trial 2, the peak of bothClostridium_sensu_stricto_1andFusobacteriumcoincided with the peak of diarrhea. It has been reported thatClostridium_sensu_stricto_1might cause epithelial inflammation in piglets [78] and stunting in infants (defined as height-for-age Z score equal to or lower than – 2, [79]). Therefore, research is needed to further investigate these two genera withrespect to their role in calf diarrhea. Our co-occurrence analysis (Fig. 7) showed that the post-D-M1 might be a driver of diarrhea recovery because of the close interaction between its constituent members and the inhibitory relationship with D-P-M.Prevotella_2andAlloprevotellaincreased their relative abundance in the post-diarrhea phase (Fig. 6C), and they dominated the post-D-M1 (Fig. 7F), which supports the importance of Prevotellaceae in resisting calf diarrhea. In particular,Prevotella_2(ASV4 and ASV26),Alloprevotella(ASV14) andPrevotella_9(ASV43), their sequences were included in Table 2, might be potential probiotics for preventing diarrhea in early stage. It should be noted that although most of the constituent members of post-D-M1 were detected before diarrhea (module 3 and 4 in Fig. 7A), their betweenness and close centrality increased in post-D-M1. This suggests that the interactions among different bacteria might play an important role in maintaining intestinal homeostasis in the gut.

Table 2 The sequences of ASVs with potentials being probiotics and biomarker of diarrhea risk

Table 2 (continued)

These potential probiotics may be supplemented in the first week after birth to prevent diarrhea, and fiber diets [80] or FSBM [71] may improve their efficacy.Clostridium_sensu_stricto_1(ASV48, Fig. 7, Table 2), was negatively related withPrevotella_2(ASV4 and ASV26),Alloprevotella(ASV14) andPrevotella_9(ASV43). All of the four ASVs established the negative relationship between the post-D-M1 and D-P-M. Thus,Clostridium_sensu_stricto_1(ASV48) might be a biomarker of diarrhea risk in the early stage.

Conclusions

In conclusion, microbial successions of the gut microbiome in dairy calves were rapid, and daily sampling is needed to capture the rapid dynamic gut microbial successions. Promoting indigenousPrevotellaand Muribaculaceae might be a new strategy to reduce the incidence of diarrhea in neonatal calves and help calves to go through the weaning transition smoothly.Prevotella_2(ASV4 and ASV26),Prevotella_9(ASV43),Alloprevotella(AVS14), unclassified Muribaculaceae (ASV28 and ASV44),UBA1819(ASV151),Ruminococcaceae_UCG-005(ASV254), andBarnesiella(ASV497) might be used as probiotics to reduce or prevent calf diarrhea;Clostridium_sensu_stricto_1(ASV48) might be a useful biomarker of diarrhea risk in this large-scale dairy farm locating in subtropical monsoon climate zone with automated milk feeder system.

Abbreviations

E. coliK99+:Escherichia coliK99+; BRV: Bovine rotavirus; BCoV: Bovine coronavirus; ASV: Amplicon sequence variant; Faith’s PD: Faith’s phylogenetic diversity; ANOSIM: Aanalysis of similarity; PCoA: Principal coordinates analysis; PERMANOVA: Permutational multivariate analysis of variance; LEfSe: Linear discriminant analysis Effect Size; D-M1: Diarrhea module 1; D-M2: Diarrhea module 2; ND-M: Non-diarrhea module; pre-D-M: Pre-diarrhea phaseassociated module; D-P-M: Diarrhea phase-associated module; post-D-M1/2: Post-diarrhea phase-associated modules.

Supplementary Information

The online version contains supplementary material available at https://doi.org/10. 1186/s40104-022-00758-4.

Additional file 1: Fig. S1.A schematic showing the experimental design, milk feeding, and fecal sample collection of the two trials.

Additional file 2: Fig. S2.Heatmap of the ASVs associated with diarrheic status in stage 4 of trial 1. The ASVs were identified based on fold change and analysis using LEfSe or DESeq2.

Additional file 3: Table S1.The ASVs in co-occurrence pattern of trial 1 stage 4.

Additional file 4: Fig. S3.Principal coordinates analysis (PCoA) plot of fecal microbiotas among ages in trial 2. The fecal microbiotas differences (P values) between diarrheic and non-diarrheic calves within the same age were compared.

Additional file 5: Fig. S4.Heatmap of the ASVs that identified to be diarrheic status transition-associated in trial 2. The ASVs were identified based on fold change and analysis using LEfSe or DESeq2.

Acknowledgements

We thank the College Experimental Teaching Center, College of Animal Sciences, Zhejiang University for experimental equipment surport and assistance. We also thank our colleagues Guangyu Zhang, Lei Cao, Weibing Shi, Shaobo Yu, Zhibo Liu and Peng Wu at the Institute of Dairy Science, College of Animal Sciences, Zhejiang University for their assistance with sampling.

Authors’ contributions

JW and HC designed the study and drafted the manuscript. HC, YL and KH conducted the animal trials and collected the fecal samples. HC performed the DNA extraction and pathogen detection. HC and BY performed the data analyses and visualization. YZ provided experiment assistance. ZY provided critical guidance in the interpretation of results and revised the manuscript. JW supervised the study and revised the manuscript. The authors read and approved the final manuscript.

Funding

This work was supported by the National Key Research and Development Program of China (2017YFD0500502).

Availability of data and material

The raw sequencing data generated in this study are publicly available in NCBI Sequence Read Archive (http://www. ncbi. nim. nih. gov/sra) under the accession number PRJNA716761.

Declarations

Ethics approval and consent to participate

All experimental protocols used in the current study were approved by the Animal Care and Use Committee of Zhejiang University (Protocol number: ZJU-20262), and all experimental procedures were performed following the approved protocols.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Author details

1Institute of Dairy Science, College of Animal Sciences, Zhejiang University, Hangzhou, China.2MoE Key Laboratory of Molecular Animal Nutrition, Zhejiang University, Hangzhou, China.3Department of Animal Sciences, The Ohio State University, Columbus, OH, USA.

Received: 15 March 2022 Accepted: 13 July 2022