Beyond magic traits: Multimodal mating cues in Heliconius butterflies

Species coexistence involves the evolution of reproductive barriers opposing gene flow. Heliconius butterflies display colorful patterns affecting mate choice and survival through warning signaling and mimicry. These patterns are called “magic traits” for speciation because divergent natural selection may promote mimicry shifts in pattern whose role as mating cue facilitates reproductive isolation. By contrast, between comimetic species, natural selection promotes pattern convergence. We addressed whether visual convergence interferes with reproductive isolation by testing for sexual isolation between two closely related species with similar patterns, H. timareta thelxinoe and H. melpomene amaryllis. Experiments with models confirmed visual attraction based on wing phenotype, leading to indiscriminate approach. Nevertheless, mate choice experiments showed assortative mating. Monitoring male behavior toward live females revealed asymmetry in male preference, H. melpomene males courting both species equally while H. timareta males strongly preferred conspecifics. Experiments with hybrid males suggested an important genetic component for such asymmetry. Behavioral observations support a key role for short‐distance cues in determining male choice in H. timareta. Scents extracts from wings and genitalia revealed interspecific divergence in chemical signatures, and hybrid female scent composition was significantly associated with courtship intensity by H. timareta males, providing candidate chemical mating cues involved in sexual isolation.     

Development of microsatellite loci from a reference genome for the Neotropical butterfly Heliconius numata and its close relatives

The Neotropical butterfly Heliconius numata (Lepidoptera: Nymphalidae: Heliconiinae) is known for its striking diversity of wing color patterns driven by the Müllerian mimicry of multiple local models and controlled by a single supergene locus. Such fine‐scale variation of traits under strong selection offers a unique opportunity for the study of the ecology and genetics of adaptation. However, little is still known of the population processes driving geographical variation in wing‐pattern phenotypes. We report the characterization of 26 microsatellite markers for the butterfly H. numata, including six located inside the wing color‐pattern supergene region. All markers are polymorphic, with allele numbers ranging from 2 to 21 per locus, an observed heterozygosity of 0.111 to 0.848 and an expected heterozygosity of 0.126 to 0.942. A subset of 18 of these markers was tested on five closely related sympatric Heliconius species with an amplification success ranging from 88% to 94%. The obtained set of microsatellite markers provides a new and useful set of tools to investigate patterns of differentiation and selection in populations of mimetic Heliconius butterflies. Moreover, markers developed within the color‐pattern supergene will facilitate characterization of the association between the genetic architecture and the functional diversity of wing patterns. Finally, the cross‐species amplification success of the described markers extends their utility to also encompass comparative population genetic studies of closely related species within a clade of rapidly diversifying species.     

Subtle variation in size and shape of the whole forewing and the red band among co‐mimics revealed by geometric morphometric analysis in Heliconius butterflies

Heliconius are unpalatable butterflies that exhibit remarkable intra‐ and interspecific variation in wing color pattern, specifically warning coloration. Species that have converged on the same pattern are often clustered in Müllerian mimicry rings. Overall, wing color patterns are nearly identical among co‐mimics. However, fine‐scale differences exist, indicating that factors in addition to natural selection may underlie wing phenotype. Here, we investigate differences in shape and size of the forewing and the red band in the Heliconius postman mimicry ring (H. erato phyllis and the co‐mimics H. besckei, H. melpomene burchelli, and H. melpomene nanna) using a landmark‐based approach. If phenotypic evolution is driven entirely by predation pressure, we expect nonsignificant differences among co‐mimics in terms of wing shape. Also, a reinforcement of wing pattern (i.e., greater similarity) could occur when co‐mimics are in sympatry. We also examined variation in the red forewing band because this trait is critical for both mimicry and sexual communication. Morphometric results revealed significant but small differences among species, particularly in the shape of the forewing of co‐mimics. Although we did not observe greater similarity when co‐mimics were in sympatry, nearly identical patterns provided evidence of convergence for mimicry. In contrast, mimetic pairs could be distinguished based on the shape (but not the size) of the red band, suggesting an “advergence” process. In addition, sexual dimorphism in the red band shape (but not size) was found for all lineages. Thus, we infer that natural selection due to predation by birds might not be the only mechanism responsible for variation in color patterns, and sexual selection could be an important driver of wing phenotypic evolution in this mimicry ring.     

