Skip to main content
Download PDF

Bioinformatics analysis identifies potential hub genes and crucial pathways in the pathogenesis of asthenozoospermia

BMC Medical Genomics volume 15, Article number: 252 (2022) Cite this article

Abstract

Background

Asthenozoospermia is a troublesome disease experienced by men in their reproductive years, but its exact etiology remains unclear. To address this problem, this study aims to identify the hub genes and crucial pathways in asthenozoospermia.

Methods

We screened two Gene Expression Omnibus (GEO) datasets (GSE92578 and GSE22331) to extract the differentially expressed genes (DEGs) between normozoospermic and asthenozoospermic men using the “Limma” package. Gene enrichment analyses of the DEGs were conducted using the “clusterProfiler” R package. The protein-protein interaction (PPI) network was then established using the STRING database. A miRNA-transcription factor-gene network was constructed based on the predicted results of hub genes using the RegNetwork database. The expression of four hub genes in asthenozoospermia and normal samples were verified using quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) and western blotting.

Results

We identified 271 DEGs, which included 218 upregulated and 53 downregulated in two asthenozoospermia datasets. These DEGs were observed to be markedly enriched in pathways with cell growth and embryonic organ development, phospholipase D signaling pathway, cGMP-PKG signaling pathway, and Wnt signaling pathway. The most significant genes were identified, including COPS7A, CUL3, KLHL7, NEDD4. We then constructed regulatory networks of these genes, miRNAs, and transcription factors. Finally, we found that the COPS7A was significantly upregulated in patients with asthenozoospermia, but CUL3, KLHL7 and NEDD4 were significantly downregulated compared with normal samples.

Conclusion

We applied bioinformatics methods to analyze the DEGs of asthenozoospermia based on the GEO database and identified the novel crucial genes and pathways in this disease. Our findings may provide novel insights into asthenozoospermia and identify new clues for the potential treatment of this disease.

Peer Review reports

Introduction

Male infertility has become a primary human health concern due to its increasing incidence and unknown etiology[1]. This disease not only introduces great psychological pressure and a substantial economic burden to the patients but also affects their overall health[2]. Asthenozoospermia accounts for a large proportion of male infertility patients, and its main clinical manifestations are decreased sperm motility. Based on standards established by the World Health Organization (WHO), asthenozoospermia is characterized by a total sperm motility lower than 40% and a progressive sperm motility lower than 32% in a semen sample[3]. The process of spermatogenesis depends on the accurate regulation of key genes, and improper regulation of this process leads to abnormalities in spermatogenesis[4]. Several factors, including smoking, drugs, chromosomal abnormalities, or mutations, are associated with asthenozoospermia, but the etiology of this disease remains unclear[5]. Several studies have found high mRNA levels in human testicular tissue and sperm, and abnormal mRNA expression is closely related to poor sperm quality[6, 7]. Although great progress toward elucidating the etiology of male infertility has been achieved, the information on the pathogenesis of asthenozoospermia remains scarce, and a large number of genes with unknown functions need further study.

With advancements in microarray and high-throughput technology, bioinformatics analyses have recently been performed to understand the molecular mechanism of asthenozoospermia. Abnormalities at the genomic level, including genes and signaling pathways related to this disease, have increasingly been identified between fertile men and patients with asthenozoospermia[8,9,10,11,12]. The large amount of gene sequencing data available has provided a convenient resource for analyzing the regulation of gene expression in asthenozoospermia using increasingly complex analytical techniques[13]. In the current manuscript, we obtained differentially expressed genes (DEGs) in patients with asthenozoospermia using the GSE92578 and GSE22331 datasets from the Gene Expression Omnibus (GEO) database. Subsequently, we conducted a series of bioinformatics analyses to identify crucial genes and pathways involved in spermatogenesis. Our findings indicated the newly discovered crucial genes and key pathways may contribute to finding the new therapeutic targets and approaches for this disease.

Methods

Sample preparation and semen collection

After obtaining informed consent from these patients and healthy volunteers, we collected 30 healthy young men and 30 patients diagnosed with asthenozoospermia in our hospital from November 2021 to January 2022. Healthy volunteers with varicocele, reproductive tract infection, and abnormal sperm parameters were excluded. All subjects were Han Chinese with no genetic association. Three ml semen samples obtained through masturbation were analyzed after three to five days of abstinence. Asthenozoospermia is defined as progressive sperm motility lower than 32% in a semen sample, which is confirmed by routine semen analysis of three semen samples collected at different time points[14]. Semen samples containing the following parameters are considered normal: semen volume ≥ 2 mL, sperm concentration ≥ 15 × 106/mL, pH ≥ 7.2, progressive motility (PR) ≥ 32%, PR + Non-progressive motility (NP) ≥ 40%. All semen samples were collected and processed by the guidelines of the World Health Organization (WHO) Laboratory Manual on Human Semen Examination and Processing (5th Edition 2010).

RNA extraction and quantitative real-time polymerase chain reaction (qRT‐PCR) analysis

