- Open Access
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
Journal of Translational Medicine volume 10, Article number: 174 (2012)
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.
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.
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.
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.
Affecting nearly one half of all Americans over the course of a lifetime and with costs exceeding $1 billion annually, 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. 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. 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.
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). (Table1) 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.
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. 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.
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 (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. 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. 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. A module of CD-HIT suite called CD-HIT-EST was employed to perform species-level read-clustering for subsequent analyses.
The orchestration and automation of steps has been achieved using a custom set of in-house utilities written in python and R programming languages. These utilities are available online as a part of the YAP package on github. 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. Further statistical analyses were accomplished using R. Heat maps were generated using the heatmap.2 function of the gplots package available on CRAN. 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. 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.
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. 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. 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, and edited using Fig Tree. 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.
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). 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. 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 file1) and Uniref90’s human database subset. 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. As described in Additional file1, 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.
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) (Table1). 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 (Additional file2) where taxonomy of species-level operational taxonomic units (OTUs) was determined using the RDP classifier. Analysis of the bacterial community in these urine samples revealed between 5 and 236 species-level (97% identity) OTUs per individual (Additional file3). Per sample taxonomic profiles were generated, showing considerable sample-to-sample variation (Additional file4). 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 (Figure1A). 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 (Figure1A). The Lactobacillales and Enterobacteriales were the two most relative abundant and changing taxonomic groups.
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 (Figure1B). 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 (Figure1B). 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 (Figure1B). 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 (Figure1C). 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 Figure1C). Lactobacillus, Klebsiella, Corynebacterium, Staphylococcus, Streptococcus, Gardnerella, Prevotella, Escherichia and Enterococcus defined statistically significant differences between healthy bladder controls and NB (“*” in Figure1C).
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 (Figure2). 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, Figure2).
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 (Figure3). 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 (Figure3). The preponderance of Enterococcus sp. in male NB observed in Figure1C was primarily within the patients with NB who emptied by spontaneous voiding (Figure3).
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 (Figure4). 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, Figure4). 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 (Figure4), 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.
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 (Figure5). 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.
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 (Table2). 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 (Table3). Flagellins are important for swarming and spread in the urogenital tract during infections. 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 (Table3). 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. Such proteases were not identified in our datasets.
Several other proteins implicated in the innate immune response were detected (Table3). Calprotectin has acute and chronic pro-inflammatory functions and recruits immune cells to the site of inflammation. 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 file1.
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. 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. 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 (Figure3 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. 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 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. 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, institutionalized Veterans utilizing intermittent catheterization, 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. 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. 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.
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.
Catheter Associated Urinary Tract Infections
Children’s National Medical Center
Comprehensive R Archive Network
Clinical and Translational Science Awards
Damage-Associated Molecular Patterns
Department of Health and Human Services
Filter-Aided Sample Preparation
High power field
Institutional Review Board
J. Craig Venter Institute
Least Common Ancestor
Linear Trap Quadrupole
National Center for Advancing Translational Sciences
National Center for Biotechnology Informatio
National Center for Medical Rehabilitation Research
National Center for Research Resources
National Institutes of Health
National Institute for Neurological Disorders and Stroke
National Rehabilitation Hospital
Oracle Grid Engine
Operational Taxonomic Unit
Phosphate Buffered Saline
Principal Component Analysis
Polymerase Chain Reaction
Ribosomal Database Project
Ribosomal Ribonucleic Acid
Spinal Cord Injury
Sodium Dodecyl Sulfate Polyacrylamide Gel Electrophoresis
Urinary Tract Infection
White Blood Cells
Yet Another Pipeline.
Infectious disease, chapter 7, urinary tract infections.http://pathmicro.med.sc.edu/infectiousdisease/UrinaryTractInfections.htm,
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. 10.1016/S1047-2797(00)00072-7.
Griebling TL: Urologic diseases in america project: trends in resource use for urinary tract infections in men. J Urol. 2005, 173: 1288-1294. 10.1097/01.ju.0000155595.98120.8e.
Litza JA, Brill JR: Urinary tract infections. Prim Care. 2010, 37: 491-507. 10.1016/j.pop.2010.04.001. viii
HHS action plan to prevent healthcare-sssociated infections.http://www.hhs.gov/ash/initiatives/hai/actionplan,
Maki DG, Tambyah PA: Engineering out the risk for infection with urinary catheters. Emerg Infect Dis. 2001, 7: 342-347. 10.3201/eid0702.010240.
Salgado CD, Karchmer RB, Farr BM: Prevention of catheter-sssociated urinary tract infections. 2003, PA: Lippincott Williams & Wilkins, Philadelphia, 4
Macejko AM, Schaeffer AJ: Asymptomatic bacteriuria and symptomatic urinary tract infections during pregnancy. Urol Clin North Am. 2007, 34: 35-42. 10.1016/j.ucl.2006.10.010.
Wald HL, Kramer AM: Nonpayment for harms resulting from medical care: catheter-associated urinary tract infections. JAMA. 2007, 298: 2782-2784. 10.1001/jama.298.23.2782.
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.
Dicianno BE, Wilson R: Hospitalizations of adults with spina bifida and congenital spinal cord anomalies. Arch Phys Med Rehabil. 2010, 91: 529-535. 10.1016/j.apmr.2009.11.023.
Filler G, Gharib M, Casier S, Lodige P: Ehrich JH. 2011, Dave S: Prevention of chronic kidney disease in spina bifida. International urology and nephrology
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. 10.1016/j.mimet.2010.02.007.
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-10.1093/nar/gkn491.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ: Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009, 75: 7537-7541. 10.1128/AEM.01541-09.
Chou HH, Holmes MH: DNA sequence quality trimming and vector removal. Bioinformatics. 2001, 17: 1093-1104. 10.1093/bioinformatics/17.12.1093.
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. 10.1038/ismej.2009.153.
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. 10.1111/j.1462-2920.2009.02051.x.
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. 10.1093/bioinformatics/btq003.
Niu B, Fu L, Sun S, Li W: Artificial and natural duplicates in pyrosequencing reads of metagenomic data. BMC Bioinformatics. 2010, 11: 187-10.1186/1471-2105-11-187.
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. 10.1093/nar/gkm864.
Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, Ciulla D, Tabbaa D, Highlander SK, Sodergren E: Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 2011, 21: 494-504. 10.1101/gr.112730.110.
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. 10.1093/nar/gkn879.
Hamady M, Knight R: Microbial community profiling for human microbiome projects: Tools, techniques, and challenges. Genome Res. 2009, 19: 1141-1152. 10.1101/gr.085464.108.
The Comprehensive R Archive Network.http://cran.r-project.org/,
YAP 16S sequence-processing pipeline.http://github.com/shpakoo/YAP,
Dray S, Dufour AB: The ade4 package: Implementing the duality diagram for ecologists. J Stat Softw. 2007, 22: 1-20.
Hollander M, Wolfe DA: Wiley series in probability and statistics. Nonparametric Statistical Methods. 1999, New York: John Wiley & Sons, Inc, 36-51. 2
Ribosome Database Project.http://rdp.cme.msu.edu,
Sonnhammer EL, Hollich V: Scoredist: a simple and robust protein sequence distance estimator. BMC Bioinformatics. 2005, 6: 108-10.1186/1471-2105-6-108.
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. 10.1007/s00705-010-0812-z.
FigTree Grapical Viewer of Phylogenetic Trees.http://tree.bio.ed.ac.uk/software/figtree/,
Wisniewski JR, Zougman A, Nagaraj N, Mann M: Universal sample preparation method for proteome analysis. Nat Methods. 2009, 6: 359-362. 10.1038/nmeth.1322.
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-10.1371/journal.pone.0026554.
Uniref90 Protein Sequence Database.http://www.ebi.ac.uk/uniprot/database/download.html,
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. 10.1038/ismej.2008.108.
Mascot Percolator online help.http://www.matrixscience.com/help/percolator_help.html,
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. 10.1128/AEM.00062-07.
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. 10.1128/JCM.41.12.5500-5510.2003.
Imirzalioglu C, Hain T, Chakraborty T, Domann E: Hidden pathogens uncovered: metagenomic analysis of urinary tract infections. Andrologia. 2008, 40: 66-71. 10.1111/j.1439-0272.2007.00830.x.
Beighton D, Carr AD, Oppenheim BA: Identification of viridans streptococci associated with bacteraemia in neutropenic cancer patients. J Med Microbiol. 1994, 40: 202-204. 10.1099/00222615-40-3-202.
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. 10.1007/s10096-005-1314-x.
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.
Hoecker JL, Pickering LK, Groschel D, Kohl S: Streptococcus salivarius sepsis in children with malignancies. J Pediatr. 1978, 92: 337-338. 10.1016/S0022-3476(78)80040-7.
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. 10.1073/pnas.0607898104.
Isaksen B, Fagerhol MK: Calprotectin inhibits matrix metalloproteinases by sequestration of zinc. Mol Pathol. 2001, 54: 289-292. 10.1136/mp.54.5.289.
Gebhardt C, Nemeth J, Angel P, Hess J: S100A8 and S100A9 in inflammation and cancer. Biochem Pharmacol. 2006, 72: 1622-1631. 10.1016/j.bcp.2006.05.017.
Bianchi ME: DAMPs, PAMPs and alarmins: all we need to know about danger. J Leukoc Biol. 2007, 81: 1-5.
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. 10.1111/j.1600-065X.2009.00790.x.
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. 10.3109/00365599.2011.599333.
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-10.1371/journal.pone.0019709.
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-10.1186/1471-2180-11-244.
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. 10.1128/JCM.05852-11.
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.
Marrie TJ, Harding GK, Ronald AR: Anaerobic and aerobic urethral flora in healthy females. J Clin Microbiol. 1978, 8: 67-72.
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. 10.1099/jmm.0.038687-0.
Lee JW, Shim YH, Lee SJ: Lactobacillus colonization status in infants with urinary tract infection. Pediatr Nephrol. 2009, 24: 135-139. 10.1007/s00467-008-0974-z.
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. 10.1093/infdis/146.6.719.
Terpenning MS, Allada R, Kauffman CA: Intermittent urethral catheterization in the elderly. J Am Geriatr Soc. 1989, 37: 411-416.
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. 10.1016/j.pmrj.2011.01.001.
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. 10.1016/S0022-5347(05)67157-1.
Wyndaele JJ: Complications of intermittent catheterization: their prevention and treatment. Spinal Cord. 2002, 40: 536-541. 10.1038/sj.sc.3101348.
Skaar EP: The battle for iron between bacterial pathogens and their vertebrate hosts. PLoS Pathog. 2010, 6: e1000949-10.1371/journal.ppat.1000949.
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.
The authors declare that they have no competing interests.
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.
Electronic supplementary material
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.(PDF 2 MB)
About this article
Cite this article
Fouts, D.E., Pieper, R., Szpakowski, S. et al. 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. J Transl Med 10, 174 (2012). https://doi.org/10.1186/1479-5876-10-174