Skip to main content

Effects of low lead exposure on sperm quality and sperm DNA methylation in adult men



Lead (Pb) exposure is a risk factor for male infertility, but the epigenetic changes in sperm DNAattributable to lead exposure is poorly defined.


In this study, we investigated whether low Pb exposure (<  10 µg/dL) affects the sperm quality. Blood, urine, and semen samples of 297 men of childbearing age were analyzed for all relevant parameters. Based on the blood Pb level (BLL), participants were allocated to RL (0–2.5 µg/dL), RM (2.5–5 µg/dL), and RH (5–10 µg/dL) groups. The 5-methylcytosine and 5-hydroxymethylcytosine patterns in the sperm DNA were identified using methylated DNA immunoprecipitation and hydroxymethylated DNA immunoprecipitation sequencing.


The non-progressive motility (NP) was significantly increased and associated with global hypomethylation of sperm DNA in the RH group compared with the RL group, indicating that aberrant sperm methylation due to low Pb exposure is possibly associated with reduced sperm motility. The hypomethylated promoter regions were primarily enriched in the calcium (Ca) homeostasis pathway. Further, the interaction between Ca and Pb was associated with sperm rapid progressive motility and asthenospermia risk, although no significant methylation abnormality was observed in those with BLL  <  5 µg/dL. When BLL was  >  5 µg/dL or when predicting NP, no significant Pb–Ca interaction was observed.


Overall, our results indicate that aberrant DNA methylation of the Ca homeostasis pathway, induced by low Pb exposure, is the potential cause for reduced sperm velocity.


Research suggests that the effect of Pb exposure on the sperm quality, an important indicator of male fertility, is dose-dependent [1]. According to the Centers for Disease Control (CDC), a blood Pb level (BLL) of  >  10 µg/dL is a public health concern [2]. However, studies have shown that the Pb exposure in men with an average BLL of  <  10 µg/dL causes reduced sperm quality. For an example, a study by Jurasovic et al. [3] in 123 Croatian volunteers with a median BLL of 5.7 µg/dL (range: 2.5–14.9 µg/dL) reported that the BLL is negatively correlated with sperm motility and morphology. Similar results were observed in another study involving 240 volunteers with a median BLL of 4.92 µg/dL (range: 1.13–14.91 µg/dL) [4]. However, it is well known that a wide concentration range may result in heteroscedasticity of data, where the variance increases with increasing concentrations [5]. The analytical results of investigations examining the biological effects of low Pb exposure may be inaccurate, especially at a lower concentration range, if individuals exposed with a high BLL are included. Notably, previous studies have failed to exclude individuals with a high BLL [3, 4]. In this study, we explored the potential effects of low Pb exposure on the sperm quality. Only men with a BLL  <  10 µg/dL were included.

Large-scale DNA methylation reprogramming occurs during gamete formation and early embryonic period in humans, and the process is susceptible to environmental factors [6, 7]. However, previous relevant research has mainly focused on determining the effect of environmental Pb exposure on DNA methylation in peripheral blood [8]. Christine et al. [9] demonstrated that maternal Pb exposure (mean BLL: 5 µg/dL) can cause changes in peripheral blood DNA methylation of important genes involved in fetal development, in line with the results of a study by Shaowei et al. [10]. In addition, DNA methylation is one of the key mechanisms to maintain Ca homeostasis. A study had shown that the Ca homeostasis pathway regulated by DNA methylation is related to the risk of coronary heart disease [11]. The maintenance of Ca homeostasis is the key factor to ensure sperm motility [12, 13]. However, only a few studies have investigated the possible effects of low-level parental Pb exposure on human sperm DNA methylation and Ca homeostasis pathway regulated by DNA methylation.

In this population-based study, we determined whether the level of Pb exposure represented by a BLL of  <  10 µg/dL is an independent risk factor for poor sperm quality. In addition, we used sperm DNA methylation data, quantified through methylated DNA immunoprecipitation (MeDIP) and hydroxymethylated DNA immunoprecipitation (hMeDIP) sequencing, to explore the mechanism associating low Pb exposure with sperm quality.

Materials and methods

Study participants

This cross-sectional study involving men of childbearing age living in the Sandu County, Guizhou Province, China, was conducted from September 2011 to April 2012 by using random cluster sampling.

After excluding those with left testis and right testis volume of  <  12 mL, sperm density  <  1  ×  106/mL abnormal semen agglutination, and blood Pb concentration  >  10 µg/dL, a total of 297 men were included in this study, with an average age of 32.71 years. The severity of smoking habit was recorded as the number of cigarettes smoked per day, and 271 participants reported of daily alcohol intake. The study was approved by the Ethics Committee and Institutional Review Board of Shanghai Institute of Planned Parenthood Research. All the study participants provided written informed consent.

Semen analysis

All participants had to abstain from sex for 2–7 days before providing the semen samples. The samples were collected by masturbation into 25-mL sterile polystyrene jars and analyzed within 1 h of ejaculation. Computer-assisted sperm analysis (CASA) (WLJY-9000, Beijing, China) was performed to obtain sperm concentration, viability, and motility after the liquefaction of sample at the clinics, in accordance with the World Health Organization (WHO) guidelines (fifth edition) [14].

An aliquot of semen was centrifuged at 3000 rpm for 5 min. The supernatant was frozen without preservatives and temporarily stored at − 20 °C. All the samples were shipped to the laboratory at the Shanghai Institute of Planned Parenthood Research (Shanghai, China) on dry ice and stored at − 80 °C for the Ca, Mg, and Zn concentration assays.

Total Zn, Ca, and Mg concentrations in semen were measured through flame atomic absorption spectrometry by using a BH5500S device (Bohui Innovation Technology Co., Ltd., Beijing, China). The calibration solution in the determinations was provided by the Bohui Company, and accuracy of the calibration solution was set according to the national standard material for Zn (GBW08620), Ca [GBW(E)080118], and Mg [GBW(E)080126].

Blood analysis

The nursing staff obtained 5 mL of venous blood into anticoagulant tubes from participants’ upper limbs. Similar to Zn, Ca, and Mg assays, the Pb and cadmium (Cd) levels in the blood were determined using a graphite furnace atomic absorption spectrophotometer (same as the detection scheme of trace elements in seminal plasma). The labeling curve was drawn by using the standard solution of Pb and Cd provided by China National Research Center for Labeling Substances. Serum FSH, LH, and total testosterone were determined using an electrochemiluminescence immunoassay with kits provided by Roche Diagnostics GmbH (Mannheim, Germany).

Urine analysis

A single-spot urine sample was collected from all the participants during their visit to the clinic. The urine samples were analyzed at the collaborative laboratory of East China University of Science and Technology (Shanghai, China). Urinary concentrations of total BPA, DEHP, MEHP and DBP were measured using modified high-performance liquid chromatography. The same method was used to detect Zn, Ca, and Mg in both urine and semen samples.

MeDIP and hMeDIP sequencing

Sperm DNA of nine men (three men randomly selected from each of the RL, RM, and RH groups) aged 20–40 years were submitted for MeDIP and hMeDIP sequencing. These men did not have a history of smoking or alcohol intake, their blood Cd concentration was  <  1 µg/dL, and BPA, DEHP, MEHP, and DBP concentrations were  <  1 µg/µmolCr. The sperm DNA was extracted from samples using the DNeasy Blood and Tissue kit (Qiagen) according to the manufacturer’s instructions. Agarose gel electrophoresis was used to detect the integrity of genome, Nano drop (Thermo Scientific) was used to detect DNA purity (OD 260/280 ratio), and the DNA concentration was measured using Qubit 2.0 fluorometer (Life Technologies). The extracted DNA was fragmented using the Covaris sonication system, and sequencing libraries were prepared from the 5-μg fragments of genomic DNA. End repair, A base addition, and adaptor ligation steps were performed using the Illumina’s Single-End DNA sample preparation kit. Adaptor-ligated DNA was immunoprecipitated by anti-5mC and anti-5hmC by using a commercial antibody. The immunoprecipitated DNA was purified and then applied for 50-bp single-end sequencing on the Illumina Hiseq2500 platform.

Bioinformatic analysis

All sequencing data passed initial quality checks for base composition through FASTQC v0.10.0 (15). For each individual,  ~  20 million reads were generated and mapped onto hg19 by using BWA [16]. After removing duplicates, we used macs2 [17] to call peaks with the following parameters: -t $path-f BAM-g hs-B-n $path_suffix-broad-broad-cutoff 0.1-fix-bimodal-extsize 300. Then, R package ChIP-seq was used to evaluate the quality of samples, and the diffbind package was used to evaluate DMRs between the samples [18, 19]. The screening threshold of DMRs was the log2(fold change)  >  1, FDR  <  0.05, and the differential binding region was annotated by CHIPseeker [20]. The R package clusterprofiler [20] was used to annotate DMRs in the promoter region of the gene and to perform Gene Ontology (GO) [21], Kyoto Encyclopedia of Genes and Genomes (KEGG) [22], and Ingenuity Pathway Analysis (IPA) analysis. The gene where the differential promoter region was located was used to create a heatmap by using the pheatmap Package.

Statistical analysis

