The dynamics of developmental system drift in the gene network underlying wing polyphenism in ants: a mathematical model

Patterns of selection across gene regulatory networks

Gene regulatory networks (GRNs) are the core engine of organismal development. If we would like to understand the origin and diversification of phenotypes, it is necessary to consider the structure of GRNs in order to reconstruct the links between genetic mutations and phenotypic change. Much of the progress in evolutionary developmental biology, however, has occurred without a nuanced consideration of the evolution of functional relationships between genes, especially in the context of their broader network interactions. Characterizing and comparing GRNs across traits and species in a more detailed way will allow us to determine how network position influences what genes drive adaptive evolution. In this perspective paper, we consider the architecture of developmental GRNs and how positive selection strength may vary across a GRN. We then propose several testable models for these patterns of selection and experimental approaches to test these models.

My road to the ants: A model clade for eco-evo-devo

This chapter is the story of how I pioneered ants as a system for studying eco-evo-devo, a field that integrates developmental biology with ecology and evolutionary biology. One aim of eco-evo-devo is to understand how the interactions between genes and their environments during development facilitates the origin and evolution of novel phenotypes. In a series of six parts, I review some of the key discoveries from my lab on how novel worker caste systems in ants--soldiers and supersoldiers--originated and evolved. I also discuss some of the ideas that emerged from these discoveries, including the role that polyphenisms, hidden developmental potentials, and rudimentary organs play in facilitating developmental and evolutionary change. As superorganisms, I argue that ants are uniquely positioned to reveal types of variation that are often difficult to observe in nature. In doing so, they have the potential to transform our view of biology and provide new perspectives in medicine, agriculture, and biodiversity conservation. With my story I hope to inspire the next generation of biologists to continue exploring the unknown regions of phenotypic space to solve some of our most pressing societal challenges.

The origin of wing polyphenism in ants: An eco-evo-devo perspective

The evolution of eusociality, where solitary individuals integrate into a single colony, is a major transition in individuality. In ants, the origin of eusociality coincided with the origin of a wing polyphenism approximately 160 million years ago, giving rise to colonies with winged queens and wingless workers. As a consequence, both eusociality and wing polyphenism are nearly universal features of all ants. Here, we synthesize fossil, ecological, developmental, and evolutionary data in an attempt to understand the factors that contributed to the origin of wing polyphenism in ants. We propose multiple models and hypotheses to explain how wing polyphenism is orchestrated at multiple levels, from environmental cues to gene networks. Furthermore, we argue that the origin of wing polyphenism enabled the subsequent evolution of morphological diversity across the ants. We finally conclude by outlining several outstanding questions for future work.

System-level analyses of keystone genes required for mammalian tooth development

When a null mutation of a gene causes a complete developmental arrest, the gene is typically considered essential for life. Yet, in most cases, null mutations have more subtle effects on the phenotype. Here we used the phenotypic severity of mutations as a tool to examine system‐level dynamics of gene expression. We classify genes required for the normal development of the mouse molar into different categories that range from essential to subtle modification of the phenotype. Collectively, we call these the developmental keystone genes. Transcriptome profiling using microarray and RNAseq analyses of patterning stage mouse molars show highly elevated expression levels for genes essential for the progression of tooth development, a result reminiscent of essential genes in single‐cell organisms. Elevated expression levels of progression genes were also detected in developing rat molars, suggesting evolutionary conservation of this system‐level dynamics. Single‐cell RNAseq analyses of developing mouse molars reveal that even though the size of the expression domain, measured in the number of cells, is the main driver of organ‐level expression, progression genes show high cell‐level transcript abundances. Progression genes are also upregulated within their pathways, which themselves are highly expressed. In contrast, a high proportion of the genes required for normal tooth patterning are secreted ligands that are expressed in fewer cells than their receptors and intracellular components. Overall, even though expression patterns of individual genes can be highly different, conserved system‐level principles of gene expression can be detected using phenotypically defined gene categories.

The phenotypic severity of mutations on mouse teeth is used to classify genes.

Genes essential for the progression of odontogenesis are highly expressed at the organ and cell level.

Many of the genes required for normal patterning are locally expressed ligands.

Interruption points in the wing gene regulatory network underlying wing polyphenism evolved independently in male and female morphs in Cardiocondyla ants

