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.

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.

The gill skeleton of cartilaginous fishes (sharks, skates, rays, and holocephalans) exhibits a striking anterior–posterior polarity, with a series of fine appendages called branchial rays projecting from the posterior margin of the gill arch cartilages. We previously demonstrated in the skate (Leucoraja erinacea) that branchial rays derive from a posterior domain of pharyngeal arch mesenchyme that is responsive to Sonic hedgehog (Shh) signaling from a distal gill arch epithelial ridge (GAER) signaling centre. However, how branchial ray progenitors are specified exclusively within posterior gill arch mesenchyme is not known. Here, we show that genes encoding several Wnt ligands are expressed in the ectoderm immediately adjacent to the skate GAER, and that these Wnt signals are transduced largely in the anterior arch environment. Using pharmacological manipulation, we show that inhibition of Wnt signalling results in an anterior expansion of Shh signal transduction in developing skate gill arches, and in the formation of ectopic anterior branchial ray cartilages. Our findings demonstrate that ectodermal Wnt signalling contributes to gill arch skeletal polarity in skate by restricting Shh signal transduction and chondrogenesis to the posterior arch environment and highlights the importance of signalling interactions at embryonic tissue boundaries for cell fate determination in vertebrate pharyngeal arches.

ABSTRACT Ageing is the accumulation of changes and overall decline of the function of cells, organs and organisms over time. At the molecular and cellular level, the concept of biological age has been established and biomarkers of biological age have been identified, notably epigenetic DNA-methylation based clocks. With the emergence of single-cell DNA methylation profiling methods, the possibility to study biological age of individual cells has been proposed, and a first proof-of-concept study, based on limited single cell datasets mostly from early developmental origin, indicated the feasibility and relevance of this approach to better understand organismal changes and cellular ageing heterogeneity. Here we generated a large single-cell DNA methylation and matched transcriptome dataset from mouse peripheral blood samples, spanning a broad range of ages (10-101 weeks of age). We observed that the number of genes expressed increased at older ages, but gene specific changes were small. We next developed a robust single cell DNA methylation age predictor (scEpiAge), which can accurately predict age in a broad range of publicly available datasets, including very sparse data and it also predicts age in single cells. Interestingly, the DNA methylation age distribution is wider than technically expected in 19% of single cells, suggesting that epigenetic age heterogeneity is present in vivo and may relate to functional differences between cells. In addition, we observe differences in epigenetic ageing between the major blood cell types. Our work provides a foundation for better single-cell and sparse data epigenetic age predictors and highlights the significance of cellular heterogeneity during ageing. Highlights - Model to estimate DNA methylation age in single cells - Large multi-omics dataset of single cells from murine blood - Epigenetic age deviations from chronological age are greater than technical expected from technical variability - Number of genes expressed increases with chronological and epigenetic age

Common variable immunodeficiency (CVID), the most prevalent symptomatic primary immunodeficiency, displays impaired terminal B-cell differentiation and defective antibody responses. Incomplete genetic penetrance and ample phenotypic expressivity in CVID suggest the participation of additional pathogenic mechanisms. Monozygotic (MZ) twins discordant for CVID are uniquely valuable for studying the contribution of epigenetics to the disease. Here, we generate a single-cell epigenomics and transcriptomics census of naïve-to-memory B cell differentiation in a CVID-discordant MZ twin pair. Our analysis identifies DNA methylation, chromatin accessibility and transcriptional defects in memory B-cells mirroring defective cell-cell communication upon activation. These findings are validated in a cohort of CVID patients and healthy donors. Our findings provide a comprehensive multi-omics map of alterations in naïve-to-memory B-cell transition in CVID and indicate links between the epigenome and immune cell cross-talk. Our resource, publicly available at the Human Cell Atlas, gives insight into future diagnosis and treatments of CVID patients.

