Skip to main content
  • Research article
  • Open access
  • Published:

Sublethal effects of acaricides and Nosema ceranae infection on immune related gene expression in honeybees


Nosema ceranae is an obligate intracellular parasite and the etiologic agent of Nosemosis that affects honeybees. Beside the stress caused by this pathogen, honeybee colonies are exposed to pesticides under beekeeper intervention, such as acaricides to control Varroa mites. These compounds can accumulate at high concentrations in apicultural matrices. In this work, the effects of parasitosis/acaricide on genes involved in honeybee immunity and survival were evaluated. Nurse bees were infected with N. ceranae and/or were chronically treated with sublethal doses of coumaphos or tau-fluvalinate, the two most abundant pesticides recorded in productive hives. Our results demonstrate the following: (1) honeybee survival was not affected by any of the treatments; (2) parasite development was not altered by acaricide treatments; (3) coumaphos exposure decreased lysozyme expression; (4) N. ceranae reduced levels of vitellogenin transcripts independently of the presence of acaricides. However, combined effects among stressors on imagoes were not recorded. Sublethal doses of acaricides and their interaction with other ubiquitous parasites in colonies, extending the experimental time, are of particular interest in further research work.


Pesticide exposure and pathogen infection are recognised as potential stressors affecting honeybee populations and consequently the profitability of beekeeping, such as parasitic mite syndrome or varroosis [1, 2], several viruses [3, 4], xenobiotics [5, 6] and microsporidiosis [7].

Nosema ceranae [810] is the predominant microsporidian pathogen, originally isolated from Apis cerana [11] but recently it has been reported as a natural parasite of A. mellifera over the world [12]. This obligate intracellular microparasite spreads among hosts via spores and is the etiologic agent of Nosemosis.

With respect to chemical stressors, honeybee colonies are exposed to pesticides used in agriculture or within bee hives by beekeeper intervention. In order to control Varroa destructor, the causative agent of varroosis, beekeepers apply synthetic acaricides, with the most commonly employed being the neurotoxic compounds coumaphos (organophosphate) and tau-fluvalinate (pyrethroid) (reviewed in [13]). These acaricides are applied directly to bee hives and accumulate in apicultural matrices, particularly beeswax [14], being detected in nearly all hives tested in several countries [1519].

It has been suggested that tau-fluvalinate increases honeybee susceptibility to deformed wing virus infection [20], and that coumaphos can act as a systemic agent inside the colonies [21] and alters some metabolic responses [22]. Nevertheless, these molecules did not alter immunity at the gene expression level in bees subject to acute exposure in laboratory conditions [23].

Insects lack an acquired immune system but have a well-developed innate response [24]. In honeybees, humoral immunity involves antimicrobial peptide (AMP) synthesis [25]. Four AMP, abaecin, apidaecin, defensin and hymenoptaecin have been identified in A. mellifera when honeybees are challenged by bacteria [26, 27]. On the contrary, cellular defence involves phagocytosis, encapsulation and nodulation mediated mostly by hemocytes [28]. These latter processes are related to melanisation, a mechanism catalysed by phenoloxidase [29] and glucose dehydrogenase (GLD) [30]. Another component of worker bee defence is vitellogenin, synthesized by the fat body and released into the hemolymph [31, 32]. This yolk protein is involved in many physiological processes, for example in the regulation of honeybee longevity, immunity, brood food synthesis or acting as a stress response mediator [33]. Antúnez et al. [34] showed evidence that honeybee colonies with different degrees of N. ceranae infection could be associated to a differential expression of vitellogenin.

Currently, the requirement of high doses of acaricides due to V. destructor resistance added to the potential of honeybee colonies to accumulate acaricide residues and their metabolites has become detrimental to honeybee health and productivity. For this reason, the evaluation of sublethal effects of acaricides is essential to gain information on the long term effects of these widely used compounds in beekeeping.

The objective of this study was to determine whether exposure of adult bees to chronic sublethal acaricide (tau-fluvalinate and coumaphos) and/or infection with N. ceranae generates significant effects on individual immune components and honeybee survival.

Materials and methods

Colony health conditions and preparation of acaricide-free combs

Local hybrid honeybees (A. mellifera mellifera/A. mellifera ligustica) were obtained from colonies located in the experimental apiary J. J. Nágera coastal station (38°10′06″S, 57°38′10″W), Mar del Plata, Argentina.

Mite infestation levels in the apiary were monitored using the natural mite fall method [13]. Also, phoretic Varroa mite infestation rates were determined according to Fries et al. [35].

In order to evaluate the presence of Nosema spp. spores in the apiary, 60 forager honeybees returning from flight from each colony were macerated and examined using light microscopy [36]. Three colonies were selected to obtain brood combs based on the lowest prevalence recorded from Varroa mites (less than 1%) and the lowest abundance from N. ceranae (an average of 5 × 105 Nosema spores/bee). Specifically, neither of the colonies used to obtain the imagoes presented any visible clinical symptoms of other diseases (i.e., American foulbrood or chalkbrood).

No chemical treatment was applied to the colonies 5 months before and during the experiment.

Four months prior to the beginning of the assays, in order to avoid the presence of long lasting acaricide residues in wax or honey [14], a plastic foundation was used allowing bees to draw out the foundation since new commercial beeswax foundation are paradoxically contaminated with residuals [37]. Plastic foundations were covered with a thin layer of virgin wax to get a better acceptance of these supports. Then they were placed in the selected colonies.

Nosema ceranae spores

Foraging bees were collected from a naturally infected colony of the experimental apiary mentioned above. The spore number was determined using a hemocytometer under light microscopy. The presence of N. ceranae spores and the absence of N. apis spores was confirmed by PCR [38].

