Identification of an embryonic differentiation stage marked by Sox1 and FoxA2 co-expression using combined cell tracking and high dimensional protein imaging.

Pluripotent mouse embryonic stem cells (ESCs) can differentiate to all germ layers and serve as an in vitro model of embryonic development. To better understand the differentiation paths traversed by ESCs committing to different lineages, we track individual differentiating ESCs by timelapse imaging followed by multiplexed high-dimensional Imaging Mass Cytometry (IMC) protein quantification. This links continuous live single-cell molecular NANOG and cellular dynamics quantification over 5-6 generations to protein expression of 37 different molecular regulators in the same single cells at the observation endpoints. Using this unique data set including kinship history and live lineage marker detection, we show that NANOG downregulation occurs generations prior to, but is not sufficient for neuroectoderm marker Sox1 upregulation. We identify a developmental cell type co-expressing both the canonical Sox1 neuroectoderm and FoxA2 endoderm markers in vitro and confirm the presence of such a population in the post-implantation embryo. RNASeq reveals cells co-expressing SOX1 and FOXA2 to have a unique cell state characterized by expression of both endoderm as well as neuroectoderm genes suggesting lineage potential towards both germ layers.

PHD2 safeguards modest mesendoderm development.

PHD2 is essential in modulating HIF-1α levels upon oxygen fluctuations. Hypoxia, a hallmark of uterus, and HIF-1α have recently emerged as opposing regulators of mesendoderm specification, suggesting a role for PHD2 therein. We found that PHD2 expression initially covered the epiblast and gradually receded from the primitive streak, which was identical to hypoxia and exclusive to HIF-1α. The investigations performed in mESCs, embryoids, and mouse embryos together demonstrated that PHD2 negatively regulated mesendoderm specification. Single-cell RNA sequencing revealed that PHD2 governed the transition from epiblast to mesendoderm. The downstream effect of PHD2 relied on the HIF-1α regulated Wnt/ß-catenin pathway, while it was regulated upstream by miR-429. In summary, our research highlights PHD2's essential role in mesendoderm specification and its interactions with hypoxia and HIF-1α.

A simple method for gene expression in endo- and ectodermal cells in mouse embryos before neural tube closure.

The lack of a widely accessible method for expressing genes of interest in wild-type embryos is a fundamental obstacle to understanding genetic regulation during embryonic development. In particular, only a few methods are available for introducing gene expression vectors into cells prior to neural tube closure, which is a period of drastic development for many tissues. In this study, we present a simple technique for injecting vectors into the amniotic cavity and allowing them to reach the ectodermal cells and the epithelia of endodermal organs of mouse embryos at E8.0 via in utero injection, using only a widely used optical fiber with an illuminator. Using this technique, retroviruses can be introduced to facilitate the labeling of cells in various tissues, including the brain, spinal cord, epidermis, and digestive and respiratory organs. We also demonstrated in utero electroporation of plasmid DNA into E7.0 and E8.0 embryos. Taking advantage of this method, we reveal the association between Ldb1 and the activity of the Neurog2 transcription factor in the mouse neocortex. This technique can aid in analyzing the roles of genes of interest during endo- and ectodermal development prior to neural tube closure.

A pre-vertebrate endodermal origin of calcitonin-producing neuroendocrine cells.

Vertebrate calcitonin-producing cells (C-cells) are neuroendocrine cells that secrete the small peptide hormone calcitonin in response to elevated blood calcium levels. Whereas mouse C-cells reside within the thyroid gland and derive from pharyngeal endoderm, avian C-cells are located within ultimobranchial glands and have been reported to derive from the neural crest. We use a comparative cell lineage tracing approach in a range of vertebrate model systems to resolve the ancestral embryonic origin of vertebrate C-cells. We find, contrary to previous studies, that chick C-cells derive from pharyngeal endoderm, with neural crest-derived cells instead contributing to connective tissue intimately associated with C-cells in the ultimobranchial gland. This endodermal origin of C-cells is conserved in a ray-finned bony fish (zebrafish) and a cartilaginous fish (the little skate, Leucoraja erinacea). Furthermore, we discover putative C-cell homologs within the endodermally-derived pharyngeal epithelium of the ascidian Ciona intestinalis and the amphioxus Branchiostoma lanceolatum, two invertebrate chordates that lack neural crest cells. Our findings point to a conserved endodermal origin of C-cells across vertebrates and to a pre-vertebrate origin of this cell type along the chordate stem.

