G9a regulates temporal preimplantation developmental program and lineage segregation in blastocyst

Version of Record

Accepted for publication after peer review and revision.

Download
Cite
Share
CommentOpen annotations (there are currently 0 annotations on this page).

Version of Record published: May 18, 2018 (This version)
Accepted Manuscript published: May 10, 2018 (Go to version)
Accepted: May 9, 2018
Received: November 14, 2017

1. Builds upon
Chromatin dynamics and the role of G9a in gene regulation and enhancer silencing during early mouse development

Jan J Zylicz, Sabine Dietmann ... M Azim Surani

Research Article Dec 9, 2015
Further reading

Abstract
Introduction
Results and discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

Early mouse development is regulated and accompanied by dynamic changes in chromatin modifications, including G9a-mediated histone H3 lysine 9 dimethylation (H3K9me2). Previously, we provided insights into its role in post-implantation development (Zylicz et al., 2015). Here we explore the impact of depleting the maternally inherited G9a in oocytes on development shortly after fertilisation. We show that G9a accumulates typically at 4 to 8 cell stage to promote timely repression of a subset of 4 cell stage-specific genes. Loss of maternal inheritance of G9a disrupts the gene regulatory network resulting in developmental delay and destabilisation of inner cell mass lineages by the late blastocyst stage. Our results indicate a vital role of this maternally inherited epigenetic regulator in creating conducive conditions for developmental progression and on cell fate choices.

https://doi.org/10.7554/eLife.33361.001

Introduction

The first developmental events in mouse are subject to regulation by information stored in the oocyte. The maternal inheritance in oocytes consists of mRNAs and proteins, which together direct rapid epigenetic reprogramming upon fertilisation, cell cleavage and activation of the zygotic genome (Ancelin et al., 2016; Li et al., 2010; Wasson et al., 2016).

The process of maternal-to-zygotic transition (MZT) in mice is first initiated in the late zygote, becoming more prominent at 2 cell stage (2C) on embryonic day (E) 1.5 (Golbus et al., 1973; Hamatani et al., 2004; Peaston et al., 2004). The initial wave of activation of many genes is followed by their repression within one or two cell cycles (Falco et al., 2007; Hamatani et al., 2004). How multiple genes are faithfully repressed during early embryogenesis to promote developmental progression remains unclear.

Extensive remodelling of the histone modifications and DNA methylation accompany division of blastomeres, which suggests a role for epigenetic regulatory mechanisms during development of the first embryonic lineages (Dahl et al., 2016; Liu et al., 2016; Smith et al., 2012; Wang et al., 2014; Wu et al., 2016; Zhang et al., 2016; Zheng et al., 2016). Outer trophectoderm (TE) cells are the first to be allocated together with the inner cell mass cells (ICM) at the blastocyst stage (E3.5). Shortly afterwards, by E4.0, the ICM segregates into primitive endoderm cells (PrE) and the pluripotent pre-Epiblast cells (Epi), which will give rise to the yolk sac and embryo proper respectively (Chazaud et al., 2006; Schrode et al., 2013).

We extend our previous study on the role of G9a-mediated H3K9me2 in mouse early post-implantation development (Zylicz et al., 2015), and examined the role of this histone methyltransferase during preimplantation development. We show that G9a (encoded by Ehmt2 gene) is maternally inherited and drives the accumulation of H3K9me2 at 4C and 8C stage, which accounts for timely repression of a subset of transcripts expressed at 4C. Severe disruption of the gene regulatory network follows upon maternal loss of G9a, resulting in developmental delay and destabilisation of ICM lineages, and frequent loss of embryos at the peri-implantation stage. Altogether, our results indicate that maternally-inherited G9a is crucial for regulating appropriate gene expression changes during preimplantation development.

Results and discussion

G9a and H3K9me2 accumulate at 4 and 8 cell stage

First, we investigated the H3K9me2 dynamics in early mouse preimplantation development. Immunofluorescence (IF) analysis of 2C (E1.5), 4C (E2.0), 8C (E2.5) and late blastocysts (E4.5) revealed progressive and significant accumulation of H3K9me2 at 4C and 8C stage; this was not the case at 2C or at E4.5 (Figure 1A,B). A more substantial enrichment follows in the epiblast of postimplantation embryos (Zylicz et al., 2015).

Figure 1 with 2 supplements see all

Download asset Open asset

H3K9me2 and G9a accumulate at 4- and 8 cell stage.

(A) Whole-mount IF staining for H3K9me2 (top panels) and G9a (bottom panels) in E1.5, E2.0, E2.5 and E4.5 embryos. DAPI intensity has been adjusted between time points for visualisation purposes (scale bar = 20 μm). IF signal is quantified (B) and visualised using box plots of median and interquartile range (IQR), with whiskers drawn 1.5xIQR away from the lower and upper quartiles. Data shows IF intensity normalised to DAPI for individual cells. At least 9 embryos were quantified for each time point. (*p<0.05 by Wilcoxon rank sum test). 2C: 2 cell stage; 4C: 4 cell stage; 8C: 8 cell stage; DAPI: 4',6-diamidino-2-phenylindole; H3K9me2: histone H3 lysine 9 dimethylation; IF: immunofluorescence; IQR: interquartile range. Also see Figure 1—figure supplements 1 and 2.

https://doi.org/10.7554/eLife.33361.002

The first wave of H3K9 dimethylation is consistent with increased level of G9a at 4C, and more significantly so at 8C stage (Figure 1A,B). What is more, G9a’s binding partner GLP also accumulates in 8C embryos (Figure 1—figure supplement 1). These results, in line with a previous immunofluorescence study (Li et al., 2013), indicate that even low levels of nuclear G9a at 4C stage are sufficient to initiate H3K9 dimethylation. Thus, following the burst of transcription at 2C stage, blastomeres accumulate substantial levels of repressive H3K9me2 mark, although its functional significance remains unclear.