In order to avoid microbial contaminants and host tissue, fresh spore suspensions were purified following a triangulation method modified by Fries et al. [39]. The obtained suspension was immediately used for experimental infection.


Technical grade chemicals were used for all trials. The acaricides coumaphos and tau-fluvalinate were obtained from Sigma-Aldrich®.

Experimental infection with N. ceranae spores

Newly emerged bees, which are free from N. ceranae infection [9, 40], from three colonies were mixed together and were placed in wooden cages with a plastic mesh (11 × 9 × 6 cm3) in groups of 150. They were also provisioned with bee bread (stored at −20 °C for 3 months in order to eliminate viable spores) manually collected from plastic combs from honeybee colonies, water and sucrose syrup ad libitum.

Individual infection was achieved on day 3 after emergence according to Porrini et al. [41]. Honeybees were starved 5 h and then more than 600 individuals received 10 μL of sucrose–water solution (2:1 w:v) with 1.5 × 105 Nosema spores. The solution was continuously vortexed to ensure a uniform suspension. Approximately 650 control honeybees were treated in the same manner using sucrose-water solution (2:1 w:v) containing no Nosema spores. Those individuals that did not consume the total amount of solution were discarded from the assay. Then, they were caged in transparent and ventilated plastic flasks of 900 cm3, these hoarding cages have inputs for feeding devices and removable sides. Four replicates of 50 imagoes per treatment were made and maintained in an incubator (28 °C; 30% HR). Honeybees were fed ad libitum with bee bread, water and sucrose syrup. A chronology of the experiment is shown in Table 1.

Table 1 Chronological description of experimental procedure

Selection of sublethal dose of acaricides

The experiment was designed in order to simulate realistic conditions of exposure to acaricides. Concentrations were based on previous data on the presence of these pesticides in pollen and honey. In the case of tau-fluvalinate, a concentration of 750 ppb was detected in honey [42] and 487.2 ppb in pollen [43], the selected dose was approximately the average between those values (666 ppb). In the case of coumaphos, 2020 ppb were found in honey samples [17] and 5800 ppb in pollen [18], so the selected dose was approximately the average (3333 ppb).

Chronic exposure bioassays

Four-day-old worker bees were used for all bioassays (Table 1). Acaricidal compounds were dissolved in ethanol 96%. The final concentration of ethanol in sugar syrup was 2%; at this percentage it has been demonstrated that behaviour components are not affected in acute exposures [44]. Additionally, control bees received the same concentration of ethanol. Honeybees were treated per os with solutions of coumaphos (3333 ppb) or tau-fluvalinate (666 ppb) added to sugar syrup. Bees were supplied with one of the following diets: sucrose syrup + tau-fluvalinate, sucrose syrup + coumaphos or control syrup. Also frozen bee bread and water were administered ad libitum and honeybees were maintained in an incubator (28 °C; 30% HR) for 9 days under the treatments mentioned above. Solutions were administered in gravity feeders and were replaced daily, recording consumption every day through differences in feeder weight. Solution evaporation was checked in order to correct the consumed volumes.

In summary, honeybees were treated as follows: non-infected without acaricide exposure (C), non-infected exposed to tau-fluvalinate (CFluv), non-infected exposed to coumaphos (CCou), Infected without acaricide exposure (I), infected exposed to tau-fluvalinate (I + Fluv), infected exposed to coumaphos (I + Cou).

Nosema ceranae development

Ten days post-infection, fifteen to thirty bees per replicate were analysed in order to individually quantify the parasite development. The digestive tract was removed by pinching the last abdominal segments and cutting the midgut. Midguts were stored at −20 °C until quantification. The number of spores in suspension (parasite intensity) was individually quantified with a hemocytometer under a light microscope [45].

Bee sample preparation for molecular analysis

For total RNA extraction, ten bees belonging to each experimental group were collected at day 13 after emergence, sliced in pieces and individually merged in 1 mL of RNA later® (Ambion, Inc. Austin, TX, USA) in order to avoid the degradation of mRNA.

RNA isolation and cDNA synthesis

Ten bees per group (3 or 4 individuals per replicate) were individually analysed. RNA later was partially removed from samples. Honeybees were homogenised using a sterile glass rod and a sterile plastic tube. An aliquot of the body homogenates obtained to evaluate gene expression was also analysed to determine the presence of Nosema spores under light microscopy. Then, 600 µL of RLT buffer (Qiagen) were added to homogenates. Total RNA was isolated using the RNeasy plus mini kit (Qiagen), according to the manufacturer’s instructions. Total RNA recovered was immediately used to generate first strand cDNA using the quanti tect reverse transcription kit (Qiagen), according to the manufacturer’s instructions. Resultant cDNA synthesized were stored at −20 °C.

Negative controls were added in parallel for each step (RNA extraction and reverse transcription reactions).

Gene and primer selection

In order to evaluate the effect of acaricides on immune-related genes of honeybees, transcript levels for the genes encoding the antimicrobial peptides abaecin (ABA), hymenoptaecin (HYM), defensin (DEF), the immune related proteins, glucose dehydrogenase (GLD), lysozyme (LYS), vitellogenin (VG) were assessed using primers previously described [25, 46, 47]. Gene transcripts were normalized relative to expression levels for the gene encoding ribosomal protein S5 (RPS5), a gene with constitutive expression across honeybee life stages [48] and disease status and β-actin, a moderately expressed housekeeping gene [47].

Primers used for quantitative PCR of honeybee transcripts are summarized in Additional file 1.

Real-time PCR

Real-time PCR reactions for the amplification of specific genes were performed in 96-well plates using a BIO-Rad CFX-96 thermal cycler (Bio-Rad Corp.). Reaction mixtures consisted of 1× QuantiTect SYBR Green PCR Master Mix, 0.5 µM of each primer (one pair of primers per reaction), RNAse free water and 5 µL of 1:10 diluted cDNA in a final volume of 25 µL.