According to a previous report, the method of TRIzol (Invitrogen) extracted RNA from sperm cells [15]and then, reverse transcribed it into cDNA based on a reverse transcription kit (Novabio, China). Subsequently, experiments were performed in an SYBR qPCR Mix kit (Novabio, China) protocol using glyceraldehyde-3-phosphate dehydrogenase (GAPDH) expression as a reference using the CFX96 system (Bio-Rad Laboratories, Inc.). Real-time RT-PCR was carried out at 95 °C for the 30s, then at 95 °C for 10s, with 40 cycles, annealing, and extension at 60 °C for 30s. All samples were normalized according to GAPDH, and the oligonucleotide primers for five genes are provided in Table 1. The relative expression of these genes was measured using the 2 ΔΔ Ct method[16].

Table 1 Primer sequences used for qRT-PCR amplification
Full size table

Sperm protein extraction and western blotting

The sperm cells in the semen sample was isolated and subsequently processed according to the previously method[17]. The samples were separated by SDS-PAGE in a 10% polyacrylamide gel and transferred to PVDF membranes. After the PVDF membrane was blocked, anti-COPS7A (cat. no.DF8375; 1:1000; Affinity Biosciences), anti-CUL3 (cat. no. DF6223; 1:1000; Affinity Biosciences), anti-KLHL7 (cat. no. A7817; 1:1000; ABclonal Technology), anti-NEDD4 (cat. no.AF4636; 1:1000; Affinity Biosciences), and anti-GAPDH (cat. no.AF7021; 1:1000; Affinity Biosciences) antibodies were incubated separately. During blotting, the blots were cut prior to hybridisation with antibodies based on the molecular weight. It was subsequently incubated with HRP-coupled goat anti-rabbit IgG (cat. no. S0001; 1:10,000;Affinity Biosciences) and goat anti-mouse IgG (cat. no. S0002; 1:10,000;Affinity Biosciences) secondary antibodies. The reaction band is detect by enhanced chemiluminescence (ECL kit, Thermo Fisher Scientific).

Data collection

In our study, two asthenozoospermia datasets were selected from the GEO database based on previous literature reports, and these data proved to be reliable[13, 15]. Two asthenozoospermia GSE profiles include GSE92578 (GPL20301 Illumina HiSeq 4000 [Homo sapiens]), containing 3 asthenozoospermia sample and 4 normal sample[9], and GSE22331 (GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array), containing one asthenozoospermia and normal sample[12]. To identify hub genes in asthenozoospermia, we filtered top genes among the DEGs based on protein-protein interaction (PPI) networks and Friends analyses.

Data processing of DEGs

