optix drives the repeated convergent evolution of butterfly wing pattern mimicry

The morphological state space revisited: what do phylogenetic patterns in homoplasy tell us about the number of possible character states?

Biological variety and major evolutionary transitions suggest that the space of possible morphologies may have varied among lineages and through time. However, most models of phylogenetic character evolution assume that the potential state space is finite. Here, I explore what the morphological state space might be like, by analysing trends in homoplasy (repeated derivation of the same character state). Analyses of ten published character matrices are compared against computer simulations with different state space models: infinite states, finite states, ordered states and an 'inertial' model, simulating phylogenetic constraints. Of these, only the infinite states model results in evolution without homoplasy, a prediction which is not generally met by real phylogenies. Many authors have interpreted the ubiquity of homoplasy as evidence that the number of evolutionary alternatives is finite. However, homoplasy is also predicted by phylogenetic constraints on the morphological distance that can be traversed between ancestor and descendent. Phylogenetic rarefaction (sub-sampling) shows that finite and inertial state spaces do produce contrasting trends in the distribution of homoplasy. Two clades show trends characteristic of phylogenetic inertia, with decreasing homoplasy (increasing consistency index) as we sub-sample more distantly related taxa. One clade shows increasing homoplasy, suggesting exhaustion of finite states. Different clades may, therefore, show different patterns of character evolution. However, when parsimony uninformative characters are excluded (which may occur without documentation in cladistic studies), it may no longer be possible to distinguish inertial and finite state spaces. Interestingly, inertial models predict that homoplasy should be clustered among comparatively close relatives (parallel evolution), whereas finite state models do not. If morphological evolution is often inertial in nature, then homoplasy (false homology) may primarily occur between close relatives, perhaps being replaced by functional analogy at higher taxonomic scales.

Different gene expressions between cattle and yak provide insights into high-altitude adaptation

DNA sequence variation has been widely reported as the genetic basis for adaptation, in both humans and other animals, to the hypoxic environment experienced at high altitudes. However, little is known about the patterns of gene expression underlying such hypoxic adaptations. In this study, we examined the differences in the transcriptomes of four organs (heart, kidney, liver and lung) between yak and cattle, a pair of closely related species distributed at high and low altitudes respectively. Of the four organs examined, heart shows the greatest differentiation between the two species in terms of gene expression profiles. Detailed analyses demonstrated that some genes associated with the oxygen supply system and the defense systems that respond to threats of hypoxia are differentially expressed. In addition, genes with significantly differentiated patterns of expression in all organs exhibited an unexpected uniformity of regulation along with an elevated frequency of nonsynonymous substitutions. This co-evolution of protein sequences and gene expression patterns is likely to be correlated with the optimization of the yak metabolic system to resist hypoxia.

Population Genomics for Understanding Adaptation in Wild Plant Species

Darwin's theory of evolution by natural selection is the foundation of modern biology. However, it has proven remarkably difficult to demonstrate at the genetic, genomic, and population level exactly how wild species adapt to their natural environments. We discuss how one can use large sets of multiple genome sequences from wild populations to understand adaptation, with an emphasis on the small herbaceous plant Arabidopsis thaliana. We present motivation for such studies; summarize progress in describing whole-genome, species-wide sequence variation; and then discuss what insights have emerged from these resources, either based on sequence information alone or in combination with phenotypic data. We conclude with thoughts on opportunities with other plant species and the impact of expected progress in sequencing technology and genome engineering for studying adaptation in nature.

Divergence with gene flow across a speciation continuum of Heliconius butterflies

Background

A key to understanding the origins of species is determining the evolutionary processes that drive the patterns of genomic divergence during speciation. New genomic technologies enable the study of high-resolution genomic patterns of divergence across natural speciation continua, where taxa pairs with different levels of reproductive isolation can be used as proxies for different stages of speciation. Empirical studies of these speciation continua can provide valuable insights into how genomes diverge during speciation.

Methods

We examine variation across a handful of genomic regions in parapatric and allopatric populations of Heliconius butterflies with varying levels of reproductive isolation. Genome sequences were mapped to 2.2-Mb of the H. erato genome, including 1-Mb across the red color pattern locus and multiple regions unlinked to color pattern variation.

Results

Phylogenetic analyses reveal a speciation continuum of pairs of hybridizing races and incipient species in the Heliconius erato clade. Comparisons of hybridizing pairs of divergently colored races and incipient species reveal that genomic divergence increases with ecological and reproductive isolation, not only across the locus responsible for adaptive variation in red wing coloration, but also at genomic regions unlinked to color pattern.

Discussion

We observe high levels of divergence between the incipient species H. erato and H. himera, suggesting that divergence may accumulate early in the speciation process. Comparisons of genomic divergence between the incipient species and allopatric races suggest that limited gene flow cannot account for the observed high levels of divergence between the incipient species.

Conclusions

Our results provide a reconstruction of the speciation continuum across the H. erato clade and provide insights into the processes that drive genomic divergence during speciation, establishing the H. erato clade as a powerful framework for the study of speciation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0486-y) contains supplementary material, which is available to authorized users.

