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.

Whole-chromosome hitchhiking driven by a male-killing endosymbiont.

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

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.

Karyotypes versus Genomes: The Nymphalid Butterflies Melitaea cinxia, Danaus plexippus, and D. chrysippus.

The number of sequenced lepidopteran genomes is increasing rapidly. However, the corresponding assemblies rarely represent whole chromosomes and generally also lack the highly repetitive W sex chromosome. Knowledge of the karyotypes can facilitate genome assembly and further our understanding of sex chromosome evolution in Lepidoptera. Here, we describe the karyotypes of the Glanville fritillary Melitaea cinxia (n = 31), the monarch Danaus plexippus (n = 30), and the African queen D. chrysippus (2n = 60 or 59, depending on the source population). We show by FISH that the telomeres are of the (TTAGG)n type, as found in most insects. M. cinxia and D. plexippus have "conventional" W chromosomes which are heterochromatic in meiotic and somatic cells. In D. chrysippus, the W is inconspicuous. Neither telomeres nor W chromosomes are represented in the published genomes of M. cinxia and D. plexippus. Representation analysis in sequenced female and male D. chrysippus genomes detected an evolutionarily old autosome-Z chromosome fusion in Danaus. Conserved synteny of whole chromosomes, so called "macro synteny", in Lepidoptera permitted us to identify the chromosomes involved in this fusion. An additional and more recent sex chromosome fusion was found in D. chrysippus by karyotype analysis and classical genetics. In a hybrid population between 2 subspecies, D. c. chrysippus and D. c. dorippus, the W chromosome was fused to an autosome that carries a wing colour locus. Thus, cytogenetics and the present state of genome data complement one another to reveal the evolutionary history of the species.

A neo-W chromosome in a tropical butterfly links colour pattern, male-killing, and speciation.

Sexually antagonistic selection can drive both the evolution of sex chromosomes and speciation itself. The tropical butterfly the African Queen, Danaus chrysippus, shows two such sexually antagonistic phenotypes, the first being sex-linked colour pattern, the second, susceptibility to a male-killing, maternally inherited mollicute, Spiroplasma ixodeti, which causes approximately 100% mortality in male eggs and first instar larvae. Importantly, this mortality is not affected by the infection status of the male parent and the horizontal transmission of Spiroplasma is unknown. In East Africa, male-killing of the Queen is prevalent in a narrow hybrid zone centred on Nairobi. This hybrid zone separates otherwise allopatric subspecies with different colour patterns. Here we show that a neo-W chromosome, a fusion between the W (female) chromosome and an autosome that controls both colour pattern and male-killing, links the two phenotypes thereby driving speciation across the hybrid zone. Studies of the population genetics of the neo-W around Nairobi show that the interaction between colour pattern and male-killer susceptibility restricts gene flow between two subspecies of D. chrysippus Our results demonstrate how a complex interplay between sex, colour pattern, male-killing, and a neo-W chromosome, has set up a genetic 'sink' that keeps the two subspecies apart. The association between the neo-W and male-killing thus provides a 'smoking gun' for an ongoing speciation process.

Penetrance interactions of colour pattern loci in the African Monarch and their implications for the evolution of dominance.

Scoring the penetrance of heterozygotes in complex phenotypes, like colour pattern, is difficult and complicates the analysis of systems in which dominance is incomplete or evolving. The African Monarch (Danaus chrysippus) represents an example where colour pattern heterozygotes, formed in the contact zone between the different subspecies, show such intermediate dominance. Colour pattern in this aposematic butterfly is controlled by three loci A, B and C. The B and C loci are closely linked in a B/C supergene and significant interaction of B and C phenotypes is therefore expected via linkage alone. The A locus, however, is not linked to B/C and is found on a different chromosome. To study interactions between these loci we generated colour pattern heterozygotes by crossing males and females bearing different A and B/C genotypes, collected from different parts of Africa. We derived a novel scoring system for the expressivity of the heterozygotes and, as predicted, we found significant interactions between the genotypes of the closely linked B and C loci. Surprisingly, however, we also found highly significant interactions between C and the unlinked A locus, modifications that generally increased the resemblance of heterozygotes to homozygous ancestors. In contrast, we found no difference in the penetrance of any of the corresponding heterozygotes from crosses conducted either in allopatry or sympatry, in reciprocal crosses of males and females, or in the presence or absence of endosymbiont mediated male-killing or its associated neoW mediated sex-linkage of colour pattern. Together, this data supports the idea that the different colour morphs of the African Monarch meet transiently in the East African contact zone and that genetic modifiers act to mask inappropriate expression of colour patterns in the incorrect environments.

