Csd Alleles Research Articles

Genetic diversity of Apis melliferacomplementary sex determiner (csd) in four Australian breeding populations

Honey bee (Apis mellifera) breeding has intensified as populations experience increasing stress and pollination demand increases. Breeding programmes risk genetic diversity losses as mating is focused on a small group of individuals. Loss of diversity at the complementary sex determiner (csd) locus results in decreased viability and reduced honey production. Bees that are homozygous at csd become inviable males rather than workers. We examined csd diversity in four Australian breeding populations: a queen bee breeder from New South Wales, another from Queensland, a Western Australian breeding programme involving 11 bee breeders, and a research population at the NSW Department of Primary Industries. We found 82 unique csd alleles across the four populations, 16 of which have not been previously reported. This study provides a snapshot of csd diversity in Australia which will be useful for the national honey bee genetic improvement programme (Plan Bee).

The Potential of Complementary Sex-Determiner Gene Allelic Diversity for Studying the Number of Patrilines within Honeybee Colonies

Polyandry, a fundamental aspect of honeybee biology, influences genetic diversity within bee colonies. The Csd gene, responsible for sex determination, exhibits a high intraspecific polymorphism, making it a promising candidate for studying patrilines. This study investigates the potential of the complementary sex-determiner gene (Csd) as a marker for genetic studies of honeybee colonies, comparing its efficacy with standard microsatellite markers. A worker bee brood from five colonies was genotyped using both Csd and microsatellite markers. The results showed that Csd alleles offer higher variability than microsatellite markers, providing a superior resolution in genotyping. The number of distinct Csd alleles in each colony surpassed microsatellite markers, demonstrating the potential of Csd for detailed genetic analyses. Despite challenges in Csd amplification efficiency, a two-step nested PCR protocol proved effective. Interestingly, Csd genotyping alone identified more patrilines than the set of five microsatellite markers, emphasizing its efficacy. Combining Csd and microsatellite genotyping enhances the resolution of genetic studies in honeybee colonies, offering valuable insights into genetic diversity, reproductive success, and social dynamics. The Csd gene emerges as a promising tool for advancing genetic studies in honeybee populations.

Genetic methods in honey bee breeding.

The honey bee Apis mellifera is a rather difficult object for selection due to the peculiarities of its biology. Breeding activities in beekeeping are aimed at obtaining bee colonies with high rates of economically useful traits, such as productivity, resistance to low temperatures and diseases, hygienic behavior, oviposition of the queen, etc. With two apiaries specializing in the breeding of A. m. mellifera and A. m. carnica as examples, the application of genetic methods in the selection of honey bees is considered. The first stage of the work was subspecies identification based on the analysis of the polymorphism of the intergenic mtDNA locus tRNAleu-COII (or COI-COII) and microsatellite nuclear DNA loci Ap243, 4a110, A24, A8, A43, A113, A88, Ap049, A28. This analysis confirmed that the studied colonies correspond to the declared subspecies. In the apiary with A. m. mellifera, hybrid colonies have been identified. A method based on the analysis of polymorphisms of the tRNAleu-COII locus and microsatellite nuclear DNA loci has been developed to identify the dark forest bee A. m. mellifera and does not allow one to differentiate subspecies from C (A. m. carnica and A. m. ligustica) and O (A. m. caucasica) evolutionary lineages from each other. The second stage was the assessment of the allelic diversity of the csd gene. In the apiary containing colonies of A. m. mellifera (N = 15), 20 csd alleles were identified. In the apiary containing colonies of A. m. carnica (N = 44), 41 alleles were identified. Six alleles are shared by both apiaries. DNA diagnostics of bee diseases showed that the studied colonies are healthy. Based on the data obtained, a scheme was developed for obtaining primary material for honey bee breeding, which can subsequently be subjected to selection according to economically useful traits. In addition, the annual assessment of the allelic diversity of the csd gene will shed light on the frequency of formation of new allelic variants and other issues related to the evolution of this gene.

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

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

Revisiting the hymenopteran diploid male vortex: a review of avoidance mechanisms and incidence

