ORIGINAL RESEARCH article

Front. Microbiol., 05 June 2025

Sec. Microorganisms in Vertebrate Digestive Systems

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1530014

Altered gut microbiota in erectile dysfunction patients: a pilot study

  • 1. Department of Urology, The First Affiliated Hospital of Dalian Medical University, Dalian, Liaoning, China

  • 2. Department of Neurology, The First Affiliated Hospital of Dalian Medical University, Dalian, Liaoning, China

Abstract

Purpose:

With the growing body of research on gut microbiota in recent years, various potential associations between gut microbiota and health or disease have been identified. However, the role of gut microbiota in Erectile dysfunction (ED) remains poorly understood. This study aimed to compare the changes in gut microbiota and metabolic pathways between ED males and healthy control group, contributing to the exploration of ED pathogenesis.

Methods:

Fecal samples were collected from 19 ED patients and 15 healthy controls (aged from 18 to 60 years), with erectile function assessed using the 5-item version of the International Index of Erectile Function (IIEF-5). Macro-genomic sequencing was performed on the NovaSeq PE 150 platform to characterize the gut microbiota distribution among the groups.

Results:

No significant differences in alpha diversity of the gut microbiota were observed between the ED and control groups. Additionally, Principal component analysis (PCA) analysis revealed no notable changes in microbiota composition between the two groups. A comparison of the abundance of key species showed that, in the ED group, species such as Ruminococcus gnavus, Thomasclavelia ramosa, Clostridium sp. AF32-12BH, Clostridium nexile, and Eubacterium siraeum were more abundant, while the abundance of Bacteroides intestinalis was decreased compared to the control group. Furthermore, pathways related to nucleotide and lipid metabolism were found to be highly expressed in the ED group.

Conclusion:

This pilot study found a decrease in the abundance of Bacteroides intestinalis and an increase in the abundance of Ruminococcus gnavus in the ED sample. These microbiota changes may contribute to ED by promoting atherosclerosis and inhibiting the degradation of branched-chain amino acids. In the future, it may be possible to achieve better outcomes for ED patients by precisely regulating the gut microbiota.

1 Introduction

Erectile dysfunction (ED) is a prevalent condition among the male population, with prevalence rates ranging from 37.2 to 48.6% (). It is characterized by the persistent inability to achieve or maintain an erection sufficient for a fulfilling sexual life, leading to significant negative impacts on a patient’s psychological well-being and marital relationships (). In response to the growing emphasis on quality of life, the therapeutic approach to ED has evolved beyond merely improving erectile function to include early diagnosis and the identification of risk factors.

Gut microbiota, a dynamic and diverse community, plays a crucial role in maintaining both the internal and external balance of the human body (; ). With the rapid advancements in sequencing technology, there has been increasing recognition of the gut microbiota’s critical role in human health, including its impact on male sexual function. Studies have established a link between intestinal endotoxins and a decline in gonadal function, affecting testosterone secretion and sperm quality, which can result in reduced fertility and quality of life (). Conversely, probiotics have been shown to enhance the intestinal barrier, reduce metabolic endotoxemia and inflammation, and subsequently improve testosterone levels (). Furthermore, a significant body of evidence highlights the close association between gut microbiota and the development of various metabolic diseases, including diabetes, hypertension, hyperlipidemia, and obesity (; ; ). These metabolic disorders are recognized as major risk factors for the development of erectile dysfunction (), suggesting a potential link between gut microbiota and the onset of ED. In addition, an imbalance in the gut microbiota can lead to the proliferation of potentially pathogenic bacteria, disrupting immune homeostasis and triggering pro-inflammatory responses. Metabolites of gut microbiota, such as lipopolysaccharides from the cell walls of Gram-negative bacteria, are considered endotoxins that cause inflammation (). In dysbiosis, lipopolysaccharides released from the death and lysis of intestinal bacteria bind to Toll-like receptor 4 (TLR4) on intestinal epithelial cells, activating the release of inflammatory cytokines (e.g., interleukins, tumor necrosis factor). The abnormal accumulation of these inflammatory factors leads to a chronic inflammatory state, causing endothelial dysfunction and potentially contributing to ED (; ). However, only a limited number of studies have utilized 16S rRNA gene sequencing to explore the microbial composition in ED patients (; ; ), and there is still a lack of comprehensive understanding of species-level information in the gut microbiota and the associated metabolic pathways.