PCR reactions were conducted using a thermal program consisting of an initial activation step at 50 °C for 2 min and 95 °C for 15 min, and 39 cycles of 94 °C for 15 s, 50 °C for 30 s and 72 °C for 30 s. Fluorescence was measured in the elongation step. Negative controls lacking template were run by excluding nucleic acids from reactions.

Specificity of the reaction was checked by analysis of the melting curve of the final amplified product (from 65 to 95 °C).

Normalization of the real-time data

Amplification results from different genes were expressed as the quantification cycle (Cq) value. According to a suggestion made by Vandesompele et al. [49], the geometric mean between the Cq of two reference genes (RPS5 and β-actin) was used for analysis. The expression ratio between each target gene and the geometric mean of reference genes was calculated according to the method described by Pfaffl [50].

Statistical analyses

Prior to analysis, graphical explanatory techniques were applied to the original data to identify outliers and collinearity [51, 52].

Abundance of parasites

The abundance of parasites observed for each honeybee specimen was initially modelled using a generalized linear mixed model (GLMM) with a Poisson distribution and log link function. Nevertheless, due to more variability (overdispersion) than expected under the assumed distribution, a negative binomial distribution was preferred, resulting in the following model:

$${\text{Parasites}}_{\text{ij}} \sim {\text{NB }}(\mu_{\text{ij}} ,{\text{ k}});{\text{ Log }}(\mu_{\text{ij}} ) \, = {\text{ Intercept}} \, + \, {\text{ factor }}\left( {\text{Treatment}} \right)_{\text{ij}} \, + \, {\text{ a}}_{\text{i}} ,{\text{ where a}}_{\text{i}} \sim {\text{ N}}\left( {0, \, \sigma^{ 2}_{\text{cage}} } \right).$$

The random intercept ai introduces a correlation structure between all observations from the same sampling cage [53, 54].

Survival analysis

A similar model was also performed using a GLMM with a binomial distribution. The response variable Y i is the number of bees that survived after treatment, of N i bees treated for a period of 13 consecutive days. To model this probability, we used the logistic link function of the form:

$${\text{logit}}\left( {\pi_{\text{ij}} } \right) \, = {\text{ Intercept }} + {\text{ factor }}\left( {\text{Treatment}} \right)_{\text{ij}} \, + {\text{ Days}}_{\text{ij}} \, + {\text{ a}}_{\text{i}} ,{\text{ where a}}_{\text{i}} \sim {\text{ N}}\left( {0, \, \sigma^{ 2} {\text{cage}}} \right).$$

The term ai is a random intercept for cagei and imposes a correlation on all observations from the same cage.

Consumption rates

The relationship between average food consumption (milligrams/hour/bee) versus day, treatment identity and their interaction, was defined using a linear mixed model (LME) [53, 55, 56]. As mentioned above, the term ai is the random effect and is assumed to be normally distributed with mean 0 and variances σ2. The term εij is the unexplained noise, and is assumed to be normally distributed with mean 0 and σ2.

The following model was applied:

$${\text{Food consumption}}_{\text{ij}} = {\text{ Intercept }} + {\text{ factor }}\left( {\text{Treatment}} \right)_{\text{ij}} *{\text{ Day}}_{\text{ij}} + {\text{ a}}_{\text{i}} + \, \varepsilon_{\text{ij}} .$$

Gene expression

Linear mixed effects models (LME) were used to model the gene expression as a response, in relation to the explanatory variable treatment. For each treatment level, gene expression was sampled from the same specimen. Hence, it is likely that individual gene expression is dependent on one another. Therefore to avoid a certain degree of pseudo-replication when fitting six separate models, an N by N correlation matrix was calculated a posteriori. The off-diagonal elements define the correlation between the residuals of each gene expression analysis and was used to assess dependency among the response variables. Thus, the following model was applied on the data obtained for each gene expression:

$${\text{Gene expression}}_{\text{ij}} = {\text{ Intercept }} + {\text{ factor }}\left( {\text{Treatment}} \right)_{\text{ij}} + {\text{ a}}_{\text{i}} + \, \varepsilon_{\text{ij}} .$$

The term ai is the random effect representing the between-cage variation and is assumed to be normally distributed with mean 0 and variances σ2. The term εij is the unexplained noise, and is assumed to be normally distributed with mean 0 and σ2. However, an initial analysis indicated large variation due to the effect of the covariate treatment, and therefore we allowed for heterogeneous residual variance structures [53, 55, 56].

The model selection followed the step-down approach described in West et al. [56] and Zuur et al. [53]. Significance of the variable treatment within the models was assessed using p values, which explains the effect of individual level of factor on the response variable as compared to the controls.

A detailed model validation was carried out by plotting residuals versus fitted values and versus each covariate to confirm that the underlying statistical assumptions were not violated [53, 54]. For the LME models, optimal model coefficients were derived using REML estimation [53, 56]. All analyses were performed with R (R Development Core Team, 2014) using the packages lme4 and nlme, respectively [57, 58].


Nosema ceranae infection on honeybees

Sublethal doses of coumaphos and tau-fluvalinate did not alter N. ceranae development. Statistical analysis revealed no significant differences between experimental groups (X2 = 1.48, df = 2, p = 0.47).

The recorded mean of spores ± standard errors per bee in infected controls was 11 534 261 ± 1 236 776; 12 104 406 ± 1 081 263 for coumaphos treated bees and 11 923 054 ± 1 114 173 for tau-fluvalinate treated ones. Less than 8% of artificial infected honeybees did not have spores in their midgut when analysed by light microscopy.

