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.     

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.     

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.     

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.     

A new mimetic species of Heliconius (Lepidoptera: Nymphalidae), from southeastern Colombia, revealed by cladistic analysis of mitochondrial DNA sequences

A new species of Heliconius and a new geographical race of Heliconius melpomene are described from the vicinity of Mocoa, Dpto. Putumayo, Colombia, based on molecular and morphological characters. The new species, H. tristero , is a close relative of H. cydno , a geographically differentiated species which lacks red coloration and engages in Müllerian mimicry with other blue and yellow Heliconius species in Central and northwestern South America. H. tristero has switched mimetic associations, instead mimicking the local, sympatric forms of two widespread mimetic species, H. erato and H. melpomene. This discovery provides evidence that the splinter species H. heurippa, H. tristero and H. timareta represent phenotypically divergent members of the H. cydno group that are endemic to successive river valleys on the eastern slope of the northern Andean Cordillera. The nominal taxon Heliconius amaryllis bellula Stichel, currently misapplied to both H. tristero and H. melpomene populations from the Mocoa region of Colombia, is considered here to represent a hybrid between H. heurippa and H. tristero. The Mocoa melpomene race is formally named Heliconius melpomene mocoa , new subspecies.     

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.     

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.     

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.     

Genomic tools and cDNA derived markers for butterflies

The Lepidoptera have long been used as examples in the study of evolution, but some questions remain difficult to resolve due to a lack of molecular genetic data. However, as technology improves, genomic tools are becoming increasingly available to tackle unanswered evolutionary questions. Here we have used expressed sequence tags (ESTs) to develop genetic markers for two Müllerian mimic species, Heliconius melpomene and Heliconius erato. In total 1363 ESTs were generated, representing 330 gene objects in H. melpomene and 431 in H. erato. User-friendly bioinformatic tools were used to construct a nonredundant database of these putative genes (available at http://www.heliconius.org), and annotate them with blast similarity searches, InterPro matches and Gene Ontology terms. This database will be continually updated with EST sequences for the Papilionideae as they become publicly available, providing a tool for gene finding in the butterflies. Alignments of the Heliconius sequences with putative homologues derived from Bombyx mori or other public data sets were used to identify conserved PCR priming sites, and develop 55 markers that can be amplified from genomic DNA in both H. erato and H. melpomene. These markers will be used for comparative linkage mapping in Heliconius and will have applications in other phylogenetic and genomic studies in the Lepidoptera.     

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.     

Phylogenetic study of heliconiine butterflies based on morphology and restriction analysis of ribosomal RNA genes

Relationships of ten heliconiine butterflies (genera Dryas and Heliconius , family Nymphalidae) were elucidated by phylogenetic analysis of characters based on ribosomal DNA restriction site variation and morphology. Agraulis vanillae , also a heliconiine, was used as the outgroup species. Although neither the morphological nor the molecular data unambiguously resolve relationships among the heliconiines, a combined analysis of both data sets results in a tree that is similar to traditional systematic arrangements and previous views of radiation in the group. Both pupal-mating and nonpupal-mating species group as clades in the combined analysis. However, the restriction site data alone do not support the monophyly of the pupal-mating clade, and the morphological data alone do not support the monophyly of the non-pupal-mating clade. Furthermore, relationships of H. melpomene, H. cydno and the silvaniform species depart from traditional arrangements based on morphology and reproductive compatibility experiments. All trees support the independent evolution of similar wing patterns of species previously suggested to be members of mimicry complexes. Several mimicry complexes appear to have a member in each of the two major monophyletic groups (pupal-mating and non-pupal-mating clades).     

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.     

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.     

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.     

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.     

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.     

Localization of Müllerian mimicry genes on a dense linkage map of Heliconius erato

We report a dense genetic linkage map of Heliconius erato, a neotropical butterfly that has undergone a remarkable adaptive radiation in warningly colored mimetic wing patterns. Our study exploited natural variation segregating in a cross between H. erato etylus and H. himera to localize wing color pattern loci on a dense linkage map containing amplified fragment length polymorphisms (AFLP), microsatellites, and single-copy nuclear loci. We unambiguously identified all 20 autosomal linkage groups and the sex chromosome (Z). The map spanned a total of 1430 Haldane cM and linkage groups varied in size from 26.3 to 97.8 cM. The average distance between markers was 5.1 cM. Within this framework, we localized two major color pattern loci to narrow regions of the genome. The first gene, D, responsible for red/orange elements, had a most likely placement in a 6.7-cM region flanked by two AFLP markers on the end of a large 87.5-cM linkage group. The second locus, Sd, affects the melanic pattern on the forewing and was found within a 6.3-cM interval between flanking AFLP loci. This study complements recent linkage analysis of H. erato's comimic, H. melpomene, and forms the basis for marker-assisted physical mapping and for studies into the comparative genetic architecture of wing-pattern mimicry in Heliconius.     

Synteny and chromosome evolution in the lepidoptera: evidence from mapping in Heliconius melpomene

The extent of conservation of synteny and gene order in the Lepidoptera has been investigated previously only by comparing a small subset of linkage groups between the moth Bombyx mori and the butterfly Heliconius melpomene. Here we report the mapping of 64 additional conserved genes in H. melpomene, which contributed 47 markers to a comparative framework of 72 orthologous loci spanning all 21 H. melpomene chromosomes and 27 of the 28 B. mori chromosomes. Comparison of the maps revealed conserved synteny across all chromosomes for the 72 loci, as well as evidence for six cases of chromosome fusion in the Heliconius lineage that contributed to the derived 21-chromosome karyotype. Comparisons of gene order on these fused chromosomes revealed two instances of colinearity between H. melpomene and B. mori, but also one instance of likely chromosomal rearrangement. B. mori is the first lepidopteran species to have its genome sequenced, and the finding that there is conserved synteny and gene order among Lepidoptera indicates that the genomic tools developed in B. mori will be broadly useful in other species.     

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.