Evidence for the evolutionary nascence of a novel sex determination pathway in honeybees

Understanding the genetics of sex determination in insects and its relevance to genetic pest management

Sex determination pathways regulate male and female-specific development and differentiation and offer potential targets for genetic pest management methods. Insect sex determination pathways are comprised of primary signals, relay genes and terminal genes. Primary signals of coleopteran, dipteran, hymenopteran and lepidopteran species are highly diverse and regulate the sex-specific splicing of relay genes based on the primary signal dosage, amino acid composition or the interaction with paternally inherited genes. In coleopterans, hymenopterans and some dipterans, relay genes are Transformer orthologs from the serine-arginine protein family that regulate sex-specific splicing of the terminal genes. Alternative genes regulate the splicing of the terminal genes in dipterans that lack Transformer orthologs and lepidopterans. Doublesex and Fruitless orthologs are the terminal genes. Doublesex and Fruitless orthologs are highly conserved zinc-finger proteins that regulate the expression of downstream proteins influencing physical traits and courtship behaviours in a sex-specific manner. Genetic pest management methods can use different mechanisms to exploit or disrupt female-specific regions of different sex determination genes. Female-specific regions of sex determination genes can be exploited to produce a lethal gene only in females or disrupted to impede female development or fertility. Reducing the number of fertile females in pest populations creates a male-biased sex ratio and eventually leads to the local elimination of the pest population. Knowledge on the genetic basis of sex determination is important to enable these sex determination pathways to be exploited for genetic pest management.

Heterozygosity at a conserved candidate sex determination locus is associated with female development in the clonal raider ant (Ooceraea biroi)

Sex determination is a developmental switch that triggers sex-specific developmental programs. This switch is "flipped" by the expression of genes that promote male- or female-specific development. Many lineages have evolved sex chromosomes that act as primary signals for sex determination. However, haplodiploidy (males are haploid and females are diploid), which occurs in ca. 12% of animal species, is incompatible with sex chromosomes. Haplodiploid taxa must, therefore, rely on other strategies for sex determination. One mechanism, "complementary sex determination" (CSD), uses heterozygosity as a proxy for diploidy. In CSD, heterozygosity at a sex determination locus triggers female development, while hemizygosity or homozygosity permits male development. CSD loci have been mapped in honeybees and two ant species, but we know little about their evolutionary history. Here, we investigate sex determination in the clonal raider ant, Ooceraea biroi. We identified a 46kb candidate CSD locus at which all females are heterozygous, but most diploid males are homozygous for either allele. As expected for CSD loci, the candidate locus has more alleles than most other loci, resulting in a peak of nucleotide diversity. This peak negligibly affects the amino acid sequences of protein-coding genes, suggesting that heterozygosity of a non-coding genomic sequence triggers female development. This locus is distinct from the CSD locus in honeybees but homologous to a CSD locus mapped in two distantly related ant species, implying that this molecular mechanism has been conserved since a common ancestor that lived approximately 112 million years ago.

Alternative double strand break repair pathways shape the evolution of high recombination in the honey bee, Apis mellifera

Social insects, particularly honey bees, have exceptionally high genomic frequencies of genetic recombination. This phenomenon and underlying mechanisms are poorly understood. To characterise the patterns of crossovers and gene conversion in the honey bee genome, a recombination map of 187 honey bee brothers was generated by whole-genome resequencing. Recombination events were heterogeneously distributed without many true hotspots. The tract lengths between phase shifts were bimodally distributed, indicating distinct crossover and gene conversion events. While crossovers predominantly occurred in G/C-rich regions and seemed to cause G/C enrichment, the gene conversions were found predominantly in A/T-rich regions. The nucleotide composition of sequences involved in gene conversions that were associated with or distant from crossovers corresponded to the differences between crossovers and gene conversions. These combined results suggest two types of DNA double-strand break repair during honey bee meiosis: non-canonical homologous recombination, leading to gene conversion and A/T enrichment of the genome, and the canonical homologous recombination based on completed double Holliday Junctions, which can result in gene conversion or crossover and is associated with G/C bias. This G/C bias may be selected for to balance the A/T-rich base composition of eusocial hymenopteran genomes. The lack of evidence for a preference of the canonical homologous recombination for double-strand break repair suggests that the high genomic recombination rate of honey bees is mainly the consequence of a high rate of double-strand breaks, which could in turn result from the life history of honey bees and their A/T-rich genome.

Dedicated developmental programing for group-supporting behaviors in eusocial honeybees

The evolutionary changes from solitary to eusocial living in vertebrates and invertebrates are associated with the diversification of social interactions and the development of queen and worker castes. Despite strong innate patterns, our understanding of the mechanisms manifesting these sophisticated behaviors is still rudimentary. Here, we show that doublesex (dsx) manifests group-supporting behaviors in the honeybee (Apis mellifera) worker caste. Computer-based individual behavioral tracking of worker bees with biallelic stop mutations in colonies revealed that the dsx gene is required for the rate and duration of group-supporting behavior that scales the relationship between bees and their work. General sensorimotor functions remained unaffected. Unexpectedly, unlike in other insects, the dsx gene is required for the neuronal wiring of the mushroom body in which the gene is spatially restricted expressed. Together, our study establishes dedicated programming for group-supporting behaviors and provides insight into the connection between development in the neuronal circuitry and behaviors regulating the formation of a eusocial society.

Honeybees' novel complementary sex-determining system: function and origin

Complementary sex determination regulates female and male development in honeybees (Apis mellifera) via heterozygous versus homo-/hemizygous genotypes of the csd (complementary sex determiner) gene involving numerous naturally occurring alleles. This lineage-specific function offers a rare opportunity to understand an undescribed regulatory mechanism and the molecular evolutionary path leading to this mechanism. We reviewed recent advances in understanding how Csd recognizes different versus identical protein variants, how these variants regulate downstream pathways and sexual differentiation, and how this mechanism has evolved and been shaped by evolutionary forces. Finally, we highlighted the shared regulatory principles of sex determination despite the diversity of primary signals and demonstrated that lineage-specific mutations are very informative for characterizing newly evolved functions.

LncRNA gene ANTSR coordinates complementary sex determination in the Argentine ant

Animals have evolved various sex determination systems. Here, we describe a newly found mechanism. A long noncoding RNA (lncRNA) transduces complementary sex determination (CSD) signal in the invasive Argentine ant. In this haplodiploid species, we identified a 5-kilobase hyper-polymorphic region underlying CSD: Heterozygous embryos become females, while homozygous and hemizygous embryos become males. Heterozygosity at the CSD locus correlates with higher expression of ANTSR, a gene that overlaps with the CSD locus and specifies an lncRNA transcript. ANTSR knockdown in CSD heterozygotes leads to male development. Comparative analyses indicated that, in Hymenoptera, ANTSR is an ancient yet rapidly evolving gene. This study reveals an lncRNA involved in genetic sex determination, alongside a previously unknown regulatory mechanism underlying sex determination based on complementarity among noncoding alleles.

Identification and functional analysis of the female determiner gene in the bean bug, Riptortus pedestris

BACKGROUND

