Diverse Routes toward Early Somites in the Mouse Embryo

A single-cell and tissue-scale analysis suite resolves Mixl1's role in heart development

Perturbation studies using gene knockouts have become a key tool for understanding the roles of regulatory genes in development. However, large-scale studies dissecting the molecular role of development master regulators in every cell type throughout the embryo are technically challenging and scarce. Here, we systematically characterize the knockout effects of the key developmental regulators T/Brachyury and Mixl1 in gastrulation and early organogenesis using single-cell profiling of chimeric mouse embryos. For the analysis of these experimental data, we present COSICC, an effective suite of statistical tools to characterize perturbation effects in complex developing cell populations. We gain insights into T's role in lateral plate mesoderm, limb development, and posterior intermediate mesoderm specification. Furthermore, we generate Mixl1 -/- embryonic chimeras and reveal the role of this key transcription factor in discrete mesoderm lineages, in particular concerning developmental dysregulation of the recently identified juxta-cardiac field.

Extended culture of 2D gastruloids to model human mesoderm development

Micropatterned human pluripotent stem cells treated with BMP4 (two-dimensional (2D) gastruloids) are among the most widely used stem cell models for human gastrulation. Due to its simplicity and reproducibility, this system is ideal for high-throughput quantitative studies of tissue patterning and has led to many insights into the mechanisms of mammalian gastrulation. However, 2D gastruloids have been studied only up to about 2 days owing to a loss of organization beyond this time with earlier protocols. Here we report an extended 2D gastruloid model to up to 10 days. We discovered a phase of highly reproducible morphogenesis between 2 and 4 days during which directed migration from the primitive streak-like region gives rise to a mesodermal layer beneath an epiblast-like layer. Multiple types of mesoderm arise with striking spatial organization including lateral plate mesoderm-like cells on the colony border and paraxial mesoderm-like cells further inside the colony. Single-cell transcriptomics showed strong similarity of these cells to mesoderm in human and nonhuman primate embryos. Our results illustrate that extended culture of 2D gastruloids provides a powerful model for human mesoderm differentiation and morphogenesis.

STAT3 signaling enhances tissue expansion during postimplantation mouse development

Signal transducer and activator of transcription (STAT)3 signaling has been studied extensively using mouse embryonic stem cells. Zygotic deletion of Stat3 enables embryo implantation, but thereafter, mutants resemble non-affected littermates from the previous day until around mid-gestation. This probably results from the loss of serine-phosphorylated STAT3, the predominant form in early postimplantation embryonic tissues associated with mitochondrial activity. Bulk RNA sequencing of isolated mouse epiblasts confirmed developmental delay transcriptionally. Single-cell RNA sequencing revealed the exclusion of derivatives of Stat3 null embryonic stem cells exclusively from the erythroid lineage of mid-gestation chimeras. We show that Stat3 null embryonic stem cells can differentiate into erythroid and hematopoietic lineages in vitro but become outcompeted when mixed with wild-type cells. Our results implicate a role for STAT3 in the temporal control of embryonic progression, particularly in tissues requiring rapid cell division, such as postimplantation epiblast and hematopoietic lineages. Interestingly, mutations in STAT3 are associated with short stature in humans.

Dual genetic tracing demonstrates the heterogeneous differentiation and function of neuromesodermal progenitors in vivo

In recent decades, the traditional paradigm of three distinct germ layers formed during gastrulation has been revised with the identification of neuromesodermal progenitors (NMPs). These progenitors emerge during gastrulation and contribute to both the neural ectoderm, particularly the spinal cord, and the adjacent paraxial mesoderm [D. Henrique et al., Development 142, 2864-2875 (2015); R. J. Garriock et al., Development 142, 1628-1638 (2015); E. Tzouanacou et al., Dev. Cell 17, 365-376 (2009)]. However, effective genetic tools for lineage tracing and functional assessments of NMPs in vivo are currently lacking. Here, we developed a dual recombinase-mediated genetic system to specifically trace and ablate Brachyury+Sox2+ NMPs. Our genetic tracing and single-cell RNA sequencing analyses revealed that NMPs consist of three distinct unipotent and bipotent progenitor populations that progressively differentiate into neural and mesodermal fates. Genetic depletion of NMPs demonstrated their critical role in trunk and tail formation. This study provides in vivo genetic evidence supporting the heterogeneity of NMPs in terms of cell fate determination and their functional roles in the developing embryo.

Progress in understanding the vertebrate segmentation clock

The segmentation clock is a molecular oscillator that regulates the periodic formation of somites from the presomitic mesoderm during vertebrate embryogenesis. Synchronous oscillatory expression of a Hairy homologue or Hairy-related basic helix-loop-helix (bHLH) transcriptional repressor in presomitic mesoderm cells regulates periodic expression of downstream factors that control somite segmentation with a periodicity that varies across species. Although many of the key components of the clock have been identified and characterized, less is known about how the clock is synchronized across cells and how species-specific periodicity is achieved. Advances in live imaging, stem cell and organoid technologies, and synthetic approaches have started to uncover the detailed mechanisms underlying these aspects of somitogenesis, providing insight into how morphogenesis is coordinated in space and time during embryonic development.

Understanding developing kidneys and Wilms tumors one cell at a time

Single-cell sequencing-based techniques are revolutionizing all fields of biomedical sciences, including normal kidney development and how this is disturbed in the development of Wilms tumor. The many different techniques and the differences between them can obscure which technique is best used to answer which question. In this review we summarize the techniques currently available, discuss which have been used in kidney development or Wilms tumor context, and which techniques can or should be combined to maximize the increase in biological understanding we can get from them.

The ratio of Wnt signaling activity to Sox2 transcription factor levels predicts neuromesodermal fate potential

Neuromesodermal progenitors (NMPs) are a vertebrate cell type that contribute descendants to both the spinal cord and the mesoderm. The undifferentiated bipotential NMP state is maintained when both Wnt signaling is active and Sox2 is present. We used transgenic reporter lines to live-image both Wnt activity and Sox2 levels in NMPs and observed a unique cellular ratio in NMPs compared to NMP-derived mesoderm or neural tissue. We used this unique signature to identify the previously unknown anatomical position of a progenitor population that gives rise to the midline tissues of the floor plate of the spinal cord and the mesodermal notochord. Thus, quantification of the active Wnt signaling to Sox2 ratio can be used to predict and identify cells with neuromesodermal potential. We also developed the auxin inducible degron 2 system for use in zebrafish to test the temporal role that Sox2 plays during midline formation. We found ectopic Sox2 in the presence of Wnt activity holds cells in the undifferentiated floor plate/notochord progenitor state, and that degradation of the ectopic Sox2 is required for cells to adopt a notochord fate.

