Mechanistic Effect Modeling of Earthworms in the Context of Pesticide Risk Assessment: Synthesis of the FORESEE Workshop

Intermediate-tier options in the environmental risk assessment of plant protection products for soil invertebrates-Synthesis of a workshop

The European environmental risk assessment (ERA) of plant protection products follows a tiered approach. The approach for soil invertebrates currently consists of two steps, starting with a Tier 1 assessment based on reproduction toxicity tests with earthworms, springtails, and predatory mites. In case an unacceptable risk is identified at Tier 1, field studies can be conducted as a higher-tier option. For soil invertebrates, intermediate tiers are not implemented. Hence, there is limited possibility to include additional information for the ERA to address specific concerns when the Tier 1 fails, as an alternative to, for example, a field study. Calibrated intermediate-tier approaches could help to address risks for soil invertebrates with less time and resources but also with sufficient certainty. A multistakeholder workshop was held on 2-4 March 2022 to discuss potential intermediate-tier options, focusing on four possible areas: (1) natural soil testing, (2) single-species tests (other than standard species), (3) assessing recovery in laboratory tests, and (4) the use of assembled soil multispecies test systems. The participants acknowledged a large potential in the intermediate-tier options but concluded that some issues need to be clarified before routine application of these approaches in the ERA is possible, that is, sensitivity, reproducibility, reliability, and standardization of potential new test systems. The definition of suitable assessment factors needed to calibrate the approaches to the protection goals was acknowledged. The aims of the workshop were to foster scientific exchange and a data-driven dialog, to discuss how the different approaches could be used in the risk assessment, and to identify research priorities for future work to address uncertainties and strengthen the tiered approach in the ERA for soil invertebrates. This article outlines the background, proposed methods, technical challenges, difficulties and opportunities in the ERA, and conclusions of the workshop. Integr Environ Assess Manag 2024;20:780-793. © 2023 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Evaluation of models to estimate the bioaccumulation of organic chemicals in earthworms

Modelling approaches to estimate the bioaccumulation of organic chemicals by earthworms are important for improving the realism in risk assessment of chemicals. However, the applicability of existing models is uncertain, partly due to the lack of independent datasets to test them. This study therefore conducted a comprehensive literature review on existing empirical and kinetic models that estimate the bioaccumulation of organic chemicals in earthworms and gathered two independent datasets from published literature to evaluate the predictive performance of these models. The Belfroid et al. (1995a) model is the best-performing empirical model, with 91.2% of earthworm body residue simulations within an order of magnitude of observation. However, this model is limited to the more hydrophobic pesticides and to the earthworm species Eisenia fetida or Eisenia andrei. The kinetic model proposed by Jager et al. (2003b) which out-performs that of Armitage and Gobas (2007), predicted uptake of PCB 153 in the earthworm E. andrei to within a factor of 10. However, the applicability of Jager et al.'s model to other organic compounds and other earthworm species is unknown due to the limited evaluation dataset. The model needs to be parameterised for different chemical, soil, and species types prior to use, which restricts its applicability to risk assessment on a broad scale. Both the empirical and kinetic models leave room for improvement in their ability to reliably predict bioaccumulation in earthworms. Whether they are fit for purpose in environmental risk assessment needs careful consideration on a case by case basis.

Ecotoxicological soil risk assessment under the new soil exposure framework - an impact assessment

For the European risk assessment (RA) for soil organisms exposed to plant protection products (PPPs) endpoints from ecotoxicological laboratory studies are compared with predicted environmental concentrations in soil (PECSOIL) at first tier. A safety margin must be met; otherwise, a higher tier RA is triggered (usually soil organism field studies). A new tiered exposure modeling guidance was published by EFSA to determine PECSOIL. This work investigates its potential impact on future soil RA. PECSOIL values for >50 active substances and metabolites were calculated and compared with the respective endpoints for soil organisms to calculate the RA failure rate. Compared to the current (FOCUS) exposure modeling, PECSOIL values for all EU regulatory zones considerably increased, e.g., resulting in active substance RA failure rates of 67%, 58% and 36% for modeling Tier-1, Tier-2 and Tier-3A, respectively. The main driving factors for elevated PECSOIL were soil bulk density, crop interception and wash-off, next to obligatory modeling and scenario adjustment factors. Spatial PECSOIL scenario selection procedures result in agronomically atypical soil characteristics (e.g., soil bulk density values in Tier-3A scenarios far below typical European agricultural areas). Consequently, exposure modeling and ecotoxicological study characteristics are inconsistent, which hinders scientifically reasonable comparison of both in the RA.

A spatiotemporally explicit modeling approach for more realistic exposure and risk assessment of off-field soil organisms

Natural and seminatural habitats of soil living organisms in cultivated landscapes can be subject to unintended exposure by active substances of plant protection products (PPPs) used in adjacent fields. Spray-drift deposition and runoff are considered major exposure routes into such off-field areas. In this work, we develop a model (xOffFieldSoil) and associated scenarios to estimate exposure of off-field soil habitats. The modular model approach consists of components, each addressing a specific aspect of exposure processes, for example, PPP use, drift deposition, runoff generation and filtering, estimation of soil concentrations. The approach is spatiotemporally explicit and operates at scales ranging from local edge-of-field to large landscapes. The outcome can be aggregated and presented to the risk assessor in a way that addresses the dimensions and scales defined in specific protection goals (SPGs). The approach can be used to assess the effect of mitigation options, for example, field margins, in-field buffers, or drift-reducing technology. The presented provisional scenarios start with a schematic edge-of-field situation and extend to real-world landscapes of up to 5 km × 5 km. A case study was conducted for two active substances of different environmental fate characteristics. Results are presented as a collection of percentiles over time and space, as contour plots, and as maps. The results show that exposure patterns of off-field soil organisms are of a complex nature due to spatial and temporal variabilities combined with landscape structure and event-based processes. Our concepts and analysis demonstrate that more realistic exposure data can be meaningfully consolidated to serve in standard-tier risk assessments. The real-world landscape-scale scenarios indicate risk hot-spots that support the identification of efficient risk mitigation. As a next step, the spatiotemporally explicit exposure data can be directly coupled to ecological effect models (e.g., for earthworms or collembola) to conduct risk assessments at biological entity levels as required by SPGs. Integr Environ Assess Manag 2024;20:263-278. © 2023 Applied Analysis Solutions LLC and WSC Scientific GmbH and Bayer AG and The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Genetic evidence behind the Cd resistance of wild Metaphire californica: The global RNA regulation rather than specific mutation of well-known gene

Heavy metal contamination presents a profound threat to terrestrial biodiversity, yet the genetic adaptation and evolution of field organisms under persistent stress are poorly understood. In this study, the Cd-resistant earthworms Metaphire californica collected from the control (Meihua, MHC) and elevated-pollution (Lupu, LPC) pairwise sites were used to elucidate the underlying genetic mechanism. A 48-h acute test showed that LPC worms exhibited 2.34 times higher LC50 (50% lethal concentration values) compared to MHC ones. The Cd bioaccumulation, metallothionein (MT) protein contents, and MT gene expression of LPC M.californica were all significantly higher than those of MHC worms. The well-known MT gene of M.californica was successfully cloned and identified, however, the encoding nucleotide and amino acids displayed non-observable mutations and the phylogenetic tree also revealed that different populations clustered together. Additionally, the results of transcriptomics sequencing demonstrated 173 differentially expressed genes between LPC and MHC worms, primarily involved in stress-response and detoxification pathways, including signal transduction, material metabolism, and protein exports. The above results confirmed that the crucial MT gene did not undergo genetic mutations but rather exhibited global mRNA regulation responsible for the Cd resistance of M.californica. The current study partially disclosed the stress adaptation and evolution of organisms under long-term in situ contamination, which provides insights into maintaining biodiversity under adverse environment.

