Impact Factor 3.518 | CiteScore 3.62
More on impact ›

Original Research ARTICLE

Front. Cell. Infect. Microbiol., 08 April 2019 | https://doi.org/10.3389/fcimb.2019.00092

The Vaginal and Urinary Microbiomes in Premenopausal Women With Interstitial Cystitis/Bladder Pain Syndrome as Compared to Unaffected Controls: A Pilot Cross-Sectional Study

Kate V. Meriwether1*, Zhenmin Lei1, Rajbir Singh2, Jeremy Gaskins3, Deslyn T. G. Hobson1 and Venkatakrishna Jala2
  • 1Department of Obstetrics and Gynecology, University of Louisville, Louisville, KY, United States
  • 2Department of Microbiology and Immunology, James Graham Brown Cancer Center, University of Louisville, Louisville, KY, United States
  • 3Department of Bioinformatics and Biostatistics, University of Louisville, Louisville, KY, United States

Interstitial cystitis/bladder pain syndrome (ICBPS) may be related to an altered genitourinary microbiome. Our aim was to compare the vaginal and urinary microbiomes between premenopausal women with ICBPS and unaffected controls. This cross-sectional study screened premenopausal women with an O'Leary-Sant questionnaire (ICBPS if score ≥6 on either index; controls <6 on both). Women completed questionnaires on health characteristics, pelvic floor symptoms (OABq, PFDI-20), body image (mBIS), and sexual function (PISQ-IR). Bacterial genomic DNA was isolated from vaginal and clean-catch urinary specimens; the bacterial 16 rRNA gene was sequenced and analyzed using the QIIME pipeline. We performed UniFrac analysis (β-diversity) and generated Chao1 estimator (richness) and Simpson index (richness and evenness) values. We analyzed 23 ICBPS and 18 non-ICBPS patients. ICBPS patients had increased vaginal deliveries, BMI, and public insurance as well as worsened OAB-q, PFDI-20, mBIS, and PISQ-IR domain scores. Lactobacilli was the most abundant genus in both cohorts, and anaerobic or fastidious predominance was similar between groups (p = 0.99). For both the urine and vagina specimens, Chao1 and Simpson indices were similar between ICBPS and unaffected women. Weighted and unweighted UniFrac analyses showed no differences between groups. A significant correlation existed between the urinary and vaginal Simpson indices in ICBPS women, but not in unaffected women. Premenopausal women with ICBPS, despite worsened socioeconomic indicators and pelvic floor function, were not found to have significantly different urinary and vaginal microbiomes compared to women without ICBPS.

Introduction

Interstitial cystitis/bladder pain syndrome (ICBPS) is a devastating and common disorder in women (Clemens et al., 2005). The pathophysiology of ICBPS is poorly understood (Patnaik et al., 2017), leading to an ineffectual treatment environment (Kanter et al., 2017). There is a need to find underlying etiologies and direct that knowledge at therapy (Hanno et al., 2011).

Microbial studies indicate that women with ICBPS may have urinary microbiome disturbances (Siddiqui et al., 2012; Nickel et al., 2016; Abernethy et al., 2017), and it is known that disturbances in the vaginal microbiome can alter the urinary microbes (Gottschick et al., 2017; Thomas-White et al., 2018). As ICBPS onset is often in the reproductive years (Clemens et al., 2012), microbial interactions with ICBPS in pre-menopausal women is an important focus of study. However, no published research integrates both urinary and vaginal niches into the microbial study of ICBPS pathology.

Our primary aim was to compare the urinary and vaginal microbiomes between two groups: pre-menopausal women with ICBPS and those women unaffected by ICBPS. Based on past study that women with urinary symptoms more altered species compared to healthy women (Pearce et al., 2014; Karstens et al., 2016), we hypothesized that women with ICBPS would have a more anaerobic and fastidious bacterial milieu than unaffected controls.

Materials and Methods

Premenopausal women were recruited from a subspecialty practice at the University of Louisville and underwent screening assessment with the O'Leary-Sant questionnaire (Lubeck et al., 2001). Women were eligible if they were hormonally pre-menopausal and excluded if they had issues that disturbed the microbiome of the genitourinary system, were hormonally post-menopausal, or unable to understand the consent process (Appendix 1).

