Biological processes are inherently dynamic, yet current methods for capturing temporal changes remain limited. Here, we present scDynaBar, a novel approach that combines CRISPR-Cas9 dynamic barcoding with single-cell sequencing. In this system, genetic barcodes gradually accumulate mutations over time; these barcodes are sequenced alongside the transcriptome of individual cells. We propose that the divergence of these barcodes from the original sequence can serve as a record of the timing of cellular events. To demonstrate the potential of this method, we track the transition from a pluripotent state to a two-cell (2C)-like state in mouse embryonic stem cells (mESCs), providing evidence for the transient nature of the 2C-like state. Additionally, our system shows consistent mutation rates across diverse cell types in a mouse gastruloid model, highlighting its applicability to other biological systems. This approach not only improves our ability to study single-cell dynamics but also opens up new possibilities for recording other temporal signals-in other words, using dynamic barcoding as a molecular clock in individual cells.
The development of complex tissues relies on the precise assignment of cell identity. At the molecular scale, this process depends on the deposition of epigenetic modifications-such as methylation-that are regulated by complex biochemical networks and occur at specific regions on the DNA and chromatin. Here we show that despite the complexity of epigenetic regulation, dynamical scaling and self-similarity of DNA methylation marks emerge in embryonic development. Drawing on single-cell multi-omics experiments, super-resolution microscopy and statistical physics, we demonstrate that these phenomena originate in dynamical feedback between DNA methylation and the formation of nanoscale dynamic chromatin aggregates. These nanoscale processes lead to genome-wide increase in DNA methylation marks following a power law and self-similar correlation functions. Using this framework, we identify methylation patterns that precede gene expression changes in embryonic symmetry breaking. Our work identifies linear sequencing measurements as a laboratory to study mesoscopic biophysical processes in vivo.
In multi-cellular eukaryotic organisms, cell type and specific functional identity are defined by the epigenetic patterning of chemical modifications to DNA and chromatin that modulate the expression and silencing of specific genes. When a cell divides, histones containing important epigenetic marks are distributed between the two daughter strands leading to a temporary dilution of epigenetic information and cell identity. In this work we introduce a physics-based model of epigenetic memory that explains how cells restore and maintain H3K9me3 and H3K27me3 histone methylation patterning after cell division. We demonstrate that emergence and maintenance of the epigenetic program is driven by an evolved mechanism that makes use of the biophysics of polymers, phase condensates and enzymatic activity. We validate our model via genome-wide epigenetic time-course simulation and comparison to experimental epigenetic data from multiple donors, multiple cell types, and for multiple epigenetic marks. Finally, we use our model as a conceptual framework to understand cellular reprogramming by hypothesizing that these processes first contend with and later utilize somatic epigenetic maintenance programs. ### Competing Interest Statement Most authors on this manuscript work at Altos labs, a company devoted to developing reprogramming based therapies.
Aging is a complex multifactorial process that affects cellular function and tissue homeostasis over time. Despite substantial research, the molecular mechanisms driving cellular aging remain poorly understood. Many studies focused on changes in DNA methylation as an indicator of aging. In particular, methylation at polycomb CpG islands was shown to be predictive of phenotypic changes associated with aging. Since many age-related pathological processes are thought to originate from single cells, we asked whether polycomb CpG island methylation occurs preferentially in a subset of cells within a population. Using single-cell whole-genome methylation data across ages and tissues, we identify polycomb CpG methylation as a hallmark of cellular aging. This revealed that aging occurs at varying rates, with faster proliferating cells showing accelerated gain of methylation. Differential gene expression analysis identified changes in immune response, translation, tumorigenesis and neurodegeneration. These results challenge traditional models of homogeneous cellular aging and suggest that aging is a highly individualized process at the single-cell level, that may be driven by programmed changes in polycomb CpG island DNA methylation.
Developmental timing differs strikingly between mammals. All embryonic and some placental lineages emerge from the pluripotent epiblast in a temporally defined sequence, lasting about two weeks in the human embryo, but only two days in mice. Moreover, the order of lineage segregation and gene expression differ between the species. We used human pluripotent stem cells to recapitulate this window of epiblast development in vitro. Simultaneous profiling of gene expression and chromatin accessibility in single cells revealed a robust, autonomous switch between cell states during this process. We reconstructed the integrative gene regulatory network (GRN) of transcription factors (TFs) and signalling molecules. Notably, this revealed a transcriptional cascade including temporally close positive and distant inhibitory connections. We suggest that this cascade acts as a transcriptional clock and governs the directionality, timing and intrinsic decisions during epiblast development. From individual gene interactions, we derived a mechanistic mathematical model of the transcriptional clock of pluripotency that closely reproduced gene expression dynamics and identified key regulatory connections. Moreover, our model revealed a bistable switch governing the transition, thus translating the GRN inference into an interpretable mechanism. Strikingly, the GRN model derived for humans predicted the acceleration of developmental timing in mice when initialised with mouse-specific expression patterns, yet still leading to a human-like expression state. Therefore, TF levels explain developmental timing, whereas the architecture of the network defines its trajectory. Together, our work provides novel conceptual insights into the intrinsic mechanisms of developmental transitions in the human embryo. ### Competing Interest Statement W.R. is an employee of Altos Labs and a consultant and shareholder of Biomodal.
During gastrulation, mouse epiblast cells form the three germ layers that establish the body plan and initiate organogenesis. While single-cell atlases have advanced our understanding of lineage diversification, spatial aspects of differentiation remain poorly defined. Here, we applied spatial transcriptomics to mouse embryos at embryonic (E) E7.25 and E7.5 days and integrated these data with existing E8.5 spatial and E6.5-E9.5 single-cell RNA-seq atlases. This resulted in a spatiotemporal atlas of over 150,000 cells with 82 refined cell-type annotations. The resource enables exploration of gene expression dynamics across anterior-posterior and dorsal-ventral axes, uncovering spatial logic guiding mesodermal fate decisions in the primitive streak. We also developed a computational pipeline to project additional single-cell datasets into this framework for comparative analysis. Freely accessible through an interactive web portal, this atlas offers a valuable tool for the developmental and stem cell biology communities to investigate mouse embryogenesis in a spatial and temporal context.
Embryonic development requires the accurate spatiotemporal execution of cell lineage-specific gene expression programs, which are controlled by transcriptional enhancers. Developmental enhancers adopt a primed chromatin state prior to their activation. How this primed enhancer state is established and maintained and how it affects the regulation of developmental gene networks remains poorly understood. Here, we use comparative multi-omic analyses of human and mouse early embryonic development to identify subsets of postgastrulation lineage-specific enhancers which are epigenetically primed ahead of their activation, marked by the histone modification H3K4me1 within the epiblast. We show that epigenetic priming occurs at lineage-specific enhancers for all three germ layers and that epigenetic priming of enhancers confers lineage-specific regulation of key developmental gene networks. Surprisingly in some cases, lineage-specific enhancers are epigenetically marked already in the zygote, weeks before their activation during lineage specification. Moreover, we outline a generalizable strategy to use naturally occurring human genetic variation to delineate important sequence determinants of primed enhancer function. Our findings identify an evolutionarily conserved program of enhancer priming and begin to dissect the temporal dynamics and mechanisms of its establishment and maintenance during early mammalian development.
Profiling combinations of histone modifications identifies gene regulatory elements in different states and discovers features controlling transcriptional and epigenetic programs. However, efforts to map chromatin states in complex, heterogeneous samples are hindered by the lack of methods that can profile multiple histone modifications together with transcriptomes in individual cells. Here, we describe single-cell multitargets and mRNA sequencing (scMTR-seq), a high-throughput method that enables simultaneous profiling of six histone modifications and transcriptome in single cells. We apply scMTR-seq to uncover dynamic and coordinated changes in chromatin states and transcriptomes during human endoderm differentiation. We also use scMTR-seq to produce lineage-resolved chromatin maps and gene regulatory networks in mouse blastocysts, revealing epigenetic asymmetries at gene regulatory regions between the three embryo lineages and identifying Trps1 as a potential repressor in epiblast cells of trophectoderm-associated enhancer networks and their target genes. Together, scMTR-seq enables investigation of combinatorial chromatin landscapes in a broad range of heterogeneous samples, providing insights into epigenetic regulatory systems.
During mouse gastrulation, extraembryonic mesoderm (ExEM) contributes to the extraembryonic yolk sac (YS) and allantois, both of which are essential for successful gestation. Although the genetic networks coordinating intra-embryonic mesodermal subtype specification are well studied, ExEM diversification remains poorly understood. Here, we identify that embryoid body (EB) in vitro differentiation generates distinct lineages of mesodermal cells, matching YS and allantois development. Combining in vitro and in vivo mouse models, we discover that Eomesodermin (Eomes) controls the formation of YS-fated ExEM but is dispensable for allantois formation. Furthermore, simultaneous disruption of Eomes and T impedes the specification of any YS or allantois mesoderm, indicating compensatory roles for T during allantois formation upon Eomes depletion. Our study highlights previously unrecognized functional and mechanistic diversity in ExEM diversification and endothelial development and introduces a tractable EB model to dissect the signaling pathways and transcriptional networks driving the formation of key extraembryonic tissues.
Mutations in DNA damage response (DDR) factors are associated with human infertility, which affects up to 15% of the population. It remains unclear if the role of DDR is solely in meiosis. One pathway implicated in human fertility is DNA translesion synthesis (TLS), which allows replication impediments to be bypassed. We find that TLS is essential for pre-meiotic germ cell development in the embryo. Loss of the central TLS component, REV1, significantly inhibits the induction of human PGC-like cells (hPGCLCs). This is recapitulated in mice, where deficiencies in TLS initiation ( Rev1 -/- or Pcna K164R/K164R ) or extension ( Rev7 -/- ) result in a >150-fold reduction in the number of primordial germ cells (PGCs) and complete sterility. In contrast, the absence of TLS does not impact the growth, function, or homeostasis of somatic tissues. Surprisingly, we find a complete failure in both activation of the germ cell transcriptional program and in DNA demethylation, a critical step in germline epigenetic reprogramming. Our findings show that for normal fertility, DNA repair is required not only for meiotic recombination but for progression through the earliest stages of germ cell development in mammals.
Early human trophoblast development has remained elusive due to the inaccessibility of the early conceptus. Non-human primate models recapitulate many features of human development and allow access to early postimplantation stages. Here, we tracked the pre- to postimplantation transition of the trophoblast lineage in superficially implanting marmoset embryos in vivo. We differentiated marmoset naive pluripotent stem cells into trophoblast stem cells (TSCs), which exhibited trophoblast-specific transcriptome, methylome, differentiation potential, and long-term self-renewal. Notably, human TSC culture conditions failed to support marmoset TSC derivation, instead inducing an extraembryonic mesoderm-like fate in marmoset cells. We show that combined MEK, TGF-β/NODAL, and histone deacetylase inhibition stabilizes a periimplantation trophoblast-like identity in marmoset TSCs. By contrast, these conditions differentiated human TSCs toward extravillous trophoblasts. Our work presents a paradigm to harness the evolutionary divergence in implantation strategies to elucidate human trophoblast development and invasion.
To implant in the uterus, mammalian embryos form blastocysts comprising trophectoderm (TE) surrounding an inner cell mass (ICM), confined to the polar region by the expanding blastocoel. The mode of implantation varies between species. Murine embryos maintain a single layered TE until they implant in the characteristic thick deciduum, whereas human blastocysts attach via polar TE directly to the uterine wall. Using immunofluorescence (IF) of rapidly isolated ICMs, blockade of RNA and protein synthesis in whole embryos, or 3D visualization of immunostained embryos, we provide evidence of multi-layering in human polar TE before implantation. This may be required for rapid uterine invasion to secure the developing human embryo and initiate formation of the placenta. Using sequential fluorescent labeling, we demonstrate that the majority of inner TE in human blastocysts arises from existing outer cells, with no evidence of conversion from the ICM in the context of the intact embryo.
Ageing is the accumulation of changes and decline of function of organisms over time. The concept and biomarkers of biological age have been established, notably DNA methylation-based clocks. The emergence of single-cell DNA methylation profiling methods opens the possibility of studying the biological age of individual cells. Here, we generate a large single-cell DNA methylation and transcriptome dataset from mouse peripheral blood samples, spanning a broad range of ages. The number of genes expressed increases with age, but gene-specific changes are small. We next develop scEpiAge, a single-cell DNA methylation age predictor, which can accurately predict age in (very sparse) publicly available datasets, and also in single cells. DNA methylation age distribution is wider than technically expected, indicating epigenetic age heterogeneity and functional differences. Our work provides a foundation for single-cell and sparse data epigenetic age predictors, validates their functionality and highlights epigenetic heterogeneity during ageing.
SUMMARYAt the onset of murine gastrulation, pluripotent epiblast cells migrate through the primitive streak, generating mesodermal and endodermal precursors, while the ectoderm arises from the remaining epiblast. Together, these germ layers establish the body plan, defining major body axes and initiating organogenesis. Although comprehensive single cell transcriptional atlases of dissociated mouse embryos across embryonic stages have provided valuable insights during gastrulation, the spatial context for cell differentiation and tissue patterning remain underexplored. In this study, we employed spatial transcriptomics to measure gene expression in mouse embryos at E6.5 and E7.5 and integrated these datasets with previously published E8.5 spatial transcriptomics1and a scRNA-seq2atlas spanning E6.5 to E9.5. This approach resulted in a comprehensive spatiotemporal atlas, comprising over 150,000 cells with 88 refined cell type annotations as well as genome-wide transcriptional imputation during mouse gastrulation and early organogenesis. The atlas facilitates exploration of gene expression dynamics along anterior-posterior and dorsal-ventral axes at cell type, tissue, and organismal scales, revealing insights into mesodermal fate decisions within the primitive streak. Moreover, we developed a bioinformatics pipeline to project additional scRNA-seq datasets into a spatiotemporal framework and demonstrate its utility by analysing cardiovascular models of gastrulation3. To maximise impact, the atlas is publicly accessible via a user-friendly web portal empowering the wider developmental and stem cell biology communities to explore mechanisms of early mouse development in a spatiotemporal context.
Ageing is the decline of physiological function over time. Statistical models termed ageing clocks can predict chronological and biological age from the longitudinal time evolution of DNA methylation. Here, we show that DNA methylation ageing is also manifest in how pairs of genomic loci evolve with respect to each other with age. Using sequencing data from a range of tissues in mouse we show that genomic correlations in DNA methylation during ageing are characterised by an enrichment of correlations between sites that are roughly 500 bp apart. We trace the origin of this behaviour to collective dynamics in the boundaries of CpG islands. We derive a simple, biophysical model that explains the origin due to a tilt in the competition between methylating and demethylating processes and an ensuing wetting-like phenomenon. Using Chip-seq data, we argue for a molecular mechanism based on a PRC2-dependent dilution of H3K27me3 during ageing. Our work gives a new perspective on epigenetic ageing, highlighting the importance of correlative as opposed to longitudinal dynamics. ### Competing Interest Statement W.R. is a consultant and shareholder of Cambridge Epigenetix. W.R. is an employee of Altos Labs. O.S. is a paid consultant of Insitro.INC. F.v.M. is a consultant and shareholder of Longevity Consultancy Group Sarl. The remaining authors declare no competing interests.
Human pluripotent stem cells (hPSCs) are of fundamental relevance in regenerative medicine. Naïve hPSCs hold promise to overcome some of the limitations of conventional (primed) hPSCs, including recurrent epigenetic anomalies. Naïve-to-primed transition (capacitation) follows transcriptional dynamics of human embryonic epiblast and is necessary for somatic differentiation from naïve hPSCs. We found that capacitated hPSCs are transcriptionally closer to postimplantation epiblast than conventional hPSCs. This prompted us to comprehensively study epigenetic and related transcriptional changes during capacitation. Our results show that CpG islands, gene regulatory elements, and retrotransposons are hotspots of epigenetic dynamics during capacitation and indicate possible distinct roles of specific epigenetic modifications in gene expression control between naïve and primed hPSCs. Unexpectedly, PRC2 activity appeared to be dispensable for the capacitation. We find that capacitated hPSCs acquire an epigenetic state similar to conventional hPSCs. Significantly, however, the X chromosome erosion frequently observed in conventional female hPSCs is reversed by resetting and subsequent capacitation.
Zygotic genome activation (ZGA) in mice takes place in two waves, a minor wave in the one-cell embryo, and a major wave at the two-cell stage, both accompanied by global transcriptional and epigenetic reprogramming. However, the orchestration of these reprogramming events by maternal factors deposited in the oocyte is not yet entirely understood. We and others have recently shown that epigenetic modifiers such as SMARCA5 (the main ATPase in ISWI complexes) can initiate the ZGA transcriptional programme in vitro . So far, the role of SMARCA5 in ZGA in vivo has not been addressed, as constitutive knock-out mice lacking SMARCA5 are not viable. We have overcome this limitation by using the targeted protein-depletion system Trim Away to degrade SMARCA5 in early zygotes. We further harnessed the power of single cell multi-omics (scNMT-seq) and showed that in the absence of SMARCA5, major ZGA genes fail to be upregulated at the two-cell stage. This is explained by the lower accessibility and disrupted nucleosome positioning at their promoters and distal regulatory regions, compared to wild-type embryos. In contrast, we show that global chromatin accessibility at the two-cell stage is higher in SMARCA5 depleted embryos compared to control embryos, and this is accompanied by other global structural changes involving heterochromatic regions. Our results show that SMARCA5 has a global regulatory role at the two-cell stage, which includes the control of ZGA gene promoters and distal regulatory regions.### Competing Interest StatementThe authors have declared no competing interest.
To enable implantation in the uterus, mammalian embryos form blastocysts comprising trophectoderm surrounding an inner cell mass, confined to the polar region by the expanding blastocyst cavity. The mode of implantation varies considerably between species. Murine embryos maintain a single layered trophectoderm until sometime after they have implanted in the characteristic thick decidual tissue, whereas human blastocysts attach directly to the uterine wall. Using immunofluorescence of rapidly isolated inner cell masses and blockade of RNA and protein synthesis in whole embryos, we provide evidence of multi-layering in human polar trophectoderm before implantation. This may be a prerequisite for rapid invasion required to secure the developing human embryo and initiate formation of the placenta. Using sequential fluorescent labelling, we demonstrate that the majority of inner trophectoderm in human blastocysts arises from existing outer cells with no evidence of conversion from the inner cell mass in the context of the intact embryo.