Homing-based gene drives targeting sex-specific lethal genes have been used for genetic control. Additionally, understanding insect sex determination provides new targets for managing insect pests. While sex determination mechanisms in holometabolous insects have been thoroughly studied and employed in pest control, the study of the sex determination pathway in hemimetabolous insects is limited to only a few species. Riptortus pedestris (Fabricius; Hemiptera: Heteroptera), commonly known as the bean bug, is a significant pest for soybeans. Nonetheless, the mechanism of its sex determination and the target gene for genetic control are not well understood.

RESULTS

We identified Rpfmd as the female determiner gene in the sex determination pathway of R. pedestris. Rpfmd encodes a female-specific serine/arginine-rich protein of 436 amino acids and one non-sex-specific short protein of 98 amino acids. Knockdown of Rpfmd in R. pedestris nymphs caused death of molting females with masculinized somatic morphology but did not affect male development. Knockdown of Rpfmd in newly emerged females inhibited ovary development, while maternal-mediated RNA interference (RNAi) knockdown of Rpfmd expression resulted in male-only offspring. Transcriptome sequencing revealed that Rpfmd regulates X chromosome dosage compensation and influences various biological processes in females but has no significant effect on males. Moreover, RNAi mediated knockdown of Rpfmd-C had no influence on the development of R. pedestris, suggesting that Rpfmd regulates sex determination through female-specific splicing isoforms. We also found that Rpfmd pre-mRNA alternative splicing regulation starts at the 24-h embryo stage, indicating the activation of sex differentiation.

CONCLUSION

Our study confirms that Rpfmd, particularly its female-specific isoform (Rpfmd-F), is the female determiner gene that regulates sex differentiation in R. pedestris. Knockdown of Rpfmd results in female-specific lethality without affecting males, making it a promising target for genetic control of this soybean pest throughout its development stages. Additionally, our findings improve the understanding of the sex-determination mechanism in hemimetabolous insects. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Recognition of polymorphic Csd proteins determines sex in the honeybee

Sex in honeybees, Apis mellifera, is genetically determined by heterozygous versus homo/hemizygous genotypes involving numerous alleles at the single complementary sex determination locus. The molecular mechanism of sex determination is however unknown because there are more than 4950 known possible allele combinations, but only two sexes in the species. We show how protein variants expressed from complementary sex determiner (csd) gene determine sex. In females, the amino acid differences between Csd variants at the potential-specifying domain (PSD) direct the selection of a conserved coiled-coil domain for binding and protein complexation. This recognition mechanism activates Csd proteins and, thus, the female pathway. In males, the absence of polymorphisms establishes other binding elements at PSD for binding and complexation of identical Csd proteins. This second recognition mechanism inactivates Csd proteins and commits male development via default pathway. Our results demonstrate that the recognition of different versus identical variants of a single protein is a mechanism to determine sex.

Masculinizer gene controls male sex determination in the codling moth, Cydia pomonella

The molecular mechanisms of sex determination in moths and butterflies (Lepidoptera) with female heterogamety (WZ/ZZ) are poorly understood, except in the silkworm Bombyx mori. However, the Masculinizer (Masc) gene that controls male development and dosage compensation in B. mori, appears to be conserved in Lepidoptera, as its masculinizing function was recently confirmed in several moth species. In this work, we investigated the role of the Masc gene in sex determination of the codling moth Cydia pomonella (Tortricidae), a globally important pest of pome fruits and walnuts. The gene structure of the C. pomonella Masc ortholog, CpMasc, is similar to B. mori Masc. However, unlike B. mori, we identified 14 splice variants of CpMasc in the available transcriptomes. Subsequent screening for sex specificity and genetic variation using publicly available data and RT-PCR revealed three male-specific splice variants. Then qPCR analysis of these variants revealed sex-biased expression showing a peak only in early male embryos. Knockdown of CpMasc by RNAi during early embryogenesis resulted in a shift from male-to female-specific splicing of the C. pomonella doublesex (Cpdsx) gene, its downstream effector, in ZZ embryos, leading to a strongly female-biased sex ratio. These data clearly demonstrate that CpMasc functions as a masculinizing gene in the sex-determining cascade of C. pomonella. Our study also showed that CpMasc transcripts are provided maternally, as they were detected in unfertilized eggs after oviposition and in mature eggs dissected from virgin females. This finding is unique, as maternal provision of mRNA has rarely been studied in Lepidoptera.

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.

Intersex is required for female sexual development in Hermetia illucens

Sex-determination pathways are extremely diverse. Understanding the mechanism of sex determination in insects is important for genetic manipulation of the pest population and for breeding of economically valuable insects. Although sex determination has been well characterized in the model species Drosophila melanogaster, little is known about this pathway in Stratiomyidae. In the present study, we first identified the Drosophila intersex (ix) homolog in Hermetia illucens, also known as the black soldier fly, which belongs to the Stratiomyidae family and which is an important insect for the conversion of various organic wastes. Phylogenetic analyses and multiple sequence alignment revealed that Hiix is conserved compared with Drosophila. We showed that Hiix is highly expressed in internal genitalia. Disruption of the Hiix gene using CRISPR/Cas9 resulted in female-specific defects in external genitalia and abnormal and undersized ovaries. Taken together, our study furthers our understanding of sex determination in insects and could facilitate breeding of H. illucens.

Transformer gene regulates feminization under two complementary sex determination loci in the ant, Vollenhovia emeryi

Organisms that reproduce sexually have evolved well-organized mechanisms to determine two sexes. Some hymenopterans (such as ants, bees, and wasps) have a complementary sex-determination system in which heterozygosity at one CSD locus induces female development, whereas hemi- or homozygosity at the locus induces male development. This system can generate a high cost of inbreeding, as individuals that are homozygous at the locus become sterile, diploid males. On the other hand, some hymenopterans have evolved a multi-locus, complementary, sex-determination system in which heterozygosity in at least one CSD locus induces female development. This system effectively reduces the proportion of sterile diploid males; however, how these multiple
primary signals based on CSD pass through a molecular cascade to regulate downstream genes has remained unclear. To clarify this matter, we used a backcross to investigate the molecular cascade in the ant, Vollenhovia emeryi, with two CSD loci. Here we show by gene disruption that transformer (tra) is necessary for proper feminization. Expression analysis of tra and doublesex (dsx) showed that heterozygosity in at least one of the two CSD loci is sufficient to promote female sex determination. Analysis of overexpression suggested that female-type Tra protein promotes splicing of tra pre-mRNA to female isoform by a positive-regulatory-feedback loop. Our data also showed that tra affects splicing of dsx. We conclude that two-loci sex determination system in V. emeryi evolved based on tra-dsx splicing cascade that is well conserved in other insect species. Finally, we suggest a cascade model to arrive at a binary determination of sex under multiple primary signals.

In Vitro Comparison of Sex-Specific Splicing Efficiencies of fem Pre-mRNA under Monoallelic and Heteroallelic Conditions of csd, a Master Sex-Determining Gene in the Honeybee