Why are there so many mimicry rings? Correlations between habitat,behaviour and mimicry in Heliconius butterflies

In the new world tropics there is an extravagant array of sympatric butterfly mimicry rings. This is puzzling under strictly coevolutionary (Müllerian) mimicry: all unpalatable species should converge as ‘co-mimics' to the same pattern. If mimicry has usually evolved in unpalatable species by one-sided (Batesian) evolution, however, it is easy to see that mimicry rings centred on different models could remain distinct. If mimicry rings were also segregated by habitat, a diversity of mimicry rings could be stabilized. In this paper we report correlations between behaviour and mimicry of nine unpalatable Heliconius species. It is already known that co-mimics fly in similar habitats, and non-mimics fly in different habitats, although there is much overlap. Contrary to a previous report, we find little difference in flight heights of heliconiine mimicry rings; all species fly from ground level to the canopy. However, co-mimics roost at night in similar habitats and at similar heights above the ground, but in different habitats and at different heights from species in other mimicry rings. Heliconius (especially the erato taxonomic group) are renowned for roosting gregariously; and co-mimics roost gregariously with each other more often than with non-mimics. Gregarious roosting is therefore common between species, as well as within species. There are thus strong links between mimicry and behavioural ecology in Heliconius. The paradoxical correlation between nocturnal roosting and visual mimicry is presumably explained by bird predation at dusk when roosts are forming, or at dawn before they have disbanded. Direct evidence of predation is lacking, but there are high rates of disturbance by birds at these times. These results, together with knowledge of the phylogeny of Heliconius, suggest that species from the melpomene-group of Heliconius have radiated to occupy mimetic niches protected by model species in the Ithomiinae and the erato-group of Heliconius. A variety of sympatric mimicry rings is apparently maintained because key models fail to converge, while more rapidly-evolving unpalatable mimics evolve towards the colour patterns of the models. The maintenance of mimetic diversity would be aided by the habitat and behavioural differences between mimicry rings revealed here, provided that different predators are found in different habitats. This explanation for the maintenance of multiple mimicry rings is more plausible for Heliconius mimicry than alternatives based on visual mating constraints, thermal ecology, or camouflage.     

Wing patterning genes and coevolution of Müllerian mimicry in Heliconius butterflies: Support from phylogeography,cophylogeny, and divergence times

Examples of long‐term coevolution are rare among free‐living organisms. Müllerian mimicry in Heliconius butterflies had been suggested as a key example of coevolution by early genetic studies. However, research over the last two decades has been dominated by the idea that the best‐studied comimics, H. erato and H. melpomene, did not coevolve at all. Recently sequenced genes associated with wing color pattern phenotype offer a new opportunity to resolve this controversy. Here, we test the hypothesis of coevolution between H. erato and H. melpomene using Bayesian multilocus analysis of five color pattern genes and five neutral genetic markers. We first explore the extent of phylogenetic agreement versus conflict between the different genes. Coevolution is then tested against three aspects of the mimicry diversifications: phylogenetic branching patterns, divergence times, and, for the first time, phylogeographic histories. We show that all three lines of evidence are compatible with strict coevolution of the diverse mimicry wing patterns, contrary to some recent suggestions. Instead, these findings tally with a coevolutionary diversification driven primarily by the ecological force of Müllerian mimicry.     

Altitude and life‐history shape the evolution of Heliconius wings

Phenotypic divergence between closely related species has long interested biologists. Taxa that inhabit a range of environments and have diverse natural histories can help understand how selection drives phenotypic divergence. In butterflies, wing color patterns have been extensively studied but diversity in wing shape and size is less well understood. Here, we assess the relative importance of phylogenetic relatedness, natural history, and habitat on shaping wing morphology in a large dataset of over 3500 individuals, representing 13 Heliconius species from across the Neotropics. We find that both larval and adult behavioral ecology correlate with patterns of wing sexual dimorphism and adult size. Species with solitary larvae have larger adult males, in contrast to gregarious Heliconius species, and indeed most Lepidoptera, where females are larger. Species in the pupal‐mating clade are smaller than those in the adult‐mating clade. Interestingly, we find that high‐altitude species tend to have rounder wings and, in one of the two major Heliconius clades, are also bigger than their lowland relatives. Furthermore, within two widespread species, we find that high‐altitude populations also have rounder wings. Thus, we reveal novel adaptive wing morphological divergence among Heliconius species beyond that imposed by natural selection on aposematic wing coloration.     

