The gene cortex controls mimicry and crypsis in butterflies and moths

The draft genome sequence of the grove snail Cepaea nemoralis

Studies on the shell color and banding polymorphism of the grove snail Cepaea nemoralis and the sister taxon Cepaea hortensis have provided compelling evidence for the fundamental role of natural selection in promoting and maintaining intraspecific variation. More recently, Cepaea has been the focus of citizen science projects on shell color evolution in relation to climate change and urbanization. C. nemoralis is particularly useful for studies on the genetics of shell polymorphism and the evolution of "supergenes," as well as evo-devo studies of shell biomineralization, because it is relatively easily maintained in captivity. However, an absence of genomic resources for C. nemoralis has generally hindered detailed genetic and molecular investigations. We therefore generated ∼23× coverage long-read data for the ∼3.5 Gb genome, and produced a draft assembly composed of 28,537 contigs with the N50 length of 333 kb. Genome completeness, estimated by BUSCO using the metazoa dataset, was 91%. Repetitive regions cover over 77% of the genome. A total of 43,519 protein-coding genes were predicted in the assembled genome, and 97.3% of these were functionally annotated from either sequence homology or protein signature searches. This first assembled and annotated genome sequence for a helicoid snail, a large group that includes edible species, agricultural pests, and parasite hosts, will be a core resource for identifying the loci that determine the shell polymorphism, as well as in a wide range of analyses in evolutionary and developmental biology, and snail biology in general.

Haplotype tagging reveals parallel formation of hybrid races in two butterfly species

Genetic variation segregates as linked sets of variants or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. Yet, genomic data often omit haplotype information due to constraints in sequencing technologies. Here, we present "haplotagging," a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabase-size haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high- and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species, the geographic clines for the major wing-pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the center of the hybrid zone. We propose that shared warning signaling (Müllerian mimicry) may couple the cline shifts seen in both species and facilitate the parallel coemergence of a novel hybrid morph in both comimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.

Deploying Big Data to Crack the Genotype to Phenotype Code

Mechanistically connecting genotypes to phenotypes is a longstanding and central mission of biology. Deciphering these connections will unite questions and datasets across all scales from molecules to ecosystems. Although high-throughput sequencing has provided a rich platform on which to launch this effort, tools for deciphering mechanisms further along the genome to phenome pipeline remain limited. Machine learning approaches and other emerging computational tools hold the promise of augmenting human efforts to overcome these obstacles. This vision paper is the result of a Reintegrating Biology Workshop, bringing together the perspectives of integrative and comparative biologists to survey challenges and opportunities in cracking the genotype to phenotype code and thereby generating predictive frameworks across biological scales. Key recommendations include promoting the development of minimum "best practices" for the experimental design and collection of data; fostering sustained and long-term data repositories; promoting programs that recruit, train, and retain a diversity of talent; and providing funding to effectively support these highly cross-disciplinary efforts. We follow this discussion by highlighting a few specific transformative research opportunities that will be advanced by these efforts.

Correlational selection in the age of genomics

Ecologists and evolutionary biologists are well aware that natural and sexual selection do not operate on traits in isolation, but instead act on combinations of traits. This long-recognized and pervasive phenomenon is known as multivariate selection, or-in the particular case where it favours correlations between interacting traits-correlational selection. Despite broad acknowledgement of correlational selection, the relevant theory has often been overlooked in genomic research. Here, we discuss theory and empirical findings from ecological, quantitative genetic and genomic research, linking key insights from different fields. Correlational selection can operate on both discrete trait combinations and quantitative characters, with profound implications for genomic architecture, linkage, pleiotropy, evolvability, modularity, phenotypic integration and phenotypic plasticity. We synthesize current knowledge and discuss promising research approaches that will enable us to understand how correlational selection shapes genomic architecture, thereby linking quantitative genetic approaches with emerging genomic methods. We suggest that research on correlational selection has great potential to integrate multiple fields in evolutionary biology, including developmental and functional biology, ecology, quantitative genetics, phenotypic polymorphisms, hybrid zones and speciation processes.

Developmental Transcriptomics Reveals a Gene Network Driving Mimetic Color Variation in a Bumble Bee

A major goal of evolutionary genetics and evo-devo is to understand how changes in genotype manifest as changes in phenotype. Bumble bees display remarkable color pattern diversity while converging onto numerous regional Müllerian mimicry patterns, thus enabling exploration of the genetic mechanisms underlying convergent phenotypic evolution. In western North America, multiple bumble bee species converge onto local mimicry patterns through parallel shifts of midabdominal segments from red to black. It was previously demonstrated that a Hox gene, Abd-B, is the key regulator of the phenotypic switch in one of these species, Bombus melanopygus, however, the mechanism by which Abd-B regulates color differentiation remains unclear. Using tissue/stage-specific transcriptomic analysis followed by qRT–PCR validation, this study reveals a suite of genes potentially involved downstream of Abd-B during color pattern differentiation. The data support differential genes expression of not only the first switch gene Abd-B, but also an intermediate developmental gene nubbin, and a whole suite of downstream melanin and redox genes that together reinforce the observed eumelanin (black)-pheomelanin (red) ratios. These include potential genes involved in the production of insect pheomelanins, a pigment until recently not thought to occur in insects and thus lacking known regulatory enzymes. The results enhance understanding of pigmentation gene networks involved in bumble bee color pattern development and diversification, while providing insights into how upstream regulators such as Hox genes interact with downstream morphogenic players to facilitate this adaptive phenotypic radiation.

