Signatures of selection among sex-determining alleles of the honey bee

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.

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.

Variability and Number of Circulating Complementary Sex Determiner (Csd) Alleles in a Breeding Population of Italian Honeybees under Controlled Mating

In Apis mellifera, csd is the primary gene involved in sex determination: haploid hemizygous eggs develop as drones, while females develop from eggs heterozygous for the csd gene. If diploid eggs are homozygous for the csd gene, diploid drones will develop, but will be eaten by worker bees before they are born. Therefore, high csd allelic diversity is a priority for colony survival and breeding. This study aims to investigate the variability of the hypervariable region (HVR) of the csd gene in bees sampled in an apiary under a selection scheme. To this end, an existing dataset of 100 whole-genome sequences was analyzed with a validated pipeline based on de novo assembly of sequences within the HVR region. In total, 102 allelic sequences were reconstructed and translated into amino acid sequences. Among these, 47 different alleles were identified, 44 of which had previously been observed, while 3 are novel alleles. The results show a high variability in the csd region in this breeding population of honeybees.

Zygosity-based sex determination in a butterfly drives hypervariability of Masculinizer

Nature has devised many ways of producing males and females. Here, we report on a previously undescribed mechanism for Lepidoptera that functions without a female-specific gene. The number of alleles or allele heterozygosity in a single Z-linked gene (BaMasc) is the primary sex-determining switch in Bicyclus anynana butterflies. Embryos carrying a single BaMasc allele develop into WZ (or Z0) females, those carrying two distinct alleles develop into ZZ males, while (ZZ) homozygotes initiate female development, have mismatched dosage compensation, and die as embryos. Consequently, selection against homozygotes has favored the evolution of spectacular allelic diversity: 205 different coding sequences of BaMasc were detected in a sample of 246 females. The structural similarity of a hypervariable region (HVR) in BaMasc to the HVR in Apis mellifera csd suggests molecular convergence between deeply diverged insect lineages. Our discovery of this primary switch highlights the fascinating diversity of sex-determining mechanisms and underlying evolutionary drivers.

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.

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.

New insights into the criteria of functional heterozygosity of the Apis mellifera complementary sex determining gene–Discovery of a functional allele pair differing by a single amino acid

The complementary sex determiner (csd) gene is responsible for controlling the sex-determination molecular switch in western honey bees (Apis mellifera): bees that are heterozygous for csd develop into females, whereas bees that are hemizygous or homozygous develop into males. The homozygous diploid males are destroyed at an early stage of their development. It has been proposed that the minimal number of amino acid differences between two csd alleles needed to fully determine femaleness is five and it has also been shown that smaller differences may result in forming an evolutionary intermediate that is not fully capable of female determination, but has increased fitness compared to the homozygous genotype. In this study, we have implemented a terminal restriction length polymorphism-based method of identifying and distinguishing paternal alleles in a given bee colony and assigning them to a particular maternal allele in order to gather information on large number of functional csd pairs and also to identify, to some extent, genotypes that are underrepresented or absent in bee colonies. The main finding of this study is the identification of a fully functional genotype consisting of csd alleles that differed from each other by a one amino acid position. The individuals carrying this genotype expressed only female-specific transcripts of feminizer and double-sex genes. By comparing the sequences differences between the csd pair identified in our study with those described earlier, we conclude that functional heterozygosity of the csd gene is dependent not only on the number of the amino acid differences but also on the sequence context and position of the change. The discovery of a functional allele pair differing by a single amino acid also implies that the generation of a new csd specificity may also occur during a single mutation step with no need for evolutionary intermediates accumulating further mutations.

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.

Highly efficient site-specific integration of DNA fragments into the honeybee genome using CRISPR/Cas9

Functional genetic studies in honeybees have been limited to transposon mediated transformation and site directed mutagenesis tools. However, site- and sequence-specific manipulations that insert DNA fragments or replace sequences at specific target sites are lacking. Such tools would enable the tagging of proteins, the expression of reporters and site-specific amino acid changes, which are all gold standard manipulations for physiological, organismal, and genetic studies. However, such manipulations must be very efficient in honeybees since screening and crossing procedures are laborious due to their social organization. Here, we report an accurate and remarkably efficient site-specific integration of DNA-sequences into the honeybee genome using clustered regularly interspaced short palindromic repeat/clustered regularly interspaced short palindromic repeat-associated protein 9-mediated homology-directed repair. We employed early embryonic injections and selected a highly efficient sgRNA in order to insert 294 and 729 bp long DNA sequences into a specific locus at the dsx gene. These sequences were locus-specifically integrated in 57% and 59% of injected bees. Most importantly, 21% and 25% of the individuals lacked the wildtype sequence demonstrating that we generated homozygous mutants in which all cells are affected (no mosaicism). The highly efficient, locus-specific insertions of nucleotide sequences generating homozygous mutants demonstrate that systematic molecular studies for honeybees are in hand that allow somatic mutation approaches via workers or studies in the next generation using queens with their worker progeny. The employment of early embryonic injections and screenings of highly efficient sgRNAs may offer the prospect of highly successful sequence- and locus-specific mutations also in other organisms.

Application of Next Generation Semiconductor-Based Sequencing for the Identification of Apis mellifera Complementary Sex Determiner (csd) Alleles from Honey DNA

Simple Summary

