Single-cell analysis of bidirectional reprogramming between early embryonic states identify mechanisms of differential lineage plasticities in mice

Two distinct lineages, pluripotent epiblast (EPI) and primitive (extra-embryonic) endoderm (PrE), arise from common inner cell mass (ICM) progenitors in mammalian embryos. To study how these sister identities are forged, we leveraged mouse embryonic stem (ES) cells and extra-embryonic endoderm (XEN) stem cells-in vitro counterparts of the EPI and PrE. Bidirectional reprogramming between ES and XEN coupled with single-cell RNA and ATAC-seq analyses showed distinct rates, efficiencies, and trajectories of state conversions, identifying drivers and roadblocks of reciprocal conversions. While GATA4-mediated ES-to-iXEN conversion was rapid and nearly deterministic, OCT4-, KLF4-, and SOX2-induced XEN-to-induced pluripotent stem (iPS) reprogramming progressed with diminished efficiency and kinetics. A dominant PrE transcriptional program, safeguarded by GATA4, alongside elevated chromatin accessibility and reduced DNA methylation of the EPI underscored the differential plasticities of the two states. Mapping in vitro to embryo trajectories tracked reprogramming cells in either direction along EPI and PrE in vivo states, without transitioning through the ICM.

Role of tristability in the robustness of the differentiation mechanism

During cell differentiation, identical pluripotent cells undergo a specification process marked by changes in the expression of key genes, regulated by transcription factors that can inhibit the transcription of a competing gene or activate their own transcription. This specification is orchestrated by gene regulatory networks (GRNs), encompassing transcription factors, biochemical reactions, and signalling cascades. Mathematical models for these GRNs have been proposed in various contexts, to replicate observed robustness in differentiation properties. This includes reproducible proportions of differentiated cells with respect to parametric or stochastic noise and the avoidance of transitions between differentiated states. Understanding the GRN components controlling these features is crucial. Our study thoroughly explored an extended version of the Toggle Switch model with auto-activation loops. This model represents cells evolving from common progenitors in one out of two fates (A or B, bistable regime) or, additionally, remaining in their progenitor state (C, tristable regime). Such a differentiation into populations with three distinct cell fates is observed during blastocyst formation in mammals, where inner cell mass cells can remain in that state or differentiate into epiblast cells or primitive endoderm. Systematic analysis revealed that the existence of a stable non-differentiated state significantly impacts the GRN's robustness against parametric variations and stochastic noise. This state reduces the sensitivity of cell populations to parameters controlling key gene expression asymmetry and prevents cells from making transitions after acquiring a new identity. Stochastic noise enhances robustness by decreasing sensitivity to initial expression levels and helping the system escape from the non-differentiated state to differentiated cell fates, making the differentiation more efficient.

Coupling of cell shape, matrix and tissue dynamics ensures embryonic patterning robustness

Tissue patterning coordinates morphogenesis, cell dynamics and fate specification. Understanding how precision in patterning is robustly achieved despite inherent developmental variability during mammalian embryogenesis remains a challenge. Here, based on cell dynamics quantification and simulation, we show how salt-and-pepper epiblast and primitive endoderm (PrE) cells pattern the inner cell mass of mouse blastocysts. Coupling cell fate and dynamics, PrE cells form apical polarity-dependent actin protrusions required for RAC1-dependent migration towards the surface of the fluid cavity, where PrE cells are trapped due to decreased tension. Concomitantly, PrE cells deposit an extracellular matrix gradient, presumably breaking the tissue-level symmetry and collectively guiding their own migration. Tissue size perturbations of mouse embryos and their comparison with monkey and human blastocysts further demonstrate that the fixed proportion of PrE/epiblast cells is optimal with respect to embryo size and tissue geometry and, despite variability, ensures patterning robustness during early mammalian development.

Divergent destinies: insights into the molecular mechanisms underlying EPI and PE fate determination

Mammalian pre-implantation development is entirely devoted to the specification of extra-embryonic lineages, which are fundamental for embryo morphogenesis and support. The second fate decision is taken just before implantation, as defined by the epiblast (EPI) and the primitive endoderm (PE) specification. Later, EPI forms the embryo proper and PE contributes to the formation of the yolk sac. The formation of EPI and PE as molecularly and morphologically distinct lineages is the final step of a multistage process, which begins when bipotent progenitor cells diverge into separate fates. Despite advances in uncovering the molecular mechanisms underlying the differential transcriptional patterns that dictate how apparently identical cells make fate decisions and how lineage integrity is maintained, a detailed overview of these mechanisms is still lacking. In this review, we dissect the EPI and PE formation process into four stages (initiation, specification, segregation, and maintenance) and we provide a comprehensive understanding of the molecular mechanisms involved in lineage establishment in the mouse. In addition, we discuss the conservation of key processes in humans, based on the most recent findings.

PI3K/AKT signaling controls ICM maturation and proper epiblast and primitive endoderm specification in mice

The inner cell mass (ICM) of early mouse embryos is specified into epiblast (Epi) and primitive endoderm (PrE) lineages during blastocyst formation. The antagonistic transcription factors (TFs) NANOG and GATA-binding protein 6 (GATA6) in combination with fibroblast growth factor (FGF)/extracellular-signal-regulated kinase (ERK) signaling are central actors in ICM fate choice. However, what initiates the specification of ICM progenitors into Epi or PrE and whether other factors are involved in this process has not been fully understood yet. Here, we show that phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT) is constitutively active during preimplantation development. Using pharmacological inhibition, we demonstrate that PI3K/AKT enables the formation of a functional ICM capable of giving rise to both the Epi and the PrE: it maintains the expression of the TF NANOG, which specifies the Epi, and confers responsiveness to FGF4, which is essential for PrE specification. Our work thus identifies PI3K/AKT signaling as an upstream regulator controlling the molecular events required for both Epi and PrE specification.

Fate specification triggers a positive feedback loop of TEAD-YAP and NANOG to promote epiblast formation in preimplantation embryos

In preimplantation embryos, epiblast (EPI) fate specification from the inner cell mass is controlled by the segregation of NANOG and GATA6 expression. TEAD-YAP interaction is activated during EPI formation and is required for pluripotency factor expression. These events occur asynchronously with similar timing during EPI formation, and their relationship remains elusive. Here, we examined the relationship between NANOG-GATA6 and TEAD-YAP. The nuclear accumulation of YAP takes place only in EPI-specified cells, and a positive feedback loop operates between NANOG and TEAD-YAP. The effects of TEAD-YAP on SOX2 upregulation in EPI-specified cells are likely indirect. EPI fate specification also alters the response of Nanog, Sox2 and Cdx2 to TEAD-YAP. These results suggest that EPI-fate specification alters the transcriptional network from the morula-like to the EPI-specified state and activates TEAD-YAP to trigger a positive feedback loop with NANOG, which stabilizes the EPI fate. The coordinated occurrence of these processes in individual cells likely supports proper EPI formation under the condition of asynchronous EPI-fate specification.

Common variants increase risk for congenital diaphragmatic hernia within the context of de novo variants

Congenital diaphragmatic hernia (CDH) is a severe congenital anomaly often accompanied by other structural anomalies and/or neurobehavioral manifestations. Rare de novo protein-coding variants and copy-number variations contribute to CDH in the population. However, most individuals with CDH remain genetically undiagnosed. Here, we perform integrated de novo and common-variant analyses using 1,469 CDH individuals, including 1,064 child-parent trios and 6,133 ancestry-matched, unaffected controls for the genome-wide association study. We identify candidate CDH variants in 15 genes, including eight novel genes, through deleterious de novo variants. We further identify two genomic loci contributing to CDH risk through common variants with similar effect sizes among Europeans and Latinx. Both loci are in putative transcriptional regulatory regions of developmental patterning genes. Estimated heritability in common variants is ∼19%. Strikingly, there is no significant difference in estimated polygenic risk scores between isolated and complex CDH or between individuals harboring deleterious de novo variants and individuals without these variants. The data support a polygenic model as part of the CDH genetic architecture.