Demographic characteristics and blood, semen, and urine parameters were statistically analyzed. The sample frequency was counted by classified variables, and the mean value, standard deviation were counted by numerical variables. Logistic regression analysis was used to classify variables and linear regression analysis was used to classify numerical variables, to explore factors affecting male fertility. By including the interaction term Pb and Ca in sperm in the full model, we investigate whether relationships of Pb and response changed based on the value of Ca in sperm. P values  <  0.05 were considered statistically significant. Apart from considering significant interaction effect, we compare two nested models that whether included interaction term or not by ANOVAs test to select the best fit. The R package stats was used to implement the regression analysis and ggpubr was used to display the results. All statistical analysis were performed using the R 3.6.3 package.


Description of study participants

Blood samples were obtained from all the participants for assessing hormone, Pb, and Cd levels. The mean hormone levels were as follows: follicle stimulating hormone (FSH): 5.14 IU/L; luteinizing hormone (LH): 5.26 U/L; and testosterone: 5.55 ng/mL. The mean BLL was 3.15 µg/dL, and the mean blood Cd level was 5.87 µg/dL. Urine analysis revealed the following plastifier levels (mean values): bisphenol A (BPA): 1.71 µg/µmolCr; dioctyl phthalate (DEHP): 6.09 µg/µmolCr; dibutyl phthalate (DBP): 2.39 µg/µmolCr; and mono-2-ethylhexyl phthalate (MEHP): 3.22 µg/µmolCr. The mean urine zinc (Zn), Ca, and magnesium (Mg) levels were 10.34, 3.06, and 1.98 mmol/L, respectively. Semen samples were evaluated for rapid progressive motility (mean, 27.35%), slow progressive motility (14.01%), non-progressive motility (13.38%), sperm concentration (54.54  ×  106/mL), and sperm viability (64.26%). Seminal plasma samples were analyzed for Zn (mean, 2.42 mmol/L), Ca (9.52 mmol/L), and Mg (4.48 mmol/L). Testicular scale was used to measure the testicular volume; the mean left testicular volume was 22.79 mL and right testicular volume was 22.6 mL. The participants were allocated to three groups based on the BLL: RH group (BLL: 5–10 µg/dL; n  =  68), RM group (BLL: 2.5–5 µg/dL; n  =  91), and RL group (BLL: 0–2.5 µg/dL; n =  138; Table 1).

Table 1 Characteristics of study participants

Risk factor assessment for sperm quality

A preliminary multiple regression model was used to analyze the association between BLL and semen parameters. The non-progressive motility was significantly increased in the RH group compared with the RL group [coefficient: 2.72, 95% confidence interval (CI) 0.075–5.37; p  =  0.044). However, no significant association was observed between Pb exposure and the proportion of rapid progressive motility, slow progressive motility, sperm viability, sperm density, and risk for asthenospermia after adjusting for age, essential metal elements, plastifiers, and smoking or alcohol status (Fig. 1). Non-progressive motility is the inability of the sperm to swim efficiently toward the egg and is considered a risk factor of asthenospermia. Our findings suggest that BLL in the range of 5–10 µg/dL is also an independent risk factor for poor sperm quality.

Fig. 1
figure 1

Forest plot showing adjusted coefficient or odds ratio (OR) and 95% confidence interval (CI) for semen parameters in the low Pb exposure group. (RL: 2.5 µL/dL  ≥  Blood lead level  ≥  0 µg/dL; RL: 5 µL/dL  ≥  Blood lead level  >  2.5 µg/dL; RL: 10 µL/dL  ≥  Blood lead level  >  5 µg/dL)

MeDIP and hMeDIP sequencing of sperm DNA exposed to low Pb level

To determine whether the poor sperm quality is due to low-level Pb exposure-induced impaired DNA methylation, two important methylation forms of sperm DNA, namely 5-methylcytosine (mC) and 5-hydroxymethylcytosine (hmC), were assessed using the MeDIP and hMeDIP sequencing. Approximately 20 million reads were generated for each sample, and the uniquely aligned ratio was about 70% (Additional file 1: Table S1). Hierarchical clustering analysis indicated extremely different sperm DNA methylation patterns in the RH and RL groups. The RM group did not display a considerable difference in the DNA methylation pattern compared with the RL group (Fig. 2A, B). According to the principal component analysis (PCA), PC1 could explain  >  80% methylation variation. Difference between the RH and RL groups in the PC1 direction was observed; however, no obvious difference in the PC1 was observed between RM and RL groups (Fig. 2C, D). We also evaluated the average signal profile across peaks, and the results revealed that the 5mC signal in the RH group decreases globally compared with the RL group and that the sperm methylation profile in the RM group is not obviously different from that in the RL group (Fig. 2E, F). These findings suggest that global sperm DNA hypomethylation induced by Pb exposure may be associated with an increased proportion of sperms with non-progressive motility. We also analyzed the distribution of different hydroxymethylation regions between the RH and RL groups. Additional file 7: Figure S1A, B indicate no obvious difference in the hydroxymethylation signal and profile between the RH and RL groups.

