Active genic machinery for epigenetic RNA modifications in bees.

Epitranscriptomics is an emerging field of investigation dedicated to the study of post-transcriptional RNA modifications. RNA methylations regulate RNA metabolism and processing, including changes in response to environmental cues. Although RNA modifications are conserved from bacteria to eukaryotes, there is little evidence of an epitranscriptomic pathway in insects. Here we identified genes related to RNA m6 A (N6-methyladenine) and m5 C (5-methylcytosine) methylation machinery in seven bee genomes (Apis mellifera, Melipona quadrifasciata, Frieseomelitta varia, Eufriesea mexicana, Bombus terrestris, Megachile rotundata and Dufourea novaeangliae). In A. mellifera, we validated the expression of methyltransferase genes and found that the global levels of m6 A and m5 C measured in the fat body and brain of adult workers differ significantly. Also, m6 A levels were differed significantly mainly between the fourth larval instar of queens and workers. Moreover, we found a conserved m5 C site in the honeybee 28S rRNA. Taken together, we confirm the existence of epitranscriptomic machinery acting in bees and open avenues for future investigations on RNA epigenetics in a wide spectrum of hymenopteran species.

miRNA-34 and miRNA-210 target hexamerin genes enhancing their differential expression during early brain development of honeybee (Apis mellifera) castes.

During the honeybee larval stage, queens develop larger brains than workers, with morphological differentiation appearing at the fourth larval phase (L4), just after a boost in nutritional difference both prospective females experience. The molecular promoters of this caste-specific brain development are already ongoing in previous larval phases. Transcriptomic analyses revealed a set of differentially expressed genes in the L3 brains of queens and workers, which represents the early molecular response to differential feeding females receive during larval development. Three genes of this set, hex70b, hex70c and hex110, are more highly transcribed in the brain of workers than in queens. The microRNAs miR-34, miR-210 and miR-317 are in higher levels in the queens' brain at the same phase of larval development. Here, we tested the hypothesis that the brain of workers expresses higher levels of hexamerins than that of queens during key phases of larval development and that this differential hexamerin genes expression is further enhanced by the repressing activity of miR-34, miR-210 and miR-317. Our transcriptional analyses showed that hex70b, hex70c and hex110 genes are differentially expressed in the brain of L3 and L4 larval phases of honeybee queens and workers. In silico reconstructed miRNA-mRNA interaction networks were validated using luciferase assays, which showed miR-34 and miR-210 negatively regulate hex70b and hex110 genes by directly and redundantly binding their 3'UTR (untranslated region) sequences. Taken together, our results suggest that miR-34 and miR-210 act together promoting differential brain development in honeybee castes by downregulating the expression of the putative antineurogenic hexamerin genes hex70b and hex110.

Caste-specific gene expression underlying the differential adult brain development in the honeybee Apis mellifera.

Apis mellifera adult workers feature more developed key brain regions than queens, which allows them to cope with the broad range of duties they need to perform in a colony. However, at the end of larval development, the brain of queens is largely more developed than that of workers. Major morphogenetic changes take place after metamorphosis that shift caste-specific brain development. Here, we tested the hypothesis that this phenomenon is hormonally governed and involves differential gene expression. Our molecular screening approach revealed a set of differentially expressed genes in Pp (first pharate-adult phase) brains between castes mainly coding for tissue remodelling and energy-converting proteins (e.g. hex 70a and ATPsynß). An in-depth qPCR analysis of the transcriptional behaviour during pupal and pharate-adult developmental stage in both castes and in response to artificially augmented hormone titres of 18 genes/variants revealed that: i. subtle differences in hormone titres between castes might be responsible for the differential expression of the EcR and insulin/insulin-like signalling (IIS) pathway genes; ii. the morphogenetic activity of the IIS in brain development must be mediated by ILP-2, iii. which together with the tum, mnb and caspase system, can constitute the molecular effectors of the caste-specific opposing brain developmental trajectories.

Hormonal control and target genes of ftz-f1 expression in the honeybee Apis mellifera: a positive loop linking juvenile hormone, ftz-f1, and vitellogenin.

Ftz-f1 is an orphan member of the nuclear hormone receptor superfamily. A 20-hydroxyecdysone pulse allows ftz-f1 gene expression, which then regulates the activity of downstream genes involved in major developmental progression events. In honeybees, the expression of genes like vitellogenin (vg), prophenoloxidase and juvenile hormone-esterase during late pharate-adult development is known to be hormonally controlled in both queens and workers by increasing juvenile hormone (JH) titres in the presence of declining levels of ecdysteroids. Since Ftz-f1 is known for mediating intracellular JH signalling, we hypothesized that ftz-f1 could mediate JH action during the pharate-adult development of honeybees, thus controlling the expression of these genes. Here, we show that ftz-f1 has caste-specific transcription profiles during this developmental period, with a peak coinciding with the increase in JH titre, and that its expression is upregulated by JH and downregulated by ecdysteroids. RNAi-mediated knock down of ftz-f1 showed that the expression of genes essential for adult development (e.g. vg and cuticular genes) depends on ftz-f1 expression. Finally, a double-repressor hypothesis-inspired vg gene knock-down experiment suggests the existence of a positive molecular loop between JH, ftz-f1 and vg.

Silencing of Apis mellifera dorsal genes reveals their role in expression of the antimicrobial peptide defensin-1.

Like all other insects, two key signalling pathways [Toll and immune deficiency (Imd)] regulate the induction of honey bee immune effectors that target microbial pathogens. Amongst these effectors are antimicrobial peptides (AMPs) that are presumed to be produced by the nuclear factors kappa B (NF-κB) Dorsal and Relish from the Toll and Imd pathways, respectively. Using in silico analysis, we previously proposed that the honey bee AMP defensin-1 was regulated by the Toll pathway, whereas hymenoptaecin was regulated by Imd and abaecin by both the Toll and Imd pathways. Here we use an RNA interference (RNAi) assay to determine the role of Dorsal in regulating abaecin and defensin-1. Honey bees have two dorsal genes (dorsal-1 and dorsal-2) and two splicing isoforms of dorsal-1 (dorsal-1A and dorsal-1B). Accordingly, we used both single and multiple (double or triple) isoform knockdown strategies to clarify the roles of dorsal proteins and their isoforms. Down-regulation of defensin-1 was observed for dorsal-1A and dorsal-2 knockdowns, but abaecin expression was not affected by dorsal RNAi. We conclude that defensin-1 is regulated by Dorsal (Toll pathway).

MicroRNA signatures characterizing caste-independent ovarian activity in queen and worker honeybees (Apis mellifera L.).

Queen and worker honeybees differ profoundly in reproductive capacity. The queen of this complex society, with 200 highly active ovarioles in each ovary, is the fertile caste, whereas the workers have approximately 20 ovarioles as a result of receiving a different diet during larval development. In a regular queenright colony, the workers have inactive ovaries and do not reproduce. However, if the queen is sensed to be absent, some of the workers activate their ovaries, producing viable haploid eggs that develop into males. Here, a deep-sequenced ovary transcriptome library of reproductive workers was used as supporting data to assess the dynamic expression of the regulatory molecules and microRNAs (miRNAs) of reproductive and nonreproductive honeybee females. In this library, most of the differentially expressed miRNAs are related to ovary physiology or oogenesis. When we quantified the dynamic expression of 19 miRNAs in the active and inactive worker ovaries and compared their expression in the ovaries of virgin and mated queens, we noted that some miRNAs (miR-1, miR-31a, miR-13b, miR-125, let-7 RNA, miR-100, miR-276, miR-12, miR-263a, miR-306, miR-317, miR-92a and miR-9a) could be used to identify reproductive and nonreproductive statuses independent of caste. Furthermore, integrative gene networks suggested that some candidate miRNAs function in the process of ovary activation in worker bees.

Rapid method for DNA extraction from the honey bee Apis mellifera and the parasitic bee mite Varroa destructor using lysis buffer and proteinase K.

We developed a rapid method for extraction of DNA from honey bees, Apis mellifera, and from the parasitic bee mite, Varroa destructor. The advantages include fast processing and low toxicity of the substances that are utilized. We used lysis buffer with nonionic detergents to lyse cell walls and proteinase K to digest proteins. We tested whole thorax, thoracic muscle mass, legs, and antennae from individual bees; the mites were processed whole (1 mite/sample). Each thorax was incubated whole, without cutting, because exocuticle color pigment darkened the extraction solution, interfering with PCR results. The procedure was performed with autoclaved equipment and laboratory gloves. For each sample, we used 100 µL lysis buffer (2 mL stock solution of 0.5 M Tris/HCl, pH 8.5, 10 mL stock solution of 2 M KCl, 500 µL solution of 1 M MgCl2, 2 mL NP40, and 27.6 g sucrose, completed to 200 mL with bidistilled water and autoclaved) and 2 µL proteinase K (10 mg/mL in bidistilled water previously autoclaved, as proteinase K cannot be autoclaved). Tissues were incubated in the solutions for 1-2 h in a water bath (62°-68 °C) or overnight at 37 °C. After incubation, the tissues were removed from the extraction solution (lysis buffer + proteinase K) and the solution heated to 92 °C for 10 min, for proteinase K inactivation. Then, the solution with the extracted DNA was stored in a refrigerator (4°-8 °C) or a freezer (-20 °C). This method does not require centrifugation or phenol/chloroform extraction. The reduced number of steps allowed us to sample many individuals/day. Whole mites and bee antennae were the most rapidly processed. All bee tissues gave the same quality DNA. This method, even using a single bee antenna or a single mite, was adequate for extraction and analysis of bee genomic and mitochondrial DNA and mite genomic DNA.

Organization, evolution and transcriptional profile of hexamerin genes of the parasitic wasp Nasonia vitripennis (Hymenoptera: Pteromalidae).

Hexamerins and prophenoloxidases (PPOs) proteins are members of the arthropod-haemocyanin superfamily. In contrast to haemocyanin and PPO, hexamerins do not bind oxygen, but mainly play a role as storage proteins that supply amino acids for insect metamorphosis. We identified seven genes encoding hexamerins, three encoding PPOs, and one hexamerin pseudogene in the genome of the parasitoid wasp Nasonia vitripennis. A phylogenetic analysis of hexamerins and PPOs from this wasp and related proteins from other insect orders suggests an essentially order-specific radiation of hexamerins. Temporal and spatial transcriptional profiles of N. vitripennis hexamerins suggest that they have physiological functions other than metamorphosis, which are arguably coupled with its lifestyle.

Morphometric and genetic changes in a population of Apis mellifera after 34 years of Africanization.

Though the replacement of European bees by Africanized honey bees in tropical America has attracted considerable attention, little is known about the temporal changes in morphological and genetic characteristics in these bee populations. We examined the changes in the morphometric and genetic profiles of an Africanized honey bee population collected near where the original African swarms escaped, after 34 years of Africanization. Workers from colonies sampled in 1968 and in 2002 were morphometrically analyzed using relative warps analysis and an Automatic Bee Identification System (ABIS). All the colonies had their mitochondrial DNA identified. The subspecies that mixed to form the Africanized honey bees were used as a comparison for the morphometric analysis. The two morphometric approaches showed great similarity of Africanized bees with the African subspecies, Apis mellifera scutellata, corroborating with other markers. We also found the population of 1968 to have the pattern of wing venation to be more similar to A. m. scutellata than the current population. The mitochondrial DNA of European origin, which was very common in the 1968 population, was not found in the current population, indicating selective pressure replacing the European with the African genome in this tropical region. Both morphometric methodologies were very effective in discriminating the A. mellifera groups; the non-linear analysis of ABIS was the most successful in identifying the bees, with more than 94% correct classifications.

Caste development and reproduction: a genome-wide analysis of hallmarks of insect eusociality.

The honey bee queen and worker castes are a model system for developmental plasticity. We used established expressed sequence tag information for a Gene Ontology based annotation of genes that are differentially expressed during caste development. Metabolic regulation emerged as a major theme, with a caste-specific difference in the expression of oxidoreductases vs. hydrolases. Motif searches in upstream regions revealed group-specific motifs, providing an entry point to cis-regulatory network studies on caste genes. For genes putatively involved in reproduction, meiosis-associated factors came out as highly conserved, whereas some determinants of embryonic axes either do not have clear orthologs (bag of marbles, gurken, torso), or appear to be lacking (trunk) in the bee genome. Our results are the outcome of a first genome-based initiative to provide an annotated framework for trends in gene regulation during female caste differentiation (representing developmental plasticity) and reproduction.

A comparative analysis of highly conserved sex-determining genes between Apis mellifera and Drosophila melanogaster

A comparison of the most conserved sex-determining genes between the fruit fly, Drosophila melanogaster, and the honey bee, Apis mellifera, was performed with bioinformatics tools developed for computational molecular biology. An initial set of protein sequences already described in the fruit fly as participants of the sex-determining cascade was retrieved from the Gene Ontology database (http://www.geneontology.org/) and aligned against a database of protein sequences predicted from the honey bee genome. The doublesex (dsx) gene is considered one of the most conserved sex-determining genes among metazoans, and a male-specific partial cDNA of putative A. mellifera dsx gene (Amdsx) was identified experimentally. The theoretical predictions were developed in the context of sequence similarity. Experimental evidence indicates that dsx is present in embryos and larvae, and that it encodes a transcription factor widely conserved in metazoans, containing a DM DNA-binding domain implicated in the regulation of the expression of genes involved in sexual phenotype formation.

Apis mellifera ultraspiracle: cDNA sequence and rapid up-regulation by juvenile hormone.

Two hormones, 20-hydroxyecdysone (20E) and juvenile hormone (JH) are key regulators of insect development including the differentiation of the alternative caste phenotypes of social insects. In addition, JH plays a different role in adult honey bees, acting as a 'behavioural pacemaker'. The functional receptor for 20E is a heterodimer consisting of the ecdysone receptor and ultraspiracle (USP) whereas the identity of the JH receptor remains unknown. We have cloned and sequenced a cDNA encoding Apis mellifera ultraspiracle (AMUSP) and examined its responses to JH. A rapid, but transient up-regulation of the AMUSP messenger is observed in the fat bodies of both queens and workers. AMusp appears to be a single copy gene that produces two transcripts ( approximately 4 and approximately 5 kb) that are differentially expressed in the animal's body. The predicted AMUSP protein shows greater sequence similarity to its orthologues from the vertebrate-crab-tick-locust group than to the dipteran-lepidopteran group. These characteristics and the rapid up-regulation by JH suggest that some of the USP functions in the honey bee may depend on ligand binding.

The vitellogenin of the honey bee, Apis mellifera: structural analysis of the cDNA and expression studies.

The cDNA of Apis mellifera vitellogenin was cloned and sequenced. It is 5440 bp long and contains an ORF of 1770 amino acids (including a putative signal peptide of 16 residues). The deduced amino acid sequence shows significant similarity with other hymenopteran vitellogenins (58% with Pimpla nipponica and 54% with Athalia rosae). The alignment with 19 insect vitellogenins shows a high number of conserved motifs; for example, close to the C-terminus there is a GL/ICG motif followed by nine cysteines, as occurs in all hymenopteran species, and, as in other insect vitellogenins, a DGXR motif is located 18 residues upstream the GL/ICG motif. Phylogenetic analysis of vitellogenin sequences available in insects gave a tree that is congruent with the currently accepted insect phylogenetic schemes. Using two fragments of the vitellogenin cDNA as probes, we analyzed by Northern blot the sex- and caste-specific patterns of vitellogenin expression in pupae and adults of A. mellifera. In queens, vitellogenin mRNA was first detected in mid-late pupal stage, whereas in workers it was first detected in late pupal stage. Vitellogenin mRNA was also observed in drones, although it was first detected not in pupae but in freshly molted adults.

Hemocyte types and total and differential counts in unparasitized and parasitized Anastrepha obliqua (Diptera, Tephritidae) larvae

The hemocyte types, in addition to total and differential hemocyte counts were studied in parasitized and unparasitized Anastrepha obliqua larvae at the beginning and at the end of the third instar. In both developmental phases, in parasitized and unparasitized larvae, prohemocytes, plasmatocytes, granulocytes, adipohemocytes, spherulocytes and oenocytoids cells were observed. Mitotic figures indicate prohemocytes as stem cells. Prohemocytes, plasmatocytes and granulocytes are the most numerous cells in the hemolymph of A. obliqua. Difference in the total number of hemocytes was observed between unparasitized and parasitized larvae at the end of the third instar, but not at the beginning

Ecdysteroid titer and reproduction in queens and workers of the honey bee and of a stingless bee: loss of ecdysteroid function at increasing levels of sociality?

Evidence from field wasps and bumblebees appoints the endocrine system as a mediator between dominance status and ovarian activity in primitively social Hymenoptera. In this comparative study on ecdysteroid titers in the highly social honey bee, Apis mellifera, and a stingless bee, Melipona quadrifasciata, we focussed on the relationship between the ecdysteroid titer, social conditions (presence or absence of the queen), and ovary activity. In contrast to bumblebees, ecdysteroid titers in honey bee and stingless bee workers were either not altered, or dropped to even lower levels after the queen was removed. We also did not detect differences between virgin queens and mated, egg laying queens. These results suggest that ecdysteroids may have lost most of their reproductive functions - yet gained functions in larval caste differentiation - as higher levels of social organization were attained in the evolution of social insects. The observation that ecdysteroid titers are transiently elevated in young workers adds a new, yet functionally still speculative facet to hormonal regulation in insect societies.

Ecdysteroid titers in pupae of highly social bees relate to distinct modes of caste development.

Modifications in endocrine programs are common mechanisms that generate alternative phenotypes. In order to understand how such changes may have evolved, we analyzed the pupal ecdysteroid titers in two closely related, highly social bees: the honey bee, Apis mellifera, and a stingless bee, Melipona quadrifasciata. In both species, the ecdysteroid titers in queens reached their peak levels earlier than in workers. Titer levels at peak maxima did not differ for the honey bee castes, but in Melipona they were twofold higher in queens than in workers. During the second half of pupal development, when the ecdysteroid titers decrease and the cuticle progressively melanizes, the titer in honey bee queens remained higher than in workers, while the reverse situation was observed in Melipona. Application of the juvenile hormone analog Pyriproxyfen((R)) to spinning-stage larvae of Melipona induced queen development. Endocrinologically this was manifest in a queen-like profile of the pupal ecdysteroid titer. Comparing these data with previous results on preimaginal hormone titers in another stingless bee, we conclude that the timing and height of the pupal ecdysteroid peak may depend on the nature of the specific stimuli that initially trigger diverging queen/worker development. In contrast, the interspecific differences in the late pupal ecdysteroid titer profiles mainly seem to be related to caste-specific programs in tissue differentiation, including cuticle pigmentation.

Hemocyte types and total and differential counts in unparasitized and parasitized Anastrepha obliqua (Diptera, Tephritidae) larvae.

The hemocyte types, in addition to total and differential hemocyte counts were studied in parasitized and unparasitized Anastrepha obliqua larvae at the beginning and at the end of the third instar. In both developmental phases, in parasitized and unparasitized larvae, prohemocytes, plasmatocytes, granulocytes, adipohemocytes, spherulocytes and oenocytoids cells were observed. Mitotic figures indicate prohemocytes as stem cells. Prohemocytes, plasmatocytes and granulocytes are the most numerous cells in the hemolymph of A. obliqua. Difference in the total number of hemocytes was observed between unparasitized and parasitized larvae at the end of the third instar, but not at the beginning.

Mispatterning in the ommatidia of Apis mellifera pupae treated with a juvenile hormone analogue.

To further understand the function of morphogenetic hormones in honeybee eye differentiation, the alterations in ommatidial patterning induced by pyriproxyfen, a juvenile hormone (JH) analogue, were studied by scanning and transmission electron microscopy. Prepupae of prospective honeybee workers were treated with pyriproxyfen and the effects on ommatidial differentiation were described at the end of the pupal development. The results show that the entire ommatidia, i.e., the dioptric as well as the receptor systems, were affected by the JH analogue. The wave of ommatidial differentiation, which progresses from the posterior to the anterior region of the pupal eyes, was arrested. In treated pupae, the rhabdomeres only differentiated at the apical axis of the retinula, the secondary and tertiary pigment cells did not develop their cytoplasm protrusions, and the cone cell quartet did not pattern correctly. Simultaneously, an intense vacuolization was observed in cells forming ommatidia. In a previous study we showed that pyriproxyfen exerts an inhibition on pupal ecdysteroid secretion. In this sense, the arrested ommatidial differentiation in pyriproxyfen-treated pupae could be due to a secondary effect resulting from an alteration in pupal ecdysteroid titers.

Hormone-dependent protein patterns in integument and cuticular pigmentation in Apis mellifera during pharate adult development.

The epidermal proteins from staged Apis mellifera pupae and pharate adults and the progress of cuticular pigmentation until adult eclosion were used as parameters to study integument differentiation under hormonal treatment. Groups of bees were treated at the beginning of the pupal stage with the juvenile hormone analog pyriproxyfen (PPN) or as pharate adults with 20-hydroxyecdysone (20E). Another group was treated with both hormones applied successively at these same developmental periods. Controls were maintained without treatment. The epidermal proteins, separated by SDS-PAGE and identified by silver staining, were studied at seven intervals during the pupal and pharate adult stages. The initiation and progress of cuticular pigmentation was also monitored and compared to controls. The results showed that PPN reduced the interval of expression of some epidermal proteins, whereas 20E had an antagonistic effect, promoting a prolongation in the time of expression of the same proteins. In PPN-treated bees, cuticular pigmentation started precociously, whereas in 20E-treated individuals this developmental event was postponed. The double hormonal treatment restored the normal progress of cuticular pigmentation and, to a large extent, the temporal epidermal protein pattern. These results are discussed in relation to the 20E titer modulation and morphogenetic hormone interaction.

Regression of the lateral oviducts during the larval-adult transformation of the reproductive system of Melipona quadrifasciata and Frieseomelitta varia.

Toward the end of the larval phase (prepupa), the reproductive systems of Melipona quadrifasciata and Frieseomelitta varia workers are anatomically similar. Scanning electron microscopy showed that during this developmental phase the right and left ovaries are fused and form a heart-shaped structure located above the midgut. Each ovary is connected to the genital chamber by a long and slender lateral oviduct. During pupal development, the lateral oviducts of workers from both species become extremely reduced due to a drastic process of cell death, as shown by transmission electron microscopy. During the lateral oviduct shortening, their simple columnar epithelial cells show some signs of apoptosis in addition to necrosis. Cell death was characterized by cytoplasmic vesiculation, peculiar accumulation of glycogen, and dilation of cytoplasmic organelles such as mitochondria and rough endoplasmic reticulum. The nuclei, at first irregularly contoured, became swollen, with chromatin flocculation and various areas of condensed chromatin next to the nuclear envelope. At the end of the pupal phase, deep recesses marked the nuclei. At emergence, worker and queen reproductive systems showed marked differences, although reduction in the lateral oviducts was an event occurring in both castes. However, in queens the ovarioles increased in length and the spermatheca was larger than that of workers. At the external anatomical level, the reproductive system of workers and queens could be distinguished in the white- and pink-eyed pupal phase. The metamorphic function of the death of lateral oviduct cells, with consequent oviduct shortening, is discussed in terms of the anatomical reorganization of the reproductive system and of the ventrolateral positioning of adult worker bee ovaries.