Honey contains traces of the DNA of the honey bees that produced it. This environmental DNA can therefore be used to investigate the genome of the honey bees. In this study, we used a next generation sequencing technology to analyze the variability of a key gene of Apis mellifera L., the complementary sex determiner (csd) gene, using honey environmental DNA as a source of honey bee DNA. This gene determines the sex of the bees. Two different alleles at this locus are needed to produce females whereas males have only one copy of this gene as they are haploid. In case two identical alleles are present in a diploid individual, the larvae are not vital and are discarded by the workers. Therefore, there is an advantage in maintaining a large csd diversity in honey bee populations. In light of the recent decline in honey bee populations, it is important to monitor the allele variability at this gene. The applied methodology provided a new strategy to disclose the genetic diversity at the csd gene at the population-wide level and identify most, if not all, csd alleles present in the colonies in a single analysis.

Abstract

The complementary sex determiner (csd) gene plays an essential role in the sex determination of Apis mellifera L. Females develop only if fertilized eggs have functional heterozygous genotypes at this gene whereas males, being haploids, are hemizygous. Two identical csd alleles produce non vital males. In light of the recent decline in honey bee populations, it is therefore important to monitor the allele variability at this gene. In this study, we tested the application of next generation semiconductor-based sequencing technology (Ion Torrent) coupled with environmental honey DNA as a source of honey bee genome information to retrieve massive sequencing data for the analysis of variability at the hypervariable region (HVR) of the csd gene. DNA was extracted from 12 honey samples collected from honeycombs directly retrieved from 12 different colonies. A specifically designed bioinformatic pipeline, applied to analyze a total of about 1.5 million reads, identified a total of 160 different csd alleles, 55% of which were novel. The average number of alleles per sample was compatible with the number of expected patrilines per colony, according to the mating behavior of the queens. Allele diversity at the csd could also provide information useful to reconstruct the history of the honey.

A simulation study of a honeybee breeding scheme accounting for polyandry, direct and maternal effects on colony performance

BACKGROUND

Efficient breeding programs are difficult to implement in honeybees due to their biological specificities (polyandry and haplo-diploidy) and complexity of the traits of interest, with performances being measured at the colony scale and resulting from the joint effects of tens of thousands of workers (called direct effects) and of the queen (called maternal effects). We implemented a Monte Carlo simulation program of a breeding plan designed specifically for Apis mellifera's populations to assess the impact of polyandry versus monoandry on colony performance, inbreeding level and genetic gain depending on the individual selection strategy considered, i.e. complete mass selection or within-family (maternal lines) selection. We simulated several scenarios with different parameter setups by varying initial genetic variances and correlations between direct and maternal effects, the selection strategy and the polyandry level. Selection was performed on colony phenotypes.

RESULTS

All scenarios showed strong increases in direct breeding values of queens after 20 years of selection. Monoandry led to significantly higher direct than maternal genetic gains, especially when a negative correlation between direct and maternal effects was simulated. However, the relative increase in these genetic gains depended also on their initial genetic variability and on the selection strategy. When polyandry was simulated, the results were very similar with either 8 or 16 drones mated to each queen. Across scenarios, polyandrous mating resulted in equivalent or higher gains in performance than monoandrous mating, but with considerably lower inbreeding rates. Mass selection conferred a ~ 20% increase in performance compared to within-family selection, but was also accompanied by a strong increase in inbreeding levels (25 to 50% higher).

CONCLUSIONS

Our study is the first to compare the long-term effects of polyandrous versus monoandrous mating in honeybee breeding. The latter is an emergent strategy to improve specific traits, such as resistance to varroa, which can be difficult or expensive to phenotype. However, if used during several generations in a closed population, monoandrous mating increases the inbreeding level of queens much more than polyandrous mating, which is a strong limitation of this strategy.

Possible Epigenetic Origin of a Recurrent Gynandromorph Pattern in Megachile Wild Bees

Simple Summary

Gynandromorphs, i.e., individuals with a mix of male and female body parts, are known for many species of insects and other animals with separate sexes. This anomaly is generally regarded as the result of localized genetic mutations in sex-determining genes. We analyzed the specific mix of male and female characters in naturally occurring gynandromorphs of 21 species of the wild bee genus Megachile and found a recurrent pattern. Based on the regularity of this pattern, and the current knowledge on sex determination and sex differentiation in the relatively closely-related honey bee, we argue that the origin of these composite phenotypes is possibly epigenetic, rather than genetic, i.e., produced by some defects in the maintenance of the regulatory signals that control sex differentiation at the level of single cell lineages, rather than triggered by genetic mutations.

Abstract

Gynandromorphs, i.e., individuals with a mix of male and female traits, are common in the wild bees of the genus Megachile (Hymenoptera, Apoidea). We described new transverse gynandromorphs in Megachile pilidens Alfkeen, 1924 and analyze the spatial distribution of body parts with male vs. female phenotype hitherto recorded in the transverse gynandromorphs of the genus Megachile. We identified 10 different arrangements, nine of which are minor variants of a very general pattern, with a combination of male and female traits largely shared by the gynandromorphs recorded in 20 out of 21 Megachile species in our dataset. Based on the recurrence of the same gynandromorph pattern, the current knowledge on sex determination and sex differentiation in the honey bee, and the results of recent gene-knockdown experiments in these insects, we suggest that these composite phenotypes are possibly epigenetic, rather than genetic, mosaics, with individual body parts of either male or female phenotype according to the locally expressed product of the alternative splicing of sex-determining gene transcripts.

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

Global allele polymorphism indicates a high rate of allele genesis at a locus under balancing selection

When selection favours rare alleles over common ones (balancing selection in the form of negative frequency-dependent selection), a locus may maintain a large number of alleles, each at similar frequency. To better understand how allelic richness is generated and maintained at such loci, we assessed 201 sequences of the complementary sex determiner (csd) of the Asian honeybee (Apis cerana), sampled from across its range. Honeybees are haplodiploid; hemizygotes at csd develop as males and heterozygotes as females, while homozygosity is lethal. Thus, csd is under strong negative frequency-dependent selection because rare alleles are less likely to end up in the lethal homozygous form. We find that in A. cerana, as in other Apis, just a few amino acid differences between csd alleles in the hypervariable region are sufficient to trigger female development. We then show that while allelic lineages are spread across geographical regions, allelic differentiation is high between populations, with most csd alleles (86.3%) detected in only one sample location. Furthermore, nucleotide diversity in the hypervariable region indicates an excess of recently arisen alleles, possibly associated with population expansion across Asia since the last glacial maximum. Only the newly invasive populations of the Austral-Pacific share most of their csd alleles. In all, the geographic patterns of csd diversity in A. cerana indicate that high mutation rates and balancing selection act together to produce high rates of allele genesis and turnover at the honeybee sex locus, which in turn leads to its exceptionally high local and global polymorphism.

An Alternative, High Throughput Method to Identify Csd Alleles of the Honey Bee

Applying instrumental insemination in closely related honey bee colonies often leads to frequent lethality of offspring causing colony collapse. This is due to the peculiarities of honey bee reproductive biology, where the complementary sex determination (csd) gene drives sex determination within a haplodiploid system. Diploid drones containing homozygous genotypes are lethal. Tracking of csd alleles using molecular markers prevents this unwanted event in closed breeding programs. Our approach described here is based on high throughput sequencing (HTS) that provides more data than traditional molecular techniques and is capable of analysing sources containing multiple alleles, including diploid individuals as the bee queen. The approach combines HTS technique and clipping wings as a minimally invasive method to detect the complementary sex determiner (csd) alleles directly from honey bee queens. Furthermore, it might also be suitable for screening alleles of honey harvested from hives of a closed breeding facility. Data on alleles of the csd gene from different honey bee subspecies are provided. It might contribute to future databases that could potentially be used to track the origin of honey. With the help of tracking csd alleles, more focused crossings will be possible, which could in turn accelerate honey bee breeding programmes targeting increase tolerance against varroosis as well.

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.

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.

A genetic switch for worker nutrition-mediated traits in honeybees

Highly social insects are characterized by caste dimorphism, with distinct size differences of reproductive organs between fertile queens and the more or less sterile workers. An abundance of nutrition or instruction via diet-specific compounds has been proposed as explanations for the nutrition-driven queen and worker polyphenism. Here, we further explored these models in the honeybee (Apis mellifera) using worker nutrition rearing and a novel mutational screening approach using the clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) method. The worker nutrition-driven size reduction of reproductive organs was restricted to the female sex, suggesting input from the sex determination pathway. Genetic screens on the sex determination genes in genetic females for size polyphenism revealed that doublesex (dsx) mutants display size-reduced reproductive organs irrespective of the sexual morphology of the organ tissue. In contrast, feminizer (fem) mutants lost the response to worker nutrition-driven size control. The first morphological worker mutants in honeybees demonstrate that the response to nutrition relies on a genetic program that is switched “ON” by the fem gene. Thus, the genetic instruction provided by the fem gene provides an entry point to genetically dissect the underlying processes that implement the size polyphenism.

In honeybees, nutrition drives dimorphic size development of reproductive organs in fertile queens and sterile workers. A study using the first induced morphological mutants in honeybees demonstrates that this developmental plasticity requires a genetic program that is switched on by the “feminizer” gene.

In honeybees, nutrition drives dimorphic size development of reproductive organs in fertile queens and sterile workers. The first induced morphological mutants in honeybees demonstrate that this developmental plasticity requires a genetic program that is switched “ON” by the feminizer (fem) gene.

Extreme polyandry aids the establishment of invasive populations of a social insect

Although monandry is believed to have facilitated the evolution of eusociality, many highly eusocial insects have since evolved extreme polyandry. The transition to extreme polyandry was likely driven by the benefits of within-colony genetic variance to task specialization and/or disease resistance, but the extent to which it confers secondary benefits, once evolved, is unclear. Here we investigate the consequences of extreme polyandry on the invasive potential of the Asian honey bee, Apis cerana. In honey bees and other Hymenoptera, small newly founded invasive populations must overcome the genetic constraint of their sex determination system that requires heterozygosity at a sex-determining locus to produce viable females. We find A. cerana queens in an invasive population mate with an average of 27 males (range 16-42) that would result in the founding queen/s carrying 75% of their source population's sex alleles in stored sperm. This mating frequency is similar to native-range Chinese A. cerana (mean 29 males, range 19-46). Simulations reveal that extreme polyandry reduces the risk, relative to monandry or moderate polyandry, that colonies produce a high incidence of inviable brood in populations that have experienced a founder event, that is, when sex allele diversity is low and/or allele frequencies are unequal. Thus, extreme polyandry aids the invasiveness of A. cerana in two ways: (1) by increasing the sex locus allelic richness carried to new populations with each founder, thereby increasing sex locus heterozygosity; and (2) by reducing the population variance in colony fitness following a founder event.

Contrasting patterns of nucleotide polymorphism suggest different selective regimes within different parts of the PgiC1 gene in Festuca ovina L

Background

Phosphoglucose isomerase (PGI, EC 5.3.1.9) is an essential metabolic enzyme in all eukaryotes. An earlier study of the PgiC1 gene, which encodes cytosolic PGI in the grass Festuca ovina L., revealed a marked difference in the levels of nucleotide polymorphism between the 5’ and 3’ portions of the gene.

Methods

In the present study, we characterized the sequence polymorphism in F. ovina PgiC1 in more detail and examined possible explanations for the non-uniform pattern of nucleotide polymorphism across the gene.

Results

Our study confirms that the two portions of the PgiC1 gene show substantially different levels of DNA polymorphism and also suggests that the peptide encoded by the 3’ portion of PgiC1 is functionally and structurally more important than that encoded by the 5’ portion. Although there was some evidence of purifying selection (d_N/d_S test) on the 5’ portion of the gene, the signature of purifying selection was considerably stronger on the 3’ portion of the gene (d_N/d_S and McDonald–Kreitman tests). Several tests support the action of balancing selection within the 5’ portion of the gene. Wall’s B and Q tests were significant only for the 5’ portion of the gene. There were also marked peaks of nucleotide diversity, Tajima’s D and the d_N/d_S ratio at or around a PgiC1 codon site (within the 5’ portion of the gene) that a previous study had suggested was subject to positive diversifying selection.

Conclusions

Our results suggest that the two portions of the gene have been subject to different selective regimes. Purifying selection appears to have been the main force contributing to the relatively low level of polymorphism within the 3’ portion of the sequence. In contrast, it is possible that balancing selection has contributed to the maintenance of the polymorphism within the 5’ portion of the gene.

Electronic supplementary material

The online version of this article (doi:10.1186/s41065-017-0032-6) contains supplementary material, which is available to authorized users.

An invasive social insect overcomes genetic load at the sex locus

Some invasive hymenopteran social insects found new populations with very few reproductive individuals. This is despite the high cost of founder effects for such insects, which generally require heterozygosity at a single locus-the complementary sex determiner, csd-to develop as females. Individuals that are homozygous at csd develop as either infertile or subfertile diploid males or not at all. Furthermore, diploid males replace the female workers that are essential for colony function. Here we document how the Asian honey bee (Apis cerana) overcame the diploid male problem during its invasion of Australia. Natural selection prevented the loss of rare csd alleles due to genetic drift and corrected the skew in allele frequencies caused by founder effects to restore high average heterozygosity. Thus, balancing selection can alleviate the genetic load at csd imposed by severe bottlenecks, and so facilitate invasiveness.

Similar but not the same: insights into the evolutionary history of paralogous sex-determining genes of the dwarf honey bee Apis florea

Studying the fate of duplicated genes provides informative insight into the evolutionary plasticity of biological pathways to which they belong. In the paralogous sex-determining genes complementary sex determiner (csd) and feminizer (fem) of honey bee species (genus Apis), only heterozygous csd initiates female development. Here, the full-length coding sequences of the genes csd and fem of the phylogenetically basal dwarf honey bee Apis florea are characterized. Compared with other Apis species, remarkable evolutionary changes in the formation and localization of a protein-interacting (coiled-coil) motif and in the amino acids coding for the csd characteristic hypervariable region (HVR) are observed. Furthermore, functionally different csd alleles were isolated as genomic fragments from a random population sample. In the predicted potential specifying domain (PSD), a high ratio of πN/πS=1.6 indicated positive selection, whereas signs of balancing selection, commonly found in other Apis species, are missing. Low nucleotide diversity on synonymous and genome-wide, non-coding sites as well as site frequency analyses indicated a strong impact of genetic drift in A. florea, likely linked to its biology. Along the evolutionary trajectory of ~30 million years of csd evolution, episodic diversifying selection seems to have acted differently among distinct Apis branches. Consistently low amino-acid differences within the PSD among pairs of functional heterozygous csd alleles indicate that the HVR is the most important region for determining allele specificity. We propose that in the early history of the lineage-specific fem duplication giving rise to csd in Apis, A. florea csd stands as a remarkable example for the plasticity of initial sex-determining signals.

Molecular tools and bumble bees: revealing hidden details of ecology and evolution in a model system

Bumble bees are a longstanding model system for studies on behaviour, ecology and evolution, due to their well-studied social lifestyle, invaluable role as wild and managed pollinators, and ubiquity and diversity across temperate ecosystems. Yet despite their importance, many aspects of bumble bee biology have remained enigmatic until the rise of the genetic and, more recently, genomic eras. Here, we review and synthesize new insights into the ecology, evolution and behaviour of bumble bees that have been gained using modern genetic and genomic techniques. Special emphasis is placed on four areas of bumble bee biology: the evolution of eusociality in this group, population-level processes, large-scale evolutionary relationships and patterns, and immunity and resistance to pesticides. We close with a prospective on the future of bumble bee genomics research, as this rapidly advancing field has the potential to further revolutionize our understanding of bumble bees, particularly in regard to adaptation and resilience. Worldwide, many bumble bee populations are in decline. As such, throughout the review, connections are drawn between new molecular insights into bumble bees and our understanding of the causal factors involved in their decline. Ongoing and potential applications to bumble bee management and conservation are also included to demonstrate how genetics- and genomics-enabled research aids in the preservation of this threatened group.

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.

Highly efficient integration and expression of piggyBac-derived cassettes in the honeybee (Apis mellifera)

Honeybees (Apis mellifera), which are important pollinators of plants, display remarkable individual behaviors that collectively contribute to the organization of a complex society. Advances in dissecting the complex processes of honeybee behavior have been limited in the recent past due to a lack of genetic manipulation tools. These tools are difficult to apply in honeybees because the unit of reproduction is the colony, and many interesting phenotypes are developmentally specified at later stages. Here, we report highly efficient integration and expression of piggyBac-derived cassettes in the honeybee. We demonstrate that 27 and 20% of queens stably transmitted two different expression cassettes to their offspring, which is a 6- to 30-fold increase in efficiency compared with those generally reported in other insect species. This high efficiency implies that an average beekeeping facility with a limited number of colonies can apply this tool. We demonstrated that the cassette stably and efficiently expressed marker genes in progeny under either an artificial or an endogenous promoter. This evidence of efficient expression encourages the use of this system to inhibit gene functions through RNAi in specific tissues and developmental stages by using various promoters. We also showed that the transgenic marker could be used to select transgenic offspring to be employed to facilitate the building of transgenic colonies via the haploid males. We present here the first to our knowledge genetic engineering tool that will efficiently allow for the systematic detection and better understanding of processes underlying the biology of honeybees.

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.

Sex determination in insects: variations on a common theme

Recent studies in a representative selection of holometabolous insects suggest that, despite diversity at the instructive level, the signal-relaying part of the sex-determining pathway is remarkably well conserved. In principle, it is composed of the transformer gene (tra), which acts as a common binary switch that transduces the selected sexual fate, female when ON, male when OFF, to the downstream effector doublesex(dsx) that controls overt sexual differentiation. An interesting recurrent feature is that tra is switched ON in the early zygote by maternally provisioned tra activity. Different male-determining signals evolved, which prevent maternal activation of zygotic tra to allow for male development. In some species, where lack of maternal activation leaves tra in the OFF state, novel female-determining signals were deployed to activate zygotic tra. It appears that both the instructive end of the pathway upstream of tra as well as the executive end downstream of dsx are primary targets of evolutionary divergence, while the transduction part seems less prone to changes. We propose that this is a feature shared with many other signaling cascades that regulate developmental fates.

Gradual molecular evolution of a sex determination switch through incomplete penetrance of femaleness

Some genes regulate phenotypes that are either present or absent. They are often important regulators of developmental switches and are involved in morphological evolution. We have little understanding of the molecular mechanisms by which these absence/presence gene functions have evolved, because the phenotype and fitness of molecular intermediate forms are unknown. Here, we studied the sex-determining switch of 14 natural sequence variants of the csd gene among 76 genotypes of the honeybee (Apis mellifera). Heterozygous genotypes (different specificities) of the csd gene determine femaleness, while hemizygous genotypes (single specificity) determine maleness. Homozygous genotypes of the csd gene (same specificity) are lethal. We found that at least five amino acid differences and length variation between Csd specificities in the specifying domain (PSD) were sufficient to regularly induce femaleness. We estimated that, on average, six pairwise amino acid differences evolved under positive selection. We also identified a natural evolutionary intermediate that showed only three amino acid length differences in the PSD relative to its parental allele. This genotype showed an intermediate fitness because it implemented lethality regularly and induced femaleness infrequently (i.e., incomplete penetrance). We suggest incomplete penetrance as a mechanism through which new molecular switches can gradually and adaptively evolve.

Nucleotide variability at its limit? Insights into the number and evolutionary dynamics of the sex-determining specificities of the honey bee Apis mellifera

Deciphering the evolutionary processes driving nucleotide variation in multiallelic genes is limited by the number of genetic systems in which such genes occur. The complementary sex determiner (csd) gene in the honey bee Apis mellifera is an informative example for studying allelic diversity and the underlying evolutionary forces in a well-described model of balancing selection. Acting as the primary signal of sex determination, diploid individuals heterozygous for csd develop into females, whereas csd homozygotes are diploid males that have zero fitness. Examining 77 of the functional heterozygous csd allele pairs, we established a combinatorical criteria that provide insights into the minimum number of amino acid differences among those pairs. Given a data set of 244 csd sequences, we show that the total number of csd alleles found in A. mellifera ranges from 53 (locally) to 87 (worldwide), which is much higher than was previously reported (20). Using a coupon-collector model, we extrapolate the presence of in total 116–145 csd alleles worldwide. The hypervariable region (HVR) is of particular importance in determining csd allele specificity, and we provide for this region evidence of high evolutionary rate for length differences exceeding those of microsatellites. The proportion of amino acids driven by positive selection and the rate of nonsynonymous substitutions in the HVR-flanking regions reach values close to 1 but differ with respect to the HVR length. Using a model of csd coalescence, we identified the high originating rate of csd specificities as a major evolutionary force, leading to an origin of a novel csd allele every 400,000 years. The csd polymorphism frequencies in natural populations indicate an excess of new mutations, whereas signs of ancestral transspecies polymorphism can still be detected. This study provides a comprehensive view of the enormous diversity and the evolutionary forces shaping a multiallelic gene.

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

The transformer (tra) gene is a key regulator in the signalling hierarchy controlling all aspects of somatic sexual differentiation in Drosophila and other insects. Here, we show that six of the seven sequenced ants have two copies of tra. Surprisingly, the two paralogues are always more similar within species than among species. Comparative sequence analyses indicate that this pattern is owing to the ongoing concerted evolution after an ancestral duplication rather than independent duplications in each of the six species. In particular, there was strong support for inter-locus recombination between the paralogues of the ant Atta cephalotes. In the five species where the location of paralogues is known, they are adjacent to each other in four cases and separated by only few genes in the fifth case. Because there have been extensive genomic rearrangements in these lineages, this suggests selection acting to conserve their synteny. In three species, we also find a signature of positive selection in one of the paralogues. In three bee species where information is available, the tra gene is also duplicated, the copies are adjacent and in at least one species there was recombination between paralogues. These results suggest that concerted evolution plays an adaptive role in the evolution of this gene family.

The Am-tra2 gene is an essential regulator of female splice regulation at two levels of the sex determination hierarchy of the honeybee

Heteroallelic and homo- or hemiallelic Complementary sex determiner (Csd) proteins determine sexual fate in the honeybee (Apis mellifera) by controlling the alternative splicing of the downstream gene fem (feminizer). Thus far, we have little understanding of how heteroallelic Csd proteins mediate the splicing of female fem messenger RNAs (mRNAs) or how Fem proteins direct the splicing of honeybee dsx (Am-dsx) pre-mRNAs. Here, we report that Am-tra2, which is an ortholog of Drosophila melanogaster tra2, is an essential component of female splicing of the fem and Am-dsx transcripts in the honeybee. The Am-tra2 transcripts are alternatively (but non-sex-specifically) spliced, and they are translated into six protein isoforms that all share the basic RNA-binding domain/RS (arginine/serine) domain structure. Knockdown studies showed that the Am-tra2 gene is required to splice fem mRNAs into the productive female and nonproductive male forms. We suggest that the Am-Tra2 proteins are essential regulators of fem pre-mRNA splicing that, together with heteroallelic Csd proteins and/or Fem proteins, implement the female pathway. In males, the Am-Tra2 proteins may enhance the switch of fem transcripts into the nonproductive male form when heteroallelic Csd proteins are absent. This dual function of Am-Tra2 proteins possibly enhances and stabilizes the binary decision process of male/female splicing. Our knockdown studies also imply that the Am-Tra2 protein is an essential regulator for Am-dsx female splice regulation, suggesting an ancestral role in holometabolous insects. We also provide evidence that the Am-tra2 gene has an essential function in honeybee embryogenesis that is unrelated to sex determination.

The sex determination gene shows no founder effect in the giant honey bee, Apis dorsata

Background

All honey bee species (Apis spp) share the same sex determination mechanism using the complementary sex determination (csd) gene. Only individuals heterogeneous at the csd allele develop into females, and the homozygous develop into diploid males, which do not survive. The honeybees are therefore under selection pressure to generate new csd alleles. Previous studies have shown that the csd gene is under balancing selection. We hypothesize that due to the long separation from the mainland of Hainan Island, China, that the giant honey bees (Apis dorsata) should show a founder effect for the csd gene, with many different alleles clustered together, and these would be absent on the mainland.

Methodology/Principal Findings

We sampled A. dorsata workers from both Hainan and Guangxi Provinces and then cloned and sequenced region 3 of the csd gene and constructed phylogenetic trees. We failed to find any clustering of the csd alleles according to their geographical origin, i.e. the Hainan and Guangxi samples did not form separate clades. Further analysis by including previously published csd sequences also failed to show any clade-forming in both the Philippines and Malaysia.

Conclusions/Significance

Results from this study and those from previous studies did not support the expectations of a founder effect. We conclude that because of the extremely high mating frequency of A. dorsata queens, a founder effect does not apply in this species.

Transcript levels of ten caste-related genes in adult diploid males of Melipona quadrifasciata (Hymenoptera, Apidae) - A comparison with haploid males, queens and workers

In Hymenoptera, homozygosity at the sex locus results in the production of diploid males. In social species, these pose a double burden by having low fitness and drawing resources normally spent for increasing the work force of a colony. Yet, diploid males are of academic interest as they can elucidate effects of ploidy (normal males are haploid, whereas the female castes, the queens and workers, are diploid) on morphology and life history. Herein we investigated expression levels of ten caste-related genes in the stingless bee Melipona quadrifasciata, comparing newly emerged and 5-day-old diploid males with haploid males, queens and workers. In diploid males, transcript levels for dunce and paramyosin were increased during the first five days of adult life, while those for diacylglycerol kinase and the transcriptional co-repressor groucho diminished. Two general trends were apparent, (i) gene expression patterns in diploid males were overall more similar to haploid ones and workers than to queens, and (ii) in queens and workers, more genes were up-regulated after emergence until day five, whereas in diploid and especially so in haploid males more genes were down-regulated. This difference between the sexes may be related to longevity, which is much longer in females than in males.

Function and evolution of sex determination mechanisms, genes and pathways in insects

Animals have evolved a bewildering diversity of mechanisms to determine the two sexes. Studies of sex determination genes – their history and function – in non-model insects and Drosophila have allowed us to begin to understand the generation of sex determination diversity. One common theme from these studies is that evolved mechanisms produce activities in either males or females to control a shared gene switch that regulates sexual development. Only a few small-scale changes in existing and duplicated genes are sufficient to generate large differences in sex determination systems. This review summarises recent findings in insects, surveys evidence of how and why sex determination mechanisms can change rapidly and suggests fruitful areas of future research.

Polymorphism analysis of csd gene in six Apis mellifera subspecies

The complementary sex determination (csd) gene is the primary gene determining the gender of honey bees (Apis spp). In this study we analyzed the polymorphism of csd gene in six Apis mellifera subspecies. The genomic region 3 of csd gene in these six A. mellifera was cloned, and identified. A total of 79 haplotypes were obtained from these six subspecies. Analysis showed that region 3 of csd gene has a high level of polymorphism in all the six A. mellifera subspecies. The A. m. anatolica subspecies has a slightly higher nucleotide diversity (π) than other subspecies, while the π values showed no significant difference among the other five subspecies. The phylogenetic tree showed that all the csd haplotypes from different A. mellifera subspecies are scattered throughout the tree, without forming six different clades. Population differentiation analysis showed that there are significant genetic differentiations among some of the subspecies. The NJ phylogenetic tree showed that the A. m. caucasica and A. m. carnica have the closest relationship, followed by A. m. ssp, A. m. ligustica, A. m. carpatica and A. m. anatolica that were gathered in the tree in turn.

A thelytokous lineage of socially parasitic honey bees has retained heterozygosity despite at least 10 years of inbreeding

The honey bee population of South Africa is divided into two subspecies: a northern population in which queenless workers reproduce arrhenotokously and a southern one in which workers reproduce thelytokously. A hybrid zone separates the two, but on at least three occasions the northern population has become infested by reproductive workers derived from the southern population. These parasitic workers lay in host colonies parthenogenetically, resulting in yet more parasites. The current infestation is 20-year old--surprising because an asexual lineage is expected to show a decline in vigor over time due to increasing homozygosity. The decline is expected to be acute in honey bees, where homozygosity at the sex locus is lethal. We surveyed colonies from the zone of infestation and genotyped putative parasites at two sets of linked microsatellite loci. We confirm that there is a single clonal lineage of parasites that shows minor variations arising from recombination events. The lineage shows high levels of heterozygosity, which may be maintained by selection against homozygotes, or by a reduction in recombination frequency within the lineage. We suggest that the clonal lineage can endure the costs of asexual reproduction because of the fitness benefits of its parasitic life history.

Origin of a function by tandem gene duplication limits the evolutionary capability of its sister copy

The most remarkable outcome of a gene duplication event is the evolution of a novel function. Little information exists on how the rise of a novel function affects the evolution of its paralogous sister gene copy, however. We studied the evolution of the feminizer (fem) gene from which the gene complementary sex determiner (csd) recently derived by tandem duplication within the honey bee (Apis) lineage. Previous studies showed that fem retained its sex determination function, whereas the rise of csd established a new primary signal of sex determination. We observed a specific reduction of nonsynonymous to synonymous substitution ratios in Apis to non-Apis fem. We found a contrasting pattern at two other genetically linked genes, suggesting that hitchhiking effects to csd, the locus under balancing selection, is not the cause of this evolutionary pattern. We also excluded higher synonymous substitution rates by relative rate testing. These results imply that stronger purifying selection is operating at the fem gene in the presence of csd. We propose that csd's new function interferes with the function of Fem protein, resulting in molecular constraints and limited evolvability of fem in the Apis lineage. Elevated silent nucleotide polymorphism in fem relative to the genome-wide average suggests that genetic linkage to the csd gene maintained more nucleotide variation in today's population. Our findings provide evidence that csd functionally and genetically interferes with fem, suggesting that a newly evolved gene and its functions can limit the evolutionary capability of other genes in the genome.

Sex determination in honeybees: two separate mechanisms induce and maintain the female pathway

Sex determination in honeybees is realized by the csd and the fem gene that establish and maintain, throughout development, sexual fates via the control of alternative splicing.

Organisms have evolved a bewildering diversity of mechanisms to generate the two sexes. The honeybee (Apis mellifera) employs an interesting system in which sex is determined by heterozygosity at a single locus (the Sex Determination Locus) harbouring the complementary sex determiner (csd) gene. Bees heterozygous at Sex Determination Locus are females, whereas bees homozygous or hemizygous are males. Little is known, however, about the regulation that links sex determination to sexual differentiation. To investigate the control of sexual development in honeybees, we analyzed the functions and the regulatory interactions of genes involved in the sex determination pathway. We show that heterozygous csd is only required to induce the female pathway, while the feminizer (fem) gene maintains this decision throughout development. By RNAi induced knockdown we show that the fem gene is essential for entire female development and that the csd gene exclusively processes the heterozygous state. Fem activity is also required to maintain the female determined pathway throughout development, which we show by mosaic structures in fem-repressed intersexuals. We use expression of Fem protein in males to demonstrate that the female maintenance mechanism is controlled by a positive feedback splicing loop in which Fem proteins mediate their own synthesis by directing female fem mRNA splicing. The csd gene is only necessary to induce this positive feedback loop in early embryogenesis by directing splicing of fem mRNAs. Finally, fem also controls the splicing of Am-doublesex transcripts encoding conserved male- and female-specific transcription factors involved in sexual differentiation. Our findings reveal how the sex determination process is realized in honeybees differing from Drosophila melanogaster.

Sexual differentiation is a fundamental process in the animal kingdom, and different species have evolved a bewildering diversity of mechanisms to generate the two sexes in the proper proportions. Sex determination in honeybees (Apis mellifera) provides an interesting and unusual system to study, as it is governed by heterozygosity of a single locus harbouring the complementary sex determiner gene (csd), in contrast to the well-studied sex chromosome system of Drosophila melanogaster. We show that the female sex determination pathway is exclusively induced by the csd gene in early embryogenesis. Later on and throughout development this inductive signal is maintained via a positive feedback loop of the feminizer (fem) gene, in which the Fem protein mediates its own synthesis. The findings reveal how the sex determination process in honeybees is realized by the regulation and function of two genes differing from Drosophila.

Experimental support for multiple-locus complementary sex determination in the parasitoid Cotesia vestalis

Despite its fundamental role in development, sex determination is highly diverse among animals. Approximately 20% of all animals are haplodiploid, with haploid males and diploid females. Haplodiploid species exhibit diverse but poorly understood mechanisms of sex determination. Some hymenopteran insect species exhibit single-locus complementary sex determination (sl-CSD), where heterozygosity at a polymorphic sex locus initiates female development. Diploid males are homozygous at the sex locus and represent a genetic load because they are inviable or sterile. Inbreeding depression associated with CSD is therefore expected to select for other modes of sex determination resulting in fewer or no diploid males. Here, we investigate an alternative, heretofore hypothetical, mode of sex determination: multiple-locus CSD (ml-CSD). Under ml-CSD, diploid males are predicted to develop only from zygotes that are homozygous at all sex loci. We show that inbreeding for eight generations in the parasitoid wasp Cotesia vestalis leads to increasing proportions of diploid males, a pattern that is consistent with ml-CSD but not sl-CSD. The proportion of diploid males (0.27 +/- 0.036) produced in the first generation of inbreeding (mother-son cross) suggests that two loci are likely involved. We also modeled diploid male production under CSD with three linked loci. Our data visually resemble CSD with linked loci because diploid male production in the second generation was lower than that in the first. To our knowledge, our data provide the first experimental support for ml-CSD.

Evolutionary genetics: changed sex determination in honeybees

In several insects and fish, and probably some mammals, the gene controlling the male-female switch has changed during evolution. It now seems that this has also happened in honeybees, where the sex-determining gene has now been shown to be a duplicate of another Hymenopteran sex-determining gene.

Thelytokous parthenogenesis in unmated queen honeybees (Apis mellifera capensis): central fusion and high recombination rates

The subspecies of honeybee indigenous to the Cape region of South Africa, Apis mellifera capensis, is unique because a high proportion of unmated workers can lay eggs that develop into females via thelytokous parthenogenesis involving central fusion of meiotic products. This ability allows pseudoclonal lineages of workers to establish, which are presently widespread as reproductive parasites within the honeybee populations of South Africa. Successful long-term propagation of a parthenogen requires the maintenance of heterozygosity at the sex locus, which in honeybees must be heterozygous for the expression of female traits. Thus, in successful lineages of parasitic workers, recombination events are reduced by an order of magnitude relative to meiosis in queens of other honeybee subspecies. Here we show that in unmated A. m. capensis queens treated to induce oviposition, no such reduction in recombination occurs, indicating that thelytoky and reduced recombination are not controlled by the same gene. Our virgin queens were able to lay both arrhenotokous male-producing haploid eggs and thelytokous female-producing diploid eggs at the same time, with evidence that they have some voluntary control over which kind of egg was laid. If so, they are able to influence the kind of second-division meiosis that occurs in their eggs post partum.

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

Sex determination in honeybees (Apis mellifera) is governed by heterozygosity at a single locus harbouring the complementary sex determiner (csd) gene, in contrast to the well-studied sex chromosome system of Drosophila melanogaster. Bees heterozygous at csd are females, whereas homozygotes and hemizygotes (haploid individuals) are males. Although at least 15 different csd alleles are known among natural bee populations, the mechanisms linking allelic interactions to switching of the sexual development programme are still obscure. Here we report a new component of the sex-determining pathway in honeybees, encoded 12 kilobases upstream of csd. The gene feminizer (fem) is the ancestrally conserved progenitor gene from which csd arose and encodes an SR-type protein, harbouring an Arg/Ser-rich domain. Fem shares the same arrangement of Arg/Ser- and proline-rich-domain with the Drosophila principal sex-determining gene transformer (tra), but lacks conserved motifs except for a 30-amino-acid motif that Fem shares only with Tra of another fly, Ceratitis capitata. Like tra, the fem transcript is alternatively spliced. The male-specific splice variant contains a premature stop codon and yields no functional product, whereas the female-specific splice variant encodes the functional protein. We show that RNA interference (RNAi)-induced knockdowns of the female-specific fem splice variant result in male bees, indicating that the fem product is required for entire female development. Furthermore, RNAi-induced knockdowns of female allelic csd transcripts result in the male-specific fem splice variant, suggesting that the fem gene implements the switch of developmental pathways controlled by heterozygosity at csd. Comparative analysis of fem and csd coding sequences from five bee species indicates a recent origin of csd in the honeybee lineage from the fem progenitor and provides evidence for positive selection at csd accompanied by purifying selection at fem. The fem locus in bees uncovers gene duplication and positive selection as evolutionary mechanisms underlying the origin of a novel sex determination pathway.