Difference in microbiota between the frozen and freeze-dried microorganisms
To confirm that the process used to prepare the fecal microorganisms separated from the donor’s feces did not alter the microbiota, frozen (F-FMT, n = 10) and freeze-dried (FD-FMT, n = 9) microorganisms were prepared for a total of 19 FMT studies (Fig. 1 and Supplementary Fig. 1). To demonstrate the FMT as a safe and virulence factor-free therapeutic, the fecal material obtained from the donor calf selected for FMT was initially confirmed to be pathogens-free (Fig. 2A). The donor’s ages did not differ between the calves used for F-FMT (50.9 ± 30.06 d old) and those used for FD-FMT (40.78 ± 15.09 d old) (Fig. 2B). The storage period starting from sample preparation to FMT was also identical between F-FMT (31.1 ± 16.68 d) and FD-FMT (39.45 ± 7.74 d) (Fig. 2C). The Procrustes test10, which visualizes the superimposition of sample coordinates for ordination analysis used to compare the microbial taxonomic profiles before and after sample preparation in each season, showed that the microbial composition of the frozen and FD microorganisms had significantly correlated with that of the intact feces (R2 = 0.6, p < 0.03 for frozen microorganisms, R2 = 0.5, p < 0.02 for FD microorganisms) (Figs. 2D and 2E). No significant differences were observed in the alpha diversity parameters, including Shannon entropy, Pielou_eveness, observed OTUs, faith_PD (Fig. 2F), and/or in the beta diversity based on weighted Uifrac distance and pairwise Permutational multivariate analysis of variance (PERMANOVA) before and after sample preparation or in both trials conducted between 2020–2021 and 2021–2022 (Fig. 2G and Supplementary Table 1).
Using BugBase11, a microbiome analysis tool that predicts the high-level phenotypes present in the microbiome samples, the proportion of each microbiome sample that includes the facultative anaerobic spp., Gram-negative spp., and gram-positive spp., and biofilm forming, and potentially pathogenic microorganisms was identical before and after sample preparation and between 2020–2021 and 2021–2022 (Fig. 2H). However, the number of aerobic microorganisms was slightly reduced after samples preparation regardless of the experimental seasons (Fig. 2H). The functional prediction of microbial communities in the isolated donor microbiota was carried out with Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2) using the 16S rRNA gene sequencing data12. The results of the principal coordinates analysis (PCoA) and pairwise PERMANOVA test conducted on Bray–Curtis distance of the predicted gene family for the enzyme commission (EC) numbers confirmed that the frozen microorganisms and their respective intact feces did not significantly differ (p = 0.585). The same phenomenon was also observed between the FD microorganisms and their respective intact feces (p = 0.960) (Fig. 2I and Supplementary Table 2). These results obtained by a metagenomic analysis indicated that the frozen and FD microorganisms were identical to their respective intact feces.
Efficacy of F-FMT and FD-FMT in treatment of diarrhea in the recipient’s calves
FMT studies were conducted with either the frozen or FD microorganisms in 2020–2021 and 2021–2022, respectively. There was no significant difference observed in ages of the recipient groups selected for two F-FMT and FD-FMT (Fig. 3A). However, in context of donors, there was a significant difference observed between the donor and recipients in both F-FMT (50.9 ± 30.06 vs 20.9 ± 19.90; p < 0.05) and FD-FMT (40.78 ± 15.9 vs 16.11 ± 12.04, p < 0.05) (Fig. 3B). In both trials of F-FMT and FD-FMT, they showed a complete reduction in the diarrheal scores (Fig. 3C) and in the fecal water content of the 19 recipients (Fig. 3D). The classical pathogen tests used for the detection of Rotavirus, coronavirus, E. coli, C. parvum, and C. perfringens in the recipient feces, revealed that C. parvum (11/19) and C. perfringens (8/19) pathogens were frequently detected regardless of the experimental seasons (Fig. 3E). However, Rotavirus (4/19) and coccidia (0/19) were rarely detected. Consistent with the previous study6, C. perfringens was detected in feces of several cases (6/19) even after recovery from diarrhea, whereas C. parvum was barely detected in most cases (3/19), regardless of the used strategy of bacterial preparation from the donor feces (Fig. 3E). In agreement with the previous study6, the blood tests used to measure the biochemical indicators demonstrated an increase in the total cholesterol level with F-FMT in 2020–2021 and FD-FMT in 2021–2022. The gamma-glutamyl transferase level decreased in both FMT treatments, whereas the other indicators [i.e., total protein, albumin, white blood cells, red blood cells, hemoglobin (Hb), hematocrit (Ht), mean corpuscular volume, mean corpuscular hemoglobin, and mean corpuscular hemoglobin concentration] remained unchanged (Fig. 3F). These results indicated that FMT was clinically effective in treating diarrhea in the recipient calves, regardless of how the donor feces were prepared (i.e., frozen or FD).
Efficacy of F-FMT and FD-FMT in altering microbiome in the recipients
F-FMT and FD-FMT showed similar clinical efficacy in treating diarrhea in the recipient calves; however, the differences between them were further evaluated to determine whether F-FMT or FD-FMT was superior in transforming the microbial community to the healthy conditions. Metagenomic analysis at the phylum level showed that the level of Bacteroidetes increased within 7 days after both F-FMT and FD-FMT, whereas the level of Proteobacteria decreased significantly (Fig. 4A and Supplementary Fig. 2). At the family level, the levels of Paraprevotellaceae and Prevotellaceae belonging to the phylum Bacteroidetes increased significantly within 7 days after both F-FMT and FD-FMT, whereas that of Enterobacteriaceae belonging to the phylum Proteobacteria decreased significantly (Fig. 4B and Supplementary Fig. 3). When considering the major genera, the levels of Fecalibacterium and Prevotella increased significantly within 7 days after both F-FMT and FD-FMT, whereas those of Bacteroides, Camphylobacter, and Veinollea decreased significantly (Fig. 4C and Supplementary Fig. 4). These results confirmed that the relative abundance of the major microbial taxa displayed similar patterns of change in both FMT treatments, regardless of the differences in samples preparation.
Linear discriminant analysis (LDA) of effect size (LEfSe) (LDA score > 2) was conducted to examine the overall bacterial features of those differentially represented before and after F-FMT and FD-FMT. Notable changes in microbial taxa were observed in FD-FMT compared with F-FMT (Fig. 4D–4G). In F-FMT recipients, the phylum Proteobacteria and gram-negative genus Prevotella were enriched before FMT and after FMT (days 0 and 7) (Figs. 4D and 4E). By contrast, FD-FMT significantly changed the microbiota. Specifically, 16 microbial taxa were detected at different concentrations before FMT (day 0), which may represent the cause of pathogenesis; at the highest score obtained from Proteobacteria, a major phylum of gram-negative bacteria (Figs. 4F and 4G). After FMT (day 7), a consortium of 19 microbial taxa, including gram-positive and -negative bacteria, such as those belonging to the phylum Bacteroidetes, and several genera of Prevotella, Blautia, and Selenomonas, were enriched in FD-FMT (Figs. 4E and 4G). The alpha diversity of the microbiota in the recipients’ calves increased close to the healthy conditions when FD-FMT (but not F-FMT) was applied (Fig. 4H). In addition, the beta diversity of the microbiota showed a lack of difference in baseline in the recipient calves used for F-FMT and FD-FMT before treatment. Interestingly, significant differences were observed when comparing the donors and recipients calves before FMT in both treatments (Fig. 4I and Supplementary Table 3). Moreover, no differences were observed between the donors and recipients’ calves 7 days after FD-FMT but not after F-FMT, confirming the possibility that the recipient calves can acquire the microbial composition of the donor within 7 days via FD-FMT (Supplementary Table 4). Further analysis of the extent of bacterial biogenesis after FMT using SourceTracker13 showed that the average contribution of the donors to the development of microbiota of the recipient calves on day 1 and 7 was 5% and 11% for F-FMT; respectively, whereas it reached 10% and 40%; respectively, for FD-FMT (Fig. 4J). Therefore, in the context of the donor’s microbiota engraftment, freeze-drying may be a superior method of microbial conditioning for FMT compared to just freezing.
Efficacy of F-FMT and FD-FMT in altering the intestinal metabolites in the recipients
To determine the exact effects of FD-FMT and F-FMT in treating diarrhea in the recipient calves, the fecal metabolites of the donors and the recipients before and 7 days after FMT were comprehensively analyzed using the capillary electrophoresis–Time-of-Flight Mass Spectrometry (CETOFMS). An interactive heatmap was provided to demonstrate all metabolites identified using Metaboanalyst14 (Fig. 5A). Principal component analysis (PCA) revealed significant differences in pretreatment (day 0) and posttreatment (day 7) recipients compared with the donors for F-FMT (Fig. 5B and Supplementary Table 5) but only between the donors and pretransplant (day 0) recipients and not between donors and posttransplant (day 7) recipients for FD-FMT (Fig. 5C and Supplementary Table 6). Procrustes analysis assessed similarity of metabolites between donors and recipients before (day 0) and after (day 7) FMT. In F-FMT, the Procrustes correlation was low (R2 = 0.2467) in donor vs. recipient on day 0 compared with donor vs. recipient on day 7 (R2 = 0.3748) (Supplementary Fig. 5A and 5B). However, in FD-FMT, the Procrustes correlation between donors and recipients was higher after FMT on day 7 than day 0 (R2 = 0.2036 vs. R2 = 0.4334) (Supplementary Fig. 5C and 5D). Because overall Procrustes correlation between donors and recipients was higher in FD-FMT than F-FMT (Supplementary Figs. 5A–5D), recipient metabolites were potentially same as donor metabolites in FD-FMT. Notably, that the number of metabolites were reduced but not increased by FMT, which were greater for FD-FMT than F-FMT. Specifically, within 7 days after F-FMT or FD-FMT, 43 or 29 metabolites were upregulated and 35 or 62 were downregulated, respectively (Supplementary Tables 7 and 8). To identify the metabolic pathways that differed between the F-FMT and FD-FMT, KEGG enrichment analysis was performed based on total number of metabolites obtained from the fold-change (FC) analysis. KEGG enrichment analysis revealed that 38 metabolic pathways were involved in F-FMT and FD-FMT (Supplementary Fig. 6A and 6B and Supplementary Table 9).
For F-FMT, the metabolites were mainly involved in purine metabolism, arginine biosynthesis, and pyrimidine metabolism (P < 0.05) (Supplementary Table 9). For FD-FMT, the metabolites were primarily involved in glycine, serine, and threonine metabolism; aminoacyl-tRNA biosynthesis; histidine metabolism; glycerophospholipid metabolism; arginine and proline metabolism; alanine, aspartate and glutamate metabolism, and glutathione metabolism (P < 0.05) (Supplementary Table 10). This study confirmed that the major metabolites in FD-FMT were involved in amino acid metabolism pathways. Next, the major metabolites detected in ATP-binding cassette (ABC) transporter, amino acid metabolism, lipid fatty acids, energy, polyamines, methylated compounds, and vitamins, short-chain fatty acids and bile acids were profiled in this study. Meanwhile, the levels of metabolites belonging to the ABC transporter metabolism, such as amino acids, were highly affected by FD-FMT rather than by F-FMT (Figs. 5D–5F and Supplementary Figs. 7–10). During the diarrheal condition, high levels of several amino acids, such as arginine, proline, and histidine were abundantly present in the feces, which may represent a microbial dysbiosis6. Thus, FD-FMT displayed a clear understanding to retain the microbial symbiosis condition by enhancing amino acid utilizing bacteria in the gut. Furthermore, by using variable importance in projection (VIP) scores obtained from the partial least-squares discriminant analysis (PLS-DA), top 15 discriminating metabolites were found, in where in FD-FMT cases amino acid histidine, phenylalanine, cysteine, leucine, valine and isoleucine were found downregulated (Figs. 5G and 5H). Consequently, these results indicated that FD-FMT was more efficient than F-FMT in changing the metabolic environment in feces of the recipient calves within 7 days after treatment.
Establishing the microbiota–metabolite correlation by FD-FMT during disease recovery
To investigate whether the microbiota and metabolites were closely associated during disease recovery post F-FMT or FD-FMT, Pearson’s correlation analysis was conducted. Specifically, the levels of the 16 selected microbiotas categorized as residents or colonizers using Microbiome Multivariable Association with Linear Models (MaAsLin2)15 in the Microbiomeanalyst16 platform (Fig. 6A and Supplementary Table 11). The 15 key metabolites shown in Figs. 5G and 5H were used for correlation analysis. Of the 240 (16 × 15) combinations, F-FMT showed five positive correlations before FMT (day 0) (Fig. 6B) and three positive and two negative correlations 7 days after FMT (day 7) (Fig. 6C). By contrast, of the 224 (16 × 14) combinations, FD-FMT displayed 14 positive and six negative correlations on day 0 (Fig. 6D) and 17 positive and two negative correlations on day 7 (Fig. 6E). For example, Campylobacteraceae and Enterobacteriaceae were positively correlated with choline (P < 0.001 and P < 0.01), N6-methyllysine (P < 0.01 and P < 0.05), and succinic acid (both P < 0.001).
Given that FD-FMT had a high number of correlations between metabolites and microbial taxa, these results indicate that compared with F-FMT, FD-FMT significantly changed the intestinal environment. To compare the predictive potential of F-FMT and FD-FMT in disease recovery through changes in microbiota and metabolites 7 days after FMT (day 7) vs. before FMT (day 0), a prediction model with area under the curve (AUC) of the receiver operating characteristic (ROC) was used17, using three feature sets: microbial taxa (n = 240), fecal metabolites (n = 596), and a combination of microbial taxa and fecal metabolites (n = 836). The AUC values of microbiota, metabolites, and microbiota + metabolites were 0.91, 0.78, and 0.76, respectively for F-FMT (Supplementary Fig. 11A), and 0.94, 0.81, and 0.81, respectively, for FD-FMT (Supplementary Fig. 11B). Thus, FD-FMT exhibited excellent discriminative power and superior performance in altering the microbiota and metabolites 7 days after FMT for diarrhea remission in calves (Fig. 6F). Overall, despite the high functional similarity between F-FMT and FD-FMT, FD-FMT exhibited broader spectral performance and was superior to F-FMT in promoting changes in the microflora and metabolites.