Hybridization and transgressive exploration of colour pattern and wing morphology in Heliconius butterflies

Hybridization can generate novel phenotypes distinct from those of parental lineages, a phenomenon known as transgressive trait variation. Transgressive phenotypes might negatively or positively affect hybrid fitness, and increase available variation. Closely related species of Heliconius butterflies regularly produce hybrids in nature, and hybridization is thought to play a role in the diversification of novel wing colour patterns despite strong stabilizing selection due to interspecific mimicry. Here, we studied wing phenotypes in first‐ and second‐generation hybrids produced by controlled crosses between either two co‐mimetic species of Heliconius or between two nonmimetic species. We quantified wing size, shape and colour pattern variation and asked whether hybrids displayed transgressive wing phenotypes. Discrete traits underlain by major‐effect loci, such as the presence or absence of colour patches, generate novel phenotypes. For quantitative traits, such as wing shape or subtle colour pattern characters, hybrids only exceed the parental range in specific dimensions of the morphological space. Overall, our study addresses some of the challenges in defining and measuring phenotypic transgression for multivariate traits and our data suggest that the extent to which transgressive trait variation in hybrids contributes to phenotypic diversity depends on the complexity and the genetic architecture of the traits.     

Evolution of novel mimicry rings facilitated by adaptive introgression in tropical butterflies

Understanding the genetic basis of phenotypic variation and the mechanisms involved in the evolution of adaptive novelty, especially in adaptive radiations, is a major goal in evolutionary biology. Here, we used whole‐genome sequence data to investigate the origin of the yellow hindwing bar in the Heliconius cydno radiation. We found modular variation associated with hindwing phenotype in two narrow noncoding regions upstream and downstream of the cortex gene, which was recently identified as a pigmentation pattern controller in multiple species of Heliconius. Genetic variation at each of these modules suggests an independent control of the dorsal and ventral hindwing patterning, with the upstream module associated with the ventral phenotype and the downstream module with the dorsal one. Furthermore, we detected introgression between H. cydno and its closely related species Heliconius melpomene in these modules, likely allowing both species to participate in novel mimicry rings. In sum, our findings support the role of regulatory modularity coupled with adaptive introgression as an elegant mechanism by which novel phenotypic combinations can evolve and fuel an adaptive radiation.     

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

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

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.     

Plasticity and genetic effects contribute to different axes of neural divergence in a community of mimetic Heliconius butterflies

Changes in ecological preference, often driven by spatial and temporal variation in resource distribution, can expose populations to environments with divergent information content. This can lead to adaptive changes in the degree to which individuals invest in sensory systems and downstream processes, to optimize behavioural performance in different contexts. At the same time, environmental conditions can produce plastic responses in nervous system development and maturation, providing an alternative route to integrating neural and ecological variation. Here, we explore how these two processes play out across a community of Heliconius butterflies. Heliconius communities exhibit multiple Mullerian mimicry rings, associated with habitat partitioning across environmental gradients. These environmental differences have previously been linked to heritable divergence in brain morphology in parapatric species pairs. They also exhibit a unique dietary adaptation, known as pollen feeding, that relies heavily on learning foraging routes, or trap-lines, between resources, which implies an important environmental influence on behavioural development. By comparing brain morphology across 133 wild-caught and insectary-reared individuals from seven Heliconius species, we find strong evidence for interspecific variation in patterns of neural investment. These largely fall into two distinct patterns of variation; first, we find consistent patterns of divergence in the size of visual brain components across both wild and insectary-reared individuals, suggesting genetically encoded divergence in the visual pathway. Second, we find interspecific differences in mushroom body size, a central component of learning and memory systems, but only among wild caught individuals. The lack of this effect in common-garden individuals suggests an extensive role for developmental plasticity in interspecific variation in the wild. Finally, we illustrate the impact of relatively small-scale spatial effects on mushroom body plasticity by performing experiments altering the cage size and structure experienced by individual H. hecale. Our data provide a comprehensive survey of community level variation in brain structure, and demonstrate that genetic effects and developmental plasticity contribute to different axes of interspecific neural variation.     

