Email updates

Keep up to date with the latest news and content from BMC Gastroenterology and BioMed Central.

Open Access Research article

Microbial community analysis reveals high level phylogenetic alterations in the overall gastrointestinal microbiota of diarrhoea-predominant irritable bowel syndrome sufferers

Lotta Krogius-Kurikka1, Anna Lyra1, Erja Malinen1, Johannes Aarnikunnas1, Jarno Tuimala27, Lars Paulin3, Harri Mäkivuokko47, Kajsa Kajander56 and Airi Palva1*

Author Affiliations

1 Department of Basic Veterinary Sciences, Faculty of Veterinary Medicine, PO Box 66, FI-00014 University of Helsinki, Helsinki, Finland

2 CSC - Scientific Computing Ltd, Espoo, Finland

3 DNA Sequencing Laboratory, Institute of Biotechnology, University of Helsinki, Helsinki, Finland

4 Danisco Innovation, Kantvik, Finland

5 Valio Ltd, Research Centre, Helsinki, Finland

6 Department of Biomedicine, Faculty of Medicine, University of Helsinki, Helsinki, Finland

7 The Finnish Red Cross, Blood Service, Helsinki, Finland

For all author emails, please log on.

BMC Gastroenterology 2009, 9:95  doi:10.1186/1471-230X-9-95

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-230X/9/95


Received:10 July 2009
Accepted:17 December 2009
Published:17 December 2009

© 2009 Krogius-Kurikka et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

A growing amount of scientific evidence suggests that microbes are involved in the aetiology of irritable bowel syndrome (IBS), and the gastrointestinal (GI) microbiota of individuals suffering from diarrhoea-predominant IBS (IBS-D) is distinguishable from other IBS-subtypes. In our study, the GI microbiota of IBS-D patients was evaluated and compared with healthy controls (HC) by using a high-resolution sequencing method. The method allowed microbial community analysis on all levels of microbial genomic guanine plus cytosine (G+C) content, including high G+C bacteria.

Methods

The collective faecal microbiota composition of ten IBS-D patients was analysed by examining sequences obtained using percent G+C (%G+C) -based profiling and fractioning combined with 16S rRNA gene clone library sequencing of 3267 clones. The IBS-D library was compared with an analogous healthy-control library of 23 subjects. Real-time PCR analysis was used to identify phylotypes belonging to the class Gammaproteobacteria and the order Coriobacteriales.

Results

Significant differences were found between clone libraries of IBS-D patients and controls. The microbial communities of IBS-D patients were enriched in Proteobacteria and Firmicutes, but reduced in the number of Actinobacteria and Bacteroidetes compared to control. In particular, 16S rDNA sequences belonging to the family Lachnospiraceae within the phylum Firmicutes were in greater abundance in the IBS-D clone library.

Conclusions

In the microbiota of IBS-D sufferers, notable differences were detected among the prominent bacterial phyla (Firmicutes, Actinobacteria, Bacteroidetes, and Proteobacteria) localized within the GI tract.

Background

Irritable bowel syndrome (IBS), a common functional disorder of the gastrointestinal (GI) tract diagnosed in 10-20% of the adult and adolescent populations, is characterized by abdominal pain or discomfort, distorted bowel habits and altered stool characteristics [1]. Although IBS does not predispose to severe ill health, it diminishes the patients' quality of life and has an economic impact on society via work absenteeism and medical leave [2]. Because of the differing symptoms experienced by IBS-diagnosed individuals, sufferers have been divided into four subtypes, i.e. constipation-predominant IBS (IBS-C), diarrhoea-predominant IBS (IBS-D), mixed-type IBS (IBS-M) and unsubtyped IBS [1]. The aetiology of IBS is currently unknown, although genetic, environmental, psychosocial or physiological factors are likely to contribute to the disorder [3,4]. Intestinal bacteria may also play a role in the onset and maintenance of IBS given that previous studies have indicated disparities in the microbiota between IBS-D sufferers and healthy individuals or, to lesser extent, those individuals diagnosed with the other IBS subtypes [5,6]. Consequently, a detailed analysis of the IBS-D associated GI microbiota present is justified.

Several studies have strengthened the argument for a role of intestinal microbes in the causation of IBS. Acute bacterial gastroenteritis often leads to lingering gastrointestinal complaints in individuals, one-third of whom develop IBS [7]. Moreover, elevated levels of serum antibodies specific for bacterial flagellins, the dominant antigens associated with Crohn's disease (CD) [8], have been detected in post-infectious IBS (PI-IBS) sufferers [9]. Evidence for a low level of mucosal inflammation within the GI tract has also emerged for all subtypes of IBS [10,11]. Reports have indicated an increased level of serine protease activity, possibly originating from bacteria [12], in faecal samples recovered from IBS-D sufferers. Although not correlated with bowel movements or present in individuals with acute infectious diarrhoea [13], serine protease activity in mice is associated with increased mucosal colonic permeability and heightened visceral hypersensitivity [12], symptoms that also occur in humans with IBS-D [14,15].

In our study, the combined 16S rRNA gene composition from the faeces of ten individuals with IBS-D, in agreement with the Rome II criteria, was analysed by %G+C profiling and fractioned DNA sequencing, compared with a similarly produced sequence library of 23 healthy individuals [16] and then examined for selected phylotypes using real-time PCR. We observed significant difference between the clone libraries constructed from the IBS-D sufferers and healthy controls samples in the abundance of the four major phyla of bacteria (Firmicutes, Actinobacteria, Bacteroidetes and Proteobacteria) found in the GI tract. To the best of our knowledge, these sequencing data represent the first comprehensive study to sufficiently describe the GI microbiota associated with IBS-D.

Methods

Study subjects

