STRONG NATURAL SELECTION IN A WARNING-COLOR HYBRID ZONE

Frequency-dependent selection on warning color can maintain narrow hybrid zones between unpalatable prey taxa. To measure such selection, we transferred marked Heliconius erato (Lepidoptera: Nymphalidae) in both directions across a 10-km-wide hybrid zone between Peruvian races differing in color pattern. These experimental H. erato were released at four sites, along with control H. erato of the phenotype native to each site. Survival of experimental butterflies was significantly lower than that of controls at two sites and overall. Most selection, measured as differences in survival, occurred soon after release. Selection against foreign morphs was 52% (confidence limits: 25–71 %) and was probably due to bird attacks on unusual warning-color morphs (more than 10% of the recaptures had beak marks). Since only three major loci determine the color-pattern differences, this suggests an average selection coefficient of 0.17 per locus, sufficient to maintain the narrow clines in H. erato.     

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.     

Movement of a Heliconius hybrid zone over 30 years: A Bayesian approach

Hybrid zones have long been of interest to biologists as natural laboratories where we can gain insight into the processes of adaptation and speciation. Repeated sampling of individual hybrid zones has been particularly useful in elucidating the dynamic balance between selection and dispersal that maintains most hybrid zones. Here, we revisit a hybrid zone between Heliconius erato butterflies in Panamá for a third time over more than 30 years. We combine a novel Bayesian extension of stepped‐cline hybrid zone models with environmental data to understand the genetic and environmental causes of cline dynamics in this species. The cline has continued to move west, likely due to dominance drive, but has slowed and broadened. Environmental analyses suggest that widespread deforestation in Panamá could be leading to decreased avian predation and relaxed selection, causing the observed changes in cline dynamics.     

GENETICS OF MIMICRY IN THE TIGER SWALLOWTAIL BUTTERFLIES,PAPILIO GLAUCUS AND P. CANADENSIS (LEPIDOPTERA: PAPILIONIDAE)

The tiger swallowtail butterfly, Papilio glaucus, exhibits a female-limited polymorphism for Batesian mimicry; the Canadian tiger swallowtail, Papilio canadensis, lacks the mimetic (dark) form entirely. The species hybridize to a limited extent where their ranges overlap. Field collections and censuses indicate that mimetic females occur throughout the range of P. glaucus but at lowest frequencies in populations at the latitudinal edges of its geographic range such as the southernmost part of Florida and along the entire northern edge of its distribution from Massachusetts to Minnesota. Frequencies of mimetic females have remained relatively stable over time. Inheritance of the mimetic form is controlled primarily by two interacting sex-linked loci. The typical matrilineal pattern of inheritance in P. glaucus can be explained by polymorphism at a Y-linked locus, b. Analysis of P. glaucus × P. canadensis crosses has also revealed an X-linked locus, s, which controls the expression of the mimetic phenotype. The P. canadensis allele, s^can, suppresses the mimetic phenotype in hybrid and backcross females. Results from more than 12 yr of rearing tiger swallowtails, including interspecies hybrids, indicate that the absence of mimetic P. canadensis females is due to both a high frequency of the “suppressing” allele s^can and low frequency of the black-pigment-determining b + allele. The frequency of s^can (or other suppressing alleles of s) in P. glaucus populations outside the hybrid zone is low. Some males heterozygous at the s locus and some suppressed mimetic females occur within the hybrid zone. A simple genetic model predicts the frequency of daughters that differ in phenotype from their mothers.     

High warning colour polymorphism in Heliconius hybrid zone roosts

1. Communal roosting behaviour has been documented among a wide range of taxa, particularly among groups of butterflies that display warning colourations. These aggregations of conspecifics and/or other species that share mimetic warning colour patterns can have a large impact on predator learning, and thus the survival of an aposematic form. Yet there has been limited investigation of communal roosting within areas where new and diverse warning colour forms are generated, such as hybrid zones. 2. Here, roosting behaviour was examined in a Heliconius erato hybrid zone in French Guiana between races with divergent warning colourations on their wings. In this hybrid zone, native individuals with nine distinct warning colourations, as well as individuals with altered forms that are not native in the French Guiana population, were marked and observed to determine if divergently coloured individuals participated in communal roosting, and if the proportions of colour pattern forms at roosts differed from the proportions that are found in the population. 3. The results demonstrated that divergently coloured individuals of the same species, including altered, non‐native forms, will readily and repeatedly participate in nocturnal communal roosting, often with extreme fidelity to specific perch locations. 4. These findings suggest that roosts composed of polymorphic warning patterns may be common in phenotypic transition zones, which could have major implications on predator training and selection dynamics in hybrid zones.     

