Skip to main content
  • Research article
  • Open access
  • Published:

Dietary inclusion of Peptiva, a peptide-based feed additive, can accelerate the maturation of the fecal bacterial microbiome in weaned pigs

Abstract

Background

Weaning is one of the most critical transition stages of the swine production cycle, as the piglet gut physiology and microbiome need to rapidly adapt to changes in diet and environmental conditions. Based on their potential for producing a vast array of bioactive molecules, peptide formulations represent a largely untapped source of compounds that could be developed into feed additives to benefit animal health and nutrition. In this context, a commercial-scale nursery trial was performed to evaluate the impact of low inclusion of a peptide-based feed additive (Peptiva, Vitech Bio-Chem Corporation) on the performance and fecal microbiome of weaned pigs.

Results

While no significant differences in body weight, daily gain, daily feed intake nor gain:feed were observed between control and treatment animals (P > 0.05), an effect of Peptiva on the fecal bacterial composition of weaned pigs was observed. The first main observation was that the fecal bacterial profiles from pigs fed Control-Phase II and Control Phase III diets were found to be very distinct, suggesting that a transition or succession stage had occurred between the two phases. Lactobacilli, represented by four main OTUs (Ssd-00002, Ssd-00019, Ssd-00025, and Ssd-00053), were more abundant at the end of Phase II (P < 0.05), while Streptococci, mostly represented by OTUs Ssd-00039 and Ssd-00048, were in higher abundance at the end of Phase III (P < 0.05). Secondly, the fecal bacterial composition from pigs fed Peptiva Phase II diets showed similarities to both Control-Phase II and Control Phase III samples, while there was no difference in fecal bacterial composition between Control-Phase III and Peptiva Phase III samples. For instance, OTUs Ssd-00019,and Ssd-00053 were in lower abundance in Peptiva Phase II samples compared to Control Phase II (P < 0.05), but no significant difference was observed in the abundance of these two OTUs when comparing Peptiva Phase II to Control Phase III (P > 0.05).

Conclusions

Together, these results suggest that Peptiva can modulate the composition of the swine microbiome during a specific window of the nursery stage, potentially by accelerating its maturation.

Background

Weaning is one of the most critical transition stages of the swine production cycle, as decreased feed intake and poor performance from sudden changes in diet and environment can result in severe economic losses [1, 2]. While a number of physiological conditions contribute to the performance and health challenges that commonly occur during the nursery phase, gastrointestinal dysfunction is generally involved. Typically, a combination of prolonged intestinal inflammation, immature immune system and transitioning gut microbial communities result in a compromised gut epithelial lining, decreased nutrient digestibility, and increased susceptibility to pathogen infection [1, 3,4,5,6,7,8,9,10,11]. Together, these conditions can lead to a higher incidence of diarrhea, resulting in higher weaned pig morbidity and mortality.

Conventional approaches to reduce the impact of weaning on nursery pig health and performance have typically combined antibiotic use to reduce the pathogen load with inclusion of high-quality protein ingredients to facilitate digestion and absorption [12]. However, implementation of stricter regulations on the prophylactic use of medically important antimicrobials, as well as higher costs of traditionally used protein sources such as fish meal, have created a need for effective substitutes and the development of innovative strategies. For instance, products such as essential oils and antimicrobial peptides are becoming more widely used as alternative antimicrobials, while modified plant ingredients with reduced levels of anti-nutritional factors (e.g. enzymatically or microbially modified soybean meal) are being included as lower cost protein-rich sources in dietary formulations [13,14,15]. In addition to these substitutes, feed additives are also developed to target other functions, such as enhancing the immune response of weaned pigs (e.g. immunoglobulin or omega-3 fatty acids), stimulating digestive functions (e.g. butyrate, glutamate, threonine or cysteine), or promoting the establishment of beneficial gut microorganisms (probiotics, prebiotics) [16,17,18].

Amongst the various products available, peptides have the unique potential to be used as multipurpose feed additives. Indeed, they are cost effective means of providing amino acids, as they are more stable, soluble, and can be absorbed at a faster rate than free amino acids [19,20,21]. In addition, certain types of peptides can control various physiological functions by acting as either antimicrobials, antioxidants, immuno-modulators or signaling molecules [20, 22, 23]. In the case of bioactive peptides supplemented in feed, they may act on either host cells and / or on the host’s microbiome [24,25,26]. As an example of peptide signaling to host cells, exorphins have been shown to modulate gastrointestinal motility, secretions, and endocrine metabolism once they have been released by digestion and absorbed by the gut epithelium [27]. Conversely, modulation of gut microbiome composition by certain antimicrobial peptides has also been reported. For instance, colicins and cecropin AD can help control the proliferation of Escherichia coli strains that can cause post weaning diarrhea in swine [25, 28, 29]. Antimicrobial peptides can also have positive effects on performance. Indeed, feeding a combination of lactoferrin, cecropin, defensin and plectasin resulted in higher average daily gain and final body weight compared to unsupplemented diets [30]. Similarly, apparent total tract digestibility of either dry matter or crude protein was found to be higher with dietary supplementation of the antimicrobial peptide-P5 [31].

Considering the importance of beneficial gut microbial communities for animal health and nutrition, manipulating the gut microbiome using peptides would represent an additional tool towards improving resistance to pathogens, optimizing the use of alternative feed ingredients or providing other benefits to the host animal. Typically, bioactive peptides remain inactive until they are released from their parent protein as a result of chemical, enzymatic, or microbial hydrolysis [32]. Since their functional characteristics would depend on their length as well as their amino acid composition and sequence [20, 23], there likely exists a wide range of potential bioactive peptides that have yet to be identified or characterized. Indeed, the search for novel bioactive peptides is still ongoing even for highly investigated sources such as milk [33]. Thus, a reasonable expectation would then be that many peptide formulations would contain bioactive peptides that can perform functions other than simply supplying dietary amino acids. However, as the effects of peptide feed additives on the gut microbiome of food animals remain largely unexplored, additional insight is required to develop further improvements in this field.

In this context, the aim of the study presented in this report was to determine the effect of a commercially-formulated peptide additive, Peptiva, on the performance and fecal bacterial community profiles of weaned pigs raised in a commercial wean-to-finish swine facility. This product has been previously reported as an acceptable protein supplement in nursery diets [34], but had not been tested at low inclusion levels. In the current study, Peptiva supplementation did not result in improved weight gains or feed efficiency of weaned pigs under the conditions tested, but it was found to affect the fecal microbiome composition of animals during the first few weeks after weaning.

Results

Effect of low inclusion of Peptiva on swine performance during the nursery phase

To test the ability of Peptiva to improve the availability of dietary amino acids in swine nursery phase diets, animals fed a Peptiva-supplemented diet that included only 90% of the recommended amino acids requirements for nursery phase diets (PEP10) were compared to animals fed the control diet (CON). After the first 3 weeks, no effect of diet on body weight was observed (P > 0.05; Table 1). Starting at week 4, however, pigs fed the CON diet tended to be heavier than PEP10-fed pigs (P = 0.07), with CON-fed pigs continuing to be heavier than PEP10-fed pigs through to week 6 (P < 0.05). While there was no difference in average daily feed intake during Phases I and II across dietary treatments, an effect of diet on daily feed intake during Phase III was observed, where CON-fed pigs had greater daily intake than PEP-fed pigs (P < 0.05). No effect of dietary treatment on average daily weight gain or gain:feed was observed. While there were 3% fewer pigs removed from the PEP10 group compared to the CON group for the entire trial period (6 wks), a statistical difference in net pig removal rate by diet was not detected. No significant differences were noted in pen weight variation amongst treatment groups.

Table 1 Growth performance of weaned pigs fed diets containing Peptiva formulated at 100 (PEP) or 90% (PEP-10) of amino acid requirements (NRC (2012)1

Effect of diet composition and Peptiva supplementation on the fecal bacterial profile of weaned pigs

To investigate the potential of Peptiva as a modulator of gut microbiome composition in weaned pigs, a comparative analysis using fecal bacterial communities as a proxy was performed on samples collected at the end of Phase II and at the end of Phase III. The average number of high-quality, non-chimeric reads for 16S rRNA gene sequences across the four sample sets (CON II, CON III, PEP II and PEP III) ranged from 14,972 ± 2792 to 26,020 ± 3191 (Table 2), with numerical differences amongst means not found to be significant (P = 0.16). Firmicutes was the most highly represented phylum, with sample set averages ranging from 77.4 to 85.3% (Table 3; Supplementary File 3). While these variations in abundance at the phylum level were found to be only numerical, the differences in representation for two families belonging to Firmicutes were found to be significant (P < 0.05). Lactobacillaceae were more abundant in CON II samples (44.8%) than in samples from pigs fed the PEP II, CON III or PEP III diets (13.0–16.0%). In contrast, Streptococcaceae were in lower abundance in CON II compared to CON III and PEP III (P < 0.05). Other well represented families belonging to Firmicutes included Lachnospiraceae (5.9–13.2%) and Clostridiaceae1 (5.9–18.9%), but the observed differences in abundance were not significant (P > 0.05). The second most abundant phylum was Bacteriodetes, with Prevotellaceae identified as its most highly represented family (11.8% - 16.0); variation across datasets was not found to be significant for either of these taxonomic groups.

Table 2 High quality sequence read yield from each sample set
Table 3 Relative abundance (%) of the main bacterial taxonomic groups in representative fecal samples from Control and Peptiva-fed pigs in each of Phase II and III diets

Comparative analysis of fecal bacterial composition by alpha and beta diversity

Community level compositional differences amongst fecal bacterial communities from CON II, CON III, PEP II and PEP III sample sets were further assessed using alpha and beta diversity analyses. A combined total of 8429 OTUs were identified across all samples analyzed (Supplementary File 2). No statistical difference was observed amongst means of the four dietary treatments for either observed OTUs, Ace, Chao1, Shannon or Simpson indices (P > 0.05; Table 4). However, principal coordinate analysis (PCoA) based on Bray-Curtis OTU composition dissimilarity revealed that samples could be clustered into three different groups according to their fecal bacterial community composition (Fig. 1). Furthermore, uneven distribution of samples from different sets amongst the three clusters of the PCoA plot suggested that distinct OTU profiles could be associated with the fecal environments of particular sets of samples.

Table 4 Alpha diversity indices and coverage from Control and Peptiva-fed pigs in each of Phase II and III diets
Fig. 1
figure 1

Comparison of fecal bacterial communities from weaned pigs under two different diets at two different time points. Principle Coordinate Analysis (PCoA) was performed using OTU composition-based Bray-Curtis dissimilarity. The x and y axes correspond to Principal Components 1 (PC1) and 2 (PC2), which explained the highest level of variation

Identification of weaned pig OTUs responding to distinct dietary treatments

As the comparative taxonomic composition analysis and PCoA both indicated differences in bacterial composition amongst sample sets, the individual profiles of major OTUs were further investigated. A total of 23 OTUs that were found to have a mean relative abundance of at least 1% in at least one sample set were designated as major OTUs. Of these most abundant OTUs, at least seven were likely to correspond to uncharacterized species, as they each showed less than 95% sequence identity to their respective closest valid taxon. Thirteen major OTUs, all affiliated to Firmicutes, were found to vary across sample types (P < 0.05) (Table 5; Supplementary File 3). Pair-wise differences between specific samples for nine of these varying OTUs were further revealed by the post hoc Nemenyi test (Fig. 2). Notably, the respective abundances of OTUs Ssd-00019 and Ssd-00053 were found to be significantly different in CON II compared to PEPII, CON III and PEP III sample sets (P < 0.05). OTUs Ssd-00002, Ssd-00025, Ssd-00039, Ssd-00048, Ssd-00061 and Ssd-00106 showed a slightly different profile, with their respective abundances being significantly different between CON II and either CON III or PEP III (P < 0.05), while no significant difference was found between PEPII and either CON II, CONIII or PEP III. Also, while Ssd-00140 was found at similar levels in CON II and PEP II, its abundance in these sample sets was significantly lower than in CON III and PEP III (P < 0.05).

Table 5 Relative abundance (%) of the most abundant OTUs in representative fecal samples from Control and Peptiva-fed pigs in each of Phase II and III diets
Fig. 2
figure 2

Main bacterial OTUs whose respective abundance was found to vary significantly amongst groups based on the post hoc Nemenyi test for multiple pairwise comparisons (P < 0.05). OTUs affiliated to the genus Lactobacillus are shown in panel (A) while OTUs affiliated to the genera Streptococcus or Roseburia are shown in panel (B). For each OTU, means with different superscripts were significantly different as determined by the post hoc Nemenyi test for multiple pairwise comparisons

Associations between main OTUs and dietary treatments

A correspondence analysis was conducted to further explore potential associations between main OTUs and dietary treatments (Fig. 3). All CON II samples clustered together with OTUs Ssd-00002, Ssd-00019, Ssd-00025, Ssd-00053 and Ssd-000106. CON III and PEP III samples were clustered into two groups, with the major group being closely associated with OTUs SSd-00048, OTUs SSd-00061 and OTUs SSd-00140, while the minor group was closely associated with OTU Ssd-00001. PEP II samples showed a very distinct distribution pattern, as half of the samples clustered with the CON II group, while the remaining samples were associated with the CON III - PEP III major cluster.

Fig. 3
figure 3

Correspondence analysis (CA) between sample type (circle) and main OTUs (star)

Discussion

Products manufactured by hydrolysis of conventional protein ingredients have the potential to include bioactive peptides that can provide other functions or benefits in addition to supplying dietary amino acids. In this study, a commercial peptide-based additive, Peptiva, was tested as a possible source of bioactive molecules using two methods. First, its ability to compensate for reduced inclusion of dietary amino acid levels in weaned pig diets, by increasing the digestibility or the efficiency of use of protein ingredients, was assessed. In the context of a commercial swine production system as used in this study, there was no difference in performance during Phases I and II post-weaning, but PEP 10-fed pigs were found to weigh significantly less than CON-fed pigs by the end of Phase III. These results would indicate that, at least in the first 6 weeks post-weaning, Peptiva supplementation at low inclusion levels was not sufficient to compensate for a 10% reduction in dietary amino acid levels. Intriguingly, the average daily feed intake for pigs fed the PEP diet was lower than for pigs fed the CON diet during Phase III, but no differences in average daily gain were observed. While future research will be required to elucidate the possible mechanisms involved, perhaps the Peptiva formulation includes bioactive peptides that can improve the efficiency of nutrient utilization by the host, acting on the absorptive capacity of enterocytes or on their metabolism. Alternatively, Peptiva could act indirectly by modulating the activity of microbial symbionts. In addition, it is worth noting that 3% fewer pigs were removed from the trial, either because of mortality or illness, in the PEP-10 group compared to the other treatments, suggesting a benefit to overall pig health. Even if this observation is not supported by statistical testing, a lower removal rate in the absence of herd-wide antibiotic administration, as was implemented in this study, represents an important finding for the swine industry, since mass administration of antibiotics in the first 7 to 10 days following weaning is not allowed in conventional production.

The second potential activity of the Peptiva product investigated in this study was the ability to change or modulate the composition of the gut microbiome in weaned pigs. Since the composition of gut microbial communities has been associated with the health status and performance of individual hosts [3, 35,36,37], compounds that can change gut symbiont profiles have the potential to be developed as tools to improve critical livestock production parameters [38]. To this end, fecal bacterial communities were used as a proxy for gut microbiome composition analysis in weaned pigs, from which two main observations were made: evidence of bacterial succession between Phase II and Phase III in control-fed animals, and a stage-specific effect in Pep-fed pigs.

Bacterial succession from Lactobacillaceae in phase II to Streptococcaceae in phase III

A comparison of the samples collected from control-fed pigs between Phase II and Phase III diets was suggestive of microbial succession, as major changes in taxonomic profiles and OTU composition were observed. For instance, members of the Lactobacillaceae family were found to be more abundant at the end of Phase II compared to the end of Phase III, which included four main OTUs (Ssd-0002, Ssd-00019, Ssd-00025, and Ssd-00053). In young animals, Lactobacilli have been reported to prevent adhesion of pathogens to the gut mucosa, inhibit growth of pathogens through production of lactate, and / or stimulate colonization of beneficial bacteria [9, 39,40,41]. Because of these types of activities, Lactobacillus species are considered beneficial to the gastrointestinal tract of animals, and are typically included in probiotic formulations. For instance, a probiotic formulation containing L. gasseri, L. reuteri, L. acidophilus and L. fermentum was reported to result in fewer incidences of diarrhea in weaned pigs and to lower E. coli counts after a pathogen challenge [39], while weaned pigs supplemented with L. reuteri were found to have higher average daily gain, longer ileal villi, as well as increased expression of the tight junction protein zonula occludens − 1 [42]. Lactobacilli have also been reported to have antimicrobial activity, as observed with L. reuteri which can inhibit the growth and mucosal adherence of enterotoxigenic E. coli [43], and L. gasseri which is known to produce a bacteriocin [44]. In the current study, three of the four most abundant Lactobacillus-affiliated OTUs were found to be closely related to L. reuteri or L. gasseri.

At the end of Phase III, members of the Streptococcaceae family became the most predominant bacterial group of the fecal microbiome in weaned pigs, while the abundance of Lactobacillus-affiliated bacteria was greatly reduced. Since the sequence identity to their respective closest Streptococcus relatives ranged between 90 and 96%, main OTUs Ssd-00039, Ssd-00048, Ssd-00061 and Ssd-00140 most likely corresponded to uncharacterized species of this genus. While the biological activities of Streptococci in the gut have not been as extensively studied as for Lactobacilli, members of this genus are also known to be lactate producers and to express bacteriocin, and thus could be involved in protection against pathogen proliferation in weaned pigs [45].

Of the factors that may be responsible for these observed changes in bacterial composition in pigs fed control diets, differences in diet formulation between Phase II and Phase III offer a reasonable explanation. Notably, three ingredients (dried whey, fish meal and zinc oxide) were included in Phase II diets, but not in Phase III diets (Supplementary File 1). As its primary use is to prevent diarrhea, zinc oxide represents a likely candidate modulator of gut microbiome composition [46,47,48,49,50]. However, its target bacterial groups in gut environments remain to be further investigated, as exemplified by two conflicting studies, one observing a decrease in Lactobacilli as a result of dietary inclusion of zinc oxide [48], while the other reported no effect [51]. Similarly, further investigations will be required to determine the effects of dried whey and fish meal, both used as high-quality protein ingredients, on the gut microbiome of weaned pigs.

Stage-specific effect of Peptiva on the microbiome of weaned pigs

The second main observation from the comparative analysis of fecal bacterial communities performed in this study was that the profiles of PEP II samples appeared to be intermediate between CON II and CON III profiles. This was well illustrated by correspondence analysis, where PEP II samples appeared to be divided into two groups, with certain samples more similar to CON II profiles while others were more similar to CON III profiles. However, the analysis of additional samples would have been beneficial in providing increased resolution and confidence in support of this observation. At the OTU level, the respective abundances of Ssd-00019 and Ssd-00053 in PEP II were found to be statistically different from CON II (P < 0.05), but not from CON III. In contrast, no difference in abundance was found for Ssd-000140 between PEP II and CON II samples, which were however both significantly lower than those observed in the CON III samples. Other OTUs, such as Ssd-00002, Ssd-00025, Ssd-00039, Ssd-00048, Ssd-00061, and Ssd-00106, were found to be statistically different between CON II and CON III, while no significant pair-wise difference was found between either CON II and PEP II or between CON III and PEP II. Finally, no major differences in fecal bacterial profiles were observed between CON III and PEP III samples, indicating that both sets of fecal bacterial communities had reached similar compositional profiles. While additional research will be required to further elucidate the mechanisms responsible for these effects, the results presented in this study would suggest that Peptiva can promote maturation of swine fecal bacterial communities during a specific period of the nursery phase.

Conclusions

Under the conditions tested in this study, Peptiva supplementation did not result in improved weight gains or feed efficiency of weaned pigs, but it was found to reduce average daily feed intake during Phase III of the nursery trial. In addition, the results presented in this report suggest that Peptiva can affect the fecal microbiome composition of animals during the first few weeks after weaning. In the context of the current understanding of gut microbiome development, early events that impact bacterial composition can have long term effects that persist in adults. For food animal production, this would suggest that development of practices or diet formulations that can establish more resistant, resilient and efficient gut microbiomes in neonates would provide lasting benefits into the growing and finishing stages. Based on their potential for producing a vast array of bioactive molecules, peptide formulations represent a largely untapped source of compounds that could be further developed into feed additives to benefit animal health and nutrition.

Methods

Animal performance trial and sample collection

The animal trial was conducted at the South Dakota State University (SDSU) Off-Site Wean-to-Finish Barn, with all procedures approved by the SDSU Institutional Animal Care and Use Committee before the start of the study (Protocol 17-035A). This swine facility is managed as a commercial-scale livestock barn to conduct nutritional and animal health research that can benefit producers in this sector. Weaned pigs (21 d of age, 5.6 ± 1.2 kg) were randomly allocated to 48 pens (24 pigs/pen), with each pen randomly assigned to one of three experimental diets (n = 16 pens/treatment): control diet (CON; formulated to meet the NRC (2012) nutrient requirements), Peptiva (PEP; control diet supplemented with Peptiva), and PEP with reduced amino acid content (PEP10; dietary amino acid content at 90% of NRC (2012) recommendations). All other dietary nutrients met or exceeded NRC (2012) recommendations for weaned pigs. Experimental diets were fed according to a standard nursery phase feeding program (Supplementary File 1): Phase I (d0-d7), Phase II (d8-d21), and Phase III (d22–42). Peptiva is a commercial product manufactured by Vitech Bio-Chem Corporation (Glendale, CA, USA) which consists of fish peptides, porcine digests and microbial peptides. In both PEP and PEP10 diets, Peptiva was included at 1, 0.5, and 0.3%, during Phases I, II, and III, respectively. The swine facility was divided into four blocks based on pen location within the barn, and each treatment was equally represented in each block (n = 4 pens/treatment/block). Pens of pigs were assigned to treatment within each block with consideration of pen weight to be equivalent between treatments as best as possible.

During the trial, one pen for the PEP group and one pen for the PEP10 group were fed the wrong diet. Consequently, these two pens were removed from the analysis, which was conducted using a total of 46 pens instead of 48 pens as originally planned (CON: n = 16; PEP: n = 15; PEP10: n = 15).

No mass antibiotic treatment via feed or water medicator was used during the course of the trial. Injectable antibiotics were only administered on an individual pig basis for treatment of scours or poor health. Individual pigs treated with injectable antibiotics that recovered from their symptoms remained in the performance trial, but they were not used to collect samples for composition analysis of fecal bacterial communities.

Body weights of the animals were measured by pen at the start of the trial, then on a weekly basis until the end of Phase III. Individual pig weights were determined at the beginning of the trial, at the end of Phase II and at the end of Phase III. Samples for microbiome analysis were collected at the end of Phase II and at the end of Phase III from ten animals fed the CON diets and ten individuals fed the PEP diet. More specifically, two representative individuals from each of five representative pens were selected for fecal sample collection for each diet. Pen weight was used to identify representative pens for each dietary treatment, and individual weight was used to identify representative animals from each selected pen. Fecal samples were collected by rectal palpation, then stored frozen (− 20 °C) until microbial genomic DNA extraction was performed.

At the conclusion of the trial, pens were randomly allotted to a separate grow finish trial, and the animals were marketed after achieving 130 kg body weight.

Microbial DNA isolation and PCR amplification of the 16S rRNA gene

Microbial genomic DNA was isolated from fecal samples using the repeated bead beating plus column method, as previously described [52]. The V1-V3 region of the bacterial 16S rRNA gene was PCR-amplified using the 27F forward [53] and 519R reverse [54] primer pair. PCR reactions were performed with the Phusion Taq DNA polymerase (Thermo Scientific, Waltham, MA, USA) under the following conditions: hot start (4 min, 98 °C), followed by 35 cycles of denaturation (10 s, 98 °C), annealing (30 s, 50 °C) and extension (30 s, 72 °C), then ending with a final extension period (10 min, 72 °C). PCR products were separated by agarose gel electrophoresis, and amplicons of the expected size (~ 500 bp) were excised for gel purification using the QiaexII Gel extraction kit (Qiagen, Hilden, Germany). For each sample, approximately 400 ng of amplified DNA were submitted to Molecular Research DNA (MRDNA, Shallowater, TX, USA) for sequencing with the Illumina MiSeq (2X300) platform to generate overlapping paired end reads.

Computational analysis of PCR generated 16S rRNA amplicon sequences

Unless specified, sequence data analysis was performed using custom written Perl scripts (Supplementary File 4). Raw bacterial 16S rRNA gene V1-V3 amplicon sequences were provided by Molecular Research DNA as assembled contigs from overlapping MiSeq (2X300) paired-end reads from the same flow cell clusters. Reads were then selected to meet the following criteria: presence of both intact 27F (forward) and 519R (reverse) primer nucleotide sequences, length between 400 and 580 nt, and a minimal quality threshold of no more than 1% of nucleotides with a Phred quality score lower than 15.

Following quality screens, sequence reads were aligned, then clustered into Operational Taxonomic Units (OTUs) at a genetic distance cutoff of 5% sequence dissimilarity [55]. While 3% is the most commonly used clustering cutoff for 16S rRNA, it was originally recommended for full length sequences, and may not be suitable for the analysis of specific subregions since nucleotide sequence variability is not constant across the entire length of the 16S rRNA gene. In this context, if 3% is a commonly accepted clustering cutoff for V4 or V4–V5 regions, which are the least variable of the hypervariable regions, then a higher cutoff should be used for the V1-V3 region, since V1 is the most variable region of the 16S rRNA gene. OTUs were screened for DNA sequence artifacts using the following methods. Chimeric sequences were first identified with the chimera.uchime and chimera.slayer commands from the MOTHUR open source software package [56]. Secondly, the integrity of the 5′ and 3′ ends of OTUs was evaluated using a database alignment search-based approach; when compared to their closest match of equal or longer sequence length from the NCBI nt database, as determined by BLAST [57], OTUs with more than five nucleotides missing from the 5′ or 3′ end of their respective alignments were discarded as artifacts. Single read OTUs were subjected to an additional screen, where only sequences that had a perfect or near perfect match to a sequence in the NCBI nt database were kept for analysis, i.e. that the alignment had to span the entire sequence of the OTU, and a maximum of 1% of dissimilar nucleotides was tolerated.

While a 5% cutoff may affect comparisons with other studies, its impact should be minimal. It would not be expected to affect the most abundant OTUs, which were the main focus of this study. While it would affect the total number of OTUs generated, the impact would be less than that of other practices that are commonly used in microbiome data analysis. For instance, the removal of all OTUs with a low number of reads would have had a greater impact, as these OTUs are by far the most abundant.

After removal of sequence chimeras and artifacts, taxonomic assignment of valid OTUs was determined using a combination of RDP Classifier [58] and BLAST [57]. The List of Prokaryotic Names with Standing in Nomenclature (LPSN - http://www.bacterio.net) was also consulted for information on valid species belonging to taxa of interest [59].

Computational analysis for alpha and beta diversity

Using custom Perl scripts, all datasets were randomly rarefied to 3000 reads, which were then used to create ‘shared’-type formatted files. Ace, Chao1, Shannon and Simpson indices, as well as observed OTUs and coverage, were determined from the shared files using summary.single in MOTHUR [56]. For Principal Coordinate Analysis (PCoA), Bray-Curtis distances were first determined using summary.shared, which were then used as input for the command pcoa, with both procedure performed using MOTHUR [56]. Principal Components 1 (PC1) and 2 (PC2), representing the highest levels of variation, were plotted using Microsoft Excel. Correspondence Analysis (CA) was conducted in R (version 3.6.1) using the command ca from the R package ‘ca’. Outputs were plotted using Microsoft Excel.

Statistical analyses

Analysis of performance data was performed using the PROC MIXED procedure of SAS (Version 9.4; SAS Inst. Inc., Cary, NC) with pen as the experimental unit and pen nested within block as the random variable. Dietary treatment was considered the fixed effect, and the effect of block was removed from the model because it was not significant. Data were a priori tested for normal distribution and homogeneity of variances. Initial body weight was used as covariate for analysis of weekly body weight. All possible comparisons between treatment means were tested using the PDiff option in SAS in combination with Tukey’s adjusted means separation when a significant main effect was observed, and data are presented as lsmeans +/− standard error of the mean. A Chi-squared test was used to evaluate the distribution of total pigs removed by treatment.

For comparison of bacterial taxonomic groups, alpha diversity indices, and OTU abundance amongst different sample groups (CON-Phase II, PEP-Phase II, CON-Phase III and PEP-Phase III), the non-parametric Friedman test (command friedman.test) and the post hoc Nemenyi test for multiple pairwise comparisons (command posthoc.friedman.nemenyi.test) were performed in R (Version R-3.2.3). Groups were considered to be significantly different when P ≤ 0.05.

Accession numbers for next generation sequencing data

Raw sequence data are available from the NCBI Sequence Read Archive under Bioproject PRJNA533644 and SRA accession SRP192997.

Availability of data and materials

Raw sequence data are available from the NCBI Sequence Read Archive under Bioproject PRJNA533644 and SRA accession SRP192997.

Abbreviations

16S rRNA:

16S ribosomal RNA

BLAST:

Basic Local Alignment Search Tool

CON II:

Samples collected at the end of Phase II from animals fed the control diet

CON III:

Samples collected at the end of Phase III from animals fed the control diet

CON:

Control diet

LPSN:

List of Prokaryotic Names with Standing in Nomenclature

NCBI:

National Center for Biotechnology Information

OTU:

Operational Taxonomic Unit

PCoA:

Principal Coordinate Analysis

PEP II:

Samples collected at the end of Phase II from animals fed the Peptiva diet

PEP III:

Samples collected at the end of Phase III from animals fed the Peptiva diet

PEP:

Control diet supplemented with Peptiva

PEP10:

Control diet containing 90% of the recommended amino acid levels, supplemented with Peptiva

RDP:

Ribosomal Database Project

SDSU:

South Dakota State University

References

  1. Campbell JM, Crenshaw JD, Polo J. The biological stress of early weaned piglets. J Anim Sci Biotechnol. 2013;4:19.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Hötzel MJ, de Souza GPP, Costa OAD, Machado Filho LCP. Disentangling the effects of weaning stressors on piglets’ behaviour and feed intake: changing the housing and social environment. Appl Anim Behav Sci. 2011;135(1):44–50.

    Article  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  4. Heo JM, Opapeju FO, Pluske JR, Kim JC, Hampson DJ, Nyachoti CM. Gastrointestinal health and function in weaned pigs: a review of feeding strategies to control post-weaning diarrhoea without using in-feed antimicrobial compounds. J Anim Physiol Anim Nutr. 2013;97(2):207–37.

    Article  CAS  Google Scholar 

  5. Pluske JR, Hampson DJ, Williams IH. Factors influencing the structure and function of the small intestine in the weaned pig: a review. Livest Sci. 1997;51(1):215–36.

    Article  Google Scholar 

  6. Moeser AJ, Pohl CS, Rajput M. Weaning stress and gastrointestinal barrier development: implications for lifelong gut health in pigs. Anim Nutr. 2017;3(4):313–21.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Pluske JR, Pethick DW, Hopwood DE, Hampson DJ. Nutritional influences on some major enteric bacterial diseases of pigs. Nut Res Rev. 2002;15:333–71.

    Article  CAS  Google Scholar 

  8. Kim JC, Hansen CF, Mullan BP, Pluske JR. Nutrition and pathology of weaner pigs: nutritional strategies to support barrier function in the gastrointestinal tract. Anim Feed Sci Tech. 2012;173(1–2):3–16.

    Article  CAS  Google Scholar 

  9. Fouhse JM, Zijlstra RT, Willing BP. The role of gut microbiota in the health and disease of pigs. Animal Frontiers. 2016;6(3):30–6.

    Article  Google Scholar 

  10. Lee YK, Mazmanian SK. Has the microbiota played a critical role in the evolution of the adaptive immune system? Science. 2010;330(6012):1768–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Brestoff JR, Artis D. Commensal bacteria at the interface of host metabolism and the immune system. Nature Immunol. 2013;14:676.

    Article  CAS  Google Scholar 

  12. Maxwell JCV, Carter SD. Feeding the weaned pig. Pages 691–716 In Swine Nutrition, Lewis, A.J. and Southern L.L. editors, CRC Press LLC., Boca Raton, FL. 2001.

  13. Franz C, Baser KHC, Windisch W. Essential oils and aromatic plants in animal feeding – a European perspective. Rev Flavour Fragr J. 2010;25:327–40.

    Article  CAS  Google Scholar 

  14. Koepke JR, Kaushik RS, Gibbons WR, Brown M, Levesque CL. Evaluation of a bio-processed soybean meal on nursery pig performance and immune status. J Anim Sci. 2017;95:5030–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Sinn S, Gibbons WR, Brown ML, DeRouchey JM, Levesque CL. Evaluation of microbially-enhanced soybean meal as an alternative to fishmeal in weaned pig diets. Animal. 2017;11:784–93.

    Article  CAS  PubMed  Google Scholar 

  16. de Lange CFM, Pluske J, Gong J, Nyachoti CM. Strategic use of feed ingredients and feed additives to stimulate gut health and development in young pigs. Livest Sci. 2010;134(1–3):124–34.

    Article  Google Scholar 

  17. Berrocoso JD, Serrano MP, Cámara L, Rebollar PG, Mateos GG. Influence of diet complexity on productive performance and nutrient digestibility of weaning pigs. Anim Feed Sci Tech. 2012;171(2):214–22.

    Article  Google Scholar 

  18. Yuan L, Chang J, Yin Q, Lu M, Di Y, Wang P, Wang Z, Wang E, Lu F. Fermented soybean meal improves the growth performance, nutrient digestibility, and microbial flora in piglets. Anim Nutr. 2017;3(1):19–24.

    Article  PubMed  Google Scholar 

  19. Lindemann MD, Cromwell GL, Monegue HJ, Cook H, Soltwedel KT, Thomas S, Easter RA. Feeding value of an enzymatically digested protein for early-weaned pigs. J Anim Sci. 2000;78:318–27.

    Article  CAS  PubMed  Google Scholar 

  20. Hou Y, Wu Z, Dai Z, Wang G, Wu G. Protein hydrolysates in animal nutrition: industrial production, bioactive peptides, and functional significance. J Anim Sci Biotechnol. 2017;8:24.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Webb KE Jr, Matthew JC, DiRienzo DB. Peptide absorption: a review of current concepts and future perspectives. J Anim Sci. 1992;70(10):3248–57.

    Article  CAS  PubMed  Google Scholar 

  22. Bhat ZF, Kumar S, Bhat HF. Bioactive peptides of animal origin: a review. J Food Sci Technol. 2015;52(9):5377–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Nasri M: Chapter Four - Protein Hydrolysates and Biopeptides: Production, Biological Activities, and Applications in Foods and Health Benefits. A Review. In: Adv Food Nutr Res. Edited by Toldrá F, vol. 81: Academic Press; 2017: 109–159.

  24. Xiong X, Yang HS, Li L, Wang YF, Huang RL, Li FN, Wang SP, Qiu W. Effects of antimicrobial peptides in nursery diets on growth performance of pigs reared on five different farms. Livest Sci. 2014;167:206–10.

    Article  Google Scholar 

  25. Wang S, Zeng X, Yang Q, Qiao S. Antimicrobial peptides as potential alternatives to antibiotics in food animal industry. Int J Mol Sci. 2016;17(5).

    Article  PubMed Central  CAS  Google Scholar 

  26. Xiao H, Shao F, Wu M, Ren W, Xiong X, Tan B, Yin Y. The application of antimicrobial peptides as growth and health promoters for swine. J Anim Sci Biotechnol. 2015;6(1):19.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Froetschel MA. Bioactive peptides in digesta that regulate gastrointestinal function and intake. J Anim Sci. 1996;74(10):2500–8.

    Article  CAS  PubMed  Google Scholar 

  28. Stahl CH, Callaway TR, Lincoln LM, Lonergan SM, Genovese KJ. Inhibitory activities of colicins against Escherichia coli strains responsible for postweaning diarrhea and edema disease in swine. Antimicrob Agents Chemother. 2004;48(8):3119–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Wu S, Zhang F, Huang Z, Liu H, Xie C, Zhang J, Thacker PA, Qiao S. Effects of the antimicrobial peptide cecropin AD on performance and intestinal health in weaned piglets challenged with Escherichia coli. Peptides. 2012;35(2):225–30.

    Article  CAS  PubMed  Google Scholar 

  30. Tang XA, Fatufe A, Yin Y, Tang Z, Wang S, Liu Z, Xinwu LT-J. Dietary supplementation with recombinant Lactoferrampin-Lactoferricin improves growth performance and affects serum parameters in piglets. J Anim Vet Adv. 2012;11:2548–55.

    Article  CAS  Google Scholar 

  31. Yoon JH, Ingale SL, Kim JS, Kim KH, Lee SH, Park YK, Kwon IK, Chae BJ. Effects of dietary supplementation of antimicrobial peptide-A3 on growth performance, nutrient digestibility, intestinal and fecal microflora and intestinal morphology in weanling pigs. Anim Feed Sci Tech. 2012;177(1–2):98–107.

    Article  CAS  Google Scholar 

  32. Korhonen H. Milk-derived bioactive peptides: from science to applications. J Funct Foods. 2009;1(2):177–87.

    Article  CAS  Google Scholar 

  33. Zanutto-Elgui MR, Vieira JCS, DZd P, MAR B, PdM P, Elgui de Oliveira D, Fleuri LF. Production of milk peptides with antimicrobial and antioxidant properties through fungal proteases. Food Chem. 2019;278:823–31.

    Article  CAS  PubMed  Google Scholar 

  34. Zhao J, Harper AF, Perkins BK, Southern LL, Shelton JL, Bidner TD, Webb KE Jr, Estienne MJ, Kuehn LA. Assessment of a marine-based hydrolyzed protein source and spray-dried plasma protein as supplements in the diet of early weaned pigs. Prof Anim Sci. 2008;24(6):604–13.

    Article  Google Scholar 

  35. Bäckhed F, Roswall J, Peng Y, Feng Q, Jia H, Kovatcheva-Datchary P, Li Y, Xia Y, Xie H, Zhong H, et al. Dynamics and stabilization of the human gut microbiome during the first year of life. Cell Host Microbe. 2015;17(5):690–703.

    Article  PubMed  CAS  Google Scholar 

  36. Kim HB, Isaacson RE. The pig gut microbial diversity: understanding the pig gut microbial ecology through the next generation high throughput sequencing. Vet Microbiol. 2015;177(3):242–51.

    Article  CAS  PubMed  Google Scholar 

  37. Richards JD, Gong J, de Lange CFM. The gastrointestinal microbiota and its role in monogastric nutrition and health with an emphasis on pigs: current understanding, possible modulations, and new technologies for ecological studies. Can J Anim Sci. 2005;85(4):421–35.

    Article  Google Scholar 

  38. Han GG, Lee J-Y, Jin G-D, Park J, Choi YH, Kang S-K, Chae BJ, Kim EB, Choi Y-J. Tracing of the fecal microbiota of commercial pigs at five growth stages from birth to shipment. Sci Rep. 2018;8(1):6012.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Huang C, Qiao S, Li D, Piao X, Ren J. Effects of Lactobacilli on the performance, diarrhea incidence, VFA concentration and gastrointestinal microbial flora of weaning pigs. Asian-Australas J Anim Sci. 2004;17(3):401–9.

    Article  Google Scholar 

  40. Valeriano VD, Balolong MP, Kang DK. Probiotic roles of Lactobacillus sp. in swine: insights from gut microbiota. J Appl Microbiol. 2017;122(3):554–67.

    Article  CAS  PubMed  Google Scholar 

  41. Yang F, Hou C, Zeng X, Qiao S. The use of lactic acid bacteria as a probiotic in swine diets. Pathogens. 2015;4(1):34–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Yi H, Wang L, Xiong Y, Wen X, Wang Z, Yang X, Gao K, Jiang Z. Effects of Lactobacillus reuteri LR1 on the growth performance, intestinal morphology, and intestinal barrier function in weaned pigs. J Anim Sci. 2018;96(6):2342–51.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Wang T, Teng K, Liu G, Liu Y, Zhang J, Zhang X, Zhang M, Tao Y, Zhong J. Lactobacillus reuteri HCM2 protects mice against enterotoxigenic Escherichia coli through modulation of gut microbiota. Sci Rep. 2018;8(1):17485.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Ritter P, Kohler C, von Ah U. Evaluation of the passage of Lactobacillus gasseri K7 and bifidobacteria from the stomach to intestines using a single reactor model. BMC Microbiol. 2009;9(1):87.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. Georgalaki MD, Van den Berghe E, Kritikos D, Devreese B, Van Beeumen J, Kalantzopoulos G, De Vuyst L, Tsakalidou E. Macedocin, a food-grade lantibiotic produced by Streptococcus macedonicus ACA-DC 198. Appl Environ Microbiol. 2002;68(12):5891–903.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Starke IC, Pieper R, Vahjen W, Zentek J. The impact of dietary zinc oxide on the bacterial diversity of the small intestinal microbiota of weaned piglets. J Vet Sci Technol. 2014;5:171.

    Google Scholar 

  47. Xia T, Lai W, Han M, Han M, Ma X, Zhang L. Dietary ZnO nanoparticles alters intestinal microbiota and inflammation response in weaned piglets. Oncotarget. 2017;8(39):64878–91.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Hojberg O, Canibe N, Poulsen HD, Hedemann MS, Jensen BB. Influence of dietary zinc oxide and copper sulfate on the gastrointestinal ecosystem in newly weaned piglets. Appl Environ Microbiol. 2005;71(5):2267–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Starke IC, Pieper R, Neumann K, Zentek J, Vahjen W. The impact of high dietary zinc oxide on the development of the intestinal microbiota in weaned piglets. FEMS Microbiol Ecol. 2014;87(2):416–27.

    Article  CAS  PubMed  Google Scholar 

  50. Shen J, Chen Y, Wang Z, Zhou A, He M, Mao L, Zou H, Peng Q, Xue B, Wang L, et al. Coated zinc oxide improves intestinal immunity function and regulates microbiota composition in weaned piglets. Br J Nutr. 2014;111(12):2123–34.

    Article  CAS  PubMed  Google Scholar 

  51. Li BT, Van Kessel AG, Caine WR, Huang SX, Kirkwood RN. Small intestinal morphology and bacterial populations in ileal digesta and feces of newly weaned pigs receiving a high dietary level of zinc oxide. Can J Anim Sci. 2001;81(4):511–6.

    Article  CAS  Google Scholar 

  52. Yu Z, Morrison M. Improved extraction of PCR-quality community DNA from digesta and fecal samples. Biotechniques. 2004;36:808–12.

    Article  CAS  PubMed  Google Scholar 

  53. Edwards U, Rogall T, Blockerl H, Emde M, Bottger EC. Isolation and direct complete nucleotide determination of entire genes. Characterization of a gene coding for 16S ribosomal RNA. Nucleic Acids Res. 1989;17(19):7843–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Lane DJ, Pace B, Olsen GJ, Stahl DA, Sogin ML, Pace NR. Rapid determination of 16S ribosomal RNA sequences for phylogenetic analyses. Proc Natl Acad Sci. 1985;82:6955–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Opdahl LJ, Gonda MG, St-Pierre B. Identification of uncultured bacterial species from Firmicutes, Bacteroidetes and CANDIDATUS Saccharibacteria as candidate cellulose utilizers from the rumen of beef cows. Microorganisms. 2018;6:17.

    Article  PubMed Central  CAS  Google Scholar 

  56. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Parte AC. LPSN--list of prokaryotic names with standing in nomenclature. Nucleic Acids Res. 2014;42(Database issue):D613–6.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the graduate and undergraduate students who helped with the collection of performance data and samples during the animal trial.

Funding

Funding was provided by Vitech Bio-Chem Corporation (Glendale, CA, USA), and by the South Dakota State University Agricultural Experiment Station. Vitech Bio-Chem Corporation approved the design of the study, but did not participate in the collection, analysis, and interpretation of data, and it was not involved in writing the manuscript. The South Dakota State University Agricultural Experiment Station was not involved in any of the above mentioned steps.

Author information

Authors and Affiliations

Authors

Contributions

CL and BS designed the study. CL, RS, and BS performed the animal performance trial. PP performed the laboratory procedures and bioinformatics for microbiome analysis. PP, CL RS, and BS analyzed and interpreted the data. PP and BS drafted the manuscript. CL and RS critically revised the manuscript. The author(s) have read and approved the manuscript.

Corresponding author

Correspondence to Benoit St-Pierre.

Ethics declarations

Ethics approval and consent to participate

All animal trial procedures were approved by the South Dakota State University Institutional Animal Care and Use Committee before the start of the study (Protocol 17-035A).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Supplementary File 1.

Ingredient composition of experimental diets. List of ingredients and their respective proportion (%) in each nursery pig phase diet used in this study.

Additional file 2: Supplementary File 2.

. Operational Taxonomic Unit (OTU) table. Relative abundance for each OTU identified in fecal samples collected from nursery pigs fed a control (CON) or Peptiva-supplemented diet (PEP) at two different stages.

Additional file 3: Supplementary File 3.

P values associated with Table 3 and Table 5. Tables showing the P values obtained from the non-parametric Friedman test and the post hoc Nemenyi test for multiple pairwise comparisons that are presented in Table 3 and Table 5.

Additional file 4: Supplementary File 4.

Custom Perl script pipeline for 16S rRNA data analysis. List of Perl scripts and their respective code that were used for analysis of 16S rRNA data.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Poudel, P., Levesque, C., Samuel, R. et al. Dietary inclusion of Peptiva, a peptide-based feed additive, can accelerate the maturation of the fecal bacterial microbiome in weaned pigs. BMC Vet Res 16, 60 (2020). https://doi.org/10.1186/s12917-020-02282-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12917-020-02282-x

Keywords