Major transcriptomic, epigenetic and metabolic changes underlie the pluripotency continuum in rabbit preimplantation embryos.

Despite the growing interest in the rabbit model for developmental and stem cell biology, the characterization of embryos at the molecular level is still poorly documented. We conducted a transcriptome analysis of rabbit preimplantation embryos from E2.7 (morula stage) to E6.6 (early primitive streak stage) using bulk and single-cell RNA-sequencing. In parallel, we studied oxidative phosphorylation and glycolysis, and analysed active and repressive epigenetic modifications during blastocyst formation and expansion. We generated a transcriptomic, epigenetic and metabolic map of the pluripotency continuum in rabbit preimplantation embryos, and identified novel markers of naive pluripotency that might be instrumental for deriving naive pluripotent stem cell lines. Although the rabbit is evolutionarily closer to mice than to primates, we found that the transcriptome of rabbit epiblast cells shares common features with those of humans and non-human primates.

Transcription of rRNA in early mouse embryos promotes chromatin reorganization and expression of major satellite repeats.

During the first cell cycles of early development, the chromatin of the embryo is highly reprogrammed while the embryonic genome starts its own transcription. The spatial organization of the genome is an important process that contributes to regulating gene transcription in time and space. It has, however, been poorly studied in the context of early embryos. To study the cause-and-effect link between transcription and spatial organization in embryos, we focused on ribosomal genes, which are silent initially but start to be transcribed in 2-cell mouse embryos. We demonstrated that ribosomal sequences and early unprocessed rRNAs are spatially organized in a very particular manner between 2-cell and 16-cell stage. By using drugs that interfere with ribosomal DNA transcription, we showed that this organization - which is totally different in somatic cells - depends on an active transcription of ribosomal genes and induces a unique chromatin environment that favors transcription of major satellite sequences once the 4-cell stage has been reached.

Maternal age affects equine day 8 embryo gene expression both in trophoblast and inner cell mass.

BACKGROUND: Breeding a mare until she is not fertile or even until her death is common in equine industry but the fertility decreases as the mare age increases. Embryo loss due to reduced embryo quality is partly accountable for this observation. Here, the effect of mare's age on blastocysts' gene expression was explored. Day 8 post-ovulation embryos were collected from multiparous young (YM, 6-year-old, N = 5) and older (OM, > 10-year-old, N = 6) non-nursing Saddlebred mares, inseminated with the semen of one stallion. Pure or inner cell mass (ICM) enriched trophoblast, obtained by embryo bisection, were RNA sequenced. Deconvolution algorithm was used to discriminate gene expression in the ICM from that in the trophoblast. Differential expression was analyzed with embryo sex and diameter as cofactors. Functional annotation and classification of differentially expressed genes and gene set enrichment analysis were also performed. RESULTS: Maternal aging did not affect embryo recovery rate, embryo diameter nor total RNA quantity. In both compartments, the expression of genes involved in mitochondria and protein metabolism were disturbed by maternal age, although more genes were affected in the ICM. Mitosis, signaling and adhesion pathways and embryo development were decreased in the ICM of embryos from old mares. In trophoblast, ion movement pathways were affected. CONCLUSIONS: This is the first study showing that maternal age affects gene expression in the equine blastocyst, demonstrating significant effects as early as 10 years of age. These perturbations may affect further embryo development and contribute to decreased fertility due to aging.

Three-dimensional analysis of nuclear heterochromatin distribution during early development in the rabbit.

Changes to the spatial organization of specific chromatin domains such as constitutive heterochromatin have been studied extensively in somatic cells. During early embryonic development, drastic epigenetic reprogramming of both the maternal and paternal genomes, followed by chromatin remodeling at the time of embryonic genome activation (EGA), have been observed in the mouse. Very few studies have been performed in other mammalian species (human, bovine, or rabbit) and the data are far from complete. During this work, we studied the three-dimensional organization of pericentromeric regions during the preimplantation period in the rabbit using specific techniques (3D-FISH) and tools (semi-automated image analysis). We observed that the pericentromeric regions (identified with specific probes for Rsat I and Rsat II genomic sequences) changed their shapes (from pearl necklaces to clusters), their nuclear localizations (from central to peripheral), as from the 4-cell stage. This reorganization goes along with histone modification changes and reduced amount of interactions with nucleolar precursor body surface. Altogether, our results suggest that the 4-cell stage may be a crucial window for events necessary before major EGA, which occurs during the 8-cell stage in the rabbit.

Different pre-implantation phenotypes of bovine blastocysts produced in vitro.

Embryo transfer in cattle is performed with blastocysts produced in vivo or in vitro using defined media. However, outdated systems such as those that use serum and co-culture remain of interest for research purposes. Here, we investigated the effect of additional culture time on in vitro-produced embryos. Specifically, we compared embryos that formed a blastocoel at different times after fertilisation to those that stayed in culture for up to two additional days with respect to their development in vivo after temporary transfer to oestrus-synchronised recipients. A pre-transfer set (D6, D6+1, D6+2, D7, D7+1, D8) was examined using microarray analyses and correlated with a post-transfer set that included two different days of transfer (D6-T6, D6+2-T8, D7+1-T8, D8-T8). All surviving conceptuses reached primitive-streak stages and filamentous sizes similarly to in vivo (D18) or in vitro controls (D7/T7). The recovery rate differed between D6 and D8 embryos that were immediately transferred (58 vs 25%). With an intermediate survival rate (33%), the D6 embryos with two additional days in culture produced nine times more IFN-tau (IFNT) at D18 than the D6 embryos that were immediately transferred. At the end of culture, D6 and D6+2 embryos displayed the highest number of gene expression differences. Despite a mortality of 40­60%, no signature was detectable in any of the transferred groups that would account for the embryos' fates. Initially reputed to be beneficial in producing more blastocysts, our culture system of B2 medium plus serum and co-culture generated blastocysts that were distinct from those developed in vivo (D7).

Eutherian mammals use diverse strategies to initiate X-chromosome inactivation during development.

X-chromosome inactivation (XCI) in female mammals allows dosage compensation for X-linked gene products between the sexes. The developmental regulation of this process has been extensively investigated in mice, where the X chromosome of paternal origin (Xp) is silenced during early embryogenesis owing to imprinted expression of the regulatory RNA, Xist (X-inactive specific transcript). Paternal XCI is reversed in the inner cell mass of the blastocyst and random XCI subsequently occurs in epiblast cells. Here we show that other eutherian mammals have very different strategies for initiating XCI. In rabbits and humans, the Xist homologue is not subject to imprinting and XCI begins later than in mice. Furthermore, Xist is upregulated on both X chromosomes in a high proportion of rabbit and human embryo cells, even in the inner cell mass. In rabbits, this triggers XCI on both X chromosomes in some cells. In humans, chromosome-wide XCI has not initiated even by the blastocyst stage, despite the upregulation of XIST. The choice of which X chromosome will finally become inactive thus occurs downstream of Xist upregulation in both rabbits and humans, unlike in mice. Our study demonstrates the remarkable diversity in XCI regulation and highlights differences between mammals in their requirement for dosage compensation during early embryogenesis.

Regulation of heat-inducible HSPA1A gene expression during maternal-to-embryo transition and in response to heat in in vitro-produced bovine embryos.

In in vitro-produced (IVP) bovine embryos, a burst in transcriptional activation of the embryonic genome (EGA) occurs at the 8-16-cell stage. To examine transcriptional regulation prior to EGA, notably in response to heat stress, we asked (1) whether the spontaneous expression of a luciferase transgene that is driven by the minimal mouse heat-shock protein 1b (hspa1b) gene promoter paralleled that of HSPA1A during EGA in IVP bovine embryo and (2) whether expression of the endogenous heat-inducible iHSPA group member HSPA1A gene and the hspa1b/luciferase transgene were induced by heat stress (HS) prior to EGA. Using two culture systems, we showed that luciferase activity levels rose during the 40-h long EGA-associated cell cycle. In contrast, iHSPA proteins were abundant in matured oocytes and in blastomeres from the two-cell to the 16-cell stages. However, normalised results detected a rise in the level of HSPA1A and luciferase mRNA during EGA, when transcription was required for their protein expression. Prior to EGA, HS-induced premature luciferase activity and transgene expression were clearly inhibited. We could not, however, establish whether this was also true for HSPA1A expression because of the decay of the abundant maternal transcripts prior to EGA. In bovine embryos, heat-induced expression of hspa1b/luciferase, and most likely of HSPA1A, was therefore strictly dependent on EGA. The level of the heat-shock transcription factor 1 molecules that were found in cell nuclei during embryonic development correlated better with the embryo's capacity for heat-shock response than with EGA-associated gene expression.

Alteration of the embryonic microenvironment and sex-specific responses of the preimplantation embryo related to a maternal high-fat diet in the rabbit model.

The maternal metabolic environment can be detrimental to the health of the offspring. In a previous work, we showed that maternal high-fat (HH) feeding in rabbit induced sex-dependent metabolic adaptation in the fetus and led to metabolic syndrome in adult offspring. As early development representing a critical window of susceptibility, in the present work we aimed to explore the effects of the HH diet on the oocyte, preimplantation embryo and its microenvironment. In oocytes from females on HH diet, transcriptomic analysis revealed a weak modification in the content of transcripts mainly involved in meiosis and translational control. The effect of maternal HH diet on the embryonic microenvironment was investigated by identifying the metabolite composition of uterine and embryonic fluids collected in vivo by biomicroscopy. Metabolomic analysis revealed differences in the HH uterine fluid surrounding the embryo, with increased pyruvate concentration. Within the blastocoelic fluid, metabolomic profiles showed decreased glucose and alanine concentrations. In addition, the blastocyst transcriptome showed under-expression of genes and pathways involved in lipid, glucose and amino acid transport and metabolism, most pronounced in female embryos. This work demonstrates that the maternal HH diet disrupts the in vivo composition of the embryonic microenvironment, where the presence of nutrients is increased. In contrast to this nutrient-rich environment, the embryo presents a decrease in nutrient sensing and metabolism suggesting a potential protective process. In addition, this work identifies a very early sex-specific response to the maternal HH diet, from the blastocyst stage.

PGE2 Supplementation of Oocyte Culture Media Improves the Developmental and Cryotolerance Performance of Bovine Blastocysts Derived From a Serum-Free in vitro Production System, Mirroring the Inner Cell Mass Transcriptome.

The culture media used throughout the in vitro production (IVP) of bovine embryos remain complex. The serum added to culture media in order to improve embryo development negatively impacts the cryotolerance of blastocysts. Periconceptional prostaglandin E2 (PGE2) signaling is known to exert prosurvival effects on in vitro-generated blastocysts. The purpose of the present study was to evaluate the effects on developmental and cryotolerance performance of a serum-free (SF) IVP system that included defined oocyte culture media supplemented or not with PGE2, versus serum-containing (SC) IVP. RNA-sequencing analysis was used to examine the gene expression of ICM derived under the different IVP conditions. We assessed the degree of cryotolerance of grade-I blastocysts during a three-day post-thaw culture by measuring survival and hatching rates, counting trophectoderm and inner cell mass (ICM) blastomere numbers. We also determined the proportion of ICM cells expressing octamer-binding transcription factor 4 protein (OCT4/POU5F1). We showed that grade-I blastocyst development rates under SF + PGE2 conditions were similar to those obtained under SC conditions, although the cleavage rate remained significantly lower. SC IVP conditions induced changes to ICM gene expression relative to several metabolic processes, catabolic activities, cell death and apoptosis. These alterations were associated with significantly higher levels of ICM cell death at day 7 post-fertilization, and lower survival and hatching rates after thawing. SF IVP conditions supplemented or not with PGE2 induced changes to ICM gene expression related to DNA replication, metabolism and double-strand break repair processes, and were associated with significantly larger ICM cell populations after thawing. SF + PGE2 IVP induced changes to ICM gene expression related to epigenetic regulation and were associated with a significantly higher proportion of ICM cells expressing OCT4. For the first time, our study thus offers a comprehensive analysis of the ICM transcriptome regulated by IVP culture conditions in terms of the cellular changes revealed during culture for three days after thawing.

Investigating the role of BCAR4 in ovarian physiology and female fertility by genome editing in rabbit.

Breast Cancer Anti-estrogen Resistance 4 (BCAR4) was previously characterised in bovine species as a gene preferentially expressed in oocytes, whose inhibition is detrimental to in vitro embryo development. But its role in oogenesis, folliculogenesis and globally fertility in vivo remains unknown. Because the gene is not conserved in mice, rabbits were chosen for investigation of BCAR4 expression and function in vivo. BCAR4 displayed preferential expression in the ovary compared to somatic organs, and within the ovarian follicle in the oocyte compared to somatic cells. The transcript was detected in follicles as early as the preantral stage. Abundance decreased throughout embryo development until the blastocyst stage. A lineage of genome-edited rabbits was produced; BCAR4 expression was abolished in follicles from homozygous animals. Females of wild-type, heterozygous and homozygous genotypes were examined for ovarian physiology and reproductive parameters. Follicle growth and the number of ovulations in response to hormonal stimulation were not significantly different between genotypes. Following insemination, homozygous females displayed a significantly lower delivery rate than their heterozygous counterparts (22 ± 7% vs 71 ± 11% (mean ± SEM)), while prolificacy was 1.8 ± 0.7 vs 6.0 ± 1.4 kittens per insemination. In conclusion, BCAR4 is not essential for follicular growth and ovulation but it contributes to optimal fertility in rabbits.

Differential correlations between changes to glutathione redox state, protein ubiquitination, and stress-inducible HSPA chaperone expression after different types of oxidative stress.

In primary bovine fibroblasts with an hspa1b/luciferase transgene, we examined the intensity of heat-shock response (HSR) following four types of oxidative stress or heat stress (HS), and its putative relationship with changes to different cell parameters, including reactive oxygen species (ROS), the redox status of the key molecules glutathione (GSH), NADP(H) NAD(H), and the post-translational protein modifications carbonylation, S-glutathionylation, and ubiquitination. We determined the sub-lethal condition generating the maximal luciferase activity and inducible HSPA protein level for treatments with hydrogen peroxide (H2O2), UVA-induced oxygen photo-activation, the superoxide-generating agent menadione (MN), and diamide (DA), an electrophilic and sulfhydryl reagent. The level of HSR induced by oxidative stress was the highest after DA and MN, followed by UVA and H2O2 treatments, and was not correlated to the level of ROS production nor to the extent of protein S-glutathionylation or carbonylation observed immediately after stress. We found a correlation following oxidative treatments between HSR and the level of GSH/GSSG immediately after stress, and the increase in protein ubiquitination during the recovery period. Conversely, HS treatment, which led to the highest HSR level, did not generate ROS nor modified or depended on GSH redox state. Furthermore, the level of protein ubiquitination was maximum immediately after HS and lower than after MN and DA treatments thereafter. In these cells, heat-induced HSR was therefore clearly different from oxidative stress-induced HSR, in which conversely early redox changes of the major cellular thiol predicted the level of HSR and polyubiquinated proteins.

Reprogramming of rabbit induced pluripotent stem cells toward epiblast and chimeric competency using Krüppel-like factors.

Rabbit induced pluripotent stem cells (rbiPSCs) possess the characteristic features of primed pluripotency as defined in rodents and primates. In the present study, we reprogrammed rbiPSCs using human Krüppel-like factors (KLFs) 2 and 4 and cultured them in a medium supplemented with fetal calf serum and leukemia inhibitory factor. These cells (designated rbEKA) were propagated by enzymatic dissociation for at least 30 passages, during which they maintained a normal karyotype. This new culturing protocol resulted in transcriptional and epigenetic reconfiguration, as substantiated by the expression of transcription factors and the presence of histone modifications associated with naïve pluripotency. Furthermore, microarray analysis of rbiPSCs, rbEKA cells, rabbit ICM cells, and rabbit epiblast showed that the global gene expression profile of the reprogrammed rbiPSCs was more similar to that of rabbit ICM and epiblast cells. Injection of rbEKA cells into 8-cell stage rabbit embryos resulted in extensive colonization of ICM in 9% early-blastocysts (E3.5), epiblast in 10% mid-blastocysts (E4.5), and embryonic disk in 1.4% pre-gastrulae (E6). Thus, these results indicate that KLF2 and KLF4 triggered the conversion of rbiPSCs into epiblast-like, embryo colonization-competent PSCs. Our results highlight some of the requirements to achieve bona fide chimeric competency.

Gene expression analysis in early embryos through reverse transcription quantitative PCR (RT-qPCR).

Real-time, reverse transcription quantitative PCR (RT-qPCR) is a highly sensitive and reproducible technology for the analysis of gene expression patterns. Its ability to detect minute quantities of nucleic acid from multifarious sources makes it an ideal technique for embryonic transcript quantification. However, complex cellular diversity and active transcriptome dynamics in early embryos necessitate particular caution to avoid erroneous results. This chapter is intended to outline basic methodology to design and execute RT-qPCR experiments in pre-implantation embryos.

Expression and localization of interleukin 1 beta and interleukin 1 receptor (type I) in the bovine endometrium and embryo.

The interleukin-1 (IL1) system likely mediates mammalian embryo-maternal communication. In cattle, we have reported that the uterine fluid of heifers carrying early embryos shows downregulated IL1 beta (IL1B), which could lead to reduced NFkB expression and dampening of maternal innate immune responses. In this work, we assessed the expression of IL 1 beta (IL1B) and its receptor, interleukin 1 receptor type I (IL1R1) in the bovine endometrium and embryos by RT-PCR, immunohistochemistry and Western blot at the time of blastocyst development. Day 8 endometrium, both collected from animals after transfer of day 5 embryos (ET) and sham transferred (ST), showed IL1B and IL1R1 mRNA transcription and protein co-localization. Similarly, day 8 blastocyst, from ET animals and entirely produced in vitro, showed IL1R1 mRNA transcription and IL1B and IL1R1 protein co-localization. IL1B mRNA was detected in the analyzed blastocysts, but at very low levels that precluded its quantification. IL1B and IL1R1 immunostaining was observed in luminal epithelial cells, glandular epithelium and stromal cells. The presence of embryos increased endometrial IL1B protein locally, while no differences regarding IL1R1 protein and IL1B and IL1R1 mRNA were detected. These results suggest that the early preimplantation bovine embryo in the maternal tract might interact with the maternal immune system through the IL1 system. Such a mechanism may allow the embryo to elicit local endometrial responses at early stages, which are required for the development of a receptive endometrium.

Bovine embryo cloning: Characterization of the recipient cytoplasts by phosphorylation patterns and kinase activities.

When in vitro-matured oocytes were enucleated, aged and kept at 10°C before reconstitution, the in vitro development of nuclear transfer embryos to the blastocyst stage did not differ from that obtained with in vitro fertilization. This suggests that these recipient cytoplasts constitute a suitable environment for the development of the nuclear transplant. The aim of the present study was to investigate, at the biochemical level, the result of the preparation of recipient oocytes, including enucleation, ageing and cooling. For this purpose the phosphorylation profiles of four groups of in vitro-matured bovine oocytes (aged oocytes, aged-cooled oocytes, enucleated-aged oocytes and enucleated-aged-cooled oocytes (recipient cytoplasts)) were analyzed. These recipient cytoplasts exhibited a phosphorylation profile similar to that of activated oocytes. Maturation promoting factor (MPF) activity, which was high in young metaphase II oocytes, in aged oocytes, in enucleated-aged oocytes and in aged-cooled oocytes, dropped to the basal level in enucleated-aged-cooled oocytes (recipient cytoplasts), while mitogen-activated protein kinase (MAPK) activity remained elevated. The combination of enucleation, ageing and cooling following oocyte in vitro maturation resulted in an interphase-like stage cytoplasm having a phosphorylation profile and low MPF activity similar to activated oocytes, but exhibiting high MAPK activity.

Sexual dimorphism of the feto-placental phenotype in response to a high fat and control maternal diets in a rabbit model.

Maternal environment during early developmental stages plays a seminal role in the establishment of adult phenotype. Using a rabbit model, we previously showed that feeding dams with a diet supplemented with 8% fat and 0.2% cholesterol (HH diet) from the prepubertal period and throughout gestation induced metabolic syndrome in adult offspring. Here, we examined the effects of the HH diet on feto-placental phenotype at 28 days post-coïtum (term = 31 days) in relation to earlier effects in the blastocyst (Day 6). At 28 days, both male and female HH fetuses were intrauterine growth retarded and dyslipidemic, with males more affected than females. Lipid droplets accumulated in the HH placentas' trophoblast, consistent with the increased concentrations in cholesteryl esters (3.2-fold), triacylglycerol (2.5-fold) and stored FA (2.12-fold). Stored FA concentrations were significantly higher in female compared to male HH placentas (2.18-fold, p<0.01), whereas triacylglycerol was increased only in HH males. Trophoblastic lipid droplet accumulation was also observed at the blastocyst stage. The expression of numerous genes involved in lipid pathways differed significantly according to diet both in term placenta and at the blastocyst stage. Among them, the expression of LXR-α in HH placentas was reduced in HH males but not females. These data demonstrate that maternal HH diet affects the blastocyst and induces sex-dependent metabolic adaptations in the placenta, which appears to protect female fetuses from developing severe dyslipidemia.