Early autonomous patterning of the anteroposterior axis in gastruloids

Minimal in vitro systems composed of embryonic stem cells (ESCs) have been shown to recapitulate the establishment of the anteroposterior (AP) axis. In contrast to the native embryo, ESC aggregates - such as gastruloids - can break symmetry, which is demarcated by polarization of the mesodermal marker T, autonomously without any localized external cues. However, associated earliest patterning events, such as the spatial restriction of cell fates and concomitant transcriptional changes, remain poorly understood. Here, we dissect the dynamics of AP axis establishment in mouse gastruloids, particularly before external Wnt stimulation. Through single-cell RNA sequencing, we identify key cell state transitions and the molecular signatures of T+ and T- populations underpinning AP polarization. We also show that this process is robust to modifications of aggregate size. Finally, transcriptomic comparison with the mouse embryo indicates that gastruloids develop similar mesendodermal cell types, despite initial differences in their primed pluripotent populations, which adopt a more mesenchymal state in lieu of an epiblast-like transcriptome. Hence, our findings suggest the possibility of alternate ESC states in vivo and in vitro that can converge onto similar cell fates.

Mapping lineage-traced cells across time points with moslin

Simultaneous profiling of single-cell gene expression and lineage history holds enormous potential for studying cellular decision-making. Recent computational approaches combine both modalities into cellular trajectories; however, they cannot make use of all available lineage information in destructive time-series experiments. Here, we present moslin, a Gromov-Wasserstein-based model to couple cellular profiles across time points based on lineage and gene expression information. We validate our approach in simulations and demonstrate on Caenorhabditis elegans embryonic development how moslin predicts fate probabilities and putative decision driver genes. Finally, we use moslin to delineate lineage relationships among transiently activated fibroblast states during zebrafish heart regeneration.

Individual variation in the emergence of anterior-to-posterior neural fates from human pluripotent stem cells

Variability between human pluripotent stem cell (hPSC) lines remains a challenge and opportunity in biomedicine. In this study, hPSC lines from multiple donors were differentiated toward neuroectoderm and mesendoderm lineages. We revealed dynamic transcriptomic patterns that delineate the emergence of these lineages, which were conserved across lines, along with individual line-specific transcriptional signatures that were invariant throughout differentiation. These transcriptomic signatures predicted an antagonism between SOX21-driven forebrain fates and retinoic acid-induced hindbrain fates. Replicate lines and paired adult tissue demonstrated the stability of these line-specific transcriptomic traits. We show that this transcriptomic variation in lineage bias had both genetic and epigenetic origins, aligned with the anterior-to-posterior structure of early mammalian development, and was present across a large collection of hPSC lines. These findings contribute to developing systematic analyses of PSCs to define the origin and consequences of variation in the early events orchestrating individual human development.

StaVia: spatially and temporally aware cartography with higher-order random walks for cell atlases

Single-cell atlases pose daunting computational challenges pertaining to the integration of spatial and temporal information and the visualization of trajectories across large atlases. We introduce StaVia, a computational framework that synergizes multi-faceted single-cell data with higher-order random walks that leverage the memory of cells' past states, fused with a cartographic Atlas View that offers intuitive graph visualization. This spatially aware cartography captures relationships between cell populations based on their spatial location as well as their gene expression and developmental stage. We demonstrate this using zebrafish gastrulation data, underscoring its potential to dissect complex biological landscapes in both spatial and temporal contexts.

Control of epiblast cell fate by mechanical cues

In amniotes, embryonic tissues originate from multipotent epiblast cells, arranged in a mosaic of presumptive territories. How these domains fated to specific lineages become segregated during body formation remains poorly understood. Using single cell RNA sequencing analysis and lineage tracing in the chicken embryo, we identify epiblast cells contributing descendants to the neural tube, somites and lateral plate after completion of gastrulation. We show that intercalation after cell division generates important movements of epiblast cells which lead to their relocation to different presumptive territories, explaining this broad spectrum of fates. This tissue remodeling phase is transient, being soon followed by the establishment of boundaries restricting cell movements therefore defining the presumptive territories of the epiblast. Finally, we find that the epiblast faces distinct mechanical constraints along the antero-posterior axis, leading to cell fate alterations when challenged. Together, we demonstrate the critical role of mechanical cues in epiblast fate determination.

A single-cell atlas of pig gastrulation as a resource for comparative embryology

Cell-fate decisions during mammalian gastrulation are poorly understood outside of rodent embryos. The embryonic disc of pig embryos mirrors humans, making them a useful proxy for studying gastrulation. Here we present a single-cell transcriptomic atlas of pig gastrulation, revealing cell-fate emergence dynamics, as well as conserved and divergent gene programs governing early porcine, primate, and murine development. We highlight heterochronicity in extraembryonic cell-types, despite the broad conservation of cell-type-specific transcriptional programs. We apply these findings in combination with functional investigations, to outline conserved spatial, molecular, and temporal events during definitive endoderm specification. We find early FOXA2 + /TBXT- embryonic disc cells directly form definitive endoderm, contrasting later-emerging FOXA2/TBXT+ node/notochord progenitors. Unlike mesoderm, none of these progenitors undergo epithelial-to-mesenchymal transition. Endoderm/Node fate hinges on balanced WNT and hypoblast-derived NODAL, which is extinguished upon endodermal differentiation. These findings emphasise the interplay between temporal and topological signalling in fate determination during gastrulation.

Mannose controls mesoderm specification and symmetry breaking in mouse gastruloids

Patterning and growth are fundamental features of embryonic development that must be tightly coordinated. To understand how metabolism impacts early mesoderm development, we used mouse embryonic stem-cell-derived gastruloids, that co-expressed glucose transporters with the mesodermal marker T/Bra. We found that the glucose mimic, 2-deoxy-D-glucose (2-DG), blocked T/Bra expression and abolished axial elongation in gastruloids. However, glucose removal did not phenocopy 2-DG treatment despite a decline in glycolytic intermediates. As 2-DG can also act as a competitive inhibitor of mannose in protein glycosylation, we added mannose together with 2-DG and found that it could rescue the mesoderm specification both in vivo and in vitro. We further showed that blocking production and intracellular recycling of mannose abrogated mesoderm specification. Proteomics analysis demonstrated that mannose reversed glycosylation of the Wnt pathway regulator, secreted frizzled receptor Frzb. Our study showed how mannose controls mesoderm specification in mouse gastruloids.

