Species specificity and intraspecific variation in the chemical profiles of Heliconius butterflies across a large geographic range

In many animals, mate choice is important for the maintenance of reproductive isolation between species. Traits important for mate choice and behavioral isolation are predicted to be under strong stabilizing selection within species; however, such traits can also exhibit variation at the population level driven by neutral and adaptive evolutionary processes. Here, we describe patterns of divergence among androconial and genital chemical profiles at inter‐ and intraspecific levels in mimetic Heliconius butterflies. Most variation in chemical bouquets was found between species, but there were also quantitative differences at the population level. We found a strong correlation between interspecific chemical and genetic divergence, but this correlation varied in intraspecific comparisons. We identified “indicator” compounds characteristic of particular species that included compounds already known to elicit a behavioral response, suggesting an approach for identification of candidate compounds for future behavioral studies in novel systems. Overall, the strong signal of species identity suggests a role for these compounds in species recognition, but with additional potentially neutral variation at the population level.     

Quantified reproductive isolation in Heliconius butterflies: Implications for introgression and hybrid speciation

Heliconius butterflies have become a model for the study of speciation with gene flow. For adaptive introgression to take place, there must be incomplete barriers to gene exchange that allow interspecific hybridization and multiple generations of backcrossing. The recent publication of estimates of individual components of reproductive isolation between several species of butterflies in the Heliconius melpomene–H. cydno clade allowed us to calculate total reproductive isolation estimates for these species. According to these estimates, the butterflies are not as promiscuous as has been implied. Differences between species are maintained by intrinsic mechanisms, while reproductive isolation of geographical races within species is mainly due to allopatry. We discuss the implications of this strong isolation for basic aspects of the hybrid speciation with introgression hypothesis.     

Reinforcement of mate preference among hybridizing Heliconius butterflies

Recent models of mate preference evolution suggest that direct selection on alleles at preference loci and correlated evolution of preference with locally adapted mating cues are more likely to drive the evolution of assortative mate preference than reinforcement. Mate preference evolution in mimetic Heliconius butterflies has been attributed to all three forms of selection, but here we show that reinforcement has been critical. By examining geographical variation in assortative mating and male mate preference among seven populations of three hybridizing Heliconius species from Costa Rica, we found pronounced character displacement of preference such that sexual isolation was enhanced in areas of interspecific contact. Of the different explanations for the evolution of assortative mate preference, only reinforcement is dependent on interspecific contact in this system. Thus, the observed pattern of reproductive character displacement of mate preference is best explained as a product of indirect selection generated by natural selection against nonmimetic hybrids.     

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.     

IDENTIFICATION OF GENETICALLY LINKED FEMALE PREFERENCE AND MALE TRAIT

Genetic variation in male traits and the female preferences for those traits allows for the evolution of sexual behavior. Trait–preference combinations are thought to improve the effectiveness of runaway sexual selection within a species, and are considered necessary for the induction of divergence between species. Novel traits, or variants of existing traits, and their associated preferences in the opposite sex are more likely to be maintained if they are genetically linked in proximity on a chromosome (the genetic coupling hypothesis), yet there is little empirical evidence that this genetic linkage occurs. Here we show for the first time that natural genetic variation at a single‐linked region can induce both species‐specific female choosiness and the male trait they are discriminating against. We found this effect in two separate regions of the genome, demonstrating that this linkage may be common. In contrast, female choosiness and male unattractiveness could not be alleviated by a single region. The close linkage of these loci and the strength of their effect provide an evolutionary means by which this preference–trait combination could arise and be maintained, thus enabling a more rapid route for runaway sexual selection, and providing empirical evidence supporting the genetic coupling hypothesis.     

Limits to runaway sexual selection: The wallflower paradox

In his mathematical treatment of Fisher's ideas on sexual selection (so-called runaway selection) Lande (1981) predicted that males may evolve increasingly elaborate sexual characters despite opposing viability selection as a consequence of the associated costs. Lande thereby assumed that female mate preferences are not subject to selection since (1) females are all inseminated and (2) the quantity and quality of their offspring are independent of the female's mate preferences. Kirkpatrick (1985) removed the latter assumption and investigated the consequences for the mean phenotype with respect to both female and male traits. He also explored the dynamics of the (co)-variance matrix by numerical methods. In this paper we consider a simpler model with just two multi-allelic loci. This enables us to derive explicit expressions for (co)-variances under steady state conditions. Rather than assume natural selection through differential fertility (as in Kirkpatrick, 1985), we take sexual selection on females into account by modelling the preference-dependent risk that females remain unmated. We argue that this wallflower effect is a realistic feature of any mating system, since it merely depends on the existence of (1) variation in mating preferences and (2) a finite mating season. Our approach provided an insight into the dynamic behaviour of the means of the phenotypes. This is because the dynamics of the means depend on the steady state (co)-variance matrix. Thus, an insight into the former requires explicit expressions for the latter. Whereas Lande and Kirkpatrick predicted runaway processes, despite opposing viability selection, our model predicts a globally stable steady state, i.e. no runaway, even without opposing viability selection (under the assumption of an asymptotically stable steady state of the (co)-variances. Admittedly, we have no analytic proof of this stability but only support for it, based on simulations.) The absence of the runaway processes in our model is caused by the wallflower effect, since it imposes constraints on the steady state of the (co)-variance matrix. When mutational input applies to female traits but not to male traits, explicit expressions for the (co)-variances under steady state conditions can be derived, and these show that: (1) both the genetic covariance and the variance of male traits are equal to zero, but (2) the variance of the female trait exceeds to zero. Should there be mutational input influencing the male trait, then these results would suggest that the male-to-female ratio of variances is much smaller than unity. This prediction is of tremendous importance for speciation through founding events.     

Male mate choice via cuticular hydrocarbon pheromones drives reproductive isolation between Drosophila species

Mate discrimination is a key mechanism restricting gene flow between species. While studied extensively with respect to female mate choice, mechanisms of male mate choice between species are far less studied. Thus, we have little knowledge of the relative frequency, importance, or overall contribution of male mate discrimination to reproductive isolation. In the present study, we estimated the relative contributions of male and female choice to reproductive isolation between Drosophila simulans and D. sechellia, and show that male mate discrimination accounts for the majority of the current isolation between these species. We further demonstrate that males discriminate based on female cuticular hydrocarbon pheromones, and collect evidence supporting the hypothesis that male mate discrimination may alleviate the costs associated with heterospecific courtship and mating. Our findings highlight the potentially significant contribution of male mate choice to the formation of reproductive isolating barriers, and thus the speciation process.     

Testing for mating isolation between ecotypes: laboratory experiments with lake, stream and hybrid stickleback

Mating isolation is a frequent contributor to ecological speciation – but how consistently does it evolve as a result of divergent selection? We tested for genetically based mating isolation between lake and stream threespine stickleback (Gasterosteus aculeatus L.) from the Misty watershed, Vancouver Island, British Columbia. We combined several design elements that are uncommon in the studies of stickleback mate choice: (i) we used second‐generation laboratory‐reared fish (to reduce environmental and maternal effects), (ii) we allowed for male–male competitive interactions (instead of the typical no‐choice trials) and (iii) we included hybrids along with pure types. Males of different types (Lake, Inlet, hybrid) were paired in aquaria, allowed to build nests and then exposed sequentially to females of all three types. We found that Lake and Inlet males differed in behaviours thought to influence stickleback mate choice (inter‐ and intra‐sexual aggression, display and nest activities), whereas hybrids were either intermediate or apparently ‘inferior’ in these behaviours. Despite these differences, Lake and Inlet fish did not mate assortatively and hybrid males did not have a mating disadvantage. Our study reinforces the noninevitability of mating isolation evolving in response to ecological differences and highlights the need to further investigate the factors promoting and constraining progress towards ecological speciation.     