Gata6 functions in zebrafish endoderm to regulate late differentiating arterial pole cardiogenesis.

Mutations in GATA6 are associated with congenital heart disease, most notably conotruncal structural defects. However, how GATA6 regulates cardiac morphology during embryogenesis is undefined. We used knockout and conditional mutant zebrafish alleles to investigate the spatiotemporal role of gata6 during cardiogenesis. Loss of gata6 specifically impacts atrioventricular valve formation and recruitment of epicardium, with a prominent loss of arterial pole cardiac cells, including those of the ventricle and outflow tract. However, there are no obvious defects in cardiac progenitor cell specification, proliferation or death. Conditional loss of gata6 starting at 24 h is sufficient to disrupt the addition of late differentiating cardiomyocytes at the arterial pole, with decreased expression levels of anterior secondary heart field (SHF) markers spry4 and mef2cb. Conditional loss of gata6 in the endoderm is sufficient to phenocopy the straight knockout, resulting in a significant loss of ventricular and outflow tract tissue. Exposure to a Dusp6 inhibitor largely rescues the loss of ventricular cells in gata6-/- larvae. Thus, gata6 functions in endoderm are mediated by FGF signaling to regulate the addition of anterior SHF progenitor derivatives during heart formation.

Medium acidosis drives cardiac differentiation during mesendoderm cell fate specification from human pluripotent stem cells.

Effective lineage-specific differentiation is essential to fulfilling the great potentials of human pluripotent stem cells (hPSCs). In this report, we investigate how modulation of medium pH and associated metabolic changes influence mesendoderm differentiation from hPSCs. We show that daily medium pH fluctuations are critical for the heterogeneity of cell fates in the absence of exogenous inducers. Acidic environment alone leads to cardiomyocyte generation without other signaling modulators. In contrast, medium alkalinization is inhibitory to cardiac fate even in the presence of classic cardiac inducers. We then demonstrate that acidic environment suppresses glycolysis to facilitate cardiac differentiation, while alkaline condition promotes glycolysis and diverts the differentiation toward other cell types. We further show that glycolysis inhibition or AMPK activation can rescue cardiac differentiation under alkalinization, and glycolysis inhibition alone can drive cardiac cell fate. This study highlights that pH changes remodel metabolic patterns and modulate signaling pathways to control cell fate.

Transcriptional and epigenetic characterization of a new in vitro platform to model the formation of human pharyngeal endoderm.

BACKGROUND: The Pharyngeal Endoderm (PE) is an extremely relevant developmental tissue, serving as the progenitor for the esophagus, parathyroids, thyroids, lungs, and thymus. While several studies have highlighted the importance of PE cells, a detailed transcriptional and epigenetic characterization of this important developmental stage is still missing, especially in humans, due to technical and ethical constraints pertaining to its early formation. RESULTS: Here we fill this knowledge gap by developing an in vitro protocol for the derivation of PE-like cells from human Embryonic Stem Cells (hESCs) and by providing an integrated multi-omics characterization. Our PE-like cells robustly express PE markers and are transcriptionally homogenous and similar to in vivo mouse PE cells. In addition, we define their epigenetic landscape and dynamic changes in response to Retinoic Acid by combining ATAC-Seq and ChIP-Seq of histone modifications. The integration of multiple high-throughput datasets leads to the identification of new putative regulatory regions and to the inference of a Retinoic Acid-centered transcription factor network orchestrating the development of PE-like cells. CONCLUSIONS: By combining hESCs differentiation with computational genomics, our work reveals the epigenetic dynamics that occur during human PE differentiation, providing a solid resource and foundation for research focused on the development of PE derivatives and the modeling of their developmental defects in genetic syndromes.