Faecal samples used for %G+C fractioning and cloning were collected from ten donors (six females, four males), with an average age of 46.5 years that suffered from IBS according to the established Rome II criteria [17]. Real-time PCR assays were performed for two additional individuals diagnosed with IBS-D (one female and one male aged 30 and 27 years, respectively) and for 22 healthy control (HC) individuals that had been used for the construction of the unfractioned DNA sample previously by Krogius-Kurikka et al. [16]. IBS-D subjects were excluded from the study if they were pregnant, lactating or unable to cooperate, had consumed antibiotics during the past two months, had undergone major or complicated abdominal surgery, were suffering from organic GI disease, severe systematic disease or endometriosis, or had been diagnosed with senile dementia. Healthy individuals with lactose intolerance, celiac disease and regular GI tract symptoms were not used as controls. All participating IBS-D individuals had undergone clinical and endoscopic GI examinations or had had a barium enema within a year prior to starting the study. Lactose-intolerant IBS-D individuals consuming low-lactose or lactose-free diets as well as individuals medicated for IBS were allowed to participate. All test subjects have been studied previously for a comparison of the GI microbiota between IBS-related and healthy individuals [18-21], and represented the placebo group of a probiotic intervention study [22]. The IBS-D faecal samples were recovered prior to onset of the intervention period and approximately two weeks to 12 months after the clinical GI examination, and the intestinal microbiota contained therein was considered typical of IBS-D sufferers.

Ethics

The study protocol was approved by the human ethics committee at the joint authority for the District Hospital of Helsinki and the Uusimaa (HUS) region. All participants provided written informed consent and were allowed to withdraw from the study at any point.

Sequencing of the %G+C profiled sample

To avoid methodological distortions and to enable comparison of IBS-D subjects' microbiota with that of healthy controls', the faecal samples from IBS-D subjects were exposed to the same processes as described in detail by Krogius-Kurikka et al. [16]. In brief, the procedure included bacterial genomic DNA isolation with the method described by Apajalahti et al. [23], followed by %G+C profiling and fractioning of an equal amount of pooled DNA samples according to their %G+C -content using 5%G+C intervals [24] (Additional file 1). The 16S rRNA was amplified from each of the %G+C fractions 30-35, 35-40, 45-50, 50-55, 60-65, 65-70 and 70-75 separately, with a low number of PCR cycles using two primer pairs with a broad bacterial range [25,26]. The amplicons were pooled and clone libraries were constructed separately for each fraction with the QIAGEN® PCR Cloning plus Kit (Qiagen, Hilden, Germany). The sequencing of the 3' end of the 16S rRNA gene clones was performed with the primer pD' [27]. For templates that were hard to sequence, 1% (v/v) dimethyl sulfoxide was used in the sequencing reaction. These templates were mainly representatives of the phylum Actinobacteria.

Additional file 1. Percent guanine plus cytosine profile of intestinal microbial genomic DNA pooled from IBS-D (n = 10) and healthy (n = 23) subjects. The amount of DNA is indicated as relative abundance (%) and the area under the curve is used for calculating the proportional amount of DNA in the separate fractions. The red line indicates IBS-D and the blue line HC. Modified from Kassinen et al. [21].

Format: PDF Size: 14KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Sequence analysis

