Skip to main content

More than royal food - Major royal jelly protein genes in sexuals and workers of the honeybee Apis mellifera



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 [1]. 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 [2]. 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 [3], 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) ([69]; for a review see [10]). 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 ([10] and references therein). Mrjps are also found in the primitively eusocial paper wasp Polistes canadensis [11], 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 [10]. 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 [14]. 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, 1625]. 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, 2631]. 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 [34] and MRJP3 [35] 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 [38] 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 [39], 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, 4044], and mrjp1 and 3 were shown to be expressed in drones (head, body, larvae) and queens (ovary, larvae) [7]. MRJP1 to 3 were also found in the haemolymph of larvae [45] with a significant higher amount than in pupae haemolymph [46] and they were shown to be down-regulated or depleted after infection with a severe pathogen of the honeybee, Paenibacillus larvae[47]. 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. [48], 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).

Figure 1
figure 1

Two-way hierarchical clustering analysis heat map and dendrogram of mrjp gene expression data over all honeybee castes, sexes 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. (a) Relative gene expression across all honeybee groups and body sections. The dendrogram for the clustering of the mrjps was generated using data of all three body sections and is exemplarily shown for the thorax. The complete heat map that served as base for the clustering of the genes is shown in Additional file 1: Figure S2. Relative gene expression values are shown in brackets. (b) Gene expression normalized to caged workers. The normalized gene expression values are shown in brackets.

Table 1 Overview of relative major royal jelly protein ( mrjp ) gene expression levels (normalized to RpS5a ) between honeybee castes, sexes and body sections
Figure 2
figure 2

Analyses of mrjp1 to 7 mRNA levels in heads of worker honeybees. The y-axis (log-scaled) indicates the relative levels of mRNA expression. Significant differences are marked by asterisks (* P < 0.05, *** P < 0.001). C, caged workers, light grey; N, hive nurses, grey; F, foragers, dark grey. Mrjp8 and 9 are not shown as these genes are not specifically up-regulated within the head.

Figure 3
figure 3

Analyses of mrjp9 mRNA levels throughout honeybee castes, sexes and body sections. The y-axis (log-scaled) indicates the relative levels of mRNA expression. Non-significant differences are marked by corresponding P values. All other between castes and sexes comparisons were significantly different (P < 0.01). H, head, light grey; T, thorax, grey; A, abdomen, dark grey.

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 [10]). 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. [49] MRJPs were shown to occur in various organs of the abdomen but not in the thoracic muscle and the thoracic salivary gland (see also [25]). Nevertheless, various MRJPs, except for MRJP4, 6 and 8, were identified in the honeybee’s nerve chord [49] that is crossing the thorax and abdomen [3]. 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 [10]. 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) [52]. 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) [9] 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. [9], especially because the protein contains a repetitive pentapeptide motive comprising several nitrogen rich amino acids [33]. MRJP3 was also associated with queen differentiation as specific post-translational modified isoforms of MRJP3 secreted into RJ lead to the development of queens [35].

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 [55] 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 [10]. 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. [44] 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 [44]. These results point to a function of MRJP7 primarily in the brain.

Hojo et al. [44] 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%, [9]) this supports the presence of an additional function of MRJP2 most likely in the brain.

MRJP1 is the most abundant protein of RJ [9] 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, 4044] and more specifically in the Kenyon cells of the brain’s mushroom bodies, cells involved in the formation and storage of memory and learning [29].

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 [10]). 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 [12]. 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

Honeybee samples

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)) [53]. 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).

Gene expression

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 [58] 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 [59] using RpS5a [57] 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 [60]. The Clustering Calculator software ( 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.



Caged workers






Major royal jelly protein




Queens with eggs


Queens with pupae


Virgin queens


Royal jelly.


  1. Wilson EO: The insect societies. 1971, Cambridge: Belknap Press of Harvard University Press

    Google Scholar 

  2. Winston ML: The biology of the honey bee. 1987, Cambridge: Harvard University Press

    Google Scholar 

  3. Snodgrass RE: Anatomy of the Honey Bee. 1910, Washington: Government Printing Office

    Book  Google Scholar 

  4. Von Planta A: Ueber den Futtersaft der Bienen. Hoppe-Seyl Zeitschr Physiol Chem. 1888, 12: 327-354.

    Google Scholar 

  5. Rembold H: Royal Jelly. Queen Rearing: Biological Basis and Technical Instruction, Chapter 2. Edited by: Ruttner F. 1983, Bucharest: Apimondia Publishing House, 35-42.

    Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Honeybee Genome Sequencing Consortium: Insights into social insects from the genome of the honeybee Apis mellifera. Nature. 2006, 443: 931-949. 10.1038/nature05260.

    Article  Google Scholar 

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

    Article  PubMed  Google Scholar 

  10. 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]

    Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  34. Kamakura M: Royalactin induces queen differentiation in honeybees. Nature. 2011, 473: 478-483. 10.1038/nature10093.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  46. 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]

    Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  49. 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]

    Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  52. Seeley TD: Adaptive significance of the age polyethism schedule in honeybee colonies. Behav Ecol Sociobiol. 1982, 11: 287-293. 10.1007/BF00299306.

    Article  Google Scholar 

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

    Article  CAS  Google Scholar 

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

    Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references


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

Author information

Authors and Affiliations


Corresponding author

Correspondence to Anja Buttstedt.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

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 [56] and [57]. 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 ( (PDF 1 MB)

Additional file 2:Statistical details for Tables S2 and S3.(XLS 142 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: