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.

Multisite imaging of neural activity using a genetically encoded calcium sensor in the honey bee.

Understanding of the neural bases for complex behaviors in Hymenoptera insect species has been limited by a lack of tools that allow measuring neuronal activity simultaneously in different brain regions. Here, we developed the first pan-neuronal genetic driver in a Hymenopteran model organism, the honey bee, and expressed the calcium indicator GCaMP6f under the control of the honey bee synapsin promoter. We show that GCaMP6f is widely expressed in the honey bee brain, allowing to record neural activity from multiple brain regions. To assess the power of this tool, we focused on the olfactory system, recording simultaneous responses from the antennal lobe, and from the more poorly investigated lateral horn (LH) and mushroom body (MB) calyces. Neural responses to 16 distinct odorants demonstrate that odorant quality (chemical structure) and quantity are faithfully encoded in the honey bee antennal lobe. In contrast, odor coding in the LH departs from this simple physico-chemical coding, supporting the role of this structure in coding the biological value of odorants. We further demonstrate robust neural responses to several bee pheromone odorants, key drivers of social behavior, in the LH. Combined, these brain recordings represent the first use of a neurogenetic tool for recording large-scale neural activity in a eusocial insect and will be of utility in assessing the neural underpinnings of olfactory and other sensory modalities and of social behaviors and cognitive abilities.

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.

Evidence for Stabilizing Selection Driving Mutational Turnover of Short Motifs in the Eukaryotic Complementary Sex Determiner (Csd) Protein.

Short linear motifs (SLiMs) can play pivotal functional roles in proteins, such as targeting proteins to specific subcellular localizations, modulating the efficiency of translation and tagging proteins for degradation. Until recently we had little knowledge about SLiM evolution. Only a few amino acids in these motifs are functionally important, making them likely to evolve ex nihilo and suggesting that they can play key roles in protein evolution. Several reports now suggest that these motifs can appear and disappear while their function in the protein is preserved, a process sometimes referred to as "turnover". However, there has been a lack of specific experiments to determine whether independently evolved motifs do indeed have the same function, which would conclusively determine whether the process of turnover actually occurs. In this study, we experimentally detected evidence for such a mutational turnover process for nuclear localization signals (NLS) during the post-duplication divergence of the Complementary sex determiner (Csd) and Feminizer (Fem) proteins in the honeybee (Apis mellifera) lineage. Experiments on the nuclear transport activity of protein segments and those of the most recent common ancestor (MRCA) sequences revealed that three new NLS motifs evolved in the Csd protein during the post-duplication divergence while other NLS motifs were lost that existed before duplication. A screen for essential and newly evolved amino acids revealed that new motifs in the Csd protein evolved by one or two missense mutations coding for lysine. Amino acids that were predating the duplication were also essential in the acquisition of the C1 motif suggesting that the ex nihilo origin was constrained by preexisting amino acids in the physical proximity. Our data support a model in which stabilizing selection maintains the constancy of nuclear transport function but allowed mutational turnover of the encoding NLS motifs.

Differential hnRNP D isoform incorporation may confer plasticity to the ESSV-mediated repressive state across HIV-1 exon 3.

Even though splicing repression by hnRNP complexes bound to exonic sequences is well-documented, the responsible effector domains of hnRNP proteins have been described for only a select number of hnRNP constituents. Thus, there is only limited information available for possible varying silencer activities amongst different hnRNP proteins and composition changes within possible hnRNP complex assemblies. In this study, we identified the glycine-rich domain (GRD) of hnRNP proteins as a unifying feature in splice site repression. We also show that all four hnRNP D isoforms can act as genuine splicing repressors when bound to exonic positions. The presence of an extended GRD, however, seemed to potentiate the hnRNP D silencer activity of isoforms p42 and p45. Moreover, we demonstrate that hnRNP D proteins associate with the HIV-1 ESSV silencer complex, probably through direct recognition of "UUAG" sequences overlapping with the previously described "UAGG" motifs bound by hnRNP A1. Consequently, this spatial proximity seems to cause mutual interference between hnRNP A1 and hnRNP D. This interplay between hnRNP A1 and D facilitates a dynamic regulation of the repressive state of HIV-1 exon 3 which manifests as fluctuating relative levels of spliced vpr- and unspliced gag/pol-mRNAs.

Behavioral and molecular studies of quantitative differences in hygienic behavior in honeybees.

BACKGROUND: Hygienic behavior (HB) enables honeybees to tolerate parasites, including infection with the parasitic mite Varroa destructor, and it is a well-known example of a quantitative genetic trait. The understanding of the molecular processes underpinning the quantitative differences in this behavior remains limited. RESULTS: We performed gene expression studies in worker bees that displayed quantitative genetic differences in HB. We established a high and low genetic source of HB performance and studied the engagements into HB of single worker bees under the same environmental conditions. We found that the percentage of worker bees that engaged in a hygienic behavioral task tripled in the high versus low HB sources, thus suggesting that genetic differences may mediate differences in stimulated states to perform HB. We found 501 differently expressed genes (DEGs) in the brains of hygienic and non-hygienic performing workers in the high HB source bees, and 342 DEGs in the brains of hygienic performing worker bees, relative to the gene expression in non-hygienic worker bees from the low HB source group. "Cell surface receptor ligand signal transduction" in the high and "negative regulation of cell communication" in the low HB source were overrepresented molecular processes, suggesting that these molecular processes in the brain may play a role in the regulation of quantitative differences in HB. Moreover, only 21 HB-associated DEGs were common between the high and low HB sources. CONCLUSIONS: The better HB colony performance is primarily achieved by a high number of bees engaging in the hygienic tasks that associate with distinct molecular processes in the brain. We propose that different gene products and pathways may mediate the quantitative genetic differences of HB.

Genomic HEXploring allows landscaping of novel potential splicing regulatory elements.

Effective splice site selection is critically controlled by flanking splicing regulatory elements (SREs) that can enhance or repress splice site use. Although several computational algorithms currently identify a multitude of potential SRE motifs, their predictive power with respect to mutation effects is limited. Following a RESCUE-type approach, we defined a hexamer-based 'HEXplorer score' as average Z-score of all six hexamers overlapping with a given nucleotide in an arbitrary genomic sequence. Plotted along genomic regions, HEXplorer score profiles varied slowly in the vicinity of splice sites. They reflected the respective splice enhancing and silencing properties of splice site neighborhoods beyond the identification of single dedicated SRE motifs. In particular, HEXplorer score differences between mutant and reference sequences faithfully represented exonic mutation effects on splice site usage. Using the HIV-1 pre-mRNA as a model system highly dependent on SREs, we found an excellent correlation in 29 mutations between splicing activity and HEXplorer score. We successfully predicted and confirmed five novel SREs and optimized mutations inactivating a known silencer. The HEXplorer score allowed landscaping of splicing regulatory regions, provided a quantitative measure of mutation effects on splice enhancing and silencing properties and permitted calculation of the mutationally most effective nucleotide.

Distinct subspecies or phenotypic plasticity? Genetic and morphological differentiation of mountain honey bees in East Africa.

Identifying the forces shaping intraspecific phenotypic and genotypic divergence are of key importance in evolutionary biology. Phenotypic divergence may result from local adaptation or, especially in species with strong gene flow, from pronounced phenotypic plasticity. Here, we examine morphological and genetic divergence among populations of the western honey bee Apis mellifera in the topographically heterogeneous East African region. The currently accepted "mountain refugia hypothesis" states that populations living in disjunct montane forests belong to a different lineage than those in savanna habitats surrounding these forests. We obtained microsatellite data, mitochondrial sequences, and morphometric data from worker honey bees collected from feral colonies in three montane forests and corresponding neighboring savanna regions in Kenya. Honey bee colonies from montane forests showed distinct worker morphology compared with colonies in savanna areas. Mitochondrial sequence data did not support the existence of the two currently accepted subspecies. Furthermore, analyses of the microsatellite data with a Bayesian clustering method did not support the existence of two source populations as it would be expected under the mountain refugia scenario. Our findings suggest that phenotypic plasticity rather than distinct ancestry is the leading cause behind the phenotypic divergence observed between montane forest and savanna honey bees. Our study thus corroborates the idea that high gene flow may select for increased plasticity.

Tra2-mediated recognition of HIV-1 5' splice site D3 as a key factor in the processing of vpr mRNA.

Small noncoding HIV-1 leader exon 3 is defined by its splice sites A2 and D3. While 3' splice site (3'ss) A2 needs to be activated for vpr mRNA formation, the location of the vpr start codon within downstream intron 3 requires silencing of splicing at 5'ss D3. Here we show that the inclusion of both HIV-1 exon 3 and vpr mRNA processing is promoted by an exonic splicing enhancer (ESE(vpr)) localized between exonic splicing silencer ESSV and 5'ss D3. The ESE(vpr) sequence was found to be bound by members of the Transformer 2 (Tra2) protein family. Coexpression of these proteins in provirus-transfected cells led to an increase in the levels of exon 3 inclusion, confirming that they act through ESE(vpr). Further analyses revealed that ESE(vpr) supports the binding of U1 snRNA at 5'ss D3, allowing bridging interactions across the upstream exon with 3'ss A2. In line with this, an increase or decrease in the complementarity of 5'ss D3 to the 5' end of U1 snRNA was accompanied by a higher or lower vpr expression level. Activation of 3'ss A2 through the proposed bridging interactions, however, was not dependent on the splicing competence of 5'ss D3 because rendering it splicing defective but still competent for efficient U1 snRNA binding maintained the enhancing function of D3. Therefore, we propose that splicing at 3'ss A2 occurs temporally between the binding of U1 snRNA and splicing at D3.

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

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.

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.

The strength of the HIV-1 3' splice sites affects Rev function.

BACKGROUND: The HIV-1 Rev protein is a key component in the early to late switch in HIV-1 splicing from early intronless (e.g. tat, rev) to late intron-containing Rev-dependent (e.g. gag, vif, env) transcripts. Previous results suggested that cis-acting sequences and inefficient 5' and 3' splice sites are a prerequisite for Rev function. However, we and other groups have shown that two of the HIV-1 5' splice sites, D1 and D4, are efficiently used in vitro and in vivo. Here, we focus on the efficiency of the HIV-1 3' splice sites taking into consideration to what extent their intrinsic efficiencies are modulated by their downstream cis-acting exonic sequences. Furthermore, we delineate their role in RNA stabilization and Rev function. RESULTS: In the presence of an efficient upstream 5' splice site the integrity of the 3' splice site is not essential for Rev function whereas an efficient 3' splice site impairs Rev function. The detrimental effect of a strong 3' splice site on the amount of Rev-dependent intron-containing HIV-1 glycoprotein coding (env) mRNA is not compensatable by weakening the strength of the upstream 5' splice site. Swapping the HIV-1 3' splice sites in an RRE-containing minigene, we found a 3' splice site usage which was variably dependent on the presence of the usual downstream exonic sequence. The most evident activation of 3' splice site usage by its usual downstream exonic sequence was observed for 3' splice site A1 which was turned from an intrinsic very weak 3' splice site into the most active 3' splice site, even abolishing Rev activity. Performing pull-down experiments with nuclear extracts of HeLa cells we identified a novel ASF/SF2-dependent exonic splicing enhancer (ESE) within HIV-1 exon 2 consisting of a heptameric sequence motif occurring twice (M1 and M2) within this short non-coding leader exon. Single point mutation of M1 within an infectious molecular clone is detrimental for HIV-1 exon 2 recognition without affecting Rev-dependent vif expression. CONCLUSION: Under the conditions of our assay, the rate limiting step of retroviral splicing, competing with Rev function, seems to be exclusively determined by the functional strength of the 3' splice site. The bipartite ASF/SF2-dependent ESE within HIV-1 exon 2 supports cross-talk between splice site pairs across exon 2 (exon definition) which is incompatible with processing of the intron-containing vif mRNA. We propose that Rev mediates a switch from exon to intron definition necessary for the expression of all intron-containing mRNAs.

Extended base pair complementarity between U1 snRNA and the 5' splice site does not inhibit splicing in higher eukaryotes, but rather increases 5' splice site recognition.

Spliceosome formation is initiated by the recognition of the 5' splice site through formation of an RNA duplex between the 5' splice site and U1 snRNA. We have previously shown that RNA duplex formation between U1 snRNA and the 5' splice site can protect pre-mRNAs from degradation prior to splicing. This initial RNA duplex must be disrupted to expose the 5' splice site sequence for base pairing with U6 snRNA and to form the active spliceosome. Here, we investigated whether hyperstabilization of the U1 snRNA/5' splice site duplex interferes with splicing efficiency in human cell lines or nuclear extracts. Unlike observations in Saccharomyces cerevisiae, we demonstrate that an extended U1 snRNA/5' splice site interaction does not decrease splicing efficiency, but rather increases 5' splice site recognition and exon inclusion. However, low complementarity of the 5' splice site to U1 snRNA significantly increases exon skipping and RNA degradation. Although the splicing mechanisms are conserved between human and S.cerevisiae, these results demonstrate that distinct differences exist in the activation of the spliceosome.