No N. ceranae spores were recorded in non-infected honeybees, demonstrating that no spore contamination in the bee bread diet and the absence of unwanted infection of bees occurred during the assays.

Honeybee survival

Model selection indicated that there was no significant effect of treatment identity on survival rates (X2 = 3.29, df = 5, p = 0.65) indicating that administrated coumaphos and tau-fluvalinate did not alter mortality of infected or non-infected honeybees. However, the probability of survival declined significantly after the start of the experiment (X2 = 534.14, df = 1, p < 0.001) (Figure 1).

Figure 1
figure 1

Honey bee survival. Predicted survival values according to each treatment and days after inoculation. Treatments: C (control); I (infected honeybees); CFluv (non-infected honeybees fed with 666 ppb of tau-fluvalinate); I + Fluv (infected honeybees fed with 666 ppb of tau-fluvalinate); CCou (non-infected honeybees fed with 3333 ppb of coumaphos); I + Cou (infected honeybees fed with 3333 ppb of coumaphos).

Food intake

The optimal model for food consumption was an LME, where food consumption was mainly affected by treatment variable (L-ratio = 22.23, df = 5, p < 0.001). There was no significant interaction between the days of consumption and treatments (L-ratio = 4.13, df = 5, p = 0.53). The effect of time (days) was neither significant (L-ratio = 0.86, df = 1, p = 0.37).

Significant differences were only detected on infected honeybees that received coumaphos treatment when compared to controls (p = 0.007). Figure 2 shows the results of LME for food consumption.

Figure 2
figure 2

LME results for food consumption. Sucrose consumption is expressed as the amount of sucrose solution (50% w/v, ad libitum delivery) consumed daily, per hour and per bee during the assay (food intake). The black dots represent the fitted values, and the lines are 95% confident intervals. Raw data were superimposed (grey dots). Treatments: C (control); I (infected honeybees); CFluv (non-infected honeybees fed with 666 ppb of tau-fluvalinate); I + Fluv (infected honeybees fed with 666 ppb of tau-fluvalinate); CCou (non-infected honeybees fed with 3333 ppb of coumaphos); I + Cou (infected honeybees fed with 3333 ppb of coumaphos). Asterisk indicates significant differences (df: 207; p = 0.0078) with control treatment (C).

Immune-related gene expression

Expression levels of genes involved in immune response were determined 10 days after infection (13 days after emergence). Gene-specific amplification was confirmed for the six genes analysed, as a single peak in the melting curve analysis and through the Tm values. No amplification occurred for negative controls (non-template controls).

Nosema ceranae spore presence was confirmed, under a light microscope, before total RNA extraction in all infected honeybee samples. In contrast, no spores were detected in non-infected honeybees.

Two of the six statistical models did not fulfil the assumptions of the LME and therefore we failed to validate them. This was mainly due to the presence of negative fitted values and/or heterogeneity of variances. Therefore no results are presented for hymenoptaecin and glucose-dehydrogenase gene expressions. With respect to the defensin gene expression, the data set contained a considerable number of missing values for a level of the factor treatment, therefore we could not include the “non-infected exposed to coumaphos” level in the analysis.

Genes that encode for antimicrobial peptides

Treatments did not show a significant effect on abaecin (L-ratio = 9.48, df = 5, p = 0.09) nor defensin transcript levels (L-ratio = 0.85, df = 5, p = 0.97).

Genes that encode lysozyme and vitellogenin

The optimal model describing vitellogenin mRNA levels was an LME that allowed for unequal variances depending on infection vs non-infection. The “non-infected exposed to coumaphos” treatment was removed from the statistical analysis due to a considerable number of missing values that were not accepted by the model. Vitellogenin gene expression was affected by treatments (L-ratio = 12.27, df = 4, p = 0.01). This is a significant, albeit weak effect (Figure 3).

Figure 3
figure 3

Vitellogenin mRNA expression. LME results. The black dots represent the fitted values, and the lines are 95% confidence intervals. Raw data were superimposed (grey dots). Gene expression was primarily affected by treatment (L-ratio = 12.27, df = 4, p = 0.01). Treatments: C (control); I (infected honeybees); CFluv (non-infected honeybees fed with 666 ppb of tau-fluvalinate); I + Fluv (infected honeybees fed with 666 ppb of tau-fluvalinate); I + Cou (infected honeybees fed with 3333 ppb of coumaphos). Treatments I, I + Fluv and I + Cou were significantly lower than control treatment (t = −2.62, t = −2.13 and t = −2.47 respectively; p = 0.01, p = 0.04 and p = 0.02 respectively). Asterisk indicates statistically significant differences with treatment (C).

Nosema ceranae infection caused a reduction in transcript abundances of vitellogenin (VG) when compared to the control treatment (p = 0.01). Also, when both stressors were applied a drop of VG was detected (p = 0.02).

The tau-fluvalinate treatment did not cause a reduction by itself in VG gene expression; however the combination of this acaricide and N. ceranae infection did (p = 0.04) (Figure 3).

In the case of lysozyme expression analysis, the optimal model was also an LME with unequal variances. The results revealed that lysozyme transcript levels were primarily affected by treatment (L-ratio = 20.92, df = 5, p < 0.001). When comparing the treatment levels with the control, we found that non-infected honeybees exposed to coumaphos expressed significantly lower relative abundances of transcript compared to the baseline (t = −3.09; p < 0.01), and that non-infected honeybees that received tau-fluvalinate treatment also had a negative but weak influence on gene expression (t = −2.09; p = 0.05) (Figure 4).

Figure 4
figure 4

Lysozyme mRNA expression. Results of the LME. The black dots represent the fitted values, and the lines are 95% confidence intervals. Raw data were superimposed (grey dots). Gene expression was primarily affected by treatment (L-ratio = 20.92, df = 5, p < 0.001). Treatments: C (control); I (infected honeybees); CFluv (non-infected honeybees fed with 666 ppb of tau-fluvalinate); I + Fluv (infected honeybees fed with 666 ppb of tau-fluvalinate); CCou (non-infected honeybees fed with 3333 ppb of coumaphos); I + Cou (infected honeybees fed with 3333 ppb of coumaphos). When compared with control, treatment CCou was significantly lower (t = −3.09; p < 0.01) and treatment CFluv also had a negative but weak influence on the gene expression (t = −2.09; p = 0.05). The asterisk indicates statistically significant differences with respect to control treatment (C).


In spite of the ubiquitous presence of N. ceranae and highly persistent acaricides in apicultural matrices, our results revealed under carefully controlled laboratory conditions, no evidence of combined effects of those stressors on honeybee survival and in most of the immune related genes analysed.

It is worth noting that coumaphos and tau-fluvalinate concentrations provided in our study have been reported in contaminated honey and pollen stored inside hives [17, 18, 45, 46], therefore acaricide-honeybee interactions are likely to occur under field-realistic scenarios. Chronic administration of these acaricides during 9 days did not affect worker survival, nor in coexposure with N. ceranae under laboratory conditions. Our test bees only received acaricide exposure as imagoes, and larvae developed in a reduced acaricide environment, therefore it will be interesting to analyse exposure to sublethal levels on survival even during larval development.

It is unsurprising that insecticides alter insect behaviour, such as feeding levels but, in this work, intake of diets was similar to control group. In other words, the observed responses were not due to poor nutrition. Only infected bees that received syrup contaminated with coumaphos showed a reduced consume when compared to controls. This result might be clarified by performing long lasting assays.

Although colonies employed during the assays did not show any clinical symptom of diseases, pathogens could be present in the administered inoculum or in the pollen diet. However, until this work was conceived, no published method have been proposed in order to avoid the possible presence of pathogens without altering the highest nutritional quality of beebread. This should be considered in future experimental designs.

Recently, an increase in Nosema spore loads in colonies treated with the insecticide imidacloprid was shown, suggesting an indirect effect of pesticides on pathogen growth in honeybees [59]. In our study, Nosema growth rate did not change in the presence of coumaphos or tau-fluvalinate at sublethal doses, in agreement with recent field work in naturally infected colonies [60].

Nosema ceranae infection by itself did not have an effect on genes that encode antimicrobial peptides nor tested enzymes at day 10 post infection. Only a significant downregulation of vitellogenin was detected in N. ceranae infected honeybees. These results agreed with those found by Antúnez et al. [61] who showed that vitellogenin transcripts decreased in N. ceranae infected bees at 7 days post infection. However, Chaimanee et al. [62] did not find differences in vitellogenin expression levels at day 13 post emergence but methodologies were dissimilar, mainly since they did not provide any protein source. In reference to the importance of the diet, a positive correlation between the level of pollen consumption and Vg expression in Nosema-infected bees has been widely demonstrated, underlining the importance of the diet [63].

Although a reduction in Vg transcripts was recorded in infected honeybees treated with tau-fluvalinate and coumaphos, this was due to the parasite action. The vitellogenin protein is a central regulator of lifespan, supporting hemocyte viability [32], recognizing cell damage [64] and acting as a metabolic regulator of stress responses [33]; altering its expression rates induced by the presence of N. ceranae, independently of acaricide exposure, could lead to an acceleration on honeybee ontogeny [33] and reduction of its antioxidant capacity. The limited knowledge of the impact of chemical pesticides on humoral immunity postulates that they do not appear to affect AMP production [65]. This agrees with our study, since no differences in defensin or abaecin transcript levels were found on bees exposed to tau-fluvalinate or coumaphos treatments.

In the present study, effects on the lysozyme gene on honeybees at sublethal chronic exposure were recorded but in a previous report, no change was observed under acute exposure (LC50) employing the same acaricides [23]. Although it should be considered that different methodologies related to acaricide exposure were employed, such as the evaluation of acute toxicity by contact, we can hypothesize that huge amounts of these molecules administrated over 24 h do not impact on the immune system as much as little doses over a longer period of time; it seems that chronic stress can lead to a dysregulation of the physiological pathways to compensate for the stress response.

The common feature of lysozymes in animals is their ability to hydrolyse bacterial cell wall [66] and to promote the synthesis of other antimicrobial peptides [67]. Chronic exposure to coumaphos reduced lysozyme transcripts but in exposed honeybees infected with N. ceranae no differences were observed with respect to control treatments. This result may be due to the absence of response by individual honeybees in the presence of both stressors or that the parasite induces an up-regulation of this gene while the acaricide downregulates it, buffering the net effect. The energetic cost could be the key reason behind this result because this parasite generates an energetic demand on its host. This could affect the mechanisms that regulate the mobilisation of energy reserves in infected individuals [68, 69].