Pervasive genetic associations between traits causing reproductive isolation in Heliconius butterflies

Ecological speciation proceeds through the accumulation of divergent traits that contribute to reproductive isolation, but in the face of gene flow traits that characterize incipient species may become disassociated through recombination. Heliconius butterflies are well known for bright mimetic warning patterns that are also used in mate recognition and cause both pre- and post-mating isolation between divergent taxa. Sympatric sister taxa representing the final stages of speciation, such as Heliconius cydno and Heliconius melpomene, also differ in ecology and hybrid fertility. We examine mate preference and sterility among offspring of crosses between these species and demonstrate the clustering of Mendelian colour pattern loci and behavioural loci that contribute to reproductive isolation. In particular, male preference for red patterns is associated with the locus responsible for the red forewing band. Two further colour pattern loci are associated, respectively, with female mating outcome and hybrid sterility. This genetic architecture in which ‘speciation genes’ are clustered in the genome can facilitate two controversial models of speciation, namely divergence in the face of gene flow and hybrid speciation.     

Cultural isolation is greater than genetic isolation across an avian hybrid zone

Elucidating the relationship between genetic and cultural evolution is important in understanding speciation, as learned premating barriers might be involved in maintaining species differences. Here, we test this relationship by examining a widely recognized premating barrier, bird song, in a hybrid zone between black‐throated green (Setophaga virens) and Townsend's warblers (S. townsendi). We use song analysis, genomic techniques and playback experiments to characterize the cultural and genetic backgrounds of individuals in this region, expecting that if song is an important reproductive barrier between these species, there should be a strong relationship between song and genotype. We show that songs in the hybrid zone correspond to the distinctly different songs found in allopatry but that song and genotype are not tightly coupled in sympatry. Allopatric individuals responded only to local songs, indicating that individuals may have learned to respond to songs they commonly hear. We observed discordance between song and genotype clines; a narrower cline suggests that cultural selection on song is stronger than natural selection on genotype. These findings indicate that song is unlikely to play a role in reproductive isolation between these species, and we suggest that spatial variation in song may nonetheless be maintained by frequency‐dependent cultural selection. This decoupling of genes and culture may contribute to hybridization in this region.     

A ROLE FOR LEARNING IN POPULATION DIVERGENCE OF MATE PREFERENCES

Learning and other forms of phenotypic plasticity have been suggested to enhance population divergence. Mate preferences can develop by learning, and species recognition might not be entirely genetic. We present data on female mate preferences of the banded demoiselle (Calopteryx splendens) that suggest a role for learning in population divergence and species recognition. Populations of this species are either allopatric or sympatric with a phenotypically similar congener (C. virgo). These two species differ mainly in the amount of wing melanization in males, and wing patches thus mediate sexual isolation. In sympatry, sexually experienced females discriminate against large melanin wing patches in heterospecific males. In contrast, in allopatric populations within the same geographic region, females show positive (“open‐ended”) preferences for such large wing patches. Virgin C. splendens females do not discriminate against heterospecific males. Moreover, physical exposure experiments of such virgin females to con‐ or hetero‐specific males significantly influences their subsequent mate preferences. Species recognition is thus not entirely genetic and it is partly influenced by interactions with mates. Learning causes pronounced population divergence in mate preferences between these weakly genetically differentiated populations, and results in a highly divergent pattern of species recognition at a small geographic scale.     

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.     

Ecological and genetic factors influencing the transition between host‐use strategies in sympatric Heliconius butterflies

