Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genome-wide alteration in DNA hydroxymethylation in the sperm from bisphenol A-exposed men

  • Huajun Zheng ,

    Contributed equally to this work with: Huajun Zheng, Xiaoyu Zhou

    Affiliations Key Laboratory of Reproduction Regulation of NPFPC, SIPPR, IRD, Fudan University, Shanghai, China, WHO Collaborating Center for Research in Human Reproduction, Shanghai, China

  • Xiaoyu Zhou ,

    Contributed equally to this work with: Huajun Zheng, Xiaoyu Zhou

    Affiliations Key Laboratory of Reproduction Regulation of NPFPC, SIPPR, IRD, Fudan University, Shanghai, China, WHO Collaborating Center for Research in Human Reproduction, Shanghai, China

  • De-kun Li,

    Affiliation Division of Research, Kaiser Foundation Research Institute, Kaiser Permanente Northern California, Oakland, and Department of Health Research and Policy, School of Medicine, Stanford University, Stanford, California, United States of America

  • Fen Yang,

    Affiliations Key Laboratory of Reproduction Regulation of NPFPC, SIPPR, IRD, Fudan University, Shanghai, China, WHO Collaborating Center for Research in Human Reproduction, Shanghai, China

  • Hongjie Pan,

    Affiliations Key Laboratory of Reproduction Regulation of NPFPC, SIPPR, IRD, Fudan University, Shanghai, China, WHO Collaborating Center for Research in Human Reproduction, Shanghai, China

  • Tianqi Li,

    Affiliations Key Laboratory of Reproduction Regulation of NPFPC, SIPPR, IRD, Fudan University, Shanghai, China, WHO Collaborating Center for Research in Human Reproduction, Shanghai, China

  • Maohua Miao ,

    runshengli2007@163.com (RSL); miaomaohua@163.com (MHM)

    Affiliations Key Laboratory of Reproduction Regulation of NPFPC, SIPPR, IRD, Fudan University, Shanghai, China, WHO Collaborating Center for Research in Human Reproduction, Shanghai, China

  • Runsheng Li ,

    runshengli2007@163.com (RSL); miaomaohua@163.com (MHM)

    Affiliations Key Laboratory of Reproduction Regulation of NPFPC, SIPPR, IRD, Fudan University, Shanghai, China, WHO Collaborating Center for Research in Human Reproduction, Shanghai, China

  • Wei Yuan

    Affiliations Key Laboratory of Reproduction Regulation of NPFPC, SIPPR, IRD, Fudan University, Shanghai, China, WHO Collaborating Center for Research in Human Reproduction, Shanghai, China

Abstract

Environmental BPA exposure has been shown to impact human sperm concentration and motility, as well as rodent spermatogenesis. However, it is unclear whether BPA exposure is associated with alteration in DNA hydroxymethylation, a marker for epigenetic modification, in human sperm. A genome-wide DNA hydroxymethylation study was performed using sperm samples of men who were occupationally exposed to BPA. Compared with controls who had no occupational BPA exposure, the total levels of 5-hydroxymethylcytosine (5hmc) increased significantly (19.37% increase) in BPA-exposed men, with 72.69% of genome regions harboring 5hmc. A total of 9,610 differential 5hmc regions (DhMRs) were revealed in BPA-exposed men relative to controls, which were mainly located in intergenic and intron regions. These DhMRs were composed of 8,670 hyper-hMRs and 940 hypo-hMRs, affecting 2,008 genes and the repetitive elements. The hyper-hMRs affected genes were enriched in pathways associated with nervous system, development, cardiovascular diseases and signal transduction. Additionally, enrichment of 5hmc was observed in the promoters of eight maternally expressed imprinted genes in BPA-exposed sperm. Some of the BPA-affected genes, for example, MLH1, CHD2, SPATA12 and SPATA20 might participate in the response to DNA damage in germ cells caused by BPA. Our analysis showed that enrichment of 5hmc both in promoters and gene bodies is higher in the genes whose expression has been detected in human sperm than those whose expression is absent. Importantly, we observed that BPA exposure affected the 5hmc level in 11.4% of these genes expressed in sperm, and in 6.85% of the sperm genome. Finally, we also observed that BPA exposure tends to change the 5hmc enrichment in the genes which was previously reported to be distributed with the trimethylated Histone 3 (H3K27me3, H3K4me2 or H3K4me3) in sperm. Thus, these results suggest that BPA exposure likely interferes with gene expression via affecting DNA hydroxymethylation in a way partially dependent on trimethylation of H3 in human spermatogenesis. Our current study reveals a new mechanism by which BPA exposure reduces human sperm quality.

Introduction

Bisphenol A (BPA) is an endocrine disrupting chemical used mainly in epoxy resin and polycarbonate plastic industry. Over 90 percent of the study subjects in the US had detectable total urinary BPA levels [1]. Human studies have shown that BPA exposure is associated with sexual function as well as sexual hormones among males [2, 3]. We identified an inverse association between total urinary BPA levels and semen concentration, total sperm count, sperm vitality and sperm motility among factory workers [4]. The inverse association was also reported in men recruited from different infertility clinics [57].

Male rodents exposed to various levels of BPA displayed several reproductive and developmental impacts, including decreased serum testosterone levels and sperm quality [8]. For example, BPA has been shown to impact sperm motility in vivo and in vitro [9, 10]. In vitro exposing sperm to BPA also impacted fertilization and early embryonic development [9]. Among the known epigenetic regulatory mechanisms, DNA methylation, via downregulating gene expression [11], has been demonstrated to play an important role during mammalian spermatogenesis [12, 13], as well as in maintenance of human sperm quality[14]. Diverse effects of BPA on different cells/tissues have been linked to abnormal pattern of DNA methylation [1517]. For example, paternal BPA exposure of mice resulted in a global alteration of DNA methylation in offspring’s sperm which is associated with an impaired spermatogenesis [17] and heart problems[18] observed in offspring.

We previously reported that sperm LINE-1 methylation level was significantly lower in BPA-exposed workers compared to unexposed workers [19]. However, little is known regarding the mechanisms by which BPA modulates DNA methylation.

5-hydroxymethylcytosine (5hmc) is oxidized from 5-methylcytosine (5mC) by TET family of proteins, and has been recently recognized as a product in intermediate process of DNA demethylation [20, 21]. Importantly, dynamics of 5hmC has been detected during mouse spermatogenesis, and enrichment of 5hmc is associated with active gene transcription [22]. Given LINE-1 methylation has been treated as a marker of global DNA methylation, therefore the hypomethylation of LINE-1 present in sperm of BPA-exposed workers [19] is consistent with the hypothesis that BPA exposure promotes DNA hydroxymethylation in sperm in these workers. In the present study, we examined whether global level of DNA hydroxymethylation is different between BPA-exposed workers and unexposed controls.

Results

DNA hydroxymethylomes of human sperm

In order to obtain a general overview of 5hmC in the matured human spermatozoa and the impact of BPA exposure to adult sperm, we collected sperm samples from 26 adults as controls, and samples from 30 workers in a factory manufacturing BPA. The sperm from control and BPA-exposed men were pooled separately for 5-hMeDIP-seq. DNA pooling approach facilitated whole-genome 5hmc analysis of sperm which usually required large quantities of DNA. Compared with data generated from individual samples, the analysis based on DNA pools was validated to provide an accurate and reliable quantitative estimate of DNA methylation, showing highly significant correlations between individual samples and pools (95% bootstrapped confidence intervals: 0.94 to 0.96) [23, 24]. Pooled samples have been widely used in epigenetic studies, including genome-wide analysis of DNA methylation in human amnion, peripheral blood mononuclear cells, cell-free DNA of cancer, etc., [2527].

Globally increased 5hmc levels in sperm from BPA-exposed men.

We profiled the 5hmc distribution in human sperm cells. Based on the hMeDIP-seq data, we calculated the short read coverage at genome-wide 100-bp bins using ‘MEDIPS’[28] and transformed them into RPKM.

DNA 5hmc in sperm from BPA-exposed men was detected in 72.69% of genome (22,500,965 regions), while it was 60.89% of genome (18,849,185 regions) in sperm from the men without occupational BPA exposure (hereafter referred to as controls), indicating that BPA exposure significantly raises global level of DNA hydroxymethylation by 19.3% in human sperm.

We then examined the 5hmc rates on each chromosome (Fig 1A), and found that in both BPA-exposed and the control samples, most chromosomes exhibited high 5hmc rates (>50%), while Chromosome Y showed the lowest 5hmc level (22.72% in control group and 27.80% in BPA-exposed group). Strikingly, we observed that BPA exposure evenly increased ~19% of 5hmc level in autosomes, while in sex chromosomes, the 5hmc level increased 25–26% (Fig 1B), suggesting that 5hmc rates in sex chromosomes are more sensitive to BPA exposure. Relative to autosomes, the increased 5hmc in sex chromosomes showed a bit bias to 5’-UTR (15.3% in sex chromosomes vs 14.1% in autosomes) and LINE (12.6% in sex chromosomes vs 10.4% in autosomes).

thumbnail
Fig 1. 5hmc distribution in sperm genome.

(a). 5hmc region ratio on the 19 autosomal chromosomes and the X and Y sex chromosomes. 100% cover ratio is defined as covering all the human genome. (b). 5hmc increase ratio on the 19 autosomal chromosomes and the X and Y sex chromosomes, where hMeDIP-seq data covered genome size in each chromosome was defined as 1. (c) The ratio of 5hmc region in each gene body. 100% cover ratio is defined as covering all of the elements of a particular type in the human genome. (d) Proportion of 5hmC regions with different CpG densities. (e) Densities of 5hmC (counts per kilobases per million reads) in each 100bp regions with different CpG densities.

https://doi.org/10.1371/journal.pone.0178535.g001

The increase of 5hmc level was also reflected in each gene body (Fig 1C and S1 Table), with a maximum increase of 26% in 5’UTR region (35.64% of 5’UTR region in control and 45.26% in BPA-exposed sperm was covered by hMeDIP-seq data) and minimum increase of 18% in intron (65.1% in control and 77.19% in BPA-exposed sperm).

5hmC and CpG density.

A total of 28,217,009 CpGs sites are present in human genome. MeDIP-seq data of human sperm had revealed that ~60% of all CpGs in the human genome were methylated, and promoters display an inverse correlation between CpG density and methylation [29]. Since Ten-eleven translocation (Tet) family proteins can enzymatically convert 5mc to 5hmC, we postulated that 5hmc levels in the sperm genome correlated with CpG density. The hMeDIP-seq data coverage analysis shows a medium CpG coverage saturation, with 16.2 million CpGs (57.3%) in the control samples and 19.8 million CpGs (70.2%) in the BPA-exposed samples, were sequenced to a depth of at least 1x coverage (S1 Fig). The BPA-exposed sperm harbored about 3.6 million (or higher by 22.5%) more hydroxymethylated CpGs than the control sperm, suggesting that BPA exposure affects spermatogenesis via disturbing DNA methylation-hydroxymethylation balance.

5hmC-enriched regions were associated with moderate CpG density (1–5 CpG per 100bp) in mouse embryonic stem cells [30]. Our study exhibited a consistent result, with 68.9% of moderate CpG density region covered by hMeDIP-seq data in normal sperm, while the coverage for low CpG density region (<1 CpG per 100bp) and high CpG density region (>5 CpG per 100bp) were 51.1% and 50.7%, respectively (Fig 1D). In low and moderate CpG density regions, the 5hmc density increased with CpG density (Fig 1E); while in high CpG density regions, the 5hmc density showed an inverse correlation with CpG density.

5hmc in CpG islands (CGI).

A total of 28,691 CGIs with average length of 761 bp and average number of CpG sites of 69 were observed in human genome (from UCSC Genome Browser “CpG Islands” track). Among them, 20,078 and 16,319 CGI were captured by hMeDIP-seq in BPA-exposed sperm (69.98%) and the control sperm (56.88%), respectively (S2 Table), indicating that the 5hmc level of CGI in sperm from BPA-exposed men is 23.0% higher than the control, in line with the result that the 5hmc level of CpGs in sperm with BPA exposure was higher by 22.5% than the control.

5hmc in promoters.

A total of 29,310 candidate promoter regions (2-kb-upstream of transcription start site (TSS)) of human genes were used to analyze the hydroxymethylation status of promoters. 10,993 promoters in sperm with BPA exposure (37.5%) with average RPKM of 28.5 and 10,538 in the control (35.9%) with average PRKM of 28.0 were covered by our hMeDIP data (S3 Table), and the difference is not statistically significant. It suggests that BPA-associated upregulated rates of 5hmc had no bias in promoter regions. In addition, 5hmc level also showed an even distribution along the whole promoter region.

Identification of differentially hydroxymethylated regions (DhMRs)

We next analyzed the hydroxymethylated regions (hMRs) whose 5hmc rates are BPA-responsive in sperm, and named them differentially 5hmc regions (DhMRs). P-values were calculated using edgeR by comparing the RPKM of the sperm with BPA exposure and control samples within each of the 100-bp windows. DhMRs were identified by filtering for windows associated with a P-value ≤0.005. This step filtered 9,610 DhMRs between sperm with BPA exposure and control, with an average RPKM 7.2 in sperm with BPA exposure and 1.8 in control (Fig 2A and S4 Table). In addition, we noticed that chromosome 19 exhibited the highest DhMRs density. Chromosome 19 is extremely gene dense relative to other human chromosomes [31], indicating that DhMRs might be associated with gene distribution. Genomic region enrichment analysis of these DhMRs indicated that they were mainly located in intron regions (Fig 2B).

thumbnail
Fig 2. Summary of DhMR location distribution in sperm genome.

a. Atlas of DhMRs distribution in different chromosomes. Red lines in inner circle represent the position of hyper-hMRs and blue lines represent the position of hypo-hMRs. The length unit of the outer circle is 100Mb. b. DhMRs in gene body. c. 5hmc density in repeats. d. DhMRs affected repeats density.

https://doi.org/10.1371/journal.pone.0178535.g002

Distribution of DhMRs in the repetitive sequences.

We analyzed the distribution of BPA-associated 5hmc in repeat regions, and found it showed a slight bias towards repeats. That is, 68.2% of repeats in the control while 80.2% of repeats in sperm with BPA exposure overlapped with 5hmc region (Fig 2C and S5 Table). We also examined the distribution of DhMRs on repetitive sequences. In general, LINEs were enriched with DhMRs (Fig 2D). Increased expression of LINE-1 is closely correlated with the decreased methylation in LINE-1 5’UTR [32]. Since the 5hmc level in 5’UTR was significantly improved in BPA-exposed sperm (Fig 3C), we speculated that the corresponding methylation level was decreased, thus BPA-related hyper-hMRs might active LINE-1, further damaging sperm DNA.

thumbnail
Fig 3. 5hmc distribution and BPA affection in promoter and gene body of sperm expressed or repressed genes.

a. Composition of 5hmc in promoter. b. Relative to control, the 5hmc level of each gene body in BPA-exposed sperm.

https://doi.org/10.1371/journal.pone.0178535.g003

Imprinted genes associated with exposure to BPA.

Imprinted genes show monoallelic or highly biased expression according to the parental origin of the allele and were controlled by DNA methylation differently in the parental germline [33]. A total of 210 human imprinted genes collected from NCBI and http://www.geneimprint.com/databases/ were analyzed in our study. We found that eight maternally expressed imprinted genes had their promoters hydroxymethylated in BPA-exposed sperm (S6 Table). Four of them (FOXF1, SALL1, CDKN1C and HOXC4) had been found hypomethylated in sperm with low motility [34] (Table 1), providing supportive evidence for our hypothesis that BPA-responsive DNA hydroxymethylation in the promoters of imprinted genes may have resulted in a low sperm motility that was previously observed in BPA-exposed men[4].

thumbnail
Table 1. Known sperm function associated genes affected by BPA.

https://doi.org/10.1371/journal.pone.0178535.t001

Functional analysis of DhMRs.

Compared with the controls, we found 8,670 hyper-hMRs and 940 hypo-hMRs in BPA-exposed sperm (S7 and S8 Tables), affecting 2,008 genes (1,870 by hyper-hMRs and 274 by hypo-hMRs). These genes showed an uneven distribution on different chromosomes, with 11.8% of genes on chr18 (46/339) while only 2.6% of genes on sex chromosomes (29/1121 on chrX and 3/113 on chrY) affected by BPA exposure. As shown in Fig 1B, BPA exposure caused more 5hmc level increase (25–26%) in sex chromosomes than in autosomes (~19%), but DhMRs affected gene ratio on sex chromosomes were much lower than on autosomes. This might be caused by low gene density on sex chromosomes, as is known that human X chromosome has roughly half the gene density of human autosomes [35]. Interestingly, the hyper-hMRs affected genes were enriched with functional annotation terms related to cell adhesion (‘cell-cell adhesion’, ‘cell adhesion’ and ‘homophilic cell adhesion’, P value<1E-4), exocytosis (GO:0006887), cell migration, protein kinase A regulatory subunit binding, etc (S9 Table). It seemed that hyper-hMRs were associated with the expression of germ cell-interaction genes. In addition, the genes with the hyper-hMRs were enriched in seven pathways, mainly associated with nervous system (four pathways), development, cardiovascular diseases and signal transduction (S10 Table). Hypo-hMRs affected genes were enriched in 12 GO terms and two KEGG pathways, including oxidative phosphorylation (S11 Table).