The sexual fate of honeybees is determined by the complementary sex determination (CSD) model: heterozygosity at a single locus (the CSD locus) determines femaleness, while hemizygosity or homozygosity at the CSD locus determines maleness. The csd gene encodes a splicing factor that regulates sex-specific splicing of the downstream target gene feminizer (fem), which is required for femaleness. The female mode of fem splicing occurs only when csd is present in the heteroallelic condition. To gain insights into how Csd proteins are only activated under the heterozygous allelic composition, we developed an in vitro assay system to evaluate the activity of Csd proteins. Consistent with the CSD model, the co-expression of two csd alleles, both of which lack splicing activity under the single-allele condition, restored the splicing activity that governs the female mode of fem splicing. RNA immunoprecipitation quantitative PCR analyses demonstrated that the CSD protein was specifically enriched in several exonic regions in the fem pre-mRNA, and enrichment in exons 3a and 5 was significantly greater under the heterozygous allelic composition than the single-allelic condition. However, in most cases csd expression under the monoallelic condition was capable of inducing the female mode of fem splicing contrary to the conventional CSD model. In contrast, repression of the male mode of fem splicing was predominant under heteroallelic conditions. These results were reproduced by real-time PCR of endogenous fem expression in female and male pupae. These findings strongly suggest that the heteroallelic composition of csd may be more important for the repression of the male splicing mode than for the induction of the female splicing mode of the fem gene.

Theme and variation in the evolution of insect sex determination

The development of dimorphic adult sexes is a critical process for most animals, one that is subject to intense selection. Work in vertebrate and insect model species has revealed that sex determination mechanisms vary widely among animal groups. However, this variation is not uniform, with a limited number of conserved factors. Therefore, sex determination offers an excellent context to consider themes and variations in gene network evolution. Here we review the literature describing sex determination in diverse insects. We have screened public genomic sequence databases for orthologs and duplicates of 25 genes involved in insect sex determination, identifying patterns of presence and absence. These genes and a 3.5 reference set of 43 others were used to infer phylogenies and compared to accepted organismal relationships to examine patterns of congruence and divergence. The function of candidate genes for roles in sex determination (virilizer, female-lethal-2-d, transformer-2) and sex chromosome dosage compensation (male specific lethal-1, msl-2, msl-3) were tested using RNA interference in the milkweed bug, Oncopeltus fasciatus. None of these candidate genes exhibited conserved roles in these processes. Amidst this variation we wish to highlight the following themes for the evolution of sex determination: (1) Unique features within taxa influence network evolution. (2) Their position in the network influences a component's evolution. Our analyses also suggest an inverse association of protein sequence conservation with functional conservation.

Transcriptome changes of Apis mellifera female embryos with fem gene knockout by CRISPR/Cas9

The sex of honey bees is decided by a regulatory cascade comprising of csd, fem and Amdsx. In order to further identify other genes involved in sex determination and differentiation of honey bees in the early stages of embryo development, the CRISPR/Cas9 method was used to knock out fem gene in the embryonic stage of diploid western honey bees, and RNA-seq was used to analyze gene expression changes in the embryo after fem knockout. Finally, we found that the bees had undergone gender changes due to fem knockout. A total of 155 differentially expressed genes (DEGs) were obtained, with 48 up-regulated and 107 down-regulated DEGs in the mutant group compared to the control group. Of them, many genes are related to sex development or differentiation. In addition, 1502 differentially expressed alternative splicing events (DEASEs) related to 1011 genes, including the main honey bee sex-determining genes csd, tra2, fem, and Amdsx, were identified between the mutant group and control group, indicating that fem regulates alternative splicing of a large number of downstream genes. Our results provide valuable clues for further investigating the molecular mechanism of sex determination and differentiation in honey bees.

The function and evolution of a genetic switch controlling sexually dimorphic eye differentiation in honeybees

Animals develop sex-specific morphological structures that are diverse between organisms. However, understanding the developmental and evolutionary mechanisms governing these traits is still limited and largely restricted to DM domain genes, which are conserved, sex-specific developmental regulators identified in genetic models. Here, we report a sex-specific developmental regulator gene, glubschauge (glu) that selectively regulates sexually dimorphic eye differentiation in honeybees. We found that the sex determination gene feminizer (fem) controls sex-specific splicing of glu transcripts, establishing a genetic switch in which Glu proteins with a zinc finger (ZnF) domain are only expressed in females. We showed that female coding sequence was essential and sufficient for partial feminization. Comparative sequence and functional studies revealed that the evolutionary origination of the genetic switch was followed by the mutational origin of the essential ZnF domain. Our results demonstrate that glu is a newly evolved sex-specific genetic switch for region-specific regulation of a dimorphic character.

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.

Sex- and caste-specific transcriptomes of larval honey bee (Apis mellifera L.) gonads: DMRT A2 and Hsp83 are differentially expressed and regulated by juvenile hormone

The gonads of honey bee, Apis mellifera, queens and drones are each composed of hundreds of serial units, the ovarioles and testioles, while the ovaries of the adult subfertile workers consist of only few ovarioles. We performed a comparative RNA-seq analysis on early fifth-instar (L5F1) larval gonads, which is a critical stage in gonad development of honey bee larvae. A total of 1834 genes were identified as differentially expressed (Padj < 0.01) among the three sex and caste phenotypes. The Gene Ontology analysis showed significant enrichment for metabolism, protein or ion binding, and oxidoreductase activity, and a KEGG analysis revealed metabolic pathways as enriched. In a principal component analysis for the total transcriptomes and hierarchical clustering of the DEGs, we found higher similarity between the queen and worker ovary transcriptomes compared to the drone testis, despite the onset of programmed cell death in the worker ovaries. Four DEGs were selected for RT-qPCR analyses, including their response to juvenile hormone (JH), which is a critical factor in the caste-specific development of the ovaries. Among these, DMRT A2 and Hsp83 were found upregulated by JH and, thus, emerged as potential molecular markers for sex- and caste-specific gonad development in honey bees.

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.

Gene dynamics of haplodiploidy favor eusociality in the Hymenoptera

The problem of whether haplodiploidy is responsible for the frequent evolution of eusociality in the Hymenoptera remains unresolved. The little-known "protected invasion hypothesis" posits that because a male will transmit a new allele for alloparental care to all his daughters under haplodiploidy, such an allele has a higher probability of spreading to fixation under haplodiploidy than under diploidy. This mechanism is investigated using the mating system and lifecycles ancestral to eusocial lineages. It is shown that although haplodiploidy increases the probability of fixation of a new allele, the effect is cancelled by a higher probability of the allele arising in a diploid population. However, the same effect of male haploidy results in a 30% lower threshold amount of reproductive help by a worker necessary to favor eusociality if the sex ratio of dispersing first-brood offspring remains even. This occurs because when first-brood daughters become workers, the sex ratio of dispersing first-brood offspring becomes male-biased, selecting for an overall female-biased first-brood sex ratio. Through this mechanism, haplodiploidy may favor eusociality in the absence of a female-biased sex ratio in dispersing reproductive offspring. The gene-centric approach used here reveals the critical role of male haploidy in structuring the social group.

Absence of complementary sex determination in two Leptopilina species (Figitidae, Hymenoptera) and a reconsideration of its incompatibility with endosymbiont-induced thelytoky

