A XEN-like State Bridges Somatic Cells to Pluripotency during Chemical Reprogramming.

Somatic cells can be reprogrammed into pluripotent stem cells (PSCs) by using pure chemicals, providing a different paradigm to study somatic reprogramming. However, the cell fate dynamics and molecular events that occur during the chemical reprogramming process remain unclear. We now show that the chemical reprogramming process requires the early formation of extra-embryonic endoderm (XEN)-like cells and a late transition from XEN-like cells to chemically-induced (Ci)PSCs, a unique route that fundamentally differs from the pathway of transcription factor-induced reprogramming. Moreover, precise manipulation of the cell fate transition in a step-wise manner through the XEN-like state allows us to identify small-molecule boosters and establish a robust chemical reprogramming system with a yield up to 1,000-fold greater than that of the previously reported protocol. These findings demonstrate that chemical reprogramming is a promising approach to manipulate cell fates.

Modelling post-implantation human development to yolk sac blood emergence.

Implantation of the human embryo begins a critical developmental stage that comprises profound events including axis formation, gastrulation and the emergence of haematopoietic system1,2. Our mechanistic knowledge of this window of human life remains limited due to restricted access to in vivo samples for both technical and ethical reasons3-5. Stem cell models of human embryo have emerged to help unlock the mysteries of this stage6-16. Here we present a genetically inducible stem cell-derived embryoid model of early post-implantation human embryogenesis that captures the reciprocal codevelopment of embryonic tissue and the extra-embryonic endoderm and mesoderm niche with early haematopoiesis. This model is produced from induced pluripotent stem cells and shows unanticipated self-organizing cellular programmes similar to those that occur in embryogenesis, including the formation of amniotic cavity and bilaminar disc morphologies as well as the generation of an anterior hypoblast pole and posterior domain. The extra-embryonic layer in these embryoids lacks trophoblast and shows advanced multilineage yolk sac tissue-like morphogenesis that harbours a process similar to distinct waves of haematopoiesis, including the emergence of erythroid-, megakaryocyte-, myeloid- and lymphoid-like cells. This model presents an easy-to-use, high-throughput, reproducible and scalable platform to probe multifaceted aspects of human development and blood formation at the early post-implantation stage. It will provide a tractable human-based model for drug testing and disease modelling.

Spatial profiling of early primate gastrulation in utero.

Gastrulation controls the emergence of cellular diversity and axis patterning in the early embryo. In mammals, this transformation is orchestrated by dynamic signalling centres at the interface of embryonic and extraembryonic tissues1-3. Elucidating the molecular framework of axis formation in vivo is fundamental for our understanding of human development4-6 and to advance stem-cell-based regenerative approaches7. Here we illuminate early gastrulation of marmoset embryos in utero using spatial transcriptomics and stem-cell-based embryo models. Gaussian process regression-based 3D transcriptomes delineate the emergence of the anterior visceral endoderm, which is hallmarked by conserved (HHEX, LEFTY2, LHX1) and primate-specific (POSTN, SDC4, FZD5) factors. WNT signalling spatially coordinates the formation of the primitive streak in the embryonic disc and is counteracted by SFRP1 and SFRP2 to sustain pluripotency in the anterior domain. Amnion specification occurs at the boundaries of the embryonic disc through ID1, ID2 and ID3 in response to BMP signalling, providing a developmental rationale for amnion differentiation of primate pluripotent stem cells (PSCs). Spatial identity mapping demonstrates that primed marmoset PSCs exhibit the highest similarity to the anterior embryonic disc, whereas naive PSCs resemble the preimplantation epiblast. Our 3D transcriptome models reveal the molecular code of lineage specification in the primate embryo and provide an in vivo reference to decipher human development.

Embryo model completes gastrulation to neurulation and organogenesis.

