Hybrid incompatibility is consistent with a hybrid origin of Heliconius heurippa Hewitson from its close relatives, Heliconius cydno Doubleday and Heliconius melpomene Linnaeus

100 years of Haldane's rule

Haldane's rule is one of the 'two rules of speciation'. It states that if one sex is 'absent, rare or sterile' in a hybrid population, then that sex will be heterogametic. Since Haldane first made this observation, 100 years have passed and still questions arise over how many independent examples exist and what the underlying causes of Haldane's rule are. This review aims to examine research that has occurred over the last century. It seeks to do so by discussing possible causes of Haldane's rule, as well as gaps in the research of these causes that could be readily addressed today. After 100 years of research, it can be concluded that Haldane's rule is a complicated one, and much current knowledge has been accrued by studying the model organisms of speciation. This has led to the primacy of dominance theory and faster-male theory as explanations for Haldane's rule. However, some of the most interesting findings of the 21st century with regard to Haldane's rule have involved investigating a wider range of taxa emphasizing the need to continue using comparative methods, including ever more taxa as new cases are discovered.

Complex basis of hybrid female sterility and Haldane's rule in Heliconius butterflies: Z-linkage and epistasis

Hybrids between species are often sterile or inviable. Hybrid unfitness usually evolves first in the heterogametic sex-a pattern known as Haldane's rule. The genetics of Haldane's rule have been extensively studied in species where the male is the heterogametic (XX/XY) sex, but its basis in taxa where the female is heterogametic (ZW/ZZ), such as Lepidoptera and birds, is largely unknown. Here, we analyse a new case of female hybrid sterility between geographic subspecies of Heliconius pardalinus. The two subspecies mate freely in captivity, but female F1 hybrids in both directions of cross are sterile. Sterility is due to arrested development of oocytes after they become differentiated from nurse cells, but before yolk deposition. We backcrossed fertile male F1 hybrids to parental females and mapped quantitative trait loci (QTLs) for female sterility. We also identified genes differentially expressed in the ovary as a function of oocyte development. The Z chromosome has a major effect, similar to the 'large X effect' in Drosophila, with strong epistatic interactions between loci at either end of the Z chromosome, and between the Z chromosome and autosomal loci on chromosomes 8 and 20. By intersecting the list of genes within these QTLs with those differentially expressed in sterile and fertile hybrids, we identified three candidate genes with relevant phenotypes. This study is the first to characterize hybrid sterility using genome mapping in the Lepidoptera and shows that it is produced by multiple complex epistatic interactions often involving the sex chromosome, as predicted by the dominance theory of Haldane's rule.

The role of climate change in pollinator decline across the Northern Hemisphere is underestimated

Pollinator biodiversity loss occurs at unprecedented rates globally, with particularly sharp declines documented in the North Temperate Zone. There is currently no consensus on the main drivers of the decline. Although climate change is expected to drive biodiversity loss in the future, current warming is often suggested to have positive impacts on pollinator assemblages in higher latitudes. Consequently, pollinator conservation initiatives in Europe and the USA tend to lack climate adaptation initiatives, an omission of which may be risky if climate change has significant negative impacts on pollinators. To gain an understanding of the impacts of climate change on pollinator biodiversity in the Northern Hemisphere, we conducted a literature review on genetic, species and community level diversity. Our findings suggest that global heating most likely causes homogenization of pollinator assemblages at all levels of pollinator biodiversity, making them less resilient to future stochasticity. Aspects of biodiversity that are rarely measured (e.g. genetic diversity, β-diversity, species evenness) tend to be most affected, while some dimensions of climate change, such as fluctuations in winter weather conditions, changes in the length of the vegetational season and increased frequency of extreme weather events, that seldom receive attention in empirical studies, tend to be particularly detrimental to pollinators. Negative effects of global heating on pollinator biodiversity are most likely exacerbated by homogenous and fragmented landscapes, widespread across Europe and the US, which limit opportunities for range-shifts and reduce micro-climatic buffering. This suggests the need for conservation initiatives to focus on increasing landscape connectivity and heterogeneity at multiple spatial scales.

Diversity in warning coloration: selective paradox or the norm?