Cell size regulates human endoderm specification through actomyosin-dependent AMOT-YAP signaling.

Cell size is a crucial physical property that significantly impacts cellular physiology and function. However, the influence of cell size on stem cell specification remains largely unknown. Here, we investigated the dynamic changes in cell size during the differentiation of human pluripotent stem cells into definitive endoderm (DE). Interestingly, cell size exhibited a gradual decrease as DE differentiation progressed with higher stiffness. Furthermore, the application of hypertonic pressure or chemical to accelerate the reduction in cell size significantly and specifically enhanced DE differentiation. By functionally intervening in mechanosensitive elements, we have identified actomyosin activity as a crucial mediator of both DE differentiation and cell size reduction. Mechanistically, the reduction in cell size induces actomyosin-dependent angiomotin (AMOT) nuclear translocation, which suppresses Yes-associated protein (YAP) activity and thus facilitates DE differentiation. Together, our study has established a novel connection between cell size diminution and DE differentiation, which is mediated by AMOT nuclear translocation. Additionally, our findings suggest that the application of osmotic pressure can effectively promote human endodermal lineage differentiation.

A Three-Step Protocol to Differentiate iPSC into Colon Organoids.

Organoids, three-dimensional, stem cell-based structures that mimic the cellular and functional architecture of tissues, have emerged as an innovative in vitro tool. They offer highly efficient models for studying both embryonic development and disease progression processes. Colon organoids can also be generated from biopsies obtained during a colonoscopy. However, the invasive nature of biopsy collection poses practical challenges and introduces biases when studying patients who are already afflicted. Therefore, the use of iPSC-derived colon organoids can be considered a more practical approach for researchers and patients alike. Numerous protocols have been published for generating colon organoids from iPSCs. While most of these protocols share a common developmental process, some are labor-intensive or require additional equipment. Taking these considerations into account, we present a cost-effective and straightforward yet functionally robust colon organoid protocol: (1) definitive endoderm differentiation, (2) hindgut endoderm differentiation, and (3) maturation of colon spheroids into mature organoids.

Cell fate decision by a morphogen-transcription factor-chromatin modifier axis.

Cell fate decisions remain poorly understood at the molecular level. Embryogenesis provides a unique opportunity to analyze molecular details associated with cell fate decisions. Works based on model organisms have provided a conceptual framework of genes that specify cell fate control, for example, transcription factors (TFs) controlling processes from pluripotency to immunity1. How TFs specify cell fate remains poorly understood. Here we report that SALL4 relies on NuRD (nucleosome-remodeling and deacetylase complex) to interpret BMP4 signal and decide cell fate in a well-controlled in vitro system. While NuRD complex cooperates with SALL4 to convert mouse embryonic fibroblasts or MEFs to pluripotency, BMP4 diverts the same process to an alternative fate, PrE (primitive endoderm). Mechanistically, BMP4 signals the dissociation of SALL4 from NuRD physically to establish a gene regulatory network for PrE. Our results provide a conceptual framework to explore the rich landscapes of cell fate choices intrinsic to development in higher organisms involving morphogen-TF-chromatin modifier pathways.

Directed differentiation of pancreatic δ cells from human pluripotent stem cells.

Dysfunction of pancreatic δ cells contributes to the etiology of diabetes. Despite their important role, human δ cells are scarce, limiting physiological studies and drug discovery targeting δ cells. To date, no directed δ-cell differentiation method has been established. Here, we demonstrate that fibroblast growth factor (FGF) 7 promotes pancreatic endoderm/progenitor differentiation, whereas FGF2 biases cells towards the pancreatic δ-cell lineage via FGF receptor 1. We develop a differentiation method to generate δ cells from human stem cells by combining FGF2 with FGF7, which synergistically directs pancreatic lineage differentiation and modulates the expression of transcription factors and SST activators during endoderm/endocrine precursor induction. These δ cells display mature RNA profiles and fine secretory granules, secrete somatostatin in response to various stimuli, and suppress insulin secretion from in vitro co-cultured ß cells and mouse ß cells upon transplantation. The generation of human pancreatic δ cells from stem cells in vitro would provide an unprecedented cell source for drug discovery and cell transplantation studies in diabetes.

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.

