Crop Pollination Exposes Honey Bees to Pesticides Which Alters Their Susceptibility to the Gut Pathogen Nosema ceranae by Jeffery S. Pettis, Elinor M. Lichtenberg, Michael Andree, Jennie Stitzinger, Robyn Rose, Dennis van Engelsdorp . Published: July 24, 2013DOI: 10.1371/journal.pone.0070182
Recent declines in honey bee populations and increasing demand for insect-pollinated crops raise concerns about pollinator shortages. Pesticide exposure and pathogens may interact to have strong negative effects on managed honey bee colonies. Such findings are of great concern given the large numbers and high levels of pesticides found in honey bee colonies. Thus it is crucial to determine how field-relevant combinations and loads of pesticides affect bee health. We collected pollen from bee hives in seven major crops to determine 1) what types of pesticides bees are exposed to when rented for pollination of various crops and 2) how field-relevant pesticide blends affect bees’ susceptibility to the gut parasite Nosema ceranae. Our samples represent pollen collected by foragers for use by the colony, and do not necessarily indicate foragers’ roles as pollinators. In blueberry, cranberry, cucumber, pumpkin and watermelon bees collected pollen almost exclusively from weeds and wildflowers during our sampling. Thus more attention must be paid to how honey bees are exposed to pesticides outside of the field in which they are placed. We detected 35 different pesticides in the sampled pollen, and found high fungicide loads. The insecticides esfenvalerate and phosmet were at a concentration higher than their median lethal dose in at least one pollen sample. While fungicides are typically seen as fairly safe for honey bees, we found an increased probability of Nosema infection in bees that consumed pollen with a higher fungicide load. Our results highlight a need for research on sub-lethal effects of fungicides and other chemicals that bees placed in an agricultural setting are exposed to.
Citation: Pettis JS, Lichtenberg EM, Andree M, Stitzinger J, Rose R, et al. (2013) Crop Pollination Exposes Honey Bees to Pesticides Which Alters Their Susceptibility to the Gut Pathogen Nosema ceranae. PLoS ONE 8(7): e70182. doi:10.1371/journal.pone.0070182
Editor: Fabio S. Nascimento, Universidade de São Paulo, Faculdade de Filosofia Ciências e Letras de Ribeirão Preto, Brazil
Received: March 25, 2013; Accepted: June 16, 2013; Published: July 24, 2013
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Funding: Funding for this study was provided by the National Honey Board (http://www.honey.com/) and the USDA-ARS Areawide Project on Bee Health (http://www.ars.usda.gov/research/projects/projects.htm?accn_no=412796). Neither the Honey Board nor USDA-ARS Program Staff had a role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: Dennis vanEngesldorp is a PLOS ONE Editor. All other authors have declared that no competing interests exist. This does not alter the authors’ adherence to all the PLOS ONE policies on sharing data and materials.
Honey bees, Apis mellifera, are one of the most important pollinators of agricultural crops . Recent declines in honey bee populations in many North American and European countries – and increasing cultivation of crops that require insects for pollination  raise concerns about pollinator shortages , . Habitat destruction, pesticide use, pathogens and climate change are thought to have contributed to these losses , , . Recent research suggests that honey bee diets, parasites, diseases and pesticides interact to have stronger negative effects on managed honey bee colonies , . Nutritional limitation ,  and exposure to sub-lethal doses of pesticides –, in particular, may alter susceptibility to or severity of diverse bee parasites and pathogens.
Recent research is uncovering diverse sub-lethal effects of pesticides on bees. Insecticides and fungicides can alter insect and spider enzyme activity, development, oviposition behavior, offspring sex ratios, mobility, navigation and orientation, feeding behavior, learning and immune function , , , –. Reduced immune functioning is of particular interest because of recent disease-related declines of bees including honey bees , . Pesticide and toxin exposure increases susceptibility to and mortality from diseases including the gut parasite Nosema spp. , . These increases may be linked to insecticide-induced alterations to immune system pathways, which have been found for several insects, including honey bees , –.
Surveys of colony food reserves and building materials (i.e. wax) have found high levels and diversity of chemicals in managed colonies , , . These mixtures have strong potential to affect individual and colony immune functioning. However, almost all research to-date on pesticides’ effects on pathogen susceptibility fed a single chemical to test bees . Because pesticides may have interactive effects on non-target organisms (e.g. ), it is crucial to determine how real world combinations and loads of pesticides affect bee health.
One pathogen of major concern to beekeepers is Nosema spp. The endoparasitic fungal infections of N. apis and N. ceranae adversely affect honey bee colony health, and can result in complete colony collapse . Infection with Nosema in the autumn leads to poor overwintering and performance the following spring , and queens can be superseded soon after becoming infected with Nosema . We chose Nosema as a model pathogen because earlier work ,  had demonstrated an interaction with pesticide exposure.
This study addresses two important questions. 1) What types of pesticides might bees be exposed to in major crops? While multiple studies have characterized the pesticide profile of various materials inside a honey bee nest , , few have looked at the pollen being brought back to the nest. 2) How do field-relevant pesticides blends affect bees’ susceptibility to infection by the Nosema parasite?
Pollen was collected from honey bees with permission of the beekeepers and the land owners.
Hive Selection and Pollen Collection
We collected pollen carried by foraging honey bees returning to the hive for nine hives in seven crops: almond, apple, blueberry, cranberry, cucumber, pumpkin, and watermelon (Table 1). For each crop, we selected three fields that were separated by at least 3.2 km. Hives were deployed in these fields for pollination services based on growers’ needs. Within each selected field, we chose the three honey bee hives with the strongest foraging forces by observing flight in the bee yard for 5–10 min, and attached plastic pollen traps (Brushy Mountain Bee Farm, Moravian Falls, NC) to these hives. Pollen traps collect the pollen pellets bees carry on their hind tibiae in flattened regions called corbiculae. Bees use this pollen to make food for larvae inside the nest. We checked traps after three days, and removed them if they contained at least 5 g of pollen. Traps with less than 5 g remained on hives until they contained 5 g of pollen or for 10 days. We placed pollen removed from traps in 50 mL centrifuge tubes and stored the samples on ice until they could be transferred to a −29°C freezer in the lab.
We determined the identity and load of pesticide residues present in pollen samples collected from all crops (except almond). For each field sampled (n = 19), we pooled pollen from the three hives for analysis. One early-season cranberry field and one cucumber field did not yield sufficient pollen in traps for pesticide analysis. Methods follow the LC/MS-MS and GC/MS methods for pollen analysis described in Mullin et al. . We used these data to determine the total number of pesticides detected in each sample, each sample’s total pesticide load, and the diversity and load of pesticides in each of 10 categories: insecticides, fungicides, herbicides, and several insecticide types (carbamates, cyclodienes, formamidines, neonicotinoids, organophosphates, oxadiazines and pyrethroids). To permit comparison between categories with different numbers of elements, we calculated diversity as the proportion of pesticides from a category found in a given sample, and load as the total load divided by the number of chemicals in that category. We only calculated diversity for categories with at least three chemicals.
The total number of pesticides present and total load did not meet parametric assumptions. We thus analyzed how these variables differ between crops using non-parametric Kruskal-Wallis tests. When separated by category and log-transformed, pesticide loads did meet parametric assumptions. We thus determined whether load varied by pesticide category using a general linear mixed model with sample as a random effect, to control for the fact that our regression included one data point per category from each sample. Insufficient degrees of freedom prevented us from expanding this model to include crop. We thus asked whether the pesticide load and diversity varied with crop for each category using one Kruskal-Wallis test per category and applying a sequential Bonferroni correction  across pesticide categories to control for multiple comparisons.
The Nosema infection experiment is similar to published methods . We obtained 210 disease-free honey bees from each of three healthy colonies at the Bee Research Laboratory. Each bee was placed into one of 21 groups upon emergence, with the ten bees in the same group and from the same colony housed together in a wooden hoarding cage (12×12×12 cm). Each group of bees was fed 1 g of pollen mixed with 0.5 mL of syrup (1:1 sucrose to water by weight), which they fully consumed in 2–4 days. These pollen cakes were placed in small petri dishes with the laboratory cages. Pollen from either one of the crop fields or one of two control diets were used. The pollen control group (“BRL”) was fed a mixed pollen diet prepared by the USDA-ARS Bee Research Laboratory. This pollen was collected in the desert Southwest (Arizona Bee Products, Tucson, AZ) and tested as pesticide-free by the USDA Agricultural Marketing Service prior to use. A protein control group was fed an artificial honey bee pollen substitute, MegaBee®. The Nosema inoculum was freshly prepared by mixing Nosema spores isolated from an infected colony (details provided in ) with 50% sucrose solution to obtain a concentration of ca. 2 million spores per 5 mL. We fed 5 mL of the Nosema inoculum to each cage during the first two days of adult life, then provided bees with ad libitum access to clean 50% (w/v) sucrose solution. We collected bees 12 days after infection and examined them for the presence or absence of N. ceranae spores by homogenizing individual abdomens in 1 mL distilled water. Here we focus only on infection prevalence, the number of individuals with Nosema spores.
To look for potential effects of individual pesticides on susceptibility to Nosema infection, we calculated the relative risk and its 95% confidence interval for bees becoming infected after consuming pollen with a specific pesticide. Relative risk measures the chance of developing a disease after a particular exposure , here each pesticide. A relative risk value of one indicates that the probability of infection is equal between exposed and non-exposed groups.
We further tested effects of pesticides in pollen on measured Nosema prevalence using a generalized linear mixed model with a bee’s Nosema status as the response variable, the source hive and pesticide variables as fixed effects, and the pollen sample fed to the bee as a random effect. Collinearity prevented developing a full model to investigate in detail how pesticides and pollen source affect bees’ susceptibility to Nosema infection. We thus selected for analysis two measures that vary with crop and are not nested: total pesticide diversity and fungal load. To graph logistic regression results in a meaningful manner, we followed recent recommendations ,  and a modification of the logi.hist.plot function in the R popbio package  that shows our mixed model output.
All pollen collected in this study contained pesticides (Table 2; mean ± se = 9.1±1.2 different chemicals, range 3–21). Pesticide loads ranged from 23.6 to 51,310.0 ppb (11,760.0±3,734.2 ppb). The maximum pesticide concentration in any single pollen sample exceeded the median lethal dose (LD50, the dose required to kill half a population within 24 or 48 h) for esfenvalerate and phosmet (Table 2). The number of pesticides detected in trapped pollen varied by the crop in which the bee hives were located (Kruskal-Wallis test: H6 = 12.96, p = 0.04), but the total pesticide load did not (H6 = 11.21, p = 0.08)(Fig. 2).
We found insecticides and fungicides in all 19, and herbicides in 23.6% of, pollen samples. Insecticides present in pollen collected by the bees came from seven categories. We found oxadiazines in 10.5%, neonicotinoids in 15.8%, carbamates in 31.6%, cyclodienes in 52.6%, formamidines in 52.6%, organophosphates in 63.2%, and pyrethroids in 100% of pollen samples. Both neonicotinoids and oxadiazines were present only in pollen collected by bees in apple orchards (Figs. 3, S1). Within a sample, pollen fungicide loads were significantly higher than loads of herbicides or any of the insecticide categories (Fig. 4; GLMM, likelihood ratio test: χ2 = 121.9, df = 8, p<0.0001).
147 of the 630 bees (23.3%) fed Nosema spores became infected. 22 of the 35 pesticides (62.9%) found in our pollen samples had relative risk values significantly different from 1 (Table 2). 8 pesticides (22.9%) were associated with increased Nosema prevalence, while the remaining 14 were associated with decreased Nosema prevalence. Two of the three detected pesticides applied by beekeepers to control hive mites (marked with a * in Table 2) had a relative risk larger than two, indicating Nosema prevalence in bees fed pollen containing those chemicals (DMPF and fluvalinate) was more than double the Nosema prevalence in bees that did not consume these chemicals. Of the seven pesticides found in pollen from over half, or at least four, of the crops, the majority were associated with higher Nosema prevalence in bees that consumed them. Both control diets had relative risk values not significantly different from one.
A pollen sample’s fungicide load significantly affected Nosema prevalence among bees fed that pollen (Fig. 5; GLMM, likelihood ratio test: χ2 = 5.8, df = 1, p = 0.02), but pesticide diversity did not (χ2 = 1.7, df = 1, p = 0.19). A bee’s source colony, included as a blocking variable, also did not affect Nosema prevalence (χ2 = 2.0, df = 2, p = 0.36). Replacing fungicide load with chlorothalonil load obtained the same result (chlorothalonil load: χ2 = 5.3, df = 1, p = 0.02; pesticide diversity: χ2 = 1.5, df = 1, p = 0.23; source colony: χ2 = 2.0, df = 2, p = 0.36; fungicide load model AIC = 612.71, chlorothalonil load model AIC = 613.15). Chlorothalonil was also the most abundant fungicide in our samples, and comprised 50.0±10.2% (mean ± se) of the per sample total fungicide load.
The results from this study highlight several patterns that merit further attention. First, despite being rented to pollinate specific crops, honey bees did not always return to the nest with corbicular pollen from those crops. These findings support other research with honey bees and native bees indicating that in some crops native bees may be more efficient pollinators . Second, fungicides were present at high levels in both crop and non-crop pollen collected by bees. Third, two fungicides (chlorothalonil and pyraclostrobin), and two miticides used by beekeepers to control varroa infestation (amitraz and fluvalinate) had a pronounced effect on bees’ ability to withstand parasite infection. Research on pesticides’ effects on bee health has focused almost exclusively on insecticides (e.g. fipronil  and the neonicotinoids imidacloprid ,  and thiacloprid ). Finally, several individual pollen samples contained loads higher than the median lethal dose for a specific pesticide. While multiple studies have shown negative effects of specific pesticides on honey bee individual and colony health , , ,  and high pesticide exposure , , ours is the first to demonstrate how real world pollen-pesticide blends affect honey bee health.
Our results show that beekeepers need to consider not only pesticide regimens of the fields in which they are placing their bees, but also spray programs near those fields that may contribute to pesticide drift onto weeds. The bees in our study collected pollen from diverse sources, often failing to collect any pollen from the target crop (Fig. 1). All of the non-target pollen that we were able to identify to genus or species was from wildflowers (Table S1), suggesting the honey bees were collecting significant amounts of pollen from weeds surrounding our focal fields. The two exceptions to this were hives placed in almond and apple orchards. Almond flowers early in the year, and almond orchards are large, thus providing honey bees with little access to other flowers. Honey bees rarely collect pollen from blueberry or cranberry flowers, which only release large quantities of pollen after being vibrated by visiting bees (buzz pollination) , . Honey bees are not capable of buzz pollination and thus are unlikely to collect large amounts of pollen from these plants to bring back to the colony. Bumble bees, which can buzz pollinate, collect mainly blueberry pollen when placed in blueberry fields . Interestingly, the two crops that saw high levels of pollen collection by honey bees are Old World crops that evolved with honey bees as natural pollinators. Crops native to the New World, where honey bees have been introduced, yielded little or no pollen in our samples.
It is possible that bees were exposed to pesticides while collecting nectar from our focal crops, even when we detected no pollen from those crops. Because pollen traps collect only corbicular pollen intended for consumption by the colony, our data indicate only flowers from which bees are actively collecting pollen and not all flowers they visited. Several studies have detected pesticides in floral nectar and pollen , , sometimes in concentrations with sublethal effects on honey and bumble bees , . Honey bees may collect nectar from blueberry and cranberry flowers via legitimate visits or “robbing” through slits cut at the base of flower corollas . However, exposure to pesticides via nectar may be unlikely in cucumber, pumpkin and watermelon. Beekeepers often report poor honey production when their hives are placed in these crops (pers. obs.).
The combination of high pesticide loads and increased Nosema infection rates in bees that consumed greater quantities of the fungicides chlorothalonil and pyraclostrobin suggest that some fungicides have stronger impacts on bee health than previously thought. Nosema infection was more than twice as likely (relative risk >2) in bees that consumed these fungicides than in bees that did not. Research on the sub-lethal effects of pesticides on honey bees has focused almost entirely on insecticides, especially neonicotinoids . In our study, neonicotinoids entered the nest only via apple pollen. However, we found fungicides at high loads in our sampled crops. While fungicides are typically less lethal to bees than insecticides (see LD50 values in Table 2), these chemicals still have potential for lethal  and sub-lethal effects. Indeed, the fungicides chlorothalonil (found at high concentrations in our pollen samples) and myclobutanil increases gut cell mortality to the same degree as imidacloprid , an insecticide with numerous sub-lethal effects (e.g. , ). Exposure to fungicides can also make bees more sensitive to acaricides, reducing medial lethal doses . In our study, consuming pollen with higher fungicide loads increased bees’ susceptibility to Nosema infection. This result is likely driven by chlorothalonil loads. The pesticide with the highest relative risk was the fungicide pyraclostrobin. Bees that consumed pollen containing pyraclostrobin were almost three times as likely (relative risk = 2.85, 95% CI 2.16–3.75; Table 2) than bees consuming pollen without this chemical to become infected after Nosema exposure. Our results show the necessity of testing for sub-lethal effects of pesticides on bees, and advocate for testing more broadly than the insecticides that are the targets of most current research.
A similarly large increased risk of Nosema infection was associated with consumption of DMPF and fluvalinate, miticides applied by beekeepers to help control the highly-destructive Varroa mite . The path from in-hive application of these miticides to pollen on foragers returning to the hive is unclear. An increasingly popular practice, rotating combs out of hives to remove accumulated pesticides, is expected to reduce miticide levels in hives, and will hopefully decrease spread of these chemicals to the environment. Potential extra-nest sources, however, would slow efforts to reduce miticide accumulation and slow the development of resistance to these chemicals.
Insecticide relative risk values showed an interesting pattern: directional separation by insecticide family. Within a family, relative risk values significantly different than one were almost all in the same direction. The formamidine (DMPF) and two of the three the pyrethroids (bifenthrin and fluvalinate, but not esfenvalerate) were associated with an increased risk of Nosema infection. The carbamate (carbaryl), all neonicotinoids (acetamiprid, imidacloprid and thiacloprid), organophosphates (coumaphos, diazinon and phosmet) and the oxadiazine (indoxacarb) were associated with reduced risk of Nosema infection. Esfenvalerate and coumaphos have previously been found to be associated with apiaries without Colony Collapse Disorder . These patterns suggest that insecticides’ modes of action have differential effects on honey bee immune functioning. Because of the relatively small number of pesticides we found in each insecticide family, however, additional sampling is necessary to determine how robust this pattern is.
The large numbers of pesticides found per sample and the high concentrations of some pesticides are concerning. First, two pollen samples contained one pesticide each at a concentration higher than the median lethal dose. Esfenvalerate (LD50 = 0.13 ppm) was measured at 0.216 ppm in pollen collected by bees in a cucumber field, and phosmet (LD50 = 8.83 ppm) at 14.7 ppm in one apple orchard. While the mean loads for these pesticides are well below their respective median lethal doses (0.0169 ppm for esfenvalerate, 0.7987 ppm for phosmet), our data indicate some bee colonies are being exposed to incredibly high levels of these chemicals. Second, research suggests that simultaneous exposure to multiple pesticides decreases lethal doses ,  or increases supersedure (queen replacement) rate . Our pollen samples contained an average of nine different pesticides, ranging as high as 21 pesticides in one cranberry field. Thus published LD50 values may not accurately indicate pesticide toxicity inside a hive containing large numbers of pesticides. Research looking at additive and synergistic effects between multiple pesticides is clearly needed. Third, pesticides can have sub-lethal effects on development, reproduction, learning and memory, and foraging behavior. The mean and maximum imidacloprid loads in our samples (0.0028 and 0.0365 ppm, respectively) are higher than some published imidacloprid concentrations with sub-lethal effects on honey and bumble bees (0.001–0.0098 ppm , , ).
It is not surprising that total pollen collection varied by crop. Bee foraging activity levels vary with weather , thus outcomes of short-term measurements may be sensitive to temperature, cloud cover or humidity during data collection. Because we collected pollen samples from different parts of the country and on different days, weather conditions undoubtedly differed between crops. Crop flowering timing and landscape-level floral availability can also affect bee activity levels. We focused our analyses on variables less affected by these factors, such as the diversity of pollen types found in samples and the proportion of a sample that was from the target crop.
Our results are consistent with previously published pesticide analyses of pollen collected by honey bees or honey bee nest material , , . The more intensive and geographically more diverse sampling of Mullin et al.  resulted in almost triple the number of pesticides we found, but the average number of pesticides per sample (7.1) is slightly lower than our 9.1. In our study and those listed above, pesticides applied by beekeepers to control hive pests were present in a large proportion of the samples, often in quantities higher than most of the pesticides that are applied to crops.
Our results combined with several recent studies of specific pesticides’ effects on Nosema infection dynamics – indicate that a detrimental interaction occurs when honey bees are exposed to both pesticides and Nosema. Specific results vary, and may depend on the pesticide or dose used. For example, bees exposed to imidacloprid and Nosema can have lower spore counts than bees only infected with the pathogen but also exhibit hindered immune functioning . Our study improves on previous methodologies by feeding pollen with real-world pesticide blends and levels that truly represents the types of exposure expected with pollination of agricultural crops. The significant increase in Nosema infection following exposure to the fungicides in pollen we found therefore indicates a pressing need for further research on lethal and sub-lethal effects of fungicides on bees. Given the diverse routes of exposure to pesticides we show, and increasing evidence that pesticide blends harm bees , , there is a pressing need for further research on the mechanisms underlying pesticide-pesticide and pesticide-disease synergistic effects on honey bee health.
Conceived and designed the experiments: JSP RR DV. Performed the experiments: JSP MA JS DV. Analyzed the data: EML DV. Wrote the paper: JP EML DV.