Next Article in Journal
Effect of Constant Temperatures on Culicoides sonorensis Midge Physiology and Vesicular Stomatitis Virus Infection
Next Article in Special Issue
Interaction between Thiamethoxam and Deformed Wing Virus Type A on Wing Characteristics and Expression of Immune and Apoptosis Genes in Apis mellifera
Previous Article in Journal
The Use of Tuning Forks for Studying Behavioural Responses in Orb Web Spiders
Previous Article in Special Issue
Sensitivity of Buff-Tailed Bumblebee (Bombus terrestris L.) to Insecticides with Different Mode of Action
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Thiamethoxam-Dressed Oilseed Rape Seeds and Nosema ceranae on Colonies of Apis mellifera iberiensis, L. under Field Conditions of Central Spain. Is Hormesis Playing a Role?

1
Unidad de Productos Fitosanitarios, Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria (INIA, CSIC), 28040 Madrid, Spain
2
Laboratorio de Mieles y Productos de las Colmenas, Centro de Investigación Apícola y Agroambiental (CIAPA), Instituto Regional de Investigación y Desarrollo Agroalimentario y Forestal (IRIAF), Consejería de Agricultura de la Junta de Comunidades de Castilla-La Mancha, 19180 Marchamalo, Spain
3
Experimentación e Investigación Hortícola Centro de Investigación Apícola y Agroambiental, Centro de Investigación Apícola y Agroambiental (CIAPA), Instituto Regional de Investigación y Desarrollo Agroalimentario y Forestal (IRIAF), Consejería de Agricultura de la Junta de Comunidades de Castilla-La Mancha, 19180 Marchamalo, Spain
4
Instituto Universitario CINQUIMA, Analytical Chemistry Group, Universidad de Valladolid, 47011 Valladolid, Spain
5
Instituto de Recursos Humanos para la Ciencia y la Tecnología (INCRECYT–ESF/EC-FSE), Fundación Parque Científico y Tecnológico de Castilla–La Mancha, 02006 Albacete, Spain
6
Laboratorio de Patología Apícola, Centro de Investigación Apícola y Agroambiental (CIAPA), Instituto Regional de Investigación y Desarrollo Agroalimentario y Forestal (IRIAF), Consejería de Agricultura de la Junta de Comunidades de Castilla-La Mancha, 19180 Marchamalo, Spain
*
Author to whom correspondence should be addressed.
Submission received: 17 March 2022 / Revised: 4 April 2022 / Accepted: 6 April 2022 / Published: 9 April 2022

Abstract

:

Simple Summary

The collapse of the honey bee colonies is a complex phenomenon in which different factors may participate in an interrelated manner (e.g., pathogen interactions, exposure to chemicals, beekeeping practices, climatology, etc.). In light of the current debate regarding the interpretation of field and monitoring studies in prospective risk assessments, here we studied how exposure to thiamethoxam affects honey bee colonies in Central Spain when applied as a seed treatment to winter oilseed rape, according to the good agricultural practice in place prior to the EU restrictions. Under the experimental conditions, exposure to thiamethoxam, alone or in combination with other stressors, did not generate and maintain sufficient chronic stress as to provoke honey bee colony collapse. The stress derived from exposure to thiamethoxam and honey bee pathogens was compensated by adjustments in the colony’s dynamics, and by an increase in the worker bee population, a behavior known as hormesis. An analysis of the factors underlying this phenomenon should be incorporated into the prospective risk assessment of plant protection products in order to improve the future interpretation of field studies and management practices.

Abstract

To study the influence of thiamethoxam exposure on colony strength and pathogen prevalence, an apiary (5 colonies) was placed in front of a plot sown with winter oilseed rape (wOSR), just before the flowering phase. Before sowing, the seeds were treated with an equivalent application of 18 g thiamethoxam/ha. For comparison, a second apiary (5 colonies) was located in front of a separate 750 m plot sown with untreated wOSR. Dead foragers at the entrance of hives were assessed every 2–3 days throughout the exposure period, while the colony strength (number of combs covered with adult honey bees and brood) and pathogens were monitored each month until the following spring. Foraging on the wOSR crop was confirmed by melissopalynology determination of the corbicular pollen collected periodically, while the chemical analysis showed that exposure to thiamethoxam was mainly through nectar. There was an increase in the accumulation of dead bees in the apiary exposed to thiamethoxam relating with the control, which was coped with an increment of bee brood surface and adult bee population. However, we did not find statistically significant differences between apiaries (α = 0.05) in terms of the evolution of pathogens. We discuss these results under hormesis perspective.

1. Introduction

In the past decade, significant losses of honey bee colonies have been reported, mainly in Europe and North America [1,2,3,4,5,6,7,8,9,10,11]. The consensus is that several factors are involved in this phenomenon: nutritional stress, pathogens, hive management practices, exposure to multiple xenobiotic residues (pesticides and veterinary drugs), invasive species, seasonal weather changes and genetic variability [7,12,13,14,15,16,17,18,19]. However, exposure to chemicals and the most prevalent pathogens are without doubt the main drivers in the phenomenon of honey bee colony loss [20].
Among the most relevant chemicals, much attention has been paid to neonicotinoids [21] in recent years due to their widespread use [22,23,24,25]. These pesticides selectively act on the insect’s central nervous system (CNS) as they are agonists of post-synaptic acetylcholine receptors [26,27,28,29]. Neonicotinoids have a low molecular weight, they are moderately-to-highly water soluble and they have a low octanol:water partitioning coefficient. These physico-chemical properties favor the efficient translocation of these chemicals through the plant xylem, such that they can systemically control piercing-sucking insects [24]. However, these characteristics also mean they have a high capacity for environmental transport [30,31,32]. Depending on the toxophore, these compounds are classified as N-nitroguanidines (imidacloprid, thiamethoxam, clothianidin and dinotefuran), nitromethylenes (nitenpyram), or N-cyano-amidines (acetamiprid and thiacloprid) neonicotinoids, with the former representing the class that has produced most concern for bees [33]. Many sub-lethal effects of N-nitroguanidine neonicotinoids have been identified in bees, influencing on orientation and olfactory learning, flight, queen performance, honey bee colony physiology and the bee’s immune system [34,35,36,37,38,39,40,41,42,43,44,45,46,47]. In particular, thiamethoxam (TMX) causes damage in the midgut, brain cells and Malpighian tubules of bees [48,49,50], affecting several biological process such as oxidative phosphorylation and tyrosine metabolism [51], as well as altering the gut microflora [52].
On the other hand, Nosema ceranae and Varroa destructor have been intensively studied in relation to the phenomenon of honey bee colony losses. N. ceranae is found worldwide [14,53,54], it has a high pathogenic capacity [55,56] and it is transmitted horizontally in the colony, through food, water, trophallaxis and cleaning duties (reviewed in [54]). It is an obligate intracellular parasite that depends on the host’s energy for reproduction, and it alters the physiology and behavior of its host to ensure conditions are optimal for it to complete its biological cycle (reviewed in [54]). As such, Apis mellifera individuals infected by N. ceranae experience: (1) a strong carbohydrate demand; (2) altered lipid and amino acid metabolism; (3) an altered vitellogenin/juvenile hormone equilibrium; (4) a reduction in hexamerins and major royal jelly proteins; (5) an upregulation of the octampamide pathway; (6) a suppressed immune response; (7) impaired apoptosis of ventricular epithelial cells [54,57,58], and (8) changes in midgut microbiota [52]. As a consequence of all these alterations, N. ceranae will perturb colony homeostasis if the infection is maintained over time [59]. In this way, infected bees generate higher Ethyl Olate (EO) titers and undertake more flight activity than non-infected bees [60]. The constant reduction in bee number due to the tissue damage, energetic stress, and altered flight behavior due to infection accelerates age polyethism in young bees to cover the energy demands of the colony [61,62]. However, precocious foragers (FBs) are less effective, so higher number of FBs will be needed to fulfil the colony’s demands, reducing the time each bee can dedicate to hive tasks. Thus, if the buffering capacity of the colony is exceeded then there is a strong risk of collapse [54].
For its part, V. destructor, one of the most important pathogens of honey bees [63,64,65], not only damages the bee by feeding on its hemolymph and fat body [66,67] but also, it acts as a vector for a number of viruses [68]. As a result, this pathogen impairs the responses of the immune system [69] and affects the nervous system [70,71].
Honey bees exposed to neonicotinoids may be more susceptible to N. ceranae and V. destructor infections and that of other viruses [72]. Indeed, N. ceranae-infected honey bees suffer increased mortality under laboratory conditions when exposed to sub-lethal concentrations of imidacloprid [73,74], thiacloprid [75,76] or TMX [77]. Similarly, the damage caused by N. ceranae in midgut cells [78] may make infected honey bees more susceptible to neonicotinoid exposure. The negative synergistic effects of Varroa and neonicotinoid insecticides have been related to immune suppression, which impairs antiviral immune barriers [47] and memory retention [79].
However, as all these effects are rarely observed under real world conditions [80,81], they remain somewhat controversial.
Here, we present the results of a field study undertaken in accordance with good beekeeping practice (prior to European restrictions) in Central Spain that set out to determine the evolution of honey bee colonies and the prevalence of pathogens following exposure to TMX through the pollen and nectar of seed-treated oilseed rape. The stress provoked by both TMX and pathogens was compensated by the dynamics of the colonies, mainly through an increase in the bee population. We discuss the possible reasons for these responses from the hormesis point of view and in the context of the prospective risk assessment of plant protection protocols.

2. Materials and Methods

2.1. Experimental Design

The experiments were designed in accordance with the EPPO 170 standards [82]. Colonies of Apis mellifera iberiensis from the same apiary located at the ‘Centro de Investigación Apícola y Agroambiental” (CIAPA, Marchamalo, Guadalajara, Spain) and with naturally mated sister queens of the same age, were used to monitor the presence of queens, and to analyze and control the sanitary status of the colonies in order to ensure the colony’s strength at the beginning of the study (10 brood chamber combs covered by bees). The experimental colonies were randomly divided into two groups of five colonies each, and they were placed in front of two plots of 2 ha (750 m apart), that were sown with winter oilseed rape (wOSR, Ginfizz variety), just before their flowering phase (Figure 1).
Sowing took place on 29 October 2014 at a density of 75 seeds m−2. In one plot, the seeds were treated with Cruiser 350 FS© before sowing (TMX 35% p/v FS, 1200 mL/100 kg seed), while both fields were treated with Butisan S© (metazachlor 50% p/v SC) at 1.5 L/ha on 31 November, receiving no further phytosanitary treatments during the study. The fields were not treated with plant protection products for at least 2 years and before sowing, the soil of the plots was sampled and characterized as described in Section 2.2.
The exposure phase lasted from 8 April to 11 May 2015, at the beginning of which an upper box was placed in each hive, which contained empty combs ready to be used by the honey bees. The stored pollen (beebread) and honey in the brood chambers was removed to stimulate foraging activity during the exposure period, and to ensure that the pollen and nectar came from the experimental plots [83]. After the exposure phase, the colonies were placed 35 km away from the area of exposure in order to monitor them until the following spring (Figure 1).
To avoid winter collapse due to nosemosis C, the colonies were treated exceptionally with Fumidil B on 17 September 2015 under controlled conditions. The treatment consisted of 4 doses of 30 mg fumagillin per colony dissolved in 250 cc of sugar syrup (50% sugar/distilled water), administered in plastic bags placed over the brood at one-week intervals to ensure full consumption. This posology was effective without leaving residues in the honey [84]. V. destructor was controlled at the beginning of September with Apitraz© (amitraz a.m.), employed at the recommended dose and in accordance with current Spanish legislation [85].

2.2. Soil Characterization and Climatological Data

The soil surface layer (0–20 cm) was characterized following ISRIC protocols [86] and using three samples gathered randomly in each plot. The anion analysis (F, Cl, NO2, NO3, PO43−, SO42−) was carried out at room temperature using a Metrohm Compact IC model 761 ion chromatograph (Metrohm, Herisau, Switzerland). Na+ and K+ were measured using a flame photometry SHERWOOD 410 apparatus, and Ca2+ and Mg2+ using atomic absorption spectrophotometry (ANALYTIKJENA NOVAA 300). Electrical conductivity (EC), pH, and humidity were measured with an EC meter CRISON micro CM2000, a pH meter CRISON GLP 21, and a thermogravimetric balance KERN DBS, respectively. Bouyoucos’ densitometry method was followed to characterize the soil texture [87], and the organic matter and carbonate content was determined with the Walkley–Black and acid neutralization methods [86], respectively.
A Walter–Lieth diagram [88] was developed with historical weather data obtained from the meteorological station at Marchamalo from the regional Irrigation Advisory Service (SIAR-CREA) network to compare it with data during the period of study. These diagrams illustrate the rainfall and temperature changes throughout the year in standardized charts which provide brief summaries of the average climatic variables and their time course. The diagrams were drawn up with the diagnostic tool of the Worldwide Bioclimatic Classification System, 1996–2021 [89].

2.3. Sampling Schedule

The number of dead worker bees was counted in box dead-bee traps [90] at the entrance of hives every 2–3 days throughout the exposure period. Colony strength was checked by determining the number of combs covered by brood and bees once every two months [59]. The honey production in each colony was evaluated in the harvesting season as the difference in the weight of the combs before and after honey extraction [91].
Corbicular pollen loads were collected using standard traps, activated for 24 h at the beginning (8 April), in the middle (15 April), and at the end of blooming (29 April). For melissopalynological and chemical analyses, at the end of the exposure period, beebread and honey samples were collected individually from one frame of the brood chambers and the upper boxes, respectively. Beebread was extracted aseptically, removing the wax from the combs as described previously [92,93,94], and storing the samples at −80 °C in the dark until they were analyzed. The honey was sampled using a sterile spatula and stored at −20 °C.
To identify viruses and quantify the Nosema spp. parasitic load in the colonies, both non-foragers, (IBs, bees sampled from no-brood combs, n > 60), and foragers (FBs, bees arriving at the colony after closing the hive entrance for 30 min, n > 60) were collected separately from each colony (Table 1) and frozen (−20 °C).

2.4. Palynological and Melissopalynological Assessment

To confirm the type of foraging flora, corbicular pollen, beebread and honey samples were analyzed according to previously described methods [95,96]. Corbicular pollen samples were weighed and stored at −20 ± 2 °C until further analysis. Pollen loads were separated based on color and texture to identify and determine the contribution of each botanical type into the samples. The pollen grains were isolated from each sample and cleaned using the Erdtam method [97]. One aliquot from each of the ten pollen loads of the same “color-texture” class was placed onto a slide in glycerin jelly mounting medium and examined microscopically in order to ensure the homogeneity of identification. The proportion of each pollen type was calculated by dividing the weight obtained for each color fraction by the overall weight of the sample analyzed [95].
Pollen was extracted from the beebread by diluting 0.5 g in 10 mL of 0.5% acidulated water (96% sulfuric acid) and centrifuging at 2500 r.p.m. for 15 min. The pellet was washed with double-distilled water and centrifuged twice, and the sediment was placed onto a slide in glycerin jelly mounting medium and examined microscopically to identify the pollen. The frequency of Brassica type pollen grains (Brassica t. pollen) was expressed as a percentage of the total pollen grains.
The honey samples were treated chemically with 0.5% acidified water (sulfuric acid 96%), and a qualitative and quantitative analysis was performed on the sediment recovered from 5 g aliquots [98]. Between 300 and 1200 pollen grains were counted in each sample, and the pollen grains were identified and classified on the basis of the identification keys [99,100]. All the reference collection pollen slides used were available at the honey laboratory at the CIAPA.

2.5. Chemical Analysis

Chemical analyses were carried out as described elsewhere [101,102,103]. FLUKA-PESTANAL analytical standards of clothianidin (CLO; purity 99.9%), TMX (purity 99.6%), and TMX-d3 (isotope-labeled standard; purity ≥ 98.0%) were purchased from Sigma-Aldrich Laborchemikalien (Seelze, Germany). LC grade methanol (MeOH), and acetonitrile (ACN) were supplied by Lab Scan Ltd. (Dublin, Ireland). Formic acid (FA), ammonium formate, and ethyl acetate (EA) were obtained from Sigma-Aldrich Chemie Gbmh (Steinheim, Germany); sodium chloride, trisodium citrate dihydrate, and trisodium citrate sesquihydrate were supplied by Panreac (Barcelona, Spain). Meanwhile, primary secondary amine (PSA) and octadecylsilane (C18) were provided by Supelco (Bellefonte, PA, USA). Finally, ultrapure water was obtained using Milipore Mili-RO plus and Mili-Q systems (Bedford, MA, USA). All chemicals used were of analytical grade.
Briefly, to quantify the levels of TMX and its metabolite clothianidin (CLO) in the seeds [101], insecticides were extracted with a mixture of acetonitrile (ACN) and sodium chloride (60:40, v/v). After centrifugation (5810 R refrigerated bench-top Eppendorf centrifuge; Hamburg, Germany), the supernatant was collected and concentrated at 60 °C (R-210/215 rotary evaporator from Buchi, Flawil, Switzerland). The dry extract was reconstituted with 1 mL of a mixture of 0.1% (v/v) FA in ACN and 0.1% (v/v) FA in water (25:75, v/v). The extract was filtered and injected onto a liquid chromatography coupled to a diode array detector (DAD; Agilent Technologies 1200 series; Palo Alto, Santa Clara, CA, USA) column (Kinetex C18, 150 × 4.6 mm, 2.6 μm, 100 Å; Phenomenex, Torrance, CA, USA) system. The limits of quantifications (LOQs) were set at 0.11 and 0.15 g/kg for TMX and CLO, respectively [101].
In relation to the analysis in bee pollen, it was employed as a sample treatment modified QuEChERS (quick, easy, cheap, effective, rugged, and safe) method. A representative amount of sample (1.0 g) was mixed with 2 mL of water and 6 mL of ACN, and the resulting mixture was shaken. Then, magnesium sulfate (1.0 g), sodium chloride (0.5 g), and tri-sodium citrate dihydrate (0.8 g) was added to the mixture, and after centrifugation (10,000 rpm, 10 °C, 5 min), 2 mL of the supernatant was transferred to an Eppendorf tube. After that, magnesium sulfate (150 mg), PSA (25 mg), and C18 (25 mg) were added to the tube, which was again centrifuged. Supernatant (1 mL) was collected and evaporated to dryness (60 °C). A further reconstitution of the dry extract was performed with 1 mL of a MeOH and water (80:20, v/v) mixture, and the resulting extract was filtered and analyzed by using an ultra-high performance liquid chromatography (UHPLC; ACQUITY, Waters, Milford, MA, USA) coupled to quadrupole time-of-flight mass spectrometry (qTOF; maXis impact, Bruker Daltonik GmbH, Bremen, Germany) system equipped with an electrospray interface (ESI), which was operated in positive mode. The obtained LOQs were 2.1 ng/g for TMX and 3.9 ng/g for CLO [102].
Finally, two different sample treatments were employed when determining TMX and CLO depending on the botanical origin of honey. Insecticides were extracted from light color honeys (multifloral and rosemary) by using a modified QuEChERS protocol. Briefly, 5.0 g of honey was mixed with 10 mL of water and 10 mL of an ACN and EA (70:30, v/v) mixture, and after shaking, different amounts of magnesium sulfate (2.0 g), sodium acetate (1.0 g), trisodium citrate dihydrate (1.5 g), and trisodium citrate sesquihydrate were added. Then, the mixture was shaken in an ultrasound device (J.P. Selecta S.A., Barcelona, Spain), centrifuged (10,000 rpm, 10 °C, 5 min); the supernatant was collected and evaporated to dryness (60 °C), while the dry extract was reconstituted with 1 mL of a MeOH and water (80:20, v/v) mixture, being the resulting extract filtered and analyzed by an ACQUITY UHPLC (Waters) coupled to a Xevo TQ-S (triple quadrupole, QqQ) mass spectrometer (Waters) equipped with an ESI interface (positive mode). The obtained LOQs were 0.06 ng/g for TMX and 0.20 ng/g for CLO [102]. On the other hand, dark honeys (heather) required a solid-phase extraction (SPE) procedure. A representative amount of dark honey (5.0 g) was diluted in 10 mL of ammonium formate (10 mmol/L) in water. The resulting mixture was loaded onto a polymeric SPE cartridge (Strata® X; Phenomenex), which was previously conditioned with 5 mL of MeOH and 5 mL of water. Then, after 5 min of drying time, the insecticides were eluted from the cartridge with 4 mL of an ACN and EA (80:20, v/v) mixture. The extract was evaporated to dryness, and the dry residue reconstituted with 1 mL of a MeOH and water (80:20, v/v) mixture, filtered and analyzed by using the same UHPLC-QqQ system that was employed for analyzing light honeys. The LOQs were 0.06 ng/g for TMX and 0.20 ng/g for CLO [103].

2.6. Identification of Pathogens

2.6.1. Detection of Varroa Mites

At each sampling time, a total of 60 FBs and 60 IBs were examined individually per colony, using sterile tweezers to detect malformations and collect parasitic mites to be identified macroscopically. A honey bee colony was considered to be infested with V. destructor when at least 1 Varroa mite was found in the sample. The rate of infestation of the bee colony was estimated by assessing the number of Varroa mites relative to the number of adult bees in each sample, and it was expressed as the number of Varroa mites/100 bees/sample.

2.6.2. Pathogen Screening

DNA and RNA were extracted at each sampling time (see Table 1) as described previously [104]. In summary, the abdomen was removed from 30 FB and 30 IB bees in each colony, placing them individually in a well of a 96-well plate containing glass beads (2 mm: Sigma®, St. Louis, MO, USA). The tissue was macerated in AL buffer (Qiagen® 19075, Hilden, Germany) and incubated with proteinase K under the same conditions as described in [104]. DNA and RNA purification were performed simultaneously in a Biosprint 96 workstation (Quiagen®, Hilden, Germany) according to the BS96 DNA Tissue extraction protocol. Total nucleic acids (RNA and DNA) were eluted in 100 µL of elution buffer.
An aliquot of 75 µL was used to detect Nosema spp. according to an earlier procedure [105]. To identify N. ceranae and N. apis, an internal control (IPC) based on mitochondrial cytochrome oxidase subunit I (COI) was used as a target in the triplex PCR assay [105]. The percentage of FBs or IBs parasitized by Nosema was estimated as 100 × number of positives/30.
The remaining 25 µL aliquots of the total nucleic acids eluted from each FB and IB sample from each colony and sampling time were pooled immediately and subjected to DNase I digestion to remove DNA (Quiagen® kit 79254), as described previously [104]. The total RNA recovered was used immediately to generate first strand cDNA with the iScriptTM cDNA Synthesis Kit (Biorad®, Hercules, CA, USA). The resultant cDNA was used to analyze the viruses present. Black Queen Cell Virus (BQCV) and Deformed Wings Virus (DWV) were detected by Real Time-PCR as detailed in [104], and using the primers and probes described by [106]. The Acute Bee Paralysis Virus-Kashmir Bee Virus-Israeli Acute Paralysis Virus (AKI) complex was detected by following the procedure described in [107]. In all cases, negative and positive controls were run in parallel for each step (sample processing, nucleic acid extraction, reverse transcription, and RT-PCR and PCR analysis).

2.7. Statistical Analysis

The statistical analysis was made according to the recommendations in [108,109]. Colony strength or pathogen parameters were compared between the two groups of colonies at the end of the exposure period (n(control) = n(TMX) = 5) with a one tailed Mann–Whitney U test (exact approach, alpha = 0.05). The Mann–Whitney U test compares the distribution of the ranks of the two populations, and to interpret this as a comparison of medians, the distributions of the two populations are assumed to be the same shape [110]. Moreover, we reported the effect size for non-parametric tests following the approach given elsewhere [111]. Mann–Whitney U tests and effect sizes were generated using the Real Statistics Resource Pack software [112] and transformed to the Cohen’s d coefficient as described previously [113]. Finally, we used G-Power 3.1.9.6 [114] to estimate the power (1-β) of the statistical tests.
To explore the evolution of the percentage of Nosema spp. infection in each honey bee colony, we used a Friedman test. This test is the non-parametric alternative to one-way ANOVA with repeated measures [115] and it was implemented with Statgraphics Centurion 18.

3. Results

3.1. Soil Characterization and Climatic Conditions

The 0–20 cm soil horizons of the control and treatment plots were classified as slightly acid clay and neutral clay loam, respectively, in both cases with low organic matter content and soil conductivity. The chemical characterization of the soil showed the CEC is dominated by Ca2+ (Table S1). The soils are derived from Miocene-tertiary aged fluvio-alluvial sediments from Henares River Basin [116], and following the classification criteria of the World Reference Base (WRB) for Soil Resources [117], they could tentatively be classified as calcic luvisols.
Based on 12 years of climatological data (period 2010–2021) from the nearest station to the experimental plots (Figure S1a), and on the bioclimatic indices estimated with the diagnostic tool of the Worldwide Bioclimatic Classification System [89], the climate in the experimental plots falls into the Upper Mesomediterranean Low Dry Blioclimatic belt, with a freezing period extending from December to March. During the period of exposure, the mean temperature was 13.0 ± 2.0 °C in April and 17.8 ± 2.5 °C in May 2015. The relative humidity ranged from 38.8 to 92.9 (mean = 61.58 ± 13.33) and the total rainfall was 29.3 mm, concentrated in 9 days and with 52.56% of the precipitation falling on 26 April 2015.
When compared to historic average climate data (Figure S1a), the Walter–Lieth diagrams show that autumn 2014 (Figure S1b) was wetter than the mean values for the period between 2010 and 2021. However, 2015 was characterized as being dryer and hotter than average, and with a wider range of temperatures (Figure S1c). By contrast, spring 2016 was more temperate and significantly wetter, especially in April (Figure S1d).

3.2. Palynological and Melissopalynological Assessments

The Brassica t. pollen content varied in the corbicular pollen sampled in April (Figure 2a). In addition to Brassica t. pollen, Prunus t. and Cistus t. pollens were also predominant in the samples at the beginning of the month. These taxa are typical of regression succession steps of the chorological series Bupleuro rigidi-Querceto rotundifoliae sigmetum in the Castellano-Maestrazgo-Manchega province [118]. In the middle of April, Brassica t. pollen was replaced by pollen from the surroundings, mainly Cistus t pollen, which was in turn replaced by Brassica t. pollen at the end of the month. At the end of the exposure period, the Brassica t. pollen content in the beebread was above 50% in all the colonies except two (C-01 and C-08), and the Brassica t. content of the honey in the cells of the frame introduced at the beginning of the exposure period ranged from 50.3% to 94.2% (Figure 2b).

3.3. Chemical Analysis

The chemical analyses revealed the level of TMX residues in treated seeds was 4660 mg/kg seed, equivalent to an application of 18 g TMX/ha. The TMX residues in bee pollen (corbicular pollen) were below the LOQ in all samples, while in the honey samples, it varied between below the limit of detection (<LOD) and 144 µg/kg (Table S13 in [103]).

3.4. Field Monitoring of the Honey Bee Colonies

During the exposure period, the accumulated honey bee death was higher in the treatment apiary (TMX) as the combs became covered with worker honey bees and brood (Figure 2, Table 2). Subsequently, both the worker honey bee population and the brood balanced out in both apiaries (Figure 3).
The mean number of dead bees was 127.6 ± 44.1 and 214.4 ± 116.5 in the control and treatment (TMX) apiaries, respectively. Likewise, the average number of combs covered with adult honey bees and honey bee brood was 12.2 ± 1.095 and 9 ± 4.96 in the controls, respectively, and 15.2 ± 2.05 and 11 ± 1.82 in the treatment apiaries. The same trend was observed in terms of production over the period of exposure, reaching 24.75 ± 15.05 (control) and 31.98 ± 12.21 kg (TMX) at the end of the exposure phase, respectively (Table 2).
Due to the high dispersion of the data within the apiaries (Table 2), at the end of the exposure period, a one-tailed Mann–Whitney U test did not find significant differences in the dead bees accumulated (Mann–Whitney U = 6; p-value = 0.111), the number of combs with bee brood (U = 6.5; p-value = 0.343) or in the honey produced (U = 9; p-value = 0.274). As such, the statistical power varied between 0.183 and 0.397 (Table 2). Regarding the number of combs covered with bees, the test did reveal a significant difference (U = 4.5; p-value = 0.048; 1-β = 0.72).

3.5. Identification of Pathogens

While N. apis was not detected in the colonies at any point in the study, the percentage of infection of N. ceranae increased in both non-forager (IB) and forager bee (FB) populations until fumagillin treatment was applied (Figure 4). The percentage of infection was significantly different between the distinct honey bee populations (Friedman Test statistic = 71.4394; p-value = 1.132 10–14) but not between apiaries. The maximum percentage infection was observed at the end of summer, with a mean of 62.6 ± 10.6% of FBs infected and a mean of 17.6 ± 5.06% of IBs.
Varroa mite was not detected in any honey bee sampled and DWV was found at a low frequency until the end of the summer, yet it was detected in all the colonies just before wintering (Table 3). By contrast, BQCV was detected throughout the experiment (Table 3), whereas the AKI complex was not identified in any sample.

4. Discussion

The foraging activity of honey bees in the experimental fields was guaranteed during the period of exposure due to the melliferous potential of the oilseed rape crop. Its attractiveness to honey bees is derived from the features of its flowers, their fragrance and nutritive nectar, and the high flower density of its plants, opening successively for 3–4 weeks [119]. The amount of Brassica t. in the bee pollen varied throughout the month of April, decreasing in the middle of April when the bulk flowering of other species took place, such as the plants of the genus Cistus [120,121], a producer of high quality pollen [122,123]. The chemical analysis confirmed that the honey bees were exposed to TMX through nectar rather than pollen [103]. The levels of residues in wOSR plants may have been influenced by the time between the sowing (Autumn 2014) and the exposure time (Spring 2015), and possibly due to soil leaching. Given their physico-chemical properties [24], TMX and its metabolite chlothianidin have high potential for lixiviation. This fact, together with the precipitation that fell in the days following drilling and the low organic matter content of the soils, may have favored the low retention of the residues of these chemicals in soils.
Several detections above LOQ were found in the honey samples of the control apiary [103], however, this does not invalidate the results. Thus, an increment on the mean accumulated deaths in the apiary located in front of the treated plot was observed when compared with the control at the end of the exposure period. Taking into account that exposure during foraging is not a discrete process but distributional [124], the mean colony exposure remained low even though effects were seen in a number of bees. The loss of forager bees can accelerate behavioral changes aimed at increasing food gathering in order to maintain the energy balance within the colony [125]. These changes bring about an increase in bee brood to compensate the interactions within the colony [126], referred to as a hormetic-like effect [127,128,129,130]. Hormesis is common in arthropods [131,132] and it is a dose-dependent phenomenon where exposure to high levels of a stressor are inhibitory, whereas low doses are stimulatory [133]. Such effects have been considered to be an example of evolutionary fitness [134] whereby organisms undergo metabolic adaptations in response to the interactions between the diverse stressors they are exposed to. Indeed, it has been postulated that such responses may be driven by antioxidants and oxidative stress [135]. Chronic stress on the colony has an important influence on colony dynamics [136], and if it surpasses a threshold, there is a risk of a destructive response, in this case, derived from the less efficient foraging of young bees [137]. This imbalance in the colony can drive the queen towards the limits of her egg-laying capacity [138], which may be accelerated by the interaction of several stressors. This behavior is similar to the Allee effect, which is related to the extinction of populations [139,140,141].
In terms of neonicotinoids, a significant increase in capped bee brood has been documented recently in small nucleus colonies exposed to low concentrations of chlothianidin over a period of 7 weeks [126]. It was reported that colonies needed to produce 1.57-fold more larvae to maintain a stable population when exposed to a concentration of 1 µg chlothianidin/L syrup, and that this effort increases as does exposure. Alternatively, effects that compensate nosemosis C [59] were evident in colonies during phase 2 of the disease, but if this does not control the stress provoked by the parasite in phase 3, it can eventually lead to colony collapse. Under the conditions of our field study, N. ceranae infection developed similarly over time in the FB and IB populations in both apiaries, with a relative infection below 5% at the end of the exposure period. The infection in the colonies had evolved to phase 2 (replacement) of nosemosis C by the end of the summer season [59]. Similar behavior has been described in colonies affected by varroosis, which might fail to survive because of their effort to redress the disturbed homeostasis [81]. Nevertheless, the pressure exerted by V. destructor was sufficiently well-controlled in the present study through the mandatory autumn treatment. In fact, according to the Spanish Ministry of Agriculture’s monitoring program, the prevalence of V. destructor was <1% in more than 80% of the colonies sampled in Castilla la Mancha in autumn 2014 and spring 2015 [142].
Immune defenses are costly to the individual, entering into a trade-off with other life-history traits such as reproduction, growth, and self-maintenance [134,143]. Suppression of the immune system may be favored by the interaction of several stressors, such as N. ceranae and TMX [77,144], and by high environmental temperatures [145], which could explain the prevalence of the BQCV virus in the colonies throughout the study. Significant increases in DWV loads have been documented in honey bees exposed to clothianidin, the metabolite of TMX, which has been related to inhibition of dorsal-1a gene transcription [36]. However, this increment was not observed here, with an increase in DWV only apparent in the autumn. These results are consistent with previous data [146] and they may be related to potential routes of transmission other than the main one through Varroa mite [147,148]. Moreover, the high environmental temperature in the spring and summer of 2015, together with a possible increase in thermogenesis due to exposure to TMX, might influence the inhibition of DWV reproduction, in part due to the host synthesis of heat shock proteins [149,150,151,152,153] and possibly by exceeding the optimal growth temperature range of this virus [154]. Nevertheless, further research is needed to understand the influence of high temperatures on the complex relationships between different traits of honey bee colonies and their pathogens.
Hormetic effects may not be statistically significant at the field level, in part because of their low magnitude and in part because of inadequate replication in experiments [155]. However, the difference between statistical significance and the biological relevance of an effect are not necessarily linked. Significance is a statistical measure reflecting whether an observed effect is likely to have occurred by chance alone [156]. EPPO standards (170) mean that statistical analysis may not be feasible in field studies due to the inherent variability of the end-points assessed and the limitations on replication [82]. Based on a Cohen’s d of 0.310, the lowest of our experimental results, it would be necessary to analyze 99 colonies/treatment to increase the statistical power up to 80%. Considering the possibility of establishing 8 hives/plot, this means a total of 12–13 plots per treatment. With this design, and according to the analysis carried out recently by the European Food Safety Authority [157], a significant effect between 7–8% could be detected. Few field studies in the literature follow such an exhaustive study design and they have generally failed to find deleterious effects at a colony level as a result of the use of neonicotinoids under good agricultural practices [158,159,160,161,162,163]. However, an increase in colony size is not always found as a result of neonicotinoid use, as observed here, suggesting that environmental conditions could play a key role in colony behavior.

5. Conclusions

