Parallel evolution of Batesian mimicry supergene in two Papilio butterflies, P. polytes and P. memnon

The molecular basis of phenotypic evolution: beyond the usual suspects

It has been well documented that mutations in coding DNA or cis-regulatory elements underlie natural phenotypic variation in many organisms. However, the development of sophisticated functional tools in recent years in a wide range of traditionally non-model systems have revealed many 'unusual suspects' in the molecular bases of phenotypic evolution, including upstream open reading frames (uORFs), cryptic splice sites, and small RNAs. Furthermore, large-scale genome sequencing, especially long-read sequencing, has identified a cornucopia of structural variation underlying phenotypic divergence and elucidated the composition of supergenes that control complex multi-trait polymorphisms. In this review article we highlight recent studies that demonstrate this great diversity of molecular mechanisms producing adaptive genetic variation and the panoply of evolutionary paths leading to the 'grandeur of life'.

A micro-RNA is the effector gene of a classic evolutionary hotspot locus

In Lepidoptera (butterflies and moths), the genomic region around the gene cortex is a 'hotspot' locus, repeatedly used to generate intraspecific melanic wing color polymorphisms across 100-million-years of evolution. However, the identity of the effector gene regulating melanic wing color within this locus remains unknown. Here, we show that none of the four candidate protein-coding genes within this locus, including cortex, serve as major effectors. Instead, a micro-RNA (miRNA), mir-193, serves as the major effector across three deeply diverged lineages of butterflies, and its function is conserved in Drosophila. In Lepidoptera, mir-193 is derived from a gigantic long non-coding RNA, ivory, and it functions by directly repressing multiple pigmentation genes. We show that a miRNA can drive repeated instances of adaptive evolution in animals.

The interplay of local adaptation and gene flow may lead to the formation of supergenes

Supergenes are genetic architectures resulting in the segregation of alternative combinations of alleles underlying complex phenotypes. The co-segregation of alleles at linked loci is often facilitated by polymorphic chromosomal rearrangements suppressing recombination locally. Supergenes are involved in many complex polymorphisms, including sexual, colour or behavioural polymorphisms in numerous plants, fungi, mammals, fish, and insects. Despite a long history of empirical and theoretical research, the formation of supergenes remains poorly understood. Here, using a two-island population genetic model, we explore how gene flow and the evolution of overdominant chromosomal inversions may jointly lead to the formation of supergenes. We show that the evolution of inversions in differentiated populations, both under disruptive selection, leads to an increase in frequency of poorly adapted, immigrant haplotypes. Indeed, rare allelic combinations, such as immigrant haplotypes, are more frequently reshuffled by recombination than common allelic combinations, and therefore benefit from the recombination suppression generated by inversions. When an inversion capturing a locally adapted haplotype spreads but is associated with a fitness cost hampering its fixation (e.g. a recessive mutation load), the maintenance of a non-inverted haplotype in the population is enhanced; under certain conditions, the immigrant haplotype persists alongside the inverted local haplotype, while the standard local haplotype disappears. This establishes a stable, local polymorphism with two non-recombining haplotypes encoding alternative adaptive strategies, that is, a supergene. These results bring new light to the importance of local adaptation, overdominance, and gene flow in the formation of supergenes and inversion polymorphisms in general.

The diversification of butterfly wing patterns: progress and prospects

Butterfly wings display rich phenotypic diversity and are associated with complex biological functions, thus serving as an important evolutionary system to address the genetic basis and evolution of phenotypic diversification. We review recent butterfly studies that revealed complex functions underlying diversified wing patterns and describe the genetic and environmental factors involved in wing pattern determinations. These factors lead to inter-specific divergence, genetic polymorphism, and phenotypic plasticity, which in many cases are decided by several key genes. We also summarize the research advances on gene co-option as an important origin of functional complexity and evolutionary novelty. These findings reveal a pattern of evolutionary innovation within a constrained developmental framework during butterfly wing morphogenesis, but further research is required to gain a systematic and comprehensive understanding.

Supergene evolution via gain of auto-regulation

The development of complex phenotypes requires the coordinated action of many genes across space and time, yet many species have evolved the ability to develop multiple discrete, alternate phenotypes1-3. Such polymorphisms are often controlled by supergenes, sets of tightly-linked mutations in one or more loci that function together to produce a complex phenotype4. Although theories of supergene evolution are well-established, the mutations that cause functional differences between supergene alleles remain essentially unknown. doublesex is the master regulator of insect sexual differentiation but functions as a supergene in multiple Papilio swallowtail butterflies, where divergent dsx alleles control development of discrete non-mimetic or mimetic female wing color patterns5-7. Here we demonstrate that the functional elements of the mimetic allele in Papilio alphenor are six new cis-regulatory elements (CREs) spread across 150 kb that are bound by DSX itself. Our findings provide experimental support to classic supergene theory and suggest that the evolution of auto-regulation may provide a simple route to supergene origination and to the co-option of pleiotropic genes into new developmental roles.

Base Composition, Codon Usage, and Patterns of Gene Sequence Evolution in Butterflies

Coding sequence evolution is influenced by both natural selection and neutral evolutionary forces. In many species, the effects of mutation bias, codon usage, and GC-biased gene conversion (gBGC) on gene sequence evolution have not been detailed. Quantification of how these forces shape substitution patterns is therefore necessary to understand the strength and direction of natural selection. Here, we used comparative genomics to investigate the association between base composition and codon usage bias on gene sequence evolution in butterflies and moths (Lepidoptera), including an in-depth analysis of underlying patterns and processes in one species, Leptidea sinapis. The data revealed significant G/C to A/T substitution bias at third codon position with some variation in the strength among different butterfly lineages. However, the substitution bias was lower than expected from previously estimated mutation rate ratios, partly due to the influence of gBGC. We found that A/T-ending codons were overrepresented in most species, but there was a positive association between the magnitude of codon usage bias and GC-content in third codon positions. In addition, the tRNA-gene population in L. sinapis showed higher GC-content at third codon positions compared to coding sequences in general and less overrepresentation of A/T-ending codons. There was an inverse relationship between synonymous substitutions and codon usage bias indicating selection on synonymous sites. We conclude that the evolutionary rate in Lepidoptera is affected by a complex interaction between underlying G/C -> A/T mutation bias and partly counteracting fixation biases, predominantly conferred by overall purifying selection, gBGC, and selection on codon usage.

Functional unit of supergene in female-limited Batesian mimicry of Papilio polytes

Supergenes are sets of genes and genetic elements that are inherited like a single gene and control complex adaptive traits, but their functional roles and units are poorly understood. In Papilio polytes, female-limited Batesian mimicry is thought to be regulated by a ∼130 kb inversion region (highly diversified region: HDR) containing 3 genes, UXT, U3X, and doublesex (dsx) which switches non-mimetic and mimetic types. To determine the functional unit, we here performed electroporation-mediated RNAi analyses (and further Crispr/Cas9 for UXT) of genes within and flanking the HDR in pupal hindwings. We first clarified that non-mimetic dsx-h had a function to form the non-mimetic trait in female and only dsx-H isoform 3 had an important function in the formation of mimetic traits. Next, we found that UXT was involved in making mimetic-type pale-yellow spots and adjacent gene sir2 in making red spots in hindwings, both of which refine more elaborate mimicry. Furthermore, downstream gene networks of dsx, U3X, and UXT screened by RNA sequencing showed that U3X upregulated dsx-H expression and repressed UXT expression. These findings demonstrate that a set of multiple genes, not only inside but also flanking HDR, can function as supergene members, which extends the definition of supergene unit than we considered before. Also, our results indicate that dsx functions as the switching gene and some other genes such as UXT and sir2 within the supergene unit work as the modifier gene.

Widespread chromosomal rearrangements preceded genetic divergence in a monitor lizard, Varanus acanthurus (Varanidae)

Chromosomal rearrangements are often associated with local adaptation and speciation because they suppress recombination, and as a result, rearrangements have been implicated in disrupting gene flow. Although there is strong evidence to suggest that chromosome rearrangements are a factor in genetic isolation of divergent populations, the underlying mechanism remains elusive. Here, we applied an integrative cytogenetics and genomics approach testing whether chromosomal rearrangements are the initial process, or a consequence, of population divergence in the dwarf goanna, Varanus acanthurus. Specifically, we tested whether chromosome rearrangements are indicators of genetic barriers that can be used to identify divergent populations by looking at gene flow within and between populations with rearrangements. We found that gene flow was present between individuals with chromosome rearrangements within populations, but there was no gene flow between populations that had similar chromosome rearrangements. Moreover, we identified a correlation between reduced genetic variation in populations with a higher frequency of homozygous submetacentric individuals. These findings suggest that chromosomal rearrangements were widespread prior to divergence, and because we found populations with higher frequencies of submetacentric chromosomes were associated with lower genetic diversity, this could indicate that polymorphisms within populations are early indicators of genetic drift.

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.

