Extrapolating Metal (Cu, Ni, Zn) Toxicity from Individuals to Populations Across Daphnia Species Using Mechanistic Models: The Roles of Uncertainty Propagation and Combined Physiological Modes of Action.

Mechanistic effect modeling is a promising tool to improve the ecological realism of environmental risk assessment. An open question for the mechanistic modeling of metal toxicity is whether the same physiological mode of action (PMoA) could be assumed for closely related species. The implications of various modeling choices, such as the use of parameter point estimates and assumption of simplistic toxicodynamic models, are largely unexplored. We conducted life-table experiments with Daphnia longispina, Daphnia magna, and Daphnia pulex exposed to the single metals Cu, Ni, and Zn, and calibrated toxicokinetic-toxicodynamic (TKTD) models based on dynamic energy budget theory. We developed TKTD models with single and combined PMoAs to compare their goodness-of-fit and predicted population-level sensitivity. We identified the PMoA reproduction efficiency as most probable in all species for Ni and Zn, but not for Cu, and found that combined-PMoA models predicted higher population-level sensitivity than single-PMoA models, which was related to the predicted individual-level sensitivity, rather than to mechanistic differences between models. Using point estimates of parameters, instead of sampling from the probability distributions of parameters, could also lead to differences in the predicted population-level sensitivity. According to model predictions, apical chronic endpoints (cumulative reproduction, survival) are conservative for single-metal population effects across metals and species. We conclude that the assumption of an identical PMoA for different species of Daphnia could be justified for Ni and Zn, but not for Cu. Single-PMoA models are more appropriate than combined-PMoA models from a model selection perspective, but propagation of the associated uncertainty should be considered. More accurate predictions of effects at low concentrations may nevertheless motivate the use of combined-PMoA models. Environ Toxicol Chem 2024;43:338-358. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Hansul S ，Fettweis A ，Smolders E ，Schamphelaere K ... -  《-》

Interactive Metal Mixture Toxicity to Daphnia magna Populations as an Emergent Property in a Dynamic Energy Budget Individual-Based Model.

Environmental risk assessment of metal mixtures is challenging due to the large number of possible mixtures and interactions. Mixture toxicity data cannot realistically be generated for all relevant scenarios. Therefore, methods for prediction of mixture toxicity from single-metal toxicity data are needed. We tested how well toxicity of Cu-Ni-Zn mixtures to Daphnia magna populations can be predicted based on the Dynamic Energy Budget theory with an individual-based model (DEB-IBM), assuming non-interactivity of metals on the physiological level. We exposed D. magna populations to Cu, Ni, and Zn and their mixture at a fixed concentration ratio. We calibrated the DEB-IBM with single-metal data and generated blind predictions of mixture toxicity (population size over time), with account for uncertainty. We compared the predictive performance of the DEB-IBM with respect to mixture effects on population density and population growth rates with that of two reference models applied on the population level, independent action and concentration addition. Our inferred physiological modes of action (pMoA) differed from literature-reported pMoAs, raising the question of whether this is a result of different model selection approaches, intraspecific variability, or whether different pMoAs might actually drive toxicity in a population context. Observed mixture effects were concentration- and endpoint-dependent. The independent action was overall more accurate than the concentration addition but concentration addition-predicted effects on population growth rate were slightly better. The DEB-IBM most accurately predicted effects on 6-week density, including antagonistic effects at high concentrations, which emerged from non-interactivity at the physiological level. Mixture effects on initial population growth rate appear to be more difficult to predict. To explain why model accuracy is endpoint-dependent, relationships between individual-level and population-level endpoints should be illuminated. Environ Toxicol Chem 2021;40:3034-3048. © 2021 SETAC.

Hansul S ，Fettweis A ，Smolders E ，De Schamphelaere K ... -  《-》

The Unexpected Absence of Nickel Effects on a Daphnia Population at 3 Temperatures is Correctly Predicted by a Dynamic Energy Budget Individual-Based Model.

Recent studies have shown that temperature affects chronic nickel (Ni) toxicity to Daphnia magna at the individual (apical) level. However, the effect of temperature on Ni toxicity to D. magna at the population level is unknown. The present study investigated whether the effect of temperature on chronic Ni toxicity to D. magna assessed on apical endpoints can be extrapolated to the population level. The results of the population experiment showed no consistent Ni effects on total D. magna population abundance at 15, 20, and 25 °C, although the Ni concentrations tested were previously reported to significantly reduce reproduction in D. magna individuals. This result supports the idea that ecological risk assessment should not extrapolate as such from apical endpoints to the population level. A dynamic energy budget individual-based model (DEB-IBM) was calibrated using apical Ni toxicity data at 15, 20, and 25 °C. The goal was to investigate whether the calibrated DEB-IBM would be able to predict the unexpected absence of effects at the population level and to further investigate the effect of temperature on Ni toxicity to a D. magna population. At the population level, the calibrated DEB-IBM correctly predicted the unexpected absence of an effect of Ni on a D. magna population. Detailed analysis of simulation output suggests that the predicted lower Ni sensitivity at the population level occurs because Ni-induced mortality is compensated by reduced starvation (less intraspecific competition). Extrapolated median effective concentration (EC50) values for population density predicted that the effect of temperature on Ni toxicity to D. magna populations was smaller (1.9-fold higher at 25 °C than at 15 °C) than on Ni toxicity to D. magna apical reproduction (the EC50 is 6.5-fold higher at 25 °C than at 15 °C). These results show that the DEB-IBM can help to replace population experiments by in silico simulations and to optimize the experimental design of population studies. Environ Toxicol Chem 2019;38:1423-1433. © 2019 SETAC.

