Open Access Highly Accessed Research

Integrated next-generation sequencing of 16S rDNA and metaproteomics differentiate the healthy urine microbiome from asymptomatic bacteriuria in neuropathic bladder associated with spinal cord injury

Derrick E Fouts1*, Rembert Pieper1, Sebastian Szpakowski1, Hans Pohl2, Susan Knoblach2, Moo-Jin Suh1, Shih-Ting Huang1, Inger Ljungberg3, Bruce M Sprague2, Sarah K Lucas1, Manolito Torralba1, Karen E Nelson1 and Suzanne L Groah34

Author Affiliations

1 J Craig Venter Institute, 9704 Medical Center Drive, Rockville, MD, 20850, USA

2 Childrens National Medical Center, 111 Michigan Avenue NW, Washington, DC, 20010, USA

3 MedStar National Rehabilitation Hospital, 102 Irving Street NW, Washington, DC, 20010, USA

4 Georgetown University Hospital, 3800 Reservoir Road, Washington, DC, 20007, USA

For all author emails, please log on.

Journal of Translational Medicine 2012, 10:174  doi:10.1186/1479-5876-10-174


The electronic version of this article is the complete one and can be found online at: http://www.translational-medicine.com/content/10/1/174


Received:3 May 2012
Accepted:2 August 2012
Published:28 August 2012

© 2012 Fouts 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

Clinical dogma is that healthy urine is sterile and the presence of bacteria with an inflammatory response is indicative of urinary tract infection (UTI). Asymptomatic bacteriuria (ABU) represents the state in which bacteria are present but the inflammatory response is negligible. Differentiating ABU from UTI is diagnostically challenging, but critical because overtreatment of ABU can perpetuate antimicrobial resistance while undertreatment of UTI can result in increased morbidity and mortality. In this study, we describe key characteristics of the healthy and ABU urine microbiomes utilizing 16S rRNA gene (16S rDNA) sequencing and metaproteomics, with the future goal of utilizing this information to personalize the treatment of UTI based on key individual characteristics.

Methods

A cross-sectional study of 26 healthy controls and 27 healthy subjects at risk for ABU due to spinal cord injury-related neuropathic bladder (NB) was conducted. Of the 27 subjects with NB, 8 voided normally, 8 utilized intermittent catheterization, and 11 utilized indwelling Foley urethral catheterization for bladder drainage. Urine was obtained by clean catch in voiders, or directly from the catheter in subjects utilizing catheters. Urinalysis, urine culture and 16S rDNA sequencing were performed on all samples, with metaproteomic analysis performed on a subsample.

Results

A total of 589454 quality-filtered 16S rDNA sequence reads were processed through a NextGen 16S rDNA analysis pipeline. Urine microbiomes differ by normal bladder function vs. NB, gender, type of bladder catheter utilized, and duration of NB. The top ten bacterial taxa showing the most relative abundance and change among samples were Lactobacillales, Enterobacteriales, Actinomycetales, Bacillales, Clostridiales, Bacteroidales, Burkholderiales, Pseudomonadales, Bifidobacteriales and Coriobacteriales. Metaproteomics confirmed the 16S rDNA results, and functional human protein-pathogen interactions were noted in subjects where host defenses were initiated.

Conclusions

Counter to clinical belief, healthy urine is not sterile. The healthy urine microbiome is characterized by a preponderance of Lactobacillales in women and Corynebacterium in men. The presence and duration of NB and method of urinary catheterization alter the healthy urine microbiome. An integrated approach of 16S rDNA sequencing with metaproteomics improves our understanding of healthy urine and facilitates a more personalized approach to prevention and treatment of infection.

Keywords:
Bacteriuria; Urine; Catheter; Neuropathic; Bladder; Microbiome; Metaproteome; Next-generation; Personalized; rRNA

Background

Affecting nearly one half of all Americans over the course of a lifetime [1] and with costs exceeding $1 billion annually [2], urinary tract infection (UTI) is a major public health problem [3,4]. It is the most common urologic disorder in the outpatient setting [3,4] and the most common health care associated infection [3,5,6]. In the health care setting alone, approximately 561777 UTIs occur annually, costing an estimated $1006 per infection, totaling more than $500 million, and being responsible for 8205 deaths [5]. This does not include the personal suffering or time lost from gainful employment.

Often considered an antecedent to UTI, asymptomatic bacteriuria (ABU) represents an asymptomatic carrier state recognized to have little impact on health or quality of life. This is in contrast to healthy urine, considered to be sterile until reaching the urethra, which is colonized by facultative anaerobic Gram-negative rods and cocci. The most common risk factor for the development of ABU and UTI is use of a urinary catheter [5-7], as catheters provide a conduit for bacterial colonization and symptomatic infection. Resolution of ABU typically occurs with removal of the urinary catheter. This is not possible in many cases, however, as the urinary catheter facilitates function and emptying in cases of bladder impairment.

The clinical distinction between symptomatic UTI and ABU is not trivial since symptomatic UTI requires treatment and perhaps further evaluation irrespective of the circumstances in which the UTI occurred, while ABU does not. Exceptions to this rule include the treatment of ABU in select at-risk populations, such as pregnant women, which has been shown to be associated with improved outcomes [8]. Distinction between these states is particularly relevant, as renewed urgency and heightened focus have been placed on UTI occurrence by national policy-makers and payers. The Centers for Medicare & Medicaid Services (CMS) has identified catheter-associated UTI, the most common hospital acquired infection, as a “never event”. Effective in 2008, this has resulted in non-reimbursement for catheter-associated UTIs that were not present on admission to acute care hospitals [9,10]. An unintended consequence of designating catheter-associated UTI as a “never event” is more aggressive screening for ABU and UTI upon admission of patients to hospitals, a strategy that may lead to increased unnecessary antibiotic treatment and emergence of antimicrobial resistance.

To achieve improved outcomes in the care of patients with ABU and UTI, improved information distinguishing states of urine in health and disease is needed. To this end, we sought to first describe states of urinary health utilizing a highly sensitive, culture independent approach to determine whether the urine microbiome of healthy people who are at risk for ABU because they utilize urinary catheters differs from that of healthy controls, and if so, to identify key factors or bacterial signatures that might ultimately lead to UTI requiring antimicrobial treatment. Subjects with neuropathic bladder due to spinal cord injury who are known to be at highest risk for ABU and UTI due to their need for catheter-assisted bladder management [4,11,12] were assessed and compared with healthy controls to achieve our goal.

Methods

The study was approved by the MedStar Institutional Review Board (IRB). All study personnel were certified in and the study protocol conformed to the ethical guidelines of the 1975 Declaration of Helsinki as reflected in approval by the MedStar IRB.

Sample acquisition and clinical urinalysis

Patients and healthy controls were recruited into this IRB-approved study (NRH IRB# 2011–019) from the outpatient clinic and inpatient ward at National Rehabilitation Hospital (Washington, DC). Following written consent, urine samples were obtained from 26 healthy, non-SCI controls and 27 with neuropathic bladder (NB) due to spinal cord injury (SCI). Patients provided urine samples by sterile collection using the means by which they customarily empty their bladder (i.e. midstream collection during voiding, or sterile catheterization if unable to void). (Table 1) The samples were coded with an anonymous research identification number and separated into two aliquots: one, for standard urine analysis and culture (Quest Diagnostics) and another for microbiome analysis. Quest Diagnostics performed analysis of urine samples for nitrite formation, leukocyte esterase, and microscopic examination for the presence and quantity of leukocytes and erythrocytes in each sample. Bacterial cultures were performed by inoculation of blood agar plates and incubation at 37°C for 48 hours.

Table 1. Patient demographics

Sample preparation and PCR

Thawed urine samples were clarified by low-speed centrifugation and bacterial genomic DNA was extracted from urine pellets by enzymatic digestion using a final concentration of 20 mg/ml Lysozyme (Invitrogen) followed by physical lysis using Lysing Matrix B tubes (QBiogene). A previous study comparing mechanical and enzymatic methods for extracting microbial genomic DNA showed that mechanical cell disruption by bead beating produced the highest bacterial diversity [13]. The samples were vortexed at maximum speed for 45 seconds using a Fastprep fp120 (MP Biomedicals) then cooled on ice. DNA was extracted from the lysate using phenol chloroform isoamyl alcohol extraction and ethanol precipitation. 16S rDNA sequences were generated by amplifying the V1-V3 region of the bacterial 16S rRNA gene using primers 27F and 534R fused with 454 adaptors and barcodes for multiplexing. Primers targeting V2 and V3 were shown to perform as well as full-length 16S rDNA sequence for community clustering and taxonomic assignments [14].

The amplicons were normalized and pooled prior to emulsion PCR and 454 sequencing (Roche, Inc.) using titanium chemistry.

DNA sequence processing

The 16S rDNA sequence-processing pipeline used for this study is composed of a selection of bioinformatics tools proven to be accurate, robust and fast. A supplementary archive contains the dot language representation of a graph depicting the entire workflow executed, outlining the specific parameters used for each command.

Initially, the SFF file, output from the sequencer, was converted into fasta and qual files using the sffinfo program included as a part of 454/Roche software package. Subsequently, the trim.seqs function in mothur [15] (version v.1.22.2) was used to de-multiplex sequencer reads. No barcode mismatches, and up to one primer mismatch were allowed past this step. The de-multiplexed reads were processed using LUCY [16-18] to filter out reads with low quality segments. At this point, the sub.sample function of mothur was used to select an equal number of reads per biological sample (n=3671 based on the biological sample with the fewest number of reads). Subsequent, the screen.seqs function of mothur was used to remove sequences shorter than 220 bases [17]. Furthermore CD-HIT-454 [19,20] was used to collapse duplicate reads, while retaining their count for subsequent enrichment statistics, (analogous to the functionality of the unique.seqs function of mothur, but orders of magnitude faster and less demanding on the computer hardware). The sequences were aligned against the SILVA database of 16S rDNA sequences [15,21] to verify 1) the orientation of noise-filtered sequences; and 2) the correct positioning of the reads with respect to the expectation of which variable regions should have been amplified and sequenced. Thereafter, the remaining sequences were subjected to mothur’s implementation of chimera slayer [15,22] to filter out chimeric reads. The processed 16S rDNA data from this study can be obtained at NCBI under BioProject ID 97505.

Taxonomical classification of OTU representative reads

Taxonomical classification of the final set of 82160 operational taxonomic unit (OTU)-representative reads down to the genus level was performed using mothur’s version of the RDP Bayesian classifier using a normalized RDP training dataset [23]. The final step of the pipeline clustered the sequences based on their similarity to produce OTUs. Customarily, a similarity threshold of 97% has been used to define OTUs at approximately the species level [24]. A module of CD-HIT suite [19] called CD-HIT-EST was employed to perform species-level read-clustering for subsequent analyses.

Statistical analyses

The orchestration and automation of steps has been achieved using a custom set of in-house utilities written in python and R [25] programming languages. These utilities are available online as a part of the YAP package on github [26]. JCVI grid infrastructure based on the Oracle Grid Engine (OGE) was used for all steps described. Relative abundance and diversity statistics were calculated within mothur [15]. Further statistical analyses were accomplished using R. Heat maps were generated using the heatmap.2 function of the gplots package available on CRAN [25]. OTU counts have been normalized to 100% per individual, to facilitate comparability. Only taxa (order or genus) with a standard deviation greater than 5% across all 52 individuals were used to generate the heat maps. Differences between subject OTU communities were assessed using the Bray-Curtis beta diversity statistic implemented in the vegdist function, a part of vegan R package available on CRAN. Clustering was accomplished using average-neighbor-joining method implemented in hclust function in the default installation of R. PCA analysis was performed using ade4[27]. P-values used to determine statistical significance in relative OTU difference plots were established using the default installation of R and kruskal.test functions implementing the Kruskal-Wallis rank-sum statistic test [28].

Phylogenetic tree building

OTU-representative sequences classified as either Lactobacillales or Enterobacteriales were aligned using tools available from Release 10 of the RDP web site [29]. Specifically, sequences from OTUs composed of reads from more than one individual were aligned to the RDP reference 16S rRNA sequence, taking into account secondary structure. At most one nearest neighbor sequence from RDP was recruited into the alignment per input sequence. The alignment was downloaded and trimmed to remove columns whose gap fractions were greater than 50%, using Belvu [30]. Based on the alignment, a bootstrapped Neighbor-Joining (NJ) tree was subsequently inferred using paupFasta, an in-house wrapper script around the PAUP* program as described [31], and edited using Fig Tree [32]. In combination with nearest neighbor taxonomy, BLASTN was used against the NCBI reference RNA database, which lacks uncultured organisms, to identify certain Lactobacillales OTU-representative branches that lacked classified RDP top matches, to the species level, using a cut-off of 97% identity.

Proteomics

