Dynamic energy budget models in ecological risk assessment: From principles to applications

Comparative temporal response of toxicity for the neonicotinoid clothianidin and organophosphate dimethoate insecticides in two species of solitary bee (Osmia bicornis and Osmia cornuta)

Solitary bees provide essential pollination services. Concerns for the decline of these wild bee species have led to calls for their inclusion in pesticide risk assessment. Solitary bees differ from honey bees in their physiology and ecology and this may affect how they respond to pesticide exposure. Here we investigate the life-time toxicity of two insecticides, the organophosphate dimethoate and neonicotinoid clothianidin, for two mason bee species, Osmia bicornis and O. cornuta using a toxicokinetic/toxicodynamic stochastic death model taken from Dynamic Energy Budget (DEBtox) theory. Both species showed concentration and exposure duration dependent effects for each chemical. LC50 values estimated from the model parameters at 48 h were ≥ 14 fold and 6 fold those at 480 h for dimethoate and clothianidin respectively. Survival modelling indicated greater sensitivity in O. bicornis than for O. cornuta to dimethoate, whilst for clothianidin, O. cornuta females but not males, were more sensitive than both sexes of O. bicornis. These sensitivity differences were not related to body size. Toxicokinetic and toxicodynamic traits derived from modelling indicated lower elimination rates in O. bicornis and higher killing rates for O. cornuta females for dimethoate and lower elimination rates for clothianidin in O. cornuta females that were related to sensitivity. This study shows the near life-time testing is possible for solitary bees and that combining adult life-time toxicity tests with toxicokinetic/toxicodynamic modelling provides a more mechanistic understanding of pesticide effects in solitary bee species.

Comparing freshwater mussel responses to stress using life-history and Dynamic Energy Budget theory

Freshwater mussels are experiencing severe population declines, affecting their critical role in freshwater ecosystems. A thorough assessment of threats posed by various stressors is needed; however, the large number of species to be considered and significant data gaps, especially for listed species, hinder the process. We combined a traits-based approach to represent multiple species grouped into three life-history categories - Equilibrium, Opportunistic, and Periodic - with the Dynamic Energy Budget modeling principles to capture the physiological mechanisms driving individual-level responses. We used the DEB model to simulate individual life cycles and explore relationships between underlying energetics and emerging individual traits of 47 freshwater mussel species and the common toxicity test surrogate, the Eastern oyster (Crassostrea virginica), under control and stressed conditions. Stress was introduced via physiological modes of action related to four key metabolic pathways: energy assimilation, maintenance, growth, and reproduction. We recorded maximum length, age at maturity, and fecundity and compared these endpoints and their stress-induced changes among life-history categories. The life-history differences among freshwater mussels directly emerged from underlying energetics, with high assimilation and maintenance supporting opportunistic traits. Stress imposed on energy assimilation had the strongest effect on all life-history traits, and a 25 % reduction in assimilation rate resulted in an average 25 % and 60 % decrease in maximum length and fecundity, respectively, and a 24 % increase in age at maturity. Equilibrium species suffered the greatest negative effects overall, indicating that this life-history strategy might be the most susceptible to stressors. The Eastern oyster displayed extreme opportunism in its life-history traits, but its responses to stress were generally within the range observed for freshwater mussels. The study provides a much-needed general understanding of stress responses across freshwater mussel life-history categories and contributes to the foundation for developing life-history-driven population models.

Bioremediation potential of the hard clam Mercenaria mercenaria as an intensive shrimp aquaculture pond polyculture condidate

Polyculture practices are important for achieving sustainable aquaculture development. Recently, hard clams polyculture in intensive shrimp ponds has been encouraged because bivalves can consume excess nutrients in aquaculture systems and sequester carbon. To evaluate the bioremediation potential of hard clams polyculture in intensive shrimp ponds, this study built an assessment model based on individual growth models and estimated the potential for nitrogen and phosphorus removal as well as CO2 fixation by hard clams. Firstly, key parameters required for model construction were obtained through field surveys and physiological experiments. Subsequently, an individual growth model for the hard clam Mercenaria mercenaria was developed based on the Dynamic Energy Budget (DEB) theory. Fitting of the growth data indicated that the model accurately replicated the growth patterns of hard clams, with relative root mean square errors of 9.87 % for shell length and 5.02 % for dry tissue weight. Finally, the assessment model for the bioremediation potential of hard clams demonstrated that, over 110 days in the intensive shrimp mariculture pond, the net removal of nitrogen and phosphorus by hard clams were 3.68 kg ha-1 and 0.81 kg ha-1, respectively, and CO2 fixation was 507.00 kg ha-1. These findings suggested that the DEB model is an effective tool for evaluating bivalve ecological remediation potential and can aid in selecting species for sustainable polyculture.

Key challenges and developments in wildlife ecological risk assessment: Problem formulation

Problem formulation (PF) is a critical initial step in planning risk assessments for chemical exposures to wildlife, used either explicitly or implicitly in various jurisdictions to include registration of new pesticides, evaluation of new and existing chemicals released to the environment, and characterization of impact when chemical releases have occurred. Despite improvements in our understanding of the environment, ecology, and biological sciences, few risk assessments have used this information to enhance their value and predictive capabilities. In addition to advances in organism-level mechanisms and methods, there have been substantive developments that focus on population- and systems-level processes. Although most of the advances have been recognized as being state-of-the-science for two decades or more, there is scant evidence that they have been incorporated into wildlife risk assessment or risk assessment in general. In this article, we identify opportunities to consider elevating the relevance of wildlife risk assessments by focusing on elements of the PF stage of risk assessment, especially in the construction of conceptual models and selection of assessment endpoints that target population- and system-level endpoints. Doing so will remain consistent with four established steps of existing guidance: (1) establish clear protection goals early in the process; (2) consider how data collection using new methods will affect decisions, given all possibilities, and develop a decision plan a priori; (3) engage all relevant stakeholders in creating a robust, holistic conceptual model that incorporates plausible stressors that could affect the targets defined in the protection goals; and (4) embrace the need for iteration throughout the PF steps (recognizing that multiple passes may be required before agreeing on a feasible plan for the rest of the risk assessment). Integr Environ Assess Manag 2024;20:658-673. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC). This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Influence of exposure scenario on the sensitivity to caffeine

The chorion acts as a protective barrier, restricting some chemical absorption into the embryo and the surrounding fluids. In this sense, larvae may only have direct contact with some chemicals after dechorionation. This study aimed to evaluate the effects of caffeine (CAF) (0, 13, 20, 44, 67, and 100 mg.L-1) under different exposure scenarios (embryos with chorion or embryos/larvae already hatched) and rank the stage sensitivity. Thus, three scenarios were investigated: from 2 to 120 hours post fertilization (hpf) (5 days of exposure- 5dE), from 72 to 120 hpf (2dE), and from 96 to 120 hpf (1dE). Heart rate (48 hpf) and energy reserves (120 hpf) were measured in the 5dE scenario, and behavior and acetylcholinesterase (AChE) activity were evaluated at 120 hpf in all scenarios (5dE, 2dE, and 1dE). At 120 hpf, some of the fish was transferred to clean medium for a 10 days depuration period (10dPE). Behavior and AChE activity were assessed after this period. In the 5dE scenario, CAF increased heartbeat (13, 20, and 30 mg.L-1) and reduced carbohydrates (67, and 100 mg.L-1), while inhibiting AChE activity (100 mg.L-1) in the 5dE, 2dE, and 1dE scenarios. CAF reduced the total distance moved in the 5dE (67, and 100 mg.L-1), 2dE (20, 30, 44, 67, and 100 mg.L-1), and 1dE fish (67, and 100 mg.L-1) and increased erratic movements. Based on the lowest observed effect concentration (LOEC) for total distance moved (20 mg.L-1) and higher inhibition of AChE activity (100 mg.L-1) (65%), 2dE fish appear to be more sensitive to CAF. After 10dPE, a recovery in behavior was detected in all scenarios (5dE, 2dE, and 1dE). AChE activity remained inhibited in the 2dE scenario while increasing in the 1dE scenario. This study demonstrated that the presence of the chorion is an important factor for the analysis of CAF toxicity. After the loss of the chorion, organisms show greater sensitivity to CAF and can be used to evaluate the toxicity of various substances, including nanomaterials or chemicals with low capacity to cross the chorion. Therefore, the use of hatched embryos in toxicity tests is suggested, as they allow a shorter and less expensive exposure scenario that provides similar outcome as the conventional scenario.

TKPlate 1.0: An Open-access platform for toxicokinetic and toxicodynamic modelling of chemicals to implement new approach methodologies in chemical risk assessment

This publication is linked to the following EFSA Supporting Publications articles: http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2023.EN-8441/full, http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2023.EN-8440/full, http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2023.EN-8437/full.

The use of new approach methodologies for the environmental risk assessment of food and feed chemicals

New Approach Methodologies (NAMs) provide tools for supporting both human and environmental risk assessment (HRA and ERA). This short review provides recent insights regarding the use of NAMs in ERA of food and feed chemicals. We highlight the usefulness of tiered methods supporting weight-of-evidence approaches in relation to problem formulation (i.e., data availability, time, and resource availability). In silico models, including quantitative structure activity relationship models, support filling data gaps when no chemical property or ecotoxicological data are available, and biologically-based models (e.g., toxicokinetic-toxicodynamic models, dynamic energy models, physiologically-based models and species sensitivity distributions) are applicable in more data rich situations, including landscape-based modelling approaches. Particular attention is given to provide practical examples to apply the approaches described in real-world settings. We conclude with future perspectives, with regards to the need for addressing complex challenges such as chemical mixtures and multiple stressors in a wide range of organisms and ecosystems.

Joint survival modelling for multiple species exposed to toxicants

In environmental risk assessment (ERA), the multitude of compounds and taxa demands cross-species extrapolation to cover the variability in sensitivity to toxicants. However, only the impact of a single compound to a single species is addressed by the general unified threshold model of survival (GUTS). The reduced GUTS is the recommended model to analyse lethal toxic effects in regulatory aquatic ERA. GUTS considers toxicokinetics and toxicodynamics. Two toxicodynamic approaches are considered: Stochastic death (SD) assumes that survival decreases with an increasing internalized amount of the toxicant. Individual tolerance (IT) assumes that individuals vary in their tolerance to toxic exposure. Existing theory suggests that the product of the threshold z_w and killing rate b_w (both SD toxicodynamic parameters) are constant across species or compounds if receptors and target sites are shared. We extend that theory and show that the shape parameter β of the loglogistic threshold distribution in IT is also constant. To verify the predicted relationships, we conducted three tests using toxicity studies for eight arthropods exposed to the insecticide flupyradifurone. We confirmed previous verifications of the relation- between SD parameters, and the newly established relation for the IT parameter β. We enhanced GUTS to jointly model survival for multiple species with shared receptors and pathways by incorporating the relations among toxicodynamic parameters described above. The joint GUTS exploits the shared parameter relations and therefore constrains parameter uncertainty for each of the separate species. Particularly for IT, the joint GUTS more precisely predicted risk to the separate species than the standard single species GUTS under environmentally realistic exposure. We suggest that joint GUTS modelling can improve cross-species extrapolation in regulatory ERA by increasing the reliability of risk estimates and reducing animal testing. Furthermore, the shared toxicodynamic response provides potential to reduce complexity of ecosystem models.

Modelling the effects of the pyrethroid insecticide cypermethrin on the life cycle of the soil dwelling annelid Enchytraeus crypticus, an original experimental design to calibrate a DEB-TKTD model

The Dynamic Energy Budget theory (DEB) enables ecotoxicologists to model the effects of chemical stressors on organism life cycles through the coupling of toxicokinetic-toxicodynamic (TK-TD) models. While good progress has been made in the application of DEB-TKTD models for aquatic organisms, applications for soil fauna are scarce, due to the lack of dedicated experimental designs suitable for collecting the required time series effect data. Enchytraeids (Annelida: Clitellata) are model organisms in soil ecology and ecotoxicology. They are recognised as indicators of biological activity in soil, and chemical stress in terrestrial ecosystems. Despite this, the application of DEB-TKTD models to investigate the impact of chemicals has not yet been tested on this family. Here we assessed the impact of the pyrethroid insecticide cypermethrin on the life cycle of Enchytraeus crypticus. We developed an original experimental design to collect the data required for the calibration of a DEB-TKTD model for this species. E. crypticus presented a slow initial growth phase that has been successfully simulated with the addition of a size-dependent food limitation for juveniles in the DEB model. The DEB-TKTD model simulations successfully agreed with the data for all endpoints and treatments over time. The highlighted physiological mode of action (pMoA) for cypermethrin was an increase of the growth energy cost. The threshold for effects on survival was estimated at 73.14 mg kg^- 1, and the threshold for effects on energy budget (i.e., sublethal effects) at 19.21 mg kg^- 1. This study demonstrates that DEB-TKTD models can be successfully applied to E. crypticus as a representative soil species, and may improve the ecological risk assessment for terrestrial ecosystems, and our mechanistic understanding of chemical effects on non-target species.

An annual profile of the impacts of simulated oil spills on the Northeast Arctic cod and haddock fisheries

We simulate the combined natural and pollutant-induced survival of early life stages of NEA cod and haddock, and the impact on the adult populations in response to the time of a major oil spill in a single year. Our simulations reveal how dynamic ocean processes, controlling both oil transport and fate and the frequency of interactions of oil with drifting fish eggs and larvae, mediate the magnitude of population losses due to an oil spill. The largest impacts on fish early life stages occurred for spills initiated in Feb-Mar, concomitant with the initial rise in marine productivity and the earliest phase of the spawning season. The reproductive health of the adult fish populations was maintained in all scenarios. The study demonstrates the application of a simulation system that provides managers with information for the planning of development activities and for the protection of fisheries resources from potential impacts.

A regression-based QSAR-model to predict acute toxicity of aromatic chemicals in tadpoles of the Japanese brown frog (Rana japonica): Calibration, validation, and future developments to support risk assessment of chemicals in amphibians

Amphibian populations are undergoing a global decline worldwide. Such decline has been attributed to their unique physiology, ecology, and exposure to multiple stressors including chemicals, temperature, and biological hazards such as fungi of the Batrachochytrium genus, viruses such as Ranavirus, and habitat reduction. There are limited toxicity data for chemicals available for amphibians and few quantitative structure-activity relationship (QSAR) models have been developed and are publicly available. Such QSARs provide important tools to assess the toxicity of chemicals particularly in a data poor context. QSARs provide important tools to assess the toxicity of chemicals particularly when no toxicological data are available. This manuscript provides a description and validation of a regression-based QSAR model to predict, in a quantitative manner, acute lethal toxicity of aromatic chemicals in tadpoles of the Japanese brown frog (Rana japonica). QSAR models for acute median lethal molar concentrations (LC50-12 h) of waterborne chemicals using the Monte Carlo method were developed. The statistical characteristics of the QSARs were described as average values obtained from five random distributions into training and validation sets. Predictions from the model gave satisfactory results for the overall training set (R² = 0.72 and RMSE = 0.33) and were even more robust for the validation set (R² = 0.96 and RMSE = 0.11). Further development of QSAR models in amphibians, particularly for other life stages and species, are discussed.

Molting enhances internal concentrations of fipronil and thereby decreases survival of two estuarine resident marine crustaceans

In aquatic arthropods, molting is a pivotal physiological process for normal development, but it may also expose them to higher risks from xenobiotics, because the organism may take up additional water during that time. This study aimed to assess the effects of molting on bioconcentration and survival after 96-h exposure to insecticide fipronil with or without oxygenase (CYP450s) inhibitor piperonyl butoxide (PBO) of two estuarine resident marine crustacean species: the sand shrimp Crangon uritai and the kuruma prawn Penaeus japonicus, with 96-h LC₅₀ value of fipronil = 2.0 µg/L and 0.2 µg/L, respectively. Two graded concentrations included group high (H) (equivalent to the 96-h LC₅₀ values) and low (L) (one-tenth of the H group concentration). Molting and survival were individually checked, and internal concentrations of fipronil and its metabolites (fipronil desulfinyl, fipronil sulfide, fipronil sulfone) were measured. The results showed that, only fipronil and fipronil sulfone were detected from organism, and that internal concentrations of these insecticides in molted specimens were higher than those of unmolted ones but comparable with those of dead ones. Accordingly, mortality was more frequent in molted specimens than those that were unmolted. Furthermore, involvement of oxygenase and higher lethal body burden threshold may confer higher tolerance to fipronil in sand shrimp than in the kuruma prawn. This study is the first to demonstrate that the body-residue-based approach is useful for deciphering the causal factors underlying fipronil toxicity, but highlights the need to consider physiological factors in arthropods, which influence and lie beyond body burden, molting and drug metabolism.

The regulating effect of growth plasticity on the dynamics of structured populations

Plasticity is the extent to which life history processes such as growth and reproduction depend on the environment. Plasticity in individual growth varies widely between taxa. Nonetheless, little is known about the effect of plasticity in individual growth on the ecological dynamics of populations. In this article, we analyse a physiologically structured population model of a consumer population in which the individual growth rate can be varied between entirely plastic to entirely non-plastic. We derive this population level model from a dynamic energy budget model to ensure an accurate energetic coupling between ingestion, somatic maintenance, growth and reproduction within an individual. We show that the consumer population is either limited by adult fecundity or juvenile survival up to maturation, depending on the level of growth plasticity and the non-plastic individual growth rate. Under these two regimes, we also find two different types of population cycles which again arise due to fluctuation in, respectively, juvenile growth rate or adult fecundity. In the end, our model not only provides insight into the effects of growth plasticity on population dynamics, but also provides a link between the dynamics found in age- and size-structured models.

In Silico Methods for Environmental Risk Assessment: Principles, Tiered Approaches, Applications, and Future Perspectives

This chapter aims to introduce the reader to the basic principles of environmental risk assessment of chemicals and highlights the usefulness of tiered approaches within weight of evidence approaches in relation to problem formulation i.e., data availability, time and resource availability. In silico models are then introduced and include quantitative structure-activity relationship (QSAR) models, which support filling data gaps when no chemical property or ecotoxicological data are available. In addition, biologically-based models can be applied in more data rich situations and these include generic or species-specific models such as toxicokinetic-toxicodynamic models, dynamic energy budget models, physiologically based models, and models for ecosystem hazard assessment i.e. species sensitivity distributions and ultimately for landscape assessment i.e. landscape-based modeling approaches. Throughout this chapter, particular attention is given to provide practical examples supporting the application of such in silico models in real-world settings. Future perspectives are discussed to address environmental risk assessment in a more holistic manner particularly for relevant complex questions, such as the risk assessment of multiple stressors and the development of harmonized approaches to ultimately quantify the relative contribution and impact of single chemicals, multiple chemicals and multiple stressors on living organisms.

(Eco)toxicological tests for assessing impacts of chemical stress to aquatic ecosystems: Facts, challenges, and future

Monitoring of chemicals in the aquatic environment by chemical analysis alone cannot completely assess and predict the effects of chemicals on aquatic species and ecosystems. This is primarily because of the increasing number of (unknown) chemical stressors and mixture effects present in the environment. In addition, the ability of ecological indices to identify underlying stressors causing negative ecological effects is limited. Therefore, additional complementary methods are needed that can address the biological effects in a direct manner and provide a link to chemical exposure, i.e. (eco)toxicological tests. (Eco)toxicological tests are defined as test systems that expose biological components (cells, individuals, populations, communities) to (environmental mixtures of) chemicals to register biological effects. These tests measure responses at the sub-organismal (biomarkers and in vitro bioassays), whole-organismal, population, or community level. We performed a literature search to obtain a state-of-the-art overview of ecotoxicological tests available for assessing impacts of chemicals to aquatic biota and to reveal datagaps. In total, we included 509 biomarkers, 207 in vitro bioassays, 422 tests measuring biological effects at the whole-organismal level, and 78 tests at the population- community- and ecosystem-level. Tests at the whole-organismal level and biomarkers were most abundant for invertebrates and fish, whilst in vitro bioassays are mostly based on mammalian cell lines. Tests at the community- and ecosystem-level were almost missing for organisms other than microorganisms and algae. In addition, we provide an overview of the various extrapolation challenges faced in using data from these tests and suggest some forward looking perspectives. Although extrapolating the measured responses to relevant protection goals remains challenging, the combination of ecotoxicological experiments and models is key for a more comprehensive assessment of the effects of chemical stressors to aquatic ecosystems.

Protein corona on biogenic silver nanoparticles provides higher stability and protects cells from toxicity in comparison to chemical nanoparticles

The development of environmentally friendly new procedures for the synthesis of metallic nanoparticles is one of the main goals of nanotechnology. Proteins and enzymes from plants, filamentous fungi, yeast, and bacteria to produce nanoparticles are both valuable and viable alternatives to conventional synthesis of nanomaterials due to their high efficiency and the low cost to scale up and generate large quantities. The aim of this work is to compare biogenic silver nanoparticles (AgNPs) obtained from cell-free filtrates from the fungus Macrophomina phaseolina to conventional chemical AgNPs, in biocidal activity and toxicity. Our results show that bio-AgNPs displayed similar bactericidal activity than chemical AgNPs, but less toxicity in the model organism Caenorhabditis elegans. We employed biochemical and proteomic techniques to profile the unique surface chemistry of the capping in the bio-AgNPs and therefore to identify the proteins involved in their synthesis and stability. These results not only suggest that the proteins involved in the synthesis of the nanoparticles and corona formation in the bio-AgNPs are responsible for keeping the silver core preserved making them more stable in time, but also masking and protecting eukaryotic cells from metal toxicity.

OpenCYP: An open source database exploring human variability in activities and frequencies of polymophisms for major cytochrome P-450 isoforms across world populations

The open source database "OpenCYP database" has been developed based on the results of extensive literature searches from the peer-reviewed literature. OpenCYP provides data on human variability on baseline of activities and polymophism frequencies for selected cytochrome P-450 isoforms (CYP1A2, CYP2A6, CYP2D6, CYP3A4/3A5 and CYP3A7) in healthy adult populations from world populations. CYP enzymatic activities were generally expressed as the metabolic ratio (MR) between an unchanged probe drug and its metabolite(s) in urine or plasma measured in healthy adults. Data on other age groups were very limited and fragmented, constituting an important data gap. Quantitative comparisons were often hampered by the different experimental conditions used. However, variability was quite limited for CYP1A2, using caffeine as a probe substrate, with a symmetrical distribution of metabolic activity values. For CYP3A4, human variability was dependent on the probe substrate itself with low variability when data considering the dextromethorphan/demethilathed metabolite MR were used and large variability when the urinary 6β-hydroxycortisol/cortisol ratio was used. The largest variability in CYP activity was shown for CYP2D6 activity, after oral dosing of dextromethorphan, for which genetic polymorphisms are well characterised and constitute a significant source of variability. It is foreseen that the OpenCYP database can contribute to promising tools to support the further development of QIVIVE and PBK models for human risk assessment of chemicals particularly when combined with information on isoform-specific content in cells using proteomic approaches.

Cross-species extrapolation of chemical sensitivity

Ecosystems are usually populated by many species. Each of these species carries the potential to show a different sensitivity towards all of the numerous chemical compounds that can be present in their environment. Since experimentally testing all possible species-chemical combinations is impossible, the ecological risk assessment of chemicals largely depends on cross-species extrapolation approaches. This review overviews currently existing cross-species extrapolation methodologies, and discusses i) how species sensitivity could be described, ii) which predictors might be useful for explaining differences in species sensitivity, and iii) which statistical considerations are important. We argue that risk assessment can benefit most from modelling approaches when sensitivity is described based on ecologically relevant and robust effects. Additionally, specific attention should be paid to heterogeneity of the training data (e.g. exposure duration, pH, temperature), since this strongly influences the reliability of the resulting models. Regarding which predictors are useful for explaining differences in species sensitivity, we review interspecies-correlation, relatedness-based, traits-based, and genomic-based extrapolation methods, describing the amount of mechanistic information the predictors contain, the amount of input data the models require, and the extent to which the different methods provide protection for ecological entities. We develop a conceptual framework, incorporating the strengths of each of the methods described. Finally, the discussion of statistical considerations reveals that regardless of the method used, statistically significant models can be found, although the usefulness, applicability, and understanding of these models varies considerably. We therefore recommend publication of scientific code along with scientific studies to simultaneously clarify modelling choices and enable elaboration on existing work. In general, this review specifies the data requirements of different cross-species extrapolation methods, aiming to make regulators and publishers more aware that access to raw- and meta-data needs to be improved to make future cross-species extrapolation efforts successful, enabling their integration into the regulatory environment.

Compartmentation and effects of lead (Pb) in the collembolan, Folsomia candida

The impact of soil lead (Pb) pollution on survival, growth, and reproduction of the collembolan, Folsomia candida, and Pb compartmentation in its gut and remaining body parts were studied by exposing animals to laboratory-spiked soil. The survival, growth, and reproduction of F. candida were significantly reduced by increasing soil Pb concentration. The LC₅₀ values of survival based on total and CaCl₂-extractable Pb concentration in soil were 2562 mg kg^-1 and 351 mg kg^-1, respectively. The EC₅₀ values of reproduction were 1244 mg kg^-1 and 48 mg kg^-1, respectively. The Pb concentration in whole body, gut, and remaining body parts was significantly increased with the increase of soil Pb concentration and followed an exponential increase when the soil Pb concentration was equal to or above a threshold (1000 mg kg^-1 for whole body and remaining body part, 500 mg kg^-1 for gut). Below this threshold, these relationships were linear. The Pb concentration in the gut was higher than whole body and remaining body part of F. candida, and the threshold of internal Pb concentration at which F. candida can compensate was in the range 7-13 mg Pb kg^-1 dry animal (corresponding to soil Pb concentration 500-1000 mg Pb kg^-1 dry soil). The results indicate that reproduction of F. candida was a more sensitive indicator of lead toxicity than survival and growth. Pb was mainly accumulated in the gut of F. candida. We discuss the internal Pb concentration as an indicator of adverse effects in the risk assessment of soil Pb pollution.

Sublethal effect modelling for environmental risk assessment of chemicals: Problem definition, model variants, application and challenges

Bioenergetic models, and specifically dynamic energy budget (DEB) theory, are gathering a great deal of interest as a tool to predict the effects of realistically variable exposure to toxicants over time on an individual animal. Here we use aquatic ecological risk assessment (ERA) as the context for a review of the different model variants within DEB and the closely related DEBkiss theory (incl. reserves, ageing, size & maturity, starvation). We propose a coherent and unifying naming scheme for all current major DEB variants, explore the implications of each model's underlying assumptions in terms of its capability and complexity and analyse differences between the models (endpoints, mathematical differences, physiological modes of action). The results imply a hierarchy of model complexity which could be used to guide the implementation of simplified model variants. We provide a decision tree to support matching the simplest suitable model to a given research or regulatory question. We detail which new insights can be gained by using DEB in toxicokinetic-toxicodynamic modelling, both generally and for the specific example of ERA, and highlight open questions. Specifically, we outline a moving time window approach to assess time-variable exposure concentrations and discuss how to account for cross-generational exposure. Where possible, we suggest valuable topics for experimental and theoretical research.

Integrating QSAR models predicting acute contact toxicity and mode of action profiling in honey bees (A. mellifera): Data curation using open source databases, performance testing and validation

Honey bees (Apis mellifera) provide key ecosystem services as pollinators bridging agriculture, the food chain and ecological communities, thereby ensuring food production and security. Ecological risk assessment of single Plant Protection Products (PPPs) requires an understanding of the exposure and toxicity. In silico tools such as QSAR models can play a major role for the prediction of structural, physico-chemical and pharmacokinetic properties of chemicals as well as toxicity of single and multiple chemicals. Here, the first integrative honey bee QSAR model has been developed for PPPs using EFSA's OpenFoodTox, US-EPA ECOTOX and Pesticide Properties DataBase i) to predict acute contact toxicity (LD₅₀) and ii) to profile the Mode of Action (MoA) of pesticides active substances. Three different classification-based and four regression-based models were developed and tested for their performance, thus identifying two models providing the most reliable predictions based on k-NN algorithm. The two-category QSAR model (toxic/non-toxic; n = 411) was validated using sensitivity (=0.93), specificity (=0.85), balanced accuracy (=0.90), and Matthews correlation coefficient (MCC = 0.78) as statistical parameters. The regression-based model (n = 113) was validated for its reliability and robustness (R² = 0.74; MAE = 0.52). Current study proposes the MoA profiling for 113 pesticides active substances and the first harmonised MoA classification scheme for acute contact toxicity in honey bees, including LD_50s data points from three different databases. The classification allows to further define MoAs and the target site of PPPs active substances, thus enabling regulators and scientists to refine chemical grouping and toxicity extrapolations for single chemicals and component-based mixture risk assessment of multiple chemicals. Relevant future perspectives are briefly addressed to integrate MoA, adverse outcome pathways (AOPs) and toxicokinetic information for the refinement of single-chemical/combined toxicity predictions and risk estimates at different levels of biological organization in the bee health context.

Invertebrate responses to microplastic ingestion: Reviewing the role of the antioxidant system

Microplastic ingestion in invertebrates reduces somatic and reproductive growth. This could be caused by energy reserves being detracted from growth processes and redistributed to maintenance processes that preserve life. A potential sink for this diverted energy is the antioxidant system, which minimises oxidative damage and reinstates redox homeostasis following disturbances caused by exposure to pollution. Several microplastic studies have used genetic and molecular redox biomarkers to assess how microplastic ingestion affects the functioning of the antioxidant system. This systematic review synthesises the current understanding of redox biomarker responses in invertebrates that have ingested microplastics. We found that biomarker response information exists for only seven invertebrate taxa, and early life stages have received little scientific attention. The microplastics used by most studies were polystyrene (45% of studies), spherical (51% of studies), and were < 10 μm in diameter (31% of studies). We found multiple examples of microplastic ingestion posing an oxidative challenge to invertebrates, which required upregulation of antioxidant system components. However, the lack of systematic experiments prevented us from clearly identifying which characteristic of microplastics caused these responses. We identify several areas for consideration when investigating biomarker responses to microplastic ingestion and offer research priorities for future studies.

Towards a spatiotemporally explicit toxicokinetic-toxicodynamic model for earthworm toxicity

The aim of the environmental risk assessment of chemicals is the prevention of unacceptable adverse effects on the environment. Therefore, the risk assessment for in-soil organisms, such as earthworms, is based on two key elements: the exposure assessment and the effect assessment. In the current risk assessment scheme, these two elements are not linked. While for the exposure assessment, advanced exposure models can take the spatial and temporal scale of substances into account, the effect assessment in the lower tiers considers only a limited temporal and spatial variability. However, for soil organisms, such as earthworms, those scales play a significant role as species move through the soil in response to environmental factors. To overcome this gap, we propose a conceptual integration of pesticide exposure, ecology, and toxicological effects on earthworms using a modular modeling approach. An essential part of this modular approach is the environment module, which utilizes exposure models to provide spatially and temporally explicit information on environmental variables (e.g., temperature, moisture, organic matter content) and chemical concentrations. The behavior module uses this information and simulates the feeding and movement of different earthworm species using a trait-based approach. The resulting exposure can be processed by a toxicokinetic-toxicodynamic (TKTD) module. TKTD models are particularly suitable to make effect predictions for time-variable exposure situations as they include the processes of uptake, elimination, internal distribution, and biotransformation of chemicals and link the internal concentration to an effect at the organism level. The population module incorporates existing population models of different earthworm species. The modular approach is illustrated using a case study with an insecticide. Our results emphasize that using a modular model approach will facilitate the integration of exposure and effects and thus enhance the risk assessment of soil organisms.

An open source physiologically based kinetic model for the chicken (Gallus gallus domesticus): Calibration and validation for the prediction residues in tissues and eggs

Xenobiotics from anthropogenic and natural origin enter animal feed and human food as regulated compounds, environmental contaminants or as part of components of the diet. After dietary exposure, a chemical is absorbed and distributed systematically to a range of organs and tissues, metabolised, and excreted. Physiologically based kinetic (PBK) models have been developed to estimate internal concentrations from external doses. In this study, a generic multi-compartment PBK model was developed for chicken. The PBK model was implemented for seven compounds (with log K_ow range -1.37-6.2) to quantitatively link external dose and internal dose for risk assessment of chemicals. Global sensitivity analysis was performed for a hydrophilic and a lipophilic compound to identify the most sensitive parameters in the PBK model. Model predictions were compared to measured data according to dataset-specific exposure scenarios. Globally, 71% of the model predictions were within a 3-fold change of the measured data for chicken and only 7% of the PBK predictions were outside a 10-fold change. While most model input parameters still rely on in vivo experiments, in vitro data were also used as model input to predict internal concentration of the coccidiostat monensin. Future developments of generic PBK models in chicken and other species of relevance to animal health risk assessment are discussed.

Ecological and human health risk assessment

The literature review presented in this paper covers the risk assessment process that is important to human health as well as the health of ecology in the form of receptors. One of the important objectives of present review is to provide summary of the scientific studies published in the year 2018. The review starts with literature published on the assessment of health risks, which are valuable to human and ecology. Most of the literature in the entire article focuses on techniques used for the analysis of scientific data and methods. In addition, review also highlights data interpretation, uncertainty, policy, and regulatory guidance associated with the management of human and ecological risks. Particularly, the review on the risk assessment related to human health and ecology is divided into two main sections. These sections provide broad state of knowledge on the risk assessment process used to health of human and ecological systems focused on investigation of polluted sites, techniques of remediation, and tools required for natural resource management.

Building and Applying Quantitative Adverse Outcome Pathway Models for Chemical Hazard and Risk Assessment

An important goal in toxicology is the development of new ways to increase the speed, accuracy, and applicability of chemical hazard and risk assessment approaches. A promising route is the integration of in vitro assays with biological pathway information. We examined how the adverse outcome pathway (AOP) framework can be used to develop pathway-based quantitative models useful for regulatory chemical safety assessment. By using AOPs as initial conceptual models and the AOP knowledge base as a source of data on key event relationships, different methods can be applied to develop computational quantitative AOP models (qAOPs) relevant for decision making. A qAOP model may not necessarily have the same structure as the AOP it is based on. Useful AOP modeling methods range from statistical, Bayesian networks, regression, and ordinary differential equations to individual-based models and should be chosen according to the questions being asked and the data available. We discuss the need for toxicokinetic models to provide linkages between exposure and qAOPs, to extrapolate from in vitro to in vivo, and to extrapolate across species. Finally, we identify best practices for modeling and model building and the necessity for transparent and comprehensive documentation to gain confidence in the use of qAOP models and ultimately their use in regulatory applications. Environ Toxicol Chem 2019;38:1850-1865. © 2019 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals, Inc. on behalf of SETAC.

Synergistic Effects of Climate Change and Marine Pollution: An Overlooked Interaction in Coastal and Estuarine Areas

Coastal areas have been increasingly affected by human activities, marine pollution and climate change are among the most important pressures affecting these environments. Human-induced pressures occur in a cumulative way and generate additive, antagonistic or synergistic effects. Knowledge on synergistic effects is crucial to coastal zone management, since they may imply a change in human uses of these systems, as well as dedicated action plans in order to reduce hazards and environmental risks. In this work, we provide an overview of the available literature on synergistic effects between climate change and chemical pollution, and discuss current knowledge, methodological approaches, and research gaps and needs. Interactions between these two pressures may be climate change dominant (climate change leads to an increase in contaminant exposure or toxicity) or contaminant-dominant (chemical exposure leads to an increase in climate change susceptibility), but the mechanistic drivers of such processes are not well known. Results from a few meta-analyses studies and reviews showed that synergistic interactions tend to be more frequent compared to additive and antagonistic ones. However, most of the studies are individual-based and assess the cumulative effects of a few contaminants individually in laboratory settings together with few climate variables, particularly temperature and pH. Nevertheless, a wide diversity of contaminants have already been individually tested, spanning from metals, persistent organic pollutants and, more recently, emergent pollutants. Population and community based approaches are less frequent but have generated very interesting and more holistic perspectives. Methodological approaches are quite diverse, from laboratory studies to mesocosm and field studies, or based on statistical or modelling tools, each with their own potential and limitations. More holistic comparisons integrating several pressures and their combinations and a multitude of habitats, taxa, life-stages, among others, are needed, as well as insights from meta-analyses and systematic reviews.

Caenorhabditis elegans as a complete model organism for biosafety assessments of nanoparticles

The number of biosafety evaluation studies of nanoparticles (NPs) using different biological models is increasing with the rapid development of nanotechnology. Thus far, nematode Caenorhabditis elegans (C. elegans), as a complete model organism, has become an important in vivo alternative assay system to assess the risk of NPs, especially at the environmental level. According to results of qualitative and quantitative analyses, it can be concluded that studies of nanoscientific research using C. elegans is persistently growing. However, the comprehensive conclusion and analysis of toxic effects of NPs in C. elegans are limited and chaotic. This review focused on the effects, especially sublethal ones, induced by NPs in C. elegans, including the development, intestinal function, immune response, neuronal function, and reproduction, as well as the underlying mechanisms of NPs causing these effects, including oxidative stress and alterations of several signaling pathways. Furthermore, we presented some factors that influence the toxic effects of NPs in C. elegans. The advantages and limitations of using nematodes in the nanotoxicology study were also discussed. Finally, we predicted that the application of C. elegans to assess long-term impacts of metal oxide NPs in the ecosystem would become a vital part of the nanoscientific research field, which provided an insight for further study.

Guidance on harmonised methodologies for human health, animal health and ecological risk assessment of combined exposure to multiple chemicals

This Guidance document describes harmonised risk assessment methodologies for combined exposure to multiple chemicals for all relevant areas within EFSA's remit, i.e. human health, animal health and ecological areas. First, a short review of the key terms, scientific basis for combined exposure risk assessment and approaches to assessing (eco)toxicology is given, including existing frameworks for these risk assessments. This background was evaluated, resulting in a harmonised framework for risk assessment of combined exposure to multiple chemicals. The framework is based on the risk assessment steps (problem formulation, exposure assessment, hazard identification and characterisation, and risk characterisation including uncertainty analysis), with tiered and stepwise approaches for both whole mixture approaches and component‐based approaches. Specific considerations are given to component‐based approaches including the grouping of chemicals into common assessment groups, the use of dose addition as a default assumption, approaches to integrate evidence of interactions and the refinement of assessment groups. Case studies are annexed in this guidance document to explore the feasibility and spectrum of applications of the proposed methods and approaches for human and animal health and ecological risk assessment. The Scientific Committee considers that this Guidance is fit for purpose for risk assessments of combined exposure to multiple chemicals and should be applied in all relevant areas of EFSA's work. Future work and research are recommended.

This publication is linked to the following EFSA Supporting Publications article: http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2019.EN-1589/full, http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2019.EN-1602/full

How do we take the pulse of an aquatic ecosystem? Current and historical approaches to measuring ecosystem integrity

Global environmental monitoring has indicated that the structure and function of some aquatic ecosystems has been significantly altered by human activities. There are many potential causes for these changes; however, one major concern is the increasing release of anthropogenic contaminants into aquatic environments. Although toxicological responses of individual organisms are typically well characterized, few studies have focused on characterizing toxicity at the ecosystem level. In fact, because of their scale and complexity, changes in ecosystem integrity are rarely considered in assessments of risks to ecosystems. This work attempts to move the conversation forward by defining integrity of ecosystems, reviewing current and historical approaches to measuring ecosystem integrity status (e.g., structural and functional measurements), and highlighting methods that could significantly contribute to the field of ecosystem toxicology (e.g., keystone species, environmental energetics, ecotoxicological modeling, and adverse outcome pathways [AOPs]). Through a critical analysis of current and historical methodologies, the present study offers a comprehensive, conceptual framework for the assessment of risks of contaminant exposure for whole ecosystems and proposes steps to facilitate better diagnoses of the integrity of aquatic systems. Environ Toxicol Chem 2019;38:289-301. © 2018 SETAC.

Calibration of parameters in Dynamic Energy Budget models using Direct-Search methods

Dynamic Energy Budget (DEB) theory aims to capture the quantitative aspects of metabolism at the individual level, for all species. The parametrization of a DEB model is based on information obtained through the observation of natural populations and experimental research. Currently the DEB toolbox estimates these parameters using the Nelder-Mead Simplex method, a derivative-free direct-search method. However, this procedure presents some limitations regarding convergence and how to address constraints. Framed in the calibration of parameters in DEB theory, this work presents a numerical comparison between the Nelder-Mead Simplex method and the SID-PSM algorithm, a Directional Direct-Search method for which convergence can be established both for unconstrained and constrained problems. A hybrid version of the two methods, named as Simplex Directional Direct-Search, provides a robust and efficient algorithm, able to solve the constrained optimization problems resulting from the parametrization of the biological models.

Current EU research activities on combined exposure to multiple chemicals

Humans and wildlife are exposed to an intractably large number of different combinations of chemicals via food, water, air, consumer products, and other media and sources. This raises concerns about their impact on public and environmental health. The risk assessment of chemicals for regulatory purposes mainly relies on the assessment of individual chemicals. If exposure to multiple chemicals is considered in a legislative framework, it is usually limited to chemicals falling within this framework and co-exposure to chemicals that are covered by a different regulatory framework is often neglected. Methodologies and guidance for assessing risks from combined exposure to multiple chemicals have been developed for different regulatory sectors, however, a harmonised, consistent approach for performing mixture risk assessments and management across different regulatory sectors is lacking. At the time of this publication, several EU research projects are running, funded by the current European Research and Innovation Programme Horizon 2020 or the Seventh Framework Programme. They aim at addressing knowledge gaps and developing methodologies to better assess chemical mixtures, by generating and making available internal and external exposure data, developing models for exposure assessment, developing tools for in silico and in vitro effect assessment to be applied in a tiered framework and for grouping of chemicals, as well as developing joint epidemiological-toxicological approaches for mixture risk assessment and for prioritising mixtures of concern. The projects EDC-MixRisk, EuroMix, EUToxRisk, HBM4EU and SOLUTIONS have started an exchange between the consortia, European Commission Services and EU Agencies, in order to identify where new methodologies have become available and where remaining gaps need to be further addressed. This paper maps how the different projects contribute to the data needs and assessment methodologies and identifies remaining challenges to be further addressed for the assessment of chemical mixtures.

Scientific Opinion on the state of the art of Toxicokinetic/Toxicodynamic (TKTD) effect models for regulatory risk assessment of pesticides for aquatic organisms

Following a request from EFSA, the Panel on Plant Protection Products and their Residues (PPR) developed an opinion on the state of the art of Toxicokinetic/Toxicodynamic (TKTD) models and their use in prospective environmental risk assessment (ERA) for pesticides and aquatic organisms. TKTD models are species- and compound-specific and can be used to predict (sub)lethal effects of pesticides under untested (time-variable) exposure conditions. Three different types of TKTD models are described, viz., (i) the 'General Unified Threshold models of Survival' (GUTS), (ii) those based on the Dynamic Energy Budget theory (DEBtox models), and (iii) models for primary producers. All these TKTD models follow the principle that the processes influencing internal exposure of an organism, (TK), are separated from the processes that lead to damage and effects/mortality (TD). GUTS models can be used to predict survival rate under untested exposure conditions. DEBtox models explore the effects on growth and reproduction of toxicants over time, even over the entire life cycle. TKTD model for primary producers and pesticides have been developed for algae, Lemna and Myriophyllum. For all TKTD model calibration, both toxicity data on standard test species and/or additional species can be used. For validation, substance and species-specific data sets from independent refined-exposure experiments are required. Based on the current state of the art (e.g. lack of documented and evaluated examples), the DEBtox modelling approach is currently limited to research applications. However, its great potential for future use in prospective ERA for pesticides is recognised. The GUTS model and the Lemna model are considered ready to be used in risk assessment.