Mimicry and the evolution of premating isolation in Heliconius melpomene Linnaeus

Spatiotemporal dynamics in butterfly hybrid zones

Evaluating whether hybrid zones are stable or mobile can provide novel insights for evolution and conservation biology. Butterflies exhibit high sensitivity to environmental changes and represent an important model system for the study of hybrid zone origins and maintenance. Here, we review the literature exploring butterfly hybrid zones, with a special focus on their spatiotemporal dynamics and the potential mechanisms that could lead to their movement or stability. We then compare different lines of evidence used to investigate hybrid zone dynamics and discuss the strengths and weaknesses of each approach. Our goal with this review is to reveal general conditions associated with the stability or mobility of butterfly hybrid zones by synthesizing evidence obtained using different types of data sampled across multiple regions and spatial scales. Finally, we discuss spatiotemporal dynamics in the context of a speciation/divergence continuum, the relevance of hybrid zones for conservation biology, and recommend key topics for future investigation.

Quantifying visual acuity in Heliconius butterflies

Heliconius butterflies are well-known for their colourful wing patterns, which advertise distastefulness to potential predators and are used during mate choice. However, the relative importance of different aspects of these signals will depend on the visual abilities of Heliconius and their predators. Previous studies have investigated colour sensitivity and neural anatomy, but visual acuity (the ability to perceive detail) has not been studied in these butterflies. Here, we provide the first estimate of visual acuity in Heliconius: from a behavioural optomotor assay, we found that mean visual acuity = 0.49 cycles-per-degree (cpd), with higher acuity in males than females. We also examined eye morphology and report more ommatidia in male eyes. Finally, we estimated how visual acuity affects Heliconius visual perception compared to a potential avian predator. Whereas the bird predator maintained high resolving power, Heliconius lost the ability to resolve detail at greater distances, though colours may remain salient. These results will inform future studies of Heliconius wing pattern evolution, as well as other aspects in these highly visual butterflies, which have emerged as an important system in studies of adaptation and speciation.

Does sexual conflict contribute to the evolution of novel warning patterns?

Why warning patterns are so diverse is an enduring evolutionary puzzle. Because predators associate particular patterns with unpleasant experiences, an individual's predation risk should decrease as the local density of its warning pattern increases, promoting pattern monomorphism. Distasteful Heliconius butterflies are known for their diversity of warning patterns. Here, we explore whether interlocus sexual conflict can contribute to their diversification. Male Heliconius use warning patterns as mating cues, but mated females may suffer costs if this leads to disturbance, favouring novel patterns. Using simulations, we show that under our model conditions drift alone is unlikely to cause pattern diversification, but that sexual conflict can assist such a process. We also find that genetic architecture influences the evolution of male preferences, which track changes in warning pattern due to sexual selection. When male attraction imposes costs on females, this affects the speed at which novel pattern alleles increase. In two experiments, females laid fewer eggs with males present. However, although males in one experiment showed less interest in females with manipulated patterns, we found no evidence that female colouration mitigates sex-specific costs. Overall, male attraction to conspecific warning patterns may impose an unrecognized cost on Heliconius females, but further work is required to determine this experimentally.

Scale-dependent environmental effects on phenotypic distributions in Heliconius butterflies

Identifying the relative importance of different mechanisms responsible for the emergence and maintenance of phenotypic diversity can be challenging, as multiple selective pressures and stochastic events are involved in these processes. Therefore, testing how environmental conditions shape the distribution of phenotypes can offer important insights on local adaptation, divergence, and speciation. The red-yellow Müllerian mimicry ring of Heliconius butterflies exhibits a wide diversity of color patterns across the Neotropics and is involved in multiple hybrid zones, making it a powerful system to investigate environmental drivers of phenotypic distributions. Using the distantly related Heliconius erato and Heliconius melpomene co-mimics and a multiscale distribution approach, we investigated whether distinct phenotypes of these species are associated with different environmental conditions. We show that Heliconius red-yellow phenotypic distribution is strongly driven by environmental gradients (especially thermal and precipitation variables), but that phenotype and environment associations vary with spatial scale. While co-mimics are usually predicted to occur in similar environments at large spatial scales, patterns at local scales are not always consistent (i.e., different variables are best predictors of phenotypic occurrence in different locations) or congruent (i.e., co-mimics show distinct associations with environment). We suggest that large-scale analyses are important for identifying how environmental factors shape broad mimetic phenotypic distributions, but that local studies are essential to understand the context-dependent biotic, abiotic, and historical mechanisms driving finer-scale phenotypic transitions.

Convergent flight morphology among Müllerian mimic mutualists

Müllerian mimicry involves a signal mutualism between prey species, shaped by visually hunting predators, and recent work has emphasized the importance of color pattern. Predators respond to more than color pattern, however, and other traits are much less studied. This article examines the hypothesis of convergent evolution in flight-related morphology among eight mimicry complexes composed of 51 butterfly species (Nymphalidae, Danainae, Ithomiini) from a single community in Ecuador. Phylogenetic comparative analyses of 14 variables indicated strong morphological differences between mimicry complexes belonging to three clusters of morphological space ("large yellow transparent," "tiger," and "transparent"), not the eight predicted based on color pattern alone. Analyses found convergence within mimicry complexes, convergence between mimicry complexes within morphospace clusters, and divergence between mimicry complexes from different morphospace clusters. These three clusters differed in size, and body and wing shape, predicting that flight biomechanics also converge (i.e., locomotor mimicry). Potential constraints on evolution of morphological mimicry related to predator discrimination, and evolutionary rates, likely e xplain why flight-related morphology differences were limited to three clusters of morphological space. Finally, the added complexity that flight-related morphology brings to signals between predator and prey indicates that evolutionary switches in color pattern are not all equally likely, potentially limiting the evolution of color patterns if they do not match morphology.

