A shark-inspired general model of tooth morphogenesis unveils developmental asymmetries in phenotype transitions.

Developmental complexity stemming from the dynamic interplay between genetic and biomechanic factors canalizes the ways genotypes and phenotypes can change in evolution. As a paradigmatic system, we explore how changes in developmental factors generate typical tooth shape transitions. Since tooth development has mainly been researched in mammals, we contribute to a more general understanding by studying the development of tooth diversity in sharks. To this end, we build a general, but realistic, mathematical model of odontogenesis. We show that it reproduces key shark-specific features of tooth development as well as real tooth shape variation in small-spotted catsharks Scyliorhinus canicula. We validate our model by comparison with experiments in vivo. Strikingly, we observe that developmental transitions between tooth shapes tend to be highly degenerate, even for complex phenotypes. We also discover that the sets of developmental parameters involved in tooth shape transitions tend to depend asymmetrically on the direction of that transition. Together, our findings provide a valuable base for furthering our understanding of how developmental changes can lead to both adaptive phenotypic change and trait convergence in complex, phenotypically highly diverse, structures.

Turing's turtles all the way down: A conserved role of EDAR in the carapacial ridge suggests a deep homology of prepatterns across ectodermal appendages.

The scutes of the turtle shell are epidermal shields that begin their formation during the early stages of shell development. Like other skin appendages, turtle scutes are hypothesized to be patterned by reaction-diffusion systems. We have previously established ex vivo and in silico systems to study these mechanisms experimentally and have further shown that mathematical models can explain the dynamics of the induction of turtle scute primordia and the generation of final scute architecture. Using these foundations, we expand our current knowledge and test the roles of ectodysplasin and activin signaling in the development of turtle scutes. We find that these molecules play important roles in the prepatterning of scute primordia along the carapacial ridge and show that blocking Edar signaling may lead to a complete loss of marginal scute primordia. We show that it is possible to reproduce these observations using simple mathematical modeling, thereby suggesting a stabilizing role for ectodysplasin within the reaction-diffusion mechanisms. Finally, we argue that our findings further entrench turtle scutes within a class of developmental systems composed of hierarchically nested reaction-diffusion mechanisms, which is conserved across ectodermal organs.

Beyond genotype-phenotype maps: Toward a phenotype-centered perspective on evolution.

Evolutionary biology is paying increasing attention to the mechanisms that enable phenotypic plasticity, evolvability, and extra-genetic inheritance. Yet, there is a concern that these phenomena remain insufficiently integrated within evolutionary theory. Understanding their evolutionary implications would require focusing on phenotypes and their variation, but this does not always fit well with the prevalent genetic representation of evolution that screens off developmental mechanisms. Here, we instead use development as a starting point, and represent it in a way that allows genetic, environmental and epigenetic sources of phenotypic variation to be independent. We show why this representation helps to understand the evolutionary consequences of both genetic and non-genetic phenotype determinants, and discuss how this approach can instigate future areas of empirical and theoretical research.

Correction: On the evolution and development of morphological complexity: A view from gene regulatory networks.

[This corrects the article DOI: 10.1371/journal.pcbi.1008570.].

On the evolution and development of morphological complexity: A view from gene regulatory networks.

How does morphological complexity evolve? This study suggests that the likelihood of mutations increasing phenotypic complexity becomes smaller when the phenotype itself is complex. In addition, the complexity of the genotype-phenotype map (GPM) also increases with the phenotypic complexity. We show that complex GPMs and the above mutational asymmetry are inevitable consequences of how genes need to be wired in order to build complex and robust phenotypes during development. We randomly wired genes and cell behaviors into networks in EmbryoMaker. EmbryoMaker is a mathematical model of development that can simulate any gene network, all animal cell behaviors (division, adhesion, apoptosis, etc.), cell signaling, cell and tissues biophysics, and the regulation of those behaviors by gene products. Through EmbryoMaker we simulated how each random network regulates development and the resulting morphology (i.e. a specific distribution of cells and gene expression in 3D). This way we obtained a zoo of possible 3D morphologies. Real gene networks are not random, but a random search allows a relatively unbiased exploration of what is needed to develop complex robust morphologies. Compared to the networks leading to simple morphologies, the networks leading to complex morphologies have the following in common: 1) They are rarer; 2) They need to be finely tuned; 3) Mutations in them tend to decrease morphological complexity; 4) They are less robust to noise; and 5) They have more complex GPMs. These results imply that, when complexity evolves, it does so at a progressively decreasing rate over generations. This is because as morphological complexity increases, the likelihood of mutations increasing complexity decreases, morphologies become less robust to noise, and the GPM becomes more complex. We find some properties in common, but also some important differences, with non-developmental GPM models (e.g. RNA, protein and gene networks in single cells).

Cell signaling stabilizes morphogenesis against noise.

Embryonic development involves gene networks, extracellular signaling, cell behaviors (cell division, adhesion, etc.) and mechanical interactions. How should these be coordinated to lead to complex and robust morphologies? To explore this question, we randomly wired genes and cell behaviors into a huge number of networks in EmbryoMaker. EmbryoMaker is a computational model of animal development that simulates how the 3D positions of cells, i.e. morphology, change over time due to such networks. We found that any gene network can lead to complex morphologies if this activates cell behaviors over large regions of the embryo. Importantly, however, for such complex morphologies to be robust to noise, gene networks should include cell signaling that compartmentalizes the embryo into small regions where cell behaviors are regulated differently. If, instead, cell behaviors are equally regulated over large regions, complex but non-robust morphologies arise. We explain how compartmentalization enhances robustness and why it is a general feature of animal development. Our results are consistent with theories proposing that robustness evolved by the co-option of gene networks and extracellular cell signaling in early animal evolution.

Environmental Causation of Turtle Scute Anomalies in ovo and in silico.

The turtle shell is often described as an evolutionary novelty that facilitated the radiation of the clade Testudines. The scutes, or keratinous plates, of the turtle shell are hypothesized to be patterned by reaction-diffusion dynamics, and this property of their development provides explanatory power to mechanisms of anomalous variation. A mathematical model of scute development predicts that anomalous variation in the phylogenetically stable pattern of scutes is achieved by environmental influence on the developmental program. We test this prediction with data on patterns of scute variation from natural nests and controlled incubation of sea turtle eggs in Florida and Western Australia. We find that high temperatures are sufficient to produce anomalous patterns in turtle scutes, and that this correlation is even stronger when conditions are dry. Furthermore, we find that the patterns of variation are not random; greater anomalous variation is found in the midline vertebral scutes and during a critical period of turtle development.

Computational modeling of development by epithelia, mesenchyme and their interactions: a unified model.

MOTIVATION: The transformation of the embryo during development requires complex gene networks, cell signaling and gene-regulated cell behaviors (division, adhesion, polarization, apoptosis, contraction, extracellular matrix secretion, signal secretion and reception, etc.). There are several models of development implementing these phenomena, but none considers at the same time the very different bio-mechanical properties of epithelia, mesenchyme, extracellular matrix and their interactions. RESULTS: Here, we present a new computational model and accompanying open-source software, EmbryoMaker, that allows the user to simulate custom developmental processes by designing custom gene networks capable of regulating cell signaling and all animal basic cell behaviors. We also include an editor to implement different initial conditions, mutations and experimental manipulations. We show the applicability of the model by simulating several complex examples of animal development. AVAILABILITY AND IMPLEMENTATION: The source code can be downloaded from: http://www.biocenter.helsinki.fi/salazar/software.html. CONTACT: isalazar@mappi.helsinki.fi SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The origin and loss of periodic patterning in the turtle shell.

The origin of the turtle shell over 200 million years ago greatly modified the amniote body plan, and the morphological plasticity of the shell has promoted the adaptive radiation of turtles. The shell, comprising a dorsal carapace and a ventral plastron, is a layered structure formed by basal endochondral axial skeletal elements (ribs, vertebrae) and plates of bone, which are overlain by keratinous ectodermal scutes. Studies of turtle development have mostly focused on the bones of the shell; however, the genetic regulation of the epidermal scutes has not been investigated. Here, we show that scutes develop from an array of patterned placodes and that these placodes are absent from a soft-shelled turtle in which scutes were lost secondarily. Experimentally inhibiting Shh, Bmp or Fgf signaling results in the disruption of the placodal pattern. Finally, a computational model is used to show how two coupled reaction-diffusion systems reproduce both natural and abnormal variation in turtle scutes. Taken together, these placodal signaling centers are likely to represent developmental modules that are responsible for the evolution of scutes in turtles, and the regulation of these centers has allowed for the diversification of the turtle shell.

Transdifferentiation of pancreatic cells by loss of contact-mediated signaling.

BACKGROUND: Replacement of dysfunctional ß-cells in the islets of Langerhans by transdifferentiation of pancreatic acinar cells has been proposed as a regenerative therapy for diabetes. Adult acinar cells spontaneously revert to a multipotent state upon tissue dissociation in vitro and can be stimulated to redifferentiate into ß-cells. Despite accumulating evidence that contact-mediated signals are involved, the mechanisms regulating acinar-to-islet cell transdifferentiation remain poorly understood. RESULTS: In this study, we propose that the crosstalk between two contact-mediated signaling mechanisms, lateral inhibition and lateral stabilization, controls cell fate stability and transdifferentiation of pancreatic cells. Analysis of a mathematical model combining gene regulation with contact-mediated signaling reveals the multistability of acinar and islet cell fates. Inhibition of one or both modes of signaling results in transdifferentiation from the acinar to the islet cell fate, either by dedifferentiation to a multipotent state or by direct lineage switching. CONCLUSIONS: This study provides a theoretical framework to understand the role of contact-mediated signaling in pancreatic cell fate control that may help to improve acinar-to-islet cell transdifferentiation strategies for ß-cell neogenesis.