BACKGROUND:Perturbation of DNA methyltransferases (DNMTs) and of the active DNA demethylation pathway via ten-eleven translocation (TET) methylcytosine dioxygenases results in severe developmental defects and embryonic lethality. Dynamic control of DNA methylation is therefore vital for embryogenesis, yet the underlying mechanisms remain poorly understood. RESULTS:Here we report a single-cell transcriptomic atlas from Dnmt and Tet mutant mouse embryos during early organogenesis. We show that both the maintenance and de novo methyltransferase enzymes are dispensable for the formation of all major cell types at E8.5. However, DNA methyltransferases are required for silencing of prior or alternative cell fates such as pluripotency and extraembryonic programmes. Deletion of all three TET enzymes produces substantial lineage biases, in particular, a failure to generate primitive erythrocytes. Single-cell multi-omics profiling moreover reveals that this is linked to a failure to demethylate distal regulatory elements in Tet triple-knockout embryos. CONCLUSIONS:This study provides a detailed analysis of the effects of perturbing DNA methylation on mouse organogenesis at a whole organism scale and affords new insights into the regulatory mechanisms of cell fate decisions.

Following gastrulation, the three primary germ layers develop into the major organs in a process known as organogenesis. Single-cell RNA sequencing has enabled the profiling of the gene expression dynamics of these cell fate decisions, yet a comprehensive map of the interplay between transcription factors and cis-regulatory elements is lacking, as are the underlying gene regulatory networks. Here we generate a multi-omics atlas of mouse early organogenesis by simultaneously profiling gene expression and chromatin accessibility from tens of thousands of single cells. We develop a computational method to leverage the multimodal readouts to predict transcription factor binding events in cis-regulatory elements, which we then use to infer gene regulatory networks that underpin lineage commitment events. Finally, we show that these models can be used to generate in silico predictions of the effect of transcription factor perturbations. We validate this experimentally by showing that Brachyury is essential for the differentiation of neuromesodermal progenitors to somitic mesoderm fate by priming cis-regulatory elements. The data set can be interactively explored at https://www.bioinformatics.babraham.ac.uk/shiny/shiny_multiome_organogenesis/

Recent breakthroughs in single-cell genomics allow probing molecular states of cells with unprecedented detail along the sequence of the DNA. Biological function relies, however, on emergent processes in the three-dimensional space of the nucleus, such as droplet formation through phase separation. Here, we use single-cell multi-omics sequencing to develop a theoretical framework to rigorously map epigenome profiling along the DNA sequence onto a description of the emergent spatial dynamics in the nucleus. Drawing on scNMT-seq multi-omics sequencing in vitro and in vivo we exemplify our approach in the context of exit from pluripotency and global de novo methylation of the genome. We show how DNA methylation patterns of the embryonic genome are established through the interplay between spatially correlated DNA methylation and topological changes to the DNA. This feedback leads to the predicted formation of 30-40nm sized condensates of methylated DNA and determines genome-scale DNA methylation rates. We verify these findings with orthogonal single cell multi-omics data that combine the methylome with HiC measurements. Notably, this scale of chromatin organization has recently been described by super-resolution microscopy. Using this framework, we identify local methylation correlations in gene bodies that precede transcriptional changes at the exit from pluripotency. Our work provides a general framework of how mechanistic insights into emergent processes underlying cell fate decisions can be gained by the combination of single-cell multi-omics and methods from theoretical physics that have not been applied in the context of genomics before. Highlights We develop methodology to infer collective spatio-temporal processes in the physical space of the nucleus from single-cell methylome sequencing experiments. We show that DNA methylation relies on a feedback between de novo methylation and nanoscale changes in DNA topology, leading to the formation of methylation condensates. Chromatin condensates at this scale have recently been described by high-resolution microscopy but have remained without mechanistic explanation. Using this framework, we identify changes in the distribution of DNA methylation marks in gene bodies that precede gene silencing at the exit from pluripotency.

