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

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

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

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

Adaptive Introgression across Species Boundaries in Heliconius Butterflies

It is widely documented that hybridisation occurs between many closely related species, but the importance of introgression in adaptive evolution remains unclear, especially in animals. Here, we have examined the role of introgressive hybridisation in transferring adaptations between mimetic Heliconius butterflies, taking advantage of the recent identification of a gene regulating red wing patterns in this genus. By sequencing regions both linked and unlinked to the red colour locus, we found a region that displays an almost perfect genotype by phenotype association across four species, H. melpomene, H. cydno, H. timareta, and H. heurippa. This particular segment is located 70 kb downstream of the red colour specification gene optix, and coalescent analysis indicates repeated introgression of adaptive alleles from H. melpomene into the H. cydno species clade. Our analytical methods complement recent genome scale data for the same region and suggest adaptive introgression has a crucial role in generating adaptive wing colour diversity in this group of butterflies

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

Conservatism and novelty in the genetic architecture of adaptation in Heliconius butterflies

10 pagesInternational audienceUnderstanding the genetic architecture of adaptive traits has been at the centre of modern evolutionary biology since Fisher; however, evaluating how the genetic architecture of ecologically important traits influences their diversification has been hampered by the scarcity of empirical data. Now, high-throughput genomics facilitates the detailed exploration of variation in the genome-to-phenotype map among closely related taxa. Here, we investigate the evolution of wing pattern diversity in Heliconius, a clade of neotropical butterflies that have undergone an adaptive radiation for wing-pattern mimicry and are influenced by distinct selection regimes. Using crosses between natural wing-pattern variants, we used genome-wide restriction site-associated DNA (RAD) genotyping, traditional linkage mapping and multivariate image analysis to study the evolution of the architecture of adaptive variation in two closely related species: Heliconius hecale and H. ismenius. We implemented a new morphometric procedure for the analysis of whole-wing pattern variation, which allows visualising spatial heatmaps of genotype-to-phenotype association for each quantitative trait locus separately. We used the H. melpomene reference genome to fine-map variation for each major wing-patterning region uncovered, evaluated the role of candidate genes and compared genetic architectures across the genus. Our results show that, although the loci responding to mimicry selection are highly conserved between species, their effect size and phenotypic action vary throughout the clade. Multilocus architecture is ancestral and maintained across species under directional selection, whereas the single-locus (supergene) inheritance controlling polymorphism in H. numata appears to have evolved only once. Nevertheless, the conservatism in the wing-patterning toolkit found throughout the genus does not appear to constrain phenotypic evolution towards local adaptive optima

Functional impact and evolution of a novel human polymorphic inversion that disrupts a gene and creates a fusion transcript

Since the discovery of chromosomal inversions almost 100 years ago, how they are maintained in natural populations has been a highly debated issue. One of the hypotheses is that inversion breakpoints could affect genes and modify gene expression levels, although evidence of this came only from laboratory mutants. In humans, a few inversions have been shown to associate with expression differences, but in all cases the molecular causes have remained elusive. Here, we have carried out a complete characterization of a new human polymorphic inversion and determined that it is specific to East Asian populations. In addition, we demonstrate that it disrupts the ZNF257 gene and, through the translocation of the first exon and regulatory sequences, creates a previously nonexistent fusion transcript, which together are associated to expression changes in several other genes. Finally, we investigate the potential evolutionary and phenotypic consequences of the inversion, and suggest that it is probably deleterious. This is therefore the first example of a natural polymorphic inversion that has position effects and creates a new chimeric gene, contributing to answer an old question in evolutionary biology