The use of Open Reading frame ESTs (ORESTES) for analysis of the honey bee transcriptome

Differential gene expression analysis following olfactory learning in honeybee (Apis mellifera L.)

Insects change their stimulus-response through the perception of associating these stimuli with important survival events such as rewards, threats, and mates. Insects develop strong associations and relate them to their experiences through several behavioral procedures. Among the insects, Apis species, Apis mellifera ligustica are known for their outstanding ability to learn with tremendous economic importance. Apis mellifera ligustica has a strong cognitive ability and promising model species for investigating the neurobiological basis of remarkable olfactory learning abilities. Here we evaluated the olfactory learning ability of A. mellifera by using the proboscis extension reflex (PER) protocol. The brains of the learner and failed-learner bees were examined for comparative transcriptome analysis by RNA-Seq to explain the difference in the learning capacity. In this study, we used an appetitive olfactory learning paradigm in the same age of A. mellifera bees to examine the differential gene expression in the brain of the learner and failed-learner. Bees that respond in 2nd and 3rd trials or only responded to 3rd trials were defined as learned bees, failed-learner individuals were those bees that did not respond in all learning trials The results indicate that the learning ability of learner bees was significantly higher than failed-learner bees for 12 days. We obtained approximately 46.7 and 46.4 million clean reads from the learner bees failed-learner bees, respectively. Gene expression profile between learners' bees and failed-learners bees identified 74 differentially expressed genes, 57 genes up-regulated in the brains of learners and 17 genes were down-regulated in the brains of the bees that fail to learn. The qRT-PCR validated the differently expressed genes. Transcriptome analyses revealed that specific genes in learner and failed-learner bees either down-regulated or up-regulated play a crucial role in brain development and learning behavior. Our finding suggests that down-regulated genes of the brain involved in the integumentary system, storage proteins, brain development, sensory processing, and neurodegenerative disorder may result in reduced olfactory discrimination and olfactory sensitivity in failed-learner bees. This study aims to contribute to a better understanding of the olfactory learning behavior and gene expression information, which opens the door for understanding of the molecular mechanism of olfactory learning behavior in honeybees.

Myosins and DYNLL1/LC8 in the honey bee (Apis mellifera L.) brain

Honey bees have brain structures with specialized and developed systems of communication that account for memory, learning capacity and behavioral organization with a set of genes homologous to vertebrate genes. Many microtubule- and actin-based molecular motors are involved in axonal/dendritic transport. Myosin-Va is present in the honey bee Apis mellifera nervous system of the larvae and adult castes and subcastes. DYNLL1/LC8 and myosin-IIb, -VI and -IXb have also been detected in the adult brain. SNARE proteins, such as CaMKII, clathrin, syntaxin, SNAP25, munc18, synaptophysin and synaptotagmin, are also expressed in the honey bee brain. Honey bee myosin-Va displayed ATP-dependent solubility and was associated with DYNLL1/LC8 and SNARE proteins in the membrane vesicle-enriched fraction. Myosin-Va expression was also decreased after the intracerebral injection of melittin and NMDA. The immunolocalization of myosin-Va and -IV, DYNLL1/LC8, and synaptophysin in mushroom bodies, and optical and antennal lobes was compared with the brain morphology based on Neo-Timm histochemistry and revealed a distinct and punctate distribution. This result suggested that the pattern of localization is associated with neuron function. Therefore, our data indicated that the roles of myosins, DYNLL1/LC8, and SNARE proteins in the nervous and visual systems of honey bees should be further studied under different developmental, caste and behavioral conditions.

Representational Difference Analysis (RDA) reveals differential expression of conserved as well as novel genes during caste-specific development of the honey bee (Apis mellifera L.) ovary