Aposematic theory has historically predicted that predators should select for warning signals to converge on a single form, as a result of frequency-dependent learning. However, widespread variation in warning signals is observed across closely related species, populations and, most problematically for evolutionary biologists, among individuals in the same population. Recent research has yielded an increased awareness of this diversity, challenging the paradigm of signal monomorphy in aposematic animals. Here we provide a comprehensive synthesis of these disparate lines of investigation, identifying within them three broad classes of explanation for variation in aposematic warning signals: genetic mechanisms, differences among predators and predator behaviour, and alternative selection pressures upon the signal. The mechanisms producing warning coloration are also important. Detailed studies of the genetic basis of warning signals in some species, most notably Heliconius butterflies, are beginning to shed light on the genetic architecture facilitating or limiting key processes such as the evolution and maintenance of polymorphisms, hybridisation, and speciation. Work on predator behaviour is changing our perception of the predator community as a single homogenous selective agent, emphasising the dynamic nature of predator-prey interactions. Predator variability in a range of factors (e.g. perceptual abilities, tolerance to chemical defences, and individual motivation), suggests that the role of predators is more complicated than previously appreciated. With complex selection regimes at work, polytypisms and polymorphisms may even occur in Müllerian mimicry systems. Meanwhile, phenotypes are often multifunctional, and thus subject to additional biotic and abiotic selection pressures. Some of these selective pressures, primarily sexual selection and thermoregulation, have received considerable attention, while others, such as disease risk and parental effects, offer promising avenues to explore. As well as reviewing the existing evidence from both empirical studies and theoretical modelling, we highlight hypotheses that could benefit from further investigation in aposematic species. Finally by collating known instances of variation in warning signals, we provide a valuable resource for understanding the taxonomic spread of diversity in aposematic signalling and with which to direct future research. A greater appreciation of the extent of variation in aposematic species, and of the selective pressures and constraints which contribute to this once-paradoxical phenomenon, yields a new perspective for the field of aposematic signalling.

Recombination rate variation shapes barriers to introgression across butterfly genomes

Hybridisation and introgression can dramatically alter the relationships among groups of species, leading to phylogenetic discordance across the genome and between populations. Introgression can also erode species differences over time, but selection against introgression at certain loci acts to maintain postmating species barriers. Theory predicts that species barriers made up of many loci throughout the genome should lead to a broad correlation between introgression and recombination rate, which determines the extent to which selection on deleterious foreign alleles will affect neutral alleles at physically linked loci. Here, we describe the variation in genealogical relationships across the genome among three species of Heliconius butterflies: H. melpomene (mel), H. cydno (cyd), and H. timareta (tim), using whole genomes of 92 individuals, and ask whether this variation can be explained by heterogeneous barriers to introgression. We find that species relationships vary predictably at the chromosomal scale. By quantifying recombination rate and admixture proportions, we then show that rates of introgression are predicted by variation in recombination rate. This implies that species barriers are highly polygenic, with selection acting against introgressed alleles across most of the genome. In addition, long chromosomes, which have lower recombination rates, produce stronger barriers on average than short chromosomes. Finally, we find a consistent difference between two species pairs on either side of the Andes, which suggests differences in the architecture of the species barriers. Our findings illustrate how the combined effects of hybridisation, recombination, and natural selection, acting at multitudes of loci over long periods, can dramatically sculpt the phylogenetic relationships among species.

Patterns of Z chromosome divergence among Heliconius species highlight the importance of historical demography

Sex chromosomes are disproportionately involved in reproductive isolation and adaptation. In support of such a “large‐X” effect, genome scans between recently diverged populations and species pairs often identify distinct patterns of divergence on the sex chromosome compared to autosomes. When measures of divergence between populations are higher on the sex chromosome compared to autosomes, such patterns could be interpreted as evidence for faster divergence on the sex chromosome, that is “faster‐X”, barriers to gene flow on the sex chromosome. However, demographic changes can strongly skew divergence estimates and are not always taken into consideration. We used 224 whole‐genome sequences representing 36 populations from two Heliconius butterfly clades (H. erato and H. melpomene) to explore patterns of Z chromosome divergence. We show that increased divergence compared to equilibrium expectations can in many cases be explained by demographic change. Among Heliconius erato populations, for instance, population size increase in the ancestral population can explain increased absolute divergence measures on the Z chromosome compared to the autosomes, as a result of increased ancestral Z chromosome genetic diversity. Nonetheless, we do identify increased divergence on the Z chromosome relative to the autosomes in parapatric or sympatric species comparisons that imply postzygotic reproductive barriers. Using simulations, we show that this is consistent with reduced gene flow on the Z chromosome, perhaps due to greater accumulation of incompatibilities. Our work demonstrates the importance of taking demography into account to interpret patterns of divergence on the Z chromosome, but nonetheless provides evidence to support the Z chromosome as a strong barrier to gene flow in incipient Heliconius butterfly species.

What shapes the continuum of reproductive isolation? Lessons from Heliconius butterflies

The process by which species evolve can be illuminated by investigating barriers that limit gene flow between taxa. Recent radiations, such as Heliconius butterflies, offer the opportunity to compare isolation between pairs of taxa at different stages of ecological, geographical, and phylogenetic divergence. Here, we report a comparative analysis of existing and novel data in order to quantify the strength and direction of isolating barriers within a well-studied clade of Heliconius Our results highlight that increased divergence is associated with the accumulation of stronger and more numerous barriers to gene flow. Wing pattern is both under natural selection for Müllerian mimicry and involved in mate choice, and therefore underlies several isolating barriers. However, pairs which share a similar wing pattern also display strong reproductive isolation mediated by traits other than wing pattern. This suggests that, while wing pattern is a key factor for early stages of divergence, it may become facultative at later stages of divergence. Additional factors including habitat partitioning, hybrid sterility, and chemically mediated mate choice are associated with complete speciation. Therefore, although most previous work has emphasized the role of wing pattern, our comparative results highlight that speciation is a multi-dimensional process, whose completion is stabilized by many factors.

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.

The functional basis of wing patterning in Heliconius butterflies: the molecules behind mimicry

Wing-pattern mimicry in butterflies has provided an important example of adaptation since Charles Darwin and Alfred Russell Wallace proposed evolution by natural selection >150 years ago. The neotropical butterfly genus Heliconius played a central role in the development of mimicry theory and has since been studied extensively in the context of ecology and population biology, behavior, and mimicry genetics. Heliconius species are notable for their diverse color patterns, and previous crossing experiments revealed that much of this variation is controlled by a small number of large-effect, Mendelian switch loci. Recent comparative analyses have shown that the same switch loci control wing-pattern diversity throughout the genus, and a number of these have now been positionally cloned. Using a combination of comparative genetic mapping, association tests, and gene expression analyses, variation in red wing patterning throughout Heliconius has been traced back to the action of the transcription factor optix. Similarly, the signaling ligand WntA has been shown to control variation in melanin patterning across Heliconius and other butterflies. Our understanding of the molecular basis of Heliconius mimicry is now providing important insights into a variety of additional evolutionary phenomena, including the origin of supergenes, the interplay between constraint and evolvability, the genetic basis of convergence, the potential for introgression to facilitate adaptation, the mechanisms of hybrid speciation in animals, and the process of ecological speciation.

The diversification of Heliconius butterflies: what have we learned in 150 years?

Research into Heliconius butterflies has made a significant contribution to evolutionary biology. Here, we review our understanding of the diversification of these butterflies, covering recent advances and a vast foundation of earlier work. Whereas no single group of organisms can be sufficient for understanding life's diversity, after years of intensive study, research into Heliconius has addressed a wide variety of evolutionary questions. We first discuss evidence for widespread gene flow between Heliconius species and what this reveals about the nature of species. We then address the evolution and diversity of warning patterns, both as the target of selection and with respect to their underlying genetic basis. The identification of major genes involved in mimetic shifts, and homology at these loci between distantly related taxa, has revealed a surprising predictability in the genetic basis of evolution. In the final sections, we consider the evolution of warning patterns, and Heliconius diversity more generally, within a broader context of ecological and sexual selection. We consider how different traits and modes of selection can interact and influence the evolution of reproductive isolation.

Evolutionary consequences of climate-induced range shifts in insects

