Selection drives divergence of eye morphology in sympatric Heliconius butterflies

When populations experience different sensory conditions, natural selection may favor sensory system divergence, affecting peripheral structures and/or downstream neural pathways. We characterized the outer eye morphology of sympatric Heliconius species from different forest types and their first-generation reciprocal hybrids to test for adaptive visual system divergence and hybrid disruption. In Panama, Heliconius cydno occurs in closed forests, whereas Heliconius melpomene resides at the forest edge. Among wild individuals, H. cydno has larger eyes than H. melpomene, and there are heritable, habitat-associated differences in the visual brain structures that exceed neutral divergence expectations. Notably, hybrids have intermediate neural phenotypes, suggesting disruption. To test for similar effects in the visual periphery, we reared both species and their hybrids in common garden conditions. We confirm that H. cydno has larger eyes and provide new evidence that this is driven by selection. Hybrid eye morphology is more H. melpomene-like despite body size being intermediate, contrasting with neural trait intermediacy. Overall, our results suggest that eye morphology differences between H. cydno and H. melpomene are adaptive, and that hybrids may suffer fitness costs due to a mismatch between the peripheral visual structures and previously described neural traits that could affect visual performance.

Transposable Element Insertions Are Associated with Batesian Mimicry in the Pantropical Butterfly Hypolimnas misippus

Hypolimnas misippus is a Batesian mimic of the toxic African Queen butterfly (Danaus chrysippus). Female H. misippus butterflies use two major wing patterning loci (M and A) to imitate three color morphs of D. chrysippus found in different regions of Africa. In this study, we examine the evolution of the M locus and identify it as an example of adaptive atavism. This phenomenon involves a morphological reversion to an ancestral character that results in an adaptive phenotype. We show that H. misippus has re-evolved an ancestral wing pattern present in other Hypolimnas species, repurposing it for Batesian mimicry of a D. chrysippus morph. Using haplotagging, a linked-read sequencing technology, and our new analytical tool, Wrath, we discover two large transposable element insertions located at the M locus and establish that these insertions are present in the dominant allele responsible for producing mimetic phenotype. By conducting a comparative analysis involving additional Hypolimnas species, we demonstrate that the dominant allele is derived. This suggests that, in the derived allele, the transposable elements disrupt a cis-regulatory element, leading to the reversion to an ancestral phenotype that is then utilized for Batesian mimicry of a distinct model, a different morph of D. chrysippus. Our findings present a compelling instance of convergent evolution and adaptive atavism, in which the same pattern element has independently evolved multiple times in Hypolimnas butterflies, repeatedly playing a role in Batesian mimicry of diverse model species.

Colour polymorphism associated with a gene duplication in male wood tiger moths

Colour is often used as an aposematic warning signal, with predator learning expected to lead to a single colour pattern within a population. However, there are many puzzling cases where aposematic signals are also polymorphic. The wood tiger moth, Arctia plantaginis, displays bright hindwing colours associated with unpalatability, and males have discrete colour morphs which vary in frequency between localities. In Finland, both white and yellow morphs can be found, and these colour morphs also differ in behavioural and life-history traits. Here, we show that male colour is linked to an extra copy of a yellow family gene that is only present in the white morphs. This white-specific duplication, which we name valkea, is highly upregulated during wing development. CRISPR targeting valkea resulted in editing of both valkea and its paralog, yellow-e, and led to the production of yellow wings. We also characterise the pigments responsible for yellow, white, and black colouration, showing that yellow is partly produced by pheomelanins, while black is dopamine-derived eumelanin. Our results add to a growing number of studies on the genetic architecture of complex and seemingly paradoxical polymorphisms, and the role of gene duplications and structural variation in adaptive evolution.

Evolutionary dynamics of genome size and content during the adaptive radiation of Heliconiini butterflies

Heliconius butterflies, a speciose genus of Müllerian mimics, represent a classic example of an adaptive radiation that includes a range of derived dietary, life history, physiological and neural traits. However, key lineages within the genus, and across the broader Heliconiini tribe, lack genomic resources, limiting our understanding of how adaptive and neutral processes shaped genome evolution during their radiation. Here, we generate highly contiguous genome assemblies for nine Heliconiini, 29 additional reference-assembled genomes, and improve 10 existing assemblies. Altogether, we provide a dataset of annotated genomes for a total of 63 species, including 58 species within the Heliconiini tribe. We use this extensive dataset to generate a robust and dated heliconiine phylogeny, describe major patterns of introgression, explore the evolution of genome architecture, and the genomic basis of key innovations in this enigmatic group, including an assessment of the evolution of putative regulatory regions at the Heliconius stem. Our work illustrates how the increased resolution provided by such dense genomic sampling improves our power to generate and test gene-phenotype hypotheses, and precisely characterize how genomes evolve.

High level of novelty under the hood of convergent evolution

Little is known about the extent to which species use homologous regulatory architectures to achieve phenotypic convergence. By characterizing chromatin accessibility and gene expression in developing wing tissues, we compared the regulatory architecture of convergence between a pair of mimetic butterfly species. Although a handful of color pattern genes are known to be involved in their convergence, our data suggest that different mutational paths underlie the integration of these genes into wing pattern development. This is supported by a large fraction of accessible chromatin being exclusive to each species, including the de novo lineage-specific evolution of a modular optix enhancer. These findings may be explained by a high level of developmental drift and evolutionary contingency that occurs during the independent evolution of mimicry.

Evolution of physical linkage between loci controlling ecological traits and mating preferences

Coupling of multiple barriers to gene-flow, such as divergent local adaptation and reproductive isolation, facilitates speciation. However, alleles at loci that contribute to barrier effects can be dissociated by recombination. Models of linkage between diverging alleles often consider elements that reduce recombination, such as chromosomal inversions and alleles that modify recombination rate between existing loci. In contrast, here, we consider the evolution of linkage due to the close proximity of loci on the same chromosome. Examples of such physical linkage exist in several species, but in other cases, strong associations are maintained without physical linkage. We use an individual-based model to study the conditions under which the physical linkage between loci controlling ecological traits and mating preferences might be expected to evolve. We modelled a single locus controlling an ecological trait that acts also as a mating cue. Mating preferences are controlled by multiple loci, formed by mutations that are randomly placed in the "genome", within varying distances from the ecological trait locus, allowing us to examine which genomic architectures spread across the population. Our model reveals that stronger physical linkage is favoured when mating preferences and selection are weaker. Under such conditions mating among divergent phenotypes is more frequent, and matching ecological trait and mating preference alleles are more likely to become dissociated by recombination, favouring the evolution of genetic linkage. While most theoretical studies on clustering of divergent loci focus on how physical linkage influences speciation, we show how physical linkage itself can arise, establishing conditions that can favour speciation.

A butterfly pan-genome reveals that a large amount of structural variation underlies the evolution of chromatin accessibility

Despite insertions and deletions being the most common structural variants (SVs) found across genomes, not much is known about how much these SVs vary within populations and between closely related species, nor their significance in evolution. To address these questions, we characterized the evolution of indel SVs using genome assemblies of three closely related Heliconius butterfly species. Over the relatively short evolutionary timescales investigated, up to 18.0% of the genome was composed of indels between two haplotypes of an individual Heliconius charithonia butterfly and up to 62.7% included lineage-specific SVs between the genomes of the most distant species (11 Mya). Lineage-specific sequences were mostly characterized as transposable elements (TEs) inserted at random throughout the genome and their overall distribution was similarly affected by linked selection as single nucleotide substitutions. Using chromatin accessibility profiles (i.e., ATAC-seq) of head tissue in caterpillars to identify sequences with potential cis-regulatory function, we found that out of the 31,066 identified differences in chromatin accessibility between species, 30.4% were within lineage-specific SVs and 9.4% were characterized as TE insertions. These TE insertions were localized closer to gene transcription start sites than expected at random and were enriched for sites with significant resemblance to several transcription factor binding sites with known function in neuron development in Drosophila We also identified 24 TE insertions with head-specific chromatin accessibility. Our results show high rates of structural genome evolution that were previously overlooked in comparative genomic studies and suggest a high potential for structural variation to serve as raw material for adaptive evolution.

Exploitation of an ancestral pheromone biosynthetic pathway contributes to diversification in Heliconius butterflies

During courtship, male butterflies of many species produce androconial secretions containing male sex pheromones (MSPs) that communicate species identity and affect female choice. MSPs are thus likely candidates as reproductive barriers, yet their role in speciation remains poorly studied. Although Heliconius butterflies are a model system in speciation, their MSPs have not been investigated from a macroevolutionary perspective. We use GC/MS to characterize male androconial secretions in 33 of the 69 species in the Heliconiini tribe. We found these blends to be species-specific, consistent with a role in reproductive isolation. We detected a burst in blend diversification rate at the most speciose genus, Heliconius; a consequence of Heliconius and Eueides species using a fatty acid (FA) metabolic pathway to unlock more complex blends than basal Heliconiini species, whose secretions are dominated by plant-like metabolites. A comparison of 10 sister species pairs demonstrates a striking positive correlation between blend dissimilarity and range overlap, consistent with character displacement or reinforcement in sympatry. These results demonstrate for the first time that MSP diversification can promote reproductive isolation across this group of butterflies, showcasing how implementation of an ancestral trait, the co-option of the FA metabolic pathway for pheromone production, can facilitate rapid speciation.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

Condition dependence in biosynthesized chemical defenses of an aposematic and mimetic Heliconius butterfly

Aposematic animals advertise their toxicity or unpalatability with bright warning coloration. However, acquiring and maintaining chemical defenses can be energetically costly, and consequent associations with other important traits could shape chemical defense evolution. Here, we have tested whether chemical defenses are involved in energetic trade‐offs with other traits, or whether the levels of chemical defenses are condition dependent, by studying associations between biosynthesized cyanogenic toxicity and a suite of key life‐history and fitness traits in a Heliconius butterfly under a controlled laboratory setting. Heliconius butterflies are well known for the diversity of their warning color patterns and widespread mimicry and can both sequester the cyanogenic glucosides of their Passiflora host plants and biosynthesize these toxins de novo. We find energetically costly life‐history traits to be either unassociated or to show a general positive association with biosynthesized cyanogenic toxicity. More toxic individuals developed faster and had higher mass as adults and a tendency for increased lifespan and fecundity. These results thus indicate that toxicity level of adult butterflies may be dependent on individual condition, influenced by genetic background or earlier conditions, with maternal effects as one strong candidate mechanism. Additionally, toxicity was higher in older individuals, consistent with previous studies indicating accumulation of toxins with age. As toxicity level at death was independent of lifespan, cyanogenic glucoside compounds may have been recycled to release resources relevant for longevity in these long‐living butterflies. Understanding the origins and maintenance of variation in defenses is necessary in building a more complete picture of factors shaping the evolution of aposematic and mimetic systems.

A large deletion at the cortex locus eliminates butterfly wing patterning