Strong reproductive isolation and narrow genomic tracts of differentiation among three woodpecker species in secondary contact

Hybrid zones allow the measurement of gene flow across the genome, producing insight into the genomic architecture of speciation. Such analysis is particularly powerful when applied to multiple pairs of hybridizing species, as patterns of genomic differentiation can then be related to age of the hybridizing species, providing a view into the build‐up of differentiation over time. We examined 33 809 single nucleotide polymorphisms (SNPs) in three hybridizing woodpecker species: Red‐breasted, Red‐naped and Yellow‐bellied sapsuckers (Sphyrapicus ruber, Sphyrapicus nuchalis and Sphyrapicus varius), two of which (ruber and nuchalis) are much more closely related than each is to the third (varius). To identify positions of SNPs on chromosomes, we developed a localization method based on comparative genomics. We found narrow clines, bimodal distributions of hybrid indices and genomic regions with decreased rates of introgression. These results suggest moderately strong reproductive isolation among species and selection against specific hybrid genotypes. We found 19 small regions of strong differentiation between species, partly shared among species pairs, but no large regions of differentiation. An association analysis revealed a single strong‐effect candidate locus associated with plumage, possibly explaining mismatch among the three species in genomic relatedness and plumage similarity. Our comparative analysis of species pairs of different age and their hybrid zones showed that moderately strong reproductive isolation can occur with little genomic differentiation, but that reproductive isolation is incomplete even with much greater genomic differentiation, implying there are long periods of time when hybridization is possible if diverging populations are in geographic contact.     

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.     

Cross‐decades stability of an avian hybrid zone

Hybrid zones are particularly valuable for understanding the evolution of partial reproductive isolation between differentiated populations. An increasing number of hybrid zones have been inferred to move over time, but in most such cases zone movement has not been tested with long‐term genomic data. The hybrid zone between Townsend's Warblers (Setophaga townsendi) and Hermit Warblers (S. occidentalis) in the Washington Cascades was previously inferred to be moving from northern S. townsendi southwards towards S. occidentalis, based on plumage and behavioural patterns as well as a 2000‐km genetic wake of hermit mitochondrial DNA (mtDNA) in coastal Townsend's Warblers. We directly tested whether hybrid zone position has changed over 2–3 decades by tracking plumage, mtDNA and nuclear genomic variation across the hybrid zone over two sampling periods (1987–94 and 2015–16). Surprisingly, there was no significant movement in genomic or plumage cline centres between the two time periods. Plumage cline widths were narrower than expected by neutral diffusion, consistent with a ‘tension zone’ model, in which selection against hybrids is balanced by movement of parental forms into the zone. Our results indicate that this hybrid zone is either stable in its location or moving at a rate that is not detectable over 2–3 decades. Despite considerable gene flow, the stable clines in multiple phenotypic and genotypic characters over decades suggest evolutionary stability of this young pair of sister species, allowing divergence to continue. We propose a novel biogeographic scenario to explain these patterns: rather than the hybrid zone having moved thousands of kilometres to its current position, inland Townsend's met coastal Hermit Warbler populations along a broad front of the British Columbia and Alaska coast and hybridization led to replacement of the Hermit Warbler plumage with Townsend's Warbler plumage patterns along this coastline. Hence, hybrid zones along British Columbia and Alaska moved only a short distance from the inland to the coast, whereas the Hermit Warbler phenotype appears stable in Washington and further south. This case provides an example of the complex biogeographic processes that have led to the distribution of current phenotypes within and among closely related species.     

Müllerian mimicry of a quantitative trait despite contrasting levels of genomic divergence and selection

Hybrid zones, where distinct populations meet and interbreed, give insight into how differences between populations are maintained despite gene flow. Studying clines in genetic loci and adaptive traits across hybrid zones is a powerful method for understanding how selection drives differentiation within a single species, but can also be used to compare parallel divergence in different species responding to a common selective pressure. Here, we study parallel divergence of wing colouration in the butterflies Heliconius erato and H. melpomene, which are distantly related Müllerian mimics which show parallel geographic variation in both discrete variation in pigmentation, and quantitative variation in structural colour. Using geographic cline analysis, we show that clines in these traits are positioned in roughly the same geographic region for both species, which is consistent with direct selection for mimicry. However, the width of the clines varies markedly between species. This difference is explained in part by variation in the strength of selection acting on colour traits within each species, but may also be influenced by differences in the dispersal rate and total strength of selection against hybrids between the species. Genotyping‐by‐sequencing also revealed weaker population structure in H. melpomene, suggesting the hybrid zones may have evolved differently in each species, which may also contribute to the patterns of phenotypic divergence in this system. Overall, we conclude that multiple factors are needed to explain patterns of clinal variation within and between these species, although mimicry has probably played a central role.     

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 zone origins, species boundaries, and the evolution of wing-pattern diversity in a polytypic species complex of North American admiral butterflies (Nymphalidae: Limenitis)