The diversification of Heliconius butterflies: what have we learned in 150 years?

Research into Heliconius butterflies has made a significant contribution to evolutionary biology. Here, we review our understanding of the diversification of these butterflies, covering recent advances and a vast foundation of earlier work. Whereas no single group of organisms can be sufficient for understanding life's diversity, after years of intensive study, research into Heliconius has addressed a wide variety of evolutionary questions. We first discuss evidence for widespread gene flow between Heliconius species and what this reveals about the nature of species. We then address the evolution and diversity of warning patterns, both as the target of selection and with respect to their underlying genetic basis. The identification of major genes involved in mimetic shifts, and homology at these loci between distantly related taxa, has revealed a surprising predictability in the genetic basis of evolution. In the final sections, we consider the evolution of warning patterns, and Heliconius diversity more generally, within a broader context of ecological and sexual selection. We consider how different traits and modes of selection can interact and influence the evolution of reproductive isolation.

Genetic Changes to a Transcriptional Silencer Element Confers Phenotypic Diversity within and between Drosophila Species

The modification of transcriptional regulation has become increasingly appreciated as a major contributor to morphological evolution. However, the role of negative-acting control elements (e.g. silencers) in generating morphological diversity has been generally overlooked relative to positive-acting “enhancer” elements. The highly variable body coloration patterns among Drosophilid insects represents a powerful model system in which the molecular alterations that underlie phenotypic diversity can be defined. In a survey of pigment phenotypes among geographically disparate Japanese populations of Drosophila auraria, we discovered a remarkable degree of variation in male-specific abdominal coloration. In testing the expression patterns of the major pigment-producing enzymes, we found that phenotypes uniquely correlated with differences in the expression of ebony, a gene required for yellow-colored cuticle. Assays of ebony’s transcriptional control region indicated that a lightly pigmented strain harbored cis-regulatory mutations that caused correlated changes in its expression. Through a series of chimeric reporter constructs between light and dark strain alleles, we localized function-altering mutations to a conserved silencer that mediates a male-specific pattern of ebony repression. This suggests that the light allele was derived through the loss of this silencer’s activity. Furthermore, examination of the ebony gene of D. serrata, a close relative of D. auraria which secondarily lost male-specific pigmentation revealed the parallel loss of this silencer element. These results demonstrate how loss-of-function mutations in a silencer element resulted in increased gene expression. We propose that the mutational inactivation of silencer elements may represent a favored path to evolve gene expression, impacting morphological traits.

One of the greatest challenges in understanding the relationship between genotype and phenotype is to discern how changes in DNA affect the normal functioning of genes. Mutations may generate a new function for a gene, yet it is frequently observed that they inactivate some aspect of a gene’s normal capacity. Investigations focused on understanding the developmental basis for the evolution of anatomical structures has found a prevalent role for mutations that alter developmental gene regulation. In animals, genes are transcriptionally activated in specific tissues during development by regulatory sequences distributed across their expansive non-protein coding regions. Regulatory elements known as silencers act to prevent genes from being expressed in certain tissues, providing a mechanism for precise control. Here, we show how a silencer that prevents expression of a pigment-producing enzyme in certain Drosophila species has repeatedly been subject to inactivating mutations that increased this gene’s expression. This example illustrates how such negative-acting regulatory sequences can represent a convenient target for increasing gene expression through the loss of a genetic element.

Gain of cis-regulatory activities underlies novel domains of wingless gene expression in Drosophila

Changes in gene expression during animal development are largely responsible for the evolution of morphological diversity. However, the genetic and molecular mechanisms responsible for the origins of new gene-expression domains have been difficult to elucidate. Here, we sought to identify molecular events underlying the origins of three novel features of wingless (wg) gene expression that are associated with distinct pigmentation patterns in Drosophila guttifera. We compared the activity of cis-regulatory sequences (enhancers) across the wg locus in D. guttifera and Drosophila melanogaster and found strong functional conservation among the enhancers that control similar patterns of wg expression in larval imaginal discs that are essential for appendage development. For pupal tissues, however, we found three novel wg enhancer activities in D. guttifera associated with novel domains of wg expression, including two enhancers located surprisingly far away in an intron of the distant Wnt10 gene. Detailed analysis of one enhancer (the vein-tip enhancer) revealed that it overlapped with a region controlling wg expression in wing crossveins (crossvein enhancer) in D. guttifera and other species. Our results indicate that one novel domain of wg expression in D. guttifera wings evolved by co-opting pre-existing regulatory sequences governing gene activity in the developing wing. We suggest that the modification of existing enhancers is a common path to the evolution of new gene-expression domains and enhancers.

New genomes clarify mimicry evolution

For over 100 years, it has been known that polymorphic mimicry is often switched by simple mendelian factors, yet the physical nature of these loci had escaped characterization. Now, the genome sequences of two swallowtail butterfly (Papilio) species have enabled the precise identification of a locus underlying mimicry, adding to unprecedented recent discoveries in mimicry genetics.

Unraveling the Tangled Skein: The Evolution of Transcriptional Regulatory Networks in Development

