Bistability in a model of mesoderm and anterior mesendoderm specification in Xenopus laevis

A gene regulatory program controlling early Xenopus mesendoderm formation: Network conservation and motifs

Germ layer formation is among the earliest differentiation events in metazoan embryos. In triploblasts, three germ layers are formed, among which the endoderm gives rise to the epithelial lining of the gut tube and associated organs including the liver, pancreas and lungs. In frogs (Xenopus), where early germ layer formation has been studied extensively, the process of endoderm specification involves the interplay of dozens of transcription factors. Here, we review the interactions between these factors, summarized in a transcriptional gene regulatory network (GRN). We highlight regulatory connections conserved between frog, fish, mouse, and human endodermal lineages. Especially prominent is the conserved role and regulatory targets of the Nodal signaling pathway and the T-box transcription factors, Vegt and Eomes. Additionally, we highlight network topologies and motifs, and speculate on their possible roles in development.

Redeployment of germ layers related TFs shows regionalized expression during two non-embryonic developments

In all non-vertebrate metazoan phyla, species that evolved non-embryonic developmental pathways as means of propagation or regeneration can be found. In this context, new bodies arise through asexual reproduction processes (such as budding) or whole body regeneration, that lack the familiar temporal and spatial cues classically associated with embryogenesis, like maternal determinants, or gastrulation. The molecular mechanisms underlying those non-embryonic developments (i.e., regeneration and asexual reproduction), and their relationship to those deployed during embryogenesis are poorly understood. We have addressed this question in the colonial ascidian Botryllus schlosseri, which undergoes an asexual reproductive process via palleal budding (PB), as well as a whole body regeneration by vascular budding (VB). We identified early regenerative structures during VB and then followed the fate of differentiating tissues during both non-embryonic developments (PB and VB) by monitoring the expression of genes known to play key functions in germ layer specification with well conserved expression patterns in solitary ascidian embryogenesis. The expression patterns of FoxA1, GATAa, GATAb, Otx, Bra, Gsc and Tbx2/3 were analysed during both PB and VB. We found that the majority of these transcription factors were expressed during both non-embryonic developmental processes, revealing a regionalization of the palleal and vascular buds. Knockdown of GATAa by siRNA in palleal buds confirmed that preventing the correct development of one of these regions blocks further tissue specification. Our results indicate that during both normal and injury-induced budding, a similar alternative developmental program operates via early commitment of epithelial regions.

Multicellular Mathematical Modelling of Mesendoderm Formation in Amphibians

The earliest cell fate decisions in a developing embryo are those associated with establishing the germ layers. The specification of the mesoderm and endoderm is of particular interest as the mesoderm is induced from the endoderm, potentially from an underlying bipotential group of cells, the mesendoderm. Mesendoderm formation has been well studied in an amphibian model frog, Xenopus laevis, and its formation is driven by a gene regulatory network (GRN) induced by maternal factors deposited in the egg. We have recently demonstrated that the axolotl, a urodele amphibian, utilises a different topology in its GRN to specify the mesendoderm. In this paper, we develop spatially structured mathematical models of the GRNs governing mesendoderm formation in a line of cells. We explore several versions of the model of mesendoderm formation in both Xenopus and the axolotl, incorporating the key differences between these two systems. Model simulations are able to reproduce known experimental data, such as Nodal expression domains in Xenopus, and also make predictions about how the positional information derived from maternal factors may be interpreted to drive cell fate decisions. We find that whilst cell-cell signalling plays a minor role in Xenopus, it is crucial for correct patterning domains in axolotl.

Two different network topologies yield bistability in models of mesoderm and anterior mesendoderm specification in amphibians

Understanding the Gene Regulatory Networks (GRNs) that underlie development is a major question for systems biology. The establishment of the germ layers is amongst the earliest events of development and has been characterised in numerous model systems. The establishment of the mesoderm is best characterised in the frog Xenopus laevis and has been well studied both experimentally and mathematically. However, the Xenopus network has significant differences from that in mouse and humans, including the presence of multiple copies of two key genes in the network, Mix and Nodal. The axolotl, a urodele amphibian, provides a model with all the benefits of amphibian embryology but crucially only a single Mix and Nodal gene required for the specification of the mesoderm. Remarkably, the number of genes within the network is not the only difference. The interaction between Mix and Brachyury, two transcription factors involved in the establishment of the endoderm and mesoderm respectively, is not conserved. While Mix represses Brachyury in Xenopus, it activates Brachyury in axolotl. Thus, whilst the topology of the networks in the two species differs, both are able to form mesoderm and endoderm in vivo. Based on current knowledge of the structure of the mesendoderm GRN we develop deterministic models that describe the time evolution of transcription factors in a single axolotl cell and compare numerical simulations with previous results from Xenopus. The models are shown to have stable steady states corresponding to mesoderm and anterior mesendoderm, with the in vitro model showing how the concentration of Activin can determine cell fate, while the in vivo model shows that β-catenin concentration can determine cell fate. Moreover, our analysis suggests that additional components must be important in the axolotl network in the specification of the full range of tissues.

•

We present models of mesendoderm specification in the urodele amphibian, the axolotl.

•

in vitro and in vivo models are simulated and compared with experimental data.

•

The model topology differs from that of the anuran amphibian, Xenopus laevis.

•

Steady states representing mesoderm and anterior mesendoderm are found in both models.

•

Both the axolotl and Xenopus topologies can account for similar qualitative data.

Wave pinning and spatial patterning in a mathematical model of Antivin/Lefty-Nodal signalling

Nodal signals are key regulators of mesoderm and endoderm development in vertebrate embryos. It has been observed experimentally that in Xenopus embryos the spatial range of Nodal signals is restricted by the signal Antivin (also known as Lefty). Nodal signals can activate both Nodal and Antivin, whereas Antivin is thought to antagonise Nodal by binding either directly to it or to its receptor. In this paper we develop a mathematical model of this signalling network in a line of cells. We consider the heterodimer and receptor-mediated inhibition mechanisms separately and find that, in both cases, the restriction by Antivin to the range of Nodal signals corresponds to wave pinning in the model. Our analysis indicates that, provided Antivin diffuses faster than Nodal, either mechanism can robustly account for the experimental data. We argue that, in the case of Xenopus development, it is wave pinning, rather than Turing-type patterning, that is underlying Nodal-Antivin dynamics. This leads to several experimentally testable predictions, which are discussed. Furthermore, for heterodimer-mediated inhibition to prevent waves of Nodal expression from propagating, the Nodal-Antivin complex must be turned over, and diffusivity of the complex must be negligible. In the absence of molecular mechanisms regulating these, we suggest that Antivin restricts Nodal signals via receptor-mediated, and not heterodimer-mediated, inhibition.