Wing polyphenism in ants, which produces a winged female queen caste and a wingless female worker caste, evolved approximately 150 million years ago and has been key to the remarkable success of ants. Approximately 20 million years ago, the myrmicine ant genus Cardiocondyla evolved an additional wing polyphenism among males producing two male morphs: wingless males that fight to enhance mating success and winged males that disperse. Here we show that interruption of rudimentary wing-disc development in larvae of the ant Cardiocondyla obscurior occurs further downstream in the network in wingless males as compared with wingless female workers. This pattern is corroborated in C. kagutsuchi, a species from a different clade within the genus, indicating that late interruption of wing development in males is conserved across Cardiocondyla. Therefore, our results show that the novel male wing polyphenism was not developmentally constrained by the pre-existing female wing polyphenism and evolved through independent alteration of interruption points in the wing gene network. Furthermore, a comparison of adult morphological characters in C. obscurior reveals that developmental trajectories lead to similar morphological trait integration between winged and wingless females, but dramatically different integration between winged and wingless males. This suggests that the alternative sex-specific developmental routes to achieve winglessness in the genus Cardiocondyla may have evolved through different selection regimes acting on wingless males and females.

Genes and genomic processes underpinning the social lives of ants

The >15000 ant species are all highly social and show great variation in colony organization, complexity and behavior. The mechanisms by which such sociality evolved, as well as those underpinning the elaboration of ant societies since their ∼140 million year old common ancestor, have long been pondered. Here, we review recent insights generated using various genomic approaches. This includes understanding the molecular mechanisms underlying caste differentiation and the diversity of social structures, studying the impact of eusociality on genomic evolutionary rates, and investigating gene expression changes associated with differences in lifespan between castes. Furthermore, functional studies involving RNAi and CRISPR have recently been successfully applied to ants, opening the door to exciting research that promises to revolutionize the understanding of the evolution and diversification of social living.

An Evo-Devo Perspective on Multicellular Development of Myxobacteria

The transition to multicellularity, recognized as one the major transitions in evolution, has occurred independently several times. While multicellular development has been extensively studied in zygotic organisms including plant and animal groups, just a few aggregative multicellular organisms have been employed as model organisms for the study of multicellularity. Studying different evolutionary origins and modes of multicellularity enables comparative analyses that can help identifying lineage-specific aspects of multicellular evolution and generic factors and mechanisms involved in the transition to multicellularity. Among aggregative multicellular organisms, myxobacteria are a valuable system to explore the particularities that aggregation confers to the evolution of multicellularity and mechanisms shared with clonal organisms. Moreover, myxobacteria species develop fruiting bodies displaying a range of morphological diversity. In this review, we aim to synthesize diverse lines of evidence regarding myxobacteria development and discuss them in the context of Evo-Devo concepts and approaches. First, we briefly describe the developmental processes in myxobacteria, present an updated comparative analysis of the genes involved in their developmental processes and discuss these and other lines of evidence in terms of co-option and developmental system drift, two concepts key to Evo-Devo studies. Next, as has been suggested from Evo-Devo approaches, we discuss how broad comparative studies and integration of diverse genetic, physicochemical, and environmental factors into experimental and theoretical models can further our understanding of myxobacterial development, phenotypic variation, and evolution.

Perspectives on Gene Regulatory Network Evolution

Animal development proceeds through the activity of genes and their cis-regulatory modules (CRMs) working together in sets of gene regulatory networks (GRNs). The emergence of species-specific traits and novel structures results from evolutionary changes in GRNs. Recent work in a wide variety of animal models, and particularly in insects, has started to reveal the modes and mechanisms of GRN evolution. I discuss here various aspects of GRN evolution and argue that developmental system drift (DSD), in which conserved phenotype is nevertheless a result of changed genetic interactions, should regularly be viewed from the perspective of GRN evolution. Advances in methods to discover related CRMs in diverse insect species, a critical requirement for detailed GRN characterization, are also described.

The evolutionary dynamics of major regulators for sexual development among Hymenoptera species

All hymenopteran species, such as bees, wasps and ants, are characterized by the common principle of haplodiploid sex determination in which haploid males arise from unfertilized eggs and females from fertilized eggs. The underlying molecular mechanism has been studied in detail in the western honey bee Apis mellifera, in which the gene complementary sex determiner (csd) acts as primary signal of the sex determining pathway, initiating female development by csd-heterozygotes. Csd arose from gene duplication of the feminizer (fem) gene, a transformer (tra) ortholog, and mediates in conjunction with transformer2 (tra2) sex-specific splicing of fem. Comparative molecular analyses identified fem/tra and its downstream target doublesex (dsx) as conserved unit within the sex determining pathway of holometabolous insects. In this study, we aim to examine evolutionary differences among these key regulators. Our main hypothesis is that sex determining key regulators in Hymenoptera species show signs of coevolution within single phylogenetic lineages. We take advantage of several newly sequenced genomes of bee species to test this hypothesis using bioinformatic approaches. We found evidences that duplications of fem are restricted to certain bee lineages and notable amino acid differences of tra2 between Apis and non-Apis species propose structural changes in Tra2 protein affecting co-regulatory function on target genes. These findings may help to gain deeper insights into the ancestral mode of hymenopteran sex determination and support the common view of the remarkable evolutionary flexibility in this regulatory pathway.

