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

Doublesex is essential for masculinization but not feminization in Lygus hesperus

In most holometabolous insects, sex differentiation occurs via a hierarchical cascade of transcription factors, with doublesex (dsx) regulating genes that control sex-specific traits. Although less is known in hemimetabolous insects, early evidence suggests that substantial differences exist from more evolutionarily advanced insects. Here, we identified and characterized dsx in Lygus hesperus (western tarnished plant bug), a hemipteran pest of many agricultural crops in western North America. The full-length transcript for L. hesperus dsx (Lhdsx) and several variants encode proteins with conserved DNA binding and oligomerization domains. Transcript profiling revealed that Lhdsx is ubiquitously expressed, likely undergoes alternative pre-mRNA splicing, and, unlike several model insects, is sex-biased rather than sex-specific. Embryonic RNA interference (RNAi) of Lhdsx only impacted sex development in adult males, which lacked both internal reproductive organs and external genitalia. No discernible impacts on adult female development or reproductivity were observed. RNAi knockdown of Lhdsx in nymphs likewise only affected adult males, which lacked the characteristic dimorphic coloration but had dramatically elevated vitellogenin transcripts. Gene knockout of Lhdsx by CRISPR/Cas9 editing yielded only females in G0 and strongly biased heterozygous G1 offspring to females with the few surviving males showing severely impaired genital development. These results indicate that L. hesperus male development requires Lhdsx, whereas female development proceeds via a basal pathway that functions independently of dsx. A fundamental understanding of sex differentiation in L. hesperus could be important for future gene-based management strategies of this important agricultural pest.

Transcriptional regulation of IAG by dsx and foxl-2 in mud crab (Scylla paramamosain)

Scylla paramamosain is an important cultured crab species on the southeast coast of China. However, the molecular regulation mechanism of its gonadal development still has not been thoroughly studied. Dsx (doublesex) and foxl-2 (forkhead transcription factor gene 2) are important transcription factors involved in gonadal development. So far, studies on the functions of dsx and foxl-2 in crustaceans are very limited. Insulin-like androgenic gland hormone (IAG) is an effector molecule that regulates the differentiation, development and sex maintenance of testes in crustaceans. In this study, the promoter region of Sp-IAG was predicted, and several potential binding sites of dsx and foxl-2 were found. Site-directed mutagenesis was performed on the predicted potential binding sites, and their promoter activity was analyzed. The results showed that there was a dsx and a foxl-2 binding site, respectively, that could regulate the expression of Sp-IAG. In order to verify the regulatory effect of these two transcription factors on Sp-IAG, we constructed the expression plasmids of dsx and foxl-2 and co-transfected them into HEK293T cell lines with the promoter of Sp-IAG, respectively. The results showed that dsx could significantly promote the expression of Sp-IAG, while foxl-2 could inhibit its expression substantially. Then we carried out in vivo RNA interference experiment on mud crabs. The expression of dsx and foxl-2 in crabs was interfered respectively. The results of qRT-PCR showed that the expression of Sp-IAG was significantly inhibited after interfering with dsx, while significantly increased after interfering with foxl-2, which was consistent with the cell experiment. In conclusion, dsx and foxl-2 transcription factors play opposite roles in regulating the expression of Sp-IAG.

Acute and Long-Term Consequences of Co-opted doublesex on the Development of Mimetic Butterfly Color Patterns

Novel phenotypes are increasingly recognized to have evolved by co-option of conserved genes into new developmental contexts, yet the process by which co-opted genes modify existing developmental programs remains obscure. Here, we provide insight into this process by characterizing the role of co-opted doublesex in butterfly wing color pattern development. dsx is the master regulator of insect sex differentiation but has been co-opted to control the switch between discrete nonmimetic and mimetic patterns in Papilio alphenor and its relatives through the evolution of novel mimetic alleles. We found dynamic spatial and temporal expression pattern differences between mimetic and nonmimetic butterflies throughout wing development. A mimetic color pattern program is switched on by a pulse of dsx expression in early pupal development that causes acute and long-term differential gene expression, particularly in Wnt and Hedgehog signaling pathways. RNAi suggested opposing, novel roles for these pathways in mimetic pattern development. Importantly, Dsx co-option caused Engrailed, a primary target of Hedgehog signaling, to gain a novel expression domain early in pupal wing development that is propagated through mid-pupal development to specify novel mimetic patterns despite becoming decoupled from Dsx expression itself. Altogether, our findings provide multiple views into how co-opted genes can both cause and elicit changes to conserved networks and pathways to result in development of novel, adaptive phenotypes.

Evidence for a single, ancient origin of a genus-wide alternative life history strategy

Understanding the evolutionary origins and factors maintaining alternative life history strategies (ALHS) within species is a major goal of evolutionary research. While alternative alleles causing discrete ALHS are expected to purge or fix over time, one-third of the ~90 species of Colias butterflies are polymorphic for a female-limited ALHS called Alba. Whether Alba arose once, evolved in parallel, or has been exchanged among taxa is currently unknown. Using comparative genome-wide association study (GWAS) and population genomic analyses, we placed the genetic basis of Alba in time-calibrated phylogenomic framework, revealing that Alba evolved once near the base of the genus and has been subsequently maintained via introgression and balancing selection. CRISPR-Cas9 mutagenesis was then used to verify a putative cis-regulatory region of Alba, which we identified using phylogenetic foot printing. We hypothesize that this cis-regulatory region acts as a modular enhancer for the induction of the Alba ALHS, which has likely facilitated its long evolutionary persistence.

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.