The long term effects of tau-fluvalinate and coumaphos and their interaction with other ubiquitous parasites are of particular concern. Thus, the possible detrimental effects of these stressors on colony health and performance cannot be discarded and warrant further research. Our interest in future studies will be to intensify sampling in order to reduce variability between individuals and to perform again those treatments or response variables that were necessary to exclude from statistical analyses. Also, since honeybee physiology is characterized by a surprising plasticity and is regulated by colony signals and resources, further field or semi-field studies are necessary.


  1. Shimanuki H, Calderone NW, Knox DA (1994) Parasitic mite syndrome: the symptoms. Am Bee J 134:827–828

    Google Scholar 

  2. Boecking O, Genersch E (2008) Varroosis—the ongoing crisis in bee keeping. J Verbr Lebensm 3:221–228

    Article  Google Scholar 

  3. Evans JD, Schwarz RS (2011) Bees brought to their knees: microbes affecting honeybee health. Trends Microbiol 12:614–620

    Article  Google Scholar 

  4. Dainat B, Evans JD, Chen YP, Gauthier L, Neumann P (2012) Predictive markers of honeybee colony collapse. PLoS One 7:e32151

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Chauzat MP, Martel AC, Blanchard P, Clément MC, Schurr F, Lair C, Ribière M, Wallner K, Rosenkranz P, Faucon JP (2010) A case report of a honeybee colony poisoning incident in France. J Apic Res 49:113–115

    Article  Google Scholar 

  6. Medrzycki P, Sgolastra F, Bortolotti L, Bogo G, Tosi S, Padovani E, Porrini C, Sabatini AG (2010) Influence of brood rearing temperature on honeybee development and susceptibility to poisoning by pesticides. J Apic Res 49:52–59

    Article  CAS  Google Scholar 

  7. Paxton RJ, Klee J, Korpela S, Fries I (2007) Nosema ceranae has infected Apis mellifera in Europe since at least 1998 and may be more virulent than Nosema apis. Apidol 38:558–565

    Article  Google Scholar 

  8. Cox-Foster DL, Conlan S, Holmes EC, Palacios G, Evans JD, Moran NA, Quan PL, Briese T, Hornig M, Geiser DM, Martinson V, vanEngelsdorp D, Kalkstein AL, Drysdale A, Hui J, Zhai J, Cui L, Hutchison SK, Simons JF, Egholm M, Pettis JS, Lipkin WI (2007) A metagenomic survey of microbes in honeybee colony collapse disorder. Science 318:283–287

    Article  CAS  PubMed  Google Scholar 

  9. Higes M, Garcia-Palencia P, Martin-Hernandez R, Meana A (2007) Experimental infection of Apis mellifera honeybees with Nosema ceranae (Microsporidia). J Invertebr Pathol 94:211–217

    Article  PubMed  Google Scholar 

  10. Paxton RJ (2010) Does infection by Nosema ceranae cause “Colony Collapse Disorder” in honeybees (Apis mellifera)? J Apic Res 49:80–84

    Article  Google Scholar 

  11. Fries I, Feng F, Silva A, Da Slemenda SB, Pieniazek NJ (1996) Nosema ceranae (Microspora, Nosematidae), morphological and molecular characterization of a microsporidian parasite of the Asian honeybee Apis cerana (Hymenoptera, Apidae). Eur J Protistol 32:356–365

    Article  Google Scholar 

  12. Klee J, Besana AM, Genersch E, Gisder S, Nanetti A, Tam DQ, Chinh TX, Puerta F, Ruz JM, Kryger P, Message D, Hatjina F, Korpela S, Fries I, Paxton RJ (2007) Widespread dispersal of the Microsporidian Nosema ceranae, an emergent pathogen of the western honeybee, Apis mellifera. J Invertebr Pathol 96:1–10

    Article  PubMed  Google Scholar 

  13. Rosenkranz P, Aumeier P, Ziegelmann B (2010) Biology and control of Varroa destructor. J Invertebr Pathol 103:96–119

    Article  Google Scholar 

  14. Bogdanov S, Kilchenmann V, Imdorf A (1996) Acaricide residues in beeswax and honey. In: Mizrahi A, Lensky Y (eds) Bee products: properties, applications and apitherapy. Plenum Press, New York, pp 239–247

    Google Scholar 

  15. Lodesani M, Costa C, Bigliardi M, Colombo R (2003) Acaricide residues in bee wax and organic beekeeping. Apiacta 38:31–33

    Google Scholar 

  16. Chauzat MP, Faucon JP (2007) Pesticide residues in beeswax samples collected from honeybee colonies (Apis mellifera L.) in France. Pest Manag Sci 63:1100–1106

    Article  CAS  PubMed  Google Scholar 

  17. Martel AC, Zeggane S, Aurieres C, Drajnudel P, Faucon JP, Aubert M (2007) Acaricide residues in honey and wax after treatment of honeybee colonies with Apivar or Asuntol 50. Apidol 38:534–544

    Article  CAS  Google Scholar 

  18. Mullin CA, Frazier M, Frazier JL, Ashcraft S, Simonds R, van Engelsdorp D, Pettis JS (2010) High levels of miticides and agrochemicals in North American apiaries: implications for honey bee health. PLoS One 5:e9754

    Article  PubMed  PubMed Central  Google Scholar 

  19. Orantes-Bermejo FJ, Gomez-Pajuelo A, Megias-Megias M, Torres Fernández- Pinar C (2010) Pesticide residues in beeswax and beebread samples collected from honeybee colonies (Apis mellifera L.) in Spain. Possible implications for bee losses. J Apic Res 49:243–250

    Article  CAS  Google Scholar 

  20. Locke B, Forsgren E, Fries I, de Miranda JR (2012) Acaricide treatment affects viral dynamics control colonies via both host physiology and mite in Varroa destructor infested honeybee. Appl Environ Microbiol 78:227–235

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. vanBuren NWM, Marien AGH, Oudejans RCHM, Velthuis HHW (1992) Perizin, an acaricide to combat the mite Varroa jacobsoni: its distribution in and influence on the honeybee Apis mellifera. Physiol Entomol 17:288–296

    Article  CAS  Google Scholar 

  22. Boncristiani H, Underwood R, Schwarz R, Evans JD, Pettis J, vanEngelsdorp D (2012) Direct effect of acaricides on pathogen loads and gene expression levels of honeybee Apis mellifera. J Insect Physiol 58:613–620

    Article  CAS  PubMed  Google Scholar 

  23. Garrido PM, Antúnez K, Martín M, Porrini MP, Zunino P, Eguaras MJ (2013) Immune- related gene expression in nurse honeybees (Apis mellifera) exposed to synthetic acaricides. J Insect Physiol 59:113–119

    Article  CAS  PubMed  Google Scholar 

  24. Lavine MD, Strand MR (2002) Insect haemocytes and their role in immunity. Insect Biochem Mol Biol 32:1295–1309

    Article  CAS  PubMed  Google Scholar 

  25. Evans JD (2006) Beepath: an ordered quantitative-PCR array for exploring honeybee immunity and disease. J Invertebr Pathol 93:135–139

    Article  CAS  PubMed  Google Scholar 

  26. Casteels P, Ampe C, Riviere L, Van Damme J, Elicone C, Fleming M, Jacobs F, Tempst P (1990) Isolation and characterization of abaecin, a major antibacterial response peptide in the honeybee (Apis mellifera). Eur J Biochem 187:381–386

    Article  CAS  PubMed  Google Scholar 

  27. Casteels P, Ampe C, Jacobs F, Tempst P (1993) Functional and chemical of Hymenoptaecin, an antibacterial polypeptide that is infection-inducible in the honeybee (Apis mellifera). J Biol Chem 268:7044–7054

    CAS  PubMed  Google Scholar 

  28. Jiravanichpaisal P, Lee BL, Söderhäll K (2006) Cell-mediated immunity in arthropods: hematopoiesis, coagulation, melanization and opsonization. Immunobiology 211:213–236

    Article  CAS  PubMed  Google Scholar 

  29. Decker H, Jaenicke E (2004) Recent findings on phenoloxidase activity and antimicrobial activity of hemocyanins. Dev Comp Immunol 28:673–687

    Article  CAS  PubMed  Google Scholar 

  30. Cox-Foster DL, Stehr JE (1994) Induction and localization of FAD-glucose dehydrogenase (GLD) during encapsulation of abiotic implants in Manduca sexta larvae. J Insect Physiol 40:235–249

    Article  CAS  Google Scholar 

  31. Wheeler DE, Kawooya JK (1990) Purification and characterization of honeybee vitellogenin. Arch Insect Biochem Physiol 14:253–267

    Article  CAS  PubMed  Google Scholar 

  32. Amdam GV, Simoes ZL, Hagen A, Norberg K, Schrøder K, Mikkelsen Ø, Kirkwood TBL, Omholt SW (2004) Hormonal control of the yolk precursor vitellogenin regulates immune function and longevity in honeybees. Exp Gerontol 39:767–773

    Article  CAS  PubMed  Google Scholar 

  33. Amdam GV (2011) Social context, stress, and plasticity of aging. Aging Cell 10:18–27

    Article  CAS  PubMed  Google Scholar 

  34. Antúnez K, Mendoza Y, Santos E, Invernizzi C (2013) Differential expression of vitellogenin in honeybees (Apis mellifera) with different degrees of Nosema ceranae infection. J Apic Res 52:227–234

    Article  Google Scholar 

  35. Fries I, Aarhus A, Hansen H, Korpela S (1991) Comparison of diagnostic methods for detection of low infestation levels of Varroa jacobsoni in honeybee (Apis mellifera) colonies. Exp Appl Acarol 10:279–287

    Article  Google Scholar 

  36. Fries I, Ekbohm G, Villumstad E (1984) Nosema apis, sampling techniques and honey yield. J Apic Res 23:102–105

    Google Scholar 

  37. Serra-Bonvehí J, Orantes-Bermejo J (2010) Acaricides and their residues in Spanish commercial beeswax. Pest Manag Sci 66:1230–1235

    Article  PubMed  Google Scholar 

  38. Martín-Hernández R, Meana A, Prieto L, Martinez Salvador A, Garrido-Bailon E, Higes M (2007) Outcome of colonization of Apis mellifera by Nosema ceranae. Appl Environ Microbiol 73:6331–6338

    Article  PubMed  PubMed Central  Google Scholar 

  39. Fries I, Chauzat MP, Chen YP, Doublet V, Genersch E, Gisder Higes M, McMahon DP, Martín-Hernández R, Natsopoulou M, Paxton RJ, Tanner G, Webster TC, Williams GR (2013) Standard methods for Nosema research. J Apic Res 52:1–28

    Article  Google Scholar 

  40. Suwannapong G, Yemor T, Boonpakdee C, Benbow ME (2010) Nosema ceranae, a new parasite in Thai honeybees. J Invertebr Pathol 106:236–241

    Article  PubMed  Google Scholar 

  41. Porrini MP, Garrido PM, Eguaras MJ (2013) Individual feeding of honeybees: modification of the Rinderer technique. J Apic Res 52:194–195

    Article  Google Scholar 

  42. Atienza J, Jimenez JJ, Bernal L, Martin MT (1993) Supercritical fluid extraction of fluvalinate residues in honey. Determination by high-performance chromatography. J Chromatogr 655:95–99

    Article  CAS  PubMed  Google Scholar 

  43. Chauzat MP, Faucon JP, Martel AC, Lachaize J, Cougoule N, Aubert M (2006) A survey on pesticide residues in pollen loads collected by honey-bees (Apis mellifera) in France. J Econ Entomol 99:253–262

    Article  CAS  PubMed  Google Scholar 

  44. Mustard JA, Edgar EA, Mazade RE, Wu C, Lillvis JL, Wright GA (2008) Acute ethanol ingestion impairs appetitive olfactory learning and odor discrimination in the honey bee. Neurobiol Learn Mem 90:633–643

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Cantwell GE (1970) Standard methods for counting Nosema spores. Am Bee J 110:222–223

    Google Scholar 

  46. Corona M, Velarde R, Remolina S, Moran-Lauter A, Wang Y, Hughes KA, Robinson GE (2007) Vitellogenin, juvenile hormone, insulin signalling, and queen honeybee longevity. Proc Natl Acad Sci U S A 104:7128–7133

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Yang X, Cox-Foster DL (2005) Impact of an ectoparasite on the immunity and pathology of an invertebrate: evidence for host immunosuppression and viral amplification. Proc Natl Acad Sci U S A 21:7470–7475

    Article  Google Scholar 

  48. Evans JD, Wheeler DE (2000) Expression profiles during honeybee caste determination. Genome Biol 2:0001

    Article  Google Scholar 

  49. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 3:0034

    Article  Google Scholar 

  50. Pfaffl M (2001) A new mathematical model for relative quantification in real time RT PCR. Nucleic Acids Res 29:2002–2007

    Article  Google Scholar 

  51. Ieno EN, Zuur AF. A Beginner’s Guide to Data Exploration and Visualisation with R. Highland Statistics Ltd. 2015

  52. Zuur AF, Ieno EN, Elphick CS (2010) A protocol for data exploration to avoid common statistical problems. Methods Ecol Evol 1:3–14

    Article  Google Scholar 

  53. Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith G. Mixed effects models and extensions in ecology with R. Springer. New York, 2009

    Book  Google Scholar 

  54. Zuur AF, Hilbe JM, Ieno EN. A Beginner’s guide to GLM and GLMM with R. A frequentist and Bayesian perspective for ecologists. Highland Statistics Ltd., Newburgh, 2013. p 256

    Google Scholar 

  55. Pinheiro J, Bates D (2000) Mixed effects models in S and S-plus. Springer, New York

    Book  Google Scholar 

  56. West BT, Welch KB, Gałecki AT (2006) Linear mixed models. A practical guide using statistical software; Chapman and Hall/CRC, Boca Raton

    Google Scholar 

  57. Bates D, Maechler M, Bolker B, Walker S. lme4: linear mixed-effects models using Eigen and S4. R package version 1.1-8. 2015.

  58. Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team nlme: linear and Nonlinear Mixed Effects Models. R package version 3.1-120. 2015.

  59. Pettis JS, van Engelsdorp D, Johnson J, Dively G (2012) Pesticide exposure in honeybees results in increased levels of the gut pathogen Nosema. Naturwissenschaften 99:153–158

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Berry JA, Hood WM, Pietravalle S, Delaplane KS (2013) Field-level sublethal effects of approved bee hive chemicals on honeybees (Apis mellifera L). PLoS One 8:e76536

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Antúnez K, Martin-Hernandez R, Prieto L, Meana A, Zunino P, Higes M (2009) Immune suppression in the honeybee (Apis mellifera) following infection by Nosema ceranae (Microsporidia). Environ Microbiol 11:2284–2290

    Article  PubMed  Google Scholar 

  62. Chaimanee V, Chantawannakul P, ChenY Evans JD, Pettis JS (2012) Differential expression of immune genes of adult honeybee (Apis mellifera) after inoculated by Nosema ceranae. J Insect Physiol 58:1090–1095

    Article  CAS  PubMed  Google Scholar 

  63. Zheng HQ, Lin ZG, Huang SK, Sohr A, Wu Z, Chen YP (2014) Spore Loads may not be used alone as a direct indicator of the severity of Nosema ceranae infection in honey bees Apis mellifera (Hymenoptera:Apidae). J Econom Entomol 107:2037–2044

    Article  Google Scholar 

  64. Havukainen H, Münch M, Baumann A, Zhong S, Halskau Ø, Krogsgaard M, Amdam GV (2013) Vitellogenin recognizes cell damage through membrane binding and shields living cells from reactive oxygen species. J Biol Chem 288:28369–28381

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. James RR, Xu J (2012) Mechanisms by which pesticides affect insect immunity. J Invertebr Pathol 109:175–182

    Article  CAS  PubMed  Google Scholar 

  66. Callewaert L, Michiels CW (2010) Lysozymes in the animal kingdom. J Biosci 35:127–160

    Article  CAS  PubMed  Google Scholar 

  67. Imler JL, Bulet P (2005) Antimicrobial peptides in Drosophila: structures, activities, and gene regulation. In: Kabelitz D, Schorder JM (eds) Mechanism of Epithelial Defense, vol 86. Karger, Basel, pp 1–21

    Chapter  Google Scholar 

  68. Mayack C, Naug D (2009) Energetic stress in the honeybee Apis mellifera from Nosema ceranae infection. J Invertebr Pathol 100:185–188

    Article  PubMed  Google Scholar 

  69. Aliferis KA, Copley T, Jabali S (2012) Gas chromatography-mass spectrometry metabolite profiling of worker honeybee (Apis mellifera L.) hemolymph for the study of Nosema ceranae infection. J Insect Physiol 58:1349–1359

    Article  CAS  PubMed  Google Scholar 

Download references

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

PMG, MPP, GMN, MJE conceived and designed the experiments. PMG, MPP, KA, BB performed the experiments and analysed the data. PMG, MPP, KA and EI performed the statistical analysis. PMG, KA, EI, PZ, GS wrote the article or revised it critically. All authors read and approved the final manuscript.


The authors would like to thank the National Research Council (CONICET), UNMdP (National University of Mar del Plata) and Bilateral Project UR1201 MinCyT-MEC. In addition, we thank the anonymous reviewers for their constructive comments which helped to improve the quality of the manuscript.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Paula Melisa Garrido.

Additional file


Additional file 1. Sequences of oligonucleotide primers used for real-time PCR quantification. Table provides primer sequences, melting temperatures and reaction efficiencies.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Garrido, P.M., Porrini, M.P., Antúnez, K. et al. Sublethal effects of acaricides and Nosema ceranae infection on immune related gene expression in honeybees. Vet Res 47, 51 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: