Pollinators are exposed to numerous parasites and pathogens when foraging on flowers. These biological stressors may affect critical cognitive abilities required for foraging. Here, we tested whether exposure to Nosema ceranae, one of the most widespread parasites of honey bees also found in wild pollinators, impacts cognition in bumblebees. We investigated different forms of olfactory learning and memory using conditioning of the proboscis extension reflex. Seven days after being exposed to parasite spores, bumblebees showed lower performance in absolute, differential and reversal learning than controls. The consistent observations across different types of olfactory learning indicate a general negative effect of N. ceranae exposure that did not specifically target particular brain areas or neural processes. We discuss the potential mechanisms by which N. ceranae impairs bumblebee cognition and the broader consequences for populations of pollinators.

Pollinators, such as bees, rely on a rich cognitive repertoire to collect pollen and nectar on flowers. These include associative learning and memories of floral traits such as odours, shapes, colours and textures, to identify the most profitable resources (Giurfa, 2015; Menzel, 2012), and spatial cues for navigation (Collett et al., 2013). Any disruption of these cognitive abilities by environmental stressors can considerably reduce the foraging performance of bees, ultimately compromising brood development and survival (Klein et al., 2017).

In particular, foraging bees are exposed to a number of parasites that can affect their physiology and behaviour (Gómez-Moracho et al., 2017). The microsporidium Nosema ceranae is one of the most prevalent parasites of bees worldwide with a large range of hosts including honey bees (Higes et al., 2006), bumblebees (Plischuk et al., 2009), solitary bees (Ravoet et al., 2014), but also other flower visitors such as wasps (Porrini et al., 2017). Insects become infected by ingesting parasite spores from contaminated water or pollen (Higes et al., 2008), or during physical contact with contaminated individuals (Smith, 2012). The spores invade the gut epithelial cells of the hosts, where they develop (Holt et al., 2013). In honey bees, N. ceranae degenerates the gut epithelium (Higes et al., 2007), alters metabolism (Mayack and Naug, 2009) and disrupts the immune response (Antúnez et al., 2009). This causes a disease (nosemosis) believed to contribute to colony collapse (Cox-Foster et al., 2007).

Nosema ceranae-infected honey bees also show impaired navigation (Wolf et al., 2014) and increased flight activity (Dussaubat et al., 2013), suggesting that their cognitive abilities are affected by the parasite. Recent studies have explored this possibility from a mechanistic point of view using Pavlovian olfactory conditioning of the proboscis extension reflex (PER; Takeda, 1961) in which harnessed bees are trained to associate an odour, or a combination of odours, with a sucrose reward (for recent reviews, see Lavond and Steinmetz, 2003; Matsumoto et al., 2012). However, their results are mixed, presumably because of important variations in parasite exposure protocols, age of bees and parasite post-exposure duration in the different studies (Bell et al., 2020; Charbonneau et al., 2016; Gage et al., 2018; Piiroinen and Goulson, 2016; Piiroinen et al., 2016). Only two studies explored these effects in bumblebees. One study suggests a slight impairment of absolute learning (Piiroinen and Goulson, 2016), and both report no effect on memory (Piiroinen et al., 2016; Piiroinen and Goulson, 2016). Note however that in these two studies less than 3% of the bumblebees exposed to N. ceranae were indeed found to be contaminated by the parasite (i.e. PCR positive) after the behavioural tests.

Given the expanding geographical distribution of N. ceranae worldwide (Klee et al., 2007), its increasing prevalence in wild bees (Plischuk et al., 2009; Porrini et al., 2017; Ravoet et al., 2014), and the potential high fitness costs incurred by bees with impaired cognition (Henry et al., 2012; Klein et al., 2017; Perry et al., 2015), clarifying its influence on host learning and memory is important for risk assessment. In particular, other critical forms of learning, such as the ability to associate one of two odours with a reward (differential learning) and reverse this association (reversal learning), have so far been unexplored. These types of learning are essential in the everyday life of bees, to discriminate flowers, olfactory landmarks and social partners, and require different brain centres [e.g. functional mushroom bodies are necessary for the acquisition of non-elemental associations but not for elemental associations (Boitard et al., 2015; Devaud et al., 2007, 2015; Giurfa and Sandoz, 2012)]. If the effects of the parasite are specific, these types of learning may be more or less impacted. Conversely, if the effects of the parasite are general, all learning types may be impacted.

Here, we built on a recently established method yielding high rates of experimental infection by N. ceranae (Gómez-Moracho et al., 2021) to study the impact of the parasite on different cognitive tasks in bumblebees. We used PER conditioning to compare the olfactory learning and memory performance of control bumblebees, bumblebees exposed to the parasite and bumblebees contaminated by the parasites (PCR positive) at 7 days post-exposure.

Bumblebees

We used bumblebee workers, Bombus terrestris (Linnaeus 1758), from 14 commercial colonies acquired from Biobest (Westerlo, Belgium). Before the experiments, we verified the absence of N. ceranae (Martín-Hernández et al., 2007), and other common parasites [Nosema bombi (Klee et al., 2006); Crithidia bombi (Schmid-Hempel and Tognazzo, 2010)] in a PCR using 15 bumblebees from each colony. We maintained bumblebees in their original colonies with ad libitum access to the syrup provided by the manufacturer and germ-free pollen (honey bee-collected pollen exposed to UV light for 12 h), in a room at 25±1°C under a 12 h:12 h light:dark photocycle, until parasite exposure.

Nosema ceranae spores

We obtained fresh spores from naturally infected honey bee colonies (Apis mellifera) maintained at our experimental apiary (University Toulouse III, France). To prepare spore solutions, we dissected the gut from 15 honey bees and crushed them in 15 ml of distilled H2O. We confirmed the presence of N. ceranae and the absence of Nosema apis (another common parasite of honey bees) in each homogenate by PCR (Martín-Hernández et al., 2007), and purified them following standard protocols (Fries et al., 2013). We centrifuged homogenates in aliquots of 1 ml at 5000 rpm for 5 min and re-suspended the pellet in 500 µl of dH2O by vortexing. This was repeated 3 times to obtain a spore solution of 85% purity (Fries et al., 2013). We counted N. ceranae spores in an improved Neubauer haemocytometer (Cantwell, 1970) in a light microscope (×400) and adjusted the spore inoculum to 15,000 spores µl−1 in 20% (w/w) of sucrose solution. Spore solutions were used within the same week they were purified.

Parasite exposure and experimental conditions

We exposed bumblebees to N. ceranae as described in Gómez-Moracho et al. (2021). Briefly, we confined individual bumblebees in a Petri dish for 5 h without food. We then exposed some bumblebees to a 20 µl drop of 20% sucrose solution containing 300,000 N. ceranae spores. Control bumblebees received 20 µl of sucrose solution (20% w/w). We only used bumblebees that consumed the entire drop of sucrose within the next 2 h. We then allocated bumblebees into microcolonies of 20–25 individuals, containing a gravity feeder with ad libitum access to food (Kraus et al., 2019). As diet can affect host–parasite relationships (Frost et al., 2008), we provided bumblebees with an artificial diet with a protein to carbohydrate ratio of 1:207, which has previously been shown to elicit highest N. ceranae prevalence in bumblebees (Gómez-Moracho et al., 2021). The diet was made with a fixed total amount of nutrients of 170 g l−1 (protein+carbohydrates) and 0.5% vitamin mixture for insects (Sigma). Carbohydrates were supplied as sucrose (Euromedex). Proteins consisted of a mixture of casein and whey (4:1) (Nutrimuscle) (Gómez-Moracho et al., 2021). We kept bumblebee microcolonies in a room at 25±1°C with a 12 h light:12 h dark photoperiod until the behavioural tests. Every day, we renewed the diet and removed dead bumblebees.

Behavioural experiments

