Can TK-TD modelling bridge the gap between in vitro and in vivo mammalian toxicity data?

Repeated dietary dose testing is used to assess longer term toxicity of chemicals, such as pesticides, to mammals. However, the internal pesticide concentration varies significantly as feeding rate relative to body size fluctuates over time. Toxicokinetic-toxicodynamic (TK-TD) models can estimate internal toxicant concentration over time and link this directly to observed effects on endpoints such as the growth rate of laboratory rats. Using TK-TD models it is therefore possible to predict the effects that would result from a constant internal concentration of a pesticide. This presents the possibility of comparison with data from in vitro experiments, potentially facilitating quantitative in vitro to in vivo extrapolation (QIVIVE). We used in vivo TK-TD models to identify relevant internal concentrations and then estimated the experimental conditions required to replicate these in cultured cells, using in vitro TK models. Cell population growth was measured, with a view to extrapolating through time and comparing effect sizes with in vivo predictions. However, observed cell proliferation was not significantly affected by the tested concentrations of any of the five pesticides in this study and so extrapolation was not possible. In light of this negative result, we highlight areas for future work towards QIVIVE of graded sublethal effects in mammals. The most pressing objective is improving the accuracy of in vivo TK predictions, which could be achieved with dietary dosing in TK studies.

Environmental Risk Assessment of Time-Variable Toxicant Exposure with Toxicokinetic-Toxicodynamic Modeling of Sublethal Endpoints and Moving Time Windows: A Case Study with Ceriodaphnia dubia.

Toxicokinetic-toxicodynamic (TKTD) modeling has received increasing attention in terms of the regulatory environmental risk assessment of chemicals. This type of mechanistic model can integrate all available data from individual-level bioassays into a single framework and enable refined risk assessments by extrapolating from laboratory results to time-variable exposure scenarios, based, for instance, on surface water exposure modeling (e.g., FOCUS). Dynamic energy budget (DEB) models coupled with TKTD modules (DEB-TKTD) constitute the leading approach to assess and predict sublethal effects of chemicals on individual organisms. However, thorough case studies are rare. We provide a state-of-the-art example with the standard aquatic test species Ceriodaphnia dubia and the fungicide azoxystrobin, including all steps, from bespoke laboratory toxicity tests to model calibration and validation, through to environmental risk assessment. Following the framework proposed in the European Food Safety Authority Scientific Opinion from 2018, we designed bespoke good laboratory practice-compliant laboratory toxicity studies based on test guideline 211 of the Organisation for Economic Co-operation and Development and then identified robust parameter values from those data for all relevant model parameters through model calibration. The DEB-TKTD model, DEBtox2019, then informed the design of the validation experiment. Once validated, the model was used to perform predictions for a time-variable exposure scenario generated by FOCUS. A moving time-window approach was used to perform the environmental risk assessment. This assessment method reduces uncertainty in the risk assessment while maintaining consistency with the traditional measures of risk. Environ Toxicol Chem 2024;00:1-13. © 2024 Syngenta Crop Protection AG. ibacon GmbH and The Author(s). Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

A User-Friendly Kinetic Model Incorporating Regression Models for Estimating Pesticide Accumulation in Diverse Earthworm Species Across Varied Soils.

Existing models for estimating pesticide bioconcentration in earthworms exhibit limited applicability across different chemicals, soils and species which restricts their potential as an alternative, intermediate tier for risk assessment. We used experimental data from uptake and elimination studies using three earthworm species (Lumbricus terrestris, Aporrectodea caliginosa, Eisenia fetida), five pesticides (log Kow 1.69-6.63) and five soils (organic matter content = 0.972-39.9 wt %) to produce a first-order kinetic accumulation model. Model applicability was evaluated against a data set of 402 internal earthworm concentrations reported from the literature including chemical and soil properties outside the data range used to produce the model. Our models accurately predict body load using either porewater or bulk soil concentrations, with at least 93.5 and 84.3% of body load predictions within a factor of 10 and 5 of corresponding observed values, respectively. This suggests that there is no need to distinguish between porewater and soil exposure routes or to consider different uptake and elimination pathways when predicting earthworm bioconcentration. Our new model not only outperformed existing models in characterizing earthworm exposure to pesticides in soil, but it could also be integrated with models that account for earthworm movement and fluctuating soil pesticide concentrations due to degradation and transport.

Fate of synthetic chemicals in the agronomic insect pest Spodoptera littoralis: experimental feeding-contact assay and toxicokinetic model.