SummaryRecent breakthroughs in single-cell genomics allow probing molecular states of cells with unprecedented detail along the sequence of the DNA. Biological function relies, however, on emergent processes in the three-dimensional space of the nucleus, such as droplet formation through phase separation. Here, we use single-cell multi-omics sequencing to develop a theoretical framework to rigorously map epigenome profiling along the DNA sequence onto a description of the emergent spatial dynamics in the nucleus. Drawing on scNMT-seq multi-omics sequencingin vitroandin vivowe exemplify our approach in the context of exit from pluripotency and globalde novomethylation of the genome. We show how DNA methylation patterns of the embryonic genome are established through the interplay between spatially correlated DNA methylation and topological changes to the DNA. This feedback leads to the predicted formation of 30-40nm sized condensates of methylated DNA and determines genome-scale DNA methylation rates. We verify these findings with orthogonal single cell multi-omics data that combine the methylome with HiC measurements. Notably, this scale of chromatin organization has recently been described by super-resolution microscopy. Using this framework, we identify local methylation correlations in gene bodies that precede transcriptional changes at the exit from pluripotency. Our work provides a general framework of how mechanistic insights into emergent processes underlying cell fate decisions can be gained by the combination of single-cell multi-omics and methods from theoretical physics that have not been applied in the context of genomics before.HighlightsWe develop methodology to infer collective spatio-temporal processes in the physical space of the nucleus from single-cell methylome sequencing experiments.We show that DNA methylation relies on a feedback betweende novomethylation and nanoscale changes in DNA topology, leading to the formation of methylation condensates.Chromatin condensates at this scale have recently been described by high-resolution microscopy but have remained without mechanistic explanation.Using this framework, we identify changes in the distribution of DNA methylation marks in gene bodies that precede gene silencing at the exit from pluripotency.

The origin of the jaw is a long-standing problem in vertebrate evolutionary biology. Classical hypotheses of serial homology propose that the upper and lower jaw evolved through modifications of dorsal and ventral gill arch skeletal elements, respectively. If the jaw and gill arches are derived members of a primitive branchial series, we predict that they would share common developmental patterning mechanisms. Using candidate and RNAseq/differential gene expression analyses, we find broad conservation of dorsoventral (DV) patterning mechanisms within the developing mandibular, hyoid, and gill arches of a cartilaginous fish, the skate (Leucoraja erinacea). Shared features include expression of genes encoding members of the ventralizing BMP and endothelin signaling pathways and their effectors, the joint markers nkx3.2 and gdf5 and prochondrogenic transcription factor barx1, and the dorsal territory marker pou3f3. Additionally, we find that mesenchymal expression of eya1/six1 is an ancestral feature of the mandibular arch of jawed vertebrates, whereas differences in notch signaling distinguish the mandibular and gill arches in skate. Comparative transcriptomic analyses of mandibular and gill arch tissues reveal additional genes differentially expressed along the DV axis of the pharyngeal arches, including scamp5 as a novel marker of the dorsal mandibular arch, as well as distinct transcriptional features of mandibular and gill arch muscle progenitors and developing gill buds. Taken together, our findings reveal conserved patterning mechanisms in the pharyngeal arches of jawed vertebrates, consistent with serial homology of their skeletal derivatives, as well as unique transcriptional features that may underpin distinct jaw and gill arch morphologies.

Germline specification in mammals occurs through an inductive process whereby competent cells in the post-implantation epiblast differentiate into primordial germ cells (PGC). The intrinsic factors that endow epiblast cells with the competence to respond to germline inductive signals remain unknown. Single-cell RNA sequencing across multiple stages of an in vitro PGC-like cells (PGCLC) differentiation system shows that PGCLC genes initially expressed in the naïve pluripotent stage become homogeneously dismantled in germline competent epiblast like-cells (EpiLC). In contrast, the decommissioning of enhancers associated with these germline genes is incomplete. Namely, a subset of these enhancers partly retain H3K4me1, accumulate less heterochromatic marks and remain accessible and responsive to transcriptional activators. Subsequently, as in vitro germline competence is lost, these enhancers get further decommissioned and lose their responsiveness to transcriptional activators. Importantly, using H3K4me1-deficient cells, we show that the loss of this histone modification reduces the germline competence of EpiLC and decreases PGCLC differentiation efficiency. Our work suggests that, although H3K4me1 might not be essential for enhancer function, it can facilitate the (re)activation of enhancers and the establishment of gene expression programs during specific developmental transitions.