Exploring the formation mechanism of short-tailed phenotypes in animals using mutant mice with the TBXT gene c.G334T developed by CRISPR/Cas9

With the change in diet structure, individuals prefer to consume mutton with less fat. However, sheep tail has a lot of fat. We identified a breed of low-fat short-tailed sheep (i.e., Hulunbuir short-tailed sheep). It is necessary to develop an animal model that can promote research on the potential mechanisms of the short-tail phenotype in sheep, which results from the TBXT gene c.G334T mutation. To create animal models, we selected mice as experimental animals. Mouse embryos lacking the TBXT protein, which crucially regulates mouse embryonic development, cannot develop normally. We utilized CRISPR/Cas9 gene editing technology to generate site-specific mutation (c.G334T) in the TBXT gene of mice, and found that the mouse TBXT mutation (c.G334T) leads to a short-tail phenotype. Furthermore, we investigated the interaction between TBXT and Wnt signaling pathways. The expressions of TBXT, Axin2, Dkk1, Wnt3, Wnt3a, and Wnt5a were discovered to be significantly different between mutant embryos and wild embryos by obtaining mouse embryos at various developmental stages and examining the expression relationship between the TBXT and Wnt signaling pathway-related components in all of these embryos. Therefore, as a transcription factor, TBXT regulates the expression of the aforementioned Wnt signaling pathway components by forming a regulatory network for the normal development of mouse embryos. This study enriches the research on the functional role of the TBXT in the development of mouse embryos and the mechanism by which the short-tailed phenotype in sheep develops.

Gastrulation-stage gene expression in Nipbl+/- mouse embryos foreshadows the development of syndromic birth defects

In animal models, Nipbl deficiency phenocopies gene expression changes and birth defects seen in Cornelia de Lange syndrome, the most common cause of which is Nipbl haploinsufficiency. Previous studies in Nipbl+/- mice suggested that heart development is abnormal as soon as cardiogenic tissue is formed. To investigate this, we performed single-cell RNA sequencing on wild-type and Nipbl+/- mouse embryos at gastrulation and early cardiac crescent stages. Nipbl+/- embryos had fewer mesoderm cells than wild-type and altered proportions of mesodermal cell subpopulations. These findings were associated with underexpression of genes implicated in driving specific mesodermal lineages. In addition, Nanog was found to be overexpressed in all germ layers, and many gene expression changes observed in Nipbl+/- embryos could be attributed to Nanog overexpression. These findings establish a link between Nipbl deficiency, Nanog overexpression, and gene expression dysregulation/lineage misallocation, which ultimately manifest as birth defects in Nipbl+/- animals and Cornelia de Lange syndrome.

The Origin and Regulation of Neuromesodermal Progenitors (NMPs) in Embryos

Neuromesodermal progenitors (NMPs), serving as the common origin of neural and paraxial mesodermal development in a large part of the trunk, have recently gained significant attention because of their critical importance in the understanding of embryonic organogenesis and the design of in vitro models of organogenesis. However, the nature of NMPs at many essential points remains only vaguely understood or even incorrectly assumed. Here, we discuss the nature of NMPs, focusing on their dynamic migratory behavior during embryogenesis and the mechanisms underlying their neural vs. mesodermal fate choice. The discussion points include the following: (1) How the sinus rhomboidals is organized; the tissue where the neural or mesodermal fate choice of NMPs occurs. (2) NMPs originating from the broad posterior epiblast are associated with Sox2 N1 enhancer activity. (3) Tbx6-dependent Sox2 repression occurs during NMP-derived paraxial mesoderm development. (4) The nephric mesenchyme, a component of the intermediate mesoderm, was newly identified as an NMP derivative. (5) The transition of embryonic tissue development from tissue-specific progenitors in the anterior part to that from NMPs occurs at the forelimb bud axial level. (6) The coexpression of Sox2 and Bra in NMPs is conditional and is not a hallmark of NMPs. (7) The ability of the NMP pool to sustain axial embryo growth depends on Wnt3a signaling in the NMP population. Current in vitro models of NMPs are also critically reviewed.

Neuromesodermal specification during head-to-tail body axis formation

The anterior-to-posterior (head-to-tail) body axis is extraordinarily diverse among vertebrates but conserved within species. Body axis development requires a population of axial progenitors that resides at the posterior of the embryo to sustain elongation and is then eliminated once axis extension is complete. These progenitors occupy distinct domains in the posterior (tail-end) of the embryo and contribute to various lineages along the body axis. The subset of axial progenitors with neuromesodermal competency will generate both the neural tube (the precursor of the spinal cord), and the trunk and tail somites (producing the musculoskeleton) during embryo development. These axial progenitors are called Neuromesodermal Competent cells (NMCs) and Neuromesodermal Progenitors (NMPs). NMCs/NMPs have recently attracted interest beyond the field of developmental biology due to their clinical potential. In the mouse, the maintenance of neuromesodermal competency relies on a fine balance between a trio of known signals: Wnt/β-catenin, FGF signalling activity and suppression of retinoic acid signalling. These signals regulate the relative expression levels of the mesodermal transcription factor Brachyury and the neural transcription factor Sox2, permitting the maintenance of progenitor identity when co-expressed, and either mesoderm or neural lineage commitment when the balance is tilted towards either Brachyury or Sox2, respectively. Despite important advances in understanding key genes and cellular behaviours involved in these fate decisions, how the balance between mesodermal and neural fates is achieved remains largely unknown. In this chapter, we provide an overview of signalling and gene regulatory networks in NMCs/NMPs. We discuss mutant phenotypes associated with axial defects, hinting at the potential significant role of lesser studied proteins in the maintenance and differentiation of the progenitors that fuel axial elongation.

Gastrulation-stage gene expression in Nipbl +/- mouse embryos foreshadows the development of syndromic birth defects

