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.

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

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.

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.

Do juvenile developmental and adult body characteristics differ among genotypes at the doublesex locus that controls female-limited Batesian mimicry polymorphism in Papilio memnon?: A test for the "cost of mimicry" hypothesis

Female-limited Batesian mimicry may have evolved because of stronger predation pressure on females than on males, but some physiological costs of mimicry may also hinder the evolution of mimicry in males. In Papilio memnon, which possesses a female-limited Batesian mimicry polymorphism, two alleles at the doublesex (dsx) locus strictly control female phenotypes. To examine whether there are physiological costs associated with mimetic genotypes in the juvenile stage, we compered mortality, juvenile growth and development, and the resultant adult characteristics among three dsx genotypes (HH, Hh, hh) at a constant temperature (25 °C) and two differing day lengths (LD 14:10 and LD 12:12; the latter might induce pupal diapause) by crossing individuals heterozygous (Hh) for the dsx allele. All pupae emerged directly without diapause irrespective of day length. The genotype frequencies of the emerged individuals were consistent with the expected 1:2:1 ratio of HH:Hh:hh. The sex ratio was significantly male-biased in one of two families, but not in the other. We found no effect of genotype on any developmental or adult characteristic, although there were sex differences in most traits. The larval development time was longer and growth rate higher in females than in males; pupal weight, forewing length, and total dry mass of the thorax and abdomen were greater in females, whereas the thoracic mass/abdominal mass ratio was greater in males. We also found that the growth rate was higher and pupal period longer with a short day than with a long day. Overall, we found no evidence for physiological costs associated with the mimetic genotypes. However, it is too early to conclude that no physiological cost of mimicry affects the evolution and maintenance of this female-limited Batesian mimicry polymorphism because we have not studied the adults of different genotypes.

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

Batesian mimicry protects animals from predators when mimics resemble distasteful models. The female-limited Batesian mimicry in Papilio butterflies is controlled by a supergene locus switching mimetic and nonmimetic forms. In Papilio polytes, recent studies revealed that a highly diversified region (HDR) containing doublesex (dsx-HDR) constitutes the supergene with dimorphic alleles and is likely maintained by a chromosomal inversion. In the closely related Papilio memnon, which exhibits a similar mimicry polymorphism, we performed whole-genome sequence analyses in 11 butterflies, which revealed a nearly identical dsx-HDR containing three genes (dsx, Nach-like, and UXT) with dimorphic sequences strictly associated with the mimetic/nonmimetic phenotypes. In addition, expression of these genes, except that of Nach-like in female hind wings, showed differences correlated with phenotype. The dimorphic dsx-HDR in P. memnon is maintained without a chromosomal inversion, suggesting that a separate mechanism causes and maintains allelic divergence in these genes. More abundant accumulation of transposable elements and repetitive sequences in the dsx-HDR than in other genomic regions may contribute to the suppression of chromosomal recombination. Gene trees for Dsx, Nach-like, and UXT indicated that mimetic alleles evolved independently in the two Papilio species. These results suggest that the genomic region involving the above three genes has repeatedly diverged so that two allelic sequences of this region function as developmental switches for mimicry polymorphism in the two Papilio species. The supergene structures revealed here suggest that independent evolutionary processes with different genetic mechanisms have led to parallel evolution of similar female-limited polymorphisms underlying Batesian mimicry in Papilio butterflies.

Temporal dynamics of the mimetic allele frequency at the doublesex locus, which controls polymorphic Batesian mimicry in Papilio memnon butterflies

Tracking allele frequencies is essential for understanding how polymorphisms of adaptive traits are maintained. In Papilio memnon butterflies, which exhibit a female-limited Batesian mimicry polymorphism (wing-pattern polymorphism), two alleles at the doublesex (dsx) locus correspond to mimetic and non-mimetic forms in females; males carry both dsx alleles but display only the non-mimetic form. This polymorphism is thought to be maintained by a negative frequency-dependent selection. By tracking dsx allele frequencies in both sexes at a Taiwanese site over four years, we found that the mimetic allele persists at intermediate frequencies even when the unpalatable model papilionid butterflies (Pachliopta and Atrophaneura species) were very rare or absent. The rates of male mate choice did not differ between the two female forms; neither did insemination number nor age composition, suggesting equivalent reproductive performance of the two forms over time. Our results characterised the temporal dynamics of the mimetic allele frequency in the field for the first time and give insights into underlying processes involved in the persistence of the female-limited Batesian mimicry polymorphism.