Linking Morphology, Toxicokinetic, and Toxicodynamic Traits of Aquatic Invertebrates to Pyrethroid Sensitivity

Improving Sediment Toxicity Testing for Very Hydrophobic Chemicals: Part 2-Exposure Duration, Upper Limit Test Concentrations, and Distinguishing Actual Toxicity from Physical Effects

Sediment toxicity testing with very hydrophobic organic chemicals (VHOCs) is challenging because of the chemicals' low aqueous solubilities and slow kinetics. The present study presents the results of experiments investigating whether the standard exposure duration of 28 days with benthic invertebrates is sufficient for VHOCs; above which concentrations in sediment VHOCs are present as "free phase," that is, crystals or non-aqueous-phase liquids (NAPLs); and whether it is possible to discriminate between actual VHOC toxicity and physical effects caused by NAPLs through fouling of the test organisms. The results suggest that the standard sediment toxicity test duration is sufficient for obtaining steady-state VHOC concentrations in Hyalella azteca and Lumbriculus variegatus, provided that spiking and equilibration are performed properly (i.e., no free phase present). Under these conditions, transient (days 3-20) peak-shaped toxicokinetics were observed, with steady-state concentrations reached at approximately 28 days. The concentration above which NAPLs are present, the so-called critical separate phase concentration (CSPC), was determined for several VHOCs by modeling and two experimental methods. Modeling resulted in unrealistic and variable data and therefore should be applied with caution. Experimentally determining CSPCs was successful and yielded values of approximately 1000 (400-2000) mg/kg dry weight, depending on the chemical. Finally, it was demonstrated that distinguishing actual toxicity from physical effects is possible by applying a well-considered test setup, combining toxicity tests with multiple invertebrates (including Lumbriculus, which serves as a negative control for fouling); a broad test concentration range, preferably up to at least 30 000 mg/kg; and passive sampling to localize the CSPC. Applying this setup, false-positive effects due to fouling, as well as false-negative results due to testing at too low concentrations (trying to stay below the CSPC), can be avoided. Environ Toxicol Chem 2024;43:1728-1739. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

The Importance of Including Variable Exposure Concentrations When Assessing Toxicity of Sediment-Associated Pharmaceuticals to an Amphipod

Pharmaceuticals have been classified as an environmental concern due to their increasing consumption globally and potential environmental impact. We examined the toxicity of sediment-associated diclofenac and citalopram administered as both single compounds and in a mixture to the sediment-living amphipod Corophium volutator. This laboratory-based study addressed the following research questions: (1) What is the toxicity of sediment-associated diclofenac and citalopram to C. volutator? (2) Can the mixture effect be described with either of the two mixture models: concentration addition (CA) or independent action (IA)? (3) What is the importance of the choice of (i) exposure measure (start concentration, time-weighted average [TWA], full exposure profile) and (ii) effect model (concentration-response vs. the toxicokinetic-toxicodynamic model general unified threshold model for survival in its reduced form [GUTS-RED]) for the derived effect concentration values? Diclofenac was more toxic than citalopram to C. volutator as a single compound (10-day exposure). Diclofenac exposure to C. volutator provided median lethal concentrations (LC50s) within the same range (11 µg g-1 dry wt sediment) using concentration-response based on TWA and both GUTS-RED models. However, concentration-response based on measured start concentrations provided an approximately 90% higher LC50 (21.6 ± 2.0 µg g-1 dry wt sediment). For citalopram, concentration-response parameters were similar regardless of model or concentration used (LC50 85-97 µg g-1 dry wt sediment), however, GUTS-RED with the assumption of individual tolerance resulted in a lower LC50 (64.9 [55.3-74.8] µg g-1 dry wt sediment). The mixture of diclofenac and citalopram followed the CA quite closely, whereas the result was synergistic when using the IA prediction. In summary, concentration-response based on TWA and GUTS-RED provided similar and reasonably good fits compared to the data set. The implications are that GUTS-RED will provide a more flexible model, which, in principle, can extend beyond the experimental period and make predictions based on variable exposure profiles (toxicity at different time frames and at different variable exposure scenarios) compared to concentration-response, which provides contaminant toxicity at one point in time. Environ Toxicol Chem 2024;43:1767-1777. © 2024 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Effects of α-olefin sulfonate (AOS) on Tubifex tubifex: toxicodynamic-toxicokinetic inferences from the general unified threshold (GUTS) model, biomarker responses and molecular docking predictions