Hormesis is common in arthropods and should be integrated into the dose response model as a continuum. However, the cost of the compensatory behavior adopted by the honey bee colonies over the long term is difficult to assess from an ecological point of view due to the complex relationships between the different traits affected. Further research is needed to understand the mechanisms underlining these compensation events, and the influence of environmental conditions on them, in order to define more appropriate management programs for this key pollinator.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/insects13040371/s1, Table S1: Characteristics of the soil at the experimental plots (mean ± sd, n = 3), Figure S1: Walter–Lieth diagrams for the period 2010–2021 (a), 2014 (b), 2015 (c), and 2016 (d), data taken from the Marchamalo meteorological station (40°40′ N 3°12′ W).

Author Contributions

Conceptualization, M.H., E.A.-P. and R.M.-H.; data curation, M.H. and E.A.-P.; formal analysis: E.A.-P.; funding acquisition, M.H. and J.B.; investigation, E.A.-P., A.V.G.-P., C.G.-V., J.A.L.-P., S.V., J.B., R.M.-H. and M.H.; methodology, E.A.-P., A.V.G.-P., C.G.-V., J.A.L.-P., S.V., J.B., R.M.-H. and M.H.; project administration, M.H. and J.B. resources, E.A.-P., A.V.G.-P., C.G.-V., J.A.L.-P., S.V., J.B., R.M.-H. and M.H; supervision, E.A.-P. and M.H.; validation, M.H., E.A.-P., R.M.-H., A.V.G.-P. and J.B. writing—original draft preparation, E.A.-P., A.V.G.-P., C.G.-V., J.A.L.-P., S.V., J.B., R.M.-H. and M.H.; writing—review and editing, E.A.-P., A.V.G.-P., C.G.-V., J.A.L.-P., S.V., J.B., R.M.-H. and M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by INIA-FEDER funds (RTA2013-00042-C10-03 and RTA2013-00042-C10-06).

Institutional Review Board Statement

The EU Directive 2010/63/EU on the protection of animals used for scientific purposes laid down the ethical framework for the use of animals in scientific experiments. The scope of this directive also includes specific invertebrate species (e.g., cephalopods) but not insects. Thus, according to European Legislation, no specific permits were required for the studies described.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on reasonable request from the corresponding author.

Acknowledgments

The authors wish to thank V. Albendea, T. Corrales, M. Gajero, and C. Uceta at the honey bee pathology laboratory for their technical support, and L. Barrios at Departamento de Estadística, CTI, CSIC for the revision and recommendations.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Brodschneider, R.; Gray, A.; Adjlane, N.; Ballis, A.; Brusbardis, V.; Charrière, J.-D.; Chlebo, R.; Coffey, M.F.; Dahle, B.; de Graaf, D.C.; et al. Multi-country loss rates of honey bee colonies during winter 2016/2017 from the COLOSS survey. J. Apic. Res. 2018, 57, 452–457. [Google Scholar] [CrossRef] [Green Version]
  2. Brodschneider, R.; Gray, A.; van der Zee, R.; Adjlane, N.; Brusbardis, V.; Charrière, J.-D.; Chlebo, R.; Coffey, M.F.; Crailsheim, K.; Dahle, B.; et al. Preliminary analysis of loss rates of honey bee colonies during winter 2015/16 from the COLOSS survey. J. Apic. Res. 2016, 55, 375–378. [Google Scholar] [CrossRef] [Green Version]
  3. Gray, A.; Adjlane, N.; Arab, A.; Ballis, A.; Brusbardis, V.; Charrière, J.-D.; Chlebo, R.; Coffey, M.F.; Cornelissen, B.; da Costa, C.A.; et al. Honey bee colony winter loss rates for 35 countries participating in the COLOSS survey for winter 2018–2019, and the effects of a new queen on the risk of colony winter loss. J. Apic. Res. 2020, 59, 744–751. [Google Scholar] [CrossRef]
  4. Gray, A.; Brodschneider, R.; Adjlane, N.; Ballis, A.; Brusbardis, V.; Charrière, J.-D.; Chlebo, R.; Coffey, M.F.; Cornelissen, B.; Amaro da Costa, C.; et al. Loss rates of honey bee colonies during winter 2017/18 in 36 countries participating in the COLOSS survey, including effects of forage sources. J. Apic. Res. 2019, 58, 479–485. [Google Scholar] [CrossRef] [Green Version]
  5. Lee, K.V.; Steinhauer, N.; Rennich, K.; Wilson, M.E.; Tarpy, D.R.; Caron, D.M.; Rose, R.; Delaplane, K.S.; Baylis, K.; Lengerich, E.J.; et al. A national survey of managed honey bee 2013–2014 annual colony losses in the USA. Apidologie 2015, 46, 292–305. [Google Scholar] [CrossRef] [Green Version]
  6. Neumann, P.; Carreck, N.L. Honey bee colony losses. J. Apic. Res. 2010, 49, 1–6. [Google Scholar] [CrossRef] [Green Version]
  7. Van der Zee, R.; Gray, A.; Pisa, L.; de Rijk, T. An Observational Study of Honey Bee Colony Winter Losses and Their Association with Varroa destructor, Neonicotinoids and Other Risk Factors. PLoS ONE 2015, 10, e0131611. [Google Scholar] [CrossRef] [Green Version]
  8. Van der Zee, R.; Pisa, L.; Andonov, S.; Brodschneider, R.; Charrière, J.-D.; Chlebo, R.; Coffey, M.F.; Crailsheim, K.; Dahle, B.; Gajda, A.; et al. Managed honey bee colony losses in Canada, China, Europe, Israel and Turkey, for the winters of 2008–2009 and 2009–2010. J. Apic. Res. 2012, 51, 100–114. [Google Scholar] [CrossRef]
  9. VanEngelsdorp, D.; Caron, D.; Hayes, J.; Underwood, R.; Henson, M.; Rennich, K.; Spleen, A.; Andree, M.; Snyder, R.; Lee, K.; et al. A national survey of managed honey bee 2010–11 winter colony losses in the USA: Results from the Bee Informed Partnership. J. Apic. Res. 2012, 51, 115–124. [Google Scholar] [CrossRef] [Green Version]
  10. VanEngelsdorp, D.; Hayes, J., Jr.; Underwood, R.M.; Pettis, J. A Survey of Honey Bee Colony Losses in the U.S., Fall 2007 to Spring 2008. PLoS ONE 2009, 3, e4071. [Google Scholar] [CrossRef]
  11. VanEngelsdorp, D.; Hayes, J.; Underwood, R.M.; Pettis, J.S. A survey of honey bee colony losses in the United States, Fall 2008 to Spring 2009. J. Apic. Res. 2010, 49, 7–14. [Google Scholar] [CrossRef] [Green Version]
  12. Carreck, N.L. A beekeeper’s perspective on the neonicotinoid ban. Pest Manag. Sci. 2017, 73, 1295–1298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Clermont, A.; Eickermann, M.; Kraus, F.; Georges, C.; Hoffmann, L.; Beyer, M. A survey on some factors potentially affecting losses of managed honey bee colonies in Luxembourg over the winters 2010/2011 and 2011/2012. J. Apic. Res. 2014, 53, 43–56. [Google Scholar] [CrossRef]
  14. Goulson, D.; Nicholls, E.; Botías, C.; Rotheray, E.L. Bee declines driven by combined stress from parasites, pesticides, and lack of flowers. Science 2015, 347, 1255957. [Google Scholar] [CrossRef] [PubMed]
  15. Hristov, P.; Shumkova, R.; Palova, N.; Neov, B. Factors Associated with Honey Bee Colony Losses: A Mini-Review. Vet. Sci. 2020, 7, 166. [Google Scholar] [CrossRef]
  16. Johnson, R.M. Honey Bee Toxicology. Annu. Rev. Entomol. 2015, 60, 415–434. [Google Scholar] [CrossRef] [Green Version]
  17. Potts, S.G.; Biesmeijer, J.C.; Kremen, C.; Neumann, P.; Schweiger, O.; Kunin, W.E. Global pollinator declines: Trends, impacts and drivers. Trends Ecol. Evol. 2010, 25, 345–353. [Google Scholar] [CrossRef]
  18. Switanek, M.; Crailsheim, K.; Truhetz, H.; Brodschneider, R. Modelling seasonal effects of temperature and precipitation on honey bee winter mortality in a temperate climate. Sci. Total Environ. 2017, 579, 1581–1587. [Google Scholar] [CrossRef]
  19. VanEngelsdorp, D.; Meixner, M.D. A historical review of managed honey bee populations in Europe and the United States and the factors that may affect them. J. Invertebr. Pathol. 2010, 103, S80–S95. [Google Scholar] [CrossRef]
  20. Neov, B.; Georgieva, A.; Shumkova, R.; Radoslavov, G.; Hristov, P. Biotic and Abiotic Factors Associated with Colonies Mortalities of Managed Honey Bee (Apis mellifera). Diversity 2019, 11, 237. [Google Scholar] [CrossRef] [Green Version]
  21. Yamamoto, I.a.C.; John, E. Nicotinoid Insecticides and the Nicotinic Acetylcholine Receptor; Springer: Tokyo, Japan, 1999. [Google Scholar]
  22. Elbert, A.; Haas, M.; Springer, B.; Thielert, W.; Nauen, R. Applied aspects of neonicotinoid uses in crop protection. Pest Manag. Sci. 2008, 64, 1099–1105. [Google Scholar] [CrossRef]
  23. Jeschke, P. Status and outlook for acaricide and insecticide discovery. Pest Manag. Sci. 2021, 77, 64–76. [Google Scholar] [CrossRef] [PubMed]
  24. Jeschke, P.; Nauen, R. Neonicotinoids—From zero to hero in insecticide chemistry. Pest Manag. Sci. 2008, 64, 1084–1098. [Google Scholar] [CrossRef] [PubMed]
  25. Jeschke, P.; Nauen, R.; Beck, M.E. Nicotinic Acetylcholine Receptor Agonists: A Milestone for Modern Crop Protection. Angew. Chem. Int. Ed. 2013, 52, 9464–9485. [Google Scholar] [CrossRef] [PubMed]
  26. Matsuda, K.; Buckingham, S.D.; Kleier, D.; Rauh, J.J.; Grauso, M.; Sattelle, D.B. Neonicotinoids: Insecticides acting on insect nicotinic acetylcholine receptors. Trends Pharmacol. Sci. 2001, 22, 573–580. [Google Scholar] [CrossRef]
  27. Matsuda, K.; Kanaoka, S.; Akamatsu, M.; Sattelle, D.B. Diverse Actions and Target-Site Selectivity of Neonicotinoids: Structural Insights. Mol. Pharmacol. 2009, 76, 1–10. [Google Scholar] [CrossRef] [Green Version]
  28. Tomizawa, M.; Casida, J.E. Neonicotinoid Insecticide Toxicology: Mechanisms of Selective Action. Annu. Rev. Pharmacol. Toxicol. 2005, 45, 247–268. [Google Scholar] [CrossRef] [Green Version]
  29. Tomizawa, M.; Lee, D.L.; Casida, J.E. Neonicotinoid Insecticides:  Molecular Features Conferring Selectivity for Insect versus Mammalian Nicotinic Receptors. J. Agric. Food Chem. 2000, 48, 6016–6024. [Google Scholar] [CrossRef]
  30. Bonmatin, J.M.; Giorio, C.; Girolami, V.; Goulson, D.; Kreutzweiser, D.P.; Krupke, C.; Liess, M.; Long, E.; Marzaro, M.; Mitchell, E.A.D.; et al. Environmental fate and exposure; Neonicotinoids and fipronil. Environ. Sci. Pollut. Res. 2015, 22, 35–67. [Google Scholar] [CrossRef]
  31. Borsuah, J.F.; Messer, T.L.; Snow, D.D.; Comfort, S.D.; Mittelstet, A.R. Literature Review: Global Neonicotinoid Insecticide Occurrence in Aquatic Environments. Water 2020, 12, 3388. [Google Scholar] [CrossRef]
  32. Wood, T.J.; Goulson, D. The environmental risks of neonicotinoid pesticides: A review of the evidence post 2013. Environ. Sci. Pollut. Res. 2017, 24, 17285–17325. [Google Scholar] [CrossRef] [PubMed]
  33. Van der Sluijs, J.P.; Simon-Delso, N.; Goulson, D.; Maxim, L.; Bonmatin, J.-M.; Belzunces, L.P. Neonicotinoids, bee disorders and the sustainability of pollinator services. Curr. Opin. Environ. Sustain. 2013, 5, 293–305. [Google Scholar] [CrossRef]
  34. Blacquière, T.; Smagghe, G.; van Gestel, C.A.M.; Mommaerts, V. Neonicotinoids in bees: A review on concentrations, side-effects and risk assessment. Ecotoxicology 2012, 21, 973–992. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Decourtye, A.; Devillers, J. Ecotoxicity of neonicotinoid insecticides to bees. Adv. Exp. Med. Biol. 2010, 683, 85–95. [Google Scholar] [CrossRef] [PubMed]
  36. Di Prisco, G.; Cavaliere, V.; Annoscia, D.; Varricchio, P.; Caprio, E.; Nazzi, F.; Gargiulo, G.; Pennacchio, F. Neonicotinoid clothianidin adversely affects insect immunity and promotes replication of a viral pathogen in honey bees. Proc. Natl. Acad. Sci. USA 2013, 110, 18466–18471. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Goulson, D. Review: An overview of the environmental risks posed by neonicotinoid insecticides. J. Appl. Ecol. 2013, 50, 977–987. [Google Scholar] [CrossRef]
  38. Hatjina, F.; Papaefthimiou, C.; Charistos, L.; Dogaroglu, T.; Bouga, M.; Emmanouil, C.; Arnold, G. Sublethal doses of imidacloprid decreased size of hypopharyngeal glands and respiratory rhythm of honeybees in vivo. Apidologie 2013, 44, 467–480. [Google Scholar] [CrossRef] [Green Version]
  39. Henry, M.; Béguin, M.; Requier, F.; Rollin, O.; Odoux, J.-F.; Aupinel, P.; Aptel, J.; Tchamitchian, S.; Decourtye, A. A Common Pesticide Decreases Foraging Success and Survival in Honey Bees. Science 2012, 336, 348–350. [Google Scholar] [CrossRef]
  40. Lundin, O.; Rundlöf, M.; Smith, H.G.; Fries, I.; Bommarco, R. Neonicotinoid Insecticides and Their Impacts on Bees: A Systematic Review of Research Approaches and Identification of Knowledge Gaps. PLoS ONE 2015, 10, e0136928. [Google Scholar] [CrossRef] [Green Version]
  41. Schneider, C.W.; Tautz, J.; Grünewald, B.; Fuchs, S. RFID Tracking of Sublethal Effects of Two Neonicotinoid Insecticides on the Foraging Behavior of Apis mellifera. PLoS ONE 2012, 7, e30023. [Google Scholar] [CrossRef] [Green Version]
  42. Tackenberg, M.C.; Giannoni-Guzmán, M.A.; Sanchez-Perez, E.; Doll, C.A.; Agosto-Rivera, J.L.; Broadie, K.; Moore, D.; McMahon, D.G. Neonicotinoids disrupt circadian rhythms and sleep in honey bees. Sci. Rep. 2020, 10, 17929. [Google Scholar] [CrossRef] [PubMed]
  43. Tavares, D.A.; Dussaubat, C.; Kretzschmar, A.; Carvalho, S.M.; Silva-Zacarin, E.C.M.; Malaspina, O.; Bérail, G.; Brunet, J.-L.; Belzunces, L.P. Exposure of larvae to thiamethoxam affects the survival and physiology of the honey bee at post-embryonic stages. Environ. Pollut. 2017, 229, 386–393. [Google Scholar] [CrossRef] [PubMed]
  44. Tavares, D.A.; Roat, T.C.; Carvalho, S.M.; Silva-Zacarin, E.C.M.; Malaspina, O. In vitro effects of thiamethoxam on larvae of Africanized honey bee Apis mellifera (Hymenoptera: Apidae). Chemosphere 2015, 135, 370–378. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Tosi, S.; Nieh, J.C. A common neonicotinoid pesticide, thiamethoxam, alters honey bee activity, motor functions, and movement to light. Sci. Rep. 2017, 7, 15132. [Google Scholar] [CrossRef]
  46. Whitehorn, P.R.; O’Connor, S.; Wackers, F.L.; Goulson, D. Neonicotinoid Pesticide Reduces Bumble Bee Colony Growth and Queen Production. Science 2012, 336, 351–352. [Google Scholar] [CrossRef] [Green Version]
  47. Annoscia, D.; Di Prisco, G.; Becchimanzi, A.; Caprio, E.; Frizzera, D.; Linguadoca, A.; Nazzi, F.; Pennacchio, F. Neonicotinoid Clothianidin reduces honey bee immune response and contributes to Varroa mite proliferation. Nat. Commun. 2020, 11, 5887. [Google Scholar] [CrossRef]
  48. Catae, A.F.; Roat, T.C.; De Oliveira, R.A.; Nocelli, R.C.; Malaspina, O. Cytotoxic effects of thiamethoxam in the midgut and malpighian tubules of Africanized Apis mellifera (Hymenoptera: Apidae). Microsc. Res. Tech. 2014, 77, 274–281. [Google Scholar] [CrossRef]
  49. Gregorc, A.; Silva-Zacarin, E.C.M.; Carvalho, S.M.; Kramberger, D.; Teixeira, E.W.; Malaspina, O. Effects of Nosema ceranae and thiametoxam in Apis mellifera: A comparative study in Africanized and Carniolan honey bees. Chemosphere 2016, 147, 328–336. [Google Scholar] [CrossRef] [Green Version]
  50. Oliveira, R.A.; Roat, T.C.; Carvalho, S.M.; Malaspina, O. Side-effects of thiamethoxam on the brain andmidgut of the africanized honeybee Apis mellifera (Hymenopptera: Apidae). Environ. Toxicol. 2014, 29, 1122–1133. [Google Scholar] [CrossRef]
  51. Shi, T.-F.; Wang, Y.-F.; Liu, F.; Qi, L.; Yu, L.-S. Sublethal Effects of the Neonicotinoid Insecticide Thiamethoxam on the Transcriptome of the Honey Bees (Hymenoptera: Apidae). J. Econ. Entomol. 2017, 110, 2283–2289. [Google Scholar] [CrossRef]
  52. Paris, L.; Peghaire, E.; Moné, A.; Diogon, M.; Debroas, D.; Delbac, F.; El Alaoui, H. Honeybee gut microbiota dysbiosis in pesticide/parasite co-exposures is mainly induced by Nosema ceranae. J. Invertebr. Pathol. 2020, 172, 107348. [Google Scholar] [CrossRef] [PubMed]
  53. Higes, M.; Martín-Hernández, R.; Meana, A. Nosema ceranae in Europe: An emergent type C nosemosis. Apidologie 2010, 41, 375–392. [Google Scholar] [CrossRef] [Green Version]
  54. Martín-Hernández, R.; Bartolomé, C.; Chejanovsky, N.; Le Conte, Y.; Dalmon, A.; Dussaubat, C.; García-Palencia, P.; Meana, A.; Pinto, M.A.; Soroker, V.; et al. Nosema ceranae in Apis mellifera: A 12 years postdetection perspective. Environ. Microbiol. 2018, 20, 1302–1329. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. García-Palencia, P.; Martín-Hernández, R.; González-Porto, A.-V.; Marin, P.; Meana, A.; Higes, M. Natural infection by Nosema ceranae causes similar lesions as in experimentally infected caged-worker honey bees (Apis mellifera). J. Apic. Res. 2010, 49, 278–283. [Google Scholar] [CrossRef]
  56. Higes, M.; García-Palencia, P.; Martín-Hernández, R.; Meana, A. Experimental infection of Apis mellifera honeybees with Nosema ceranae (Microsporidia). J. Invertebr. Pathol. 2007, 94, 211–217. [Google Scholar] [CrossRef]
  57. Goblirsch, M. Nosema ceranae disease of the honey bee (Apis mellifera). Apidologie 2018, 49, 131–150. [Google Scholar] [CrossRef] [Green Version]
  58. Paris, L.; El Alaoui, H.; Delbac, F.; Diogon, M. Effects of the gut parasite Nosema ceranae on honey bee physiology and behavior. Curr. Opin. Insect Sci. 2018, 26, 149–154. [Google Scholar] [CrossRef]
  59. 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. How natural infection by Nosema ceranae causes honeybee colony collapse. Environ. Microbiol. 2008, 10, 2659–2669. [Google Scholar] [CrossRef]
  60. Dussaubat, C.; Maisonnasse, A.; Alaux, C.; Tchamitchan, S.; Brunet, J.-L.; Plettner, E.; Belzunces, L.P.; Le Conte, Y. Nosema spp. Infection Alters Pheromone Production in Honey Bees (Apis mellifera). J. Chem. Ecol. 2010, 36, 522–525. [Google Scholar] [CrossRef]
  61. Goblirsch, M.; Huang, Z.Y.; Spivak, M. Physiological and Behavioral Changes in Honey Bees (Apis mellifera) Induced by Nosema ceranae Infection. PLoS ONE 2013, 8, e58165. [Google Scholar] [CrossRef]
  62. Ares, A.M.; Nozal, M.J.; Bernal, J.L.; Martín-Hernández, R.; Higes, M.; Bernal, J. Liquid chromatography coupled to ion trap-tandem mass spectrometry to evaluate juvenile hormone III levels in bee hemolymph from Nosema spp. infected colonies. J. Chromatogr. B 2012, 899, 146–153. [Google Scholar] [CrossRef] [PubMed]
  63. Vilarem, C.; Piou, V.; Vogelweith, F.; Vétillard, A. Varroa destructor from the Laboratory to the Field: Control, Biocontrol and IPM Perspectives—A Review. Insects 2021, 12, 800. [Google Scholar] [CrossRef] [PubMed]
  64. Traynor, K.S.; Mondet, F.; de Miranda, J.R.; Techer, M.; Kowallik, V.; Oddie, M.A.Y.; Chantawannakul, P.; McAfee, A. Varroa destructor: A Complex Parasite, Crippling Honey Bees Worldwide. Trends Parasitol. 2020, 36, 592–606. [Google Scholar] [CrossRef] [PubMed]
  65. Rosenkranz, P.; Aumeier, P.; Ziegelmann, B. Biology and control of Varroa destructor. J. Invertebr. Pathol. 2010, 103, S96–S119. [Google Scholar] [CrossRef]
  66. Ramsey, S.D.; Ochoa, R.; Bauchan, G.; Gulbronson, C.; Mowery, J.D.; Cohen, A.; Lim, D.; Joklik, J.; Cicero, J.M.; Ellis, J.D.; et al. Varroa destructor feeds primarily on honey bee fat body tissue and not hemolymph. Proc. Natl. Acad. Sci. USA 2019, 116, 1792–1801. [Google Scholar] [CrossRef] [Green Version]
  67. Weinberg, K.P.; Madel, G. The influence of the mite Varroa jacobsoni Oud. on the protein concentration and the haemolymph volume of the brood of worker bees and drones of the honey bee Apis mellifera L. Apidologie 1985, 16, 421–436. [Google Scholar] [CrossRef] [Green Version]
  68. Sammataro, D.; Gerson, U.; Needham, G. Parasitic Mites of Honey Bees: Life History, Implications, and Impact. Annu. Rev. Entomol. 2000, 45, 519–548. [Google Scholar] [CrossRef]
  69. Yang, X.; Cox-Foster, D.L. Impact of an ectoparasite on the immunity and pathology of an invertebrate: Evidence for host immunosuppression and viral amplification. Proc. Natl. Acad. Sci. USA 2005, 102, 7470–7475. [Google Scholar] [CrossRef] [Green Version]
  70. Kralj, J.; Brockmann, A.; Fuchs, S.; Tautz, J. The parasitic mite Varroa destructor affects non-associative learning in honey bee foragers, Apis mellifera L. J. Comp. Physiol. A 2007, 193, 363–370. [Google Scholar] [CrossRef]
  71. Kralj, J.; Fuchs, S. Parasitic Varroa destructor mites influence flight duration and homing ability of infested Apis mellifera foragers. Apidologie 2006, 37, 577–587. [Google Scholar] [CrossRef] [Green Version]
  72. Coulon, M.; Dalmon, A.; Di Prisco, G.; Prado, A.; Arban, F.; Dubois, E.; Ribière-Chabert, M.; Alaux, C.; Thiéry, R.; Le Conte, Y. Interactions Between Thiamethoxam and Deformed Wing Virus Can Drastically Impair Flight Behavior of Honey Bees. Front. Microbiol. 2020, 11, 722. [Google Scholar] [CrossRef] [PubMed]
  73. Alaux, C.; Brunet, J.-L.; Dussaubat, C.; Mondet, F.; Tchamitchan, S.; Cousin, M.; Brillard, J.; Baldy, A.; Belzunces, L.P.; Le Conte, Y. Interactions between Nosema microspores and a neonicotinoid weaken honeybees (Apis mellifera). Environ. Microbiol. 2010, 12, 774–782. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Pettis, J.S.; vanEngelsdorp, D.; Johnson, J.; Dively, G. Pesticide exposure in honey bees results in increased levels of the gut pathogen Nosema. Naturwissenschaften 2012, 99, 153–158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Doublet, V.; Labarussias, M.; de Miranda, J.R.; Moritz, R.F.A.; Paxton, R.J. Bees under stress: Sublethal doses of a neonicotinoid pesticide and pathogens interact to elevate honey bee mortality across the life cycle. Environ. Microbiol. 2015, 17, 969–983. [Google Scholar] [CrossRef]
  76. Vidau, C.; Diogon, M.; Aufauvre, J.; Fontbonne, R.; Viguès, B.; Brunet, J.-L.; Texier, C.; Biron, D.G.; Blot, N.; El Alaoui, H.; et al. Exposure to Sublethal Doses of Fipronil and Thiacloprid Highly Increases Mortality of Honeybees Previously Infected by Nosema ceranae. PLoS ONE 2011, 6, e21550. [Google Scholar] [CrossRef] [Green Version]
  77. Tesovnik, T.; Zorc, M.; Ristanić, M.; Glavinić, U.; Stevanović, J.; Narat, M.; Stanimirović, Z. Exposure of honey bee larvae to thiamethoxam and its interaction with Nosema ceranae infection in adult honey bees. Environ. Pollut. 2020, 256, 113443. [Google Scholar] [CrossRef]
  78. Dussaubat, C.; Brunet, J.-L.; Higes, M.; Colbourne, J.K.; Lopez, J.; Choi, J.-H.; Martín-Hernández, R.; Botías, C.; Cousin, M.; McDonnell, C.; et al. Gut pathology and responses to the microsporidium Nosema ceranae in the honey bee Apis mellifera. PLoS ONE 2012, 7, e37017. [Google Scholar] [CrossRef] [Green Version]
  79. Morfin, N.; Goodwin, P.H.; Guzman-Novoa, E. The combined effects of Varroa destructor parasitism and exposure to neonicotinoids affects honey bee (Apis mellifera L.) Memory and Gene Expression. Biology 2020, 9, 237. [Google Scholar] [CrossRef]
  80. Odemer, R. Effects of Chronic Pesticide and Pathogen Exposure on Honey Bee (Apis mellifera L.) Health at the Colony Level; Kommunikations, Informations und Medienzentrum der Universität Hohenheim: Stuttgart, Germany, 2018. [Google Scholar]
  81. Van Dooremalen, C.; Cornelissen, B.; Poleij-Hok-Ahin, C.; Blacquière, T. Single and interactive effects of Varroa destructor, Nosema spp., and imidacloprid on honey bee colonies (Apis mellifera). Ecosphere 2018, 9, e02378. [Google Scholar] [CrossRef] [Green Version]
  82. European and Mediterranean Plant Protection Organization. PP 1/170 (4): Side-effects on honeybees. EPPO Bull. 2010, 40, 313–319. [Google Scholar] [CrossRef]
  83. Schulz, D.J.; Huang, Z.-Y.; Robinson, G.E. Effects of colony food shortage on behavioral development in honey bees. Behav. Ecol. Sociobiol. 1998, 42, 295–303. [Google Scholar] [CrossRef]
  84. Higes, M.; Nozal, M.J.; Alvaro, A.; Barrios, L.; Meana, A.; Martín-Hernández, R.; Bernal, J.L.; Bernal, J. The stability and effectiveness of fumagillin in controlling Nosema ceranae (Microsporidia) infection in honey bees (Apis mellifera) under laboratory and field conditions. Apidologie 2011, 42, 364–377. [Google Scholar] [CrossRef] [Green Version]
  85. Ministerio de Agricultura, Pesca y Alimentación. Real Decreto 608/2006, de 19 de Mayo, por el que se Establece y Regula un Programa Nacional de Lucha y Control de las Enfermedades de las Abejas de la Miel; Ministerio de Agricultura: Madrid, Spain, 2006. [Google Scholar]
  86. International Soil Reference and Information Centre. Procedures for Soil Analysis, 6th ed.; van Reeuwijk, L., Ed.; ISRIC: Wageningen, The Netherlands, 2002. [Google Scholar]
  87. Bouyoucos, G.J. A Recalibration of the Hydrometer Method for Making Mechanical Analysis of Soils1. Agron. J. 1951, 43, 434–438. [Google Scholar] [CrossRef] [Green Version]
  88. Walter, H.; Lieth, H. Klimadiagramm—Weltatlas; VEB Gustav Fisher: Jena, Germany, 1967. [Google Scholar]
  89. Rivas-Martinez, S. Worldwide Bioclimatic Classification System, 1996–2021. Available online: http://www.globalbioclimatics.org/default.htm (accessed on 7 July 2021).
  90. Pérez, J.L.; Higes, M.; Suárez, M.; Llorente, J.; Meana, A. Easy ways to determine honey bee mortality using dead-bee traps. J. Apic. Res. 2001, 40, 25–28. [Google Scholar] [CrossRef]
  91. Botías, C.; Martín-Hernández, R.; Días, J.; García-Palencia, P.; Matabuena, M.; Juarranz, Á.; Barrios, L.; Meana, A.; Nanetti, A.; Higes, M. The effect of induced queen replacement on Nosema spp. infection in honey bee (Apis mellifera iberiensis) colonies. Environ. Microbiol. 2012, 14, 845–859. [Google Scholar] [CrossRef] [PubMed]
  92. Bernal, J.; Garrido-Bailón, E.; Del Nozal, M.J.; González-Porto, A.V.; Martín-Hernández, R.; Diego, J.C.; Jiménez, J.J.; Bernal, J.L.; Higes, M. Overview of pesticide residues in stored pollen and their potential effect on bee colony (Apis mellifera) losses in Spain. J. Econ. Entomol. 2010, 103, 1964–1971. [Google Scholar] [CrossRef] [PubMed]
  93. Cepero, A.; Ravoet, J.; Gómez-Moracho, T.; Bernal, J.L.; Del Nozal, M.J.; Bartolomé, C.; Maside, X.; Meana, A.; González-Porto, A.V.; de Graaf, D.C.; et al. Holistic screening of collapsing honey bee colonies in Spain: A case study. BMC Res. Notes 2014, 7, 649. [Google Scholar] [CrossRef] [Green Version]
  94. Higes, M.; Martín-Hernández, R.; Garrido-Bailón, E.; González-Porto, A.V.; García-Palencia, P.; Meana, A.; Del Nozal, M.J.; Mayo, R.; Bernal, J.L. Honeybee colony collapse due to Nosema ceranae in professional apiaries. Environ. Microbiol. Rep. 2009, 1, 110–113. [Google Scholar] [CrossRef]
  95. Hidalgo, M.I.B.; Pacheco, M.L.; Pacheco, J. Origen floral de las cargas de polen recogidas por Apis mellifera L. en Álora (Málaga, España). Acta Bot. Malacit. 1990, 15, 33–44. [Google Scholar] [CrossRef]
  96. Louveaux, J.; Maurizio, A.; Vorwohl, G. Methods of Melissopalynology. Bee World 1978, 59, 139–157. [Google Scholar] [CrossRef]
  97. Erdtman, G. Morphology—Taxonomy—Ecology. An Introduction to the Study of Pollen Grains and Spores. In Handbook of Palynology; Munksgaard: Copenhagen, Denmark, 1969. [Google Scholar]
  98. Ohe, W.V.D.; Oddo, L.P.; Piana, M.L.; Morlot, M.; Martin, P. Harmonized methods of melissopalynology. Apidologie 2004, 35, S18–S25. [Google Scholar]
  99. Faegri, K.; Iversen, J. Textbook of Pollen Analysis, 4th ed.; John Wiley & Sons: Chichester, UK, 1989. [Google Scholar]
  100. Valdés, B.; Fernández, I.D.M.J. Atlas Polínico de Andalucía Occidental; Instituto de Desarrollo Regional de la Universidad de Sevilla y Excma, Diputación de Cádiz: Sevilla, Spain, 1987. [Google Scholar]
  101. Valverde, S.; Ares, A.M.; Bernal, J.L.; Nozal, M.J.; Bernal, J. Effect of the Storage Conditions (Light and Temperature) on the Detection of Thiamethoxam and Clothianidin Content in Rapeseeds by LC-DAD. Food Anal. Methods 2018, 11, 161–169. [Google Scholar] [CrossRef]
  102. Valverde, S.; Bernal, J.L.; Martín, M.T.; Nozal, M.J.; Bernal, J. Fast determination of neonicotinoid insecticides in bee pollen using QuEChERS and ultra-high performance liquid chromatography coupled to quadrupole time-of-flight mass spectrometry. Electrophoresis 2016, 37, 2470–2477. [Google Scholar] [CrossRef] [PubMed]
  103. Valverde, S.; Ibáñez, M.; Bernal, J.L.; Nozal, M.J.; Hernández, F.; Bernal, J. Development and validation of ultra high performance-liquid chromatography–tandem mass spectrometry based methods for the determination of neonicotinoid insecticides in honey. Food Chem. 2018, 266, 215–222. [Google Scholar] [CrossRef]
  104. Jara, L.; Ruiz, C.; Martín-Hernández, R.; Muñoz, I.; Higes, M.; Serrano, J.; De la Rúa, P. The Effect of Migratory Beekeeping on the Infestation Rate of Parasites in Honey Bee (Apis mellifera) Colonies and on Their Genetic Variability. Microorganisms 2021, 9, 22. [Google Scholar] [CrossRef]
  105. Martín-Hernández, R.; Botías, C.; Bailón, E.G.; Martínez-Salvador, A.; Prieto, L.; Meana, A.; Higes, M. Microsporidia infecting Apis mellifera: Coexistence or competition. Is Nosema ceranae replacing Nosema apis? Environ. Microbiol. 2012, 14, 2127–2138. [Google Scholar] [CrossRef]
  106. Chantawannakul, P.; Ward, L.; Boonham, N.; Brown, M. A scientific note on the detection of honeybee viruses using real-time PCR (TaqMan) in Varroa mites collected from a Thai honeybee (Apis mellifera) apiary. J. Invertebr. Pathol. 2006, 91, 69–73. [Google Scholar] [CrossRef]
  107. Francis, R.; Kryger, P. Single Assay Detection of Acute Bee Paralysis Virus, Kashmir Bee Virus and Israeli Acute Paralysis Virus. J. Apic. Sci. 2012, 56, 137–146. [Google Scholar] [CrossRef] [Green Version]
  108. Medrzycki, P.; Giffard, H.; Aupinel, P.; Belzunces, L.P.; Chauzat, M.-P.; Claßen, C.; Colin, M.E.; Dupont, T.; Girolami, V.; Johnson, R.; et al. Standard methods for toxicology research in Apis mellifera. J. Apic. Res. 2013, 52, 1–60. [Google Scholar] [CrossRef] [Green Version]
  109. Pirk, C.W.W.; de Miranda, J.R.; Kramer, M.; Murray, T.E.; Nazzi, F.; Shutler, D.; van der Steen, J.J.M.; van Dooremalen, C. Statistical guidelines for Apis mellifera research. J. Apic. Res. 2013, 52, 1–24. [Google Scholar] [CrossRef] [Green Version]
  110. Divine, G.W.; Norton, H.J.; Barón, A.E.; Juarez-Colunga, E. The Wilcoxon–Mann–Whitney Procedure Fails as a Test of Medians. Am. Stat. 2018, 72, 278–286. [Google Scholar] [CrossRef] [Green Version]
  111. Fritz, C.O.; Morris, P.E.; Richler, J.J. Effect size estimates: Current use, calculations, and interpretation. J. Exp. Psychol. Gen. 2012, 141, 2–18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  112. Zaiontz, C. Real Statistics Using Excel. Available online: www.real-statistics.com (accessed on 15 March 2021).
  113. Lenhard, W.; Lenhard, A. Calculation of Effect Sizes. Available online: https://www.psychometrica.de/effect_size.html (accessed on 6 June 2021).
  114. Faul, F.; Erdfelder, E.; Lang, A.-G.; Buchner, A. G*Power 3: A flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav. Res. Methods 2007, 39, 175–191. [Google Scholar] [CrossRef] [PubMed]
  115. Friedman, M. The Use of Ranks to Avoid the Assumption of Normality Implicit in the Analysis of Variance. J. Am. Stat. Assoc. 1937, 32, 675–701. [Google Scholar] [CrossRef]
  116. Instituto Geológico y Minero de España. IGME_Geode_50 (MapServer): Continuous Digital Geological Map of Spain, Scale 1:50.000; Instituto Geológico y Minero de España: Madrid, Spain, 2014. [Google Scholar]
  117. World Reference Base, I.W.G. World Reference Base for Soil Resources 2014, Update 2015 International Soil Classification System for Naming Soils and Creating Legends for Soil Maps; FAO: Rome, Italy, 2015; Volume 106. [Google Scholar]
  118. Rivas-Martinez, S. Memoria del Mapa de Series de Vegetación de España; Serie Técnica; Ministerio de Agricultura, Pesca y Alimentación; I.C.O.N.A: Madrid, Spain, 1987. [Google Scholar]
  119. Nedić, N.; Mačukanović-Jocić, M.; Rančić, D.; Rørslett, B.; Šoštarić, I.; Stevanović, Z.D.; Mladenović, M. Melliferous potential of Brassica napus L. subsp. napus (Cruciferae). Arthropod-Plant Interact. 2013, 7, 323–333. [Google Scholar] [CrossRef]
  120. Talavera, S.; Gibbs, P.E.; Herrera, J. Reproductive biology of Cistus ladanifer (Cistaceae). Plant Syst. Evol. 1993, 186, 123–134. [Google Scholar] [CrossRef]
  121. Ortiz, P.L. The Cistaceae as food resources for honey bees in SW Spain. J. Apic. Res. 1994, 33, 136–144. [Google Scholar] [CrossRef]
  122. Atsalakis, E.; Chinou, I.; Makropoulou, M.; Karabournioti, S.; Graikou, K. Evaluation of Phenolic Compounds in Cistus creticus Bee Pollen from Greece. Antioxidant and Antimicrobial Properties. Nat. Prod. Commun. 2017, 12, 1934578X1701201141. [Google Scholar] [CrossRef] [Green Version]
  123. Papaefthimiou, D.; Papanikolaou, A.; Falara, V.; Givanoudi, S.; Kostas, S.; Kanellis, A.K. Genus Cistus: A model for exploring labdane-type diterpenes’ biosynthesis and a natural source of high value products with biological, aromatic, and pharmacological properties. Front. Chem. 2014, 2, 35. [Google Scholar] [CrossRef] [Green Version]
  124. Sponsler, D.B.; Johnson, R.M. Mechanistic modeling of pesticide exposure: The missing keystone of honey bee toxicology. Environ. Toxicol. Chem. 2017, 36, 871–881. [Google Scholar] [CrossRef]
  125. Huang, Z.-Y.; Robinson, G.E. Regulation of honey bee division of labor by colony age demography. Behav. Ecol. Sociobiol. 1996, 39, 147–158. [Google Scholar] [CrossRef]
  126. Schott, M.; Sandmann, M.; Cresswell, J.E.; Becher, M.A.; Eichner, G.; Brandt, D.T.; Halitschke, R.; Krueger, S.; Morlock, G.; Düring, R.-A.; et al. Honeybee colonies compensate for pesticide-induced effects on royal jelly composition and brood survival with increased brood production. Sci. Rep. 2021, 11, 62. [Google Scholar] [CrossRef]
  127. Calabrese, E.J.; Baldwin, L.A. The Hormetic Dose-Response Model Is More Common than the Threshold Model in Toxicology. Toxicol. Sci. 2003, 71, 246–250. [Google Scholar] [CrossRef]
  128. Luckey, T.D. Insecticide Hormoligosis. J. Econ. Entomol. 1968, 61, 7–12. [Google Scholar] [CrossRef] [PubMed]
  129. Southam, C.M.; Ehrlich, J. Effects of Extract of western red-cedar heartwood on certain wood-decaying fungi in culture. Phytopathology 1943, 33, 517–524. [Google Scholar]
  130. Kendig, E.L.; Le, H.H.; Belcher, S.M. Defining hormesis: Evaluation of a complex concentration response phenomenon. Int. J. Toxicol. 2010, 29, 235–246. [Google Scholar] [CrossRef] [PubMed]
  131. Cohen, E. Pesticide-mediated homeostatic modulation in arthropods. Pestic. Biochem. Physiol. 2006, 85, 21–27. [Google Scholar] [CrossRef]
  132. Guedes, R.N.C.; Cutler, G.C. Insecticide-induced hormesis and arthropod pest management. Pest Manag. Sci. 2014, 70, 690–697. [Google Scholar] [CrossRef] [Green Version]
  133. Calabrese, E.J.; Blain, R. The occurrence of hormetic dose responses in the toxicological literature, the hormesis database: An overview. Toxicol. Appl. Pharm. 2005, 202, 289–301. [Google Scholar] [CrossRef]
  134. Costantini, D.; Metcalfe, N.B.; Monaghan, P. Ecological processes in a hormetic framework. Ecol. Lett. 2010, 13, 1435–1447. [Google Scholar] [CrossRef]
  135. Berry, R.; López-Martínez, G. A dose of experimental hormesis: When mild stress protects and improves animal performance. Comp. Biochem. Physiol. Part A Mol. Integr. Physiol. 2020, 242, 110658. [Google Scholar] [CrossRef] [PubMed]
  136. Bryden, J.; Gill, R.J.; Mitton, R.A.A.; Raine, N.E.; Jansen, V.A.A. Chronic sublethal stress causes bee colony failure. Ecol. Lett. 2013, 16, 1463–1469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  137. Perry, C.J.; Søvik, E.; Myerscough, M.R.; Barron, A.B. Rapid behavioral maturation accelerates failure of stressed honey bee colonies. Proc. Natl. Acad. Sci. USA 2015, 112, 3427–3432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  138. Higes, M.; Meana, A.; Bartolomé, C.; Botías, C.; Martín-Hernández, R. Nosema ceranae (Microsporidia), a controversial 21st century honey bee pathogen. Environ. Microbiol. Rep. 2013, 5, 17–29. [Google Scholar] [CrossRef]
  139. Saitanis, C.J.; Agathokleous, E. Stress response and population dynamics: Is Allee effect hormesis? Sci. Total Environ. 2019, 682, 623–628. [Google Scholar] [CrossRef]
  140. Dennis, B.; Kemp, W.P. How Hives Collapse: Allee Effects, Ecological Resilience, and the Honey Bee. PLoS ONE 2016, 11, e0150055. [Google Scholar] [CrossRef] [Green Version]
  141. Courchamp, F.; Clutton-Brock, T.; Grenfell, B. Inverse density dependence and the Allee effect. Trends Ecol. Evol. 1999, 14, 405–410. [Google Scholar] [CrossRef]
  142. Ministerio de Agricultura, Pesca y Alimentación. Informes de Resultados del Programa de Vigilancia Sobre la Pérdida de Colonias de Abejas. Available online: https://www.mapa.gob.es/es/ganaderia/temas/sanidad-animal-higiene-ganadera/sanidad-animal/enfermedades/otras-enfermedades-abejas/otras-enf-abejas.aspx (accessed on 11 November 2021).
  143. Sadd, B.M.; Schmid-Hempel, P. Perspective: Principles of ecological immunology. Evol. Appl. 2009, 2, 113–121. [Google Scholar] [CrossRef]
  144. Liu, Z.; Li, S.; Li, H. Interactions between sublethal doses of thiamethoxam and Nosema ceranae in the honey bee, Apis mellifera. J. Apic. Res. 2021, 60, 717–725. [Google Scholar] [CrossRef]
  145. McKinstry, M.; Chung, C.; Truong, H.; Johnston, B.A.; Snow, J.W. The heat shock response and humoral immune response are mutually antagonistic in honey bees. Sci. Rep. 2017, 7, 8850. [Google Scholar] [CrossRef] [Green Version]
  146. Straub, L.; Williams, G.R.; Vidondo, B.; Khongphinitbunjong, K.; Retschnig, G.; Schneeberger, A.; Chantawannakul, P.; Dietemann, V.; Neumann, P. Neonicotinoids and ectoparasitic mites synergistically impact honeybees. Sci. Rep. 2019, 9, 8159. [Google Scholar] [CrossRef] [PubMed]
  147. Locke, B.; Semberg, E.; Forsgren, E.; de Miranda, J.R. Persistence of subclinical deformed wing virus infections in honeybees following Varroa mite removal and a bee population turnover. PLoS ONE 2017, 12, e0180910. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  148. Möckel, N.; Gisder, S.; Genersch, E. Horizontal transmission of deformed wing virus: Pathological consequences in adult bees (Apis mellifera) depend on the transmission route. J. Gen. Virol. 2011, 92, 370–377. [Google Scholar] [CrossRef] [PubMed]
  149. Potts, R.; Clarke, R.M.; Oldfield, S.E.; Wood, L.K.; de Ibarra, N.H.; Cresswell, J.E. The effect of dietary neonicotinoid pesticides on non-flight thermogenesis in worker bumble bees (Bombus terrestris). J. Insect Physiol. 2018, 104, 33–39. [Google Scholar] [CrossRef] [PubMed]
  150. Tosi, S.; Démares, F.J.; Nicolson, S.W.; Medrzycki, P.; Pirk, C.W.W.; Human, H. Effects of a neonicotinoid pesticide on thermoregulation of African honey bees (Apis mellifera scutellata). J. Insect Physiol. 2016, 93, 56–63. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  151. Bordier, C.; Dechatre, H.; Suchail, S.; Peruzzi, M.; Soubeyrand, S.; Pioz, M.; Pélissier, M.; Crauser, D.; Conte, Y.L.; Alaux, C. Colony adaptive response to simulated heat waves and consequences at the individual level in honeybees (Apis mellifera). Sci. Rep. 2017, 7, 3760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  152. Dalmon, A.; Peruzzi, M.; Le Conte, Y.; Alaux, C.; Pioz, M. Temperature-driven changes in viral loads in the honey bee Apis mellifera. J. Invertebr. Pathol. 2019, 160, 87–94. [Google Scholar] [CrossRef]
  153. McMenamin, A.J.; Daughenbaugh, K.F.; Flenniken, M.L. The Heat Shock Response in the Western Honey Bee (Apis mellifera) is Antiviral. Viruses 2020, 12, 245. [Google Scholar] [CrossRef] [Green Version]
  154. Ryan, S.J.; Carlson, C.J.; Mordecai, E.A.; Johnson, L.R. Global expansion and redistribution of Aedes-borne virus transmission risk with climate change. PLoS Negl. Trop. Dis. 2019, 13, e0007213. [Google Scholar] [CrossRef] [Green Version]
  155. Cutler, G.C.; Rix, R.R. Can poisons stimulate bees? Appreciating the potential of hormesis in bee–pesticide research. Pest Manag. Sci. 2015, 71, 1368–1370. [Google Scholar] [CrossRef]
  156. EFSA Scientific Committee. Statistical Significance and Biological Relevance. EFSA J. 2011, 9, 2372. [Google Scholar] [CrossRef]
  157. Authority, E.F.S.; Ippolito, A.; Focks, A.; Rundlöf, M.; Arce, A.; Marchesi, M.; Neri, F.M.; Rortais, A.; Szentes, C.; Auteri, D. Analysis of background variability of honey bee colony size. EFSA Supporting Publ. 2021, 18, 6518E. [Google Scholar] [CrossRef]
  158. Cutler, G.C.; Scott-Dupree, C.D.; Sultan, M.; McFarlane, A.D.; Brewer, L. A large-scale field study examining effects of exposure to clothianidin seed-treated canola on honey bee colony health, development, and overwintering success. PeerJ 2014, 2, e652. [Google Scholar] [CrossRef] [PubMed]
  159. Flores, J.M.; Gámiz, V.; Gil-Lebrero, S.; Rodríguez, I.; Navas, F.J.; García-Valcárcel, A.I.; Cutillas, V.; Fernández-Alba, A.R.; Hernando, M.D. A three-year large scale study on the risk of honey bee colony exposure to blooming sunflowers grown from seeds treated with thiamethoxam and clothianidin neonicotinoids. Chemosphere 2021, 262, 127735. [Google Scholar] [CrossRef] [PubMed]
  160. Hernando, M.D.; Gámiz, V.; Gil-Lebrero, S.; Rodríguez, I.; García-Valcárcel, A.I.; Cutillas, V.; Fernández-Alba, A.R.; Flores, J.M. Viability of honeybee colonies exposed to sunflowers grown from seeds treated with the neonicotinoids thiamethoxam and clothianidin. Chemosphere 2018, 202, 609–617. [Google Scholar] [CrossRef]
  161. Rolke, D.; Fuchs, S.; Grünewald, B.; Gao, Z.; Blenau, W. Large-scale monitoring of effects of clothianidin-dressed oilseed rape seeds on pollinating insects in Northern Germany: Effects on honey bees (Apis mellifera). Ecotoxicology 2016, 25, 1648–1665. [Google Scholar] [CrossRef] [Green Version]
  162. Woodcock, B.A.; Bullock, J.M.; Shore, R.F.; Heard, M.S.; Pereira, M.G.; Redhead, J.; Ridding, L.; Dean, H.; Sleep, D.; Henrys, P.; et al. Country-specific effects of neonicotinoid pesticides on honey bees and wild bees. Science 2017, 356, 1393–1395. [Google Scholar] [CrossRef] [Green Version]
  163. Osterman, J.; Wintermantel, D.; Locke, B.; Jonsson, O.; Semberg, E.; Onorati, P.; Forsgren, E.; Rosenkranz, P.; Rahbek-Pedersen, T.; Bommarco, R.; et al. Clothianidin seed-treatment has no detectable negative impact on honeybee colonies and their pathogens. Nat. Commun. 2019, 10, 692. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Site of the experimental study: (a) Map of Spain with the Province of Guadalajara highlighted in green; (b) Location of the apiaries during (A: Marchamalo) and after (B: Fuentelahiguera de Albatages) the exposure phase; (c) Detailed location of the apiaries during the exposure period in Marchamalo (C: control, T: treatment); (d) Detailed location of the apiaries after the exposure period in Fuentelahiguera de Albatages (B).
