Phylogeography of Heliconius cydno and its closest relatives: disentangling their origin and diversification

The origins of the extraordinary diversity within the Neotropics have long fascinated biologists and naturalists. Yet, the underlying factors that have given rise to this diversity remain controversial. To test the relative importance of Quaternary climatic change and Neogene tectonic and paleogeographic reorganizations in the generation of biodiversity, we examine intraspecific variation across the Heliconius cydno radiation and compare this variation to that within the closely related Heliconius melpomene and Heliconius timareta radiations. Our data, which consist of both mtDNA and genome‐scan data from nearly 2250 amplified fragment length polymorphism (AFLP) loci, reveal a complex history of differentiation and admixture at different geographic scales. Both mtDNA and AFLP phylogenies suggest that H. timareta and H. cydno are probably geographic extremes of the same radiation that probably diverged from H. melpomene prior to the Pliocene–Pleistocene boundary, consistent with hypotheses of diversification that rely on geological events in the Pliocene. The mtDNA suggests that this radiation originated in Central America or the northwestern region of South America, with a subsequent colonization of the eastern and western slopes of the Andes. Our genome‐scan data indicate significant admixture among sympatric H. cydno/H. timareta and H. melpomene populations across the extensive geographic ranges of the two radiations. Within H. cydno, both mtDNA and AFLP data indicate significant population structure at local scales, with strong genetic differences even among adjacent H. cydno colour pattern races. These genetic patterns highlight the importance of past geoclimatic events, intraspecific gene flow, and local population differentiation in the origin and establishment of new adaptive forms.     

A Narrow Heliconius cydno (Nymphalidae; Heliconiini) Hybrid Zone With Differences in Morph Sex Ratios1

In the western cordillera of the Cauca valley, Colombia, there is a narrow hybrid zone between two races of Heliconius cydno, one of which is polymorphic. Three large broods show that most of the phenotypic variation observed can be explained by four loci of major effect, named Sb, Yb, L, and K. Two further loci, G and Wo, were identified that control more minor phenotypic variation. Sb, Yb, and Wo are linked and the latter is differentially expressed between the sexes. The transition between H. c. cydnides in the north and H. c. weymeri to the south occurs over approximately 15 km. Collections from a single site near the center of the hybrid zone show that gene frequencies have been stable over the 10 yr from 1991 to 2001 and that color‐pattern genes are in Hardy–Weinberg equilibrium, with little evidence for linkage disequilibrium. This suggests that mating is random between color‐pattern forms. Nonetheless, there is evidence for variation in the sex ratio, with parental phenotypes and the most melanic hybrid phenotypes showing a fairly even sex ratio compared to a strong male bias observed in the remaining hybrids. We hypothesize that this might be explained by differential selection between sexes and phenotypes, phenomena that could have important implications for hybrid zone analysis and the evolution of Müllerian mimicry.     

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.     

Ecological and genetic associations across a Heliconius hybrid zone

Differences in habitat use can bridge early and late stages of speciation by initiating assortative mating. Heliconius colour pattern races might select habitats over which each pattern confers a relative fitness advantage because signal efficacy of wing patterns can vary by environment. Thus habitat preferences could serve to promote the evolution of mimetic colour patterns for mate choice. Here I compare colour pattern genotype and phenotype frequencies to environmental variation across the H. erato hydara x H. erato erato hybrid zone in French Guiana to determine whether races exhibit habitat preferences. I found that genotype and phenotype frequencies correspond to differences in land cover moreso than to other environmental factors. Temporal shifts in colour pattern genotypes, phenotypes and land cover also were associated at individual sample sites, which further suggests that H. erato races differ in habitat use and that habitat preferences may promote speciation among Heliconius butterflies.     

Multiple sources of reproductive isolation in a bimodal butterfly hybrid zone

An important evolutionary question concerns whether one or many barriers are involved in the early stages of speciation. We examine pre‐ and post‐zygotic reproductive barriers between two species of butterflies (Heliconius erato chestertonii and H. e. venus) separated by a bimodal hybrid zone in the Cauca Valley, Colombia. We show that there is both strong pre‐ and post‐mating reproductive isolation, together leading to a 98% reduction in gene flow between the species. Pre‐mating isolation plays a primary role, contributing strongly to this isolation (87%), similar to previous examples in Heliconius. Post‐mating isolation was also strong, with absence of Haldane’s rule, but an asymmetric reduction in fertility (< 11%) in inter‐specific crosses depending on maternal genotype. In summary, this is one of the first examples of post‐zygotic reproductive isolation playing a significant role in early stages of parapatric speciation in Heliconius and demonstrates the importance of multiple barriers to gene flow in the speciation process.     

Hybrid zones and the speciation continuum in Heliconius butterflies

Tropical butterflies in the genus Heliconius have long been models in the study of the stages of speciation. Heliconius are unpalatable to predators, and many species are notable for multiple geographic populations with striking warning colour pattern differences associated with Müllerian mimicry. A speciation continuum is evident in Heliconius hybrid zones. Examples range from hybrid zones across which (a) there is little genetic differentiation other than at mimicry loci, but where hybrids are common, (b) to ‘bimodal‘ hybrid zones with strong genetic divergence and few hybrids, (c) through to ‘good’ sympatric species, with hybridization extremely rare or absent. Now, in this issue of Molecular Ecology, Arias et al. ( 2012 ) have found an intermediate case in Colombian Heliconius cydno showing evidence for assortative mating and molecular differences, but where hybrids are abundant.     

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.     

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

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

Polyphyly and gene flow between non-sibling Heliconius species