Sequences were checked manually and the primer sequences were removed with the Staden Package pregap4 version 1.5 and gap version 4.10 assembly programs [28]. Sequences present in more than one clone library were considered to be non-chimeric. Potential chimeras were also searched by manually browsing the ClustalW 1.83 sequence alignment [29] with Bio Edit version 7.0.5.3 [30]. The sequences from %G+C fractions 25-30, 40-45 and 55-60 (AM276372 to AM277303 [21], comprising one-third of the total sample DNA, were included in any further analyses to encompass the whole sample.

Determination of operational taxonomic units and diversity measurements

MAFFT 6.603b [31] available at the CSC-IT Center for Science (Espoo, Finland) was used for 16S rRNA gene sequence alignments. The MAFFT FFT-NS-2 alignment option was used to align the 16S rRNA gene clone sequences of the IBS-D and healthy control (HC) libraries separately and together. The sequences were cut from the Escherichia coli position 430 (totally conserved GTAAA), resulting in an alignment including the 16S rRNA variable regions V1 and V2. The alignments were visually inspected, but they were not edited manually to avoid subjectivity and to maintain reproducibility of the alignments. Distance matrices were created from the cut alignments with Phylip 3.66 Dnadist [32] using the F84 evolution model. The sequences were assigned into phylotypes (operational taxonomic units, OTUs) with DOTUR [33] by applying the furthest neighbor rule in which all sequences within an OTU fulfil the similarity criterion. The 98% cut-off for sequence similarity was used to delimit an OTU. The OTU representatives determined separately for IBS-D and HC libraries were used in the RDP library compare [34] and UniFrac [35] analyses. Common OTU representatives were used in the SONS analysis [36] and in constructing a phylogenetic tree of the family Lachnospiraceae. From each OTU with less than 95% similarity to any EMBL nucleotide sequence database entry a representative clone was sequenced to near full-length as described by Krogius-Kurikka et al. [16], but the sequence analyses were preformed with the 16S rRNA region covering approximately 450 bp from the 5'-end of the 16S rRNA gene.

The indices for diversity and richness estimates were calculated using DOTUR [33]. Simpson's 1/D and Shannon's indices take into account the number of species present and the abundance of each species, with the value of the indices increasing with greater diversity [37,38]. The Chao estimator for richness [39] considers singletons and doubletons as rare species and the ACE richness estimator considers OTUs represented with less than ten sequences as rare species [40]. Coverage of the clone libraries was calculated with the formula of Good [41], which takes into account the number of singletons in the clone library. The Fasta EMBL Environmental and EMBL Prokaryote database searches [42] and the Ribosomal Database Project II (RDP II) Classifier Tool [43] were used to affiliate the phylotypes.

Comparison of IBS-D and HC libraries

The 16S rRNA gene clone libraries of the IBS-D patients and healthy controls [16] were compared using SONS [36], RDP library compare [34] and UniFrac [35]. SONS calculates the observed fraction of sequences in shared OTUs in each library and the observed fraction of shared OTUs in each library. A distance matrix calculated based on the common MAFFT alignment of IBS-D and HC sequences, as described above in the Determination of operational taxonomic units and diversity measurements -section, was used in the SONS analysis. The RDP library compare online analysis tool was used to make microbial community comparisons and the UniFrac program was applied to compare the microbial communities in a phylogenetic context. In both analyses, a comparison between IBS-D and HC libraries was performed for separately determined OTUs, as described in the section Determination of operational taxonomic units and diversity measurements. The RDP analysis was also performed at the sequence level.

In constructing a phylogenetic tree for the UniFrac analysis, a representative sequence of each OTU was aligned with MAFFT [31] using the E-INS-i alignment algorithm. Thereafter, 16S rRNA reference sequences from the European ribosomal RNA database [44] (Additional file 2) were selected and similarly aligned. The MAFFT-profile alignment option was used for constructing a combined profile alignment from the above mentioned alignments. The alignment was cut from E. coli position 430 (totally conserved GTAAA) and reference sequences, except for Methanobrevibacter smithii, were then deleted from the alignment with BioEdit version 7.0.5.3 [30]. A F84-corrected distance matrix was created using Phylip 3.66 dnadist [32]. The OTU representatives in the tree were labelled with taxonomic information from the Ribosomal Database Project II Classifier Tool [43] to identify the sequence affiliations. UniFrac Significance using abundance weights and P Test Significance analyses were used to describe whether the communities were significantly different overall. Lineage-Specific Analysis was used to test whether significant differences were present between the libraries within separate lineages at a specified distance from the root.

Additional file 2. RDP reference sequences. The RDP reference sequences [44] used in the profile alignments for UniFrac analysis [35] and in construction of the phylogenetic tree for the family Lachnospiraceae. Roman numerals indicate Clostridium rRNA clusters.

Format: PDF Size: 8KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Phylogenetic tree of the family Lachnospiraceae

A phylogenetic tree was constructed for the common OTU representatives from the IBS-D and HC libraries, with over 50% confidence threshold to the family Lachnospiraceae (Firmicutes) according to the RDP Classifier [43] (Additional file 3). The sequences were aligned with the reference sequences representing the Clostridium rRNA XIV group [45] and selected additional reference sequences (Additional files 2 and 3). An MAFFT profile alignment and an F84 distance matrix were constructed as for the UniFrac analysis in the Comparison of IBS-D and HC libraries -section, with the exception that sequences from the European ribosomal RNA database representing Clostridium rRNA Cluster XIV and Clostridium leptum AF262239 were used as references. Bootstrapping the data with one hundred replicates and construction of the tree were performed with the Seqboot and Consense programs of Phylip 3.66 [32]. A phylogenetic tree was generated with a neighbor-joining algorithm from the distance matrix using Phylip 3.66 neighbor [32]. The tree was visualized with MEGA4 [46].

Additional file 3. Phylogenetic tree of the family Lachnospiraceae. A neighbor-joining tree containing 201 common Lachnospiraceae OTUs for IBS-D and HC libraries. The number of sequences within an OTU is denoted after the abbreviation IBS-D or HC. Reference sequences for real-time PCR analyses from the studies by Kassinen et al. [21] and Lyra et al. [5] and the sequence for bacterium A4 (DQ789118) associated with CD [51] are denoted with red and blue font, respectively. Reference sequences presenting the Clostridium rRNA XIV group are denoted with green font. Bootstrap values are percentages of 100 resamplings and the scale bar represents 0.06 substitutions per nucleotide position.

Format: PDF Size: 47KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Real-time PCR assays

The real-time PCR assays for Enterobacteriaceae and Eggerthella lenta -like phylotypes were developed since the common OTUs for HC and IBS-D, presenting these phylotypes, were notably more abundant in sequences from the IBS-D library (Figure 1) and the genera were significantly more abundant in the IBS-D library according to the RDP library compare tool (Additional file 4). The primers and assays were designed as described in Lyra et al. [5]. The forward and reverse primers and positive control clones for real-time PCR assays were 5'-CATAACGTCGCAAGACCAAAGA-3', 5'-GAGTCTGGACCGTGTCTCAGTTC-3' and AM276420 for the first and 5'-GTGACCAACCTGCCCCTTG-3', 5'-GACCCCATCCCTTGCCGT-3' and AM276078 for the latter. The clones AM275716 and AM693283 were used as negative controls for the assays targeting Enterobacteriaceae and Eggerthella lenta -like phylotypes, respectively. The negative controls were highly similar in sequence with the positive target clones, but included mismatches within the primer annealing sites.

Additional file 4. RDP library compare results for sequences and OTUs. Significantly differing (p-values < 0.01) groups of sequences and OTUs in the IBS-D and HC libraries and their phylogenetic affiliation according to RDP library compare [34]. The more abundant group is indicated in boldface.

Format: PDF Size: 6KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 1. Distribution of sequences and OTUs among IBS-D and HC clone libraries. The unique sequences or OTUs (cut-off-level of 98%) in the IBS-D and HC libraries are indicated in red and blue, respectively. The grey area in the group-wise bars is a mirror-image of shared sequences or OTUs and it is presented on both sides of the y-axis. The number of singleton OTUs is given in parentheses. The roman numerals XIV and IV within the Firmicutes indicate the corresponding Clostridium rRNA clusters, 1) Other Firmicutes, 2) Acidobacteria, Cyanobacteria, TM7 and Verrucomicrobia.