Alteration of DNA demethylation dynamics by in vitro culture conditions in rabbit pre-implantation embryos.

Alterations to DNA methylation have been attributed to in vitro culture and may affect normal embryo development. We chose to analyze DNA methylation reprogramming in the rabbit which, of the species with delayed transcriptional activation of the embryonic genome, allows easy comparisons between in vivo-developed (IVD) and in vitro-cultured (IVC) embryos. In this species, variations in DNA methylation had not previously been quantified, even in IVD embryos. IVD and IVC embryos were recovered at the 2, 4, 8 and 16-cell, morula and blastocyst stages. Immunostaining for 5-methyl-cytidine and normalization of the quantity of methylated DNA vs. the total DNA content were then performed. Our quantitative results evidenced DNA demethylation during pre-implantation development in both IVD and IVC embryos, but with different kinetics. Demethylation occurred earlier in vitro than in vivo between the 2 and 8-cell stages in IVC embryos, reaching its lowest level, while it only started at the 4-cell stage and ended at the 16-cell stage in IVD embryos. We also showed that an absence of serum from the culture medium significantly altered the degree of DNA demethylation. Finally, at the blastocyst stage, ICM was more methylated than the trophectoderm in all cases. Despite a morphological delay observed in in vitro cultured blastocysts, the difference in DNA methylation between ICM and trophectoderm cells appeared at the same time post-fertilization in IVD and IVC embryos, which may reflect another difference in the dynamics of DNA methylation during blastocyst formation. Our data thus clearly establish an effect of embryonic environment on DNA methylation reprogramming during pre-implantation development in a non-rodent species.

Expression of pluripotency master regulators during two key developmental transitions: EGA and early lineage specification in the bovine embryo.

Pluripotency genes are implicated in mouse embryonic genome activation (EGA) and pluripotent lineage specification. Moreover, their expression levels have been correlated with embryonic term development. In bovine, however, little information is available about dynamics of pluripotency genes during these processes. In this study, we charted quantitative and/or qualitative spatio-temporal expression patterns of transcripts and proteins of pluripotency genes (OCT4, SOX2 and NANOG) and mRNA levels of some of their downstream targets in bovine oocytes and early embryos. Furthermore, to correlate expression patterns of these genes with term developmental potential, we used cloned embryos, having similar in vitro but different full term development rates. Our findings affirm: firstly, the core triad of pluripotency genes is probably not implicated in bovine EGA since their proteins were not detected during pre-EGA phase, despite the transcripts for OCT4 and SOX2 were present. Secondly, an earlier ICM specification of transcripts and proteins of SOX2 and NANOG makes them pertinent candidates of bovine pluripotent lineage specification than OCT4. Thirdly, embryos with low term development potential have higher transcription rates; nevertheless, precarious balance between pluripotency genes is maintained. This balance presages normal in vitro development but, probably higher transcription rate disturbs it at later stage that abrogates term development.

Dynamics of DNA methylation levels in maternal and paternal rabbit genomes after fertilization.

The reprogramming of DNA methylation in early embryos has been considered to be essential for the reprogramming of differentiated parental genomes to totipotency, the transcription of embryonic genome activation (EGA) and subsequent development. However, its degree appears to differ as a function of species and it may be altered by the in vitro environment. While the rabbit is a pertinent model for species with a delayed EGA because both in vivo and in vitro developed embryos are easily available, the status of DNA methylation levels in both parental genomes after fertilization remains controversial. In order to generate precise data on the DNA methylation status in rabbit zygotes, we first of all defined five pronuclear (PN) stages during the first cell cycle and then classified in vivo and in vitro developed rabbit zygotes according to these PN stages. Using this classification we precisely quantified both methylated DNA and the total DNA content during the one cell stage. The quantification of methylated DNA, normalized for the total DNA content, showed that both pronuclei displayed distinct patterns of DNA methylation reprogramming. In the maternal pronucleus (MP) the methylation level remained constant throughout the one cell stage, thanks to maintenance methylation during the S phase. Conversely, in the paternal pronucleus (PP) partial demethylation occurred before replication, probably as a result of active DNA demethylation, while maintenance methylation subsequently occurred during the S phase. Interestingly, we showed that PP DNA methylation reprogramming was partially altered by the in vitro environment. Taken together, our approach evidenced that rabbit is one of the species displaying partial DNA demethylation in the PP, and for the first time demonstrated maintenance methylation activity in both pronuclei during the first S phase.