As the genetic basis of natural and domesticated variation has been described in recent years, a number of hotspot genes have been repeatedly identified as the targets of selection, Heliconius butterflies display a spectacular diversity of pattern variants in the wild and the genetic basis of these patterns has been well-described. Here, we sought to identify the mechanism behind an unusual pattern variant that is instead found in captivity, the ivory mutant, in which all scales on both the wings and body become white or yellow. Using a combination of autozygosity mapping and coverage analysis from 37 captive individuals, we identify a 78-kb deletion at the cortex wing patterning locus, a gene which has been associated with wing pattern evolution in H. melpomene and 10 divergent lepidopteran species. This deletion is undetected among 458 wild Heliconius genomes samples, and its dosage explains both homozygous and heterozygous ivory phenotypes found in captivity. The deletion spans a large 5′ region of the cortex gene that includes a facultative 5′UTR exon detected in larval wing disk transcriptomes. CRISPR mutagenesis of this exon replicates the wing phenotypes from coding knock-outs of cortex, consistent with a functional role of ivory-deleted elements in establishing scale color fate. Population demographics reveal that the stock giving rise to the ivory mutant has a mixed origin from across the wild range of H. melpomene, and supports a scenario where the ivory mutation occurred after the introduction of cortex haplotypes from Ecuador. Homozygotes for the ivory deletion are inviable while heterozygotes are the targets of artificial selection, joining 40 other examples of allelic variants that provide heterozygous advantage in animal populations under artificial selection by fanciers and breeders. Finally, our results highlight the promise of autozygosity and association mapping for identifying the genetic basis of aberrant mutations in captive insect populations.

Functional genomics of supergene-controlled behavior in the white-throated sparrow

Supergenes are regions of suppressed recombination that may span hundreds of genes and can control variation in key ecological phenotypes. Since genetic analysis is made impossible by the absence of recombination between genes, it has been difficult to establish how individual genes within these regions contribute to supergene-controlled phenotypes. The white-throated sparrow is a classic example in which a supergene controls behavioral differences as well as distinct coloration that determines mate choice. A landmark study now demonstrates that differences between supergene variants in the promoter sequences of a hormone receptor gene change its expression and control changes in behavior. To unambiguously establish the link between genotype and phenotype, the authors used antisense oligonucleotides to alter the level of gene expression in a focal brain region targeted through a cannula. The study showcases a powerful approach to the functional genomic manipulation of a wild vertebrate species.

Identification and Composition of Clasper Scent Gland Components of the Butterfly Heliconius erato and Its Relation to Mimicry

The butterfly Heliconius erato occurs in various mimetic morphs. The male clasper scent gland releases an anti‐aphrodisiac pheromone and additionally contains a complex mixture of up to 350 components, varying between individuals. In 114 samples of five different mimicry groups and their hybrids 750 different compounds were detected by gas chromatography/mass spectrometry (GC/MS). Many unknown components occurred, which were identified using their mass spectra, gas chromatography/infrared spectroscopy (GC/IR)‐analyses, derivatization, and synthesis. Key compounds proved to be various esters of 3‐oxohexan‐1‐ol and (Z)‐3‐hexen‐1‐ol with (S)‐2,3‐dihydrofarnesoic acid, accompanied by a large variety of other esters with longer terpene acids, fatty acids, and various alcohols. In addition, linear terpenes with up to seven uniformly connected isoprene units occur, e. g. farnesylfarnesol. A large number of the compounds have not been reported before from nature. Discriminant analyses of principal components of the gland contents showed that the iridescent mimicry group differs strongly from the other, mostly also separated, mimicry groups. Comparison with data from other species indicated that Heliconius recruits different biosynthetic pathways in a species‐specific manner for semiochemical formation.

Cortex cis-regulatory switches establish scale colour identity and pattern diversity in Heliconius

In Heliconius butterflies, wing colour pattern diversity and scale types are controlled by a few genes of large effect that regulate colour pattern switches between morphs and species across a large mimetic radiation. One of these genes, cortex, has been repeatedly associated with colour pattern evolution in butterflies. Here we carried out CRISPR knockouts in multiple Heliconius species and show that cortex is a major determinant of scale cell identity. Chromatin accessibility profiling and introgression scans identified cis-regulatory regions associated with discrete phenotypic switches. CRISPR perturbation of these regions in black hindwing genotypes recreated a yellow bar, revealing their spatially limited activity. In the H. melpomene/timareta lineage, the candidate CRE from yellow-barred phenotype morphs is interrupted by a transposable element, suggesting that cis-regulatory structural variation underlies these mimetic adaptations. Our work shows that cortex functionally controls scale colour fate and that its cis-regulatory regions control a phenotypic switch in a modular and pattern-specific fashion.

Genomics of altitude-associated wing shape in two tropical butterflies

Understanding how organisms adapt to their local environment is central to evolution. With new whole-genome sequencing technologies and the explosion of data, deciphering the genomic basis of complex traits that are ecologically relevant is becoming increasingly feasible. Here, we studied the genomic basis of wing shape in two Neotropical butterflies that inhabit large geographical ranges. Heliconius butterflies at high elevations have been shown to generally have rounder wings than those in the lowlands. We reared over 1,100 butterflies from 71 broods of H. erato and H. melpomene in common-garden conditions and showed that wing aspect ratio, that is, elongatedness, is highly heritable in both species and that elevation-associated wing aspect ratio differences are maintained. Genome-wide associations with a published data set of 666 whole genomes from across a hybrid zone, uncovered a highly polygenic basis to wing aspect ratio variation in the wild. We identified several genes that have roles in wing morphogenesis or wing aspect ratio variation in Drosophila flies, making them promising candidates for future studies. There was little evidence for molecular parallelism in the two species, with only one shared candidate gene, nor for a role of the four known colour pattern loci, except for optix in H. erato. Thus, we present the first insights into the heritability and genomic basis of within-species wing aspect ratio in two Heliconius species, adding to a growing body of evidence that polygenic adaptation may underlie many ecologically relevant traits.

Rampant Genome-Wide Admixture across the Heliconius Radiation

