# Design and Investigation of PolyFermS In Vitro Continuous Fermentation Models Inoculated with Immobilized Fecal Microbiota Mimicking the Elderly Colon

Sophie Fehlbaum, Christophe Chassard, Martina C. Haug, Candice Fourmestraux, Muriel Derrien, Christophe Lacroix

► To cite this version:

Sophie Fehlbaum, Christophe Chassard, Martina C. Haug, Candice Fourmestraux, Muriel Derrien, et al.. Design and Investigation of PolyFermS In Vitro Continuous Fermentation Models Inoculated with Immobilized Fecal Microbiota Mimicking the Elderly Colon. PLoS ONE, 2015, 10 (11), pp.e0142793. 10.1371/journal.pone.0142793. hal-04658070

HAL Id: hal-04658070

https://hal.science/hal-04658070v1

Submitted on 23 Jul 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Design and Investigation of PolyFermS *In Vitro* Continuous Fermentation Models Inoculated with Immobilized Fecal Microbiota Mimicking the Elderly Colon

Sophie Fehlbaum¹, Christophe Chassard¹, Martina C. Haug¹, Candice Fourmestraux², Muriel Derrien², Christophe Lacroix¹*

1 Laboratory of Food Biotechnology, Institute of Food, Nutrition and Health, ETH Zurich, Zurich, Switzerland,

2 Danone Nutricia Research, Palaiseau, France

* christophe.lacroix@hest.ethz.ch

<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Fehlbaum S, Chassard C, Haug MC, Fourmestraux C, Derrien M, Lacroix C (2015) Design and Investigation of PolyFermS *In Vitro* Continuous Fermentation Models Inoculated with Immobilized Fecal Microbiota Mimicking the Elderly Colon. PLoS ONE 10(11): e0142793. doi:10.1371/journal. pone.0142793

**Editor:** Tiffany L. Weir, Colorado State University, UNITED STATES

**Received:** April 8, 2015

**Accepted:** October 27, 2015

**Published:** November 11, 2015

**Copyright:** © 2015 Fehlbaum et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All 454 pyrosequencing files have been deposited to the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under accession number SRP053000.

**Funding:** Funding for this study was provided by Danone Research (Palaiseau, France). The funder supported the salary for SF. The authors MD and CF are employees of Danone Research. The specific roles of these authors are articulated in the 'author contributions' section.

## Abstract

*In vitro* gut modeling is a useful approach to investigate some factors and mechanisms of the gut microbiota independent of the effects of the host. This study tested the use of immobilized fecal microbiota to develop different designs of continuous colonic fermentation models mimicking elderly gut fermentation. Model 1 was a three-stage fermentation mimicking the proximal, transverse and distal colon. Models 2 and 3 were based on the new PolyFermS platform composed of an inoculum reactor seeded with immobilized fecal microbiota and used to continuously inoculate with the same microbiota different second-stage reactors mounted in parallel. The main gut bacterial groups, microbial diversity and metabolite production were monitored in effluents of all reactors using quantitative PCR, 16S rRNA gene 454-pyrosequencing, and HPLC, respectively. In all models, a diverse microbiota resembling the one tested in donor's fecal sample was established. Metabolic stability in inoculum reactors seeded with immobilized fecal microbiota was shown for operation times of up to 80 days. A high microbial and metabolic reproducibility was demonstrated for downstream control and experimental reactors of a PolyFermS model. The PolyFermS models tested here are particularly suited to investigate the effects of environmental factors, such as diet and drugs, in a controlled setting with the same microbiota source.

## Introduction

The human colon harbors a large number of microbes forming a complex ecosystem responsible for various processes in the host. Under normal conditions, the gut microbiota acts as a barrier against enteropathogens, contributes to the development of the immune system and exerts important metabolic functions; which includes the production of short chain fatty acids (SCFA; such as acetate, propionate and butyrate) by breaking down complex carbohydrates that provide energy to epithelial cells and to the host [1, 2]. Each human harbors a unique gut
**Competing Interests:** MD and CF are employees of Danone Research. This does not alter the authors' adherence to PLOS ONE policies on sharing data and materials.

microbiota composition consisting of bacteria belonging mainly to the phyla Firmicutes or Bacteroidetes and, to a lesser extent, to Actinobacteria, Proteobacteria and Verrucomicrobia [3, 4]. Colonization of the gut occurs first during birth, and throughout the first 2–3 years of life the microbial composition becomes established towards an adult-like microbiota. Recent studies indicate that the gut microbiota remains stable in adulthood, except for temporary alterations due to diet, disease and antibiotic treatment. However, an important shift in the microbial composition occurs during old age that is associated with a reduction in stability and often in biodiversity [3–7]. Additionally, a large inter-individual variability of the gut microbiota composition was reported for elderly Irish subjects of a large-scale *in vivo* study, with pyrosequencing reads assigned to the phyla *Bacteroidetes* and *Firmicutes* ranging from 3 to 94% [8]. To date, no common core microbiota for the elderly was defined, partly due to the various physiological factors, including lifestyle, diet and need for medications that change in old age. Thus, establishing the changes in the composition with ageing still require further investigations [4, 7].

Intestinal fermentation models allow the *in vitro* cultivation of gut microbiota to study their composition and function, uncoupled from the host. As such, models provide greater control, easier manipulation, and no ethical restrictions relative to *in vivo* studies, and are very complementary to *in vivo* strategies for elucidating mechanisms of gut microbiota [9]. Intestinal models have developed from batch for short-term fermentation studies to single or multistage continuous models that allow long-term studies due to substrate replenishment and toxic product removal [10]. However, one of the main challenges of the continuous culture models is the reproduction of the biofilm-associated microbes of the gut that is important to prevent washout of the less competitive bacteria. Immobilization of gut microbiota in gellan-xanthan gel beads has shown to reproduce the free and biofilm associated states of bacterial populations and to maintain the bacterial diversity at high cell densities in continuous intestinal reactors over periods of up to 71 days [10–12]. Furthermore, reproducibility and biological replication of continuous intestinal models was recently improved by the introduction of the PolyFermS model that allows the parallel testing of treatments with the same gut microbiota, and which has been validated for the child and the swine proximal colon [13, 14].

*In vitro* intestinal fermentations models have been developed and validated [10] to investigate factors of microbiota composition and metabolism of infants to adults while the elderly gut microbiota was only scarcely analyzed. Several studies were performed in continuous three-stage models for investigating the effects of antibiotics on *Clostridium difficile* infection. For these studies *C. difficile* was inoculated with mixed fecal samples from multiple elder donors and the system was challenged with antibiotics to promote the germination and growth of the sporulated bacteria, while microbiota analysis was only done with cultivation [15–18]. In a recent study batch cultures and continuous three-stage models inoculated with microbiota from single fecal samples of elder donors were used to investigate probiotics, prebiotics and synbiotics. Fluorescent *in situ* hybridization methods were used to monitor gut microbiota composition [19]. To date, no study has reported a detailed analysis of gut microbiota establishment and diversity in *in vitro* fermentation models reproducing the gut of aged (over 65 years) people.

The aim of this study was to investigate the use of immobilized fecal microbiota to develop different designs of continuous colonic fermentation models mimicking elderly gut fermentation. Immobilization of fecal microbiota obtained from three different donors was performed independently. Fecal beads were used to inoculate an immobilized cell reactor (IR) operated with conditions selected to mimic the proximal colon section of an elder. Three model designs, all starting with an IR used to generate a constant gut microbiota composition in proximal colon conditions, were tested for different experimental questions. These models were set in
sequential order with the aim to investigate colonization of *Clostridium difficile* (data not shown in this paper). Model 1 was based on the three-stage design, with immobilized microbiota inoculated in a first proximal colon reactor connected to a transverse and a distal colon reactor, previously validated for infant and child microbiota fermentation [12, 20, 21]. Models 2 and 3 were developed based on the PolyFermS platform, which recently has been validated with child [13] and swine [14] microbiota, with adjusted conditions for the elderly microbiota. In model 2 IR containing immobilized fecal microbiota was used to inoculate (10% v/v) two parallel sets of 2-stage reactors mimicking the proximal and distal colon. Because *C. difficile* growth was only detected in distal colon reactors, in model 3 IR was used to feed (100% v/v) five reactors mounted in parallel, and mimicking conditions of the distal colon. This design allowed to test in parallel four treatments compared to a control in distal colon reactors. The microbiota composition in reactor effluents was monitored and compared to that of the corresponding fecal donor, and temporal stability of the models and reproducibility of downstream reactors within a PolyFermS model were demonstrated. Microbiota composition, diversity (qPCR and pyrosequencing) and activity (HPLC) were monitored in reactor effluents over operation periods of up to 80 days.

# Materials and Methods

## Ethics Statement

The Ethics Committee of ETH Zurich exempted this study from review because sample collection was not in terms of intervention. An informed written consent was, however, obtained from the fecal donors.

## Fecal inoculum and immobilization

For each fermentation experiment a fresh fecal sample from a different donor was used for the immobilization procedure. Fecal samples were collected from three healthy women, aged 71 (fermentation 1), 72 (fermentation 2) and 78 years (fermentation 3), who did not receive anti-biotic treatment for at least three months prior to sample collection, and who did not consume probiotics on a regular basis.

Immediately after defecating, the fecal sample was transferred to a tube containing 5 mL of sterile, pre-reduced peptone water (0.1%, pH 7), placed in an anaerobic jar (Anaerojar, Oxoid, Hampshire, England), and transported and processed within three hours. Handling and encapsulation of the fecal microbiota into 1–2 mm gel beads composed of gellan (2.5% w/v), xanthan (0.25% w/v), and sodium citrate (0.2% w/v, Sigma-Aldrich Chemie GmbH, Buchs, Switzerland) was performed in an anaerobic chamber as previously described [21].

## Fermentation medium

The fermentation medium was based on the composition described by MacFarlane *et al.* [22] for simulation of adult chyme entering the colon. It contained (g L⁻¹ of distilled water): pectin from citrus (2), xylan from oat spelt (2), arabinogalactan from larch wood (2), guar gum (1), inulin (1), soluble potato starch (5), mucin (4), casein acid hydrolysate (3), peptone water (5), tryptone (5), yeast extract (4.5), cysteine (0.8), bile salts (0.4), KH₂PO₄ (0.5), NaHCO₃ (1.5), NaCl (4.5), KCl (4.5), MgSO₄ anhydrous (0.6), CaCl₂ x 2H₂O (0.1), MnCl₂ x 4H₂O (0.2), FeSO₄ x 7H₂O (0.005), hemin (0.05) and Tween 80 (1). One mL of a filter-sterilized (0.2 µm pore-size) vitamin solution [23] was added to 1 L of autoclaved (20 min, 120°C) and cooled medium. All components of the nutritive medium were purchased from Sigma-Aldrich Chemie, except for inulin (Orafti®, BENEO kindly provided by RPN Foodtechnology AG, Sursee,
Switzerland), peptone water (Oxoid AG, Pratteln, Switzerland), bile salts (Oxoid AG), tryptone (Becton Dickinson AG, Allschwil, Switzerland) and KH₂PO₄ (VWR International AG).

## Experimental setup

Various reactor set-ups were applied for the three continuous intestinal fermentation experiments (Fig 1). The inoculum reactor (IR) seeded with donor's microbiota immobilized in polysaccharide gel beads, and operated with proximal colon conditions, was common to all models. Model 1 was a classical three-stage system consisting of three reactors placed in series and operated under conditions of the proximal (PC corresponding to IR), transverse (TC) and distal colon (DC) [21]. Model 2 was based on the recently developed PolyFermS platform [13, 14] and consisted of an IR with immobilized fecal microbiota in proximal colon conditions used to continuously inoculate two parallel systems (10% v/v of the feed), each composed of a proximal (PC1 and PC2) and a transverse-distal reactor (DC1 and DC2). For model 3, the chyme medium fermented in IR with immobilized microbiota and operated with proximal colon conditions was used to continuously feed (100% v/v of the feed) five reactors mounted in parallel and mimicking conditions of a transverse-distal colon. In models 2 and 3, one system or reactor downstream to IR was used as control while the other system (PC2-DC2) or reactors (TR1-TR4) were used to comparatively test treatments, respectively.

## Fermentation procedures

The IR inoculated with 30% v/v gel beads was used to first colonize the beads with repeated-batch fermentations. The fresh medium was replaced every 12 h, for a total fermentation time of 60 to 72 h, depending on the model. Temperature was set at 37°C, stirring speed at 120 rpm and the pH was controlled at 5.5 or 5.7 by addition of 2.5 M NaOH. Sterile nutritive medium (4°C) was pumped continuously via a peristaltic pump (Reglo analog, Ismatec, Glattbrugg, Switzerland). Total mean retention times between 28 to 38.5 h were used in models 1–3, which is in the range of previously measured values for total colonic transit time in healthy elderly subjects that were between 25 and 66 h [24–28].

The reactors and nutritive media were continuously flushed with a low flow of CO₂ to maintain anaerobic conditions during fermentation. The conditions of the models along with the model specific design and trials are briefly presented below and summarized in Fig 1.

**Model 1.** Six repeated-batch fermentations (72 h) were performed to colonize beads in PC (fermentation volume of 200 mL, pH of 5.5). Then two reactors mimicking conditions of transverse (TC) and distal colon (DC) (400 mL, pH 6.2 and 6.8, respectively) were connected in series to PC. The feed flow rate of the nutritive medium was set at 26 mL h⁻¹, giving mean retention times of 7.7 h for IR and 15.4 h for TC and DC, for a total system retention time of 38.5 h. The model was stabilized for 14 days.

**Model 2.** Five repeated-batch fermentations (60 h) were performed to colonize beads in IR (225 mL and pH 5.7). The model was then switched to continuous mode with a medium flow rate of 25 mL h⁻¹ for an additional five days. Two proximal colon reactors (PC1 and PC2, 250 mL and pH 5.7) were attached to IR, and each PC was connected to distal colon reactors (DC1 and DC2, 450 mL and pH 6.8) used to mimic transverse and distal colon conditions. Compared with model 1, the pH control set-point was increased from 5.5 to 5.7. This modification was implemented in order to better match the pH of the proximal colon section *in vivo*, which has been reported to be in the range from 5.5 to 5.9 [29] and to enhance metabolic activity in IR and PC's. PC1 and PC2 were continuously inoculated with 2.5 mL h⁻¹ (10%) fermented medium from IR and 22.5 mL h⁻¹ (90%) fresh nutritive medium, for total flow rate of 25 mL h⁻¹. Mean retention times in PC and DC reactors were 10 and 18 h, respectively, for a
A

<figure><img src="image_4.png" /></figure>

B

<figure><img src="image_5.png" /></figure>

C

<figure><img src="image_6.png" /><figcaption>**Fig 1.** Set-up of the continuous fermentation models with immobilized gut microbiota. (A) **Model 1.** 3-stage model consisting of a proximal, transverse and distal colon reactor (B) **Model 2.** 2-stage model with an inoculum reactor connected to two parallel test systems consisting of a proximal and distal colon reactor (C) **Model 3.** 2-stage model with an inoculum reactor (proximal colon conditions) feeding 5 distal colon reactors connected in parallel; RT: Retention time, V: Volume.</figcaption></figure>

doi:10.1371/journal.pone.0142793.g001
total retention time of 28 h (PC1 and DC1; PC2 and DC2). The model was continuously operated for an additional 18 days to reach stability and used for testing, for a total fermentation time of 55 days.

**Model 3.** Beads colonization was carried out in five repeated-batch cultures (60 h) in IR (450 mL and pH 5.7). Five DC reactors (225 mL and pH 6.8) mounted in parallel and mimicking conditions of transverse-distal colon were connected to IR. The fresh nutritive medium flow rate in IR was set at 45 mL h⁻¹, while each DC reactors was fed with 9 mL h⁻¹ medium fermented in IR. The mean retention times in IR and DC reactors of model 3 were 10 and 25 h, respectively, for a total retention time of 35 h. The model was stabilized for 14 days and used for testing, for a total fermentation of 80 days.

## Sampling and analysis

During continuous fermentations, effluent samples (10 mL) were collected daily from each reactor. Because the different models were also used for experimental trials with *Clostridium difficile* and antibiotic and probiotic treatments, only samples obtained during periods of control conditions are reported for model assessment. Analyses of microbial composition by quantitative polymerase chain reaction (qPCR) and pyrosequencing (model 2 and 3) were performed on samples from three days at the end of stabilization: days 9, 11 and 13 of model 1, days 16, 17 and 18 of model 2 and days 10, 11 and 12 of model 3. Metabolite concentrations in effluents of all reactors were tested daily during the entire fermentation by high-performance liquid chromatography (HPLC). Data corresponding to control condition periods (no treatment applied) are reported. Long term temporal stability was tested with IR data from models 2 and 3 since these reactors was not subjected to any manipulation over the entire culture period.

## DNA extraction

For qPCR and pyrosequencing analyses total microbial DNA of 200 mg feces and 2 mL effluent samples was extracted using the FastDNA® SPIN Kit for Soil (MP Biomedicals, Illkirch, France) and a final elution volume of 100 µL. DNA concentrations were determined using a Nanodrop® ND-1000 Spectrophotometer (Witec AG, Littau, Switzerland).

### qPCR analysis

Total bacteria and predominant bacterial groups were enumerated using specific primers (S1 Table). One µL of 10- or 100-fold diluted DNA was amplified in a total volume of 25 µL as described in [21], using 2 x SYBR Green PCR Master Mix (Applied Biosystems, Zug, Switzerland). Each reaction was run in duplicate on an ABI PRISM 7500-PCR sequence detection system (Applied Biosystems). For quantification, standard curves were produced by amplification of the DNA of the reference strain of the respective target group [30].

### HPLC analysis

SCFA (acetate, propionate, butyrate, valerate, isobutyrate and isovalerate) as well as formate and lactate concentrations) in fermentation effluent samples from all reactors were determined by HPLC analysis (Thermo Fisher Scientific Inc. Accela, Wohlen, Switzerland) in duplicate [31]. Effluent supernatants were 2-fold diluted with sterile ultra-pure water and filtered directly into vials through a 0.45 µm nylon HPLC filter (Infochroma AG, Zug, Switzerland). The analysis was run at a flow rate of 0.4 mL min⁻¹ using an Aminex HPX-87H column (Bio-Rad Laboratories AG, Reinach, Switzerland) and 10 mM H₂SO₄ as eluent.
Microbiota profiling by 454 pyrosequencing

