Sublethal acetamiprid exposure induces immunity, suppresses pathways linked to juvenile hormone synthesis in queens and affects cycle-related signaling in emerging bees

Acetamiprid is the only neonicotinoid registered in the European Union because the risks of neonicotinoids to honey bees and other pollinators are strictly regulated. Herein, we orally exposed honey bee colonies to sublethal concentrations of acetamiprid (20 μg/L) under isolated conditions. After one month of continuous exposure, the emerging bees and queens were collected and analyzed via high-throughput label-free quantitative proteomics using a data-independent acquisition strategy. Six and 34 significantly differentially expressed proteins (DEPs) were identified in the emerging bees and queens, respectively. Mrjp3 was the only DEP found in both sample types/castes, and its opposite regulation illustrated a differential response. The DEPs in the emerging bees (H/ACA RNP, Rap1GAP, Mrjp3, and JHE) suggested that sublethal exposure to acetamiprid affected cell cycle-related signaling, which may affect the life history of workers in the colony. The DEPs with increased levels in queens, such as Mrjps 1-4 and 6-7, hymenoptaecin, and apidaecin 22, indicated an activated immune response. Additionally, the level of farnesyl pyrophosphate synthase (FPPS), which is essential for the mevalonate pathway and juvenile hormone biosynthesis, was significantly decreased in queens. The impaired utilization of juvenile hormone in queens supported the identification of additional DEPs. Furthermore, the proteome changes suggested the existence of increased neonicotinoid detoxification by UDP-glucuronosyltransferase and increased amino acid metabolism. The results suggest that the continuous exposure of bee colonies to acetamiprid at low doses (nanograms per gram in feed) may pose a threat to the colonies. The different exposure routes and durations for the emerging bees and queens in our experiment must be considered, i.e., the emerging bees were exposed as larvae via feeding royal jelly and beebread provided by workers (nurse bees), whereas the queens were fed royal jelly throughout the experiment. The biological consequences of the proteomic changes resulting from sublethal/chronic exposure require future determination.

Systemic pesticides in a solitary bee pollen food store affect larval development and increase pupal mortality

Solitary bees are often exposed to various pesticides applied for pest control on farmland while providing pollination services to food crops. Increasing evidence suggests that sublethal toxicity of agricultural pesticides affects solitary bees differently than the social bees used to determine regulatory thresholds, such as honey bees and bumblebees. Studies on solitary bees are challenging because of the difficulties in obtaining large numbers of eggs or young larvae for bioassays. Here we show the toxic and sublethal developmental effects of four widely used plant systemic pesticides on the Japanese orchard bee (Osmia cornifrons). Pollen food stores of this solitary bee were treated with different concentrations of three insecticides (acetamiprid, flonicamid, and sulfoxaflor) and a fungicide (dodine). Eggs were transplanted to the treated pollen and larvae were allowed to feed on the pollen stores after egg hatch. The effects of chronic ingestion of contaminated pollen were measured until adult eclosion. This year-long study revealed that chronic exposure to all tested pesticides delayed larval development and lowered larval and adult body weights. Additionally, exposure to the systemic fungicide resulted in abnormal larval defecation and increased mortality at the pupal stage, indicating potential risk to bees from fungicide exposure. These findings demonstrate potential threats to solitary bees from systemic insecticides and fungicides and will help in making policy decisions to mitigate these effects.

Environmental Risk Assessment with Energy Budget Models: A Comparison Between Two Models of Different Complexity

The extrapolation of effects from controlled standard laboratory tests to real environmental conditions is a major challenge facing ecological risk assessment (ERA) of chemicals. Toxicokinetic-toxicodynamic (TKTD) models, such as those based on dynamic energy budget (DEB) theory, can play an important role in filling this gap. Through the years, different practical TKTD models have been derived from DEB theory, ranging from the full "standard" DEB animal model to simplified "DEBtox" models. It is currently unclear what impact a different level of model complexity can have on the regulatory risk assessment. In the present study, we compare the performance of two DEB-TKTD models with different levels of complexity, focusing on model calibration on standard test data and on forward predictions for untested time-variable exposure profiles. The first model is based on the standard DEB model with primary parameters, whereas the second is a reduced version with compound parameters, based on DEBkiss. After harmonization of the modeling choices, we demonstrate that these two models can achieve very similar performances both in the calibration step and in the forward prediction step. With the data presented in the present study, selection of the most suitable TKTD model for ERA therefore cannot be based alone on goodness-of-fit or on the precision of model predictions (within current ERA procedures for pesticides) but would likely be based on the trade-off between ease of use and model flexibility. We also stress the importance of modeling choices, such as how to fill gaps in the information content of experimental toxicity data and how to accommodate differences in growth and reproduction between different data sets for the same chemical-species combination. Environ Toxicol Chem 2024;43:440-449. © 2023 ibacon GmbH. Bayer AG and The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Permethrin Contamination of Sawgrass Marshes and Potential Risk for the Imperiled Klot's Skipper Butterfly (Euphyes pilatka klotsi)

Nontarget effects from mosquito control operations are possible in habitats adjacent to areas targeted by ultra-low-volume (ULV) sprays of permethrin for adult mosquito control. We assessed the risks of permethrin exposure to butterflies, particularly the imperiled Klot's skipper, when exposed to ground-based ULV sprays. Samples of larval host plant leaves (sawgrass) were collected in June (in mosquito season) and January (outside mosquito season) of 2015 from sawgrass marsh habitats of the National Key Deer Wildlife Refuge (Big Pine Key, FL, USA) and analyzed for permethrin. Permethrin detection was higher in June (detected on 70% of samples) than in January (30%), and concentrations were significantly higher in June (geomean = 2.1 ng/g, median = 2.4) relative to January (0.4 ng/g, median = 0.2). Dietary risk for 4th to 5th-instar larvae was low based on the measured residues. The AGricultural DISPersal model (Ver. 8.26) was used to estimate permethrin residues on sawgrass following ULV sprays (deposited residues) to estimate immediate postspray risk. Estimated deposited residues (33-543 ng/g) were much higher than measured residues, which leads to a higher risk likelihood for butterfly larvae immediately after ULV sprays. The difference between estimated and measured residues, and between the two risk estimations, reflects uncertainty in risk estimates based on the measured residues. Research on modeling deposited pesticide residues following ground-based ULV spray is limited. More research on estimating deposited pesticide residues from truck-mounted ULV sprayers could help reduce uncertainty in the risk predictions for nontarget insects like butterflies. Environ Toxicol Chem 2024;43:267-278. Published 2023. This article is a U.S. Government work and is in the public domain in the USA. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Trisiloxane Surfactants Negatively Affect Reproductive Behaviors and Enhance Viral Replication in Honey Bees

Trisiloxane surfactants are often applied in formulated adjuvant products to blooming crops, including almonds, exposing the managed honey bees (Apis mellifera) used for pollination of these crops and persisting in colony matrices, such as bee bread. Despite this, little is known regarding the effects of trisiloxane surfactants on important aspects of colony health, such as reproduction. In the present study, we use laboratory assays to examine how exposure to field-relevant concentrations of three trisiloxane surfactants found in commonly used adjuvant formulations affect queen oviposition rates, worker interactions with the queen, and worker susceptibility to endogenous viral pathogens. Trisiloxane surfactants were administered at 5 mg/kg in pollen supplement diet for 14 days. No effects on worker behavior or physiology could be detected, but our results demonstrate that hydroxy-capped trisiloxane surfactants can negatively affect queen oviposition and methyl-capped trisiloxane surfactants cause increased replication of Deformed Wing Virus in workers, suggesting that trisiloxane surfactant use while honey bees are foraging may negatively impact colony longevity and growth. Environ Toxicol Chem 2024;43:222-233. © 2023 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Further Limitations of Synthetic Fungicide Use and Expansion of Organic Agriculture in Europe Will Increase the Environmental and Health Risks of Chemical Crop Protection Caused by Copper-Containing Fungicides

Copper-containing fungicides have been used in agriculture since 1885. The divalent copper ion is a nonbiodegradable multisite inhibitor that has a strictly protective, nonsystemic effect on plants. Copper-containing plant protection products currently approved in Germany contain copper oxychloride, copper hydroxide, and tribasic copper sulfate. Copper is primarily used to control oomycete pathogens in grapevine, hop, potato, and fungal diseases in fruit production. In the environment, copper is highly persistent and toxic to nontarget organisms. The latter applies for terrestric and aquatic organisms such as earthworms, insects, birds, fish, Daphnia, and algae. Hence, copper fungicides are currently classified in the European Union as candidates for substitution. Pertinently, copper also exhibits significant mammalian toxicity (median lethal dose oral = 300-2500 mg/kg body wt in rats). To date, organic production still profoundly relies on the use of copper fungicides. Attempts to reduce doses of copper applications and the search for copper substitutes have not been successful. Copper compounds compared with modern synthetic fungicides with similar areas of use display significantly higher risks for honey bees (3- to 20-fold), beneficial insects (6- to 2000-fold), birds (2- to 13-fold), and mammals (up to 17-fold). These data contradict current views that crop protection in organic farming is associated with lower environmental or health risks. Further limitations in the range and use of modern single-site fungicides may force conventional production to fill the gaps with copper fungicides to counteract fungicide resistance. In contrast to the European Union Green Deal goals, the intended expansion of organic farming in Europe would further enhance the use of copper fungicides and hence increase the overall risks of chemical crop protection in Europe. Environ Toxicol Chem 2024;43:19-30. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Pesticides bioassays using neotropical aquatic species: Trends during the last twenty years and future challenges in Argentina

