Qi Sun1†
Leqian Li1†Hao Zhou2†
Ying Wu3Yubo Gao1Bingyi Wu4
Yifeng Qiu1Zhipeng Zhou1
Qixiang Song5*
Jie Zhao6*
Peng Wu1*- 1Department of Urology, Nanfang Hospital, Southern Medical University, Guangzhou, China
- 2Department of Hospital Infection Management, Nanfang Hospital, Southern Medical University, Guangzhou, China
- 3School of Public Health, Southern Medical University, Guangzhou, China
- 4Medical Research Center, Nanfang Hospital, Southern Medical University, Guangzhou, China
- 5Department of Urology, Renji Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
- 6School of Pharmaceutical Sciences, Southern Medical University, Guangzhou, China
Although it is known that changes in bacterial components of the urinary microbiome are associated with overactive bladder (OAB), the specific role of viruses is still insufficiently investigated. The aim of the present study is to evaluate the role of urinary viruses in woman with OAB, and analyze the potential relationship between viruses, bacteria and disease. Catheterized urine samples were collected from 55 women with OAB and 18 control individuals. OAB patients fulfilling the following criteria were considered eligible for this study: female, 18 years of age or older; presented with classic OAB symptoms defined by the International Continence Society; and OAB Symptom Score (OABSS) total score ≥ 3 points and question 3 (urgency) score ≥ 2 points. Based on results of metagenomic next-generation sequencing (mNGS), all participants were divided into virus-infected and virus-uninfected groups for analysis. The results of mNGS showed that the diversity of the OAB group was lower than that of the control group when focused on bacterial sequences, which was consistent with our previous study. According to the questionnaire filled out by the patients, OABSS and 8-item OAB questionnaire, female OAB patients who had viruses detected in their urine had more severe symptoms. In parallel, John Cunningham virus (mainly subtype 7 and subtype 2) was the most frequently detected virus in urine. Correlation analysis indicated that risk factors for virus infection in OAB patients include age, habit of holding urine and pelvic surgery history. Given our preliminary data, viral infection can aggravate OAB severity and affect the composition of bacterial. Further research is required to explain how viral infections can aggravate OAB patient symptoms and cause bacterial changes.
Introduction
Overactive bladder (OAB) is characterized as “urgency, with or without urgency urinary incontinence (UUI), usually with increased daytime frequency and nocturia” defined by the International Continence Society (ICS; Abrams et al., 2003). The pioneering works of Brading suggested that abnormal electrical couping of smooth muscle cells resulting in detrusor contraction while other researches suggesting sensation of urgency may be provoked by urothelial/suburothelial dysfunction (Drake et al., 2001; Lee et al., 2010; Chen et al., 2017). It has been proposed that metabolic syndrome, affective disorders, sex hormone deficiency, urinary microbiota, gastrointestinal dysfunction and subclinical autonomic nervous system dysfunction are potential pathophysiological cofactors of OAB (Peyronnet et al., 2019). However, there are still a number of patients failed to respond to the medications, including antimuscarinics and the β3-adrenoreceptor agonist, requiring invasive interventions (Wadie, 2014; De Wachter et al., 2019). Therefore, the multifactorial nature of OAB and the effective pharmacotherapy of refractory cases compose the two main obstacles in basic science research and clinical management.
Emerging evidence indicated that alteration to urinary microbiome is closely related to the pathogenesis of OAB using 16S rRNA gene sequencing and expanded quantitative urine culture (EQUC) (Price et al., 2016; Gasiorek et al., 2019; Lundy et al., 2021). Wolfe et al. identified urinary bacteria community with 16S rRNA sequencing in female patients undergoing urogynecologic surgery with negative standard urine cultures, and developed its linkage with diverse urinary diseases (Wolfe et al., 2012). Pearce et al. (2014) used EQUC and 16S rRNA sequencing to compare the urinary microbiome of women with and without wet OAB (UUI) and their data suggest a potentially important difference in the urinary microbiome of women in the two groups. Combined, with the in-depth research on the urinary microbiome, from sterile to fertile, a growing number of scholars realize that the current research results are just the tip of the iceberg (Hilt et al., 2014; Price et al., 2016; Brubaker et al., 2021).
The urinary microbiome, however, contains not just bacteria, but also viruses, fungi, parasites, and other species (Lloyd-Price et al., 2016). Incidentally, excluding the bacteria, studies on other microorganisms, particularly viruses and their association with OAB are limited (Brubaker et al., 2021). Among them, most frequently detected viruses in urinary tract are human polyomavirus type 1 (BK Virus) and 2 (JC Virus), both of which are mainly studied in immunocompromised patients, such as kidney and bone marrow transplant recipients (Jhang et al., 2018). Meanwhile, both JC and BK are double-stranded DNA viruses and belong to polyomavirus (Salabura et al., 2021). Detection of JC Virus is reported in recent studies of urinary diseases such as lower urinary tract symptoms (LUTS) and OAB, but its relationships with urinary bacteria and OAB symptoms require further research (Garretto et al., 2018; Thomas et al., 2021). Additionally, human papillomavirus, adenovirus, and cytomegalovirus infections have been associated with LUTS (Paduch, 2007; Brubaker et al., 2021). Although urinary tract infections (UTIs) are mainly caused by bacteria such as Escherichia coli and viral infections are rare pathogens, viral infections of the urinary tract are linked to considerable morbidity and distress, including increased mortality in immunocompromised patients (Khonsari et al., 2021; Salabura et al., 2021). However, there is an important link between bacteria and viruses, such as increased antibiotic resistance in UTI patients, which leads to an increase in the global medical burden, and bacteriophages, as the most abundant group of human viruses, can be used as a feasible treatment for UTI (Behzadi et al., 2020; Ahmadi, 2021). It is important to note that, especially in female patients, misdiagnosis of UTI and genital tract infections (GTIs) may lead to imbalances in bacterial and viral populations and worsen disease (Behzadi et al., 2019).
Next-generation sequencing (NGS) of bacterial 16S rRNA is widely applied in clinical pathogen detection including urology (Derosa et al., 2020; Salabura et al., 2021). One such application is the unbiased next-generation metagenomic sequencing (mNGS), which overcomes the limitations of current diagnostic tests and allows untargeted, culture-independent pathogen detection directly from clinical samples (Neugent et al., 2020). At the same time, the development of bioinformatics for the analysis of sequencing data can help reduce possible biases in the results to obtain accurate and reliable analysis results (Behzadi and Ranjbar, 2018). Here, we seek to evaluate the role of urinary viruses in the onset of OAB using a mNGS approach, and attempt to reveal the potential relationship between viruses, bacteria, and OAB patients.
Materials and methods
Clinical procedures and sample acquisition
A total of 73 eligible individuals were identified from the outpatient clinic, including 55 OAB patients and 18 asymptomatic controls between January 2021 and August 2021. OAB patients fulfilling the following criteria were considered eligible for this study: (1) female, 18 years of age or older; (2) presented with classic OAB symptoms defined by the ICS; and (3) OAB Symptom Score (OABSS) total score ≥ 3 points and question 3 (urgency) score ≥ 2 points. The exclusion criteria include current urinary tract infection, with indwelling catheter, with antibiotic exposure in the past 30 days, chemical cystitis, cyclophosphamide cystitis, genitourinary cancer, urinary stones, neurogenic bladder, pelvic organ prolapse, pelvic radiation and pregnancy. The inclusion/exclusion criteria were consistent for the control subjects, except for the absence of OAB symptoms. The degree of OAB was appraised with both OABSS (mild ≤ 5, 6 < moderate < 11, severe ≥ 12) and 8-item OAB questionnaire (OAB-V8) during the recruitment. OAB patients were separated into the OAB virus-infected (OAB-VI) group and OAB virus-uninfected (OAB-VU) group based on the mNGS results. Similarly, the control group was divided into the control virus-infected group (CON-VI) and control virus-uninfected group (CON-VU).
As previously described, urine specimens were collected by transurethral catheterization to avoid contamination from distal urethra (Wolfe et al., 2012). In brief, with the patient in a lithotomy position, the perineum is thoroughly disinfected, followed by introducing a 16Fr catheter into the bladder to collect urine (Karstens et al., 2016; Yoo et al., 2021). Finally, 50 ml of urine was collected from each participant in this study, of which 20 ml was tested using standard urine culture methods to exclude UTI. The remaining 30 ml urine sample was temporarily stored in a 4°C refrigerator for 1 h and transferred to a −80°C freezer with DNA/RNA Shield (Zymo Research Corporation, Irvine, CA, United States) added until further processing. The study protocol was approved by the Institutional Review Board of Nanfang hospital (reference number: NFEC-2020-123).
DNA extraction and metagenomic next generation sequencing
The DNA extraction was extracted using the TIANamp Micro DNA Kit (DP316, TIANGEN BIOTECH) sonicated to a size of 200–300 bps fragments (Bioruptor Pico protocols). Then, DNA libraries were constructed through DNA-fragmentation, end-repair, adapter-ligation and PCR amplification.
Quality qualified libraries were pooled. DNA Nanoball was made and sequenced by BGISEQ-50 platform. High-quality sequencing data was generated by removing low-quality reads, followed by computational subtraction of human host sequences mapped to the human reference genome (hg19) using Burrows–Wheeler Alignment tool1 (Li and Durbin, 2009). Low-quality reads refer to the basic characteristics that do not meet the following indicators: (1) the number of Q30 bases is more than 80%; (2) joint contamination ratio should not exceed 1%; (3) the effective sequence length is not <50 bp; and (4) the effective alignment rate of the data should be greater than 70% (Gu et al., 2019). The remaining data by removal of low-complexity reads were classified by simultaneously aligning to pathogens metagenomics Database (PMDB, BGI-Shenzhen), consisting of bacteria, fungi, viruses and parasites. We used a threshold of 10% minimum coverage of the viral genome or a viral stringently mapped reads number >100 for diagnosis (Jerome et al., 2019). The raw data of mNGS have been submitted to the Short Read Archive (PRJNA814640).
Clinical data and bioinformatics analysis
Clinical characteristics between or among OAB and control groups were compared by the t-test or Mann–Whitney U-test for continuous variables, and by the Chi-square test for categorical variables. The Spearman correlation matrix was calculated to examine the pairwise between any two of all variables. Binary logistic regression was preformed and a predictive nomogram was drawn. All p values were two tailed, and the statistical significance level was set as 0.05. All statistical analyses were performed by R software (Version 4.1.1) and SPSS software (Version 21).
Kraken2 was used to annotate and classify all valid sequences of each sample to analyze microbial composition and diversity information (Wood and Salzberg, 2014). After using Kraken2 to classify mNGS data, Bracken 2.0 (default parameters) was used to estimate the relative abundance of metagenomic samples at different classification levels. Chao1 index, Shannon index and Simpson index calculated by QIIME2 were used to evaluate Alpha diversity. Through linear discriminant analysis (LDA) effect size (LEfSe) analysis, the differences of bacteria characteristics at all classification levels were identified (Segata et al., 2011).
Phylogenetic tree analysis
We searched all available JC Virus complete genomes from the NCBI database (http://www.ncbi.nlm.nih.gov/, last accessed September 14, 2021) and combined with previous research databases that identified JC Virus subtypes and population distribution (Shackelton et al., 2006; Garretto et al., 2018; Forni et al., 2020). Next, 672 available JC Virus complete genomes and 24 JC Virus genomes identified through mNGS were compared, using MUSCLE to align the sequences through Geneious (Version 10.1.2). Subsequently, the maximum likelihood multi-model inference method was used to assign each query sequence to a reference sequence with known genotypes, followed by using the IQ-TREE to obtain a phylogenetic tree and the online tool Interactive Tree Of Life2 for annotation (Nguyen et al., 2015).
Results
The characteristics and microbiota composition of subjects: OAB vs. control
In total, 55 OAB patients and 18 asymptomatic controls were recruited in the study. There was no significant difference in the demographic characteristics between the two groups, with the exception of UTI history (60% vs. 17%, p = 0.002), pelvic surgery history (43% vs. 0%, p = 0.000) and the habit of holding urine (36% vs. 6%, p = 0.015, Table 1).
Table 1. The comparisons of demographic characteristics and symptom scores between OAB patients and Asymptomatic controls.
A total of 5,883,488 bacterial reads were obtained. At the phylum level, Proteobacteria, Firmicutes, Actinobacteria and Bacteroidetes were the most common types in both groups, with Bacteroidetes decreased significantly in patients with OAB (p = 0.0005, Figures 1A,B). Although the Simpson index (p = 0.056, Figure 1C) was similar, the Shannon index and chao1 index were significantly reduced in the OAB group at the species level (p = 0.043, Figure 1D; p = 0.022, Figure 1E). The clinical information corresponding to the urine sample and microbial composition and relative abundance are shown in Figure 1F, which is intuitive to show that Lactobacillus inner, Burkholderia vietnamiensis and Escherichia coli are enriched in the OAB group. To identify specific OAB related bacteria, the LEfSe algorithm was used, indicating significantly elevated compositional abundance of 14 species in OAB patients and other 85 species in control subjects (Figure 2).
Figure 1. The urinary microbiome profile of participants. (A) The bar graph showing the relative abundance of urine microbiota at phylum level. (B) Composition Changes of Actinobacteria, Bacteroidetes, Firmicutes and Proteobacteria between the control and OAB groups. Alpha diversity analysis for the control and OAB urinary microbiomes including Simpson index (C), Shannon index (D) and Chao1 index (E). (F) Heatmap demonstrates the relative abundance of the microbiota at different levels of the control and OAB groups and clinical data.
Figure 2. Microbial taxa associated with Overactive bladder. (A) Cladogram representation of the urinary microbial taxa associated with OAB (yellow red) and control (blue). (B) Association of specific microbiota taxa with OAB group and control group by linear discriminant analysis effect size (LEfSe). Yellow red indicates taxa enriched in OAB group and blue indicates taxa enriched in control group. Only species meeting a linear discriminant analysis score threshold >2 are shown.
The characteristics and microbiota composition of OAB patients with or without urinary virus
Twenty-six OAB patients (47.3%) identified with positive urinary virus were classified as OAB-VI group and the remaining 29 patients (52.7%) without detectable virus were allocated to the OAB-VU group (Supplementary Table 1). Notably, patients in the OAB-VI group demonstrated much worse OAB symptoms than that in the OAB-VU group as evidenced by both the OABSS and OAB-V8 questionnaires (Table 2). Plus, there were significantly more moderate-to-severe cases in the OAB-VI group, while mild OAB symptoms were predominant if without virus detection (Supplementary Figure S1). Crucially, at species level, the relative abundance of Staphylococcus warneri (p = 0.0092), Staphylococcus hominis (p = 0.0017) and Staphylococcus epidermidis (p = 0.0041) was significantly increased in the OAB-VI group (Figures 3A–C; Supplementary Figure S2). In contrast, there was no significant difference in Lactobacillus inner richness in OAB-VI compared with OAB-VU, but there was a downward trend (Figure 3D).
Table 2. The comparison of clinical features of OAB patients with or without viral infection.
Figure 3. Bacterial average relative abundance in OAB-VI and OAB-VU. Relative abundance of (A) Staphylococcus warneri, (B) Staphylococcus hominis, (C) Staphylococcus epidermidis and (D) Lactobacillus inner. *p≤0.05, **p≤0.01.
Virus genomes were detected in 32 of 73 samples, of which 26 samples were from the OAB group and 6 samples were from the control group (Supplementary Table 1). JC virus was the most common of all viral genomes (n = 25), BK virus, EB virus and human beta-herpes virus 6A were also detected (Supplementary Table 1). However, the JC virus genome detected in urine sample of OAB-VI25 was not included because the coverage rate was too low to be classified (Supplementary Table 1). The main subtype of JC Virus strains isolated in the samples were mainly subtype 7B (OAB, n = 8; control, n = 1) and subtype 7A (OAB, n = 7; control, n = 2) followed by subtype 7C (OAB, n = 2; control, n = 0), subtype 2D (OAB, n = 2; control, n = 0) and subtype 2A (OAB, n = 1; control, n = 1). Based on the geographic distribution of these virus subtypes, we found that they were mainly from the Asian region (n = 19, Figure 4).
Figure 4. Phylogenetic analysis JC virus. Phylogenetic tree of 672 JC Virus complete genomes inferred using a maximum likelihood approach, including strains in this study (a red square at the end of a branch indicates). The genomes from OAB-VI (red squares) and control (black squares) subjects were inserted into the circular dendrite map. The tree branch color represents different types and peripheral color bands represent geographic distribution. Samples analyzed in the study are labeled with the sample name and JC Virus subtype.
Finally, by comparing alpha diversity of the CON-VI, CON-VU, OAB-VI and OAB-VU groups, we only found that the Chao1 index of the OAB-VI group was significantly lower than that of the CON-VU group (Supplementary Figures S3A–C).
Correlation analysis between OAB symptoms and risk factors for viral infection in OAB patients
Spearman correlation analysis was used to identify intra-variable linear correlation between variables, results indicated existence of high correlations between viral infection and OABSS (r = 0.58), age and pausimenia (r = 0.68), hypertension and age (r = 0.53) respectively (Figure 5A). In addition, viral infection is associated with over 10 times more likely to have a worse OAB symptom (Table 3).
Figure 5. Correlation analyze with OAB symptoms and risk factors for viral infection. (A) Heatmap for pearson correlation coefficient (PCCs) matrix for virus infection, age, BMI, pregnancy, HP, DM, pausimenia, estrogen use, UTI history, habit of holding urine history, OABSS, age group, BMI group, Pelvic surgery, OAB severity (n = 55). (B) Forest plot for screened risk factor and protective factor of viral infection (p < 0.05). (C) Receiver operating characteristic curve (ROC curve): area under the ROC curve is 0.915 (p < 0.001). (D) Nomogram for the prediction of viral infection. Instruction for utilization: locate the patient’s age on the corresponding axis. Draw a line straight down to the point axis to determine how many points toward the probability of viral infection. Repeat the process for each additional predictor and add up all the points. Locate the final sum on the total points axis. Draw a line straight up to find the patient’s risk of viral infection so as to provide guidance of mNGS application. (E) Calibration curve for the nomogram to predict probability of viral infection.
Table 3. OR value calculation of risk factors for OAB symptom severity.
Estimates from binary logistic regression model indicated risk factors for virus infection in OAB patients including age (OR = 1.99; 95%CI: 0.02, 2.61), holding urine habit (OR = 2.16; 95%CI: 0.18, 3.85) and pelvic surgery (OR = 2.53; 95%CI: 0.54, 4.27). Protective factors include estrogen use history (OR = −0.74; 95% CI: −4.37, 1.97) (Table 4; Figure 5B). Five clinical variables with values greater than 1 are: virus infection (VIP = 1.94), UTI history (VIP = 1.21), BMI (VIP = 1.13), pausimenia (VIP = 1.02, Supplementary Figure S4). The area under ROC curve is 0.915 (p < 0.0001) for our binary logistic regression predictive model (Figures 5C,D). Then, nomogram predictive model based on visualized binary logistic regression result was built for virus infection prediction, including patient history of hypertension, estrogen usage, UTI history, pelvic surgery history, the habit of holding urine as independent predictive factors for virus infection (Figure 5E).
Table 4. Risk factors for virus infection identified by binary logistic regression.
Discussion
The correlations between alteration of urine microbiome and occurrence of OAB have not been revealed until recent years (Wolfe et al., 2012; Pearce et al., 2014; Wu et al., 2017; Brubaker et al., 2021). We have previously demonstrated that at the genus level, the abundance of Aerococcus and Staphylococcus increased, while Prevotella and Lactobacillus decreased in OAB patients (Wu et al., 2017). By comparing the urinary microbiome of women with UUI and asymptomatic controls, Pearce et al. found that the microbiome from UUI showed a decreased in Lactobacillus and an increased in Grdnerella abundance from the UUI cohort by EQUC (Pearce et al., 2014). Additionally, Karstens et al. (2016) determined that an increased severity of UUI symptoms was associated with decreased microbial diversity in female UUI patients. However, most studies on urinary microbiome were measured by 16S rRNA gene sequencing, and investigations were limited to the bacteria. In the present study, we took the lead in unveiling the effects of urinary viral infection on the presence of OAB, and providing evidence on the connections between modified microbial features and degree of OAB symptoms.
Consistent with our previous report, the analysis of Chao1 and Shannon index indicated that the OAB group had lower diversity than the control group, which was consistent from both 16S rRNA gene sequencing and mNGS results (Wu et al., 2017). Since most studies on the urinary microbiome used 16S rRNA amplicon sequencing and focused on the genus level analysis, few studies have similar results at the species level. Recently, Wolfe et al. have provided new insights into taxonomic schemes for species-level bacterial annotations in 16S amplicon data (Hoffman et al., 2021). In the results obtained by mNGS, at species level, the OAB-VI group demonstrated significantly decreased Lactobacillus iners, but increased coagulase-negative Staphylococci (CoNS), including Staphylococcus epidermidis, Staphylococcus hominis and Staphylococcus warner. Previous studies of Staphylococcus warneri, Staphylococcus epidermidis and Staphylococcus hominis are mainly elaborated in human skin microbiome (Byrd et al., 2018). Staphylococcus epidermidis may act both as an opportunistic pathogen and a probiotics, which prevents bacterial biofilm formation via expression of a serine type protease, termed Esp (Otto, 2009; Sugimoto et al., 2013). On the contrary, known as the most sensitive microbiome marker in urinary tract diseases, a decreased Lactobacillus in the urine was associated with the presence of OAB/UUI (Wu et al., 2017; Peyronnet et al., 2019). In a study, whole-genome phylogenetic analysis of bacterial strains isolated from the vagina and bladder of the same female identified highly similar Lactobacillus iners and Lactobacillus crispatus, both of these Lactobacillus have been linked to health in the vagina, as well as the bladder in asymptomatic women (Thomas-White, 2018). These conclusions may indicate its potential role in prevention and treatment of OAB as a probiotics.
Viral genomes, especially JC Virus and BK Virus, were detected predominantly more frequent in patients with OAB. Similarly, Garretto et al. (2018) found that urine samples from OAB patients contain complete genomes of JC Virus while no viral genes were detected in the control group. In reality, JC Virus may infect healthy persons latently, resulting in asymptomatic replication and sporadic detection of the virus in the urine in up to 30% of the general population (White and Khalili, 2011). Previous studies have determined the subtype of JC Virus and global distribution (Forni et al., 2020). The main subtype of JC Virus strains detected in urine of OAB patients were subtype 7B (n = 8) and subtype 7A (n = 7) and these two subtypes of JC Viruses are mainly prevalent in Asia, which was obviously inconsistent with the subtype measured by Wolfe et al., since their samples were mainly from Caucasians and African Americans (Garretto et al., 2018). Interactions between virus and bacteria are widely reported in gastrointestinal tract, mostly focusing on the detrimental effects of bacteriophage on bacteria and end up with reduce in specific bacteria (Santiago-Rodriguez and Hollister, 2019). The bacterial microbiota is required for many enteric viruses during the infection of eukaryotic cells (Monedero et al., 2018). The effects of bacteria on eucaryotic viruses include direct binding with virions to enhance stability and indirect effects led by immune activation through the attachment of host cells (Rowe et al., 2019). Nevertheless, the evidence on whether the virus interacts directly or indirectly with urinary tract bacteria remains to be discovered.
In a cross-sectional study, women with unhealthy habit of holding urine suffer from more frequent LUTS, and this feature was also found in our clinical practice (Reynolds et al., 2020). Here, we found that history of pelvic surgery, holding urine habit and UTI history are risk factors for viral infection, which may make the symptoms of OAB worse. However, GTIs was not taken into account in this study. The association between history of GTIs and OAB and viral infection in women deserves further exploration (Behzadi et al., 2019, 2020). Estrogen usage history whose protective impact of vaginal estrogen therapy has been reported on postmenopausal OAB patients associated with increase of Lactobacillus is found to be a protective factor of viral infection (Thomas-White et al., 2020). Patients in the virus infected group present more severe symptoms and received higher OABSS and OAB-V8 score, which could be explained by the bladder epithelium damage caused by JC or BK viral replication (Li et al., 2013). Given the multifactorial nature of OAB, further research and clinical observations to identify the triple relationship among OAB, virus, bacteria as well as the importance of viral detection and targeted treatment in OAB patients would greatly benefit clinical practice.
Conclusion
In conclusion, we focused at the bacterial and viral composition of urine from OAB and healthy individuals, and found that more severe symptoms of OAB patients seems to be related to viral infection. This is the first published research on the function of human urine bacteria and viruses in OAB female patients that we are aware of. It is our hope that this work lays the foundation for future research on female urine microbiota and OAB disease.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
PW, JZ, and QSo: study concept and design. PW, HZ, YG, BW, YQ, and ZZ: investigation. PW, HZ, and BW: supervision. HZ, QSu, LL, and YW: analysis and interpretation of data. PW, QSo, QSu, and LL: writing – original draft. PW, JZ, QSo, and YW: writing – review and editing. PW: funding acquisition. All authors contributed to the article and approvedthe submitted version.
Funding
This current study was supported by funding from the National Natural Science Foundation of China (grant nos. 81870522 and 82173304) and the Natural Science Foundation of Guangdong Province (grant nos. 2018A030313148 and 2021A1515012262).
Acknowledgments
We would like to thank all the participants in the study from Nanfang Hospital for their assistance to this study.
Conflict of interest
The 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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.984234/full#supplementary-material
Footnotes
References
Abrams, P., Cardozo, L., Fall, M., Griffiths, D., Rosier, P., Ulmsten, U., et al. (2003). The standardisation of terminology in lower urinary tract function: report from the standardisation sub-committee of the International Continence Society. Urology 61, 37–49. doi: 10.1016/S0090-4295(02)02243-4
Ahmadi, M. (2021). Virulence factors, antibiotic resistance patterns, and molecular types of clinical isolates of Klebsiella pneumoniae. Expert Rev. Anti-Infect. Ther. 20, 463–472. doi: 10.1080/14787210.2022.1990040
Behzadi, P., Behzadi, E., Agnieszka, E., and Adamska, P. (2019). Urinary tract infections (UTIs) or genital tract infections (GTIs)? It’s the diagnostics that count. GMS Hyg. Infect. Control 14, Doc14. doi: 10.3205/dgkh000320
Behzadi, P., and Ranjbar, R. (2018). DNA microarray technology and bioinformatic web services. Acta Microbiol. Immunol. Hung. 66, 19–30. doi: 10.1556/030.65.2018.028
Behzadi, P., Urbán, E., Matuz, M., Benkő, R., and Gajdács, M. (2020). “The role of Gram-negative bacteria in urinary tract infections: current concepts and therapeutic options,” in Advances in microbiology, infectious diseases and public health advances in experimental medicine and biology. ed. G. Donelli (Cham: Springer International Publishing), 35–69.
Brubaker, L., Putonti, C., Dong, Q., and Wolfe, A. J. (2021). The human urobiome. Mamm. Genome 32, 232–238. doi: 10.1007/s00335-021-09862-8
Byrd, A. L., Belkaid, Y., and Segre, J. A. (2018). The human skin microbiome. Nat. Rev. Microbiol. 16, 143–155. doi: 10.1038/nrmicro.2017.157
Chen, S.-L., Ng, S.-C., Huang, Y.-H., and Chen, G.-D. (2017). Are patients with bladder oversensitivity different from those with urodynamically proven detrusor overactivity in female overactive bladder syndrome? J. Chin. Med. Assoc. 80, 644–650. doi: 10.1016/j.jcma.2017.03.009
De Wachter, S., Benson, K. D., Dmochowski, R. R., Rovner, E. S., Versi, E., Miller, L. E., et al. (2019). Six-month results of selective bladder denervation in women with refractory overactive bladder. J. Urol. 201, 573–580. doi: 10.1016/j.juro.2018.09.043
Derosa, L., Routy, B., Fidelle, M., Iebba, V., Alla, L., Pasolli, E., et al. (2020). Gut bacteria composition drives primary resistance to cancer immunotherapy in renal cell carcinoma patients. Eur. Urol. 78, 195–206. doi: 10.1016/j.eururo.2020.04.044
Drake, M., Mills, I., and Gillespie, J. (2001). Model of peripheral autonomous modules and a myovesical plexus in normal and overactive bladder function. Lancet 358, 401–403. doi: 10.1016/S0140-6736(01)05549-0
Forni, D., Cagliani, R., Clerici, M., Pozzoli, U., and Sironi, M. (2020). You will never walk alone: codispersal of JC polyomavirus with human populations. Mol. Biol. Evol. 37, 442–454. doi: 10.1093/molbev/msz227
Garretto, A., Thomas-White, K., Wolfe, A. J., and Putonti, C. (2018). Detecting viral genomes in the female urinary microbiome. J. Gen. Virol. 99, 1141–1146. doi: 10.1099/jgv.0.001097
Gasiorek, M., Hsieh, M. H., and Forster, C. S. (2019). Utility of DNA next-generation sequencing and expanded quantitative urine culture in diagnosis and management of chronic or persistent lower urinary tract symptoms. J. Clin. Microbiol. 58, e00204-19. doi: 10.1128/JCM.00204-19
Gu, W., Miller, S., and Chiu, C. Y. (2019). Clinical metagenomic next-generation sequencing for pathogen detection. Annu. Rev. Pathol. Mech. Dis. 14, 319–338. doi: 10.1146/annurev-pathmechdis-012418-012751
Hilt, E. E., McKinley, K., Pearce, M. M., Rosenfeld, A. B., Zilliox, M. J., Mueller, E. R., et al. (2014). Urine is not sterile: use of enhanced urine culture techniques to detect resident bacterial Flora in the adult female bladder. J. Clin. Microbiol. 52, 871–876. doi: 10.1128/JCM.02876-13
Hoffman, C., Siddiqui, N. Y., Fields, I., Gregory, W. T., Simon, H. M., Mooney, M. A., et al. (2021). Species-level resolution of female bladder microbiota from 16S rRNA amplicon sequencing. mSystems 6, e0051821. doi: 10.1128/mSystems.00518-21
Jerome, H., Taylor, C., Sreenu, V. B., Klymenko, T., Filipe, A. D. S., Jackson, C., et al. (2019). Metagenomic next-generation sequencing aids the diagnosis of viral infections in febrile returning travellers. J. Infect. 79, 383–388. doi: 10.1016/j.jinf.2019.08.003
Jhang, J.-F., Hsu, Y.-H., Peng, C.-W., Jiang, Y.-H., Ho, H.-C., and Kuo, H.-C. (2018). Epstein-Barr virus as a potential etiology of persistent bladder inflammation in human interstitial cystitis/bladder pain syndrome. J. Urol. 200, 590–596. doi: 10.1016/j.juro.2018.03.133
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
Khonsari, M. S., Behzadi, P., and Foroohi, F. (2021). The prevalence of type 3 fimbriae in Uropathogenic Escherichia coli isolated from clinical urine samples. Meta Gene 28:100881. doi: 10.1016/j.mgene.2021.100881
Lee, S. R., Kim, H. J., Kim, A., and Kim, J. H. (2010). Overactive bladder is not only overactive but also hypersensitive. Urology 75, 1053–1059. doi: 10.1016/j.urology.2009.10.045
Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754–1760. doi: 10.1093/bioinformatics/btp324
Li, R., Sharma, B. N., Linder, S., Gutteberg, T. J., Hirsch, H. H., and Rinaldo, C. H. (2013). Characteristics of polyomavirus BK (BKPyV) infection in primary human urothelial cells. Virology 440, 41–50. doi: 10.1016/j.virol.2013.01.024
Lloyd-Price, J., Abu-Ali, G., and Huttenhower, C. (2016). The healthy human microbiome. Genome Med. 8, 51. doi: 10.1186/s13073-016-0307-y
Lundy, S. D., Sangwan, N., Parekh, N. V., Selvam, M. K. P., Gupta, S., McCaffrey, P., et al. (2021). Functional and taxonomic dysbiosis of the gut, urine, and semen microbiomes in male infertility. Eur. Urol. 79, 826–836. doi: 10.1016/j.eururo.2021.01.014
Monedero, V., Collado, M. C., and Rodríguez-Díaz, J. (2018). Therapeutic opportunities in intestinal microbiota–virus interactions. Trends Biotechnol. 36, 645–648. doi: 10.1016/j.tibtech.2017.12.009
Neugent, M. L., Hulyalkar, N. V., Nguyen, V. H., Zimmern, P. E., and De Nisco, N. J. (2020). Advances in understanding the human urinary microbiome and its potential role in urinary tract infection. mBio 11, e00218-20. doi: 10.1128/mBio.00218-20
Nguyen, L.-T., Schmidt, H. A., von Haeseler, A., and Minh, B. Q. (2015). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274. doi: 10.1093/molbev/msu300
Otto, M. (2009). Staphylococcus epidermidis – the “accidental” pathogen. Nat. Rev. Microbiol. 7, 555–567. doi: 10.1038/nrmicro2182
Paduch, D. A. (2007). Viral lower urinary tract infections. Curr. Urol. Rep. 8, 324–335. doi: 10.1007/s11934-007-0080-y
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-14. doi: 10.1128/mBio.01283-14
Peyronnet, B., Mironska, E., Chapple, C., Cardozo, L., Oelke, M., Dmochowski, R., et al. (2019). A comprehensive review of overactive bladder pathophysiology: on the way to tailored treatment. Eur. Urol. 75, 988–1000. doi: 10.1016/j.eururo.2019.02.038
Price, T. K., Dune, T., Hilt, E. E., Thomas-White, K. J., Kliethermes, S., Brincat, C., et al. (2016). The clinical urine culture: enhanced techniques improve detection of clinically relevant microorganisms. J. Clin. Microbiol. 54, 1216–1222. doi: 10.1128/JCM.00044-16
Reynolds, W. S., Kowalik, C., Kaufman, M. R., Dmochowski, R. R., and Fowke, J. H. (2020). Women’s perceptions of public restrooms and the relationships with toileting behaviors and bladder symptoms: a cross-sectional study. J. Urol. 204, 310–315. doi: 10.1097/JU.0000000000000812
Rowe, H. M., Meliopoulos, V. A., Iverson, A., Bomme, P., Schultz-Cherry, S., and Rosch, J. W. (2019). Direct interactions with influenza promote bacterial adherence during respiratory infections. Nat. Microbiol. 4, 1328–1336. doi: 10.1038/s41564-019-0447-0
Salabura, A., Łuniewski, A., Kucharska, M., Myszak, D., Dołęgowska, B., Ciechanowski, K., et al. (2021). Urinary tract virome as an urgent target for metagenomics. Life (Basel) 11, 1264. doi: 10.3390/life11111264
Santiago-Rodriguez, T. M., and Hollister, E. B. (2019). Human virome and disease: high-throughput sequencing for virus discovery, identification of phage-bacteria dysbiosis and development of therapeutic approaches with emphasis on the human gut. Viruses 11, 656. doi: 10.3390/v11070656
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12, R60. doi: 10.1186/gb-2011-12-6-r60
Shackelton, L. A., Rambaut, A., Pybus, O. G., and Holmes, E. C. (2006). JC virus evolution and its association with human populations. J. Virol. 80, 9928–9933. doi: 10.1128/JVI.00441-06
Sugimoto, S., Iwamoto, T., Takada, K., Okuda, K., Tajima, A., Iwase, T., et al. (2013). Staphylococcus epidermidis Esp degrades specific proteins associated with Staphylococcus aureus biofilm formation and host-pathogen interaction. J. Bacteriol. 195, 1645–1655. doi: 10.1128/JB.01672-12
Thomas, S., Dunn, C. D., Campbell, L. J., Strand, D. W., Vezina, C. M., Bjorling, D. E., et al. (2021). A multi-omic investigation of male lower urinary tract symptoms: potential role for JC virus. PLoS One 16:e0246266. doi: 10.1371/journal.pone.0246266
Thomas-White, K. (2018). Culturing of female bladder bacteria reveals an interconnected urogenital microbiota. Nat. Commun. 9, 1557. doi: 10.1038/s41467-018-03968-5
Thomas-White, K., Taege, S., Limeira, R., Brincat, C., Joyce, C., Hilt, E. E., et al. (2020). Vaginal estrogen therapy is associated with increased lactobacillus in the urine of postmenopausal women with overactive bladder symptoms. Am. J. Obstet. Gynecol. 223, 727.e1–727.e11. doi: 10.1016/j.ajog.2020.08.006
Wadie, B. S. (2014). Management of refractory OAB in the non-neurogenic patient. Curr. Urol. Rep. 15, 438. doi: 10.1007/s11934-014-0438-x
White, M. K., and Khalili, K. (2011). Pathogenesis of progressive multifocal Leukoencephalopathy—revisited. J. Infect. Dis. 203, 578–586. doi: 10.1093/infdis/jiq097
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
Wood, D. E., and Salzberg, S. L. (2014). Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 15, R46. doi: 10.1186/gb-2014-15-3-r46
Wu, P., Chen, Y., Zhao, J., Zhang, G., Chen, J., Wang, J., et al. (2017). Urinary microbiome and psychological factors in women with overactive bladder. Front. Cell. Infect. Microbiol. 7:488. doi: 10.3389/fcimb.2017.00488
Keywords: female urinary microbiome, overactive bladder, metagenome, virome, microbiota
Citation: Sun Q, Li L, Zhou H, Wu Y, Gao Y, Wu B, Qiu Y, Zhou Z, Song Q, Zhao J and Wu P (2022) The detection of urinary viruses is associated with aggravated symptoms and altered bacteriome in female with overactive bladder. Front. Microbiol. 13:984234. doi: 10.3389/fmicb.2022.984234
Edited by:
Seil Kim, Korea Research Institute of Standards and Science, South KoreaReviewed by:
Payam Behzadi, Islamic Azad University, IranMangesh Vasant Suryavanshi, Lerner Research Institute, Cleveland Clinic, United States
Copyright © 2022 Sun, Li, Zhou, Wu, Gao, Wu, Qiu, Zhou, Song, Zhao and Wu. 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: Peng Wu, ZG9jdG9yd3VwZW5nQGdtYWlsLmNvbQ==; Jie Zhao, emhhb2ppZV8wNDEyQDE2My5jb20=; Qixiang Song, c29uZ3FpeGlhbmdAcmVuamkuY29t
†These authors have contributed equally to this work