In highly eusocial insects, such as the honey bee, Apis mellifera, the reproductive bias has become embedded in morphological caste differences. These are most expressively denoted in ovary size, with adult queens having large ovaries consisting of 150-200 ovarioles each, while workers typically have only 1-20 ovarioles per ovary. This morphological differentiation is a result of hormonal signals triggered by the diet change in the third larval instar, which eventually generate caste-specific gene expression patterns. To reveal these we produced differential gene expression libraries by Representational Difference Analysis (RDA) for queen and worker ovaries in a developmental stage when cell death is a prominent feature in the ovarioles of workers, whereas all ovarioles are maintained and extend in length in queens. In the queen library, 48% of the gene set represented homologs of known Drosophila genes, whereas in the worker ovary, the largest set (59%) were ESTs evidencing novel genes, not even computationally predicted in the honey bee genome. Differential expression was confirmed by quantitative RT-PCR for a selected gene set, denoting major differences for two queen and two worker library genes. These included two unpredicted genes located in chromosome 11 (Group11.35 and Group11.31, respectively) possibly representing long non-coding RNAs. Being candidates as modulators of ovary development, their expression and functional analysis should be a focal point for future studies.

Suppression subtractive hybridization analysis reveals expression of conserved and novel genes in male accessory glands of the ant Leptothorax gredleri

Background

During mating, insect males eject accessory gland proteins (Acps) into the female genital tract. These substances are known to affect female post-mating behavior and physiology. In addition, they may harm the female, e.g., in reducing its lifespan. This is interpreted as a consequence of sexual antagonistic co-evolution. Whereas sexual conflict abounds in non-social species, the peculiar life history of social insects (ants, bees, wasps) with lifelong pair-bonding and no re-mating aligns the reproductive interests of the sexes. Harming the female during mating would negatively affect male fitness and sexual antagonism is therefore not expected. Indeed, mating appears to increase female longevity in at least one ant species. Acps are presumed to play a role in this phenomenon, but the underlying mechanisms are unknown. In this study, we investigated genes, which are preferentially expressed in male accessory glands of the ant Leptothorax gredleri, to determine which proteins might be transferred in the seminal fluid.

Results

By a suppression subtractive hybridization protocol we obtained 20 unique sequences (USs). Twelve had mutual best matches with genes predicted for Apis mellifera and Nasonia vitripennis. Functional information (Gene Ontology) was available only for seven of these, including intracellular signaling, energy-dependent transport and metabolic enzyme activities. The remaining eight USs did not match sequences from other species. Six genes were further analyzed by quantitative RT-PCR in three life cycle stages of male ants. A gene with carboxy-lyase activity and one of unpredicted function were significantly overexpressed in accessory glands of sexually mature males.

Conclusions

Our study is the first one to investigate differential gene expression in ants in a context related to mating. Our findings indicate that male accessory glands of L. gredleri express a series of genes that are unique to this species, possibly representing novel genes, in addition to conserved ones for which functions can be predicted. Identifying differentially expressed genes might help to better understand molecular mechanisms involved in reproductive processes in eusocial Hymenoptera. While the novel genes could account for rapidly evolving ones driven by intra-sexual conflict between males, conserved genes imply that rather beneficial traits might get fixed by a process described as inter-sexual cooperation between males and females.

Identification of major royal jelly proteins in the brain of the honeybee Apis mellifera

The consumption of royal jelly (RJ) determines the differences between castes and behavioral development in the honeybee Apis mellifera. However, it is not known whether the proteins of RJ are related to these differences, or which proteins are responsible for the changes. To understand the functions of RJ proteins that are present in other tissues of the bee, in addition to hypopharyngeal gland, we used a polyclonal antibody anti-MRJP1 to investigate the presence of this protein in nervous system of honeybee. This study showed the presence of three polypeptides (p57, p70 and p128) in specific tissues of bee brain. Mushroom body, optic lobe and antennal lobe neuropils all contained proteins recognized by anti-MRJP1. Proteomic analysis showed that the three polypeptides are correlated with proteins of the MRJP family. p57 is correlated with MRJP1, p70 with MRJP3, while p128 may be an oligomeric form or a new polypeptide. Immunostaining of the brain and hypopharyngeal gland revealed differential expression of MRJPs in various brain regions and in different honeybee castes and subcastes. The identification and localization of these MRJPs contribute to the elucidation of the biological roles of this protein family.

Transcriptome analysis of Taenia solium cysticerci using Open Reading Frame ESTs (ORESTES)

BACKGROUND

Human infection by the pork tapeworm Taenia solium affects more than 50 million people worldwide, particularly in underdeveloped and developing countries. Cysticercosis which arises from larval encystation can be life threatening and difficult to treat. Here, we investigate for the first time the transcriptome of the clinically relevant cysticerci larval form.

RESULTS

Using Expressed Sequence Tags (ESTs) produced by the ORESTES method, a total of 1,520 high quality ESTs were generated from 20 ORESTES cDNA mini-libraries and its analysis revealed fragments of genes with promising applications including 51 ESTs matching antigens previously described in other species, as well as 113 sequences representing proteins with potential extracellular localization, with obvious applications for immune-diagnosis or vaccine development.

CONCLUSION

The set of sequences described here will contribute to deciphering the expression profile of this important parasite and will be informative for the genome assembly and annotation, as well as for studies of intra- and inter-specific sequence variability. Genes of interest for developing new diagnostic and therapeutic tools are described and discussed.

A non-invasive method for silencing gene transcription in honeybees maintained under natural conditions

In the Apis mellifera post-genomic era, RNAi protocols have been used in functional approaches. However, sample manipulation and invasive methods such as injection of double-stranded RNA (dsRNA) can compromise physiology and survival. To circumvent these problems, we developed a non-invasive method for honeybee gene knockdown, using a well-established vitellogenin RNAi system as a model. Second instar larvae received dsRNA for vitellogenin (dsVg-RNA) in their natural diet. For exogenous control, larvae received dsRNA for GFP (dsGFP-RNA). Untreated larvae formed another control group. Around 60% of the treated larvae naturally developed until adult emergence when 0.5 microg of dsVg-RNA or dsGFP-RNA was offered while no larvae that received 3.0 microg of dsRNA reached pupal stages. Diet dilution did not affect the removal rates. Viability depends not only on the delivered doses but also on the internal conditions of colonies. The weight of treated and untreated groups showed no statistical differences. This showed that RNAi ingestion did not elicit drastic collateral effects. Approximately 90% of vitellogenin transcripts from 7-day-old workers were silenced compared to controls. A large number of samples are handled in a relatively short time and smaller quantities of RNAi molecules are used compared to invasive methods. These advantages culminate in a versatile and a cost-effective approach.

Farnesoic acid O-methyl transferase (FAMeT) isoforms: conserved traits and gene expression patterns related to caste differentiation in the stingless bee, Melipona scutellaris

Farnesoic acid O-methyl transferase (FAMeT) is the enzyme that catalyzes the formation of methyl farnesoate (MF) from farnesoic acid (FA) in the biosynthetic pathway of juvenile hormone (JH). This work reports the cloning, sequencing, and expression of FAMeT gene from the stingless bee Melipona scutellaris (MsFAMeT). The MsFAMeT in silico analysis showed that greatest sequence similarity is found in Apis mellifera and other insects, while relatively less similarity is shown in crustaceans. Evidence of alternative splicing of a 27 nucleotide (nt) microexon explains the presence of the detected isoforms, 1 and 2. The expression analysis of the two isoforms showed a marked difference when castes were compared, suggesting that they could be involved differently in the JH metabolism in M. scutellaris, providing new insights for the comprehension of female plasticity.

Wasp Gene Expression Supports an Evolutionary Link Between Maternal Behavior and Eusociality

The presence of workers that forgo reproduction and care for their siblings is a defining feature of eusociality and a major challenge for evolutionary theory. It has been proposed that worker behavior evolved from maternal care behavior. We explored this idea by studying gene expression in the primitively eusocial wasp Polistes metricus. Because little genomic information existed for this species, we used 454 sequencing to generate 391,157 brain complementary DNA reads, resulting in robust hits to 3017 genes from the honey bee genome, from which we identified and assayed orthologs of 32 honey bee behaviorally related genes. Wasp brain gene expression in workers was more similar to that in foundresses, which show maternal care, than to that in queens and gynes, which do not. Insulin-related genes were among the differentially regulated genes, suggesting that the evolution of eusociality involved major nutritional and reproductive pathways.

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.

A highly divergent gene cluster in honey bees encodes a novel silk family

The pupal cocoon of the domesticated silk moth Bombyx mori is the best known and most extensively studied insect silk. It is not widely known that Apis mellifera larvae also produce silk. We have used a combination of genomic and proteomic techniques to identify four honey bee fiber genes (AmelFibroin1-4) and two silk-associated genes (AmelSA1 and 2). The four fiber genes are small, comprise a single exon each, and are clustered on a short genomic region where the open reading frames are GC-rich amid low GC intergenic regions. The genes encode similar proteins that are highly helical and predicted to form unusually tight coiled coils. Despite the similarity in size, structure, and composition of the encoded proteins, the genes have low primary sequence identity. We propose that the four fiber genes have arisen from gene duplication events but have subsequently diverged significantly. The silk-associated genes encode proteins likely to act as a glue (AmelSA1) and involved in silk processing (AmelSA2). Although the silks of honey bees and silkmoths both originate in larval labial glands, the silk proteins are completely different in their primary, secondary, and tertiary structures as well as the genomic arrangement of the genes encoding them. This implies independent evolutionary origins for these functionally related proteins.

Gene expression profiles underlying alternative caste phenotypes in a highly eusocial bee, Melipona quadrifasciata

To evaluate caste-biased gene expression in Melipona quadrifasciata, a stingless bee, we generated 1278 ESTs using Representational Difference Analysis. Most annotated sequences were similar to honey bee genes of unknown function. Only few queen-biased sequences had their putative function assigned by sequence comparison, contrasting with the worker-biased ESTs. The expression of six annotated genes connected to caste specificity was validated by real time PCR. Interestingly, queens that were developmentally induced by treatment with a juvenile hormone analogue displayed an expression profile clearly different from natural queens for this set of genes. In summary, this study represents an important first step in applying a comparative genomic approach to queen/worker polyphenism in the bee.

Molecular cloning and expression of a hexamerin cDNA from the honey bee, Apis mellifera

A cDNA encoding a hexamerin subunit of the Africanized honey bee (Apis mellifera) was isolated and completely sequenced. In the deduced translation product we identified the N-terminal sequence typical of the honey bee HEX 70b hexamerin. The genomic sequence consists of seven exons flanked by GT/AT exon/intron splicing sites, which encode a 683 amino acid polypeptide with an estimated molecular mass of 79.5 kDa, and pI value of 6.72. Semi-quantitative RT-PCR revealed high levels of Hex 70b message in larval stages, followed by an abrupt decrease during prepupal-pupal transition. This coincides with decaying titers of juvenile hormone (JH) and ecdysteroids that is the signal for the metamorphic molt. To verify whether the high Hex 70b expression is dependent on high hormone levels, we treated 5th instar larvae with JH or 20-hydroxyecdysone (20E). In treated larvae, Hex 70b expression was maintained at high levels for a prolonged period of time than in the respective controls, thus indicating a positive hormone regulation at the transcriptional level. Experiments designed to verify the influence of the diet on Hex 70b expression showed similar transcript amounts in adult workers fed on a protein-enriched diet or fed exclusively on sugar. However, sugar-fed workers responded to the lack of dietary proteins by diminishing significantly the amount of HEX 70b subunits in hemolymph. Apparently, they use HEX 70b to compensate the lack of dietary proteins.