STRONG NATURAL SELECTION IN A WARNING-COLOR HYBRID ZONE

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.

Polymorphic mimicry, microhabitat use, and sex-specific behaviour

In order to assess the adaptive importance of microhabitat segregation for the maintenance of mimetic diversity, I explore how flight height varies between the sympatric forms of the polymorphic butterfly Heliconius numata and their respective models in the genus Melinaea. There is no evidence for vertical stratification of mimicry rings in these tiger-patterned butterflies, but males of H. numata tend to fly significantly higher than females and the Melinaea models. This difference in microhabitat preference likely results from females searching for host plants whereas males are patrolling for mates. I then present an extension of Muller's mimicry model for the case of partial behavioural or spatial segregation of sexes. The analysis suggests that sex-specific behaviours can make mimicry more beneficial, simply by reducing the effective population size participating in mimicry. The interaction between mimicry and sex-specific behaviours may therefore facilitate the evolution of polymorphism via enhanced, fine-scale local adaptation.

Reproductive isolation caused by visual predation on migrants between divergent environments

In theory, natural selection can drive adaptation within species while simultaneously promoting the formation of new species by causing the evolution of reproductive isolation. Cryptic coloration is widespread in nature and is generally considered to be a clear visual example of adaptation. I provide evidence that population divergence in cryptic coloration can also cause reproductive isolation. First, a manipulative field experiment using walking-stick insects demonstrates that the relative survival of different colour-pattern morphs depends on the host-plant species on which they are resting, but only in the presence of avian predation. Second, natural populations adapted to different host plants have diverged in colour-pattern-morph frequencies such that between-host migrants are more likely to be the locally less-cryptic morph than are residents. Collectively, these data indicate that high rates of visual predation on less-cryptic migrants are likely to reduce encounters, and thus interbreeding, between host-associated populations. Comparison with previous estimates of sexual isolation reveals that the contribution of selection against between-host migrants to total premating isolation is as strong as, or stronger than, that of sexual isolation. These findings highlight the potential role of natural selection against migrants between divergent environments in the formation of new species.

The evolution of müllerian mimicry in multispecies communities

Prey species that are unprofitable to attack often share conspicuous colours and patterns with other coexisting defended species. This phenomenon, termed müllerian mimicry, has long been explained as a consequence of selection on defended prey to adopt a common way of advertising their unprofitability. However, studies using two unpalatable prey types have not always supported this theory. Here we show, using a system of humans hunting for computer-generated prey, that predators do not always generate strong selection for mimicry when there are two unprofitable prey types. By contrast, we demonstrate that when predators are faced with a range of different prey species, selection on unprofitable prey to resemble one another can be intense. Here the primary selective force is not one in which predators evaluate the profitabilities of distinct prey types independently, but one in which predators learn better to avoid unprofitable phenotypes that share traits distinguishing them from profitable prey. This need to simplify decision making readily facilitates the spread of imperfect mimetic forms from rarity, and suggests that müllerian mimicry is more likely to arise in multispecies communities.

Historical demography of Mullerian mimicry in the neotropical Heliconius butterflies

We compare the historical demographies of two Müllerian comimetic butterfly species: Heliconius erato and Heliconius melpomene. These species show an extensive parallel geographic divergence in their aposematic wing phenotypes. Recent studies suggest that this coincident mosaic results from simultaneous demographic processes shaped by extrinsic forces over Pleistocene climate fluctuations. However, DNA sequence variation at two rapidly evolving unlinked nuclear loci, Mannose phosphate isomerase (Mpi) and Triose phosphate isomerase (Tpi), show that the comimetic species have quite different quaternary demographies. In H. erato, despite ongoing lineage sorting across the Andes, nuclear genealogical estimates showed little geographical structure, suggesting high historical gene flow. Coalescent-based demographic analysis revealed population growth since the Pliocene period. Although these patterns suggest vicariant population subdivision associated with the Andean orogeny, they are not consistent with hypotheses of Pleistocene population fragmentation facilitating allopatric wing phenotype radiation in H. erato. In contrast, nuclear genetic diversity, theta, in H. melpomene was reduced relative to its comimic and revealed three phylogeographical clades. The pattern of coalescent events within regional clades was most consistent with population growth in relatively isolated populations after a recent period of restricted population size. These different demographic histories suggest that the wing-pattern radiations were not coincident in the two species. Instead, larger effective population size (N(e)) in H. erato, together with profound population change in H. melpomene, supports an earlier hypothesis that H. erato diversified first as the model species of this remarkable mimetic association.

The importance of pattern similarity between Müllerian mimics in predator avoidance learning

Müllerian mimicry, where unpalatable prey share common warning patterns, has long fascinated evolutionary biologists. It is commonly assumed that Müllerian mimics benefit by sharing the costs of predator education, thus reducing per capita mortality, although there has been no direct test of this assumption. Here, we specifically measure the selection pressure exerted by avian predators on unpalatable prey with different degrees of visual similarity in their warning patterns. Using wild-caught birds foraging on novel patterned prey in the laboratory, we unexpectedly found that pattern similarity did not increase the speed of avoidance learning, and even dissimilar mimics shared the education of naive predators. This was a consistent finding across two different densities of unpalatable prey, although mortalities were lower at the higher density as expected. Interestingly, the mortalities of Müllerian mimics were affected by pattern similarity in the predicted way by the end of our experiment, although the result was not quite significant. This suggests that the benefits to Müllerian mimics may emerge only later in the learning process, and that predator experience of the patterns may affect the degree to which pattern similarity is important. This highlights the need to measure the behaviour of real predators if we are to understand fully the evolution of mimicry systems.

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.

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.

Predators favour mimicry in a tropical reef fish

Batesian mimicry evolves when the 'umbrella' of protection provided by resemblance to a conspicuous unpalatable model species is sufficient to overcome increased predation risk associated with greater conspicuousness. However, the shape and extent of this umbrella, that is, how the level of protection provided by mimicry changes with degree of resemblance between model and mimic, is poorly known. We investigated the response of wild predatory fishes to plastic replicas of a model-mimic species pair of tropical reef fishes, Canthigaster valentini (a toxic pufferfish, the model) and Paraluteres prionurus (the putative mimic), and additional replicas with progressively lower degrees of resemblance to the mimic species. Our results reveal a relatively broad region of protection, indicated by a reduced approach rate by piscivorous fishes, surrounding the colour pattern of the model species. Protection increased with increasing resemblance. By contrast, the response of non-piscivorous fishes was unrelated to degree of resemblance of replicas to the model. Our results suggest that piscivorous fishes on the reef are educated regarding the toxicity of C. valentini, and that avoidance of fish having the pufferfish colour pattern has generated selection favouring mimetic resemblance by the palatable P. prionurus. The relatively broad protective umbrella has probably facilitated the initial evolution of resemblance in the palatable prey species despite the potential hazards of greater conspicuousness.

Phylogenetic discordance at the species boundary: comparative gene genealogies among rapidly radiating Heliconius butterflies

Recent adaptive radiations provide excellent model systems for understanding speciation, but rapid diversification can cause problems for phylogenetic inference. Here we use gene genealogies to investigate the phylogeny of recent speciation in the heliconiine butterflies. We sequenced three gene regions, intron 3 ( approximately 550 bp) of sex-linked triose-phosphate isomerase (Tpi), intron 3 ( approximately 450 bp) of autosomal mannose-phosphate isomerase (Mpi), and 1,603 bp of mitochondrial cytochrome oxidase subunits I and II (COI and COII), for 37 individuals from 25 species of Heliconius and related genera. The nuclear intron sequences evolved at rates similar to those of mitochondrial coding sequences, but the phylogenetic utility of introns was restricted to closely related geographic populations and species due to high levels of indel variation. For two sister species pairs, Heliconius erato-Heliconius himera and Heliconius melpomene-Heliconius cydno, there was highly significant discordance between the three genes. At mtDNA and Tpi, the hypotheses of reciprocal monophyly and paraphyly of at least one species with respect to its sister could not be distinguished. In contrast alleles sampled from the third locus, Mpi, showed polyphyletic relationships between both species pairs. In all cases, recent coalescence of mtDNA lineages within species suggests that polyphyly of nuclear genes is not unexpected. In addition, very similar alleles were shared between melpomene and cydno, implying recent gene flow. Our finding of discordant genealogies between genes is consistent with models of adaptive speciation with ongoing gene flow and highlights the need for multiple locus comparisons to resolve phylogeny among closely related species.

Correlational selection and the evolution of genomic architecture

