- Open Access
Ovarian dysfunction following prenatal exposure to an insecticide, chlordecone, associates with altered epigenetic features
- Louis Legoff†1,
- Ouzna Dali†1, 2,
- Shereen Cynthia D’Cruz†1,
- Antonio Suglia1,
- Aurore Gely-Pernot1,
- Chloé Hémery1,
- Pierre-Yves Kernanec1,
- Abbassia Demmouche2,
- Christine Kervarrec1,
- Sergei Tevosian3,
- Luc Multigner1 and
- Fatima Smagulova1Email authorView ORCID ID profile
© The Author(s) 2019
- Received: 21 December 2018
- Accepted: 2 May 2019
- Published: 13 May 2019
Chlordecone (CD) is an insecticide that was used in the French West Indies for several years to control the banana root borer pest. Given its nonsignificant degradation, it persists in the environment. CD is a carcinogenic compound with reproductive and developmental toxicity and is a recognized endocrine-disrupting chemical. In this study, we examined the effects of CD on female reproductive system of mice with the focus on epigenetic features in ovary. Our data show that gestational exposure to low dose of CD affects meiotic double-strand breaks repair in female embryos. In adult mice derived from CD-treated pregnant females, we observed delayed puberty, decreased number of primordial and increased number of atretic follicles. Gene expression analysis revealed that Rcbtb2 and Rbpms genes were not expressed in embryonic gonads. Estrogen signaling- and oocyte maturation-associated genes were downregulated in adult ovaries. The morphological changes were associated with altered epigenetic features: increased H2Aub and increased H3K27me3 and decreased H4ac and H3K4me3 in embryonic oocytes. The DNA damage-associated, γH2AX marks were detected in the follicles of treated but not control adult ovaries. We also found reduced H3K4me3 and H4ac in fully grown oocytes of the treated ovaries. The ChIP-seq analysis of H3K4me3 in adult ovaries showed that target genes of ZFP57 and TRIM28, which regulate pluripotency and imprinting, were significantly enriched in altered regions. Our study clearly demonstrates that gestational exposure to a low dose of CD impairs the function of female reproductive system and the changes are associated with altered epigenetic features.
Environmental factors are known to affect numerous biological processes, particularly development, due to high proliferative activity of cells during organogenesis. While the ability of environmental compounds to cause changes in the DNA makeup of the developing cell is well known, the ability of the chemicals to alter epigenetic processes in the cell is just beginning to be understood. During mammalian development, two rounds of epigenetic reprogramming occur to promote distinct cellular fates that would give rise to germ cells. The first event occurs in preimplantation embryo and the second one in the primordial germ cells when the cells undergo the somatic-to-germline transition (SGT) (e.g., reviewed in ). During SGT, germ cells of both sexes reset their epigenetic marks and become committed to oocyte or male gonocyte fates. Consequently, the developmental window of SGT is very sensitive to exogenous environmental factors.
In male animals, environmentally induced changes have been shown to induce stable modifications in the mammalian germline and contribute to transgenerational effects [2–5]. In some cases, toxicant-induced effects were detected up to the fourth generation of males in the absence of the initial stimulus . In contrast to male gonocytes, the embryonic ovarian germ cells after SGT undergo meiosis, a specialized division process that gives rise to haploid gametes. Subsequently, meiosis gets arrested in females just before birth. During meiosis, DNA is particularly vulnerable to environmental exposure because meiosis entails the formation of double-strand DNA breaks (DSBs) throughout the entire genome. DSBs formation is an integral part of the homologous recombination process, with the major goal of bringing together the parental chromosomes before first meiotic division. DSBs formation normally occurs at discrete regions of the genome known as recombination hotspots; up to 25000 hotspots were found in the mouse genome in meiotic cells . Hence, exposure to toxicants at this stage is likely to result in chromosome missegregation and the formation of abnormal gametes.
Meiosis in females remains dormant until puberty and is completed only after fertilization. Successful meiosis requires the orchestrated actions of a comprehensive histone modification ensemble, which include histone phosphorylation, trimethylation, acetylation, ubiquitination, among others [8, 9]. Histone H2A and H2B ubiquitination are critical for meiosis. For instance, H2A ubiquitination (H2Aub) in meiotic prophase cells was observed at the heterochromatic X and Y chromosomes in males , as well as at unsynapsed chromosomes of both sexes . Ubiquitination of H2A is also important for DNA repair; these marks are abundant at DNA damage sites . It is suggested that E3 ubiquitin ligase RNF8 catalyzes the ubiquitination of H2A and H2AX at damaged sites [13, 14]. On the other hand, H2A ubiquitination is also critical for polycomb repressor complex 1 (PRC1)-dependent repression of developmental genes, which is required to maintain embryonic stem (ES) cell identity . During this process, polycomb repressor complex 2 (PRC2) is first recruited to the developmental genes where it catalyzes the trimethylation of histone H3 at lysine 27 (H3K27me3) [16, 17]. Genes with H3K27me3 marks are subsequently targeted by PRC1, and their chromatin acquires ubiquitination at histone H2A. This PRC1/PRC2 combination regulates, for example, Hox gene family and targets of Oct4/Nanog/Sox2 genes. Monoubiquitination of H2A at developmental genes is catalyzed by E3 ligases, such as RING1A and RING1B . Thus, histone H2Aub appears to play an essential role in gene regulation and maintain DNA integrity during meiosis.
Recent studies have shown that growing oocytes in arrested meiosis undergo another specialized reprogramming event . Specifically, the fully grown oocytes acquire the non-canonical broad domains of histones H3K4me3 [19, 20]. Notably, DNA methylation in growing oocytes negatively correlates with these broad peaks, which suggests that blocking DNA methylation at these domains is essential for the early histone-mediated developmental program after fertilization . Both the establishment of the broad H3K4me3 domains in the oocyte and their timely removal past fertilization are critical for zygote genome activation (ZGA) and early embryogenesis [19, 20]. Some of the broad peaks, however, are preserved in late embryos which points to their potential role in later development .
During early mammalian development, reversal of histone modifications plays an equally important role as their addition. Many cellular factors including nuclear receptors are implicated in the regulation of histone demethylation and thereby could influence the preservation/removal of histone marks. For example, estrogen receptors (ERs) regulate the activity of histone demethylase [21, 22]. ERs, particularly estrogen receptor beta (ESR2), are expressed in many developing tissues, indicating multiple important functions of estrogen signaling during development (reviewed in ). Experimental evidence shows that exposure to environmental toxicants such as bisphenol A (BPA) affects developmental processes, at least partially, through binding to estrogen receptors . BPA also modulates meiosis by affecting chromosomal synapsing, DSBs repair or changing the expression of meiotic genes that rely on the ER function. The organochlorine insecticide, chlordecone (CD, also known as Kepone), has also been shown to bind both ESR1 and ESR2 [25, 26]. Our previous studies on the effects of CD on the male reproductive system showed that CD causes transgenerational effects . We showed that transgenerational effects are associated with alterations in histone H3K4me3 methylation marks residing in DNA locations that are strongly enriched for ESR1 binding motifs in directly exposed F1 and non-exposed F3 generation males . Our study demonstrated that histone epigenetic marks could mediate the transmission of CD-promoted effects to subsequent generations via male germline. Since CD is a very stable compound in the environment, it was predicted that it would persist in the French West Indies for many centuries due to its intensive use from 1973 to 1993 to control the banana root borer pest in the region; therefore, the population there is constantly exposed to this toxicant . Hence, it is very important to reveal the mechanistic effects of CD on organ development including the reproductive system to identify new biomarkers of exposure.
In this new study, we examined the effects of CD on murine female reproductive system exposed in utero, with the focus on epigenetic histones marks that are important for oocyte development and maturation. Our results show that gestational exposure to low dose of CD leads to oxidative stress and affects meiosis and gene expression in embryonic ovaries. In the adult animals from the same lineage, we observed delayed puberty and decreased number of primordial (but increased the number of atretic) follicles. The morphological changes promoted by CD are associated with epigenetic changes in both embryonic and adult ovaries. In summary, CD causes deleterious effects on the female reproductive system, and these effects are associated with modified epigenetic features.
The increase in 8-oxo-G in embryonic oocytes after CD exposure
Meiotic defects are increased in CD-exposed ovaries
To determine whether ovarian DNA repair is associated with altered activity of DNA repair proteins, we analyzed the RNA expression of the genes that encode major DNA repair proteins (DMC1, RAD51, MRE11A). The mRNA expressions of Rad51, Mre11a and Dmc1 have increased, but the changes were not significant (Fig. 2e). In summary, DSB repair efficiency was altered in CD-exposed oocytes.
Increase in H2Aub and H3K27me3 and a decrease in H3K4me3 and H4 acetylation after CD exposure
H4 acetylation is important both for DNA repair and for chromatin-remodeling events . We also analyzed histone acetylation using immunostaining of meiotic spreads against H4ac (Penta, acetylated at residues K5, K8, K12, K16) (Fig. 5c). The quantitative analysis showed that H4ac intensity has decreased 2.1 times in the CD-exposed oocytes (Fig. 5d), but no significant changes were observed in somatic cells. Our data suggest that epigenetic marks associated with transcriptional activities and chromatin remodeling were also altered upon CD exposure in embryonic oocytes.
The effects of CD on embryonic gene expression
In our previous work that examined the effects of CD on males, we found that a large number of genes were deregulated . To identify whether these genes are affected in female mice, we analyzed the expression of some of these genes in E15.5 embryonic ovaries. We found that the expression of Rcbtb2 (regulator of chromosome condensation ) and Rbpms (RNA-binding protein critical repressor of AP-1 signaling ) was totally absent in E15.5-treated ovaries (Additional file 1: Fig. S3), which implies that CD targets similar genes in ovaries as well.
We also analyzed the expression of genes encoding for transcription factors important for ovarian development, FOXL2 and FOXO3. We determined that expression of Foxo3 RNA increased 2.3-fold. In contrast, no significant changes in the expression of Foxl2 were noted. Since CD has estrogenic effects, we tested whether exposure to CD affects the expression level of nuclear receptors such as Esr1 and Esr2, but no significant alterations could be observed (Additional file 1: Fig. S3). Overall, our data demonstrate that gestational exposure to CD affects genes associated with important functions including development and chromatin remodeling.
Exposure to CD leads to delayed puberty and affects the expression of multiple genes important for follicle development
To explore the molecular mechanisms involved in delayed puberty, we performed RT-qPCR gene expression analysis in the ovaries dissected at first estrus. For the analysis, we selected several candidate genes, including germ cell markers (Ddx4, Kit), genes controlling the estrogen signaling pathway (Star, Esr1, Esr2, Inha, Inhba) and the genes involved in transcriptional regulation of ovarian development (Foxl2, Foxo3, Smad4). Because many toxic substances cause oxidative stress, we verified the expression of the genes that are involved in oxidative stress response (Gpx1, Sod1 and Nudt1). We also analyzed the expression of genes that encode for DNA repair proteins (Mre11a, Rad51, Rcbtb2) and the enzymes that introduce histone H3 methylation (Wdr5). We found a decrease in the expression of most of these genes (Fig. 6b). For example, Esr2, Inha, Foxl2 and Foxo3 were downregulated in the treated ovaries (Fig. 6b). Reanalysis of gene expression in the ovaries of three-month-old animals (Additional file 1: Fig. S4) showed a dramatic decrease in the expression of Inhba (10 times), which was associated with a strong decrease in Esr2 (3 times) compared to the control.
Our data suggest that gestational exposure to CD leads to delayed puberty, which is associated with modified expression of genes implicated in the estrogenic signaling pathway, oxidative stress, DNA repair and histone modification.
CD exposure affects the ovarian cell population
In summary, our data demonstrate that exposure to CD leads to a decrease in primordial ovarian follicle population and an increase in atretic follicles, suggesting that CD affects the constitution of ovarian follicle pool and folliculogenesis.
Gestational exposure to CD reduces H3K4me3 and H4ac marks in adult oocytes
H4 acetylation is important for DNA repair . To determine whether CD exposure had any impact on H4 acetylation, we immunostained ovarian sections using an antibody against hyperacetylated histone H4 (Fig. 8c), and the signal was detectable only in fully grown oocytes. We performed analysis of this mark only in fully growing oocytes and in surrounding granulosa cells. The quantitative analysis of immunofluorescence showed 2.6 times decrease in H4 acetylation in treated oocytes compared to control with no change in surrounding granulosa cells (Fig. 8d).
To sum up, our analysis reveals that gestational exposure to CD associates with the decreased level of histone H3 trimethylation and H4 acetylation in fully grown oocytes.
The genome-wide analysis of H3K4me3 revealed changes in genes important for regulating imprinting and pluripotency
We also assigned regions to the genes located within differential H3K4me3 peaks using GREAT and identified 185 genes that are located near the altered H3K4me3 peaks. Our analysis revealed that genes associated with mitotic spindle orientation (Dynlt1a, Dynlt1b) were most enriched in differential peaks (Fig. 10d).
To explore the role for transcriptional regulation in altering epigenetic modifications in the CD-exposed ovaries, we analyzed the altered H3K4me3 peaks for motif enrichment using MEME-ChIP as described in Methods section. We identified two significantly enriched motifs in altered H3K4me3 peaks; the part of the first motif is significantly similar to ZBTB17 and WT1 binding sites (p value = 5.94e−06 and 1.50e−05 and q value = 5.42e−03 and 7.62e−03, respectively) (Fig. 10e). For example, differential peak in Cop1 gene has ZBTB17 binding site (Additional file 1: Fig. S7). In the second identified motif, part of the sequence is similar to ONECUT2 binding site (p value = 2.42e−05, q value = 1.54e−02) (Fig. 10f).
In summary, our data suggest that changes in limited regions in the mouse genome, such as genes associated with mitotic division function, could contribute to observed phenotype, such as abnormal mitosis of granulosa cells, which could lead to atresia of ovarian follicles.
The embryonic effects of CD exposure
In this study, we examined the effects of CD on ovarian function. Given the potential of CD to cause transgenerational effects and the fact that humans remain exposed to CD despite discontinued use, understanding the mechanisms of CD actions on female germline is of great importance. In this study, we used a low dose of CD (100 ug/kg/day), which was previously considered to be safe . This dose was used in our previous study, and we showed that it leads to transgenerational effects in males . So, we asked the question whether this dose could promote changes in female reproductive system. Following exposure to CD, we analyzed the ovaries of mice at several developmental time points: at prophase I of meiosis (embryonic day 15 and 17), at first estrus (day 27 ± 1, in control) and at the age of three and five months. We found that CD causes oxidative stress in ovaries, which is manifested as increased production of 8-oxo-G in the nuclei of oocytes. It has been established that exposure to chemicals leads to the formation of excess reactive oxygen species (ROS) via several mechanisms, including the mitochondrial respiration and the detoxification pathways (reviewed in ). Guanine nucleotide is particularly vulnerable to oxidation and could lead to the formation of 8-oxo-derivative . The increase in the expression of Nudt1 mRNA confirms that cells have responded to increased production of 8-oxo-G so that they could be eliminated. We believe that formation of 8-oxo-G affected meiosis because it could delay meiotic DNA repair. The higher number of DMC1 foci at 15.5 that persists until E17.5 in CD-treated oocytes could be explained by the complexity of the regulation of meiotic DSB repair. It has been found that RPA regulates the loading of DMC1 on DNA . RPA-mutant mice do not form DMC1 foci in spite of the presence of DMC1 protein in the cells , suggesting that recombinase activity of DMC1 is regulated by protein factors. Up to date, many accessory proteins were found to cooperate with DMC1 to facilitate the homologous search and strand exchange activity, and among them are HOP2, MND1, SWI5, SFR1 and MEIOB (reviewed in ).
On the other hand, increased DSBs that we observed in the treated samples could also be the result of inefficient DSBs repair due to decreased levels of cellular ATP that arise from elevated oxidative stress. Most DNA repair proteins are ATP-dependent. Moreover, chromosome mobility that is essential for repair of DSBs [50, 51] as well as homologous search  are both ATP-dependent processes. However, another reason for the DSBs increase could be related to estrogenic properties of CD; it has been shown that estrogens are capable of inducing DSBs at gene promoters . The dramatic increase in H2Aub in zygotene stage (where DSBs started to repair) of CD-exposed ovaries suggests that some breaks are induced via DNA damage and they are persistent. We found that changes in histone H2Aub correlate with an increase in H3K27me3. This could be due to the functional link between enzymes, which introduce these two marks [28, 29]. The increase in mRNA expression of histone methyltransferase EZH2, a catalytic unit of the polycomb repressive complex 2 (PRC2) that introduces H3K27me3 modification, is in consistence with the global increase in H3K27me3 in meiotic and somatic cells.
We also found a dramatic decrease in the expression of genes such as Rcbtb2 and Rbpms in embryonic ovaries. Rcbtb2 is essential for the formation of primary cilia in mammalian cells . The presence of primary cilia has been documented in granulosa cells, antral follicles and ovarian surface epithelium [54, 55]. The major function of primary cilia is to sense a variety of external stimuli; several developmentally important signaling pathways are active in primary cilia. The absence or defects in the formation of primary cilia has been observed in epithelial ovarian cancer  that comprises the vast majority of the human ovarian cancers . Therefore, the absence of Rcbtb2 mRNA expression in the embryonic ovaries of CD-treated mice could have deleterious effects on ovarian signaling pathways. Among the downregulated genes, Rbpms has a broad expression pattern in adult ovaries  and is an important regulator of oocyte development . Studies using Zebrafish have shown that Rbpms mutants fail to undergo oogenesis during sexual differentiation . In humans, moderate-to-negligible protein expression of RBPMS has been observed in serous papillary ovarian cancer tissues when compared to control ovaries . Moreover, proliferation and invasion of ovarian cancer cells could be impeded by inhibiting microRNA-21-associated target genes that included RBPMS, RCBT1 and ZNF608 . Therefore, CD-induced decreased expression of the genes essential for ovarian functions could have detrimental effects on ovarian development.
The effects of CD on follicle population
We found several changes in the ovaries of adult mice that had undergone in utero CD exposure. The delay in puberty and perturbed gene expression suggest that CD exposure affects signaling pathways in the ovary. Specifically, we observed a dramatic reduction in the expression of Esr2, Inhba, Inha and Smad4 in the ovaries of three-month-old animals that reflect an abnormal granulosa cell state. The Inhba, Inha and Smad4 genes implicated in signaling (TGF-beta) pathway play an essential role for granulosa cell growth and differentiation as well as ovarian steroidogenesis . Inha knockout mice exhibited numerous multilayered follicles, and they are in an advanced stage compared to controls . Activin B is a regulator of estrogens in the mouse ovary . Decreased Esr2 expression observed in our study could be due to reduced activin B levels. On the other hand, changes in Esr2 are associated with changes in Foxl2 expression. It has been shown that FOXL2 is required for efficient gene regulation by steroid receptors in follicular cells, with ESR2 being the main vector of estradiol signaling . The dramatic decrease in Esr2 in three-month-old ovaries and delayed puberty suggests that CD causes changes in the female reproductive system, and it seems to occur via estrogenic signaling pathways. In summary, profound changes in gene expression suggest underlying defects in granulosa cell development and function, which could explain ovarian pathology (e.g., the presence of a large number of atretic follicles and the decrease in primordial follicles).
Epigenetic effects in adult ovary
We found a decrease in H3K4me3 staining in fully grown oocytes. As discussed in Introduction, growing oocytes gain broad domains of H3K4me3 during maturation, and these domains are essential for the normal development of the organism. These regions are also candidates for transgenerational inheritance because some of the broad domains are preserved in embryos suggesting their importance in development. Therefore, changes in the epigenetic state could perturb the future organism. However, the exact mechanism through which this occurs is yet to be identified. We found that expression of Wdr5 has decreased. WDR5 is a part of the MLL2 histone methyltransferase complex and it gets recruited to the promoters of ESR1 target genes upon estrogen stimulation . The exact role of WDR5 in the ovary is not investigated. However, an abundance of WDR5 transcript was observed in porcine oocytes at the germinal vesicle stage, followed by a stage-specific expression pattern in early embryos . Knockdown of WDR5 using siRNA had an impact on the expression of several pluripotent genes along with DNA repair defects in blastocysts, indicating that WDR5 is possibly a maternal-effect gene that is crucial for early development [65–67]. Therefore, it is conceivable that decreased expression of Wdr5 caused by a reduction in ERs contributed to the global decrease in H3K4me3 levels and could possibly have an effect on development.
We found that altered H3K4me3 marks are enriched in gene targets of ZFP57 and TRIM28 (also known as KAP1). ZFP57 and TRIM28 are both required for maintaining the repressive DNA methylation and H3-lysine-9-trimethylation (H3K9me3) at imprint control regions (ICRs) [68, 69]. TRIM28 is involved in heterochromatin formation during the process of reprogramming of mouse somatic cells into induced pluripotent stem cells, suggesting its role in reprogramming . In our previous study on transgenerational inheritance in males exposed to CD, we found that altered H3K4me3 regions were also enriched in ZFP57 and TRIM28 targets and the changes were detected up to third generations . It is likely that the changes were triggered during somatic-to-germline transition, and it seems to be the mechanism of CD action on female lineage as well. Importantly, haploinsufficiency of Samd9 l gene, another target of ZFP57, causes myeloid malignancies in mice mimicking human diseases with monosomy 7 . Chd9 is also a target of ZFP57, and it is highly expressed in oocytes where it plays an important role in the acquisition of a highly loosened chromatin structure . We found epigenetic changes in genes important for mitotic function (e.g., Dynlt1f), which could be implicated in abnormal granulosa cell division. It is conceivable that decreased histone H3K4me3 intensity in oocyte nuclei could be affected by chromatin structure. Further work is required to confirm this.
How embryonic effects of CD on ovary could affect oocyte maturation?
The decrease in primordial follicle pool suggests that some of these follicles were eliminated during the late prophase I. It was proposed that checkpoint activity (CHK2) gets triggered in oocytes that carry more than 10 DSBs in late prophase I, and their repair is inhibited by HORMAD1/2 by preventing intersister recombination . In addition, low levels of H3K4me3 in embryonic oocytes and in adult growing oocytes suggest that the epigenetic state changes might have got altered in utero which sustained up to adulthood. We believe that a decrease in primordial follicles could be explained by prophase I meiotic defects; these oocytes could be eliminated and leads to lower primordial follicle population. The increase in atretic follicles could be explained by the fact that gestational exposure during SGT transition introduces epigenetic changes in important genes, which are implicated in cell division. Indeed, our ChIP-seq data showed that only limited regions have altered H3K4me3 marks, but these regions contain genes important for mitotic spindle orientation, chromatin structure, regulation of pluripotency and reprogramming. We believe that these changes occurred during early development and preserved at limited regions in the adults. We identified binding sites for important developmental genes such Wt1, Zbtb17 and Onecut2. It is conceivable that exposure to CD affects these transcriptional factors and these factors, in turn, could change the transcription activities of their target genes.
In summary, exposure to a low dose of CD leads to meiosis defects and decreases the oocytes stockpile. Our results demonstrate that histone epigenetic marks could be one of the responsible factors for the effects of CD on adult ovaries.
We used the same dose of CD as reported in our previous study . The dose that we used is relatively a low dose (100 µkg/kg/day) and does not cause general toxic effects in mice. The exposure window was embryonic day E6.5 to E15.5, which corresponds to somatic-to-germline transition. Briefly, CD was suspended in olive oil and administered via oral gavage at a dose of 100 μg/kg/day in a volume of 150 μl for each mouse. The control mice were treated with the same volume of oil. Ovaries were analyzed at embryonic day E15 and E17.5, and at adult stage at first estrus (day 27 ± 1, in control) and in three- and five-month-old animals.
Preparation and immunostaining of structurally preserved nuclei
Structurally preserved nuclei (SPN) for three-dimensional analysis were prepared by cutting fresh E15.5 ovarian tissues in DMEM medium (Life Technologies, GIBCO) with 0.5% protease inhibitor. The ovary suspensions were mixed with equal amounts of 3.7% (vol/vol) paraformaldehyde and 0.1 M sucrose and spread on glass slides; the slides were air-dried and kept at −80 °C. Slides were washed several times with PBS before use and 2 min with 0.1 M glycine in PBS to remove the traces of paraformaldehyde. Slides were permeabilized during 30 min in PBS/0.5% Triton at RT, washed with PBS and blocked for 30 min in solution containing 0.1% (v/v) donkey serum, 0.03% (w/v) BSA, and 0.005% (v/v) Triton X-100 in PBS. SPN were immunostained with mouse monoclonal 8-oxo-G (1:100, Merck Millipore, MAB3560) and rabbit polyclonal SYCP3 (1:100, Santa Cruz, SC-33195) antibodies at 37C for 1 h followed by several washes and incubation with fluorescent Alexa antibodies. Z-stacks were acquired with 500 nm step; 21 individual planes were taken for each individual channel, for DAPI (Blue, 350 nm), SYCP3 (Green, 488 nm) and 8-oxo-G (Red, 594 nm) using Zen Pro (version 2.3) program. All images for control and exposed samples were taken with fixed exposure time. Deconvolution was performed using the “Fast Iterative” algorithm provided by Zen Pro. The mean intensity images were generated for each z-stack, and the resulting images were analyzed in ImageJ v1.52n. We used the lasso tool for nucleus and cytoplasm contouring, and the integrated density immunofluorescence for each nucleus was calculated. The similar area with background was used, and the background was subtracted. We analyzed 4 independent biological replicates for control and treated groups, and at least 12 cells for each replicate at zygotene stage were analyzed. The data were plotted in excel and presented as corrected total cell fluorescence (CTCF) of normalized fluorescence for nucleus or cytoplasm compared to control ± SEM, *p < 0.05, **p < 0.01, t test.
Meiotic chromosome spreads
The ovaries were dissected at E15.5 and E17.5 and kept in ice-cold PBS. The ovaries were placed in a drop of 100 mM sucrose (20 µl), disrupted carefully with tweezers and pipetted until a cell suspension was formed. The cells were transferred onto a glass slide with 100 µl of 1% paraformaldehyde containing 0.1% (v/v) Triton X-100 and were kept for 2–4 h in a humidified chamber at RT and then dried. The slides were washed 4 times for 1 min each in 0.4% (v/v) Kodak Photo-Flo solution, air-dried and kept at − 80 °C until use.
Immunolabeling of ovarian spreads
The ovarian spreads were blocked with blocking buffer (0.1% (v/v) donkey serum, 0.03% (w/v) BSA, and 0.005% (v/v) Triton X-100 in PBS) for 20 min at 37 °C in a humidified chamber. The slides were incubated with a primary antibody diluted in blocking solution for 1 h at 37 °C in a humidified chamber, washed three times with 0.4% (v/v) Kodak Photo-Flo/PBS and incubated with a secondary fluorescent Alexa antibody diluted 1:500 for 20 min at RT. The slides were washed three times with 0.4% (v/v) Kodak Photo-Flo/PBS and mounted with Vectashield mounting media containing 5 µg/ml DAPI (Vector Laboratories, Burlingame, CA). The slides were analyzed with an epifluorescence Axio Imager M1 microscope (ZEISS, France), and the pictures were taken using a 63X or 100X objective with AxioVision Rel 4.7.1 (Zeiss, France). For DMC1 foci visualization, the slides containing E15.5 and E17.5 ovaries were incubated with primary rabbit anti-DMC1 (1:100, Santa Cruz, SC-22768) and mouse anti-SYPC3 (1:100, Santa Cruz, SC-33195) antibodies. DMC1 foci were counted manually using Photoshop, and the average number of foci per cell were compared between control and treated groups. For H2Aub visualization, the slides containing E15.5 ovaries were incubated with primary mouse anti-H2Aub (1:100, Merck Millipore, 05-678) and rabbit anti-SYPC3 (1:100, Santa Cruz, SC-33195) antibodies. For H3K27me3 visualization, the slides containing E15.5 ovaries were incubated with primary mouse anti-H3K27me3 (1:100, Abcam, ab6002) and rabbit anti-SYPC3 (1:100, Santa Cruz, SC-33195) antibodies. For H3K4me3 visualization, the slides containing E15.5 ovaries were incubated with primary rabbit anti-H3K4me3 (1:100, Millipore, 07-473) and rabbit anti-SYPC3 (1:100, Santa Cruz, SC-33195) antibodies, and for H4Ac visualization, the slides containing E15.5 ovaries were incubated with primary rabbit anti-H4penta (1:100, Millipore, 06-946) and mouse anti-SYPC3 (1:100, Santa Cruz, SC-74569) antibodies. We used minimum 4 biological replicates for control and treated groups derived from 4 different treated litters. We analyzed at least 20 cells for each biological replicate. We used ImageJ for quantitative analysis of immunofluorescence of histone marks. We measured the intensity of fluorescence by contouring the nuclei with lasso tool of ImageJ. The signals for histone marks, SYCP3, DAPI and background were measured. We calculated CTCF (corrected total cell fluorescence) as multiplication of mean value of fluorescence to cell area with background subtraction. The immunofluorescence signal was normalized to SYCP3 or DAPI signal and data presented as normalized fluorescence compared to control.
Immunohistochemistry of sections
For immunostaining, the ovaries from control and CD-treated groups were fixed in 4% (w/v) PFA solution for 16 h, dehydrated and embedded in paraffin. The sections were deparaffinized and rehydrated, and the epitopes were unmasked in 0.01 M citrate buffer, pH 6 at 80 °C for 45 min. After washing in 1X PBS-0.05% Tween (PBS-T), the sections were incubated with rabbit anti-H3K4me3 (1:500, Merck Millipore, 07-473) and mouse anti-MSY2 (1:500, Santa Cruz, SC-3938440) antibody or rabbit polyclonal hyperacetylated histone 4 (1:500, Merck Millipore, 06-956) and mouse anti-MSY2, rabbit polyclonal anti-phosphorylated Histone H2AX (1:500, γH2AX, Trevigen, 4411-PC-100) or goat anti-AMH antibody (1:500, Santa Cruz, SC-6886). The sections with primary antibodies were incubated in PBS-T overnight at 4 °C in a humidified chamber. After washing in PBS-T, the sections were incubated with appropriate fluorescent secondary antibody (1:100) for 1 h in a humidified chamber at room temperature. The sections were counterstained with 0.001% (v/v) 4,6-diamidino-2-phenylindole dihydrochloride (DAPI) and mounted using Vectashield solution. The images were taken using an AxioImager microscope equipped with an AxioCam MRc5 camera and AxioVision software version 4.8.2 (Zeiss, Le Pecq, France) with a 5X or 40X objective lens. We analyzed intensity of IF only in antral follicles. For the analysis, we chose oocytes with visibly similar size nucleus. We analyzed minimum 4 follicles per replicate and quantitated signal intensity in oocyte and surrounding granulosa cells separately using ImageJ v1.52n. We subtracted background using the regions with no cells. The mean intensity for H3K4me3 or H4Ac signals was counted, normalized to oocyte marker MSY2 and presented as normalized mean fluorescence compared to control.
For γH2AX analysis, we analyzed sections with 5X magnification. We considered the follicle as impaired if the number of γH2AX-positive cells exceeded 5. We counted on average at least 4 sections for each biological replicate and presented the data as a number of “damaged” follicles per section.
Quantification of germ cells and follicles
For the analysis of follicles, ovaries were fixed in Bouin solution, washed in PBS and 70% ethanol, dehydrated and embedded in paraffin. Five-micrometer sections of the entire ovaries were cut and every 5th section was placed on the slide. The sections were deparaffinized, rehydrated and stained with PAS hematoxylin. The follicles were quantified using the NDPview software. The follicles were counted when the nuclei of the oocytes were visible in the follicle. The follicles were classified as primordial when the oocyte was surrounded by a single layer of flattened granulosa cells; as primary, when at least one of the granulosa cells of the single layer became cuboidal; as secondary/preantral, when the number of granulosa cell layers exceeded one; and as antral, when the antrum cavity appeared. The follicles were classified as atretic when they contained at least 2 pyknotic granulosa cells. We analyzed on average 91 sections per control and 80 sections per CD-treated samples.
RNA extraction and quantitative PCR
Total RNA was extracted using the RNeasy plus mini kit (Qiagen) according to the manufacturer’s instructions. For embryonic ovary analysis, RNA was extracted from pool of 4 embryonic ovaries. We used 5 biological replicates for control and 6 for treatment group analysis. For RNA analysis in first estrus mice, 8 individual ovaries from 4 different litters from each group were used. For RNA analysis in three-month-old mice, 6 individual ovaries from 6 litters of each group were used. Reverse transcription was performed with 1 µg of RNA using the iScript™ cDNA Synthesis Kit (Bio-Rad). The resulting cDNA was diluted 10 times and used for quantitative PCR. The primer sequences used for qPCR are shown in Additional file 1: Table S4. QPCR was performed using the iTaq Universal SYBR Green Supermix (Bio-Rad) according to the manufacturer’s instructions on a CFX384 Touch Real-Time PCR Detection system (Bio-Rad). The quantification cycle (Cq) values were calculated with Bio-Rad CFX Manager 3.1. Cq values of Hprt and Dazl (embryonic gonad) or Rpl37a (adult ovaries) cDNA were used for normalization. The data were analyzed and are presented as mean values of Fold Change (FC) compared to control ± SEM, *p < 0.05, **, P < 0.01, *** p < 0.001, t test.
Chromatin immunoprecipitation and library preparation
We used two ovaries from different animals for each biological replicate; we used two replicates for control and treated groups. Chromatin immunoprecipitation in adult tissue was performed as previously described , except for small modifications. Specifically, chromatin was not dialyzed after sonication but instead diluted 6 times with ChIP-dilution buffer. We performed ChIP using 1 µg of H3K4me3 (Merck Millipore, 07-473). Immunoprecipitated DNA was purified with Qiagen MinElute kit, and 10 µg of DNA was used for library construction. Library was prepared using NEBNext® Ultra™ II DNA Library Prep Kit for Illumina (E7645) and submitted to Genomic platform to perform massively parallel 50-bp sequencing in single-end mode.
Peak calling and identification of differential peaks
All analyses were performed as previously described . Approximately 77.6 million tags were derived from the anti-H3K4me3 ChIP. The resulting sequences were mapped back to the mouse mm10/Ensembl genome using Bowtie 1.0.0 with a seed length of 40. H3K4me3 mark peaks were identified using MACS 2.1.1 with two biological replicate samples, including the corresponding input, using a shift-size window of 73 bp, no model, and a p value threshold < 1e−05. Statistical significance was calculated using a Limma test.
Functional annotation of ChIP-seq data
Functional annotation was performed using GREAT version 3.0.0 or EnrichR with default parameters . We used ChEA to functionally annotate differential peaks. ChEA is a database of genes potentially regulated by specific transcription factors, where the data were extracted manually from different ChIP-seq experiments for each transcription factor. Motif finding was performed using MEME-ChIP  with the default parameters. Identified motifs were compared with known motifs using TomTom  with the default parameters.
All sequencing data from this study are publicly available and have been deposited in the National Center for Biotechnology Information Gene Expression Omnibus (GEO); GSE number is GSE129897.
The authors express their gratitude to Severine Mazaud-Guittot for her help with ovarian section staining and classification of follicles. We are thankful for H2P2 platform of Rennes for assistance with tissue preparation for histology. Sequencing was performed by the GenomEast platform, a member of the “France Génomique” consortium (ANR-10-INBS-0009). We also would like to thank the French Institute of Bioinformatics (IFB, ANR-11-INBS-0013) for providing storage and computing resources on its national life science Cloud.
This work was supported by Atip-Avenir program (R13139NS) to FS. LL was supported by ARED/INSERM fellowship.
FS designed the research. FS, PYK, AGP and CK did animal work. OD and SCD performed the histological section, immunohistochemistry and RT-qPCR analysis. LL participated in the follicle analysis, performed quantitative analysis of z-stacks, prepared libraries for the genome-wide sequencing, analyzed the sequencing data and prepared figures for the manuscript. AS performed quantitative analysis of histone marks images. AD, LM, SCD and ST participated in discussion and commented on the manuscript. LL, ST, SCD and FS wrote the manuscript. All authors read and approved the final manuscript.
The animal facility used for the present study is licensed by the French Ministry of Agriculture (agreement D35–238–19). All animal procedures were performed according to the Ethics Committee of the Ministry of the Research of France (agreement number: 01861.02). All experimental procedures followed the ethical principles outlined in the Ministry of Research Guide for Care and Use of the Laboratory Animals and were approved by the local Animal Experimentation Ethics Committee (C2EA-07).
Consent for publication
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- 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(1609):20110330.PubMedPubMed CentralGoogle Scholar
- Gely-Pernot A, Hao C, Legoff L, Multigner L, D’Cruz SC, Kervarrec C, Jegou B, Tevosian S, Smagulova F. Gestational exposure to chlordecone promotes transgenerational changes in the murine reproductive system of males. Sci Rep. 2018;8(1):10274.PubMedPubMed CentralGoogle Scholar
- Hao C, Gely-Pernot A, Kervarrec C, Boudjema M, Becker E, Khil P, Tevosian S, Jegou B, Smagulova F. Exposure to the widely used herbicide atrazine results in deregulation of global tissue-specific RNA transcription in the third generation and is associated with a global decrease of histone trimethylation in mice. Nucleic Acids Res. 2016;44(20):9784–802.PubMedPubMed CentralGoogle Scholar
- Ben Maamar M, Sadler-Riggleman I, Beck D, McBirney M, Nilsson E, Klukovich R, Xie Y, Tang C, Yan W, Skinner MK. Alterations in sperm DNA methylation, non-coding RNA expression, and histone retention mediate vinclozolin-induced epigenetic transgenerational inheritance of disease. Environ Epigenetics. 2018;4(2):dvy010.Google Scholar
- McBirney M, King SE, Pappalardo M, Houser E, Unkefer M, Nilsson E, Sadler-Riggleman I, Beck D, Winchester P, Skinner MK. Atrazine induced epigenetic transgenerational inheritance of disease, lean phenotype and sperm epimutation pathology biomarkers. PLoS ONE. 2017;12(9):e0184306.PubMedPubMed CentralGoogle Scholar
- Anway MD, Cupp AS, Uzumcu M, Skinner MK. Epigenetic transgenerational actions of endocrine disruptors and male fertility. Science. 2005;308(5727):1466–9.PubMedGoogle Scholar
- Brick K, Smagulova F, Khil P, Camerini-Otero RD, Petukhova GV. Genetic recombination is directed away from functional genomic elements in mice. Nature. 2012;485(7400):642–5.PubMedPubMed CentralGoogle Scholar
- Gu L, Wang Q, Sun QY. Histone modifications during mammalian oocyte maturation: dynamics, regulation and functions. Cell Cycle. 2010;9(10):1942–50.PubMedGoogle Scholar
- Kota SK, Feil R. Epigenetic transitions in germ cell development and meiosis. Dev Cell. 2010;19(5):675–86.Google Scholar
- Baarends WM, Hoogerbrugge JW, Roest HP, Ooms M, Vreeburg J, Hoeijmakers JH, Grootegoed JA. Histone ubiquitination and chromatin remodeling in mouse spermatogenesis. Dev Biol. 1999;207(2):322–33.PubMedGoogle Scholar
- Baarends WM, Wassenaar E, van der Laan R, Hoogerbrugge J, Sleddens-Linkels E, Hoeijmakers JH, de Boer P, Grootegoed JA. Silencing of unpaired chromatin and histone H2A ubiquitination in mammalian meiosis. Mol Cell Biol. 2005;25(3):1041–53.PubMedPubMed CentralGoogle Scholar
- Bergink S, Salomons FA, Hoogstraten D, Groothuis TA, de Waard H, Wu J, Yuan L, Citterio E, Houtsmuller AB, Neefjes J, et al. DNA damage triggers nucleotide excision repair-dependent monoubiquitylation of histone H2A. Genes Dev. 2006;20(10):1343–52.PubMedPubMed CentralGoogle Scholar
- Mailand N, Bekker-Jensen S, Faustrup H, Melander F, Bartek J, Lukas C, Lukas J. RNF8 ubiquitylates histones at DNA double-strand breaks and promotes assembly of repair proteins. Cell. 2007;131(5):887–900.PubMedPubMed CentralGoogle Scholar
- Huen MS, Grant R, Manke I, Minn K, Yu X, Yaffe MB, Chen J. RNF8 transduces the DNA-damage signal via histone ubiquitylation and checkpoint protein assembly. Cell. 2007;131(5):901–14.PubMedPubMed CentralGoogle Scholar
- Endoh M, Endo TA, Endoh T, Isono K, Sharif J, Ohara O, Toyoda T, Ito T, Eskeland R, Bickmore WA, et al. Histone H2A mono-ubiquitination is a crucial step to mediate PRC1-dependent repression of developmental genes to maintain ES cell identity. PLoS Genet. 2012;8(7):e1002774.PubMedPubMed CentralGoogle Scholar
- Czermin B, Melfi R, McCabe D, Seitz V, Imhof A, Pirrotta V. Drosophila enhancer of Zeste/ESC complexes have a histone H3 methyltransferase activity that marks chromosomal Polycomb sites. Cell. 2002;111(2):185–96.Google Scholar
- Cao R, Wang L, Wang H, Xia L, Erdjument-Bromage H, Tempst P, Jones RS, Zhang Y. Role of histone H3 lysine 27 methylation in Polycomb-group silencing. Science. 2002;298(5595):1039–43.PubMedGoogle Scholar
- Wang H, Wang L, Erdjument-Bromage H, Vidal M, Tempst P, Jones RS, Zhang Y. Role of histone H2A ubiquitination in Polycomb silencing. Nature. 2004;431(7010):873–8.PubMedGoogle Scholar
- Dahl JA, Jung I, Aanes H, Greggains GD, Manaf A, Lerdrup M, Li G, Kuan S, Li B, Lee AY, et al. Broad histone H3K4me3 domains in mouse oocytes modulate maternal-to-zygotic transition. Nature. 2016;537(7621):548–52.PubMedPubMed CentralGoogle Scholar
- Zhang B, Zheng H, Huang B, Li W, Xiang Y, Peng X, Ming J, Wu X, Zhang Y, Xu Q, et al. Allelic reprogramming of the histone modification H3K4me3 in early mammalian development. Nature. 2016;537(7621):553–7.PubMedGoogle Scholar
- Garcia-Bassets I, Kwon YS, Telese F, Prefontaine GG, Hutt KR, Cheng CS, Ju BG, Ohgi KA, Wang J, Escoubet-Lozach L, et al. Histone methylation-dependent mechanisms impose ligand dependency for gene activation by nuclear receptors. Cell. 2007;128(3):505–18.PubMedPubMed CentralGoogle Scholar
- Carnesecchi J, Cerutti C, Vanacker JM, Forcet C. ERRalpha protein is stabilized by LSD1 in a demethylation-independent manner. PLoS ONE. 2017;12(11):e0188871.PubMedPubMed CentralGoogle Scholar
- Bondesson M, Hao R, Lin CY, Williams C, Gustafsson JA. Estrogen receptor signaling during vertebrate development. Biochem Biophys Acta. 2015;1849(2):142–51.PubMedGoogle Scholar
- Ozgyin L, Erdos E, Bojcsuk D, Balint BL. Nuclear receptors in transgenerational epigenetic inheritance. Prog Biophys Mol Biol. 2015;118(1–2):34–43.PubMedGoogle Scholar
- Kuiper GG, Lemmen JG, Carlsson B, Corton JC, Safe SH, van der Saag PT, van der Burg B, Gustafsson JA. Interaction of estrogenic chemicals and phytoestrogens with estrogen receptor beta. Endocrinology. 1998;139(10):4252–63.PubMedGoogle Scholar
- Lemaire G, Mnif W, Mauvais P, Balaguer P, Rahmani R. Activation of alpha- and beta-estrogen receptors by persistent pesticides in reporter cell lines. Life Sci. 2006;79(12):1160–9.PubMedGoogle Scholar
- Multigner L, Kadhel P, Rouget F, Blanchet P, Cordier S. Chlordecone exposure and adverse effects in French West Indies populations. Environ Sci Pollut Res Int. 2016;23(1):3–8.PubMedGoogle Scholar
- Cooper S, Grijzenhout A, Underwood E, Ancelin K, Zhang T, Nesterova TB, Anil-Kirmizitas B, Bassett A, Kooistra SM, Agger K, et al. Jarid2 binds mono-ubiquitylated H2A lysine 119 to mediate crosstalk between Polycomb complexes PRC1 and PRC2. Nat Commun. 2016;7:13661.PubMedPubMed CentralGoogle Scholar
- Kalb R, Latwiel S, Baymaz HI, Jansen PW, Muller CW, Vermeulen M, Muller J. Histone H2A monoubiquitination promotes histone H3 methylation in Polycomb repression. Nat Struct Mol Biol. 2014;21(6):569–71.PubMedGoogle Scholar
- Guenther MG, Levine SS, Boyer LA, Jaenisch R, Young RA. A chromatin landmark and transcription initiation at most promoters in human cells. Cell. 2007;130(1):77–88.PubMedPubMed CentralGoogle Scholar
- Xu J, Ma H, Jin J, Uttam S, Fu R, Huang Y, Liu Y. Super-resolution imaging of higher-order chromatin structures at different epigenomic states in single mammalian cells. Cell Rep. 2018;24(4):873–82.PubMedPubMed CentralGoogle Scholar
- Dhar S, Gursoy-Yuzugullu O, Parasuram R, Price BD. The tale of a tail: histone H4 acetylation and the repair of DNA breaks. Philos Trans R Soc Lond Ser B Biol Sci. 2017;372(1731):1–12.Google Scholar
- Zhang L, Li W, Ni J, Wu J, Liu J, Zhang Z, Zhang Y, Li H, Shi Y, Teves ME, et al. RC/BTB2 is essential for formation of primary cilia in mammalian cells. Cytoskeleton. 2015;72(4):171–81.PubMedGoogle Scholar
- Fu J, Cheng L, Wang Y, Yuan P, Xu X, Ding L, Zhang H, Jiang K, Song H, Chen Z, et al. The RNA-binding protein RBPMS1 represses AP-1 signaling and regulates breast cancer cell proliferation and migration. Biochem Biophys Acta. 2015;1853(1):1–13.PubMedGoogle Scholar
- Iwase A, Nakamura T, Osuka S, Takikawa S, Goto M, Kikkawa F. Anti-Mullerian hormone as a marker of ovarian reserve: what have we learned, and what should we know? Reprod Med Biol. 2016;15(3):127–36.PubMedGoogle Scholar
- Kouzarides T. Chromatin modifications and their function. Cell. 2007;128(4):693–705.Google Scholar
- Vermeulen M, Eberl HC, Matarese F, Marks H, Denissov S, Butter F, Lee KK, Olsen JV, Hyman AA, Stunnenberg HG, et al. Quantitative interaction proteomics and genome-wide profiling of epigenetic histone marks and their readers. Cell. 2010;142(6):967–80.PubMedGoogle Scholar
- Sharma GG, So S, Gupta A, Kumar R, Cayrou C, Avvakumov N, Bhadra U, Pandita RK, Porteus MH, Chen DJ, et al. MOF and histone H4 acetylation at lysine 16 are critical for DNA damage response and double-strand break repair. Mol Cell Biol. 2010;30(14):3582–95.PubMedPubMed CentralGoogle Scholar
- Antognelli C, Cecchetti R, Riuzzi F, Peirce MJ, Talesa VN. Glyoxalase 1 sustains the metastatic phenotype of prostate cancer cells via EMT control. J Cell Mol Med. 2018;22(5):2865–83.PubMedPubMed CentralGoogle Scholar
- Tamori S, Nozaki Y, Motomura H, Nakane H, Katayama R, Onaga C, Kikuchi E, Shimada N, Suzuki Y, Noike M, et al. Glyoxalase 1 gene is highly expressed in basal-like human breast cancers and contributes to survival of ALDH1-positive breast cancer stem cells. Oncotarget. 2018;9(92):36515–29.PubMedPubMed CentralGoogle Scholar
- Jandial R, Neman J, Lim PP, Tamae D, Kowolik CM, Wuenschell GE, Shuck SC, Ciminera AK, De Jesus LR, Ouyang C et al. Inhibition of GLO1 in glioblastoma multiforme increases DNA-AGEs, stimulates RAGE expression, and inhibits brain tumor growth in orthotopic mouse models. Int J Mol Sci. 2018;19(2):1–21.Google Scholar
- Hosoda F, Arai Y, Okada N, Shimizu H, Miyamoto M, Kitagawa N, Katai H, Taniguchi H, Yanagihara K, Imoto I, et al. Integrated genomic and functional analyses reveal glyoxalase I as a novel metabolic oncogene in human gastric cancer. Oncogene. 2015;34(9):1196–206.PubMedPubMed CentralGoogle Scholar
- Lachmann A, Xu H, Krishnan J, Berger SI, Mazloom AR, Ma’ayan A. ChEA: transcription factor regulation inferred from integrating genome-wide ChIP-X experiments. Bioinformatics. 2010;26(19):2438–44.PubMedPubMed CentralGoogle Scholar
- Cohen-Dvashi H, Ben-Chetrit N, Russell R, Carvalho S, Lauriola M, Nisani S, Mancini M, Nataraj N, Kedmi M, Roth L, et al. Navigator-3, a modulator of cell migration, may act as a suppressor of breast cancer progression. EMBO Mol Med. 2015;7(3):299–314.PubMedPubMed CentralGoogle Scholar
- Tayama S, Motohara T, Narantuya D, Li C, Fujimoto K, Sakaguchi I, Tashiro H, Saya H, Nagano O, Katabuchi H. The impact of EpCAM expression on response to chemotherapy and clinical outcomes in patients with epithelial ovarian cancer. Oncotarget. 2017;8(27):44312–25.PubMedPubMed CentralGoogle Scholar
- ATSDR. Toxicological profile for Mirex and Chlordecone. Atlanta: US Department of Health and Human Services, Public Health Service; 1995.Google Scholar
- Markkanen E. Not breathing is not an option: how to deal with oxidative DNA damage. DNA Repair. 2017;59:82–105.PubMedGoogle Scholar
- Shi B, Xue J, Yin H, Guo R, Luo M, Ye L, Shi Q, Huang X, Liu M, Sha J, et al. Dual functions for the ssDNA-binding protein RPA in meiotic recombination. PLoS Genet. 2019;15(2):e1007952.PubMedPubMed CentralGoogle Scholar
- Brown MS, Bishop DK. DNA strand exchange and RecA homologs in meiosis. Cold Spring Harb Perspect Biol. 2014;7(1):a016659.PubMedGoogle Scholar
- Mine-Hattab J, Rothstein R. DNA in motion during double-strand break repair. Trends Cell Biol. 2013;23(11):529–36.PubMedGoogle Scholar
- Silva N, Jantsch V. UNC-84: “LINC-ing” chromosome movement and double strand break repair. J Cell Biol. 2016;215(6):753–6.PubMedPubMed CentralGoogle Scholar
- Woglar A, Jantsch V. Chromosome movement in meiosis I prophase of Caenorhabditis elegans. Chromosoma. 2014;123(1–2):15–24.PubMedGoogle Scholar
- Williamson LM, Lees-Miller SP. Estrogen receptor alpha-mediated transcription induces cell cycle-dependent DNA double-strand breaks. Carcinogenesis. 2011;32(3):279–85.PubMedGoogle Scholar
- Teilmann SC, Byskov AG, Pedersen PA, Wheatley DN, Pazour GJ, Christensen ST. Localization of transient receptor potential ion channels in primary and motile cilia of the female murine reproductive organs. Mol Reprod Dev. 2005;71(4):444–52.PubMedGoogle Scholar
- Egeberg DL, Lethan M, Manguso R, Schneider L, Awan A, Jorgensen TS, Byskov AG, Pedersen LB, Christensen ST. Primary cilia and aberrant cell signaling in epithelial ovarian cancer. Cilia. 2012;1(1):15.PubMedPubMed CentralGoogle Scholar
- Lheureux S, Gourley C, Vergote I, Oza AM. Epithelial ovarian cancer. Lancet. 2019;393(10177):1240–53.PubMedGoogle Scholar
- Fagerberg L, Hallstrom BM, Oksvold P, Kampf C, Djureinovic D, Odeberg J, Habuka M, Tahmasebpoor S, Danielsson A, Edlund K, et al. Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics. Mol Cell Proteomics MCP. 2014;13(2):397–406.PubMedGoogle Scholar
- Kaufman OH, Lee K, Martin M, Rothhamel S, Marlow FL. rbpms2 functions in Balbiani body architecture and ovary fate. PLoS Genet. 2018;14(7):e1007489.PubMedPubMed CentralGoogle Scholar
- Baez-Vega PM, Echevarria Vargas IM, Valiyeva F, Encarnacion-Rosado J, Roman A, Flores J, Marcos-Martinez MJ, Vivas-Mejia PE. Targeting miR-21-3p inhibits proliferation and invasion of ovarian cancer cells. Oncotarget. 2016;7(24):36321–37.PubMedPubMed CentralGoogle Scholar
- Fang L, Chang HM, Cheng JC, Leung PC, Sun YP. TGF-beta1 downregulates StAR expression and decreases progesterone production through Smad3 and ERK1/2 signaling pathways in human granulosa cells. J Clin Endocrinol Metab. 2014;99(11):E2234–43.PubMedGoogle Scholar
- Myers M, Middlebrook BS, Matzuk MM, Pangas SA. Loss of inhibin alpha uncouples oocyte-granulosa cell dynamics and disrupts postnatal folliculogenesis. Dev Biol. 2009;334(2):458–67.PubMedPubMed CentralGoogle Scholar
- Kipp JL, Kilen SM, Bristol-Gould S, Woodruff TK, Mayo KE. Neonatal exposure to estrogens suppresses activin expression and signaling in the mouse ovary. Endocrinology. 2007;148(5):1968–76.PubMedGoogle Scholar
- Georges A, L’Hote D, Todeschini AL, Auguste A, Legois B, Zider A, Veitia RA. The transcription factor FOXL2 mobilizes estrogen signaling to maintain the identity of ovarian granulosa cells. eLife. 2014;3:1–19.Google Scholar
- Mo R, Rao SM, Zhu YJ. Identification of the MLL2 complex as a coactivator for estrogen receptor alpha. J Biol Chem. 2006;281(23):15714–20.PubMedGoogle Scholar
- Ding B, Cao Z, Hong R, Li H, Zuo X, Luo L, Li Y, Huang W, Li W, Zhang K, et al. WDR5 in porcine preimplantation embryos: expression, regulation of epigenetic modifications and requirement for early developmentdagger. Biol Reprod. 2017;96(4):758–71.PubMedGoogle Scholar
- Kim JY, Banerjee T, Vinckevicius A, Luo Q, Parker JB, Baker MR, Radhakrishnan I, Wei JJ, Barish GD, Chakravarti D. A role for WDR5 in integrating threonine 11 phosphorylation to lysine 4 methylation on histone H3 during androgen signaling and in prostate cancer. Mol Cell. 2014;54(4):613–25.PubMedPubMed CentralGoogle Scholar
- Mungamuri SK, Murk W, Grumolato L, Bernstein E, Aaronson SA. Chromatin modifications sequentially enhance ErbB2 expression in ErbB2-positive breast cancers. Cell Rep. 2013;5(2):302–13.PubMedPubMed CentralGoogle Scholar
- Riso V, Cammisa M, Kukreja H, Anvar Z, Verde G, Sparago A, Acurzio B, Lad S, Lonardo E, Sankar A, et al. ZFP57 maintains the parent-of-origin-specific expression of the imprinted genes and differentially affects non-imprinted targets in mouse embryonic stem cells. Nucleic Acids Res. 2016;44(17):8165–78.PubMedPubMed CentralGoogle Scholar
- Anvar Z, Cammisa M, Riso V, Baglivo I, Kukreja H, Sparago A, Girardot M, Lad S, De Feis I, Cerrato F, et al. ZFP57 recognizes multiple and closely spaced sequence motif variants to maintain repressive epigenetic marks in mouse embryonic stem cells. Nucleic Acids Res. 2016;44(3):1118–32.PubMedGoogle Scholar
- Klimczak M, Czerwinska P, Mazurek S, Sozanska B, Biecek P, Mackiewicz A, Wiznerowicz M. TRIM28 epigenetic corepressor is indispensable for stable induced pluripotent stem cell formation. Stem cell Res. 2017;23:163–72.PubMedGoogle Scholar
- Nagamachi A, Matsui H, Asou H, Ozaki Y, Aki D, Kanai A, Takubo K, Suda T, Nakamura T, Wolff L, et al. Haploinsufficiency of SAMD9L, an endosome fusion facilitator, causes myeloid malignancies in mice mimicking human diseases with monosomy 7. Cancer Cell. 2013;24(3):305–17.Google Scholar
- Ooga M, Funaya S, Hashioka Y, Fujii W, Naito K, Suzuki MG, Aoki F. Chd9 mediates highly loosened chromatin structure in growing mouse oocytes. Biochem Biophys Res Commun. 2018;500(3):583–8.PubMedGoogle Scholar
- Rinaldi VD, Bolcun-Filas E, Kogo H, Kurahashi H, Schimenti JC. The DNA damage checkpoint eliminates mouse oocytes with chromosome synapsis failure. Mol Cell 2017;67(6):1026–1036 e1022.PubMedPubMed CentralGoogle Scholar
- McLean CY, Bristor D, Hiller M, Clarke SL, Schaar BT, Lowe CB, Wenger AM, Bejerano G. GREAT improves functional interpretation of cis-regulatory regions. Nat Biotechnol. 2010;28(5):495–501.PubMedPubMed CentralGoogle Scholar
- Machanick P, Bailey TL. MEME-ChIP: motif analysis of large DNA datasets. Bioinformatics. 2011;27(12):1696–7.PubMedPubMed CentralGoogle Scholar
- Tanaka E, Bailey T, Grant CE, Noble WS, Keich U. Improved similarity scores for comparing motifs. Bioinformatics. 2011;27(12):1603–9.PubMedPubMed CentralGoogle Scholar