Zacytuj

Introduction

The healthy creature intestine is home to microorganisms (Eckburg et al. 2005; Ley et al. 2006; Lozupone et al. 2012). Although microbiota resides in the intestines, it plays a critical role in the digestion and absorption of nutrients, maturation of immune system, anti-colonization, and stimulation of diverse host functions (Turnbaugh et al. 2006; Levy et al. 2017). The jejunum is a significant site for nutrient absorption (Martinez-Guryn et al. 2019). The jejunal microbiota is closely related to amino acid metabolism (Dai et al. 2010) and lipid deposition (Li et al. 2019).

Previous studies indicated that the host’s genetics shapes the microbial repertoire (Goodrich et al. 2014; Goodrich et al. 2016). It was discovered that the intestinal microbiota in exotic pig breeds varies from Chinese indigenous pig breeds (Yang et al. 2014). To further explore this observation, two pig breeds with different host genetics (Saba and Landrace) were selected as the subjects in this study. The Saba pig is an indigenous breed in Chuxiong of Yunnan Province, China, and it is on the list of National Conservation Program for Chinese Indigenous Livestock Germplasm. Saba pigs grow slow, but this breed is characteristic of a high propensity for meat quality, ability to adapt to the environment, and disease resistance (Jeong et al. 2014; Diao et al. 2019). In contrast, the Landrace breed was commercially selected over generations for rapid growth and enhanced carcass yield (Briggs 1983).

Before birth the intestine of newborns is believed to be free of microbes (Turnbaugh and Turnbaugh 2008). Due to contact with sows and exposure to the surrounding environment, a complex microbial community rapidly colonizes the newborn mammal (Frese et al. 2015). The balanced microbiological system (diverse intestinal microbes) is a significant hallmark of piglet health (Patil et al. 2019). Suckling piglets are an essential stage in the life of pigs, and thus more attention should be paid to the intestinal microbiota of piglets. The 35-day-old piglets easy to cause any diseases or dramatic internal environmental changes, they are about to wean; therefore, we selected 35-day-old Saba and Landrace suckling piglets as the subjects of this study. A comparison of their jejunal microbiota diversity will help comprehend the composition and functionality of gut microbiota in Chinese indigenous pigs.

Experimental
Materials and Methods

Animals and samples collection. All Saba and Landrace pigs were raised on a commercial farm in Chuxiong of Yunnan Province, China. Three Saba and three Landrace sows of third parity were selected for this study. They lived in six enclosures in an environmentally controlled room and were fed with the National Research Council (NRC) diet without antibiotics. After parturition, all piglets from every sow were placed in a single enclosure and fed by sow’s milk until day 35 (35 d). From each sow, three piglets were randomly selected, and their birth weight and 35 d body weight were recorded, and the average daily gain of both groups was calculated. Piglets were sacrificed, and the content from the middle of the jejunum was collected for 16S rRNA sequencing analysis.

DNA extraction and PCR amplification. Based on the manufacturer’s instructions, the QIAamp® Fast DNA Stool Mini Kit (Qiagen, Cat No.: 19593) was used to extracted Genomic DNA from 18 samples. The V3-V4 region of the bacterial 16S ribosomal RNA genes was amplified following the method of Fadrosh and coworkers (Fadrosh et al. 2014).

Illumina MiSeq PE250 sequencing. The Qubit® 2.0 (Invitrogen, USA) was used to quantify DNA in the samples for library preparation. During the amplification, the barcodes were introduced by the ligated primers, which included sequencing adaptor, barcode, and sequence binding to V3-V4 region. The libraries were sequenced on the MiSeq platform (Illumina, Inc., CA, USA). All jejunal content samples from 18 piglets were subjected to 16S rRNA sequencing; however, one sample from Saba piglets and two samples from Landrace piglets failed to build a database.