The molecular and genetic basis for the evolution of anatomical diversity is a major question that has inspired evolutionary and developmental biologists for decades. Because morphology takes form during development, a true comprehension of how anatomical structures evolve requires an understanding of the evolutionary events that alter developmental genetic programs. Vast gene regulatory networks (GRNs) that connect transcription factors to their target regulatory sequences control gene expression in time and space and therefore determine the tissue-specific genetic programs that shape morphological structures. In recent years, many new examples have greatly advanced our understanding of the genetic alterations that modify GRNs to generate newly evolved morphologies. Here, we review several aspects of GRN evolution, including their deep preservation, their mechanisms of alteration, and how they originate to generate novel developmental programs.

Towards the identification of the loci of adaptive evolution

1. Establishing the genetic and molecular basis underlying adaptive traits is one of the major goals of evolutionary geneticists in order to understand the connection between genotype and phenotype and elucidate the mechanisms of evolutionary change. Despite considerable effort to address this question, there remain relatively few systems in which the genes shaping adaptations have been identified. 2. Here, we review the experimental tools that have been applied to document the molecular basis underlying evolution in several natural systems, in order to highlight their benefits, limitations and suitability. In most cases, a combination of DNA, RNA and functional methodologies with field experiments will be needed to uncover the genes and mechanisms shaping adaptation in nature.

Conservatism and novelty in the genetic architecture of adaptation in Heliconius butterflies

Understanding the genetic architecture of adaptive traits has been at the centre of modern evolutionary biology since Fisher; however, evaluating how the genetic architecture of ecologically important traits influences their diversification has been hampered by the scarcity of empirical data. Now, high-throughput genomics facilitates the detailed exploration of variation in the genome-to-phenotype map among closely related taxa. Here, we investigate the evolution of wing pattern diversity in Heliconius, a clade of neotropical butterflies that have undergone an adaptive radiation for wing-pattern mimicry and are influenced by distinct selection regimes. Using crosses between natural wing-pattern variants, we used genome-wide restriction site-associated DNA (RAD) genotyping, traditional linkage mapping and multivariate image analysis to study the evolution of the architecture of adaptive variation in two closely related species: Heliconius hecale and H. ismenius. We implemented a new morphometric procedure for the analysis of whole-wing pattern variation, which allows visualising spatial heatmaps of genotype-to-phenotype association for each quantitative trait locus separately. We used the H. melpomene reference genome to fine-map variation for each major wing-patterning region uncovered, evaluated the role of candidate genes and compared genetic architectures across the genus. Our results show that, although the loci responding to mimicry selection are highly conserved between species, their effect size and phenotypic action vary throughout the clade. Multilocus architecture is ancestral and maintained across species under directional selection, whereas the single-locus (supergene) inheritance controlling polymorphism in H. numata appears to have evolved only once. Nevertheless, the conservatism in the wing-patterning toolkit found throughout the genus does not appear to constrain phenotypic evolution towards local adaptive optima.

How do regulatory networks evolve and expand throughout evolution?

Throughout evolution, regulatory networks need to expand and adapt to accommodate novel genes and gene functions. However, the molecular details explaining how gene networks evolve remain largely unknown. Recent studies demonstrate that changes in transcription factors contribute to the evolution of regulatory networks. In particular, duplication of transcription factors followed by specific mutations in their DNA-binding or interaction domains propels the divergence and emergence of new networks. The innate promiscuity and modularity of regulatory networks contributes to their evolvability: duplicated promiscuous regulators and their target promoters can acquire mutations that lead to gradual increases in specificity, allowing neofunctionalization or subfunctionalization.

Inter and intra island introgression in a wolf spider radiation from the Galápagos, and its implications for parallel evolution

Parallel radiations within island systems are often assumed to follow a simple scenario in which single colonization events are followed by in situ adaptive divergence. However, subsequent gene exchange after the initial colonization and during the divergence process might have important evolutionary impacts on species radiations. Gene exchange among ecologically similar species from different islands may lead to introgression of adaptive genetic variation and influence the parallel divergence process. In this study, we estimate levels of gene exchange within a wolf spider radiation of the genus Hogna Simon, 1885, from the Galápagos, wherein habitat specialization into 'high elevation' and 'coastal dry' species apparently evolved repeatedly on two islands. By using a multilocus approach we show that low levels of inter-island and relatively higher levels of intra island introgression shaped genetic variation in this species complex. Using these estimates, we demonstrate by means of a coalescence simulation that under these inter- and intra-island migration rates parallel evolution most likely evolves by introgression of adaptive alleles among islands, rather than through independent mutations despite the close genetic relationship of species within islands. As species phylogenies within radiations are frequently used to infer the divergence pattern, even relatively low levels of interspecific gene flow should not be neglected when interpreting parallel trait evolution.

[The genetic walk of evolution]

