Understanding the origin and maintenance of adaptive phenotypic novelty is a central goal of evolutionary biology. However, both hybridization and incomplete lineage sorting can lead to genealogical discordance between the regions of the genome underlying adaptive traits and the remainder of the genome, decoupling inferences about character evolution from population history. Here, to disentangle these effects, we investigated the evolutionary origins and maintenance of Batesian mimicry between North American admiral butterflies (Limenitis arthemis) and their chemically defended model (Battus philenor) using a combination of de novo genome sequencing, whole-genome resequencing, and statistical introgression mapping. Our results suggest that balancing selection, arising from geographic variation in the presence or absence of the unpalatable model, has maintained two deeply divergent color patterning haplotypes that have been repeatedly sieved among distinct mimetic and nonmimetic lineages of Limenitis via introgressive hybridization.

Despite more than a century of biological research on the evolution and maintenance of mimetic signals, the relative frequencies of models and mimics necessary to establish and maintain Batesian mimicry in natural populations remain understudied. Here we investigate the frequency-dependent dynamics of imperfect Batesian mimicry, using predation experiments involving artificial butterfly models. We use two geographically distinct populations of Adelpha butterflies that vary in their relative frequencies of a putatively defended model (Adelpha iphiclus) and Batesian mimic (Adelpha serpa). We found that in Costa Rica, where both species share similar abundances, Batesian mimicry breaks down, and predators more readily attack artificial butterfly models of the presumed mimic, A. serpa By contrast, in Ecuador, where A. iphiclus (model) is significantly more abundant than A. serpa (mimic), both species are equally protected from predation. Our results provide compelling experimental evidence that imperfect Batesian mimicry is frequency-dependent on the relative abundance of models and mimics in natural populations, and contribute to the growing body of evidence that complex dynamics, such as seasonality or the availability of alternative prey, influence the evolution of mimetic traits.

Rapid diversification is often associated with morphological or ecological adaptations that allow organisms to radiate into novel niches. Neotropical Adelpha butterflies, which comprise over 200 species and subspecies, are characterized by extraordinary breadth in host plant use and wing colour patterns compared to their closest relatives. To examine the relationship between phenotypic and species diversification, we reconstructed the phylogenetic history of Adelpha and its temperate sister genus Limenitis using genomewide restriction-site-associated DNA (RAD) sequencing. Despite a declining fraction of shared markers with increasing evolutionary distance, the RAD-Seq data consistently generated well-supported trees using a variety of phylogenetic methods. These well-resolved phylogenies allow the identification of an ecologically important relationship with a toxic host plant family, as well as the confirmation of widespread, convergent wing pattern mimicry throughout the genus. Taken together, our results support the hypothesis that evolutionary innovations in both larvae and adults have permitted the colonization of novel host plants and fuelled adaptive diversification within this large butterfly radiation.