Processing of sequencing data. The sequencing data analysis referred to the method of Li and coworkers (Li et al. 2019). Trimming of barcodes and primers was performed using Pandaseq (https://github.com/neufeld/pandaseq/releases/tag/v2.8.1), followed by the quality control (e.g., the lengths of reads and an average base quality) using Fastqc (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). 16S rDNA tags between 220 bp and 500 bp, with no more than three ambiguous N, were kept, and the average Phred score of bases was no worse than 20 (Q20). The copy number of tags was enumerated, and the redundant tags were removed. Only the tags with a frequency higher than 1, which are more reliable in general, were clustered into Operational Taxonomic Units (OTUs). Each OTU had a representative tag. OTUs were clustered with a criterion of 97% similarity using the Uparse (http://drive5.com/uparse/), with chimeric sequences identified and removed using the Userach (version 7.0). Each representative tag was assigned to taxa by RDP Classifier (http://rdp.cme.msu.edu/) against the RDP database (http://rdp.cme.msu.edu/) using a confidence threshold of 0.8.

The OTU profiling table and alpha/beta diversity were also achieved by Python scripts of QIIME. Alpha diversity was the species diversity in each sample, including community abundance (Chao1 index), the diversity (Shannon and Simpson index), the phylogenetic diversity index (PD whole tree), and coverage (Good’s coverage values). QIIME software was used to calculate the samples’ alpha diversity index based on the OTU results and to generate the corresponding dilution curve. The Bray-Curtis distance was calculated to estimate the dissimilarity in the community structure, which was visualized using principal coordinates analysis (PCoA). Analysis of similarities (ANOSIM) was performed in the Mothur v1.380. We determined the strength of these groups using multiresponse permutation procedures (MRPP). Both analyses were performed in the PC-ORD. In addition to p-values, PC-ORD generated T and A values for all comparisons in the MRPP. T was a measure of separation between groups, with more negative values indicating a stronger separation. Group homogeneity was described by A and was scaled between 0 and 1.

The linear discriminant analysis (LDA) effect size (LEfSe) method (p < 0.05, LDA > 2) was used to identify the most differentially abundant OTUs between groups, with the LDA obtained by a pair-wise computation. The Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) based on a closed-reference operational taxonomic unit (OTU) was used to predict the abundances of functional categories in the Kyoto Encyclopedia of Genes and Genomes (KEGG) ortholog (KO). The correlation coefficients between KEGG pathways and bacterial compositions were calculated using Pearson’s correlation test in GraphPad Prism 7.

Statistical analysis. The experimental data, including growth performance and microbiota abundances, were analyzed with the SPSS 22.0 software. Grade and quantitative data were compared with the t-test between the two groups.

We used Spearman’s test to estimate the correlation between KEGG pathway and jejunal microbial composition and host growth performance. P < 0.05 was deemed to statistical significance.

Results

Growth performance of Saba and Landrace piglets. The growth performance of Saba and Landrace piglets is shown in Table I. The birth weight of Landrace piglets was significantly higher than Saba piglets (p < 0.001). On day 35, the body weight and average daily gain (p < 0.001) of Landrace piglets were higher than Saba piglets (p < 0.001).

Growth performance of Saba and Landrace piglets.

SBLA
Birth weight (kg)0.76 ± 0.20B1.99 ± 0.14A
Body weight (kg)4.69 ± 1.14B10.22 ± 0.57A
Average daily gain (kg)0.11 ± 0.03B0.24 ± 0.02A
Stem length (cm)37.00 ± 4.84B50.00 ± 2.27A
Height at withers (cm)22.75 ± 1.83B29.00 ± 1.6A
Chest measurement (cm)37.88 ± 3.48B48.38 ± 1.41A
Chest depth (cm)10.75 ± 0.71B14.38 ± 2.07A
Abdominal girth (cm)38.75 ± 3.85B49.50 ± 1.60A
Cannon circumference (cm)8.06 ± 0.56B10.13 ± 0.35A

Different superscripts in the same column indicate significant difference (p < 0.001)

SB – Saba piglets, LA – Landrace piglets

Gut microbiota DNA sequence data and quality control. Sequencing of the amplicons of the 16S rRNA gene at MiSeq generated 884,982 clean reads (mean length of 415 bp) with 368,167,415 base pairs in total, yielding an average of 58,999 clean reads (55,547– 63,409), and 24,544,494 base pairs (22,667,528 bp – 26,990,140 bp) per sample (Table II). Out of the high-quality sequences, about 99.53% were between 420 and 460 bp for these two breeds.

Description of the assembly results of jejunum microbiota from piglets.