The Hypolimnas misippus Genome Supports a Common Origin of the W Chromosome in Lepidoptera.

Moths and butterflies (Lepidoptera) have a heterogametic sex chromosome system with females carrying ZW chromosomes and males ZZ. The lack of W chromosomes in early-diverging lepidopteran lineages has led to the suggestion of an ancestral Z0 system in this clade and a B chromosome origin of the W. This contrasts with the canonical model of W chromosome evolution in which the W would have originated from the same homologous autosomal pair as the Z chromosome. Despite the distinct models proposed, the rapid evolution of the W chromosome has hindered the elucidation of its origin. Here, we present high-quality, chromosome-level genome assemblies of 2 Hypolimnas species (Hypolimnas misippus and Hypolimnas bolina) and use the H. misippus assembly to explore the evolution of W chromosomes in butterflies and moths. We show that in H. misippus, the W chromosome has higher similarity to the Z chromosome than any other chromosome, which could suggest a possible origin from the same homologous autosome pair as the Z chromosome. However, using genome assemblies of closely related species (ditrysian lineages) containing assembled W chromosomes, we present contrasting evidence suggesting that the W chromosome might have evolved from a B chromosome instead. Crucially, by using a synteny analysis to infer homology, we show that W chromosomes are likely to share a common evolutionary origin in Lepidoptera. This study highlights the difficulty of studying the evolution of W chromosomes and contributes to better understanding its evolutionary origins.

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.

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'.

Hybrid effects in field populations of the African monarch butterfly, Danaus chrysippus (L.) (Lepidoptera: Nymphalidae).

Heterosis, Haldane and Bateson-Dobzhansky-Muller effects have been widely documented amongst a range of plants and animals. However, typically these effects are shown by taking parents of known genotype into the laboratory and measuring components of the F1 progeny under laboratory conditions. This leaves in doubt the real significance of such effects in the field. Here we use the well-known colour pattern genotypes of the African monarch or queen (Danaus chrysippus), which also control wing length, to test these effects both in the laboratory and in a contact zone in the field. By measuring the wing lengths in animals of known colour pattern genotype we show clear evidence for all three hybrid effects at the A and BC colour patterning loci, and importantly, that these same effects persist in the same presumptive F1s when measured in hybrid populations in the field. This demonstrates the power of a system in which genotypes can be directly inferred in the field and highlights that all three hybrid effects can be seen in the East African contact zone of this fascinating butterfly.

Neo Sex Chromosomes, Colour Polymorphism and Male-Killing in the African Queen Butterfly, Danaus chrysippus (L.).

Danaus chrysippus (L.), one of the world's commonest butterflies, has an extensive range throughout the Old-World tropics. In Africa it is divided into four geographical subspecies which overlap and hybridise freely in the East African Rift: Here alone a male-killing (MK) endosymbiont, Spiroplasma ixodetis, has invaded, causing female-biased populations to predominate. In ssp. chrysippus, inside the Rift only, an autosome carrying a colour locus has fused with the W chromosome to create a neo-W chromosome. A total of 40-100% of Rift females are neo-W and carry Spiroplasma, thus transmitting a linked, matrilineal neo-W, MK complex. As neo-W females have no sons, half the mother's genes are lost in each generation. Paradoxically, although neo-W females have no close male relatives and are thereby forced to outbreed, MK restricts gene flow between subspecies and may thus promote speciation. The neo-W chromosome originated in the Nairobi region around 2.2 k years ago and subsequently spread throughout the Rift contact zone in some 26 k generations, possibly assisted by not having any competing brothers. Our work on the neo-W chromosome, the spread of Spiroplasma and possible speciation is ongoing.