Analysis of association of DhMRs with sperm-expressed and repressed genes

The studies based on different techniques including microarray [36], SAGE [37] and RNA-seq [38, 39] have proved that human spermatozoa contain a complex repertoire of mRNAs and small RNAs, though their expression level is much lower than testes cells. By using RNA-seq, Sendler et al. identified 6,709 expressed genes (RPKM>1) in human sperm [39], with 726 genes considered as highly expressed genes (RPKM>50) in human sperm (S12 Table), while 523 abundantly expressed genes in testes (PRKM>25) were validated not expressed in sperm (RPKM<1) (S13 Table), which is defined as sperm-repressed genes.

Given 5hmC in the mouse genome is crucial for regulation of gene expression during the differentiation of spermatogenic cells[22], we speculated that it also modifies gene expression before human sperm maturation, considering that sperm are known to be transcriptionally silent. Therefore, we compared DhMRs in the genes with different expression in sperm. 678 genes sperm-expressed were observed to have hyper-hMRs, while 84 sperm-expressed genes were observed to have hypo-hMRs (Fig 3), indicating that 11.4% (762/6,709) of sperm-expressed genes have altered status of DNA hydroxymethylation as a result of BPA exposure. In contrast, status of DNA hydroxymethylation of only 6.85% (2,008/29,310) of total genes is responsible to BPA exposure, as mentioned above (S5, S7 and S8 Tables). Therefore, our analysis strongly suggests the presence of trend that BPA preferentially affects level of DNA hydroxymethylation of sperm-expressed genes (FDR = 2.05E-56).

Association of the 5hmc rates in the promoters and gene bodies with level of gene expression was next analyzed. We checked 5hmc level in promoter and gene body of sperm-expressed and repressed genes[39]. The 3.0k-bp fragment of each gene (from -2.5k to 0.5k, relative to TSS) was subdivided into six bins, with 500-bp for each bin to calculate the 5hmc level (Fig 3A). The 5hmc levels of promoter regions in sperm-expressed genes were significantly higher than in sperm-repressed genes (Fig 3A), indicating that enrichment of 5hmc in promoter plays a critical role for gene expression. This is in accordance with reported observations that 5hmC levels in the promoters are positively correlated with gene expression [22, 40].

We observed the presence of BPA-associated enrichment of 5hmc in every analyzed part of genes. Compared to sperm-expressed genes, the sperm-repressed genes showed a lower enrichment of 5hmc in gene bodies including 5’UTR, exon, intron and 3’UTR (Fig 3B), consistent with the above-mentioned observation that sperm-expressed genes are more sensitive to BPA exposure.

As mentioned above, BPA-associated alteration in 5hmc rate had no bias to the promoter region. However, we did find hyper-hMRs in the promoters of 58 genes. 24 out of the 58 genes harbored no detectable 5hmc in their promoters in the control, but significantly increased 5hmc level after BPA exposure (S14 Table). For example, the 5hmc rate in the promoter of AChE gene, which encodes acetylcholine hydrolyzing enzyme acetylcholinesterase, significantly raised from 0 in the control sperm to 5–6 (RPKM) in BPA-exposed sperm. AChE is sperm-expressed gene[39]. Recently, we detected an inverse correlation between the 5hmc rate in the promoter of AChE gene and sperm motility in the subfertile male population (submitted). Considering that sperm AChE activity was negatively related to the development of sperm motility during sperm maturation[41], the BPA-associated enrichment of 5hmc in the promoter of AChE gene suggests the presence of upregulation of AChE expression in the sperm from BPA-exposed men, eventually making an contribution to development of sperm with poor motility observed in BPA-exposed men[41].

We observed that 38 sperm-repressed genes [39] were enriched with hyper-hMRs, with seven of them in the promoter, 5’UTR, exon or TTS (S13 Table). One of them is PPP2R3C, which encodes Protein phosphatase 2A, regulatory subunit B, gamma. It had been found to be associated with globozoospermia [42].

BPA exposure was previously shown to cause DNA damage in human sperm[7] and murine sperm[43]. The transcripts of some DNA damage-associated genes are abundant in human sperm [39]. We found that their 5hmc levels are upregulated by BPA exposure (S12 Table). For example, a hyper-hMRs was present in the promoter of CHD2 (chromodomain helicase DNA-binding protein 2), which was proposed to play an important role in the DNA damage response and genome stability maintenance [44]. SPATA12 (spermatogenesis associated protein 12) has been revealed interacting with CHD2 and playing a role in DNA damage signaling [45], and we also found the hyper-hMRs in the 1st intron of SPATA12 (S12 Table). Additionally, the hyper-hMRs was observed to be present in its 5’UTR of DNA mismatch repair gene MLH1 (S13 Table). Taken together, the results suggest upregulated expression of the sperm-expressed genes via raised DNA hydroxymethylation participates in BPA-associated DNA damage and repairing during human spermatogenesis (Table 1).

On the contrary, nine sperm-expressed genes were observed to have hypo-hMRs in BPA-exposed sperm (S12 Table), and their expression might be accordingly down-regulated. RAD23B, highly expressed in the human adult testis and spermatozoa, plays an important role in spermatogenesis and sperm maturation by inhibiting Ub-dependent proteolysis [46]. CUL3 (Cullin-3) was expressed in male germline and required for male fertility [47]. The testis-specific HSPA1L was associated with the human sperm plasma membrane during spermatogenesis and its expression peaks in spermatids [48]. The list of genes also includes UBAP2 (ubiquitin associated protein 2) (Table 1), whose expression was reduced by 52-fold reduction expression in spermatozoa of oligozoospermic infertile men [49].

Our analysis identified hyper-hMRs in the gene bodies of 131 genes (S15 Table) whose expression was reported to be absent both in human testis and sperm[39]. One of them is NOS. Testicular NO synthase (NOS) activity is higher in old rats than younger one, and were proposed to be involved in NO-mediated intrinsic pathway signaling in age-related increase in germ-cell apoptosis in male rats[50]. Considering the hyper-hMRs in gene bodies is correlated with upregulated gene expression, therefore, BPA-associated hyper-hMR in the NOS gene could affect spermatogenesis via a raised rate in germ-cell apoptosis.

Analysis of association of DhMRs with H3 methylation in sperm-expressed genes

The vast majority of mammalian sperm DNA is bound protamines, while nucleosome-binding sites constituting 2.9% of the paternal genome had been revealed in the human sperm genome [51]. We next analyzed whether BPA-associated alterations of DNA hydroxymethylation is associated with Histone H3 methylation. Based on ChIP-Seq analysis, two types of H3 Lysine 4 methylation (H3K4me3 and H3K4me2: two activating histone modification) were reported to exist in 5,455 and 9,073 genes, respectively in human sperm, while H3 Lysine 27 trimethylation (H3K27me 3: a repressive histone modification) is only present in 3,911 gene [52, 53]. Keeping in line with the reports, our analysis showed 43.9% (2,390/5,455) and 40.9% (3,708/9,073) of genes binding with H3K4me2 and H3K4me3, respectively, are expressed in sperm (Fig 4A and 4B). In contrast, only 20.1% (788/3911) of H3K27me3-binding genes is expressed in sperm (Fig 4C).

thumbnail
Fig 4. Relationship of 5hmc, histone binding and expression of sperm genes.

a. relationship of H3K4me2 binding genes with DhMR genes and gene expression. b. relationship of H3K4me3 binding genes with DhMR genes and gene expression. c. relationship of H3K27me3 binding genes with DhMR genes and gene expression. d. relationship of H3K4me3/H3K27me3 binding genes with DhMR genes and gene expression.

https://doi.org/10.1371/journal.pone.0178535.g004

