Hybrid incompatibility is consistent with a hybrid origin of Heliconius heurippa Hewitson from its close relatives, Heliconius cydno Doubleday and Heliconius melpomene Linnaeus

Abstract Shared ancestral variation and introgression complicates the reconstruction of phylogenetic relationships among closely related taxa. Here we use overall genomic compatibility as an alternative estimate of species relationships in a group where divergence is rapid and genetic exchange is common. Heliconius heurippa, a butterfly species endemic to Colombia, has a colour pattern genetically intermediate between H. cydno and H. melpomene: its hindwing is nearly indistinguishable from that of H. melpomene and its forewing band is an intermediate phenotype between both species. This observation has lead to the suggestion that the pattern of H. heurippa arose through hybridization. We present a genetic analysis of hybrid compatibility in crosses between the three taxa. Heliconius heurippa x H. cydno and female H. melpomene x male H. heurippa yield fertile and viable F1 hybrids, but male H. melpomene x female H. heurippa crosses yield sterile F1 females. In contrast, Haldane's rule has previously been detected between H. melpomene and H cydno in both directions. Therefore, H. heurippa is most closely related to H. cydno, with some evidence for introgression of genes from H. melpomene. The results are compatible with the hypothesis of a hybrid origin for H. heurippa. In addition, backcrosses using F1 hybrid males provide evidence for a large Z(X)-chromosome effect on sterility and for recessive autosomal sterility factors as predicted by Dominance Theory.     

CASE STUDIES AND MATHEMATICAL MODELS OF ECOLOGICAL SPECIATION. 4. HYBRID SPECIATION IN BUTTERFLIES IN A JUNGLE

We build a spatial individual-based multilocus model of homoploid hybrid speciation tailored for a tentative case of hybrid origin of Heliconius heurippa from H. melpomene and H. cydno in South America. Our model attempts to account for empirical patterns and data on genetic incompatibility, mating preferences and selection by predation (both based on coloration patterns), habitat preference, and local adaptation for all three Heliconius species. Using this model, we study the likelihood of recombinational speciation and identify the effects of various ecological and genetic parameters on the dynamics, patterns, and consequences of hybrid ecological speciation. Overall, our model supports the possibility of hybrid origin of H. heurippa under certain conditions. The most plausible scenario would include hybridization between H. melpomene and H. cydno in an area geographically isolated from the rest of both parental species with subsequent long-lasting geographic isolation of the new hybrid species, followed by changes in the species ranges, the secondary contact, and disappearance of H. melpomene -type ecomorph in the hybrid species. However, much more work (both empirical and theoretical) is necessary to be able to make more definite conclusions on the importance of homoploid hybrid speciation in animals.     

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.     

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.     

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.     

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.     

Multilocus analyses of admixture and introgression among hybridizing Heliconius butterflies

Introgressive hybridization is an important evolutionary process and new analytical methods provide substantial power to detect and quantify it. In this study we use variation in the frequency of 657 AFLP fragments and DNA sequence variation from 15 genes to measure the extent of admixture and the direction of interspecific gene flow among three Heliconius butterfly species that diverged recently as a result of natural selection for Miillerian mimicry, and which continue to hybridize. Bayesian clustering based on AFLP genotypes correctly delineated the three species and identified four H. cydno, three H. pachinus, and three H. melpomene individuals that were of mixed ancestry. Gene genealogies revealed substantial shared DNA sequence variation among all three species and coalescent simulations based on the Isolation with Migration (IM) model pointed to interspecific gene flow as its cause. The IM simulations further indicated that interspecific gene flow was significantly asymmetrical, with greater gene flow from H. pachinus into H. cydno (2Nm = 4.326) than the reverse (2Nm = 0.502), and unidirectional gene flow from H. cydno and H. pachinus into H. melpomene (2Nm = 0.294 and 0.252, respectively). These asymmetries are in the directions expected based on the genetics of wing patterning and the probability that hybrids of various phenotypes will survive and reproduce in different mimetic environments. This empirical demonstration of extensive interspecific gene flow is in contrast to a previous study which found little evidence of gene flow between another pair of hybridizing Heliconius species, H. himera and H. erato, and it highlights the critical role of natural selection in maintaining species diversity. Furthermore, these results lend support to the hypotheses that phenotypic diversification in the genus Heliconius has been fueled by introgressive hybridization and that reinforcement has driven the evolution of assortative mate preferences.     

Comparative population genetics of a mimicry locus among hybridizing Heliconius butterfly species

The comimetic Heliconius butterfly species pair, H. erato and H. melpomene, appear to use a conserved Mendelian switch locus to generate their matching red wing patterns. Here we investigate whether H. cydno and H. pachinus, species closely related to H. melpomene, use this same switch locus to generate their highly divergent red and brown color pattern elements. Using an F2 intercross between H. cydno and H. pachinus, we first map the genomic positions of two novel red/brown wing pattern elements; the G locus, which controls the presence of red vs brown at the base of the ventral wings, and the Br locus, which controls the presence vs absence of a brown oval pattern on the ventral hind wing. The results reveal that the G locus is tightly linked to markers in the genomic interval that controls red wing pattern elements of H. erato and H. melpomene. Br is on the same linkage group but approximately 26 cM away. Next, we analyze fine-scale patterns of genetic differentiation and linkage disequilibrium throughout the G locus candidate interval in H. cydno, H. pachinus and H. melpomene, and find evidence for elevated differentiation between H. cydno and H. pachinus, but no localized signature of association. Overall, these results indicate that the G locus maps to the same interval as the locus controlling red patterning in H. melpomene and H. erato. This, in turn, suggests that the genes controlling red pattern elements may be homologous across Heliconius, supporting the hypothesis that Heliconius butterflies use a limited suite of conserved genetic switch loci to generate both convergent and divergent wing patterns.     

Mimicry and the evolution of premating isolation in Heliconius melpomene Linnaeus

Ecological divergence can cause speciation if adaptive traits have pleiotropic effects on mate choice. In Heliconius butterflies, mimetic patterns play a role in mate detection between sister species, as well as signalling to predators. Here we show that male butterflies from four recently diverged parapatric populations of Heliconius melpomene are more likely to approach and court their own colour patterns as compared with those of other races. A few exceptions, where males were more attracted to patterns other than their own, suggest that some mimetic patterns are sub-optimal in mate choice. Genotype frequencies in hybrid zones between races of H. melpomene suggest that mating is random, so reinforcement is unlikely to have played a role in intra-specific divergence. In summary, co-evolved divergence of colour pattern and mate preference occurs rapidly and is likely the first step in Heliconius speciation.     

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.     

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.     

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.     

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.     

Disruptive sexual selection against hybrids contributes to speciation between Heliconius cydno and Heliconius melpomene

Understanding the fate of hybrids in wild populations is fundamental to understanding speciation. Here we provide evidence for disruptive sexual selection against hybrids between Heliconius cydno and Heliconius melpomene. The two species are sympatric across most of Central and Andean South America, and coexist despite a low level of hybridization. No-choice mating experiments show strong assortative mating between the species. Hybrids mate readily with one another, but both sexes show a reduction in mating success of over 50% with the parental species. Mating preference is associated with a shift in the adult colour pattern, which is involved in predator defence through Müllerian mimicry, but also strongly affects male courtship probability. The hybrids, which lie outside the curve of protection afforded by mimetic resemblance to the parental species, are also largely outside the curves of parental mating preference. Disruptive sexual selection against F(1) hybrids therefore forms an additional post-mating barrier to gene flow, blurring the distinction between pre-mating and post-mating isolation, and helping to maintain the distinctness of these hybridizing species.     

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.     

A genetic linkage map of the mimetic butterfly Heliconius melpomene

Heliconius melpomene is a mimetic butterfly that exhibits great geographic variation in color pattern. We present here a genetic linkage map based on analysis of genetic markers in 73 individuals from a single F(2) family, offspring of a cross between H. m. cythera from western Ecuador and H. m. melpomene from French Guiana. A novel "three-step method" is described for the analysis of dominant markers in an F(2) cross, using outbred parental strains and taking advantage of the lack of crossing over in female Lepidoptera. This method is likely to prove useful for future mapping studies in outbred species with crossing over restricted to one sex, such as the Lepidoptera and Drosophila. The resulting linkage map has 21 linkage groups corresponding to the 21 chromosomes of H. melpomene and includes 219 AFLP markers, 23 microsatellites, 19 single-copy nuclear genes, and the color pattern switch genes Yb and Sb. The marker density is high, averaging >1/7 cM. The total map length is 1616 cM and the average chromosome length is 77 cM. The genome size of H. melpomene was estimated to be 292 Mb, giving a relationship of physical-to-map distance of 180 kb/cM. This map forms the basis for future comparative linkage analysis of color pattern evolution in Heliconius.     