The analyses were performed using the iCycler iQ Real-Time Detection System (Bio-Rad, Hercules, CA, USA) associated with the iCycler Optical System Interface software (version 2.3; Bio-Rad). The 12 IBS-D and 22 HC patients' individual faecal DNA samples were run as triplicates using the following optimized reaction conditions: 25 ng of faecal DNA, 1:75 000 dilution of SYBR Green I (Lonza Biosciences, Basel, Switzerland), 10 mM Tris-HCl (pH 8.8), 50 mM KCl, 0.1% Triton X-100, 2 mM and 3 mM of MgCl2 for Enterobacteriaceae and E. lenta (respectively), 100 μM each dNTP, 0.5 μM each primer, 0.024 U Dynazyme II polymerase (Finnzymes, Espoo, Finland) and 5 μl of either template or water. Standards ranging from 102 to 107 16S rRNA gene copies, amplified from the positive controls, were applied. After an initial denaturation at 95°C for 5 min the real-time PCR amplification proceeded with 40 cycles of denaturation at 95°C for 20 s, primer annealing at 64°C for Enterobacteriaceae and at 68°C for E. lenta for 20 s, extension at 72°C for 30 s and a fluorescence detection step at 85°C for 30 s. A melt curve analysis was preformed after amplification by slow cooling from 95°C to 60°C, with fluorescence collection at 0.3°C intervals and a hold of 10 s at each decrement to check the specificity of the real-time PCR assay. The raw data were transformed to log10 ratios of the number of 16S rRNA gene copies in one gram of faeces. The R software environment for statistical computing and graphics [47] was used for performing non-parametric Mann-Whitney U-tests.

Nucleotide sequence accession numbers

The 16S rRNA gene sequences reported in this study were deposited in the EMBL Nucleotide Sequence Database under accession numbers AM691850 to AM694184.

Results

Library characteristics

Approximately 300 sequences were recovered from each %G+C fraction of the IBS-D patients' (n = 10) pooled sample, resulting in a total of 3267 sequences. The IBS-D clone library constructions and sequencing were performed similarly to those of the healthy controls (n = 23) that comprised of 3199 sequences [16] with the exception of amount of subjects pooled to construct the sample (Table 1). According to Good's formula [41], the coverage of clone libraries was above 95%. Less phylotypes were present in the IBS-D library (n = 302) than in the HC library (n = 428), with a 98% cut-off level for OTUs [16] (Table 1). The Shannon and Simpson indices for diversity and Chao and ACE richness estimates were lower for the IBS-D library. The rank abundance curves of the libraries showed highly similar OTU evenness (Figure 2).

thumbnailFigure 2. Rank-abundance plots of the IBS-D and HC libraries. The curves for IBS-D and HC libraries are indicated in red and blue, respectively. To make the image more compact, the OTUs with the highest number of sequences are deleted (an OTU of 395 sequences for IBS-D and an OTU of 160 sequences for HC).

Table 1. Characteristics of the IBS-D and HC libraries.

Microbial community comparisons

The microbial community comparison at the phylum level using RDP classifier revealed that the IBS-D library had significantly more representatives of Proteobacteria and Firmicutes than the HC library, whereas the condition was the opposite with Actinobacteria and Bacteroidetes (Figure 3). When the comparison was made with OTUs, the IBS-D library was significantly richer in Firmicutes than the HC library (Figure 3). In the UniFrac Lineage-Specific Analysis, the phylum Actinobacteria differed significantly (p = 0.0013), and the phylum Bacteroidetes, the Firmicutes families Lachnospiraceae and Ruminococcaeae and the Proteobacteria classes Gammaproteobacteria and Alphaproteobacteria differed highly significantly (p < 0.001) between the IBS-D and HC libraries. Overall, however, the libraries did not differ significantly.

thumbnailFigure 3. Relative abundance of phyla in the IBS-D and HC libraries in the RDP library compare analysis [34]. Significantly differing (p < 0.01) abundances of sequences (s) or OTUs with a cut-off-level of 98% for sequence similarity (o) between the IBS-D and HC libraries are indicated. The "unclassified" phyla have a bootstrap value below 80%.

A substantial proportion of sequences were members of the family Lachnospiraceae in both libraries; 45% and 33% in the IBS-D and HC libraries, respectively (Additional file 3). The Lachnospiraceae diverged significantly between IBS-D and HC according to the RDP library compare sequence and OTU-based analyses (Additional file 4). Among significantly differing groups of Actinobacteria, Eggerthella was the only genus more abundant in the IBS-D library than in the HC library (Additional file 4). The IBS-D library contained significantly less Bacteroidetes sequences than the HC library (29 vs. 96) in both the RDP library compare and UniFrac analyses. Presence of GI pathogens was not demonstrated by sequencing.

The combined sequence pool from both community samples comprised 578 OTUs, 30.4% of which were shared (Figure 4). The IBS-D library had only half the number of unique OTUs found in the HC library. The majority of the sequences were, however, shared (81.0%), and the proportions of unique sequences in both libraries were similar. The proportion of shared OTUs in the IBS-D library was greater than in the HC library at all OTU cut-off levels, with the highest difference being observed at the (cut-off) level of 90% (Figure 5). The HC library harboured more unique sequences and OTUs than the IBS-D library, with a few exceptions; the portion of unique sequences in the classes Bacilli and Proteobacteria, and the portion of unique OTUs affiliating with the Clostridium XIV cluster and the order Bifidobacteriales were more abundant in the IBS-D library (Figure 1).

thumbnailFigure 4. Percentage of shared and unique OTUs (n = 578) and sequences (n = 6466) in the combined IBS-D and HC libraries. Number of sequences and OTUs (cut-off-level of 98%) are given in parentheses. The number of singletons is indicated after a semicolon. A common alignment was used for the determination of shared OTUs.