Association of DhMRs with the three types of H3 methylation was analyzed accordingly. Contrast to 6.85% (2,008/29,310) genes in sperm genome was associated with BPA exposure in sperm, a significantly higher percentage (from 9.86–12.1%) of H3 (bearing one of three types of lysine methylation) binding genes has DhMRs (Fig 4A–4C), indicating that H3 methylation raises efficiency of BPA-modifying DNA hydroxymethylation in sperm. Our analysis showed that 12.1% (474/3,911) of H3K27me3-binding genes were affected by DhMRs, while they were detected in 10.5% (573/5,455) H3K4me2- and 9.86% (895/9,073) of H3K4me3-binding genes, respectively (Fig 4A–4C). The difference is statistically significant (P-value between H3K4me2 and H3K27me3: 0.0038; P-value between H3K4me3 and H3K27me3: 2.529e-05), Similarly, in the 2,478 genes enriched both with H3K4me3 and H3K27me3 in human sperm[52] (or defined as bivalent genes), we observed that 12.3% (306/2,478) of bivalent genes harbored DhMRs (Fig 4D). These results strongly suggested that BPA has a potential to preferentially raise level of 5hmc in the genes occupied with H3K27me3.

We further performed an analysis similar to the one above, focusing on whether the three types of H3 methylation differentially affect association of DhMRs with sperm profile of gene expression. Our analysis revealed that 44.6% (895/2,008) of the genes with BPA-responsive DhMRs bind H3K4me3, while only 28.5% (573/2,008), 23.6% (474/2,008) and 15.2% (306/2,008) of them binds with H3K4me2, H3K27me3 and H3K4me3/H3K27me3, respectively (The p-value between H3K4me3 and H3K4me2 was 2.2E-16, between H3K4me3 and H3K27me3 was 2.2E-163H) (Fig 4A–4D). The analysis demonstrated that the sperm-expressed genes are more sensitive to BPA for their alterations in DNA hydroxymethylation when they are enriched with H3K4me3 than other examined styles of H3 methylation.

Analysis of BPA-associated 5hmc via qPCR

We also compared the rates of 5hmc in sperm from men exposed to BPA with those from the men without significant BPA exposure. Four genes (AChE, ATP2A1, PPP2R3C and LINE-1) were detected to have hyper-DhMRs in the regions of promoter, intron or 5’UTR, as indicated in Fig 5A, and were selected for the analysis.

thumbnail
Fig 5. BPA raised 5hmc rate in sperm genes analyzed by qPCR.

a. Schematic representation of location of the CCGG loci in the 5UTR /promoters/ gene body of LINE-1, AChE, ATP2A1 and PPP2R3C, respectively. The red line indicates the region where enrichment of BPA-associated 5hmc was detected through high-throughput sequencing, and it is 200-bp long in LINE-1, AChE, ATP2A1, and 300-bp long in PPP2R3C, respectively. The green line represents the location of qPCR product. b and d. The 5hmc rates of four genes were measured using the pooled sperm samples from BPA-exposed men and the control men. The sperm DNA was treated with glucosylation followed by CCGG loci-dependent MspI/HpaII digestion (b), or was immuno-precipitated with the 5hmc antibody (d), and then used as the template of qPCR. The data were obtained from at least three independent experiments, and the values are shown as the mean ± SEM; * p< 0.05. c. The 5hmc rates were measured using individual samples. Five BPA-exposed samples and five control samples were randomly chosen for measurement of 5hmc rate. The values are shown as the mean ± SEM; * p< 0.05.

https://doi.org/10.1371/journal.pone.0178535.g005

The 5hmC rates in the CpG sites were measured by qPCR after sperm DNA was applied to treatment with glucosylation followed by CCGG loci-dependent MspI/HpaII digestion[54]. We first used the pooled DNA samples for the measurement. qPCR assays across the selected CCGG sites showed a varied ranges of 5hmC rates, i.e., from 0.85% to 5.38% in the four genes, with the lowest rate detected in the AChE promoter. Moreover, we observed significant higher 5hmc rates in all the four genes in sperm from men exposed to BPA than the control (Fig 5B).

Considering the variations of 5hmc rates among different individual samples were not reflected in the pooled samples-based analysis, we repeated the above measurement in five randomly-selected sperm samples with BPA exposure and five randomly-selected control samples (Fig 5C). Again, our results showed BPA exposure raised the 5hmc rates in these genes with an efficiency comparable to that observed in the pooled samples.

Given our data were generated based on hMeDIP assay in the present study, we finally compared levels of 5hmc in these genes in BPA-exposed samples with the control by using the pooled DNA samples that were immuno-precipitated with the 5hmc antibody. As expected, a BPA-associated rise in the 5hmc rate was detected in all the four genes (Fig 5D). Importantly, by using two different techniques, we observed BPA exposure raised the 5hmc rate with a similar efficiency in the four genes, respectively (Fig 5B and 5D).

Discussion

We reported in the present study that global level of DNA hydroxymethylation is higher by 19.6% in sperm from BPA-exposed men than from men without occupational exposure to BPA. The BPA-responsive alteration in DNA hydroxymethylation is characterized by the following two features: (1), 5hmc rates in sex chromosomes are more sensitive than autosomes to BPA exposure. It is currently unknown whether it is associated with BPA’s weak estrogen activity; and (2), LINES has been identified as the repetitive elements that are more sensitive to BPA exposure than any other genes and the repetitive elements. Of the genes affected by BPA exposure, the retrotransposons are particularly intriguing. They are thought to be maintained in a transcriptionally silent state by DNA methylation[55], but they are activated via BPA-induced DNA demethylation [56]. Positive correlation between expression and 5hmC density of retrotransposons during spermatogenesis was recently reported, highlighting the important role of 5hmc in regulating of retrotransposon activity [22]. It is noteworthy, therefore, that BPA raises 5hmc rates in retrotransposons including LINE-1 of human sperm, consistent with the hypothesis that BPA-responsive hypomethylation of LINE-1 detected in human sperm[19] is mediated by enhanced DNA hydroxymethylation. Given more than 100 LINE-1-mediated insertions which result in genetic diseases have been reported [57] and derepression of LINE-1 along with complete male sterility is linked to azoospermia [58], upregulated 5hmc rate in LINE-1 likely raises genomic instability, eventually making a contribution to reduced sperm concentration in the men exposed to BPA.

Besides the repetitive elements, approximately 2,000 sperm genes were observed to have changed their 5hmc rates as a result of BPA exposure. In a rational manner, we hypothesized that BPA exposure modifies the epigenetic setting of genes at different stages of spermatogenesis, and altered expression of some of these genes impairs spermatogenensis via distinct mechanisms. It has been recently recognized that 5hmC distribution undergoes an interesting dynamic change in the entire genome during mouse spermatogenesis[22]. Consistent with the report, levels of expression of Tet1-3 are positively correlated with sperm quality [59]. Importantly, concerning ICSI-outcome, the lowest levels of TET1-3 mRNAs in sperm were found in the non-pregnant group, while increased TET2 in sperm was significantly associated with pregnancy[59]. Globally altered expression of 5hmc has been reported to abnormal sperm of humans [60]. These reports indicated that an optimal pattern of genomic DNA hydroxymethylation is pivotal for maintenance of seminal quality and male fertility. Thus, genome-wide rise of the 5hmC marker is a new epigenetic feature of reproductive toxicity of BPA, and could be potentially regarded as a biomarker for BPA exposure.

The mechanism(s) underlying globally raised 5hmC in sperm responding to BPA exposure needs to be investigated and understood. The global level of 5hmC is reduced as male mouse meiosis starts [22], and BPA exposure disrupted meiosis progression [61]. These results suggest the BPA represses expression of TET1-3, which are the key regulators of global level of 5hmc during spermatogenesis. The mechanistic process may be investigated using in vitro approaches in future studies. However, multiple factors may be involved in the control of sperm genome-wide DNA hydroxymethylation in sperm, as suggested by an opposite changes in level of 5hmc [62] and TETs [59] was observed in sperm from aging men. Methyl-CpG-binding protein 2 (MeCP2) has been identified as the major 5hmC-binding protein in the brain [63]. Its expression also improves the development of mouse somatic cell nuclear transfer embryos, via a 5hmc-upregulating mechanism[64]. Interestingly, BPA exposure obviously upregulated expression of MeCP2, given abundance of MeCP2-expressing neurons significantly increased[65]. Therefore, it may be speculated that BPA raises genome-wide level of sperm 5hmc in the way depending on MeCP2. It is worthy for future study.

