Potential of Pectins to Beneficially Modulate the Gut Microbiota Depends on Their Structural Properties (2024)

Pectins

Pectins (nine in total) were produced by CP Kelco (Denmark) from orange (P1 and P8), lemon (P2, P3, P6, and P9), lime (P5 and P7) and sugar beet (P4) using different extraction methods (harsh, mild, and differentially extracted) as shown in Table 1. Pectins P5 and P6, respectively, were chemically and enzymatically deesterified, and pectins P7 and P8 were amidated. Pectic derivative RGI (P10) was purified from lime.

High Performance Anion Exchange Chromatography (HPAEC)

Monosaccharide composition of pectins (P1–P10) was determined by HPAEC equipped with a PA20 column and Pulsed Amperometric Detection (PAD, Dionex, CA, United States) as reported previously (Larsen et al., 2018). Monosaccharides, GalA, Ara, Rha, Gal, Glc, and Xyl were quantified from two independent HPAEC-PAD analyses. DBr of the rhamnosyl residues with neutral sugars side chains was calculated by the equation DBr(%) = 100% × Rha(%)/[(Ara(%) + Gal(%)]. The molar content of hom*ogalacturonan (HG) and rhamnogalacturonan (RG) was estimated using equations: HG(%) = GalA(%) – Rha(%); RG(%) = 2Rha(%) + Ara(%) + Gal(%) (Yapo and Koffi, 2014). Differences between the pectins in the composition of monosaccharides were tested by the one-way ANOVA, Tukey’s post hoc test (Supplementary Table S1).

TIM-2 Fermentations and Sample Collection

In vitro fermentations were performed using TIM-2 model of proximal colon model, providing fully anaerobic conditions, peristaltic movements and removal of metabolites during fermentation. The TIM-2 model has been validated and used as a predictive model for clinical trials (Venema, 2015). The microbiota inoculum consisted of an active, pooled fecal samples from 8 healthy Caucasian adults (male: n = 4, female: n = 4; age of 25–42 years). Fecal samples were maintained under anaerobiosis by using anaerobic packs (AnaeroGen^TM, Oxoid, Cambridge, United Kingdom). The fecal microbiota was hom*ogenized under anaerobic conditions, snap-frozen in liquid nitrogen and stored at -80°C before inoculation in TIM-2 (Aguirre et al., 2014). The TIM-2 colon model has been described previously, including operational units, maintenance, conditions of fermentation, and sample loading (Aguirre et al., 2014). Briefly, all units of the system were flushed with nitrogen prior to inoculation and throughout the experiments. Fermentations were performed at 37°C with the pH kept at 5.8 (or slightly above) by automatic titration with 2M NaOH. A 30 ml portion of culture hom*ogenate was used to inoculate the units for each experiment. Following inoculation the microbiota was incubated for 16 h for adaptation to the new environment. During this period the basal SIEM, composed of complex carbohydrates, protein, ox-bile, Tween 80, vitamins, and minerals (Venema, 2015), was gradually introduced into the system in a total volume of 40 ml. After the adaptation, the culture was deprived from any medium for 2 h (starvation period). Afterward, the units were fed with SIEM containing 7.5 g pectin per day as the only carbohydrate source. The doses of pectin in TIM-2 fermentations was based on the amount of carbohydrates present in SIEM, and corresponded to approximately 250 g citrus fruit (grapefruit, lemon, and orange) with the typical content of pectins of 2.3–4.5% (Thakur et al., 1997). Luminal and dialysate samples were removed at the start of fermentation (0 h) and after 24, 48, 56, and 72 h fermentation for analyses of microbiota, SCFA and BCFA. Samples were snap-frozen in liquid nitrogen and stored (-80°C) until analysis. Fermentation of each substrate was performed in two independent TIM-2 experiments.

DNA Purification

Microbiota samples from TIM-2 experiments were centrifuged (10,000 ×g, 10 min) and the fecal water discarded. The total bacterial DNA was isolated from the pellet using the PowerLyzer@PowerSoil DNA Isolation Kit (Qiagen Nordic, Denmark) according to manufacturer’s protocol with few modifications. Modifications included resuspension of fecal slurries in the Bead Solution, transfer to the PowerLyser^® Glass bead tubes and heating the tubes at 65°C for 10 min before the bead beating step. The hom*ogenization was performed at speed setting of 6.5 m/s for 3 cycles, 30 s each with 30 sec interval (FastPrep-24^®, MP Biomedicals, Solon OH, United States). The concentration and purity of DNA was determined using NanoDrop 2000 Spectrophotometer (Thermofisher Scientific, Denmark). The DNA concentration was typically in the range of 20–100 ng/μl and A260/A280 ratios of 1.8–2.1.

16S rRNA Gene Amplicon Sequencing and Data Processing

The DNA library for amplicon sequencing was prepared according to Williams et al. (2017). In brief, the V3 region (∼190 bp) of the 16S rRNA gene was PCR amplified using AccuPrime SuperMix II (Life Technologies, CA, United States) and primers NXt_388F and NXt_518R compatible with Nextera Index Kit (Illumina). Amplification steps included initial denaturation at 95°C for 2 min, 33 cycles of denaturation at 95°C for 15 s, annealing of primer at 55°C for 15 s and elongation at 68°C for 30 s. The second PCR (PCR II) was performed using Phusion High-Fidelity PCR Master Mix (ThermoFisher Scientific, MA, United States) and primers P5 and P7 (Nextera Index Kit) to incorporate adapters and tags in the PCR product. The PCRII setup included 13 cycles of denaturation at 98°C for 10 s, annealing of primer at 55°C for 20 s and elongation at 72°C for 20 s, followed by an extension at 72°C for 5 min. Amplifications were performed in a Veriti 96-well thermal cycler (Applied Biosystems, United States). The amplified fragments with adapters and tags were purified using AMPure XP beats (Beckman Coulter Genomic, CA, United States). Tag-encoded 16S rRNA gene sequencing was performed using Illumina NextSeq 550 Sequencing System using the 2 × 150 cycles MID V2 kit (Illumina, CA, United States). The raw dataset, containing pair-ended reads with corresponding quality scores, were trimmed, merged, clustered [operational taxonomic units (OTU) with 97% similarity], filtered from chimeric sequences using UPARSE and taxonomically assigned using the GreenGenes database (version 12.10) (Williams et al., 2017).

Analysis of Sequencing Data and Statistic

Quantitative insight into microbial ecology (QIIME) software package (version 1.8.0) was used for subsequent analysis steps (Caporaso et al., 2010). The sequences were classified using Greengenes database (version 13.8) as a reference 16S rRNA gene database (Mcdonald et al., 2012). Alpha-diversity was evaluated by the number of observed species and Chao1 (97% similarity OTUs) computed for each OTU table rarefied to 10,000 sequences per sample, based on the lowest number of sequences produced per sample. Good’s coverage estimates were over 98.2% across the samples (data not shown), suggesting that sampling depth of 10,000 reads was sufficient to capture majority of the OTUs. Differences in alpha-diversity over time were tested using Student’s t-test, employing the non-parametric (Monte Carlo) method (999 permutations). Differences between the pectins in Chao1 indices were assessed by a one-way ANOVA with LSD post hoc test using data obtained in two independent experiments. Community differences (beta-diversity) were revealed by weighted and unweighted UniFrac distances matrices computed from the rarefied OTU tables and evaluated by analysis of similarities (ANOSIM). PCoA plots were generated with Jackknifed Beta Diversity workflow using 10 subsampled OTU tables. Group differences in microbial community structure were assessed by hierarchical clustering of OTUs (97% similarity), applying UPGMA (unweighted pair group method using arithmetic averages) algorithm and visualized using R-studio software heatmap2 (cut-off 0.01%). Statistical differences in bacterial abundances were evaluated for pooled data set (48, 56, and 72 h from two independent experiments, n = 6 for each substrate) compared to time 0 h (baseline, n = 20) by the non-parametric Wilcoxon Rank Sum test combined with Bonferroni multiplicity correction, using 0.05 as significance level. Correlation between the relative abundances of OTUs and the structural characteristics of nine pectins (P1–P9) was evaluated with the Spearman rank correlation test implemented in the otu_category_significance.py script (QIIME 1.8.0). Statistical significance was evaluated from the conservative FDR corrected p-values for multiple comparisons (FDR-adjusted, p < 0.05).

Analysis of Short Chain Fatty and Branched Chain Fatty Acids

The amounts of SCFA and BCFA were determined in the lumen and dialysate TIM-2 samples by Brightlabs (Venlo, Netherlands). Samples were centrifuged (10,000 ×g, 20 min) and the supernatants were diluted with 1.5 mM sulfuric acid. The samples (10 μl) were injected onto an 883 Basic IC plus ion-chromatography system with suppressed conductivity detection used in positive mode (Metrohm, Herisa, Switzerland). The acids were separated by isocratic elution on ICsep ION300 Ion exclusion column (300 mm × 7.8 mm, 7 μm particle size) and Metrosep RP2 Guard column (Transgenomic, New Haven, CT, United States) using 1.5 mM sulfuric acid as a mobile phase. Samples were eluted at a flow rate of 0.4 ml/min with a pressure of 5.5 MPa at 65°C. Acetic, propionic, butyric, iso-butyric, and iso-valeric acids were used as calibration standards (Sigma Chemical, St. Louis, MO, United States). Statistical differences between the pectins in SCFA and BCFA production over time were assessed by one-way ANOVA, Tukey’s post hoc test. Differences between the high methoxyl pectins (P1, P2, P3, and P4) and LM pectins (P5, P6, P7, P8, and P9) in fatty acids production were evaluated by the two-sided Wilcoxon Rank Sum test. Correlation between the relative abundances of OTUs and cumulative production of fatty acids was assessed with the Spearman rank correlation test (QIIME 1.8.0) and statistical significance was evaluated from the FDR corrected p-values. In all tests differences were assumed to be statistically significant at p-values < 0.05.

Composition and Structural Properties of Pectins

In total nine pectins (P1–P9) and RGI (P10), were used as substrates for microbiota fermentation in TIM-2 colon model (Table 1). The structural and physicochemical properties of pectins have been described previously (Larsen et al., 2018). The composition of monosaccharides referred to in this study is presented in Supplementary Table S1. Pectins P1, P2, P3, and P4 were high methoxyl pectins with highest DE of 70.0 and 74.7% in P2 and P3, respectively; others were LM pectins with DE between 11.4 and 35.6%, lowest in P5. The major neutral sugars in pectin molecules were galactose (Gal, 9.1–30.7%), arabinose (Ara, 0.8–16.7%), rhamnose (Rha, 1.7–3.7%) and Glc (Glc, 0.7–8.7%), while Xyl was found in lesser amounts (<1%). Content of GalA and HG was largest in pectins P2, P3, P6, and P9 from lemon (GalA of 70.7–83.0% and HG of 69.6–81.2%) and lowest in sugar beet pectin (P4) (46.9% of GalA and 43.9% of HG). Additionally, sugar beet pectin was characterized by relatively higher fractions of RG (50.7%) and neutral sugars (Gal, Ara, and Rha). RGI (P10) differed from pectins by low percentage of GalA (11.0%) and Gal (7.9%), and high fraction of Ara (69.9%), Rha (6.6%), and Xyl (1.8%). Estimated DBr of rhamnosyl residues with neutral sugars side chains varied between 7.4% (P3) and 19.2% (P6).

Alpha-Diversity of Microbiota

Changes in alpha-diversity of microbiota over time are presented by Chao1 indices (Figure 1A,B) and observed species number (Supplementary Figure S1). Both Chao1 and the number of observed species increased significantly after 24 h fermentation (Chao1of 746 ± 30) compared to baseline (Chao1of 723 ± 22) and, afterward, decreased gradually reaching the lowest values at the end of fermentation (Chao1of 640 ± 27 after 72 h) (Figure 1A). Fermentations of pectins P7 and P8 generally produced more diverse microbiota communities after 24–48 h compared to other substrates (Chao 1 of 724 ± 30 and 765 ± 12 after 48 h, respectively), whereas diversity in P6 and P9 fermentations was lowest (Chao 1 of 653 ± 36 and 651 ± 22 after 48 h, respectively) (Figure 1B). Differences between the substrates in alpha-diversity indices after 72 h fermentation were insignificant.

Open in a separate window

FIGURE 1

Analysis of alpha-diversity over fermentation time and between the pectins presented by Chao1 index. Microbiota was analyzed using Illumina NextSeq 550 Sequencing System. (A) Chao1 at baseline (0 h) and after 24, 48, and 72 h fermentation computed for combined samples (n = 20 for each time point). Asterisks denote the values significantly different from time 0 h (non-parametric Student’s t-test, ^∗p < 0.05, ^∗∗p < 0.01). (B) Chao1 after 24 h (brown columns), 48 h (green columns), and 72 h (purple columns) fermentation of pectins P1–P10, presented by means and standard deviation (bars) from two biological replicates. Superscripts (a, b, c, and d) show significant differences between the substrates at time-points 24 and 48 h (one-way ANOVA, LSD post hoc t test, p < 0.05). Chao1 indices at 72 h were not significantly different (letters not shown).

Group Diversity of Microbiota

Analysis of group diversity (beta-diversity) across the sampling times (Figure 2A,B) and the pectins (Figure 3A,B) is presented by the PCoA plots using unweighted and weighted UniFrac distance matrices. Unweighted UniFrac matrix is based on the presence/absence of bacterial species, accounting for both abundant and rare lineages; while weighted matrix relies on the absolute proportions, and it is most sensitive to detect the changes in dominant taxa. Major shifts in microbial diversity occurred after 24 h fermentation, as seen from both unweighted and weighted PCoA and confirmed by ANOSIM (p = 0.001) (Figure 2A,B). Differences in group diversity between the time-points 48, 56, and 72 h were insignificant, especially based on weighted UniFrac (R = 0.002, p = 0.774; Figure 2B), indicating stabilization of microbial communities after 48 h fermentation. Variation between the pectins could be best explained by the weighted PCoA (PC1 66% and PC2 27%), suggesting the major effect of abundant bacterial taxa for group differentiation (Figure 3B). Pectins P4 and P10 (group 1) were clearly separated from other substrates along the PC1 (weighted PCoA, Figure 3B). Other distinct groups combined pectins P1, P2, and P3 (group 2), and pectins P7 and P8 (group 3) (Figure 3A,B).

Open in a separate window

FIGURE 2

Principal coordinate analysis (PCoA) plots, showing group diversity over time, 0–72 h (middle plots), 0–24 h (plots to the left), and 48–72 h (plots to the right), assessed by unweighted (A) and weighted (B) UniFrac distance matrices. Samples are colored according to the time-points: blue dots – 0 h, red dots – 24 h, green dots – 48 h, yellow dots – 56 h, and pink dots – 72 h. Samples were collected from two independent fermentations of each substrate (P1–P10) and analyzed for microbiota composition using Illumina NextSeq 550 Sequencing System. Group differences were tested by ANOSIM as shown by R-values and p-values. The samples are plotted on the first two principal coordinates PC1 and PC2. The ellipse center indicates group means for each time-point.

Open in a separate window

FIGURE 3

Principal coordinate analysis, showing group diversity between the pectins (P1–P10), as determined by unweighted (A) and weighted (B) UniFrac distance matrices. PCoA plots were constructed for the samples collected after 48, 56, and 72 h fermentation from two independent experiments. Microbiota composition was analyzed using Illumina NextSeq 550 Sequencing System. The ellipse center indicates group means. The samples are plotted on the first two principal coordinates PC1 and PC2.

Effect of Pectins on Microbiota Composition

In total 1218 OTUs were assigned to more than 150 bacterial taxa at species-level. Significant changes in bacterial abundances after 48 h fermentation of pectins (P1–P10) as shown by the heatmap in Figure 4. Regarding microbiota similarity between 48 and 72 h fermentation (Figure 2A,B), statistical differences between the pectins were evaluated for the combined data. The hierarchical clustering of individual samples and the relative abundances of species-level OTUs (cut-off 0.01%) at baseline (0 h) and after 48, 56, and 72 h fermentation are presented in Supplementary Figure S2. Four bacterial taxa were prevalent at fermentation start (64.5% in total), including family Lachnospiraceae (19.3%), family Ruminococcaceae (16.5%) and genus Ruminococcus (10.4%) within phylum Firmicutes, and species Prevotella copri (18.3%) within phylum Bacteroidetes. Other 26 common bacterial populations, belonging to phyla Firmicutes and Bacteroidetes, were found in moderate numbers (0.1–10%); among them, genera Blautia, Dorea, Coprococcus, Ruminococcus, Bacteroides, Prevotella and species F. prausnitzii. Less abundant phyla at baseline were Actinobacteria (0.41%), presented mainly by genera Bifidobacterium and Collinsella, phylum Proteobacteria (0.09%), presented by Desulfovibrio spp., Sutterella and Enterobacteriaceae, and phylum Tenericutes (0.02%), comprising Anaeroplasma and order RF39.

Open in a separate window

FIGURE 4

Changes in relative abundances of species level OTUs in TIM-2 fermentations of pectins (P1–P10). Microbiota composition was analyzed using Illumina NextSeq 550 Sequencing System. The heatmap and dendograms were generated using UPGMA algorithm and R-studio software (heatmap2). Abundances at fermentation start (0 h) are shown by numbers (means from two independent experiments, n = 20). Statistical differences after 48 h fermentation were evaluated for combined data (48, 56, and 72 h) obtained for each pectin in two independent experiments (n = 6), using the Wilcoxon Rank Sum test and Bonferroni multiplicity correction (p < 0.05). Significant increases in relative abundances compared to baseline (0 h) are marked by dots (∙) and decreases by asterisks (^•). Taxa denoted as “Other” indicate ambiguity in the assignment; taxa in square brackets indicate a proposed taxonomy. The color key presents the log-transformed mean values.

Cluster analysis revealed grouping of pectins into several clusters after similarity in microbiota profiles (Figure 4 and Supplementary Figure S2). A distinct cluster combined pectins P4 and P10, which separation was likely driven by comparatively higher numbers of P. copri, Bacteroides plebeius, Ruminococcus gnavus and lower proportions of Ruminococcaceae (Figure 4). Pectins P7 and P8 were merged into another cluster, having increased levels of order Clostridiales, Coprococcus, Bacteroides, family Christensenellaceae and Proteobacteria (Disulfovibrio D168, Klebsiella, and Enterobacter), along with reduction in F. prausnitzii. Grouping of pectins P1, P2, P3, and P9 was characterized by lower abundances of Coprococcus and Paraprevotellaceae, and relatively high numbers of F. prausnitzii. Clustering of P5, P6 can be partially explained by higher proportions of Lachnospira, Butyrivibrio, Phascolarctobacterium and Coprococcus, along with lower fraction of P. copri. Bacterial populations generally reduced by all pectins, included Ruminococcus spp., genera Blautia, Roseburia, Catenibacterium, Bifidobacterium, Paraprevotella, and families Lachnospiraceae and Erysipelotrichaceae. Concurrently, abundances of Lachnospira, Oscillospira, Clostridium, Butyrivibrio, Prevotella, family Ruminococcaceae and order Bacteroidales, were mostly increased. Additionally, the rare taxa, i.e., genus Mogibacterium, family Enterobacteriaceae (genus Citrobacter) and order RF39 (phylum Tenericutes) were stimulated in all fermentations, though to a different extent (Figure 4).

Correlation Between Bacterial Abundances and the Structural Properties of Pectins

Figure 5 shows correlation between the relative abundances of species-level OTUs and characteristics of pectins, i.e., content of monosaccharides, GalA, RG fraction, DE and DBr. The levels of Oscillospira, Blautia, Dorea, Ruminococcus, Coprococcus, R. torques, Lachnospiraceae and Clostridiales within phylum Firmicutes, and Paraprevotella, B. uniformis, B. ovatus, P. distasonis and Prevotella within phylum Bacteroidetes, correlated significantly (r = 0.35–0.71, p < 0.05) or showed a tendency to increase with higher fractions of simple sugars (Gal, Rha, Xyl, and Glc) and lower content of GalA and DE values. Species F. prausnitzii and family Ruminococcaceae showed positive correlation with the sugar residues, GalA and DE (r = 0.42–0.72, p < 0.01), in contrast to bacterial taxa above. Abundance of P. copri correlated positively with Ara (r = 0.61, p < 0.001) and DE (r = 0.77, p < 0.001) and negatively with DBr (r = 0.49, p < 0.01). Additionally, genera Coprococcus and Lachnospira correlated positively with DBr (r = 0.48–0.68, p < 0.01) and negatively (r = 0.48–0.82, p < 0.01) with DE and Ara.

Open in a separate window

FIGURE 5

Correlation between the relative abundances of bacterial taxa and the structural characteristics of pectins (P1–P9), i.e., content of galactose (Gal), rhamnose (Rha), xylose (Xyl), glucose (Glc), arabinose (Ara), galacturonic acid (GalA), rhamnogalacturonan (RG), degree of branching (DBr), and degree of esterification (DE). Analysis was performed by the Spearman rank correlation test using combined data (48, 56, and 72 h), obtained for each pectin from two independent experiments. Significant correlation is indicated by the correlation coefficient (FDR-adjusted, p < 0.05).

Production of Short Chain and Branched Chain Fatty Acids

The cumulative production of SCFA (acetate, propionate and n-butyrate) after 72 h fermentation is shown in Table 2. The amounts of acetate and total SCFA were not different between the pectins, comprising 81–95 mmol and 131–162 mmol, respectively. Production of propionate was largest for sugar beet pectin (P4) (43.3 mmol) and RGI (P10) (47.6 mmol), and differed significantly from the citrus pectins P2, P5, P6, P7, and P9 (24.2–30.1 mmol). Generally, the HM pectins (P1, P2, P3, and P4) generated higher amounts (p < 0.05) of propionate and total SCFA compared to LM pectins (P5, P6, P7, P8, and P9). Production of butyrate was lowest in fermentations of P10 (19.3 mmol) and statistically different from pectins P1, P2, P3, and P9 (27.9–32.9 mmol). The lowest acetate to propionate ratios (2.0–2.4) were determined for P4, P1, and P10 and highest (3.2–3.6) for P5 and P6 fermentations. Production of total BCFA, i-butyrate and i-valerate were slightly lower in RGI fermentations, though not significantly different from other treatments (Supplementary Table S2). The amount of propionate throughout fermentation (24–72 h) correlated positively with P. copri, R. torques and unclassified Clostridiales, and, negatively, with genera Coprococcus, Butyrivibrio and Lachnospira (Table 3). Production of butyrate decreased with increasing abundances of B. ovatus and Bacteroides after 48 h fermentation. Correlation between other bacterial taxa and SCFA was insignificant.

Increase in Bacterial Abundances Was Related to Species-Specific Utilization of Pectins

In this study we conducted TIM-2 fermentations with diverse pectins to investigate the relationship between the structural properties of pectins and the changes in microbiota composition. Enrichment or reduction of bacterial populations during fermentation would depend on their ability to degrade pectins, and/or utilize POS and other metabolites (cross-feeding interactions). Among bacterial taxa increased in this study, the members of Bacteroidales, genera Prevotella, F. prausnitzii and family Enterobacteriaceae, have been reported as fiber- or pectin-degraders. In addition to the primary pectin-degraders Bacteroidetes and Prevotella, described in introduction, structure and activity of a large variety of pectinases have been characterized within the members of Enterobacteriaceae (Erwinia, Yersinia, and E. coli) and Clostridium species of intestinal origin (Nakajima et al., 2002; Abbott and Boraston, 2008). Ability of intestinal isolates of Clostridium, Bacteroides, Enterococcus, and F. prausnitzii to utilize pectins was previously confirmed in culture growth experiments (Olano-Martin et al., 2002; Lopez-Siles et al., 2012). Besides, Oscillospira and Clostridium species were found to be associated with degradation of plant material and enhanced by plant fiber content in human diet (Mackie et al., 2003; Tims et al., 2013). Interestingly, shifts in related species of Bacteroides ovatus, B. plebeius and B. uniformis in this study were substrate dependent, which can be explained by the species-specific activity of pectin-degrading enzymes and the hierarchical preference of substrate utilization. Accordingly, Tuncil et al. (2017) demonstrated that human gut symbionts B. thetaiotamicron and B. ovatus, grown together on POS, had inverse growth profile, concurrently with the different expression pattern of glycan utilization genes. Preferential utilization of metabolites from pectin degradation can be expected, regarding the highly competitive environment in fecal fermentations. High increase in genera Butyrivibrio and Lachnospira, previously shown incapable to degrade pectins (van Laere et al., 2000; Flint et al., 2012), might be linked to metabolic cross-feeding interactions between the members of the fecal microbiota.

Bacterial Populations Associated With Human Health Were Selectively Stimulated by Pectins

Pectin fermentations affected abundances of various bacterial taxa associated with microbiota dysbiosis in human diseases, such as obesity and IBD. Thus, it was previously shown that obesity-related bias in the gut microbiota included higher levels of Blautia, Eubacterium, Roseburia, Dorea and Ruminococcus, along with reduction of F. prausnitzii, Oscillospira, Christensenellaceae, Prevotella, Bacteroides and genera within Proteobacteria (Nadal et al., 2009; Furet et al., 2010; Tims et al., 2013; Verdam et al., 2013; Goodrich et al., 2014; Kasai et al., 2015). Besides, R. torques and R. gnavus, decreased in this study, were found in higher numbers in IBD patients (Png et al., 2010). Notably, shifts in F. prausnitzii, Coprococcus, B. ovatus, B. plebeius, P. copri and Sutterella were strongly dependent on the feed, suggesting that these species can be modulated by specific pectins. Among them, F. prausnitzii, is commonly referred to as a marker for intestinal health, exhibiting anti-inflammatory effects in the gut (Walters et al., 2014). This study indicated that stimulation of F. prausnitzii could be achieved by fermentation of high methoxyl pectins (P1–P3) rather than LM pectins (P7 and P8). In their turn, LM citrus pectins were efficient to decrease proportions of P. copri, a microbe associated with induced insulin resistance in mice and rheumatoid arthritis in humans (Pedersen et al., 2016; Pianta et al., 2017); and to increase the levels of Coprococcus linked to reduced severity of IBS in humans (Tap et al., 2017).

Production of Propionate, Butyrate and Total SCFA Differed Between the Pectins

The overall production of propionate was highest in RGI fermentations and in fermentations of high methoxyl pectins (especially P4) and correlated positively with the relative abundances of P. copri, Ruminococcus spp, and unidentified Clostridiales, suggesting that these species were able to generate propionate from pectin fermentation products. Supporting our results, Gulfi and coworkers reported that HM pectins had a tendency to produce larger amounts of propionic acid in batch fecal fermentations, compared to LM pectins (Gulfi et al., 2005). Additionally, the ability of Prevotella and Ruminococcus to produce propionate was recently confirmed by genomic analysis of colonic anaerobes combined with growth experiments (Qin et al., 2010; Reichardt et al., 2014; Chen et al., 2017). Production of propionate by Ruminococcus is known to be enhanced in the presence of Rha and fucose (Reichardt et al., 2014), which are the common structural units of pectin molecules. Consequently, increased production of propionate in fermentations of sugar beet pectin (P4) and RGI (P10) might be related to the relatively higher content of Rha and other neutral sugars. Interestingly, the acetate to propionate ratios were found to be reduced in RGI and the high methoxyl pectin group, lowest for P4 and harsh extracted pectin from orange (P1). Lower acetate/propionate ratios have been associated with an anti-cholesterolemic effect and reduction of cardiovascular disease risk, and generally considered as beneficial (den Besten et al., 2013). The lowest amount of butyrate was detected in RGI fermentations, which might be explained by decreased levels of F. prausnitzii, a predominant butyrate producer in the gut (Louis and Flint, 2017). The negative correlation between butyrate and the numbers of Bacteroides spp. was, probably, related to their inability to produce butyrate (Louis and Flint, 2017). Production of SCFA in fermentations of citrus pectins (P1–P3 and P5–P9) was not statistically different, despite the differences in microbiota composition. This effect was presumably caused by the functional redundancy of microbial populations and the metabolic cross-feeding interactions (Reichardt et al., 2018).

Changes in Microbiota Were Related to the Structural Features of Pectins

Shifts in bacterial abundances in TIM-2 fermentations were related to the structural features of pectins. We identified at least five factors essential for microbiota shaping: (i) DE of polygalacturonic acid, (ii) composition of neutral sugars, (iii) distribution of HG and RG fractions, (iv) DBr, and (v) modification of pectic backbone, e.g., by amidation. DE is probably the most important parameter, as seen from the beta-diversity analysis and the correlation of bacterial taxa with DE. So far, only few studies indicated a link between the DE in pectins and microbiota composition. In agreement with our results, enrichment in Prevotella spp. was found in colonic microbiota of pigs fed with LM pectins, while genus Bacteroides was increased in fecal batch fermentations of LM pectin contrary to HM pectin (Dongowski et al., 2002; Tian et al., 2017). Results of the correlation analysis (Figure 5) pointed out on the possibility of differential stimulation of bacterial populations using pectins with different sugar content, especially Rha and Gal. Thus, correlation of F. prausnitzii with the major sugars (negative) and GalA/DE (positive) suggested, that HM pectins with high fraction of HG over RG would be preferable for growth of this beneficial microbe. Interestingly, the residue Ara differed from other sugars, correlating positively with P. copri and negatively with Coprococcus and Lachnospira. Both Lachnospira spp. and P. copri possess α-arabinofuranosidase activity (Dodd et al., 2011), having potential to degrade arabinan side chains in pectin molecules. Thus, inverse correlation is most probably linked to the differences between bacterial species in activity of arabinanolytic enzymes and the metabolic hierarchy as discussed above. Unexpectedly, we observed similarity in microbiota profiling between the structurally different substrates, i.e., sugar beet pectin and RGI. This similarity might be associated with HG/RG distribution, as both sugar beet pectin and RGI were distinguished by the low HG (GalA) and high RG content compared to other pectins. Likewise, grouping of the LM pectin P9 together with the high methoxyl pectins (P1 and P2) might be related to high percentage of HG in these substrates. Similarity between pectins P7 and P8 is most probably attributed to amidation of the C-6 uronate groups in addition to the equally low DE values. It is well-known that structural features of pectins affect their rheological behavior and functional properties, e.g., gelling capacity, viscosity, molecular conformation, and solubility (Sila et al., 2009). We suppose that these properties might have an impact on substrate-bacterial interactions and pectin utilization in microbiota fermentations. Furthermore, carboxyl groups in the polygalacturonic chains can be involved in electrostatic interactions and, together with amide groups, in hydrogen bonding. In the previous study, we observed that the surface charge (or zeta-potential) of the LM pectins, was lower than that of HM pectins, due to the higher fraction of non-esterified carboxyl groups (Larsen et al., 2018). Lower zeta-potential in LM pectins indicated stronger electrostatic repulsion upon interactions with bacterial cell wall groups (hydroxyl, carboxyl, etc.) negatively charged at neutral pH (Jiang et al., 2004). High DBr and RG fraction can either provide steric hindrance or contribute to hydrophobic interactions with non-polar groups in bacterial cell walls. Involvement in the different types of interactions can, probably, explain the opposite correlation found for DE (and GalA) compared to DBr and neutral sugars with bacterial taxa in this study.

In conclusion, this study provided evidence that modulation of the gut microbiota by pectins depended on their structural features. We identified specific bacterial taxa which abundances were differentially affected by pectins and proposed the main factors, linked to differences in microbiota composition. Understanding of the interplay between the gut commensals and the structural properties of pectins is essential to predict physiological effects of ingested pectins and to provide ideas for development of pectin-containing dietary fibers, targeting beneficial bacteria to facilitate more balanced microbiota profiles. It would be relevant to verify our findings in relation to other gut microbiotas, e.g., from diseased subjects. Additionally, further comparative in vitro and in vivo studies with structurally diverse pectins and their derivatives are needed to achieve the detailed knowledge of structure-function relationship of pectins in the gut.

Funding. This research was funded by the Danish Council for Strategic Research (DSR, Project BioSyn, No. 3050-00005B) within the frame of the project “Strategic Research Collaboration in Food Science in the State of São Paulo, Brazil and Denmark” and supported by the Brazilian National Council for Scientific and Technological Development (CNPq-Brazil), the Program “Science without Borders.” From the CNPq-Brazil CBS received a scholarship 246027/2012-6, and TC – scholarship 233521/2014-3. This research has additionally been supported by the Dutch Province of Limburg.

• Abbott D. W., Boraston A. B. (2008). Structural biology of pectin degradation by Enterobacteriaceae.Microbiol. Mol. Biol. Rev.72301–316. 10.1128/MMBR.00038-07 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Aguirre M., Jonkers D. M. A. E., Troost F. J., Roeselers G., Venema K. (2014). In vitro characterization of the impact of different substrates on metabolite production, energy extraction and composition of gut microbiota from lean and obese subjects.PLoS One9:e113864. 10.1371/journal.pone.0113864 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Caporaso J. G., Kuczynski J., Stombaugh J., Bittinger K., Bushman F. D., Costello E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data.Nat. Methods7335–336. 10.1038/nmeth.f.303 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Chan S. Y., Choo W. S., Young D. J., Loh X. J. (2017). Pectin as a rheology modifier: origin, structure, commercial production and rheology.Carbohydr. Polym.161118–139. 10.1016/j.carbpol.2016.12.033 [PubMed] [CrossRef] [Google Scholar]
• Chen J., Liang R. H., Liu W., Li T., Liu C. M., Wu S. S., et al. (2013). Pectic-oligosaccharides prepared by dynamic high-pressure microfluidization and their in vitro fermentation properties.Carbohydr. Polym.91175–182. 10.1016/j.carbpol.2012.08.021 [PubMed] [CrossRef] [Google Scholar]
• Chen T., Long W., Zhang C., Liu S., Zhao L., Hamaker B. R. (2017). Fiber-utilizing capacity varies in Prevotella- versus Bacteroides-dominated gut microbiota.Sci. Rep.7:2594. 10.1038/s41598-017-02995-4 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• den Besten G., van Eunen K., Groen A. K., Venema K., Reijngoud D. J., Bakker B. M. (2013). The role of short-chain fatty acids in the interplay between diet, gut microbiota, and host energy metabolism.J. Lipid Res.542325–2340. 10.1194/jlr.R036012 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Dodd D., Mackie R. I., Cann I. K. O. (2011). Xylan degradation, a metabolic property shared by rumen and human colonic Bacteroidetes.Mol. Microbiol.79292–304. 10.1111/j.1365-2958.2010.07473.x [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Dongowski G., Lorenz A., Proll J. (2002). The degree of methylation influences the degradation of pectin in the intestinal tract of rats and in vitro.J. Nutr.1321935–1944. 10.1093/jn/132.7.1935 [PubMed] [CrossRef] [Google Scholar]
• Flint H. J., Scott K. P., Duncan S. H., Louis P., Forano E. (2012). Microbial degradation of complex carbohydrates in the gut.Gut Microbes3286–306. 10.4161/gmic.19897 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• f*ckunaga T., Sasaki M., Araki Y., Okamoto T., Yasuoka T., Tsujikawa T., et al. (2003). Effects of the soluble fibre pectin on intestinal cell proliferation, fecal short chain fatty acid production and microbial population.Digestion6742–49. 10.1159/000069705 [PubMed] [CrossRef] [Google Scholar]
• Furet J.-P., Kong L.-C., Tap J., Poitou C., Basdevant A., Bouillot J.-L., et al. (2010). Differential adaptation of human gut microbiota to bariatric surgery–induced weight loss: links with metabolic and low-grade inflammation markers.Diabetes123049–3057. 10.2337/db10-0253 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Gómez B., Gullón B., Remoroza C., Schols H. A., Parajó J. C., Alonso J. L. (2014). Purification, characterization, and prebiotic properties of pectic oligosaccharides from orange peel wastes.J. Agric. Food Chem.629769–9782. 10.1021/jf503475b [PubMed] [CrossRef] [Google Scholar]
• Gómez B., Gullón B., Yáñez R., Schols H., Alonso J. L. (2016). Prebiotic potential of pectins and pectic oligosaccharides derived from lemon peel wastes and sugar beet pulp: a comparative evaluation.J. Funct. Foods20108–121. 10.1016/j.jff.2015.10.029 [CrossRef] [Google Scholar]
• Goodrich J. K., Waters J. L., Poole A. C., Sutter J. L., Koren O., Blekhman R., et al. (2014). Human genetics shape the gut microbiome.Cell159789–799. 10.1016/j.cell.2014.09.053 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Grondin J. M., Tamura K., Déjean G., Abbott D. W., Brumer H. (2017). Polysaccharide utilization loci: fueling microbial communities.J. Bacteriol.199e860–e816. 10.1128/JB.00860-16 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Gulfi M., Arrigoni E., Amadò R. (2005). Influence of structure on in vitro fermentability of commercial pectins and partially hydrolysed pectin preparations.Carbohydr. Polym.59247–255. 10.1016/j.carbpol.2004.09.018 [CrossRef] [Google Scholar]
• Jiang W., Saxena A., Song B., Ward B. B., Beveridge T. J., Myneni S. C. B. (2004). Elucidation of functional groups on gram-positive and gram-negative bacterial surfaces using Infrared spectroscopy.Langmuir2011433–11442. 10.1021/la049043+ [PubMed] [CrossRef] [Google Scholar]
• Kasai C., Sugimoto K., Moritani I., Tanaka J., Oya Y., Inoue H., et al. (2015). Comparison of the gut microbiota composition between obese and non-obese individuals in a Japanese population, as analyzed by terminal restriction fragment length polymorphism and next-generation sequencing.BMC Gastroenterol.15:100. 10.1186/s12876-015-0330-2 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Larsen N., Cahú T. B., Marta S., Saad I., Blennow A., Jespersen L. (2018). The effect of pectins on survival of probiotic Lactobacillus spp. in gastrointestinal juices is related to their structure and physical properties.Food Microbiol.7411–20. 10.1016/j.fm.2018.02.015 [PubMed] [CrossRef] [Google Scholar]
• Leijdekkers A. G. M., Aguirre M., Venema K., Bosch G., Gruppen H., Schols H. A. (2014). In vitro fermentability of sugar beet pulp derived oligosaccharides using human and pig fecal inocula.J. Agric. Food Chem.621079–1087. 10.1021/jf4049676 [PubMed] [CrossRef] [Google Scholar]
• Li P., Xia J., Nie Z., Shan Y. (2016). Pectic oligosaccharides hydrolyzed from orange peel by fungal multi-enzyme complexes and their prebiotic and antibacterial potentials.LWT Food Sci. Technol.69203–206. 10.1016/j.lwt.2016.01.042 [CrossRef] [Google Scholar]
• Lopez-Siles M., Khan T. M., Duncan S. H., Harmsen H. J. M., Garcia-Gil L. J., Flint H. J. (2012). Cultured representatives of two major phylogroups of human colonic Faecalibacterium prausnitzii can utilize pectin, uronic acids, and host-derived substrates for growth.Appl. Environ. Microbiol.78420–428. 10.1128/AEM.06858-11 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Louis P., Flint H. J. (2017). Formation of propionate and butyrate by the human colonic microbiota.Environ. Microbiol.1929–41. 10.1111/1462-2920.13589 [PubMed] [CrossRef] [Google Scholar]
• Mackie R. I., Aminov R. I., Hu W., Klieve A. V., Ouwerkerk D., Sundset M. A., et al. (2003). Ecology of uncultivated Oscillospira species in the rumen of cattle, sheep, and reindeer as assessed by microscopy and molecular approaches.Appl. Environ. Microbiol.696808–6815. 10.1128/AEM.69.11.6808-6815.2003 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Manderson K., Pinart M., Tuohy K. M., Grace W. E., Hotchkiss A. T., Widmer W., et al. (2005). In vitro determination of prebiotic properties of oligosaccharides derived from an orange juice manufacturing by-product stream.Appl. Environ. Microbiol.718383–8389. 10.1128/AEM.71.12.8383-8389.2005 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Martens E. C., Lowe E. C., Chiang H., Pudlo N. A., Wu M., McNulty N. P., et al. (2011). Recognition and degradation of plant cell wall polysaccharides by two human gut symbionts.PLoS Biol.9:e1001221. 10.1371/journal.pbio.1001221 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Mcdonald D., Price M. N., Goodrich J., Nawrocki E. P., Desantis T. Z., Probst A., et al. (2012). An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea.ISME J.6610–618. 10.1038/ismej.2011.139 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Nadal I., Santacruz A., Marcos A., Warnberg J., Garagorri M., Moreno L., et al. (2009). Shifts in clostridia, Bacteroides and immunoglobulin- coating fecal bacteria associated with weight loss in obese adolescents.Int. J. Obes.33758–767. 10.1038/ijo.2008.260 [PubMed] [CrossRef] [Google Scholar]
• Nakajima N., Ishihara K., Matsuura Y. (2002). Dietary-fiber-degrading enzymes from a human intestinal Clostridium and their application to oligosaccharide production from nonstarchy polysaccharides using immobilized cells.Appl. Microbiol. Biotechnol.59182–189. 10.1007/s00253-002-1015-7 [PubMed] [CrossRef] [Google Scholar]
• Olano-Martin E., Gibson G. R., Rastall R. A. (2002). Comparison of the in vitro bifidogenic properties of pectins and pectic-oligosaccharides.J. Appl. Microbiol.93505–511. 10.1046/j.1365-2672.2002.01719.x [PubMed] [CrossRef] [Google Scholar]
• Onumpai C., Kolida S., Bonnin E., Rastall R. A. (2011). Microbial utilization and selectivity of pectin fractions with various structures.Appl. Environ. Microbiol.775747–5754. 10.1128/AEM.00179-11 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Parkar S. G., Redgate E. L., Wibisono R., Luo X., Koh E. T. H., Schröder R. (2010). Gut health benefits of kiwifruit pectins: comparison with commercial functional polysaccharides.J. Funct. Foods2210–218. 10.1016/j.jff.2010.04.009 [CrossRef] [Google Scholar]
• Pedersen H. K., Gudmundsdottir V., Nielsen H. B., Hyotylainen T., Nielsen T., Jensen B. A. H., et al. (2016). Human gut microbes impact host serum metabolome and insulin sensitivity.Nature535376–381. 10.1038/nature18646 [PubMed] [CrossRef] [Google Scholar]
• Pianta A., Arvikar S., Strle K., Drouin E. E., Wang Q., Costello C. E., et al. (2017). Evidence of the immune relevance of prevotella copri, a gut microbe, in patients with rheumatoid arthritis.Arthritis Rheumatol.69964–975. 10.1002/art.40003 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Png C. W., Lindén S. K., Gilshenan K. S., Zoetendal E. G., McSweeney C. S., Sly L. I., et al. (2010). Mucolytic bacteria with increased prevalence in IBD mucosa augment in vitro utilization of mucin by other bacteria.Am. J. Gastroenterol.1052420–2428. 10.1038/ajg.2010.281 [PubMed] [CrossRef] [Google Scholar]
• Qin J., Li R., Raes J., Arumugam M., Burgdorf K. S., Manichanh C., et al. (2010). A human gut microbial gene catalogue established by metagenomic sequencing.Nature46459–65. 10.1038/nature08821 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Regulation (EC) (2008). Regulation (EC) No 1333/2008 of the European Parliament, and of the Council of 16 December 2008 on food additives. Available at: http://data.europa.eu/eli/reg/2008/1333/oj. [Google Scholar]
• Reichardt N., Duncan S. H., Young P., Belenguer A., Leitch C. M., Scott K. P., et al. (2014). Phylogenetic distribution of three pathways for propionate production within the human gut microbiota.ISME J.81323–1335. 10.1038/ismej.2014.14 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Reichardt N., Vollmer M., Holtrop G., Farquharson F. M., Wefers D., Bunzel M., et al. (2018). Specific substrate-driven changes in human faecal microbiota composition contrast with functional redundancy in short-chain fatty acid production.ISME J.12610–622. 10.1038/ismej.2017.196 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Schwartz S. E., Levine R. A., Weinstock R. S., Petokas S., Mills C. A., Thomas F. D. (1988). Sustained pectin ingestion: effect on gastric emptying and glucose tolerance in non-insulin-dependent diabetic patients.Am. J. Clin. Nutr.481413–1417. 10.1093/ajcn/48.6.1413 [PubMed] [CrossRef] [Google Scholar]
• Shinohara K., Ohashi Y., Kawasumi K., Terada A., Fujisawa T. (2010). Effect of apple intake on fecal microbiota and metabolites in humans.Anaerobe16510–515. 10.1016/j.anaerobe.2010.03.005 [PubMed] [CrossRef] [Google Scholar]
• Sila D. N., Van Buggenhout S., Duvetter T., Van Loey A., Hendrickx M. (2009). Pectins in processed fruits and vegetables: part II-structure function relationships.Compr. Rev. Food Sci. Food Saf.8105–117. 10.1111/j.1541-4337.2009.00071.x [CrossRef] [Google Scholar]
• Sulek K., Vigsnaes L. K., Schmidt L. R., Holck J., Frandsen H. L., Smedsgaard J., et al. (2014). A combined metabolomic and phylogenetic study reveals putatively prebiotic effects of high molecular weight arabino-oligosaccharides when assessed by invitro fermentation in bacterial communities derived from humans.Anaerobe2868–77. 10.1016/j.anaerobe.2014.05.007 [PubMed] [CrossRef] [Google Scholar]
• Tap J., Derrien M., Törnblom H., Brazeilles R., Cools-portier S., Doré J., et al. (2017). Identification of an intestinal microbiota signature associated with severity of irritable bowel syndrome.Gastroenterology152111–123. 10.1053/j.gastro.2016.09.049 [PubMed] [CrossRef] [Google Scholar]
• Thakur B. R., Singh R. K., Handa A. K. (1997). Chemistry and Uses of Pectin - A Review.Crit. Rev. Food Sci. Nutr.3747–73. 10.1080/10408399709527767 [PubMed] [CrossRef] [Google Scholar]
• Tian L., Bruggeman G., van den Berg M., Borewicz K., Scheurink A. J. W., Bruininx E., et al. (2017). Effects of pectin on fermentation characteristics, carbohydrate utilization, and microbial community composition in the gastrointestinal tract of weaning pigs.Mol. Nutr. Food Res.61:1600186. 10.1002/mnfr.201600186 [PubMed] [CrossRef] [Google Scholar]
• Tims S., Derom C., Jonkers D. M., Vlietinck R., Saris W. H., Kleerebezem M., et al. (2013). Microbiota conservation and BMI signatures in adult monozygotic twins.ISME J.7707–717. 10.1038/ismej.2012.146 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Tuncil Y. E., Nakatsu C. H., Kazem A. E., Arioglu-Tuncil S., Reuhs B., Martens E. C., et al. (2017). Delayed utilization of some fast-fermenting soluble dietary fibers by human gut microbiota when presented in a mixture.J. Funct. Foods32347–357. 10.1016/j.jff.2017.03.001 [CrossRef] [Google Scholar]
• van Laere K. M. J., Hartemink R., Bosveld M., Schols H. A., Voragen A. G. J. (2000). Fermentation of plant cell wall derived polysaccharides and their corresponding oligosaccharides by intestinal bacteria.J. Agric. Food Chem.481644–1652. 10.1021/jf990519i [PubMed] [CrossRef] [Google Scholar]
• Venema K. (2015). “The TNO in vitro model of the colon ( TIM-2),” inThe Impact of Food Bioactives on Health: In Vitro and Ex Vivo Models, edsVerhoeckx K., Cotter P., Kleiveland I. C., Lea T., Mackie A.et al. (Cham: Springer; ), 293–304. 10.1007/978-3-319-16104-4 [CrossRef] [Google Scholar]
• Verdam F. J., Fuentes S., De Jonge C., Zoetendal E. G., Erbil R., Greve J. W., et al. (2013). Human intestinal microbiota composition is associated with local and systemic inflammation in obesity.Obesity21E607–E615. 10.1002/oby.20466 [PubMed] [CrossRef] [Google Scholar]
• Walters W. A., Xu Z., Knight R. (2014). Meta-analyses of human gut microbes associated with obesity and IBD.FEBS Lett.5884223–4233. 10.1016/j.febslet.2014.09.039 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Wicker L., Kim Y., Kim M. J., Thirkield B., Lin Z., Jung J. (2014). Pectin as a bioactive polysaccharide - extracting tailored function from less.Food Hydrocoll.42251–259. 10.1016/j.foodhyd.2014.01.002 [CrossRef] [Google Scholar]
• Williams A. R., Hansen T. V. A., Krych L., Ahmad H. F. B., Nielsen D. S., Skovgaard K., et al. (2017). Dietary cinnamaldehyde enhances acquisition of specific antibodies following helminth infection in pigs.Vet. Immunol. Immunopathol.18943–52. 10.1016/j.vetimm.2017.06.004 [PubMed] [CrossRef] [Google Scholar]
• Yapo B. M., Koffi K. L. (2014). Extraction and characterization of highly gelling low methoxy pectin from cashew apple pomace.Foods31–12. 10.3390/foods3010001 [PMC free article] [PubMed] [CrossRef] [Google Scholar]