Mimicry--whereby warning signals in different species evolve to look similar--has long served as a paradigm of convergent evolution. Little is known, however, about the genes that underlie the evolution of mimetic phenotypes or to what extent the same or different genes drive such convergence. Here, we characterize one of the major genes responsible for mimetic wing pattern evolution in Heliconius butterflies. Mapping, gene expression, and population genetic work all identify a single gene, optix, that controls extreme red wing pattern variation across multiple species of Heliconius. Our results show that the cis-regulatory evolution of a single transcription factor can repeatedly drive the convergent evolution of complex color patterns in distantly related species, thus blurring the distinction between convergence and homology.

One of the most striking examples of sexual dimorphism is sex-limited mimicry in butterflies, a phenomenon in which one sex--usually the female--mimics a toxic model species, whereas the other sex displays a different wing pattern. Sex-limited mimicry is phylogenetically widespread in the swallowtail butterfly genus Papilio, in which it is often associated with female mimetic polymorphism. In multiple polymorphic species, the entire wing pattern phenotype is controlled by a single Mendelian 'supergene'. Although theoretical work has explored the evolutionary dynamics of supergene mimicry, there are almost no empirical data that address the critical issue of what a mimicry supergene actually is at a functional level. Using an integrative approach combining genetic and association mapping, transcriptome and genome sequencing, and gene expression analyses, we show that a single gene, doublesex, controls supergene mimicry in Papilio polytes. This is in contrast to the long-held view that supergenes are likely to be controlled by a tightly linked cluster of loci. Analysis of gene expression and DNA sequence variation indicates that isoform expression differences contribute to the functional differences between dsx mimicry alleles, and protein sequence evolution may also have a role. Our results combine elements from different hypotheses for the identity of supergenes, showing that a single gene can switch the entire wing pattern among mimicry phenotypes but may require multiple, tightly linked mutations to do so.

We used 20 de novo genome assemblies to probe the speciation history and architecture of gene flow in rapidly radiating Heliconius butterflies. Our tests to distinguish incomplete lineage sorting from introgression indicate that gene flow has obscured several ancient phylogenetic relationships in this group over large swathes of the genome. Introgressed loci are underrepresented in low-recombination and gene-rich regions, consistent with the purging of foreign alleles more tightly linked to incompatibility loci. Here, we identify a hitherto unknown inversion that traps a color pattern switch locus. We infer that this inversion was transferred between lineages by introgression and is convergent with a similar rearrangement in another part of the genus. These multiple de novo genome sequences enable improved understanding of the importance of introgression and selective processes in adaptive radiation.

Sexual isolation is a critical form of reproductive isolation in the early stages of animal speciation, yet little is known about the genetic basis of divergent mate preferences and preference cues in young species. Heliconius butterflies, well known for their diversity of wing color patterns, mate assortatively as a result of divergence in male preference for wing patterns. Here we show that the specific cue used by Heliconius cydno and Heliconius pachinus males to recognize conspecific females is the color of patches on the wings. In addition, male mate preference segregates with forewing color in hybrids, indicating a genetic association between the loci responsible for preference and preference cue. Quantitative trait locus mapping places a preference locus coincident with the locus that determines forewing color, which itself is perfectly linked to the wing patterning candidate gene, wingless. Furthermore, yellow-colored males of the polymorphic race H. cydno alithea prefer to court yellow females, indicating that wing color and color preference are controlled by loci that are located in an inversion or are pleiotropic effects of a single locus. Tight genetic associations between preference and preference cue, although rare, make divergence and speciation particularly likely because the effects of natural and sexual selection on one trait are transferred to the other, leading to the coordinated evolution of mate recognition. This effect of linkage on divergence is especially important in Heliconius because differentiation of wing color patterns in the genus has been driven and maintained by natural selection for Mü llerian mimicry.Heliconius ͉ Lepidoptera ͉ sexual isolation ͉ speciation