Advancing maternal age causes a progressive reduction in fertility. The decline in developmental competence of the oocyte with age is likely to be a consequence of multiple contributory factors. Loss of epigenetic quality of the oocyte could impair early developmental events or programme adverse outcomes in offspring that manifest only later in life. Here, we undertake joint profiling of the transcriptome and DNA methylome of individual oocytes from reproductively young and old mice undergoing natural ovulation. We find reduced complexity as well as increased variance in the transcriptome of oocytes from aged females. This transcriptome heterogeneity is reflected in the identification of discrete sub-populations. Oocytes with a transcriptome characteristic of immature chromatin configuration (NSN) clustered into two groups: one with reduced developmental competence, as indicated by lower expression of maternal effect genes, and one with a young-like transcriptome. Oocytes from older females had on average reduced CpG methylation, but the characteristic bimodal methylation landscape of the oocyte was preserved. Germline differentially methylated regions of imprinted genes were appropriately methylated irrespective of age. For the majority of differentially expressed transcripts, the absence of correlated methylation changes suggests a post-transcriptional basis for most age-related effects on the transcriptome. However, we did find differences in gene body methylation at which there were corresponding changes in gene expression, indicating age-related effects on transcription that translate into methylation differences. Interestingly, oocytes varied in expression and methylation of these genes, which could contribute to variable competence of oocytes or penetrance of maternal age-related phenotypes in offspring.

Background Alternative splicing is a key regulatory mechanism in eukaryotic cells and increases the effective number of functionally distinct gene products. Using bulk RNA sequencing, splicing variation has been studied across human tissues and in genetically diverse populations. This has identified disease-relevant splicing events, as well as associations between splicing and genomic features, including sequence composition and conservation. However, variability in splicing between single cells from the same tissue or cell type and its determinants remains poorly understood. Results We applied parallel DNA methylation and transcriptome sequencing to differentiating human induced pluripotent stem cells to characterize splicing variation (exon skipping) and its determinants. Our results show that variation in single-cell splicing can be accurately predicted based on local sequence composition and genomic features. We observe moderate but consistent contributions from local DNA methylation profiles to splicing variation across cells. A combined model that is built based on genomic features as well as DNA methylation information accurately predicts different splicing modes of individual cassette exons. These categories include the conventional inclusion and exclusion patterns, but also more subtle modes of cell-to-cell variation in splicing. Finally, we identified and characterized associations between DNA methylation and splicing changes during cell differentiation. Conclusions Our study yields new insights into alternative splicing at the single-cell level and reveals a previously underappreciated link between DNA methylation variation and splicing.

Bivalent chromatin marks developmental promoters in pluripotent cells, yet their targeting and precise impact on lineage commitment remains unclear. We uncover Developmental Pluripotency Associated 2 (Dppa2) and 4 (Dppa4) as epigenetic priming factors, establishing chromatin bivalency. Single-cell transcriptomics and differentiation assays reveal Dppa2/4 double knockout embryonic stem cells fail to exit pluripotency and differentiate efficiently. Dppa2/4 associate with COMPASS and Polycomb complexes and are required to recruit and maintain their binding at a subset of developmentally important bivalent promoters which are characterised by low expression and poised RNA polymerase. Consequently, upon Dppa2/4 knockout, these dependent promoters gain DNA methylation and are unable to be activated upon differentiation. Our findings uncover a novel targeting principle for bivalency to developmental promoters, poising them for future lineage specific activation.

