Distorted bacterial network in male factor infertility
To compare the seminal microbiome in men with normal or abnormal spermiogram, we recruited 88 patients seeking assisted reproduction with either abnormal (male-factor infertility, n=58) or normal spermiogram (partners of female-factor infertility or idiopathic infertility, n=30) (Supplemental Table S1). Patients with abnormal spermiogram include oligozoospermia, asthenozoospermia, teratozoospermia or combined cases. The age of the patient, the age of the partner, the treatment outcome and semen volume were comparable. However, the sperm concentration, total motility, progressive motility and morphology were significantly lowered in the male-factor infertile group.
We then performed 16S rRNA gene sequencing on seminal plasma samples from this cohort of patients. After sequencing, a sequence curation pipeline optimized for analyses of amplicon libraries was performed for quality control with a low sequencing error rate [29]. In total, 2591 amplicon sequence variants (ASVs) were identified across all the seminal plasma samples. After removing all 178 contaminants in the negative controls, 2413 ASVs were maintained. Consistent with previous reports [2], Firmicutes, Proteobacteria, Actinobacteria and Bacteroidetes were the major phyla in the seminal microbiome that constituted at least 80% of the phyla identified (Fig 1A). We then compared the richness and evenness of the bacterial community using multiple indices. In our cohort, we observed a significant difference in Weighted Unifrac distance between the normozoospermia and male-factor infertile group but no significant difference in other α- and β- diversities between the normozoospermia and male-factor infertile group were observed (Fig 1B-D, Supplemental Fig S1-2). At the genus level, we observed an enrichment of Veillonella, Gardnerella and Lactobacillus in normozoospermia patients. However, the differential abundance of these genera was not statistically significant.
To further investigate the structure and interaction of seminal microbial communities, we performed community network analysis using the Weighted Correlation Network Analysis (WGCNA) algorithm [30], which has not been applied in the seminal microbiome. In this analysis, the nodes represent ASVs and the edges that connect these nodes represent correlations between ASVs. Notably, while the number of nodes was comparable, we observed a marked decrease in the number of edges in male-factor infertile patients (Fig 1E). In microbial community networks, the ASVs clustered into independent modules, a property known as modularity, with a small group of ASVs serving as module connectors. We observed the absence of a module hub and the disappearance of peripheral nodes in three phyla, Fusobacteria, TM7, and Spirochaetes in male-factor infertile patients which lead to the shrinkage of network diameter and heterogeneity (Fig 1E-F). These data suggest that the bacterial network is distorted in patients with abnormal spermiogram despite the comparable richness and evenness of the metacommunities.
Elevation of DEFB119 stratifies the dysbiosis in male-factor infertility
Previous studies have utilized sperm parameters as stratification factors in the analysis of seminal microbiomes. We speculated that a host factor potentially involved in the host-microbe interactions would provide a better stratification of the seminal microbiome and a higher resolution of the dysbiosis associated with male infertility. Therefore, we explored the involvement of β-defensins in regulating the seminal microbiome, we focused on DEFB119, a β-defensin that plays pivotal roles in sperm production and functions [18, 19]. We determined the protein level of DEFB119 in seminal plasma by enzyme-linked immunoassay with an antibody against the C-terminus of the protein. We observed a range of DEFB119 levels in subjects with normal spermiogram (mean 291.30, CI 210.50 – 372.00 ng/ml). Intriguingly, while the patients with abnormal spermiogram expressed a subtle increased level of DEFB119 in the seminal plasma (mean 393.60, CI 245·70 - 541·50 ng/ml), a subgroup of patients demonstrated a marked elevation of DEFB119 in seminal plasma above the 100th centile (>900 ng/ml) of subjects with normal spermiogram (Fig 2A & Table 1).
To examine the metacommunity structure associated with elevated DEFB119, we categorized the patients according to the following grouping: G1 - Normal spermiogram and low DEFB119 level (n=30); G2 - Abnormal spermiogram and low DEFB119 (n=52); and G3 - Abnormal spermiogram and elevated level of DEFB119 (n=5). We observed a lowered sperm concentration, motility and morphology in G2 and G3 as compared to G1. However, the spermiogram parameters were comparable in male-factor infertile patients with normal (G2) or elevated levels of DEFB119 (G3)(Table 1). The phylotypes in G3 were dominated by Firmicutes and Proteobacteria while the abundance of Actinobacteria and Bacteroidetes was diminished (Fig 2B). Comparing the top forty most abundant genera, the abundance of Bactobacillus, Gardnerella, and Prevotella decreased (Fig 2C). We have also compared the richness and evenness of the bacterial community among the three groups using multiple indices as in our previous analysis. The α-diversity was significantly reduced in patients with elevated levels of DEFB119 as compared to G1 and G2 groups (Shannon index p < 0.05, Fig 2D and Supplemental Fig S3). Similarly, the β-diversity as measured by Jaccard distance, Bray Curtis Distance and Weighted Unifrac distance showed that the bacterial community in patients with high levels of DEFB119 (G3) was significantly different from those in G1 and G2 groups (p < 0.05, Fig 2E-F, Supplemental Fig S4). No significant difference in α- and β- diversities except Weighted Unifrac distance were observed between G1 and G2 groups.
Differential abundance analysis at the genus level revealed a decrease in Clostridium and increases in eight genera in G3 (p < 0.05), including Sporosarcina, AF12, Helicobacter, Desulfovibrio, Phyllobacterium, Enterobacter, Anaerobacillus and Carnobacterium, despite their rare occurrence and low relative abundance (Supplemental Fig S5). WGCNA community network analysis showed that patients with elevated levels of DEFB119 (G3) showed a marked decrease in the numbers of both nodes and edges (Fig 2G). The network diameter and heterogeneity were further diminished in G3 as compared to G2. These data suggest that the elevated level of DEFB119 is associated with dysbiosis of the seminal microbiome and severe distortion of bacterial networks in male-factor infertile patients.
Mediators of abnormal sperm parameters in male infertility
Members of the β-defensin family are known to regulate sperm functions required for the migration in the female reproductive tract and successful fertilization [14, 15]. While our data showed that the elevated level of DEFB119 was associated with the dysbiosis of the seminal microbiome in male infertility, the infertile outcome could be attributed to the effect of DEFB119 on sperm functions per se, the indirect effect from the dysbiosis of the seminal microbiome or both. To study this, we examined the correlation of DEFB119 level with various sperm parameters in normozoospermia and male-factor infertile patients. Among the sperm parameters, a significant negative correlation was observed in progressive motility (Supplemental Fig S6, p < 0.05). In line with this, redundancy analysis (RDA) also revealed a negative correlation of DEFB119 with sperm motility and progressive motility (Fig 3A). Although genera such as Gardnerella and Campylobacter were positively associated with sperm concentration, Veillonella was positively associated with sperm motility and Prevotella was positively associated with sperm morphology and concentration, these associations were not significant (Fig 3A). Next, we set out to test if the abnormal sperm parameter was mediated by the elevation of DEFB119 and the dysbiosis of the seminal microbiome. Mediation analysis revealed a notable increase in the number of significant mediations in the male-factor infertile patients but not the normozoospermia control cohort (Fig 3B). Unexpectedly, significant mediations were observed in 4 out of 6 tested models, regardless of the initiators and mediators. These results suggest that the elevation of DEFB119 and the dysbiosis of the seminal microbiome could be the cause or the effect and vice versa. Nonetheless, the RDA and mediation analysis suggests the possible involvement of elevated DEFB119 and dysbiotic seminal microbiome in mediating the abnormal sperm parameters.
Elevated DEFB119 decreases sperm motility in male infertility
To validate the mediation analysis, we set out to examine the effect of elevated DEFB119 on sperm parameters. We mimicked the elevated level of DEFB119, as observed in G3, by recombinant DEFB119 (rDEFB119) treatment in a separate cohort of patients with normal or abnormal spermiogram profiles and normal level of DEFB119 in the seminal plasma i.e. being classified as G1 or G2. Since the decrease in sperm count and morphology are spermatogenic factors that would not be altered in a short incubation period of rDEFB119 treatment theoretically and both progressive motility and total motility demonstrated a negative correlation with DEFB119 in RDA (Fig 3A), we examined if the elevated DEFB119 affects the motility of sperm. Our results showed that rDEFB119 significantly decreased both total and progressive motility when compared with the vehicle control (Fig 4A-B, p<0·001). Interestingly, when the sample cohort was further categorized into normozoospermic and male-factor infertile groups based on the spermiogram profile, the decrease in progressive motility was only observed only in male-factor infertile patients (Fig 4C-D). These results, in corroboration with the mediation analysis, suggest that elevated DEFB119 level has detrimental effects on sperm motility, and may contribute to the infertile outcome of the male-factor infertile patients with disturbed seminal microbial networks.
Species-specific antimicrobial activity of DEFB119 shapes the seminal microbiome
It is well established that the β-defensins possess antimicrobial activity that contributes to the host defence against pathogens [31–35]. Of note, reproductive-tract-specific β-defensins play dual roles in host defence and sperm functions in the reproductive tract, their expressions are known to be induced by both physiological and pathological stimuli including hormone and bacterial toxin lipopolysaccharides [36, 37]. Therefore, the elevated DEFB119 in seminal plasma could be a cause or an effect of the dysbiosis of the seminal microbiome. To investigate if DEFB119 plays an active role in shaping the seminal microbiome and provoking dysbiosis, we performed bactericidal assay against two dominant genera observed in G1 and G2, Prevotella and Streptococcus, and the genus found in G3, Pseudomonas (Fig 5A-C). We include species that can be propagated in vitro and have either been reported in the semen, including Prevotella copri, Streptococcus mitis, Streptococcus anginosus and Pseudomonas aeruginosa [38–41], or found in other organ systems such as Prevotella stercorea in the gut microbiome and Streptococcus intermedius in the central nervous system [42, 43].
We treated the bacterial culture with various doses of rDEFB119 and observed the bacterial growth after culture. We observed a dose-dependent decrease in the amount of Prevotella stercorea, Streptococcus mitis, Streptococcus anginosus and Streptococcus intermedius after rDEFB119 treatments as compared to the vehicle controls, suggesting the bactericidal effect of DEFB119 on these species (Fig 5D-I). Intriguingly, rDEFB119 significantly promoted the growth of Streptococcus mitis at 3 ug/ml while higher doses demonstrate significant inhibitory effects. This result suggests that reproductive tract β-defensins can shape the seminal microbiome by exerting both promoting and inhibitory effects on specific species in a dose-dependent manner. Recombinant DEFB119 exerted a negligible effect on the growth of Pseudomonas aeruginosa and Prevotella copri regardless of the dosage used (Fig 5D-I), suggesting that these bacteria were resistant to DEFB119. Taken together, our results suggest that the elevation of DEFB119 in semen plays ab active role in provoking the dysbiosis of the seminal microbiome.