Our analysis demonstrated that the 5hmc levels in gene bodies and promoters were correlated with their levels of expression in human sperm, respectively, consistent with the previously reported observation that the 5hmc level positively regulated gene expression in eight consecutive types of mouse spermatogenic cells[22]. In a rational manner, we speculate that BPA exposure modifies the epigenetic setting of 2,008 genes at different stages of spermatogenesis. We observed that BPA preferentially affects the 5hmc levels of sperm-expressed genes, and such a modification probably depends on H3 methylation, given that it raised BPA-altering 5hmc in sperm. Additionally, studies have revealed a close association of 5hmc with histone modifications. For example, enrichment of 5hmc exists in the genes bound with active (H3K4me2 and 3) and repressive (H3K27me3) H3 modification [66]. Via the analysis based on co-immunolabeling of 5hmC and H3K27me3, Haffner et al. also reported that global level of 5hmc and H3K27me3 are tightly co-regulated during hierarchical differentiation in adult tissues [67]. In the present study, we observed that BPA has a potential to preferentially raise the enrichment of 5hmc in the genes occupied with H3K27me3. Enhancer of Zeste Homolog 2 (EZH2) is known to be the catalytic subunit of the polycomb repressive complex 2 (PRC2), a histone methyltransferase catalyzing trimethylation of H3K27. The distribution of Ezh2 was detected in spermatogonial stem cells [68] and the round spermatids [69]. Although it remains unknown whether BPA affects EZH2 expression during spermatogenesis, it did so in both cultured breast cancer cells and mouse mammary [70, 71]. Taken together, these results suggest BPA raises 5hmc rates during human spermatogenesis in the way dependent on EZH2-based H3K27me3.

Altered pattern of DNA methylation in sperm has been proposed to play an important role in a transgenerational toxic effect caused by paternal exposure of BPA [17, 72] and other environmental factors. Reprogramming of the paternal genome upon fertilization involves genome-wide oxidation of 5-methylcytosine-mediated by manteral TET3 [41]. Aberration in status of paternal DNA hydroxymethylation has been shown to be linked to poor preimplantation embryo development [73, 74]. For example, human [60] and mouse [22] round spermatids have a globally-altered pattern of 5hmc, compared with mature sperm, while active DNA demethylation normal-sized embryos were detected in a much reduced rate in the embryo with injection of round spermatid[73]. These reports suggest a close relationship between paternal alteration of 5hmc and impaired development [75]. Accordingly, our present study also indicates a potential risk that paternal alteration in DNA hydroxymethylation participates in a transgenerational toxic effect of BPA.

It should be noted that urine BPA level in the BPA-exposed men in the present study is quite high compared with those who were exposed through environmental sources. Further investigation is needed to examine whether the above observed association also exists at lower BPA exposure levels.

Conclusions

We reported a genome-wide rise in 5hmc rate in the sperm from BPA-exposed men in the present study. To our knowledge, this is the first report regarding the effect of BPA on DNA hydroxymethylation in human tissue. Our results suggest that BPA may affect gene expression and reduce genomic stability via altered rate of 5hmc, eventually resulting a previously reported low concentration and motility of sperm observed in BPA-exposed men [4]. Therefore, our study provides new evidence for a potential underlying mechanism of sperm toxicity of BPA exposure.

Materials and methods

Ethics statement

The study was reviewed and approved by the ethics committee board of Shanghai Institute of Planned Parenthood Research (IRB00008297). All participants gave written informed consent before participation in the study.

Semen samples

Collection of semen samples from the workers with BPA exposure and from the control men without significant BPA exposure was previously described [4]. The semen samples of exposed group were from 30 workers with an average age of 33.6 years, while those of the unexposed group were from 26 workers with an average age of 34.4 years. To be eligible as a subjects in the present study, the exposed workers had to be: 1) working in a BPA manufacturer or epoxy resin manufacturer; 2)aged between 20 to 49 years; 3) with a job title of packaging, chemical reaction, or material feeding, which were supposed to be exposed to high doses of BPA. While the unexposed workers are recruited if they were: 1) working in a factory in the same city as the exposed factory with no known occupational BPA exposure, nor known reproductive toxicants; 2) aged between 20 to 49 years. The median of urine BPA levels in the exposed and unexposed group were 251.33 ug/g creatine and 0 ug/g creatine, respectively. The percentages of smoking for the exposed and unexposed group were 43.3% and 42.3%.

hMeDIP

hmeDIP assays were performed on sperm DNA as described for 5mC as previously described by Mohn et al. [76] using commercial antibodies specific to 5hmC (Active Motif). Genomic DNA was extracted from the sperm from 30 BPA-exposed men and 26 controls, respectively. The BPA-exposed DNA and the control DNA were generated by the pooling of each sample with an equal amount of DNA. The pooled DNA was then sonicated to an average size of 200bp (range 100-500bp). DNA fragments were end-repaired, A-tailed and custom TruSeq adapters (non-methylated) were ligated. 1μg of adapter-ligated DNA was denatured and incubated with 1μl of 5hmC antibody at 4°C overnight. Antibody-DNA complexes were captured by protein A/G beads. The immunoprecipitated DNA was purified and then applied for sequencing analysis.

DNA hydroxymethylomes

The genome-wide maps of 5hmC were produced by performing hydroxymethyl-DNA immunoprecipitation followed by massively parallel sequencing with an Illumina Genome Analyzer (hMeDIP-seq). A total of 35,995,280 reads of 51 bp was obtained from control, and 44,263,854 reads of 51 bp was obtained from sperm with BPA exposure. Sequence quality control was performed using FASTQC.

SAMtools [77] was used to remove reads with low alignment score (Q < 20). Alignment to reference genome hg19 build downloaded from UCSC was performed using BWA (v0.5.8) [78] with default parameters. Then hMeDIP-seq data were processed using R-based Package ‘MEDIPS’ [28] with genomic window size set to 100 bp. For each sample, the aligned reads count at each window is transformed into reads per kilobase per million (RPKM) format. The region was defined as being covered if any part of the region was covered by a minimum of one sequencing reads. The saturation analyses resulted in genome-wide coverage saturation of 0.93 for control and 0.92 for case (S2 Fig).

Moreover, we tested the enrichment of CpG-rich short reads derived from the immunoprecipitation step, and found a relative enrichment for CpG-rich short reads from the control sample (1.05) and case sample (1.03) compared with the reference genome. Unexpectedly, the relative CpG enrichment is close to one for both samples (S16 Table).

DhMR analysis

P-values were calculated using edgeR by comparing the RPKM of the sperm with BPA exposure and control samples within each of the 100-bp windows. Differentially 5hmc regions (DhMRs) were identified by filtering for windows associated with a P-value ≤0.005. Then GO and pathway enrichment analysis of genes affected by hyper-hMRs or hypo-hMRs were separately performed based on EntrezID using DAVID [79].

Supporting information

S1 Fig. CpG coverage analysis.

The figure shows the results of hMeDIPS coverage analysis for control (a) and case (b).

https://doi.org/10.1371/journal.pone.0178535.s001

(TIF)

S2 Fig. Saturation analysis.

The figure shows the results of the saturation analysis of the MEDIPS package analyzing hMeDIP-seq data from control (a) and case (b).

https://doi.org/10.1371/journal.pone.0178535.s002

(TIF)

S1 Table. Distribution of 5hmc in gene body.

https://doi.org/10.1371/journal.pone.0178535.s003

(XLS)

S2 Table. Distribution of CpG islands covered by hMeDIP reads.

https://doi.org/10.1371/journal.pone.0178535.s004

(XLS)

S3 Table. Distribution of promoter region covered by hMeDIP reads.

https://doi.org/10.1371/journal.pone.0178535.s005

(XLS)

S4 Table. Differentially hydroxymethylated regions (DhMRs) between BPA-exposed men and controls.

https://doi.org/10.1371/journal.pone.0178535.s006

(XLS)

S5 Table. Distribution of hydroxymethylated regions in repeats.

https://doi.org/10.1371/journal.pone.0178535.s007

(XLS)

S6 Table. 5hmc level in promoter of human imprinted genes.

https://doi.org/10.1371/journal.pone.0178535.s008

(XLS)

S7 Table. Genes affected by hyper-hMRs in BPA-exposed population.

https://doi.org/10.1371/journal.pone.0178535.s009

(XLS)

S8 Table. Genes affected by hypo-hMRs in BPA-exposed population.