In light of the above, the present study utilized macro-genomic sequencing to investigate the changes in gut microbiota and metabolic pathways between Chinese ED patients and healthy individuals. The aim was to compare the changes in gut microbiota and metabolic pathways between ED males and healthy control group.

2 Matirials and methods

2.1 Study design

This was a case-control study. A total of 19 male patients diagnosed with ED were recruited from the Department of Reproductive Men’s Medicine at the First Affiliated Hospital of Dalian Medical University between May and July 2023. The diagnostic criteria for ED were based on the 5-item version of the International Index of Erectile Function (IIEF-5). Additionally, 15 healthy men with normal sexual function were recruited as controls. Inclusion criteria for participants were as follows: (1) aged between 18 and 60 years; (2) heterosexual, with a regular sexual partner and a stable sexual life for at least 3 months; and (3) normal external genitalia development with no history of trauma. Participants were excluded if they met any of the following criteria: (1) use of antimicrobials, probiotics, or any gastrointestinal disease treatment in the past month; (2) history of diabetes mellitus, hypertension, or thyroid dysfunction; (3) previous radical prostatectomy, pelvic trauma, or surgery; or (4) a history of severe psychiatric disorders such as anxiety, depression, or other serious mental illnesses.

The study was reviewed and approved by the Ethics Committee of the First Affiliated Hospital of Dalian Medical University (PJ-KS-KY-2023-167), and all participants provided written informed consent.

2.2 Clinical information collection and analysis

Demographic and clinical information was collected through questionnaires and electronic medical records. The data collected included age, body mass index (BMI), International Prostate Symptom Score (IPSS), testosterone, estradiol, follicle-stimulating hormone (FSH), luteinizing hormone (LH), prolactin, glucose, HDL, LDL, triglycerides, total cholesterol, homocysteine, and lipoproteins.

2.3 Fecal specimen collection and pretreatment

Fecal samples were collected from all participants using sterile plastic spoons and stored in sterile plastic tubes. The samples were promptly placed in a –20°C refrigerator and subsequently transferred to a –80°C freezer within 24 h for long-term storage.

2.4 Gene sequencing and data analysis

(1) Genomic DNA was extracted from fecal samples using the Tengen Magnetic Bead Kit. (2) DNA purity and integrity were assessed via 1% agarose gel electrophoresis, and DNA concentration was quantified using the Qubit® dsDNA Assay Kit in the Qubit® 2.0 Fluorometer (Life Technologies, CA, United States). The appropriate amount of sample was placed in a centrifuge tube, diluted with sterile water, and adjusted to an OD value between 1.8 and 2.0. (3) Library construction involved using 1 μg of genomic DNA, which was randomly fragmented into ∼350 bp pieces using a Covaris Ultrasonic Crusher. Following end repair, A-tail addition, sequencing junction addition, purification, PCR amplification, and other steps, the library was completed. After library construction, preliminary quantification was performed using Qubit 2.0, and the library was diluted to 2 ng/μL. The insert size was assessed using the Agilent 2100 system. After confirming the insert size met expectations, Q-PCR was employed to accurately determine the effective concentration of the library (with a required concentration > 3 nM) to ensure high-quality output. Once the library passed quality control, different libraries were pooled according to their effective concentration and target sequencing volume for PE150 sequencing.

2.5 Bioinformatics analysis