The presence of pesticides in aquatic ecosystems is one of the most relevant stressors which biota usually face. Laboratory tests using model organisms for pesticides toxicity assessment are employed worldwide. The use of these species has been encouraged in the scientific community due to their advantageous features and their acceptation by regulatory and standardization organizations. However, non-model species as well as those belonging particular ecosystems could contribute in the laboratory-field toxicity extrapolation. In this context, this work aims on exploring the state of the ecotoxicological studies of pesticides in neotropical aquatic species, focusing on bioassays performed in Argentina over the last 20 years as a case of study. Furthermore, we analyzed the possible advantages and disadvantages of these studies, possible differential sensitivities among native and model species, and future challenges to be faced. The analysis of more than 150 publications allowed identify the chemical identity of tested compounds, organisms used for the bioassays, characteristics of the experimental designs, and the toxicity endpoints. Particularly, the studied cases showed that the tested chemicals are related to those most used in the agricultural activity in Argentina, the predilection for particular species in some taxonomic groups (e.g. amphibians), and the wide election of biochemical biomarkers in the studies. Regarding the sensitivity comparison between native and non-native species, the amount of data available indicates that there is not a clear difference beyond some particular cases. However, deeper understanding of toxic effects of pesticides on non-model species could help in a more comprehensive ecological risk assessment in different ecosystems.

Real-time monitoring of honeybee colony daily activity and bee loss rates can highlight the risk posed by a pesticide

Information on honeybee foraging performance and especially bee loss rates at the colony level are crucial for evaluating the magnitude of effects due to pesticide exposure, thereby ensuring that protection goals for honeybee colonies are met (i.e. threshold of acceptable effects). However, current methods for monitoring honeybee foraging activity and mortality are very approximate (visual records) or are time-limited and mostly based on single cohort analysis. We therefore assess the potential of bee counters, that enable a colony-level and continuous monitoring of bee flight activity and mortality, in pesticide risk assessment. After assessing the background activity and bee loss rates, we exposed colonies to two concentrations of sulfoxaflor (a neurotoxic insecticide) in sugar syrup: a concentration that was considered to be field realistic (0.59 μg/ml) and a higher concentration (2.36 μg/ml) representing a worst-case exposure scenario. We did not find any effect of the field-realistic concentration on flight activity and bee loss rates. However, a two-fold decrease in daily flight activity and a 10-fold increase in daily bee losses were detected in colonies exposed to the highest sulfoxaflor concentration as compared to before exposure. When compared to the theoretical trigger values associated with the specific protection goal of 7 % colony-size reduction, the observed fold changes in daily bee losses were often found to be at risk for colonies. In conclusion, the real-time and colony-level monitoring of bee loss rates, combined with threshold values indicating at which levels bee loss rates threaten the colony, have great potential for improving regulatory pesticide risk assessments for honeybees under field conditions.

Field Availability and Avoidance of Imidacloprid-Treated Soybean Seeds and Cotyledons by Birds

Treated seeds and their cotyledons can present a toxicological risk to seed-eating birds. To assess whether avoidance behavior limits exposure and consequently the risk to birds, three fields were sown with soybeans. Half of the surface of each field was sown with seeds treated with 42 g/100 kg seed of insecticide imidacloprid (T plot, treated) and the other half with seeds without imidacloprid (C plot, control). Unburied seeds were surveyed in C and T plots at 12 and 48 h post-sowing. Damaged seedlings were surveyed in C and T plots at 12 days post-sowing. The abundance and richness of birds was surveyed at the field level (without distinguishing between C and T plots) before, during, and after sowing, and 12 days post-sowing. Unburied seed density was higher in the headlands of the T plots than in the C plots, but did not differ between 12 and 48 h. The damage to cotyledons of seedlings was 15.4% higher in C plots than in T plots. The abundance and richness/ha of birds that eat seeds and cotyledons were lower after sowing, indicating a deterrent effect on birds by sowing imidacloprid-treated seeds. Although the variation in seed density over time does not allow solid conclusions to be drawn about the avoidance of seeds treated by birds, the seedling results suggest an aversive effect of imidacloprid-treated soybeans on birds. The dominant species was the eared dove (Zenaida auriculata), whose risk of acute poisoning by imidacloprid in soybean seeds and cotyledons was low, according to its toxicity exposure ratio, foraged area of concern, and foraged time of concern. Environ Toxicol Chem 2023;42:1049-1060. © 2023 SETAC.

Nematode Community of a Natural Grassland Responds Sensitively to the Broad-Spectrum Fungicide Mancozeb in Soil Microcosms

Fungicides make up the largest part of total pesticide use, with the dithiocarbamate mancozeb being widely applied as a nonsystemic contact pesticide to protect a wide range of field crops against fungal diseases. Although nematodes are key drivers of soil functioning, data on effects of fungicides, and especially mancozeb, on these nontarget organisms are scarce. Therefore, the effects of mancozeb on a soil nematode community from a natural grassland were assessed in small-scale soil microcosms. Nematodes were exposed to mancozeb-spiked soil in six nominal concentrations (7-133 mg/kg dry soil) and analyzed after 14, 56, and 84 days in terms of densities, genus composition, and functional traits. Because this fungicide is known to quickly degrade in soils (50% degradation time <1 day), mancozeb concentrations were analyzed for all sampling occasions. Chemical analysis revealed considerably lower measured concentrations compared with the aimed nominal soil concentrations at the beginning of the exposure (1-18 mg/kg dry soil), suggesting fast degradation during the spiking process. Nevertheless, the native nematode community responded sensitively to the fungicide mancozeb, revealing lower no-observed-effect concentration and 10% effect concentration (EC10) values than reported for other soil invertebrates such as springtails and earthworms. Using the EC10 for the most sensitive nematode community endpoint (percentage of predators and omnivores: 1.2 mg/kg dry soil), the risk assessment exhibited a toxicity exposure ratio of 0.66 and, thus, a high risk of mancozeb for soil nematodes. Keeping in mind their abundance and their central roles in soil food-web functioning, the demonstrated sensitivity to a widely applied fungicide underscores the relevance of the inclusion of nematodes into routine risk-assessment programs for pesticides. Environ Toxicol Chem 2022;41:2420-2430. © 2022 SETAC.

Honey bee toxicological responses do not accurately predict environmental risk of imidacloprid to a solitary ground-nesting bee species

Honey bees (Apis mellifera) are the current model species for pesticide risk assessments, but considering bee diversity, their life histories, and paucity of non-eusocial bee data, this approach could underestimate risk. We assessed whether honey bees were an adequate risk predictor to non-targets. We conducted oral and contact bioassays for Leioproctus paahaumaa, a solitary ground-nesting bee, and A. mellifera, using imidacloprid (neonicotinoid) and dimethoate (organophosphate). The bees responded inconsistently; L. paahaumaa were 36 and 194 times more susceptible to oral and topically applied imidacloprid than A. mellifera, but showed comparable sensitivity to dimethoate. Furthermore, the proposed safety factor of ten applied to honey bee endpoints did not cover the interspecific sensitivity difference. Our standard-setting study highlights the urgent need for more comparative inter-species toxicity studies and the development of standardized toxicity protocols to ensure regulatory pesticide risk assessment frameworks are protective of diverse pollinators.

Chronic and Acute Effects of Imidacloprid on a Simulated BEEHAVE Honeybee Colony

Honeybees (Apis mellifera) are important pollinators for wild plants as well as for crops, but honeybee performance is threatened by several stressors including varroa mites, gaps in foraging supply, and pesticides. The consequences of bee colony longtime exposure to multiple stressors are not well understood. The vast number of possible stressor combinations and necessary study duration require research comprising field, laboratory, and simulation experiments. We simulated long-term exposure of a honeybee colony to the insecticide imidacloprid and to varroa mites carrying the deformed wing virus in landscapes with different temporal gaps in resource availability as single stressors and in combinations. Furthermore, we put a strong emphasis on chronic lethal, acute sublethal, and acute lethal effects of imidacloprid on honeybees. We have chosen conservative published values to parameterize our model (e.g., highest reported imidacloprid contamination). As expected, combinations of stressors had a stronger negative effect on bee performance than each single stressor alone, and effect sizes were larger after 3 years of exposure than after the first year. Imidacloprid-caused reduction in bee performance was almost exclusively due to chronic lethal effects because the thresholds for acute effects were rarely met in simulations. In addition, honeybee colony extinctions were observed by the last day of the first year but more pronounced on the last days of the second and third simulation year. In conclusion, our study highlights the need for more long-term studies on chronic lethal effects of pesticides on honeybees. Environ Toxicol Chem 2022;41:2318-2327. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

A systematic review of research conducted by pioneer groups in ecotoxicological studies with bees in Brazil: advances and perspectives

Brazil presents the most threatened endemic or rare species among neotropical regions, with the Hymenoptera order, to which bees belong, classified as a high-risk category. In Brazil, the main cause of bee death is the indiscriminate use of pesticides. In this context, groups such as Bee Ecotoxicology and Conservation Laboratory (LECA in Portuguese) and Bees and Environmental Services (ASAs in Portuguese) have become a reference in studies evaluating the impacts of pesticides on bees since 1976. Thus, the objective of this review was to conduct a quantitative and qualitative review of the studies conducted by these groups to evaluate and compile the advances made over the years, identify potential knowledge gaps for future studies, and support the sensitivities of stingless bees when compared to the species Apis mellifera. The quantitative analyses showed that most studies were carried out in the genus Apis, under laboratory conditions. However, more recently (since 2003), studies have also focused on stingless bees and the neonicotinoid class of insecticides. The most relevant gaps identified were the lack of studies under field conditions and on bee biology. The qualitative analyses indicated that Brazilian stingless bees are more susceptible to pesticides than A. mellifera and require a much lower average dose, concentration, or lethal time to display morphological and behavioral damage or decreased lifespan. Thus, future studies should work towards establishing more representative protocols for stingless bees. Furthermore, public policies must be created for the protection and conservation of bees native to Brazil.

Mix-Tool: An Edge-of-Field Approach to Predict Pesticide Mixtures of Concern in Surface Water From Agricultural Crops

Current regulation on the authorization of plant protection products (PPPs) in the European Union is limited to the evaluation of ecological risks for the single active substances they contain. However, plant protection treatments in agriculture often consist of PPPs already containing more than one active substance; moreover, each cropped field receives multiple applications per year, leading to complex pesticide mixtures in the environment. Different transport processes lead to a multitude of heterogeneous and potentially toxic substances that, for example, may reach water bodies and act simultaneously on natural freshwater ecosystems. In this context, the development of methodologies and tools to manage risks of pesticides mixtures is imperative to improve the current ecological risk assessment procedures and to avoid further deterioration of ecological quality of natural resources. The present study suggests new procedures for identifying pesticide mixtures of potential concern released from agricultural crops in surface water. The approach follows the European Union regulatory context for the authorization of PPPs in the market (edge-of field risk assessment) and requires the use of Forum for the Co-ordination of pesticide fate models and their Use (FOCUS) models (Step 3 and 4) for calculating the concentrations in surface water of mixture components on a daily basis. Moreover, it uses concentration addition models to calculate the toxic potency of the pesticide mixtures released by a treated crop. To implement this procedure, we developed a simple Microsoft-Excel-based tool. We also considered two case studies (maize and apple tree), representative of Italian agricultural scenarios for annual and perennial crops. Moreover, we compared results with 3 years of monitoring data of surface water bodies of the Lombardia region (northern Italy) where the two crops are largely present. Environ Toxicol Chem 2022;41:2028-2038. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Sensitivity of the stripe-faced dunnart, Sminthopsis macroura (Gould 1845), to the insecticide, fipronil; implications for pesticide risk assessments in Australia

A lack of toxicity data quantifying responses of Australian native mammals to agricultural pesticides prompted an investigation into the sensitivity of the stripe-faced dunnart, Sminthopsis macroura (Gould 1845) to the insecticide, fipronil (5-amino-3-cyano-1-(2,6-dichloro-4-trifluoromethylphenyl)-4-trifluoromethylsulfinyl pyrazole, CAS No. 120068-37-3). Using the Up-And-Down method for determining acute oral toxicity in mammals (OECD) median lethal dose estimates of 990 mg kg^-1 (95% confidence interval (CI) = 580.7-4770.0 mg kg^-1) and 270.4 mg kg^-1 (95% CI = 0.0->20,000.0 mg kg^-1) were resolved for male and female S. macroura, respectively. The difference between median lethal dose estimates for males and females may have been influenced by the older ages of two female dunnarts. Consequently, further modelling of female responses to fipronil doses used the following assumptions: (a) death at 2000 mg kg^-1, (b) survival at 500 mg kg^-1 and (c) a differential response (both survival and death) at 990 mg kg^-1. This modelling revealed median lethal dose estimates for female S. macroura of 669.1 mg kg^-1 (95% CI = 550-990 mg kg^-1; assuming death at 990 mg kg^-1) and 990 mg kg^-1 (95% CI = 544.7-1470 mg kg^-1; assuming survival at 990 mg kg^-1). These median lethal dose estimates are 3-10-fold higher than available LD₅₀ values of 94 mg kg^-1 for a similarly sized eutherian mammal, Mus musculus (L. 1758) and 97 mg kg^-1 for Rattus norvegicus (Birkenhout 1769). Implications for pesticide risk assessments in Australia are discussed.

TWAc-Check: A New Approach to Determine the Appropriate Use of Time-Weighted Average Concentration in Aquatic Risk Assessment

In pesticide risk assessment, regulatory acceptable concentrations for surface water bodies (RACsw,ch) are used that are derived from standard studies with continuous exposure of organisms to a test compound for days or months. These RACsw,ch are compared with the maximum tested concentration of more realistic exposure scenarios. However, the actual exposure duration could be notably shorter (e.g., hours) than the standard study, which intentionally leads to an overly conservative Tier 1 risk assessment. This discrepancy can be addressed in a risk assessment using the time-weighted average concentration (TWAc). In Europe, the applicability of TWAc for a particular risk assessment is evaluated using a complex decision scheme, which has been controversial; thus we propose an alternative approach: We used TWAc-check (which is based on the idea that the TWAc concept is just a model for aquatic risk assessment) to test whether the use of a TWAc is appropriate for such assessment. The TWAc-check method works by using predicted-measured diagrams to test how well the TWAc model predicts experimental data from peak exposure experiments. Overestimated effects are accepted because the conservatism of the TWAc model is prioritized over the goodness of fit. We illustrate the applicability of TWAc-check by applying it to various data sets for different species and substances. We demonstrate that the applicability is case dependent. Specifically, TWAc-check correctly identifies that the use of TWAc is not appropriate for early onset of effects or delayed effects. The proposed concept shows that the time window is a decisive factor as to whether or not the model is acceptable and that this concept can be used as a potential refinement option prior to the use of toxicokinetic-toxicodynamic models. Environ Toxicol Chem 2022;41:1778-1787. © 2022 Bayer AG. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Probabilistic co-occurrence assessment for suites of listed species

Section 7 of the Endangered Species Act requires the US Environmental Protection Agency (US EPA) to consult with the Services (US Fish and Wildlife Service and National Marine Fisheries Service) over potential pesticide impacts on federally listed species. Consultation is complicated by the large number of pesticide products and listed species, as well as by lack of consensus on best practices for conducting co-occurrence analyses. Previous work demonstrates that probabilistic estimates of species' ranges and pesticide use patterns improve these analyses. Here we demonstrate that such estimates can be made for suites of sympatric listed species. Focusing on two watersheds, one in Iowa and the other in Mississippi, we obtained distribution records for 13 species of terrestrial and aquatic listed plants and animals occurring therein. We used maximum entropy modeling and bioclimatic, topographic, hydrographic, and land cover variables to predict species' ranges at high spatial resolution. We constructed probabilistic spatial models of use areas for two pesticides based on the US Department of Agriculture Cropland Data Layer and reduced classification errors by incorporating information on the relationships between individual pixels and their neighbors using object-based images analysis. We then combined species distribution and crop footprint models to derive overall probability of co-occurrence of listed species and pesticide use. For aquatic species, we also integrated an estimate of downstream residue transport. We report each separate species-by-use-area co-occurrence estimate and also combine these modeled co-occurrence probabilities across species within watersheds to produce an overall metric of potential pesticide exposure risk for these listed species at the watershed level. We propose that the consultation process between US EPA and the Services be based on such batched estimation of probabilistic co-occurrence for multiple listed species at a regional scale. Integr Environ Assess Manag 2022;18:1088-1100. © 2021 SETAC.

Acute oral toxicity and risks of four classes of systemic insecticide to the Common Eastern Bumblebee (Bombus impatiens)

The Common Eastern Bumblebee (Bombus impatiens) is native to North America with an expanding range across Eastern Canada and the USA. This species is commercially produced primarily for greenhouse crop pollination and is a common and abundant component of the wild bumblebee fauna in agricultural, suburban and urban landscapes. However, there is a dearth of pesticide toxicity information about North American bumblebees. The present study determined the acute oral lethal toxicity (48-h LD50) of: the butenolide, flupyradifurone (>1.7 μg/bee); the diamide, cyantraniliprole (>0.54 μg/bee); the neonicotinoid, thiamethoxam (0.0012 μg/bee); and the sulfoximine, sulfoxaflor (0.0177 μg/bee). Compared with published honey bee (Apis mellifera) LD50 values, the present study shows that sulfoxaflor and thiamethoxam are 8.3× and 3.3× more acutely toxic to B. impatiens, whereas flupyradifurone is more acutely toxic to A. mellifera. The current rule of thumb for toxicity extrapolation beyond the honey bee as a model species, termed 10× safety factor, may be sufficient for bumblebee acute oral toxicity. A comparison of five risk assessment equations suggested that the Standard Risk Approach (SRA) and Fixed Dose Risk Approach (FDRA) provide more nuanced levels of risk evaluation compared to the Exposure Toxicity Ratio (ETR), Hazard Quotient (HQ), and Risk Quotient (RQ), primarily because the SRA and FDRA take into account real world variability in pollen and nectar pesticide residues and the chances that bees may be exposed to them.

Delayed effects of a single dose of a neurotoxic pesticide (sulfoxaflor) on honeybee foraging activity