https://doi.org/10.1371/journal.pone.0178535.s010

(XLS)

S9 Table. Enrichment analysis of hyper-hMRs affected genes.

https://doi.org/10.1371/journal.pone.0178535.s011

(XLS)

S10 Table. Pathway enrichment analysis of hyper-hMRs affected genes.

https://doi.org/10.1371/journal.pone.0178535.s012

(XLS)

S11 Table. Enrichment analysis of hypo-hMRs affected genes.

https://doi.org/10.1371/journal.pone.0178535.s013

(XLS)

S12 Table. Highly expressed genes in sperm was affected by hMRs.

https://doi.org/10.1371/journal.pone.0178535.s014

(XLS)

S13 Table. Lowly expressed genes in sperm was affected by hyper-hMRs.

https://doi.org/10.1371/journal.pone.0178535.s015

(XLS)

S14 Table. Genes with promoter affected by DhMRs in the BPA-exposed population.

https://doi.org/10.1371/journal.pone.0178535.s016

(XLS)

S15 Table. Genes with gene body affected by hyper-hMRs in the BPA-exposed population.

https://doi.org/10.1371/journal.pone.0178535.s017

(XLS)

S16 Table. CpG enrichment analysis of the MEDIPS package.

https://doi.org/10.1371/journal.pone.0178535.s018

(XLS)

Acknowledgments

This work was supported by grants from the National Natural Science Foundation of China [Grant numbers. 81571495, 81501318 and 81270760], the National Basic Research Program of China [Grant numbers. 2014CB943104], Shanghai Municipal Committee of Science and Technology [Grant No. 15431902800] and Key Laboratory of Reproduction Regulation of NPFPC.

Author Contributions

  1. Conceptualization: RSL WY.
  2. Data curation: HJZ XYZ MHM.
  3. Formal analysis: HJZ XYZ DKL FY HJP.
  4. Funding acquisition: RSL MHM.
  5. Investigation: HJZ XYZ HJP TQL.
  6. Methodology: HJZ XYZ DKL MHM RSL.
  7. Project administration: RSL WY.
  8. Resources: MHM RSL WY.
  9. Software: HJZ XYZ TQL.
  10. Supervision: MHM RSL WY.
  11. Validation: HJZ XYZ DKL MHM RSL WY.
  12. Visualization: HJZ XYZ RSL.
  13. Writing – original draft: HJZ XYZ RSL.
  14. Writing – review & editing: HJZ DKL MHM RSL WY.