Contrasting genomic and phenotypic outcomes of hybridization between pairs of mimetic butterfly taxa across a suture zone

Hybrid zones, whereby divergent lineages come into contact and eventually hybridize, can provide insights on the mechanisms involved in population differentiation and reproductive isolation, and ultimately speciation. Suture zones offer the opportunity to compare these processes across multiple species. In this paper we use reduced‐complexity genomic data to compare the genetic and phenotypic structure and hybridization patterns of two mimetic butterfly species, Ithomia salapia and Oleria onega (Nymphalidae: Ithomiini), each consisting of a pair of lineages differentiated for their wing colour pattern and that come into contact in the Andean foothills of Peru. Despite similarities in their life history, we highlight major differences, both at the genomic and phenotypic level, between the two species. These differences include the presence of hybrids, variations in wing phenotype, and genomic patterns of introgression and differentiation. In I. salapia, the two lineages appear to hybridize only rarely, whereas in O. onega the hybrids are not only more common, but also genetically and phenotypically more variable. We also detected loci statistically associated with wing colour pattern variation, but in both species these loci were not over‐represented among the candidate barrier loci, suggesting that traits other than wing colour pattern may be important for reproductive isolation. Our results contrast with the genomic patterns observed between hybridizing lineages in the mimetic Heliconius butterflies, and call for a broader investigation into the genomics of speciation in Ithomiini ‐ the largest radiation of mimetic butterflies.     

PARALLEL RACE FORMATION AND THE EVOLUTION OF MIMICRY IN HELICONIUS BUTTERFLIES: A PHYLOGENETIC HYPOTHESIS FROM MITOCHONDRIAL DNA SEQUENCES

Mimicry has been a fundamental focus of research since the birth of evolutionary biology yet rarely has been studied from a phylogenetic perspective beyond the simple recognition that mimics are not similar due to common descent. The difficulty of finding characters to discern relationships among closely related and convergent taxa has challenged systematists for more than a century. The phenotypic diversity of wing pattens among mimetic Heliconius adds an additional twist to the problem, because single species contain more than a dozen radically different-looking geographical races even though the mimetic advantage is theoretically highest when all individuals within and between species appear the same. Mitochondrial DNA (mtDNA) offers an independent way to address these issues. In this study, Cytochrome Oxidase I and II sequences from multiple, parallel races of Heliconius erato and Heliconius melpomene are examined, to estimate intraspecific phylogeny and gauge sequence divergence and ages of clades among races within each species. Although phenotypes of sympatric races exhibit remarkable concordance between the two species, the mitochondrial cladograms show that the species have not shared a common evolutionary history. H. erato exhibits a basal split between trans- and cis-Andean groups of races, whereas H. melpomene originates in the Guiana Shield. Diverse races in either species appear to have evolved within the last 200,000 yr, and convergent phenotypes have evolved independently within as well as between species. These results contradict prior theories of the evolution of mimicry based on analysis of wing-pattern genetics.     

Primers for the amplification of nuclear introns in Heliconius butterflies

I described the development of degenerate polymerase chain reaction (PCR) primers for intron‐containing regions of nine candidate wing patterning and pigmentation genes in Heliconius butterflies. Primers were developed by comparing sequence data from Drosophila melanogaster, Precis coenia, and a variety of other insects so they are likely to be applicable widely among the butterfly family Nymphalidae and perhaps Lepidoptera in general. The amplified regions are highly variable and should be useful for inferring relationships among closely related species and estimating the phylogeographical and population genetic structure of individual species.     

Evaluating female remating rates in light of spermatophore degradation in Heliconius butterflies: pupal‐mating monandry versus adult‐mating polyandry