A urinary pellet specimen equivalent to 5–10 ml voided urine washed two times with ~10 ml ice-cold PBS was re-suspended in 1 ml of 10 mM ammonium bicarbonate containing 0.1% Triton-X100, 0.5% octylglucoside, 5 μg/ml leupeptin, 10 mM EDTA and 2 mM benzamidine. The suspension was heated to 85°C for 5 min followed by sonication (amplitude 4, Misonex 3000 sonicator) in 30s on/15s off cycles 10 times on ice. The suspension was centrifuged for 15 min at 16,100 x g and the supernatant recovered. Following protein quantity estimates based on SDS-PAGE analysis, an aliquot with ~10 μg protein was digested at a 1:50 ratio (trypsin/protein) using Filter-Aided Sample Preparation (FASP) [33]. The digestion mixture was reconstituted in 50 μl 0.1% formic acid. For shotgun proteomic analysis, peptides in a 20 μl sample aliquot were separated on a capillary C18 LC column in 122 min binary gradient runs from 97% solvent A (0.1% formic acid) to 80% solvent B (0.1% formic acid, 90% acetonitrile) at a flow rate of 350 nl/min. Nano-electrospray into the source was followed by mass analysis and spectral acquisition in automated MS/MS mode, with the top five parent ions selected for fragmentation in scans of the m/z range 350–2,000 and a dynamic exclusion setting of 90 sec. The LC-MS/MS workflow using the LTQ-XL ion trap system (Thermo Fisher Scientific) was previously described in more detail [34]. The instrument was calibrated prior to performance of LC-MS/MS experiments with 200 nmol human [Glu1-fibrinopeptide B (M.W. 1570.57), verifying that peaks representing ion counts had widths at half-height of <0.25 min, signal/noise ratios >200 and peak heights >107. The LTQ search parameters (+1 to +3 ions) included mass error tolerances of ± 1.4 Da for peptide precursor ions and ± 0.5 Da for peptide fragment ions, selecting monoisotopic m/z values. The search engine used to select these parameters and identify peptides and proteins was Mascot v.2.3 (Matrix Science). The protein sequence database was comprised of 19 bacterial genomes (details in Additional file 1) and Uniref90’s human database subset [35]. We limited this database to the human Uniref90 subset and 19 bacterial genomes frequently associated with ABU and UTI, because significant computational challenges for peptide-spectral match (PSM) assignments are encountered when very large protein databases are used [36]. As described in Additional file 1, we used stringent criteria for PSMs (q-value <=0.01; PEP-value <=10-4) using the Mascot Percolator algorithm that improves discrimination between correct and incorrect PSMs, particularly when the searched sequence space in the database is large [37].

Additional file 1. Proteomic analyses of urinary pellet samples using the LTQ XL ion trap instrument (Thermo-Electron) and the Mascot search engine version 2.3 (Matrix Science) for spectral matches with a 19-species database.

Format: PDF Size: 1.6MB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Results

Urine samples were obtained from 26 healthy, non-SCI controls and 27 individuals with NB due to SCI at the National Rehabilitation Hospital in Washington, DC. Among those with SCI and NB, 8 voided spontaneously without a catheter, 8 emptied by clean intermittent catheterization (IC), and 11 used indwelling urethral Foley catheterization (FC) (Table 1). Frozen, resuspended pellets from the urine samples were analyzed with culture-independent surveys of bacterial 16S rDNA and urinary proteins.

A total of 589454 quality-filtered 454 16S rDNA sequence reads were processed through a NextGen 16S rDNA analysis pipeline [26] (Additional file 2) where taxonomy of species-level operational taxonomic units (OTUs) was determined using the RDP classifier [38]. Analysis of the bacterial community in these urine samples revealed between 5 and 236 species-level (97% identity) OTUs per individual (Additional file 3). Per sample taxonomic profiles were generated, showing considerable sample-to-sample variation (Additional file 4). To better visualize sample-to-sample taxonomic profiles, a heat map was generated, clustering the distribution of OTU taxonomy at the level of bacterial order using the Bray-Curtis index (Figure 1A). The top ten bacterial taxa showing the most relative abundance and change among samples were Lactobacillales, Enterobacteriales, Actinomycetales, Bacillales, Clostridiales, Bacteroidales, Burkholderiales, Pseudomonadales, Bifidobacteriales and Coriobacteriales (Figure 1A). The Lactobacillales and Enterobacteriales were the two most relative abundant and changing taxonomic groups.

Additional file 2. YAP 16S rDNA sequence-processing and analysis pipeline workflow diagram.

Format: PDF Size: 53KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 3. A table summarizing the results of bacterial OTU-based microbiome analysis.

Format: PDF Size: 174KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 4. Bacterial OTU taxonomic profiles at the genus level per individual.

Format: PDF Size: 3.9MB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 1. Differences in relative bacterial OTU counts between neuropathic and healthy bladder in males and females. For every individual, the OTU counts were normalized to the individual's total OTU count. A heat map of the clustered distribution of OTU taxonomy at the level of bacterial order (A) and genus (B) was constructed using the Bray-Curtis index. For panels A and B, only the taxa with a standard deviation > 5% across all individuals are shown. Differences in the average OTU count between females and males are plotted in light and dark gray, respectively (C). In panel (C), the top 15 (8%) most abundant bacterial OTUs are represented. The X-axis indicates the difference in relative OTU counts per bacterial genus indicated on the Y-axis. Statistical significance was established using Kruskal-Wallis test. Significant differences (P < 0.05) between the relative OTU counts are indicated by an asterisk (*) for bladder function, and a plus sign (+) for gender.

Urinary microbiome differed by gender and bladder function

To determine whether distinct microbial signatures were associated with gender (male versus female) and/or bladder function (healthy control versus NB), samples were grouped by these variables. Taxonomic counts were then normalized by total number of OTUs per sample, and visualized by a heatmap, clustering the distribution of OTUs at the level of bacterial genus using the Bray-Curtis index (Figure 1B). The samples were colored by relative abundance (red/warm most abundant to blue/cool least abundant). Six main taxonomic profile clusters emerged with distinct patterns when grouped by gender and bladder function. Cluster 1 was composed of almost an equal proportion of healthy and NB samples, clusters 2 and 3 were dominated by healthy controls, whereas clusters 4–6 were entirely composed of patients with NB (Figure 1B). The 2 “healthy” clusters, (2 and 3) were distinguished by gender, with females in cluster 2 and males in cluster 3, and by bacterial genus, Lactobacillus grouping with females in cluster 2 and completely absent in the male-dominated cluster 3. Cluster 3 was composed of different Gram-positive organisms. Cluster 1 had the most diverse bacterial genus profile, composed largely of Lactobacillus, but not as abundant as cluster 2 with elements of cluster 3 and a few potentially pathogenic genera (Enterococcus, Salmonella, and Peptoniphilus). Closer inspection of clusters 4–6 showed a very different pattern of bacteria, with known UTI pathogens dominating the profiles (Figure 1B). Cluster 4 was primarily composed of Enterococcus, Escherichia and Salmonella. Cluster 5 was dominated by Klebsiella sp. and was the only cluster comprised of all males. Cluster 6 had the most Enterococcus counts of any cluster and also contained Aerococcus and Proteus sp.

By taking the difference in normalized relative abundance between controls and SCI groups, the most relative abundant bacterial taxa per group as confirmed (Figure 1C). Statistical significance was established using the Kruskal-Wallis test. Significant differences (P < 0.05) between the top relative OTU counts suggest that Lactobacillus, Corynebacterium, Gardnerella, Prevotella and Enterococcus define gender differences (“+” in Figure 1C). Lactobacillus, Klebsiella, Corynebacterium, Staphylococcus, Streptococcus, Gardnerella, Prevotella, Escherichia and Enterococcus defined statistically significant differences between healthy bladder controls and NB (“*” in Figure 1C).

Urinary microbiome differed with duration of NB

To further investigate the dynamics of the urinary microbiome in NB, the NB group was divided into bins based on the time post SCI. The bins consisted of 0–2 months (n=5), 3–12 months (n=5), 13–48 months (n=3), and 49+ months (n=14). Principal component analysis (PCA) revealed that the healthy control group and the 0–2 month duration were very similar (Figure 2). Likewise, the 13–48 and 49+ month groups were similar to each other, while the 3–12 month group was too distorted by one outlier to provide meaningful information. Enterococcus and Escherichia emerged as major contributors to these profiles as illustrated in the vector diagram (inset, Figure 2).

thumbnailFigure 2. OTU differences among individuals by duration of neuropathic bladder. A PCA analysis of the OTU counts of 52 individuals. The points are circled and colored based on the duration (in months) of neuropathic bladder (see key). The inset depicts a vector plot indicating the most influential principal component (bacterial genus).

Urinary microbiome differed with bladder management

To get an overview of bacterial species variation within the NB group, differences in microbiome taxonomic profiles between males and females and among bladder management groups were plotted (Figure 3). Lactobacillus species were significant (P < 0.05) contributors to female healthy controls and the void NB groups, while Corynebacterium sp. defined the healthy bladder male urinary microbiome (Figure 3). The preponderance of Enterococcus sp. in male NB observed in Figure 1C was primarily within the patients with NB who emptied by spontaneous voiding (Figure 3).

thumbnailFigure 3. Differences in the relative OTU counts between males and females stratified by bladder management. For every individual, the OTU counts were normalized to the individual’s total OTU count. Differences between relative OTU counts were calculated by subtracting the average OTU counts from females and males per bladder management category; healthy control, void (SCI patient with no catheter usage), indwelling catheter (IC), and Foley catheter (FC) (see key for color coding and sample sizes). The X-axis indicates the difference in relative OTU counts per bacterial genus indicated on the Y-axis. Significant differences (P < 0.05) between the relative OTU counts are indicated by plus sign (+) for gender, and an asterisk (*) for bladder management. The inset depicts the mean and standard deviation of OTU counts of the indicated genera for each catheter management group.

Phylogenetic tree of the urinary Lactobacillales revealed potentially pathogenic species

A phylogenetic tree was constructed with OTU representatives of all OTUs with matches to Lactobacillales and nearest neighbors in the RDP database in order to provide an environmental context and determine species-level taxonomy (Figure 4). Labels of the tree were colored based on NB composition. Red indicates OTUs composed only of samples from NB, light red indicates OTUs with a majority composition of NB samples, while dark blue denotes OTUs only solely of healthy controls, and light blue illustrates OTUs with a majority of reads from healthy controls. There is much greater phylogenetic diversity among the Lactobacillus and Streptococcus branches (green and gray shading, respectively, Figure 4). Aerococcus and Enterococcus branches are dominated by OTUs (red leaves) from urine of NB patients, suggesting qualitatively that these two genera may be important indicators of bacteriuria, or of catheter usage. In contrast to these two branches, the Lactobacillus and Streptococcus branches are largely comprised of OTUs from healthy controls with a few exceptions (red leaves). One such exception was Lactobacillus iners (Figure 4), which was previously shown to significantly contribute to UTI [39,40]. Another exception was determined to be Streptococcus salivarius, a lactic acid-producing Gram-positive organism typically found in the oral cavity, and an opportunistic pathogen implicated in bacteremia [41-43] and septicemia [44].

thumbnailFigure 4. Phylogenetic diversity of Lactobacillales 16S rDNA sequences in human urine. NJ tree clustering of Lactobacillales OTU representatives labeled based on similarity to known RDP database sequences (gray), and OTU composition. Leaves are colored as follows: OTUs consisting of only healthy individuals (dark blue), mostly healthy (light blue), only NB (red), mostly NB (pink/salmon). Branches were highlighted and labeled by identifiable bacterial genera. Genus-level classification was based on the OTU representative RDP classification and the classification of nearest neighbors the RDP alignment. The nodes show SequenceID_#male/#female_#SCI/#healthy subjects.

Phylogenetic diversity of the urinary enterobacteriales

In contrast to the Lactobacillales, the Enterobacteriales had no OTUs composed solely of healthy controls; indicating this group of bacteria may be a potential indicator of future UTI (Figure 5). Shaded regions of the tree noted unambiguous genus taxonomy as follows: Escherichia, Enterobacter, Klebsiella, Proteus, Morganella, and Providencia. All of these genera have been associated with UTI, and dominated the tree. Klebsiella stood out from the others, having been divided into three distinct branches. Top matches to N2-fixing, plant-associated Klebsiella pneumoniae and K. variicola isolates rather than known K. pneumoniae clinical strains suggests that N2-fixing, plant-associated Klebsiella may exist in urine. Further work is needed to confirm this result.

thumbnailFigure 5. Phylogenetic diversity of Enterobacteriales 16S rDNA sequences in human urine. NJ tree clustering of Enterobacteriales OTU representatives labeled based on similarity to known RDP database sequences (gray), and OTU composition. Leaves are colored as follows: OTUs consisting of only healthy individuals (dark blue), mostly healthy (light blue), only NB (red), mostly NB (pink/salmon). Branches were highlighted and labeled by identifiable bacterial genera. Genus-level classification was based on the OTU representative RDP classification and the classification of nearest neighbors the RDP alignment. The nodes show SequenceID_#male/#female_#SCI/#healthy subjects.

Metaproteomics reveals anti-microbial and pro-inflammatory responses in the absence of diagnosed UTI

Shotgun proteomic data for a subset of nine urinary samples searching a protein sequence databases of 19 bacterial genomes were compared with the OTU analysis derived from 16S rDNA sequencing. In addition to protein identifications targeting the 19 species frequently associated with UTIs or urethra colonization, a non-redundant human protein database was searched to assess the detection of host responses towards bacterial colonization. While proteomic analysis identified more species than those diagnosed by culture methods, it did not identify the fastidious anaerobic or microaerophilic bacteria frequently assigned by the RDP classifier in the 16S rDNA analysis with the exception of Lactobacillus. We assume this to be due to species abundance issues. Proteins encoded by Escherichia, Klebsiella and Enterobacter species were identified in half of the examined samples, generally in agreement with 16S rDNA data (Table 2). Metaproteomic profiles allowed preliminary insights into the production of bacterial factors interacting with the urinary tract environment. Identifications of such proteins as limited to those urine donor samples with evidence for the initiation of pro-inflammatory and microbicidal host defenses (high spectral counts for calprotectin subunits; identification of lactotransferrin, myeloperoxidase and eosinophil peroxidase). We emphasize that the sample sizes are too small to predict whether these observations will be reproducible using larger-scale proteomic surveys. Pseudomonas aeruginosa, Enterobacter hormaechei and E. coli (known opportunistic pathogens in the urinary tract) produced proteins for iron/siderophore acquisition and high mobility (flagellins) in subjects 36 and 51 (Table 3). Flagellins are important for swarming and spread in the urogenital tract during infections [45]. The iron-sequestering protein lactotransferrin may be expressed and released into the urinary tract by neutrophils to sequester iron, which in turn, E. coli and E. hormaechei may counteract by addition of copious amounts of iron/siderophore receptors to the outer membrane protein repertoire (Table 3). While the S100-A8 and S100-A9 calprotectins appear to have numerous physiological functions, one of them is to sequester zinc in response to sensing the presence of bacterial pathogens. This sequestration can inactivate metalloproteinases that the bacteria produce and secrete to invade the host tissue [46]. Such proteases were not identified in our datasets.

Table 2. Bacterial profiles of urinary samples

Table 3. Human and bacterial proteins potentially contributing directly or indirectly to host-pathogen interactions in the urinary tract

Several other proteins implicated in the innate immune response were detected (Table 3). Calprotectin has acute and chronic pro-inflammatory functions and recruits immune cells to the site of inflammation [47]. Antimicrobial peroxidases, such as eosinophil and myeloperoxidase, produce reactive oxygen species during the respiratory burst of neutrophils and are directly microbicidal. Like calprotectin, annexin A1 is a calcium-binding protein that was detected in high abundance primarily in those donors where respiratory burst of neutrophils did seem to be muted (subjects 1, 16, 31, 34, 37, 39). These proteins are implicated in innate immunity, influence cell apoptosis and function as damage-associated molecular patterns (DAMPs) [48,49]. For a complete list of bacterial and human proteins profiled in this study, see Additional file 1.

Discussion

In this study we describe the healthy urine microbiome in a heterogeneous population of men and women with and without NB, using both 16S rDNA sequencing and metaproteomic analysis. Based on other reports [50-53] and including our data, this is further confirmation that the commonly held clinical belief that healthy urine should be sterile is false. Specifically, our data indicate that (1) when collected by the routine method of “clean, midstream catch”, healthy urine is not aseptic; (2) the healthy and NB urine microbiomes differs by gender; (3) the asymptomatic bacteriuria urine microbiome of people with NB differs from that of healthy controls; and (4) the asymptomatic bacteriuria urine microbiome of people with NB differs depending on duration of exposure to and type of urinary catheter.

This is the first report comparing the healthy urine microbiome in male and female subjects. Historically, and based on cultivation results, clinicians have assumed urine to be ‘sterile’. However, Wolfe et al. recently used 16S rDNA sequencing to identify uncultivated bacteria in the urine of healthy women [53]. Our data confirm these results in women, and further show that uncultivated bacteria are also present in the urine of healthy men. Moreover, we demonstrate that the healthy urine microbiome of males and females differs, with men having significantly more relative abundance of Corynebacterium, a common component of the superficial skin flora, and women having significantly greater relative abundance of Lactobacillus. Critical to this discussion is an understanding that because all samples from healthy subjects and those of subjects with NB who voided were collected by midstream clean catch, it is not possible to distinguish whether the microbes identified originated in the bladder, urethra, or both. Therefore, the possibility exists that the identified urine microbiomes are populated by species that exist in the bladder, urethra, or both.

Our finding of predominance of Corynebacterium in healthy males suggests that this microbe may contribute to the healthy urine microbiome. Not only was Corynebacterium identified to a significantly greater degree in healthy males compared to those with NB, males with NB who voided and were sampled by clean catch had the lowest abundance of this species. Dong et al. compared clean catch urine and distal urethral swabs in healthy volunteers and similarly found a predominance of Corynebacterium in both types of sampling but in significantly greater amounts in the clean catch urine samples [51]. Taken together, the data suggest that Corynebacterium may reside in the proximal urethra and/or bladder in addition to the distal urethra, and may play a role in the healthy urine microbiome.

While it is well established that in most healthy women of childbearing age the vaginal tract is colonized by Lactobacillus species [54-56], this has not been investigated in women with NB. We too found a clear preponderance of Lactobacillus in healthy control females; however, in addition to the greatest relative abundance of Lactobacillus in healthy control females, there was a progressive reduction in abundance of Lactobacillus in females with NB who void (clean catch sample) or who use intermittent catheterization, and females with NB who use indwelling (Foley) catheters (Figure 3 inset). This may suggest that either increasing exposure to a urinary catheter and/or increasing severity of NB can influence the ability of Lactobacillus to colonize the urinary tract. Alternatively, Lactobacillus may merely be a contaminant of the external urethra that arises from proximity to the vaginal flora. However, this appears less likely since Wolfe et al. showed the presence of Lactobacillus in urine collected by transurethral catheters and suprapubic aspirate, which samples the bladder directly [53]. Because lactic acid-producing Lactobacillus species contribute to controlling the growth of more virulent bacteria that cannot survive in a more acidic environment, the presence of Lactobacillus within the urethra and/or bladder may be protective. This has been shown to be the case in infants [57] and may also be true for males as Lactobacillus has been shown to be present in clean catch urine samples of healthy males by Dong et al. [51] and in our study. Our findings suggest that Lactobacillus may be a commensal organism present during states of health, more in females than in males, and that the microbiome of at-risk populations may be characterized by a distinct lack of Lactobacillus, which perhaps creates a better environment for the growth of pathogenic microorganisms. Together, these findings suggest that the clinical objective of ‘sterile’, microbe-free urine may not be optimal for the patient.

This is the first report of 16S rDNA sequencing of urine in people with NB, providing much more detail about the ABU state than has previously been available through cultivation-based evidence. Standard cultivation diagnostics show populations vulnerable to bacteriuria include nursing home residents utilizing long-term indwelling catheterization [58], institutionalized Veterans utilizing intermittent catheterization [59], and patients with SCI who utilize urinary catheters [60-62]. Our analysis not only confirms the cultivation-based evidence, but also shows that a significantly different microbiome was present in the NB group, and that 16S rDNA sequencing specifically identifies microorganisms, such as Enterobacteriales, as potential major contributors to a pathogenic microbiome. These results supplement those from Bank et al., where urine specimens of 142 consecutive patients with varied genitourinary pathology (kidney stones, indwelling catheters, or suspected UTI) were analyzed with standard cultivation and screened for Actinobaculum schaalii using PCR [50]. Those authors found that the most heavily colonized patients were those who were older and who utilized indwelling urinary catheters, while the younger patients who typically use intermittent catheterization (and may have utilized urinary catheters for a shorter amount of time) were colonized with bacteria to a lesser degree. The Enterobacteriaceae was most commonly isolated in the catheterized specimens in that study. Significant variance in medical comorbidities, underlying genitourinary pathology, and method of urine collection limit any further comparisons to that study. Moreover, because our NB population was reportedly asymptomatic (i.e. infection free), our findings demonstrate that in the catheterized population, the microbiome is intrinsically different than in controls, even in the absence of illness. This distinction is important, since in the clinical setting ABU is often inappropriately treated with antibiotics, which may further disrupt the NB microbiome.

The present data also indicate that the urine microbiome of healthy subjects with NB became increasingly abundant with Enterobacteriales with increasing duration of NB, whereas Lactobacillus decreased over time, both in men and women. While the urinary microbiome of men and women with NB remained similar to that of healthy controls during the first several months after NB diagnosis, by one year the urine microbiome was nearly devoid of Lactobacillus and dominated by Enterococcus. This further suggests a change in the microbiome with duration of NB that may place patients at increased risk of UTI.

Fundamental to these discoveries is the diverse sample population and our novel analytic approach of utilizing a combination of 16S rDNA sequencing in all subjects and metaproteomics in a subsample. Clinical gold standard diagnostic testing when a patient presents with signs and/or symptoms of UTI includes (1) urinalysis to confirm urinary tract inflammation and (2) urine culture to identify, quantify and predict antimicrobial resistance to a given pathogen(s). We believe that 16S rDNA sequencing has the potential for translation to the clinic, offering significant clinical advancement over diagnostic urine culture because it provides a greater depth of understanding and sensitivity pertaining to the composition of commensal and potentially pathogenic microbes present in urine. Furthermore, prospective assessments during varying periods of health and disease may allow personalization of care that has not been possible to date with our current diagnostic methods. Urinary metaproteomic profiles in parallel may contribute to the identification of a host inflammatory response utilizing urinary biomarkers with greater sensitivity and specificity for UTI than traditional measures of urinary leukocyte esterase production or white blood cell count detected by urinalysis. We hypothesize that protein profiles with distinct abundance ratios of immunomodulatory versus pro-inflammatory and microbiocidal molecules, are indicative of either UTI or reflect asymptomatic colonization. For instance, the presence of lactotransferrin in urine has been used to support the notion of a “battle for iron” being waged between the host and the pathogen, involving E. coli, particularly in the case of UTI by E. coli given the abundance of its triad of iron acquisition receptors [63]. While metaproteomics holds promise as a diagnostic method to discriminate between symptomatic UTI and ABU, more in-depth fractionation of samples is needed to evaluate whether this method can reach the sensitivity of 16S rDNA profiling methods. Larger patient cohorts, including those diagnosed with UTI symptoms, are required to assess if such ‘omics methods’ more accurately differentiate UTI from ABU than current diagnostic standards. In addition, they suggest new considerations that may impact preventive and treatment options for people at-risk for UTI.

Several limitations of this study are to be considered when interpreting the results. The major limiting factor is that some subjects (healthy controls and NB subjects who voided spontaneously) had urine collection via midstream clean catch sampling while the subjects with NB who managed their bladder with intermittent catheterization or indwelling Foley (urethral) catheters were sampled directly from the catheter. Therefore, microbes identified in the former groups could be representatives from the bladder, urethra, skin, or a combination of these microbiomes, whereas urine from subjects in the latter groups represents bladder microbiota. This distinction raises other questions, such as to what degree do differing microbiomes in the urethra and bladder influence each other, and does a changing urethral microbiome represent a potential antecedent UTI state. Further, given that clinical urine sampling is unlikely to become more invasive (via direct sampling from the bladder), how do clinicians interpret the presence of bacteria in clean catch midstream urine samples that could potentially originate from multiple anatomic locations? Lastly, future studies will employ a larger sample size as our results indicate that the urine microbiome differs by a number of clinical factors requiring multiple stratification points.

Conclusions

Contrary to clinical dogma that healthy urine is sterile; these results suggest that the state of healthy urine is, in fact, one of ‘asymptomatic bacteriuria’. Utilizing high throughput sequencing and metaproteomics, we have described the healthy urine microbiome of a number of populations: male and female healthy controls and healthy subjects with NB. Differing urinary microbiomes for males and females were described. We have demonstrated that NB and/or urinary catheterization impacts the healthy urine microbiome in both genders and this varies by type of bladder management and duration of NB. Furthermore, the presence of a variety of urine microbiomes differing on key, clinical characteristics suggests the benefit of a more personalized approach to UTI care. Clearly, DNA sequencing techniques allow for more specific assessment of the contributing microorganisms than do current clinical diagnostic standards, offering the potential for significant clinical advancement of diagnostic methods for UTI, which have otherwise remained relatively unchanged for decades. Longitudinal differentiation of the urine microbiome at the time of, prior to, and after infection also will be necessary to fully describe the course of disease and its antecedents. These findings advance clinical translational science toward improved diagnostics and more targeted use of therapeutics.

Abbreviations

ABU: Asymptomatic Bacteriuria; CAUTI: Catheter Associated Urinary Tract Infections; CNMC: Children’s National Medical Center; CRAN: Comprehensive R Archive Network; CTSA: Clinical and Translational Science Awards; DAMPs: Damage-Associated Molecular Patterns; DHHS: Department of Health and Human Services; DNA: Deoxyribonucleic Acid; Ec: Escherichia coli; EDTA: Ethylenediaminetetraacetic Acid; Ef: Enterococcus faecalis; Eh: Enterobacter hormaechei; ESBL: Extended-Spectrum-Beta-Lactamases; FASP: Filter-Aided Sample Preparation; FC: Foley Catheter; HC: Healthy Control; Hpf: High power field; IC: Intermittent Catheter; IRB: Institutional Review Board; JCVI: J. Craig Venter Institute; Kp: Klebsiella pneumoniae; LCA: Least Common Ancestor; Lj: Lactobacillus jensenii; LTQ: Linear Trap Quadrupole; NB: Neuropathic Bladder; NCATS: National Center for Advancing Translational Sciences; NCBI: National Center for Biotechnology Informatio; NCMRR: National Center for Medical Rehabilitation Research; NCRR: National Center for Research Resources; NEG: Negative; NIH: National Institutes of Health; NINDS: National Institute for Neurological Disorders and Stroke; NJ: Neighbor-Joining; No: Number; NRH: National Rehabilitation Hospital; NT: Nucleotide; OGE: Oracle Grid Engine; OTU: Operational Taxonomic Unit; Pa: Pseudomonas aeruginosa; PBS: Phosphate Buffered Saline; PCA: Principal Component Analysis; PCR: Polymerase Chain Reaction; Pm: Proteus mirabilis; Ps: Providencia stuartii; RDP: Ribosomal Database Project; RNA: Ribonucleic Acid; rRNA: Ribosomal Ribonucleic Acid; SCI: Spinal Cord Injury; SDS-PAGE: Sodium Dodecyl Sulfate Polyacrylamide Gel Electrophoresis; Sp: Streptococcus pneumoniae; SSU: Small Subunit; TRA: Trace; UTI: Urinary Tract Infection; WBC: White Blood Cells; YAP: Yet Another Pipeline.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

DEF, RP, SLG conceived and organized the study; MS, BS, SKL, MT performed laboratory experiments; DEF, RP, SS, HP, SK, MS, IL, BS, SLG processed and/or analyzed data; and DEF, RP, SS, HP, SK, MT, KEN, SLG wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank Shih-Ting Huang for her contributions to the MS database assembly and computational MS data processing and Anthony Scott Durkin for registering the NCBI project and for assistance with data submission. We are grateful to J. Craig Venter at JCVI and Eric Hoffman at CNMC for internal support of this project. This work was partially supported by NIH NCMRR/NINDS 2R24HD050846-06 (NCMRR-DC Core Molecular and Functional Outcome Measures in Rehabilitation Medicine). In addition, this project has been funded in whole or in part with Federal funds (Grant # UL1TR000101) from the National Center for Advancing Translational Sciences (NCATS), National Institutes of Health (NIH), through the Clinical and Translational Science Awards Program (CTSA), a trademark of DHHS, part of the Roadmap Initiative, “Re-Engineering the Clinical Research Enterprise.

References

  1. Infectious disease, chapter 7, urinary tract infections.

    http://pathmicro.med.sc.edu/infectiousdisease/UrinaryTractInfections.htm webcite

    OpenURL

  2. Foxman B, Barlow R, D'Arcy H, Gillespie B, Sobel JD: Urinary tract infection: self-reported incidence and associated costs.

    Ann Epidemiol 2000, 10:509-515. PubMed Abstract | Publisher Full Text OpenURL

  3. Griebling TL: Urologic diseases in america project: trends in resource use for urinary tract infections in men.

    J Urol 2005, 173:1288-1294. PubMed Abstract | Publisher Full Text OpenURL

  4. Litza JA, Brill JR: Urinary tract infections.

    Prim Care 2010, 37:491-507.

    viii

    PubMed Abstract | Publisher Full Text OpenURL

  5. HHS action plan to prevent healthcare-sssociated infections.

    http://www.hhs.gov/ash/initiatives/hai/actionplan webcite

    OpenURL

  6. Maki DG, Tambyah PA: Engineering out the risk for infection with urinary catheters.

    Emerg Infect Dis 2001, 7:342-347. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Salgado CD, Karchmer RB, Farr BM: Prevention of catheter-sssociated urinary tract infections. 4th edition. PA: Lippincott Williams & Wilkins, Philadelphia; 2003. OpenURL

  8. Macejko AM, Schaeffer AJ: Asymptomatic bacteriuria and symptomatic urinary tract infections during pregnancy.

    Urol Clin North Am 2007, 34:35-42. PubMed Abstract | Publisher Full Text OpenURL

  9. Wald HL, Kramer AM: Nonpayment for harms resulting from medical care: catheter-associated urinary tract infections.

    JAMA 2007, 298:2782-2784. PubMed Abstract | Publisher Full Text OpenURL

  10. Saint S, Meddings JA, Calfee D, Kowalski CP, Krein SL: Catheter-associated urinary tract infection and the medicare rule changes.

    Ann Intern Med 2009, 150:877-884. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Dicianno BE, Wilson R: Hospitalizations of adults with spina bifida and congenital spinal cord anomalies.

    Arch Phys Med Rehabil 2010, 91:529-535. PubMed Abstract | Publisher Full Text OpenURL

  12. Filler G, Gharib M, Casier S, Lodige P: Ehrich JH. Dave S: Prevention of chronic kidney disease in spina bifida. International urology and nephrology; 2011. OpenURL

  13. Salonen A, Nikkila J, Jalanka-Tuovinen J, Immonen O, Rajilic-Stojanovic M, Kekkonen RA, Palva A, de Vos WM: Comparative analysis of fecal DNA extraction methods with phylogenetic microarray: effective recovery of bacterial and archaeal DNA using mechanical cell lysis.

    J Microbiol Methods 2010, 81:127-134. PubMed Abstract | Publisher Full Text OpenURL

  14. Liu Z, DeSantis TZ, Andersen GL, Knight R: Accurate taxonomy assignments from 16S rRNA sequences produced by highly parallel pyrosequencers.

    Nucleic Acids Res 2008, 36:e120. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al.: Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities.

    Appl Environ Microbiol 2009, 75:7537-7541. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Chou HH, Holmes MH: DNA sequence quality trimming and vector removal.

    Bioinformatics 2001, 17:1093-1104. PubMed Abstract | Publisher Full Text OpenURL

  17. Engelbrektson A, Kunin V, Wrighton KC, Zvenigorodsky N, Chen F, Ochman H, Hugenholtz P: Experimental factors affecting PCR-based estimates of microbial species richness and evenness.

    ISME J 2010, 4:642-647. PubMed Abstract | Publisher Full Text OpenURL

  18. Kunin V, Engelbrektson A, Ochman H, Hugenholtz P: Wrinkles in the rare biosphere: pyrosequencing errors can lead to artificial inflation of diversity estimates.

    Environ Microbiol 2010, 12:118-123. PubMed Abstract | Publisher Full Text OpenURL

  19. Huang Y, Niu B, Gao Y, Fu L, Li W: CD-HIT Suite: a web server for clustering and comparing biological sequences.

    Bioinformatics 2010, 26:680-682. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Niu B, Fu L, Sun S, Li W: Artificial and natural duplicates in pyrosequencing reads of metagenomic data.

    BMC Bioinformatics 2010, 11:187. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  21. Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glockner FO: SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB.

    Nucleic Acids Res 2007, 35:7188-7196. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, Ciulla D, Tabbaa D, Highlander SK, Sodergren E, et al.: Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons.

    Genome Res 2011, 21:494-504. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. 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:D141-D145. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Hamady M, Knight R: Microbial community profiling for human microbiome projects: Tools, techniques, and challenges.

    Genome Res 2009, 19:1141-1152. PubMed Abstract | Publisher Full Text OpenURL

  25. The Comprehensive R Archive Network.

    http://cran.r-project.org/ webcite

    OpenURL

  26. YAP 16S sequence-processing pipeline.

    http://github.com/shpakoo/YAP webcite

    OpenURL

  27. Dray S, Dufour AB: The ade4 package: Implementing the duality diagram for ecologists.

    J Stat Softw 2007, 22:1-20. OpenURL

  28. Hollander M, Wolfe DA: Wiley series in probability and statistics. In Nonparametric Statistical Methods. 2nd edition. New York: John Wiley & Sons, Inc; 1999:36-51. OpenURL

  29. Ribosome Database Project.

    http://rdp.cme.msu.edu webcite

    OpenURL

  30. Sonnhammer EL, Hollich V: Scoredist: a simple and robust protein sequence distance estimator.

    BMC Bioinformatics 2005, 6:108. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  31. Seal BS, Fouts DE, Simmons M, Garrish JK, Kuntz RL, Woolsey R, Schegg KM, Kropinski AM, Ackermann HW, Siragusa GR: Clostridium perfringens bacteriophages PhiCP39O and PhiCP26F: genomic organization and proteomic analysis of the virions.

    Arch Virol 2011, 156:25-35. PubMed Abstract | Publisher Full Text OpenURL

  32. FigTree Grapical Viewer of Phylogenetic Trees.

    http://tree.bio.ed.ac.uk/software/figtree/ webcite

    OpenURL

  33. Wisniewski JR, Zougman A, Nagaraj N, Mann M: Universal sample preparation method for proteome analysis.

    Nat Methods 2009, 6:359-362. PubMed Abstract | Publisher Full Text OpenURL

  34. Pieper R, Zhang Q, Clark DJ, Huang ST, Suh MJ, Braisted JC, Payne SH, Fleischmann RD, Peterson SN, Tzipori S: Characterizing the Escherichia coli O157:H7 Proteome Including Protein Associations with Higher Order Assemblies.

    PLoS One 2011, 6:e26554. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Uniref90 Protein Sequence Database.

    http://www.ebi.ac.uk/uniprot/database/download.html webcite

    OpenURL

  36. Verberkmoes NC, Russell AL, Shah M, Godzik A, Rosenquist M, Halfvarson J, Lefsrud MG, Apajalahti J, Tysk C, Hettich RL, Jansson JK: Shotgun metaproteomics of the human distal gut microbiota.

    ISME J 2009, 3:179-189. PubMed Abstract | Publisher Full Text OpenURL

  37. Mascot Percolator online help. http://www.matrixscience.com/help/percolator_help.html

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

    Appl Environ Microb 2007, 73:5261-5267. Publisher Full Text OpenURL

  39. Domann E, Hong G, Imirzalioglu C, Turschner S, Kuhle J, Watzel C, Hain T, Hossain H, Chakraborty T: Culture-independent identification of pathogenic bacteria and polymicrobial infections in the genitourinary tract of renal transplant recipients.

    J Clin Microbiol 2003, 41:5500-5510. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Imirzalioglu C, Hain T, Chakraborty T, Domann E: Hidden pathogens uncovered: metagenomic analysis of urinary tract infections.

    Andrologia 2008, 40:66-71. PubMed Abstract | Publisher Full Text OpenURL

  41. Beighton D, Carr AD, Oppenheim BA: Identification of viridans streptococci associated with bacteraemia in neutropenic cancer patients.

    J Med Microbiol 1994, 40:202-204. PubMed Abstract | Publisher Full Text OpenURL

  42. Corredoira JC, Alonso MP, Garcia JF, Casariego E, Coira A, Rodriguez A, Pita J, Louzao C, Pombo B, Lopez MJ, Varela J: Clinical characteristics and significance of Streptococcus salivarius bacteremia and Streptococcus bovis bacteremia: a prospective 16-year study.

    Eur J Clin Microbiol Infect Dis 2005, 24:250-255. PubMed Abstract | Publisher Full Text OpenURL

  43. Ruoff KL, Miller SI, Garner CV, Ferraro MJ, Calderwood SB: Bacteremia with Streptococcus bovis and Streptococcus salivarius: clinical correlates of more accurate identification of isolates.

    J Clin Microbiol 1989, 27:305-308. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Hoecker JL, Pickering LK, Groschel D, Kohl S: Streptococcus salivarius sepsis in children with malignancies.

    J Pediatr 1978, 92:337-338. PubMed Abstract | Publisher Full Text OpenURL

  45. Lane MC, Alteri CJ, Smith SN, Mobley HL: Expression of flagella is coincident with uropathogenic Escherichia coli ascension to the upper urinary tract.

    Proc Natl Acad Sci USA 2007, 104:16669-16674. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Isaksen B, Fagerhol MK: Calprotectin inhibits matrix metalloproteinases by sequestration of zinc.

    Mol Pathol 2001, 54:289-292. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. Gebhardt C, Nemeth J, Angel P, Hess J: S100A8 and S100A9 in inflammation and cancer.

    Biochem Pharmacol 2006, 72:1622-1631. PubMed Abstract | Publisher Full Text OpenURL

  48. Bianchi ME: DAMPs, PAMPs and alarmins: all we need to know about danger.

    J Leukoc Biol 2007, 81:1-5. PubMed Abstract OpenURL

  49. Sato S, St-Pierre C, Bhaumik P, Nieminen J: Galectins in innate immunity: dual functions of host soluble beta-galactoside-binding lectins as damage-associated molecular patterns (DAMPs) and as receptors for pathogen-associated molecular patterns (PAMPs).

    Immunol Rev 2009, 230:172-187. PubMed Abstract | Publisher Full Text OpenURL

  50. Bank S, Hansen TM, KM SB, Lund L, Prag J: Actinobaculum schaalii in urological patients, screened with real-time polymerase chain reaction.

    Scand J Urol Nephrol 2011, 45:406-410. PubMed Abstract | Publisher Full Text OpenURL

  51. Dong Q, Nelson DE, Toh E, Diao L, Gao X, Fortenberry JD, Van der Pol B: The microbial communities in male first catch urine are highly similar to those in paired urethral swab specimens.

    PLoS One 2011, 6:e19709. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Siddiqui H, Nederbragt AJ, Lagesen K, Jeansson SL, Jakobsen KS: Assessing diversity of the female urine microbiota by high throughput sequencing of 16S rDNA amplicons.

    BMC Microbiol 2011, 11:244. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  53. Wolfe AJ, Toh E, Shibata N, Rong R, Kenton K, Fitzgerald M, Mueller ER, Schreckenberger P, Dong Q, Nelson DE, Brubaker L: Evidence of Uncultivated Bacteria in the Adult Female Bladder.

    J Clin Microbiol 2012, 50:1376-1383. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Pfau A, Sacks T: The bacterial flora of the vaginal vestibule, urethra and vagina in the normal premenopausal woman.

    J Urol 1977, 118:292-295. PubMed Abstract OpenURL

  55. Marrie TJ, Harding GK, Ronald AR: Anaerobic and aerobic urethral flora in healthy females.

    J Clin Microbiol 1978, 8:67-72. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  56. Zhang R, Daroczy K, Xiao B, Yu L, Chen R, Liao Q: Qualitative and semi-quantitative analysis ofLactobacillusspecies in vagina of healthy fertile and postmenopausal Chinese women.

    J Med Microbiol 2012, 61:729-739. PubMed Abstract | Publisher Full Text OpenURL

  57. Lee JW, Shim YH, Lee SJ: Lactobacillus colonization status in infants with urinary tract infection.

    Pediatr Nephrol 2009, 24:135-139. PubMed Abstract | Publisher Full Text OpenURL

  58. Warren JW, Tenney JH, Hoopes JM, Muncie HL, Anthony WC: A prospective microbiologic study of bacteriuria in patients with chronic indwelling urethral catheters.

    J Infect Dis 1982, 146:719-723. PubMed Abstract | Publisher Full Text OpenURL

  59. Terpenning MS, Allada R, Kauffman CA: Intermittent urethral catheterization in the elderly.

    J Am Geriatr Soc 1989, 37:411-416. PubMed Abstract OpenURL

  60. Cardenas DD, Moore KN, Dannels-McClure A, Scelza WM, Graves DE, Brooks M, Busch AK: Intermittent catheterization with a hydrophilic-coated catheter delays urinary tract infections in acute spinal cord injury: a prospective, randomized, multicenter trial.

    Pm R 2011, 3:408-417. PubMed Abstract | Publisher Full Text OpenURL

  61. Esclarin De Ruz A, Garcia Leoni E, Herruzo Cabrera R: Epidemiology and risk factors for urinary tract infection in patients with spinal cord injury.

    J Urol 2000, 164:1285-1289. PubMed Abstract | Publisher Full Text OpenURL

  62. Wyndaele JJ: Complications of intermittent catheterization: their prevention and treatment.

    Spinal Cord 2002, 40:536-541. PubMed Abstract | Publisher Full Text OpenURL

  63. Skaar EP: The battle for iron between bacterial pathogens and their vertebrate hosts.

    PLoS Pathog 2010, 6:e1000949. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL