Duplication and concerted evolution in a master sex determiner under balancing selection

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.

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

Genomic insights into mite phylogeny, fitness, development, and reproduction

BACKGROUND

Predatory mites (Acari: Phytoseiidae) are the most important beneficial arthropods used in augmentative biological pest control of protected crops around the world. However, the genomes of mites are far less well understood than those of insects and the evolutionary relationships among mite and other chelicerate orders are contested, with the enigmatic origin of mites at one of the centres in discussion of the evolution of Arachnida.

RESULTS

We here report the 173 Mb nuclear genome (from 51.75 Gb pairs of Illumina reads) of the predatory mite, Neoseiulus cucumeris, a biocontrol agent against pests such as mites and thrips worldwide. We identified nearly 20.6 Mb (~ 11.93% of this genome) of repetitive sequences and annotated 18,735 protein-coding genes (a typical gene 2888 bp in size); the total length of protein-coding genes was about 50.55 Mb (29.2% of this assembly). About 37% (6981) of the genes are unique to N. cucumeris based on comparison with other arachnid genomes. Our phylogenomic analysis supported the monophyly of Acari, therefore rejecting the biphyletic origin of mites advocated by other studies based on limited gene fragments or few taxa in recent years. Our transcriptomic analyses of different life stages of N. cucumeris provide new insights into genes involved in its development. Putative genes involved in vitellogenesis, regulation of oviposition, sex determination, development of legs, signal perception, detoxification and stress-resistance, and innate immune systems are identified.

CONCLUSIONS

Our genomics and developmental transcriptomics analyses of N. cucumeris provide invaluable resources for further research on the development, reproduction, and fitness of this economically important mite in particular and Arachnida in general.

Evidence for involvement of a transformer paralogue in sex determination of the wasp Leptopilina clavipes

Transformer (tra) is the central gear in many insect sex determination pathways and transduces a wide range of primary signals. Mediated by transformer-2 (tra2) it directs sexual development into the female or male mode. Duplications of tra have been detected in numerous Hymenoptera, but a function in sex determination has been confirmed only in Apis mellifera. We identified a tra2 orthologue (Lc-tra2), a tra orthologue (Lc-tra) and a tra paralogue (Lc-traB) in the genome of Leptopilina clavipes (Hymenoptera: Cynipidae). We compared the sequence and structural conservation of these genes between sexual (arrhenotokous) and asexual all-female producing (thelytokous) individuals. Lc-tra is sex-specifically spliced in adults consistent with its orthologous function. The male-specific regions of Lc-tra are conserved in both reproductive modes. The paralogue Lc-traB lacks the genomic region coding for male-specific exons and can only be translated into a full-length TRA-like peptide sequence. Furthermore, unlike LC-TRA, the LC-TRAB interstrain sequence variation is not differentiated into a sexual and an asexual haplotype. The LC-TRAB protein interacts with LC-TRA as well as LC-TRA2. This suggests that Lc-traB functions as a conserved element in sex determination of sexual and asexual individuals.

The doublesex gene integrates multi-locus complementary sex determination signals in the Japanese ant, Vollenhovia emeryi

A female diploid, male haploid sex determination system (haplodiploidy) is found in hymenopteran taxa, such as ants, wasps, bees and sawflies. In this system, a single, complementary sex-determination (sl-CSD) locus functions as the primary sex-determination signal. In the taxa that has evolved this system, females and males are heterozygous and hemi/homozygous at the CSD locus, respectively. While the sl-CSD system enables females to alter sex ratios in the nest, it carries a high cost in terms of inbreeding, as individuals that are homozygous at the CSD locus become sterile diploid males. To counter this risk, some of hymenopteran species have evolved a multi-locus CSD (ml-CSD) system, which effectively reduces the proportion of sterile males. However, the mechanism by which these multiple primary signals are integrated and how they affect the terminal sex-differentiation signal of the molecular cascade have not yet been clarified. To resolve these questions, we examined the molecular cascade in the Japanese ant Vollenhovia emeryi, which we previously confirmed has two CSD loci. Here, we showed that the sex-determination gene, doublesex (dsx), which is highly conserved among phylogenetically distant taxa, is responsible for integrating two CSD signals in V. emeryi. After identifying and characterizing dsx, genotypes containing two CSD loci and splicing patterns of dsx were found to correspond to the sexual phenotype, suggesting that two primary signals are integrated into dsx. These findings will facilitate future molecular and functional studies of the sex determination cascade in V. emeryi, and shed light on the evolution and diversification of sex determination systems in insects.

Fire ant social chromosomes: Differences in number, sequence and expression of odorant binding proteins

Variation in social behavior is common yet our knowledge of the mechanisms underpinning its evolution is limited. The fire ant Solenopsis invicta provides a textbook example of a Mendelian element controlling social organization: alternate alleles of a genetic element first identified as encoding an odorant binding protein (OBP) named Gp-9 determine whether a colony accepts one or multiple queens. The potential roles of such a protein in perceiving olfactory cues and evidence of positive selection on its amino acid sequence made it an appealing candidate gene. However, we recently showed that recombination is suppressed between Gp-9 and hundreds of other genes as part of a >19 Mb supergene-like region carried by a pair of social chromosomes. This finding raises the need to reassess the potential role of Gp-9. We identify 23 OBPs in the fire ant genome assembly, including nine located in the region of suppressed recombination with Gp-9. For six of these, the alleles carried by the two variants of the supergene-like region differ in protein-coding sequence and thus likely in function, with Gp-9 showing the strongest evidence of positive selection. We identify an additional OBP specific to the Sb variant of the region. Finally, we find that 14 OBPs are differentially expressed between single- and multiple-queen colonies. These results are consistent with multiple OBPs playing a role in determining social structure.

Identification and functional characterization of the sex-determining gene doublesex in the sawfly, Athalia rosae (Hymenoptera: Tenthredinidae)

Sexual fate of the sawfly, Athalia rosae (Hymenoptera: Tenthredinidae) is determined by the complementary sex determination (CSD) mechanism as is the case in honeybees. However, to date, genes involved in sex determination have not been identified in this species. In this study, we attempted to identify orthologs of complementary sex-determiner (csd), feminizer (fem), and doublesex (dsx) from the A. rosae genome, all of which are crucial components of the sex determination cascade in the honeybee. As a result, we identified a sawfly ortholog of dsx (designated as Ardsx). Rapid amplification of cDNA ends (RACE) using total RNA extracted from male and female larvae identified three male-specific variants and three female-specific variants. Comparison between the full-length Ardsx cDNAs and the genomic sequence revealed that exon 5 was differentially spliced between the male- and female-specific variants. RT-PCR analysis demonstrated that Ardsx pre-mRNA was spliced alternatively in a sex-dependent manner at almost all the developmental stages. RNAi-mediated knockdown of Ardsx in males caused severe defects in the reproductive organs and, notably, induced development of the ovipository apparatus containing the dorsal pair of blades and the sheath. These males also showed abnormalities in testes and seminal vesicles and lacked mature sperm. The present study provides the first direct evidence that dsx is essential for sexual development in hymenopteran species.

The transformer genes in the fig wasp Ceratosolen solmsi provide new evidence for duplications independent of complementary sex determination

Transformer (tra) is the key gene that turns on the sex-determination cascade in Drosophila melanogaster and in some other insects. The honeybee Apis mellifera has two duplicates of tra, one of which (complementary sex determiner, csd) is the primary signal for complementary sex-determination (CSD), regulating the other duplicate (feminizer). Two tra duplicates have been found in some other hymenopteran species, resulting in the assumption that a single ancestral duplication of tra took place in the Hymenoptera. Here, we searched for tra homologues and pseudogenes in the Hymenoptera, focusing on five newly published hymenopteran genomes. We found three tra copies in the fig wasp Ceratosolen solmsi. Further evolutionary and expression analyses also showed that the two duplicates (Csoltra-B and Csoltra-C) are under positive selection, and have female-specific expression, suggesting possible sex-related functions. Moreover, Aculeata species exhibit many pseudogenes generated by lineage-specific duplications. We conclude that phylogenetic reconstruction and pseudogene screening provide novel evidence supporting the hypothesis of independent duplications rather an ancestral origin of multiple tra paralogues in the Hymenoptera. The case of C. solmsi is the first example of a non-CSD species with duplicated tra, contrary to the previous assumption that derived tra paralogues function as the CSD locus.

Evolution of Social Insect Polyphenism Facilitated by the Sex Differentiation Cascade

The major transition to eusociality required the evolution of a switch to canalize development into either a reproductive or a helper, the nature of which is currently unknown. Following predictions from the ‘theory of facilitated variation’, we identify sex differentiation pathways as promising candidates because of their pre-adaptation to regulating development of complex phenotypes. We show that conserved core genes, including the juvenile hormone-sensitive master sex differentiation gene doublesex (dsx) and a krüppel homolog 2 (kr-h2) with putative regulatory function, exhibit both sex and morph-specific expression across life stages in the ant Cardiocondyla obscurior. We hypothesize that genes in the sex differentiation cascade evolved perception of alternative input signals for caste differentiation (i.e. environmental or genetic cues), and that their inherent switch-like and epistatic behavior facilitated signal transfer to downstream targets, thus allowing them to control differential development into morphological castes.

Division of labor into reproductive queens and helper workers in the societies of ants, bees and wasps is achieved by phenotypic plasticity, which allows individuals to embark on discrete developmental trajectories in response to variable signals. These signals can be genetic, epigenetic or environmental, thereby resembling the extreme variation in signals for sex determination across multicellular animals. We show that common developmental pathways downstream of these input signals, including the conserved sex differentiation gene doublesex, regulate sex and caste-specific phenotypic differentiation in the ant species Cardiocondyla obscurior. Many different mechanisms of gene regulation have been implicated in controlling caste-specific development in social insects but these all depend on a higher-level genetic switch. We propose that highly conserved hub genes such as dsx, which can translate variable input signals into large transcription differences using intermediate-level regulators, are tightly linked with the repeated evolutionary transition to eusociality and caste polyphenism.

QTL Mapping of Sex Determination Loci Supports an Ancient Pathway in Ants and Honey Bees

Sex determination mechanisms play a central role in life-history characteristics, affecting mating systems, sex ratios, inbreeding tolerance, etc. Downstream components of sex determination pathways are highly conserved, but upstream components evolve rapidly. Evolutionary dynamics of sex determination remain poorly understood, particularly because mechanisms appear so diverse. Here we investigate the origins and evolution of complementary sex determination (CSD) in ants and bees. The honey bee has a well-characterized CSD locus, containing tandemly arranged homologs of the transformer gene [complementary sex determiner (csd) and feminizer (fem)]. Such tandem paralogs appear frequently in aculeate hymenopteran genomes. However, only comparative genomic, but not functional, data support a broader role for csd/fem in sex determination, and whether species other than the honey bee use this pathway remains controversial. Here we used a backcross to test whether csd/fem acts as a CSD locus in an ant (Vollenhovia emeryi). After sequencing and assembling the genome, we computed a linkage map, and conducted a quantitative trait locus (QTL) analysis of diploid male production using 68 diploid males and 171 workers. We found two QTLs on separate linkage groups (CsdQTL1 and CsdQTL2) that jointly explained 98.0% of the phenotypic variance. CsdQTL1 included two tandem transformer homologs. These data support the prediction that the same CSD mechanism has indeed been conserved for over 100 million years. CsdQTL2 had no similarity to CsdQTL1 and included a 236-kb region with no obvious CSD gene candidates, making it impossible to conclusively characterize it using our data. The sequence of this locus was conserved in at least one other ant genome that diverged >75 million years ago. By applying QTL analysis to ants for the first time, we support the hypothesis that elements of hymenopteran CSD are ancient, but also show that more remains to be learned about the diversity of CSD mechanisms.

Whether a developing embryo becomes male or female has significant downstream consequences. Depending on the species, sex can be determined by a wide variety of mechanisms. Sex determination systems can evolve rapidly, but how this occurs, and even how widespread the same mechanism is within a given taxonomic group, remains largely unknown. By experimentally mapping the sex determination architecture in the ant, Vollenhovia emeryi, we found that the well-characterized honey bee sex determination locus originated more than 100 million years ago. However, we also found an additional locus that has no homology to the first. Currently uncharacterized, this locus suggests that different species may use a variety of complementary sex determination mechanisms. Yet, core elements of the complementary sex determination machinery appear to be ancient.

Social evolution. Genomic signatures of evolutionary transitions from solitary to group living

The evolution of eusociality is one of the major transitions in evolution, but the underlying genomic changes are unknown. We compared the genomes of 10 bee species that vary in social complexity, representing multiple independent transitions in social evolution, and report three major findings. First, many important genes show evidence of neutral evolution as a consequence of relaxed selection with increasing social complexity. Second, there is no single road map to eusociality; independent evolutionary transitions in sociality have independent genetic underpinnings. Third, though clearly independent in detail, these transitions do have similar general features, including an increase in constrained protein evolution accompanied by increases in the potential for gene regulation and decreases in diversity and abundance of transposable elements. Eusociality may arise through different mechanisms each time, but would likely always involve an increase in the complexity of gene networks.

The genomes of two key bumblebee species with primitive eusocial organization

Background

The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats.

Results

We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits.

Conclusions

These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0623-3) contains supplementary material, which is available to authorized users.

The evolutionary dynamics of major regulators for sexual development among Hymenoptera species

All hymenopteran species, such as bees, wasps and ants, are characterized by the common principle of haplodiploid sex determination in which haploid males arise from unfertilized eggs and females from fertilized eggs. The underlying molecular mechanism has been studied in detail in the western honey bee Apis mellifera, in which the gene complementary sex determiner (csd) acts as primary signal of the sex determining pathway, initiating female development by csd-heterozygotes. Csd arose from gene duplication of the feminizer (fem) gene, a transformer (tra) ortholog, and mediates in conjunction with transformer2 (tra2) sex-specific splicing of fem. Comparative molecular analyses identified fem/tra and its downstream target doublesex (dsx) as conserved unit within the sex determining pathway of holometabolous insects. In this study, we aim to examine evolutionary differences among these key regulators. Our main hypothesis is that sex determining key regulators in Hymenoptera species show signs of coevolution within single phylogenetic lineages. We take advantage of several newly sequenced genomes of bee species to test this hypothesis using bioinformatic approaches. We found evidences that duplications of fem are restricted to certain bee lineages and notable amino acid differences of tra2 between Apis and non-Apis species propose structural changes in Tra2 protein affecting co-regulatory function on target genes. These findings may help to gain deeper insights into the ancestral mode of hymenopteran sex determination and support the common view of the remarkable evolutionary flexibility in this regulatory pathway.

Independent evolutionary origin of fem paralogous genes and complementary sex determination in hymenopteran insects

The primary signal of sex determination in the honeybee, the complementary sex determiner (csd) gene, evolved from a gene duplication event from an ancestral copy of the fem gene. Recently, other paralogs of the fem gene have been identified in several ant and bumblebee genomes. This discovery and the close phylogenetic relationship of the paralogous gene sequences led to the hypothesis of a single ancestry of the csd genetic system of complementary sex determination in the Hymenopteran insects, in which the fem and csd gene copies evolved as a unit in concert with the mutual transfers of sequences (concerted evolution). Here, we show that the paralogous gene copies evolved repeatedly through independent gene duplication events in the honeybee, bumblebee, and ant lineage. We detected no sequence tracts that would indicate a DNA transfer between the fem and the fem1/csd genes between different ant and bee species. Instead, we found tracts of duplication events in other genomic locations, suggesting that gene duplication was a frequent event in the evolution of these genes. These and other evidences suggest that the fem1/csd gene originated repeatedly through gene duplications in the bumblebee, honeybee, and ant lineages in the last 100 million years. Signatures of concerted evolution were not detectable, implicating that the gene tree based on neutral synonymous sites represents the phylogenetic relationships and origins of the fem and fem1/csd genes. Our results further imply that the fem1 and csd gene in bumblebees, honeybees, and ants are not orthologs, because they originated independently from the fem gene. Hence, the widely shared and conserved complementary sex determination mechanism in Hymenopteran insects is controlled by different genes and molecular processes. These findings highlight the limits of comparative genomics and emphasize the requirement to study gene functions in different species and major hymenopteran lineages.

Phylogenetic distribution and evolutionary dynamics of the sex determination genes doublesex and transformer in insects

Sex determination in insects is characterized by a gene cascade that is conserved at the bottom but contains diverse primary signals at the top. The bottom master switch gene doublesex is found in all insects. Its upstream regulator transformer is present in the orders Hymenoptera, Coleoptera and Diptera, but has thus far not been found in Lepidoptera and in the basal lineages of Diptera. transformer is presumed to be ancestral to the holometabolous insects based on its shared domains and conserved features of autoregulation and sex-specific splicing. We interpret that its absence in basal lineages of Diptera and its order-specific conserved domains indicate multiple independent losses or recruitments into the sex determination cascade. Duplications of transformer are found in derived families within the Hymenoptera, characterized by their complementary sex determination mechanism. As duplications are not found in any other insect order, they appear linked to the haplodiploid reproduction of the Hymenoptera. Further phylogenetic analyses combined with functional studies are needed to understand the evolutionary history of the transformer gene among insects.