The Comet Cometh: Evolving Developmental Systems

In a recent opinion piece, Denis Duboule has claimed that the increasing shift towards systems biology is driving evolutionary and developmental biology apart, and that a true reunification of these two disciplines within the framework of evolutionary developmental biology (EvoDevo) may easily take another 100 years. He identifies methodological, epistemological, and social differences as causes for this supposed separation. Our article provides a contrasting view. We argue that Duboule’s prediction is based on a one-sided understanding of systems biology as a science that is only interested in functional, not evolutionary, aspects of biological processes. Instead, we propose a research program for an evolutionary systems biology, which is based on local exploration of the configuration space in evolving developmental systems. We call this approach—which is based on reverse engineering, simulation, and mathematical analysis—the natural history of configuration space. We discuss a number of illustrative examples that demonstrate the past success of local exploration, as opposed to global mapping, in different biological contexts. We argue that this pragmatic mode of inquiry can be extended and applied to the mathematical analysis of the developmental repertoire and evolutionary potential of evolving developmental mechanisms and that evolutionary systems biology so conceived provides a pragmatic epistemological framework for the EvoDevo synthesis.

Gene duplication and the environmental regulation of physiology and development

When different life stages have different environmental tolerances, development needs to be regulated so that each life stage experiences environmental conditions that are suitable for it, if fitness is to be maintained. Restricting the timing of developmental transitions to occur under specific combinations of environmental conditions is therefore adaptively important. However, impeding development can itself incur demographic and fitness costs. How do organisms regulate development and physiological processes so that they occur under the broadest range of permissive conditions? Gene duplication offers one solution: Multiple genes contribute to the same downstream process, but do so under distinct combinations of environmental conditions. We present a simple model to examine how environmental sensitivities of genes and how gene duplication influence the distribution of environmental conditions under which an end process will proceed. The model shows that the duplication of genes that retain their downstream function but diverge in environmental sensitivities can allow an end process to proceed under more than one distinct combination of environmental conditions. The outcomes depend on how upstream genes regulate downstream components, which genes in the pathway have diversified in their sensitivities, and the structure of the pathway.

Cryptic genetic variation: evolution's hidden substrate

Cryptic genetic variation (CGV) is invisible under normal conditions, but it can fuel evolution when circumstances change. In theory, CGV can represent a massive cache of adaptive potential or a pool of deleterious alleles that are in need of constant suppression. CGV emerges from both neutral and selective processes, and it may inform about how human populations respond to change. CGV facilitates adaptation in experimental settings, but does it have an important role in the real world? Here, we review the empirical support for widespread CGV in natural populations, including its potential role in emerging human diseases and the growing evidence of its contribution to evolution.

The wing-patterning network in the wingless castes of Myrmicine and Formicine ant species is a mix of evolutionarily labile and non-labile genes

Wing polyphenism in ants is the ability of a single genome to produce winged or wingless castes in a colony in response to environmental cues. Although wing polyphenism is a universal and homologous feature of ants, the gene network underlying wing polyphenism is conserved in the winged castes, but is labile in the wingless castes, that is, the network is interrupted at different points in the wingless castes of different ant species. Because the expression of all genes sampled so far in this network in the wingless castes is evolutionarily labile across species, an important question is whether all "interruption points" in the network are evolutionarily labile or are there interruption points that are evolutionarily non-labile. Here we show that in the wingless castes, the expression of the gene brinker (brk), which mediates growth, patterning, and apoptosis in the Drosophila wing disc, is non-labile; it is absent in vestigial wing discs of four ants species. In contrast, the expression of engrailed (en), a gene upstream of brk is labile; it is present in some species but absent in others. In the winged castes, both brk and en expression are conserved relative to their expression in Drosophila wing discs. The differential lability of genes in the network in wingless castes may be a general feature of networks underlying polyphenic traits. This raises the possibility that some genes, like brk, may be under stabilizing selection while most others, like en, may be evolving via directional selection or neutral drift.

