Next Article in Journal
The First Cretaceous Epyrine Wasp (Hymenoptera: Bethylidae): A New Genus and Species from Early Cenomanian Kachin Amber
Next Article in Special Issue
The Impact of Wildflower Habitat on Insect Functional Group Abundance in Turfgrass Systems
Previous Article in Journal
Evaluation of Cotton Fleahopper (Pseudatomoscelis seriatus (Reuter)) Feeding on Mpp51Aa2-Traited Cotton Utilizing Electrical Penetration Graph (EPG) Waveforms
Previous Article in Special Issue
Six Steps towards a Spatial Design for Large-Scale Pollinator Surveillance Monitoring
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 15
Issue 5
10.3390/insects15050317
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Different Sensitivity of Flower-Visiting Diptera to a Neonicotinoid Insecticide: Expanding the Base for a Multiple-Species Risk Assessment Approach

by
Cátia Ariana Henriques Martins
1,
Celeste Azpiazu
2,3,
Jordi Bosch
2,
Giovanni Burgio
1,
Maria Luisa Dindo
1,
Santolo Francati
1,
Daniele Sommaggio
4,5 and
Fabio Sgolastra
1,*
1
Dipartimento di Scienze e Tecnologie Agro-Alimentari, Alma Mater Studiorum Università di Bologna, 40127 Bologna, Italy
2
CREAF, Centre de Recerca Ecològica i Aplicacions Forestals, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
3
Universidad Politécnica de Madrid, 28040 Madrid, Spain
4
Dipartimento di Scienze della Vita, Università di Modena e Reggio Emilia, 41121 Modena, Italy
5
National Biodiversity Future Center (NBFC), Piazza Marina 61, 90133 Palermo, Italy
*
Author to whom correspondence should be addressed.
Insects 2024, 15(5), 317; https://doi.org/10.3390/insects15050317
Submission received: 26 March 2024 / Revised: 18 April 2024 / Accepted: 24 April 2024 / Published: 29 April 2024
(This article belongs to the Special Issue Pollinator Biodiversity and Ecosystem Services)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Insect pollinators play an essential service in agricultural systems, but are commonly exposed to pesticides. Although pollinators are present in several insect orders, above all dipterans, information on pesticide sensitivity is mostly restricted to bees. We assessed the sensitivity of two hoverflies (Sphaerophoria rueppellii, Eristalinus aeneus) and one tachinid fly (Exorista larvarum) to a neonicotinoid insecticide (Confidor®, imidacloprid). We adapted the standardized methodology of acute contact exposure in honey bees to build dose–response curves and calculate median lethal doses (LD50) for the three species. S. rueppelli was the most sensitive, E. aeneus the least. Results were compared with those available in the literature for other pollinator species using a species sensitivity distribution (SSD) approach: as a result, the 95th percentile of pollinator species would be protected by a safety factor of 100 times the Apis mellifera endpoint. Dipterans were less sensitive to imidacloprid than most bee species. We measured the number of eggs laid following exposure to different insecticide doses and assessed the potential trade-off between oviposition and survival through the sublethal sensitivity index (SSI). Pesticide exposure had a significant effect on fecundity, and SSI values indicated that oviposition is a sensitive endpoint for the three dipteran species tested.

Abstract

Insects play an essential role as pollinators of wild flowers and crops. At the same time, pollinators in agricultural environments are commonly exposed to pesticides, compromising their survival and the provision of pollination services. Although pollinators include a wide range of species from several insect orders, information on pesticide sensitivity is mostly restricted to bees. In addition, the disparity of methodological procedures used for different insect groups hinders the comparison of toxicity data between bees and other pollinators. Dipterans are a highly diverse insect order that includes some important pollinators. Therefore, in this study, we assessed the sensitivity of two hoverflies (Sphaerophoria rueppellii, Eristalinus aeneus) and one tachinid fly (Exorista larvarum) to a neonicotinoid insecticide (Confidor®, imidacloprid) following a comparative approach. We adapted the standardized methodology of acute contact exposure in honey bees to build dose–response curves and calculate median lethal doses (LD50) for the three species. The methodology consisted in applying 1 µL of the test solution on the thorax of each insect. Sphaerophoria rueppelli was the most sensitive species (LD50 = 10.23 ng/insect), and E. aeneus (LD50 = 18,176 ng/insect) the least. We then compared our results with those available in the literature for other pollinator species using species sensitivity distribution (SSD). Based on the SSD curve, the 95th percentile of pollinator species would be protected by a safety factor of 100 times the Apis mellifera endpoint. Overall, dipterans were less sensitive to imidacloprid than most bee species. As opposed to most bee species, oviposition and fecundity of many dipteran species can be reliably assessed in the laboratory. We measured the number of eggs laid following exposure to different insecticide doses and assessed the potential trade-off between oviposition and survival through the sublethal sensitivity index (SSI). Exposure to imidacloprid had a significant effect on fecundity, and SSI values indicated that oviposition is a sensitive endpoint for the three dipteran species tested. Future studies should integrate this information related to population dynamics in simulation models for environmental risk assessment.

Graphical Abstract

1. Introduction

Insect pollinators play an essential role in ecosystem functioning and crop pollination [1,2]. Although bees usually take most of the credit for this ecosystem service, awareness of the importance of other insect pollinators is growing [3]. Dipterans, in particular, are one of the most species-rich insect orders, comprising more than 160,000 species grouped in ca. 160 families, of which at least 75 are flower visitors [4,5,6]. Different groups of flies, including hoverflies (Syrphidae), bee flies (Bombyliidae) and many groups of muscoid flies, visit flowers mostly to feed on pollen and nectar, but also to find prey and hosts for their larvae, oviposit, mate and take shelter [7]. Several studies show that flies may be highly effective pollinators [3,8,9,10] and highlight their contribution to wild flower and crop pollination by visiting at least 555 plant species, including over 100 cultivated plants [3,11,12]. Excluding Hymenoptera, most of which are bees, Diptera is the most frequent flower-visiting order (72% of crops) [13]. Based on their dual role as pollinators [5,14,15] and biological control agents [16,17] some studies advocate for the use of fly populations in integrated pest management [10].
Declines in pollinator abundance and diversity over the last century have been mostly documented for bees [18,19] and to a lesser extent hoverflies [20,21,22]. Importantly, these declines are partly associated with agricultural intensification and its strong reliance on systemic and non-systemic pesticides, including seed treatments [23,24,25]. To safeguard the well-being and survival of pollinators and other non-target organisms, pesticides undergo a process of environmental risk assessment [26,27]. However, non-bee pollinators are poorly represented in current environmental risk assessment. Bee risk assessment schemes rely on a surrogate species, the western honey bee, Apis mellifera L., and are based on the calculation of the dose of pesticide that is lethal to 50% of the population (LD50) following acute exposure [28,29]. This approach does not account for potential sensitivity differences among bee species [30,31,32,33,34,35]. Recent studies have compared the sensitivity of honey bees and other bees, including bumblebees, mason bees and stingless bees [30,32,36,37,38]. Non-bee pollinators are also overlooked in risk assessment procedures to assess the effects of pesticides on non-target arthropods [39]. For these reasons, the European Food Safety Authority recently published a road map to expand environmental risk assessment to additional pollinator species [40]. In addition, the few studies addressing pesticide sensitivity in non-bee pollinator taxa use different methodological procedures, hindering direct comparison with results obtained on bees [41,42,43].
In this study, we assessed the toxicity of a neonicotinoid insecticide (imidacloprid) in two hoverflies, Sphaerophoria rueppellii (Wiedemann) (Diptera: Syrphidae) and Eristalinus aeneus (Scopoli) (Diptera: Syrphidae), and a parasitoid fly, Exorista larvarum (L.) (Diptera: Tachinidae). The three tested species play important roles as pollinators for several plant species and in various ecological contexts. Exorista larvarum is a polyphagous gregarious larval parasitoid of several Lepidoptera, such as the gypsy and cabbage moths, which are pests in forest and agricultural environments, respectively [44]. It is also a common visitor of Apiaceae [45]. Sphaerophoria rueppellii is an aphidophagous hoverfly [46,47]. Populations of this species are commercially available due to their high effectiveness as biological control agents of pest aphids and as a pollinator of some crops [48]. Eristalinus aeneus is a saprophagous hoverfly. The adults are effective pollinators of hybrid seed crops [49].
Due to their foraging activities in crop plants during adult stage, all these species can be routinely exposed to pesticides. Imidacloprid is a neurotoxic insecticide with a high affinity for the nicotinic acetylcholine receptor (nAChR). Due to its physicochemical characteristics, imidacloprid can easily translocate to all plant tissues, thus providing protection from a wide range of herbivorous (mainly sap-feeding) insects. Imidacloprid has been registered in 120 countries for use on over 140 crops against several insect pests, including aphids, whiteflies, and the Colorado potato beetle [50]. It was introduced in the market in 1991 and quickly became the most worldwide used insecticide, accounting for 41.5% of the neonicotinoid market in 2009 [50]. The widespread use of this compound, together with its systemic properties and long persistence in the environment, raised concern about its impacts on biodiversity, ecosystem functioning and services [51,52,53]. Due to its high acute toxicity to bees, its use in the EU was banned by the European Commission in 2018 [54]. Despite this ban, many EU Member States provide emergency authorizations for its use, and it is still widely used in many parts of the world [55].
Historically, pesticide acute toxicity to honey bees has been tested following the protocols developed by the Organisation for Economic Cooperation and Development (OECD). These protocols are currently used for honey bees and, with some minor adaptations, for bumblebees [56,57,58,59]. However, similar standard protocols for non-bee pollinators are lacking. To cover this gap and gain comparable data to assess potential interspecific differences among pollinators, in this study, we adapted the standardized acute contact test for honey bees [56] to expose individual flies to a range of insecticide doses, build dose–response curves, and calculate LD50 values. We then followed a species sensitivity distribution curve (SSD) approach to compare the sensitivity of the three dipterans tested to that of other pollinators, including honey bees. This distribution statistically describes the variation among a set of species to a particular chemical or mixture, and attributes a fifth percentile hazard concentration/dose (HC5/HD5) that ensures a proper level of protection [60]. The SSD approach is a probabilistic approach to establish adequate protection exposure levels for most or all species when limited single-species toxicity data are available [61]. Based on previous studies showing large variability in pesticide sensitivity among different bee species [31], we expected large variability among the three dipteran species and between bees and dipterans.
The population effects of pesticide exposure not only depend on the intrinsic sensitivity (i.e., mortality) of the target species but also on their ability to reproduce during sublethal exposure, which is constrained by life history traits. Under stress conditions, some species prioritize investment in physiological maintenance and survival at the expense of reproduction. On the other hand, other species, often called opportunistic, maintain high levels of investment in reproduction until they reach levels of lethal exposure [62]. The population effects of a toxicant can be measured by the sublethal sensitivity index (SSI), which is the ratio between the LD50 or LC50 (median lethal dose or concentration) and the NOEL (no observed effect level) for reproduction [62]. For a given LD50, population effects will manifest earlier (at lower doses) in species with high SSI, that is, in species prioritizing energetic investment in somatic maintenance, including stress resistance and immunity [63]. Measuring reproductive output in the laboratory is not feasible for most bee species, which require flying space and large amounts of pollen and nectar resources [64]. However, by providing adequate oviposition substrates, fecundity (number of eggs laid) can be reasonably assessed in many dipterans under laboratory conditions. In this study, we assessed the effects of dosage on fecundity (number of eggs laid) and used the results to calculate the SSI for the three species.
The objectives of our study were: (1) to measure the contact toxicity of imidacloprid in three dipteran species using a methodology comparable to that used for bees; (2) to assess interspecific variation in sensitivity to imidacloprid using a sensitivity distribution curve approach; and (3) to assess the sublethal effects of imidacloprid exposure on fecundity and assess the vulnerability of the three species using the SSI approach.

2. Materials and Methods

2.1. Species, Populations, and Test Conditions

A population of E. larvarum was maintained at the Department of Agricultural and Food Sciences (DISTAL—University of Bologna, Italy) following the standard rearing procedures described by Dindo et al. [65]. Galleria mellonella L. (Lepidoptera: Pyralidae) larvae, reared at 30 ± 1 °C, 65 ± 5% relative humidity, and complete darkness, were used as a factitious host [66]. Mature E. larvarum puparia weighing 35–55 mg were selected for the experiment. Newly emerged individuals were kept in Plexiglas cages (40 cm × 30 cm × 30 cm) in a 1:1 sex ratio and left to mate for 3 days. Mated females were individually exposed to imidacloprid. Following exposure, six groups of five females of the same treatment group were transferred to Plexiglas cages (20 cm × 20 cm × 20 cm) within each of the eight treatments.
Sphaerophoria rueppellii were obtained from a population reared in Belgium (Spharophoria-System, Biobest N.V.) and shipped as pupae to DISTAL in February–April 2021. Upon arrival, pupae were introduced to standard mesh cages (60 cm × 40 cm × 40 cm) and transferred to a 26 °C chamber for incubation. Upon emergence, we let adults (1:1 sex ratio) mate for 2 days. Following the procedure previously described and after the exposure to imidacloprid, females were individually transferred to cages (V = 150 cc) with the lid perforated for air circulation and maintained there until the end of the test.
Eristalinus aeneus pupae were obtained from a population reared in Spain (Goldfly®, Polyfly S.L., Spain) and shipped to DISTAL in November 2021 and January 2022. Pupae were incubated as described above. Emerged individuals were kept in standard mesh cages (60 cm × 40 cm × 40 cm) for 5 days in a 1:1 sex ratio. Mated females were treated as described for S. rueppellii.
All flies were kept under light:dark conditions (16:8 h) in a climatic chamber at 26 ± 1 °C and 70–80% relative humidity during emergence and mating and for the duration of the test. Upon emergence, all fly species were fed ad libitum, with honey bee pollen pellets from organic beekeeping (Bona Mel®, Spain and sugar cubes. Pollen was supplied to promote ovary development [7]. Distilled water was provided in plastic 50 mL drinking straws (E. larvarum) and in 1 mL Eppendorf tubes with the tip cut off and plugged with dampened cotton (S. rueppellii and E. aeneus). Water was refilled as needed. Initial sample sizes were ca. 30 females per species and dose. Each female represented a replicate in the statistical analyses. A randomly selected group of mated females of each species (n = 26–55) were weighted to obtain an average measure of fresh body weight.

2.2. Imidacloprid Solutions

We used the commercially available formulation Confidor®, CropScience S.r.l., Italy (imidacloprid 20% w/v). Following preliminary range-finding tests, we started by diluting Confidor® in HPLC-grade acetone at the highest tested concentration for each species (1563 ng·µL−1 for E. larvarum, 16 ng·µL−1 for S. rueppellii, and 43,750 ng·µL−1 for E. aeneus) and then we applied serial dilutions to achieve the appropriate doses.

2.3. Experimental Design, Exposure, and Mortality Assessment

For the three fly species, two negative control groups were included: negative control (untreated) and solvent control (pure acetone). To obtain imidacloprid dose–response curves, we exposed flies to six/seven doses in a geometric series. Following the preliminary range-finding trials, we used a dilution factor of 2.5 (5 imidacloprid doses from 40 to 1563 ng insect−1) in E. larvarum, 2.5 (5 doses from 0.41 to 16 ng insect−1) in S. rueppellii, and 3 (6 doses from 540 to 43,750 ng insect−1) in E. aeneus. After meeting, adult female flies were individually transferred to petri dishes and anesthetized on ice blocks for 2–3 min at 4–5 °C, and then exposed to the imidacloprid doses by topical application with 1 µL per female. The dose was applied on the dorsum of the thorax between the neck and the wing base by using a Gilson® Repetman automatic precision micropipette (Middleton, WI USA). Mated females of the solvent control were exposed to pure acetone. In the negative control, females were anesthetized as in the other groups, but the exposure phase was only simulated Mortality was recorded at 4, 24, 48, 72 and 96 h after exposure and compared with solvent control values. The number of replicates (single females) for the estimation of the dose–response curves was ca. 30 per dose.

2.4. Sublethal Effect Assessment: Reproduction (Oviposition Rate and Fecundity)

In the same live females used for the assessment of mortality, fecundity was evaluated 24 h after the application of the insecticide. Fecundity assessment was adjusted to the life history of each species. In E. larvarum, we counted the number of eggs laid on host larvae. Mated females in groups of five were introduced in Plexiglas cages (20 cm × 20 cm × 20 cm) with last instar G. mellonella larvae (three larvae per female) for 1 h [67]. After the oviposition period, the overall number of eggs counted on the host larvae was divided by the number of females in each cage to determine fecundity, expressed as the average number of eggs laid per female in 1 h. Contrary to the other two species, oviposition rate (proportion of females that laid eggs) in E. larvarum was assessed at the cage level (five females per cage) rather than at the individual level. In S. rueppellii, eggs laid by each individual female were counted daily until the end of the test (96 h). In each cage, we added two sprouts of pea (Pisum sativum L.) infested with aphids (M. persicae) and one non-infested sprout to provide a free surface for egg laying. Aphid populations were maintained on pea plants in a climate chamber at 20 ± 1 °C, 60–80% relative humidity and 16:8 h (light:dark) photoperiod. Sprout roots were wrapped in humid (distilled water) cotton and aluminum foil to keep the plant alive. Sprouts and aphids were substituted as needed. Eggs were counted daily and hoverfly larvae were removed as detected to avoid egg predation [68]. Fecundity of individual E. aeneus females was assessed daily for 7 days. To promote egg laying, we provided a substrate of decaying soaked oat grains [69,70] for each individual female. The substrate was prepared with organic oat and distilled water (200 g per 175 mL) and left for 24 h at room temperature in complete darkness. Each female was provided with ca. 9.3 g of substrate, which was renewed every two days to avoid mold formation. Fecundity in S. rueppellii and E. aeneus was expressed as the daily number of eggs laid per female. Simple sizes (N) were 16–28 for S. rueppellii, 6–17 (groups of 5 flies) for E. larvarum, and 21–29 for E. aeneus.

2.5. Ecotoxicological Data from Literature

Median lethal doses (LD50; expressed in ng/insect and in µg/g of insect) of imidacloprid for several bee and beetle species were obtained from the literature (see Supplementary Materials, Table S1). Unfortunately, we could not find data for other pollinator groups (butterflies and wasps). To increase the number of species available for the analysis, LD50 values from both active ingredient and commercial products were included. Although the presence of co-formulants may influence the LD50 values [71,72], inspection of the data gathered did not show any consistent trend when both LD50 estimates based on formulated products and on the active ingredient were available for the same species (Table S1). LD50 active ingredient estimates for Apis mellifera (n = 6) show a high level of variability, indicating that other factors such as population origin, test conditions and laboratory procedures may also influence the results. Where more than one value per species was available, we used the mean to build the SSD curve.

2.6. Statistical Analysis

Dose–response models were fitted to mortality data at the different exposure doses for each species using the drc package [73] in R ver. 4.1.2 [74]. Mortality was corrected with Abbott’s formula [75] using the untreated control as a reference. Differences in mortality between the untreated and solvent control were analyzed with chi-squared tests, and minimal detectable differences (MDDs) were calculated with the propTestMdd function of the EnvStats package (considering a 0.80 associated power and a 0.05 p-value). The mselect function of drc was used to determine which model function was most appropriate based on the Akaike information criterion (AIC). Due to the non-monotonic response of E. aeneus, in this species, the dose–response curve was divided into the two ascending parts [76]. Then, the median lethal dose (LD50) values at 48 h were calculated with the ED function of the drc package. Weight-normalized LD50 values were calculated based on the mean fresh weight of each species (Figure S1).
The species sensitivity distribution (SSD) was fitted to a log-normal dataset with the values of 48 h LD50 using the package ssdtools [77]. From the resulting curve, we obtained the 5% (HD5, as the lower limit of the distribution) hazardous dose, and calculated the 95% confidence intervals (CIs, 1000 interactions). Using the approach of Arena and Sgolastra [28], we calculated the sensitivity ratio (R) between A. mellifera and the other pollinator species:
R = L D 50   A .   m e l l i f e r a L D 50   p o l l i n a t o r
The reproduction endpoints were assessed only in individuals that survived at least 24 h after insecticide exposure. For E. larvarum, we assessed the number of eggs laid during 1 h. For S. rueppellii and E. aeneus, eggs were counted daily for 4 and 7 days, respectively. Firstly, we verified for both oviposition and fecundity that there were no significant differences between the two controls (untreated and solvent control). We use chi-squared and Mann–Whitney U tests for oviposition and reproduction, respectively. Minimal detectable differences (MDDs) on parametric tests were assessed as explained above. When no significant differences were found, solvent controls were used in the comparison with treated groups [78]. Then, the effect of treatment on oviposition rate (proportion of females that laid eggs) was examined with generalized linear models (GLMs) by fitting a binomial error distribution using the function glm. Gaussian error distribution was used to detect differences among treatments in fecundity (number of eggs laid per female). Data were log-transformed prior to analysis to achieve normal error distribution. All results are reported with a significance level of 5%.
The sublethal sensitivity index (SSI) is the ratio between the median lethal dose (LD50) and the no-observed-effect level (NOEL) for specific endpoints [62]. However, NOEL is highly dependent on sample size and dose selection [79]. For this reason, the lower confidence bound of the benchmark dose (BMDL) (i.e., the estimated lowest dose that produces an adverse response compared to the negative control) is often used instead of NOEL [79]. The BMD is the dose, derived from the estimated dose–response curve, associated with a given BMR (benchmark response, also known as critical effect size). The BMR was set at 10%, as recommended for quantal data analysis. In our study, BMDL values were estimated with the US EPA’s Benchmark Dose Software (BMDS 3.3.2 online version) [80] using the number of females that laid eggs (oviposition rate) for each treatment, and SSI was calculated for each species. High SSI values (>1) indicate that reproduction is a sensitive endpoint because it is inhibited at doses much lower than the LD50. Low SSI values (<1) indicate that the organism maintains its investment in reproduction until death.

3. Results

3.1. Species Sensitivity Distribution: Mortality

Mortality of the untreated control at 48 h, used to correct mortality of the other treatments, was 17% in E. larvarum, 0% in S. rueppellii, and 21% in E. aeneus. No significant differences were observed between the untreated and solvent control groups (E. larvarum: χ2 = 0.15, d.f. = 1, p = 0.70, MDD: 38%; S. rueppellii: χ2 = 1.48, d.f. = 1, p = 0.22, MDD: 29%; E. aeneus: χ2 = 0.083, d.f. = 1, p = 0.78, MDD: 36%). Imidacloprid dose–response curves varied substantially between species and the LD50 at 48 h increased in the following order: S. rueppellii < E. larvarum < E. aeneus (Table 1). This sensitivity ranking was maintained when LD50 values were corrected by the fresh body weight of each species (LD50 expressed in µg divided the average body weight expressed in grams). The dose–response curve of E. aeneus was non-monotonic and we could only calculate the LD50 associated with the first ascending part of the curve. Dose–response model parameters and LD50 values of the three dipteran species are reported in Table 1 and plots in the Supplementary Material (Figure S2).
We fitted the 48 h LD50 of our three target species and that of eight bee and two beetle species (obtained from the literature) into the SSD curve. The resulting HD5s (5% hazardous doses) were 0.615 ng/insect (95% CI [0.087; 7.416]) and 0.0105 µg/g of insect body weight (95% CI [0.0016; 0.118]) (Figure 1). The SSD curve shows that the sensitivity of A. mellifera is intermediate. That is, 59% (41% when body weight is accounted for) of the species are not protected by the honey bee LD50. Based on sensitivity ratio (R), the 10-fold safety factor recommended by the EFSA [27] to extrapolate the sensitivity of A. mellifera to other pollinators, was protective for 75% of the tested species, including the three dipterans, but not for three bee species (Apis cerana Fab., Melipona scutellaris Latreille and Leioproctus paahaumaa Donovan; Table S2).

3.2. Sublethal Effects: Reproduction (Oviposition Rate and Fecundity)

Oviposition rate (number of females that laid eggs) showed no differences between untreated and solvent control females in any of the species tested (S. rueppellii: χ2 = 3.4; d.f. = 1, p = 0.064, MDD: 37%; E. aeneus: χ2 = 0.42, d.f. = 1, p = 0.518, MDD: 39%). Analyses for E. larvarum were not conducted because the data were identical in both treatments. Fecundity (number of eggs laid per female) did not differ between untreated and solvent control females in any of the species tested either (E. larvarum: U = 43.0, d.f. = 1, p = 0.86; S. rueppellii: U = 276.5, d.f. = 1, p = 0.090; E. aeneus: U = 503.5, d.f. = 1, p = 0.214). For these reasons, only the solvent control was included in subsequent analyses. In E. larvarum, no differences were observed between treatments in oviposition (GLM: χ2 = 9.35 d.f. = 5, p = 0.090; Figure 2a). However, average eggs per female differed across treatments (GLM: χ2 = 122.34, d.f. = 5, p < 0.001; Figure 2a). All doses were significantly different from the control group (Figure 2a). In S. rueppellii, differences were observed among treatments in both oviposition (GLM: χ2 = 13.04, d.f. = 5, p = 0.023; Figure 2b) and fecundity (GLM: χ2 = 12.85, d.f. = 5, p = 0.025; Figure 2b). At the highest dose (i.e., 16 ng/insect) both endpoints were reduced compared to the control (Figure 2b). In E. aeneus, differences between treatments in oviposition were detected (GLM, χ2 = 12.68, d.f. = 6, p = 0.048, Figure 2c). The proportion of females that laid eggs was reduced both at 1620 and 29,166 ng/insects respect to the control (Figure 2c). In this species, fecundity was also different between treatments (GLM: χ2 = 18.17, d.f. = 6, p = 0.005; Figure 2c). All doses, with the exception of 540 ng/insect, were different from the control (Figure 2c). The BMDL calculated on oviposition (proportion of females laying eggs) for the three target species are shown in Table 2. All SSIs were > 1, indicating that oviposition was inhibited by exposure to the insecticide (Table 3).

4. Discussion

Most ecotoxicological studies on pollinators focus on the western honey bee, A. mellifera, and to a lesser extent bumblebees and solitary bees. Information on the sensitivity of other pollinator groups is mostly lacking, and when available, results are often difficult to compare with results of bee studies due to the use of different methodologies (i.e., oral bioassay [41], test on larvae [43]). Here, we provide relevant ecotoxicological endpoints for three dipteran pollinator species covering a range of life history traits and larval feeding habits, including a parasitoid, a predator, and a saprophagous species. As far as we know, this is the first study testing acute contact toxicity in dipteran pollinators using a risk assessment standard methodology. This is important because environmental risk assessment for pollinators is increasingly shifting towards a system-based approach integrating multiple species [40].
Among the three species tested, the aphidophagous hoverfly S. rueppellii exhibited the highest sensitivity to imidacloprid. High sensitivity of this species to this neonicotinoid was also reported in a previous study [97], in which adults were fed honeydew laced with imidacloprid at concentrations (15.6 ng/mL) detected in treated plants. On the other hand, the saprophagous hoverfly E. aeneus was by far the least sensitive species. This result is in line with other studies that observed a relatively high tolerance of another Eristalini, Eristalis tenax (L.), to imidacloprid and another neonicotinoid, thiamethoxam [41,43]. Interestingly, we obtained a non-monotonic dose response in E. aeneus. At intermediate doses, the toxicity dropped and started to increase again at higher doses. A similar biphasic hormetic dose–response pattern has been previously documented in several insects following exposure to neonicotinoids, including imidacloprid [76,98,99,100,101]. Hormetic response emerges when organisms mitigate stressor effects at low levels, often at the expense of other physiological processes. However, this balance becomes challenging to sustain at higher doses [102,103]. Numerous specific biochemical or molecular pathways could potentially contribute to hormesis, such as those associated with detoxification enzyme networks, oxidative stress response, and hormone signaling cascades [104,105,106,107]. Nonetheless, the precise reasons for the non-linear variation in neonicotinoid toxicity remain unclear and should be further explored.
The SSD curve shows a wide range of sensitivity to imidacloprid, which can be explained by differences in binding properties of the nicotinic acetylcholine receptors among species [108,109,110]. Interspecific variation in pesticide sensitivity is often phylogenetically constrained [111] and can also be partially explained by differences in body size [31,89] and other biological traits (e.g., hemolymph pH, lipid content, cuticle thickness). Especially when accounting for body weight, Coleoptera and Diptera species appear to be more tolerant to imidacloprid than bees. In this regard, toxicogenomics has strong potential to elucidate the molecular mechanisms evolved by different phylogenetic groups to respond to different toxicants. Sphaerophoria rueppellii, the most sensitive non-bee species in our SSD analysis, has a greater number of genes in various detoxification families, such as GSTs and CYP450s, in comparison with bees [112]. Some studies indicate that this expansion of detoxification genes in hoverflies may have evolved as an adaptation to exposure to environments with highly heterogeneous chemical backgrounds [112,113,114,115].
Our findings have important implications for environmental risk assessment. Bee pesticide risk assessment schemes rely almost exclusively on the western honey bee, A. mellifera, even though various studies have shown that this species is often less sensitive to pesticides than other bees [31,116]. For this reason, attempts to include other bee species in risk assessment schemes [27] have recently resulted in a roadmap towards a more holistic environmental risk assessment inclusive of all insect pollinator groups [40]. Species sensitivity distribution (SSD) curves are a useful approach to obtain estimates of pesticide sensitivity across species and to predict the proportion of potentially affected species at different concentrations using a limited set of ecotoxicity data [117]. Our results show that 60% of the species for which data are available, including the hoverfly S. rueppellii, fall below the honey bee LD50 threshold. When accounting for body weight, the level of protection derived from the honey bee LD50 increases, but still 40% of the species, all of them bees, remain unprotected. SSD curves can also be used to establish dose/concentration thresholds that would protect the 95th percentile of the species (HD5). In our study, the HD5 was approximately 100-fold lower than the LD50 of the honey bee. This value is 10 times higher than the safety factor of 10 recommended by the European Food Safety Authority [27] to estimate the sensitivity of untested bee species using honey bee data. Three (A. cerana, M. scutellaris and L. paahaumaa) of the 13 tested species in our dataset would have LD50 values lower than the honey bee LD50 even when applying a correction factor of 10. This outcome is in line with other studies emphasizing the need to apply safety factors higher than 10 to protect non-Apis bees [37,38]. The potential of the SSD approach is limited by the number and phylogenetic composition of species for which reliable LD50 values are available. Our study underscores the need to include a range of model species in pollinator and non-target arthropod risk assessment, including underrepresented insect orders such as Diptera and Lepidoptera [40]. Importantly, similar experimental procedures should be used across taxa so as to reduce the potential impact of confounding factors in the assessment of pesticide sensitivity. This is also important at the species level, as evidenced by the high variability found in the literature for honey bee LD50 estimates. To overcome this potential shortcoming, it is important to include a toxic reference in the analysis. Toxic references in risk assessment are used as a mortality reference and to validate the response of the test system by ensuring exposure consistency and the sensitivity of the test. Dimethoate LD50 is used as a toxic reference in bee risk assessment schemes, but no information is available for other pollinator species. Therefore, in our study, the results of a positive control could not be compared to a reference standard. Future studies should address the establishment of toxic reference values for non-bee pollinators.
Regulatory bodies also need repeatable methods to generate reliable data on intrinsic sensitivity. However, the ultimate impact of a pesticide will also depend on population resilience [118]. The latter may be defined as the capacity of the population to recover following disturbance, which is directly linked to oviposition and fecundity. Declines in these parameters can be considered early warning signals (EWSs) before population collapse. We assessed the effect of imidacloprid exposure on oviposition and fecundity in the three dipteran species. Again, S. rueppellii was the most sensitive species, followed by E. larvarum and E. aeneus. The lower confidence bound of the benchmark dose (BMDL) for oviposition in this species was 28 times lower than that of E. larvarum and 210 times lower than that of E. aeneus. Different species show different strategies and/or levels of adaptive plasticity to allocate resources between individual maintenance (survival) and reproduction. The sublethal sensitivity index (SSI) is a useful approach to predict the ecological effects of a pesticide, where a high SSI value indicates a high level of vulnerability. In agreement with results obtained in other invertebrate groups [62], SSI values for the three dipterans analyzed in our study were higher than 1, indicating that oviposition is a sensitive endpoint in these species, and therefore represents a useful EWS. The sublethal (BMDL) doses were 6 to 50 times lower than the lethal (LD50) doses, indicating that these species reduce the investment in oviposition in the presence of a pesticide stressor.

5. Conclusions

In conclusion, our study provides new information on the sensitivity of dipteran species and demonstrates that it is possible to apply standardized bee ecotoxicological procedures to other pollinator groups. Future studies should explore the mechanisms that cause non-monotonic dose–response relationships and create standardized toxicity testing techniques for dipteran pollinators. The integration of our results into SSD curves represents a first step towards an advanced environmental risk assessment, covering an expanded pool of species and including non-bee pollinators. The possibility of measuring oviposition and fecundity (and fertility in future studies) in some dipteran species under laboratory conditions, together with the use of SSI, widens our ability to predict the impact of pesticides at the population level and provides a crucial first step for the integration of population dynamics in simulation models for environmental risk assessment [119].

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/insects15050317/s1. Figure S1: Boxplots of fresh weight for the three fly species.; Table S1: Toxicity data used to build the species sensitivity distribution (SSD) curve. LD50 values: median lethal dose of imidacloprid at 48 h. Where no weight standardization was performed in the original study, we indicate the reference for the mean weight used to transform the data.; Table S2: Species sensitivity ratio (R). Values in bold exceed the range of 10-fold safety factor from the endpoint of A. mellifera, recommended by the EFSA (2013); Figure S2: Fitted dose–response curves in (A) Exorista larvarum (L.), (B) Sphaerophoria rueppellii (Wiedemann) (Diptera: Syrphidae) and (C) Eristalinus aeneus (Scopoli) (Diptera: Syrphidae). Ref. [120] is cited in Supplementary Materials.

Author Contributions

F.S., J.B., M.L.D., G.B. and D.S. conceived the experiments. C.A.H.M., C.A., F.S., J.B., M.L.D., G.B. and D.S. designed the experiments. C.A.H.M., C.A. and S.F. collected the data. C.A.H.M. and C.A. analyzed the data. C.A.H.M. and F.S. took the lead in writing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by a grant from the University of Bologna (RFO2021_SGOLASTRA) and by the European Union NextGenerationEU program under the National Recovery and Resilience Plan (NRRP), project title “Agritech, National Research Center,” project code CN_00000022 to Fabio Sgolastra and Maria Luisa Dindo. A doctoral scholarship from the Italian Ministry of University and Research supported Cátia A.H. Martins, and Celeste Azpiazu was supported by a Margarita Salas postdoctoral fellowship from the Spanish Ministry of Universities through the EU NextGeneration program.

Data Availability Statement

Data will be made available on request.

Acknowledgments

We thank several students, especially Gioia De Iulis, Giovanni Dentale and Federico Agostini, for their help with laboratory data collection. We are grateful to Marc Vaez-Olivera and Manuela Sanchez Fernandez (POLIFLY) and Apostolos Pekas (BioBest) for providing the hoverfly populations. We thank R. Molowny-Horas for statistical advice.

Conflicts of Interest

We have no competing interests.

References

  1. Klein, A.-M.; Vaissière, B.E.; Cane, J.H.; Steffan-Dewenter, I.; Cunningham, S.A.; Kremen, C.; Tscharntke, T. Importance of Pollinators in Changing Landscapes for World Crops. Proc. R. Soc. B Biol. Sci. 2007, 274, 303–313. [Google Scholar] [CrossRef] [PubMed]
  2. Ollerton, J.; Winfree, R.; Tarrant, S. How Many Flowering Plants Are Pollinated by Animals? Oikos 2011, 120, 321–326. [Google Scholar] [CrossRef]
  3. Rader, R.; Bartomeus, I.; Garibaldi, L.A.; Garratt, M.P.D.; Howlett, B.G.; Winfree, R.; Cunningham, S.A.; Mayfield, M.M.; Arthur, A.D.; Andersson, G.K.S.; et al. Non-Bee Insects Are Important Contributors to Global Crop Pollination. Proc. Natl. Acad. Sci. USA 2016, 113, 146–151. [Google Scholar] [CrossRef] [PubMed]
  4. Skevington, J.H.; Dang, P.T. Exploring the Diversity of Flies (Diptera). Biodiversity 2002, 3, 3–27. [Google Scholar] [CrossRef]
  5. Ssymank, A.; Kearns, C.A.; Pape, T.; Thompson, F.C. Pollinating Flies (Diptera): A Major Contribution to Plant Diversity and Agricultural Production. Biodiversity 2008, 9, 86–89. [Google Scholar] [CrossRef]
  6. Courtney, G.W.; Pape, T.; Skevington, J.H.; Sinclair, B.J. Biodiversity of Diptera. In Insect Biodiversity; Wiley: Hoboken, NJ, USA, 2017; pp. 229–278. [Google Scholar]
  7. Woodcock, T.S.; Larson, B.M.H.; Kevan, P.G.; Inouye, D.W.; Lunau, K. Flies and Flowers II: Floral Attractants and Rewards. J. Pollinat. Ecol. 2014, 12, 63–94. [Google Scholar] [CrossRef]
  8. Inouye, D.W.; Larson, B.M.H.; Ssymank, A.; Kevan, P.G. Flies and Flowers III: Ecology of Foraging and Pollination. J. Pollinat. Ecol. 2015, 16, 115–133. [Google Scholar] [CrossRef]
  9. Roquer-Beni, L.; Arnan, X.; Rodrigo, A.; Bosch, J. What Makes a Good Pollinator? Relationship between Pollinator Traits and Pollination Effectiveness in Apple Flowers. Entomol. Gen. 2022, 42, 875–882. [Google Scholar] [CrossRef]
  10. Rodríguez-Gasol, N.; Alins, G.; Veronesi, E.R.; Wratten, S. The Ecology of Predatory Hoverflies as Ecosystem-Service Providers in Agricultural Systems. Biol. Control 2020, 151, 104405. [Google Scholar] [CrossRef]
  11. Larson, B.M.H.; Kevan, P.G.; Inouye, D.W. Flies and Flowers: Taxonomic Diversity of Anthophiles and Pollinators. Can. Entomol. 2001, 133, 439–465. [Google Scholar] [CrossRef]
  12. Raguso, R.A. Don’t Forget the Flies: Dipteran Diversity and Its Consequences for Floral Ecology and Evolution. Appl. Entomol. Zool. 2020, 55, 1–7. [Google Scholar] [CrossRef]
  13. Rader, R.; Cunningham, S.A.; Howlett, B.G.; Inouye, D.W. Non-Bee Insects as Visitors and Pollinators of Crops: Biology, Ecology, and Management. Annu. Rev. Entomol. 2020, 65, 391–407. [Google Scholar] [CrossRef] [PubMed]
  14. Al-Dobai, S.; Reitz, S.; Sivinski, J. Tachinidae (Diptera) Associated with Flowering Plants: Estimating Floral Attractiveness. Biol. Control 2012, 61, 230–239. [Google Scholar] [CrossRef]
  15. Orford, K.A.; Vaughan, I.P.; Memmott, J. The Forgotten Flies: The Importance of Non-Syrphid Diptera as Pollinators. Proc. R. Soc. B Biol. Sci. 2015, 282, 20142934. [Google Scholar] [CrossRef] [PubMed]
  16. Dunn, L.; Lequerica, M.; Reid, C.R.; Latty, T. Dual Ecosystem Services of Syrphid Flies (Diptera: Syrphidae): Pollinators and Biological Control Agents. Pest Manag. Sci. 2020, 76, 1973–1979. [Google Scholar] [CrossRef]
  17. Dindo, M.L.; Grenier, S. Production of Dipteran Parasitoids. In Mass Production of Beneficial Organisms; Elsevier: Amsterdam, The Netherlands, 2023; pp. 71–100. [Google Scholar]
  18. Biesmeijer, J.C.; Roberts, S.P.M.; Reemer, M.; Ohlemuller, R.; Edwards, M.; Peeters, T.; Schaffers, A.P.; Potts, S.G.; Kleukers, R.; Thomas, C.D.; et al. Parallel Declines in Pollinators and Insect-Pollinated Plants in Britain and the Netherlands. Science 2006, 313, 351–354. [Google Scholar] [CrossRef]
  19. Zattara, E.E.; Aizen, M.A. Worldwide Occurrence Records Suggest a Global Decline in Bee Species Richness. One Earth 2021, 4, 114–123. [Google Scholar] [CrossRef]
  20. Hallmann, C.A.; Sorg, M.; Jongejans, E.; Siepel, H.; Hofland, N.; Schwan, H.; Stenmans, W.; Müller, A.; Sumser, H.; Hörren, T.; et al. More than 75 Percent Decline over 27 Years in Total Flying Insect Biomass in Protected Areas. PLoS ONE 2017, 12, e0185809. [Google Scholar] [CrossRef]
  21. Hallmann, C.A.; Ssymank, A.; Sorg, M.; de Kroon, H.; Jongejans, E. Insect Biomass Decline Scaled to Species Diversity: General Patterns Derived from a Hoverfly Community. Proc. Natl. Acad. Sci. USA 2021, 118, e2002554117. [Google Scholar] [CrossRef]
  22. Barendregt, A.; Zeegers, T.; van Steenis, W.; Jongejans, E. Forest Hoverfly Community Collapse: Abundance and Species Richness Drop over Four Decades. Insect Conserv. Divers. 2022, 15, 510–521. [Google Scholar] [CrossRef]
  23. 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]
  24. Rundlöf, M.; Andersson, G.K.S.; Bommarco, R.; Fries, I.; Hederström, V.; Herbertsson, L.; Jonsson, O.; Klatt, B.K.; Pedersen, T.R.; Yourstone, J.; et al. Seed Coating with a Neonicotinoid Insecticide Negatively Affects Wild Bees. Nature 2015, 521, 77–80. [Google Scholar] [CrossRef]
  25. 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]
  26. U.S. Environmental Protection Agency (USEPA); Health Canada Pest Management Regulatory Agency (PMRA); California Departament of Pesticide Regulatio. USEPA Guidance for Assessing Pesticide Risks to Bees, USEPA: Washington, DC, USA, 2014.
  27. EFSA Guidance on the Risk Assessment of Plant Protection Products on Bees (Apis Mellifera, Bombus Spp. and Solitary Bees). EFSA J. 2013, 11, 3295. [CrossRef]
  28. Sgolastra, F.; Medrzycki, P.; Bortolotti, L.; Maini, S.; Porrini, C.; Simon-Delso, N.; Bosch, J. Bees and Pesticide Regulation: Lessons from the Neonicotinoid Experience. Biol. Conserv. 2020, 241, 108356. [Google Scholar] [CrossRef]
  29. Sanchez-Bayo, F.; Goka, K. Pesticide Residues and Bees–A Risk Assessment. PLoS ONE 2014, 9, e94482. [Google Scholar] [CrossRef]
  30. Azpiazu, C.; Bosch, J.; Bortolotti, L.; Medrzycki, P.; Teper, D.; Molowny-Horas, R.; Sgolastra, F. Toxicity of the Insecticide Sulfoxaflor Alone and in Combination with the Fungicide Fluxapyroxad in Three Bee Species. Sci. Rep. 2021, 11, 6821. [Google Scholar] [CrossRef]
  31. Arena, M.; Sgolastra, F. A Meta-Analysis Comparing the Sensitivity of Bees to Pesticides. Ecotoxicology 2014, 23, 324–334. [Google Scholar] [CrossRef]
  32. Sgolastra, F.; Medrzycki, P.; Bortolotti, L.; Renzi, M.T.; Tosi, S.; Bogo, G.; Teper, D.; Porrini, C.; Molowny-Horas, R.; Bosch, J. Synergistic Mortality between a Neonicotinoid Insecticide and an Ergosterol-Biosynthesis-Inhibiting Fungicide in Three Bee Species. Pest Manag. Sci. 2017, 73, 1236–1243. [Google Scholar] [CrossRef]
  33. Biddinger, D.J.; Robertson, J.L.; Mullin, C.; Frazier, J.; Ashcraft, S.A.; Rajotte, E.G.; Joshi, N.K.; Vaughn, M. Comparative Toxicities and Synergism of Apple Orchard Pesticides to Apis Mellifera (L.) and Osmia Cornifrons (Radoszkowski). PLoS ONE 2013, 8, e72587. [Google Scholar] [CrossRef]
  34. Heard, M.S.; Baas, J.; Dorne, J.-L.; Lahive, E.; Robinson, A.G.; Rortais, A.; Spurgeon, D.J.; Svendsen, C.; Hesketh, H. Comparative Toxicity of Pesticides and Environmental Contaminants in Bees: Are Honey Bees a Useful Proxy for Wild Bee Species? Sci. Total Environ. 2017, 578, 357–365. [Google Scholar] [CrossRef] [PubMed]
  35. Robinson, A.; Hesketh, H.; Lahive, E.; Horton, A.A.; Svendsen, C.; Rortais, A.; Dorne, J.L.; Baas, J.; Heard, M.S.; Spurgeon, D.J. Comparing Bee Species Responses to Chemical Mixtures: Common Response Patterns? PLoS ONE 2017, 12, e0176289. [Google Scholar] [CrossRef] [PubMed]
  36. Linguadoca, A.; Jürison, M.; Hellström, S.; Straw, E.A.; Šima, P.; Karise, R.; Costa, C.; Serra, G.; Colombo, R.; Paxton, R.J.; et al. Intra-Specific Variation in Sensitivity of Bombus Terrestris and Osmia Bicornis to Three Pesticides. Sci. Rep. 2022, 12, 17311. [Google Scholar] [CrossRef] [PubMed]
  37. de Assis, J.C.; Tadei, R.; Menezes-Oliveira, V.B.; Silva-Zacarin, E.C.M. Are Native Bees in Brazil at Risk from the Exposure to the Neonicotinoid Imidacloprid? Environ. Res. 2022, 212, 113127. [Google Scholar] [CrossRef] [PubMed]
  38. Lourencetti, A.P.S.; Azevedo, P.; Miotelo, L.; Malaspina, O.; Nocelli, R.C.F. Surrogate Species in Pesticide Risk Assessments: Toxicological Data of Three Stingless Bees Species. Environ. Pollut. 2023, 318, 120842. [Google Scholar] [CrossRef] [PubMed]
  39. Beneficial Arthropod Regulatory Testing Group; European and Mediterranean Plant Protection Organisation; Council of Europe; Organisation for Economic Co-operation and Development; International Organisation for Biological and Integrated Control of Noxious Animals and Plants; Society of Environmental Toxicology and Chemistry-Europe; Commission of the European Communities. Guidance Document on Regulatory Testing and Risk Assessment Procedures for Plant Protection Products with Non-Target Arthropods: From the ESCORT 2 Workshop (European Standard Characteristics of Non-Target Arthropod Regulatory Testing), Candolfi, M.P., Barrett, K.L., Campbell, P.J., Forster, R., Grandy, N., Huet, M.-C., Lewis, G., Oomen, P.A., Schmuck, R., Vog, H., Eds.; SETAC: Brussels, Belgium, 2001; ISBN 9781880611524.
  40. Williams, J.H.; Bordoni, A.; Bednarska, A.; Pinto, A.; Martins, C.A.H.; Henriques, D.; Sgolastra, F.; Knapp, J.; Loureiro, J.; Sousa, J.P.; et al. Roadmap for Action on the Environmental Risk Assessment of Chemicals for Insect Pollinators (IPol-ERA). EFSA Support. Publ. 2023, 20, 8431E. [Google Scholar] [CrossRef]
  41. Basley, K.; Davenport, B.; Vogiatzis, K.; Goulson, D. Effects of Chronic Exposure to Thiamethoxam on Larvae of the Hoverfly Eristalis Tenax (Diptera, Syrphidae). PeerJ 2018, 6, e4258. [Google Scholar] [CrossRef] [PubMed]
  42. Moens, J.; De Clercq, P.; Tirry, L. Side Effects of Pesticides on the Larvae of the Hoverfly Episyrphus Balteatus in the Laboratory. Phytoparasitica 2011, 39, 1–9. [Google Scholar] [CrossRef]
  43. Nagloo, N.; Rigosi, E.; O’Carroll, D.C. Acute and Chronic Toxicity of Imidacloprid in the Pollinator Fly, Eristalis Tenax L., Assessed Using a Novel Oral Bioassay. Ecotoxicol. Environ. Saf. 2023, 251, 114505. [Google Scholar] [CrossRef]
  44. Benelli, M.; Tóth, F.; Dindo, M.L. Low-temperature Storage of Exorista Larvarum Puparia as a Tool for Assisting Parasitoid Production. Entomol. Exp. Appl. 2018, 166, 914–924. [Google Scholar] [CrossRef]
  45. Mellini, E.; Coulibaly, A.K. Un Decennio Di Sperimentazione Sul Sistema Ospite-Parassita Galleria Mellonella L.-Pseudogonia rufifrons Wied: Sintesi dei risultati. Boll. Ist. Ent. “G. Grandi” Univ. Bologna 1991, 45, 191–249. [Google Scholar]
  46. Amorós-Jiménez, R.; Pineda, A.; Fereres, A.; Marcos-García, M.Á. Prey Availability and Abiotic Requirements of Immature Stages of the Aphid Predator Sphaerophoria Rueppellii. Biol. Control 2012, 63, 17–24. [Google Scholar] [CrossRef]
  47. Amorós-Jiménez, R.; Pineda, A.; Fereres, A.; Marcos-García, M.Á. Feeding Preferences of the Aphidophagous Hoverfly Sphaerophoria Rueppellii Affect the Performance of Its Offspring. BioControl 2014, 59, 427–435. [Google Scholar] [CrossRef]
  48. Pekas, A.; De Craecker, I.; Boonen, S.; Wäckers, F.L.; Moerkens, R. One Stone; Two Birds: Concurrent Pest Control and Pollination Services Provided by Aphidophagous Hoverflies. Biol. Control 2020, 149, 104328. [Google Scholar] [CrossRef]
  49. Sánchez, M.; Belliure, B.; Montserrat, M.; Gil, J.; Velásquez, Y. Pollination by the Hoverfly Eristalinus Aeneus (Diptera: Syrphidae) in Two Hybrid Seed Crops: Celery and Fennel (Apiaceae). J. Agric. Sci. 2022, 160, 194–206. [Google Scholar] [CrossRef]
  50. Simon-Delso, N.; Amaral-Rogers, V.; Belzunces, L.P.; Bonmatin, J.M.; Chagnon, M.; Downs, C.; Furlan, L.; Gibbons, D.W.; Giorio, C.; Girolami, V.; et al. Systemic Insecticides (Neonicotinoids and Fipronil): Trends, Uses, Mode of Action and Metabolites. Environ. Sci. Pollut. Res. 2015, 22, 5–34. [Google Scholar] [CrossRef]
  51. van der Sluijs, J.P.; Amaral-Rogers, V.; Belzunces, L.P.; Bijleveld van Lexmond, M.F.I.J.; Bonmatin, J.-M.; Chagnon, M.; Downs, C.A.; Furlan, L.; Gibbons, D.W.; Giorio, C.; et al. Conclusions of the Worldwide Integrated Assessment on the Risks of Neonicotinoids and Fipronil to Biodiversity and Ecosystem Functioning. Environ. Sci. Pollut. Res. 2015, 22, 148–154. [Google Scholar] [CrossRef]
  52. Maini, S.; Medrzycki, P.; Porrini, C. The Puzzle of Honey Bee Losses: A Brief Review. Bull. Insectology 2010, 63, 153–160. [Google Scholar]
  53. 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]
  54. OJEU Commission Implementing Regulation (EU) 2018/783 of 29 May 2018 Amending Implementing Regulation (EU) No 540/2011 as Regards the Conditions of Approval of the Active Substance Imidacloprid. Off. J. Eur. Union 2018, 132, 31–34.
  55. Goulson, D. Pesticides, Corporate Irresponsibility, and the Fate of Our Planet. One Earth 2020, 2, 302–305. [Google Scholar] [CrossRef]
  56. OECD. Test No. 214: Honeybees, Acute Contact Toxicity Test, OECD Guidelines for the Testing of Chemicals, Section 2; OECD: Paris, France, 1998; ISBN 9789264070189.
  57. OECD. Test No. 213: Honeybees, Acute Oral Toxicity Test, OECD Guidelines for the Testing of Chemicals, Section 2; OECD: Paris, France, 1998; ISBN 9789264070165.
  58. OECD. Test No. 246: Bumblebee, Acute Contact Toxicity Test, OECD Guidelines for the Testing of Chemicals, Section 2; OECD: Paris, France, 2017; ISBN 9789264284104.
  59. OECD. Test No. 247: Bumblebee, Acute Oral Toxicity Test, OECD Guidelines for the Testing of Chemicals, Section 2; OECD: Paris, France, 2017; ISBN 9789264284128.
  60. Posthuma, L.; Suter II, G.W.; Traas, T.P. Species Sensitivity Distributions in Ecotoxicology; CRC Press: Boca Raton, FL, USA, 2001; ISBN 1420032313. [Google Scholar]
  61. Wheeler, J.; Grist, E.P.; Leung, K.M.; Morritt, D.; Crane, M. Species Sensitivity Distributions: Data and Model Choice. Mar. Pollut. Bull. 2002, 45, 192–202. [Google Scholar] [CrossRef] [PubMed]
  62. van Straalen, N.M. Biodiversity of Ecotoxicological Responses in Animals. Netherlands J. Zool. 1994, 44, 112–129. [Google Scholar] [CrossRef]
  63. Rodrigues, M.A.; Flatt, T. Endocrine Uncoupling of the Trade-off between Reproduction and Somatic Maintenance in Eusocial Insects. Curr. Opin. Insect Sci. 2016, 16, 1–8. [Google Scholar] [CrossRef] [PubMed]
  64. Sgolastra, F.; Hinarejos, S.; Pitts-Singer, T.L.; Boyle, N.K.; Joseph, T.; Lūckmann, J.; Raine, N.E.; Singh, R.; Williams, N.M.; Bosch, J. Pesticide Exposure Assessment Paradigm for Solitary Bees. Environ. Entomol. 2019, 48, 22–35. [Google Scholar] [CrossRef] [PubMed]
  65. Dindo, M.L.; Modesto, M.; Rossi, C.; Di Vito, M.; Burgio, G.; Barbanti, L.; Mattarelli, P. Monarda Fistulosa Hydrolate as Antimicrobial Agent in Artificial Media for the in Vitro Rearing of the Tachinid Parasitoid Exorista Larvarum. Entomol. Exp. Appl. 2021, 169, 79–89. [Google Scholar] [CrossRef]
  66. Dindo, M.L.; Marchetti, E.; Baronio, P. In Vitro Rearing of the Parasitoid Exorista Larvarum (Diptera: Tachinidae) from Eggs Laid Out of Host. J. Econ. Entomol. 2007, 100, 26–30. [Google Scholar] [CrossRef] [PubMed]
  67. Dindo, M.L.; Rezaei, M.; De Clercq, P. Improvements in the Rearing of the Tachinid Parasitoid Exorista Larvarum (Diptera: Tachinidae): Influence of Adult Food on Female Longevity and Reproduction Capacity. J. Insect Sci. 2019, 19, 6. [Google Scholar] [CrossRef] [PubMed]
  68. Branquart, E.; Hemptinne, J.-L.; Bauffe, C.; Benfekih, L. Cannibalism In Episyrphus Balteatus (Ditp: Syrphidae). BioControl 1997, 42, 145–152. [Google Scholar] [CrossRef]
  69. Campoy, A.; Lutsyk, M.; Pérez-Bañón, C.; Rojo, S. Age-Stage Two-Sex Life Table Analysis of Eristalinus Aeneus (Diptera, Syrphidae) Reared with Two Different Larval Media. Bull. Entomol. Res. 2022, 112, 13–20. [Google Scholar] [CrossRef]
  70. Gladis, T. Laborzucht Einiger Eristalinen (Diptera, Syrphidae) Und Möglichkeiten Für Ihren Einsatz in Der Pflanzenzüchtung [Laboratory Rearing of Some Eristalines (Diptera, Syrphidae) and the Possibility of Their Use in Plant Cultures]. Verhandlungen der Westdtsch. Entomol. 1994, 1993, 139–152. [Google Scholar]
  71. Straw, E.A.; Carpentier, E.N.; Brown, M.J.F. Roundup Causes High Levels of Mortality Following Contact Exposure in Bumble Bees. J. Appl. Ecol. 2021, 58, 1167–1176. [Google Scholar] [CrossRef]
  72. Straw, E.A.; Brown, M.J.F. Co-Formulant in a Commercial Fungicide Product Causes Lethal and Sub-Lethal Effects in Bumble Bees. Sci. Rep. 2021, 11, 21653. [Google Scholar] [CrossRef] [PubMed]
  73. Ritz, C.; Baty, F.; Streibig, J.C.; Gerhard, D. Dose-Response Analysis Using R. PLoS ONE 2015, 10, e0146021. [Google Scholar] [CrossRef] [PubMed]
  74. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria; Available online: https://www.r-project.org/ (accessed on 30 March 2022).
  75. Abbott, W.S. A Method of Computing the Effectiveness of an Insecticide. J. Econ. Entomol 1925, 18, 265–267. [Google Scholar] [CrossRef]
  76. Suchail, S.; Guez, D.; Belzunces, L.P. Characteristics of Imidacloprid Toxicity in Two Apis Mellifera Subspecies. Environ. Toxicol. Chem. 2000, 19, 1901–1905. [Google Scholar] [CrossRef]
  77. Thorley, J.; Schwarz, C. Ssdtools: An R Package to Fit Species Sensitivity Distributions. J. Open Source Softw. 2018, 3, 1082. [Google Scholar] [CrossRef]
  78. OECD. Current Approaches in the Statistical Analysis of Ecotoxicity Data: A Guidance to Application, OECD Environment Health and Safety Publications Series on Testing and Assessment; OECD: Paris, France, 2016.
  79. Laskowski, R. Some Good Reasons to Ban the Use of NOEC, LOEC and Related Concepts in Ecotoxicology. Oikos 1995, 73, 140–144. [Google Scholar] [CrossRef]
  80. USEPA Benchmark Dose Tools (BMDS) Online. Available online: https://bmdsonline.epa.gov/ (accessed on 12 April 2023).
  81. Yasuda, M.; Sakamoto, Y.; Goka, K.; Nagamitsu, T.; Taki, H. Insecticide Susceptibility in Asian Honey Bees (Apis Cerana (Hymenoptera: Apidae)) and Implications for Wild Honey Bees in Asia. J. Econ. Entomol. 2017, 110, 447–452. [Google Scholar] [CrossRef]
  82. Soares, H.M.; Jacob, C.R.O.; Carvalho, S.M.; Nocelli, R.C.F.; Malaspina, O. Toxicity of Imidacloprid to the Stingless Bee Scaptotrigona Postica Latreille, 1807 (Hymenoptera: Apidae). Bull. Environ. Contam. Toxicol. 2015, 94, 675–680. [Google Scholar] [CrossRef]
  83. da Costa, L.M.; Grella, T.C.; Barbosa, R.A.; Malaspina, O.; Nocelli, R.C.F. Determination of Acute Lethal Doses (LD50 and LC50) of Imidacloprid for the Native Bee Melipona Scutellaris Latreille, 1811 (Hymenoptera: Apidae). Sociobiology 2015, 62. [Google Scholar] [CrossRef]
  84. Kueh Tai, F.; Pattemore, D.E.; Jochym, M.; Beggs, J.R.; Northcott, G.L.; Mortensen, A.N. Honey Bee Toxicological Responses Do Not Accurately Predict Environmental Risk of Imidacloprid to a Solitary Ground-Nesting Bee Species. Sci. Total Environ. 2022, 839, 156398. [Google Scholar] [CrossRef]
  85. Youn, Y.N.; Seo, M.J.; Shin, J.G.; Jang, C.; Yu, Y.M. Toxicity of Greenhouse Pesticides to Multicolored Asian Lady Beetles, Harmonia Axyridis (Coleoptera: Coccinellidae). Biol. Control 2003, 28, 164–170. [Google Scholar] [CrossRef]
  86. Lucas, É.; Giroux, S.; Demougeot, S.; Duchesne, R.-M.; Coderre, D. Compatibility of a Natural Enemy, Coleomegilla Maculata Lengi (Col., Coccinellidae) and Four Insecticides Used against the Colorado Potato Beetle (Col., Chrysomelidae). J. Appl. Entomol. 2004, 128, 233–239. [Google Scholar] [CrossRef]
  87. Hartfelder, K.; Engels, W. Allometric and Multivariate Analysis of Sex and Caste Polymorphism in the Neotropical Stingless Bee, Scaptotrigona Postica. Insectes Soc. 1992, 39, 251–266. [Google Scholar] [CrossRef]
  88. Lourenço, C.T.; Carvalho, S.M.; Malaspina, O.; Nocelli, R.C.F. Oral Toxicity of Fipronil Insecticide Against the Stingless Bee Melipona Scutellaris (Latreille, 1811). Bull. Environ. Contam. Toxicol. 2012, 89, 921–924. [Google Scholar] [CrossRef] [PubMed]
  89. Thompson, H. Extrapolation of Acute Toxicity across Bee Species. Integr. Environ. Assess. Manag. 2016, 12, 622–626. [Google Scholar] [CrossRef] [PubMed]
  90. Beadle, K.; Singh, K.S.; Troczka, B.J.; Randall, E.; Zaworra, M.; Zimmer, C.T.; Hayward, A.; Reid, R.; Kor, L.; Kohler, M.; et al. Genomic Insights into Neonicotinoid Sensitivity in the Solitary Bee Osmia Bicornis. PLoS Genet. 2019, 15, e1007903. [Google Scholar] [CrossRef] [PubMed]
  91. Phan, N.T.; Joshi, N.K.; Rajotte, E.G.; López-Uribe, M.M.; Zhu, F.; Biddinger, D.J. A New Ingestion Bioassay Protocol for Assessing Pesticide Toxicity to the Adult Japanese Orchard Bee (Osmia Cornifrons). Sci. Rep. 2020, 10, 9517. [Google Scholar] [CrossRef] [PubMed]
  92. Uhl, P.; Awanbor, O.; Schulz, R.S.; Brühl, C.A. Osmia Bicornis Is Rarely an Adequate Regulatory Surrogate Species. Comparing Its Acute Sensitivity towards Multiple Insecticides with Regulatory Apis Mellifera Endpoints. PLoS ONE 2019, 14, e0201081. [Google Scholar] [CrossRef]
  93. ECOTOX Curated Toxicity Data Were Retrieved from the ECOTOXicology Knowledgebase. U.S. Environmental Protection Agency. Available online: http://www.epa.gov/ecotox/ (accessed on 1 November 2022).
  94. Bortolotti, L.; Porrini, C.; Sbrenna, G. Effetti Dell’imidacloprid Nei Confronti Di Bombus Terrestris (L.). Prove Di Laboratorio. Inf. Fitopatol. 2002, 3, 66–71. [Google Scholar]
  95. Hagen, M.; Wikelski, M.; Kissling, W.D. Space Use of Bumblebees (Bombus Spp.) Revealed by Radio-Tracking. PLoS ONE 2011, 6, e19997. [Google Scholar] [CrossRef] [PubMed]
  96. Suchail, S.; Debrauwer, L.; Belzunces, L.P. Metabolism of Imidacloprid in Apis Mellifera. Pest Manag. Sci. 2004, 60, 291–296. [Google Scholar] [CrossRef]
  97. Calvo-Agudo, M.; González-Cabrera, J.; Picó, Y.; Calatayud-Vernich, P.; Urbaneja, A.; Dicke, M.; Tena, A. Neonicotinoids in Excretion Product of Phloem-Feeding Insects Kill Beneficial Insects. Proc. Natl. Acad. Sci. USA 2019, 116, 16817–16822. [Google Scholar] [CrossRef] [PubMed]
  98. 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–94, 56–63. [Google Scholar] [CrossRef] [PubMed]
  99. Anderson, N.L.; Harmon-Threatt, A.N. Chronic Contact with Realistic Soil Concentrations of Imidacloprid Affects the Mass, Immature Development Speed, and Adult Longevity of Solitary Bees. Sci. Rep. 2019, 9, 3724. [Google Scholar] [CrossRef] [PubMed]
  100. Potts, R.; Clarke, R.M.; Oldfield, S.E.; Wood, L.K.; Hempel de Ibarra, N.; 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]
  101. Haddi, K.; Mendes, M.V.; Barcellos, M.S.; Lino-Neto, J.; Freitas, H.L.; Guedes, R.N.C.; Oliveira, E.E. Sexual Success after Stress? Imidacloprid-Induced Hormesis in Males of the Neotropical Stink Bug Euschistus Heros. PLoS ONE 2016, 11, e0156616. [Google Scholar] [CrossRef] [PubMed]
  102. Agathokleous, E.; Calabrese, E.J. Environmental Toxicology and Ecotoxicology: How Clean Is Clean? Rethinking Dose-Response Analysis. Sci. Total Environ. 2020, 746, 138769. [Google Scholar] [CrossRef]
  103. Cutler, G.C.; Amichot, M.; Benelli, G.; Guedes, R.N.C.; Qu, Y.; Rix, R.R.; Ullah, F.; Desneux, N. Hormesis and Insects: Effects and Interactions in Agroecosystems. Sci. Total Environ. 2022, 825, 153899. [Google Scholar] [CrossRef]
  104. Calabrese, E.J. Hormetic Mechanisms. Crit. Rev. Toxicol. 2013, 43, 580–606. [Google Scholar] [CrossRef] [PubMed]
  105. De Smet, L.; Hatjina, F.; Ioannidis, P.; Hamamtzoglou, A.; Schoonvaere, K.; Francis, F.; Meeus, I.; Smagghe, G.; de Graaf, D.C. Stress Indicator Gene Expression Profiles, Colony Dynamics and Tissue Development of Honey Bees Exposed to Sub-Lethal Doses of Imidacloprid in Laboratory and Field Experiments. PLoS ONE 2017, 12, e0171529. [Google Scholar] [CrossRef] [PubMed]
  106. Derecka, K.; Blythe, M.J.; Malla, S.; Genereux, D.P.; Guffanti, A.; Pavan, P.; Moles, A.; Snart, C.; Ryder, T.; Ortori, C.A. Transient Exposure to Low Levels of Insecticide Affects Metabolic Networks of Honeybee Larvae. PLoS ONE 2013, 8, e68191. [Google Scholar] [CrossRef] [PubMed]
  107. Erofeeva, E.A. Environmental Hormesis: From Cell to Ecosystem. Curr. Opin. Environ. Sci. Heal. 2022, 29, 100378. [Google Scholar] [CrossRef]
  108. Crossthwaite, A.J.; Bigot, A.; Camblin, P.; Goodchild, J.; Lind, R.J.; Slater, R.; Maienfisch, P. The Invertebrate Pharmacology of Insecticides Acting at Nicotinic Acetylcholine Receptors. J. Pestic. Sci. 2017, 42, 67–83. [Google Scholar] [CrossRef]
  109. Jones, A.K.; Brown, L.A.; Sattelle, D.B. Insect Nicotinic Acetylcholine Receptor Gene Families: From Genetic Model Organism to Vector, Pest and Beneficial Species. Invertebr. Neurosci. 2007, 7, 67–73. [Google Scholar] [CrossRef]
  110. Maloney, E.M.; Taillebois, E.; Gilles, N.; Morrissey, C.A.; Liber, K.; Servent, D.; Thany, S.H. Binding Properties to Nicotinic Acetylcholine Receptors Can Explain Differential Toxicity of Neonicotinoid Insecticides in Chironomidae. Aquat. Toxicol. 2021, 230, 105701. [Google Scholar] [CrossRef]
  111. Pamminger, T. Extrapolating Acute Contact Bee Sensitivity to Insecticides Based on Body Weight Using a Phylogenetically Informed Interspecies Scaling Framework. Environ. Toxicol. Chem. 2021, 40, 2042–2050. [Google Scholar] [CrossRef]
  112. Bailey, E.; Field, L.; Rawlings, C.; King, R.; Mohareb, F.; Pak, K.-H.; Hughes, D.; Williamson, M.; Ganko, E.; Buer, B.; et al. A Near-Chromosome Level Genome Assembly of the European Hoverfly, Sphaerophoria Rueppellii (Diptera: Syrphidae), Provides Comparative Insights into Insecticide Resistance-Related Gene Family Evolution. BMC Genom. 2022, 23, 198. [Google Scholar] [CrossRef]
  113. Doyle, T.; Jimenez-Guri, E.; Hawkes, W.L.S.; Massy, R.; Mantica, F.; Permanyer, J.; Cozzuto, L.; Hermoso Pulido, T.; Baril, T.; Hayward, A.; et al. Genome-wide Transcriptomic Changes Reveal the Genetic Pathways Involved in Insect Migration. Mol. Ecol. 2022, 31, 4332–4350. [Google Scholar] [CrossRef]
  114. Yuan, H.; Gao, B.; Wu, C.; Zhang, L.; Li, H.; Xiao, Y.; Wu, K. Genome of the Hoverfly Eupeodes Corollae Provides Insights into the Evolution of Predation and Pollination in Insects. BMC Biol. 2022, 20, 157. [Google Scholar] [CrossRef] [PubMed]
  115. Haas, J.; Hayward, A.; Buer, B.; Maiwald, F.; Nebelsiek, B.; Glaubitz, J.; Bass, C.; Nauen, R. Phylogenomic and Functional Characterization of an Evolutionary Conserved Cytochrome P450-Based Insecticide Detoxification Mechanism in Bees. Proc. Natl. Acad. Sci. USA 2022, 119, e2205850119. [Google Scholar] [CrossRef] [PubMed]
  116. Uhl, P.; Franke, L.A.; Rehberg, C.; Wollmann, C.; Stahlschmidt, P.; Jeker, L.; Brühl, C.A. Interspecific Sensitivity of Bees towards Dimethoate and Implications for Environmental Risk Assessment. Sci. Rep. 2016, 6, 34439. [Google Scholar] [CrossRef] [PubMed]
  117. Belanger, S.; Barron, M.; Craig, P.; Dyer, S.; Galay-Burgos, M.; Hamer, M.; Marshall, S.; Posthuma, L.; Raimondo, S.; Whitehouse, P. Future Needs and Recommendations in the Development of Species Sensitivity Distributions: Estimating Toxicity Thresholds for Aquatic Ecological Communities and Assessing Impacts of Chemical Exposures. Integr. Environ. Assess. Manag. 2017, 13, 664–674. [Google Scholar] [CrossRef] [PubMed]
  118. Schmolke, A.; Galic, N.; Feken, M.; Thompson, H.; Sgolastra, F.; Pitts-Singer, T.; Elston, C.; Pamminger, T.; Hinarejos, S. Assessment of the Vulnerability to Pesticide Exposures across Bee Species. Environ. Toxicol. Chem. 2021, 40, 2640–2651. [Google Scholar] [CrossRef]
  119. Topping, C.J.; Aldrich, A.; Berny, P. Overhaul Environmental Risk Assessment for Pesticides. Science 2020, 367, 360–363. [Google Scholar] [CrossRef]
  120. Hätönen, M.; Kantner, C.; Lopez Losada, R.; Ludwig, N.; Benavent González, A.; Riedhammer, C.; Kunz, P.; Panico, S.C.; Laakkonen, E.; Parramon Dolcet, L.; et al. European Arthropods and Their Role in Pollination: Scientific Report of Their Biodiversity, Ecology and Sensitivity to Biocides; European Chemicals Agency: Helsinki, Finland, 2022; ISBN 978-92-9468-131-7. [Google Scholar]
Insects 15 00317 g001
Figure 1. Distribution of sensitivity to imidacloprid calculated from 48 h LD50 values of three fly, eight bee and two beetle species, not accounting (a) and accounting (b) for fresh body weight. The gray area shows the parametric 95% CI (1000 interactions). The red dashed line indicates the LD50 for Apis mellifera. The gray dashed line indicates the HD5 hazardous dose (the dose protecting 95% of the species). LD50 values for bee and beetle species were obtained from the literature [81,82,83,84,85,86,87,88,89]. Values for A. mellifera, Bombus terrestris, Osmia bicornis and Osmia cornifrons are means from different studies [27,33,76,84,89,90,91,92,93,94,95,96] (Table S1).
Figure 1. Distribution of sensitivity to imidacloprid calculated from 48 h LD50 values of three fly, eight bee and two beetle species, not accounting (a) and accounting (b) for fresh body weight. The gray area shows the parametric 95% CI (1000 interactions). The red dashed line indicates the LD50 for Apis mellifera. The gray dashed line indicates the HD5 hazardous dose (the dose protecting 95% of the species). LD50 values for bee and beetle species were obtained from the literature [81,82,83,84,85,86,87,88,89]. Values for A. mellifera, Bombus terrestris, Osmia bicornis and Osmia cornifrons are means from different studies [27,33,76,84,89,90,91,92,93,94,95,96] (Table S1).
Insects 15 00317 g001
Insects 15 00317 g002aInsects 15 00317 g002b
Figure 2. Reproduction comparison between control and treated females at different doses of Confidor® (a.i. imidacloprid) in three dipteran species: Exorista larvarum (a), Sphaerophoria rueppellii (b), and Eristalinus aeneus (c). Bar plots indicate oviposition (the number of mated females that laid eggs), and box plots indicate the fecundity rate per female (E. larvarum number of eggs per female assessed during 1 h; S. rueppellii number of eggs per day per female assessed during 72 h; E. aeneus number of eggs per day per female assessed during 168 h. The line in the center of each box is the median, and the box spans the interquartile range (IQR), containing 50% of the data. Whiskers extend to the minimum and maximum values within 1.5 times the IQR. Outliers are also indicated. Bar plots and box plots with asterisks are significantly different from the control (* p <0.05; ** p < 0.01, *** p < 0.001). Simple size: S. rueppellii (N = 16–28), E. larvarum (N = 6–17 †), and E. aeneus (N = 21–29). † Groups of 5 flies per cage.
Figure 2. Reproduction comparison between control and treated females at different doses of Confidor® (a.i. imidacloprid) in three dipteran species: Exorista larvarum (a), Sphaerophoria rueppellii (b), and Eristalinus aeneus (c). Bar plots indicate oviposition (the number of mated females that laid eggs), and box plots indicate the fecundity rate per female (E. larvarum number of eggs per female assessed during 1 h; S. rueppellii number of eggs per day per female assessed during 72 h; E. aeneus number of eggs per day per female assessed during 168 h. The line in the center of each box is the median, and the box spans the interquartile range (IQR), containing 50% of the data. Whiskers extend to the minimum and maximum values within 1.5 times the IQR. Outliers are also indicated. Bar plots and box plots with asterisks are significantly different from the control (* p <0.05; ** p < 0.01, *** p < 0.001). Simple size: S. rueppellii (N = 16–28), E. larvarum (N = 6–17 †), and E. aeneus (N = 21–29). † Groups of 5 flies per cage.
Insects 15 00317 g002aInsects 15 00317 g002b
Table 1. LD50 values at 48 h (in ng/insect and µg/g of insect body weight) and dose–response model parameters for three dipteran species following acute topical exposure to Confidor® (a.i. imidacloprid). * Values for Eristalinus aeneus were split for the two ascending parts of the dose–effect curve (LD50 could only be calculated considering the first ascending part of the curve. The upper doses were not included in the model).
Table 1. LD50 values at 48 h (in ng/insect and µg/g of insect body weight) and dose–response model parameters for three dipteran species following acute topical exposure to Confidor® (a.i. imidacloprid). * Values for Eristalinus aeneus were split for the two ascending parts of the dose–effect curve (LD50 could only be calculated considering the first ascending part of the curve. The upper doses were not included in the model).
SpeciesnModelSlopeLog-Likelihood Residual Standard Error (df)AICp-Value48 h LD50 95% CI48 h LD50 95% CI
(ng/Insect)(µg/g Insect)
Exorista larvarum233Log-logistic−0.765−14.874.89 (3)37.750.0087467.46302.28–632.6511.667.54–15.79
Sphaerophoria rueppellii205Log-logistic−1.668−19.355.11 (3)46.850.01210.237.81–12.651.350.83–1.86
Eristalinus aeneus *193Log-logistic−1.062−9.862.75 (2)27.740.03618,176.208005.6–28,346.9344.77151.85–537.69
Table 2. BMDL (lower confidence bound of the benchmark dose), BMD (benchmark dose) and BMDU (upper confidence bound of the benchmark dose) that affected oviposition rate (number of females that laid eggs) in ng/insect in three dipteran species. A BMR (benchmark response) of 10% change after acute topical exposure to Confidor® (a.i. imidacloprid) was considered. Best-fit models and their related parameters for the three dipteran species are provided; n = sample size of females that survived >24 h after exposure.
Table 2. BMDL (lower confidence bound of the benchmark dose), BMD (benchmark dose) and BMDU (upper confidence bound of the benchmark dose) that affected oviposition rate (number of females that laid eggs) in ng/insect in three dipteran species. A BMR (benchmark response) of 10% change after acute topical exposure to Confidor® (a.i. imidacloprid) was considered. Best-fit models and their related parameters for the three dipteran species are provided; n = sample size of females that survived >24 h after exposure.
SpeciesnModelBMDLBMDBMDUp-ValueAIC
(ng/Insect)(ng/Insect)(ng/Insect)
Exorista larvarum178Multistage 1°47.0875.11127.120.63933.535
Sphaerophoria rueppellii174Logistic1.702.504.690.649190.797
Eristalinus aeneus209Log-Logistic356.791080.918498.470.178165.298
Table 3. Comparison of lethal (LD50, median lethal dose at 48 h) and sublethal (BMDL, lower confidence bound of the benchmark dose that affects oviposition) doses after acute topical exposure to Confidor® (a.i. imidacloprid) in three dipteran species. SSI = sublethal sensitivity index (LD50/BMDL).
Table 3. Comparison of lethal (LD50, median lethal dose at 48 h) and sublethal (BMDL, lower confidence bound of the benchmark dose that affects oviposition) doses after acute topical exposure to Confidor® (a.i. imidacloprid) in three dipteran species. SSI = sublethal sensitivity index (LD50/BMDL).
SpeciesLD50 BMDLSSI
(ng/Insect)(ng/Insect)
Exorista larvarum467.4647.089.93
Sphaerophoria rueppellii10.231.706.01
Eristalinus aeneus18,176.20356.7950.94
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Henriques Martins, C.A.; Azpiazu, C.; Bosch, J.; Burgio, G.; Dindo, M.L.; Francati, S.; Sommaggio, D.; Sgolastra, F. Different Sensitivity of Flower-Visiting Diptera to a Neonicotinoid Insecticide: Expanding the Base for a Multiple-Species Risk Assessment Approach. Insects 2024, 15, 317. https://doi.org/10.3390/insects15050317

AMA Style

Henriques Martins CA, Azpiazu C, Bosch J, Burgio G, Dindo ML, Francati S, Sommaggio D, Sgolastra F. Different Sensitivity of Flower-Visiting Diptera to a Neonicotinoid Insecticide: Expanding the Base for a Multiple-Species Risk Assessment Approach. Insects. 2024; 15(5):317. https://doi.org/10.3390/insects15050317

Chicago/Turabian Style

Henriques Martins, Cátia Ariana, Celeste Azpiazu, Jordi Bosch, Giovanni Burgio, Maria Luisa Dindo, Santolo Francati, Daniele Sommaggio, and Fabio Sgolastra. 2024. "Different Sensitivity of Flower-Visiting Diptera to a Neonicotinoid Insecticide: Expanding the Base for a Multiple-Species Risk Assessment Approach" Insects 15, no. 5: 317. https://doi.org/10.3390/insects15050317

APA Style

Henriques Martins, C. A., Azpiazu, C., Bosch, J., Burgio, G., Dindo, M. L., Francati, S., Sommaggio, D., & Sgolastra, F. (2024). Different Sensitivity of Flower-Visiting Diptera to a Neonicotinoid Insecticide: Expanding the Base for a Multiple-Species Risk Assessment Approach. Insects, 15(5), 317. https://doi.org/10.3390/insects15050317

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