Gene expression dynamics underlying cell fate emergence in 2D micropatterned human embryonic stem cell gastruloids

Human embryonic stem cells cultured in 2D micropatterns with BMP4 differentiate into a radial arrangement of germ layers and extraembryonic cells. Single-cell transcriptomes demonstrate generation of cell types transcriptionally similar to their in vivo counterparts in Carnegie stage 7 human gastrula. Time-course analyses indicate sequential differentiation, where the epiblast arises by 12 h between the prospective ectoderm in the center and the cells initiating differentiation toward extraembryonic fates at the edge. Extraembryonic and mesendoderm precursors arise from the epiblast by 24 h, while nascent mesoderm, endoderm, and primordial germ cell-like cells form by 44 h. Dynamic changes in transcripts encoding signaling components support a BMP, WNT, and Nodal hierarchy underlying germ-layer specification conserved across mammals, and FGF and HIPPO pathways being active throughout differentiation. This work also provides a resource for mining genes and pathways expressed in a stereotyped 2D gastruloid model, common with other species or unique to human gastrulation

Dissecting cell identity via network inference and in silico gene perturbation

Cell identity is governed by the complex regulation of gene expression, represented as gene-regulatory network

Stability analysis of non-autonomous reaction-diffusion systems: the effects of growing domains

By using asymptotic theory, we generalise the Turing diffusively-driven instability conditions for reaction-diffusion systems with slow, isotropic domain growth. There are two fundamental biological differences between the Turing conditions on fixed and growing domains, namely: (i) we need not enforce cross nor pure kinetic conditions and (ii) the restriction to activator-inhibitor kinetics to induce pattern formation on a growing biological system is no longer a requirement. Our theoretical findings are confirmed and reinforced by numerical simulations for the special cases of isotropic linear, exponential and logistic growth profiles. In particular we illustrate an example of a reaction-diffusion system which cannot exhibit a diffusively-driven instability on a fixed domain but is unstable in the presence of slow growth

High-resolution transcriptional and morphogenetic profiling of cells from micropatterned human ESC gastruloid cultures

During mammalian gastrulation, germ layers arise and are shaped into the body plan while extraembryonic layers sustain the embryo. Human embryonic stem cells, cultured with BMP4 on extracellular matrix micro-discs, reproducibly differentiate into gastruloids, expressing markers of germ layers and extraembryonic cells in radial arrangement. Using single-cell RNA sequencing and cross-species comparisons with mouse, cynomolgus monkey gastrulae, and post-implantation human embryos, we reveal that gastruloids contain cells transcriptionally similar to epiblast, ectoderm, mesoderm, endoderm, primordial germ cells, trophectoderm, and amnion. Upon gastruloid dissociation, single cells reseeded onto micro-discs were motile and aggregated with the same but segregated from distinct cell types. Ectodermal cells segregated from endodermal and extraembryonic but mixed with mesodermal cells. Our work demonstrates that the gastruloid system models primate-specific features of embryogenesis, and that gastruloid cells exhibit evolutionarily conserved sorting behaviors. This work generates a resource for transcriptomes of human extraembryonic and embryonic germ layers differentiated in a stereotyped arrangement

Atp7a determines a hierarchy of copper metabolism essential for notochord development

SummaryThe critical developmental and genetic requirements of copper metabolism during embryogenesis are unknown. Utilizing a chemical genetic screen in zebrafish, we identified small molecules that perturb copper homeostasis. Our findings reveal a role for copper in notochord formation and demonstrate a hierarchy of copper metabolism within the embryo. To elucidate these observations, we interrogated a genetic screen for embryos phenocopied by copper deficiency, identifying calamity, a mutant defective in the zebrafish ortholog of the Menkes disease gene (atp7a). Copper metabolism in calamity is restored by human ATP7A, and transplantation experiments reveal that atp7a functions cell autonomously, findings with important therapeutic implications. The gene dosage of atp7a determines the sensitivity to copper deprivation, revealing that the observed developmental hierarchy of copper metabolism is informed by specific genetic factors. Our data provide insight into the developmental pathophysiology of copper metabolism and suggest that suboptimal copper metabolism may contribute to birth defects

Gα12/13 regulate epiboly by inhibiting E-cadherin activity and modulating the actin cytoskeleton

Epiboly spreads and thins the blastoderm over the yolk cell during zebrafish gastrulation, and involves coordinated movements of several cell layers. Although recent studies have begun to elucidate the processes that underlie these epibolic movements, the cellular and molecular mechanisms involved remain to be fully defined. Here, we show that gastrulae with altered Gα12/13 signaling display delayed epibolic movement of the deep cells, abnormal movement of dorsal forerunner cells, and dissociation of cells from the blastoderm, phenocopying e-cadherin mutants. Biochemical and genetic studies indicate that Gα12/13 regulate epiboly, in part by associating with the cytoplasmic terminus of E-cadherin, and thereby inhibiting E-cadherin activity and cell adhesion. Furthermore, we demonstrate that Gα12/13 modulate epibolic movements of the enveloping layer by regulating actin cytoskeleton organization through a RhoGEF/Rho-dependent pathway. These results provide the first in vivo evidence that Gα12/13 regulate epiboly through two distinct mechanisms: limiting E-cadherin activity and modulating the organization of the actin cytoskeleton

Characterization of Turing diffusion-driven instability on evolving domains

In this paper we establish a general theoretical framework for Turing diffusion-driven instability for reaction-diffusion systems on time-dependent evolving domains. The main result is that Turing diffusion-driven instability for reaction-diffusion systems on evolving domains is characterised by Lyapunov exponents of the evolution family associated with the linearised system (obtained by linearising the original system along a spatially independent solution). This framework allows for the inclusion of the analysis of the long-time behavior of the solutions of reaction-diffusion systems. Applications to two special types of evolving domains are considered: (i) time-dependent domains which evolve to a final limiting fixed domain and (ii) time-dependent domains which are eventually time periodic. Reaction-diffusion systems have been widely proposed as plausible mechanisms for pattern formation in morphogenesis

Aberrant behaviours of reaction diffusion self-organisation models on growing domains in the presence of gene expression time delays

Turing’s pattern formation mechanism exhibits sensitivity to the details of the initial conditions suggesting that, in isolation, it cannot robustly generate pattern within noisy biological environments. Nonetheless, secondary aspects of developmental self-organisation, such as a growing domain, have been shown to ameliorate this aberrant model behaviour. Furthermore, while in-situ hybridisation reveals the presence of gene expression in developmental processes, the influence of such dynamics on Turing’s model has received limited attention. Here, we novelly focus on the Gierer–Meinhardt reaction diffusion system considering delays due the time taken for gene expression, while incorporating a number of different domain growth profiles to further explore the influence and interplay of domain growth and gene expression on Turing’s mechanism. We find extensive pathological model behaviour, exhibiting one or more of the following: temporal oscillations with no spatial structure, a failure of the Turing instability and an extreme sensitivity to the initial conditions, the growth profile and the duration of gene expression. This deviant behaviour is even more severe than observed in previous studies of Schnakenberg kinetics on exponentially growing domains in the presence of gene expression (Gaffney and Monk in Bull. Math. Biol. 68:99–130, 2006). Our results emphasise that gene expression dynamics induce unrealistic behaviour in Turing’s model for multiple choices of kinetics and thus such aberrant modelling predictions are likely to be generic. They also highlight that domain growth can no longer ameliorate the excessive sensitivity of Turing’s mechanism in the presence of gene expression time delays. The above, extensive, pathologies suggest that, in the presence of gene expression, Turing’s mechanism would generally require a novel and extensive secondary mechanism to control reaction diffusion patterning

Zebrafish TRPA1 Channels Are Required for Chemosensation But Not for Thermosensation or Mechanosensory Hair Cell Function

Transient receptor potential (TRP) ion channels have been implicated in detecting chemical, thermal, and mechanical stimuli in organisms ranging from mammals to Caenorhabditis elegans. It is well established that TRPA1 detects and mediates behavioral responses to chemical irritants. However, the role of TRPA1 in detecting thermal and mechanical stimuli is controversial. To further clarify the functions of TRPA1 channels in vertebrates, we analyzed their roles in zebrafish. The two zebrafish TRPA1 paralogs are expressed in sensory neurons and are activated by several chemical irritants in vitro. High-throughput behavioral analyses of trpa1a and trpa1b mutant larvae indicate that TRPA1b is necessary for behavioral responses to these chemical irritants. However, TRPA1 paralogs are not required for behavioral responses to temperature changes or for mechanosensory hair cell function in the inner ear or lateral line. These results support a role for zebrafish TRPA1 in chemical but not thermal or mechanical sensing, and establish a high-throughput system to identify genes and small molecules that modulate chemosensation, thermosensation, and mechanosensation