Malnutrition is a major factor affecting animal health, resistance to disease and survival. In honey bees (Apis mellifera), pollen, which is the main dietary source of proteins, amino acids and lipids, is essential to adult bee physiological development while reducing their susceptibility to parasites and pathogens. However, the molecular mechanisms underlying pollen's nutritive impact on honey bee health remained to be determined. For that purpose, we investigated the influence of pollen nutrients on the transcriptome of worker bees parasitized by the mite Varroa destructor, known for suppressing immunity and decreasing lifespan. The 4 experimental groups (control bees without a pollen diet, control bees fed with pollen, varroa-parasitized bees without a pollen diet and varroa-parasitized bees fed with pollen) were analyzed by performing a digital gene expression (DGE) analysis on bee abdomens.
Around 36, 000 unique tags were generated per DGE-tag library, which matched about 8, 000 genes (60% of the genes in the honey bee genome). Comparing the transcriptome of bees fed with pollen and sugar and bees restricted to a sugar diet, we found that pollen activates nutrient-sensing and metabolic pathways. In addition, those nutrients had a positive influence on genes affecting longevity and the production of some antimicrobial peptides. However, varroa parasitism caused the development of viral populations and a decrease in metabolism, specifically by inhibiting protein metabolism essential to bee health. This harmful effect was not reversed by pollen intake.
The DGE-tag profiling methods used in this study proved to be a powerful means for analyzing transcriptome variation related to nutrient intake in honey bees. Ultimately, with such an approach, applying genomics tools to nutrition research, nutrigenomics promises to offer a better understanding of how nutrition influences body homeostasis and may help reduce the susceptibility of bees to (less virulent) pathogens.
The availability of environmental nutrients is a key determinant of organism growth and survival. Therefore, signal transduction and the transcriptional network linking extracellular nutrients to the expression of metabolic genes have been widely studied in model organisms, like single cell organisms  and the fruit fly . In addition, nutrients in the diet play a crucial role in developing optimal immune response, such that a deficient or improper diet can have negative consequences on the susceptibility to a variety of pathogens [3, 4]. This relationship between health and nutrition has led to the development of a new scientific discipline, called nutrigenomics, that integrates high-throughput genomics tools with nutrition research [5, 6]. One of the main goals is to understand how nutrition influences metabolic pathways and homeostasis of the organism and ultimately develop "dietary-intervention strategies to recover normal homeostasis and to prevent diet-related diseases" .
The honey bee (Apis mellifera) is an important model for such studies, because its dietary requirements are well known  and the sequencing of its genome has been completed . If floral nectar (carbohydrates) represents the energetic fuel of honey bees, pollen provides the nutrients required for their internal organ development . Protein constitutes 2.5 to 61% of the dry mass of pollen , which is virtually the only source of protein naturally available. Nutrients reported also include lipids (1-20%), amino acids, starch, sterols, vitamins and minerals [10, 11]. Together these nutrients combine to make pollen nutrition one of the most important factors influencing the longevity of newly-emerged bees . In addition, proper quantity and quality of pollen not only reduce sensitivity to pesticides , but also help the bees to fight against pathogens, like microsporidia , bacteria  and virus , as pollen nutrition actually enhances honey bee immune function . However, the influence of pollen nutrition on honey bee health and resistance to disease cannot be fully understood without determining how nutrients act at the molecular level. This research gap coupled with the availability of new genomic tools in honey bees makes the development of a nutrigenomics approach in honey bees relevant.
Using this approach, we attempted to determine the pollen-induced molecular mechanisms that influence honey bee health. For that purpose, we characterized the transcriptome changes caused by pollen feeding in normal and sick bees. Sick bees were parasitized by the mite Varroa destructor, which is an ectoparasite that feeds on the hemolymph of immature and adult bees  and in this way infects bees with viruses . The consequences are dramatic since varroa parasitism provokes the deficiency of the honey bee physiological metabolism and immune system [20, 21]. Previous studies analyzed the transcriptional response to mite parasitism in prepupae  and pupae . Here, the transcriptome analysis was performed on whole abdomen of adult bees, which encompasses the midgut, where pollen digestion takes place, and the fat body, the main site of antimicrobial peptide synthesis and nutrient storage. Here, we used the digital gene expression (DGE) method to perform a deep transcriptome analysis of bees fed with pollen. This method combines the sequencing serial analysis of gene expression (SAGE) principle with the sequencing technology for generating a digital output proportional to the number of transcripts per mRNA . The benefit of not requiring presynthesized oligonucleotide probes (as in microarrays) allows the direct enumeration of transcripts, which is highly replicable, accurate and comparable across experiments [25, 26]. In addition, the quantification of weakly expressed genes can be analyzed which cannot be assessed using microarray. After generating the DGE-tag libraries we determined how pollen feeding influences metabolic pathways, genes affecting lifespan, immune functions and virus prevalence.
Results and discussion
Statistics of tag sequencing
Four DGE-tag libraries were generated from our experimental groups: control bees without a pollen diet (V-P-), control bees fed with pollen (V-P+), varroa-parasitized bees without a pollen diet (V+P-) and varroa-parasitized bees fed with pollen (V+P+). In total, around 14 × 106 tags (excluding the adaptor sequences) were sequenced from each library (see Table 1 for details). The libraries from varroa-parasitized bees contained more aligned tags (~7.3 × 106 tags) than non-parasitized bees (~4.5 × 106 tags). This higher number of aligned tags is due to the increase in tags matching virus sequences, which reached 1, 417, 481 in the V-P- library, 824, 212 in V-P+, 4, 258, 724 in V+P- and 5, 061, 614 in V+P+. These virus tags represent 32% of the total number of aligned tags in the V-P- library, 18% in V-P+, 60% in V+P- and 67% in V+P+. This clearly shows that varroa triggers a boost of pathogen loads in bees. This was likely induced by repeated varroa feeding on bees hemolymph that not only affect bee physiology but also contribute to the transmission of different viruses [19, 21, 27, 28]. In addition, we noticed that healthy (non-parasitized) bees, which were fed with pollen, had an overall lower load of pathogen tags than bees restricted to a sugar diet. However, few viruses were present at significant different levels between each treatment group (see the Effects of pollen feeding and varroa parasitism on virus prevalence).
Statistics of DGE sequencing
20, 397, 219
17, 130, 032
19, 857, 085
20, 633, 649
Total filtered tags
14, 817, 859
13, 508, 333
14, 346, 766
13, 867, 965
Total tags excluding adaptor sequences
14, 471, 228
13, 316, 644
14, 138, 414
13, 598, 573
Total aligned tags
4, 373, 417
4, 616, 713
7, 195, 133
7, 566, 291
Distinct aligned tags
Number of unique honey bee gene hit
Number of digital tags generated from the four libraries: control bees without a pollen diet (V-P-), control bees fed with pollen (V-P+), varroa-parasitized bees without a pollen diet (V+P-), varroa-parasitized bees fed with pollen (V+P+).
These distinct tags and their genomic count as well as raw data are available from NCBI Gene Expression Omnibus (GEO) database with the accession number: GSE25161. The number of distinct aligned tags was similar between each library (Table 1) and their distribution shows a consistent pattern across the 4 DGE libraries, with most of the distinct tags (34 to 40%) appearing only once and around 40% appearing between two and ten times (Additional file 1). The percentage of tags with a high count number dropped dramatically and only 1% of distinct tags have a count above 1, 000. Such pattern of DGE-Tag profiling was previously observed in honey bee heads  and in others insects and vertebrates [30, 31]. In total, around 8, 000 honey bee genes were maps from the four abdominal libraries (Table 1), which represents approximately 60% of the genes so far identified from the honey bee genome . As a comparison, 55 to 60% of the genes were also found to be expressed in the bee brain using DGE and microarray analysis [29, 32].
Differentially expressed genes and profiling of marker gene expression
The number of genes, whose expression level was significantly affected by pollen feeding or varroa-parasitism at an adjusted P-value for multiple comparison (Q-value) <0.01 is given in Table 2 and the lists are provided in Additional file 2. While a similar number of genes were up- and downregulated by pollen nutrition in healthy bees, varroa-parasitism had a major inhibitory effect on gene expression. This inhibition at the transcription level might explain the loss of weight and the decreased metabolism induced by varroa in emerging bees .
Genes differentially transcribed after pollen feeding or varroa parasitism
Genes were considered significantly differentially expressed between two treatments if their associated Q-value < 0.01.
To assess the reliability of our DGE-tag profiling, we checked whether the regulation of some candidates genes was consistent with the pattern of expression previously reported in earlier studies and performed reverse-transcriptase PCR (RT-PCR) on those genes. Vitellogenin (Vg) is a yolk protein taken up by developing oocytes and associated to egg production in queens [33, 34] but also has antioxidant functions that protect bees from oxidative stress and enhance longevity [35, 36]. Therefore, Vg level can be assimilated to an indicator of longevity. Since, pollen nutrition promotes the development of fat body [12, 17] and Vg is synthesized in the fat body , we expected an increase in Vg level in bees fed with pollen as previously found . That was both confirmed by the DGE and RT-PCR analysis (Additional file 2 and Figure 1). This suggests that one way by which pollen feeding increase longevity  is via the increase of Vg level. Interestingly, varroa parasitism decreased the Vg level in bees fed with pollen (see also ), but this level was still significantly higher than in parasitized bees restricted to a sugar diet.
Another interesting gene is malvolio (mvl), which encodes a manganese transmembrane transporter involved in sucrose responsiveness . In honey bees, foraging workers that are more responsive to sucrose have a higher level of mvl mRNA in the brain than young nurses taking care of the brood inside the hive . Pollen feeding induced an expression pattern opposite to Vg (Figure 1). Bees fed only with sugar had a higher level of mvl in the abdomen. Since old foragers feed mainly on carbohydrates and young nurses feed on both carbohydrates and pollen , emerging bees that are restricted to a sugar diet tend to develop a mvl expression profile similar to foragers as compared to bees fed with pollen, which are more similar to nurse bees. Indeed, as opposed to poor diet, a high protein diet might stimulate nurse tasks via the development of hypopharyngeal glands used for secreting larval food . Similarly, varroa parasitism induced an increase in mvl level (although not significant in bees fed with pollen), suggesting also a forager profile.
The expression level of the gene prophenoloxidase (PPO), involved in insect immunity , was also confirmed by RT-PCRs (Figure 1) and previous data. Indeed, varroa parasitism inhibited PPO expression, which was consistent with the immune suppression previously reported , but pollen feeding did not affect PPO level (see also  at the enzyme activity level).
Finally, we checked the expression level of the Drosophila ortholog spaetzle, which is a member of the Toll immune signaling pathway against fungi and bacteria . The DGE-tag profiling of this gene, confirmed by RT-PCR, showed that pollen feeding significantly increased the number of spaetzle transcripts in healthy bees (Figure 1), which suggests that pollen nutrients can enhance immune function in bees. The expression of this gene was reduced by varroa parasitism in bees fed with pollen, supporting again the immune suppression induced by varroa. Altogether, these results validated the DGE-tag profiling of the four libraries.
Functional analysis of differentially expressed genes
Gene Ontology analysis was performed to explore which functional components are regulated by pollen feeding and affected by varroa parasitism. Although the pollen wall, composed of sporopollenin and cellulose, is poorly digested by honey bee digestive enzymes, 50 to 98% of pollen contents can be extracted from ingested pollen, depending on pollen type and age of the bee . Bees of 9 days-old, close to the age of bees that we analyzed, can digest as much as 95% of pollen contents . The upregulation of many metabolic functions in bees fed with pollen highlights the high digestibility of pollen (Figure 2A). Indeed, in response to pollen nutrients, composed mostly of proteins and lipids, genes involved in proteolysis, lipid metabolism and coding for peptidases were upregulated. Carbohydrate metabolism was also upregulated by pollen feeding, suggesting a digestive response to pollen starches. Besides the metabolic categories, "oxidative phosphorylation" and "generation of precursor metabolites and energy" showed significant enrichment among the tested categories. This might not be surprising since nutrients are used to produce energy, like ATP via oxidative phosphorylation.
Groups of genes implicated in translational processes and ribosome biogenesis were downregulated by pollen feeding (Figure 2B). In all eukaryotes, ribosome biogenesis and translation regulation are essential for tissue growth and are regulated by food intake. Proper regulation of those processes are required to ensure cell and tissue homeostasis upon changes in nutrition and are inhibited during starvation for dietary nutrients . Therefore, an upregulation of ribosome biogenesis and translation was expected upon pollen feeding. However, the majority of pollen consumption and tissue development (e.g. hypopharyngeal glands, ovaries and fat body) occur within the first week after adult emergence . The prior completion of physiological development before our analysis could explain the downregulation of translational processes and ribosome biogenesis.
The effect of pollen nutrition on honey bee physiology was dramatically altered by varroa parasitism. Whereas groups of genes involved in ribosome biogenesis and translation were also downregulated, metabolic categories were no longer upregulated (Figure 2A). Only small sets of genes involved in proteolysis, peptidase activity and carbohydrates metabolism were still upregulated by pollen nutrition. In contrast most of the genes implicated in protein metabolism were downregulated in varroa-parasitized bees (Figure 2B), demonstrating that varroa mites inhibited the digestion and/or the use of proteins. In addition, the "lipid metabolic process", "oxidative phosphorylation" and "generation of precursor metabolites and energy" categories were no longer enriched. Metabolic processes were also found to be mostly downregulated by varroa parasitism in prepupae . Altogether, these results indicate that parasitized bees could not correctly assimilate and use the pollen nutrients required for their physiological development.
This was further confirmed by the effect of varroa on the bee transcriptome. Functional analysis showed that all of the previous categories, such as ribosome biogenesis, translation, metabolic functions and energy production, were inhibited by the varroa mite (Figure 2D). The inhibition by varroa of ribosome biogenesis and translation, essential to cell growth, might lead to lower tissue development. The lower concentrations of proteins and carbohydrates induced by varroa  might be explained by our results on the inhibition of metabolic genes and genes coding for endopeptidase. This reduction of metabolic reserves induced by varroa was reported at the emergence , therefore the downregulation of metabolic genes lasting 8 days after the emergence in bees fed with pollen, demonstrates that pollen nutrients are not able to compensate those losses.
A large set of genes were regulated by pollen feeding both in parasitized and non-parasitized bees and presented a similar expression pattern (Table 3 part 1). However, a subset of genes was upregulated by pollen feeding in healthy bees and downregulated in parasitized bees. A similar pattern was observed regarding the effect of varroa on bees fed with or without pollen (Table 3 part 2). In order to decipher the biological meaning of those opposite regulations, we performed a functional analysis of those subsets. Among the subset upregulated by pollen feeding in healthy bees but downregulated in parasitized bees, groups of genes are involved in ribosome biogenesis, translation, protein metabolism and energy production (Additional file 3 part A). The subset upregulated by varroa parasitism in bees fed without pollen but downregulated in pollen-fed bees, showed significant enrichment for genes coding for oxidoreductase and involved in amino acid metabolism (Additional file 3 part B). Those results confirmed the negative effect of varroa on those previously reported functions. Indeed, the energy metabolism and nutrient demand of the mite is high, utilizing up to 25% of the nutritional reserves of pupae . The decrease in metabolic reserves lasting in adults might come from the development of viruses transmitted by mites at the larvae or pupae stage.
Overlap of genes regulated in the abdomen by pollen feeding or varroa parasitism
χ2 = 190.1
χ2 = 74.6
P < 0.001
P < 0.001
1) Number of genes regulated by pollen feeding in bees parasitized (Varroa +) or not by varroa (Varroa -). 2) Number of genes regulated by varroa parasitism in bees fed with (Pollen +) or without pollen (Pollen -). ↑ and ↓ indicate genes that are up- or downregulated. Chi-square tests with Yates correction were performed to determine whether overlapping genes follow the same pattern of expression between the two conditions.
Similar results to the functional analysis were found with the pathways analysis (Additional file 4). Using the later analysis, we also studied the effect of pollen and varroa on amino acids. Indeed, honey bees require ten essential amino acids for their adult development: arginine, histidine, lysine, tryptophan, phenylalanine, methionine, threonine, leucine, isoleucine, and valine . This means those amino acids cannot be synthesized de novo by the bees and therefore must be supplied directly in the diet either as free amino acids or digested from proteins. We therefore looked at the amino acid metabolism pathways that were significantly enriched after pollen feeding and found that four of them (arginine, leucine, isoleucine, and valine) were upregulated (Additional file 4). Interestingly, five of those amino acid pathways were downregulated by varroa (lysine, tryptophan, leucine, isoleucine, and valine), confirming that varroa significantly reduces the development of bees. However, none of these essential amino acid pathways were significantly affected in varroa-parasitized bees fed with pollen compared to parasitized bees restricted to a sugar diet, suggesting that pollen nutrients compensated for the negative impacts of varroa on the metabolism of those amino acids.
Effects of pollen feeding and varroa parasitism on the nutrient sensing pathway
The insulin/TOR pathway is a conserved signaling cascade that functions as a nutrient sensing pathway by linking food-intake to animal growth and metabolism, including reproduction and lifespan [46, 47]. In honey bees, this pathway plays a major role in the regulation of aging of individuals [35, 48] but also at a social level by regulating the division of labor between worker bees [49, 50]. Less is known about the nutrient regulation of this pathway in honey bees. However, a recent study showed that the inhibition of Insulin receptor substrate expression in peripheral tissues (as opposed to central brain tissues) made forager bees collect pollen rather than nectar . By analyzing the DGE-tag profiling of peripheral tissues (abdomen) we found that most of the genes from the mTOR signaling pathways were upregulated by pollen nutrition in healthy bees (Figure 3). The insulin-like peptide 2 (ILP-2) expression notably increased in response to pollen nutrition, as it was found in response to amino acids , which might promote nurse physiology that extends bee longevity . Altogether these results show that the insulin/TOR pathway responds positively to pollen nutrition and might trigger the pollen-induced physiological development. The opposite effect was observed regarding ILP-1, with a higher expression under poor nutrition conditions . However, those inverse patterns might actually be explained by the functional differences of ILP-1 and ILP-2 
. In varroa-parasitized bees, pollen nutrition tended to inhibit this pathway. This could be caused by the higher prevalence of viruses in parasitized bees fed with pollen compared to parasitized bees fed only with sugar (see Effects of pollen feeding and varroa parasitism on virus prevalence). Indeed, varroa parasitism (coupled with a higher viral load) induced a similar pattern of the TOR pathway (Figure 3). The parasite inhibits protein metabolism and decreases the nutrient reserve of bees (see above), which might be explained by the downregulation of the TOR pathway mediating tissue growth in response to nutrient intake. This is further confirmed by the inhibition of slimfast by varroa parasitism in bees fed with pollen. The downregulation of this gene, which encodes an amino acid transporter and is a downstream target of TOR pathway, causes a deficient global growth in Drosophila similar to that seen in individuals raised under poor nutritional conditions [46, 53].
TOR is known to respond to the presence of amino acids and induces an upregulation of ribosome biogenesis, translation [2, 54–57] and energy metabolism [58, 59] required for tissue growth. Interestingly, the downregulation of TOR pathway was followed by the downregulation of each biological process (Figure 3), indicating that the same mechanisms of the nutrient sensing pathway are found in honey bees. The upregulation of the pathway in pollen-fed bees was followed by an increase in energy metabolism but the ribosome biogenesis and translation were both inhibited. As explained above, it is possible that tissue development that requires the activation of translational processes was already completed.
The expression level of genes did not always match well to the activator/inhibitor functions established in the Drosophila TOR pathway. The expression patterns correspond to global DGE-tag profilings of whole bee abdomens but many genes possess tissue-specific functions. Therefore, targeting specific tissues like the fat body, which has important storage functions associated with nutrition, should provide more accurate expression data for pathway analysis.
Effects of pollen feeding and varroa parasitism on genes affecting lifespan and immune genes
Since pollen nutrition and varroa increases and reduces the longevity of bees [12, 60], respectively, we looked at the expression pattern of genes that have been found to affect longevity in Drosophila. Those genes were found by using the functional analysis from DAVID 6.7 bioinformatic resources  and GoToolbox . We found 20 genes including Vg, known to affect longevity, to be differentially expressed between the different treatments. A mutation at the loci of each of those genes causes a decrease in longevity except for Tor, whose inhibition increase longevity . Among those genes, we expected to find a larger set to be upregulated by pollen feeding than downregulated and inversely regarding varroa effect. According to our expectation, 7 out of 10 genes regulated by pollen intake were upregulated in non-parasitized bees (Figure 4A). For example, Sod and Trxr-1, which encodes a cytoplasmic Cu-Zn superoxide dismutase and a thioredoxin reductase, respectively , are known to increase lifespan [65, 66], thanks to the degradation of superoxide radicals and H2O2 or the neutralization of reactive oxygen species. The positive effect of pollen feeding on the expression pattern of both genes was confirmed by reverse-transcriptase PCR (RT-PCR) analysis performed on a different set of experimental bees (n = 8 pools of bees per treatment, Additional file 5). Those genes might be involved in the pollen-induced molecular mechanisms that enhance life expectancy of honey bees. However, Sod 2 (mitochondrial Mn superoxide dismutase), which is also involved in degradation of superoxide radicals and H2O2 , was not affected by pollen feeding (Additional file 5), contrary to the downregulation found with the DGE analysis. Due to the higher level of Vg (Figure 1), a longer lifespan is expected in parasitized bees supplied with pollen compared to parasitized bees fed only with sugar (Figure 4A). However, except Vg, genes that were regulated by pollen intake were all downregulated. Whether pollen nutrition can affect varroa resistance remains to be tested but it might be possible that the effect of pollen on the genes affecting longevity is compensated by the increase in Vg. The picture was clearer with regards to the effect of varroa on the 19 genes affecting longevity (Figure 4A and Additional file 5). All of them were downregulated (including Sod 2 in presence of pollen, Additional file 5) supporting a decrease of longevity induced by varroa.
Next, we examined the effect of pollen intake and varroa on the expression pattern of immune genes, that have been previously annotated . Few immune genes were upregulated by pollen feeding in healthy bees (Figure 4B). Three of them encoded antimicrobial peptides (AMPs) (Lyzozyme-2 and -3, Defensin-1), which might be linked to a higher content of fat body, the major site of synthesis of AMPs. Interestingly, in the absence of immune challenge (healthy bees), Spaetzle and the gene coding for the peptidoglycan recognition protein (PGRP)-LC, which are activators of the Toll and Imd pathways respectively, were both upregulated by protein feeding. The upregulation of those key immune genes confirmed by RT-PCR analysis (Figure 1 for Spaetzle and Additional file 5 for Defensin-1 and PGRP-LC) suggests that bees might be able to mount an efficient immune response to pathogen infections. In parasitized bees, only Imd, a specific gene for antimicrobial defense, was upregulated by pollen feeding. A slight increase of Imd expression was also observed with RT-PCR analysis although not significant (p = 0.058, Additional file 5). In healthy bees, Imd expression did not change after pollen feeding (Additional file 5). Since, protein and amino acid nutrition can enhance immune functions in humans [4, 67] and insects [68, 69], an upregulation of immune genes was expected. The functional analysis showed that in varroa-parasitized bees fed with pollen, a significant number of genes involved in protein metabolism were downregulated compared to parasitized bees fed only with sugar (Figure 2B). A deficiency in protein metabolism might have affected the immune system of those bees. This immunosupression could be explained by the higher virus load in parasitized bees fed with pollen compared to parasitized bees fed only with sugar (see Effects of pollen feeding and varroa parasitism on virus prevalence). However, when fed with pollen, it is possible that parasitized bees invest more energy in other defense mechanisms than immunity which is inhibited by varroa, e.g. the antioxidant system (higher Vg level in bees fed with pollen).
Varroa induced the inhibition of most of the immune genes in bees fed with or without pollen (including PGRP-LC, Additional file 5). However, a few immune genes were upregulated, like the AMPs (Defensin-1, Hymenoptaecin, Apidaecin-1 and Lyzozyme-3) and cact-1, which controls the production on AMPs . The upregulation of Defensin-1 by varroa parasitism was confirmed by RT-PCR analysis (Additional file 5). Those genes, except Lyzozyme-3 have already been found to be upregulated by a immune challenge . The mechanisms triggering the immunosupression are still misunderstood but one explanation might be that the deficient protein metabolism induced by varroa and associated viruses (see above) block the development of immune function.
Immune genes that were upregulated by pollen feeding and varroa (Figure 4B) belonged to the antibacterial peptide family or were key activators of immune pathways. Since the characterization of immune defense in honey bee is still at his infancy, it is possible that the others immune genes don't actually play a major role in immune response in bees. Indeed, genes from immune signaling pathway have usually different roles, like developmental functions.
Effects of pollen feeding and varroa parasitism on virus prevalence
Using the DGE-tag libraries from honey bee abdomen samples, we performed a metagenomic analysis to determine the effects of pollen feeding and varroa on the prevalence of 10 honey bee viruses: Chronic bee paralysis virus RNA 1, Chronic bee paralysis virus RNA 2, Sacbrood virus, Deformed wing virus, Black queen cell virus, Acute bee paralysis virus, Kashmir bee virus, Kakugo virus, Varroa destructor virus 1 and Israel acute paralysis virus. All of them were identified in our DGE-tag libraries but very low numbers of tag copies of the Chronic bee paralysis virus RNA 1, Chronic bee paralysis virus RNA 2, Sacbrood virus, Black queen cell virus, Acute bee paralysis virus, Kashmir bee virus and Israel acute paralysis virus were found (from 1 to 65, see data deposited in GEO database: GSE25161) and no difference in their abundance was detected between treatments. However, tags mapping the deformed wing virus (DWV), Varroa destructor virus (VDV) and Kakugo virus (KV) were highly prevalent and significantly affected by the different treatments (Table 4). Varroa destructor is an efficient vector of bee viruses by promoting their transmission and development (DWV [19, 21], VDV , KV ) but the presence of viruses in control bees can be explained by a vertical transmission from queens and drones to worker offspring through fertilized eggs .
Effects of pollen feeding and varroa parasitism on the prevalence of honey bee virus
Number of tag copies
Tag odds ratio
Deformed winged virus-9489
1, 653, 968
2, 400, 089
Deformed winged virus-4761
Varroa destructor virus-5865
Varroa destructor virus-9462
Varroa destructor virus-9753
1, 202, 745
2, 514, 691
2, 556, 286
Virus tags that were differentially expressed between two conditions (Q-value < 0.01) in pools of 12 bees. Numbers of tag copies and odds ratio are indicated. ns stands for "not significant". Each tag ID is composed of the corresponding virus name and tag position in the sequence (indicated by the last numbers).
As expected, varroa provoked a boost of virus development (Table 4). Since AMPs have a broad-spectrum antimicrobial activity against bacteria, yeast, fungi, parasites but also viruses , the upregulation of AMPs in varroa-parasitized bees (Figure 4B) might represent the immune response to the high viral infection. Different studies showed that pollen nutrients can help the bees to fight against parasites and pathogens [14, 15]. For example, a recent study found that pollen nutrition reduce the prevalence of DWV in normal adult bees . We obtained similar results in healthy (non-parasitized) bees with pollen nutrition reducing the load of VDV. However, we observed the inverse pattern with DWV and KV. Viruses utilize cell machinery to multiply, so higher number of cells available could contribute also to a better development of viral population. For example, DWV has been found to multiply in the bee fat body  and this tissue is more developed in pollen-fed bees. We also found that pollen feeding increased viruses prevalence in varroa-parasitized bees, indicating pollen nutrition is not effective for limiting virus loads in heavily infected bees. Larger physiological machinery might, to the contrary, favor their multiplication, but it might also help bees to withstand viral infection.
With this study, we showed that nutrigenomics, originally developed for human health purposes, is also promising for deciphering the nutritive influence on metabolic pathways in insects. More specifically, this work is the first step towards understanding how pollen nutrients shapes honey bee health. While pollen nutrition enhances macromolecule metabolism and activates the TOR pathway required for the tissue growth and development, it also contains nutrients that stimulate the expression of genes involved in longevity, like genes coding for antioxidants (vitellogenin and superoxide dismutase). However, the negative impacts of varroa on the bee metabolism and immune functions could not be reversed by pollen feeding. By inhibiting the macromolecule metabolism (notably proteins), varroa parasitism prevented bees from accessing the beneficial effects of pollen. This demonstrates that pathology associated to varroa is extremely virulent and difficult to reverse, probably due to the multiplication of multiple viruses transmitted by the mite.
In the future, studying more specific tissues (like fat body) and particular dietary components will hopefully provide a more complete and comprehensible picture of the beneficial or adversarial effects of pollen components in bees. Finally, given that honey bees are victims of a wide range of pathogens that might be less virulent than varroa, developing dietary intervention based on knowledge of nutritional requirement and status can be useful to prevent, mitigate or cure chronic disease.
The experiment was performed with local hybrid honey bees (A. m. ligustica/A. m. mellifera). In order to obtain a large number of varroa-infected bees, an infested colony that was not treated with miticide was selected. The multiple-mated queen was placed in a queen cage to stop egg-laying for 10 days. After 10 days, all the brood produced before the queen confinement was sealed and the mites had no brood cell to parasitize. Then, the queen was transferred to a queen-excluder for 3 days with an empty frame. The brood produced by the queen on this frame was then uniformly parasitized by varroa mites. Three weeks after the queen laid eggs in the queen-excluder, the frame was removed from the colony and the sealed cells were opened carefully to remove the adult workers just before emergence and to check for the presence of reproducing varroa mites. Varroa-parasitized bees with deformed wings were discarded due to their extremely short lifespan.
Varroa-parasitized and non-parasitized bees (control) were placed in distinct cages in the dark at 32°C and 70% RH and fed ad libitum with candy (30% honey, 70% powdered sugar) and water. To simulate as much as possible the colony rearing conditions, cages were exposed to a Beeboost® (Pherotech), releasing one queen-equivalent of queen mandibular pheromone per day. Half of the groups were also fed with a polyfloral pollen diet, which was replaced every day and prepared with fresh pollen obtained from our colonies in Avignon (France) and 1/10 water. Thus we obtained 4 experimental groups: control bees without a pollen diet (V-P-), control bees fed with pollen (V-P+), varroa-parasitized bees without a pollen diet (V+P-) and varroa-parasitized bees fed with pollen (V+P+). Since their pollen consumption diminishes when bees are 8-10 days old , bees from each group were flash frozen in liquid nitrogen at the age of 8 days. One group (pool of 12 bees) per treatment was used for comparing the 'average' responses to the different treatments with the digital transcriptome analysis. Single replicate analyses have generated consistent results in previous studies [30, 31, 74], but confirmation of some of the results obtained from digital expression analysis was performed with RT-PCR on four new experimental replicates (see below). The quantities of pollen consumed in the V-P+ and V+P+ groups were similar (8.2 ± 2 and 7.7 ± 2.6 mg per bee and per day, respectively).
Twelve abdomens per group were homogenized in Trizol (Invitrogen cat. No. 15596-026) and RNA extraction was carried out as indicated in the Qiagen RNeasy kit (Qiagen cat. No. 74104) for total RNA with on-column DNase I treatment (Qiagen cat. No. 79254).
Sequencing libraries were constructed using Illumina's DGE Tag Profiling DpnII Sample Prep kit (Illumina cat. No. FC-102-1007) according to the manufacturer's instructions (cat. No. 1004241 rev.A). Briefly, polyA+ RNA were isolated from 2 μg of total RNA using Sera-magnetic oligo-dT beads. Reverse transcription was performed on immobilized polyA+ RNA using SuperScript II Reverse Transcriptase (Invitrogen cat. No. 18064-14). Second strand was synthesized using DNA Polymerase I. Immobilized cDNAs were digested with DpnII. Digested fragments were washed off, and the retained cDNA fragments were ligated to GEX DpnII Adapter 1. Tags were released by MmeI digestion generating 20 bp tags that are dephosphorylated using CIAP and ligated to GEX Adapter 2. The adapter-ligated cDNA tags were enriched using 15 cycles of PCR, and purified on a 6% Novex TBE PAGE gel (Invitrogen cat. No. EC6265Box). 85-bp fragments were excised from the gel, eluted and precipitated. Tags were quantified on an Agilent Bioanalyzer using DNA1000 chip (Agilent cat. No. 5067-1504). Each library was diluted to 10 nM, denatured, and diluted again to 8 pM. 100 μl of the diluted library were hybridized on a lane of an Illumina's Flow Cell. Clustering and 36 cycles sequencing were performed according to Illumina's instructions.
Analysis and mapping of DGE tags
We first filtered the raw data via the default Illumina pipeline quality filter, which calculates a Chastity score as the ratio of the highest of the four (base type) intensities to the sum of highest two. This score is used to remove clusters with low signal to noise ratio (often caused by clusters being too close to each other so their signals bleed into one another). The sequences were passed through the chastity filter of value > = 0.6. Then, adapters were trimmed from tag sequences by an in-house script. Reads displaying sequence ambiguity were excluded from further analysis. For each library, read counts were recorded for every unique sequence. Tag sequences were mapped to the honey bee genome (BeeBase Genome - release 2 http://genomes.arc.georgetown.edu/drupal/beebase/?q=home) with Eland v1.4, and only perfect tag matches of coding sequences were kept for statistical analysis. Mapping was also performed on different sequences of honey bee virus genomes (Chronic bee paralysis virus RNA 1: GenBank EU122229, Chronic bee paralysis virus RNA 2: EU122230, Sacbrood virus: AF092924, Deformed wing virus: AJ489744, Black queen cell virus: AF183905, Acute bee paralysis virus: AF150629, Kashmir bee virus: AY275710, Kakugo virus: AB070959, Varroa destructor virus 1: AY251269 and Israel acute paralysis virus: EF219380).
In the case of tag census, the probability distribution of the measured parameter 'tag abundance' is known a priori and follows a hypergeometric distribution. Accordingly, the probability that 2 census of the same transcript in 2 different conditions are significantly different can be calculated without replicate. Differentially expressed tags were therefore identified by a Fisher's exact test with Benjamini-Hochberg method, which controls the expected proportion of incorrectly rejected null hypotheses (type I errors), to compute the False Discovery Rate and a cut-off Q-value of 0.01 as implemented in the SAGE Genie resource http://cgap.nci.nih.gov/SAGE. The expression pattern of the tags was based on their sequence odds ratio between two pools A and B: (sequences in A/total sequences in pool A)/(sequences in B/total sequences in pool B).
We explored whether any particular molecular functions or biological processes from the differentially expressed gene lists were represented by larger numbers of genes than expected on the basis of chance. Drosophila melanogaster orthologs were identified by reciprocal best BLASTX match to bee genes, and Gene Ontology (GO) terms were assigned based on annotation of Drosophila genes. GOToolBox  was used to identify overrepresented terms (hypergeometric tests followed by the Benjamini Hochberg correction for multiple testing). Enrichment analysis of bio-pathways were performed using DAVID 6.7 bioinformatic resources .
Verification by reverse-transcriptase PCR
RT-PCR was performed on selected genes known to be affected by nutrition and varroa parasitism (vitellogenin (Vg), malvolio (mvl), Prophenoloxidase (PPO) and Spaetzle (spz)), three genes involved in lifespan (Superoxide dismutase (Sod), Superoxide dismutase 2 (Sod2), Thioredoxin reductase 1(Trxr-1)) and three immune genes (Peptidoglycan recognition protein LC (PGRP-LC), defensin1, Immune deficiency (Imd)). The RT-PCR analysis was carried out on a different set of bees obtained by using the same experimental procedure (see above). Four experimental replicates per treatment were performed and 2 pools of 3 bees per replicate (cage) were analyzed (n = 8 pools/treatment). The transcript abundance was measured for with an M × 3000 P QPCR Systems (Agilent) and the SYBR green detection method (Agilent cat. No. 600828). RT-PCR values of the selected genes were normalized by using an exogenous human RNA control (β-actin) spiked during cDNA synthesis performed with the SuperScript III (Invitrogen, France). Relative expression was calculated by raising 2 to the power of the difference in Ct values. Primer sequences are listed in Additional file 6.
We thank L. Journot for his advice in the design and analysis of the molecular studies and members of the lab for comments that improve the manuscript. This work was funded by a BEEDOC grant (FP7, RTD REG/E.4(2009)D/561221) and C. Alaux was supported by an INRA young researcher position (INRA SPE department).
INRA, UMR 406 Abeilles et Environnement, Domaine Saint-Paul
Rinderer TE, Rothenbuhler WC, Gochnauer TA: The influence of pollen on the susceptibility of honey-bee larvae toBacilluslarvae.J Invertebr Pathol 1974, 23:347–350.PubMedView Article
Degrandi-Hoffman G, Chen Y, Huang E, Huang MH: The effect of diet on protein concentration, hypopharyngeal gland development and virus load in worker honey bees (Apis melliferaL.).J Insect Physiol 2010, 56:1184–1191.PubMedView Article
Alaux C, Ducloz F, Crauser D, Le Conte Y: Diet effects on honeybee immunocompetence.Biol Lett 2010, 6:562–565.PubMedView Article
Sammataro D, Gerson U, Needham G: Parasitic mites of honey bees: life history, implications, and impact.Annu Rev Entomol 2000, 45:519–548.PubMedView Article
Bowen-Walker PL, Martin SJ, Gunn A: The transmission of deformed wing virus between honeybees (Apis melliferaL.) by the ectoparasitic mite Varroa jacobsoni Oud.J Invertebr Pathol 1999, 73:101–106.PubMedView Article
Bowen-Walker PL, Gunn A: The effect of the ectoparasitic mite, Varroa destructor on adult worker honeybee (Apis mellifera) emergence weights, water, protein, carbohydrate, and lipid levels.Entomol Exp Appl 2001, 101:207–217.View Article
Yang X, Cox-Foster DL: Impact of an ectoparasite on the immunity and pathology of an invertebrate: evidence for host immunosuppression and viral amplification.Proc Natl Acad Sci USA 2005,102(21):7470–7475.PubMedView Article
Zhang Y, Liu X, Zhang W, Han R: Differential gene expression of the honey beesApis melliferaandA. ceranainduced byVarroa destructorinfection.J Insect Physiol 2010,56(9):1207–1218.PubMedView Article
Navajas M, Migeon A, Alaux C, Martin-Magniette M, Robinson G, Evans J, Cros-Arteil S, Crauser D, Le Conte Y: Differential gene expression of the honey beeApis melliferaassociated with Varroa destructor infection.BMC Genomics 2008, 9:301.PubMedView Article
Anisimov SV: Serial Analysis of Gene Expression (SAGE): 13 years of application in research.Curr Pharm Biotechnol 2008,9(5):338–350.PubMedView Article
Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y: RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays.Genome Res 2008,18(9):1509–1517.PubMedView Article
't Hoen PA, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RH, de Menezes RX, Boer JM, van Ommen GJ, den Dunnen JT: Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms.Nucleic Acids Res 2008, 36:e141.PubMedView Article
Ongus JR, Peters D, Bonmatin JM, Bengsch E, Vlak JM, van Oers MM: Complete sequence of a picorna-like virus of the genus Iflavirus replicating in the miteVarroa destructor.J Gen Virol 2004,85(Pt 12):3747–3755.PubMedView Article
Fujiyuki T, Ohka S, Takeuchi H, Ono M, Nomoto A, Kubo T: Prevalence and phylogeny of Kakugo virus, a novel insect picorna-like virus that infects the honeybee (Apis melliferaL.), under various colony conditions.J Virol 2006,80(23):11528–11538.PubMedView Article
Liu F, Li W, Li Z, Zhang S, Chen S, Su S: High-abundance mRNAs inApis mellifera: comparison between nurses and foragers.J Insect Physiol 2011, 57:(274–279).PubMedView Article
Wang XW, Luan JB, Li JM, Bao YY, Zhang CX, Liu SS: De novo characterization of a whitefly transcriptome and analysis of its gene expression during development.BMC Genomics 2010, 11:400.PubMedView Article
Wang Z, Dong D, Ru B, Young RL, Han N, Guo T, Zhang S: Digital gene expression tag profiling of bat digits provides robust candidates contributing to wing formation.BMC Genomics 2010, 11:619.PubMedView Article
Alaux C, Le Conte Y, Adams HA, Rodriguez-Zas S, Grozinger CM, Sinha S, Robinson GE: Regulation of brain gene expression in honey bees by brood pheromone.Genes Brain Behav 2009,8(3):309–319.PubMedView Article
Engels E: Occurrence and significance of vitellogenins in female castes of social Hymenoptera.Am Zool 1974, 14:1229–1237.
Tanaka ED, Hartfelder K: The initial stages of oogenesis and their relation to differential fertility in the honey bee (Apis mellifera) castes.Arthropod Struct Dev 2004,33(4):431–442.PubMedView Article
Corona M, Velarde RA, Remolina S, Moran-Lauter A, Wang Y, Hughes KA, Robinson GE: Vitellogenin, juvenile hormone, insulin signaling, and queen honey bee longevity.Proc Natl Acad Sci USA 2007,104(17):7128–7133.PubMedView Article
Seehuus SC, Norberg K, Gimsa U, Krekling T, Amdam GV: Reproductive protein protects functionally sterile honey bee workers from oxidative stress.Proc Natl Acad Sci USA 2006,103(4):962–967.PubMedView Article
Bitondi MMG, Simões ZLP: The relationship between level of pollen in the diet, vitellogenin and juvenile hormone titres in AfricanizedApis melliferaworkers.J Apic Res 1996, 35:27–36.
Amdam GV, Hartfelder K, Norberg K, Hagen A, Omholt SW: Altered physiology in worker honey bees (Hymenoptera: Apidae) infested with the miteVarroa destructor(Acari: Varroidae): A factor in colony loss during overwintering?J Econ Entomol 2004, 97:741–747.PubMedView Article
Orgad S, Nelson H, Segal D, Nelson N: Metal ions suppress the abnormal taste behavior of theDrosophilamutant malvolio.J Exp Biol 1998,201(Pt 1):115–120.PubMed
Ben-Shahar Y, Dudek NL, Robinson GE: Phenotypic deconstruction reveals involvement of manganese transporter malvolio in honey bee division of labor.J Exp Biol 2004,207(Pt 19):3281–3288.PubMedView Article
Crailsheim K, Schneider LHW, Hrassnigg N, Bühlmann G, Brosch U, Gmeinbauer R, Schöffmann B: Pollen consumption and utilization in worker honeybees (Apis melliferacarnica): dependence on individual age and function.J Insect Physiol 1992, 38:409–419.View Article
Evans JD, Aronstein K, Chen YP, Hetru C, Imler JL, Jiang H, Kanost M, Thompson GJ, Zou Z, Hultmark D: Immune pathways and defence mechanisms in honey beesApis mellifera.Insect Mol Biol 2006, 15:645–656.PubMedView Article
Pernal SF, Currie RW: Pollen quality of fresh and 1-year-old single pollen diets for worker honey bees (Apis melliferaL.).Apidologie 2000, 31:387–409.View Article
Garedew A, Schmolz E, Lamprecht I: The energy and nutritional demand of the parasitic life of the miteVarroa destructor.Apidologie 2004, 35:419–430.View Article
de Groot AP: Protein and amino acid requirements of the hone y bee (Apis mellificaL.).Physiol Comp Oecol 1953, 3:197–285.
Colombani J, Raisin S, Pantalacci S, Radimerski T, Montagne J, Leopold P: A nutrient sensor mechanism controlsDrosophilagrowth.Cell 2003,114(6):739–749.PubMedView Article
Oldham S, Hafen E: Insulin/IGF and target of rapamycin signaling: a TOR de force in growth control.Trends Cell Biol 2003,13(2):79–85.PubMedView Article
Münch D, Amdam GV: The curious case of aging plasticity in honey bees.FEBS Lett 2010,584(12):2496–2503.PubMedView Article
Wang Y, Amdam GV, Rueppell O, Wallrichs M, Fondrk MK, Kaftanoglu O, Page REJ: PDK1 and HR46 gene homologs tie social behavior to ovary signals.PLoS ONE 2009, 4:e4899.PubMedView Article
Ament SA, Corona M, Pollock HS, Robinson GE: Insulin signaling is involved in the regulation of worker division of labor in honey bee colonies.Proc Natl Acad Sci USA 2008, 105:4226–4231.PubMedView Article
Wang Y, Mutti NS, Ihle KE, Siegel A, Dolezal AG, Kraftanoglu O, Amdam GV: Down-regulation of honeybee IRS gene biases behavior toward food rich in protein.PLoS Genet 2010, 6:e1000896.PubMedView Article
Amdam GV: Social context, stress, and plasticity of aging.Aging Cell 2011,10(1):18–27.PubMedView Article
Hennig KM, Colombani J, Neufeld TP: TOR coordinates bulk and targeted endocytosis in theDrosophila melanogasterfat body to regulate cell growth.J Cell Biol 2006,173(6):963–974.PubMedView Article
Li L, Edgar BA, Grewal SS: Nutritional control of gene expression inDrosophilalarvae via TOR, Myc and a novel cis-regulatory element.BMC Cell Biol 2010, 11:7.PubMedView Article
Hay N, Sonenberg N: Upstream and downstream of mTOR.Genes Dev 2004,18(16):1926–1945.PubMedView Article
Grewal SS, Evans JR, Edgar BA: DrosophilaTIF-IA is required for ribosome synthesis and cell growth and is regulated by the TOR pathway.J Cell Biol 2007, 179:1105–1113.PubMedView Article
Guertin DA, Guntur KV, Bell GW, Thoreen CC, Sabatini DM: Functional genomics identifies TOR-regulated genes that control growth and division.Curr Biol 2006,16(10):958–970.PubMedView Article
Tiefenbock SK, Baltzer C, Egli NA, Frei C: TheDrosophilaPGC-1 homologue Spargel coordinates mitochondrial activity to insulin signalling.Embo J 2009,29(1):171–183.PubMedView Article
Baltzer C, Tiefenbock SK, Marti M, Frei C: Nutrition controls mitochondrial biogenesis in theDrosophilaadipose tissue through Delg and cyclin D/Cdk4.PLoS ONE 2009,4(9):e6935.PubMedView Article
Schneider P, Drescher W: The influence ofVarroa jacobsoniOud. on weight, development of weight and hypopharyngeal glands, and longevity ofApis melliferaL.Apidologie 1987, 18:101–110.View Article
Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources.Nat Protoc 2009, 4:44–57.View Article
Martin D, Brun C, Remy E, Mouren P, Thieffry D, Jacq B: GOToolBox: functional analysis of gene datasets based on Gene Ontology.Genome Biol 2004,5(12):R101.PubMedView Article
Kapahi P, Zid BM, Harper T, Koslover D, Sapin V, Benzer S: Regulation of lifespan inDrosophilaby modulation of genes in the TOR signaling pathway.Curr Biol 2004,14(10):885–890.PubMedView Article
Corona M, Hughes KA, Weaver DB, Robinson GE: Gene expression patterns associated with queen honey bee longevity.Mech Ageing Dev 2005,126(11):1230–1238.PubMedView Article
Sun J, Tower J: FLP recombinase-mediated induction of Cu/Zn-superoxide dismutase transgene expression can extend the life span of adultDrosophila melanogasterflies.Mol Cell Biol 1999,19(1):216–228.PubMed
Orr WC, Mockett RJ, Sohal RS: Overexpression of glutathione reductase extends survival in transgenicDrosophila melanogasterunder hyperoxia but not normoxia.Faseb J 1999,13(13):1733–1742.PubMed
Lee KP, Simpson SJ, Wilson K: Dietary protein-quality influences melanization and immune function in an insect.Funct Ecol 2008, 22:1052–1061.View Article
Lee KP, Cory JS, Wilson K, Raubenheimer D, Simpson SJ: Flexible diet choice offsets protein costs of pathogen resistance in a caterpillar.Proc R Biol Sci 2006,273(1588):823–829.View Article
Lemaitre B, Nicolas E, Michaut L, Reichhart JM, Hoffmann JA: The dorsoventral regulatory gene cassette spatzle/Toll/cactus controls the potent antifungal response inDrosophilaadults.Cell 1996,86(6):973–983.PubMedView Article
Yue C, Schroder M, Gisder S, Genersch E: Vertical-transmission routes for deformed wing virus of honeybees (Apis mellifera).J Gen Virol 2007,88(Pt 8):2329–2336.PubMedView Article
Vizioli J, Salzet M: Antimicrobial peptides from animals: focus on invertebrates.Trends Pharmacol Sci 2002,23(11):494–496.PubMedView Article
Fievet J, Tentcheva D, Gauthier L, de Miranda J, Cousserans F, Colin ME, Bergoin M: Localization of deformed wing virus infection in queen and droneApis melliferaL.Virol J 2006, 3:16.PubMedView Article
Bradford JR, Hey Y, Yates T, Li Y, Pepper SD, Miller CJ: A comparison of massively parallel nucleotide sequencing with oligonucleotide microarrays for global transcription profiling.BMC Genomics 2010, 11:282.PubMedView Article
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.