Genome-Wide Analysis of Odorant and Gustatory Receptors in Six Papilio Butterflies (Lepidoptera: Papilionidae)

The chemical interactions of insects and host plants are shaping the evolution of chemosensory receptor gene families. However, the correlation between host range and chemoreceptor gene repertoire sizes is still elusive in Papilionidae. Here, we addressed the issue of whether host plant diversities are correlated with the expansions of odorant (ORs) or gustatory (GRs) receptors in six Papilio butterflies. By combining genomics, transcriptomics and bioinformatics approaches, 381 ORs and 328 GRs were annotated in the genomes of a generalist P. glaucus and five specialists, P. xuthus, P. polytes, P. memnon, P. machaon and P. dardanus. Orthologous ORs or GRs in Papilio had highly conserved gene structure. Five Papilio specialists exhibited a similar frequency of intron lengths for ORs or GRs, but which was different from those in the generalist. Phylogenetic analysis revealed 60 orthologous OR groups, 45 of which shared one-to-one relationships. Such a single gene in each butterfly also occurred in 26 GR groups. Intriguingly, bitter GRs had fewer introns than other GRs and clustered into a large clade. Focusing on the two chemoreceptor gene families in P. xuthus, most PxutORs (52/58) were expressed in antennae and 31 genes in reproductive tissues. Eleven out of 28 foretarsus-expressed PxutGRs were female-biased genes, as strong candidates for sensing oviposition stimulants. These results indicate that the host range may not shape the large-scale expansions of ORs and GRs in Papilio butterflies and identify important molecular targets involved in olfaction, oviposition or reproduction in P. xuthus.

Genomic architecture and functional unit of mimicry supergene in female limited Batesian mimic Papilio butterflies

It has long been suggested that dimorphic female-limited Batesian mimicry of two closely related Papilio butterflies, Papilio memnon and Papilio polytes, is controlled by supergenes. Whole-genome sequencing, genome-wide association studies and functional analyses have recently identified mimicry supergenes, including the doublesex (dsx) gene. Although supergenes of both the species are composed of highly divergent regions between mimetic and non-mimetic alleles and are located at the same chromosomal locus, they show critical differences in genomic architecture, particularly with or without an inversion: P. polytes has an inversion, but P. memnon does not. This review introduces and compares the detailed genomic structure of mimicry supergenes in two Papilio species, including gene composition, repetitive sequence composition, breakpoint/boundary site structure, chromosomal inversion and linkage disequilibrium. Expression patterns and functional analyses of the respective genes within or flanking the supergene suggest that dsx and other genes are involved in mimetic traits. In addition, structural comparison of the corresponding region for the mimicry supergene among further Papilio species suggests three scenarios for the evolution of the mimicry supergene between the two Papilio species. The structural features revealed in the Papilio mimicry supergene provide insight into the formation, maintenance and evolution of supergenes. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Genomic architecture of supergenes: connecting form and function

Supergenes are tightly linked sets of loci that are inherited together and control complex phenotypes. While classical supergenes-governing traits such as wing patterns in Heliconius butterflies or heterostyly in Primula-have been studied since the Modern Synthesis, we still understand very little about how they evolve and persist in nature. The genetic architecture of supergenes is a critical factor affecting their evolutionary fate, as it can change key parameters such as recombination rate and effective population size, potentially redirecting molecular evolution of the supergene in addition to the surrounding genomic region. To understand supergene evolution, we must link genomic architecture with evolutionary patterns and processes. This is now becoming possible with recent advances in sequencing technology and powerful forward computer simulations. The present theme issue brings together theoretical and empirical papers, as well as opinion and synthesis papers, which showcase the architectural diversity of supergenes and connect this to critical processes in supergene evolution, such as polymorphism maintenance and mutation accumulation. Here, we summarize those insights to highlight new ideas and methods that illuminate the path forward for the study of supergenes in nature. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

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

doublesex Controls Both Hindwing and Abdominal Mimicry Traits in the Female-Limited Batesian Mimicry of Papilio memnon

