Profiling the Urinary Microbiota in Male Patients With Bladder Cancer in China

Introduction

Bladder cancer is diagnosed in more than 430,000 patients worldwide every year, making it the ninth most common malignancy (Kamat et al., 2016). In addition, most of bladder cancer patients are male, five to six times more common than female in China (Li et al., 2015). In the past decades, bladder cancer has aroused scientists' attention for its high morbidity and mortality rates. Unfortunately, the etiology and pathophysiology of bladder cancer remain unknown. It may be caused by genetic mutations and external risk factors, including tobacco smoking, carcinogen exposure, the chlorination of drinking water and possibly cyclophosphamide (Babjuk et al., 2017).

Notably, all above factors are known to affect the composition of microbiota, the ensemble of symbiotic bacteria, fungi, parasites, and viruses that inhabit the epithelial barrier surfaces of our body (Costello et al., 2012). The microbiome affects human physiological functions, such as metabolism, immunity and haematopoiesis (Dzutsev et al., 2015). In addition, the microbiome also plays a role in the development of malignancies both at epithelial barriers and in tissues (Roy and Trinchieri, 2017). Studies suggest that microbial dysbiosis at various body sites may promote disease progression, such as periodontitis, inflammatory bowel disease, colorectal cancer and breast cancer (Darveau, 2010; Irrazábal et al., 2014; Chu et al., 2016; Urbaniak et al., 2016). Recently, emerging evidence overturns the dogma that urine in healthy individuals must be sterile (Whiteside et al., 2015). Furthermore, dysbiosis of the urinary microbiome has been revealed responsible for various urological disorders, such as urgency urinary incontinence, interstitial cystitis, overactive bladder and prostate cancer (Siddiqui et al., 2012; Pearce et al., 2014; Wu et al., 2017; Shrestha et al., 2018).

Accordingly, it is conceivable that alteration of urinary microbiome may be associated with bladder cancer. Although the link between specific pathogens and cancer is well established, such as Helicobacter pylori and gastric cancer (El-Omar et al., 1997), there is currently no hard evidence linking microbiome and bladder cancer. Intravesical instillation of Mycobacterium bovis bacillus Calmette–Guérin and oral administration of Lactobacillus after removal of the bladder tumor could reduce the probability of recurrence (Aso et al., 1995; Zitvogel et al., 2017). Dysbiosis caused by repeated antibiotic use can increase the incidence of cancers, including bladder cancer (Boursi et al., 2015). In addition, a case-control study showed that regular probiotic intake reduced the risk of bladder cancer in the healthy population (Ohashi et al., 2002). Taken together, these results strongly support the hypothesis that microbiome might be involved in bladder carcinogenesis, progression and relapse.

However, despite mounting researches on the human microbiome have yielded multiple insights into health and disease including cancers (Thomas et al., 2016), a detailed and comprehensive analysis of microbiota in urine of bladder cancer has not been yet reported. A recent study suggested that microbiome may be a factor in bladder cancer pathology and further studies on the urinary microbiota of bladder cancer would direct urologists to new therapeutic and prognostic options (Bucevic Popovic et al., 2017). Our primary purpose was to characterize urinary microbiota associated with bladder cancer in China and to explore the role of microbiome in bladder carcinogenesis.

Materials and Methods

Subject Recruitment and Specimen Collection

Urine specimens were collected from male patients with bladder cancer and non-neoplastic patients admitted to Nanfang Hospital in China between March 2017 and September 2017. All cancer cases were histologically confirmed as urothelial carcinoma and male controls were cases for a wide spectrum of non-neoplastic conditions, such as renal cyst. Subjects with prior known sexually transmitted infection or a recent history of urinary tract infections or antibiotic usage for any indication (within 1 month) were excluded. All subjects were required to finish a structured questionnaire to collect information on socio-demographic characteristics. Data collection followed the principles outlined in the Declaration of Helsinki. All participants had signed a written informed consent to contribute their own anonymous information to this study. Our study was approved by the Medicine Institutional Review Board of Southern Medical University.

Mid-stream urine specimens were collected by the clean catch method under the guidance of urotherapy nurses, then centrifuged at 16,000 g for 10 min immediately and stored at −80°C until further processing.

DNA Isolation and 16s rRNA Gene Sequencing

To avoid contamination, DNA isolation was performed using the cultured cells protocol supplied with the DNeasy Blood and Tissue Kit (Qiagen, Germany) in a laminar flow hood. The concentration of extracted DNA was determined through a Nanodrop ND-1000 spectrophotometer (Thermo Electron Corporation, USA). The genomic DNA isolated from the clinical samples was amplified using primer sets specific for V4 regions (515F: GTGCCAGCMGCCGCGGTAA; and 806R:GGACTACHVGGGTWTCTAAT). In order to evaluate contribution of extraneous DNA from reagents, extraction negative controls (no urine) and PCR negative controls (no template) were included. The resultant PCR products were purified by Qiaquick PCR purification kit (Qiagen, Valencia, CA). Finally, purified samples were normalized to equal DNA concentration and sequenced using the Illumina Miseq sequencer (Illumina, Inc., USA). The 16S rRNA gene sequences have been submitted to the Short Read Archive (SRA) under accession number SUB3915640.