Fig. 2
figure 2

5mC methylation partners in different environments Pb exposure groups. A Correlation heatmap showing 5 mC methylation correlation between the RH group (blood Pb, 5–10 µg/dL) and the RL group (blood Pb, 0–2.5 µg/dL). B Correlation heatmap showing 5 mC methylation correlation between the RM group (blood Pb, 2.5–5 µg/dL) and RL group (blood Pb, 0–2.5 µg/dL). C Principal component analysis (PCA) of methylation measurements in the RH and RL groups. D PCA of methylation measurements in the RM and RL groups. EComparison of the average coverage profile at  ±  200 bp near the 5 mC methylation peaks summit between the RH and RL groups. F Comparison of the average coverage profile at  ±  200 bp near the 5 mC methylation peaks summit between the RM and RL groups

Distribution of differentially methylated regions and hydroxymethylated regions

MACS2 was used to identify the distribution of different methylation regions between the RH and RL groups. A total of 18,259 significantly different peaks were identified in the whole genome, of which 17,558 peaks were downregulated and 701 were upregulated (Fig. 3A; Additional file 2: Table S2). These differential peaks had a high reads alignment (Fig. 3B). The 5 mC methylation of the promoter region can cause a low expression of related genes. We screened 1401 genes (1373 hypomethylation genes and 28 hypermethylation genes; Fig. 3C, Additional file 3: Table S3) related to environmental Pb exposure.The promoter regions of these genes contain at least one differentially methylated peaks. Because the methylation variation between the RM and RL groups was small, we found only one significant differential peak (data not shown). Finally, we obtained 12 differentiated peaks, and interestingly, most of them were found to be located on the Y chromosome (Additional file 8: Figure S2).

Fig. 3
figure 3

A Volcano plot showing differentially methylated regions between the RH and RL groups. B MA plots of the log average of the methylation levels on the X-axis and log ratio of the methylation levels on the Y-axis, RH and RL groups. C Pie charts depicting genes with differential methylation in the promoter regions

Enrichment analysis of DMR-related genes

We performed the enrichment analysis of differential genes. The enriched GO terms were related to biological processes, as shown in Fig. 4A; Additional file 4: Table S4. Exposure to Pb seemed to alter several pathways in the KEGG enrichment analysis, as is shown in Fig. 4B; Additional file 5: Table S5. Genes with methylation differences in the promoter regions were submitted to the IPA analysis. The analysis showed that differentially methylated genes are mainly related to the neurological system, and cytoskeleton and Ca pathways (Fig. 4D; Additional file 6: Table S6), in agreement with the enriched terms identified using the GO and KEGG. Moreover, the promoter regions of the Ca pathway genes enriched in GO, KEGG, and IPA were all hypomethylated. We further tried to understand the biological relevance of these genes associated with differentially methylated regions (DMRs) in the Ca pathway homeostasis regulation.

Fig. 4
figure 4

A Gene Ontology enrichment analysis for biological processes. B KEGG signaling pathways analysis. C Visualization of the previously determined calcium signaling pathway using the Pathview package based on KEGG pathways. D Ingenuity pathway analysis (IPA) of differentially regulated genes. E Visualization of genes enriched in the calcium signaling pathway in IPA analysis. F Heatmap summarizing significantly differentially methylated genes enriched in the calcium homeostasis pathway in all three databases (GO, KEGG, and IPA)

Pathway analysis revealed 17 genes with promoter hypomethylation in the KEGG Ca signaling pathway and 15 genes with promoter hypomethylation in the IPA Ca signaling pathway (Fig. 4C, E). Eleven differentially methylated genes simultaneously enriched in the GO (Ca ion transport), KEGG (Ca signaling pathway), and IPA (Ca signaling pathway) were highlighted in the heatmap (Fig. 4F). These findings indicate that the mechanism of Pb toxicity is not only limited to essential metals and metal-dependent biological enzymes but also affects the methylation level of promoters related to the Ca homeostasis pathway genes.

Effect modification of Ca and Pb on sperm motility

Although no significant association of low Pb exposure or Ca concentration with sperm progressive motility and density was observed in the preliminary model, DNA methylation results showed a close association between Pb and Ca homeostasis pathways, even in the population with low Pb exposure. Therefore, we further evaluated the presence of interactions between Pb and Ca in these individuals. A test for interaction revealed a two-way effect modification between Pb exposure and seminal Ca concentration for predicting rapid progressive motility. However, this effect modification was found only in the RM group (coefficient: 1.2, 95% CI 0.12–2.2; p  =  0.029) compared with the RL group. Results of the nested model test indicated that the full model that involves the interaction effect is superior to the preliminary model in predicting rapid progressive motility of sperms (Table 2; p  =  0.057). A relationship between low Pb exposure and sperm motility was not observed in the preliminary model possibly because the interaction between Ca and Pb masked the influence of Pb exposure on sperm motility. We also did not observe a significant interaction between urinary Ca and low Pb exposure, and therefore, we excluded it from the full model (data not shown). In the full model of other semen parameters, we did not observe a significant interaction between Pb and Cd.

