Ecological and genetic associations across a Heliconius hybrid zone

Differences in habitat use can bridge early and late stages of speciation by initiating assortative mating. Heliconius colour pattern races might select habitats over which each pattern confers a relative fitness advantage because signal efficacy of wing patterns can vary by environment. Thus habitat preferences could serve to promote the evolution of mimetic colour patterns for mate choice. Here I compare colour pattern genotype and phenotype frequencies to environmental variation across the H. erato hydara x H. erato erato hybrid zone in French Guiana to determine whether races exhibit habitat preferences. I found that genotype and phenotype frequencies correspond to differences in land cover moreso than to other environmental factors. Temporal shifts in colour pattern genotypes, phenotypes and land cover also were associated at individual sample sites, which further suggests that H. erato races differ in habitat use and that habitat preferences may promote speciation among Heliconius butterflies.     

Interspecific sexual attraction because of convergence in warning colouration: is there a conflict between natural and sexual selection in mimetic species?

When species converge in their colour patterns because of mimicry, and those patterns are also used in mate recognition, there is a probability of conflicting selection pressures. Closely related species that mimic one another are particularly likely to face such confusion because of similarities in their courtship behaviour and ecology. We conducted experiments in greenhouse conditions to study interspecific attraction between two mimetic butterfly species, Heliconius erato and Heliconius melpomene. Both species spent considerable time approaching and courting females of the co-mimic species. Experiments using wing models demonstrated the importance of colour pattern in this interspecific attraction. Although males of H. melpomene were attracted to their co-mimics as much as to their own females, H. erato males were more efficient at distinguishing conspecifics, possibly using wing odours. Although preliminary, these results suggest that the use of additional cues may have evolved in H. erato to reduce the cost of convergence in visual signals with H. melpomene. Overall, our results showed that there might be a cost of mimetic convergence because of a reduction in the efficiency of species recognition. Such cost may contribute to explain the apparently stable diversity in Müllerian mimetic patterns in many tropical butterfly assemblages.     

A hybrid zone provides evidence for incipient ecological speciation in Heliconius butterflies

In Heliconius butterflies, it has been proposed that speciation occurs through a combination of divergence in ecological habitat preferences and mimetic colour patterns. Here we test this hypothesis by investigating a parapatric form of the widespread species Heliconius erato. Mendelian (colour patterns) and molecular genetic data permit us to address hypotheses about introgression and genetic differentiation between different populations. Combined analysis of colour pattern, microsatellite loci and mitochondrial DNA showed that Heliconius erato venus and Heliconius erato chestertonii form a bimodal hybrid zone implying partial reproductive isolation. In a sample of 121 individuals collected in sympatry, 25% were hybrids representing a significant deficit of heterozygotes compared to the Hardy-Weinberg expectation. Seven microsatellite loci, analysed for a subset of these individuals, showed marked differentiation between the parental taxa, and unambiguously identified two genotypic clusters concordant with our phenotypic classification of individuals. Mitochondrial DNA analysis showed H. erato venus as a monophyletic group well differentiated from H. erato chestertonii, implying a lack of historical introgression between the populations. Heliconius erato chestertonii is therefore an incipient species that maintains its integrity despite high levels of hybridization. Moreover, H. erato chestertonii is found at higher altitudes than other races of H. erato and has a distinct colour pattern and mimetic relationship. Hence, there are now two examples of parapatric incipient species related to H. erato, H. himera and H. erato chestertonii, both of which are associated with higher altitudes, more arid habitats and distinct mimetic relationships. This implies that parapatric habitat adaptation is a likely cause of speciation in this group.     

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.     

Patterns of pollen feeding and habitat preference among Heliconius species

Abstract 1. The ecological circumstances that precipitate speciation remain poorly understood. Here, a community of Heliconius butterflies in lowland Panama was studied to investigate patterns of pollen use, and more specifically the ecological changes associated with the recent divergence of Heliconius melpomene (Linnaeus) and H. cydno (Doubleday).
2. Considering the seven commonest Heliconius species in the community, 32 types of pollen or spore were encountered in pollen loads but only five pollen species were common. Systematic exploitation of pollen was therefore confined to a small proportion of the flowers visited.
3. Most of the variation in pollen load composition between individuals was explained by differences in collecting locality. The exception was Psiguria , which was used in all habitats by the melpomene / hecale clade far more than by the erato / sapho clade. This may suggest an ancestral switch within Heliconius towards increased reliance on Psiguria pollen.
4. Heliconius cydno and H. melpomene differed significantly in pollen load composition for three of the five most commonly collected pollen species. This is most probably explained by differences in habitat preference; H. melpomene and its co-mimic H. erato are found in open habitat while H. cydno and its co-mimic H. sapho are found in closed-canopy forest.
5. As melpomene and cydno are known to hybridise occasionally, such differences in adult microhabitat contribute to pre-mating isolation. Habitat divergence between H. cydno and H. melpomene , which is associated with changes in mimicry, must have played a role in their recent speciation.     

Phylogenetic codivergence supports coevolution of mimetic Heliconius butterflies

The unpalatable and warning-patterned butterflies Heliconius erato and Heliconius melpomene provide the best studied example of mutualistic Müllerian mimicry, thought-but rarely demonstrated-to promote coevolution. Some of the strongest available evidence for coevolution comes from phylogenetic codivergence, the parallel divergence of ecologically associated lineages. Early evolutionary reconstructions suggested codivergence between mimetic populations of H. erato and H. melpomene, and this was initially hailed as one of the most striking known cases of coevolution. However, subsequent molecular phylogenetic analyses found discrepancies in phylogenetic branching patterns and timing (topological and temporal incongruence) that argued against codivergence. We present the first explicit cophylogenetic test of codivergence between mimetic populations of H. erato and H. melpomene, and re-examine the timing of these radiations. We find statistically significant topological congruence between multilocus coalescent population phylogenies of H. erato and H. melpomene. Cophylogenetic historical reconstructions support repeated codivergence of mimetic populations, from the base of the sampled radiations. Pairwise distance correlation tests, based on our coalescent analyses plus recently published AFLP and wing colour pattern gene data, also suggest that the phylogenies of H. erato and H. melpomene show significant topological congruence. Divergence time estimates, based on a Bayesian coalescent model, suggest that the evolutionary radiations of H. erato and H. melpomene occurred over the same time period, and are compatible with a series of temporally congruent codivergence events. Our results suggest that differences in within-species genetic divergence are the result of a greater overall effective population size for H. erato relative to H. melpomene and do not imply incongruence in the timing of their phylogenetic radiations. Repeated codivergence between Müllerian co-mimics, predicted to exert mutual selection pressures, strongly suggests coevolution. Our results therefore support a history of reciprocal coevolution between Müllerian co-mimics characterised by phylogenetic codivergence and parallel phenotypic change.     

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.     

Review. Hybrid trait speciation and Heliconius butterflies

Homoploid hybrid speciation (HHS) is the establishment of a novel species through introgressive hybridization without a change in chromosome number. We discuss different routes by which this might occur and propose a novel term, 'hybrid trait speciation', which combines the idea that hybridization can generate adaptive novelty with the 'magic trait' model of ecological speciation. Heliconius butterflies contain many putative examples of hybrid colour patterns, but only recently has the HHS hypothesis been tested explicitly in this group. Molecular data has shown evidence for gene flow between many distinct species. Furthermore, the colour pattern of Heliconius heurippa can be recreated in laboratory crosses between Heliconius melpomene and Heliconius cydno and, crucially, plays a role in assortative mating between the three species. Nonetheless, although the genome of H. heurippa shows evidence for hybridization, it is not a mosaic of the two parental species. Instead, ongoing hybridization has likely blurred any signal of the original speciation event. We argue that where hybridization leads to novel adaptive traits that also cause reproductive isolation, it is likely to trigger speciation.     

Pervasive genetic associations between traits causing reproductive isolation in Heliconius butterflies

Ecological speciation proceeds through the accumulation of divergent traits that contribute to reproductive isolation, but in the face of gene flow traits that characterize incipient species may become disassociated through recombination. Heliconius butterflies are well known for bright mimetic warning patterns that are also used in mate recognition and cause both pre- and post-mating isolation between divergent taxa. Sympatric sister taxa representing the final stages of speciation, such as Heliconius cydno and Heliconius melpomene, also differ in ecology and hybrid fertility. We examine mate preference and sterility among offspring of crosses between these species and demonstrate the clustering of Mendelian colour pattern loci and behavioural loci that contribute to reproductive isolation. In particular, male preference for red patterns is associated with the locus responsible for the red forewing band. Two further colour pattern loci are associated, respectively, with female mating outcome and hybrid sterility. This genetic architecture in which ‘speciation genes’ are clustered in the genome can facilitate two controversial models of speciation, namely divergence in the face of gene flow and hybrid speciation.     

Adaptive introgression across species boundaries in Heliconius butterflies

It is widely documented that hybridisation occurs between many closely related species, but the importance of introgression in adaptive evolution remains unclear, especially in animals. Here, we have examined the role of introgressive hybridisation in transferring adaptations between mimetic Heliconius butterflies, taking advantage of the recent identification of a gene regulating red wing patterns in this genus. By sequencing regions both linked and unlinked to the red colour locus, we found a region that displays an almost perfect genotype by phenotype association across four species, H. melpomene, H. cydno, H. timareta, and H. heurippa. This particular segment is located 70 kb downstream of the red colour specification gene optix, and coalescent analysis indicates repeated introgression of adaptive alleles from H. melpomene into the H. cydno species clade. Our analytical methods complement recent genome scale data for the same region and suggest adaptive introgression has a crucial role in generating adaptive wing colour diversity in this group of butterflies.     

Evolution of the wave: aerodynamic and aposematic functions of butterfly wing motion

Many unpalatable butterfly species use coloration to signal their distastefulness to birds, but motion cues may also be crucial to ward off predatory attacks. In previous research, captive passion-vine butterflies Heliconius mimetic in colour pattern were also mimetic in motion. Here, I investigate whether wing motion changes with the flight demands of different behaviours. If birds select for wing motion as a warning signal, aposematic butterflies should maintain wing motion independently of behavioural context. Members of one mimicry group (Heliconius cydno and Heliconius sapho) beat their wings more slowly and their wing strokes were more asymmetric than their sister-species (Heliconius melpomene and Heliconius erato, respectively), which were members of another mimicry group having a quick and steady wing motion. Within mimicry groups, wing beat frequency declined as its role in generating lift also declined in different behavioural contexts. In contrast, asymmetry of the stroke was not associated with wing beat frequency or behavioural context-strong indication that birds process and store the Fourier motion energy of butterfly wings. Although direct evidence that birds respond to subtle differences in butterfly wing motion is lacking, birds appear to generalize a motion pattern as much as they encounter members of a mimicry group in different behavioural contexts.     

ASSORTATIVE MATING PREFERENCES AMONG HYBRIDS OFFERS A ROUTE TO HYBRID SPECIATION

Homoploid speciation generates species without a change in chromosome number via introgressive hybridization and has been considered rare in animals. Heliconius butterflies exhibit bright aposematic color patterns that also act as cues in assortative mating. Heliconius heurippa has a color pattern that can be recreated by introgression of the H. melpomene red band into an H. cydno genetic background. Wild H. heurippa males show assortative mating based on color pattern and we here investigate the origin of this preference by studying first-generation backcross hybrids between H. melpomene and H. cydno that resemble H. heurippa . These hybrids show assortative mating preferences, showing a strong preference for their own color pattern over that of either parental species. This is consistent with a genetic basis to wing pattern preference and implies, first, that assortative mating preferences would facilitate the initial establishment of a homozygous hybrid color pattern by increasing the likelihood that early generation hybrids mate among themselves. Second, once established such a lineage would inherit assortative mating preferences that would lead to partial reproductive isolation from parental lineages.     

Shared and divergent expression domains on mimetic Heliconius wings