G9a promotes developmental progression and primitive endoderm (PrE) segregation

To determine the role of G9a in preimplantation development, we first induced depletion of the maternal pool of both the Ehmt2 transcript and G9a protein in the oocyte (Figure 1—figure supplement 2). To do so, we used Zp3-Cre-expressing Ehmt2^F/- females (Ehmt2^F/+ for controls), in which conditional allele is recombined under Cre recombinase expression during oogenesis (de Vries et al., 2000), and crossed them with Ehmt2^+/- males. Whereas loss of both maternal and zygotic G9a (Ehmt2^M/Z) does not appear to grossly affect early development up to the E2.5 stage, these embryos show slight developmental delay and lack of substantial accumulation of H3K9me2 (Figure 2A,B, Figure 2—figure supplement 1A). Upon further development in vitro, the maternally depleted embryos do form blastocysts (E4.5), however with fewer Pou5f1-positive ICM cells (Figure 2—figure supplement 1B,C).

Figure 2 with 1 supplement see all

Download asset Open asset

Lack of maternal G9a leads to smaller blastocysts with fewer PrE cells.

(A) Whole-mount IF staining for G9a and H3K9me2 in E2.5 *Ehmt2^M/Z* and *Ehmt2^M/+* embryos. (Scale bar = 20 μm). (B) IF signal quantification for G9a and H3K9me2 from Figure 2A. Box plots show median and interquartile range (IQR), with whiskers drawn 1.5xIQR away from the lower and upper quartiles. Data shows IF intensity normalised to DAPI. At least six embryos were quantified for each genotype. (*p<0.05 in Wilcoxon rank sum test). (C) Whole-mount IF staining of E4.5 *Ehmt2^Mat* and *Ehmt2^Cntr* blastocysts using anti-CDX2 (TE, Blue), anti-SOX2 (Epi, Green) and anti-SOX17 (PrE, Red) antibodies. White arrows point towards nuclei devoid of staining for any lineage marker. (Scale bar = 20 μm). (**D–F**) Dot plots showing IF quantification from Figure 2C in relation to embryo genotypes. Each dot represents one embryo, 10 *Ehmt2^Mat* and 8 *Ehmt2^Cntr* embryos were quantified. (D) Total number of cells in an embryo. (E) Percentage of cells within the ICM (SOX2⁺ or SOX17⁺). (F) Percentage of cells within each lineage or showing no marker gene expression (SOX2^-SOX17^-CDX2^-). Line shows the median. (*p<0.05 in Wilcoxon rank sum test). CDX2: Caudal Type Homeobox 2; DAPI: 4',6-diamidino-2-phenylindole; *Ehmt2^Cntr*: control embryos with maternally inherited G9a; *Ehmt2^M/+* embryos maternally depleted but with zygotic expression of G9a; *Ehmt2^M/Z* embryos without both maternal and zygotic expression of G9a; *Ehmt2^Mat* embryos maternally depleted of G9a;Epi: pre-epiblast; H3K9me2: histone H3 lysine 9 dimethylation; ICM: inner cell mass; IF: immunofluorescence; IQR: interquartile range; PrE: primitive endoderm; SOX2: SRY box 2; SOX17: SRY box 17; TE: trophectoderm. Also see Figure 2—figure supplement 1.

https://doi.org/10.7554/eLife.33361.005

Next, we performed a detailed IF analysis of mutant and control embryos recovered at E4.5 for SOX2, SOX17 and CDX2, the critical markers of epiblast (Epi), primitive endoderm (PrE) and trophectoderm (TE), respectively (Figure 2C). We found no significant difference concerning lineage allocation between Ehmt2^M/Z and Ehmt2^M/+ embryos and thus decided to group them together as maternally depleted embryos (Ehmt2^Mat). Similarly, control embryos (Ehmt2^Cntr) of different genotypes (e.g. Ehmt2^+/- vs Ehmt2^-/-, data not shown) also do not reveal detectable phenotypic differences at E4.5. In contrast, Ehmt2^Mat compared to Ehmt2^Cntr blastocysts show significantly reduced cell numbers (80 vs 166 cells, respectively), and with fewer ICM (SOX2^+ve or SOX17^+ve) cells (25.0% vs 39.1% respectively) in line with our initial observation (Figure 2D,E, Figure 2—figure supplement 1B,C). We attribute the changes primarily to a relative reduction in PrE cells (7.9% vs 25.3%), as well as to a relative increase in the ICM cells that do not display expression of any of the three lineage markers (Figure 2F).

This observation suggests that maternal loss of G9a results in delayed development and lineage segregation within the ICM. The observed phenotype is however not entirely due to a developmental delay since we also see defects in lineage stabilisation within the ICM, potentially resulting in cell loss. Consistently, when comparing Ehmt2^Mat and Ehmt2^Cntr blastocysts we observed a significant increase in the levels γH2A.X, a hallmark of genetic stress and misregulation of the cell cycle (Figure 2—figure supplement 1B,C) (Turinetto and Giachino, 2015). Accumulation of γH2A.X also occurs in cells with decondensed chromatin potentially due to loss of H3K9me2 (Banáth et al., 2009). Notably, ‘knockdown’ of GLP, a binding partner of G9a, also results in defects in blastocyst, increased cell death, and loss of cells in the ICM (Huang et al., 2015). Thus, G9a-GLP complex likely promotes developmental progression in vivo, with a timely specification of the PrE and stabilisation of lineage choices. The significant reduction in the number of cells in mutant embryos as well as accumulation of γH2A.X indicate potential misregulation of the cell cycle as we reported in the postimplantation epiblast (Zylicz et al., 2015). Nevertheless, mutant embryos show defects in timely lineage segregation and stabilisation at the blastocyst stage, which indicates that epigenetic programming might create conditions conducive to cell fate choices. Previous studies have similarly indicated that arginine methylation by CARM1 and PRDM14, might also be important for the establishment of the Epi (Burton et al., 2013; Torres-Padilla et al., 2007). We thus hypothesise that maternal G9a is involved in setting up a stable transcriptional and epigenetic network allowing for timely developmental progression and lineage segregation at the blastocyst stage.