Engaging the community in pollinator research: the effect of wing pattern and weather on butterfly behavior

Community science, which engages students and the public in data collection and scientific inquiry, is often integrated into conservation and long-term monitoring efforts. However, it has the potential to also introduce the public to, and be useful for, sensory ecology and other fields of study. Here we describe a community science project that exposes participants to animal behavior and sensory ecology using the rich butterfly community of Northwest Arkansas, USA. Butterflies use visual signals to communicate and to attract mates. Brighter colors can produce stronger signals for mate attraction but can also unintentionally attract negative attention from predators. Environmental conditions such as weather can affect visual signaling as well, by influencing the wavelengths of light available and subsequent signal detection. However, we do not know whether the signals butterflies present correlate broadly with how they behave. In this study, we collaborated with hundreds of students and community members at the University of Arkansas (UARK) and the Botanical Gardens of the Ozarks (BGO) for over 3.5 years to examine relationships among wing pattern, weather, time of day, behavior, and flower choice. We found that both weather and wing color influenced general butterfly behavior. Butterflies were seen feeding more on cloudy days than on sunny or partly cloudy days. Brown butterflies fed or sat more often, while white butterflies flew more often relative to other butterfly colors. We also found that there was an interaction between the effects of weather and wing color on butterfly behavior. Furthermore, butterfly color predicted the choice of flower colors that butterflies visited, though this effect was influenced by observer group (UARK student or BGO participant). These results suggest that flower choice may be associated with butterfly wing pattern, and that different environmental conditions may influence butterfly behavior in wing-pattern-specific ways. They also illustrate one way that public involvement in behavioral studies can facilitate the identification of coarse-scale, community-wide behavioral patterns, and lay the groundwork for future studies of sensory niches.

Visual mate preference evolution during butterfly speciation is linked to neural processing genes

Many animal species remain separate not because their individuals fail to produce viable hybrids but because they "choose" not to mate. However, we still know very little of the genetic mechanisms underlying changes in these mate preference behaviours. Heliconius butterflies display bright warning patterns, which they also use to recognize conspecifics. Here, we couple QTL for divergence in visual preference behaviours with population genomic and gene expression analyses of neural tissue (central brain, optic lobes and ommatidia) across development in two sympatric Heliconius species. Within a region containing 200 genes, we identify five genes that are strongly associated with divergent visual preferences. Three of these have previously been implicated in key components of neural signalling (specifically an ionotropic glutamate receptor and two regucalcins), and overall our candidates suggest shifts in behaviour involve changes in visual integration or processing. This would allow preference evolution without altering perception of the wider environment.

Transparency improves concealment in cryptically coloured moths

Predation is a ubiquitous and strong selective pressure on living organisms. Transparency is a predation defence widespread in water but rare on land. Some Lepidoptera display transparent patches combined with already cryptic opaque patches. A recent study showed that transparency reduced detectability of aposematic prey with conspicuous patches. However, whether transparency has any effect at reducing detectability of already cryptic prey is still unknown. We conducted field predation experiments with free avian predators where we monitored and compared survival of a fully opaque grey artificial form (cryptic), a form including transparent windows and a wingless artificial butterfly body. Survival of the transparent forms was similar to that of wingless bodies and higher than that of fully opaque forms, suggesting a reduction of detectability conferred by transparency. This is the first evidence that transparency decreases detectability in cryptic terrestrial prey. Future studies should explore the organization of transparent and opaque patches in animals and their interplay on survival, as well as the costs and other potential benefits associated with transparency on land.

Genetic dissection of assortative mating behavior

The evolution of new species is made easier when traits under divergent ecological selection are also mating cues. Such ecological mating cues are now considered more common than previously thought, but we still know little about the genetic changes underlying their evolution or more generally about the genetic basis for assortative mating behaviors. Both tight physical linkage and the existence of large-effect preference loci will strengthen genetic associations between behavioral and ecological barriers, promoting the evolution of assortative mating. The warning patterns of Heliconius melpomene and H. cydno are under disruptive selection due to increased predation of nonmimetic hybrids and are used during mate recognition. We carried out a genome-wide quantitative trait locus (QTL) analysis of preference behaviors between these species and showed that divergent male preference has a simple genetic basis. We identify three QTLs that together explain a large proportion (approximately 60%) of the difference in preference behavior observed between the parental species. One of these QTLs is just 1.2 (0-4.8) centiMorgans (cM) from the major color pattern gene optix, and, individually, all three have a large effect on the preference phenotype. Genomic divergence between H. cydno and H. melpomene is high but broadly heterogenous, and admixture is reduced at the preference-optix color pattern locus but not the other preference QTLs. The simple genetic architecture we reveal will facilitate the evolution and maintenance of new species despite ongoing gene flow by coupling behavioral and ecological aspects of reproductive isolation.

Recombination rate variation shapes barriers to introgression across butterfly genomes

