Wenzhi Du
Tianyi Wang
Wenxiu Zhang
Yu Xiao
Xinghuan Wang- 1Department of Urology, Zhongnan Hospital of Wuhan University, Wuhan University, Wuhan, China
- 2Department of Internal Medicine, Zhongnan Hospital of Wuhan University, Wuhan University, Wuhan, China
- 3Department of Pediatrics, Maternal and Child Health Care Hospital of Shandong Province, Jinan, China
- 4Department of Biological Repositories, Zhongnan Hospital of Wuhan University, Wuhan, China
Background: Observational studies have suggested a possible association between benign prostate hyperplasia (BPH) and bladder cancer (BLCA). However, these studies are prone to errors and limitations or confounding factors, making them unsuitable for assessing the causal relationship between BPH and BLCA.
Objective: Two-sample Mendelian randomization (MR) was performed to determine a possible association between genetically predicted BPH and the risk of BLCA.
Methods: A two-sample MR analysis was performed utilizing the Integrative Epidemiology Unit genome-wide association (GWAS) database of the Medical Research Council, United Kingdom A series of control steps, including five primary methods, were performed to identify the most suitable instrumental variables (IVs) for MR analysis. Sensitivity analysis was conducted to avoid statistical errors, including heterogeneity and pleiotropic bias.
Results: Genetic variants associated with BPH (P < 5 × 10–8) and BLCA (P < 5 × 10–6) were identified as instrumental variables and assessed using GWAS summary data (BPH, 4,670 cases vs. 458,340 controls; BLCA, 1,279 cases vs. 372,016 controls). BPH exhibited a positive effect on the occurrence of BLCA (inverse variance weighted (IVW), odds ratio (OR) = 1.095, 95% confidence interval (CI) = 1.030–1.165, p = 0.003), but there was no causal effect for BLCA on BPH (IVW, OR = 1.092, 95% CI = 0.814–1.465, p = 0.554).
Conclusion: Genetically predicted BPH was associated with a higher risk of BLCA in all histological subtypes. In contrast, the evidence was not significant to back the causality of genetically induced BLCA on BPH. These findings indicate that BPH plays a key role in developing BLCA in the European population. Further studies are needed to uncover the underlying mechanisms.
Introduction
Benign prostate hyperplasia (BPH), obstructive disease of the lower urinary system, is more prevalent among men of age more than 50 years, accompanied mainly by lower urinary tract symptoms (LUTS) (Chughtai et al., 2016), subsequently leading to bladder outflow obstruction (BOO), impairing urine emptying, producing residual urine (PVR) (Roehrborn, 2008; Foo, 2019), and reducing the quality of life. Long-term BOO can lead to bladder fibrosis, causing decreased bladder compliance and changes in postganglionic cholinergic fibers (Chertin et al., 2003; Wang et al., 2021a). In the early stage, the symptoms are primarily atypical because of forced urinary muscle function compensation. As the lower urinary tract obstruction gradually aggravates, the clinical symptoms become increasingly apparent, with the main manifestations being frequent urination, urgency, dysuria, increased nocturia, and hematuria (Irwin et al., 2011). Several studies have uncovered genetic determinants or biomarkers of BPH risk to provide essential guidance for novel pharmacotherapies or risk prediction (Sengupta et al., 2008; Zeng et al., 2009; Desaulniers et al., 2017; Hellwege et al., 2019).
Bladder cancer (BLCA) is the most common cancer of the urinary system, ranking 13th in mortality among all malignant tumors worldwide (Mao et al., 2021; Xiong et al., 2021). The number of new cases of BLCA in the US has increased to more than 80,000, with global mortality of 200,000 (Lenis et al., 2020). The development of BLCA is closely associated with genetic abnormalities and epigenetic changes. Tobacco use (Freedman et al., 2011), chemical exposure (Cumberbatch et al., 2018), radiation therapy, and chemotherapy are the leading causes of BLCA. Although both BPH and BLCA are independent diseases with evidently different pathophysiological changes, epidemiology, and risk factors (Chughtai et al., 2016; Lenis et al., 2020), the coexistence of BPH and BLCA has often been observed in many patients (Kang et al., 2007; Tseng, 2013). Nonetheless, there is little evidence to establish a definite causality between the two diseases. According to some studies, the risk of BLCA increases with BPH (Dai et al., 2016; Ham et al., 2021), and patients with a family history of BLCA may be more likely to have BPH (Negri et al., 2005). However, as an unavoidable inherent disadvantage, traditional observational studies have several limitations that weaken the reliability of their findings regarding causal effects (Thiese, 2014). Thus, there is a need to establish the causal connection between the two diseases to understand the etiology of disease processes and identify a novel therapeutic modality.
Mendelian randomization (MR) can estimate the exposure and outcome of a disease by estimating its genetics, especially in cases with insufficient evidence to conduct a randomized controlled trial (Duan et al., 2021). For MR, researchers use instrumental variables (IVs) comprising single-nucleotide polymorphisms (SNPs) (Noordam et al., 2016) to perform causal analysis, making MR more resistant to reverse causality, bias, and confounders compared to traditional observational studies. This method has been widely used in epidemiological studies to show the causal link between conditions and diseases (Wang et al., 2021b). The application of MR should be based on critical assumptions that—despite confounding factors—genetic variants are highly correlated with exposure and associated with the outcome only after exposure (Wang et al., 2022). Instead of relying on genetic variants from individuals, MR estimates the exposure-outcome relationship using a summary of genome-wide association (GWAS) data (Zhao et al., 2021). No MR study has examined the causality between BPH and BLCA. Based on all the above-mentioned information, we hypothesized that while BPH can increase the risk of BLCA, but the latter does not exacerbate the condition of BPH. In this study, we attempted to identify the causal relationships and direction between the two diseases by performing a two-sample bidirectional MR analysis.
Methods
Source of data
For this study, data were collected from two large GWAS datasets. When there is a variation between the allele frequencies and exposure distributions of different ethnic groups, the result of population stratification is an association between the outcome and exposure (Larsson et al., 2019). To avoid population stratification bias, we selected genetic variants only from European ancestries, which ensured that the results did not reflect the racial or ethnic composition of the population. The Integrative Epidemiology Unit (IEU) has created the IEU GWAS database at the University of Bristol (Duan et al., 2021). It comprises over 40,000 publicly available GWAS datasets, continuously updated and used for academic research, including studies concerning prostate (Lv et al., 2022) and bladder (Chen et al., 2021). To assess the causal relationship between BPH and BLCA, we downloaded publicly available datasets from IEU summary-level GWAS data (https://gwas.mrcieu.ac.uk/) and performed two-sample bidirectional MR estimates, including BPH (ukb-b-11601) and BLCA (ieu-b-4874). The BPH dataset derives from the IEU consortium of the Medical Research Council (MRC-IEU), including 463,010 individuals (4,670 cases vs. 458,340 controls) and 9,851,867 SNPs of European ancestry. The GWAS dataset for BLCA comprises 373,295 (1,279 cases vs. 372,016 controls) and 9,904,926 SNPs of European ancestry. As mentioned above, we only adopted European ancestry datasets to prevent deviations from population stratification. (Table 1).
TABLE 1. Basic information of datasets.
Single nucleotide polymorphism selection
A range of quality control standards was established using the above-mentioned datasets, and the specific IVs having solid associations with exposure were identified. First, the identified IVs needed to have a genome-wide significance level. We selected IVs with P < 5 × 10−8 from exposures, which can guarantee a significant association of IVs with exposures (BPH or BLCA). In our BLCA data analysis, we detected only two genome-wide significant SNPs (rs146190477 and rs74903734), which may result in weak instruments (Choi et al., 2019). Therefore, we applied a set of genetic instruments to reach a more relaxed threshold (P < 5 × 10−6). Then, to meet the above-mentioned basic IV assumptions, we removed linkage disequilibria (LD) —a measure of the association between alleles at related loci (Slatkin, 2008) —to prevent the genetic variation that LD caused. Accordingly, the TwoSampleMR package default parameters were set (physical distance threshold, 10,000 kb; r2 < 0.001) for clumping data and shearing the genetic variants that could be possible LD. Then, we used the PhenoSCanner database, which includes a large amount of information about SNPs with possible linkages to diseases and phenotypes (Staley et al., 2016), and found no association of SNPs between BPH and BLCA. Finally, the selected SNPs were not considered biased due to their minor allele frequencies (MAF) (not < 0.01) in the original datasets. From the outcome study, if no proxy variant was identified for a target variant, it was excluded. Similarly, when using the TwoSampleMR package, a palindromic SNP with intermediate allele frequencies must be eliminated to harmonize the outcome and exposure. We evaluated the F statistics to analyze the intensity of the association between exposure and SNPs (set threshold value > 10; formula: F = β2/SE2; the effect on the risk of exposure, β; the standard error, SE) (Burgess and Thompson, 2011). One of the SNPs was considered a “weak instrument” (value < 10), and was thus removed from the analysis (Carter et al., 2021). The explained genetic variance was also calculated for each instrument (formula:
Analysis tool
Data were analyzed using R software (4.1.0) and the R packages by MR and MR-Pleiotropy RESidual Sum and Outlier (PRESSO) (Hemani et al., 2018; Verbanck et al., 2018). The code used is openly accessible.
Mendelian randomization analysis
For statistical analysis, we implemented five approaches, including inverse variance weighted (IVW), MR-Egger, weighted median, simple mode, and weighted mode. First, all selected SNPs were analyzed by the IVW method, which can efficiently find whether directional pleiotropy exists or not. Then, the slope of the weighted linear regression against the horizontal pleiotropy was used in the MR-Egger method to obtain a causal estimate (Burgess and Thompson, 2017). Then, the weighted median estimate was used when 50% of the weight was derived from valid IVs (Bowden et al., 2016a). To further reflect the authenticity of our MR analysis, we adopted the simple mode, weighted modes, and the maximum likelihood ratio test. Lastly, considering that exposure and outcome can interact as both cause and effect, we reversed the role of exposure and outcome in bidirectional MR studies and further established the causal relationship between BPH and BLCA.
Sensitivity analysis
As a barometer of the MR assumption violations (Bowden et al., 2016b), the degree of heterogeneity of the IVs was evaluated using the IVW and ER-Egger methods. We implemented Cochran’s Q test to quantify the heterogeneity when p < 0.05 reflected significant heterogeneity among the causal appraisal for the selected SNPs (Cohen et al., 2015). Moreover, to exclude the possibility that any specific SNP may change the effect proportionately, we performed a leave-one-out plot analysis (Cheng et al., 2017) to determine whether no such distinct SNP exists; this method remarkably decides the confidence and scientificalness of the MR analysis. We employed funnel plots, MR-PRESSO, and MR-Egger to find the potential unbalanced horizontal pleiotropy between the IVs. MR-Egger intercept tool measures the likelihood of a directional pleiotropy of a genetic variant (p < 0.05), and the MR-PRESSO test can detect the potential pleiotropy even in a tiny number of loci (Burgess and Thompson, 2017; Hemani et al., 2018). In addition, compared with the MR-Egger regression test, the MR-PRESSO outliner test has higher accuracy (Hemani et al., 2018). In our bidirectional MR analysis, all sensitivity tests mentioned above are two-sided.
Results
Benign prostate hyperplasia increases the risk of bladder cancer
From the GWAS dataset for BPH without LD (r2 < 0.001), 12 IVs (Table 2 and Supplementary Figure S3) were selected; these were within the 10,000 kb physical distance threshold and achieved genome-wide significance (P < 5 × 10−8). To avoid SNPs with palindromic nature and intermediate allele frequencies, we abolished one SNP (rs12131120) with the help of the two-sample MR function of the R package. We conducted the F statistic to ensure that no weak instrument exists, and the results were in the range of 37.1–417.4.
TABLE 2. MR effect size of BPH on BLCA with details of selected SNPs.
We selected 12 compelling IVs above, we use IVW (p = 0.003 < 0.05), MR-Egger (p = 0.103 > 0.05), weighted median (p = 0.099 > 0.05), simple mode (p = 0.497 > 0.05), and weighted mode (p = 0.521 > 0.05) methods to analyze the causal effects of BPH and BLCA (Table 3 and Supplementary Figures S1B,C). The IVW method indicated a symbolic association; nonetheless, the results of the five above-mentioned approaches showed no consistency. Therefore, to further establish the credibility of our conclusion, we performed the maximum likelihood method (p = 0.003 < 0.05), which complemented the outcomes of the IVW approach and verified the persuasiveness of our study.
TABLE 3. BPH & BLCA: Bidirectional MR results.
Assessment of the heterogeneity of the selected IVs using the MR-Egger (Q = 5.508; p = 0.598 > 0.05) and IVW (Q = 7.613; p = 0.472 > 0.05) methods (Table 4) showed that no apparent heterogeneity. Leave-one-out cross-validation was conducted, and no specific SNP was observed to influence the effects disproportionately (Supplementary Figure S1A). The MR-Egger horizontal pleiotropy test suggested the absence of any bias in the results (Intercept = −0.00052, p = 0.190 > 0.05). Besides, the funnel plot test revealed the absence of horizontal pleiotropy (Table 4 and Supplementary Figure S1D). Moreover, MR-PRESSO global test and outliner test did not detect any heterogeneity or outliners.
TABLE 4. Heterogeneity and horizontal pleiotropy check.
Bladder cancer does not increase the risk of benign prostate hyperplasia
We reversed the outcome and exposure location in the bi-direction MR analysis and echoed the above-mentioned causal effect MR analysis approaches of BPH to BLCA. From the GWAS dataset for BLCA with no LD (r2 < 0.001), we selected 25 IVs (Table 5 and Supplementary Figure S3) that were within the 10,000 kb physical distance threshold, had acquired genome-wide significance (P < 5 × 10–6) and had no potential pleiotropy as tested by the PhenoScanner tool. These IVs were found in BLCA within the MAF threshold (> 0.01). There were no palindromic SNPs with intermediate allele frequencies. Individual SNPs had F statistics in the range of 36.4 and 401.3.
TABLE 5. MR effect size of BLCA on BPH with details of selected SNPs.
We repeated IVW (p = 0.554 > 0.05), MR-Egger (p = 0.732 > 0.05), weighted median (p = 0.517 > 0.05), simple mode (p = 0.464 > 0.05), weighted mode (p = 0.500 > 0.05), and found no persuasive clue of BLCA increasing the risk of BPH (Table 3 and Supplementary Figures S2B–C).
Subsequently, although the heterogeneity was detected in the selected SNPs through the IVW (Q = 34.512, p = 0.0005 < 0.05) and MR-Egger (Q = 33.852, p = 0.0003 < 0.05) estimation (Table 4), the leave-one-out cross-validation revealed no SNPs that disproportionally influence the effects (Supplementary Figure S2A); besides, the MR-Egger horizontal pleiotropy test (p = 0.652 > 0.05) and funnel plot test showed no bias in the results (Table 4 and Supplementary Figure S2D). Furthermore, the MR-PRESSO test also showed no heterogeneity or outliners. In summary, we have reached a stable and reliable conclusion.
Discussion
Males suffering from LUTS often have BPH, characterized by proliferation features in the central zone of epithelial components and stroma of the prostate (Jahan et al., 2021). The overall annual incidence of LUTS induced by BPH is estimated to be 15/1,000 person-years (95% CI = 14.8–16.1). In individuals aged 45–49 years, the incidence is 3/1,000 person-years, and in those aged 75–79 years, the incidence increases to 38/1,000 person-years (Naderi et al., 2004). More than 30% of patients with severe LUTS need medical or surgical care (Madersbacher et al., 2019). Long-term BOO can cause pathophysiological changes in the bladder. Many experts speculate a possible causal relationship between BPH and BLCA, a highly malignant disease with high recurrence and progression rates. BLCA has a low 5-year cancer-specific survival rate (Yang et al., 2018), causing 550,000 cases worldwide and 17,000 deaths in the US alone annually (Siegel et al., 2019; Richters et al., 2020).
Nonetheless, the published evidence showed no consistency in the causality of BPH and BLCA. According to a few studies, BPH may increase the likelihood of developing BLCA (Fang et al., 2018; Lenis et al., 2020) and is a crucial risk factor for BLCA, especially for patients with type 2 diabetes (Tseng, 2013). In one meta-analysis that systematically reviewed 16 case-control studies and 10 cohort studies, BPH was associated with a higher incidence of BLCA (18). To reduce prostate volume and lower the recurrence rate of BLCA, some experts used 5-α reductase inhibitors or other androgen deprivation therapies (Shiota et al., 2017; Zhu et al., 2021). Besides, simultaneous surgical treatment to cure BPH and BLCA is safe in patients having both diseases; the process adds no greater risk of BLCA recurrence or metastasis rates, and patients enjoy an improved quality of life after the procedure (Zhou et al., 2020; Wang et al., 2021c). However, some researchers found that BPH has no positive effect on the occurrence of BLCA. In a study on a cohort of almost 80,000 Swedish men with BPH, no association was found between BPH and BLCA, except for those patients who underwent transurethral prostate resection (Kang et al., 2007). Some researchers have also found that transurethral resection of the prostate cannot lower the incidence of bladder cancer (Lin et al., 2018). Nonetheless, these conventional observational studies showed the influence of BPH on BLCA and remained vulnerable to the threat of potential confounders, such as smoking, occupation, and heredity, that are unavoidable in crowd screening. Besides, researchers may find chance results in multiple comparisons and subgroup analyses in a cohort study. As these two common diseases afflict middle-aged and older men, it is crucial to understand the causal relationship between them for better clinical management.
Genetic information is always transmitted unidirectionally from parent to offspring, limiting the impact of reverse causation (Swerdlow et al., 2016). At this point, traditional epidemiological studies are equally unavoidable. To further clarify the inconsistencies and provide a valuable guideline for clinical management, we performed the first two-sample and bidirectional MR analysis to assess the causal direction between BPH and BLCA. To establish our hypothesis, we performed an array of quality control approaches and MR analytical methods and tested heterogeneity; any potential pleiotropy (P > 0.05) was not detected in this MR-based bidirectional causality study. The results demonstrate a causal effect between BPH and BLCA, wherein BPH can increase the risk of BLCA, but the opposite was not observed. According to the findings of this study, we suggest some recommendations for patients with BPH or BLCA. First, patients, especially those with long-term BPH and severe LUTS, need to be examined timely and periodically to prevent the occurrence of BLCA and acquire a better living status. Second, patients who suffer from both BPH and BLCA are suggested to adopt early intervention to avert the development of BPH and slow down the malignant process of BLCA. Third, patients with BLCA but not BPH need to be attentive to any prostate condition during diagnosis and treatment to prevent the aggravation of BLCA due to the later occurrence of BPH. Finally, patients with BLCA experiencing urinary tract irritation should seek medical attention to distinguish whether the source of this symptom comes from BLCA alone or combined BPH.
Our research findings are consistent with the conclusions of previous studies. There are many mechanisms by which BPH leads to BLCA. The ratio of acid phosphatase to β-glucuronidase in the urine of patients with prostate cancer was only 10% effective in those having BPH (Muntzing and Nilsson, 1972). BPH causes lower urinary tract obstruction, poor urination, urinary retention, and increased β-glucuronidase activity in urine, which decomposes β-naphthylamine, 4-aminobiphenyl, and benzidine, that do not easily and directly cause carcinogenesis by being converted into 2-amino-1 naphthol that is a carcinogen, thus significantly prolonging the contact time between some carcinogens and bladder mucosa (Melicow et al., 1961). Matsumoto et al. performed an animal study and observed that BOO prolonged exposure to carcinogens and caused a higher risk of BLCA (Matsumoto et al., 2009). Besides, long-term urinary retention can breed bacteria and increase the risk of cystitis and urinary tract infections, especially for patients with diabetes mellitus (Selius and Subedi, 2008; Odabasi and Mert, 2020). With the deposition of crystalline materials and chemical substances (Coe et al., 2016; Barghouthy et al., 2021), bladder stones are easily formed, consequently aggravating bladder inflammation (Ma et al., 2016; Cicione et al., 2018). These complications induce the proliferation or metaplasia of bladder transitional epithelial cells, causing the formation of transitional cell carcinoma and increasing the chance of recurrence. In one review on BLCA diagnosis and management, researchers reported that the bladder diverticulum has a higher risk of developing BLCA, which might be caused by prolonged exposure of the bladder mucosa to carcinogenic substances due to urinary stasis (Walker et al., 2014). Diluting urinary metabolites by taking more water, removing BOO, and increasing the frequency of urination is associated with a lower incidence of BLCA (Zhou et al., 2014). Chronic kidney disease (CKD) is significantly associated with BPH; patients with CKD have a higher risk of developing BLCA than either patients with BPH or CKD only (Chung et al., 2007; Chen et al., 2008). The potential mechanism may be that CKD tends to be accompanied by obstructive urinary tract disease (Rassweiler et al., 2006). In addition, androgens have a high affinity to prostate tissue through the androgen receptors (AR), which can significantly influence the physiological state of prostate tissues and facilitate BPH development (Heinlein and Chang, 2004; Tong and Zhou, 2020), and trigger the occurrence of various bladder carcinomas (Wang et al., 2020a; Martinez-Rojo et al., 2021). Patients with BLCA who received androgen suppression medicine were more likely to survive than those who did not (Wang et al., 2020b; Wu et al., 2021). Both female and male BLCA cells are highly dependent on AR (Sottnik et al., 2021). AR expression or dihydrotestosterone (DHT) treatment can decrease the forkhead box O1 (FOXO1) expression and cause FOXO1 phosphorylation through the protein kinase B signaling pathway, which works as a tumor suppressor to BLCA cells metastasis, migration, invasion, and proliferation (Ide et al., 2020). The AR-FOXO1-AKT axis might be an essential pathway to the development of BLCA.
Our study had several strengths which are listed as follows. First, this is the first MR study to systematically determine the causal relationship between BPH and the risk of developing BLCA in the European population. This study revealed that the occurrence of BPH caused by specific genetic variants could be a potential cause of BLCA and indicates that there exists a unidirectional genetic link between BPH and BLCA, compensating for the shortcomings of previous studies. Second, our two-sample MR analyses are consistent with the conclusions of previous observational studies, making the conclusions more reliable and trustworthy. However, a few limitations of our study warrant our attention. First, the lack of information regarding the degree of BPH and the BLCA stage prevented a subgroup analysis. Second, although the study supported a link between BPH and BLCA, the sample was only derived from the European population, which is not large enough to generalize various ancestries. Third, in correlation analysis with BPH, we could not exclude the sample of females in the BLCA group; this is caused by the limitations of the GWAS database and may cause statistical bias. Nevertheless, although such a potential effect exists, it is relatively small and unlikely to significantly influence our conclusion. Fourth, overlapping participants should not be involved in the exposure and outcome studies of two-sample MR analysis. In this study, we could not assess the extent of overlap between the two samples. However, the bias due to the overlap between the two samples can be minimized using more powerful tools estimated by the F statistic (much more significant than 10) (Pierce and Burgess, 2013; Wu et al., 2020). Fifth, we extracted data from the largest available GWAS of BLCA, and a few genome-wide significant SNPs were identified, which may have led to weak genetic instruments. Therefore, we applied statistical criteria to meet a more relaxed threshold to include more SNPs as instruments; this method has been performed in many MR studies to solve the problem of insufficient genome-wide significant SNPs (Gage et al., 2017; Hartwig et al., 2017). Last, we found positive results only in IVW but not through other MR statistical methods; nonetheless, the conclusion is still reliable because IVW is the most used, primary, and persuasive MR analysis method (Bowden et al., 2016b). This method combined the Wald ratios of exposure and outcome and provides the most accurate estimation for all selected SNPs (Burgess et al., 2013; Emdin et al., 2017). Meanwhile, we also performed a maximum likelihood ratio analysis, whose outcome was consistent with the IVW result.
Conclusion
This MR study demonstrated that BPH increased the risk of BLCA, but BLCA did not influence the overall risk of BPH. More studies are needed to explore the underlying biochemical mechanism and find potential therapeutic targets for the prevention of these diseases.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics statement
No additional ethical approval was required because publicly available summary-level data were adopted in the MR analysis. Each original publication can identify ethical approval and informed consent for each included study. All detailed studies followed the Declaration of Helsinki.
Author contributions
WD browsed and gathered the data, performed the MR analysis process using R packages, and wrote and edited the manuscript. WD, TW, and WZ supervised the MR analysis process and checked the data. XW and YX offered instructive and practical suggestions during the process of topic designing. All authors reviewed the manuscript, and their approval was obtained.
Funding
This study was supported by the Improvement Project for Theranostic Ability on Difficulty Miscellaneous Disease from the National Health Commission of China (ZLYNXM202006), the Science and Technology Department of Hubei Province Key Project (2022EJD001 and YYXKNL2022001) and the Fundamental Research Funds for the Central Universities (4206-413100049).
Acknowledgments
We gratefully acknowledge the exceptional technical assistance provided by Wei Liu and Zhiwen He from Zhongnan Hospital of Wuhan University. We want to acknowledge the GWAS database for providing free access to data.
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/fgene.2022.1016696/full#supplementary-material
Abbreviations
AR, androgen receptors; BLCA, bladder cancer; BOO, obstruction of bladder outflow; BPH, benign prostate hyperplasia; CI, confidence interval; CKD, chronic kidney disease; DHT, dihydrotestosterone; FOXO1, forkhead box O1; GWAS, genome-wide association study; IVW, inverse variance weighted; IVs, instrumental variables; LD, linkage disequilibria; LUTS, lower urinary tract symptoms; MR, mendelian randomization; MAF, minor allele frequencies; OR, odds ratio; PVR, producing residual urine; PV, prostate volume; SNPs, single nucleotide polymorphisms; VE, variance explained.
References
Barghouthy, Y., Corrales, M., Doizi, S., Somani, B. K., and Traxer, O. (2021). Tea and coffee consumption and pathophysiology related to kidney stone formation: A systematic review. World J. Urol. 39 (7), 2417–2426. doi:10.1007/s00345-020-03466-8
Bowden, J., Davey Smith, G., Haycock, P. C., and Burgess, S. (2016). Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet. Epidemiol. 40 (4), 304–314. doi:10.1002/gepi.21965
Bowden, J., Del Greco, M. F., Minelli, C., Davey Smith, G., Sheehan, N. A., and Thompson, J. R. (2016). Assessing the suitability of summary data for two-sample mendelian randomization analyses using MR-egger regression: The role of the I2 statistic. Int. J. Epidemiol. 45 (6), 1961–1974. doi:10.1093/ije/dyw220
Burgess, S., Butterworth, A., and Thompson, S. G. (2013). Mendelian randomization analysis with multiple genetic variants using summarized data. Genet. Epidemiol. 37 (7), 658–665. doi:10.1002/gepi.21758
Burgess, S., and Thompson, S. G. (2011). Avoiding bias from weak instruments in Mendelian randomization studies. Int. J. Epidemiol. 40 (3), 755–764. doi:10.1093/ije/dyr036
Burgess, S., and Thompson, S. G. (2017). Interpreting findings from Mendelian randomization using the MR-Egger method. Eur. J. Epidemiol. 32 (5), 377–389. doi:10.1007/s10654-017-0255-x
Carter, A. R., Sanderson, E., Hammerton, G., Richmond, R. C., Davey Smith, G., Heron, J., et al. (2021). Mendelian randomisation for mediation analysis: Current methods and challenges for implementation. Eur. J. Epidemiol. 36 (5), 465–478. doi:10.1007/s10654-021-00757-1
Chen, C. H., Shun, C. T., Huang, K. H., Huang, C. Y., Yu, H. J., and Pu, Y. S. (2008). Characteristics of female non-muscle-invasive bladder cancer in taiwan: Association with upper tract urothelial carcinoma and end-stage renal disease. Urology 71 (6), 1155–1160. doi:10.1016/j.urology.2007.11.140
Chen, M., Xu, Y., Xu, J., Chancoco, H., and Gu, J. (2021). Leukocyte telomere length and bladder cancer risk: A large case-control study and mendelian randomization analysis. Cancer Epidemiol. Biomarkers Prev. 30 (1), 203–209. doi:10.1158/1055-9965.EPI-20-0351
Cheng, H., Garrick, D. J., and Fernando, R. L. (2017). Efficient strategies for leave-one-out cross validation for genomic best linear unbiased prediction. J. Anim. Sci. Biotechnol. 8, 38. doi:10.1186/s40104-017-0164-6
Chertin, B., Rolle, U., Cascio, S., and Puri, P. (2003). Alterations in cholinergic and neuropeptide innervation of urinary bladder following partial bladder outlet obstruction. Pediatr. Surg. Int. 19 (6), 427–431. doi:10.1007/s00383-002-0937-6
Choi, K. W., Chen, C. Y., Stein, M. B., Klimentidis, Y. C., Wang, M. J., Koenen, K. C., et al. (2019). Assessment of bidirectional relationships between physical activity and depression among adults: A 2-sample mendelian randomization study. JAMA Psychiatry 76 (4), 399–408. doi:10.1001/jamapsychiatry.2018.4175
Chughtai, B., Forde, J. C., Thomas, D. D., Laor, L., Hossack, T., Woo, H. H., et al. (2016). Benign prostatic hyperplasia. Nat. Rev. Dis. Prim. 2, 16031. doi:10.1038/nrdp.2016.31
Chung, S. D., Huang, K. H., Lai, M. K., Huang, C. Y., Chen, C. H., Pu, Y. S., et al. (2007). CKD as a risk factor for bladder recurrence after nephroureterectomy for upper urinary tract urothelial carcinoma. Am. J. Kidney Dis. 50 (5), 743–753. doi:10.1053/j.ajkd.2007.08.007
Cicione, A. C. D. E. N., Manno, S., Damiano, R., Posti, A., Lima, E., et al. (2018). Bladder stone management: An update. Minerva Urol. Nefrol. 70 (1), 53–65. doi:10.23736/S0393-2249.17.02972-1
Coe, F. L., Worcester, E. M., and Evan, A. P. (2016). Idiopathic hypercalciuria and formation of calcium renal stones. Nat. Rev. Nephrol. 12 (9), 519–533. doi:10.1038/nrneph.2016.101
Cohen, J. F., Chalumeau, M., Cohen, R., Korevaar, D. A., Khoshnood, B., and Bossuyt, P. M. (2015). Cochran's Q test was useful to assess heterogeneity in likelihood ratios in studies of diagnostic accuracy. J. Clin. Epidemiol. 68 (3), 299–306. doi:10.1016/j.jclinepi.2014.09.005
Cumberbatch, M. G. K., Jubber, I., Black, P. C., Esperto, F., Figueroa, J. D., Kamat, A. M., et al. (2018). Epidemiology of bladder cancer: A systematic review and contemporary update of risk factors in 2018. Eur. Urol. 74 (6), 784–795. doi:10.1016/j.eururo.2018.09.001
Dai, X., Fang, X., Ma, Y., and Xianyu, J. (2016). Benign prostatic hyperplasia and the risk of prostate cancer and bladder cancer: A meta-analysis of observational studies. Med. Baltim. 95 (18), e3493. doi:10.1097/MD.0000000000003493
Desaulniers, A. T., Cederberg, R. A., Lents, C. A., and White, B. R. (2017). Expression and role of gonadotropin-releasing hormone 2 and its receptor in mammals. Front. Endocrinol. 8, 269. doi:10.3389/fendo.2017.00269
Duan, C., Shi, J., Yuan, G., Shou, X., Chen, T., Zhu, X., et al. (2021). Causal association between heart failure and alzheimer's disease: A two-sample bidirectional mendelian randomization study. Front. Genet. 12, 772343. doi:10.3389/fgene.2021.772343
Emdin, C. A., Khera, A. V., and Kathiresan, S. (2017). Mendelian randomization. JAMA 318 (19), 1925–1926. doi:10.1001/jama.2017.17219
Fang, C. W., Liao, C. H., Wu, S. C., and Muo, C. H. (2018). Association of benign prostatic hyperplasia and subsequent risk of bladder cancer: An asian population cohort study. World J. Urol. 36 (6), 931–938. doi:10.1007/s00345-018-2216-8
Foo, K. T. (2019). What is a disease? What is the disease clinical benign prostatic hyperplasia (BPH)? World J. Urol. 37 (7), 1293–1296. doi:10.1007/s00345-019-02691-0
Freedman, N. D., Silverman, D. T., Hollenbeck, A. R., Schatzkin, A., and Abnet, C. C. (2011). Association between smoking and risk of bladder cancer among men and women. JAMA 306 (7), 737–745. doi:10.1001/jama.2011.1142
Gage, S. H., Jones, H. J., Burgess, S., Bowden, J., Davey Smith, G., Zammit, S., et al. (2017). Assessing causality in associations between cannabis use and schizophrenia risk: A two-sample mendelian randomization study. Psychol. Med. 47 (5), 971–980. doi:10.1017/S0033291716003172
Ham, W. S., Park, J. S., Jang, W. S., Choi, Y. D., and Kim, J. (2021). Prognostic value of prostate volume in non-muscle invasive bladder cancer. Sci. Rep. 11 (1), 18784. doi:10.1038/s41598-021-98045-1
Hartwig, F. P., Borges, M. C., Horta, B. L., Bowden, J., and Davey Smith, G. (2017). Inflammatory biomarkers and risk of schizophrenia: A 2-sample mendelian randomization study. JAMA Psychiatry 74 (12), 1226–1233. doi:10.1001/jamapsychiatry.2017.3191
Heinlein, C. A., and Chang, C. (2004). Androgen receptor in prostate cancer. Endocr. Rev. 25 (2), 276–308. doi:10.1210/er.2002-0032
Hellwege, J. N., Stallings, S., Torstenson, E. S., Carroll, R., Borthwick, K. M., Brilliant, M. H., et al. (2019). Heritability and genome-wide association study of benign prostatic hyperplasia (BPH) in the eMERGE network. Sci. Rep. 9 (1), 6077. doi:10.1038/s41598-019-42427-z
Hemani, G., Zheng, J., Elsworth, B., Wade, K. H., Haberland, V., Baird, D., et al. (2018). The MR-Base platform supports systematic causal inference across the human phenome. Elife 7, e34408. doi:10.7554/eLife.34408
Ide, H., Mizushima, T., Jiang, G., Goto, T., Nagata, Y., Teramoto, Y., et al. (2020). FOXO1 as a tumor suppressor inactivated via AR/ERβ signals in urothelial cells. Endocr. Relat. Cancer 27 (4), 231–244. doi:10.1530/ERC-20-0004
Irwin, D. E., Kopp, Z. S., Agatep, B., Milsom, I., and Abrams, P. (2011). Worldwide prevalence estimates of lower urinary tract symptoms, overactive bladder, urinary incontinence and bladder outlet obstruction. BJU Int. 108 (7), 1132–1138. doi:10.1111/j.1464-410X.2010.09993.x
Jahan, N., Chowdhury, A., Li, T., Xu, K., Wei, F., and Wang, S. (2021). Neferine improves oxidative stress and apoptosis in benign prostate hyperplasia via Nrf2-ARE pathway. Redox Rep. 26 (1), 1–9. doi:10.1080/13510002.2021.1871814
Kang, D., Chokkalingam, A. P., Gridley, G., Nyren, O., Johansson, J. E., Adami, H. O., et al. (2007). Benign prostatic hyperplasia and subsequent risk of bladder cancer. Br. J. Cancer 96 (9), 1475–1479. doi:10.1038/sj.bjc.6603730
Larsson, S. C., Michaëlsson, K., and Burgess, S. (2019). Mendelian randomization in the bone field. Bone 126, 51–58. doi:10.1016/j.bone.2018.10.011
Lenis, A. T., Lec, P. M., Chamie, K., and Mshs, M. D. (2020). Bladder cancer: A review. JAMA 324 (19), 1980–1991. doi:10.1001/jama.2020.17598
Lin, Y. H., Hou, C. P., Chen, T. H., Juang, H. H., Chang, P. L., Yang, P. S., et al. (2018). Transurethral resection of the prostate provides more favorable clinical outcomes compared with conservative medical treatment in patients with urinary retention caused by benign prostatic obstruction. BMC Geriatr. 18 (1), 15. doi:10.1186/s12877-018-0709-3
Lv, L., Ye, D., Chen, J., Qian, Y., Fu, A. N., Song, J., et al. (2022). Circulating phosphorus concentration and risk of prostate cancer: A mendelian randomization study. Am. J. Clin. Nutr. 115 (2), 534–543. doi:10.1093/ajcn/nqab342
Ma, C., Lu, B., and Sun, E. (2016). Giant bladder stone in a male patient: A case report. Med. Baltim. 95 (30), e4323. doi:10.1097/MD.0000000000004323
Madersbacher, S., Sampson, N., and Culig, Z. (2019). Pathophysiology of benign prostatic hyperplasia and benign prostatic enlargement: A mini-review. Gerontology 65 (5), 458–464. doi:10.1159/000496289
Mao, W., Wang, K., Sun, S., Wu, J., Chen, M., Geng, J., et al. (2021). ID2 inhibits bladder cancer progression and metastasis via PI3K/AKT signaling pathway. Front. Cell Dev. Biol. 9, 738364. doi:10.3389/fcell.2021.738364
Martinez-Rojo, E., Berumen, L. C., Garcia-Alcocer, G., and Escobar-Cabrera, J. (2021). The role of androgens and androgen receptor in human bladder cancer. Biomolecules 11 (4), 594. doi:10.3390/biom11040594
Matsumoto, S., Shimizu, N., Hanai, T., Uemura, H., and Levin, R. (2009). Bladder outlet obstruction accelerates bladder carcinogenesis. BJU Int. 103 (10), 1436–1439. doi:10.1111/j.1464-410X.2008.08261.x
Melicow, M. M., Uson, A. C., and Lipton, R. (1961). beta-Glucuronidase activity in the urine of patients with bladder cancer and other conditions. J. Urol. 86, 89–94. doi:10.1016/s0022-5347(17)65114-0
Muntzing, J., and Nilsson, T. (1972). Enzyme activity and distribution in the hyperplastic and cancerous human prostate. Scand. J. Urol. Nephrol. 6 (2), 107–111. doi:10.3109/00365597209133620
Naderi, N., Mochtar, C. A., and de la Rosette, J. J. (2004). Real life practice in the management of benign prostatic hyperplasia. Curr. Opin. Urol. 14 (1), 41–44. doi:10.1097/00042307-200401000-00009
Negri, E., Pelucchi, C., Talamini, R., Montella, M., Gallus, S., Bosetti, C., et al. (2005). Family history of cancer and the risk of prostate cancer and benign prostatic hyperplasia. Int. J. Cancer 114 (4), 648–652. doi:10.1002/ijc.20755
Noordam, R., Smit, R. A., Postmus, I., Trompet, S., and van Heemst, D. (2016). Assessment of causality between serum gamma-glutamyltransferase and type 2 diabetes mellitus using publicly available data: A mendelian randomization study. Int. J. Epidemiol. 45 (6), 1953–1960. doi:10.1093/ije/dyw306
Odabasi, Z., and Mert, A. (2020). Candida urinary tract infections in adults. World J. Urol. 38 (11), 2699–2707. doi:10.1007/s00345-019-02991-5
Pierce, B. L., and Burgess, S. (2013). Efficient design for mendelian randomization studies: Subsample and 2-sample instrumental variable estimators. Am. J. Epidemiol. 178 (7), 1177–1184. doi:10.1093/aje/kwt084
Rassweiler, J., Teber, D., Kuntz, R., and Hofmann, R. (2006). Complications of transurethral resection of the prostate (TURP)-incidence, management, and prevention. Eur. Urol. 50 (5), 969–979. doi:10.1016/j.eururo.2005.12.042
Richters, A., Aben, K. K. H., and Kiemeney, L. (2020). The global burden of urinary bladder cancer: An update. World J. Urol. 38 (8), 1895–1904. doi:10.1007/s00345-019-02984-4
Roehrborn, C. G. (2008). BPH progression: Concept and key learning from MTOPS, ALTESS, COMBAT, and ALF-ONE. BJU Int. 101 , 17–21. doi:10.1111/j.1464-410X.2008.07497.x
Selius, B. A., and Subedi, R. (2008). Urinary retention in adults: Diagnosis and initial management. Am. Fam. Physician 77 (5), 643–650.
Sengupta, S., Amling, C., D'Amico, A. V., and Blute, M. L. (2008). Prostate specific antigen kinetics in the management of prostate cancer. J. Urol. 179 (3), 821–826. doi:10.1016/j.juro.2007.10.023
Shiota, M., Kiyoshima, K., Yokomizo, A., Takeuchi, A., Kashiwagi, E., Dejima, T., et al. (2017). Suppressed recurrent bladder cancer after androgen suppression with androgen deprivation therapy or 5α-reductase inhibitor. J. Urol. 197 (2), 308–313. doi:10.1016/j.juro.2016.08.006
Siegel, R. L., Miller, K. D., and Jemal, A. (2019). Cancer statistics. Ca. Cancer J. Clin. 69 (1), 7–34. doi:10.3322/caac.21551
Slatkin, M. (2008). Linkage disequilibrium-understanding the evolutionary past and mapping the medical future. Nat. Rev. Genet. 9 (6), 477–485. doi:10.1038/nrg2361
Sottnik, J. L., Vanderlinden, L., Joshi, M., Chauca-Diaz, A., Owens, C., Hansel, D. E., et al. (2021). Androgen receptor regulates CD44 expression in bladder cancer. Cancer Res. 81 (11), 2833–2846. doi:10.1158/0008-5472.CAN-20-3095
Staley, J. R., Blackshaw, J., Kamat, M. A., Ellis, S., Surendran, P., Sun, B. B., et al. (2016). PhenoScanner: A database of human genotype-phenotype associations. Bioinformatics 32 (20), 3207–3209. doi:10.1093/bioinformatics/btw373
Swerdlow, D. I., Kuchenbaecker, K. B., Shah, S., Sofat, R., Holmes, M. V., White, J., et al. (2016). Selecting instruments for Mendelian randomization in the wake of genome-wide association studies. Int. J. Epidemiol. 45 (5), 1600–1616. doi:10.1093/ije/dyw088
Thiese, M. S. (2014). Observational and interventional study design types; an overview. Biochem. Med. 24 (2), 199–210. doi:10.11613/BM.2014.022
Tong, Y., and Zhou, R. Y. (2020). Review of the roles and interaction of androgen and inflammation in benign prostatic hyperplasia. Mediat. Inflamm. 2020, 7958316–7958317. doi:10.1155/2020/7958316
Tseng, C. H. (2013). Benign prostatic hyperplasia is a significant risk factor for bladder cancer in diabetic patients: A population-based cohort study using the national Health insurance in taiwan. BMC Cancer 13, 7. doi:10.1186/1471-2407-13-7
Verbanck, M., Chen, C. Y., Neale, B., and Do, R. (2018). Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet. 50 (5), 693–698. doi:10.1038/s41588-018-0099-7
Walker, N. F., Gan, C., Olsburgh, J., and Khan, M. S. (2014). Diagnosis and management of intradiverticular bladder tumours. Nat. Rev. Urol. 11 (7), 383–390. doi:10.1038/nrurol.2014.131
Wang, C. S., Li, C. C., Juan, Y. S., Wu, W. J., and Lee, H. Y. (2020). 5α-reductase inhibitors impact prognosis of urothelial carcinoma. BMC Cancer 20 (1), 872. doi:10.1186/s12885-020-07373-4
Wang, D., Li, M., Li, J., Wan, X., Huang, Y., Wang, C., et al. (2020). Comprehensive characterization of androgen-responsive lncRNAs mediated regulatory network in hormone-related cancers. Dis. Markers 2020, 8884450. doi:10.1155/2020/8884450
Wang, M., Jian, Z., Gao, X., Yuan, C., Jin, X., Li, H., et al. (2021). Causal associations between educational attainment and 14 urological and reproductive Health outcomes: A mendelian randomization study. Front. Public Health 9, 742952. doi:10.3389/fpubh.2021.742952
Wang, N., Lu, L., Cao, Q. F., Qian, S., Ding, J., Wang, C., et al. (2021). Partial inhibition of activin receptor-like kinase 4 alleviates bladder fibrosis caused by bladder outlet obstruction. Exp. Cell Res. 406 (1), 112724. doi:10.1016/j.yexcr.2021.112724
Wang, T., Cheng, J., and Wang, Y. (2022). Genetic support of a causal relationship between iron status and atrial fibrillation: A mendelian randomization study. Genes Nutr. 17 (1), 8. doi:10.1186/s12263-022-00708-9
Wang, Y., Zhang, Z., Shao, J., Lu, Y., Li, X., and Li, R. (2021). Outcome of simultaneous thulium laser enucleation of bladder tumor and prostate in patients with non-muscle invasive bladder tumor and benign prostatic hyperplasia: A matched-pair comparison with a long-term follow-up. World J. Urol. 39 (1), 105–111. doi:10.1007/s00345-020-03158-3
Wu, F., Huang, Y., Hu, J., and Shao, Z. (2020). Mendelian randomization study of inflammatory bowel disease and bone mineral density. BMC Med. 18 (1), 312. doi:10.1186/s12916-020-01778-5
Wu, S. C., Kwon, D., Jue, J. S., Chen, F. V., Velasquez Escobar, M. C., Punnen, S., et al. (2021). Androgen suppression therapy is associated with lower recurrence of non-muscle-invasive bladder cancer. Eur. Urol. Focus 7 (1), 142–147. doi:10.1016/j.euf.2019.04.021
Xiong, Y., Ju, L., Yuan, L., Chen, L., Wang, G., Xu, H., et al. (2021). KNSTRN promotes tumorigenesis and gemcitabine resistance by activating AKT in bladder cancer. Oncogene 40 (9), 1595–1608. doi:10.1038/s41388-020-01634-z
Yang, C., Yuan, W., Yang, X., Li, P., Wang, J., Han, J., et al. (2018). Circular RNA circ-ITCH inhibits bladder cancer progression by sponging miR-17/miR-224 and regulating p21, PTEN expression. Mol. Cancer 17 (1), 19. doi:10.1186/s12943-018-0771-7
Zeng, Y., Kakehi, Y., Nouh, M. A., Tsunemori, H., Sugimoto, M., and Wu, X. X. (2009). Gene expression profiles of lysophosphatidic acid-related molecules in the prostate: Relevance to prostate cancer and benign hyperplasia. Prostate 69 (3), 283–292. doi:10.1002/pros.20879
Zhao, Y., Fan, Y., Wang, M., Yu, C., Zhou, M., Jiang, D., et al. (2021). Kidney stone disease and cardiovascular events: A study on bidirectional causality based on mendelian randomization. Transl. Androl. Urol. 10 (12), 4344–4352. doi:10.21037/tau-21-899
Zhou, J., Kelsey, K. T., Giovannucci, E., and Michaud, D. S. (2014). Fluid intake and risk of bladder cancer in the Nurses' Health Studies. Int. J. Cancer 135 (5), 1229–1237. doi:10.1002/ijc.28764
Zhou, L., Liang, X., and Zhang, K. (2020). Assessment of the clinical efficacy of simultaneous transurethral resection of both bladder cancer and the prostate: A systematic review and meta-analysis. Aging Male 23 (5), 1182–1193. doi:10.1080/13685538.2020.1718637
Keywords: mendelian randomization, genome-wide association study, benign prostate hyperplasia, bladder cancer, causal inference
Citation: Du W, Wang T, Zhang W, Xiao Y and Wang X (2022) Genetically supported causality between benign prostate hyperplasia and urinary bladder neoplasms: A mendelian randomization study. Front. Genet. 13:1016696. doi: 10.3389/fgene.2022.1016696
Received: 11 August 2022; Accepted: 31 October 2022;
Published: 17 November 2022.
Edited by:
Giuseppe Morgia, University of Catania, ItalyReviewed by:
Jacklyn Hellwege, Vanderbilt University Medical Center, United StatesZeyan Li, Shandong University, China
Copyright © 2022 Du, Wang, Zhang, Xiao and Wang. 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: Yu Xiao, eXUueGlhb0B3aHUuZWR1LmNu; Xinghuan Wang, d2FuZ3hpbmdodWFuQHdodS5lZHUuY24=