Multilocus species trees show the recent adaptive radiation of the mimetic heliconius butterflies

Genome-Wide Gene Birth-Death Dynamics Are Associated with Diet Breadth Variation in Lepidoptera

Comparative analyses of gene birth-death dynamics have the potential to reveal gene families that played an important role in the evolution of morphological, behavioral, or physiological variation. Here, we used whole genomes of 30 species of butterflies and moths to identify gene birth-death dynamics among the Lepidoptera that are associated with specialist or generalist feeding strategies. Our work advances this field using a uniform set of annotated proteins for all genomes, investigating associations while correcting for phylogeny, and assessing all gene families rather than a priori subsets. We discovered that the sizes of several important gene families (e.g. those associated with pesticide resistance, xenobiotic detoxification, and/or protein digestion) are significantly correlated with diet breadth. We also found 22 gene families showing significant shifts in gene birth-death dynamics at the butterfly (Papilionoidea) crown node, the most notable of which was a family of pheromone receptors that underwent a contraction potentially linked with a shift to visual-based mate recognition. Our findings highlight the importance of uniform annotations, phylogenetic corrections, and unbiased gene family analyses in generating a list of candidate genes that warrant further exploration.

Genomic evidence reveals three W-autosome fusions in Heliconius butterflies

Sex chromosomes are evolutionarily labile in many animals and sometimes fuse with autosomes, creating so-called neo-sex chromosomes. Fusions between sex chromosomes and autosomes have been proposed to reduce sexual conflict and to promote adaptation and reproductive isolation among species. Recently, advances in genomics have fuelled the discovery of such fusions across the tree of life. Here, we discovered multiple fusions leading to neo-sex chromosomes in the sapho subclade of the classical adaptive radiation of Heliconius butterflies. Heliconius butterflies generally have 21 chromosomes with very high synteny. However, the five Heliconius species in the sapho subclade show large variation in chromosome number ranging from 21 to 60. We find that the W chromosome is fused with chromosome 4 in all of them. Two sister species pairs show subsequent fusions between the W and chromosomes 9 or 14, respectively. These fusions between autosomes and sex chromosomes make Heliconius butterflies an ideal system for studying the role of neo-sex chromosomes in adaptive radiations and the degeneration of sex chromosomes over time. Our findings emphasize the capability of short-read resequencing to detect genomic signatures of fusion events between sex chromosomes and autosomes even when sex chromosomes are not explicitly assembled.

Hybrid speciation driven by multilocus introgression of ecological traits

Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.

A genetic cause of male mate preference

A gene for mate preference has been shared between hybridizing butterfly species.

Pervasive mimicry in flight behavior among aposematic butterflies

Flight was a key innovation in the adaptive radiation of insects. However, it is a complex trait influenced by a large number of interacting biotic and abiotic factors, making it difficult to unravel the evolutionary drivers. We investigate flight patterns in neotropical heliconiine butterflies, well known for mimicry of their aposematic wing color patterns. We quantify the flight patterns (wing beat frequency and wing angles) of 351 individuals representing 29 heliconiine and 9 ithomiine species belonging to ten color pattern mimicry groupings. For wing beat frequency and up wing angles, we show that heliconiine species group by color pattern mimicry affiliation. Convergence of down wing angles to mimicry groupings is less pronounced, indicating that distinct components of flight are under different selection pressures and constraints. The flight characteristics of the Tiger mimicry group are particularly divergent due to convergence with distantly related ithomiine species. Predator-driven selection for mimicry also explained variation in flight among subspecies, indicating that this convergence can occur over relatively short evolutionary timescales. Our results suggest that the flight convergence is driven by aposematic signaling rather than shared habitat between comimics. We demonstrate that behavioral mimicry can occur between lineages that have separated over evolutionary timescales ranging from <0.5 to 70 My.

Major patterns in the introgression history of Heliconius butterflies

Gene flow between species, although usually deleterious, is an important evolutionary process that can facilitate adaptation and lead to species diversification. It also makes estimation of species relationships difficult. Here, we use the full-likelihood multispecies coalescent (MSC) approach to estimate species phylogeny and major introgression events in Heliconius butterflies from whole-genome sequence data. We obtain a robust estimate of species branching order among major clades in the genus, including the 'melpomene-silvaniform' group, which shows extensive historical and ongoing gene flow. We obtain chromosome-level estimates of key parameters in the species phylogeny, including species divergence times, present-day and ancestral population sizes, as well as the direction, timing, and intensity of gene flow. Our analysis leads to a phylogeny with introgression events that differ from those obtained in previous studies. We find that Heliconius aoede most likely represents the earliest-branching lineage of the genus and that 'silvaniform' species are paraphyletic within the melpomene-silvaniform group. Our phylogeny provides new, parsimonious histories for the origins of key traits in Heliconius, including pollen feeding and an inversion involved in wing pattern mimicry. Our results demonstrate the power and feasibility of the full-likelihood MSC approach for estimating species phylogeny and key population parameters despite extensive gene flow. The methods used here should be useful for analysis of other difficult species groups with high rates of introgression.

Variation of Clasper Scent Gland Composition of Heliconius Butterflies from a Biodiversity Hotspot