Sample nameClean ReadsBases (bp)Q20 (%)Q30 (%)GC (%)Average length (bp)
LA-15827824288798    0.96020.88730.5229416
LA-25597523383577    0.95910.88550.5188417
LA-36239625557801    0.96960.90810.527409
LA-45656523778906    0.96140.89220.5365420
LA-55554722929087    0.96710.90540.5221412
LA-65633623469885    0.96120.89280.5234416
LA-757139237718480.960.89020.5332416
SB-16339326990140    0.96090.89170.5154425
SB-25595922667528    0.96970.91190.5271405
SB-36248425749078    0.96610.90350.5185412
SB-46152726100871    0.95940.88670.5117424
SB-55690724136231    0.96020.88970.5123424
SB-66340926890638  0.9580.88480.5496424
SB-75892023815087    0.96560.90290.5306404
SB-86014724637940    0.96260.89660.5253409

SB – Saba piglets, LA – Landrace piglets

Diversity in jejunal microbiota of Saba and Landrace piglets. We revealed that the jejunal microbiome was different in Saba and Landrace piglets. The USEARCH algorithm was used to cluster at a 0.97 similarity level, and the clustered sequences were filtered by a chimera. We obtained 489 and 365 OTUs (Fig. 1) from Saba and Landrace piglets, respectively. In total, 254 OTUs were shared by Saba and Landrace piglets. The alpha diversity index of the samples is shown in Table III. The PD whole tree in Saba piglets (17.76) was significantly higher (p < 0.05) than in Landrace piglets (13.31). The Chao1 index (242.85 vs. 229.04), the observed species index (139.75 vs. 124.75), and Shannon (3.06 vs. 2.84) and Simpson (0.79 vs. 0.76) indexes for microbiota from Saba piglets are higher than those from Landrace piglets, but the statistical significance (p > 0.05) was not noticed. Besides, we compared the beta diversities between all the samples (Table IV). The test statistic (R) of Multi Response Permutation Procedure (MRPP) was 0.031 (p = 0.102) on the weighted UniFrac, and 0.034 (p = 0.005) on the unweighted Uni-Frac. Also, the test statistic (R) of ANOSIM was 0.131 (p = 0.01; Fig. 2) on the unweighted UniFrac. Using the unweighted UniFrac metric, the Principal Coordinates Analysis (PCoA) showed a clear separation between Saba and Landrace piglet samples (Fig. 2).

Fig. 1.

Venn diagram of OTUs clustered at 97% sequence identity of microbiotas from Saba and Landrace piglets. The number of overlapping parts is the total number of OTUs between the groups, while the numbers in non-overlapping parts indicate the number of unique OTUs for each group. SB – Saba piglets, LA – Landrace piglets.

Alpha diversity in jejunal microbiota between Saba and Landrace piglets.

LASBp value
Chao1 index229.04 ± 38.23242.85 ± 10.920.531
The observed_species index124.75 ± 24.27139.75 ± 36.570.520
PD_whole_tree13.31 ± 1.9817.76 ± 2.170.023
Shannon index2.84 ± 0.423.06 ± 0.890.671
Simpson index0.76 ± 0.0770.79 ± 0.0960.665
Goods_coverage0.998 ± 0.000150.999 ± 0.000240.149

SB – Saba piglets, LA – Landrace piglets

MRPP of the 16S rRNA gene between Saba and Landrace piglets.

AObserve DeltaExpect DeltaSignificance
The weighted_unifrac0.03126998502417340.2766829870044870.285614136784429  0.102
The unweighted_unifrac0.03383530543520210.5515816682596760.5708981826417620.00

SB – Saba piglets, LA – Landrace piglets

Fig. 2.

Principal coordinate analysis (PCoA) illustrated bacterial community structures based on Bray-Curtis distances. On the PCoA plot, each color represents one group. Unweighted and weighted PCoA of β-diversity measures of all samples. PCOA1 (19.67%) and PCOA2 (13.63%).

