SummaryDNA methylation is highly dynamic during mammalian embryogenesis. It is broadly accepted that the paternal genome is actively depleted of 5-methyl cytosine at fertilization, followed by passive loss that reaches a minimum at the blastocyst stage. However, this model is based on limited data, and to date no base-resolution maps exist to support and refine it. Here, we generated genome-scale DNA methylation maps in mouse gametes and through post-implantation embryogenesis. We find that the oocyte already exhibits global hypomethylation, most prominently at specific families of long interspersed element-1 and long terminal repeat retro-elements, which are disparate between gametes and resolve to lower methylation values in zygote. Surprisingly, the oocyte contributes a unique set of Differentially Methylated Regions (DMRs), including many CpG Island promoter regions, that are maintained in the early embryo but are lost upon specification and absent from somatic cells. In contrast, sperm-contributed DMRs are largely intergenic and resolve to hypermethylation after the blastocyst stage. Our data provide a complete genome-scale, base-resolution timeline of DNA methylation in the pre-specified embryo, when this epigenetic modification is most dynamic, before returning to the canonical somatic pattern.

In mammals, cytosine methylation is predominantly restricted to CpG dinucleotides and stably distributed across the genome, with local, cell type-specific regulation directed by DNA binding factors1-3. This comparatively static landscape dramatically contrasts the events of fertilization, where the paternal genome is globally reprogrammed. Paternal genome demethylation includes the majority of CpGs, though methylation is maintained at several notable features4-7. While these dynamics have been extensively characterized in the mouse, only limited observations are available in other mammals, and direct measurements are required to understand the extent to which early embryonic landscapes are conserved8-10. We present genome-scale DNA methylation maps of human preimplantation development and embryonic stem cell (ESC) derivation, confirming a transient state of global hypomethylation that includes most CpGs, while sites of persistent maintenance are primarily restricted to gene bodies. While most features share similar dynamics to mouse, maternally contributed methylation is divergently targeted to species-specific sets of CpG island (CGI) promoters that extend beyond known Imprint Control Regions (ICRs). Retrotransposon regulation is also highly diverse and transitions from maternally to embryonically expressed, species-specific elements. Together, our data confirm that paternal genome demethylation is a general attribute of early mammalian development that is characterized by distinct modes of epigenetic regulation.

A transcribed, multi-channel, and continuously evolving molecular recorder 77 To achieve our goal of a tunable, high information content molecular recorder, we 78 utilized Cas9 to generate insertions or deletions (indels) upon repair of double-stranded breaks, 79 which are inherited in the next generation of cells 11-16. We record within a 205 base pair, synthetic DNA "target site" containing three "cut sites" and a static 8 base pair "integration barcode" (intBC), which are delivered in multiple copies via piggyBac transposition (Fig. 1a, b). We embedded this sequence into the 3'UTR of a constitutively transcribed fluorescent protein to enable profiling from the transcriptome. A second cassette encodes three independently transcribed and complementary guide RNAs to permit recording of multiple, distinct signals (Fig. 1a, b) 18. Our system is capable of high information storage due to the diversity of heritable repair outcomes, and the large number of targeted sites, which can be distinguished by the intBC (Fig. 1c). DNA repair generates hundreds of unique indels, and the distribution for each cut site is different and nonuniform: some produce highly biased outcomes while others create a diverse series (Fig. 1c, Extended Data Fig. 1) 19-21. To identify sequences that can tune the mutation rate of our recorder for timescales that are not pre-defined, and may extend from days to months, we screened several guide RNA series containing mismatches to their targets 22 by monitoring their activity on a GFP reporter over a 20-day timecourse and selected those that demonstrated a broad dynamic range (Fig. 1d). Slower cutting rates may improve viability in vivo, as frequent Cas9mediated double-strand breaks can cause cellular toxicity 23,24. To demonstrate information recovery from single cell transcriptomes, we stably transduced K562 cells with our technology and generated a primary, cell-barcoded cDNA pool via the 10x Genomics platform, allowing us * * *