Male Heliconius butterflies possess two pheromone emitting structures, wing androconia and abdominal clasper scent glands. The composition of the clasper scent gland of males of 17 Heliconius and Eueides species from an overlapping area in Ecuador, comprising three mimicry groups, was investigated by GC/MS. The chemical signal serves as an anti-aphrodisiac signal that is transferred from males to females during mating, indicating the mating status of the female to prevent them from harassment by other males. In addition, the odour may also serve in predator defence. There is potential for convergence driven by mimicry, although, such convergence might be detrimental for species recognition of the butterflies within the mimicry ring, making mating more difficult. More than 500 compounds were detected, consisting of volatile, semi-volatile or non-volatile compounds, including terpenes, fatty acid esters or aromatic compounds. Several novel esters were identified by GC/MS and GC/IR data, microderivatisation and synthesis, including butyl (Z)-3-dodecenoate and other (Z)-3-alkenoates, 3-oxohexyl citronellate and 5-methylhexa-3,5-dienyl (E)-2,3-dihydrofarnesoate. The secretions were found to be species specific, potentially allowing for species differentiation. Statistical analysis of the compounds showed differentiation by phylogenetic clade and species, but not by mimicry group.

A modified method to analyse cell proliferation using EdU labelling in large insect brains

The study of neurogenesis is critical to understanding of the evolution of nervous systems. Within invertebrates, this process has been extensively studied in Drosophila melanogaster, which is the predominant model thanks to the availability of advanced genetic tools. However, insect nervous systems are extremely diverse, and by studying a range of taxa we can gain additional information about how nervous systems and their development evolve. One example of the high diversity of insect nervous system diversity is provided by the mushroom bodies. Mushroom bodies have critical roles in learning and memory and vary dramatically across species in relative size and the type(s) of sensory information they process. Heliconiini butterflies provide a useful snapshot of this diversity within a closely related clade. Within Heliconiini, the genus Heliconius contains species where mushroom bodies are 3-4 times larger than other closely related genera, relative to the rest of the brain. This variation in size is largely explained by increases in the number of Kenyon cells, the intrinsic neurons which form the mushroom body. Hence, variation in mushroom body size is the product of changes in cell proliferation during Kenyon cell neurogenesis. Studying this variation requires adapting labelling techniques for use in less commonly studied organisms, as methods developed for common laboratory insects often do not work. Here, we present a modified protocol for EdU staining to examine neurogenesis in large-brained insects, using Heliconiini butterflies as our primary case, but also demonstrating applicability to cockroaches, another large-brained insect.

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.

Sex-linked gene traffic underlies the acquisition of sexually dimorphic UV color vision in Heliconius butterflies

The acquisition of novel sexually dimorphic traits poses an evolutionary puzzle: How do new traits arise and become sex-limited? Recently acquired color vision, sexually dimorphic in animals like primates and butterflies, presents a compelling model for understanding how traits become sex-biased. For example, some Heliconius butterflies uniquely possess UV (ultraviolet) color vision, which correlates with the expression of two differentially tuned UV-sensitive rhodopsins, UVRh1 and UVRh2. To discover how such traits become sexually dimorphic, we studied Heliconius charithonia, which exhibits female-specific UVRh1 expression. We demonstrate that females, but not males, discriminate different UV wavelengths. Through whole-genome shotgun sequencing and assembly of the H. charithonia genome, we discovered that UVRh1 is present on the W chromosome, making it obligately female-specific. By knocking out UVRh1, we show that UVRh1 protein expression is absent in mutant female eye tissue, as in wild-type male eyes. A PCR survey of UVRh1 sex-linkage across the genus shows that species with female-specific UVRh1 expression lack UVRh1 gDNA in males. Thus, acquisition of sex linkage is sufficient to achieve female-specific expression of UVRh1, though this does not preclude other mechanisms, like cis-regulatory evolution from also contributing. Moreover, both this event, and mutations leading to differential UV opsin sensitivity, occurred early in the history of Heliconius. These results suggest a path for acquiring sexual dimorphism distinct from existing mechanistic models. We propose a model where gene traffic to heterosomes (the W or the Y) genetically partitions a trait by sex before a phenotype shifts (spectral tuning of UV sensitivity).

Plasticity and genetic effects contribute to different axes of neural divergence in a community of mimetic Heliconius butterflies

Changes in ecological preference, often driven by spatial and temporal variation in resource distribution, can expose populations to environments with divergent information content. This can lead to adaptive changes in the degree to which individuals invest in sensory systems and downstream processes, to optimize behavioural performance in different contexts. At the same time, environmental conditions can produce plastic responses in nervous system development and maturation, providing an alternative route to integrating neural and ecological variation. Here, we explore how these two processes play out across a community of Heliconius butterflies. Heliconius communities exhibit multiple Mullerian mimicry rings, associated with habitat partitioning across environmental gradients. These environmental differences have previously been linked to heritable divergence in brain morphology in parapatric species pairs. They also exhibit a unique dietary adaptation, known as pollen feeding, that relies heavily on learning foraging routes, or trap-lines, between resources, which implies an important environmental influence on behavioural development. By comparing brain morphology across 133 wild-caught and insectary-reared individuals from seven Heliconius species, we find strong evidence for interspecific variation in patterns of neural investment. These largely fall into two distinct patterns of variation; first, we find consistent patterns of divergence in the size of visual brain components across both wild and insectary-reared individuals, suggesting genetically encoded divergence in the visual pathway. Second, we find interspecific differences in mushroom body size, a central component of learning and memory systems, but only among wild caught individuals. The lack of this effect in common-garden individuals suggests an extensive role for developmental plasticity in interspecific variation in the wild. Finally, we illustrate the impact of relatively small-scale spatial effects on mushroom body plasticity by performing experiments altering the cage size and structure experienced by individual H. hecale. Our data provide a comprehensive survey of community level variation in brain structure, and demonstrate that genetic effects and developmental plasticity contribute to different axes of interspecific neural variation.

First chromosome scale genomes of ithomiine butterflies (Nymphalidae: Ithomiini): Comparative models for mimicry genetic studies

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies.

Does sexual conflict contribute to the evolution of novel warning patterns?

Why warning patterns are so diverse is an enduring evolutionary puzzle. Because predators associate particular patterns with unpleasant experiences, an individual's predation risk should decrease as the local density of its warning pattern increases, promoting pattern monomorphism. Distasteful Heliconius butterflies are known for their diversity of warning patterns. Here, we explore whether interlocus sexual conflict can contribute to their diversification. Male Heliconius use warning patterns as mating cues, but mated females may suffer costs if this leads to disturbance, favouring novel patterns. Using simulations, we show that under our model conditions drift alone is unlikely to cause pattern diversification, but that sexual conflict can assist such a process. We also find that genetic architecture influences the evolution of male preferences, which track changes in warning pattern due to sexual selection. When male attraction imposes costs on females, this affects the speed at which novel pattern alleles increase. In two experiments, females laid fewer eggs with males present. However, although males in one experiment showed less interest in females with manipulated patterns, we found no evidence that female colouration mitigates sex-specific costs. Overall, male attraction to conspecific warning patterns may impose an unrecognized cost on Heliconius females, but further work is required to determine this experimentally.

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.

Transposons and non-coding regions drive the intrafamily differences of genome size in insects

Genome size (GS) can vary considerably between phylogenetically close species, but the landscape of GS changes in insects remain largely unclear. To better understand the specific evolutionary factors that determine GS in insects, we examined flow cytometry-based published GS data from 1,326 insect species, spanning 700 genera, 155 families, and 21 orders. Model fitting showed that GS generally followed an Ornstein-Uhlenbeck adaptive evolutionary model in Insecta overall. Ancestral reconstruction indicated a likely GS of 1,069 Mb, suggesting that most insect clades appeared to undergo massive genome expansions or contractions. Quantification of genomic components in 56 species from nine families in four insect orders revealed that the proliferation of transposable elements contributed to high variation in GS between close species, such as within Coleoptera. This study sheds lights on the pattern of GS variation in insects and provides a better understanding of insect GS evolution.

One's trash is someone else's treasure: sequence read archives from Lepidoptera genomes provide material for genome reconstruction of their endosymbionts

BACKGROUND

Maternally inherited bacterial symbionts are extremely widespread in insects. They owe their success to their ability to promote their own transmission through various manipulations of their hosts' life-histories. Many symbionts however very often go undetected. Consequently, we have only a restricted idea of the true symbiont diversity in insects, which may hinder our understanding of even bigger questions in the field such as the evolution or establishment of symbiosis.

RESULTS

In this study, we screened publicly available Lepidoptera genomic material for two of the most common insect endosymbionts, namely Wolbachia and Spiroplasma, in 1904 entries, encompassing 106 distinct species. We compared the performance of two screening software, Kraken2 and MetaPhlAn2, to identify the bacterial infections and using a baiting approach we reconstruct endosymbiont genome assemblies. Of the 106 species screened, 20 (19%) and nine (8.5%) were found to be infected with either Wolbachia or Spiroplasma, respectively. Construction of partial symbiotic genomes and phylogenetic analyses suggested the Wolbachia strains from the supergroup B were the most prevalent type of symbionts, while Spiroplasma infections were scarce in the Lepidoptera species screened here.

CONCLUSIONS

Our results indicate that many of the host-symbiont associations remain largely unexplored, with the majority of associations we identify never being recorded before. This highlights the usefulness of public databases to explore the hidden diversity of symbiotic entities, allowing the development of hypotheses regarding host-symbiont associations. The ever-expanding genomic databases provide a diverse databank from which one can characterize and explore the true diversity of symbiotic entities.

Repeated genetic adaptation to altitude in two tropical butterflies

Repeated evolution can provide insight into the mechanisms that facilitate adaptation to novel or changing environments. Here we study adaptation to altitude in two tropical butterflies, Heliconius erato and H. melpomene, which have repeatedly and independently adapted to montane habitats on either side of the Andes. We sequenced 518 whole genomes from altitudinal transects and found many regions differentiated between highland (~ 1200 m) and lowland (~ 200 m) populations. We show repeated genetic differentiation across replicate populations within species, including allopatric comparisons. In contrast, there is little molecular parallelism between the two species. By sampling five close relatives, we find that a large proportion of divergent regions identified within species have arisen from standing variation and putative adaptive introgression from high-altitude specialist species. Taken together our study supports a role for both standing genetic variation and gene flow from independently adapted species in promoting parallel local adaptation to the environment.

Here, the authors study adaptation to altitude in 518 whole genomes from two species of tropical butterflies. They find repeated genetic differentiation within species, little molecular parallelism between these species, and introgression from closely related species, concluding that standing genetic variation promotes parallel local adaptation.

Validation of reference-assisted assembly using existing and novel Heliothine genomes

Chloridea subflexa and Chloridea virescens are a pair of closely related noctuid species exhibiting pheromone-based sexual isolation and divergent host plant preferences. We produced a novel Illumina short read C. subflexa genome assembly and an improved C. virescens genome assembly, which offer opportunities to study the genomic basis for evolutionarily important traits in this lepidopteran family with few genomic resources. We then examined the feasibility of reference-assisted assembly, an approach that leverages existing high quality genomic resources for genome improvement in closely related taxa and applied it to our Heliothine genomes. Our work demonstrates that reference-assisted assembly has the potential to enhance contiguity and completeness of existing insect genomic resources with minimal additional laboratory costs. We conclude by discussing both the potential and pitfalls of reference-assisted assembly according to the intended downstream assembly application.

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.

A peculiar new species of Dione (Agraulis) Boisduval & Le Conte (Lepidoptera, Nymphalidae, Heliconiinae) associated with Malesherbia Ruiz & Pavón (Passifloraceae) in xeric western slopes of the Andes

Butterflies associated with xerophytic environments of the Andes have been little studied, and they exhibit high levels of endemism. Herein Dione (Agraulis) dodona Lamas & Farfán, sp. nov. (Nymphalidae; Heliconiinae) is described, distributed on the western slopes of the Andes of Peru and northern Chile, between 800 and 3,000 m elevation. Adults of both sexes, and the immature stages, are described and illustrated based on light and scanning electron microscopy. The immature stages are associated with MalesherbiatenuifoliaD. Don (Passifloraceae) found in xeric environments, representing a new record of this genus as a host plant for the subfamily Heliconiinae. Conspicuous morphological differences are presented for all stages at the generic level. Based on a phylogenetic analysis of the COI barcode mitochondrial gene fragment, D. (A.) dodona Lamas & Farfán, sp. nov. is distinguished as an independent lineage within the Agraulis clade of Dione, with ca. 5% difference to congeneric species.

The evolution of adult pollen feeding did not alter postembryonic growth in Heliconius butterflies

For many animals, the availability and provision of dietary resources can vary markedly between juvenile and adult stages, often leading to a temporal separation of nutrient acquisition and use. Juvenile developmental programs are likely limited by the energetic demands of many adult tissues and processes with early developmental origins. Enhanced dietary quality in the adult stage may, therefore, alter selection on life history and growth patterns in juvenile stages. Heliconius are unique among butterflies in actively collecting and digesting pollen grains, which provide an adult source of essential amino acids. The origin of pollen feeding has therefore previously been hypothesized to lift constraints on larval growth rates, allowing Heliconius to spend less time as larvae when they are most vulnerable to predation. By measuring larval and pupal life-history traits across three pollen-feeding and three nonpollen-feeding Heliconiini, we provide the first test of this hypothesis. Although we detect significant interspecific variation in larval and pupal development, we do not find any consistent shift associated with pollen feeding. We discuss how this result may fit with patterns of nitrogen allocation, the benefits of nitrogenous stores, and developmental limitations on growth. Our results provide a framework for studies aiming to link innovations in adult Heliconius to altered selection regimes and developmental programs in early life stages.

The genetic basis of structural colour variation in mimetic Heliconius butterflies

Structural colours, produced by the reflection of light from ultrastructures, have evolved multiple times in butterflies. Unlike pigmentary colours and patterns, little is known about the genetic basis of these colours. Reflective structures on wing-scale ridges are responsible for iridescent structural colour in many butterflies, including the Müllerian mimics Heliconius erato and Heliconius melpomene. Here, we quantify aspects of scale ultrastructure variation and colour in crosses between iridescent and non-iridescent subspecies of both of these species and perform quantitative trait locus (QTL) mapping. We show that iridescent structural colour has a complex genetic basis in both species, with offspring from crosses having a wide variation in blue colour (both hue and brightness) and scale structure measurements. We detect two different genomic regions in each species that explain modest amounts of this variation, with a sex-linked QTL in H. erato but not H. melpomene. We also find differences between species in the relationships between structure and colour, overall suggesting that these species have followed different evolutionary trajectories in their evolution of structural colour. We then identify genes within the QTL intervals that are differentially expressed between subspecies and/or wing regions, revealing likely candidates for genes controlling structural colour formation.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Hybrid enrichment of adaptive variation revealed by genotype-environment associations in montane sedges

The role of hybridization in diversification is complex and may result in many possible outcomes. Not only can hybridization produce new lineages, but those lineages may contain unique combinations of adaptive genetic variation derived from parental taxa that allow hybrid‐origin lineages to occupy unique environmental space relative to one (or both) parent(s). We document such a case of hybridization between two sedge species, Carex nova and Carex nelsonii (Cyperaceae), that occupy partially overlapping environmental space in the southern Rocky Mountains, USA. In the region hypothesized to be the origin of the hybrid lineage, one parental taxon (C. nelsonii) is at the edge of its environmental tolerance. Hybrid‐origin individuals display mixed ancestry between the parental taxa—of nearly 7000 unlinked loci sampled, almost 30% showed evidence of excess ancestry from one parental lineage—approximately half displayed a genomic background skewed towards one parent, and half skewed towards the other. To test whether excess ancestry loci may have conferred an adaptive advantage to the hybrid‐origin lineage, we conducted genotype–environment association analyses on different combinations of loci—with and without excess ancestry—and with multiple contrasts between the hybrids and parental taxa. Loci with skewed ancestry showed significant environmental associations distinguishing the hybrid lineage from one parent (C. nelsonii), whereas loci with relatively equal representation of parental ancestries showed no such environmental associations. Moreover, the overwhelming majority of candidate adaptive loci with respect to environmental gradients also had excess ancestry from a parental lineage, implying these loci have facilitated the persistence of the hybrid lineage in an environment unsuitable to at least one parent.

Multiple mechanisms of photoreceptor spectral tuning in Heliconius butterflies

The evolution of color vision is often studied through the lens of receptor gain relative to an ancestor with fewer spectral classes of photoreceptor. For instance, in Heliconius butterflies, a genus-specific UVRh opsin duplication led to the evolution of UV color discrimination in Heliconius erato females, a rare trait among butterflies. However, color vision evolution is not well understood in the context of loss. In Heliconius melpomene and Heliconius ismenius lineages, the UV2 receptor subtype has been lost, which limits female color vision in shorter wavelengths. Here, we compare the visual systems of butterflies that have either retained or lost the UV2 photoreceptor using intracellular recordings, ATAC-seq, and antibody staining. We identify several ways these butterflies modulate their color vision. In H. melpomene, chromatin reorganization has downregulated an otherwise intact UVRh2 gene, whereas in H. ismenius, pseudogenization has led to the truncation of UVRh2. In species that lack the UV2 receptor, the peak sensitivity of the remaining UV1 photoreceptor cell is shifted to longer wavelengths. Across Heliconius, we identify the widespread use of filtering pigments and co-expression of two opsins in the same photoreceptor cells. Multiple mechanisms of spectral tuning, including the molecular evolution of blue opsins, have led to the divergence of receptor sensitivities between species. The diversity of photoreceptor and ommatidial subtypes between species suggests that Heliconius visual systems are under varying selection pressures for color discrimination. Modulating the wavelengths of peak sensitivities of both the blue- and remaining UV-sensitive photoreceptor cells suggests that Heliconius species may have compensated for UV receptor loss.

OUP accepted manuscript

Introgressive hybridization plays a key role in adaptive evolution and species diversification in many groups of species. However, frequent hybridization and gene flow between species make estimation of the species phylogeny and key population parameters challenging. Here, we show that by accounting for phasing and using full-likelihood methods, introgression histories and population parameters can be estimated reliably from whole-genome sequence data. We employ the multispecies coalescent (MSC) model with and without gene flow to infer the species phylogeny and cross-species introgression events using genomic data from six members of the erato-sara clade of Heliconius butterflies. The methods naturally accommodate random fluctuations in genealogical history across the genome due to deep coalescence. To avoid heterozygote phasing errors in haploid sequences commonly produced by genome assembly methods, we process and compile unphased diploid sequence alignments and use analytical methods to average over uncertainties in heterozygote phase resolution. There is robust evidence for introgression across the genome, both among distantly related species deep in the phylogeny and between sister species in shallow parts of the tree. We obtain chromosome-specific estimates of key population parameters such as introgression directions, times and probabilities, as well as species divergence times and population sizes for modern and ancestral species. We confirm ancestral gene flow between the sara clade and an ancestral population of Heliconius telesiphe, a likely hybrid speciation origin for Heliconius hecalesia, and gene flow between the sister species Heliconius erato and Heliconius himera. Inferred introgression among ancestral species also explains the history of two chromosomal inversions deep in the phylogeny of the group. This study illustrates how a full-likelihood approach based on the MSC makes it possible to extract rich historical information of species divergence and gene flow from genomic data. [3s; bpp; gene flow; Heliconius; hybrid speciation; introgression; inversion; multispecies coalescent]

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.

An accurate assignment test for extremely low-coverage whole-genome sequence data

Genomic assignment tests can provide important diagnostic biological characteristics, such as population of origin or ecotype. Yet, assignment tests often rely on moderate- to high-coverage sequence data that can be difficult to obtain for fields such as molecular ecology and ancient DNA. We have developed a novel approach that efficiently assigns biologically relevant information (i.e., population identity or structural variants such as inversions) in extremely low-coverage sequence data. First, we generate databases from existing reference data using a subset of diagnostic single nucleotide polymorphisms (SNPs) associated with a biological characteristic. Low-coverage alignment files are subsequently compared to these databases to ascertain allelic state, yielding a joint probability for each association. To assess the efficacy of this approach, we assigned haplotypes and population identity in Heliconius butterflies, Atlantic herring, and Atlantic cod using chromosomal inversion sites and whole-genome data. We scored both modern and ancient specimens, including the first whole-genome sequence data recovered from ancient Atlantic herring bones. The method accurately assigns biological characteristics, including population membership, using extremely low-coverage data (as low as 0.0001x) based on genome-wide SNPs. This approach will therefore increase the number of samples in evolutionary, ecological and archaeological research for which relevant biological information can be obtained.

Genetic Causes and Consequences of Sympatric Morph Divergence in Salmonidae: A Search for Mechanisms

Repeatedly and recently evolved sympatric morphs exhibiting consistent phenotypic differences provide natural experimental replicates of speciation. Because such morphs are observed frequently in Salmonidae, this clade provides a rare opportunity to uncover the genomic mechanisms underpinning speciation. Such insight is also critical for conserving salmonid diversity, the loss of which could have significant ecological and economic consequences. Our review suggests that genetic differentiation among sympatric morphs is largely nonparallel apart from a few key genes that may be critical for consistently driving morph differentiation. We discuss alternative levels of parallelism likely underlying consistent morph differentiation and identify several factors that may temper this incipient speciation between sympatric morphs, including glacial history and contemporary selective pressures. Our synthesis demonstrates that salmonids are useful for studying speciation and poses additional research questions to be answered by future study of this family. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Environmental Drivers of Diversification and Hybridization in Neotropical Butterflies

Studying how the environment shapes current biodiversity patterns in species rich regions is a fundamental issue in biogeography, ecology, and conservation. However, in the Neotropics, the study of the forces driving species distribution and richness, is mostly based on vertebrates and plants. In this study, we used 54,392 georeferenced records for 46 species and 1,012 georeferenced records for 38 interspecific hybrids of the Neotropical Heliconius butterflies to investigate the role of the environment in shaping their distribution and richness, as well as their geographic patterns of phylogenetic diversity and phylogenetic endemism. We also evaluated whether niche similarity promotes hybridization in Heliconius. We found that these insects display five general distribution patterns mostly explained by precipitation and isothermality, and to a lesser extent, by altitude. Interestingly, altitude plays a major role as a predictor of species richness and phylogenetic diversity, while precipitation explains patterns of phylogenetic endemism. We did not find evidence supporting the role of the environment in facilitating hybridization because hybridizing species do not necessarily share the same climatic niche despite some of them having largely overlapping geographic distributions. Overall, we confirmed that, as in other organisms, high annual temperature, a constant supply of water, and spatio-topographic complexity are the main predictors of diversity in Heliconius. However, future studies at large scale need to investigate the effect of microclimate variables and ecological interactions.

Conserved ancestral tropical niche but different continental histories explain the latitudinal diversity gradient in brush-footed butterflies

The global increase in species richness toward the tropics across continents and taxonomic groups, referred to as the latitudinal diversity gradient, stimulated the formulation of many hypotheses to explain the underlying mechanisms of this pattern. We evaluate several of these hypotheses to explain spatial diversity patterns in a butterfly family, the Nymphalidae, by assessing the contributions of speciation, extinction, and dispersal, and also the extent to which these processes differ among regions at the same latitude. We generate a time-calibrated phylogeny containing 2,866 nymphalid species (~45% of extant diversity). Neither speciation nor extinction rate variations consistently explain the latitudinal diversity gradient among regions because temporal diversification dynamics differ greatly across longitude. The Neotropical diversity results from low extinction rates, not high speciation rates, and biotic interchanges with other regions are rare. Southeast Asia is also characterized by a low speciation rate but, unlike the Neotropics, is the main source of dispersal events through time. Our results suggest that global climate change throughout the Cenozoic, combined with tropical niche conservatism, played a major role in generating the modern latitudinal diversity gradient of nymphalid butterflies.

Comparative transcriptome analysis at the onset of speciation in a mimetic butterfly-The Ithomiini Melinaea marsaeus

Ecological speciation entails divergent selection on specific traits and ultimately on the developmental pathways responsible for these traits. Selection can act on gene sequences but also on regulatory regions responsible for gene expression. Mimetic butterflies are a relevant system for speciation studies because wing colour pattern (WCP) often diverges between closely related taxa and is thought to drive speciation through assortative mating and increased predation on hybrids. Here, we generate the first transcriptomic resources for a mimetic butterfly of the tribe Ithomiini, Melinaea marsaeus, to examine patterns of differential expression between two subspecies and between tissues that express traits that likely drive reproductive isolation; WCP and chemosensory genes. We sequenced whole transcriptomes of three life stages to cover a large catalogue of transcripts, and we investigated differential expression between subspecies in pupal wing discs and antennae. Eighteen known WCP genes were expressed in wing discs and 115 chemosensory genes were expressed in antennae, with a remarkable diversity of chemosensory protein genes. Many transcripts were differentially expressed between subspecies, including two WCP genes and one odorant receptor. Our results suggest that in M. marsaeus the same genes as in other mimetic butterflies are involved in traits causing reproductive isolation, and point at possible candidates for the differences in those traits between subspecies. Differential expression analyses of other developmental stages and body organs and functional studies are needed to confirm and expand these results. Our work provides key resources for comparative genomics in mimetic butterflies, and more generally in Lepidoptera.

Chromosome Fusion Affects Genetic Diversity and Evolutionary Turnover of Functional Loci but Consistently Depends on Chromosome Size

Major changes in chromosome number and structure are linked to a series of evolutionary phenomena, including intrinsic barriers to gene flow or suppression of recombination due to chromosomal rearrangements. However, chromosome rearrangements can also affect the fundamental dynamics of molecular evolution within populations by changing relationships between linked loci and altering rates of recombination. Here, we build chromosome-level assembly Eueides isabella and, together with a recent chromosome-level assembly of Dryas iulia, examine the evolutionary consequences of multiple chromosome fusions in Heliconius butterflies. These assemblies pinpoint fusion points on 10 of the 20 autosomal chromosomes and reveal striking differences in the characteristics of fused and unfused chromosomes. The ten smallest autosomes in D. iulia and E. isabella, which have each fused to a longer chromosome in Heliconius, have higher repeat and GC content, and longer introns than predicted by their chromosome length. When fused, these characteristics change to become more in line with chromosome length. The fusions also led to reduced diversity, which likely reflects increased background selection and selection against introgression between diverging populations, following a reduction in per-base recombination rate. We further show that chromosome size and fusion impact turnover rates of functional loci at a macroevolutionary scale. Together these results provide further evidence that chromosome fusion in Heliconius likely had dramatic effects on population level processes shaping rates of neutral and adaptive divergence. These effects may have impacted patterns of diversification in Heliconius, a classic example of an adaptive radiation.

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.

The Dryas iulia Genome Supports Multiple Gains of a W Chromosome from a B Chromosome in Butterflies

In butterflies and moths, which exhibit highly variable sex determination mechanisms, the homogametic Z chromosome is deeply conserved and is featured in many genome assemblies. The evolution and origin of the female W sex chromosome, however, remains mostly unknown. Previous studies have proposed that a ZZ/Z0 sex determination system is ancestral to Lepidoptera, and that W chromosomes may originate from sex-linked B chromosomes. Here, we sequence and assemble the female Dryas iulia genome into 32 highly contiguous ordered and oriented chromosomes, including the Z and W sex chromosomes. We then use sex-specific Hi-C, ATAC-seq, PRO-seq, and whole-genome DNA sequence data sets to test if features of the D. iulia W chromosome are consistent with a hypothesized B chromosome origin. We show that the putative W chromosome displays female-associated DNA sequence, gene expression, and chromatin accessibility to confirm the sex-linked function of the W sequence. In contrast with expectations from studies of homologous sex chromosomes, highly repetitive DNA content on the W chromosome, the sole presence of domesticated repetitive elements in functional DNA, and lack of sequence homology with the Z chromosome or autosomes is most consistent with a B chromosome origin for the W, although it remains challenging to rule out extensive sequence divergence. Synteny analysis of the D. iulia W chromosome with other female lepidopteran genome assemblies shows no homology between W chromosomes and suggests multiple, independent origins of the W chromosome from a B chromosome likely occurred in butterflies.

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.

CARACTERIZACIÓN GENÉTICA DE LA POBLACIÓN DE Heliconius sara (Nymphalidae) EN LA ISLA GORGONA, COLOMBIA

La estructura genética de poblaciones de mariposas con distribución en islas y sus pares continentales ha sido poco documentada para el neotrópico. Este estudio presenta la caracterización de una población de Heliconius sara con distribución en la Isla Gorgona, ubicada en la región del Pacífico Oriental Colombiano. Para esto se examinaron secuencias parciales de un marcador mitocondrial incluyendo información obtenida del GenBank. Se comparó la diversidad y estructura genética con sus conespecíficos continentales y también con congéneres, con los que comparte un ancestro común cercano en el clado Sapho-Sara. Para el análisis de diversidad y estructura genética se realizó un análisis molecular de varianza. Este análisis muestra que la distancia entre la población de la isla y sus pares en el continente es consistente con la variación intraespecífica observada en otras especies del género Heliconius. Para la reconstrucción de la genealogía y datación reciente en el Pleistoceno superior del grupo monofilético de secuencias de H. sara, se realizó un análisis de inferencia bayesiana, así como una de máxima verosimilitud. Del análisis demográfico se seleccionó un modelo histórico de flujo asimétrico desde la isla hacia el continente que sugiere baja resistencia de la discontinuidad geográfica a la dispersión de esta mariposa diurna desde la isla. Este es el primer estudio en examinar un posible evento de aislamiento de una población insular de mariposas en Colombia.

Synteny-Based Genome Assembly for 16 Species of Heliconius Butterflies, and an Assessment of Structural Variation across the Genus

Heliconius butterflies (Lepidoptera: Nymphalidae) are a group of 48 neotropical species widely studied in evolutionary research. Despite the wealth of genomic data generated in past years, chromosomal level genome assemblies currently exist for only two species, Heliconius melpomene and Heliconius erato, each a representative of one of the two major clades of the genus. Here, we use these reference genomes to improve the contiguity of previously published draft genome assemblies of 16 Heliconius species. Using a reference-assisted scaffolding approach, we place and order the scaffolds of these genomes onto chromosomes, resulting in 95.7-99.9% of their genomes anchored to chromosomes. Genome sizes are somewhat variable among species (270-422 Mb) and in one small group of species (Heliconius hecale, Heliconius elevatus, and Heliconius pardalinus) expansions in genome size are driven mainly by repetitive sequences that map to four small regions in the H. melpomene reference genome. Genes from these repeat regions show an increase in exon copy number, an absence of internal stop codons, evidence of constraint on nonsynonymous changes, and increased expression, all of which suggest that at least some of the extra copies are functional. Finally, we conducted a systematic search for inversions and identified five moderately large inversions fixed between the two major Heliconius clades. We infer that one of these inversions was transferred by introgression between the lineages leading to the erato/sara and burneyi/doris clades. These reference-guided assemblies represent a major improvement in Heliconius genomic resources that enable further genetic and evolutionary discoveries in this genus.

Habitat generalist species constrain the diversity of mimicry rings in heterogeneous habitats

How evolution creates and maintains trait patterns in species-rich communities is still an unsolved topic in evolutionary ecology. One classical example of community-level pattern is the unexpected coexistence of different mimicry rings, each of which is a group of mimetic species with the same warning signal. The coexistence of different mimicry rings in a community seems paradoxical because selection among unpalatable species should favor convergence to a single warning pattern. We combined mathematical modeling based on network theory and numerical simulations to explore how different types of selection, such as mimetic and environmental selections, and habitat use by mimetic species influence the formation of coexisting rings. We show that when habitat and mimicry are strong sources of selection, the formation of multiple rings takes longer due to conflicting selective pressures. Moreover, habitat generalist species decrease the distinctiveness of different mimicry rings’ patterns and a few habitat generalist species can generate a “small-world effect”, preventing the formation of multiple mimicry rings. These results may explain why the coexistence of mimicry rings is more common in groups of animals that tend towards habitat specialism, such as butterflies.

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.

Convergence in sympatry: Evolution of blue-banded wing pattern in Morpho butterflies

Species interactions such as mimicry can promote trait convergence but disentangling this effect from those of shared ecology, evolutionary history, and niche conservatism is often challenging. Here by focusing on wing colour pattern variation within and between three butterfly species living in sympatry in a large proportion of their range, we tested the effect of species interactions on trait diversification. These butterflies display a conspicuous iridescent blue coloration on the dorsal side of their wings and a cryptic brownish colour on the ventral side. Combined with an erratic and fast flight, these colour patterns increase the difficulty of capture by predators and contribute to the high escape abilities of these butterflies. We hypothesize that, beyond their direct contribution to predator escape, these wing patterns can be used as signals of escape abilities by predators, resulting in positive frequency-dependent selection favouring convergence in wing pattern in sympatry. To test this hypothesis, we quantified dorsal wing pattern variations of 723 butterflies from the three species sampled throughout their distribution, including sympatric and allopatric situations and compared the phenotypic distances between species, sex and localities. We detected a significant effect of localities on colour pattern, and higher inter-specific resemblance in sympatry as compared to allopatry, consistent with the hypothesis of local convergence of wing patterns. Our results provide support to the existence of escape mimicry in the wild and stress the importance of estimating trait variation within species to understand trait variation between species, and to a larger extent, trait diversification at the macro-evolutionary scale.

Heliconius Butterflies Host Characteristic and Phylogenetically Structured Adult-Stage Microbiomes

Lepidoptera (butterflies and moths) are diverse and ecologically important, yet we know little about how they interact with microbes as adults. Due to metamorphosis, the form and function of their adult-stage microbiomes might be very different from those of microbiomes in the larval stage (caterpillars). We studied adult-stage microbiomes of Heliconius and closely related passion-vine butterflies (Heliconiini), which are an important model system in evolutionary biology. To characterize the structure and dynamics of heliconiine microbiomes, we used field collections of wild butterflies, 16S rRNA gene sequencing, quantitative PCR, and shotgun metagenomics. We found that Heliconius butterflies harbor simple and abundant bacterial communities that are moderately consistent among conspecific individuals and over time. Heliconiine microbiomes also exhibited a strong signal of the host phylogeny, with a major distinction between Heliconius and other butterflies. These patterns were largely driven by differing relative abundances of bacterial phylotypes shared among host species and genera, as opposed to the presence or absence of host-specific phylotypes. We suggest that the phylogenetic structure in heliconiine microbiomes arises from conserved host traits that differentially filter microbes from the environment. While the relative importance of different traits remains unclear, our data indicate that pollen feeding (unique to Heliconius) is not a primary driver. Using shotgun metagenomics, we also discovered trypanosomatids and microsporidia to be prevalent in butterfly guts, raising the possibility of antagonistic interactions between eukaryotic parasites and colocalized gut bacteria. Our discovery of characteristic and phylogenetically structured microbiomes provides a foundation for tests of adult-stage microbiome function, a poorly understood aspect of lepidopteran biology.IMPORTANCE Many insects host microbiomes with important ecological functions. However, the prevalence of this phenomenon is unclear because in many insect taxa, microbiomes have been studied in only part of the life cycle, if at all. A prominent example is butterflies and moths, in which the composition and functional role of adult-stage microbiomes are largely unknown. We comprehensively characterized microbiomes in adult passion-vine butterflies. Butterfly-associated bacterial communities are generally abundant in guts, consistent within populations, and composed of taxa widely shared among hosts. More closely related butterflies harbor more similar microbiomes, with the most dramatic shift in microbiome composition occurring in tandem with a suite of ecological and life history traits unique to the genus Heliconius Butterflies are also frequently infected with previously undescribed eukaryotic parasites, which may interact with bacteria in important ways. These findings advance our understanding of butterfly biology and insect-microbe interactions generally.

Pollen feeding in Heliconius butterflies: the singular evolution of an adaptive suite

Major evolutionary transitions can be triggered by behavioural novelty, and are often associated with 'adaptive suites', which involve shifts in multiple co-adapted traits subject to complex interactions. Heliconius butterflies represent one such example, actively feeding on pollen, a behaviour unique among butterflies. Pollen feeding permits a prolonged reproductive lifespan, and co-occurs with a constellation of behavioural, neuroanatomical, life history, morphological and physiological traits that are absent in closely related, non-pollen-feeding genera. As a highly tractable system, supported by considerable ecological and genomic data, Heliconius are an excellent model for investigating how behavioural innovation can trigger a cascade of adaptive shifts in multiple diverse, but interrelated, traits. Here, we synthesize current knowledge of pollen feeding in Heliconius, and explore potential interactions between associated, putatively adaptive, traits. Currently, no physiological, morphological or molecular innovation has been explicitly linked to the origin of pollen feeding, and several hypothesized links between different aspects of Heliconius biology remain poorly tested. However, resolving these uncertainties will contribute to our understanding of how behavioural innovations evolve and subsequently alter the evolutionary trajectories of diverse traits impacting resource acquisition, life history, senescence and cognition.

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.

Phylomitogenomics provides new perspectives on the Euphasmatodea radiation (Insecta: Phasmatodea)

Phasmatodea species diversity lies almost entirely within its suborder Euphasmatodea, which exhibits a pantropical distribution and is considered to derive from a recent and rapid evolutionary radiation. To shed light on Euphasmatodea origins and diversification, we assembled the mitogenomes of 17 species from transcriptomic sequencing data and analysed them along with 22 already available Phasmatodea mitogenomes and 33 mitogenomes representing most of the Polyneoptera lineages. Maximum Likelihood and Bayesian Inference approaches retrieved consistent topologies, both showing the widespread conflict between phylogenetic approaches and traditional systematics. We performed a divergence time analysis leveraging ten fossil specimens representative of most polyneopteran lineages: the time tree obtained supports an older radiation of the clade with respect to previous hypotheses. Euphasmatodea diversification is inferred to have started ~ 187 million years ago, suggesting that the Triassic-Jurassic mass extinction and the breakup of Pangea could have contributed to the process. We also investigated Euphasmatodea mitogenomes patterns of dN, dS and dN/dS ratio throughout our time-tree, trying to characterize the selective regime which may have shaped the clade evolution.

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.

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.