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

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.

Hierarchical morphogenesis of swallowtail butterfly wing scale nanostructures

The study of color patterns in the animal integument is a fundamental question in biology, with many lepidopteran species being exemplary models in this endeavor due to their relative simplicity and elegance. While significant advances have been made in unraveling the cellular and molecular basis of lepidopteran pigmentary coloration, the morphogenesis of wing scale nanostructures involved in structural color production is not well understood. Contemporary research on this topic largely focuses on a few nymphalid model taxa (e.g., Bicyclus, Heliconius), despite an overwhelming diversity in the hierarchical nanostructural organization of lepidopteran wing scales. Here, we present a time-resolved, comparative developmental study of hierarchical scale nanostructures in Parides eurimedes and five other papilionid species. Our results uphold the putative conserved role of F-actin bundles in acting as spacers between developing ridges, as previously documented in several nymphalid species. Interestingly, while ridges are developing in P. eurimedes, plasma membrane manifests irregular mesh-like crossribs characteristic of Papilionidae, which delineate the accretion of cuticle into rows of planar disks in between ridges. Once the ridges have grown, disintegrating F-actin bundles appear to reorganize into a network that supports the invagination of plasma membrane underlying the disks, subsequently forming an extruded honeycomb lattice. Our results uncover a previously undocumented role for F-actin in the morphogenesis of complex wing scale nanostructures, likely specific to Papilionidae.

Convergence in sympatric swallowtail butterflies reveals ecological interactions as a key driver of worldwide trait diversification

Ecological interactions can promote phenotypic diversification in sympatric species. While competition can enhance trait divergence, other ecological interactions may promote convergence in sympatric species. Within butterflies, evolutionary convergences in wing color patterns have been reported between distantly related species, especially in females of palatable species, where mimetic color patterns are promoted by predator communities shared with defended species living in sympatry. Wing color patterns are also often involved in species recognition in butterflies, and divergence in this trait has been reported in closely related species living in sympatry as a result of reproductive character displacement. Here, we investigate the effect of sympatry between species on the convergence vs. divergence of their wing color patterns in relation to phylogenetic distance, focusing on the iconic swallowtail butterflies (family Papilionidae). We developed an unsupervised machine learning-based method to estimate phenotypic distances among wing color patterns of 337 species, enabling us to finely quantify morphological diversity at the global scale among species and allowing us to compute pairwise phenotypic distances between sympatric and allopatric species pairs. We found phenotypic convergence in sympatry, stronger among distantly related species, while divergence was weaker and restricted to closely related males. The convergence was stronger among females than males, suggesting that differential selective pressures acting on the two sexes drove sexual dimorphism. Our results highlight the significant effect of ecological interactions driven by predation pressures on trait diversification in Papilionidae and provide evidence for the interaction between phylogenetic proximity and ecological interactions in sympatry, acting on macroevolutionary patterns of phenotypic diversification.

Spatially differential regulation of ATF2 phosphorylation contributes to warning coloration of gregarious locusts

Warning coloration are common defense strategies used by animals to deter predators. Pestilential gregarious locusts exhibit a notable black-brown pattern as a form of warning coloration. However, the mechanisms regulating this distinctive pattern remain largely unknown. Here, we revealed that the black and brown integuments of locusts are governed by varying amounts of β-carotene and β-carotene-binding protein (βCBP) complexes. βCBP expression is regulated by the bZIP transcription factor activation transcription factor 2 (ATF2), which is activated by protein kinase C alpha in response to crowding. Specifically, ATF2 is phosphorylated at Ser327 and translocates to the nucleus, where it binds to the βCBP promoter and stimulates overexpression. Differential phosphorylation of ATF2 leads to the divergent black and brown coloration in gregarious locusts. The accumulation of red pigments vital for creating the brown sternum depends on βCBP overexpression. The spatial variation in ATF2 phosphorylation enables locusts to rapidly adapt to changing environment for aposematism.

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.

Genomic knockout of hsp23 both decreases and increases fitness under opposing thermal extremes in Drosophila melanogaster

Under exposure to harmful environmental stresses, organisms exhibit a general stress response involving upregulation of the expression of heat shock proteins (HSPs) which is thought to be adaptive. Small heat shock proteins (sHSPs) are key components of this response, although shsp genes may have other essential roles in development. However, the upregulation of expression of a suite of genes under stress may not necessarily be evidence of an adaptive response to stress that involves those genes. To explore this issue, we used the CRISPR/Cas9 system to investigate pleiotropic effects of the hsp23 gene in Drosophila melanogaster. Transgenic flies carrying a pCFD5 plasmid containing sgRNAs were created to generate a complete knockout of the hsp23 gene. The transgenic line lacking hsp23 showed an increased hatch rate and no major fitness costs under an intermediate temperature used for culturing the flies. In addition, hsp23 knockout affected tolerance to hot and cold temperature extremes but in opposing directions; knockout flies had reduced tolerance to cold, but increased tolerance to heat. Despite this, hsp23 expression (in wild type flies) was increased under both hot and cold conditions. The hsp23 gene was required for heat hardening at the pupal stage, but not at the 1st-instar larval stage, even though the gene was upregulated in wild type controls at that life stage. The phenotypic effects of hsp23 were not compensated for by expression changes in other shsps. Our study shows that the fitness consequences of an hsp gene knockout depends on environmental conditions, with potential fitness benefits of gene loss even under conditions when the gene is normally upregulated.