Common modes of ERK induction resolve into context-specific signalling via emergent networks and cell-type-specific transcriptional repression

Fibroblast Growth Factor signalling via ERK exerts diverse roles in development and disease. In mammalian preimplantation embryos and naïve pluripotent stem cells ERK promotes differentiation, whereas in primed pluripotent states closer to somatic differentiation ERK sustains self-renewal. How can the same pathway produce different outcomes in two related cell types? To explore context-dependent ERK signalling we generated cell and mouse lines that allow for tissue- and time-specific ERK activation. Using these tools, we find that specificity in ERK response is mostly mediated by repression of transcriptional targets that occur in tandem with reductions in chromatin accessibility at regulatory regions. Furthermore, immediate early ERK responses are largely shared by different cell types but produce cell-specific programmes as these responses interface with emergent networks in the responding cells. Induction in naïve pluripotency is accompanied by chromatin changes, whereas in later stages it is not, suggesting that chromatin context does not shape signalling response. Altogether, our data suggest that cell-type-specific responses to ERK signalling exploit the same immediate early response, but then sculpt it to specific lineages via repression of distinct cellular programmes.

AI-powered simulation-based inference of a genuinely spatial-stochastic gene regulation model of early mouse embryogenesis

Understanding how multicellular organisms reliably orchestrate cell-fate decisions is a central challenge in developmental biology, particularly in early mammalian development, where tissue-level differentiation arises from seemingly cell-autonomous mechanisms. In this study, we present a multi-scale, spatial-stochastic simulation framework for mouse embryogenesis, focusing on inner cell mass (ICM) differentiation into epiblast (EPI) and primitive endoderm (PRE) at the blastocyst stage. Our framework models key regulatory and tissue-scale interactions in a biophysically realistic fashion, capturing the inherent stochasticity of intracellular gene expression and intercellular signaling, while efficiently simulating these processes by advancing event-driven simulation techniques. Leveraging the power of Simulation-Based Inference (SBI) through the AI-driven Sequential Neural Posterior Estimation (SNPE) algorithm, we conduct a large-scale Bayesian inferential analysis to identify parameter sets that faithfully reproduce experimentally observed features of ICM specification. Our results reveal mechanistic insights into how the combined action of autocrine and paracrine FGF4 signaling coordinates stochastic gene expression at the cellular scale to achieve robust and reproducible ICM patterning at the tissue scale. We further demonstrate that the ICM exhibits a specific time window of sensitivity to exogenous FGF4, enabling lineage proportions to be adjusted based on timing and dosage, thereby extending current experimental findings and providing quantitative predictions for both mutant and wild-type ICM systems. Notably, FGF4 signaling not only ensures correct EPI-PRE lineage proportions but also enhances ICM resilience to perturbations, reducing fate-proportioning errors by 10-20% compared to a purely cell-autonomous system. Additionally, we uncover a surprising role for variability in intracellular initial conditions, showing that high gene-expression heterogeneity can improve both the accuracy and precision of cell-fate proportioning, which remains robust when fewer than 25% of the ICM population experiences perturbed initial conditions. Our work offers a comprehensive, spatial-stochastic description of the biochemical processes driving ICM differentiation and identifies the necessary conditions for its robust unfolding. It also provides a framework for future exploration of similar spatial-stochastic systems in developmental biology.

EnhancerNet: a predictive model of cell identity dynamics through enhancer selection

Understanding how cell identity is encoded by the genome and acquired during differentiation is a central challenge in cell biology. I have developed a theoretical framework called EnhancerNet, which models the regulation of cell identity through the lens of transcription factor-enhancer interactions. I demonstrate that autoregulation in these interactions imposes a constraint on the model, resulting in simplified dynamics that can be parameterized from observed cell identities. Despite its simplicity, EnhancerNet recapitulates a broad range of experimental observations on cell identity dynamics, including enhancer selection, cell fate induction, hierarchical differentiation through multipotent progenitor states and direct reprogramming by transcription factor overexpression. The model makes specific quantitative predictions, reproducing known reprogramming recipes and the complex haematopoietic differentiation hierarchy without fitting unobserved parameters. EnhancerNet provides insights into how new cell types could evolve and highlights the functional importance of distal regulatory elements with dynamic chromatin in multicellular evolution.

Naive pluripotent stem cell-based models capture FGF-dependent human hypoblast lineage specification

The hypoblast is an essential extraembryonic tissue set aside within the inner cell mass in the blastocyst. Research with human embryos is challenging. Thus, stem cell models that reproduce hypoblast differentiation provide valuable alternatives. We show here that human naive pluripotent stem cell (PSC) to hypoblast differentiation proceeds via reversion to a transitional ICM-like state from which the hypoblast emerges in concordance with the trajectory in human blastocysts. We identified a window when fibroblast growth factor (FGF) signaling is critical for hypoblast specification. Revisiting FGF signaling in human embryos revealed that inhibition in the early blastocyst suppresses hypoblast formation. In vitro, the induction of hypoblast is synergistically enhanced by limiting trophectoderm and epiblast fates. This finding revises previous reports and establishes a conservation in lineage specification between mice and humans. Overall, this study demonstrates the utility of human naive PSC-based models in elucidating the mechanistic features of early human embryogenesis.

Living in a noisy world-origins of gene expression noise and its impact on cellular decision-making

The expression level of a gene can vary between genetically identical cells under the same environmental condition-a phenomenon referred to as gene expression noise. Several studies have now elucidated a central role of transcription factors in the generation of expression noise. Transcription factors, as the key components of gene regulatory networks, drive many important cellular decisions in response to cellular and environmental signals. Therefore, a very relevant question is how expression noise impacts gene regulation and influences cellular decision-making. In this Review, we summarize the current understanding of the molecular origins of expression noise, highlighting the role of transcription factors in this process, and discuss the ways in which noise can influence cellular decision-making. As advances in single-cell technologies open new avenues for studying expression noise as well as gene regulatory circuits, a better understanding of the influence of noise on cellular decisions will have important implications for many biological processes.

Exploiting transcription factors to target EMT and cancer stem cells for tumor modulation and therapy

Transcription factors (TFs) are essential in controlling gene regulatory networks that determine cellular fate during embryogenesis and tumor development. TFs are the major players in promoting cancer stemness by regulating the function of cancer stem cells (CSCs). Understanding how TFs interact with their downstream targets for determining cell fate during embryogenesis and tumor development is a critical area of research. CSCs are increasingly recognized for their significance in tumorigenesis and patient prognosis, as they play a significant role in cancer initiation, progression, metastasis, and treatment resistance. However, traditional therapies have limited effectiveness in eliminating this subset of cells, allowing CSCs to persist and potentially form secondary tumors. Recent studies have revealed that cancer cells and tumors with CSC-like features also exhibit genes related to the epithelial-to-mesenchymal transition (EMT). EMT-associated transcription factors (EMT-TFs) like TWIST and Snail/Slug can upregulate EMT-related genes and reprogram cancer cells into a stem-like phenotype. Importantly, the regulation of EMT-TFs, particularly through post-translational modifications (PTMs), plays a significant role in cancer metastasis and the acquisition of stem cell-like features. PTMs, including phosphorylation, ubiquitination, and SUMOylation, can alter the stability, localization, and activity of EMT-TFs, thereby modulating their ability to drive EMT and stemness properties in cancer cells. Although targeting EMT-TFs holds potential in tackling CSCs, current pharmacological approaches to do so directly are unavailable. Therefore, this review aims to explore the role of EMT- and CSC-TFs, their connection and impact in cellular development and cancer, emphasizing the potential of TF networks as targets for therapeutic intervention.

A dynamical perspective: moving towards mechanism in single-cell transcriptomics

As the field of single-cell transcriptomics matures, research is shifting focus from phenomenological descriptions of cellular phenotypes to a mechanistic understanding of the gene regulation underneath. This perspective considers the value of capturing dynamical information at single-cell resolution for gaining mechanistic insight; reviews the available technologies for recording and inferring temporal information in single cells; and explores whether better dynamical resolution is sufficient to adequately capture the causal relationships driving complex biological systems. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Using extended growth of cattle embryos in culture to gain insights into bovine developmental events on embryonic days 8 to 10

We have previously described an extended embryo culture system, based on uterine composition, growth factors and the cell culture additive B27, for growing cattle embryos in vitro beyond embryonic day 7. Here, extended in vitro embryos are compared to embryos developed in the uterus and are used to establish a developmental staging framework useful for understanding developmental events occurring until Day 10. Immunofluorescence or mRNA expression of the ICM/epiblast markers OCT4, SOX2 and NANOG, hypoblast markers GATA6, SOX17 and GATA4 and trophoblast genes CDX2, GATA3, ASCL2 and IFNT revealed the presence of four stages during this period that can be molecularly distinguished. These are expanded blastocyst, hatched blastocyst, hypoblast layering and early hypoblast migration. Interestingly NANOG and SOX17 show reciprocal expression at the expanded blastocyst stage, well before SOX2 and GATA6 expression refines to a similar so-called "salt and pepper" mutually exclusive expression in the ICM at the hatched blastocyst stage. GATA4 expression is only seen from stages when the hypoblast starts migrating around the blastocyst cavity. Intriguingly, trophoblast still expresses GATA6 and OCT4 in all cells during the expanded blastocyst phase, while SOX2 and SOX17 are seen in only some trophoblast cells. By the hypoblast-epiblast layering stage no trophoblast expression remains except for that of OCT4 protein, which starts waning in trophoblast once the hypoblast begins migrating. Lastly, it is shown that cultured embryos exhibit increased expression of the stress marker TP53 in the epiblast and hypoblast at late stages in comparison to embryos produced in the uterine environment.

NANOG expression in parthenogenetic porcine blastocysts is required for intact lineage specification and pluripotency

OBJECTIVE

Nanog homeobox (NANOG) is a core transcription factor that contributes to pluripotency along with octamer binding transcription factor-4 (OCT4) and sex determining region-Y box-2 (SOX2). It is an epiblast lineage marker in mammalian pre-implantation embryos and exhibits a species-specific expression pattern. Therefore, it is important to understand the lineage of NANOG, the trophectoderm, and the primitive endoderm in the pig embryo.

METHODS

A loss- and gain-of-function analysis was done to determine the role of NANOG in lineage specification in parthenogenetic porcine blastocysts. We analyzed the relationship between NANOG and pluripotent core transcription factors and other lineage makers.

RESULTS

In NANOG-null late blastocysts, OCT4-, SOX2-, and SOX17-positive cells were decreased, whereas GATA binding protein 6 (GATA6)-positive cells were increased. Quantitative real-time polymerase chain reaction revealed that the expression of SOX2 was decreased in NANOG-null blastocysts, whereas that of primitive endoderm makers, except SOX17, was increased. In NANOG-overexpressing blastocysts, caudal type homeobox 2 (CDX2-), SOX17-, and GATA6-positive cells were decreased. The results indicated that the expression of primitive endoderm markers and trophectoderm-related genes was decreased.

CONCLUSION

Taken together, the results demonstrate that NANOG is involved in the epiblast and primitive endoderm differentiation and is essential for maintaining pluripotency within the epiblast.

In vitro generation of mouse morula-like cells

Generating cells with the molecular and functional properties of embryo cells and with full developmental potential is an aim with fundamental biological significance. Here we report the in vitro generation of mouse transient morula-like cells (MLCs) via the manipulation of signaling pathways. MLCs are molecularly distinct from embryonic stem cells (ESCs) and cluster instead with embryo 8- to 16-cell stage cells. A single MLC can generate a blastoid, and the efficiency increases to 80% when 8-10 MLCs are used. MLCs make embryoids directly, efficiently, and within 4 days. Transcriptomic analysis shows that day 4-5 MLC-derived embryoids contain the cell types found in natural embryos at early gastrulation. Furthermore, MLCs introduced into morulae segregate into epiblast (EPI), primitive endoderm (PrE), and trophectoderm (TE) fates in blastocyst chimeras and have a molecular signature indistinguishable from that of host embryo cells. These findings represent the generation of cells that are molecularly and functionally similar to the precursors of the first three cell lineages of the embryo.

The salt-and-pepper pattern in mouse blastocysts is compatible with signaling beyond the nearest neighbors

Embryos develop in a concerted sequence of spatiotemporal arrangements of cells. In the preimplantation mouse embryo, the distribution of the cells in the inner cell mass evolves from a salt-and-pepper pattern to spatial segregation of two distinct cell types. The exact properties of the salt-and-pepper pattern have not been analyzed so far. We investigate the spatiotemporal distribution of NANOG- and GATA6-expressing cells in the ICM of the mouse blastocysts with quantitative three-dimensional single-cell-based neighborhood analyses. A combination of spatial statistics and agent-based modeling reveals that the cell fate distribution follows a local clustering pattern. Using ordinary differential equations modeling, we show that this pattern can be established by a distance-based signaling mechanism enabling cells to integrate information from the whole inner cell mass into their cell fate decision. Our work highlights the importance of longer-range signaling to ensure coordinated decisions in groups of cells to successfully build embryos.

Local cellular interactions during the self-organization of stem cells

Stem cell models for early mammalian development offer new experimental opportunities to access spatio-temporal details of the cell-cell interactions that govern cell differentiation and tissue patterning. This review summarizes recent studies that have used stem cell models to investigate the spatial range of developmental cell-cell communication systems. A key message from these works is that important biochemical signals for cell differentiation in these systems, such as Nodal and fibroblast growth factors (FGFs), often act over short distances of only a few cell diameters. The formation of long-range patterns at the tissue scale associated with these signals then results from signal relays and cell rearrangements. The modular view of differentiation and patterning emerging from research on stem cell models can offer a fresh perspective on the corresponding processes in the embryo.

Constructing maps between distinct cell fates and parametric conditions by systematic perturbations

MOTIVATION

Cell fate transitions are common in many developmental processes. Therefore, identifying the mechanisms behind them is crucial. Traditionally, due to complexity of networks and existence of plenty of kinetic parameters, dynamical analysis of biomolecular networks can only be performed by simultaneously perturbing a small number of parameters. Although many efforts have focused on how cell states change under specific perturbations, conversely, how to infer parametric conditions underlying distinct cell fates by systematic perturbations is less clear and needs to be further investigated.

RESULTS

In this article, we present a general computational method by integrating systematic perturbations, unsupervised clustering, principal component analysis, and fitting analysis. The method can be used to to construct maps between distinct cell fates and parametric conditions by systematic perturbations. In particular, there are no needs of accurate parameter measurements and occurrence of bifurcations to establish the maps. To validate feasibility and inference performance of the method, we use toggle switch, inner cell mass, and epithelial mesenchymal transition as model systems to show how the maps are constructed and how system parameters encode essential information on cell fates. The maps tell us how systematic perturbations drive cell fate decisions and transitions, and allow us to purposefully predict, manipulate, and even control cell states. The approach is especially helpful in understanding crucial roles of certain parameter combinations during fate transitions. We hope that the approach can provide us valuable information on parametric or perturbation conditions so some specific targets, e.g. directional differentiation, can be realized.

AVAILABILITY AND IMPLEMENTATION

No public data are used. The data we used are generated by randomly chosen values of model parameters in certain ranges, and the corresponding parameters are already attached in supplementary materials.

Transcriptional network governing extraembryonic endoderm cell fate choice

