Effects of immunostimulation on social behavior, chemical communication and genome-wide gene expression in honey bee workers (Apis mellifera)
© Richard et al.; licensee BioMed Central Ltd. 2012
Received: 24 April 2012
Accepted: 8 October 2012
Published: 16 October 2012
Social insects, such as honey bees, use molecular, physiological and behavioral responses to combat pathogens and parasites. The honey bee genome contains all of the canonical insect immune response pathways, and several studies have demonstrated that pathogens can activate expression of immune effectors. Honey bees also use behavioral responses, termed social immunity, to collectively defend their hives from pathogens and parasites. These responses include hygienic behavior (where workers remove diseased brood) and allo-grooming (where workers remove ectoparasites from nestmates). We have previously demonstrated that immunostimulation causes changes in the cuticular hydrocarbon profiles of workers, which results in altered worker-worker social interactions. Thus, cuticular hydrocarbons may enable workers to identify sick nestmates, and adjust their behavior in response. Here, we test the specificity of behavioral, chemical and genomic responses to immunostimulation by challenging workers with a panel of different immune stimulants (saline, Sephadex beads and Gram-negative bacteria E. coli).
While only bacteria-injected bees elicited altered behavioral responses from healthy nestmates compared to controls, all treatments resulted in significant changes in cuticular hydrocarbon profiles. Immunostimulation caused significant changes in expression of hundreds of genes, the majority of which have not been identified as members of the canonical immune response pathways. Furthermore, several new candidate genes that may play a role in cuticular hydrocarbon biosynthesis were identified. Effects of immune challenge expression of several genes involved in immune response, cuticular hydrocarbon biosynthesis, and the Notch signaling pathway were confirmed using quantitative real-time PCR. Finally, we identified common genes regulated by pathogen challenge in honey bees and other insects.
These results demonstrate that honey bee genomic responses to immunostimulation are substantially broader than the previously identified canonical immune response pathways, and may mediate the behavioral changes associated with social immunity by orchestrating changes in chemical signaling. These studies lay the groundwork for future research into the genomic responses of honey bees to native honey bee parasites and pathogens.
KeywordsHoney bees Immunity Genomics Social behavior Social immunity Chemical ecology Cuticular hydrocarbons Nestmate recognition
Honey bees are an outstanding model system for studying the molecular, physiological and social basis of disease transmission and resistance. Honey bees are plagued by a number of parasites and pathogens (reviewed in), and the social colony environment (with up to 50,000 densely packed worker bees) provides excellent conditions for disease transmission. While the innate immune systems of invertebrates and vertebrates are surprisingly conserved, there can be large differences in the numbers of genes involved in the different molecular arms of the immune response system across species. As is the case with many insects with recently sequenced genomes[4–8], bees have a much smaller number of known canonical immune response genes relative to Drosophila, which is one of the best characterized models of insect immunity. Thus, bees (and other insects) may utilize alternative genetic and physiological mechanisms to respond to infections. Furthermore, honey bees can resist pathogens and parasites by employing sophisticated behavioral defense mechanisms, termed “social immunity” (reviewed in). Here, we examine the effects of a panel of general immune elicitors (injection with saline, Sephadex beads or bacteria) on worker-worker social interactions and chemical communication, and use whole-genome microarrays to characterize global gene expression responses to these immune elicitors.
Honey bee populations have been in decline worldwide, with beekeepers recently reporting massive annual losses (>30%) (reviewed in). This decline is undoubtedly due in part to the multitude of parasites and pathogens that target honey bees, several of which have only recently been identified. Honey bees are host to over 20 viruses[12–14], as well as a number of bacterial and fungal pathogens[15–18] including the gut microspordian parasites Nosema apis and Nosema ceranae. Honey bees are also severely impacted by Varroa mites (Varroa destructor), tracheal mites (Acarapis woodi) and other ectoparasites).
Sequencing of the honey bee genome identified 177 genes associated with the canonical immune response pathways in insects. Innate immune pathways in insects, mostly obtained from studies in Drosophila, consists of both cellular and humoral responses, which can be systemic or local (reviewed in[10, 20–24]). Cellular immune responses involve a number of differentiated hemocytes. Pathogens activate phenoloxidase and associated immune cascades, resulting in phagocytosis, encapsulation, and/or melanization of invading organisms or wounds. Humoral responses include cytotoxic molecules (such as reactive oxygen and nitrogen species), lysozymes, cytokines, and antimicrobial peptides (AMPs). AMPs are primarily produced by the fat bodies in insects. A number of signal transduction pathways are involved in moderating immune responses, including the JAK/STAT (which seems to respond primarily to tissue damage), JNK (which mediates wound repair), Imd (which primarily regulates responses to Gram-negative bacteria), and Spaetzle/Toll (which generally regulates responses to Gram-positive bacteria and fungi) pathways. Activation of one of these signal transduction pathways may lead to up- or down-regulation of the others depending on the host-parasite system. For example, studies in Drosophila have shown that the relationships between immune pathways and effectors can be highly pathogen specific and much still remains to be characterized. Furthermore, canonical pathways derived from Drosophila experiments, may not be fully generalizable to other systems.
Honey bees also employ a number of behavioral mechanisms – termed “social immunity” – to reduce the impacts of parasites and diseases (reviewed in[11, 26]). Grooming can remove ectoparasites, and hygienic behavior can reduce levels of brood diseases and depress Varroa mite populations (mites feed on developing pupae). Workers also coat the interior of their colonies with propolis (derived from plant sap), which has antimicrobial properties and results in lowered levels of immune response genes in individual bees. Furthermore, diseased or otherwise stressed bees undergo accelerated behavioral maturation, moving from in-hive tasks to foraging, which removes them from the brood nest and colony, thereby potentially reducing transmission[28, 29]. Finally, studies in termites suggest that healthy nestmates can become “socially vaccinated” by exposure to infected nestmates, thereby improving their resistance to a pathogen. Future investigations may uncover similar defense mechanisms in honey bees.
Honey bees can also use cuticular hydrocarbons as chemical cues to distinguish between healthy and immuno-stimulated adult nestmates, and respond differently to them. Cuticular hydrocarbons are synthesized in oenocytes, which are embedded in the fat body tissue under the epithelia, and deposited on the cuticular surface. Cuticular hydrocarbon patterns can be modified by genotype, physiological state, and environmental context (reviewed in[33, 34]), including social status in honey bees. Immunostimulation of honey bee workers with lipopolysaccharides derived from bacterial cell walls caused significant changes in cuticular hydrocarbon profiles in worker bees after four hours, and resulted in altered social interactions. Healthy bees were more aggressive towards their nestmates that were coated with extracts of immunostimulated nestmates than to nestmates that were coated with extracts of healthy nestmates, indicating that these changes in chemical profile could alter worker-worker interactions. Other studies have demonstrated that Varroa parasitized pupae and adults have modified cuticular hydrocarbon profiles, while infection with a virus can elicit aggressive reactions from healthy nestmates in a colony.
Here, we examined the behavioral, chemical and molecular responses in honey bee workers to injection with saline, Sephadex beads, and Gram-negative bacteria (freeze-killed E. coli cells), six hours after injection. We determined if these different immune elicitors stimulate unique responses at these three levels; for example, theoretically only E. coli injection should activate the Imd pathway, which could lead to distinct chemical profiles and behavioral responses. We also determined if immunostimulation resulted in significantly altered expression of previously annotated honey bee immune genes, if there was overlap with other studies of the effects of immunochallenge or parasitization in honey bees and Drosophila[38–41], and if there were changes in expression of genes associated with cuticular hydrocarbon biosynthesis pathways.
Pilot study on the effects of immune elicitors on locomotion and mortality
Lethality of our treatments was assessed 8 hours after treatment, in a pilot study using bees from Colony 3. After treatment, 1 out of 10 of the control (handled and CO2 anesthetized) and 1 of the 10 bead- or bacteria-injected individuals died, respectively. None of the saline-injected workers died.
The behavioral effects of treatments (N = 10 for each treatment) were individually assessed on isolated bees placed in circular arenas for 5 minutes as in. Statistical analyses revealed no significant effects of treatments on general activity (total time spent self-grooming: control: 63 ± 23 sec; saline-injected: 56 ± 15 sec; bead-injected: 86 ± 24 sec; bacteria-injected: 83 ± 31 sec; Kruskal-Wallis: H(3)=0.58; p=0.8) or locomotion (number of lines crossed: control: 50 ± 16 sec; saline-injected: 29 ± 8.5 sec; bead-injected: 39 ± 12 sec; bacteria-injected: 43 ± 15.5 sec; Kruskal-Wallis: H(3)=0.507; p = 0.9).
Effects of immune elicitors on social interactions
To monitor social interactions, we used a nestmate recognition assay in which an individual, 6 hours after treatment, was returned to her cage with her original nestmates and social interactions were monitored. 975 bees from Colony 1, separated in 65 petri dishes, were used for this experiment. This assay has been used regularly in bees and ants to examine nestmate recognition[31, 35]. As in the pilot study with isolated individual bees, there was no effect of treatment on locomotion (number of lines crossed: control: 6 ± 4.75; saline-injected: 4 ± 4; bead-injected: 5.5 ± 3; bacteria-injected: 2.5 ± 2.5; Kruskal-Wallis: H(2)=−4.6; p>0.05).
Aggression level index of honey bee behavioral interactions
Allo-grooming(antennal contact, grooming or licking)
Biting with gaster flexion
Furthermore, unlike the pilot study with isolated bees, the frequency of self-grooming behavior was significantly lower in bacteria-injected workers compared to control and saline-injected workers, while self-grooming by bead-injected workers was intermediate (Kruskal-Wallis: H (3, N= 65) = 15,82 p =0.001 followed by post-hoc pairwise comparisons p<0.01, Figure1B). However, bacteria-injected workers were the targets of significantly higher allo-grooming behavior than control workers, while saline- and bead-injected workers were intermediate (Kruskal-Wallis: H (3, N= 65) =7,97 p =0.046 followed by pairwise comparisons p<0.05, Figure1B).
Effects of immune elicitors on cuticular hydrocarbon profiles
Non-polar compounds were extracted from honey bee workers using pentane, and chemical profiles of cuticular extracts were analyzed using gas chromatography (GC). In our experimental conditions, all the chemical compounds we identified by mass spectrometry (MS) analysis were hydrocarbons.
Effects of immunostimulation on cuticular hydrocarbon profiles of bees in Colony 1
Median ± SIQR
Median ± SIQR
Median ± SIQR
Median ± SIQR
Effects of immunostimulation on cuticular hydrocarbon profiles of bees in Colony 2
Median ± SIQR
Median ± SIQR
Median ± SIQR
Median ± SIQR
Hentriacontene Isomere 2
Global gene expression responses to immunostimulation
Effects of individual immune elicitors on gene expression
Analysis of overlap among treatment groups
# of significantly, differentially regulated transcripts
p-value (Fisher's Exact Test)
Saline x Control
p < 0.001
Bead x Control
p < 0.001
Bacteria x Control
p < 0.001
Functional analysis of regulated genes
Gene ontology analysis of the significantly regulated 302 transcripts found in both colonies (207 of which had unique Drosophila orthologs with Flybase annotations and thus were used in the analysis) revealed an overrepresentation of genes involved in immune response, spermatogenesis, wing disc dorsal/ventral pattern formation, tissue development, post-transcriptional regulation of gene expression and protein polymerization (p<0.05, see Additional file1: Table S6 for a listing of the genes in these GO categories). However, none of these categories survived the Benjamini correction. Immune genes corresponding to the major immune pathways were significantly regulated (defensin-1, relish, domeless, cactus, melanization protease 1, death related ced-3/DREDD, PGRP-SC2, kayak, spirit among others). Pale (ple), a tyrosine hydroxylase involved in melanization and wound repair, was also significantly up-regulated by bacteria, bead and saline injection in Colony 1 and bacteria and bead injection in Colony 2. Other genes included those involved in cell growth and proliferation (insulin-like receptor), cytoskeleton structure (basigin, chickadee, twinstar, annexin ix, and isoforms of tubulin), extracellular matrix components (pericardin and laminin A), Notch signaling (apterous, pebbled, groucho), phagocytosis (draper), and cabut, which encodes a transcription factor that is regulated by the JNK cascade.
Of the 22 transcripts (corresponding to 17 unique Flybase genes) significantly regulated by injection with saline, beads, and bacteria relative to controls, two categories were significantly overrepresented: organ morphogenesis, p<0.005, and developmental process, p<0.001 (see Additional file1: Table S7 for a listing of genes in these GO categories). Again, neither of these categories survived the Benjamini correction. Several of the previously discussed genes (domeless, insulin-like receptor, basigin, twinstar, and groucho) are part of this group, as well a serine protease immune response integrator (spirit) which functions in Toll pathway activation. Notably, two genes involved in lipid metabolism were also found in this group, which is of particular interest since fatty acids are the precursors of cuticular hydrocarbon synthesis in insects (reviewed in). Bubblegum (bgm) encodes a very long chain fatty acid CoA ligase, which plays a role in fatty acid metabolism, and was down-regulated in the treatment groups. Lipid-storage droplet 2 (Lsd-2) is involved in lipid storage and accumulation, and was up-regulated in the treatment groups. Lipid storage droplet 1 (Lsd-1), a component of lipid droplets in fat bodies, was found to be significantly down-regulated in both colonies, but not by all treatments.
GO analysis of the 68 transcripts whose expression was specifically altered in bacteria-injected workers (corresponding to 50 unique Flybase genes) yielded only one significantly overrepresented cluster (immune response; p=0.001). These genes included cactus, kayak, draper, defensin-1, relish, PGRP-SC2. Other regulated genes included those that may be involved in metabolism, such as Gr28b, a gustatory receptor that seems to mediate diet-related changes in immune response, sorbital dehydrogenase-2 (Sodh-2), an alcohol dehydrogenase of sugars, and trehalose-6-phospate synthetase 1 (Tps1) which may play a role in protection from hypoxia and anoxic injury. See Additional file1: Tables S8 and S9 for a listing of the bacteria-regulated genes and GO categories.
Comparisons of gene expression patterns with previous studies
We examined overlap in gene expression patterns between our study and previous studies in honey bees and Drosophila. It is important to note that there were large differences in study design (including pathogen challenge, timecourse and tissue) as well as gene expression platform and statistical analyses, which makes direct comparisons challenging. However we sought to identify common regulated genes that could indicate conserved immune response mechanisms across different host organisms and pathogen challenges.
Navajas et al. performed a microarray analysis of control and Varroa mite-parasitized honey bee pupae from two strains (resistant and sensitive to Varroa mite infestation). Expression levels of 32 genes (with 22 matching unique GB identifiers) were significantly altered by Varroa parasitization. There was no significant overlap with the 302 transcripts (matched to 239 GB identifiers) in the current study (Fisher’s Exact Test, p>0.09; Table5). However, two genes, including ple, were regulated in both studies, though directionality was not necessarily conserved between studies.
In a separate study, Alaux and colleagues examined the effects of Varroa parasitization and nutrition on honey bee worker gene expression using an RNA-seq approach. We compared our 302 regulated transcripts (matched to 239 GB identifiers) with the list of genes that were significantly, differentially regulated by mite parasitization in pollen-fed bees and found 117 overlapping genes, including genes involved in immune and metabolic processes. Though directionality was not necessarily conserved between studies, the overlap was greater than expected by chance (Fisher’s Exact Test, p<0.0001; Table5, see Additional file1: Table S11 for a listing of overlapping genes). Bgm, Lsd-2, and domeless were among the genes significantly regulated in both studies. A GO analysis of the overlapping transcripts (corresponding to 102 unique Flybase genes) found one significantly overrepresented category which did not survive Benjamini correction (metabolic process, p<0.04).
Finally, we compared the 302 significantly regulated transcripts in our study (matched to 207 unique Flybase identifiers) with significantly regulated genes in two separate studies examining immune response in fruit flies. De Gregario and colleagues identified 400 significantly regulated genes after septic injury of male Drosophila with bacteria contaminated needles or feeding with fungal spores (corresponding to 497 unique Flybase genes). Roxstrom-Linquist and colleagues identified approximately 390 upregulated genes (corresponding to 464 unique Flybase genes) in Drosophila melanogaster orally infected with protozoa, viruses, bacteria or fungi. Eight genes were found in all three studies (Additional file1: Table S12), including a number of genes with immune related functions: ple, cactus, defensin-1, relish and PGRP-SC2. Serpin 28D (CG7219) was also regulated in all three studies, and is involved in melanization.
Quantitative real-time PCR validation of expression of candidate genes
We have demonstrated that immunostimulation of adult honey bee workers results in altered social interactions only in the case of bacteria-injected bees. However, all the treatment groups have significantly different cuticular hydrocarbon chemical profiles, and all treatments caused large scale changes in gene expression patterns in abdominal fat body and epithelial tissue, which includes oenocytes. It is important to note that these dramatic changes occurred in a relatively short timeframe, only six hours after immunostimulation. Importantly, this study demonstrates that expression of a large number of genes, not just those canonically associated with immune response pathways, are modulated by immunostimulation, as has been demonstrated by several other studies of genome-wide responses to immune challenges (for example,[38–41]). However, comparisons with studies in Drosophila found limited overlap of gene expression patterns, and the handful of commonly regulated genes were key members of canonical immune response pathways.
It is somewhat perplexing that only bacteria-injected bees were subjected to significantly increased allo-grooming and aggression. All treatments caused changes in chemical profiles, and the changes did not appear to be substantially larger in bacteria-injected bees, though there was a trend for lower total quantity of branched alkanes in this group. However, previous studies have suggested that alkenes, rather than alkanes, play a significant role in nestmate recognition. Based on the hierarchical clustering analysis, the gene expression patterns of bacteria-injected bees were clearly distinct from saline- and bead-injected bees. However, overall it did not appear that bacterial injection altered expression of an entirely unique set of genes (out of 302 transcripts, expression of 68 transcripts were significantly regulated only by bacterial injection, relative to controls, and expression of many of these genes changed non-significantly in the other treatments as well). Thus, perhaps bacteria injection did not cause change in expression of a large, unique subset genes, but rather caused more substantial changes in the magnitude of expression levels of some genes, as indicated by hierarchical clustering. Alternatively, it may be that expression of a few key genes and/or chemical profile components in bacteria-treated bees provoked behavioral changes in nestmates.
It remains to be determined if infection with honey bee specific parasites and pathogens causes similar changes in gene expression and behavioral responses, though there is some indication that infection with Nosema microsporidia does stimulate expression changes in a significantly overlapping set of genes (Holt, Aronstein and Grozinger, unpublished data). It also remains to be determined if these behavioral changes are adaptive: it is possible that bees have evolved to “match” the strength of the signal to the virulence of the pathogen, and thus infection with bacteria would elicit a greater change in the behavioral responses of nestmates than simple cuticular wounding. However, colony-level assays would need to be performed in order to determine if these changes in social interactions actually result in differences in the spread of pathogens through the colony and their impacts on infected individuals. For example, the increased grooming behavior we observed towards bacteria-injected individuals could facilitate the spread of pathogens through the colony (as observed in carpenter ants,), slow the spread by causing isolation or removal of infected individuals[56, 57] or reduce the impact on the infected individual[58, 59].
Our studies also demonstrate that immunostimulation elicits complex gene expression changes in the epithelial tissues in honey bee workers. Immunostimulation resulted in significant expression changes of 302 common transcripts in worker bees from the two colonies examined. Only 14 of these genes corresponded to previously annotated immune genes identified from the honey bee genome. Several other biological processes were modified, including cell growth and proliferation, cytoskelatal structure, metabolism and components of the Notch signaling pathway. Cell growth, proliferation, and migration, particularly involving actin-mediated cytoskeletal changes, are required for repairing epithelial wounds. The Notch signaling pathway has not yet been linked to immune response or wound repair, but wound repair uses many of the same developmental pathways that function during dorsal closure in Drosophila development, which does involve Notch signaling. Alternatively, since the insect fat body is involved in regulating many key processes, including metabolism (which was commonly regulated in our study and by Varroa parasitization), these changes may reflect general physiological changes after immunostimulation or stress. Expression changes in three genes of the Notch signaling pathway (apterous, groucho, and pebbled) were confirmed using quantitative real-time PCR. Interestingly, expression was significantly higher in saline-injected bees relative to controls for all three genes, but only expression of apterous was affected in bacteria-injected bees, suggesting that changes in Notch signaling may be modulated temporally or by other signaling pathways.
We found significant expression changes of a number of key immune response genes (for a review of the function of these genes, see[10, 23, 24]). The JAK/STAT pathway is regulated by the Domeless receptor; domeless expression was significantly regulated by immune stimulation in both genotypes of bees in our study. Activation of the IMD pathway requires cleavage of Relish by the caspase DREDD; both Dredd and Relish were significantly regulated in our study. We also observed significant regulation of PGRP-SC2 which suppresses activation the IMD pathway, with bacteria-injected bees showing the highest levels of PGRP-SC2 expression (data not shown). In fruit flies, PGRP-SC1 and PGR-SC2 may function in preventing over-activation of the IMD pathway. The IMD pathway triggers the JNK pathway, which activates the transcriptional regulator AP-1 (which contains Kayak/D-fos), and AP-1 in turn negatively regulates Relish-dependent transcription. We found significant regulation of kayak in our study. Interestingly, we also found significant changes in expression of cabut, which is regulated by the JNK pathway but has not yet been linked to immune function. The Toll pathway operates through transduction factors including spirit, which acts extracellularly and upstream of Spaetzle and the NF-κB protein Dorsal (which is negatively regulated by cactus). We found significant regulation of spirit and cactus in our study. Pale, which plays an important role in melanization and wound repair, and draper, which functions in phagocytosis, were also significantly regulated in our study. As a Gram-negative bacteria, it would be expected that E. coli would primarily stimulate activation of the IMD pathway. However, we observed changes in gene expression of members of the Toll pathway, including cactus, spirit and defensin-1, which exhibits Gram-positive antimicrobial activity and is not regulated by the IMD pathway. Thus, there is likely considerable cross-talk between the pathways.
Our studies also identified several genes which may play a role in altering cuticular hydrocarbon patterns. Cuticular hydrocarbons are synthesized primarily in the oenocytes (reviewed in), which are embedded in the fat body of adult honey bees. Cuticular hydrocarbon biosynthesis involves activation of fatty acids by an acyl-CoA synthetase, chain elongations of fatty-acyl-CoAs to produce very long chain fatty acids, and subsequent conversion to a hydrocarbon, likely by a p450 enzyme. Fatty acids are stored in lipid droplets in the adipoctye cells of the insect fat body. Fatty acids can be released from droplets in the adipocytes and accumulate in the oenocytes; this occurs under starvation conditions in particular. This process is mediated in part by lipid storage droplet-2 (Lsd-2): increased Lsd-2 expression in the fat bodies decreases lipid movement to the oenocytes. We found increased expression of Lsd-2 in immunostimulated bees (see Figure6), suggesting reduced movement of lipids to oenocytes, and perhaps reduced levels of cuticular hydrocarbons. We did observe a decrease in the total relative quantity of all branched alkanes in bacteria-injected workers in both colonies but this difference was not significant.
Bubblegum (bgm) activates long chain fatty acids to form acyl-CoAs (reviewed in), a key step in cuticular hydrocarbon biosynthesis. Bgm was originally described as a Drosophila mutant that resulted in elevated levels of very long chain fatty acids and neurodegeneration. Bgm homologs have been identified in numerous species, including humans and mice, and have been demonstrated to activate long chain (C16) and very long chain (C24) fatty acids. In our study, bgm expression was significantly decreased relative to controls (see Figure6).
Despite large differences in study designs and analysis methods, we found some overlap in gene expression with previous studies examining the effects of Varroa mite parasitization on honey bees[38, 40]. Varroa-responsive genes were significantly associated with basic cellular processes, including cell organization, biogenesis and metabolism. We also found significant changes in functional categories associated with basic cellular processes, such as cell growth, proliferation and cytoskeletal structure. Varroa parasitization also caused changes in expression of pale. As discussed above, this may represent cellular mechanisms for wound-healing. Expression of potential hydrocarbon synthesis genes, namely bgm and Lsd-2, were also regulated by Varroa parasitization. Indeed, Varroa parasitized pupae and adults have modified cuticular hydrocarbon profiles. These differences are likely responsible for stimulating hygienic behavior, in which diseased larvae are removed by adult worker bees, a key component of Varroa resistance.
Comparison with two previous studies[39, 41] examining the effects of immunostimulation on Drosophila global gene expression patterns revealed conserved changes in expression of key immune genes in Drosophila and honey bees (including relish, cactus, defensin-1, spirit, PGRP-SC2, and pale), but otherwise limited overlap in the significantly regulated genes. The lack of similarity could represent species-specific immune responses or simply technical differences – for example, Roxstrom-Lindquist orally infected young male flies with bacteria, fungi and microsporidia, and measured whole-body gene expression changes in only two replicates using Affymetrix microarrays.
Our results suggest that immunostimulation of honey bee workers causes significant changes in gene expression patterns, cuticular hydrocarbon chemical profiles, and, in the case of bacteria-injection, social behavior. While bacteria-injection did not cause expression changes in an entirely unique set of genes or production of unique cuticular hydrocarbons compared to other treatment groups, it does appear that the magnitude of the expression and production changes was greater, and this may have resulted in the altered behavioral responses. As demonstrated by other studies examining genome-wide expression changes associated with immune challenges[38–41], we found that even a short-term immune challenge can result in dramatic changes in gene expression that encompass far more genes than those represented by canonical immune response pathways. For example, we found several genes associated with the Notch signaling pathway, which suggests this pathway may also play a role in mediating immune responses. Furthermore, our study has highlighted potentially new candidate genes for regulating cuticular hydrocarbon synthesis in insects. These chemicals serve many functions in insects, including operating as sex and caste-specific pheromones, and most research on their biosynthesis has focused on desaturase enzymes and p450s. These studies lay the groundwork for future work examining the molecular pathways that mediate immune responses to acute stimulators and chronic infections in bees and other insects, and has paved the way for research into the genes that regulate social immunity.
Honey bee stocks
Honey bee colonies were maintained according to standard practices at the Lake Wheeler Honey Bee Research Facility at North Carolina State University in 2008. Workers for these studies were obtained from three source colonies, headed by queens each instrumentally inseminated with semen from a single, different, male (Glenn Apiaries, Fallbrook, CA). The specific source colonies used for the different experiments were Colonies 1, 2 and 3, and are listed in the experimental details below. Since honey bees are haplodiploid, the coefficient of relatedness among the workers in each colony was 0.75, thereby reducing variation in gene expression responses, chemical profiles, and presumably behavioral responses among the nestmates. Honeycomb frames of late-stage pupae were removed from the colonies and placed in incubators overnight at 33°C, 50% relative humidity (RH).
Honey bee rearing
Newly emerged bees (<12 hours old) were brushed from the frames and placed into modified 10 cm Petri dishes in groups of 15. Dishes were maintained under red light in a temperature and humidity-controlled environmental room (33°C, 50% RH, Phytotron Facility, NCSU). Bees were fed 50% sucrose/water solutions and 45% honey/45% pollen/10% water paste ad libitum. Food was replaced every two days. Bees were also exposed to 0.1 queen bee equivalents of queen mandibular pheromone (QMP) (Pherotech, Vancouver, Canada). Every day, 10 μl of QMP (0.01 queen equivalents/μL in an isopropanol/1% water solution) was placed on a microscope slide and allowed to evaporate before being placed in the cage. This amount of QMP mimics a live queen in assays of worker behavior and physiology[68, 69] and thus should help simulate normal rearing conditions.
When the bees were 10 days old, five individual bees were removed from each cage (leaving 10 in the cage). Three of these were then subjected to one of four treatments. The first group of bees was handled and anesthetized with CO2 for 1 minute (control treatment). A second group of bees was anesthetized and injected with 8 μl of sterile bee saline (130 mM NaCl, 6 mM KCl, 4 mM MgCl2, 5 mM CaCl2, 160 mM sucrose, 25 mM glucose, 10 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid in distilled water, pH 6.7, 500 mOsmol, as in). A third group of bees was anesthetized and injected with 8 μl of a CM-25 Sephadex beads (Sigma-Aldrich, Steinheim, Germany) solution (0.01 g of beads mixed with 500 μl sterile bee saline and vortexed, resulting in approximately 110 beads/8 μl). The fourth group of bees was anesthetized and injected with 8 μL of an E. coli (JM101, Sigma, St Louis, MO) solution suspended in sterile bee saline (3.8*105 cells/bee, following a protocol modified from). The bacteria were grown in LB media (Sigma, St Louis, MO), collected and resuspended in the sterile bee saline at the appropriate concentration (4.75*104 cells/μL), and then stored at −80°C until it was thawed for the injections. Injections were performed into the abdominal cavity through tergites with a nano-injector equipped with a glass needle (Schley Compact Model II Instrument; Honey Bee Insemination Services, Davis, CA, US), using a binocular microscope (Leica MZ6 stereomicroscope, Leica Microsystems, Buffalo Grove, IL). Treated bees were marked with a dot of Testor’s paint on their thorax and maintained individually in 10x10x7 cm3 Plexigas cages in an incubator under red light at 33°C, 50%RH for 6 hours, with 50% sucrose. After 6 hours, two treated bees per cage were collected onto dry ice and stored at −80°C for microarray or chemical analysis, while the third was reintroduced (individually) to her original Petri dish cages for behavioral analysis.
Behavioral analysis was performed on bees from Colony 1 (15 control, 13 saline-, 18 bead-, and 19 bacteria-injected individuals), while chemical analysis was performed on bees from Colony 1 (13 control, 13 saline-, 12 bead- and 16 bacteria-injected individuals) and Colony 2 (12 control, 11 saline-, 13 bead- and 12 bacteria-injected individuals). Microarray analysis was performed on 6 individuals/treatment in Colony 1, and 4 individuals/treatment in Colony 2. Colony 3 was used for pilot studies of the deleterious effects of the treatments (see below).
Behavioral assessment of the deleterious effects of treatment
Assessment of grooming behavior and locomotion
To determine whether the treatments induced any deleterious effects, we conducted a preliminary study in 10-day-old honey bees (N = 10 per treatment, from Colony 3). Immediately following the injection, individual workers were placed into the Petri dish cages, where they were later observed. Grooming behavior and locomotion were assessed after four hours. To assess locomotor activity, a cross was drawn under the circular arena and the number of lines crossed during 5 minutes was used as a locomotion index. All studies were conducted blind to the treatment group. Mortality was assessed 8 hours after treatment and only two workers died.
Assessment of social interactions
Six hours after being exposed to one of the four treatments, individual workers were reintroduced to their original Petri dish cage. Assays were conducted blind to the treatment group. Each cage was used only once. Beginning 30 seconds after the reintroduction, the number and type of social behaviors of the non-treated bees toward the treated bee were recorded through a scanning method, every 10 seconds for 5 minutes. Aggression levels were defined using standard practices, as described in Table1[31, 35]. For each treated bee, the aggression score was calculated by multiplying the number of non-treated bees observed performing a particular behavior during the observation period by the aggression level index (Table1) for that behavior, and summing across all interactions observed.
Hydrocarbon extraction and analysis
Cuticular hydrocarbons were extracted by submerging each frozen bee individually into 1 ml n-pentane with 10 μg of hexadecane (Sigma, St Louis, MO, U.S.A.) as an internal standard, and agitating gently for 10 min. Whole bees were used for all samples. The solvent was reduced to dryness with a gentle stream of N2, and the sample resuspended in 100 μl. Hydrocarbons were quantified by gas chromatography and identified by gas-chromatography–mass spectrometry (GC–MS). We analyzed 2 μl of solution on a HP5890II GC (Agilent, Palo Alto, CA, U.S.A.) equipped with a flame-ionization detector and interfaced with a HP ChemStation (Rev. A.09.03). Splitless injection was made into a 30 m x 0.32 mm x1 μm HP-5 capillary column operated at 150°C for 2 min, increased at 15°C/min to 250°C then at 3°C/min to 300°C and held at this temperature for 15 min. The injector and detector were held at 300°C and 310°C, respectively. Compound identification was achieved by splitless capillary GC–MS using a HP6890 GC and a model 5973A MSD with an electron impact ion source and a HP-5 capillary column (30 m x 0.25 mm ID x 0.25 μm film thickness). This protocol was modified from.
Statistical analyses of behavioral and chemical studies
For behavioral assays, owing to sample size, nonparametric procedures were used to assess significant variations in our data. Non-parametric Kruskal–Wallis ANOVAs on ranks for global comparison were performed, followed by corrected Mann–Whitney U tests for multiple comparisons to compare significant differences in aggression indices between treatment groups. For statistical analyses of the chemical profiles, all compounds present above 5% relative levels were used. To assess the similarities of cuticular hydrocarbon profiles, a stepwise discriminant analysis of the relative proportions of the compounds was employed using Statistica 6.0 (StatSoft, Inc., Tulsa, OK, U.S.A.) as in. Prior to analysis, each peak area was standardized according to. Multivariate analysis of variance is known to be resistant to non-normality in the data because of skewness. However, before conducting the analyses, we tested for normality and equality of variance over all groups, both of which were non-significant.
Individual bees were thawed and dissected under cold RNAlater (Qiagen, Valencia, CA). Abdomens were eviscerated and the cuticles, with the associated fat bodies and oenocytes, were collected. RNA was extracted using an RNeasy kit (Qiagen). RNA was quantified using a Nanodrop 1000 spectrophotometer (Thermo Fischer Scientific, Wilmington, DE) and sample integrity and quality was monitored using agarose gel electrophoresis to confirm the presence of ribosomal RNA bands. 500 ng of RNA/individual was amplified using the Ambion MessageAmp II aRNA Amplification kit (AM1751, Life Technologies, Grand Island, NY). 5 μg of amplified RNA from each sample were labeled independently with Cy3 and Cy5 dyes using the ULS aRNA fluorescent labeling kit (EA-006, Kreatech, Amesterdam, Netherlands). Samples were hybridized to microarrays (two samples/array) in a loop design with dye swaps incorporated, using 24 microarrays for Colony 1 (6 biological replicates, 2 technical replicates per sample) and 16 arrays for Colony 2 (4 biological replicates, 2 technical replicates per sample). Whole genome microarrays (<http://www.ebi.ac.uk/arrayexpress/arrays/A-MEXP-755>) were purchased from the W.M. Keck Center for Functional Genomics at the University of Illinois, Urbana-Champaign. Arrays were scanned using the Axon Genepix 4000B scanner (Molecular Devices, Sunnyvale, CA) using GENEPIX software (Agilent Technologies, Santa Clara, CA).
Analysis of gene expression
where Y represents the residual from the previous model. Treatment, spot and dye are fixed effects and array is a random effect. p-values were corrected for multiple testing using a false discovery rate < 0.01 (proc MULTTEST, SAS). Hierarchical clustering, using the ward method was performed in JMP 9.0.2 (SAS, Cary, NC). Gene ontology analysis was performed using DAVID version 6.7[73, 74] with a cutoff of p <0.05. For all gene ontology (GO) analyses, array transcripts were matched to their Drosophila orthologs in Flybase (http://flybase.org/). All of the array transcripts with Drosophila orthologs were used as a background list. Analysis of overlap between the significantly regulated transcripts in our study or between significantly regulated genes in our study and other gene expression studies was performed using Fisher’s Exact Tests (Dr. Oyvind Langsrud, Statistics Norway, <http://www.langsrud.com/fisher.htm>). Common, differentially expressed transcripts or genes between studies were identified with Venny. For overlap comparisons between our study and other studies examining gene expression in honey bees ([9, 38, 40]), we converted microarray transcript identifiers (AM numbers) to GB identifiers from the honey bee genome annotation (see BeeBasehttp://hymenopteragenome.org/beebase/). For comparisons between our study and other studies examining gene expression in fruit flies ([39, 41]), we converted all genes to Flybase identifiers to find common overlap. The array data has been deposited on the ArrayExpress website (E-MEXP-3708).
Validation of candidate gene expression patterns using quantitative real-time PCR
Primer sequences used to quantify expression of candidate genes from a third biological replicate using quantitative real-time PCR
Accession number/ reference
Lipid Storage Droplet 2
(Huising & Flik, 2005)
False Discovery Rate
- Grh :
Immune deficiency pathway
- Lsd :
Lipid storage droplet
Janus protein kinase (JAK)/Signal transducers and activators of transcription (STAT)
Jun N-terminal Kinase
- K W:
- LB media:
Luria Broth media
Queen Mandibular Pheromone
Trehalose-6-phospate synthetase 1.
We would like to thank Davis Murphy for assistance with the cage experiments, Joe Flowers for expert beekeeping assistance, Coby Schal for use of his GC instruments, Mike Roe for the use of his GC-MS, Sarah Kocher for assistance with statistical analysis of the data, Elina Lastro Niño and Tracey Baumgarten for valuable assistance with the qRT-PCR validation, and members of the Grozinger lab for helpful discussions and critical reading of the manuscript. This research was funded by a USDA-NIFA-AFRI grant to CMG and FJR and an NSF Predoctoral Fellowship to HLH.
- Vanengelsdorp D, Meixner MD: A historical review of managed honey bee populations in Europe and the United States and the factors that may affect them. J Invertebr Pathol. 2010, 103 (Suppl 1): S80-S95.View ArticlePubMed
- Winston ML: The Biology of the Honey Bee. 1987, Cambridge: Harvard University Press
- Brown P: Cinderella goes to the ball. Nature. 2001, 410 (6832): 1018-1020. 10.1038/35074279.View ArticlePubMed
- Bonasio R, Zhang G, Ye C, Mutti NS, Fang X, Qin N, Donahue G, Yang P, Li Q, Li C, et al: Genomic comparison of the ants Camponotus floridanus and Harpegnathos saltator. Science. 2010, 329 (5995): 1068-1071. 10.1126/science.1192428.PubMed CentralView ArticlePubMed
- Gerardo NM, Altincicek B, Anselme C, Atamian H, Barribeau SM, de Vos M, Duncan EJ, Evans JD, Gabaldon T, Ghanim M, et al: Immunity and other defenses in pea aphids, Acyrthosiphon pisum. Genome Biol. 2010, 11 (2): R21-10.1186/gb-2010-11-2-r21.PubMed CentralView ArticlePubMed
- Smith CD, Zimin A, Holt C, Abouheif E, Benton R, Cash E, Croset V, Currie CR, Elhaik E, Elsik CG, et al: Draft genome of the globally widespread and invasive Argentine ant (Linepithema humile). Proc Natl Acad Sci USA. 2011, 108 (14): 5673-5678. 10.1073/pnas.1008617108.PubMed CentralView ArticlePubMed
- Smith CR, Smith CD, Robertson HM, Helmkampf M, Zimin A, Yandell M, Holt C, Hu H, Abouheif E, Benton R, et al: Draft genome of the red harvester ant Pogonomyrmex barbatus. Proc Natl Acad Sci USA. 2011, 108 (14): 5667-5672. 10.1073/pnas.1007901108.PubMed CentralView ArticlePubMed
- Werren JH, Richards S, Desjardins CA, Niehuis O, Gadau J, Colbourne JK, Beukeboom LW, Desplan C, Elsik CG, Grimmelikhuijzen CJ, et al: Functional and evolutionary insights from the genomes of three parasitoid Nasonia species. Science. 2010, 327 (5963): 343-348. 10.1126/science.1178028.View ArticlePubMed
- 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 (5): 645-656. 10.1111/j.1365-2583.2006.00682.x.PubMed CentralView ArticlePubMed
- Lemaitre B, Hoffmann J: The host defense of Drosophila melanogaster. Annu Rev Immunol. 2007, 25: 697-743. 10.1146/annurev.immunol.25.022106.141615.View ArticlePubMed
- Wilson-Rich N, Spivak M, Fefferman NH, Starks PT: Genetic, individual, and group facilitation of disease resistance in insect societies. Annu Rev Entomol. 2009, 54: 405-423. 10.1146/annurev.ento.53.103106.093301.View ArticlePubMed
- Chen YP, Siede R: Honey Bee Viruses. Advances in Virus Research. Edited by: Maramorosh K, Shatkin AJ, Murphy FA. 2007, San Diego, California: Elsevier Academic Press, 34-80.
- Maori E, Lavi S, Mozes-Koch R, Gantman Y, Peretz Y, Edelbaum O, Tanne E, Sela I: Isolation and characterization of Israeli acute paralysis virus, a dicistrovirus affecting honeybees in Israel: evidence for diversity due to intra- and inter-species recombination. J Gen Virol. 2007, 88 (Pt 12): 3428-3438.View ArticlePubMed
- Runckel C, Flenniken ML, Engel JC, Ruby JG, Ganem D, Andino R, DeRisi JL: Temporal Analysis of the Honey Bee Microbiome Reveals Four Novel Viruses and Seasonal Prevalence of Known Viruses, Nosema, and Crithidia. PLoS One. 2011, 6 (6): e20656-10.1371/journal.pone.0020656.PubMed CentralView ArticlePubMed
- Aronstein KA, Murray KD: Chalkbrood disease in honey bees. J Invertebr Pathol. 2010, 103 (Suppl 1): S20-S29.View ArticlePubMed
- Forsgren E: European foulbrood in honey bees. J Invertebr Pathol. 2010, 103 (Suppl 1): S5-S9.View ArticlePubMed
- Genersch E: Paenibacillus larvae and American Foulbrood - long since known and still surprising. Journal fur Verbrauchershutz und Lebensmittelsicherheit. 2008, 3: 429-434. 10.1007/s00003-008-0379-8.View Article
- Chen YP, Huang ZY: Nosema ceranae, a newly identified pathogen of Apis mellifera in the USA and Asia. Apidologie. 2010, 41 (3): 364-374. 10.1051/apido/2010021.View Article
- 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 ArticlePubMed
- Broderick NA, Welchman DP, Lemaitre B: Recognition and response to microbial infection in Drosophila. Insect Infection and Immunity: Evolution, Ecology and Mechanisms. Edited by: Rolff J, Reynolds S. 2009, Oxford: Oxford University Press, 13-33.View Article
- Rolff J, Reynolds SE: Introducing insect infection and immunity. Insect Infection and Immunity: Evolution, Ecology and Mechanisms. Edited by: Rolff J, Reynolds S. 2009, Oxford:Oxford University Press, 1-12.View Article
- Schmid-Hempel P: Evolutionary ecology of insect immune defenses. Annu Rev Entomol. 2005, 50: 529-551. 10.1146/annurev.ento.50.071803.130420.View ArticlePubMed
- Ganesan S, Aggarwal K, Paquette N, Silverman N: NF-kappa B/Rel Proteins and the Humoral Immune Responses of Drosophila melanogaster. Nf-Kb in Health and Disease. 2011, 349: 25-60.View Article
- Martin P, Parkhurst SM: Parallels between tissue repair and embryo morphogenesis. Development. 2004, 131: 3021-3034. 10.1242/dev.01253.View ArticlePubMed
- Dionne MS, Schneider DS: Models of infectious diseases in the fruit fly Drosophila melanogaster. Dis Model Mech. 2008, 1 (1): 43-49. 10.1242/dmm.000307.PubMed CentralView ArticlePubMed
- Cremer S, Armitage SA, Schmid-Hempel P: Social immunity. Curr Biol. 2007, 17 (16): R693-R702. 10.1016/j.cub.2007.06.008.View ArticlePubMed
- Simone M, Evans JD, Spivak M: Resin collection and social immunity in honey bees. Evolution. 2009, 63 (11): 3016-3022. 10.1111/j.1558-5646.2009.00772.x.View ArticlePubMed
- Rueppell O, Hayworth MK, Ross NP: Altruistic self-removal of health-compromised honey bee workers from their hive. J Evol Biol. 2010, 23 (7): 1538-1546. 10.1111/j.1420-9101.2010.02022.x.View ArticlePubMed
- Hassanein MH: The Influence of Infection with Nosema-Apis on the Activities and Longevity of the Worker Honeybee. Ann Appl Biol. 1953, 40 (2): 418-423. 10.1111/j.1744-7348.1953.tb01093.x.View Article
- Traniello JF, Rosengaus RB, Savoie K: The development of immunity in a social insect: evidence for the group facilitation of disease resistance. Proc Natl Acad Sci USA. 2002, 99 (10): 6838-6842. 10.1073/pnas.102176599.PubMed CentralView ArticlePubMed
- Richard FJ, Aubert A, Grozinger CM: Modulation of social interactions by immune stimulation in honey bee, Apis mellifera, workers. BMC Biol. 2008, 6: 50-10.1186/1741-7007-6-50.PubMed CentralView ArticlePubMed
- Fan Y, Zurek L, Dykstra MJ, Schal C: Hydrocarbon synthesis by enzymatically dissociated oenocytes of the abdominal integument of the German Cockroach, Blattella germanica. Naturwissenschaften. 2003, 90 (3): 121-126.PubMed
- Howard RW, Blomquist GJ: Ecological, behavioral, and biochemical aspects of insect hydrocarbons. Annu Rev Entomol. 2005, 50: 371-393. 10.1146/annurev.ento.50.071803.130359.View ArticlePubMed
- Blomquist GJ, Bagneres AG: Insect hydrocarbons: biology, biochemistry and chemical ecology. 2010, Cambridge: Cambridge University PressView Article
- Fan YL, Richard FJ, Rouf N, Grozinger CM: Effects of queen mandibular pheromone on nestmate recognition in worker honeybees, Apis mellifera. Anim Behav. 2010, 79 (3): 649-656. 10.1016/j.anbehav.2009.12.013.View Article
- Salvy M, Martin C, Bagneres AG, Provost E, Roux M, Le Conte Y, Clement JL: Modifications of the cuticular hydrocarbon profile of Apis mellifera worker bees in the presence of the ectoparasitic mite Varroa jacobsoni in brood cells. Parasitology. 2001, 122 (Pt 2): 145-159.PubMed
- Kulincevic JM, Stairs GR, Rothenbuhler WC: A disease of the honey bee causing behavioral changes and mortality. J Invertebr Pathol. 1969, 14 (1): 13-17. 10.1016/0022-2011(69)90004-4.View ArticlePubMed
- Alaux C, Dantec C, Parrinello H, Le Conte Y: Nutrigenomics in honey bees: digital gene expression analysis of pollen's nutritive effects on healthy and varroa-parasitized bees. BMC Genomics. 2011, 12: 496-10.1186/1471-2164-12-496.PubMed CentralView ArticlePubMed
- De Gregorio E, Spellman PT, Rubin GM, Lemaitre B: Genome-wide analysis of the Drosophila immune response by using oligonucleotide microarrays. Proc Natl Acad Sci USA. 2001, 98 (22): 12590-12595. 10.1073/pnas.221458698.PubMed CentralView ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- Roxstrom-Lindquist K, Terenius O, Faye I: Parasite-specific immune response in adult Drosophila melanogaster: a genomic study. EMBO Rep. 2004, 5 (2): 207-212. 10.1038/sj.embor.7400073.PubMed CentralView ArticlePubMed
- Decanini LI, Collins AM, Evans JD: Variation and heritability in immune gene expression by diseased honeybees. J Hered. 2007, 98 (3): 195-201. 10.1093/jhered/esm008.View ArticlePubMed
- Evans JD: Transcriptional immune responses by honey bee larvae during invasion by the bacterial pathogen, Paenibacillus larvae. J Invertebr Pathol. 2004, 85 (2): 105-111. 10.1016/j.jip.2004.02.004.View ArticlePubMed
- Wittkopp PJ, Beldade P: Development and evolution of insect pigmentation: genetic mechanisms and the potential consequences of pleiotropy. Semin Cell Dev Biol. 2009, 20 (1): 65-71. 10.1016/j.semcdb.2008.10.002.View ArticlePubMed
- Munoz-Descalzo S, Terol J, Paricio N: Cabut, aC(2)H(2) zinc finger transcription factor, is required during Drosophila dorsal closure downstream of JNK signaling. Dev Biol. 2005, 287 (1): 168-179. 10.1016/j.ydbio.2005.08.048.View ArticlePubMed
- Kambris Z, Brun S, Jang IH, Nam HJ, Romeo Y, Takahashi K, Lee WJ, Ueda R, Lemaitre B: Drosophila immunity: a large-scale in vivo RNAi screen identifies five serine proteases required for Toll activation. Curr Biol. 2006, 16 (8): 808-813. 10.1016/j.cub.2006.03.020.View ArticlePubMed
- Blomquist GJ: Biosynthesis of cuticular hydrocarbons. Insect hydrocarbons: biology, biochemistry and chemical ecology. Edited by: Blomquist GJ, Bagneres AG. 2010, Cambridge: Cambridge University Press, 35-52.View Article
- Min KT, Benzer S: Preventing neurodegeneration in the Drosophila mutant bubblegum. Science. 1999, 284 (5422): 1985-1988. 10.1126/science.284.5422.1985.View ArticlePubMed
- Teixeira L, Rabouille C, Rorth P, Ephrussi A, Vanzo NF: Drosophila perilipin/ADRP homologue Lsd2 regulates lipid metabolism. Mech Dev. 2003, 120 (9): 1071-1081. 10.1016/S0925-4773(03)00158-8.View ArticlePubMed
- Beller M, Riedel D, Jansch L, Dieterich G, Wehland J, Jackle H, Kuhnlein RP: Characterization of the Drosophila lipid droplet subproteome. Mol Cell Proteomics. 2006, 5 (6): 1082-1094. 10.1074/mcp.M600011-MCP200.View ArticlePubMed
- Ayres JS, Schneider DS: The Role of Anorexia in Resistance and Tolerance to Infections in Drosophila. PLoS Biol. 2009, 7 (7): e1000150-10.1371/journal.pbio.1000150.PubMed CentralView ArticlePubMed
- Luque T, Hjelmqvist L, Marfany G, Danielsson O, El-Ahmad M, Persson B, Jornvall H, Gonzalez-Duarte R: Sorbitol dehydrogenase of Drosophila - Gene, protein, and expression data show a two-gene system. J Biol Chem. 1998, 273 (51): 34293-34301. 10.1074/jbc.273.51.34293.View ArticlePubMed
- Chen QF, Haddad GG: Role of trehalose phosphate synthase and trehalose during hypoxia: from flies to mammals. J Exp Biol. 2004, 207 (18): 3125-3129. 10.1242/jeb.01133.View ArticlePubMed
- Scherfer C, Tang H, Kambris Z, Lhocine N, Hashimoto C, Lemaitre B: Drosophila Serpin-28D regulates hemolymph phenoloxidase activity and adult pigmentation. Dev Biol. 2008, 323 (2): 189-196. 10.1016/j.ydbio.2008.08.030.View ArticlePubMed
- Dani FR, Jones GR, Corsi S, Beard R, Pradella D, Turillazzi S: Nestmate recognition cues in the honey bee: differential importance of cuticular alkanes and alkenes. Chem Senses. 2005, 30 (6): 477-489. 10.1093/chemse/bji040.View ArticlePubMed
- Hamilton C, Lejeune BT, Rosengaus RB: Trophallaxis and prophylaxis: social immunity in the carpenter ant Camponotus pennsylvanicus. Biol Lett. 2011, 7 (1): 89-92. 10.1098/rsbl.2010.0466.PubMed CentralView ArticlePubMed
- Aubert A, Richard FJ: Social management of LPS-induced inflammation in Formica polyctena ants. Brain Behav Immun. 2008, 22 (6): 833-837. 10.1016/j.bbi.2008.01.010.View ArticlePubMed
- Hughes WOH, Eilenberg J, Boomsma JJ: Trade-offs in group living: transmission and disease resistance in leaf-cutting ants. Proc R Soc Lond B. 2002, 269: 1811-1819. 10.1098/rspb.2002.2113.View Article
- Traniello JFA, Rosengaus RB, Savoie K: The development of immunity in a social evidence for the group facilitation of disease resistance. Proc Natl Acad Sci USA. 2002, 99: 6838-6842. 10.1073/pnas.102176599.PubMed CentralView ArticlePubMed
- Zecchini V, Brennan K, Martinez-Arias A: An activity of Notch regulates JNK signalling and affects dorsal closure in Drosophila. Curr Biol. 1999, 9 (9): 460-469. 10.1016/S0960-9822(99)80211-5.View ArticlePubMed
- Arrese EL, Soulages JL: Insect Fat Body: Energy, Metabolism, and Regulation. Annu Rev Entomol. 2010, 55: 207-225. 10.1146/annurev-ento-112408-085356.PubMed CentralView ArticlePubMed
- Bischoff V, Vignal C, Duvic B, Boneca IG, Hoffmann JA, Royet J: Downregulation of the Drosophila immune response by peptidoglycan-recognition proteins SC1 and SC2. PLoS Pathog. 2006, 2 (2): e14-10.1371/journal.ppat.0020014.PubMed CentralView ArticlePubMed
- Schluns H, Crozier RH: Relish regulates expression of antimicrobial peptide genes in the honeybee, Apis mellifera, shown by RNA interference. Insect Mol Biol. 2007, 16 (6): 753-759. 10.1111/j.1365-2583.2007.00768.x.View ArticlePubMed
- Dade HA: Anatomy and Dissection of the Honeybee. 1994, Cardiff: International Bee Research Association
- Gutierrez E, Wiggins D, Fielding B, Gould AP: Specialized hepatocyte-like cells regulate Drosophila lipid metabolism. Nature. 2007, 445 (7125): 275-280. 10.1038/nature05382.View ArticlePubMed
- Mashek DG, Li LO, Coleman RA: Long-chain acyl-CoA synthetases and fatty acid channeling. Future Lipidology. 2007, 2 (4): 465-476. 10.2217/174608220.127.116.115.PubMed CentralView ArticlePubMed
- Pei ZT, Oey NA, Zuidervaart MM, Jia ZZ, Li YY, Steinberg SJ, Smith KD, Watkins PA: The acyl-CoA synthetase "bubblegum" (lipidosin) - Further characterization and role in neuronal fatty acid beta-oxidation. J Biol Chem. 2003, 278 (47): 47070-47078. 10.1074/jbc.M310075200.View ArticlePubMed
- Grozinger CM, Sharabash NM, Whitfield CW, Robinson GE: Pheromone-mediated gene expression in the honey bee brain. Proc Natl Acad Sci USA. 2003, 100 (Suppl 2): 14519-14525.PubMed CentralView ArticlePubMed
- Hoover SE, Keeling CI, Winston ML, Slessor KN: The effect of queen pheromones on worker honey bee ovary development. Naturwissenschaften. 2003, 90 (10): 477-480. 10.1007/s00114-003-0462-z.View ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- Richard FJ, Poulsen M, Drijfhout F, Jones G, Boomsma JJ: Specificity in chemical profiles of workers, brood and mutualistic fungi in Atta, Acromyrmex, and Sericomyrmex fungus-growing ants. J Chem Ecol. 2007, 33 (12): 2281-2292. 10.1007/s10886-007-9385-z.View ArticlePubMed
- Reyment RA: Compositional data analysis. Terra Nova. 1989, 1 (1): 29-34. 10.1111/j.1365-3121.1989.tb00322.x.View Article
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4 (5): P3-10.1186/gb-2003-4-5-p3.View ArticlePubMed
- da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57.View ArticlePubMed
- VENNY: An interactive tool for comparing lists with Venn Diagrams.http://bioinfogp.cnb.csic.es/tools/venny/index.html,
- Fussnecker BL, McKenzie AM, Grozinger CM: cGMP modulates responses to queen mandibular pheromone in worker honey bees. J Comp Physiol A Neuroethol Sens Neural Behav Physiol. 2011, 197 (9): 939-948. 10.1007/s00359-011-0654-5.PubMed CentralView ArticlePubMed
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.