Aquatic Risks at the Landscape Scale: A Case Study for Pyrethroid Use in Pome Fruit Orchards in Belgium

Procedures for environmental risk assessment for pesticides are under continuous development and subject to debate, especially at higher tier levels. Spatiotemporal dynamics of both pesticide exposure and effects at the landscape scale are largely ignored, which is a major flaw of the current risk assessment system. Furthermore, concrete guidance on risk assessment at landscape scales in the regulatory context is lacking. In this regard, we present an integrated modular simulation model system that includes spatiotemporally explicit simulation of pesticide application, fate, and effects on aquatic organisms. As a case study, the landscape model was applied to the Rummen, a river catchment in Belgium with a high density of pome fruit orchards. The application of a pyrethroid to pome fruit and the corresponding drift deposition on surface water and fate dynamics were simulated. Risk to aquatic organisms was quantified using a toxicokinetic/toxicodynamic model for individual survival at different levels of spatial aggregation, ranging from the catchment scale to individual stream segments. Although the derivation of landscape-scale risk assessment end points from model outputs is straightforward, a dialogue within the community, building on concrete examples as provided by this case study, is urgently needed in order to decide on the appropriate end points and on the definition of representative landscape scenarios for use in risk assessment.

Variation in the Chemical Sensitivity of Earthworms from Field Populations to Imidacloprid and Copper

The chemical risk of pesticides for nontarget soil macroorganisms has mainly been assessed using the compost earthworm Eisenia fetida. However, E. fetida does not occur in agroecosystems, and it is generally less sensitive than other earthworm species. Thus, the extrapolation of its response to pesticides to other earthworm species may lead to uncertainties in risk assessment. Because toxicity data for other earthworms are scarce, we assessed the chemical sensitivity of five species (Allolobophora chlorotica, Aporrectodea caliginosa, Aporrectodea longa, Aporrectodea rosea, and Lumbricus rubellus) from different habitats (forests, wetlands, and grasslands), as well as E. fetida, to imidacloprid and copper in single-species acute toxicity tests. In addition, we examined the relationship between earthworm traits (ecotype and weight), habitat characteristics (ecosystem type and soil pH), and chemical sensitivity. The lower limits of the hazardous concentration affecting 5% (HC5) of species were 178.99 and 0.32 mg active ingredient/kg dry weight for copper and imidacloprid, respectively. Some concentrations that have been measured in European agroecosystems for both pesticides were above the HC5s, indicating toxic risks for these organisms. Furthermore, soil pH from the sampling habitat played a significant role, with earthworms sampled from extremely acidic soils being less sensitive to copper than earthworms from neutral soils. In addition, endogeic earthworms were more sensitive to imidacloprid than epigeic earthworms. This may translate to changes in soil functions such as bioturbation, which is mainly carried out by endogeic earthworms. Our results suggest that risk assessment should include a wider range of earthworms covering different habitats and ecosystem functions to achieve a better protection of the biological functions carried out by these key soil organisms. Environ Toxicol Chem 2023;42:939-947. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Predicting the bioremediation potential of earthworms of different ecotypes through a multi-biomarker approach

Earthworms are attracting the attention of bioremediation research because of their short-term impact on pollutant fate. However, earthworm-assisted bioremediation largely depends on the earthworm sensitivity to target pollutants and its metabolic capacity to break down contaminants. The most studied species in soil bioremediation has been Eisenia fetida, which inhabits the soil surface feeding on decomposing organic residues. Therefore, its bioremediation potential may be limited to organic matter-rich topsoil. We compared the detoxification potential against organophosphate (OP) pesticides of three earthworm species representative of the main ecotypes: epigeic, anecic, and endogeic. Selected biomarkers of pesticide detoxification (esterases, cytochrome P450-dependent monooxygenase, and glutathione S-transferase) and oxidative homeostasis (total antioxidant capacity, glutathione levels, and glutathione reductase [GR] and catalase activities) were measured in the muscle wall and gastrointestinal tract of E. fetida (epigeic), Lumbricus terrestris (anecic) and Aporrectodea caliginosa (endogeic). Our results show that L. terrestris was the most suitable species to bioremediate OP-contaminated soil for the following reasons: 1) Gut carboxylesterase (CbE) activity of L. terrestris was higher than that of E. fetida, whereas muscle CbE activity was more sensitivity to OP inhibition than that of E. fetida, which means a high capacity to inactivate the toxic oxon metabolites of OPs. 2) Muscle and gut phosphotriesterase activities were significantly higher in L. terrestris than in the other species. 3) Enzymatic (catalase and GR) and molecular mechanisms of free radical inactivation (glutathione) were 3- to 4-fold higher in L. terrestris concerning E. fetida and A. caliginosa, which reveals a higher potential to keep the cellular oxidative homeostasis against reactive metabolites formed during OP metabolism. Together with biological and ecological traits, these toxicological traits suggest L. terrestris a better candidate for soil bioremediation than epigeic earthworms.

Revealing microcystin-LR ecotoxicity to earthworm (Eisenia fetida) at the intestinal cell level

Potential adverse effects of microcystin-LR (MC-LR) on soil invertebrates have not been studied. Here we investigated the mechanism of MC-LR toxicity to earthworm (Eisenia fetida) intestine at the individual level and at the cellular level. The results showed an inverse relationship between the bodyweight and survival rate of earthworms over exposure time- and MC-LR doses in soil. Dose-dependent intestinal lesions and disturbances of enzymatic activities (e.g., cellulase, Na⁺/K⁺-ATPase, and AChE) were observed, which resulted in intestinal dysfunction. Excessive reactive oxygen species generation led to DNA damage and lipid peroxidation of intestinal cells. The oxidative damage to DNA prolonged cell cycle arrest at the G2/M-phase transition in mitosis, thus stimulating and accelerating apoptosis in earthworm intestine. MC-LR target earthworm intestine tissue. MC-LR at low concentrations can damage earthworm intestine regardless of exposure routes (oral or contact). High toxicity of MC-LR to earthworms delineates its ecological risks to terrestrial ecosystems.

Three perspectives on the prediction of chemical effects in ecosystems

The increasing production, use and emission of synthetic chemicals into the environment represents a major driver of global change. The large number of synthetic chemicals, limited knowledge on exposure patterns and effects in organisms and their interaction with other global change drivers hamper the prediction of effects in ecosystems. However, recent advances in biomolecular and computational methods are promising to improve our capacity for prediction. We delineate three idealised perspectives for the prediction of chemical effects: the suborganismal, organismal and ecological perspective, which are currently largely separated. Each of the outlined perspectives includes essential and complementary theories and tools for prediction but captures only part of the phenomenon of chemical effects. Links between the perspectives may foster predictive modelling of chemical effects in ecosystems and extrapolation between species. A major challenge for the linkage is the lack of data sets simultaneously covering different levels of biological organisation (here referred to as biological levels) as well as varying temporal and spatial scales. Synthesising the three perspectives, some central aspects and associated types of data seem particularly necessary to improve prediction. First, suborganism- and organism-level responses to chemicals need to be recorded and tested for relationships with chemical groups and organism traits. Second, metrics that are measurable at many biological levels, such as energy, need to be scrutinised for their potential to integrate across levels. Third, experimental data on the simultaneous response over multiple biological levels and spatiotemporal scales are required. These could be collected in nested and interconnected micro- and mesocosm experiments. Lastly, prioritisation of processes involved in the prediction framework needs to find a balance between simplification and capturing the essential complexity of a system. For example, in some cases, eco-evolutionary dynamics and interactions may need stronger consideration. Prediction needs to move from a static to a real-world eco-evolutionary view.

The BEEHAVEecotox Model-Integrating a Mechanistic Effect Module into the Honeybee Colony Model