The mechanism by which transcription factor (TF) network instructs cell-type-specific transcriptional programs to drive primitive endoderm (PrE) progenitors to commit to parietal endoderm (PE) versus visceral endoderm (VE) cell fates remains poorly understood. To address the question, we analyzed the single-cell transcriptional signatures defining PrE, PE, and VE cell states during the onset of the PE-VE lineage bifurcation. By coupling with the epigenomic comparison of active enhancers unique to PE and VE cells, we identified GATA6, SOX17, and FOXA2 as central regulators for the lineage divergence. Transcriptomic analysis of cXEN cells, an in vitro model for PE cells, after the acute depletion of GATA6 or SOX17 demonstrated that these factors induce Mycn, imparting the self-renewal properties of PE cells. Concurrently, they suppress the VE gene program, including key genes like Hnf4a and Ttr, among others. We proceeded with RNA-seq analysis on cXEN cells with FOXA2 knockout, in conjunction with GATA6 or SOX17 depletion. We found FOXA2 acts as a potent suppressor of Mycn while simultaneously activating the VE gene program. The antagonistic gene regulatory activities of GATA6/SOX17 and FOXA2 in promoting alternative cell fates, and their physical co-bindings at the enhancers provide molecular insights to the plasticity of the PrE lineage. Finally, we show that the external cue, BMP signaling, promotes the VE cell fate by activation of VE TFs and repression of PE TFs including GATA6 and SOX17. These data reveal a putative core gene regulatory module that underpins PE and VE cell fate choice.

Recognition and reconstruction of cell differentiation patterns with deep learning

Cell lineage decisions occur in three-dimensional spatial patterns that are difficult to identify by eye. There is an ongoing effort to replicate such patterns using mathematical modeling. One approach uses long ranging cell-cell communication to replicate common spatial arrangements like checkerboard and engulfing patterns. In this model, the cell-cell communication has been implemented as a signal that disperses throughout the tissue. On the other hand, machine learning models have been developed for pattern recognition and pattern reconstruction tasks. We combined synthetic data generated by the mathematical model with spatial summary statistics and deep learning algorithms to recognize and reconstruct cell fate patterns in organoids of mouse embryonic stem cells. Application of Moran's index and pair correlation functions for in vitro and synthetic data from the model showed local clustering and radial segregation. To assess the patterns as a whole, a graph neural network was developed and trained on synthetic data from the model. Application to in vitro data predicted a low signal dispersion value. To test this result, we implemented a multilayer perceptron for the prediction of a given cell fate based on the fates of the neighboring cells. The results show a 70% accuracy of cell fate imputation based on the nine nearest neighbors of a cell. Overall, our approach combines deep learning with mathematical modeling to link cell fate patterns with potential underlying mechanisms.

A bipartite function of ESRRB can integrate signaling over time to balance self-renewal and differentiation

Cooperative DNA binding of transcription factors (TFs) integrates the cellular context to support cell specification during development. Naive mouse embryonic stem cells are derived from early development and can sustain their pluripotent identity indefinitely. Here, we ask whether TFs associated with pluripotency evolved to directly support this state or if the state emerges from their combinatorial action. NANOG and ESRRB are key pluripotency factors that co-bind DNA. We find that when both factors are expressed, ESRRB supports pluripotency. However, when NANOG is absent, ESRRB supports a bistable culture of cells with an embryo-like primitive endoderm identity ancillary to pluripotency. The stoichiometry between NANOG and ESRRB allows quantitative titration of this differentiation, and in silico modeling of bipartite ESRRB activity suggests it safeguards plasticity in differentiation. Thus, the concerted activity of cooperative TFs can transform their effect to sustain intermediate cell identities and allow ex vivo expansion of immortal stem cells. A record of this paper's transparent peer review process is included in the supplemental information.

insideOutside: an accessible algorithm for classifying interior and exterior points, with applications in embryology

A crucial aspect of embryology is relating the position of individual cells to the broader geometry of the embryo. A classic example of this is the first cell-fate decision of the mouse embryo, where interior cells become inner cell mass and exterior cells become trophectoderm. Fluorescent labelling, imaging, and quantification of tissue-specific proteins have advanced our understanding of this dynamic process. However, instances arise where these markers are either not available, or not reliable, and we are left only with the cells' spatial locations. Therefore, a simple, robust method for classifying interior and exterior cells of an embryo using spatial information is required. Here, we describe a simple mathematical framework and an unsupervised machine learning approach, termed insideOutside, for classifying interior and exterior points of a three-dimensional point-cloud, a common output from imaged cells within the early mouse embryo. We benchmark our method against other published methods to demonstrate that it yields greater accuracy in classification of nuclei from the pre-implantation mouse embryos and greater accuracy when challenged with local surface concavities. We have made MATLAB and Python implementations of the method freely available. This method should prove useful for embryology, with broader applications to similar data arising in the life sciences.

Adjusting the range of cell-cell communication enables fine-tuning of cell fate patterns from checkerboard to engulfing

During development, spatio-temporal patterns ranging from checkerboard to engulfing occur with precise proportions of the respective cell fates. Key developmental regulators are intracellular transcriptional interactions and intercellular signaling. We present an analytically tractable mathematical model based on signaling that reliably generates different cell type patterns with specified proportions. Employing statistical mechanics, We derived a cell fate decision model for two cell types. A detailed steady state analysis on the resulting dynamical system yielded necessary conditions to generate spatially heterogeneous patterns. This allows the cell type proportions to be controlled by a single model parameter. Cell-cell communication is realized by local and global signaling mechanisms. These result in different cell type patterns. A nearest neighbor signal yields checkerboard patterns. Increasing the signal dispersion, cell fate clusters and an engulfing pattern can be generated. Altogether, the presented model allows us to reliably generate heterogeneous cell type patterns of different kinds as well as desired proportions.

Young LINE-1 transposon 5' UTRs marked by elongation factor ELL3 function as enhancers to regulate naïve pluripotency in embryonic stem cells

LINE-1s are the major clade of retrotransposons with autonomous retrotransposition activity. Despite the potential genotoxicity, LINE-1s are highly activated in early embryos. Here we show that a subset of young LINE-1s, L1Md_Ts, are marked by the RNA polymerase II elongation factor ELL3, and function as enhancers in mouse embryonic stem cells. ELL3 depletion dislodges the DNA hydroxymethylase TET1 and the co-repressor SIN3A from L1Md_Ts, but increases the enrichment of the Bromodomain protein BRD4, leading to loss of 5hmC, gain of H3K27ac, and upregulation of the L1Md_T nearby genes. Specifically, ELL3 occupies and represses the L1Md_T-based enhancer located within Akt3, which encodes a key regulator of AKT pathway. ELL3 is required for proper ERK activation and efficient shutdown of naïve pluripotency through inhibiting Akt3 during naïve-primed transition. Our study reveals that the enhancer function of a subset of young LINE-1s controlled by ELL3 in transcription regulation and mouse early embryo development.

Single-cell analysis of bidirectional reprogramming between early embryonic states reveals mechanisms of differential lineage plasticities

Two distinct fates, pluripotent epiblast (EPI) and primitive (extra-embryonic) endoderm (PrE), arise from common progenitor cells, the inner cell mass (ICM), in mammalian embryos. To study how these sister identities are forged, we leveraged embryonic (ES) and eXtraembryonic ENdoderm (XEN) stem cells - in vitro counterparts of the EPI and PrE. Bidirectional reprogramming between ES and XEN coupled with single-cell RNA and ATAC-seq analyses uncovered distinct rates, efficiencies and trajectories of state conversions, identifying drivers and roadblocks of reciprocal conversions. While GATA4-mediated ES-to-iXEN conversion was rapid and nearly deterministic, OCT4, KLF4 and SOX2-induced XEN-to-iPS reprogramming progressed with diminished efficiency and kinetics. The dominant PrE transcriptional program, safeguarded by Gata4, and globally elevated chromatin accessibility of EPI underscored the differential plasticities of the two states. Mapping in vitro trajectories to embryos revealed reprogramming in either direction tracked along, and toggled between, EPI and PrE in vivo states without transitioning through the ICM.

