Model-based experimental design for assessing effects of mixtures of chemicals

Predicting the Metal Mixture Toxicity with a Toxicokinetic-Toxicodynamic Model Considering the Time-Dependent Adverse Outcome Pathways

Chemicals mainly exist in ecosystems as mixtures, and understanding and predicting their effects are major challenges in ecotoxicology. While the adverse outcome pathway (AOP) and toxicokinetic-toxicodynamic (TK-TD) models show promise as mechanistic approaches in chemical risk assessment, there is still a lack of methodology to incorporate the AOP into a TK-TD model. Here, we describe a novel approach that integrates the AOP and TK-TD models to predict mixture toxicity using metal mixtures (specifically Cd-Cu) as a case study. We preliminarily constructed an AOP of the metal mixture through temporal transcriptome analysis together with confirmatory bioassays. The AOP revealed that prolonged exposure time activated more key events and adverse outcomes, indicating different modes of action over time. We selected a potential key event as a proxy for damage and used it as a measurable parameter to replace the theoretical parameter (scaled damage) in the TK-TD model. This refined model, which connects molecular responses to organism outcomes, effectively predicts Cd-Cu mixture toxicity over time and can be extended to other metal mixtures and even multicomponent mixtures. Overall, our results contribute to a better understanding of metal mixture toxicity and provide insights for integrating the AOP and TK-TD models to improve risk assessment for chemical mixtures.

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.

Evaluation of uncertainties in environmental impact assessment of naturally occurring radiation exposure situations on example of undisturbed and legacy NORM sites in the Fen Complex, Norway

The risk from naturally occurring radioactive materials (NORM) has been extensively assessed, and this has led to the integration of specific NORM radiation protection requirements within the latest EU Directive 2013/59. Nevertheless, it has been internationally recognised that remaining NORM knowledge gaps and uncertainties now present similarly significant issues in addressing recent regulatory requirements. The multi-tiered nature of environmental impact assessment (EIA) implies per se possibility for uncertainties, but when EIA at radiation exposure sites includes consideration of sites with multiple radiation and contamination sources, different ecosystem transport pathways, effects and risks by applying different parameters and models, level of overall uncertainty increases. The results of EIA study in the Fen area in Norway, comprised of undisturbed and legacy NORM sites, have been evaluated in this analysis, in order to identify all existing input uncertainties and how they may impact the final conclusions, and thus, influence any subsequent decision-making. The main uncertainties have been identified in the measurement and exposure analysis tier, and were related to the heterogeneous distribution of radionuclides, radionuclide speciation, as well as to generic variability issues in the concepts used for mobility and biota uptake analysis (such as distribution coefficient, transfer factors and concentration ratios) as well as radioecological modelling. The uncertainties in the input values to the calculation of the dose arising from radon exposure in the Fen area led to an overall elevated uncertainty of the magnitude of the radiation exposure dose of humans. It has been concluded that an integrated, ecosystem-based approach with consideration of complexity of prevailing environmental conditions and interconnections must be applied to fully understand possible radiation effects and risks.

Effect of fluoranthene on antioxidative defense in different tissues of Lymantria dispar and Euproctis chrysorrhoea larvae

This study examined the effect of long-term exposure to environmentally relevant concentrations of dietary fluoranthene (6.7 and 67 ng / g dry food weight) on defense mechanisms of the polyphagous forest insects Lymantria dispar L. and Euproctis chrysorrhoea L. The activities and expression of isoforms of superoxide dismutase (SOD) and catalase (CAT), the activities of glutathione S-transferase (GST) and glutathione reductase (GR), and total glutathione content (GSH) were determined in the whole midgut and midgut tissue, while SOD and CAT activities were assessed in hemolymph of the larvae. The results showed significant changes of enzyme activities, with more pronounced responses in larval midgut tissues, and between-species differences in patterns of response. Significantly increased activity of SOD was recorded in the whole midgut and midgut tissue of L. dispar larvae, as well as in midgut tissue of E. chrysorrhoea larvae. Fluoranthene increased CAT activity in midgut tissue of L. dispar larvae, and in the whole midgut and midgut tissue of E. chrysorrhoea larvae. Different expression patterns were detected for enzyme isoforms in tissues of larvae exposed to dietary fluoranthene. Total GSH content and GST activity increased in E. chrysorrhoea larval midgut tissue. Significantly decreased SOD activity in hemolymph of L. dispar larvae, and opposite changes in CAT activity were recorded in the hemolymph of larvae of two insect species. The tissue-specific responses of enzymes to dietary fluoranthene, recorded in each species, enabled the larvae to overcome the pollutant induced oxidative stress, and suggest further assessment of their possible use as early-warning signals of environmental pollution.

