Skip to main content

Trace amounts of antibiotic altered metabolomic and microbial profiles of weaned pigs infected with a pathogenic E. coli

Abstract

Background

Our previous study has shown that supplementation of trace amounts of antibiotic exacerbated the detrimental effects of enterotoxigenic E. coli (ETEC) infection and delayed the recovery of pigs that may be associated with modified metabolites and metabolic pathways. Therefore, the objective of this study was to explore the impacts of trace levels of antibiotic (carbadox) on host metabolic profiles and colon microbiota of weaned pigs experimentally infected with ETEC F18.

Results

The multivariate analysis highlighted a distinct metabolomic profile of serum and colon digesta between trace amounts of antibiotic (TRA; 0.5 mg/kg carbadox) and label-recommended dose antibiotic (REC; 50 mg/kg carbadox) on d 5 post-inoculation (PI). The relative abundance of metabolomic markers of amino acids, carbohydrates, and purine metabolism were significantly differentiated between the TRA and REC groups (q < 0.2). In addition, pigs in REC group had the highest (P < 0.05) relative abundance of Lactobacillaceae and tended to have increased (P < 0.10) relative abundance of Lachnospiraceae in the colon digesta on d 5 PI. On d 11 PI, pigs in REC had greater (P < 0.05) relative abundance of Clostridiaceae compared with other groups, whereas had reduced (P < 0.05) relative abundance of Prevotellaceae than pigs in control group.

Conclusions

Trace amounts of antibiotic resulted in differential metabolites and metabolic pathways that may be associated with its slow responses against ETEC F18 infection. The altered gut microbiota profiles by label-recommended dose antibiotic may contribute to the promotion of disease resistance in weaned pigs.

Background

Trace amounts of antibiotics are emerging contaminants of environmental concern due to their potential risks on non-target organisms and the spread of bacterial resistance [1]. The excessive and imprudent use of human and veterinary antibiotics significantly contributes to a continuous release of antibiotics into the environment, thus, a variety of trace concentrations of antibiotics have been detected in surface water, wastewater, soil, air, and dust [2,3,4,5]. Exposure to trace levels of antibiotics may cause unexpected adverse effects on the host, such as allergic reactions, disruption of digestive system function, and shaping the metabolites and microbial community [6, 7]. In particular, exposure to trace amounts of antibiotics at early life can result in the alteration of microbiota and metabolic networks, which further accelerate the possible development of resistant strains [8, 9]. Previous studies in mice have shown that subtherapeutic concentrations of penicillin exposure at early life interfered with the development of the intestinal immune system [10] and induced metabolic changes due to altered intestinal microbiota [11].

Our recent study has shown that trace amounts of antibiotic in feed exacerbated the detrimental effects of enterotoxigenic E. coli (ETEC) infection by increasing diarrhea and systemic inflammation in weanling pigs [12]. The exact mechanisms for the exacerbation of ETEC infection by trace amounts of antibiotic has not been determined, but it has been suggested that trace concentrations of antibiotics can act as signaling molecules to trigger specific bacterial responses [13]. Therefore, the objective of this study was to explore the impacts of feeding trace levels of antibiotic on host metabolic profiles and colon microbiota of weaned pigs experimentally infected with ETEC F18.

Materials and methods

Animals, housing, experimental design, and diet

The protocol for this study was reviewed and approved by the Institutional Animal Care and Use Committee at the University of California, Davis (IACUC #19322). Samples analyzed here were obtained as described in a previously published study by Kim et al. [12]. Briefly, 34 weanling pigs (crossbred; initial BW: 6.88 ± 1.03 kg) with an equal number of gilts and barrows were randomly assigned to one of three treatments in a randomized complete block design with body weight within sex and litter as the blocks and pig as the experimental unit. The 3 dietary treatments included: 1) the complex nursery basal diet (control; CON), 2) addition of 0.5 mg/kg carbadox (trace amounts of antibiotic; TRA) to the basal diet, or 3) addition of 50 mg/kg carbadox (label-recommended dose of antibiotic; REC) to the basal diet. All diets were formulated to meet pig nutritional requirements and provided as mash form throughout the experiment.

After 7 d of adaptation, all pigs were orally inoculated with 3 mL of ETEC F18 using a disposable Luer-lock syringe for 3 consecutive days from d 0 post-inoculation (PI). The ETEC F18 was originally isolated from a field disease outbreak by the University of Illinois Veterinary Diagnostic Lab (isolate number: U.IL-VDL # 05–27,242) and expresses heat-labile toxin (LT), heat-stable toxin b (STb), and Shiga-like toxin (SLT-2). The ETEC F18 inoculums were provided at 1010 colony-forming unit (CFU) per 3 mL dose in phosphate-buffered saline (PBS). This dose caused mild diarrhea in the current study as reported in Kim et al. [12], which is consistent with our previous published researches [14,15,16]. Growth performance, blood profiles, and immune responses were also reported in previous work [12].

Sample collections

Sixteen pigs (6 pigs in CON, 4 pigs in TRA, and 6 pigs in REC) were randomly selected and euthanized on d 5 PI near the peak of infection, and the remaining pigs (6 pigs in CON, 5 pigs in TRA, and 7 pigs in REC) were euthanized at the end of the experiment (d 11 PI) that was the recovery period of the infection. The selection of necropsy time was based on the results of clinical observations and immune response parameters that were reported in previously published research using the same ETEC strain and inoculation dose [15, 16]. Before euthanasia, pigs were anesthetized with a 1-mL mixture of 100 mg telazol, 50 mg ketamine, and 50 mg xylazine (2:1:1) by intramuscular injection. After anesthesia, intracardiac injection with 78 mg sodium pentobarbital (Vortech Pharmaceuticals, Ltd., Dearborn, MI, USA) per 1 kg of BW was used to euthanize each pig. Blood samples were collected from the jugular vein of all pigs without EDTA to yield serum before ETEC challenge (d 0), and on d 5, and 11 PI. Serum samples were collected by centrifuging approximately 5 mL of whole blood samples at 20 °C at 1500 × g for 15 min and immediately stored at − 80 °C until untargeted metabolomics analysis. Colon digesta were collected from the distal colon of pigs on d 5 and 11 PI and immediately snap-frozen in liquid nitrogen and stored at − 80 °C for untargeted metabolomics and microbiome analysis.

Untargeted metabolomics analysis

The untargeted metabolomics analysis was performed by the NIH West Coast Metabolomics Center using gas chromatography (Agilent 6890 gas chromatograph controlled using Leco ChromaTOF software version 2.32, Agilent, Santa Clara, CA, USA) coupled with time-of-flight mass spectrometry (GC/TOF-MS) (Leco Pegasus IV time-of-flight mass spectrometer controlled using Leco ChromaTOF software version3 2.32, Leco, Joseph, MI, USA). Metabolite extraction was performed following procedures described previously [17]. Briefly, frozen colon digesta samples (approximately 10 mg) were homogenized using a Retsch ball mill (Retsch, Newtown, PA, USA) for 30 s at 25 times/s. After homogenization, a prechilled (− 20 °C) extraction solution (isopropanol/acetonitrile/water at the volume ratio 3:3:2, degassed with liquid nitrogen) was added at a volume of 1 mL/20 mg of sample. Samples were then vortexed and shaken for metabolite extraction. After centrifugation at 12,800 × g for 2 min, the supernatant was collected and divided into two equal aliquots and concentrated at room temperature for 4 h in a cold-trap vacuum concentrator (Labconco Centrivap, Kansas City, MO, USA). To separate complex lipids and waxes, the residue was re-suspended in 500 μL of 50% aqueous acetonitrile and centrifuged at 12,800 × g for 2 min. The resultant supernatant was collected and concentrated in the vacuum concentrator. Dried sample extracts were derivatized and mixed with internal retention index markers (fatty acid methyl esters with the chain length of C8 to C30). The samples were injected for GC/TOF analysis, and all samples were analyzed in a single batch. Data acquisition by mass spectrometry and mass calibration using FC43 (perfluorotributylamine) before starting analysis sequences. Metabolite identifications was performed based on the two parameters: 1) Retention index window ± 2000 U (around ± 2 sec retention time deviation), and 2) Mass spectral similarity plus additional confidence criteria as detailed below (Data analysis). A detailed methodology for data acquisition and metabolite identification described in a previously published article by Fiehn et al. [17].