Table 2 Effect modification of seminal calcium and Pb of sperm quality and nested model comparisons

Furthermore, we analyzed whether low Pb exposure is related to the risk of asthenospermia, as stated in the fifth edition of the WHO guidelines. The results demonstrated that there is a two-way effect modification between Pb exposure and seminal Ca concentration for predicting asthenospermia, and the effect was found to exist only in the RM group when compared with the RL group (OR: 0.83, 95% CI 0.68–0.99; p  =  0.05). Results of the nested model test supported the superiority of the full model that involves the interaction effect to the preliminary model (p  =  0.030; Table 2).


Our findings showed that a BLL of 5–10 µg/dL can still reduce the sperm motility in men of childbearing age, which is inconsistent with the results of a previous Mexican cohort study [23]. The primary reason for this difference may be the lack in previous studies of assessment of other risk factors, such as Cd, trace elements, or plastifiers, that are associated with male reproductive toxicity. Unlike a recent study in Southwest China, which reported that low Pb exposure adversely effects the sperm movement and density, our research shows that low Pb exposure is related only to the sperm motility. The difference may be that the individuals with BLL  >  10 µg/dL were not excluded in the study in Southwest China [24].

Majority of studies have focused on determining the effects of environmental Pb exposure on DNA methylation in human peripheral blood cells [9, 25, 26]. Moreover, analysis in some studies that focused on determining the epigenetic effects of Pb exposure on the male reproductive system had been limited to sperm DNA damage and changes in the sperm chromatin [23, 27]. Our research confirms that the effect of BLL as low as 5 µg/dL on decreased sperm motility can be similar to that of the occupational Pb exposure. We found that low Pb exposure (BLL of 5–10 vs 0–2.5 µg/dL) changes the DNA methylation pattern of human sperm and causes global hypomethylation, which established the mechanism linking the decreased sperm motility with low Pb exposure and DNA methylation. However, we did not observe a global change in sperm DNA methylation in the 2.5–5 µg/dL group. The finding indicates that the effect of Pb exposure on sperm DNA methylation is dose-dependent and that a BLL of 5 µg/dL could be used as a threshold value to indicate the safety of DNA methylation in male gametes. In addition, we detected the DNA hydroxymethylation level in the RH and RL groups and observed no global difference in the sperm DNA hydroxymethylation, which imply that hydroxymethylation may be insensitive to low Pb exposure. Compared with other studies, this study expands the perception of the effect of environmental Pb exposure on DNA methylation of male gametes [28, 29].

We also identified 1401 differentially methylated genes (DMGs) in promoter regions. Among them, the promoter regions of 1373 genes were hypomethylated, suggesting that environmental Pb exposure may be associated with abnormally high expression of genes in sperm cells. This finding is consistent with a study that reported a correlation between Pb exposure and decreased DNA methylation [30]. We analyzed the enrichment of these DMGs through IPA, GO, and KEGG analysis. Results were consistent and indicated that the functions of DMGs are mainly enriched in the cytoskeleton and neurological system pathways. Abnormal cytoskeleton has been reported to be an important factor for sperm motility [31]. Our result confirms that the reduced sperm motility reduction is associated with the DNA methylation pattern caused by low Pb exposure. Moreover, significant enrichment of DMGs in the neurological system necessitates further study to understand the association between paternal Pb exposure and offspring nervous system disorders [32].

Our findings revealed that low Pb exposure causes hypomethylation and hence impairs the Ca homeostasis pathway genes in sperm cells. However, in the preliminary model, we did not observe a significant correlation between seminal Ca or urinary Ca and sperm progressive motility. Previous studies have reported an antagonistic interaction between Pb and Ca; therefore, Ca supplementation is considered an effective antidote for chronic Pb poisoning [33,34,35]. We assume that Ca influx may antagonize the effect of low Pb exposure on sperm motility. In the full model with the interaction factors of blood Pb and seminal Ca, we observed statistically significant modification effects of low Pb exposure and seminal Ca on the rapid progressive motility of sperms and risk of asthenospermia as well.

Genomic regulation through epigenetic modifications in Ca homeostasis-related pathways is an important mechanism to maintain human health [11]. Our research showed that Ca may be a protective factor for men of childbearing age with a BLL of  <  5 µg/dL and may help to avoid the deterioration of sperm motility induced by Pb exposure. However, in those with a BLL of  >  5 µg/dL, we observed no significant interaction between Pb and Ca accompanied by abnormal methylation of Ca homeostasis regulatory genes in rapid progressive motility sperm cells. This finding suggests that the interaction between Ca and Pb is epigenetically controlled under the condition of low Pb exposure. Moreover, Ca does not protect against abnormal DNA methylation caused by high Pb exposure beyond 5 µg/dL. The potential reason may be Pb exposure beyond 5 µg/dL interfered methylation-dependent Ca homeostasis regulation. Similarly, we did not observe the interaction between Pb and Ca in non-progressive motility sperm cells which have been reported carrying inherent methylation abnormalities [36]. This potential mechanism may account for the ineffectiveness of Ca supplementation in alleviating Pb poisoning.

Our study has certain limitations. The results must be interpreted in the context of the Chinese population and limitations inherent to human data. Moreover, although data were adjusted for multiple known confounding factors such as age, smoking and alcohol status, endocrine disruptors, and Cd, individuals are typically exposed to factors, such as chemical stressors (exposure to nearly 200 chemicals) and nonchemical stressors (anxiety and stress), that can affect the sperm quality. We, therefore, cannot exclude the possibility that the reported results may be related to other exposure factors. These limitations notwithstanding our data suggest that Pb exposure is associated with epigenetic perturbations, including Ca pathway gene hypomethylation, which has been linked to various diseases. Currently, another cohort study in the eastern China is investigating the impact of environmental Pb exposure on the Ca pathway.

Availability of data and materials

All underlying data are provided in the manuscript or available upon request.



Blood lead level


Non-progressive motility


Methylated DNA immunoprecipitation


Hydroxymethylated DNA immunoprecipitation


Follicle stimulating hormone


Luteinizing hormone


Dioctyl phthalate


Bisphenol A


Dibutyl phthalate


Mono-2-ethylhexyl phthalate


Differentially methylated genes


  1. Onul NM, Biletska EM, Stus VP, Polion MY. The role of lead in the etiopathogenesis of male fertility reduction. Wiad Lek. 2018;71:1155–60.

    PubMed  Google Scholar 

  2. Frumkin H, Buchanan S, Brown MJ. Interpreting and managing blood lead levels. MMWR Recomm Rep. 2007.

  3. Jurasovic J, Cvitkovic P, Pizent A, Colak B, Telisman S. Semen quality and reproductive endocrine function with regard to blood cadmium in Croatian male subjects. Biometals. 2004;17:735–43.

    Article  CAS  Google Scholar 

  4. Telisman S, Colak B, Pizent A, Jurasovic J, Cvitkovic P. Reproductive toxicity of low-level lead exposure in men. Environ Res. 2007;105:256–66.

    Article  CAS  Google Scholar 

  5. Sandeep S, Sonawane SSC, Sana S, Sanjay JK. An approach to select linear regression model in bioanalytical method validation. J Anal Sci Technol. 2019.

    Article  Google Scholar 

  6. Messerschmidt DM, Knowles BB, Solter D. DNA methylation dynamics during epigenetic reprogramming in the germline and preimplantation embryos. Genes Dev. 2014;28:812–28.

    Article  CAS  Google Scholar 

  7. Seisenberger S, Peat JR, Hore TA, Santos F, Dean W, Reik W. Reprogramming DNA methylation in the mammalian life cycle: building and breaking epigenetic barriers. Philos Trans R Soc Lond B Biol Sci. 2013;368:20110330.

    Article  Google Scholar 

  8. Hanna CW, Bloom MS, Robinson WP, Kim D, Parsons PJ, vom Saal FS, Taylor JA, Steuerwald AJ, Fujimoto VY. DNA methylation changes in whole blood is associated with exposure to the environmental contaminants, mercury, lead, cadmium and bisphenol A, in women undergoing ovarian stimulation for IVF. Hum Reprod. 2012;27:1401–10.

    Article  CAS  Google Scholar 

  9. Rygiel CA, Dolinoy DC, Perng W, Jones TR, Solano M, Hu H, Tellez-Rojo MM, Peterson KE, Goodrich JM. Trimester-specific associations of prenatal lead exposure with infant cord blood DNA methylation at birth. Epigenet Insights. 2020;13:2516865720938669.

    Article  Google Scholar 

  10. Wu S, Hivert MF, Cardenas A, Zhong J, Rifas-Shiman SL, Agha G, Colicino E, Just AC, Amarasiriwardena C, Lin X, Litonjua AA, DeMeo DL, Gillman MW, Wright RO, Oken E, Baccarelli AA. Exposure to low levels of lead in utero and umbilical cord blood DNA methylation in project viva: an epigenome-wide association study. Environ Health Perspect. 2017;125:087019.

    Article  Google Scholar 

  11. Agha G, Mendelson MM, Ward-Caviness CK, Joehanes R, Huan T, Gondalia R, Salfati E, Brody JA, Fiorito G, Bressler J, Chen BH, Ligthart S, Guarrera S, Colicino E, Just AC, Wahl S, Gieger C, Vandiver AR, Tanaka T, Hernandez DG, Pilling LC, Singleton AB, Sacerdote C, Krogh V, Panico S, Tumino R, Li Y, Zhang G, Stewart JD, Floyd JS, Wiggins KL, Rotter JI, Multhaup M, Bakulski K, Horvath S, Tsao PS, Absher DM, Vokonas P, Hirschhorn J, Fallin MD, Liu C, Bandinelli S, Boerwinkle E, Dehghan A, Schwartz JD, Psaty BM, Feinberg AP, Hou L, Ferrucci L, Sotoodehnia N, Matullo G, Peters A, Fornage M, Assimes TL, Whitsel EA, Levy D, Baccarelli AA. Blood leukocyte DNA methylation predicts risk of future myocardial infarction and coronary heart disease. Circulation. 2019;140:645–57.

    Article  CAS  Google Scholar 

  12. Darszon A, Nishigaki T, Wood C, Trevino CL, Felix R, Beltran C. Calcium channels and Ca2+ fluctuations in sperm physiology. Int Rev Cytol. 2005;243:79–172.

    Article  CAS  Google Scholar 

  13. Xu Y, Fan Y, Fan W, Jing J, Xue K, Zhang X, Ye B, Ji Y, Liu Y, Ding Z. RNASET2 impairs the sperm motility via PKA/PI3K/calcium signal pathways. Reproduction. 2018;155:383–92.

    Article  CAS  Google Scholar 

  14. Cooper TG. WHO laboratory manual for the examination and processing of human semen, 5th ed. Research, WHO D. o. R. H. a; 2010.

  15. Andrews S. A Quality Control Tool for High Throughput Sequence Data. FastQC: [Online]. 2010.

  16. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.

    Article  CAS  Google Scholar 

  17. Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9:R137.

    Article  Google Scholar 

  18. Sarkar D, Gentlemen R, Lawrence M, Yao Z. chipseq: chipseq: a package for analyzing chipseq data. R package version 1.42.0. Berlin: Springer; 2021.

    Google Scholar 

  19. Ross-Innes CS, Stark R, Teschendorff AE, Holmes KA, Ali HR, Dunning MJ, Brown GD, Gojis O, Ellis IO, Green AR, Ali S, Chin SF, Palmieri C, Caldas C, Carroll JS. Differential oestrogen receptor binding is associated with clinical outcome in breast cancer. Nature. 2012;481:389–93.

    Article  CAS  Google Scholar 

  20. Yu G, Wang LG, He QY. ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization. Bioinformatics. 2015;31:2382–3.

    Article  CAS  Google Scholar 

  21. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G. Gene ontology: tool for the unification of biology. The gene ontology. Consortium Nat Genet. 2000;25:25–9.

    Article  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  23. Hernandez-Ochoa I, Garcia-Vargas G, Lopez-Carrillo L, Rubio-Andrade M, Moran-Martinez J, Cebrian ME, Quintanilla-Vega B. Low lead environmental exposure alters semen quality and sperm chromatin condensation in northern Mexico. Reprod Toxicol. 2005;20:221–8.

    Article  CAS  Google Scholar 

  24. Ren J, Cui J, Chen Q, Zhou N, Zhou Z, Zhang GH, Wu W, Yang H, Cao J. Low-level lead exposure is associated with aberrant sperm quality and reproductive hormone levels in Chinese male individuals: Results from the MARHCS study low-level lead exposure is associated with aberrant sperm quality. Chemosphere. 2020;244:125402.

    Article  CAS  Google Scholar 

  25. Wright RO, Schwartz J, Wright RJ, Bollati V, Tarantini L, Park SK, Hu H, Sparrow D, Vokonas P, Baccarelli A. Biomarkers of lead exposure and DNA methylation within retrotransposons. Environ Health Perspect. 2010;118:790–5.

    Article  CAS  Google Scholar 

  26. Kovatsi L, Georgiou E, Ioannou A, Haitoglou C, Tzimagiorgis G, Tsoukali H, Kouidou S. p16 promoter methylation in Pb2+ -exposed individuals. Clin Toxicol (Phila). 2010;48:124–8.

    Article  Google Scholar 

  27. Li C, Zhao K, Zhang H, Liu L, Xiong F, Wang K, Chen B. Lead exposure reduces sperm quality and DNA integrity in mice. Environ Toxicol. 2018;33:594–602.

    Article  CAS  Google Scholar 

  28. Pacchierotti F, Spano M. Environmental impact on DNA methylation in the germline: state of the art and gaps of knowledge. Biomed Res Int. 2015;2015:123484.

    Article  Google Scholar 

  29. Khalid M, Abdollahi M. Epigenetic modifications associated with pathophysiological effects of lead exposure. J Environ Sci Health C Environ Carcinog Ecotoxicol Rev. 2019;37:235–87.

    Article  CAS  Google Scholar 

  30. Ruiz-Hernandez A, Kuo CC, Rentero-Garrido P, Tang WY, Redon J, Ordovas JM, Navas-Acien A, Tellez-Plaza M. Environmental chemicals and DNA methylation in adults: a systematic review of the epidemiologic evidence. Clin Epigenetics. 2015;7:55.

    Article  Google Scholar 

  31. Godinez-Solis Y, Solis-Heredia MJ, Roa-Espitia A, Parra-Forero LY, Hernandez-Gonzalez EO, Hernandez-Ochoa I, Quintanilla-Vega B. Low concentrations of lead decrease the sperm fertilization ability by altering the acrosome reaction in mice. Toxicol Appl Pharmacol. 2019;380:114694.

    Article  CAS  Google Scholar 

  32. Brady K, Herrera Y, Zenick H. Influence of parental lead exposure on subsequent learning ability of offspring. Pharmacol Biochem Behav. 1975;3:561–5.

    Article  CAS  Google Scholar 

  33. Vigeh M, Smith DR, Hsu PC. How does lead induce male infertility? Iran J Reprod Med. 2011;9:1–8.

    PubMed  PubMed Central  Google Scholar 

  34. Ettinger AS, Lamadrid-Figueroa H, Tellez-Rojo MM, Mercado-Garcia A, Peterson KE, Schwartz J, Hu H, Hernandez-Avila M. Effect of calcium supplementation on blood lead levels in pregnancy: a randomized placebo-controlled trial. Environ Health Perspect. 2009;117:26–31.

    Article  CAS  Google Scholar 

  35. He Y, Zou Q, Chen H, Weng S, Luo T, Zeng X. Lead inhibits human sperm functions by reducing the levels of intracellular calcium, cAMP, and tyrosine phosphorylation. Tohoku J Exp Med. 2016;238:295–303.

    Article  CAS  Google Scholar 

  36. Capra E, Lazzari B, Turri F, Cremonesi P, Portela AMR, Ajmone-Marsan P, Stella A, Pizzi F. Epigenetic analysis of high and low motile sperm populations reveals methylation variation in satellite regions within the pericentromeric position and in genes functionally related to sperm DNA organization and maintenance in Bos taurus. BMC Genomics. 2019;20:940.

    Article  CAS  Google Scholar 