Hybrid zones present opportunities to study the effects of gene flow, selection, and recombination in natural populations and, thus, provide insights into the genetic and phenotypic changes that occur early in speciation. Here we investigate a hybrid zone between mimetic (Limenitis arthemis astyanax) and nonmimetic (Limenitis arthemis arthemis) populations of admiral butterflies using DNA sequence variation from mtDNA and seven nuclear gene loci. We find three distinct mitochondrial clades within this complex, and observe a strong overall concordance between wing-pattern phenotypes and mitochondrial variation. Nuclear gene genealogies, in contrast, revealed no evidence of exclusivity for either wing-pattern phenotype, suggesting incomplete barriers to gene exchange and/or insufficient time for lineage sorting. Coalescent simulations indicate that gene flow between these two subspecies is highly asymmetric, with the majority of migration occurring from mimetic into nonmimetic populations. Selective sweeps of alleles responsible for mimetic phenotypes may have occurred more than once when mimetic and nonmimetic Limenitis occurred together in the presence of the model (Battus philenor).     

Cryptic differences in colour among Müllerian mimics: how can the visual capacities of predators and prey shape the evolution of wing colours?

Antagonistic interactions between predators and prey often lead to co‐evolution. In the case of toxic prey, aposematic colours act as warning signals for predators and play a protective role. Evolutionary convergence in colour patterns among toxic prey evolves due to positive density‐dependent selection and the benefits of mutual resemblance in spreading the mortality cost of educating predators over a larger prey assemblage. Comimetic species evolve highly similar colour patterns, but such convergence may interfere with intraspecific signalling and recognition in the prey community, especially for species involved in polymorphic mimicry. Using spectrophotometry measures, we investigated the variation in wing coloration among comimetic butterflies from distantly related lineages. We focused on seven morphs of the polymorphic species Heliconius numata and the seven corresponding comimetic species from the genus Melinaea. Significant differences in the yellow, orange and black patches of the wing were detected between genera. Perceptions of these cryptic differences by bird and butterfly observers were then estimated using models of animal vision based on physiological data. Our results showed that the most strikingly perceived differences were obtained for the contrast of yellow against a black background. The capacity to discriminate between comimetic genera based on this colour contrast was also evaluated to be higher for butterflies than for birds, suggesting that this variation in colour, likely undetectable to birds, might be used by butterflies for distinguishing mating partners without losing the benefits of mimicry. The evolution of wing colour in mimetic butterflies might thus be shaped by the opposite selective pressures exerted by predation and species recognition.     

REINFORCING SELECTION IS EFFECTIVE UNDER A RELATIVELY BROAD SET OF CONDITIONS IN A MOSAIC HYBRID ZONE

Many hybrid zones have a mosaic structure, yet we know of no theoretical work that examines the impact of mosaicism on the outcome of evolution. We developed a computer simulation model designed to test whether the outcome of reinforcing selection differs in a mosaic and a clinal hybrid zone. Our model was a one-dimensional stepping-stone model. The mosaic and clinal hybrid zones that we modeled were, respectively, a mosaic maintained by differential fitness of the interacting taxa in patchy habitats and a tension zone. We modeled changes in gene frequency at two biallelic loci, A and B. Hybrids at the A locus were selected against. An allele at the B locus caused assortative mating at the A locus, which promoted reinforcement; there was a selective cost to this allele. In a mosaic hybrid zone, spatial variation in the fitness of A-locus homozygotes in different patches caused gene and genotype frequencies at the A and B loci to differ greatly from those in a tension zone. Compared to a tension zone, a mosaic hybrid zone had a broader region in which hybrids could be formed and, thus, a broader region in which the assortative-mating allele provided a net selective advantage (via decreased production of the less fit A-locus hybrids). This caused the assortative-mating allele to be favored under a broader set of conditions in a mosaic hybrid zone than in a tension zone. In mosaic and tension hybrid zones, both low and high levels of migration could prevent the establishment of the allele that promoted reinforcement, but the allele could establish under a wider range of migration rates in a mosaic than in a tension zone. In a tension zone, both low and high levels of selection against A-locus hybrids could prevent the establishment of the assortative-mating allele. In a mosaic hybrid zone, the assortative-mating allele established under lower levels of selection against hybrids than in a tension zone, and high levels of selection did not impede the establishment of this allele. Overall, our work illustrates how the structure of a hybrid zone can alter the outcome of an important evolutionary process, in this case, reinforcement.     