AbstractAll hymenopterans (bees, ants, wasps, and sawflies) have haplodiploid sex determination. Generally, this involves haploid males developing from unfertilized eggs and diploid females developing from fertilized eggs, but diploid male production (DMP) frequently occurs. Some species have complementary sex determination (CSD), in which heterozygotes of a csd locus or loci are diploid females, hemizygotes are haploid males, and homozygotes are typically sterile diploid males. These diploid males underlie the diploid male vortex (DMV), the phenomenon of inbreeding reducing csd allelic diversity and causing increasingly higher production of sterile diploid males until the population dies out. The DMV has been cited as an extinction risk for many species across the order and a danger to both human‐controlled and natural hymenopteran populations important in biological control, pollination, and endangered species conservation. However, it has been unclear how frequently it occurs. We review the known mechanisms for DMV avoidance. Many of these mechanisms are linked to lifestyle, and so we structure our investigation around eusocial lifestyle and non‐eusocial (solitary or gregarious) species. We also review documented DMV incidence to make an inference about its prevalence. There are many means to avoid inbreeding or diploid male production, including eusocial‐exclusive polygyny, polyandry, and diploid queen execution. For both eusocial and non‐eusocial species these include biological mechanisms such as multi‐locus CSD or non‐CSD sex determination, diploid male inviability or brood removal, partial diploid male fecundity, and mating avoidance of kin or diploid males; or population genetic mechanisms such as gene flow and dispersal, balancing selection of csd alleles (with some being more prominent in eusocial vs. non‐eusocial species). Documented cases of DMV are uncommon, and incidence is often tied to exacerbatory conditions such as habitat fragmentation and host declines. With this review we suggest that due to numerous avoidance mechanisms, DMV risk may not be as high as previously believed, and requires specific circumstances. In doing so we aim to better inform hymenopteran breeding and conservation efforts.

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.

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

Simple SummaryHoney 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.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.

The effects of functional response and host abundance fluctuations on genetic rescue in parasitoids with single-locus sex determination.

Many parasitoids have single‐locus complementary sex determination (sl‐CSD), which produces sterile or inviable males when homozygous at the sex determining locus. A previous study theoretically showed that small populations have elevated risks of extinction due to the positive feedback between inbreeding and small population size, referred to as the diploid male vortex. A few modeling studies have suggested that the diploid male vortex may not be as common because balancing selection at sex determining loci tends to maintain high allelic diversity in spatially structured populations. However, the generality of the conclusion is yet uncertain, as they were drawn either from models developed for particular systems or from a general‐purpose competition model. To attest the conclusion, we study several well‐studied host–parasitoid models that incorporate functional response specifying the number of attacked hosts given a host density and derive the conditions for a diploid male vortex in a single population. Then, we develop spatially structured individual‐based versions of the models to include female behavior, diploid male fertility, and temporal fluctuations. The results show that producing a handful of successful offspring per female parasitoid could enable parasitoid persistence when a typical number of CSD alleles are present. The effect of functional response depends on the levels of fluctuations in host abundance, and inviable or partially fertile diploid males and a small increase in dispersal can alleviate the risk of a diploid male vortex. Our work supports the generality of effective genetic rescue in spatially connected parasitoid populations with sl‐CSD. However, under more variable climate, the efficacy of the CSD mechanism may substantially decline.

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.

Uneven distribution of complementary sex determiner (csd) alleles in Apis mellifera population

The complementary sex determiner (csd) gene determines the sex of the western honey bee (Apis mellifera L.). Bees that are heterozygous at the csd locus develop into females; whereas hemizygous bees develop into males. The co-occurrence of two identical csd alleles in a single diploid genome leads to the genetic death of the bee. Thus, the maintenance of csd diversity in the population is favoured. The number and distribution of csd alleles is particularly interesting in light of the recent decline in the honey bee population. In this study, we analysed the distribution of csd alleles in two Polish populations separated by about 100 km. We analysed the maternal alleles of 193 colonies and found 121 different alleles. We also analysed the distribution and frequency of the alleles, and found that they are distributed unevenly. We show that the methods that have been used so far to estimate the total worldwide number of csd alleles have significantly underestimated their diversity. We also show that the uneven distribution of csd alleles is caused by a large number of infrequent alleles, which most likely results from the fact that these alleles are generated very frequently.

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.

Everything You Wanted to Know about Honeybee Sex

There is an exquisite genetic control behind a honeybee’s fate in the hive—from the lowly drone to the almighty queen— which literally represents the bee’s knees for evolutionary scientists exploring how multiple mutations, or alleles, of a single gene called the complementary sex determiner (csd )c an have ap rofound inf luence on honeybee society. Unlike people, there are no X and Y sex chromosomes for bees. Rather, sex is determined by csd gene and its allelic composition and whether a queen bee chooses to fertilize her eggs. Female bees (queens or workers) come from fertilized eggs, always receiving two different (heterozygote) copies of csd. Fertile males always come from unfertilized eggs, receiving only one copy of csd. Two identical (homozygote) copies ofcsd in fertilized eggs is always lethal; these individuals are being killed at the early larval stage by worker bees as they would develop into diploid males which do not contribute to colony fitness. Lechner et al. (2013) have now examined the exquisite molecular control behind the sex determination, finely identifying and tracing back a comprehensive number of csd alleles to create a richer understanding of the variability of the csd gene over evolutionary time. They looked at a data set of 244 csd sequences from queens, worker bees, and drones and showed that the total number of csd alleles found in bees ranges from at least 53 (locally) to 87 (worldwide), which is much higher than previously reported. Using an evolutionary model, they also extrapolated the presence of a total 116–145 csd alleles worldwide, a great example of the enormous sequence variability within csd. They were able to finely decipher the minimum number of mutations leading to heterozygous csd, identify faster evolving hot spots within the csd gene, and how these may contribute to variability. “Comprehensive insights into the sequence variability of the sex determining gene csd in honeybees elucidate the evolutionary processes that lead to the enormous number of csd alleles found worldwide,” said Hasselmann. Finally, they traced the data back over evolutionary time and found that a novel csd function affecting sex determination arises about every 400,000 years. The study provides one of the most comprehensive views of the enormous genetic diversity and the evolutionary forces shaping sex determination in bees, as well as how changes in csd affect honeybee colony fitness.

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.

Genetic tests for alleles of complementary-sex-determiner to support honeybee breeding programmes

The honeybee haplodiploid sex determination system depends on genetic variation at the complementary sex-determiner (csd) locus. In closed populations of honeybees, especially those undergoing selective breeding, the number of csd alleles can drop such that brood viability is affected. Here we present two polymerase chain reaction tests that allow the discrimination of csd alleles. Such tests should find utility in bee breeding programmes allowing the tracking and maintenance of csd alleles through successive generations.

The Sex Determination Gene Shows No Founder Effect in the Giant Honey Bee, Apis dorsata

BackgroundAll 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 FindingsWe 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/SignificanceResults 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.

Csd alleles in the red dwarf honey bee (Apis florea, Hymenoptera: Apidae) show exceptionally high nucleotide diversity

Abstract The single locus complementary sex determination (sl-csd) gene is the primary gene determining the gender of honey bees (Apis spp.). While the csd gene has been well studied in the Western honey bee (Apis mellifera), and comparable data exist in both the Eastern honey bee (Apis cerana) and the giant honey bee (Apis dorsata), no studies have been conducted in the red dwarf honey bee, Apis florea. In this study we cloned the genomic region 3 of the A. florea csd gene from 60 workers, and identified 12 csd alleles. Analysis showed that similar to A. mellifera, region 3 of the csd gene contains a RS domain at the N terminal, a proline-rich domain at the C terminal, and a hypervariable region in the middle. However, the A. florea csd gene possessed a much higher level of nucleotide diversity, compared to A. mellifera, A. cerana and Apis dorsata. We also show that similar to the other three Apis species, in A. florea, nonsynonymous mutations in the csd gene are selectively favored in young alleles.

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.