Skip to main content
Download PDF

DNA methylation patterns in patients with asthenospermia and oligoasthenospermia

BMC Genomics volume 25, Article number: 602 (2024) Cite this article

Abstract

Background

Spermatogenesis is a highly regulated and complex process in which DNA methylation plays a crucial role. This study aimed to explore the differential methylation profiles in sperm DNA between patients with asthenospermia (AS) and healthy controls (HCs), those with oligoasthenospermia (OAS) and HCs, and patients with AS and those with OAS.

Results

Semen samples and clinical data were collected from five patients with AS, five patients with OAS, and six age-matched HCs. Reduced representation bisulfite sequencing (RRBS) was performed to identify differentially methylated regions (DMRs) in sperm cells among the different types of patients and HCs. A total of 6520, 28,019, and 16,432 DMRs were detected between AS and HC, OAS and HC, and AS and OAS groups, respectively. These DMRs were predominantly located within gene bodies and mapped to 2868, 9296, and 9090 genes in the respective groups. Of note, 12, 9, and 8 DMRs in each group were closely associated with spermatogenesis and male infertility. Furthermore, BDNF, SMARCB1, PIK3CA, and DDX27; RBMX and SPATA17; ASZ1, CDH1, and CHDH were identified as strong differentially methylated candidate genes in each group, respectively. Meanwhile, the GO analysis of DMR-associated genes in the AS vs. HC groups revealed that protein binding, cytoplasm, and transcription (DNA-templated) were the most enriched terms in the biological process (BP), cellular component (CC), and molecular function (MF), respectively. Likewise, in both the OAS vs. HC and AS vs. OAS groups, GO analysis revealed protein binding, nucleus, and transcription (DNA-templated) as the most enriched terms in BP, CC, and MF, respectively. Finally, the KEGG analysis of DMR-annotated genes and these genes at promoters suggested that metabolic pathways were the most significantly associated across all three groups.

Conclusions

The current study results revealed distinctive sperm DNA methylation patterns in the AS vs. HC and OAS vs. HC groups, particularly between patients with AS and those with OAS. The identification of key genes associated with spermatogenesis and male infertility in addition to the differentially enriched metabolic pathways may contribute to uncovering the potential pathogenesis in different types of abnormal sperm parameters.

Peer Review reports

Introduction

Infertility affects approximately 10–20% of couples globally, with male factors contributing to 30–50% of these cases [1]. Asthenospermia (AS) is characterized by a decreased percentage of motile spermatozoa (sperm motility defined as below 32% progressive motility [PR] or below 40% total motility), whereas oligoasthenospermia (OAS) is characterized by a low concentration of spermatozoa and low percentage of motile spermatozoa (sperm concentration defined as below 15 million sperm per milliliter [or below 39 million sperm in the ejaculate] and sperm motility defined as below 32% PR or below 40% total motility), both of which fall below the 5th World Health Organization (WHO) manual for examination of human semen [1, 2]. Although both conditions represent the principal manifestations of impaired spermatogenesis, their exact etiology remains unclear.

Spermatogenesis is a highly regulated and complex process that is influenced by various genetic and environmental factors [3, 4]. Increasing evidence suggests that increased prevalence of genetic abnormalities was associated with decreased sperm concentration and motility [5,6,7]. In the past decades, several novel gene mutations have been discovered in AS and OAS, such as DNAH9, KCNU1, MNS1, GSTM1 and T1, etc [8,9,10,11,12]. Epigenetics, another layer of gene regulation, could strongly affect gene expression and function in both sperm function and fertilizing ability [5, 13]. Furthermore, epigenetic regulation was also changed by environmental factors [5]; the dynamic epigenome would undergo substantial changes over years, affecting disease susceptibility and determining the resulting male infertility phenotype [5, 14]. Therefore, it is important to identify epigenome-based biomarkers to assess spermatogenesis capacity in unexplained infertility, complementary to the traditional genetic biomarkers. However, differences at the genetic and epigenetic levels between the different types of impaired spermatogenesis, particularly AS and OAS, remain poorly understood.

Several data have shown that epigenetic alterations, particularly DNA methylation, play a significant role in spermatogenesis, might resulting in decreased sperm number and motility [4, 5, 15]; for instance, increased methylation of the DAZL promoter CpG island has been linked to defective human sperm in oligoasthenoteratozoospermic (OAT) men [16]. More recent study has also reported that AS patients exhibited significant hypermethylation at 2 CpG sites of IGF-2, and significant decreases in methylation at one CpG site of KCNQ1, along with three CpG sites of MEST compared to healthy controls (HCs) [17]. In OAS patients, DYDC2 and NXF3 were hypermethylated [18]. However, to date, no study has assessed epigenetic differences between AS and OAS and only a few studies have specifically investigated the differences between patients with OAS and HCs.

Previous studies have primarily focused on the aberrant methylation signature of imprinted and nonimprinted genes using the 450 K BeadChip or the bisulfite sequencing polymerase chain reaction (PCR); these studies reported the abnormally methylated genes, such as PTPRN2, RPH3AL, and APCS, etc. in oligo-/oligoastheno-zoospermic, asthenozoospermic, teratozoospermic and asthenoteratozoospermic patients [19,20,21]. Despite the advancements made in these studies, the complex correlation between epigenetics and sperm abnormalities remains to be comprehensively illustrated. Reduced representation bisulfite sequencing (RRBS), a cost-effective and efficient method for genome-wide methylation profiling, can offer enhanced coverage and depth for analyzing epigenetics [22].

In the present study, RRBS was performed to comprehensively investigate the global DNA methylation profiles in three distinct groups (AS vs. HCs, OAS vs. HCs, and AS vs. OAS). The study aimed to elucidate the underlying pathogenesis of spermatogenesis in AS and OAS patients, and to provide novel insights into the distinct characteristics of different spermatogenesis abnormalities.

Methods

Participants

Five patients each with AS and OAS and six age-matched HCs were recruited from Peking Union Medical College Hospital (PUMCH) in April 2021. The inclusion criteria for AS participants were as follows: patients with no chromosomal abnormalities, varicocele, congenital hypogonadotropic hypogonadism, chronic illnesses, history of genitourinary infections impacting fertility, or exposure to heavy metals and those who did not smoke. Any diagnosed or suspected factors that might affect reproductive health were excluded through a thorough medical history conducted by experienced urologists. The same inclusion criteria were followed for patients with OAS. These patients we enrolled had no history of hormone substance abuse or medication affecting hormone balance, and no recent acute infections. The inclusion criteria for the HC group involved no history of AS or OAS, absence of diseases known to cause male infertility, and recent childbirth. Moreover, all individuals in the HC group exhibited normal sperm morphology and sperm parameters (Table 1). The included patients with AS and OAS received the same basic traditional medical treatment before collection; however, they showed poor response to the treatment. The PUMCH ethics committee approved this study (No. JS-3112).