(1) The raw data obtained from the NovaSeq sequencing platform were preprocessed using fastp to generate clean data for subsequent analysis. The clean data were then assembled and analyzed using MEGAHIT software (); (2) MetaGeneMark was employed to predict the open reading frames (ORFs) for each sample, and the results were de-redundant using CD-HIT software to generate a non-redundant initial gene catalog (). Bowtie2 was used to align the clean data from each sample to the initial gene catalog, and the final gene catalog (non-redundant) was obtained after filtering. The final gene catalog (unigenes) was used for further analysis (); (3) The unigenes were annotated against Micro_NR sequences and the KEGG database using DIAMOND software (). Alpha diversity was assessed using Shannon and Simpson indices. Principal component analysis (PCA) was performed, and differences between groups were examined using ANOSIM analysis. MetaGenomeSeq and LEfSe analyses were then conducted to identify species that differed between groups (). MetaGenomeSeq analysis was used to perform hypothesis testing between groups for each taxonomic stratum to obtain p-values and Q-values, while LEfSe analysis was performed using LEfSe software (LDA score set to 2 by default) (). Random forest algorithms were applied for regression analyses (R-pROC and randomforest packages) (), and random forest models were constructed by selecting species by gradient at the species level. Important species were screened by MeanDecreaseAccurance, and then each model was cross-validated (default 10-fold) to plot ROC curves.

2.6 Statistical analysis

Statistical analysis was performed using SPSS26 and version 3.6.1 R software. Comparisons between the two groups were made using the independent samples t-test and were expressed as ( ± s). p < 0.05 was considered statistically significant.

3 Results

3.1 Clinical characteristics of the participants

A total of 19 ED patients and 15 healthy controls participated in this study. As shown in Table 1, there were no significant differences between the two groups in terms of age and body mass index. Regarding sex hormone levels, no significant differences were observed between the groups in testosterone, estrogen, or prolactin levels. Similarly, for biochemical indices, no significant differences were found in the levels of total cholesterol (TC) and triglycerides (TG) between the two groups.

TABLE 1

Cohort characteristicED groupControl grouptP-value
Age (y)31.16 ± 3.1129.73 ± 5.120.9480.35
BMI (kg/m2)27.37 ± 7.2123.47 ± 3.251.9440.06
IIEF-512.11 ± 4.922.53 ± 2–8.436<0.001
IPSS4.21 ± 5.262.4 ± 4.871.0300.31
Testosterone (ng/mL)5.33 ± 2.65.85 ± 1.88–0.6450.52
Estradiol (pg/mL)27.67 ± 7.3524.85 ± 3.821.4430.160
Follicle stimulating hormone (mIU/mL)5.53 ± 2.324.53 ± 1.631.4170.17
Luteinizing hormone (mIU/mL)4.81 ± 1.864.62 ± 1.330.3230.75
Prolactin (μIU/mL)219.5 ± 95.1279.59 ± 137.65–1.5040.14
Blood sugar (mmol/L)7.41 ± 4.165.37 ± 0.652.1090.05
HDL (mmol/L)1.27 ± 0.261.25 ± 0.20.1880.85
LDL (mmol/L)2.79 ± 0.512.67 ± 0.820.5130.61
Triglyceride (mmol/L)2.25 ± 2.331.68 ± 0.970.8790.39
Total cholesterol (mmol/L)5.06 ± 0.734.8 ± 1.060.8520.40
Homocysteine (mmol/L)17.86 ± 10.7323.62 ± 22.41–0.9160.37
Lipoprotein (mmol/L)162.21 ± 79.37156.52 ± 150.790.1320.90

The demographics and serum characteristics of patients with ED and healthy controls.

3.2 Alterations in the diversity and structure of gut microbiota in ED patients

We investigated changes in the composition of the gut microbiota at both the “phylum” and “genus” levels between ED patients and controls. At the phylum level, the two most abundant microbiota in the gut of both groups were Bacillota (51.1% vs. 47.6%) and Bacteroidota (26.6% vs. 28.0%). At the genus level, Prevotella (9.8%) and Bacteroides (5.1%) were the two most abundant genera in ED patients, while in the control group, Bacteroides (7.8%) had the highest abundance, followed by Prevotella (3.3%) (Figures 1a,b).

FIGURE 1

