Chemical signals act as the main reproductive barrier between sister and mimetic Heliconius butterflies.

Colour pattern is the main trait that drives mate recognition between Heliconius species that are phylogenetically close. However, when this cue is compromised such as in cases of mimetic, sympatric and closely related species, alternative mating signals must evolve to ensure reproductive isolation and species integrity. The closely related species Heliconius melpomene malleti and H. timareta florencia occur in the same geographical region, and despite being co-mimics, they display strong reproductive isolation. In order to test which cues differ between species, and potentially contribute to reproductive isolation, we quantified differences in the wing phenotype and the male chemical profile. As expected, the wing colour pattern was indistinguishable between the two species, while the chemical profile of the androconial and genital males' extracts showed marked differences. We then conducted behavioural experiments to study the importance of these signals in mate recognition by females. In agreement with our previous results, we found that chemical blends and not wing colour pattern drive the preference of females for conspecific males. Also, experiments with hybrid males and females suggested an important genetic component for both chemical production and preference. Altogether, these results suggest that chemicals are the major reproductive barrier opposing gene flow between these two sister and co-mimic 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.

The comparative landscape of duplications in Heliconius melpomene and Heliconius cydno.

Gene duplications can facilitate adaptation and may lead to interpopulation divergence, causing reproductive isolation. We used whole-genome resequencing data from 34 butterflies to detect duplications in two Heliconius species, Heliconius cydno and Heliconius melpomene. Taking advantage of three distinctive signals of duplication in short-read sequencing data, we identified 744 duplicated loci in H. cydno and H. melpomene and evaluated the accuracy of our approach using single-molecule sequencing. We have found that duplications overlap genes significantly less than expected at random in H. melpomene, consistent with the action of background selection against duplicates in functional regions of the genome. Duplicate loci that are highly differentiated between H. melpomene and H. cydno map to four different chromosomes. Four duplications were identified with a strong signal of divergent selection, including an odorant binding protein and another in close proximity with a known wing colour pattern locus that differs between the two species.

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.

Supergenes and their role in evolution.

Adaptation is commonly a multidimensional problem, with changes in multiple traits required to match a complex environment. This is epitomized by balanced polymorphisms in which multiple phenotypes co-exist and are maintained in a population by a balance of selective forces. Consideration of such polymorphisms led to the concept of the supergene, where alternative phenotypes in a balanced polymorphism segregate as if controlled by a single genetic locus, resulting from tight genetic linkage between multiple functional loci. Recently, the molecular basis for several supergenes has been resolved. Thus, major chromosomal inversions have been shown to be associated with polymorphisms in butterflies, ants and birds, offering a mechanism for localised reduction in recombination. In several examples of plant self-incompatibility, the functional role of multiple elements within the supergene architecture has been demonstrated, conclusively showing that balanced polymorphism can be maintained at multiple coadapted and tightly linked elements. Despite recent criticism, we argue that the supergene concept remains relevant and is more testable than ever with modern molecular methods.

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.

Hybridization and speciation.

Hybridization has many and varied impacts on the process of speciation. Hybridization may slow or reverse differentiation by allowing gene flow and recombination. It may accelerate speciation via adaptive introgression or cause near-instantaneous speciation by allopolyploidization. It may have multiple effects at different stages and in different spatial contexts within a single speciation event. We offer a perspective on the context and evolutionary significance of hybridization during speciation, highlighting issues of current interest and debate. In secondary contact zones, it is uncertain if barriers to gene flow will be strengthened or broken down due to recombination and gene flow. Theory and empirical evidence suggest the latter is more likely, except within and around strongly selected genomic regions. Hybridization may contribute to speciation through the formation of new hybrid taxa, whereas introgression of a few loci may promote adaptive divergence and so facilitate speciation. Gene regulatory networks, epigenetic effects and the evolution of selfish genetic material in the genome suggest that the Dobzhansky-Muller model of hybrid incompatibilities requires a broader interpretation. Finally, although the incidence of reinforcement remains uncertain, this and other interactions in areas of sympatry may have knock-on effects on speciation both within and outside regions of hybridization.

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.