Insecticides prevent or reduce insect crop damage, maintaining crop quality and quantity. Physiological traits, such as an insect's feeding behavior, influence the way insecticides are absorbed and processed in the body (toxicokinetics), which can be exploited to improve species selectivity. To fully understand the uptake of insecticides, it is essential to study their total uptake and toxicokinetics independent of their toxic effects on insects. We studied the toxicokinetics (TK) of insecticidally inactive test compounds incorporating agro-like structural motifs in larvae of the Egyptian cotton leafworm (Spodoptera littoralis, Lepidoptera), and their distribution across all biological matrices, using laboratory experiments and modeling. We measured Spodoptera larval behavior and temporal changes of whole-body concentrations of test compounds during feeding on treated soybean leaf disks and throughout a subsequent depuration period. Differences in the distribution of the total quantities of compounds were found between the biological matrices leaf, larva, and feces. Rate constants for uptake and elimination of test compounds were derived by calibrating a toxicokinetic model to the whole-body concentrations. Uptake and elimination rate constants depended on the physicochemical properties of the test compounds. Increasing hydrophobicity increased the bioaccumulation potential of test compounds. Incomplete quantities in larval matrices indicated that some compounds may undergo biotransformation. As fecal excretion was a major elimination pathway, the variable time of release and number of feces pellets led to a high variability in the body burden. We provide quantitative models to predict the toxicokinetics and bioaccumulation potential of inactive insecticide analogs (parent compounds) in Spodoptera.

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.

Earthworm lipid content and size help account for differences in pesticide bioconcentration between species.

The uptake and elimination kinetics of pesticides from soil to earthworms are important in characterising the risk of pesticides to soil organisms and the risk from secondary poisoning. However, the understanding of the relative importance of chemical, soil, and species differences in determining pesticide bioconcentration into earthworms is limited. Furthermore, there is insufficient independent data in the literature to fully evaluate existing predictive bioconcentration models. We conducted kinetic uptake and elimination experiments for three contrasting earthworm species (Lumbricus terrestris, Aporrectodea caliginosa, Eisenia fetida) in five soils using a mixture of five pesticides (log Kow 1.69 - 6.63). Bioconcentration increased with pesticide hydrophobicity and decreased with soil organic matter. Bioconcentration factors were comparable between earthworm species for hydrophilic pesticides due to the similar water content of earthworm species. Inter-species variations in bioconcentration of hydrophobic pesticides were primarily accounted for by earthworm lipid content and specific surface area (SSA). Existing bioconcentration models either failed to perform well across earthworm species and for more hydrophilic compounds (log Kow < 2) or were not parameterised for a wide range of compounds and earthworm species. Refined models should incorporate earthworm properties (lipid content and SSA) to account for inter-species differences in pesticide uptake from soil.

Characterization of patterns and variability in the dynamics of outdoor aquatic mesocosms: exploring the capabilities and challenges in data supporting aquatic system models.

Aquatic mesocosms are complex test systems used within regulatory risk assessment of plant protection products. These model ecosystems allow researchers to capture interactions of multiple species under realistic environmental conditions. They enable assessment of direct and indirect effects of stressors at all trophic levels (i.e., from primary producers to secondary consumers) and impacts on ecosystem functions. Due to the limited ability to test the multitude of potential exposure scenarios, cross-linking aquatic mesocosm studies with virtual mesocosms, i.e., aquatic system models (ASMs), can serve to meet the demand for more environmental realism and ecological relevance in risk assessment. In this study, full control data sets from seven aquatic mesocosm studies conducted at a single test facility under GLP were analysed graphically and using descriptive statistics. Thereby, not only a comprehensive data base but also an insight into the species present, their dynamics over time, and variability in unchallenged mesocosms was observed. While consistency in dynamics could be discerned for physical and chemical parameters, variability was evident for several biological endpoints. This variability points to amplification of small differences over time as well as to stochastic processes. The outline of existing gaps and uncertainties in data leads to the estimation of what can be expected to be captured and predicted by ASMs.

A Framework for Algae Modeling in Regulatory Risk Assessment.