Comparison of jejunal microbiota of Saba and Landrace piglets. The jejunal bacterial taxa were diversified between Saba and Landrace piglets at the phylum level (Fig. 3A). Among these taxa, Firmicutes occurred with the highest abundance within the jejunal microbiota of Landrace piglets (95.82%), followed by Proteobacteria (1.67%), and Bacteroidetes (0.052%). Similarly, the higher relative abundance of Firmicutes (97.6%, p > 0.05) than Bacteroidetes (1.14%, p < 0.05), and a lower proportion of Proteobacteria (0.71%, p > 0.05) were observed in microbiota of Saba piglets.

Fig. 3.

Community composition of the jejunum microbial of Saba and Landrace piglets at the phylum (A) and genus (B) levels, respectively. Data are expressed as means + MSE, *p < 0.05.

The abundance of bacterial species within jejunal taxa is shown in Fig. 3B. The two most abundant genera were Lactobacilli and Clostridium XI, accounting for 28.0% and 42.24% of the jejunal species in Landrace piglets, respectively. Compared with the microbiota of Landrace piglets, the jejunal microbial of Saba piglets had a higher abundance of Lactobacilli (36.81%) and a lower abundance of Clostridium XI (40.02%), but the difference was not statistically significant (p > 0.05). Moreover, Veillonella (0.58% vs. 2.34%), Streptococcus (0.23% vs. 1.32%), and Saccharibacteria genera incertae sedis (0.19% vs. 1.04%) of Saba piglets were remarkably lower than in microbiota of Landrace piglets (p < 0.05).

Differences of bacterial taxa between Saba and Landrace piglets. The different number of OTUs was observed between the jejunal microbiota of Saba and Landrace piglets (Fig. 4). There was one main phylum (Firmicutes) and two genera (Coprococcus and Parabacteroides) significantly enriched in jejunal microbiota of Saba piglets. Also, multiple biomarkers were significantly enriched in jejunal microbiota of Landrace piglets, including two phyla (Candidatus Saccharibacteria and Proteobacteria), two classes (Epsilonproteobacteria and Gammaproteobacteria), two orders (Campylobacterales and Pasteurellales), eight families (Fusobacteriaceae, Leuconostocaceae, Actinomycetaceae, Enterococcaceae, Campylobacteraceae, Dermatophilaceae, Streptococcaceae, and Pasteurellaceae), thirteen genera (Enterococcus, Actinomyces, Fusobacterium, Weissella, Pediococcus, Campylobacter, Oribacterium, Sharpea, Tonsilliphilus, Pasteurella, Saccharibacteria genera incertae sedis, Streptococcus, and Actinobacillus). Furthermore, the increase in the abundance of the phylum Candidatus Saccharibacteria was represented by an increased abundance of the genus Saccharibacteria genera incertae sedis (Fig. 5).

Fig. 4.

Alteration of the relative abundance of bacteria in the Saba and Landrace piglets using linear discriminant analysis effect size (LEfSe). Each bar represents the log 10 effect size (LDA score) for a specific taxon. A longer bar represents a higher LDA score. Only taxa meeting an LDA significant threshold of 2 are shown. These taxa showed a statistically significant difference between the Saba and Landrace piglets (p < 0.05 by the Wilcoxon test); each color represents one group.

p – phylum, c – class, o – order, f – family, and g – genus.

Fig. 5.

A cladogram showed a comparison of the bacterial microbial profiles from Saba and Landrace piglets.

p – phylum, c – class, o – order, f – family, and g – genus.

