ABSTRACT Single cell genomics has become a popular approach to uncover the cellular heterogeneity of progenitor and terminally differentiated cell types with great precision. This approach can also delineate lineage hierarchies and identify molecular programmes of cell-fate acquisition and segregation. Nowadays, tens of thousands of cells are routinely sequenced in single cell-based methods and even more are expected to be analysed in the future. However, interpretation of the resulting data is challenging and requires computational models at multiple levels of abstraction. In contrast to other applications of single cell sequencing, where clustering approaches dominate, developmental systems are generally modelled using continuous structures, trajectories and trees. These trajectory models carry the promise of elucidating mechanisms of development, disease and stimulation response at very high molecular resolution. However, their reliable analysis and biological interpretation requires an understanding of their underlying assumptions and limitations. Here, we review the basic concepts of such computational approaches and discuss the characteristics of developmental processes that can be learnt from trajectory models.

ABSTRACT Deciphering mechanisms of endocrine cell induction, specification and lineage allocation in vivo will provide valuable insights into how the islets of Langerhans are generated. Currently, it is ill defined how endocrine progenitors segregate into different endocrine subtypes during development. Here, we generated a novel neurogenin 3 (Ngn3)-Venus fusion (NVF) reporter mouse line, that closely mirrors the transient endogenous Ngn3 protein expression. To define an in vivo roadmap of endocrinogenesis, we performed single cell RNA sequencing of 36,351 pancreatic epithelial and NVF + cells during secondary transition. This allowed Ngn3 low endocrine progenitors, Ngn3 high endocrine precursors, Fev + endocrine lineage and hormone + endocrine subtypes to be distinguished and time-resolved, and molecular programs during the step-wise lineage restriction steps to be delineated. Strikingly, we identified 58 novel signature genes that show the same transient expression dynamics as Ngn3 in the 7260 profiled Ngn3 -expressing cells. The differential expression of these genes in endocrine precursors associated with their cell-fate allocation towards distinct endocrine cell types. Thus, the generation of an accurately regulated NVF reporter allowed us to temporally resolve endocrine lineage development to provide a fine-grained single cell molecular profile of endocrinogenesis in vivo .

ABSTRACT Mutations in Hes1 , a target gene of the Notch signalling pathway, lead to ectopic pancreas by a poorly described mechanism. Here, we use genetic inactivation of Hes1 combined with lineage tracing and live imaging to reveal an endodermal requirement for Hes1, and show that ectopic pancreas tissue is derived from the dorsal pancreas primordium. RNA-seq analysis of sorted E10.5 Hes1 +/+ and Hes1 −/− Pdx1-GFP + cells suggested that upregulation of endocrine lineage genes in Hes1 −/− embryos was the major defect and, accordingly, early pancreas morphogenesis was normalized, and the ectopic pancreas phenotype suppressed, in Hes1 −/− Neurog3 −/− embryos. In Mib1 mutants, we found a near total depletion of dorsal progenitors, which was replaced by an anterior Gcg + extension. Together, our results demonstrate that aberrant morphogenesis is the cause of ectopic pancreas and that a part of the endocrine differentiation program is mechanistically involved in the dysgenesis. Our results suggest that the ratio of endocrine lineage to progenitor cells is important for morphogenesis and that a strong endocrinogenic phenotype without complete progenitor depletion, as seen in Hes1 mutants, provokes an extreme dysgenesis that causes ectopic pancreas.

ABSTRACT The pancreas is an endoderm-derived glandular organ that participates in the regulation of systemic glucose metabolism and food digestion through the function of its endocrine and exocrine compartments, respectively. While intensive research has explored the signaling pathways and transcriptional programs that govern pancreas development, much remains to be discovered regarding the cellular processes that orchestrate pancreas morphogenesis. Here, we discuss the developmental mechanisms and principles that are known to underlie pancreas development, from induction and lineage formation to morphogenesis and organogenesis. Elucidating such principles will help to identify novel candidate disease genes and unravel the pathogenesis of pancreas-related diseases, such as diabetes, pancreatitis and cancer.

Current understanding infers a neural crest origin of thyroid C cells, the major source of calcitonin in mammals and ancestors to neuroendocrine thyroid tumors. The concept is primarily based on investigations in quail–chick chimeras involving fate mapping of neural crest cells to the ultimobranchial glands that regulate Ca 2+ homeostasis in birds, reptiles, amphibians and fishes, but whether mammalian C cell development involves a homologous ontogenetic trajectory has not been experimentally verified. With lineage tracing, we now provide direct evidence that Sox17+ anterior endoderm is the only source of differentiated C cells and their progenitors in mice. Like many gut endoderm derivatives, embryonic C cells were found to coexpress pioneer factors forkhead box (Fox) a1 and Foxa2 before neuroendocrine differentiation takes place. In the ultimobranchial body epithelium emerging from pharyngeal pouch endoderm in early organogenesis, differential Foxa1/Foxa2 expression distinguished two spatially separated pools of C cell precursors with different growth properties. A similar expression pattern was recapitulated in medullary thyroid carcinoma cells in vivo , consistent with a growth-promoting role of Foxa1. In contrast to embryonic precursor cells, C cell-derived tumor cells invading the stromal compartment downregulated Foxa2, foregoing epithelial-to-mesenchymal transition designated by loss of E-cadherin; both Foxa2 and E-cadherin were re-expressed at metastatic sites. These findings revise mammalian C cell ontogeny, expand the neuroendocrine repertoire of endoderm and redefine the boundaries of neural crest diversification. The data further underpin distinct functions of Foxa1 and Foxa2 in both embryonic and tumor development.

Transcription factors (TFs) pattern developing tissues and determine cell fates; however, how spatio-temporal TF gradients are generated is ill defined. Here we show that miR-335 fine-tunes TF gradients in the endoderm and promotes mesendodermal lineage segregation. Initially, we identified miR-335 as a regulated intronic miRNA in differentiating embryonic stem cells (ESCs). miR-335 is encoded in the mesoderm-specific transcript ( Mest ) and targets the 3′-UTRs of the endoderm-determining TFs Foxa2 and Sox17. Mest and miR-335 are co-expressed and highly accumulate in the mesoderm, but are transiently expressed in endoderm progenitors. Overexpression of miR-335 does not affect initial mesendoderm induction, but blocks Foxa2- and Sox17-mediated endoderm differentiation in ESCs and ESC-derived embryos. Conversely, inhibition of miR-335 activity leads to increased Foxa2 and Sox17 protein accumulation and endoderm formation. Mathematical modeling predicts that transient miR-335 expression in endoderm progenitors shapes a TF gradient in the endoderm, which we confirm by functional studies in vivo . Taken together, our results suggest that miR-335 targets endoderm TFs for spatio-temporal gradient formation in the endoderm and to stabilize lineage decisions during mesendoderm formation.

Several signalling cascades are implicated in the formation and patterning of the three principal germ layers, but their precise temporal-spatial mode of action in progenitor populations remains undefined. We have used conditional gene deletion of mouse β-catenin in Sox17-positive embryonic and extra-embryonic endoderm as well as vascular endothelial progenitors to address the function of canonical Wnt signalling in cell lineage formation and patterning. Conditional mutants fail to form anterior brain structures and exhibit posterior body axis truncations, whereas initial blood vessel formation appears normal. Tetraploid rescue experiments reveal that lack of β-catenin in the anterior visceral endoderm results in defects in head organizer formation. Sox17 lineage tracing in the definitive endoderm (DE) shows a cell-autonomous requirement for β-catenin in midgut and hindgut formation. Surprisingly, wild-type posterior visceral endoderm (PVE) in midgut- and hindgut-deficient tetraploid chimera rescues the posterior body axis truncation, indicating that the PVE is important for tail organizer formation. Upon loss of β-catenin in the visceral endoderm and DE lineages, but not in the vascular endothelial lineage, Sox17 expression is not maintained, suggesting downstream regulation by canonical Wnt signalling. Strikingly, Tcf4/β-catenin transactivation complexes accumulated on Sox17 cis-regulatory elements specifically upon endoderm induction in an embryonic stem cell differentiation system. Together, these results indicate that the Wnt/β-catenin signalling pathway regulates Sox17 expression for visceral endoderm pattering and DE formation and provide the first functional evidence that the PVE is necessary for gastrula organizer gene induction and posterior axis development.

The thymus and parathyroid glands are derived from the third pharyngeal pouch endoderm. The mechanisms that establish distinct molecular domains in the third pouch and control the subsequent separation of these organ primordia from the pharynx are poorly understood. Here, we report that mouse embryos that lack two FGF feedback antagonists, Spry1 and Spry2, display parathyroid and thymus hypoplasia and a failure of these organ primordia to completely separate from the pharynx. We show that FGF ligands and downstream reporter genes are expressed in highly regionalised patterns in the third pouch and that sprouty gene deletion results in upregulated FGF signalling throughout the pouch endoderm. As a consequence, the initiation of markers of parathyroid and thymus fate is altered. In addition, a normal apoptotic programme that is associated with the separation of the primordia from the pharynx is disrupted, resulting in the maintenance of a thymus-pharynx attachment and a subsequent inability of the thymus to migrate to its appropriate position above the heart. We demonstrate that the sprouty genes function in the pharyngeal endoderm itself to control these processes and that the defects in sprouty-deficient mutants are, at least in part, due to hyper-responsiveness to Fgf8. Finally, we provide evidence to suggest that parathyroid hypoplasia in these mutants is due to early gene expression defects in the third pouch, whereas thymus hypoplasia is caused by reduced proliferation of thymic epithelial cells in the thymus primordium.

In the mouse, one of the earliest events in the determination of cell fate is the segregation of cells into germ layers during gastrulation; however, the cellular and molecular details are not well defined due to intrauterine development. We were able to visualize a clear sequence of events occurring in the process of germ-layer formation, using immunohistochemistry and time-lapse confocal imaging. The T-box transcription factor brachyury (T) and the Forkhead transcription factor Foxa2 specify mesoderm and endoderm in the posterior epiblast. Fate-specified epiblast cells lose their polarity and undergo epithelial-mesenchymal transition to invade into the primitive streak region, where these cell populations quickly separate and differentiate into morphologically and molecularly distinct Foxa2-positive endoderm and T-positive mesoderm populations. The endoderm cells flatten and acquire apical-basal polarity during intercalation into the outside epithelium in order to establish proper intracellular junctions with pre-existing cells. By contrast, the mesodermal cells become spherical during migration and acquire a mesenchymal fate. Interestingly, axial mesodermal cells are descended from Foxa2-positive epiblast cells that upregulate T protein in the anterior primitive streak region. These cells, as well as Foxa2-positive endoderm cells, are highly polarized and epithelialized, suggesting that Foxa2 promotes an epithelial fate and suppresses a mesenchymal fate. This observation is supported by the fact that Foxa2 mutant endodermal cells fail to maintain polarity and do not establish proper cellular junctions, and are thus unable to functionally integrate into the endoderm epithelium. We propose that Foxa2 regulates a molecular program that induces an epithelial cellular phenotype.

Differential gene regulation integrated in time and space drives developmental programs during embryogenesis. To understand how the program of gastrulation is regulated by Wnt/β-catenin signaling, we have used genome-wide expression profiling of conditional β-catenin mutant embryos. Known Wnt/β-catenin target genes, known components of other signaling pathways, as well as a number of uncharacterized genes were downregulated in these mutants. To further narrow down the set of differentially expressed genes, we used whole-mount in situ screening to associate gene expression with putative domains of Wnt activity. Several potential novel target genes were identified by this means and two, Grsf1 and Fragilis2 , were functionally analyzed by RNA interference (RNAi) in completely embryonic stem(ES) cell-derived embryos. We show that the gene encoding the RNA-binding factor Grsf1 is important for axial elongation, mid/hindbrain development and axial mesoderm specification, and that Fragilis2 ,encoding a transmembrane protein, regulates epithelialization of the somites and paraxial mesoderm formation. Intriguingly, the knock-down phenotypes recapitulate several aspects of Wnt pathway mutants, suggesting that these genes are components of the downstream Wnt response. This functional genomic approach allows the rapid identification of functionally important components of embryonic development from large datasets of putative targets.