Comparative genomics of the mimicry switch in Papilio dardanus

The African Mocker Swallowtail, Papilio dardanus, is a textbook example in evolutionary genetics. Classical breeding experiments have shown that wing pattern variation in this polymorphic Batesian mimic is determined by the polyallelic H locus that controls a set of distinct mimetic phenotypes. Using bacterial artificial chromosome (BAC) sequencing, recombination analyses and comparative genomics, we show that H co-segregates with an interval of less than 500 kb that is collinear with two other Lepidoptera genomes and contains 24 genes, including the transcription factor genes engrailed (en) and invected (inv). H is located in a region of conserved gene order, which argues against any role for genomic translocations in the evolution of a hypothesized multi-gene mimicry locus. Natural populations of P. dardanus show significant associations of specific morphs with single nucleotide polymorphisms (SNPs), centred on en. In addition, SNP variation in the H region reveals evidence of non-neutral molecular evolution in the en gene alone. We find evidence for a duplication potentially driving physical constraints on recombination in the lamborni morph. Absence of perfect linkage disequilibrium between different genes in the other morphs suggests that H is limited to nucleotide positions in the regulatory and coding regions of en. Our results therefore support the hypothesis that a single gene underlies wing pattern variation in P. dardanus

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.This work was supported by ERC Starting Grant Stg-243179 (MimEvol) and French Research Agency grant ANR-12-JSV7-0005 (HybEvol) to M.J. C.S. was funded by the Universidad del Rosario FIUR grant QDN-DG001 and COLCIENCIAS (Grant FP44842-5-2017). R.M.M. was funded by a Research Fellowship at King's College, Cambridge, UK. R.M.M. and C.D.J. were funded by an ERC grant (Speciation Genetics 339873)

Wing Patterns in the Mist

Arnaud Martin is with University of California Irvine, Durrell D. Kapan is with University of Hawaii at Manoa, Lawrence E. Gilbert is with UT Austin.The aesthetic appeal of butterfly wing patterns has been costly to their status as a tool of fundamental scientific inquiry. Thus, while mimetic convergence in wing patterns between edible “Batesian” mimics and distasteful models, or between different distasteful “Müllerian” mimics (species that cooperate to educate predators) has long been the subject of genetic analysis [1] and field experiments [2], most biology text books confine mimicry to sections on striking adaptations without applying these examples to broader topics of evolution. Meanwhile, the study of color patterns in animals, often tucked into the same sections of texts, is undergoing a revolution in this age of evo-devo and genomics [3]. Among insect models for studying color pattern, the genus Heliconius is gaining the attention of an ever-widening audience.Biological Sciences, School o

Major improvements to the Heliconius melpomene genome assembly used to confirm 10 chromosome fusion events in 6 million years of butterfly evolution

The Heliconius butterflies are a widely studied adaptive radiation of 46 species spread across Central and South America, several of which are known to hybridize in the wild. Here, we present a substantially improved assembly of the Heliconius melpomene genome, developed using novel methods that should be applicable to improving other genome assemblies produced using short read sequencing. First, we whole-genome-sequenced a pedigree to produce a linkage map incorporating 99% of the genome. Second, we incorporated haplotype scaffolds extensively to produce a more complete haploid version of the draft genome. Third, we incorporated ~20x coverage of Pacific Biosciences sequencing, and scaffolded the haploid genome using an assembly of this long-read sequence. These improvements result in a genome of 795 scaffolds, 275 Mb in length, with an N50 length of 2.1 Mb, an N50 number of 34, and with 99% of the genome placed, and 84% anchored on chromosomes. We use the new genome assembly to confirm that the Heliconius genome underwent 10 chromosome fusions since the split with its sister genus Eueides, over a period of about 6 million yr

Selective sweeps on novel and introgressed variation shape mimicry loci in a butterfly adaptive radiation.