A combined RAD-Seq and WGS approach reveals the genomic basis of yellow color variation in bumble bee Bombus terrestris

Bumble bees exhibit exceptional diversity in their segmental body coloration largely as a result of mimicry. In this study we sought to discover genes involved in this variation through studying a lab-generated mutant in bumble bee Bombus terrestris, in which the typical black coloration of the pleuron, scutellum, and first metasomal tergite is replaced by yellow, a color variant also found in sister lineages to B. terrestris. Utilizing a combination of RAD-Seq and whole-genome re-sequencing, we localized the color-generating variant to a single SNP in the protein-coding sequence of transcription factor cut. This mutation generates an amino acid change that modifies the conformation of a coiled-coil structure outside DNA-binding domains. We found that all sequenced Hymenoptera, including sister lineages, possess the non-mutant allele, indicating different mechanisms are involved in the same color transition in nature. Cut is important for multiple facets of development, yet this mutation generated no noticeable external phenotypic effects outside of setal characteristics. Reproductive capacity was reduced, however, as queens were less likely to mate and produce female offspring, exhibiting behavior similar to that of workers. Our research implicates a novel developmental player in pigmentation, and potentially caste, thus contributing to a better understanding of the evolution of diversity in both of these processes.

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 Genomic Architecture and Evolutionary Fates of Supergenes

Supergenes are genomic regions containing sets of tightly linked loci that control multi-trait phenotypic polymorphisms under balancing selection. Recent advances in genomics have uncovered significant variation in both the genomic architecture as well as the mode of origin of supergenes across diverse organismal systems. Although the role of genomic architecture for the origin of supergenes has been much discussed, differences in the genomic architecture also subsequently affect the evolutionary trajectory of supergenes and the rate of degeneration of supergene haplotypes. In this review, we synthesize recent genomic work and historical models of supergene evolution, highlighting how the genomic architecture of supergenes affects their evolutionary fate. We discuss how recent findings on classic supergenes involved in governing ant colony social form, mimicry in butterflies, and heterostyly in flowering plants relate to theoretical expectations. Furthermore, we use forward simulations to demonstrate that differences in genomic architecture affect the degeneration of supergenes. Finally, we discuss implications of the evolution of supergene haplotypes for the long-term fate of balanced polymorphisms governed by supergenes.

Mutation load at a mimicry supergene sheds new light on the evolution of inversion polymorphisms

Chromosomal inversions are ubiquitous in genomes and often coordinate complex phenotypes, such as the covariation of behavior and morphology in many birds, fishes, insects or mammals^1-11. However, why and how inversions become associated with polymorphic traits remains obscure. Here we show that despite a strong selective advantage when they form, inversions accumulate recessive deleterious mutations that generate frequency-dependent selection and promote their maintenance at intermediate frequency. Combining genomics and in vivo fitness analyses in a model butterfly for wing-pattern polymorphism, Heliconius numata, we reveal that three ecologically advantageous inversions have built up a heavy mutational load from the sequential accumulation of deleterious mutations and transposable elements. Inversions associate with sharply reduced viability when homozygous, which prevents them from replacing ancestral chromosome arrangements. Our results suggest that other complex polymorphisms, rather than representing adaptations to competing ecological optima, could evolve because chromosomal rearrangements are intrinsically prone to carrying recessive harmful mutations.

The Fractal Geometry of the Nymphalid Groundplan: Self-Similar Configuration of Color Pattern Symmetry Systems in Butterfly Wings

Simple Summary

Highly diverse color patterns of butterfly wings can be explained as modifications of an archetypical color pattern of nymphalid butterflies called the nymphalid groundplan. The nymphalid groundplan contains three major symmetry systems and a discal symmetry system, but their relationships have been elusive. Here, the morphological and spatial relationships among these symmetry systems were studied based on cross-species color-pattern comparisons of the hindwings in nymphalid butterflies. It was shown that all symmetry systems can be expressed as various structures, suggesting the equivalence (homology) of these systems in developmental potential. In some cases, the discal symmetry system is circularly surrounded by the central symmetry system, which may then be surrounded by the border and basal symmetry systems, indicating a unified supersymmetry system covering the entire wing. These results suggest that butterfly color patterns are hierarchically constructed; one system is nested within another system, which is a self-similar relationship that achieves the fractal geometry. This self-similarity is likely mediated by the serial induction of organizers during development, and a possible mechanism is proposed for symmetry breaking of the system morphology, which contributes to the diversity of butterfly wing color patterns.

Abstract

The nymphalid groundplan is an archetypical color pattern of nymphalid butterflies involving three major symmetry systems and a discal symmetry system, which share the basic morphogenesis unit. Here, the morphological and spatial relationships among these symmetry systems were studied based on cross-species comparisons of nymphalid hindwings. Based on findings in Neope and Symbrenthia, all three major symmetry systems can be expressed as bands, spots, or eyespot-like structures, suggesting equivalence (homology) of these systems in developmental potential. The discal symmetry system can also be expressed as various structures. The discal symmetry system is circularly surrounded by the central symmetry system, which may then be surrounded by the border and basal symmetry systems, based mainly on findings in Agrias, indicating a unified supersymmetry system covering the entire wing. The border symmetry system can occupy the central part of the wing when the central symmetry system is compromised, as seen in Callicore. These results suggest that butterfly color patterns are hierarchically constructed in a self-similar fashion, as the fractal geometry of the nymphalid groundplan. This self-similarity is likely mediated by the serial induction of organizers, and symmetry breaking of the system morphology may be generated by the collision of opposing signals during development.