References

  1. 1. Cantonwine DE, Hauser R, Meeker JD. Bisphenol A and Human Reproductive Health. Expert review of obstetrics & gynecology. 2013;8(4). Epub 2013/11/05. PubMed Central PMCID: PMC3811157.
  2. 2. Hanaoka T, Kawamura N, Hara K, Tsugane S. Urinary bisphenol A and plasma hormone concentrations in male workers exposed to bisphenol A diglycidyl ether and mixed organic solvents. Occupational and environmental medicine. 2002;59(9):625–8. Epub 2002/09/03. PubMed Central PMCID: PMC1740362. pmid:12205237
  3. 3. Meeker JD, Calafat AM, Hauser R. Urinary bisphenol A concentrations in relation to serum thyroid and reproductive hormone levels in men from an infertility clinic. Environmental science & technology. 2010;44(4):1458–63. Epub 2009/12/25. PubMed Central PMCID: PMC2823133.
  4. 4. Li DK, Zhou Z, Miao M, He Y, Wang J, Ferber J, et al. Urine bisphenol-A (BPA) level in relation to semen quality. Fertility and sterility. 2011;95(2):625–30 e1-4. Epub 2010/11/03. pmid:21035116
  5. 5. Knez J, Kranvogl R, Breznik BP, Voncina E, Vlaisavljevic V. Are urinary bisphenol A levels in men related to semen quality and embryo development after medically assisted reproduction? Fertility and sterility. 2014;101(1):215–21 e5. Epub 2013/11/05. pmid:24182411
  6. 6. Lassen TH, Frederiksen H, Jensen TK, Petersen JH, Joensen UN, Main KM, et al. Urinary bisphenol A levels in young men: association with reproductive hormones and semen quality. Environmental health perspectives. 2014;122(5):478–84. Epub 2014/05/03. PubMed Central PMCID: PMC4014766. pmid:24786630
  7. 7. Meeker JD, Ehrlich S, Toth TL, Wright DL, Calafat AM, Trisini AT, et al. Semen quality and sperm DNA damage in relation to urinary bisphenol A among men from an infertility clinic. Reprod Toxicol. 2010;30(4):532–9. Epub 2010/07/27. PubMed Central PMCID: PMC2993767. pmid:20656017
  8. 8. Richter CA, Birnbaum LS, Farabollini F, Newbold RR, Rubin BS, Talsness CE, et al. In vivo effects of bisphenol A in laboratory rodent studies. Reprod Toxicol. 2007;24(2):199–224. Epub 2007/08/09. PubMed Central PMCID: PMC2151845. pmid:17683900
  9. 9. Rahman MS, Kwon WS, Lee JS, Yoon SJ, Ryu BY, Pang MG. Bisphenol-A affects male fertility via fertility-related proteins in spermatozoa. Scientific reports. 2015;5:9169. Epub 2015/03/17. PubMed Central PMCID: PMC4360475. pmid:25772901
  10. 10. Wang HF, Liu M, Li N, Luo T, Zheng LP, Zeng XH. Bisphenol A impairs Mature Sperm Functions by a CatSper-Relevant Mechanism. Toxicological sciences: an official journal of the Society of Toxicology. 2016. Epub 2016/04/30.
  11. 11. Kundu TK, Rao MR. CpG islands in chromatin organization and gene expression. Journal of biochemistry. 1999;125(2):217–22. Epub 1999/02/17. pmid:9990116
  12. 12. Calicchio R, Doridot L, Miralles F, Mehats C, Vaiman D. DNA methylation, an epigenetic mode of gene expression regulation in reproductive science. Current pharmaceutical design. 2014;20(11):1726–50. Epub 2013/07/31. pmid:23888966
  13. 13. Di Giacomo M, Comazzetto S, Saini H, De Fazio S, Carrieri C, Morgan M, et al. Multiple epigenetic mechanisms and the piRNA pathway enforce LINE1 silencing during adult spermatogenesis. Molecular cell. 2013;50(4):601–8. Epub 2013/05/28. pmid:23706823
  14. 14. Marques CJ, Costa P, Vaz B, Carvalho F, Fernandes S, Barros A, et al. Abnormal methylation of imprinted genes in human sperm is associated with oligozoospermia. Molecular human reproduction. 2008;14(2):67–74. Epub 2008/01/08. pmid:18178607
  15. 15. Yaoi T, Itoh K, Nakamura K, Ogi H, Fujiwara Y, Fushiki S. Genome-wide analysis of epigenomic alterations in fetal mouse forebrain after exposure to low doses of bisphenol A. Biochemical and biophysical research communications. 2008;376(3):563–7. Epub 2008/09/23. pmid:18804091
  16. 16. Doshi T, Mehta SS, Dighe V, Balasinor N, Vanage G. Hypermethylation of estrogen receptor promoter region in adult testis of rats exposed neonatally to bisphenol A. Toxicology. 2011;289(2–3):74–82. Epub 2011/08/11. pmid:21827818
  17. 17. Manikkam M, Tracey R, Guerrero-Bosagna C, Skinner MK. Plastics derived endocrine disruptors (BPA, DEHP and DBP) induce epigenetic transgenerational inheritance of obesity, reproductive disease and sperm epimutations. PloS one. 2013;8(1):e55387. Epub 2013/01/30. PubMed Central PMCID: PMC3554682. pmid:23359474
  18. 18. Lombo M, Fernandez-Diez C, Gonzalez-Rojo S, Navarro C, Robles V, Herraez MP. Transgenerational inheritance of heart disorders caused by paternal bisphenol A exposure. Environ Pollut. 2015;206:667–78. Epub 2015/09/01. pmid:26322593
  19. 19. Miao M, Zhou X, Li Y, Zhang O, Zhou Z, Li T, et al. LINE-1 hypomethylation in spermatozoa is associated with Bisphenol A exposure. Andrology. 2014;2(1):138–44. Epub 2013/12/03. pmid:24293158
  20. 20. Hill PW, Amouroux R, Hajkova P. DNA demethylation, Tet proteins and 5-hydroxymethylcytosine in epigenetic reprogramming: an emerging complex story. Genomics. 2014;104(5):324–33. Epub 2014/09/01. pmid:25173569
  21. 21. Branco MR, Ficz G, Reik W. Uncovering the role of 5-hydroxymethylcytosine in the epigenome. Nature reviews Genetics. 2012;13(1):7–13. Epub 2011/11/16.
  22. 22. Gan H, Wen L, Liao S, Lin X, Ma T, Liu J, et al. Dynamics of 5-hydroxymethylcytosine during mouse spermatogenesis. Nature communications. 2013;4:1995. pmid:23759713
  23. 23. Docherty SJ, Davis OS, Haworth CM, Plomin R, Mill J. Bisulfite-based epityping on pooled genomic DNA provides an accurate estimate of average group DNA methylation. Epigenetics & chromatin. 2009;2(1):3. PubMed Central PMCID: PMC2657899.
  24. 24. Gallego-Fabrega C, Carrera C, Muino E, Montaner J, Krupinski J, Fernandez-Cadenas I, et al. DNA methylation levels are highly correlated between pooled samples and averaged values when analysed using the Infinium HumanMethylation450 BeadChip array. Clinical epigenetics. 2015;7(1):78. PubMed Central PMCID: PMC4521379.
  25. 25. Kim J, Pitlick MM, Christine PJ, Schaefer AR, Saleme C, Comas B, et al. Genome-wide analysis of DNA methylation in human amnion. TheScientificWorldJournal. 2013;2013:678156. PubMed Central PMCID: PMC3590748. pmid:23533356
  26. 26. Legendre C, Gooden GC, Johnson K, Martinez RA, Liang WS, Salhia B. Whole-genome bisulfite sequencing of cell-free DNA identifies signature associated with metastatic breast cancer. Clinical epigenetics. 2015;7(1):100. PubMed Central PMCID: PMC4573288.
  27. 27. Sun L, Lin J, Du H, Hu C, Huang Z, Lv Z, et al. Gender-specific DNA methylome analysis of a Han Chinese longevity population. BioMed research international. 2014;2014:396727. PubMed Central PMCID: PMC4009103. pmid:24822201
  28. 28. Chavez L, Jozefczuk J, Grimm C, Dietrich J, Timmermann B, Lehrach H, et al. Computational analysis of genome-wide DNA methylation during the differentiation of human embryonic stem cells along the endodermal lineage. Genome research. 2010;20(10):1441–50. PubMed Central PMCID: PMC2945193. pmid:20802089
  29. 29. Down TA, Rakyan VK, Turner DJ, Flicek P, Li H, Kulesha E, et al. A Bayesian deconvolution strategy for immunoprecipitation-based DNA methylome analysis. Nature biotechnology. 2008;26(7):779–85. PubMed Central PMCID: PMC2644410. pmid:18612301
  30. 30. Wu H, D'Alessio AC, Ito S, Wang Z, Cui K, Zhao K, et al. Genome-wide analysis of 5-hydroxymethylcytosine distribution reveals its dual function in transcriptional regulation in mouse embryonic stem cells. Genes & development. 2011;25(7):679–84. PubMed Central PMCID: PMC3070931.
  31. 31. Grimwood J, Gordon LA, Olsen A, Terry A, Schmutz J, Lamerdin J, et al. The DNA sequence and biology of human chromosome 19. Nature. 2004;428(6982):529–35. pmid:15057824
  32. 32. He ZM, Li J, Hwa YL, Brost B, Fang Q, Jiang SW. Transition of LINE-1 DNA methylation status and altered expression in first and third trimester placentas. PloS one. 2014;9(5):e96994. PubMed Central PMCID: PMC4018393. pmid:24821186
  33. 33. Ishida M, Moore GE. The role of imprinted genes in humans. Molecular aspects of medicine. 2013;34(4):826–40. pmid:22771538
  34. 34. Pacheco SE, Houseman EA, Christensen BC, Marsit CJ, Kelsey KT, Sigman M, et al. 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. PubMed Central PMCID: PMC3107223. pmid:21674046
  35. 35. Bellott DW, Skaletsky H, Pyntikova T, Mardis ER, Graves T, Kremitzki C, et al. Convergent evolution of chicken Z and human X chromosomes by expansion and gene acquisition. Nature. 2010;466(7306):612–6. PubMed Central PMCID: PMC2943333. pmid:20622855
  36. 36. Ostermeier GC, Dix DJ, Miller D, Khatri P, Krawetz SA. Spermatozoal RNA profiles of normal fertile men. Lancet. 2002;360(9335):772–7. pmid:12241836
  37. 37. Zhao Y, Li Q, Yao C, Wang Z, Zhou Y, Wang Y, et al. Characterization and quantification of mRNA transcripts in ejaculated spermatozoa of fertile men by serial analysis of gene expression. Human reproduction. 2006;21(6):1583–90. pmid:16501037
  38. 38. Krawetz SA, Kruger A, Lalancette C, Tagett R, Anton E, Draghici S, et al. A survey of small RNAs in human sperm. Human reproduction. 2011;26(12):3401–12. PubMed Central PMCID: PMC3212879. pmid:21989093
  39. 39. Sendler E, Johnson GD, Mao S, Goodrich RJ, Diamond MP, Hauser R, et al. Stability, delivery and functions of human sperm RNAs at fertilization. Nucleic acids research. 2013;41(7):4104–17. PubMed Central PMCID: PMC3627604. pmid:23471003
  40. 40. Chen K, Zhang J, Guo Z, Ma Q, Xu Z, Zhou Y, et al. Loss of 5-hydroxymethylcytosine is linked to gene body hypermethylation in kidney cancer. Cell research. 2016;26(1):103–18. PubMed Central PMCID: PMC4816137. pmid:26680004
  41. 41. Wossidlo M, Nakamura T, Lepikhov K, Marques CJ, Zakhartchenko V, Boiani M, et al. 5-Hydroxymethylcytosine in the mammalian zygote is linked with epigenetic reprogramming. Nature communications. 2011;2:241. Epub 2011/03/17. pmid:21407207
  42. 42. Karaouzene T, El Atifi M, Issartel JP, Grepillat M, Coutton C, Martinez D, et al. Comparative testicular transcriptome of wild type and globozoospermic Dpy19l2 knock out mice. Basic and clinical andrology. 2013;23:7. PubMed Central PMCID: PMC4346239. pmid:25780569
  43. 43. De Flora S, Micale RT, La Maestra S, Izzotti A, D'Agostini F, Camoirano A, et al. Upregulation of clusterin in prostate and DNA damage in spermatozoa from bisphenol A-treated rats and formation of DNA adducts in cultured human prostatic cells. Toxicological sciences: an official journal of the Society of Toxicology. 2011;122(1):45–51. Epub 2011/05/04.
  44. 44. Nagarajan P, Onami TM, Rajagopalan S, Kania S, Donnell R, Venkatachalam S. Role of chromodomain helicase DNA-binding protein 2 in DNA damage response signaling and tumorigenesis. Oncogene. 2009;28(8):1053–62. PubMed Central PMCID: PMC2648865. pmid:19137022
  45. 45. Zhang Y, Yang L, Lin Y, Rong Z, Liu X, Li D. SPATA12 and its possible role in DNA damage induced by ultraviolet-C. PloS one. 2013;8(10):e78201. PubMed Central PMCID: PMC3799742. pmid:24205157
  46. 46. Huang X, Wang H, Xu M, Lu L, Xu Z, Li J, et al. Expression of a novel RAD23B mRNA splice variant in the human testis. Journal of andrology. 2004;25(3):363–8. pmid:15064313
  47. 47. Arama E, Bader M, Rieckhof GE, Steller H. A ubiquitin ligase complex regulates caspase activation during sperm differentiation in Drosophila. PLoS biology. 2007;5(10):e251. PubMed Central PMCID: PMC1976628. pmid:17880263
  48. 48. Naaby-Hansen S, Herr JC. Heat shock proteins on the human sperm surface. Journal of reproductive immunology. 2010;84(1):32–40. PubMed Central PMCID: PMC2898571. pmid:19962198
  49. 49. Montjean D, De La Grange P, Gentien D, Rapinat A, Belloc S, Cohen-Bacrie P, et al. Sperm transcriptome profiling in oligozoospermia. Journal of assisted reproduction and genetics. 2012;29(1):3–10. PubMed Central PMCID: PMC3252406. pmid:21989496
  50. 50. Molaro A, Hodges E, Fang F, Song Q, McCombie WR, Hannon GJ, et al. Sperm methylation profiles reveal features of epigenetic inheritance and evolution in primates. Cell. 2011;146(6):1029–41. PubMed Central PMCID: PMC3205962. pmid:21925323
  51. 51. Samans B, Yang Y, Krebs S, Sarode GV, Blum H, Reichenbach M, et al. Uniformity of nucleosome preservation pattern in Mammalian sperm and its connection to repetitive DNA elements. Developmental cell. 2014;30(1):23–35. pmid:24998597
  52. 52. Hammoud SS, Nix DA, Zhang H, Purwar J, Carrell DT, Cairns BR. Distinctive chromatin in human sperm packages genes for embryo development. Nature. 2009;460(7254):473–8. PubMed Central PMCID: PMC2858064. pmid:19525931
  53. 53. Brykczynska U, Hisano M, Erkek S, Ramos L, Oakeley EJ, Roloff TC, et al. Repressive and active histone methylation mark distinct promoters in human and mouse spermatozoa. Nature structural & molecular biology. 2010;17(6):679–87.
  54. 54. Davis T, Vaisvila R. High sensitivity 5-hydroxymethylcytosine detection in Balb/C brain tissue. Journal of visualized experiments: JoVE. 2011;(48). Epub 2011/02/11. PubMed Central PMCID: PMC3339836.
  55. 55. Walsh CP, Chaillet JR, Bestor TH. Transcription of IAP endogenous retroviruses is constrained by cytosine methylation. Nature genetics. 1998;20(2):116–7. Epub 1998/10/15. pmid:9771701
  56. 56. Dolinoy DC, Huang D, Jirtle RL. Maternal nutrient supplementation counteracts bisphenol A-induced DNA hypomethylation in early development. Proceedings of the National Academy of Sciences of the United States of America. 2007;104(32):13056–61. Epub 2007/08/03. PubMed Central PMCID: PMC1941790. pmid:17670942
  57. 57. Hancks DC, Kazazian HH Jr. Roles for retrotransposon insertions in human disease. Mobile DNA. 2016;7:9. Epub 2016/05/10. PubMed Central PMCID: PMC4859970. pmid:27158268
  58. 58. Soper SF, van der Heijden GW, Hardiman TC, Goodheart M, Martin SL, de Boer P, et al. Mouse maelstrom, a component of nuage, is essential for spermatogenesis and transposon repression in meiosis. Developmental cell. 2008;15(2):285–97. Epub 2008/08/13. PubMed Central PMCID: PMC2546488. pmid:18694567
  59. 59. Ni K, Dansranjavin T, Rogenhofer N, Oeztuerk N, Deuker J, Bergmann M, et al. TET enzymes are successively expressed during human spermatogenesis and their expression level is pivotal for male fertility. Hum Reprod. 2016. Epub 2016/05/04.
  60. 60. Wang XX, Sun BF, Jiao J, Chong ZC, Chen YS, Wang XL, et al. Genome-wide 5-hydroxymethylcytosine modification pattern is a novel epigenetic feature of globozoospermia. Oncotarget. 2015;6(9):6535–43. PubMed Central PMCID: PMC4466632. pmid:25762640
  61. 61. Xie M, Bu P, Li F, Lan S, Wu H, Yuan L, et al. Neonatal bisphenol A exposure induces meiotic arrest and apoptosis of spermatogenic cells. Oncotarget. 2016;7(9):10606–15. Epub 2016/02/11. pmid:26863571
  62. 62. Jenkins TG, Aston KI, Cairns BR, Carrell DT. Paternal aging and associated intraindividual alterations of global sperm 5-methylcytosine and 5-hydroxymethylcytosine levels. Fertility and sterility. 2013;100(4):945–51. Epub 2013/07/03. pmid:23809503
  63. 63. Mellen M, Ayata P, Dewell S, Kriaucionis S, Heintz N. MeCP2 binds to 5hmC enriched within active genes and accessible chromatin in the nervous system. Cell. 2012;151(7):1417–30. Epub 2012/12/25. PubMed Central PMCID: PMC3653293. pmid:23260135
  64. 64. Wang ZD, Duan L, Zhang ZH, Song SH, Bai GY, Zhang N, et al. Methyl-CpG-Binding Protein 2 Improves the Development of Mouse Somatic Cell Nuclear Transfer Embryos. Cellular reprogramming. 2016;18(2):78–86. Epub 2016/03/17. pmid:26982160
  65. 65. Yeo M, Berglund K, Hanna M, Guo JU, Kittur J, Torres MD, et al. Bisphenol A delays the perinatal chloride shift in cortical neurons by epigenetic effects on the Kcc2 promoter. Proceedings of the National Academy of Sciences of the United States of America. 2013;110(11):4315–20. Epub 2013/02/27. PubMed Central PMCID: PMC3600491. pmid:23440186
  66. 66. Berger SL. Transgenerational inheritance of longevity: epigenetic mysteries abound. Cell metabolism. 2012;15(1):6–7. Epub 2012/01/10. pmid:22225870
  67. 67. Haffner MC, Pellakuru LG, Ghosh S, Lotan TL, Nelson WG, De Marzo AM, et al. Tight correlation of 5-hydroxymethylcytosine and Polycomb marks in health and disease. Cell Cycle. 2013;12(12):1835–41. Epub 2013/05/17. PubMed Central PMCID: PMC3735697. pmid:23676216
  68. 68. Zhou Q, Guo Y, Zheng B, Shao B, Jiang M, Wang G, et al. Establishment of a proteome profile and identification of molecular markers for mouse spermatogonial stem cells. Journal of cellular and molecular medicine. 2015;19(3):521–34. Epub 2014/10/30. PubMed Central PMCID: PMC4369810. pmid:25352495
  69. 69. Lambrot R, Jones S, Saint-Phar S, Kimmins S. Specialized distribution of the histone methyltransferase Ezh2 in the nuclear apical region of round spermatids and its interaction with the histone variant H1t2. Journal of andrology. 2012;33(5):1058–66. Epub 2012/02/11. pmid:22323620
  70. 70. Bhan A, Hussain I, Ansari KI, Bobzean SA, Perrotti LI, Mandal SS. Histone methyltransferase EZH2 is transcriptionally induced by estradiol as well as estrogenic endocrine disruptors bisphenol-A and diethylstilbestrol. Journal of molecular biology. 2014;426(20):3426–41. Epub 2014/08/05. pmid:25088689
  71. 71. Doherty LF, Bromer JG, Zhou Y, Aldad TS, Taylor HS. In utero exposure to diethylstilbestrol (DES) or bisphenol-A (BPA) increases EZH2 expression in the mammary gland: an epigenetic mechanism linking endocrine disruptors to breast cancer. Hormones & cancer. 2010;1(3):146–55. Epub 2011/07/16. PubMed Central PMCID: PMC3140020.
  72. 72. Salian S, Doshi T, Vanage G. Perinatal exposure of rats to Bisphenol A affects the fertility of male offspring. Life sciences. 2009;85(21–22):742–52. Epub 2009/10/20. pmid:19837096
  73. 73. Kurotaki YK, Hatanaka Y, Kamimura S, Oikawa M, Inoue H, Ogonuki N, et al. Impaired active DNA demethylation in zygotes generated by round spermatid injection. Hum Reprod. 2015;30(5):1178–87. Epub 2015/03/06. pmid:25740879
  74. 74. Jia G, Fu X, Cheng K, Yue M, Jia B, Hou Y, et al. Spermatozoa cryopreservation alters pronuclear formation and zygotic DNA demethylation in mice. Theriogenology. 2015;83(6):1000–6. Epub 2014/12/31. pmid:25547286
  75. 75. Gu TP, Guo F, Yang H, Wu HP, Xu GF, Liu W, et al. The role of Tet3 DNA dioxygenase in epigenetic reprogramming by oocytes. Nature. 2011;477(7366):606–10. Epub 2011/09/06. pmid:21892189
  76. 76. Mohn F, Weber M, Schubeler D, Roloff TC. Methylated DNA immunoprecipitation (MeDIP). Methods Mol Biol. 2009;507:55–64. Epub 2008/11/07. pmid:18987806
  77. 77. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9. PubMed Central PMCID: PMC2723002. pmid:19505943
  78. 78. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60. PubMed Central PMCID: PMC2705234. pmid:19451168
  79. 79. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature protocols. 2009;4(1):44–57. pmid:19131956