thumbnailFigure 5. Abundance of shared OTUs. Proportion of shared OTUs in the IBS-D and HC libraries according to the SONS analysis [36].

Real-time PCR analysis

The average PCR efficiencies were 89% for the Enterobacteriaceae and 96% for E. lenta -like real-time PCR assays. The differences between the IBS-D and HC libraries were not significant according to the Mann-Whitney U-tests (Figure 6). The result remained similar when the amounts of detected phylotypes were proportioned to the total amount of bacteria; Mann-Whitney U-test p-values 0.15 and 0.63 for Enterobacteriaceae and Eggerthella lenta -like phylotypes, respectively (data not shown). According to the dry weights, the moisture content of the samples did not affect the real-time PCR results (data not shown).

thumbnailFigure 6. Real-time PCR results from assays for Enterobacteriaceae (A) and Eggerthella lenta -like (B) phylotypes. Samples from IBS-D and HC subjects are in red and blue, respectively. The values are 16S rRNA gene copy numbers per gram of faeces (log10 values), and the detection limit is set to 104. Vertical lines are medians from the Mann-Whitney U-test.

Discussion

To determine whether diarrhoea-predominant IBS is linked to particular changes in the GI microbiota, 16S rRNA gene sequence data constructed from pooled faecal DNA samples of ten IBS patients with diarrhoea were compared with the similarly constructed data of 23 healthy controls [16]. The %G+C fractioning enhances the distinction of a wide spectrum of bacterial phylotypes during the subsequent 16S rRNA gene library construction and sequence analysis. The same extent of cloning was applied to %G+C fractions with varying amount of DNA, both within a sample and between samples, possibly skewing the relative abundance of OTUs detected within the libraries. However, the method allowed thorough comparison of disturbed (IBS-D) and healthy GI microbiota also in the high G+C region.

Lower number of OTUs in the IBS-D library

The lower diversity and number of OTUs in the IBS-D library is probably mostly due to the smaller number of subjects used to construct the library. On the other hand, a decreased diversity in the GI microbiota of IBS-D patients compared with healthy individuals may exist, as has been observed with CD patients [48]. Acute diarrhoea has been shown to reduce the overall microbiota composition, which could simply be caused by washing out of the commensal bacteria [49]. At the approximate family/genus level (90% OTU cut-off), the number of shared OTUs was above 85% in the IBS-D library compared with 50% in the HC library. Thus the lower diversity of IBS-D patients' GI microbiota seems more apparent at a higher taxonomic level, possibly indicating dysbiosis. A likely cause for these observations is the subject number difference in the pooled samples. On the other hand, the Simpson diversity indexes were very similar below approximately 90% OTU cut-off levels, suggesting little difference between the two libraries (data not shown). However, comparative analysis at higher taxonomic level is likely to be less affected by individuality and thus less prone to such bias.

Abundance of Lachnospiraceae in IBS-D library

Some data suggest an overlap between the aetiology of inflammatory bowel disease (IBD) and IBS [50]. Duck and colleagues [51] isolated a bacterium (A4) with a highly similar flagellin (A4-Fla2) to Fla-X, which has been shown to be an important antigen in CD [8]. The A4 bacterium was classified as a member of the family Lachnospiraceae and the Clostridium cluster XIVa. Interestingly, elevated antibody concentrations towards the flagellins A4-Fla2 and Fla-X have also been associated with IBS, especially with the PI-IBS subgroup of patients [9], which is most often diarrhoea-predominant [52]. In our study, the family Lachnospiraceae was significantly more abundant in the IBS-D library than in the HC library. A large part of the sequences and OTUs in IBS-D (45% and 41%) and HC (33% and 30%) libraries are affiliated with the family Lachnospiraceae, including the largest OTU in the IBS-D (n = 395) and HC (n = 160) libraries, which is affiliated with Eubacterium rectale, a prevalent member of this family in the gut. Contradictory results on the abundance of Clostridium coccoides -Eubacterium rectale group bacteria among IBS patients have been obtained in previous studies, which may be due to the different analysis methods used and the broadness of the group [18,20].

Rajilić-Stojanović et al. [6] observed a significantly lower level of Bacteroides spp. in IBS patients than in healthy controls, but at the Bacteroidetes phylum level no difference was detected between the groups. In this study, the number of Bacteroidetes phylum sequences was lower in the IBS-D library. However, the number of sequences assigned to this phylum was overall low, and thus, extra caution should be employed in interpreting the data. The group Bacteroides-Prevotella-Porphyromonas has previously been quantified from the faeces of healthy controls and IBS patients, with no significant difference between the groups [18]. The maximum four hour delay due to transportation from the study subjects to the laboratory in the freezing of the samples in -70°C may have lowered the proportion of Bacteroides detected (Salonen et al., personal communications). All samples were, however, stored and processed in the same manner to minimize any technical bias into the comparative analysis.

Real-time PCR quantification compared with community comparison

The sequencing results were in accordance with the earlier real-time PCR quantification of a Collinsella aerofaciens -like phylotype [21], Bifidobacterium spp. and Lactobacillus spp. [18], all being less abundant in the IBS-D library than in the HC library. Several OTUs comprising at least 100 sequences were not commonly found in the IBS-D library. These OTUs were affiliated with Bifidobacterium longum (Bifidobacteriaceae), Bifidobacterium pseudocatenulatum (Bifidobacteriaceae), Eggerthella lenta (Coriobacteriaceae), a phylotype with 93% similarity to E. lenta (Coriobacteriaceae), Enterobacter ludwigii (Gammaproteobacteria) and Streptococcus bovis (Bacillus). However, real-time PCR assays performed on individual samples for targeting these OTUs, e.g. B. longum [18], a B. catenulatum/B. pseudocatenulatum-like phylotype [21] and a S. bovis-like phylotype [21], have not revealed significant differences between IBS subjects and healthy controls. Additionally, phylotypes affiliating with Lachnospiraceae resembling Ruminococcus torques with less than 95% 16S rRNA gene sequence similarity have been detected with an association to either IBS-D or healthy depending on the phylotype [5]. Caution should be taken when comparing the real-time PCR results and the sequencing data since the targets of the real-time PCR assays and the taxonomic grouping of the sequence data may not be congruent.