Recombinant hybrids retain heterozygosity at many loci: new insights into the genomics of reproductive isolation in Populus

The maintenance of species barriers in the face of gene flow is often thought to result from strong selection against intermediate genotypes, thereby preserving genetic differentiation. Most speciation genomic studies thus aim to identify exceptionally divergent loci between populations, but divergence will be affected by many processes other than reproductive isolation (RI) and speciation. Through genomic studies of recombinant hybrids sampled in the wild, genetic variation associated with RI can be observed in situ, because selection against incompatible genotypes will leave detectable patterns of variation in the hybrid genomes. To better understand the mechanisms directly involved in RI, we investigated three natural ‘replicate’ hybrid zones between two divergent Populus species via locus‐specific patterns of ancestry across recombinant hybrid genomes. As expected, genomic patterns in hybrids and their parental species were consistent with the presence of underdominant selection at several genomic regions. Surprisingly, many loci displayed greatly increased between‐species heterozygosity in recombinant hybrids despite striking genetic differentiation between the parental genomes, the opposite of what would be expected with selection against intermediate genotypes. Only a limited, reproducible set of genotypic combinations was present in hybrid genomes across localities. In the absence of clearly delimited ‘hybrid habitats’, our results suggest that complex epistatic interactions within genomes play an important role in advanced stages of RI between these ecologically divergent forest trees. This calls for more genomic studies that test for unusual patterns of genomic ancestry in hybridizing species.     

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.     

First chromosome scale genomes of ithomiine butterflies (Nymphalidae: Ithomiini): Comparative models for mimicry genetic studies

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies.     

Spatial mosaic formation through frequency-dependent selection in Müllerian mimicry complexes

Although contemporary models of Müllerian mimicry have considered the movement of interfacial boundaries between two distinct mimetic forms, and even the possibility of polymorphisms in two patch systems, no model has considered how multiple forms of Müllerian mimics might evolve and be maintained over large geographical areas. A spatially explicit individual-based model for the evolution of Müllerian mimicry is presented, in which two unpalatable species are distributed over discrete cells within a regular lattice. Populations in each cell are capable of genetic drift and experience localized dispersal as well as frequency-dependent selection by predators. When each unpalatable prey species was introduced into a random cell and allowed to spread, then mimicry evolved throughout the system in the form of a spatial mosaic of phenotypes, separated by narrow "hybrid zones". The primary mechanism generating phenotypic diversity was the occasional establishment of new mutant forms in unoccupied cells and their subsequent maintenance (and spread) through frequency-dependent selection. The mean number of discrete clusters of the same morph that formed in the lattice was higher the higher the intensity of predation, and higher the lower the dispersal rate of unpalatable prey. Under certain conditions the hybrid zones moved, in a direction dependent on the curvature of their interfacial boundaries. However, the mimetic mosaics were highly stable when the intensity of predation was high and the rate of prey dispersal was low. Overall, this model highlights how a stable mosaic of different mimetic forms can evolve from a range of starting conditions through a combination of chance effects and localized frequency-dependent selection.     

EVOLUTIONARY IMPLICATIONS OF DIVERGENT CLINES IN AN AVIAN (MANACUS: AVES) HYBRID ZONE

Abstract A previous study of the hybrid zone in western Panama between white‐collared (Manacus candei) and golden‐collared manakins (M. vitellinus) documented the unidirectional introgression of vitellinus male secondary sexual traits across the zone. Here, we examine the hybrid zone in greater genetic and morphological detail. Statistical comparisons of clines are performed using maximum‐likelihood and nonparametric bootstrap methods. Our results demonstrate that an array of six molecular and two morphometric markers agree in cline position and width. Clines for male collar and belly color are similar in width to the first eight clines, but are shifted in position by at least five cline widths. The result is that birds in intervening populations are genetically and morphometrically very like parental candei, but males have the plumage color of parental vitellinus. Neither neutral diffusion nor nonlinearity of color scales appear to be viable explanations for the large cline shifts. Genetic dominance of vitellinus plumage traits is another potential explanation that will require breeding experiments to test. Sexual selection remains a plausible explanation for the observed introgression of vitellinus color traits in these highly dimorphic, polygynous, lek‐mating birds. Two other clines, including a nondiagnostic isozyme locus, are similar in position to the main cluster of clines, but are broader in width. Thus, introgression at some loci is greater than that detected with diagnostic markers. Assuming that narrow clines are maintained by selection, variation in cline width indicates that selection is not uniform throughout the genome and that diagnostic markers are under more intense selective pressure. The traditional focus on diagnostic markers in studies of hybrid zones may therefore lead to underestimates of average introgression. This effect may be more pronounced in organisms with low levels of genetic divergence between hybridizing taxa.     

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.     

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.