Papilio butterflies are known to possess female-limited Batesian mimicry polymorphisms. In Papilio memnon, females have mimetic and non-mimetic forms, whereas males are monomorphic and non-mimetic. Mimetic females are characterized by color patterns and tails in the hindwing and yellow abdomens. Recently, an analysis of whole-genome sequences has shown that an approximately 160 kb region of chromosome 25 is responsible for mimicry and has high diversity between mimetic (A) and non-mimetic (a) alleles (highly diversified region: HDR). The HDR includes three genes, UXT, doublesex (dsx), and Nach-like, but the functions of these genes are unknown. Here, we investigated the function of dsx, a gene involved in sexual differentiation, which is expected to be functionally important for hindwing and abdominal mimetic traits in P. memnon. Expression analysis by reverse transcription quantitative PCR (RT-qPCR) and RNA sequencing showed that mimetic dsx (dsx-A) was highly expressed in the hindwings in the early pupal stage. In the abdomen, both dsx-A and dsx-a were highly expressed during the early pupal stage. When dsx was knocked down using small interfering RNAs (siRNAs) designed in the common region of dsx-A and dsx-a, a male-like pattern appeared on the hindwings of mimetic and non-mimetic females. Similarly, when dsx was knocked down in the abdomen, the yellow scales characteristic of mimetic females changed to black. Furthermore, when dsx-a was specifically knocked down, the color pattern of the hindwings changed, as in the case of dsx knockdown in non-mimetic females but not mimetic females. These results suggest that dsx-a is involved in color pattern formation on the hindwings of non-mimetic females, whereas dsx-A is involved in hindwing and abdominal mimetic traits. dsx was involved in abdominal and hindwing mimetic traits, but dsx expression patterns in the hindwing and abdomen were different, suggesting that different regulatory mechanisms may exist. Our study is the first to show that the same gene (dsx) regulates both the hindwing and abdominal mimetic traits. This is the first functional analysis of abdominal mimicry in butterflies.

On the evolution of mimicry in avian nestlings

Batesian mimicry (BM), where a nontoxic species resembles a toxic species with aposematic coloring, has been recently described for a Neotropical species of the suboscine passerine (Laniocera hypopyrra). Understanding the order and series in which these characteristics evolved is unknown and requires character information from closely related taxa. Here, we trace the origin of mimetic traits and how they evolved by examining antipredator characteristics using images and other field‐collected trait data from nest and nestlings along with data available in the literature for the Laniisominae clade and closely related taxa. We found that morphological modifications of the downy feathers appeared first in the broader clade leading to the Laniisominae clade followed by further morphological and behavioral characteristics within the Laniisominae clade leading to the full BM. Images of nestlings in the Laniisominae and closely related clades demonstrated the extent of antipredator and camouflage characteristics. We found a complex set of behavioral and morphological traits in this clade for reducing predation from hiding to camouflage to mimicry. We further propose the evolution of two distinctive mimicry strategies in the Laniisominae clade: (1) Batesian Mimicry, as described above and (2) Masquerade, resemblance to inedible objects commonly found in their local environment. This complex set of antipredator traits shed light on the diversity of antipredator characteristics in avian nestlings, particularly in neotropical areas where the avian diversity is highest. Unfortunately, there are hundreds of species in the neotropics that lack basic natural history information on nesting traits, and therefore, we are likely missing critical information on the diversity of antipredator characteristics across avian nestlings.

The effect of the doublesex gene in body colour masculinization of the damselfly Ischnura senegalensis

Odonata species display a remarkable diversity of colour patterns, including intrasexual polymorphisms. In the damselfly (Ischnura senegalensis), the expression of a sex-determining transcription factor, the doublesex (Isdsx) gene is reportedly associated with female colour polymorphism (CP) (gynomorph for female-specific colour and andromorph for male-mimicking colour). Here, the function of Isdsx in thoracic coloration was investigated by electroporation-mediated RNA interference (RNAi). RNAi of the Isdsx common region in males and andromorphic females reduced melanization and thus changed the colour pattern into that of gynomorphic females, while the gynomorphic colour pattern was not affected. By contrast, RNAi against the Isdsx long isoform produced no changes, suggesting that the Isdsx short isoform is important for body colour masculinization in both males and andromorphic females. When examining the expression levels of five genes with differences between sexes and female morphs, two melanin-suppressing genes, black and ebony, were expressed at higher levels in the Isdsx RNAi body area than a control area. Therefore, the Isdsx short isoform may induce thoracic colour differentiation by suppressing black and ebony, thereby generating female CP in I. senegalensis. These findings contribute to the understanding of the molecular and evolutionary mechanisms underlying female CP in Odonata.

Parallel Genetic Origin of Foot Feathering in Birds

Understanding the genetic basis of similar phenotypes shared between lineages is a long-lasting research interest. Even though animal evolution offers many examples of parallelism, for many phenotypes little is known about the underlying genes and mutations. We here use a combination of whole-genome sequencing, expression analyses, and comparative genomics to study the parallel genetic origin of ptilopody (Pti) in chicken. Ptilopody (or foot feathering) is a polygenic trait that can be observed in domesticated and wild avian species and is characterized by the partial or complete development of feathers on the ankle and feet. In domesticated birds, ptilopody is easily selected to fixation, though extensive variation in the type and level of feather development is often observed. By means of a genome-wide association analysis, we identified two genomic regions associated with ptilopody. At one of the loci, we identified a 17-kb deletion affecting PITX1 expression, a gene known to encode a transcription regulator of hindlimb identity and development. Similarly to pigeon, at the second loci, we observed ectopic expression of TBX5, a gene involved in forelimb identity and a key determinant of foot feather development. We also observed that the trait evolved only once as foot-feathered birds share the same haplotype upstream TBX5. Our findings indicate that in chicken and pigeon ptilopody is determined by the same set of genes that affect similar molecular pathways. Our study confirms that ptilopody has evolved through parallel evolution in chicken and pigeon.

Evolution of sexual development and sexual dimorphism in insects

Most animal species consist of two distinct sexes. At the morphological, physiological, and behavioral levels the differences between males and females are numerous and dramatic, yet at the genomic level they are often slight or absent. This disconnect is overcome because simple genetic differences or environmental signals are able to direct the sex-specific expression of a shared genome. A canonical picture of how this process works in insects emerged from decades of work on Drosophila. But recent years have seen an explosion of molecular-genetic and developmental work on a broad range of insects. Drawing these studies together, we describe the evolution of sexual dimorphism from a comparative perspective and argue that insect sex determination and differentiation systems are composites of rapidly evolving and highly conserved elements.

Disentangling Population History and Character Evolution among Hybridizing Lineages

Understanding the origin and maintenance of adaptive phenotypic novelty is a central goal of evolutionary biology. However, both hybridization and incomplete lineage sorting can lead to genealogical discordance between the regions of the genome underlying adaptive traits and the remainder of the genome, decoupling inferences about character evolution from population history. Here, to disentangle these effects, we investigated the evolutionary origins and maintenance of Batesian mimicry between North American admiral butterflies (Limenitis arthemis) and their chemically defended model (Battus philenor) using a combination of de novo genome sequencing, whole-genome resequencing, and statistical introgression mapping. Our results suggest that balancing selection, arising from geographic variation in the presence or absence of the unpalatable model, has maintained two deeply divergent color patterning haplotypes that have been repeatedly sieved among distinct mimetic and nonmimetic lineages of Limenitis via introgressive hybridization.

The Genomic Architecture and Evolutionary Fates of Supergenes

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

Population genetic structure and evolution of Batesian mimicry in Papilio polytes from the Ryukyu Islands, Japan, analyzed by genotyping-by-sequencing

Batesian mimicry is a striking example of Darwinian evolution, in which a mimetic species resembles toxic or unpalatable model species, thereby receiving protection from predators. In some species exhibiting Batesian mimicry, nonmimetic individuals coexist as polymorphism in the same population despite the benefits of mimicry. In a previous study, we proposed that the abundance of mimics is limited by that of the models, leading to polymorphic Batesian mimicry in the swallowtail butterfly, Papilio polytes, on the Ryukyu Islands in Japan. We found that their mimic ratios (MRs), which varied among the Islands, were explained by the model abundance of each habitat, rather than isolation by distance or phylogenetic constraint based on the mitochondrial DNA (mtDNA) analysis. In the present study, this possibility was reexamined based on hundreds of nuclear single nucleotide polymorphisms (SNPs) of 93 P. polytes individuals from five Islands of the Ryukyus. We found that the population genetic and phylogenetic structures of P. polytes largely corresponded to the geographic arrangement of the habitat Islands, and the genetic distances among island populations show significant correlation with the geographic distances, which was not evident by the mtDNA-based analysis. A partial Mantel test controlling for the present SNP-based genetic distances revealed that the MRs of P. polytes were strongly correlated with the model abundance of each island, implying that negative frequency-dependent selection interacting with model species shaped and maintained the mimetic polymorphism. Taken together, our results support the possibility that predation pressure, not isolation by distance or other neutral factors, is a major driving force of evolution of the Batesian mimicry in P. polytes from the Ryukyus.

Genetic switch in UV response of mimicry-related pale-yellow colors in Batesian mimic butterfly, Papilio polytes

In a Batesian mimic butterfly Papilio polytes, mimetic females resemble an unpalatable model, Pachliopta aristolochiae, but exhibit a different color pattern from nonmimetic females and males. In particular, the pale-yellow region on hind wings, which correspondingly sends important putative signals for mimicry and mate preference, is different in shape and chemical features between nonmimetic and mimetic morphs. Recently, we found that mimetic-type doublesex [dsx (H)] causes mimetic traits; however, the control of dimorphic pale-yellow colors remains unclear. Here, we revealed that dsx (H) switched the pale-yellow colors from UV-excited fluorescent type (nonmimetic) to UV-reflecting type (mimetic), by repressing the papiliochrome II synthesis genes and nanostructural changes in wing scales. Photoreceptor reactivities showed that some birds and butterflies could effectively recognize mimetic and nonmimetic pale-yellow colors, suggesting that a genetic switch in the UV response of pale-yellow colors may play essential roles in establishing the dimorphic female-limited Batesian mimicry.

Batesian mimicry has evolved with deleterious effects of the pleiotropic gene doublesex

Dimorphic female-limited Batesian mimicry in the swallowtail butterfly Papilio polytes is regulated by the supergene locus H, harbouring the mimetic (H) and non-mimetic (h) doublesex (dsx) gene. In the present study, we demonstrated that dsx-H negatively affects the number of eggs laid, hatching rate, larval survival rate, and adult lifespan. When crossed with hh males, the number of eggs laid of mimetic females (genotype HH) was lower than that of non-mimetic females (hh). Moreover, hh and Hh females laid fewer eggs when crossed with HH males. The hatching and larval survival rates were lower when both female and male parents harboured dsx-H. The adult lifespan of HH females was shorter than that of hh females, while it was similar in males regardless of the genotype. These findings suggest the presence of a cost-benefit balance of Batesian mimicry, which is evolved to avoid predation but is accompanied by physiological deficits, in this species.

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

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

Genome assembly of the basket willow, Salix viminalis, reveals earliest stages of sex chromosome expansion

BACKGROUND

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear.

RESULTS

Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes.

CONCLUSIONS

Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

Genome assembly of the basket willow, Salix viminalis, reveals earliest stages of sex chromosome expansion

BACKGROUND

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear.

RESULTS

Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes.

CONCLUSIONS

Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

Mimicry diversification in Papilio dardanus via a genomic inversion in the regulatory region of engrailed-invected

Polymorphic Batesian mimics exhibit multiple protective morphs that each mimic a different noxious model. Here, we study the genomic transitions leading to the evolution of different mimetic wing patterns in the polymorphic Mocker Swallowtail Papilio dardanus. We generated a draft genome (231 Mb over 30 chromosomes) and re-sequenced individuals of three morphs. Genome-wide single nucleotide polymorphism (SNP) analysis revealed elevated linkage disequilibrium and divergence between morphs in the regulatory region of engrailed, a developmental gene previously implicated in the mimicry switch. The diverged region exhibits a discrete chromosomal inversion (of 40 kb) relative to the ancestral orientation that is associated with the cenea morph, but not with the bottom-recessive hippocoonides morph or with non-mimetic allopatric populations. The functional role of this inversion in the expression of the novel phenotype is currently unknown, but by preventing recombination, it allows the stable inheritance of divergent alleles enabling geographic spread and local coexistence of multiple adaptive morphs.

The genomics of coloration provides insights into adaptive evolution

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

Butterfly Mimicry Polymorphisms Highlight Phylogenetic Limits of Gene Reuse in the Evolution of Diverse Adaptations