Correlation between microbiota and growth performance. The correlation between jejunal microbiota and host growth performance was shown in Fig. 6. The Coprococcus was negatively correlated with the body weight (p = 0.046, R = –0.53), average daily gain (p = 0.046, R = –0.53), stem length (p = 0.018, R = –0.61), height at withers (p = 0.0063, R = –0.68), chest measurement (p = 0.046, R = –0.53), chest depth (p = 0.012, R = –0.64), abdominal girth (p = 0.045, R = –0.53), and cannon circumference (p = 0.011, R = –0.65). The Parabacteroides was negatively correlated with the cannon circumference (p = 0.048, R = –0.52). The Tonsilliphilus was positively correlated with the body weight (p = 0.025, R = 0.58), average daily gain (p = 0.022, R = 0.59), stem length (p = 0.026, R = 0.58), height at withers (p = 0.013, R = 0.64), chest measurement (p = 0.0031, R = 0.72), chest depth (p = 0.0015, R = 0.74), abdominal girth (p = 0.0053, R = 0.69), and cannon circumference (p = 0.022, R = 0.72). The Saccharibacteria genera incertae sedis was positively correlated with the stem length (p = 0.036, R = 0.55), chest measurement (p = 0.038, R = 0.54), and abdominal girth (p = 0.031, R = 0.56). The Enterococcus was positively correlated with the body weight (p = 0.0097, R = 0.65), average daily gain (p = 0.015, R = 0.62), stem length (p = 0.049, R = 0.52), height at withers (p = 0.0048, R = 0.70), chest measurement (p = 0.0066, R = 0.68), chest depth (p = 0.0069, R = 0.68), abdominal girth (p = 0.0070, R = 0.67), and cannon circumference (p = 0.0052, R = 0.71). The Pediococcus was positively correlated with the body weight (p = 0.023, R = 0.59), average daily gain (p = 0.021, R = 0.59), stem length (p = 0.026, R = 0.58), height at withers (p = 0.013, R = 0.64), chest measurement (p = 0.0025, R = 0.72), chest depth (p = 0.0021, R = 0.74), abdominal girth (p = 0.0048, R = 0.69), and cannon circumference (p = 0.018, R = 0.63). The Weissella was positively correlated with the body weight (p = 0.033, R = 0.56), average daily gain (p = 0.033, R = 0.56), stem length (p = 0.038, R = 0.54), height at withers (p = 0.011, R = 0.66), chest measurement (p = 0.0039, R = 0.70), chest depth (p = 0.0029, R = 0.73), abdominal girth (p = 0.0078, R = 0.66), and cannon circumference (p = 0.022, R = 0.62). The Streptococcus was positively correlated with the body weight (p = 0.0024, R = 0.74), average daily gain (p = 0.0020, R = 0.75), stem length (p = 0.0037, R = 0.71), height at withers (p = 0.0073, R = 0.67), chest measurement (p = 0.0011, R = 0.77), chest depth (p = 0.0099, R = 0.65), abdominal girth (p = 0.00097, R = 0.78), and cannon circumference (p = 0.0038, R = 0.72). The Sharpea was positively correlated with the body weight (p = 0.023, R = 0.59), average daily gain (p = 0.023, R = 0.59), height at withers (p = 0.0060, R = 0.69), chest measurement (p = 0.0096, R = 0.66), chest depth (p = 0.0079, R = 0.67), abdominal girth (p = 0.010, R = 0.65), and cannon circumference (p = 0.0099, R = 0.66). The Campylobacter was positively correlated with the body weight (p = 0.0034, R = 0.72), average daily gain (p = 0.0031, R = 0.73), stem length (p = 0.0012, R = 0.78), height at withers (p = 0.011, R = 0.65), chest measurement (p = 0.0022, R = 0.75), chest depth (p = 0.0012, R = 0.78), abdominal girth (p = 0.0013, R = 0.77), and cannon circumference (p = 0.0066, R = 0.69).

Fig. 6.

Heatmap analysis of the correlation between microbiota and growth performance.

KEGG pathway and their correlation with microbiota. To assess the jejunal microbiota’s metabolic potential, we performed Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) (Langille et al. 2013; Javurek et al. 2016). KEGG pathway (L3 hierarchy) analysis is shown in Fig. 7. The “d-glutamine and d-glutamate metabolism”, and “one carbon pool by folate” pathway were enriched in Saba piglets. The “limonene and pinene degradation”, “tryptophan metabolism”, and “sulfur relay system” were enriched in Landrace piglets.

Fig. 7.

KEGG enrichment analysis of the difference within groups at the L3 hierarchy.