Download references


We would like to thank TopEdit ( for English language editing of this manuscript.


This work was supported by National Key R&D Program of China (2018YFC1005001), The Foundation of Science and Technology Commission of Shanghai Municipality (17JC1420103),Shanghai Municipal Science and Technology Major Project (2017SHZDZX01),Youth Program of National Natural Science Foundation of China (31801252),Youth Program of Shanghai Municipal Health Commission (20204Y0276),Shanghai Municipal Health Commission Scientific Research Project (201640369).

Author information

Authors and Affiliations



JQ, MM and YG take primary responsibility for the paper. TZ and YR collected the data. TZ, BW and HD did the statistical analysis. LC and JZ coordinated the research. JL, XW, ZW, XW, XS and JW participated in the data analysis. HS drafted the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Jun Qian, Maohua Miao or Yihua Gu.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethics Committee and Institutional Review Board of Shanghai Institute of Planned Parenthood Research. All the study participants provided written informed consent.

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: Table S1.

Reads alignment ratio.

Additional file 2: Table S2.

Significant differential 5 mC peaks between the RH and RL groups.

Additional file 3: Table S3.

Genes differentially 5 mC methylated in the promoter regions between the RH and RL groups.

Additional file 4: Table S4.

Gene Ontology analysis result.

Additional file 5: Table S5.

Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis results.

Additional file 6: Table S6.

Ingenuity pathway analysis results.

Additional file 7: Figure S1.

(A) Correlation heatmap showing the 5hmC correlation between the RH and RL groups. (B) Principal component analysis (PCA) of 5hmC measurements in the RH and RL groups. (C)Comparison of the average coverage profile at ±200 bp near 5hmC peaks between the RH and RL groups.

Additional file 8: Figure S2.

Heatmap summarizing the significant differential hydroxymethylcytosine regions between the RH and RL groups.

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 The Creative Commons Public Domain Dedication waiver ( 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

Verify currency and authenticity via CrossMark

Cite this article

Zhang, T., Ru, Y.F., Wu, B. et al. Effects of low lead exposure on sperm quality and sperm DNA methylation in adult men. Cell Biosci 11, 150 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Low lead
  • Sperm motility
  • DNA methylation
  • Calcium and lead interaction
  • Calcium homeostasis regulation pathway