Range shifts can rapidly create new areas of geographic overlap between formerly allopatric taxa and evidence is accumulating that this can affect species persistence. We review the emerging literature on the short- and long-term consequences of these geographic range shifts. Specifically, we focus on the evolutionary consequences of novel species interactions in newly created sympatric areas by describing the potential (i) short-term processes acting on reproductive barriers between species and (ii) long-term consequences of range shifts on the stability of hybrid zones, introgression and ultimately speciation and extinction rates. Subsequently, we (iii) review the empirical literature on insects to evaluate which processes have been studied, and (iv) outline some areas that deserve increased attention in the future, namely the genomics of hybridisation and introgression, our ability to forecast range shifts and the impending threat from insect vectors and pests on biodiversity, human health and crop production. Our review shows that species interactions in de novo sympatric areas can be manifold, sometimes increasing and sometimes decreasing species diversity. A key issue that emerges is that climate-induced hybridisations in insects are much more widespread than anticipated and that rising temperatures and increased anthropogenic disturbances are accelerating the process of species mixing. The existing evidence only shows the tip of the iceberg and we are likely to see many more cases of species mixing following range shifts in the near future.

An introgressed wing pattern acts as a mating cue

Heliconius butterflies provide good examples of both homoploid hybrid speciation and ecological speciation. In particular, examples of adaptive introgression have been detected among the subspecies of Heliconius timareta, which acquired red color pattern elements from H. melpomene. We tested whether the introgression of red wing pattern elements into H. timareta florencia might also be associated with incipient reproductive isolation (RI) from its close relative, H. timareta subsp. nov., found in the eastern Andes. No choice experiments show a 50% reduction in mating between females of H. t. subsp. nov. and males of H .t. florencia, but not in the reciprocal direction. In choice experiments using wing models, males of H. timareta subsp. nov. approach and court red phenotypes less than their own, whereas males of H. t. florencia prefer models with a red phenotype. Intrinsic postzygotic isolation was not detected in crosses between these H. timareta races. These results suggest that a color pattern trait gained by introgression is triggering RI between H. timareta subsp. nov. and H. t. florencia.

Introgression of wing pattern alleles and speciation via homoploid hybridization in Heliconius butterflies: a review of evidence from the genome

The diverse Müllerian mimetic wing patterns of neotropical Heliconius (Nymphalidae) have been proposed to be not only aposematic signals to potential predators, but also intra- and interspecific recognition signals that allow the butterflies to maintain their specific identities, and which perhaps drive the process of speciation, as well. Adaptive features under differential selection that also serve as cues for assortative mating have been referred to as 'magic traits', which can drive ecological speciation. Such traits are expected to exhibit allelic differentiation between closely related species with ongoing gene flow, whereas unlinked neutral traits are expected to be homogenized to a greater degree by introgression. However, recent evidence suggests that interspecific hybridization among Heliconius butterflies may have resulted in adaptive introgression of these very same traits across species boundaries, and in the evolution of new species by homoploid hybrid speciation. The theory and data supporting various aspects of the apparent paradox of 'magic trait' introgression are reviewed, with emphasis on population genomic comparisons of Heliconius melpomene and its close relatives.

Hybrid speciation in Heliconius butterflies? A review and critique of the evidence

The evidence supporting the recent hypothesis of a homoploid hybrid origin for the butterfly species Heliconius heurippa is evaluated. Data from selective breeding experiments, mate-choice studies, and a wide variety of DNA markers are reviewed, and an alternative hypothesis for the origin of the species and its close relatives is proposed. A scenario of occasional red wing-pattern mutations in peripheral populations of Heliconius cydno with subsequent adaptive convergence towards sympatric mimicry rings involving H. melpomene and H. erato is offered as an alternative to the HHS hypothesis. Recent twists of this tale are addressed in a postscript.

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.

Genetic evidence for hybrid trait speciation in heliconius butterflies

Homoploid hybrid speciation is the formation of a new hybrid species without change in chromosome number. So far, there has been a lack of direct molecular evidence for hybridization generating novel traits directly involved in animal speciation. Heliconius butterflies exhibit bright aposematic color patterns that also act as cues in assortative mating. Heliconius heurippa has been proposed as a hybrid species, and its color pattern can be recreated by introgression of the H. m. melpomene red band into the genetic background of the yellow banded H. cydno cordula. This hybrid color pattern is also involved in mate choice and leads to reproductive isolation between H. heurippa and its close relatives. Here, we provide molecular evidence for adaptive introgression by sequencing genes across the Heliconius red band locus and comparing them to unlinked wing patterning genes in H. melpomene, H. cydno, and H. heurippa. 670 SNPs distributed among 29 unlinked coding genes (25,847bp) showed H. heurippa was related to H. c. cordula or the three species were intermixed. In contrast, among 344 SNPs distributed among 13 genes in the red band region (18,629bp), most showed H. heurippa related with H. c. cordula, but a block of around 6,5kb located in the 3' of a putative kinesin gene grouped H. heurippa with H. m. melpomene, supporting the hybrid introgression hypothesis. Genealogical reconstruction showed that this introgression occurred after divergence of the parental species, perhaps around 0.43Mya. Expression of the kinesin gene is spatially restricted to the distal region of the forewing, suggesting a mechanism for pattern regulation. This gene therefore constitutes the first molecular evidence for adaptive introgression during hybrid speciation and is the first clear candidate for a Heliconius wing patterning locus.

