New Insights to Compare and Choose TKTD Models for Survival Based on an Interlaboratory Study for Lymnaea stagnalis Exposed to Cd

Application of toxicokinetic-toxicodynamic models in the aquatic ecological risk assessment of metals: A review

The issue of toxic metal pollution is a considerable environmental concern owing to its complex nature, spatial and temporal variability, and susceptibility to environmental factors. Current water quality criteria and ecological risk assessments of metals are based on single-metal toxicity data from short-term, simplified indoor exposure conditions, ignoring the complexity of actual environmental conditions. This results in increased uncertainty in predicting toxic metal toxicity and risk assessment. Using appropriate bioavailability and effect modeling of metals is critical for establishing environmental quality standards and performing risk assessments for metals. Traditional dose-effect models are based on a static statistical relationship and fall short of revealing the bioavailability and effect processes of metals and do not effectively assess ecological impacts under complex exposure conditions. This paper summarizes the toxicokinetic-toxicodynamic (TK-TD) model, which is gaining interest in environmental and ecotoxicological research. The key concepts, and theories of its construction theories, are discussed and the application of the TK-TD model in toxicity prediction and risk assessment of different metals in the aquatic environment, and trends in the development of the TK-TD model are highlighted. The findings of our review prove that the TK-TD model can effectively predict toxic metal toxicity in real time and under complex exposure conditions in the future.

In Vitro Toxicity and Modeling Reveal Nanoplastic Effects on Marine Bivalves

Nanoplastics (NPs) represent a growing concern for global environmental health, particularly in marine ecosystems where they predominantly accumulate. The impact of NPs on marine benthic organisms, such as bivalves, raises critical questions regarding ecological integrity and food safety. Traditional methods for assessing NP toxicity are often limited by their time-intensive nature and ethical considerations. Herein, we explore the toxicological effects of NPs on the marine bivalve Ruditapes philippinarum, employing a combination of in vitro cellular assays and advanced modeling techniques. Results indicate a range of adverse effects at the organismal level, including growth inhibition (69.5-108%), oxidative stress, lipid peroxidation, and DNA damage in bivalves, following exposure to NPs at concentrations in the range of 1.6 × 109-1.6 × 1011 particles/mL (p/mL). Interestingly, the growth inhibition predicted by models (54.7-104%), based on in vitro cellular proliferation assays, shows strong agreement with the in vivo outcomes of NP exposure. Furthermore, we establish a clear correlation between cytotoxicity observed in vitro and the toxicological responses at the organismal level. Taken together, this work suggests that the integration of computational modeling with in vitro toxicity assays can predict the detrimental effects of NPs on bivalves, offering insightful references for assessing the environmental risk assessment of NPs in marine benthic ecosystems.

Mechanistic analysis of the sub chronic toxicity of La and Gd in Daphnia magna based on TKTD modelling and synchrotron X-ray fluorescence imaging

The release of lanthanides (Ln) into the environment has increased in recent decades due to their expanding applications in society. Studying their toxicity in aquatic ecosystems is urgent and challenging, with contradictory evidence presented in the literature. This study compared the biodistribution of La and Gd in Daphnia magna exposed to sub-chronic conditions and developed the first Toxicokinetic-Toxicodynamic (TKTD) model for these lanthanides with this model crustacean. D. magna were initially exposed for 7 days to concentrations close to the LC50 of La (2.10 mg L-1) and Gd (1.70 mg L-1). After exposure, half of the live daphnids were introduced in a clean media to allow depuration over 24 h, while the other organisms were directly prepared for synchrotron imaging measurements. Synchrotron X-ray fluorescence analysis revealed that metal distribution in the organisms was similar for both La and Gd, predominantly localized in the intestinal tract, even after the depuration process. These results indicate that ingested metal can adversely affect organisms under sub-chronic exposure conditions, highlighting the importance of using nominal concentrations as a more suitable indicator of metal bioavailability for risk assessment. The General Unified Threshold Model of Survival (GUTS) TKTD framework, in its reduced form (GUTS-RED), was developed for La and Gd using dissolved and nominal concentrations. D. magna were exposed for 7 days to concentrations from 0.5 to 5 mg L-1 of La or Gd and mortality monitored daily. The mechanistic model revealed a faster toxicokinetics for La than Gd and a higher toxicity for Gd than La in the organism. This study confirmed, despite similar chemical properties, the variation in both toxicity and toxicokinetics between these two metals.

