Dominant and diet-responsive groups of bacteria within the human colonic microbiota (original) (raw)
Journal Article
,
Pathogen Genomics, Wellcome Trust Sanger Institute
,
Cambridge
,
UK
Search for other works by this author on:
,
Rowett Institute of Nutrition and Health, University of Aberdeen
,
Aberdeen
,
UK
Search for other works by this author on:
,
Rowett Institute of Nutrition and Health, University of Aberdeen
,
Aberdeen
,
UK
Search for other works by this author on:
,
Rowett Institute of Nutrition and Health, University of Aberdeen
,
Aberdeen
,
UK
Search for other works by this author on:
,
Biomathematics and Statistics Scotland
,
Aberdeen
,
UK
Search for other works by this author on:
,
Rowett Institute of Nutrition and Health, University of Aberdeen
,
Aberdeen
,
UK
Search for other works by this author on:
,
Rowett Institute of Nutrition and Health, University of Aberdeen
,
Aberdeen
,
UK
Search for other works by this author on:
,
Pathogen Genomics, Wellcome Trust Sanger Institute
,
Cambridge
,
UK
Search for other works by this author on:
,
Pathogen Genomics, Wellcome Trust Sanger Institute
,
Cambridge
,
UK
Search for other works by this author on:
,
Rowett Institute of Nutrition and Health, University of Aberdeen
,
Aberdeen
,
UK
Search for other works by this author on:
Received:
25 January 2010
Revision received:
11 May 2010
Published:
05 August 2010
Cite
Alan W Walker, Jennifer Ince, Sylvia H Duncan, Lucy M Webster, Grietje Holtrop, Xiaolei Ze, David Brown, Mark D Stares, Paul Scott, Aurore Bergerat, Petra Louis, Freda McIntosh, Alexandra M Johnstone, Gerald E Lobley, Julian Parkhill, Harry J Flint, Dominant and diet-responsive groups of bacteria within the human colonic microbiota, The ISME Journal, Volume 5, Issue 2, February 2011, Pages 220–230, https://doi.org/10.1038/ismej.2010.118
Close
Navbar Search Filter Mobile Enter search term Search
Abstract
The populations of dominant species within the human colonic microbiota can potentially be modified by dietary intake with consequences for health. Here we examined the influence of precisely controlled diets in 14 overweight men. Volunteers were provided successively with a control diet, diets high in resistant starch (RS) or non-starch polysaccharides (NSPs) and a reduced carbohydrate weight loss (WL) diet, over 10 weeks. Analysis of 16S rRNA sequences in stool samples of six volunteers detected 320 phylotypes (defined at >98% identity) of which 26, including 19 cultured species, each accounted for >1% of sequences. Although samples clustered more strongly by individual than by diet, time courses obtained by targeted qPCR revealed that ‘blooms’ in specific bacterial groups occurred rapidly after a dietary change. These were rapidly reversed by the subsequent diet. Relatives of Ruminococcus bromii (R-ruminococci) increased in most volunteers on the RS diet, accounting for a mean of 17% of total bacteria compared with 3.8% on the NSP diet, whereas the uncultured Oscillibacter group increased on the RS and WL diets. Relatives of Eubacterium rectale increased on RS (to mean 10.1%) but decreased, along with Collinsella aerofaciens, on WL. Inter-individual variation was marked, however, with >60% of RS remaining unfermented in two volunteers on the RS diet, compared to <4% in the other 12 volunteers; these two individuals also showed low numbers of R-ruminococci (<1%). Dietary non-digestible carbohydrate can produce marked changes in the gut microbiota, but these depend on the initial composition of an individual's gut microbiota.
Introduction
The remarkable diversity of the human colonic microbiota at the level of bacterial species and phylotypes has become apparent from 16S rRNA-based analyses. Not only are hundreds of phylotypes typically estimated to be present in the human colonic microbiota from a given faecal sample (Suau et al., 1999; Eckburg et al., 2005), but samples from different individuals have been reported to show limited overlap in the phylotypes present (Ley et al., 2006; Turnbaugh et al., 2008). A recent study (Tap et al., 2009), however, has indicated that certain phylotypes occur more commonly than others among the dominant faecal bacteria of different individuals. Identifying the dominant bacterial species that colonize the large intestine and the extent to which these are influenced by diet and host factors is of key importance in uncovering the impact of the colonic microbiota upon human health (Flint et al., 2007; Sokol et al., 2008).
A few studies to date have examined temporal changes, and these suggest a degree of stability in the colonic microbiota of individuals consuming their normal diets (Franks et al., 1998; Zoetendal et al., 1998; Costello et al., 2009). In contrast, little is known about the impact of dietary change upon microbial community composition. There is evidence that dietary supplementation with prebiotics such as fructo-oligosaccharides and inulin can promote specific groups of bacteria, including bifidobacteria (Bouhnik et al., 2004; Ramirez-Farias et al., 2009). It has also been shown that reductions in total carbohydrate content, in weight loss (WL) diets for obese subjects, have major effects upon the composition and metabolic outputs of the bacterial community in the colon (Duncan et al., 2007, 2008; Brinkworth et al., 2009). These changes are assumed to reflect the fermentation of non-digestible (ND) carbohydrate components (mainly non-starch polysaccharides (NSPs), resistant starch (RS) and certain oligosaccharides) that reach the large intestine. The impact upon the colonic microbiota of controlled changes in the main types of ND carbohydrate normally present in the diet (RS and NSP) has not however been examined in any detail. The ND carbohydrate content of the human diet is considered to influence health. For example, diets high in RS have been shown to benefit insulin sensitivity, possibly mediated by bacterial fermentative activity in the colon (Robertson et al., 2005). Diets containing RS and NSP offer potential benefits in prevention of colorectal cancer through the delivery of fermentation acids, in particular butyrate, to the distal colon (McIntyre et al., 1993; Duncan et al., 2007). Microbial breakdown of NSP also releases bound phytochemicals into the colon (Gill and Rowland, 2002). These health benefits may be particularly important in obese and overweight subjects who are at increased risk of developing colorectal cancer and diabetes (Polednak, 2003; Cani et al., 2007).
This study asks to what extent changes in the major type of ND carbohydrate in the diet influence the dominant bacterial phylotypes present in an individual. It also explores the dynamics of such changes, including their reversibility, for the first time through the use of detailed time courses and a cross-over design. A necessary supplementary question is whether the same dietary change elicits similar changes in the gut microbiota of different individuals. Our results reveal rapid and marked changes in the colonic microbiota of a group of overweight volunteers after a shift in the main type of ingested fermentable carbohydrate (RS or NSP). We also find that these changes can be highly specific to the individual, with potentially important implications for the design of fibre-enriched diets.
Materials and methods
Volunteer recruitment
Volunteers were initially recruited by newspaper advertisement. Inclusion criteria for men were a body mass index >27 kg m−2, waist circumference >102 cm and fasting plasma glucose concentration >6.0 mmol l−1, indicating metabolic syndrome. They underwent a medical examination and their general practitioner was contacted to confirm medical and medication history. None of the subjects had a history of gastrointestinal disease or disturbance, or took antibiotics, either in the 6 months leading up to the study or during the study. Two of the original 16 subjects left the study for reasons unconnected with the protocol. The remaining 14 volunteers had a mean age of 54 years (range 27–73 years) with a mean body mass index (kg m−2) of 39.4 (range 27.9–51.3) (Supplementary Table S1). The study was approved by the North of Scotland Research Ethics Service and all volunteers provided informed, written consent.
Experimental design and dietary regime
The volunteers were weight stable (<3 kg change in the past 4 months) before entry on the trial. Over the first 7 weeks, intakes were provided at energy maintenance as 1.4–1.5 × measured resting metabolic rate for each individual. The final intervention period consisted of high protein WL diet for an additional 3 weeks (Supplementary Table S2, Supplementary Figure S1). The initial maintenance diet (M diet) comprised protein/carbohydrate/fat % as 13:52:35 of metabolizable energy and 27.7 g per day NSP, provided for 7 days. Subjects were then provided with fixed intakes of two diets, which consisted of either a high RS diet or a high NSP diet each supplied for 3 weeks in a randomized cross-over design (Supplementary Figure S1). The RS diet contained added type III RS whereas the NSP diet intake contained added wheat bran. All meals were of the same energy density (5.5 MJ kg−1) and daily intakes were recorded by weigh-back after each meal (Supplementary Table S2). Daily macronutrient intakes were calculated using the Windiet software program (Robert Gordon University, Aberdeen, UK), based on the type and quantity of each ingredient consumed and published food composition tables (Food Standards Agency 2002). Faecal samples were collected on average twice each week. In addition, one 24 h collection from the final week of each dietary period was used for chemical analysis of digestibility.
Chemical analysis of diet composition and digestibility
Diets were analysed for total gross energy, RS and insoluble and soluble NSPs (Supplementary Table S2). Estimation of RS content of diets and 24 h faecal sample collections was as described by Englyst et al. (1992) whereas NSPs were quantified as described by Englyst and Cummings (1988).
DNA extraction from faecal samples
Faecal samples were kept at 4 °C and processed within 5 h of collection. This short period of storage is not expected to influence molecular estimation of microbial community composition (Lauber et al., 2010). Each sample was mixed and 5 g dispersed (3500 r.p.m. for 1 min using a Dispomix Drive, Medic Tools, Lussiwag, Switzerland) in 10 ml sterile phosphate-buffered saline buffer before aliquoting. One aliquot was used immediately for DNA extraction using the FastDNA Spin for soil kit following manufacturer's (Qbiogene, MP Biomedicals, Illkirch, France) instructions.
Phylogenetic analysis of sequences derived from 16S rRNA clone libraries
PCR amplification of 16S rRNA genes from the extracted DNA involved initial denaturation at 94 °C for 5 min; 20 cycles of denaturation (30 s at 95 °C), annealing (30 s at 52 °C), extension at 72 °C for 2 min with a final cycle at 72 °C for 8 min. Taq1 polymerase from Promega (Southampton, UK) was used at 0.025 U μl−1 in the presence of 2.5 mm MgCl2. The forward primer comprised a 4:1:1:1:1 mixture of 7f, 27fChl, 27fBor, 27fBif and 27fAto and the reverse primer was 1510r (Supplementary Table S3). PCR products were cleaned using the Wizard PCR product purification kit (Promega) and were then cloned and sequenced as described (Lawley et al., 2009). The sequences (spanning variable regions V2–V5) were aligned using the SILVA-derived reference alignment in the mothur software package (Schloss et al., 2009) and were tested for the presence of chimeras by Mallard (Ashelford et al., 2006), Pintail (Ashelford et al., 2005) and BLAST (Johnson et al., 2008). After removal of chimeric sequences, the alignment was subjected to extensive manual curation using the editor function in the ARB package (Ludwig et al., 2004). Using the curated alignment, a distance matrix, with Felsenstein correction, was created using ARB. This matrix was then used as an input for DOTUR (Schloss and Handelsman, 2005) using a 98% identity cut-off under the default furthest-neighbour setting. Sequences with >98% phylogenetic similarity were regarded as belonging to the same phylotype. This generated 320 phylotypes (Supplementary Table S4), which were assigned an identity from the phylum to genus level by the RDP Classifier (Cole et al., 2009) and to the closest match in the NCBI database using MegaBLAST (Johnson et al., 2008). Rarefaction analysis, Shannon and Simpson diversity index calculations, community structure comparisons with ∫-Libshuff, parsimony and weighted and unweighted versions of UniFrac, Jaccard and Yue and Clayton theta tree clustering analysis were all performed using the mothur software package (Schloss et al., 2009).
The final 16S rRNA clone library data set from six volunteers (v16, v19, v20, v22, v23 and v24) contained 5915 sequences (median length of 678 bases) (GenBank accession numbers GU238434–GU244348). These sequences were from one faecal sample from each of the four diets per individual, except for volunteer v16 when two samples were analysed from the M and NSP diets (see Supplementary Table S4). An additional 768 sequences were subsequently obtained from the RS and NSP diets only of v14 and v25 (Supplementary Table S5) (GenBank accession numbers HM191774–HM192541).
qPCR and DGGE analysis of 16S rRNA genes
Quantitative real-time PCR (qPCR) was performed on all available samples from all 14 volunteers. Primer sequences and conditions are reported in Supplementary Table S3, based largely on the study by Ramirez-Farias et al. (2009). Primers for _Oscillibacter_-related bacteria were designed in this study.
For denaturing gradient gel electrophoresis (DGGE), bacterial DNA was amplified by PCR using primers for the hyper-variable V3 region of the 16S gene (Supplementary Table S3). Reaction mixtures contained 1.5 mm MgCl2 and 0.025 U Taq1 polymerase per μl. Products were separated in an 8% acrylamide denaturing gradient gel containing urea and formamide (35–60% gradient) at 80 V for 16 h at 60 °C and visualized by staining with SybrGold dye (Invitrogen, Paisley, UK).
Statistical analysis
For phylotypes with large abundance (more than 1% and found in 5 or 6 out of 6 volunteers, which was the case for 8 out of 320 phylotypes) the Sanger sequences (clone counts) were analysed as Hierarchical Generalized Linear Models using a binomial model with logistic link and estimated overdispersion. This was carried out for one phylotype at a time, with the total number of clones per sample regarded as the total count. Volunteer and diet were regarded as random (normally distributed) and fixed effects, respectively. Similarities between faecal samples were displayed graphically by means of principal coordinate analysis, where the Canberra distance matrix was applied to the numbers of clones assigned to each phylotype (based on all 320 phylotypes) for each faecal sample. Principal coordinate analysis was also used to obtain a graphical representation of the phylogenetic distances (obtained from the DOTUR software package) between the 320 phylotypes.
The effect of diet was tested by means of analysis of variance with volunteer regarded as random and diet regarded as fixed effect. When the effect of diet was significant, diet means were compared by _post hoc t_-test. Data containing many zeros (that is, below limit of detection) were analysed by Friedman's non-parametric analysis of variance instead. For all analyses, the effect of order was included initially as a fixed effect but it and its interaction with diet were found not to be significant and was removed from subsequent analyses. All statistical analyses were performed by using Genstat 11th Edition Release 11.1 (VSN International Ltd, Hemel Hempstead, Hertfordshire, UK). Significance was set at _P_-value of <0.05.
Results
Study design
Fourteen overweight male volunteers completed this study (see Materials and Methods section) (Supplementary Table S1). Complete diets were provided daily resulting in precise control over all dietary intake throughout the 10-week study period (Supplementary Table S2). Each volunteer received a control weight maintenance (M) diet for the first week, followed by three weeks each of diets in high type III resistant starch (RS diet) or high in wheat bran (NSP) and a WL diet reduced in carbohydrate but high in protein (WL diet). The first three diets were provided at energy maintenance and were closely matched for total carbohydrate, protein and fat content; the order of the RS and NSP diets was reversed for half of the volunteers (Supplementary Figure S1, Supplementary Table S2). Faecal samples were processed as soon as possible after collection, without prior freezing, for the extraction of nucleic acids (see Materials and Methods section).
Dominant bacterial phylotypes within the intestinal bacterial community
16S rRNA profiles from all available faecal samples were analysed first using DGGE, which revealed that profiles were consistent over time within an individual for a given diet (Supplementary Figure S2). A switch between the RS and NSP diets, however, altered the profile substantially, indicating a shift in bacterial community composition. Sequences of amplified 16S rRNA genes were analysed for 6 of the 14 volunteers (3 from each diet order) using a sample from the final week of each dietary regime. This yielded 5915 non-chimeric sequences (averaging 227 sequences per sample). Phylogenetic analysis using a 98% sequence identity cut-off revealed 320 phylotypes across all samples (Supplementary Table S4). The proportion of cultured strains increased with increasing abundance of phylotypes, with all phylotypes present at >2% of the total being identified with cultured species (Figure 1). None of the 320 phylotypes was detected in every sample, but 32 phylotypes were present in all six volunteers accounting for 47.1% of all sequences, whereas 77 phylotypes occurred in at least four of the six individuals (Figure 2, Supplementary Figure S3). Figure 2 suggests a bimodal frequency distribution for phylotype occurrence, indicating that there may be a group of widespread and highly abundant phylotypes within the human gut microbiota (Tap et al., 2009). The most abundant phylotype was Eubacterium rectale although three phylotypes corresponding to three divergent strains of Faecalibacterium prausnitzii together comprised 8.0% of total sequences (Table 1).
Figure 1
Culturability of 16S rRNA phylotypes in relation to their abundance. 16S rRNA sequences (5915) obtained from faecal samples of six volunteers (for all four diets) were classified into 320 phylotypes (defined at >98% sequence identity) (Supplementary Table S4). The % of phylotypes showing >98% sequence identity to a cultured bacterium is seen to increase with increasing phylotype abundance. The phylotype frequency distribution is shown as a percentage; actual numbers of phylotypes were 9 (>2% of all sequences), 17 (>1%<2%), 24 (>0.5%<1%) and 270 (<0.5%) (total 320). Data refer to volunteers 16, 20, 22 (diet order RS-NSP) and 19, 23, 24 (diet order NSP-RS).
Figure 2
Incidence of phylotypes in different individuals. The distribution of all 320 16S rRNA phylotypes was determined across the six volunteers (see Figure 1, Supplementary Table S4); the numbers of phylotypes found in all six individuals (32) and in five (25), four (20), three (36), two (62) or one (145) of the six are shown here according to bacterial phylum.
Table 1
Impact of diet upon the 10 most abundant 16S rRNA phylotypes detected in six overweight male volunteers
Phylotype, speciesa | Phylum | Mean % clones | Diets | P_-value_b | |||
---|---|---|---|---|---|---|---|
% total clones | M | NSP | RS | WL | |||
Eubacterium rectale | F (L) | 4.43 | 2.9x | 4.1x | 8.0y | 1.7x | <0.001 |
Collinsella aerofaciens | A | 3.67 | 4.3x | 3.4x | 5.2x | 0.6y | 0.032 |
Clostridium clostridioforme | F (L) | 3.83 | 4.8 | 2.6 | 2.5 | 3.1 | 0.217 |
Bacteroides vulgatus | B | 3.21 | 2.6 | 4.1 | 1.3 | 4.9 | 0.057 |
Faecalibacterium prausnitzii L2-6 | F (R) | 2.96 | 2.5 | 3.4 | 1.6 | 3.8 | 0.224 |
Faecalibacterium prausnitzii A2-165 | F (R) | 2.55 | 1.8 | 2.8 | 2.9 | 3.4 | 0.338 |
Faecalibacterium prausnitzii M21/2 | F (R) | 2.47 | 2.1 | 1.8 | 2.2 | 3.1 | 0.771 |
Anaerostipes colic SSC/2 | F (L) | 2.25 | 3.0 | 2.5 | 1.1 | 1.8 | 0.091 |
Ruminococcus bromii | F (R) | 2.11 | 1.5x | 0.4x | 5.0y | 1.5x | <0.001 |
Eubacterium hallii | F (L) | 2.00 | 2.4 | 2.5 | 1.2 | 1.4 | 0.251 |
Phylotype, speciesa | Phylum | Mean % clones | Diets | P_-value_b | |||
---|---|---|---|---|---|---|---|
% total clones | M | NSP | RS | WL | |||
Eubacterium rectale | F (L) | 4.43 | 2.9x | 4.1x | 8.0y | 1.7x | <0.001 |
Collinsella aerofaciens | A | 3.67 | 4.3x | 3.4x | 5.2x | 0.6y | 0.032 |
Clostridium clostridioforme | F (L) | 3.83 | 4.8 | 2.6 | 2.5 | 3.1 | 0.217 |
Bacteroides vulgatus | B | 3.21 | 2.6 | 4.1 | 1.3 | 4.9 | 0.057 |
Faecalibacterium prausnitzii L2-6 | F (R) | 2.96 | 2.5 | 3.4 | 1.6 | 3.8 | 0.224 |
Faecalibacterium prausnitzii A2-165 | F (R) | 2.55 | 1.8 | 2.8 | 2.9 | 3.4 | 0.338 |
Faecalibacterium prausnitzii M21/2 | F (R) | 2.47 | 2.1 | 1.8 | 2.2 | 3.1 | 0.771 |
Anaerostipes colic SSC/2 | F (L) | 2.25 | 3.0 | 2.5 | 1.1 | 1.8 | 0.091 |
Ruminococcus bromii | F (R) | 2.11 | 1.5x | 0.4x | 5.0y | 1.5x | <0.001 |
Eubacterium hallii | F (L) | 2.00 | 2.4 | 2.5 | 1.2 | 1.4 | 0.251 |
Abbreviations: F, Firmicutes (L, Lachnospiraceae; R, Ruminococcaceae); B, Bacteroidetes; A, Actinobacteria.
a
Detected in all six volunteers, except for R. bromii (four of six volunteers).
b
Within a row values not sharing superscripts (x and y) differ significantly (P< 0.05). Clone counts were analysed with Hierarchical Generalized Linear Model with Volunteer as random and Diet as fixed effect.
c
New species name soon to be formally proposed (Allen-Vercoe et al., in preparation).
Table 1
Impact of diet upon the 10 most abundant 16S rRNA phylotypes detected in six overweight male volunteers
Phylotype, speciesa | Phylum | Mean % clones | Diets | P_-value_b | |||
---|---|---|---|---|---|---|---|
% total clones | M | NSP | RS | WL | |||
Eubacterium rectale | F (L) | 4.43 | 2.9x | 4.1x | 8.0y | 1.7x | <0.001 |
Collinsella aerofaciens | A | 3.67 | 4.3x | 3.4x | 5.2x | 0.6y | 0.032 |
Clostridium clostridioforme | F (L) | 3.83 | 4.8 | 2.6 | 2.5 | 3.1 | 0.217 |
Bacteroides vulgatus | B | 3.21 | 2.6 | 4.1 | 1.3 | 4.9 | 0.057 |
Faecalibacterium prausnitzii L2-6 | F (R) | 2.96 | 2.5 | 3.4 | 1.6 | 3.8 | 0.224 |
Faecalibacterium prausnitzii A2-165 | F (R) | 2.55 | 1.8 | 2.8 | 2.9 | 3.4 | 0.338 |
Faecalibacterium prausnitzii M21/2 | F (R) | 2.47 | 2.1 | 1.8 | 2.2 | 3.1 | 0.771 |
Anaerostipes colic SSC/2 | F (L) | 2.25 | 3.0 | 2.5 | 1.1 | 1.8 | 0.091 |
Ruminococcus bromii | F (R) | 2.11 | 1.5x | 0.4x | 5.0y | 1.5x | <0.001 |
Eubacterium hallii | F (L) | 2.00 | 2.4 | 2.5 | 1.2 | 1.4 | 0.251 |
Phylotype, speciesa | Phylum | Mean % clones | Diets | P_-value_b | |||
---|---|---|---|---|---|---|---|
% total clones | M | NSP | RS | WL | |||
Eubacterium rectale | F (L) | 4.43 | 2.9x | 4.1x | 8.0y | 1.7x | <0.001 |
Collinsella aerofaciens | A | 3.67 | 4.3x | 3.4x | 5.2x | 0.6y | 0.032 |
Clostridium clostridioforme | F (L) | 3.83 | 4.8 | 2.6 | 2.5 | 3.1 | 0.217 |
Bacteroides vulgatus | B | 3.21 | 2.6 | 4.1 | 1.3 | 4.9 | 0.057 |
Faecalibacterium prausnitzii L2-6 | F (R) | 2.96 | 2.5 | 3.4 | 1.6 | 3.8 | 0.224 |
Faecalibacterium prausnitzii A2-165 | F (R) | 2.55 | 1.8 | 2.8 | 2.9 | 3.4 | 0.338 |
Faecalibacterium prausnitzii M21/2 | F (R) | 2.47 | 2.1 | 1.8 | 2.2 | 3.1 | 0.771 |
Anaerostipes colic SSC/2 | F (L) | 2.25 | 3.0 | 2.5 | 1.1 | 1.8 | 0.091 |
Ruminococcus bromii | F (R) | 2.11 | 1.5x | 0.4x | 5.0y | 1.5x | <0.001 |
Eubacterium hallii | F (L) | 2.00 | 2.4 | 2.5 | 1.2 | 1.4 | 0.251 |
Abbreviations: F, Firmicutes (L, Lachnospiraceae; R, Ruminococcaceae); B, Bacteroidetes; A, Actinobacteria.
a
Detected in all six volunteers, except for R. bromii (four of six volunteers).
b
Within a row values not sharing superscripts (x and y) differ significantly (P< 0.05). Clone counts were analysed with Hierarchical Generalized Linear Model with Volunteer as random and Diet as fixed effect.
c
New species name soon to be formally proposed (Allen-Vercoe et al., in preparation).
Overall community structure was analysed for each sample using the mothur software package (Schloss et al., 2009). ∫-Libshuff, parsimony and UniFrac (weighted and unweighted) analysis revealed a high degree of variation between individuals and showed that each sample was significantly different from all others. Similarity in community membership and structure was estimated using the Jaccard and the Yue and Clayton theta similarity coefficients (Yue and Clayton, 2005) (Supplementary Figure S4). These results and those from principal coordinate analysis (Figure 3) indicate that samples predominantly clustered by volunteer, suggesting that responses to dietary change are influenced by the initial species composition of the individual's colonic microbiota. Rarefaction analysis revealed that the observed species richness was similar for all four diets (Supplementary Figure S5). Simpson (1/D) and Shannon diversity indices were slightly lower on the RS and NSP diets compared with the M and WL diets, but did not differ significantly (Supplementary Figure S5; see also Supplementary Table S4 for individual values).
Figure 3
Impact of diet and individual variation upon faecal microbiota composition. A principal coordinates analysis (using Canberra distance matrix) based on 16S rRNA clone libraries from 26 faecal samples (obtained from six donors under four dietary conditions) (Supplementary Table S4). Colour code is based on donor (v16, v19, v20, v22, v23 and v24). Diets are indicated as M, maintenance; NSP, non-starch polysaccharide; RS, resistant starch; WL, weight loss. Two samples were analysed for v16 from the M and NSP diets.
Response of the colonic microbial community to dietary change
At the phylum level there was no significant effect of diet upon the proportions of Bacteroidetes (mean 21.5%), Firmicutes (mean 70.6%), Actinobacteria (mean 4.9%) or Proteobacteria (mean 3.0%) within the faecal microbiota. At finer taxonomic levels, however, two individual phylotypes, E. rectale and Ruminococcus bromii, showed increased proportions on the RS diet whereas Collinsella aerofaciens showed decreased proportions on the WL diet (P<0.001) (Table 1). Because the sequence data also suggested that groups of related phylotypes might be affected by diet, we decided to use qPCR for further analysis.
PCR primer sets (Supplementary Table S3) were used to target groups for which sequence analysis indicated a response to the RS diet (Rum, Ros, Osc), or for which previous in vitro evidence has shown the ability to use starch (Bac, Bif, Fprau) (Table 2). Methanogens were also monitored because of their potential impact on the polysaccharide-utilizing community (Robert and Bernalier-Donadille, 2003). All available time points were analysed. Mean data for the targeted groups are shown in Table 2 and complete time courses for all 14 subjects are shown in Figure 4 for those bacterial groups that showed responses to the RS diet. The Rum group includes relatives of R. bromii, R. flavefaciens and R. albus that belong to the Ruminococcaceae, and will be referred to here as ‘R-ruminococci’. There was a significant increase of around 4.5-fold (P<0.001) in this group on the RS diet, compared to the NSP diet (Table 2, Figure 4). A new primer set was designed to recognize another group of Ruminococcaceae, related to Oscillibacter valericigenes (Iino et al., 2007) that has not been cultured from the human gut. This group increased significantly on average on the RS and WL relative to the M and NSP diets (P<0.001) (Table 2, Figure 4). A previously reported primer set (Mackie et al., 2003) targeting relatives of this group from the rumen was also used here, but detected <0.1% of human faecal bacteria (results not shown). Populations of a third abundant group of Ruminococcaceae, F. prausnitzii, showed no significant response to changes of diet (Table 2).
Table 2
Populations estimated by qPCR of bacteria and methanogenic archaea expressed as % relative to total bacteria
Bacterial group (%) | Diet | P dieta | ||||
---|---|---|---|---|---|---|
M | NSP | RS | WL | SED | ||
Bac | 27.8 | 25.7 | 20.2 | 24.0 | 3.58 | 0.387 |
Ros | 7.3x | 6.5x | 10.1y | 3.3z | 1.381 | <0.001 |
Fprau | 11.2 | 14.4 | 12.1 | 12.5 | 1.622 | 0.278 |
Rum | 6.5x | 3.8x | 17.0y | 7.5x | 2.253 | <0.001 |
Osc | 0.74x | 0.77x | 2.0y | 1.6y | 0.34 | <0.001 |
Bif | 1.9 | 1.8 | 2.4 | 0.8 | 0.632 | 0.059 |
Met | 0.080 | 0.034 | 0.124 | 0.123 | 0.042 | 0.155b |
Bacterial group (%) | Diet | P dieta | ||||
---|---|---|---|---|---|---|
M | NSP | RS | WL | SED | ||
Bac | 27.8 | 25.7 | 20.2 | 24.0 | 3.58 | 0.387 |
Ros | 7.3x | 6.5x | 10.1y | 3.3z | 1.381 | <0.001 |
Fprau | 11.2 | 14.4 | 12.1 | 12.5 | 1.622 | 0.278 |
Rum | 6.5x | 3.8x | 17.0y | 7.5x | 2.253 | <0.001 |
Osc | 0.74x | 0.77x | 2.0y | 1.6y | 0.34 | <0.001 |
Bif | 1.9 | 1.8 | 2.4 | 0.8 | 0.632 | 0.059 |
Met | 0.080 | 0.034 | 0.124 | 0.123 | 0.042 | 0.155b |
Abbreviations: Bac, Bacteroides/Prevotella; Ros, E. rectale/Roseburia spp.; Fprau, F. prausnitzii; Rum, R-ruminococci; Osc, Oscillibacter relatives; Bif, Bifidobacterium spp.; Met, methanogens.
a
Significance based on analysis of variance (ANOVA). Within each row, values not superscripts (x, y and z) differ significantly (P<0.05, post hoc _t_-test).
b
Friedman non-parametric ANOVA.
Primer sequences and amplification conditions are shown in Supplementary Table S3. One mean value is shown per diet for the 14 volunteers. For the M diet, means are for all samples from the 1 week of this diet, whereas for the NSP, RS and WL diets means are for all samples from the last 2 weeks on each diet.
Table 2
Populations estimated by qPCR of bacteria and methanogenic archaea expressed as % relative to total bacteria
Bacterial group (%) | Diet | P dieta | ||||
---|---|---|---|---|---|---|
M | NSP | RS | WL | SED | ||
Bac | 27.8 | 25.7 | 20.2 | 24.0 | 3.58 | 0.387 |
Ros | 7.3x | 6.5x | 10.1y | 3.3z | 1.381 | <0.001 |
Fprau | 11.2 | 14.4 | 12.1 | 12.5 | 1.622 | 0.278 |
Rum | 6.5x | 3.8x | 17.0y | 7.5x | 2.253 | <0.001 |
Osc | 0.74x | 0.77x | 2.0y | 1.6y | 0.34 | <0.001 |
Bif | 1.9 | 1.8 | 2.4 | 0.8 | 0.632 | 0.059 |
Met | 0.080 | 0.034 | 0.124 | 0.123 | 0.042 | 0.155b |
Bacterial group (%) | Diet | P dieta | ||||
---|---|---|---|---|---|---|
M | NSP | RS | WL | SED | ||
Bac | 27.8 | 25.7 | 20.2 | 24.0 | 3.58 | 0.387 |
Ros | 7.3x | 6.5x | 10.1y | 3.3z | 1.381 | <0.001 |
Fprau | 11.2 | 14.4 | 12.1 | 12.5 | 1.622 | 0.278 |
Rum | 6.5x | 3.8x | 17.0y | 7.5x | 2.253 | <0.001 |
Osc | 0.74x | 0.77x | 2.0y | 1.6y | 0.34 | <0.001 |
Bif | 1.9 | 1.8 | 2.4 | 0.8 | 0.632 | 0.059 |
Met | 0.080 | 0.034 | 0.124 | 0.123 | 0.042 | 0.155b |
Abbreviations: Bac, Bacteroides/Prevotella; Ros, E. rectale/Roseburia spp.; Fprau, F. prausnitzii; Rum, R-ruminococci; Osc, Oscillibacter relatives; Bif, Bifidobacterium spp.; Met, methanogens.
a
Significance based on analysis of variance (ANOVA). Within each row, values not superscripts (x, y and z) differ significantly (P<0.05, post hoc _t_-test).
b
Friedman non-parametric ANOVA.
Primer sequences and amplification conditions are shown in Supplementary Table S3. One mean value is shown per diet for the 14 volunteers. For the M diet, means are for all samples from the 1 week of this diet, whereas for the NSP, RS and WL diets means are for all samples from the last 2 weeks on each diet.
Figure 4
Diet-driven changes in four groups of human colonic bacteria detected by qPCR. Abundance for each targeted group is expressed as a percentage of the signal obtained with a general bacterial primer set (see Supplementary Table S3 for primers and conditions used). All available time points are shown for the 14 volunteers: left-hand panels show the diet order M-NSP-RS-WL and right-hand panels M-RS-NSP-WL, (a and b) R-ruminococci (relatives of R. bromii); (c and d) relatives of E. rectale and Roseburia spp.; (e and f) Bifidobacterium spp.; (g and h) relatives of O. valericigenes.
Firmicutes bacteria related to Roseburia and E. rectale also increased significantly (P<0.048) on the RS diet, and decreased (approximately threefold compared to RS) (P<0.029) on the reduced carbohydrate WL diet. (Ros, Table 2). Individual differences were evident for all groups shown in Figure 4. Notably, Bifidobacterium spp. did not show a significant mean response to RS, but their population increased markedly on RS in one volunteer (v23) (Figures 4 and 5). Bifidobacteria are often poorly represented in 16S rRNA clone libraries (Eckburg et al., 2005; Hold et al., 2002) and this was also the case in our analysis in spite of using a broad bacterial primer set reported to include this group (Frank et al., 2008). There was no significant change in the % Bacteroides in response to diet (Table 2).
Figure 5
Populations of three groups of potentially amylolytic bacteria on the RS (high-resistant starch) and NSP (low-resistant starch) diets. Populations, estimated by qPCR, are shown for R-ruminococci, E. rectale/Roseburia spp. and Bifidobacterium spp. for all 14 volunteers, according to diet order. Data are the mean values for each volunteer during the second and third week of the RS and NSP diets (overall means are given in Table 2).
Methanogenic archaea were detected by qPCR in >90% of samples from nine volunteers, but were detected in fewer than 10% of samples from v11, 14, 17, 24 and 25. This group showed no significant effect of diet (Table 2).
Digestibility of RS and NSP
Whole tract digestibility of RS and NSP was estimated from chemical analysis of the diet and of 24 h faecal collections. For 12 of the 14 volunteers <4% of ingested RS was recovered in faecal samples for all four diets, reflecting almost complete microbial fermentation. For two individuals (v14 and v25), however, faecal RS recovery was 69% and 65%, respectively (31% and 35% digestibility) on the RS diet (Figure 6). These two individuals both harboured low R-ruminococci numbers as assessed by qPCR, although one of them showed an increase in the Roseburia and E. rectale group on the RS diet (Figures 4 and 5). Additional clone library analyses were subsequently performed on the final samples from the RS and NSP dietary periods for these two individuals. No R-ruminococci were detected whereas 36% and 46.5% of clones belonged to the Bacteroidetes among the 381 and 387 clones analysed for v14 and v25, respectively (Supplementary Table S5).
Figure 6
Starch digestibilities. The whole tract % digestibility of resistant starch was determined after estimating resistant starch present in the supplied diet and in 24 h faecal samples (see Materials and Methods section). Results are shown only for the NSP and RS diet periods. The markedly reduced resistant starch digestibilities in v14 and v25 correspond with low R-ruminococcal numbers in these individuals (see text and Supplementary Table S5).
Soluble NSP digestibility averaged 90% and was slightly lower (P<0.05) for RS and WL diets (88%) than for the M diet (92%). Insoluble NSP digestibility averaged 66% and showed no difference between diets (mean range 58–71%).
Discussion
This is the first study to combine precise control and monitoring of human dietary intake and digestibility with detailed analysis of changes in the faecal microbiota at the level of bacterial phylotypes. As a result it has been possible to reveal the interplay between diet and inter-individual differences that influence the human colonic community. Analysis of amplified 16S rRNA sequences suggested that samples from six overweight male volunteers clustered more strongly by individual than by diet. This is in line with previously reported evidence of inter-individual variation in the faecal microbiota (Franks et al., 1998; Ley et al., 2006). In contrast, diet composition had very substantial effects on specific groups of bacteria that were detected both through clone libraries and extensive qPCR analysis on 279 faecal samples. The time courses show that most diet-driven changes occurred rapidly, being detectable within 3–4 days, and were reversed equally rapidly. These kinetics appear consistent with immediate effects of dietary residue upon relative bacterial growth rates in the colon that are subsequently reflected in faecal samples as colonic contents turnover, assuming mean transit colonic times of around 60 h (Stephen et al., 1987).
Firmicutes bacteria related to R. bromii (R-ruminococci) and E. rectale were commonly stimulated by the RS diet. In most individuals, qPCR analysis revealed a surge in the population of R-ruminococci with values exceeding 25% of total bacteria in some samples. R-ruminococci were the only group of human gut bacteria previously found to be preferentially associated with particulate material in human faecal samples (Walker et al., 2008), suggesting that they have an important role in the breakdown of particulate substrates. R. bromii isolates show amylolytic activity (Salyers et al., 1977) but our results suggest that many related, but uncultured, bacteria may also possess this activity. One previous study, using non-quantitative DGGE analysis, reported that _R. bromii_-related bacteria were prominent in faecal samples from humans on a diet high in RS and NSP (Abell et al., 2008). The E. rectale and Roseburia group is also known to contain amylolytic species (Ramsay et al., 2006) but has not previously been shown to respond to RS in vivo. An earlier study failed to detect an increase in specific Eubacterium spp. in human subjects on a diet high in RS type III but, importantly, detection methods were not available for E. rectale (Schwiertz et al., 2002). In this work we also detected a dietary response for an uncultured group of Ruminococcaceae related to Oscillibacter. It is not known whether these bacteria are starch degraders, and their increase both on the RS and WL diets suggests that other factors must be involved in their response to diet.
Early studies indicated a broad distribution of amylase activity among cultured human gut bacteria with >50% of strains, including many Bacteroides, able to grow on amylose or amylopectin starch (Salyers et al., 1977). The potential role of Bifidobacterium spp. has been emphasized particularly in the degradation of high amylose RS (Macfarlane and Englyst, 1986; Wang et al., 1999). In this study, bifidobacterial numbers showed a strong response to RS only in one individual (v23), whereas there was no evidence for a response of the Bacteroides group to RS. Bacteroides spp. may perhaps be better adapted to the use of solubilized starch molecules (Flint et al., 2008). These findings in vivo however correspond better with more recent work in vitro that identified R. bromii, E. rectale and Bifidobacterium spp. as the major species in human faeces that colonize starch particles (Leitch et al., 2007) and that use 13C-labelled starch (Kovatcheva-Datchary et al., 2009).
Different groups of amylolytic bacteria differ markedly in their metabolic products and potential host interactions. Bifidobacterium spp., which produce lactate and acetate, are widely used as probiotics and as targets for prebiosis (Furrie et al., 2005). Members of the E. rectale group are flagellated bacteria that are major producers of butyrate in the large intestine, and may therefore contribute to the butyrogenic effect of RS (Aminov et al., 2006; Duncan et al., 2007; Louis et al., 2010). R. bromii is a producer of acetate, ethanol and hydrogen that is likely to contribute to gas production, but otherwise little is currently known about the impact of this group upon the host. It will be important in future to establish whether different types of RS select for different groups of amylolytic bacteria in view of their differing effects on fermentation and host responses (Le Leu et al., 2009).
Starch is considered to show high digestibility across the whole gastrointestinal tract, with most RSs being completely fermented in the large intestine (Bird et al., 2000). Remarkably, however, significant amounts of starch survived fermentation to be recovered in the faeces in 2 of the 14 individuals studied here. This difference in fermentation cannot be ascribed to dietary intake, which was standardized for all 14 subjects, and seems likely to lie with the strain composition of each individual's colonic microbiota. These two individuals showed very low numbers of R-ruminococci, and the relationship between community composition and starch fermentation will clearly warrant further investigation. Although these two individuals were also non-methanogenic, no simple relationship was evident between methanogens and R-ruminococci. Methanogen populations are known to vary between individuals and are influenced by a variety of factors (Florin et al., 2000).
In contrast to these responses to RS, there was little evidence that the high NSP diet resulted in major alterations in the composition of the faecal microbiota. In part, however, this may reflect the fact that a smaller increase was achieved in NSP intake (1.5-fold) than with RS intake (4.8-fold) when compared with the M diet. It is possible that larger changes in specific NSP components would affect the populations of specific groups of colonic bacteria, as was observed with RS. Significant decreases were observed for C. aerofaciens and for the E. rectale group on the WL diet. The WL changes do not show a simple relationship with RS and NSP intakes, and it is possible that the increased dietary protein content of this diet might have a role in altering microbiota composition.
The distribution of major bacterial phyla observed here did not depart dramatically from that reported in non-obese subjects (Eckburg et al., 2005; Walker et al., 2008). The data that we obtained by qPCR and clone library analysis are consistent with recent reports on obese subjects using fluorescence in situ hybridization microscopy (Duncan et al., 2007, 2008; Schwiertz et al., 2010) or 16S rRNA sequencing (Zhang et al., 2008) although lower % Bacteroides have been reported in another study (Ley et al., 2006). More subtle differences may occur between the gut microbiota of obese and non-obese individuals at the species level, and indeed the dietary responses reported here make this likely. Nevertheless 5 of the 10 most abundant phylotypes identified in this study group were also among the 10 most abundant phylotypes described by Tap et al. (2009) in non-obese volunteers.
A high proportion of the most abundant 16S rRNA phylotypes that we detected here corresponded to cultured bacteria, with 66% of the 50 phylotypes that accounted for >0.5% of sequences having close cultured relatives. This suggests that the limited coverage of the human gut microbiota through cultivation may be because of insufficient anaerobic isolation work, rather than to intrinsic non-culturability of human colonic bacteria. Four of the five most abundant phylotypes detected here (E. rectale, F. prausnitzii, C. aerofaciens and B. vulgatus) corresponded with the most abundant species reported by anaerobic cultivation in Japanese–Hawaiians, North American Caucasians and polyp patients (Moore and Moore, 1995). Overall, however, 33.4% of phylotypes showed <98% identity with cultured bacteria in this study, reflecting the relatively poor coverage of the less abundant bacterial groups by cultured strains. Many abundant phylotypes were found to be widespread, with 32 of the 320 phylotypes detected being present in all six individuals surveyed and accounting for almost half (47%) of all sequences.
In conclusion, we find that increased intake of RS, an important ND carbohydrate in the human diet, can substantially alter the species composition of the colonic microbiota. Such responses are likely to occur also with other fermentable ND dietary components. It follows that the colonic microbial community must typically be in a state of continuous change over time, driven by short-term changes in dietary intake. Thus, only the most successful and versatile organisms will be found commonly among the dominant microbiota at different sampling times and in different individuals. In addition, however, our evidence shows that the bacterial strain composition of the colonic microbiota is subject to inter-individual variation. Furthermore, we show that this can be associated with profound inter-individual differences in the response of the microbial community to dietary change, and in microbial fermentation of dietary substrates in the colon. This suggests that dietary advice on the consumption of ND carbohydrates might need to be personalized in the future.
Acknowledgements
We acknowledge support for this project from the World Cancer Research Fund. The Rowett Institute of Nutrition and Health (University of Aberdeen) and Biomathematics and Statistics Scotland received support from the Scottish Government Rural Environment Research and Analysis Directorate. We thank Claire Fyfe and the RINH Human Studies Unit. We also thank David Harris and the rest of his team at the Wellcome Trust Sanger Institute for generating the 16S rRNA gene sequences. Funding for AWW, MDS, PS, JP and sequencing was provided by The Wellcome Trust (grant no. WT 76964).
Data availability
Accession codes
Accessions
GenBank/EMBL/DDBJ
GenBank/EMBL/DDBJ
References
Abell
GCJ
,
Cooke
CM
,
Bennett
CN
,
Conlon
MA
,
McOrist
AL
. (
2008
).
Phylotypes related to Ruminococcus bromii are abundant in the large bowel of humans and increase in response to a diet high in resistant starch
.
FEMS Microbiol Ecol
66
:
505
–
515
.
Aminov
RI
,
Walker
AW
,
Duncan
SH
,
Harmsen
HJM
,
Welling
GW
,
Flint
HJ
. (
2006
).
Molecular diversity, cultivation, and improved FISH detection of a dominant group of human gut bacteria related to Roseburia spp. or Eubacterium rectale
.
Appl Environ Microbiol
72
:
6371
–
6376
.
Ashelford
KE
,
Chuzhanova
NA
,
Fry
JC
,
Jones
AJ
,
Weightman
AJ
. (
2005
).
At least 1 in 20 sequence records currently held in public repositories is estimated to contain substantial anomalies
.
Appl Environ Microbiol
71
:
7724
–
7736
.
Ashelford
KE
,
Chuzhanova
NA
,
Fry
JC
,
Jones
AJ
,
Weightman
AJ
. (
2006
).
New screening software shows that most recent large 16S rRNA gene clone libraries contain chimeras
.
Appl Environ Microbiol
72
:
5734
–
5741
.
Bird
AR
,
Brown
IL
,
Topping
DL
. (
2000
).
Starches, resistant starches, the gut microflora and human heath
.
Curr Issues Intest Microbiol
1
:
25
–
37
.
Bouhnik
Y
,
Raskine
L
,
Simoneau
G
,
Vicaut
E
,
Neut
C
,
Flourie
B
et al. (
2004
).
The capacity of nondigestible carbohydrates to stimulate faecal bifidobacteria in healthy humans: a double blind, randomized, placebo-controlled, parallel-group, dose response relation study
.
Am J Clin Nutr
80
:
1658
–
1664
.
Brinkworth
GD
,
Noakes
M
,
Clifton
PM
,
Bird
AR
. (
2009
).
Comparative effects of very low carbohydrate, high fat and high carbohydrate, low-fat weight loss diets on bowel habit and faecal short chain fatty acids and bacterial populations
.
Brit J Nutr
101
:
1493
–
1502
.
Cani
PD
,
Amar
J
,
Iglesias
MA
,
Poggi
M
,
Knauf
C
,
Bastelica
D
et al. (
2007
).
Metabolic endotoxemia initiates obesity and insulin resistance
.
Diabetes
56
:
1761
–
1772
.
Cole
JR
,
Wang
Q
,
Cardenas
E
,
Fish
J
,
Chai
B
,
Farris
RJ
et al. (
2009
).
The Ribosomal Database Project: improved alignments and new tools for rRNA analysis
.
Nucleic Acids Res
37
:
D141
–
D145
.
Costello
EK
,
Lauber
CL
,
Hamady
M
,
Fierer
N
,
Gordon
JI
,
Knight
R
. (
2009
).
Bacterial community variation in human body habitats across time and space
.
Science
326
:
1694
–
1697
.
Duncan
SH
,
Belenguer
A
,
Holtrop
G
,
Johnstone
AM
,
Flint
HJ
,
Lobley
GE
. (
2007
).
Reduced dietary intake of carbohydrates by obese subjects results in decreased concentrations of butyrate and butyrate-producing bacteria in feces
.
Appl Environ Microbiol
73
:
1073
–
1078
.
Duncan
SH
,
Lobley
GE
,
Holtrop
G
,
Ince
J
,
Johnstone
AM
,
Louis
P
et al. (
2008
).
Human colonic microbiota associated with diet, obesity and weight loss
.
Int J Obesity
32
:
1720
–
1724
.
Eckburg
PB
,
Bernstein
CN
,
Purdom
E
,
Dethlefsen
L
,
Sargent
M
,
Gill
SR
et al. (
2005
).
Diversity of the human intestinal microbial flora
.
Science
308
:
1635
–
1638
.
Englyst
HN
,
Cummings
JH
. (
1988
).
Improved method for measurement of dietary fiber as non-starch polysaccharides in plant foods
.
J Assoc Off Anal Chem
1
:
247
–
286
.
Englyst
HN
,
Kingman
SM
,
Cummings
JE
. (
1992
).
Classification and measurement of nutritionally important starch fractions
.
Eur J Clin Nutr
46
:
Suppl 2
S33
–
S50
.
Food Standards Agency (2002). McCance and Widdowson's The Composition of Foods, 6th summary edn. Royal Society of Chemistry: Cambridge..
Flint
HJ
,
Bayer
EA
,
Rincon
MT
,
Lamed
R
,
White
BA
. (
2008
).
Polysaccharide utilization by gut bacteria: potential for new insights from genomic analysis
.
Nat Rev Microbiol
6
:
121
–
131
.
Flint
HJ
,
Duncan
SH
,
Scott
KP
,
Louis
P
. (
2007
).
Interactions and competition within the microbial community of the human colon: links between diet and health
.
Environ Microbiol
9
:
1101
–
1111
.
Florin
THF
,
Zhu
G
,
Kirk
KM
,
Martin
NG
. (
2000
).
Shared and unique environmental factors determine the ecology of methanogens in humans and rats
.
Am J Gastroenterol
95
:
2872
–
2879
.
Frank
JA
,
Reich
CI
,
Sharma
S
,
Weisbaum
JS
,
Wilson
BA
,
Olsen
GJ
. (
2008
).
Critical evaluation of two primers commonly used for amplification of bacterial 16S rRNA genes
.
Appl Environ Microbiol
74
:
2461
–
2470
.
Franks
AH
,
Harmsen
HJM
,
Raangs
GC
,
Jansen
GJ
,
Schut
F
,
Welling
GW
. (
1998
).
Variations of bacterial populations in human feces measured by fluorescent in situ hybridisation with group-specific 16S rRNA-targeted oligonucleotide probes
.
Appl Environ Microbiol
64
:
3336
–
3345
.
Furrie
E
,
Macfarlane
S
,
Kennedy
A
,
Cummings
JH
,
Walsh
SV
,
O’Neill
DA
et al. (
2005
).
Synbiotic therapy (Bifidobacterium longum/Synergy 1) initiates resolution of inflammation in patients with active ulcerative colitis: a randomised controlled pilot trial
.
Gut
54
:
242
–
249
.
Gill
CIR
,
Rowland
IR
. (
2002
).
Diet and cancer: assessing the risk
.
Br J Nutr
88
:
Suppl 1
S73
–
S87
.
Hold
GL
,
Pryde
SE
,
Russell
VJ
,
Furrie
E
,
Flint
HJ
. (
2002
).
Assessment of microbial diversity in human colonic samples by 16S rDNA sequence analysis
.
FEMS Microbiol Ecol
39
:
33
–
39
.
Johnson
M
,
Zaretskaya
I
,
Raytselis
Y
,
Merezhuk
Y
,
McGinnis
S
,
Madden
TL
. (
2008
).
NCBI BLAST: a better web interface
.
Nucleic Acids Res
36
:
W5
–
W9
.
Kovatcheva-Datchary
P
,
Egert
M
,
Maathuis
A
,
Rajilic-Stojanovic
M
,
de Graaf
AA
,
Smidt
H
et al. (
2009
).
Linking phylogenetic identities of bacteria to starch fermentation in an in vitro model of the large intestine by RNA-based stable isotope probing
.
Environ Microbiol
11
:
914
–
926
.
Iino
T
,
Mori
K
,
Tanaka
K
,
Suzuki
K-I
,
Harayama
S
. (
2007
).
Oscillibacter valericigenes gen nov., a valerate-producing anaerobic bacterium isolated from the alimentary canal of a Japanese corbicula clam
.
Int J Syst Evol Microbiol
57
:
1840
–
1845
.
Lauber
CL
,
Zhou
N
,
Gordon
JI
,
Knight
R
,
Fierer
N
. (
2010
).
Effect of storage conditions on the assessment of bacterial community structure in soil and human associated samples
.
FEMS Microbiol Lett
307
:
80
–
85
.
Lawley
TD
,
Clare
S
,
Walker
AW
,
Goulding
D
,
Stabler
RA
,
Croucher
N
et al. (
2009
).
Antibiotic treatment of Clostridium difficile carrier mice triggers a supershedder state, spore-mediated transmission, and severe disease in immunocompromised hosts
.
Infect Immun
77
:
3661
–
3669
.
Leitch
ECM
,
Walker
AW
,
Duncan
SH
,
Holtrop
G
,
Flint
HJ
. (
2007
).
Selective colonization of insoluble substrates by human colonic bacteria
.
Environ Microbiol
72
:
667
–
679
.
Le Leu
RK
,
Hu
Y
,
Brown
IL
,
Young
GP
. (
2009
).
Effect of high amylose maize starches on colonic fermentation and apoptotic response to DNA-damage in the colon of rats
.
BMC Nutr Metab
6
:
11
–
19
.
Ley
RE
,
Turnbaugh
PJ
,
Klein
S
,
Gordon
JI
. (
2006
).
Microbial ecology—human gut microbes associated with obesity
.
Nature
444
:
1022
–
1023
.
Louis
P
,
Young
P
,
Holtrop
G
,
Flint
HJ
. (
2010
).
Diversity of human colonic butyrate-producing bacteria revealed by analysis of the butyryl-CoA:acetate CoA transferase gene
.
Environ Microbiol
12
:
304
–
314
.
Ludwig
W
,
Strunk
O
,
Westram
R
,
Richter
L
,
Meier
H
,
Yadhukumar
et al. (
2004
).
ARB: a software environment for sequence data
.
Nucleic Acids Res
32
:
1363
–
1371
.
Macfarlane
GT
,
Englyst
HN
. (
1986
).
Starch utilization by the human large intestinal microflora
.
J Appl Bacteriol
60
:
195
–
201
.
Mackie
RI
,
Aminov
RI
,
Hu
W
,
Klieve
AV
,
Ouwekerk
D
,
Sundset
MA
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
69
:
6808
–
6815
.
McIntyre
AP
,
Gibson
P
,
Young
GP
. (
1993
).
Butyrate production from dietary fibre and protection against large bowel cancer in a gut model
.
Gut
34
:
386
–
391
.
Moore
WEC
,
Moore
LH
. (
1995
).
Intestinal floras of populations that have a high risk of colon cancer
.
Appl Environ Microbiol
61
:
3202
–
3207
.
Polednak
AP
. (
2003
).
Trends in incidence rates for obesity-associated cancers in the U.S
.
Cancer Detect Prev
27
:
415
–
421
.
Ramirez-Farias
C
,
Slezak
K
,
Fuller
Z
,
Duncan
A
,
Holtrop
G
,
Louis
P
. (
2009
).
Effect of inulin on the human gut microbiota: stimulation of Bifidobacterium adolescentis and Faecalibacterium prausnitzii
.
Br J Nutr
101
:
541
–
550
.
Ramsay
AG
,
Scott
KP
,
Martin
CJ
,
Rincon
MT
,
Flint
HJ
. (
2006
).
Cell-associated α-amylases of butyrate-producing Firmicute bacteria from the human colon
.
Microbiology
152
:
3281
–
3290
.
Robert
C
,
Bernalier-Donadille
A
. (
2003
).
The cellulolytic microflora of the human colon: evidence of microcrystalline cellulose-degrading bacteria in methane-excreting subjects
.
FEMS Microbiol Ecol
46
:
81
–
89
.
Robertson
MD
,
Bickerton
AS
,
Dennis
AL
,
Vidal
H
,
Frayn
KN
. (
2005
).
Insulin-sensitizing effects of dietary resistant starch and effects on skeletal muscle and adipose tissue metabolism
.
Am J Clin Nutr
82
:
559
–
567
.
Salyers
AA
,
West
SEH
,
Vercellotti
JR
,
Wilkins
TD
. (
1977
).
Fermentation of mucins and plant polysaccharides by anaerobic bacteria from the human colon
.
Appl Environ Microbiol
34
:
529
–
533
.
Schloss
PD
,
Handelsman
J
. (
2005
).
Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness
.
Appl Environ Microbiol
71
:
1501
–
1506
.
Schloss
PD
,
Westcott
SL
,
Ryabin
T
,
Hall
JR
,
Hartmann
M
,
Hollister
EB
et al. (
2009
).
Introducing mothur: open source, platform-independent, community-supported software for describing and comparing microbial communities
.
Appl Environ Microbiol
75
:
7537
–
7541
.
Schwiertz
A
,
Lehmann
U
,
Jacobasch
G
,
Blaut
M
. (
2002
).
Influence of resistant starch on the SCFA production and cell counts of butyrate-producing Eubacterium spp. in the human intestine
.
J Appl Microbiol
93
:
157
–
162
.
Schwiertz
A
,
Tara
D
,
Schafer
K
,
Beijer
S
,
Bos
NA
,
Donus
C
et al. (
2010
).
Microbiota and SCFA in lean and overweight healthy subjects
.
Obesity
18
:
190
–
195
.
Sokol
H
,
Pigneur
B
,
Watterlot
L
,
Lakhdari
O
,
Bermúdez-Humarán
LG
,
Gratadoux
J
et al. (
2008
).
Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn's disease patients
.
Proc Natl Acad Sci USA
105
:
16731
–
16736
.
Stephen
AM
,
Wiggins
HS
,
Cummings
JH
. (
1987
).
Effect of changing transit time on colonic microbial metabolism in man
.
Gut
28
:
601
–
609
.
Suau
A
,
Bonnet
R
,
Sutren
M
,
Godon
JJ
,
Gibson
GR
,
Collins
MD
et al. (
1999
).
Direct analysis of genes encoding 16S rRNA from complex communities reveals many novel molecular species within the human gut
.
Appl Environ Microbiol
24
:
4799
–
4807
.
Tap
J
,
Mondot
S
,
Levenez
F
,
Pelletier
E
,
Caron
C
,
Furet
JP
et al. (
2009
).
Towards the human intestinal microbiota phylogenetic core
.
Environ Microbiol
11
:
2574
–
2584
.
Turnbaugh
PJ
,
Backhed
F
,
Fulton
L
,
Gordon
JI
. (
2008
).
Diet-induced obesity is linked to marked but reversible alterations in the mouse distal gut microbiome
.
Cell Host Microbe
3
:
213
–
223
.
Walker
AW
,
Duncan
SH
,
Harmsen
HJM
,
Holtrop
G
,
Welling
GW
,
Flint
HJ
. (
2008
).
The species composition of the human intestinal microbiota differs between particle-associated and liquid phase communities
.
Environ Microbiol
10
:
3275
–
3283
.
Wang
X
,
Conway
PL
,
Brown
IL
,
Evans
AJ
. (
1999
).
In vitro utilization of amylopectin and high amylose maize (amylomaize) starch granules by human colonic bacteria
.
Appl Environ Microbiol
65
:
4848
–
4854
.
Yue
JC
,
Clayton
MK
. (
2005
).
A similarity measure based on species proportions
.
Commun Stat Theo Methods
34
:
2123
–
2131
.
Zoetendal
EG
,
Akkermans
ADL
,
de Vos
WM
. (
1998
).
Temperature gradient gel electrophoresis analysis of 16S rRNA from human faecal samples reveals stable and host-specific communities of active bacteria
.
Appl Environ Microbiol
64
:
3854
–
3859
.
Zhang
H
,
DiBaise
JK
,
Zuccolo
A
,
Kudma
D
,
Braidotti
M
,
Yu
Y
et al. (
2008
).
Human gut microbiota in obesity and after gastric bypass
.
Proc Natl Acad Sci USA
106
:
2365
–
2370
.
Supplementary Information accompanies the paper on The ISME Journal website
Supplementary information The online version of this article (doi:10.1038/ismej.2010.118) contains supplementary material, which is available to authorized users.
© 2011 International Society for Microbial Ecology
Supplementary data
Citations
Views
Altmetric
Metrics
Total Views 1,028
751 Pageviews
277 PDF Downloads
Since 2/1/2024
Month: | Total Views: |
---|---|
February 2024 | 92 |
March 2024 | 2 |
April 2024 | 1 |
May 2024 | 149 |
June 2024 | 137 |
July 2024 | 154 |
August 2024 | 132 |
September 2024 | 135 |
October 2024 | 165 |
November 2024 | 61 |
×
Email alerts
Email alerts
Citing articles via
More from Oxford Academic