Women who scored <6 points on both indices of the O'Leary-Sant questionnaire were considered unaffected, non-ICBPS participants, and women with a score ≥6 on either index were considered ICBPS participants. The study was approved by the Institutional Review Board, and women underwent an informed consent prior to participation. All participants had a midstream, clean-catch urinary specimen sent for microscopic urinalysis (mUA) with reflex urine culture if indicated. The urinary specimen was collected after a front to back wipe with a moistened gynecological cleaning towelette, with the patient instructed to void for 1–2 s prior to collection in the sterile cup. If the women were found to have urinary infection or other exclusion criteria after the study visit (Appendix 1), they were excluded.

Participating patients completed questionnaires on their characteristics, health history, use of tobacco, vaginal medications, hormones/contraceptives, or probiotics. If possible, we calculated the days since their last menstrual cycle (LMP). Patients completed the pelvic floor distress inventory (PFDI-20), overactive bladder questionnaire (OAB-q), modified body image scale (mBIS), and pelvic organ prolapse/incontinence sexual questionnaire: IUGA revised (PISQ-IR). Women underwent collection of a BD ProbeTecTM swab (Becton Dickinson, Franklin Lakes, NJ) from the mid-vaginal walls and a clean-catch, midstream urinary specimen in a sterile cup (separate from mUA specimen). Vaginal and urinary samples were both sealed and stored in refrigeration (4°C) for <12 h prior to laboratory processing.

At the laboratory, urine specimens were sampled after agitation by pipetting three aliquots from each specimen (1.0 mL each) into three 1.5 mL nucleic acid-free tubes. Tubes were centrifuged at 12,000 g for 10 min and the supernatant discarded (Pearce et al., 2014; Abernethy et al., 2017). The pellet was suspended in a PBS buffer, centrifuged again for 5 min at 12,000 g, and the remaining supernatant was removed.

At the laboratory, vaginal samples were also placed in −80°C until pellet extraction, where thawed swabs were agitated for 30 s in a tube containing 1 mL of sterile PBS and pressed against wall of the tube multiple times. The pellet was isolated by centrifuging the solution at 12,000 g for 5 min, and the supernatant was discarded.

Urine and vaginal pellets were stored at −80°C until DNA extraction. A bead based extraction method was used for isolation of the total bacterial genomic DNA (gDNA) from the pellets using a QIAmp BiOstic Bacteremia DNA kit (Qiagen Inc, Germantown, MD).

The hypervariable V3-V4 regions of the 16S rRNA gene were amplified using preexisting primers (Integrated DNA Technologies, Skokie, Illinois) as follows:

Forward Primer: 5′TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG3′

Reverse Primer: 5′GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGACTACHVGGGTATCTAATCC3′

We utilized 12.5 ngmicrobial genomic DNA with appropriate PCR conditions, and an extraction-negative controlwas run with amplification and sequencing steps for vaginal and urinary specimens to ensure lack of laboratory contamination. Each PCR reaction was performed in triplicate and then pooled with respective samples. The amplified bacterial gDNA wasquantified using NanoDrop8000(ThermoScientific, Waltham, MA). ThePCR products were cleaned using AMPure XP beads. These PCR products weresubjected to index PCR to attach dual indices and Illumina (San Diego, CA) sequencing adapters usingthe Nextera XT kit (FC-121-1012). The resulting products were again cleaned using AMPure XP beads and quantified by usingQubit dsDNA BR Assay kit (ThermoScientific, Waltham, MA). The samples (4 nM) were pooled and denatured to perform sequencing usingMiSeq protocols alongside a control library. The sequencing was performed usingthe Illumina MiSeq Reagents kit v3(600 cycles; MS-102-3003) on the MiSeqinstrument (Illumina, San Diego, CA).

Microbiome data were pre-processed with Microbiome Helper v1.0.2 (Comeau et al., 2017). Reads were stitched and filtered for quality reads. Paired-end reads were stitched using PEAR v0.0.10 (Zhang et al., 2014). Reads with >10% of bases having a quality score of ≤30 were removed using FASTX-Toolkit. Reads ≥400 base-pairs that contained the matching primer sequences were retained. Following this, chimeric sequences were removed using the implementation of the UCHIME algorithm (Edgar et al., 2011) in VSEARCH v2.4.4 (Rognes et al., 2016).

Sequences were grouped into Operational Taxonomic Units (OTU's) using QIIME v1.9.1 (Navas-Molina et al., 2013). We grouped reads into OTUs at 97% similarity using an open-reference approach in QIIME. Reads were first clustered with the algorithm SortMeRNA (Kopylova et al., 2012) against the Greengenes database (v13_8) (DeSantis et al., 2006). This was followed by de novo clustering with SUMACLUST (v1.0.00) (Mercier and Coissac, 2013). The OTUs were normalized by random subsampling (rarefying) to 5,000 and 10,000 reads per sample. We analyzed the data that was rarefied to a depth of 5,000 sequences/sample, as 10,000 sequences/sample was not found more informative and eliminated viable data (Figure 2).

The primary outcome was the presence of an anaerobic or fastidious organisms as the predominant genera of the sample (>50% of reads), based on past study indicating an increase in these microorganisms to the urinary microbiome of women with urinary disease (Pearce et al., 2014; Karstens et al., 2016). We defined anaerobic or fastidious bacteria as those including Prevotella, Mobiculus, Gardnerella vaginalis, Mycoplasma, Ureaplasma, Bacteroides, and aerobic bacteria that are not classically cultivated in culture1. This was a pilot study without prior data on the topic for power calculation; we aimed to recruit 20 patients per group for a convenience sample.

We calculated two measures of alpha diversity for each sample (urine and vagina of each patient) using QIIME: the Chao1 estimator (richness) and the Simpson index (evenness and richness). The McIntosh index of evenness and the Shannon index (evenness and richness) were also calculated for each sample as additional measures of alpha diversity. Beta diversity principle co-ordinate plots (PCOs) were generated using phylogenetic UniFrac distances (Lozupone and Knight, 2005). Permanova analyses (weighted and unweighted) calculated the significance of Unifrac distances between the two groups. The percentile of reads of a certain genus for the sample (out of the total reads in the sample) was calculated for the top 5 most predominant genera, as beyond the top 5 genera the percentiles were so minute (<1%) that comparisons between them would not provide meaningful data.

To compare between ICBPS and control groups, we utilized chi-square (Fisher's test for n<5) for categorical variables and t-tests for continuous, parametric variables (Mann-Whitney for non-parametric). The multivariate Hotelling's t-test was used to compare measures of alpha diversity (Chao1 estimator and Simpson index of both urine and vaginal samples) between the groups. If data for the Chao1 estimator and Simpson index were skewed (violated parametric assumptions), these were analyzed on the logarithmic scale. The ANOSIM (analysis of similarity) test was utilized on the data as well without logarithmic transformation to account for the non-parametric nature of these outcomes.

Kendall correlation measures were used to analyze interactions between different sample locations (urine and vagina) for alpha diversity measures within the separate groups (ICBPS and unaffected). The relationships between baseline patient characteristics or questionnaire scores (OABq, PFDI, mBIS, and PISQ-IR) and microbial diversity were analyzed utilizing marginal univariate regression models, considering the pairs of urinary and vaginal alpha measures separately due to low power. To correct for multiple comparisons, we determined significance by applying an FDR correction across all regressions with the same measure (Chao or Simpson) and sample site (urinary or vagina). For covariates affecting the alpha measures at the 0.20 FDR level, we prepared to create a multiple regression model to explore the effect of group (ICBPS or unaffected) controlling for discovered confounders. All statistical analyses were performed using the R statistical software, version 3.3.1.

Stacked bar plots for the top 10 genus classifications for the urinary and vaginal samples of ICBPS women and non-ICBPS women were created using cluster 3.0 following Euclidean distance hierarchical matrices. Images were generated with Java Tree Software2.

Results

We screened 136 women for participation and 46 patients enrolled (Figure 1). Of women who enrolled, their presenting pelvic floor complaints included urinary incontinence (n = 12 ICBPS; n = 5 control), pelvic pain (n = 8 ICBPS; n = 4 control), vaginal prolapse (n = 5 ICBPS; n = 2 control), recurrent urinary infections (n = 2 ICBPS; n = 0 control), and other pelvic complaints (n = 0 ICBPS; n = 4 control)., Four women without pelvic floor complaints volunteered from the community and were enrolled in the control cohort. Five women were excluded after enrollment due to urinary tract infections (n = 2) or other exclusion criteria discovered after the study visit (n = 3), leaving 41 women (23 ICBPS; 18 unaffected) in the analyses. One woman in the non-ICBPS group had an uninterpretable vaginal swab and urine microbiome measures that were extreme outliers, so her samples were excluded from alpha and beta diversity analyses.

FIGURE 1
www.frontiersin.org

Figure 1. Summary of screening and inclusion of study population.

The mean days since the LMP was 22.2 days in the ICBPS group and 18.3 days in the unaffected group. However, this was based on only 12 women in each group, due lack of data on LMP (n = 1 ICBPS; n = 4 control) or with past procedures (hysterectomy, endometrial ablation, or hormonal implant) that precluded menstruation (n = 10 ICBPS; n = 2 control) in other patients. This variable of mean days since LMP did not have a significant impact on any microbiome outcomes (p-values 0.39–0.72).

The ICBPS patients and non-ICBPS patients had similar age, racial distribution, smoking status, sexual activity, hormone use, history of pelvic surgery/treatment, and Charleston health index scores (Table 1). However, ICBPS patients had increased vaginal deliveries, BMI, and public insurance as well as worsened scores on the OAB-q, PFDI-20, mBIS, and several PISQ-IR domains in sexually active women.

TABLE 1
www.frontiersin.org

Table 1. Patient characteristics in the study groups and relationship of characteristics and outcomes to the vaginal and urinary microbiome.

We obtained 5,755,661 total quality sequence reads from 81 samples (average 71,058 reads/sample; 53,387 for urinary specimens and 89,170 for vaginal specimens). There were significantly more reads/sample in the vaginal space as opposed to the urinary space in this population (p<0.01). The average number of reads/sample for ICBPS women was 75,794, and for normal women was 64,832, which was not different between the groups (p = 0.31). The average reads/sample for urinary samples was not different between groups (ICBPS 56,070 and non-ICBPS 49,959, p = 0.62), and the average reads/sample for vaginal samples was also similar between groups (ICBPS 95,18 and non-ICBPS 80,580, p = 0.36).

One vaginal and 3 urine pellet samples contained <10,000 reads, and increasing rarefaction to 10,000 eliminated these data without substantially changing results (Figure 2), so rarefaction at 5,000 reads/sample was maintained. Urinary and vaginal microbiomes for observed species was similar between ICBPS women to non-ICBPS women, but urinary and vaginal samples were distinct from one another in both groups with rarefaction up to 5,000 reads/sample, and this did not change at up to 10,000 reads per sample (Figure 2). We did not detect significant differences between ICBPS and non-ICBPS women in beta diversity (multivariate Hotelling test, p = 0.58; ANOSIM p = 0.90). UniFrac analyses, both weighted and unweighted, yielded no significance in beta diversity between the ICBPS and control groups for the urine or the vagina (all p > 0.05, Figure 3).

FIGURE 2
www.frontiersin.org

Figure 2. Plotting of rarefaction measure (y-axis) against sequences per sample (x-axis), either for observed species (top graphs) or the entire phylogenetic tree (bottom graphs), for each sample location (urinary or vaginal) and each group (ICBPS or non-ICBPS women) up to 5,000 sequences per sample (left graphs) or up to 10,000 sequences per sample (right graphs).

FIGURE 3
www.frontiersin.org

Figure 3. Unweighted and weighted UniFrac analyses with 3D plotting by primary component (PC) distances of normal vaginal (yellow) and normal urine (purple) samples as well as ICBPS women vaginal (blue) and urine samples (red).

Anaerobic predominance in the microbiome (the primary outcome) was rare, with only 12% (5/41; n = 2 ICBPS patients, n = 3 unaffected patients; n = 2 ICBPS) of women in the study having an anaerobic or fastidious genus as the primary contributor (>50% of reads) to either the vaginal or urinary microbiome. Lactobacilli was the dominant genus in the majority of urine (40/41, 98%) and vaginal (36/40, 85%) samples. The frequency of Lactobacilli dominance between urinary and vaginal samples trended toward the urine have more Lactobacilli dominance, albeit not statistically significant (p = 0.057). For vaginal samples, 2/23 ICBPS patients (9%) and 2/17 (12%) non-ICBPS patients (p = 0.99) had predominance of a non-Lactobacilli genus. For these 4 patients where Lactobacilli did not dominate the vaginal microbiome, dominant genera included Prevotella (n = 2 unaffected patients; n = 1 ICBPS) or Shuttleworthia (n = 1 ICBPS). For the urine, no ICBPS patients and only one non-ICBPS patient (1/18, 6%) had dominance of a non-Lactobacilli genus (44% Prevotella). Stacked bar plots of the vagina (Figure 4A) and the urinary samples (Figure 4B) are shown. At the OTU level, microbial species were also similar between groups for urine and vagina (Supplementary Table 1 and Supplementary Figure 1).

FIGURE 4
www.frontiersin.org

Figure 4. Stacked bar plots to the genus level of the top 10 genera for the vagina (A) and urinary samples (B) for the entire study population, with all genera not in the top 10 listed as “other.” Samples are sorted into unaffected (left) and ICBPS women (right) and labeled below with the individual woman's study identification number. Each color corresponds to a different genus, with Lactobacilli (orange) being predominant in most specimens.

In comparing measures of alpha diversity between ICBPS and non-ICBPS groups on the logarithmic scale, no significant differences were seen in the urinary (p = 0.11) or vaginal Simpson indices (p = 0.70). For the Chao1 estimator of richness, the urine (p = 0.91) and the vagina (p = 0.64) were also similar. Shannon index and McIntosh evenness comparison between groups also showed no significant differences (Supplementary Table 2). ICBPS specimens demonstrated Kendall correlation between the vaginal and urinary Simpson indices (tau = 0.49, 95% CI 0.23 to 0.68, p < 0.01, but this was not seen for non-ICBPS women (tau = 0.27, 95% CI−0.03 to 0.58, p = 0.06). No other significant correlations were found between the urine and vagina for Chao1 or Simpson indices in either group.

Without FDR correction, we detected no significant relationships (Table 1) between the urinary Simpson index and the factors of smoking status (p = 0.11), OABq symptom index (p = 0.19) and quality of life (p = 0.20) scores, PISQ-IR condition impact score (p = 0.07), and mBIS score (p = 0.06). For vaginal samples, again without FDR correction, the PISQ-IR condition impact score (p = 0.15) and the mBIS score (p = 0.06) were non-significantly related to the Simpson index, and the CRADI-8 (p = 0.04) score was related to the vaginal Chao1 estimator (Table 1). However, after applying the FDR correction for marginal regression models at each site, none of these effects were significant at either the 0.05 or 0.20 FDR level.

Discussion

In this cross-sectional pilot study of the urinary and vaginal microbiome comparing premenopausal women with ICBPS to those unaffected by ICBPS, we did not detect any significant differences in the vaginal or urinary microbiome. We did, however, note severe discrepancies in patient characteristics and pelvic and sexual function between the groups. These results indicate that ICBPS women have many risk factors for compromised pelvic health, and ICBPS is not a microbiome pathology.

Of note, the correlation between the diversity (Simpson index) of the vaginal and urinary samples was stronger in the ICBPS group than in unaffected women despite the use of clean-catch urinary specimens. ICBPS women are known to have altered epithelial function (Liu et al., 2012), so it is possible that this finding may indicate some difference in the natural barrier between the urinary and vaginal spaces in ICBPS women that this study was not designed to determine.

Unfortunately, this study lacks data indicating that alpha diversity correlates with validated urinary symptom scores, so ability to channel this finding at an intervention for ICBPS is limited. However, given that no prior study has explored the vaginal microenvironment in ICBPS, this data is novel and suggests a hypothesis that ICBPS allows more communication between the urinary and vaginal microbiomes.

This study demonstrates that the urinary and vaginal microbiomes in pre-menopausal women seen to have an overwhelming predominance of Lactobacilli, regardless of the presence of ICBPS. It is known that urine is not sterile and holds microbes that correspond greatly to the vaginal microbiome (Wolfe et al., 2012), but urine still holds less bacterial genomic material than the genitalia (Gottschick et al., 2017). Lactobacilli are associated with health of the vaginal microbiome (Stapleton, 2016; Lewis et al., 2017) and are predominant in the premenopausal urinary microbiome as well (Thomas-White et al., 2017). It appears that the vast majority of bacteria communicating between the urine from the vaginal environments are Lactobacilli, which is both statistically probable and physiologically sensible (the urethral epithelium has inherent defenses against more pathologic bacteria) (Stapleton, 2016). Nevertheless, the Lactobacilli composition in the urine or vagina was similar between ICBPS and non-ICBPS groups in this detailed study, so we cannot conclude that vaginal Lactobacilli disturbance is related to ICBPS pathology.

Strengths of the study include the use of detailed, reliable methods to determine the microbial composition of the genitourinary system. Genetic analyses are known to be more reliable than simple culture or enzyme-testing, particularly in urinary samples (Wolfe et al., 2012; Mouraviev and McDonald, 2018). This study also included the vaginal microbiomes in the analyses, as opposed to former studies that only investigate the urinary microbiome (Pearce et al., 2014; Thomas-White et al., 2017). As the urine and vaginal microbiomes are likely to be closely related, ignoring the vagina may fail to capture the entire pathophysiology. Finally, this study collected and investigated multiple possible confounders of the genitourinary microbiome, an intricate system that cannot be considered without patient factors.

Limitations of this study are mostly due to the smaller size of the study population, limiting power and ability to analyze subpopulations. The vaginal microbial environment in the premenopausal woman is extremely complex and changes with factors such as diet, the menstrual cycle, and vaginal medication use (Ravel et al., 2011; Nunn and Forney, 2016). We attempted to account for the days since a woman's LMP, but many women lacked these data, and the cycle day may be inaccurate or not correlate to the hormonal phase. Another limitation was our use of clean-catch urine. Past studies have found microbiome differences between clean-catch and catheterized urine specimens (Wolfe et al., 2012), but we cannot know what proportion of the urinary sample is vaginal contamination. Also, we did not explore viral or fungal data in these analyses, which focused only on the bacterial microbiome. Lastly, this population included a variety of genitourinary complaints within the sample group, making it difficult to isolate ICBPS pathology. Furthermore, the control population was mostly recruited from a subspecialty clinic where gynecological/urinary pathology is common regardless of the presence of ICBPS. ICBPS has varied presentation, and the “gold standard” used to validate some screening tools, the potassium challenge test, is known to be unreliable (Hanno et al., 2011). It is well known that ICBPS is a clinical diagnosis, but is extremely hard to diagnose with high specificity and sensitivity with current tools, given the large amoung of overlap with other pelvic pain syndromes and irritative voiding symptoms (Hanno et al., 2011; Patnaik et al., 2017). We gave preference to a questionnaire proven to correspond to disease severity and treatment response (Lubeck et al., 2001), which are most meaningful to patients and clinicians.

Conclusion

The current study highlights that ICBPS is a disease with multiple layers of dysfunction. This pilot study did not detect alterations in vaginal or urinary microbiome associated with the presence of ICBPS. However, novel findings in this study could guide future investigation of altered urinary-vaginal microbial relationships in this disease state.

Ethics Statement

University of Louisville Internal Review Board approved this study prior to any patient recruitment (IRB#13.0301), and all participating patients engaged in an informed consent conversation and signed a written consent document prior to participation or the collection of any patient data and samples.

Author Contributions

KM and VJ contributed to the protocol and project development, data collection and management, data analysis, and the writing of the manuscript. ZL contributed to the protocol and project development, data collection and management, and the editing of the manuscript. RS contributed to the data collection and management, and data analysis. JG contributed to the data analysis and the editing of the manuscript. DH contributed to the data collection and management, and the editing of the manuscript.

Funding

University of Louisville School of Medicine Basic Grant, Speed Type 57124.

Conflict of Interest Statement

KM is a textbook editor for Elsevier publishing, and in the future will receive royalties for this textbook. The content of the textbook is unrelated to the content of this research article.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2019.00092/full#supplementary-material

Supplementary Table 1. Operational taxonomic unit (OTU) level data showing the percentile of reads for each individual OUT for each sample, and the mean percentile for the sample group (ICBPS urine, unaffected urine, ICBPS vagina, unaffected vagina) for each OTU.

Supplementary Table 2. Table of alpha diversity measures (Chao1 richness measures, McIntosh evenness index, Simpson index, and Shannon index) at 5,000 sequences per sample for all specimens, with means displayed for each group (ICBPS and unaffected) for each group of samples (urinary and vaginal).

Supplementary Figure 1. OTU level data for ICBPS women for the urinary (ICBPS Image A; unaffected Image B) and vaginal spaces and unaffected women for the urinary and vagina (ICBPS Image C; unaffected Image D) spaces. Only the top 10 species for each group/space are shown. The y-axes are the percentiles of reads in the sample, and the x-axes display the individual patients and, on the far right, the mean percentile of reads in the sample for each species across the sample group. Note the large percentile of Lactobacilli species seen in the majority of samples.

Appendix 1. Details of study exclusion criteria.

Footnotes

1. ^Bacterial Vaginosis. [cited 2018 April 4] Available online at: https://www.cdc.gov/std/tg2015/bv.htm.

2. ^Castle T. JHeatChart. [cited 2018 Feb 23] Available online at: www.javaheatmap.com/.

References

Abernethy, M. G., Rosenfeld, A., White, J. R., Mueller, M. G., Lewicky-Gaupp, C., and Kenton, K. (2017). Urinary microbiome and cytokine levels in women with interstitial cystitis. Obstet. Gynecol. 129, 500–506. doi: 10.1097/AOG.0000000000001892

PubMed Abstract | CrossRef Full Text | Google Scholar

Clemens, J. Q., Elliott, M. N., Suttorp, M., and Berry, S. H. (2012). Temporal ordering of interstitial cystitis/bladder pain syndrome and non-bladder conditions. Urology 80, 1227–1231. doi: 10.1016/j.urology.2012.06.059

PubMed Abstract | CrossRef Full Text | Google Scholar

Clemens, J. Q., Meenan, R. T., O'Keeffe Rosetti, M. C., Brown, S. O., Gao, S. Y., and Calhoun, E. A. (2005). Prevalence of interstitial cystitis symptoms in a managed care population. J. Urol. 174, 576–580. doi: 10.1097/01.ju.0000165170.43617.be

PubMed Abstract | CrossRef Full Text | Google Scholar

Comeau, A. M., Douglas, G. M., and Langille, M. G. (2017). Microbiome helper: a custom and streamlined workflow for microbiome research. mSystems 2:e00127-16. doi: 10.1128/mSystems.00127-16

PubMed Abstract | CrossRef Full Text | Google Scholar

DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/AEM.03006-05

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381

PubMed Abstract | CrossRef Full Text | Google Scholar

Gottschick, C., Deng, Z. L., Vital, M., Masur, C., Abels, C., Pieper, D. H., et al. (2017). The urinary microbiota of men and women and its changes in women during bacterial vaginosis and antibiotic treatment. Microbiome 5:99. doi: 10.1186/s40168-017-0305-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanno, P. M., Burks, D. A., Clemens, J. Q., Dmochowski, R. R., Erickson, D., Fitzgerald, M. P., et al. (2011). AUA guideline for the diagnosis and treatment of interstitial cystitis/bladder pain syndrome. J. Urol. 185, 2162–2170. doi: 10.1016/j.juro.2011.03.064

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanter, G., Volpe, K. A., Dunivan, G. C., Cichowski, S. B., Jeppson, P. C., Rogers, R. G., et al. (2017). Important role of physicians in addressing psychological aspects of interstitial cystitis/bladder pain syndrome (IC/BPS): a qualitative analysis. Int. Urogynecol. J. 28, 249–256. doi: 10.1007/s00192-016-3109-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Karstens, L., Asquith, M., Davin, S., Stauffer, P., Fair, D., Gregory, W. T., et al. (2016). Does the urinary microbiome play a role in urgency urinary incontinence and its severity? Front. Cell. Infect. Microbiol. 6:78. doi: 10.3389/fcimb.2016.00078

PubMed Abstract | CrossRef Full Text | Google Scholar

Kopylova, E., Noe, L., and Touzet, H. (2012). SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics 28, 3211–3217. doi: 10.1093/bioinformatics/bts611

PubMed Abstract | CrossRef Full Text | Google Scholar

Lewis, F. M., Bernstein, K. T., and Aral, S. O. (2017). Vaginal microbiome and its relationship to behavior, sexual health, and sexually transmitted diseases. Obstet. Gynecol. 129, 643–654. doi: 10.1097/AOG.0000000000001932

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H. T., Shie, J. H., Chen, S. H., Wang, Y. S., and Kuo, H. C. (2012). Differences in mast cell infiltration, E-cadherin, and zonula occludens-1 expression between patients with overactive bladder and interstitial cystitis/bladder pain syndrome. Urology 80:225. doi: 10.1016/j.urology.2012.01.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Lozupone, C., and Knight, R. (2005). UniFrac: a new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 71, 8228–8235. doi: 10.1128/AEM.71.12.8228-8235.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Lubeck, D. P., Whitmore, K., Sant, G. R., Alvarez-Horine, S., and Lai, C. (2001). Psychometric validation of the O'leary-Sant interstitial cystitis symptom index in a clinical trial of pentosan polysulfate sodium. Urology 57(6 Suppl. 1):62–66. doi: 10.1016/S0090-4295(01)01126-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Mercier, C. B. A., and Coissac, E. (2013). SUMATRA and SUMACLUST: Fast and Exact Comparison and Clustering of Sequences. Available online at: http://metabarcoding.org/sumatra/ (accessed April 17, 2017).

Mouraviev, V., and McDonald, M. (2018). An implementation of next generation sequencing for prevention and diagnosis of urinary tract infection in urology. Can. J. Urol. 25, 9349–9356.

PubMed Abstract | Google Scholar

Navas-Molina, J. A., Peralta-Sanchez, J. M., Gonzalez, A., McMurdie, P. J., Vazquez-Baeza, Y., Xu, Z., et al. (2013). Advancing our understanding of the human microbiome using QIIME. Meth. Enzymol. 531, 371–444. doi: 10.1016/B978-0-12-407863-5.00019-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Nickel, J. C., Stephens, A., Landis, J. R., Mullins, C., van Bokhoven, A., Lucia, M. S., et al. (2016). Assessment of the lower urinary tract microbiota during symptom flare in women with urologic chronic pelvic pain syndrome: a MAPP network study. J. Urol. 195, 356–362. doi: 10.1016/j.juro.2015.09.075

PubMed Abstract | CrossRef Full Text | Google Scholar

Nunn, K. L., and Forney, L. J. (2016). Unraveling the dynamics of the human vaginal microbiome. Yale J. Biol. Med. 89, 331–337.

PubMed Abstract | Google Scholar

Patnaik, S. S., Lagana, A. S., Vitale, S. G., Buttice, S., Noventa, M., Gizzo, S., et al. (2017). Etiology, pathophysiology and biomarkers of interstitial cystitis/painful bladder syndrome. Arch. Gynecol. Obstet. 295, 1341–1359. doi: 10.1007/s00404-017-4364-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Pearce, M. M., Hilt, E. E., Rosenfeld, A. B., Zilliox, M. J., Thomas-White, K., Fok, C., et al. (2014). The female urinary microbiome: a comparison of women with and without urgency urinary incontinence. MBio 5, e01283–e01214. doi: 10.1128/mBio.01283-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Ravel, J., Gajer, P., Abdo, Z., Schneider, G. M., Koenig, S. S., McCulle, S. L., et al. (2011). Vaginal microbiome of reproductive-age women. Proc. Natl. Acad. Sci. U S A. 108 (Suppl. 1), 4680–4687. doi: 10.1073/pnas.1002611107

PubMed Abstract | CrossRef Full Text | Google Scholar

Rognes, T., Flouri, T., Nichols, B., Quince, C., and Mahe, F. (2016). VSEARCH: a versatile open source tool for metagenomics. PeerJ 4:e2584. doi: 10.7717/peerj.2584

PubMed Abstract | CrossRef Full Text | Google Scholar

Siddiqui, H., Lagesen, K., Nederbragt, A. J., Jeansson, S. L., and Jakobsen, K. S. (2012). Alterations of microbiota in urine from women with interstitial cystitis. BMC Microbiol. 12:205. doi: 10.1186/1471-2180-12-205

PubMed Abstract | CrossRef Full Text | Google Scholar

Stapleton, A. E. (2016). The vaginal microbiota and urinary tract infection. Microbiol. Spectr. 4:UTI-0025-2016. doi: 10.1128/microbiolspec.UTI-0025-2016

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomas-White, K., Forster, S. C., Kumar, N., Van Kuiken, M., Putonti, C., Stares, M. D., et al. (2018). Culturing of female bladder bacteria reveals an interconnected urogenital microbiota. Nat. Commun. 9:1557. doi: 10.1038/s41467-018-03968-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomas-White, K. J., Kliethermes, S., Rickey, L., Lukacz, E. S., Richter, H. E., Moalli, P., et al. (2017). Evaluation of the urinary microbiota of women with uncomplicated stress urinary incontinence. Am. J. Obstet. Gynecol. 216:55. doi: 10.1016/j.ajog.2016.07.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolfe, A. J., Toh, E., Shibata, N., Rong, R., Kenton, K., Fitzgerald, M., et al. (2012). Evidence of uncultivated bacteria in the adult female bladder. J. Clin. Microbiol. 50, 1376–1383. doi: 10.1128/JCM.05852-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Kobert, K., Flouri, T., and Stamatakis, A. (2014). PEAR: a fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics 30, 614–620. doi: 10.1093/bioinformatics/btt593

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: interstitial cystitis, painful bladder, bladder pain syndrome, microbiome, urine, vagina, genitourinary

Citation: Meriwether KV, Lei Z, Singh R, Gaskins J, Hobson DTG and Jala V (2019) The Vaginal and Urinary Microbiomes in Premenopausal Women With Interstitial Cystitis/Bladder Pain Syndrome as Compared to Unaffected Controls: A Pilot Cross-Sectional Study. Front. Cell. Infect. Microbiol. 9:92. doi: 10.3389/fcimb.2019.00092

Received: 14 December 2018; Accepted: 15 March 2019;
Published: 08 April 2019.

Edited by:

Shai Bel, Bar-Ilan University, Israel

Reviewed by:

Alan J. Wolfe, Loyola University Chicago, United States
Gena D. Tribble, University of Texas Health Science Center at Houston, United States

Copyright © 2019 Meriwether, Lei, Singh, Gaskins, Hobson and Jala. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Kate V. Meriwether, meriwet2@gmail.com