Distinct developmental mechanisms influence sexual dimorphisms in the milkweed bug Oncopeltus fasciatus

Sexual dimorphism is common in animals. The most complete model of sex determination comes from Drosophila melanogaster, where the relative dosage of autosomes and X chromosomes leads indirectly to sex-specific transcripts of doublesex (dsx). Female Dsx interacts with a mediator complex protein encoded by intersex (ix) to activate female development. In males, the transcription factor encoded by fruitless (fru) promotes male-specific behaviour. The genetics of sex determination have been examined in a small number of other insects, yet several questions remain about the plesiomorphic state. Is dsx required for female and male development? Is fru conserved in male behaviour or morphology? Are other components such as ix functionally conserved? To address these questions, we report expression and functional tests of dsx, ix and fru in the hemipteran Oncopeltus fasciatus, characterizing three sexual dimorphisms. dsx prevents ix phenotypes in all sexes and dimorphic traits in the milkweed bug. ix and fru are expressed across the body, in females and males. fru and ix also affect the genitalia of both sexes, but have effects limited to different dimorphic structures in different sexes. These results reveal roles for ix and fru distinct from other insects, and demonstrate distinct development mechanisms in different sexually dimorphic structures.

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

Phenotypic Plasticity of the Mimetic Swallowtail Butterfly Papilio polytes: Color Pattern Modifications and Their Implications in Mimicry Evolution

Simple Summary

Diverse butterfly wing color patterns are evolutionary products in response to environmental changes in the past. Environmental stress, such as temperature shock, is known to induce color pattern modifications in various butterfly species, and this phenotypic plasticity plays an important role in color pattern evolution. However, the potential contributions of phenotypic plasticity to mimicry evolution have not been evaluated. Here, we focused on the swallowtail butterfly Papilio polytes, which has nonmimetic and mimetic forms in females, to examine its plastic phenotypes. Cold shock and heat shock treatments in the nonmimetic form induced color pattern modifications that were partly similar to those of the mimetic form, and nonmimetic females were more sensitive than males and mimetic females. These results suggest that phenotypic plasticity in nonmimetic females might have provided a basis of natural selection for mimetic color patterns during evolution.

Abstract

Butterfly wing color patterns are sensitive to environmental stress, such as temperature shock, and this phenotypic plasticity plays an important role in color pattern evolution. However, the potential contributions of phenotypic plasticity to mimicry evolution have not been evaluated. Here, we focused on the swallowtail butterfly Papilio polytes, which has nonmimetic and mimetic forms in females, to examine its plastic phenotypes. In the nonmimetic form, medial white spots and submarginal reddish spots in the ventral hindwings were enlarged by cold shock but were mostly reduced in size by heat shock. These temperature-shock-induced color pattern modifications were partly similar to mimetic color patterns, and nonmimetic females were more sensitive than males and mimetic females. Unexpectedly, injection of tungstate, a known modification inducer in nymphalid and lycaenid butterflies, did not induce any modification, but fluorescent brightener 28, another inducer discovered recently, induced unique modifications. These results suggest that phenotypic plasticity in nonmimetic females might have provided a basis of natural selection for mimetic color patterns during evolution.

Alternative splicing in seasonal plasticity and the potential for adaptation to environmental change

Seasonal plasticity is accomplished via tightly regulated developmental cascades that translate environmental cues into trait changes. Little is known about how alternative splicing and other posttranscriptional molecular mechanisms contribute to plasticity or how these mechanisms impact how plasticity evolves. Here, we use transcriptomic and genomic data from the butterfly Bicyclus anynana, a model system for seasonal plasticity, to compare the extent of differential expression and splicing and test how these axes of transcriptional plasticity differ in their potential for evolutionary change. Between seasonal morphs, we find that differential splicing affects a smaller but functionally unique set of genes compared to differential expression. Further, we find strong support for the novel hypothesis that spliced genes are more susceptible than differentially expressed genes to erosion of genetic variation due to selection on seasonal plasticity. Our results suggest that splicing plasticity is especially likely to experience genetic constraints that could affect the potential of wild populations to respond to rapidly changing environments.

Novel doublesex duplication associated with sexually dimorphic development of dogface butterfly wings

Sexually dimorphic development is responsible for some of the most remarkable phenotypic variation found in nature. Alternative splicing of the transcription factor gene doublesex (dsx) is a highly conserved developmental switch controlling the expression of sex-specific pathways. Here, we leverage sex-specific differences in butterfly wing color pattern to characterize the genetic basis of sexually dimorphic development. We use RNA-seq, immunolocalization, and motif binding site analysis to test specific predictions about the role of dsx in the development of structurally based ultraviolet (UV) wing patterns in Zerene cesonia (Southern Dogface). Unexpectedly, we discover a novel duplication of dsx that shows a sex-specific burst of expression associated with the sexually dimorphic UV coloration. The derived copy consists of a single exon that encodes a DNA binding but no protein-binding domain and has experienced rapid amino-acid divergence. We propose the novel dsx paralog may suppress UV scale differentiation in females, which is supported by an excess of Dsx-binding sites at cytoskeletal and chitin-related genes with sex-biased expression. These findings illustrate the molecular flexibility of the dsx gene in mediating the differentiation of secondary sexual characteristics.