The nuclear and mitochondrial genomes of Frieseomelitta varia - a highly eusocial stingless bee (Meliponini) with a permanently sterile worker caste.

BACKGROUND: Most of our understanding on the social behavior and genomics of bees and other social insects is centered on the Western honey bee, Apis mellifera. The genus Apis, however, is a highly derived branch comprising less than a dozen species, four of which genomically characterized. In contrast, for the equally highly eusocial, yet taxonomically and biologically more diverse Meliponini, a full genome sequence was so far available for a single Melipona species only. We present here the genome sequence of Frieseomelitta varia, a stingless bee that has, as a peculiarity, a completely sterile worker caste. RESULTS: The assembly of 243,974,526 high quality Illumina reads resulted in a predicted assembled genome size of 275 Mb composed of 2173 scaffolds. A BUSCO analysis for the 10,526 predicted genes showed that these represent 96.6% of the expected hymenopteran orthologs. We also predicted 169,371 repetitive genomic components, 2083 putative transposable elements, and 1946 genes for non-coding RNAs, largely long non-coding RNAs. The mitochondrial genome comprises 15,144 bp, encoding 13 proteins, 22 tRNAs and 2 rRNAs. We observed considerable rearrangement in the mitochondrial gene order compared to other bees. For an in-depth analysis of genes related to social biology, we manually checked the annotations for 533 automatically predicted gene models, including 127 genes related to reproductive processes, 104 to development, and 174 immunity-related genes. We also performed specific searches for genes containing transcription factor domains and genes related to neurogenesis and chemosensory communication. CONCLUSIONS: The total genome size for F. varia is similar to the sequenced genomes of other bees. Using specific prediction methods, we identified a large number of repetitive genome components and long non-coding RNAs, which could provide the molecular basis for gene regulatory plasticity, including worker reproduction. The remarkable reshuffling in gene order in the mitochondrial genome suggests that stingless bees may be a hotspot for mtDNA evolution. Hence, while being just the second stingless bee genome sequenced, we expect that subsequent targeting of a selected set of species from this diverse clade of highly eusocial bees will reveal relevant evolutionary signals and trends related to eusociality in these important pollinators.

Molecular determinants of caste differentiation in the highly eusocial honeybee Apis mellifera.

BACKGROUND: In honeybees, differential feeding of female larvae promotes the occurrence of two different phenotypes, a queen and a worker, from identical genotypes, through incremental alterations, which affect general growth, and character state alterations that result in the presence or absence of specific structures. Although previous studies revealed a link between incremental alterations and differential expression of physiometabolic genes, the molecular changes accompanying character state alterations remain unknown. RESULTS: By using cDNA microarray analyses of >6,000 Apis mellifera ESTs, we found 240 differentially expressed genes (DEGs) between developing queens and workers. Many genes recorded as up-regulated in prospective workers appear to be unique to A. mellifera, suggesting that the workers' developmental pathway involves the participation of novel genes. Workers up-regulate more developmental genes than queens, whereas queens up-regulate a greater proportion of physiometabolic genes, including genes coding for metabolic enzymes and genes whose products are known to regulate the rate of mass-transforming processes and the general growth of the organism (e.g., tor). Many DEGs are likely to be involved in processes favoring the development of caste-biased structures, like brain, legs and ovaries, as well as genes that code for cytoskeleton constituents. Treatment of developing worker larvae with juvenile hormone (JH) revealed 52 JH responsive genes, specifically during the critical period of caste development. Using Gibbs sampling and Expectation Maximization algorithms, we discovered eight overrepresented cis-elements from four gene groups. Graph theory and complex networks concepts were adopted to attain powerful graphical representations of the interrelation between cis-elements and genes and objectively quantify the degree of relationship between these entities. CONCLUSION: We suggest that clusters of functionally related DEGs are co-regulated during caste development in honeybees. This network of interactions is activated by nutrition-driven stimuli in early larval stages. Our data are consistent with the hypothesis that JH is a key component of the developmental determination of queen-like characters. Finally, we propose a conceptual model of caste differentiation in A. mellifera based on gene-regulatory networks.

Vitellogenin regulates hormonal dynamics in the worker caste of a eusocial insect.

Functionally sterile honey bee workers synthesize the yolk protein vitellogenin while performing nest tasks. The subsequent shift to foraging is linked to a reduced vitellogenin and an increased juvenile hormone (JH) titer. JH is a principal controller of vitellogenin expression and behavioral development. Yet, we show here that silencing of vitellogenin expression causes a significant increase in JH titer and its putative receptor. Mathematically, the increase corresponds to a dynamic dose-response. This role of vitellogenin in the tuning of the endocrine system is uncommon and may elucidate how an ancestral pathway of fertility regulation has been remodeled into a novel circuit controlling social behavior.

Neuroligin-associated microRNA-932 targets actin and regulates memory in the honeybee.

Increasing evidence suggests small non-coding RNAs (ncRNAs) such as microRNAs (miRNAs) control levels of mRNA expression during experience-related remodelling of the brain. Here we use an associative olfactory learning paradigm in the honeybee Apis mellifera to examine gene expression changes in the brain during memory formation. Brain transcriptome analysis reveals a general downregulation of protein-coding genes, including asparagine synthetase and actin, and upregulation of ncRNAs. miRNA-mRNA network predictions together with PCR validation suggest miRNAs including miR-210 and miR-932 target the downregulated protein-coding genes. Feeding cholesterol-conjugated antisense RNA to bees results in the inhibition of miR-210 and of miR-932. Loss of miR-932 impairs long-term memory formation, but not memory acquisition. Functional analyses show that miR-932 interacts with Act5C, providing evidence for direct regulation of actin expression by an miRNA. An activity-dependent increase in miR-932 expression may therefore control actin-related plasticity mechanisms and affect memory formation in the brain.

Developmental regulation of ecdysone receptor (EcR) and EcR-controlled gene expression during pharate-adult development of honeybees (Apis mellifera).

Major developmental transitions in multicellular organisms are driven by steroid hormones. In insects, these, together with juvenile hormone (JH), control development, metamorphosis, reproduction and aging, and are also suggested to play an important role in caste differentiation of social insects. Here, we aimed to determine how EcR transcription and ecdysteroid titers are related during honeybee postembryonic development and what may actually be the role of EcR in caste development of this social insect. In addition, we expected that knocking-down EcR gene expression would give us information on the participation of the respective protein in regulating downstream targets of EcR. We found that in Apis mellifera females, EcR-A is the predominantly expressed variant in postembryonic development, while EcR-B transcript levels are higher in embryos, indicating an early developmental switch in EcR function. During larval and pupal stages, EcR-B expression levels are very low, while EcR-A transcripts are more variable and abundant in workers compared to queens. Strikingly, these transcript levels are opposite to the ecdysteroid titer profile. 20-hydroxyecdysone (20E) application experiments revealed that low 20E levels induce EcR expression during development, whereas high ecdysteroid titers seem to be repressive. By means of RNAi-mediated knockdown (KD) of both EcR transcript variants we detected the differential expression of 234 poly-A(+) transcripts encoding genes such as CYPs, MRJPs and certain hormone response genes (Kr-h1 and ftz-f1). EcR-KD also promoted the differential expression of 70 miRNAs, including highly conserved ones (e.g., miR-133 and miR-375), as well honeybee-specific ones (e.g., miR-3745 and miR-3761). Our results put in evidence a broad spectrum of EcR-controlled gene expression during postembryonic development of honeybees, revealing new facets of EcR biology in this social insect.

Non-Target Effects of Green Fluorescent Protein (GFP)-Derived Double-Stranded RNA (dsRNA-GFP) Used in Honey Bee RNA Interference (RNAi) Assays.

RNA interference has been frequently applied to modulate gene function in organisms where the production and maintenance of mutants is challenging, as in our model of study, the honey bee, Apis mellifera. A green fluorescent protein (GFP)-derived double-stranded RNA (dsRNA-GFP) is currently commonly used as control in honey bee RNAi experiments, since its gene does not exist in the A. mellifera genome. Although dsRNA-GFP is not expected to trigger RNAi responses in treated bees, undesirable effects on gene expression, pigmentation or developmental timing are often observed. Here, we performed three independent experiments using microarrays to examine the effect of dsRNA-GFP treatment (introduced by feeding) on global gene expression patterns in developing worker bees. Our data revealed that the expression of nearly 1,400 genes was altered in response to dsRNA-GFP, representing around 10% of known honey bee genes. Expression changes appear to be the result of both direct off-target effects and indirect downstream secondary effects; indeed, there were several instances of sequence similarity between putative siRNAs generated from the dsRNA-GFP construct and genes whose expression levels were altered. In general, the affected genes are involved in important developmental and metabolic processes associated with RNA processing and transport, hormone metabolism, immunity, response to external stimulus and to stress. These results suggest that multiple dsRNA controls should be employed in RNAi studies in honey bees. Furthermore, any RNAi studies involving these genes affected by dsRNA-GFP in our studies should use a different dsRNA control.

Functional and evolutionary insights from the genomes of three parasitoid Nasonia species.

We report here genome sequences and comparative analyses of three closely related parasitoid wasps: Nasonia vitripennis, N. giraulti, and N. longicornis. Parasitoids are important regulators of arthropod populations, including major agricultural pests and disease vectors, and Nasonia is an emerging genetic model, particularly for evolutionary and developmental genetics. Key findings include the identification of a functional DNA methylation tool kit; hymenopteran-specific genes including diverse venoms; lateral gene transfers among Pox viruses, Wolbachia, and Nasonia; and the rapid evolution of genes involved in nuclear-mitochondrial interactions that are implicated in speciation. Newly developed genome resources advance Nasonia for genetic research, accelerate mapping and cloning of quantitative trait loci, and will ultimately provide tools and knowledge for further increasing the utility of parasitoids as pest insect-control agents.