In animal models, Nipbl-deficiency phenocopies gene expression changes and birth defects seen in Cornelia de Lange Syndrome (CdLS), the most common cause of which is Nipbl-haploinsufficiency. Previous studies in Nipbl+/- mice suggested that heart development is abnormal as soon as cardiogenic tissue is formed. To investigate this, we performed single-cell RNA-sequencing on wildtype (WT) and Nipbl+/- mouse embryos at gastrulation and early cardiac crescent stages. Nipbl+/- embryos had fewer mesoderm cells than WT and altered proportions of mesodermal cell subpopulations. These findings were associated with underexpression of genes implicated in driving specific mesodermal lineages. In addition, Nanog was found to be overexpressed in all germ layers, and many gene expression changes observed in Nipbl+/- embryos could be attributed to Nanog overexpression. These findings establish a link between Nipbl-deficiency, Nanog overexpression, and gene expression dysregulation/lineage misallocation, which ultimately manifest as birth defects in Nipbl+/- animals and CdLS.

Teaser

Gene expression changes during gastrulation of Nipbl-deficient mice shed light on early origins of structural birth defects.

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.

Notch signalling influences cell fate decisions and HOX gene induction in axial progenitors

The generation of the post-cranial embryonic body relies on the coordinated production of spinal cord neurectoderm and presomitic mesoderm cells from neuromesodermal progenitors (NMPs). This process is orchestrated by pro-neural and pro-mesodermal transcription factors that are co-expressed in NMPs together with Hox genes, which are essential for axial allocation of NMP derivatives. NMPs reside in a posterior growth region, which is marked by the expression of Wnt, FGF and Notch signalling components. Although the importance of Wnt and FGF in influencing the induction and differentiation of NMPs is well established, the precise role of Notch remains unclear. Here, we show that the Wnt/FGF-driven induction of NMPs from human embryonic stem cells (hESCs) relies on Notch signalling. Using hESC-derived NMPs and chick embryo grafting, we demonstrate that Notch directs a pro-mesodermal character at the expense of neural fate. We show that Notch also contributes to activation of HOX gene expression in human NMPs, partly in a non-cell-autonomous manner. Finally, we provide evidence that Notch exerts its effects via the establishment of a negative-feedback loop with FGF signalling.

Stabilized mosaic single-cell data integration using unshared features

Currently available single-cell omics technologies capture many unique features with different biological information content. Data integration aims to place cells, captured with different technologies, onto a common embedding to facilitate downstream analytical tasks. Current horizontal data integration techniques use a set of common features, thereby ignoring non-overlapping features and losing information. Here we introduce StabMap, a mosaic data integration technique that stabilizes mapping of single-cell data by exploiting the non-overlapping features. StabMap first infers a mosaic data topology based on shared features, then projects all cells onto supervised or unsupervised reference coordinates by traversing shortest paths along the topology. We show that StabMap performs well in various simulation contexts, facilitates 'multi-hop' mosaic data integration where some datasets do not share any features and enables the use of spatial gene expression features for mapping dissociated single-cell data onto a spatial transcriptomic reference.

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.

Gastrulation: Its Principles and Variations

As epiblast cells initiate development into various somatic cells, they undergo a large-scale reorganization, called gastrulation. The gastrulation of the epiblast cells produces three groups of cells: the endoderm layer, the collection of miscellaneous mesodermal tissues, and the ectodermal layer, which includes the neural, epidermal, and associated tissues. Most studies of gastrulation have focused on the formation of the tissues that provide the primary route for cell reorganization, that is, the primitive streak, in the chicken and mouse. In contrast, how gastrulation alters epiblast-derived cells has remained underinvestigated. This chapter highlights the regulation of cell and tissue fate via the gastrulation process. The roles and regulatory functions of neuromesodermal progenitors (NMPs) in the gastrulation process, elucidated in the last decade, are discussed in depth to resolve points of confusion. Chicken and mouse embryos, which form a primitive streak as the site of mesoderm precursor ingression, have been investigated extensively. However, primitive streak formation is an exception, even among amniotes. The roles of gastrulation processes in generating various somatic tissues will be discussed broadly.

From signalling to form: the coordination of neural tube patterning

The development of the vertebrate spinal cord involves the formation of the neural tube and the generation of multiple distinct cell types. The process starts during gastrulation, combining axial elongation with specification of neural cells and the formation of the neuroepithelium. Tissue movements produce the neural tube which is then exposed to signals that provide patterning information to neural progenitors. The intracellular response to these signals, via a gene regulatory network, governs the spatial and temporal differentiation of progenitors into specific cell types, facilitating the assembly of functional neuronal circuits. The interplay between the gene regulatory network, cell movement, and tissue mechanics generates the conserved neural tube pattern observed across species. In this review we offer an overview of the molecular and cellular processes governing the formation and patterning of the neural tube, highlighting how the remarkable complexity and precision of vertebrate nervous system arises. We argue that a multidisciplinary and multiscale understanding of the neural tube development, paired with the study of species-specific strategies, will be crucial to tackle the open questions.

The control of transitions along the main body axis

Although vertebrates display a large variety of forms and sizes, the mechanisms controlling the layout of the basic body plan are substantially conserved throughout the clade. Following gastrulation, head, trunk, and tail are sequentially generated through the continuous addition of tissue at the caudal embryonic end. Development of each of these major embryonic regions is regulated by a distinct genetic network. The transitions from head-to-trunk and from trunk-to-tail development thus involve major changes in regulatory mechanisms, requiring proper coordination to guarantee smooth progression of embryonic development. In this review, we will discuss the key cellular and embryological events associated with those transitions giving particular attention to their regulation, aiming to provide a cohesive outlook of this important component of vertebrate development.

Mapping oto-pharyngeal development in a human inner ear organoid model

Inner ear development requires the coordination of cell types from distinct epithelial, mesenchymal and neuronal lineages. Although we have learned much from animal models, many details about human inner ear development remain elusive. We recently developed an in vitro model of human inner ear organogenesis using pluripotent stem cells in a 3D culture, fostering the growth of a sensorineural circuit, including hair cells and neurons. Despite previously characterizing some cell types, many remain undefined. This study aimed to chart the in vitro development timeline of the inner ear organoid to understand the mechanisms at play. Using single-cell RNA sequencing at ten stages during the first 36 days of differentiation, we tracked the evolution from pluripotency to various ear cell types after exposure to specific signaling modulators. Our findings showcase gene expression that influences differentiation, identifying a plethora of ectodermal and mesenchymal cell types. We also discern aspects of the organoid model consistent with in vivo development, while highlighting potential discrepancies. Our study establishes the Inner Ear Organoid Developmental Atlas (IODA), offering deeper insights into human biology and improving inner ear tissue differentiation.

