The growth performance does not statistically differ between DLY and NX piglets
To compare the growth performance between NX and DLY piglets, we measured the average daily feed intake (ADFI), average daily gain (ADG), feed conversion ratio (FCR) and organ index in the trial. As shown in Fig S1A, there was no significantly difference in ADFI between NX and DLY piglets, despite the DLY piglets exhibited the trend of higher ADFI than NX piglets. Regarding the ADG between DLY and NX piglets, DLY piglets showed a significantly increased ADG than NX piglets in the last week of the trial (Fig S1B). While the FCR and the organ index (organ weight/ body weight) between DLY and NX piglets did not show statistically significant p value (Fig S1C and S1D).
Pig breeds exhibit distinct colonic microbial characteristics
To compare the structure and function of gut microorganisms between NX and DLY piglets, we performed metagenomic analyses in the colon of these two breeds. Alpha-diversity analyses indicated that the DLY piglets consist of the higher total number of gut microorganisms in the colon when compared to NX piglets (Fig. 1A). The species diversity of gut microorganisms did not significantly differ between DLY and NX piglets which were measured by using Simpson and Shannon indices. While NX piglets manifested the trend of higher gut microbial species and richness than DLY piglets in the colon (Fig. 1A). There was a significantly difference between DLY and NX piglets in compositional similarity measured by Beta-diversity as shown in Fig. 1B.
To further uncover the distinct gut microbial composition between DLY and NX piglets, we analyzed the taxonomic composition of the colonic microbiota at different levels. DLY piglets had much more archaea distribution in the colon than NX piglets at the Kingdom level (Fig S2A and S3A). The phylum-level bacteria analyses found that the abundance of Bacteroidetes was greater in NX piglets than DLY piglets, while DLY piglets showed an increased abundance of Spirochaetes, Proteobacteria, Euryarchaeota and Fimicutes in the colon when compared to NX piglets (Fig S2B and S3B). At the Class level, NX piglets displayed a higher abundance of colonic Clostridia, Erysipelotrichia and Bacteroidia than DLY piglets (Fig S2C and S3C). Furthermore, the abundance of Clostridiales, Bacillales, Erysipelotrichales, Bacteroidales, Flavobacteriales and Enterobacterales were significantly increased in the NX piglets than DLY piglets at the order level (Fig S2D and S3D). Similarly, the increased abundance of Lactobacillaceae, Bacteroidaceae, Lachnospiraceae, Prevotellaceae, Ruminococcaceae, Clostridiaceae, Peptostreptococcaceae, Enterobacteriaceae, Odoribacteraceae, Rikenellaceae and Erysipelotrichaeae were observed in the colon of NX piglets at the Family level (Fig. S2E and S3E). Regarding the differences in the colonic microbiota between DLY and NX piglets at the Genus level, we found that the abundance of Bacteroides, Roseburia, Prevotella, Phocaeicola, Anaerrobutyricum, Lacobacillus were significantly enriched in NX piglets than that in DLY piglets (Fig. 1C). In a sharp contrast, the abundance of Megasphaera, Streptococcus and Treponema were significantly increased in DLY piglets than that in NX piglets.
NX and DLY piglets colonize different microbial species in the colon
Next, we discriminated the differentially colonized microbes (DCMs) at the species level between DLY and NX piglets, and identified 665 of gut microbial species that are specifically enriched in DLY piglets, and 371 of gut microbial species that are specifically enriched in NX piglets (Fig. 1D and Fig. 2A). Then we sought to compare the dissimilarity of colonic bacteria, archaea and fungi between these two breeds, respectively. Firstly, we evaluated the top 30 DCMs between DLY and NX piglets and found that the abundance of Clostridium species (C. butyricum, C. botulinum and C. perfringens), Bacteroides species (B. thetaiotaomicron and B. fragilis), Prevotella species (P. enoeca, P. intermedia and P. dentalis) and Lactobacillus_crispatus, are dominantly enriched in NX piglets. In contrast, DLY piglets showed the high level of colonic Streptococcus species (S. lutetiensis, S. infantarius, S. equninus, S. pasteurianus and S. lutetiensis) and Selenomonas_sputigena (Fig. 2B). Next, we also noticed the varied colonization of archaea between DLY and NX piglets. The abundance of Pyrodictium_delaneyi and Methanohalopilus_halophilus were significantly increased in NX piglet than DLY piglet, while the Halopenitus_persicus, Halolamina_sp_CBA1230, Halorubrum_sp.PV6 Halarchaeum_sp.CBA1220 and Methanosphaera_stadtmanae were preferentially colonized in the colon of DLY piglets (Fig. 2C). Lastly the abundance of some fungal species like Tetrapisispora_blattae and Zygosaccharomyces_rouxii were significantly increased in the colon of NX piglets than DLY piglets (Fig. 2D).
To further identify specific colonic bacterial genera that may act as the unique characteristics of these two breeds, LEFSe was used to evaluate the bacterial composition between DLY and NX piglets. Our analyses found that the species of Streptococcus (S. lutetiensis, S. pasteurianus, S. thermophilus, S. gallolyticus and S. infantarius) and Selenomonas_sputigens were specifically enriched in DLY piglets when compared with NX piglets. While the species of Bacteroides (B. fragilis, B. thetaiotaomicron, B. intestinalis and B. sp_PHL2737), Prevotella (P. dentalis, P. ruminicola, P. oris) were exclusively colonialized in NX piglets (Fig. 2E). Together, our analyses suggest that DLY and NX piglets display the vast differences in the colonization of gut microbial species.
NX piglets-enriched gut microbes modulate several metabolic functions
To further examine the function of the gut microbiome between these two breeds, we first assessed DEMs by conducting Kyoto Encyclopedia of Genes and Genomes (KEGG), Cazy, BacMet and Swissprot analyses. As shown in Fig. 3A-3D, all four analyses displayed statistically significant p value (p < 0.05) between NX and DLY piglets. KEGG data showed that the functional categories related to several metabolic pathways such as methane metabolism, fructose and mannose metabolism, glycerolipid metabolism, selenocompound metabolism, purines metabolism and D-amino acid metabolism were significantly enriched in the colon of NX piglets. While the function of DLY piglets-enriched gut microbes were dominant in phenylalanine metabolism, pentose and glucuronate interconversions, and oxidative phosphorylation pathway (Fig. 3E).
Then the Cazy analysis identified a subset of carbohydrate-active enzymes including glycoside hydrolases (GHs) and polysaccharide lyases (PLs) that differ between DLY and NX piglets. There were 20 upregulated GHs and PLs in NX piglets when compared to DLY piglets, among which GH43, GH43_10, PL1, PL1_2 and carbohydrate esterases 12 (CE12) represented the most significantly upregulated ones. By contrast, GH25, GH13_31, GH13_26, GH13_3 and PL0 were slightly up-regulated in DLY piglets than NX piglets (Fig. 3F).
Next, we performed LDA analyses to further unravel the correlations between KEGG and DEMs. Our results showed that the Prevotella species (P.dentalls, P.oris, P.ruminicola) and Bacteroides species (B.fragilis, B.intestinalis, B.sp.PHL_2737 ) were positively correlated with several nutritional metabolic processes, such as amino acid metabolism, energy metabolism, carbohydrate and lipid metabolism. While the Streptococcus species and Acidaminocossus_fementans were positively correlated with xenobiotics biodegeneration and metabolism and nucleotide metabolism (Fig. 3G). Similarly, we also determined the correlations between carbohydrate-active enzymes and DEMs, and found that the Streptococcus species (S. equinus, S. infantanius, S. spuligena, S. pasteurianus) were positively correlated with GT112, GH4 and GH13_26 enzymes. And the activity of GH43, GH43_10, GH43_19, GH127, GH146, PL1 and CE0 were significantly correlated with the Prevotella species (P. dentails, P. ruminiccla) and Bacteroides species (B. thetaiotaomiccron, B. fragilis, B. intestinalis) (Fig. 3H).
Gut microbiota correlates with growth performance both in DLY and NX piglets
The gut microbiota has the ability to modulate nutrition level and physical state, affecting the growth performance of the host. Therefore, we sought to determine whether the gut microbiota correlates with ADG and ADFI in DLY and NX piglets by using Spearman correlation coefficient and Pearson correlation coefficient analyses. Our data found that the Bacteroides species (B.xylanisolvens, B.fragilis) and Clostridium species ( C.butyricum and C.botulinum) were negatively correlated with ADFI and ADG (Fig. 4A-4B), suggesting that the high abundance of Bacterides and Clostridium species may decrease the ADFI and ADG in DLY and NX piglets. Then, we further analyzed the correlations between the gut microbial function and ADFI and ADG. As shown in Fig. 4C-4F, despite the low correlations between the KEGG and ADFI and ADG, or the Cazy and ADFI and ADG were observed, there were still some small categories showing statistically significant p vulues. For instance, the small subunit ribosomal protein S12 (K02950), putative transposase (K07491), large subunit ribosomal protein L16 (K02878), ABC-2 type transport system permease protein (K01992) were negatively correlated with ADG and ADFI. While the ferrous iron transport protein A (K04758), transposase, IS30 family (K07482), large subunit ribosomal protein L33 (K02913), DNA-damage-inducible protein J (K07473) and small subunit ribosomal protein S16 (K02959) were positively correlated with ADG and ADFI (Fig. 4D). In addition, the PLs, such as PL1 and PL1_2, were negatively correlated with ADG, but not ADFI. The group of GH43_10, GH35, GH43, GT14 and GH127 showed negatively correlation with ADFI and ADG. But the PL0, GH29, GH20, GH29 and GH92 were positively correlated with ADG and ADFI, respectively (Fig. 4F). Lastly, the BacMet (a database for antibacterial biocide and metal resistance genes) and Card (the antibiotic resistance database) analyses showed that copper, hydrogen peroxide, penam and lincosamide antibiotic were positively correlated with ADFI, and the manganese was negatively correlated with ADG and ADFI, respectively (Fig. S4A-4D).
Intestinal genes expression pattern differs between DLY and NX piglets
To reveal the global mRNA gene expression pattern in these two breeds, we conducted mRNA sequencing in the colon of DLY and NX piglets. We identified 433 differentially expressed genes (DEGs), with 205 upregulated and 228 downregulated genes in NX piglets than DLY piglets (Fig. S5A). KEGG pathway analyses showed that DEGs are mainly involved in regulating cytokine-cytokine receptor interaction, intestinal immune activity, neuroactive ligand-receptor interaction, and chemokine signaling pathway (Fig S5B). Gene ontology analyses according to “Biological process “of significantly up- or downregulated genes revealed the enrichment of genes involved in carbohydrate binding activity, C-C chemokine activity and cytokine binding and transmembrane transporter activity (Fig S5C).
NX piglets exhibit improved intestinal barrier function than DLY piglets
Intestinal epithelium barrier not only contributes to the gut development but also acts as the core immune system to defend against microbial infections. Then, we compared the colonic barrier function between DLY and NX piglets. There were 30 DEGs related to intestinal barrier function that have been identified in our mRNA sequencing data, with 17 upregulated and 13 downregulated genes in NX piglets when compared to DLY piglets (Fig. 5A). As expected, qPCR experiment confirmed the significantly upregulation of a subset of intestinal barrier related genes like E-cadeherin, WAS (WASP actin nucleation promoting factor), TAD1 (tRNA-specific adenosine deaminase) and HCLS (hematopoietic cell-specific Lyn substrate 1) in the colon of NX piglets, in consistent with our mRNA sequencing data (Fig. 5A). These results indicate that NX piglets have better intestinal barrier function than DLY piglet.
Next, GO and KEGG analyses were further performed to examine the potential function of differentially expressed intestinal barrier genes, respectively. GO data showed that these upregulated genes are enriched in the binding activity of toll-like receptor and lipoteichoic acid receptor, oxidized low-density lipoprotein particle and short-chain fatty acid transmembrane transporter activity (Fig. 5B). While the group of downregulated genes showed dominant binding activity in cell adhesion molecule and protein domain specific factors (Fig. 5D). Moreover, KEGG analyses revealed the enrichment of upregulated genes in tight junction formation, fat digestion and absorption, PPAR (peroxisome proliferators-activated receptors) signaling pathway and cholesterol metabolism and lipid metabolism (Fig. 5C). Downregulated genes were involved in the formation of tight junction and defense against pathogenic Escherichia coli infection (Fig. 5E).
Varied intestinal immune reactions between NX and DLY piglets
Indigenous pig breeds have been shown to be associated with higher capacity of disease resistance. We first compared the expression of antimicrobial peptide genes (AMPs) in NX and DLY piglets. There were 22 differentially expressed AMPs, with 13 upregulated and 9 downregulated genes in NX piglets (Fig. 6A). Furthermore, RT-qPCR confirmed the significantly upregulated expression of regenerating islet-derived protein Ⅲ-γ (Reg-Ⅲ-γ) as well as two AMPs regulators Toll like receptor 1 (TLR1) and myeloid differentiation primary response gene 88 (Myd88)in NX piglets. While the mRNA expressions of beta-defensin 1 (PBD-1) and proline-arginine-rich (PR39) were drastically down-regulated in the colon of NX piglets when compared to DLY piglets (Fig. 6C). Our data suggest that these two breeds have different intestinal AMPs expression pattern and may exert distinct intestinal immune function. Additionally, we also examined the expression of inflammatory genes between NX and DLY piglets and identified 29 differentially expressed genes, with 20 upregulated and 9 downregulated inflammatory genes in the colon of NX piglets (Fig. 6B). Similarly, RT-qPCR validated the increased mRNA expressions of colonic interleukin-2 (IL-2), IL-17F, IL-18 and decreased mRNA expression of IL-8 and IL-22 in NX piglets when compared to DLY piglets (Fig. 6D). By conducting GO and KEGG analyses, we found that these differentially expressed immune responsive genes are mainly involved in modulating cytokine and chemokine activities, CCR chemokine receptor binding activity, and immune responsive signaling pathways including NF-kappaB inhibitor alpha (NF- κB) and TLRs signaling pathways (Fig. 6E-6H).