We investigated the potential ecological risks and harm to aquatic organisms posed by anionic surfactants such as α-olefin sulfonate (AOS), which are commonly found in industrial and consumer products, including detergents. This study assessed acute (96-h) and subchronic (14-day) responses using antioxidant activity, protein levels, and histopathological changes in Tubifex tubifex exposed to different AOS concentrations (10% of the LC50, 20% of the LC50, and a control). Molecular docking was used to investigate the potential interactions between the key stress biomarker enzymes (superoxide dismutase, catalase, and cytochrome c oxidase) of Tubifex tubifex. Acute AOS exposure showed a concentration-dependent decrease in survival, and the general unified threshold (GUTS) model revealed that survivorship is linked to individual response patterns rather than random (stochastic) fluctuations. The GUTS model also revealed dose-dependent toxicity patterns in Tubifex tubifex exposed to α-olefin sulfonate (AOS), with adaptive mechanisms at lower concentrations but significant increases in mortality beyond a certain threshold, emphasizing the role of the AOS concentration in shaping its toxicological impact. Exposure to AOS disrupted antioxidant activity, inducing oxidative stress, with GST and GPx showing positive associations with surfactant concentration and increased lipid peroxidation (elevated MDA levels); moreover, AOS exposure decreased protein concentration, signifying disturbances in vital cellular processes. Histopathological examinations revealed various tissue-level alterations, including cellular vacuolation, cytoplasmic swelling, inflammation, necrosis, and apoptosis. Molecular docking analysis demonstrated interactions between AOS and enzymes (-catalase, superoxide dismutase, and cytochrome c oxidase) in Tubifex tubifex, including hydrophobic and hydrogen bond interactions, with the potential to disrupt enzyme structures and activities, leading to cellular process disruptions, oxidative stress, and tissue damage. According to the species sensitivity distribution (SSD), the difference in toxicity between Tilapia melanopleura (higher sensitivity) and Daphnia magna (low sensitivity) to AOS suggests distinct toxicokinetic and toxicodynamic mechanisms attributable to more complex physiology in Tilapia and efficient detoxification in Daphnia due to its smaller size.

Spatial variation in the sensitivity of freshwater macroinvertebrate assemblages to chemical stressors

Assessing spatial variation in the chemical sensitivity of natural assemblages will enhance ecological relevance and reduce uncertainty in ecological risk assessments and the derivation of environmental quality standards (EQSs). However, the majority of species in natural communities have not undergone toxicity testing for any chemical, which poses a major challenge when assessing their sensitivity. We investigated spatial variation and patterns in the sensitivity of 4084 freshwater macroinvertebrate assemblages across England to 5 general-acting chemicals (heavy metals) and 13 specifically acting chemicals (insecticides) using a novel hierarchical species sensitivity distribution method based on taxonomic relatedness. Furthermore, we explored how river typology relates to spatial variation in assemblage sensitivity to chemicals and the potential impacts of such variation on current EQSs. Our findings revealed that, whereas assemblages with similar taxonomic compositions exhibit comparable sensitivity distributions, assemblages with different taxonomic compositions could have very similar or very different sensitivity distributions. The variation in assemblage sensitivity was greater for specifically acting chemicals than for general-acting chemicals and exhibited spatial clustering patterns. These spatial clustering patterns varied depending on the chemical, and the regions where assemblages were most sensitive to metals were generally not the same as the regions where assemblages were most sensitive to insecticides. Spatial variation in assemblage sensitivity was related to river typology with sensitive assemblages being more common than expected in lowland calcareous (or mixed geology) rivers within very small to small catchments. Comparing spatial variation in assemblage-specific chemical sensitivity to EQSs, we found that the operational EQSs in England would protect most study assemblages (i.e., > 99.5 %), although a small proportion of assemblages may face potential risks associated with azinphos-methyl, copper, and malathion. In many cases the EQSs were very precautionary, potentially requiring expensive control measures or restricting beneficial chemical use with no additional environmental benefit. The development of spatially defined EQSs, possibly based on river types, could be developed to target areas that require the highest level of protection and thus strike a balance between the benefits of chemical use and environmental protection.