Eomes restricts Brachyury functions at the onset of mouse gastrulation

Mammalian specification of mesoderm and definitive endoderm (DE) is instructed by the two related Tbx transcription factors (TFs) Eomesodermin (Eomes) and Brachyury sharing partially redundant functions. Gross differences in mutant embryonic phenotypes suggest specific functions of each TF. To date, the molecular details of separated lineage-specific gene regulation by Eomes and Brachyury remain poorly understood. Here, we combine mouse embryonic and stem-cell-based analyses to delineate the non-overlapping, lineage-specific transcriptional activities. On a genome-wide scale, binding of both TFs overlaps at promoters of target genes but shows specificity for distal enhancer regions that is conferred by differences in Tbx DNA-binding motifs. The unique binding to enhancer sites instructs the specification of anterior mesoderm (AM) and DE by Eomes and caudal mesoderm by Brachyury. Remarkably, EOMES antagonizes BRACHYURY gene regulatory functions in coexpressing cells during early gastrulation to ensure the proper sequence of early AM and DE lineage specification followed by posterior mesoderm derivatives.

A unified pipeline for FISH spatial transcriptomics

High-throughput spatial transcriptomics has emerged as a powerful tool for investigating the spatial distribution of mRNA expression and its effects on cellular function. There is a lack of standardized tools for analyzing spatial transcriptomics data, leading many groups to write their own in-house tools that are often poorly documented and not generalizable. To address this, we have expanded and improved the starfish library and used those tools to create PIPEFISH, a semi-automated and generalizable pipeline that performs transcript annotation for fluorescence in situ hybridization (FISH)-based spatial transcriptomics. We used this pipeline to annotate transcript locations from three real datasets from three different common types of FISH image-based experiments, MERFISH, seqFISH, and targeted in situ sequencing (ISS), and verified that the results were high quality using the internal quality metrics of the pipeline and also a comparison with an orthogonal method of measuring RNA expression. PIPEFISH is a publicly available and open-source tool.

Oct1 cooperates with the Smad family of transcription factors to promote mesodermal lineage specification

The transition between pluripotent and tissue-specific states is a key aspect of development. Understanding the pathways driving these transitions will facilitate the engineering of properly differentiated cells for experimental and therapeutic uses. Here, we showed that during mesoderm differentiation, the transcription factor Oct1 activated developmental lineage-appropriate genes that were silent in pluripotent cells. Using mouse embryonic stem cells (ESCs) with an inducible knockout of Oct1, we showed that Oct1 deficiency resulted in poor induction of mesoderm-specific genes, leading to impaired mesodermal and terminal muscle differentiation. Oct1-deficient cells exhibited poor temporal coordination of the induction of lineage-specific genes and showed inappropriate developmental lineage branching, resulting in poorly differentiated cell states retaining epithelial characteristics. In ESCs, Oct1 localized with the pluripotency factor Oct4 at mesoderm-associated genes and remained bound to those loci during differentiation after the dissociation of Oct4. Binding events for Oct1 overlapped with those for the histone lysine demethylase Utx, and an interaction between Oct1 and Utx suggested that these two proteins cooperate to activate gene expression. The specificity of the ubiquitous Oct1 for the induction of mesodermal genes could be partially explained by the frequent coexistence of Smad and Oct binding sites at mesoderm-specific genes and the cooperative stimulation of mesodermal gene transcription by Oct1 and Smad3. Together, these results identify Oct1 as a key mediator of mesoderm lineage-specific gene induction.

A single-cell transcriptional timelapse of mouse embryonic development, from gastrula to pup

The house mouse, Mus musculus, is an exceptional model system, combining genetic tractability with close homology to human biology. Gestation in mouse development lasts just under three weeks, a period during which its genome orchestrates the astonishing transformation of a single cell zygote into a free-living pup composed of >500 million cells. Towards a global framework for exploring mammalian development, we applied single cell combinatorial indexing (sci-*) to profile the transcriptional states of 12.4 million nuclei from 83 precisely staged embryos spanning late gastrulation (embryonic day 8 or E8) to birth (postnatal day 0 or P0), with 2-hr temporal resolution during somitogenesis, 6-hr resolution through to birth, and 20-min resolution during the immediate postpartum period. From these data (E8 to P0), we annotate dozens of trajectories and hundreds of cell types and perform deeper analyses of the unfolding of the posterior embryo during somitogenesis as well as the ontogenesis of the kidney, mesenchyme, retina, and early neurons. Finally, we leverage the depth and temporal resolution of these whole embryo snapshots, together with other published data, to construct and curate a rooted tree of cell type relationships that spans mouse development from zygote to pup. Throughout this tree, we systematically nominate sets of transcription factors (TFs) and other genes as candidate drivers of the in vivo differentiation of hundreds of mammalian cell types. Remarkably, the most dramatic shifts in transcriptional state are observed in a restricted set of cell types in the hours immediately following birth, and presumably underlie the massive changes in physiology that must accompany the successful transition of a placental mammal to extrauterine life.

An interview with Valerie Wilson

Valerie Wilson is Professor of Early Embryo Development and a group leader at the Centre for Regenerative Medicine at the University of Edinburgh, UK. Val's lab works on the development of the vertebrate head-to-tail axis, focusing on the differentiation potential of axial progenitors. In 2022, she was awarded the British Society for Developmental Biology's (BSDB) Waddington Medal, which recognises outstanding individuals who have made major contributions to UK developmental biology. After her award lecture, we spoke to Val about her career and heard her thoughts on the past, present and future of the field.

Quantitative Experimental Embryology: A Modern Classical Approach

Experimental Embryology is often referred to as a classical approach of developmental biology that has been to some extent replaced by the introduction of molecular biology and genetic techniques to the field. Inspired by the combination of this approach with advanced techniques to uncover core principles of neural crest development by the laboratory of Roberto Mayor, we review key quantitative examples of experimental embryology from recent work in a broad range of developmental biology questions. We propose that quantitative experimental embryology offers essential ways to explore the reaction of cells and tissues to targeted cell addition, removal, and confinement. In doing so, it is an essential methodology to uncover principles of development that remain elusive such as pattern regulation, scaling, and self-organisation.