Embryonic stem (ES) cells can undergo many aspects of mammalian embryogenesis in vitro1-5, but their developmental potential is substantially extended by interactions with extraembryonic stem cells, including trophoblast stem (TS) cells, extraembryonic endoderm stem (XEN) cells and inducible XEN (iXEN) cells6-11. Here we assembled stem cell-derived embryos in vitro from mouse ES cells, TS cells and iXEN cells and showed that they recapitulate the development of whole natural mouse embryo in utero up to day 8.5 post-fertilization. Our embryo model displays headfolds with defined forebrain and midbrain regions and develops a beating heart-like structure, a trunk comprising a neural tube and somites, a tail bud containing neuromesodermal progenitors, a gut tube, and primordial germ cells. This complete embryo model develops within an extraembryonic yolk sac that initiates blood island development. Notably, we demonstrate that the neurulating embryo model assembled from Pax6-knockout ES cells aggregated with wild-type TS cells and iXEN cells recapitulates the ventral domain expansion of the neural tube that occurs in natural, ubiquitous Pax6-knockout embryos. Thus, these complete embryoids are a powerful in vitro model for dissecting the roles of diverse cell lineages and genes in development. Our results demonstrate the self-organization ability of ES cells and two types of extraembryonic stem cells to reconstitute mammalian development through and beyond gastrulation to neurulation and early organogenesis.

A geometrical model of cell fate specification in the mouse blastocyst.

The lineage decision that generates the epiblast and primitive endoderm from the inner cell mass (ICM) is a paradigm for cell fate specification. Recent mathematics has formalized Waddington's landscape metaphor and proven that lineage decisions in detailed gene network models must conform to a small list of low-dimensional stereotypic changes called bifurcations. The most plausible bifurcation for the ICM is the so-called heteroclinic flip that we define and elaborate here. Our re-analysis of recent data suggests that there is sufficient cell movement in the ICM so the FGF signal, which drives the lineage decision, can be treated as spatially uniform. We thus extend the bifurcation model for a single cell to the entire ICM by means of a self-consistently defined time-dependent FGF signal. This model is consistent with available data and we propose additional dynamic experiments to test it further. This demonstrates that simplified, quantitative and intuitively transparent descriptions are possible when attention is shifted from specific genes to lineages. The flip bifurcation is a very plausible model for any situation where the embryo needs control over the relative proportions of two fates by a morphogen feedback.

CellRank 2: unified fate mapping in multiview single-cell data.

Single-cell RNA sequencing allows us to model cellular state dynamics and fate decisions using expression similarity or RNA velocity to reconstruct state-change trajectories; however, trajectory inference does not incorporate valuable time point information or utilize additional modalities, whereas methods that address these different data views cannot be combined or do not scale. Here we present CellRank 2, a versatile and scalable framework to study cellular fate using multiview single-cell data of up to millions of cells in a unified fashion. CellRank 2 consistently recovers terminal states and fate probabilities across data modalities in human hematopoiesis and endodermal development. Our framework also allows combining transitions within and across experimental time points, a feature we use to recover genes promoting medullary thymic epithelial cell formation during pharyngeal endoderm development. Moreover, we enable estimating cell-specific transcription and degradation rates from metabolic-labeling data, which we apply to an intestinal organoid system to delineate differentiation trajectories and pinpoint regulatory strategies.

Single-cell transcriptomic characterization of a gastrulating human embryo.

Gastrulation is the fundamental process in all multicellular animals through which the basic body plan is first laid down1-4. It is pivotal in generating cellular diversity coordinated with spatial patterning. In humans, gastrulation occurs in the third week after fertilization. Our understanding of this process in humans is relatively limited and based primarily on historical specimens5-8, experimental models9-12 or, more recently, in vitro cultured samples13-16. Here we characterize in a spatially resolved manner the single-cell transcriptional profile of an entire gastrulating human embryo, staged to be between 16 and 19 days after fertilization. We use these data to analyse the cell types present and to make comparisons with other model systems. In addition to pluripotent epiblast, we identified primordial germ cells, red blood cells and various mesodermal and endodermal cell types. This dataset offers a unique glimpse into a central but inaccessible stage of our development. This characterization provides new context for interpreting experiments in other model systems and represents a valuable resource for guiding directed differentiation of human cells in vitro.

Control of embryonic stem cell lineage commitment by core promoter factor, TAF3.