We review and discuss the importance of correlational selection (selection for optimal character combinations) in natural populations. If two or more traits subject to multivariate selection are heritable, correlational selection builds favourable genetic correlations through the formation of linkage disequilibrium at underlying loci governing the traits. However, linkage disequilibria built up by correlational selection are expected to decay rapidly (ie, within a few generations), unless correlational selection is strong and chronic. We argue that frequency-dependent biotic interactions that have 'Red Queen dynamics' (eg, host-parasite interactions, predator-prey relationships or intraspecific arms races) often fuel chronic correlational selection, which is strong enough to maintain adaptive genetic correlations of the kind we describe. We illustrate these processes and phenomena using empirical examples from various plant and animal systems, including our own recent work on the evolutionary dynamics of a heritable throat colour polymorphism in the side-blotched lizard Uta stansburiana. In particular, male and female colour morphs of side-blotched lizards cycle on five- and two-generation (year) timescales under the force of strong frequency-dependent selection. Each morph refines the other morph in a Red Queen dynamic. Strong correlational selection gradients among life history, immunological and morphological traits shape the genetic correlations of the side-blotched lizard polymorphism. We discuss the broader evolutionary consequences of the buildup of co-adapted trait complexes within species, such as the implications for speciation processes.

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.

The genetics and evo-devo of butterfly wing patterns

Understanding how the spectacular diversity of colour patterns on butterfly wings is shaped by natural selection, and how particular pattern elements are generated, has been the focus of both evolutionary and developmental biologists. The growing field of evolutionary developmental biology has now begun to provide a link between genetic variation and the phenotypes that are produced by developmental processes and that are sorted by natural selection. Butterfly wing patterns are set to become one of the few examples of morphological diversity to be studied successfully at many levels of biological organization, and thus to yield a more complete picture of adaptive morphological evolution.

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

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

Reproductive isolation caused by colour pattern mimicry

Speciation is facilitated if ecological adaptation directly causes assortative mating, but few natural examples are known. Here we show that a shift in colour pattern mimicry was crucial in the origin of two butterfly species. The sister species Heliconius melpomene and Heliconius cydno recently diverged to mimic different model taxa, and our experiments show that their mimetic coloration is also important in choosing mates. Assortative mating between the sister species means that hybridization is rare in nature, and the few hybrids that are produced are non-mimetic, poorly adapted intermediates. Thus, the mimetic shift has caused both pre-mating and post-mating isolation. In addition, individuals from a population of H. melpomene allopatric to H. cydno court and mate with H. cydno more readily than those from a sympatric population. This suggests that assortative mating has been enhanced in sympatry.

Three-butterfly system provides a field test of müllerian mimicry

In 1879, Müller proposed that two brightly coloured distasteful butterfly species (co-models) that share a single warning-colour pattern would benefit by spreading the selective burden of educating predators. The mutual benefit of sharing warning signals among distasteful species, so-called müllerian mimicry, is supported by comparative evidence, theoretical studies and laboratory simulations; however, to date, this key exemplar of adaptive evolution has not been experimentally tested in the field. To measure natural selection generated by müllerian mimicry, I exploited the unusual polymorphism of Heliconius cydno (Lepidoptera: Nymphalidae). Here I show increased survival of H. cydno morphs that match locally abundant monomorphic co-model species. This study demonstrates müllerian mimicry in the field. It also shows that müllerian mimicry with several co-models generates geographically divergent selection, which explains the existence of polymorphism in distasteful species with warning coloration.

What initiates speciation in passion-vine butterflies?

Studies of the continuum between geographic races and species provide the clearest insights into the causes of speciation. Here we report on mate choice and hybrid viability experiments in a pair of warningly colored butterflies, Heliconius erato and Heliconius himera, that maintain their genetic integrity in the face of hybridization. Hybrid sterility and inviability have been unimportant in the early stages of speciation of these two Heliconius. We find no evidence of reduced fecundity, egg hatch, or larval survival nor increases in developmental time in three generations of hybrid crosses. Instead, speciation in this pair appears to have been catalyzed by the association of strong mating preferences with divergence in warning coloration and ecology. In mate choice experiments, matings between the two species are a tenth as likely as matings within species. F1 hybrids of both sexes mate frequently with both pure forms. However, male F1 progeny from crosses between H. himera mothers and H. erato fathers have somewhat reduced mating success. The strong barrier to gene flow provided by divergence in mate preference is probably enhanced by frequency-dependent predation against hybrids similar to the type known to occur across interracial hybrid zones of H. erato. In addition, the transition between this pair falls at the boundary between wet and dry forest, and rare hybrids may also be selected against because they are poorly adapted to either biotope. These results add to a growing body of evidence that challenge the importance of genomic incompatibilities in the earliest stages of speciation.