Hybridisation and introgression can dramatically alter the relationships among groups of species, leading to phylogenetic discordance across the genome and between populations. Introgression can also erode species differences over time, but selection against introgression at certain loci acts to maintain postmating species barriers. Theory predicts that species barriers made up of many loci throughout the genome should lead to a broad correlation between introgression and recombination rate, which determines the extent to which selection on deleterious foreign alleles will affect neutral alleles at physically linked loci. Here, we describe the variation in genealogical relationships across the genome among three species of Heliconius butterflies: H. melpomene (mel), H. cydno (cyd), and H. timareta (tim), using whole genomes of 92 individuals, and ask whether this variation can be explained by heterogeneous barriers to introgression. We find that species relationships vary predictably at the chromosomal scale. By quantifying recombination rate and admixture proportions, we then show that rates of introgression are predicted by variation in recombination rate. This implies that species barriers are highly polygenic, with selection acting against introgressed alleles across most of the genome. In addition, long chromosomes, which have lower recombination rates, produce stronger barriers on average than short chromosomes. Finally, we find a consistent difference between two species pairs on either side of the Andes, which suggests differences in the architecture of the species barriers. Our findings illustrate how the combined effects of hybridisation, recombination, and natural selection, acting at multitudes of loci over long periods, can dramatically sculpt the phylogenetic relationships among species.

Aristaless Controls Butterfly Wing Color Variation Used in Mimicry and Mate Choice

Neotropical Heliconius butterflies display a diversity of warningly colored wing patterns, which serve roles in both Müllerian mimicry and mate choice behavior. Wing pattern diversity in Heliconius is controlled by a small number of unlinked, Mendelian "switch" loci [1]. One of these, termed the K locus, switches between yellow and white color patterns, important mimicry signals as well as mating cues [2-4]. Furthermore, mate preference behavior is tightly linked to this locus [4]. K controls the distribution of white versus yellow scales on the wing, with a dominant white allele and a recessive yellow allele. Here, we combine fine-scale genetic mapping, genome-wide association studies, gene expression analyses, population and comparative genomics, and genome editing with CRISPR/Cas9 to characterize the molecular basis of the K locus in Heliconius and to infer its evolutionary history. We show that white versus yellow color variation in Heliconius cydno is due to alternate haplotypes at a putative cis-regulatory element (CRE) downstream of a tandem duplication of the homeodomain transcription factor aristaless. Aristaless1 (al1) and aristaless2 (al2) are differentially regulated between white and yellow wings throughout development with elevated expression of al1 in developing white wings, suggesting a role in repressing pigmentation. Consistent with this, knockout of al1 causes white wings to become yellow. The evolution of wing color in this group has been marked by retention of the ancestral yellow color in many lineages, a single origin of white coloration in H. cydno, and subsequent introgression of white color from H. cydno into H. melpomene.

The appearance of mimetic Heliconius butterflies to predators and conspecifics

Adaptive coloration is under conflicting selection pressures: choosing potential mates and warning signaling against visually guided predators. Different elements of the color signal may therefore be tuned by evolution for different functions. We investigated how mimicry in four pairs of Heliconius comimics is potentially seen both from the perspective of butterflies and birds. Visual sensitivities of eight candidate avian predators were predicted through genetic analysis of their opsin genes. Using digital image color analysis, combined with bird and butterfly visual system models, we explored how predators and conspecifics may visualize mimetic patterns. Ultraviolet vision (UVS) birds are able to discriminate between the yellow and white colors of comimics better than violet vision (VS) birds. For Heliconius vision, males and females differ in their ability to discriminate comimics. Female vision and red filtering pigments have a significant effect on the perception of the yellow forewing band and the red ventral forewing pattern. A behavioral experiment showed that UV cues are used in mating behavior; removal of such cues was associated with an increased tendency to approach comimics as compared to conspecifics. We have therefore shown that visual signals can act to both reduce the cost of confusion in courtship and maintain the advantages of mimicry.

Testing the adaptive hypothesis of Batesian mimicry among hybridizing North American admiral butterflies

Batesian mimicry is characterized by phenotypic convergence between an unpalatable model and a palatable mimic. However, because convergent evolution may arise via alternative evolutionary mechanisms, putative examples of Batesian mimicry must be rigorously tested. Here, we used artificial butterfly facsimiles (N = 4000) to test the prediction that (1) palatable Limenitis lorquini butterflies should experience reduced predation when in sympatry with their putative model, Adelpha californica, (2) protection from predation on L. lorquini should erode outside of the geographical range of the model, and (3) mimetic color pattern traits are more variable in allopatry, consistent with relaxed selection for mimicry. We find support for these predictions, implying that this convergence is the result of selection for Batesian mimicry. Additionally, we conducted mark-recapture studies to examine the effect of mimicry and found that mimics survive significantly longer at sites where the model is abundant. Finally, in contrast to theoretical predictions, we found evidence that the Batesian model (A. californica) is protected from predation outside of its geographic range. We discuss these results considering the ongoing hybridization between L. lorquini and its sister species, L. weidemeyerii, and growing evidence that selection for mimicry predictably leads to a reduction in gene flow between nascent species.

Male sex pheromone components in Heliconius butterflies released by the androconia affect female choice

Sex-specific pheromones are known to play an important role in butterfly courtship, and may influence both individual reproductive success and reproductive isolation between species. Extensive ecological, behavioural and genetic studies of Heliconius butterflies have made a substantial contribution to our understanding of speciation. Male pheromones, although long suspected to play an important role, have received relatively little attention in this genus. Here, we combine morphological, chemical and behavioural analyses of male pheromones in the Neotropical butterfly Heliconius melpomene. First, we identify putative androconia that are specialized brush-like scales that lie within the shiny grey region of the male hindwing. We then describe putative male sex pheromone compounds, which are largely confined to the androconial region of the hindwing of mature males, but are absent in immature males and females. Finally, behavioural choice experiments reveal that females of H. melpomene, H. erato and H. timareta strongly discriminate against conspecific males which have their androconial region experimentally blocked. As well as demonstrating the importance of chemical signalling for female mate choice in Heliconius butterflies, the results describe structures involved in release of the pheromone and a list of potential male sex pheromone compounds.

Estimating the age of Heliconius butterflies from calibrated photographs

Mating behaviour and predation avoidance in Heliconius involve visual colour signals; however, there is considerable inter-individual phenotypic variation in the appearance of colours. In particular, the red pigment varies from bright crimson to faded red. It has been thought that this variation is primarily due to pigment fading with age, although this has not been explicitly tested. Previous studies have shown the importance of red patterns in mate choice and that birds and butterflies might perceive these small colour differences. Using digital photography and calibrated colour images, we investigated whether the hue variation in the forewing dorsal red band of Heliconius melpomene rosina corresponds with age. We found that the red hue and age were highly associated, suggesting that red colour can indeed be used as a proxy for age in the study of wild-caught butterflies.

Maintaining mimicry diversity: optimal warning colour patterns differ among microhabitats in Amazonian clearwing butterflies

Mimicry is one of the best-studied examples of adaptation, and recent studies have provided new insights into the role of mimicry in speciation and diversification. Classical Müllerian mimicry theory predicts convergence in warning signal among protected species, yet tropical butterflies are exuberantly diverse in warning colour patterns, even within communities. We tested the hypothesis that microhabitat partitioning in aposematic butterflies and insectivorous birds can lead to selection for different colour patterns in different microhabitats and thus help maintain mimicry diversity. We measured distribution across flight height and topography for 64 species of clearwing butterflies (Ithomiini) and their co-mimics, and 127 species of insectivorous birds, in an Amazon rainforest community. For the majority of bird species, estimated encounter rates were non-random for the two most abundant mimicry rings. Furthermore, most butterfly species in these two mimicry rings displayed the warning colour pattern predicted to be optimal for anti-predator defence in their preferred microhabitats. These conclusions were supported by a field trial using butterfly specimens, which showed significantly different predation rates on colour patterns in two microhabitats. We therefore provide the first direct evidence to support the hypothesis that different mimicry patterns can represent stable, community-level adaptations to differing biotic environments.

Poison frog color morphs express assortative mate preferences in allopatry but not sympatry

The concurrent divergence of mating traits and preferences is necessary for the evolution of reproductive isolation via sexual selection, and such coevolution has been demonstrated in diverse lineages. However, the extent to which assortative mate preferences are sufficient to drive reproductive isolation in nature is less clear. Natural contact zones between lineages divergent in traits and preferences provide exceptional opportunities for testing the predicted evolutionary consequences of such divergence. The strawberry poison frog (Oophaga pumilio) displays extreme color polymorphism in and around the young Bocas del Toro archipelago. In a transition zone between red and blue allopatric lineages, we asked whether female preferences diverged along with coloration, and whether any divergent preferences persist in a zone of sympatry. When choosing among red, blue and phenotypically intermediate males, females from monomorphic red and monomorphic blue populations both expressed assortative preferences. However, red, blue, and intermediate females from the contact zone all preferred red males, suggesting that divergent preferences may be insufficient to effect behavioral isolation. Our results highlight the complexity of behavioral isolation, and the need for studies that can reveal the circumstances under which divergent preferences do and do not contribute to speciation.

Population differentiation at a regional scale in spadefoot toads: contributions of distance and divergent selective environments

The causes of population differentiation can provide insight into the origins of early barriers to gene flow. Two key drivers of population differentiation are geographic distance and local adaptation to divergent selective environments. When reproductive isolation arises because some populations of a species are under selection to avoid hybridization while others are not, population differentiation and even speciation can result. Spadefoot toad populations Spea multiplicata that are sympatric with a congener have undergone reinforcement. This reinforcement has resulted not only in increased reproductive isolation from the congener, but also in the evolution of reproductive isolation from nearby and distant conspecific allopatric populations. We used multiple approaches to evaluate the contributions of geographic distance and divergent selective environments to population structure across this regional scale in S. multiplicata, based on genotypes from six nuclear microsatellite markers. We compared groups of populations varying in both geographic location and in the presence of a congener. Hierarchical F-statistics and results from cluster analyses and discriminant analyses of principal components all indicate that geographic distance is the stronger contributor to genetic differentiation among S. multiplicata populations at a regional scale. However, we found evidence that adaptation to divergent selective environments also contributes to population structure. Our findings highlight how variation in the balance of evolutionary forces acting across a species' range can lead to variation in the relative contributions of geographic distance and local adaptation to population differentiation across different spatial scales.

The functional basis of wing patterning in Heliconius butterflies: the molecules behind mimicry

Wing-pattern mimicry in butterflies has provided an important example of adaptation since Charles Darwin and Alfred Russell Wallace proposed evolution by natural selection >150 years ago. The neotropical butterfly genus Heliconius played a central role in the development of mimicry theory and has since been studied extensively in the context of ecology and population biology, behavior, and mimicry genetics. Heliconius species are notable for their diverse color patterns, and previous crossing experiments revealed that much of this variation is controlled by a small number of large-effect, Mendelian switch loci. Recent comparative analyses have shown that the same switch loci control wing-pattern diversity throughout the genus, and a number of these have now been positionally cloned. Using a combination of comparative genetic mapping, association tests, and gene expression analyses, variation in red wing patterning throughout Heliconius has been traced back to the action of the transcription factor optix. Similarly, the signaling ligand WntA has been shown to control variation in melanin patterning across Heliconius and other butterflies. Our understanding of the molecular basis of Heliconius mimicry is now providing important insights into a variety of additional evolutionary phenomena, including the origin of supergenes, the interplay between constraint and evolvability, the genetic basis of convergence, the potential for introgression to facilitate adaptation, the mechanisms of hybrid speciation in animals, and the process of ecological speciation.

Mimetic Divergence and the Speciation Continuum in the Mimic Poison Frog Ranitomeya imitator

While divergent ecological adaptation can drive speciation, understanding the factors that facilitate or constrain this process remains a major goal in speciation research. Here, we study two mimetic transition zones in the poison frog Ranitomeya imitator, a species that has undergone a Müllerian mimetic radiation to establish four morphs in Peru. We find that mimetic morphs are strongly phenotypically differentiated, producing geographic clines with varying widths. However, distinct morphs show little neutral genetic divergence, and landscape genetic analyses implicate isolation by distance as the primary determinant of among-population genetic differentiation. Mate choice experiments suggest random mating at the transition zones, although certain allopatric populations show a preference for their own morph. We present evidence that this preference may be mediated by color pattern specifically. These results contrast with an earlier study of a third transition zone, in which a mimetic shift was associated with reproductive isolation. Overall, our results suggest that the three known mimetic transition zones in R. imitator reflect a speciation continuum, which we have characterized at the geographic, phenotypic, behavioral, and genetic levels. We discuss possible explanations for variable progress toward speciation, suggesting that multifarious selection on both mimetic color pattern and body size may be responsible for generating reproductive isolation.

Prey from the eyes of predators: Color discriminability of aposematic and mimetic butterflies from an avian visual perspective

Predation exerts strong selection on mimetic butterfly wing color patterns, which also serve other functions such as sexual selection. Therefore, specific selection pressures may affect the sexes and signal components differentially. We tested three predictions about the evolution of mimetic resemblance by comparing wing coloration of aposematic butterflies and their Batesian mimics: (a) females gain greater mimetic advantage than males and therefore are better mimics, (b) due to intersexual genetic correlations, sexually monomorphic mimics are better mimics than female-limited mimics, and (c) mimetic resemblance is better on the dorsal wing surface that is visible to predators in flight. Using a physiological model of avian color vision, we quantified mimetic resemblance from predators' perspective, which showed that female butterflies were better mimics than males. Mimetic resemblance in female-limited mimics was comparable to that in sexually monomorphic mimics, suggesting that intersexual genetic correlations did not constrain adaptive response to selection for female-limited mimicry. Mimetic resemblance on the ventral wing surface was better than that on the dorsal wing surface, implying stronger natural and sexual selection on ventral and dorsal surfaces, respectively. These results suggest that mimetic resemblance in butterfly mimicry rings has evolved under various selective pressures acting in a sex- and wing surface-specific manner.

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.

Warning signals are seductive: relative contributions of color and pattern to predator avoidance and mate attraction in Heliconius butterflies

Visual signaling in animals can serve many uses, including predator deterrence and mate attraction. In many cases, signals used to advertise unprofitability to predators are also used for intraspecific communication. Although aposematism and mate choice are significant forces driving the evolution of many animal phenotypes, the interplay between relevant visual signals remains little explored. Here, we address this question in the aposematic passion-vine butterfly Heliconius erato by using color- and pattern-manipulated models to test the contributions of different visual features to both mate choice and warning coloration. We found that the relative effectiveness of a model at escaping predation was correlated with its effectiveness at inducing mating behavior, and in both cases wing color was more predictive of presumptive fitness benefits than wing pattern. Overall, however, a combination of the natural (local) color and pattern was most successful for both predator deterrence and mate attraction. By exploring the relative contributions of color versus pattern composition in predation and mate preference studies, we have shown how both natural and sexual selection may work in parallel to drive the evolution of specific animal color patterns.

Phylogeography of Heliconius cydno and its closest relatives: disentangling their origin and diversification

The origins of the extraordinary diversity within the Neotropics have long fascinated biologists and naturalists. Yet, the underlying factors that have given rise to this diversity remain controversial. To test the relative importance of Quaternary climatic change and Neogene tectonic and paleogeographic reorganizations in the generation of biodiversity, we examine intraspecific variation across the Heliconius cydno radiation and compare this variation to that within the closely related Heliconius melpomene and Heliconius timareta radiations. Our data, which consist of both mtDNA and genome-scan data from nearly 2250 amplified fragment length polymorphism (AFLP) loci, reveal a complex history of differentiation and admixture at different geographic scales. Both mtDNA and AFLP phylogenies suggest that H. timareta and H. cydno are probably geographic extremes of the same radiation that probably diverged from H. melpomene prior to the Pliocene-Pleistocene boundary, consistent with hypotheses of diversification that rely on geological events in the Pliocene. The mtDNA suggests that this radiation originated in Central America or the northwestern region of South America, with a subsequent colonization of the eastern and western slopes of the Andes. Our genome-scan data indicate significant admixture among sympatric H. cydno/H. timareta and H. melpomene populations across the extensive geographic ranges of the two radiations. Within H. cydno, both mtDNA and AFLP data indicate significant population structure at local scales, with strong genetic differences even among adjacent H. cydno colour pattern races. These genetic patterns highlight the importance of past geoclimatic events, intraspecific gene flow, and local population differentiation in the origin and establishment of new adaptive forms.

Divergent warning patterns contribute to assortative mating between incipient Heliconius species

Theoretical models suggest that traits under divergent ecological selection, which also contribute to assortative mating, will facilitate speciation with gene flow. Evidence for these so-called "magic traits" now exists across a range of taxa. However, their importance during speciation will depend on the extent to which they contribute to reproductive isolation. Addressing this requires experiments to determine the exact cues involved as well as estimates of assortative mating in the wild. Heliconius butterflies are well known for their diversity of bright warning color patterns, and their amenability to experimental manipulation has provided an excellent opportunity to test their role in reproductive isolation. Here, we reveal that divergent color patterns contribute to mate recognition between the incipient species Heliconius himera and H. erato, a taxon pair for which assortative mating by color pattern has been demonstrated among wild individuals: First, we demonstrate that males are more likely to attempt to mate conspecific females; second, we show that males are more likely to approach pinned females that share their own warning pattern. These data are valuable as these taxa likely represent the early stages of speciation, but unusually also allow comparisons with rates of interbreeding between divergent ecologically relevant phenotypes measured in the wild.

Sexual Size Dimorphism in the Color Pattern Elements of Two Mimetic Heliconius Butterflies

Sexual dichromatism and sexual dimorphism of body size are reasonably well studied in butterflies. Sexual size dimorphism of color pattern elements, however, is much less explored. The object of this study is Heliconius, a genus of butterflies well known for the coevolution between mate color preferences and mimicry. Given the sexual role of wing coloration, we investigated the existence of sexual size dimorphism in the wing color elements of a mimetic pair-Heliconius erato phyllis Fabricius and Heliconius besckei Ménétriés-and analyzed the allometric patterns of these traits. Correlation between size of elements in the dorsal and ventral wing surfaces were also estimated. In both species, three out of four elements were larger in males, but the non-dimorphic element was not the same. With regard to the allometric patterns, our most important finding was that smaller males of one species have proportionally larger yellow bars. This is the first study specifically concerning quantitative sexual dimorphism in the coloration of this well-known genus of butterflies and it opens new prospects to investigate sex-related natural selection and sexual selection of color pattern elements.

Incomplete behavioural isolation and asymmetric female preference in darter sister species (Percidae: Etheostoma)

In this study, female preferences and behavioural isolation were estimated in a pair of allopatric sister species, Etheostoma duryi and Etheostoma flavum. Dichotomous mate preference trials were conducted to determine whether females prefer to associate with conspecific over heterospecific males and free-spawning assays were conducted to determine whether those preferences translated into behavioural isolation. Dichotomous mate choice trials revealed asymmetric female preference, as female E. flavum preferred conspecific males, whereas female E. duryi showed no preference. Free-spawning assays indicated that behavioural isolation remains incomplete between E. duryi and E. flavum (IB = 0·19). In addition to female mating preferences, male behaviour also appeared to influence mating outcomes as male E. flavum consistently courted conspecific females more often in free-spawning assays whereas male E. duryi did not. The data therefore suggest that despite marked divergence in male nuptial colour, divergence in female preferences between these species may not be sufficient to maintain species boundaries upon secondary contact. These results contrast with similar work in a sympatric darter species pair and may be explained by considering the contributions of reinforcement and differences in colour pattern as well as colour value.

Male courtship preferences demonstrate discrimination against allopatric colour morphs in a cichlid fish

Whether premating isolation is achieved by male-specific, female-specific or sex-independent assortative preferences often depends on the underlying evolutionary processes. Here we test mate preferences of males presented with females of different allopatric colour variants of the cichlid fish Tropheus sp., a Lake Tanganyika endemic with rich geographical colour pattern variation, in which the strength of sexual isolation varies between populations. We conducted two-way mate choice experiments to compare behaviour of males of a red-bodied morph (population Moliro) towards females from their own population with behaviour towards females from four allopatric populations at different stages of phylogenetic and phenotypic divergence. Males courted same-population females significantly more intensely than females of other populations, and reduced their heteromorphic courtship efforts both with increasing genetic and increasing phenotypic distinctness of the females. In particular, females of a closely related red-bodied population received significantly more courtship than either genetically distinct, similarly coloured females (‘Kirschfleck’ morph) or genetically related, differently coloured females (‘yellow-blotch’ morph), both of which were courted similarly. Genetically and phenotypically distinct females (Tropheus polli) were not courted at all. Consistent with previous female-choice experiments, female courtship activity also decreased with increasing genetic distance from the males’ population. Given successful experimental and natural introgression between colour morphs and the pervasive allopatry of related variants, we consider it unlikely that assortative preferences of both sexes were driven by direct selection during periods of secondary contact or, in turn, drove colour pattern differentiation in allopatry. Rather, we suggest that sexual isolation evolved as by-product of allopatric divergence.

Disruptive ecological selection on a mating cue

Adaptation to divergent ecological niches can result in speciation. Traits subject to disruptive selection that also contribute to non-random mating will facilitate speciation with gene flow. Such 'magic' or 'multiple-effect' traits may be widespread and important for generating biodiversity, but strong empirical evidence is still lacking. Although there is evidence that putative ecological traits are indeed involved in assortative mating, evidence that these same traits are under divergent selection is considerably weaker. Heliconius butterfly wing patterns are subject to positive frequency-dependent selection by predators, owing to aposematism and Müllerian mimicry, and divergent colour patterns are used by closely related species to recognize potential mates. The amenability of colour patterns to experimental manipulation, independent of other traits, presents an excellent opportunity to test their role during speciation. We conducted field experiments with artificial butterflies, designed to match natural butterflies with respect to avian vision. These were complemented with enclosure trials with live birds and real butterflies. Our experiments showed that hybrid colour-pattern phenotypes are attacked more frequently than parental forms. For the first time, we demonstrate disruptive ecological selection on a trait that also acts as a mating cue.

Biased learning affects mate choice in a butterfly

Early acquisition of mate preferences or mate-preference learning is associated with signal diversity and speciation in a wide variety of animal species. However, the diversity of mechanisms of mate-preference learning across taxa remains poorly understood. Using the butterfly Bicyclus anynana we uncover a mechanism that can lead to directional sexual selection via mate-preference learning: a bias in learning enhanced ornamentation, which is independent of preexisting mating biases. Naïve females mated preferentially with wild-type males over males with enhanced wing ornamentation, but females briefly exposed to enhanced males mated significantly more often with enhanced males. In contrast, females exposed to males with reduced wing ornamentation did not learn to prefer drab males. Thus, we observe both a learned change of a preexisting mating bias, and a bias in ability to learn enhanced male ornaments over reduced ornaments. Our findings demonstrate that females are able to change their preferences in response to a single social event, and suggest a role for biased learning in the evolution of visual sexual ornamentation.

Mate choice and the genetic basis for colour variation in a polymorphic dart frog: inferences from a wild pedigree

Understanding how reproductive barriers evolve during speciation remains an important question in evolution. Divergence in mating preferences may be a common first step in this process. The striking colour pattern diversity of strawberry dart frog (Dendrobates pumilio) populations has likely been shaped by sexual selection. Previous laboratory studies have shown that females attend to male coloration and prefer to court with males of their own colour, suggesting that divergent morphs may be reproductively isolated. To test this hypothesis, we used molecular data to estimate pedigree relationships from a polymorphic population. Whereas in the laboratory both red and yellow females preferred to court with males of their own phenotype, our pedigree shows a pattern of assortative mating only for red females. In the wild, yellow females appear to be less choosy about their mates, perhaps because they incur higher costs associated with searching than females of the more common red phenotype. We also used our pedigree to investigate the genetic basis for colour-pattern variation. The phenotype frequencies we observed were consistent with those expected if dorsal background coloration is controlled by a single locus, with complete dominance of red over yellow. Our results not only help clarify the role of sexual selection in reducing gene flow, but also shed light on the mechanisms underlying colour-pattern variation among sympatric colour morphs. The difference we observed between mating preferences measured under laboratory conditions and the pattern of mate choice observed in the wild highlight the importance of field studies for understanding behavioural reproductive isolation.

Genomic islands of divergence in hybridizing Heliconius butterflies identified by large-scale targeted sequencing

Heliconius butterflies represent a recent radiation of species, in which wing pattern divergence has been implicated in speciation. Several loci that control wing pattern phenotypes have been mapped and two were identified through sequencing. These same gene regions play a role in adaptation across the whole Heliconius radiation. Previous studies of population genetic patterns at these regions have sequenced small amplicons. Here, we use targeted next-generation sequence capture to survey patterns of divergence across these entire regions in divergent geographical races and species of Heliconius. This technique was successful both within and between species for obtaining high coverage of almost all coding regions and sufficient coverage of non-coding regions to perform population genetic analyses. We find major peaks of elevated population differentiation between races across hybrid zones, which indicate regions under strong divergent selection. These 'islands' of divergence appear to be more extensive between closely related species, but there is less clear evidence for such islands between more distantly related species at two further points along the 'speciation continuum'. We also sequence fosmid clones across these regions in different Heliconius melpomene races. We find no major structural rearrangements but many relatively large (greater than 1 kb) insertion/deletion events (including gain/loss of transposable elements) that are variable between races.

Introgressive hybridization between color morphs in a population of cichlid fishes twelve years after human-induced secondary admixis

In the extremely species-rich haplochromine cichlid fishes of the East African Great Lakes, prezygotic isolation between closely related species is often maintained by color-assortative mating. In 1998, local fisherman working for the ornamental fish trade released different color morphs of the cichlid genus Tropheus into a small harbor basin in the southern part of Lake Tanganyika. This artificial amalgamation of color morphs provides a unique possibility to study mating patterns in cichlids in a natural environment over time. In a precursor study, we analyzed genotypes and phenotypes of almost 500 individuals sampled between 1999 and 2001 and uncovered a marked degree of color-assortative mating, which depended on the level of color pattern dissimilarity between morphs. Twelve years after introduction of nonindigenous morphs, we again sampled Tropheus individuals from the harbor basin and an adjacent, originally pure population and analyzed phenotypes (coloration) and genotypes (mitochondrial control region and 9 microsatellite loci) to assess the current status of the admixed population. Principal component analyses of color score data and population assignment tests demonstrate an increasing level of introgressive hybridization between morphs but also some ongoing color-assortative mating within morphs. The observed mating pattern might have been influenced by fluctuating environmental conditions such as periodic algal blooms or increased sedimentation causing turbid conditions in an otherwise clear lake.

Genes versus phenotypes in the study of speciation

Despite persistent debate on the nature of species, the widespread adoption of Mayr's biological species concept has led to a heavy emphasis on the importance of reproductive isolation to the speciation process. Equating the origin of species with the evolution of reproductive isolation has become common practice in the study of speciation, coincident with an increasing focus on elucidating the specific genetic changes (i.e.-speciation genes) underlying intrinsic reproductive barriers between species. In contrast, some have recognized that reproductive isolation is usually a signature effect rather than a primary cause of speciation. Here we describe a research paradigm that shifts emphasis from effects to causes in order to resolve this apparent contradiction and galvanize the study of speciation. We identify major elements necessary for a balanced and comprehensive investigation of the origin of species and place the study of so-called "speciation genes" into its appropriate context. We emphasize the importance of characterizing diverging phenotypes, identifying relevant evolutionary forces acting on those phenotypes and their role in the causal origins of reduced gene flow between incipient species, and the nature of the genetic and phenotypic boundaries that results from such processes. This approach has the potential to unify the field of speciation research, by allowing us to make better "historical" predictions about the fate of diverging populations regardless of taxon.

Characterisation and expression of microRNAs in developing wings of the neotropical butterfly Heliconius melpomene

BACKGROUND

Heliconius butterflies are an excellent system for studies of adaptive convergent and divergent phenotypic traits. Wing colour patterns are used as signals to both predators and potential mates and are inherited in a Mendelian manner. The underlying genetic mechanisms of pattern formation have been studied for many years and shed light on broad issues, such as the repeatability of evolution. In Heliconius melpomene, the yellow hindwing bar is controlled by the HmYb locus. MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression that have key roles in many biological processes, including development. miRNAs could act as regulators of genes involved in wing development, patterning and pigmentation. For this reason we characterised miRNAs in developing butterfly wings and examined differences in their expression between colour pattern races.

RESULTS

We sequenced small RNA libraries from two colour pattern races and detected 142 Heliconius miRNAs with homology to others found in miRBase. Several highly abundant miRNAs were differentially represented in the libraries between colour pattern races. These candidates were tested further using Northern blots, showing that differences in expression were primarily due to developmental stage rather than colour pattern. Assembly of sequenced reads to the HmYb region identified hme-miR-193 and hme-miR-2788; located 2380 bp apart in an intergenic region. These two miRNAs are expressed in wings and show an upregulation between 24 and 72 hours post-pupation, indicating a potential role in butterfly wing development. A search for miRNAs in all available H. melpomene BAC sequences (~2.5 Mb) did not reveal any other miRNAs and no novel miRNAs were predicted.

CONCLUSIONS

Here we describe the first butterfly miRNAs and characterise their expression in developing wings. Some show differences in expression across developing pupal stages and may have important functions in butterfly wing development. Two miRNAs were located in the HmYb region and were expressed in developing pupal wings. Future work will examine the expression of these miRNAs in different colour pattern races and identify miRNA targets among wing patterning genes.

Mate preference across the speciation continuum in a clade of mimetic butterflies

Premating behavioral isolation is increasingly recognized as an important part of ecological speciation, where divergent natural selection causes the evolution of reproductive barriers. A number of studies have now demonstrated that traits under divergent natural selection also affect mate preferences. However, studies of single species pairs only capture a snapshot of the speciation process, making it difficult to assess the role of mate preferences throughout the entire process. Heliconius butterflies are well known for their brightly colored mimetic warning patterns, and previous studies have shown that these patterns are also used as mate recognition cues. Here, we present mate preference data for four pairs of sister taxa, representing different stages of divergence, which together allow us to compare diverging mate preferences across the continuum of Heliconius speciation. Using a novel Bayesian approach, our results support a model of ecological speciation in which strong premating isolation arises early, but continues to increase throughout the continuum from polymorphic populations through to "good," sympatric ecologically divergent species.

A guide to the genomics of ecological speciation in natural animal populations

Interest in ecological speciation is growing, as evidence accumulates showing that natural selection can lead to rapid divergence between subpopulations. However, whether and how ecological divergence can lead to the buildup of reproductive isolation remains under debate. What is the relative importance of natural selection vs. neutral processes? How does adaptation generate reproductive isolation? Can ecological speciation occur despite homogenizing gene flow? These questions can be addressed using genomic approaches, and with the rapid development of genomic technology, will become more answerable in studies of wild populations than ever before. In this article, we identify open questions in ecological speciation theory and suggest useful genomic methods for addressing these questions in natural animal populations. We aim to provide a practical guide for ecologists interested in incorporating genomic methods into their research programs. An increased integration between ecological research and genomics has the potential to shed novel light on the origin of species.