Bioinformatics Analysis

Raw data were filtered to eliminate reads with adapter pollution and low quality to obtain clean reads by using QIIME (Caporaso et al., 2010). Filtered sequences were clustered by 97% identity into operational taxonomic units (OTUs) using UPARSE (Edgar, 2013), and subsequently, a single representative sequence from each clustered OTU was used to align to the SILVA database (Quast et al., 2013) and the Greengenes database (DeSantis et al., 2006) by Ribosomal Database Project Classifier (Wang et al., 2007).

QIIME was used to evaluate alpha diversity, which is composed of the Observed Species, Chao1, Shannon, Simpson and Ace indexes. Among them, the Observed Species, Chao1 and Ace indexes are indicators of species richness, while Shannon and Simpson indexes are indicators of species diversity. The difference of alpha diversity between groups was evaluated by Wilcoxon Rank-Sum Test (group number = 2) and Kruskal–Wallis test (n > 2) using SPSS (version 22).

To compare microbial composition between groups, beta diversity was evaluated by calculating the Bray Curtis, weighted UniFrac and unweighted UniFrac distances. Principal coordinate analysis (PCoA) was applied to generate three-dimensional plots in QIIME based on these distance matrices. The PERMANOVA was performed to test for statistical significance between groups using 999 permutations in QIIME.

Differential abundance analysis between groups was performed using Metastats and P-values were adjusted for multiple hypothesis testing using the False Discovery Rate based on the Benjamini-Hochberg (White et al., 2009). To identify significantly different bacteria between groups, taxa summaries were reformatted and input into Linear discriminant analysis effect size (LEfSe) via the Huttenhower Lab Galaxy Server (Segata et al., 2011). The Kruskal-Wallis rank sum test and Wilcoxon test were used to identify biomarkers, and linear discriminant analysis (LDA) was used to score them. Only taxa with logarithmic LDA score greater than 2 at a P < 0.05 were considered significantly enriched. To predict the functional pathways from microbiota composition data, Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) was performed for reconstruction of metagenome (Langille et al., 2013). Predicted functional genes were categorized into Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology and compared across patient groups using STAMP (version 2.1.3, Parks et al., 2014).

Statistical Analysis

Data are presented as median (first quartile to the third quartile) for continuous variables or number of cases (%) for counts data. The statistical significance of differences between groups were evaluated using Mann-Whitney U-test for continuous variables and Pearson's chi-square test or Fisher's Exact Test for count data through SPSS software (Version 22.0). All tests were two sided and P < 0.05 were considered statistically significant. PASS programme (PASS 11, NCSS, Kaysville UT, USA) provides estimates of power by simulation. Estimates were obtained 2000 simulations. We estimated that with 31 patients in cancer group and 18 subjects in control group, we would have 90% power to detect differences at the 0.05 significance level (alpha) using a two-sided Mann-Whitney Test.

Results

Subjects and Samples Characteristics

A total of 60 samples were analyzed, including 35 from male cancer patients and 25 from male non-neoplastic controls, while 4 samples in cancer group and 7 samples in control group were excluded for samples with too little sequencing reads (Supplementary Table 1). Detailed information on clinical characteristics and pathological parameters can be found in Table 1 and Supplementary Table 1. No significant difference was observed in the demographic characteristics between cancer and non-cancer group, except cigarette smoking (Table 1).

TABLE 1

Table 1. Comparisons of demographic characteristics and parameter of alpha diversity between cancer patients and non-cancer controls.

Cancer group was composed of 26 patients with non-muscle-invasive bladder cancer (NMIBC) and 5 patients with muscle-invasive bladder cancer (MIBC). The 2004 WHO grading system categorized cancer group as papillary urothelial neoplasm of low malignant potential (PUNLMP, n = 5) or low-grade papillary urothelial carcinoma group (LG, n = 11) or high-grade papillary urothelial carcinoma group (HG, n = 15) (Supplementary Table 2). Based on European Organization for Research and Treatment of Cancer (EORTC) scoring system (Lughezzani et al., 2011) (Supplementary Table 3), NMIBC group were stratified into lower risk of recurrence group (LER, recurrence score of EORTC ≤ 4, n = 16) and higher risk of recurrence group (HER, recurrence score of EORTC ≥ 5, n = 10), lower risk of progression group (LEP, progression score of EORTC ≤ 6, n = 15) and higher risk of progression group (HEP, progression score of EORTC ≥ 7, n = 11), respectively (Supplementary Table 2).

Sequencing Data, Alpha, and Beta Diversity