Determining the probability of hemiplasy in the presence of incomplete lineage sorting and introgression

The incongruence of character states with phylogenetic relationships is often interpreted as evidence of convergent evolution. However, trait evolution along discordant gene trees can also generate these incongruences - a phenomenon known as hemiplasy. Classic comparative methods do not account for discordance, resulting in incorrect inferences about the number, timing, and direction of trait transitions. Biological sources of discordance include incomplete lineage sorting (ILS) and introgression, but only ILS has received theoretical consideration in the context of hemiplasy. Here, we present a model that shows introgression makes hemiplasy more likely, such that methods that account for ILS alone will be conservative. We also present a method and software (HeIST) for making statistical inferences about the probability of hemiplasy and homoplasy in large datasets that contain both ILS and introgression. We apply our methods to two empirical datasets, finding that hemiplasy is likely to contribute to the observed trait incongruences in both.

The roles of wing color pattern and geography in the evolution of Neotropical Preponini butterflies

Diversification rates and evolutionary trajectories are known to be influenced by phenotypic traits and the geographic history of the landscapes that organisms inhabit. One of the most conspicuous traits in butterflies is their wing color pattern, which has been shown to be important in speciation. The evolution of many taxa in the Neotropics has also been influenced by major geological events. Using a dated, species-level molecular phylogenetic hypothesis for Preponini, a colorful Neotropical butterfly tribe, we evaluated whether diversification rates were constant or varied through time, and how they were influenced by color pattern evolution and biogeographical events. We found that Preponini originated approximately 28 million years ago and that diversification has increased through time consistent with major periods of Andean uplift. Even though some clades show evolutionarily rapid transitions in coloration, contrary to our expectations, these shifts were not correlated with shifts in diversification. Involvement in mimicry with other butterfly groups might explain the rapid changes in dorsal color patterns in this tribe, but such changes have not increased species diversification in this group. However, we found evidence for an influence of major Miocene and Pliocene geological events on the tribe's evolution. Preponini apparently originated within South America, and range evolution has since been dynamic, congruent with Andean geologic activity, closure of the Panama Isthmus, and Miocene climate variability.

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.

Visual mate preference evolution during butterfly speciation is linked to neural processing genes

Many animal species remain separate not because their individuals fail to produce viable hybrids but because they "choose" not to mate. However, we still know very little of the genetic mechanisms underlying changes in these mate preference behaviours. Heliconius butterflies display bright warning patterns, which they also use to recognize conspecifics. Here, we couple QTL for divergence in visual preference behaviours with population genomic and gene expression analyses of neural tissue (central brain, optic lobes and ommatidia) across development in two sympatric Heliconius species. Within a region containing 200 genes, we identify five genes that are strongly associated with divergent visual preferences. Three of these have previously been implicated in key components of neural signalling (specifically an ionotropic glutamate receptor and two regucalcins), and overall our candidates suggest shifts in behaviour involve changes in visual integration or processing. This would allow preference evolution without altering perception of the wider environment.

Ancient incomplete lineage sorting of Hyles and Rhodafra (Lepidoptera: Sphingidae)

The hawkmoth genus Rhodafra comprises two African species with unclear relationships, as their wing patterns are markedly different, with one species closely resembling species of a related genus, Hyles. The present paper aims to investigate the monophyly and phylogenetic position of Rhodafra in relation to Hyles and other genera of the subtribe Choerocampina (Sphingidae: Macroglossinae: Macroglossini) using mitochondrial and nuclear sequence data from more species and individuals than have hitherto been studied. As no fresh tissue of Rhodafra was available, ancient-DNA methodology was applied. All data corroborate the genus as monophyletic and that a similar wing pattern is not a good indicator of close phylogenetic relationship in this group of moths. Phylogenetic trees based on mitochondrial data agree in placing Rhodafra within Hyles. In contrast, analysis of nuclear EF1alpha sequences produces a topology in which Rhodafra is placed as the sister clade to Hyles. Although multispecies coalescent analyses suggest a polytomy between Rhodafra, Hyles lineata and the remaining Hyles, total evidence analyses corroborate Rhodafra as sister to Hyles. This relationship is interpreted as the favoured topology. For a more robust result, the question should be re-examined using genomic approaches.

Tissue-specific developmental regulation and isoform usage underlie the role of doublesex in sex differentiation and mimicry in Papilio swallowtails