Tissue-intrinsic beta-catenin signals antagonize Nodal-driven anterior visceral endoderm differentiation.

The anterior-posterior axis of the mammalian embryo is laid down by the anterior visceral endoderm (AVE), an extraembryonic signaling center that is specified within the visceral endoderm. Current models posit that AVE differentiation is promoted globally by epiblast-derived Nodal signals, and spatially restricted by a BMP gradient established by the extraembryonic ectoderm. Here, we report spatially restricted AVE differentiation in bilayered embryo-like aggregates made from mouse embryonic stem cells that lack an extraembryonic ectoderm. Notably, clusters of AVE cells also form in pure visceral endoderm cultures upon activation of Nodal signaling, indicating that tissue-intrinsic factors can restrict AVE differentiation. We identify ß-catenin activity as a tissue-intrinsic factor that antagonizes AVE-inducing Nodal signals. Together, our results show how an AVE-like population can arise through interactions between epiblast and visceral endoderm alone. This mechanism may be a flexible solution for axis patterning in a wide range of embryo geometries, and provide robustness to axis patterning when coupled with signal gradients.

A robust model for cell type-specific interindividual variation in single-cell RNA sequencing data.

Single-cell RNA sequencing (scRNA-seq) has been widely used to characterize cell types based on their average gene expression profiles. However, most studies do not consider cell type-specific variation across donors. Modelling this cell type-specific inter-individual variation could help elucidate cell type-specific biology and inform genes and cell types underlying complex traits. We therefore develop a new model to detect and quantify cell type-specific variation across individuals called CTMM (Cell Type-specific linear Mixed Model). We use extensive simulations to show that CTMM is powerful and unbiased in realistic settings. We also derive calibrated tests for cell type-specific interindividual variation, which is challenging given the modest sample sizes in scRNA-seq. We apply CTMM to scRNA-seq data from human induced pluripotent stem cells to characterize the transcriptomic variation across donors as cells differentiate into endoderm. We find that almost 100% of transcriptome-wide variability between donors is differentiation stage-specific. CTMM also identifies individual genes with statistically significant stage-specific variability across samples, including 85 genes that do not have significant stage-specific mean expression. Finally, we extend CTMM to partition interindividual covariance between stages, which recapitulates the overall differentiation trajectory. Overall, CTMM is a powerful tool to illuminate cell type-specific biology in scRNA-seq.

Extraembryonic gut endoderm cells undergo programmed cell death during development.

Despite a distinct developmental origin, extraembryonic cells in mice contribute to gut endoderm and converge to transcriptionally resemble their embryonic counterparts. Notably, all extraembryonic progenitors share a non-canonical epigenome, raising several pertinent questions, including whether this landscape is reset to match the embryonic regulation and if extraembryonic cells persist into later development. Here we developed a two-colour lineage-tracing strategy to track and isolate extraembryonic cells over time. We find that extraembryonic gut cells display substantial memory of their developmental origin including retention of the original DNA methylation landscape and resulting transcriptional signatures. Furthermore, we show that extraembryonic gut cells undergo programmed cell death and neighbouring embryonic cells clear their remnants via non-professional phagocytosis. By midgestation, we no longer detect extraembryonic cells in the wild-type gut, whereas they persist and differentiate further in p53-mutant embryos. Our study provides key insights into the molecular and developmental fate of extraembryonic cells inside the embryo.

TGF-ß modulates cell fate in human ES cell-derived foregut endoderm by inhibiting Wnt and BMP signaling.