Single-cell multi-omics profiling links dynamic DNA methylation to cell fate decisions during mouse early organogenesis

BACKGROUND

Perturbation of DNA methyltransferases (DNMTs) and of the active DNA demethylation pathway via ten-eleven translocation (TET) methylcytosine dioxygenases results in severe developmental defects and embryonic lethality. Dynamic control of DNA methylation is therefore vital for embryogenesis, yet the underlying mechanisms remain poorly understood.

RESULTS

Here we report a single-cell transcriptomic atlas from Dnmt and Tet mutant mouse embryos during early organogenesis. We show that both the maintenance and de novo methyltransferase enzymes are dispensable for the formation of all major cell types at E8.5. However, DNA methyltransferases are required for silencing of prior or alternative cell fates such as pluripotency and extraembryonic programmes. Deletion of all three TET enzymes produces substantial lineage biases, in particular, a failure to generate primitive erythrocytes. Single-cell multi-omics profiling moreover reveals that this is linked to a failure to demethylate distal regulatory elements in Tet triple-knockout embryos.

CONCLUSIONS

This study provides a detailed analysis of the effects of perturbing DNA methylation on mouse organogenesis at a whole organism scale and affords new insights into the regulatory mechanisms of cell fate decisions.

Specification of the haematopoietic stem cell lineage: From blood-fated mesodermal angioblasts to haemogenic endothelium

Haematopoietic stem and progenitor cells emerge from specialized haemogenic endothelial cells in select vascular beds during embryonic development. Specification and commitment to the blood lineage, however, occur before endothelial cells are endowed with haemogenic competence, at the time of mesoderm patterning and production of endothelial cell progenitors (angioblasts). Whilst early blood cell fate specification has long been recognized, very little is known about the mechanisms that induce endothelial cell diversification and progressive acquisition of a blood identity by a subset of these cells. Here, we review the endothelial origin of the haematopoietic system and the complex developmental journey of blood-fated angioblasts. We discuss how recent technological advances will be instrumental to examine the diversity of the embryonic anatomical niches, signaling pathways and downstream epigenetic and transcriptional processes controlling endothelial cell heterogeneity and blood cell fate specification. Ultimately, this will give essential insights into the ontogeny of the cells giving rise to haematopoietic stem cells, that may aid in the development of novel strategies for their in vitro production for clinical purposes.

Mesp1 controls the chromatin and enhancer landscapes essential for spatiotemporal patterning of early cardiovascular progenitors

The mammalian heart arises from various populations of Mesp1-expressing cardiovascular progenitors (CPs) that are specified during the early stages of gastrulation. Mesp1 is a transcription factor that acts as a master regulator of CP specification and differentiation. However, how Mesp1 regulates the chromatin landscape of nascent mesodermal cells to define the temporal and spatial patterning of the distinct populations of CPs remains unknown. Here, by combining ChIP-seq, RNA-seq and ATAC-seq during mouse pluripotent stem cell differentiation, we defined the dynamic remodelling of the chromatin landscape mediated by Mesp1. We identified different enhancers that are temporally regulated to erase the pluripotent state and specify the pools of CPs that mediate heart development. We identified Zic2 and Zic3 as essential cofactors that act with Mesp1 to regulate its transcription-factor activity at key mesodermal enhancers, thereby regulating the chromatin remodelling and gene expression associated with the specification of the different populations of CPs in vivo. Our study identifies the dynamics of the chromatin landscape and enhancer remodelling associated with temporal patterning of early mesodermal cells into the distinct populations of CPs that mediate heart development.

Mesoderm induction and patterning: Insights from neuromesodermal progenitors

The discovery of mesoderm inducing signals helped usher in the era of molecular developmental biology, and today the mechanisms of mesoderm induction and patterning are still intensely studied. Mesoderm induction begins during gastrulation, but recent evidence in vertebrates shows that this process continues after gastrulation in a group of posteriorly localized cells called neuromesodermal progenitors (NMPs). NMPs reside within the post-gastrulation embryonic structure called the tailbud, where they make a lineage decision between ectoderm (spinal cord) and mesoderm. The majority of NMP-derived mesoderm generates somites, but also contributes to lateral mesoderm fates such as endothelium. The discovery of NMPs provides a new paradigm in which to study vertebrate mesoderm induction. This review will discuss mechanisms of mesoderm induction within NMPs, and how they have informed our understanding of mesoderm induction more broadly within vertebrates as well as animal species outside of the vertebrate lineage. Special focus will be given to the signaling networks underlying NMP-derived mesoderm induction and patterning, as well as emerging work on the significance of partial epithelial-mesenchymal states in coordinating cell fate and morphogenesis.

Generation and characterization of hair-bearing skin organoids from human pluripotent stem cells

Human skin uses millions of hairs and glands distributed across the body surface to function as an external barrier, thermoregulator and stimuli sensor. The large-scale generation of human skin with these appendages would be beneficial, but is challenging. Here, we describe a detailed protocol for generating hair-bearing skin tissue entirely from a homogeneous population of human pluripotent stem cells in a three-dimensional in vitro culture system. Defined culture conditions are used over a 2-week period to induce differentiation of pluripotent stem cells to surface ectoderm and cranial neural crest cells, which give rise to the epidermis and dermis, respectively, in each organoid unit. After 60 d of incubation, the skin organoids produce hair follicles. By day ~130, the skin organoids reach full complexity and contain stratified skin layers, pigmented hair follicles, sebaceous glands, Merkel cells and sensory neurons, recapitulating the cell composition and architecture of fetal skin tissue at week 18 of gestation. Skin organoids can be maintained in culture using this protocol for up to 150 d, enabling the organoids to be used to investigate basic skin biology, model disease and, further, reconstruct or regenerate skin tissue.

Spatial transcriptomic profiles of mouse uterine microenvironments at pregnancy day 7.5†

Uterine dysfunctions lead to fertility disorders and pregnancy complications. Normal uterine functions at pregnancy depend on crosstalk among multiple cell types in uterine microenvironments. Here, we performed the spatial transcriptomics and single-cell RNA-seq assays to determine local gene expression profiles at the embryo implantation site of the mouse uterus on pregnancy day 7.5 (D7.5). The spatial transcriptomic annotation identified 11 domains of distinct gene signatures, including a mesometrial myometrium, an anti-mesometrial myometrium, a mesometrial decidua enriched with natural killer cells, a vascular sinus zone for maternal vessel remodeling, a fetal-maternal interface, a primary decidual zone, a transition decidual zone, a secondary decidual zone, undifferentiated stroma, uterine glands, and the embryo. The scRNA-Seq identified 12 types of cells in the D7.5 uterus including three types of stromal fibroblasts with differentiated and undifferentiated markers, one cluster of epithelium including luminal and glandular epithelium, mesothelium, endothelia, pericytes, myelomonocytic cell, natural killer cells, and lymphocyte B. These single-cell RNA signatures were then utilized to deconvolute the cell-type compositions of each individual uterine microenvironment. Functional annotation assays on spatial transcriptomic data revealed uterine microenvironments with distinguished metabolic preferences, immune responses, and various cellular behaviors that are regulated by region-specific endocrine and paracrine signals. Global interactome among regions is also projected based on the spatial transcriptomic data. This study provides high-resolution transcriptome profiles with locality information at the embryo implantation site to facilitate further investigations on molecular mechanisms for normal pregnancy progression.

Regulation of Myogenesis by a Na/K-ATPase α1 Caveolin-Binding Motif

The N-terminal caveolin-binding motif (CBM) in Na/K-ATPase (NKA) α1 subunit is essential for cell signaling and somitogenesis in animals. To further investigate the molecular mechanism, we have generated CBM mutant human-induced pluripotent stem cells (iPSCs) through CRISPR/Cas9 genome editing and examined their ability to differentiate into skeletal muscle (Skm) cells. Compared with the parental wild-type human iPSCs, the CBM mutant cells lost their ability of Skm differentiation, which was evidenced by the absence of spontaneous cell contraction, marker gene expression, and subcellular myofiber banding structures in the final differentiated induced Skm cells. Another NKA functional mutant, A420P, which lacks NKA/Src signaling function, did not produce a similar defect. Indeed, A420P mutant iPSCs retained intact pluripotency and ability of Skm differentiation. Mechanistically, the myogenic transcription factor MYOD was greatly suppressed by the CBM mutation. Overexpression of a mouse Myod cDNA through lentiviral delivery restored the CBM mutant cells' ability to differentiate into Skm. Upstream of MYOD, Wnt signaling was demonstrated from the TOPFlash assay to have a similar inhibition. This effect on Wnt activity was further confirmed functionally by defective induction of the presomitic mesoderm marker genes BRACHYURY (T) and MESOGENIN1 (MSGN1) by Wnt3a ligand or the GSK3 inhibitor/Wnt pathway activator CHIR. Further investigation through immunofluorescence imaging and cell fractionation revealed a shifted membrane localization of β-catenin in CBM mutant iPSCs, revealing a novel molecular component of NKA-Wnt regulation. This study sheds light on a genetic regulation of myogenesis through the CBM of NKA and control of Wnt/β-catenin signaling.

Epha1 is a cell-surface marker for the neuromesodermal competent population

The vertebrate body is built during embryonic development by the sequential addition of new tissue as the embryo grows at its caudal end. During this process, progenitor cells within the neuromesodermal competent (NMC) region generate the postcranial neural tube and paraxial mesoderm. Here, we have applied a genetic strategy to recover the NMC cell population from mouse embryonic tissues and have searched their transcriptome for cell-surface markers that would give access to these cells without previous genetic modifications. We found that Epha1 expression is restricted to the axial progenitor-containing areas of the mouse embryo. Epha1-positive cells isolated from the mouse tailbud generate neural and mesodermal derivatives when cultured in vitro. This observation, together with their enrichment in the Sox2+/Tbxt+ molecular phenotype, indicates a direct association between Epha1 and the NMC population. Additional analyses suggest that tailbud cells expressing low Epha1 levels might also contain notochord progenitors, and that high Epha1 expression might be associated with progenitors entering paraxial mesoderm differentiation. Epha1 could thus be a valuable cell-surface marker for labeling and recovering physiologically active axial progenitors from embryonic tissues.

Integration of spatial and single-cell transcriptomic data elucidates mouse organogenesis

Molecular profiling of single cells has advanced our knowledge of the molecular basis of development. However, current approaches mostly rely on dissociating cells from tissues, thereby losing the crucial spatial context of regulatory processes. Here, we apply an image-based single-cell transcriptomics method, sequential fluorescence in situ hybridization (seqFISH), to detect mRNAs for 387 target genes in tissue sections of mouse embryos at the 8-12 somite stage. By integrating spatial context and multiplexed transcriptional measurements with two single-cell transcriptome atlases, we characterize cell types across the embryo and demonstrate that spatially resolved expression of genes not profiled by seqFISH can be imputed. We use this high-resolution spatial map to characterize fundamental steps in the patterning of the midbrain-hindbrain boundary (MHB) and the developing gut tube. We uncover axes of cell differentiation that are not apparent from single-cell RNA-sequencing (scRNA-seq) data, such as early dorsal-ventral separation of esophageal and tracheal progenitor populations in the gut tube. Our method provides an approach for studying cell fate decisions in complex tissues and development.

Diversification of muscle types in Drosophila embryos

Drosophila embryonic somatic muscles represent a simple and tractable model system to study the gene regulatory networks that control diversification of cell types. Somatic myogenesis in Drosophila is initiated by intrinsic action of the mesodermal master gene twist, which activates a cascade of transcriptional outputs including myogenic differentiation factor Mef2, which triggers all aspects of the myogenic differentiation program. In parallel, the expression of a combinatorial code of identity transcription factors (iTFs) defines discrete particular features of each muscle fiber, such as number of fusion events, and specific attachment to tendon cells or innervation, thus ensuring diversification of muscle types. Here, we take the example of a subset of lateral transverse (LT) muscles and discuss how the iTF code and downstream effector genes progressively define individual LT properties such as fusion program, attachment and innervation. We discuss new challenges in the field including the contribution of posttranscriptional and epitranscriptomic regulation of gene expression in the diversification of cell types.

Building consensus in neuromesodermal research: Current advances and future biomedical perspectives

The development of the vertebrate body axis relies on the activity of different populations of axial progenitors, including neuromesodermal progenitors. Currently, the term 'Neuromesodermal progenitors' is associated with various definitions. Here, we use distinct terminologies to highlight advances in our understanding of this cell type at both the single-cell and population levels. We discuss how these recent insights prompt new opportunities to address a range of biomedical questions spanning cancer metastasis, congenital disorders, cellular metabolism, regenerative medicine, and evolution. Finally, we outline some of the major unanswered questions and propose future directions at the forefront of neuromesodermal research.