Pesticide risk-assessment guidelines for honeybees (Apis mellifera) generally require determining the acute toxicity of a chemical over the short-term through fix-duration tests. However, potential long-lasting or delayed effects resulting from an acute exposure (e.g. a single dose) are often overlooked, although the modification of a developmental process may have life-long consequences. To investigate this question, we exposed young honeybee workers to a single sublethal field-realistic dose of a neurotoxic pesticide, sulfoxaflor, at one of two amounts (16 or 60 ng), at the moment when they initiated orientation flights (preceding foraging activity). We then tracked in the field their flight activity and lifespan with automated life-long monitoring devices. Both amounts of sulfoxaflor administered reduced the total number of flights but did not affect bee survival and flight duration. When looking at the time series of flight activity, effects were not immediate but delayed until foraging activity with a decrease in the daily number of foraging flights and consequently in their total number (24 and 33% less for the 16 and 60 ng doses, respectively). The results of our study therefore blur the general assumption in honeybee toxicology that acute exposure results in immediate and rapid effects and call for long-term recording and/or time-to-effect measurements, even upon exposure to a single dose of pesticide.

Bird exposure to fungicides through the consumption of treated seeds: A study of wild red-legged partridges in central Spain

Sown seeds are a key component of many farmland birds' diets due to natural food shortages in autumn and winter. Because these seeds are often treated with pesticides, their ingestion by birds can result in toxic effects. For risk assessment, data on treated seed toxicity should be combined with information about exposure risk for wild birds and the factors that modulate it. We characterized the exposure of red-legged partridges to pesticide-treated seeds through the analysis of digestive contents of birds shot by hunters (n = 194) in an agricultural region in central Spain. We measured the contribution of sown seeds to the partridges' diet and how it related to pesticide exposure. Moreover, we evaluated the influence of landscape composition on the intake of sown seeds and pesticides by partridges. During peak sowing time, seeds constituted half (50.7%) of the fresh biomass ingested by partridges, which consumed mostly winter cereal seeds (42.3% of biomass). Residues of seven fungicides and one insecticide (active ingredients) were detected in 33.0% of birds. The presence of pesticides in digestive contents was linked to the ingestion of cereal sown seeds. Moreover, dietary exposure of birds to pesticides was modulated by landscape characteristics, being lower in areas with heterogeneous landscapes, greater habitat mosaic and more natural vegetation. The estimated dietary intake of pesticides resulting from our field observations, in combination with experimental data on pesticide toxicity, raise concerns about the risks that pesticide-treated cereal seeds pose to granivorous bird populations. Our results highlight the importance of farming landscape composition and diversification, which should be considered as a priority in the agricultural policy to mitigate pesticide risks to farmland birds through the consumption of treated seeds.

The Role of Source-Sink Dynamics in the Assessment of Risk to Nontarget Arthropods from the Use of Plant Protection Products

The concept of source-sink dynamics as a potentially important component of metapopulation dynamics was introduced in the 1980s. The objective of the present review was to review the considerable body of work that has been developed, to consider its theoretical implications as well as to understand how source-sink dynamics may manifest under field conditions in the specific case of nontarget arthropods in the agricultural environment. Our review concludes that metapopulation dynamics based on field observations are often far more complex than existing theoretical source-sink models would indicate, because they are dependent on numerous population processes and influencing factors. The difficulty in identifying and measuring these factors likely explains why empirical studies assessing source-sink dynamics are scarce. Furthermore, we highlight the importance of considering the spatial and temporal heterogeneity of agricultural landscapes when assessing the population dynamics of nontarget arthropods in the context of the risk from the use of plant protection products. A need is identified to further develop and thoroughly validate predictive population models, which can incorporate all factors relevant to a specific system. Once reliable predictive models for a number of representative nontarget arthropod species are available, they could provide a meaningful tool for refined risk evaluations (higher tier level risk assessment), addressing specific concerns identified at the initial evaluation stages (lower tier level risk assessment). Environ Toxicol Chem 2021;40:2667-2679. © 2021 ERM, FMC, Syngenta, Bayer AG, BASF SE, Corteva agriscience. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Applying a Hybrid Modeling Approach to Evaluate Potential Pesticide Effects and Mitigation Effectiveness for an Endangered Fish in Simulated Oxbow Habitats

The occurrence of some species listed under the United States' Endangered Species Act in agricultural landscapes suggests that their habitats could potentially be exposed to pesticides. However, the potential effects from such exposures on populations are difficult to estimate. Mechanistic models can provide an avenue to estimating the potential impacts on populations, considering realistic assumptions about the ecology of the species, the ecosystem it is part of, and the potential exposures within the habitat. In the present study, we applied a hybrid model of the Topeka shiner (Notropis topeka), a small endangered cyprinid fish endemic to the US Midwest, to assess the potential population-level effects of realistic exposures to a fungicide (benzovindiflupyr). The Topeka shiner populations were simulated in the context of the food web found in oxbow habitats that are the focus of ongoing habitat restoration efforts for the species. We applied realistic, time-variable exposure scenarios and represented lethal and sublethal effects to individual Topeka shiners using toxicokinetic-toxicodynamic models. With fish in general showing the highest sensitivity to the compound, direct effects on simulated Topeka shiner populations governed the population-level effects. We characterized the population-level effects of different exposure scenarios with exposure multiplication factors (EMFs) applied. The introduction of a vegetative filter strip (VFS; 15 ft; 4.6 m) between the treated area and the oxbow habitat was shown to be effective as mitigation because EMFs were 2 to 3 times higher than for the exposure scenario without VFS. Environ Toxicol Chem 2021;40:2615-2628. © 2021 SETAC.

Assessment of the Vulnerability to Pesticide Exposures Across Bee Species

In many countries, the western honey bee is used as surrogate in pesticide risk assessments for bees. However, uncertainty remains in the estimation of pesticide risk to non-Apis bees because their potential routes of exposure to pesticides, life histories, and ecologies differ from those of honey bees. We applied the vulnerability concept in pesticide risk assessment to 10 bee species including the honey bee, 2 bumble bee species, and 7 solitary bee species with different nesting strategies. Trait-based vulnerability considers the evaluation of a species at the level of both the organism (exposure and effect) and the population (recovery), which goes beyond the sensitivity of individuals to a toxicant assessed in standard laboratory toxicity studies by including effects on populations in the field. Based on expert judgment, each trait was classified by its relationship to the vulnerability to pesticide exposure, effects (intrinsic sensitivity), and population recovery. The results suggested that the non-Apis bees included in our approach are potentially more vulnerable to pesticides than the honey bee due to traits governing exposure and population recovery potential. Our analysis highlights many uncertainties related to the interaction between bee ecology and the potential exposures and population-level effects of pesticides, emphasizing the need for more research to identify suitable surrogate species for higher tier bee risk assessments. Environ Toxicol Chem 2021;40:2640-2651. © 2021 SETAC.

Human Exposures to Neonicotinoids in Kumasi, Ghana

Neonicotinoid insecticides (NNIs) are now popular in many agricultural systems across Africa; however, the extent of human exposures to NNIs in African countries is scarcely reported. The present study evaluates neonicotinoid exposures in the consumer population of Kumasi, a cosmopolitan city in Ghana. A total of 75 human urine samples were collected from healthy volunteers (nonfarmers, aged 13-80 yr) and analyzed with a liquid chromatography electrospray ionization tandem mass spectrometry system. Seven NNIs and 3 NNI metabolites were detected in the following pattern (frequency, median concentration, maximum concentration): N-dm-acetamiprid (94.7%, 0.41 µg/L, 8.79 µg/L) > imidacloprid (70.7%, 0.15 µg/L, 211.62 µg/L) > N-(6-chloro-3-pyridylmethyl)-N-ethyl-N'-methylformamidine (62.2%, 0.43 µg/L, 53.85 µg/L) > 2-[N-(6-chloro-3-pyridylmethyl)-N-ethylamino]-2-(methylimino)acetic acid (56.8%, 0.10 µg/L, 3.53 µg/L) > clothianidin (40%, >limit of quantification [LOQ], 0.45 µg/L) > nitenpyram (18.7%, >LOQ, 0.14 µg/L) ≈ thiamethoxam (18.7%, >LOQ, 0.21 µg/L) > dinotefuran (12.0%, >LOQ, 1.01 µg/L) > acetamiprid (2.7%, >LOQ, 0.08 µg/L) ≈ thiacloprid (2.7%, >LOQ, 0.14 µg/L). Approximately 92% of the subjects were found to be exposed to multiple neonicotinoids simultaneously. The mean, median, and maximum imidacloprid equivalent of the relative potency factor of NNIs were found to be 1.6, 0.5, and 22.52, respectively. The median estimated daily intakes of acetamiprid, imidacloprid, and nitenpyram were 0.47, 1.27, and 0.02 µg/kg/d for females and 0.91, 0.66, and 0.08 µg/kg/d for males, respectively. The maximum daily intakes of all the NNIs were <1% of their chronic reference doses (cRfDs), except for imidacloprid and thiacloprid which recorded maximum daily intakes corresponding to 17.97 and 8.28% of cRfDs, respectively. To the best of our knowledge, the present study is the first biomonitoring report on neonicotinoid insecticides in Africa. Environ Toxicol Chem 2021;40:2306-2318. © 2021 SETAC.

Effects of glyphosate spray-drift on plant flowering

Recent studies have shown that sub-lethal doses of herbicides may affect plant flowering, however, no study has established a direct relationship between the concentrations of deposited herbicide and plant flowering. Here the aim was to investigate the relationship between herbicide spray drift deposited on non-target plants and plant flowering in a realistic agro-ecosystem setting. The concentrations of the herbicide glyphosate deposited on plants were estimated by measuring the concentration of a dye tracer applied together with the herbicide. The estimated maximal and average deposition of glyphosate within the experimental area corresponded to 30 g glyphosate/ha (2.08% of the label rate of 1440 g a.i./ha) and 2.4 g glyphosate/ha (0.15% label rate), respectively, and the concentrations decreased rapidly with increasing distance from the spraying track. However, there were not a unique relation between distance and deposition, which indicate that heterogeneities of turbulence, wind speed and/or direction can strongly influence the deposition from 1 min to another during spraying. The effects of glyphosate on cumulative flower numbers and flowering time were modelled using Gompertz growth models on four non-target species. Glyphosate had a significantly negative effect on the cumulative number of flowers on Trifolium pratense and Lotus corniculatus, whereas there were no significant effects on Trifolium repens, and a positive, but non-significant, effect on number of flowers on Cichorium intybus. Glyphosate did not affect the flowering time of any of the four species significantly. Lack of floral resources is known to be of major importance for pollinator declines. The implications of the presented results for pesticide risk assessment are discussed.

Larval oral exposure to thiacloprid: Dose-response toxicity testing in solitary bees, Osmia spp. (Hymenoptera: Megachilidae)

Risk assessment of pesticides involves ecotoxicological testing. In case pesticide exposure to bees is likely, toxicity tests are performed with honey bees (Apis mellifera), with a tiered approach, for which validated and internationally accepted test protocols exist. However, concerns have grown regarding the protection of non-Apis bees [bumble bees (Bombus spp.), solitary and stingless bees], given their different life cycles and therefore distinct exposure routes. Larvae of solitary bees of the genus Osmia feed on unprocessed pollen during development, yet no toxicity test protocol is internationally accepted or validated to assess the impact of pesticide exposure during this stage of their life cycle. Therefore, the purpose of this study is to further validate a test protocol with two solitary bee species (O. cornuta and O. bicornis) to assess lethal and sublethal effects of pesticide exposure on larval development. Larvae were exposed to thiacloprid (neonicotinoid insecticide) mixed in a new, artificial pollen provision. Both lethal (developmental and winter mortality) and sublethal endpoints (larval development time, pollen provision consumption, cocoon weight, emergence time and adult longevity) were recorded. Effects of lower, more environmentally realistic doses were only reflected in sublethal endpoints. In both bee species, thiacloprid treatment was associated with increased developmental mortality and larval development time, and decreased pollen provision consumption and cocoon weight. The test protocol proved valid and robust and showed that for higher doses of thiacloprid the acute endpoint (larval mortality) is sufficient. In addition, new insights needed to develop a standardized test protocol were acquired, such as testing of a positive control for the first time and selection of male and female individuals at egg level.

Contact toxicity of three insecticides for use in tier I pesticide risk assessments with Megachile rotundata (Hymenoptera: Megachilidae)

The current pesticide risk assessment paradigm may not adequately protect solitary bees as it focuses primarily on the honey bee (Apis mellifera). The alfalfa leafcutting bee (Megachile rotundata) is a potential surrogate species for use in pesticide risk assessment for solitary bees in North America. However, the toxicity of potential toxic reference standards to M. rotundata will need to be determined before pesticide risk assessment tests (tier I trials) can be implemented. Therefore, we assessed the acute topical toxicity and generated LD50 values for three insecticides: dimethoate (62.08 ng a.i./bee), permethrin (50.01 ng a.i./bee), and imidacloprid (12.82 ng a.i/bee). The variation in the mass of individual bees had a significant but small effect on these toxicity estimates. Overall, the toxicity of these insecticides to M. rotundata were within the 10-fold safety factor currently used with A. mellifera toxicity estimates from tier I trials to estimate risk to other bee species. Therefore, tier I pesticide risk assessments with solitary bees may not be necessary, and efforts could be directed to developing more realistic, higher-tier pesticide risk assessment trials for solitary bees.

Application of General Unified Threshold Models of Survival Models for Regulatory Aquatic Pesticide Risk Assessment Illustrated with an Example for the Insecticide Chlorpyrifos

Mathematical models within the General Unified Threshold models of Survival (GUTS) framework translate time-variable chemical exposure information into expected survival of animals. The GUTS models are species and compound specific and explicitly describe the internal exposure dynamics in an organism (toxicokinetics) and the related damage and effect dynamics (toxicodynamics), thereby connecting the external exposure concentration dynamics with the simulated mortality or immobility over time. In a recent scientific opinion on toxicokinetic-toxicodynamic (TKTD) models published by the European Food Safety Authority (EFSA), the GUTS modeling framework was considered ready for use in the aquatic risk assessment for pesticides and aquatic fauna. The GUTS models are suggested for use in risk assessment, if they are sufficiently validated for a specific substance-species combination. This paper aims to illustrate how they can be used in the regulatory environmental risk assessment for pesticides for a specific type of refinement, that is, when risks are triggered by lower tiers in acute as well as in chronic risk assessment and mortality or immobility is the critical endpoint. This approach involves the evaluation of time-variable exposure regimes in a so-called "Tier-2C" assessment. The insecticide chlorpyrifos was selected as an example compound because a large data set was available. The GUTS models for 13 different freshwater arthropods and 8 different theoretical aquatic exposure profiles were used to calculate a series of GUTS-based risk estimates, including exposure profile-specific multiplication factors leading to 50% mortality or immobility at the end of the tested profile (LP50/EP50) as "margins of safety." To put the use of GUTS models within the tiered aquatic risk assessment into perspective, GUTS models for the 13 aquatic arthropods were also used to predict the environmental risks of a measured chlorpyrifos exposure profile from an experimental ditch study (Tier-3 approach), and the results are discussed in the context of calibration of the tiered approach. Integr Environ Assess Manag 2021;17:243-258. © 2020 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Multi-level analysis of exposure to triazole fungicides through treated seed ingestion in the red-legged partridge

Triazole fungicides are the most widely used products to treat cereal seeds. Granivorous birds, such as red-legged partridges (Alectoris rufa), which consume seeds left on the surface of fields after sowing, have a high risk of exposure. As triazole fungicides can affect sterol synthesis, we tested the hypothesis that treated seed consumption could alter the synthesis of sex hormones and reduce the reproductive capacity of partridges. We exposed adult partridges to seeds treated with four different formulations containing triazoles as active ingredients (flutriafol, prothioconazole, tebuconazole, and a mixture of the latter two) simulating a field exposure during the late autumn sowing season. All treatments produced biochemical changes and an overexpression of genes encoding for enzymes involved in the biosynthesis of sterols and steroid hormones, such as PMVK, ABCA1, MVD, PSCK9, DHCR7 and HSD17B7. Plasma levels of oestradiol were reduced in partridges exposed to tebuconazole. We also monitored reproduction 3 months after exposure (laying date, egg fertilization and hatching rates). We observed a 14-day delay in the laying onset of partridges that had been exposed to flutriafol as compared to controls. These results show that the consumption of seeds treated with triazole fungicides has the potential to affect granivorous bird reproduction. We recommend the evaluation of lagged reproductive effects as part of the protocols of environmental risk assessment of pesticides in wild birds in light of the effects resulting from the exposure to triazole-treated seeds.

Prediction of Soil Adsorption Coefficient in Pesticides Using Physicochemical Properties and Molecular Descriptors by Machine Learning Models

The soil adsorption coefficient (K_OC ) plays an important role in environmental risk assessment of pesticide registration. Based on this risk assessment, applied and registered pesticides can be allowed in the European Union. Almost 1 yr is required to study and obtain the K_OC value of a pesticide. Furthermore, acquiring the K_OC requires a large cost. It is necessary to efficiently estimate the K_OC value in the early stages of pesticide development. In the present study, the experimental values of physicochemical properties and molecular descriptors of chemical structures were collected to develop a quantitative structure-property relationship (QSPR) model, and the prediction performance of the model was evaluated. More specifically, we compared the accuracies of models based on a gradient boosting decision tree, multiple linear regression, and support vector machine. The experimental results suggest that it is possible to develop a QSPR model with high accuracy using both the molecular descriptors calculated from the structural formula and experimental values of physicochemical properties from open literature and databases. Comparing to the previously established models, we achieved high prediction accuracy, fitness, and robustness by only using freeware. Therefore, our developed QSPR models can be useful preliminary risk assessment in the early developmental stages of pesticides. Environ Toxicol Chem 2020;39:1451-1459. © 2020 SETAC.

A Probabilistic Co-Occurrence Approach for Estimating Likelihood of Spatial Overlap Between Listed Species Distribution and Pesticide Use Patterns

