Weight Gain, Diarrhea Incidence, and mRNA Expression of Inflammatory Cytokines and Barrier Proteins in Suckling Bamei Pigs
As can be summarized from Table 1, piglets fed the supplement (ES) demonstrated significantly increased (p<0.05) final body weight and ADG compared with piglets in the Con group, while diarrhea incidence was significantly decreased for piglets fed ES compared with Con piglets (p<0.05). Early dietary supplementation had a substantial effect on the expression of inflammatory cytokines and tight junction protein genes in the jejunum (Figure 1). In the ES group, TLR4, ZO-1, Occludin and Claudin-1 mRNA levels were higher compared to the Con group (p<0.05). For Bamei piglets fed ES, jejunum TNFα and IL-8 mRNA concentration were significantly lower compared to Bamei piglets in the Con group (p<0.05).
ES Induced Changes to the Jejunum Microbiome of Bamei Suckling Piglets
An average of 61,438 clean tags per sample were obtained from the 35 jejunal content samples analyzed. A total of 672 OTUs based on a 97% sequence similarity were identified from these sequences. Of these 672 OTUs, 498 OTUs were common among both groups (Figure 2A). These 498 OTUs mapped to 20 phyla, 48 classes, 106 orders, 182 families, 350 genera, and 377 species. The alpha diversity was estimated through the diversity index (Shannon and Simpson) and richness estimate (Chao1 and Ace). As can be seen in Figure 2B, the richness estimate (ACE and Chao1) increased significantly for ES group compared to the Con group (p<0.05), whereas the diversity indices (Simpson and Shannon) were similar among both groups (p>0.05). The NMDS plot (Stress = 0.125<0.2), which is used to illustrate the dissimilarity of the microbial community, revealed distinct structures between the Con and ES groups (Figure 2C). Similarly, the jackknifed beta diversity and hierarchical clustering analysis via the UPGMA demonstrated that different groups were clustered into their own groups (Figure 2D). These results suggest the gut microbiota composition in suckling Bamei pigs was being altered by early dietary supplementation.
Here, a total of 20 phylum and 350 genera were identified within the jejunum microbiome. Figures 3A, B report the top 10 most abundant microbes in the piglet’s jejunum intestinal contents for both groups. Firmicutes and Proteobacteria were observed as the predominant phyla in the jejunal microbiota, followed by Bacteroidetes, Chlamydiae, and Actinobacteria. The relative jejunal microbiota abundances at the genus level are presented in Figure 3B. The results demonstrate that Lactobacillus, Clostridium_sensu_stricto_1, Buchnera, Actinobacillus, and Acinetobacter were the predominant genera. Linear discriminant analysis (LDA) effect size taxonomic cladogram (Figure 3C) and LDA value distribution histogram (Figure 3D) were used to determine the jejunal microbiota structure and the predominant bacteria present with the greatest taxa differences between the two groups displayed (i.e. LDA score > 4). The relative abundances of Proteobacteria, Buchnera, Enterobacteriales, Enterobacteriaceae, Alphaproteobacteria, Rickettsiales, Rickettsiaceae, Rickettsia, Bacteroidia, Bacteroidetes, and Bacteroidales were significantly higher for piglets fed ES, while the relative abundances of Veillonellaceae, Selenomonadales, Negativicutes, Romboutsia, Peptostreptococcaceae, Actinobacillus, Pasteurellaceae, Pasteurellales, Firmicutes, Bacilli, Lactobacillales, Lactobacillus, and Lactobacillaceae were reduced compared to piglet’s fed control. The data from the Kruskal-Wallis rank sum test analysis histogram of bacterial phyla and genera data is presented in Figures 3E-F. The relative abundances of Bacteroidetes, Proteobacteria and Firmicutes accounted for more than 1% of the total microbiome. For piglet’s fed ES the relative abundance of Bacteroidetes and Proteobacteria were significantly decreased (p<0.05), and Firmicutes was significantly lower compared with piglets fed Con (p<0.05). Seven (7) genera demonstrated relative abundances of more than 1%, whereas piglet’s fed ES demonstrated lower (p<0.05) Romboutsia, Actinobacillus, Bacteroides and Lactobacillus were lower than piglet’s fed Con.
In order to evaluate the functional capacity of the jejunal bacterial communities for both piglet groups, PICRUSt was used to further analyze the KEGG pathway compositions. The second level KEGG pathway analysis showed that Lipid metabolism, Energy metabolism, Metabolism of terpenoids and polyketides, Infectious diseases: Viral, Cardiovascular diseases, Infectious diseases: Parasitic, Neurodegenerative diseases, Circulatory system, Transport and catabolism, Cancers: Specific types, Endocrine system, Substance dependence, Endocrine and metabolic diseases, and Signal transduction were enriched (p<0.05), while Carbohydrate metabolism, Nucleotide metabolism and Membrane transport were decreased (p<0.05) in the ES group (Figure 4).
Metabolites and Metabolic Pathways Within the Jejunum
To further explore early dietary supplementation influence on jejunal microbiota, the jejunum metabolite concentrations for both groups were analyzed. A total of 283 metabolites were found. These metabolites, including amino acids, carbohydrates, organic acids, lipids, nucleotides, and others, are involved in multiple jejunum biochemical processes for the Bamei piglets. The PCA score plots from were derived from the LC-TOF/MS metabolic profiles of jejunal contents showed separation between piglets fed ES and Con (Figure. 5A, B). As shown in the OPLS-DA score plot (Figure. 5C, D), piglets fed ES compared to groups and piglets fed Con could be separated into distinct clusters according to their metabolic differences (OPLS-DA models +: R2Y = 0.869 and Q2 = 0.471; -: R2Y = 0.836 and Q2 = 0.321). In addition, the permutation test for OPLS-DA demonstrated the Q2 regression line had a negative intercept. Additionally, all R2 and Q2 values on the left were lower than the original points on the right (OPLS-DA validate models +: R2Y = 0.729 and Q2 = -0.669; -: R2Y = 0.5448 and Q2 = -0.6042) (Figure. 5E, F), demonstrating that the OPLS-DA model in the present study is valid.
To identify compound differences between the two treatment groups, the parameters of VIP>1 and p<0.05 were used as a criterion. In the jejunum, 23 compounds (L-Citrulline, Betaine, 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine, 5-Methylcytosine, Cytosine, Glycitein, Daidzein, N-Oleoylethanolamine, L-Histidine, Acetylcarnitine, 1-Myristoyl-sn-glycero-3-phosphocholine, 1-Oleoyl-L-.alpha.-lysophosphatidic acid, 2'-Deoxyinosine, 1-Palmitoyl-sn-glycero-3-phosphocholine, Thioetheramide-PC, PC(16:0/16:0), Linoleoyl ethanolamide, Thymine, Arachidonic Acid (peroxide free), 2'-Deoxyuridine, Genistein, Chenodeoxycholate, 4-Androsten-17.beta.-ol-3-one glucosiduronate) were increased and 35 compounds (Guanosine, Pro-Glu, 2-Hydroxyadenine, N-Acetylmannosamine, N-Acetylneuraminic acid, Pro-Ala, Uridine, N-Acetyl-D-glucosamine, L-Arginine, Pro-Thr, Pro-Phe, MG(18:2(9Z,12Z)/0:0/0:0)[rac], Allopurinol riboside, Ile-Pro, Asp-Leu, Thr-Ala, Riboflavin, Cholic acid, Hypoxanthine, S-Methyl-5'-thioadenosine, Trimethylamine N-oxide, Adenine, 3-Methoxy-4-Hydroxyphenylglycol Sulfate, Muramic acid, all cis-(6,9,12)-Linolenic acid, Lumichrome, Alpha-D-Glucose, D-Mannose, Pantothenate, L-Asparagine, D-Lyxose, L-Threonine, L-Aspartate, Inosine) were decreased for piglets fed ES compared to piglets fed Con (Table S1). Further metabolic pathway enrichment analysis demonstrated that piglets fed ES significantly altered their arginine biosynthesis, pyrimidine metabolism, primary bile acid biosynthesis, valine, leucine and isoleucine biosynthesis, alanine, aspartate and glutamate metabolism, linoleic acid metabolism, taurine and hypotaurine metabolism, pantothenate and CoA biosynthesis, and riboflavin metabolism (p<0.05, rich factor>0.1, Figure 6) compared to piglets fed Con.
Correlation analysis
Metabolites with VIP>1 (p<0.01) and genera with significantly different abundances (p<0.05) between piglets fed ES and Con were used for the Pearson’s correlation coefficient analysis (r > 0.4 or < -0.4, p<0.05). As is shown in Figure 7, the relative abundance of 3-Methoxy-4-Hydroxyphenylglycol Sulfate was positively correlated with Actinobacillus, Romboutsia, Terrisporobacter, and Veillonella, while negatively correlated with Buchnera. The relative abundance of N-Acetylneuraminic acid showed positive correlations with Terrisporobacter, Pseudomonas, Terrisporobacter, and Veillonella. It was also observed that Pro-Thr and Uridine levels were positively correlated with Bacteroides and Tyzzerella. Furthermore, cis-(6,9,12)-Linolenic acid was positively correlated with Romboutsia and Terrisporobacter. Thymine was positively correlated with Enterobacter. Buchnera was positively correlated with Betaine, Glycitein, and L-Citrulline, and negatively correlated with Muramic acid and Urea. Levels of TLR4 were positively correlated with the relative abundances of Buchnera, but was negatively correlated with the relative abundance of Terrisporobacter. Conversely, TNFα was positively correlated with the relative abundances of Romboutsia and Terrisporobacter, and was negatively correlated with the relative abundance of Buchnera and Enterobacter. Similarly, IL-8 was positively correlated with the relative abundances of Terrisporobacter, and was negatively correlated with the relative abundance of Buchnera and Enterobacter. It was also observed that ZO-1 was positively correlated with the relative abundances of Buchnera and Enterobacter, and negative correlations were observed with the relative abundance of Romboutsia, Terrisporobacter and Veillonella. Occludin was positively correlated with the relative abundances of Buchnera. Claudin-1 was positively correlated with the relative abundances of Buchnera and Enterobacter, although it was negatively correlated with the relative abundance of Romboutsia and Terrisporobacter. The ADG was positively correlated with the relative abundances of TLR4, ZO-1, and Claudin-1, and negative correlations were observed with the relative abundance of 3-Methoxy-4-Hydroxyphenylglycol Sulfate, L-Arginine, N-Acetylneuraminic acid, Pro-Ala, Uridine, IL-8, TNFα, and Veillonella.