The Wnt/TCF7L1 transcriptional repressor axis drives primitive endoderm formation by antagonizing naive and formative pluripotency

Early during preimplantation development and in heterogeneous mouse embryonic stem cells (mESC) culture, pluripotent cells are specified towards either the primed epiblast or the primitive endoderm (PE) lineage. Canonical Wnt signaling is crucial for safeguarding naive pluripotency and embryo implantation, yet the role and relevance of canonical Wnt inhibition during early mammalian development remains unknown. Here, we demonstrate that transcriptional repression exerted by Wnt/TCF7L1 promotes PE differentiation of mESCs and in preimplantation inner cell mass. Time-series RNA sequencing and promoter occupancy data reveal that TCF7L1 binds and represses genes encoding essential naive pluripotency factors and indispensable regulators of the formative pluripotency program, including Otx2 and Lef1. Consequently, TCF7L1 promotes pluripotency exit and suppresses epiblast lineage formation, thereby driving cells into PE specification. Conversely, TCF7L1 is required for PE specification as deletion of Tcf7l1 abrogates PE differentiation without restraining epiblast priming. Taken together, our study underscores the importance of transcriptional Wnt inhibition in regulating lineage specification in ESCs and preimplantation embryo development as well as identifies TCF7L1 as key regulator of this process.

Inhibition of FGF receptor impairs primitive endoderm differentiation in bovine embryos

The first cellular differentiation event in the pre-implantation embryo results in the trophectoderm (TE) and the inner cell mass (ICM). A second event occurs in the latter, resulting in the epiblast and the primitive endoderm (PE). This second differentiation is still not fully characterized in bovine development, although it is likely to involve FGF signalling. Thus, in this study, we tested the hypothesis that stimulation or inhibition of the FGF pathway during bovine embryo in vitro culture would only interfere with PE differentiation if maintained until later blastocyst stages. At first, we characterized the expression of PE marker SOX17 at different blastocyst stages. Then, we treated in vitro produced embryos during different windows of time: days 5.0-7.0 (D5-D7), D7-D9, and D5-D9 with 1 μg/ml FGF4 and 1 μg/ml heparin or 1 mM FGFR inhibitor, AZD4547. We observed that the SOX17-positive cell number only increases in late-stage blastocysts compared to early stages. Treatment of embryos with FGF4 did not change the number of SOX17-positive cells, while inhibition of FGFR signalling reduced SOX17-positive cells from D5-D7 and completely ablated SOX17 expression when kept until D9. In conclusion, FGFR inhibition repressed PE differentiation in bovine embryos at all time points, although stimulation with FGF4 did not interfere with PE cell numbers.

Biophysical models of early mammalian embryogenesis

Embryo development is a critical and fascinating stage in the life cycle of many organisms. Despite decades of research, the earliest stages of mammalian embryogenesis are still poorly understood, caused by a scarcity of high-resolution spatial and temporal data, the use of only a few model organisms, and a paucity of truly multidisciplinary approaches that combine biological research with biophysical modeling and computational simulation. Here, we explain the theoretical frameworks and biophysical processes that are best suited to modeling the early mammalian embryo, review a comprehensive list of previous models, and discuss the most promising avenues for future work.

Early Cell Specification in Mammalian Fertilized and Somatic Cell Nuclear Transfer Embryos

Early cell specification in mammalian preimplantation embryos is an intricate cellular process that leads to coordinated spatial and temporal expression of specific genes. Proper segregation into the first two cell lineages, the inner cell mass (ICM) and the trophectoderm (TE), is imperative for developing the embryo proper and the placenta, respectively. Somatic cell nuclear transfer (SCNT) allows the formation of a blastocyst containing both ICM and TE from a differentiated cell nucleus, which means that this differentiated genome must be reprogrammed to a totipotent state. Although blastocysts can be generated efficiently through SCNT, the full-term development of SCNT embryos is impaired mostly due to placental defects. In this review, we examine the early cell fate decisions in fertilized embryos and compare them to observations in SCNT-derived embryos, in order to understand if these processes are affected by SCNT and could be responsible for the low success of reproductive cloning.

Journey of the mouse primitive endoderm: from specification to maturation

The blastocyst is a conserved stage and distinct milestone in the development of the mammalian embryo. Blastocyst stage embryos comprise three cell lineages which arise through two sequential binary cell fate specification steps. In the first, extra-embryonic trophectoderm (TE) cells segregate from inner cell mass (ICM) cells. Subsequently, ICM cells acquire a pluripotent epiblast (Epi) or extra-embryonic primitive endoderm (PrE, also referred to as hypoblast) identity. In the mouse, nascent Epi and PrE cells emerge in a salt-and-pepper distribution in the early blastocyst and are subsequently sorted into adjacent tissue layers by the late blastocyst stage. Epi cells cluster at the interior of the ICM, while PrE cells are positioned on its surface interfacing the blastocyst cavity, where they display apicobasal polarity. As the embryo implants into the maternal uterus, cells at the periphery of the PrE epithelium, at the intersection with the TE, break away and migrate along the TE as they mature into parietal endoderm (ParE). PrE cells remaining in association with the Epi mature into visceral endoderm. In this review, we discuss our current understanding of the PrE from its specification to its maturation. This article is part of the theme issue 'Extraembryonic tissues: exploring concepts, definitions and functions across the animal kingdom'.

The physiological and pathological mechanisms of early embryonic development

Early embryonic development is a complex process. The zygote undergoes several rounds of division to form a blastocyst, and during this process, the zygote undergoes the maternal-to-zygotic transition to gain control of embryonic development and makes two cell fate decisions to differentiate into an embryonic and two extra-embryonic lineages. With the use of new molecular biotechnologies and animal models, we can now further study the molecular mechanisms of early embryonic development and the pathological causes of early embryonic arrest. Here, we first summarize the known molecular regulatory mechanisms of early embryonic development in mice. Then we discuss the pathological factors leading to the early embryonic arrest. We hope that this review will give researchers a relatively complete view of the physiology and pathology of early embryonic development.

Paracrine interactions through FGFR1 and FGFR2 receptors regulate the development of preimplantation mouse chimaeric embryo

The preimplantation mammalian embryo has the potential to self-organize, allowing the formation of a correctly patterned embryo despite experimental perturbation. To better understand the mechanisms controlling the developmental plasticity of the early mouse embryo, we used chimaeras composed of an embryonic day (E)3.5 or E4.5 inner cell mass (ICM) and cleaving 8-cell embryo. We revealed that the restricted potential of the ICM can be compensated for by uncommitted 8-cell embryo-derived blastomeres, thus leading to the formation of a normal chimaeric blastocyst that can undergo full development. However, whether such chimaeras maintain developmental competence depends on the presence or specific orientation of the polarized primitive endoderm layer in the ICM component. We also demonstrated that downregulated FGFR1 and FGFR2 expression in 8-cell embryos disturbs intercellular interactions between both components and results in an inverse proportion of primitive endoderm and epiblast within the resulting ICM and abnormal embryo development. This finding suggests that FGF signalling is a key part of the regulatory mechanism that assigns cells to a given lineage and ensures the proper composition of the blastocyst, which is a prerequisite for its successful implantation in the uterus and for further development.

A pendulum of induction between the epiblast and extra-embryonic endoderm supports post-implantation progression