How frequent is gene flow between species? The pattern of evolution is typically portrayed as a phylogenetic tree, yet gene flow between good species may be an important mechanism in diversification, spreading adaptive traits and leading to a complex pattern of phylogenetic incongruence. This process has thus far been studied mainly among a few closely related species, or in geographically restricted areas such as islands, but not on the scale of a continental radiation. Using a genomic representation of 40 out of 47 species in the genus, we demonstrate that admixture has played a role throughout the evolution of the charismatic Neotropical butterflies Heliconius. Modeling of phylogenetic networks based on the exome uncovers up to 13 instances of interspecific gene flow. Admixture is detected among the relatives of Heliconius erato, as well as between the ancient lineages leading to modern clades. Interspecific gene flow played a role throughout the evolution of the genus, although the process has been most frequent in the clade of Heliconius melpomene and relatives. We identify Heliconius hecalesia and relatives as putative hybrids, including new evidence for introgression at the loci controlling the mimetic wing patterns. Models accounting for interspecific gene flow yield a more complete picture of the radiation as a network, which will improve our ability to study trait evolution in a realistic comparative framework.

Evolutionary and ecological processes influencing chemical defense variation in an aposematic and mimetic Heliconius butterfly

Chemical defences against predators underlie the evolution of aposematic coloration and mimicry, which are classic examples of adaptive evolution. Surprisingly little is known about the roles of ecological and evolutionary processes maintaining defence variation, and how they may feedback to shape the evolutionary dynamics of species. Cyanogenic Heliconius butterflies exhibit diverse warning color patterns and mimicry, thus providing a useful framework for investigating these questions. We studied intraspecific variation in de novo biosynthesized cyanogenic toxicity and its potential ecological and evolutionary sources in wild populations of Heliconius erato along environmental gradients, in common-garden broods and with feeding treatments. Our results demonstrate substantial intraspecific variation, including detectable variation among broods reared in a common garden. The latter estimate suggests considerable evolutionary potential in this trait, although predicting the response to selection is likely complicated due to the observed skewed distribution of toxicity values and the signatures of maternal contributions to the inheritance of toxicity. Larval diet contributed little to toxicity variation. Furthermore, toxicity profiles were similar along steep rainfall and altitudinal gradients, providing little evidence for these factors explaining variation in biosynthesized toxicity in natural populations. In contrast, there were striking differences in the chemical profiles of H. erato from geographically distant populations, implying potential local adaptation in the acquisition mechanisms and levels of defensive compounds. The results highlight the extensive variation and potential for adaptive evolution in defense traits for aposematic and mimetic species, which may contribute to the high diversity often found in these systems.

Population structure, adaptation and divergence of the meadow spittlebug, Philaenus spumarius (Hemiptera, Aphrophoridae), revealed by genomic and morphological data

Understanding patterns of population differentiation and gene flow in insect vectors of plant diseases is crucial for the implementation of management programs of disease. We investigated morphological and genome-wide variation across the distribution range of the spittlebug Philaenus spumarius (Linnaeus, 1758) (Hemiptera, Auchenorrhyncha, Aphrophoridae), presently the most important vector of the plant pathogenic bacterium Xylella fastidiosa Wells et al., 1987 in Europe. We found genome-wide divergence between P. spumarius and a very closely related species, P. tesselatus Melichar, 1899, at RAD sequencing markers. The two species may be identified by the morphology of male genitalia but are not differentiated at mitochondrial COI, making DNA barcoding with this gene ineffective. This highlights the importance of using integrative approaches in taxonomy. We detected admixture between P. tesselatus from Morocco and P. spumarius from the Iberian Peninsula, suggesting gene-flow between them. Within P. spumarius, we found a pattern of isolation-by-distance in European populations, likely acting alongside other factors restricting gene flow. Varying levels of co-occurrence of different lineages, showing heterogeneous levels of admixture, suggest other isolation mechanisms. The transatlantic populations of North America and Azores were genetically closer to the British population analyzed here, suggesting an origin from North-Western Europe, as already detected with mitochondrial DNA. Nevertheless, these may have been produced through different colonization events. We detected SNPs with signatures of positive selection associated with environmental variables, especially related to extremes and range variation in temperature and precipitation. The population genomics approach provided new insights into the patterns of divergence, gene flow and adaptation in these spittlebugs and led to several hypotheses that require further local investigation.

Insights into invasive species from whole-genome resequencing

Studies of invasive species can simultaneously inform management strategies and quantify rapid evolution in the wild. The role of genomics in invasion science is increasingly recognised, and the growing availability of reference genomes for invasive species is paving the way for whole-genome resequencing studies in a wide range of systems. Here, we survey the literature to assess the application of whole-genome resequencing data in invasion biology. For some applications, such as the reconstruction of invasion routes in time and space, sequencing the whole genome of many individuals can increase the accuracy of existing methods. In other cases, population genomic approaches such as haplotype analysis can permit entirely new questions to be addressed and new technologies applied. To date whole-genome resequencing has only been used in a handful of invasive systems, but these studies have confirmed the importance of processes such as balancing selection and hybridization in allowing invasive species to reuse existing adaptations and rapidly overcome the challenges of a foreign ecosystem. The use of genomic data does not constitute a paradigm shift per se, but by leveraging new theory, tools, and technologies, population genomics can provide unprecedented insight into basic and applied aspects of invasion science.

A high-quality, chromosome-level genome assembly of the Black Soldier Fly (Hermetia illucens L.)

Hermetia illucens L. (Diptera: Stratiomyidae), the Black Soldier Fly (BSF) is an increasingly important species for bioconversion of organic material into animal feed. We generated a high-quality chromosome-scale genome assembly of the BSF using Pacific Bioscience, 10X Genomics linked read and high-throughput chromosome conformation capture sequencing technology. Scaffolding the final assembly with Hi-C data produced a highly contiguous 1.01 Gb genome with 99.75% of scaffolds assembled into pseudochromosomes representing seven chromosomes with 16.01 Mb contig and 180.46 Mb scaffold N50 values. The highly complete genome obtained a Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness of 98.6%. We masked 67.32% of the genome as repetitive sequences and annotated a total of 16,478 protein-coding genes using the BRAKER2 pipeline. We analyzed an established lab population to investigate the genomic variation and architecture of the BSF revealing six autosomes and an X chromosome. Additionally, we estimated the inbreeding coefficient (1.9%) of the lab population by assessing runs of homozygosity. This provided evidence for inbreeding events including long runs of homozygosity on chromosome 5. The release of this novel chromosome-scale BSF genome assembly will provide an improved resource for further genomic studies, functional characterization of genes of interest and genetic modification of this economically important species.

Global population genetic structure and demographic trajectories of the black soldier fly, Hermetia illucens

Background

The black soldier fly (Hermetia illucens) is the most promising insect candidate for nutrient-recycling through bioconversion of organic waste into biomass, thereby improving sustainability of protein supplies for animal feed and facilitating transition to a circular economy. Contrary to conventional livestock, genetic resources of farmed insects remain poorly characterised. We present the first comprehensive population genetic characterisation of H. illucens. Based on 15 novel microsatellite markers, we genotyped and analysed 2862 individuals from 150 wild and captive populations originating from 57 countries on seven subcontinents.

Results

We identified 16 well-distinguished genetic clusters indicating substantial global population structure. The data revealed genetic hotspots in central South America and successive northwards range expansions within the indigenous ranges of the Americas. Colonisations and naturalisations of largely unique genetic profiles occurred on all non-native continents, either preceded by demographically independent founder events from various single sources or involving admixture scenarios. A decisive primarily admixed Polynesian bridgehead population serially colonised the entire Australasian region and its secondarily admixed descendants successively mediated invasions into Africa and Europe. Conversely, captive populations from several continents traced back to a single North American origin and exhibit considerably reduced genetic diversity, although some farmed strains carry distinct genetic signatures. We highlight genetic footprints characteristic of progressing domestication due to increasing socio-economic importance of H. illucens, and ongoing introgression between domesticated strains globally traded for large-scale farming and wild populations in some regions.

Conclusions

We document the dynamic population genetic history of a cosmopolitan dipteran of South American origin shaped by striking geographic patterns. These reflect both ancient dispersal routes, and stochastic and heterogeneous anthropogenic introductions during the last century leading to pronounced diversification of worldwide structure of H. illucens. Upon the recent advent of its agronomic commercialisation, however, current human-mediated translocations of the black soldier fly largely involve genetically highly uniform domesticated strains, which meanwhile threaten the genetic integrity of differentiated unique local resources through introgression. Our in-depth reconstruction of the contemporary and historical demographic trajectories of H. illucens emphasises benchmarking potential for applied future research on this emerging model of the prospering insect-livestock sector.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01029-w.

Phenotypic plasticity in chemical defence of butterflies allows usage of diverse host plants

Host plant specialization is a major force driving ecological niche partitioning and diversification in insect herbivores. The cyanogenic defences of Passiflora plants keep most herbivores at bay, but not the larvae of Heliconius butterflies, which can both sequester and biosynthesize cyanogenic compounds. Here, we demonstrate that both Heliconius cydno chioneus and H. melpomene rosina have remarkable plasticity in their chemical defences. When feeding on Passiflora species with cyanogenic compounds that they can readily sequester, both species downregulate the biosynthesis of these compounds. By contrast, when fed on Passiflora plants that do not contain cyanogenic glucosides that can be sequestered, both species increase biosynthesis. This biochemical plasticity comes at a fitness cost for the more specialist H. m. rosina, as adult size and weight for this species negatively correlate with biosynthesis levels, but not for the more generalist H. c. chioneus. By contrast, H. m rosina has increased performance when sequestration is possible on its specialized host plant. In summary, phenotypic plasticity in biochemical responses to different host plants offers these butterflies the ability to widen their range of potential hosts within the Passiflora genus, while maintaining their chemical defences.

A novel terpene synthase controls differences in anti-aphrodisiac pheromone production between closely related Heliconius butterflies

Plants and insects often use the same compounds for chemical communication, but not much is known about the genetics of convergent evolution of chemical signals. The terpene (E)-β-ocimene is a common component of floral scent and is also used by the butterfly Heliconius melpomene as an anti-aphrodisiac pheromone. While the biosynthesis of terpenes has been described in plants and microorganisms, few terpene synthases (TPSs) have been identified in insects. Here, we study the recent divergence of 2 species, H. melpomene and Heliconius cydno, which differ in the presence of (E)-β-ocimene; combining linkage mapping, gene expression, and functional analyses, we identify 2 novel TPSs. Furthermore, we demonstrate that one, HmelOS, is able to synthesise (E)-β-ocimene in vitro. We find no evidence for TPS activity in HcydOS (HmelOS ortholog of H. cydno), suggesting that the loss of (E)-β-ocimene in this species is the result of coding, not regulatory, differences. The TPS enzymes we discovered are unrelated to previously described plant and insect TPSs, demonstrating that chemical convergence has independent evolutionary origins.

Plants and insects often use the same compounds for chemical communication, but little is known about the convergent evolution of such chemical signals. This study identifies a novel terpene synthase involved in production of an anti-aphrodisiac pheromone by the butterfly Heliconius melpomene. This enzyme is unrelated to other insect terpene synthases, providing evidence that the ability to synthesise terpenes has arisen multiple times independently within the insects.

Clustering of loci controlling species differences in male chemical bouquets of sympatric Heliconius butterflies

The degree to which loci promoting reproductive isolation cluster in the genome-that is, the genetic architecture of reproductive isolation-can influence the tempo and mode of speciation. Tight linkage between these loci can facilitate speciation in the face of gene flow. Pheromones play a role in reproductive isolation in many Lepidoptera species, and the role of endogenously produced compounds as secondary metabolites decreases the likelihood of pleiotropy associated with many barrier loci. Heliconius butterflies use male sex pheromones to both court females (aphrodisiac wing pheromones) and ward off male courtship (male-transferred antiaphrodisiac genital pheromones), and it is likely that these compounds play a role in reproductive isolation between Heliconius species. Using a set of backcross hybrids between H. melpomene and H. cydno, we investigated the genetic architecture of putative male pheromone compound production. We found a set of 40 significant quantitative trait loci (QTL) representing 33 potential pheromone compounds. QTL clustered significantly on two chromosomes, chromosome 8 for genital compounds and chromosome 20 for wing compounds, and chromosome 20 was enriched for potential pheromone biosynthesis genes. There was minimal overlap between pheromone QTL and known QTL for mate choice and color pattern. Nonetheless, we did detect linkage between a QTL for wing androconial area and optix, a color pattern locus known to play a role in reproductive isolation in these species. This tight clustering of putative pheromone loci might contribute to coincident reproductive isolating barriers, facilitating speciation despite ongoing gene flow.

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.

Adaptive Introgression across Semipermeable Species Boundaries between Local Helicoverpa zea and Invasive Helicoverpa armigera Moths

Hybridization between invasive and native species has raised global concern, given the dramatic increase in species range shifts and pest outbreaks due to anthropogenic dispersal. Nevertheless, secondary contact between sister lineages of local and invasive species provides a natural laboratory to understand the factors that determine introgression and the maintenance or loss of species barriers. Here, we characterize the early evolutionary outcomes following secondary contact between invasive Helicoverpa armigera and native H. zea in Brazil. We carried out whole-genome resequencing of Helicoverpa moths from Brazil in two temporal samples: during the outbreak of H. armigera in 2013 and 2017. There is evidence for a burst of hybridization and widespread introgression from local H. zea into invasive H. armigera coinciding with H. armigera expansion in 2013. However, in H. armigera, the admixture proportion and the length of introgressed blocks were significantly reduced between 2013 and 2017, suggesting selection against admixture. In contrast to the genome-wide pattern, there was striking evidence for adaptive introgression of a single region from the invasive H. armigera into local H. zea, including an insecticide resistance allele that increased in frequency over time. In summary, despite extensive gene flow after secondary contact, the species boundaries are largely maintained except for the single introgressed region containing the insecticide-resistant locus. We document the worst-case scenario for an invasive species, in which there are now two pest species instead of one, and the native species has acquired resistance to pyrethroid insecticides through introgression.

A haplotype-resolved, de novo genome assembly for the wood tiger moth (Arctia plantaginis) through trio binning

BACKGROUND

Diploid genome assembly is typically impeded by heterozygosity because it introduces errors when haplotypes are collapsed into a consensus sequence. Trio binning offers an innovative solution that exploits heterozygosity for assembly. Short, parental reads are used to assign parental origin to long reads from their F1 offspring before assembly, enabling complete haplotype resolution. Trio binning could therefore provide an effective strategy for assembling highly heterozygous genomes, which are traditionally problematic, such as insect genomes. This includes the wood tiger moth (Arctia plantaginis), which is an evolutionary study system for warning colour polymorphism.

FINDINGS

We produced a high-quality, haplotype-resolved assembly for Arctia plantaginis through trio binning. We sequenced a same-species family (F1 heterozygosity ∼1.9%) and used parental Illumina reads to bin 99.98% of offspring Pacific Biosciences reads by parental origin, before assembling each haplotype separately and scaffolding with 10X linked reads. Both assemblies are contiguous (mean scaffold N50: 8.2 Mb) and complete (mean BUSCO completeness: 97.3%), with annotations and 31 chromosomes identified through karyotyping. We used the assembly to analyse genome-wide population structure and relationships between 40 wild resequenced individuals from 5 populations across Europe, revealing the Georgian population as the most genetically differentiated with the lowest genetic diversity.

CONCLUSIONS

We present the first invertebrate genome to be assembled via trio binning. This assembly is one of the highest quality genomes available for Lepidoptera, supporting trio binning as a potent strategy for assembling heterozygous genomes. Using our assembly, we provide genomic insights into the geographic population structure of A. plantaginis.

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.

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.

Chemical signals act as the main reproductive barrier between sister and mimetic Heliconius butterflies

Colour pattern is the main trait that drives mate recognition between Heliconius species that are phylogenetically close. However, when this cue is compromised such as in cases of mimetic, sympatric and closely related species, alternative mating signals must evolve to ensure reproductive isolation and species integrity. The closely related species Heliconius melpomene malleti and H. timareta florencia occur in the same geographical region, and despite being co-mimics, they display strong reproductive isolation. In order to test which cues differ between species, and potentially contribute to reproductive isolation, we quantified differences in the wing phenotype and the male chemical profile. As expected, the wing colour pattern was indistinguishable between the two species, while the chemical profile of the androconial and genital males' extracts showed marked differences. We then conducted behavioural experiments to study the importance of these signals in mate recognition by females. In agreement with our previous results, we found that chemical blends and not wing colour pattern drive the preference of females for conspecific males. Also, experiments with hybrid males and females suggested an important genetic component for both chemical production and preference. Altogether, these results suggest that chemicals are the major reproductive barrier opposing gene flow between these two sister and co-mimic species.

Species specificity and intraspecific variation in the chemical profiles of Heliconius butterflies across a large geographic range

In many animals, mate choice is important for the maintenance of reproductive isolation between species. Traits important for mate choice and behavioral isolation are predicted to be under strong stabilizing selection within species; however, such traits can also exhibit variation at the population level driven by neutral and adaptive evolutionary processes. Here, we describe patterns of divergence among androconial and genital chemical profiles at inter- and intraspecific levels in mimetic Heliconius butterflies. Most variation in chemical bouquets was found between species, but there were also quantitative differences at the population level. We found a strong correlation between interspecific chemical and genetic divergence, but this correlation varied in intraspecific comparisons. We identified "indicator" compounds characteristic of particular species that included compounds already known to elicit a behavioral response, suggesting an approach for identification of candidate compounds for future behavioral studies in novel systems. Overall, the strong signal of species identity suggests a role for these compounds in species recognition, but with additional potentially neutral variation at the population level.

Hybridization and transgressive exploration of colour pattern and wing morphology in Heliconius butterflies

Hybridization can generate novel phenotypes distinct from those of parental lineages, a phenomenon known as transgressive trait variation. Transgressive phenotypes might negatively or positively affect hybrid fitness, and increase available variation. Closely related species of Heliconius butterflies regularly produce hybrids in nature, and hybridization is thought to play a role in the diversification of novel wing colour patterns despite strong stabilizing selection due to interspecific mimicry. Here, we studied wing phenotypes in first- and second-generation hybrids produced by controlled crosses between either two co-mimetic species of Heliconius or between two nonmimetic species. We quantified wing size, shape and colour pattern variation and asked whether hybrids displayed transgressive wing phenotypes. Discrete traits underlain by major-effect loci, such as the presence or absence of colour patches, generate novel phenotypes. For quantitative traits, such as wing shape or subtle colour pattern characters, hybrids only exceed the parental range in specific dimensions of the morphological space. Overall, our study addresses some of the challenges in defining and measuring phenotypic transgression for multivariate traits and our data suggest that the extent to which transgressive trait variation in hybrids contributes to phenotypic diversity depends on the complexity and the genetic architecture of the traits.

Microclimate buffering and thermal tolerance across elevations in a tropical butterfly

Microclimatic variability in tropical forests plays a key role in shaping species distributions and their ability to cope with environmental change, especially for ectotherms. Nonetheless, currently available climatic datasets lack data from the forest interior and, furthermore, our knowledge of thermal tolerance among tropical ectotherms is limited. We therefore studied natural variation in the microclimate experienced by tropical butterflies in the genus Heliconius across their Andean range in a single year. We found that the forest strongly buffers temperature and humidity in the understorey, especially in the lowlands, where temperatures are more extreme. There were systematic differences between our yearly records and macroclimate databases (WorldClim2), with lower interpolated minimum temperatures and maximum temperatures higher than expected. We then assessed thermal tolerance of 10 Heliconius butterfly species in the wild and found that populations at high elevations had significantly lower heat tolerance than those at lower elevations. However, when we reared populations of the widespread H. erato from high and low elevations in a common-garden environment, the difference in heat tolerance across elevations was reduced, indicating plasticity in this trait. Microclimate buffering is not currently captured in publicly available datasets, but could be crucial for enabling upland shifting of species sensitive to heat such as highland Heliconius Plasticity in thermal tolerance may alleviate the effects of global warming on some widespread ectotherm species, but more research is needed to understand the long-term consequences of plasticity on populations and species.

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.

A major locus controls a biologically active pheromone component in Heliconius melpomene

Understanding the production, response, and genetics of signals used in mate choice can inform our understanding of the evolution of both intraspecific mate choice and reproductive isolation. Sex pheromones are important for courtship and mate choice in many insects, but we know relatively little of their role in butterflies. The butterfly Heliconius melpomene uses a complex blend of wing androconial compounds during courtship. Electroantennography in H. melpomene and its close relative Heliconius cydno showed that responses to androconial extracts were not species specific. Females of both species responded equally strongly to extracts of both species, suggesting conservation of peripheral nervous system elements across the two species. Individual blend components provoked little to no response, with the exception of octadecanal, a major component of the H. melpomene blend. Supplementing octadecanal on the wings of octadecanal-rich H. melpomene males led to an increase in the time until mating, demonstrating the bioactivity of octadecanal in Heliconius. Using quantitative trait locus (QTL) mapping, we identified a single locus on chromosome 20 responsible for 41% of the parental species' difference in octadecanal production. This QTL does not overlap with any of the major wing color or mate choice loci, nor does it overlap with known regions of elevated or reduced FST . A set of 16 candidate fatty acid biosynthesis genes lies underneath the QTL. Pheromones in Heliconius carry information relevant for mate choice and are under simple genetic control, suggesting they could be important during speciation.

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.

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.

Suppression of Wolbachia-mediated male-killing in the butterfly Hypolimnas bolina involves a single genomic region

Background

Sex ratio distorting agents (maternally inherited symbionts and meiotically-driving sex chromosomes) are common in insects. When these agents rise to high frequencies they create strong population sex ratio bias and selection then favours mutations that act to restore the rare sex. Despite this strong selection pressure, the evolution of mutations that suppress sex ratio distorting elements appears to be constrained in many cases, where sex-biased populations persist for many generations. This scenario has been observed in the butterfly Hypolimnas bolina, where Wolbachia-mediated male killing endured for 800–1,000 generations across multiple populations before the evolution of suppression. Here we test the hypothesis that this evolutionary lag is the result of suppression being a multilocus trait requiring multiple mutations.

Methods

We developed genetic markers, based on conservation of synteny, for each H. bolina chromosome and verified coverage using recombinational mapping. We then used a Wolbachia-infected mapping family to assess each chromosome for the presence of loci required for male survival, as determined by the presence of markers in all surviving sons.

Results

Informative markers were obtained for each of the 31 chromosomes in H. bolina. The only marker that cosegregated with suppression was located on chromosome 25. A genomic region necessary for suppression has previously been located on this chromosome. We therefore conclude that a single genomic region of the H. bolina genome is necessary for male-killing suppression.

Discussion

The evolutionary lag observed in our system is not caused by a need for changes at multiple genomic locations. The findings favour hypotheses in which either multiple mutations are required within a single genomic region, or the suppressor mutation is a singularly rare event.

Can genomics shed light on the origin of species?

Evolutionary biologists are increasingly using population genetic variation across genomes to address questions around the origin and ongoing evolution of species. Patterns of differentiation between closely related species are highly variable across the genome, and a wide variety of processes contribute to that variation. There is an emerging pattern of parallelism, whereby different species pairs in groups of related species show similar differentiation patterns across their genomes, offering an opportunity to test hypotheses regarding the processes underlying species differentiation. A recent study used both simulations and empirical data to investigate different forms of selection in a radiation of monkeyflowers. The parallel patterns emerged very rapidly after divergence and could not be readily explained by selection for removal of deleterious mutations but instead likely results from some combination of adaptive evolution, species incompatibilities, and ongoing gene flow. Overall, an emerging pattern is that there may be a surprising degree of predictability in the genetic architecture of species differences across groups of related species.

Genetic variation across genomes can be used to address long-standing questions over the origin and ongoing evolution of species. This Primer explores the implications of one such study that is suggesting a surprising repeatability in patterns of genetic divergence across genomes and between species.

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.

Male pheromone composition depends on larval but not adult diet in Heliconius melpomene

1. Condition-dependent traits can act as honest signals of mate quality, with fitter individuals being able to display preferred phenotypes. Nutrition is known to be an important determinant of individual condition, with diet known to affect many secondary sexual traits. 2. In Heliconius butterflies, male chemical signalling plays an important role in female mate choice. Potential male sex pheromone components have been identified previously, although it is unclear what information they convey to the female. 3. In the present study, the effect of diet on androconial and genital compound production is tested in male Heliconius melpomene rosina. To manipulate larval diet, larvae are reared on three different Passiflora host plants: Passiflora menispermifolia, the preferred host plant, Passiflora vitifolia and Passiflora platyloba. To manipulate adult diet, adult butterflies are reared with and without access to pollen, a key component of their diet. 4. No evidence is found to suggest that adult pollen consumption affects compound production in the first 10 days after eclosion. There is also a strong overlap in the chemical profiles of individuals reared on different larval host plants. The most abundant compounds produced by the butterflies do not differ between host plant groups. However, some compounds found in small amounts differ both qualitatively and quantitatively. Some of these compounds are predicted to be of plant origin and the others synthesised by the butterfly. Further electrophysiological and behavioural experiments will be needed to determine the biological significance of these differences.

Genetic dissection of assortative mating behavior

The evolution of new species is made easier when traits under divergent ecological selection are also mating cues. Such ecological mating cues are now considered more common than previously thought, but we still know little about the genetic changes underlying their evolution or more generally about the genetic basis for assortative mating behaviors. Both tight physical linkage and the existence of large-effect preference loci will strengthen genetic associations between behavioral and ecological barriers, promoting the evolution of assortative mating. The warning patterns of Heliconius melpomene and H. cydno are under disruptive selection due to increased predation of nonmimetic hybrids and are used during mate recognition. We carried out a genome-wide quantitative trait locus (QTL) analysis of preference behaviors between these species and showed that divergent male preference has a simple genetic basis. We identify three QTLs that together explain a large proportion (approximately 60%) of the difference in preference behavior observed between the parental species. One of these QTLs is just 1.2 (0-4.8) centiMorgans (cM) from the major color pattern gene optix, and, individually, all three have a large effect on the preference phenotype. Genomic divergence between H. cydno and H. melpomene is high but broadly heterogenous, and admixture is reduced at the preference-optix color pattern locus but not the other preference QTLs. The simple genetic architecture we reveal will facilitate the evolution and maintenance of new species despite ongoing gene flow by coupling behavioral and ecological aspects of reproductive isolation.

Recombination rate variation shapes barriers to introgression across butterfly genomes

Hybridisation and introgression can dramatically alter the relationships among groups of species, leading to phylogenetic discordance across the genome and between populations. Introgression can also erode species differences over time, but selection against introgression at certain loci acts to maintain postmating species barriers. Theory predicts that species barriers made up of many loci throughout the genome should lead to a broad correlation between introgression and recombination rate, which determines the extent to which selection on deleterious foreign alleles will affect neutral alleles at physically linked loci. Here, we describe the variation in genealogical relationships across the genome among three species of Heliconius butterflies: H. melpomene (mel), H. cydno (cyd), and H. timareta (tim), using whole genomes of 92 individuals, and ask whether this variation can be explained by heterogeneous barriers to introgression. We find that species relationships vary predictably at the chromosomal scale. By quantifying recombination rate and admixture proportions, we then show that rates of introgression are predicted by variation in recombination rate. This implies that species barriers are highly polygenic, with selection acting against introgressed alleles across most of the genome. In addition, long chromosomes, which have lower recombination rates, produce stronger barriers on average than short chromosomes. Finally, we find a consistent difference between two species pairs on either side of the Andes, which suggests differences in the architecture of the species barriers. Our findings illustrate how the combined effects of hybridisation, recombination, and natural selection, acting at multitudes of loci over long periods, can dramatically sculpt the phylogenetic relationships among species.

Sexually dimorphic gene expression and transcriptome evolution provide mixed evidence for a fast-Z effect in Heliconius

Sex chromosomes have different evolutionary properties compared to autosomes due to their hemizygous nature. In particular, recessive mutations are more readily exposed to selection, which can lead to faster rates of molecular evolution. Here, we report patterns of gene expression and molecular evolution for a group of butterflies. First, we improve the completeness of the Heliconius melpomene reference annotation, a neotropical butterfly with a ZW sex determination system. Then, we analyse RNA from male and female whole abdomens and sequence female ovary and gut tissue to identify sex‐ and tissue‐specific gene expression profiles in H. melpomene. Using these expression profiles, we compare (a) sequence divergence and polymorphism; (b) the strength of positive and negative selection; and (c) rates of adaptive evolution, for Z and autosomal genes between two species of Heliconius butterflies, H. melpomene and H. erato. We show that the rate of adaptive substitutions is higher for Z than autosomal genes, but contrary to expectation, it is also higher for male‐biased than female‐biased genes. Additionally, we find no significant increase in the rate of adaptive evolution or purifying selection on genes expressed in ovary tissue, a heterogametic‐specific tissue. Our results contribute to a growing body of literature from other ZW systems that also provide mixed evidence for a fast‐Z effect where hemizygosity influences the rate of adaptive substitutions.