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.

Palatal segment contributions to midfacial growth in three inbred mouse strains.

Following facial prominence fusion, anterior-posterior (A-P) elongation of the palate is a critical aspect of palatogenesis and integrated midfacial elongation. Reciprocal epithelial-mesenchymal interactions drive secondary palate elongation and periodic signaling center formation within the rugae growth zone (RGZ). However, the relationship between RGZ dynamics and the morphogenetic behavior of underlying palatal bone mesenchymal precursors has remained enigmatic. Our results indicate that cellular activity at the RGZ simultaneously drives rugae formation and elongation of the maxillary bone primordium within the anterior secondary palate, which more than doubles in length prior to palatal shelf fusion. The first formed palatal ruga, found just posterior to the RGZ, represents a consistent morphological boundary between anterior and posterior secondary palate bone precursors, being found at the future maxillary-palatine suture. These results suggest a model where changes in RGZ-driven A-P growth of the anterior secondary palate may produce interspecies and intraspecies differences in facial prognathism and differences in the proportional contribution of palatal segment-associated bones to total palate length. An ontogenetic comparison of three inbred mouse strains indicated that while RGZ-driven growth of the anterior secondary palate is critical for early midfacial outgrowth, subtle strain-specific bony contributions to adult palate length are not present until after this initial palatal growth period. This multifaceted illustration of normal midfacial growth dynamics confirms a one-to-one relationship between palatal segments and upper jaw bones during the earliest stages of palatal growth, which may serve as the basis for evolutionary change in upper jaw morphology. Additionally, identified mouse strain-specific differences in palatal segment elongation provide a useful foundation for understanding the impact of background genetic effects on facial morphogenesis.

The spontaneous mouse mutant low set ears (Lse) is caused by tandem duplication of Fgf3 and Fgf4.

The external ear develops from an organized convergence of ventrally migrating neural crest cells into the first and second branchial arches. Defects in external ear position are often symptomatic of complex syndromes such as Apert, Treacher-Collins, and Crouzon Syndrome. The low set ears (Lse) spontaneous mouse mutant is characterized by the dominant inheritance of a ventrally shifted external ear position and an abnormal external auditory meatus (EAM). We identified the causative mutation as a 148 Kb tandem duplication on Chromosome 7, which includes the entire coding sequences of Fgf3 and Fgf4. Duplications of FGF3 and FGF4 occur in 11q duplication syndrome in humans and are associated with craniofacial anomalies, among other features. Intercrosses of Lse-affected mice revealed perinatal lethality in homozygotes, and Lse/Lse embryos display additional phenotypes including polydactyly, abnormal eye morphology, and cleft secondary palate. The duplication results in increased Fgf3 and Fgf4 expression in the branchial arches and additional discrete domains in the developing embryo. This ectopic overexpression resulted in functional FGF signaling, demonstrated by increased Spry2 and Etv5 expression in overlapping domains of the developing arches. Finally, a genetic interaction between Fgf3/4 overexpression and Twist1, a regulator of skull suture development, resulted in perinatal lethality, cleft palate, and polydactyly in compound heterozygotes. These data indicate a role for Fgf3 and Fgf4 in external ear and palate development and provide a novel mouse model for further interrogation of the biological consequences of human FGF3/4 duplication.

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.

Systematic reconstruction of cellular trajectories across mouse embryogenesis.

Mammalian embryogenesis is characterized by rapid cellular proliferation and diversification. Within a few weeks, a single-cell zygote gives rise to millions of cells expressing a panoply of molecular programs. Although intensively studied, a comprehensive delineation of the major cellular trajectories that comprise mammalian development in vivo remains elusive. Here, we set out to integrate several single-cell RNA-sequencing (scRNA-seq) datasets that collectively span mouse gastrulation and organogenesis, supplemented with new profiling of ~150,000 nuclei from approximately embryonic day 8.5 (E8.5) embryos staged in one-somite increments. Overall, we define cell states at each of 19 successive stages spanning E3.5 to E13.5 and heuristically connect them to their pseudoancestors and pseudodescendants. Although constructed through automated procedures, the resulting directed acyclic graph (TOME (trajectories of mammalian embryogenesis)) is largely consistent with our contemporary understanding of mammalian development. We leverage TOME to systematically nominate transcription factors (TFs) as candidate regulators of each cell type's specification, as well as 'cell-type homologs' across vertebrate evolution.

Loss of Extreme Long-Range Enhancers in Human Neural Crest Drives a Craniofacial Disorder.

Non-coding mutations at the far end of a large gene desert surrounding the SOX9 gene result in a human craniofacial disorder called Pierre Robin sequence (PRS). Leveraging a human stem cell differentiation model, we identify two clusters of enhancers within the PRS-associated region that regulate SOX9 expression during a restricted window of facial progenitor development at distances up to 1.45 Mb. Enhancers within the 1.45 Mb cluster exhibit highly synergistic activity that is dependent on the Coordinator motif. Using mouse models, we demonstrate that PRS phenotypic specificity arises from the convergence of two mechanisms: confinement of Sox9 dosage perturbation to developing facial structures through context-specific enhancer activity and heightened sensitivity of the lower jaw to Sox9 expression reduction. Overall, we characterize the longest-range human enhancers involved in congenital malformations, directly demonstrate that PRS is an enhanceropathy, and illustrate how small changes in gene expression can lead to morphological variation.

FaceBase 3: analytical tools and FAIR resources for craniofacial and dental research.