Gut microbiota in distal Colon

Bacterial DNA was extracted from digesta samples using the Quick-DNA Fecal/Soil Microbe Kit (Zymo Research, Irvine, CA, USA) following the manufacturer’s instructions. Extracted bacterial DNA was amplified with PCR, targeting the V4 region of the 16S rRNA gene with primers 515 F (5′- XXXXXXXXGTGTGCCAGCMGCCGCGGTAA-3′) with an 8 bp barcode (X) and Illumina adapter (GT) and 806 R (5′-GGACTACHVGGGTWTCTAAT-3′) [18]. Amplification included thermocycling conditions of 94 °C for 3 min for denaturation, 35 cycles of 94 °C for 45 s, 50 °C for 1 min, 72 °C for 1.5 min, and 72 °C for 10 min (final elongation). To reduce PCR bias, each sample was amplified in triplicate. Each PCR reaction included 2 μL of template DNA, 0.5 μL (10 μmol/L) of barcoded forward primer, 0.5 μL (10 μmol/L) of reverse primer, 12.5 μL of GoTaq 2X Green Master Mix (Promega, Madison, WI, USA), and 9.5 μL of nuclease-free water. The triplicate PCR products were pooled and subjectively quantified based on the brightness of the bands on a 2% agarose gel with SYBR safe (Invitrogen Co., Carlsbad, CA, USA). All amplicons were then pooled at equal amounts and further purified using the QIAquick PCR Purification Kit (Qiagen, Hilden, Germany). The purified library was submitted to the UC Davis Genome Center DNA Technologies Core for 250 bp paired-end sequencing on the Illumina MiSeq platform (Illumina, Inc. San Diego, CA, USA).