We used Spearman’s correlation heatmap (Fig. 8) to study the correlation between the jejunal microbiota and the KEGG pathway. The “d-glutamine and d-glutamate metabolism” pathway was positively correlated with the presence of Firmicutes (family, p = 0.0080, R = 0.67), while negatively correlated with Fusobacteriaceae (family, p = 0.020, R = –0.60), and Fusobacterium (genus, p = 0.017, R = –0.62). The pathway of “one carbon pool by folate” was negatively correlated with Fusobacteriaceae (family, p = 0.049, R = –0.52) and Fusobacterium (genus, p = 0.048, R = –0.52). The pathway “tryptophan metabolism” was negatively correlated with Coprococcus (genus, p = 0.023, R = –0.59), but positively correlated with Gammaproteobacteria (class, p = 0.024, R = 0.59), Enterococcaceae (family, p = 0.0089, R = 0.66), Proteobacteria (family, p = 0.022, R = 0.59), and Enterococcus (genus, p = 0.031, R = 0.56).

Fig. 8.

Pearson’s correlation analysis of microorganisms and signal pathways in Saba and Landrace piglets. Heatmap analysis of the correlation between microorganisms and signal pathways. Correlations with p < 0.05 are shown. Blue represents a significant negative correlation (p < 0.05), red represents a significant positive correlation (p < 0.05), and white represents no significant correlation (p > 0.05). The number represents the value of R (p < 0.05).

Discussion

The Landrace pig from Denmark is a typical commercial pig breed of fast growth and high carcass yield (Briggs 1983). The previous research reported that the body weight of Landrace piglets was 1.68 kg and 6.52 kg on day 1 and day 27, respectively (Li et al. 2013). In contrast, the Saba pig is an indigenous breed from China, with a relatively slow growth rate. In our study, the birth weight, body weight (day 35), and average daily gain of Landrace piglets (1.99 kg, 10.22 kg, and 0.24 kg/d, respectively) were higher than those of Saba piglets (0.76 kg, 4.69 kg, and 0.11 kg/d, respectively). The data indicated that the growth performance of Landrace piglets was higher than Saba piglets. Previous studies have shown that growth performance and intestinal microbes were different in Jinhua pigs and Landrace pigs of the same age (Xiao et al. 2018). In addition, our results show that the Coprococcus and Parabacteroides were negatively correlated with the growth performance, while the Tonsilliphilus, Enterococcus, Pediococcus, Weissella, Streptococcus, Campylobacter, Saccharibacteria genera incertae sedis, and Sharpea were positively associated with the growth performance. The above results suggested that the composition of intestinal microbiota was significantly and closely connected with the pig breed.

It is generally believed that intestinal microorganisms have abundant metabolic profiles to maintain their basic life and have a considerable impact on host growth and health (Turnbaugh et al. 2006). The accumulating evidence suggested that diet (Pluske 2013), environment (Thompson et al. 2008), and host’s genetics (Büsing and Zeyner 2015; Hancox et al. 2015) can affect the composition of intestinal microbiota. The previous research (Yang et al. 2014) revealed that the percentages of Firmicutes and Bacteroidetes in the Chinese indigenous pig breeds (Xiaomeishan, Meishan, and Bama sows) were higher than those of exotic breeds (Landrace, Yorkshire, and Duroc sows). It is consistent with our finding that the Firmicutes, Proteobacteria, and Bacteroidetes dominated in the jejunum of both pig breeds. Furthermore, the percentage of Bacteroidetes in Saba piglets was significantly higher than in Landrace piglets. Saba piglets are obese, and Landrace piglets are lean. A previous study has shown that fat deposition is positively correlated with the presence of Bacteroidetes and Firmicutes within the intestinal microbiota (Turnbaugh et al. 2006). Nevertheless, the mechanisms between intestinal microbiota and fat deposition are still unclear, and further study is needed.

At the genus level, the two most numerous genera in the pig’s small intestine were Lactobacillus and Clostridium (Crespo-Piazuelo et al. 2018). Our study demonstrated that Clostridium XI is the most numerous bacteria in jejune of both Saba and Landrace piglets. We speculate that this was principally due to the adaptation of the microbial system for milk nutrition, as sucking lambs have a large proportion of Clostridium XI in the intestinal microbial community (Bi et al. 2019). The Lactobacilli were usually considered beneficial bacteria responsible for more effective anti-inflammation and out-competing microbiota competences (Etzold et al. 2014). Therefore, some Lactobacillus species have been used as substitutes for antibiotics for growth promotion. The average daily gain (ADG) of weaned piglets fed with Lactobacillus plantarum PFM 105 was significantly improved after three weeks (Wang et al. 2019). Again, two Lactobacillus strains (Lactobacillus frumenti and Lactobacillus gasseri LA39), when taken orally, could significantly prevent stress-induced diarrhea caused by early weaning of piglets (Hu et al. 2018). In line with that, in this study, the abundance of Lactobacillus in Saba piglets was higher than in Landrace piglets.

