Single-Cell Landscape of Transcriptional Heterogeneity and Cell Fate Decisions during Mouse Early Gastrulation

EpiCHAOS: a metric to quantify epigenomic heterogeneity in single-cell data

Epigenetic heterogeneity is a fundamental property of biological systems and is recognized as a potential driver of tumor plasticity and therapy resistance. Single-cell epigenomics technologies have been widely employed to study epigenetic variation between-but not within-cellular clusters. We introduce epiCHAOS: a quantitative metric of cell-to-cell heterogeneity, applicable to any single-cell epigenomics data type. After validation in synthetic datasets, we apply epiCHAOS to investigate global and region-specific patterns of epigenetic heterogeneity across diverse biological systems. EpiCHAOS provides an excellent approximation of stemness and plasticity in development and malignancy, making it a valuable addition to single-cell cancer epigenomics analyses.

Genetic and epigenetic alterations in aging and rejuvenation of human

All the information essential for life is encoded within our genome and epigenome, which orchestrates diverse cellular states spatially and temporally. In particular, the epigenome interacts with internal and external stimuli, encoding and preserving cellular experiences, and it serves as the regulatory base of the transcriptome across diverse cell types. The emergence of single-cell transcriptomic and epigenomic data collection has revealed unique omics signatures in diverse tissues, highlighting cellular heterogeneity. Recent research has documented age-related epigenetic changes at the single-cell level, alongside the validation of cellular rejuvenation through partial reprogramming, which involves simultaneous epigenetic modifications. These dynamic shifts, primarily fueled by stem cell plasticity, have catalyzed significant interest and cross-disciplinary research endeavors. This review explores the genomic and epigenomic alterations with aging, elucidating their reciprocal interactions. Additionally, it seeks to discuss the evolving landscape of rejuvenation research, with a particular emphasis on dissecting stem cell behavior through the lens of single-cell analysis. Moreover, it proposes potential research methodologies for future studies.

The transcription factor GABPA is a master regulator of naïve pluripotency

The establishment of naïve pluripotency is a continuous process starting with the generation of inner cell mass (ICM) which then differentiating into epiblast (EPI). Recent studies have revealed key transcription factors (TFs) for ICM formation, but which TFs initiate EPI specification remains unknown. Here, using a targeted rapid protein degradation system, we show that GABPA is not only a regulator of major ZGA, but also a master EPI specifier required for naïve pluripotency establishment by regulating 47% of EPI genes during E3.5 to E4.5 transition. Chromatin binding dynamics analysis suggests that GABPA controls EPI formation at least partly by binding to the ICM gene promoters occupied by the pluripotency regulators TFAP2C and SOX2 at E3.5 to establish naïve pluripotency at E4.5. Our study not only uncovers GABPA as a master pluripotency regulator, but also supports the notion that mammalian pluripotency establishment requires a dynamic and stepwise multi-TFs regulatory network.

Deep learning-based models for preimplantation mouse and human embryos based on single-cell RNA sequencing

The rapid growth of single-cell transcriptomic technology has produced an increasing number of datasets for both embryonic development and in vitro pluripotent stem cell-derived models. This avalanche of data surrounding pluripotency and the process of lineage specification has meant it has become increasingly difficult to define specific cell types or states in vivo, and compare these with in vitro differentiation. Here we utilize a set of deep learning tools to integrate and classify multiple datasets. This allows the definition of both mouse and human embryo cell types, lineages and states, thereby maximizing the information one can garner from these precious experimental resources. Our approaches are built on recent initiatives for large-scale human organ atlases, but here we focus on material that is difficult to obtain and process, spanning early mouse and human development. Using publicly available data for these stages, we test different deep learning approaches and develop a model to classify cell types in an unbiased fashion at the same time as defining the set of genes used by the model to identify lineages, cell types and states. We used our models trained on in vivo development to classify pluripotent stem cell models for both mouse and human development, showcasing the importance of this resource as a dynamic reference for early embryogenesis.

A Single Cell Transcriptomic Fingerprint of Stressed Premature, Imbalanced Differentiation of Embryonic Stem Cells

BACKGROUND

Miscarriages cause a greater loss-of-life than cardiovascular diseases, but knowledge about environmentally induced miscarriages is limited. Cultured naïve pluripotent embryonic stem cells (ESC) differentiate into extra-embryonic endoderm/extraembryonic endoderm (XEN) or formative pluripotent ESC, during the period emulating maximal miscarriage of peri-implantation development. In previous reports using small marker sets, hyperosmotic sorbitol, or retinoic acid (RA) decreased naïve pluripotency and increased XEN by FACS quantitation.

METHODS

Bulk and single cell (sc)RNAseq analyses of two cultured ESC lines was done, corroborated by qPCR. Transcriptomic responses were analyzed of cultured ESC stressed by Sorbitol, with Leukemia inhibitory factor (LIF + ; stemness growth factor), RA without LIF to control for XEN induction, and compared with normal differentiation (LIF - , ND).

RESULTS

Sorbitol and RA increase subpopulations of 2-cell embryo-like (2CEL) and XEN sub-lineages; primitive, parietal, and visceral endoderm (VE) cells and suppress formative pluripotency, imbalancing alternate lineage choices of initial naïve pluripotent cultured ESC compared with ND. Although bulk RNAseq and gene ontology (GO) group analyses suggest that stress induces anterior VE-head organizer and placental markers, scRNAseq reveals relatively few cells. But VE and placental markers/cells were in adjacent stressed cell clusters in the UMAP, like recent, normal UMAP of conceptuses. UMAPs show that dose-dependent stress overrides stemness to force premature lineage imbalance.

CONCLUSIONS

Hyperosmotic stress, and other toxicological stresses, like drugs with active ingredient RA, may cause premature, lineage imbalance, resulting in miscarriages or birth defects.

AI-powered simulation-based inference of a genuinely spatial-stochastic gene regulation model of early mouse embryogenesis

Understanding how multicellular organisms reliably orchestrate cell-fate decisions is a central challenge in developmental biology, particularly in early mammalian development, where tissue-level differentiation arises from seemingly cell-autonomous mechanisms. In this study, we present a multi-scale, spatial-stochastic simulation framework for mouse embryogenesis, focusing on inner cell mass (ICM) differentiation into epiblast (EPI) and primitive endoderm (PRE) at the blastocyst stage. Our framework models key regulatory and tissue-scale interactions in a biophysically realistic fashion, capturing the inherent stochasticity of intracellular gene expression and intercellular signaling, while efficiently simulating these processes by advancing event-driven simulation techniques. Leveraging the power of Simulation-Based Inference (SBI) through the AI-driven Sequential Neural Posterior Estimation (SNPE) algorithm, we conduct a large-scale Bayesian inferential analysis to identify parameter sets that faithfully reproduce experimentally observed features of ICM specification. Our results reveal mechanistic insights into how the combined action of autocrine and paracrine FGF4 signaling coordinates stochastic gene expression at the cellular scale to achieve robust and reproducible ICM patterning at the tissue scale. We further demonstrate that the ICM exhibits a specific time window of sensitivity to exogenous FGF4, enabling lineage proportions to be adjusted based on timing and dosage, thereby extending current experimental findings and providing quantitative predictions for both mutant and wild-type ICM systems. Notably, FGF4 signaling not only ensures correct EPI-PRE lineage proportions but also enhances ICM resilience to perturbations, reducing fate-proportioning errors by 10-20% compared to a purely cell-autonomous system. Additionally, we uncover a surprising role for variability in intracellular initial conditions, showing that high gene-expression heterogeneity can improve both the accuracy and precision of cell-fate proportioning, which remains robust when fewer than 25% of the ICM population experiences perturbed initial conditions. Our work offers a comprehensive, spatial-stochastic description of the biochemical processes driving ICM differentiation and identifies the necessary conditions for its robust unfolding. It also provides a framework for future exploration of similar spatial-stochastic systems in developmental biology.

PI3K/AKT signaling controls ICM maturation and proper epiblast and primitive endoderm specification in mice

The inner cell mass (ICM) of early mouse embryos is specified into epiblast (Epi) and primitive endoderm (PrE) lineages during blastocyst formation. The antagonistic transcription factors (TFs) NANOG and GATA-binding protein 6 (GATA6) in combination with fibroblast growth factor (FGF)/extracellular-signal-regulated kinase (ERK) signaling are central actors in ICM fate choice. However, what initiates the specification of ICM progenitors into Epi or PrE and whether other factors are involved in this process has not been fully understood yet. Here, we show that phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT) is constitutively active during preimplantation development. Using pharmacological inhibition, we demonstrate that PI3K/AKT enables the formation of a functional ICM capable of giving rise to both the Epi and the PrE: it maintains the expression of the TF NANOG, which specifies the Epi, and confers responsiveness to FGF4, which is essential for PrE specification. Our work thus identifies PI3K/AKT signaling as an upstream regulator controlling the molecular events required for both Epi and PrE specification.

Cpt1a Drives primed-to-naïve pluripotency transition through lipid remodeling

Metabolism has been implicated in cell fate determination, particularly through epigenetic modifications. Similarly, lipid remodeling also plays a role in regulating cell fate. Here, we present comprehensive lipidomics analysis during BMP4-driven primed to naive pluripotency transition or BiPNT and demonstrate that lipid remodeling plays an essential role. We further identify Cpt1a as a rate-limiting factor in BiPNT, driving lipid remodeling and metabolic reprogramming while simultaneously increasing intracellular acetyl-CoA levels and enhancing H3K27ac at chromatin open sites. Perturbation of BiPNT by histone acetylation inhibitors suppresses lipid remodeling and pluripotency transition. Together, our study suggests that lipid remodeling promotes pluripotency transitions and further regulates cell fate decisions, implicating Cpt1a as a critical regulator between primed-naive cell fate control.

An integrated approach identifies the molecular underpinnings of murine anterior visceral endoderm migration

The anterior visceral endoderm (AVE) differs from the surrounding visceral endoderm (VE) in its migratory behavior and ability to restrict primitive streak formation to the opposite side of the mouse embryo. To characterize the molecular bases for the unique properties of the AVE, we combined single-cell RNA sequencing of the VE prior to and during AVE migration with phosphoproteomics, high-resolution live-imaging, and short-term lineage labeling and intervention. This identified the transient nature of the AVE with attenuation of "anteriorizing" gene expression as cells migrate and the emergence of heterogeneities in transcriptional states relative to the AVE's position. Using cell communication analysis, we identified the requirement of semaphorin signaling for normal AVE migration. Lattice light-sheet microscopy showed that Sema6D mutants have abnormalities in basal projections and migration speed. These findings point to a tight coupling between transcriptional state and position of the AVE and identify molecular controllers of AVE migration.

FoxO transcription factors actuate the formative pluripotency specific gene expression programme

Naïve pluripotency is sustained by a self-reinforcing gene regulatory network (GRN) comprising core and naïve pluripotency-specific transcription factors (TFs). Upon exiting naïve pluripotency, embryonic stem cells (ESCs) transition through a formative post-implantation-like pluripotent state, where they acquire competence for lineage choice. However, the mechanisms underlying disengagement from the naïve GRN and initiation of the formative GRN are unclear. Here, we demonstrate that phosphorylated AKT acts as a gatekeeper that prevents nuclear localisation of FoxO TFs in naïve ESCs. PTEN-mediated reduction of AKT activity upon exit from naïve pluripotency allows nuclear entry of FoxO TFs, enforcing a cell fate transition by binding and activating formative pluripotency-specific enhancers. Indeed, FoxO TFs are necessary and sufficient for the activation of the formative pluripotency-specific GRN. Our work uncovers a pivotal role for FoxO TFs in establishing formative post-implantation pluripotency, a critical early embryonic cell fate transition.

Integrin signaling in pluripotent cells acts as a gatekeeper of mouse germline entry

Primordial germ cells (PGCs) are the precursors of gametes and the sole mechanism by which animals transmit genetic information across generations. In the mouse embryo, the transcriptional and epigenetic regulation of PGC specification has been extensively characterized. However, the initial event that triggers the soma-germline segregation remains poorly understood. Here, we uncover a critical role for the basement membrane in regulating germline entry. We show that PGCs arise in a region of the mouse embryo that lacks contact with the basement membrane, and the addition of exogenous extracellular matrix (ECM) inhibits both PGC and PGC-like cell (PGCLC) specification in mouse embryos and stem cell models, respectively. Mechanistically, we demonstrate that the engagement of β1 integrin with laminin blocks PGCLC specification by preventing the Wnt signaling-dependent down-regulation of the PGC transcriptional repressor, Otx2. In this way, the physical segregation of cells away from the basement membrane acts as a morphogenetic fate switch that controls the soma-germline bifurcation.

Evolutionary origin of the chordate nervous system revealed by amphioxus developmental trajectories

The common ancestor of all vertebrates had a highly sophisticated nervous system, but questions remain about the evolution of vertebrate neural cell types. The amphioxus, a chordate that diverged before the origin of vertebrates, can inform vertebrate evolution. Here we develop and analyse a single-cell RNA-sequencing dataset from seven amphioxus embryo stages to understand chordate cell type evolution and to study vertebrate neural cell type origins. We identified many new amphioxus cell types, including homologues to the vertebrate hypothalamus and neurohypophysis, rooting the evolutionary origin of these structures. On the basis of ancestor-descendant reconstruction of cell trajectories of the amphioxus and other species, we inferred expression dynamics of transcription factor genes throughout embryogenesis and identified three ancient developmental routes forming chordate neurons. We characterized cell specification at the mechanistic level and generated mutant lines to examine the function of five key transcription factors involved in neural specification. Our results show three developmental origins for the vertebrate nervous system: an anterior FoxQ2-dependent mechanism that is deeply conserved in invertebrates, a less-conserved route leading to more posterior neurons in the vertebrate spinal cord and a mechanism for specifying neuromesoderm progenitors that is restricted to chordates. The evolution of neuromesoderm progenitors may have led to a dramatic shift in posterior neural and mesodermal cell fate decisions and the body elongation process in a stem chordate.

Non-apical mitoses contribute to cell delamination during mouse gastrulation

During the epithelial-mesenchymal transition driving mouse embryo gastrulation, cells divide more frequently at the primitive streak, and half of those divisions happen away from the apical pole. These observations suggest that non-apical mitoses might play a role in cell delamination. We aim to uncover and challenge the molecular determinants of mitosis position in different regions of the epiblast through computational modeling and pharmacological treatments of embryos and stem cell-based epiblast spheroids. Blocking basement membrane degradation at the streak has no impact on the asymmetry in mitosis frequency and position. By contrast, disturbance of the actomyosin cytoskeleton or cell cycle dynamics elicits ectopic non-apical mitosis and shows that the streak region is characterized by local relaxation of the actomyosin cytoskeleton and less stringent regulation of cell division. These factors are essential for normal dynamics at the streak and favor cell delamination from the epiblast.

Dynamic intrauterine crosstalk promotes porcine embryo implantation during early pregnancy

Dynamic crosstalk between the embryo and mother is crucial during implantation. Here, we comprehensively profile the single-cell transcriptome of pig peri-implantation embryos and corresponding maternal endometrium, identifying 4 different lineages in embryos and 13 cell types in the endometrium. Cell-specific gene expression characterizes 4 distinct trophectoderm subpopulations, showing development from undifferentiated trophectoderm to polar and mural trophectoderm. Dynamic expression of genes in different types of endometrial cells illustrates their molecular response to embryos during implantation. Then, we developed a novel tool, ExtraCellTalk, generating an overall dynamic map of maternal-foetal crosstalk using uterine luminal proteins as bridges. Through cross-species comparisons, we identified a conserved RBP4/STRA6 pathway in which embryonic-derived RBP4 could target the STRA6 receptor on stromal cells to regulate the interaction with other endometrial cells. These results provide insight into the maternal-foetal crosstalk during embryo implantation and represent a valuable resource for further studies to improve embryo implantation.

Polycomb function in early mouse development

Epigenetic factors are crucial for ensuring proper chromatin dynamics during the initial stages of embryo development. Among these factors, the Polycomb group (PcG) of proteins plays a key role in establishing correct transcriptional programmes during mouse embryogenesis. PcG proteins are classified into two complexes: Polycomb repressive complex 1 (PRC1) and PRC2. Both complexes decorate histone proteins with distinct post-translational modifications (PTMs) that are predictive of a silent transcriptional chromatin state. In recent years, a critical adaptation of the classical techniques to analyse chromatin profiles and to study biochemical interactions at low-input resolution has allowed us to deeply explore PcG molecular mechanisms in the very early stages of mouse embryo development- from fertilisation to gastrulation, and from zygotic genome activation (ZGA) to specific lineages differentiation. These advancements provide a foundation for a deeper understanding of the fundamental role Polycomb complexes play in early development and have elucidated the mechanistic dynamics of PRC1 and PRC2. In this review, we discuss the functions and molecular mechanisms of both PRC1 and PRC2 during early mouse embryo development, integrating new studies with existing knowledge. Furthermore, we highlight the molecular functionality of Polycomb complexes from ZGA through gastrulation, with a particular focus on non-canonical imprinted and bivalent genes, and Hox cluster regulation.

Protocol for generating mouse morula-like cells resembling 8- to 16-cell stage embryo cells

Generating cell types with properties of embryo cells with full developmental potential is of great biological importance. Here, we present a protocol for generating mouse morula-like cells (MLCs) resembling 8- to 16-cell stage embryo cells. We describe steps for induction, via increasing Stat3 activation, and the isolation of MLCs. We then detail procedures for segregating MLCs into blastocyst cell fates and how to create embryo-like structures from them. This system provides a stem-cell-based embryo model to study early embryo development. For complete details on the use and execution of this protocol, please refer to Li et al.1.

Lineage regulators TFAP2C and NR5A2 function as bipotency activators in totipotent embryos

During the first lineage segregation, a mammalian totipotent embryo differentiates into the inner cell mass (ICM) and trophectoderm (TE). However, how transcription factors (TFs) regulate this earliest cell-fate decision in vivo remains elusive, with their regulomes primarily inferred from cultured cells. Here, we investigated the TF regulomes during the first lineage specification in early mouse embryos, spanning the pre-initiation, initiation, commitment, and maintenance phases. Unexpectedly, we found that TFAP2C, a trophoblast regulator, bound and activated both early TE and inner cell mass (ICM) genes at the totipotent (two- to eight-cell) stages ('bipotency activation'). Tfap2c deficiency caused downregulation of early ICM genes, including Nanog, Nr5a2, and Tdgf1, and early TE genes, including Tfeb and Itgb5, in eight-cell embryos. Transcription defects in both ICM and TE lineages were also found in blastocysts, accompanied by increased apoptosis and reduced cell numbers in ICMs. Upon trophoblast commitment, TFAP2C left early ICM genes but acquired binding to late TE genes in blastocysts, where it co-bound with CDX2, and later to extra-embryonic ectoderm (ExE) genes, where it cooperatively co-occupied with the former ICM regulator SOX2. Finally, 'bipotency activation' in totipotent embryos also applied to a pluripotency regulator NR5A2, which similarly bound and activated both ICM and TE lineage genes at the eight-cell stage. These data reveal a unique transcription circuity of totipotency underpinned by highly adaptable lineage regulators.

3D reconstruction of a gastrulating human embryo

Progress in understanding early human development has been impeded by the scarcity of reference datasets from natural embryos, particularly those with spatial information during crucial stages like gastrulation. We conducted high-resolution spatial transcriptomics profiling on 38,562 spots from 62 transverse sections of an intact Carnegie stage (CS) 8 human embryo. From this spatial transcriptomic dataset, we constructed a 3D model of the CS8 embryo, in which a range of cell subtypes are identified, based on gene expression patterns and positional register, along the anterior-posterior, medial-lateral, and dorsal-ventral axis in the embryo. We further characterized the lineage trajectories of embryonic and extra-embryonic tissues and associated regulons and the regionalization of signaling centers and signaling activities that underpin lineage progression and tissue patterning during gastrulation. Collectively, the findings of this study provide insights into gastrulation and post-gastrulation development of the human embryo.

Basal delamination during mouse gastrulation primes pluripotent cells for differentiation

The blueprint of the mammalian body plan is laid out during gastrulation, when a trilaminar embryo is formed. This process entails a burst of proliferation, the ingression of embryonic epiblast cells at the primitive streak, and their priming toward primitive streak fates. How these different events are coordinated remains unknown. Here, we developed and characterized a 3D culture of self-renewing mouse embryonic cells that captures the main transcriptional and architectural features of the early gastrulating mouse epiblast. Using this system in combination with microfabrication and in vivo experiments, we found that proliferation-induced crowding triggers delamination of cells that express high levels of the apical polarity protein aPKC. Upon delamination, cells become more sensitive to Wnt signaling and upregulate the expression of primitive streak markers such as Brachyury. This mechanistic coupling between ingression and differentiation ensures that the right cell types become specified at the right place during embryonic development.

Unsolved morphogenesis problems and the hidden order

In this work, the morphogenesis mechanisms are considered from the complexity perspective. It is shown that both morphogenesis and the functioning of organs should be unstable in the case of short-range interaction potentials. The repeatability of forms during evolution is a strong argument for its directionality. The formation of organs during evolution can occur only in the presence of a priori information about the structure of such an organ. The focus of the discussion is not merely on constraining potential possibilities but on the concept of directed evolution itself. A morphogenesis model was constructed based on nontrivial quantum effects. These interaction effects between biologically important molecules ensure the accurate synthesis of cells, tissues, and organs.

Single-cell RNA sequencing of mutant whole mouse embryos: from the epiblast to the end of gastrulation

Over the last decade, single-cell approaches have become the gold standard for studying gene expression dynamics, cell heterogeneity, and cell states within samples. Before single-cell advances, the feasibility of capturing the dynamic cellular landscape and rapid cell transitions during early development was limited. In this paper, we designed a robust pipeline to perform single-cell and nuclei analysis on mouse embryos from E6.5 to E8, corresponding to the onset and completion of gastrulation. Gastrulation is a fundamental process during development that establishes the three germinal layers: mesoderm, ectoderm, and endoderm, which are essential for organogenesis. Extensive literature is available on single-cell omics applied to WT perigastrulating embryos. However, single-cell analysis of mutant embryos is still scarce and often limited to FACS-sorted populations. This is partially due to the technical constraints associated with the need for genotyping, timed pregnancies, the count of embryos with desired genotypes per pregnancy, and the number of cells per embryo at these stages. Here, we present a methodology designed to overcome these limitations. This method establishes breeding and timed pregnancy guidelines to achieve a higher chance of synchronized pregnancies with desired genotypes. Optimization steps in the embryo isolation process coupled with FAST genotyping protocol (3 hours) allow for microdroplet-based single-cell to be performed on the same day, ensuring the high viability of cells and robust results. We also include guidelines for optimal nuclei isolations from embryos. Thus, these approaches increase the feasibility of single-cell approaches of mutant embryos at the gastrulation stage. We anticipate this method will facilitate the analysis of how mutations shape the cellular landscape of the gastrula.

SUMMARY

We establish a pipeline for high-quality single-cell and nuclei suspensions of gastrulating mouse embryos for sequencing of single cells and nuclei.

The effect of data transformation on low-dimensional integration of single-cell RNA-seq

BACKGROUND

Recent developments in single-cell RNA sequencing have opened up a multitude of possibilities to study tissues at the level of cellular populations. However, the heterogeneity in single-cell sequencing data necessitates appropriate procedures to adjust for technological limitations and various sources of noise when integrating datasets from different studies. While many analysis procedures employ various preprocessing steps, they often overlook the importance of selecting and optimizing the employed data transformation methods.

RESULTS

This work investigates data transformation approaches used in single-cell clustering analysis tools and their effects on batch integration analysis. In particular, we compare 16 transformations and their impact on the low-dimensional representations, aiming to reduce the batch effect and integrate multiple single-cell sequencing data. Our results show that data transformations strongly influence the results of single-cell clustering on low-dimensional data space, such as those generated by UMAP or PCA. Moreover, these changes in low-dimensional space significantly affect trajectory analysis using multiple datasets, as well. However, the performance of the data transformations greatly varies across datasets, and the optimal method was different for each dataset. Additionally, we explored how data transformation impacts the analysis of deep feature encodings using deep neural network-based models, including autoencoder-based models and proto-typical networks. Data transformation also strongly affects the outcome of deep neural network models.

CONCLUSIONS

Our findings suggest that the batch effect and noise in integrative analysis are highly influenced by data transformation. Low-dimensional features can integrate different batches well when proper data transformation is applied. Furthermore, we found that the batch mixing score on low-dimensional space can guide the selection of the optimal data transformation. In conclusion, data preprocessing is one of the most crucial analysis steps and needs to be cautiously considered in the integrative analysis of multiple scRNA-seq datasets.

Annelid adult cell type diversity and their pluripotent cellular origins

Many annelids can regenerate missing body parts or reproduce asexually, generating all cell types in adult stages. However, the putative adult stem cell populations involved in these processes, and the diversity of cell types generated by them, are still unknown. To address this, we recover 75,218 single cell transcriptomes of the highly regenerative and asexually-reproducing annelid Pristina leidyi. Our results uncover a rich cell type diversity including annelid specific types as well as novel types. Moreover, we characterise transcription factors and gene networks that are expressed specifically in these populations. Finally, we uncover a broadly abundant cluster of putative stem cells with a pluripotent signature. This population expresses well-known stem cell markers such as vasa, piwi and nanos homologues, but also shows heterogeneous expression of differentiated cell markers and their transcription factors. We find conserved expression of pluripotency regulators, including multiple chromatin remodelling and epigenetic factors, in piwi+ cells. Finally, lineage reconstruction analyses reveal computational differentiation trajectories from piwi+ cells to diverse adult types. Our data reveal the cell type diversity of adult annelids by single cell transcriptomics and suggest that a piwi+ cell population with a pluripotent stem cell signature is associated with adult cell type differentiation.

Integrative analysis of transcriptomic and epigenomic data reveals distinct patterns for developmental and housekeeping gene regulation

BACKGROUND

Regulation of transcription is central to the emergence of new cell types during development, and it often involves activation of genes via proximal and distal regulatory regions. The activity of regulatory elements is determined by transcription factors (TFs) and epigenetic marks, but despite extensive mapping of such patterns, the extraction of regulatory principles remains challenging.

RESULTS

Here we study differentially and similarly expressed genes along with their associated epigenomic profiles, chromatin accessibility and DNA methylation, during lineage specification at gastrulation in mice. Comparison of the three lineages allows us to identify genomic and epigenomic features that distinguish the two classes of genes. We show that differentially expressed genes are primarily regulated by distal elements, while similarly expressed genes are controlled by proximal housekeeping regulatory programs. Differentially expressed genes are relatively isolated within topologically associated domains, while similarly expressed genes tend to be located in gene clusters. Transcription of differentially expressed genes is associated with differentially open chromatin at distal elements including enhancers, while that of similarly expressed genes is associated with ubiquitously accessible chromatin at promoters.

CONCLUSION

Based on these associations of (linearly) distal genes' transcription start sites (TSSs) and putative enhancers for developmental genes, our findings allow us to link putative enhancers to their target promoters and to infer lineage-specific repertoires of putative driver transcription factors, within which we define subgroups of pioneers and co-operators.

Mapping single-cell developmental potential in health and disease with interpretable deep learning

Single-cell RNA sequencing (scRNA-seq) has transformed our understanding of cell fate in developmental systems. However, identifying the molecular hallmarks of potency - the capacity of a cell to differentiate into other cell types - has remained challenging. Here, we introduce CytoTRACE 2, an interpretable deep learning framework for characterizing potency and differentiation states on an absolute scale from scRNA-seq data. Across 31 human and mouse scRNA-seq datasets encompassing 28 tissue types, CytoTRACE 2 outperformed existing methods for recovering experimentally determined potency levels and differentiation states covering the entire range of cellular ontogeny. Moreover, it reconstructed the temporal hierarchy of mouse embryogenesis across 62 timepoints; identified pan-tissue expression programs that discriminate major potency levels; and facilitated discovery of cellular phenotypes in cancer linked to survival and immunotherapy resistance. Our results illuminate a fundamental feature of cell biology and provide a broadly applicable platform for delineating single-cell differentiation landscapes in health and disease.

Distinct pathways drive anterior hypoblast specification in the implanting human embryo

Development requires coordinated interactions between the epiblast, which generates the embryo proper; the trophectoderm, which generates the placenta; and the hypoblast, which forms both the anterior signalling centre and the yolk sac. These interactions remain poorly understood in human embryogenesis because mechanistic studies have only recently become possible. Here we examine signalling interactions post-implantation using human embryos and stem cell models of the epiblast and hypoblast. We find anterior hypoblast specification is NODAL dependent, as in the mouse. However, while BMP inhibits anterior signalling centre specification in the mouse, it is essential for its maintenance in human. We also find contrasting requirements for BMP in the naive pre-implantation epiblast of mouse and human embryos. Finally, we show that NOTCH signalling is important for human epiblast survival. Our findings of conserved and species-specific factors that drive these early stages of embryonic development highlight the strengths of comparative species studies.

A single-cell time-lapse of mouse prenatal development from gastrula to birth

The house mouse (Mus musculus) is an exceptional model system, combining genetic tractability with close evolutionary affinity to humans1,2. Mouse gestation lasts only 3 weeks, during which the genome orchestrates the astonishing transformation of a single-cell zygote into a free-living pup composed of more than 500 million cells. Here, to establish a global framework for exploring mammalian development, we applied optimized single-cell combinatorial indexing3 to profile the transcriptional states of 12.4 million nuclei from 83 embryos, precisely staged at 2- to 6-hour intervals spanning late gastrulation (embryonic day 8) to birth (postnatal day 0). From these data, we annotate hundreds of cell types and explore the ontogenesis of the posterior embryo during somitogenesis and of kidney, mesenchyme, retina and early neurons. We leverage the temporal resolution and sampling depth of these whole-embryo snapshots, together with published data4-8 from earlier timepoints, to construct a rooted tree of cell-type relationships that spans the entirety of prenatal development, from zygote to birth. Throughout this tree, we systematically nominate genes encoding transcription factors and other proteins as candidate drivers of the in vivo differentiation of hundreds of cell types. Remarkably, the most marked temporal shifts in cell states are observed within one hour of birth and presumably underlie the massive physiological adaptations that must accompany the successful transition of a mammalian fetus to life outside the womb.

Tracking early mammalian organogenesis - prediction and validation of differentiation trajectories at whole organism scale

Early organogenesis represents a key step in animal development, during which pluripotent cells diversify to initiate organ formation. Here, we sampled 300,000 single-cell transcriptomes from mouse embryos between E8.5 and E9.5 in 6-h intervals and combined this new dataset with our previous atlas (E6.5-E8.5) to produce a densely sampled timecourse of >400,000 cells from early gastrulation to organogenesis. Computational lineage reconstruction identified complex waves of blood and endothelial development, including a new programme for somite-derived endothelium. We also dissected the E7.5 primitive streak into four adjacent regions, performed scRNA-seq and predicted cell fates computationally. Finally, we defined developmental state/fate relationships by combining orthotopic grafting, microscopic analysis and scRNA-seq to transcriptionally determine cell fates of grafted primitive streak regions after 24 h of in vitro embryo culture. Experimentally determined fate outcomes were in good agreement with computationally predicted fates, demonstrating how classical grafting experiments can be revisited to establish high-resolution cell state/fate relationships. Such interdisciplinary approaches will benefit future studies in developmental biology and guide the in vitro production of cells for organ regeneration and repair.

OCT4 is expressed in extraembryonic endoderm stem (XEN) cell progenitors during somatic cell reprogramming

During development, progenitors of embryonic stem (ES) and extraembryonic endoderm stem (XEN) cells are concomitantly specified within the inner cell mass (ICM) of the mouse blastocyst. Similarly, XEN cells are induced (iXEN cells) alongside induced pluripotent stem (iPS) cells following overexpression of Oct4, Sox2, Klf4 and Myc (OSKM) during somatic cell reprogramming. It is unclear how or why this cocktail produces both stem cell types, but OCT4 has been associated with non-pluripotent outcomes. In this report, we show that, during OSKM reprogramming, many individual Oct4-GFP-expressing cells are fated to become iXEN cells. Interestingly, SKM alone was also sufficient to induce iXEN cell formation, likely via activation of endogenous Oct4. These observations indicate that iXEN cell formation is not strictly an artifact of Oct4 overexpression. Moreover, our results suggest that a pathway to XEN may be an integral feature of establishing pluripotency during reprogramming, as in early embryo development.

The Nucleosome Remodelling and Deacetylation complex coordinates the transcriptional response to lineage commitment in pluripotent cells

As cells exit the pluripotent state and begin to commit to a specific lineage they must activate genes appropriate for that lineage while silencing genes associated with pluripotency and preventing activation of lineage-inappropriate genes. The Nucleosome Remodelling and Deacetylation (NuRD) complex is essential for pluripotent cells to successfully undergo lineage commitment. NuRD controls nucleosome density at regulatory sequences to facilitate transcriptional responses, and also has been shown to prevent unscheduled transcription (transcriptional noise) in undifferentiated pluripotent cells. How these activities combine to ensure cells engage a gene expression program suitable for successful lineage commitment has not been determined. Here, we show that NuRD is not required to silence all genes. Rather, it restricts expression of genes primed for activation upon exit from the pluripotent state, but maintains them in a transcriptionally permissive state in self-renewing conditions, which facilitates their subsequent activation upon exit from naïve pluripotency. We further show that NuRD coordinates gene expression changes, which acts to maintain a barrier between different stable states. Thus NuRD-mediated chromatin remodelling serves multiple functions, including reducing transcriptional noise, priming genes for activation and coordinating the transcriptional response to facilitate lineage commitment.

A relay velocity model infers cell-dependent RNA velocity

RNA velocity provides an approach for inferring cellular state transitions from single-cell RNA sequencing (scRNA-seq) data. Conventional RNA velocity models infer universal kinetics from all cells in an scRNA-seq experiment, resulting in unpredictable performance in experiments with multi-stage and/or multi-lineage transition of cell states where the assumption of the same kinetic rates for all cells no longer holds. Here we present cellDancer, a scalable deep neural network that locally infers velocity for each cell from its neighbors and then relays a series of local velocities to provide single-cell resolution inference of velocity kinetics. In the simulation benchmark, cellDancer shows robust performance in multiple kinetic regimes, high dropout ratio datasets and sparse datasets. We show that cellDancer overcomes the limitations of existing RNA velocity models in modeling erythroid maturation and hippocampus development. Moreover, cellDancer provides cell-specific predictions of transcription, splicing and degradation rates, which we identify as potential indicators of cell fate in the mouse pancreas.

Evolution of hierarchy and irreversibility in theoretical cell differentiation model

The process of cell differentiation in multicellular organisms is characterized by hierarchy and irreversibility in many cases. However, the conditions and selection pressures that give rise to these characteristics remain poorly understood. By using a mathematical model, here we show that the network of differentiation potency (differentiation diagram) becomes necessarily hierarchical and irreversible by increasing the number of terminally differentiated states under certain conditions. The mechanisms generating these characteristics are clarified using geometry in the cell state space. The results demonstrate that the hierarchical organization and irreversibility can manifest independently of direct selection pressures associated with these characteristics, instead they appear to evolve as byproducts of selective forces favoring a diversity of differentiated cell types. The study also provides a new perspective on the structure of gene regulatory networks that produce hierarchical and irreversible differentiation diagrams. These results indicate some constraints on cell differentiation, which are expected to provide a starting point for theoretical discussion of the implicit limits and directions of evolution in multicellular organisms.

Dissecting embryonic and extraembryonic lineage crosstalk with stem cell co-culture

Embryogenesis necessitates harmonious coordination between embryonic and extraembryonic tissues. Although stem cells of both embryonic and extraembryonic origins have been generated, they are grown in different culture conditions. In this study, utilizing a unified culture condition that activates the FGF, TGF-β, and WNT pathways, we have successfully derived embryonic stem cells (FTW-ESCs), extraembryonic endoderm stem cells (FTW-XENs), and trophoblast stem cells (FTW-TSCs) from the three foundational tissues of mouse and cynomolgus monkey (Macaca fascicularis) blastocysts. This approach facilitates the co-culture of embryonic and extraembryonic stem cells, revealing a growth inhibition effect exerted by extraembryonic endoderm cells on pluripotent cells, partially through extracellular matrix signaling. Additionally, our cross-species analysis identified both shared and unique transcription factors and pathways regulating FTW-XENs. The embryonic and extraembryonic stem cell co-culture strategy offers promising avenues for developing more faithful embryo models and devising more developmentally pertinent differentiation protocols.

Multifaceted SOX2-chromatin interaction underpins pluripotency progression in early embryos

Pioneer transcription factors (TFs), such as OCT4 and SOX2, play crucial roles in pluripotency regulation. However, the master TF-governed pluripotency regulatory circuitry was largely inferred from cultured cells. In this work, we investigated SOX2 binding from embryonic day 3.5 (E3.5) to E7.5 in the mouse. In E3.5 inner cell mass (ICM), SOX2 regulates the ICM-trophectoderm program but is dispensable for opening global enhancers. Instead, SOX2 occupies preaccessible enhancers in part opened by early-stage expressing TFs TFAP2C and NR5A2. SOX2 then widely redistributes when cells adopt naive and formative pluripotency by opening enhancers or poising them for rapid future activation. Hence, multifaceted pioneer TF-enhancer interaction underpins pluripotency progression in embryos, including a distinctive state in E3.5 ICM that bridges totipotency and pluripotency.

Gestational ozone inhalation elicits maternal cardiac dysfunction and transcriptional changes to placental pericytes and endothelial cells

Ozone (O3) is a criteria air pollutant with the most frequent incidence of exceeding air quality standards. Inhalation of O3 is known to cause lung inflammation and consequent systemic health effects, including endothelial dysfunction. Epidemiologic data have shown that gestational exposure to air pollutants correlates with complications of pregnancy, including low birth weight, intrauterine growth deficiency, preeclampsia, and premature birth. Mechanisms underlying how air pollution may facilitate or exacerbate gestational complications remain poorly defined. The current study sought to uncover how gestational O3 exposure impacted maternal cardiovascular function, as well as the development of the placenta. Pregnant mice were exposed to 1PPM O3 or a sham filtered air (FA) exposure for 4 h on gestational day (GD) 10.5, and evaluated for cardiac function via echocardiography on GD18.5. Echocardiography revealed a significant reduction in maternal stroke volume and ejection fraction in maternally exposed dams. To examine the impact of maternal O3 exposure on the maternal-fetal interface, placentae were analyzed by single-cell RNA sequencing analysis. Mid-gestational O3 exposure led to significant differential expression of 4021 transcripts compared with controls, and pericytes displayed the greatest transcriptional modulation. Pathway analysis identified extracellular matrix organization to be significantly altered after the exposure, with the greatest modifications in trophoblasts, pericytes, and endothelial cells. This study provides insights into potential molecular processes during pregnancy that may be altered due to the inhalation of environmental toxicants.

Transcriptomic analysis of PADI4 target genes during multi-lineage differentiation of embryonic stem cells

During mammalian embryo development, pluripotent epiblast cells diversify into the three primary germ layers, which will later give rise to all fetal and adult tissues. These processes involve profound transcriptional and epigenetic changes that require precise coordination. Peptidylarginine deiminase IV (PADI4) is a transcriptional regulator that is strongly associated with inflammation and carcinogenesis but whose physiological roles are less well understood. We previously found that Padi4 expression is associated with pluripotency. Here, we examined the role of PADI4 in maintaining the multi-lineage differentiation potential of mouse embryonic stem (ES) cells. Using bulk and single-cell transcriptomic analyses of embryoid bodies (EBs) derived from Padi4 knock-out (Padi4-KO) mouse ES cells, we find that PADI4 loss impairs mesoderm diversification and differentiation of cardimyocytes and endothelial cells. Additionally, Padi4 deletion leads to concerted downregulation of genes associated with polarized growth, sterol metabolism and the extracellular matrix (ECM). This study indicates a requirement for Padi4 in the specification of the mesodermal lineage and reports the Padi4 associated transcriptome, providing a platform for understanding the physiological functions of Padi4 in development and homeostasis. This article is part of the Theo Murphy meeting issue 'The virtues and vices of protein citrullination'.

In vitro generation of mouse morula-like cells

Generating cells with the molecular and functional properties of embryo cells and with full developmental potential is an aim with fundamental biological significance. Here we report the in vitro generation of mouse transient morula-like cells (MLCs) via the manipulation of signaling pathways. MLCs are molecularly distinct from embryonic stem cells (ESCs) and cluster instead with embryo 8- to 16-cell stage cells. A single MLC can generate a blastoid, and the efficiency increases to 80% when 8-10 MLCs are used. MLCs make embryoids directly, efficiently, and within 4 days. Transcriptomic analysis shows that day 4-5 MLC-derived embryoids contain the cell types found in natural embryos at early gastrulation. Furthermore, MLCs introduced into morulae segregate into epiblast (EPI), primitive endoderm (PrE), and trophectoderm (TE) fates in blastocyst chimeras and have a molecular signature indistinguishable from that of host embryo cells. These findings represent the generation of cells that are molecularly and functionally similar to the precursors of the first three cell lineages of the embryo.

Single-cell, whole-embryo phenotyping of mammalian developmental disorders

Mouse models are a critical tool for studying human diseases, particularly developmental disorders1. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse2. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing3 to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions4,5. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be 'decomposable' through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

Transcriptional network governing extraembryonic endoderm cell fate choice

The mechanism by which transcription factor (TF) network instructs cell-type-specific transcriptional programs to drive primitive endoderm (PrE) progenitors to commit to parietal endoderm (PE) versus visceral endoderm (VE) cell fates remains poorly understood. To address the question, we analyzed the single-cell transcriptional signatures defining PrE, PE, and VE cell states during the onset of the PE-VE lineage bifurcation. By coupling with the epigenomic comparison of active enhancers unique to PE and VE cells, we identified GATA6, SOX17, and FOXA2 as central regulators for the lineage divergence. Transcriptomic analysis of cXEN cells, an in vitro model for PE cells, after the acute depletion of GATA6 or SOX17 demonstrated that these factors induce Mycn, imparting the self-renewal properties of PE cells. Concurrently, they suppress the VE gene program, including key genes like Hnf4a and Ttr, among others. We proceeded with RNA-seq analysis on cXEN cells with FOXA2 knockout, in conjunction with GATA6 or SOX17 depletion. We found FOXA2 acts as a potent suppressor of Mycn while simultaneously activating the VE gene program. The antagonistic gene regulatory activities of GATA6/SOX17 and FOXA2 in promoting alternative cell fates, and their physical co-bindings at the enhancers provide molecular insights to the plasticity of the PrE lineage. Finally, we show that the external cue, BMP signaling, promotes the VE cell fate by activation of VE TFs and repression of PE TFs including GATA6 and SOX17. These data reveal a putative core gene regulatory module that underpins PE and VE cell fate choice.

Omics Views of Mechanisms for Cell Fate Determination in Early Mammalian Development

During mammalian preimplantation development, a totipotent zygote undergoes several cell cleavages and two rounds of cell fate determination, ultimately forming a mature blastocyst. Along with compaction, the establishment of apicobasal cell polarity breaks the symmetry of an embryo and guides subsequent cell fate choice. Although the lineage segregation of the inner cell mass (ICM) and trophectoderm (TE) is the first symbol of cell differentiation, several molecules have been shown to bias the early cell fate through their inter-cellular variations at much earlier stages, including the 2- and 4-cell stages. The underlying mechanisms of early cell fate determination have long been an important research topic. In this review, we summarize the molecular events that occur during early embryogenesis, as well as the current understanding of their regulatory roles in cell fate decisions. Moreover, as powerful tools for early embryogenesis research, single-cell omics techniques have been applied to both mouse and human preimplantation embryos and have contributed to the discovery of cell fate regulators. Here, we summarize their applications in the research of preimplantation embryos, and provide new insights and perspectives on cell fate regulation.

SGAE: single-cell gene association entropy for revealing critical states of cell transitions during embryonic development

The critical point or pivotal threshold of cell transition occurs in early embryonic development when cell differentiation culminates in its transition to specific cell fates, at which the cell population undergoes an abrupt and qualitative shift. Revealing such critical points of cell transitions can track cellular heterogeneity and shed light on the molecular mechanisms of cell differentiation. However, precise detection of critical state transitions proves challenging when relying on single-cell RNA sequencing data due to their inherent sparsity, noise, and heterogeneity. In this study, diverging from conventional methods like differential gene analysis or static techniques that emphasize classification of cell types, an innovative computational approach, single-cell gene association entropy (SGAE), is designed for the analysis of single-cell RNA-seq data and utilizes gene association information to reveal critical states of cell transitions. More specifically, through the translation of gene expression data into local SGAE scores, the proposed SGAE can serve as an index to quantitatively assess the resilience and critical properties of genetic regulatory networks, consequently detecting the signal of cell transitions. Analyses of five single-cell datasets for embryonic development demonstrate that the SGAE method achieves better performance in facilitating the characterization of a critical phase transition compared with other existing methods. Moreover, the SGAE value can effectively discriminate cellular heterogeneity over time and performs well in the temporal clustering of cells. Besides, biological functional analysis also indicates the effectiveness of the proposed approach.

Ribosomal RNA 2'-O-methylation dynamics impact cell fate decisions

Translational regulation impacts both pluripotency maintenance and cell differentiation. To what degree the ribosome exerts control over this process remains unanswered. Accumulating evidence has demonstrated heterogeneity in ribosome composition in various organisms. 2'-O-methylation (2'-O-me) of rRNA represents an important source of heterogeneity, where site-specific alteration of methylation levels can modulate translation. Here, we examine changes in rRNA 2'-O-me during mouse brain development and tri-lineage differentiation of human embryonic stem cells (hESCs). We find distinct alterations between brain regions, as well as clear dynamics during cortex development and germ layer differentiation. We identify a methylation site impacting neuronal differentiation. Modulation of its methylation levels affects ribosome association of the fragile X mental retardation protein (FMRP) and is accompanied by an altered translation of WNT pathway-related mRNAs. Together, these data identify ribosome heterogeneity through rRNA 2'-O-me during early development and differentiation and suggest a direct role for ribosomes in regulating translation during cell fate acquisition.

Inhibition of Wnt activity improves peri-implantation development of somatic cell nuclear transfer embryos

Somatic cell nuclear transfer (SCNT) can reprogram differentiated somatic cells into totipotency. Although pre-implantation development of SCNT embryos has greatly improved, most SCNT blastocysts are still arrested at the peri-implantation stage, and the underlying mechanism remains elusive. Here, we develop a 3D in vitro culture system for SCNT peri-implantation embryos and discover that persistent Wnt signals block the naïve-to-primed pluripotency transition of epiblasts with aberrant H3K27me3 occupancy, which in turn leads to defects in epiblast transformation events and subsequent implantation failure. Strikingly, manipulating Wnt signals can attenuate the pluripotency transition and H3K27me3 deposition defects in epiblasts and achieve up to a 9-fold increase in cloning efficiency. Finally, single-cell RNA-seq analysis reveals that Wnt inhibition markedly enhances the lineage developmental trajectories of SCNT blastocysts during peri-implantation development. Overall, these findings reveal diminished potentials of SCNT blastocysts for lineage specification and validate a critical peri-implantation barrier for SCNT embryos.

Optogenetic control of Wnt signaling models cell-intrinsic embryogenic patterning using 2D human pluripotent stem cell culture

In embryonic stem cell (ESC) models for early development, spatially and temporally varying patterns of signaling and cell types emerge spontaneously. However, mechanistic insight into this dynamic self-organization is limited by a lack of methods for spatiotemporal control of signaling, and the relevance of signal dynamics and cell-to-cell variability to pattern emergence remains unknown. Here, we combine optogenetic stimulation, imaging and transcriptomic approaches to study self-organization of human ESCs (hESC) in two-dimensional (2D) culture. Morphogen dynamics were controlled via optogenetic activation of canonical Wnt/β-catenin signaling (optoWnt), which drove broad transcriptional changes and mesendoderm differentiation at high efficiency (>99% cells). When activated within cell subpopulations, optoWnt induced cell self-organization into distinct epithelial and mesenchymal domains, mediated by changes in cell migration, an epithelial to mesenchymal-like transition and TGFβ signaling. Furthermore, we demonstrate that such optogenetic control of cell subpopulations can be used to uncover signaling feedback mechanisms between neighboring cell types. These findings reveal that cell-to-cell variability in Wnt signaling is sufficient to generate tissue-scale patterning and establish a hESC model system for investigating feedback mechanisms relevant to early human embryogenesis.

Plakoglobin is a mechanoresponsive regulator of naive pluripotency

Biomechanical cues are instrumental in guiding embryonic development and cell differentiation. Understanding how these physical stimuli translate into transcriptional programs will provide insight into mechanisms underlying mammalian pre-implantation development. Here, we explore this type of regulation by exerting microenvironmental control over mouse embryonic stem cells. Microfluidic encapsulation of mouse embryonic stem cells in agarose microgels stabilizes the naive pluripotency network and specifically induces expression of Plakoglobin (Jup), a vertebrate homolog of β-catenin. Overexpression of Plakoglobin is sufficient to fully re-establish the naive pluripotency gene regulatory network under metastable pluripotency conditions, as confirmed by single-cell transcriptome profiling. Finally, we find that, in the epiblast, Plakoglobin was exclusively expressed at the blastocyst stage in human and mouse embryos - further strengthening the link between Plakoglobin and naive pluripotency in vivo. Our work reveals Plakoglobin as a mechanosensitive regulator of naive pluripotency and provides a paradigm to interrogate the effects of volumetric confinement on cell-fate transitions.

RBBP4 is an epigenetic barrier for the induced transition of pluripotent stem cells into totipotent 2C-like cells

Cellular totipotency is critical for whole-organism generation, yet how totipotency is established remains poorly illustrated. Abundant transposable elements (TEs) are activated in totipotent cells, which is critical for embryonic totipotency. Here, we show that the histone chaperone RBBP4, but not its homolog RBBP7, is indispensable for maintaining the identity of mouse embryonic stem cells (mESCs). Auxin-induced degradation of RBBP4, but not RBBP7, reprograms mESCs to the totipotent 2C-like cells. Also, loss of RBBP4 enhances transition from mESCs to trophoblast cells. Mechanistically, RBBP4 binds to the endogenous retroviruses (ERVs) and functions as an upstream regulator by recruiting G9a to deposit H3K9me2 on ERVL elements, and recruiting KAP1 to deposit H3K9me3 on ERV1/ERVK elements, respectively. Moreover, RBBP4 facilitates the maintenance of nucleosome occupancy at the ERVK and ERVL sites within heterochromatin regions through the chromatin remodeler CHD4. RBBP4 depletion leads to the loss of the heterochromatin marks and activation of TEs and 2C genes. Together, our findings illustrate that RBBP4 is required for heterochromatin assembly and is a critical barrier for inducing cell fate transition from pluripotency to totipotency.

Genetic control of the dynamic transcriptional response to immune stimuli and glucocorticoids at single-cell resolution

Synthetic glucocorticoids, such as dexamethasone, have been used as a treatment for many immune conditions, such as asthma and, more recently, severe COVID-19. Single-cell data can capture more fine-grained details on transcriptional variability and dynamics to gain a better understanding of the molecular underpinnings of inter-individual variation in drug response. Here, we used single-cell RNA-seq to study the dynamics of the transcriptional response to glucocorticoids in activated peripheral blood mononuclear cells from 96 African American children. We used novel statistical approaches to calculate a mean-independent measure of gene expression variability and a measure of transcriptional response pseudotime. Using these approaches, we showed that glucocorticoids reverse the effects of immune stimulation on both gene expression mean and variability. Our novel measure of gene expression response dynamics, based on the diagonal linear discriminant analysis, separated individual cells by response status on the basis of their transcriptional profiles and allowed us to identify different dynamic patterns of gene expression along the response pseudotime. We identified genetic variants regulating gene expression mean and variability, including treatment-specific effects, and showed widespread genetic regulation of the transcriptional dynamics of the gene expression response.

A single cell transcriptomic fingerprint of stressed premature, imbalanced differentiation of embryonic stem cells

Cultured naïve pluripotent ESC differentiate into first lineage, XEN or second lineage, formative pluripotency. Hyperosmotic stress (sorbitol), like retinoic acid, decreases naive pluripotency and increases XEN in two ESC lines, as reported by bulk and scRNAseq, analyzed by UMAP. Sorbitol overrides pluripotency in two ESC lines as reported by bulk and scRNAseq, analyzed by UMAP. UMAP analyzed the effects of 5 stimuli - three stressed (200-300mM sorbitol with leukemia inhibitory factor +LIF) and two unstressed (+LIF, normal stemness-NS and -LIF, normal differentiation-ND). Sorbitol and RA decrease naive pluripotency and increase subpopulations of 2-cell embryo-like and XEN sub-lineages; primitive, parietal, and visceral endoderm (VE). Between the naïve pluripotency and primitive endoderm clusters is a stress-induced cluster with transient intermediate cells with higher LIF receptor signaling, with increased Stat3, Klf4, and Tbx3 expression. Sorbitol, like RA, also suppresses formative pluripotency, increasing lineage imbalance. Although bulk RNAseq and gene ontology group analyses suggest that stress induces head organizer and placental markers, scRNAseq reveals few cells. But VE and placental markers/cells were in adjacent clusters, like recent reports. UMAPs show that dose-dependent stress overrides stemness to force premature lineage imbalance. Hyperosmotic stress induces lineage imbalance, and other toxicological stresses, like drugs with RA, may cause lineage imbalance, resulting in miscarriages or birth defects.

Esrrb guides naive pluripotent cells through the formative transcriptional programme

During embryonic development, naive pluripotent epiblast cells transit to a formative state. The formative epiblast cells form a polarized epithelium, exhibit distinct transcriptional and epigenetic profiles and acquire competence to differentiate into all somatic and germline lineages. However, we have limited understanding of how the transition to a formative state is molecularly controlled. Here we used murine embryonic stem cell models to show that ESRRB is both required and sufficient to activate formative genes. Genetic inactivation of Esrrb leads to illegitimate expression of mesendoderm and extra-embryonic markers, impaired formative expression and failure to self-organize in 3D. Functionally, this results in impaired ability to generate formative stem cells and primordial germ cells in the absence of Esrrb. Computational modelling and genomic analyses revealed that ESRRB occupies key formative genes in naive cells and throughout the formative state. In so doing, ESRRB kickstarts the formative transition, leading to timely and unbiased capacity for multi-lineage differentiation.

Cross-activation of FGF, NODAL, and WNT pathways constrains BMP-signaling-mediated induction of the totipotent state in mouse embryonic stem cells

Embryonic stem cells (ESCs) are an attractive model to study the relationship between signaling and cell fates. Cultured mouse ESCs can exist in multiple states resembling distinct stages of early embryogenesis, such as totipotent, pluripotent, primed, and primitive endoderm. The signaling mechanisms regulating the totipotent state and coexistence of these states are poorly understood. Here we identify bone morphogenetic protein (BMP) signaling as an inducer of the totipotent state. However, we discover that BMP's role is constrained by the cross-activation of FGF, NODAL, and WNT pathways. We exploit this finding to enhance the proportion of totipotent cells by rationally inhibiting the cross-activated pathways. Single-cell mRNA sequencing reveals that induction of the totipotent state is accompanied by suppression of primed and primitive endoderm states. Furthermore, reprogrammed totipotent cells we generate in culture resemble totipotent cells of preimplantation embryo. Our findings reveal a BMP signaling mechanism regulating both the totipotent state and heterogeneity of ESCs.