The monarch butterfly, Danaus plexippus, is famous for its spectacular annual migration across North America, recent worldwide dispersal, and orange warning coloration. Despite decades of study and broad public interest, we know little about the genetic basis of these hallmark traits. By sequencing 101 monarch genomes from around the globe, we uncover the history of the monarch's evolutionary origin and global dispersal, characterize the genes and pathways associated with migratory behavior, and identify the discrete genetic basis of warning coloration. The results rewrite our understanding of this classic system, showing that D. plexippus was ancestrally migratory and dispersed out of North America to occupy its broad distribution. We find the strongest signatures of selection associated with migration center on flight muscle function, resulting in greater flight efficiency among migratory monarchs, and that variation in monarch warning coloration is controlled by a single myosin gene not previously implicated in insect pigmentation.

Although animals display a rich variety of shapes and patterns, the genetic changes that explain how complex forms arise are still unclear. Here we take advantage of the extensive diversity of Heliconius butterflies to identify a gene that causes adaptive variation of black wing patterns within and between species. Linkage mapping in two species groups, gene-expression analysis in seven species, and pharmacological treatments all indicate that cis-regulatory evolution of the WntA ligand underpins discrete changes in color pattern features across the Heliconius genus. These results illustrate how the direct modulation of morphogen sources can generate a wide array of unique morphologies, thus providing a link between natural genetic variation, pattern formation, and adaptation.Mü llerian mimicry | Wnt pathway | Mendelian genetics | evolutionary-developmental biology

We here pioneer a low-cost assembly strategy for 20 Heliconiini genomes to characterize the evolutionary history of the rapidly radiating genus Heliconius. A bifurcating tree provides a poor fit to the data, and we therefore explore a reticulate phylogeny for Heliconius. We probe the genomic architecture of gene flow, and develop a new method to distinguish incomplete lineage sorting from introgression. We find that most loci with non-canonical histories arose through introgression, and are strongly underrepresented in regions of low recombination and high gene density. This is expected if introgressed alleles are more likely to be purged in such regions due to tighter linkage with incompatibility loci. Finally, we identify a hitherto unrecognized inversion, and show it is a convergent structural rearrangement that captures a known color pattern switch locus within the genus. Our multi-genome assembly approach enables an improved understanding of adaptive radiation.

The mimetic butterflies Heliconius erato and Heliconius melpomene have undergone parallel radiations to form a near-identical patchwork of over 20 different wing-pattern races across the Neotropics. Previous molecular phylogenetic work on these radiations has suggested that similar but geographically disjunct color patterns arose multiple times independently in each species. The neutral markers used in these studies, however, can move freely across color pattern boundaries, and therefore might not represent the history of the adaptive traits as accurately as markers linked to color pattern genes. To assess the evolutionary histories across different loci, we compared relationships among races within H. erato and within H. melpomene using a series of unlinked genes, genes linked to color pattern loci, and optix, a gene recently shown to control red colorpattern variation. We found that although unlinked genes partition populations by geographic region, optix had a different history, structuring lineages by red color patterns and supporting a single origin of red-rayed patterns within each species. Genes closely linked (80-250 kb) to optix exhibited only weak associations with color pattern. This study empirically demonstrates the necessity of examining phenotype-determining genomic regions to understand the history of adaptive change in rapidly radiating lineages. With these refined relationships, we resolve a long-standing debate about the origins of the races within each species, supporting the hypothesis that the red-rayed Amazonian pattern evolved recently and expanded, causing disjunctions of more ancestral patterns.Müllerian mimicry | population genetics | phylogeography R esearchers typically rely on neutrally evolving loci to generate a phylogenetic and population genetic history of adaptive divergence. The rationale is that these markers provide an unbiased view of the relationships among divergent phenotypes and a better understanding of the evolutionary processes generating variation. However, the genome is a complicated mosaic shaped by an interplay of mutation, drift, selection, and recombination. Recombination allows different regions of the genome to experience alternative restrictions to gene flow, and thus develop different evolutionary trajectories. The closer a genetic marker is to the alleles responsible for adaptive differences, the more likely that it will trace the history of phenotypic change.Understanding how phenotypic variation is generated in nature is greatly enhanced by studying groups that are actively undergoing diversification. By deciphering the history of such diverse phenotypes we gain a clearer understanding of the evolutionary process, including the tempo and mode of phenotypic change. Heliconius butterflies present one of the most striking examples of a recent phenotypic radiation. The 40 species in the genus exhibit hundreds of wing patterns that are involved in Müllerian mimicry complexes, where distasteful species converge on a shared warning signal to avoid predation. Th...