Locomotor mimicry in butterflies? The associations of positions of centres of mass among groups of mimetic, unprofitable prey

With detailed measurements of flight-related morphological parameters of 18 species within the tribe Heliconiini and 10 of their non-heliconiine comimics, I found that morphological parameters relevant to flight biomechanics are associated with three escape tactics of these Neotropical butterflies: evasive flight, distastefulness, and mimicry. Two distinct character suites, one pertaining to the position of centre of body mass and another pertaining to wing shape, were identified using principal-components analysis. A third component correlated most highly with the position of centre of wing mass. Mimicry groups were best discriminated by the positions of centres of body and wing mass, suggesting that these features converged within mimicry groups. Centres of mass were positioned nearer to the wing base in mimicry groups composed of more palatable species, presumably decreasing the radial moments of inertia of the body and wings and increasing flight speed and turning performance. Predation has selected for body morphology that increases flight speed and manoeuvrability in palatable butterflies, whereas the morphology of distasteful species compromises flight performance. Convergent selection may reduce morphological differences of species within mimicry groups arising from distantly related lineages, and hence it may have enhanced the morphological diversification of palatable and unpalatable butterflies within lineages that have more recently evolved distastefulness.

Comparisons among morphological characters and between localities in the Chorthippus parallelus hybrid zone (Orthoptera: Acrididae)

Two subspecies of the grasshopper, Chorthippus parallelus , meet and hybridize in the Pyrenees. The hybrid zone between the two taxa is believed to have formed following range expansion at the end of the last glaciation and to be maintained by a balance between gene flow and selection against hybrids. Laboratory FI hybrid males are sterile. We have examined morphological characters in two cols about 200 km apart and compared the positions and widths of dines for these characters both within and between cols. Clines are neither coincident (centred in the same position), nor concordant (equal in width) in either col. Furthermore clines are wider and more dispersed in the western col. The lack of concordance can be explained on several models for the maintenance of the hybrid zone and is a common observation. There are many fewer examples of lack of coincidence and this result is more difficult to interpret. We suggest that it is due to a difference in population structure between the cols: a patchy distribution of grasshoppers in the western col produces wider, more dispersed dines while a sharp break in distribution in the eastern col tends to concentrate morphological change over a shorter distance.

Did forest islands drive the diversity of warningly coloured butterflies? Biotic drift and the shifting balance

Species of the South American butterfly genus Heliconius have undergone remarkably wide racial divergence in their patterns, and most of the resulting races are muellerian mimics. As warning coloration normally imposes stabilizing selection on the pattern, this divergence is much in need of explanation. Two models have been suggested. Brown, Sheppard and Turner proposed that the divergence results from ‘mimetic capture’, the switching of patterns between adaptive peaks generated by changes in the overall composition of the local biota (‘biotic drift’) and hence of the mimicry rings to which each species belongs; these changes have in turn been generated by long term patterns of species extinction in island refuges as biota became progressively isolated and contiguous during contraction and expansion of the rain forest during the Pleistocene. An alternative model, proposed by Mallet, is that truly novel colour patterns became established by mutation and random drift, then spreading to become predominant in local areas; subsequently the novel patterns spread over wide areas by the migration of dines. Under this application of Wright’s shifting balance model, refuges are not necessary for divergence, and muellerian mimicry evolves after divergence rather than being the driving force for race formation. Although our respective models appear diametrically opposed, the hypotheses are difficult to distinguish and there are broad areas of agreement; in both models there is an initial stochastic event, followed by natural selection for mimicry, and both will operate either in parapatry or allopatry. The diversity of warning patterns is better explained by the shifting balance model, but there are alternative selectionist explanations such as sexual selection.

Adaptation, speciation and hybrid zones

Many species are divided into a mosaic of genetically distinct populations, separated by narrow zones of hybridization. Studies of hybrid zones allow us to quantify the genetic differences responsible for speciation, to measure the diffusion of genes between diverging taxa, and to understand the spread of alternative adaptations.