Understanding how gene regulatory networks control the progressive restriction of cell fates is a long-standing challenge. Recent advances in measuring single cell gene expression are providing new insights into lineage commitment. However, the regulatory events underlying these changes remain elusive. Here we investigate the dynamics of chromatin regulatory landscapes during embryogenesis at single cell resolution. Using single cell combinatorial indexing assay for transposase accessible chromatin (sci-ATAC-seq)1, we profiled chromatin accessibility in over 20,000 single nuclei from fixed Drosophila embryos spanning three landmark embryonic stages: 2-4 hours (hrs) after egg laying (predominantly stage 5 blastoderm nuclei), when each embryo comprises ~6,000 multipotent cells; 6-8hrs (predominantly stage 10-11), to capture a midpoint in embryonic development when major lineages in the mesoderm and ectoderm are specified; and 10-12hrs (predominantly stage 13), when each of the embryo’s >20,000 cells are undergoing terminal differentiation. Our results reveal spatial heterogeneity in the usage of the regulatory genome prior to gastrulation, a feature that aligns with future cell fate, and nuclei can be temporally ordered along developmental trajectories. During mid-embryogenesis, tissue granularity emerges such that individual cell types can be inferred by their chromatin accessibility, while maintaining a signature of their germ layer of origin. The data reveal overlapping usage of regulatory elements between cells of the endoderm and non-myogenic mesoderm, suggesting a common developmental program reminiscent of the mesendoderm lineage in other species2–4. Altogether, we identify over 30,000 distal regulatory elements exhibiting tissue-specific accessibility. We validated the germ layer specificity of a subset of these predicted enhancers in transgenic embryos, achieving 90% accuracy. Overall, our results demonstrate the power of shotgun single cell profiling of embryos to resolve dynamic changes in the chromatin landscape during development, and to uncover the cis-regulatory programs of metazoan germ layers and cell types.

The discovery of substantial amounts of 5-hydroxymethylcytosine (5hmC), formed by the oxidation of 5-methylcytosine (5mC), in various mouse tissues and human embryonic stem (ES) cells has necessitated a reevaluation of our knowledge of 5mC/5hmC patterns and functions in mammalian cells. Here, we investigate the tissue specificity of both the global levels and locus-specific distribution of 5hmC in several human tissues and cell lines. We find that global 5hmC content of normal human tissues is highly variable, does not correlate with global 5mC content, and decreases rapidly as cells from normal tissue adapt to cell culture. Using tiling microarrays to map 5hmC levels in DNA from normal human tissues, we find that 5hmC patterns are tissue specific; unsupervised hierarchical clustering based solely on 5hmC patterns groups independent biological samples by tissue type. Moreover, in agreement with previous studies, we find 5hmC associated primarily, but not exclusively, with the body of transcribed genes, and that within these genes 5hmC levels are positively correlated with transcription levels. However, using quantitative 5hmC-qPCR, we find that the absolute levels of 5hmC for any given gene are primarily determined by tissue type, gene expression having a secondary influence on 5hmC levels. That is, a gene transcribed at a similar level in several different tissues may have vastly different levels of 5hmC (>20-fold) dependent on tissue type. Our findings highlight tissue type as a major modifier of 5hmC levels in expressed genes and emphasize the importance of using quantitative analyses in the study of 5hmC levels.

BackgroundDNA methylation and the Polycomb repression system are epigenetic mechanisms that play important roles in maintaining transcriptional repression. Recent evidence suggests that DNA methylation can attenuate the binding of Polycomb protein components to chromatin and thus plays a role in determining their genomic targeting. However, whether this role of DNA methylation is important in the context of transcriptional regulation is unclear.ResultsBy genome-wide mapping of the Polycomb Repressive Complex 2-signature histone mark, H3K27me3, in severely DNA hypomethylated mouse somatic cells, we show that hypomethylation leads to widespread H3K27me3 redistribution, in a manner that reflects the local DNA methylation status in wild-type cells. Unexpectedly, we observe striking loss of H3K27me3 and Polycomb Repressive Complex 2 from Polycomb target gene promoters in DNA hypomethylated cells, including Hox gene clusters. Importantly, we show that many of these genes become ectopically expressed in DNA hypomethylated cells, consistent with loss of Polycomb-mediated repression.ConclusionsAn intact DNA methylome is required for appropriate Polycomb-mediated gene repression by constraining Polycomb Repressive Complex 2 targeting. These observations identify a previously unappreciated role for DNA methylation in gene regulation and therefore influence our understanding of how this epigenetic mechanism contributes to normal development and disease.

SUMMARYMouse primordial germ cells (PGCs) erase global DNA methylation (5mC) as part of the comprehensive epigenetic reprogramming that occurs during PGC development. 5mC plays an important role in maintaining stable gene silencing and repression of transposable elements (TE) but it is not clear how the extensive loss of DNA methylation impacts on gene expression and TE repression in developing PGCs. Using a novel epigenetic disruption and recovery screen and genetic analyses, we identified a core set of germline-specific genes that are dependent exclusively on promoter DNA methylation for initiation and maintenance of developmental silencing. These gene promoters appear to possess a specialised chromatin environment that does not acquire any of the repressive H3K27me3, H3K9me2, H3K9me3 or H4K20me3 histone modifications when silenced by DNA methylation. Intriguingly, this methylation-dependent subset is highly enriched in genes with roles in suppressing TE activity in germ cells. We show that the mechanism for developmental regulation of the germline genome-defence genes involves DNMT3B-dependent de novo DNA methylation. These genes are then activated by lineage-specific promoter demethylation during distinct global epigenetic reprogramming events in migratory (~E8.5) and post-migratory (E10.5-11.5) PGCs. We propose that genes involved in genome defence are developmentally regulated primarily by promoter DNA methylation as a sensory mechanism that is coupled to the potential for TE activation during global 5mC erasure, thereby acting as a failsafe to ensure TE suppression and maintain genomic integrity in the germline.