Somite development and regionalisation of the vertebral axial skeleton

A critical stage in the development of all vertebrate embryos is the generation of the body plan and its subsequent patterning and regionalisation along the main anterior-posterior axis. This includes the formation of the vertebral axial skeleton. Its organisation begins during early embryonic development with the periodic formation of paired blocks of mesoderm tissue called somites. Here, we review axial patterning of somites, with a focus on studies using amniote model systems - avian and mouse. We summarise the molecular and cellular mechanisms that generate paraxial mesoderm and review how the different anatomical regions of the vertebral column acquire their specific identity and thus shape the body plan. We also discuss the generation of organoids and embryo-like structures from embryonic stem cells, which provide insights regarding axis formation and promise to be useful for disease modelling.

Cell-to-cell heterogeneity in Sox2 and Bra expression guides progenitor motility and destiny

Although cell-to-cell heterogeneity in gene and protein expression within cell populations has been widely documented, we know little about its biological functions. By studying progenitors of the posterior region of bird embryos, we found that expression levels of transcription factors Sox2 and Bra, respectively involved in neural tube (NT) and mesoderm specification, display a high degree of cell-to-cell heterogeneity. By combining forced expression and downregulation approaches with time-lapse imaging, we demonstrate that Sox2-to-Bra ratio guides progenitor's motility and their ability to stay in or exit the progenitor zone to integrate neural or mesodermal tissues. Indeed, high Bra levels confer high motility that pushes cells to join the paraxial mesoderm, while high levels of Sox2 tend to inhibit cell movement forcing cells to integrate the NT. Mathematical modeling captures the importance of cell motility regulation in this process and further suggests that randomness in Sox2/Bra cell-to-cell distribution favors cell rearrangements and tissue shape conservation.

Dynamics of primitive streak regression controls the fate of neuromesodermal progenitors in the chicken embryo

In classical descriptions of vertebrate development, the segregation of the three embryonic germ layers completes by the end of gastrulation. Body formation then proceeds in a head to tail fashion by progressive deposition of lineage-committed progenitors during regression of the primitive streak (PS) and tail bud (TB). The identification by retrospective clonal analysis of a population of neuromesodermal progenitors (NMPs) contributing to both musculoskeletal precursors (paraxial mesoderm) and spinal cord during axis formation challenged these notions. However, classical fate mapping studies of the PS region in amniotes have so far failed to provide direct evidence for such bipotential cells at the single-cell level. Here, using lineage tracing and single-cell RNA sequencing in the chicken embryo, we identify a resident cell population of the anterior PS epiblast, which contributes to neural and mesodermal lineages in trunk and tail. These cells initially behave as monopotent progenitors as classically described and only acquire a bipotential fate later, in more posterior regions. We show that NMPs exhibit a conserved transcriptomic signature during axis elongation but lose their epithelial characteristicsin the TB. Posterior to anterior gradients of convergence speed and ingression along the PS lead to asymmetric exhaustion of PS mesodermal precursor territories. Through limited ingression and increased proliferation, NMPs are maintained and amplified as a cell population which constitute the main progenitors in the TB. Together, our studies provide a novel understanding of the PS and TB contribution through the NMPs to the formation of the body of amniote embryos.

Coordinated changes in gene expression kinetics underlie both mouse and human erythroid maturation

BACKGROUND

Single-cell technologies are transforming biomedical research, including the recent demonstration that unspliced pre-mRNA present in single-cell RNA-Seq permits prediction of future expression states. Here we apply this RNA velocity concept to an extended timecourse dataset covering mouse gastrulation and early organogenesis.

RESULTS

Intriguingly, RNA velocity correctly identifies epiblast cells as the starting point, but several trajectory predictions at later stages are inconsistent with both real-time ordering and existing knowledge. The most striking discrepancy concerns red blood cell maturation, with velocity-inferred trajectories opposing the true differentiation path. Investigating the underlying causes reveals a group of genes with a coordinated step-change in transcription, thus violating the assumptions behind current velocity analysis suites, which do not accommodate time-dependent changes in expression dynamics. Using scRNA-Seq analysis of chimeric mouse embryos lacking the major erythroid regulator Gata1, we show that genes with the step-changes in expression dynamics during erythroid differentiation fail to be upregulated in the mutant cells, thus underscoring the coordination of modulating transcription rate along a differentiation trajectory. In addition to the expected block in erythroid maturation, the Gata1-chimera dataset reveals induction of PU.1 and expansion of megakaryocyte progenitors. Finally, we show that erythropoiesis in human fetal liver is similarly characterized by a coordinated step-change in gene expression.

CONCLUSIONS

By identifying a limitation of the current velocity framework coupled with in vivo analysis of mutant cells, we reveal a coordinated step-change in gene expression kinetics during erythropoiesis, with likely implications for many other differentiation processes.

Axial elongation of caudalized human organoids mimics aspects of neural tube development

Axial elongation of the neural tube is crucial during mammalian embryogenesis for anterior-posterior body axis establishment and subsequent spinal cord development, but these processes cannot be interrogated directly in humans as they occur post-implantation. Here, we report an organoid model of neural tube extension derived from human pluripotent stem cell (hPSC) aggregates that have been caudalized with Wnt agonism, enabling them to recapitulate aspects of the morphological and temporal gene expression patterns of neural tube development. Elongating organoids consist largely of neuroepithelial compartments and contain TBXT+SOX2+ neuro-mesodermal progenitors in addition to PAX6+NES+ neural progenitors. A critical threshold of Wnt agonism stimulated singular axial extensions while maintaining multiple cell lineages, such that organoids displayed regionalized anterior-to-posterior HOX gene expression with hindbrain (HOXB1) regions spatially distinct from brachial (HOXC6) and thoracic (HOXB9) regions. CRISPR interference-mediated silencing of TBXT, a Wnt pathway target, increased neuroepithelial compartmentalization, abrogated HOX expression and disrupted uniaxial elongation. Together, these results demonstrate the potent capacity of caudalized hPSC organoids to undergo axial elongation in a manner that can be used to dissect the cellular organization and patterning decisions that dictate early human nervous system development.