Sharp genetic discontinuity across a unimodal Heliconius hybrid zone

Hybrid zones are powerful natural systems to study evolutionary processes to gain an understanding of adaptation and speciation. In the Cauca Valley (Colombia), two butterfly races, Heliconius cydno cydnides and Heliconius cydno weymeri, meet and hybridize. We characterized this hybrid zone using a combination of mitochondrial DNA (mtDNA) sequences, amplified fragment length polymorphisms (AFLPs), microsatellites and sequences for nuclear loci within and outside of the genomic regions that cause differences in wing colour pattern. The hybrid zone is largely composed of individuals of mixed ancestry. However, there is strong genetic discontinuity between the hybridizing races in mtDNA and, to a lesser extent, in all nuclear markers surveyed. The mtDNA clustering of H. c. cydnides with the H. cydno race from the Magdalena Valley and H. c. weymeri with the H. cydno race from the pacific coast suggests that H. c. cydnides colonized the Cauca Valley from the north, whereas H. c. weymeri did so by crossing the Andes in the southern part, implying a secondary contact origin. Colonization of the valley by H. cydno was accompanied by mimicry shift. Strong ecological isolation, driven by locally adaptive differences in mimetic wing patterns, is playing an important role in maintaining the hybrid zone. However, selection on wing pattern alone is not sufficient to explain the genetic discontinuity observed. There is evidence for differences in male mating preference, but the contribution of additional barriers needs further investigation. Overall, our results support the idea that speciation is a cumulative process, where the combination of multiple isolation barriers, combined with major phenotypic differences, facilitates population divergence in face of gene flow.     

Hybrid sterility,Haldane's rule and speciation in Heliconius cydno and H. melpomene

Most genetic studies of Haldane's rule, in which hybrid sterility or inviability affects the heterogametic sex preferentially, have focused on Drosophila. It therefore remains unclear to what extent the conclusions of that work apply more generally, particularly in female-heterogametic taxa such as birds and Lepidoptera. Here we present a genetic analysis of Haldane's rule in Heliconius butterflies. Female F(1) hybrids between Heliconius melpomene and H. cydno are completely sterile, while males have normal to mildly reduced fertility. In backcrosses of male F(1) hybrids, female offspring range from completely sterile to fully fertile. Linkage analysis using the Z-linked triose-phosphate isomerase locus demonstrates a "large X" (Z) effect on sterility. Expression of female sterility varies among crosses in this and a previous study of Heliconius. Sterility may result from the production of normal but infertile eggs, production of small infertile eggs, or from a complete failure to develop ovarioles, which suggests multiple routes to the evolution of hybrid sterility in these Heliconius species. These results conform to the expectations of the "dominance" rather than "faster male" theories of Haldane's rule and suggest that relatively few loci are responsible. The two species are broadly sympatric and hybridize in the wild, so that female hybrid sterility forms one of several strong but incomplete barriers to gene flow in nature. The effect of female sterility is comparable to that of selection against non-mimetic hybrids, while mate choice forms a much stronger barrier to gene transfer.     

Some early works on heliconiine butterflies and their biology (Lepidoptera, Nymphalidae)

In the course of a survey of the history of the study of Heliconiinae, particularly their biology, in the eighteenth and nineteenth centuries it is shown that Heliconius nigromarginatus (Goeze) and H. pallescens (Goeze) are nomina oblita , and that H. cinereofuscus (Goeze) is not, as always supposed, a species in its own right. A new name is proposed for H. melpomene cybele (Cramer). The origin and application of some of Linnaeus' names is discussed.
Illustrations of Heliconius by Seba, Clerck, Petiver and Cramer show that three species, melpomene, erato and doris have been polymorphic at least since the middle of the eighteenth century.
The Indo-Australian genus Cethosia should probably be included in the Heliconiinae.