Review. Hybrid trait speciation and Heliconius butterflies

Homoploid hybrid speciation (HHS) is the establishment of a novel species through introgressive hybridization without a change in chromosome number. We discuss different routes by which this might occur and propose a novel term, 'hybrid trait speciation', which combines the idea that hybridization can generate adaptive novelty with the 'magic trait' model of ecological speciation. Heliconius butterflies contain many putative examples of hybrid colour patterns, but only recently has the HHS hypothesis been tested explicitly in this group. Molecular data has shown evidence for gene flow between many distinct species. Furthermore, the colour pattern of Heliconius heurippa can be recreated in laboratory crosses between Heliconius melpomene and Heliconius cydno and, crucially, plays a role in assortative mating between the three species. Nonetheless, although the genome of H. heurippa shows evidence for hybridization, it is not a mosaic of the two parental species. Instead, ongoing hybridization has likely blurred any signal of the original speciation event. We argue that where hybridization leads to novel adaptive traits that also cause reproductive isolation, it is likely to trigger speciation.

Gene flow and the genealogical history of Heliconius heurippa

Background

The neotropical butterfly Heliconius heurippa has a hybrid colour pattern, which also contributes to reproductive isolation, making it a likely example of hybrid speciation. Here we used phylogenetic and coalescent-based analyses of multilocus sequence data to investigate the origin of H. heurippa.

Results

We sequenced a mitochondrial region (CoI and CoII), a sex-linked locus (Tpi) and two autosomal loci (w and sd) from H. heurippa and the putative parental species, H. cydno and H. melpomene. These were analysed in combination with data from two previously sequenced autosomal loci, Dll and Inv. H. heurippa was monophyletic at mtDNA and Tpi, but showed a shared distribution of alleles derived from both parental lineages at all four autosomal loci. Estimates of genetic differentiation showed that H. heurippa is closer to H. cydno at mtDNA and three autosomal loci, intermediate at Tpi, and closer to H. melpomene at Dll. Using coalescent simulations with the Isolation-Migration model (IM), we attempted to establish the incidence of gene flow in the origin of H. heurippa. This analysis suggested that ongoing introgression is frequent between all three species and variable in extent between loci.

Conclusion

Introgression, which is a necessary precursor of hybrid speciation, seems to have also blurred the coalescent history of these species. The origin of Heliconius heurippa may have been restricted to introgression of few colour pattern genes from H. melpomene into the H. cydno genome, with little evidence of genomic mosaicism.

Do pollen feeding, pupal-mating and larval gregariousness have a single origin in Heliconius butterflies? Inferences from multilocus DNA sequence data

Phylogenetic information is useful in understanding the evolutionary history of adaptive traits. Here, we present a well-resolved phylogenetic hypothesis for Heliconius butterflies and related genera. We use this tree to investigate the evolution of three traits, pollen feeding, pupal-mating behaviour and larval gregariousness. Phylogenetic relationships among 60 Heliconiina species (86% of the subtribe) were inferred from partial DNA sequences of the mitochondrial genes cytochrome oxidase I, cytochrome oxidase II and 16S rRNA, and fragments of the nuclear genes elongation factor-1α, apterous, decapentaplegic and wingless (3834 bp in total). The results corroborate previous hypotheses based on sequence data in showing that Heliconius is paraphyletic, with Laparus doris and Neruda falling within the genus, demonstrating a single origin for pollen feeding but with a loss of the trait in Neruda. However, different genes are not congruent in their placement of Neruda; therefore, monophyly of the pollen feeding species cannot be ruled out. There is also a highly supported monophyletic ‘pupal-mating clade’ suggesting that pupal mating behaviour evolved only once in the Heliconiina. Additionally, we observed at least three independent origins for larval gregariousness from a solitary ancestor, showing that gregarious larval behaviour arose after warning coloration.