The inheritance of process: a dynamical systems approach

A central unresolved problem of evolutionary biology concerns the way in which evolution at the genotypic level relates to the evolution of phenotypes. This genotype-phenotype map involves developmental and physiological processes, which are complex and not well understood. These processes co-determine the rate and direction of adaptive change by shaping the distribution of phenotypic variability on which selection can act. In this study, we argue-expanding on earlier ideas by Goodwin, Oster, and Alberch-that an explicit treatment of this map in terms of dynamical systems theory can provide an integrated understanding of evolution and development. We describe a conceptual framework, which demonstrates how development determines the probability of possible phenotypic transitions-and hence the evolvability of a biological system. We use a simple conceptual model to illustrate how the regulatory dynamics of the genotype-phenotype map can be passed on from generation to generation, and how heredity itself can be treated as a dynamic process. Our model yields explanations for punctuated evolutionary dynamics, the difference between micro- and macroevolution, and for the role of the environment in major phenotypic transitions. We propose a quantitative research program in evolutionary developmental systems biology-combining experimental methods with mathematical modeling-which aims at elaborating our conceptual framework by applying it to a wide range of evolving developmental systems. This requires a large and sustained effort, which we believe is justified by the significant potential benefits of an extended evolutionary theory that uses dynamic molecular genetic data to reintegrate development and evolution.

Phenotypic models of evolution and development: geometry as destiny

Quantitative models of development that consider all relevant genes typically are difficult to fit to embryonic data alone and have many redundant parameters. Computational evolution supplies models of phenotype with relatively few variables and parameters that allows the patterning dynamics to be reduced to a geometrical picture for how the state of a cell moves. The clock and wavefront model, that defines the phenotype of somitogenesis, can be represented as a sequence of two discrete dynamical transitions (bifurcations). The expression-time to space map for Hox genes and the posterior dominance rule are phenotypes that naturally follow from computational evolution without considering the genetics of Hox regulation.

Developmental evolution in social insects: regulatory networks from genes to societies

The evolution and development of complex phenotypes in social insect colonies, such as queen-worker dimorphism or division of labor, can, in our opinion, only be fully understood within an expanded mechanistic framework of Developmental Evolution. Conversely, social insects offer a fertile research area in which fundamental questions of Developmental Evolution can be addressed empirically. We review the concept of gene regulatory networks (GRNs) that aims to fully describe the battery of interacting genomic modules that are differentially expressed during the development of individual organisms. We discuss how distinct types of network models have been used to study different levels of biological organization in social insects, from GRNs to social networks. We propose that these hierarchical networks spanning different organizational levels from genes to societies should be integrated and incorporated into full GRN models to elucidate the evolutionary and developmental mechanisms underlying social insect phenotypes. Finally, we discuss prospects and approaches to achieve such an integration.

Evolutionary history of the recruitment of conserved developmental genes in association to the formation and diversification of a novel trait

Background

The origin and modification of novel traits are important aspects of biological diversification. Studies combining concepts and approaches of developmental genetics and evolutionary biology have uncovered many examples of the recruitment, or co-option, of genes conserved across lineages for the formation of novel, lineage-restricted traits. However, little is known about the evolutionary history of the recruitment of those genes, and of the relationship between them -for example, whether the co-option involves whole or parts of existing networks, or whether it occurs by redeployment of individual genes with de novo rewiring. We use a model novel trait, color pattern elements on butterfly wings called eyespots, to explore these questions. Eyespots have greatly diversified under natural and sexual selection, and their formation involves genetic circuitries shared across insects.

Results

We investigated the evolutionary history of the recruitment and co-recruitment of four conserved transcription regulators to the larval wing disc region where circular pattern elements develop. The co-localization of Antennapedia, Notch, Distal-less, and Spalt with presumptive (eye)spot organizers was examined in 13 butterfly species, providing the largest comparative dataset available for the system. We found variation between families, between subfamilies, and between tribes. Phylogenetic reconstructions by parsimony and maximum likelihood methods revealed an unambiguous evolutionary history only for Antennapedia, with a resolved single origin of eyespot-associated expression, and many homoplastic events for Notch, Distal-less, and Spalt. The flexibility in the (co-)recruitment of the targeted genes includes cases where different gene combinations are associated with morphologically similar eyespots, as well as cases where identical protein combinations are associated with very different phenotypes.

Conclusions