1. Butterflies are frequently used in comparative studies of sexual selection because of their diverse mating systems. In Heliconius, the two major clades in the genus are characterised by contrasting pupal‐mating and adult‐mating strategies. Adult‐mating females are considered to be promiscuous whereas pupal‐mating females are thought to be monandrous. 2. Counting spermatophores in female Lepidoptera is a common method for assessing patterns of female remating. However, in pupal‐mating Heliconius butterflies spermatophores can become completely degraded, potentially leading to underestimation of female remating rates. 3. We qualitatively characterised the different states of spermatophore degradation, and showed that complete degradation takes approximately 3 weeks in captive‐bred H. erato females. 4. We counted spermatophores and/or assayed spermatophore degradation in > 500 Heliconius females across 28 species sampled from natural populations. Among pupal‐maters these observations yielded a few rare observations of double mating by recently eclosed females, but generally indicated a lack of rematings. In contrast, approximately 25% of sampled adult‐mating females remated at least once. 5. Using a novel statistical analysis, we estimated remating rates from patterns of spermatophore degradation or from counts stratified by age, as indicated by wing‐wear. This analysis showed no statistically significant evidence for remating for the pupal‐mating H. erato whereas significant remating rates were detected for adult‐mating species. 6. The present results support the established view of Heliconius mating systems in which pupal‐maters are largely monandrous, whereas adult‐maters are polyandrous.     

Molecular patterns of differentiation in canyon treefrogs (Hyla arenicolor): evidence for introgressive hybridization with the Arizona treefrog (H. wrightorum) and correlations with advertisement call differences

Detection of genetic and behavioural diversity within morphologically similar species has led to the discovery of cryptic species complexes. We tested the hypothesis that US populations of the canyon treefrog (Hyla arenicolor) may consist of cryptic species by examining mate‐attraction signals among three divergent clades defined by mtDNA. Using a multi‐locus approach, we re‐analysed phylogenetic relationships among the three clades and a closely related, but morphologically and behaviourally dissimilar species, the Arizona treefrog (H. wrightorum). We found evidence for introgression of H. wrightorum’s mitochondrial genome into H. arenicolor. Additionally, the two‐clade topology based on nuclear data is more congruent with patterns of call variation than the three‐clade topology from the mitochondrial dataset. The magnitude of the call divergence is probably insufficient to promote isolation of the nuclear DNA‐defined clades should they become sympatric, but further divergence in call properties significant in species identification could promote speciation in the future.     

Variable Selection and the Coexistence of Multiple mimetic forms of the Butterfly Heliconius numata

