Ancient homology underlies adaptive mimetic diversity across butterflies

A micro-RNA is the effector gene of a classic evolutionary hotspot locus

In Lepidoptera (butterflies and moths), the genomic region around the gene cortex is a 'hotspot' locus, repeatedly used to generate intraspecific melanic wing color polymorphisms across 100-million-years of evolution. However, the identity of the effector gene regulating melanic wing color within this locus remains unknown. Here, we show that none of the four candidate protein-coding genes within this locus, including cortex, serve as major effectors. Instead, a micro-RNA (miRNA), mir-193, serves as the major effector across three deeply diverged lineages of butterflies, and its function is conserved in Drosophila. In Lepidoptera, mir-193 is derived from a gigantic long non-coding RNA, ivory, and it functions by directly repressing multiple pigmentation genes. We show that a miRNA can drive repeated instances of adaptive evolution in animals.

Potential and progress of studying mountain biodiversity by means of butterfly genetics and genomics

Mountains are rich in biodiversity, and butterflies are species-rich and have a good ecological and evolutionary research foundation. This review addresses the potential and progress of studying mountain biodiversity using butterflies as a model. We discuss the uniqueness of mountain ecosystems, factors influencing the distribution of mountain butterflies, representative genetic and evolutionary models in butterfly research, and evolutionary studies of mountain biodiversity involving butterfly genetics and genomics. Finally, we demonstrate the necessity of studying mountain butterflies and propose future perspectives. This review provides insights for studying the biodiversity of mountain butterflies as well as a summary of research methods for reference.

The diversification of butterfly wing patterns: progress and prospects

Butterfly wings display rich phenotypic diversity and are associated with complex biological functions, thus serving as an important evolutionary system to address the genetic basis and evolution of phenotypic diversification. We review recent butterfly studies that revealed complex functions underlying diversified wing patterns and describe the genetic and environmental factors involved in wing pattern determinations. These factors lead to inter-specific divergence, genetic polymorphism, and phenotypic plasticity, which in many cases are decided by several key genes. We also summarize the research advances on gene co-option as an important origin of functional complexity and evolutionary novelty. These findings reveal a pattern of evolutionary innovation within a constrained developmental framework during butterfly wing morphogenesis, but further research is required to gain a systematic and comprehensive understanding.

Frizzled2 receives WntA signaling during butterfly wing pattern formation

Butterfly color patterns provide visible and biodiverse phenotypic readouts of the patterning processes. Although the secreted ligand WntA has been shown to instruct the color pattern formation in butterflies, its mode of reception remains elusive. Butterfly genomes encode four homologs of the Frizzled-family of Wnt receptors. Here, we show that CRISPR mosaic knockouts of frizzled2 (fz2) phenocopy the color pattern effects of WntA loss of function in multiple nymphalids. Whereas WntA mosaic clones result in intermediate patterns of reduced size, fz2 clones are cell-autonomous, consistent with a morphogen function. Shifts in expression of WntA and fz2 in WntA crispant pupae show that they are under positive and negative feedback, respectively. Fz1 is required for Wnt-independent planar cell polarity in the wing epithelium. Fz3 and Fz4 show phenotypes consistent with Wnt competitive-antagonist functions in vein formation (Fz3 and Fz4), wing margin specification (Fz3), and color patterning in the Discalis and Marginal Band Systems (Fz4). Overall, these data show that the WntA/Frizzled2 morphogen-receptor pair forms a signaling axis that instructs butterfly color patterning and shed light on the functional diversity of insect Frizzled receptors.

A histological analysis of coloration in the Peruvian mimic poison frog (Ranitomeya imitator)

Aposematism continues to be a phenomenon of central interest in evolutionary biology. The life history of the mimic poison frog, Ranitomeya imitator, relies heavily on aposematism. In order for aposematic signals to be effective, predators must be able to learn to avoid the associated phenotype. However, in R. imitator, aposematism is associated with four different color phenotypes that mimic a complex of congeneric species occurring across the mimic frog's geographic range. Investigations of the underlying mechanics of color production in these frogs can provide insights into how and why these different morphs evolved. We used histological samples to examine divergence in the color production mechanisms used by R. imitator to produce effective aposematic signals across its geographic range. We measured the coverage of melanophores and xanthophores (the area covered by chromatophores divided by total area of the skin section) in each color morph. We find that morphs that produce orange skin exhibit a higher coverage of xanthophores and lower coverage of melanophores than those that produce yellow skin. In turn, morphs that produce yellow skin exhibit a higher coverage of xanthophores and lower coverage of melanophores than those that produce green skin. Generally, across the morphs, a high ratio of xanthophores to melanophores is associated with colors of brighter spectral reflectance. Together, our results contribute to the understanding of color production in amphibians and document divergence in the histology of a species that is subject to divergent selection associated with aposematism.

The influence of transposable elements on animal colouration

Transposable elements (TEs) are mobile genetic sequences present within host genomes. TEs can contribute to the evolution of host traits, since transposition is mutagenic and TEs often contain host regulatory and protein coding sequences. We review cases where TEs influence animal colouration, reporting major patterns and outstanding questions. TE-induced colouration phenotypes typically arise via introduction of novel regulatory sequences and splice sites, affecting pigment cell development or pigment synthesis. We discuss if particular TE types may be more frequently involved in the evolution of colour variation in animals, given that examples involving long terminal repeat (LTR) elements appear to dominate. Currently, examples of TE-induced colouration phenotypes in animals mainly concern model and domesticated insect and mammal species. However, several influential recent examples, coupled with increases in genome sequencing, suggest cases reported from wild species will increase considerably.

Widespread chromosomal rearrangements preceded genetic divergence in a monitor lizard, Varanus acanthurus (Varanidae)

Chromosomal rearrangements are often associated with local adaptation and speciation because they suppress recombination, and as a result, rearrangements have been implicated in disrupting gene flow. Although there is strong evidence to suggest that chromosome rearrangements are a factor in genetic isolation of divergent populations, the underlying mechanism remains elusive. Here, we applied an integrative cytogenetics and genomics approach testing whether chromosomal rearrangements are the initial process, or a consequence, of population divergence in the dwarf goanna, Varanus acanthurus. Specifically, we tested whether chromosome rearrangements are indicators of genetic barriers that can be used to identify divergent populations by looking at gene flow within and between populations with rearrangements. We found that gene flow was present between individuals with chromosome rearrangements within populations, but there was no gene flow between populations that had similar chromosome rearrangements. Moreover, we identified a correlation between reduced genetic variation in populations with a higher frequency of homozygous submetacentric individuals. These findings suggest that chromosomal rearrangements were widespread prior to divergence, and because we found populations with higher frequencies of submetacentric chromosomes were associated with lower genetic diversity, this could indicate that polymorphisms within populations are early indicators of genetic drift.

How butterfly wings got their pattern

Gene regulatory elements play a crucial role in the pattern formation of butterfly wings.

Deep cis-regulatory homology of the butterfly wing pattern ground plan

Butterfly wing patterns derive from a deeply conserved developmental ground plan yet are diverse and evolve rapidly. It is poorly understood how gene regulatory architectures can accommodate both deep homology and adaptive change. To address this, we characterized the cis-regulatory evolution of the color pattern gene WntA in nymphalid butterflies. Comparative assay for transposase-accessible chromatin using sequencing (ATAC-seq) and in vivo deletions spanning 46 cis-regulatory elements across five species revealed deep homology of ground plan-determining sequences, except in monarch butterflies. Furthermore, noncoding deletions displayed both positive and negative regulatory effects that were often broad in nature. Our results provide little support for models predicting rapid enhancer turnover and suggest that deeply ancestral, multifunctional noncoding elements can underlie rapidly evolving trait systems.

Aculeate Hymenopterans as Aposematic and Mimetic Models

Aposematism and mimicry are complex phenomena which have been studied extensively; however, much of our knowledge comes from just a few focal groups, especially butterflies. Aposematic species combine a warning signal with a secondary defense that reduces their profitability as prey. Aculeate hymenopterans are an extremely diverse lineage defined by the modification of the ovipositor into a stinger which represents a potent defense against predators. Aculeates are often brightly colored and broadly mimicked by members of other arthropod groups including Diptera, Lepidoptera, Coleoptera, and Araneae. However, aculeates are surprisingly understudied as aposematic and mimetic model organisms. Recent studies have described novel pigments contributing to warning coloration in insects and identified changes in cis-regulatory elements as potential drivers of color pattern evolution. Many biotic and abiotic factors contribute to the evolution and maintenance of conspicuous color patterns. Predator distribution and diversity seem to influence the phenotypic diversity of aposematic velvet ants while studies on bumble bees underscore the importance of intermediate mimetic phenotypes in transition zones between putative mimicry rings. Aculeate hymenopterans are attractive models for studying sex-based intraspecific mimicry as male aculeates lack the defense conferred by the females’ stinger. In some species, evolution of male and female color patterns appears to be decoupled. Future studies on aposematic aculeates and their associated mimics hold great promise for unraveling outstanding questions about the evolution of conspicuous color patterns and the factors which determine the composition and distribution of mimetic communities.

Ecological basis and genetic architecture of crypsis polymorphism in the desert clicker grasshopper (Ligurotettix coquilletti)

Color polymorphic species can offer exceptional insight into the ecology and genetics of adaptation. Although the genetic architecture of animal coloration is diverse, many color polymorphisms are associated with large structural variants and maintained by biotic interactions. Grasshoppers are notably polymorphic in both color and karyotype, which makes them excellent models for understanding the ecological drivers and genetic underpinnings of color variation. Banded and uniform morphs of the desert clicker grasshopper (Ligurotettix coquilletti) are found across the western deserts of North America. To address the hypothesis that predation maintains local color polymorphism and shapes regional crypsis variation, we surveyed morph frequencies and tested for covariation with two predation environments. Morphs coexisted at intermediate frequencies at most sites, consistent with local balancing selection. Morph frequencies covaried with the appearance of desert substrate-an environment used only by females-suggesting that ground-foraging predators are major agents of selection on crypsis. We next addressed the hypothesized link between morph variation and genome structure. To do so, we designed an approach for detecting inversions and indels using only RADseq data. The banded morph was perfectly correlated with a large putative indel. Remarkably, indel dominance differed among populations, a rare example of dominance evolution in nature.

Convergent flight morphology among Müllerian mimic mutualists

Müllerian mimicry involves a signal mutualism between prey species, shaped by visually hunting predators, and recent work has emphasized the importance of color pattern. Predators respond to more than color pattern, however, and other traits are much less studied. This article examines the hypothesis of convergent evolution in flight-related morphology among eight mimicry complexes composed of 51 butterfly species (Nymphalidae, Danainae, Ithomiini) from a single community in Ecuador. Phylogenetic comparative analyses of 14 variables indicated strong morphological differences between mimicry complexes belonging to three clusters of morphological space ("large yellow transparent," "tiger," and "transparent"), not the eight predicted based on color pattern alone. Analyses found convergence within mimicry complexes, convergence between mimicry complexes within morphospace clusters, and divergence between mimicry complexes from different morphospace clusters. These three clusters differed in size, and body and wing shape, predicting that flight biomechanics also converge (i.e., locomotor mimicry). Potential constraints on evolution of morphological mimicry related to predator discrimination, and evolutionary rates, likely e xplain why flight-related morphology differences were limited to three clusters of morphological space. Finally, the added complexity that flight-related morphology brings to signals between predator and prey indicates that evolutionary switches in color pattern are not all equally likely, potentially limiting the evolution of color patterns if they do not match morphology.

Reinventing the wheel? Reassessing the roles of gene flow, sorting and convergence in repeated evolution

Biologists have long been intrigued by apparently predictable and repetitive evolutionary trajectories inferred across a variety of lineages and systems. In recent years, high-throughput sequencing analyses have started to transform our understanding of such repetitive shifts. While researchers have traditionally categorized such shifts as either "convergent" or "parallel," based on relatedness of the lineages involved, emerging genomic insights provide an opportunity to better describe the actual evolutionary mechanisms at play. A synthesis of recent genomic analyses confirms that convergence is the predominant driver of repetitive evolution among species, whereas repeated sorting of standing variation is the major driver of repeated shifts within species. However, emerging data reveal numerous notable exceptions to these expectations, with recent examples of de novo mutations underpinning convergent shifts among even very closely related lineages, while repetitive sorting processes have occurred among even deeply divergent taxa, sometimes via introgression. A number of very recent analyses have found evidence for both processes occurring on different scales within taxa. We suggest that the relative importance of convergent versus sorting processes depends on the interplay between gene flow among populations, and phylogenetic relatedness of the lineages involved.

Porcupine/Wntless-dependent trafficking of the conserved WntA ligand in butterflies

