|Home | About | Journals | Submit | Contact Us | Français|
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.
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.
Volunteers were initially recruited by newspaper advertisement. Inclusion criteria for men were a body mass index >27kgm−2, waist circumference >102cm and fasting plasma glucose concentration >6.0mmoll−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 (kgm−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.
The volunteers were weight stable (<3kg 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.7g 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.5MJkg−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 24h collection from the final week of each dietary period was used for chemical analysis of digestibility.
Diets were analysed for total gross energy, RS and insoluble and soluble NSPs (Supplementary Table S2). Estimation of RS content of diets and 24h faecal sample collections was as described by Englyst et al. (1992) whereas NSPs were quantified as described by Englyst and Cummings (1988).
Faecal samples were kept at 4°C and processed within 5h 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 5g dispersed (3500r.p.m. for 1min using a Dispomix Drive, Medic Tools, Lussiwag, Switzerland) in 10ml 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.
PCR amplification of 16S rRNA genes from the extracted DNA involved initial denaturation at 94°C for 5min; 20 cycles of denaturation (30s at 95°C), annealing (30s at 52°C), extension at 72°C for 2min with a final cycle at 72°C for 8min. Taq1 polymerase from Promega (Southampton, UK) was used at 0.025Uμl−1 in the presence of 2.5m 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).
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.5m MgCl2 and 0.025U Taq1 polymerase per μl. Products were separated in an 8% acrylamide denaturing gradient gel containing urea and formamide (35–60% gradient) at 80V for 16h at 60°C and visualized by staining with SybrGold dye (Invitrogen, Paisley, UK).
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.
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).
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).
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).
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).
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 and5).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).
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).
Whole tract digestibility of RS and NSP was estimated from chemical analysis of the diet and of 24h 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 and5).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).
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%).
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 60h (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.
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).
Supplementary Information accompanies the paper on The ISME Journal website (http://www.nature.com/ismej)