Nutrigenomics in honey bees: digital gene expression analysis of pollen's nutritive effects on healthy and varroa-parasitized bees
© Alaux et al; licensee BioMed Central Ltd. 2011
Received: 7 January 2011
Accepted: 10 October 2011
Published: 10 October 2011
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
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
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
Genes differentially transcribed after pollen feeding or varroa parasitism
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, v arroa 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
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.
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
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
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
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
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
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).
- Lopez-Maury L, Marguerat S, Bahler J: Tuning gene expression to changing environments: from rapid responses to evolutionary adaptation. Nat Rev Genet. 2008, 9 (8): 583-593.View ArticlePubMedGoogle Scholar
- Grewal SS: Insulin/TOR signaling in growth and homeostasis: a view from the fly world. Int J Biochem Cell Biol. 2009, 41 (5): 1006-1010. 10.1016/j.biocel.2008.10.010.View ArticlePubMedGoogle Scholar
- Field CJ, Johnson IR, Schley PD: Nutrients and their role in host resistance to infection. J Leukoc Biol. 2002, 71: 16-32.PubMedGoogle Scholar
- Li P, Yin YL, Li D, Kim SW, Wu G: Amino acids and immune function. Brit J Nutr. 2007, 98: 237-252. 10.1017/S000711450769936X.View ArticlePubMedGoogle Scholar
- Kaput J, Rodriguez RL: Nutritional genomics: the next frontier in the postgenomic era. Physiol Genomics. 2004, 16 (2): 166-177.View ArticlePubMedGoogle Scholar
- Müller M, Kersten S: Nutrigenomics: goals and strategies. Nat Rev Genet. 2003, 4 (4): 315-322. 10.1038/nrg1047.View ArticlePubMedGoogle Scholar
- Brodschneider R, Crailsheim K: Nutrition and health in honey bees. Apidologie. 2010, 41: 278-294. 10.1051/apido/2010012.View ArticleGoogle Scholar
- Honey Bee Genome Sequencing Consortium: Insights into social insects from the genome of the honeybee Apis mellifera. Nature. 2006, 443: 931-949. 10.1038/nature05260.View ArticleGoogle Scholar
- Roulston TH, Buchmann SL: A phylogenetic reconsideration of the pollen starch-pollination correlation. Evol Ecol Res. 2000, 2: 627-643.Google Scholar
- Roulston TH, Cane JH: Pollen nutritional content and digestibility for animals. Plant Syst Evol. 2000, 222: 187-209. 10.1007/BF00984102.View ArticleGoogle Scholar
- Stanley RG, Linskens HF: Pollen: Biology, biochemistry, management. 1974, Heidelberg, Germany: SpringerView ArticleGoogle Scholar
- Haydak MH: Honey bee nutrition. Annu Rev Entomol. 1970, 15: 143-156. 10.1146/annurev.en.15.010170.001043.View ArticleGoogle Scholar
- Wahl O, Ulm K: Influence of pollen feeding and physiological condition on pesticide sensitivity of the honey bee Apis mellifera carnica. Oecologia. 1983, 59: 106-128. 10.1007/BF00388082.View ArticlePubMedGoogle Scholar
- Rinderer TE, Elliott KD: Worker honey bee response to infection with Nosema apis. J Econ Entomol. 1977, 70: 431-433.View ArticleGoogle Scholar
- Rinderer TE, Rothenbuhler WC, Gochnauer TA: The influence of pollen on the susceptibility of honey-bee larvae to Bacillus larvae. J Invertebr Pathol. 1974, 23: 347-350. 10.1016/0022-2011(74)90100-1.View ArticlePubMedGoogle Scholar
- 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 mellifera L.). J Insect Physiol. 2010, 56: 1184-1191. 10.1016/j.jinsphys.2010.03.017.View ArticlePubMedGoogle Scholar
- Alaux C, Ducloz F, Crauser D, Le Conte Y: Diet effects on honeybee immunocompetence. Biol Lett. 2010, 6: 562-565. 10.1098/rsbl.2009.0986.View ArticlePubMedPubMed CentralGoogle Scholar
- Sammataro D, Gerson U, Needham G: Parasitic mites of honey bees: life history, implications, and impact. Annu Rev Entomol. 2000, 45: 519-548. 10.1146/annurev.ento.45.1.519.View ArticlePubMedGoogle Scholar
- Bowen-Walker PL, Martin SJ, Gunn A: The transmission of deformed wing virus between honeybees (Apis mellifera L.) by the ectoparasitic mite Varroa jacobsoni Oud. J Invertebr Pathol. 1999, 73: 101-106. 10.1006/jipa.1998.4807.View ArticlePubMedGoogle Scholar
- 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. 10.1046/j.1570-7458.2001.00905.x.View ArticleGoogle Scholar
- 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. 10.1073/pnas.0501860102.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Y, Liu X, Zhang W, Han R: Differential gene expression of the honey bees Apis mellifera and A. cerana induced by Varroa destructor infection. J Insect Physiol. 2010, 56 (9): 1207-1218. 10.1016/j.jinsphys.2010.03.019.View ArticlePubMedGoogle Scholar
- 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 bee Apis mellifera associated with Varroa destructor infection. BMC Genomics. 2008, 9: 301-10.1186/1471-2164-9-301.View ArticlePubMedPubMed CentralGoogle Scholar
- Anisimov SV: Serial Analysis of Gene Expression (SAGE): 13 years of application in research. Curr Pharm Biotechnol. 2008, 9 (5): 338-350. 10.2174/138920108785915148.View ArticlePubMedGoogle Scholar
- 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. 10.1101/gr.079558.108.View ArticlePubMedPubMed CentralGoogle Scholar
- '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-10.1093/nar/gkn705.View ArticlePubMedPubMed CentralGoogle Scholar
- 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 mite Varroa destructor. J Gen Virol. 2004, 85 (Pt 12): 3747-3755.View ArticlePubMedGoogle Scholar
- 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 mellifera L.), under various colony conditions. J Virol. 2006, 80 (23): 11528-11538. 10.1128/JVI.00754-06.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu F, Li W, Li Z, Zhang S, Chen S, Su S: High-abundance mRNAs in Apis mellifera: comparison between nurses and foragers. J Insect Physiol. 2011, 57: (274-279). 10.1016/j.jinsphys.2010.11.015.View ArticlePubMedGoogle Scholar
- 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-10.1186/1471-2164-11-400.View ArticlePubMedPubMed CentralGoogle Scholar
- 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-10.1186/1471-2164-11-619.View ArticlePubMedPubMed CentralGoogle Scholar
- 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. 10.1111/j.1601-183X.2009.00480.x.View ArticlePubMedGoogle Scholar
- Engels E: Occurrence and significance of vitellogenins in female castes of social Hymenoptera. Am Zool. 1974, 14: 1229-1237.View ArticleGoogle Scholar
- 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. 10.1016/j.asd.2004.06.006.View ArticlePubMedGoogle Scholar
- 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. 10.1073/pnas.0701909104.View ArticlePubMedPubMed CentralGoogle Scholar
- 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. 10.1073/pnas.0502681103.View ArticlePubMedPubMed CentralGoogle Scholar
- Bitondi MMG, Simões ZLP: The relationship between level of pollen in the diet, vitellogenin and juvenile hormone titres in Africanized Apis mellifera workers. J Apic Res. 1996, 35: 27-36.View ArticleGoogle Scholar
- Amdam GV, Hartfelder K, Norberg K, Hagen A, Omholt SW: Altered physiology in worker honey bees (Hymenoptera: Apidae) infested with the mite Varroa destructor (Acari: Varroidae): A factor in colony loss during overwintering?. J Econ Entomol. 2004, 97: 741-747. 10.1603/0022-0493(2004)097[0741:APIWHB]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Orgad S, Nelson H, Segal D, Nelson N: Metal ions suppress the abnormal taste behavior of the Drosophila mutant malvolio. J Exp Biol. 1998, 201 (Pt 1): 115-120.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Crailsheim K, Schneider LHW, Hrassnigg N, Bühlmann G, Brosch U, Gmeinbauer R, Schöffmann B: Pollen consumption and utilization in worker honeybees (Apis mellifera carnica): dependence on individual age and function. J Insect Physiol. 1992, 38: 409-419. 10.1016/0022-1910(92)90117-V.View ArticleGoogle Scholar
- 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 bees Apis mellifera. Insect Mol Biol. 2006, 15: 645-656. 10.1111/j.1365-2583.2006.00682.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Pernal SF, Currie RW: Pollen quality of fresh and 1-year-old single pollen diets for worker honey bees (Apis mellifera L.). Apidologie. 2000, 31: 387-409. 10.1051/apido:2000130.View ArticleGoogle Scholar
- Garedew A, Schmolz E, Lamprecht I: The energy and nutritional demand of the parasitic life of the mite Varroa destructor. Apidologie. 2004, 35: 419-430. 10.1051/apido:2004032.View ArticleGoogle Scholar
- de Groot AP: Protein and amino acid requirements of the hone y bee (Apis mellifica L.). Physiol Comp Oecol. 1953, 3: 197-285.Google Scholar
- Colombani J, Raisin S, Pantalacci S, Radimerski T, Montagne J, Leopold P: A nutrient sensor mechanism controls Drosophila growth. Cell. 2003, 114 (6): 739-749. 10.1016/S0092-8674(03)00713-X.View ArticlePubMedGoogle Scholar
- 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. 10.1016/S0962-8924(02)00042-9.View ArticlePubMedGoogle Scholar
- Münch D, Amdam GV: The curious case of aging plasticity in honey bees. FEBS Lett. 2010, 584 (12): 2496-2503. 10.1016/j.febslet.2010.04.007.View ArticlePubMedGoogle Scholar
- 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-10.1371/journal.pone.0004899.View ArticlePubMedPubMed CentralGoogle Scholar
- 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. 10.1073/pnas.0800630105.View ArticlePubMedPubMed CentralGoogle Scholar
- 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-10.1371/journal.pgen.1000896.View ArticlePubMedPubMed CentralGoogle Scholar
- Amdam GV: Social context, stress, and plasticity of aging. Aging Cell. 2011, 10 (1): 18-27. 10.1111/j.1474-9726.2010.00647.x.View ArticlePubMedGoogle Scholar
- Hennig KM, Colombani J, Neufeld TP: TOR coordinates bulk and targeted endocytosis in the Drosophila melanogaster fat body to regulate cell growth. J Cell Biol. 2006, 173 (6): 963-974. 10.1083/jcb.200511140.View ArticlePubMedPubMed CentralGoogle Scholar
- Li L, Edgar BA, Grewal SS: Nutritional control of gene expression in Drosophila larvae via TOR, Myc and a novel cis-regulatory element. BMC Cell Biol. 2010, 11: 7-10.1186/1471-2121-11-7.View ArticlePubMedPubMed CentralGoogle Scholar
- Hay N, Sonenberg N: Upstream and downstream of mTOR. Genes Dev. 2004, 18 (16): 1926-1945. 10.1101/gad.1212704.View ArticlePubMedGoogle Scholar
- Grewal SS, Evans JR, Edgar BA: Drosophila TIF-IA is required for ribosome synthesis and cell growth and is regulated by the TOR pathway. J Cell Biol. 2007, 179: 1105-1113. 10.1083/jcb.200709044.View ArticlePubMedPubMed CentralGoogle Scholar
- 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. 10.1016/j.cub.2006.03.084.View ArticlePubMedGoogle Scholar
- Tiefenbock SK, Baltzer C, Egli NA, Frei C: The Drosophila PGC-1 homologue Spargel coordinates mitochondrial activity to insulin signalling. Embo J. 2009, 29 (1): 171-183.View ArticlePubMedPubMed CentralGoogle Scholar
- Baltzer C, Tiefenbock SK, Marti M, Frei C: Nutrition controls mitochondrial biogenesis in the Drosophila adipose tissue through Delg and cyclin D/Cdk4. PLoS ONE. 2009, 4 (9): e6935-10.1371/journal.pone.0006935.View ArticlePubMedPubMed CentralGoogle Scholar
- Schneider P, Drescher W: The influence of Varroa jacobsoni Oud. on weight, development of weight and hypopharyngeal glands, and longevity of Apis mellifera L. Apidologie. 1987, 18: 101-110. 10.1051/apido:19870108.View ArticleGoogle Scholar
- 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 ArticleGoogle Scholar
- 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-10.1186/gb-2004-5-12-r101.View ArticlePubMedPubMed CentralGoogle Scholar
- Kapahi P, Zid BM, Harper T, Koslover D, Sapin V, Benzer S: Regulation of lifespan in Drosophila by modulation of genes in the TOR signaling pathway. Curr Biol. 2004, 14 (10): 885-890. 10.1016/j.cub.2004.03.059.View ArticlePubMedPubMed CentralGoogle Scholar
- 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. 10.1016/j.mad.2005.07.004.View ArticlePubMedGoogle Scholar
- Sun J, Tower J: FLP recombinase-mediated induction of Cu/Zn-superoxide dismutase transgene expression can extend the life span of adult Drosophila melanogaster flies. Mol Cell Biol. 1999, 19 (1): 216-228.View ArticlePubMedPubMed CentralGoogle Scholar
- Orr WC, Mockett RJ, Sohal RS: Overexpression of glutathione reductase extends survival in transgenic Drosophila melanogaster under hyperoxia but not normoxia. Faseb J. 1999, 13 (13): 1733-1742.PubMedGoogle Scholar
- Woodward B: Protein, calories, and immune defenses. Nutr Rev. 1998, 56 (1 Pt 2): S84-92.PubMedGoogle Scholar
- Lee KP, Simpson SJ, Wilson K: Dietary protein-quality influences melanization and immune function in an insect. Funct Ecol. 2008, 22: 1052-1061. 10.1111/j.1365-2435.2008.01459.x.View ArticleGoogle Scholar
- 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. 10.1098/rspb.2005.3385.View ArticleGoogle Scholar
- Lemaitre B, Nicolas E, Michaut L, Reichhart JM, Hoffmann JA: The dorsoventral regulatory gene cassette spatzle/Toll/cactus controls the potent antifungal response in Drosophila adults. Cell. 1996, 86 (6): 973-983. 10.1016/S0092-8674(00)80172-5.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Vizioli J, Salzet M: Antimicrobial peptides from animals: focus on invertebrates. Trends Pharmacol Sci. 2002, 23 (11): 494-496. 10.1016/S0165-6147(02)02105-3.View ArticlePubMedGoogle Scholar
- Fievet J, Tentcheva D, Gauthier L, de Miranda J, Cousserans F, Colin ME, Bergoin M: Localization of deformed wing virus infection in queen and drone Apis mellifera L. Virol J. 2006, 3: 16-10.1186/1743-422X-3-16.View ArticlePubMedPubMed CentralGoogle Scholar
- 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-10.1186/1471-2164-11-282.View ArticlePubMedPubMed CentralGoogle Scholar
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.