Wnt ligands are key signaling molecules in animals, but little is known about the evolutionary dynamics and mode of action of the WntA orthologs, which are not present in the vertebrates or in Drosophila. Here we show that the WntA subfamily evolved at the base of the Bilateria + Cnidaria clade, and conserved the thumb region and Ser209 acylation site present in most other Wnts, suggesting WntA requires the core Wnt secretory pathway. WntA proteins are distinguishable from other Wnts by a synapomorphic Iso/Val/Ala216 amino-acid residue that replaces the otherwise ubiquitous Thr216 position. WntA embryonic expression is conserved between beetles and butterflies, suggesting functionality, but the WntA gene was lost three times within arthropods, in podoplean copepods, in the cyclorrhaphan fly radiation, and in ensiferan crickets and katydids. Finally, CRISPR mosaic knockouts (KOs) of porcupine and wntless phenocopied the pattern-specific effects of WntA KOs in the wings of Vanessa cardui butterflies. These results highlight the molecular conservation of the WntA protein across invertebrates, and imply it functions as a typical Wnt ligand that is acylated and secreted through the Porcupine/Wntless secretory pathway.

Disentangling Population History and Character Evolution among Hybridizing Lineages

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.

The evolution and genetics of sexually dimorphic 'dual' mimicry in the butterfly Elymnias hypermnestra

Sexual dimorphism is a major component of morphological variation across the tree of life, but the mechanisms underlying phenotypic differences between sexes of a single species are poorly understood. We examined the population genomics and biogeography of the common palmfly Elymnias hypermnestra, a dual mimic in which female wing colour patterns are either dark brown (melanic) or bright orange, mimicking toxic Euploea and Danaus species, respectively. As males always have a melanic wing colour pattern, this makes E. hypermnestra a fascinating model organism in which populations vary in sexual dimorphism. Population structure analysis revealed that there were three genetically distinct E. hypermnestra populations, which we further validated by creating a phylogenomic species tree and inferring historical barriers to gene flow. This species tree demonstrated that multiple lineages with orange females do not form a monophyletic group, and the same is true of clades with melanic females. We identified two single nucleotide polymorphisms (SNPs) near the colour patterning gene WntA that were significantly associated with the female colour pattern polymorphism, suggesting that this gene affects sexual dimorphism. Given WntA's role in colour patterning across Nymphalidae, E. hypermnestra females demonstrate the repeatability of the evolution of sexual dimorphism.

Deep Convergence, Shared Ancestry, and Evolutionary Novelty in the Genetic Architecture of Heliconius Mimicry

Convergent evolution can occur through different genetic mechanisms in different species. It is now clear that convergence at the genetic level is also widespread, and can be caused by either (i) parallel genetic evolution, where independently evolved convergent mutations arise in different populations or species, or (ii) collateral evolution in which shared ancestry results from either ancestral polymorphism or introgression among taxa. The adaptive radiation of Heliconius butterflies shows color pattern variation within species, as well as mimetic convergence between species. Using comparisons from across multiple hybrid zones, we use signals of shared ancestry to identify and refine multiple putative regulatory elements in Heliconius melpomene and its comimics, Heliconius elevatus and Heliconius besckei, around three known major color patterning genes: optix, WntA, and cortex While we find that convergence between H. melpomene and H. elevatus is caused by a complex history of collateral evolution via introgression in the Amazon, convergence between these species in the Guianas appears to have evolved independently. Thus, we find adaptive convergent genetic evolution to be a key driver of regulatory changes that lead to rapid phenotypic changes. Furthermore, we uncover evidence of parallel genetic evolution at some loci around optix and WntA in H. melpomene and its distant comimic Heliconius erato Ultimately, we show that all three of convergence, conservation, and novelty underlie the modular architecture of Heliconius color pattern mimicry.

Perfect mimicry between Heliconius butterflies is constrained by genetics and development

Müllerian mimicry strongly exemplifies the power of natural selection. However, the exact measure of such adaptive phenotypic convergence and the possible causes of its imperfection often remain unidentified. Here, we first quantify wing colour pattern differences in the forewing region of 14 co-mimetic colour pattern morphs of the butterfly species Heliconius erato and Heliconius melpomene and measure the extent to which mimicking colour pattern morphs are not perfectly identical. Next, using gene-editing CRISPR/Cas9 KO experiments of the gene WntA, which has been mapped to colour pattern diversity in these butterflies, we explore the exact areas of the wings in which WntA affects colour pattern formation differently in H. erato and H. melpomene. We find that, while the relative size of the forewing pattern is generally nearly identical between co-mimics, the CRISPR/Cas9 KO results highlight divergent boundaries in the wing that prevent the co-mimics from achieving perfect mimicry. We suggest that this mismatch may be explained by divergence in the gene regulatory network that defines wing colour patterning in both species, thus constraining morphological evolution even between closely related species.

Gene Flow in the Müllerian Mimicry Ring of a Poisonous Papuan Songbird Clade (Pitohui; Aves)

Müllerian mimicry rings are remarkable symbiotic species assemblages in which multiple members share a similar phenotype. However, their evolutionary origin remains poorly understood. Although gene flow among species has been shown to generate mimetic patterns in some Heliconius butterflies, mimicry is believed to be due to true convergence without gene flow in many other cases. We investigated the evolutionary history of multiple members of a passerine mimicry ring in the poisonous Papuan pitohuis. Previous phylogenetic evidence indicates that the aposematic coloration shared by many, but not all, members of this genus is ancestral and has only been retained by members of the mimicry ring. Using a newly assembled genome and thousands of genomic DNA markers, we demonstrate gene flow from the hooded pitohui (Pitohui dichrous) into the southern variable pitohui (Pitohui uropygialis), consistent with shared patterns of aposematic coloration. The vicinity of putatively introgressed loci is significantly enriched for genes that are important in melanin pigment expression and toxin resistance, suggesting that gene flow may have been instrumental in the sharing of plumage patterns and toxicity. These results indicate that interspecies gene flow may be a more general mechanism in generating mimicry rings than hitherto appreciated.

Butterfly Mimicry Polymorphisms Highlight Phylogenetic Limits of Gene Reuse in the Evolution of Diverse Adaptations

Some genes have repeatedly been found to control diverse adaptations in a wide variety of organisms. Such gene reuse reveals not only the diversity of phenotypes these unique genes control but also the composition of developmental gene networks and the genetic routes available to and taken by organisms during adaptation. However, the causes of gene reuse remain unclear. A small number of large-effect Mendelian loci control a huge diversity of mimetic butterfly wing color patterns, but reasons for their reuse are difficult to identify because the genetic basis of mimicry has primarily been studied in two systems with correlated factors: female-limited Batesian mimicry in Papilio swallowtails (Papilionidae) and non-sex-limited Müllerian mimicry in Heliconius longwings (Nymphalidae). Here, we break the correlation between phylogenetic relationship and sex-limited mimicry by identifying loci controlling female-limited mimicry polymorphism Hypolimnas misippus (Nymphalidae) and non-sex-limited mimicry polymorphism in Papilio clytia (Papilionidae). The Papilio clytia polymorphism is controlled by the genome region containing the gene cortex, the classic P supergene in Heliconius numata, and loci controlling color pattern variation across Lepidoptera. In contrast, female-limited mimicry polymorphism in Hypolimnas misippus is associated with a locus not previously implicated in color patterning. Thus, although many species repeatedly converged on cortex and its neighboring genes over 120 My of evolution of diverse color patterns, female-limited mimicry polymorphisms each evolved using a different gene. Our results support conclusions that gene reuse occurs mainly within ∼10 My and highlight the puzzling diversity of genes controlling seemingly complex female-limited mimicry polymorphisms.

Studying convergent evolution to relate genotype to behavioral phenotype

Neuroscience has a long, rich history in embracing unusual animals for research. Over the past several decades, there has been a technology-driven bottleneck in the species used for neuroscience research. However, an oncoming wave of technologies applicable to many animals hold promise for enabling researchers to address challenging scientific questions that cannot be solved using traditional laboratory animals. Here, we discuss how leveraging the convergent evolution of physiological or behavioral phenotypes can empower research mapping genotype to phenotype interactions. We present two case studies using electric fish and poison frogs and discuss how comparative work can teach us about evolutionary constraint and flexibility at various levels of biological organization. We also offer advice on the potential and pitfalls of establishing novel model systems in neuroscience research. Finally, we end with a discussion on the use of charismatic animals in neuroscience research and their utility in public outreach. Overall, we argue that convergent evolution frameworks can help identify generalizable principles of neuroscience.

Comparative analysis of integument transcriptomes identifies genes that participate in marking pattern formation in three allelic mutants of silkworm, Bombyx mori

The diversity markings and pigment patterns in insects are outcomes of adaptive evolution. The elucidation of the molecular mechanism underlying variations in pigment patterns may improve our understanding of the origin and evolution of these spectacular diverse phenotypes. Melanin, ommochrome, and pteridine are the three main types of insect pigments, and the genes that directly participate in pigment biosynthesis have been extensively studied. However, available information on gene interactions and the whole pigment regulatory network is limited. In this study, we performed integument transcriptome sequencing to analyze three larval marking allelic mutants, namely, multi lunar (L), L^C, and L^Ca, which have similar twin-spot markings on the dorsal side of multiple segments. Further analysis identified 336 differentially expressed genes (DEGs) between L and Dazao (wild type which exhibits normal markings), 68 DEGs between L^C/+ and +^LC/+^LC, and 188 DEGs between L^Ca/+ and +^LCa/+^LCa. Gene Ontology (GO) analysis indicated a significant DEG enrichment of the functional terms catalytic activity, binding, metabolic process, and cellular process. Furthermore, three mutants share six common enriched KEGG pathways. We finally identified eight common DEGs among three pairwise comparisons, including Krueppel-like factor, TATA-binding protein, protein patched, UDP-glycosyltransferase, an unknown secreted protein, and three cuticular proteins. Microarray-based gene expression analysis revealed that the eight genes are upregulated during molting, which coincides with marking formation, and are significantly differentially expressed between marking and non-marking regions. The results suggest that the eight common genes are involved in the construction of the multiple twin-spot marking patterns in the three mutants.

Modular cis-regulatory logic of yellow gene expression in silkmoth larvae

Colour patterns in butterflies and moths are crucial traits for adaptation. Previous investigations have highlighted genes responsible for pigmentation (ie yellow and ebony). However, the mechanisms by which these genes are regulated in lepidopteran insects remain poorly understood. To elucidate this, molecular studies involving dipterans have largely analysed the cis-regulatory regions of pigmentation genes and have revealed cis-regulatory modularity. Here, we used well-developed transgenic techniques in Bombyx mori and demonstrated that cis-regulatory modularity controls tissue-specific expression of the yellow gene. We first identified which body parts are regulated by the yellow gene via black pigmentation. We then isolated three discrete regulatory elements driving tissue-specific gene expression in three regions of B. mori larvae. Finally, we found that there is no apparent sequence conservation of cis-regulatory regions between B. mori and Drosophila melanogaster, and no expression driven by the regulatory regions of one species when introduced into the other species. Therefore, the trans-regulatory landscapes of the yellow gene differ significantly between the two taxa. The results of this study confirm that lepidopteran species use cis-regulatory modules to control gene expression related to pigmentation, and represent a powerful cadre of transgenic tools for studying evolutionary developmental mechanisms.

The genetic architecture of adaptation: convergence and pleiotropy in Heliconius wing pattern evolution

Unravelling the genetic basis of adaptive traits is a major challenge in evolutionary biology. Doing so informs our understanding of evolution towards an adaptive optimum, the distribution of locus effect sizes, and the influence of genetic architecture on the evolvability of a trait. In the Müllerian co-mimics Heliconius melpomene and Heliconius erato some Mendelian loci affecting mimicry shifts are well known. However, several phenotypes in H. melpomene remain to be mapped, and the quantitative genetics of colour pattern variation has rarely been analysed. Here we use quantitative trait loci (QTL) analyses of crosses between H. melpomene races from Peru and Suriname to map, for the first time, the control of the broken band phenotype to WntA and identify a ~100 kb region controlling this variation. Additionally, we map variation in basal forewing red-orange pigmentation to a locus centred around the gene ventral veins lacking (vvl). The locus also appears to affect medial band shape variation as it was previously known to do in H. erato. This adds to the list of homologous regions controlling convergent phenotypes between these two species. Finally we show that Heliconius wing-patterning genes are strikingly pleiotropic among wing pattern traits. Our results demonstrate how genetic architecture can shape, aid and constrain adaptive evolution.

Conservation and flexibility in the gene regulatory landscape of heliconiine butterfly wings

BACKGROUND