We tested the cognitive performance of bumblebees of unknown age using PER at day 7 after parasite exposure (Fig. 1A). The day before the behavioural tests, we kept diets to low levels (∼200 µl per bumblebee) to keep bumblebees motivated for the PER experiments. Three hours before the behavioural tests, we collected bumblebees from the microcolonies, chilled them on ice for 5 min and restrained them in a modified 2 ml Eppendorf tube (hereafter, capsule) that we cut in length to fit each bumblebee (adapted from Toda et al., 2009; Fig. 1B). Bumblebees were tested in the horizontal position and could move forward and backward inside the tube, and therefore retract their head (Toda et al., 2009). We found these conditions better suited to perform PER experiments with bumblebees than the classical vertical harnessing used for honey bees (Giurfa and Sandoz, 2012), in which bumblebees appeared paralysed (T.G.-M., T.D. and M.L., unpublished data). In this approach, we obtained comparable to better learning performance than in previous studies [e.g. 58.1% in our study versus 44% in Laloi et al. (1999) or 57% in Piiroinen et al. (2016)]. Once the bumblebees were in the capsule, we kept them in the dark, in an incubator at 28°C, with no access to food. Bumblebees were left in the capsules for 3 h before the experiments, and for the whole duration of each conditioning protocol [i.e. a total of 4 h for sucrose responsiveness test, 5 h for absolute learning with short-term memory (STM) test, 6 h for reversal learning, and 28 h for absolute learning with 24 h memory (long-term memory, LTM) test]. All bumblebees that finished the conditioning protocols were kept at −20°C for later analysis of their infection status through PCR.

Fig. 1.

Conditioning protocols and odour delivery setup. (A) Schematic representation of the proboscis extension reflex (PER) protocols used in cognitive assays using two odorants (A and B). (i) Sequences used in the absolute learning and memory tasks. (ii) Sequences used for rewarded and unrewarded trials in the differential and reversal tasks. (iii) Sequence of events used in every trial. White bars represent odourless air flow before and after conditioning. Odour (CS) alone is represented by right diagonal lines. Sucrose (US) alone is represented by left diagonal lines. Hatching indicates the overlap of CS and US presentation. (B) Odour delivery setup. The bumblebee is placed inside the capsule in front of the air delivery setup. After the odour is delivered, sucrose is presented with a toothpick to the bumblebee, which extends its proboscis to drink the reward.

Fig. 1.

Conditioning protocols and odour delivery setup. (A) Schematic representation of the proboscis extension reflex (PER) protocols used in cognitive assays using two odorants (A and B). (i) Sequences used in the absolute learning and memory tasks. (ii) Sequences used for rewarded and unrewarded trials in the differential and reversal tasks. (iii) Sequence of events used in every trial. White bars represent odourless air flow before and after conditioning. Odour (CS) alone is represented by right diagonal lines. Sucrose (US) alone is represented by left diagonal lines. Hatching indicates the overlap of CS and US presentation. (B) Odour delivery setup. The bumblebee is placed inside the capsule in front of the air delivery setup. After the odour is delivered, sucrose is presented with a toothpick to the bumblebee, which extends its proboscis to drink the reward.

Sucrose responsiveness

We tested the sucrose responsiveness of bumblebees to control for potential influences of N. ceranae on their reward perception or feeding motivation. We presented seven sucrose solutions to each bumblebee, from concentrations of 0% (pure water) to 60% (w/w), with increments of 10% (Graystock et al., 2013). For each concentration, we touched the antennae of the bumblebee with a toothpick soaked in the corresponding sucrose solution to elicit the PER. We presented solutions in an increasing concentration gradient with an inter-trial interval of 5 min between concentrations. We discarded all the bumblebees responding to water (i.e. 0% sucrose solution) to avoid the effect of thirst on sucrose responsiveness (Baracchi et al., 2018), and those showing an inconsistent response (i.e. bumblebees responding to lower but not higher sucrose concentrations; Table 1) (Scheiner et al., 2003). As sucrose concentrations were systematically presented in the same increasing order to the bumblebees, we calculated a gustatory score for each bumblebee based on the lowest concentration at which they responded. Bumblebees whose first response was observed following exposure to 10% sucrose had a score of 1, whereas bumblebees that responded for the first time to 60% sucrose had a score of 6. Therefore, the lower the score, the lower the sucrose sensitivity threshold of the bumblebee.

Table 1.

Sample sizes for each experiment 7 days post-exposure

Sample sizes for each experiment 7 days post-exposure
Sample sizes for each experiment 7 days post-exposure

Conditioning experiments

All experiments shared the same general protocol (Fig. 1A). An encapsulated bumblebee (Fig. 1B) was placed 1.5 cm in front of an automated conditioning setup (described in Aguiar et al., 2018) delivering a continuous stream of odourless air at 1.2 ml s−1 to which specific odours were selectively added (Raiser et al., 2017). We used two odorants as conditioned stimulus (CS): nonanal and phenylacetaldehyde (Palottini et al., 2018; Sommerlandt et al., 2014), at a 1:100 dilution in mineral oil. Before conditioning, we tested the responsiveness of bumblebees to sucrose by touching both antennae with a toothpick soaked in 50% (w/w) sucrose solution without allowing them to lick. Bumblebees extending their proboscis were considered motivated and kept for the experiments. Conditioning trials (Fig. 1Aiii) consisted in 15 s of odourless airflow, followed by 6 s of CS, and 3 s of unconditioned stimulus (US, i.e. 50% sucrose solution applied with a toothpick on the bumblebee's antennae), with 2 s of overlap between CS+US, and 20 s of odourless airflow (Aguiar et al., 2018). The inter-trial interval was 10 min. An air extractor was placed behind the bumblebee to prevent odorant accumulation during CS delivery. Bumblebees extending their proboscis within 3 s of US presentation (i.e. 2 s CS+US and 1 s US) were allowed to lick the toothpick soaked in sucrose (50% w/w). In unrewarded trials (see reversal learning and memory tests, below) no US was applied. We scored a conditioned response if the bumblebee extended its proboscis to the odour delivery before sucrose presentation. Bumblebees that responded to the odour in the first conditioning trial were discarded from the analysis. We used conditioned responses to calculate individual scores for each bumblebee describing its performance during conditioning (i.e. acquisition score) (Monchanin et al., 2021; see details below). Exposed and control bumblebees were always conditioned in parallel.

Absolute learning, STM and LTM

We tested the effects of N. ceranae on the ability of bumblebees to associate an odour with a reward. This form of learning primarily requires peripheral brain centres, (i.e. the antennal lobes) as it can be observed in bees with a non-functional central brain (Giurfa and Sandoz, 2012). We trained bumblebees in a spaced 3-trial absolute conditioning learning protocol (Fig. 1Ai) that was shown to generate robust LTM in bees (Menzel et al., 2001). We used the same rewarded odour (CS+) during training of a given bumblebee, but both nonanal and phenylacetaldehyde were used as a CS+ for different bumblebees. For each of these bees, we calculated an acquisition score (sum of PER responses, i.e. 0–2), and compared the learning curve to assess the increase of PER responses over trials. Bumblebees that did not respond to the US in at least one trial were considered not motivated and were removed from the learning analysis (i.e. 30.3% control and 29.6% exposed bumblebees; Table 1).

We tested memory retention of bumblebee responders that showed at least one conditioned response in either of the last two trials (Simcock et al., 2018; Wright et al., 2015) (i.e. 41.3% of bumblebees conditioned for STM and 30–68% of bumblebees conditioned for LTM; see details in Table 1). We performed tests either 1 h (i.e. STM) or 24 h (i.e. LTM) after the last acquisition trial. Bumblebees were tested for either 1 h or 24 h memory, but never for both because of the risk of reconsolidation or extinction of memory when the CS is presented several times without a reward to bees (Bouton and Moody, 2004). For tests performed after 24 h, bumblebees were fed until satiation with 50% (w/w) sucrose solution right after conditioning; unfortunately, 9–20% of responder bumblebees died before the test (Table 1). In bees, LTM is dependent on protein synthesis whereas STM is not (Menzel, 2001). Studying these two types of memory was thus a means to explore whether exposure to N. ceranae interfered with protein synthesis. We presented bumblebees with the two odorants without any reward: one odour used as a CS+ to test for memory formation, and the second odour as a novel odorant (NOd), to control for potential generalization (Matsumoto et al., 2012). For example, when nonanal was used as the CS+, phenylacetaldehyde was used as a NOd, and vice versa. Bumblebees responding only to the CS+, and not to the NOd, were considered to have generated a specific memory to the rewarded odour. Other responses registered were those to the NOd alone (inverted response), both odours (general response) or neither of the odours (no memory). Just after the memory test, we tested the motivation of bumblebees by touching their antennae with a toothpick soaked with 50% (w/w) sucrose solution. Bumblebees that did not respond to the US were discarded from the analysis (i.e. 14.06% of alive responder bumblebees in LTM; Table 1). Sample sizes dropped between learning and memory tests because of the selection protocol of responders, mortality in tubes and loss of motivation by the bumblebees (see details in Table 1).

Reversal learning

We tested the effects of N. ceranae on the ability of bumblebees to learn to discriminate two odours and reverse the task. This form of learning involves two phases. The differential conditioning phase (phase 1) requires the antennal lobes but is not dependent on central brain centres (i.e. mushroom bodies) (Boitard et al., 2015; Devaud et al., 2007). The reversal learning phase (phase 2) requires both functional antennal lobes and mushroom bodies to be observed (Boitard et al., 2015; Devaud et al., 2007). In phase 1, we trained bumblebees to discriminate between two odours. This consisted in 10 trials (Fig. 1Aii), five with each odour that was either paired with the US (A+) or unpaired (B−), presented in a pseudo-random order. The rewarded and unrewarded odours were randomized on different training days. Bumblebees that did not respond to the US in two or more trials were discarded from the analysis and for the reversal phase. In phase 2, we trained bumblebees to invert the first learnt contingency. This reversal phase started 1 h after the end of the differential phase. Here, we trained bumblebees in 12 trials, 6 with each odour presented in a pseudo-random order. The previously rewarded odour was not associated with a reward anymore (A−) while the previously unrewarded odour became rewarded (B+). To start from the same level of learning, in the first trial of the reversal phase we presented the odour A−, and kept for analysis only those bumblebees that extended their proboscis (Table 1). We analysed the performance of bumblebees in each phase separately by attributing them an acquisition score (sum of all trials where the bee responded to CS+ but not to CS−) for each phase.

Infection status

We assessed the infection status of bumblebees that finished the tests by PCR using 218MITOC primers (Martín-Hernández et al., 2007). These primers are 100% specific for N. ceranae and have a 0% error rate (false negatives) for reactions containing DNA equivalent to 2000 spores (Martín-Hernández et al., 2007). For each bumblebee, the entire gut was extracted, homogenized in sterile dH2O and vortexed with 2 mm glass beads (Labbox Labware). Genomic DNA was extracted using proteinase K (20 mg ml−1; Euromedex) and 1 mmol l−1 of Tris-EDTA buffer (pH 8). A sample with N. ceranae spores was included in each round of extraction as a positive control. PCR was performed with the Taq Polymerase Direct Loading Buffer (5 U μl−1; MP Biomedicals) following the manufacturer's instructions. We used a final volume of 25 μl with 0.4 μmol l−1 of each pair of primers (Martín-Hernández et al., 2007), 200 μmol l−1 of dNTPs (Jena Biosciences), 0.48 μg μl−1 of BSA (Sigma) and 2.5 μl of DNA sample. PCR reactions were carried out in a S1000™ Thermal Cycler (Bio-Rad). Thermal conditions were 94°C for 2 min, 35 cycles of 94°C for 30 s, 61.8°C for 45 s and 72°C for 2 min, and a final step of 72°C for 7 min. The length of the PCR products (i.e. 218 bp) was checked by electrophoresis with a 1.2% agarose gel stained with SYBR Safe DNA Stain (Edvotek, Washington, DC, USA). Positive and negative PCR controls were run in parallel. Based on the PCR results, we classified bumblebees into three different infection statuses: control, Nosema-exposed negative (NE−) or Nosema-exposed positive (NE+). Bumblebees that were not exposed to the parasite but nevertheless showed a positive result in PCR were excluded from the analysis (i.e. 6.26%; 23 out of 367 control bumblebees). These positives may be due to the fact that despite our precautions before starting the experiments (PCR screening of about 10% of the workers in the colonies, use of UV-treated pollens), it is possible that the commercial colonies we used were not fully free of parasites. Additionally, we cannot completely exclude low levels of cross-contamination between treatments during manipulations.

Statistical analysis

All analyses were conducted in RStudio (version 4.1.0). We evaluated the effects of parasite exposure and infection on learning curves, and gustatory, acquisition, learning and memory scores. The proportion of responses to the different sucrose concentrations was analysed in a generalized linear mixed model (GLMM) (package lme4; Bates et al., 2015), with infection status and concentration as fixed factors. Learning curves of absolute and reversal learning experiments were analysed with a binomial GLMM with infection status, trial and the interaction between them as fixed factors. Whenever the interaction was not significant it was removed from the model. The responses to rewarded and unrewarded odours during reversal learning were analysed separately. To determine the ability of bumblebees to learn to differentiate a rewarded from an unrewarded odour (i.e. differential phase) and the opposite (i.e. reversal phase) in the reversal learning experiment, we studied the interaction of infection status, trial and reward in a binomial GLMM, followed by a Tukey pairwise comparison applying the function lsmeans (package emmeans). All these models included colony and bee identity as random factors. Gustatory scores for sucrose responsiveness, and acquisition scores during learning experiments were analysed with a linear model (LM). These models included bee infection status as a fixed factor and colony of origin as a random factor. We performed Tukey post hoc pairwise comparisons (package multcomp; Hothorn et al., 2008) to assess the relationship between the three bee infection statuses. The effect of infection status on the ability of bumblebees to generate specific memory (i.e. response to CS+ only) or no memory (i.e. response to any odour) was compared with a Chi-square test. Raw data are available from Zenodo (doi:10.5281/zenodo.4376362).

Parasite exposure did not influence sucrose responsiveness

We tested responsiveness to different sucrose concentrations in 63 bumblebees showing a consistent response (37 controls, 16 NE−, 10 NE+; Table 1). Bumblebees increased their response to sucrose solution with concentration (Fig. 2A; GLMM, estimate=0.098, s.e.=0.009, P<0.001) in the three infection statuses (GLMM, infection status: χ2=4.754, d.f.=2, P=0.09). Gustatory scores ranged from 1.32±1.22 (mean±s.e.m.) in controls to 1.43±0.24 in NE− and 1.80±0.41 in NE+ bumblebees, but did not differ between infection status (Fig. 2B; LM, F=1.077, P=0.347). Therefore, exposure to N. ceranae affected neither the sucrose sensitivity nor the feeding motivation of bumblebees.

Fig. 2.

Sucrose responsiveness. (A) Proportion of control (blue, n=37), Nosema-exposed negative (NE−; yellow, n=16) and Nosema-exposed positive (NE+; red, n=10) bumblebees responding to an increase in sucrose concentration. Pie chart represents the percentage of NE− and NE+ bees (total n=30). (B) Violin plots showing the gustatory score of bumblebees as the sum of all responses for each bumblebee. Black diamonds represent the mean score for each infection status. White circles represent the score for each individual. n is the sample size. Gustatory scores did not differ between infection statuses (GLMM, P>0.05).

Fig. 2.

Sucrose responsiveness. (A) Proportion of control (blue, n=37), Nosema-exposed negative (NE−; yellow, n=16) and Nosema-exposed positive (NE+; red, n=10) bumblebees responding to an increase in sucrose concentration. Pie chart represents the percentage of NE− and NE+ bees (total n=30). (B) Violin plots showing the gustatory score of bumblebees as the sum of all responses for each bumblebee. Black diamonds represent the mean score for each infection status. White circles represent the score for each individual. n is the sample size. Gustatory scores did not differ between infection statuses (GLMM, P>0.05).

Parasite exposure reduced absolute learning but not memory

We analysed absolute conditioning in 420 bumblebees (141 controls, 228 NE−, 51 NE+; Table 1). The proportion of bumblebees showing a conditioned response to CS+ increased with trials (GLMM; trial: χ2=58.432, d.f.=2, P<0.001), but this trend was lower for NE− (estimate=−1.056, s.e.=0.2950, P=0.00034) and NE+ (estimate=−1.388, s.e.=0.474, P=0.003) bumblebees than for controls (Fig. 3). Likewise, exposed bumblebees (either NE− or NE+) had significantly lower acquisition scores (Fig. 3B; LM, infection status: χ2=15.897, d.f.=2, P<0.001) than controls. Acquisition scores were similar for NE− and NE+ (Tukey P>0.05; Table S1).

Fig. 3.

Absolute learning and memory. (A) Learning curves show the percentage of control (blue, n=141), NE− (yellow, n=228) and NE+ (red, n=51) bumblebees extending their proboscis to the odour during conditioning. Pie chart shows the percentage of NE− and NE+ bumblebees that finished conditioning (total n=279). (B) Violin plots for acquisition score (sum of correct responses). Black diamonds represent the mean score for each infection status. White circles represent the score of each individual. n is the sample size. (C,D) Short-term memory (STM) and long-term memory (LTM) tests. Bar plots show the proportion of bumblebees responding to both the CS+ and a novel odorant (NOd; i.e. general memory), the CS+ only (specific memory), the NOd only (inverted memory) or no odours (no memory), 1 h (STM; C) and 24 h (LTM; D) after training. Numbers inside the bars represent the sample size. Different letters represent significant differences between infection statuses (GLMM; P<0.05).

Fig. 3.

Absolute learning and memory. (A) Learning curves show the percentage of control (blue, n=141), NE− (yellow, n=228) and NE+ (red, n=51) bumblebees extending their proboscis to the odour during conditioning. Pie chart shows the percentage of NE− and NE+ bumblebees that finished conditioning (total n=279). (B) Violin plots for acquisition score (sum of correct responses). Black diamonds represent the mean score for each infection status. White circles represent the score of each individual. n is the sample size. (C,D) Short-term memory (STM) and long-term memory (LTM) tests. Bar plots show the proportion of bumblebees responding to both the CS+ and a novel odorant (NOd; i.e. general memory), the CS+ only (specific memory), the NOd only (inverted memory) or no odours (no memory), 1 h (STM; C) and 24 h (LTM; D) after training. Numbers inside the bars represent the sample size. Different letters represent significant differences between infection statuses (GLMM; P<0.05).

We analysed STM and LTM formation of bumblebee responders (i.e. those that showed at least one conditioned response in trial 2 and/or trial 3: 41.3% for STM and 50% for LTM; Table 1). One hour after training, bumblebees either responded to the CS+ only (i.e. specific memory; 82.02±0.011%, mean±s.e.m.) or did not respond to any odour (i.e. no memory), with no bumblebees responding to NOd (i.e. neither generalizers nor inverters; Fig. 3C). The proportion of bumblebees with specific memory was similar in the three infection statuses (χ2=0.013, d.f.=2, P=0.993), indicating that N. ceranae did not affect STM (Fig. 3C). A much lower proportion of the bumblebees (33.03±6.35%, mean±s.e.m.) showed specific memory for the conditioned odour after 24 h (Fig. 3D). Two bumblebees generalized their responses (one control and one NE−), and only one NE− bumblebee inverted its response (i.e. response to NOd only). The proportion of bumblebees showing specific memory to the CS+ did not differ between infection statuses (Fig. 3D; χ2=2.349, d.f.=2, P=0.310), presumably because of the low number of NE+ individuals (9 bumblebees).

Parasite exposure reduced differential and reversal learning

We analysed differential learning in 125 bumblebees (64 controls, 41 NE−, 20 NE+; Table 1). The proportion of bumblebees responding to A+ was affected by the interaction between the infection status and trial (GLMM; infection status×trial: χ2=7.987, d.f.=2, P=0.018). Responses increased with trials in control and NE− bumblebees (estimate=1.220, s.e.=0.1879, P<0.001), but not in NE+ (estimate=−0.732, s.e.=0.299, P=0.014). Control bumblebees discriminated the two odours (i.e. higher proportion of responses to A+ than to B−) at trial 2 (Tukey; z=3.674, P=0.044) and NE− bumblebees discriminated the two odours at trial 4 (Tukey: z=4.083, P=0.009), but NE+ bumblebees did not show this ability after 5 trials (Tukey: z=2.514, P=0.665) (Fig. 4A). Parasite exposure affected acquisition scores (GLMM, infection status, χ2=29.978, d.f.=2, P<0.001; Fig. 4B). NE− and NE+ bumblebees showed similar acquisition scores (Tukey test: P>0.05; Table S1) and these scores were significantly lower than those of controls (Tukey test: P<0.05; Table S1). Thus overall, exposure to N. ceranae reduced differential learning performance. Exposed bumblebees were slower (i.e. NE−) or unable (i.e. NE+) to solve the task.

Fig. 4.

Reversal learning. (A,B) Differential learning phase. (A) Percentage of PER responses to rewarded (A+, circles) and unrewarded (B−, triangles) odours by control (blue, n=64), NE− (yellow, n=41) and NE+ (red, n=20) bumblebees. Pie chart shows the percentage of NE− and NE+ bumblebees that finished conditioning (total n=61). (B) Violin plots of acquisition scores (i.e. sum of correct responses). Black diamonds represent the mean score of each infection status. White circles are the scores for each individual. n is the sample size. (C,D) Reversal learning phase. (C) Curves show the increase in the percentage of PER responses to B+ over A− for trials in control (n=56), NE− (n=32) and NE+ (n=15) bumblebees. (D) Violin plots of acquisition scores. Different letters above violin plots represent significant differences between infection statuses (GLMM, P<0.05).

Fig. 4.

Reversal learning. (A,B) Differential learning phase. (A) Percentage of PER responses to rewarded (A+, circles) and unrewarded (B−, triangles) odours by control (blue, n=64), NE− (yellow, n=41) and NE+ (red, n=20) bumblebees. Pie chart shows the percentage of NE− and NE+ bumblebees that finished conditioning (total n=61). (B) Violin plots of acquisition scores (i.e. sum of correct responses). Black diamonds represent the mean score of each infection status. White circles are the scores for each individual. n is the sample size. (C,D) Reversal learning phase. (C) Curves show the increase in the percentage of PER responses to B+ over A− for trials in control (n=56), NE− (n=32) and NE+ (n=15) bumblebees. (D) Violin plots of acquisition scores. Different letters above violin plots represent significant differences between infection statuses (GLMM, P<0.05).

We analysed reversal learning in 103 bumblebees that finished the differential phase, and that responded to A− in the first trial of the reversal phase (56 controls, 32 NE−, 15 NE+; Table 1). All bumblebees reduced their response to A− over trials (GLMM, trial: estimate=−1.187, s.e.=0.134; P<0.001) and increased their response to B+ (GLMM, trial; estimate=1.215, s.e.=0.124; P<0.001) (Fig. 4D). Infection status did not affect the response to A− (GLMM; infection status: χ2=9, d.f.=2, P=0.389). The proportion of bumblebees responding to B+ was not different between controls and NE− (estimate=−0.579, s.e.=0.348, P=0.096). However, it was significantly reduced in NE+ (estimate=−1.597, s.e.=0.450, P<0.001). Control bumblebees reversed their response to odours earlier, at trial 5 (i.e. higher proportion of bees responding to B+ than A−; Tukey: z=4.478, P=0.002), while NE− did so at trial 6 (Tukey: z=4.281, P=0.006), and NE+ never did so (Tukey; trial 6: z=2.103, P=0.960). This was also reflected in the acquisition scores, which were significantly different between bumblebees of different infection statuses (Fig. 4E; GLMM, infection status, χ2=6.783, P=0.033). NE+ bumblebees had lower acquisition scores than controls (Tukey test: P=0.053; Table S1) and NE− bumblebees (Tukey test: P=0.044; Table S1), suggesting they had a lower ability to reverse the task. Thus, overall, exposure to N. ceranae also impaired the reversal phase of reversal learning. We found no effect of N. ceranae exposure on either phase of reversal learning when bumblebees were tested at 2 days post-exposure, suggesting that stress due to parasite exposure or parasite infection requires a longer time to be established (see SupplementaryMaterials and Methods, Table S2 and Fig. S1).

Bees are exposed to a number of parasites that can affect cognitive ability supporting crucial behaviour (Koch et al., 2017; Schmid-Hempel, 2013). Previous studies exploring the effect of N. ceranae on absolute olfactory learning and memory in bees reported contrasting results, presumably because of differences in conditioning protocols and infection rates across studies (Bell et al., 2020; Charbonneau et al., 2016; Gage et al., 2018; Piiroinen and Goulson, 2016; Piiroinen et al., 2016). Here, we ran a suite of standard olfactory cognitive assays showing that feeding bumblebees spores of this parasite consistently impaired different types of olfactory learning but not memory 7 days after exposure.

Exposure to N. ceranae in food impaired the ability of bumblebees to associate an odour with a reward (absolute learning), discriminate two odours (differential learning) and learn an opposite association (reversal learning). These are fundamental cognitive operations a bee must display to efficiently forage on flowers (Giurfa and Sandoz, 2012). Our results agree with a previous study reporting reduced absolute learning in N. ceranae exposed bumblebees (Piiroinen and Goulson, 2016). We also found that N. ceranae did not affect sucrose responsiveness, contrary to observations in honey bees in which it was found to increase their hunger (Mayack and Naug, 2009). The parasite may thus not produce the same energetic stress observed previously in honey bees (Mayack and Naug, 2009), where N. ceranae seems to be a more specific parasite (van der Steen et al., 2022). However, further experiments are needed to confirm this as protocols in these previous studies differed slightly from ours. For instance, under our conditions we cannot discard the possibility that response to high sucrose concentrations was due to sensitization as no water was used between sucrose concentrations. In bees, sucrose perception through the antennae and olfactory learning require processing of olfactory information through the antennal lobes (Giurfa and Sandoz, 2012). While simple forms of associations can be acquired only with functional antennal lobes, other types of learning (i.e. reversal learning, configural learning) also require information processing in the mushroom bodies (Devaud et al., 2007, 2015; Boitard et al., 2015). In our experiments, the fact that all types of learning were impaired and that sucrose sensitivity was not suggests that N. ceranae did not specifically target the antennal lobes or the mushroom bodies. Rather, it probably impacted the learning processes in general.

By contrast, we found no evidence that N. ceranae influenced memory or sucrose sensitivity. During training, animals learn and form short-term memories that are later consolidated and transformed into stable long-term memories (Menzel and Müller, 1996) after protein synthesis (Menzel, 2001). In our experiments, N. ceranae impaired neither STM nor LTM.

It has recently been questioned whether bumblebees are natural hosts of N. ceranae based on the lack of evidence of parasitic forms inside host cells (Gisder et al., 2020). Several studies have nevertheless reported N. ceranae in wild bumblebees at low (e.g. 4.76%; Sinpoo et al., 2018) and high prevalence (e.g. 72%; Arbulo et al., 2015). Whether or not bumblebees are suitable hosts for N. ceranae replication, our results imply they are impacted by an acute exposure to the parasite. Such exposure may be extremely frequent in nature because of the high prevalence of N. ceranae in honey bees (Runckel et al., 2011) that contaminate flowers with spores through physical contact or in their faeces (Graystock et al., 2015). Our protocol of parasite exposure significantly increased the infection rate of bumblebees to 28% in comparison to previous studies (Piiroinen and Goulson, 2016; Piiroinen et al., 2016), which allowed the evaluation of cognitive traits in bumblebees in the three infection statuses. Bumblebees that tested positive to N. ceranae showed a tendency for lower cognitive performance than those that were exposed but tested negative. They reached the lowest learning during absolute conditioning and did not discriminate odours, suggesting that infection may interfere with some aspects of cognition. However further experiments are needed to tackle this question, as the lower performance of exposed bumblebees positive to N. ceranae could also be related to a worse health status and therefore higher susceptibility to become infected.

Through which mechanism does the parasite impair learning? Ingestion of N. ceranae spores exerts a stress that can reduce cognition. Nosema ceranae is known to alter the immune system of bees; for example, by modulating the expression of antimicrobial peptides (Antúnez et al., 2009; Botías et al., 2021; Sinpoo et al., 2018). Stimulation of the immune system with non-pathogenic elicitors, such as lipopolysaccharides (LPS), was shown to reduce learning ability in honey bees, which were less able to associate an odour with a reward (Mallon et al., 2003), and bumblebees, which showed a lower performance in odour (Mobley and Gegear, 2018) and colour differential learning tasks (Mobley and Gegear, 2018). It is thus possible that the observed effects of N. ceranae exposure on bumblebee cognition were caused by activation of the immune response. Exposed positive bumblebees showed lower learning performance than exposed negative bumblebees during differential and reversal learning tasks, suggesting a further effect of parasite infection, rather than just exposure. In honey bees, infection with N. ceranae was shown to downregulate the expression of genes in the brain (Doublet et al., 2016), some of which are linked to olfaction (Badaoui et al., 2017; Doublet et al., 2016), potentially leading to changes in behaviour and cognition. Whether N. ceranae downregulates gene expression in the brain of bumblebees needs to be addressed. Interestingly, none of these effects were observed at 2 days post-exposure (see Supplementary Materials and Methods, Table S2 and Fig. S1). So far it is unknown whether N. ceranae triggers the immune response at this time. In honey bees, the earliest effects have been observed after 3 days (Chaimanee et al., 2012).

As olfactory learning is essential for foraging, this sublethal effect of N. ceranae exposure on bumblebee cognition can compromise colony foraging success, as well as chemical communication between bees, ultimately leading to colony collapse. Parasite load in the field can range from a few to thousands spores (Meana et al., 2010). Here, we used a substantially higher spore load of N. ceranae to infect commercially reared bumblebees than the actual infection rates found in wild bumblebees (e.g. 6800 spores per individual; Graystock et al., 2013). Commercial and wild bumblebee colonies exhibit physiological and behavioural differences as a result of different selective pressures (Velthuis and van Doorn, 2006), and may, therefore, show different susceptibility to parasites. Therefore, further studies are needed to analyse the effects of different concentrations of N. ceranae spores and their possible interactions with other stressors in the field in order to assess their real impact on wild pollinators. Beyond bees, these effects may also have broader fundamental consequences for plants and parasites. From the plant perspective, an impaired flower constancy by pollinators may increase pollen transfer between incompatible flowers of different species, and therefore reduce pollination efficiency. From the parasite perspective, foraging errors due to impaired learning by bees may decrease they tendency for flower constancy. This would be beneficial for the parasite as it may favour its spread across flower species and thus possibly increase the range of hosts.

We thank Philipp Heeb, Cristian Pasquaretta and Maria Eugenia Villar Damiani for constructive discussions, Charline Schartier for helping running some experiments, and Martin Giurfa for lending the automated PER conditioning setup.

Author contributions

Conceptualization: T.G.-M., M.L.; Methodology: T.G.-M., T.D.; Validation: T.G.-M., M.L.; Formal analysis: T.G.-M.; Investigation: T.G.-M., M.L.; Data curation: T.G.-M., T.D.; Writing - original draft: T.G.-M., T.D., M.L.; Writing - review & editing: T.G.-M., T.D., M.L.; Funding acquisition: M.L.

Funding

T.G.-M. was funded by a postdoctoral grant of the Fondation Fyssen. During writing, T.G.-M., T.D. and M.L. were funded by grants from the Agence Nationale de la Recherche (POLLINET ANR-16-CE02-0002-01, 3DNaviBee ANR-19-CE37-0024, BEE-MOVE ANR-20-ERC8-0004-01), the European Regional Development Fund (ECONECT), and the Agence de la Transition Ecologique (LOTAPIS).

Data availability

Raw data are available from Zenodo (doi:10.5281/zenodo.4376362).

Aguiar
,
J. M. R. B. V.
,
Roselino
,
A. C.
,
Sazima
,
M.
and
Giurfa
,
M.
(
2018
).
Can honey bees discriminate between floral-fragrance isomers?
J. Exp. Biol.
221
,
jeb180844
.
Antúnez
,
K.
,
Martín-Hernández
,
R.
,
Prieto
,
L.
,
Meana
,
A.
,
Zunino
,
P.
and
Higes
,
M.
(
2009
).
Immune suppression in the honey bee (Apis mellifera) following infection by Nosema ceranae (Microsporidia)
.
Environ. Microbiol.
11
,
2284
-
2290
.
Arbulo
,
N.
,
Antúnez
,
K.
,
Salvarrey
,
S.
,
Santos
,
E.
,
Branchiccela
,
B.
,
Martín-Hernández
,
R.
,
Higes
,
M.
and
Invernizzi
,
C.
(
2015
).
High prevalence and infection levels of Nosema ceranae in bumblebees Bombus atratus and Bombus bellicosus from Uruguay
.
J. Invertebr. Pathol.
130
,
165
-
168
.
Badaoui
,
B.
,
Fougeroux
,
A.
,
Petit
,
F.
,
Anselmo
,
A.
,
Gorni
,
C.
,
Cucurachi
,
M.
,
Cersini
,
A.
,
Granato
,
A.
,
Cardeti
,
G.
,
Formato
,
G.
et al. 
(
2017
).
RNA-sequence analysis of gene expression from honeybees (Apis mellifera) infected with Nosema ceranae
.
PLoS One
12
,
e0173438
.
Baracchi
,
D.
,
Rigosi
,
E.
,
de Brito Sanchez
,
G.
,
Giurfa
,
M.
(
2018
).
Lateralization of sucrose responsiveness and non-associative learning in honeybees
.
Front. Psychol.
9
,
425
.
Bates
,
D.
,
Mächler
,
M.
,
Bolker
,
B.
and
Walker
,
S.
(
2015
).
Fitting linear mixed-effects models using lme4
.
J. Stat. Softw.
67
,
1
-
48
.
Bell
,
H. C.
,
Montgomery
,
C. N.
,
Benavides
,
J. E.
and
Nieh
,
J. C.
(
2020
).
Effects of Nosema ceranae (Dissociodihaplophasida: Nosematidae) and flupyradifurone on olfactory learning in honey bees, Apis mellifera (Hymenoptera: Apidae)
.
J. Insect Sci.
20
,
29
.
Boitard
,
C.
,
Devaud
,
J.-M.
,
Isabel
,
G.
and
Giurfa
,
M.
(
2015
).
GABAergic feedback signaling into the calyces of the mushroom bodies enables olfactory reversal learning in honey bees
.
Front. Behav. Neurosci.
9
,
198
.
Botías
,
C.
,
Jones
,
J. C.
,
Pamminger
,
T.
,
Bartomeus
,
I.
,
Hughes
,
W. O. H.
and
Goulson
,
D.
(
2021
).
Multiple stressors interact to impair the performance of bumblebee Bombus terrestris colonies
.
J. Anim. Ecol.
90
,
415
-
431
.
Bouton
,
M. E.
and
Moody
,
E. W.
(
2004
).
Memory processes in classical conditioning
.
Neurosci. Biobehav. Rev.
28
,
663
-
674
.
Cantwell
,
G.
(
1970
).
Standard methods for counting Nosema spores
.
Am. Bee J.
110
,
222
-
223
.
Chaimanee
,
V.
,
Chantawannakul
,
P.
,
Chen
,
Y.
,
Evans
,
J. D.
and
Pettis
,
J. S.
(
2012
).
Differential expression of immune genes of adult honey bee (Apis mellifera) after inoculated by Nosema ceranae
.
J. Insect Physiol.
58
,
1090
-
1095
.
Charbonneau
,
L. R.
,
Hillier
,
N. K.
,
Rogers
,
R. E. L.
,
Williams
,
G. R.
and
Shutler
,
D.
(
2016
).
Effects of Nosema apis, N. ceranae, and coinfections on honey bee (Apis mellifera) learning and memory
.
Sci. Rep.
6
,
22626
.
Collett
,
M.
,
Chittka
,
L.
and
Collett
,
T. S.
(
2013
).
Spatial memory in insect navigation
.
Curr. Biol.
23
,
R789
-
R800
.
Cox-Foster
,
D. L.
,
Conlan
,
S.
,
Holmes
,
E. C.
,
Palacios
,
G.
,
Evans
,
J. D.
,
Moran
,
N. A.
,
Quan
,
P.-L.
,
Briese
,
T.
,
Hornig
,
M.
,
Geiser
,
D. M.
et al. 
(
2007
).
A metagenomic survey of microbes in honey bee Colony Collapse Disorder
.
Science
318
,
283
-
287
.
Devaud
,
J.-M.
,
Blunk
,
A.
,
Podufall
,
J.
,
Giurfa
,
M.
and
Grünewald
,
B.
(
2007
).
Using local anaesthetics to block neuronal activity and map specific learning tasks to the mushroom bodies of an insect brain: local anaesthetics and odour learning in honeybees
.
Eur. J. Neurosci.
26
,
3193
-
3206
.
Devaud
,
J.-M.
,
Papouin
,
T.
,
Carcaud
,
J.
,
Sandoz
,
J.-C.
,
Grünewald
,
B.
and
Giurfa
,
M.
(
2015
).
Neural substrate for higher-order learning in an insect: mushroom bodies are necessary for configural discriminations
.
Proc. Natl. Acad. Sci. USA
112
,
E5854
-
E5862
.
Doublet
,
V.
,
Paxton
,
R. J.
,
McDonnell
,
C. M.
,
Dubois
,
E.
,
Nidelet
,
S.
,
Moritz
,
R. F. A.
,
Alaux
,
C.
and
Le Conte
,
Y.
(
2016
).
Brain transcriptomes of honey bees (Apis mellifera) experimentally infected by two pathogens: black queen cell virus and Nosema ceranae
.
Genomics Data
10
,
79
-
82
.
Dussaubat
,
C.
,
Maisonnasse
,
A.
,
Crauser
,
D.
,
Beslay
,
D.
,
Costagliola
,
G.
,
Soubeyrand
,
S.
,
Kretzchmar
,
A.
and
Le Conte
,
Y.
(
2013
).
Flight behavior and pheromone changes associated to Nosema ceranae infection of honey bee workers (Apis mellifera) in field conditions
.
J. Invertebr. Pathol.
113
,
42
-
51
.
Fries
,
I.
,
Chauzat
,
M.-P.
,
Chen
,
Y.-P.
,
Doublet
,
V.
,
Genersch
,
E.
,
Gisder
,
S.
,
Higes
,
M.
,
McMahon
,
D. P.
,
Martín-Hernández
,
R.
,
Natsopoulou
,
M.
et al. 
(
2013
).
Standard methods for Nosema research
.
J. Apic. Res.
52
,
1
-
28
.
Frost
,
P. C.
,
Ebert
,
D.
and
Smith
,
V. H.
(
2008
).
Bacterial infection changes the elemental composition of Daphnia magna
.
J. Anim. Ecol.
77
,
1265
-
1272
.
Gage
,
S. L.
,
Kramer
,
C.
,
Calle
,
S.
,
Carroll
,
M.
,
Heien
,
M.
and
DeGrandi-Hoffman
,
G.
(
2018
).
Nosema ceranae parasitism impacts olfactory learning and memory and neurochemistry in honey bees (Apis mellifera)
.
J. Exp. Biol.
221
,
jeb161489
.
Gisder
,
S.
,
Horchler
,
L.
,
Pieper
,
F.
,
Schüler
,
V.
,
Šima
,
P.
and
Genersch
,
E.
(
2020
).
Rapid gastrointestinal passage may protect Bombus terrestris from becoming a true host for Nosema ceranae
.
Appl. Environ. Microbiol.
86
,
e00629-20
.
Giurfa
,
M.
(
2015
).
Learning and cognition in insects: learning and insect cognition
.
Wiley Interdiscip. Rev. Cogn. Sci.
6
,
383
-
395
.
Giurfa
,
M.
and
Sandoz
,
J.-C.
(
2012
).
Invertebrate learning and memory: fifty years of olfactory conditioning of the proboscis extension response in honeybees
.
Learn. Mem.
19
,
54
-
66
.
Gómez-Moracho
,
T.
,
Heeb
,
P.
and
Lihoreau
,
M.
(
2017
).
Effects of parasites and pathogens on bee cognition: bee parasites, pathogens and cognition
.
Ecol. Entomol.
42
,
51
-
64
.
Gómez-Moracho
,
T.
,
Durand
,
T.
,
Pasquaretta
,
C.
,
Heeb
,
P.
and
Lihoreau
,
M.
(
2021
).
Artificial diets modulate infection rates by Nosema ceranae in bumblebees
.
Microorganisms
9
,
158
.
Graystock
,
P.
,
Yates
,
K.
,
Darvill
,
B.
,
Goulson
,
D.
and
Hughes
,
W. O. H.
(
2013
).
Emerging dangers: deadly effects of an emergent parasite in a new pollinator host
.
J. Invertebr. Pathol.
114
,
114
-
119
.
Graystock
,
P.
,
Goulson
,
D.
and
Hughes
,
W.
(
2015
).
Parasites in bloom: flowers aid dispersal and transmission of pollinator parasites within and between bee species
.
Proc. R. Soc. B Biol. Sci.
282
,
20151371
.
Henry
,
M.
,
Beguin
,
M.
,
Requier
,
F.
,
Rollin
,
O.
,
Odoux
,
J.-F.
,
Aupinel
,
P.
,
Aptel
,
J.
,
Tchamitchian
,
S.
and
Decourtye
,
A.
(
2012
).
A common pesticide decreases foraging success and survival in honey bees
.
Science
336
,
348
-
350
.
Higes
,
M.
,
Martín
,
R.
and
Meana
,
A.
(
2006
).
Nosema ceranae, a new microsporidian parasite in honeybees in Europe
.
J. Invertebr. Pathol.
92
,
93
-
95
.
Higes
,
M.
,
García-Palencia
,
P.
,
Martín-Hernández
,
R.
and
Meana
,
A.
(
2007
).
Experimental infection of Apis mellifera honeybees with Nosema ceranae (Microsporidia)
.
J. Invertebr. Pathol.
94
,
211
-
217
.
Higes
,
M.
,
Martín-Hernández
,
R.
,
Botías
,
C.
,
Bailón
,
E. G.
,
González-Porto
,
A. V.
,
Barrios
,
L.
,
del Nozal
,
M. J.
,
Bernal
,
J. L.
,
Jiménez
,
J. J.
,
Palencia
,
P. G.
et al. 
(
2008
).
How natural infection by Nosema ceranae causes honeybee colony collapse
.
Environ. Microbiol.
10
,
2659
-
2669
.
Holt
,
H. L.
,
Aronstein
,
K. A.
and
Grozinger
,
C. M.
(
2013
).
Chronic parasitization by Nosema microsporidia causes global expression changes in core nutritional, metabolic and behavioral pathways in honey bee workers (Apis mellifera)
.
BMC Genom.
14
,
799
.
Hothorn
,
T.
,
Bretz
,
F.
and
Westfall
,
P.
(
2008
).
Simultaneous inference in general parametric models
.
Biom. J.
50
,
346
-
363
.
Klee
,
J.
,
Tek Tay
,
W.
and
Paxton
,
R. J.
(
2006
).
Specific and sensitive detection of Nosema bombi (Microsporidia: Nosematidae) in bumble bees (Bombus spp.; Hymenoptera: Apidae) by PCR of partial rRNA gene sequences
.
J. Invertebr. Pathol.
91
,
98
-
104
.
Klee
,
J.
,
Besana
,
A. M.
,
Genersch
,
E.
,
Gisder
,
S.
,
Nanetti
,
A.
,
Tam
,
D. Q.
,
Chinh
,
T. X.
,
Puerta
,
F.
,
Ruz
,
J. M.
,
Kryger
,
P.
et al. 
(
2007
).
Widespread dispersal of the microsporidian Nosema ceranae, an emergent pathogen of the western honey bee, Apis mellifera
.
J. Invertebr. Pathol.
96
,
1
-
10
.
Klein
,
S.
,
Cabirol
,
A.
,
Devaud
,
J.-M.
,
Barron
,
A. B.
and
Lihoreau
,
M.
(
2017
).
Why bees are so vulnerable to environmental stressors
.
Trends Ecol. Evol.
32
,
268
-
278
.
Koch
,
H.
,
Brown
,
M. J.
and
Stevenson
,
P. C.
(
2017
).
The role of disease in bee foraging ecology
.
Curr. Opin. Insect Sci.
21
,
60
-
67
.
Kraus
,
S.
,
Gómez-Moracho
,
T.
,
Pasquaretta
,
C.
,
Latil
,
G.
,
Dussutour
,
A.
and
Lihoreau
,
M.
(
2019
).
Bumblebees adjust protein and lipid collection rules to the presence of brood
.
Curr. Zool.
65
,
437
-
446
.
Laloi
,
D.
,
Sandoz
,
J. c.
,
Picard-Nizou
,
A. l.
,
Marchesi
,
A.
,
Pouvreau
,
A.
,
Taséi
,
J. n.
,
Poppy
,
G.
and
Pham-delègue
,
M. h.
(
1999
).
Olfactory conditioning of the proboscis extension in bumble bees
.
Entomol. Exp. Appl.
90
,
123
-
129
.
Lavond
,
D. G.
and
Steinmetz
,
J. E.
(
2003
).
Handbook of Classical Conditioning
.
Boston, MA
:
Springer US
.
Mallon
,
E. B.
,
Brockmann
,
A.
and
Schmid-Hempel
,
P.
(
2003
).
Immune response inhibits associative learning in insects
.
Proc. R. Soc. Lond. B Biol. Sci.
270
,
2471
-
2473
.
Martín-Hernández
,
R.
,
Meana
,
A.
,
Prieto
,
L.
,
Salvador
,
A. M.
,
Garrido-Bailón
,
E.
and
Higes
,
M.
(
2007
).
Outcome of colonization of Apis mellifera by Nosema ceranae
.
Appl. Environ. Microbiol.
73
,
6331
-
6338
.
Matsumoto
,
Y.
,
Menzel
,
R.
,
Sandoz
,
J.-C.
and
Giurfa
,
M.
(
2012
).
Revisiting olfactory classical conditioning of the proboscis extension response in honey bees: a step toward standardized procedures
.
J. Neurosci. Methods
211
,
159
-
167
.
Mayack
,
C.
and
Naug
,
D.
(
2009
).
Energetic stress in the honeybee Apis mellifera from Nosema ceranae infection
.
J. Invertebr. Pathol.
100
,
185
-
188
.
Meana
,
A.
,
Martín-Hernández
,
R.
and
Higes
,
M.
(
2010
).
The reliability of spore counts to diagnose Nosema ceranae infections in honey bees
.
J. Apic. Res.
49
,
212
-
214
.
Menzel
,
R.
(
2001
).
Searching for the memory trace in a mini-brain, the honeybee
.
Learn. Mem.
8
,
53
-
62
.
Menzel
,
R.
(
2012
).
The honeybee as a model for understanding the basis of cognition
.
Nat. Rev. Neurosci.
13
,
758
-
768
.
Menzel
,
R.
and
Müller
,
U.
(
1996
).
Learning and memory in honeybees: from behavior to neural substrates
.
Annu. Rev. Neurosci.
19
,
379
-
404
.
Menzel
,
R.
,
Manz
,
G.
,
Menzel
,
R.
and
Greggers
,
U.
(
2001
).
Massed and spaced learning in honeybees: the role of cs, us, the intertrial interval, and the test interval
.
Learn. Mem.
8
,
198
-
208
.
Mobley
,
M. W.
and
Gegear
,
R. J.
(
2018
).
Immune-cognitive system connectivity reduces bumblebee foraging success in complex multisensory floral environments
.
Sci. Rep.
8
,
5953
.
Monchanin
,
C.
,
Drujont
,
E.
,
Devaud
,
J.-M.
,
Lihoreau
,
M.
and
Barron
,
A. B.
(
2021
).
Metal pollutants have additive negative effects on honey bee cognition
.
J. Exp. Biol.
224
, jeb
241869
.
Palottini
,
F.
,
Estravis Barcala
,
M. C.
and
Farina
,
W. M.
(
2018
).
Odor learning and its experience-dependent modulation in the South American native bumblebee Bombus atratus (hymenoptera: apidae)
.
Front. Psychol.
9
,
603
.
Perry
,
C. J.
,
Søvik
,
E.
,
Myerscough
,
M. R.
and
Barron
,
A. B.
(
2015
).
Rapid behavioral maturation accelerates failure of stressed honey bee colonies
.
Proc. Natl. Acad. Sci USA
112
,
3427
-
3432
.
Piiroinen
,
S.
and
Goulson
,
D.
(
2016
).
Chronic neonicotinoid pesticide exposure and parasite stress differentially affects learning in honeybees and bumblebees
.
Proc. R. Soc. B Biol. Sci.
283
,
20160246
.
Piiroinen
,
S.
,
Botías
,
C.
,
Nicholls
,
E.
and
Goulson
,
D.
(
2016
).
No effect of low-level chronic neonicotinoid exposure on bumblebee learning and fecundity
.
PeerJ
4
,
e1808
.
Plischuk
,
S.
,
Martín-Hernández
,
R.
,
Prieto
,
L.
,
Lucía
,
M.
,
Botías
,
C.
,
Meana
,
A.
,
Abrahamovich
,
A. H.
,
Lange
,
C.
and
Higes
,
M.
(
2009
).
South American native bumblebees (Hymenoptera: Apidae) infected by Nosema ceranae (Microsporidia), an emerging pathogen of honeybees (Apis mellifera)
.
Environ. Microbiol. Rep.
1
,
131
-
135
.
Porrini
,
M. P.
,
Porrini
,
L. P.
,
Garrido
,
P. M.
,
de Melo E Silva Neto
,
C.
,
Porrini
,
D. P.
,
Muller
,
F.
,
Nuñez
,
L. A.
,
Alvarez
,
L.
,
Iriarte
,
P. F.
and
Eguaras
,
M. J.
(
2017
).
Nosema ceranae in South American native stingless bees and social wasp
.
Microb. Ecol.
74
,
761
-
764
.
Raiser
,
G.
,
Galizia
,
C. G.
and
Szyszka
,
P.
(
2017
).
A high-bandwidth dual-channel olfactory stimulator for studying temporal sensitivity of olfactory processing
.
Chem. Senses
42
,
141
-
151
.
Ravoet
,
J.
,
De Smet
,
L.
,
Meeus
,
I.
,
Smagghe
,
G.
,
Wenseleers
,
T.
and
de Graaf
,
D. C.
(
2014
).
Widespread occurrence of honey bee pathogens in solitary bees
.
J. Invertebr. Pathol.
122
,
55
-
58
.
Runckel
,
C.
,
Flenniken
,
M. L.
,
Engel
,
J. C.
,
Ruby
,
J. G.
,
Ganem
,
D.
,
Andino
,
R.
and
DeRisi
,
J. L.
(
2011
).
Temporal analysis of the honey bee microbiome reveals four novel viruses and seasonal prevalence of known viruses, Nosema, and Crithidia
.
PLoS One
6
,
e20656
.
Scheiner
,
R.
,
Barnert
,
M.
and
Erber
,
J.
(
2003
).
Variation in water and sucrose responsiveness during the foraging season affects proboscis extension learning in honey bees
.
Apidologie
34
,
67
-
72
.
Schmid-Hempel
,
P.
(
2013
).
Evolutionary Parasitology: The Integrated Study of Infections, Immunology, Ecology, and Genetics
.
Oxford University Press
.
Schmid-Hempel
,
R.
and
Tognazzo
,
M.
(
2010
).
Molecular divergence defines two distinct lineages of Crithidia bombi (Trypanosomatidae), parasites of bumblebees
.
J. Eukaryot. Microbiol.
57
,
337
-
345
.
Simcock
,
N. K.
,
Gray
,
H.
,
Bouchebti
,
S.
and
Wright
,
G. A.
(
2018
).
Appetitive olfactory learning and memory in the honeybee depend on sugar reward identity
.
J. Insect Physiol.
106
,
71
-
77
.
Sinpoo
,
C.
,
Paxton
,
R. J.
,
Disayathanoowat
,
T.
,
Krongdang
,
S.
and
Chantawannakul
,
P.
(
2018
).
Impact of Nosema ceranae and Nosema apis on individual worker bees of the two host species (Apis cerana and Apis mellifera) and regulation of host immune response
.
J. Insect Physiol.
105
,
1
-
8
.
Smith
,
M. L.
(
2012
).
The honey bee parasite Nosema ceranae: transmissible via food exchange?
PLoS One
7
,
e43319
.
Sommerlandt
,
F. M. J.
,
Rössler
,
W.
and
Spaethe
,
J.
(
2014
).
Elemental and non-elemental olfactory learning using PER conditioning in the bumblebee, Bombus terrestris
.
Apidologie
45
,
106
-
115
.
Takeda
,
K.
(
1961
).
Classical conditioned response in the honey bee
.
J. Insect Physiol.
6
,
168
-
179
.
Toda
,
N. R. T.
,
Song
,
J.
and
Nieh
,
J. C.
(
2009
).
Bumblebees exhibit the memory spacing effect
.
Naturwissenschaften
96
,
1185
-
1191
.
van der Steen
,
J. J. M.
,
Hendriks
,
M. J. A.
,
van Diepeningen
,
A. D.
,
van Gent-Pelzer
,
M. P. E.
and
van der Lee
,
T. A. J.
(
2022
).
Live and dead qPCR detection demonstrates that feeding of Nosema ceranae results in infection in the honey bee but not the bumble bee
.
J. Apic. Res.
61
,
352
-
364
.
Velthuis
,
H. H. W.
and
van Doorn
,
A.
(
2006
).
A century of advances in bumblebee domestication and the economic and environmental aspects of its commercialization for pollination
.
Apidologie
37
,
421
-
451
.
Wolf
,
S.
,
McMahon
,
D. P.
,
Lim
,
K. S.
,
Pull
,
C. D.
,
Clark
,
S. J.
,
Paxton
,
R. J.
and
Osborne
,
J. L.
(
2014
).
So near and yet so far: harmonic radar reveals reduced homing ability of Nosema infected honeybees
.
PLoS One
9
,
e103989
.
Wright
,
G. A.
,
Softley
,
S.
and
Earnshaw
,
H.
(
2015
).
Low doses of neonicotinoid pesticides in food rewards impair short-term olfactory memory in foraging-age honeybees
.
Sci. Rep.
5
,
15322
.

Competing interests

The authors declare no competing or financial interests.

Supplementary information