The process used in this study is shown in Fig. 1. From the GSE22331 dataset, DEGs were screened using R software (v 4.1.1) based on the “Limma” package based on preset thresholds of |logFC| > 1.0 and P-value < 0.05[18]. The DEG results from the GSE92578 dataset were downloaded from the supplementary material of a previously published study (https://academic.oup.com/biolreprod/article/100/4/982/5230853#134234708)[9]. DEGs were defined as genes that met the criteria Q-value < 0.001 and |logFC | > 1.0 in the GSE92578 dataset. Intersecting genes were obtained by Venn analyses[19]. A visual hierarchical cluster analysis was also conducted to display a volcano plot of the DEGs.

Fig. 1
figure 1

Study design and data preprocessing

Full size image

GO annotation and KEGG pathway enrichment analyses of DEGs

We performed GO annotation and KEGG pathway of DEGs[20,21,22]. The GO annotation (including biological processes, cellular components, and molecular functions) and KEGG pathway (P-value < 0.05) were enriched using the “clusterProfiler” R language package[23].

Top genes based on the semantic similarities of GO terms

Subsequently, we analyzed the functional similarity between the common DEGs using the “GOSemSim” package[24], which analyses the semantic similarity of GO terms obtained for the genes (also called Friends analysis). The top 24 DEGs with high average functional similarity were considered hub genes in the Friends analysis.

Construction of PPI network and screening of hub genes

To explore the correlation among DEGs, the PPI network was built by the STRING database (https://www.string-db.org/). We then visualized the results using Cytoscape software (v 3.7.2), and the top 24 DEGs were identified using the cytoHubba based on node degree[25]. To identify hub genes in asthenozoospermia, we filtered top genes among the DEGs based on protein-protein interaction (PPI) networks and Friends analyses.

Construction of the miRNA-TFs-genes network

To study the regulatory mechanisms of four crucial genes, miRNAs and TFs that bind to hub genes were predicted using RegNetwork (http://www.regnetworkweb.org/), which is a target gene prediction database. We then constructed regulatory networks of these hub genes, including TF-miRNA-gene regulatory networks, using Cytoscape software (v 3.7.2).

Statistical analysis

Experimental data were reported by x ± SD, and the differences among the experimental groups were compared by Student’s t-test. All experimental data analysis was performed using SPSS 23.0 (IBM Corp.), and P < 0.05 indicated statistical significance.

Results

DEGs in asthenozoospermia

The analysis of the GSE22331 dataset successfully identified 4,621 DEGs, which included 3,019 upregulated and 1,602 downregulated genes (Additional file 1: Tables S1 and S2). The DEGs identified from the GSE22331 database are shown in Fig. 2 A, where the red and green dots represent the up-and down-regulated DEGs, respectively. From the GSE92578 dataset, we identified 4,259 DEGs, which included 3,061 upregulated and 1,198 downregulated genes (Additional file 1: Tables S3 and S4), and the DEGs identified from the GSE92578 database are shown in Fig. 2B. Venn analysis was then performed to compare the DEGs identified from the GSE22331 dataset with those obtained from the GSE92578 dataset. In this study, we identified 271 DEGs from two datasets involving 218 up-regulated and 53 down-regulated genes (Fig. 2 C and Fig. 2D; Additional file 1: Tables S5 and S6).

Table 2 Top 24 genes based on the Friends analysis and PPI network analysis
Full size table
Fig. 2
figure 2

Identification of common DEGs from GSE22331 and GSE92578 datasets. A Volcano plot of the 4,621 DEGs in GSE22331 dataset. B Volcano plot of the 4,259 DEGs in GSE92578 dataset.Venn diagram of C upregulated and D downregulated DEGs based on the two GEO datasets. Red, upregulation; green, downregulation. DEG, differentially expressed gene; GEO, Gene Expression Omnibus

Full size image
Table 3 The basic characteristics and sperm parameters in patients with asthenozoospermia and healthy volunteers
Full size table

GO annotation and KEGG pathway enrichment analysis

Functional enrichment analysis was performed on two datasets from the GEO database, respectively. In the GSE22331 dataset, the crucial biological functions were mainly enriched in related processes of cell signal transduction (Additional file 1: Fig. S1A-C), and the crucial pathways were mainly enriched in Cytokine-cytokine receptor interaction, Herpes simplex virus 1 infection, Neuroactive ligand-eceptor interaction, Toll − like receptor signaling pathway, and Viral protein interaction with cytokine and cytokine receptor (Additional file 1: Fig.S1D-E). In the GSE92587 dataset, the crucial biological functions were mainly enriched in cell growth and signal transduction (Additional file 1: Fig. S2A-C), and the crucial pathways were mainly enriched in Axon guidance, Focal adhesion, JAK − STAT signaling pathway, PD − L1 expression and PD − 1 checkpoint pathway in cancer, Ubiquitin mediated proteolysis (Additional file 1: Fig. S1D-E).

To clarify the mechanism of these DEGs in the pathogenesis of asthenozoospermia, the biological roles of the 271 DEGs were enriched by GO and KEGG pathway. We present the top 20 significant terms of GO and pathways. The crucial biological processes of these DEGs were mainly enriched in cell growth, embryonic organ development, regulation of cell growth, sodium ion transport, organic acid transmembrane transport, and carboxylic acid transmembrane transport (Fig. 3 A). The cellular components analysis of these DEGs were significantly enriched in postsynaptic specialization, postsynaptic density, neuron-to-neuron synapse, asymmetric synapse, focal adhesion, cell-substrate junction, neuron spine, dendritic spine, and fibrillar center (Fig. 3B). These DEGs were significantly enriched in the following molecular functions: active transmembrane transporter activity, secondary active transmembrane transporter activity, organic anion transmembrane transporter activity, organic acid transmembrane transporter activity, carboxylic acid transmembrane transporter activity, and amino acid transmembrane transporter activity (Fig. 3 C). In addition, the pathways in which the DEGs were significantly enriched were the phospholipase D signaling pathway, cGMP-PKG signaling pathway, Wnt signaling pathway, mTOR signaling pathway, apelin signaling pathway, glucagon signaling pathway, and VEGF signaling pathway (Fig. 4 A). We also identified the underlying connection between the DEGs and the top 5 KEGG pathways, namely, herpes simplex virus 1 infection, phospholipase D signaling pathway, cGMP-PKG signaling pathway, Wnt signaling pathway, and platelet activation (Fig. 4B).

Fig. 3
figure 3

GO annotation enrichment analysis of DEGs. The top 20 enriched GO A BP, B CC and C MF terms. GO, gene ontology; DEG, differentially expressed gene; BP, biological process; CC, cellular component; MF, molecular function

Full size image
Fig. 4
figure 4

KEGG pathway enrichment analysis of DEGs. A The top 20 enriched KEGG pathways. B The underlying connection between DEGs and pathways in the top 5 KEGG pathway. KEGG, Kyoto encyclopedia of genes and genomes; DEG, differentially expressed gene

Full size image

Identification of top genes based on the Friends analysis

To identify interacting proteins in asthenozoospermia, proteins were ranked according to the average functional similarity relationship between proteins in the interactoms. The distribution of the functional similarities is summarized as box plots. The Friends analysis results showed that the top 24 genes were ZNF490, ZNF17, SAP30BP, ZNF783, COPS7B, COPS7A, ZNF415, MLLT1, DAB2, MRGBP, TTI1, OLFM2, TRAF3, CUL3, JADE1, CTNNBIP1, STRN4, PRKCH, NEDD4, RBM4, GLO1, MIA3, HDX, and KLHL (Fig. 5).

Fig. 5
figure 5

Summary of functional similarities of the interactome in asthenozoospermia. The box plot represents the middle 50% similarity (the upper and lower boundaries represent the 75th and 25th percentiles). The lines in the boxes represent the average of the functional similarities. Proteins with higher mean functional similarities (cutoff > 0.75) are considered to be central proteins in interacting groups in asthenospermia. The dashed line indicates the cutoff value

Full size image

Construction of protein-protein interaction (PPI) network and screening of hub genes

To understand the interaction between these DEGs, we performed PPI network analysis to screen of hub genes in asthenozoospermia. The PPI network included 175 nodes and 490 edges (Fig. 6 and Table S7 ). Subsequently, the most valuable module was identified from the PPI network using MCODE, and the top 24 genes were screened by maximal clique centrality algorithm (MCC) in the cytoHubba (Fig. 7 A). A Venn analysis was used to identify four key genes based on two methods, including COPS7A, CUL3, KLHL7, NEDD4, which may be potentially important genes in asthenozoospermia (Fig. 7B; Table 2). The mRNA expression of COPS7A was upregulated in patients with asthenozoospermia, but the expression of CUL3, KLHL7, and NEDD4 was downregulated in these patients (Fig. 7 C).

Fig. 6
figure 6

PPI network was constructed by all the 271 DEGs using STRING database. PPI, protein-protein interaction; DEG, differentially expressed gene; STRING, search tool for the retrieval of interacting genes

Full size image
Fig. 7
figure 7

The hub genes identification. A Module analysis of PPI network. B The 4 common hub genes between the PPI network and Friends analysis were screened. C The heat map showed the expression of four hub genes in the GSE22331 and GSE92578 datasets. PPI, protein-protein interaction.*P < 0.05, ****P < 0.001

Full size image

Construction of the miRNA-TF-genes network

Then, we predicted miRNAs and transcription factors (TFs) upstream of the four hub genes in the RegNetwork database. Regulatory networks of the four hub genes were subsequently constructed using Cytoscape (Fig. 8 and Additional file 1: Table S8). In this network, miR-570 and four TF (SPI1, MAX, GABPA, and MYC), were verified to be coregulators of COPS7A and CUL3, and miR-577 was confirmed as a coregulator of CUL3 and KLHL7. In addition, the TF of CTCF was identified as a coregulator of KLHL7 and NEDD4, the TF of PAX2 and TBP were proven to be coregulators of NEDD4 and CUL3, and the TF of GATA1 and EGR1 were found to be coregulators of COPS7A and NEDD4.

Fig. 8
figure 8

The miRNAs-TFs-genes networks. microRNAs, miRNAs; TFs, transcription factor

Full size image

Validation expression of the hub genes by qRT-PCR and western blotting.

The basic characteristics and sperm parameters from asthenozoospermia and normal sample were listed in Table 3. Expressions of four hub genes were validated by qRT-PCR in asthenozoospermia and normal samples, which was consistent with bioinformatical results (Fig. 9 A). A western blotting assay using 16 semen samples (including 8 asthenospermia and 8 normal samples) revealed similar results, indicating COPS7A was significantly upregulated in patients with asthenozoospermia, but CUL3, KLHL7 and NEDD4 were significantly downregulated compared with normal samples (Fig. 9B and Additional file 1: Fig. S3).

Fig. 9
figure 9

Hub genes expression in normozoospermia and asthenozoospermia samples

A Expressions of four hub genes were validated by qRT-PCR in spermatozoa from normozoospermia and asthenozoospermia. B Expressions of four hub genes were validated by western blotting experiments in spermatozoa from normozoospermia and asthenozoospermia. *P < 0.05. Asth, asthenozoospermia; Noth, normozoospermia

Full size image

Discussion

Approximately 15% of human couples worldwide are reportedly infertile, and more than half of these cases have been linked to male infertility[26]. Clinically, the main manifestation of this disease is poor sperm motility. For example, sperm functional or structural defects and seminal plasma abnormalities may impair sperm motility[1]. Therefore, exploring the molecular mechanisms of asthenozoospermia will help us identify breakthrough treatments for male infertility. In this study, we screened novel hub genes and crucial pathways in asthenozoospermia by comprehensive bioinformatics analysis.

Screening DEGs in patients with asthenozoospermia can be used to discover the pathogenesis of this disease. These DEGs are related to various molecular functions, especially the up-regulated DEGs, including cell growth, embryonic organ development, metabolism and transportation. Our results describe a general trend to up-regulation of the DEGs in asthenozoospermia, but previous investigations concur that asthenozoospermia is mostly characterized by gene down-regulation[8, 11, 15]. The differences may be associated with different data processing and diverse molecular functions. Together, these results suggest that dysregulation of key gene expression is an important pathogenesis of asthenozoospermia, as these genes are essential for sperm growth and development and energy metabolism. Sperm capacitation depends on a series of biochemical processes that occur in the female reproductive tract, which include the activation or inhibition of key signaling pathways[27]. Activation of the phospholipase D signaling pathway is thought to promote sperm capacitation[28]. A C-type natriuretic peptide secreted by the female reproductive tract may stimulate cGMP/PKG signaling in sperm cells, promote sperm activation, and induce acrosome reactions[29]. The Wnt signaling pathway has been proved to be involved in the regulation of GSK3 activity in sperm and affects sperm activity[30]. The Wnt signaling pathway has been proved to be involved in the regulation of GSK3 activity in sperm and affects sperm activity[31]. Although it has been reported that the Apelin signaling pathway, Glucagon signaling pathway, and VEGF signaling pathway are involved in the cell energy metabolism, proliferation, and apoptosis, their role in sperm function is unclear[32,33,34]. Studies have shown that sperm motility is regulated by the flagellum[35]. Defects in the function or structure of the flagellum are associated with abnormal energy metabolism and impairments in epididymal maturation and capacitation, and these effects lead to a decline in sperm motility and asthenozoospermia[36]. Taken together, our findings suggest that there is a potential correlation between these abnormal pathways and asthenozoospermia.

Sperm are highly functional cells that exhibit complex functions during reproduction[37]. The orderly expression of key genes during sperm development is crucial, and abnormal gene expression can lead to poor sperm quality or functional defects[36, 38, 39]. Moreover, the mRNA expression detected in testicular tissue was found to be consistent with that detected in sperm cells, which suggests that sperm mRNA may have potential roles to predict spermatogenesis and elucidate sperm function[7]. In this study, we studied mRNAs that exhibit significant dysregulation in patients with asthenozoospermia to explore the potential pathogenesis of this disease. In this manuscript, we identified novel key genes in asthenozoospermia based on two different algorithms. The intersection among two groups of hub genes included COPS7A, CUL3, KLHL7, and NEDD4, which can be considered the main key genes in asthenozoospermia. Moreover, evidence from the protein expression profiles of human sperm indicated that sperm protein exhibits a significantly higher or lower expression in asthenozoospermia[40]. In our study, the expression of these key genes was verified by collecting samples from asthenospermia patients, and the results were consistent with this hypothesis.

During spermatogenesis, protein ubiquitination has been shown to be involved in the stability of sperm cell structure and function[41,42,43]. The four hub genes are related to ubiquitin conjugation and protein degradation. COPS7A, the most significant gene, belongs to the COP9 signalosome (CSN) complex, which is an evolutionarily conserved multisubunit protease that regulates the activity of the ubiquitin conjugation pathway[44]. The interaction of COPS7A and polyamine regulator 1 (PMF-1) is reportedly involved in the transcriptional expression of growth-related genes[45]. We observed that COPS7A expression was upregulated in asthenospermia, but the correlation between its expression and spermatogenesis remains to be understood. CUL3 encodes a member of the cullin protein family and is involved in the regulation of multiple ubiquitin pathways and cell cycles[46, 47]. The expression of CUL3 in testicular tissue is reportedly involved in mediating protein ubiquitination in spermatogenesis by binding to KLHL10 to form the ubiquitin E3 ligase[48]. Moreover, KHLH7 has been reported to be involved in regulating nucleolar integrity and promoting terminal uryltransferase 1 ubiquitination[49]. Zhou et al. found that NEDD4 expression is critical for maintaining spermatogonial stem cell homeostasis and stress response[50]. Although the bioinformatics evidence indicates that the expression disorder of these genes can affect the development and function of male sperm, it is noteworthy that our results still need to be verified by further in vivo experiments.

Increasing lines of evidence show that sperm RNAs, including various mRNAs, miRNAs, and lncRNAs, play a crucial role in regulating male fertility[51]. Seminal plasma contains large amounts of miRNAs, and their epigenetic signatures in spermatogenesis remain unclear. These miRNAs participate in the spermatogenesis process, can reflect the pathophysiological changes of sperm cells, and serve as a potential biomarker[52]. Moreover, the synergistic effect of TFs and miRNA in the epigenetic regulatory network may be involved in the functional regulation mechanism of target genes[53]. In the current study, we constructed a miRNA-mRNA-TF network, and the interaction of all the centers provided new clues for exploring the potential mechanisms between miRNAs and their target mRNAs in asthenozoospermia. This network provides a new method for elucidating the potential mechanisms and molecular targets of asthenozoospermia.

Our study has some limitations. The two asthenozoospermia datasets retrieved from the GEO database have very small sample sizes, which could result in poor repeatability of the results. Due to insufficient clinical samples and funds, the pathogenic mechanism of four hub genes in the asthenozoospermia could not be verified. Moreover, the mechanism through which miRNAs and TFs participate in the regulation of hub genes has not been determined experimentally. We will conduct experimental verification in the future to confirm the potential mechanisms of the key genes and signaling pathways identified in this study.

Conclusion

In summary, we firstly identified 271 DEGs that are potentially involved in asthenozoospermia based on the GEO database. Importantly, the interaction between the identified DEGs identified a regulatory pathway for asthenozoospermia, which includes the phospholipase D signaling pathway, the cGMP-PKG signaling pathway, and the Wnt signaling pathway. Among the identified DEGs, four novel key DEGs were observed in the pathogenesis of asthenozoospermia. Understanding the key genes and signaling pathways in asthenozoospermia can help determine the causes of male infertility, which will aid the search for new therapies and strategies for male contraception.

Availability of data and materials

Publicly available datasets were analyzed in this study, and these datasets can be found here: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE92578; https://academic.oup.com/biolrepr.

od/article/100/4/982/5,230,853#134,234,708.

Abbreviations

GEO:

Gene Expression Omnibus

DEGs:

Differentially expressed genes

PPI:

Protein-protein interaction

qRT-PCR:

quantitative real-time reverse transcription-polymerase chain reaction

WHO:

World Health Organization

PR:

Progressive Motility

NP:

Non-progressive motility

MCC:

maximal clique centrality algorithm

References

  1. Krausz C, Riera-Escamilla A. Genetics of male infertility. Nat Rev Urol. 2018;15(6):369–84.

    Article  CAS  Google Scholar 

  2. Jensen TK, Jacobsen R, Christensen K, Nielsen NC, Bostofte E. Good semen quality and life expectancy: a cohort study of 43,277 men. Am J Epidemiol. 2009;170(5):559–65.

    Article  Google Scholar 

  3. World Health Organization. WHO laboratory manual for the examination and processing of human semen. 2010.

  4. Neto FT, Bach PV, Najari BB, Li PS, Goldstein M. Spermatogenesis in humans and its affecting factors. Semin Cell Dev Biol. 2016;59:10–26.

    Article  Google Scholar 

  5. Gilany K, Moazeni-Pourasil RS, Jafarzadeh N, Savadi-Shiraz E. Metabolomics fingerprinting of the human seminal plasma of asthenozoospermic patients. Mol Reprod Dev. 2014;81(1):84–6.

    Article  CAS  Google Scholar 

  6. Zhang Y, Shi J, Rassoulzadegan M, Tuorto F, Chen Q. Sperm RNA code programmes the metabolic health of offspring. Nat Rev Endocrinol. 2019;15(8):489–98.

    Article  Google Scholar 

  7. Ostermeier GC, Dix DJ, Miller D, Khatri P, Krawetz SA. Spermatozoal RNA profiles of normal fertile men. Lancet. 2002;360(9335):772–7.

    Article  CAS  Google Scholar 

  8. Caballero-Campo P, Lira-Albarrán S, Barrera D, Borja-Cacho E, Godoy-Morales HS, Rangel-Escareño C, Larrea F, Chirinos M. Gene transcription profiling of astheno-and normo-zoospermic sperm subpopulations. Asian J Androl. 2020;22(6):608–15.

    Article  CAS  Google Scholar 

  9. Zhang X, Zhang P, Song D, Xiong S, Zhang H, Fu J, Gao F, Chen H, Zeng X. Expression profiles and characteristics of human lncRNA in normal and asthenozoospermia sperm†. Biol Reprod. 2019;100(4):982–93.

    Article  Google Scholar 

  10. Bansal SK, Gupta N, Sankhwar SN, Rajender S. Differential genes expression between fertile and infertile spermatozoa revealed by transcriptome analysis. PLoS ONE. 2015;10(5):e127007.

    Article  Google Scholar 

  11. Jodar M, Kalko S, Castillo J, Ballescà JL, Oliva R. Differential RNAs in the sperm cells of asthenozoospermic patients. Hum Reprod. 2012;27(5):1431–8.

    Article  CAS  Google Scholar 

  12. Chen J, Wang Y, Wei B, Lai Y, Yan Q, Gui Y, Cai Z. Functional expression of ropporin in human testis and ejaculated spermatozoa. J Androl. 2011;32(1):26–32.

    Article  Google Scholar 

  13. Li L, Chen S. Screening, identification and interaction analysis of key MicroRNAs and genes in Asthenozoospermia. Int J Med Sci. 2021;18(7):1670–9.

    Article  CAS  Google Scholar 

  14. Lu JC, Huang YF, Lü NQ. WHO laboratory manual for the examination and processing of human semen: its applicability to andrology laboratories in China. Zhonghua Nan Ke Xue. 2010;16(10):867–71.

    Google Scholar 

  15. Liu XX, Cai L, Liu FJ. An in silico analysis of human sperm genes associated with asthenozoospermia and its implication in male infertility. Medicine (Baltim). 2018;97(49):e13338.

    Article  CAS  Google Scholar 

  16. Pfaffl MW. A new mathematical model for relative quantification in real-time RT--PCR. Nucleic Acids Res. 2001;29(9):e45.

    Article  CAS  Google Scholar 

  17. Garbarino AD, Saucedo L, Giordana S, Magri ML, Buffone MG, Neuspiller F, Vazquez-Levin MH, Marín-Briggiler CI. Fibroblast growth factor 2 (FGF2) is present in human spermatozoa and is related with sperm motility. The use of recombinant FGF2 to improve motile sperm recovery. Androl US. 2017;5(5):990–8.

    Article  Google Scholar 

  18. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

    Article  Google Scholar 

  19. Hur B, Kang D, Lee S, Moon JH, Lee G, Kim S. Venn-diaNet: venn diagram based network propagation analysis framework for comparing multiple biological experiments. BMC Bioinform. 2019;20(Suppl 23):667.

    Article  Google Scholar 

  20. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545-51.

    Article  CAS  Google Scholar 

  21. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    Article  CAS  Google Scholar 

  22. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  Google Scholar 

  23. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.

    Article  CAS  Google Scholar 

  24. Yu G, Li F, Qin Y, Bo X, Wu Y, Wang S. GOSemSim: an R package for measuring semantic similarity among GO terms and gene products. Bioinformatics. 2010;26(7):976–8.

    Article  CAS  Google Scholar 

  25. Chin CH, Chen SH, Wu HH, Ho CW, Ko MT, Lin CY. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8(Suppl 4):11.

    Article  Google Scholar 

  26. Katz DJ, Teloken P, Shoshany O. Male infertility-the other side of the equation. Aust Fam Phys. 2017;46(9):641–6.

    Google Scholar 

  27. Freitas MJ, Vijayaraghavan S, Fardilha M. Signaling mechanisms in mammalian sperm motility. Biol Reprod. 2017;96(1):2–12.

    Google Scholar 

  28. Cohen G, Rubinstein S, Gur Y, Breitbart H. Crosstalk between protein kinase A and C regulates phospholipase D and F-actin formation during sperm capacitation. Dev Biol. 2004;267(1):230–41.

    Article  CAS  Google Scholar 

  29. Wu K, Mei C, Chen Y, Guo L, Yu Y, Huang D. C-type natriuretic peptide regulates sperm capacitation by the cGMP/PKG signalling pathway via ca(2+) influx and tyrosine phosphorylation. Reprod Biomed Online. 2019;38(3):289–99.

    Article  CAS  Google Scholar 

  30. Koch S, Acebron SP, Herbst J, Hatiboglu G, Niehrs C. Post-transcriptional wnt signaling governs epididymal sperm maturation. Cell. 2015;163(5):1225–36.

    Article  CAS  Google Scholar 

  31. Silva JV, Cabral M, Correia BR, Carvalho P, Sousa M, Oliveira PF, Fardilha M. mTOR signaling pathway regulates sperm quality in older men. Cells Basel. 2019;8(6):629.

    Article  CAS  Google Scholar 

  32. Rago V, De Rose D, Santoro M, Panza S, Malivindi R, Andò S, D’Agata R, Aquila S. Human sperm express the receptor for glucagon-like peptide-1 (GLP-1), which affects sperm function and metabolism. Endocrinology. 2020. https://doi.org/10.1210/endocr/bqaa031.

    Article  Google Scholar 

  33. Ku JY, Park MJ, Park HJ, Park NC, Joo BS. Combination of korean Red Ginseng extract and hydrogen-rich water improves spermatogenesis and sperm motility in male mice. Chin J Integr Med. 2020;26(5):361–9.

    Article  CAS  Google Scholar 

  34. Kurowska P, Barbe A, Różycka M, Chmielińska J, Dupont J, Rak A. Apelin in reproductive physiology and pathology of different species: a critical review. Int J Endocrinol. 2018;2018:9170480.

    Article  Google Scholar 

  35. Lehti MS, Sironen A. Formation and function of sperm tail structures in association with sperm motility defects. Biol Reprod. 2017;97(4):522–36.

    Article  Google Scholar 

  36. Zhang B, Ma H, Khan T, Ma A, Li T, Zhang H, Gao J, Zhou J, Li Y, Yu C, et al: A DNAH17 missense variant causes flagella destabilization and asthenozoospermia. J Exp Med 2020, 217(2).

  37. Chu DS, Shakes DC. Spermatogenesis. Adv Exp Med Biol. 2013;757:171–203.

    Article  CAS  Google Scholar 

  38. Coutton C, Vargas AS, Amiri-Yekta A, Kherraf ZE, Ben MS, Le Tanno P, Wambergue-Legrand C, Karaouzène T, Martinez G, Crouzy S, et al. Mutations in CFAP43 and CFAP44 cause male infertility and flagellum defects in Trypanosoma and human. Nat Commun. 2018;9(1):686.

    Article  Google Scholar 

  39. Ben KM, Coutton C, Zouari R, Karaouzène T, Rendu J, Bidart M, Yassine S, Pierre V, Delaroche J, Hennebicq S, et al. Mutations in DNAH1, which encodes an inner arm heavy chain dynein, lead to male infertility from multiple morphological abnormalities of the sperm flagella. Am J Hum Genet. 2014;94(1):95–104.

    Article  Google Scholar 

  40. Cao X, Cui Y, Zhang X, Lou J, Zhou J, Bei H, Wei R. Proteomic profile of human spermatozoa in healthy and asthenozoospermic individuals. Reprod Biol Endocrinol. 2018;16(1):16.

    Article  Google Scholar 

  41. Foot N, Henshall T, Kumar S. Ubiquitination and the regulation of membrane proteins. Physiol Rev. 2017;97(1):253–81.

    Article  Google Scholar 

  42. Sidjanin DJ, Park AK, Ronchetti A, Martins J, Jackson WT. TBC1D20 mediates autophagy as a key regulator of autophagosome maturation. Autophagy. 2016;12(10):1759–75.

    Article  CAS  Google Scholar 

  43. Chen HY, Sun JM, Zhang Y, Davie JR, Meistrich ML. Ubiquitination of histone H3 in elongating spermatids of rat testes. J Biol Chem. 1998;273(21):13165–9.

    Article  CAS  Google Scholar 

  44. Wei N, Serino G, Deng XW. The COP9 signalosome: more than a protease. Trends Biochem Sci. 2008;33(12):592–600.

    Article  CAS  Google Scholar 

  45. Wang Y, Devereux W, Stewart TM, Casero RJ. Polyamine-modulated factor 1 binds to the human homologue of the 7a subunit of the Arabidopsis COP9 signalosome: implications in gene expression. Biochem J. 2002;366(Pt 1):79–86.

    Article  CAS  Google Scholar 

  46. Sumara I, Quadroni M, Frei C, Olma MH, Sumara G, Ricci R, Peter M. A Cul3-based E3 ligase removes Aurora B from mitotic chromosomes, regulating mitotic progression and completion of cytokinesis in human cells. Dev Cell. 2007;12(6):887–900.

    Article  CAS  Google Scholar 

  47. Jiang J. Regulation of Hh/Gli signaling by dual ubiquitin pathways. Cell Cycle. 2006;5(21):2457–63.

    Article  CAS  Google Scholar 

  48. Wang S, Zheng H, Esaki Y, Kelly F, Yan W. Cullin3 is a KLHL10-interacting protein preferentially expressed during late spermiogenesis. Biol Reprod. 2006;74(1):102–8.

    Article  CAS  Google Scholar 

  49. Kim J, Tsuruta F, Okajima T, Yano S, Sato B, Chiba T. KLHL7 promotes TUT1 ubiquitination associated with nucleolar integrity: implications for retinitis pigmentosa. Biochem Biophys Res Commun. 2017;494(1–2):220–6.

    Article  CAS  Google Scholar 

  50. Zhou Z, Kawabe H, Suzuki A, Shinmyozu K, Saga Y. NEDD4 controls spermatogonial stem cell homeostasis and stress response by regulating messenger ribonucleoprotein complexes. Nat Commun. 2017;8:15662.

    Article  CAS  Google Scholar 

  51. Corral-Vazquez C, Salas-Huetos A, Blanco J, Vidal F, Sarrate Z, Anton E. Sperm microRNA pairs: new perspectives in the search for male fertility biomarkers. Fertil Steril. 2019;112(5):831–41.

    Article  CAS  Google Scholar 

  52. Vojtech L, Woo S, Hughes S, Levy C, Ballweber L, Sauteraud RP, Strobl J, Westerberg K, Gottardo R, Tewari M, et al. Exosomes in human semen carry a distinctive repertoire of small non-coding RNAs with potential regulatory functions. Nucleic Acids Res. 2014;42(11):7290–304.

    Article  CAS  Google Scholar 

  53. Catalanotto C, Cogoni C, Zardo G. MicroRNA in control of gene expression: an overview of nuclear functions. Int J Mol Sci. 2016;17(10):1712.

    Article  Google Scholar 

Download references

Acknowledgements

We would like to thank American Journal Experts (https://www.aje.com) for English language editing.

Funding

No funding.

Author information

Authors and Affiliations

  1. Department of Urology, The Second Affiliated Hospital of Anhui Medical of University, 230601, Hefei, China

    Ci Zou, Shen Xu, Hao Geng, Enlai Li, Wei Sun & Dexin Yu

  2. Department of Urology, The Second Affiliated Hospital of Anhui Medical University, NO.678 Furong Road, 230601, Hefei, China

    Dexin Yu

Authors
  1. Ci Zou
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Shen Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Hao Geng
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Enlai Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Wei Sun
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Dexin Yu
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

DY and CZ conceived the project, designed the study, and drafted the manuscript. CZ and SX directed the study. HG, EL and WS collected the public data. DY revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Dexin Yu.

Ethics declarations

Ethics approval and consent to participate

This study followed the principles of the Helsinki Declaration and was approved by the Research Ethics Committee of the Second Affiliated Hospital of Anhui Medical University (YX2022-046 [F1]), and written informed consent was obtained from all patients.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Supplementary results of gene differential analysis and functional enrichment analysis for asthenozoospermia.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zou, C., Xu, S., Geng, H. et al. Bioinformatics analysis identifies potential hub genes and crucial pathways in the pathogenesis of asthenozoospermia. BMC Med Genomics 15, 252 (2022). https://doi.org/10.1186/s12920-022-01407-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-022-01407-5

Keywords

  • Asthenozoospermia
  • Bioinformatics analysis
  • Hub genes
  • Pathways
  • miRNAs

BMC Medical Genomics

ISSN: 1755-8794

Contact us