hb or not hb - when and why accounting for background mortality in toxicological survival models matters?

Decisions in Environmental Risk Assessment (ERA) about impacts of chemical compounds on different species are based on critical effect indicators such as the 50% lethal concentration (LC₅₀). Regulatory documents recommend concentration-response (or concentration-effect) model fitting on standard toxicity test data to get LC₅₀ values. However, toxicokinetic-toxicodynamic (TKTD) models proved their efficiency to better exploit toxicity test data, at Tier-2 but also at Tier-1, delivering time-independent indicators. In particular, LC₅₀ values can be obtained from the reduced General Unified Threshold model of Survival (GUTS-RED) with both variants, Stochastic Death and Individual Tolerance, that include parameter h_b, the background mortality. Estimating h_b during the fitting process or not depends on studies and fitting habits, while it may strongly influence the other GUTS-RED parameters, and consequently the LC₅₀ estimate. We hypothesized that estimating h_b from all data in all replicates over time should provide more precise LC₅₀ estimates. We then explored how estimating h_b impacted: (i) GUTS-RED model parameters; (ii) goodness-of-fit criteria (fitting plot, posterior predictive check, parameter correlations); (iii) LC₅₀ accuracy and precision. We finally show that estimating h_b does not impact the LC₅₀ precision while providing more accurate and precise GUTS parameter estimates. Hence, estimating h_b would lead to a more protective ERA.

Modeling Time-Dependent Aquatic Toxicity of Hydrocarbons: Role of Organism Weight, Temperature, and Substance Hydrophobicity

Oil spill exposures are highly dynamic and are not comparable to laboratory exposures used in standard toxicity tests. Toxicokinetic-toxicodynamic (TKTD) models allow translation of effects observed in the laboratory to the field. To improve TKTD model calibration, new and previously published data from 148 tests were analyzed to estimate rates characterizing the time course of toxicity for 10 fish and 42 invertebrate species across 37 hydrocarbons. A key parameter in the TKTD model is the first-order rate that incorporates passive elimination, biotransformation, and damage repair processes. The results indicated that temperature (4-26 °C), organism size (0.0001-10 g), and substance log octanol-water partition coefficient (2-6) had limited influence on this parameter, which exhibited a 5th to 95th percentile range of 0.2-2.5 day-1 (median 0.7 day-1 ). A species sensitivity distribution approach is proposed to quantify the variability of this parameter across taxa, with further studies needed for aliphatic hydrocarbons and plant species. Study findings allow existing oil spill models to be refined to improve effect predictions. Environ Toxicol Chem 2022;41:3070-3083. © 2022 ExxonMobil Biomedical Science Inc. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

The effects of different temperatures in mercury toxicity to the terrestrial isopod Porcellionides pruinosus

Climate changes and metal contamination are pervasive stressors for soil ecosystems. Mercury (Hg), one of the most toxic metals, has been reported to interact with temperature. However, compared to aquatic biota, little is known about how temperature affects Hg toxicity and bioaccumulation to soil organisms. Here, toxicity and bioaccumulation experiments were replicated at 15 °C, 20 °C, and 25 °C to understand how sub-optimal temperatures affect the toxicokinetics and toxicodynamics of Hg via soil. Genotoxicity and energy reserves were also assessed to disclose potential trade-offs in life-history traits. Results underpin the complexity of temperature-Hg interactions. Survival was determined mainly by toxicokinetics, but toxicodynamics also played a significant role in defining survival probability during early stages. The processes determining survival probability were faster at 25 °C: General Unified Threshold of Survival (GUTS) model identified an earlier/steeper decline in survival, compared to 20 °C or 15 °C, but it also approached the threshold faster. Despite potentiation of Hg genotoxicity, temperature promoted faster detoxification, either increasing toxicokinetics rates or damage repair mechanisms. This metabolism-driven increase in detoxification led to higher depletion of energy reserves and likely triggered stress response pathways. This work emphasized the need for comprehensive experimental approaches that can integrate the multiple processes involved in temperature-metal interactions.

The toxicity of cadmium-copper mixtures on daphnids and microalgae analyzed using the Biotic Ligand Model

