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

Establishing a Male-Positive Genetic Sexing Strain in the Asian Malaria Vector Anopheles stephensi

Genetic biocontrol interventions targeting mosquito-borne diseases require the release of male mosquitoes exclusively, as only females consume blood and transmit human pathogens. This reduces the risk of spreading pathogens while enabling effective population control. Robust sex sorting methods to enable early larval sorting in mosquitoes need to be developed to allow for scalable sex sorting for genetic biocontrol interventions. This study applies the SEPARATOR (Sexing Element Produced by Alternative RNA-splicing of A Transgenic Observable Reporter) system, previously developed for Aedes aegypti, to the Asian malaria vector Anopheles stephensi. We hypothesized that the intron from the doublesex gene in Anopheles gambiae would function in An. stephensi due to evolutionary conservation. Our results confirm that the splicing module from An. gambiae operates effectively in An. stephensi, demonstrating evolutionary conservation in sex-specific splicing events between these species. This system enables reliable positive male selection from first instar larval to pupal stages. RT-PCR analysis demonstrates that male-specific EGFP expression is dependent on doublesex sex-specific splicing events. The SEPARATOR system's independence from sex-chromosome linkage confers resistance to meiotic recombination and chromosomal rearrangements. This approach may facilitate the mass release of males, and the cross-species portability of SEPARATOR establishes it as a valuable tool for genetic biocontrol interventions across various pest species.

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.

Bivalves as Emerging Model Systems to Study the Mechanisms and Evolution of Sex Determination: A Genomic Point of View

Bivalves are a diverse group of molluscs that have recently attained a central role in plenty of biological research fields, thanks to their peculiar life history traits. Here, we propose that bivalves should be considered as emerging model systems also in sex-determination (SD) studies, since they would allow to investigate: 1) the transition between environmental and genetic SD, with respect to different reproductive backgrounds and sexual systems (from species with strict gonochorism to species with various forms of hermaphroditism); 2) the genomic evolution of sex chromosomes (SCs), considering that no heteromorphic SCs are currently known and that homomorphic SCs have been identified only in a few species of scallops; 3) the putative role of mitochondria at some level of the SD signaling pathway, in a mechanism that may resemble the cytoplasmatic male sterility of plants; 4) the evolutionary history of SD-related gene (SRG) families with respect to other animal groups. In particular, we think that this last topic may lay the foundations for expanding our understanding of bivalve SD, as our current knowledge is quite fragmented and limited to a few species. As a matter of fact, tracing the phylogenetic history and diversity of SRG families (such as the Dmrt, Sox, and Fox genes) would allow not only to perform more targeted functional experiments and genomic analyses, but also to foster the possibility of establishing a solid comparative framework.

The identification and expression pattern of the sex determination genes and their sex-specific variants in the egg parasitoid Trichogramma dendrolimi Matsumura (Hymenoptera: Trichogrammatidae)

Introduction: Trichogramma wasps are egg parasitoids of agricultural lepidopteran pests. The sex of Trichogramma is determined by its ploidy as well as certain sex ratio distorters, such as the endosymbiotic bacteria Wolbachia spp. and the paternal sex ratio (PSR) chromosome. The sex determination systems of hymenopterans, such as Trichogramma spp., involve cascades of the genes transformer (tra), transformer-2 (tra2), and doublesex (dsx) and are associated with sex-specific tra and dsx splicing. First, these genes and their sex-specific variants must be identified to elucidate the interactions between the sex ratio disorders and the sex determination mechanism of Trichogramma. Methods: Here, we characterized the sex determination genes tra, tra2, and dsx in Trichogramma dendrolimi. Sex-specific tra and dsx variants were detected in cDNA samples obtained from both male and female Trichogramma wasps. They were observed in the early embryos (1-10 h), late embryos (12-20 h), larvae (32 h and 48 h), pre-pupae (96 h), and pupae (144 h, 168 h, 192 h, and 216 h) of both male and female T. dendrolimi offspring. Results: We detected female-specific tra variants throughout the entire early female offspring stage. The male-specific variant began to express at 9-10 h as the egg was not fertilized. However, we did not find any maternally derived, female-specific tra variant in the early male embryo. This observation suggests that the female-specific tra variant expressed in the female embryo at 1-9 h may not have originated from the maternal female wasp. Discussion: The present study might be the first to identify the sex determination genes and sex-specific gene splicing in Trichogramma wasps. The findings of this study lay the foundation for investigating the sex determination mechanisms of Trichogramma and other wasps. They also facilitate sex identification in immature T. dendrolimi and the application of this important egg parasitoid in biological insect pest control programs.

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.

A history of the genetic and molecular identification of genes and their functions controlling insect sex determination