Genetic differentiation and natural hybridization between the Sardinian endemic Maniola nurag and the European Maniola jurtina

The Mediterranean island of Sardinia is known for its multitude of unique genetic lineages. We view one of them in a larger phylogeographic context. The endemic Sardinian Meadow Brown butterfly, Maniola nurag, is restricted to the mountainous areas of the island, whereas its widespread close relative, Maniola jurtina, also occurs on the coast. At intermediate altitudes the species' distributions overlap. There, a number of individuals exhibit phenotypic characteristics intermediate between the two species. We examined patterns of intra- and interpopulation variation in 10 M. nurag populations from Sardinia and 16 M. jurtina populations from Sardinia and continental Europe, as well as 17 intermediate individuals, sampled in 1999-2002, by means of allozyme markers, combining it with a morphometric analysis based on 18 wing-characters of 52 males. At the 15 loci studied (aldolase, aat-1, aat-2, g6pdh, gpd, idh-1, idh-2, mdh-1, mdh-2, mpi, me, leu-ala, pgi, pgm, and 6pgdh), 76 different alleles were detected, 63 of which were shared by M. nurag and M. jurtina. None of the loci was found to be alternatively fixed between the two species. In that respect, this study testifies to the difficulties that may arise when trying to identify hybrids from genotypic data. Levels of genetic variation in island populations (M. jurtina: H(O) = 0.137-0.189; M. nurag: H(O) = 0.141-0.270) were comparable to those of mainland M. jurtina (H(O) = 0.141-0.236). A Bayesian admixture analysis supported the hypothesis of mixed (hybrid) ancestry of individuals occurring at intermediate altitudes. Similarly, neighbour-joining and unweighted pair-group method with arithmetic averaging (UPGMA) analyses, as well as morphometrics hinted at the existence of a Maniola-hybrid zone in Sardinia at intermediate altitudes. We discuss the results in the light of the phylogeography of other Sardinian taxa with the aim to reach a general understanding of the biogeographic history of this island's endemic species.

No genomic mosaicism in a putative hybrid butterfly species

Recent descriptions of hybrid animal species have spurred interest in this phenomenon, but little genomic data exist to support it. Here, we use frequency variation for 657 amplified fragment length polymorphism (AFLP) markers and DNA sequence variation from 16 genes to determine whether the genome of Heliconius pachinus, a suspected hybrid butterfly species, is a mixture of the putative parental species, Heliconius cydno and Heliconius melpomene. Despite substantial shared genetic variation among all three species, we show that the genome of H. pachinus is not a mosaic; both AFLP and DNA sequence data overwhelmingly associate H. pachinus with just one of the potential parents, H. cydno. This pattern also applies to the gene wingless, which is tightly linked to the locus that determines forewing colour-one specific H. pachinus trait that has been hypothesized to have originated from H. melpomene. As a whole, the data support a traditional, bifurcating model of speciation in which H. pachinus split from a common ancestor with H. cydno without a genetic contribution from H. melpomene. However, comparison of our data to DNA sequence data for another putative hybrid Heliconius species, Heliconius heurippa, suggests that the H. heurippa genome may be a mosaic.

Asymmetric postmating isolation: Darwin's corollary to Haldane's rule

Asymmetric postmating isolation, where reciprocal interspecific crosses produce different levels of fertilization success or hybrid sterility/inviability, is very common. Darwin emphasized its pervasiveness in plants, but it occurs in all taxa assayed. This asymmetry often results from Dobzhansky-Muller incompatibilities (DMIs) involving uniparentally inherited genetic factors (e.g., gametophyte-sporophyte interactions in plants or cytoplasmic-nuclear interactions). Typically, unidirectional (U) DMIs act simultaneously with bidirectional (B) DMIs between autosomal loci that affect reciprocal crosses equally. We model both classes of two-locus DMIs to make quantitative and qualitative predictions concerning patterns of isolation asymmetry in parental species crosses and in the hybrid F(1) generation. First, we find conditions that produce expected differences. Second, we present a stochastic analysis of DMI accumulation to predict probable levels of asymmetry as divergence time increases. We find that systematic interspecific differences in relative rates of evolution for autosomal vs. nonautosomal loci can lead to different expected F(1) fitnesses from reciprocal crosses, but asymmetries are more simply explained by stochastic differences in the accumulation of U DMIs. The magnitude of asymmetry depends primarily on the cumulative effects of U vs. B DMIs (which depend on heterozygous effects of DMIs), the average number of DMIs required to produce complete reproductive isolation (more asymmetry occurs when fewer DMIs are required), and the shape of the function describing how fitness declines as DMIs accumulate. Comparing our predictions to data from diverse taxa indicates that unidirectional DMIs, specifically involving sex chromosomes, cytoplasmic elements, and maternal effects, are likely to play an important role in postmating isolation.

Reconstructing the history of selection during homoploid hybrid speciation

This study aims to identify selection pressures during the historical process of homoploid hybrid speciation in three Helianthus (sunflower) hybrid species. If selection against intrinsic genetic incompatibilities (fertility selection) or for important morphological/ecological traits (phenotypic selection) were important in hybrid speciation, we would expect this selection to have influenced the parentage of molecular markers or chromosomal segments in the hybrid species' genomes. To infer past selection, we compared the parentage of molecular markers in high-density maps of the three hybrid species with predicted marker parentage from an analysis of fertility selection in artificial hybrids and from the directions of quantitative trait loci effects with respect to the phenotypes of the hybrid species. Multiple logistic regression models were consistent with both fertility and phenotypic selection in all three species. To further investigate traits under selection, we used a permutation test to determine whether marker parentage predicted from groups of functionally related traits differed from neutral expectations. Our results suggest that trait groups associated with ecological divergence were under selection during hybrid speciation. This study presents a new method to test for selection and supports earlier claims that fertility selection and phenotypic selection on ecologically relevant traits have operated simultaneously during sunflower hybrid speciation.

Natural hybridization in heliconiine butterflies: the species boundary as a continuum

BACKGROUND

To understand speciation and the maintenance of taxa as separate entities, we need information about natural hybridization and gene flow among species.

RESULTS

Interspecific hybrids occur regularly in Heliconius and Eueides (Lepidoptera: Nymphalidae) in the wild: 26-29% of the species of Heliconiina are involved, depending on species concept employed. Hybridization is, however, rare on a per-individual basis. For one well-studied case of species hybridizing in parapatric contact (Heliconius erato and H. himera), phenotypically detectable hybrids form around 10% of the population, but for species in sympatry hybrids usually form less than 0.05% of individuals. There is a roughly exponential decline with genetic distance in the numbers of natural hybrids in collections, both between and within species, suggesting a simple "exponential failure law" of compatibility as found in some prokaryotes.

CONCLUSION

Hybridization between species of Heliconius appears to be a natural phenomenon; there is no evidence that it has been enhanced by recent human habitat disturbance. In some well-studied cases, backcrossing occurs in the field and fertile backcrosses have been verified in insectaries, which indicates that introgression is likely, and recent molecular work shows that alleles at some but not all loci are exchanged between pairs of sympatric, hybridizing species. Molecular clock dating suggests that gene exchange may continue for more than 3 million years after speciation. In addition, one species, H. heurippa, appears to have formed as a result of hybrid speciation. Introgression may often contribute to adaptive evolution as well as sometimes to speciation itself, via hybrid speciation. Geographic races and species that coexist in sympatry therefore form part of a continuum in terms of hybridization rates or probability of gene flow. This finding concurs with the view that processes leading to speciation are continuous, rather than sudden, and that they are the same as those operating within species, rather than requiring special punctuated effects or complete allopatry. Although not qualitatively distinct from geographic races, nor "real" in terms of phylogenetic species concepts or the biological species concept, hybridizing species of Heliconius are stably distinct in sympatry, and remain useful groups for predicting morphological, ecological, behavioural and genetic characteristics.

Speciation by hybridization in Heliconius butterflies

Speciation is generally regarded to result from the splitting of a single lineage. An alternative is hybrid speciation, considered to be extremely rare, in which two distinct lineages contribute genes to a daughter species. Here we show that a hybrid trait in an animal species can directly cause reproductive isolation. The butterfly species Heliconius heurippa is known to have an intermediate morphology and a hybrid genome, and we have recreated its intermediate wing colour and pattern through laboratory crosses between H. melpomene, H. cydno and their F1 hybrids. We then used mate preference experiments to show that the phenotype of H. heurippa reproductively isolates it from both parental species. There is strong assortative mating between all three species, and in H. heurippa the wing pattern and colour elements derived from H. melpomene and H. cydno are both critical for mate recognition by males.