Mechanistic effect models are powerful tools for extrapolating from laboratory studies to field conditions. For bees, several good models are available that can simulate colony dynamics. Controlled and reliable experimental systems are also available to estimate the inherent toxicity of pesticides to individuals. However, there is currently no systematic and mechanistic way of linking the output of experimental ecotoxicological testing to bee models for bee risk assessment. We introduce an ecotoxicological module that mechanistically links exposure with the hazard profile of a pesticide for individual honeybees so that colony effects emerge. This mechanistic link allows the translation of results from standard laboratory studies to relevant parameters and processes for simulating bee colony dynamics. The module was integrated into the state-of-the-art honeybee model BEEHAVE. For the integration, BEEHAVE was adapted to mechanistically link the exposure and effects on different cohorts to colony dynamics. The BEEHAVEecotox model was tested against semifield (tunnel) studies, which were deemed the best study type to test whether BEEHAVEecotox predicted realistic effect sizes under controlled conditions. Two pesticides used as toxic standards were chosen for this validation to represent two different modes of action: acute mortality of foragers and chronic brood effects. The ecotoxicological module was able to predict effect sizes in the tunnel studies based on information from standard laboratory tests. In conclusion, the BEEHAVEecotox model is an excellent tool to be used for honeybee risk assessment, interpretation of field and semifield studies, and exploring the efficiency of different mitigation measures. The principles for exposure and effect modules are portable and could be used for any well-constructed honeybee model. Environ Toxicol Chem 2022;41:2870-2882. © 2022 Bayer AG & Sygenta, et al. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Integrating earthworm movement and life history through dynamic energy budgets

Earthworms are considered ecosystem engineers and, as such, they are an integral part of the soil ecosystem. The movement of earthworms is significantly influenced by environmental factors such as temperature and soil properties. As movement may directly be linked to food ingestion, especially of endogeic species like Aporrectodea caliginosa, changes in those environmental factors also affect life history traits such as growth and reproduction. In our laboratory studies, earthworms showed a decrease in burrowing activity with decreasing moisture levels and, to some extent, the organic matter content. The burrowing activity of earthworms was also affected by temperature, for which the casts produced per earthworm was used as a proxy in laboratory experiments. We integrated changes in earthworm movement and life histories in response to temperature, soil organic matter content and the moisture level, as observed in our experiment and reported in the literature, through dynamic energy budget (DEB) modelling. The joint parametrization of a DEB model for A. caliginosa based on movement and life history data revealed that food ingestion via movement is an integral part of the earthworms' energy budgets. Our findings highlight the importance of soil properties to be considered in the model development for earthworms. Furthermore, by understanding and incorporating the effect of environmental factors on the physiology, this mechanistic approach can help assess the impact of environmental changes such as temperature rise or drought.

A critical review of effect modeling for ecological risk assessment of plant protection products

A wide diversity of plant protection products (PPP) is used for crop protection leading to the contamination of soil, water, and air, which can have ecotoxicological impacts on living organisms. It is inconceivable to study the effects of each compound on each species from each compartment, experimental studies being time consuming and cost prohibitive, and animal testing having to be avoided. Therefore, numerous models are developed to assess PPP ecotoxicological effects. Our objective was to provide an overview of the modeling approaches enabling the assessment of PPP effects (including biopesticides) on the biota. Six categories of models were inventoried: (Q)SAR, DR and TKTD, population, multi-species, landscape, and mixture models. They were developed for various species (terrestrial and aquatic vertebrates and invertebrates, primary producers, micro-organisms) belonging to diverse environmental compartments, to address different goals (e.g., species sensitivity or PPP bioaccumulation assessment, ecosystem services protection). Among them, mechanistic models are increasingly recognized by EFSA for PPP regulatory risk assessment but, to date, remain not considered in notified guidance documents. The strengths and limits of the reviewed models are discussed together with improvement avenues (multigenerational effects, multiple biotic and abiotic stressors). This review also underlines a lack of model testing by means of field data and of sensitivity and uncertainty analyses. Accurate and robust modeling of PPP effects and other stressors on living organisms, from their application in the field to their functional consequences on the ecosystems at different scales of time and space, would help going toward a more sustainable management of the environment. Graphical Abstract Combination of the keyword lists composing the first bibliographic query. Columns were joined together with the logical operator AND. All keyword lists are available in Supplementary Information at https://doi.org/10.5281/zenodo.5775038 (Larras et al. 2021).

Ecotoxicological Effects of Bimetallic PdNi/MWCNT and PdCu/MWCNT Nanoparticles onto DNA Damage and Oxidative Stress in Earthworms

Bimetallic nanoparticles are synthesized using two different metal elements and used recently in many fields. However, limited studies related to the ecotoxic effects of nanoparticles available in the literature. The purpose of this study is to synthesize and characterize bimetallic PdCu/MWCNT and PdNi/MWCNT NPs and investigate their ecotoxic effects on earthworms. For this purpose, we injected approximately 20 µL of various concentrations of bimetallic PdCu/MWCNT and PdNi/MWCNT NPs (1, 10, 100, 1000, and 2000 mg/L) into the coelomic space of earthworms. We evaluated survival rate, malformations, reactive oxygen species (ROS) level, 8-OHdG content, and histopathological changes in earthworms at the 48th hour after exposure. PdCu/MWCNT and PdNi/MWCNT NPs were characterized by scanning electron microscope (SEM), transmission electron microscope (TEM), X-ray diffraction (XRD) pattern, and Raman-scattering spectroscopy. Toxicological examinations showed that PdCu/MWCNT NPs reduced the survival rate of earthworms (2000 mg/L, 84%) and caused various malformations (various lesions, thinning, swelling, and rupture), but nonsignificant effects of survival rate and malformations were observed in earthworms using PdNi/MWCNT NPs. The histopathological examinations of earthworm tissues exposed with PdNi/MWCNT determined that tissues in all treatment groups had a normal histological appearance. However, at a concentration of 2000 mg/L of PdCu/MWCNT NPs, atrophy in the longitudinal muscle layer and less degenerative cells in the epidermis layer were observed in earthworm tissues. It was determined that PdNi/MWCNT and PdCu/MWCNT NPs caused significant increases in ROS levels and 8-OHdG activity in earthworm tissues after 48 h. Finally, our results demonstrated that the toxicity of PdNi/MWCNT NPs was detected to be lower than PdCu/MWCNT NPs. However, both nanoparticles may pose a toxicological risk at high concentrations (1000 and 2000 mg/L). These findings will provide valuable information to studies on the use of PdNi/MWCNT NPs in wastewater treatment systems, industrial and medical fields, which have been determined to have less ecotoxicological risk.

Reveal the metal handling and resistance of earthworm Metaphire californica with different exposure history through toxicokinetic modeling

Toxicokinetic (TK) model provides a new approach to mechanistically elucidate the natural variation of metal handling strategy by adaptive and sensitive earthworm populations. Here, TK model was applied to explore the metal handling and resistance strategy of wild Metaphire californica with different historical exposure history through a 12-day re-exposure and another 12-day elimination incubation. M. californica populations showed different kinetic strategies for non-essential metals (Cd and Pb) and essential metals (Zn and Cu), which were closely related to their exposure history. M. californica from the most serious Cd-contaminated soil showed the fastest kinetic rates of both Cd uptake (K₁ = 0.78 g_soil/g_worm/day) and elimination (K₂ = 0.23 day^-1), and also had the lowest Cd half-life (t_1/2 = 3.01 day), which demonstrated the potential Cd-resistance of wild M. californica from Cd-contaminated soils. Besides, the comparative experiment showed totally different metal kinetics of laboratory Eisenia fetida from field M. californica, suggesting the impacts of distinct exposure history and species-specifical sensitivities. These findings provide a novel approach to identify and quantify resistance using TK model and also imply the risk of overlooking existing exposure background and interspecies extrapolation in eco-toxicological studies and risk assessments.