Characterizing potential spatial overlap between federally threatened and endangered ("listed") species distributions and registered pesticide use patterns is important for accurate risk assessment of threatened and endangered species. Because accurate range information for such rare species is often limited and agricultural pesticide use patterns are dynamic, simple spatial co-occurrence methods may overestimate or underestimate overlap and result in decisions that benefit neither listed species nor the regulatory process. Here, we demonstrate a new method of co-occurrence analysis that employs probability theory to estimate spatial distribution of rare species populations and areas of pesticide use to determine the likelihood of potential exposure. Specifically, we 1) describe a probabilistic method to estimate pesticide use based on crop production patterns; 2) construct species distribution models for 2 listed insect species whose ranges were previously incompletely described, the rusty-patched bumble bee (Bombus affinis) and the Poweshiek skipperling (Oarisma poweshiek); and 3) develop a probabilistic co-occurrence methodology and assessment framework. Using the principles of the Bayes' theorem, we constructed probabilistic spatial models of pesticide use areas by integrating information from land-cover spatial data, agriculture statistics, and remote-sensing data. We used maximum entropy methods to build species distribution models for 2 listed insects based on species collection and observation records and predictor variables relevant to the species' biogeography and natural history. We further developed novel methods for refinement of these models at spatial scales relevant to US Fish and Wildlife Service (FWS) regulatory priorities (e.g., critical habitat areas). Integrating both probabilistic assessments and focusing on USFWS priority management areas, we demonstrate that spatial overlap (i.e., potential for exposure) is not deterministic but instead a function of both species distribution and land use patterns. Our work serves as a framework to enhance the accuracy and efficiency of threatened and endangered species assessments using a data-driven likelihood analysis of species co-occurrence. Integr Environ Assess Manag 2019;00:1-12. © 2019 SETAC.

Time-Variable Exposure Experiments in Conjunction with Higher Tier Population and Effect Modeling to Assess the Risk of Chlorotoluron to Green Algae

An algae population model was applied to describe measured effects of pulsed exposure to chlorotoluron on populations of Pseudokirchneriella subcapitata in 2 laboratory flow-through chemostat tests with different exposure regimes. Both tests enabled evaluation of adverse effects on algae during the exposure and population recovery afterward. Impacts on population densities after chlorotoluron exposure were directly visible as biomass loss in the chemostats. Recovery was observed after each exposure peak. The test results indicate that P. subcapitata is unlikely to show an increased sensitivity to chlorotoluron after pulsed exposure. No altered response or adaptation of the algae to chlorotoluron was observed, with the exception of the last high peak in flow-through test 2. Therefore, an adaptation to the test substance cannot be excluded after long-term exposure. However, recovery to the steady-state level after this peak indicates that the growth rate (fitness) was not significantly reduced in the population with higher tolerance. No differences in chlorotoluron impact on the populations over time in terms of growth were detected. Model predictions agreed well with the measured data. The tests and modeling results validate the model to simulate population dynamics of P. subcapitata after pulsed exposure to chlorotoluron. Model predictions and extrapolations with different exposure patterns are considered reliable for chlorotoluron. The good reproducibility of the population behavior in the test systems supports this conclusion. An example modeled extrapolation of the experimental results to other (untested) exposure scenarios shows a potential approach to using the validated model as a supportive tool in risk assessment. Environ Toxicol Chem 2019;38:2520-2534. © 2019 SETAC.

Selecting appropriate focal species for assessing the risk to birds from newly drilled pesticide-treated winter cereal fields in France

Identifying focal bird species appropriate to the situation in which a plant protection product is used is important for refined risk assessment (EFSA ). We analyzed the results of extensive field observations of newly drilled cereal fields in France in autumn over 2 seasons to determine real bird focal species. In 2011, birds were observed before and after drilling on wheat and barley fields drilled with imidacloprid-treated seeds (i.e., "treatment" fields) or seeds treated with compounds other than imidacloprid (i.e., "alternative treatment" fields). Bird abundance, species richness, and diversity were significantly higher in wheat fields than barley fields; these findings led us to monitor only wheat fields in 2012. Statistical analyses did not show a significant effect of the drilling itself or between the treatment fields and the alternative treatment fields on the number and type of bird species. These results led to the pooling of 2011 data on all fields for focal species determination. Similarly, all bird monitoring data generated in 2012 before and after drilling were pooled and analyzed. Rules for determination of candidate focal species detailed in the EFSA () guidance were followed. Carrion crow, wood pigeon, gray partridge, skylark, common starling, and pied wagtail were the bird species most frequently observed on wheat fields. This list of candidate species was processed to determine the most relevant focal species according to the method of Dietzen et al. (); this process resulted in the selection of skylark, gray partridge, wood pigeon, and pied wagtail as focal species to assess risks to birds for pesticides applied during drilling of winter cereals in France (September through November). Integr Environ Assess Manag 2019;00:000-000. © 2018 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals, Inc. on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Thiamethoxam: Long-term effects following honey bee colony-level exposure and implications for risk assessment

Neonicotinoid insecticides have been used in a wide range of crops through seed treatment, soil and foliar applications and a large database exists on both their lethal and sub-lethal effects on honey bees under controlled laboratory conditions. However, colony-level studies on the effects of neonicotinoids in field studies are limited, primarily due to their complexity and the resources required. This paper reports the combined results of two large-scale colony-feeding studies, each with 6 weeks of continuous dosing of 12 colonies per treatment (24 control) to 12.5, 25, 37.5, 50 or 100 ng thiamethoxam/g sucrose solution. Exposure continued beyond dosing with residues present in stored nectar and bee-bread. The studies were conducted in an area with limited alternative forage and colonies were required to forage for pollen and additional nectar The studies provide colony-level endpoints: significant effects (reductions in bees, brood) were observed after exposure to the two highest dose rates, colony loss occurred at the highest dose rate, but colonies were able to recover (2-3 brood cycles after the end of dosing) after dosing with 50 ng thiamethoxam/g sucrose. No significant colony-level effects were observed at lower dose rates. The data reported here support the conclusions of previous colony-level crop-based field studies with thiamethoxam, in which residues in pollen and nectar were an order of magnitude below the colony-level NOEC of 37.5 ng thiamethoxam/g sucrose. The feeding study data are also compared to the outcomes of regulatory Tier 1 risk assessments conducted using guidance provided by the USA, Canada, Brazil and the EU regulatory authorities. We propose an adaptation of the European chronic adult bee risk assessment that takes into account the full dataset generated in laboratory studies while still providing an order of magnitude of safety compared with the colony feeding study NOEC.

Target site model: Predicting mode of action and aquatic organism acute toxicity using Abraham parameters and feature-weighted k-nearest neighbors classification

A database of 1480 chemicals with 47 associated modes of action compiled from the literature encompasses a wide range of chemical classes (alkanes, polycyclic aromatic hydrocarbons, pesticides, and polar compounds) and includes toxicity data for 79 different aquatic genera. The data were split into a calibration group and a validation group (80/20) to apply k-nearest neighbors (k-NN) methodology to predict the toxic mode of action for the compound. Other approaches were tested (support vector machines and linear discriminant analysis) as well as variations in the k-NN technique (distance weighting, feature weighting). Best-prediction results were found with k = 3, in a voting platform with optimized feature weighting. Using the predicted mode of action, the appropriate polyparameter target site model for that mode of action is applied to calculate the 50% lethal concentration (LC50). Predicted LC50s for the validation database resulted in a root-mean squared error (RMSE) of 0.752. This can be compared to an RMSE of 0.655 for the same validation set using the reference mode of action labels. The complete database resulted in an RMSE of 0.793 for reference mode of action labels. This confirms that the classification model has sufficient accuracy for predicting the mode of action and for determining toxicity using the target site model. Environ Toxicol Chem 2019;38:375-386. © 2018 SETAC.

Honey bee colony-level exposure and effects in realistic landscapes: An application of BEEHAVE simulating clothianidin residues in corn pollen

Discerning potential effects of insecticides on honey bee colonies in field studies conducted under realistic conditions can be challenging because of concurrent interactions with other environmental conditions. Honey bee colony models can control exposures and other environmental factors, as well as assess links among pollen and nectar residues in the landscape, their influx into the colony, and the resulting exposures and effects on bees at different developmental stages. We extended the colony model BEEHAVE to represent exposure to the insecticide clothianidin via residues in pollen from treated cornfields set in real agricultural landscapes in the US Midwest. We assessed their potential risks to honey bee colonies over a 1-yr cycle. Clothianidin effects on colony strength were only observed if unrealistically high residue levels in the pollen were simulated. The landscape composition significantly impacted the collection of pollen (residue exposure) from the cornfields, resulting in higher colony-level effects in landscapes with lower proportions of semi-natural land. The application of the extended BEEHAVE model with a pollen exposure-effects module provides a case study for the application of a mechanistic honey bee colony model in pesticide risk assessment integrating the impact of a range of landscape compositions. Environ Toxicol Chem 2019;38:423-435. © 2018 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals, Inc. on behalf of SETAC.

Target site model: Application of the polyparameter target lipid model to predict aquatic organism acute toxicity for various modes of action

A database of 2049 chemicals with 47 associated modes of action (MoA) was compiled from the literature. The database includes alkanes, polycyclic aromatic hydrocarbons, pesticides, inorganic, and polar compounds. Brief descriptions of some critical MoA classification groups are provided. The MoA from the 14 sources were assigned using a variety of reliable experimental and modeling techniques. Toxicity information, chemical parameters, and solubility limits were combined with the MoA label information to create the data set used for model development. The model database was used to generate linear free energy relationships for each specific MoA using multilinear regression analysis. The model uses chemical-specific Abraham solute parameters estimated from AbSolv to determine MoA-specific solvent parameters. With this procedure, critical target site concentrations are determined for each genus. Statistical analysis showed a wide range in values of the solvent parameters for the significant MoA. Environ Toxicol Chem 2019;38:222-239. © 2018 SETAC.

Exposures of children to neonicotinoids in pine wilt disease control areas

Neonicotinoid insecticides that have been on the market since 1992 have been used globally including in Japan. Because they are sprayed over forests and agricultural areas, inadvertent toxicity in nontarget insects (especially honey bees) and humans is a matter of public concern. However, information on exposure levels and potential health impacts of neonicotinoids in children living around sprayed areas is scarce. Thus, we determined neonicotinoid exposure levels in children living in communities where thiacloprid was used to control pine wilt disease. A total of 46 children (23 males and 23 females) were recruited for the present study, and informed written consent was obtained from their guardians. Urine specimens were collected before, during, and after insecticide spraying events; and atmospheric particulate matter was also collected. Concentrations of thiacloprid and 6 other neonicotinoid compounds were determined in urine samples and in atmospheric particulate matter specimens using liquid chromatography-electrospray ionization-tandem mass spectrometry. In urine specimens, thiacloprid concentrations were <0.13 μg/L and were detectable in approximately 30% of all samples. Concentrations of the other neonicotinoids, N-dm-acetamiprid, thiamethoxam, dinotefuran, and clothianidin, were 18.7, 1.92, 72.3, and 6.02 µg/L, respectively. Estimated daily intakes of these neonicotinoids were then calculated from urinary levels; although the estimated daily intakes of the neonicotinoids were lower than current acceptable daily intake values, the children were found to be exposed to multiple neonicotinoids on a daily basis. Environ Toxicol Chem 2019;38:71-79. © 2018 SETAC.

Coupling toxicokinetic-toxicodynamic and population models for assessing aquatic ecological risks to time-varying pesticide exposures

Population modeling evaluations of pesticide exposure time series were compared with aspects of a currently used risk assessment process. The US Environmental Protection Agency's Office of Pesticide Programs models daily aquatic 30-yr pesticide exposure distributions in its risk assessments, but does not routinely make full use of the information in such time series. We used mysid shrimp Americamysis bahia toxicity and demographic data to demonstrate the value of a toxicokinetic-toxicodynamic model coupled with a series of matrix population models in risk assessment refinements. This species is a small epibenthic marine crustacean routinely used in regulatory toxicity tests. We demonstrate how the model coupling can refine current risk assessments using only existing standard regulatory toxicity test results. Several exposure scenarios (each with the same initial risk characterization as determined by a more traditional organism-based approach) were created within which population modeling documented risks different from those of assessments based on the traditional approach. We also present different acute and chronic toxicity data scenarios by which toxicokinetic-toxicodynamic coupled with population modeling can distinguish responses that traditional risk evaluations are not designed to detect. Our results reinforce the benefits of this type of modeling in risk evaluations, especially related to time-varying exposure concentrations. Environ Toxicol Chem 2018;37:2633-2644. Published 2018 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

Focal Species Candidates for Pesticide Risk Assessment in European Rice Fields: A Review

An assessment of potential risks of pesticides on wildlife is required during the process of product registration within Europe because of the importance of agricultural landscapes as wildlife habitats. Despite their peculiarity and their specific role as artificial wetlands, rice paddies are to date pooled with cereals in guidance documents on how to conduct risk assessments for birds and mammals in Europe. Hence, the focal species currently considered in risk assessments for rice paddies are those known from cereal fields and can therefore be expected to differ significantly from the species actually occurring in the wet environments of rice paddies. We present results of a comprehensive review on bird and mammal species regularly occurring in rice paddies during a time of potential pesticide exposure to identify appropriate focal species candidates for ecotoxicological pesticide risk assessment according to the European Food Safety Authority (EFSA). In addition, we present data on rice cultivation areas and agricultural practices in Europe to give background information supporting the species selection process. Our literature search identified a general scarcity of relevant data, particularly for mammals, which highlights the need for crop-specific focal species studies. However, our results clearly indicate that the relevant bird and mammal species in rice fields indeed differ strongly from the focal species used for the cereal risk assessment. They can thus be used as a baseline for more realistic wildlife risk assessments specific to rice and the development of a revised guidance document to bridge the gap for regulatory decision makers. Integr Environ Assess Manag 2018;14:537-551. © 2018 SETAC.

Assessing and mitigating simulated population-level effects of 3 herbicides to a threatened plant: Application of a species-specific population model of Boltonia decurrens

Extrapolating from organism-level endpoints, as generated from standard pesticide toxicity tests, to populations is an important step in threatened and endangered species risk assessments. We apply a population model for a threatened herbaceous plant species, Boltonia decurrens, to estimate the potential population-level impacts of 3 herbicides. We combine conservative exposure scenarios with dose-response relationships for growth and survival of standard test species and apply those in the species-specific model. Exposure profiles applied in the B. decurrens model were estimated using exposure modeling approaches. Spray buffer zones were simulated by using corresponding exposure profiles, and their effectiveness at mitigating simulated effects on the plant populations was assessed with the model. From simulated exposure effects scenarios that affect plant populations, the present results suggest that B. decurrens populations may be more sensitive to exposures from herbicide spray drift affecting vegetative stages than from runoff affecting early seedling survival and growth. Spray application buffer zones were shown to be effective at reducing effects on simulated populations. Our case study demonstrates how species-specific population models can be applied in pesticide risk assessment to bring organism-level endpoints, exposure assumptions, and species characteristics together in an ecologically relevant context. Environ Toxicol Chem 2018;37:1545-1555. © 2018 SETAC.

Uptake and dissipation of neonicotinoid residues in nectar and foliage of systemically treated woody landscape plants

Systemic neonicotinoid insecticides used in urban arboriculture could pose a risk to bees and other pollinators foraging on treated plants. We measured uptake and dissipation of soil-applied imidacloprid and dinotefuran in nectar and leaves of 2 woody plant species, a broadleaf evergreen tree (Ilex × attenuata) and a deciduous shrub (Clethra alnifolia), to assess concentrations to which pollinators and pests might be exposed in landscape settings. Three application timings, autumn (postbloom), spring (prebloom), and summer (early postbloom), were evaluated to see if taking advantage of differences in the neonicotinoids' systemic mobility and persistence might enable pest control while minimizing transference into nectar. Nectar and tissue samples were collected from in-ground plants and analyzed for residues by high-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) in 2 successive years. Concentrations found in nectar following autumn or spring applications ranged from 166 to 515 ng/g for imidacloprid and from 70 to 1235 ng/gg for dinotefuran, depending on plant and timing. These residues exceed concentrations shown to adversely affect individual- and colony-level traits of bees. Summer application mitigated concentrations of imidacloprid (8-31 ng/g), but not dinotefuran (235-1191 ng/g), in nectar. Our data suggest that dinotefuran may be more persistent than is generally believed. Implications for integrated pest and pollinator management in urban landscapes are discussed. Environ Toxicol Chem 2018;37:860-870. © 2017 SETAC.

Thiamethoxam honey bee colony feeding study: Linking effects at the level of the individual to those at the colony level

Neonicotinoid insecticides have been used globally on a wide range of crops through seed treatment as well as soil and foliar applications and have been increasingly studied in relation to the potential risk to bees because of their detection in pollen and nectar of bee-attractive crops. The present article reports the results of laboratory studies (10-d adult and 22-d larval toxicity studies assessing the chronic toxicity of thiamethoxam to adult honey bees and larvae, respectively) and a colony feeding study, with 6 wk of exposure in an area with limited alternative forage, to provide a prewintering colony-level endpoint. The endpoints following exposure of individuals in the laboratory (10-d adult chronic no-observed-effect concentration [NOEC] for mortality 117 μg thiamethoxam/kg sucrose solution, 141 μg thiamethoxam/L sucrose solution; 22-d larval chronic NOEC 102 μg thiamethoxam/kg diet) are compared with those generated at the colony level, which incorporates sublethal effects (no-observed-adverse-effect concentration [NOAEC] 50 μg thiamethoxam/L sucrose solution, 43 μg thiamethoxam/kg sucrose solution). The data for sucrose-fed honey bee colonies support the lack of effects identified in previous colony-level field studies with thiamethoxam. However, unlike these field studies demonstrating no effects, colony feeding study data also provide a threshold level of exposure likely to result in adverse effects on the colony in the absence of alternative forage, and a basis by which to evaluate the potential risk of thiamethoxam residues detected in pollen, nectar, or water following treatment of bee-attractive crops. Environ Toxicol Chem 2018;37:816-828. © 2017 SETAC.

Developing population models: A systematic approach for pesticide risk assessment using herbaceous plants as an example

Population models are used as tools in species management and conservation and are increasingly recognized as important tools in pesticide risk assessments. A wide variety of population model applications and resources on modeling techniques, evaluation and documentation can be found in the literature. In this paper, we add to these resources by introducing a systematic, transparent approach to developing population models. The decision guide that we propose is intended to help model developers systematically address data availability for their purpose and the steps that need to be taken in any model development. The resulting conceptual model includes the necessary complexity to address the model purpose on the basis of current understanding and available data. We provide specific guidance for the development of population models for herbaceous plant species in pesticide risk assessment and demonstrate the approach with an example of a conceptual model developed following the decision guide for herbicide risk assessment of Mead's milkweed (Asclepias meadii), a species listed as threatened under the US Endangered Species Act. The decision guide specific to herbaceous plants demonstrates the details, but the general approach can be adapted for other species groups and management objectives. Population models provide a tool to link population-level dynamics, species and habitat characteristics as well as information about stressors in a single approach. Developing such models in a systematic, transparent way will increase their applicability and credibility, reduce development efforts, and result in models that are readily available for use in species management and risk assessments.

Response and recovery of the macrophytes Elodea canadensis and Myriophyllum spicatum following a pulse exposure to the herbicide iofensulfuron-sodium in outdoor stream mesocosms

Interest in stream mesocosms has recently revived for higher tier aquatic macrophyte risk assessment of plant protection products mainly because 1) the highest predicted environmental concentrations for the assessment of effects are frequently derived from stream scenarios, and 2) they allow an effect assessment using stream-typical pulse exposures. Therefore, the present stream mesocosm study used an herbicide pulse exposure and evaluated the responses of Elodea canadensis and Myriophyllum spicatum. Macrophytes were exposed for 24 h to 1 μg/L, 3 μg/L, 10 μg/L, and 30 μg/L of the herbicide iofensulfuron-sodium with a subsequent recovery period of 42 d. Biological endpoints were growth rates of the main, side, and total shoot length, the shoot number, the maximum root length, and the dry weight. The total shoot length was identified as the most sensitive endpoint; the growth rate of the total shoot length was inhibited by up to 66% and 45% in M. spicatum and E. canadensis, respectively. The lowest no observed effect concentrations (NOECs) were observed at day 7 and/or day 14 after herbicide treatment and were 1 μg/L for M. spicatum and 3 μg/L for E. canadensis. The no-observed-ecologically-adverse-effect concentrations (NOEAECs) were 10 μg/L and 30 μg/L for M. spicatum and E. canadensis, respectively. Such or similar mesocosm designs are useful to simulate typical stream exposures and estimate herbicide effects on aquatic macrophytes in stream systems. Environ Toxicol Chem 2017;36:1090-1100. © 2016 SETAC.

Mechanistic modeling of pesticide exposure: The missing keystone of honey bee toxicology

The role of pesticides in recent honey bee losses is controversial, partly because field studies often fail to detect effects predicted by laboratory studies. This dissonance highlights a critical gap in the field of honey bee toxicology: there exists little mechanistic understanding of the patterns and processes of exposure that link honey bees to pesticides in their environment. The authors submit that 2 key processes underlie honey bee pesticide exposure: 1) the acquisition of pesticide by foraging bees, and 2) the in-hive distribution of pesticide returned by foragers. The acquisition of pesticide by foraging bees must be understood as the spatiotemporal intersection between environmental contamination and honey bee foraging activity. This implies that exposure is distributional, not discrete, and that a subset of foragers may acquire harmful doses of pesticide while the mean colony exposure would appear safe. The in-hive distribution of pesticide is a complex process driven principally by food transfer interactions between colony members, and this process differs importantly between pollen and nectar. High priority should be placed on applying the extensive literature on honey bee biology to the development of more rigorously mechanistic models of honey bee pesticide exposure. In combination with mechanistic effects modeling, mechanistic exposure modeling has the potential to integrate the field of honey bee toxicology, advancing both risk assessment and basic research. Environ Toxicol Chem 2017;36:871-881. © 2016 SETAC.

Population modeling for pesticide risk assessment of threatened species-A case study of a terrestrial plant, Boltonia decurrens

Although population models are recognized as necessary tools in the ecological risk assessment of pesticides, particularly for species listed under the Endangered Species Act, their application in this context is currently limited to very few cases. The authors developed a detailed, individual-based population model for a threatened plant species, the decurrent false aster (Boltonia decurrens), for application in pesticide risk assessment. Floods and competition with other plant species are known factors that drive the species' population dynamics and were included in the model approach. The authors use the model to compare the population-level effects of 5 toxicity surrogates applied to B. decurrens under varying environmental conditions. The model results suggest that the environmental conditions under which herbicide applications occur may have a higher impact on populations than organism-level sensitivities to an herbicide within a realistic range. Indirect effects may be as important as the direct effects of herbicide applications by shifting competition strength if competing species have different sensitivities to the herbicide. The model approach provides a case study for population-level risk assessments of listed species. Population-level effects of herbicides can be assessed in a realistic and species-specific context, and uncertainties can be addressed explicitly. The authors discuss how their approach can inform the future development and application of modeling for population-level risk assessments of listed species, and ecological risk assessment in general. Environ Toxicol Chem 2017;36:480-491. © 2016 SETAC.

Use pattern of pesticides and their predicted mobility into shallow groundwater and surface water bodies of paddy lands in Mahaweli river basin in Sri Lanka

Pesticides applied on agricultural lands reach groundwater by leaching, and move to offsite water bodies by direct runoff, erosion and spray drift. Therefore, an assessment of the mobility of pesticides in water resources is important to safeguard such resources. Mobility of pesticides on agricultural lands of Mahaweli river basin in Sri Lanka has not been reported to date. In this context, the mobility potential of 32 pesticides on surface water and groundwater was assessed by widely used pesticide risk indicators, such as Attenuation Factor (AF) index and the Pesticide Impact Rating Index (PIRI) with some modifications. Four surface water bodies having greater than 20% land use of the catchment under agriculture, and shallow groundwater table at 3.0 m depth were selected for the risk assessment. According to AF, carbofuran, quinclorac and thiamethoxam are three most leachable pesticides having AF values 1.44 × 10^-2, 1.87 × 10^-3 and 5.70 × 10^-4, respectively. Using PIRI, offsite movement of pesticides by direct runoff was found to be greater than with the erosion of soil particles for the study area. Carbofuran and quinclorac are most mobile pesticides by direct runoff with runoff fractions of 0.01 and 0.08, respectively, at the studied area. Thiamethoxam and novaluron are the most mobile pesticides by erosion with erosion factions of 1.02 × 10^-4 and 1.05 × 10^-4, respectively. Expected pesticide residue levels in both surface and groundwater were predicted to remain below the USEPA health advisory levels, except for carbofuran, indicating that pesticide pollution is unlikely to exceed the available health guidelines in the Mahaweli river basin in Sri Lanka.

Using BEEHAVE to explore pesticide protection goals for European honeybee (Apis melifera L.) worker losses at different forage qualities

Losses of honeybee colonies are intensely debated and although honeybees suffer multiple stressors, the main focus has been on pesticides. As a result, the European Food Safety Authority (EFSA) revised the guidance for pesticide risk assessment for honeybees. The European Food Safety Authority reported a protection goal of negligible effect at 7% of colony size and then used the Khoury honeybee colony model to set trigger values for forager losses. However, the Khoury model is very simplistic and simulates colonies in an idealized state. In the present study, the authors demonstrate how a more realistic published honeybee model, BEEHAVE, with a few simple changes, can be used to explore pesticide risks. The results show that forage availability interacts with pesticide-induced worker losses, and colony resilience increases with forage quality. Adding alternative unexposed forage to the landscape also substantially mitigates the effects of pesticide exposure. The results indicate that EFSA's reported protection goal of 7% of colony size and triggers for daily worker losses are overly conservative. The authors conclude that forage availability is critical for colony resilience and that with adequate forage the colonies are resilient to even high levels of worker losses. However, the authors recommend setting protection goals using suboptimal forage conditions to ensure conservatism and for such suboptimal forage, a total of 20% reduction in colony size was safe. Environ Toxicol Chem 2017;36:254-264. © 2016 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals, Inc. on behalf of SETAC.

Comparison of buckwheat, red clover, and purple tansy as potential surrogate plants for use in semi-field pesticide risk assessments with Bombus impatiens

Background. Bumble bees (Bombus spp.) are important wild and managed pollinators. There is increased interest in incorporating data on bumble bees into risk assessments for pesticides, but standardized methods for assessing hazards of pesticides in semi-field and field settings have not yet been established for bumble bees. During semi-field studies, colonies are caged with pesticide-treated flowering surrogate plants, which must be attractive to foragers to ensure colony exposure to the test compound, and must produce an ample nectar and pollen to sustain colonies during testing. However, it is not known which plant(s) are suitable for use in semi-field studies with bumble bees.

Materials and Methods. We compared B. impatiens foraging activity and colony development on small plots of flowering buckwheat (Fagopyrum esculentum, var. common), red clover (Trifolium pratense), and purple tansy (Phacelia tanacetifolia) under semi-field conditions to assess their suitability as surrogate plants for pesticide risk assessment studies with bumble bees. We also compared the growth characteristics and input requirements of each plant type.

Results. All three plant types generally established and grew well. Red clover and purple tansy experienced significant weed pressure and/or insect pest damage. In contrast, pest pressure was extremely low in buckwheat. Overall, B. impatiens foraging activity was significantly greater on buckwheat plots than red clover or purple tansy, but plant type had no effect on number of individuals produced per colony or colony weight.

Discussion. Because of the consistently high foraging activity and successful colony development observed on buckwheat plots, combined with its favourable growth characteristics and low maintenance requirements, we recommend buckwheat as a surrogate plant for use in semi-field pesticide toxicity assessments with B. impatiens.