The FaceBase Consortium was established by the National Institute of Dental and Craniofacial Research in 2009 as a 'big data' resource for the craniofacial research community. Over the past decade, researchers have deposited hundreds of annotated and curated datasets on both normal and disordered craniofacial development in FaceBase, all freely available to the research community on the FaceBase Hub website. The Hub has developed numerous visualization and analysis tools designed to promote integration of multidisciplinary data while remaining dedicated to the FAIR principles of data management (findability, accessibility, interoperability and reusability) and providing a faceted search infrastructure for locating desired data efficiently. Summaries of the datasets generated by the FaceBase projects from 2014 to 2019 are provided here. FaceBase 3 now welcomes contributions of data on craniofacial and dental development in humans, model organisms and cell lines. Collectively, the FaceBase Consortium, along with other NIH-supported data resources, provide a continuously growing, dynamic and current resource for the scientific community while improving data reproducibility and fulfilling data sharing requirements.

Pbx loss in cranial neural crest, unlike in epithelium, results in cleft palate only and a broader midface.

Orofacial clefting represents the most common craniofacial birth defect. Cleft lip with or without cleft palate (CL/P) is genetically distinct from cleft palate only (CPO). Numerous transcription factors (TFs) regulate normal development of the midface, comprising the premaxilla, maxilla and palatine bones, through control of basic cellular behaviors. Within the Pbx family of genes encoding Three Amino-acid Loop Extension (TALE) homeodomain-containing TFs, we previously established that in the mouse, Pbx1 plays a preeminent role in midfacial morphogenesis, and Pbx2 and Pbx3 execute collaborative functions in domains of coexpression. We also reported that Pbx1 loss from cephalic epithelial domains, on a Pbx2- or Pbx3-deficient background, results in CL/P via disruption of a regulatory network that controls apoptosis at the seam of frontonasal and maxillary process fusion. Conversely, Pbx1 loss in cranial neural crest cell (CNCC)-derived mesenchyme on a Pbx2-deficient background results in CPO, a phenotype not yet characterized. In this study, we provide in-depth analysis of PBX1 and PBX2 protein localization from early stages of midfacial morphogenesis throughout development of the secondary palate. We further establish CNCC-specific roles of PBX TFs and describe the developmental abnormalities resulting from their loss in the murine embryonic secondary palate. Additionally, we compare and contrast the phenotypes arising from PBX1 loss in CNCC with those caused by its loss in the epithelium and show that CNCC-specific Pbx1 deletion affects only later secondary palate morphogenesis. Moreover, CNCC mutants exhibit perturbed rostro-caudal organization and broadening of the midfacial complex. Proliferation defects are pronounced in CNCC mutants at gestational day (E)12.5, suggesting altered proliferation of mutant palatal progenitor cells, consistent with roles of PBX factors in maintaining progenitor cell state. Although the craniofacial skeletal abnormalities in CNCC mutants do not result from overt patterning defects, osteogenesis is delayed, underscoring a critical role of PBX factors in CNCC morphogenesis and differentiation. Overall, the characterization of tissue-specific Pbx loss-of-function mouse models with orofacial clefting establishes these strains as unique tools to further dissect the complexities of this congenital craniofacial malformation. This study closely links PBX TALE homeodomain proteins to the variation in maxillary shape and size that occurs in pathological settings and during evolution of midfacial morphology.

Chromatin Architecture of the Pitx2 Locus Requires CTCF- and Pitx2-Dependent Asymmetry that Mirrors Embryonic Gut Laterality.

Expression of Pitx2 on the left side of the embryo patterns left-right (LR) organs including the dorsal mesentery (DM), whose asymmetric cell behavior directs gut looping. Despite the importance of organ laterality, chromatin-level regulation of Pitx2 remains undefined. Here, we show that genes immediately neighboring Pitx2 in chicken and mouse, including a long noncoding RNA (Pitx2 locus-asymmetric regulated RNA or Playrr), are expressed on the right side and repressed by Pitx2. CRISPR/Cas9 genome editing of Playrr, 3D fluorescent in situ hybridization (FISH), and variations of chromatin conformation capture (3C) demonstrate that mutual antagonism between Pitx2 and Playrr is coordinated by asymmetric chromatin interactions dependent on Pitx2 and CTCF. We demonstrate that transcriptional and morphological asymmetries driving gut looping are mirrored by chromatin architectural asymmetries at the Pitx2 locus. We propose a model whereby Pitx2 auto-regulation directs chromatin topology to coordinate LR transcription of this locus essential for LR organogenesis.

The left-right Pitx2 pathway drives organ-specific arterial and lymphatic development in the intestine.

The dorsal mesentery (DM) is the major conduit for blood and lymphatic vessels in the gut. The mechanisms underlying their morphogenesis are challenging to study and remain unknown. Here we show that arteriogenesis in the DM begins during gut rotation and proceeds strictly on the left side, dependent on the Pitx2 target gene Cxcl12. Although competent Cxcr4-positive angioblasts are present on the right, they fail to form vessels and progressively emigrate. Surprisingly, gut lymphatics also initiate in the left DM and arise only after-and dependent on-arteriogenesis, implicating arteries as drivers of gut lymphangiogenesis. Our data begin to unravel the origin of two distinct vascular systems and demonstrate how early left-right molecular asymmetries are translated into organ-specific vascular patterns. We propose a dual origin of gut lymphangiogenesis in which prior arterial growth is required to initiate local lymphatics that only subsequently connect to the vascular system.