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.

Demographically explicit scans for barriers to gene flow using gIMble.

Identifying regions of the genome that act as barriers to gene flow between recently diverged taxa has remained challenging given the many evolutionary forces that generate variation in genetic diversity and divergence along the genome, and the stochastic nature of this variation. Progress has been impeded by a conceptual and methodological divide between analyses that infer the demographic history of speciation and genome scans aimed at identifying locally maladaptive alleles i.e. genomic barriers to gene flow. Here we implement genomewide IM blockwise likelihood estimation (gIMble), a composite likelihood approach for the quantification of barriers, that bridges this divide. This analytic framework captures background selection and selection against barriers in a model of isolation with migration (IM) as heterogeneity in effective population size (Ne) and effective migration rate (me), respectively. Variation in both effective demographic parameters is estimated in sliding windows via pre-computed likelihood grids. gIMble includes modules for pre-processing/filtering of genomic data and performing parametric bootstraps using coalescent simulations. To demonstrate the new approach, we analyse data from a well-studied pair of sister species of tropical butterflies with a known history of post-divergence gene flow: Heliconius melpomene and H. cydno. Our analyses uncover both large-effect barrier loci (including well-known wing-pattern genes) and a genome-wide signal of a polygenic barrier architecture.

Do chromosome rearrangements fix by genetic drift or natural selection? Insights from Brenthis butterflies.

Large-scale chromosome rearrangements, such as fissions and fusions, are a common feature of eukaryote evolution. They can have considerable influence on the evolution of populations, yet it remains unclear exactly how rearrangements become established and eventually fix. Rearrangements could fix by genetic drift if they are weakly deleterious or neutral, or they may instead be favoured by positive natural selection. Here, we compare genome assemblies of three closely related Brenthis butterfly species and characterize a complex history of fission and fusion rearrangements. An inferred demographic history of these species suggests that rearrangements became fixed in populations with large long-term effective size (Ne ), consistent with rearrangements being selectively neutral or only very weakly underdominant. Using a recently developed analytic framework for characterizing hard selective sweeps, we find that chromosome fusions are not enriched for evidence of past sweeps compared to other regions of the genome. Nonetheless, we do infer a strong and recent selective sweep around one chromosome fusion in the B. daphne genome. Our results suggest that rearrangements in these species likely have weak absolute fitness effects and fix by genetic drift. However, one putative selective sweep raises the possibility that natural selection may sometimes play a role in the fixation of chromosome fusions.

Morphological changes in female reproductive organs in the African monarch butterfly, host to a male-killing Spiroplasma.

Background: Sexual selection and conflicts within and between sexes promote morphological diversity of reproductive traits within species. Variation in the morphology of diagnostic reproductive characters within species offer an excellent opportunity to study these evolutionary processes as drivers of species diversification. The African monarch, Danaus chrysippus (Linnaeus, 1758), is widespread across Africa. The species is polytypic, with the respective geographical ranges of the four colour morphs only overlapping in East Africa. Furthermore, some of the populations host an endosymbiotic bacterium, Spiroplasma, which induces son-killing and distorts the local host population sex-ratio, creating sexual conflicts between the females seeking to optimize their fecundity and the limited mating capacity of the rare males. Methods: We dissected females from Kenya, Rwanda and South Africa, where Spiroplasma vary in presence and prevalence (high, variable and absent, respectively), and conducted microscopy imaging of their reproductive organs. We then characterized the effect of population, female body size, and female mating status, on the size and shape of different genitalia characters of the D. chrysippus female butterflies. Results: We showed that although the general morphology of the organs is conserved in D. chrysippus, female genitalia vary in size and shape between and within populations. The virgin females have smaller organs, while the same organs were expanded in mated females. Females from highly female-biased populations, where the male-killing Spiroplasma is prevalent, also have a larger area of their corpus bursae covered with signa structures. However, this pattern occurs because a larger proportion of the females remains virgin in the female-biased populations rather than because of male depletion due to the symbiont, as males from sex-ratio distorted populations did not produce significantly smaller nutritious spermatophores.

Space-for-time substitutions in climate change ecology and evolution.

In an epoch of rapid environmental change, understanding and predicting how biodiversity will respond to a changing climate is an urgent challenge. Since we seldom have sufficient long-term biological data to use the past to anticipate the future, spatial climate-biotic relationships are often used as a proxy for predicting biotic responses to climate change over time. These 'space-for-time substitutions' (SFTS) have become near ubiquitous in global change biology, but with different subfields largely developing methods in isolation. We review how climate-focussed SFTS are used in four subfields of ecology and evolution, each focussed on a different type of biotic variable - population phenotypes, population genotypes, species' distributions, and ecological communities. We then examine the similarities and differences between subfields in terms of methods, limitations and opportunities. While SFTS are used for a wide range of applications, two main approaches are applied across the four subfields: spatial in situ gradient methods and transplant experiments. We find that SFTS methods share common limitations relating to (i) the causality of identified spatial climate-biotic relationships and (ii) the transferability of these relationships, i.e. whether climate-biotic relationships observed over space are equivalent to those occurring over time. Moreover, despite widespread application of SFTS in climate change research, key assumptions remain largely untested. We highlight opportunities to enhance the robustness of SFTS by addressing key assumptions and limitations, with a particular emphasis on where approaches could be shared between the four subfields.

Subtle Introgression Footprints at the End of the Speciation Continuum in a Clade of Heliconius Butterflies.

Quantifying gene flow between lineages at different stages of the speciation continuum is central to understanding speciation. Heliconius butterflies have undergone an adaptive radiation in wing color patterns driven partly by natural selection for local mimicry. Color patterns are also known to be used as assortative mating cues. Therefore, wing pattern divergence is considered to play a role in speciation. A corollary is that mimicry between closely related species may be associated with hybridization and interfere with reproductive isolation. Here, we take a multifaceted approach to explore speciation history, species boundaries, and traits involved in species differentiation between the two closely related species, Heliconius hecale and Heliconius ismenius. We focus on geographic regions where the two species mimic each other and contrast this with geographic regions where they do not mimic each other. To examine population history and patterns of gene flow, we tested and compared a four-population model accounting for linked selection. This model suggests that the two species have remained isolated for a large part of their history, yet with a small amount of gene exchange. Accordingly, signatures of genomic introgression were small except at a major wing pattern allele and chemosensing genes and stronger in the mimetic populations compared with nonmimetic populations. Behavioral assays confirm that visual confusion exists but that short-range cues determine strong sexual isolation. Tests for chemical differentiation between species identified major differences in putative pheromones which likely mediate mate choice and the maintenance of species differences.

Recent Evolution of a Maternally Acting Sex-Determining Supergene in a Fly with Single-Sex Broods.

Sex determination is a key developmental process, yet it is remarkably variable across the tree of life. The dipteran family Sciaridae exhibits one of the most unusual sex determination systems in which mothers control offspring sex through selective elimination of paternal X chromosomes. Whereas in some members of the family females produce mixed-sex broods, others such as the dark-winged fungus gnat Bradysia coprophila are monogenic, with females producing single-sex broods. Female-producing females were previously found to be heterozygous for a large X-linked paracentric inversion (X'), which is maternally inherited and absent from male-producing females. Here, we assembled and characterized the X' sequence. As close sequence homology between the X and X' made identification of the inversion challenging, we developed a k-mer-based approach to bin genomic reads before assembly. We confirmed that the inversion spans most of the X' chromosome (∼55 Mb) and encodes ∼3,500 genes. Analysis of the divergence between the inversion and the homologous region of the X revealed that it originated very recently (<0.5 Ma). Surprisingly, we found that the X' is more complex than previously thought and is likely to have undergone multiple rearrangements that have produced regions of varying ages, resembling a supergene composed of evolutionary strata. We found functional degradation of ∼7.3% of genes within the region of recombination suppression, but no evidence of accumulation of repetitive elements. Our findings provide an indication that sex-linked inversions are driving turnover of the strange sex determination system in this family of flies.

Genome sequence of Ophryocystis elektroscirrha, an apicomplexan parasite of monarch butterflies: cryptic diversity and response to host-sequestered plant chemicals.

Apicomplexa are ancient and diverse organisms which have been poorly characterized by modern genomics. To better understand the evolution and diversity of these single-celled eukaryotes, we sequenced the genome of Ophryocystis elektroscirrha, a parasite of monarch butterflies, Danaus plexippus. We contextualize our newly generated resources within apicomplexan genomics before answering longstanding questions specific to this host-parasite system. To start, the genome is miniscule, totaling only 9 million bases and containing fewer than 3,000 genes, half the gene content of two other sequenced invertebrate-infecting apicomplexans, Porospora gigantea and Gregarina niphandrodes. We found that O. elektroscirrha shares different orthologs with each sequenced relative, suggesting the true set of universally conserved apicomplexan genes is very small indeed. Next, we show that sequencing data from other potential host butterflies can be used to diagnose infection status as well as to study diversity of parasite sequences. We recovered a similarly sized parasite genome from another butterfly, Danaus chrysippus, that was highly diverged from the O. elektroscirrha reference, possibly representing a distinct species. Using these two new genomes, we investigated potential evolutionary response by parasites to toxic phytochemicals their hosts ingest and sequester. Monarch butterflies are well-known to tolerate toxic cardenolides thanks to changes in the sequence of their Type II ATPase sodium pumps. We show that Ophryocystis completely lacks Type II or Type 4 sodium pumps, and related proteins PMCA calcium pumps show extreme sequence divergence compared to other Apicomplexa, demonstrating new avenues of research opened by genome sequencing of non-model Apicomplexa.

African Queens find mates when males are rare.

In butterflies and moths, male-killing endosymbionts are transmitted from infected females via their eggs, and the male progeny then perish. This means that successful transmission of the parasite relies on the successful mating of the host. Paradoxically, at the population level, parasite transmission also reduces the number of adult males present in the final population for infected females to mate with. Here we investigate if successful female mating when males are rare is indeed a likely rate-limiting step in the transmission of male-killing Spiroplasma in the African Monarch, Danaus chrysippus. In Lepidoptera, successful pairings are hallmarked by the transfer of a sperm-containing spermatophore from the male to the female during copulation. Conveniently, this spermatophore remains detectable within the female upon dissection, and thus, spermatophore counts can be used to assess the frequency of successful mating in the field. We used such spermatophore counts to examine if altered sex ratios in the D. chrysippus do indeed affect female mating success. We examined two different field sites in East Africa where males were often rare. Surprisingly, mated females carried an average of 1.5 spermatophores each, regardless of male frequency, and importantly, only 10-20% remained unmated. This suggests that infected females will still be able to mate in the face of either Spiroplasma-mediated male killing and/or fluctuations in adult sex ratio over the wet-dry season cycle. These observations may begin to explain how the male-killing mollicute can still be successfully transmitted in a population where males are rare.

Chromosome Fissions and Fusions Act as Barriers to Gene Flow between Brenthis Fritillary Butterflies.

Chromosome rearrangements are thought to promote reproductive isolation between incipient species. However, it is unclear how often, and under what conditions, fission and fusion rearrangements act as barriers to gene flow. Here we investigate speciation between two largely sympatric fritillary butterflies, Brenthis daphne and Brenthis ino. We use a composite likelihood approach to infer the demographic history of these species from whole-genome sequence data. We then compare chromosome-level genome assemblies of individuals from each species and identify a total of nine chromosome fissions and fusions. Finally, we fit a demographic model where effective population sizes and effective migration rate vary across the genome, allowing us to quantify the effects of chromosome rearrangements on reproductive isolation. We show that chromosomes involved in rearrangements experienced less effective migration since the onset of species divergence and that genomic regions near rearrangement points have a further reduction in effective migration rate. Our results suggest that the evolution of multiple rearrangements in the B. daphne and B. ino populations, including alternative fusions of the same chromosomes, have resulted in a reduction in gene flow. Although fission and fusion of chromosomes are unlikely to be the only processes that have led to speciation between these butterflies, this study shows that these rearrangements can directly promote reproductive isolation and may be involved in speciation when karyotypes evolve quickly.

Rapid and predictable genome evolution across three hybrid ant populations.

Hybridization is frequent in the wild but it is unclear when admixture events lead to predictable outcomes and if so, at what timescale. We show that selection led to correlated sorting of genetic variation rapidly after admixture in 3 hybrid Formica aquilonia × F. polyctena ant populations. Removal of ancestry from the species with the lowest effective population size happened in all populations, consistent with purging of deleterious load. This process was modulated by recombination rate variation and the density of functional sites. Moreover, haplotypes with signatures of positive selection in either species were more likely to fix in hybrids. These mechanisms led to mosaic genomes with comparable ancestry proportions. Our work demonstrates predictable evolution over short timescales after admixture in nature.

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.

Global biogeography of warning coloration in the butterfly Danaus chrysippus.

Warning coloration provides a textbook example of natural selection, but the frequent observation of polymorphism in aposematic species presents an evolutionary puzzle. We investigated biogeography and polymorphism of warning patterns in the widespread butterfly Danaus chrysippus using records from citizen science (n = 5467), museums (n = 8864) and fieldwork (n = 2586). We find that polymorphism in three traits controlled by known mendelian loci is extensive. Broad allele frequency clines, hundreds of kilometres wide, suggest a balance between long-range dispersal and predation of unfamiliar morphs. Mismatched clines for the white hindwing and forewing tip in East Africa are consistent with a previous finding that the black wingtip allele has spread recently in the region through hitchhiking with a heritable endosymbiont. Light/dark background coloration shows more extensive polymorphism. The darker genotype is more common in cooler regions, possibly reflecting a trade-off between thermoregulation and predator warning. Overall, our findings show how studying local adaptation at the global scale provides a more complete picture of the evolutionary forces involved.

Stepwise evolution of a butterfly supergene via duplication and inversion.

Supergenes maintain adaptive clusters of alleles in the face of genetic mixing. Although usually attributed to inversions, supergenes can be complex, and reconstructing the precise processes that led to recombination suppression and their timing is challenging. We investigated the origin of the BC supergene, which controls variation in warning coloration in the African monarch butterfly, Danaus chrysippus. By generating chromosome-scale assemblies for all three alleles, we identified multiple structural differences. Most strikingly, we find that a region of more than 1 million bp underwent several segmental duplications at least 7.5 Ma. The resulting duplicated fragments appear to have triggered four inversions in surrounding parts of the chromosome, resulting in stepwise growth of the region of suppressed recombination. Phylogenies for the inversions are incongruent with the species tree and suggest that structural polymorphisms have persisted for at least 4.1 Myr. In addition to the role of duplications in triggering inversions, our results suggest a previously undescribed mechanism of recombination suppression through independent losses of divergent duplicated tracts. Overall, our findings add support for a stepwise model of supergene evolution involving a variety of structural changes. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Genome assembly of Danaus chrysippus and comparison with the Monarch Danaus plexippus.

Milkweed butterflies in the genus Danaus are studied in a diverse range of research fields including the neurobiology of migration, biochemistry of plant detoxification, host-parasite interactions, evolution of sex chromosomes, and speciation. We have assembled a nearly chromosomal genome for Danaus chrysippus (known as the African Monarch, African Queen, and Plain Tiger) using long-read sequencing data. This species is of particular interest for the study of genome structural change and its consequences for evolution. Comparison with the genome of the North American Monarch Danaus plexippus reveals generally strong synteny but highlights 3 inversion differences. The 3 chromosomes involved were previously found to carry peaks of intraspecific differentiation in D. chrysippus in Africa, suggesting that these inversions may be polymorphic and associated with local adaptation. The D. chrysippus genome is over 40% larger than that of D. plexippus, and nearly all of the additional ∼100 Megabases of DNA comprises repeats. Future comparative genomic studies within this genus will shed light on the evolution of genome architecture.

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.

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.

Signatures of Introgression across the Allele Frequency Spectrum.

The detection of introgression from genomic data is transforming our view of species and the origins of adaptive variation. Among the most widely used approaches to detect introgression is the so-called ABBA-BABA test or D-statistic, which identifies excess allele sharing between nonsister taxa. Part of the appeal of D is its simplicity, but this also limits its informativeness, particularly about the timing and direction of introgression. Here we present a simple extension, D frequency spectrum or DFS, in which D is partitioned according to the frequencies of derived alleles. We use simulations over a large parameter space to show how DFS carries information about various factors. In particular, recent introgression reliably leads to a peak in DFS among low-frequency derived alleles, whereas violation of model assumptions can lead to a lack of signal at low frequencies. We also reanalyze published empirical data from six different animal and plant taxa, and interpret the results in the light of our simulations, showing how DFS provides novel insights. We currently see DFS as a descriptive tool that will augment both simple and sophisticated tests for introgression, but in the future it may be usefully incorporated into probabilistic inference frameworks.

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.

Genetic diversity, demographic history and neo-sex chromosomes in the Critically Endangered Raso lark.

Small effective population sizes could expose island species to inbreeding and loss of genetic variation. Here, we investigate factors shaping genetic diversity in the Raso lark, which has been restricted to a single islet for approximately 500 years, with a population size of a few hundred. We assembled a reference genome for the related Eurasian skylark and then assessed diversity and demographic history using RAD-seq data (75 samples from Raso larks and two related mainland species). We first identify broad tracts of suppressed recombination in females, indicating enlarged neo-sex chromosomes. We then show that genetic diversity across autosomes in the Raso lark is lower than in its mainland relatives, but inconsistent with long-term persistence at its current population size. Finally, we find that genetic signatures of the recent population contraction are overshadowed by an ancient expansion and persistence of a very large population until the human settlement of Cape Verde. Our findings show how genome-wide approaches to study endangered species can help avoid confounding effects of genome architecture on diversity estimates, and how present-day diversity can be shaped by ancient demographic events.