Influence of soil organic matter content on the toxicity of pesticides to soil invertebrates: A review

A better understanding of how soil organic matter (OM) content influences pesticide toxicity to soil invertebrates is needed to improve the ecological relevance of risk assessment approaches. In the current study, soil invertebrate toxicity data (LC50 and EC50 values) were collected from studies determining the toxicity of organic chemicals in soils with varying OM content. Relevant studies were identified by performing a literature search and through the use of toxicity databases. The data were used to address the following questions: (1) Can the relationship between toxicity and soil OM content be quantified? (2) Does soil OM content influence different toxicity endpoints in a similar way? (3) Is the influence of soil OM content on sensitivity to pesticides different between species? The results indicate that toxicity-OM relationships are chemical dependent, differ between endpoints, and are species-specific. Hence, the grouping of chemicals based solely on their lipophilicity, as well as having only one correction factor for multiple species, may not be an appropriate approach to risk assessment. Integr Environ Assess Manag 2023;19:1457-1472. © 2023 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Short-term insecticide exposure amid co-occurring stressors reduces diversity and densities in north-east Indian experimental aquatic invertebrate communities

Globally, river pesticide concentrations are associated with regional and local stream invertebrate diversity declines. Pesticides often co-occur with elevated nutrients (e.g. nitrogen and phosphorus) and sediments related to agriculture, making their individual effects difficult to disentangle. These effects are also less well studied in Asia, than in other geographic regions. Within Asia, India is one of the largest producers and users of pesticides and has approximately 60% of total land mass used for agriculture. Here we examine the responses of Indian river invertebrate communities subjected to malathion, nutrients, and sediment additions in a semi-orthogonal design, in three sequential (through time) short-term (120 h) mesocosm experiments. Additionally, a series of single-species toxicity tests were run that used 24 h exposure and 72 h recovery to examine the sensitivity of 13 local invertebrate taxa to malathion, and 9 taxa to cypermethrin, comparing these results to those from other biogeographic regions. Mesocosm results indicate that malathion exposure had a major effect compared to other stressors on communities, with a lesser effect of nutrients and/or sediments. In mesocosms, taxa richness, total abundance and the abundance of sensitive species all declined associated with malathion concentrations. Comparisons of organism sensitivities from other geographic locations and those in the current paper suggest taxa in India are relatively tolerant to malathion and cypermethrin. Our results further reinforce that the high observed aquatic pesticide concentrations known to occur in Asian freshwater ecosystems are likely to be negatively affecting biodiversity, homogenising biota towards those most stress tolerant.

Two-Compartmental Toxicokinetic Model Predicts Interspecies Sensitivity Variation of Imidacloprid to Aquatic Invertebrates

Interspecies sensitivity to the same chemical can be several orders of magnitude different. Quantifying toxicologically internal levels and toxicokinetic (TK) parameters is critical in elucidating the interspecies sensitivity. Herein, a two-compartmental TK model was constructed to characterize the uptake, distribution, and elimination kinetics toward interspecies sensitivity to an insecticide, imidacloprid. Imidacloprid exhibited the highest lethality to the insect Chironomus dilutus, followed by Lumbriculus variegatus, Hyalella azteca, and Daphnia magna. Interspecies sensitivity of imidacloprid to these invertebrates varied by ∼1000 folds based on water concentrations (LC50). Remarkably, the sensitivity variation decreased to ∼50 folds based on the internal residues (LR50), highlighting the critical role of TK in interspecies sensitivity. A one-compartmental TK model failed to simulate the bioaccumulation of imidacloprid in these invertebrates except for D. magna. Instead, a two-compartmental model successfully simulated the slow elimination of imidacloprid in the remaining three species by internally distinguishing the highly dynamic (C1) and toxicologically available (C2) fractions. We further showed that the species sensitivity of the invertebrates to imidacloprid was significantly related to C2, demonstrating that C2 was toxicologically available and responsible for the toxicity of imidacloprid. This mechanistic-based model bridged the internal distribution of organic contaminants in small invertebrates and the associated toxic potency.