Deciphering the molecular basis of pluripotency is fundamental to our understanding of development and embryonic stem cell function. Here, we report that TAF3, a TBP-associated core promoter factor, is highly enriched in ES cells. In this context, TAF3 is required for endoderm lineage differentiation and prevents premature specification of neuroectoderm and mesoderm. In addition to its role in the core promoter recognition complex TFIID, genome-wide binding studies reveal that TAF3 localizes to a subset of chromosomal regions bound by CTCF/cohesin that are selectively associated with genes upregulated by TAF3. Notably, CTCF directly recruits TAF3 to promoter distal sites and TAF3-dependent DNA looping is observed between the promoter distal sites and core promoters occupied by TAF3/CTCF/cohesin. Together, our findings support a new role of TAF3 in mediating long-range chromatin regulatory interactions that safeguard the finely-balanced transcriptional programs underlying pluripotency.

HBO1 determines SMAD action in pluripotency and mesendoderm specification.

TGF-ß signaling family plays an essential role to regulate fate decisions in pluripotency and lineage specification. How the action of TGF-ß family signaling is intrinsically executed remains not fully elucidated. Here, we show that HBO1, a MYST histone acetyltransferase (HAT) is an essential cell intrinsic determinant for TGF-ß signaling in human embryonic stem cells (hESCs). HBO1-/- hESCs fail to response to TGF-ß signaling to maintain pluripotency and spontaneously differentiate into neuroectoderm. Moreover, HBO1 deficient hESCs show complete defect in mesendoderm specification in BMP4-triggered gastruloids or teratomas. Molecularly, HBO1 interacts with SMAD4 and co-binds the open chromatin labeled by H3K14ac and H3K4me3 in undifferentiated hESCs. Upon differentiation, HBO1/SMAD4 co-bind and maintain the mesoderm genes in BMP4-triggered mesoderm cells while lose chromatin occupancy in neural cells induced by dual-SMAD inhibition. Our data reveal an essential role of HBO1, a chromatin factor to determine the action of SMAD in both human pluripotency and mesendoderm specification.

ECM-integrin signalling instructs cellular position sensing to pattern the early mouse embryo.

Development entails patterned emergence of diverse cell types within the embryo. In mammals, cells positioned inside the embryo give rise to the inner cell mass (ICM), which eventually forms the embryo itself. Yet, the molecular basis of how these cells recognise their 'inside' position to instruct their fate is unknown. Here, we show that provision of extracellular matrix (ECM) to isolated embryonic cells induces ICM specification and alters the subsequent spatial arrangement between epiblast (EPI) and primitive endoderm (PrE) cells that emerge within the ICM. Notably, this effect is dependent on integrin ß1 activity and involves apical-to-basal conversion of cell polarity. We demonstrate that ECM-integrin activity is sufficient for 'inside' positional signalling and is required for correct EPI/PrE patterning. Thus, our findings highlight the significance of ECM-integrin adhesion in enabling position sensing by cells to achieve tissue patterning.

Leveling Waddington: the emergence of direct programming and the loss of cell fate hierarchies.

For decades, Waddington's concept of the 'epigenetic landscape' has served as an educative hierarchical model to illustrate the progressive restriction of cell differentiation potential during normal development. While still being highly valuable in the context of normal development, the Waddington model falls short of accommodating recent breakthroughs in cell programming. The advent of induced pluripotent stem (iPS) cells and advances in direct cell fate conversion (also known as transdifferentiation) suggest that somatic and pluripotent cell fates can be interconverted without transiting through distinct hierarchies. We propose a non-hierarchical 'epigenetic disc' model to explain such cell fate transitions, which provides an alternative landscape for modelling cell programming and reprogramming.

The emergent landscape of the mouse gut endoderm at single-cell resolution.

Here we delineate the ontogeny of the mammalian endoderm by generating 112,217 single-cell transcriptomes, which represent all endoderm populations within the mouse embryo until midgestation. We use graph-based approaches to model differentiating cells, which provides a spatio-temporal characterization of developmental trajectories and defines the transcriptional architecture that accompanies the emergence of the first (primitive or extra-embryonic) endodermal population and its sister pluripotent (embryonic) epiblast lineage. We uncover a relationship between descendants of these two lineages, in which epiblast cells differentiate into endoderm at two distinct time points-before and during gastrulation. Trajectories of endoderm cells were mapped as they acquired embryonic versus extra-embryonic fates and as they spatially converged within the nascent gut endoderm, which revealed these cells to be globally similar but retain aspects of their lineage history. We observed the regionalized identity of cells along the anterior-posterior axis of the emergent gut tube, which reflects their embryonic or extra-embryonic origin, and the coordinated patterning of these cells into organ-specific territories.