Joint effects of gamma radiation and cadmium on subcellular-, individual- and population-level endpoints of the green microalga Raphidocelis subcapitata

Interpreting and predicting the combined effects of toxicants in the environment is an important challenge in ecotoxicology. How such effects are connected across different levels of biological organisation is an additional matter of uncertainty. Such knowledge gaps are particularly prominent with regards to how ionising radiation interacts with contaminants. We assessed the response of twelve endpoints at the subcellular, individual and population level in a green microalga when exposed singly and jointly to gamma radiation and cadmium (Cd). We used a fully factorial experimental design where observed effects were compared to those predicted by the Independent Action (IA) model for mixture toxicity to determine whether they deviated from additivity. Subcellular endpoints (e.g., catalase, thiamine diphosphate, xanthophyll cycle pigments) showed an increased antioxidant and/or photoprotective response. However, our results indicate that this protection was not sufficient to prevent lipid peroxidation, which also increased with dose. At ecologically relevant doses, most interactions between gamma radiation and Cd regarding subcellular-, individual- and population-level endpoints were additive as predicted by the IA model. However, exposure to binary mixtures displayed antagonistic interactions between gamma radiation and Cd at the higher end of the tested dose spectrum. No correlations were observed between subcellular endpoints and higher-level endpoints, but there were linkages between individual and population endpoints. Our results suggest that antagonistic interactions between gamma radiation and Cd can occur at higher doses and that these interactions seem to disseminate from subcellular and individual to population level. Possible consequences for aquatic primary production and food-web interactions are discussed.

DNA damage exerted by mixtures of commercial formulations of glyphosate and imazethapyr herbicides in Rhinella arenarum (Anura, Bufonidae) tadpoles

Glyphosate (GLY) and imazethapyr (IMZT) are two herbicides commonly used worldwide, either alone or in mixtures. They represent key pesticides in modern agricultural management. The toxicity that results when employed as mixtures has not been characterized so far. Acute toxicity of the 48% GLY-based herbicide (GBH) Credit® and the 10.59% IMZT-based herbicide (IBH) Pivot® H alone and their binary combinations was analyzed in Rhinella arenarum tadpoles exposed in a semi-static renewal test. Lethal effects were determined using mortality as the end-point, whereas sublethal effects were determined employing the single-cell gel electrophoresis (SCGE) bioassay. Based on mortality experiments, results revealed LC5096 h values of 78.18 mg/L GBH and 0.99 mg/L IBH for Credit® and Pivot® H, respectively. An increase in the genetic damage index (GDI) was found after exposure to Credit® or Pivot® H at 5 and 10% of LC5096 h values. The combinations of 5% Credit®-5% Pivot® H LC5096 h and 10% Credit®-10% Pivot® H LC5096 h concentrations significantly enhanced the GDI in comparison with tadpoles exposed only to Credit® or Pivot® H. Thus, the effect of interaction between GBH and IBH inducing DNA damage in R. arenarum blood cells can be considered to be synergistic.

Environmental occurrence and hazard of organic UV stabilizers and UV filters in the sediment of European North and Baltic Seas