Figure 1. Site of the experimental study: (a) Map of Spain with the Province of Guadalajara highlighted in green; (b) Location of the apiaries during (A: Marchamalo) and after (B: Fuentelahiguera de Albatages) the exposure phase; (c) Detailed location of the apiaries during the exposure period in Marchamalo (C: control, T: treatment); (d) Detailed location of the apiaries after the exposure period in Fuentelahiguera de Albatages (B).
Insects 13 00371 g001
Figure 2. Melissopalynological analysis: (a) The % (w/w) Brassica t. loads in corbicular pollen at the beginning (8 April 2015), middle (13 April 2015), and end (29 April 2015) of winter oilseed rape (wOSR) flowering; (b) The frequency (%) of Brassica t. pollen in beebread Insects 13 00371 i001 and honey Insects 13 00371 i002: C-01 to C-05 = Control hives; C-06 to C-10 = hives exposed to TMX.
Figure 2. Melissopalynological analysis: (a) The % (w/w) Brassica t. loads in corbicular pollen at the beginning (8 April 2015), middle (13 April 2015), and end (29 April 2015) of winter oilseed rape (wOSR) flowering; (b) The frequency (%) of Brassica t. pollen in beebread Insects 13 00371 i001 and honey Insects 13 00371 i002: C-01 to C-05 = Control hives; C-06 to C-10 = hives exposed to TMX.
Insects 13 00371 g002
Figure 3. Colony strength parameters throughout the study in control and treatment (TMX) apiaries: mean accumulated dead bees at the entrance of the colonies (  control; TMX) and mean number of combs covered with bees ( control;TMX) or brood ( control; TMX).
Figure 3. Colony strength parameters throughout the study in control and treatment (TMX) apiaries: mean accumulated dead bees at the entrance of the colonies (  control; TMX) and mean number of combs covered with bees ( control;TMX) or brood ( control; TMX).
Insects 13 00371 g003
Figure 4. Box and Whisker plots of N. ceranae infection in non-forager honey bees (IB; Insects 13 00371 i003 control; Insects 13 00371 i004 TMX) and forager honey bees (FB; Insects 13 00371 i005 control; Insects 13 00371 i006 TMX).
Figure 4. Box and Whisker plots of N. ceranae infection in non-forager honey bees (IB; Insects 13 00371 i003 control; Insects 13 00371 i004 TMX) and forager honey bees (FB; Insects 13 00371 i005 control; Insects 13 00371 i006 TMX).
Insects 13 00371 g004
Table 1. Field schedule.
Table 1. Field schedule.
Sampling DateForager DeathColony StrengthPathogensBee PollenBeebreadHoney
30 March 2015XXX
8 April 2015X X
10 April 2015X
13 April 2015XX X
15 April 2015X
17 April 2015X
20 April 2015X
22 April 2015X
24 April 2015X
27 April 2015X
29 April 2015X X
4 May 2015X
6 May 2015X
8 May 2015X
11 May 2015XXX XX
15 July 2015 XX
15 September 2015 XX
17 September 2015
15 November 2015 XX
15 March 2016 XX
Table 2. Descriptive statistics of the variables of strength in control and treatment (TMX) apiaries at the end of the exposure period.
Table 2. Descriptive statistics of the variables of strength in control and treatment (TMX) apiaries at the end of the exposure period.
VariableDescriptive StatisticsControl (n = 5)TMX (n = 5)
Accumulated deathsMin8092
Max192394
Median129214
Mean127.6214.4
Standard deviation (SD)44.106116.468
Variation coefficient (%)34.56554.323
Mean 95% confidence interval[72.836; 182.364][69.785; 359.015]
Power (1-β)0.379
Cohen’s d0.863
Number of combs with worker beesMin1111
Max1416
Median1216
Mean12.215
Standard deviation (SD)1.0952.236
Mean 95% confidence interval[10.84; 13.56][12.22; 17.77]
Variation coefficient (%)8.97914.907
Power (1-β)0.541
Cohen’s d1.245
Number of combs with broodMin29
Max1313
Median10.511.0
Mean9.011.0
Standard deviation (SD)4.9661.826
Variation coefficient (%)55.18416.598
Mean 95% confidence interval[1.097; 16.903][8.095; 13.905]
Power (1-β)0.103
Cohen’s d0.310
Honey production (kg)Min10.2413.02
Max45.2245.58
Median16.20432.034
Mean24.75531.977
Standard deviation (SD)15.07512.21
Variation coefficient (%)60.89438.184
Mean 95% confidence interval[6.038; 43.473][16.816; 47.137]
Power (1-β)0.165
Cohen’s d0.475
Table 3. Detections of Black queen cell virus (BQCV) and Deformed wing virus (DWV) during the study in non-forager (IB) and forager bees (FB).
Table 3. Detections of Black queen cell virus (BQCV) and Deformed wing virus (DWV) during the study in non-forager (IB) and forager bees (FB).
Hive30 March
2015
11 May
2015
15 July
2015
15 September
2015
17 November
2015
15 March 2016
IBFBIBFBIBFBIBFBIBFBIBFB
C 01BQCV-BQCVBQCVBQCVBQCV--DWVDWVDWVDWV
C 02BQCV-BQCV-BQCVBQCVBQCV-BQCV
DWV
DWV--
C 03BQCVBQCV-DWVBQCVBQCVBQCVBQCV-BQCVBQCV-DWVDWV--
C 04BQCVBQCV-DWVBQCVDWVBQCVBQCVBQCVBQCVBQCV-DWVDWVBQCV-
C 05DWVBQCVBQCV-DWV-DWVBQCVBQCV-BQCV-DWVDWVBQCV-
C 06BQCVBQCVBQCV-BQCVBQCV--BQCV-DWVDWVBQCV-
C 07BQCV-BQCV-BQCVBQCV-BQCVBQCV-DWVDWVBQCV-
C 08BQCV-BQCVBQCVBQCVBQCV--BQCV-DWVDWV-DWV
C 09-BQCV--BQCVBQCVBQCV-BQCV-DWVDWV--
C 10-BQCV-DWV-BQCVBQCVBQCVBQCV-BQCV-DWVDWV--
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alonso-Prados, E.; González-Porto, A.V.; García-Villarubia, C.; López-Pérez, J.A.; Valverde, S.; Bernal, J.; Martín-Hernández, R.; Higes, M. Effects of Thiamethoxam-Dressed Oilseed Rape Seeds and Nosema ceranae on Colonies of Apis mellifera iberiensis, L. under Field Conditions of Central Spain. Is Hormesis Playing a Role? Insects 2022, 13, 371. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13040371

AMA Style

Alonso-Prados E, González-Porto AV, García-Villarubia C, López-Pérez JA, Valverde S, Bernal J, Martín-Hernández R, Higes M. Effects of Thiamethoxam-Dressed Oilseed Rape Seeds and Nosema ceranae on Colonies of Apis mellifera iberiensis, L. under Field Conditions of Central Spain. Is Hormesis Playing a Role? Insects. 2022; 13(4):371. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13040371

Chicago/Turabian Style

Alonso-Prados, Elena, Amelia Virginia González-Porto, Carlos García-Villarubia, José Antonio López-Pérez, Silvia Valverde, José Bernal, Raquel Martín-Hernández, and Mariano Higes. 2022. "Effects of Thiamethoxam-Dressed Oilseed Rape Seeds and Nosema ceranae on Colonies of Apis mellifera iberiensis, L. under Field Conditions of Central Spain. Is Hormesis Playing a Role?" Insects 13, no. 4: 371. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13040371

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop