- Open Access
Metagenomic reconstructions of gut microbial metabolism in weanling pigs
- Weilan Wang†1, 2,
- Huifeng Hu†2,
- Ruurd T. Zijlstra1,
- Jinshui Zheng†2Email author and
- Michael G. Gänzle†1, 3Email authorView ORCID ID profile
© The Author(s). 2019
- Received: 28 November 2018
- Accepted: 11 March 2019
- Published: 26 March 2019
The piglets’ transition from milk to solid feed induces a succession of bacterial communities, enhancing the hosts’ ability to harvest energy from dietary carbohydrates. To reconstruct microbial carbohydrate metabolism in weanling pigs, this study combined 16S rRNA gene sequencing (n = 191) and shotgun metagenomics (n = 72).
Time and wheat content in feed explained most of the variation of the microbiota as assessed by 16S rRNA gene sequencing in weanling pigs. De novo metagenomic binning reconstructed 360 high-quality genomes that represented 11 prokaryotic and 1 archaeal phylum. Analysis of carbohydrate metabolism in these genomes revealed that starch fermentation is carried out by a consortium of Firmicutes expressing extracellular α-(1 → 4)-glucan branching enzyme (GH13) and Bacteroidetes expressing periplasmic neopullulanase (GH13) and α-glucosidase (GH97). Fructans were degraded by extracellular GH32 enzymes from Bacteriodetes and Lactobacillus. Lactose fermentation by β-galactosidases (GH2 and GH42) was identified in Firmicutes. In conclusion, the assembly of 360 high-quality genomes as the first metagenomic reference for swine intestinal microbiota allowed identification of key microbial contributors to degradation of starch, fructans, and lactose.
Microbial consortia that are responsible for degradation of these glycans differ substantially from the microbial consortia that degrade the same glycans in humans. Our study thus enables improvement of feeding models with higher feed efficiency and better pathogen control for weanling pigs.
- Weanling pigs
- Microbial degradation
Culture-dependent and culture-independent approaches have advanced our understanding of the assembly of intestinal microbiota and their importance for their host [1–4]. In general, the development of gut microbiota is influenced by host genetic variation [2, 5, 6], environmental factors, and stochastic events [5, 7, 8]. The association between the host genetics and gut microbiome is mediated by immunity-related pathways and the secretion of antimicrobial compounds . Multiple environmental factors, such as antibiotics, social contacts, and the environment, also shape the architecture of gut microbiota [5, 7, 9]. The diet and particularly dietary carbohydrates are key determinants for the composition and activity of the intestinal microbiome [1, 5, 7, 10, 11].
At weaning, mammals gradually transition from lactose to plant carbohydrates as main source of dietary carbohydrates; this dietary shift also induces a major shift of the intestinal microbiome. In contrast, current swine production systems impose an abrupt transition from sow’s milk to solid food in piglets; this also induces an accelerated succession of microbial communities in weanling piglets [7, 12]. Gut microbial communities play a pivotal role in facilitating adaption of weanling piglets to fibrous feed and in minimizing the risk of colonization by pathogens after weaning . Carbohydrates are the main energy source for pigs; in commercial pig production, carbohydrates account for more than 60% of the dry matter and 60–70% of the dietary energy intake [14, 15]. However, digestive enzymes secreted by the host do not degrade dietary polysaccharides other than starch. Pigs and other mammalian hosts rely on symbiotic gut microorganisms to ferment this abundance feed ingredient and provide metabolic energy . Because specific microbial taxa are specialized for degradation of specific dietary carbohydrates, the composition of the diet alters the composition and activity of intestinal microbiota .
Current analyses of intestinal microbiota are largely based on sequence-based methodology, avoiding the time-consuming culture-based analysis of intestinal microbiota. When assessing the function of intestinal microbiota on the basis of high-throughput sequencing data, metagenomic binning and genome-scale metabolic reconstructions has bridged the gap between the taxonomic analysis of microbial communities on the basis of 16S rRNA sequences and the description of the metabolic repertoire of individual members of gut microbiome by analysis of the abundance and distribution of metabolic enzymes . Metagenomic binning is also an essential tool to understand metabolic cooperativity between different representatives of the microbiome. Metabolic binning of metagenomic sequence data and the assignment of bacterial taxonomy to metabolic activity is thus an important tool to substitute untargeted microbiome modulation with targeted or predictable modulation of gut microbiome .
Swine are an important livestock species. Particularly at weaning, dietary management of the microbiome of piglets is an important tool to reduce the susceptibility to pathogens, and to improve feed efficiency. However, past studies of the interactions between gut microbiome and diets in pigs were limited to 16S rRNA gene sequencing, or metagenomic analyses without metagenomic binning [7, 18–20]. This study therefore aimed to unravel the adaptation of the swine microbiome to the dietary shift after weaning, and to establish a metagenomic reference by binning of genomes of swine gut bacteria from 72 samples from 18 animals. The metagenomic reference was used to predict the metabolic capacity of the fecal microbiome for metabolism of dietary carbohydrates by CAZy annotation, and by detailed analysis of metabolic pathways of major substrates present in wheat.
Diets, animals, and samples
Six experimental diets were prepared by mixing 98% basal diet (Additional file 1: Table S1) with 2% unfermented wheat, acidified wheat (pH 3.8), wheat fermented with Lactobacillus casei K9-1 and Lactobacillus fermentum K9-2 (CanBiocin Inc., Edmonton, AB, Canada), or unfermented wheat with freeze-dried cultures of L. casei K9–1 and L. fermentum K9-2 (approximately 109 CFU /g), wheat fermented with Lactobacillus reuteri TMW1.656, and wheat fermented L. reuteri TMW1.656ΔrtcN (Additional file 1: Table S2) . Feed fermentation was performed as previously described [12, 20].
The six dietary treatments were randomly allocated to 48 crossbred castrated male piglets (21 days of age) with randomized block design to provide 8 replicates per dietary treatment. Pigs were raised in a temperature-controlled room (28 ± 2 .5 °C) with one pig per pen and divided into six blocks. Pigs had access to ad libitum feed and clean water.
A total of 191 fecal samples were collected from the pen floors at days at weaning (0 day) and 7, 14, and 21 days after weaning. To avoid contamination of fecal samples, the pen floor was cleaned before sample collection and samples were collected during the peak hours of defecation 1–2 h after morning feeding. Fecal samples were immediately stored at − 20 °C after sampling. Subsampled samples (2–3 g) were stored at − 80 °C after thawing and mixing of frozen samples. Total bacterial DNA was extracted from fecal samples using QIAamp Fast DNA stool mini kit (Qiagen, Inc., Valencia, CA, USA) pre-treated with bead-beating (BioSpec Products, Inc., Bartlesville, USA) for 30 s × 8 times. Purified DNA with an A260/280 ratio higher than 1.8 were selected for sequencing analysis. PicoGreen quantification assay was included by default as quality control for 16S rRNA gene sequencing (University of Minnesota Genomics Center, Minneapolis, MN, USA) and shotgun metagenomic sequencing (McGill University and Genome Quebec Innovation Center, Montreal, Canada).
Intestinal microbial community analysis using 16S rRNA gene sequencing
Genomic DNAs were sequenced on Illumina MiSeq (2 × 300 bp reads) by amplifying the V5-V6 domain of the 16S rRNA gene using TruSeq Nano protocol (University of Minnesota Genomics Center, Minneapolis, MN, USA). A total of 6,647,893 sequences with an average length of 266 bp, corresponding to 34,805 16S rRNA sequences for each of the 191 samples, were retained for downstream analysis after the quality filtering from QIIME pipeline (MacQIIME 1.9.1) . Sequences with 97% similarity were clustered into operational taxonomic units (OTUs) by UCLUST  after de-replication and de-multiplexing. The GreenGenes database was used for taxonomy assignment with the default cut-off of 97% average nucleotide identity (ANI) at the genus level . OTUs that were represented by only one or two sequences (relative abundance < 0.005%) were discarded. Principle coordinates analysis (PCoA) and analysis of similarities (ANOSIM) were performed using weighted UniFrac distance matrix calculated by beta_diversity.py .
Metagenomic sequencing, assembly, binning, and genome annotation
Samples from 18 pigs were randomly selected to include 6 samples from each of the following dietary groups: acidified wheat, wheat fermented with Lactobacillus casei K9-1 and L. fermentum K9-2, or unfermented wheat with freeze-dried cultures of L. casei K9-1 and L. fermentum K9-2. Feces taken from these 18 piglets at weaning (0 day) and 7, 14, and 21 days after weaning were sequenced on Illumina Hiseq 2500 PE125 platform with low input shotgun metagenomic library protocol (McGill University and Genome Quebec Innovation Center, Montreal, Canada). A total of 399.25 Gb of sequence data were obtained, corresponding to 2.22 × 1011 reads for each of the 72 samples. After quality check by FastQC, adapters were trimmed from raw reads by Trimmomatic  using a local adapter database. Trimmed reads were assembled into contigs using IDBA_UD with default parameters.
Binning was performed with MaxBin2  using contigs longer than 3000 bp. After a two-step de-replication with dRep , 596 bins were obtained from the sample clusters pooled by pigs (four time points for each pig). CheckM  assessment indicated that all 596 bins were ≥ 50% complete; of these 596 bins, 458 bins were substantially complete (completeness ≥ 70%) and 240 bins were nearly complete (completeness ≥ 90%)  (Additional file 1: Table S3). Of the 458 substantially complete bins, 360 bins with contamination < 5% were regarded as high-quality assembled genomes and selected for further analyses.
Open reading frames (ORFs) were identified by prodigal v.2.6.1 . ORFs were annotated with BLAST against Clusters of Orthologous Groups (COG) database and CAZy database with an e value ≤ 1e-5 .
Phylogenetic identification and calculation of the relative abundance
The taxonomy of 360 high-quality bins was assigned by Phylophlan with 3737 reference genomes  on the basis of 400 proteins. The bins were assigned at the species, genus, and family level when average amino acid identity of encoded proteins to the reference genome was greater than ≥ 90%, 60%, and 45%, respectively [33, 34], in at least 50 proteins . The average coverage of bins were determined using MaxBin2  by recruiting reads (from each sample) to scaffolds. The average coverage normalized to the total number of reads in each sample corresponds to the relative abundances of bins.
Reconstruction of metabolic pathways for carbohydrate fermentation
CAZy were clustered into five categories based on the substrate specificity of glycoside hydrolases (GHs) and carbohydrate esterases (CEs). Enzymes from the GH families GH13, GH31, GH97, GH4, GH14, GH15, GH57, and GH63 were assigned to starch-degrading enzymes. GH families GH32, GH91, and GH68 were assigned as fructan-hydrolyzing enzymes. GH families containing β-glucanases including licheninase, β-glucan endohydrolase, endo-(1, 4) β-glucanase are GH8, GH16, GH26, GH5, GH6, GH9, GH10, GH12, GH44, GH48, GH45, GH51. GH and CE families harboring xylanase, arabinofuranosidase, α-glucuronsidase, and acetyl-xylan esterase were regarded as arabinoxylan-specific and include GH5, GH10, GH11, GH8, GH43, GH51, GH67, GH115, CE1, CE2, CE4, CE6, and CE7. Enzymes degrading O-linked and N-linked host glycans include GH20, GH84, GH110, GH89, GH125, GH109, CE14, GH123, and CE9. The degradation capacity of each bin corresponds to the sum of positive hits of GHs or CEs under each category.
Metabolic pathways of starch, fructan, and lactose were studied by blasting sequences of key enzymes that were characterized biochemically (Additional file 1: Table S4) against 360 assembled genomes. An amino acid identity of ≥ 40% and e value ≤ 1e-5 were used as threshold values. The relative abundance of enzymes over time was calculated by sum of all positive hits normalizing with corresponding abundance of target bins at four time points.
The data for average daily gain, feed intake of pigs, feed efficiency, and relative abundance of bins were analyzed using linear mixed-effects (LME) models in R (version 3.4.3). In the model, time was treated as fixed factor; pig was considered as experimental unit and random factor. P values < 0.05 with Bonferroni-adjustment were considered significant. Results are presented as means ± standard deviation. Alpha-diversity parameters between time points were compared using Kruskal-Wallis rank-sum test. Two-way ANOVA was applied to investigate the longitudinal differences of UniFrac distances between littermates and between all pigs from different litters. Statistical significance of ANOSIM was determined through permutations between dietary categories. The R value was calculated by the following formula: R = difference of mean rank (all distances between groups − all distances within groups) / (N(N-1)/4). The larger R value between 0 and 1 reflects the higher dissimilarity between the groups.
Growth performance and gut health of pigs
Growth performance of weanling pigs during the first 3 weeks after weanling
Average feed intakea (g DMb/day)
266.56 ± 2.90 C
468.98 ± 4.11 B
749.07 ± 16.85 A
Average daily gain (g/day)
177.31 ± 11.48 C
375.00 ± 11.33 B
615.38 ± 17.49 A
Feed efficiency (G/F)
0.65 ± 0.04 B
0.80 ± 0.02 A
0.84 ± 0.03 A
Bacterial community composition analysis by 16S rRNA gene sequencing
ANOSIM of fecal microbiota after weaning based on weighted UniFrac distance matrix calculated with partial 16S rRNA sequences
Metagenomic reconstruction of fecal microbiome in weanling pigs
The capacity of 360 bins for glycan degradation was initially predicted by identification of glycoside hydrolases (GHs) and carbohydrate esterases with similar substrate preference (Fig. 2, Additional file 1: Table S3). Starch-degrading enzymes were widely found in genomes of Firmicutes and Bacteroidetes. These metagenomic bins were relatively abundant and increased over time, particularly in genomes of Faecalibacterium spp. Only few genomes harbored fructan-degrading enzymes, examples include genomes of Lactobacillus, Escherichia coli, and several unclassified Bacteroidales. The distribution of bins carrying enzymes for β-glucan, arabinoxylan, and host-glycan metabolism overlapped; these GH families were widely distributed in genomes of Bacteroidetes, Ruminococcus, and Lachnospiraceae.
Reconstruction of metabolic pathways for starch, fructan, and lactose metabolism in weanling pigs
Firmicutes and Bacteroidetes harbored distinct pathways for starch degradation. Firmicutes convert starch by an extracellular α-(1 → 4)-glucan branching enzyme (GlgB) and pullulanases (Amy12); these enzymes occurred only in Firmicutes. GlgB was detected in 191 of 216 bins assigned to Firmicutes. Remarkably, most of the starch-degrading Firmicutes carry GlgB but no other starch-hydrolyzing enzymes; only five members of Faecalibacterium prausnitzii and one Eubacterium rectale additional carry pullulanases. Bin 254 was the only member of Firmicutes with an extracellular neopullulanase2. The periplasmic starch utilization by susA (GH13, 54.2%) and susB (GH97, 74.0%) was present in Bacteroidetes. Multiple metagenomics bins assigned to Bacteroidetes contained multiple genes for starch digestion, including extracellular and periplasmic enzymes. Among the 80 starch-degrading Bacteroidetes, only 2 Bacteroidetes were without susA or susB, whereas 23 Bacteroidetes additionally harbored extracellular amylase (Amy1, 4) or neopullulanase (SusG).
Firmicutes and Bacteroidetes also harbored distinct enzymes for fructan hydrolysis. Firmicutes catabolize fructan by intracellular β-fructofuranosidase (ScrA/ScrB, GH32) and extracellular fructansucrases (Inu, GH68, and FruA, GH32). FruA, ScrB, and Inu were present only in lactobacilli. ScrA was present in other Firmicutes, including Subdoligranulum variabile, Faecalibacterium prausnitzii, Eubacterium, and unclassified Clostridiales. Bacteroidetes metabolize fructan by β-(2 → 6) endo-fructanases (GH32) including BT_1760 (extracellular), BT_3082 (periplasmic), and BT_1765/1754 (intracellular). Among 49 fructan-degrading Bacteroidetes, only three members did not carry intracellular β-(2 → 6) endo-fructanase, nine members of Bacteroidetes additional harbored extracellular BT_1760, another two members additional carried periplasmic BT_3082.
Lactose hydrolysis was identified only in Firmicutes with exception of bin20 representing a member of Coriobacteriaceae. Most lactose-degrading bacteria (27 out of 31) hydrolyze lactose by intracellular GH2 β-galactosidase (BbgI, LacM, and BbgIV) or GH42 β-galactosidase LacA, including members of Lactobacillus, Subdoligranulum, and Ruminococcus. The abundance of Lactobacillus delbrueckii, a species that is specialized for lactose conversion, decreased over time (Fig. 2). The other lactose decomposers were capable of lactose hydrolysis by extracellular BbgIII (GH2).
Microbial composition in piglets differed by host-derived factors and reshaped by diet
The fecal microbiota of weanling pigs is unstable during early life and stabilizes within 2–3 weeks after weaning [12, 18]. In this study, we analyzed the structure and function of the fecal microbial community in piglet during the first 3 weeks after weaning by partial 16S rRNA gene sequencing and metagenomics analysis. Analysis of fecal microbiota enables sampling from the same animal over time. Swine fecal microbiota closely resemble microbiota of the colon, where complex carbohydrates are fermented, but differ substantially from ileal microbiota that are specialized on fermentation of mono- and disaccharides [43–46].
Diet together with age were the most significant factors shaping community assembly in weanling piglets; litter effects were transient and minor. Age alters the physiology of the gastrointestinal tract including immune and metabolic functions during weaning in piglets; these changes occur rapidly in response to the transition to solid food [47, 48]. Initial differences in microbiota structure between animals and litters were altered by the uniform post-weaning diet. Inclusion of probiotic bacteria modulated gastric communities of lactobacilli but did not significantly alter the structure of fecal microbiota . The use of the same probiotic L. reuteri at a tenfold higher dose significantly altered the abundance of only very few bacterial taxa . The limited impact of probiotic bacteria on the overall composition of intestinal microbiota matches observations in humans [47, 48]. Most lactobacilli used in probiotic studies are allochthonous to the hosts’ gastrointestinal tract. These allochthones are less likely to alter autochthonous microbiota than host-adapted probiotic strains [21, 48, 49].
Novel reference with 596 genomes were reconstructed for swine fecal microorganisms
The importance of the intestinal microbiome for host physiology highlights the need for comprehensive analysis based on genomic and phenotypic assays. Culture-dependent analysis of intestinal microbiota, however, lags the identification of bacterial taxa by high-throughput sequencing approaches [3, 50]. This study reconstructed 596 genomes including 360 high-quality and substantially completed genomes from swine fecal microbiota. Owing to the high error rates in sequencing and assembly of regions with repetitive sequences including rRNA operons , bacterial taxa were identified on the basis of the average amino acid identity (AAI) of multiple conserved proteins sequences distributed across the genomes . About 30% of the genomes—106 of 360 genomes—were assigned to the genus or species level. This greatly expands the taxonomic assignment of swine microbiota when compared to previous metagenomic studies conducted without metagenomic binning, which assigned 7.6% of taxa to the species or genus level , but also emphasizes the need for further metagenomic and culture-based studies to characterize previously unidentified bacterial taxa. For the remaining 254 high-quality genomes in this study, we provide novel genome sequence data and expand current knowledge on metabolic diversity of swine intestinal microorganisms [52, 53]. The phylogenetic analysis of 360 reconstructed genomes and the annotation of open reading frames serves as the first reference for metagenomes of swine microbiota. Different from past swine metagenomics studies [7, 18, 54], this reference allows the assignment of the metabolic activity of intestinal organisms to their taxonomic identification.
Phylogeny and functions of high-quality genomes reconstructed for swine fecal microorganisms
The increasing abundance of Bacteroidales and Clostridiales over time reflects their ability to derive metabolic energy from diverse plant polysaccharides [41, 55–57]. CAZy annotation indicated the ability of multiple Bacteroidetes species to degrade starch, β-glucan, arabinoxylan, and host glycans . Ruminococcaceae, a family in the phylum Firmicutes, also includes species with the capability to hydrolyze a broad range of polysaccharides, matching the identification of Ruminococcus bromii as a keystone species for starch degradation in the human colon . The relative abundance of metagenomic bins representing bacterial taxa with multiple polysaccharide-degrading enzymes increased in response to the inclusion of plant carbohydrates after weaning. Examples include members of Bacteroidetes and Ruminococcaceae. In contrast, some early colonizers, mostly Proteobacteria including E. coli and members of the genus Clostridium, decreased dramatically after weaning. The reduced abundance of Proteobacteria may relate to a lower protein intake concomitant with a higher fiber intake and a more developed immune function [47, 59].
The relative abundance of the species Lactobacillus amylovorus, Lactobacillus johnsonii, and Lactobacillus reuteri that was determined by metagenomic binning (this study) matched previous analyses of the same samples that were based on 16S rRNA sequencing and species-specific quantitative PCR . Metagenomic binning also identified Lactobacillus delbrueckii in the microbiome of piglets; this organism decreased rapidly after weaning. L. delbrueckii has not been considered a member of animal intestinal microbiota . Different from other species of the L. delbrueckii group, which maintain enzymes for metabolism of a relatively broad array of carbohydrates, the genome of L. delbrueckii underwent reductive evolution that silenced most carbohydrate metabolic enzymes . This metabolic focus of L. delbrueckii on lactose as main source of metabolic energy was interpreted as adaptation to the milk environment or dairy fermentations ; however, our data suggests that the metabolic focus on lactose may alternatively represent adaptation to the intestine of suckling mammals. Re-analysis of the intestinal microbiome of weanling piglets  indeed revealed that L. delbrueckii was also detected in piglets on the day of weaning but no longer detectable 2 or 3 weeks after weaning (Additional file 1: Figure S2).
Microbial degradation of starch, fructans, and lactose in weanling pigs
Even though a large panel of CAZymes have been cataloged based on substrate specificity, CAZy family-based classification of enzymes needs to be complemented by a more detailed analysis that is based on reference sequences of enzymes that were biochemically characterized [37, 62]. Moreover, classification of proteins in GH or CE families not always allows an unambiguous prediction of their substrate specificity or cellular location. Many substrates are degraded by enzymes from several families, and enzymes in many GH or CE families are active on more than one substrate .
In humans, starch entering the large intestine is degraded by microbial consortia contributing diverse extracellular, periplasmatic, and intracellular starch-converting and -hydrolyzing enzymes  while lactose, GOS, and dietary fructans are degraded by few bacterial groups, particularly Bifidobacterium spp. [64, 65]. Our analysis revealed that microbial consortia and species that degrade starch, fructans, and lactose in weanling piglets differ substantially from microorganisms or microbial consortia that are responsible for the corresponding metabolic activities in human intestinal microbiota.
Bacterial degradation of starch is mediated by amylases and pullulanases, which hydrolyze α-(1 → 4)- and α-(1 → 6)-glucosidic bonds, respectively . Members from GH13 families were identified as the principal starch degrading enzymes but Firmicutes and Bacteroidetes used distinct starch-utilization systems. Extracellular glycosidases were identified mainly in Firmicutes while periplasmic enzymes were only found in Bacteriodetes. The high abundance of an extracellular α-(1 → 4)-glucan branching enzyme suggests this enzyme is important for the primary degradation of starch in the swine GIT. The α-(1 → 4)-glucanotransferase GlgB catalyzes glucan chain transfer to form α-(1 → 6)-glucosidic linkages; this enzyme was found in Firmicutes only. The enzyme was suggested to improve accessibility to insoluble starch [40, 66] and is broadly distributed in intestinal microbiota of different hosts including humans, chicken, cattle, and swine . Following starch hydrolysis by extracellular enzymes, the α-glucosidases SusA and SusB further degrade gluco-oligosaccharides in the periplasm. Disruption of SusA and SusB from Bacteroides thetaiotaomicron reduced the rate of growth but did not eliminate the growth of the strain ; however, periplasmic starch-degrading enzymes may reduce access of competitors to the products of hydrolysis. With exception of F. prausnitzii, starch-degrading Firmicutes carried a single glycosidase; in contrast, redundant enzymes with different activities or different locations were commonly detected within a single genome of Bacteroidetes. The high abundance of starch-utilizing enzymes in fecal microbiota demonstrates that wheat starch, despite its hydrolysis by pancreatic amylases and brush border enzymes, is a major carbohydrate source for colonic microbiota. The distribution of extracellular and periplasmic enzymes for starch degradation highlights a high level of metabolic cooperativity that was also noted in human starch and cellulose-degrading microbial communities [37, 58]. In humans, Ruminococcus bromii plays a key role in fermentation of type 3-resistant starch and enhanced the growth of B. thetaiotaomicron, Bacteroides adolescentis, or E. rectale on resistant starch [42, 58]; however, this study indicates that this species does not fulfill a comparable role in swine microbiota.
Low molecular weight fructans are among the major non-starch polysaccharides in wheat . Fructans were degraded by Bacteriodetes and lactobacilli. The linear structure of fructans allows hydrolysis by single enzyme that is classified in the GH32 or GH68 families . In contrast to the complex and partially redundant starch-degrading enzymes in Bacteriodetes, fructans degraders carried fewer fructanases. Fructan utilization is not conserved within members of a specific species [41, 42, 70]. Three GH32 enzymes, including BT1760 (extracellular), BT3082 (periplasmic), and BT1765 (intracellular), as well as hybrid two-component (HTC) signaling system, BT1754 are required for fructan utilization in B. thetaiotaomicron and related Bacteroides spp. . B. thetaiotaomicron utilized levan while Bacteroides caccae ferments inulin [41, 69]. GH68 family enzymes in lactobacilli are extracellular levansucrases which are necessary for biofilm formation on non-secretary epithelia of the upper GI tract; these enzymes synthesize levan but do not contribute to fructan hydrolysis . Intracellular GH32 β-fructofuranosidases of lactobacilli (ScrB) utilize only di- and trisaccharides that are transported across the membrane [36, 70, 72]. The metagenomic analysis is the first to report the presence of extracellular fructanases (FruA) in intestinal lactobacilli. FruA is common in oral streptococci; however, its presence is exceptional in lactobacilli and was previously identified only in type II sourdough microbiota . The exclusive presence of FruA in Lactobacillus species representing swine intestinal communities may reflect specific nutritional requirements in pigs. The identification of Bacteroides and Lactobacillus spp. as major fructan-degrading organisms also differentiates human and swine microbiota; in humans, bifidobacteria are the main organisms that degrade fructans .
Lactose accounts for about 26.7% of sow milk solids ; transition diets contain 10–15% of lactose. Lactose is a major dietary carbohydrate in suckling and weanling pigs. Only Firmicutes fermented lactose with the LacS/LacLM pathway widely distributed in lactobacilli [75–77]. Lactose is transported into the cytoplasm by lactose permease and hydrolyzed by intracellular GH2 β-galactosidases common in Firmicutes, and GH42 β-galactosidases of lactobacilli. Lactobacilli colonize the stomach of swine where dietary lactose is available; in contrast, Bacteroidetes are dominant only in hindgut microbiota after full or partial digestion of lactose in the small intestine. Accordingly, their extracellular enzymes may target β-glycosidic linkages in host or plant glycans rather than lactose. Microbial fermentation of lactose in the terminal ileum and the large intestine contributes to lactose digestion particularly in lactase-non-persistent humans . The distribution of β-galactosidases in human microbiota remains poorly characterized; bifidobacteria are considered to be the main organisms involved in metabolism of lactose and related β-galacto-oligosaccharides .
In conclusion, we present a metagenomic reference for swine fecal microbiome by assigning taxonomies and metabolic functions to the 360 high-quality assembled genomes. Along with the clear evidence for dietary carbohydrates acting as the most significant drivers for diversification of microbiota, we further determined the key microbial contributors to degradation of major substrates in starter diet, including starch, fructans, and lactose. Starch is a substrate for colonic microbiota and its metabolism is dependent on metabolic cooperativity between Firmicutes and Bacteroidetes. Fructans and lactose are fermented by simple enzyme systems present in Bacteroides and Lactobacillus spp., respectively. Our genome-based functional analysis thus improves the understanding of carbohydrate fermentation in the swine intestine when compared to previous studies that report gene annotation without metagenomic binning [20, 44]. It will enable future studies linking composition and function of piglet microbiota to establish feeding systems that improve feed efficiency and animal health while reducing microbial resistance to antibiotics. It may also facilitate the design and interpretation of swine as a highly suitable animal model to understand carbohydrate digestion in the human intestine [51, 79–81].
Canbiocin Inc. is acknowledged for providing bacterial strains.
Ruurd Zijlstra and Michael Gänzle received funding from Alberta Agriculture and Forestry (Grant No. 2013R002R); Jinshui Zheng is supported by the National Natural Science Foundation of China (NSFC) (31770003; 31500003). Weilan Wang and Michael Gänzle acknowledge support from the China Scholarship Council and the Canada Research Chairs Program, respectively.
Availability of data and materials
The datasets during and/or analyzed during the current study are available from the corresponding author on reasonable request; sequence data will be submitted to the NCBI database.
WW and HFH conducted bioinformatics analyses, RTZ and MGG designed the animal experiment, MGG and JSZ planned the project, and WW and MGG were major contributors to the manuscript. All authors read and approved the final manuscript.
This animal study was performed at the University of Alberta Swine Research and Technology Centre (Edmonton, AB, Canada), approved by the University of Alberta Animal Care and Use Committee with guidelines of the Canadian Council on Animal Care.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis 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.
- Bäckhed F, Roswall J, Peng Y, Feng Q, Jia H, Kovatcheva-Datchary P, et al. Dynamics and stabilization of the human gut microbiome during the first year of life. Cell Host Microbe. 2015;17:690–703.View ArticleGoogle Scholar
- Blekhman R, Goodrich JK, Huang K, Sun Q, Bukowski R, Bell JT, et al. Host genetic variation impacts microbiome composition across human body sites. Genome Biol. 2015;16:191.View ArticleGoogle Scholar
- Browne HP, Forster SC, Anonye BO, Kumar N, Neville BA, Stares MD, et al. Culturing of ‘unculturable’ human microbiota reveals novel taxa and extensive sporulation. Nature. 2016;533:543–6.View ArticleGoogle Scholar
- Dethlefsen L, Eckburg PB, Bik EM, Relman DA. Assembly of the human intestinal microbiota. Trends Ecol Evol. 2006;21:517–23.View ArticleGoogle Scholar
- Benson AK, Kelly SA, Legge R, Ma F, Low SJ, Kim J, et al. Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors. PNAS. 2010;107:18933–8.View ArticleGoogle Scholar
- Camarinha-Silva A, Maushammer M, Wellmann R, Vital M, Preuss S, Bennewitz J. Host genome influence on gut microbial composition and microbial prediction of complex traits in pigs. Genetics. 2017;206:1637–44.View ArticleGoogle Scholar
- Frese SA, Parker K, Calvert CC, Mills DA. Diet shapes the gut microbiome of pigs during nursing and weaning. Microbiome. 2015;3:28.View ArticleGoogle Scholar
- Koropatkin NM, Cameron EA, Martens EC. How glycan metabolism shapes the human gut microbiota. Nat Rev Microbiol. 2012;10:323–35.View ArticleGoogle Scholar
- Chung WSF, Walker AW, Louis P, Parkhill J, Vermeiren J, Bosscher D, et al. Modulation of the human gut microbiota by dietary fibres occurs at the species level. BMC Biol. 2016;14:1–14.View ArticleGoogle Scholar
- Sharon I, Morowitz MJ, Thomas BC, Costello EK, Relman DA, Banfield JF. Time series community genomics analysis reveals rapid shifts in bacterial species, strains, and phage during infant gut colonization. Genome Res. 2013;23:111–20.View ArticleGoogle Scholar
- Sullam KE, Rubin BER, Dalton CM, Kilham SS, Flecker AS, Russell JA. Divergence across diet, time and populations rules out parallel evolution in the gut microbiomes of Trinidadian guppies. ISME J. 2015;9:1508–22.View ArticleGoogle Scholar
- Yang Y, Galle S, Le MHA, Zijlstra RT, Gänzle MG. Feed fermentation with reuteran- and levan-producing Lactobacillus reuteri reduces colonization of weanling pigs by enterotoxigenic Escherichia coli. Appl Env Microbiol. 2015;81:5743–52.View ArticleGoogle Scholar
- Molist F, van Oostrum M, Pérez JF, Mateos GG, Nyachoti CM, van der Aar PJ. Relevance of functional properties of dietary fibre in diets for weanling pigs. Anim Feed Sci Technol. 2014;189:1–10.View ArticleGoogle Scholar
- Bach Knudsen KE, Hedemann MS, Lærke HN. The role of carbohydrates in intestinal health of pigs. Anim Feed Sci Technol. 2012;173:41–53.View ArticleGoogle Scholar
- Jha R, Berrocoso JD. Review: dietary fiber utilization and its effects on physiological functions and gut health of swine. Animal. 2015;9:1441–52.View ArticleGoogle Scholar
- Magnúsdóttir S, Heinken A, Kutt L, Ravcheev DA, Bauer E, Noronha A, et al. Generation of genome-scale metabolic reconstructions for 773 members of the human gut microbiota. Nat Biotechnol. 2017;35:81–9.View ArticleGoogle Scholar
- Schmidt TSB, Raes J, Bork P. The human gut microbiome: from association to modulation. Cell. 2018;172:1198–215.View ArticleGoogle Scholar
- 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–12.View ArticleGoogle Scholar
- Yang Y, Zhao X, Le MHA, Zijlstra RT, Gänzle MG. Reutericyclin producing Lactobacillus reuteri modulates development of fecal microbiota in weanling pigs. Front Microbiol. 2015;6:762.PubMedPubMed CentralGoogle Scholar
- Xiao L, Estellé J, Kiilerich P, Ramayo-Caldas Y, Xia Z, Feng Q, et al. A reference gene catalogue of the pig gut microbiome. Nat Microbiol. 2016;1:1–6.Google Scholar
- Zhao X, Wang W, Blaine A, Kane ST, Zijlstra RT, Gänzle MG. Impact of probiotic Lactobacillus sp.on autochthonous lactobacilli in weaned piglets. J Appl Microbiol. 2018; in pressGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.View ArticleGoogle Scholar
- Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.View ArticleGoogle Scholar
- Edgar RC. Accuracy of taxonomy prediction for 16S rRNA and fungal ITS sequences. PeerJ. 2018;6:e4652.View ArticleGoogle Scholar
- Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.View ArticleGoogle Scholar
- Wu Y, Tang Y, Tringe SG, Simmons BA, Singer SW. MaxBin: an automated binning method to recover individual genomes from metagenomes using an expectation-maximization algorithm. Microbiome. 2014;2:1–18.View ArticleGoogle Scholar
- Olm MR, Brown CT, Brooks B, Banfield JF. DRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. ISME J. 2017;11:2864–8.View ArticleGoogle Scholar
- Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.View ArticleGoogle Scholar
- Svartström O, Alneberg J, Terrapon N, Lombard V, de Bruijn I, Malmsten J, et al. Ninety-nine de novo assembled genomes from the moose (Alces alces) rumen microbiome provide new insights into microbial plant biomass degradation. ISME J. 2017;11:2538–51.View ArticleGoogle Scholar
- Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:119.View ArticleGoogle Scholar
- Finn RD, Clements J, Arndt W, Miller BL, Wheeler TJ, Schreiber F, Bateman A, Eddy SR. HMMER web server: 2015 update. Nucleic Acids Res. 2015;43:W30–8.View ArticleGoogle Scholar
- Segata N, Börnigen D, Morgan XC, Huttenhower C. PhyloPhlAn is a new method for improved phylogenetic and taxonomic placement of microbes. Nat Commun. 2013;4:2304Google Scholar
- Konstantinidis KT, Rosselló-Móra R, Amann R. Uncultivated microbes in need of their own taxonomy. ISME J. 2017;11:2399–406.View ArticleGoogle Scholar
- Rodriguez-R LM, Konstantinidis KT. Bypassing cultivation to identify bacterial species culture-independent genomic approaches identify credibly distinct clusters, avoid cultivation bias, and provide true insights into microbial species. Microbe. 2014;9:111–7.Google Scholar
- Gupta RS, Lo B, Son J. Phylogenomics and comparative genomic studies robustly support division of the genus Mycobacterium into an emended genus Mycobacterium and four novel genera. Front Microbiol. 2018;9:1–41.View ArticleGoogle Scholar
- Chen C, Zhou F, Ren J, Ai L, Dong Y, Wu Z, et al. Cloning, expression and functional validation of a β-fructofuranosidase from Lactobacillus plantarum. Process Biochem. 2014;49:758–67.View ArticleGoogle Scholar
- Flint HJ, Scott KP, Duncan SH, Louis P, Forano E. Microbial degradation of complex carbohydrates in the gut. Gut Microbes. 2012;3:289–306.View ArticleGoogle Scholar
- Martens EC, Koropatkin NM, Smith TJ, Gordon JI. Complex glycan catabolism by the human gut microbiota: the Bacteroidetes Sus-like paradigm. J Biol Chem. 2009;284:24673–7.View ArticleGoogle Scholar
- Ryan SM, Fitzgerald GF, Van Sinderen D. Screening for and identification of starch-, amylopectin-, and pullulan-degrading activities in bifidobacterial strains. Appl Env Microbiol. 2006;72:5289–96.View ArticleGoogle Scholar
- Scott KP, Martin JC, Chassard C, Clerget M, Potrykus J, Campbell G, et al. Substrate-driven gene expression in Roseburia inulinivorans: importance of inducible enzymes in the utilization of inulin and starch. PNAS. 2011;108:4672–9.View ArticleGoogle Scholar
- Sonnenburg ED, Zheng H, Joglekar P, Higginbottom SK, Firbank SJ, Bolam DN, et al. Specificity of polysaccharide use in intestinal bacteroides species determines diet-induced microbiota alterations. Cell. 2010;141:1241–52.View ArticleGoogle Scholar
- Ze X, Ben DY, Laverde-Gomez JA, Dassa B, Sheridan PO, Duncan SH, et al. Unique organization of extracellular amylases into amylosomes in the resistant starch-utilizing human colonic firmicutes bacterium Ruminococcus bromii. MBio. 2015;6:1–11.View ArticleGoogle Scholar
- Loponen J, Gänzle M, Loponen J, Gänzle MG. Use of sourdough in low FODMAP baking. Foods. 2018;7:96.View ArticleGoogle Scholar
- Looft T, Allen HK, Cantarel BL, Levine UY, Bayles DO, Alt DP, et al. Bacteria, phages and pigs: the effects of in-feed antibiotics on the microbiome at different gut locations. ISME J. 2014;8:1566–76.View ArticleGoogle Scholar
- Metzler-Zebeli BU, Hooda S, Pieper R, Zijlstra RT, Van Kessel AG, Mosenthin R, et al. Nonstarch polysaccharides modulate bacterial microbiota, pathways for butyrate production, and abundance of pathogenic escherichia coli in the pig gastrointestinal tract. Appl Environ Microbiol. 2010;76:3692–701.View ArticleGoogle Scholar
- Zhang Q, Widmer G, Tzipori S. A pig model of the human gastrointestinal tract. Gut Microbes. 2013;4:193–200.View ArticleGoogle Scholar
- Pluske JR. Feed- and feed additives-related aspects of gut health and development in weanling pigs. J Anim Sci Biotechnol. 2013;4:1.View ArticleGoogle Scholar
- Walter J, Maldonado-Gomez MX, Martinez I. To engraft or not to engraft: an ecological framework for gut microbiome modulation with live microbes. Curr Opin Biotechnol. 2018;49:129–39.View ArticleGoogle Scholar
- Kristensen NB, Bryrup T, Allin KH, Nielsen T, Hansen TH, Pedersen O. Alterations in fecal microbiota composition by probiotic supplementation in healthy adults: a systematic review of randomized controlled trials. Genome Med. 2016;8:52.View ArticleGoogle Scholar
- Hugerth LW, Larsson J, Alneberg J, Lindh MV, Legrand C, Pinhassi J, et al. Metagenome-assembled genomes uncover a global brackish microbiome. Genome Biol. 2015;16:1–18.View ArticleGoogle Scholar
- Peng Y, Leung HCM, Yiu SM, Chin FYL. IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics. 2012;28:1420–8.View ArticleGoogle Scholar
- Hug LA, Baker BJ, Anantharaman K, Brown CT, Probst AJ, Castelle CJ, et al. A new view of the tree of life. Nat Microbiol. 2016;1:1–6.View ArticleGoogle Scholar
- Mukherjee S, Seshadri R, Varghese NJ, Eloe-Fadrosh EA, Meier-Kolthoff JP, Göker M, et al. 1,003 reference genomes of bacterial and archaeal isolates expand coverage of the tree of life. Nat Biotechnol. 2017;35:676–83.View ArticleGoogle Scholar
- Di Cesare A, Fontaneto D, Doppelbauer J, Corno G. Fitness and recovery of bacterial communities and antibiotic resistance genes in urban wastewaters exposed to classical disinfection treatments. Env Sci Technol. 2016;50:10153–61.View ArticleGoogle Scholar
- Hamaker BR, Tuncil YE. A perspective on the complexity of dietary fiber structures and their potential effect on the gut microbiota. J Mol Biol. 2014;426:3838–50.View ArticleGoogle Scholar
- Metzler-Zebeli BU, Schmitz-Esser S, Mann E, Grüll D, Molnar T, Zebeli Q. Adaptation of the cecal bacterial microbiome of growing pigs in response to resistant starch type 4. Appl Environ Microbiol. 2015;81:8489–99.View ArticleGoogle Scholar
- Robert C, Bernalier-Donadille A. The cellulolytic microflora of the human colon: evidence of microcrystalline cellulose-degrading bacteria in methane-excreting subjects. FEMS Microbiol Ecol. 2003;46:81–9.View ArticleGoogle Scholar
- Ze X, Duncan SH, Louis P, Flint HJ. Ruminococcus bromii is a keystone species for the degradation of resistant starch in the human colon. ISME J. 2012;6:1535–43.View ArticleGoogle Scholar
- Haenen D, Zhang J, Souza da Silva C, Bosch G, van der Meer IM, van Arkel J, et al. A diet high in resistant starch modulates microbiota composition, SCFA concentrations, and gene expression in pig intestine. J Nutr. 2013;143:274–83.View ArticleGoogle Scholar
- Duar RM, Lin XB, Zheng J, Martino ME, Grenier T, Perez-Munoz ME, et al. Lifestyles in transition: evolution and natural history of the genus Lactobacillus. FEMS Microbiol Rev. 2017;41:S27–48.View ArticleGoogle Scholar
- van de Guchte M, Penaud S, Grimaldi C, Barbe V, Bryson K, Nicolas P, et al. The complete genome sequence of Lactobacillus bulgaricus reveals extensive and ongoing reductive evolution. PNAS. 2006;103:9274–9.View ArticleGoogle Scholar
- Sørensen HR, Pedersen S, Jørgensen CT, Meyer AS. Enzymatic hydrolysis of wheat arabinoxylan by a recombinant “minimal” enzyme cocktail containing β-xylosidase and novel endo-1,4-β-xylanase and α-L-arabinofuranosidase activities. Biotechnol Prog. 2007;23:100–7.View ArticleGoogle Scholar
- Cockburn DW, Orlovsky NI, Foley MH, Kwiatkowski KJ, Bahr CM, Maynard M, et al. Molecular details of a starch utilization pathway in the human gut symbiont Eubacterium rectale. Mol Microbiol. 2015;95:209–30.View ArticleGoogle Scholar
- Reichardt N, Vollmer M, Holtrop G, Farquharson FM, Wefers D, Bunzel M, et al. Specific substrate-driven changes in human faecal microbiota composition contrast with functional redundancy in short-chain fatty acid production. ISME J. 2018;12:610–22.View ArticleGoogle Scholar
- Krumbeck JA, Rasmussen HE, Hutkins RW, Clarke J, Shawron K, Keshavarzian A, et al. Probiotic Bifidobacterium strains and galactooligosaccharides improve intestinal barrier function in obese adults but show no synergism when used together as synbiotics. Microbiome. 2018;6:121.View ArticleGoogle Scholar
- Rumbak E, Rawlings DE, Lindsey GG, Woods DR. Characterization of the Butyrivibrio fibrisolvens glgB gene, which encodes a glycogen-branching enzyme with starch-clearing activity. J Bacteriol. 1991;173:6732–41.View ArticleGoogle Scholar
- Lee S, Cantarel B, Henrissat B, Gevers D, Birren BW, Huttenhower C, et al. Gene-targeted metagenomic analysis of glucan-branching enzyme gene profiles among human and animal fecal microbiota. ISME J. 2014;8:493–503.View ArticleGoogle Scholar
- D’Elia JN, Salyers AA. Contribution of a neopullulanase, a pullulanase, and an alpha-glucosidase to growth of Bacteroides thetaiotaomicron on starch. J Bacteriol. 1996;178:7173–9.View ArticleGoogle Scholar
- Wilson B, Whelan K. Prebiotic inulin-type fructans and galacto-oligosaccharides: definition, specificity, function, and application in gastrointestinal disorders. J Gastroenterol Hepatol. 2017;32:64–8.View ArticleGoogle Scholar
- Petrov K, Popova L, Petrova P. High lactic acid and fructose production via Mn2+−mediated conversion of inulin by Lactobacillus paracasei. Appl Microbiol Biotechnol. 2017;101:4433–45.View ArticleGoogle Scholar
- Schwab C, Vogel R, Gänzle MG. Influence of oligosaccharides on the viability and membrane properties of lactobacillus reuteri TMW1.106 during freeze-drying. Cryobiology. 2007;55:108–14.View ArticleGoogle Scholar
- Pridmore RD, Berger B, Desiere F, Vilanova D, Barretto C, Pittet A-C, et al. The genome sequence of the probiotic intestinal bacterium Lactobacillus johnsonii NCC 533. Proc Natl Acad Sci U S A. 2004;101:2512–7.View ArticleGoogle Scholar
- Flint HJ, Duncan SH, Louis P. The impact of nutrition on intestinal bacterial communities. Curr Opin Microbiol. 2017;38:59–65.View ArticleGoogle Scholar
- Klobasa F, Werhahn E, Butler JE. Composition of sow milk during lactation. J Anim Sci. 1987;64:1458–66.View ArticleGoogle Scholar
- Gänzle MG, Follador R. Metabolism of oligosaccharides and starch in Lactobacilli: a review. Front Microbiol. 2012;3:340.View ArticleGoogle Scholar
- Honda H, Kataoka F, Nagaoka S, Kawai Y, Kitazawa H, Itoh H, et al. β-Galactosidase, phospho-β-galactosidase and phospho-β- glucosidase activities in lactobacilli strains isolated from human faeces. Lett Appl Microbiol. 2007;45:461–6.View ArticleGoogle Scholar
- Schwab C, Sørensen KI, Gänzle MG. Heterologous expression of glycoside hydrolase family 2 and 42 β-galactosidases of lactic acid bacteria in Lactococcus lactis. Syst Appl Microbiol. 2010;33:300–7.View ArticleGoogle Scholar
- Arola H, Tamm A. Metabolism of lactose in the human body. Scand J Gastroenterol. 1994;29:21–5.View ArticleGoogle Scholar
- Krumbeck JA, Rasmussen HE, Hutkins RW, Clarke J, Shawron K, Keshavarzian A WJ. Probiotic Bifidobacterium strains and galactooligosaccharides improve intestinal barrier function in obese adults but show no synergism when used together as synbiotics. Microbiome. 2018;6(1):121.Google Scholar
- Guilloteau P, Zabielski R, Hammon HM, Metges CC. Nutritional programming of gastrointestinal tract development. Is the pig a good model for man? Nutr Res Rev. 2010;23:4–22.View ArticleGoogle Scholar
- Heinritz SN, Mosenthin R, Weiss E. Use of pigs as a potential model for research into dietary modulation of the human gut microbiota. Nutr Res Rev. 2013;26:191–209.View ArticleGoogle Scholar