Metamodels and phylogenetic replication: a systematic approach to the evolution of developmental pathways

Multiple recent co-options of Optix associated with novel traits in adaptive butterfly wing radiations

Background

While the ecological factors that drive phenotypic radiations are often well understood, less is known about the generative mechanisms that cause the emergence and subsequent diversification of novel features. Heliconius butterflies display an extraordinary diversity of wing patterns due in part to mimicry and sexual selection. Identifying the genetic drivers of this crucible of evolution is now within reach, as it was recently shown that cis-regulatory variation of the optix transcription factor explains red pattern differences in the adaptive radiations of the Heliconius melpomene and Heliconius erato species groups.

Results

Here, we compare the developmental expression of the Optix protein across a large phylogenetic sample of butterflies and infer that its color patterning role originated at the base of the neotropical passion-vine butterfly clade (Lepidoptera, Nymphalidae, Tribe: Heliconiini), shortly predating multiple Optix-driven wing pattern radiations in the speciose Heliconius and Eueides genera. We also characterize novel Optix and Doublesex expression in the male-specific pheromone wing scales of the basal heliconiines Dryas and Agraulis, thus illustrating that within the Heliconinii lineage, Optix has been evolutionarily redeployed in multiple contexts in association with diverse wing features.

Conclusions

Our findings reveal that the repeated co-option of Optix into various aspects of wing scale specification was associated with multiple evolutionary novelties over a relatively short evolutionary time scale. In particular, the recruitment of Optix expression in colored scale cell precursors was a necessary condition to the explosive diversification of passion-vine butterfly wing patterns. The novel deployment of a gene followed by spatial modulation of its expression in a given cell type could be a common mode of developmental innovation for triggering phenotypic radiations.

A perspective on micro-evo-devo: progress and potential

The term "micro-evo-devo" refers to the combined study of the genetic and developmental bases of natural variation in populations and the evolutionary forces that have shaped this variation. It thus represents a synthesis of the fields of evolutionary developmental biology and population genetics. As has been pointed out by several others, this synthesis can provide insights into the evolution of organismal form and function that have not been possible within these individual disciplines separately. Despite a number of important successes in micro-evo-devo, however, it appears that evo devo and population genetics remain largely separate spheres of research, limiting their ability to address evolutionary questions. This also risks pushing contemporary evo devo to the fringes of evolutionary biology because it does not describe the causative molecular changes underlying evolution or the evolutionary forces involved. Here we reemphasize the theoretical and practical importance of micro-evo-devo as a strategy for understanding phenotypic evolution, review the key recent insights that it has provided, and present a perspective on both the potential and the remaining challenges of this exciting interdisciplinary field.

Genetic architecture of parallel pelvic reduction in ninespine sticklebacks

Teleost fish genomes are known to be evolving faster than those of other vertebrate taxa. Thus, fish are suited to address the extent to which the same vs. different genes are responsible for similar phenotypic changes in rapidly evolving genomes of evolutionary independent lineages. To gain insights into the genetic basis and evolutionary processes behind parallel phenotypic changes within and between species, we identified the genomic regions involved in pelvic reduction in Northern European ninespine sticklebacks (Pungitius pungitius) and compared them to those of North American ninespine and threespine sticklebacks (Gasterosteus aculeatus). To this end, we conducted quantitative trait locus (QTL) mapping using 283 F₂ progeny from an interpopulation cross. Phenotypic analyses indicated that pelvic reduction is a recessive trait and is inherited in a simple Mendelian fashion. Significant QTL for pelvic spine and girdle lengths were identified in the region of the Pituitary homeobox transcription factor 1 (Pitx1) gene, also responsible for pelvic reduction in threespine sticklebacks. The fact that no QTL was observed in the region identified in the mapping study of North American ninespine sticklebacks suggests that an alternative QTL for pelvic reduction has emerged in this species within the past 1.6 million years after the split between Northern European and North American populations. In general, our study provides empirical support for the view that alternative genetic mechanisms that lead to similar phenotypes can evolve over short evolutionary time scales.

Flower color as a model system for studies of plant evo-devo

Even though pigmentation traits have had substantial impacts on the field of animal evolutionary developmental biology, they have played only relatively minor roles in plant evo-devo. This is surprising given the often direct connection between flower color and fitness variation mediated through the effects of pollinators. At the same time, ecological and evolutionary genetic studies have utilized the molecular resources available for the anthocyanin pathway to generate several examples of the molecular basis of putatively adaptive transitions in flower color. Despite this opportunity to synthesize experimental approaches in ecology, evolution, and developmental biology, the investigation of many fundamental questions in evo-devo using this powerful model is only at its earliest stages. For example, a long-standing question is whether predictable genetic changes accompany the repeated evolution of a trait. Due to the conserved nature of the biochemical and regulatory control of anthocyanin biosynthesis, it has become possible to determine whether, and under what circumstances, different types of mutations responsible for flower color variation are preferentially targeted by natural selection. In addition, because plants use anthocyanin and related compounds in vegetative tissue for other important physiological functions, the identification of naturally occurring transitions from unpigmented to pigmented flowers provides the opportunity to examine the mechanisms by which regulatory networks are co-opted into new developmental domains. Here, we review what is known about the ecological and molecular basis of anthocyanic flower color transitions in natural systems, focusing on the evolutionary and developmental features involved. In doing so, we provide suggestions for future work on this trait and suggest that there is still much to be learned from the evolutionary development of flower color transitions in nature.

Combining genome-wide methods to investigate the genetic complexity of courtship song variation in Drosophila melanogaster

Little is currently known about the genetic complexity of quantitative behavioral variation, the types of genes involved, or their effects on intermediate phenotypes. Here, we conduct a genome-wide association study of Drosophila melanogaster courtship song variation using 168 sequenced inbred lines, and fail to find highly significant associations. However, by combining these data with results from a well-powered Evolve and Resequence (E&R) study on the same trait, we provide statistical evidence that some power to associate genotype and phenotype is available. Genes that are significant in both analyses are enriched for expression in the nervous system, and affect neural development and synaptic growth when perturbed. Quantitative complementation at one of these loci, Syntrophin-like 1, supports a hypothesis that variation at this locus affects variation in the inter-pulse interval of courtship song. These results suggest that experimental evolution may provide an approach for genome-scale replication in Drosophila.

The Loci of repeated evolution: a catalog of genetic hotspots of phenotypic variation

What is the nature of the genetic changes underlying phenotypic evolution? We have catalogued 1008 alleles described in the literature that cause phenotypic differences among animals, plants, and yeasts. Surprisingly, evolution of similar traits in distinct lineages often involves mutations in the same gene ("gene reuse"). This compilation yields three important qualitative implications about repeated evolution. First, the apparent evolution of similar traits by gene reuse can be traced back to two alternatives, either several independent causative mutations or a single original mutational event followed by sorting processes. Second, hotspots of evolution-defined as the repeated occurrence of de novo mutations at orthologous loci and causing similar phenotypic variation-are omnipresent in the literature with more than 100 examples covering various levels of analysis, including numerous gain-of-function events. Finally, several alleles of large effect have been shown to result from the aggregation of multiple small-effect mutations at the same hotspot locus, thus reconciling micromutationist theories of adaptation with the empirical observation of large-effect variants. Although data heterogeneity and experimental biases prevented us from extracting quantitative trends, our synthesis highlights the existence of genetic paths of least resistance leading to viable evolutionary change.

Investigating the role of the melanocortin-1 receptor gene in an extreme case of microgeographical variation in the pattern of melanin-based plumage pigmentation

The Réunion grey white-eye (Zosterops borbonicus) is a single-island endemic passerine bird that exhibits striking geographically structured melanic polymorphism at a very small spatial scale. We investigated the genetic basis of this color polymorphism by testing whether the melanocortin-1 receptor (MC1R), a gene often involved in natural melanic polymorphism in birds, was associated with the observed plumage variation. Although we found three non-synonymous mutations, we detected no association between MC1R variants and color morphs, and the main amino-acid variant found in the Réunion grey white-eye was also present at high frequency in the Mauritius grey white-eye (Zosterops mauritianus), its sister species which shows no melanic polymorphism. In addition, neutrality tests and analysis of population structure did not reveal any obvious pattern of positive or balancing selection acting on MC1R. Altogether these results indicate that MC1R does not play a role in explaining the melanic variation observed in the Réunion grey white-eye. We propose that other genes such as POMC, Agouti or any other genes involved in pigment synthesis will need to be investigated in future studies if we are to understand how selection shapes complex patterns of melanin-based plumage pigmentation.

TRIAL REGISTRATION

All sequences submitted to Genbank. Accession number: JX914505 to JX914564.

The ontogeny of color: developmental origins of divergent pigmentation in Drosophila americana and D. novamexicana

Pigmentation is a model trait for evolutionary and developmental analysis that is particularly amenable to molecular investigation in the genus Drosophila. To better understand how this phenotype evolves, we examined divergent pigmentation and gene expression over developmental time in the dark-bodied D. americana and its light-bodied sister species D. novamexicana. Prior genetic analysis implicated two enzyme-encoding genes, tan and ebony, in pigmentation divergence between these species, but questions remain about the underlying molecular mechanisms. Here, we describe stages of pupal development in both species and use this staging to determine when pigmentation develops and diverges between D. americana and D. novamexicana. For the developmental stages encompassing pigment divergence, we compare mRNA expression of tan and ebony over time and between species. Finally, we use allele-specific expression assays to determine whether interspecific differences in mRNA abundance have a cis-regulatory basis and find evidence of cis-regulatory divergence for both tan and ebony. cis-regulatory divergence affecting tan had a small effect on mRNA abundance and was limited to a few developmental stages, yet previous data suggests that this divergence is likely to be biologically meaningful. Our study suggests that small and developmentally transient expression changes may contribute to phenotypic diversification more often than commonly appreciated. Recognizing the potential phenotypic impact of such changes is important for a scientific community increasingly focused on dissecting quantitative variation, but detecting these types of changes will be a major challenge to elucidating the molecular basis of complex traits.

Unraveling the thread of nature's tapestry: the genetics of diversity and convergence in animal pigmentation

Animals display incredibly diverse color patterns yet little is known about the underlying genetic basis of these phenotypes. However, emerging results are reshaping our view of how the process of phenotypic evolution occurs. Here, we outline recent research from three particularly active areas of investigation: melanin pigmentation in Drosophila, wing patterning in butterflies, and pigment variation in lizards. For each system, we highlight (i) the function and evolution of color variation, (ii) various approaches that have been used to explore the genetic basis of pigment variation, and (iii) conclusions regarding the genetic basis of convergent evolution which have emerged from comparative analyses. Results from these studies indicate that natural variation in pigmentation is a particularly powerful tool to examine the molecular basis of evolution, especially with regard to convergent or parallel evolution. Comparison of these systems also reveals that the molecular basis of convergent evolution is heterogeneous, sometimes involving conserved mechanisms and sometimes not. In the near future, additional work in other emerging systems will substantially expand the scope of available comparisons.

Lessons from flower colour evolution on targets of selection

The genetic basis of flower colour evolution provides a useful system to address the debate over the relative contribution of regulatory vs. functional mutations in evolution. The relative importance of these two categories depends on the type of flower colour transition and the genes involved in those transitions. These differences reflect differences in the degree of deleterious pleiotropy associated with functional inactivation of various anthocyanin pathway genes. Our findings illustrate how generalized statements regarding the contributions of regulatory and functional mutations to broad categories of traits, such as morphological vs. physiological, ignore differences among traits within categories and in doing so overlook important factors determining the relative importance of regulatory and functional mutations.

Diversification of complex butterfly wing patterns by repeated regulatory evolution of a Wnt ligand

Although animals display a rich variety of shapes and patterns, the genetic changes that explain how complex forms arise are still unclear. Here we take advantage of the extensive diversity of Heliconius butterflies to identify a gene that causes adaptive variation of black wing patterns within and between species. Linkage mapping in two species groups, gene-expression analysis in seven species, and pharmacological treatments all indicate that cis-regulatory evolution of the WntA ligand underpins discrete changes in color pattern features across the Heliconius genus. These results illustrate how the direct modulation of morphogen sources can generate a wide array of unique morphologies, thus providing a link between natural genetic variation, pattern formation, and adaptation.

Repeated origin and loss of adhesive toepads in geckos

Geckos are well known for their extraordinary clinging abilities and many species easily scale vertical or even inverted surfaces. This ability is enabled by a complex digital adhesive mechanism (adhesive toepads) that employs van der Waals based adhesion, augmented by frictional forces. Numerous morphological traits and behaviors have evolved to facilitate deployment of the adhesive mechanism, maximize adhesive force and enable release from the substrate. The complex digital morphologies that result allow geckos to interact with their environment in a novel fashion quite differently from most other lizards. Details of toepad morphology suggest multiple gains and losses of the adhesive mechanism, but lack of a comprehensive phylogeny has hindered efforts to determine how frequently adhesive toepads have been gained and lost. Here we present a multigene phylogeny of geckos, including 107 of 118 recognized genera, and determine that adhesive toepads have been gained and lost multiple times, and remarkably, with approximately equal frequency. The most likely hypothesis suggests that adhesive toepads evolved 11 times and were lost nine times. The overall external morphology of the toepad is strikingly similar in many lineages in which it is independently derived, but lineage-specific differences are evident, particularly regarding internal anatomy, with unique morphological patterns defining each independent derivation.

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.

L'enfant terrible at 30: the maturation of evolutionary developmental biology

The recent Keystone Symposium on Evolutionary Developmental Biology at Tahoe City in February 2011 provided an opportunity to take stock of where the past three decades have brought this interdisciplinary field. It revealed maturation on several fronts, including increased experimental rigor, the softening of dichotomies that were crucial to its founding and growth, and its growing relevance to both basic and biomedical biology.

Involvement of the conserved Hox gene Antennapedia in the development and evolution of a novel trait

BACKGROUND

Hox proteins specify segment identity during embryogenesis and have typical associated expression patterns. Changes in embryonic expression and activity of Hox genes were crucial in the evolution of animal body plans, but their role in the post-embryonic development of lineage-specific traits remains largely unexplored. Here, we focus on the insect Hox genes Ultrabithorax (Ubx) and Antennapedia (Antp), and implicate the latter in the formation and diversification of novel, butterfly-specific wing patterns.

RESULTS

First, we describe a conserved pattern of Ubx expression and a novel pattern of Antp expression in wing discs of Bicyclus anynana butterflies. The discrete, reiterated domains of Antp contrast with the typical expression of Hox genes in single continuous regions in arthropod embryos. Second, we show that this pattern is associated with the establishment of the organizing centres of eyespots. Antp upregulation is the earliest event in organizer development described to date, and in contrast to all genes implicated in eyespot formation, is exclusive to those centres. Third, our comparative analysis of gene expression across nymphalids reveals unexpected differences in organizer determination.

CONCLUSIONS

We show that the Antp's recruitment for the formation of novel traits in butterfly wing discs involved the evolution of new expression domains, and is restricted to a particular lineage. This study contributes novel insights into the evolution of Antp expression, as well as into the genetic mechanisms underlying morphological diversification. Our results also underscore how a wider representation of morphological and phylogenetic diversity is essential in evolutionary developmental biology.

Convergence, adaptation, and constraint

Convergent evolution of similar phenotypic features in similar environmental contexts has long been taken as evidence of adaptation. Nonetheless, recent conceptual and empirical developments in many fields have led to a proliferation of ideas about the relationship between convergence and adaptation. Despite criticism from some systematically minded biologists, I reaffirm that convergence in taxa occupying similar selective environments often is the result of natural selection. However, convergent evolution of a trait in a particular environment can occur for reasons other than selection on that trait in that environment, and species can respond to similar selective pressures by evolving nonconvergent adaptations. For these reasons, studies of convergence should be coupled with other methods-such as direct measurements of selection or investigations of the functional correlates of trait evolution-to test hypotheses of adaptation. The independent acquisition of similar phenotypes by the same genetic or developmental pathway has been suggested as evidence of constraints on adaptation, a view widely repeated as genomic studies have documented phenotypic convergence resulting from change in the same genes, sometimes even by the same mutation. Contrary to some claims, convergence by changes in the same genes is not necessarily evidence of constraint, but rather suggests hypotheses that can test the relative roles of constraint and selection in directing phenotypic evolution.

The developmental role of Agouti in color pattern evolution

Animal color patterns can affect fitness in the wild; however, little is known about the mechanisms that control their formation and subsequent evolution. We took advantage of two locally camouflaged populations of Peromyscus mice to show that the negative regulator of adult pigmentation, Agouti, also plays a key developmental role in color pattern evolution. Genetic and functional analyses showed that ventral-specific embryonic expression of Agouti establishes a prepattern by delaying the terminal differentiation of ventral melanocytes. Moreover, a skin-specific increase in both the level and spatial domain of Agouti expression prevents melanocyte maturation in a regionalized manner, resulting in a novel and adaptive color pattern. Thus, natural selection favors late-acting, tissue-specific changes in embryonic Agouti expression to produce large changes in adult color pattern.

Population genetics, pleiotropy, and the preferential fixation of mutations during adaptive evolution

Ongoing debate centers on whether certain types of mutations are fixed preferentially during adaptive evolution. Although there has been much discussion, no quantitative framework currently exists to test for these biases. Here, we describe a method for distinguishing between the two processes that likely account for biased rates of substitution: variation in mutation rates and variation in the probability that a mutation becomes fixed once it arises. We then use this approach to examine the type and magnitude of these biases during evolutionary transitions across multiple scales: those involving repeated origins of individual traits (flower color change), and transitions involving broad suites of traits (morphological and physiological trait evolution in plants and animals). We show that fixation biases can be strong at both levels of comparison, likely due to differences in the magnitude of deleterious pleiotropy associated with alternative mutation categories. However, we also show that the scale at which these comparisons are made greatly influences the results, as broad comparisons that simultaneously analyze multiple traits obscure heterogeneity in the direction and magnitude of these biases. We conclude that preferential fixation of mutations likely is common in nature, but should be studied on a trait-by-trait basis.

Convergence in pigmentation at multiple levels: mutations, genes and function

Convergence--the independent evolution of the same trait by two or more taxa--has long been of interest to evolutionary biologists, but only recently has the molecular basis of phenotypic convergence been identified. Here, we highlight studies of rapid evolution of cryptic coloration in vertebrates to demonstrate that phenotypic convergence can occur at multiple levels: mutations, genes and gene function. We first show that different genes can be responsible for convergent phenotypes even among closely related populations, for example, in the pale beach mice inhabiting Florida's Gulf and Atlantic coasts. By contrast, the exact same mutation can create similar phenotypes in distantly related species such as mice and mammoths. Next, we show that different mutations in the same gene need not be functionally equivalent to produce similar phenotypes. For example, separate mutations produce divergent protein function but convergent pale coloration in two lizard species. Similarly, mutations that alter the expression of a gene in different ways can, nevertheless, result in similar phenotypes, as demonstrated by sister species of deer mice. Together these studies underscore the importance of identifying not only the genes, but also the precise mutations and their effects on protein function, that contribute to adaptation and highlight how convergence can occur at different genetic levels.

Transcriptional network structure has little effect on the rate of regulatory evolution in yeast

Studies in evolutionary developmental biology suggest that the structure of genetic pathways may bias the fixation of natural variation toward particular nodes in these pathways. In an attempt to test this trend genome wide, we integrated several previously published data sets to examine whether the position of genes in the whole-genome transcriptional network of Saccharomyces cerevisiae is associated with the amount of cis-regulatory expression divergence between S. cerevisiae and its sibling species Saccharomyces paradoxus. We find little evidence for an association between connectivity and divergence in the global network that combines data from multiple conditions. However, relationships between connectivity and divergence are apparent in some of the smaller subnetworks. Despite a slight tendency for genes with more transcriptional interactions to show greater divergence, these differences explain no more than a small fraction of variation in evolutionary rates. These results suggest that the systems biology focus on large interactomes may miss some critical details of local interactions. More detailed experimental analysis will be needed to define the genetic pathways that control specific phenotypic traits and quantify the rate of regulatory changes at different points in these pathways.

Reciprocal insights into adaptation from agricultural and evolutionary studies in tomato

Although traditionally separated by different aims and methodologies, research on agricultural and evolutionary problems shares a common goal of understanding the mechanisms underlying functionally important traits. As such, research in both fields offers potential complementary and reciprocal insights. Here, we discuss adaptive stress responses (specifically to water stress) as an example of potentially fruitful research reciprocity, where agricultural research has clearly produced advances that could benefit evolutionary studies, while evolutionary studies offer approaches and insights underexplored in crop studies. We focus on research on Solanum species that include the domesticated tomato and its wild relatives. Integrated approaches to understanding ecological adaptation are particularly attractive in tomato and its wild relatives: many presumptively adaptive phenotypic differences characterize wild species, and the physiological and mechanistic basis of many relevant traits and environmental responses has already been examined in the context of cultivated tomato and some wild species. We highlight four specific instances where these reciprocal insights can be combined to better address questions that are fundamental both to agriculture and evolution.

Parallel adaptation: one or many waves of advance of an advantageous allele?

Models for detecting the effect of adaptation on population genomic diversity are often predicated on a single newly arisen mutation sweeping rapidly to fixation. However, a population can also adapt to a new environment by multiple mutations of similar phenotypic effect that arise in parallel, at the same locus or different loci. These mutations can each quickly reach intermediate frequency, preventing any single one from rapidly sweeping to fixation globally, leading to a "soft" sweep in the population. Here we study various models of parallel mutation in a continuous, geographically spread population adapting to a global selection pressure. The slow geographic spread of a selected allele due to limited dispersal can allow other selected alleles to arise and start to spread elsewhere in the species range. When these different selected alleles meet, their spread can slow dramatically and so initially form a geographic patchwork, a random tessellation, which could be mistaken for a signal of local adaptation. This spatial tessellation will dissipate over time due to mixing by migration, leaving a set of partial sweeps within the global population. We show that the spatial tessellation initially formed by mutational types is closely connected to Poisson process models of crystallization, which we extend. We find that the probability of parallel mutation and the spatial scale on which parallel mutation occurs are captured by a single compound parameter, a characteristic length, which reflects the expected distance a spreading allele travels before it encounters a different spreading allele. This characteristic length depends on the mutation rate, the dispersal parameter, the effective local density of individuals, and to a much lesser extent the strength of selection. While our knowledge of these parameters is poor, we argue that even in widely dispersing species, such parallel geographic sweeps may be surprisingly common. Thus, we predict that as more data become available, many more examples of intraspecies parallel adaptation will be uncovered.