Pereira CMS ，Vlaeminck K ，Viaene K ，De Schamphelaere KAC ... -  《-》

Development and Validation of a Mixture Toxicity Implementation in the Dynamic Energy Budget-Individual-Based Model: Effects of Copper and Zinc on Daphnia magna Populations.

Mechanistic population models are gaining considerable interest in ecological risk assessment. The dynamic energy budget approach for toxicity (DEBtox) and the general unified threshold model for survival (GUTS) are well-established theoretical frameworks that describe sublethal and lethal effects of a chemical stressor, respectively. However, there have been limited applications of these models for mixtures of chemicals, especially to predict long-term effects on populations. We used DEBtox and GUTS in an individual-based model (IBM) framework to predict both single and combined effects of copper and zinc on Daphnia magna populations. The model was calibrated based on standard chronic toxicity test results with the single substances. A mixture toxicity implementation based on the general independent action model for mixtures was developed and validated with data from a population experiment with copper and zinc mixtures. Population-level effects of exposure to individual metals were accurately predicted by DEB-IBM. The DEB-IBM framework also allowed us to identify the potential mechanisms underlying these observations. Under independent action the DEB-IBM was able to predict the population dynamics observed in populations exposed to the single metals and their mixtures (R2  > 65% in all treatments). Our modeling shows that it is possible to extrapolate from single-substance effects at the individual level to mixture toxicity effects at the population level, without the need for mixture toxicity data at the individual level from standard mixture toxicity tests. The application of such modeling techniques can increase the ecological realism in risk assessment. Environ Toxicol Chem 2021;40:513-527. © 2020 SETAC.

Vlaeminck K ，Viaene KPJ ，Van Sprang P ，De Schamphelaere KAC ... -  《-》

The Use of Mechanistic Population Models in Metal Risk Assessment: Combined Effects of Copper and Food Source on Lymnaea stagnalis Populations.

Environmental risk assessment (ERA) of chemicals aims to protect populations, communities, and ecosystems. Population models are considered more frequent in ERA because they can bridge the gap between the individual and the population level. Lymnaea stagnalis (the great pond snail) is an organism that is particularly sensitive to various metals, including copper (Cu). In addition, the sensitivity of this species to Cu differs between food sources. The first goal of the present study was to investigate whether we could explain the variability in sensitivity between food sources (lettuce and fish flakes) at the individual level with a dynamic energy budget (DEB) model. By adapting an existing DEB model and calibrating it with Cu toxicity data, thereby combining information from 3 studies and 2 endpoints (growth and reproduction), we put forward inhibition of energy assimilation as the most plausible physiological mode of action (PMoA) of Cu. Furthermore, the variation in Cu sensitivity between both food sources was considerably lower at the PMoA level than at the individual level. Higher Cu sensitivity at individual level under conditions of lower food quality or availability appears to emerge from first DEB principles when inhibition of assimilation is the PMoA. This supports the idea that DEB explained Cu sensitivity variation between food sources. Our second goal was to investigate whether this food source effect propagated to the population level. By incorporating DEB in an individual-based model (IBM), population-level effects were predicted. Based on our simulations, the food source effect was still present at the population level, albeit less prominently. Finally, we compared predicted population-level effect concentration, x% (ECx) values with individual-level ECx values for different studies. Using the DEB-IBM, the range of effect concentrations decreased significantly: at the individual level, the difference in chronic EC10 values between studies was a factor of 70 (1.13-78 µg dissolved Cu/L), whereas at the population level the difference was a factor of 15 (2.9-44.6 µg dissolved Cu/L). To improve interstudy comparability, a bioavailability correction for differences in water chemistry was performed with a biotic ligand model. This further decreased the variation, down to a factor of 7.4. Applying the population model in combination with a bioavailability correction thus significantly decreased the variability of chronic effect concentrations of Cu for L. stagnalis. Overall, the results of the present study illustrate the potential usefulness of transitioning to a more modeling-based environmental risk assessment. Environ Toxicol Chem 2019;00:1-16. © 2019 SETAC.

Vlaeminck K ，Viaene KPJ ，Van Sprang P ，Baken S ，De Schamphelaere KAC ... -  《-》