Genetic differences between pluripotent stem cell lines cause variable activity of extracellular signaling pathways, limiting reproducibility of directed differentiation protocols. Here we used human embryonic stem cells (hESCs) to interrogate how exogenous factors modulate endogenous signaling events during specification of foregut endoderm lineages. We find that transforming growth factor ß1 (TGF-ß1) activates a putative human OTX2/LHX1 gene regulatory network which promotes anterior fate by antagonizing endogenous Wnt signaling. In contrast to Porcupine inhibition, TGF-ß1 effects cannot be reversed by exogenous Wnt ligands, suggesting that induction of SHISA proteins and intracellular accumulation of Fzd receptors render TGF-ß1-treated cells refractory to Wnt signaling. Subsequently, TGF-ß1-mediated inhibition of BMP and Wnt signaling suppresses liver fate and promotes pancreas fate. Furthermore, combined TGF-ß1 treatment and Wnt inhibition during pancreatic specification reproducibly and robustly enhance INSULIN+ cell yield across hESC lines. This modification of widely used differentiation protocols will enhance pancreatic ß cell yield for cell-based therapeutic applications.

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

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

The primitive endoderm supports lineage plasticity to enable regulative development.

Mammalian blastocyst formation involves the specification of the trophectoderm followed by the differentiation of the inner cell mass into embryonic epiblast and extra-embryonic primitive endoderm (PrE). During this time, the embryo maintains a window of plasticity and can redirect its cellular fate when challenged experimentally. In this context, we found that the PrE alone was sufficient to regenerate a complete blastocyst and continue post-implantation development. We identify an in vitro population similar to the early PrE in vivo that exhibits the same embryonic and extra-embryonic potency and can form complete stem cell-based embryo models, termed blastoids. Commitment in the PrE is suppressed by JAK/STAT signaling, collaborating with OCT4 and the sustained expression of a subset of pluripotency-related transcription factors that safeguard an enhancer landscape permissive for multi-lineage differentiation. Our observations support the notion that transcription factor persistence underlies plasticity in regulative development and highlight the importance of the PrE in perturbed development.

Coordination between endoderm progression and mouse gastruloid elongation controls endodermal morphotype choice.

Embryonic development is highly robust. Morphogenetic variability between embryos (under ideal conditions) is largely quantitative. This robustness stands in contrast to in vitro embryo-like models, which, like most organoids, can display a high degree of tissue morphogenetic variability. The source of this difference is not fully understood. We use the mouse gastruloid model to study the morphogenetic progression of definitive endoderm (DE) and its divergence. We first catalog the different morphologies and characterize their statistics. We then learn predictive models for DE morphotype based on earlier expression and morphology measurements. Finally, we analyze these models to identify key drivers of morphotype variability and devise gastruloid-specific and global interventions that can lower this variability and steer morphotype choice. In the process, we identify two types of coordination lacking in the in vitro model but required for robust gut-tube formation. This approach can help improve the quality and usability of 3D embryo-like models.

An integrated approach identifies the molecular underpinnings of murine anterior visceral endoderm migration.

The anterior visceral endoderm (AVE) differs from the surrounding visceral endoderm (VE) in its migratory behavior and ability to restrict primitive streak formation to the opposite side of the mouse embryo. To characterize the molecular bases for the unique properties of the AVE, we combined single-cell RNA sequencing of the VE prior to and during AVE migration with phosphoproteomics, high-resolution live-imaging, and short-term lineage labeling and intervention. This identified the transient nature of the AVE with attenuation of "anteriorizing" gene expression as cells migrate and the emergence of heterogeneities in transcriptional states relative to the AVE's position. Using cell communication analysis, we identified the requirement of semaphorin signaling for normal AVE migration. Lattice light-sheet microscopy showed that Sema6D mutants have abnormalities in basal projections and migration speed. These findings point to a tight coupling between transcriptional state and position of the AVE and identify molecular controllers of AVE migration.