Furthermore, some taxa are recognized as negatively correlated with host health. Veillonella is present in piglets infected with porcine epidemic diarrhea virus (PEDV), caused by the disorder of intestinal amino acid metabolism and energy metabolism (Huang et al. 2018). Also, Veillonella and Streptococcus can increase the secretion of inflammatory cytokine and the production of anti-microbial peptides (AMP), resulting in improved mucosal thickness and epithelial barrier function in 3D-reconstructed human gingiva (Shang et al. 2018). In the present study, a higher proportion of Lactobacilli and a lower proportion of Streptococcus and Veillonella in the microbiome of Saba piglets, suggested that Saba piglets might have stronger disease resistance. Overall, our data indicated that the pig breed could influence the microbiome composition.

The gut microbiome is beneficial for pigs, contributing to improved vitamin K production, cellulose fermentation, and increased resistance to pathogens (Kim and Isaacson 2015; Stokes 2017; Yang et al. 2017). PICRUSt analysis of jejunal microbiotas’ metabolic potential showed that metabolic pathways were significantly different in Saba and Landrace piglets. It is noteworthy that the different amino acid metabolism pathways were enriched in Saba and Landrace piglets. The “d-glutamine and d-glutamate metabolism” pathway was enriched in Saba piglets, while “tryptophan metabolism” pathway was enriched in Landrace piglets. The “d-glutamine and d-glutamate metabolism” pathway participates in the inhibition of lipid peroxidation and quenches free radicals during oxidative stress (Qu et al. 2020). It has been reported that tryptophan is related to the immune response regulation, inflammation, and oxidative stress (Anesi et al. 2019; Liu et al. 2019). Therefore, different pig breeds might have regulated their health through different metabolic pathways.

The central part of amino acids absorption is the small intestine (Wu 1998). The catabolism of arginine and lysine in the jejunum could exceed their transport into intestinal cells (Dai et al. 2010). This phenomenon may be due to the role of intestinal microorganisms. Besides, owing to the deficiency of several key enzymes, the threonine, tryptophan, histidine, lysine, and methionine cannot be metabolized by porcine intestinal cells in the presence of amino acids at physiological concentrations (Chen et al. 2007). However, the histidine, glutamate, threonine, and lysine were utilized by microbiota of the porcine small intestine (Dai et al. 2010). The above results indicated that jejunal microorganisms could participate in amino acid metabolism. In our study, the Firmicutes enriched in the jejunal microbiome of Saba piglets were positively associated with the “d-glutamine and d-glutamate metabolism” pathway. In contrast Fusobacteriaceae and Fusobacterium that enriched the jejunal microbiome of Landrace pigs were negatively correlated with the above pathway. Finally, Coprococcus from the jejunal microbiome of Saba piglets was negatively correlated with the “tryptophan metabolism” pathway, and Enterococcaceae, Enterococcus, Proteobacteria, and Gammaproteobacteria from Landrace piglets were positively correlated with tryptophan metabolism. Therefore, these taxa may be related to the metabolism of amino acids in jejunum.

Conclusions

In summary, the growth performance was higher for Landrace piglets compared to Saba piglets due to their different genetic characteristics. The rich diversity and fewer infection-associated taxa were observed in Saba piglets, partially accounting for their high adaptability to environmental perturbations compared to Landrace piglets. Several taxa in the jejunum of Saba and Landrace piglets were associated with “d-glutamine and d-glutamate metabolism” and “tryptophan metabolism” respectively, suggesting that pig breeds may regulate their health through different metabolic pathways. Although the interaction between pig and microbiota needs further extensive investigations, our study would shed more light on the functional exploration and resource development of local pig intestinal microbiota in China.

eISSN:
2544-4646
Język:
Angielski
Częstotliwość wydawania:
4 razy w roku
Dziedziny czasopisma:
Life Sciences, Microbiology and Virology