In this study, an Eggerthella lenta -like phylotype and an Enterobacteriaceae phylotype were quantified with real-time PCR. The median of the 16S rRNA gene quantities for the E. lenta -like phylotype was lower for IBS- than for HC (Figure 6), similarly as previously reported for the Atopobíum group, which contains the E. lenta -like pylotype [18]. However, these differences were not significant. A minor tendency for Enterobacteriaceae to be more abundant among IBS-D patients than controls was detected. Consistent with this, elevated number of Enterobacteria have previously been detected among CD patients [53]. Evidence of higher quantities of aerobic bacteria in the clone library of IBS-D subjects was seen among Gammaproteobacteria and Bacilli. An elevated aerobe:anaerobe ratio [19] and higher amounts of Enterobacteriaceae [54] have earlier been detected in association with IBS using culture-based techniques. Furthermore, in association with acute diarrhoea, increases in the number of enteric (aerobic) bacteria have been reported [49,55], while anaerobic bacterial counts have decreased [56].

Conclusions

Based on the 16S rRNA gene sequencing approach applied here, the faecal sample of IBS-D patients showed indications of a dysbiotic microbiota, even though the overall structure was similar to that of healthy controls. Notable differences between IBS-D and controls were identified in all dominant bacterial phyla of the GI microbiota; Firmicutes, Actinobacteria, Bacteroidetes and Proteobacteria. Differences on such a high taxonomic level are not as prone to biases in sequence library comparisons as lower level comparisons. Within Firmicutes, the family Lachnospiraceae was significantly increased in the IBS-D group regarding the number of sequences and OTUs. In future research, the role of this family in IBS-D should receive more attention.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

LK-K and JA edited the sequence data. LK-K and AL conducted the community analyses and together with EM prepared the manuscript. JT acted as a bioinformatics consultant, LP supervised the sequencing process and HM performed the %G+C profiling and fractioning. KK recruited the IBS-D subjects and planned and coordinated the collection of samples. AP coordinated and supervised the study. All authors made corrections and approved the final manuscript.

Acknowledgements

This study was supported by the Finnish Funding Agency for Technology and Innovation (Grant no. 40160/05), the Academy of Finland (Grant no. 214 157) and the Finnish Graduate School on Applied Bioscience. This work was performed at the Centre of Excellence on Microbial Food Safety Research, Academy of Finland. We are grateful to Sinikka Ahonen, Anu Suoranta, Matias Rantanen, Laura Mäkelä and Annemari Wickström for technical assistance and Professor Jukka Corander and Doctor Janne Nikkilä for statistical consultation. Doctors Maria Saarela and Jaana Mättö are gratefully acknowledged for recruiting HC study subjects and management of sample collection. Doctor Ingemar von Ossowski is kindly acknowledged for language revision. Kyösti Kurikka, MSc, is thanked for drawing the figures.