Genetic induction and mechanochemical propagation of a morphogenetic wave.

Tissue morphogenesis arises from coordinated changes in cell shape driven by actomyosin contractions. Patterns of gene expression regionalize cell behaviours by controlling actomyosin contractility. Here we report two modes of control over Rho1 and myosin II (MyoII) activation in the Drosophila endoderm. First, Rho1-MyoII are induced in a spatially restricted primordium via localized transcription of the G-protein-coupled receptor ligand Fog. Second, a tissue-scale wave of Rho1-MyoII activation and cell invagination progresses anteriorly away from the primordium. The wave does not require sustained gene transcription, and is not governed by regulated Fog delivery. Instead, MyoII inhibition blocks Rho1 activation and propagation, revealing a mechanical feedback driven by MyoII. We find that MyoII activation and invagination in each row of cells drives adhesion to the vitelline membrane mediated by integrins, apical spreading, MyoII activation and invagination in the next row. Endoderm morphogenesis thus emerges from local transcriptional initiation and a mechanically driven cycle of cell deformation.

Molecular control of macroscopic forces drives formation of the vertebrate hindgut.

The embryonic gut tube is a cylindrical structure from which the respiratory and gastrointestinal tracts develop1. Although the early emergence of the endoderm as an epithelial sheet2,3 and later morphogenesis of the definitive digestive and respiratory organs4-6 have been investigated, the intervening process of gut tube formation remains relatively understudied7,8. Here we investigate the molecular control of macroscopic forces underlying early morphogenesis of the gut tube in the chick embryo. The gut tube has been described as forming from two endodermal invaginations-the anterior intestinal portal (AIP) towards the rostral end of the embryo and the caudal intestinal portal (CIP) at the caudal end-that migrate towards one another, internalizing the endoderm until they meet at the yolk stalk (umbilicus in mammals)1,6. Migration of the AIP to form foregut has been descriptively characterized8,9, but the hindgut is likely to form by a distinct mechanism that has not been fully explained10. We find that the hindgut is formed by collective cell movements through a stationary CIP, rather than by movement of the CIP itself. Further, combining in vivo imaging, biophysics and mathematical modelling with molecular and embryological approaches, we identify a contractile force gradient that drives cell movements in the hindgut-forming endoderm, enabling tissue-scale posterior extension of the forming hindgut tube. The force gradient, in turn, is established in response to a morphogenic gradient of fibroblast growth factor signalling. As a result, we propose that an important positive feedback arises, whereby contracting cells draw passive cells from low to high fibroblast growth factor levels, recruiting them to contract and pull more cells into the elongating hindgut. In addition to providing insight into the early gut development, these findings illustrate how large-scale tissue level forces can be traced to developmental signals during vertebrate morphogenesis.

Multi-omics profiling of mouse gastrulation at single-cell resolution.

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

Modelling human hepato-biliary-pancreatic organogenesis from the foregut-midgut boundary.

Organogenesis is a complex and interconnected process that is orchestrated by multiple boundary tissue interactions1-7. However, it remains unclear how individual, neighbouring components coordinate to establish an integral multi-organ structure. Here we report the continuous patterning and dynamic morphogenesis of hepatic, biliary and pancreatic structures, invaginating from a three-dimensional culture of human pluripotent stem cells. The boundary interactions between anterior and posterior gut spheroids differentiated from human pluripotent stem cells enables retinoic acid-dependent emergence of hepato-biliary-pancreatic organ domains specified at the foregut-midgut boundary organoids in the absence of extrinsic factors. Whereas transplant-derived tissues are dominated by midgut derivatives, long-term-cultured microdissected hepato-biliary-pancreatic organoids develop into segregated multi-organ anlages, which then recapitulate early morphogenetic events including the invagination and branching of three different and interconnected organ structures, reminiscent of tissues derived from mouse explanted foregut-midgut culture. Mis-segregation of multi-organ domains caused by a genetic mutation in HES1 abolishes the biliary specification potential in culture, as seen in vivo8,9. In sum, we demonstrate that the experimental multi-organ integrated model can be established by the juxtapositioning of foregut and midgut tissues, and potentially serves as a tractable, manipulatable and easily accessible model for the study of complex human endoderm organogenesis.

A single-cell molecular map of mouse gastrulation and early organogenesis.

Across the animal kingdom, gastrulation represents a key developmental event during which embryonic pluripotent cells diversify into lineage-specific precursors that will generate the adult organism. Here we report the transcriptional profiles of 116,312 single cells from mouse embryos collected at nine sequential time points ranging from 6.5 to 8.5 days post-fertilization. We construct a molecular map of cellular differentiation from pluripotency towards all major embryonic lineages, and explore the complex events involved in the convergence of visceral and primitive streak-derived endoderm. Furthermore, we use single-cell profiling to show that Tal1-/- chimeric embryos display defects in early mesoderm diversification, and we thus demonstrate how combining temporal and transcriptional information can illuminate gene function. Together, this comprehensive delineation of mammalian cell differentiation trajectories in vivo represents a baseline for understanding the effects of gene mutations during development, as well as a roadmap for the optimization of in vitro differentiation protocols for regenerative medicine.

Developmental single-cell transcriptomics in the Lytechinus variegatus sea urchin embryo.

Using scRNA-seq coupled with computational approaches, we studied transcriptional changes in cell states of sea urchin embryos during development to the larval stage. Eighteen closely spaced time points were taken during the first 24 h of development of Lytechinus variegatus (Lv). Developmental trajectories were constructed using Waddington-OT, a computational approach to 'stitch' together developmental time points. Skeletogenic and primordial germ cell trajectories diverged early in cleavage. Ectodermal progenitors were distinct from other lineages by the 6th cleavage, although a small percentage of ectoderm cells briefly co-expressed endoderm markers that indicated an early ecto-endoderm cell state, likely in cells originating from the equatorial region of the egg. Endomesoderm cells also originated at the 6th cleavage and this state persisted for more than two cleavages, then diverged into distinct endoderm and mesoderm fates asynchronously, with some cells retaining an intermediate specification status until gastrulation. Seventy-nine out of 80 genes (99%) examined, and included in published developmental gene regulatory networks (dGRNs), are present in the Lv-scRNA-seq dataset and are expressed in the correct lineages in which the dGRN circuits operate.

Krüppel-like factor gene function in the ctenophore Mnemiopsis leidyi assessed by CRISPR/Cas9-mediated genome editing.

The Krüppel-like factor (Klf) gene family encodes transcription factors that play an important role in the regulation of stem cell proliferation, cell differentiation and development in bilaterians. Although Klf genes have been shown to specify functionally various cell types in non-bilaterian animals, their role in early-diverging animal lineages has not been assessed. Thus, the ancestral activity of these transcription factors in animal development is not well understood. The ctenophore Mnemiopsis leidyi has emerged as an important non-bilaterian model system for understanding early animal evolution. Here, we characterize the expression and functional role of Klf genes during M. leidyi embryogenesis. Zygotic Klf gene function was assessed with both CRISPR/Cas9-mediated genome editing and splice-blocking morpholino oligonucleotide knockdown approaches. Abrogation of zygotic Klf expression during M. leidyi embryogenesis resulted in abnormal development of several organs, including the pharynx, tentacle bulbs and apical organ. Our data suggest an ancient role for Klf genes in regulating endodermal patterning, possibly through regulation of cell proliferation.

Cell-cell communication through FGF4 generates and maintains robust proportions of differentiated cell types in embryonic stem cells.

During embryonic development and tissue homeostasis, reproducible proportions of differentiated cell types are specified from populations of multipotent precursor cells. Molecular mechanisms that enable both robust cell-type proportioning despite variable initial conditions in the precursor cells, and the re-establishment of these proportions upon perturbations in a developing tissue remain to be characterized. Here, we report that the differentiation of robust proportions of epiblast-like and primitive endoderm-like cells in mouse embryonic stem cell cultures emerges at the population level through cell-cell communication via a short-range fibroblast growth factor 4 (FGF4) signal. We characterize the molecular and dynamical properties of the communication mechanism and show how it controls both robust cell-type proportioning from a wide range of experimentally controlled initial conditions, as well as the autonomous re-establishment of these proportions following the isolation of one cell type. The generation and maintenance of reproducible proportions of discrete cell types is a new function for FGF signaling that might operate in a range of developing tissues.