Marianne Otte

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.

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.

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.

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.

Tra2-Mediated Recognition of HIV-1 5′ Splice Site D3 as a Key Factor in the Processing of vpr mRNA

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

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.

Honey bee promoter sequences for targeted gene expression

The honey bee, Apis mellifera, displays a rich behavioural repertoire, social organization and caste differentiation, and has an interesting mode of sex determination, but we still know little about its underlying genetic programs. We lack stable transgenic tools in honey bees that would allow genetic control of gene activity in stable transgenic lines. As an initial step towards a transgenic method, we identified promoter sequences in the honey bee that can drive constitutive, tissue‐specific and cold shock‐induced gene expression. We identified the promoter sequences of Am‐actin5c, elp2l, Am‐hsp83 and Am‐hsp70 and showed that, except for the elp2l sequence, the identified sequences were able to drive reporter gene expression in Sf21 cells. We further demonstrated through electroporation experiments that the putative neuron‐specific elp2l promoter sequence can direct gene expression in the honey bee brain. The identification of these promoter sequences is an important initial step in studying the function of genes with transgenic experiments in the honey bee, an organism with a rich set of interesting phenotypes.

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

BackgroundHygienic 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.ResultsWe 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.ConclusionsThe 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.

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.

HEXploring of the HIV-1 genome allows landscaping of new potential splicing regulatory elements.

Effective selection between true and decoy splice sites is critically controlled by flanking splicing regulatory elements (SREs), which can enhance or repress splice site use. Recent experimental evidence suggests that the entire regional context of SREs rather than a single enhancer/silencer hexamer jointly contribute to splicing. Extending the hexamer score concept [Fairbrother et al.: Science 2002, 297:1007-13], we represent the splicing regulatory property of an entire 5’ss neighborhood by a weighted average of normalized Z-scores for all hexamers overlapping with the target region. These novel “HEXplorer” scores describe the degrees of exon- or intron-likeness (ZEI) and enhancer-likeness (ZWS) for a given region upstream of a 5’ss. These scores can be graphically represented by positive or negative oriented areas along the sequence. Mutation effects on an entire 5’ss neighborhood are then captured by comparing the HEXplorer areas of wild type and mutant sequences upstream of a 5’ss. The fundamental datasets of weak and strong 5’ss used in the definition of these HEXplorer scores were derived based on the HBond score that measures the 5’ss complementarity to U1 snRNA. In a first test, we scanned the small non-coding HIV-1 leader exon 3 for regions enriched in SREs. Here, HEXplorer scores correctly indicated both the well-known exonic splicing silencer ESSV and the recently discovered exonic splicing enhancer ESEvpr upstream of 5’ss D3. Next, we tested the HEXplorer’s capability to predict mutations’ potency to modify 5’ss D3 usage. We systematically examined this ESE region using various single and double mutations predicted to either alter 5’ss usage or act neutrally. In 20 tested mutations, the HEXplorer prediction correlated well with the experimentally detected level of exon inclusion. Extending the HEXplorer approach to all HIV-1 exons, we were able to identify three novel exonic splicing enhancers that contribute to the inclusion of the viral exons 2, 2b and exon 4. All three novel ESEs were experimentally confirmed by HEXplorer predicted point-mutations. Beyond application to HIV-1 5’ss usage, the HEXplorer may also prove particularly useful as a method for assessing pathogenic human exonic mutations.