Exposure to cypermethrin pesticide disturbs the microbiome and disseminates antibiotic resistance genes in soil and the gut of Enchytraeus crypticus

Worldwide, pyrethroids, such as cypermethrin, are the second most applied group of insecticides, however, their effects on the soil microbiome and non-target soil fauna remain largely unknown. Herein, we assessed the change of bacterial communities and antibiotic resistance genes (ARGs) of soil and in the gut of the model soil species Enchytraeus crypticus using a combination of 16S rRNA gene amplicon sequencing, and high-throughput qPCR of ARGs. Results indicate that cypermethrin exposure enriches potential pathogens (e.g. Bacillus anthracis) in the soil and gut microbiome of E. crypticus, heavily disrupting the latter's microbiome structure, and even disrupts activities of the E. crypticus immune system. The co-occurrence of potential pathogens (e.g. Acinetobacter baumannii), ARGs, and mobile genetic elements (MGEs) revealed the increased risk of pathogenicity as well as antibiotic resistance in potential pathogens. Moreover, structural equation modeling demonstrated that the dissemination of ARGs was not only promoted by MGEs, but also by the ratio of the core to non-core bacterial abundance. Collectively, these results provide an in-depth view of the previously unappreciated environmental risk of cypermethrin on the dissemination of ARGs in the soil and non-target soil fauna.

Species sensitivity distributions for two widely used anti-sea lice chemotherapeutants in the salmon aquaculture industry

The main objective of the present study is to construct acute aquatic species sensitivity distributions (SSD) and generate proposed HC₅ values (i.e. the hazardous concentration for which 5 % of species are affected or potentially affected) for two aquaculture anti-sea lice bath pesticides, azamethiphos, and hydrogen peroxide. These values could be used as the basis for the establishment of environmental quality standards (EQS). We have generated SSDs and inferred HC₅ values for mortality and sublethal endpoints using LC50, EC50, and NOEC/LOEC data points separately and for each bath pesticide. Through the examination of literature data on the toxicity of both compounds, we opted to use tests with limited exposure times to ensure environmental relevance for bath pesticides. We also separated life stages for some of the sensitive taxa to account for differences in sensitivities and risk of exposure. The resulting threshold concentrations in environmental seawater are 0.10 μg/L for azamethiphos and 0.15 mg/L for hydrogen peroxide. These suggested azamethiphos and hydrogen peroxide thresholds are comparable to some previously reported EQS values. Further considerations need to be included in how to better use these thresholds in a regulatory context in relation to dispersion patterns. It is also clear that delayed mortality and sublethal effects documented in the literature require further study to fully anticipate the environmental risks of using these two bath pesticides.

Joint survival modelling for multiple species exposed to toxicants

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

Elucidating the Differences in Metal Toxicity by Quantitative Adverse Outcome Pathways

Numerous studies have reported that the toxicity differences among metals are widespread; however, little is known about the mechanism of differences in metal toxicity to aquatic organisms due to the lack of quantitative understanding of their adverse outcome pathway. Here, we investigated the effects of Cd and Cu on bioaccumulation, gene expression, physiological responses, and apical effects in zebrafish larvae. RNA sequencing was conducted to provide supplementary mechanistic information for the effects of Cd and Cu exposure. On this basis, we proposed a quantitative adverse outcome pathway (qAOP) suitable for metal risk assessment of aquatic organisms. Our work provides a mechanistic explanation for the differences in metal toxicity where the strong bioaccumulation of Cu enables the newly accumulated Cu to reach the threshold that causes different adverse effects faster than Cd in zebrafish larvae, resulting in a higher toxicity of Cu than that of Cd. Furthermore, we proposed a parameter C_IT/BCF (the ratio of internal threshold concentration and bioaccumulation factor) that helps to understand the toxicity differences by combining the information of bioaccumulation and internal threshold of adverse effects. This work demonstrated that qAOP is an effective quantitative tool for understanding the toxicity mechanism and highlight the importance of toxicokinetics and toxicodynamics at different biological levels in determining the metal toxicity.

