- 1Department of Reproductive Medicine, Zhongda Hospital Affiliated to Southeast University, Nanjing, China
- 2School of Medicine, Southeast University, Nanjing, China
Background: Recent researches have shown a correlation between the gut microbiota (GM) and various diseases. However, it remains uncertain whether the relationship between GM and benign prostatic hyperplasia (BPH) is causal.
Methods: We carried out a two-sample Mendelian randomization (MR) analysis, utilizing data from the most extensive GM-focused genome-wide association study by the MiBioGen consortium, with a sample size of 13,266. Data for BPH, encompassing 26,358 cases and 110,070 controls, were obtained from the R8 release of the FinnGen consortium. We employed multiple techniques, such as inverse variance weighted (IVW), constrained maximum likelihood and model averaging methods, maximum likelihood, MR-Pleiotropy RESidual Sum and Outlier (MRPRESSO),MR-Egger, and weighted median methods, to investigate the causal relationship between GM and BPH. To evaluate the heterogeneity among the instrumental variables, Cochran’s Q statistics were employed. Additionally, the presence of horizontal pleiotropy was assessed through the application of both MR-Egger and MR-PRESSO tests. The direction of causality was scrutinized for robustness using the MR-Steiger directionality test. A reverse MR analysis examined the GM previously linked to BPH through a causal relationship in the forward MR assessment.
Results: According to the analysis conducted using IVW,Eisenbergiella (odds ratio [OR]=0.92, 95% confidence interval [CI]: 0.85–0.99,P=0.022) and Ruminococcaceae (UCG009) (OR=0.88, 95% CI: 0.79–0.99, P=0.027) were found to reduce the risk of BPH, while Escherichia shigella (OR=1.19, 95% CI: 1.05–1.36, P=0.0082) appeared to increase it. The subsequent reverse MR analysis revealed that the three GM were not significantly influenced by BPH, and there was no noticeable heterogeneity or horizontal pleiotropy among the instrumental variables.Conclusion: These results indicated a causal relationship between Eisenbergiella, Ruminococcaceae (UCG009), and Escherichia shigella and BPH. Further randomized controlled trials are needed to explore more comprehensively the roles and operational mechanisms of these GM in relation to BPH.
1 Introduction
Benign prostatic hyperplasia (BPH) is a common disease affecting middle-aged and elderly men; the clinical symptoms include frequent urination, urgent urination, and increased nocturia (Berry et al., 1984). The incidence of BPH gradually increases with age. The clinical prevalence of BPH in men aged >50 years is 50–75%, which increases with age. The prevalence of BPH in men aged > 70 years is >80% (Wang et al., 2020). The disease course seriously affects the quality of life of middle-aged and elderly men. Many studies have been conducted on the pathogenesis of BPH; however, its etiology has not yet been clarified. Currently, the accepted risk factors for the pathogenesis of BPH include functional testis and aging (Braeckman and Denis, 2017). Additionally, there are other theories, such as the role of androgen and its receptors (Carson and Rittmaster, 2003), imbalance between cell proliferation and apoptosis (Rho et al., 2019; Yuan et al., 2020), role of growth factor neurotransmitters (La Vignera et al., 2016), prostate interstitial glandular epithelial interaction (Timme et al., 1995), and inflammatory factors (Kohnen and Drach, 1979). Smoking, obesity, alcoholism, family history, race, and geographical environment are also associated with BPH occurrence (Bushman, 2009; Vuichoud and Loughlin, 2015). In recent years, many studies have shown that changes in the gut microbiota (GM) are of great significance for the diagnosis, prevention, and treatment of BPH (An et al., 2023; Russo et al., 2023).
The GM comprises a genome of microbiota that inhabits the gastrointestinal tract and works symbiotically with the host to maintain health (Connelly et al., 2019). According to recent publications, the human GM contains trillions of microorganisms, including bacteria, fungi, protists, archaea, and viruses (Yatsunenko et al., 2012). The human GM is viewed across age and geography. The composition of each individual’s GM is influenced by several factors, including age, diet, lifestyle, and drug intake (Yatsunenko et al., 2012; Goodrich et al., 2014). The GM has become a research hotspot in recent years. There have been many reports on the relationship between GM and various diseases, confirming that GM plays an important role in the occurrence and progression of diseases. The relationship between GM dysbiosis and several diseases has been established; some examples include periodontal disease, cancer, obesity, diabetes, and chronic fatigue syndrome (Ditto et al., 2021). However, a growing number of studies have shown that differences between healthy and bio-imbalanced microbiota may be associated with prostate diseases, such as chronic prostatitis, BPH, prostate cancer, and other prostate diseases (Fernandes et al., 2022). Tsai et al. observed significant differences in the microbial composition of urine samples between patients with BPH and healthy controls (Tsai et al., 2022). Li et al. used intestinal samples from BPH rats and healthy control rats, and the GM composition was analyzed using 16S rDNA sequencing and liquid chromatography-tandem mass spectrometry. BPH may affect the GM and gut metabolites in a variety of ways, including changing the proportion of GM through the nervous and psychological systems. This regulates the body’s metabolism and hormone synthesis and mediates the inflammatory response of the host and immune system to regulate the changes in GM. Consequently, GM and gut metabolites are affected (Li et al., 2022). BPH can cause complications such as prostatitis, renal damage, and metabolic syndrome, which are associated with GM (Tang et al., 2015; Bishehsari et al., 2020; Liu et al., 2021). In the present study, we explored the relationship between GM and BPH.
Currently, a few high-quality randomized controlled trials (RCTs) have investigated the relationship between the GM and BPH. Furthermore, observational studies cannot completely eliminate the influence of confounding factors and are susceptible to recall bias and reverse causation. Mendelian randomization (MR) is a contemporary method for examining the correlation between GM and BPH. MR employs instrumental variables (IVs), in which genetic variants serve as instruments. During pregnancy, alleles are randomly assigned, a process akin to the randomization of treatment and control groups in RCTs (Li et al., 2023), and common confounding factors have no impact on the correlation between genetic variations and outcomes (Davey and Hemani, 2014). To date, MR has been extensively employed for assessing potential causal links between exposure and disease, as mentioned in reference (Shu et al., 2019). In this research, we conducted a two-sample MR analysis to evaluate the potential connection between GM and BPH. Summary statistics from genome-wide association study (GWAS) conducted by the MiBioGen and FinnGen consortiums were utilized for this purpose.
2 Methods
2.1 Data source
In the MiBioGen consortium study, 18,340 participants underwent analysis for host genotypes and 16S fecal microbiome rRNA gene sequencing profiles. The cohorts included participants from a wide array of countries: USA, Canada, Israel, South Korea, Germany, Denmark, the Netherlands, Belgium, Sweden, Finland, and the UK. Among these cohorts, 20 consisted of individuals from single ancestries, namely European (16 cohorts, N=13,266), Middle-Eastern (1 cohort, N=481), East Asian (1 cohort, N=811), American Hispanic/Latin (1 cohort, N=1,097), and African American (1 cohort, N=114). Additionally, four cohorts had mixed ancestry, totaling 2,571 participants. More comprehensive details on these 24 cohorts are available in a previously published study (Kurilshikov et al., 2021). The GWAS identified 211 GM taxa (from genus to phylum levels) and 122,110 associated single nucleotide polymorphisms (SNPs) (Kurilshikov et al., 2021). Genetic variants correlated with 9 phyla, 16 classes, 20 orders, 35 families, and 131 genera have been identified (van der Velde et al., 2019; Zeng et al., 2023). In the present study, the lowest taxonomic classification was at the genus level. The findings indicated that 131 genera were identified with a mean abundance exceeding 1%, of which 12 were unidentifiable (Kurilshikov et al., 2021). Accordingly, 119 genera were incorporated in the current investigation for analytical purposes. Summary statistics for the BPH GWAS were acquired through the R8 release dataset provided by the FinnGen consortium (Dai et al., 2021), encompassing 26,358 cases and 110,070 controls. The FinnGen research initiative harmonizes genomic information in Finnish biobanks with health-related data from the country’s healthcare databases, research endpoints in this research were defined using ICD codes (Kurki et al., 2023). In the analysis, adjustments were made for factors such as gender, age, genotyping batch, and the initial ten principal elements. Comprehensive details about the utilized GWAS are provided in Table 1.
2.2 Instrumental variable
To ensure accuracy of the outcome regarding the causal influence of GM on BPH, a series of quality control procedures were implemented for the selection of genetic predictors related to microbiome characteristics. First, the IVs chosen for the analysis must demonstrate a significant correlation with the exposure factors. In this study, SNPs with p-values below the locus-wide significance threshold (1 × 10-5) were selected for adequate IVs. Second, only SNPs with the lowest p-values were included in the calculation of linkage disequilibrium (LD) among SNPs with R2 <0.1 (clumping window size=500 kb) (Ni et al., 2021) to ensure the independence of all IVs. This approach minimized the effects of LD, which violates randomized allele allocation. Third, allele frequency data was employed to deduce forward-strand alleles when dealing with palindromic SNPs. Finally, after harmonizing the effects of exposure and outcomes, we employed the F-statistic of SNPs to ascertain the strength and stability of the link between IVs and exposure factors. IVs with F values ≤10, indicating bias toward weak IVs, were excluded. The F-statistic was calculated as F=β2 exposure/SE2 exposure.
2.3 Statistical analysis
The basic principle of MR is based on the Mendelian law of inheritance. During gamete formation, offspring alleles are randomly assigned by the parents. Before an individual is born, genes and the phenotype determined by the genes are randomly assigned to the population. MR study can cleverly avoid the influence of confounding factors brought by postnatal environmental changes and interference of reverse causality; therefore, it is referred to as “nature’s randomized trial” (Hingorani and Humphries, 2005).
In this research, the MR approach was employed to investigate the potential causal link between GM and BPH, while rigorously minimizing the influence of confounding factors. We performed a two-sample MR analysis, synthesizing data from GM and BPH GWAS, to further scrutinize this hypothesized causal association. Various MR methods were uesd in our analyses, comprising inverse variance weighting (IVW), constrained maximum likelihood and model averaging (cML-MA), maximum likelihood (ML), MR-Egger regression, weighted median, MR Pleiotropy Residual Sum and Outlier (MR-PRESSO), and MR Steiger. The IVW method is a weighted linear regression model proposed by Burgess et al. in 2013, which is the most commonly used calculation method and most important method in MR studies (Burgess et al., 2013). A consistent estimate of the causal effect of exposure on the outcome was obtained by combining the Wald ratio of the causal effect of each SNP. The regression did not consider the existence of an intercept term. The ML method is a conventional approach that exhibits a low standard error, similar to that of the IVW method. In the absence of horizontal pleiotropy or heterogeneity, the outcomes will remain unbiased, with standard errors smaller compared to IVW (Pierce and Burgess, 2013). The largest difference between the MR-Egger regression and IVW analysis methods is the addition of intercept term, which is mainly used to determine whether there is horizontal pleiotropy (Slob et al., 2017). Even if all genetic variants are invalid IVs, the MR-Egger test yields consistent causal effect estimates, as long as the association of each genetic variant with exposure is independent of the pleiotropic effect of the variant (rather than through exposure) (Bowden et al., 2015). When only half of the IVs are valid, applying the weighted median method can provide robust estimates (Bowden et al., 2016). The cML-MA method is a MR technique rooted in model averaging that does not rely on the InSIDE assumptions when managing pleiotropic impacts, whether they are linked or independent (Xue et al., 2021).
The Cochran’s Q statistic served as the measurement tool for the heterogeneity among IVs. The MR-PRESSO is commonly used for verifying horizontal pleiotropy. It consists of three parts: a global, outlier, and distortion tests. The MR-PRESSO method can calculate whether there is horizontal pleiotropy and screen for abnormal SNPs. If abnormal SNPs are screened, they must be deleted before calculating the MR results to ensure their stability (Verbanck et al., 2018). In this study, SNPs exhibiting pleiotropic outliers (P<0.05) were excluded using MR-PRESSO. Furthermore, a “leave-one-out” analysis was conducted to visually identify potentially heterogeneous SNPs by skipping each SNP in turn. The MR Steiger directionality test provided an assessment of the robustness in determining directionality. In order to establish the connection between GM and BPH, a reverse MR analysis was performed on the bacterial strains that were significantly linked with BPH in the forward MR analysis. The settings and procedures corresponded to those of the forward MR.
The statistical power of the MR estimates was determined using an online calculator tool provided by Stephen Burgess. Computation of MR estimates requires the inclusion of R2 summations for each SNP (Burgess, 2014).
R2 denotes the fraction of variability in exposure that can be accounted for by genetic variation:
It was postulated that a significant causal association between GM and BPH existed and is subject to the fulfilment of the following criteria: First, the IVW method revealed a notable distinction (P<0.05); Second, the five approaches were in harmony regarding the direction of their estimates; Third, both the MR-Egger intercept and the MR-PRESSO global test did not achieve statistical relevance (P>0.05); Fourth, the MR Steiger directionality tests showed TRUE. The analysis was carried out using the TwoSampleMR (version 0.5.6) (Hemani et al., 2018), MRcML (version 0.0.0.9000) (Xue et al., 2021), and MR-PRESSO (version 1.0) (Verbanck et al., 2018) packages, all within the R software environment (version 4.2.2).
3 Results
3.1 Instrumental variable selection
A total of 1,269 SNPs were chosen as IVs for 119 GM, based on predefined criteria. The analysis revealed F-statistics greater than 10 for these SNPs, indicating their robustness as IVs. Consequently, the study’s outcomes show no signs of weak instrument bias, affirming the reliability of the results. Detailed information is available in Supplementary File 1: Table S1.
3.2 Two-sample MR analysis
After applying the standard, three GM genera—Eisenbergiella, Escherichia shigella, and Ruminococcaceae (UCG009)—were identified as having a connection to BPH. The IVW method highlighted this relationship by showing a marked difference (P<0.05), with a uniform direction observed in estimates from all five methods.
The IVW analysis identified that Eisenbergiella (odds ratio [OR]=0.92, 95% confidence interval [CI]:0.85–0.99, P=0.022) and Ruminococcaceae (UCG009) (OR=0.88, 95% CI: 0.79–0.99, P=0.027) appeared to confer a protective influence against BPH. Conversely, the IVW analysis also showed that Escherichia shigella (OR=1.19, 95% CI: 1.05–1.36, P=0.0082) was implicated in elevating BPH risk (Table 2). The Cochran’s Q test demonstrated uniformity across the IVs, with no detected heterogeneity. This is further validated by MR-Egger and MR-PRESSO analyses, which largely ruled out horizontal pleiotropic effects in these IVs. In scrutinizing the data plots (Figures 1, 2), no significant deviations or anomalies were apparent, and the MR-PRESSO analysis did not identify any SNPs as outliers. The MR Steiger directionality tests further reinforced the robust connection from the GM to BPH for all evaluated results (Table 3). Comprehensive statistics for the influence of the 119 GM on BPH can be found in Supplementary Tables S2–7.
Figure 1 Scatter plots of significant causality between GM and BPH. Scatter plot of the effect size and 95% CI of each SNP on GM and BPH risk. The horizontal axis reflects genetic effect of each SNP on GM. The vertical axis represents the genetic effect of each SNP on BPH risk. MR, Mendelian randomization; BPH, Benign prostatic hyperplasia; GM, gut microbiota; SNP, Single nucleotide polymorphism; CI, confidence interval.
Figure 2 Leave-one-out analysis for the impact of individual SNPs on the association between GM and BPH risk. By leaving out exactly one SNP, it shows how each individual SNP influences the overall estimate. SNPs, Single nucleotide polymorphisms; GM, gut microbiota; BPH, Benign prostatic hyperplasia.
The reverse MR analysis provided no definitive proof of a causal link between BPH and the three GM, as detailed in Table 3 and Supplementary File 1: Table S9. The Cochran’s Q test, MR-Egger, and MR-PRESSO analyses (Table 3: Supplementary File 1: Table S10–12) revealed no significant heterogeneity or horizontal pleiotropy. The MR Steiger directionality tests confirmed the relationship between BPH and the three GM, as indicated as TRUE (Table 3: Supplementary File 1: Table S13).
4 Discussion
In recent years, the human GM has been extensively studied. GM evolved with the evolution of the host, and the flora obtained at birth developed with the development of the host (Adak and Khan, 2019). GM dynamically colonizes the intestinal mucous layer of the human body. The total number of GM is approximately 10 times that of human cells, and the total number of genomes is approximately 150 times that of human genes. It is often called the “virtual organ” of the human body (Mu et al., 2017) and the “second genome” that affects human health. There is a two-way interaction between the host and GM, which is affected by the diversity and abundance of GM, and it has a certain effect on the health of the host body (Dominguez-Bello et al., 2010). It plays important roles in the regulation of digestion, nutrient absorption, energy metabolism, fat metabolism, and immunity (Farre et al., 2020). Studies have shown that an imbalance in GM is related to cardiovascular (Wang and Zhao, 2018), neurological (Socala et al., 2021), respiratory (Chunxi et al., 2020) and metabolic (Fan and Pedersen, 2021) diseases, malignant tumors (Long et al., 2023), and other diseases (Gomaa, 2020).
BPH may affect GM and gut metabolites in many ways. BPH can cause complications such as prostatitis, renal damage, and metabolic syndrome, which have been found to be associated with GM (Li et al., 2022). To date, the effect of GM on the occurrence and development of BPH is unclear, and there are few related studies. Jinho et al. found that changes in the GM induced by BPH and finasteride treatment were closely related to the regulation of prostate morphology, hormones, and prostate cell apoptosis. Therefore, intestinal microorganisms can be used as indicators and therapeutic measures (An et al., 2023).
In our study, a two-sample MR analysis was conducted to investigate the potential causal relationship between GM and BPH. This study marks an initial extensive MR study into the causal association between GM and BPH at the gene prediction level, employing a significant volume of GWAS data. The findings of this study may provide valuable insights for future management of BPH by highlighting the relevance of GM to the condition. Specifically, this study found that Eisenbergiella and Ruminococcaceae (UCG009) exhibited protective effects against BPH, whereas Escherichia shigella showed the opposite effects.
Eisenbergiella is a proposed genus of Spirillaceae, which includes Eisenbergiella tayi, an anaerobic bacterium. Butyric acid is the primary fermentation end product and energy source of colonic cells (Jousimies-Somer, 2002) and is involved in the maintenance of colonic mucosal health (Amir et al., 2014). Although the biological function information of this genus is very limited, Eisenbergiella belongs to Lachnospiraceae. Studies have shown that an increase in Lachnospiraceae may be related to the occurrence of metabolic diseases, such as diabetes (Kameyama and Itoh, 2014). However, few studies have examined the relationship between Eisenbergiella and BPH. Escherichia shigella is a conditional pathogen. Its DNA is closely related to that of Escherichia coli and can cause acute gastrointestinal infections. It is one of the most important pathogenic bacteria in the human body and is balanced with other bacteria in the intestinal tract. Escherichia shigella has two important roles in the human body: first, it plays an important role in inducing intestinal inflammation; second, in an inflammatory environment, it can gain a stronger survival advantage by secreting enterosin and aggravating the disorder of intestinal microecology (Singh et al., 2015). Lee et al. recruited 77 BPH patients and 30 controls who had not recently taken antibiotics. They found that there was a significant difference in the abundance of Escherichia shigella between the BPH and control groups (Lee et al., 2021). Ruminococcaceae (UCG009) is an anti-inflammatory bacterium found in the cecum and colon. As one of the producers of short-chain fatty acids, it is responsible for the degradation of various indigestible polysaccharides and fibers (Hooda et al., 2012; Donaldson et al., 2016). However, there are limited studies on the effects of Ruminococcaceae (UCG009) on BPH.
This study had several strengths. In this research, MR analysis was conducted to explore the causal link between GM and BPH in humans, eliminating the influence of confounding variables. The GM summary statistics were sourced from the most extensive GWAS meta-analysis conducted by the MiBioGen consortium, while the BPH summary statistics came from the FinnGen consortium’s R8 release data. This provided robustness to the instruments and incorporated the most up-to-date information in the MR analysis. The MR-PRESSO and MR-Egger regression intercept-term tests were employed to identify and rule out horizontal pleiotropy. Furthermore, cML-MA was used to exclude bias from both linked and independent pleiotropies. Non-overlapping exposure and result summary-level data were employed in this MR study to prevent bias.
Nonetheless, the interpretation of the study’s outcomes must take several limitations into account. First, we could not study the non-linear correlations because the analysis was conducted on aggregated numbers instead of raw data. Second, the GWAS analysis of GM did not confined only to male subjects. Although this analysis considers sex differences and ignores any genetic variation in male and female chromosomes, it cannot avoid any potential bias caused by sex (Verbanck et al., 2018). Third, our study was limited to European patients; therefore, the exclusive focus on European patients raises questions regarding the applicability of the findings to non-European populations, due to inherent genetic distinctions among ethnic groups. Finally, because only generic-level exposure data were available, it was impossible to further explore the connection between GM and BPH at the species level.
5 Conclusion
The outcomes of this two-sample MR study indicate that Eisenbergiella and Ruminococcaceae (UCG009) may confer protective effects against BPH, whereas Escherichia shigella may have a positive relationship with the risk of BPH. Further research is needed to obtain a more comprehensive understanding of the possible beneficial or harmful impacts of these GM on BPH and the mechanisms behind them. Although there is a dearth of evidence demonstrating the impact of BPH on GM, it is possible that an effect may exist, this possibility calls for more research to ascertain the truth.
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
Ethical approval was not provided for this study on human participants because this research has been conducted using published studies and consortia providing publicly available summary statistics. All original studies have been approved by the corresponding ethical review board. In addition, no individual-level data was used in this study. Therefore, no new ethical review board approval was required.
Author contributions
Conception and design: DX, JW, and XZ. Data collection and assembly: DX, JW, XZ, TS, and LL. Data analysis and interpretation: DX, JW, XZ, TS, and LL. Manuscript writing: DX. Final approval of the manuscript: All authors. DX and JW are the co-first authors who made equal contributions to this study. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by grants from the Self-Select Fund of Zhongda Hospital [ZDYYZXKT2019001].
Acknowledgments
We would like to thank all contributors who shared the data from this study. We would like to thank all participants and investigators in the FinnGen study, all participants and investigators in the gut microbiota study, and investigators and participants in the original GWASs.
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/fcimb.2023.1248381/full#supplementary-material
Abbreviations
BPH, benign prostatic hyperplasia; CI, confidence interval; cML-MA, constrained maximum likelihood and model averaging; GM, gut microbiota; GWAS, genome-wide association study; IV, Instrumental Variable; IVW, inverse variance weighted; LD, linkage disequilibrium; MR, Mendelian randomization; MR-PRESSO, MR Pleiotropy Residual Sum and Outlier; OR, odds ratio; RCTs, randomized controlled trials; SNPs, Single nucleotide polymorphisms.
References
Adak, A., Khan, M. R. (2019). An insight into gut microbiota and its functionalities. Cell Mol. Life Sci. 76 (3), 473–493. doi: 10.1007/s00018-018-2943-4
Amir, I., Bouvet, P., Legeay, C., Gophna, U., Weinberger, A. (2014). Eisenbergiella tayi gen. nov., sp. nov., isolated from human blood. Int. J. Syst. Evol. Microbiol. 64 (Pt 3), 907–914. doi: 10.1099/ijs.0.057331-0
An, J., Song, Y., Kim, S., Kong, H., Kim, K. (2023). Alteration of gut microbes in benign prostatic hyperplasia model and finasteride treatment model. Int. J. Mol. Sci. 24 (6). doi: 10.3390/ijms24065904
Berry, S. J., Coffey, D. S., Walsh, P. C., Ewing, L. L. (1984). The development of human benign prostatic hyperplasia with age. J. Urol. 132 (3), 474–479. doi: 10.1016/S0022-5347(17)49698-4
Bishehsari, F., Voigt, R. M., Keshavarzian, A. (2020). Circadian rhythms and the gut microbiota: from the metabolic syndrome to cancer. Nat. Rev. Endocrinol. 16 (12), 731–739. doi: 10.1038/s41574-020-00427-4
Bowden, J., Davey, S. G., Burgess, S. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int. J. Epidemiol. 44 (2), 512–525. doi: 10.1093/ije/dyv080
Bowden, J., Davey, S. G., Haycock, P. C., 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
Braeckman, J., Denis, L. (2017). Management of BPH then 2000 and now 2016 - From BPH to BPO. Asian J. Urol. 4 (3), 138–147. doi: 10.1016/j.ajur.2017.02.002
Burgess, S. (2014). Sample size and power calculations in Mendelian randomization with a single instrumental variable and a binary outcome. Int. J. Epidemiol. 43 (3), 922–929. doi: 10.1093/ije/dyu005
Burgess, S., Butterworth, A., 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
Bushman, W. (2009). Etiology, epidemiology, and natural history of benign prostatic hyperplasia. Urol Clin. North Am. 36 (4), 403–415. doi: 10.1016/j.ucl.2009.07.003
Carson, C. R., and Rittmaster, R. (2003). The role of dihydrotestosterone in benign prostatic hyperplasia. Urology 61 (4 Suppl 1), 2–7. doi: 10.1016/S0090-4295(03)00045-1
Chunxi, L., Haiyue, L., Yanxia, L., Jianbing, P., Jin, S. (2020). The gut microbiota and respiratory diseases: new evidence. J. Immunol. Res. 2020, 2340670. doi: 10.1155/2020/2340670
Connelly, S., Fanelli, B., Hasan, N. A., Colwell, R. R., Kaleko, M. (2019). Oral beta-lactamase protects the canine gut microbiome from oral amoxicillin-mediated damage. Microorganisms 7 (5). doi: 10.3390/microorganisms7050150
Dai, J., Zhang, T., Guo, J., Zhou, Q., Gu, Y., Zhang, J., et al. (2021). Homozygous pathogenic variants in ACTL9 cause fertilization failure and male infertility in humans and mice. Am. J. Hum. Genet. 108 (3), 469–481. doi: 10.1016/j.ajhg.2021.02.004
Davey, S. G., Hemani, G. (2014). Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum. Mol. Genet. 23 (R1), R89–R98. doi: 10.1093/hmg/ddu328
Ditto, M. C., Parisi, S., Landolfi, G., Borrelli, R., Realmuto, C., Finucci, A., et al. (2021). Intestinal microbiota changes induced by TNF-inhibitors in IBD-related spondyloarthritis. RMD Open 7 (3). doi: 10.1136/rmdopen-2021-001755
Dominguez-Bello, M. G., Costello, E. K., Contreras, M., Magris, M., Hidalgo, G., Fierer, N., et al. (2010). Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proc. Natl. Acad. Sci. U.S.A. 107 (26), 11971–11975. doi: 10.1073/pnas.1002601107
Donaldson, G. P., Lee, S. M., Mazmanian, S. K. (2016). Gut biogeography of the bacterial microbiota. Nat. Rev. Microbiol. 14 (1), 20–32. doi: 10.1038/nrmicro3552
Fan, Y., Pedersen, O. (2021). Gut microbiota in human metabolic health and disease. Nat. Rev. Microbiol. 19 (1), 55–71. doi: 10.1038/s41579-020-0433-9
Farre, R., Fiorani, M., Abdu, R. S., Matteoli, G. (2020). Intestinal permeability, inflammation and the role of nutrients. Nutrients 12 (4). doi: 10.3390/nu12041185
Fernandes, A., Oliveira, A., Guedes, C., Fernandes, R., Soares, R., Barata, P.. (2022). Effect of radium-223 on the gut microbiota of prostate cancer patients: A pilot case series study. Curr. Issues Mol. Biol. 44 (10), 4950–4959. doi: 10.3390/cimb44100336
Gomaa, E. Z. (2020). Human gut microbiota/microbiome in health and diseases: a review. Antonie Van Leeuwenhoek 113 (12), 2019–2040. doi: 10.1007/s10482-020-01474-7
Goodrich, J. K., Waters, J. L., Poole, A. C., Sutter, J. L., Koren, O., Blekhman, R., et al. (2014). Human genetics shape the gut microbiome. Cell 159 (4), 789–799. doi: 10.1016/j.cell.2014.09.053
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. doi: 10.7554/eLife.34408.012
Hingorani, A., Humphries, S. (2005). Nature’s randomised trials. Lancet 366 (9501), 1906–1908. doi: 10.1016/S0140-6736(05)67767-7
Hooda, S., Boler, B. M., Serao, M. C., Brulc, J. M., Staeger, M. A., Boileau, T. W., et al. (2012). 454 pyrosequencing reveals a shift in fecal microbiota of healthy adult men consuming polydextrose or soluble corn fiber. J. Nutr. 142 (7), 1259–1265. doi: 10.3945/jn.112.158766
Jousimies-Somer, H. R. (2002). Wadsworth-KTL Anaerobic Bacteriology Manual (Belmont: Star Publisher).
Kameyama, K., Itoh, K. (2014). Intestinal colonization by a Lachnospiraceae bacterium contributes to the development of diabetes in obese mice. Microbes Environ. 29 (4), 427–430. doi: 10.1264/jsme2.ME14054
Kohnen, P. W., Drach, G. W. (1979). Patterns of inflammation in prostatic hyperplasia: a histologic and bacteriologic study. J. Urol. 121 (6), 755–760. doi: 10.1016/S0022-5347(17)56980-3
Kurilshikov, A., Medina-Gomez, C., Bacigalupe, R., Radjabzadeh, D., Wang, J., Demirkan, A., et al. (2021). Large-scale association analyses identify host factors influencing human gut microbiome composition. Nat. Genet. 53 (2), 156–165. doi: 10.1038/s41588-020-00763-1
Kurki, M. I., Karjalainen, J., Palta, P., Sipila, T. P., Kristiansson, K., Donner, K. M., et al. (2023). FinnGen provides genetic insights from a well-phenotyped isolated population. Nature 613 (7944), 508–518. doi: 10.1038/s41586-022-05473-8
La Vignera, S., Condorelli, R. A., Russo, G. I., Morgia, G., Calogero, A. E. (2016). Endocrine control of benign prostatic hyperplasia. Andrology 4 (3), 404–411. doi: 10.1111/andr.12186
Lee, H. Y., Wang, J. W., Juan, Y. S., Li, C. C., Liu, C. J., Cho, S. Y., et al. (2021). The impact of urine microbiota in patients with lower urinary tract symptoms. Ann. Clin. Microbiol. Antimicrob. 20 (1), 23. doi: 10.1186/s12941-021-00428-9
Li, L. Y., Han, J., Wu, L., Fang, C., Li, W. G., Gu, J. M., et al. (2022). Alterations of gut microbiota diversity, composition and metabonomics in testosterone-induced benign prostatic hyperplasia rats. Mil Med. Res. 9 (1), 12. doi: 10.1186/s40779-022-00373-4
Li, B. H., Yan, S. Y., Luo, L. S., Zeng, X. T., Wang, Y. B., Wang, X. H. (2023). Ten interleukins and risk of prostate cancer. Front. Oncol. 13, 1108633. doi: 10.3389/fonc.2023.1108633
Liu, J., Liu, L., Zhang, G., Peng, X. (2021). Poria cocos polysaccharides attenuate chronic nonbacterial prostatitis by targeting the gut microbiota: Comparative study of Poria cocos polysaccharides and finasteride in treating chronic prostatitis. Int. J. Biol. Macromol 189, 346–355. doi: 10.1016/j.ijbiomac.2021.08.139
Long, Y., Tang, L., Zhou, Y., Zhao, S., Zhu, H. (2023). Causal relationship between gut microbiota and cancers: a two-sample Mendelian randomisation study. BMC Med. 21 (1), 66. doi: 10.1186/s12916-023-02761-6
Mu, Q., Kirby, J., Reilly, C. M., Luo, X. M. (2017). Leaky gut as a danger signal for autoimmune diseases. Front. Immunol. 8, 598. doi: 10.3389/fimmu.2017.00598
Ni, J. J., Xu, Q., Yan, S. S., Han, B. X., Zhang, H., Wei, X. T., et al. (2021). Gut microbiota and psychiatric disorders: A two-sample mendelian randomization study. Front. Microbiol. 12, 737197. doi: 10.3389/fmicb.2021.737197
Pierce, B. L., 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
Rho, J., Seo, C. S., Park, H. S., Wijerathne, C. U., Jeong, H. Y., Moon, O. S., et al. (2019). Ulmus macrocarpa Hance improves benign prostatic hyperplasia by regulating prostatic cell apoptosis. J. Ethnopharmacol 233, 115–122. doi: 10.1016/j.jep.2018.11.042
Russo, G. I., Bongiorno, D., Bonomo, C., Musso, N., Stefani, S., Sokolakis, I., et al. (2023). The relationship between the gut microbiota, benign prostatic hyperplasia, and erectile dysfunction. Int. J. Impot Res. 35 (4), 350–355. doi: 10.1038/s41443-022-00569-1
Shu, X., Wu, L., Khankari, N. K., Shu, X. O., Wang, T. J., Michailidou, K., et al. (2019). Associations of obesity and circulating insulin and glucose with breast cancer risk: a Mendelian randomization analysis. Int. J. Epidemiol. 48 (3), 795–806. doi: 10.1093/ije/dyy201
Singh, V., Yeoh, B. S., Xiao, X., Kumar, M., Bachman, M., Borregaard, N., et al. (2015). Interplay between enterobactin, myeloperoxidase and lipocalin 2 regulates E. coli survival in the inflamed gut. Nat. Commun. 6, 7113. doi: 10.1038/ncomms8113
Slob, E., Groenen, P., Thurik, A. R., Rietveld, C. A. (2017). A note on the use of Egger regression in Mendelian randomization studies. Int. J. Epidemiol. 46 (6), 2094–2097. doi: 10.1093/ije/dyx191
Socala, K., Doboszewska, U., Szopa, A., Serefko, A., Wlodarczyk, M., Zielinska, A., et al. (2021). The role of microbiota-gut-brain axis in neuropsychiatric and neurological disorders. Pharmacol. Res. 172, 105840. doi: 10.1016/j.phrs.2021.105840
Tang, W. H., Wang, Z., Kennedy, D. J., Wu, Y., Buffa, J. A., Agatisa-Boyle, B., et al. (2015). Gut microbiota-dependent trimethylamine N-oxide (TMAO) pathway contributes to both development of renal insufficiency and mortality risk in chronic kidney disease. Circ. Res. 116 (3), 448–455. doi: 10.1161/CIRCRESAHA.116.305360
Timme, T. L., Truong, L. D., Slawin, K. M., Kadmon, D., Park, S. H., Thompson, T. C. (1995). Mesenchymal-epithelial interactions and transforming growth factor-beta 1 expression during normal and abnormal prostatic growth. Microsc Res. Tech 30 (4), 333–341. doi: 10.1002/jemt.1070300408
Tsai, K. Y., Wu, D. C., Wu, W. J., Wang, J. W., Juan, Y. S., Li, C. C., et al. (2022). Exploring the Association between Gut and Urine Microbiota and Prostatic Disease including Benign Prostatic Hyperplasia and Prostate Cancer Using 16S rRNA Sequencing. Biomedicines 10 (11). doi: 10.3390/biomedicines10112676
van der Velde, K. J., Imhann, F., Charbon, B., Pang, C., van Enckevort, D., Slofstra, M., et al. (2019). MOLGENIS research: advanced bioinformatics data software for non-bioinformaticians. Bioinformatics 35 (6), 1076–1078. doi: 10.1093/bioinformatics/bty742
Verbanck, M., Chen, C. Y., Neale, B., 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
Vuichoud, C., Loughlin, K. R. (2015). Benign prostatic hyperplasia: epidemiology, economics and evaluation. Can. J. Urol. 22 Suppl 1, 1–6.
Wang, B., Zhang, S., Sun, C., Du, C., He, G., Wen, J. (2020). Comparison between a transurethral prostate split and transurethral prostate resection for benign prostatic hyperplasia treatment in a small prostate volume: a prospective controlled study. Ann. Transl. Med. 8 (16), 1016. doi: 10.21037/atm-20-5462
Wang, Z., Zhao, Y. (2018). Gut microbiota derived metabolites in cardiovascular health and disease. Protein Cell 9 (5), 416–431. doi: 10.1007/s13238-018-0549-0
Xue, H., Shen, X., Pan, W. (2021). Constrained maximum likelihood-based Mendelian randomization robust to both correlated and uncorrelated pleiotropic effects. Am. J. Hum. Genet. 108 (7), 1251–1269. doi: 10.1016/j.ajhg.2021.05.014
Yatsunenko, T., Rey, F. E., Manary, M. J., Trehan, I., Dominguez-Bello, M. G., Contreras, M., et al. (2012). Human gut microbiome viewed across age and geography. Nature 486 (7402), 222–227. doi: 10.1038/nature11053
Yuan, Y. F., Zhu, W.X., Liu, T., He, J. Q., Zhou, Q., Zhou, X., et al. (2020). Cyclopamine functions as a suppressor of benign prostatic hyperplasia by inhibiting epithelial and stromal cell proliferation via suppression of the Hedgehog signaling pathway. Int. J. Mol. Med. 46 (1), 311–319. doi: 10.3892/ijmm.2020.4569
Keywords: benign prostate hyperplasia, gut microbiota, causal inference, mendelian randomization study, genome wide association study
Citation: Xia D, Wang J, Zhao X, Shen T, Ling L and Liang Y (2023) Association between gut microbiota and benign prostatic hyperplasia: a two-sample mendelian randomization study. Front. Cell. Infect. Microbiol. 13:1248381. doi: 10.3389/fcimb.2023.1248381
Received: 02 July 2023; Accepted: 30 August 2023;
Published: 20 September 2023.
Edited by:
Maria De Lourdes Pereira, University of Aveiro, PortugalCopyright © 2023 Xia, Wang, Zhao, Shen, Ling and Liang. 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: Yuanjiao Liang, WXVhbmppYW8xOTY1QDEyNi5jb20=
†These authors have contributed equally to this work