Age-based dynamic changes of phylogenetic composition and interaction networks of health pig gut microbiome feeding in a uniformed condition
BMC Veterinary Research volume 15, Article number: 172 (2019)
The gut microbiota impacts on a range of host biological processes, and the imbalances in its composition are associated with pathology. Though the understanding of contribution of the many factors, e.g. gender, diet and age, in the development of gut microbiota has been well established, the dynamic changes of the phylogenetic composition and the interaction networks along with the age remain unclear in pigs.
Here we applied 16S ribosomal RNA gene sequencing, enterotype-like clustering (Classification of the gut microbiome into distinct types) and phylogenetic co-occurrence network to explore the dynamic changes of pig gut microbiome following the ages with a successive investigation at four ages in a cohort of 953 pigs. We found that Firmicutes and Bacteroidetes are two predominant phyla throughout the experimental period. The richness of gut microbiota was significantly increased from 25 to 240 days of age. Principal coordinates analysis showed a clear difference in the gut microbial community compositions between pre-weaning piglets and the pigs at the other three age groups. The gut microbiota of pre-weaning piglets was clearly classified into two enterotypes, which were dominated by Fusobacterium and p-75-a5, respectively. However, Prevotella and Treponema were the main drivers of the enterotypes for pigs at the age of 80, 120 and 240 days. Besides the piglets, even some adult pigs switched putative enterotypes between ages. We confirmed that the topological features of phylogenetic co-occurrence networks, including scale, stability and complexity were increased along with the age. The biological significance for modules in the network of piglets were mainly associated with the utilization of simple carbohydrate and lactose, whereas the sub-networks identified at the ages of 80, 120 and 240 days may be involved in the digestion of complex dietary polysaccharide. The modules related to the metabolism of protein and amino acids could be identified in the networks at 120 and 240 days. This dynamic change of the functional capacities of gut microbiome was further supported by functional prediction analysis.
The present study provided meaningful biological insights into the age-based dynamic shifts of ecological community of porcine gut microbiota.
The microbial community of the mammalian gastrointestinal tract harbors a complex and dynamic ecosystem, which is populated with as many as trillions of microbes, including viruses, bacteria, archaea, fungi and protists . Gut microbiota executes numerous vital functions for homeostasis ranging from harvesting essential nutrients, regulating host metabolism and maintaining the immune system . The establishment and maintenance of the gut microbial ecosystem are crucial to host health because they ensure the persistent occurrence of functions associated with those beneficial microbes . The gut microbiota can be influenced by diet, age, host genetics, and many environmental factors [4,5,6,7]. A previous study indicated clear differences in the composition of gut microbiota among infants, toddlers, adults and the elderly . The gut microbiota of newborns undergoes substantial modulation, in which facultative and strict anaerobes gradually dominate the community by replacing aerobes . The age-dependent shifts in gut microbiome reflect the change from the primary lactate metabolism to the enhancement of plant polysaccharide metabolism ability . The bacterial diversity of human gut microbiota increases with age and the phylogenetic composition evolves towards an adult-like configuration within the first three years after birth . The most noticeable feature in the gut microbial composition between young and elderly individuals is an alteration in the ratio of Firmicutes to Bacteroidetes, in which the young adults have a higher proportion of Firmicutes and the older adults show a higher proportion of Bacteroidetes . In pigs, previous reports indicated that the alpha-diversity and the taxonomic composition of intestinal microbiota were significantly altered with the growth of pigs [12,13,14,15]. A study in commercially raised pigs found a significant association between aging and an increasing measure of richness and diversity as well as distinct changes of the core microbiota .
The complexity of intestinal flora is reflected in both of community structure and function capacity due to their dynamic nature and compositional variability. Microbes in gastrointestinal tract form a complex ecological system with various symbiotic relationships, rather than a simplified collection of independent individuals . The interactions between or among organisms in a microecological system can be characterized into the relationships of mutualism, commensalism, amensalism, predation and competition . The co-occurrence patterns of microorganisms and the microbial relationships are essential for community assembly and stability , and can deduce the various effects on the host health. The application of system and network theory can facilitate the analyze of the composition and interaction of gut microbiota, enhancing our understanding of its complex ecological characteristics . Although the classical longitudinal analysis captures the changes of gut microbial composition in outcomes over time, the information about dynamic changes of the interaction networks between bacteria remains unknown. To our knowledge, there has been no systematic study about longitudinal dynamics of the phylogenetic composition and the interaction network of gut microbiota in pigs, especially, in an experimental pig cohort collected fecal samples at multiple ages. Furthermore, due to the similarity of the pig and human in gastrointestinal system and metabolism physiology , swine has been considered as one of the ideal models for studying the dynamic natures of gut microbiome.
In this study, we employed enterotype-like clustering analysis and co-occurrence network analysis to investigate the age-based dynamic shifts of pig gut microbiota in an experimental pig population collected the fecal samples at four different ages. Our results provided biological insights into the organization, function and evolution of pig gut microbial community following the age.
The experimental cohort was sourced from one pig farm of Jiangxi Agricultural University in Nanchang, Jiangxi Province and comprised of 953 pigs from a F6 population of heterogeneous pig cross, which had been constructed with eight founder breeds, including four Western commercial pig breeds of Landrace, Large White, Duroc and Pietran, and four Chinese indigenous pig breeds of Erhualian, Bamaxiang, Tibetan and Laiwu. All pigs were raised at the experimental pig farm of Jiangxi Agricultural University in Nanchang, Jiangxi. Piglets were housed with their mothers during suckling period and creep feed was provided during the last week of lactation (Additional file 1: Table S1). After the weaning at 28 days, all experimental pigs were transferred to a uniform fattening house and fed two times a day with formula diets (Additional file 1: Table S1). Water was available ad libitum from nipple drinkers. All boars were castrated at 80 days. All experimental animals used in this study were healthy and did not receive any antibiotic treatments within two months before fecal sample collection. The experimental pigs were slaughtered at a commercial slaughterhouse at 240 days of age by bleeding after electrical stunning.
Collection of fecal samples and DNA extraction
Fecal samples were collected at the ages of 25, 120 and 240 days, which represented the time of preweaning, mid-stage of fattening and slaughtering. A total of 1417 fecal samples were harvested from the rectum of experimental pigs in the cohort, including 175 samples collected at the age of 25 days (from 92 females and 83 entire males), 551 samples at the age of 120 days (from 302 females and 249 castrated males) and 691 samples at the age of 240 days (from 377 females and 314 castrated males). In another study, to investigate the effect of castration on the phylogenetic composition of gut microbiota, we collected additional 64 fecal samples (from 36 females and 28 entire males) in the same experimental cohort at the age of 80 days (boars were castrated immediately after collected the fecal samples). To analyze dynamic changes of gut microbiota at more age times, the bacterial composition data of these 64 samples were also included in this study. The detailed information about sampling is shown in the Additional file 2: Figure S1.
All samples were collected in sterile plastic centrifuge tubes and deep frozen in liquid nitrogen. After transported to the laboratory, the samples were stored at − 80 °C freezer until use. DNA extraction from the fecal samples was performed using QIAamp DNA Stool Mini Kit (Qiagen, Germany) following manufacturer’s instructions (McOrist et al., 2002). The quality and integrity of the DNA samples were checked by electrophoresis with 0.8% agarose gel, and DNA concentrations were measured with a ND-1000 spectrophotometer (Nanodrop Technologies, USA).
Amplification and sequencing of bacterial 16S rRNA gene
The V3-V4 region of 16S rRNA gene was amplified by PCR with the barcode fusion forward primer 338F [ACTCCTACGGGAGGCAGCAG] and the reverse primer 806R [GGACTACHVGGGTWTCTAAT] under the melting temperature of 55 °C with 28 cycles. For the bacterial composition data of 64 samples harvested at the age of 80 days mentioned above, only V4 hypervariable region of 16S rRNA gene was amplified by the fusion primers 515F [GTGCCAGCMGCCGCGGTAA] and 806R [GGACTACHVGGGTWTCTAAT] under the melting temperature of 56 °C with 30 cycles. The sequencing of the PCR amplicons was performed with a 250-bp paired-end procedure on an Illumina MiSeq platform (Illumina, USA) according to the manufacturer’s manuals.
Bioinformatics processing of sequence reads
To obtain the clean sequence reads, the primer and barcode sequences, and the low-quality reads were excluded from further analysis. FLASH (version 1.2.11) was used to assemble the paired-end clean reads into tags . To avoid the bias generated by the different sequencing depth , we rarefied the library size to 10,000 tags per sample. Unique bacterial sequences with 97% of the sequence similarity were clustered into operational taxonomic unit (OTU) using the QIIME software (the toolbox for Quantitative Insights Into Microbial Ecology) , which uses UCLUST (an algorithm to cluster sequence reads based on similarity) to perform the clustering . OTUs were matched to bacteria by using a primer-specific version of the GreenGenes reference database (version 13.5) . We filtered out those OTUs which had relative abundance less than 0.01% and were presented in less than 1% of the experimental pigs from further analyses. The alpha-diversity indexes of chao1, ACE, observed species, Simpson and Shannon were calculated by Mothur software (Version 1.42.0) . The comparison of the relative abundances of bacterial taxa and the alpha-diversity indexes between two age groups was performed by Wilcoxon t-test. The P values were adjusted the multiple tests by Benjamini-Hochberg (BH) method with the threshold of false discovery rate (FDR) ≤ 5%. In order to evaluate the effects of host age and sex on microbial composition of fecal samples, a Bray-Curtis similarity matrix was constructed based on the relative abundances of OTUs. We split the sequences of V3-V4 region of 16S rRNA gene into the sequences containing only V4 region for the samples at 25, 120 and 240 days of age. And then, all sequence data of the V4 region were incorporated together for further calculation of Bray-Curtis matrix. Principal coordinated analysis (PCoA) of Bray-Curtis distances was performed using the R function pcoa.
The enterotype-like clustering analysis of fecal microbiota at four ages was separately performed with the method described previously . Briefly, we calculated Jensen-Shannon divergence (JSD) distance for the relative abundances of bacterial taxa at the genus level by applying the Partitioning Around Medoids (PAM) method. The optimal number of clusters and the groups’ robustness were evaluated with CalinskiHarabasz (CH) index and silhouette value, respectively. To identify the genera with different abundances between the two enterotypes, a linear discriminant analysis (LDA) effect size (LEfSe) analysis was performed under the condition ɑ = 0.01, with an LDA score of at least 3 . Sparse Correlations for Compositional data (SparCC) was employed to determine co-abundance (positive) and co-exclusion (negative) relationships between genera based on their relative abundances . Significant correlations between bacterial genera were detected using the partial correlation and information theory (PCIT) algorithm . The main absolute correlations (top 15%) were transformed into links between two genera in the genus network and the networks were visualized in Cytoscape (version 3.4.0).
Construction of phylogenetic interaction networks
We filtered out those OTUs which had less than 0.05% of relative abundance and were presented in less than 10% of the tested samples from further interaction network analysis. Phylogenetic co-occurrence networks of porcine fecal microbiota at the ages of 25, 80, 120 and 240 days were separately inferred based on the SparCC algorithm . Significant correlations due to the compositional structures of OTUs were detected using the PCIT algorithm . The confidence of the interactions between nodes was established with > 0.55 of absolute sparse correlation coefficient. We used Cytoscape (version 3.4.0) to visualize co-occurrence networks and calculate their topological characteristics including densities, clustering coefficients and scale-free properties . The clustering of sub-modules was based on the vertex weighting by local neighborhood density and outward traversal from a locally dense seed node according to given parameters. Assessment of biologically important modules in OTU networks was performed using Molecular Complex Detection (MCODE) plugin in the Cytoscape software . To further compare the topological characteristics of main modules in the interaction networks, we used Network Analysis Profler (NAP) to numerically calculate the topological features and metrics of networks (e.g. average eccentricity, average number of neighbors, centralization betweenness, centralization closeness and centralization degree) .
Prediction and comparison of functional capacities of fecal microbiome among four age groups
To compare the potential function capacities of fecal microbiome among four age groups, we used the phylogenetic investigation of communities by reconstruction of unobserved states (PICRUST) algorithm to obtain the function profiles of bacterial community . Functional genes were categorized into KEGG pathways. The relative abundance of each subclass term of the KEGG pathways was calculated by summing the abundances of functional genes that were annotated to the functional subsystem. Tukey-Kramer post-hoc test was used for pairwise comparisons of potential function capacities of fecal microbiomes among four age groups. Story’s FDR was used to correct the multiple tests. All tests were performed using Statistical Analysis of Metagenomic Profiles (STAMP) software .
Longitudinal dynamics of fecal microbial composition and richness
The total numbers of high-quality sequence tags for samples at the ages of 25, 80, 120 and 240 days were 5,746,871 (an average of 32,839 tags per sample), 1,082,852 (an average of 16,919 tags per sample), 16,825,154 (an average of 30,536 tags per sample) and 19,749,079 (an average of 28,580 tags per sample), respectively. A total of 3392, 2344, 3401 and 3609 OTUs were obtained. To investigate the dynamic shifts of fecal microbiota structure along with the age, the alpha-diversity of fecal microbiota was compared among four ages using the Chao1, ACE, observed species, Shannon and Simpson index. Interestingly, the experimental pigs showed the continuously increased indexes of Chao1, ACE and observed species following the age (Fig. 1A-C). Shannon index, which reflects the richness and evenness of bacterial taxa, was increased markedly from the age of 25 to 120 days, but showed no significant difference between the ages of 120 and 240 days (Fig. 1D). The Simpson index was significantly decreased from 25 to 120 days, but slightly increased from 120 to 240 days (Fig. 1E). This result suggested that the richness of gut microbiota was continuously increased from 25 to 240 days, but the highest diversity was achieved at 120 days. PCoA analysis based on the relative abundances of OTUs revealed the distinctly microbial compositions between pre-weaning piglets and the pigs at the other three ages (Fig. 1F). The bacterial composition of fecal samples at the age of 80 days was more similar to that at 120 and 240 days. However, no distinct difference was observed in fecal bacterial structures between the ages of 120 and 240 days.
We further evaluated the age-based dynamic changes of microbial composition at the taxonomic level. The phylogenetic compositions of fecal microbiota at the phylum level at four ages are shown in Fig. 2A. Firmicutes was the most predominant phylum at all four age stages, but had slightly less abundance in piglets (25 days) compared with that in pigs at the other three ages (P > 0.05). The relative abundance of Bacteroidetes was increased from the age of 25 to 80 days (P < 0.0001), but subsequently decreased after 80 days (P < 0.01). The relative abundances of Proteobacteria and Fusobacteria were decreased dramatically after the age of 25 days (P < 0.0001). At the genus level, the fecal microbiota of piglets was dominated by Bacteroides, Fusobacterium, p-75-a5 and Prevotella, whereas Lactobacillus, Prevotella, Ruminococcus and Treponema were the core bacteria for fecal samples at 80, 120 and 240 days of age (Fig. 2B). The relative abundances of three predominant genera in pre-weaning pigs, including Bacteroides, Fusobacterium and p-75-a5 were significantly decreased after 25 days of age (P < 0.0001; Fig. 2B). The relative abundances of Lactobacillus, Prevotella and Ruminococcus were significantly increased from 25 to 80 days of age (P < 0.0001), but subsequently decreased from the age of 80 to 120 days (P < 0.05). However, there was no significant difference between 120 and 240 days. The proportion of Treponema was significantly increased from 25 to 120 days (P < 0.05), but showed no significant difference between 120 and 240 days (P > 0.05).
Besides age, we also evaluated the effect of host sex on gut microbial community structures at all four age stages. Although the samples harvested at 80 days of age showed a slight difference in bacterial compositions between gilts and entire boars, we did not observe distinct sex-biased bacterial compositions for samples at the other three age stages in PCoA analysis (Additional file 3: Figure S2). This should be due to the reasons of sexual immuration for piglets at the age of 25 days and the castration for boars at 120 and 240 days. However, this discrepancy may also be caused by the different amplicons in 16S rRNA gene sequencing analysis for samples at 80 days.
Longitudinal dynamics of enterotype-like clustering of porcine fecal microbiota
We separately performed the enterotype-like clustering analysis at each of four ages. As the result, two enterotype-like clusters were identified for each age. Distinct from the enterotypes for pre-weaning piglets, which were dominated by p-75-a5 and Fusobacterium, respectively (enterotype 1 and 2), the enterotype clusters for the pigs at 80, 120 and 240 days showed a high degree of similarity and could be grouped into the other two enterotypes that were overrepresented by Prevotella and Treponema, respectively (enterotype 3 and 4) (Fig. 3). The differential taxa of bacteria between the two enterotypes at each age are shown in Additional file 4: Figure S3. To further investigate the discrepancy of the enterotypes between pre-weaning piglets and the pigs at the other three ages, we constructed the SparCC network at the genus level in all four age groups (Additional file 5: Figure S4). We observed that Fusobacterium and p-75-a5 were not only the main drivers for the enterotypes at 25 days (Fig. 3A), but also the hub nodes in the network (Additional file 5: Figure S4). However, Treponema was also overrepresented in the enterotype 1, whereas Prevotella was more abundant in the enterotype 2 (Additional file 4: Fig. S3), and these two genera were the other two hubs of the network at 25 days. Furthermore, Fusobacterium showed co-abundance with Prevotella and co-exclusion with Treponema. This result suggested that the gut microbiota in pre-weaning piglets should not change completely with age, but only increased the diversity. Some bacterial taxa gradually dominated the bacterial composition by replacing the initial predominant taxa following the age.
To evaluate the stability of enterotypes following the age, we particularly checked the enterotype distribution in pigs with fecal samples at more than one age. Of the 57 pigs having fecal samples at both 25 and 120 days (Additional file 2: Fig. S1), twenty-three (40.35%) were assigned to the Fusobacterium enterotype, and the others (59.65%) were occupied by p-75-a5 enterotype at 25 days of age. However, at the age of 120 days, among those pigs with Fusobacterium enterotype at 25 days of age, sixteen (69.57%) switched their putative enterotypes to Prevotella, and the others (30.43%) changed their enterotypes to Treponema. A similar result was also obtained in pre-weaning pigs with p-75-a5 enterotype. Twenty-six (76.47%) out of these piglets shifted the enterotypes to Prevotella, whereas the left 23.53% of piglets switched the enterotypes to Treponema at 120 days. There were 377 pigs involved in the enterotype analysis at both 120 and 240 days (Additional file 2: Figure S1). Eighty-four (22.3%) pigs with Prevotella enterotype at 120 days switched their putative enterotypes to Treponema at 240 days. Inversely, 142 (37.7%) pigs having Treponema enterotype changed the enterotypes to Prevotella at 240 days.
Dynamic changes of interaction networks of porcine fecal bacteria with age
To evaluate the dynamic changes of potential interaction networks among bacterial taxa at different ages, 273, 334, 351 and 322 key OTUs at the age of 25, 80, 120 and 240 days were used for constructing the phylogenetic interaction networks, respectively. Overall, in agreement with the previous report , OTUs exhibited a various degree of connectivity at each age. The complexity of the phylogenetic interaction networks reflecting in the average number of edges per node was increased from 25 (1.48) to 120 days (2.4), but no significant difference between 120 and 240 days (2.29). The stability (i.e. percentage of negative interactions) was obviously increased from 1.0% at 25 days to 16.9% at 120 days. However, this percentage was decreased to 8.9% at 240 days. The highest percentage of negative interactions was identified at 80 days (21.4%). This continuously declined percentage from the ages of 80 and 120 days to 240 days suggested that the highest stability of gut microbial ecology was achieved before 120 days. However, we could not exclude the influence of the smaller sample size and the different amplicon in 16S rRNA gene sequencing on the percentage of negative interactions for samples at 80 days. In more details, the interaction network of the gut microbiota at the age of 25 days was comprised of 66 nodes and 98 edges (Fig. 4A). This phylogenetic interaction network was mostly divided into three distinct modules, including two piglet-specific modules (module 1 and 3). The nodes in the module 1 were mainly annotated to Leptonema, PSB-M-3, Shewanella, Pseudomonas and Sulfurospirillum (Additional file 6: Table S2). These bacterial genera are mainly involved in the utilization of carbohydrate in diet milk [38,39,40,41]. The module 2 was comprised of Sphaerochaeta (Additional file 6: Table S2) which can produce a diverse set of saccharolytic enzymes that participate in the glycolytic and pentose phosphate pathways . Based on the overlapping structures and topological features (Fig. 5 and Additional file 7: Figure S5), this module was also presented at 120 (module 3) and 240 (module 3) days of age. Module 3 contained the bacteria that mainly takes part in the utilization of the lactose and glutamic acid (Additional file 6: Table S2) [43,44,45,46]. The interaction network at the age of 80 days was comprised of 129 nodes and 350 interactions. Three modules were identified in this network (Fig. 4B). Most of the nodes in this co-occurrence network belonged to the genera Blautia, Faecalibacterium, Prevotella, Ruminococcus and Treponema (Additional file 6: Table S2). These microbes are mainly involved in the digestion of dietary polysaccharides (e.g. fiber, cellulose and lignin) and the production of significant amounts of short chain fatty acids (SCFAs) [5, 47,48,49].
The interaction network at the age of 120 days was comprised of 150 nodes and 360 edges, and mainly divided into three modules (Fig. 4C). Based on the topological features, the module 1 and 3 were highly similar to the module 2 at the age of 240 days and the module 2 at the age of 25 days, respectively (Fig. 5 and Additional file 7: Figure S5). The module 2 consisted of Prevotella and Bacteroides (Additional file 6: Table S2), which play essential roles in the process of the degradation of complex dietary polysaccharides and amino acid metabolism [47, 50, 51]. The interaction network at the age of 240 days was comprised of 142 nodes and 325 edges. Three distinct modules were identified in this co-occurrence network (Fig. 4D), and only the module 1 which was primarily comprised of the bacteria related to dietary protein metabolism (e.g. Cetobacterium somerae, Coprococcus eutactus, Mycobacterium, CF231, Comamonas, Ralstonia, Sporosarcina and Nitrincola) was unique to this network [52,53,54,55,56,57,58,59]. The modules 2 and 3 were highly similar to the modules identified at 120 (module 1) and 25 days (module 2), respectively.
Longitudinal dynamics of predicted function capacity of fecal microbiome
To investigate the dynamic changes of the potential function capacity of fecal microbiome with age, the relative abundances of KEGG pathways were predicted by PICRUST based on 16S rRNA gene sequences. A total of 39 KEGG pathways showed significantly different enrichments at different ages (P < 0.05, Fig. 6). Eleven out of these 39 KEGG pathways were significantly enriched in pre-weaning piglets (Additional file 8: Table S3), including carbohydrate metabolism, energy metabolism, and xenobiotics biodegradation and metabolism. There were four KEGG pathways showing significant enrichments at the age of 80 days, including enzyme families, metabolism, nervous system and transcription (P < 0.001). We also detected one (folding, sorting and degradation, P < 0.01) and four (amino acid metabolism, endocrine system, excretory system, and transport and catabolism, P < 0.001) pathways that were significantly enriched at 120 and 240 days, respectively. Furthermore, the relative abundances of cell motility and amino acid metabolism were significantly increased following the age (P < 0.001), whereas the relative abundance of membrane transport was decreased with age (P < 0.001).
Several studies have reported the age-related gut microbial shifts in pigs [16, 60,61,62]. However, to our knowledge, there are few studies about the dynamic shifts of enterotype-like clusters and interaction networks of gut microbiota following the age in an experimental pig cohort with large sample size. In this study, we explored the longitudinal changes of gut microbial composition, enterotypes, co-occurrence interaction network and potential function capacity following the age in an experimental cohort comprised of 953 pigs. Especially, most of the experimental pigs were continuously collected fecal samples at more than one age. Our study provided an integrated overview of longitudinal changes of the porcine gut microbiome.
In accordance with the previous studies [63, 64], the dominant bacterial phyla in porcine fecal microbiota throughout the experimental period were Firmicutes and Bacteroidetes. The age-based dynamic changes of the relative abundance of Bacteroidetes were mainly contributed by Bacteroides and Prevotella. The decreased percentage of Fusobacteria was mainly caused by reduced abundance of Fusobacterium, while the variation of Spirochaetes was a result of increased abundance of Treponema. In agreement with the previous reports [16, 62, 65, 66], Proteobacteria was relatively more abundant in the gut of piglets, but its relative abundance was considerably decreased in the adulthood. The relative abundances of opportunistic pathogens (e.g. Actinobacillus and Helicobacter) might be reduced following the maturation of pig gut microbiota. It has also been demonstrated that a high diversity of gut microbiota is generally considered related to the host health and fitness, and it is also regard as the sign of the maturation of gut microbiota [67, 68]. We found that the richness of fecal microbiota was continuously increased following the age in this study.
Consistent with the findings from previous report , two enterotype-like clusters were identified at each age. Fusobacterium and p-75-a5 were the main drivers of the enterotypes in pre-weaning piglet. The genus p-75-a5 was a member within the family Erysipelotrichaceae, which was highly correlated with the digestion of protein and fat . As we have well known, milk contained high concentration of protein and fat is the main diet for piglets. It has been reported that Fusobacterium was associated with some diseases in animals , indicating that opportunistic pathogenic bacteria were commonly presented in suckling piglets. Our data coincides with the fact that the piglet has an immature immune system and thus is unable to manage and control the parasitic invaders . March et al. reported that Ruminococcaceae was the main driver for the only enterotype-like cluster in the 14-day-old suckling pigs . This discrepancy may be caused by the different environments, diets and pig breed [72, 73]. Prevotella and Treponema were the main drivers of the enterotypes at 80, 120 and 240 days. The previous report also found that Prevotella and Treponema were the main drivers for the enterotypes of pigs at 60 days of age . Prevotella is capable of metabolizing dietary polysaccharide and producing large amounts of SCFAs . Treponema plays an essential role in cellulose and lignin degradation . This result suggested that the gut microbiota was mature when the age of pigs was above 80 days, and its ability digesting diet polysaccharide was improved. The enterotypes in pre-weaning pigs did not remain stable (see Results). Surprisingly, the enterotypes did not keep stable even in some adult pigs, e.g. more than 22.3% of pigs with Prevotella enterotype at 120 days showed Treponema enterotype at 240 days. It has been reported that Prevotella enterotype may be better adapted to the diet for growing pigs which contains more plant polysaccharides . In this study, the diets provided to the pigs at 120 and 240 days of age contained slight difference of coarse fiber levels (increased from 5 to 8%) (Additional file 1: Table S1). Costea et al. also reported that although all three enterotypes of human gut microbiota kept stable in overall, 16% of individuals switched putative enterotypes between visits .
The scale, complexity and stability of the phylogenetic co-occurrence network at 25 days of age were lower than that at the other three ages, but the networks at 120 and 240 days were highly similar. This was consistent with the result of enterotype analysis, which suggested an unstable structure of fecal microbial community in pre-weaning piglets. The gut microbial ecology would remain stable once the gut microbiota achieves maturation although it may be fluctuations in populations . The main modules in the co-occurrence network of piglets were associated with the utilization of simple carbohydrate (module 1) and lactose (module 3). It is well known that sow milk was the main source of nutrients for suckling piglets. The module commonly identified at 25 (module 2), 120 (module 3) and 240 (module 3) days of age was comprised of Sphaerochaeta which has been reported to take part in glucose metabolism through glycolytic and pentose phosphate pathways [42, 76]. This result suggested that glycolytic and pentose phosphate pathways of gut microbiome were essential for pigs throughout the whole-life span. However, this module was absent in the network at 80 days. This discrepancy may be caused by the smaller sample size and the different sequencing region of 16S rRNA gene. Several modules identified at 80, 120 (module 1 and 2) and 240 (module 2) days were related to the digestion of nondigestible carbohydrates and the production of SCFAs [5, 49], suggesting that the gut microbiota above the age of 80 days obtained the functional capacity of digesting dietary fiber. The modules involved in the metabolism of dietary protein and amino acids were identified at 120 (module 2) and 240 (module 1) days, indicating that the functional capacity of gut microbiome digesting dietary protein had been significantly improved since the age of 120 days.
Consistent with the result observed in the interaction network analysis, the KEGG pathways enriched in piglets were related to carbohydrate and energy metabolism. This may be due to the high fat and carbohydrate concentration of sow milk diet. The functional capacities of cell motility and amino acid metabolism were increased with ages. Cellular motility is essential for biological processes, such as the development and immune response . As mentioned above, the abilities of immune response and amino acid metabolism of gut microbiome were increased following the age.
The present study demonstrated the dynamic changes of the phylogenetic composition, enterotypes, interaction networks and potential functional capacity of porcine gut microbiome following the age. These longitudinal changes of gut microbiome contribute to the metabolism of diet nutrients and the maturation of immune system in pigs. The results provide meaningful biological insights into the age-based dynamic shifts of ecological community of porcine gut microbiota and provide valuable knowledge about establishment and maturation of gut microbiome.
Kyoto encyclopedia of genes and genomes
Operational taxonomic units
Ribosomal robonucleic acid
Short-chain fatty acid
Gill SR, Pop M, Deboy RT, Eckburg PB, Turnbaugh PJ, Samuel BS, Gordon JI, Relman DA, Fraser-Liggett CM, Nelson KE. Metagenomic analysis of the human distal gut microbiome. Science. 2006;312(5778):1355–9.
Nicholson JK, Holmes E, Kinross J, Burcelin R, Gibson G, Jia W, Pettersson S. Host-gut microbiota metabolic interactions. Science. 2012;336(6086):1262–7.
Lozupone CA, Stombaugh JI, Gordon JI, Jansson JK, Knight R. Diversity, stability and resilience of the human gut microbiota. Nature. 2012;489(7415):220–30.
Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, Magris M, Hidalgo G, Baldassano RN, Anokhin AP, et al. Human gut microbiome viewed across age and geography. Nature. 2012;486(7402):222–7.
De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, Collini S, Pieraccini G, Lionetti P. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc Natl Acad Sci U S A. 2010;107(33):14691–6.
O'Toole PW, Jeffery IB. Gut microbiota and aging. Science. 2015;350(6265):1214–5.
Clark RI, Salazar A, Yamada R, Fitz-Gibbon S, Morselli M, Alcaraz J, Rana A, Rera M, Pellegrini M, Ja WW, et al. Distinct shifts in microbiota composition during Drosophila aging impair intestinal function and drive mortality. Cell Rep. 2015;12(10):1656–67.
Conlon MA, Bird AR. The impact of diet and lifestyle on gut microbiota and human health. Nutrients. 2014;7(1):17–44.
Costello EK, Stagaman K, Dethlefsen L, Bohannan BJ, Relman DA. The application of ecological theory toward an understanding of the human microbiome. Science. 2012;336(6086):1255–62.
Langille MG, Meehan CJ, Koenig JE, Dhanani AS, Rose RA, Howlett SE, Beiko RG. Microbial shifts in the aging mouse gut. Microbiome. 2014;2(1):50.
Mariat D, Firmesse O, Levenez F, Guimaraes V, Sokol H, Dore J, Corthier G, Furet JP. The Firmicutes/Bacteroidetes ratio of the human microbiota changes with age. BMC Microbiol. 2009;9:123.
Chen L, Xu Y, Chen X, Fang C, Zhao L, Chen F. The maturing development of gut microbiota in commercial piglets during the weaning transition. Front Microbiol. 2017;8:1688.
Kim J, Nguyen SG, Guevarra RB, Lee I, Unno T. Analysis of swine fecal microbiota at various growth stages. Arch Microbiol. 2015;197(6):753–9.
Han GG, Lee JY, Jin GD, Park J, Choi YH, Kang SK, Chae BJ, Kim EB, Choi YJ. Tracing of the fecal microbiota of commercial pigs at five growth stages from birth to shipment. Sci Rep. 2018;8(1):6012.
Hu J, Nie Y, Chen J, Zhang Y, Wang Z, Fan Q, Yan X. Gradual changes of gut microbiota in weaned miniature piglets. Front Microbiol. 2016;7:1727.
Slifierz MJ, Friendship RM, Weese JS. Longitudinal study of the early-life fecal and nasal microbiotas of the domestic pig. BMC Microbiol. 2015;15(1):184.
Young VB, Schmidt TM. Overview of the gastrointestinal microbiota. Adv Exp Med Biol. 2008;635:29–40.
Faust K, Raes J. Microbial interactions: from networks to models. Nat Rev Microbiol. 2012;10(8):538–50.
Faust K, Sathirapongsasuti JF, Izard J, Segata N, Gevers D, Raes J, Huttenhower C. Microbial co-occurrence relationships in the human microbiome. PLoS Comput Biol. 2012;8(7):e1002606.
Layeghifard M, Hwang DM, Guttman DS. Disentangling interactions in the microbiome: a network perspective. Trends Microbiol. 2017;25(3):217–28.
Wang M, Donovan SM. Human microbiota-associated swine: current progress and future opportunities. ILAR J. 2015;56(1):63–73.
Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63.
Hughes JB, Hellmann JJ. The application of rarefaction techniques to molecular inventories of microbial diversity. Methods Enzymol. 2005;397:292–308.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1.
DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72(7):5069–72.
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.
Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, Fernandes GR, Tap J, Bruls T, Batto JM, et al. Enterotypes of the human gut microbiome. Nature. 2011;473(7346):174–80.
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60.
Friedman J, Alm EJ. Inferring correlation networks from genomic survey data. PLoS Comput Biol. 2012;8(9):e1002687.
Reverter A, Chan EK. Combining partial correlation and an information theory approach to the reversed engineering of gene co-expression networks. Bioinformatics. 2008;24(21):2491–7.
Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011;27(3):431–2.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC bioinformatics. 2003;4(2).
Theodosiou T, Efstathiou G, Papanikolaou N, Kyrpides NC, Bagos PG, Iliopoulos I, Pavlopoulos GA. NAP: the network analysis profiler, a web tool for easier topological analysis and comparison of medium-scale biological networks. BMC Res Notes. 2017;10(1):278.
Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Vega Thurber RL, Knight R, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31(9):814–21.
Parks DH, Tyson GW, Hugenholtz P, Beiko RG. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics. 2014;30(21):3123–4.
Ramayo-Caldas Y, Mach N, Lepage P, Levenez F, Denis C, Lemonnier G, Leplat JJ, Billon Y, Berri M, Dore J, et al. Phylogenetic network analysis applied to pig gut microbiota identifies an ecosystem structure linked with growth traits. ISME J. 2016;10(12):2973–7.
Achouak W, Sutra L, Heulin T, Meyer JM, Fromin N, Degraeve S, Christen R, Gardan L. Pseudomonas brassicacearum sp. nov. and Pseudomonas thivervalensis sp. nov., two root-associated bacteria isolated from Brassica napus and Arabidopsis thaliana. Int J Syst Evol Microbiol. 2000;50(Pt 1):9–18.
Rodionov DA, Yang C, Li X, Rodionova IA, Wang Y, Obraztsova AY, Zagnitko OP, Overbeek R, Romine MF, Reed S, et al. Genomic encyclopedia of sugar utilization pathways in the Shewanella genus. BMC Genomics. 2010;11:494.
Huntemann M, Stackebrandt E, Held B, Nolan M, Lucas S, Hammon N, Deshpande S, Cheng JF, Tapia R, Goodwin LA, et al. Genome sequence of the phylogenetically isolated spirochete Leptonema illini type strain (3055(T)). Stand Genomic Sci. 2013;8(2):177–87.
Miura T, Uchino Y, Tsuchikane K, Ohtsubo Y, Ohji S, Hosoyama A, Ito M, Takahata Y, Yamazoe A, Suzuki K, et al. Complete genome sequences of Sulfurospirillum strains UCH001 and UCH003 isolated from groundwater in Japan. Genome announcements. 2015;3(2).
Caro-Quintero A, Ritalahti KM, Cusick KD, Loffler FE, Konstantinidis KT. The chimeric genome of Sphaerochaeta: nonspiral spirochetes that break with the prevalent dogma in spirochete biology. MBio. 2012;3(3).
Downes J, Olsvik B, Hiom SJ, Spratt DA, Cheeseman SL, Olsen I, Weightman AJ, Wade WG. Bulleidia extructa gen. Nov., sp. nov., isolated from the oral cavity. Int J Syst Evol Microbiol. 2000;50(Pt 3):979–83.
Jumas-Bilak E, Carlier JP, Jean-Pierre H, Mory F, Teyssier C, Gay B, Campos J, Marchandin H. Acidaminococcus intestini sp. nov., isolated from human clinical samples. Int J Syst Evol Microbiol. 2007;57(Pt 10:2314–9.
Lee CC, Lo WC, Lai SM, Chen YP, Tang CY, Lyu PC. Metabolic classification of microbial genomes using functional probes. BMC Genomics. 2012;13:157.
Tian L, Bruggeman G, van den Berg M, Borewicz K, Scheurink AJ, Bruininx E, de Vos P, Smidt H, Schols HA, Gruppen H. Effects of pectin on fermentation characteristics, carbohydrate utilization, and microbial community composition in the gastrointestinal tract of weaning pigs. Mol Nutr Food Res. 2017;61(1).
Ellekilde M, Selfjord E, Larsen CS, Jakesevic M, Rune I, Tranberg B, Vogensen FK, Nielsen DS, Bahl MI, Licht TR, et al. Transfer of gut microbiota from lean and obese mice to antibiotic-treated mice. Sci Rep. 2014;4:5922.
Niu Q, Li P, Hao S, Zhang Y, Kim SW, Li H, Ma X, Gao S, He L, Wu W, et al. Dynamic distribution of the gut microbiota and the relationship with apparent crude fiber digestibility and growth stages in pigs. Sci Rep. 2015;5:9938.
Claesson MJ, Jeffery IB, Conde S, Power SE, O'Connor EM, Cusack S, Harris HM, Coakley M, Lakshminarayanan B, O'Sullivan O, et al. Gut microbiota composition correlates with diet and health in the elderly. Nature. 2012;488(7410):178–84.
Chen T, Long W, Zhang C, Liu S, Zhao L, Hamaker BR. Fiber-utilizing capacity varies in Prevotella- versus Bacteroides-dominated gut microbiota. Sci Rep. 2017;7(1):2594.
Rios-Covian D, Salazar N, Gueimonde M, de Los Reyes-Gavilan CG. shaping the metabolism of intestinal Bacteroides population through diet to improve human health. Front Microbiol. 2017;8:376.
Zhao L, Meng Q, Ren L, Liu W, Zhang X, Huo Y, Zhou Z. Effects of nitrate addition on rumen fermentation, bacterial biodiversity and abundance. Asian Australas J Anim Sci. 2015;28(10):1433–41.
Ma YF, Zhang Y, Zhang JY, Chen DW, Zhu Y, Zheng H, Wang SY, Jiang CY, Zhao GP, Liu SJ. The complete genome of Comamonas testosteroni reveals its genetic adaptations to changing environments. Appl Environ Microbiol. 2009;75(21):6812–9.
Reichardt N, Duncan SH, Young P, Belenguer A, McWilliam Leitch C, Scott KP, Flint HJ, Louis P. Phylogenetic distribution of three pathways for propionate production within the human gut microbiota. The ISME journal. 2014;8(6):1323–35.
Finegold SM, Vaisanen ML, Molitoris DR, Tomzynski TJ, Song Y, Liu C, Collins MD, Lawson PA. Cetobacterium somerae sp. nov. from human feces and emended description of the genus Cetobacterium. Syst Appl Microbiol. 2003;26(2):177–81.
Cook GM, Berney M, Gebhard S, Heinemann M, Cox RA, Danilchanka O, Niederweis M. Physiology of mycobacteria. Adv Microb Physiol. 2009;55:81–182 318-189.
Zhang Y, Qiu S. Phylogenomic analysis of the genus Ralstonia based on 686 single-copy genes. Antonie Van Leeuwenhoek. 2016;109(1):71–82.
Joshi A, Thite S, Kulkarni G, Dhotre D, Joseph N, Venkata Ramana V, Polkade A, Shouche Y. Nitrincola alkalisediminis sp. nov., an alkaliphilic bacterium isolated from an alkaline lake. Int J Syst Evol Microbiol. 2016;66(3):1254–9.
McCoy DD, Cetin A, Hausinger RP. Characterization of urease from Sporosarcina ureae. Arch Microbiol. 1992;157(5):411–6.
Kim HB, Borewicz K, White BA, Singer RS, Sreevatsan S, Tu ZJ, Isaacson RE. Longitudinal investigation of the age-related bacterial diversity in the feces of commercial pigs. Vet Microbiol. 2011;153(1–2):124–33.
Mach N, Berri M, Estelle J, Levenez F, Lemonnier G, Denis C, Leplat JJ, Chevaleyre C, Billon Y, Dore J, et al. Early-life establishment of the swine gut microbiome and impact on host phenotypes. Environ Microbiol Rep. 2015;7(3):554–69.
Alain BPE, Chae JP, Balolong MP, Bum Kim H, Kang DK. Assessment of fecal bacterial diversity among healthy piglets during the weaning transition. J Gen Appl Microbiol. 2014;60(4):140–6.
Isaacson R, Kim HB. The intestinal microbiome of the pig. Anim Health Res Rev. 2012;13(1):100–9.
Looft T, Johnson TA, Allen HK, Bayles DO, Alt DP, Stedtfeld RD, Sul WJ, Stedtfeld TM, Chai B, Cole JR, et al. In-feed antibiotic effects on the swine intestinal microbiome. Proc Natl Acad Sci U S A. 2012;109(5):1691–6.
Hansen R, Thomson JM, Fox JG, El-Omar EM, Hold GL. Could helicobacter organisms cause inflammatory bowel disease? FEMS Immunol Med Microbiol. 2011;61(1):1–14.
Thompson AB, Postey RC, Snider T, Pasma T. Actinobacillus equuli as a primary pathogen in breeding sows and piglets. The Canadian veterinary journal = La revue veterinaire canadienne. 2010;51(11):1223–5.
Le Chatelier E, Nielsen T, Qin J, Prifti E, Hildebrand F, Falony G, Almeida M, Arumugam M, Batto JM, Kennedy S, et al. Richness of human gut microbiome correlates with metabolic markers. Nature. 2013;500(7464):541–6.
Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457(7228):480–4.
Bermingham EN, Maclean P, Thomas DG, Cave NJ, Young W. Key bacterial families (Clostridiaceae, Erysipelotrichaceae and Bacteroidaceae) are related to the digestion of protein and energy in dogs. PeerJ. 2017;5:e3019.
De Witte C, Flahou B, Ducatelle R, Smet A, De Bruyne E, Cnockaert M, Taminiau B, Daube G, Vandamme P, Haesebrouck F. Detection, isolation and characterization of fusobacterium gastrosuis sp. nov. colonizing the stomach of pigs. Syst Appl Microbiol. 2017;40(1):42–50.
Stokes CR. The development and role of microbial-host interactions in gut mucosal immune development. Journal of animal science and biotechnology. 2017;8:12.
Thompson CL, Wang B, Holmes AJ. The immediate environment during postnatal development has long-term impact on gut community structure in pigs. The ISME journal. 2008;2(7):739–48.
Pluske JR. Feed- and feed additives-related aspects of gut health and development in weanling pigs. J animal science and biotechnology. 2013;4(1):1.
Liu H, Ivarsson E, Dicksved J, Lundh T, Lindberg JE. Inclusion of chicory (Cichorium intybus L.) in pigs' diets affects the intestinal microenvironment and the gut microbiota. Appl Environ Microbiol. 2012;78(12):4102–9.
Costea PI, Hildebrand F, Arumugam M, Backhed F, Blaser MJ, Bushman FD, de Vos WM, Ehrlich SD, Fraser CM, Hattori M, et al. Enterotypes in the landscape of gut microbial community composition. Nat Microbiol. 2018;3(1):8–16.
Patra KC, Hay N. The pentose phosphate pathway and cancer. Trends Biochem Sci. 2014;39(8):347–54.
Ananthakrishnan R, Ehrlicher A. The forces behind cell movement. Int J Biol Sci. 2007;3(5):303–17.
We thank the colleagues in the State Key Laboratory of Pig Genetic Improvement and Production Technology, Jiangxi Agricultural University for their efforts in the management of experimental pig populations.
This work was supported by National Natural Science Foundation of China (31772579, 31472071). The funder was not involved in the design of the study and collection, analysis and interpretation of data and in writing the manuscript.
Availability of data and materials
The 16S rRNA gene sequencing data obtained in this study were deposited in NCBI Sequence Read Archive (SRA) under accession number SRP151335.
Ethics approval and consent to participate
All procedures involving animals were conducted according to the guidelines for the care and use of experimental animals established by the Ministry of Agricultural of China. The owner of the experimental pigs used in this study gave verbal consent for the sample collections. There were no vulnerable populations involved, and no endangered species was used in the experiments. The project was also approved by Animal Care and Use Committee (ACUC) in Jiangxi Agricultural University.
Consent for publication
The authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Summary description of dietary components for pigs at different ages. (XLSX 9 kb)
Figure S1. The venn diagram showing the sample distribution of experimental pig cohort among four ages. (TIF 269 kb)
Figure S2. PCoA showing the effect of sex on fecal bacterial composition based on the Bray-Curtis similarity matrix. (TIF 3183 kb)
Figure S3. Histogram of the linear discriminant analysis (LDA) score for differentially abundant genera between the enterotypes. (A) 25 days. (B) 80 days. (C) 120 days. (D) 240 days. Genera with LDA scores > 3 were presented. (TIF 823 kb)
Figure S4. The SparCC network at the genus level showing that the main driver of each enterotype was also the hub node. (A) 25 days. (B) 80 days. (C) 120 days. (D) 240 days. The size of the nodes represents the relative abundance of each genus. Node colors indicate the phylum that each genus belongs to. Edge colors represent positive (green) and negative (red) correlations, and the thickness of the edges indicates the values of the correlations. (TIF 1973 kb)
Table S2. Summary of annotation information of node OTUs from biologically important modules in the phylogenetic co-occurrence networks at four ages. (XLSX 15 kb)
Figure S5. Comparison of the topological features for main modules in the networks of porcine gut microbiota at different ages. (TIF 146 kb)
Table S3. Comparison of predicted functional capacities of fecal microbiome among four ages. (XLSX 15 kb)
About this article
Cite this article
Ke, S., Fang, S., He, M. et al. Age-based dynamic changes of phylogenetic composition and interaction networks of health pig gut microbiome feeding in a uniformed condition. BMC Vet Res 15, 172 (2019). https://doi.org/10.1186/s12917-019-1918-5