Shifts in host‐plant use by phytophagous insects have played a central role in their diversification. Evolving host‐use strategies will reflect a trade‐off between selection pressures. The ecological niche of herbivorous insects is partitioned along several dimensions, and if populations remain in contact, recombination will break down associations between relevant loci. As such, genetic architecture can profoundly affect the coordinated divergence of traits and subsequently the ability to exploit novel habitats. The closely related species Heliconius cydno and H. melpomene differ in mimetic colour pattern, habitat and host‐plant use. We investigate the selection pressures and genetic basis underlying host‐use differences in these two species. Host‐plant surveys reveal that H. melpomene specializes on a single species of Passiflora. This is also true for the majority of other Heliconius species in secondary growth forest at our study site, as expected under a model of interspecific competition. In contrast, H. cydno, which uses closed‐forest habitats where both Heliconius and Passiflora are less common, appears not to be restricted by competition and uses a broad selection of the available Passiflora. However, other selection pressures are likely involved, and field experiments reveal that early larval survival of both butterfly species is highest on Passiflora menispermifolia, but most markedly so for H. melpomene, the specialist on that host. Finally, we demonstrate an association between host‐plant acceptance and colour pattern amongst interspecific hybrids, suggesting that major loci underlying these important ecological traits are physically linked in the genome. Together, our results reveal ecological and genetic associations between shifts in habitat, host use and mimetic colour pattern that have likely facilitated both speciation and coexistence.     

Geographic contrasts between pre‐ and postzygotic barriers are consistent with reinforcement in Heliconius butterflies

Identifying the traits causing reproductive isolation and the order in which they evolve is fundamental to understanding speciation. Here, we quantify prezygotic and intrinsic postzygotic isolation among allopatric, parapatric, and sympatric populations of the butterflies Heliconius elevatus and Heliconius pardalinus. Sympatric populations from the Amazon (H. elevatus and H. p. butleri) exhibit strong prezygotic isolation and rarely mate in captivity; however, hybrids are fertile. Allopatric populations from the Amazon (H. p. butleri) and Andes (H. p. sergestus) mate freely when brought together in captivity, but the female F1 hybrids are sterile. Parapatric populations (H. elevatus and H. p. sergestus) exhibit both assortative mating and sterility of female F1s. Assortative mating in sympatric populations is consistent with reinforcement in the face of gene flow, where the driving force, selection against hybrids, is due to disruption of mimicry and other ecological traits rather than hybrid sterility. In contrast, the lack of assortative mating and hybrid sterility observed in allopatric populations suggests that geographic isolation enables the evolution of intrinsic postzygotic reproductive isolation. Our results show how the types of reproductive barriers that evolve between species may depend on geography.     

Divergent patterns of diversification in courtship and genitalic characters of Timema walking‐sticks

Understanding the patterns of diversification in sexual traits and the selection underlying such diversification represents a major unresolved question in evolutionary biology. We examined the phylogenetic diversification for courtship and external genitalic characters across ten species of Timema walking‐sticks, to infer the tempos and modes of character change in these sexual traits and to draw inferences regarding the selective pressures underlying speciation and diversification in this clade. Rates of inferred change in male courtship behaviours were proportional to speciation events, but male external genitalic structures showed a pattern of continuous change across evolutionary time, with divergence proportional to branch lengths. These findings suggest that diversification of courtship behaviour is mediated by processes that occur in association with speciation, whereas diversification of genitalia occurs more or less continuously, most likely driven by forces of sexual selection.     

Independent inheritance of preference and performance in hybrids between host races of Mitoura butterflies (Lepidoptera: Lycaenidae)

Divergent natural selection contributes to reproductive isolation among populations adapting to different habitats or resources if hybrids between populations are intermediate in phenotype and suffer an associated, environmentally dependent reduction in fitness. This prediction was tested using two host races of Mitoura butterflies. Thirty-five F1 hybrid and parental lines were created, larvae were raised on the two host plants, and oviposition preferences were assayed in choice arenas. Larvae from both reciprocal hybrid crosses suffered a host-specific reduction in performance: when reared on incense cedar, hybrid survival was approximately 30% less than the survival of pure lines of the cedar-associated host race. The performance of hybrid larvae reared on the other host, MacNab cypress, was not reduced relative to parental genotypes. Females from both reciprocal hybrid crosses preferred to oviposit on incense cedar, the same host that resulted in the reduced survival of hybrid larvae. Thus, dominance is implicated in the inheritance of traits involved in both preference and performance, which do not appear to be genetically linked in Mitoura butterflies. Gene flow between host races may be reduced because the correlation between preference and performance that was previously described in parental populations is essentially broken by hybridization.     

Differentiation in putative male sex pheromone components across and within populations of the African butterfly Bicyclus anynana as a potential driver of reproductive isolation

Sexual traits are often the most divergent characters among closely related species, suggesting an important role of sexual traits in speciation. However, to prove this, we need to show that sexual trait differences accumulate before or during the speciation process, rather than being a consequence of it. Here, we contrast patterns of divergence among putative male sex pheromone (pMSP) composition and the genetic structure inferred from variation in the mitochondrial cytochrome oxidase 1 and nuclear CAD loci in the African butterfly Bicyclus anynana (Butler, 1879) to determine whether the evolution of “pheromonal dialects” occurs before or after the differentiation process. We observed differences in abundance of some shared pMSP components as well as differences in the composition of the pMSP among B. anynana populations. In addition, B. anynana individuals from Kenya displayed differences in the pMSP composition within a single population that appeared not associated with genetic differences. These differences in pMSP composition both between and within B. anynana populations were as large as those found between different Bicyclus species. Our results suggest that “pheromonal dialects” evolved within and among populations of B. anynana and may therefore act as precursors of an ongoing speciation process.     

Secondary contact between Lycaeides idas and L. melissa in the Rocky Mountains: extensive admixture and a patchy hybrid zone

Studies of hybridization have increased our understanding of the nature of species boundaries, the process of speciation, and the effects of hybridization on the evolution of populations and species. In the present study we use genetic and morphological data to determine the outcome and consequences of secondary contact and hybridization between the butterfly species Lycaeides idas and L. melissa in the Rocky Mountains. Admixture proportions estimated from structure and geographical cline analysis indicate L. idas and L. melissa have hybridized extensively in the Rocky Mountains and that reproductive isolation was insufficient to prevent introgression for much of the genome. Geographical patterns of admixture suggest that hybridization between L. idas and L. melissa has led to the formation of a hybrid zone. The hybrid zone is relatively wide, given estimates of dispersal for Lycaeides butterflies, and does not show strong evidence of cline concordance among characters. We believe the structure of the Lycaeides hybrid zone might be best explained by the patchy distribution of Lycaeides, local extinction and colonization of habitat patches, environmental variation and weak overall selection against hybrids. We found no evidence that hybridization in the Rocky Mountains has resulted in the formation of independent hybrid species, in contrast to the outcome of hybridization between L. idas and L. melissa in the Sierra Nevada. Finally, our results suggest that differences in male morphology between L. idas and L. melissa might contribute to isolation, or perhaps even that selection has favoured the spread of L. melissa male genitalia alleles.     

Prenuptial agreements: mating frequency predicts gift-giving in Heliconius species

Theory predicts that when males provision females with nuptial gifts that include nutrients, the degree of polyandry should be positively correlated with the size or quality of the gift. This is because larger and more nutritious gifts tend to increase female refractory period, reducing the chances the female will remate soon. This decreases the likelihood of sperm competiton and consequently increases the donor male fitness. Butterflies in the genus Heliconius Kluk (Lepidoptera: Nymphalidae: Heliconini) exhibit variable mating systems that include monandry and polyandry. In addition to protein in the spermatophore, males increase gift quality by providing females with cyanide, which may contribute to protection of the female or her eggs. We tested whether degree of polyandry and gift quality (spermatophore weight and cyanide content) were correlated in nine Heliconius species from greenhouse populations. As predicted, both spermatophore weight and cyanide content were correlated with mating frequency. This is the first report to show that degree of polyandry correlates with allocation of defensive chemical as part of a nuptial gift.