A total of 4,427,184 clean reads were obtained from the 49 samples. The median number of reads in cancer patients was 105,440, and in the non-cancer patients was 98,289 (Table 1, P > 0.05). The reads were classified into 1,653 OTUs that were used for downstream analysis. More OTUs were identified in urine from cancer patients, with an average of 120 OTUs in cancer group and 56 OTUs in control group (P = 0.008). Higher Observed Species index (richness, P = 0.008), Chao1 index (richness, P = 0.008), Ace index (richness, P = 0.003), Shannon index (diversity, P > 0.05) and lower Simpson index (diversity, P > 0.05) were presented in cancer group, which indicates that bacterial richness significantly increased in cancer patients while difference of species diversity was not significant (Figures 1A–E, Table 1).

FIGURE 1

Figure 1. Alpha and principal coordinate analysis (PCoA) for bladder-cancer samples and non-cancer samples. (A–E) Box-plot showing alpha diversity in samples using different metrics (A, observed species index; B, Chao1 index; C, Ace index; D, Shannon index; E, Simpson index). (F) PCoA plots of unweighed UniFrac distances in which samples were colored by clinical outcome. The PERMANOVA performed on the unweighted UniFrac distances showed that the observed differences were statistically significant (999 permutations; F = 2.53; P < 0.001).

Of interest was that we found that bacterial richness increased in HER group and HEP group, compared to LER group and LEP group, respectively (Supplementary Figures 1A–E, 2A–E; Supplementary Table 4), though no significant association was found between alpha diversity and tumor grade (Supplementary Figures 3A–E; Supplementary Table 5).

A comparison of urine from cancer patients and that from non-cancer controls showed significantly different bacterial profiles on unweighted UniFrac PCoA plots (Figure 1F). The PERMANOVA performed on the data set showed that the observed differences were statistically significant (F = 1.97, P < 0.05; F = 2.53, P < 0.001; and F = 1.61, P < 0.05, for weighted UniFrac, unweighted UniFrac and Bray Curtis distances, respectively). In addition, a clear hierarchical clustering of cancer samples was observed on dendrogram based on unweighted UniFrac distance metric (Figure 2). We next studied whether the microbial profile was different among different risk level and tumor grade. The results showed that the microbiota composition of patients with high risk of recurrence and progression was significantly different from that of patients with low risk of recurrence and progression (PERMANOA, F = 2.31, P < 0.01, LER vs. HER; F = 2.40, P < 0.01, LEP vs. HEP, for weighted UniFrac distances). However, no significant difference in the microbiota profiles of bladder cancer was observed for tumor grade (Supplementary Figures 1F, 2F, 3F, Supplementary Table 6).

FIGURE 2

Figure 2. The urinary microbiota profile of participants. A clear hierarchical clustering of cancer samples was observed in the dendrogram at genus level (left; based on the unweighted UniFrac distance metric). In the histogram (right), each colored box represents a bacterial taxon and each bar, a subject. The height of a colored box represents the relative abundance of that organism within the sample. Bacterial genera with a relative abundance <0.5% and unclassified genera are grouped as “Other”.

Relative Abundance of Urinary Bacteria in Cancer and Control Samples

At phylum level, the urinary microbiota was dominated by Proteobacteria (39.7% cancer, 49.0% control) and Firmicutes (32.8% cancer, 28.1% control), followed by Actinobacteria (7.0% cancer, 6.2% control) and Bacteroidetes(3.9% cancer, 9.4% control) (Figure 3A, Table 2). The microbial composition of all samples at class, order and family level were demonstrated in Figures 3B–D. In addition, the genera compositions of all samples were demonstrated in Figure 2. Though bacterial relative abundance at phylum, class or order level differed between cancer and non-cancer group, no significant difference was found using Metastats algorithm after False Discovery Rate adjustment (Table 2). However, it is notable that Sphingobacteriaceae (P = 0.047) was significantly more abundant in patients and Thermoactinomycetaceae (P = 0.005) in control group at family level, while Acinetobacter (P = 0.048) were significantly more abundant in cancer patients and Serratia (P = 0.003), Proteus (P = 0.003), Laceyella (P = 0.003) in non-cancer group at genus level (Table 2).

FIGURE 3

Figure 3. Bacterial average relative abundance in bladder cancer and non-cancer samples. Average distribution of major taxa is represented by bar graphs. (A) phylum; (B) class; (C) order; (D) family. Each colored box represents a bacterial taxon and the height of a colored box represents the relative abundance of that organism within the sample. Bacterial genera with a relative abundance <1% and unclassified genera are grouped as “Other”.

TABLE 2

Table 2. Comparison of relative abundance of urinary microbiome between cancer group and control group at all taxonomic levels.

Specific Genera Associated With Bladder Cancer