Background  

The view that gene flow between related animal species is rare and evolutionarily unimportant largely antedates sensitive molecular techniques. Here we use DNA sequencing to investigate a pair of morphologically and ecologically divergent, non-sibling butterfly species, Heliconius cydno and H. melpomene (Lepidoptera: Nymphalidae), whose distributions overlap in Central and Northwestern South America.     

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

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

Gene flow and the genealogical history of <Emphasis Type="Italic">Heliconius heurippa</Emphasis>

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.     

Cryptic genetic and wing pattern diversity in a mimetic Heliconius butterfly

Despite rampant colour pattern diversity in South America, Heliconius erato exhibits a ‘postman’ wing pattern throughout most of Central America. We examined genetic variation across the range of H. erato, including dense sampling in Central America, and discovered a deep genetic break, centred on the mountain range that runs through Costa Rica. This break is characterized by a novel mitochondrial lineage, which is nearly fixed in northern Central America, that branches basal to all previously described mitochondrial diversity in the species. Strong genetic differentiation also appears in Z‐linked and autosomal markers, and it is further associated with a distinct, but subtle, shift in wing pattern phenotype. Comparison of clines in wing phenotype, mtDNA and nuclear markers indicate they are all centred on the mountains dividing Costa Rica, but that cline width differs among data sets. Phylogeographical analyses, accounting for this new diversity, rewrite our understanding of mimicry evolution in this system. For instance, these results suggest that H. erato originated west of the Andes, perhaps in Central America, and as many as 1 million years before its co‐mimic, H. melpomene. Overall our data indicate that neutral genetic markers and colour pattern loci are congruent and converge on the same hypothesis—H. erato originated in northwest South America or Central America with a ‘postman’ phenotype and then radiated into the wealth of colour patterns present today.     

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.     

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.     

Estimating the relative fitness of local adaptive peaks: the aerodynamic costs of flight in mimetic passion-vine butterflies Heliconius

In a related paper, we demonstrated that mimetic Heliconius butterflies have converged in wing-beat frequency and degree of asymmetry in the wing motion, whereas sister species are dissimilar in these same traits. Warning signals of sympatric, distasteful species converge in evolutionary models in order to educate their predators more efficiently that the signal is associated with unprofitable prey. Barring other constraints, the behaviours of the different co-mimetic pairs should ultimately converge on that behaviour which minimizes the energetic cost of flight. We estimated the energetic cost of each mimic''s flight behaviour in order to predict the difference in height of each fitness peak and the direction of convergent selection qualitatively. Following adjustments for body mass, mimetic Heliconius melpomene and Heliconius erato required more aerodynamic power than Heliconius cydno and Heliconius sapho. This difference was attributed to the slower flight speeds and higher wing-beat frequencies of H. melpomene and H. erato. Consequently, H. melpomene and H. erato expended more energy per unit distance per unit body mass than H. cydno and H. sapho. However, differences in body mass may equalize energy budgets and stabilize the sympatric coexistence of the two pairs of co-mimics.     

Comparative population genetics of a mimicry locus among hybridizing Heliconius butterfly species

The comimetic Heliconius butterfly species pair, H. erato and H. melpomene, appear to use a conserved Mendelian switch locus to generate their matching red wing patterns. Here we investigate whether H. cydno and H. pachinus, species closely related to H. melpomene, use this same switch locus to generate their highly divergent red and brown color pattern elements. Using an F2 intercross between H. cydno and H. pachinus, we first map the genomic positions of two novel red/brown wing pattern elements; the G locus, which controls the presence of red vs brown at the base of the ventral wings, and the Br locus, which controls the presence vs absence of a brown oval pattern on the ventral hind wing. The results reveal that the G locus is tightly linked to markers in the genomic interval that controls red wing pattern elements of H. erato and H. melpomene. Br is on the same linkage group but approximately 26 cM away. Next, we analyze fine-scale patterns of genetic differentiation and linkage disequilibrium throughout the G locus candidate interval in H. cydno, H. pachinus and H. melpomene, and find evidence for elevated differentiation between H. cydno and H. pachinus, but no localized signature of association. Overall, these results indicate that the G locus maps to the same interval as the locus controlling red patterning in H. melpomene and H. erato. This, in turn, suggests that the genes controlling red pattern elements may be homologous across Heliconius, supporting the hypothesis that Heliconius butterflies use a limited suite of conserved genetic switch loci to generate both convergent and divergent wing patterns.     

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

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

Evidence for concerted movement of nuclear and mitochondrial clines in a lizard hybrid zone

Moving hybrid zones provide compelling examples of evolution in action, yet long‐term studies that test the assumptions of hybrid zone stability are rare. Using replicated transect samples collected over a 10‐year interval from 2002 to 2012, we find evidence for concerted movement of genetic clines in a plateau fence lizard hybrid zone (Sceloporus tristichus) in Arizona. Cline‐fitting analyses of SNP and mtDNA data both provide evidence that the hybrid zone shifted northward by approximately 2 km during the 10‐year interval. For each sampling period, the mtDNA cline centre is displaced from the SNP cline centre and maintaining an introgression distance of approximately 3 km. The northward expansion of juniper trees into the Little Colorado River Basin in the early 1900s provides a plausible mechanism for hybrid zone formation and movement, and a broadscale quantification of recent land cover change provides support for increased woody species encroachment at the southern end of the hybrid zone. However, population processes can also contribute to hybrid zone movement, and the current stability of the ecotone habitats in the centre of the hybrid zone suggests that movement could decelerate in the future.