Comparing the acute and chronic toxicity of flupyradifurone and imidacloprid to non-target aquatic arthropod species

Flupyradifurone (FPF) is a new type of butenolide insecticide. It was launched on the market in 2015 and is considered an alternative to the widely used neonicotinoids, like imidacloprid (IMI), some of which are banned from outdoor use in the European Union. FPF is claimed to be safe for bees, but its safety for aquatic organisms is unknown. Its high water solubility, persistence in the environment, and potential large-scale use make it urgent to evaluate possible impacts on aquatic systems. The current study assessed the acute and chronic toxicity of FPF for aquatic arthropod species and compared these results with those of imidacloprid. Besides, toxicokinetics and toxicokinetic-toxicodynamic models were used to understand the mechanisms of the toxicity of FPF. The present study results showed that organisms take up FPF slower than IMI and eliminate it faster. In addition, the hazardous concentration 5th percentiles (HC₀₅) value of FPF derived from a species sensitivity distribution (SSD) based on acute toxicity was found to be 0.052 µmol/L (corresponding to 15 µg/L), which was 37 times higher than IMI (0.0014 µmol/L, corresponding to 0.36 µg/L). The chronic 28 days EC₁₀ of FPF for Cloeon dipterum and Gammarus pulex were 7.5 µg/L and 2.9 µg/L, respectively. For G. pulex, after 28 days of exposure, the no observed effect concentration (NOEC) of FPF for food consumption was 0.3 µg/L. A toxicokinetic-toxicodynamic (TKTD) model parameterised on the acute toxicity data well predicted the observed chronic effects of FPF on G. pulex, indicating that toxicity mechanisms of FPF did not change with prolonged exposure time, which is not the case for IMI.

Exposures to deltamethrin on immature Chironomus columbiensis drive sublethal and transgenerational effects on their reproduction and wing morphology

Sublethal exposure to insecticides can trigger unintended responses in non-target insects that may disrupt reproductive and developmental performances of these organisms. Here, we assessed whether sublethal exposure to the pyrethroid insecticide deltamethrin in early life had sublethal and transgenerational effects on the reproduction (i.e., fecundity and fertility) and wing morphology of Chironomus columbiensis, an aquatic insect used as a water quality indicator. We first conducted concentration-response bioassays to evaluate the susceptibility of C. columbiensis larvae to deltamethrin. Our results revealed that deltamethrin toxicity was approximately 7-fold higher when C. columbiensis larvae where exposed to 96 h (LC₅₀ = 0.17 [0.15-0.20] μg/L) than to 24 h (LC₅₀ = 1.17 [0.97-1.43] μg/L). Furthermore, the sublethal exposures (at LC₁ = 0.02 μg/L or LC₁₀ = 0.05 μg/L) of immature C. columbiensis resulted in lower fecundity (e.g., reduced eggs production) and morphometric variation wing shapes. Further reduction in fertility rates (quantity of viable eggs) occurred at deltamethrin LC₁₀ (0.05 μg/L). Almost 80% of the fecundity was recovered with only a single recovery generation; however, two subsequent recovery generations were not sufficient to fully recover fecundity in C. columbiensis. Specimens recovered from 98.5% of wing morphometric variation after two consecutive generations without deltamethrin exposure. Collectively, our findings demonstrates that sublethal exposure to synthetic pyrethroids such as deltamethrin detrimentally affect the reproduction and wing shape of C. columbiensis, but also indicate that proper management of these compounds (e.g., concentration and frequency of application) would suffice for these insects' population recovery.

A critical review of effect modeling for ecological risk assessment of plant protection products