454-pyrosequencing analysis of total genomic DNA of fecal and effluent samples was carried
out at DNAVision (Gosselies, Belgium). The V5-V6 hypervariable 16S RNA region was ampli-
fied using specific primers 784F (5'- AGGATTAGATACCCTKGTA-3') and 1061R (5'-CRR
CACGAGCTGACGAC-3') [32]. The forward primer contained the sequence of the Titanium
A adaptor and a unique barcode sequence. Pyrosequencing was carried out using primer A on
a 454 Life Sciences Genome Sequencer FLX instrument (Roche Applied Science, Vilvoorde,
Belgium) following Titanium chemistry. The data obtained was analyzed using the open source
software package Quantitative Insights Into Microbial Ecology (QIIME), v1.7 [33]. Raw
sequencing reads were filtered based on selected quality criteria such as: (1) no mismatch with
the primer sequences and barcode tags; (2) no ambiguous bases (Ns); (3) read-lengths not
shorter than 200 base pairs (bp) or longer than 1000 bp; (4) the average quality score in a
sliding window of 50 bp not to fall below 25; (5) excluding homopolymer runs higher than
6 nt. Sequences that passed quality filtering were clustered into OTUs at 97% identity level
using cd-hit [34]. Representative sequences (the most abundant) for each OTU were aligned
using PyNAST and taxonomically assigned using Greengenes v_13_08 database. Chimera-
Slayer was used to discard chimeric sequences, based on a reference data set of sequences [35].
This led to 8245 +/- 1924 (mean +/- SD) reads per sample. These phylogenies were combined
with absence/presence or abundance information for each OTU to calculate unweighted or
weighted UniFrac distances, respectively, using rarefaction of 7000 sequences per samples.
Unifrac measures the phylogenetic distance between sets of taxa in a phylogenetic tree as the
fraction of the branch length of the tree that leads to descendants from either one environment
or the other, but not both [36, 37]. Weighted and unweighted Unifrac metrics were used to
build phylogenetic distance matrices. Principal coordinates analysis (PCoA) was applied to the
distance matrices for visualization. Alpha diversity (diversity within sample) was calculated
using Shannon (evenness) indexes. All 454-pyrosequencing files have been deposited to the
National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under
accession number SRP053000.

Statistical analysis

Statistical analyses of HPLC and qPCR data (log10-transformed) were performed using JMP
8.0 (SAS Institute Inc., Cary, NC). Data are expressed as means ± SD of three days at the end
of the stabilization period of each fermentation model. For every model the qPCR and HPLC
data were compared between the reactors using the nonparametric Kruskal-Wallis test.
P values < 0.05 were considered significant. Monte Carlo permutation procedure was used to
determine difference between proximal and distal colon using 999 permutations. Correlation
between genus-level phylotypes and metabolites (acetate, propionate, butyrate, isobutyrate, iso-
valerate and valerate) were done in fermentation models 2 and 3. Analysis was done using R
package "Microbiome" [38] using Spearman correlation. P-values were corrected for multiple
testing using Benjamini-Hochberg. Resulting q values < 0.05 were considered as significant.

Results

Microbial composition of fecal microbiota

The composition of dominant bacterial groups in fecal donor samples was assessed by analyz-
ing the 16S rRNA gene copy numbers of total and selected bacterial groups using qPCR. All
bacterial populations tested were detected in fecal samples, except for the *Roseburia* spp./*E. rec-
tale* group and *Methanobacteriales* that were below the detection limit in the fecal inoculum of
<table><caption>Table 1. qPCR enumeration of bacterial groups in fecal inocula and effluent samples of models' reactors at the end of the stabilization period.</caption><thead><tr><th></th><th>Total 16S rRNA gene</th><th>Bacteroides spp.</th><th>Entero-bacteriaceae</th><th>Lactobacillus spp.</th><th>Bifido-bacterium spp.</th><th>F. prausnitzii</th><th>Clostridium Cluster IV</th><th>Roseburia spp./ E. rectale</th><th>Methano-bacteriales</th></tr></thead><tbody><tr><td colspan="10"><strong>Model 1</strong></td></tr><tr><td>Donor<sup>1</sup></td><td>11.4</td><td>9.3</td><td>7.7</td><td>7.6</td><td>6.6</td><td>8.7</td><td>8.7</td><td>ND</td><td>ND</td></tr><tr><td>PC<sup>2</sup></td><td>10.6 ± 0.02<sup>A</sup></td><td>7.1 ± 0.1<sup>A</sup></td><td>9.4 ± 0.1<sup>A</sup></td><td>7.5 ± 0.2<sup>A</sup></td><td>8.2 ± 0.1<sup>A</sup></td><td>5.4 ± 0.2<sup>A</sup></td><td>6.7 ± 0.4<sup>A</sup></td><td>ND</td><td>ND</td></tr><tr><td>DC<sup>2</sup></td><td>10.5 ± 0.1<sup>A</sup></td><td>9.5 ± 0.01<sup>B</sup></td><td>8.2 ± 0.1<sup>B</sup></td><td>8.0 ± 0.2<sup>B</sup></td><td>8.3 ± 0.1<sup>A</sup></td><td>7.4 ± 0.1<sup>B</sup></td><td>8.1 ± 0.1<sup>B</sup></td><td>ND</td><td>ND</td></tr><tr><td colspan="10"><strong>Model 2</strong></td></tr><tr><td>Donor<sup>1</sup></td><td>11.1</td><td>9.6</td><td>7.2</td><td>8.3</td><td>8.5</td><td>10.1</td><td>10.3</td><td>9.4</td><td>8.6</td></tr><tr><td>IR<sup>2</sup></td><td>10.2 ± 0.2<sup>A</sup></td><td>10.0 ± 0.2<sup>A</sup></td><td>8.9 ± 0.1<sup>A</sup></td><td>6.1 ± 0.1<sup>A</sup></td><td>6.7 ± 0.1<sup>A</sup></td><td>9.1 ± 0.6<sup>A</sup></td><td>10.1 ± 0.2<sup>A</sup></td><td>8.6 ± 0.3<sup>A</sup></td><td>8.8 ± 0.1<sup>A</sup></td></tr><tr><td>PC1<sup>2</sup></td><td>10.3 ± 0.05<sup>A</sup></td><td>10.0 ± 0.2<sup>A</sup></td><td>9.1 ± 0.3<sup>A</sup></td><td>6.1 ± 0.02<sup>A</sup></td><td>7.7 ± 0.1<sup>B</sup></td><td>9.2 ± 0.4<sup>A</sup></td><td>10.1 ± 0.1<sup>A</sup></td><td>8.6 ± 0.3<sup>A</sup></td><td>7.0 ± 0.4<sup>B</sup></td></tr><tr><td>DC1<sup>2</sup></td><td>10.3 ± 0.2<sup>A</sup></td><td>9.9 ± 0.1<sup>A</sup></td><td>8.9 ± 0.1<sup>AB</sup></td><td>6.2 ± 0.1<sup>AB</sup></td><td>7.6 ± 0.1<sup>B</sup></td><td>9.2 ± 0.2<sup>A</sup></td><td>9.9 ± 0.1<sup>A</sup></td><td>8.6 ± 0.2<sup>A</sup></td><td>8.9 ± 0.04<sup>A</sup></td></tr><tr><td>PC2<sup>2</sup></td><td>10.3 ± 0.1<sup>A</sup></td><td>10.0 ± 0.1<sup>A</sup></td><td>8.8 ± 0.3<sup>AB</sup></td><td>6.2 ± 0.1<sup>A</sup></td><td>6.6 ± 0.1<sup>A</sup></td><td>9.0 ± 0.6<sup>A</sup></td><td>10.0 ± 0.2<sup>A</sup></td><td>8.4 ± 0.3<sup>A</sup></td><td>7.3 ± 0.3<sup>B</sup></td></tr><tr><td>DC2<sup>2</sup></td><td>10.1 ± 0.2<sup>A</sup></td><td>9.8 ± 0.2<sup>A</sup></td><td>8.7 ± 0.1<sup>B</sup></td><td>6.3 ± 0.001<sup>B</sup></td><td>6.8 ± 0.2<sup>A</sup></td><td>9.0 ± 0.2<sup>A</sup></td><td>9.9 ± 0.1<sup>A</sup></td><td>8.3 ± 0.2<sup>A</sup></td><td>8.9 ± 0.1<sup>A</sup></td></tr><tr><td colspan="10"><strong>Model 3</strong></td></tr><tr><td>Donor<sup>1</sup></td><td>11.5</td><td>10.6</td><td>6.3</td><td>7.5</td><td>6.2</td><td>10.6</td><td>10.3</td><td>9.7</td><td>7.7</td></tr><tr><td>IR<sup>2</sup></td><td>11.2 ± 0.03<sup>A</sup></td><td>9.9 ± 0.1<sup>A</sup></td><td>9.6 ± 0.1<sup>A</sup></td><td>6.7 ± 0.1<sup>A</sup></td><td>7.8 ± 0.2<sup>AB</sup></td><td>10.2 ± 0.1<sup>A</sup></td><td>9.9 ± 0.1<sup>A</sup></td><td>8.8 ± 0.1<sup>A</sup></td><td>7.2 ± 1.0<sup>A</sup></td></tr><tr><td>CR<sup>2</sup></td><td>11.1 ± 0.1<sup>A</sup></td><td>10.0 ± 0.1<sup>AB</sup></td><td>9.5 ± 0.1<sup>A</sup></td><td>6.9 ± 0.1<sup>AB</sup></td><td>7.6 ± 0.1<sup>AB</sup></td><td>10.1 ± 0.1<sup>A</sup></td><td>9.8 ± 0.1<sup>A</sup></td><td>8.5 ± 0.2<sup>B</sup></td><td>8.3 ± 0.04<sup>A</sup></td></tr><tr><td>TR1<sup>2</sup></td><td>11.2 ± 0.2<sup>A</sup></td><td>10.1 ± 0.03<sup>B</sup></td><td>9.5 ± 0.1<sup>A</sup></td><td>7.0 ± 0.01<sup>B</sup></td><td>7.5 ± 0.1<sup>A</sup></td><td>10.1 ± 0.04<sup>A</sup></td><td>9.8 ± 0.1<sup>A</sup></td><td>8.5 ± 0.2<sup>AB</sup></td><td>8.3 ± 0.04<sup>A</sup></td></tr><tr><td>TR2<sup>2</sup></td><td>11.2 ± 0.05<sup>A</sup></td><td>10.1 ± 0.03<sup>B</sup></td><td>9.5 ± 0.1<sup>A</sup></td><td>7.0 ± 0.2<sup>AB</sup></td><td>7.6 ± 0.2<sup>AB</sup></td><td>10.2 ± 0.02<sup>A</sup></td><td>9.9 ± 0.02<sup>A</sup></td><td>8.6 ± 0.1<sup>AB</sup></td><td>8.3 ± 0.2<sup>A</sup></td></tr><tr><td>TR3<sup>2</sup></td><td>11.2 ± 0.1<sup>A</sup></td><td>10 ± 0.2<sup>AB</sup></td><td>9.5 ± 0.1<sup>A</sup></td><td>7.0 ± 0.2<sup>AB</sup></td><td>7.5 ± 0.2<sup>AB</sup></td><td>10.1 ± 0.1<sup>A</sup></td><td>9.8 ± 0.1<sup>A</sup></td><td>8.5 ± 0.2<sup>B</sup></td><td>8.3 ± 0.2<sup>A</sup></td></tr><tr><td>TR4<sup>2</sup></td><td>11.2 ± 0.1<sup>A</sup></td><td>10.0 ± 0.1<sup>AB</sup></td><td>9.5 ± 0.1<sup>A</sup></td><td>6.8 ± 0.4<sup>AB</sup></td><td>7.7 ± 0.2<sup>B</sup></td><td>10.2 ± 0.1<sup>A</sup></td><td>9.9 ± 0.2<sup>A</sup></td><td>8.6 ± 0.3<sup>AB</sup></td><td>7.9 ± 1.0<sup>A</sup></td></tr></tbody></table>

PC, proximal colon reactor; DC, distal colon reactor; IR, inoculum reactor; CR, control reactor; TR, test reactor; ND, not detected

doi:10.1371/journal.pone.0142793.t001

model 1 (Table 1). Predominant bacterial groups of all three fecal samples were *Bacteroides* spp. and *Clostridium* Cluster IV within which *Faecalibacterium prausnitzii* was dominant. *Enterobacteriaceae*, *Lactobacillus* spp. and *Bifidobacterium* spp. belonged to the subdominant populations in all three fecal inocula.

The microbial profile of fecal inocula of model 2 and 3 was additionally analyzed by pyrosequencing (Fig 2). Phyla of both fecal samples were mainly assigned to Firmicutes and Bacteroidetes, with a Bacteroidetes:Firmicutes ratio of 0.15 and 0.30 for donor 2 and 3, respectively, and followed by Actinobacteria, Proteobacteria and Tenericutes (Fig 2A). The dominant families, *Lachnospiraceae*, *Bacteroidaceae* and *Ruminococcaceae* were similar for both microbiota (Fig 2B). At the genus level an unassigned genus belonging to the family of *Lachnospiraceae* was most abundant in fecal sample 2, followed by *Blautia* and *Bacteroides* (Fig 2C). In fecal sample 3 the same unassigned genus at similar abundance (~ 20%) to fecal sample 2 was most abundant, closely followed by *Bacteroides*.

## Microbial composition of effluents determined by qPCR

Reactor effluents at the end of stabilization of all models were analyzed by qPCR to compare the microbial composition to the fecal donor and between reactors of a model (Table 1). Total bacterial numbers were, in general, high (> 10 log₁₀ gene copies mL⁻<fnref n="1" />) but between 0.3 to 1.0 log₁₀ lower compared to the corresponding donor's fecal sample while no differences were

<footnote n="1">¹Data are mean log₁₀ copies 16S rRNA gene g⁻¹ feces; samples were analyzed in duplicate.</footnote>

<footnote n="2">²Data are mean log₁₀ copies 16S rRNA gene g⁻¹ fermentation effluent ± SD of three last days at the end of the stabilization period; samples were analyzed in duplicate. Values with different letters are significantly different within one model (P < 0.05)</footnote>
A

<figure><img src="image_7.png" /></figure>

B

<figure><img src="image_8.png" /></figure>

C

<figure><img src="image_9.png" /></figure>

**Fig 2.** Microbial composition of fecal samples and reactors of PolyFermS models measured by 454 pyrosequencing. Relative abundance at (A) phylum level of fecal samples of donors 2 and 3, IR and DCI of model 2 and IR and CR of model 3 (B) family level and (C) genus level of fecal samples of donors 2 and 3, all
reactors of model 2 and reactors IR, CR, TR3 and TR4 of model 3 identified by pyrosequencing of the V5-V6 hypervariable regions of the 16S rRNA gene. Effluent samples are average values of three last days at the end of the stabilization period. Parentheses indicate an unknown family belonging to an order or an unknown genus belonging to a family or order. Values < 1% are summarized in the group "others".

doi:10.1371/journal.pone.0142793.g002

observed between reactors within a model. In effluent samples of model 1, the *Roseburia* spp./ *E. rectale* group and *Methanobacteriales* were not detected which was consistent with the lack of these groups in the corresponding fecal donor. *Bacteroides* spp. and *Enterobacteriaceae* were predominant in the PC and DC reactors of model 1, respectively. However copy numbers of specific population groups were significantly different between PC and DC effluents of model 1, with the exception of the total 16S rRNA gene and *Bifidobacterium* spp. The microbial composition of the DC reactor was more similar to the fecal donor than for the PC reactor, except for total bacteria, *Lactobacillus* spp. and *Bifidobacterium* spp.

Similar to the corresponding fecal samples, the predominant bacterial groups in effluents from model 2 and 3 comprised *Bacteroides* spp., *Clostridium* Cluster IV and *Faecalibacterium prausnitzii*. In model 2, no significant difference among all reactors was found for copy numbers of total 16S rRNA gene, *Bacteroides* spp., *Faecalibacterium prausnitzii*, *Clostridium* Cluster IV and *Roseburia* spp./*E. rectale* group. Only small (less than 0.4 log₁₀) but significant differences were detected in DC2 for *Enterobacteriaceae* compared to IR and PC1, and for *Lactobacillus* spp. compared to IR, PC1 and PC2. *Bifidobacterium* spp. gene copy numbers were approx. 1 log higher in test system 1 compared to test system 2 and IR of model 2. *Methanobacteriales* numbers were significantly lower at PC conditions compared to IR and DC reactors.

In model 3, no significant difference was observed between reactors for most tested populations. Only small (≤ 0.3 log₁₀ gene copies mL⁻¹), but significant differences were measured for *Bifidobacterium* spp. numbers in TR4 compared to TR1, for *Bacteroides* and *Lactobacillus* spp. between IR and test reactors or CR, and for *Roseburia* spp. in CR and TR3 compared to IR.

## Microbiota profile and diversity in effluents determined by pyrosequencing

To assess the microbial diversity sequencing of the V5-V6 region of 16S rRNA gene was performed by 454 FLX pyrosequencing of all reactors effluent samples of model 2 and selected reactors (IR, CR, TR3 and TR4) of model 3 (Fig 2) at the end of stabilization phase and compared to diversity of corresponding feces. The main phyla in IR of model 2 and 3 were Firmicutes and Bacteroidetes followed by Actinobacteria, Proteobacteria and Tenericutes (Fig 2A). The ratio Bacteroidetes:Firmicutes was increased in IR and distal colon reactors of model 2 (ratios of 1.7 in IR and 1.2 in DC1) and model 3 (ratios of 0.8 in IR and 1.4 in CR) relative to the corresponding fecal inoculum (ratios 0.2 and 0.3, respectively). Bacteroidetes:Firmicutes ratios in DC2 of model 2 and TR3 and TR4 of model 3 were similar to DC1 and CR, respectively (data not shown). On the family and genus levels *Bacteroidaceae* and *Bacteroides* were dominant in all reactors of both models 2 and 3 (Fig 2B and 2C). *Bacteroidaceae* abundances increased from 10% to approx. 60% and from 36% to approx. 53% in the effluent samples of model 2 and 3, respectively, compared to the fecal donor samples. In general, very similar microbial patterns (family and genus level) were obtained for all reactors within a model. In both DC reactors of model 2 the abundance of *Bacteroidaceae* (54.8% for DC1 and 52.0% for DC2, Fig 2) decreased compared to PC reactors (61.2% for PC1 and 70.1% for PC2). Other small differences at family and genus levels between the PC and DC reactors of each test system were observed. In model 3, minor differences between composition in IR and DC reactors were observed while microbial patterns were highly comparable between the DC reactors.
Weighted Unifrac

A

<figure><img src="image_10.png" /></figure>

Unweighted Unifrac

<figure><img src="image_11.png" /></figure>

B

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

**Fig 3.** PCoA analysis of PolyFermS models based on weighted and unweighted UniFrac analysis. Each symbol is representing a different reactor. (A) Three last days at the end of the stabilization period of all reactors of model 2 (IR, PC1, DC1, PC2 and DC2) and (B) three last days at the end of the stabilization period of model 3 (IR, CR, TR3 and TR4).

doi:10.1371/journal.pone.0142793.g003

Beta diversity (that measures the diversity between samples) of bacterial populations at the end of stabilization phase of model 2 and 3 was analyzed using Principal coordinate analysis (PCoA) (Fig 3). Significant differences between DC and PC were observed using both Unifrac
<figure><img src="image_14.png" /><figcaption>**Fig 4. Shannon diversity index of fecal samples and reactors of PolyFermS models.** The Shannon index was assessed in fecal donors 2 and 3, all reactors of model 2 and reactors IR, CR, TR3 and TR4 of model 3 of three last days at the end of stabilization phase. A higher Shannon index reflects a more diverse community (in abundance and evenness).</figcaption></figure>

doi:10.1371/journal.pone.0142793.g004

distances (p<0.005). In models 2 and 3 a clear separation of reactors operated with proximal colon conditions (IRs, PC1 and PC2) and distal colon conditions (DC1- DC2, CR-TR3-TR4, respectively) was observed.

The Shannon diversity index was assessed for fecal donor samples and effluent samples of models 2 and 3 (Fig 4). A lower diversity was measured in model 2 (mean Shannon index of 5.3 ± 0.4 calculated for all reactors) and model 3 reactors (mean Shannon index of 6.4 ± 0.1 calculated for IR, CR, TR3 and TR4) compared to that of the corresponding fecal samples (Shannon index of 7.5 and 7.4 for model 2 and 3, respectively). In model 2, a higher diversity was obtained for DC (5.6 ± 0.2) compared to PC reactors (4.9 ± 0.2); while in model 3 the Shannon diversity was similar for all tested reactors with values between 6.2 and 6.6.

## Metabolic activity

SCFA were measured by HPLC in fermentation effluents at the end of the stabilization period of all reactors of models 1, 2 and 3 to assess the metabolic activity and intra model stability (Table 2). After the initial stabilization periods, high and stable metabolic activities were
<table><caption>Table 2. Metabolites concentration (mM) and ratios (%) measured by HPLC in effluent samples of models’ reactors at the end of the stabilization period.</caption><thead><tr><th rowspan="2"></th><th colspan="7">Concentrations (mM)</th><th colspan="5">Ratios (%)</th></tr><tr><th>Acetate</th><th>Butyrate</th><th>Propionate</th><th>Valerate</th><th>Isobutyrate</th><th>Isovalerate</th><th>Total SCFA</th><th>Acetate</th><th>Butyrate</th><th>Propionate</th><th>Valerate</th><th>BCFA</th></tr></thead><tbody><tr><td colspan="13"><strong>Model 1</strong></td></tr><tr><td>PC</td><td>74.1 ± 2.4<sup>A</sup></td><td>5.4 ± 1.0<sup>A</sup></td><td>5.4 ± 0.7<sup>A</sup></td><td>ND</td><td>ND</td><td>ND</td><td>84.9 ± 2.7<sup>A</sup></td><td>87.3</td><td>6.4</td><td>6.4</td><td>-</td><td>-</td></tr><tr><td>TC</td><td>108.6 ± 7.1<sup>B</sup></td><td>15.4 ± 2.1<sup>B</sup></td><td>15.3 ± 0.3<sup>B</sup></td><td>ND</td><td>1.9 ± 2.1<sup>A</sup></td><td>4.9 ± 1.2<sup>A</sup></td><td>146.1 ± 7.4<sup>B</sup></td><td>74.3</td><td>10.5</td><td>10.5</td><td>4.7</td><td>4.7</td></tr><tr><td>DC</td><td>110.7 ± 7.8<sup>B</sup></td><td>15.4 ± 0.8<sup>B</sup></td><td>15.3 ± 2.0<sup>B</sup></td><td>ND</td><td>1.6 ± 0.3<sup>A</sup></td><td>4.6 ± 1.5<sup>A</sup></td><td>147.6 ± 8.1<sup>B</sup></td><td>75.0</td><td>10.4</td><td>10.4</td><td>4.2</td><td>4.2</td></tr><tr><td colspan="13"><strong>Model 2</strong></td></tr><tr><td>IR</td><td>66.3 ± 2.7<sup>A</sup></td><td>25 ± 2.9<sup>A</sup></td><td>19.1 ± 0.8<sup>A</sup></td><td>ND</td><td>ND</td><td>ND</td><td>110.4 ± 4.0<sup>A</sup></td><td>60.1</td><td>22.6</td><td>17.3</td><td>-</td><td>-</td></tr><tr><td>PC1</td><td>69.1 ± 2.1<sup>A</sup></td><td>20.8 ± 1.4<sup>AB</sup></td><td>19.9 ± 0.6<sup>A</sup></td><td>ND</td><td>ND</td><td>ND</td><td>109.8 ± 2.6<sup>A</sup></td><td>62.9</td><td>18.9</td><td>18.1</td><td>-</td><td>-</td></tr><tr><td>DC1</td><td>84.7 ± 2.5<sup>B</sup></td><td>20.5 ± 1.3<sup>B</sup></td><td>24.6 ± 0.2<sup>B</sup></td><td>8.2 ± 0.2<sup>A</sup></td><td>6.3 ± 0.2<sup>A</sup></td><td>7.2 ± 0.1<sup>A</sup></td><td>151.5 ± 2.8<sup>B</sup></td><td>55.9</td><td>13.5</td><td>16.2</td><td>4.2</td><td>10.2</td></tr><tr><td>PC2</td><td>70.8 ± 2.9<sup>A</sup></td><td>20.2 ± 3.3<sup>AB</sup></td><td>26.5 ± 3.3<sup>BC</sup></td><td>ND</td><td>ND</td><td>ND</td><td>117.5 ± 5.5<sup>C</sup></td><td>60.3</td><td>17.2</td><td>22.6</td><td>-</td><td>-</td></tr><tr><td>DC2</td><td>81.6 ± 1.1<sup>B</sup></td><td>19.7 ± 2.6<sup>B</sup></td><td>26.5 ± 1.8<sup>C</sup></td><td>5.7 ± 0.2<sup>B</sup></td><td>6.4 ± 0.2<sup>A</sup></td><td>6.9 ± 0.4<sup>B</sup></td><td>146.8 ± 3.3<sup>D</sup></td><td>55.6</td><td>13.4</td><td>18.1</td><td>4.4</td><td>8.6</td></tr><tr><td colspan="13"><strong>Model 3</strong></td></tr><tr><td>IR</td><td>56.5 ± 5.0<sup>A</sup></td><td>40.2 ± 5.6<sup>A</sup></td><td>23.7 ± 3.7<sup>A</sup></td><td>ND</td><td>ND</td><td>ND</td><td>120.4 ± 8.4<sup>A</sup></td><td>46.9</td><td>33.4</td><td>19.7</td><td>-</td><td>-</td></tr><tr><td>CR</td><td>80.7 ± 4.5<sup>B</sup></td><td>42.8 ± 1.8<sup>A</sup></td><td>31.5 ± 2.4<sup>BC</sup></td><td>ND</td><td>ND</td><td>0.8 ± 0.1<sup>A</sup></td><td>155.8 ± 5.4<sup>B</sup></td><td>51.8</td><td>27.5</td><td>20.2</td><td>-</td><td>0.5</td></tr><tr><td>TR1</td><td>76.8 ± 4.1<sup>BC</sup></td><td>42.7 ± 5.0<sup>A</sup></td><td>31.0 ± 1.7<sup>BC</sup></td><td>ND</td><td>ND</td><td>ND</td><td>150.5 ± 6.7<sup>BC</sup></td><td>51.0</td><td>28.4</td><td>20.6</td><td>-</td><td>-</td></tr><tr><td>TR2</td><td>73.2 ± 1.9<sup>C</sup></td><td>40.9 ± 1.3<sup>A</sup></td><td>29.0 ± 0.9<sup>B</sup></td><td>ND</td><td>ND</td><td>0.4 ± 0.1<sup>A</sup></td><td>143.5 ± 2.5<sup>C</sup></td><td>51.0</td><td>28.5</td><td>20.2</td><td>-</td><td>0.3</td></tr><tr><td>TR3</td><td>78.1 ± 1.1<sup>B</sup></td><td>39.7 ± 1.4<sup>A</sup></td><td>30.9 ± 0.9<sup>C</sup></td><td>ND</td><td>ND</td><td>1.0 ± 0.8<sup>A</sup></td><td>149.7 ± 2.0<sup>B</sup></td><td>52.2</td><td>26.5</td><td>20.6</td><td>-</td><td>0.7</td></tr><tr><td>TR4</td><td>73.7 ± 0.7<sup>C</sup></td><td>42.4 ± 1.9<sup>A</sup></td><td>28.8 ± 1.1<sup>B</sup></td><td>ND</td><td>ND</td><td>ND</td><td>144.9 ± 2.3<sup>C</sup></td><td>50.9</td><td>29.3</td><td>19.9</td><td>-</td><td>-</td></tr></tbody></table>

PC, proximal colon reactor; DC, distal colon reactor; IR, inoculum reactor; CR, control reactor; TR, test reactor

Data are means ± SD of three last days at the end of the stabilization period; samples were analyzed in duplicate. ND, not detected

Values with different letters are significantly different within one model (P < 0.05)

doi:10.1371/journal.pone.0142793.t002

measured over the entire fermentation in IR's of models 2 and 3 which were operated under constant conditions and used to demonstrate temporal stability of the PolyFermS models (Fig 5).

The concentrations of SCFA tested in reactors were model and reactor (-proximal and distal colon conditions) dependent, while intermediate products lactate and formate remained undetected in fermentation effluents. Acetate was the main metabolite in reactor effluents of all models, followed by butyrate and propionate which were generally produced at similar levels within a reactor. Butyrate concentrations were around 10 mM higher than propionate in CR and TR reactors of model 3. The molar ratios of acetate, butyrate and propionate in IR of the three models were different. In IR of model 1, operated at pH 5.5, a higher acetate fraction was produced (87/6/6) compared to IR of model 2 (60/23/17) and 3 (47/33/20) which were operated at a higher pH of 5.7. Higher concentrations of acetate and propionate were measured in TC and DC reactors of model 1, and DC reactors of models 2 and 3 compared to reactors IR and PCs of the same models, operated with proximal colon conditions. Butyrate concentrations increased along the reactors of the 3-stage model 1, but remained unchanged between PC and DC reactors of models 2 and 3. The mean concentrations of acetate, butyrate and propionate in the 5 distal reactors of model 3, were 76.5 ± 6.5, 41.7 ± 6.0 and 30.2 ± 3.4 mM, respectively, with small (less than 8 mM) but significant differences among reactors for acetate and propionate. Valerate was only detected in DC reactors from model 2.

Branched-chain fatty acids (BCFA) were not detected in any PC reactors. Isovalerate and isobutyrate were present in effluents samples from DC reactors of model 1 and model 2, at higher concentration in the latter. For model 3 isovalerate was only measured at low concentrations (≤ 1.0 mM) close to the detection limit in some DC reactors.
<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /><figcaption>**Fig 5.** Daily mean SCFA concentrations in fermentation effluents of IR of PolyFermS models measured by HPLC. Initial stabilization: stabilization period in continuous mode to reach pseudo steady-state.; stabilized operation mode: continuous operation mode during pseudo steady-state conditions. (A) Model 2 and (B) model 3; (♦) acetate, (●) butyrate, and (■) propionate.</figcaption></figure>

doi:10.1371/journal.pone.0142793.g005

## Correlations between microbiota composition and metabolite production

Pyrosequencing data on the genus level and metabolite concentrations measured by HPLC were investigated to test significant correlations between phylogenetic groups and metabolic activity. For model 2, significant negative correlations were calculated between isobutyrate, isovalerate and valerate concentrations and the dominant genera *Ruminococcus* and *Bacteroides* (Fig 6). Butyrate was positively correlated with the dominant genus *Roseburia* and unclassified members of *Ruminococcaceae* and *Lachnospiraceae*. In contrast, a dominant unclassified member of *Enterococcaceae* was negatively correlated with butyrate and positively correlated with all other metabolites detected. Furthermore many genera detected at less than 1% (*Dialister*, *Anaerococcus* and unclassified genera of *Rikenellaceae* and *Mogibacteriaceae*) showed positive correlations with isobutyrate, isovalerate and valerate with the exception of *Oscillospira*, *Peptostreptococcaceae* (with abundances above 1% but in distal colon reactors only). Correlations between metabolites (only acetate and propionate) and phylogenetic groups were also found for model 3 (S1 Fig).
<figure><img src="image_17.png" /><figcaption>**Fig 6.** Correlations between genus-level phylogenetic groups and metabolites (SCFA, BCFA) of three last days at the end of stabilization period of model 2. The correlations, assessed by Spearman are indicated by either red (positive) and blue (negative), the significant correlations (q < 0.05) are indicated by '+'. Only genus related phylotypes > 0.1% and with at least one significant correlation with metabolites are depicted. Parentheses indicate an unknown genus belonging to a family or order.</figcaption></figure>

doi:10.1371/journal.pone.0142793.g006

## Discussion

Colonic fermentation models are useful tools to investigate factors that can influence the composition and metabolism of the gut microbiota, such as diet, antibiotic treatment, and bacterial infections *in vitro* and independent of the host [9, 10]. An important aspect for *in vitro* studies is the rational design of models and conditions, considering host target, model characteristics and limits, and the recognition that models are not perfect representation of reality. Therefore differences are often observed between fermentation samples and donor's feces. A major discriminatory factor of *in vitro* models is the technique used for fecal inoculation. For most models, a fresh liquid fecal suspension is inoculated whereas this has been shown to lead to limited stability (washout of less competitive or slow-growing bacteria), cell density and difficulty to reproduce both the planktonic (lumen) and sessile (food particule and mucus associated) microbiota of the colon.

In the present study, we report the first-time investigation of continuous fermentation models with fecal microbiota obtained from different healthy volunteers aged between 71 and 78 years using in-depth characterization methods of the microbial diversity. We immobilized the fecal microbiota and inoculated the fecal biocatalysts in the inoculum reactor of the tested
models with different designs. A dense and diverse microbiota could be established in Poly-
FermS models, with reproducible microbial composition and metabolic activity for down-
stream test and control reactors within a model.

During collection and immobilization of fecal microbiota from each elderly donor special attention was paid to keep anaerobic conditions from donor to reactor, in order to reproduce both the planktonic and sessile forms of bacteria in the colon, as previously suggested [10]. Gel beads can provide a protective microenvironment for the bacteria and allow the growth of complex and stable gut ecosystems at high cell densities of up to ca. log 11 cells per mL effluent as observed in the present study with elder gut microbiota, preventing the loss of slow growing bacteria. As expected from the lack of water reabsorption total bacteria numbers in reactor effluents of models were up to 1.0 log₁₀ lower compared to the corresponding donor's fecal sample. All bacterial groups tested in the fecal inoculum with qPCR were present in the corresponding models. The main differences in the bacterial composition and metabolic activity amongst models can be assigned to the different fecal inoculum used. In particular, *Roseburia* spp. was not detected in the feces and effluents of model 1 and this may explain the high acetate and low butyrate concentrations in this model (Table 2) since *Roseburia* is a main contributor for the conversion of acetate into butyrate [39, 40].

The pH of IR of model 1 was set to 5.5 in order to replicate the pH set in previous elder gut
fermentation models [15, 17]. However, this pH is in the low range for the human proximal
colon *in vivo*, [29, 41] and for models 2 and 3, the pH in IR and PCs was set to 5.7. This pH ele-
vation induced an increase in total metabolites by approximately 20%, in agreement with previ-
ous observations of pH effect made in the PolyFermS model with child microbiota [13]. In
contrast, butyrate concentrations did not increase in model 1 (pH 5.5) relative to model 2 and
3 (pH 5.7), as would have been expected from the stimulation of butyrate production at the
lower pH, as observed in the previous study [13]. This is likely due to the lack of *Roseburia* spp.
and the low *F. prausnitzii* numbers in model 1 which are the main butyrate producers in the
human gut microbiota [42]. In models 2 and 3, the microbiota composition tested with qPCR
was very similar between IR or PC reactors and distal reactors, while some limited changes
were measured with pyrosequencing. In contrast, most targeted populations by qPCR signifi-
cantly increased from PC (IR) to DC of model 1, suggesting that the low pH of 5.5 limited the
growth of the targeted groups.

With qPCR we detected high *Enterobacteriaceae* copy numbers in reactor effluents com-
pared to feces for all three models. This was observed in previous gut fermentation models [14,
21] and may be due to competitive advantage of these fast growing and robust bacteria that
allows them to occupy niches during the immobilization process and the succeeding batch fer-
mentation. The low levels of SCFA in the beginning of batch fermentation may further explain
the increase in *Enterobacteriaceae* in reactors, as SCFA have inhibitory effects against *Entero-
bacteriaceae*, such as *Escherichia coli* [43].

The microbial composition of models 2 and 3 and corresponding fecal inocula may be con-
sidered more representative of the elderly population than the fecal inoculum of model 1,
which did not harbor *Roseburia* spp., although the genus *Roseburia* was assigned at approxi-
mately 3% in fecal samples from elderly Irish subjects [8].

454 pyrosequencing was performed using the V5-V6 hypervariable region that was previ-
ously used to profile gut microbiota [44–46]. Sun *et al.* [47] recently reported that intragenomic
heterogeneity for the V6 region may introduce overestimation of prokaryotes diversity. How-
ever in our study pyrosequencing data were used to compare of composition of donor and
reactor samples within a model which should not be affected by this possible bias. In general,
similar microbial profiles between effluent samples of model 2 and 3 and corresponding fecal
donors were obtained. However, in both models the ratio of Firmicutes:Bacteroidetes was
decreased when compared to the fecal sample, likely due to host-related factors including water and metabolite absorption and intestinal cells and host interaction, both of which are lacking in the fermentation models [13, 48]. Changes in microbiota composition and diversity may also reflect a possible loss of bacteria from donor to reactor, even though great care was taken to protect viability (strict anaerobiosis, mild conditions, short time), and the use of fecal microbiota to inoculate the inoculum reactor run in proximal colon conditions. Furthermore, the strictly controlled environmental factors, such as pH, transit time and medium composition in the *in vitro* models do no fully represent the specific donor conditions, thereby further contributing to *in vitro* and *in vivo* variations [10]. Despite the increase in Bacteroidetes in models 2 and 3, the Firmicutes:Bacteroidetes ratios of the fecal donor and the models were all in the range of previously reported data recorded during a large-scale *in vivo* study with elderly Irish people [8]. Indeed large inter-individual variations were observed in this study; however, the Bacteroidetes:Firmicutes ratio was shown to be higher in the elder, relative to the adult population. The majority of the reads was assigned to Firmicutes and Bacteroidetes while only low levels of Proteobacteria and Actinobacteria were detected in concordance with *in vivo* findings [8, 49]. Many of the predominant genera (abundance > 1%) including *Bacteroides*, *Faecalibacterium*, *Roseburia* and *Ruminococcus* were also found above 1% in elderly Irish subjects [8]. No in-depth characterization of the microbiota in donor and effluent samples was reported in previous investigations of fermentation models of the elderly microbiota, in which only traditional plating methods [15, 18] or FISH [19] were used.

SCFA are mainly produced from carbohydrate fermentation and to a lesser extent via degradation of proteins and amino acids; the effects of SCFA on the host are well documented [1]. It was previously found that the major SCFA found in stools of healthy volunteers between 14–74 years of age were: acetate, propionate and butyrate at an approx. ratio of 3:1:1 [50]. In our study, similar ratios were found in distal colon reactors of models 2 (4:1.5:1) and 3 (3:1.5:2), whereas in model 1 the acetate fraction was considerably higher (7:1:1), likely due to the low pH of 5.5 in IR and the fecal microbiota used in this model as discussed above. *In vivo* investigations are, however, hampered by the continuous absorption of metabolic products, which results in less than 5% of total production excreted in the feces [51] along with the difficulty associated with obtaining samples from different regions of the colon. Therefore, metabolite concentrations and ratios in feces are not indicative for the colonic microbiota activity. In contrast *in vitro* modeling allows accurate measuring of the metabolic activity of the gut microbiota for the tested model conditions. Stable SCFA concentrations were obtained throughout the fermentation in the untreated IR's of the PolyFermS models demonstrating maintenance of gut microbial activity over the entire fermentation of 55 and 80 days for model 2 and 3, respectively (Fig 5).

BCFA are products of protein and amino acid fermentation but the formation of BCFA and associated species is not well studied [52]. In the colon of elders, an increase in proteolytic activity and a decrease in concentrations of SCFA were reported [53, 54]. Metabolites of amino acid fermentation can have toxic effects on the colonic lumen and were associated with several gut disorders [52]. In the tested models, BCFA were solely detected in significant levels within the distal reactors of model 1 and 2. This observation is consistent with the understanding that the distal colon is the major site for proteolysis whereas carbohydrate fermentation is the main energy yielding process in the proximal colon, resulting in a lower pH in this section [43, 55]. In model 2, genera with abundances of less than 1% were positively correlated with isobutyrate, valerate and isovalerate, suggesting that the dominant bacteria were mainly responsible for saccharolytic fermentation while proteins were degraded by the subdominant populations. Very low or no BCFA were detected in CR and TR reactors of model 3 which was set to mimic fermentation of transverse-distal colon sections within one reactor. This may be explained by the
microbiota composition of model 3 that was different from model 1 and 2. The microbiota-
dependent production of BCFA suggests the importance of using individual microbiota for
inoculating intestinal fermentation models instead of pooling microbiota from different
donors, as done in many studies for inoculation of gut fermentation models.

A major feature of the PolyFermS models over the three-stage model (model 1) is that sev-
eral treatments can be investigated simultaneously and compared to a control inoculated with
the same microbiota, thus generating reproducible and accurate data rather than when treat-
ments are applied during consecutive periods. In both PolyFermS models microbial diversity
and metabolic activity was very similar between control and test reactors. Model 2 built with
two sets of proximal and distal colon reactors can be used for a broad range of studies, in proxi-
mal and distal colon conditions, such as the effect of an altered diet and administration of anti-
biotics on the gut microbiota in old age. Furthermore, the model is applicable for the *in vitro*
investigation of the elderly microbiota in combination with health-related questions such as
the manipulation of the gut microbiota using pro- and prebiotics [56]. PolyFermS model 3
built with multiple parallel distal colon reactors can be especially useful to study the effect of
factors related to age on microbial metabolism in the lower colon, such as promotion of putre-
faction due to low fiber intake. The PolyFermS intestinal platform has potential to be scaled
down and adapted with multi-reactors to enhance screening efficiency.

To conclude, in the present study we showed the stability and reproducibility of PolyFermS
continuous colonic fermentation models inoculated with immobilized elderly microbiota.
Immobilization requires only small amounts of high quality fecal material to prime a gut
model that can be stably operated over several months for testing parallel treatments in conse-
cutive blocks [56]. The PolyFermS platform should be suitable for a range of *in vitro* gut micro-
biota investigations, from classical microbe interaction studies to complex ecological studies of
the elderly gut microbiome investigated by in-depth analysis of the microbial diversity.

**Supporting Information**

S1 Fig. Correlations between genus-level phylogenetic groups and acetate and propionate of three last days at the end of the stabilization period of model 3. The correlations, assessed by Spearman are indicated by either red (positive) and blue (negative), the significant correlations (q < 0.05) are indicated by '+'. Only genus related phylotypes > 0.1% and with at least one significant correlation with metabolites are depicted.
(TIF)

S1 Table. Primers used for enumeration of bacterial groups by qPCR.
(DOCX)

**Acknowledgments**

We thank Liselotte Wallquist, Sophie Poeker, Lukas Meile, Simon Galenda and Markus Reich-
lin for technical assistance.

**Author Contributions**

Conceived and designed the experiments: SF CC MCH CL. Performed the experiments: SF MCH CC. Analyzed the data: SF CC MCH MD CL. Contributed reagents/materials/analysis tools: CL CF. Wrote the paper: SF CC MD CL.
## References

1. Chassard C, Lacroix C. Carbohydrates and the human gut microbiota. Curr Opin Clin Nutr Metab Care. 2013; 16(4):453-60. Epub 2013/05/31. doi: 10.1097/MCO.0b013e3283619e63 PMID: 23719143.

2. Robles Alonso V, Guarner F. Linking the gut microbiota to human health. Brit J Nutr. 2013; 109 Suppl 2: S21-6. Epub 2013/02/13. doi: 10.1017/S0007114512005235 S0007114512005235 [pii]. PMID: 23360877.

3. Duncan SH, Flint HJ. Probiotics and prebiotics and health in ageing populations. Maturitas. 2013; 75 (1):44-50. Epub 2013/03/16. S0378-5122(13)00049-2 [pii] doi: 10.1016/j.maturitas.2013.02.004 PMID: 23489554.

4. Power SE, O'Toole PW, Stanton C, Ross RP, Fitzgerald GF. Intestinal microbiota, diet and health. Brit J Nutr. 2014; 111(3):387-402. Epub 2013/08/13. doi: 10.1017/S0007114513002560 S0007114513002560 [pii]. PMID: 23931069.

5. O'Toole PW, Claesson MJ. Gut microbiota: Changes throughout the lifespan from infancy to elderly. Int Dairy J. 2010; 20(4):281-91. doi: 10.1016/j.idairyj.2009.11.010 PMID: ISI:000275924800012.

6. Biagi E, Candela M, Turroni S, Garagnani P, Franceschi C, Brigidi P. Ageing and gut microbes: Perspectives for health maintenance and longevity. Pharmacol Res. 2013; 69(1):11-20. Epub 2012/10/20. doi: 10.1016/j.phrs.2012.10.005 S1043-6618(12)00192-2 [pii]. PMID: 23079287.

7. Ottman N, Smidt H, de Vos WM, Belzer C. The function of our microbiota: who is out there and what do they do? Front Cell Infect Mi. 2012; 2:104. Unsp 104 doi: 10.3389/Fcimb.2012.00104 PMID: WOS:000324630500002; PubMed Central PMCID: PMC3417542.

8. Claesson MJ, Cusack S, O'Sullivan O, Greene-Diniz R, de Weerd H, Flannery E, et al. Composition, variability, and temporal stability of the intestinal microbiota of the elderly. Proc Natl Acad Sci U S A. 2011; 108 Suppl 1:4586-91. Epub 2010/06/24. doi: 10.1073/pnas.1000097107 1000097107 [pii]. PMID: 20571116; PubMed Central PMCID: PMC3063589.

9. Lacroix C, de Wouters T, Chassard C. Integrated multi-scale strategies to investigate nutritional compounds and their effect on the gut microbiota. Curr Opin Biotechnol. 2015; 32C:149-55. doi: 10.1016/j.copbio.2014.12.009 PMID: 25562815.

10. Payne AN, Zihler A, Chassard C, Lacroix C. Advances and perspectives in *in vitro* human gut fermentation modeling. Trends Biotechnol. 2012; 30(1):17-25. Epub 2011/07/19. doi: 10.1016/j.tibtech.2011.06.011 S0167-7799(11)00115-6 [pii]. PMID: 21764163.

11. Cinquin C, Le Blay G, Fliss I, Lacroix C. Immobilization of infant fecal microbiota and utilization in an *in vitro* colonic fermentation model. Microb Ecol. 2004; 48(1):128-38. Epub 2004/04/16. doi: 10.1007/s00248-003-2022-7 PMID: 15085302.

12. Cinquin C, Le Blay G, Fliss I, Lacroix C. New three-stage *in vitro* model for infant colonic fermentation with immobilized fecal microbiota. FEMS Microbiol Ecol. 2006; 57(2):324-36. Epub 2006/07/27. FEM117 [pii] doi: 10.1111/j.1574-6941.2006.00117.x PMID: 16867149.

13. Zihler Berner A, Fuentes S, Dostal A, Payne AN, Vazquez Gutierrez P, Chassard C, et al. Novel Poly-fermentor intestinal model (PolyFermS) for controlled ecological studies: validation and effect of pH. PLoS One. 2013; 8(10):e77772. Epub 2013/11/10. doi: 10.1371/journal.pone.0077772 PONE-D-13-16596 [pii]. PMID: 24204958; PubMed Central PMCID: PMC3813750.

14. Tanner SA, Zihler Berner A, Rigozzi E, Grattepanche F, Chassard C, Lacroix C. *In vitro* continuous fermentation model (PolyFermS) of the swine proximal colon for simultaneous testing on the same gut microbiota. PLoS One. 2014; 9(4):e94123. doi: 10.1371/journal.pone.0094123 PMID: 24709947; PubMed Central PMCID: PMC3978012.

15. Baines SD, Freeman J, Wilcox MH. Effects of piperacillin/tazobactam on *Clostridium difficile* growth and toxin production in a human gut model. J Antimicrob Chemother. 2005; 55(6):974-82. Epub 2005/04/30. doi:10.1093/jac/dki120 [pii] doi: 10.1093/jac/dki120 PMID: 15860551.

16. Baines SD, O'Connor R, Saxton K, Freeman J, Wilcox MH. Activity of vancomycin against epidemic *Clostridium difficile* strains in a human gut model. J Antimicrob Chemother. 2009; 63(3):520-5. Epub 2008/12/30. doi:10.1093/jac/dkn502 [pii] doi: 10.1093/jac/dkn502 PMID: 19112083.

17. Freeman J, Baines SD, Saxton K, Wilcox MH. Effect of metronidazole on growth and toxin production by epidemic *Clostridium difficile* PCR ribotypes 001 and 027 in a human gut model. J Antimicrob Chemother. 2007; 60(1):83-91. Epub 2007/05/08. doi:10.1093/jac/dkm113 [pii] doi: 10.1093/jac/dkm113 PMID: 17483547.

18. Crowther GS, Chilton CH, Todhunter SL, Nicholson S, Freeman J, Baines SD, et al. Development and validation of a chemostat gut model to study both planktonic and biofilm modes of growth of *Clostridium difficile* and human microbiota. PLoS One. 2014; 9(2):e88396. Epub 2014/02/12. doi: 10.1371/journal.pone.0088396 PONE-D-13-46790 [pii]. PMID: 24516647.
19. Likotrafiti E, Tuohy KM, Gibson GR, Rastall RA. An *in vitro* study of the effect of probiotics, prebiotics and synbiotics on the elderly faecal microbiota. Anaerobe. 2014; 27:50–5. doi: 10.1016/j.anaerobe. 2014.03.009 PMID: 24685554.

20. Payne AN, Chassard C, Banz Y, Lacroix C. The composition and metabolic activity of child gut micro- biota demonstrate differential adaptation to varied nutrient loads in an *in vitro* model of colonic fermen- tation. FEMS Microbiol Ecol. 2012; 80(3):608–23. Epub 2012/02/14. doi: 10.1111/j.1574-6941.2012. 01330.x PMID: 22324938.

21. Zihler A, Gagnon M, Chassard C, Hegland A, Stevens MJ, Braegger CP, et al. Unexpected conse- quences of administering bacteriocinogenic probiotic strains for *Salmonella* populations, revealed by an *in vitro* colonic model of the child gut. Microbiology. 2010; 156(Pt 11):3342–53. Epub 2010/08/07. mic.0.042036–0 [pii] doi: 10.1099/mic.0.042036–0 PMID: 20688827.

22. Macfarlane GT, Macfarlane S, Gibson GR. Validation of a three-stage compound continuous culture system for investigating the effect of retention time on the ecology and metabolism of bacteria in the human colon. Microb Ecol. 1998; 35(2):180–7. Epub 1998/12/16. MECOJKF97-5 [pii]. PMID: 9541554.

23. Michel C, Kravtchenko TP, David A, Gueneau S, Kozlowski F, Cherbut C. *In vitro* prebiotic effects of Acacia gums onto the human intestinal microbiota depends on both botanical origin and environmental pH. Anaerobe. 1998; 4(6):257–66. doi: 10.1006/anae.1998.0178 PMID: 16887651.

24. Nagengast FM, van der Werf SD, Lamers HL, Hectors MP, Buys WC, van Tongeren JM. Influence of age, intestinal transit time, and dietary composition on fecal bile acid profiles in healthy subjects. Digest Dis Sci. 1988; 33(6):673–8. PMID: 3371139.

25. Madsen JL. Effects of gender, age, and body mass index on gastrointestinal transit times. Digest Dis Sci. 1992; 37(10):1548–53. Epub 1992/10/01. PMID: 1396002.

26. Madsen JL, Graff J. Effects of ageing on gastrointestinal motor function. Age Ageing. 2004; 33(2):154–9. doi: 10.1093/ageing/afh040 PMID: 14960431.

27. Merkel IS, Locher J, Burgio K, Towers A, Wald A. Physiologic and psychologic characteristics of an elderly population with chronic constipation. Am J Gastroenterol. 1993; 88(11):1854–9. Epub 1993/11/ 01. PMID: 8237932.

28. Evans JM, Fleming KC, Talley NJ, Schleck CD, Zinsmeister AR, Melton LJ 3rd. Relation of colonic tran- sit to functional bowel disease in older people: a population-based study. J Am Geriatr Soc. 1998; 46 (1):83–7. Epub 1998/01/22. PMID: 9434670.

29. Lawley TD, Walker AW. Intestinal colonization resistance. Immunology. 2013; 138(1):1–11. doi: 10. 1111/j.1365-2567.2012.03616.x PMID: 23240815; PubMed Central PMCID: PMC3533696.

30. Dostal A, Fehlbaum S, Chassard C, Zimmermann MB, Lacroix C. Low iron availability in continuous *in vitro* colonic fermentations induces strong dysbiosis of the child gut microbial consortium and a decrease in main metabolites. FEMS Microbiol Ecol. 2013; 83(1):161–75. Epub 2012/08/01. doi: 10. 1111/j.1574-6941.2012.01461.x PMID: 22845175; PubMed Central PMCID: PMC3511601.

31. Cleusix V, Lacroix C, Vollenweider S, Le Blay G. Glycerol induces reuterin production and decreases *Escherichia coli* population in an *in vitro* model of colonic fermentation with immobilized human feces. FEMS Microbiol Ecol. 2008; 63(1):56–64. Epub 2007/11/22. FEM412 [pii] doi: 10.1111/j.1574-6941. 2007.00412.x PMID: 18028400.

32. Andersson AF, Lindberg M, Jakobsson H, Backhed F, Nyren P, Engstrand L. Comparative analysis of human gut microbiota by barcoded pyrosequencing. PLoS One. 2008; 3(7):e2836. Epub 2008/07/31. doi: 10.1371/journal.pone.0002836 PMID: 18665274; PubMed Central PMCID: PMC2475661.

33. 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(5):335–6. doi: 10.1038/ nmeth.f.303 PMID: 20383131; PubMed Central PMCID: PMC3156573.

34. Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006; 22(13):1658–9. doi: 10.1093/bioinformatics/btl158 PMID: 16731699.

35. Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, et al. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 2011; 21(3):494–504. doi: 10.1101/gr.112730.110 PMID: 21212162; PubMed Central PMCID: PMC3044863.

36. Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microb. 2005; 71(12):8228–35. doi: 10.1128/AEM.71.12.8228–8235.2005 PMID: 16332807; PubMed Central PMCID: PMC1317376.

37. Lozupone CA, Hamady M, Kelley ST, Knight R. Quantitative and qualitative beta diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microb. 2007; 73(5):1576–85. doi: 10.1128/AEM.01996-06 PMID: 17220268; PubMed Central PMCID: PMC1828774.
38. Lahti L, Salonen A, Kekkonen RA, Salojarvi J, Jalanka-Tuovinen J, Palva A, et al. Associations between the human intestinal microbiota, *Lactobacillus rhamnosus* GG and serum lipids indicated by integrated analysis of high-throughput profiling data. PeerJ. 2013; 1:e32. doi: 10.7717/peerj.32 PMID: 23638368; PubMed Central PMCID: PMC3628737.

39. Duncan SH, Barcenilla A, Stewart CS, Pryde SE, Flint HJ. Acetate utilization and butyryl coenzyme A (CoA):acetate-CoA transferase in butyrate-producing bacteria from the human large intestine. Appl Environ Microb. 2002; 68(10):5186-90. PMID: 12324374; PubMed Central PMCID: PMC126392.

40. Belenguer A, Duncan SH, Calder AG, Holtrop G, Louis P, Lobley GE, et al. Two routes of metabolic cross-feeding between *Bifidobacterium adolescentis* and butyrate-producing anaerobes from the human gut. Appl Environ Microb. 2006; 72(5):3593-9. doi: 10.1128/AEM.72.5.3593-3599.2006 PMID: 16672507; PubMed Central PMCID: PMC1472403.

41. Nugent SG, Kumar D, Rampton DS, Evans DF. Intestinal luminal pH in inflammatory bowel disease: possible determinants and implications for therapy with aminosalicylates and other drugs. Gut. 2001; 48(4):571-7. PMID: 11247905; PubMed Central PMCID: PMC1728243.

42. Walker AW, Duncan SH, McWilliam Leitch EC, Child MW, Flint HJ. pH and peptide supply can radically alter bacterial populations and short-chain fatty acid ratios within microbial communities from the human colon. Appl Environ Microb. 2005; 71(7):3692-700. Epub 2005/07/08. 71/7/3692 [pii] doi: 10.1128/AEM.71.7.3692-3700.2005 PMID: 16000778; PubMed Central PMCID: PMC1169066.

43. Duncan SH, Louis P, Thomson JM, Flint HJ. The role of pH in determining the species composition of the human colonic microbiota. Environ Microbiol. 2009; 11(8):2112-22. Epub 2009/04/29. EMI1931 [pii] doi: 10.1111/j.1462-2920.2009.01931.x PMID: 19397676.

44. Zhang J, Guo Z, Xue Z, Sun Z, Zhang M, Wang L, et al. A phylo-functional core of gut microbiota in healthy young Chinese cohorts across lifestyles, geography and ethnicities. ISME J. 2015; 9(9):1979-90. doi: 10.1038/ismej.2015.11 PMID: 25647347.

45. Montassier E, Batard E, Massart S, Gastinne T, Carton T, Caillon J, et al. 16S rRNA gene pyrosequencing reveals shift in patient faecal microbiota during high-dose chemotherapy as conditioning regimen for bone marrow transplantation. Microb Ecol. 2014; 67(3):690-9. doi: 10.1007/s00248-013-0355-4 PMID: 24402367.

46. De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, et al. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc Natl Acad Sci U S A. 2010; 107(33):14691-6. doi: 10.1073/pnas.1005963107 PMID: 20679230; PubMed Central PMCID: PMC2930426.

47. Sun DL, Jiang X, Wu QL, Zhou NY. Intragenomic heterogeneity of 16S rRNA genes causes overestimation of prokaryotic diversity. Appl Environ Microb. 2013; 79(19):5962-9. doi: 10.1128/AEM.01282-13 PMID: 23872556; PubMed Central PMCID: PMC3811346.

48. Van den Abbeele P, Grootaert C, Marzorati M, Possemiers S, Verstraete W, Gerard P, et al. Microbial community development in a dynamic gut model is reproducible, colon region specific, and selective for Bacteroidetes and *Clostridium* cluster IX. Appl Environ Microb. 2010; 76(15):5237-46. doi: 10.1128/AEM.00759-10 PMID: 20562281; PubMed Central PMCID: PMC2916472.

49. Biagi E, Nylund L, Candela M, Ostan R, Bucci L, Pini E, et al. Through ageing, and beyond: gut microbiota and inflammatory status in seniors and centenarians. PLoS One. 2010; 5(5):e10667. Epub 2010/05/26. doi: 10.1371/journal.pone.0010667 PMID: 20498852; PubMed Central PMCID: PMC2871786.

50. Schwiertz A, Taras D, Schafer K, Beijer S, Bos NA, Donus C, et al. Microbiota and SCFA in lean and overweight healthy subjects. Obesity. 2010; 18(1):190-5. doi: 10.1038/oby.2009.167 PMID: 19498350.

51. den Besten G, van Eunen K, Groen AK, Venema K, Reijingoud DJ, Bakker BM. The role of short-chain fatty acids in the interplay between diet, gut microbiota, and host energy metabolism. J Lipid Res. 2013; 54(9):2325-40. doi: 10.1194/jlr.R036012 PMID: 23821742; PubMed Central PMCID: PMC3735932.

52. Nyangale EP, Mottram DS, Gibson GR. Gut microbial activity, implications for health and disease: the potential role of metabolite analysis. J Proteome Res. 2012; 11(12):5573-85. Epub 2012/11/03. doi: 10.1021/pr300637d PMID: 23116228.

53. Woodmansey EJ. Intestinal bacteria and ageing. J Appl Microbiol. 2007; 102(5):1178-86. Epub 2007/04/24. JAM3400 [pii] doi: 10.1111/j.1365-2672.2007.03400.x PMID: 17448153.

54. Tiihonen K, Tynkkynen S, Ouwehand A, Ahlroos T, Rautonen N. The effect of ageing with and without non-steroidal anti-inflammatory drugs on gastrointestinal microbiology and immunology. Brit J Nutr. 2008; 100(1):130-7. Epub 2008/02/19. doi: 10.1017/S000711450888871X S000711450888871X [pii]. PMID: 18279548.

55. Cummings JH, Macfarlane GT. The control and consequences of bacterial fermentation in the human colon. J Appl Bacteriol. 1991; 70(6):443-59. PMID: 1938669.
56. Tanner SA, Chassard C, Zihler Berner A, Lacroix C. Synergistic effects of *Bifidobacterium thermophilum* RBL67 and selected prebiotics on inhibition of *Salmonella* colonization in the swine proximal colon. PolyFermS model. Gut Pathog. 2014; 6(1):44. doi: 10.1186/s13099-014-0044-y PMID: 25364390; PubMed Central PMCID: PMC4215022.