Complementary sex determination (CSD) is a widespread sex determination mechanism in haplodiploid Hymenoptera. Under CSD, sex is determined by the allelic state of one or multiple CSD loci. Heterozygosity at one or more loci leads to female development, whereas hemizygosity of haploid eggs and homozygosity of diploid eggs results in male development. Sexual (arrhenotokous) reproduction normally yields haploid male and diploid female offspring. Under asexual reproduction (thelytoky), diploidized unfertilized eggs develop into females. Thelytoky is often induced by bacterial endosymbionts that achieve egg diploidization by gamete duplication. As gamete duplication leads to complete homozygosity, endosymbiont-induced thelytokous reproduction is presumed to be incompatible with CSD, which relies on heterozygosity for female development. Previously, we excluded CSD in four Asobara (Braconidae) species and proposed a two-step mechanism for Wolbachia-induced thelytoky in Asobara japonica. Here, we conclusively reject CSD in two cynipid wasp species, Leptopilina heterotoma and Leptopilina clavipes. We further show that thelytoky in L. clavipes depends on Wolbachia titer but that diploidization and feminization steps cannot be separated, unlike in A. japonica. We discuss what these results reveal about the sex determination mechanism of L. clavipes and the presumed incompatibility between CSD and endosymbiont-induced thelytoky in the Hymenoptera.

Analysis of Complementary Sex-Determiner (csd) Allele Diversity in Different Honeybee Subspecies from Italy Based on NGS Data

Sexual regulation in Apis mellifera is controlled by the complementary sex-determiner (csd) gene: females (queens and workers) are heterozygous at this locus and males (drones) are hemizygous. When homozygous diploid drones develop, they are eaten by worker bees. High csd allelic diversity in honeybee populations is a priority for colony survival. The focus of this study is to investigate csd variability in the genomic sequence of the hypervariable region (HVR) of the csd gene in honeybee subspecies sampled in Italy. During the summer of 2017 and 2018, worker bees belonging to 125 colonies were sampled. The honeybees belonged to seven different A. mellifera subspecies: A. m. ligustica, A. m. sicula, A. m cecropia, A. m. carnica, A. m. mellifera, Buckfast and hybrid Carnica. Illumina genomic resequencing of all samples was performed and used for the characterization of global variability among colonies. In this work, a pipeline using existing resequencing data to explore the csd gene allelic variants present in the subspecies collection, based on de novo assembly of sequences falling within the HVR region, is described. On the whole, 138 allelic sequences were successfully reconstructed. Among these, 88 different alleles were identified, 68 of which match with csd alleles present in the NCBI GenBank database.

Abundant small RNAs in the reproductive tissues and eggs of the honey bee, Apis mellifera

Background

Polyandrous social insects such as the honey bee are prime candidates for parental manipulation of gene expression in offspring. Although there is good evidence for parent-of-origin effects in honey bees the epigenetic mechanisms that underlie these effects remain a mystery. Small RNA molecules such as miRNAs, piRNAs and siRNAs play important roles in transgenerational epigenetic inheritance and in the regulation of gene expression during development.

Results

Here we present the first characterisation of small RNAs present in honey bee reproductive tissues: ovaries, spermatheca, semen, fertilised and unfertilised eggs, and testes. We show that semen contains fewer piRNAs relative to eggs and ovaries, and that piRNAs and miRNAs which map antisense to genes involved in DNA regulation and developmental processes are differentially expressed between tissues. tRNA fragments are highly abundant in semen and have a similar profile to those seen in the semen of other animals. Intriguingly we also find abundant piRNAs that target the sex determination locus, suggesting that piRNAs may play a role in honey bee sex determination.

Conclusions

We conclude that small RNAs may play a fundamental role in honey bee gametogenesis and reproduction and provide a plausible mechanism for parent-of-origin effects on gene expression and reproductive physiology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08478-9.

Sexual diversification of splicing regulation during embryonic development in honeybees (Apis mellifera), A haplodiploid system

The honeybee is a haplodiploid organism in which sexual development is determined by the complementary sex determiner (csd) gene and realized by sex-specific splicing processes involving the feminizer (fem) gene. We used high throughput transcriptome sequencing (RNA-Seq) to characterize the transcriptional differences between the sexes caused by the fertilization and sex determination processes in honeybee (Apis mellifera) embryos. We identified 758, 372 and 43 differentially expressed genes (DEGs) and 58, 176 and 233 differentially spliced genes (DSGs) in 10-15-h-old, 25-40-h-old and 55-70-h-old female and male embryos, respectively. The early difference in male and female embryos in response to the fertilization and non-fertilization processes resulted mainly in differential expression of genes (758 DEGs vs. 58 DSGs). In the latest sampled embryonic stage, the transcriptional differences between the sexes were dominated by alternative splicing of transcripts (43 DEGs vs. 233 DSGs). Interestingly, differentially spliced transcripts that encode RNA-binding properties were overrepresented in 55-70-h-old embryos, indicating a more diverse regulation via alternative splicing than previous work on the sex determination pathway suggested. These stage- and sex-specific transcriptome data from honeybee embryos provide a comprehensive resource for examining the roles of fertilization and sex determination in developmental programming in a haplodiploid system.

Manipulating Insect Sex Determination Pathways for Genetic Pest Management: Opportunities and Challenges

Sex determination pathways in insects are generally characterised by an upstream primary signal, which is highly variable across species, and that regulates the splicing of a suite of downstream but highly-conserved genes (transformer, doublesex and fruitless). In turn, these downstream genes then regulate the expression of sex-specific characteristics in males and females. Identification of sex determination pathways has and continues to be, a critical component of insect population suppression technologies. For example, "first-generation" transgenic technologies such as fsRIDL (Female-Specific Release of Insects carrying Dominant Lethals) enabled efficient selective removal of females from a target population as a significant improvement on the sterile insect technique (SIT). Second-generation technologies such as CRISPR/Cas9 homing gene drives and precision-guided SIT (pgSIT) have used gene editing technologies to manipulate sex determination genes in vivo. The development of future, third-generation control technologies, such as Y-linked drives, (female to male) sex-reversal, or X-shredding, will require additional knowledge of aspects of sexual development, including a deeper understanding of the nature of primary signals and dosage compensation. This review shows how knowledge of sex determination in target pest species is fundamental to all phases of the development of control technologies.

Alternative splicing patterns of doublesex reveal a missing link between Nix and doublesex in the sex determination cascade of Aedes albopictus

Sexual development in insects is regulated by a complicated hierarchical cascade of sex determination. The primary signals are diverse, whereas the central nexus doublesex (dsx) gene is relatively conserved within the pathway. Aedes (Stegomyia) albopictus is an important vector with an extensive worldwide distribution. We previously reported that Ae. albopictus dsx (Aalbdsx) yields one male- (Aalbdsx^M ) and three female-specific isoforms (Aalbdsx^F1-3 ); however, the spatiotemporal expression profiles and mechanisms regulating sex-specific alternative splicing require further investigation. In this study, we demonstrated that the Aalbdsx^M messenger RNA (mRNA) represents the default pattern when analyzed in human foreskin fibroblasts and HeLa cells. We combined reverse transcription polymerase chain reaction with RNA immunoprecipitation using specific antibodies against tagged Ae. albopictus male-determining factor AalNix and confirmed that AalNix indirectly regulates dsx pre-mRNA and regulates its alternative splicing. During the early embryo stage (0-2 and 4-8 h), maternal dsx^F and default splicing dsx^M were detected in both sexes; the expression of dsx^M then decreased until sufficient AalNix transcripts accumulated in male embryos at 20-24 h. These findings suggest that one or more potential dsx splicing enhancers can shift dsx^M to dsx^F in both sexes; however, the presence of Nix influences the function of this unknown splicing enhancer and ultimately leads to the formation of dsx^M in males. Finally, our results provide important insight into the regulatory mechanism of dsx alternative splicing in the mosquito.

A feminizing switch in a hemimetabolous insect

The mechanism of sex determination remains poorly understood in hemimetabolous insects. Here, in the brown planthopper (BPH), Nilaparvata lugens, a hemipteran rice pest, we identified a feminizing switch or a female determiner (Nlfmd) that encodes a serine/arginine-rich protein. Knockdown of Nlfmd in female nymphs resulted in masculinization of both the somatic morphology and doublesex splicing. The female-specific isoform of Nlfmd, Nlfmd-F, is maternally deposited and zygotically transcribed. Depletion of Nlfmd by maternal RNAi or CRISPR-Cas9 resulted in female-specific embryonic lethality. Knockdown of an hnRNP40 family gene named female determiner 2 (Nlfmd2) also conferred masculinization. In vitro experiments showed that an Nlfmd2 isoform, NlFMD2340, bound the RAAGAA repeat motif in the Nldsx pre-mRNA and formed a protein complex with NlFMD-F to modulate Nldsx splicing, suggesting that NlFMD2 may function as an RNA binding partner of the feminizing switch NlFMD. Our results provide novel insights into the diverse mechanisms of insect sex determination.

Effects of Functional Depletion of Doublesex on Male Development in the Sawfly, Athalia rosae

Simple Summary

The sawfly, Athalia rosae, exploits a haplodiploid mode of reproduction, in which fertilized eggs develop into diploid females, whereas unfertilized eggs parthenogenetically develop into haploid males. The doublesex (dsx) gene is a well-conserved transcription factor that regulates sexual differentiation in insects. In the present study, we knocked down the A. rosae ortholog of dsx (Ardsx) during several developmental stages with repeated double-stranded RNA (dsRNA) injections. As a result, knockdown of Ardsx in haploid males caused almost complete male-to-female sex reversal, but the resulting eggs were infertile. The same knockdown approach using diploid males caused complete male-to-female sex reversal; they were able to produce fertile eggs and exhibited female behaviors. The same RNAi treatment did not affect female differentiation. These results demonstrated that dsx in the sawfly is essential for male development and its depletion caused complete male-to-female sex reversal. This is the first demonstration of functional depletion of dsx not causing intersexuality but inducing total sex reversal in males instead.

Abstract

The doublesex (dsx) gene, which encodes a transcription factor, regulates sexual differentiation in insects. Sex-specific splicing of dsx occurs to yield male- and female-specific isoforms, which promote male and female development, respectively. Thus, functional disruption of dsx leads to an intersexual phenotype in both sexes. We previously identified a dsx ortholog in the sawfly, Athalia rosae. Similar to dsx in other insects, dsx in the sawfly yields different isoforms in males and females as a result of alternative splicing. The sawfly exploits a haplodiploid mode of reproduction, in which fertilized eggs develop into diploid females, whereas unfertilized eggs parthenogenetically develop into haploid males. In the present study, we knocked down the A. rosae ortholog of dsx (Ardsx) during several developmental stages with repeated double-stranded RNA (dsRNA) injections. Knockdown of Ardsx via parental RNA interference (RNAi), which enables knockdown of genes in offspring embryos, led to a lack of internal and external genitalia in haploid male progeny. Additional injection of dsRNA targeting Ardsx in these animals caused almost complete male-to-female sex reversal, but the resulting eggs were infertile. Notably, the same knockdown approach using diploid males obtained by sib-crossing caused complete male-to-female sex reversal; they were morphologically and behaviorally females. The same RNAi treatment did not affect female differentiation. These results indicate that dsx in the sawfly is essential for male development and its depletion caused complete male-to-female sex reversal. This is the first demonstration of functional depletion of dsx not causing intersexuality but inducing total sex reversal in males instead.

Sex-Specific Development in Haplodiploid Honeybee Is Controlled by the Female-Embryo-Specific Activation of Thousands of Intronic LncRNAs

Embryonic development depends on a highly coordinated shift in transcription programs known as the maternal-to-zygotic transition (MZT). It remains unclear how haploid and diploid embryo coordinate their genomic activation and embryonic development during MZT in haplodiploid animals. Here, we applied a single-embryo RNA-seq approach to characterize the embryonic transcriptome dynamics in haploid males vs. diploid females of the haplodiploid insect honeybee (Apis mellifera). We observed typical zygotic genome activation (ZGA) occurred in three major waves specifically in female honeybee embryos; haploid genome activation was much weaker and occurred later. Strikingly, we also observed three waves of transcriptional activation for thousands of long non-coding transcripts (lncRNA), 73% of which are transcribed from intronic regions and 65% were specific to female honeybee embryos. These findings support a model in which introns encode thousands of lncRNAs that are expressed in a diploid-embryo-specific and ZGA-triggered manner that may have potential functions to regulate gene expression during early embryonic development in the haplodiploid insect honeybee.

Identification and Characterization of the Masculinizing Function of the Helicoverpa armigera Masc Gene

The Masculinizer (Masc) gene has been known to control sex development and dosage compensation in lepidopterans. However, it remains unclear whether its ortholog exists and plays the same roles in distantly related lepidopterans such as Helicoverpa armigera. To address this question, we cloned Masc from H. armigera (HaMasc), which contains all essential functional domains of BmMasc, albeit with less than 30% amino acid sequence identity with BmMasc. Genomic PCR and qPCR analyses showed that HaMasc is a Z chromosome-linked gene since its genomic content in males (ZZ) was two times greater than that in females (ZW). RT-PCR and RT-qPCR analyses revealed that HaMasc expression was sex- and stage-biased, with significantly more transcripts in males and eggs than in females and other stages. Transfection of a mixture of three siRNAs of HaMasc into a male embryonic cell line of H. armigera led to the appearance of female-specific mRNA splicing isoforms of H. armigeradoublesex (Hadsx), a downstream target gene of HaMasc in the H. armigera sex determination pathway. The knockdown of HaMasc, starting from the third instar larvae resulted in a shift of Hadsx splicing from male to female isoforms, smaller male pupa and testes, fewer but larger/longer spermatocytes and sperm bundles, delayed pupation and internal fusion of the testes and follicles. These data demonstrate that HaMasc functions as a masculinizing gene in the H. armigera sex-determination cascade.

Transcriptional Profiles of Diploid Mutant Apis mellifera Embryos after Knockout of csd by CRISPR/Cas9

Simple Summary

In honey bees, males are haploid while females are diploid, leading to a fundamental difference in genetic materials between the sexes. In order to better control the comparison of gene expression between males and females, diploid mutant males were generated by knocking out the sex-determining gene, complementary sex determiner (csd), in fertilized embryos. The diploid mutant drones had male external morphological features, as well as male gonads. RNA sequencing was performed on the diploid mutant embryos and one-day-old larvae. The transcriptome analysis showed that several female-biased genes, such as worker-enriched antennal (Wat), vitellogenin (Vg), and some venom-related genes, were down-regulated in the diploid mutant males. In contrast, some male-biased genes, like takeout and apolipophorin-III-like protein (A4), were up-regulated. Moreover, the co-expression gene networks suggested that csd might interact very closely with fruitless (fru), feminizer (fem) might have connections with hexamerin 70c (hex70c), and transformer-2 (tra2) might play roles with troponin T (TpnT). Foundational information about the differences in the gene expression caused by sex differentiation was provided in this study. It is believed that this study will pave the ground for further research on the different mechanisms between males and females in honey bees.

Abstract

In honey bees, complementary sex determiner (csd) is the primary signal of sex determination. Its allelic composition is heterozygous in females, and hemizygous or homozygous in males. To explore the transcriptome differences after sex differentiation between males and females, with genetic differences excluded, csd in fertilized embryos was knocked out by CRISPR/Cas9. The diploid mutant males at 24 h, 48 h, 72 h, and 96 h after egg laying (AEL) and the mock-treated females derived from the same fertilized queen were investigated through RNA-seq. Mutations were detected in the target sequence in diploid mutants. The diploid mutant drones had typical male morphological characteristics and gonads. Transcriptome analysis showed that several female-biased genes, such as worker-enriched antennal (Wat), vitellogenin (Vg), and some venom-related genes, were down-regulated in the diploid mutant males. In contrast, some male-biased genes, such as takeout and apolipophorin-III-like protein (A4), had higher expressions in the diploid mutant males. Weighted gene co-expression network analysis (WGCNA) indicated that there might be interactions between csd and fruitless (fru), feminizer (fem) and hexamerin 70c (hex70c), transformer-2 (tra2) and troponin T (TpnT). The information provided by this study will benefit further research on the sex dimorphism and development of honey bees and other insects in Hymenoptera.

A conserved role of the duplicated Masculinizer gene in sex determination of the Mediterranean flour moth, Ephestia kuehniella

Sex determination in the silkworm, Bombyx mori, is based on Feminizer (Fem), a W-linked Fem piRNA that triggers female development in WZ individuals, and the Z-linked Masculinizer (Masc), which initiates male development and dosage compensation in ZZ individuals. While Fem piRNA is missing in a close relative of B. mori, Masc determines sex in several representatives of distant lepidopteran lineages. We studied the molecular mechanisms of sex determination in the Mediterranean flour moth, Ephestia kuehniella (Pyralidae). We identified an E. kuehniella Masc ortholog, EkMasc, and its paralog resulting from a recent duplication, EkMascB. Both genes are located on the Z chromosome and encode a similar Masc protein that contains two conserved domains but has lost the conserved double zinc finger domain. We developed PCR-based genetic sexing and demonstrated a peak in the expression of EkMasc and EkMascB genes only in early male embryos. Simultaneous knock-down experiments of both EkMasc and EkMascB using RNAi during early embryogenesis led to a shift from male- to female-specific splicing of the E. kuehniella doublesex gene (Ekdsx), their downstream effector, in ZZ embryos and resulted in a strong female-biased sex-ratio. Our results thus confirmed the conserved role of EkMasc and/or EkMascB in masculinization. We suggest that the C-terminal proline-rich domain, we have identified in all functionally confirmed Masc proteins, in conjunction with the masculinizing domain, is important for transcriptional regulation of sex determination in Lepidoptera. The function of the Masc double zinc finger domain is still unknown, but appears to have been lost in E. kuehniella.

The sex-determining cascade in the silkworm, Bombyx mori, differs greatly from those of other insects. In B. mori, female development is initiated by Fem piRNA expressed from the W chromosome during early embryogenesis. Fem piRNA silences Masculinizer (Masc) thereby blocking the male pathway resulting in female development. It is currently unknown whether this cascade is conserved across Lepidoptera. In the Mediterranean flour moth, Ephestia kuehniella, we identified an ortholog of Masc and discovered its functional duplication on the Z chromosome, which has not yet been found in any other lepidopteran species. We provide two lines of evidence that the EkMasc and/or EkMascB genes play an essential role in masculinization: (i) they show a peak of expression during early embryogenesis in ZZ but not in WZ embryos and (ii) their simultaneous silencing by RNAi results in female-specific splicing of the E. kuehniella doublesex gene (Ekdsx) in ZZ embryos and in a female-biased sex ratio. Our results suggest a conserved role of the duplicated Masc gene in sex determination of E. kuehniella.

Modulation of fatty acid elongation in cockroaches sustains sexually dimorphic hydrocarbons and female attractiveness

Insect cuticular hydrocarbons (CHCs) serve as important intersexual signaling chemicals and generally show variation between the sexes, but little is known about the generation of sexually dimorphic hydrocarbons (SDHCs) in insects. In this study, we report the molecular mechanism and biological significance that underlie the generation of SDHC in the German cockroach Blattella germanica. Sexually mature females possess more C29 CHCs, especially the contact sex pheromone precursor 3,11-DimeC29. RNA interference (RNAi) screen against the fatty acid elongase family members combined with heterologous expression of the genes in yeast revealed that both BgElo12 and BgElo24 were involved in hydrocarbon (HC) production, but BgElo24 is of wide catalytic activities and is able to provide substrates for BgElo12, and only the female-enriched BgElo12 is responsible for sustaining female-specific HC profile. Repressing BgElo12 masculinized the female CHC profile, decreased contact sex pheromone level, and consequently reduced the sexual attractiveness of female cockroaches. Moreover, the asymmetric expression of BgElo12 between the sexes is modulated by sex differentiation cascade. Specifically, male-specific BgDsx represses the transcription of BgElo12 in males, while BgTra is able to remove this effect in females. Our study reveals a novel molecular mechanism responsible for the formation of SDHCs and also provide evidences on shaping of the SDHCs by sexual selection, as females use them to generate high levels of contact sex pheromone.

Sexual dimorphism of body waxes is prevalent in insects; this study reveals that the sex-differentiation pathway regulates fatty acid elongation, ensuring production of the sexually dimorphic cuticular hydrocarbons needed for high levels of sex pheromone and sexual attractiveness in female cockroaches.

Diversity of Modes of Reproduction and Sex Determination Systems in Invertebrates, and the Putative Contribution of Genetic Conflict

About eight million animal species are estimated to live on Earth, and all except those belonging to one subphylum are invertebrates. Invertebrates are incredibly diverse in their morphologies, life histories, and in the range of the ecological niches that they occupy. A great variety of modes of reproduction and sex determination systems is also observed among them, and their mosaic-distribution across the phylogeny shows that transitions between them occur frequently and rapidly. Genetic conflict in its various forms is a long-standing theory to explain what drives those evolutionary transitions. Here, we review (1) the different modes of reproduction among invertebrate species, highlighting sexual reproduction as the probable ancestral state; (2) the paradoxical diversity of sex determination systems; (3) the different types of genetic conflicts that could drive the evolution of such different systems.

Genus-Wide Characterization of Bumblebee Genomes Provides Insights into Their Evolution and Variation in Ecological and Behavioral Traits

Bumblebees are a diverse group of globally important pollinators in natural ecosystems and for agricultural food production. With both eusocial and solitary life-cycle phases, and some social parasite species, they are especially interesting models to understand social evolution, behavior, and ecology. Reports of many species in decline point to pathogen transmission, habitat loss, pesticide usage, and global climate change, as interconnected causes. These threats to bumblebee diversity make our reliance on a handful of well-studied species for agricultural pollination particularly precarious. To broadly sample bumblebee genomic and phenotypic diversity, we de novo sequenced and assembled the genomes of 17 species, representing all 15 subgenera, producing the first genus-wide quantification of genetic and genomic variation potentially underlying key ecological and behavioral traits. The species phylogeny resolves subgenera relationships, whereas incomplete lineage sorting likely drives high levels of gene tree discordance. Five chromosome-level assemblies show a stable 18-chromosome karyotype, with major rearrangements creating 25 chromosomes in social parasites. Differential transposable element activity drives changes in genome sizes, with putative domestications of repetitive sequences influencing gene coding and regulatory potential. Dynamically evolving gene families and signatures of positive selection point to genus-wide variation in processes linked to foraging, diet and metabolism, immunity and detoxification, as well as adaptations for life at high altitudes. Our study reveals how bumblebee genes and genomes have evolved across the Bombus phylogeny and identifies variations potentially linked to key ecological and behavioral traits of these important pollinators.

Comparative Transcriptome Analysis Reveals Sex-Based Differences during the Development of the Adult Parasitic Wasp Cotesia vestalis (Hymenoptera: Braconidae)

The endoparasitic wasp Cotesia vestalis is an important biological agent for controlling the population of Plutella xylostella, a major pest of cruciferous crops worldwide. Though the genome of C. vestalis has recently been reported, molecular mechanisms associated with sexual development have not been comprehensively studied. Here, we combined PacBio Iso-Seq and Illumina RNA-Seq to perform genome-wide profiling of pharate adult and adult development of male and female C. vestalis. Taking advantage of Iso-Seq full-length reads, we identified 14,466 novel transcripts as well as 8770 lncRNAs, with many lncRNAs showing a sex- and stage-specific expression pattern. The differentially expressed gene (DEG) analyses showed 2125 stage-specific and 326 sex-specific expressed genes. We also found that 4819 genes showed 11,856 alternative splicing events through combining the Iso-Seq and RNA-Seq data. The results of comparative analyses showed that most genes were alternatively spliced across developmental stages, and alternative splicing (AS) events were more prevalent in females than in males. Furthermore, we identified six sex-determining genes in this parasitic wasp and verified their sex-specific alternative splicing profiles. Specifically, the characterization of feminizer and doublesex splicing between male and female implies a conserved regulation mechanism of sexual development in parasitic wasps.

Conservation insights from wild bee genetic studies: Geographic differences, susceptibility to inbreeding, and signs of local adaptation

Conserving bees are critical both ecologically and economically. Genetic tools are valuable for monitoring these vital pollinators since tracking these small, fast-flying insects by traditional means is difficult. By surveying the current state of the literature, this review discusses how recent advances in landscape genetic and genomic research are elucidating how wild bees respond to anthropogenic threats. Current literature suggests that there may be geographic differences in the vulnerability of bee species to landscape changes. Populations of temperate bee species are becoming more isolated and more genetically depauperate as their landscape becomes more fragmented, but tropical bee species appear unaffected. These differences may be an artifact of historical differences in land-use, or it suggests that different management plans are needed for temperate and tropical bee species. Encouragingly, genetic studies on invasive bee species indicate that low levels of genetic diversity may not lead to rapid extinction in bees as once predicted. Additionally, next-generation sequencing has given researchers the power to identify potential genes under selection, which are likely critical to species' survival in their rapidly changing environment. While genetic studies provide insights into wild bee biology, more studies focusing on a greater phylogenetic and life-history breadth of species are needed. Therefore, caution should be taken when making broad conservation decisions based on the currently few species examined.

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.

A chimeric gene paternally instructs female sex determination in the haplodiploid wasp Nasonia

Various primary signals direct insect sex determination. In hymenopteran insects, the presence of a paternal genome is needed to initiate female development. When absent, uniparental haploid males develop. We molecularly and functionally identified the instructor sex-determination gene, wasp overruler of masculinization (wom), of the haplodiploid wasp Nasonia vitripennis This gene contains a P53-like domain coding region and arose by gene duplication and genomic rearrangements. Maternal silencing of wom results in male development of haploid embryos. Upon fertilization, early zygotic transcription from the paternal wom allele is initiated, followed by a timely zygotic expression of transformer (tra), leading to female development. Wom is an instructor gene with a parent-of-origin effect in sex determination.

Bracon brevicornis Genome Showcases the Potential of Linked-Read Sequencing in Identifying a Putative Complementary Sex Determiner Gene

Bracon brevicornis is an ectoparasitoid of a wide range of larval-stage Lepidopterans, including several pests of important crops, such as the corn borer, Ostrinia nubilalis. It is also one of the earliest documented cases of complementary sex determination in Hymenoptera. Here, we present the linked-read-based genome of B. brevicornis, complete with an ab initio-derived annotation and protein comparisons with fellow braconids, Fopius arisanus and Diachasma alloeum. We demonstrate the potential of linked-read assemblies in exploring regions of heterozygosity and search for structural and homology-derived evidence of the complementary sex determiner gene (csd).

Nix alone is sufficient to convert female Aedes aegypti into fertile males and myo-sex is needed for male flight

The presence of a dominant male-determining locus (M-locus) in one of a pair of autosomes establishes the male sex in the dengue fever mosquito Aedes aegypti. The Ae. aegypti M-locus contains 30 genes, including Nix, a previously reported male-determining factor. Here we show that the Nix transgene alone was sufficient to convert females into fertile males, which continued to produce sex-converted progeny. We also show that a second M-locus gene named myo-sex was needed for male flight. Nix-mediated sex conversion was 100% penetrant, heritable, and stable, indicating great potential for developing mosquito-control strategies to reduce vector populations by female-to-male conversion. This work also sheds lights into the molecular basis of the function of the M-locus.

A dominant male-determining locus (M-locus) establishes the male sex (M/m) in the yellow fever mosquito, Aedes aegypti. Nix, a gene in the M-locus, was shown to be a male-determining factor (M factor) as somatic knockout of Nix led to feminized males (M/m) while transient expression of Nix resulted in partially masculinized females (m/m), with male reproductive organs but retained female antennae. It was not clear whether any of the other 29 genes in the 1.3-Mb M-locus are also needed for complete sex-conversion. Here, we report the generation of multiple transgenic lines that express Nix under the control of its own promoter. Genetic and molecular analyses of these lines provided insights unattainable from previous transient experiments. We show that the Nix transgene alone, in the absence of the M-locus, was sufficient to convert females into males with all male-specific sexually dimorphic features and male-like gene expression. The converted m/m males are flightless, unable to perform the nuptial flight required for mating. However, they were able to father sex-converted progeny when presented with cold-anesthetized wild-type females. We show that myo-sex, a myosin heavy-chain gene also in the M-locus, was required for male flight as knockout of myo-sex rendered wild-type males flightless. We also show that Nix-mediated female-to-male conversion was 100% penetrant and stable over many generations. Therefore, Nix has great potential for developing mosquito control strategies to reduce vector populations by female-to-male sex conversion, or to aid in a sterile insect technique that requires releasing only non-biting males.

Genome of the Parasitoid Wasp Diachasma alloeum, an Emerging Model for Ecological Speciation and Transitions to Asexual Reproduction

Parasitoid wasps are among the most speciose animals, yet have relatively few available genomic resources. We report a draft genome assembly of the wasp Diachasma alloeum (Hymenoptera: Braconidae), a host-specific parasitoid of the apple maggot fly Rhagoletis pomonella (Diptera: Tephritidae), and a developing model for understanding how ecological speciation can “cascade” across trophic levels. Identification of gene content confirmed the overall quality of the draft genome, and we manually annotated ∼400 genes as part of this study, including those involved in oxidative phosphorylation, chemosensation, and reproduction. Through comparisons to model hymenopterans such as the European honeybee Apis mellifera and parasitoid wasp Nasonia vitripennis, as well as a more closely related braconid parasitoid Microplitis demolitor, we identified a proliferation of transposable elements in the genome, an expansion of chemosensory genes in parasitoid wasps, and the maintenance of several key genes with known roles in sexual reproduction and sex determination. The D. alloeum genome will provide a valuable resource for comparative genomics studies in Hymenoptera as well as specific investigations into the genomic changes associated with ecological speciation and transitions to asexuality.

Mapping of Multiple Complementary Sex Determination Loci in a Parasitoid Wasp

Sex determination has evolved in a variety of ways and can depend on environmental and genetic signals. A widespread form of genetic sex determination is haplodiploidy, where unfertilized, haploid eggs develop into males and fertilized diploid eggs into females. One of the molecular mechanisms underlying haplodiploidy in Hymenoptera, the large insect order comprising ants, bees, and wasps, is complementary sex determination (CSD). In species with CSD, heterozygosity at one or several loci induces female development. Here, we identify the genomic regions putatively underlying multilocus CSD in the parasitoid wasp Lysiphlebus fabarum using restriction-site associated DNA sequencing. By analyzing segregation patterns at polymorphic sites among 331 diploid males and females, we identify up to four CSD candidate regions, all on different chromosomes. None of the candidate regions feature evidence for homology with the csd gene from the honey bee, the only species in which CSD has been characterized, suggesting that CSD in L. fabarum is regulated via a novel molecular mechanism. Moreover, no homology is shared between the candidate loci, in contrast to the idea that multilocus CSD should emerge from duplications of an ancestral single-locus system. Taken together, our results suggest that the molecular mechanisms underlying CSD in Hymenoptera are not conserved between species, raising the question as to whether CSD may have evolved multiple times independently in the group.

Nix is a male-determining factor in the Asian tiger mosquito Aedes albopictus

The initial signal that governs sex determination is highly variable among insects. A homolog of Nix, the male-determining factor in Aedes aegypti, was previously found in the Asian tiger mosquito Ae. albopictus. Here we show that the Ae. albopictus Nix (AalNix) is more complex in gene structure and splice isoforms than its Ae. aegypti homolog (AaeNix). AalNix shows a similar transcription profile compared to AaeNix. CRISPR/Cas9-mediated knockouts of AalNix in vivo and in the Ae. albopictus C6/36 cells lead to a shift of dsx and fru splicing towards the female isoforms. G₀ knockout males showed feminization and deformities including feminized antennae, absence or partial absence of gonocoxites, gonostyli, testes and accessory glands, and the formation of ovaries. Despite ~70 MY of divergence, Nix functions as a conserved male-determining factor in the two most important arboviral vectors, namely Ae. aegypti and Ae. albopictus.

Identification of potential candidate genes involved in the sex determination cascade in an aquatic firefly, Sclerotia aquatilis (Coleoptera, Lampyridae)

Sexual differentiation, dimorphism, and courtship behavior are the downstream developmental programs of the sex determination cascade. The sex determination cascade in arthropods often involves key genes, transformer (tra), doublesex (dsx), transformer-2 (tra2), and fruitless (fru). These genes are conserved among insect taxa; however, they have never been reported in fireflies. In this study, the candidate genes for these key genes were identified for the first time in an aquatic firefly, Sclerotia aquatilis using transcriptome analysis. A comparative protein-protein interaction (PPI) network of sex determination cascade was reconstructed for S. aquatilis based on a network of a model insect, Drosophila melanogaster. Subsequently, a sex determination cascade in S. aquatilis was proposed based on the amino acid sequence structures and expression profiles of these candidates. This study describes the first efforts toward understanding the molecular control of sex determination cascade in fireflies.

Hemimetabolous insects elucidate the origin of sexual development via alternative splicing

Insects are the only known animals in which sexual differentiation is controlled by sex-specific splicing. The doublesex transcription factor produces distinct male and female isoforms, which are both essential for sex-specific development. dsx splicing depends on transformer, which is also alternatively spliced such that functional Tra is only present in females. This pathway has evolved from an ancestral mechanism where dsx was independent of tra and expressed and required only in males. To reconstruct this transition, we examined three basal, hemimetabolous insect orders: Hemiptera, Phthiraptera, and Blattodea. We show that tra and dsx have distinct functions in these insects, reflecting different stages in the changeover from a transcription-based to a splicing-based mode of sexual differentiation. We propose that the canonical insect tra-dsx pathway evolved via merger between expanding dsx function (from males to both sexes) and narrowing tra function (from a general splicing factor to dedicated regulator of dsx).

Substantial Heritable Variation in Recombination Rate on Multiple Scales in Honeybees and Bumblebees

Meiotic recombination shuffles genetic variation and promotes correct segregation of chromosomes. Rates of recombination vary on several scales, both within genomes and between individuals, and this variation is affected by both genetic and environmental factors. Social insects have extremely high rates of recombination, although the evolutionary causes of this are not known. Here, we estimate rates of crossovers and gene conversions in 22 colonies of the honeybee, Apis mellifera, and 9 colonies of the bumblebee, Bombus terrestris, using direct sequencing of 299 haploid drone offspring. We confirm that both species have extremely elevated crossover rates, with higher rates measured in the highly eusocial honeybee than the primitively social bumblebee. There are also significant differences in recombination rate between subspecies of honeybee. There is substantial variation in genome-wide recombination rate between individuals of both A. mellifera and B. terrestris and the distribution of these rates overlap between species. A large proportion of interindividual variation in recombination rate is heritable, which indicates the presence of variation in trans-acting factors that influence recombination genome-wide. We infer that levels of crossover interference are significantly lower in honeybees compared to bumblebees, which may be one mechanism that contributes to higher recombination rates in honeybees. We also find a significant increase in recombination rate with distance from the centromere, mirrored by methylation differences. We detect a strong transmission bias due to GC-biased gene conversion associated with noncrossover gene conversions. Our results shed light on the mechanistic causes of extreme rates of recombination in social insects and the genetic architecture of recombination rate variation.