A wide diversity of plant protection products (PPP) is used for crop protection leading to the contamination of soil, water, and air, which can have ecotoxicological impacts on living organisms. It is inconceivable to study the effects of each compound on each species from each compartment, experimental studies being time consuming and cost prohibitive, and animal testing having to be avoided. Therefore, numerous models are developed to assess PPP ecotoxicological effects. Our objective was to provide an overview of the modeling approaches enabling the assessment of PPP effects (including biopesticides) on the biota. Six categories of models were inventoried: (Q)SAR, DR and TKTD, population, multi-species, landscape, and mixture models. They were developed for various species (terrestrial and aquatic vertebrates and invertebrates, primary producers, micro-organisms) belonging to diverse environmental compartments, to address different goals (e.g., species sensitivity or PPP bioaccumulation assessment, ecosystem services protection). Among them, mechanistic models are increasingly recognized by EFSA for PPP regulatory risk assessment but, to date, remain not considered in notified guidance documents. The strengths and limits of the reviewed models are discussed together with improvement avenues (multigenerational effects, multiple biotic and abiotic stressors). This review also underlines a lack of model testing by means of field data and of sensitivity and uncertainty analyses. Accurate and robust modeling of PPP effects and other stressors on living organisms, from their application in the field to their functional consequences on the ecosystems at different scales of time and space, would help going toward a more sustainable management of the environment. Graphical Abstract Combination of the keyword lists composing the first bibliographic query. Columns were joined together with the logical operator AND. All keyword lists are available in Supplementary Information at https://doi.org/10.5281/zenodo.5775038 (Larras et al. 2021).

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.

Using TKTD Models in Combination with In Vivo Enzyme Inhibition Assays to Investigate the Mechanisms behind Synergistic Interactions across Two Species

The aim of this study is to compare the azole synergy across an insect, Chironomus riparius, and a crustacean species, Daphnia magna. We use a combination of in vivo measurements of cytochrome P450 monooxygenase (CYP) biotransformation potential and toxicokinetic (TK) and toxicodynamic (TD) modeling to understand the mechanism behind the synergy of two azole fungicides: the imidazole prochloraz and the triazole propiconazole on the pyrethroid insecticide α-cypermethrin. For both species, the synergistic effect of prochloraz was well-described by its effect on in vivo CYP activity, which corresponded to the biotransformation rate of the TK model parameterized on the survival data of the mixture experiment. For propiconazole, however, there were 100-fold and 50-fold differences between the 50% effect concentration of in vivo CYP activity and the modeled biotransformation rate for C. riparius and D. magna, respectively. Propiconazole, therefore, seems to induce synergy through a mechanism that cannot be quantified solely by the CYP activity assay used in this study in either of the two species. We discuss the differences between prochloraz and propiconazole as synergists across the two species in the light of the type and time dynamics of affected biotransformation processes.

Integrating Bioavailability of Metals in Fish Population Models

Population models are increasingly being used to extrapolate individual-level effects of chemicals, including metals, to population-level effects. For metals, it is also important to take into account their bioavailability to correctly predict metal toxicity in natural waters. However, to our knowledge, no models exist that integrate metal bioavailability into population modeling. Therefore, our main aims were to 1) incorporate the bioavailability of copper (Cu) and zinc (Zn) into an individual-based model (IBM) of rainbow trout (Oncorhynchus mykiss), and 2) predict how survival-time concentration data translate to population-level effects. For each test water, reduced versions of the general unified threshold model of survival (GUTS-RED) were calibrated using the complete survival-time concentration data. The GUTS-RED individual tolerance (IT) showed the best fit in the different test waters. Little variation between the different test waters was found for 2 GUTS-RED-IT parameters. The GUTS-RED-IT parameter "median of distribution of thresholds" (m_w ) showed a strong positive relation with the Ca²⁺ , Mg²⁺ , Na⁺ , and H⁺ ion activities. Therefore, m_w formed the base of the calibrated GUTS bioavailability model (GUTS-BLM), which predicted 30-d x% lethal concentration (LCx) values within a 2-fold error. The GUTS-BLM was combined with an IBM, inSTREAM-Gen, into a GUTS-BLM-IBM. Assuming that juvenile survival was the only effect of Cu and Zn exposure, population-level effect concentrations were predicted to be 1.3 to 6.2 times higher than 30-d laboratory LCx values, with the larger differences being associated with higher interindividual variation of metal sensitivity. The proposed GUTS-BLM-IBM model can provide insight into metal bioavailability and effects at the population level and could be further improved by incorporating sublethal effects of Cu and Zn. Environ Toxicol Chem 2021;40:2764-2780. © 2021 SETAC.

Dynamic interaction processes of rare earth metal mixtures in terrestrial organisms interpreted by toxicokinetic and toxicodynamic model

Despite the progress in explanation of mixture toxicity of rare earth elements (REEs), a large knowledge gap still exists in interpreting their mixed effects from a dynamic perspective. Here, we investigated the effects of La-Ce mixtures in Enchytraeus crypticus at different exposure times. The single and mixture toxicity of La and Ce increased with time, as reflected by the reduced LC50/MT50 values. With concentration addition as the reference model, the interactions between La and Ce were quantified by MIXTOX modelling tool, showing a time-dependent pattern with antagonistic effect after 1 and 2 d but additive effects afterwards. The dynamic accumulation and toxicity of La/Ce in organisms exposed to REE mixtures was fitted using a process-based toxicokinetic and toxicodynamic (TK-TD) model to unravel how the elements interacted. Generally, the estimated uptake, elimination, and damage rate constants of La/Ce declined with increasing level of each other, suggesting inhibited uptake and subsequently reduced toxicity of La/Ce due to competition effect. The interplay of La and Ce in TK and TD processes seemed responsible for the observed antagonism. Our study showed that mixture toxicity and interaction of REEs are time-dependent processes and application of TK-TD model may provide more insight into this dynamic effect.

Ozonolysis of Permethrin in the Atmosphere: Mechanism, Kinetics, and Evaluation of Toxicity

Pyrethroid, a pesticide widely used worldwide, could mimic, block, or synergize the effects of endogenous hormones in humans or mammals after entering into the atmosphere and after being sprayed and applied in large quantities. This research aims to study the mechanism, kinetics, and eco-toxicity evaluation of the ozonolysis of permethrin (PER)-one of the typical pyrethroid (type I) pesticides. Existing experimental studies only predicted that ozonolysis of PER could generate a cycloperoxy analogue of PER (IM13-1-11), and the reaction mechanism has not yet been completed. To make up for the lack of experimental results, the 13 primary reaction pathways of PER and ozone, as well as the subsequent reactions of Criegee intermediates with small molecules such as NOx, COx, SO2, and O2, have been studied to propose new reaction paths by quantum chemical calculations in this work. We calculated the total reaction rate constant of PER and ozone at 298 K and 1 atm based on the calculated thermodynamic data and the transition state theory (TST), which was compared with the experimental values to prove the reliability of our results. Based on the quantitative structure and activity relationship, we predicted the acute and chronic toxicity of PER and its products of ozonolysis to three representative organisms-fish, daphnia, and green algae to avoid animal experiments. The results show that ozonolysis products of PER are still extremely harmful to the environment and should be taken seriously, although the products have less toxicity than PER.

Role of Toxicokinetic and Toxicodynamic Parameters in Explaining the Sensitivity of Zebrafish Larvae to Four Metals

Given the persistence and toxic potencies of metal contaminants in ecosystems, animals, and human beings, they are considered to be hazardous global pollutants. While the lethality of metal toxicities (e.g., LC₅₀) can significantly vary, even within the same species, the underlying mechanisms are less well-understood. In this study, we developed a subcellular two-compartment toxicokinetic-toxicodynamic (TK-TD) model for zebrafish larvae when exposed to four metals (cadmium, lead, copper, and zinc) to reveal whether differences in metal toxicity (LC₅₀ values) were dominated by the TK or TD processes. Results showed that the subcellular TK and TD parameters of the four metals were significantly different, and the bioconcentration factor (BCF) value of copper was higher than those of the other metals. We also found that the TD parameter internal threshold concentration (C_IT) was significantly positively correlated to the LC₅₀ values (R² = 0.7), suggesting a dominant role of TD processes in metal toxicity. Furthermore, the combined parameter C_IT/BCF for a metal-sensitive fraction (BCF_MSF), which linked exposure to effects through the TK-TD approach, explained up to 89% of the variation in toxicity to the four metals. The present study suggests that the observed variation in toxicity of these four metals was mainly determined by TD processes but that TK processes should not be ignored, especially for copper.

Off-Target Stoichiometric Binding Identified from Toxicogenomics Explains Why Some Species Are More Sensitive than Others to a Widely Used Neonicotinoid

Neonicotinoids are currently licensed for use in 120 countries, making accurate nontarget species sensitivity predictions critical. Unfortunately, such predictions are fraught with uncertainty, as sensitivity is extrapolated from only a few test species and neonicotinoid sensitivities can differ greatly between closely related taxa. Combining classical toxicology with de novo toxicogenomics could greatly improve sensitivity predictions and identify unexpectedly susceptible species. We show that there is a >30-fold differential species sensitivity (DSS) for the neonicotinoid imidacloprid between five earthworm species, a critical nontarget taxon. This variation could not be explained by differential toxicokinetics. Furthermore, comparing key motif expression in subunit genes of the classical nicotinic acetylcholine receptor (nAChR) target predicts only minor differences in the ligand binding domains (LBDs). In contrast, predicted dissimilarities in LBDs do occur in the highly expressed but nonclassical targets, acetylcholine binding proteins (AChBPs). Critically, the predicted AChBP divergence is capable of explaining DSS. We propose that high expression levels of putative nonsynaptic AChBPs with high imidacloprid affinities reduce imidacloprid binding to critical nAChRs involved in vital synaptic neurotransmission. This study provides a clear example of how pragmatic interrogation of key motif expression in complex multisubunit receptors can predict observed DSS, thereby informing sensitivity predictions for essential nontarget species.

Can Organophosphates and Carbamates Cause Synergisms by Inhibiting Esterases Responsible for Biotransformation of Pyrethroids?

Hydrolysis catalyzed by general esterases (GEs) is the most efficient route for hydrolyzation of pyrethroid insecticides. Organophosphate (OP) and carbamate (CB) insecticides are known to inhibit GEs in addition to acetylcholinesterase (AChE), which is their main target. We hypothesize that synergies can be induced by OPs and CBs when mixed with pyrethroids, due to their inhibition of GE-dependent detoxification of pyrethroids. To test this hypothesis, we conducted mixture toxicity experiments with Daphnia magna using α-cypermethrin (α-cyp) in combination with the noninsecticidal OP tetraisopropyl pyrophosphoramide (iso-OMPA) and five AChE inhibitors diazinon, chlorpyrifos, chlorfenviphos, parathion, and aldicarb. In addition, the in vivo GE activity inhibition was measured for all compounds. Up to 10-fold synergy was found between α-cyp and iso-OMPA, and the degree of synergy correlated linearly with the inhibition of the GE activity. No synergy, however, was found in any of the insecticide mixtures nor was the GE activity inhibited within the nonlethal concentration range tested. It was concluded that the effect of the insecticides on AChE occurred at lower concentrations than their effect on GEs, making the daphnids become immobilized before any synergistic effects on mortality could be observed. The implications of the findings are discussed from a risk assessment perspective.

Effect Modeling Quantifies the Difference Between the Toxicity of Average Pesticide Concentrations and Time-Variable Exposures from Water Quality Monitoring

Synthetic chemicals are frequently detected in water bodies, and their concentrations vary over time. Water monitoring programs typically employ either a sequence of grab samples or continuous sampling, followed by chemical analysis. Continuous time-proportional sampling yields the time-weighted average concentration, which is taken as proxy for the real, time-variable exposure. However, we do not know how much the toxicity of the average concentration differs from the toxicity of the corresponding fluctuating exposure profile. We used toxicokinetic-toxicodynamic models (invertebrates, fish) and population growth models (algae, duckweed) to calculate the margin of safety in moving time windows across measured aquatic concentration time series (7 pesticides) in 5 streams. A longer sampling period (14 d) for time-proportional sampling leads to more deviations from the real chemical stress than shorter sampling durations (3 d). The associated error is a factor of 4 or less in the margin of safety value toward underestimating and an error of factor 9 toward overestimating chemical stress in the most toxic time windows. Under- and overestimations occur with approximate equal frequency and are very small compared with the overall variation, which ranged from 0.027 to 2.4 × 10¹⁰ (margin of safety values). We conclude that continuous, time-proportional sampling for a period of 3 and 14 d for acute and chronic assessment, respectively, yields sufficiently accurate average concentrations to assess ecotoxicological effects. Environ Toxicol Chem 2020;39:2158-2168. © 2020 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.