Formation of the three primary germ layers during gastrulation is an essential step in the establishment of the vertebrate body plan. Recent studies employing single cell RNA-sequencing have identified major transcriptional changes associated with germ layer specification. Global epigenetic reprogramming accompanies these changes, but the role of the epigenome in regulating early cell fate choice remains unresolved, and the coordination between different epigenetic layers is unclear. Here we describe the first single cell triple-omics map of chromatin accessibility, DNA methylation and RNA expression during the exit from pluripotency and the onset of gastrulation in mouse embryos. We find dynamic dependencies between the different molecular layers, with evidence for distinct modes of epigenetic regulation. The initial exit from pluripotency coincides with the establishment of a global repressive epigenetic landscape, followed by the emergence of local lineage-specific epigenetic patterns during gastrulation. Notably, cells committed to mesoderm and endoderm undergo widespread coordinated epigenetic rearrangements, driven by loss of methylation in enhancer marks and a concomitant increase of chromatin accessibility. In striking contrast, the epigenetic landscape of ectodermal cells is already established in the early epiblast. Hence, regulatory elements associated with each germ layer are either epigenetically primed or epigenetically remodelled prior to overt cell fate decisions during gastrulation, providing the molecular logic for a hierarchical emergence of the primary germ layers. Highlights First map of mouse gastrulation using comprehensive single cell triple-omic analysis. Exit from pluripotency is associated with a global repressive epigenetic landscape, driven by a sharp gain of DNA methylation and a gradual decrease of chromatin accessibility. DNA methylation and chromatin accessibility changes in enhancers, but not in promoters, are associated with germ layer formation. Mesoderm and endoderm enhancers become open and demethylated upon lineage commitment. Ectoderm enhancers are primed in the early epiblast and protected from the global repressive dynamics, supporting a default model of ectoderm commitment in vivo .

Formation of the three primary germ layers during gastrulation is an essential step in the establishment of the vertebrate body plan and is associated with major transcriptional changes1-5. Global epigenetic reprogramming accompanies these changes6-8, but the role of the epigenome in regulating early cell-fate choice remains unresolved, and the coordination between different molecular layers is unclear. Here we describe a single-cell multi-omics map of chromatin accessibility, DNA methylation and RNA expression during the onset of gastrulation in mouse embryos. The initial exit from pluripotency coincides with the establishment of a global repressive epigenetic landscape, followed by the emergence of lineage-specific epigenetic patterns during gastrulation. Notably, cells committed to mesoderm and endoderm undergo widespread coordinated epigenetic rearrangements at enhancer marks, driven by ten-eleven translocation (TET)-mediated demethylation and a concomitant increase of accessibility. By contrast, the methylation and accessibility landscape of ectodermal cells is already established in the early epiblast. Hence, regulatory elements associated with each germ layer are either epigenetically primed or remodelled before cell-fate decisions, providing the molecular framework for a hierarchical emergence of the primary germ layers.

Pluripotency is accompanied by the erasure of parental epigenetic memory, with naïve pluripotent cells exhibiting global DNA hypomethylation both in vitro and in vivo. Exit from pluripotency and priming for differentiation into somatic lineages is associated with genome-wide de novo DNA methylation. We show that during this phase, co-expression of enzymes required for DNA methylation turnover, DNMT3s and TETs, promotes cell-to-cell variability in this epigenetic mark. Using a combination of single-cell sequencing and quantitative biophysical modeling, we show that this variability is associated with coherent, genome-scale oscillations in DNA methylation with an amplitude dependent on CpG density. Analysis of parallel single-cell transcriptional and epigenetic profiling provides evidence for oscillatory dynamics both in vitro and in vivo. These observations provide insights into the emergence of epigenetic heterogeneity during early embryo development, indicating that dynamic changes in DNA methylation might influence early cell fate decisions.