Alterations of microbiota in urine from women with interstitial cystitis
© Siddiqui et al.; licensee BioMed Central Ltd. 2012
Received: 25 June 2012
Accepted: 5 September 2012
Published: 13 September 2012
Skip to main content
© Siddiqui et al.; licensee BioMed Central Ltd. 2012
Received: 25 June 2012
Accepted: 5 September 2012
Published: 13 September 2012
Interstitial Cystitis (IC) is a chronic inflammatory condition of the bladder with unknown etiology. The aim of this study was to characterize the microbial community present in the urine from IC female patients by 454 high throughput sequencing of the 16S variable regions V1V2 and V6. The taxonomical composition, richness and diversity of the IC microbiota were determined and compared to the microbial profile of asymptomatic healthy female (HF) urine.
The composition and distribution of bacterial sequences differed between the urine microbiota of IC patients and HFs. Reduced sequence richness and diversity were found in IC patient urine, and a significant difference in the community structure of IC urine in relation to HF urine was observed. More than 90% of the IC sequence reads were identified as belonging to the bacterial genus Lactobacillus, a marked increase compared to 60% in HF urine.
The 16S rDNA sequence data demonstrates a shift in the composition of the bacterial community in IC urine. The reduced microbial diversity and richness is accompanied by a higher abundance of the bacterial genus Lactobacillus, compared to HF urine. This study demonstrates that high throughput sequencing analysis of urine microbiota in IC patients is a powerful tool towards a better understanding of this enigmatic disease.
Interstitial Cystitis or Painful Bladder Syndrome (IC/PBS) is a chronic condition characterized by frequent urination and bladder pain, which often results in reduced quality of life. Clinicians experience that this disease is becoming more prevalent . While evidence suggests that about 90% of those affected are female, some urologists consider chronic bacterial prostatitis to be the male equivalent of IC.
The traditional definition of this illness is chronic sterile bladder inflammation of unknown etiology and it has not been possible to prove any causative pathogenic agent for this syndrome [2, 3]. Currently there are four major hypotheses of pathogenesis: 1) autoimmunity, 2) deficiency of the glycosaminoglycan layer causing increased bladder wall permeability, 3) neurogenic inflammation and 4) chronic infection .
While several features of IC have suggested an infective etiology, numerous studies using traditional culture techniques have failed to provide consistent evidence that IC is associated with infection. It has been proposed that possible microbial agents causing this disease could be difficult to cultivate or are present in numbers too low to be confirmed in the laboratory . Advances in molecular-based diagnostics have made it possible to overcome the limitations of culture-based detection. Investigators have used PCR, cloning and 16S ribosomal DNA (rDNA) sequencing to search for pathogenic agents in bladder biopsies and urine specimens of IC patients [6–11], but with conflicting results. However, some of these studies have indicated that women with IC may have a higher prevalence of bacteria in the urine than those without IC [6, 8, 9].
Furthermore, clinical studies have demonstrated that administration of antibiotics may sometimes be correlated with decreased symptoms in patients [12–14]. This can be due to both inhibition of bacterial growth or as a conventional anti-inflammatory effect of doxycycline. A study by Zhang et al. (2010)  not only demonstrates improvement in symptoms, but also a decreased level of nanobacteria after antibiotic treatment, strongly suggesting a microbial association of IC in some cases.
We recently developed approaches to assess the major microbial populations in female human urine, based on 16S rDNA PCR followed by 454 pyrosequencing and analyses using a suite of bioinformatics tools (Siddiqui et al. (2011) ) . We have shown that healthy female (HF) urine is a complex milieu with many different bacteria present. The normal human urine microbiota includes numerous fastidious and anaerobic microbes, which are potentially pathogenic [16–19]. In this work we applied these techniques in a prospective study to describe the microbial community present in urine from IC patients. We also performed a comparative analysis between the IC sequence dataset and the HF dataset previously generated  to determine to what extent the bacterial profiles differ. Our analyses indicate important differences between the two microbiota. We observe a lower complexity and variation between urine from IC individuals in relation to HF individuals.
This study was approved by the Regional Committee for Medical Research Ethics East –Norway (REK Øst Prosjekt 110-08141c 1.2008.367), and the samples were taken with informed consent. 8 female patients of age from 27 to 67 years (P1 = 59, P2 = 40, P3 = 27, P4 = 47, P5 = 31, P6 = 35, P7 = 32, P8 = 67) underwent thorough clinical examination including cystoscopy and fulfilled the criteria of European Society for the Study of Interstitial Cystitis (ESSIC) . All patients had an established diagnosis of IC for more than four years. Midstream urine (30 ml) was collected by the clean catch method with labial separation supervised by an urotherapy nurse. Specimens were kept at 4°C, and within an hour processed for DNA isolation. All specimens used were culture-negative, as tested by the Urological Clinic at the University Hospital HF Aker-Oslo. None of the patients was receiving antibiotics at the time samples were taken, nor prior to that according to hospital records.
Sample processing and DNA extraction was performed as previously described in Siddiqui et al. (2011) . Briefly, urine aliquots (30 ml) were pelleted by centrifugation and total DNA was isolated from sediments using DNeasy Blood & Tissue kit (QIAGEN, Germany), preceded by incubation with POWERlyse (lysis buffer) (NorDiag ASA, Oslo, Norway). Finally, the DNA was eluted in 100 μl of AE buffer from the kit. The DNA concentrations in the samples (P1-P8) were measured by Quant-iT PicoGreeen dsDNA assay kit (Molecular Probes, Invitrogen USA) and ranged from 0.22 ng/μl to 4.36 ng/μl.
For each IC urine sample, we amplified 16S rDNA sequences using two different primer sets specific for the V1V2 and V6 hypervariable regions followed by 454 pyrosequencing as described in Siddiqui et al. (2011) . Each of the primers consisted of a target specific region at their 3’ end (V1V2 or V6) and an adapter sequence (Primer A or Primer B) at their 5’ end as needed for GS FLX amplicon sequencing (454 Life Sciences, USA). Equal amounts of the two different amplicons (both V1V2- and V6-region) for a single subject were pooled and sequenced using GS-FLX chemistry in the same lane of a Pico-Titer plate divided into 16 lanes, except for samples P1, P2 and, P3, for which each amplicon (V1V2 and V6) was sequenced in a separate lane. 454 pyrosequencing was performed by the Norwegian Sequencing Centre (NSC) at the Department of Biology, University of Oslo, Norway.
Sampling depth and biodiversity found by amplicon 454 pyrosequencing V1V2 and V6 region from urine
Combined sequence pool from HF urine 1
Combined sequence pool from IC urine 2
Length cutoff 3
OTU and Diversity indices
Silva 16S alignment7
Shannon index10 (3%)
Inverse Simpson index11 (3%)
The bacterial identification technique of broad range 16S rDNA PCR is highly sensitive to environmental contamination. To control for this the IC urine sample sequence sets were stripped for sequences that could stem from contamination sources. This was done by using contamination control sequences (total = 25,246) from negative control extractions (buffer and kit reagents) followed by PCR and pyrosequencing, as reported in Siddiqui et al. (2011) . A complete linkage clustering at 1% genetic difference of each sample together with its respective control was performed using ESPRIT ( http://www.biotech.ufl.edu/people/sun/esprit.html). Any sample sequences found in clusters where ≥ 50% of the sequences belonged to the contamination control were excluded from subsequent analyses (details in Siddiqui et al. (2011) ).
Sequence data generated in this study were submitted to the Sequence Read Archive with the study accession number ERP001705. The dataset is available at http://www.ebi.ac.uk/ena/data/view/ERP001705.
For taxonomic grouping of the sequence reads, MEGAN V3.4 http://www-ab.informatik.uni-tuebingen.de/software/megan/welcome.html[23, 24] was used. First, the sequence reads were compared to a curated version of the SSUrdp database  using blastn with a maximum expectation value (E) of 10-5. To reflect the actual abundance behind every denoised sequence cluster, each entry in the blast result file was replicated as many times as the total number of reads that mapped to that query sequence (for detailed procedure and parameters see Siddiqui et al. (2011) ). When comparing the individual datasets using MEGAN, numbers of reads were normalized up to 100,000 for every dataset.
Metastats, statistical methods ( http://metastats.cbcb.umd.edu/, [26, 27]) for detecting differentially abundant taxa, was used to reveal significant differences between IC urine microbiota and HF urine microbiota (taxonomy assessed in Siddiqui et al. 2011 ). This method employs a false discovery rate to improve specificity in high-complexity environments, and in addition handles sparsely sampled features using Fisher’s exact test. The Metastats p - values at different taxon levels, which were assigned using MEGAN, are listed in Additional file 1: Table S1. A p - value ≤ 0.05 was considered significant.
Numbers of operational taxonomical units (OTUs), rarefaction curves and diversity indices were calculated using MOTHUR v1.22.2 [28, 29] (see Table 1). To enable comparisons, the HF sequences generated in Siddiqui et al. (2011)  were reanalyzed along with the IC dataset from this study. Briefly, the sequences were aligned to the Silva 16S alignment as recommended by MOTHUR  – sequences not aligned or aligned outside of where 95% of all of the sequences aligned were removed from the datasets. For an improved OTU clustering single linkage preclustering  was performed, allowing two nucleotides to differ between sequences, before clustering using average linkage. The processing was done both on each separate sample and on pooled V1V2 and V6 sequences for both IC and HF samples. We also calculated the OTUs and Shannon index for normalized numbers of sequences for each separate sample . A random number of reads, corresponding to the lowest number of sequences in a sample group, i.e. 2,720 for V1V2 and 2,988 for V6, was picked 100 times from each sequence set. These new sequence sets were processed through MOTHUR in the same fashion as the full sequence sets and the average of the resulting OTUs and Shannon values are shown in Additional file 2: Table S2.
The differences in Shannon indices and inverse Simpson indices for the two communities (IC and HF urine) were also statistically evaluated by Wilcox rank sum test in R, for both the V1V2 and V6 amplicon datasets.
Venn diagrams were generated for both data sets using MOTHUR to calculate how many OTUs were shared between the two communities. To further explore the relationships between the two microbial communities, samples were clustered into Newick-formatted trees (using the UPGMA algorithm) with distance between communities calculated with θYC coefficient as a measurement of dissimilarity between community structures  in MOTHUR. In addition, weighted UniFrac testing  was performed to determine the statistical significance of clustering within the tree. A non-metric multidimensional scaling (NMDS) plot was generated in R for the distances calculated using θYC measures for each sequence dataset (V1V2 and V6), knowing that θYC weighs rare and abundant OTUs more evenly than other metrics such as Jaccard.
After preprocessing of the raw IC 454 reads as described in Siddiqui et al. (2011) , we obtained a total of 46, 138 and 62,032 16S rDNA sequences for V1V2 and V6 regions, respectively, see Table 1. For comparison purposes, the preprocessing information for the HF urine sequences reported in Siddiqui et al. (2011)  is also listed in the table. Average number of reads per IC sample was 5,767 and 7,754 for V1V2 and V6, respectively (range: V1V2 3035–9506; V6 4900–14602) see Additional file 2: Table S2. 97% of the preprocessed sequences were classified to phylum, order and family level, and 95% of the sequences were identified down to genus level.
In comparison to HF urine (Siddiqui et al. (2011) ), IC urine has a significantly higher proportion of Firmicutes (p ≤ 0.05, p value from Metastats for V1V2) (65% vs 93%, respectively) and reduced proportions of the other 5 common phyla (Figure 1A). Interestingly, the phylum Nitrospirae was only detected in IC urine. Five additional phyla present in HF urine (Siddiqui et al. (2011) ) were not identified in IC urine at all (Figure 1A). The distribution of major phyla in IC urine was similar for both the V1V2 and V6 sequence dataset, although Fusobacteria and Nitrospirae were only identified by the V6 sequence dataset.
Sequence reads for all phyla but one (Nitrospirae 0.003% of the reads) were further classified to order level. 16 of the 22 orders identified in healthy urine (Siddiqui et al. (2011) ) were also detected in IC urine. A significant shift in the bacterial composition was observed as a result of an increase of Lactobacillales (Figure 1B and C) (p ≤ 0.05, p value from Metastats for V1V2) in the IC urine microbial community relative to HF urine. 92% and 91% of the reads for V1V2 and V6 respectively, were assigned to this order. In HF urine only 53% of the reads for V1V2 and 55% for V6 were assigned to Lactobacillales. The abundance of other major orders seen in HF urine is reduced in IC samples (Figure 1B and Additional file 1: Table S1).
All sequence reads assigned to the order level were additionally assigned to family level. Among the 26 families identified, only 21 were assigned to different genera. Four of those families that were not further assigned (Pasteurelacae, Neisseriacae, Methyliphilaceae, and Micrococcaceae) were also detected in the HF urine study. Saprospiraceae, on the other hand was only found in IC urine.
At the genus level, the pooled sequences were assigned to 31 different genera, with 23 and 25 different genera for V1V2 and V6 analysis, respectively. Lactobacillus was the most abundant genus in both datasets and comprised a total of 92% of the sequences. Gardnerella and Corynebacterium were the two other major genera identified with 2% sequence abundance each. Prevotella and Ureaplasma were each represented by 1% of the sequences assigned. The other 26 genera determined in IC urine constituted only 2% of the total IC urine bacterial community.
In contrast to HF urine, there was a considerable reduction in total numbers of genera identified in IC urine (45 genera vs. 31 genera, respectively) (Additional file 1: Table S1). Additionally, the abundance of common genera was found to differ between IC patients and healthy females. The significant increase of Lactobacillus (p ≤ 0.05, p values from Metastats for both V1V2 and V6) in IC urine compared to HF urine again suggested a structural shift in the microbiota of IC patients. Enterococcus, Atopobium, Proteus and Cronobacter are 4 genera identified in IC urine that were not detected in our previous HF urine study, while a group of 17 genera were only associated with HF urine.
Sequences from 16 of the genera identified in the IC samples were further assigned to 22 different species (Additional file 3: Table S3). When comparing to our previous study, 13 of these species are already found in asymptomatic HF urine. However, nine of these species were not identified in our previous study, nor associated with IC according to literature.
In all but two IC urine samples, Lactobacillus accounted for more than ~95% of the sequences for both V1V2 and V6 data. Lactobacillus was not only the most abundant genus, but also the most frequent genus among all IC urine specimens with its rRNA sequences present in all eight samples, in contrast to urine samples from HF (6/8). Sequences assigned to Prevotella, Peptoniphilus and Anaerococcus were also frequently detected (5/8), followed by Staphylococcus and Finegoldia (4/8), and Gardnerella, Streptococcus and Dialister (3/8) in IC urine. Including Ureaplasma, 7 genera were identified by reads belonging to 2 urine samples and another 15 genera were only detected in 1 out of the 8 samples.
We have characterized the urine microbiota of IC patients using high throughput 454 pyrosesequencing of 16S rDNA amplicons. These results were compared to HF data from our previous study (Siddiqui et al. (2011) ). Our results did not reveal any single potential pathogenic bacterium common to all IC patients. However, important differences were detected between the IC and HF microbiota.
The use of primers for both V1V2 and V6 regions yielded complementary results for IC urine in line with the previous study of HF urine (Siddiqui et al. (2011) ), and thus maximized the detection of bacterial diversity.
Knowing that urine samples are at risk of contamination by bacterial flora of the female urogential system [34, 35], mid-stream urine sampling was performed under guidance of an experienced urotherapy nurse. Suprapubic puncture was suggested as an alternative method, but the method was considered to be too invasive. Interestingly, comparing results from our previous microbiome study on female mid-stream urine (Siddiqui et al. 2011)  with recent results from suprapubic aspirate by Wolfe et al. (2012) , the major findings are the same; a strong indication that mid-stream urine will give comparable results in a urine microbiome analysis.
A decrease in overall richness and ecological diversity (as indicated by rarefaction analysis, number of OTUs, Shannon index and inverse Simpson index estimations) of IC urine microbiota was detected in contrast to HF urine (Table 1 and Figure 3). In addition, the ß-diversity analysis (θYC distances between all urine samples) suggested that the microbiota of HF samples are more dissimilar from each other than the microbiota of IC individuals.
The taxonomical analysis indicated a shift in composition of urine microbiota of IC patients, with changes in bacterial groups spanning from genus to phyla level and a reduction in microbial complexity compared to HF. More importantly, a significant increase in Lactobacillus in IC patients was revealed. Lactobacillus was detected in every IC urine sample and a total sequence abundance of 92% in contrast to 57% in HF urine was observed. This shift was also clearly displayed both at the order and phylum level (Lactobacillales and Firmicutes, respectively). In contrast, Prevotella, - a genus belonging to the phylum Bacteroidetes (order Bacteriodales) – was present only at 1%, significantly lower than in HF urine, where it was previously reported as one of the major genera with an abundance of 19%. Gardnerella, another dominant genus in female urine, was present with the same frequency in IC urine but with a general lower abundance.
A reduction in bacterial diversity and shift in the microbiota as observed in this chronic inflammatory state has also been reported for other clinical conditions such as obesity, irritable bowel syndrome, and inflammatory bowel disease including Crohn’s disease [36–38].
Attempts to identify an infectious etiology for IC have not yet found any evidence for a specific pathogen. However, previous culture-dependent studies of samples from IC patients (i.e. bladder biopsy, midstream urine) have reported organisms such as Gardnerella, Lactobacillus sp., Streptococcus ssp., Escherichia coli, Proteus mirabilis, Corynebacterium ssp., Klebsiella sp., Enterococcus sp., Propionbacterium, Prevotella, Bacteroides sp., and Peptostreptococcus[6, 9, 39]. Lactobacillus, Gardnerella and Streptococcus were repeatedly detected in these studies and were also seen in our study. Haarala et al. (1999)  using culture techniques concluded that bacterial flora of midstream urine from patients with IC clearly differs from that of healthy women, in line with our findings. A study by Zhang et al. (2010)  suggested nanobacteria as a possible causative agent for IC. The two latter studies also reported a reduction in bacterial levels and urinary symptoms upon antibiotic treatment of the IC patients. The primer pairs both for V1V2 and V6 amplicons used in our study would be expected to amplify 16S rDNA regions of all of the organisms mentioned above. Nevertheless we did not identify Klebsiella, E.coli, Peptostreptococcus or nanobacteria in any of our IC urine samples.
Studies reporting results from culture-independent 16S rDNA PCR approaches on samples (i.e. bladder biopsy, midstream urine) from IC patients, have yielded somewhat conflicting results both in terms of positive PCRs and the resulting bacterial profiles [7, 8, 10, 11, 40]. While two of the reports [11, 40] found no evidence of bacterial DNA in biopsy and urine specimens from IC patients, Dominique et al. (1995)  demonstrated bacterial DNA in bladder tissues in 29% of patients with IC. The 4 sequences retrieved showed homology to E. coli (2) and Pseudomonas (2), however neither of these bacteria was found in our study. Heritz et al. (1997)  also reported bacterial DNA in both biopsies and urines from IC patients (53% and 46%, respectively). They concluded that there is a difference in the bacterial profile between the patients and controls, and further suggested a link between one or more bacterial species and IC. However, in their study only 11 bacterial clones from 3 different IC patients were analyzed and the bacterial sequences were related to E.coli, Abiotrophia defectivus, Veillonella and Rothia dentocariosa. Except for Veillonella, these bacteria were not detected in our study. All these 4 previously reported studies used different primer sets (likely to explain some of the differences in the results) and classical cloning strategies (explaining the very few sequences analyzed). In contrast, our study represents the first 16S rDNA amplicon high throughput sequencing approach on IC urine, increasing both the sensitivity and resolution of the investigation.
Lactobacillus has not only been indicated or shown in IC urine samples from females (100% of the cases in this study and as shown by others [6, 9, 39]) but also demonstrated in IC urine from a male subject . In our study we also detected a significant increase in abundance of this genus, considering its supposedly commensal presence in human urine from healthy subjects [16, 18, 19].
Lactobacillus is generally considered to be of low virulence, rarely causing infections in humans. It is best known for its presence in vaginal microflora, where it normally generates and maintains a physiological acidic environment, which prevents infections. Because of these properties, Lactobacillus has been used in probiotics, and is thought to prevent or even treat urinary tract infection (UTI) . However, there are increasing indications that specific Lactobacillus spp are of pathogenic relevance and may be involved in urinary tract infections [43, 44].
Many female patients with symptoms suggestive of UTI, but with culture-negative urines are often treated with antibacterial agents since their symptoms may be indistinguishable from those with a proven UTI . It has been proposed that Lactobacillus, resistant to widely used antibiotics, may multiply during treatment, giving the genus an advantage over antibiotic-sensitive commensals, and allowing it to invade the proximal urethra and paraurethral tissues causing inflammatory changes . This organism has also been related to the presence of UTI symptoms in otherwise culture-negative urines [43, 44, 46]. In a study by Maskell et al. (1983)  antibacterial treatment was withheld over the course of 2 years from symptomatic women with culture-negative urine. During the course of the study Lactobacilli (detected by special culture techniques) gradually disappeared from the urine of most of the patients who also became symptom free. A similar association of Lactobacillus and urinary symptoms was reported by Darbro et al. (2009) . Over a six-month period a female patient had constant symptoms and high counts of Lactobacillus spp (> 50,000 CFU/ml) in the urine; she became symptom-free with culture-negative urine after treatment with Lactobacillus-targeted medication. These results also suggest that a shift in the microbial community towards Lactobacillus in IC urine samples may be an important etiological factor for the severe symptoms reported by the patients. Since additional culture techniques such as 48 h incubation in an atmosphere containing 7% CO2 are needed for detection of Lactobacillus, this may be the reason why IC urine samples have not yet been associated with bacterial growth in routine clinical investigations. However, in our study this problem was overcome by a culture-independent approach.
This investigation did not reveal any obvious putative causative bacterial agents of IC. However, the greater abundance of Lactobacillus in IC urine and its lower occurrence in HF urine is an important finding that requires further study to establish whether these microbial changes play a part in the development of IC. To this end, whole genome sequencing of Lactobacillus from IC patients may be a possible approach. Even if an increased presence of Lactobacillus is merely a secondary marker, understanding its IC associated genomics could aid in diagnosis and therapeutic assessment.
The authors would like to thank Hege Junita Gaup for technical assistance, William Ryan Easterday for critical reading of the manuscript and the Norwegian Sequencing Centre (NSC, www.sequencing.uio.no), Department of Biology, University of Oslo, for sequencing services. We are very grateful to Professor Lars M Eri and urotherapists Turid H Hoel and Bodil Svendsen at Aker University Hospital HF, Urological Clinic for specimen collection. Additionally we thank two anonymous reviewers, whose comments helped to improve the manuscript.
Financial support for this research was provided by grants from the Research Council of Norway to KSJ and from CEES to HS.
This article is published under license to 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.