Embryogenesis is supported by dynamic loops of cellular interactions. Here, we create a partial mouse embryo model to elucidate the principles of epiblast (Epi) and extra-embryonic endoderm co-development (XEn). We trigger naive mouse embryonic stem cells to form a blastocyst-stage niche of Epi-like cells and XEn-like cells (3D, hydrogel free and serum free). Once established, these two lineages autonomously progress in minimal medium to form an inner pro-amniotic-like cavity surrounded by polarized Epi-like cells covered with visceral endoderm (VE)-like cells. The progression occurs through reciprocal inductions by which the Epi supports the primitive endoderm (PrE) to produce a basal lamina that subsequently regulates Epi polarization and/or cavitation, which, in return, channels the transcriptomic progression to VE. This VE then contributes to Epi bifurcation into anterior- and posterior-like states. Similarly, boosting the formation of PrE-like cells within blastoids supports developmental progression. We argue that self-organization can arise from lineage bifurcation followed by a pendulum of induction that propagates over time.

A NANOG‐pERK reciprocal regulatory circuit regulates Nanog autoregulation and ERK signaling dynamics

The self‐renewal and differentiation potential of embryonic stem cells (ESCs) is maintained by the regulated expression of core pluripotency factors. Expression levels of the core pluripotency factor Nanog are tightly regulated by a negative feedback autorepression loop. However, it remains unclear how ESCs perceive NANOG levels and execute autorepression. Here, we show that a dose‐dependent induction of Fgfbp1 and Fgfr2 by NANOG activates autocrine‐mediated ERK signaling in Nanog‐high cells to trigger autorepression. pERK recruits NONO to the Nanog locus to repress transcription by preventing POL2 loading. This Nanog autorepression process establishes a self‐perpetuating reciprocal NANOG‐pERK regulatory circuit. We further demonstrate that this reciprocal regulatory circuit induces pERK heterogeneity and ERK signaling dynamics in pluripotent stem cells. Collectively our data suggest that NANOG induces Fgfr2 and Fgfbp1 to activate ERK signaling in Nanog‐high cells to establish a NANOG‐pERK reciprocal regulatory circuit. This circuit regulates ERK signaling dynamics and Nanog autoregulation in pluripotent cells.

Initial source of heterogeneity in a model for cell fate decision in the early mammalian embryo

During development, cells from a population of common progenitors evolve towards different fates characterized by distinct levels of specific transcription factors, a process known as cell differentiation. This evolution is governed by gene regulatory networks modulated by intercellular signalling. In order to evolve towards distinct fates, cells forming the population of common progenitors must display some heterogeneity. We applied a modelling approach to obtain insights into the possible sources of cell-to-cell variability initiating the specification of cells of the inner cell mass into epiblast or primitive endoderm cells in early mammalian embryo. At the single-cell level, these cell fates correspond to three possible steady states of the model. A combination of numerical simulations and bifurcation analyses predicts that the behaviour of the model is preserved with respect to the source of variability and that cell-cell coupling induces the emergence of multiple steady states associated with various cell fate configurations, and to a distribution of the levels of expression of key transcription factors. Statistical analysis of these time-dependent distributions reveals differences in the evolutions of the variance-to-mean ratios of key variables of the system, depending on the simulated source of variability, and, by comparison with experimental data, points to the rate of synthesis of the key transcription factor NANOG as a likely initial source of heterogeneity.

TEAD4 regulates trophectoderm differentiation upstream of CDX2 in a GATA3-independent manner in the human preimplantation embryo

STUDY QUESTION

What is the role of transcriptional-enhanced associate (TEA) domain family member 4 (TEAD4) in trophectoderm (TE) differentiation during human embryo preimplantation development in comparison to mouse?

SUMMARY ANSWER

TEAD4 regulates TE lineage differentiation in the human preimplantation embryo acting upstream of caudal-type homeobox protein 2 (CDX2), but in contrast to the mouse in a GATA-binding protein 3 (GATA3)-independent manner.

WHAT IS KNOWN ALREADY

Tead4 is one of the earliest transcription factors expressed during mouse embryo preimplantation development and is required for the expression of TE-associated genes. Functional knock-out studies in mouse, inactivating Tead4 by site-specific recombination, have shown that Tead4-targeted embryos have compromised development and expression of the TE-specific Cdx2 and Gata3 is downregulated. Cdx2 and Gata3 act in parallel pathways downstream of Tead4 to induce successful TE differentiation. Downstream loss of Cdx2 expression, compromises TE differentiation and subsequent blastocoel formation and leads to the ectopic expression of inner cell mass (ICM) genes, including POU Class 5 homeobox 1 (Pou5f1) and SRY-box transcription factor (Sox2). Cdx2 is a more potent regulator of TE fate in mouse as loss of Cdx2 expression induces more severe phenotypes compared with loss of Gata3 expression. The role of TEAD4 and its downstream effectors during human preimplantation embryo development has not been investigated yet.

STUDY DESIGN, SIZE, DURATION

The clustered regularly interspaced short palindromic repeats-clustered regularly interspaced short palindromic repeats (CRISPR)-associated genes (CRISPR-Cas9) system was first introduced in pronuclei (PN)-stage mouse zygotes aiming to identify a guide RNA (gRNA), yielding high editing efficiency and effective disruption of the Tead4 locus. Three guides were tested (gRNA1-3), each time targeting a distinct region of Exon 2 of Tead4. The effects of targeting on developmental capacity were studied in Tead4-targeted embryos (n = 164-summarized data from gRNA1-3) and were compared with two control groups; sham-injected embryos (n = 26) and non-injected media-control embryos (n = 51). The editing efficiency was determined by next-generation sequencing (NGS). In total, n = 55 (summarized data from gRNA1-3) targeted mouse embryos were analysed by NGS. Immunofluorescence analysis to confirm successful targeting by gRNA1 was performed in Tead4-targeted embryos, and non-injected media-control embryos. The downregulation of secondary TE-associated markers Cdx2 and Gata3 was used as an indirect confirmation of successful Tead4-targeting (previously shown to be expressed downstream of Tead4). Additional groups of gRNA1 Tead4-targeted (n = 45) and media control (n = 36) embryos were cultured for an extended period of 8.5 days, to further assess the developmental capacity of the Tead4-targeted group to develop beyond implantation stages. Following the mouse investigation, human metaphase-II (MII) oocytes obtained by IVM were microinjected with gRNA-Cas9 during ICSI (n = 74) to target TEAD4 or used as media-control (n = 33). The editing efficiency was successfully assessed in n = 25 TEAD4-targeted human embryos. Finally, immunofluorescence analysis for TEAD4, CDX2, GATA3 and the ICM marker SOX2 was performed in TEAD4-targeted (n = 10) and non-injected media-control embryos (n = 29).

PARTICIPANTS/MATERIALS, SETTING, METHODS

A ribonucleoprotein complex consisting of a gRNA-Cas9 mixture, designed to target Exon 2 of Tead4/TEAD4, was microinjected in mouse PN stage zygotes or human IVM MII oocytes along with sperm. Generated embryos were cultured in vitro for 4 days in mouse or 6.5 days in human. In mouse, an additional group of Tead4-targeted and media-control embryos was cultured in vitro for an extended period of 8.5 days. Embryonic development and morphology were assessed daily, during culture in vitro of mouse and human embryos and was followed by a detailed scoring at late blastocyst stage. Targeting efficiency following gRNA-Cas9 introduction was assessed via immunostaining and NGS analysis.

MAIN RESULTS AND THE ROLE OF CHANCE