Sample collection and sperm cell isolation

Semen samples were collected by masturbation after sexual abstinence for a period of 2–7 days, stored in sterile containers, and the semen analysis were conducted as per the 5th WHO guidelines [2]. The sperm cell isolation was carried out using density gradient centrifugation following the 5th WHO standard. Percoll (Sigma-Aldrich, Dorset, UK) and ×10 Human Tubal Fluid (HTF) were used to prepare an iso-osmotic gradient [23]. Briefly, discontinuous double-density gradients with two layers (2 mL 40% and 2 mL 80%) were prepared. Then, 1 mL of semen mixed with the density gradient fluid was centrifuged for 20 min at 300 g. After removing the supernatan, the collected sperm cells were washed with 5 mL of x1 HTF and centrifuged for 6 min at 200 ×g. Then this procedure was repeated once more. Finally, the supernatant sperm was stored at − 80 °C for DNA extraction.

DNA extraction and library construction

We performed DNA isolation from frozen sperm cells using a kit (FineMag Universal Genomic DNA Extraction Kit by Magnetic Bead Method, Genfine, Beijing, China, M202). Spermatozoa precipitates with a cell count between 1 × 106 and 107 were resuspended in 300 μL of Buffer MDA and 20 μL of Proteinase K. Next, 300 μL of Buffer GHL was added, and the mixture was inverted 3–5 times for complete dissolution. Then, 200–500 μL of the treated sample was processed with magnetic beads to isolate high-purity DNA. The quality control of the isolated DNA was performed to ensure the purity and concentration. Finally, spectrophotometry and fluorescence were used to quantify the isolated DNA samples.

For DNA library construction, all DNA samples were subjected to RRBS library preparation using the Acegen Rapid RRBS Library Prep Kit (Acegen, Cat. No. AG0422) according to the manufacturer’s protocol. In brief, 100 ng of input DNA was digested with MspI. The DNA fragments underwent end repair, A-tailing, and cytosine ligation. Fragments ranging from 150 to 300 base pairs were gel-extracted and subjected to bisulfite treatment using the EZ DNA Methylation Gold Kit (Zymo Research). After bisulfite treatment, these fragments were then PCR-amplified to construct the final DNA libraries. The library quality was assessed using Qubit 2.0, Agilent 2100, and q-PCR for estimating insert size and concentration.

RRBS and data processing

After comprehensive DNA library analysis, the fragments were sequenced using the Illumina NovaSeq 6000 system (Illumina, Inc., San Diego, CA, USA), with the resulting raw data being stored in the FASTQ format. Base calling was performed using Illumina Casava 1.8, and data quality was assessed and refined using FastQC and Trimmomatic. In brief, we first used FastQC v0.11.7 to perform basic statistics on the raw data. Then, the raw data was processed with the software Trimmomatic [24] v0.36. For methylation analysis, Bsmap [25] v 2.7.3 was used to aligne the clean data to the reference genome downloaded from genome website directly. After mapping to genome and generating SAM file, SAMTOOLS [26] v1.9 was used to convert SAM file to BAM file, sort and index for subsequent analysis. Then the measurement and reliability assessment of differentially methylated sites were performed using the Bsmap v 2.7.3 software. Finally, the Metilene [27] v0.2-8 software identified DMRs and calculated the average methylation levels. P-value was adjusted using the false discovery rate (FDR). DMRs meeting the criteria of FDR < 0.05 were considered significant differences.

Function enrichment analysis and statistical analysis

