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

The guinea pig serves as an alternative model to study human preimplantation development

Preimplantation development is an important window of human embryogenesis. However, ethical constraints and the limitations involved in studying human embryos often necessitate the use of alternative model systems. Here we identify the guinea pig as a promising small animal model to study human preimplantation development. Using single-cell RNA-sequencing, we generated an atlas of guinea pig preimplantation development, revealing its close resemblance to early human embryogenesis in terms of the timing of compaction, early-, mid- and late-blastocyst formation, and implantation, and the spatio-temporal expression of key lineage markers. We also show conserved roles of Hippo, MEK-ERK and JAK-STAT signalling. Furthermore, multi-species analysis highlights the spatio-temporal expression of conserved and divergent genes during preimplantation development and pluripotency. The guinea pig serves as a valuable animal model for advancing preimplantation development and stem cell research, and can be leveraged to better understand the longer-term impact of early exposures on offspring outcomes.

Origin, fate and function of extraembryonic tissues during mammalian development

Extraembryonic tissues have pivotal roles in morphogenesis and patterning of the early mammalian embryo. Developmental programmes mediated through signalling pathways and gene regulatory networks determine the sequence in which fate determination and lineage commitment of extraembryonic tissues take place, and epigenetic processes allow the memory of cell identity and state to be sustained throughout and beyond embryo development, even extending across generations. In this Review, we discuss the molecular and cellular mechanisms necessary for the different extraembryonic tissues to develop and function, from their initial specification up until the end of gastrulation, when the body plan of the embryo and the anatomical organization of its supporting extraembryonic structures are established. We examine the interaction between extraembryonic and embryonic tissues during early patterning and morphogenesis, and outline how epigenetic memory supports extraembryonic tissue development.

Chromatin accessibility landscape of mouse early embryos revealed by single-cell NanoATAC-seq2

In mammals, fertilized eggs undergo genome-wide epigenetic reprogramming to generate the organism. However, our understanding of epigenetic dynamics during preimplantation development at single-cell resolution remains incomplete. Here, we developed scNanoATAC-seq2, a single-cell assay for transposase-accessible chromatin using long-read sequencing for scarce samples. We present a detailed chromatin accessibility landscape of mouse preimplantation development, revealing distinct chromatin signatures in the epiblast, primitive endoderm, and trophectoderm during lineage segregation. Differences between zygotes and two-cell embryos highlight reprogramming in chromatin accessibility during the maternal-to-zygotic transition. Single-cell long-read sequencing enables in-depth analysis of chromatin accessibility in noncanonical imprinting, imprinted X chromosome inactivation, and low-mappability genomic regions, such as repetitive elements and paralogs. Our data provide insights into chromatin dynamics during mammalian preimplantation development and lineage differentiation.

Metabolic characterizations of PFOS-induced disruptions in early embryonic development

BACKGROUND

Perfluorooctane sulfonates (PFOS) are persistent environmental pollutants linked to developmental toxicity, but the mechanisms remain unclear. This study investigates the metabolic changes induced by PFOS exposure during early embryonic development and integrates metabolomic, transcriptomic, and molecular docking analyses to explore underlying mechanisms.

METHODS

Mouse embryoid bodies (mEBs) were exposed to PFOS for 2 days, 4 days and 6 days. Metabolomic profiling was conducted to identify differential metabolites. Transcriptomic data were integrated with metabolomics using Cytoscape to map metabolic pathway alterations. Molecular docking simulations were performed to assess PFOS binding to key enzymes.

RESULTS

PFOS exposure resulted in significant alterations in lipid (Erucic acid, L-carnitine), amino acid (L-methionine, creatine, hippuric acid, and spermine), and nucleotide metabolism (e.g., hypoxanthine). Integrated transcriptomic and metabolomic analysis revealed disrupted pathways included SLC25A20 regulated L-carnitine metabolism. Molecular docking simulations indicated that PFOS binds to methionine synthase and hypoxanthine guanine phosphoribosyl transferase, potentially inhibiting their function and disrupting metabolic homeostasis for L-methionine and hypoxanthine during embryonic development.

CONCLUSION

PFOS exposure disrupts key metabolic pathways critical for embryogenesis, including lipid, amino acid, and nucleotide metabolism. Molecular docking and transcriptomic integration highlight enzyme targeting as a potential mechanism of PFOS-induced developmental toxicity. These findings provide novel insights into the molecular and metabolic disruptions caused by PFOS, with implications for understanding its developmental toxicity.

Aging activates escape of the silent X chromosome in the female mouse hippocampus

Women live longer than men and exhibit less cognitive aging. The X chromosome contributes to sex differences, as females harbor an inactive X (Xi) and active X (Xa), in contrast to males with only an Xa. Thus, reactivation of silent Xi genes may contribute to sex differences. We use allele-specific, single-nucleus RNA sequencing to show that aging remodels transcription of the Xi and Xa across hippocampal cell types. Aging preferentially changed gene expression on the X's relative to autosomes. Select genes on the Xi underwent activation, with new escape across cells including in the dentate gyrus, critical to learning and memory. Expression of the Xi escapee Plp1, a myelin component, was increased in the aging hippocampus of female mice and parahippocampus of women. AAV-mediated Plp1 elevation in the dentate gyrus of aging male and female mice improved cognition. Understanding how the Xi may confer female advantage could lead to novel targets that counter brain aging and disease in both sexes.

Mouse totipotent blastomere-like cells model embryogenesis from zygotic genome activation to post implantation

Embryo development begins with zygotic genome activation (ZGA), eventually generating blastocysts for implantation. However, in vitro systems modeling the pre-implantation development are still absent and challenging. Here, we used mouse totipotent blastomere-like cells (TBLCs) to develop spontaneous differentiation and blastoid formation systems, respectively. We found Wnt signaling enabled the rapid expansion of TBLCs and the optimization of their culture medium. We successfully developed a TBLC-spontaneous differentiation system in which mouse TBLCs (mTBLCs) firstly converted into two types of ZGA-like cells (ZLCs) distinguished by Zscan4 expression. Surprisingly, Zscan4-, but not Zscan4+, ZLCs further passed through intermediate 4-cell and then 8-cell/morula stages to produce epiblast, primitive endoderm, and trophectoderm lineages. Significantly, single TBLCs underwent expansion, compaction, and polarization to efficiently generate blastocyst-like structures and even post-implantation egg-cylinder-like structures. Conclusively, we established TBLC-based differentiation and embryo-like structure formation systems to model early embryonic development, offering criteria for evaluating and understanding totipotency.

Generation of transient totipotent blastomere-like stem cells by short-term high-dose Pladienolide B treatment

As an alternative model for studying the dynamic process of early mammalian embryonic development, much progress has been made in using mouse embryonic stem cells (mESCs) to generate embryo-like structures, especially by modifying the starting cells. A previous study has demonstrated that totipotent blastomere-like cells (TBLCs) can be obtained by continuous treatment of mESCs with a low-dose splicing inhibitor, Pladienolide B (PlaB). However, these totipotent mESCs have limited proliferative capacity. Here, we report that short-term high-dose PlaB treatment can also induce mESCs to acquire totipotency. This treatment equips this novel type of stem cells with the ability to self-organize into blastoids and recapitulate key preimplantation developmental processes. Therefore, the stem cells are termed transient totipotent blastomere-like stem cells (tTBLCs). Transcriptome analysis showed that tTBLC blastoids bore similarities to mouse E3.5 blastocysts, E4.5 blastocysts, and TBLC blastoids. Additionally, we found that tTBLC blastoids could develop beyond the implantation stage, forming egg-cylinder-like structures both in vitro and in vivo. In summary, our research provides an alternative rapid and convenient method to generate the starting cells capable of developing into blastoids, which have immense application in various fields, not only in the basic study of early mouse embryogenesis but also in high-throughput drug screening.

Mapping cells through time and space with moscot

Single-cell genomic technologies enable the multimodal profiling of millions of cells across temporal and spatial dimensions. However, experimental limitations hinder the comprehensive measurement of cells under native temporal dynamics and in their native spatial tissue niche. Optimal transport has emerged as a powerful tool to address these constraints and has facilitated the recovery of the original cellular context1-4. Yet, most optimal transport applications are unable to incorporate multimodal information or scale to single-cell atlases. Here we introduce multi-omics single-cell optimal transport (moscot), a scalable framework for optimal transport in single-cell genomics that supports multimodality across all applications. We demonstrate the capability of moscot to efficiently reconstruct developmental trajectories of 1.7 million cells from mouse embryos across 20 time points. To illustrate the capability of moscot in space, we enrich spatial transcriptomic datasets by mapping multimodal information from single-cell profiles in a mouse liver sample and align multiple coronal sections of the mouse brain. We present moscot.spatiotemporal, an approach that leverages gene-expression data across both spatial and temporal dimensions to uncover the spatiotemporal dynamics of mouse embryogenesis. We also resolve endocrine-lineage relationships of delta and epsilon cells in a previously unpublished mouse, time-resolved pancreas development dataset using paired measurements of gene expression and chromatin accessibility. Our findings are confirmed through experimental validation of NEUROD2 as a regulator of epsilon progenitor cells in a model of human induced pluripotent stem cell islet cell differentiation. Moscot is available as open-source software, accompanied by extensive documentation.

Capture of Totipotency in Mouse Embryonic Stem Cells in the Absence of Pdzk1

Totipotent cells can differentiate into three lineages: the epiblast, primitive endoderm, and trophectoderm. Naturally, only early fertilized embryos possess totipotency, and they lose this ability as they develop. The expansion of stem cell differentiation potential has been a hot topic in developmental biology for years, particularly with respect to the generation totipotent-like stem cells. Here, the study describes the establishment of totipotency in embryonic stem cells (ESCs) via the deletion of a single gene, Pdzk1. Pdzk1-knockout (KO) ESCs substantially contribute to the fetus, placenta, and yolk sac in chimera assays but can also self-organize to form standard blastocyst-like structures containing the three lineages efficiently; thus, they exhibit full developmental potential as early blastomeres. Single-cell transcriptome and bulk RNA-seq comprehensive analyses revealed that Pdzk1-KO activates several lineage inducers (C1qa, C1qb, Fgf5, and Cdx2) to break down barriers between embryonic and extraembryonic tissues, making these lineages switch smoothly and resulting in a totipotent-like state. This versatile and scalable system provides a robust experimental model for differentiation potency and cell fate studies.

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.

Multimodal screen identifies noise-regulatory proteins

Gene-expression noise can influence cell-fate choices across pathology and physiology. However, a crucial question persists: do regulatory proteins or pathways exist that control noise independently of mean expression levels? Our integrative approach, combining single-cell RNA sequencing with proteomics and regulator enrichment analysis, identifies 32 putative noise regulators. SON, a nuclear speckle-associated protein, alters transcriptional noise without changing mean expression levels. Furthermore, SON's noise control can propagate to the protein level. Long-read and total RNA sequencing shows that SON's noise control does not significantly change isoform usage or splicing efficiency. Moreover, SON depletion reduces state switching in pluripotent mouse embryonic stem cells and impacts their fate choice during differentiation. Collectively, we demonstrate a class of proteins that control noise orthogonally to mean expression levels. This work serves as a proof of concept that can identify other functional noise regulators throughout development and disease progression.

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.

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.

The transcription factor GABPA is a master regulator of naive pluripotency

The establishment of naive pluripotency is a continuous process starting with the generation of inner cell mass (ICM) that then differentiates 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 naive 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 naive 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-TF regulatory network.

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.

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.

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.

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.

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, a robust pipeline was designed to perform single-cell and nuclei analysis on mouse embryos from embryonic day 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 wild-type 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, a methodology is presented 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 a same-day genotyping protocol (3 h) allow for microdroplet-based single-cell to be performed on the same day, ensuring the high viability of cells and robust results. This method further includes 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 that this method will facilitate the analysis of how mutations shape the cellular landscape of the gastrula.

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.