The use of toxicokinetic-toxicodynamic (TKTD) modeling in regulatory risk assessment of plant protection products is increasingly popular, especially since the 2018 European Food Safety Authority (EFSA) opinion on TKTD modeling announced that several established models are ready for use in risk assessment. With careful adherence to the guidelines laid out by EFSA, we present a stepwise approach to validation and use of the Simple Algae Model Extended (SAM-X) for regulatory submission in Tier 2C. We demonstrate how the use of moving time windows across time-variable exposure profiles can generate thousands of virtual laboratory mimic simulations that seamlessly predict the effects of time-variable exposures across a full exposure profile while maintaining the laboratory conditions of the standard Organisation for Economic Co-operation and Development (OECD) growth inhibition test. Thus, every virtual laboratory test has a duration of 72 h, with OECD medium and constant light and temperature conditions. The only deviation from the standard test setup is the replacement of constant exposure conditions for time-variable concentrations. The present study demonstrates that for simulation of 72-h toxicity tests, the nutrient dynamics in the SAM-X model are not required, and we propose the alternative use of a simplified model version. For risk assessment, in accordance with the EFSA guidelines we use a median exposure profile of 10 as a threshold, meaning that if a time window within the exposure profile causes 50% growth inhibition when magnified by a factor of 10, the threshold will have been exceeded. We present a simplified example for chlorotoluron and isoproturon. The present case study brings to life our proposed framework for TKTD modeling of algae to establish whether a given exposure can be considered to be of low risk. Environ Toxicol Chem 2023;42:1823-1838. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

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.

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.

How to analyse and account for interactions in mixture toxicity with toxicokinetic-toxicodynamic models.

The assessment of chemical mixture toxicity is one of the major challenges in ecotoxicology. Chemicals can interact, leading to more or less effects than expected, commonly named synergism and antagonism respectively. The classic ad hoc approach for the assessment of mixture effects is based on dose-response curves at a single time point, and is limited to identifying a mixture interaction but cannot provide predictions for untested exposure durations, nor for scenarios where exposure varies in time. We here propose a new approach using toxicokinetic-toxicodynamic modelling: The General Unified Threshold model of Survival (GUTS) framework, recently extended for mixture toxicity assessment. We designed a dedicated mechanistic interaction module coupled with the GUTS mixture model to i) identify interactions, ii) test hypotheses to identify which chemical is likely responsible for the interaction, and finally iii) simulate and predict the effect of synergistic and antagonistic mixtures. We tested the modelling approach experimentally with two species (Enchytraeus crypticus and Mamestra brassicae) exposed to different potentially synergistic mixtures (composed of: prochloraz, imidacloprid, cypermethrin, azoxystrobin, chlorothalonil, and chlorpyrifos). Furthermore, we also tested the model with previously published experimental data on two other species (Bombus terrestris and Daphnia magna) exposed to pesticide mixtures (clothianidin, propiconazole, dimethoate, imidacloprid and thiacloprid) found to be synergistic or antagonistic with the classic approach. The results showed an accurate simulation of synergistic and antagonistic effects for the different tested species and mixtures. This modelling approach can identify interactions accounting for the entire time of exposure, and not only at one time point as in the classic approach, and provides predictions of the mixture effect for untested mixture exposure scenarios, including those with time-variable mixture composition.

Interactive effects of multiple stressors vary with consumer interactions, stressor dynamics and magnitude.

Predicting the impacts of multiple stressors is important for informing ecosystem management but is impeded by a lack of a general framework for predicting whether stressors interact synergistically, additively or antagonistically. Here, we use process-based models to study how interactions generalise across three levels of biological organisation (physiological, population and consumer-resource) for a two-stressor experiment on a seagrass model system. We found that the same underlying processes could result in synergistic, additive or antagonistic interactions, with interaction type depending on initial conditions, experiment duration, stressor dynamics and consumer presence. Our results help explain why meta-analyses of multiple stressor experimental results have struggled to identify predictors of consistently non-additive interactions in the natural environment. Experiments run over extended temporal scales, with treatments across gradients of stressor magnitude, are needed to identify the processes that underpin how stressors interact and provide useful predictions to management.

Fish Species Sensitivity Ranking Depends on Pesticide Exposure Profiles.

In the regulatory environmental risk assessment of plant protection products, the exposure tested in standard toxicity tests assumes simple exposure dynamics, such as constant exposure at the first stage of testing. However, environmental exposure can be highly dynamic. A species response to exposure is governed by toxicokinetics (TK) and toxicodynamics (TD). Therefore, it can be expected that the sensitivity of a species to a substance is dependent on the interplay of TKTD processes with the dynamics of the exposure. We investigated whether exposure dynamics affects species sensitivity of five fish species and if their sensitivity rankings differ among exposure profiles. We analyzed individual survival under projected surface water exposure to benzovindiflupyr. For this purpose, we calibrated compound- and species-specific reduced general unified threshold models of survival (GUTS-RED) models from standard laboratory toxicity data with the assumptions of stochastic death and individual tolerance. Using the calibrated models, we generated species sensitivity distributions based on median lethal profile multiplication factors for three characteristic exposure profiles. The analysis was performed using different GUTS-RED implementations: openGUTS (MATLAB® and Windows® versions) and the R package morse. The sensitivity rankings of the fish species changed as a function of exposure profile. For a multiple-peak scenario, rainbow trout was the most sensitive species. For a single peak followed by a slow concentration decline the most sensitive species was the fathead minnow (GUTS-RED-stochastic death) or the common carp (GUTS-RED-individual tolerance). Our results suggest that a single most sensitive species cannot be defined for all situations, all exposure profiles, and both GUTS-RED variants. Environ Toxicol Chem 2022;41:1732-1741. © 2022 Syngenta. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

The application and limitations of exposure multiplication factors in sublethal effect modelling.

Thanks to growing interest and research in the field, toxicokinetic-toxicodynamic (TKTD) models are close to realising their potential in environmental risk assessment (ERA) of chemicals such as plant protection products. A fundamental application is to find a multiplicative scale factor which-when applied to an exposure profile-results in some specified effect relative to a control. The approach is similar to applying assessment factors to experimental results, common in regulatory frameworks. It also relies on the same core assumption: that increasing the scaling always produces more extreme effects. Unlike experimental approaches, TKTD models offer an opportunity to interrogate this assumption in a mathematically rigorous manner. For four well-known TKTD models we seek to prove that the approach guarantees a unique scale factor for any percentage effect. Somewhat surprisingly, certain model configurations may have multiple scale factors which result in the same percentage effect. These cases require a more cautious regulatory approach and generate open biological and mathematical questions. We provide examples of the violations and suggest how to deal with them. Mathematical proofs provide the strongest possible backing for TKTD modelling approaches in ERA, since the applicability of the models can be determined exactly.

Modelling the effects of variability in feeding rate on growth - a vital step for DEB-TKTD modelling.

A major limitation of dietary toxicity studies on rodents is that food consumption often differs between treatments. The control treatment serves as a reference of how animals would have grown if not for the toxicant in their diet, but this comparison unavoidably conflates the effects of toxicity and feeding rate on body weight over time. A key advantage of toxicity models based on dynamic energy budget theory (DEB) is that chemical stress and food consumption are separate model inputs, so their effects on growth rate can be separated. To reduce data requirements, DEB convention is to derive a simplified feeding input, f, from food availability; its value ranges from zero (starvation) to one (food available ad libitum). Observed food consumption in dietary toxicity studies shows that, even in the control treatment, rats limit their food consumption, contradicting DEB assumptions regarding feeding rate. Relatively little work has focused on addressing this mismatch, but accurately modelling the effects of food intake on growth rate is essential for the effects of toxicity to be isolated. This can provide greater insight into the results of chronic toxicity studies and allows accurate extrapolation of toxic effects from laboratory data. Here we trial a new method for calculating f, based on the observed relationships between food consumption and body size in laboratory rats. We compare model results with those of the conventional DEB method and a previous effort to calculate f using observed food consumption data. Our results showed that the new method improved model accuracy while modelled reserve dynamics closely followed observed body fat percentage over time. The new method assumes that digestive efficiency increases with body size. Verifying this relationship through data collection would strengthen the basis of DEB theory and support the case for its use in ecological risk assessment.

Considerations for using reproduction data in toxicokinetic-toxicodynamic modeling.

Toxicokinetic-toxicodynamic (TKTD) modeling is essential to make sense of the time dependence of toxic effects, and to interpret and predict consequences of time-varying exposure. These advantages have been recognized in the regulatory arena, especially for environmental risk assessment of pesticides, where time-varying exposure is the norm. We critically evaluate the link between the modeled variables in TKTD models and the observations from laboratory ecotoxicity tests. For the endpoint reproduction, this link is far from trivial. The relevant TKTD models for sublethal effects are based on dynamic energy budget (DEB) theory, which specifies a continuous investment flux into reproduction. In contrast, experimental tests score egg or offspring release by the mother. The link between model and data is particularly troublesome when a species reproduces in discrete clutches and, even more so, when eggs are incubated in the mother's brood pouch (and release of neonates is scored in the test). This situation is quite common among aquatic invertebrates (e.g., cladocerans, amphipods, mysids), including many popular test species. In this discussion paper, we treat these and other issues with reproduction data, reflect on their potential impact on DEB-TKTD analysis, and provide preliminary recommendations to correct them. Both modelers and users of model results need to be aware of these complications, as ignoring them could easily lead to unnecessary failure of DEB-TKTD models during calibration, or when validating them against independent data for other exposure scenarios. Integr Environ Assess Manag 2022;18:479-487. © 2021 SETAC.