For the prediction of metal mixture ecotoxicity, the BLM approach is promising since it evaluates the amount of metals accumulated on the biotic ligand on the basis of water chemistry, i.e., species (major cations) competing with metals, and related toxicity. Based on previous work by Farley et al. (2015) (MMME research project), this study aimed at modeling toxicity of Cd:Cu mixtures (0:1 - 1:1 - 1:0 - 1:2 - 1:3 - 2:1 - 3:1 - 4:1 - 5:1 - 6:1) to the crustacean Daphnia magna(48-h immobilization tests) and the microalga Pseudokirchneriella subcapitata (72-h growth inhibition tests). The 2012 version of the USGS model was chosen, assuming additivity of effects and accumulation of metals on a single site. The assumption that EDTA could contribute to toxicity through metals complexing was also tested, and potential effects due to reduction of ions Ca²⁺ absorption by metals were considered. Modeling started with parameter values of Farley et al. (2015) and some of these parameters were adjusted to fit modeled data on observed data. The results show that toxicity can be correctly predicted for the microalgae and that the hypothesis of additivity is verified. For daphnids, the prediction was roughly correct, but taking into account CuEDTA led to more realistic parameter values close to that reported by Farley et al. (2015). However, It seems that, for daphnids responses, metals interact either antagonistically or synergistically depending on the Cu:Cd ratio. Furthermore, synergy could not be explained by additional effects linked to a reduction of Ca absorption since this reduction, mainly due to Cd, increased inversely to synergy. Finally, the USGS model applied to our data was able to predict Cu:Cd mixture toxicity to microalgae and daphnids, giving rise to estimated EC50s roughly reflecting EC50s calculated from observed toxicity.

Taking full advantage of modelling to better assess environmental risk due to xenobiotics-the all-in-one facility MOSAIC

In the European Union, more than 100,000 man-made chemical substances are awaiting an environmental risk assessment (ERA). Simultaneously, ERA of these chemicals has now entered a new era requiring determination of risks for physiologically diverse species exposed to several chemicals, often in mixtures. Additionally, recent recommendations from regulatory bodies underline a crucial need for the use of mechanistic effect models, allowing assessments that are not only ecologically relevant, but also more integrative, consistent and efficient. At the individual level, toxicokinetic-toxicodynamic (TKTD) models are particularly encouraged for the regulatory assessment of pesticide-related risks on aquatic organisms. In this paper, we first briefly present a classical dose-response model to showcase the on-line MOSAIC tool, which offers all necessary services in a turnkey web platform, whatever the type of data analyzed. Secondly, we focus on the necessity to account for the time-dimension of the exposure by illustrating how MOSAIC can support a robust calculation of bioaccumulation metrics. Finally, we show how MOSAIC can be of valuable help to fully complete the EFSA workflow regarding the use of TKTD models, especially with GUTS models, providing a user-friendly interface for calibrating, validating and predicting survival over time under any time-variable exposure scenario of interest. Our conclusion proposes a few lines of thought for an easier use of modelling in ERA.

Terrestrial organisms react differently to nano and non-nano Cu(OH)2 forms

The application of nanomaterials in agriculture is deemed as a promising strategy to increase the sector's sustainability. Nanopesticides are expected to improve solutions for pest/disease management, while reducing adverse effects to ecosystems, in accordance with the European Green Deal strategy. Hitherto, nanopesticide research has mostly focused on assessing effects to target species or crops, and less attention has been devoted to non-target soil species. In this study, we explored whether three copper hydroxide-based forms (nano and non-nano) show different toxicity and accumulation patterns in two terrestrial invertebrates, the isopod Porcellionides pruinosus and the mealworm larvae Tenebrio molitor. Toxicity and bioaccumulation experiments were performed and time-course toxicity and toxicokinetics analyzed and modelled. We found important differences in copper sensitivity and accumulation: T. molitor was more resistant and able to eliminate copper compared to P. pruinosus. The nanopesticide mostly elicited lower effects for both species. Slower toxicokinetics by the nanopesticide and lower bioaccumulative potential to P. pruinosus must be partially responsible for the differences in toxicity. Although no toxicokinetics differences were found in T. molitor between Cu forms, distinct negative effects on growth must reflect different energetic costs for copper regulation. Even though effects on toxicity and bioaccumulation are species-specific, copper hydroxide-based nanopesticides may cause less adverse effects to soil organisms than conventional products, being a good solution for reducing the environmental impact. CAPSULE: Copper hydroxide-based nanopesticide was less toxic and bioaccumulative to soil organisms than conventional copper hydroxide products, being a good solution for reducing the environmental impact.

Assessing the Relative Toxicity of Different Road Salts and Effect of Temperature on Salinity Toxicity: LCx Values versus No-Effect Concentration (NEC) Values

Freshwater biota are at risk globally from increasing salinity, including increases from deicing salts in cold regions. A variety of metrics of toxicity are used when estimating the toxicity of substances and comparing the toxicity between substances. However, the implications of using different metrics are not widely appreciated. Using the mayfly Colobruscoides giganteus (Ephemeroptera: Colobruscoidea), we compare the toxicity of seven different salts where toxicity was estimated using two metrics: (1) the no-effect concentrations (NEC) and (2) the lethal concentrations for 10, 25 and 50% of the test populations (LC_x). The LC_x values were estimated using two different models, the classic log-logistic model and the newer toxicokinetic-toxicodynamic (TKTD) model. The NEC and both types of LC_x values were estimated using Bayesian statistics. We also compared the toxicity of two salts (NaCl and CaCl₂) for C. giganteus at water temperatures of 4 °C, 7 °C and 15 °C using the same metrics of toxicity. Our motivation for using a mayfly to assess salinity toxicity was because mayflies are generally salt sensitive, are ecologically important and are common in Australian (sub-)alpine streams. The temperature ranges were chosen to mimic winter, spring and summer water temperatures for Australian (sub-)alpine streams. Considering 144-h classical LC_x values, we found toxicity differed between various salts, i.e., the lowest 144-h LC₅₀ (8 mS/cm) for a salt used by a ski resort was half that of the highest 144-h LC₅₀ from artificial marine salts and CaCl₂ applied to roads (16 mS/cm). The analytical grade NaCl (as shown by 144-h LC₅₀ value at 7 °C) was substantially more toxic (7.3 mS/cm) compared to analytical grade CaCl₂ (12.5 mS/cm). Yet for NEC values, there were comparably fewer differences in toxicity between salts and none between the same salts at different temperatures. We conclude that LC_x values are better suited to compare the difference in toxicity between substances or between the same substance at different test temperatures, while NEC values are better suited to estimating concentrations of substances that have no effect to the test species and endpoint measured under laboratory conditions.

Assessing the Relative Toxicity of Different Road Salts and Effect of Temperature on Salinity Toxicity: LCx Values versus No-Effect Concentration (NEC) Values

Freshwater biota are at risk globally from increasing salinity, including increases from deicing salts in cold regions. A variety of metrics of toxicity are used when estimating the toxicity of substances and comparing the toxicity between substances. However, the implications of using different metrics are not widely appreciated. Using the mayfly Colobruscoides giganteus (Ephemeroptera: Colobruscoidea), we compare the toxicity of seven different salts where toxicity was estimated using two metrics: (1) the no-effect concentrations (NEC) and (2) the lethal concentrations for 10, 25 and 50% of the test populations (LC_x). The LC_x values were estimated using two different models, the classic log-logistic model and the newer toxicokinetic-toxicodynamic (TKTD) model. The NEC and both types of LC_x values were estimated using Bayesian statistics. We also compared the toxicity of two salts (NaCl and CaCl₂) for C. giganteus at water temperatures of 4 °C, 7 °C and 15 °C using the same metrics of toxicity. Our motivation for using a mayfly to assess salinity toxicity was because mayflies are generally salt sensitive, are ecologically important and are common in Australian (sub-)alpine streams. The temperature ranges were chosen to mimic winter, spring and summer water temperatures for Australian (sub-)alpine streams. Considering 144-h classical LC_x values, we found toxicity differed between various salts, i.e., the lowest 144-h LC₅₀ (8 mS/cm) for a salt used by a ski resort was half that of the highest 144-h LC₅₀ from artificial marine salts and CaCl₂ applied to roads (16 mS/cm). The analytical grade NaCl (as shown by 144-h LC₅₀ value at 7 °C) was substantially more toxic (7.3 mS/cm) compared to analytical grade CaCl₂ (12.5 mS/cm). Yet for NEC values, there were comparably fewer differences in toxicity between salts and none between the same salts at different temperatures. We conclude that LC_x values are better suited to compare the difference in toxicity between substances or between the same substance at different test temperatures, while NEC values are better suited to estimating concentrations of substances that have no effect to the test species and endpoint measured under laboratory conditions.

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.

Issues, challenges, directives, and limitations concerning the improvement of environmental risk assessment of pharmaceutically active compounds

Nowadays, when tons of different chemicals including pharmaceutically active compounds (PhACs) are known to be released into the environment, applying adequate risk assessment relating to the protection of the ecosystem is very important. A broad body of scientific papers expresses a growing demand for improvement of the method(s) for ecological/environmental risk assessment (ERA). Although certain issues about ERA often emerge in the community, most of them cannot be considered as real problems and its methodology was developed keeping several limitations in mind. Nevertheless, the current approaches can be improved in order to better serve the intended purposes. For example, there is a lack of an integrated, manageable ecotoxicological database. It is not uncommon for basic, but extremely important, influencing factors such as time of exposure, interactions between different compounds, and characteristics of different habitats to be ignored. Discussing under the basic regulatory framework used in the EU, this correspondence paper deals with these and other examples to present the current features of ERA, identify gaps in process and application, and propose possible improvements/directives.

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

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

Improving Risk Assessment by Predicting the Survival of Field Gammarids Exposed to Dynamic Pesticide Mixtures

Exposure assessment of pesticides has substantially improved over time, with methods that now include a combination of advanced analytical techniques and fate/transport models to evaluate their spatiotemporal distribution. However, the current regulatory environmental risk assessment considers thresholds from laboratory studies completed under standardized conditions that do not reflect environmental dynamics. Using the General Unified Threshold model for Survival (GUTS) model framework, we predicted the impact of time-varying pesticide exposures on the survival of gammarids in a small agricultural stream. The LP₅₀ values were used as an additional metric for assessing risks (defined in GUTS as a multiplication factor applied to the concentration time series to induce 50% mortality by the end of exposure). Although real-case exposures to individual pesticides were predicted to produce little to no impact on survival, the LP₅₀ values indicate acute (LP₅₀ ≤ 100) and/or chronic (LP₅₀ ≤ 10) toxicities for azoxystrobin, chlorpyrifos, diazinon, and imidacloprid, while risk to propiconazole exposure was considered very low (LP₅₀ ≫ 100). Finally, the model was extended to reflect mixture toxicity via concentration addition. It predicted risks under acute and chronic exposures to organophosphates and neonicotinoids. Given that gammarids are simultaneously exposed to multiple chemicals and other stressors throughout their lifetime, a decline in survival probabilities due to chemical stress can likely influence their overall fitness. We recognize that some assumptions require validation, but our work included a level of realism that can assist risk managers when evaluating the cumulative consequences of chemical exposure.

The unlimited potential of the great pond snail, Lymnaea stagnalis

Only a limited number of animal species lend themselves to becoming model organisms in multiple biological disciplines: one of these is the great pond snail, Lymnaea stagnalis. Extensively used since the 1970s to study fundamental mechanisms in neurobiology, the value of this freshwater snail has been also recognised in fields as diverse as host-parasite interactions, ecotoxicology, evolution, genome editing and 'omics', and human disease modelling. While there is knowledge about the natural history of this species, what is currently lacking is an integration of findings from the laboratory and the field. With this in mind, this article aims to summarise the applicability of L. stagnalis and points out that this multipurpose model organism is an excellent, contemporary choice for addressing a large range of different biological questions, problems and phenomena.

Comparison of the General Threshold Model of Survival and Dose-Response Models in Simulating the Acute Toxicity of Metals to Danio rerio

We exposed zebrafish (Danio rerio) to different concentrations of lead and cadmium, and monitored them for survival at 24, 48, 72, and 96 h. Metal toxicity was predicted and compared using the dose-response and general threshold survival models in terms of required data sets, fit performance, and applicability. Environ Toxicol Chem 2019;38:2169-2177. © 2019 SETAC.

Applications of dynamic models in predicting the bioaccumulation, transport and toxicity of trace metals in aquatic organisms

This review evaluates the three dynamic models (biokinetic model: BK, physiologically based pharmacokinetic model: PBPK, and toxicokinetic-toxicodynamic model: TKTD) in our understanding of the key questions in metal ecotoxicology in aquatic systems, i.e., bioaccumulation, transport and toxicity. All the models rely on the first-order kinetics principle of metal uptake and elimination. The BK model basically treats organisms as a single compartment, and is both physiologically and geochemically based. With a good understanding of each kinetic parameter, bioaccumulation of metals in any aquatic organisms can be studied holistically and mechanistically. Modeling efforts are not merely restrained from the prediction of metal accumulation in the tissues, but instead provide the direction of the key processes that need to be addressed. PBPK is more physiologically based since it mainly addresses the transportation, transformation and distribution of metals in the organisms. It can be treated conceptually as a multi-compartmental kinetic model, whereas the physiology is driving the development of any good PBPK model which is no generic for aquatic animals and contaminants. There are now increasingly applications of the PBPK modeling specifically in metal studies, which reveal many important processes that are impossible to be teased out by direct experimental measurements without adequate modeling. TKTD models further focus on metal toxicity in addition to metal bioaccumulation. The TK part links exposure and bioaccumulation, while the TD part links bioaccumulation and toxic effects. The separation of TK and TD makes it possible to model processes, e.g., toxicity modification by environmental factors, interaction between different metals, at both the toxicokinetic and toxicodynamic levels. TKTD models provide a framework for making full use of metal toxicity data, and thus provide more information for environmental risk assessments. Overall, the three models reviewed here will continue to provide guiding principles in our further studies of metal bioaccumulation and toxicity in aquatic organisms.

Recommendations to address uncertainties in environmental risk assessment using toxicokinetic-toxicodynamic models

Providing reliable environmental quality standards (EQSs) is a challenging issue in environmental risk assessment (ERA). These EQSs are derived from toxicity endpoints estimated from dose-response models to identify and characterize the environmental hazard of chemical compounds released by human activities. These toxicity endpoints include the classical x% effect/lethal concentrations at a specific time t (EC/LC(x, t)) and the new multiplication factors applied to environmental exposure profiles leading to x% effect reduction at a specific time t (MF(x, t), or denoted LP(x, t) by the EFSA). However, classical dose-response models used to estimate toxicity endpoints have some weaknesses, such as their dependency on observation time points, which are likely to differ between species (e.g., experiment duration). Furthermore, real-world exposure profiles are rarely constant over time, which makes the use of classical dose-response models difficult and may prevent the derivation of MF(x, t). When dealing with survival or immobility toxicity test data, these issues can be overcome with the use of the general unified threshold model of survival (GUTS), a toxicokinetic-toxicodynamic (TKTD) model that provides an explicit framework to analyse both time- and concentration-dependent data sets as well as obtain a mechanistic derivation of EC/LC(x, t) and MF(x, t) regardless of x and at any time t of interest. In ERA, the assessment of a risk is inherently built upon probability distributions, such that the next critical step is to characterize the uncertainties of toxicity endpoints and, consequently, those of EQSs. With this perspective, we investigated the use of a Bayesian framework to obtain the uncertainties from the calibration process and to propagate them to model predictions, including LC(x, t) and MF(x, t) derivations. We also explored the mathematical properties of LC(x, t) and MF(x, t) as well as the impact of different experimental designs to provide some recommendations for a robust derivation of toxicity endpoints leading to reliable EQSs: avoid computing LC(x, t) and MF(x, t) for extreme x values (0 or 100%), where uncertainty is maximal; compute MF(x, t) after a long period of time to take depuration time into account and test survival under pulses with different periods of time between them.

Fit Reduced GUTS Models Online: From Theory to Practice

Mechanistic modeling approaches, such as the toxicokinetic-toxicodynamic (TKTD) framework, are promoted by international institutions such as the European Food Safety Authority and the Organization for Economic Cooperation and Development to assess the environmental risk of chemical products generated by human activities. TKTD models can encompass a large set of mechanisms describing the kinetics of compounds inside organisms (e.g., uptake and elimination) and their effect at the level of individuals (e.g., damage accrual, recovery, and death mechanism). Compared to classical dose-response models, TKTD approaches have many advantages, including accounting for temporal aspects of exposure and toxicity, considering data points all along the experiment and not only at the end, and making predictions for untested situations as realistic exposure scenarios. Among TKTD models, the general unified threshold model of survival (GUTS) is within the most recent and innovative framework but is still underused in practice, especially by risk assessors, because specialist programming and statistical skills are necessary to run it. Making GUTS models easier to use through a new module freely available from the web platform MOSAIC (standing for MOdeling and StAtistical tools for ecotoxIClogy) should promote GUTS operability in support of the daily work of environmental risk assessors. This paper presents the main features of MOSAIC_GUTS: uploading of the experimental data, GUTS fitting analysis, and LCx estimates with their uncertainty. These features will be exemplified from literature data. Integr Environ Assess Manag 2018;14:625-630. © 2018 SETAC.

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.