UV absorbing compounds are of emerging concern due to their large production volumes, their persistence or pseudo-persistence, and their potential adverse effects. This is the first study investigating the environmental occurrence and potential hazard of organic UV stabilizers and UV filters in the North and Baltic Sea surface sediments, including the connecting Skagerrak and Kattegat straits. In total, nineteen substances were identified over the entire study area, including the rarely studied compounds ethylhexyl triazone (EHT) and bisoctrizole (UV-360). Octocrylene (OC) was the predominant compound in this study with regard to detection frequency (79%) and concentrations (up to 9.7 ng/g dw). OC accounted for more than 65% of UV stabilizer contamination in the German Bight. The triazine derivative EHT was quantified in the Rhine-Meuse-Delta and the German Bight in concentrations up to 2.0 ng/g dw. In the Baltic Sea, benzotriazole UV stabilizers accounted for 60% of the contamination, with UV-360 as the main substance. The estimated environmental hazard quotients indicated a negligible impact on benthic and sediment-dwelling organisms in the North and Baltic Seas. Region-specific contamination pattern and riverine influences were revealed. The results suggest that both direct and indirect sources contribute to the UV stabilizer and UV filter contamination in the study area.

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

In ecological risk assessment of chemicals, hazard identification and hazard characterisation are most often based on ecotoxicological tests and expressed as summary statistics such as No Observed Effect Concentrations or Lethal Concentration values and No Effect Concentrations. Considerable research is currently ongoing to further improve methodologies to take into account toxico kinetic aspects in toxicological assessments, extrapolations of toxic effects observed on individuals to population effects and combined effects of multiple chemicals effects. In this context, the principles of the Dynamic Energy Budget (DEB), namely the conserved allocation of energy to different life-supporting processes in a wide variety of different species, have been applied successfully to the development of a number of DEB models. DEB models allow the incorporation of effects on growth, reproduction and survival within one consistent framework. This review aims to discuss the principles of the DEB theory together with available DEB models, databases available and applications in ecological risk assessment of chemicals for a wide range of species and taxa. Future perspectives are also discussed with particular emphasis on ongoing research efforts to develop DEB models as open source tools to further support the research and regulatory community to integrate quantitative biology in ecotoxicological risk assessment.

Interaction patterns and toxicities of binary and ternary pesticide mixtures to Daphnia magna estimated by an accelerated failure time model

Organisms in natural environments are often exposed to a broad variety of chemicals, and the multi-chemical mixtures exposure may produce significant toxic effects, even though the individual chemicals are present at concentrations below their no-observed-effect concentrations. This study represents the first attempt that uses the accelerated failure time (AFT) model to quantify the interaction and toxicity of multi-chemical mixtures in environmental toxicology. We firstly conducted the acute immobilization tests with Daphnia magna exposed to mixtures of diazinon (DZN), fenitrothion (MEP); and thiobencarb (TB) in single, binary, and ternary formulations, and then fitted the results to the AFT model. The 48-h EC₅₀ (concentration required to immobilize 50% of the daphnids at 48h) values for each pesticide obtained from the AFT model are within a factor of 2 of the corresponding values calculated from the single pesticide exposure tests, indicating the methodology is able to provide credible toxicity values. The AFT model revealed either significant synergistic (DZN and MEP; DZN and TB) or antagonistic (MEP and TB) interactions in binary mixtures, while the interaction pattern of ternary mixture depended on both the concentration levels and concentration ratios of pesticides. With a factor of 2, the AFT model accurately estimated the toxicities for 78% of binary mixture formulations that exhibited significant synergistic effects, and the toxicities for all the ternary formulations. Our results showed that the AFT model can provide a simple and efficient way to quantify the interactions between pesticides and to assess the toxicity of their mixtures. This ability may greatly facilitate the ecotoxicological risk assessment of exposure to multi-chemical mixtures.

Comparative toxicity of pesticides and environmental contaminants in bees: Are honey bees a useful proxy for wild bee species?