Some genes have repeatedly been found to control diverse adaptations in a wide variety of organisms. Such gene reuse reveals not only the diversity of phenotypes these unique genes control but also the composition of developmental gene networks and the genetic routes available to and taken by organisms during adaptation. However, the causes of gene reuse remain unclear. A small number of large-effect Mendelian loci control a huge diversity of mimetic butterfly wing color patterns, but reasons for their reuse are difficult to identify because the genetic basis of mimicry has primarily been studied in two systems with correlated factors: female-limited Batesian mimicry in Papilio swallowtails (Papilionidae) and non-sex-limited Müllerian mimicry in Heliconius longwings (Nymphalidae). Here, we break the correlation between phylogenetic relationship and sex-limited mimicry by identifying loci controlling female-limited mimicry polymorphism Hypolimnas misippus (Nymphalidae) and non-sex-limited mimicry polymorphism in Papilio clytia (Papilionidae). The Papilio clytia polymorphism is controlled by the genome region containing the gene cortex, the classic P supergene in Heliconius numata, and loci controlling color pattern variation across Lepidoptera. In contrast, female-limited mimicry polymorphism in Hypolimnas misippus is associated with a locus not previously implicated in color patterning. Thus, although many species repeatedly converged on cortex and its neighboring genes over 120 My of evolution of diverse color patterns, female-limited mimicry polymorphisms each evolved using a different gene. Our results support conclusions that gene reuse occurs mainly within ∼10 My and highlight the puzzling diversity of genes controlling seemingly complex female-limited mimicry polymorphisms.

A shared genetic basis of mimicry across swallowtail butterflies points to ancestral co-option of doublesex

Uncovering whether convergent adaptations share a genetic basis is consequential for understanding the evolution of phenotypic diversity. This information can help us understand the extent to which shared ancestry or independent evolution shape adaptive phenotypes. In this study, we first ask whether the same genes underlie polymorphic mimicry in Papilio swallowtail butterflies. By comparing signatures of genetic variation between polymorphic and monomorphic species, we then investigate how ancestral variation, hybridization, and independent evolution contributed to wing pattern diversity in this group. We report that a single gene, doublesex (dsx), controls mimicry across multiple taxa, but with species-specific patterns of genetic differentiation and linkage disequilibrium. In contrast to widespread examples of phenotypic evolution driven by introgression, our analyses reveal distinct mimicry alleles. We conclude that mimicry evolution in this group was likely facilitated by ancestral polymorphism resulting from early co-option of dsx as a mimicry locus, and that evolutionary turnover of dsx alleles may underlie the wing pattern diversity of extant polymorphic and monomorphic lineages.

An Ancient and Eroded Social Supergene Is Widespread across Formica Ants

Supergenes, clusters of tightly linked genes, play a key role in the evolution of complex adaptive variation [1, 2]. Although supergenes have been identified in many species, we lack an understanding of their origin, evolution, and persistence [3]. Here, we uncover 20-40 Ma of evolutionary history of a supergene associated with polymorphic social organization in Formica ants [4]. We show that five Formica species exhibit homologous divergent haplotypes spanning 11 Mbp on chromosome 3. Despite the supergene's size, only 142 single nucleotide polymorphisms (SNPs) consistently distinguish alternative supergene haplotypes across all five species. These conserved trans-species SNPs are localized in a small number of disjunct clusters distributed across the supergene. This unexpected pattern of divergence indicates that the Formica supergene does not follow standard models of sex chromosome evolution, in which distinct evolutionary strata reflect an expanding region of suppressed recombination [5]. We propose an alternative "eroded strata model" in which clusters of conserved trans-species SNPs represent functionally important areas maintained by selection in the face of rare recombination between ancestral haplotypes. The comparison of whole-genome sequences across 10 additional Formica species reveals that the most conserved region of the supergene contains a transcription factor essential for motor neuron development in Drosophila [6]. The discovery that a very small portion of this large and ancient supergene harbors conserved trans-species SNPs linked to colony social organization suggests that the ancestral haplotypes have been eroded by recombination, with selection preserving differentiation at one or a few genes generating alternative social organization.

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

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

Chromosomal-level reference genome of Chinese peacock butterfly (Papilio bianor) based on third-generation DNA sequencing and Hi-C analysis

BACKGROUND

Papilio bianor Cramer, 1777 (commonly known as the Chinese peacock butterfly) (Insecta, Lepidoptera, Papilionidae) is a widely distributed swallowtail butterfly with a wide number of geographic populations ranging from the southeast of Russia to China, Japan, India, Vietnam, Myanmar, and Thailand. Its wing color consists of both pigmentary colored scales (black, reddish) and structural colored scales (iridescent blue or green dust). A high-quality reference genome of P. bianor is an important foundation for investigating iridescent color evolution, phylogeography, and the evolution of swallowtail butterflies.

FINDINGS

We obtained a chromosome-level de novo genome assembly of the highly heterozygous P. bianor using long Pacific Biosciences sequencing reads and high-throughput chromosome conformation capture technology. The final assembly is 421.52 Mb on 30 chromosomes (29 autosomes and 1 Z sex chromosome) with 13.12 Mb scaffold N50. In total, 15,375 protein-coding genes and 233.09 Mb of repetitive sequences were identified. Phylogenetic analyses indicated that P. bianor separated from a common ancestor of swallowtails ∼23.69-36.04 million years ago. Demographic history suggested that the population expansion of this species from the last interglacial period to the last glacial maximum possibly resulted from its decreased natural enemies and its adaptation to climate change during the glacial period.

CONCLUSIONS

We present a high-quality chromosome-level reference genome of P. bianor using long-read single-molecule sequencing and Hi-C-based chromatin interaction maps. Our results lay the foundation for exploring the genetic basis of special biological features of P. bianor and also provide a useful data source for comparative genomics and phylogenomics among butterflies and moths.

Genomic evidence of human selection on Vavilovian mimicry

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

Conservation and flexibility in the gene regulatory landscape of heliconiine butterfly wings

BACKGROUND

Many traits evolve by cis-regulatory modification, by which changes to noncoding sequences affect the binding affinity for available transcription factors and thus modify the expression profile of genes. Multiple examples of cis-regulatory evolution have been described at pattern switch genes responsible for butterfly wing pattern polymorphism, including in the diverse neotropical genus Heliconius, but the identities of the factors that can regulate these switch genes have not been identified.

RESULTS

We investigated the spatial transcriptomic landscape across the wings of three closely related butterfly species, two of which have a convergently evolved co-mimetic pattern and the other having a divergent pattern. We identified candidate factors for regulating the expression of wing patterning genes, including transcription factors with a conserved expression profile in all three species, and others, including both transcription factors and Wnt pathway genes, with markedly different profiles in each of the three species. We verified the conserved expression profile of the transcription factor homothorax by immunofluorescence and showed that its expression profile strongly correlates with that of the selector gene optix in butterflies with the Amazonian forewing pattern element 'dennis.'

CONCLUSION

Here we show that, in addition to factors with conserved expression profiles like homothorax, there are also a variety of transcription factors and signaling pathway components that appear to vary in their expression profiles between closely related butterfly species, highlighting the importance of genome-wide regulatory evolution between species.

The mimetic wing pattern of Papilio polytes butterflies is regulated by a doublesex-orchestrated gene network

The swallowtail butterfly Papilio polytes is sexually dimorphic and exhibits female-limited Batesian mimicry. This species also has two female forms, a non-mimetic form with male-like wing patterns, and a mimetic form resembling an unpalatable model, Pachliopta aristolochiae. The mimicry locus H constitutes a dimorphic Mendelian 'supergene', including a transcription factor gene doublesex (dsx). However, how the mimetic-type dsx (dsx-H) orchestrates the downstream gene network and causes the mimetic traits remains unclear. Here we performed RNA-seq-based gene screening and found that Wnt1 and Wnt6 are up-regulated by dsx-H during the early pupal stage and are involved in the red/white pigmentation and patterning of mimetic female wings. In contrast, a homeobox gene abdominal-A is repressed by dsx-H and involved in the non-mimetic colouration pattern. These findings suggest that dual regulation by dsx-H, induction of mimetic gene networks and repression of non-mimetic gene networks, is essential for the switch from non-mimetic to mimetic pattern in mimetic female wings.

Genetic architecture and sex-specific selection govern modular, male-biased evolution of doublesex

doublesex regulates early embryonic sex differentiation in holometabolous insects, along with the development of species-, sex-, and morph-specific adaptations during pupal stages. How does a highly conserved gene with a critical developmental role also remain functionally dynamic enough to gain ecologically important adaptations that are divergent in sister species? We analyzed patterns of exon-level molecular evolution and protein structural homology of doublesex from 145 species of four insect orders representing 350 million years of divergence. This analysis revealed that evolution of doublesex was governed by a modular architecture: Functional domains and female-specific regions were highly conserved, whereas male-specific sequences and protein structures evolved up to thousand-fold faster, with sites under pervasive and/or episodic positive selection. This pattern of sex bias was reversed in Hymenoptera. Thus, highly conserved yet dynamic master regulators such as doublesex may partition specific conserved and novel functions in different genic modules at deep evolutionary time scales.