Predicting Mixture Effects over Time with Toxicokinetic-Toxicodynamic Models (GUTS): Assumptions, Experimental Testing, and Predictive Power.

Current methods to assess the impact of chemical mixtures on organisms ignore the temporal dimension. The General Unified Threshold model for Survival (GUTS) provides a framework for deriving toxicokinetic-toxicodynamic (TKTD) models, which account for effects of toxicant exposure on survival in time. Starting from the classic assumptions of independent action and concentration addition, we derive equations for the GUTS reduced (GUTS-RED) model corresponding to these mixture toxicity concepts and go on to demonstrate their application. Using experimental binary mixture studies with Enchytraeus crypticus and previously published data for Daphnia magna and Apis mellifera, we assessed the predictive power of the extended GUTS-RED framework for mixture assessment. The extended models accurately predicted the mixture effect. The GUTS parameters on single exposure data, mixture model calibration, and predictive power analyses on mixture exposure data offer novel diagnostic tools to inform on the chemical mode of action, specifically whether a similar or dissimilar form of damage is caused by mixture components. Finally, observed deviations from model predictions can identify interactions, e.g., synergism or antagonism, between chemicals in the mixture, which are not accounted for by the models. TKTD models, such as GUTS-RED, thus offer a framework to implement new mechanistic knowledge in mixture hazard assessments.

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

Earthworms are important ecosystem engineers, and assessment of the risk of plant protection products toward them is part of the European environmental risk assessment (ERA). In the current ERA scheme, exposure and effects are represented simplistically and are not well integrated, resulting in uncertainty when the results are applied to ecosystems. Modeling offers a powerful tool to integrate the effects observed in lower tier laboratory studies with the environmental conditions under which exposure is expected in the field. This paper provides a summary of the (In)Field Organism Risk modEling by coupling Soil Exposure and Effect (FORESEE) Workshop held 28-30 January 2020 in Düsseldorf, Germany. This workshop focused on toxicokinetic-toxicodynamic (TKTD) and population modeling of earthworms in the context of ERA. The goal was to bring together scientists from different stakeholder groups to discuss the current state of soil invertebrate modeling and to explore how earthworm modeling could be applied to risk assessments, in particular how the different model outputs can be used in the tiered ERA approach. In support of these goals, the workshop aimed at addressing the requirements and concerns of the different stakeholder groups to support further model development. The modeling approach included 4 submodules to cover the most relevant processes for earthworm risk assessment: environment, behavior (feeding, vertical movement), TKTD, and population. Four workgroups examined different aspects of the model with relevance for risk assessment, earthworm ecology, uptake routes, and cross-species extrapolation and model testing. Here, we present the perspectives of each workgroup and highlight how the collaborative effort of participants from multidisciplinary backgrounds helped to establish common ground. In addition, we provide a list of recommendations for how earthworm TKTD modeling could address some of the uncertainties in current risk assessments for plant protection products. Integr Environ Assess Manag 2021;17:352-363. © 2020 SETAC.

Bioenergetics modelling to analyse and predict the joint effects of multiple stressors: Meta-analysis and model corroboration.

Understanding the consequences of the combined effects of multiple stressors-including stress from man-made chemicals-is important for conservation management, the ecological risk assessment of chemicals, and many other ecological applications. Our current ability to predict and analyse the joint effects of multiple stressors is insufficient to make the prospective risk assessment of chemicals more ecologically relevant because we lack a full understanding of how organisms respond to stress factors alone and in combination. Here, we describe a Dynamic Energy Budget (DEB) based bioenergetics model that predicts the potential effects of single or multiple natural and chemical stressors on life history traits. We demonstrate the plausibility of the model using a meta-analysis of 128 existing studies on freshwater invertebrates. We then validate our model by comparing its predictions for a combination of three stressors (i.e. chemical, temperature, and food availability) with new, independent experimental data on life history traits in the daphnid Ceriodaphnia dubia. We found that the model predictions are in agreement with observed growth curves and reproductive traits. To the best of our knowledge, this is the first time that the combined effects of three stress factors on life history traits observed in laboratory studies have been predicted successfully in invertebrates. We suggest that a re-analysis of existing studies on multiple stressors within the modelling framework outlined here will provide a robust null model for identifying stressor interactions, and expect that a better understanding of the underlying mechanisms will arise from these new analyses. Bioenergetics modelling could be applied more broadly to support environmental management decision making.