Variable extent of sex-biased dispersal in a strongly polygynous mammal.

For mammals with a polygynous mating system, dispersal is expected to be male-biased. However, with the increase in empirical studies, discrepancies are arising between the expected and observed direction/extent of the bias in dispersal. In this study, we assessed sex-biased dispersal in red deer (Cervus elaphus) on 13 estates from the Scottish Highlands. A total of 568 adult individuals were genotyped at 21 microsatellite markers and sequenced for 821 bp of the mitochondrial control region. Estimates of population structure with mitochondrial sequences were eight times larger than that obtained with microsatellite data (F(st'-mtDNA) = 0.831; F(st'-micros) = 0.096) indicating overall male-biased dispersal in the study area. Comparisons of microsatellite data between the sexes indicated a predominance of male-biased dispersal in the study area but values of F(ST) and relatedness were only slighter larger for females. Individual-based spatial autocorrelation analysis generated a similar pattern of relatedness across geographical distances for both sexes, with differences only significant at two distance intervals (25-30 and 70-112 km). Patterns of relatedness differed between estates, male biased-dispersal was detected in eight estates but no sex-biased dispersal was found in the remaining five. Neither population density nor landscape cover was found to be associated with the patterns of relatedness found across the estates. Differences in management strategies that could influence age structure, sex ratio and dispersal behaviour are proposed as potential factors influencing the relatedness patterns observed. This study provides new insights on dispersal of a strongly polygynous mammal at geographical scales relevant for management and conservation.

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.

Out of the Andes: patterns of diversification in clearwing butterflies.

Global biodiversity peaks in the tropical forests of the Andes, a striking geological feature that has likely been instrumental in generating biodiversity by providing opportunities for both vicariant and ecological speciation. However, the role of these mountains in the diversification of insects, which dominate biodiversity, has been poorly explored using phylogenetic methods. Here we study the role of the Andes in the evolution of a diverse Neotropical insect group, the clearwing butterflies. We used dated species-level phylogenies to investigate the time course of speciation and to infer ancestral elevation ranges for two diverse genera. We show that both genera likely originated at middle elevations in the Andes in the Middle Miocene, contrasting with most published results in vertebrates that point to a lowland origin. Although we detected a signature of vicariance caused by the uplift of the Andes at the Miocene-Pliocene boundary, most sister species were parapatric without any obvious vicariant barrier. Combined with an overall decelerating speciation rate, these results suggest an important role for ecological speciation and adaptive radiation, rather than simple vicariance.

Genetic diversity and population structure of Scottish Highland red deer (Cervus elaphus) populations: a mitochondrial survey.

The largest population of red deer (Cervus elaphus) in Europe is found in Scotland. However, human impacts through hunting and introduction of foreign deer stock have disturbed the population's genetics to an unknown extent. In this study, we analysed mitochondrial control region sequences of 625 individuals to assess signatures of human and natural historical influence on the genetic diversity and population structure of red deer in the Scottish Highlands. Genetic diversity was high with 74 haplotypes found in our study area (115 x 87 km). Phylogenetic analyses revealed that none of the individuals had introgressed mtDNA from foreign species or subspecies of deer and only suggested a very few localized red deer translocations among British localities. A haplotype network and population analyses indicated significant genetic structure (Phi(ST)=0.3452, F(ST)=0.2478), largely concordant with the geographical location of the populations. Mismatch distribution analysis and neutrality tests indicated a significant population expansion for one of the main haplogroups found in the study area, approximately dated c. 8200 or 16 400 years ago when applying a fast or slow mutation rate, respectively. Contrary to general belief, our results strongly suggest that native Scottish red deer mtDNA haplotypes have persisted in the Scottish Highlands and that the population retains a largely natural haplotype diversity and structure in our study area.

Butterfly speciation and the distribution of gene effect sizes fixed during adaptation.

Mimicry has had a significant historical influence as a tractable system for studying adaptation and is known to play a role in speciation. Here, we discuss recent theoretical treatment of adaptive walks to local adaptive peaks and contrast this with the adaptive landscape of mimicry. Evolution of novel Müllerian mimicry patterns almost certainly involves substitution of a major mutation to provide an initial similarity to the model, such that major gene effects are expected to an even greater degree than for other adaptive traits. The likelihood of large adaptive peak shifts in mimicry evolution may therefore promote speciation. In addition, mimicry adaptive peaks are determined by the local abundance of particular patterns and may be more fluid than the case for other traits. It will therefore be of considerable interest to test empirically the distribution of effect sizes fixed during mimicry evolution. Here, we show the feasibility of this by presenting a preliminary quantitative trait locus (QTL) analysis of Heliconius colour patterns. This shows that a number of modifier loci of different effect sizes influence forewing band morphology. We also show multiple pleiotropic effects of major Heliconius patterning loci and discuss the likelihood of multiple substitutions at the same loci in pattern evolution, which would inflate the importance of major loci in QTL analysis of the gene effect sizes. Analyses such as these have the potential to uncover the genetic architecture of both within and between species adaptive differences.

Interspecific sexual attraction because of convergence in warning colouration: is there a conflict between natural and sexual selection in mimetic species?

When species converge in their colour patterns because of mimicry, and those patterns are also used in mate recognition, there is a probability of conflicting selection pressures. Closely related species that mimic one another are particularly likely to face such confusion because of similarities in their courtship behaviour and ecology. We conducted experiments in greenhouse conditions to study interspecific attraction between two mimetic butterfly species, Heliconius erato and Heliconius melpomene. Both species spent considerable time approaching and courting females of the co-mimic species. Experiments using wing models demonstrated the importance of colour pattern in this interspecific attraction. Although males of H. melpomene were attracted to their co-mimics as much as to their own females, H. erato males were more efficient at distinguishing conspecifics, possibly using wing odours. Although preliminary, these results suggest that the use of additional cues may have evolved in H. erato to reduce the cost of convergence in visual signals with H. melpomene. Overall, our results showed that there might be a cost of mimetic convergence because of a reduction in the efficiency of species recognition. Such cost may contribute to explain the apparently stable diversity in Müllerian mimetic patterns in many tropical butterfly assemblages.

Landscape features affect gene flow of Scottish Highland red deer (Cervus elaphus).

Landscape features have been shown to strongly influence dispersal and, consequently, the genetic population structure of organisms. Studies quantifying the effect of landscape features on gene flow of large mammals with high dispersal capabilities are rare and have mainly been focused at large geographical scales. In this study, we assessed the influence of several natural and human-made landscape features on red deer gene flow in the Scottish Highlands by analysing 695 individuals for 21 microsatellite markers. Despite the relatively small scale of the study area (115 x 87 km), significant population structure was found using F-statistics (F(ST) = 0.019) and the program structure, with major differentiation found between populations sampled on either side of the main geographical barrier (the Great Glen). To assess the effect of landscape features on red deer population structure, the ArcMap GIS was used to create cost-distance matrices for moving between populations, using a range of cost values for each of the landscape features under consideration. Landscape features were shown to significantly affect red deer gene flow as they explained a greater proportion of the genetic variation than the geographical distance between populations. Sea lochs were found to be the most important red deer gene flow barriers in our study area, followed by mountain slopes, roads and forests. Inland lochs and rivers were identified as landscape features that might facilitate gene flow of red deer. Additionally, we explored the effect of choosing arbitrary cell cost values to construct least cost-distance matrices and described a method for improving the selection of cell cost values for a particular landscape feature.

Heliconius wing patterns: an evo-devo model for understanding phenotypic diversity.

Evolutionary Developmental Biology aims for a mechanistic understanding of phenotypic diversity, and present knowledge is largely based on gene expression and interaction patterns from a small number of well-known model organisms. However, our understanding of biological diversification depends on our ability to pinpoint the causes of natural variation at a micro-evolutionary level, and therefore requires the isolation of genetic and developmental variation in a controlled genetic background. The colour patterns of Heliconius butterflies (Nymphalidae: Heliconiinae) provide a rich suite of naturally occurring variants with striking phenotypic diversity and multiple taxonomic levels of variation. Diversification in the genus is well known for its dramatic colour-pattern divergence between races or closely related species, and for Müllerian mimicry convergence between distantly related species, providing a unique system to study the development basis of colour-pattern evolution. A long history of genetic studies has showed that pattern variation is based on allelic combinations at a surprisingly small number of loci, and recent developmental evidence suggests that pattern development in Heliconius is different from the eyespot determination of other butterflies. Fine-scale genetic mapping studies have shown that a shared toolkit of genes is used to produce both convergent and divergent phenotypes. These exciting results and the development of new genomic resources make Heliconius a very promising evo-devo model for the study of adaptive change.

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

Abstract Shared ancestral variation and introgression complicates the reconstruction of phylogenetic relationships among closely related taxa. Here we use overall genomic compatibility as an alternative estimate of species relationships in a group where divergence is rapid and genetic exchange is common. Heliconius heurippa, a butterfly species endemic to Colombia, has a colour pattern genetically intermediate between H. cydno and H. melpomene: its hindwing is nearly indistinguishable from that of H. melpomene and its forewing band is an intermediate phenotype between both species. This observation has lead to the suggestion that the pattern of H. heurippa arose through hybridization. We present a genetic analysis of hybrid compatibility in crosses between the three taxa. Heliconius heurippa x H. cydno and female H. melpomene x male H. heurippa yield fertile and viable F1 hybrids, but male H. melpomene x female H. heurippa crosses yield sterile F1 females. In contrast, Haldane's rule has previously been detected between H. melpomene and H cydno in both directions. Therefore, H. heurippa is most closely related to H. cydno, with some evidence for introgression of genes from H. melpomene. The results are compatible with the hypothesis of a hybrid origin for H. heurippa. In addition, backcrosses using F1 hybrid males provide evidence for a large Z(X)-chromosome effect on sterility and for recessive autosomal sterility factors as predicted by Dominance Theory.

First-generation linkage map of the warningly colored butterfly Heliconius erato.

We report the first genetic linkage map of Heliconius erato, a species that shows remarkable variation in its warningly colored wing patterns. We use crosses between H. erato and its sister species, H. himera, to place two major color pattern genes, D and Cr, on a linkage map containing AFLP, allozyme, microsatellite and single-copy nuclear loci. We identified all 21 linkage groups in an initial genetic screen of 22 progeny from an F1 female x male H. himera family. Of the 229 markers, 87 used to identify linkage groups were also informative in 35 progeny from a sibling backcross (H. himera female x F1 male). With these, and an additional 33 markers informative in the second family, we constructed recombinational maps for 19 of the 21 linkage groups. These maps varied in length from 18.1 to 431.1 centimorgans (cM) and yielded an estimated total length of 2400 cM. The average distance between markers was 23 cM, and eight of the 19 linkage groups, including the sex chromosome (Z) and the chromosome containing the Cr locus, contained two or more codominant anchor loci. Of the three potential candidate genes mapped here, Cubitus interruptus (Ci), Decapentaplegic (Dpp) and Wingless (Wg), only Ci was linked, although loosely, to a known Heliconius color pattern locus. This work is an important first step for constructing a denser genetic map of the H. erato color pattern radiation and for a comparative genomic study of the architecture of mimicry in Heliconius butterflies.