The genetics of the sex determination regulatory cascade in Drosophila melanogaster has a fascinating history, interlinked with the foundation of the Genetics discipline itself. The discovery that alternative splicing rather than differential transcription is the molecular mechanism underlying the upstream control of sex differences in the Drosophila model system was surprising. This notion is now fully integrated into the scientific canon, appearing in many genetics textbooks and online education resources. In the last three decades, it was a key reference point for starting evolutionary studies in other insect species by using homology-based approaches. This review will introduce a very brief history of Drosophila genetics. It will describe the genetic and molecular approaches applied for the identifying and cloning key genes involved in sex determination in Drosophila and in many other insect species. These comparative analyses led to supporting the idea that sex-determining pathways have evolved mainly by recruiting different upstream signals/genes while maintaining widely conserved intermediate and downstream regulatory genes. The review also provides examples of the link between technological advances and research achievements, to stimulate reflections on how science is produced. It aims to hopefully strengthen the related historical and conceptual knowledge of general readers of other disciplines and of younger geneticists, often focused on the latest technical-molecular approaches.

A male-specific doublesex isoform reveals an evolutionary pathway of sexual development via distinct alternative splicing mechanisms

The doublesex/mab-3 related transcription factor (Dmrt) genes regulate sexual development in metazoans. Studies of the doublesex (dsx) gene in insects, in particular Drosophila melanogaster, reveal that alternative splicing of dsx generates sex-specific Dsx isoforms underlying sexual differentiation. Such a splicing-based mechanism underlying sex-specific Dmrt function is thought to be evolved from a transcription-based mechanism used in non-insect species, but how such transition occurs during evolution is not known. Here we identified a male-specific dsx transcript (dsx^M2) through intron retention (IR), in addition to previously identified dsx^M and dsx^F transcripts through alternative polyadenylation (APA) with mutually exclusive exons. We found that Dsx^M2 had similarly masculinizing function as Dsx^M. We also found that the IR-based mechanism generating sex-specific dsx transcripts was conserved from flies to cockroaches. Further analysis of these dsx transcripts suggested an evolutionary pathway from sexually monomorphic to sex-specific dsx via the sequential use of IR-based and APA-based alternative splicing.

An ancestral male-specific doublesex isoform, dsxM2, is identified via intron retention, with a masculinizing function weaker than the modern dsxM

Genome analysis of Phrixothrix hirtus (Phengodidae) railroad worm shows the expansion of odorant-binding gene families and positive selection on morphogenesis and sex determination genes

Among bioluminescent beetles of the Elateroidea superfamily, Phengodidae is the third largest family, with 244 bioluminescent species distributed only in the Americas, but is still the least studied from the phylogenetic and evolutionary points of view. The railroad worm Phrixothrix hirtus is an essential biological model and symbolic species due to its bicolor bioluminescence, being the only organism that produces true red light among bioluminescent terrestrial species. Here, we performed partial genome assembly of P. hirtus, combining short and long reads generated with Illumina sequencing, providing the first source of genomic information and a framework for comparative analyses of the bioluminescent system in Elateroidea. This is the largest genome described in the Elateroidea superfamily, with an estimated size of ∼3.4 Gb, displaying 32 % GC content, and 67 % transposable elements. Comparative genomic analyses showed a positive selection of genes and gene family expansion events of growths and morphogenesis gene products, which could be associated with the atypical anatomical development and morphogenesis found in paedomorphic females and underdeveloped males. We also observed gene family expansion among distinct odorant-binding receptors, which could be associated with the pheromone communication system typical of these beetles, and retrotransposable elements. Common genes putatively regulating bioluminescence production and control, including two luciferase genes corresponding to lateral lanterns green-emitting and head lanterns red-emitting luciferases with 7 exons and 6 introns, and genes potentially involved in luciferin biosynthesis were found, indicating that there are no clear differences about the presence or absence of gene families associated with bioluminescence in Elateroidea.

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.

Evolutionary history of sexual differentiation mechanism in insects

Alternative splicing underpins functional diversity in proteins and the complexity and diversity of eukaryotes. An example is the doublesex gene, the key transcriptional factor in arthropod sexual differentiation. doublesex is controlled by sex-specific splicing and promotes both male and female differentiation in holometabolan insects, whereas in hemimetabolan species, doublesex has sex-specific isoforms but is not required for female differentiation. How doublesex evolved to be essential for female development remains largely unknown. Here, we investigate ancestral states of doublesex using Thermobia domestica belonging to Zygentoma, the sister group of Pterygota, that is, winged insects. We find that, in T. domestica, doublesex expresses sex-specific isoforms but is only necessary for male differentiation of sexual morphology. This result supports the hypothesis that doublesex initially promoted male differentiation during insect evolution. However, T. domestica doublesex has a short female-specific region and upregulates the expression of vitellogenin homologs in females, suggesting that doublesex may already play some role in female morphogenesis of the common ancestor of Pterygota. Reconstruction of the ancestral sequence and prediction of protein structures show that the female-specific isoform of doublesex has an extended C-terminal disordered region in holometabolan insects but not in nonholometabolan species. We propose that doublesex acquired its function in female morphogenesis through a change in the protein motif structure rather than the emergence of the female-specific exon.

Evolutionary transition of doublesex regulation from sex-specific splicing to male-specific transcription in termites

The sex determination gene doublesex (dsx) encodes a transcription factor with two domains, oligomerization domain 1 (OD1) and OD2, and is present throughout insects. Sex-specific Dsx splicing isoforms regulate the transcription of target genes and trigger sex differentiation in all Holometabola examined to date. However, in some hemimetabolous insects, dsx is not spliced sexually and its sequence is less conserved. Here, to elucidate evolutionary changes in dsx in domain organisation and regulation in termites, we searched genome and/or transcriptome databases for the dsx OD1 and OD2 in seven termite species and their sister group (Cryptocercus woodroaches). Molecular phylogenetic and synteny analyses identified OD1 sequences of termites and C. punctulatus that clustered with dsx of Holometabola and regarded them as dsx orthologues. The Cryptocercus dsx orthologue containing OD2 was spliced sexually, as previously shown in other insects. However, OD2 was not found in all termite dsx orthologues. These orthologues were encoded by a single exon in three termites for which genome information is available; they were not alternatively spliced but transcribed in a male-specific manner in two examined species. Evolution of dsx regulation from sex-specific splicing to male-specific transcription may have occurred at an early stage of social evolution in termites.

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.

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.

Genome assembly, sex-biased gene expression and dosage compensation in the damselfly Ischnura elegans

The evolution of sex chromosomes, and patterns of sex-biased gene expression and dosage compensation, are poorly known among early winged insects such as odonates. We assembled and annotated the genome of Ischnura elegans (blue-tailed damselfly), which, like other odonates, has a male-hemigametic sex-determining system (X0 males, XX females). By identifying X-linked genes in I. elegans and their orthologs in other insect genomes, we found homologies between the X chromosome in odonates and chromosomes of other orders, including the X chromosome in Coleoptera. Next, we showed balanced expression of X-linked genes between sexes in adult I. elegans, i.e. evidence of dosage compensation. Finally, among the genes in the sex-determining pathway only fruitless was found to be X-linked, while only doublesex showed sex-biased expression. This study reveals partly conserved sex chromosome synteny and independent evolution of dosage compensation among insect orders separated by several hundred million years of evolutionary history.

The direct regulation of Aalbdsx on AalVgR is indispensable for ovarian development in Aedes albopictus

BACKGROUND

Aedes albopictus is an important vector with an extensive worldwide distribution. Only female mosquitoes play a significant role in the transmission of pathogens. Doublesex (dsx) is a central nexus gene in the insect somatic sex determination hierarchy.

RESULTS

In this study, we characterized the full-length sex-specific splicing forms of the Ae. albopictus dsx (Aalbdsx) gene. Then, we identified 15 direct target genes of DSX in adult females using digital gene expression combined with quantitative real-time polymerase chain reaction (qPCR) by performing a chromatin immunoprecipitation (ChIP) assay with specific DSX antibodies. Knockdown of Aalbdsx suppressed ovarian development and decreased the transcript levels of the Aalbdsx target vitellogenin receptor (VgR) gene, whereas vitellogenin (Vg) expression showed an increase in the fat body. Genes in the major Vg regulatory pathway were also up-regulated.

CONCLUSION

Our results suggest that both Vg and VgR are direct target genes of Aalbdsx and that direct regulation of Aalbdsx on VgR is indispensable for ovarian development in Ae. albopictus, which not only provides a reference for the further elucidation of the evolutionarily conserved role of dsx in Ae. albopictus sexual differentiation but also reveals potential molecular targets for application to the development of sterile male mosquitoes to be released for vector control. © 2020 Society of Chemical Industry.

Doublesex Mediates the Development of Sex-Specific Pheromone Organs in Bicyclus Butterflies via Multiple Mechanisms

The Bicyclus lineage of satyrid butterflies exhibits male-specific traits, the scent organs, used for chemical communication during courtship. These organs consist of tightly packed brush-like scales (hair-pencils) that rub against scent patches to disperse pheromones, but the evolution and molecular basis of these organ’s male-limited development remains unknown. Here, we examine the evolution of the number and location of the scent patches and hair-pencils within 53 species of Bicyclus butterflies, and the involvement of the sex determinant gene doublesex (dsx) in scent organ development in Bicyclus anynana using CRISPR/Cas9. We show that scent patches and hair-pencils arose via multiple, independent gains, in a correlated manner. Further, an initially nonsex-specific Dsx protein expression pattern in developing wing discs becomes male-specific and spatially refined to areas that develop the scent patches. Functional perturbations of dsx show that this gene activates patch development in males whereas hair-pencils develop in both sexes without Dsx input. Dsx in females is, instead, required to repress hair-pencils whereas Dsx in males regulates minor aspects of its development. These findings suggest that the patches and hair-pencils evolve as correlated composite organs presumably due to their functional integration. Divergence in the function of dsx isoforms occurred across the sexes, where the male isoform promotes patch development in males and the female isoform represses hair-pencil development in females, both leading to the development of male-limited traits. Furthermore, evolution in number of patches in males is due to the evolution of spatial regulation of dsx.

The Sex Determination Cascade in the Silkworm

In insects, sex determination pathways involve three levels of master regulators: primary signals, which determine the sex; executors, which control sex-specific differentiation of tissues and organs; and transducers, which link the primary signals to the executors. The primary signals differ widely among insect species. In Diptera alone, several unrelated primary sex determiners have been identified. However, the doublesex (dsx) gene is highly conserved as the executor component across multiple insect orders. The transducer level shows an intermediate level of conservation. In many, but not all examined insects, a key transducer role is performed by transformer (tra), which controls sex-specific splicing of dsx. In Lepidoptera, studies of sex determination have focused on the lepidopteran model species Bombyx mori (the silkworm). In B. mori, the primary signal of sex determination cascade starts from Fem, a female-specific PIWI-interacting RNA, and its targeting gene Masc, which is apparently specific to and conserved among Lepidoptera. Tra has not been found in Lepidoptera. Instead, the B. mori PSI protein binds directly to dsx pre-mRNA and regulates its alternative splicing to produce male- and female-specific transcripts. Despite this basic understanding of the molecular mechanisms underlying sex determination, the links among the primary signals, transducers and executors remain largely unknown in Lepidoptera. In this review, we focus on the latest findings regarding the functions and working mechanisms of genes involved in feminization and masculinization in Lepidoptera and discuss directions for future research of sex determination in the silkworm.

Genomics Reveals Widespread Ecological Speciation in Flightless Insects

Recent genomic analyses have highlighted parallel divergence in response to ecological gradients, but the extent to which altitude can underpin such repeated speciation remains unclear. Wing reduction and flight loss have apparently evolved repeatedly in montane insect assemblages and have been suggested as important drivers of hexapod diversification. We test this hypothesis using genomic analyses of a widespread wing-polymorphic stonefly species complex in New Zealand. We identified over 50,000 polymorphic genetic markers generated across almost 200 Zelandoperla fenestrata stonefly specimens using a newly generated plecopteran reference genome, to reveal widespread parallel speciation between sympatric full-winged and wing-reduced ecotypes. Rather than the existence of a single, widespread, flightless taxon (Zelandoperla pennulata), evolutionary genomic data reveal that wing-reduced upland lineages have speciated repeatedly and independently from full-winged Z. fenestrata. This repeated evolution of reproductive isolation between local ecotype pairs that lack mitochondrial DNA differentiation suggests that ecological speciation has evolved recently. A cluster of outlier single-nucleotide polymorphisms detected in independently wing-reduced lineages, tightly linked in an approximately 85 kb genomic region that includes the developmental “supergene” doublesex, suggests that this “island of divergence” may play a key role in rapid ecological speciation. [Ecological speciation; genome assembly; genomic island of differentiation; genotyping-by-sequencing; incipient species; plecoptera; wing reduction.]

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.

Copulatory function and development shape modular architecture of genitalia differently in males and females

Genitalia are multitasking structures whose development is mediated by numerous regulatory pathways. This multifactorial nature provides an avenue for multiple sources of selection. As a result, genitalia tend to evolve as modular systems comprising semi-independent subsets of structures, yet the processes that give rise to those patterns are still poorly understood. Here, we ask what are the relative roles of development and function in shaping modular patterns of genitalia within populations and across species of stink-bugs. We found that male genitalia are less integrated, more modular, and primarily shaped by functional demands. In contrast, females show higher integration, lower modularity, and a predominant role of developmental processes. Further, interactions among parts of each sex are more determinant to modularity than those between the sexes, and patterns of modularity are equivalent between and within species. Our results strongly indicate that genitalia have been subjected to sex-specific selection, although male and female genitalia are homologous and functionally associated. Moreover, modular patterns are seemingly constant in the evolutionary history of stink-bugs, suggesting a scenario of multivariate stabilizing selection within each sex. Our study demonstrates that interactions among genital parts of the same sex may be more fundamental to genital evolution than previously thought.