Threats to wild and managed insect pollinators in Europe are cause for both ecological and socio-economic concern. Multiple anthropogenic pressures may be exacerbating pollinator declines. One key pressure is exposure to chemicals including pesticides and other contaminants. Historically the honey bee (Apis mellifera spp.) has been used as an 'indicator' species for 'standard' ecotoxicological testing but it has been suggested that it is not always a good proxy for other types of eusocial and solitary bees because of species differences in autecology and sensitivity to various stressors. We developed a common toxicity test system to conduct acute and chronic exposures of up to 240h of similar doses of seven chemicals, targeting different metabolic pathways, on three bee species (Apis mellifera spp., Bombus terrestris and Osmia bicornis). We compared the relative sensitivity between species in terms of potency between the chemicals and the influence of exposure time on toxicity. While there were significant interspecific differences that varied through time, overall the magnitude of these differences (in terms of treatment effect ratios) was generally comparable (<2 fold) although there were some large divergences from this pattern. Our results suggest that A. mellifera spp. could be used as a proxy for other bee species provided a reasonable assessment factor is used to cover interspecific variation. Perhaps more importantly our results show significant and large time dependency of toxicity across all three tested species that greatly exceeds species differences (>25 fold within test). These are rarely considered in standard regulatory testing but may have severe environmental consequences, especially when coupled with the likelihood of differential species exposures in the wild. These insights indicate that further work is required to understand how differences in toxicokinetics vary between species and mixtures of chemicals.

The Effects of Heavy Metals and Total Petroleum Hydrocarbons on Soil Bacterial Activity and Functional Diversity in the Upper Silesia Industrial Region (Poland)

Various inorganic and organic pollutants in industrial soils may adversely affect soil microorganisms and terrestrial ecosystem functioning. The aim of the study was to explore the relationship between the microbial activity, microbial biomass, and functional diversity of soil bacteria and the metals and total petroleum hydrocarbons (TPHs) in the Upper Silesian Industrial Region (Poland). We collected soil samples in pine-dominated forest stands and analyzed them according to a range of soil physicochemical properties, including metal content (cadmium, lead, and zinc) and TPH content. Metal concentrations were normalized to their toxicity to soil microorganisms and integrated in a toxicity index (TI). Soil microbial activity measured as soil respiration rate, microbial biomass measured as substrate-induced respiration rate, and the bacterial catabolic activity (area under the curve, AUC) assessed using Biolog® ECO plates were negatively related to TPH pollution as shown in multiple regressions. The canonical correspondence analysis (CCA) showed that both TPH and TI affected the community-level physiological profiles (CLPPs) of soil bacteria and the pollutants' effects were much stronger than the effects of other soil properties, including nutrient content.

Joint Toxicity of Cadmium and Ionizing Radiation on Zooplankton Carbon Incorporation, Growth and Mobility

The risk of exposure to radioactive elements is seldom assessed considering mixture toxicity, potentially over- or underestimating biological and ecological effects on ecosystems. This study investigated how three end points, carbon transfer between phytoplankton and Daphnia magna, D. magna mobility and growth, responded to exposure to γ-radiation in combination with the heavy metal cadmium (Cd), using the MIXTOX approach. Observed effects were compared with mixture effects predicted by concentration addition (CA) and independent action (IA) models and with deviations for synergistic/antagonistic (S/A), dose-level (DL), and dose-ratio (DR) dependency interactions. Several patterns of response were observed depending on the end point tested. DL-dependent deviation from the IA model was observed for carbon incorporation with antagonism switching to synergism at higher doses, while the CA model indicated synergism, mainly driven by effects at high doses of γ-radiation. CA detected antagonism regarding acute immobilization, while IA predicted DR-dependency. Both CA and IA also identified antagonism for daphnid growth. In general, effects of combinations of γ-radiation and Cd seem to be antagonistic at lower doses, but synergistic at the higher range of the doses tested. Our results highlight the importance of investigating the effects of exposure to γ-radiation in a multistressor context.

Effects of fluoranthene on the fitness-related traits and antioxidative defense in Lymantria dispar L

This study aimed to examine the effects of ubiquitous polycyclic aromatic hydrocarbon fluoranthene, supplemented to an artificial diet, on the fitness-related traits and activity of midgut antioxidative enzymes-superoxide dismutase (SOD) and catalase (CAT), and expression of their isoforms in the fifth-instar gypsy moth Lymantria dispar L. Prolonged duration of development and reduced weight and relative growth rate were recorded in larvae reared on the diets supplemented with different concentrations of fluoranthene. SOD and CAT activities were significantly higher in the midguts of fluoranthene-treated larvae, compared to that of the control group. Different expression patterns were detected for SOD as well as for CAT isoforms, depending on the supplemented concentration of fluoranthene. Obtained results suggest that the activity of these enzymes in gypsy moth larvae may be used as biomarkers for assessing pollution, even at low concentrations of the pollutant.

Tools and perspectives for assessing chemical mixtures and multiple stressors

The present paper summarizes the most important insights and findings of the EU NoMiracle project with a focus on (1) risk assessment of chemical mixtures, (2) combinations of chemical and natural stressors, and (3) the receptor-oriented approach in cumulative risk assessment. The project aimed at integration of methods for human and ecological risk assessment. A mechanistically based model, considering uptake and toxicity as a processes in time, has demonstrated considerable potential for predicting mixture effects in ecotoxicology, but requires the measurement of toxicity endpoints at different moments in time. Within a novel framework for risk assessment of chemical mixtures, the importance of environmental factors on toxicokinetic processes is highlighted. A new paradigm for applying personal characteristics that determine individual exposure and sensitivity in human risk assessment is suggested. The results are discussed in the light of recent developments in risk assessment of mixtures and multiple stressors.

A review of multiple stressor studies that include ionising radiation

Studies were reviewed that investigated the combined effects of ionising radiation and other stressors on non-human biota. The aim was to determine the state of research in this area of science, and determine if a review of the literature might permit a gross generalization as to whether the combined effects of multi-stressors and radiation are fundamentally additive, synergistic or antagonistic. A multiple stressor database was established for different organism groups. Information was collected on species, stressors applied and effects evaluated. Studies were mostly laboratory based and investigated two-component mixtures. Interactions declared positive occurred in 58% of the studies, while 26% found negative interactions. Interactions were dependent on dose/concentration, on organism's life stage and exposure time and differed among endpoints. Except for one study, none of the studies predicted combined effects following Concentration Addition or Independent Action, and hence, no justified conclusions can be made about synergism or antagonism.

Ecotoxicological mechanisms and models in an impact analysis tool for oil spills

In an international collaborative effort, an impact analysis tool is being developed to predict the effect of accidental oil spills on recruitment and production of Atlantic cod (Gadus morhua) in the Barents Sea. The tool consisted of three coupled ecological models that describe (1) plankton biomass dynamics, (2) cod larvae growth, and (3) fish stock dynamics. The discussions from a series of workshops are presented in which variables and parameters of the first two ecological models were listed that may be affected by oil-related compounds. In addition, ecotoxicological algorithms are suggested that may be used to quantify such effects and what the challenges and opportunities are for algorithm parameterization. Based on model exercises described in the literature, survival and individual growth of cod larvae, survival and reproduction of zooplankton, and phytoplankton population growth are denoted as variables and parameters from the ecological models that might be affected in case of an oil spill. Because toxicity databases mostly (67%) contain data for freshwater species in temperate environments, parameterization of the ecotoxicological algorithms describing effects on these endpoints in the subarctic marine environment is not straightforward. Therefore, it is proposed that metadata analyses be used to estimate the sensitivity of subarctic marine species from available databases. To perform such analyses and reduce associated uncertainty and variability, mechanistic models of varying complexity, possibly aided by new experimental data, are proposed. Lastly, examples are given of how seasonality in ecosystems may influence chemical effects, in particular in the subarctic environment. Food availability and length of day were identified as important characteristics as these determine nutritional status and phototoxicity, respectively.