Adaptive phenotypes often arise by rewiring existing developmental networks. Co-option of transcription factors in novel contexts has facilitated the evolution of ecologically important adaptations. doublesex (dsx) governs fundamental sex differentiation during embryonic stages and has been co-opted to regulate diverse secondary sexual dimorphisms during pupal development of holometabolous insects. In Papilio polytes, dsx regulates female-limited mimetic polymorphism, resulting in mimetic and non-mimetic forms. To understand how a critical gene such as dsx regulates novel wing patterns while maintaining its basic function in sex differentiation, we traced its expression through metamorphosis in P. polytes using developmental transcriptome data. We found three key dsx expression peaks: (i) eggs in pre- and post-ovisposition stages; (ii) developing wing discs and body in final larval instar; and (iii) 3-day pupae. We identified potential dsx targets using co-expression and differential expression analysis, and found distinct, non-overlapping sets of genes-containing putative dsx-binding sites-in developing wings versus abdominal tissue and in mimetic versus non-mimetic individuals. This suggests that dsx regulates distinct downstream targets in different tissues and wing colour morphs and has perhaps acquired new, previously unknown targets, for regulating mimetic polymorphism. Additionally, we observed that the three female isoforms of dsx were differentially expressed across stages (from eggs to adults) and tissues and differed in their protein structure. This may promote differential protein-protein interactions for each isoform and facilitate sub-functionalization of dsx activity across its isoforms. Our findings suggest that dsx employs tissue-specific downstream effectors and partitions its functions across multiple isoforms to regulate primary and secondary sexual dimorphism through insect development.

Many functionally connected loci foster adaptive diversification along a neotropical hybrid zone

Characterizing the genetic complexity of adaptation and trait evolution is a major emphasis of evolutionary biology and genetics. Incongruent findings from genetic studies have resulted in conceptual models ranging from a few large-effect loci to massively polygenic architectures. Here, we combine chromatin immunoprecipitation sequencing, Hi-C, RNA sequencing, and 40 whole-genome sequences from Heliconius butterflies to show that red color pattern diversification occurred via many genomic loci. We find that the red wing pattern master regulatory transcription factor Optix binds dozens of loci also under selection, which frequently form three-dimensional adaptive hubs with selection acting on multiple physically interacting genes. Many Optix-bound genes under selection are tied to pigmentation and wing development, and these loci collectively maintain separation between adaptive red color pattern phenotypes in natural populations. We propose a model of trait evolution where functional connections between loci may resolve much of the disparity between large-effect and polygenic evolutionary models.

Large-scale mutation in the evolution of a gene complex for cryptic coloration

The types of mutations affecting adaptation in the wild are only beginning to be understood. In particular, whether structural changes shape adaptation by suppressing recombination or by creating new mutations is unresolved. Here, we show that multiple linked but recombining loci underlie cryptic color morphs of Timema chumash stick insects. In a related species, these loci are found in a region of suppressed recombination, forming a supergene. However, in seven species of Timema, we found that a megabase-size "supermutation" has deleted color loci in green morphs. Moreover, we found that balancing selection likely contributes more to maintaining this mutation than does introgression. Our results show how suppressed recombination and large-scale mutation can help to package gene complexes into discrete units of diversity such as morphs, ecotypes, or species.

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.

Divergence of chemosensing during the early stages of speciation

Chemosensory communication is essential to insect biology, playing indispensable roles during mate-finding, foraging, and oviposition behaviors. These traits are particularly important during speciation, where chemical perception may serve to establish species barriers. However, identifying genes associated with such complex behavioral traits remains a significant challenge. Through a combination of transcriptomic and genomic approaches, we characterize the genetic architecture of chemoperception and the role of chemosensing during speciation for a young species pair of Heliconius butterflies, Heliconius melpomene and Heliconius cydno We provide a detailed description of chemosensory gene-expression profiles as they relate to sensory tissue (antennae, legs, and mouthparts), sex (male and female), and life stage (unmated and mated female butterflies). Our results untangle the potential role of chemical communication in establishing barriers during speciation and identify strong candidate genes for mate and host plant choice behaviors. Of the 252 chemosensory genes, HmOBP20 (involved in volatile detection) and HmGr56 (a putative synephrine-related receptor) emerge as strong candidates for divergence in pheromone detection and host plant discrimination, respectively. These two genes are not physically linked to wing-color pattern loci or other genomic regions associated with visual mate preference. Altogether, our results provide evidence for chemosensory divergence between H. melpomene and H. cydno, two rarely hybridizing butterflies with distinct mate and host plant preferences, a finding that supports a polygenic architecture of species boundaries.

Speciation through chromosomal fusion and fission in Lepidoptera

Changes in chromosome numbers may strongly affect reproductive barriers, because individuals heterozygous for distinct karyotypes are typically expected to be at least partially sterile or to show reduced recombination. Therefore, several classic speciation models are based on chromosomal changes. One import mechanism generating variation in chromosome numbers is fusion and fission of existing chromosomes, which is particularly likely in species with holocentric chromosomes, i.e. chromosomes that lack a single centromere. Holocentric chromosomes evolved repeatedly across the tree of life, including in Lepidoptera. Although changes in chromosome numbers are hypothesized to be an important driver of the spectacular diversification of Lepidoptera, comparative studies across the order are lacking. We performed the first comprehensive literature survey of karyotypes for Lepidoptera species since the 1970s and tested if, and how, chromosomal variation might affect speciation. Even though a meta-analysis of karyological differences between closely related taxa did not reveal an effect on the degree of reproductive isolation, phylogenetic diversification rate analyses across the 16 best-covered genera indicated a strong, positive association of rates of chromosome number evolution and speciation. These findings suggest a macroevolutionary impact of varying chromosome numbers in Lepidoptera and likely apply to other taxonomic groups, especially to those with holocentric chromosomes. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

The roles of hybridization and habitat fragmentation in the evolution of Brazil's enigmatic longwing butterflies, Heliconius nattereri and H. hermathena

BACKGROUND

Heliconius butterflies are widely distributed across the Neotropics and have evolved a stunning array of wing color patterns that mediate Müllerian mimicry and mating behavior. Their rapid radiation has been strongly influenced by hybridization, which has created new species and allowed sharing of color patterning alleles between mimetic species pairs. While these processes have frequently been observed in widespread species with contiguous distributions, many Heliconius species inhabit patchy or rare habitats that may strongly influence the origin and spread of species and color patterns. Here, we assess the effects of historical population fragmentation and unique biology on the origins, genetic health, and color pattern evolution of two rare and sparsely distributed Brazilian butterflies, Heliconius hermathena and Heliconius nattereri.

RESULTS

We assembled genomes and re-sequenced whole genomes of eight H. nattereri and 71 H. hermathena individuals. These species harbor little genetic diversity, skewed site frequency spectra, and high deleterious mutation loads consistent with recent population bottlenecks. Heliconius hermathena consists of discrete, strongly isolated populations that likely arose from a single population that dispersed after the last glacial maximum. Despite having a unique color pattern combination that suggested a hybrid origin, we found no genome-wide evidence that H. hermathena is a hybrid species. However, H. hermathena mimicry evolved via introgression, from co-mimetic Heliconius erato, of a small genomic region upstream of the color patterning gene cortex.

CONCLUSIONS

Heliconius hermathena and H. nattereri population fragmentation, potentially driven by historical climate change and recent deforestation, has significantly reduced the genetic health of these rare species. Our results contribute to a growing body of evidence that introgression of color patterning alleles between co-mimetic species appears to be a general feature of Heliconius evolution.

The genomics of coloration provides insights into adaptive evolution

Coloration is an easily quantifiable visual trait that has proven to be a highly tractable system for genetic analysis and for studying adaptive evolution. The application of genomic approaches to evolutionary studies of coloration is providing new insight into the genetic architectures underlying colour traits, including the importance of large-effect mutations and supergenes, the role of development in shaping genetic variation and the origins of adaptive variation, which often involves adaptive introgression. Improved knowledge of the genetic basis of traits can facilitate field studies of natural selection and sexual selection, making it possible for strong selection and its influence on the genome to be demonstrated in wild populations.

De Novo Genome Assembly of the Meadow Brown Butterfly, Maniola jurtina

Meadow brown butterflies (Maniola jurtina) on the Isles of Scilly represent an ideal model in which to dissect the links between genotype, phenotype and long-term patterns of selection in the wild - a largely unfulfilled but fundamental aim of modern biology. To meet this aim, a clear description of genotype is required. Here we present the draft genome sequence of M. jurtina to serve as a founding genetic resource for this species. Seven libraries were constructed using pooled DNA from five wild caught spotted females and sequenced using Illumina, PacBio RSII and MinION technology. A novel hybrid assembly approach was employed to generate a final assembly with an N50 of 214 kb (longest scaffold 2.9 Mb). The sequence assembly described here predicts a gene count of 36,294 and includes variants and gene duplicates from five genotypes. Core BUSCO (Benchmarking Universal Single-Copy Orthologs) gene sets of Arthropoda and Insecta recovered 90.5% and 88.7% complete and single-copy genes respectively. Comparisons with 17 other Lepidopteran species placed 86.5% of the assembled genes in orthogroups. Our results provide the first high-quality draft genome and annotation of the butterfly M. jurtina.

Evolution of wing pigmentation in Drosophila: Diversity, physiological regulation, and cis-regulatory evolution

Fruit flies (Drosophila and its close relatives, or “drosophilids”) are a group that includes an important model organism, Drosophila melanogaster, and also very diverse species distributed worldwide. Many of these species have black or brown pigmentation patterns on their wings, and have been used as material for evo‐devo research. Pigmentation patterns are thought to have evolved rapidly compared with body plans or body shapes; hence they are advantageous model systems for studying evolutionary gains of traits and parallel evolution. Various groups of drosophilids, including genus Idiomyia (Hawaiian Drosophila), have a variety of pigmentations, ranging from simple black pigmentations around crossveins to a single antero‐distal spot and a more complex mottled pattern. Pigmentation patterns are sometimes obviously used for sexual displays; however, in some cases they may have other functions. The process of wing formation in Drosophila, the general mechanism of pigmentation formation, and the transport of substances necessary for pigmentation, including melanin precursors, through wing veins are summarized here. Lastly, the evolution of the expression of genes regulating pigmentation patterns, the role of cis‐regulatory regions, and the conditions required for the evolutionary emergence of pigmentation patterns are discussed. Future prospects for research on the evolution of wing pigmentation pattern formation in drosophilids are presented, particularly from the point of view of how they compare with other studies of the evolution of new traits.

Contrasting genomic and phenotypic outcomes of hybridization between pairs of mimetic butterfly taxa across a suture zone

Hybrid zones, whereby divergent lineages come into contact and eventually hybridize, can provide insights on the mechanisms involved in population differentiation and reproductive isolation, and ultimately speciation. Suture zones offer the opportunity to compare these processes across multiple species. In this paper we use reduced-complexity genomic data to compare the genetic and phenotypic structure and hybridization patterns of two mimetic butterfly species, Ithomia salapia and Oleria onega (Nymphalidae: Ithomiini), each consisting of a pair of lineages differentiated for their wing colour pattern and that come into contact in the Andean foothills of Peru. Despite similarities in their life history, we highlight major differences, both at the genomic and phenotypic level, between the two species. These differences include the presence of hybrids, variations in wing phenotype, and genomic patterns of introgression and differentiation. In I. salapia, the two lineages appear to hybridize only rarely, whereas in O. onega the hybrids are not only more common, but also genetically and phenotypically more variable. We also detected loci statistically associated with wing colour pattern variation, but in both species these loci were not over-represented among the candidate barrier loci, suggesting that traits other than wing colour pattern may be important for reproductive isolation. Our results contrast with the genomic patterns observed between hybridizing lineages in the mimetic Heliconius butterflies, and call for a broader investigation into the genomics of speciation in Ithomiini - the largest radiation of mimetic butterflies.

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.

Whole-chromosome hitchhiking driven by a male-killing endosymbiont

Neo-sex chromosomes are found in many taxa, but the forces driving their emergence and spread are poorly understood. The female-specific neo-W chromosome of the African monarch (or queen) butterfly Danaus chrysippus presents an intriguing case study because it is restricted to a single 'contact zone' population, involves a putative colour patterning supergene, and co-occurs with infection by the male-killing endosymbiont Spiroplasma. We investigated the origin and evolution of this system using whole genome sequencing. We first identify the 'BC supergene', a broad region of suppressed recombination across nearly half a chromosome, which links two colour patterning loci. Association analysis suggests that the genes yellow and arrow in this region control the forewing colour pattern differences between D. chrysippus subspecies. We then show that the same chromosome has recently formed a neo-W that has spread through the contact zone within approximately 2,200 years. We also assembled the genome of the male-killing Spiroplasma, and find that it shows perfect genealogical congruence with the neo-W, suggesting that the neo-W has hitchhiked to high frequency as the male-killer has spread through the population. The complete absence of female crossing-over in the Lepidoptera causes whole-chromosome hitchhiking of a single neo-W haplotype, carrying a single allele of the BC supergene and dragging multiple non-synonymous mutations to high frequency. This has created a population of infected females that all carry the same recessive colour patterning allele, making the phenotypes of each successive generation highly dependent on uninfected male immigrants. Our findings show how hitchhiking can occur between the physically unlinked genomes of host and endosymbiont, with dramatic consequences.

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.

The Genomic Complexity of a Large Inversion in Great Tits

Chromosome inversions have clear effects on genome evolution and have been associated with speciation, adaptation, and the evolution of the sex chromosomes. In birds, these inversions may play an important role in hybridization of species and disassortative mating. We identified a large (≈64 Mb) inversion polymorphism in the great tit (Parus major) that encompasses almost 1,000 genes and more than 90% of Chromosome 1A. The inversion occurs at a low frequency in a set of over 2,300 genotyped great tits in the Netherlands with only 5% of the birds being heterozygous for the inversion. In an additional analysis of 29 resequenced birds from across Europe, we found two heterozygotes. The likely inversion breakpoints show considerable genomic complexity, including multiple copy number variable segments. We identified different haplotypes for the inversion, which differ in the degree of recombination in the center of the chromosome. Overall, this remarkable genetic variant is widespread among distinct great tit populations and future studies of the inversion haplotype, including how it affects the fitness of carriers, may help to understand the mechanisms that maintain it.

Symmetry systems on the wings of Dichromodes Guenée (Lepidoptera: Geometridae) are unconstrained by venation

The nymphalid groundplan, an idealized schematic illustrating the essential elements of butterfly wing patterns, predicts a consistent relationship between color pattern and wing venation. Moths in the family Geometridae have wing shapes and patterns that often resemble those of butterflies, and until recently, this family was believed to be among butterflies’ closest relatives. However, an examination of the geometrid genus Dichromodes Guenée, 1858 shows no consistent relationship between the central symmetry system and wing venation. Whereas the distal edge of the central symmetry system is predicted to reach the costal margin proximal to the Subcostal vein in butterflies and acronictine moths, it has no consistent relationship with the Subcostal, Radius, or Radial Sector 1 veins in Dichromodes. This finding highlights developmental diversity that was previously overlooked due to the overwhelming preference for butterflies in studies of lepidopteran wing patterns.

Transcriptional Landscapes of Divergent Sporophyte Development in Two Mosses, Physcomitrium (Physcomitrella) patens and Funaria hygrometrica

Understanding the molecular basis of morphological shifts is a fundamental question of evolutionary biology. New morphologies may arise through the birth/death of genes (gene gain/loss) or by reutilizing existing gene sets. Yet, the relative contribution of these two processes to radical morphological shifts is still poorly understood. Here, we use the model system of two mosses, Funaria hygrometrica and Physcomitrium (Physcomitrella) patens, to investigate the molecular mechanisms underlying contrasting sporophyte architectures. We used comparative analysis of time-series expression data for four stages of sporophyte development in both species to address this question in detail. We found that large-scale differences in sporophytic architecture are mainly governed by orthologous (i.e., shared) genes frequently experiencing temporal gene expression shifts between the two species. While the absolute number of species-specific genes expressed during sporophyte development is somewhat smaller, we observed a significant increase of their proportion in preferentially sporophyte expressed genes, suggesting a fundamental role in the sporophyte phase. However, further functional studies are necessary to determine their contribution to diverging sporophyte morphologies. Our results add to the growing set of studies suggesting that radical changes in morphology may rely on the heterochronic expression of conserved regulators.

A transposable element insertion is associated with an alternative life history strategy

Tradeoffs affect resource allocation during development and result in fitness consequences that drive the evolution of life history strategies. Yet despite their importance, we know little about the mechanisms underlying life history tradeoffs. Many species of Colias butterflies exhibit an alternative life history strategy (ALHS) where females divert resources from wing pigment synthesis to reproductive and somatic development. Due to this reallocation, a wing color polymorphism is associated with the ALHS: either yellow/orange or white. Here we map the locus associated with this ALHS in Colias crocea to a transposable element insertion located downstream of the Colias homolog of BarH-1, a homeobox transcription factor. Using CRISPR/Cas9 gene editing, antibody staining, and electron microscopy we find white-specific expression of BarH-1 suppresses the formation of pigment granules in wing scales and gives rise to white wing color. Lipid and transcriptome analyses reveal physiological differences associated with the ALHS. Together, these findings characterize a mechanism for a female-limited ALHS.

Tradeoffs are central to life history theory and evolutionary biology, yet almost nothing is known about their mechanistic basis. Here the authors characterize one such mechanism and find a transposable element insertion is associated with the switch between alternative life history strategies.

Evolutionary genetics in insect phenotypic radiations: the value of a comparative genomic approach

Expanding genome sequencing and transgenic technologies are enabling the discovery of genes driving phenotypic diversity across insect taxa. Limitations in downstream functional genetic approaches, however, have been an obstacle for developing non-model systems for evolutionary genetics. Phenotypically diverse radiations, such as those exhibiting convergence and divergence as a result of mimicry, are ideal for evolutionary genetics as they can lead to insights using comparative genomic approaches alone. The varied and repeated instances of phenotypes in highly polymorphic systems allow assessment of whether similar loci are repeatedly targeted by selection and can inform how alleles sort across lineages. Comparative genomics of these taxa can be used to decipher components of gene regulatory networks, dissect regulatory regions, and validate genes.

Altitude and life-history shape the evolution of Heliconius wings

Phenotypic divergence between closely related species has long interested biologists. Taxa that inhabit a range of environments and have diverse natural histories can help understand how selection drives phenotypic divergence. In butterflies, wing color patterns have been extensively studied but diversity in wing shape and size is less well understood. Here, we assess the relative importance of phylogenetic relatedness, natural history, and habitat on shaping wing morphology in a large dataset of over 3500 individuals, representing 13 Heliconius species from across the Neotropics. We find that both larval and adult behavioral ecology correlate with patterns of wing sexual dimorphism and adult size. Species with solitary larvae have larger adult males, in contrast to gregarious Heliconius species, and indeed most Lepidoptera, where females are larger. Species in the pupal-mating clade are smaller than those in the adult-mating clade. Interestingly, we find that high-altitude species tend to have rounder wings and, in one of the two major Heliconius clades, are also bigger than their lowland relatives. Furthermore, within two widespread species, we find that high-altitude populations also have rounder wings. Thus, we reveal novel adaptive wing morphological divergence among Heliconius species beyond that imposed by natural selection on aposematic wing coloration.

Parallel evolution of ancient, pleiotropic enhancers underlies butterfly wing pattern mimicry

Color pattern mimicry in Heliconius butterflies is a classic case study of complex trait adaptation via selection on a few large effect genes. Association studies have linked color pattern variation to a handful of noncoding regions, yet the presumptive cis-regulatory elements (CREs) that control color patterning remain unknown. Here we combine chromatin assays, DNA sequence associations, and genome editing to functionally characterize 5 cis-regulatory elements of the color pattern gene optix We were surprised to find that the cis-regulatory architecture of optix is characterized by pleiotropy and regulatory fragility, where deletion of individual cis-regulatory elements has broad effects on both color pattern and wing vein development. Remarkably, we found orthologous cis-regulatory elements associate with wing pattern convergence of distantly related comimics, suggesting that parallel coevolution of ancestral elements facilitated pattern mimicry. Our results support a model of color pattern evolution in Heliconius where changes to ancient, multifunctional cis-regulatory elements underlie adaptive radiation.

Genomic architecture and introgression shape a butterfly radiation

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.

Genetic convergence of industrial melanism in three geometrid moths

The rise of dark (melanic) forms of many species of moth in heavily coal-polluted areas of nineteenth- and twentieth-century Britain, and their post-1970s fall, point to a common selective pressure (camouflage against bird predators) acting at the community level. The extent to which this convergent phenotypic response relied on similar genetic and developmental mechanisms is unknown. We examine this problem by testing the hypothesis that the locus controlling melanism in Phigalia pilosaria and Odontopera bidentata, two species of geometrid moth that showed strong associations between melanism and coal pollution, is the same as that controlling melanism in Biston betularia, previously identified as the gene cortex. Comparative linkage mapping using family material supports the hypothesis for both species, indicating a deeply conserved developmental mechanism for melanism involving cortex. However, in contrast to the strong selective sweep signature seen in British B. betularia, no significant association was detected between cortex-region markers and melanic morphs in wild-caught samples of P. pilosaria and O. bidentata, implying much older, or diverse, origins of melanic morph alleles in these latter species.

Wing scale ultrastructure underlying convergent and divergent iridescent colours in mimetic Heliconius butterflies

Iridescence is an optical phenomenon whereby colour changes with the illumination and viewing angle. It can be produced by thin film interference or diffraction. Iridescent optical structures are fairly common in nature, but relatively little is known about their production or evolution. Here we describe the structures responsible for producing blue-green iridescent colour in Heliconius butterflies. Overall the wing scale structures of iridescent and non-iridescent Heliconius species are very similar, both having longitudinal ridges joined by cross-ribs. However, iridescent scales have ridges composed of layered lamellae, which act as multilayer reflectors. Differences in brightness between species can be explained by the extent of overlap of the lamellae and their curvature as well as the density of ridges on the scale. Heliconius are well known for their Müllerian mimicry. We find that iridescent structural colour is not closely matched between co-mimetic species. Differences appear less pronounced in models of Heliconius vision than models of avian vision, suggesting that they are not driven by selection to avoid heterospecific courtship by co-mimics. Ridge profiles appear to evolve relatively slowly, being similar between closely related taxa, while ridge density evolves faster and is similar between distantly related co-mimics.

Genomic evidence of human selection on Vavilovian mimicry

Vavilovian mimicry is an evolutionary process by which weeds evolve to resemble domesticated crop plants and is thought to be the result of unintentional selection by humans. Unravelling its molecular mechanisms will extend our knowledge of mimicry and contribute to our understanding of the origin and evolution of agricultural weeds, an important component of crop biology. To this end, we compared mimetic and non-mimetic populations of Echinochloa crus-galli from the Yangtze River basin phenotypically and by genome resequencing, and we show that this weed in rice paddies has evolved a small tiller angle, allowing it to phenocopy cultivated rice at the seedling stage. We demonstrate that mimetic lines evolved from the non-mimetic population as recently as 1,000 yr ago and were subject to a genetic bottleneck, and that genomic regions containing 87 putative plant architecture-related genes (including LAZY1, a key gene controlling plant tiller angle) were under selection during the mimicry process. Our data provide genome-level evidence for the action of human selection on Vavilovian mimicry.

The population genetics of crypsis in vertebrates: recent insights from mice, hares, and lizards

By combining well-established population genetic theory with high-throughput sequencing data from natural populations, major strides have recently been made in understanding how, why, and when vertebrate populations evolve crypsis. Here, we focus on background matching, a particular facet of crypsis that involves the ability of an organism to conceal itself through matching its color to the surrounding environment. While interesting in and of itself, the study of this phenotype has also provided fruitful population genetic insights into the interplay of strong positive selection with other evolutionary processes. Specifically, and predicated upon the findings of previous candidate gene association studies, a primary focus of this recent literature involves the realization that the inference of selection from DNA sequence data first requires a robust model of population demography in order to identify genomic regions which do not conform to neutral expectations. Moreover, these demographic estimates provide crucial information about the origin and timing of the onset of selective pressures associated with, for example, the colonization of a novel environment. Furthermore, such inference has revealed crypsis to be a particularly useful phenotype for investigating the interplay of migration and selection—with examples of gene flow constraining rates of adaptation, or alternatively providing the genetic variants that may ultimately sweep through the population. Here, we evaluate the underlying evidence, review the strengths and weaknesses of the many population genetic methodologies used in these studies, and discuss how these insights have aided our general understanding of the evolutionary process.

Unravelling the genes forming the wing pattern supergene in the polymorphic butterfly Heliconius numata

Background

Unravelling the genetic basis of polymorphic characters is central to our understanding of the origins and diversification of living organisms. Recently, supergenes have been implicated in a wide range of complex polymorphisms, from adaptive colouration in butterflies and fish to reproductive strategies in birds and plants. The concept of a supergene is now a hot topic in biology, and identification of its functional elements is needed to shed light on the evolution of highly divergent adaptive traits. Here, we apply different gene expression analyses to study the supergene P that controls polymorphism of mimetic wing colour patterns in the neotropical butterfly Heliconius numata.

Results

We performed de novo transcriptome assembly and differential expression analyses using high-throughput Illumina RNA sequencing on developing wing discs of different H. numata morphs. Within the P interval, 30 and 17 of the 191 transcripts were expressed differentially in prepupae and day-1 pupae, respectively. Among these is the gene cortex, known to play a role in wing pattern formation in Heliconius and other Lepidoptera. Our in situ hybridization experiments confirmed the relationship between cortex expression and adult wing patterns.

Conclusions

This study found the majority of genes in the P interval to be expressed in the developing wing discs during the critical stages of colour pattern formation, and detect drastic changes in expression patterns in multiple genes associated with structural variants. The patterns of expression of cortex only partially recapitulate the variation in adult phenotype, suggesting that the remaining phenotypic variation could be controlled by other genes within the P interval. Although functional studies on cortex are now needed to determine its exact developmental role, our results are in accordance with the classical supergene hypothesis, whereby several genes inherited together due to tight linkage control a major developmental switch.

Electronic supplementary material

The online version of this article (10.1186/s13227-019-0129-2) contains supplementary material, which is available to authorized users.

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.