Natural selection leaves distinct signatures in the genome that can reveal the targets and history of adaptive evolution. By analysing high-coverage genome sequence data from 4 major colour pattern loci sampled from nearly 600 individuals in 53 populations, we show pervasive selection on wing patterns in the Heliconius adaptive radiation. The strongest signatures correspond to loci with the greatest phenotypic effects, consistent with visual selection by predators, and are found in colour patterns with geographically restricted distributions. These recent sweeps are similar between co-mimics and indicate colour pattern turn-over events despite strong stabilising selection. Using simulations, we compare sweep signatures expected under classic hard sweeps with those resulting from adaptive introgression, an important aspect of mimicry evolution in Heliconius butterflies. Simulated recipient populations show a distinct 'volcano' pattern with peaks of increased genetic diversity around the selected target, characteristic of sweeps of introgressed variation and consistent with diversity patterns found in some populations. Our genomic data reveal a surprisingly dynamic history of colour pattern selection and co-evolution in this adaptive radiation

Phylogenetic Codivergence Supports Coevolution of Mimetic Heliconius Butterflies

The unpalatable and warning-patterned butterflies _Heliconius erato_ and _Heliconius melpomene_ provide the best studied example of mutualistic Müllerian mimicry, thought – but rarely demonstrated – to promote coevolution. Some of the strongest available evidence for coevolution comes from phylogenetic codivergence, the parallel divergence of ecologically associated lineages. Early evolutionary reconstructions suggested codivergence between mimetic populations of _H. erato_ and _H. melpomene_, and this was initially hailed as the most striking known case of coevolution. However, subsequent molecular phylogenetic analyses found discrepancies in phylogenetic branching patterns and timing (topological and temporal incongruence) that argued against codivergence. We present the first explicit cophylogenetic test of codivergence between mimetic populations of _H. erato_ and _H. melpomene_, and re-examine the timing of these radiations. We find statistically significant topological congruence between multilocus coalescent population phylogenies of _H. erato_ and _H. melpomene_, supporting repeated codivergence of mimetic populations. Divergence time estimates, based on a Bayesian coalescent model, suggest that the evolutionary radiations of _H. erato_ and _H. melpomene_ occurred over the same time period, and are compatible with a series of temporally congruent codivergence events. This evidence supports a history of reciprocal coevolution between Müllerian co-mimics characterised by phylogenetic codivergence and parallel phenotypic change

Genomic hotspots for adaptation: the population genetics of Müllerian mimicry in Heliconius erato

This is the final version of the article. Available from the publisher via the DOI in this record.Wing pattern evolution in Heliconius butterflies provides some of the most striking examples of adaptation by natural selection. The genes controlling pattern variation are classic examples of Mendelian loci of large effect, where allelic variation causes large and discrete phenotypic changes and is responsible for both convergent and highly divergent wing pattern evolution across the genus. We characterize nucleotide variation, genotype-by-phenotype associations, linkage disequilibrium (LD), and candidate gene expression patterns across two unlinked genomic intervals that control yellow and red wing pattern variation among mimetic forms of Heliconius erato. Despite very strong natural selection on color pattern, we see neither a strong reduction in genetic diversity nor evidence for extended LD across either patterning interval. This observation highlights the extent that recombination can erase the signature of selection in natural populations and is consistent with the hypothesis that either the adaptive radiation or the alleles controlling it are quite old. However, across both patterning intervals we identified SNPs clustered in several coding regions that were strongly associated with color pattern phenotype. Interestingly, coding regions with associated SNPs were widely separated, suggesting that color pattern alleles may be composed of multiple functional sites, conforming to previous descriptions of these loci as "supergenes." Examination of gene expression levels of genes flanking these regions in both H. erato and its co-mimic, H. melpomene, implicate a gene with high sequence similarity to a kinesin as playing a key role in modulating pattern and provides convincing evidence for parallel changes in gene regulation across co-mimetic lineages. The complex genetic architecture at these color pattern loci stands in marked contrast to the single casual mutations often identified in genetic studies of adaptation, but may be more indicative of the type of genetic changes responsible for much of the adaptive variation found in natural populations.Funding: Funding for this study was provided by National Science Foundation grants to WOM (DEB-0715096 and IBN-0344705) and BAC (DEB-0513424). Funding for work on H. melpomene came from a BBSRC grant to CDJ and RHf-C (011845). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript

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.RMM is funded by a Junior Research Fellowship at King’s College, Cambridge. KMK is supported by the Balfour Studentship, University of Cambridge, SHMa by a Research Fellowship at St John's College, Cambridge, and SHMo by a Research Fellowship from the Royal Commission for the Exhibition of 1851. Our work on Heliconius has been additionally supported by the Agence Nationale de la Recherche (France), the Biology and Biotechnology Research Council (UK), the British Ecological Society, the European Research Council, the Natural Environment Research Council (UK), and the Smithsonian Tropical Research Institute.This is the author accepted manuscript. The final version is available from Wiley via http://dx.doi.org/10.1111/jeb.1267

Optimal-Foraging Predator Favors Commensalistic Batesian Mimicry

BACKGROUND:Mimicry, in which one prey species (the Mimic) imitates the aposematic signals of another prey (the Model) to deceive their predators, has attracted the general interest of evolutionary biologists. Predator psychology, especially how the predator learns and forgets, has recently been recognized as an important factor in a predator-prey system. This idea is supported by both theoretical and experimental evidence, but is also the source of a good deal of controversy because of its novel prediction that in a Model/Mimic relationship even a moderately unpalatable Mimic increases the risk of the Model (quasi-Batesian mimicry). METHODOLOGY/PRINCIPAL FINDINGS:We developed a psychology-based Monte Carlo model simulation of mimicry that incorporates a "Pavlovian" predator that practices an optimal foraging strategy, and examined how various ecological and psychological factors affect the relationships between a Model prey species and its Mimic. The behavior of the predator in our model is consistent with that reported by experimental studies, but our simulation's predictions differed markedly from those of previous models of mimicry because a more abundant Mimic did not increase the predation risk of the Model when alternative prey were abundant. Moreover, a quasi-Batesian relationship emerges only when no or very few alternative prey items were available. Therefore, the availability of alternative prey rather than the precise method of predator learning critically determines the relationship between Model and Mimic. Moreover, the predation risk to the Model and Mimic is determined by the absolute density of the Model rather than by its density relative to that of the Mimic. CONCLUSIONS/SIGNIFICANCE:Although these predictions are counterintuitive, they can explain various kinds of data that have been offered in support of competitive theories. Our model results suggest that to understand mimicry in nature it is important to consider the likely presence of alternative prey and the possibility that predation pressure is not constant

Genomic Hotspots for Adaptation: The Population Genetics of Müllerian Mimicry in the Heliconius melpomene Clade

Wing patterning in Heliconius butterflies is a longstanding example of both Müllerian mimicry and phenotypic radiation under strong natural selection. The loci controlling such patterns are “hotspots” for adaptive evolution with great allelic diversity across different species in the genus. We characterise nucleotide variation, genotype-by-phenotype associations, linkage disequilibrium, and candidate gene expression at two loci and across multiple hybrid zones in Heliconius melpomene and relatives. Alleles at HmB control the presence or absence of the red forewing band, while alleles at HmYb control the yellow hindwing bar. Across HmYb two regions, separated by ∼100 kb, show significant genotype-by-phenotype associations that are replicated across independent hybrid zones. In contrast, at HmB a single peak of association indicates the likely position of functional sites at three genes, encoding a kinesin, a G-protein coupled receptor, and an mRNA splicing factor. At both HmYb and HmB there is evidence for enhanced linkage disequilibrium (LD) between associated sites separated by up to 14 kb, suggesting that multiple sites are under selection. However, there was no evidence for reduced variation or deviations from neutrality that might indicate a recent selective sweep, consistent with these alleles being relatively old. Of the three genes showing an association with the HmB locus, the kinesin shows differences in wing disc expression between races that are replicated in the co-mimic, Heliconius erato, providing striking evidence for parallel changes in gene expression between Müllerian co-mimics. Wing patterning loci in Heliconius melpomene therefore show a haplotype structure maintained by selection, but no evidence for a recent selective sweep. The complex genetic pattern contrasts with the simple genetic basis of many adaptive traits studied previously, but may provide a better model for most adaptation in natural populations that has arisen over millions rather than tens of years