References

  1. Longstreth GF, Thompson WG, Chey WD, Houghton LA, Mearin F, Spiller RC: Functional bowel disorders.

    Gastroenterology 2006, 130(5):1480-1491. PubMed Abstract | Publisher Full Text OpenURL

  2. Talley NJ: Functional gastrointestinal disorders as a public health problem.

    Neurogastroenterol Motil 2008, 20(Suppl 1):121-129. PubMed Abstract | Publisher Full Text OpenURL

  3. Öhman L, Simrén M: New insights into the pathogenesis and pathophysiology of irritable bowel syndrome.

    Dig Liver Dis 2007, 39(3):201-215. PubMed Abstract | Publisher Full Text OpenURL

  4. Drossman DA: The functional gastrointestinal disorders and the Rome III process.

    Gastroenterology 2006, 130(5):1377-1390. PubMed Abstract | Publisher Full Text OpenURL

  5. Lyra A, Rinttilä T, Nikkilä J, Krogius-Kurikka L, Kajander K, Malinen E, Mättö J, Mäkelä L, Palva A: Diarrhoea-predominant irritable bowel syndrome distinguishable by 16S rRNA gene phylotype quantification.

    World J Gastroenterol 15(47):5936-45. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Rajilić-Stojanović M: Diversity of the Human Gastrointestinal Microbiota - Novel Perspectives from High Troughput Analyses. University of Wageningen, Wageningen, The Netherlands; 2007.

  7. Dupont AW: Post-infectious irritable bowel syndrome.

    Curr Gastroenterol Rep 2007, 9(5):378-384. PubMed Abstract | Publisher Full Text OpenURL

  8. Lodes MJ, Cong Y, Elson CO, Mohamath R, Landers CJ, Targan SR, Fort M, Hershberg RM: Bacterial flagellin is a dominant antigen in Crohn disease.

    J Clin Invest 2004, 113(9):1296-1306. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Schoepfer AM, Schaffer T, Seibold-Schmid B, Muller S, Seibold F: Antibodies to flagellin indicate reactivity to bacterial antigens in IBS patients.

    Neurogastroenterol Motil 2008, 20(10):1110-1118. PubMed Abstract | Publisher Full Text OpenURL

  10. Quigley EM: Changing face of irritable bowel syndrome.

    World J Gastroenterol 2006, 12(1):1-5. PubMed Abstract | Publisher Full Text OpenURL

  11. Parkes GC, Brostoff J, Whelan K, Sanderson JD: Gastrointestinal microbiota in irritable bowel syndrome: their role in its pathogenesis and treatment.

    Am J Gastroenterol 2008, 103(6):1557-1567. PubMed Abstract | Publisher Full Text OpenURL

  12. Gecse K, Róka R, Ferrier L, Leveque M, Eutamene H, Cartier C, Ait-Belgnaoui A, Rosztóczy A, Izbéki F, Fioramonti J, Wittmann T, Bueno L: Increased faecal serine protease activity in diarrhoeic IBS patients: a colonic lumenal factor impairing colonic permeability and sensitivity.

    Gut 2008, 57(5):591-599. PubMed Abstract | Publisher Full Text OpenURL

  13. Róka R, Rosztóczy A, Leveque M, Izbéki F, Nagy F, Molnár T, Lonovics J, Garcia-Villar R, Fioramonti J, Wittmann T, Bueno L: A pilot study of fecal serine-protease activity: a pathophysiologic factor in diarrhea-predominant irritable bowel syndrome.

    Clin Gastroenterol Hepatol 2007, 5(5):550-555. PubMed Abstract | Publisher Full Text OpenURL

  14. Dunlop SP, Hebden J, Campbell E, Naesdal J, Olbe L, Perkins AC, Spiller RC: Abnormal intestinal permeability in subgroups of diarrhea-predominant irritable bowel syndromes.

    Am J Gastroenterol 2006, 101(6):1288-1294. PubMed Abstract | Publisher Full Text OpenURL

  15. Annaházi A, Gecse K, Dabek M, Ait-Belgnaoui A, Rosztóczy A, Róka R, Molnár T, Theodorou V, Wittmann T, Bueno L, Eutamene H: Fecal proteases from diarrheic-IBS and ulcerative colitis patients exert opposite effect on visceral sensitivity in mice.

    Pain 2009, 144(1-2>):209-217. PubMed Abstract | Publisher Full Text OpenURL

  16. Krogius-Kurikka L, Kassinen A, Paulin L, Corander J, Mäkivuokko H, Tuimala J, Palva A: Sequence Analysis of Percent G+C Fraction Libraries of Human Faecal Bacterial DNA Reveals a High Number of Actinobacteria.

    BMC Microbiol 2009, (9):68. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  17. Thompson WG, Longstreth GF, Drossman DA, Heaton KW, Irvine EJ, Muller-Lissner SA: Functional bowel disorders and functional abdominal pain.

    Gut 1999, 45(Suppl 2):II43-7. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Malinen E, Rinttilä T, Kajander K, Mättö J, Kassinen A, Krogius L, Saarela M, Korpela R, Palva A: Analysis of the fecal microbiota of irritable bowel syndrome patients and healthy controls with real-time PCR.

    Am J Gastroenterol 2005, 100(2):373-382. PubMed Abstract | Publisher Full Text OpenURL

  19. Mättö J, Maunuksela L, Kajander K, Palva A, Korpela R, Kassinen A, Saarela M: Composition and temporal stability of gastrointestinal microbiota in irritable bowel syndrome--a longitudinal study in IBS and control subjects.

    FEMS Immunol Med Microbiol 2005, 43(2):213-222. PubMed Abstract | Publisher Full Text OpenURL

  20. Maukonen J, Satokari R, Mättö J, Söderlund H, Mattila-Sandholm T, Saarela M: Prevalence and temporal stability of selected clostridial groups in irritable bowel syndrome in relation to predominant faecal bacteria.

    J Med Microbiol 2006, 55(Pt 5):625-633. PubMed Abstract | Publisher Full Text OpenURL

  21. Kassinen A, Krogius-Kurikka L, Mäkivuokko H, Rinttilä T, Paulin L, Corander J, Malinen E, Apajalahti J, Palva A: The fecal microbiota of irritable bowel syndrome patients differs significantly from that of healthy subjects.

    Gastroenterology 2007, 133(1):24-33. PubMed Abstract | Publisher Full Text OpenURL

  22. Kajander K, Krogius-Kurikka L, Rinttilä T, Karjalainen H, Palva A, Korpela R: Effects of multispecies probiotic supplementation on intestinal microbiota in irritable bowel syndrome.

    Aliment Pharmacol Ther 2007, 26(3):463-473. PubMed Abstract | Publisher Full Text OpenURL

  23. Apajalahti JH, Särkilahti LK, Mäki BR, Heikkinen JP, Nurminen PH, Holben WE: Effective recovery of bacterial DNA and percent-guanine-plus-cytosine-based analysis of community structure in the gastrointestinal tract of broiler chickens.

    Appl Environ Microbiol 1998, 64(10):4084-4088. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Holben WE, Harris D: DNA-based monitoring of total bacterial community structure in environmental samples.

    Mol Ecol 1995, 4(5):627-631. PubMed Abstract | Publisher Full Text OpenURL

  25. Hicks RE, Amann RI, Stahl DA: Dual staining of natural bacterioplankton with 4',6-diamidino-2-phenylindole and fluorescent oligonucleotide probes targeting kingdom-level 16S rRNA sequences.

    Appl Environ Microbiol 1992, 58(7):2158-2163. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Kane MD, Poulsen LK, Stahl DA: Monitoring the enrichment and isolation of sulfate-reducing bacteria by using oligonucleotide hybridization probes designed from environmentally derived 16S rRNA sequences.

    Appl Environ Microbiol 1993, 59(3):682-686. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Edwards U, Rogall T, Blöcker H, Emde M, Böttger EC: Isolation and direct complete nucleotide determination of entire genes. Characterization of a gene coding for 16S ribosomal RNA.

    Nucleic Acids Res 1989, 17(19):7843-7853. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Staden R, Beal KF, Bonfield JK: The Staden package.

    Methods Mol Biol 2000, 132:115-130. PubMed Abstract OpenURL

  29. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice.

    Nucleic Acids Res 1994, 22(22):4673-4680. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Hall T: BioEdit. Biological sequence alignment editor for Windows. [http://www.mbio.ncsu.edu/BioEdit/bioedit.html] webcite

    North Carolina State University, NC, USA; 1998.

  31. Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program.

    Brief Bioinform 2008, 9(4):286-298. PubMed Abstract | Publisher Full Text OpenURL

  32. Felsenstein J: PHYLIP - Phylogeny Inference package (Version 3.2).

    Cladistics 1989, (17):164-166. OpenURL

  33. Schloss PD, Handelsman J: Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness.

    Appl Environ Microbiol 2005, 71(3):1501-1506. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ, Kulam-Syed-Mohideen AS, McGarrell DM, Marsh T, Garrity GM, Tiedje JM: The Ribosomal Database Project: improved alignments and new tools for rRNA analysis.

    Nucleic Acids Res 2009, (37 Database):D141-5. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Lozupone C, Hamady M, Knight R: UniFrac--an online tool for comparing microbial community diversity in a phylogenetic context.

    BMC Bioinformatics 2006, 7:371. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  36. Schloss PD, Handelsman J: Introducing SONS, a tool for operational taxonomic unit-based comparisons of microbial community memberships and structures.

    Appl Environ Microbiol 2006, 72(10):6773-6779. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Simpson EH: Measurement of diversity.

    Nature 1949, 163:688. Publisher Full Text OpenURL

  38. Shannon CE: A mathematical theory of communication.

    Bell Syst Tech J 1948., 27(379-423 and 623-656) OpenURL

  39. Chao A: Nonparametric estimation of the number of classes in a population.

    Scand J Statist 1984, (11):265-270. OpenURL

  40. Chao A, Lee S: Estimating the number of classes via sample coverage.

    J Amer Stat Assoc 1992, 87(417):210-217. Publisher Full Text OpenURL

  41. Good IJ: The population frequencies of species and the estimation of population parameters.

    Biometrika 1953, 40:237-264. OpenURL

  42. Pearson WR: Rapid and sensitive sequence comparison with FASTP and FASTA.

    Methods Enzymol 1990, 183:63-98. PubMed Abstract | Publisher Full Text OpenURL

  43. Wang Q, Garrity GM, Tiedje JM, Cole JR: Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy.

    Appl Environ Microbiol 2007, 73(16):5261-5267. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Wuyts J, Perriere G, Peer Y: The European ribosomal RNA database.

    Nucleic Acids Res 2004, (32 Database):D101-3. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Collins MD, Lawson PA, Willems A, Cordoba JJ, Fernandez-Garayzabal J, Garcia P, Cai J, Hippe H, Farrow JA: The phylogeny of the genus Clostridium: proposal of five new genera and eleven new species combinations.

    Int J Syst Bacteriol 1994, 44(4):812-826. PubMed Abstract | Publisher Full Text OpenURL

  46. Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment.

    Brief Bioinform 2004, 5(2):150-163. PubMed Abstract | Publisher Full Text OpenURL

  47. R Development Core Team: R: A Language and Environment for Statistical Computing. [http://www.r-project.org] webcite

    Vienna, Austria 2007. OpenURL

  48. Manichanh C, Rigottier-Gois L, Bonnaud E, Gloux K, Pelletier E, Frangeul L, Nalin R, Jarrin C, Chardon P, Marteau P, Roca J, Doré J: Reduced diversity of faecal microbiota in Crohn's disease revealed by a metagenomic approach.

    Gut 2006, 55(2):205-211. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  49. Mai V, Braden CR, Heckendorf J, Pironis B, Hirshon JM: Monitoring of stool microbiota in subjects with diarrhea indicates distortions in composition.

    J Clin Microbiol 2006, 44(12):4550-4552. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  50. Quigley EM: Irritable bowel syndrome and inflammatory bowel disease: interrelated diseases?

    Chin J Dig Dis 2005, 6(3):122-132. PubMed Abstract | Publisher Full Text OpenURL

  51. Duck LW, Walter MR, Novak J, Kelly D, Tomasi M, Cong Y, Elson CO: Isolation of flagellated bacteria implicated in Crohn's disease.

    Inflamm Bowel Dis 2007, 13(10):1191-1201. PubMed Abstract | Publisher Full Text OpenURL

  52. Dunlop SP, Jenkins D, Spiller RC: Distinctive clinical, psychological, and histological features of postinfective irritable bowel syndrome.

    Am J Gastroenterol 2003, 98(7):1578-1583. PubMed Abstract | Publisher Full Text OpenURL

  53. Seksik P, Rigottier-Gois L, Gramet G, Sutren M, Pochart P, Marteau P, Jian R, Doré J: Alterations of the dominant faecal bacterial groups in patients with Crohn's disease of the colon.

    Gut 2003, 52(2):237-242. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Si JM, Yu YC, Fan YJ, Chen SJ: Intestinal microecology and quality of life in irritable bowel syndrome patients.

    World J Gastroenterol 2004, 10(12):1802-1805. PubMed Abstract | Publisher Full Text OpenURL

  55. Young VB, Schmidt TM: Antibiotic-associated diarrhea accompanied by large-scale alterations in the composition of the fecal microbiota.

    J Clin Microbiol 2004, 42(3):1203-1206. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  56. Barman M, Unold D, Shifley K, Amir E, Hung K, Bos N, Salzman N: Enteric salmonellosis disrupts the microbial ecology of the murine gastrointestinal tract.

    Infect Immun 2008, 76(3):907-915. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-230X/9/95/prepub