The software sabre (https://github.com/najoshi/sabre) was used to demultiplex and remove barcodes from raw sequences. Sequences were then imported into Quantitative Insights Into Microbial Ecology 2 (QIIME2; version 2018.6) for downstream filtering and bioinformatics analysis [19, 20]. Plugin q2-dada2 [21] was used for quality control and constructing features. Taxonomic classification was assigned using the feature-classifier plugin trained with the SILVA rRNA database 99% operational taxonomic units (OTU), version 132 [22, 23].

Data analysis

The metabolomics data were analyzed using different modules of a web-based platform, MetaboAnalyst 5.0 (https://www.metaboanalyst.ca) [24]. Data were filtered for peaks with detection rates less than 30% of missing abundances and normalized using logarithmic transformation and auto-scaling. Mass univariate analysis was performed using one-way ANOVA followed by Tukey’s post hoc test (adjusted P ≤ 0.05). Fold change analysis and t-tests were also conducted to determine the fold change and significance of each identified metabolite. Statistical significance was declared at a false discovery rate (FDR, Benjamini and Hochberg correction; q) q < 0.2 and fold change > 2.0. Partial least squares discriminant analysis (PLS-DA) was carried out to further identify discriminative variables (metabolites) among the treatment groups. Pathway analysis and metabolite set enrichment analysis were performed on identified metabolites that had a Variable Importance in Projection (VIP) score > 1.

Data visualization and statistical analysis for colon microbiota were conducted using R (version 3.6.1). Two alpha diversity indices, Chao1 and Shannon, were calculated using the phyloseq package. Relative abundance was calculated using the phyloseq package and visualized using the ggplot2 package in R. Relative abundance data were aggregated at various taxonomical levels. Shapiro-Wilk normality test and Bartlett test were used to verify normality and constant variance, respectively, in alpha diversity and relative abundance. Shannon index was analyzed using ANOVA with the statistical model, including sample collection days within treatment as fixed effects. Significance in Chao1 index and relative abundance was observed using Kruskal-Wallis rank-sum test followed by a Conover test for multiple pairwise comparisons using the agricolae package. Beta diversity was calculated based on the Bray-Curtis dissimilarity for principal coordinates analysis (PCoA). The homogeneity of multivariate dispersions was tested by the vegan package using the betadisper function, before the adonis function was used to calculate PERMANOVA with 999 replicate permutations. Statistical significance and tendency were considered at P < 0.05 and 0.05 ≤ P < 0.10, respectively.

Results

Metabolite profiles in serum

A total of 354 (134 identified and 220 unidentified) metabolites were detected in serum samples. Based on VIP scores and relative abundance, 3 metabolites (fructose, mannonic acid, and propyleneglycol) were up-regulated by TRA, compared with the pigs in REC on d 0 before ETEC challenge (Table 1). Supplementation of REC changed the abundances of 6 metabolites (2 up-regulated and 4 down-regulated) compared with CON, while REC changed 16 metabolites (6 up-regulated and 10 down-regulated) in comparison with TRA on d 5 PI. On d 11 PI, chenodeoxycholic acid was enriched, while glycerol and inositol-4-monophosphate were reduced in the CON group compared with REC. Pigs in TRA had greater chenodeoxycholic acid than pigs in REC, but 5 metabolites (p-tolyl glucuronide, glycerol, mannitol, 2-ketoisocaproic acid, and inositol-4-monophosphate) were decreased in pigs supplemented with TRA compared with pigs in REC. No differential metabolites were identified when comparing CON vs. TRA throughout the experiment. Based on the identified metabolites, a PLS-DA score plot with 95% confidence ranges (shaded areas) showed a clear separation between the TRA and REC groups throughout the experiment (Fig. 1). To further explore the metabolic profile differences among two dietary treatments, PLS-DA was performed for the following comparisons: (1) TRA vs. REC, and (2) CON vs. REC on d 0 before ETEC challenge, d 5 and d 11 PI. The score plot again distinguished the TRA from the REC, and also revealed the metabolic profile differences between CON and REC (Fig. S1).

Table 1 Serum metabolites that differed among the dietary treatment groups1
Fig. 1
figure 1

Partial Least Squares Discriminant Analysis (PLS-DA) 2D score plot of the metabolites in serum showed separated clusters between the TRA and REC groups on d 0 before ETEC challenge (A), d 5 (B) and 11 (C) PI.  = CON (Control);  = TRA (Trace amounts of antibiotic);  = REC (Label-recommended dose of antibiotic). Shaded areas in different colors represent in 95% confidence interval

Pathway analysis and metabolite set enrichment analysis were performed on metabolites in serum with VIP > 1. On d 0 before ETEC challenge, inositol phosphate metabolism, glyoxylate and discarboxylate metabolism, glycine, serine and threonine metabolism, citrate cycle, and ascorbate and aldarate metabolism were the most affected metabolic pathways when comparing CON with REC (Fig. S2A, C). Citrate (TCA) cycle, arginine biosynthesis, and alanine, aspartate, and glutamate metabolism were the most affected metabolic pathways when TRA was compared with REC (Fig. S2B, D). On d 5 PI, aminoacyl-tRNA biosynthesis, glycine, serine, and threonine metabolism, and phenylalanine, tyrosine, and tryptophan biosynthesis were the most affected metabolic pathways when comparing CON vs. REC (Fig. 2 A, C), while aminoacyl-tRNA biosynthesis, alanine, aspartate, and glutamate metabolism, and glycolysis and gluconeogenesis were the most affected metabolic pathways in a comparison of TRA vs. REC (Fig. 2B, D). On d 11 PI, arginine biosynthesis, alanine, aspartate and glutamate metabolism, D-glutamine and D-glutamate metabolism, pyrimidine metabolism, and citrate cycle were the most affected metabolic pathways when comparing CON with REC (Fig. S3A, C). Arginine biosynthesis, aminoacyl-tRNA biosynthesis, alanine, aspartate, and glutamate metabolism, and D-glutamine and D-glutamate metabolism were the most affected metabolic pathways in a comparison of TRA with REC (Fig. S3 B, D).

Fig. 2
figure 2

Significantly changed pathways in serum between the control (CON) and label-recommended dose of antibiotic (REC) groups (A), and trace amounts of antibiotic (TRA) and REC groups (B) on d 5 post-inoculation. The x-axis represents the pathway impact values and the y-axis represents the -log(P) values from the pathway enrichment analysis. Metabolite set enrichment analysis (C, D) shows the metabolic pathways were enriched in CON compared to REC, and TRA compared to REC on d 5 post-inoculation, respectively. Both pathway analysis and metabolite set enrichment analysis were performed using identified metabolites with VIP > 1

Metabolite profiles in distal colon digesta

A total of 398 (178 identified and 220 unidentified) metabolites were detected in colon digesta. Based on VIP score and relative abundance, 12 metabolites (9 up-regulated and 3 down-regulated) were differentiated on d 5 PI, and one metabolite, inosine, was decreased on d 11 PI in pigs fed with TRA when compared with pigs in the REC group (Table 2). No differential metabolites were identified when comparing CON vs. TRA, and CON vs. REC at d 5 and 11 PI. Based on the identified metabolites, a PLS-DA score plot with 95% confidence ranges (shaded areas) showed a clear separation between the TRA and REC groups at both PI time points (Fig. 3). The PLS-DA score plots in a pairwise manner also clearly separated TRA from REC on d 5 and 11 PI (Fig. S4).

Table 2 Colon digesta metabolites that differed among the dietary treatment groups1
Fig. 3
figure 3

Partial Least Squares Discriminant Analysis (PLS-DA) 2D score plot of the metabolites in colon digesta showed separated clusters between the TRA and REC groups on d 5 (A) and 11 (B) post-inoculation.  = CON (Control);  = TRA (Trace amounts of antibiotic);  = REC (Label-recommended dose of antibiotic). Shaded areas in different colors represent in 95% confidence interval

Pathway analysis and metabolite set enrichment analysis were performed on metabolites in colon digesta with VIP > 1. Starch and sucrose metabolism, purine metabolism, arginine biosynthesis, and arginine and proline metabolism were the most affected metabolic pathways when TRA group was compared with the REC group on d 5 PI (Fig. 4 A, C). Aminoacyl-tRNA biosynthesis, arginine biosynthesis, pentose and glucuronate interconversions, arginine and proline metabolism, alanine, aspartate, and glutamate metabolism, glutathione metabolism, and glyoxylate and dicarboxylate metabolism were the most affected metabolic pathways on d 11 PI when TRA group was compared with the REC group (Fig. 4 B, D).

Fig. 4
figure 4

Significantly changed pathways in colon digesta between the trace amounts of antibiotic (TRA) and label-recommended dose of antibiotic (REC) groups on d 5 (A) and 11 post-inoculation (B). The x-axis represents the pathway impact values and the y-axis represents the -log(P) values from the pathway enrichment analysis. Metabolite set enrichment analysis (C, D) shows the metabolic pathways were enriched in TRA compared to REC group on d 5 and 11 post-inoculation, respectively. Both pathway analysis and metabolite set enrichment analysis were performed using identified metabolites with VIP > 1

Microbial profiles in distal colon digesta

A total of 481,102 qualified reads were obtained with a mean of 15,034 reads per sample. A total of 3561 OTUs were identified in the current experiment. No differences were observed in the alpha diversity of distal colon content among dietary treatments on d 5 and d 11 PI. Both Chao1 and Shannon indices of the distal colon content were lower (P < 0.05) on d 11 PI than d 5 PI for pigs in the CON group (Fig. S5). Beta diversity (Adonis analysis based on the Bray-Curtis distance) indicated that day (days PI) was a significant factor associated with composition distance (R2 = 0.11, P < 0.05; Fig. S6). Compositional differences of the distal colon microbiota were also observed between CON vs. REC and TRA vs. REC on d 5 and d 11 PI (Pairwise-Adonis, P < 0.05; Fig. S6).

The dominant phyla in distal colon content were Firmicutes, Bacteroidetes, Proteobacteria, and Actinobacteria, regardless of treatment or sampling day (Fig. S7). Pigs in the TRA or REC group had a lower (P < 0.05) relative abundance of Actinobacteria than pigs in the CON group on d 11 PI. Within the Firmicutes phylum (Fig. 5), pigs in the TRA group had lower (P < 0.05) relative abundance of Lactobacillaceae (8.91% vs. 21.33%) than pigs in REC on d 5 PI, whereas REC had lower (P < 0.05) relative abundance of Lactobacillaceae (5.82% vs. 23.90% or 27.69%) than pigs in the CON or TRA groups on d 11 PI. Pigs in the REC group had higher (P < 0.05) relative abundance of Clostridiaceae (17.14% vs. 1.45%) and Streptococcaceae (10.09% vs. 0.21%), but lower (P < 0.05) relative abundance of Lachnospiraceae (20.25% vs. 27.44%) in the distal colon on d 11 PI than on d 5 PI. Within the Bacteroidetes phylum (Fig. 6), pigs in the TRA group had reduced (P < 0.05) relative abundance of Muribaculaceae (0.60% vs. 2.46%) and Rikenellaceae (0.61% vs. 3.09%) in distal colon on d 11 PI than on d 5 PI. On d 11 PI, pigs in the CON group had higher (P < 0.05) relative abundance of Prevotellaceae (13.78% vs. 9.32%) in distal colon content, compared with pigs in the REC group.

Fig. 5
figure 5

Stacked bar plot showing the relative abundance of Firmicutes family in colon digesta of enterotoxigenic E. coli F18 challenged pigs fed diets supplemented with different dose of antibiotic on d 5 and 11 post-inoculation (A). Violin plot showing the relative abundance of individual bacterial phylum (B). a-cMeans without a common superscript are different across both time points (Diet × day, P < 0.05). Each least squares mean represents 4 to 7 observations. CON = Control; TRA = Trace amount of antibiotic; REC = Label-recommended dose of antibiotic

Fig. 6
figure 6

Stacked bar plot showing the relative abundance of Bacteroidetes family in colon digesta of enterotoxigenic E. coli F18 challenged pigs fed diets supplemented with different dose of antibiotic on d 5 and 11 post-inoculation (A). Violin plot showing the relative abundance of individual bacterial phylum (B). a-bMeans without a common superscript are different across both time points (Diet × day, P < 0.05). Each least squares mean represents 4 to 7 observations. CON = Control; TRA = Trace amount of antibiotic; REC = Label-recommended dose of antibiotic

Discussion

In-feed antibiotics can influence nutrient metabolism and many biological processes in pigs by altering microbiota and metabolites [25, 26]. Antimicrobial agent used in present study, carbadox, is one of the most common antibiotics and widely used in the U.S. swine industry to control enteric diseases and to promote the growth of nursery pigs [27]. However, little is known about the impacts of trace amounts of antibiotics on metabolic and microbial changes in piglets, especially under disease-challenged conditions. The present study investigated the alteration of metabolic pathways in the serum and colon digesta by using an untargeted metabolomics approach when pigs were supplemented with different levels of the antibiotic carbadox. Results from the current study highlighted that supplementing label-recommended doses of antibiotics altered metabolomic markers related to nutrient metabolism in the serum and colon digesta. Moreover, supplementation of different levels of antibiotic modified microbial community composition and diversity to a different extent in the colon digesta of pigs challenged with ETEC F18. Our previous research reported that supplementing the label-recommended dose of antibiotic enhanced disease resistance and promoted growth, whereas supplementing trace amounts of antibiotic exacerbated the detrimental effects of ETEC F18 infection on performance and diarrhea, and systemic inflammation of weaned pigs [12]. Results from the current study will help us to understand the negative impacts of trace amounts of antibiotic on performance and health of young pigs by focusing on the gut microbiome and their metabolites and the host metabolism.

The metabolomics approach exploits high-throughput analytical measurements to identify host and microbiota metabolites and associated biological changes that are affected by internal or external factors to maintain homeostasis [28]. In the present study, differences in the metabolic profiles of serum and colon digesta were found predominately between pigs supplemented with trace amounts of antibiotic and label-recommended dose of antibiotic, especially during the peak infection period (d 5 PI). These findings suggest the comparative dose-response metabolic effects of antibiotics during ETEC infection in weaned pigs.

In-feed antibiotics mediate growth enhancement as a result of improved nutrient utilization in pigs. Growing evidence suggests that the administration of in-feed antibiotics can enhance nutrient digestibility and regulate the nutrient metabolism of the host [29]. The bacteriostatic activity of in-feed antibiotics may also impact the intestinal microbial metabolites by reducing growth depressing microbiota [30]. It was reported that in-feed antibiotics at a subtherapeutic concentration could enhance amino acid availability in piglets, as indicated by increased serum metabolomic markers that are associated with amino acid metabolism [25]. Amino acid metabolism is extremely important to support animal growth, maintain homeostasis, and regulate other biological processes in the host and intestinal microbiota [31, 32]. In the present study, metabolites related to amino acid metabolism (2-hydroxyvaleric acid, pipecolinic acid, histidine, and creatine) were enriched in the serum of pigs fed with the label-recommended antibiotic dose compared with pigs in the trace amounts of antibiotic group. This was likely due to the reduced peptide catabolism initiated by microbial protease activities when feeding label-recommended dose of antibiotic [33]. However, 2-hydroxyvaleric acid, a metabolomic marker of branched-chain amino acid catabolism, was observed to be reduced in the colon digesta of pigs fed with label-recommended dose antibiotic compared with the trace amounts of antibiotic group. These observations are in agreement with a previous research, in which Mu et al. [25] also reported that increased serum metabolites related to amino acid metabolism were concomitant with a decrease in jejunal metabolites associated with amino acid metabolism in pigs fed with a mixture of antibiotics at a growth-promoting dose. Thus, these results suggest that the systemic interplay between microbiota and metabolite profiles was promoted by feeding label-recommended dose of antibiotics. A previous study using metagenomic analysis also observed that antibiotics at subtherapeutic doses reduced the abundance of clusters of orthologous groups involved in protein metabolism in the fecal microbiota of pigs [34]. Consistent with performance data and clinical signs [12], it is not surprising to observe that trace amount of antibiotics had different impacts on serum and colon digesta metabolites that are associated with amino acid metabolism when compared with label-recommended dose of antibiotic. Previous in vitro research suggested that E. coli cells stimulated cellular functions and metabolic modifications of amino acid catabolism upon exposure to the antibiotic ampicillin at below the minimal inhibitory concentrations (sub-MIC) [35]. More specifically, E. coli cells treated with sub-MIC ampicillin resulted in increased amino acid depletion in Luria-Bertani (LB) media due to stress responses, which provided amino acids as a major energy source for cultured cells. This finding indicates that the alteration of metabolomic markers of amino acid metabolism caused by trace amounts of antibiotic in the current study may be related to the depletion of amino acids during the host response to ETEC infection. Subsequently, less amounts of amino acids might be available to support the growth of the pigs when they were challenged with ETEC and supplemented with trace amount of antibiotics.

Carbohydrate metabolism is essential to support the virulence of pathogenic enterobacteria [36]. It has been reported that the colonization of pathogenic E. coli in the mouse intestine was supported by the catabolism of several carbohydrates, including galactose, fucose, mannose, and maltose [37, 38]. In the present study, metabolomic markers related to galactose metabolism (glycerol and myo-inositol) and carbohydrate digestion and absorption (maltotriose) were enriched in the serum or colon digesta from pigs supplemented with trace amounts of antibiotic. This finding suggests that trace amounts of antibiotic may assist in constitution of an ecological niche for ETEC F18 colonization in the intestine of pigs, rather than exhibit its antibacterial activity. Besides the carbohydrate utilization by pathogens to colonize, carbohydrate metabolism is also vital for the systemic inflammatory response [39]. Baurhoo et al. [40] reported that a significant mobilization and catabolism of carbohydrates were observed in the intestine of chickens during LPS-induced systemic inflammation. In the present animal study, trace amounts of antibiotic exacerbated the intestinal and systemic inflammatory status of ETEC F18 challenged pigs [12]. Thus, the increased metabolites associated with carbohydrate metabolism in pigs supplemented with trace amount of antibiotics during the peak of ETEC infection indicates that these pigs may utilize more carbohydrates as energy sources to support their immune responses and recovery processes against ETEC F18 instead of growth.

Interestingly, supplementation of trace amounts of antibiotic also enriched serum metabolomic markers of purine metabolism (hypoxanthine and guanine) during the peak of ETEC infection. A previous in vitro study reported that Pasteurella multocida significantly increased the expression of proteins involved in purine synthesis and metabolism, in response to sub-MIC antibiotics, including amoxicillin, chlortetracycline, and enrofloxacin [41]. Ng et al. [42] also demonstrated that extremely low concentrations of antibiotics, such as tetracycline and macrolide, upregulated the expression of genes associated with purine metabolism in Streptococcus pneumoniae. The metabolites involved in purine metabolism are often upregulated in the activated immune cells as important immune signaling molecules [43]. For instance, previous research reported that mice infected with E. coli had enriched plasma metabolites that are linked to the purine metabolic pathway [44]. Likewise, growing evidence also suggests that trace concentrations of antibiotics may perform as signaling agents and trigger special bacterial responses, such as increased purine metabolism, following an infection [43, 45], Thus, our results indicate that purine metabolism might contribute to the elevated systemic inflammation in pigs fed with trace amounts of antibiotic [12].

The composition and diversity of gut microbial communities in pigs are greatly impacted by their age, health status, and nutrient components in feed [46,47,48]. To test the impacts of trace amounts of antibiotic on gut microbiota diversity, distal colon contents were collected, and 16S rRNA gene sequencing was performed. Consistent with previously published research, antibiotics-treated at recommended concentrations clustered separately from non-treated groups [34, 49], indicating that antibiotics administration at label-recommended dose altered colon microbiota composition and diversity. However, there was no clear separation in distal colon microbiota between pigs supplemented with trace amounts of antibiotic and pigs in the control group.

The three most abundant phyla found in the colon digesta of pigs in the present study were Firmicutes, Bacteroidetes, and Proteobacteria, which was consistent with previously published research [50, 51]. Within the Firmicutes and Bacteroidetes phyla, the relative abundance of Lachnospiraceae and Lactobacillaceae were enriched in the distal colon of pigs supplemented with label-recommended dose of antibiotic, but the relative abundance of Lactobacillaceae was reduced in colon digesta of trace amounts of antibiotic pigs during the peak infection period. Lachnospiraceae family contain numerous genera involved in producing butyric acid, which provides energy for other microbes and host epithelial cells and prevents the growth of other microbes in the digestive tract [52, 53]. Moreover, Lactobacillaceae were reported to be positively correlated with feed efficiency [54] and nitrogen, energy, cellulose, and hemicellulose digestibility in pigs [55]. Although the exact mechanism of antimicrobial effects is not yet clear, Lactobacillaceae are known for their health-promoting effects and for their ability to inhibit intestinal pathogens such as E. coli and Salmonella [56]. Thus, Lachnospiraceae and Lactobacillaceae have been proposed and investigated as biomarkers to predict the health status of pigs [57, 58]. Rhouma et al. [59] demonstrated that the ETEC F4 challenge suppressed the relative abundance of Lachnospiraceae and Lactobacillaceae in fecal contents of pigs, compared with unchallenged pigs. In addition, Dou et al. [60] also reported that diarrheic pigs, in natural post-weaning diarrhea, had a lower abundance of Lachnospiraceae and Lactobacillaceae in feces, compared with healthy pigs. Therefore, the modified intestinal microbial environment, including the enhanced presence of Lachnospiraceae and Lactobacillaceae, may contribute to the accelerated recovery from ETEC 18 infection in pigs supplemented with label-recommended dose of antibiotics.

Previous studies have also reported the contribution of intestinal microbiota to weight gain in pigs. For example, Kim et al. [61] observed that Clostridiaceae in the distal gut of pigs were positively correlated with weight gain, while Unno et al. [62] reported a negative correlation between weight and Prevotellaceae in feces when pigs were supplemented with different types of antibiotics. In the present study, pigs fed label-recommended dose antibiotic had increased relative abundance of Clostridiaceae but reduced relative abundance of Prevotellaceae compared to pigs in the control group on d 11 PI. These observations are consistent with the literature and confirmed the effectiveness of label-recommended dose of antibiotic for growth-promoting purposes [12].

In conclusion, the metabolomics and microbiome approaches in the present study identified the differential metabolites and their pathways in the serum and distal colon digesta of ETEC F18 challenged pigs. Compared with label-recommended dose of antibiotic, trace amounts of antibiotic oppositely affected metabolomic markers related to the metabolisms of amino acids, carbohydrates, and purine. Pigs administered label-recommended dose of antibiotic had marked modulation of microbial composition, which may be highly correlated with their enhanced growth performance and disease resistance in weaned pigs. The impacted metabolic pathways and colonic microbial shift may also be closely associated with the slow growth and delayed recovery from ETEC infection of weaned pigs supplemented with trace amounts of antibiotic. Future studies will consider incorporating targeted metabolomics and metagenomics to provide more insights into the potential risk of trace amounts of antibiotic on the host response to ETEC infection. The exploration of metabolomic markers and gut microbiome interaction will be important to decipher the mechanisms of how trace amounts of antibiotic negatively impact the health of young animals.

Availability of data and materials

All data generated or analyzed during this study are available from the corresponding author upon reasonable request.

Abbreviations

FDR:

False discovery rate

ETEC:

Enterotoxigenic E. coli

MIC:

Minimal inhibitory concentrations

OTU:

Operational taxonomic units

PI:

Post-inoculation

PLS-DA:

Partial least squares discriminant analysis

PCoA:

Principal coordinates analysis

VIP:

Variable importance in the projection

References

  1. Ben Y, Fu C, Hu M, Liu L, Wong MH, Zheng C. Human health risk assessment of antibiotic resistance associated with antibiotic residues in the environment: a review. Environ Res. 2019;169:483–93. https://doi.org/10.1016/j.envres.2018.11.040.

    Article  CAS  PubMed  Google Scholar 

  2. Jechalke S, Heuer H, Siemens J, Amelung W, Smalla K. Fate and effects of veterinary antibiotics in soil. Trends Microbiol. 2014;22(9):536–45. https://doi.org/10.1016/j.tim.2014.05.005.

    Article  CAS  PubMed  Google Scholar 

  3. Wei R, Ge F, Huang S, Chen M, Wang R. Occurrence of veterinary antibiotics in animal wastewater and surface water around farms in Jiangsu Province, China. Chemosphere. 2011;82(10):1408–14. https://doi.org/10.1016/j.chemosphere.2010.11.067.

    Article  CAS  PubMed  Google Scholar 

  4. Hamscher G, Hartung J. Veterinary antibiotics in dust: sources, environmental concentrations, and possible health hazards. In: Kümmerer K, editor. Pharm environ. Berlin: Springer; 2008. p. 95–102. https://doi.org/10.1007/978-3-540-74664-5_7.

    Chapter  Google Scholar 

  5. Manzetti S, Ghisi R. The environmental release and fate of antibiotics. Mar Pollut Bull. 2014;79(1-2):7–15. https://doi.org/10.1016/j.marpolbul.2014.01.005.

    Article  CAS  PubMed  Google Scholar 

  6. Sarmah AK, Meyer MT, Boxall ABA. A global perspective on the use, sales, exposure pathways, occurrence, fate and effects of veterinary antibiotics (VAs) in the environment. Chemosphere. 2006;65(5):725–59. https://doi.org/10.1016/j.chemosphere.2006.03.026.

    Article  CAS  PubMed  Google Scholar 

  7. Chen J, Ying GG, Deng WJ. Antibiotic residues in food: extraction, analysis, and human health concerns. J Agric Food Chem. 2019;67(27):7569–86. https://doi.org/10.1021/acs.jafc.9b01334.

    Article  CAS  PubMed  Google Scholar 

  8. Nisha AR. Antibiotic residues – a global health hazard. Vet World. 2008;1:375–7. https://doi.org/10.5455/vetworld.2008.375-377.

    Article  Google Scholar 

  9. Davies J, Spiegelman GB, Yim G. The world of subinhibitory antibiotic concentrations. Curr Opin Microbiol. 2006;9(5):445–53. https://doi.org/10.1016/j.mib.2006.08.006.

    Article  CAS  PubMed  Google Scholar 

  10. Jin S, Zhao D, Cai C, Song D, Shen J, Xu A, et al. Low-dose penicillin exposure in early life decreases Th17 and the susceptibility to DSS colitis in mice through gut microbiota modification. Sci Rep. 2017;7(1):43662. https://doi.org/10.1038/srep43662.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Cox LM, Yamanishi S, Sohn J, Alekseyenko AV, Leung JM, Cho I, et al. Altering the intestinal microbiota during a critical developmental window has lasting metabolic consequences. Cell. 2014;158(4):705–21. https://doi.org/10.1016/j.cell.2014.05.052.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Kim K, He Y, Jinno C, Kovanda L, Li X, Song M, et al. Trace amounts of antibiotic exacerbated diarrhea and systemic inflammation of weaned pigs infected with a pathogenic Escherichia coli. J Anim Sci. 2021;99(3):1–13. https://doi.org/10.1093/jas/skab073.

    Article  Google Scholar 

  13. Fajardo A, Martínez JL. Antibiotics as signals that trigger specific bacterial responses. Curr Opin Microbiol. 2008;11(2):161–7. https://doi.org/10.1016/j.mib.2008.02.006.

    Article  CAS  PubMed  Google Scholar 

  14. Liu Y, Song M, Che TM, Almeida JAS, Lee JJ, Bravo D, et al. Dietary plant extracts alleviate diarrhea and alter immune responses of weaned pigs experimentally infected with a pathogenic Escherichia coli. J Anim Sci. 2013;91(11):5294–306. https://doi.org/10.2527/jas.2012-6194.

    Article  CAS  PubMed  Google Scholar 

  15. Kim K, Ehrlich A, Perng V, Chase JA, Raybould H, Li X, et al. Algae-derived β-glucan enhanced gut health and immune responses of weaned pigs experimentally infected with a pathogenic E. coli. Anim Feed Sci Technol. 2019;248:114–25. https://doi.org/10.1016/j.anifeedsci.2018.12.004.

    Article  CAS  Google Scholar 

  16. Kim K, He Y, Xiong X, Ehrlich A, Li X, Raybould H, et al. Dietary supplementation of Bacillus subtilis influenced intestinal health of weaned pigs experimentally infected with a pathogenic E coli. J Anim Sci Biotechnol. 2019;10(1):52. https://doi.org/10.1186/s40104-019-0364-3.

    Article  CAS  Google Scholar 

  17. Fiehn O, Wohlgemuth G, Scholz M, Kind T, Lee DY, Lu Y, et al. Quality control for plant metabolomics: reporting MSI-compliant studies. Plant J. 2008;53(4):691–704. https://doi.org/10.1111/j.1365-313X.2007.03387.x.

    Article  CAS  PubMed  Google Scholar 

  18. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6(8):1621–4. https://doi.org/10.1038/ismej.2012.8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7. https://doi.org/10.1038/s41587-019-0209-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6. https://doi.org/10.1038/nmeth.f.303.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2012;41(D1):D590–6. https://doi.org/10.1093/nar/gks1219.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Bokulich NA, Kaehler BD, Rideout JR, Dillon M, Bolyen E, Knight R, et al. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome. 2018;6(1):1–17. https://doi.org/10.1186/s40168-018-0470-z.

    Article  Google Scholar 

  24. Pang Z, Chong J, Zhou G, de Lima Morais DA, Chang L, Barrette M, et al. MetaboAnalyst 5.0: narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 2021;49(W1):W388–96. https://doi.org/10.1093/nar/gkab382.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Mu C, Yang Y, Yu K, Yu M, Zhang C, Su Y, et al. Alteration of metabolomic markers of amino-acid metabolism in piglets with in-feed antibiotics. Amino Acids. 2017;49(4):771–81. https://doi.org/10.1007/s00726-017-2379-4.

    Article  CAS  PubMed  Google Scholar 

  26. Pi Y, Gao K, Peng Y, Mu CL, Zhu WY. Antibiotic-induced alterations of the gut microbiota and microbial fermentation in protein parallel the changes in host nitrogen metabolism of growing pigs. Animal. 2019;13(2):262–72. https://doi.org/10.1017/S1751731118001416.

    Article  CAS  PubMed  Google Scholar 

  27. USDA. 2015. Swine 2012; part II: reference of swine health and health management in the United States, 2012. USDA–APHIS–VS–CEAH–NAHMS. #676.0216.

  28. Weckwerth W. Metabolomics in systems biology. Annu Rev Plant Biol. 2003;54(1):669–89. https://doi.org/10.1146/annurev.arplant.54.031902.135014.

    Article  CAS  PubMed  Google Scholar 

  29. Gaskins HR, Collier CT, Anderson DB. Antibiotics as growth promotants: mode of action. Anim Biotechnol. 2002;13(1):29–42. https://doi.org/10.1081/ABIO-120005768.

    Article  CAS  PubMed  Google Scholar 

  30. Brüssow H. Growth promotion and gut microbiota: insights from antibiotic use. Environ Microbiol. 2015;17(7):2216–27. https://doi.org/10.1111/1462-2920.12786.

    Article  PubMed  Google Scholar 

  31. Christensen HN. Role of amino acid transport and countertransport in nutrition and metabolism. Physiol Rev. 1990;70(1):43–77. https://doi.org/10.1152/physrev.1990.70.1.43.

    Article  CAS  PubMed  Google Scholar 

  32. Wu G. Amino acids: metabolism, functions, and nutrition. Amino Acids. 2009;37(1):1–17. https://doi.org/10.1007/s00726-009-0269-0.

    Article  CAS  PubMed  Google Scholar 

  33. Romick-Rosendale LE, Goodpaster AM, Hanwright PJ, Patel NB, Wheeler ET, Chona DL, et al. NMR-based metabonomics analysis of mouse urine and fecal extracts following oral treatment with the broad-spectrum antibiotic enrofloxacin (Baytril). Magn Reson Chem. 2009;47(S1):S36–46. https://doi.org/10.1002/mrc.2511.

    Article  CAS  PubMed  Google Scholar 

  34. Looft T, Johnson TA, Allen HK, Bayles DO, Alt DP, Stedtfeld RD, et al. In-feed antibiotic effects on the swine intestinal microbiome. Proc Natl Acad Sci. 2012;109(5):1691–6. https://doi.org/10.1073/pnas.1120238109.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Mathieu A, Fleurier S, Frénoy A, Dairou J, Bredeche MF, Sanchez-Vizuete P, et al. Discovery and function of a general core hormetic stress response in E. coli induced by sublethal concentrations of antibiotics. Cell Rep. 2016;17(1):46–57. https://doi.org/10.1016/j.celrep.2016.09.001.

    Article  CAS  PubMed  Google Scholar 

  36. Le Bouguénec C, Schouler C. Sugar metabolism, an additional virulence factor in enterobacteria. Int J Med Microbiol. 2011;301(1):1–6. https://doi.org/10.1016/j.ijmm.2010.04.021.

    Article  CAS  PubMed  Google Scholar 

  37. Fabich AJ, Jones SA, Chowdhury FZ, Cernosek A, Anderson A, Smalley D, et al. Comparison of carbon nutrition for pathogenic and commensal Escherichia coli strains in the mouse intestine. Infect Immun. 2008;76(3):1143–52. https://doi.org/10.1128/IAI.01386-07.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Jones SA, Jorgensen M, Chowdhury FZ, Rodgers R, Hartline J, Leatham MP, et al. Glycogen and maltose utilization by Escherichia coli O157:H7 in the mouse intestine. Infect Immun. 2008;76(6):2531–40. https://doi.org/10.1128/IAI.00096-08.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Wolowczuk I, Verwaerde C, Viltart O, Delanoye A, Delacre M, Pot B, et al. Feeding our immune system: impact on metabolism. Clin Dev Immunol. 2008;2008:639803. https://doi.org/10.1155/2008/639803.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Baurhoo B, Ferket P, Ashwell CM, de Oliviera J, Zhao X. Cell walls of Saccharomyces cerevisiae differentially modulated innate immunity and glucose metabolism during late systemic inflammation. PLoS One. 2012;7(1):1–16. https://doi.org/10.1371/journal.pone.0030323.

    Article  CAS  Google Scholar 

  41. Nanduri B, Lawrence ML, Peddinti DS, Burgess SC. Effects of subminimum inhibitory concentrations of antibiotics on the Pasteurella multocida proteome: a systems approach. Comp Funct Genomics. 2008;2008:254836. https://doi.org/10.1155/2008/254836.

    Article  CAS  PubMed Central  Google Scholar 

  42. Ng WL, Kazmierczak KM, Robertson GT, Gilmour R, Winkler ME. Transcriptional regulation and signature patterns revealed by microarray analyses of Streptococcus pneumoniae R6 challenged with sublethal concentrations of translation inhibitors. Am Soc Microbiol. 2003;185(1):359–70. https://doi.org/10.1128/JB.185.1.359-370.2003.

    Article  CAS  Google Scholar 

  43. Cekic C, Linden J. Purinergic regulation of the immune system. Nat Rev Immunol. 2016;16(3):177–92. https://doi.org/10.1038/nri.2016.4.

    Article  CAS  PubMed  Google Scholar 

  44. Yang JH, Bhargava P, McCloskey D, Mao N, Palsson BO, Collins JJ. Antibiotic-induced changes to the host metabolic environment inhibit drug efficacy and alter immune function. Cell Host Microbe. 2017;22:757–65.e3. https://doi.org/10.1016/j.chom.2017.10.020.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Romero D, Traxler MF, López D, Kolter R. Antibiotics as signal molecules. Chem Rev. 2011;111(9):5492–505. https://doi.org/10.1021/cr2000509.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Gresse R, Chaucheyras-Durand F, Fleury MA, Van de Wiele T, Forano E, Blanquet-Diot S. Gut microbiota dysbiosis in postweaning piglets: understanding the keys to health. Trends Microbiol. 2017;25(10):851–73. https://doi.org/10.1016/j.tim.2017.05.004.

    Article  CAS  PubMed  Google Scholar 

  47. Frese SA, Parker K, Calvert CC, Mills DA. Diet shapes the gut microbiome of pigs during nursing and weaning. Microbiome. 2015;3(1):28. https://doi.org/10.1186/s40168-015-0091-8.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Guevarra RB, Hong SH, Cho JH, Kim BR, Shin J, Lee JH, et al. The dynamics of the piglet gut microbiome during the weaning transition in association with health and nutrition. J Anim Sci Biotechnol. 2018;9(1):54. https://doi.org/10.1186/s40104-018-0269-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Schokker D, Zhang J, Vastenhouw SA, Heilig HGHJ, Smidt H, Rebel JMJ, et al. Long-lasting effects of early-life antibiotic treatment and routine animal handling on gut microbiota composition and immune system in pigs. PLoS One. 2015;10(2):1–18. https://doi.org/10.1371/journal.pone.0116523.

    Article  CAS  Google Scholar 

  50. Crespo-Piazuelo D, Estellé J, Revilla M, Criado-Mesas L, Ramayo-Caldas Y, Óvilo C, et al. Characterization of bacterial microbiota compositions along the intestinal tract in pigs and their interactions and functions. Sci Rep. 2018;8(1):12727. https://doi.org/10.1038/s41598-018-30932-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. He Y, Jinno C, Kim K, Wu Z, Tan B, Li X, et al. Dietary Bacillus spp enhanced growth and disease resistance of weaned pigs by modulating intestinal microbiota and systemic immunity. J Anim Sci Biotechnol. 2020;11:101. https://doi.org/10.1186/s40104-020-00498-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Meehan CJ, Beiko RG. A Phylogenomic view of ecological specialization in the Lachnospiraceae, a family of digestive tract-associated bacteria. Genome Biol Evol. 2014;6(3):703–13. https://doi.org/10.1093/gbe/evu050.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Guarner F, Malagelada JR. Gut flora in health and disease. Lancet. 2003;361(9356):512–9. https://doi.org/10.1016/S0140-6736(03)12489-0.

    Article  PubMed  Google Scholar 

  54. Quan J, Cai G, Ye J, Yang M, Ding R, Wang X, et al. A global comparison of the microbiome compositions of three gut locations in commercial pigs with extreme feed conversion ratios. Sci Rep. 2018;8(1):4536–46. https://doi.org/10.1038/s41598-018-22692-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Le Sciellour M, Labussière E, Zemb O, Renaudeau D. Effect of dietary fiber content on nutrient digestibility and fecal microbiota composition in growing-finishing pigs. PLoS One. 2018;13(10):e0206159. https://doi.org/10.1371/journal.pone.0206159.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Stecher B, Hardt W-D. The role of microbiota in infectious disease. Trends Microbiol. 2008;16(3):107–14. https://doi.org/10.1016/j.tim.2007.12.008.

    Article  CAS  PubMed  Google Scholar 

  57. Guevarra RB, Lee JH, Lee SH, Seok M-J, Kim DW, Kang BN, et al. Piglet gut microbial shifts early in life: causes and effects. J Anim Sci Biotechnol. 2019;10(1):1. https://doi.org/10.1186/s40104-018-0308-3.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Arfken AM, Frey JF, Summers KL. Temporal dynamics of the gut bacteriome and mycobiome in the weanling pig. Microorganisms. 2020;8(6):868. https://doi.org/10.3390/microorganisms8060868.

    Article  CAS  PubMed Central  Google Scholar 

  59. Rhouma M, Braley C, Th W, Thibodeau A, Quessy S, Fravalo P. Evolution of pig fecal microbiota composition and diversity in response to enterotoxigenic Escherichia coli infection and colistin treatment in weaned piglets. Microorganisms. 2021;7(7):1459–74. https://doi.org/10.3390/microorganisms9071459.

    Article  CAS  Google Scholar 

  60. Dou S, Gadonna-Widehem P, Rome V, Hamoudi D, Rhazi L, Lakhal L, et al. Characterisation of early-life fecal microbiota in susceptible and healthy pigs to post-weaning diarrhoea. PLoS One. 2017;12(1):e0169851. https://doi.org/10.1371/journal.pone.0169851.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Kim HB, Borewicz K, White BA, Singer RS, Sreevatsan S, Tu ZJ, et al. Microbial shifts in the swine distal gut in response to the treatment with antimicrobial growth promoter, tylosin. Proc Natl Acad Sci. 2012;109(38):15485–90. https://doi.org/10.1073/pnas.1205147109.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Unno T, Kim J, Guevarra RB, Nguyen SG. Effects of antibiotic growth promoter and characterization of ecological succession in swine gut microbiota. J Microbiol Biotechnol. 2015;25(4):431–8. https://doi.org/10.4014/jmb.1408.08063.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

Not applicable.

Funding

This project was supported by the United States Department of Agriculture (USDA) National Institute of Food and Agriculture (NIFA), multistate projects W4002 and NC1202.

Author information

Authors and Affiliations

Authors

Contributions

The contributions of the authors were as follows: KK conducted the experiment and wrote the manuscript. CJ and PJ helped to conduct animal trial, part of the laboratory work, data analysis, and helped to revise the manuscript. YL was the principal investigator. YL designed the experiment, oversaw the development of the study and wrote the last version of the manuscript. The authors declare no conflicts of interest. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Yanhong Liu.

Ethics declarations

Ethics approval and consent to participate

The protocol for this study was reviewed and approved by the Institutional Animal Care and Use Committee at the University of California, Davis (IACAC #19322). The study was conducted at the Cole Facility at the University of California, Davis.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Supplementary Information

Additional file 1: Fig. S1.

Partial Least Squares Discriminant Analysis (PLS-DA) 2D score plot of the metabolites in serum revealed significant differences on d 0 before E. coli challenge, d 5 and 11 post-inoculation between the TRA and REC groups (A-C) and CON and REC (D-E), respectively.  = CON (Control);  = TRA (Trace amounts of antibiotic);  = REC (Label-recommended dose of antibiotic). Shaded areas in different colors represent in 95% confidence interval. Fig. S2. Significantly changed pathways in serum between the control (CON) and label-recommended dose of antibiotic (REC) groups (A), and trace amounts of antibiotic (TRA) and REC groups (B) on d 0 before E. coli challenge. The x-axis represents the pathway impact values and the y-axis represents the -log(P) values from the pathway enrichment analysis. Metabolite set enrichment analysis (C, D) shows the metabolic pathways were enriched in CON compared to REC, and TRA compared to REC on d 0 before E. coli challenge, respectively. Both pathway analysis and metabolite set enrichment analysis were performed using identified metabolites with VIP > 1. Fig. S3. Significantly changed pathways in serum between the control (CON) and label-recommended dose of antibiotic (REC) groups (A), and trace amounts of antibiotic (TRA) and REC groups (B) on d 11 post-inoculation. The x-axis represents the pathway impact values and the y-axis represents the -log(P) values from the pathway enrichment analysis. Metabolite set enrichment analysis (C, D) shows the metabolic pathways were enriched in CON compared to REC, and TRA compared to REC on d 11 post-inoculation, respectively. Both pathway analysis and metabolite set enrichment analysis were performed using identified metabolites with VIP > 1. Fig. S4. Partial Least Squares Discriminant Analysis (PLS-DA) 2D score plot of the metabolites in colon digesta revealed significant differences between the TRA and REC groups on d 5 (A) and 11 (B) post-inoculation.  = TRA (Trace amounts of antibiotic);  = REC (Label-recommended dose of antibiotic). Shaded areas in different colors represent in 95% confidence interval. Fig. S5. Alpha diversity as indicated by Chao 1 (A) and Shannon (B) indices in colon digesta of enterotoxigenic E. coli F18 challenged pigs fed diets supplemented with different dose of antibiotic on d 5 and 11 post-inoculation. a-bMeans without a common superscript are different across both time points (Diet × day, P < 0.05). Each least squares mean represents 4 to 7 observations. CON = Control; TRA = Trace amount of antibiotic; REC = Label-recommended dose of antibiotic. Fig. S6. Beta diversity of colon digesta of enterotoxigenic E. coli F18 challenged pigs fed diets supplemented with different dose of antibiotic on d 5 and 11 post-inoculation. Data were analyzed by principal coordinate analysis (PCoA) based on the Bray-Curtis dissimilarity. Symbols indicate dietary treatments and colors indicate different sampling dates. CON = Control; TRA = Trace amount of antibiotic; REC = Label-recommended dose of antibiotic. Fig. S7. Stacked bar plot showing the relative abundance of bacterial phyla in colon digesta of enterotoxigenic E. coli F18 challenged pigs fed diets supplemented with different dose of antibiotic on d 5 and 11 post-inoculation (A). Violin plot showing the relative abundance of individual bacterial phylum (B). a-cMeans without a common superscript are different across both time points (Diet × day, P < 0.05). Each least squares mean represents 4 to 7 observations. CON = Control; TRA = Trace amount of antibiotic; REC = Label-recommended dose of antibiotic.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kim, K., Jinno, C., Ji, P. et al. Trace amounts of antibiotic altered metabolomic and microbial profiles of weaned pigs infected with a pathogenic E. coli. J Animal Sci Biotechnol 13, 59 (2022). https://doi.org/10.1186/s40104-022-00703-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40104-022-00703-5

Keywords