Imidacloprid is a neonicotinoid pesticide heavily used by the agricultural industry and shown to have negative impacts on honey bees above certain concentrations. We evaluated the effects of different imidacloprid concentrations in sugar syrup using cage and field studies, and across different environments. Honey bee colonies fed sublethal concentrations of imidicloprid (0, 5, 20 and 100 ppb) over 6 weeks in field trials at a desert site (Arizona), a site near intensive agriculture (Arkansas) and a site with little nearby agriculture but abundant natural forage (Mississippi) were monitored with respect to colony metrics, such as adult bee and brood population sizes, as well as pesticide residues. Hive weight and internal hive temperature were monitored continuously over two trials in Arizona. Colonies fed 100 ppb imidacloprid in Arizona had significantly lower adult bee populations, brood surface areas and average frame weights, and reduced temperature control, compared to colonies in one or more of the other treatment groups, and consumption rates of those colonies were lower compared to other colonies in Arizona and Arkansas, although no differences in capped brood or average frame weight were observed among treatments in Arkansas. At the Mississippi site, also rich in alternative forage, colonies fed 5 ppb imidacloprid had less capped brood than control colonies, but contamination of control colonies was detected. In contrast, significantly higher daily hive weight variability among colonies fed 5 ppb imidacloprid in Arizona suggested greater foraging activity during a nectar flow post treatment, than any other treatment group. Imidacloprid concentrations in stored honey corresponded well with the respective syrup concentrations fed to the colonies and remained stable within the hive for at least 7 months after the end of treatment.
Citation: Meikle WG, Adamczyk JJ, Weiss M, Gregorc A, Johnson DR, Stewart SD, et al. (2016) Sublethal Effects of Imidacloprid on Honey Bee Colony Growth and Activity at Three Sites in the U.S. PLoS ONE 11(12): e0168603. https://doi.org/10.1371/journal.pone.0168603
Editor: Gadi V.P. Reddy, Montana State University Bozeman, UNITED STATES
Received: July 21, 2016; Accepted: December 2, 2016; Published: December 28, 2016
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.
Data Availability: All relevant data are within the paper and its Supporting Information files. All data used in inferential statistics are included in Supporting Information.
Funding: The authors received no specific funding for this work.
Competing interests: The authors have declared that no competing interests exist.
Beekeepers and the public have raised concerns about the effects of neonicotinoids, a class of neurotoxic insecticides used as systemic pesticides, on honey bee health. The use of neonicotinoids in agricultural and urban environments has grown dramatically since the introduction of the first, imidacloprid, in 1994; at present, neonicotinoids account for more than 15% of global pesticide sales [1, 2] and imidacloprid is now a top-selling pesticide globally . Neonicotinoids are applied as field sprays and seed treatments , so honey bees may be exposed via direct application, contact with treated surfaces, seed dusts and plant products including pollen, nectar and exudate from extra-floral nectaries [3–8]. Acute toxicity is clearly detrimental to honey bees, but the LC50 of imidacloprid, 1760 ng/L, is above what is considered a “field realistic” range . However, sublethal effects of neonicotinoids have been observed in laboratory experiments  and field studies [9, 10], although few studies have examined colony-level behaviors such as foraging activity or temperature control. Sublethal neonicotinoid pesticides can directly impair honey bee learning and sensory capabilities [11, 12], decrease foraging success and survivorship [8, 10], and can have indirect effects on bees by acting as repellents . Dively et al.  reported effects of sublethal imidacloprid concentrations in pollen diet, which is mainly consumed by newly-emerged bees and nurse bees on brood production, queen replacement, foraging activity and winter survivorship. Pesticides can interact with pests and pathogens; bees exposed to neonicotinoids have been found to have higher Varroa and Nosema densities [2, 14–16]. Synergistic effects of imidacloprid with a pyrethroid on pollinators has also been observed . While size of the worker bee population is often the main criterion used by beekeepers for judging colony health, sublethal pesticide exposure may affect many aspects of honey bee ecology and social organization such as temperature regulation . Sublethal pesticide exposure may have delayed effects as the colony exploits nutritional reserves in honey, pollen, and the bodies of the workers themselves. Few colony-level effects of field-relevant concentrations of imidacloprid exposure via sugar syrup, which is consumed throughout the life of the adult bee , have been reported. Identifying symptoms of sublethal exposure on colony-level behavior would help beekeepers detect such exposure.
Two bee cage studies were conducted to measure survivorship and food consumption by newly-emerged adult bees exposed to various concentrations of imidacloprid in sugar syrup. Four field studies were conducted in which honey bee colonies were exposed to imidacloprid concentrations in three different landscapes in the southern half of the U.S. for a qualitative examination of the role of environment. Two studies were conducted at sites in Arkansas and Mississippi with abundant forage, both agricultural and natural, and two studies conducted at a relatively isolated site in southern Arizona with no commercial agriculture and little natural forage during the treatment period. Capped brood levels and other metrics of colony phenology and size were measured at all sites, and continuous weight and temperature data, which have been shown to reflect honey bee colony growth, adult bee population size, foraging activity and brood production [20, 21] were monitored at the Arizona sites. Coumaphos, often used against bee pests and a common contaminant of hive products [22, 23], was applied to colonies in Mississippi before the imidacloprid application and those data were evaluated for interactions with imidacloprid. Proper apiary management at each site was given primacy over consistency across sites. Imidacloprid concentrations were chosen to include both a low (5 ppb) concentration, observed in field samples  and a concentration high enough (100 ppb) was used here to increase the chances of provoking a response in terms of hive growth and activity, and to harmonize to some extent this study with the experimental designs of other workers (e.g., ). A summary of the experiments is provided (S1 Table).
Materials and Methods
Two cage studies were conducted: the first from 11 Sept. to 31 Oct. 2014 in which groups of 100 bees were fed sugar syrup (see S2 Table) containing either 100, 5 or 0 ppb imidacloprid (10 replicate cages per treatment) and the second from 27 July to 16 Sept 2015 in which groups were fed syrup containing either 100, 20, 5 or 0 ppb imidacloprid (7–8 replicate cages per treatment). In each experiment, one frame of mature brood was removed from each of four 6-month-old colonies with Cordovan Italian queens (C.F. Koehnen & Sons, Inc) established at the Carl Hayden Bee Research Center, Tucson, AZ and placed in an incubator (Percival model I36VL) at 30°C and 50% r.h. Adult bees emerging over the following 48 h were distributed among 30 Plexiglas® cages (internal volume 785 cm3) until each cage had 100 bees. Each cage had plastic feeding bottles on top containing 15 mL sugar solution and 50 mL water, and a 4x4 cm square of wax foundation hung vertically in the center. A mixture of 10g each of pollen (Natural Foods Inc., Toledo, OH), sucrose and inverted sucrose was placed inside a rubber gasket accessed via a hole in the side of the cage. Dead bees were removed and counted daily. Syrup consumption was measured weekly by weighing bottles of syrup before and after use; vials were emptied and refilled with fresh syrup. Consumption per bee was calculated as the observed consumption for a given cage divided by the number of “bee-days” for that time period, in which a bee-day represents one bee alive for one day in that cage. Water consumption was measured the same way but less often and only total water consumption was analyzed. Five to eight live bees were removed from each cage after 4 and 6 weeks; those data were censored. Bee samples were submitted to the Laboratory Approval and Testing Division, Agricultural Marketing Service, USDA (LATD) for residue analysis. Adult bee survivorship in cage studies was analyzed using Proc LifeReg (SAS Inc. 2002). An appropriate distribution was first chosen to model survivorship, survivorship curves were generated for each replicate cage based on that model, and treatments compared using ANOVA (Proc Glimmix, SAS Inc. 2002) with respect to the 30th and 50th percentiles, and shape was estimated by subtracting the 40th percentile from the 30th percentile. Weekly syrup consumption per bee was analyzed using repeated measures MANOVA.
Two studies were conducted. In the first study, from May 2014 to March 2015 (hereafter 2014 experiment), bee colonies were fed sugar syrup containing either 100, 5 or 0 ppb imidacloprid. In the second study, from June to September 2015 (hereafter the 2015 experiment), colonies were fed syrup containing either 100, 20, 5 or 0 ppb imidacloprid. In each study, N = 4 bee colonies per treatment group. All colonies were established from packages (C.F. Koehnen & Sons, Inc., Glenn, CA 95943) three months prior to treatment in painted, 10-frame, wooden Langstroth boxes (43.7 l capacity) (Mann Lake Ltd,) with migratory wooden lids at the Santa Rita Experimental Range south of Tucson, AZ (31°46'39"N, 110°51'46"W). At establishment each colony was given three frames of drawn comb with six frames of plastic comb foundation and fed 2 kg sugar syrup (1:1 w:w) and 250 g pollen patty (see Supplementary material 1). The apiary was surrounded by native, unmanaged plants, particularly mesquite (Prosopis spp.). After 4–6 weeks, hives were placed on stainless steel electronic scales (TEKFA® model B-2418 and Avery Weigh-Tronix model BSAO1824-200) (max. capacity 100 kg) connected to 12-bit dataloggers (Hobo® U-12, Onset Computer Corporation), set to record weight every 15 minutes. The system had an overall precision of approximately ±20 g. On the same day, a temperature sensor (iButton Thermochron, precision ±0.06°C) enclosed in brass mesh was stapled to the center of the top bar on the 5th frame in each hive and set to record every 30 min. Hives were inspected at 4-week intervals starting 5 weeks before treatment, using a published protocol . For each inspection, frames were gently shaken to dislodge adult bees, then weighed, photographed using a 16.3 megapixel digital camera (Pentax K-01, Ricoh Imaging Co., Ltd.) and replaced in the hive. The area of sealed brood per frame was estimated from the photographs using ImageJ version 1.47 software (W. Rasband, National Institutes of Health, USA). The total weight of the adult bee population was calculated by subtracting the combined weights of hive components (i.e. lid, inner cover, box, bottom board, frames, entrance reducer, internal feeder) obtained at the start of the experiment (model EC15, OHaus) from the total hive weight recorded the midnight prior to the inspection. At each inspection, 3–5 g of adult bees, wax and honey were each collected from each hive into 50 ml centrifuge tubes and stored at -20°C. Wax was collected from all frames with comb pre-treatment, and from recently built comb during and post treatment. Adult bees were collected from a frame next to the brood cluster, and nectar was collected from capped cells pre-treatment and from recently filled cells during and post-treatment. Pooled samples collected prior to treatment were analyzed for residues of 174 compounds by LATD; later samples were analyzed only for neonicotinoid pesticides and breakdown products.
Hives were separated into four groups and randomly assigned to treatment group. Hives within a group were 0.5–1 m apart and groups were >3 m apart. Just prior to treatment all broodless frames containing honey and/or pollen were replaced with frames of empty drawn comb collected earlier from the same apiary, leaving colonies with an average (±s.e.) of 1620 (±246) g food in 2014 and 1561 (±145) g food in 2015. Disposable latex gloves, hive tools and sampling equipment used in hive inspections were changed for each treatment group; hives were inspected from lowest to highest concentration treatments. Starting 17 July 2014 for the first experiment and 9 July 2015 for the second, colonies were given 2 kg of treated syrup twice per week for the first 2 weeks and 3 kg twice per week for the last 4 weeks, and not fed thereafter. Infestation levels of Varroa destructor were monitored on four occasions during and post treatment in 2014 and once during treatment in 2015 by placing adhesive cardboard strips (Mann Lake Ltd) under the brood box for 3 d and counting the mites. Colonies in the first experiment were treated with amitraz (Apivar®) on 4 Nov. The 2015 study was terminated shortly after the end of treatment owing to visits by a large animal, probably a bear.
Continuous data were divided into daily average data and within-day detrended data. Detrended data were calculated as the difference between the 25 hour running average and the raw data . Sine curves were fit 3-day subsamples of detrended data taken sequentially by day (see ) and curve amplitudes, representing estimates of daily variability, were used as a response variables. Weight and temperature amplitude datasets were reduced to a data point every 5 days for repeated measures analysis to ensure no overlap between 3 d samples. For consistency, running average data were treated in the same fashion. Ambient weather data for 2014 was obtained for comparison (AmeriFlux US-SRM Santa Rita Mesquite, doi:10.17190/AMF/1246104).
In May 2015, 16 honey bee colonies were established from packages with Italian queens from the same breeding line obtained from a commercial breeder (Diamond Lakes Apiaries, Murfreesboro, AR) in new 10-frame Langstroth deep boxes (Dadant & Sons, Inc., IL) with plastic comb foundation near Lonoke, AR (34°40'38"N, 91°55'21"W) and fed sugar syrup (1:2 w:w). Forage consisted of cotton, soybeans and rice in addition to trees and wildflowers. Colonies were inspected prior to treatment on 7 July. At each inspection, adult bee, honey and wax were sample, and hive frames weighed and photographed, as described above for the Arizona experiments; samples were analyzed for all pesticide residues at LATD. Honey frames were not removed prior to treatment. Amitraz (Apivar) strips were installed in colonies prior to treatment, and the strips left in place during syrup feeding. Hives were weighed 14 July and every 14–18 days thereafter until 29 Sept. Colonies were randomly assigned to treatment groups and starting 15 July were fed 2 kg sugar syrup with imidacloprid concentrations of 0, 5, 20 or 100 ppb twice per week, until 3 Aug., after which the quantity of syrup was increased to 3 kg twice per week until 24 Aug. Hives were inspected on 6 Aug., 28–29 Aug. and 29–30 Sept, and samples taken on 31 Aug. and 23 Sept.
In April 2015, 15 honey bee colonies, (Group A) were established from packages with Italian queens from the same breeding line (Gunter Honey, Inc., Lumberton, MS) in new 10-frame Langstroth deep boxes (Dadant & Sons, Inc.) with plastic comb foundations at the Mississippi State University Coastal Experiment Station in McNeill, MS (30°39'46"N, 89°38'01"W). Forage consisted of Chinese tallow trees (Triadica sebifera) in June and goldenrod (Solidago sp.) in Oct. In June a second group of 13 bee colonies (Group B) was established at the same site. Group A colonies were inspected 21 May and 23 June. At each inspection, hive strength was estimated by counting “frame spaces” (spaces between frames occupied by adult bees, with the spaces between the frame and box on each side being 0.5 spaces each; the values are whole numbers between 0 and 10 inclusive). Capped brood surface area was measured by placing a transparent plastic sheet marked with regular hexagons (23.4 cm2 each) on the comb and counting occupied hexagons, rounding up if more than half filled (see ). Adult bees were sampled near brood frames and honey was sampled from brood frames; samples were stored at -20°C. Varroa densities were assessed in Aug. and Oct. by collecting ≥230 adult bees into 1 L Mason jars with 70% ethanol, agitating the jar, and counting mites and bees.
Eight randomly-selected colonies from Group A received 1 kg carbohydrate patty (Pro Winter, Mann Lake Ltd) mixed with coumaphos at 5.8 ppm (see S2 Table), a concentration reported from pollen , on 21 May, and again four weeks later. The remaining seven colonies received coumaphos-free carbohydrate patties. On 16 July, all 28 Group A and Group B colonies were inspected and divided into seven treatment groups of four colonies each. Frames containing only stored honey were replaced with empty drawn comb. On 17 July, the eight Group A colonies exposed to coumaphos were split into two groups, with four colonies getting coumaphos at 5.8 ppm in addition to imidacloprid at 5 ppb, and the other four getting 5.8 ppm coumaphos plus 20 ppb imidacloprid. Of the seven colonies in Group A that did not get coumaphos, four were given 5.8 ppm coumaphos syrup and the remaining three colonies given syrup with no additives. The same day, the 13 colonies in Group B were divided into three groups of four colonies, and the colonies received syrup with imidacloprid at 5, 20 or 100 ppb, respectively, with the remaining colony included as an untreated control (see S1 Fig). During treatment, each colony received 1 kg syrup three times per week until the end of treatment period on 19 Aug. for a cumulative total of 12 kg. All 28 colonies were inspected on 17 Aug.,17 Sept. and 19 Oct. Colonies received a further five 1 kg feedings of 2:1 sugar syrup from 20 Aug. to 14 Oct. Adult bee and honey samples were pooled within treatment group for each sampling date and analyzed at LATD for residues of imidacloprid, coumaphos and breakdown products. The effects of treatment on brood surface area were evaluated using mixed-model ANOVA, with frame space estimates prior to treatment used as covariates. Frame spaces were analyzed using non-parametric tests (Proc Npar1way, SAS Inc. 2002).
Use of the apiary site in Arizona was granted through Range Use Agreement with the University of Arizona Agricultural Experiment Station, Santa Rita Experimental Range. Use of the apiary site in Arkansas was granted through an agreement with a private landowner and the site is registered with the Arkansas State Plant Board, in accordance with Arkansas apiary law. The use of the apiary site in Mississippi was granted through the USDA Agricultural Research Service. The land use did not involve endangered or protected species.
Experiments were terminated after 50 d with 25–70 remaining alive per cage in 2014 and 23–64 per cage in 2015. Bees collected before treatment in both trials tested negative for neonicotinoid pesticides or breakdown products, and no imidacloprid was detected in any samples from 2015. Only bees collected from the 100 ppb group in 2014 during treatment tested positive (< 2 ppb), possibly due to contaminated syrup in their digestive tracts. A Weibull distribution with censoring was selected to model survivorship curves for each replicate. The intercept and scale parameters from the fit for each replicate were significant (P<0.001). The standard error of the intercept as a percentage of the intercept ranged from 1.1 to 10.9% for the 2014 data and from 1.1 to 6.8% for the 2015 data. For the scale parameter those values ranged from 5.2 to 24.9% for the 2014 data and from 10.7 to 20.7% for the 2015 data (see Figs A-D in S1 File). The 0, 5 and 100 ppb treatments were analyzed across both experiments, and the analysis was also conducted using the 0, 5, 20 and 100 ppb data from the second experiment. Treatment did not have a significant effect on the 30th percentile, 50th percentile, or shape response variables for either analysis (S3 Table).
Syrup and water consumption were likewise compared 1) across experiments for the 0, 5 and 100 ppb groups; and 2) within the second experiment for the 0, 5, 20 and 100 ppb groups. In the first analysis, treatment (F2,347 = 19.68, P<0.0001) and experiment (F1,349 = 182.40, P<0.0001) were significant (Fig 1). Bees in the 100 ppb group consumed less syrup on average than those in either the 5 ppb or control groups (P<0.0001 for both). In the second analysis, no significant differences were observed (P = 0.14). Water consumption was significantly different among groups in the first analysis (F2,47 = 3.79, P = 0.0298) and between years (F1,47 = 36.07, P<0.0001) (S2 Fig). Bees in the 100 ppb group drank more water than bees in the control (P = 0.0419) and bees in 2015 drank more than those in 2014. In the second analysis, water consumption was significantly different among groups (F3,26 = 3.20, P = 0.0400) but no pairwise contrasts were significant.
Adult honey bees were kept in cages and fed sugar syrup with imidacloprid concentrations of either 100 ppb, 20 ppb, 5 ppb or 0 ppb. Each cage was stocked with 100 bees (7–8 cages per group). A) Experiment conducted Sept-Oct. 2014 (no 20 ppb treatment); B) Experiment conducted Aug.-Sept. 2015.
No agrochemical residues were detected in adult bees or honey in 2014, but samples of wax pooled among hives contained miticides fluvalinate (32–34 ppb) and coumaphos (9 ppb) and the fungicide carbendazim (16–65 ppb). One sample also had an insect growth regulator (IGR) hydroprene, at 144 ppb. Wax samples from the 2015 experiment contained four varroacides: fluvalinate (up to 781 ppb), thymol (up to 648 ppb), fenpyroximate (up to 33 ppb) and trace amounts of coumaphos, as well as carbendazim (up to 27 ppb) and trace amounts of hexythiazox, an IGR. Honey from 2015 showed trace carbendazim. Thymol was detected in adult bees (194–392 ppb) but not applied. Imidacloprid was detected in adult bees only in the 100 ppb treatment in both field experiments (Table 1). In 2014, imidacloprid concentrations remained largely stable in stored honey through the following March. Imidacloprid was detected in low levels in wax in the100-ppb group post treatment in 2014 but not in 2015.
Adult bee population (Fig 2), brood area (Fig 3) and average frame weights from the 2014 study were divided into 4 time periods: 1) pre-treatment (22 May to 16 July); 2) during treatment (13 Aug.); 3) post treatment (9 Sept. to 4 Nov.); and 4) post winter (26 Feb. to 23 Mar.) (Table A in S2 File). Pre-treatment no differences were found among groups. During treatment, brood area was significantly lower on average in the 100 ppb group than the 5 ppb and control groups (a weight error prevented adult bee population calculation in Aug.) (see Table 2 and Table B in S2 File). Post treatment, adult bee populations and brood areas were significantly lower among 100 ppb colonies, after controlling for pre-treatment values, than control colonies. Brood area and average frame weight, after controlling for pre-treatment frame weight, were lower in the 100 ppb treatment than the 5 ppb treatment. Post winter, capped brood areas were significantly lower in the 100 ppb group than either the 5 ppb or control groups. Average daily Varroa mite fall over four sampling occasions (Aug., Oct., Nov. and Mar.) were not different across treatments (P = 0.08) but were across sampling occasions (F3,10 = 42.82, P<0.0001) with the highest being an average of 19.0 mites per day in Oct. Visual assessment of the colonies revealed no disease symptoms.
The syrup had imidacloprid concentrations of either 100 ppb, 5 ppb or 0 ppb (four colonies per group), in southern Arizona. A) Experiment initiated in 2014; B) Experiment initiated in 2015. Gray zone indicates treatment period. * indicates that masses were estimated using detrended weight amplitudes (see Meikle et al. 2016) rather than inspection data.
The syrup had imidacloprid concentrations of either 100 ppb, 5 ppb or 0 ppb (four colonies per group), in southern Arizona. A) Experiment initiated in 2014; B) Experiment initiated in 2015. Gray zone indicates time period during which the syrup was given.
Hourly hive weight data from May 2014 to March 2015 showed similar trends (Fig 4). Daily weight amplitudes were divided into three periods: 1) pre-treatment (22 May to 16 July); 2) during treatment (17 July to 26Aug.); and 3) a post treatment nectar flow (3 Sept. to 28 Oct.), with a nectar flow being defined as > 7 d of largely monotonic increases in hive weight (see ). No differences were observed among treatments pre- or during treatment (Tables A and B in S3 File), but during a nectar flow starting mid-September, weight amplitudes were significantly higher for colonies in the 5 ppb treatment than for colonies either the 100 ppb or control groups (Fig 5), after controlling for the adult bee population, suggesting greater foraging activity. The data also revealed the robbing of at least one hive in the 100 ppb concentration group. A colony fed syrup is expected to lose weight as moisture in the syrup is driven off by the bees, but a loss rate, even in a desert environment, of > 1.5 kg per day is excessive. A loss rate of 2 kg in < 6 h (Fig 6) indicated robbing.
The syrup had imidacloprid concentrations of either 100 ppb, 5 ppb or 0 ppb (four colonies per group), in southern Arizona. A) Experiment initiated in 2014; B) Experiment initiated in 2015. Gray zone indicates treatment period.
Data were collected during a nectar flow in Sept.-Oct. 2014 from bee colonies given sugar syrup containing imidacloprid concentrations of either 100 ppb, 5 ppb or 0 ppb (four colonies per group), in southern Arizona.
Data collected from honey bee colonies involved in a feeding experiment Aug.- Sept. 2014. “High concentration” colonies were fed sugar syrup containing imidacloprid at 100 ppb and “Control” colonies fed sugar syrup alone
Temperature data for the 2014 experiment were also divided into three periods: 1) 1) pre-treatment (22 May to 16 July); 2) post treatment (27 Aug. to 30 Nov.); and 3) winter (1 Dec. to 26 Feb.); some data were missing due to equipment problems. Running average temperatures did not differ among treatments during any of those periods (S3 Fig) nor did detrended amplitudes pre- or post-treatment. However, amplitudes in colonies in the 100 ppb treatment group were significantly higher in winter, indicating reduced temperature control, than those in the 5 ppb treatment (Fig 7) but not the control. One hive (100 ppb treatment) died during winter and those data were removed from that analysis.
Data were collected during winter (Dec. 2014 to Feb. 2015) from honey bee colonies fed sugar syrup containing either 100 ppb, 5 ppb or 0 ppb imidacloprid during July-Aug. 2014.
Discrete data from the 2015 study in Arizona were divided into two time periods, owing to the early termination: 1) pre-treatment (3 June to 8 July); and 2) during treatment (5 Aug. to 2 Sept., which included the final inspection). As in 2014, no significant pre-treatment differences were observed and adult bee populations in the 100 ppb group during treatment were significantly lower than other groups after controlling for pre-treatment values (Tables A and B in S4 File). Similar results were observed for brood surface area. Average frame weights were not affected, but differences may have been obscured by syrup storage from feeding. Average daily Varroa fall data were not different among treatments (P = 0.81).
All syrup was consumed in the 2014 treatment period but not in 2015. Unconsumed syrup was measured on 7 and 24 Aug.; the syrup was replaced with fresh syrup to maintain similar imidacloprid exposure between the two experiments. Colonies in the 100 ppb treatment left on average 4.21 kg of syrup while colonies in the other groups left 0.50–0.80 kg. These data were analyzed with similar data from Arkansas (see below).
Continuous weight and temperature data in 2015 were grouped into the same pre-treatment and during treatment time periods as the hive inspection data. Weight amplitudes were not different among treatments, but the lack of a nectar flow would have made any differences more difficult to detect. Average temperatures, with adult bee population as a covariate, were not different pre-treatment but were significantly lower for the 100 ppb treatment than for the 20 ppb, 5 ppb and control groups during treatment. Temperature amplitudes in the 100 ppb treatment were again significantly higher than either the 5 ppb treatment or control but not the 20 ppb treatment (Fig 8).
The syrup contained either 100 ppb, 20 ppb, 5 ppb or 0 ppb imidacloprid during July-Aug. 2015. A) Amplitudes of sine curves fit to detrended hourly internal temperature data; B) Average daily hive temperature.
Neonicotinoid residue analysis differed in two important respects from the Arizona studies: 1) imidacloprid levels were lower in Arkansas post treatment, particularly in honey (Table 3) than in the Arizona studies, and were also lower in the syrup itself for the 100 ppb treatment; and 2) measurable amounts of another neonicotinoid pesticide, clothianidin, were detected with imidacloprid in some samples of bees, wax and nectar, always at concentrations of <7 ppb.
No differences among groups with respect to brood surface area, average frame weight or total hive weight were detected, after controlling for pre-existing differences, after the start of treatment. The exposure of colonies to clothianidin may have affected differences among treatment groups. Unconsumed syrup was measured on 24 July, 20 and 28 Aug. Colonies in the 100 ppb treatment left on average over 5.42 kg of syrup while colonies in the other treatment groups left on average 1.05 to 1.74 kg. Treatment group had a significant effect on the log amount of consumed syrup (F3,27 = 7.92, P = 0.0017) but site did not and, as noted above, no pre-existing size differences among groups were detected. Colonies in the 100 ppb group consumed less syrup than those in other groups (P = 0.0042, 0.0013 and 0.0047 for the 20 ppb, 5 ppb and control groups, respectively).
Coumaphos residues were detected in honey and adult bees sampled before treatment (Table 4). Both coumaphos and imidacloprid were detected in honey and adult bees post treatment. Imidacloprid levels were consistent with those observed in Arizona; it was also detected in honey sampled from colonies that had not been fed imidacloprid, indicating robbing, so comparisons between the 0 and 5 ppb imidacloprid treatments should be interpreted cautiously.
Exposure to coumaphos in patties did not have a significant effect for the May, June or July datasets on capped brood area (P = 0.35) or on the number of frame spaces (P = 0.22, 0.66 and 0.10, respectively, using nonparametric ANOVA). Analysis of all groups after syrup feeding showed a significant effect of imidacloprid, but not coumaphos, and a significant interaction of imidacloprid and coumaphos (Tables A and B in S5 File.) (Fig 9). Only one post hoc contrast was significant, showing that colonies not given coumaphos but fed 20 ppb imidacloprid syrup had significantly less brood than those given both coumaphos and 5 ppb imidacloprid. Excluding colonies fed coumaphos, control colonies had more capped brood on average than colonies in other groups and colonies fed 5 ppb had more brood than those fed 100 ppb (S4 Table). Varroa densities, measured as mites per bee, were not different across treatments (P = 0.51) or sampling occasions (P = 0.052); densities reached 1.9 mites per bee in Oct.
Colonies were fed either 100, 20, 5 or 0 ppb imidacloprid in sugar syrup over approximately 6 weeks at two locations. A) Mississippi apiary, which also included some hives given coumaphos at 5.8 ppm; B) Arkansas apiary. Gray area indicates treatment period.
Sublethal concentrations of neonicotinoid pesticides can affect several aspects of honey bee behavior and ecology. The focus of this work was on colony-level behaviors such as brood rearing, foraging activity and temperature management over time, across different environments. Cage studies in controlled conditions were conducted to focus on adult bee survivorship and consumption. No differences were found in survivorship. Bees in cages fed 100 ppb imidacloprid consumed less syrup on average than those in other groups but the differences were not large. Interestingly, caged bees in the 2015 study consumed less syrup on average, 18.5 mg per bee per day, than those in the 2014 study, 22.6 mg; the reason for this difference is not known. Caged bees ate on average 420 mg (100 ppb group) to 480 mg (control group) sugar from emergence to 50 d. Sugar consumption by an adult honey bee has been estimated to be 534 to 1005 mg (pollen and nectar foragers, respectively) from emergence to 31 d . The lower consumption rate of syrup by caged bees would result in a lower imidacloprid exposure compared to actively-foraging field bees. This was born out in residue analyses. Caged bees were restricted to treatment syrup but imidacloprid was detected only once in two studies (at <2 ppb in the 100 ppb treatment) while bees from that treatment group in Arizona field studies had 6.8 to 19.3 ppb imidacloprid. Imidacloprid is metabolized within hours by bees  so those concentrations probably represent mostly undigested syrup in the gut.
Determination of colony-level effects of sublethal exposure is difficult, largely due to natural variation among colonies and to uncontrolled factors such as differential exposure of experimental colonies to exogenous agrochemicals, and up to 24 to 80 colonies per treatment have been suggested . We observed differences with far fewer colonies per treatment although not at all sites. In Arizona, colonies in the 100 ppb imidacloprid treatment group had significantly less brood than other groups during and post treatment. Colonies in the 2014 experiment also had fewer adult bees and lighter frames post treatment, indicating less stored food. Temperature variability, as measured by daily temperature amplitudes, were higher on average in the 100 ppb group than those in the 5 ppb group during winter for the 2014 study and during treatment for the 2015 study. Temperature data can be affected by the distance of the bee cluster to the sensor, which can vary among colonies and over time ; we assumed that distance was random with respect to treatment. Low temperature variability has been associated with brood rearing and high variability with little or no brood .
High coumaphos concentrations, such as 100 ppm or greater, affect queens by reducing oviposition and ovarian weight, causing early supercedure and queen cell rejection , and increasing worker mortality . In the Mississippi study presented here, coumaphos levels were lower (5.8 ppm) and coumaphos alone did not have a significant impact on the amount of brood, although its interaction with imidacloprid was significant. Among colonies not fed coumaphos, imidacloprid exposure had a negative effect on brood production at all concentrations. In interpreting these results, the effects of unknown local conditions cannot be ruled out for two reasons: 1) differences in capped brood area in coumaphos-free groups fed syrup with 5 and 0 ppb imidacloprid were not observed at any other sites; and 2) imidacloprid was detected in honey from both the coumaphos positive and coumaphos negative control groups by the end of treatment, indicating robbing. No significant differences in brood or frame weights were observed in Arkansas. Experimental execution and environment between Arkansas and the other sites differed in several respects: 1) Arkansas colonies were not deprived of stored honey; 2) syrup in the 100 ppb treatment in Arkansas apparently had less imidacloprid than planned; and 3) Arkansas colonies were exposed to other agrochemicals, including another neonicotinoid, clothianidin.
Syrup consumption differed among treatment groups, sites and between years. Colonies in the 100 ppb treatment groups in Arizona in 2015 (but not 2014) and in Arkansas left more syrup than colonies in other groups, on average by a factor of at least three. Colonies in Arizona had few alternatives for nectar; the lack of consumption may have been due to lower activity levels (reduced foraging activity was observed) or to a dislike for the treated syrup . Reduced consumption was observed in the 100 ppb group in the cage studies, and newly-emerged honey bees can detect imidacloprid at < 3 ppb . Colonies in Arkansas may have sought other nectar sources. Some colonies in Mississippi may have avoided high imidacloprid exposure by absconding: hives exposed to 0, 5, and 20 ppb retained their original queens but three of four hives fed 100 ppb syrup absconded by September, one of which was recovered later with the original queen and overwintered successfully in new equipment. Interestingly, consumption rates were also different across years for caged bees, which would not have had the option of alternative forage or absconding.
Colonies in the 5 ppb and control groups had more brood area, heavier frames and greater temperature control than colonies in the 100 ppb groups. In addition, colonies in the 5 ppb group had higher daily weight amplitudes, which is associated with flight activity, during a nectar flow post treatment (in 2014) than colonies in other groups, including the control. Average frame weights, which would reflect nectar storage, were not different between the 5 ppb and control groups, so either the flight activity did not increase foraging success, or consumption by the colonies was higher. Imidacloprid is a nicotine analog and affects insect nicotinic acetylcholine receptors . Honey bees have been found to have a tolerance for the “trace” concentrations of nicotine, 0.1 to 5 ppm, found naturally in some nectars  and the survival of bees from “weak” colonies (colonies exhibiting a lack of vigor for unknown reasons) improved when placed in cages and given a 300 μM solution of nicotine, or about 48.6 ppm . That concentration is several magnitudes higher than a 5 ppb imidacloprid solution.
These experiments involved a particular kind of exposure to imidacloprid (via syrup) with exposure to constant concentrations over 5–6 weeks. Honey bees can be exposed to pesticides via other means, including in pollen, seed treatment dusts and water and means of exposure may play a role in the impact on colony health. Exposing bees via contaminated supplemental pollen feed (e.g., [2, 32]) may impact nurse bees and larvae more than it would foragers because those life stages feed more on pollen than foragers [33, 34]. While an effort was made to harmonize colony management among sites in this study, colonies had site-specific pesticide exposures, due to agriculture, pre-existing contaminants in drawn comb, and bee pest treatments that should be taken into account when interpreting the data across sites. Regarding the treatments selected here, even in agriculturally intensive areas, such as the southeastern U.S., concentrations as high as 100 ppb have seldom been reported in pollen or nectar. Stewart et al.  found that neonicotinoid levels in soils and wild flowers after seed treatment applications in experimental and commercial plots seldom exceeded 50 ppb. Bees exposed to high levels of imidacloprid in nectar also may simply reduce or avoid feeding on contaminated nectar, whether or not alternative forage is available . Detection of effects of low concentrations of a neonicotinoid pesticide on colony level behavior suggests that the interaction between pollinators and commercial agriculture is more complex than previously thought.
- Imidacloprid concentrations were stable in stored nectar or honey in the hive environment for at least 7 months;
- Bee colonies with little stored honey and fed an imidacloprid concentration of 100 ppb in syrup generally had lower adult bee populations, less capped brood and worse winter thermoregulation compared to controls;
- Syrup consumption differed with respect to imidacloprid concentration as well as year in cage studies, and with respect to imidacloprid concentration, year and site in field studies; caged bees or bee colonies fed 100 ppb syrup often consumed less than other groups;
- Within-day hive weight changes were significantly higher, indicating increased foraging activity, in the 5 ppb treatment than either the 100 ppb treatment or the control during a nectar flow.
S1 Fig. Diagram of treatments for Mississippi trial.
S2 Fig. Total water and syrup consumption for bees kept in cages and fed sugar syrup with imidacloprid concentrations of either 100 ppb, 20 ppb, 5 ppb or 0 ppb.
S3 Fig. Average daily internal temperatures for bee colonies subjected to one of three treatments: fed syrup containing 100 ppb imidacloprid, fed syrup containing 5 ppb imidacloprid, and fed syrup with no imidacloprid.
Figs. A-D. Intercept and scale parameters, and associated distribution scale and shape parameters, for a regression-based fit of a Weibull distribution with censoring to bee survivorship in cage studies over two trials.
TablesA and B. Analysis and post hoc contrast results for hive inspection data from a field experiment conducted in Arizona 2014.
Tables A and B. Analysis and post hoc contrast results for continuous weight and temperature data from a field experiment conducted in Arizona 2014.
Tables A and B. Analysis and post hoc contrast results for hive inspection and continuous monitoring data from a field experiment conducted in Arizona 2015.
Tables A and B. Analysis and post hoc contrast results from a field experiment conducted in Mississippi 2015.
S6 File. Experimental data.
S1 Table. Summary of experimental designs.
S2 Table. Preparation of pesticide treatments.
S3 Table. Analysis of select parameters from Weibull distributions to cage survivorship data.
S4 Table. Analysis and post hoc contrast results showing the effect of exposure to imidacloprid on capped brood area for hives not exposed to coumaphos from a field experiment conducted in Mississippi 2015.
The authors would like to warmly thank M. Heitlinger at the Santa Rita Experimental Range and M. Mcclaran at the University of Arizona, and E. Willman, Jr in Arkansas, for permission for land use; D. Boykin for crucial statistical guidance, R. Scott for weather data, and A. Barron, R. Johnson, K. Kaplan and four anonymous reviewers for their tremendous help improving the manuscript.
- Conceptualization: WGM MW JJA DRJ SDS MJC GML.
- Data curation: WGM MW JJA AG DRJ.
- Formal analysis: WGM JJA AG DRJ.
- Funding acquisition: WGM JJA SDS GML.
- Investigation: WGM MW JJA AG DRJ JZ GML.
- Methodology: WGM MW DRJ JZ.
- Project administration: WGM JJA.
- Resources: WGM JJA GML.
- Software: WGM.
- Supervision: WGM JJA SDS GML.
- Validation: JJA AG DRJ.
- Visualization: WGM MW JJA MJC.
- Writing – original draft: WGM JJA.
- Writing – review & editing: WGM MW JJA MJC.
- 1. Blacquière T, Smagghe G, van Gestel CAM, Mommaerts V (2012) Neonicotinoids in bees: a review on concentrations, side-effects and risk assessment. Ecotoxicology 21:973–992. pmid:22350105
- 2. Dively GP, Embrey MS, Kamel A, Hawthorne DJ, Pettis JS (2015) Assessment of chronic sublethal effects of imidacloprid on honey bee colony health. PLoS ONE 10(3): e0118748. pmid:25786127
- 3. Stoner KA, Eitzer BD (2012) Movement of soil-applied imidacloprid and thiamethoxam into nectar and pollen of squash (Cucurbita pepo). PLoS ONE 7(6): e39114. pmid:22761727
- 4. Krupke C.H, Hunt GJ, Eitzer BD, Andino G, Given K (2012) Multiple routes of pesticide exposure for honey bees living near agricultural fields. PLoS ONE 7(1): e29268. pmid:22235278
- 5. Girolami V, Mazzon L, Squartini A, Mori N, Marzaro M, DiBernardo A, et al. (2009) Translocation of neonicotinoid insecticides from coated seeds to seedling guttation drops: A novel way of intoxication for bees. J Econ Entomol 102:1808–1815. pmid:19886445
- 6. Reetz J, Zühlke S, Spiteller M, Wallner K (2011) Neonicotinoid insecticides translocated in guttated droplets of seed-treated maize and wheat: a threat to honeybees? Apidologie 42:596–606.
- 7. Pochi D, Biocca M, Fanigliulo R, Pulcini P, Conte E (2012) Potential exposure of bees, Apis mellifera L., to particulate matter and pesticides derived from seed dressing during maize sowing. Bull Environ Contam Toxicol 89:354–361. pmid:22562269
- 8. Zhu YC, Adamczyk J, Rinderer T, Yao J, Danka R, Luttrell R, Gore J (2015) Spray toxicity and risk potential of 42 commonly used formulations of row crop pesticides to adult honey bees (Hymenoptera: Apidae). J Econ Entomol 1–8
- 9. Cresswell JE (2011) A meta-analysis of experiments testing the effects of a neonicotinoid insecticide (imidacloprid) on honey bees. Ecotoxicology 20:149–157. pmid:21080222
- 10. Henry M, Béguin M, Requier F, Rollin O, Odoux J-F, Aupinel P, et al. (2012) A common pesticide decreases foraging success and survival in honey bees. Science 336: 348–350. pmid:22461498
- 11. El Hassani AK, Dacher M, Gary V, Lambin M, Gauthier M, Armengaud C (2008) Effects of sublethal doses of acetamiprid and thiamethoxam on the behavior of the Honeybee (Apis mellifera). Arch Environ Contam Toxicol 54:653–661. pmid:18026773
- 12. Yang E-C, Chang H-C, Wu W-Y, Chen Y-W (2012) Impaired olfactory associative behavior of honeybee workers due to contamination of imidacloprid in the larval stage. PLoS ONE 7(11): e49472. pmid:23166680
- 13. Easton AH, Goulson D (2013) The neonicotinoid insecticide imidacloprid repels pollinating flies and beetles at field-realistic concentrations. PLoS ONE 8(1): e54819. pmid:23382980
- 14. Alaux C, Brunet J-L, Dussaubat C, Mondet F, Tchamitchan S, Cousin M, et al. (2010) Interactions between Nosema microspores and a neonicotinoid weaken honeybees (Apis mellifera). Env Microbiol 12: 774–782.
- 15. Collison E, Hird H, Cresswell J, Tyler C (2015) Interactive effects of pesticide exposure and pathogen infection on bee health–a critical analysis. Biol Rev (2015),
- 16. Pettis JS, vanEngelsdorp D, Johnson J, Dively G (2012) Pesticide exposure in honey bees results in increased levels of the gut pathogen Nosema. Naturwissenschaften 99: 153–158. pmid:22246149
- 17. Gill RJ, Ramos-Rodriguez O, Raine NE (2012) Combined pesticide exposure severely affects individual- and colony-level traits in bees. Nature 491: 105–109. pmid:23086150
- 18. Basile R, Pirk CWW, Tautz J (2008) Trophallactic activities in the honeybee brood nest—Heaters get supplied with high performance fuel. Zoology 111: 433–44. pmid:18565745
- 19. Rortais A, Arnold G, Halm M- P, Touffet-Briens F (2005) Modes of honeybees exposure to systemic insecticides: estimated amounts of contaminated pollen and nectar consumed by different categories of bees. Apidologie 36: 71–83
- 20. Meikle W G, Rector BG, Mercadier G, Holst N (2008) Within-day variation in continuous hive weight data as a measure of honey bee colony activity. Apidologie 39: 694–707.
- 21. Meikle WG, Weiss M, Stilwell AR (2016) Monitoring colony phenology using within-day variability in continuous weight and temperature of honey bee hives. Apidologie 47:1–14
- 22. Berry JA, Hood WM, Pietravalle S, Delaplane KS (2013) Field-level sublethal effects of approved bee hive chemicals on honey bees (Apis mellifera L). PLoS ONE 8(10): e76536. pmid:24204638
- 23. Mullin CA, Frazier M, Frazier JL, Ashcraft S, Simonds R, VanEngelsdorp D, et al. (2010) High levels of miticides and agrochemicals in North American apiaries: implications for honey bee health. PLoS ONE 5:e9754. pmid:20333298
Delaplane KS, van der Steen J, Guzman-Novoa E (2013) Standard methods for estimating strength parameters of Apis mellifera colonies. In Dietemann V.,Ellis JD, Neumann P (Eds) The COLOSS BEEBOOK, Volume I: standard methods for Apis mellifera research. J Apic Res 52(1)
- 25. Johnson R (2015) Honey bee toxicology. Annu Rev Entomol 60:415–34. pmid:25341092
Owens CD (1971) The thermology of wintering honey bee colonies. United States Department of Agriculture, Agricultural Research Service, Technical Bulletin no. 1429.
- 27. Pettis JS, Collins AM, Wilbanks R, Feldlaufer MF (2004) Effects of coumaphos on queen rearing in the honey bee, Apis mellifera. Apidologie 35:605–610.
- 28. Kessler SC, Tiedeken EJ, Simcock KL, Derveau S, Mitchell J, Softley S, et al. (2015) Bees prefer foods containing neonicotinoid pesticides. Nature 521 (7550): 74–76. pmid:25901684
- 29. Tomizawa M, Casida JE (2005) Neonicotinoid insecticide toxicology: Mechanisms of selective action. Annu Rev Pharmicol Toxicol 45: 247–268.
- 30. Singaravelan N, Inbar M, Ne'eman G, Distl M, Wink M, Izhaki I (2006) The effects of nectar-nicotine on colony fitness of caged honeybees. J Chem Ecol 32(1): 49–59. pmid:16525869
- 31. Köhler A, Pirk CWW, Nicolson SW (2012) Honeybees and nectar nicotine: Deterrence and reduced survival versus potential health benefits. J Ins Physiol 58(2): 286–292.
- 32. Sandrock C, Tanadini M, Tanadini LG, Fauser-Misslin A, Potts SG, Neumann P (2014) Impact of chronic neonicotinoid exposure on honeybee colony performance and queen supersedure. PLoS One
- 33. Toth AL, Robinson GE (2005) Worker nutrition and division of labour in honeybees. Anim Behav 69: 427–435.
Toth P (2009) M.Sc. Thesis. Lethal and sublethal effects of imidacloprid and amitraz on Apis mellifera (Linneaus) (Hymenoptera: Apidae) larvae and pupae. The University of Florida. Available at http://ufdc.ufl.edu/UFE0024135/00001.
- 35. Stewart SD, Lorenz GM, Catchot AL, Gore J, Cook D, Skinner J, et al. (2014) Potential exposure of pollinators to neonicotinoid insecticides from the use of insecticide seed treatments in the mid-southern United States. Environ Sci Technol 48: 9762−9769. pmid:25010122