G9a represses a subset of 4C upregulated genes

To address the functional relevance of G9a in transcriptional regulation, we focused on the 8C stage (E2.5), when we detect accumulation of the highest levels of G9a and H3K9me2. Loss of maternal G9a at 8C has yet to result in an overt phenotype except for a mild developmental delay. To determine if there are underlying consequences of loss of G9a already at 8C, we performed single-embryo RNA-seq on ten control (Ehmt2^Cntr), and ten maternally depleted embryos (Ehmt2^Mat), which were morphologically indistinguishable (Figure 3—source data 1). Of note, we observed differential expression of only six genes between Ehmt2^M/Z and Ehmt2^M/+; we decided therefore to group them together. Maternally depleted G9a embryos however are largely transcriptionally distinct from the controls as shown by principal component analysis (PCA) (Figure 3—figure supplement 1A). The differential expression of 1467 genes accounts for the differences between the two groups, of which 44% become upregulated (Figure 3A, Figure 3—source data 2). Lack of enrichment for derepressed genes upon maternal loss of a transcriptional repressor implies that the transcriptome is affected prior to 8C, by which stage we observe substantial accumulation of secondary transcriptional effects. Interestingly, there were only four genes upregulated in both E2.5 Ehmt2^Mat and E6.25 Ehmt2^-/- embryos (Zylicz et al., 2015)(Chi2 p>0.77; Figure 3B), indicating that G9a plays a distinct role in pre- and post-implantation development.

Figure 3 with 3 supplements see all

Download asset Open asset

Maternal G9a represses a subset of genes induced at 4 cell stage.

(A) Scatter plot showing transcript expression levels in *Ehmt2^Cntr* and *Ehmt2^Mat* 8C (E2.5) stage embryos. Blue points are differentially expressed genes (adjusted p<0.05 in DEseq2). Shown is the average from ten biological replicates. (B) Venn diagram showing the overlap between upregulated genes upon maternal loss of G9a (*Ehmt2^Mat*) at 8C (E2.5) and those upregulated in zygotic *Ehmt2* deletion in E6.25 epiblast (Zylicz et al., 2015). (C) Expression profiles of genes within specific clusters during wildtype development from 2C to 8C (left panels). Dotted line represents mean expression within the cluster. Data from GSE22182 (Tang et al., 2011) was used and within all expressed genes 14 specific clusters were identified. Shown are clusters most enriched for genes becoming derepressed (Clusters 1 and 6) or downregulated (Cluster 11) at 8C in Ehmt2^Mat. Right panel of bar plots shows the percentage of genes upregulated (Green) or downregulated (Red) at 8C in Ehmt2^Mat, which belong to identified clusters. Significance of enrichments was calculated using Chi2 test (* 10⁻⁵ < p < 0.05; ** 10⁻¹⁰ < p < 10⁻⁵; ***p<10⁻¹⁰). (D) Fold expression changes of PrE (*Fgfr1, Gata6*, *Gata4*) and TE (*Tead4, Tfap2c*) marker genes in *Ehmt2^Mat*compared to *Ehmt2^Cntr* 8C embryos. Box plots show median and interquartile range (IQR), with whiskers drawn 1.5xIQR away from the lower and upper quartiles. Each dot is a single embryo and data was normalised to median expression levels in *Ehmt2^Cntr*. (*adjusted p<0.05 in DEseq2). 2C: 2 cell embryo, 4C: 4 cell embryo; 8C: 8 cell embryo; *Ehmt2^Cntr*: control embryos with maternally inherited G9a; *Ehmt2^Mat* embryos maternally depleted of G9a; *Fgfr1*: fibroblast growth factor receptor 1; *Gata4*: GATA binding protein 4; *Gata6*: GATA binding protein 6; IQR: interquartile range; *Tcfap2c*: Transcription factor AP-2 gamma; TE: trophectoderm; *Tead4*: TEA domain transcription factor 4. Also see Figure 3—figure supplements 1–3 and Figure 3—source data 1–3.

https://doi.org/10.7554/eLife.33361.007

Figure 3—source data 1 List of all single-embryo RNAseq samples sequenced. Shown are numbers of uniquely mapped reads as well as the genotype of each embryo.: https://doi.org/10.7554/eLife.33361.011
Download elife-33361-fig3-data1-v2.xlsx
Figure 3—source data 2 List of differentially expressed genes between Ehmt2^Cntr and Ehmt2^Mat at 8C stage embryos. Differentially expressed genes were identified using DEseq2 with adjusted p-value<0.05.: https://doi.org/10.7554/eLife.33361.012
Download elife-33361-fig3-data2-v2.xlsx
Figure 3—source data 3 List of enriched GO terms in genes upregulated or downregulated in 8C Ehmt2^Mat embryos. GO term enrichment for biological processes was calculated using DAVID software with minimum five genes in a category and EASE p-value<0.05.: https://doi.org/10.7554/eLife.33361.013
Download elife-33361-fig3-data3-v2.xlsx

For further insight into the role of G9a in regulating early development, we performed Gene Ontology (GO) enrichment analysis (Figure 3—source data 3). Consistent with the reduced blastocyst size and loss of ICM cells at E4.5, the GO terms indicated down regulation of genes involved in cell cycle (p<10⁻⁸), transcription (p<10⁻⁴) and stem cell population maintenance (p<10⁻⁴), in Ehmt2^Mat 8C embryos (Figure 3—figure supplement 1B). However, the GO term most strikingly enriched was chromatin silencing (p<10⁻¹¹), which is entirely due to reduced levels of 18 transcripts coding for H2A, MacroH2A and H2A.J (Figure 3—figure supplement 1B,C). Invariably, these histones are already upregulated in a subset of Ehmt2^Cntr embryos by E2.5, which is not the case in the mutants (Figure 3—figure supplement 1D,E)(Wu et al., 2014). Thus, despite being morphologically indistinguishable from controls, Ehmt2^Mat 8C embryos deviate from them concerning developmental progression. On the other hand, there was high enrichment of genes for GO terms such as rRNA processing (p<10⁻²³) and mRNA processing (p<10⁻¹⁷) indicating that G9a preferentially represses genes linked to RNA metabolism (Figure 3—figure supplement 1F).

To put these results in the context of the developmental progression, we have integrated our results with an independent single-cell RNAseq dataset of wildtype mouse preimplantation embryos (Tang et al., 2011). Firstly, we performed clustering of genes using Gaussian mixture fitting based on their dynamic expression between 2C and 8C stages (Wang et al., 2012). To describe the expression dynamics, we found an optimal number of 14 clusters of genes. Next, we overlapped genes within clusters with those upregulated or downregulated at the 8C stage in Ehmt2^Mat embryos (Figure 3C). Strikingly, of the upregulated genes, 28% and 30% were part of cluster 1 and 6, respectively (Chi² <10⁻⁶⁸ and Chi² <10⁻⁴⁸). Genes within those clusters are transiently upregulated at the 4C stage. On the other hand, cluster 11 of stably expressed genes, is most enriched for genes downregulated upon loss of maternal G9a (12% of genes, Chi2 <10⁻¹⁵, Figure 3C). To further validate the specific derepression of 4C genes in mutants, we have also analysed gene expression changes between 2C and 4C, as well as 4C and 8C stages in wildtype embryos (Figure 3—figure supplement 2A). Our analysis revealed that genes upregulated in Ehmt2^Mat 8C embryos (green box), show a significant increase in expression at 4C stage, which are silenced later in wildtype embryos. Our findings indicate that during preimplantation development, G9a mediates silencing of a specific set of genes, which are transiently upregulated at 4C stages. Whereas there is significant accumulation of H3K9me2 globally as judged by IF analysis at the 4C stage (Figure 1A,B), the precise distribution of the mark within the genome cannot be ascertained by this observation. We propose that it is only at the 8C stage that the increasing levels of G9a levels allow H3K9me2 deposition at the 4C-specific genes.

Next, to explore the roots of lineage destabilisation in maternally depleted embryos, we focused on the misexpression of transcriptional factors (TFs) in our RNAseq dataset (n = 71) (Figure 3—figure supplement 2B). Consistent with fewer PrE and more TE cells, there was a reduction in the levels of PrE specifiers (Gata6, Fgfr1), and an increase in the expression of TE markers (Tead4, Tcfap2c) in maternally deleted mutants compared to controls (Figure 3D). Intriguingly, expression of Gata4, normally seen in late blastocysts, already shows upregulation in mutant embryos (Plusa et al., 2008) (Figure 3D). These results indicate that maternal loss of G9a results in destabilisation of specific gene regulatory networks involved in pre-implantation development. Efficient setting up of such transcriptional circuitry as early as the 8C stage seems necessary for subsequent timely lineage segregation. However, it remains unclear which subset of the differentially expressed genes directly results in PrE defects. Indeed, this might be linked to increased expression of a critical trophectoderm specifier Tead4 (Yagi et al., 2007), or downregulation of vital PrE regulators Fgfr1 (Kang et al., 2017; Molotkov et al., 2017) and Gata6 (Schrode et al., 2014). Alternatively, ICM cells might show increased sensitivity to the deregulation of such processes as RNA metabolism and cell cycle. Importantly, derepression of transposable elements (TE) might also result in increased DNA damage and cell death (Ancelin et al., 2016).

To address the latter hypothesis, we performed differential TE expression analysis on our RNA-seq dataset. Of particular interest are murine endogenous retrovirus-like elements (MuERV-L), which were shown to be derepressed upon loss of G9a in embryonic stem cells (ESCs)(Maksakova et al., 2013). This class of TEs is transiently expressed at the 2C stage (Kigami et al., 2003; Macfarlan et al., 2012; Peaston et al., 2004) and regulates zygotic genome activation, transcriptional networks and developmenal progression (Huang et al., 2017; Kigami et al., 2003). Our TE analysis of mutant embryos however revealed only a very mild reactivation of the ERV-L family of TEs, which contains MuERV-L (Figure 3—figure supplement 3). Thus, G9a is dispensible for timely repression of these retrotransposons. We have however observed persistent H3K9me2 staning at DAPI-dense foci in Ehmt2^M/Z at the 8C stage (Figure 2A), suggesting that another histone methyltransferase is likely to deposit silent chromatin marks at repetetive elements in the absence of G9a. Indeed, a recent study suggests that a loss of SETDB1, an H3K9 methyltransferase, results in upregulation of both ERV and LINE1 elements in oocytes (Kim et al., 2016).

In conclusion, the role of G9a during preimplantation development is distinct from that observed during postimplantation development. Notably, G9a is required for the repression of a specific set of transiently upregulated 4C genes in preimplantation embryos. Moreover, maternal G9a allows for the appropriate setting of transcriptional circuits, which promote developmental progression, timely PrE specification and stabilisation of ICM lineages. Together with PRDM14 and arginine methylation reported previously (Burton et al., 2013; Torres-Padilla et al., 2007), our findings suggest that the earliest epigenetic programming events at the onset of development are involved in creating a competent environment for cell fate choices that ensue. These epigenetic regulators, amongst others, ensure establishment of stable transcriptional circuitry, which instructs lineage segregation later in development.

Materials and methods

Key resources table

Reagent type (species) or resource	Designation	Source or reference
Strain (mouse)	Ehmt2Flox/Flox	PMID:17707231
Strain (mouse)	Zp3-Cre	PMID: 10686600

GO term enrichment analysis were performed exactly as in (Zylicz et al., 2015).

Mouse breading, embryo collection

Request a detailed protocol

Timed natural matings were used for all experiments unless otherwise stated. Noon of the day when the vaginal plugs of mated females were identified was scored as E0.5. For Ehmt2 matings a published conditional allele was used (Sampath et al., 2007). To obtain Ehmt2^Cntr or Ehmt2^Mat embryos, Ehmt2^F/+ Zp3-Cre^+ve or Ehmt2^F/- Zp3-Cre^+ve females were used respectively (de Vries et al., 2000). When stated a ΔPE-Pou5f1-EGFP reporter line was crossed in (GGOF) (Yeom et al., 1996). All husbandry and experiments involving mice were carried out according to the local ethics committee and were performed in a facility designated by the Home Office.

Immunofluorescence

Request a detailed protocol

Embryos were treated as previously described (Nichols and Smith, 2009). Primary antibodies used are as follows: anti-CDX2 (Biogenex, clone CDX2-88), anti-H3K9me2 (Abcam, UK, ab1220), anti-GFP (Nacalai tesque, Japan, GF090R), anti-G9a (Cell Signaling, MA, 68851T), anti-GLP (Research and Diagnostic Systems, MN, PP-B0422-00), anti-SOX2 (Abcam, UK, ab92494), anti-SOX17(Research and Diagnostic Systems, MN, AF1924). Mean nuclear intensities of IF and DAPI signal were quantified using mageJ and corrected for the staining background. As nuclear size is changing between stages all IF measurements were normalised to DAPI signal as a proxy for DNA content.

Single-Embryo RNAseq

Request a detailed protocol

The embryos used were from natural matings and were morphologically assessed to ensure only viable samples were collected. cDNA was prepared and amplified as earlier described (Tang et al., 2010). Illumina libraries were prepared as published (Huang et al., 2017). Single-end 50 bp sequencing was performed with HiSeq4000 (Illumina, San Diego, CA). RNA-seq reads were adapter- and quality-trimmed, and aligned with Tophat2 (Kim et al., 2013) against the mouse reference (GRCm38/mm10) genome. Read counts per ENSEMBL transcript were obtained by SeqMonk. Differential expression was evaluated with the DESeq2 package (Love et al., 2014). Gene was deemed differentially expressed when p-value<0.05 after Benjamini and Hochberg correction. For clustering the dynamics of gene expression in distinct environments we have used an R package GeneClusterNet based on Gaussian mixture fitting (Wang et al., 2012). For this analysis data was downloaded from GEO (GSE22182) (Tang et al., 2011). Optimal numbers of clusters (14) was identified by finding minimal Bayesian Information Criterion. For TE analysis, RNA-seq reads were aligned with bowtie (options: ‘-m 1 – v1 --best --strata’) selecting for uniquely mapping reads only. RepeatMasker annotations of individual TE elements were downloaded from the UCSC Table Browser. Read counts per TE elements were obtained by featureCounts (http://bioinf.wehi.edu.au/featureCounts). Data is available under GSE106790.

Data availability

Sequencing data have been deposited in GEO under accession codes GSE106790.

The following data sets were generated

1. Zylicz JJ
2. Borensztein M
3. Dietmann S
4. Huang Y
5. Lee C
6. Surani MA
(2017) G9a regulates temporal preimplantation developmental program and lineage segregation in blastocyst
Publicly available at the NCBI Gene Expression Omnibus (accession no: GSE106790).

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE106790

The following previously published data sets were used

1. Tang F
2. Barbacioru C
3. Lao K
4. Surani MA
(2010) Global Deterministic and Stochastic Allelic Specific Gene Expression in Single Blastomeres of Mouse Early Embryos
Publicly available at the NCBI Gene Expression Omnibus (accession no: GSE22182).

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE22182
1. Wu J
2. Huang B
3. Chen H
4. Xie W
(2016) The landscape of accessible chromatin in mammalian pre-implantation embryos (RNA-Seq)
Publicly available at the NCBI Gene Expression Omnibus (accession no: GSE66582).

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE66582

References

1. Ancelin K
2. Syx L
3. Borensztein M
4. Ranisavljevic N
5. Vassilev I
6. Briseño-Roa L
7. Liu T
8. Metzger E
9. Servant N
10. Barillot E
11. Chen CJ
12. Schüle R
13. Heard E
(2016) Maternal LSD1/KDM1A is an essential regulator of chromatin and transcription landscapes during zygotic genome activation
eLife 5:e08851.
https://doi.org/10.7554/eLife.08851
- PubMed
- Google Scholar
(2009) Explanation for excessive DNA single-strand breaks and endogenous repair foci in pluripotent mouse embryonic stem cells
Experimental Cell Research 315:1505–1520.
https://doi.org/10.1016/j.yexcr.2008.12.007
- PubMed
- Google Scholar
(2013) Single-cell profiling of epigenetic modifiers identifies PRDM14 as an inducer of cell fate in the mammalian embryo
Cell Reports 5:687–701.
https://doi.org/10.1016/j.celrep.2013.09.044
- PubMed
- Google Scholar
(2006) Early lineage segregation between epiblast and primitive endoderm in mouse blastocysts through the Grb2-MAPK pathway
Developmental Cell 10:615–624.
https://doi.org/10.1016/j.devcel.2006.02.020
- PubMed
- Google Scholar
1. Dahl JA
2. Jung I
3. Aanes H
4. Greggains GD
5. Manaf A
6. Lerdrup M
7. Li G
8. Kuan S
9. Li B
10. Lee AY
11. Preissl S
12. Jermstad I
13. Haugen MH
14. Suganthan R
15. Bjørås M
16. Hansen K
17. Dalen KT
18. Fedorcsak P
19. Ren B
20. Klungland A
(2016) Broad histone H3K4me3 domains in mouse oocytes modulate maternal-to-zygotic transition
Nature 537:548–552.
https://doi.org/10.1038/nature19360
- PubMed
- Google Scholar
1. de Vries WN
2. Binns LT
3. Fancher KS
4. Dean J
5. Moore R
6. Kemler R
7. Knowles BB
(2000) Expression of Cre recombinase in mouse oocytes: a means to study maternal effect genes
Genesis 26:110–112.
https://doi.org/10.1002/(SICI)1526-968X(200002)26:2<110::AID-GENE2>3.0.CO;2-8
- PubMed
- Google Scholar
1. Falco G
2. Lee SL
3. Stanghellini I
4. Bassey UC
5. Hamatani T
6. Ko MS
(2007) Zscan4: a novel gene expressed exclusively in late 2-cell embryos and embryonic stem cells
Developmental Biology 307:539–550.
https://doi.org/10.1016/j.ydbio.2007.05.003
- PubMed
- Google Scholar
(1973) The effects of inhibitors of RNA synthesis (alpha-amanitin and actinomycin D) on preimplantation mouse embryogenesis
Journal of Experimental Zoology 186:207–216.
https://doi.org/10.1002/jez.1401860211
- PubMed
- Google Scholar
1. Hamatani T
2. Carter MG
3. Sharov AA
4. Ko MS
(2004) Dynamics of global gene expression changes during mouse preimplantation development
Developmental Cell 6:117–131.
https://doi.org/10.1016/S1534-5807(03)00373-3
- PubMed
- Google Scholar
1. Huang XJ
2. Ma X
3. Wang X
4. Zhou X
5. Li J
6. Sun SC
7. Liu H
(2015) Involvement of G9A-like protein (GLP) in the development of mouse preimplantation embryos in vitro
Reproduction, fertility, and development.
https://doi.org/10.1071/RD14341
- PubMed
- Google Scholar
1. Huang Y
2. Kim JK
3. Do DV
4. Lee C
5. Penfold CA
6. Zylicz JJ
7. Marioni JC
8. Hackett JA
9. Surani MA
(2017) Stella modulates transcriptional and endogenous retrovirus programs during maternal-to-zygotic transition
eLife 6:e22345.
https://doi.org/10.7554/eLife.22345
- PubMed
- Google Scholar
(2017) Lineage establishment and progression within the inner cell mass of the mouse blastocyst requires FGFR1 and FGFR2
Developmental Cell 41:496–510.
https://doi.org/10.1016/j.devcel.2017.05.003
- PubMed
- Google Scholar
1. Kigami D
2. Minami N
3. Takayama H
4. Imai H
(2003) MuERV-L is one of the earliest transcribed genes in mouse one-cell embryos
Biology of Reproduction 68:651–654.
https://doi.org/10.1095/biolreprod.102.007906
- PubMed
- Google Scholar
1. Kim D
2. Pertea G
3. Trapnell C
4. Pimentel H
5. Kelley R
6. Salzberg SL
(2013) TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions
Genome Biology 14:R36.
https://doi.org/10.1186/gb-2013-14-4-r36
- PubMed
- Google Scholar
1. Kim J
2. Zhao H
3. Dan J
4. Kim S
5. Hardikar S
6. Hollowell D
7. Lin K
8. Lu Y
9. Takata Y
10. Shen J
11. Chen T
(2016) Maternal Setdb1 is required for meiotic progression and preimplantation development in mouse
PLOS Genetics 12:e1005970.
https://doi.org/10.1371/journal.pgen.1005970
- PubMed
- Google Scholar
1. Li B
2. Tang N
3. Chen S
4. Li X
5. Huang X
6. Wang X
7. Sun F
(2013) G9a co-localized with histone H3 lysine 9 monomethylation but not dimethylation in a nuclear membrane-dependent manner during mouse preimplantation embryo development
Journal of Assisted Reproduction and Genetics 30:441–448.
https://doi.org/10.1007/s10815-012-9911-y
- PubMed
- Google Scholar
1. Li L
2. Zheng P
3. Dean J
(2010) Maternal control of early mouse development
Development 137:859–870.
https://doi.org/10.1242/dev.039487
- PubMed
- Google Scholar
1. Liu X
2. Wang C
3. Liu W
4. Li J
5. Li C
6. Kou X
7. Chen J
8. Zhao Y
9. Gao H
10. Wang H
11. Zhang Y
12. Gao Y
13. Gao S
(2016) Distinct features of H3K4me3 and H3K27me3 chromatin domains in pre-implantation embryos
Nature 537:558–562.
https://doi.org/10.1038/nature19362
- PubMed
- Google Scholar
1. Love MI
2. Huber W
3. Anders S
(2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2
Genome Biology 15:550.
https://doi.org/10.1186/s13059-014-0550-8
- PubMed
- Google Scholar
1. Macfarlan TS
2. Gifford WD
3. Driscoll S
4. Lettieri K
5. Rowe HM
6. Bonanomi D
7. Firth A
8. Singer O
9. Trono D
10. Pfaff SL
(2012) Embryonic stem cell potency fluctuates with endogenous retrovirus activity
Nature 487:57–63.
https://doi.org/10.1038/nature11244
- PubMed
- Google Scholar
1. Maksakova IA
2. Thompson PJ
3. Goyal P
4. Jones SJ
5. Singh PB
6. Karimi MM
7. Lorincz MC
(2013) Distinct roles of KAP1, HP1 and G9a/GLP in silencing of the two-cell-specific retrotransposon MERVL in mouse ES cells
Epigenetics & Chromatin 6:15.
https://doi.org/10.1186/1756-8935-6-15
- PubMed
- Google Scholar
1. Molotkov A
2. Mazot P
3. Brewer JR
4. Cinalli RM
5. Soriano P
(2017) Distinct requirements for FGFR1 and FGFR2 in primitive endoderm development and exit from pluripotency
Developmental Cell 41:511–526.
https://doi.org/10.1016/j.devcel.2017.05.004
- PubMed
- Google Scholar
1. Nichols J
2. Smith A
(2009) Naive and primed pluripotent states
Cell Stem Cell 4:487–492.
https://doi.org/10.1016/j.stem.2009.05.015
- PubMed
- Google Scholar
(2004) Retrotransposons regulate host genes in mouse oocytes and preimplantation embryos
Developmental Cell 7:597–606.
https://doi.org/10.1016/j.devcel.2004.09.004
- PubMed
- Google Scholar
(2008) Distinct sequential cell behaviours direct primitive endoderm formation in the mouse blastocyst
Development 135:3081–3091.
https://doi.org/10.1242/dev.021519
- PubMed
- Google Scholar
1. Sampath SC
2. Marazzi I
3. Yap KL
4. Sampath SC
5. Krutchinsky AN
6. Mecklenbräuker I
7. Viale A
8. Rudensky E
9. Zhou MM
10. Chait BT
11. Tarakhovsky A
(2007) Methylation of a histone mimic within the histone methyltransferase G9a regulates protein complex assembly
Molecular Cell 27:596–608.
https://doi.org/10.1016/j.molcel.2007.06.026
- PubMed
- Google Scholar
(2014) GATA6 levels modulate primitive endoderm cell fate choice and timing in the mouse blastocyst
Developmental Cell 29:454–467.
https://doi.org/10.1016/j.devcel.2014.04.011
- PubMed
- Google Scholar
(2013) Anatomy of a blastocyst: cell behaviors driving cell fate choice and morphogenesis in the early mouse embryo
Genesis 51:219–233.
https://doi.org/10.1002/dvg.22368
- PubMed
- Google Scholar
1. Smith ZD
2. Chan MM
3. Mikkelsen TS
4. Gu H
5. Gnirke A
6. Regev A
7. Meissner A
(2012) A unique regulatory phase of DNA methylation in the early mammalian embryo
Nature 484:339–344.
https://doi.org/10.1038/nature10960
- PubMed
- Google Scholar
1. Tang F
2. Barbacioru C
3. Nordman E
4. Bao S
5. Lee C
6. Wang X
7. Tuch BB
8. Heard E
9. Lao K
10. Surani MA
(2011) Deterministic and stochastic allele specific gene expression in single mouse blastomeres
PLoS One 6:e21208.
https://doi.org/10.1371/journal.pone.0021208
- PubMed
- Google Scholar
1. Tang F
2. Barbacioru C
3. Nordman E
4. Li B
5. Xu N
6. Bashkirov VI
7. Lao K
8. Surani MA
(2010) RNA-Seq analysis to capture the transcriptome landscape of a single cell
Nature Protocols 5:516–535.
https://doi.org/10.1038/nprot.2009.236
- PubMed
- Google Scholar
(2007) Histone arginine methylation regulates pluripotency in the early mouse embryo
Nature 445:214–218.
https://doi.org/10.1038/nature05458
- PubMed
- Google Scholar
1. Turinetto V
2. Giachino C
(2015) Multiple facets of histone variant H2AX: a DNA double-strand-break marker with several biological functions
Nucleic Acids Research 43:2489–2498.
https://doi.org/10.1093/nar/gkv061
- PubMed
- Google Scholar
1. Wang L
2. Zhang J
3. Duan J
4. Gao X
5. Zhu W
6. Lu X
7. Yang L
8. Zhang J
9. Li G
10. Ci W
11. Li W
12. Zhou Q
13. Aluru N
14. Tang F
15. He C
16. Huang X
17. Liu J
(2014) Programming and inheritance of parental DNA methylomes in mammals
Cell 157:979–991.
https://doi.org/10.1016/j.cell.2014.04.017
- PubMed
- Google Scholar
1. Wang Y
2. Xu M
3. Wang Z
4. Tao M
5. Zhu J
6. Wang L
7. Li R
8. Berceli SA
9. Wu R
(2012) How to cluster gene expression dynamics in response to environmental signals
Briefings in Bioinformatics 13:162–174.
https://doi.org/10.1093/bib/bbr032
- PubMed
- Google Scholar
1. Wasson JA
2. Simon AK
3. Myrick DA
4. Wolf G
5. Driscoll S
6. Pfaff SL
7. Macfarlan TS
8. Katz DJ
(2016) Maternally provided LSD1/KDM1A enables the maternal-to-zygotic transition and prevents defects that manifest postnatally
eLife 5:e08848.
https://doi.org/10.7554/eLife.08848
- PubMed
- Google Scholar
1. Wu BJ
2. Dong FL
3. Ma XS
4. Wang XG
5. Lin F
6. Liu HL
(2014) Localization and expression of histone H2A variants during mouse oogenesis and preimplantation embryo development
Genetics and Molecular Research 13:5929–5939.
https://doi.org/10.4238/2014.August.7.8
- PubMed
- Google Scholar
1. Wu J
2. Huang B
3. Chen H
4. Yin Q
5. Liu Y
6. Xiang Y
7. Zhang B
8. Liu B
9. Wang Q
10. Xia W
11. Li W
12. Li Y
13. Ma J
14. Peng X
15. Zheng H
16. Ming J
17. Zhang W
18. Zhang J
19. Tian G
20. Xu F
21. Chang Z
22. Na J
23. Yang X
24. Xie W
(2016) The landscape of accessible chromatin in mammalian preimplantation embryos
Nature 534:652–657.
https://doi.org/10.1038/nature18606
- PubMed
- Google Scholar
(2007) Transcription factor TEAD4 specifies the trophectoderm lineage at the beginning of mammalian development
Development 134:3827–3836.
https://doi.org/10.1242/dev.010223
- PubMed
- Google Scholar
1. Yeom YI
2. Fuhrmann G
3. Ovitt CE
4. Brehm A
5. Ohbo K
6. Gross M
7. Hübner K
8. Schöler HR
(1996)
Germline regulatory element of Oct-4 specific for the totipotent cycle of embryonal cells

Development 122:881–894.
- PubMed
- Google Scholar
1. Zhang B
2. Zheng H
3. Huang B
4. Li W
5. Xiang Y
6. Peng X
7. Ming J
8. Wu X
9. Zhang Y
10. Xu Q
11. Liu W
12. Kou X
13. Zhao Y
14. He W
15. Li C
16. Chen B
17. Li Y
18. Wang Q
19. Ma J
20. Yin Q
21. Kee K
22. Meng A
23. Gao S
24. Xu F
25. Na J
26. Xie W
(2016) Allelic reprogramming of the histone modification H3K4me3 in early mammalian development
Nature 537:553–557.
https://doi.org/10.1038/nature19361
- PubMed
- Google Scholar
1. Zheng H
2. Huang B
3. Zhang B
4. Xiang Y
5. Du Z
6. Xu Q
7. Li Y
8. Wang Q
9. Ma J
10. Peng X
11. Xu F
12. Xie W
(2016) Resetting epigenetic memory by reprogramming of histone modifications in mammals
Molecular Cell 63:1066–1079.
https://doi.org/10.1016/j.molcel.2016.08.032
- PubMed
- Google Scholar
1. Zylicz JJ
2. Dietmann S
3. Günesdogan U
4. Hackett JA
5. Cougot D
6. Lee C
7. Surani MA
(2015) Chromatin dynamics and the role of G9a in gene regulation and enhancer silencing during early mouse development
eLife 4:e09571.
https://doi.org/10.7554/eLife.09571
- PubMed
- Google Scholar

Article and author information

Author details

Jan J Zylicz
1. Wellcome Trust/Cancer Research United Kingdom Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
2. Department of Physiology, Development and Neuroscience, University of Cambridge, Cambridge, United Kingdom
3. Wellcome Trust/Medical Research Council Stem Cell Institute, University of Cambridge, Cambridge, United Kingdom
Contribution
Conceptualization, Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing—original draft

Contributed equally with
Maud Borensztein

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-9622-5658
Maud Borensztein
1. Wellcome Trust/Cancer Research United Kingdom Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
2. Department of Physiology, Development and Neuroscience, University of Cambridge, Cambridge, United Kingdom
Contribution
Validation, Investigation, Writing—review and editing

Contributed equally with
Jan J Zylicz

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-4378-5018
Frederick CK Wong
1. Wellcome Trust/Cancer Research United Kingdom Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
2. Department of Physiology, Development and Neuroscience, University of Cambridge, Cambridge, United Kingdom
Contribution
Validation, Investigation

Competing interests
No competing interests declared
Yun Huang
1. Wellcome Trust/Cancer Research United Kingdom Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
2. Department of Physiology, Development and Neuroscience, University of Cambridge, Cambridge, United Kingdom
Contribution
Resources, Data curation, Formal analysis, Investigation, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-7843-9126
Caroline Lee
1. Wellcome Trust/Cancer Research United Kingdom Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
2. Department of Physiology, Development and Neuroscience, University of Cambridge, Cambridge, United Kingdom
Contribution
Investigation

Competing interests
No competing interests declared
Sabine Dietmann

Wellcome Trust/Medical Research Council Stem Cell Institute, University of Cambridge, Cambridge, United Kingdom

Contribution
Resources, Data curation, Formal analysis, Investigation, Writing—review and editing

Competing interests
No competing interests declared
M Azim Surani
1. Wellcome Trust/Cancer Research United Kingdom Gurdon Institute, University of Cambridge, Cambridge, United Kingdom
2. Department of Physiology, Development and Neuroscience, University of Cambridge, Cambridge, United Kingdom
Contribution
Conceptualization, Supervision, Funding acquisition, Project administration, Writing—review and editing

For correspondence
a.surani@gurdon.cam.ac.uk

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-8640-4318

Funding

Wellcome (096738)

Jan J Zylicz
Maud Borensztein
Yun Huang
Caroline Lee
Sabine Dietmann
M Azim Surani

Wellcome (RG44593)

Jan J Zylicz

H2020 Marie Skłodowska-Curie Actions (706144)

Maud Borensztein

Cancer Research UK (C6946/A14492)

Jan J Zylicz
Maud Borensztein
Yun Huang
Caroline Lee
Sabine Dietmann
M Azim Surani

James Baird Fund, University of Cambridge

Yun Huang

Wellcome (092096)

Jan J Zylicz
Maud Borensztein
Yun Huang
Caroline Lee
Sabine Dietmann
M Azim Surani

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We are grateful to Alexander Tarakhovsky and Dónal O'Carroll for sharing G9a conditional knockout mice. We thank Dang Vinh Do for critical input into the project and members of the Surani Lab for helpful discussions.

Ethics

Animal experimentation: Animal experimentation: All husbandry and experiments involving mice were authorised by a UK Home Office Project Licenses 80/2637 and PE596D1FE and carried out in a Home Office-designated facility.

Version history

Received: November 14, 2017
Accepted: May 9, 2018
Accepted Manuscript published: May 10, 2018 (version 1)
Version of Record published: May 18, 2018 (version 2)

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.