Genetic mutations are the main fuel of evolution. In each generation, they produce new variations, which may be sorted out by natural or sexual selection. Mutations are generated by chance; yet which are the mutations actually sorted out by evolution, and why? This review presents some recent advances regarding this question. First, we gather results obtained at molecular and cellular levels, through synthetic experiments and under artificial selection paradigms. Next, we highlight studies at the multi-cellular level, especially studies of repeated evolution, whereby independent lineages acquire similar traits. Recent meta-analysis and quantifications are being presented; together they suggest that evolutionary relevant mutations accumulate around hotspots, spanning different levels of genetic organization. Pioneering work suggests that many causes, corresponding to many biological contexts, may explain the existence of these genetic hotspots. We finally discuss methodological limits, empirical challenges and a few future potential directions for this domain of research dedicated to the genetic path of evolution.

Gradual and contingent evolutionary emergence of leaf mimicry in butterfly wing patterns

Background

Special resemblance of animals to natural objects such as leaves provides a representative example of evolutionary adaptation. The existence of such sophisticated features challenges our understanding of how complex adaptive phenotypes evolved. Leaf mimicry typically consists of several pattern elements, the spatial arrangement of which generates the leaf venation-like appearance. However, the process by which leaf patterns evolved remains unclear.

Results

In this study we show the evolutionary origin and process for the leaf pattern in Kallima (Nymphalidae) butterflies. Using comparative morphological analyses, we reveal that the wing patterns of Kallima and 45 closely related species share the same ground plan, suggesting that the pattern elements of leaf mimicry have been inherited across species with lineage-specific changes of their character states. On the basis of these analyses, phylogenetic comparative methods estimated past states of the pattern elements and enabled reconstruction of the wing patterns of the most recent common ancestor. This analysis shows that the leaf pattern has evolved through several intermediate patterns. Further, we use Bayesian statistical methods to estimate the temporal order of character-state changes in the pattern elements by which leaf mimesis evolved, and show that the pattern elements changed their spatial arrangement (e.g., from a curved line to a straight line) in a stepwise manner and finally establish a close resemblance to a leaf venation-like appearance.

Conclusions

Our study provides the first evidence for stepwise and contingent evolution of leaf mimicry.　 Leaf mimicry patterns evolved in a gradual, rather than a sudden, manner from a non-mimetic ancestor. Through a lineage of Kallima butterflies, the leaf patterns evolutionarily originated through temporal accumulation of orchestrated changes in multiple pattern elements.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0229-5) contains supplementary material, which is available to authorized users.

Stable Heliconius butterfly hybrid zones are correlated with a local rainfall peak at the edge of the Amazon basin

Multilocus clines between Müllerian mimetic races of Heliconius butterflies provide a classic example of the maintenance of hybrid zones and their importance in speciation. Concordant hybrid zones in the mimics Heliconius erato and H. melpomene in northern Peru were carefully documented in the 1980s, and this prior work now permits a historical analysis of the movement or stasis of the zones. Previous work predicted that these zones might be moving toward the Andes due to selective asymmetry. Extensive deforestation and climate change might also be expected to affect the positions and widths of the hybrid zones. We show that the positions and shapes of these hybrid zones have instead remained remarkably stable between 1985 and 2012. The stability of this interaction strongly implicates continued selection, rather than neutral mixing following secondary contact. The stability of cline widths and strong linkage disequilibria (gametic correlation coefficients Rmax = 0.35-0.56 among unlinked loci) over 25 years suggest that mimetic selection pressures on each color pattern locus have remained approximately constant (s ≈ 0.13-0.40 per locus in both species). Exceptionally high levels of precipitation at the edge of the easternmost Andes may act as a population density trough for butterflies, trapping the hybrid zones at the foot of the mountains, and preventing movement. As such, our results falsify one prediction of the Pleistocene Refugium theory: That the ranges of divergent species or subspecies should be centered on regions characterized by maxima of rainfall, with hybrid zones falling in more arid regions between them.

Müllerian mimicry as a result of codivergence between velvet ants and spider wasps

Recent studies have delineated a large Nearctic Müllerian mimicry complex in Dasymutilla velvet ants. Psorthaspis spider wasps live in areas where this mimicry complex is found and are phenotypically similar to Dasymutilla. We tested the idea that Psorthaspis spider wasps are participating in the Dasymutilla mimicry complex and that they codiverged with Dasymutilla. We performed morphometric analyses and human perception tests, and tabulated distributional records to determine the fit of Psorthaspis to the Dasymutilla mimicry complex. We inferred a dated phylogeny using nuclear molecular markers (28S, elongation factor 1-alpha, long-wavelength rhodopsin and wingless) for Psorthaspis species and compared it to a dated phylogeny of Dasymutilla. We tested for codivergence between the two groups using two statistical analyses. Our results show that Psorthaspis spider wasps are morphologically similar to the Dasymutilla mimicry rings. In addition, our tests indicate that Psorthaspis and Dasymutilla codiverged to produce similar color patterns. This study expands the breadth of the Dasymutilla Müllerian mimicry complex and provides insights about how codivergence influenced the evolution of mimicry in these groups.

Contrasting modes and tempos of venom expression evolution in two snake species

Selection is predicted to drive diversification within species and lead to local adaptation, but understanding the mechanistic details underlying this process and thus the genetic basis of adaptive evolution requires the mapping of genotype to phenotype. Venom is complex and involves many genes, but the specialization of the venom gland toward toxin production allows specific transcripts to be correlated with specific toxic proteins, establishing a direct link from genotype to phenotype. To determine the extent of expression variation and identify the processes driving patterns of phenotypic diversity, we constructed genotype-phenotype maps and compared range-wide toxin-protein expression variation for two species of snake with nearly identical ranges: the eastern diamondback rattlesnake (Crotalus adamanteus) and the eastern coral snake (Micrurus fulvius). We detected significant expression variation in C. adamanteus, identified the specific loci associated with population differentiation, and found that loci expressed at all levels contributed to this divergence. Contrary to expectations, we found no expression variation in M. fulvius, suggesting that M. fulvius populations are not locally adapted. Our results not only linked expression variation at specific loci to divergence in a polygenic, complex trait but also have extensive conservation and biomedical implications. C. adamanteus is currently a candidate for federal listing under the Endangered Species Act, and the loss of any major population would result in the irrevocable loss of a unique venom phenotype. The lack of variation in M. fulvius has significant biomedical application because our data will assist in the development of effective antivenom for this species.

Predictable transcriptome evolution in the convergent and complex bioluminescent organs of squid

Despite contingency in life's history, the similarity of evolutionarily convergent traits may represent predictable solutions to common conditions. However, the extent to which overall gene expression levels (transcriptomes) underlying convergent traits are themselves convergent remains largely unexplored. Here, we show strong statistical support for convergent evolutionary origins and massively parallel evolution of the entire transcriptomes in symbiotic bioluminescent organs (bacterial photophores) from two divergent squid species. The gene expression similarities are so strong that regression models of one species' photophore can predict organ identity of a distantly related photophore from gene expression levels alone. Our results point to widespread parallel changes in gene expression evolution associated with convergent origins of complex organs. Therefore, predictable solutions may drive not only the evolution of novel, complex organs but also the evolution of overall gene expression levels that underlie them.

A nonsynonymous mutation in the transcriptional regulator lbh is associated with cichlid craniofacial adaptation and neural crest cell development

Since the time of Darwin, biologists have sought to understand the origins and maintenance of life's diversity of form. However, the nature of the exact DNA mutations and molecular mechanisms that result in morphological differences between species remains unclear. Here, we characterize a nonsynonymous mutation in a transcriptional coactivator, limb bud and heart homolog (lbh), which is associated with adaptive variation in the lower jaw of cichlid fishes. Using both zebrafish and Xenopus, we demonstrate that lbh mediates migration of cranial neural crest cells, the cellular source of the craniofacial skeleton. A single amino acid change that is alternatively fixed in cichlids with differing facial morphologies results in discrete shifts in migration patterns of this multipotent cell type that are consistent with both embryological and adult craniofacial phenotypes. Among animals, this polymorphism in lbh represents a rare example of a coding change that is associated with continuous morphological variation. This work offers novel insights into the development and evolution of the craniofacial skeleton, underscores the evolutionary potential of neural crest cells, and extends our understanding of the genetic nature of mutations that underlie divergence in complex phenotypes.

Duplication of a promiscuous transcription factor drives the emergence of a new regulatory network

The emergence of new genes throughout evolution requires rewiring and extension of regulatory networks. However, the molecular details of how the transcriptional regulation of new gene copies evolves remain largely unexplored. Here we show how duplication of a transcription factor gene allowed the emergence of two independent regulatory circuits. Interestingly, the ancestral transcription factor was promiscuous and could bind different motifs in its target promoters. After duplication, one paralogue evolved increased binding specificity so that it only binds one type of motif, whereas the other copy evolved a decreased activity so that it only activates promoters that contain multiple binding sites. Interestingly, only a few mutations in both the DNA-binding domains and in the promoter binding sites were required to gradually disentangle the two networks. These results reveal how duplication of a promiscuous transcription factor followed by concerted cis and trans mutations allows expansion of a regulatory network.

Ancient homology underlies adaptive mimetic diversity across butterflies

Convergent evolution provides a rare, natural experiment with which to test the predictability of adaptation at the molecular level. Little is known about the molecular basis of convergence over macro-evolutionary timescales. Here we use a combination of positional cloning, population genomic resequencing, association mapping and developmental data to demonstrate that positionally orthologous nucleotide variants in the upstream region of the same gene, WntA, are responsible for parallel mimetic variation in two butterfly lineages that diverged >65 million years ago. Furthermore, characterization of spatial patterns of WntA expression during development suggests that alternative regulatory mechanisms underlie wing pattern variation in each system. Taken together, our results reveal a strikingly predictable molecular basis for phenotypic convergence over deep evolutionary time.

Little is known about the genetic basis of convergent evolution in deeply diverged species. Here, the authors show that variation in the WntA gene is associated with parallel wing pattern variation in two butterflies that diverged more than 65 million years ago.

Wnt signaling underlies evolution and development of the butterfly wing pattern symmetry systems

Most butterfly wing patterns are proposed to be derived from a set of conserved pattern elements known as symmetry systems. Symmetry systems are so-named because they are often associated with parallel color stripes mirrored around linear organizing centers that run between the anterior and posterior wing margins. Even though the symmetry systems are the most prominent and diverse wing pattern elements, their study has been confounded by a lack of knowledge regarding the molecular basis of their development, as well as the difficulty of drawing pattern homologies across species with highly derived wing patterns. Here we present the first molecular characterization of symmetry system development by showing that WntA expression is consistently associated with the major basal, discal, central, and external symmetry system patterns of nymphalid butterflies. Pharmacological manipulations of signaling gradients using heparin and dextran sulfate showed that pattern organizing centers correspond precisely with WntA, wingless, Wnt6, and Wnt10 expression patterns, thus suggesting a role for Wnt signaling in color pattern induction. Importantly, this model is supported by recent genetic and population genomic work identifying WntA as the causative locus underlying wing pattern variation within several butterfly species. By comparing the expression of WntA between nymphalid butterflies representing a range of prototypical symmetry systems, slightly deviated symmetry systems, and highly derived wing patterns, we were able to infer symmetry system homologies in several challenging cases. Our work illustrates how highly divergent morphologies can be derived from modifications to a common ground plan across both micro- and macro-evolutionary time scales.

Genetic mapping of horizontal stripes in Lake Victoria cichlid fishes: benefits and pitfalls of using RAD markers for dense linkage mapping

The genetic dissection of naturally occurring phenotypes sheds light on many fundamental and longstanding questions in speciation and adaptation and is a central research topic in evolutionary biology. Until recently, forward-genetic approaches were virtually impossible to apply to nonmodel organisms, but the development of next-generation sequencing techniques eases this difficulty. Here, we use the ddRAD-seq method to map a colour trait with a known adaptive function in cichlid fishes, well-known textbook examples for rapid rates of speciation and astonishing phenotypic diversification. A suite of phenotypic key innovations is related to speciation and adaptation in cichlids, among which body coloration features prominently. The focal trait of this study, horizontal stripes, evolved in parallel in several cichlid radiations and is associated with piscivorous foraging behaviour. We conducted interspecific crosses between Haplochromis sauvagei and H. nyererei and constructed a linkage map with 867 SNP markers distributed on 22 linkage groups and total size of 1130.63 cM. Lateral stripes are inherited as a Mendelian trait and map to a single genomic interval that harbours a paralog of a gene with known function in stripe patterning. Dorsolateral and mid-lateral stripes were always coinherited and are thus under the same genetic control. Additionally, we directly quantify the genotyping error rates in RAD markers and offer guidelines for identifying and dealing with errors. Uncritical marker selection was found to severely impact linkage map construction. Fortunately, by applying appropriate quality control steps, a genotyping accuracy of >99.9% can be reached, thus allowing for efficient linkage mapping of evolutionarily relevant traits.

Embryogenesis of Heliconius erato (Lepidoptera, Nymphalidae): a contribution to the anatomical development of an evo-devo model organism

This study reports on the embryogenesis of Heliconius erato phyllis between blastoderm formation and the prehatching larval stage. Syncytial blastoderm formation occurred approximately 2 h after egg laying (AEL) and at about 4 h, the cellular blastoderm was formed. The germ band arose from the entire length of the blastoderm, and rapidly became compacted occupying approximately two-thirds of the egg length. At about 7 h AEL, protocephalon and protocorm differentiation occurred. Continued proliferation of the germ band was followed by penetration into the yolk mass, forming a C-shaped embryo at about 10 h. Approximately 12 h AEL, the gnathal, thoracic and abdominal segments became visible. The primordium of the mouthparts and thoracic legs formed as paired evaginations, while the prolegs formed as paired lobes. At about 30 h, the embryo reversed dorsoventrally. Approximately 32 h AEL, the protocephalon and gnathal segments fused, shifting the relative position of the rudimentary appendages in this region. At about 52 h, the embryo was U-shaped in lateral view and at approximately 56 h, the bristles began evagination from the larval cuticle. Larvae hatched at about 72 h. We found that H. erato phyllis followed an embryonic pattern consistent with long-germ embryogenesis. Thus, we believe that H. erato phyllis should be classified as a long-germ lepidopteran. The study of H. erato phyllis embryogenesis provided a structural glimpse into the morphogenetic events that occur in the Heliconius egg period. This study could help future molecular approaches to understanding the evolution of Heliconius development.

Different tradeoffs result from alternate genetic adaptations to a common environment

Fitness tradeoffs are often assumed by evolutionary theory, yet little is known about the frequency of fitness tradeoffs during stress adaptation. Even less is known about the genetic factors that confer these tradeoffs and whether alternative adaptive mutations yield contrasting tradeoff dynamics. We addressed these issues using 114 clones of Escherichia coli that were evolved independently for 2,000 generations under thermal stress (42.2 °C). For each clone, we measured their fitness relative to the ancestral clone at 37 °C and 20 °C. Tradeoffs were common at 37 °C but more prevalent at 20 °C, where 56% of clones were outperformed by the ancestor. We also characterized the upper and lower thermal boundaries of each clone. All clones shifted their upper boundary to at least 45 °C; roughly half increased their lower niche boundary concomitantly, representing a shift of thermal niche. The remaining clones expanded their thermal niche by increasing their upper limit without a commensurate increase of lower limit. We associated these niche dynamics with genotypes and confirmed associations by engineering single mutations in the rpoB gene, which encodes the beta subunit of RNA polymerase, and the rho gene, which encodes a termination factor. Single mutations in the rpoB gene exhibit antagonistic pleiotropy, with fitness tradeoffs at 18 °C and fitness benefits at 42.2 °C. In contrast, a mutation within the rho transcriptional terminator, which defines an alternative adaptive pathway from that of rpoB, had no demonstrable effect on fitness at 18 °C. This study suggests that two different genetic pathways toward high-temperature adaptation have contrasting effects with respect to thermal tradeoffs.

Pigmentary and photonic coloration mechanisms reveal taxonomic relationships of the Cattlehearts (Lepidoptera: Papilionidae: Parides)

BACKGROUND

The colorful wing patterns of butterflies, a prime example of biodiversity, can change dramatically within closely related species. Wing pattern diversity is specifically present among papilionid butterflies. Whether a correlation between color and the evolution of these butterflies exists so far remained unsolved.

RESULTS

We here investigate the Cattlehearts, Parides, a small Neotropical genus of papilionid butterflies with 36 members, the wings of which are marked by distinctly colored patches. By applying various physical techniques, we investigate the coloration toolkit of the wing scales. The wing scales contain two different, wavelength-selective absorbing pigments, causing pigmentary colorations. Scale ridges with multilayered lamellae, lumen multilayers or gyroid photonic crystals in the scale lumen create structural colors that are variously combined with these pigmentary colors.

CONCLUSIONS

The pigmentary and structural traits strongly correlate with the taxonomical distribution of Parides species. The experimental findings add crucial insight into the evolution of butterfly wing scales and show the importance of morphological parameter mapping for butterfly phylogenetics.

Gene co-expression modules underlying polymorphic and monomorphic zooids in the colonial hydrozoan, Hydractinia symbiolongicarpus

Advances in sequencing technology have forced a quantitative revolution in Evolutionary Biology. One important feature of this renaissance is that comprehensive genomic resources can be obtained quickly for almost any taxon, thus speeding the development of new model organisms. Here, we analyze 20 RNA-seq libraries from morphologically, sexually, and genetically distinct polyp types from the gonochoristic colonial hydrozoan, Hydractinia symbiolongicarpus (Cnidaria). Analyses of these data using weighted gene co-expression networks highlight deeply conserved genetic elements of animal spermatogenesis and demonstrate the utility of these methods in identifying modules of genes that correlate with different zooid types across various statistical contrasts. RNA-seq data and analytical scripts described here are deposited in publicly available databases.

Dynamics of F-actin prefigure the structure of butterfly wing scales

The wings of butterflies and moths consist of dorsal and ventral epidermal surfaces that give rise to overlapping layers of scales and hairs (Lepidoptera, "scale wing"). Wing scales (average length ~200 µm) are homologous to insect bristles (macrochaetes), and their colors create the patterns that characterize lepidopteran wings. The topology and surface sculpture of wing scales vary widely, and this architectural complexity arises from variations in the developmental program of the individual scale cells of the wing epithelium. One of the more striking features of lepidopteran wing scales are the longitudinal ridges that run the length of the mature (dead) cell, gathering the cuticularized scale cell surface into pleats on the sides of each scale. While also present around the periphery of other insect bristles and hairs, longitudinal ridges in lepidopteran wing scales gain new significance for their creation of iridescent color through microribs and lamellae. Here we show the dynamics of the highly organized F-actin filaments during scale cell development, and present experimental manipulations of actin polymerization that reveal the essential role of this cytoskeletal component in wing scale elongation and the positioning of longitudinal ribs.

The expansion of body coloration involves coordinated evolution in cis and trans within the pigmentation regulatory network of Drosophila prostipennis

The generation of complex morphological features requires the precisely orchestrated expression of numerous genes during development. While several traits have been resolved to evolutionary changes within a single gene, the evolutionary path by which genes derive co-localized or mutually excluded expression patterns is currently a mystery. Here we investigate how the Drosophila pigmentation gene network was altered in Drosophila prostipennis, a species in the Drosophila melanogaster subgroup, that evolved expanded abdominal pigmentation. We show that this expansion involved broadened expression of the melanin-promoting enzyme genes tan and yellow, and a reciprocal withdrawn pattern of the melanin-suppressing enzyme gene ebony. To examine whether these coordinated changes to the network were generated through mutations in the cis-regulatory elements (CREs) of these genes, we cloned and tested CREs of D. prostipennis tan, ebony, and yellow in transgenic reporter assays. Regulatory regions of both tan and ebony failed to recapitulate the derived D. prostipennis expression phenotype, implicating the modification of a factor or factors upstream of both genes. However, the D. prostipennis yellow cis-regulatory region recapitulated the expanded expression pattern observed in this species, implicating causative mutations in cis to yellow. Our results provide an example in which a coordinated expression program evolved through independent changes at multiple loci, rather than through changes to a single "master regulator" directing a suite of downstream target genes. This implies a complex network structure in which each gene may be subject to a unique set of inputs, and resultantly may require individualized evolutionary paths to yield correlated gene expression patterns.

The evolutionary genomics of cichlid fishes: explosive speciation and adaptation in the postgenomic era

With more than 1,500 species, cichlid fishes provide textbook examples of recent and diverse adaptive radiations, rapid rates of speciation, and the parallel evolution of adaptive phenotypes among both recently and distantly related lineages. This extraordinary diversity has attracted considerable interest from researchers across several biological disciplines. Their broad phenotypic variation coupled with recent divergence makes cichlids an ideal model system for understanding speciation, adaptation, and phenotypic diversification. Genetic mapping, genome-wide analyses, and genome projects have flourished in the past decade and have added new insights on the question of why there are so many cichlids. These recent findings also show that the sharing of older DNA polymorphisms is extensive and suggest that linage sorting is incomplete and that adaptive introgression played a role in the African radiation. Here, we review the results of genetic and genomic research on cichlids in the past decade and suggest some potential avenues to further exploit the potential of the cichlid model system to provide a better understanding of the genomics of adaptation and speciation.

Comparative functional analyses of ultrabithorax reveal multiple steps and paths to diversification of legs in the adaptive radiation of semi-aquatic insects

Invasion of new ecological habitats is often associated with lineage diversification, yet the genetic changes underlying invasions and radiations are poorly understood. Over 200 million years ago, the semi-aquatic insects invaded water surface from a common terrestrial ancestor and diversified to exploit a wide array of niches. Here, we uncover the changes in regulation and function of the gene Ultrabithorax associated with both the invasion of water surface and the subsequent diversification of the group. In the common ancestor of the semi-aquatic insects, a novel deployment of Ubx protein in the mid-legs increased their length, thereby enhancing their role in water surface walking. In derived lineages that specialize in rowing on the open water, additional changes in the timing of Ubx expression further elongated the mid-legs thereby facilitating their function as oars. In addition, Ubx protein function was selectively reversed to shorten specific rear-leg segments, thereby enabling their function as rudders. These changes in Ubx have generated distinct niche-specialized morphologies that account for the remarkable diversification of the semi-aquatic insects. Therefore, changes in the regulation and function of a key developmental gene may facilitate both the morphological change necessary to transition to novel habitats and fuel subsequent morphological diversification.

Population genomics of parallel hybrid zones in the mimetic butterflies, H. melpomene and H. erato

Hybrid zones can be valuable tools for studying evolution and identifying genomic regions responsible for adaptive divergence and underlying phenotypic variation. Hybrid zones between subspecies of Heliconius butterflies can be very narrow and are maintained by strong selection acting on color pattern. The comimetic species, H. erato and H. melpomene, have parallel hybrid zones in which both species undergo a change from one color pattern form to another. We use restriction-associated DNA sequencing to obtain several thousand genome-wide sequence markers and use these to analyze patterns of population divergence across two pairs of parallel hybrid zones in Peru and Ecuador. We compare two approaches for analysis of this type of data—alignment to a reference genome and de novo assembly—and find that alignment gives the best results for species both closely (H. melpomene) and distantly (H. erato, ∼15% divergent) related to the reference sequence. Our results confirm that the color pattern controlling loci account for the majority of divergent regions across the genome, but we also detect other divergent regions apparently unlinked to color pattern differences. We also use association mapping to identify previously unmapped color pattern loci, in particular the Ro locus. Finally, we identify a new cryptic population of H. timareta in Ecuador, which occurs at relatively low altitude and is mimetic with H. melpomene malleti.

Parallel developmental genetic features underlie stickleback gill raker evolution

Background

Convergent evolution, the repeated evolution of similar phenotypes in independent lineages, provides natural replicates to study mechanisms of evolution. Cases of convergent evolution might have the same underlying developmental and genetic bases, implying that some evolutionary trajectories might be predictable. In a classic example of convergent evolution, most freshwater populations of threespine stickleback fish have independently evolved a reduction of gill raker number to adapt to novel diets. Gill rakers are a segmentally reiterated set of dermal bones important for fish feeding. A previous large quantitative trait locus (QTL) mapping study using a marine × freshwater F2 cross identified QTL on chromosomes 4 and 20 with large effects on evolved gill raker reduction.

Results

By examining skeletal morphology in adult and developing sticklebacks, we find heritable marine/freshwater differences in gill raker number and spacing that are specified early in development. Using the expression of the Ectodysplasin receptor (Edar) gene as a marker of raker primordia, we find that the differences are present before the budding of gill rakers occurs, suggesting an early change to a lateral inhibition process controlling raker primordia spacing. Through linkage mapping in F2 fish from crosses with three independently derived freshwater populations, we find in all three crosses QTL overlapping both previously identified QTL on chromosomes 4 and 20 that control raker number. These two QTL affect the early spacing of gill raker buds.

Conclusions

Collectively, these data demonstrate that parallel developmental genetic features underlie the convergent evolution of gill raker reduction in freshwater sticklebacks, suggesting that even highly polygenic adaptive traits can have a predictable developmental genetic basis.

Genomics and the origin of species

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.