To compare the microbial community richness and homogeneity between the two groups, we analyzed the samples for alpha diversity. The results showed no significant differences in gut microbiota diversity (Shannon and Simpson indices) between the ED and control groups (p = 0.8 vs. p = 0.4). However, there was a trend toward decreasing alpha diversity in the ED patients (Figures 2A,B). Furthermore, we performed β-diversity analysis to assess the differences in microbial community composition and distribution between the groups. The results revealed no significant separation between the two groups after clustering the ED group (black) and the control group (red). The unweighted UniFrac distance also showed no statistically significant difference (R = 0.017, p = 0.26), indicating that there were no significant structural changes in the gut microbiota between the two groups (Figures 2C,D).

FIGURE 2

3.3 Identification of key gut microbiota between ED and control groups

To identify the species that most significantly contributed to the differences between the two groups, LefSe analysis was conducted on specimens from both the ED and control groups. A total of 61 gut microorganisms were found to be significantly different, with 35 species being more abundant in the ED group and 26 in the control group (Figure 3A). Furthermore, a comparison of the abundance of the main species revealed that in the ED group, Ruminococcus gnavus, Thomasclavelia ramosa, Clostridium sp. AF32-12BH, Clostridium nexile, and Eubacterium siraeum were more abundant, while the abundance of Bacteroides intestinalis was significantly reduced compared to the control group (Figures 3B–G).

FIGURE 3

3.4 Differential analysis of gut microbial metabolic pathways

Metabolic pathways of the gut microbiota were predicted, and those with differential expression between the controls and ED patients were identified through LEfSe analysis. As shown in Table 2, pathways associated with nucleotide and lipid metabolism, such as purine metabolism (ko00230) and glycerophospholipid metabolism (ko00564), were more highly expressed in the gut microbiota of ED patients. In contrast, pathways related to amino acid and vitamin metabolism, such as valine, leucine, and isoleucine degradation (ko00280) and lipoic acid metabolism (ko00785), were more abundant in the control group (P < 0.05).

TABLE 2

CharacteristicIDNameClassLDAP-value
EDko00230Purine metabolismMetabolism; nucleotide metabolism2.6970.006
ko00564Glycerophospholipid metabolismMetabolism; lipid metabolism2.3540.030
ko00983Drug metabolism - other enzymesMetabolism; xenobiotics biodegradation and metabolism2.3390.042
ko00770Pantothenate and CoA biosynthesisMetabolism; metabolism of cofactors and vitamins2.2220.036
ko05012Parkinson diseaseHuman diseases; neurodegenerative disease2.1730.033
ko05418Fluid shear stress and atherosclerosisHuman diseases; cardiovascular disease2.1570.033
Controlko00280Valine, leucine and isoleucine degradationMetabolism; amino acid metabolism2.3650.004
ko00785Lipoic acid metabolismMetabolism; metabolism of cofactors and vitamins2.0640.023
ko05230Central carbon metabolism in cancerHuman diseases; cancer: overview2.1560.021

Linear discriminant analysis effect size (lefse) analysis of two groups of pathway abundance.

3.5 Predictive modeling of gut microbes

To explore the potential of gut microbes as biomarkers for ED, we developed a diagnostic model using the random forest algorithm. The model, based on 10 species, demonstrated high performance, with AUC values of 95.45 and 100% for the training and validation sets, respectively (Figures 4A,B). These results indicate that the model can accurately differentiate between ED patients and healthy controls. Additionally, the species identified as significant by the model were highlighted using mean decreasing accuracy (Figure 4C).

FIGURE 4

4 Discussion

The human gut microbiota consists primarily of six major groups: thick-walled bacilli, anaplasma, actinobacteria, Mycobacterium, and Clostridium, with the dominant groups being anaplasma and thick-walled bacilli (). These microbiota play a key role in regulating the body’s immunity, metabolism, and endocrine processes, primarily through their own activities and metabolic products ().

We observed differences in the abundance of certain gut microbiota between the control and ED groups through macro gene sequencing. Specifically, the abundance of Ruminococcus gnavus, Thomasclavelia ramosa, Clostridium sp. AF32-12BH, Clostridium nexile, and Eubacterium siraeum species was increased in ED patients, while the abundance of Bacteroides intestinalis decreased. Previous studies have linked Thomasclavelia ramosa, a species associated with ecological disorders, to non-alcoholic fatty liver disease (). Clostridium ramosa, a Gram-positive bacterium belonging to the phylum Firmicutes that forms spores and grows in anaerobic environments (), has been found to be closely associated with fructose and bile acid metabolism, suggesting it may be a potential target for treating non-alcoholic steatohepatitis (). Zhu et al. conducted a Mendelian randomization analysis exploring the causal relationship between gut microbiota and ED, finding that Tyzzerella3 species were linked to an increased risk of ED (). In a study by Smadar Shilo et al., LefSe analysis revealed that Eubacterium siraeum was significantly enriched in individuals with type 1 diabetes (). Hu et al. compared the gut microbiome at the species level and found that Ruminococcus gnavus was enriched in the obese group (). Summarizing, the gut microbes identified in our study, which showed increased abundance in ED patients, appear to influence glucose and lipid metabolism, but their precise biological roles in ED development remain to be further explored.

Notably, combining the results from LEfSe and Metastat analyses, only the abundance of Bacteroides intestinalis was significantly reduced in the ED group (). Previous studies have shown that Enterobacter aerogenes degrades complex arabinoxylans and xylans found in dietary fibers such as wheat, rye, oats, and barley (). The degradation products, butyrate and ferulic acid, have demonstrated various functional activities, including antioxidant, anti-inflammatory, and antimicrobial effects (; ). These metabolites can interact with the immune system, modulating the development and function of almost all types of immune cells in the intestinal immune cell pool, thereby influencing immune function and preventing overactive immune responses (; ). Therefore, we speculate that the reduced abundance of Bacteroides intestinalis may contribute to a prolonged inflammatory state in the body, impairing the vascular endothelium and promoting the development of ED. However, the exact mechanism remains to be confirmed. Additionally, previous studies have found that supplementing with polysaccharides from Lyophyllum mushrooms can reduce obesity and hyperlipidemia in mice by increasing the abundance of Bacteroides intestinalis (). Thus, Bacteroides intestinalis may represent a potential therapeutic target for ED in the future.

There is some controversy regarding changes in the diversity of the gut microbiota in ED patients. Qiao et al.’s study found higher alpha diversity (Shannon index) in the ED group (), which contrasts with Geng et al.’s findings (). In contrast, Kang et al.’s study presented a third conclusion, reporting no significant difference between the ED and control groups across several indices, including species richness (Chao1 and observed OTUs) and diversity of gut microbiota (Shannon and Simpson indices) (). In the present study, the ED group showed a trend toward lower alpha diversity, but this difference was not statistically significant. Notably, unlike previous studies, no statistically significant difference in β-diversity was observed between the two groups, suggesting that there may be no significant difference in the microbiota composition between the groups. We have conducted a comprehensive analysis to explore the reasons for the differences in research results. Firstly, a small sample size may be a key factor contributing to the discrepancies. Secondly, the different stages of ED could also impact the characteristics of the gut microbiota. In the early stages of the disease, the abundance of gut microbiota decreases, leading to a significant reduction in diversity. However, as the disease progresses, the diversity of the gut microbiota may stabilize through compensatory regulation, changes in dietary structure, and drug treatments. Additionally, differences in dietary habits among populations may contribute to variability in findings (). Newman et al. noted that diet is a critical factor influencing gut bacterial diversity (). The diversity results observed in our study are similar to those of , which may be due to the fact that our study sample is from northern China, where dietary patterns are similar. Given these factors, particularly the controversial results regarding microbiota diversity in ED patients, it is essential to expand the sample size in future studies and minimize the impact of regional differences. This approach will help yield more reliable and comprehensive conclusions.

Our study compared the functional genes of the gut microbiota using the KEGG database, and the results indicated that the functions of microbiota with significant changes in ED were primarily reflected in the upregulation of genes related to atherosclerosis, purine metabolism, and glycerophospholipid metabolism, along with the downregulation of genes involved in lipoic acid metabolism, and valine, leucine, and isoleucine degradation. Branched-chain amino acids (BCAAs) leucine, isoleucine, and valine are essential amino acids that cannot be produced by the human body and must be obtained through the diet (). These amino acids play a crucial role in various metabolic processes, including energy generation and utilization, protein synthesis and breakdown, fat conversion, and sugar metabolism. As early as 1960, studies noted an increase in BCAA levels in obese patients, and with the advent of metabolomics techniques, the correlation between obesity, insulin resistance, and BCAAs has been further confirmed (; ). Vanessa et al. found that transplanting fecal microbiota from obese twins into germ-free mice resulted in higher circulating BCAA levels in the mice compared to those receiving microbiota from lean twins (). This suggests that a reduction in gut microbiota capable of biodegrading BCAAs could lead to elevated circulating BCAA levels, potentially contributing to obesity and glucose metabolism disorders. Additionally, the ED group showed significant upregulation of functional genes related to lipid metabolism, particularly glycerophospholipid metabolism. This finding is consistent with previous literature, which reports that diabetes, hypertension, dyslipidemia, and obesity are independent risk factors for ED ().

There are some limitations in this study that should be addressed in future research. First, as a pilot study, the sample size was small, and additional cases are needed to confirm our findings. Second, ED can be categorized into cardiac, organic, and mixed types based on underlying pathogenesis. However, the ED patients in this study were not specifically categorized, and thus, we could not compare microbiota changes between these subgroups. Finally, as mentioned earlier, the study sample has certain geographical characteristics. Although we established strict inclusion criteria, it is difficult to fully exclude the impact of dietary factors or other drugs on the gut microbiota.

5 Conclusion

This study found no significant difference in the diversity and structure of the gut microbiota between the ED population and healthy individuals. However, changes were observed in the abundance of certain gut microbiota, notably a significant decrease in Bacteroides intestinalis and an increase in Ruminococcus gnavus in the ED group. These microbial alterations may contribute to the increased risk of ED by promoting atherosclerosis, enhancing lipid metabolism, and inhibiting the degradation of branched-chain amino acids. Improving the intestinal microecology could potentially become a novel approach for the prevention and treatment of ED in the future.

Statements

Data availability statement

The original contributions presented in the study are publicly available. This data can be at: https://www.ncbi.nlm.nih.gov/sra/PRJNA1261102.

Ethics statement

The studies involving humans were approved by the Ethics Committee of the First Affiliated Hospital of Dalian Medical University (PJ-KS-KY-2023-167). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

QS: Formal Analysis, Methodology, Software, Writing – original draft. KW: Data curation, Investigation, Supervision Writing – original draft. YL: Formal Analysis, Methodology, Visualization, Writing – review & editing. QT: Methodology, Project administration, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

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.

Generative AI statement

The authors declare that no Generative AI was used in the creation of this manuscript.

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.

References

Summary

Keywords

Bacteroides intestinalis, erectile dysfunction, gut microbiota, metabolic pathways, macrogenomics

Citation

Su Q, Wang K, Luo Y and Tang Q (2025) Altered gut microbiota in erectile dysfunction patients: a pilot study. Front. Microbiol. 16:1530014. doi: 10.3389/fmicb.2025.1530014

Received

19 November 2024

Accepted

15 May 2025

Published

05 June 2025

Volume

16 - 2025

Edited by

George Tsiamis, University of Patras, Greece

Reviewed by

Muawuz Ijaz, University of Veterinary and Animal Sciences, Pakistan

Daniela Elena Serban, University of Medicine and Pharmacy Iuliu Hatieganu, Romania

Updates

Copyright

*Correspondence: Qizhen Tang,

†These authors have contributed equally to this work and share first authorship

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics