More than royal food - Major royal jelly protein genes in sexuals and workers of the honeybee Apis mellifera
Frontiers in Zoology volume 10, Article number: 72 (2013)
In the honeybee Apis mellifera, female larvae destined to become a queen are fed with royal jelly, a secretion of the hypopharyngeal glands of young nurse bees that rear the brood. The protein moiety of royal jelly comprises mostly major royal jelly proteins (MRJPs) of which the coding genes (mrjp1-9) have been identified on chromosome 11 in the honeybee’s genome.
We determined the expression of mrjp1-9 among the honeybee worker caste (nurses, foragers) and the sexuals (queens (unmated, mated) and drones) in various body parts (head, thorax, abdomen). Specific mrjp expression was not only found in brood rearing nurse bees, but also in foragers and the sexuals.
The expression of mrjp1 to 7 is characteristic for the heads of worker bees, with an elevated expression of mrjp1-4 and 7 in nurse bees compared to foragers. Mrjp5 and 6 were higher in foragers compared to nurses suggesting functions in addition to those of brood food proteins. Furthermore, the expression of mrjp9 was high in the heads, thoraces and abdomen of almost all female bees, suggesting a function irrespective of body section. This completely different expression profile suggests mrjp9 to code for the most ancestral major royal jelly protein of the honeybee.
The cornerstone of eusociality in insects is the caste differentiation . Female individuals are polymorphic with functional and/or morphological different castes resulting in division of labour particular for reproduction. The basics of caste determination are best elucidated for the western honeybee Apis mellifera. Its female caste system with up to 80,000 sterile worker bees and a single queen in the colony is a textbook classic . Female larvae either develop into a queen or a worker bee depending on the quality and the quantity of the diet provided by the nurse bees. Royal larvae exclusively receive royal jelly (RJ), a secretion of the hypopharyngeal glands of nurse bees that rear the brood , whereas pollen and honey is added to the diet of worker destined larvae.
Besides water (60-70%) RJ consists of 10-16% sugar, 12-15% crude protein, 3-6% lipids and traces of salts, free amino acids and vitamins [4, 5]. The protein moiety comprises primarily major royal jelly proteins (MRJPs) of which nine different encoding genes have been identified (mrjp1-9) ([6–9]; for a review see ). Genes encoding MRJPs are not only found within the genus Apis, and are common in other Hymenopteran species including the solitary, parasitoid jewel wasp Nasonia vitripennis, the alfalfa leafcutter bee Megachile rotundata, as well as several bumble bees and ants ( and references therein). Mrjps are also found in the primitively eusocial paper wasp Polistes canadensis , representing the first mrjp- exhibiting species in the family Vespidae. The ancestral state of the genes seems to be the single copy condition [10, 12]. Nevertheless, few cases of copy number radiations have occurred in Hymenopteran species including the parasitic wasp N. vitripennis, the fungus-growing ant Acromyrmex echinator, the argentine ant Linepithema humile and the leaf cutter ant Atta cephalotes . The evolution of mrjps is well understood and they originated by duplication from yellow genes, a group of genes common throughout arthropods [7, 13]. The yellow gene was originally identified as a Drosophila melanogaster mutation interfering with the melanic cuticle pigmentation . Subsequently yellow genes were shown to be associated with sex-specific reproductive maturation and behaviour [13, 15].
All MRJPs, except for MRJP8, have been detected in RJ or the hypopharyngeal glands of A. mellifera [9, 16–25]. The genes mrjp1 to 5 have been suggested to be primarily expressed in the hypopharyngeal glands of nurses but not in those of foragers [7, 26–31]. This matched well with the notion that MRJPs are primarily a nutrient for the developing larvae [9, 32, 33]. In addition, MRJPs have been repeatedly studied because of their potential role as a queen determinator in RJ. However, many of these studies are inconclusive and not in line with each other. MRJP1  and MRJP3  have been claimed to induce larvae growth and function as queen determinators. However, also changing glucose and fructose concentrations [36, 37] or the unsaturated fatty acid 10-hydroxy-2-decenoic acid  interfered with caste determination. In general it seems clear that the ancestral MRJP did not evolve to a caste determination protein in general in social insects. For example the genome of the solitary wasp Nasonia vitripennis which has no workers comprises a suite of 10 mrjp genes. Furthermore, given that mrjp genes are also present in e.g., the red harvester ant Pogonomyrmex barbatus which shows genetic caste determination based on an individual's genotype , it seems to be unlikely that the ancient MRJP functions as ‘queen determinator’. So if any MRJP is involved in environmental queen determination, this must be a very specific function that exclusively evolved in Apis.
The expression of mrjp1-8 has been repeatedly shown also in the brain of nurse bees [29, 40–44], and mrjp1 and 3 were shown to be expressed in drones (head, body, larvae) and queens (ovary, larvae) . MRJP1 to 3 were also found in the haemolymph of larvae  with a significant higher amount than in pupae haemolymph  and they were shown to be down-regulated or depleted after infection with a severe pathogen of the honeybee, Paenibacillus larvae. Hence, MRJPs seem to have important functions for honeybee physiology, development and colonial organization in general [10, 44] and not just as a food compound in RJ.
In this study we provide a comprehensive analysis on major royal jelly protein gene expression in all honeybee castes and both sexes (workers (nurses, foragers), queens (unmated, mated) and drones) in various body parts (head, thorax, abdomen) to assess to what extent expression of these genes is affected by caste and sex.
Confirmation of primer specificity and validation of reference genes
The specificity of qPCR products was documented with the high-resolution automatic capillary electrophoresis system QIAxcel® (Qiagen, Hilden, Germany). No primer-dimers or unspecific products were generated and all mrjp transcripts resulted in single products with the predicted size (Additional file 1: Figure S1). In addition, melting curve analyses resulted in the corresponding specific melting temperatures (Additional file 1: Table S1).
Two genes, RpS5a and arp1, were analyzed in regard to their suitability to serve as reference genes. According to Pfaffl et al. , any studied gene with a Ct value standard deviation higher than 1 can be considered as inconsistent. Whereas RpS5a showed a standard deviation of 0.86, arp1 had with 1.55 a much higher value. Thus, arp1 was coregulated within the different castes, sexes and body sections and excluded as reference gene from further analyses.
Tissue specific gene expression
Mrjp transcript abundance was highly variable among the various groups and tissues (Figure 1, Table 1). Within the drones and queens mrjp transcript abundance was generally low in all tissues, particulary for mrjp1 to 5 (Figure 1a, Table 1). Mrjps1 to 7 showed significantly increased transcript abundances in heads of nurses (N) and foragers (F), compared to caged workers (C), drones (D) and queens (Q) (Figure 2 and Additional file 1: Table S2). This drastic different gene expression leads to the division of the dendrogram into two different clusters, one comprising foragers and nurses and the other one the sexuals (Figure 1). Mrjp8 was very evenly expressed in almost all analyzed groups and body sections (mean ± SE: 0.25 ± 0.03) suggesting a function independent of caste, sex and tissue. Mrjp9 showed a completely different expression pattern from all of the other mrjps, with a significant increase in the expression in all worker bees and virgin queens compared to drones and mated queens (Figure 3), leading to a separate cluster in the dendrogram (Figure 1a). Except for drones, mrjp9 was strongly expressed in the thorax and abdomen of all groups and always among the genes with the highest transcript abundances (Figure 1a and Additional file 1: Table S3).
Expression in workers
To better visualize the different regulation of mrjps especially among the worker caste, Figure 1b shows the expression data normalized to caged workers, which had a protein poor diet and could not activate food glands. Those workers showed only a basic expression of mrjp1 to 7 in the head compared to nurse bees in the hive (Figure 1b and Additional file 1: Table S2). The observed up-regulation in nurses was fair for mrjp 5 and 6 (30- and 6-fold, respectively; P < 0.001, one-way ANOVA, Bonferroni post-hoc test), higher for mrjp4 and 7 (both 600-fold, P < 0.001) and tremendous for mrjp3 (3,000-fold, P < 0.001), mrjp2 (15,000-fold, P < 0.001) and mrjp1 (13,000-fold, P < 0.001) (Figure 2).
Surprisingly, this general increase in the expression of the seven mrjp genes was not only found in nurses but also in foragers (Figure 2 and Additional file 1: Table S2). Mrjps 1, 2, 5 and 7 were significantly up-regulated in foragers compared to caged workers (P < 0.001) but showed no significant difference to the up-regulation in nurses. The expression level of mrjp6 was even 10-fold higher in foragers compared to nurses (P < 0.001). The increase in mrjp3 expression in foragers (25-fold, P < 0.001) was significantly lower compared to the 3,000-fold increase in nurses (P < 0.001). This reduced increase was also found for mrjp4, which was only 28-fold higher in foragers, but 600-fold higher in nurses compared to caged workers (P < 0.001). In all three types of worker bees, mrjp1 to 7 were only slightly expressed in thorax and abdomen, with an apparent minimal expression of mrjp1 to 4 (relative gene expression below 0.01, Figure 1a, Table 1). The expression of mrjp8, similar in head and abdomen among the three worker types, differed only in the thorax with a 4- to 5-fold increased expression level in foragers and nurses compared to caged workers (P < 0.01) (Figure 1b). Mrjp9 showed no difference between caged workers and nurses but the expression levels in forager head and thorax were 2- to 6-fold higher (P < 0.001).
Expression in sexuals - queens and drones
There were no significant differences in transcript abundance for any of the mrjp genes between the two types of mated queens (Figure 1 and Additional file 1: Table S2). Virgin queens differed from mated queens primarily in mrjp9, with a 5- to 30-fold increased expression in all body parts (P < 0.001, one-way ANOVA, Bonferroni post-hoc test, pooled mated queens) (Figure 3). Mrjp1 to 5 were only marginally expressed in all body sections of the three queen types with a relative expression value ≤ 0.01 (Figure 1a, Table 1).
The drones generally had the lowest mrjp expression levels. The highest expressed mrjp is with a relative expression of 0.29 ± 0.05 (mean ± SE) mrjp8, the one that is barely regulated and consistently expressed over all groups (Figure 1, Table 1). All of the other mrjps are even lower expressed with a relative expression ≤ 0.1. The analysis of mrjp expression profiles shows that drones cluster together with the three types of queens in one cluster far from nurses and foragers (Figure 1).
Mrjps are expressed in all body sections of both sexes and all castes of the honeybee Apis mellifera suggesting that the functions of these proteins are wide-ranging and not only constrained to the food glands. Our results show that the focus of the expression of mrjp1 to 7 lies clearly in the worker heads (Figure 1). This is in accordance with previous studies confirming the association of some MRJPs with the activated food glands (for a recent review see ). Since the food glands are absent in drones and queens and not developed in the caged workers it is not surprising to see mrjp 1 to 7 only minimal expressed (Figure 1). In the thorax and abdomen, the expression of the genes is similar low but not absent, showing that mrjps are involved in more than just producing nutritional protein.
In the recently compiled honeybee protein atlas of Chan et al.  MRJPs were shown to occur in various organs of the abdomen but not in the thoracic muscle and the thoracic salivary gland (see also ). Nevertheless, various MRJPs, except for MRJP4, 6 and 8, were identified in the honeybee’s nerve chord  that is crossing the thorax and abdomen . Therefore, the basic expression of some mrjps in the thorax might be caused by expression in the nerve chord ganglia.
The ancestors – MRJP8 and 9
The two genes mrjp 8 and 9 showed a very different expression profile compared to the other seven mrjps. Mrjp8 and 9 are not up-regulated in worker heads and expressed fair to high in the thorax and abdomen (Figure 1). Thus, both proteins are not tissue specific. Already Peiren et al. [50, 51] showed that both proteins are not constrained to the food glands but also present in bee venom. In fact, MRJP8 has never been identified as component of RJ in A. mellifera and also MRJP9 was only identified in three out of numerous studies on the RJ proteome [20, 24, 25] alluding to a low amount of the protein in RJ. This is well in line with the previous notion that mrjp1-7 derived from mrjp8 and 9 during honeybee evolution . Whereas mrjp8 was consistently expressed in almost all analyzed samples, mrjp9 was up-regulated in workers and virgin queens, compared to drones and mated queens (Figure 3, Table 1). Furthermore, mrjp9 was the only mrjp with an expression higher than the reference gene in the thorax and the abdomen. As mrjp9 shows a different expression profile, it clusters also in the dendrogram independent from all of the other mrjps (Figure 1). In addition, the MRJPs of species that possess just one MRJP (e.g., Bombus terrestris, Camponotus floridanus, Harpegnathos saltator and Megachile rotundata) show always the highest similarity (44 to 56%) to MRJP9 of the honeybee (Additional file 1: Table S4). This defines mrjp9 as the most ancestral mrjp in A. mellifera, suggesting a completely different function for MRJP9 from the other MRJPs.
Comparison across caste – nurses and foragers
There was a tremendous up-regulation of mrjp1-7 in both the heads of nurses and foragers compared to the caged controls. The transcript abundances of mrjp1-4 and 7 in nurses were higher than in foragers, although not always significantly higher (Figure 2). The lack of significance may have been due to the low sample size but also to the lack of age control in the foragers. Indeed, the variance for transcript abundance was significantly higher among the foragers than the nurse bees. We cannot exclude, that some foragers still had partially activated food glands which might have diluted the difference to the nurses. Although workers are typically considered as mature foragers if they carry pollen, there will be age variance among the analyzed foragers with potential residual food gland activity. Some workers may start foraging as early as day 10 (~10% of the time) . Thus, foragers arriving with pollen at the hive may still be also performing nursing tasks within the hive. Since the caged bees were kept on a protein free diet and had no contact to brood, their hypopharyngeal glands were not activated [53, 54] and did produce no or heavily reduced MRJPs 1–7. Clearly these bees showed that the hive context and the diet have a massive impact on mrjp transcripts and caged workers are important controls to visualize this effect.
In RJ, MRJP1, 2, 3 and 5 amount to 82% of total RJ proteins (31%, 16%, 26% and 9%, respectively)  and those four MRJPs are commonly identified with 1D-SDS-PAGE. This suggests that primarily MRJP1, 2, 3 and 5 have a function in RJ.
In this study, mrjp3 was most strongly up-regulated in nurse bees (3,000-fold compared to caged workers and 130-fold compared to foragers). As the expression in all of the other analysed groups was only marginal, the protein seems to be particularly specific for young nurses and serves as important RJ ingredient. The nutritional function as food protein was already suggested for MRJP3 by Schmitzová et al. , especially because the protein contains a repetitive pentapeptide motive comprising several nitrogen rich amino acids . MRJP3 was also associated with queen differentiation as specific post-translational modified isoforms of MRJP3 secreted into RJ lead to the development of queens .
Since also MRJP5 contains a repetitive motive with nitrogen rich amino acids, it has been considered as nutritional protein as well [9, 32] and we found a 30-fold upregulation of mrjp5 in nurse heads compared to caged workers. However, forager heads revealed a further 2-fold upregulation. Although this was not significant compared to the nurse bees, it nevertheless suggests that MRJP5 may have yet unknown functions in foragers beyond a plain food protein. Indeed, MRJP5 was detected in a significant higher amount in the brain proteome of foragers compared to nurses, whereas MRJP1-4 and 7 were higher in nurse brains  similar to our study.
Mrjp6 was the only mrjp transcript which was significantly more abundant in foragers than in nurses (10-fold), suggesting again a function more specific to older worker bees. Both MRJP5 and 6 show the highest sequence identity (74%) among all MRJPs . Since both genes also cluster together in the dendrogram based on mrjp expression levels (Figure 1a), it seems highly suggestive that both proteins are functionally important for foraging bees.
Mrjp7 was 600-fold up-regulated in nurse heads compared to caged workers but is not a major component of RJ. Hojo et al.  found mrjp7 to be one of the most abundant transcripts of the mushroom bodies in the honeybee brain. Furthermore, it was only marginally expressed in the hypopharyngeal glands . These results point to a function of MRJP7 primarily in the brain.
Hojo et al.  also found mrjp2 to be expressed in the mushroom bodies. In this study the expression level of mrjp2 was very similar to the expression level of mrjp3 in nurse bee heads. Since the amounts of these both proteins are very different in RJ (MRJP2 – 16%, MRJP3 – 26%, ) this supports the presence of an additional function of MRJP2 most likely in the brain.
MRJP1 is the most abundant protein of RJ  and shows with a relative expression of 226 in nurse bee heads by far the highest expression of all mrjps in any analyzed group. Since mrjp1 expression is still elevated in foragers it is unlikely that the nutritional role and/or the possible role in queen differentiation are the only functions of MRJP1. The expression of mrjp1 was already shown also in the brain [29, 40–44] and more specifically in the Kenyon cells of the brain’s mushroom bodies, cells involved in the formation and storage of memory and learning .
Comparison across sexuals – drones and queens
The expression patterns of mrjps are very similar among all sexuals and across all body sections (Figure 1). Most striking is that none of the mrjps showed any up-regulation in the heads. So whatever brain function MRJPs have in workers, the sexuals must do without it. Among queens, expression of mrjp1 to 8 seems to be independent of queen age and mating. Virgin queens show a 5- to 30-fold up-regulation of mrjp9 in all body parts. Nevertheless, the similarity of the expression levels among the queens is high and they form a tight cluster in the dendrogram, supported by high bootstrap values (Figure 1). Drones show barely any differences to mated queens but have a significant down-regulation of mrjp9 compared to virgin queens. But as long as the functions of MRJP9 are not known, any functional interpretation remains pure speculation. Despite this small difference in mrjp9 expression, drones show with a similarity distance of 3.76 a higher similarity to queens than to any analysed group of the worker caste (similarity distances above 4.4).
To date, mrjp genes have been identified with different copy numbers in several species of the Hymenoptera, including the solitary parasitoid wasp Nasonia vitripennis, the leaf cutter ant Atta cephalotes and the buff-tailed bumblebee Bombus terrestris (for review see ). The high copy numbers in the genus Apis resulted from duplications of an ancient mrjp gene [7, 10] and the results of the present study further confirm that this ancestor is mrjp9.
Although the influence of caste is much higher on mrjp expression than the influence of sex, we can show that mrjps are widely expressed in all castes and both sexes and not limited to the food glands of nurse bees. Queens and drones show almost the same expression patterns for all mrjps. Thus, both sexes cluster together in a discrete branch in the group-tree based on mrjp expression (Figure 1). Furthermore, the sexuals are more similar to the caged controls (similarity distance 4.44) than to the other two worker castes (similarity distance 6.66), due to the head expression of mrjp1-7 in foragers and hive nurses. So clearly MRJP1-7 are typical worker caste proteins connected to labour in the colony. But even in the worker caste, functions of the proteins differ and some MRJPs, like MRJP1 and 2, are obviously polyfunctional.
In spite of this dominating and clear cut differentiation in the female caste, all of our results concerning mrjp expression in workers suggest that MRJPs not only have a nutritional function in royal jelly, but also a more general physiological one for all organisms in the colony. In particular, the evolutionary old mrjp9 and mrjp8 cannot serve as food proteins since neither queens nor drones feed anybody in the colony. Furthermore, the consistent expression across all body sections indicates a more general physiological role. For the sole MRJP in B. terrestris a non-nutritive function was already suggested before . The picture of mrjp s in honeybees remains complex and albeit feeding larvae represents without doubt a function for some MRJPs others have clearly profound roles in the brain and fulfill tasks that still have to be elucidated.
Material and methods
Honeybees (Apis mellifera) were sampled in June and July 2013 from the University apiary. Drones (D) and pollen foragers (F) were caught directly at the flight entrance of one colony and frozen in liquid nitrogen.
To rear nurse bees, a brood frame was removed from the hive and incubated at 34°C and ~60% relative humidity until the bees hatched. Freshly emerged workers were paint marked and returned to the hive (nurse bees (N)). To test for the effect of age, we also reared freshly emerged workers in hoarding cages in the incubator fed with a non-protein diet (79% powdered sugar (w/w), 20% honey (w/w)) ad libitum to obtain equally aged workers but without developed food glands (caged workers (C)) . After four days, caged workers and hive nurses were freeze killed in liquid nitrogen.
Queen bees were raised from queen cells with the help of nurse bees in incubators at 34°C and ~60% humidity. Virgin queens (QV) were freeze killed a few days after hatching. For each group (drones, caged workers, nurses, foragers and virgin queens) nine bees were collected and stored at -80°C until further processing. Queens destined to be mated were introduced into mating hives with ~2000 worker bees and allowed to perform mating flights. Three of the mated queens were freeze killed directly after initiation of oviposition and the production of first eggs (QE) whereas three other queens (QP) were freeze killed until the first pupae had developed (~ 15 d after first oviposition - confirmed by white worker pupae with red eyes).
Total RNA was extracted from the head, thorax and abdomen of nine individuals of each group (three for the two mated queen groups) using the RNeasy Mini Kit (Qiagen, Hilden, Germany) according to the manufacture’s protocol. Quality and quantity of the total RNA were photometrically determined with a NanoDrop 1000 (Thermo Fisher Scientific, Wilmington, DE, USA). 500 ng total RNA were reverse transcribed using 0.4 μg Oligo (dT)15 Primer (Promega, Mannheim, Germany), 0.8 μl dNTPs (10 mM) and 80 U M-MLV reverse transcriptase (Promega, Mannheim, Germany). cDNA was purified with the QIAquick PCR Purification Kit (Qiagen, Hilden, Germany) as described in the manufacture’s protocol and the concentration set to 15 ng/μl. cDNA of mated queens was used directly for quantitative real-time PCR (qPCR) analyses whereas the cDNA of all other groups was pooled to minimize individual variation from three individuals, i.e. body sections, to one pool. Finally three pools per body section and group, except for mated queens, were used.
For qPCR reactions, 1 μl cDNA (15 ng/μl) was mixed with 5 μl SensiMixPlus SYBR & Fluorescein Kit (Bioline, Luckenwalde, Germany), 0.3 μM of each gene specific primer and 3.4 μl DEPC-water. Gene specific primers were designed to span at least one intron using sequences of Apis mellifera mrjp1-9 that were available on GenBank in May 2013 and the programme Primer-BLAST of the National Center for Biotechnology Information (NCBI) (Additional file 1: Table S1). Ribosomal protein S5a (RpS5a) and actin related protein 1 (arp1, also known as Actin) were initially chosen in order to standardize expression levels between individuals and groups [56, 57].
The same qPCR protocol was used for all primer pairs. An initial denaturation step of 10 min at 95°C was followed by 40 amplification cycles (95°C, 15 sec; 57°C, 30 sec; 72°C, 30 sec), and a subsequent melting curve analysis between 55°C and 98°C, reading the fluorescence at 1°C increments. Two technical replicates were run for each sample using Chromo4™ (Bio-Rad, Munich, Germany) and repeated if necessary, to reach at maximum an in between replicate threshold cycle (Ct) difference of 0.5.
LinRegPCR version 12.10  was used to determine the Ct values after baseline subtraction. PCR efficiency for each target gene was estimated by serial dilution qPCR (Additional file 1: Table S1) and relative target gene expression was determined as described  using RpS5a  as honeybee reference gene (Table 1). Heat maps combined with dendrograms were used to show gene and caste specific expression differences. We performed a hierarchical clustering analysis using the software MultiExperiment Viewer (MeV) version 4.9. Within this analysis, a dendrogram tree was built for each group (samples, genes) based on the raw gene expression data. Furthermore, we used the option ‘optimized gene and caste leaf order’ pending on the data set. To build the heat map and dendrogram the Euclidean distance between samples and genes and single linkage clustering as clustering method were used to build the final figure . The Clustering Calculator software (http://www2.biology.ualberta.ca/jbrzusto/cluster.php) was used with the same settings as mentioned above to verify the topology of the inferred dendrogram (bootstrap resampling, 1000 replications), and to calculate similarity distances.
All statistical analyses were performed with STATISTICA 8.0 (StatSoft, Tulsa, OK, USA). Data were tested for deviations form a normal distribution by Kolmogorov-Smirnov tests and log- or Box-Cox-transformed in the case of significant deviations from normality and homoscedasticity. Mrjp gene specific comparisons between and within groups (e.g. caste and/or body section) were done using one-way analysis of variance (ANOVA) with Bonferroni post-hoc test. Statistical details for all groups and genes are given in the Additional file 2.
Major royal jelly protein
Queens with eggs
Queens with pupae
Wilson EO: The insect societies. 1971, Cambridge: Belknap Press of Harvard University Press
Winston ML: The biology of the honey bee. 1987, Cambridge: Harvard University Press
Snodgrass RE: Anatomy of the Honey Bee. 1910, Washington: Government Printing Office
Von Planta A: Ueber den Futtersaft der Bienen. Hoppe-Seyl Zeitschr Physiol Chem. 1888, 12: 327-354.
Rembold H: Royal Jelly. Queen Rearing: Biological Basis and Technical Instruction, Chapter 2. Edited by: Ruttner F. 1983, Bucharest: Apimondia Publishing House, 35-42.
Albert Š, Klaudiny J: The MRJP/YELLOW protein family of Apis mellifera: identification of new members in the EST library. J Insect Physiol. 2004, 50: 51-59. 10.1016/j.jinsphys.2003.09.008.
Drapeau MD, Albert Š, Kucharski R, Prusko C, Maleszka R: Evolution of the Yellow/Major Royal Jelly Protein family and the emergence of social behavior in honey bees. Genome Res. 2006, 16: 1385-1394. 10.1101/gr.5012006.
Honeybee Genome Sequencing Consortium: Insights into social insects from the genome of the honeybee Apis mellifera. Nature. 2006, 443: 931-949. 10.1038/nature05260.
Schmitzová J, Klaudiny J, Albert Š, Schröder W, Schreckengost W, Hanes J, Júdová J, Simúth J: A family of major royal jelly proteins of the honeybee Apis mellifera L. Cell Mol Life Sci. 1998, 54: 1020-1030. 10.1007/s000180050229.
Buttstedt A, Moritz RFA, Erler S: Origin and function of the major royal jelly proteins of the honeybee (Apis mellifera) as members of the yellow gene family. Biol Rev. 2013, (doi: 10.1111/brv.12052) [epub ahead of print]
Sumner S, Pereboom JJ, Jordan WC: Differential gene expression and phenotypic plasticity in behavioural castes of the primitively eusocial wasp, Polistes canadensis. Proc Biol Sci. 2006, 273: 19-26. 10.1098/rspb.2005.3291.
Kupke J, Spaethe J, Mueller MJ, Rössler W, Albert Š: Molecular and biochemical characterization of the major royal jelly protein in bumblebees suggest a non-nutritive function. Insect Biochem Mol Biol. 2012, 42: 647-54. 10.1016/j.ibmb.2012.05.003.
Ferguson LC, Green J, Surridge A, Jiggins CD: Evolution of the insect yellow gene family. Mol Biol Evol. 2011, 28: 257-272. 10.1093/molbev/msq192.
Nash WG: Pattern of pigmentation color states regulated by the y locus in Drosophila melanogaster. Develop Biol. 1976, 48: 336-343. 10.1016/0012-1606(76)90095-6.
Drapeau MD: A novel hypothesis on the biochemical role of the Drosophila Yellow protein. Biochem Biophys Res Commun. 2003, 311: 1-3. 10.1016/j.bbrc.2003.09.106.
Sano O, Kunikata T, Kohno K, Iwaki K, Ikeda M, Kurimoto M: Characterization of royal jelly proteins in both Africanized and European honeybees (Apis mellifera) by two-dimensional gel electrophoresis. J Agric Food Chem. 2004, 52: 15-20. 10.1021/jf030340e.
Scarselli R, Donadio E, Giuffrida MG, Fortunato D, Conti A, Balestreri E, Felicioli R, Pinzauti M, Sabatini AG, Felicioli A: Towards royal jelly proteome. Proteomics. 2005, 5: 769-776. 10.1002/pmic.200401149.
Li J, Wang T, Zhang Z, Pan Y: Proteomic analysis of royal jelly from three strains of western honeybees (Apis mellifera). J Agric Food Chem. 2007, 55: 8411-8422. 10.1021/jf0717440.
Li JK, Feng M, Zhang L, Zhang ZH, Pan YH: Proteomics analysis of major royal jelly protein changes under different storage conditions. J Proteome Res. 2008, 7: 3339-3353. 10.1021/pr8002276.
Schönleben S, Sickmann A, Mueller MJ, Reinders J: Proteome analysis of Apis mellifera royal jelly. Anal Bioanal Chem. 2007, 389: 1087-1093. 10.1007/s00216-007-1498-2.
Furusawa T, Rakwal R, Nam HW, Shibato J, Agrawal GK, Kim YS, Ogawa Y, Yoshida Y, Kouzuma Y, Masuo Y, Yonekura M: Comprehensive royal jelly (RJ) proteomics using one- and two-dimensional proteomics platforms reveals novel RJ proteins and potential phospho/glycoproteins. J Proteome Res. 2008, 7: 3194-3229. 10.1021/pr800061j.
Yu F, Mao F, Jianke L: Royal jelly proteome comparison between A. mellifera ligustica and A. cerana cerana. J Proteome Res. 2010, 9: 2207-2215. 10.1021/pr900979h.
Han B, Li C, Zhang L, Fang Y, Feng M, Li J: Novel royal jelly proteins identified by gel-based and gel-free proteomics. J Agric Food Chem. 2011, 59: 10346-10355. 10.1021/jf202355n.
Zhang L, Fang Y, Li R, Feng M, Han B, Zhou T, Li J: Towards posttranslational modification proteome of royal jelly. J Proteomics. 2012, 75: 5327-5341. 10.1016/j.jprot.2012.06.008.
Fujita T, Kozuka-Hata H, Ao-Kondo H, Kunieda T, Oyama M, Kubo T: Proteomic analysis of the royal jelly and characterization of the functions of its derivation glands in the honeybee. J Proteome Res. 2013, 12: 404-11. 10.1021/pr300700e.
Klaudiny J, Kulifajová J, Crailsheim K, Simúth J: New approach to the study of division of labour in the honeybee colony (Apis mellifera L.). Apidologie. 1994, 25: 596-600. 10.1051/apido:19940610.
Kubo T, Sasaki M, Nakamura J, Sasagawa H, Ohashi K, Takeuchi H, Natori S: Change in the expression of hypopharyngeal-gland proteins of the worker honeybees (Apis mellifera L.) with age and/or role. J Biochem. 1996, 119: 291-295. 10.1093/oxfordjournals.jbchem.a021237.
Ohashi K, Natori S, Kubo T: Change in the mode of gene expression of the hypopharyngeal gland cells with an age-dependent role change of the worker honeybee Apis mellifera L. Eur J Biochem. 1997, 249: 797-802. 10.1111/j.1432-1033.1997.t01-1-00797.x.
Kucharski R, Maleszka R, Hayward DC, Ball EE: A royal jelly protein is expressed in a subset of Kenyon cells in the mushroom bodies of the honey bee brain. Naturwissenschaften. 1998, 85: 343-346. 10.1007/s001140050512.
Feng M, Fang Y, Li J: Proteomic analysis of honeybee worker (Apis mellifera) hypopharyngeal gland development. BMC Genomics. 2009, 10: 645-10.1186/1471-2164-10-645.
Ueno T, Nakaoka T, Takeuchi H, Kubo T: Differential gene expression in the hypopharyngeal glands of worker honeybees (Apis mellifera L.) associated with an age-dependent role change. Zoolog Sci. 2009, 26: 557-563. 10.2108/zsj.26.557.
Albert Š, Bhattacharya D, Klaudiny J, Schmitzová J, Šimúth J: The family of major royal jelly proteins and its evolution. J Mol Evol. 1999, 49: 290-297. 10.1007/PL00006551.
Albert Š, Klaudiny J, Šimúth J: Molecular characterization of MRJP3, highly polymorphic protein of honeybee (Apis mellifera) royal jelly. Insect Biochem Mol Biol. 1999, 29: 427-434. 10.1016/S0965-1748(99)00019-3.
Kamakura M: Royalactin induces queen differentiation in honeybees. Nature. 2011, 473: 478-483. 10.1038/nature10093.
Huang CY, Chi LL, Huang WJ, Chen YW, Chen WJ, Kuo YC, Yuan CM, Chen CN: Growth stimulating effect on queen bee larvae of histone deacetylase inhibitors. J Agric Food Chem. 2012, 60: 6139-6149. 10.1021/jf300815b.
Asencot M, Lensky Y: The effect of sugars and juvenile hormone on the differentiation of the female honeybee larvae (Apis mellifera L.) to queens. Life Sci. 1976, 18: 693-699. 10.1016/0024-3205(76)90180-6.
Kaftanoglu O, Linksvayer TA, Page RE: Rearing honey bees, Apis mellifera, in vitro 1: effects of sugar concentrations on survival and development. J Insect Sci. 2011, 11: 96-
Spannhoff A, Kim YK, Raynal NJ, Gharibyan V, Su MB, Zhou YY, Li J, Castellano S, Sbardella G, Issa JP, Bedford MT: Histone deacetylase inhibitor activity in royal jelly might facilitate caste switching in bees. EMBO Rep. 2011, 12: 238-243. 10.1038/embor.2011.9.
Volny VP, Gordon DM: Genetic basis for queen–worker dimorphism in a social insect. Proc Natl Acad Sci USA. 2002, 99: 6108-6111. 10.1073/pnas.092066699.
Whitfield CW, Band MR, Bonaldo MF, Kumar CG, Liu L, Pardinas JR, Robertson HM, Soares MB, Robinson GE: Annotated expressed sequence tags and cDNA microarrays for studies of brain and behavior in the honey bee. Genome Res. 2002, 12: 555-566. 10.1101/gr.5302.
Thompson GJ, Kucharski R, Maleszka R, Oldroyd BP: Towards a molecular definition of worker sterility: differential gene expression and reproductive plasticity in honey bees. Insect Mol Biol. 2006, 15: 637-644.
Garcia L, Saraiva Garcia CH, Calábria LK, da Cruz Costa Nunes G, Sánchez Puentes A, Báo SN, Fontes W, Ricart CA, Salmen Espindola F, de Sousa Valle M: Proteomic analysis of honey bee brain upon ontogenetic and behavioral development. J Proteome Res. 2009, 8: 1464-1473. 10.1021/pr800823r.
Peixoto LG, Calábria LK, Garcia L, Capparelli FE, Goulart LR, de Sousa MV, Espindola FS: Identification of major royal jelly proteins in the brain of the honeybee Apis mellifera. J Insect Physiol. 2009, 55: 671-677. 10.1016/j.jinsphys.2009.05.005.
Hojo M, Kagami T, Sasaki T, Nakamura J, Sasaki M: Reduced expression of major royal jelly protein 1 gene in the mushroom bodies of worker honeybees with reduced learning ability. Apidologie. 2010, 41: 194-202. 10.1051/apido/2009075.
Randolt K, Gimple O, Geissendörfer J, Reinders J, Prusko C, Mueller MJ, Albert S, Tautz J, Beier H: Immune-related proteins induced in the hemolymph after aseptic and septic injury differ in honey bee worker larvae and adults. Arch Insect Biochem Physiol. 2008, 69: 155-67. 10.1002/arch.20269.
Woltedji D, Fang Y, Han B, Feng M, Li R, Lu X, Li J: Proteome analysis of hemolymph changes during larval to pupal development stage of honeybee worker (Apis mellifera ligustica). J Proteome Res. 2013, (doi: 10.1021/pr400519d) [epub ahead of print]
Chan QW, Melathopoulos AP, Pernal SF, Foster LJ: The innate immune and systemic response in honey bees to a bacterial pathogen, Paenibacillus larvae. BMC Genomics. 2009, 10: 387-10.1186/1471-2164-10-387.
Pfaffl MW, Tichopád A, Prgomet C, Neuvians TP: Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper – Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004, 26: 509-515.
Chan QW, Chan MY, Logan M, Fang Y, Higo H, Foster LJ: Honey bee protein atlas at organ-level resolution. Genome Res. 2013, (doi: 10.1101/gr.155994.113) [epub ahead of print]
Peiren N, de Graaf DC, Vanrobaeys F, Danneels EL, Devreese B, Van Beeumen J, Jacobs FJ: Proteomic analysis of the honey bee worker venom gland focusing on the mechanisms of protection against tissue damage. Toxicon. 2008, 52: 72-83. 10.1016/j.toxicon.2008.05.003.
Peiren N, Vanrobaeys F, de Graaf DC, Devreese B, Van Beeumen J, Jacobs FJ: The protein composition of honeybee venom reconsidered by a proteomic approach. Biochim Biophys Acta. 2005, 1752: 1-5. 10.1016/j.bbapap.2005.07.017.
Seeley TD: Adaptive significance of the age polyethism schedule in honeybee colonies. Behav Ecol Sociobiol. 1982, 11: 287-293. 10.1007/BF00299306.
Huang ZY, Otis GW, Teal PEA: Nature of brood signal activating the protein synthesis of hypopharyngeal gland in honey bees, Apis mellifera (Apidae: Hymenoptera). Apidologie. 1989, 20: 455-464. 10.1051/apido:19890601.
Naiem ES, Hrassnigg N, Crailsheim K: Nurse bees support the physiological development of young bees (Apis mellifera L.). J Com. Physiol B. 1999, 169: 271-279.
Hernández LG, Lu B, da Cruz GC, Calábria LK, Martins NF, Togawa R, Espindola FS, Yates JR, Cunha RB, de Sousa MV: Worker honeybee brain proteome. J Proteome Res. 2012, 11: 1485-93. 10.1021/pr2007818.
Gregorc A, Evans JD, Scharf M, Ellis JD: Gene expression in honey bee (Apis mellifera) larvae exposed to pesticides and Varroa mites (Varroa destructor). J Insect Physiol. 2012, 58: 1042-1049. 10.1016/j.jinsphys.2012.03.015.
Evans JD: Beepath: an ordered quantitative-PCR array for exploring honey bee immunity and disease. J Invertebr Pathol. 2006, 93: 135-139. 10.1016/j.jip.2006.04.004.
Ruijter JM, Ramakers C, Hoogaars WMH, Karlen Y, Bakker O, van den Hoff MJB, Moorman AFM: Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucl Acids Res. 2009, 37: e45-10.1093/nar/gkp045.
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucl Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.
We are particularly grateful to Holger Scharpenberg for handling the bees required for this study. AB was supported by a fellowship of the Prorectorate for Research and Young Academics of the Martin-Luther-University Halle-Wittenberg. The DFG provided financial support for chemicals and supplies (RFAM).
The authors declare that they have no competing interests.
AB, RFAM and SE conceived the study and designed the experimental set-up. AB executed the experiments. AB and SE analysed the data. AB wrote the first draft of the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Analysis of qPCR product specificity; mrjp1 – 89 bp, mrjp2 – 90 bp, mrjp3 – 102 bp, mrjp4 – 91 bp, mrjp5 – 80 bp, mrjp6 – 146 bp, mrjp7 – 95 bp, mrjp8 – 130 bp, mrjp9 – 103 bp, RpSa5 – 115 bp, arp1 – 120 bp. Figure S2. Two-way hierarchical clustering analysis heat map and dendrogram of mrjp gene expression data over all honeybee castes and groups. Branch lengths in dendrograms produced from cluster analysis correspond to the relative degree of similarity between branches. Differential gene expression is represented for all genes as a color gradient across all samples from deep blue (lowest) to light yellow (highest). The expression values were log transformed and visualized using MultiExperiment Viewer (MeV) version 4.9. Table S1. Primer characteristics for qPCR. Primer sequences for RpS5a and arp1 were adopted from  and . Table S2. Expression profile of a specific major royal jelly protein over all analyzed groups within honeybee body sections. Table S3. Expression profiles of all major royal jelly proteins in a specific group within honeybee body sections. Table S4. Amino acid sequence identities of all Apis mellifera MRJPs and the single MRJPs of Bombus terrestris (ADW82102.1), Megachile rotundata (XP_003708472.1), Camponotus floridanus (EFN61808.1) and Harpegnathos saltator (EFN81209.1). Sequence identities were determined using ClustalW2 (http://www.ebi.ac.uk/Tools/msa/clustalw2/). (PDF 1 MB)
About this article
Cite this article
Buttstedt, A., Moritz, R.F. & Erler, S. More than royal food - Major royal jelly protein genes in sexuals and workers of the honeybee Apis mellifera. Front Zool 10, 72 (2013). https://doi.org/10.1186/1742-9994-10-72
- Social insect
- Queen determination
- Caste determination
- Caste system