IBM SPSS Statistics version 27.0 (IBM, Armonk, NY, USA) was used for statistical analyses. The protein–protein interaction (PPI) network was constructed using the Search Tool for the Retrieval of Interacting Genes (STRING) database (https://string-db.org/) [28], with a threshold score of 0.150. The STRING database systematically collects and integrates protein–protein interactions, encompassing both physical interactions and functional associations. Its primary aim is to provide approximate confidence levels for the associations between proteins [28]. Thereafter, Cytoscape software and its plug-in Cyto-Hubba were utilized to attain a higher degree of connectivity, in order to identify the hub gene based on maximal clique centrality (MCC) [29]. MCC is a precise predictive method for ranking essential genes or gene products, with better performance compared to other methods [29, 30]. Then we selected the top 10 node genes ranked by MCC in CytoHubba analysis. These nodes are considered to be the most important in the interaction network.

Gene Ontology (GO, http://www.geneontology.org/) [31] and Kyoto Encyclopedia of Genes and Genomes (KEGG) [32] analyses were performed using DAVID [33] (P value < 0.05 was considered as significant enrichment) to provide potential biological functions and pathways.

Results

Clinical characteristics of participants

The average age of participants in the AS, OAS, and HC groups was 33.75 ± 3.94, 30.75 ± 1.80, and 34.00 ± 1.08 years, respectively. All participants were of Han Chinese ethnicity. In the AS and OAS groups, even after administering treatment, the sperm count, sperm concentration, and sperm motility of patients in these groups remained lower than those of HCs. The clinical characteristics of the five patients with AS, five with OAS, and six HCs included in this study are shown in Table 1.

Table 1 Clinical characteristics of the participants
Full size table

Global differential DNA methylation profiling

To elucidate the global DNA methylation differences across various types of abnormal sperm parameters, RRBS was performed for DNA methylation profiling between three distinct pairs: AS and HCs, OAS and HCs, and AS and OAS. A total of 6520, 28,019, and 16,432 DMRs were identified between each group, respectively. Of these DMRs, the number of hypermethylated and hypomethylated DMRs was 4251 and 2269, 25,747 and 2272, and 3179 and 23,253, respectively. According to Fig. 1A, B and C, the average DNA methylation between each group was significantly different. Furthermore, the genomic zone distribution analysis of these DMRs indicated a consistent pattern across all groups, with gene bodies (44%, 39%, and 40%) being the most frequently affected region, followed by introns (40%, 36%, and 36%) and CGI shores (12%, 21%, and 22%) (Fig. 1D). The distribution of these DMRs in chromosomes is shown in Fig. 1E, which illustrates the distribution of hyper-DMRs and hypo-DMRs.

Fig. 1
figure 1

Global differential DNA methylation profiling. A–C DMR methylation level distribution and clustering heat map of AS vs. HC, OAS vs. HC and AS vs. OAS group. D DMR anchoring region distribution. E Circos plots of the significance of DMRs on specific chromosomes. From the outside to the inside, the distribution of hyper-DMR is shown as a scatter plot with the direction toward the outside, i.e., higher significance = more dots toward the outside; the GC content of sequences is shown as a heatmap, darker color = higher GC content; the gene density of sequences is shown as a heatmap, darker color = higher GC content; and the distribution of hypo-DMR is shown as a scatter plot with the direction toward the inside, i.e., higher significance = more dots toward the outside

Full size image

For AS and HC, OAS and HC, and AS and OAS groups, the 6520, 28,019 and 16,432 identified DMRs were mapped to 2868, 9296, and 9090 genes, respectively. The top 100 DMR-associated genes for each group are shown in additional file 1: Tables S1, S2, and S3. Among these, 56 genes were hypermethylated and 9 were hypomethylated in the AS group compared with the HC group. In the OAS and HC groups, all 58 annotated genes were hypermethylated, whereas in the AS and OAS groups, 44 genes were hypomethylated. In summary, these results suggest that the top 100 genes associated with DMRs show consistent patterns of methylation levels, both high and low, which mirror the overall methylation profiles observed across the entire genome.

Integrative analysis and PPI network construction of DMR-associated genes

Genes that play pivotal roles in spermatogenesis, male infertility, asthenozoospermia, and oligoasthenospermia were comprehensively assessed using the National Center for Biotechnology Information (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/) [34], GeneCards (https://www.genecards.org/) [35], the Online Mendelian Inheritance in Man (https://www.omim.org/) database [36], and an extensive review of current literature. The three groups included 12, 8, and 9 differential methylated genes, respectively, among the DMRs (Table 2). These genes were hypermethylated, hypermethylated, and hypomethylated in the three groups, respectively. Next, the node genes in each group were identified. In the AS vs. HC group, BDNF, SMARCB1, and PIK3A emerged as the most prominently ranked genes, whereas DDX27, RBMX, and SPATA17 and ASZ1, CDH1, and CHDH held the highest rankings in the OAS vs. HC and AS vs. OAS groups, respectively (Fig. 2A-C). Overall, the findings indicated that these differential methylated genes may be promising candidates for future research on AS and OAS.

Table 2 DMR-harboring genes involved in spermatogenesis or male infertility or asthenospermia or oligoasthenospermia in the AS vs. HC, OAS vs. HC, and AS vs. OAS groups
Full size table
Fig. 2
figure 2

Identification of key hub genes using Cytoscape’s CytoHubba. A-C Top 10 node genes in each group based on maximum clique centrality

Full size image

GO and KEGG analysis of DMR-associated genes

The above 2868, 9296, and 9090 genes identified from the AS vs. HC, OAS vs. HC, and AS vs. OAS groups, respectively, were subjected to GO and KEGG analyses. For the AS and HC group, GO displayed protein binding as the most enriched pathway in the biological process (BP), with cytoplasm the most enriched in the cellular component (CC) and transcription (DNA-templated) the most enriched in the molecular function (MF). Of note, similar results were obtained for both the OAS vs. HC and AS vs. OAS groups, with protein binding, nucleus, and transcription (DNA-templated) being the most enriched pathways in BP, CC, and MF, respectively (Fig. 3D-F)

Fig. 3
figure 3

GO and KEGG annotation of DMR-associated genes. A-C GO analysis. D-E KEGG analysis

Full size image

GO and KEGG analysis of DMR-associated genes in promoters

The DMR-associated genes among the AS and HC, OAS and HC, and AS and OAS groups were subjected to GO and KEGG analysis. GO results were consistent across the three groups, with protein binding emerging as the most enriched pathway in BP, nucleus in CC, and transcription (DNA-templated) in MF (Fig. 4A-C). Likewise, Fig. 4D-F shows that the metabolic pathway was the most affected pathway for all groups, in line with the results of the whole-genome analysis.

Fig. 4
figure 4

GO and KEGG annotation of DMR-associated genes in promoters. A-C GO analysis. D-E KEGG analysis

Full size image

Discussion

In our data, all HCs exhibited normal semen parameters and sperm morphology. This consistency underscores the reliability of our findings. Patients with AS and OAS demonstrated impaired sperm motility (and number), further indicating the presence of underlying mechanisms contributing to reproductive dysfunction in these patients. A comprehensive analysis also unveiled 6520, 28,019, and 16,432 DMRs differentiating patients with AS from HCs, patients with OAS from HCs, and patients with AS from those with OAS, respectively. Of note, the identified DMRs exhibited a prevalent hypermethylation pattern in both patients with AS and OAS compared with HCs. These findings align with a previous report on oligo-/oligoasthenozoospermic individuals, wherein 1436/1680 DMRs were hypermethylated [19]. The present study also revealed that these hypermethylated DMRs were primarily located within gene bodies, in line with earlier research on male infertility and sperm methylation [37,38,39]. Taken together, these results indicated that the DNA methylation pattern of spermatozoa may be distinct compared with that of somatic cells, challenging the conventional belief that methylation is principally enriched in the promoter region. The underlying mechanism, however, remains unknown, potentially attributed to significant differences in biological function between germ cells and somatic cells [40]. It should also be noted that sperm may undergo a complex process of epigenetic reprogramming distinct from that in somatic cells, encompassing the erasure and reconstruction of sperm epigenetic profiles [41, 42]. This intricate reprogramming may involve alterations in the dynamic equilibrium between methylation and demethylation processes throughout spermatogenesis and development, ultimately leading to distinct methylation patterns. Thus, further research focusing on DMR-associated genes located in gene bodies is needed. Of note, when patients with OAS and AS were compared, the identified DMRs contained increased hypermethylation in the former group. This indicated that a subset of samples with poorer semen parameters may harbor increased methylated alterations. A similar study reported that decreased sperm count, known as oligospermia, may be linked to modified sperm DNA methylation patterns across both imprinted and nonimprinted genes [43]. However, conflicting findings have emerged from studies investigating the correlation between DNA methylation and sperm counts. Some studies indicated no significant association, whereas others reported only select regional alterations to be associated with sperm motility [44, 45]. This discrepancy might be due to the significant heterogeneity between and within samples and owing to the small sample size analyzed in the current study.

GO analysis revealed that the identified DMR-associated genes are linked to “protein binding” in BP, in the “cytoplasm” and “nucleus” in CC, and in “transcription (DNA-templated)” in MF across the three experimental groups. An accumulation of cytoplasm has the potential to disrupt the balance of intracellular molecules, thereby triggering an elevation in the level of reactive oxygen species (ROS) [46]. Of interest, excessive levels of ROS have been reported to lead to genome-wide hypomethylation [47, 48]. They are positively correlated with IGF2 and MLH1 hypermethylation and negatively correlated with H19 hypomethylation in idiopathic male infertility, all of which are associated with semen quality and spermatogenesis interruption [48]. Furthermore, oxidative stress and DNA methylation share a common denominator: the one-carbon cycle, suggesting a significant correlation between oxidative stress, methylation processes, and epigenesis [49]. The identification of ROS-induced epigenetic damage warrants further investigation. In addition, changes in cytoplasm in the human sperm nucleus are associated with male infertility and sperm abnormalities [50]. Furthermore, the nucleus and mitochondria interact intricately, with mitochondria playing diverse roles in epigenetic regulation, gametogenesis, and reproduction [51]. In summary, these findings may explain the results observed in the GO analysis of patients with AS and OAS. Moreover, the epigenetic alterations in genes involved in both cytoplasmic and nuclear components may signify changes in the intracellular milieu in these patients. Nevertheless, this hypothesis remains speculative, and the relationship among cytoplasmic and nuclear genes and ROS levels, mitochondrial function, and sperm quality needs further study. In addition, genes participating in DNA transcription that play important roles in regulating gene expression and function may be implicated in the critical regulation of spermatogenesis. This coincides with the finding of “transcription (DNA-templated)” in the MF GO terms. Meanwhile, KEGG analysis highlighted that DMR-associated genes were associated with “metabolic pathways” in all three groups. Studies have reported metabolic alterations in the seminal plasma of patients with AS and OAS [52,53,54,55]. In line, the current study data indicated that distinct metabolic alterations may be involved in the pathology of AS and OAS. However, despite the centrality of energy metabolism in the underlying etiology of AS, few studies have identified the associated metabolic pathways. Thus, the most affected metabolic pathways associated with these DMR-annotated genes in both AS and OAS need to be explored.

A total of 12, 9, and 8 differentially methylated genes were identified in AS and HC, OAS and HC, and AS and OAS groups, respectively. These genes are closely associated with spermatogenesis and male infertility. Among these, BDNF, SMARCB1, PIK3A, and DDX27; RBMX and SPATA17; ASZ1, CDH1, and CHDH were identified as strong differentially methylated candidate genes in the three pair groups, respectively. BDNF, a member of the nerve growth factor family located in the human sperm, has been reported to participate in mitigating oxidative stress, which is important for maintaining sperm motility and facilitating sperm–egg. Therefore, alterations in BDNF may be associated with an impaired capacity of sperm to defend against oxidative stress, contributing to asthenospermia and related disorders [56, 57]. Of note, a recent investigation conducted in zebrafish models revealed that paternal exposure to microcystin–leucine arginine, a substance known to adversely affect male reproduction and disrupt the developmental trajectory of the offspring, resulted in a notable increase in the DNA methylation of the BDNF promoter within sperm, indicating that hypermethylated BDNF in AS observed in the present study might explain a potential connection with male reproductive outcomes [58]. SMARCB1, known as Snf5, is a key constituent of the SWI/SNF complex that regulates diverse cellular processes in model eukaryotes. The SWI/SNF complexes have been illustrated to potentially collaborate with certain specific transcription factors, thereby assisting in the activation of male-specific germline genes [59]. PIK3CA, a subtype of the regulatory subunit of phosphoinositide 3-kinase (PI3K), phosphorylates phosphatidylinositol and participates in the PI3K-Akt signaling pathway, the activation of which enhances the proliferation and confers anti-apoptotic properties to Sertoli cells [60]. In addition, the PI3K/Akt/mTOR signaling pathway is involved in the regulation and induction of hypoxia-inducible factor 1 alpha (HIF1α) expression; the overexpression of HIF1α has been linked to detrimental effects on testicular spermatogenesis [61]. Therefore, the identified hypermethylated genes in AS likely play a complex role in spermatogenic cells.

Compared with HCs, RBMX was hypermethylated in OAS. RBMX serves as the primary source for several autosomal retrotransposed genes, one of which is the human HNRNPGT gene (encoding the hnRNP G-T protein). The haploinsufficiency of HNRNPGT resulted in abnormal sperm production in mice [62]. In addition, a recent case report revealed that deleterious variants within the RBMXL2 gene have the potential to induce male infertility by contributing to a complete meiotic arrest, potentially by disrupting the normal trajectory of gene expression [63]. However, further research is needed to understand its role in spermatogenesis.

Between the AS and OAS groups, studies have shown that spermatogenesis-associated protein 17 (SPATA17) promoted apoptosis in the GC-1 cell line and that its overexpression mediated apoptosis in transgenic male mouse germ cells [64,65,66,67]. Regarding ASZ1, research on mice revealed that nuclear respiratory factor (NRF)-1, a key transcription factor that regulates mitochondrial biogenesis in germ cell development, exhibited a precise crosstalk with epigenetic modulation. This regulatory crosstalk directly influenced the expression of various germ cell-specific genes, particularly ASZ1, thereby interfering with germ cell proliferation and inducing male sterility [68]. ASZ1 also exhibited the ability to suppress retrotransposons and retrotransposon-driven aberrant chimeric transcripts within adenocytes via DNA methylation in the piRNA pathway, which is critical for the integrity of spermatogenesis [69]. Thus, the NRF-1–ASZ1–piRNA axis may be a promising regulatory pathway in spermatogenesis. In bovine studies, sperm-specific HMRs have been identified to target various spermatogenesis-related genes, including ASZ1 [70]. Meanwhile, CDH1 is considered a definitive marker delineating undifferentiated spermatogonia. It plays an indispensable role in orchestrating gonadotropin production within the testes through its facilitation of the migration and aggregation of primordial germ cells and somatic cells and establishing a critical regulatory nexus in the intricate processes of spermatogenesis [71]. CHDH deficiency in mice was shown to reduce sperm motility [72]. Moreover, the CHDH polymorphism rs12676 in males was associated with reduced ATP concentration in sperm [73].

Taken together, these studies highlight the roles of the identified genes in spermatogenesis and male infertility. They could act as promising candidates for in-depth studies to uncover the mechanisms of AS and OAS and explore their association with alterations in DNA methylation. However, it is not clear why these genes exhibit a more pronounced hypermethylation pattern in patients with OAS compared to those with AS. Thus, further studies are needed to explore the association between the DNA methylation of DMRs and gene expression.

There are several limitations in this study. First, the sample size was small, with no an independent validation cohort, thus the sampling bias was inevitable. While our study identified certain differentially methylated genes associated with spermatogenesis, and male infertility, the generalizability of these findings might be limited due to the restricted sample size. Our study inclined towards identifying trends in change rather than drawing final conclusions. Additionally, the small sample size limits subgroup analyses to potentially reveal different spermatogenic dysfunction with different demographic or clinical characteristics. Future studies with larger sample size are warranted to validate our findings. Second, this study lacked oligospermia controls. Therefore, the specific impact of decreased sperm count on sperm methylation status could not be clarified. Third, genes associated with DMRs located at gene bodies were not investigated. Future research using deep sequencing is essential to validate the candidate genes identified in the current study. Last, we utilized frozen sperm samples for DNA extraction. We acknowledge that frozen sperm may exhibit lower DNA integrity, viability, and motility values compared to fresh sperm [74, 75]; however, it is challenging to preserve and manage fresh sperm over the long term; and there is limited information available regarding the effects of cryopreservation on epigenetic modulation in sperm, and existing studies present contradictory results [75]. In the future studies, we intend to investigate the impact of the cryopreservation process on sperm DNA methylation.

Conclusions

The current study findings revealed distinct sperm DNA methylation patterns among patients with AS and HCs, patients with OAS and HCs, and patients with AS and those with OAS. These DMRs could be the promising biomarkers of impaired spermatogenesis. The key DMR-associated genes identified in the study were associated with spermatogenesis and male infertility, suggesting their involvement in the regulation of sperm motility and quantity via aberrant DNA methylation alterations. These DMR-associated genes were enriched in metabolic pathways, shedding light on the different metabolism-related pathogenesis in different dys-spermatism types. These findings collectively provide novel sights into the impact of DNA methylation on sperm motility, sperm number, and male infertility. The validation of these candidate genes and differentially methylated loci using deep sequencing in large populations is necessary to explore their effect on sperm motility and number.

Data availability

The methylation sequencing data of healthy controls is available from NCBI and have been assigned BioProject accession PRJNA1117580(https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/?term=PRJNA1117580). The methylation sequencing data of patients with asthenospermia and oligoasthenospermia is available from NCBI and have been assigned BioProject accession PRJNA1116754(https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/?term=PRJNA1116754).

Abbreviations

AS:

Asthenospermia

BP:

Biological process

CC:

Cellular component

DMRs:

Differentially methylated regions

GO:

Gene Ontology

HCs:

Healthy controls

KEGG:

Kyoto Encyclopedia of Genes and Genomes

MF:

Molecular function

NP:

Non-progressive motility

OAS:

Oligoasthenospermia

PPI:

Protein–protein interaction

PR:

Progressive motility

RRBS:

Reduced representation bisulfite sequencing

References

  1. Eisenberg ML, Esteves SC, Lamb DJ, Hotaling JM, Giwercman A, Hwang K, Cheng YS. Male infertility. Nat Reviews Disease Primers. 2023;9(1):49.

    Article  PubMed  Google Scholar 

  2. Cooper TG, Noonan E, von Eckardstein S, Auger J, Baker HW, Behre HM, Haugen TB, Kruger T, Wang C, Mbizvo MT, et al. World Health Organization reference values for human semen characteristics. Hum Reprod Update. 2010;16(3):231–45.

    Article  PubMed  Google Scholar 

  3. Sharma A, Minhas S, Dhillo WS, Jayasena CN. Male infertility due to testicular disorders. J Clin Endocrinol Metab. 2021;106(2):e442–59.

    Article  PubMed  Google Scholar 

  4. Du L, Chen W, Cheng Z, Wu S, He J, Han L, He Z, Qin W. Novel gene regulation in normal and abnormal spermatogenesis. Cells. 2021;10(3).

  5. Gunes S, Esteves SC. Role of genetics and epigenetics in male infertility. Andrologia. 2021;53(1):e13586.

    Article  PubMed  Google Scholar 

  6. Li WN, Zhu L, Jia MM, Yin SL, Lu GX, Liu G. Missense mutation in DNAJB13 gene correlated with male fertility in asthenozoospermia. Andrology. 2020;8(2):299–306.

    Article  CAS  PubMed  Google Scholar 

  7. Cavarocchi E, Whitfield M, Saez F, Touré A. Sperm ion transporters and channels in human asthenozoospermia: genetic etiology, lessons from animal models, and clinical perspectives. International journal of molecular sciences. 2022;23(7).

  8. Li Q, Shen Y, Zhao LJ, Wang JB, Huang X. Mutations in CFTR genes are associated with oligoasthenospermia in infertile men undergoing IVF. Andrologia. 2022;54(3):e14355.

    Article  CAS  PubMed  Google Scholar 

  9. Zhang H, He J, Zhao Y, Wu Q, Zou T, Sun J, Zhu H, Wang X, Sun F, Xing J, et al. Effect of glutathione S-transferase gene polymorphisms on semen quality in patients with idiopathic male infertility. J Int Med Res. 2021;49(12):3000605211061045.

    Article  CAS  PubMed  Google Scholar 

  10. Li Y, Wang WL, Tu CF, Meng LL, Hu TY, Du J, Lin G, Nie HC, Tan YQ. A novel homozygous frameshift mutation in MNS1 associated with severe oligoasthenoteratozoospermia in humans. Asian J Androl. 2021;23(2):197–204.

    Article  CAS  PubMed  Google Scholar 

  11. Liu R, Yan Z, Fan Y, Qu R, Chen B, Li B, Wu L, Wu H, Mu J, Zhao L, et al. Bi-allelic variants in KCNU1 cause impaired acrosome reactions and male infertility. Hum Reprod (Oxford England). 2022;37(7):1394–405.

    Article  CAS  Google Scholar 

  12. Tang D, Sha Y, Gao Y, Zhang J, Cheng H, Zhang J, Ni X, Wang C, Xu C, Geng H, et al. Novel variants in DNAH9 lead to nonsyndromic severe asthenozoospermia. Reproductive Biology Endocrinology: RB&E. 2021;19(1):27.

    Article  CAS  Google Scholar 

  13. Lande-Diner L, Zhang J, Ben-Porath I, Amariglio N, Keshet I, Hecht M, Azuara V, Fisher AG, Rechavi G, Cedar H. Role of DNA methylation in stable gene repression. J Biol Chem. 2007;282(16):12194–200.

    Article  CAS  PubMed  Google Scholar 

  14. Sarkar S, Sujit KM, Singh V, Pandey R, Trivedi S, Singh K, Gupta G, Rajender S. Array-based DNA methylation profiling reveals peripheral blood differential methylation in male infertility. Fertil Steril. 2019;112(1):61–e721.

    Article  CAS  PubMed  Google Scholar 

  15. Rotondo JC, Lanzillotti C, Mazziotta C, Tognon M, Martini F. Epigenetics of male infertility: the role of DNA methylation. Front cell Dev Biology. 2021;9:689624.

    Article  Google Scholar 

  16. Navarro-Costa P, Nogueira P, Carvalho M, Leal F, Cordeiro I, Calhaz-Jorge C, Gonçalves J, Plancha CE. Incorrect DNA methylation of the DAZL promoter CpG island associates with defective human sperm. Hum Reprod (Oxford England). 2010;25(10):2647–54.

    Article  CAS  Google Scholar 

  17. Song B, Chen Y, Wang C, Li G, Wei Z, He X, Cao Y. Poor semen parameters are associated with abnormal methylation of imprinted genes in sperm DNA. Reproductive Biology Endocrinology: RB&E. 2022;20(1):155.

    Article  CAS  Google Scholar 

  18. Ma RH, Zhang ZG, Zhang YT, Jian SY, Li BY. Detection of aberrant DNA methylation patterns in sperm of male recurrent spontaneous abortion patients. Zygote (Cambridge England). 2023;31(2):163–72.

    Article  CAS  PubMed  Google Scholar 

  19. Sujit KM, Sarkar S, Singh V, Pandey R, Agrawal NK, Trivedi S, Singh K, Gupta G, Rajender S. Genome-wide differential methylation analyses identifies methylation signatures of male infertility. Hum Reprod (Oxford England). 2018;33(12):2256–67.

    Article  CAS  Google Scholar 

  20. Camprubí C, Salas-Huetos A, Aiese-Cigliano R, Godo A, Pons MC, Castellano G, Grossmann M, Sanseverino W, Martin-Subero JI, Garrido N, et al. Spermatozoa from infertile patients exhibit differences of DNA methylation associated with spermatogenesis-related processes: an array-based analysis. Reprod Biomed Online. 2016;33(6):709–19.

    Article  PubMed  Google Scholar 

  21. Tang Q, Pan F, Yang J, Fu Z, Lu Y, Wu X, Han X, Chen M, Lu C, Xia Y, et al. Idiopathic male infertility is strongly associated with aberrant DNA methylation of imprinted loci in sperm: a case-control study. Clin Epigenetics. 2018;10(1):134.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Beck D, Ben Maamar M, Skinner MK. Genome-wide CpG density and DNA methylation analysis method (MeDIP, RRBS, and WGBS) comparisons. Epigenetics. 2022;17(5):518–30.

    Article  PubMed  Google Scholar 

  23. Gooch J, Tungsirisurp S, Costanzo H, Napier R, Frascione N. Generating aptamers towards human sperm cells using massively parallel sequencing. Anal Bioanal Chem. 2021;413(23):5821–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for illumina sequence data. Bioinf (Oxford England). 2014;30(15):2114–20.

    CAS  Google Scholar 

  25. Xi Y, Li W. BSMAP: whole genome bisulfite sequence MAPping program. BMC Bioinformatics. 2009;10:232.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The sequence alignment/map format and sAMtools. Bioinf (Oxford England). 2009;25(16):2078–9.

    Google Scholar 

  27. Jühling F, Kretzmer H, Bernhart SH, Otto C, Stadler PF, Hoffmann S. Metilene: fast and sensitive calling of differentially methylated regions from bisulfite sequencing data. Genome Res. 2016;26(2):256–62.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.

    Article  CAS  PubMed  Google Scholar 

  29. 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):S11.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36(Database issue):D480–4.

    CAS  PubMed  Google Scholar 

  33. Dennis G Jr., Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4(5):P3.

    Article  PubMed  Google Scholar 

  34. Sayers EW, Bolton EE, Brister JR, Canese K, Chan J, Comeau DC, Farrell CM, Feldgarden M, Fine AM, Funk K, et al. Database resources of the National Center for Biotechnology Information in 2023. Nucleic Acids Res. 2023;51(D1):D29–38.

    Article  CAS  PubMed  Google Scholar 

  35. Stelzer G, Rosen N, Plaschkes I, Zimmerman S, Twik M, Fishilevich S, Stein TI, Nudel R, Lieder I, Mazor Y et al. The GeneCards suite: from gene data mining to disease genome sequence analyses. Current protocols in bioinformatics. 2016;54:1.30.1-1.3.

  36. Amberger JS, Bocchini CA, Scott AF, Hamosh A. OMIM.org: leveraging knowledge across phenotype-gene relationships. Nucleic Acids Res. 2019;47(D1):D1038–43.

    Article  CAS  PubMed  Google Scholar 

  37. He W, Sun Υ, Zhang S, Feng X, Xu M, Dai J, Ni X, Wang X, Wu Q. Profiling the DNA methylation patterns of imprinted genes in abnormal semen samples by next-generation bisulfite sequencing. J Assist Reprod Genet. 2020;37(9):2211–21.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Du Y, Li M, Chen J, Duan Y, Wang X, Qiu Y, Cai Z, Gui Y, Jiang H. Promoter targeted bisulfite sequencing reveals DNA methylation profiles associated with low sperm motility in asthenozoospermia. Hum Reprod (Oxford England). 2016;31(1):24–33.

    Article  CAS  Google Scholar 

  39. Molaro A, Hodges E, Fang F, Song Q, McCombie WR, Hannon GJ, Smith AD. Sperm methylation profiles reveal features of epigenetic inheritance and evolution in primates. Cell. 2011;146(6):1029–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Bajusz I, Henry S, Sutus E, Kovács G, Pirity MK. Evolving role of RING1 and YY1 binding protein in the regulation of germ-cell-specific transcription. Genes. 2019;10(11).

  41. McSwiggin HM, O’Doherty AM. Epigenetic reprogramming during spermatogenesis and male factor infertility. Reprod (Cambridge England). 2018;156(2):R9–21.

    Article  CAS  Google Scholar 

  42. Hernandez Mora JR, Buhigas C, Clark S, Del Gallego Bonilla R, Daskeviciute D, Monteagudo-Sánchez A, Poo-Llanillo ME, Medrano JV, Simón C, Meseguer M, et al. Single-cell multi-omic analysis profiles defective genome activation and epigenetic reprogramming associated with human pre-implantation embryo arrest. Cell Rep. 2023;42(2):112100.

    Article  CAS  PubMed  Google Scholar 

  43. Pacheco SE, Houseman EA, Christensen BC, Marsit CJ, Kelsey KT, Sigman M, Boekelheide K. Integrative DNA methylation and gene expression analyses identify DNA packaging and epigenetic regulatory genes associated with low motility sperm. PLoS ONE. 2011;6(6):e20280.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Xu A, Hua Y, Zhang J, Chen W, Zhao K, Xi W, Wang H, Fang J, Su S, Tang M, et al. Abnormal hypermethylation of the VDAC2 promoter is a potential cause of idiopathic asthenospermia in men. Sci Rep. 2016;6:37836.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Botezatu A, Socolov R, Socolov D, Iancu IV, Anton G. Methylation pattern of methylene tetrahydrofolate reductase and small nuclear ribonucleoprotein polypeptide N promoters in oligoasthenospermia: a case-control study. Reprod Biomed Online. 2014;28(2):225–31.

    Article  CAS  PubMed  Google Scholar 

  46. Rengan AK, Agarwal A, van der Linde M, du Plessis SS. An investigation of excess residual cytoplasm in human spermatozoa and its distinction from the cytoplasmic droplet. Reproductive biology and endocrinology: RB&E. 2012;10:92.

  47. Chianese R, Pierantoni R. Mitochondrial reactive oxygen species (ROS) production alters sperm quality. Antioxid (Basel Switzerland). 2021;10(1).

  48. Darbandi M, Darbandi S, Agarwal A, Baskaran S, Dutta S, Sengupta P, Khorram Khorshid HR, Esteves S, Gilany K, Hedayati M, et al. Reactive oxygen species-induced alterations in H19-Igf2 methylation patterns, seminal plasma metabolites, and semen quality. J Assist Reprod Genet. 2019;36(2):241–53.

    Article  PubMed  Google Scholar 

  49. Menezo YJ, Silvestris E, Dale B, Elder K. Oxidative stress and alterations in DNA methylation: two sides of the same coin in reproduction. Reprod Biomed Online. 2016;33(6):668–83.

    Article  CAS  PubMed  Google Scholar 

  50. Zhu WJ. Transmission electron microscopy analysis of the origin and incidence of sperm intranuclear cytoplasmic retention in fertile and teratozoospermia men. Andrology. 2018;6(2):317–24.

    Article  CAS  PubMed  Google Scholar 

  51. Vertika S, Singh KK, Rajender S. Mitochondria, spermatogenesis, and male infertility - an update. Mitochondrion. 2020;54:26–40.

    Article  CAS  PubMed  Google Scholar 

  52. Ma P, Zhang Z, Zhou X, Luo J, Lu H, Wang Y. Characterizing semen abnormality male infertility using non-targeted blood plasma metabolomics. PLoS ONE. 2019;14(7):e0219179.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Hou Y, He Y, Wang J, Zhou Y, Ma B, Li X, Shi L, Zhang Q. Comprehensive metabolomics profiling of seminal plasma in asthenozoospermia caused by different etiologies. Clin Chim Acta. 2023;548:117530.

    Article  CAS  PubMed  Google Scholar 

  54. Chen L, Wen CW, Deng MJ, Ping L, Zhang ZD, Zhou ZH, Wang X. Metabolic and transcriptional changes in seminal plasma of asthenozoospermia patients. Biomedical Chromatography: BMC. 2020;34(3):e4769.

    Article  CAS  PubMed  Google Scholar 

  55. Boguenet M, Bocca C, Bouet PE, Serri O, Chupin S, Tessier L, Blanchet O, El Hachem H, de la Chao JM, Reynier P, et al. Metabolomic signature of the seminal plasma in men with severe oligoasthenospermia. Andrology. 2020;8(6):1859–66.

    Article  CAS  PubMed  Google Scholar 

  56. Najafi A, Asadi E, Moawad AR, Mikaeili S, Amidi F, Adutwum E, Safa M, Sobhani AG. Supplementation of freezing and thawing media with brain-derived neurotrophic factor protects human sperm from freeze-thaw-induced damage. Fertil Steril. 2016;106(7):1658–e654.

    Article  CAS  PubMed  Google Scholar 

  57. Najafi A, Amidi F, Sedighi Gilani MA, Moawad AR, Asadi E, Khanlarkhni N, Fallah P, Rezaiian Z, Sobhani A. Effect of brain-derived neurotrophic factor on sperm function, oxidative stress and membrane integrity in human. Andrologia. 2017;49(2).

  58. Xu J, Zhang W, Zhong S, Xie X, Che H, Si W, Tuo X, Xu D, Zhao S. Microcystin-leucine-arginine affects brain gene expression programs and behaviors of offspring through paternal epigenetic information. Sci Total Environ. 2023;857(Pt 1):159032.

    Article  CAS  PubMed  Google Scholar 

  59. Ito T, Osada A, Ohta M, Yokota K, Nishiyama A, Niikura Y, Tamura T, Sekita Y, Kimura T. SWI/SNF chromatin remodeling complex is required for initiation of sex-dependent differentiation in mouse germline. Sci Rep. 2021;11(1):24074.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Riera MF, Regueira M, Galardo MN, Pellizzari EH, Meroni SB, Cigorraga SB. Signal transduction pathways in FSH regulation of rat sertoli cell proliferation. Am J Physiol Endocrinol Metabolism. 2012;302(8):E914–23.

    Article  CAS  Google Scholar 

  61. Lv D, Ji Y, Zhang Q, Shi Z, Chen T, Zhang C, Wang X, Ren T, Gao Z, Zhong C. Mailuoshutong pill for varicocele-associated male infertility-phytochemical characterisation and multitarget mechanism. Front Pharmacol. 2022;13:961011.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Ehrmann I, Dalgliesh C, Tsaousi A, Paronetto MP, Heinrich B, Kist R, Cairns P, Li W, Mueller C, Jackson M, et al. Haploinsufficiency of the germ cell-specific nuclear RNA binding protein hnRNP G-T prevents functional spermatogenesis in the mouse. Hum Mol Genet. 2008;17(18):2803–18.

    Article  CAS  PubMed  Google Scholar 

  63. Ghieh F, Izard V, Poulain M, Fortemps J, Kazdar N, Mandon-Pepin B, Ferlicot S, Ayoubi JM, Vialard F. Cryptic splice site poisoning and meiotic arrest caused by a homozygous frameshift mutation in RBMXL2: a case report. Andrologia. 2022;54(11):e14595.

    Article  CAS  PubMed  Google Scholar 

  64. Nie DS, Liu Y, Juan H, Yang X. Overexpression of human SPATA17 protein induces germ cell apoptosis in transgenic male mice. Mol Biol Rep. 2013;40(2):1905–10.

    Article  CAS  PubMed  Google Scholar 

  65. Deng Y, Hu LS, Lu GX. Expression and identification of a novel apoptosis gene Spata17 (MSRG-11) in mouse spermatogenic cells. Acta Biochim Biophys Sin. 2006;38(1):37–45.

    Article  CAS  PubMed  Google Scholar 

  66. Nie D, Liu Y, Xiang Y. Overexpression a novel zebra fish spermatogenesis-associated gene 17 (SPATA17) induces apoptosis in GC-1 cells. Mol Biol Rep. 2011;38(6):3945–52.

    Article  CAS  PubMed  Google Scholar 

  67. Miyamoto T, Tsujimura A, Miyagawa Y, Koh E, Sakugawa N, Miyakawa H, Sato H, Namiki M, Okuyama A, Sengoku K. A single nucleotide polymorphism in SPATA17 may be a genetic risk factor for Japanese patients with meiotic arrest. Asian J Androl. 2009;11(5):623–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Wang J, Tang C, Wang Q, Su J, Ni T, Yang W, Wang Y, Chen W, Liu X, Wang S, et al. NRF1 coordinates with DNA methylation to regulate spermatogenesis. FASEB Journal: Official Publication Federation Am Soc Experimental Biology. 2017;31(11):4959–70.

    Article  CAS  Google Scholar 

  69. Ikeda S, Tanaka K, Ohtani R, Kanda A, Sotomaru Y, Kono T, Obata Y. Disruption of piRNA machinery by deletion of ASZ1/GASZ results in the expression of aberrant chimeric transcripts in gonocytes. J Reprod Dev. 2022;68(2):125–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Zhou Y, Connor EE, Bickhart DM, Li C, Baldwin RL, Schroeder SG, Rosen BD, Yang L, Van Tassell CP, Liu GE. Comparative whole genome DNA methylation profiling of cattle sperm and somatic tissues reveals striking hypomethylated patterns in sperm. GigaScience. 2018;7(5).

  71. Nguyen HT, Martin LJ. Classical cadherins in the testis: how are they regulated? Reproduction, fertility, and development. 2023;35(14):641–60.

  72. Johnson AR, Craciunescu CN, Guo Z, Teng YW, Thresher RJ, Blusztajn JK, Zeisel SH. Deletion of murine choline dehydrogenase results in diminished sperm motility. FASEB Journal: Official Publication Federation Am Soc Experimental Biology. 2010;24(8):2752–61.

    Article  CAS  Google Scholar 

  73. Johnson AR, Lao S, Wang T, Galanko JA, Zeisel SH. Choline dehydrogenase polymorphism rs12676 is a functional variation and is associated with changes in human sperm cell function. PLoS ONE. 2012;7(4):e36047.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Juanpanich T, Suttirojpattana T, Parnpai R, Vutyavanich T. The relationship between reactive oxygen species, DNA fragmentation, and sperm parameters in human sperm using simplified sucrose vitrification with or without triple antioxidant supplementation. Clin Experimental Reproductive Med. 2022;49(2):117–26.

    Article  Google Scholar 

  75. Khosravizadeh Z, Khodamoradi K, Rashidi Z, Jahromi M, Shiri E, Salehi E, Talebi A. Sperm cryopreservation and DNA methylation: possible implications for ART success and the health of offspring. J Assist Reprod Genet. 2022;39(8):1815–24.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank all the patients and healthy controls who participated in our study.

Funding

National Key Research and Development Program of China (2018YFE0207300), Beijing Natural Science Foundation (M23008), the National High Level Hospital Clinical Research Funding (2022-PUMCH-B-124), the National High Level Hospital Clinical Research Funding (2022-PUMCH-A-057).

Author information

Author notes

  1. Jingdi Zhang and Xiaogang Li contributed equally to this work.

Authors and Affiliations

  1. Department of Clinical Laboratory, State Key Laboratory of Complex Severe and Rare Diseases, Peking Union Medical College Hospital, Chinese Academy of Medical Science and Peking Union Medical College, No.1 Shuaifuyuan, Dongcheng District, Beijing, 100730, China

    Jingdi Zhang, Rongrong Wang, Xinxin Feng, Siyu Wang, Yongzhe Li & Ye Guo

  2. Department of Urology, State Key Laboratory of Complex Severe and Rare Diseases, Peking Union Medical College, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences, Beijing, China

    Hai Wang, Yutao Wang & Hongjun Li

  3. Medical Science Research Center, State Key Laboratory of Complex Severe and Rare Diseases, Peking Union Medical College Hospital, Chinese Academy of Medical Science and Peking Union Medical College, Beijing, China

    Xiaogang Li

Authors
  1. Jingdi Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Xiaogang Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Rongrong Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Xinxin Feng
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Siyu Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Hai Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Yutao Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Hongjun Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Yongzhe Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Ye Guo
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

J.D.Z: Data analysis and interpretation, Visualization, Writing - Original Draft. X.G.L: Conception and design, Writing - Review & Editing. R.R.W, X.X.F, S.Y.W, H.W, Y.T.W, H.J.L.: Collection and assembly of data. Y.Z.L: Funding acquisition, Revising and examining the manuscript. Y.G：Conception and design, Provision of study materials or patients, Collection and assembly of data.

Corresponding authors

Correspondence to Yongzhe Li or Ye Guo.

Ethics declarations

Ethics approval and consent to participate

The participates in our study written informed consent. This study was approved by the Peking Union Medical College Hospital ethics committee (No.JS-3112).

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Supplementary Material 3

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

Zhang, J., Li, X., Wang, R. et al. DNA methylation patterns in patients with asthenospermia and oligoasthenospermia. BMC Genomics 25, 602 (2024). https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-024-10491-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-024-10491-z

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us