SUMMARY Heliconius butterfly wing patterns show repeated convergence between species and have adaptive value in mimicry and mate choice, offering an opportunity to connect adaptive changes in phenotype with their underlying genotypes. Here we study forewing ommochrome pigmentation in Heliconius melpomene . We clone two new ommochrome pathway genes for the Lepidoptera, karmoisin and kynurenine formamidase ( kf  ), and analyze the expression patterns of all known ommochrome genes across pupal wing development. In combination with published work, this generates the first comparative gene expression data for the co-mimics Heliconius erato and H. melpomene . In both species cinnabar expression correlates with the forewing band, but the expression pattern of vermillion differs significantly between the mimics. This demonstrates that both shared and divergent expression patterns are associated with mimetic phenotypes between Heliconius species. Two genes not studied in H. erato, scarlet and possibly kf , also show enhanced expression in the forewing band of H. melpomene , implying co-ordinated upregulation of several members of this biosynthetic pathway during pattern formation.     

A conserved supergene locus controls colour pattern diversity in Heliconius butterflies

We studied whether similar developmental genetic mechanisms are involved in both convergent and divergent evolution. Mimetic insects are known for their diversity of patterns as well as their remarkable evolutionary convergence, and they have played an important role in controversies over the respective roles of selection and constraints in adaptive evolution. Here we contrast three butterfly species, all classic examples of Müllerian mimicry. We used a genetic linkage map to show that a locus, Yb, which controls the presence of a yellow band in geographic races of Heliconius melpomene, maps precisely to the same location as the locus Cr, which has very similar phenotypic effects in its co-mimic H. erato. Furthermore, the same genomic location acts as a "supergene", determining multiple sympatric morphs in a third species, H. numata. H. numata is a species with a very different phenotypic appearance, whose many forms mimic different unrelated ithomiine butterflies in the genus Melinaea. Other unlinked colour pattern loci map to a homologous linkage group in the co-mimics H. melpomene and H. erato, but they are not involved in mimetic polymorphism in H. numata. Hence, a single region from the multilocus colour pattern architecture of H. melpomene and H. erato appears to have gained control of the entire wing-pattern variability in H. numata, presumably as a result of selection for mimetic "supergene" polymorphism without intermediates. Although we cannot at this stage confirm the homology of the loci segregating in the three species, our results imply that a conserved yet relatively unconstrained mechanism underlying pattern switching can affect mimicry in radically different ways. We also show that adaptive evolution, both convergent and diversifying, can occur by the repeated involvement of the same genomic regions.     

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.     

Genetic analysis of a wild-caught hybrid between non-sister Heliconius butterfly species

Interspecific hybridization occurs regularly in wild Heliconius butterflies, although hybrid individuals are usually very rare. However, hybridization generally occurs only between the most closely related species. We report a rare naturally occurring hybrid between non-sister species and carry out the first genetic analysis of such distant hybridization. Mitochondrial and nuclear genes indicate that the specimen is an F1 hybrid between a female Heliconius ethilla and a male Heliconius melpomene, originating from a group of 13 species estimated to have diverged over 2.5 Myr ago. The presence of such distant natural hybrids, together with evidence for backcrossing, suggests that gene flow across species boundaries can take place long after speciation. Adaptive genes such as those involved in wing coloration could thus be widely shared among members of this highly mimetic genus.     

Host choice promotes reproductive isolation between host races of the larch budmoth Zeiraphera diniana

The chances for sympatric speciation are improved if ecological divergence leads to assortative mating as a by-product. This effect is known in parasites that find mates using host cues, but studies of larch- and pine-feeding races of the larch budmoth (Zeiraphera diniana, Lepidoptera: Tortricidae) suggest it may also occur when mate attraction is via sex pheromones that are independent of habitat. We have previously shown that females releasing pheromones on or near their own host attract more males of their own race than if placed on the alternative host. This host effect would enhance assortative mating provided adults preferentially alight on their native hosts. Here we investigate alighting preferences in natural mixed forest using a novel likelihood analysis of genotypic clusters based on three semidiagnostic allozyme loci. Both larch and pine females show a realized alighting preference for their own host of 86%. The equivalent preferences of males were 79% for the larch race and 85% for the pine race. These preferences are also detectable in small-scale laboratory experiments, where alighting preferences of larch and pine races towards their own hosts were, respectively, 67 and 66% in females and 69 and 63% in males. Pure larch race moths reared in the laboratory had alighting choice similar to moths from natural populations, while hybrids were intermediate, showing that alighting preferences were heritable and approximately additive. The field estimates of alighting preference, coupled with earlier work on mate choice, yield an estimated rate of natural hybridization between sympatric host races of 2.2-3.8% per generation. Divergent alighting choice enhances pheromone-mediated assortative mating today, and is likely to have been an important cause of assortative mating during initial divergence in host use. Because resources are normally 'coarse-grained' in space and time, assortative mating due to ecological divergence may be a more important catalyst of sympatric speciation than generally realized.     

Reinforcement of mate preference among hybridizing Heliconius butterflies

Recent models of mate preference evolution suggest that direct selection on alleles at preference loci and correlated evolution of preference with locally adapted mating cues are more likely to drive the evolution of assortative mate preference than reinforcement. Mate preference evolution in mimetic Heliconius butterflies has been attributed to all three forms of selection, but here we show that reinforcement has been critical. By examining geographical variation in assortative mating and male mate preference among seven populations of three hybridizing Heliconius species from Costa Rica, we found pronounced character displacement of preference such that sexual isolation was enhanced in areas of interspecific contact. Of the different explanations for the evolution of assortative mate preference, only reinforcement is dependent on interspecific contact in this system. Thus, the observed pattern of reproductive character displacement of mate preference is best explained as a product of indirect selection generated by natural selection against nonmimetic hybrids.     

What can hybrid zones tell us about speciation? The case of Heliconius erato and H. himera (Lepidoptera: Nymphalidae)

To understand speciation we need to study the genetics and ecology of intermediate cases where interspecific hybridization still occurs. Two closely related species of Heliconius butterflies meet this criterion: Heliconius himera is endemic to dry forest and thorn scrub in southern Ecuador and northern Peru, while its sister species, H. erato , is ubiquitous in wet forest throughout south and central America. In three known zones of contact, the two species remain distinct, while hybrids are found at low frequency. Collections in southern Ecuador show that the contact zone is about 5 km wide, half the width of the narrowest clines between colour pattern races of H. erato. The narrowness of this dine argues that very strong selection (s ≅ 1) is maintaining the parapatric distributions of these two species. The zone is closely related with a habitat transition from wet to dry forest, which suggests that the narrow zone of parapatry is maintained primarily by ecological adaptation. Selection on colour pattern loci, assortative mating and hybrid inviability may also be important. The genetics of hybrids between the two species shows that the major gene control of pattern elements is similar to that found in previous studies of H. erato races, and some of the loci are homologous. This suggests that similar genetic processes are involved in the morphological divergence of species and races. Evidence from related Heliconius supports a hypothesis that ecological adaptation is the driving force for speciation in the group.     

Mimicry: developmental genes that contribute to speciation

Despite renewed interest in the role of natural selection as a catalyst for the origin of species, the developmental and genetic basis of speciation remains poorly understood. Here we describe the genetics of Müllerian mimicry in Heliconius cydno and H. melpomene (Lepidoptera: Nymphalidae), sister species that recently diverged to mimic other Heliconius. This mimetic shift was a key step in their speciation, leading to pre- and postmating isolation. We identify 10 autosomal loci, half of which have major effects. At least eight appear to be homologous with genes known to control pattern differences within each species. Dominance has evolved under the influence of identifiable "modifier" loci rather than being a fixed characteristic of each locus. Epistasis is found at many levels: phenotypic interaction between specific pairs of genes, developmental canalization due to polygenic modifiers so that patterns are less sharply defined in hybrids, and overall fitness through ecological selection against nonmimetic hybrid genotypes. Most of the loci are clustered into two genomic regions or "supergenes," suggesting color pattern evolution is constrained by preexisting linked elements that may have arisen via tandem duplication rather than having been assembled by natural selection. Linkage, modifiers, and epistasis affect the strength of mimicry as a barrier to gene flow between these naturally hybridizing species and may permit introgression in genomic regions unlinked to those under disruptive selection. Müllerian mimics in Heliconius use different genetic architectures to achieve the same mimetic patterns, implying few developmental constraints. Therefore, although developmental and genomic constraints undoubtedly influence the evolutionary process, their effects are probably not strong in comparison with natural selection.