The LEfSe, which allows for identifying specific taxa associated with cancer, showed significantly higher compositional abundances of Acinetobacter, Anaerococcus, Rubrobacter, Sphingobacterium, Atopostipes, Geobacillus in cancer patients and Serratia, Proteus, Roseomonas, Ruminiclostridium-6, and Eubacterium–xylanophilum in control group at genus level (Figures 4A,B). Further analysis showed that Acinetobacter (cancer vs. control; 31 of 31 vs. 15 of 18, P = 0.044), Anaerococcus (cancer vs. control; 19 of 31 vs. 5 of 18, P = 0.035), Rubrobacter (cancer vs. control; 30 of 31 vs. 11 of 18, P = 0.002), Sphingobacterium (cancer vs. control; 13 of 31 vs. 2 of 18, P = 0.024), Atopostipes (13 of 31 vs. 2 of 18, P = 0.024) and Geobacillus (cancer vs. control; 12 of 31 vs. 1 of 18, P = 0.017) were detected in more cancer group samples than control group samples.

FIGURE 4

Figure 4. Microbial taxa associated with bladder cancer. (A) Cladogram representation of the urinary microbial taxa associated with bladder cancer (red) and non-cancer (green). (B) Association of specific microbiota taxa with cancer group and non-cancer group by linear discriminant analysis effect size (LEfSe). Red indicates taxa enriched in cancer group and green indicates taxa enriched in non-cancer group.

Since distinct difference of microbial profile was observed among groups with different risk of recurrence and progression, we next identified the specific taxa associated with high risk of recurrence and progression by using LEfSe analysis. The results showed that 6 genera were overrepresented in patients with high risk of recurrence and 4 genera in patients with high risk of progression, including Herbaspirillum, Gemella, Bacteroides, Porphyrobacter, Faecalibacterium, Aeromonas in HER group and Herbaspirillum, Porphyrobacter, Bacteroides, Marmoricola in HEP group (Supplementary Figures 4, 5).

Potential Functional Pathways Associated With Bladder Cancer

To infer the functional pathways based on the microbial community profiles we utilized PICRUSt. Overall, the microbial profiles present in patients with bladder cancer and non-cancer could not be distinguished clearly based on their functions (Figure 5A). The predicted KEGG pathways significantly enriched in bladder cancer included Staphylococcus aureus infection, glycerolipid metabolism, retinol metabolism, ethylbenzene degradation and carotenoid biosynthesis (Figure 5B).

FIGURE 5

Figure 5. Potential functional pathways associated with bladder cancer. (A) Principal component analysis (PCA) plot comparing the metagenome predictions (for KEGG orthology using PICRUSt) of the microbiota of patients with bladder cancer and non-cancer. (B) Microbial pathways that were significantly differentially enriched between the cases and controls. Red indicates taxa enriched in cancer group and blue indicates taxa enriched in non-cancer group.

Discussion

In this study, we have characterized the urinary microbial profile of bladder cancer by using 16S rRNA gene sequencing and the results showed that bacterial richness significantly increased in bladder cancer patients. In addition, higher bacterial richness was presented in urine samples with higher risk of recurrence and progression, which suggests that higher bacterial richness may be a potential indicator of high risk of recurrence and progression of NMIBC.

A reduction in microbial diversity has now been considered as a feature of gut disease, such as ulcerative colitis, Crohn's disease and colorectal cancer (Lepage et al., 2011; Ahn et al., 2013; Gevers et al., 2014). However, no consistent changes in microbial diversity were found among urinary tract disorders. Increased microbial diversity was observed in urgency urinary incontinence (Pearce et al., 2015), reduced diversity was found in interstitial cystitis (Siddiqui et al., 2012) and overactive bladder (Wu et al., 2017), while no significant difference in microbial diversity was found in prostate cancer (Shrestha et al., 2018).

Similar to microbiome inhabiting the gut, inter-individual microbial community heterogeneity of the urinary tract is influenced by lots of genetic and environmental factors, such as spatial distribution or lifestyle. However, despite the distinct inter-individual differences between samples, a common microbial feature appears to emerge, as shown in the PCoA analysis that clustered cancer group and control group separately (Figure 1F), indicating a possibility on common dysbiosis associated to bladder cancer. Studies suggest that changes in microbial composition and function may contribute to carcinogenesis, tumor progression and dissemination at these sites (Schwabe and Jobin, 2013). Alteration of urinary microbiome might be also involved in the development and progression of bladder cancer.

Bucevic Popovic et al. reported that no significant difference was observed in overall microbial profiles between bladder cancer (n = 12) and control group (n = 11) (Bucevic Popovic et al., 2017). Number of cases, ethnicity, gender and age may be reasons of different results between our study and Bucevic Popovic's study.

Acinetobacter and Anaerococcus were two genera found in higher abundances in bladder cancer patients than in non-cancer group. Acinetobacter spp. was reported one of the most abundant Gram-negative bacteria isolated from the urine of cattle affected by urothelial tumors of the urinary bladder (Roperto et al., 2012). Acinetobacter is a complex genus associated with nosocomial infections, including urinary tract infections. The virulence factors of Acinetobacter baumannii, a species of Acinetobacter, identified to date were confirmed to be involved in biofilm formation, adherence and invasion of epithelial cells, bacterial dissemination by degrading phospholipids present at mucosal barrier, and escape from the host immune response (McConnell et al., 2013). As for Anaerococcus, it was reported as a member of the Gram-positive anaerobic cocci, which was able to induce inflammation, remodeling of extracellular matrix (ECM) and re-epithelialization (Murphy and Frick, 2013). Based on above analysis, we raise an intriguing possibility that the interplay of ECM and microbiome and concomitant inflammation might play a role in bladder carcinogenesis.

The link between chronic inflammation, microbiome and the initiation and progression of solid tumors has been established for various neoplastic diseases, especially colorectal cancer (Irrazábal et al., 2014). Transient inflammation is considered as part of body's immune defense against pathogen, but persistent inflammation could potentially contribute to the development of cancer (Ainsworth, 2017). One large epidemiological study reported that repeated, regular bouts of cystitis were associated with increased risk of bladder cancer (Vermeulen et al., 2015). Alfano et al. reported that bacteria produce proteases, acting as virulence factors with crucial roles in host tissue degradation, as well as immune system evasion and destruction of host physical barriers, further promoting inflammation, remodeling of ECM and the generation of oxygen radicals, which results in mutagenesis that may promote the onset and progression of cancer (Alfano et al., 2016).

In our study, the LEfSe analysis showed that bacterial taxa along Sphingobacteriaceae-to-Sphingobacterium lineage (Figures 4A,B) were enriched in cancer patients. Sphingobacterium are recognized as etiological agents of cystic fibrosis and urinary tract infections. Besides, ceramides and sphingophospholipids of S. spiritivorum, one species of the genus Sphingobacterium, could induce DNA fragmentation, caspase-3 activation, changes in morphology and cell cycle shortening (Lambiase, 2014). The PICRUSt showed that in comparison with controls, Staphylococcus aureus infection was increased in patients with bladder cancer. Staphylococcus aureus could produce various enzymes such as alkaline protease, elastase and phospholipase C, which could degrade ECM components, break down elastin, disrupt tight junction, damage tissue, and cleaves various bonds in phospholipids. Above analysis also raises the possibility that microbiome-mediated modifications of ECM and concomitant inflammation might play an role in the initiation and development of bladder cancer. However, this study cannot answer the question whether changes in the microbiome contribute to cancer or vice versa. Numerous studies are required to determine whether the described profile is associated to, correlated with, or even responsible for bladder carcinogenesis and development.

Xu et al. reported enrichment of Streptococcus spp. in urine from urothelial carcinoma patients (n = 8) compared to healthy individuals (n = 6) (Xu et al., 2014). Bucevic Popovic et al. reported that genus Fusobacterium was significantly enriched in urine of bladder cancer patients. Streptococcus and Fusobacterium were also detected in our study. However, despite relative abundance of both genera was higher in bladder cancer patients, no significant difference was found between cancer and non-cancer group using Metastats algorithm after False Discovery Rate adjustment (Table 2).

Alfano et al. reported that specific pathogenic bacteria might promote the initiation and development of malignancies and tumor-associated microenvironments probably selectively facilitate the growth of specific bacteria in turn (Alfano et al., 2016). Enrichment of Herbaspirillum, Porphyrobacter, and Bacteroides was observed in bladder cancer patients with high risk of recurrence and progression (Supplementary Figures 4, 5), which suggests that these genera maybe potential biomarkers for risk stratification and to predict cancer prognosis. In addition, given the promising results of fecal microbiota transplantation, which was used to “re-establish the balance of nature” within the intestinal environment for the treatment of several refractory gastrointestinal disorders, strategies to restore normal bladder-associated microenvironment might be a potential option to reduce incidence or recurrence of bladder cancer.

Immune checkpoint targets have become the focus of investigation for the treatment of bladder cancer, including programmed death ligand-1 (PD-L1). Atezolizumab was the first PD-L1 inhibitor confirmed active in bladder cancer and is currently the only PD-L1 inhibitor approved by the FDA for patients with metastatic or locally advanced urothelial carcinoma based on promising overall response rates in clinical trials (Bellmunt et al., 2017). Though checkpoint blockade therapies have had remarkable results in management for bladder cancer, only small part of patients responds to PD-L1 blockers (Inman et al., 2017). It has become evident that microbiome plays a crucial role in modulating the response to cancer therapy, including immunotherapy (Roy and Trinchieri, 2017). Zitvogel et al reported that the microbiome could affect the therapeutic efficacy of PD-L1 blockade (Zitvogel et al., 2016). If the true association between microbial community and bladder carcinogenesis is confirmed in the future, urinary microbiome might be a target for enhancing bladder cancer responses PD-L1 inhibitor.

Since gender can affect microbiome studies (Kim et al., 2017) and men are at considerably higher risk of developing bladder cancer, this study included only male patients (Dobruch et al., 2016). On one hand the restriction to male patients limited the generalizability of the data, on the other hand the restriction minimized confounding factor of gender. In addition, to avoid potential bias caused by hospital noise and environmental changes, we recruited non-neoplastic patients as controls. However, urinary microbiota may be influenced by various disorders, despite excluding confounding factors of neoplasms.

Our study is not devoid of limitations. Firstly, it is not possible to determine the cause-effect relationship between microbiome and bladder cancer for retrospective study and low number of cases. Thus, prospective follow-up studies with a larger sample number and animal experiment studies will be needed to clarify the role of microbiome in development and progression of bladder cancer. Secondly, despite the fact that potential harm brought by catheterization can be avoided, mid-stream urine sample collected by the clean catch method is potentially contaminated with microbiota surrounding the urethral orifice. Thirdly, bladder mucosa would provide the best samples to characterize the bladder-associated microbiome, which could rule out sample contamination with microorganisms present in the urethra. However, obtaining bladder biopsies or suprapubic aspirates in healthy individuals is unethical. Hence, profiling microbiome of tumor, peri-tumor and non-tumor tissues of bladder cancer patients would be crucial to ascertain the association between pathogenic effect of microbial dysbiosis and the onset, progression and relapse of cancer. Another limitation is that we did not have relevant data on urinary incontinence or retention profiles. The effect of profiles on urination may contribute to the results in the study. Last but not least, the 16S rRNA gene sequencing-based approach enabled us to detect bacteria that present even in low numbers on one hand, on the other hand this method cannot identify bacteria well at species level and detect non-bacterial microorganisms, such as viruses and fungi (Ainsworth, 2017).

In conclusion, we have profiled the urinary microbiome associated with bladder cancer in the most comprehensive study to date. Our study suggests that urinary microbiota may be associated with bladder cancer, but the cause-effect relationship remains unclear. When a true association between urinary microbiota and bladder cancer is determined in the future, a better understanding of the role of microbiome in the development and progression of bladder cancer could provide novel diagnostic and prognostic biomarkers, as well as more microbiome-targeted therapeutic options.

Author Contributions

PW, GZ, JZ: conception and design; GZ, YC, JC: acquisition of data; GZ, YC, JZ: analysis and interpretation of data; PW, GZ, JZ: drafting of the manuscript; PW, WH, JRZ: critical revision of the manuscript for important intellectual content; PW, GZ, JZ, JLZ: statistical analysis; PW: obtaining funding; PW: supervision.

Funding

Supported by the Science and the Technology Planning Project of Guangdong Province (grant no. 2014A020212195), the Natural Science Foundation of Guangdong Province (grant no. 2016A030313605) and the National Natural Science Foundation of China (grant No.31401499).

Conflict of Interest Statement

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.

Supplementary Material

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

References

Ahn, J., Sinha, R., Pei, Z., Dominianni, C., Wu, J., Shi, J., et al. (2013). Human gut microbiome and risk for colorectal cancer. J. Natl. Cancer Inst. 105, 1907–1911. doi: 10.1093/jnci/djt300

PubMed Abstract | CrossRef Full Text | Google Scholar

Ainsworth, C. (2017). Microbiome: a bag of surprises. Nature 551, S40–S41. doi: 10.1038/551S40a

PubMed Abstract | CrossRef Full Text | Google Scholar

Alfano, M., Canducci, F., Nebuloni, M., Clementi, M., Montorsi, F., and Salonia, A. (2016). The interplay of extracellular matrix and microbiome in urothelial bladder cancer. Nat. Rev. Urol. 13, 77–90. doi: 10.1038/nrurol.2015.292

PubMed Abstract | CrossRef Full Text | Google Scholar

Aso, Y., Akaza, H., Kotake, T., Tsukamoto, T., Imai, K., and Naito, S. (1995). Preventive effect of a Lactobacillus casei preparation on the recurrence of superficial bladder cancer in a double-blind trial. BLP Stud. Group. Eur. Urol. 27, 104–109. doi: 10.1159/000475138

PubMed Abstract | CrossRef Full Text | Google Scholar

Babjuk, M., Böhle, A., Burger, M., Capoun, O., Cohen, D., Compérat, E. M., et al. (2017). EAU guidelines on non–muscle-invasive urothelial carcinoma of the bladder: update 2016. Eur. Urol. 71, 447–461. doi: 10.1016/j.eururo.2016.05.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Bellmunt, J., Powles, T., and Vogelzang, N. J. (2017). A review on the evolution of PD-1/PD-L1 immunotherapy for bladder cancer: the future is now. Cancer Treat. Rev. 54, 58–67. doi: 10.1016/j.ctrv.2017.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Boursi, B., Mamtani, R., Haynes, K., and Yang, Y. X. (2015). Recurrent antibiotic exposure may promote cancer formation–another step in understanding the role of the human microbiota? Eur. J. Cancer 51, 2655–2664. doi: 10.1016/j.ejca.2015.08.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Bucevic Popovic, V., Situm, M., Chow, C. T., Chan, L. S., Roje, B., and Terzic, J. (2017). The urinary microbiome associated with bladder cancer. BioRxiv doi: 10.1101/224667

CrossRef Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Chu, H., Khosravi, A., Kusumawardhani, I. P., Kwon, A. H. K., Vasconcelos, A. C., Cunha, L. D., et al. (2016). Gene-microbiota interactions contribute to the pathogenesis of inflammatory bowel disease. Science 352, 1116–1120. doi: 10.1126/science.aad9948

PubMed Abstract | CrossRef Full Text | Google Scholar

Costello, E. K., Stagaman, K., Dethlefsen, L., Bohannan, B. J., and Relman, D. A. (2012). The application of ecological theory toward an understanding of the human microbiome. Science 336, 1255–1262. doi: 10.1126/science.1224203

PubMed Abstract | CrossRef Full Text | Google Scholar

Darveau, R. P. (2010). Periodontitis: a polymicrobial disruption of host homeostasis. Nat. Rev. Microbiol 8, 481–490. doi: 10.1038/nrmicro2337

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

Dobruch, J., Daneshmand, S., Fisch, M., Lotan, Y., Noon, A. P., Resnick, M. J., et al. (2016). Gender and bladder cancer: a collaborative review of etiology, biology, and outcomes. Eur. Urol. 69, 300–310. doi: 10.1016/j.eururo.2015.08.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Dzutsev, A., Goldszmid, R. S., Viaud, S., Zitvogel, L., and Trinchieri, G. (2015). The role of the microbiota in inflammation, carcinogenesis, and cancer therapy. Eur. J. Immunol. 45, 17–31. doi: 10.1002/eji.201444972

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C (2013).UPARSE:highly accurate OUT sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604

PubMed Abstract | CrossRef Full Text

El-Omar, E. M., Oien, K., El-Nujumi, A., Gillen, D., Wirz, A., Dahill, S., et al. (1997). Helicobacter pylori infection and chronic gastric acid hyposecretion. Gastroenterology 113, 15–24. doi: 10.1016/S0016-5085(97)70075-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Gevers, D., Kugathasan, S., Denson, L. A., Vázquez-Baeza, Y., Van Treuren, W., Ren, B., et al. (2014). The treatment-naive microbiome in new-onset crohn's disease. Cell Host Microbe 15, 382–392. doi: 10.1016/j.chom.2014

PubMed Abstract | CrossRef Full Text | Google Scholar

Inman, B. A., Longo, T. A., Ramalingam, S., and Harrison, M. R. (2017). Atezolizumab: a PD-L1-blocking antibody for bladder cancer. Clin. Cancer Res. 23, 1886–1890. doi: 10.1158/1078-0432.CCR-16-1417

PubMed Abstract | CrossRef Full Text | Google Scholar

Irrazábal, T., Belcheva, A., Girardin, S. E., Martin, A., and Philpott, D. J. (2014). The multifaceted role of the intestinal microbiota in colon cancer. Mol. Cell 54, 309–320. doi: 10.1016/j.molcel.2014.03.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Kamat, A. M., Hahn, N. M., Efstathiou, J. A., Lerner, S. P., Malmstrom, P. U., Choi, W., et al. (2016). Bladder cancer. Lancet 388, 2796–2810. doi: 10.1016/S0140-6736(16)30512-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D., Hofstaedter, C. E., Zhao, C., Mattei, L., Tanes, C., Clarke, E., et al. (2017). Optimizing methods and dodging pitfalls in microbiome research. Microbiome 5:52. doi: 10.1186/s40168-017-0267-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Lambiase, A. (2014). “The family sphingobacteriaceae,” in The Prokaryotes: Other Major Lineages of Bacteria and The Archaea, eds E. Rosenberg, E.F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin, Heidelberg: Springer), 907–914.

Google Scholar

Langille, M. G., Zaneveld, J., Caporaso, J. G., McDonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. doi: 10.1038/nbt.2676

PubMed Abstract | CrossRef Full Text | Google Scholar

Lepage, P., Häsler, R., Spehlmann, M. E., Rehman, A., Zvirbliene, A., Begun, A., et al. (2011). Twin study indicates loss of interaction between microbiota and mucosa of patients with ulcerative colitis. Gastroenterology 141, 227–236. doi: 10.1053/j.gastro.2011.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, K., Lin, T., Xue, W., Mu, X., Xu, E., Yang, X., et al. (2015). Current status of diagnosis and treatment of bladder cancer in china – analyses of chinese bladder cancer consortium database. Asian J. Urol. 2, 63–69. doi: 10.1016/j.ajur.2015.04.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Lughezzani, G., Sun, M., Shariat, S. F., Budaus, L., Thuret, R., Jeldres, C., et al. (2011). A population-based competing-risks analysis of the survival of patients treated with radical cystectomy for bladder cancer. Cancer Am. Cancer Soc. 117, 103–109. doi: 10.1002/cncr.25345

PubMed Abstract | CrossRef Full Text | Google Scholar

McConnell, M. J., Actis, L., and Pachón, J. (2013). Acinetobacter baumannii: human infections, factors contributing to pathogenesis and animal models. FEMS Microbiol. Rev. 37, 130–155. doi: 10.1111/j.1574-6976.2012.00344.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Murphy, E. C., and Frick, I. M. (2013). Gram-positive anaerobic cocci–commensals and opportunistic pathogens. FEMS Microbiol. Rev. 37, 520–553. doi: 10.1111/1574-6976.12005

PubMed Abstract | CrossRef Full Text | Google Scholar

Ohashi, Y., Nakai, S., Tsukamoto, T., Masumori, N., Akaza, H., Miyanaga, N., et al. (2002). Habitual intake of lactic acid bacteria and risk reduction of bladder cancer. Urol. Int. 68, 273–280. doi: 10.1159/000058450

PubMed Abstract | CrossRef Full Text | Google Scholar

Parks, D. H., Tyson, G. W., Hugenholtz, P., and Beiko, R. G. (2014). STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics 30, 3123–3124. doi: 10.1093/bioinformatics/btu494

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-14. doi: 10.1128/mBio.01283-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Pearce, M. M., Zilliox, M. J., Rosenfeld, A. B., Thomas-White, K. J., Richter, H. E., Nager, C. W., et al. (2015). The female urinary microbiome in urgency urinary incontinence. Am. J. Obstet. Gynecol. 213, e1–11. doi: 10.1016/j.ajog.2015.07.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219

PubMed Abstract | CrossRef Full Text | Google Scholar

Roperto, S., Di Guardo, G., Leonardi, L., Pagnini, U., Manco, E., Paciello, O., et al. (2012). Bacterial isolates from the urine of cattle affected by urothelial tumors of the urinary bladder. Res. Vet. Sci. 93, 1361–1366. doi: 10.1016/j.rvsc.2012.06.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Roy, S., and Trinchieri, G. (2017). Microbiota: a key orchestrator of cancer therapy. Nat. Rev. Cancer 17, 271–285. doi: 10.1038/nrc.2017.13

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwabe, R. F., and Jobin, C. (2013). The microbiome and cancer. Nat. Rev. Cancer 13, 800–812. doi: 10.1038/nrc3610

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Shrestha, E., White, J. R., Yu, S. H., Kulac, I., Ertunc, O., De Marzo, A. M., et al. (2018). Profiling the urinary microbiome in men with positive versus negative biopsies for prostate cancer. J. Urol. 199, 161–171. doi: 10.1016/j.juro.2017.08.001

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

Thomas, A. M., Jesus, E. C., Lopes, A., Aguiar, S. J., Begnami, M. D., Rocha, R. M., et al. (2016). Tissue-associated bacterial alterations in rectal carcinoma patients revealed by 16S rRNA community profiling. Front. Cell. Infect. Microbiol. 6:179. doi: 10.3389/fcimb.2016.00179

PubMed Abstract | CrossRef Full Text | Google Scholar

Urbaniak, C., Gloor, G. B., Brackstone, M., Scott, L., Tangney, M., and Reid, G. (2016). The Microbiota of breast tissue and its association with breast cancer. Appl. Environ. Microbiol. 82, 5039–5048. doi: 10.1128/AEM.01235-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Vermeulen, S. H., Hanum, N., Grotenhuis, A. J., Castaño-Vinyals, G., van der Heijden, A. G., Aben, K. K., et al. (2015). Recurrent urinary tract infection and risk of bladder cancer in the Nijmegen bladder cancer study. Br. J. Cancer 112, 594–600. doi: 10.1038/bjc.2014.601

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07

PubMed Abstract | CrossRef Full Text | Google Scholar

White, J. R., Nagarajan, N., and Pop, M. (2009). Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput. Biol. 5:e1000352. doi: 10.1371/journal.pcbi.1000352

PubMed Abstract | CrossRef Full Text | Google Scholar

Whiteside, S. A., Razvi, H., Dave, S., Reid, G., and Burton, J. P. (2015). The microbiome of the urinary tract–a role beyond infection. Nat. Rev. Urol. 12, 81–90. doi: 10.1038/nrurol.2014.361

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, W., Yang, L., Lee, P., Huang, W. C., Nossa, C., Ma, Y., et al. (2014). Mini-review: perspective of the microbiome in the pathogenesis of urothelial carcinoma. Am. J. Clin. Exp. Urol. 2, 57–61.

PubMed Abstract | Google Scholar

Zitvogel, L., Ayyoub, M., Routy, B., and Kroemer, G. (2016). Microbiome and anticancer immunosurveillance. Cell 165, 276–287. doi: 10.1016/j.cell.2016.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Zitvogel, L., Daillère, R., Roberti, M. P., Routy, B., and Kroemer, G. (2017). Anticancer effects of the microbiome and its products. Nat. Rev. Microbiol. 15, 465–478. doi: 10.1038/nrmicro.2017.44

PubMed Abstract | CrossRef Full Text | Google Scholar