Polymorphism in aposematic animals and coexistence of multiple mimicry rings within a habitat are not predicted by classical Müllerian mimicry. The butterfly Heliconius numata Cramer (Lepidoptera: Nymphalidae; Heliconiinae) is both polymorphic and aposematic. The polymorphism is due to variation at a single locus (or `supergene') which determines colour patterns involved in Müllerian mimicry. We sampled 11 sites in a small area (approx. 60×30km) of North-eastern Peru for H. numata and its co-mimics in the genus Melinaea and Athyrtis (Ithomiinae), and examined the role of temporal and spatial heterogeneity in the maintenance of polymorphism. Colour-patterns of Melinaea communities, which constitute the likely `mimetic environment' for H. numata, are differentiated on a more local scale than morphs of H. numata, but the latter do show a strong and significant response to local selection for colour-pattern. In contrast, analysis of enzyme polymorphism in H. numata across the region revealed no spatial structure, which is consistent with a high mobility of this species. Differences in spatial variability in the two taxa may have caused H. numata to become polymorphic, while temporal variability, not significant in this study, probably has a lesser effect. The mimetic polymorphism is therefore explained by means of multiple selection-migration clines at a single locus, a similar process to that which explains narrow hybrid zones between geographic races of other Heliconius butterflies. This revised version was published online in July 2006 with corrections to the Cover Date.     

Mimetic colour pattern evolution in the highly polymorphic Bombus trifasciatus (Hymenoptera: Apidae) species complex and its comimics

Müllerian mimetic systems have uncovered some of the dynamic processes by which natural selection can drive the radiation of convergent and divergent phenotypes. We examined evolution involving Müllerian mimicry in bumble bees by documenting the distribution and evolution of colour patterns amongst three colour‐polymorphic lineages –Bombus trifasciatus Smith, Bombus haemorrhoidalis Smith, and Bombus breviceps Smith – that mimic each other across ～14 colour groups in South‐East Asia. Using mitochondrial DNA sequence data, we estimated relationships within each lineage to infer the processes that gave rise to the colour diversity and develop hypotheses on species recognition. We expanded on our assessment of species delineation in the B. trifasciatus lineage using three nuclear gene fragments and morphometrics. Comparison of colour patterns amongst georeferenced specimens showed considerable variation in the degree and geographical range of mimicry amongst mimicry groups. Phylogenetic estimates show high rates of colour pattern evolution, with colour variation often exceeding variation within the fast‐evolving mitochondrial genes. The molecular data, and to some degree the morphometric data, support unique histories for several taxa recognized previously within the B. trifasciatus lineage, which may include several species. Early vicariant events within the B. trifasciatus lineage are likely to have occurred ～2.2 Mya in the mountains of south‐west China. © 2012 The Linnean Society of London, Zoological Journal of the Linnean Society, 2012, 166 , 805–826.     

Clustering of loci controlling species differences in male chemical bouquets of sympatric Heliconius butterflies

The degree to which loci promoting reproductive isolation cluster in the genome—that is, the genetic architecture of reproductive isolation—can influence the tempo and mode of speciation. Tight linkage between these loci can facilitate speciation in the face of gene flow. Pheromones play a role in reproductive isolation in many Lepidoptera species, and the role of endogenously produced compounds as secondary metabolites decreases the likelihood of pleiotropy associated with many barrier loci. Heliconius butterflies use male sex pheromones to both court females (aphrodisiac wing pheromones) and ward off male courtship (male‐transferred antiaphrodisiac genital pheromones), and it is likely that these compounds play a role in reproductive isolation between Heliconius species. Using a set of backcross hybrids between H. melpomene and H. cydno, we investigated the genetic architecture of putative male pheromone compound production. We found a set of 40 significant quantitative trait loci (QTL) representing 33 potential pheromone compounds. QTL clustered significantly on two chromosomes, chromosome 8 for genital compounds and chromosome 20 for wing compounds, and chromosome 20 was enriched for potential pheromone biosynthesis genes. There was minimal overlap between pheromone QTL and known QTL for mate choice and color pattern. Nonetheless, we did detect linkage between a QTL for wing androconial area and optix, a color pattern locus known to play a role in reproductive isolation in these species. This tight clustering of putative pheromone loci might contribute to coincident reproductive isolating barriers, facilitating speciation despite ongoing gene flow.     

Inter-species variation in unpalatability does not explain polymorphism in a mimetic species

Conspicuous colouration in unpalatable organisms acts as a warning signal of their unprofitability, a phenomenon known as aposematism. The protection conferred by such colouration can lead to evolutionary convergence in warning signals between aposematic species, because sharing warning signals reduces the per capita cost of predator learning. Consequently, most aposematic species display a single colour pattern and participate in a single mimetic community (i.e. mimicry ring) at any given locality. However, some, like the Amazonian butterfly Heliconius numata, are polymorphic and participate in several mimicry rings within the same locality. We tested whether the unexpected polymorphism of H. numata could be due to a weak defence against predators. Poorly defended species participating in a mimicry ring are subject to negative frequency dependent selection, because their presence weakens the protection provided by the shared signal. This could promote polymorphism and participation in multiple mimicry rings. Using wild caught great tits (Parus major), we compared the palatability of H. numata to one of its locally monomorphic co-mimics (Mechanitis polymnia) and to two other locally monomorphic Heliconius species (H. melpomene and H. erato). The tested birds strongly rejected the polymorphic species H. numata, as well as the two other Heliconius species. Unexpectedly, a significantly weaker rejection was found towards M. polymnia, which relies on different toxic compounds to Heliconius. Our study demonstrates that the origin of polymorphic mimicry in H. numata is unlikely to stem from low unpalatability and raises new questions on defence variation within mimetic communities.