NGS analysis of the Tead4-targeted locus revealed very high editing efficiencies for all three guides, with 100% of the mouse embryos (55 out of 55) carrying genetic modifications resulting from CRISPR-Cas9 genome editing. More specifically, 65.22% (15 out 23) of the PN zygotes microinjected with gRNA1-Cas9, which exhibited the highest efficiency, carried exclusively mutated alleles. The developmental capacity of targeted embryos was significantly reduced (data from gRNA1), as 44.17% of the embryos arrested at the morula stage (2.5 days post coitum), coincident with the initiation of TE lineage differentiation, compared with 8.51% in control and 12.50% in sham control groups. High-quality blastocyst formation rates (Grade 3) were 8.97% in the gRNA1-targeted group, compared with 87.23% in the media-control and 87.50% in the sham group. Immunofluorescence analysis in targeted embryos confirmed downregulation of Tead4, Cdx2, and Gata3 expression, which resulted from successful targeting of the Tead4 locus. Tead4-targeted mouse embryos stained positive for the ICM markers Pou5f1 and Sox2, indicating that expression of ICM lineage markers is not affected. Tead4-targeted embryos were able to cavitate and form a blastocoel without being able to hatch. Extended embryo culture following zona pellucida removal, revealed that the targeted embryos can attach and form egg-cylinder-like structures in the absence of trophoblast giant cells. In human embryos, Exon 2 of TEAD4 was successfully targeted by CRISPR-Cas9 (n = 74). In total, 25 embryos from various developmental stages were analysed by NGS and 96.00% (24 out of 25) of the embryos carried genetic modifications because of gRNA-Cas9 editing. In the subgroup of the 24 edited embryos, 17 (70.83%) carried only mutant alleles and 11 out of these 17 (64.70%) carried exclusively frameshift mutations. Six out of 11 embryos reached the blastocyst stage. In contrast to mice, human-targeted embryos formed blastocysts at a rate (25.00%) that did not differ significantly from the control group (23.81%). However, blastocyst morphology and TE quality were significantly compromised following TEAD4-targeting, showing grade C TE scores, with TE containing very few cells. Immunofluorescence analysis of TEAD4-targeted embryos (n = 10) confirmed successful editing by the complete absence of TEAD4 and its downstream TE marker CDX2, but the embryos generated retained expression of GATA3, which is in contrast to what we have observed and has previously been reported in mouse. In this regard, our results indicate that GATA3 acts in parallel with TEAD4/CDX2 towards TE differentiation in human.

LARGE SCALE DATA

N/A.

LIMITATIONS, REASONS FOR CAUTION

CRISPR-Cas9 germline genome editing, in some cases, induces mosaic genotypes. These genotypes are a result of inefficient and delayed editing, and complicate the phenotypic analysis and developmental assessment of the injected embryos. We cannot exclude the possibility that the observed differences between mouse and human are the result of variable effects triggered by the culture conditions, which were however similar for both mouse and human embryos in this study. Furthermore, this study utilized human oocytes obtained by IVM, which may not fully recapitulate the developmental behaviour of in vivo matured oocytes.

WIDER IMPLICATIONS OF THE FINDINGS

Elucidation of the evolutionary conservation of molecular mechanisms that regulate the differentiation and formation of the trophoblast lineage can give us fundamental insights into early implantation failure, which accounts for ∼15% of human conceptions.

STUDY FUNDING/COMPETING INTEREST(S)

The research was funded by the FWO-Vlaanderen (Flemish fund for scientific research, Grant no. G051516N), and Hercules funding (FWO.HMZ.2016.00.02.01) and Ghent University (BOF.BAS.2018.0018.01). G.C. is supported by FWO-Vlaanderen (Flemish fund for scientific research, Grant no. 11L8822N). A.B. is supported by FWO-Vlaanderen (Flemish fund for scientific research, Grant no. 1298722 N). We further thank Ferring Pharmaceuticals (Aalst, Belgium) for their unrestricted educational grant. The authors declare no competing interests.

TRIAL REGISTRATION NUMBER

N/A.

Extensive co-binding and rapid redistribution of NANOG and GATA6 during emergence of divergent lineages

Fate-determining transcription factors (TFs) can promote lineage-restricted transcriptional programs from common progenitor states. The inner cell mass (ICM) of mouse blastocysts co-expresses the TFs NANOG and GATA6, which drive the bifurcation of the ICM into either the epiblast (Epi) or the primitive endoderm (PrE), respectively. Here, we induce GATA6 in embryonic stem cells-that also express NANOG-to characterize how a state of co-expression of opposing TFs resolves into divergent lineages. Surprisingly, we find that GATA6 and NANOG co-bind at the vast majority of Epi and PrE enhancers, a phenomenon we also observe in blastocysts. The co-bound state is followed by eviction and repression of Epi TFs, and quick remodeling of chromatin and enhancer-promoter contacts thus establishing the PrE lineage while repressing the Epi fate. We propose that co-binding of GATA6 and NANOG at shared enhancers maintains ICM plasticity and promotes the rapid establishment of Epi- and PrE-specific transcriptional programs.

NANOG initiates epiblast fate through the coordination of pluripotency genes expression

The epiblast is the source of all mammalian embryonic tissues and of pluripotent embryonic stem cells. It differentiates alongside the primitive endoderm in a “salt and pepper” pattern from inner cell mass (ICM) progenitors during the preimplantation stages through the activity of NANOG, GATA6 and the FGF pathway. When and how epiblast lineage specification is initiated is still unclear. Here, we show that the coordinated expression of pluripotency markers defines epiblast identity. Conversely, ICM progenitor cells display random cell-to-cell variability in expression of various pluripotency markers, remarkably dissimilar from the epiblast signature and independently from NANOG, GATA6 and FGF activities. Coordination of pluripotency markers expression fails in Nanog and Gata6 double KO (DKO) embryos. Collectively, our data suggest that NANOG triggers epiblast specification by ensuring the coordinated expression of pluripotency markers in a subset of cells, implying a stochastic mechanism. These features are likely conserved, as suggested by analysis of human embryos.

Pluripotent epiblast cells segregate from primitive endoderm in the blastocyst inner cell mass (ICM). Here the authors show that mosaic epiblast differentiation during mouse and human preimplantation development initiates stochastically in ICM progenitors, independently of the FGF pathway, and requires NANOG activity

Dynamical landscapes of cell fate decisions

The generation of cellular diversity during development involves differentiating cells transitioning between discrete cell states. In the 1940s, the developmental biologist Conrad Waddington introduced a landscape metaphor to describe this process. The developmental path of a cell was pictured as a ball rolling through a terrain of branching valleys with cell fate decisions represented by the branch points at which the ball decides between one of two available valleys. Here we discuss progress in constructing quantitative dynamical models inspired by this view of cellular differentiation. We describe a framework based on catastrophe theory and dynamical systems methods that provides the foundations for quantitative geometric models of cellular differentiation. These models can be fit to experimental data and used to make quantitative predictions about cellular differentiation. The theory indicates that cell fate decisions can be described by a small number of decision structures, such that there are only two distinct ways in which cells make a binary choice between one of two fates. We discuss the biological relevance of these mechanisms and suggest the approach is broadly applicable for the quantitative analysis of differentiation dynamics and for determining principles of developmental decisions.

Multi-synchronization and other patterns of multi-rhythmicity in oscillatory biological systems

While experimental and theoretical studies have established the prevalence of rhythmic behaviour at all levels of biological organization, less common is the coexistence between multiple oscillatory regimes (multi-rhythmicity), which has been predicted by a variety of models for biological oscillators. The phenomenon of multi-rhythmicity involves, most commonly, the coexistence between two (birhythmicity) or three (trirhythmicity) distinct regimes of self-sustained oscillations. Birhythmicity has been observed experimentally in a few chemical reactions and in biological examples pertaining to cardiac cell physiology, neurobiology, human voice patterns and ecology. The present study consists of two parts. We first review the mechanisms underlying multi-rhythmicity in models for biochemical and cellular oscillations in which the phenomenon was investigated over the years. In the second part, we focus on the coupling of the cell cycle and the circadian clock and show how an additional source of multi-rhythmicity arises from the bidirectional coupling of these two cellular oscillators. Upon bidirectional coupling, the two oscillatory networks generally synchronize in a unique manner characterized by a single, common period. In some conditions, however, the two oscillators may synchronize in two or three different ways characterized by distinct waveforms and periods. We refer to this type of multi-rhythmicity as ‘multi-synchronization’.

The role of FGF-4 and FGFR-2 on preimplantation embryo development in experimental maternal diabetes

OBJECTIVE

Diabetes mellitus can cause spontaneous abortion, neonatal diseases, congenital malformations, and death. There are many studies related to the damage of in vitro hyperglycemia on embryogenesis in literature, but not enough studies on in vivo hyperglycemia effects on embryogenesis. Fibroblast growth factor (FGF) molecules play an essential role in pre-implantation embryo development and diabetes pathogenesis. In our study, we researched whether FGF-4 and FGFR-2 were playing a role in maternal diabetes' effects on embryo development.

MATERIAL AND METHODS

Thirty adult virgin female BALB/c mice were randomly divided into two groups: control and diabetic. The experimental diabetes model was generated by streptozotocin (55 mg/kg, once, intraperitoneally). The control and the diabetic group were mated. Embryos were collected at the morula and blastocyte stages corresponding to the third and fourth days of pregnancy. Embryo's FGF-4 and FGFR-2 molecules were evaluated by their immunofluorescence staining and immunoreactivity score.

RESULT

The results clearly showed that the FGF-4 and FGFR-2 immunofluorescence reactivity was higher in the diabetes group.

CONCLUSION

We concluded that FGF-4 and FGFR-2 overexpression might impair mouse pre-implantation embryo development in maternal diabetes and suggest investigating whether they have crucial effects on human embryo development and infertility in maternal diabetes.

Transcriptional and epigenetic control of early life cell fate decisions

PURPOSE OF REVIEW

Global epigenetic reprogramming of the parental genomes after fertilization ensures the establishment of genome organization permissive for cell specialization and differentiation during development. In this review, we highlight selected, well-characterized relationships between epigenetic factors and transcriptional cell fate regulators during the initial stages of mouse development.

RECENT FINDINGS

Blastomeres of the mouse embryo are characterized by atypical and dynamic histone modification arrangements, noncoding RNAs and DNA methylation profiles. Moreover, asymmetries in epigenomic patterning between embryonic cells arise as early as the first cleavage, with potentially instructive roles during the first lineage allocations in the mouse embryo. Although it is widely appreciated that transcription factors and developmental signaling pathways play a crucial role in cell fate specification at the onset of development, it is increasingly clear that their function is tightly connected to the underlying epigenetic status of the embryonic cells in which they act.

SUMMARY

Findings on the interplay between genetic, epigenetic and environmental factors during reprogramming and differentiation in the embryo are crucial for understanding the molecular underpinnings of disease processes, particularly tumorigenesis, which is characterized by global epigenetic rewiring and progressive loss of cellular identity.

Gaining confidence in inferred networks

Network inference is a notoriously challenging problem. Inferred networks are associated with high uncertainty and likely riddled with false positive and false negative interactions. Especially for biological networks we do not have good ways of judging the performance of inference methods against real networks, and instead we often rely solely on the performance against simulated data. Gaining confidence in networks inferred from real data nevertheless thus requires establishing reliable validation methods. Here, we argue that the expectation of mixing patterns in biological networks such as gene regulatory networks offers a reasonable starting point: interactions are more likely to occur between nodes with similar biological functions. We can quantify this behaviour using the assortativity coefficient, and here we show that the resulting heuristic, functional assortativity, offers a reliable and informative route for comparing different inference algorithms.

Synthetic multistability in mammalian cells

In multicellular organisms, gene regulatory circuits generate thousands of molecularly distinct, mitotically heritable states through the property of multistability. Designing synthetic multistable circuits would provide insight into natural cell fate control circuit architectures and would allow engineering of multicellular programs that require interactions among distinct cell types. We created MultiFate, a naturally inspired, synthetic circuit that supports long-term, controllable, and expandable multistability in mammalian cells. MultiFate uses engineered zinc finger transcription factors that transcriptionally self-activate as homodimers and mutually inhibit one another through heterodimerization. Using a model-based design, we engineered MultiFate circuits that generate as many as seven states, each stable for at least 18 days. MultiFate permits controlled state switching and modulation of state stability through external inputs and can be expanded with additional transcription factors. These results provide a foundation for engineering multicellular behaviors in mammalian cells.

Specification and role of extraembryonic endoderm lineages in the periimplantation mouse embryo

During mammalian embryo development, the correct formation of the first extraembryonic endoderm lineages is fundamental for successful development. In the periimplantation blastocyst, the primitive endoderm (PrE) is formed, which gives rise to the parietal endoderm (PE) and visceral endoderm (VE) during further developmental stages. These PrE-derived lineages show significant differences in both their formation and roles. Whereas differentiation of the PE as a migratory lineage has been suggested to represent the first epithelial-to-mesenchymal transition (EMT) in development, organisation of the epithelial VE is of utmost importance for the correct axis definition and patterning of the embryo. Despite sharing a common origin, the striking differences between the VE and PE are indicative of their distinct roles in early development. However, there is a significant disparity in the current knowledge of each lineage, which reflects the need for a deeper understanding of their respective specification processes. In this review, we will discuss the origin and maturation of the PrE, PE, and VE during the periimplantation period using the mouse model as an example. Additionally, we consider the latest findings regarding the role of the PrE-derived lineages and early embryo morphogenesis, as obtained from the most recent in vitro models.

How a cell decides its own fate: a single-cell view of molecular mechanisms and dynamics of cell-type specification

On its path from a fertilized egg to one of the many cell types in a multicellular organism, a cell turns the blank canvas of its early embryonic state into a molecular profile fine-tuned to achieve a vital organismal function. This remarkable transformation emerges from the interplay between dynamically changing external signals, the cell's internal, variable state, and tremendously complex molecular machinery; we are only beginning to understand. Recently developed single-cell omics techniques have started to provide an unprecedented, comprehensive view of the molecular changes during cell-type specification and promise to reveal the underlying gene regulatory mechanism. The exponentially increasing amount of quantitative molecular data being created at the moment is slated to inform predictive, mathematical models. Such models can suggest novel ways to manipulate cell types experimentally, which has important biomedical applications. This review is meant to give the reader a starting point to participate in this exciting phase of molecular developmental biology. We first introduce some of the principal molecular players involved in cell-type specification and discuss the important organizing ability of biomolecular condensates, which has been discovered recently. We then review some of the most important single-cell omics methods and relevant findings they produced. We devote special attention to the dynamics of the molecular changes and discuss methods to measure them, most importantly lineage tracing. Finally, we introduce a conceptual framework that connects all molecular agents in a mathematical model and helps us make sense of the experimental data.

Capturing Pluripotency and Beyond

During the development of a multicellular organism, the specification of different cell lineages originates in a small group of pluripotent cells, the epiblasts, formed in the preimplantation embryo. The pluripotent epiblast is protected from premature differentiation until exposure to inductive cues in strictly controlled spatially and temporally organized patterns guiding fetus formation. Epiblasts cultured in vitro are embryonic stem cells (ESCs), which recapitulate the self-renewal and lineage specification properties of their endogenous counterparts. The characteristics of totipotency, although less understood than pluripotency, are becoming clearer. Recent studies have shown that a minor ESC subpopulation exhibits expanded developmental potential beyond pluripotency, displaying a characteristic reminiscent of two-cell embryo blastomeres (2CLCs). In addition, reprogramming both mouse and human ESCs in defined media can produce expanded/extended pluripotent stem cells (EPSCs) similar to but different from 2CLCs. Further, the molecular roadmaps driving the transition of various potency states have been clarified. These recent key findings will allow us to understand eutherian mammalian development by comparing the underlying differences between potency network components during development. Using the mouse as a paradigm and recent progress in human PSCs, we review the epiblast's identity acquisition during embryogenesis and their ESC counterparts regarding their pluripotent fates and beyond.