- Research Article
- Open Access
- Published:
Enteroendocrine peptides, growth, and the microbiome during the porcine weaning transition
Animal Microbiome volume 4, Article number: 56 (2022)
Abstract
Background
Growth rate in pigs can be affected by numerous factors that also affect feeding behavior and the microbiome. Recent studies report some communication between the microbiome and the enteroendocrine system. The present study examined if changes in the piglet microbiome between birth and during the weaning transition can be correlated either positively or negatively with growth rate and plasma concentrations of enteroendocrine peptides.
Results
During the post-weaning transition, a 49% reduction in average daily gain was observed at day 24 (P < 0.05) relative to day 21. Pigs recovered by day 28 with body weight and average daily gain increases of 17% and 175%, respectively relative to day 24 and the highest rate of gain was measured at day 35 (462 g/day). The time interval between day 21–24 had the highest number of correlations (n = 25) between the relative abundance differences in taxa over time and corresponding percent weight gain. Amplicon sequence variants with the greatest correlation with percent weight gain between day 21–24 belonged to families Prevotellaceae NK3B31 (ρ = 0.65, P < 0.001), Veillonellaceae (ρ = 0.63, P < 0.001) and Rikenellaceae RC9 (ρ = 0.62, P < 0.001). Seven taxa were positively correlated with percent weight gain between day 24–28. Eight taxa were positively correlated with percent weight gain between day 28–35, of which four were Clostridia. Only Lactobacillus reuteri was positively correlated across both day 24–28 and day 28–35 analyses. Insulin-like growth factor 1 (IGF-1; R2 = 0.61, P < 0.001), glucose-dependent insulinotropic polypeptide (GIP; R2 = 0.20, P < 0.001), glucagon-like peptide 1 (GLP-1; R2 = 0.51, P < 0.001), and glucagon-like peptide 2 (GLP-2; R2 = 0.21, P < 0.001) were significantly associated with the piglet fecal community NMDS, while serotonin showed no significant association (R2 = 0.03, P = 0.15). Higher concentrations of GLP-1 and GLP-2 characterized day 1 fecal communities, while GIP levels had the strongest relationship primarily with samples ordinated with the day 21 cluster.
Conclusions
Demonstration of an association of certain taxa with individual gut peptides at specific ages suggests the potential for the microbiome to elicit changes in the gut enteroendocrine system during early postnatal development in the pig.
Introduction
The gut of the healthy neonatal pig experiences an initial microbial colonization at birth with the adaptation to the external environment in combination with the introduction of a milk-based diet. At weaning, the subsequent transition to solid food imparts another transition for the gut microbiome. Several studies in swine have demonstrated that the physical and metabolic/biochemical interactions of the microbiome with the early postnatal diet result in changes in the relative contribution of various bacterial families to the overall microbial ecology of the host environment [1,2,3,4,5].
How the microbiome interacts with the physiology of the neonatal animal to elicit changes in growth rate is of current interest. Several studies have demonstrated that changes in the microbiome elicit changes in short chain fatty acid production and affect intestinal growth and barrier function [4, 6, 7]. Additionally, interactions of the microbial community with the gut wall result in changes in cytokine and growth factor expression by intestinal cells [8, 9]. These cytokines and growth factors may potentially impact feeding behavior, peripheral metabolism and consequently growth rate.
Growth rate can be affected by environment, diet, behavior and stress (psychological and physiological). Mediators for many of these factors that can impact growth rate can also affect feeding behavior, directly or indirectly, and thus potentially affect the microbiome [2, 10, 11]. Studies with post-weaning pigs have demonstrated a relationship between the microbiome and feed efficiency, therefore growth rate [12]. The interaction of the microbiome with growth rate during the weaning transition has not been thoroughly evaluated. The interaction of the microbiome with feed efficiency and growth rate is most likely through short chain fatty acids, biogenic amines, neuroendocrine secretions, and the enteroendocrine system [13, 14]. Recent work has demonstrated interactions exist between the microbiome and the enteroendocrine system [15, 16]. For example, the bacterial metabolite and quorum sensing molecule indole can stimulate glucagon-like peptide 1 (GLP-1) secretion into the vasculature [17]. Germ free mice have reduced enterochromaffin cell (EC) seratonin levels, while colonization of the gut with donor microbiota restores seratonin levels in EC cells [18]. Butyrate, the bacterial metabolites and short chain fatty acid has been reported to stimulate glucose-dependent insulinotropic polypeptide (GIP) secretion following oral administration [19].
The enteroendocrine system produces a variety of gut peptides. Several of these have been reported to impact feed intake, peripheral tissue metabolism and tissue accretion. Many of these actions are through receptors present on vagal afferents; however, some responses are direct effects of these molecules in the periphery. Gut serotonin enters the bloodstream and promotes insulin secretion in the pancreas and stimulates lipogenesis in adipose tissue and liver, while inhibiting lipolysis [20]. GIP has been demonstrated to stimulate lipid uptake and storage in adipocytes while inhibiting lipolysis and is widely considered an adipogenic hormone [21]. GLP-1 is an incretin hormone that stimulates insulin secretion and inhibits glucagon release from the pancreas [22]. Recently, GLP-1 was reported to alter skeletal muscle remodeling by altering glucose uptake, mitochondrial biogenesis and type 1 fiber formation [23]. GIP, GLP-1 and glucagon-like peptide-2 (GLP-2) inhibit bone resorption and thus may impact bone remodeling [24]. These various actions of the gut peptides suggest they may have important roles in growth.
The present study examined if changes in the piglet microbiome between birth and the adaptation to solid feed during the weaning transition can be correlated either positively or negatively with growth rate and plasma concentrations of the enteroendocrine peptides GIP, GLP-1, GLP-2 and serotonin.
Materials and methods
Animals
Large White × Landrace piglets from 8 litters (L.1-8) were assessed from birth through day 35 of age and were weaned at day 21. Piglets were not provided milk replacer/supplement or creep feed at any point throughout the experiment. The diets were formulated to meet the National Research Council estimate of nutrient requirements. From days 21–28, piglets received Nursery Diet 1 exclusively, followed by Nursery Diet 2 exclusively from days 29–35 (Additional file 1). Piglets were given free access to feed and water, evaluated daily for health, and weights were recorded for each sample time point. Animal health was evaluated by veterinary staff. A healthy state was defined by weight gain trends, lack of respiratory issues, absence of scabbing or wounds, alertness. All piglets used in this study were observed to be healthy. No antibiotics, antifungals, or supplemental additives were administered to the piglets at any time during the experiment. Care and treatment of all pigs were approved by the USDA-ARS Institutional Animal Care and Use Committee of the Beltsville Agricultural Research Center.
Fecal sample collection
Fresh fecal samples were collected into sterile cryovial tubes from the rectum of individual piglets (n = 48) on days 1, 21 (weaning), 24, 28 and 35. Samples were initially flash frozen in liquid nitrogen and stored at − 80 °C until further processing. A total of 243 fecal samples were collected.
Blood collection
Blood (0.5–1.5 mL) was collected from the cranial vena cava, and plasma was harvested following centrifugation of EDTA-treated syringes at days 1, 21, and 35 of age with S-Monovette K3E (Sarstedt, Nümbrecht, GE). Plasma samples were stored at − 20 °C until analyzed.
ELISA
Plasma insulin-like growth factor 1 (IGF-I, ALPCO, Salem, NH, 03079), serotonin (5HT, Enzo Life Sciences Inc., Farmingdale, NY 11735), gastric inhibitory peptide (GIP, Millipore, Billerica, MA 01821), glucagon-like peptide 1 (GLP-1, Millipore, Billerica, MA 01821) and glucagon-like peptide 2 (GLP-2, AssayPro LLC, St. Charles, MO, 63301) were analyzed with commercial ELISA kits according to the manufacturer’s instructions. A standard curve was prepared according to manufacturer’s recommendations using a PL4 log-transformed sigmoidal algorithm (SigmaPlot 14.0, Systat Software Inc., Chicago, IL); unknowns were extrapolated from that standard curve.
DNA extraction and sequencing
DNA was isolated from 0.25 g feces using the MagAttract Power Microbiome Kit (Qiagen, Hilden, Germany) by the Microbial Systems Molecular Biology Laboratory at the University of Michigan. Cells were lysed to isolate DNA using mechanical bead beating for 20 total minutes with 20 frequency/second and extracted using magnetic bead technology according to the Qiagen protocol. The V4 region of the 16S rRNA-encoding gene was amplified from extracted DNA using the barcoded dual-index primers developed previously [25]. 16S regions were sequenced with the Illumina MiSeq Sequencing platform generating 2 × 250 bp paired-end reads. All sequences are publicly available in the NCBI’s sequence read archive (SRA) under BioProjects PRJNA613280 and PRJNA673841.
Sequence processing
Quality filtering, pairing, denoising, amplicon sequence variants (ASVs) determination, and chimera removal was conducted with the DADA2 v. 2020.8.0 [26] plugin in QIIME2 v. 2020.4 [27]. For quality trimming, reverse reads were truncated to 200 bp. Taxonomic classification of the ASVs was performed using the Qiime2 pretrained 16S 515F/806R [28,29,30] from the Silva 138 database [31]. Samples with < 5000 sequences (n = 7) and ASVs identified as archaea, chloroplast, mitochondria, or unassigned were removed. A total of 10,039,464 sequences with an average of 41,315 ± 1919 reads per sample (n = 243) were selected for downstream analysis.
Statistical analysis
Assumptions of samples’ independence, normal distribution and equal variances have been verified for all variables. Equality of variances was tested using Levene’s test in R 4.1.3 (https://r-project.org/). When the data were not normally distributed and the assumption of equal variance has been rejected (serotonin and GLP-1) two-way non-parametric ANOVA (Scheirer–Ray–Hare test) was performed with Dunn’s test as post-hoc analysis (Kruskal–Wallis multiple comparison with p values adjusted with Benjamin-Hochberg method) in R. Data characterized by normal distribution and equal variances were analyzed by two-way ANOVA using the GLM procedure in SAS (SAS, Inc., Cary, NC). Least significant differences were used to determine which means differed significantly. The model included the main effects of age, sex, pig and sow. Sex and pig had no effect (P > 0.05) for any of the measured parameters, so the only two-way interaction tested was age (or age range for ADG) x sow.
Non-metric multidimensional scaling (NMDS) of piglet fecal communities was performed using the R-packages vegan (https://CRAN.R-project.org/package=vegan) and phyloseq [32] with the metaMDS function (k = 2, autotransform = False, trymax = 1000) on cumulative sum scaling transformed (CSS) ASVs and Bray–Curtis dissimilarity distances. Permutational analysis of variance (PERMANOVA) was used to determine the main effect of growth time points on the piglet fecal microbial communities using the adonis function (permutations = 1000) in vegan. The strata option (strata = piglet) was used to account for repeated sampling of individual piglets. ASVs present < 1% of samples were removed prior to analysis to reduce noise and the occurrences of rare ASVs. Concentrations of serotonin, GIP, GLP-1, and IGF-1 concentrations were log transformed; a pseudocount of 1 was added to GLP-2 log transformations to prevent negative values. The envfit function (perm = 999) in vegan was then used to fit peptide vectors to the piglet fecal community ordination plot using regression and was visualized using ggordiplots (http://github.com/jfq3/ggordiplots) and ggplot2 packages [33]. Interaction network analyses among centered-log-ratio transformed (CLR) ASVs and log transformed peptide concentrations at each time point were conducted with the Julia based FlashWeave program (sensitive = true, heterogenous = false) [34] which uses a local-to-global learning approach to test for conditional independence and predict direct relationships among microbial taxa and metadata variables. Prior to transformations and network analysis, ASVs in < 30% of samples were filtered from samples to remove rare ASVs and avoid overabundance of zeroes in log ratio transformations, which can lead to spurious results. Network plots were visualized using the Cytoscape v. 3.8.0 [35]. Spearman rank correlations between relative abundance changes in ASVs and percent weight gain in piglets over consecutive time points were performed on total sum scale normalized (TSS) ASVs using the function rcorr (type = “spearman”) in the R-package corrplot (https://github.com/taiyun/corrplot). To reduce correlation bias of zero-inflated data, ASVs present < 50% of the samples at each time point comparison were removed from analysis. Corresponding P values were calculated using the function cor.test and corrected for the false discovery rate (FDR). Correlation plots were visualized using the corrplot function. To visualize relationships between taxa, a neighbor-joining phylogenetic tree [36] based on the multiple sequence comparison by the log-expectation [37] aligned V4 region of 16S rRNA genes were constructed in MEGAX [38]. All analyses in this study used a significance threshold of P < 0.05.
Results
Animal weight gain with age
Pigs used in this study were born at an average birth weight of 1.66 ± 0.04 kg (Fig. 1). Weaning at day 21 of age was accompanied by a mean body weight of 6.94 ± 0.16 kg and an average daily gain (ADG) of 248 ± 8 g/d between day 1 and day 21 of age (Table 1). Weaning resulted in a 49% reduction in average daily gain at day 24 (P < 0.05) relative to day 21 and no change in body weight (P > 0.05). Body weight and average daily gain increased 17% and 175%, respectively, at day 28 relative to day 24 (P < 0.05). The heaviest body weight was observed at day 35 and accompanied by a further increase of 32% in average daily gain between day 28 and day 35 relative to day 24 to day 28 (P < 0.05).
Plasma peptide concentrations
Analysis of plasma peptides demonstrated an age-related increase in IGF-1 from day 1 to day 35 (Fig. 2a, P < 0.001). In contrast, no change in plasma seratonin was detected with age (Fig. 2b, P > 0.05). Plasma GIP did not differ at day 1 and day 21 of age (P > 0.05) but declined by 46% at day 35 (Fig. 2c, P < 0.01). Plasma GLP-1 concentration was 274 ng/mL at d1 in these pigs and dropped by 87–35 ng/mL at day 21, but with no further change at day 35 of age (Fig. 2d, P < 0.05). The plasma concentration of GLP-2 also declined between day 1 and day 21 (Fig. 2e, P < 0.001), but with no further change at day 35 (P > 0.05). Only plasma IGF-1 demonstrated a correlation with body weight in these pigs between day 1 and day 21 (R2 = 0.422, P < 0.0001) or day 1 and day 35 of age (R2 = 0.633; P < 0.001) of the plasma peptides quantified. Average daily gain was also correlated with plasma IGF-1 between day 1 and day 21 (R2 = 0.247, P < 0.001).
Plasma peptide concentrations in pigs at day 1, day 21 and day 35 of age. a Plasma Insulin-like growth factor 1 (IGF-1). Significant age effect (P < 0.0001); significant sow effect (P < 0.0001). b Plasma Seratonin. Significant sow effect (P < 0.0001) and age × sow effect (P < 0.0001). c Plasma Gastric inhibitory peptide (GIP). Significant age effect (P < 0.0001), sow effect (P < 0.0001) and age × sow effect (P < 0.0001). d Plasma Glucagon-like peptide 1 (GLP-1). Significant age effect (P < 0.05), sow effect (P < 0.01) and age × sow effect (P < 0.01). e Plasma Glucagon-like peptide 2 (GLP-2). Significant age effect (P < 0.0001) and age × sow effect (P < 0.001). Legend: Bars represent means and error bars represent standard errors. Means not sharing a common superscript letter are significantly different for age (n = 48; P < 0.05)
Correlation between ASVs and weight gain over the weaning transition
ASV analysis did not detect any significant changes in the microbiome between day 1 and day 21 of age. However, Spearman Rank correlations between weight gain and the microbiome revealed that Clostridium scindens was unique (among those bacteria detected by 16S rRNA analysis) in being associated with growth in the preweaning pig (r = 0.604, P = 0.00004, FDR = 0.00119).
The time interval from day 21 to day 24 had the highest number of significant correlations (n = 25) between the relative abundance differences in taxa over time and corresponding percent weight gain compared to time periods from day 24 to day 28 (n = 8) and from day 28 to day 35 (n = 16; Fig. 3). Between day 21 and day 24, increased relative abundances of ASVs belonging to families Prevotellaceae NK3B31 (\(\rho\) = 0.65, P < 0.001), Veillonellaceae (\(\rho\) = 0.63, P < 0.001) and Rikenellaceae RC9 (\(\rho\) = 0.62, P < 0.001) were significantly correlated with an increase in percent weight gain. In contrast, increased relative abundances of ASVs Parabacteroides merdae (\(\rho\) = − 0.58, P < 0.001), Bacteroides (\(\rho\) = − 0.58, P < 0.001) and Cloacibacillus porcorum (\(\rho\) = − 0.57, P < 0.001) were significantly correlated with a decrease in percent weight gain. Bacteroidia, Clostridia and Negativicutes accounted for 19 of the 25 identified genera or species that affected ADG between day 21 and day 24. Throughout the weaning transition, increased relative abundances of ASVs belonging to the Lactobacillus genera all showed significant and positive correlations with weight gain, while increased relative abundances of ASVs belonging to the Christenellaceae R7 group were significantly and negatively correlated with percent weight gain. ASV UCG-002 showed a significant positive correlation with percent weight gain between day 21 and day 24 (\(\rho\) = 0.55, P < 0.001) but a significant negative correlation with percent weight gain between day 28 and day 35 (\(\rho\) = − 0.52, P = 0.01). Seven bacterial genera or species were positively correlated with percent weight gain between day 24 and day 28, with only one genus correlating with percent weight gain for both day 21 to day 24 and day 24 to day 28: Megaspaera. Analysis of the correlations from day 28 to day 35 indicated the increased relative abundances of eight taxa that were positively correlated with percent weight gain, of which four were Clostridia. Only Lactobacillus reuteri was positively correlated across both day 24 to day 28 and day 28 day 35 analyses. Among the eight taxa that were negatively correlated with percent weight gain in the day 28 to day 35 analysis, all but Succiniaclasticum were of the Clostridia class.
Correlations between weight gain and ASVs. Legend: Spearman rank correlation plot showing significant correlations (P < 0.05, FDR corrected) between change in percent relative abundance and percent relative weight gain between time points (D21-D24, D24-D28 and D28-D35). Red circles indicate negative and blue circles indicate positive correlations. The size of the circle indicates the strength of the correlation. Neighbor-joining tree with bootstrap values (bootstrap = 100, values < 50 not shown) based on MUSCLE alignment of 16S rRNA v4 gene region depict phylogenetic relationship between ASVs
Piglet fecal community structure and peptides associated with age
As expected, age had an effect on the piglet fecal microbiomes at day 1, day 21 and day 35 (PERMANOVA F = 38.94, R2 = 0.37, P < 0.001) and samples (n = 136) showed distinct clusters by age (Fig. 4, stress = 0.093, k = 2). Peptides related to the sample plot distributions were fitted to the ordination using regression. IGF-1 (R2 = 0.61, P < 0.001), GIP (R2 = 0.20, P < 0.001), GLP-1 (R2 = 0.51, P < 0.001), and GLP-2 (R2 = 0.21, P < 0.001) were significantly associated with the piglet fecal community NMDS, while serotonin showed no significant association (R2 = 0.03, P = 0.15). Higher plasma concentrations of GLP-1 and GLP-2 were associated with day 1 fecal communities, while day 35 fecal samples were more closely associated with higher plasma levels of IGF-1. Plasma GIP levels had the strongest relationship primarily with fecal communities ordinated with the day 21 cluster.
Non-metric multidimensional scaling (NMDS) plot of the piglet fecal community structure at different time points based on Bray–Curtis dissimilarity matrix. Legend: The effect of age on community was tested using permutational analysis of variance (PERMANOVA). The envfit function was used to fit peptide variables to the ordination plot. Ellipses indicate 1 standard error from centroid. Arrow lengths are scaled at 0.5. Only peptides that showed significant correlation (P < 0.05) with the unconstrained ordination axes were plotted
ASVs associated with plasma peptides and weight gain
Direct associations between ASVs in the piglet fecal microbiome and peptide concentrations or ADG were determined using a network-based approach for each time point at day 1, day 21 and day 35 (Fig. 5). Weights for each association are found in Additional file 2. Only a limited number of relationships were discovered, therefore the significant relationships between gut peptides and the microbiome were extracted from the overall networks at each age for ease of visualization. At day 1, no significant associations between fecal taxa and gut peptide concentrations or weight gain were determined; however, there was a positive association between IGF-1 and weight gain at day 1 of age. Day 21 microbiomes showed the most direct associations between taxa and peptides or weight gain (2 negative, 5 positive) including a positive association between an unclassified species from genus UCG-002 (family Oscillospiraceae) and IGF-1, positive associations between an unclassified Akkermansia species and GIP as well as between Clostridiales and GIP. Negative associations were identified between an unclassified Desulfovibrio species and weight gain, while GLP-1 was negatively associated with Actinomyces at day 21 of age.
A–C Network plots of microbiota associations with peptides and average daily gain (ADG). Legend: Network plots showing direct associations among ASVs, peptides and ADG at time points (A) day 1, (B) day 21 and (C) day 35 using the FlashWeave program. Only significant associations between ASVs and metadata variables are shown (P < 0.01, FDR corrected). Edge color indicates association: negative (red), positive (blue); edge thickness indicates association weight. ASV node size is proportional to the mean centered log ratio abundance
At day 35, 3 positive associations and 2 negative associations between taxa and peptides were found. GIP was positively associated with L. reuteri and negatively with Clostridium sensu stricto 1. GLP-1 was positively associated with Acidaminococcus at day 35 of age. For all time points, peptide IGF-1 and weight gain were positively associated with each other.
Discussion
Growth in the pig is characterized by a moderate preweaning growth rate with competition for nutrients amongst littermates between birth and weaning and dependent on the lactational capacity of the sow. This is often followed by a post weaning lag in growth rate, which was seen in this study between day 21 and day 24 of age when the average daily gain dropped by almost half during a period of major changes in the gut microbiome and gut permeability associated with a transition to a corn/soy based diet from the milk based diet. Successful transition to the weaning diet between day 24 and day 28 resulted in a rapid increase in average daily weight gain to compensate for the reduced energy intake and consequent reduced growth rate immediately following weaning (day 21 to day 24).
As mentioned, C. scindens was unique in the preweaning microbiome with association with preweaning growth. The high correlation of C. scindens with preweaning weight gain suggests a beneficial role for C. scindens in net nutrient uptake from the gut. C. scindens may interact with other components of the microbiome to promote net nutrient uptake or alter the gut ecology to augment gut health and consequent animal growth. Previous research has indicated that C. scindens can either directly inhibit the growth of C. difficile or indirectly inhibit C. difficile growth through secondary transformation of bile acids by hydroxylation [39, 40], and thus, promote overall gut health. In addition, these 7-hyroxylated bile acids can function as signaling molecules [41], for example as agonists for the Takeda G-protein receptor 5 (TGR5) and may indirectly contribute to the regulation of energy metabolism and inflammation [42]. C. scindens has also been demonstrated to induce GLP-1 expression; GLP-1 has an insulinotropic effect which suggests a promotion of net energy accretion in the pig [22]. Therefore, these data suggest C. scindens may serve as a potential marker for a gut microbial community beneficial to growth of the preweaning pig. The correlation between C. scindens and weight gain was gone after day 21 as the microbiome adapted to the dietary and physiological changes that occur with weaning. This study utilized the Illumina MiSeq platform, targeting the V4 region of the bacterial 16S rRNA gene. This region has been accepted as the most accurate target for bacterial taxonomy [25]. A caveat of this study is that the microbiome field continues to evolve through advances in sequencing technology, such as full-length 16S rRNA gene sequencing and metagenomics. While full-length 16S rRNA gene sequencing can allow for better taxonomic resolution, and thus, more accurate compositional data, the cost remains restrictive [43, 44]. Metagenomics, which investigates the total nucleotide sequences in a sample, often misses numerically inferior species and cannot investigate important members of the rare biosphere without additional targeted techniques [45, 46]. As costs continue to decrease for sequencing techniques, and updates are made to bacterial genome and pathway databases, microbiome data will become even more informative.
The changes in weight gain accompanied very large changes in microbiome dynamics following weaning. The reduction in weight gain between day 21 and day 24 following weaning was accompanied by both negative and positive correlations between weight gain and specific bacterial taxa. The high correlation between Prevotellaceae and weight gain was relatively unexpected due to the large reduction in growth rate between these ages and a presumed reduction in feed intake. The changes in the bacterial community as a result of the transition from a diet comprised of milk oligosaccharides used by Bacteroidaceae, prior to weaning, to the complex polysaccharides utilized by Prevotellaceae post-weaning has been previously reported [47,48,49] and may account for the large negative and positive correlations with weight gain observed for these two members of the bacterial community, respectively, between day 21 and day 24 of age. The positive association with Prevotellaceae may also be the consequence of the unique mucin glycoprotein degrading sulfoglycosidase that permits it to grow during the early weaning transition [50]. Following the transition to the complex carbohydrate based post-weaning diet, no correlations between Prevotellaceae and weight gain were detected at day 28 when weight gain had improved, suggesting that Prevotellaceae may be an early adapter to changes in diet in these swine. However, other members of the bacterial community also contribute in the present study during the transition to a more complex grain based diet as demonstrated by the positive correlations of Veillonellaceae, Butyricicoccaceae, Spirochaetaceae, Muribaculaceae and Oscillospiraceae with weight gain between day 21 and day 24. Veillonellaceae, Butyricicoccaceae, Spirochaetaceae and Oscillospiraceae are all short chain fatty acid producing bacterial families [51, 52] and in the present study may be important in the early adaption to metabolizing complex polysaccharides from grain starches. Muribaculaceae has been a largely uncharacterized clade, but whose abundance has recently been associated with an increase in fecal proprionate concentrations in mice [53] and swine [54], and genomes assembled from metagenomes indicate Muribaculaceae may express the necessary metabolic pathways to ferment polysaccharides into propionate, acetate and succinate [55]. Rikenellaceae RC9 has been associated with negative impacts on health and growth in several studies in mice [56,57,58], while the present study suggests a positive association with growth during the weaning adaptation.
The time from day 24 to day 28 demonstrated a rapid increase in weight gain, indicating the pigs’ successful adaption to the postweaning diet. This timespan was characterized by a smaller number of bacterial families correlating with weight gain; these correlations with weight gain trended weaker with only Lactobacillaceae and Ruminococceae correlations improved with ADG between day 24 and day 28 from day 21 to day 24. A different milieu of bacteria genera and species accounted for positive and negative correlations with body weight gain between day 24 and day 28 than those during the period between day 21 and day 24, with Blautia and L. reuteri demonstrating the highest positive correlations with body weight gain. Blautia has previously been correlated with increases in postweaning feed intake in swine [59], which can contribute to an increase in ADG; while specific Lactobacillus reuteri strains have been demonstrated to have probiotic properties in vitro [60, 61] and by oral administration to alter the composition of the microbial community within the gut of the newly weaned pig [62] as well as promote an increase in ADG in newly weaned pigs by enhancing intestinal barrier function [63]. However, the L. reuteri strains in the present study have not been identified. Christensenellaceae R-7 were the only genus to have a negative correlation with ADG between day 24 and day 28. Previous research [64, 65] has indicated that unclassified members of Christensenellaceae are associated with higher feed efficiency in swine post-weaning and function to maintain gastrointestinal tract structure and function, which suggests that further study is necessary to elucidate the specific role of Christensenellaceae R-7 relative to other genera of Christensenellaceae in gut health and animal growth. The large reduction in the number of negative correlations between taxa from day 21 to day 24 relative to day 24 to day 28 suggests that the pig successfully transitions to the Phase 1 diet by no later than day 28 of age as only one taxa is negatively associated with growth rate in that analysis group.
Weight gain continued to increase from day 28 to day 35 as the pigs grew more rapidly than earlier ages, suggesting the swine had successfully adapted to the dietary and environmental changes associated with weaning. Pigs were transitioned from nursery diet 1 to nursery diet 2 at day 28. This change in nursery diets follows standard practice in the swine industry to reduce feed costs by eliminating the highly digestible sprayed dried plasma proteins and reducing the amount of calf milk replacer found in the weanling diet, substituting soybean meal as a protein source and increasing the amount of corn in the diet as an energy source. Correlations between ADG and the microbial community experienced further changes during this timeframe, coinciding with the transition from nursery diet 1 to nursery diet 2. Previous research has suggested that unidentified Prevotella and Faecalibacterium can respond to dietary protein [66] and thus may account for their positive association with weight gain. Further research is needed to identify if soy protein can affect microbiome composition with phased nursery diet feeding.
The pathway analysis provided a few insights into the relationship of the microbiome to plasma peptides that can originate from the gastrointestinal tract and associated organs of the digestive system. The newborn gut ecology is undergoing significant changes with the introduction of dietary milk and exposure to the extrauterine environment, so relationships between the gut peptides and the microbiome at day 1 of age were not expected. However, the analysis demonstrated that GLP-1 and GLP-2 were both associated with the day 1 piglet fecal community microbiome. Both peptides were at their highest plasma concentrations at day 1 when comparing the three sampling ages. GLP-1 and GLP-2 both respond to nutrient ingestion and absorption [67, 68] and the response to the initial nutrient presentation during colostrum feeding may induce a maximal secretory response. Therefore, the association of plasma GLP-1 and GLP-2 with the day 1 fecal microbiome may be the consequence of the initial interactions of the adapting and developing microbiome with the digesta from colostrum and the subsequent absorption of simple sugars and lipids through the intestine or these nutrients may be bypassing the microbiome for direct absorption. IGF-1 was not associated with the microbiome at day 1 but was correlated with body weight. Analysis of plasma IGF-1 was included in this study due to its known associations with changes in body weight with age [69]. Previous research has shown that IGF-1 in newborn healthy humans is correlated to body weight [70], however the only study in swine to examine this relationship compared normal and growth retarded newborn pigs, finding a negative association [71], and no analysis was performed within the normal sized piglet group.
Plasma GIP concentration was associated with the day 21 fecal microbiome community. The microbiome at day 21 is much more complex that in the newborn pig [1, 4, 5, 11] and is fully capable of metabolic functions that produce numerous microbial metabolites that can communicate with the enteroendocrine system [15, 16]. Butyrate is one microbial metabolite that has been demonstrated to stimulate GIP release [19]. GIP functions to stimulate insulin secretion by the pancreas and to promote lipid accretion in adipose tissue. The neonatal pig is born with limited fat stores and rapidly accumulates adipose tissue lipids through weaning [72]. Thus, GIP may function in the preweaning pig to stimulate lipogenesis by adipose tissue in response to interactions with the microbiome.
Association of plasma IGF-1 with the microbiome has been previously reported [73]. Yan et al. [74] colonized the gut of germ-free mice and found an increase in serum IGF-1 and bone growth, while providing short chain fatty acids to antibiotic treated mice restored growth rates and increased serum IGF-1. If probiotics can be developed to specifically target microbial signaling for IGF-1 release, the potential for improving efficiency of production is significant as promotion of lean, fast growth is always a priority for the swine industry. Demonstration of an association of specific taxa with specific gut peptides in the present study suggests the potential for developing such probiotics that may result in a specific enteroendocrine response. However, further studies are necessary to elucidate these specific relationships.
Conclusions
Demonstration of an association of certain taxa with individual gut peptides at specific ages suggests the potential for the microbiome to elicit changes in the gut enteroendocrine system during early postnatal development in the pig.
Availability of data and materials
All sequences are publicly available in the NCBI’s Sequence Read Archive (SRA) under BioProjects PRJNA613280 (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA613280) and PRJNA673841 (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA673841).
Abbreviations
- ADG:
-
Average daily gain
- ASV:
-
Amplicon sequence variant
- CLR:
-
Centered-log-ratio
- CSS:
-
Cumulative sum scaling
- FDR:
-
False discovery rate
- g:
-
Grams
- GIP:
-
Glucose-dependent insulinotropic polypeptide
- GLP-1:
-
Glucagon-like peptide 1
- GLP-2:
-
Glucagon-like peptide 2
- IGF-1:
-
Insulin-like growth factor 1
- kg:
-
Kilogram
- NMDS:
-
Non-metric multidimensional scaling
- PERMANOVA:
-
Permutational analysis of variance
- TGR5:
-
Takeda G-protein receptor 5
References
Wang X, et al. Longitudinal investigation of the swine gut microbiome from birth to market reveals stage and growth performance associated bacteria. Microbiome. 2019;7(1):109.
Holman DB, et al. Meta-analysis to define a core microbiota in the swine gut. mSystems. 2017;2(3):e00004-17.
Wu J, et al. Dietary nutrients shape gut microbes and intestinal mucosa via epigenetic modifications. Crit Rev Food Sci Nutr. 2020;62:1–15.
Qi R, et al. Changes of gut microbiota and its correlation with short chain fatty acids and bioamine in piglets at the early growth stage. Front Vet Sci. 2020;7:617259.
Arfken AM, Frey JF, Summers KL. Temporal dynamics of the gut bacteriome and mycobiome in the weanling pig. Microorganisms. 2020;8(6):868.
LeBlanc JG, et al. Beneficial effects on host energy metabolism of short-chain fatty acids and vitamins produced by commensal and probiotic bacteria. Microb Cell Fact. 2017;16(1):79.
Xiong X, et al. Nutritional intervention for the intestinal development and health of weaned pigs. Front Vet Sci. 2019;6:46.
Broom LJ, Kogut MH. Gut immunity: its development and reasons and opportunities for modulation in monogastric production animals. Anim Health Res Rev. 2018;19(1):46–52.
Xin J, et al. Probiotic Lactobacillus johnsonii BS15 promotes growth performance, intestinal immunity, and gut microbiota in piglets. Probiotics Antimicrob Proteins. 2020;12(1):184–93.
Isaacson R, Kim HB. The intestinal microbiome of the pig. Anim Health Res Rev. 2012;13(1):100–9.
Chen L, et al. The maturing development of gut microbiota in commercial piglets during the weaning transition. Front Microbiol. 2017;8:1688.
Gardiner GE, Metzler-Zebeli BU, Lawlor PG. Impact of intestinal microbiota on growth and feed efficiency in pigs: a review. Microorganisms. 2020;8(12):1886.
Silver M, Fowden AL. Pituitary-adrenocortical activity in the fetal pig in the last third of gestation. Q J Exp Physiol. 1989;74(2):197–206.
Steinert RE, et al. Digestive physiology of the pig symposium: secretion of gastrointestinal hormones and eating control. J Anim Sci. 2013;91(5):1963–73.
Martin AM, et al. The influence of the gut microbiome on host metabolism through the regulation of gut hormone release. Front Physiol. 2019;10:428.
Lupien-Meilleur J, et al. Interplay between gut microbiota and gastrointestinal peptides: potential outcomes on the regulation of glucose control. Can J Diabetes. 2020;44(4):359–67.
Kuhre RE, et al. Bile acids are important direct and indirect regulators of the secretion of appetite- and metabolism-regulating hormones from the gut and pancreas. Mol Metab. 2018;11:84–95.
Yano JM, et al. Indigenous bacteria from the gut microbiota regulate host serotonin biosynthesis. Cell. 2015;161(2):264–76.
Lin HV, et al. Butyrate and propionate protect against diet-induced obesity and regulate gut hormones via free fatty acid receptor 3-independent mechanisms. PLoS ONE. 2012;7(4):e35240.
Yabut JM, et al. Emerging roles for serotonin in regulating metabolism: new implications for an ancient molecule. Endocr Rev. 2019;40(4):1092–107.
Gogebakan O, et al. Glucose-dependent insulinotropic polypeptide reduces fat-specific expression and activity of 11beta-hydroxysteroid dehydrogenase type 1 and inhibits release of free fatty acids. Diabetes. 2012;61(2):292–300.
Muller TD, et al. Glucagon-like peptide 1 (GLP-1). Mol Metab. 2019;30:72–130.
Wu L, et al. GLP-1 regulates exercise endurance and skeletal muscle remodeling via GLP-1R/AMPK pathway. Biochim Biophys Acta Mol Cell Res. 2022;1869(9):119300.
Schiellerup SP, et al. Gut hormones and their effect on bone metabolism. Potential drug therapies in future osteoporosis treatment. Front Endocrinol (Lausanne). 2019;10:75.
Kozich JJ, et al. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79(17):5112–20.
Callahan BJ, et al. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581–3.
Bolyen E, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7.
Bokulich A, Parker W. Data models, representation and adequacy-for-purpose. Eur J Philos Sci. 2021;11(1):31.
Bokulich NA, et al. Optimizing taxonmic classification of marker-gen amplicon sequences with QIIME 2's q2-feature-classifier plugin. Microbiome 2018;6.
Bokulich NA, et al. bokulich-lab/RESCRIPt: 2020.11. 2020.
Yilmaz P, et al. The SILVA and “all-species living tree project (LTP)” taxonomic frameworks. Nucleic Acids Res. 2014;42(Database issue):D643–8.
McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8(4):e61217.
Wickham H. ggplot2: elegant graphics for data analysis. 2nd ed. Berlin: Springer; 2016.
Tackmann J, Matias Rodrigues JF, von Mering C. Rapid inference of direct interactions in large-scale ecological networks from heterogeneous microbial sequencing data. Cell Syst. 2019;9(3):286–96.
Shannon P, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25.
Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 2004;5:113.
Kumar S, et al. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9.
Buffie CG, et al. Precision microbiome reconstitution restores bile acid mediated resistance to Clostridium difficile. Nature. 2015;517(7533):205–8.
Marion S, et al. In vitro and in vivo characterization of Clostridium scindens bile acid transformations. Gut Microbes. 2019;10(4):481–503.
Keitel V, Stindt J, Haussinger D. Bile acid-activated receptors: GPBAR1 (TGR5) and other G protein-coupled receptors. Handb Exp Pharmacol. 2019;256:19–49.
Schaap FG, Trauner M, Jansen PL. Bile acid receptors as targets for drug development. Nat Rev Gastroenterol Hepatol. 2014;11(1):55–67.
Katiraei S, et al. Evaluation of full-length versus V4-region 16S rRNA sequencing for phylogenetic analysis of mouse intestinal microbiota after a dietary intervention. Curr Microbiol. 2022;79(9):276.
Gotschlich EC, Colbert RA, Gill T. Methods in microbiome research: past, present, and future. Best Pract Res Clin Rheumatol. 2019;33(6):101498.
Huffnagle GB, Noverr MC. The emerging world of the fungal microbiome. Trends Microbiol. 2013;21(7):334–41.
Liu S, et al. Opportunities and challenges of using metagenomic data to bring uncultured microbes into cultivation. Microbiome. 2022;10(1):76.
Guevarra RB, 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:54.
Frese SA, et al. Diet shapes the gut microbiome of pigs during nursing and weaning. Microbiome. 2015;3:28.
Summers KL, et al. The piglet mycobiome during the weaning transition: a pilot study. J Anim Sci. 2019;97(7):2889–900.
Rho JH, et al. A novel mechanism for desulfation of mucin: identification and cloning of a mucin-desulfating glycosidase (sulfoglycosidase) from Prevotella strain RS2. J Bacteriol. 2005;187(5):1543–51.
Vital M, Howe AC, Tiedje JM. Revealing the bacterial butyrate synthesis pathways by analyzing (meta)genomic data. MBio. 2014;5(2):e00889.
Beaumont M, et al. Gut microbiota-derived metabolite signature in suckling and weaned piglets. J Proteome Res. 2021;20(1):982–94.
Smith BJ, et al. Changes in the gut microbiome and fermentation products concurrent with enhanced longevity in acarbose-treated mice. BMC Microbiol. 2019;19(1):130.
Wu Y, et al. Short administration of combined prebiotics improved microbial colonization, gut barrier, and growth performance of neonatal piglets. ACS Omega. 2020;5(32):20506–16.
Ormerod KL, et al. Genomic characterization of the uncultured Bacteroidales family S24–7 inhabiting the guts of homeothermic animals. Microbiome. 2016;4(1):36.
Gao X, et al. Correlations between alpha-linolenic acid-improved multitissue homeostasis and gut microbiota in mice fed a high-fat diet. mSystems. 2020;5(6):e0039120.
Sun L, et al. Cecal gut microbiota and metabolites might contribute to the severity of acute myocardial ischemia by impacting the intestinal permeability, oxidative stress, and energy metabolism. Front Microbiol. 2019;10:1745.
Zheng Y, et al. Lycium barbarum mitigates radiation injury via regulation of the immune function, gut microbiota, and related metabolites. Biomed Pharmacother. 2021;139:111654.
Renganathan V, Madyastha KM. Linalyl acetate is metabolized by pseudomonas incognita with the acetoxy group intact. Appl Environ Microbiol. 1983;45(1):6–15.
Simpson JM, et al. Denaturing gradient gel electrophoresis analysis of 16S ribosomal DNA amplicons to monitor changes in fecal bacterial populations of weaning pigs after introduction of Lactobacillus reuteri strain MM53. Appl Environ Microbiol. 2000;66(11):4705–14.
Zhang WM, et al. Lactobacillus reuteri glyceraldehyde-3-phosphate dehydrogenase functions in adhesion to intestinal epithelial cells. Can J Microbiol. 2015;61(5):373–80.
Zhang D, et al. Changes in the diversity and composition of gut microbiota of weaned piglets after oral administration of Lactobacillus or an antibiotic. Appl Microbiol Biotechnol. 2016;100(23):10081–93.
Tian Z, et al. Effect of long-term dietary probiotic Lactobacillus reuteri 1 or antibiotics on meat quality, muscular amino acids and fatty acids in pigs. Meat Sci. 2021;171:108234.
McCormack UM, et al. Exploring a possible link between the intestinal microbiota and feed efficiency in pigs. Appl Environ Microbiol. 2017;83(15):e00380-17.
Jenkins S, et al. Relationships between diets different in fibre type and content with growth, Escherichia coli shedding, and faecal microbial diversity after weaning. Anim Prod Sci. 2015;55:1451.
Hou L, et al. Effects of protein restriction and subsequent realimentation on body composition, gut microbiota and metabolite profiles in weaned piglets. Animals (Basel). 2021;11(3):686.
Knapper JM, Morgan LM, Fletcher JM. Nutrient-induced secretion and metabolic effects of glucose-dependent insulinotropic polypeptide and glucagon-like peptide-1. Proc Nutr Soc. 1996;55(1B):291–305.
Burrin DG, et al. Glucagon-like peptide 2 dose-dependently activates intestinal cell survival and proliferation in neonatal piglets. Endocrinology. 2005;146(1):22–32.
Li JQ, et al. Genetic effects of IGF-1 gene on the performance in Landrace x Lantang pig resource population. Yi Chuan Xue Bao. 2003;30(9):835–9.
Ong K, et al. Size at birth and cord blood levels of insulin, insulin-like growth factor I (IGF-I), IGF-II, IGF-binding protein-1 (IGFBP-1), IGFBP-3, and the soluble IGF-II/mannose-6-phosphate receptor in term human infants. The ALSPAC Study Team. Avon Longitudinal Study of Pregnancy and Childhood. J Clin Endocrinol Metab. 2000;85(11):4266–9.
Bauer R, et al. Body weight distribution and organ size in newborn swine (sus scrofa domestica)—a study describing an animal model for asymmetrical intrauterine growth retardation. Exp Toxicol Pathol. 1998;50(1):59–65.
Mersmann HJ, et al. Adipose tissue composition and lipogenic capacity in developing swine. Am J Physiol. 1973;224(5):1130–5.
Jensen EA, et al. Crosstalk between the growth hormone/insulin-like growth factor-1 axis and the gut microbiome: a new frontier for microbial endocrinology. Growth Horm IGF Res. 2020;53–54:101333.
Yan J, et al. Gut microbiota induce IGF-1 and promote bone formation and growth. Proc Natl Acad Sci U S A. 2016;113(47):E7554–63.
Acknowledgements
The authors wish to thank the USDA Beltsville Agricultural Research Center farm employees for their assistance with animal care and handling. The authors would like to thank Craig Storuzuk, DVM, and Emily Beall for veterinary care of the experimental animals in this study. The authors thank Monika Proszkowiec-Weglarz, Ph.D. for assistance with the statistical analysis of animal performance and plasma peptides. This research was supported by work performed by the University of Michigan Microbial Systems Molecular Biology Laboratory. The authors wish to thank Ms. J. Foster-Frey for assistance with sample processing. Dr. A.M. Arfken is currently affiliated with the CDC Environmental and Microbiology Engineering Lab in the Waterborne Disease Prevention Branch through the Oak Ridge Institute for Science and Education, Oak Ridge Associated Universities, Oak Ridge TN.
Funding
No grants or outside funding were used for this study.
Author information
Authors and Affiliations
Contributions
TGR, AMA and KLS contributed to the design of the study. TGR, AMA and KLS performed animal research and sample collection. TGR performed the ELISA analysis. AMA performed statistical analysis; TGR and AMA wrote the first draft of the manuscript, all authors contributed to manuscript revision, read, and approved the submitted version.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Care and treatment of all pigs were approved by the USDA-ARS Institutional Animal Care and Use Committee of the Beltsville Agricultural Research Center.
Consent for publication
Not applicable.
Competing interests
The authors declare no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Mention of trade name, proprietary product or vendor does not constitute a guarantee or warranty of the product by the U.S. Department of Agriculture or imply its approval to the exclusion of other products or vendors that also may be suitable
Supplementary Information
Additional file 1.
Nursery diet composition.
Additional file 2.
Analyses of ASVs, time points, and animal weight by Spearman rank correlation.
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/.
About this article
Cite this article
Ramsay, T.G., Arfken, A.M. & Summers, K.L. Enteroendocrine peptides, growth, and the microbiome during the porcine weaning transition. anim microbiome 4, 56 (2022). https://doi.org/10.1186/s42523-022-00206-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s42523-022-00206-8
Keywords
- Piglet
- Weaning
- Microbiome
- Bacteriome
- Enteroendocrine
- GLP-1
- GLP-2
- GIP
- IGF-1