Many traits evolve by cis-regulatory modification, by which changes to noncoding sequences affect the binding affinity for available transcription factors and thus modify the expression profile of genes. Multiple examples of cis-regulatory evolution have been described at pattern switch genes responsible for butterfly wing pattern polymorphism, including in the diverse neotropical genus Heliconius, but the identities of the factors that can regulate these switch genes have not been identified.

RESULTS

We investigated the spatial transcriptomic landscape across the wings of three closely related butterfly species, two of which have a convergently evolved co-mimetic pattern and the other having a divergent pattern. We identified candidate factors for regulating the expression of wing patterning genes, including transcription factors with a conserved expression profile in all three species, and others, including both transcription factors and Wnt pathway genes, with markedly different profiles in each of the three species. We verified the conserved expression profile of the transcription factor homothorax by immunofluorescence and showed that its expression profile strongly correlates with that of the selector gene optix in butterflies with the Amazonian forewing pattern element 'dennis.'

CONCLUSION

Here we show that, in addition to factors with conserved expression profiles like homothorax, there are also a variety of transcription factors and signaling pathway components that appear to vary in their expression profiles between closely related butterfly species, highlighting the importance of genome-wide regulatory evolution between species.

The mimetic wing pattern of Papilio polytes butterflies is regulated by a doublesex-orchestrated gene network

The swallowtail butterfly Papilio polytes is sexually dimorphic and exhibits female-limited Batesian mimicry. This species also has two female forms, a non-mimetic form with male-like wing patterns, and a mimetic form resembling an unpalatable model, Pachliopta aristolochiae. The mimicry locus H constitutes a dimorphic Mendelian 'supergene', including a transcription factor gene doublesex (dsx). However, how the mimetic-type dsx (dsx-H) orchestrates the downstream gene network and causes the mimetic traits remains unclear. Here we performed RNA-seq-based gene screening and found that Wnt1 and Wnt6 are up-regulated by dsx-H during the early pupal stage and are involved in the red/white pigmentation and patterning of mimetic female wings. In contrast, a homeobox gene abdominal-A is repressed by dsx-H and involved in the non-mimetic colouration pattern. These findings suggest that dual regulation by dsx-H, induction of mimetic gene networks and repression of non-mimetic gene networks, is essential for the switch from non-mimetic to mimetic pattern in mimetic female wings.

A homeotic shift late in development drives mimetic color variation in a bumble bee

Natural phenotypic radiations, with their high diversity and convergence, are well-suited for informing how genomic changes translate to natural phenotypic variation. New genomic tools enable discovery in such traditionally nonmodel systems. Here, we characterize the genomic basis of color pattern variation in bumble bees (Hymenoptera, Apidae, Bombus), a group that has undergone extensive convergence of setal color patterns as a result of Müllerian mimicry. In western North America, multiple species converge on local mimicry patterns through parallel shifts of midabdominal segments from red to black. Using genome-wide association, we establish that a cis-regulatory locus between the abdominal fate-determining Hox genes, abd-A and Abd-B, controls the red-black color switch in a western species, Bombus melanopygus Gene expression analysis reveals distinct shifts in Abd-B aligned with the duration of setal pigmentation at the pupal-adult transition. This results in atypical anterior Abd-B expression, a late developmental homeotic shift. Changing expression of Hox genes can have widespread effects, given their important role across segmental phenotypes; however, the late timing reduces this pleiotropy, making Hox genes suitable targets. Analysis of this locus across mimics and relatives reveals that other species follow independent genetic routes to obtain the same phenotypes.

Development and evolution of color patterns in ladybird beetles: A case study in Harmonia axyridis

Many organisms show various geometric color patterns on their bodies, and the developmental, evolutionary, genetic, and ecological bases of these patterns have been intensely studied in various organisms. Ladybird beetles display highly diverse patterns of wing (elytral) color and are one of the most attractive model organisms for studying these characteristics. In this study, we reviewed the genetic history of elytral color patterns in the Asian multicolored ladybird beetle Harmonia axyridis from the classical genetic studies led by the pupils of Thomas Hunt Morgan and Theodosius Dobzhansky to recent genomic studies that revealed that a single GATA transcription factor gene, pannier, regulates the highly diverse elytral color patterns in this species. We also reviewed and discussed the developmental and evolutionary mechanisms driven by the pannier locus in H. axyridis. In the development sections, we focused on the following two topics: (a) how the red (carotenoid) and black (melanin) pigmentation of elytra is regulated by the pannier and pigmentation genes and (b) how the diverse color patterns are formed by integrating regulatory inputs from other genes involved in wing development. In the evolution section, we subsequently focused on the highly diversified DNA sequences within the first intron of pannier that are 56-76 kb long and that were generated through recurrent multiple inversions. Furthermore, we discussed how these recurrent inversions have driven the diversification of color patterns throughout evolution.

Insights into eyespot color-pattern formation mechanisms from color gradients, boundary scales, and rudimentary eyespots in butterfly wings

Butterfly eyespot color patterns are traditionally explained by the gradient model, where positional information is stably provided by a morphogen gradient from a single organizer and its output is a set of non-graded (or graded) colors based on pre-determined threshold levels. An alternative model is the induction model, in which the outer black ring and the inner black core disk of an eyespot are specified by graded signals from the primary and secondary organizers that also involve lateral induction. To examine the feasibility of these models, we analyzed eyespot color gradients, boundary scales, and rudimentary eyespots in various nymphalid butterflies. Most parts of eyespots showed color gradients with gradual or fluctuating changes with sharp boundaries in many species, but some species had eyespots that were composed of a constant color within a given part. Thus, a plausible model should be flexible enough to incorporate this diversity. Some boundary scales appeared to have two kinds of pigments, and others had "misplaced" colors, suggesting an overlapping of two signals and a difficulty in assuming sharp threshold boundaries. Rudimentary eyespots of three Junonia species revealed that the outer black ring is likely determined first and the inner yellow or red ring is laterally induced. This outside-to-inside determination together with the lateral induction may favor the induction model, in which dynamic signal interactions play a major role. The implications of these results for the ploidy hypothesis and color-pattern rules are discussed.

How a growing organismal perspective is adding new depth to integrative studies of morphological evolution

Over the past half century, the field of Evolutionary Developmental Biology, or Evo-devo, has integrated diverse fields of biology into a more synthetic understanding of morphological diversity. This has resulted in numerous insights into how development can evolve and reciprocally influence morphological evolution, as well as generated several novel theoretical areas. Although comparative by default, there remains a great gap in our understanding of adaptive morphological diversification and how developmental mechanisms influence the shape and pattern of phenotypic variation. Herein we highlight areas of research that are in the process of filling this void, and areas, if investigated more fully, that will add new insights into the diversification of morphology. At the centre of our discussion is an explicit awareness of organismal biology. Here we discuss an organismal framework that is supported by three distinct pillars. First, there is a need for Evo-devo to adopt a high-resolution phylogenetic approach in the study of morphological variation and its developmental underpinnings. Secondly, we propose that to understand the dynamic nature of morphological evolution, investigators need to give more explicit attention to the processes that generate evolutionarily relevant variation at the population level. Finally, we emphasize the need to address more thoroughly the processes that structure variation at micro- and macroevolutionary scales including modularity, morphological integration, constraint, and plasticity. We illustrate the power of these three pillars using numerous examples from both invertebrates and vertebrates to emphasize that many of these approaches are already present within the field, but have yet to be formally integrated into many research programs. We feel that the most exciting new insights will come where the traditional experimental approaches to Evo-devo are integrated more thoroughly with the principles of this organismal framework.

Aristaless Controls Butterfly Wing Color Variation Used in Mimicry and Mate Choice

Neotropical Heliconius butterflies display a diversity of warningly colored wing patterns, which serve roles in both Müllerian mimicry and mate choice behavior. Wing pattern diversity in Heliconius is controlled by a small number of unlinked, Mendelian "switch" loci [1]. One of these, termed the K locus, switches between yellow and white color patterns, important mimicry signals as well as mating cues [2-4]. Furthermore, mate preference behavior is tightly linked to this locus [4]. K controls the distribution of white versus yellow scales on the wing, with a dominant white allele and a recessive yellow allele. Here, we combine fine-scale genetic mapping, genome-wide association studies, gene expression analyses, population and comparative genomics, and genome editing with CRISPR/Cas9 to characterize the molecular basis of the K locus in Heliconius and to infer its evolutionary history. We show that white versus yellow color variation in Heliconius cydno is due to alternate haplotypes at a putative cis-regulatory element (CRE) downstream of a tandem duplication of the homeodomain transcription factor aristaless. Aristaless1 (al1) and aristaless2 (al2) are differentially regulated between white and yellow wings throughout development with elevated expression of al1 in developing white wings, suggesting a role in repressing pigmentation. Consistent with this, knockout of al1 causes white wings to become yellow. The evolution of wing color in this group has been marked by retention of the ancestral yellow color in many lineages, a single origin of white coloration in H. cydno, and subsequent introgression of white color from H. cydno into H. melpomene.

From pattern to process: studies at the interface of gene regulatory networks, morphogenesis, and evolution

The development of anatomical structures is complex, beginning with patterning of gene expression by multiple gene regulatory networks (GRNs). These networks ultimately regulate the activity of effector molecules, which in turn alter cellular behavior during development. Together these processes biomechanically produce the three-dimensional shape that the anatomical structure adopts over time. However, the interfaces between these processes are often overlooked and also include counter-intuitive feedback mechanisms. In this review, we examine each step in this extraordinarily complex process and explore how evolutionary developmental biology model systems, such as butterfly scales, vertebrate teeth, and the Drosophila dorsal appendage offer a complementary approach to expose the multifactorial integration of genetics and morphogenesis from an alternative perspective.

Genomics of coloration in natural animal populations

Animal coloration has traditionally been the target of genetic and evolutionary studies. However, until very recently, the study of the genetic basis of animal coloration has been mainly restricted to model species, whereas research on non-model species has been either neglected or mainly based on candidate approaches, and thereby limited by the knowledge obtained in model species. Recent high-throughput sequencing technologies allow us to overcome previous limitations, and open new avenues to study the genetic basis of animal coloration in a broader number of species and colour traits, and to address the general relevance of different genetic structures and their implications for the evolution of colour. In this review, we highlight aspects where genome-wide studies could be of major utility to fill in the gaps in our understanding of the biology and evolution of animal coloration. The new genomic approaches have been promptly adopted to study animal coloration although substantial work is still needed to consider a larger range of species and colour traits, such as those exhibiting continuous variation or based on reflective structures. We argue that a robust advancement in the study of animal coloration will also require large efforts to validate the functional role of the genes and variants discovered using genome-wide tools.This article is part of the themed issue 'Animal coloration: production, perception, function and application'.

Learning of salient prey traits explains Batesian mimicry evolution

Batesian mimicry evolution involves an initial major mutation that produces a rough resemblance to the model, followed by smaller improving changes. To examine the learning psychology of this process, we applied established ideas about mimicry in Papilio polyxenes asterius of the model Battus philenor. We performed experiments with wild birds as predators and butterfly wings as semiartificial prey. Wings of hybrids of P. p. asterius and Papilio machaon were used to approximate the first mutant, with melanism as the hypothesized first mimetic trait. Based on previous results about learning psychology and imperfect mimicry, we predicted that: melanism should have high salience (i.e., being noticeable and prominent), meaning that predators readily discriminate a melanistic mutant from appearances similar to P. machaon; the difference between the first mutant and the model should have intermediate salience to allow further improvement of mimicry; and the final difference in appearance between P. p. asterius and B. philenor should have very low salience, causing improvement to level off. Our results supported both the traditional hypothesis and all our predictions about relative salience. We conclude that there is good agreement between long-held ideas about how Batesian mimicry evolves and recent insights from learning psychology about the role of salience in mimicry evolution.

Waiting in the wings: what can we learn about gene co-option from the diversification of butterfly wing patterns?

A major challenge is to understand how conserved gene regulatory networks control the wonderful diversity of form that we see among animals and plants. Butterfly wing patterns are an excellent example of this diversity. Butterfly wings form as imaginal discs in the caterpillar and are constructed by a gene regulatory network, much of which is conserved across the holometabolous insects. Recent work in Heliconius butterflies takes advantage of genomic approaches and offers insights into how the diversification of wing patterns is overlaid onto this conserved network. WntA is a patterning morphogen that alters spatial information in the wing. Optix is a transcription factor that acts later in development to paint specific wing regions red. Both of these loci fit the paradigm of conserved protein-coding loci with diverse regulatory elements and developmental roles that have taken on novel derived functions in patterning wings. These discoveries offer insights into the 'Nymphalid Ground Plan', which offers a unifying hypothesis for pattern formation across nymphalid butterflies. These loci also represent 'hotspots' for morphological change that have been targeted repeatedly during evolution. Both convergent and divergent evolution of a great diversity of patterns is controlled by complex alleles at just a few genes. We suggest that evolutionary change has become focused on one or a few genetic loci for two reasons. First, pre-existing complex cis-regulatory loci that already interact with potentially relevant transcription factors are more likely to acquire novel functions in wing patterning. Second, the shape of wing regulatory networks may constrain evolutionary change to one or a few loci. Overall, genomic approaches that have identified wing patterning loci in these butterflies offer broad insight into how gene regulatory networks evolve to produce diversity.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Macroevolutionary shifts of WntA function potentiate butterfly wing-pattern diversity

Butterfly wing patterns provide a rich comparative framework to study how morphological complexity develops and evolves. Here we used CRISPR/Cas9 somatic mutagenesis to test a patterning role for WntA, a signaling ligand gene previously identified as a hotspot of shape-tuning alleles involved in wing mimicry. We show that WntA loss-of-function causes multiple modifications of pattern elements in seven nymphalid butterfly species. In three butterflies with a conserved wing-pattern arrangement, WntA is necessary for the induction of stripe-like patterns known as symmetry systems and acquired a novel eyespot activator role specific to Vanessa forewings. In two Heliconius species, WntA specifies the boundaries between melanic fields and the light-color patterns that they contour. In the passionvine butterfly Agraulis, WntA removal shows opposite effects on adjacent pattern elements, revealing a dual role across the wing field. Finally, WntA acquired a divergent role in the patterning of interveinous patterns in the monarch, a basal nymphalid butterfly that lacks stripe-like symmetry systems. These results identify WntA as an instructive signal for the prepatterning of a biological system of exuberant diversity and illustrate how shifts in the deployment and effects of a single developmental gene underlie morphological change.

Mimicry in butterflies: co-option and a bag of magnificent developmental genetic tricks

Butterfly wing patterns are key adaptations that are controlled by remarkable developmental and genetic mechanisms that facilitate rapid evolutionary change. With swift advancements in the fields of genomics and genetic manipulations, identifying the regulators of wing development and mimetic wing patterns has become feasible even in nonmodel organisms such as butterflies. Recent mapping and gene expression studies have identified single switch loci of major effects such as transcription factors and supergenes as the main drivers of adaptive evolution of mimetic and polymorphic butterfly wing patterns. We highlight several of these examples, with emphasis on doublesex, optix, WntA and other dynamic, yet essential, master regulators that control critical color variation and sex-specific traits. Co-option emerges as a predominant theme, where typically embryonic and other early-stage developmental genes and networks have been rewired to regulate polymorphic and sex-limited mimetic wing patterns in iconic butterfly adaptations. Drawing comparisons from our knowledge of wing development in Drosophila, we illustrate the functional space of genes that have been recruited to regulate butterfly wing patterns. We also propose a developmental pathway that potentially results in dorsoventral mismatch in butterfly wing patterns. Such dorsoventrally mismatched color patterns modulate signal components of butterfly wings that are used in intra- and inter-specific communication. Recent advances-fuelled by RNAi-mediated knockdowns and CRISPR/Cas9-based genomic edits-in the developmental genetics of butterfly wing patterns, and the underlying biological diversity and complexity of wing coloration, are pushing butterflies as an emerging model system in ecological genetics and evolutionary developmental biology. WIREs Dev Biol 2018, 7:e291. doi: 10.1002/wdev.291 This article is categorized under: Gene Expression and Transcriptional Hierarchies > Regulatory Mechanisms Comparative Development and Evolution > Regulation of Organ Diversity Comparative Development and Evolution > Evolutionary Novelties.

Investigating the genomic basis of discrete phenotypes using a Pool-Seq-only approach: New insights into the genetics underlying colour variation in diverse taxa

While large-scale genomic approaches are increasingly revealing the genetic basis of polymorphic phenotypes such as colour morphs, such approaches are almost exclusively conducted in species with high-quality genomes and annotations. Here, we use Pool-Seq data for both genome assembly and SNP frequency estimation, followed by scanning for F_ST outliers to identify divergent genomic regions. Using paired-end, short-read sequencing data from two groups of individuals expressing divergent phenotypes, we generate a de novo rough-draft genome, identify SNPs and calculate genomewide F_ST differences between phenotypic groups. As genomes generated by Pool-Seq data are highly fragmented, we also present an approach for super-scaffolding contigs using existing protein-coding data sets. Using this approach, we reanalysed genomic data from two recent studies of birds and butterflies investigating colour pattern variation and replicated their core findings, demonstrating the accuracy and power of a Pool-Seq-only approach. Additionally, we discovered new regions of high divergence and new annotations that together suggest novel parallels between birds and butterflies in the origins of their colour pattern variation.

Deep-Time Convergence in Rove Beetle Symbionts of Army Ants

Recent adaptive radiations provide striking examples of convergence [1-4], but the predictability of evolution over much deeper timescales is controversial, with a scarcity of ancient clades exhibiting repetitive patterns of phenotypic evolution [5, 6]. Army ants are ecologically dominant arthropod predators of the world's tropics, with large nomadic colonies housing diverse communities of socially parasitic myrmecophiles [7]. Remarkable among these are many species of rove beetle (Staphylinidae) that exhibit ant-mimicking "myrmecoid" body forms and are behaviorally accepted into their aggressive hosts' societies: emigrating with colonies and inhabiting temporary nest bivouacs, grooming and feeding with workers, but also consuming the brood [8-11]. Here, we demonstrate that myrmecoid rove beetles are strongly polyphyletic, with this adaptive morphological and behavioral syndrome having evolved at least 12 times during the evolution of a single staphylinid subfamily, Aleocharinae. Each independent myrmecoid clade is restricted to one zoogeographic region and highly host specific on a single army ant genus. Dating estimates reveal that myrmecoid clades are separated by substantial phylogenetic distances-as much as 105 million years. All such groups arose in parallel during the Cenozoic, when army ants diversified into modern genera [12] and rose to ecological dominance [13, 14]. This work uncovers a rare example of an ancient system of complex morphological and behavioral convergence, with replicate beetle lineages following a predictable phenotypic trajectory during their parasitic adaptation to host colonies.

Complex modular architecture around a simple toolkit of wing pattern genes

Identifying the genomic changes that control morphological variation and understanding how they generate diversity is a major goal of evolutionary biology. In Heliconius butterflies, a small number of genes control the development of diverse wing color patterns. Here, we used full genome sequencing of individuals across the Heliconius erato radiation and closely related species to characterize genomic variation associated with wing pattern diversity. We show that variation around color pattern genes is highly modular, with narrow genomic intervals associated with specific differences in color and pattern. This modular architecture explains the diversity of color patterns and provides a flexible mechanism for rapid morphological diversification.

Genes controlling mimetic colour pattern variation in butterflies

Butterfly wing patterns are made up of arrays of coloured scales. There are two genera in which within-species variation in wing patterning is common and has been investigated at the molecular level, Heliconius and Papilio. Both of these species have mimetic relationships with other butterfly species that increase their protection from predators. Heliconius have a 'tool-kit' of five genetic loci that control colour pattern, three of which have been identified at the gene level, and which have been repeatedly used to modify colour pattern by different species in the genus. By contrast, the three Papilio species that have been investigated each have different genetic mechanisms controlling their polymorphic wing patterns.

Advances in finding Alba: the locus affecting life history and color polymorphism in a Colias butterfly

Although alternative life-history strategies exist within many populations, very little is known about their genetic basis and mechanistic insight into these traits could greatly advance the understanding of eco-evolutionary dynamics. Many species of butterfly within the genus Colias exhibit a sex-limited wing colour polymorphism, called Alba, which is correlated with an alternative life-history strategy. Here, we have taken the first steps in localizing the region carrying Alba in Colias croceus, a species with no genomic resources, by generating whole genome sequence of a single Alba mother and two sequencing pools, one for her Alba and another for her orange, offspring. These data were used in a bulk-segregant analysis wherein SNPs fulfilling the Mendelian inheritance expectations of Alba were identified. Then, using the conserved synteny in Lepidoptera, the Alba locus was assigned to chromosome 15 in Bombyx mori. We then identified candidate regions within the chromosome by investigating the distribution of Alba SNPs along the chromosome and the difference in nucleotide diversity in exons between the two pools. A region spanning ~ 5.7 Mbp at the 5' end of the chromosome was identified as likely to contain the Alba locus. These insights set the stage for more detailed genomic scans and mapping of the Alba phenotype, and demonstrate an efficient use of genomic resources in a novel species.

The phenotypic plasticity of developmental modules

Background

Organisms develop and evolve in a modular fashion, but how individual modules interact with the environment remains poorly understood. Phenotypically plastic traits are often under selection, and studies are needed to address how traits respond to the environment in a modular fashion. In this study, tissue-specific plasticity of melanic spots was examined in the large milkweed bug, Oncopeltus fasciatus.

Results

Although the size of the abdominal melanic bands varied according to rearing temperatures, wing melanic bands were more robust. To explore the regulation of abdominal pigmentation plasticity, candidate genes involved in abdominal melanic spot patterning and biosynthesis of melanin were analyzed. While the knockdown of dopa decarboxylase (Ddc) led to lighter pigmentation in both the wings and the abdomen, the shape of the melanic elements remained unaffected. Although the knockdown of Abdominal-B (Abd-B) partially phenocopied the low-temperature phenotype, the abdominal bands were still sensitive to temperature shifts. These observations suggest that regulators downstream of Abd-B but upstream of DDC are responsible for the temperature response of the abdomen. Ablation of wings led to the regeneration of a smaller wing with reduced melanic bands that were shifted proximally. In addition, the knockdown of the Wnt signaling nuclear effector genes, armadillo1 and armadillo 2, altered both the melanic bands and the wing shape. Thus, the pleiotropic effects of Wnt signaling may constrain the amount of plasticity in wing melanic bands.

Conclusions

We propose that when traits are regulated by distinct pre-patterning mechanisms, they can respond to the environment in a modular fashion, whereas when the environment impacts developmental regulators that are shared between different modules, phenotypic plasticity can manifest as a developmentally integrated system.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-016-0053-7) contains supplementary material, which is available to authorized users.

Genome-wide introgression among distantly related Heliconius butterfly species

Background

Although hybridization is thought to be relatively rare in animals, the raw genetic material introduced via introgression may play an important role in fueling adaptation and adaptive radiation. The butterfly genus Heliconius is an excellent system to study hybridization and introgression but most studies have focused on closely related species such as H. cydno and H. melpomene. Here we characterize genome-wide patterns of introgression between H. besckei, the only species with a red and yellow banded ‘postman’ wing pattern in the tiger-striped silvaniform clade, and co-mimetic H. melpomene nanna.

Results

We find a pronounced signature of putative introgression from H. melpomene into H. besckei in the genomic region upstream of the gene optix, known to control red wing patterning, suggesting adaptive introgression of wing pattern mimicry between these two distantly related species. At least 39 additional genomic regions show signals of introgression as strong or stronger than this mimicry locus. Gene flow has been on-going, with evidence of gene exchange at multiple time points, and bidirectional, moving from the melpomene to the silvaniform clade and vice versa. The history of gene exchange has also been complex, with contributions from multiple silvaniform species in addition to H. besckei. We also detect a signature of ancient introgression of the entire Z chromosome between the silvaniform and melpomene/cydno clades.

Conclusions

Our study provides a genome-wide portrait of introgression between distantly related butterfly species. We further propose a comprehensive and efficient workflow for gene flow identification in genomic data sets.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-0889-0) contains supplementary material, which is available to authorized users.

The functional basis of wing patterning in Heliconius butterflies: the molecules behind mimicry

Wing-pattern mimicry in butterflies has provided an important example of adaptation since Charles Darwin and Alfred Russell Wallace proposed evolution by natural selection >150 years ago. The neotropical butterfly genus Heliconius played a central role in the development of mimicry theory and has since been studied extensively in the context of ecology and population biology, behavior, and mimicry genetics. Heliconius species are notable for their diverse color patterns, and previous crossing experiments revealed that much of this variation is controlled by a small number of large-effect, Mendelian switch loci. Recent comparative analyses have shown that the same switch loci control wing-pattern diversity throughout the genus, and a number of these have now been positionally cloned. Using a combination of comparative genetic mapping, association tests, and gene expression analyses, variation in red wing patterning throughout Heliconius has been traced back to the action of the transcription factor optix. Similarly, the signaling ligand WntA has been shown to control variation in melanin patterning across Heliconius and other butterflies. Our understanding of the molecular basis of Heliconius mimicry is now providing important insights into a variety of additional evolutionary phenomena, including the origin of supergenes, the interplay between constraint and evolvability, the genetic basis of convergence, the potential for introgression to facilitate adaptation, the mechanisms of hybrid speciation in animals, and the process of ecological speciation.