The evolutionary history of gene (co-)recruitment is consistent with both divergence from a recruited putative ancestral network, and with independent co-option of individual genes. The diversity in the combinations of genes expressed in association with eyespot formation does not parallel diversity in characteristics of the adult phenotype. We discuss these results in the context of inferring homology. Our study underscores the importance of widening the representation of phylogenetic, morphological, and genetic diversity in order to establish general principles about the mechanisms behind the evolution of novel traits.

Life's attractors : understanding developmental systems through reverse engineering and in silico evolution

We propose an approach to evolutionary systems biology which is based on reverse engineering of gene regulatory networks and in silico evolutionary simulations. We infer regulatory parameters for gene networks by fitting computational models to quantitative expression data. This allows us to characterize the regulatory structure and dynamical repertoire of evolving gene regulatory networks with a reasonable amount of experimental and computational effort. We use the resulting network models to identify those regulatory interactions that are conserved, and those that have diverged between different species. Moreover, we use the models obtained by data fitting as starting points for simulations of evolutionary transitions between species. These simulations enable us to investigate whether such transitions are random, or whether they show stereotypical series of regulatory changes which depend on the structure and dynamical repertoire of an evolving network. Finally, we present a case study-the gap gene network in dipterans (flies, midges, and mosquitoes)-to illustrate the practical application of the proposed methodology, and to highlight the kind of biological insights that can be gained by this approach.

Functional data analysis for identifying nonlinear models of gene regulatory networks

Background

A key problem in systems biology is estimating dynamical models of gene regulatory networks. Traditionally, this has been done using regression or other ad-hoc methods when the model is linear. More detailed, realistic modeling studies usually employ nonlinear dynamical models, which lead to computationally difficult parameter estimation problems. Functional data analysis methods, however, offer a means to simplify fitting by transforming the problem from one of matching modeled and observed dynamics to one of matching modeled and observed time derivatives–a regression problem, albeit a nonlinear one.

Results

We formulate a functional data analysis approach for estimating the parameters of nonlinear dynamical models and evaluate this approach on data from two real systems, the gap gene system of Drosophila melanogaster and the synthetic IRMA network, which was created expressly as a test case for genetic network inference. We also evaluate the approach on simulated data sets generated by the GeneNetWeaver program, the basis for the annual DREAM reverse engineering challenge. We assess the accuracy with which the correct regulatory relationships within the networks are extracted, and consider alternative methods of regularization for the purpose of overfitting avoidance. We also show that the computational efficiency of the functional data analysis approach, and the decomposability of the resulting regression problem, allow us to explicitly enumerate and evaluate all possible regulator combinations for every gene. This gives deeper insight into the the relevance of different regulators or regulator combinations, and lets one check for alternative regulatory explanations.

Conclusions

Functional data analysis is a powerful approach for estimating detailed nonlinear models of gene expression dynamics, allowing efficient and accurate estimation of regulatory architecture.

Mechanisms driving neural crest induction and migration in the zebrafish and Xenopus laevis

The neural crest is an evolutionary adaptation, with roots in the formation of mesoderm. Modification of neural crest behavior has been is critical for the evolutionary diversification of the vertebrates and defects in neural crest underlie a range of human birth defects. There has been a tremendous increase in our knowledge of the molecular, cellular, and inductive interactions that converge on defining the neural crest and determining its behavior. While there is a temptation to look for simple models to explain neural crest behavior, the reality is that the system is complex in its circuitry. In this review, our goal is to identify the broad features of neural crest origins (developmentally) and migration (cellularly) using data from the zebrafish (teleost) and Xenopus laevis (tetrapod amphibian) in order to illuminate where general mechanisms appear to be in play, and equally importantly, where disparities in experimental results suggest areas of profitable study.

Robustness of positional specification by the Hedgehog morphogen gradient

Spatial gradients of Hedgehog signalling play a central role in many patterning events during animal development, regulating cell fate determination and tissue growth in a variety of tissues and developmental stages. Experimental evidence suggests that many of the proteins responsible for regulating Hedgehog signalling and transport are themselves targets of Hedgehog signalling, leading to multiple levels of feedback within the system. We use mathematical modelling to analyse how these overlapping feedbacks combine to regulate patterning and potentially enhance robustness in the Drosophila wing imaginal disc. Our results predict that the regulation of Hedgehog transport and stability by glypicans, as well as multiple overlapping feedbacks in the Hedgehog response network, can combine to enhance the robustness of positional specification against variability in Hedgehog levels. We also discuss potential trade-offs between robustness and additional features of the Hedgehog gradient, such as signalling range and size regulation.