Body size-dependent toxicokinetics and toxicodynamics could explain intra- and interspecies variability in sensitivity

Diterpenoid and C20 diterpenoid alkaloid as a potent inhibitor of SARS-CoV-2 main protease (Mpro): from Piper barberi Gamble, an endemic and endangered species of Southern Western Ghats

The present study investigated the phytochemicals and in silico anti-nCoV properties of Piper barberi, an endangered and endemic species of Southern Western Ghats. Using conventional soxhlet extraction method, the leaf and stem were extracted separately with methanol (PBLM and PBSM). The bioactive compounds from the extracts were identified using HR-LCMS/MS-qTOF analysis. These compounds were subjected to various in silico analyses to identify potential drug candidates against nCoV. The HR LCMS/MS analysis of PBLM and PBSM revealed the presence of phenols, flavonoids, alkaloids, and terpenoids in it and this is the first report of the phytoconstituents present in the species P. barberi. All the identified bioactive compounds were subjected to predict ADMET. Out of 49 identified compounds, only 31 passed drug-likeness properties and toxicity tests. Molecular interaction studies were conducted using the AutoDockTools 4.2.6., which showed that only 13 compounds exhibited acceptable binding affinity with the nCoV target Mpro. Structural stability and binding free energy analyses of the five compounds with the higher binding affinity indicated that the bioactive compounds Hetisine and Ajaconine are stable with both hydrogen bonds and hydrophobic interactions. Hetisine shows stable binding among these two compounds with two hydrogen bond interactions with the crucial catalytic dyad residue (His41). Thus, this study concludes that these compounds might potentially be used as an alternative drug candidate for managing nCoV. However, further experimental validation, including in vitro and in vivo assays, is required to substantiate the results.Communicated by Ramaswamy H. Sarma.

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.

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

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

Comparability of comparative toxicity: insect sensitivity to imidacloprid reveals huge variations across species but also within species

Pesticides have been identified as major drivers of insect biodiversity loss. Thus, the study of their effects on non-pest insect species has attracted a lot of attention in recent decades. In general toxicology, the 'gold standard' to assess the toxicity of a substance is to measure mass-specific LD50 (i.e. median lethal dose per unit body mass). In entomology, reviews attempting to compare these data across all available studies are lacking. To fill this gap in knowledge, we performed a systematic review of the lethality of imidacloprid for adult insects. Imidacloprid is possibly the most extensively studied insecticide in recent times, yet we found that little is comparable across studies, owing to both methodological divergence and missing estimates of body mass. By accounting for body mass whenever possible, we show how imidacloprid sensitivity spans across an apparent range of approximately six orders of magnitude across insect species. Very high variability within species can also be observed owing to differences in exposure methods and observation time. We suggest that a more comparable and comprehensive approach has both biological and economic relevance. Ultimately, this would help to identify differences that could direct research towards preventing non-target species from being negatively affected.

A Toxicokinetic-Toxicodynamic Modeling Workflow Assessing the Quality of Input Mortality Data

Toxicokinetic-toxicodynamic (TKTD) models simulate organismal uptake and elimination of a substance (TK) and its effects on the organism (TD). The Reduced General Unified Threshold model of Survival (GUTS-RED) is a TKTD modeling framework that is well established for aquatic risk assessment to simulate effects on survival. The TKTD models are applied in three steps: parameterization based on experimental data (calibration), comparing predictions with independent data (validation), and prediction of endpoints under environmental scenarios. Despite a clear understanding of the sensitivity of GUTS-RED predictions to the model parameters, the influence of the input data on the quality of GUTS-RED calibration and validation has not been systematically explored. We analyzed the performance of GUTS-RED calibration and validation based on a unique, comprehensive data set, covering different types of substances, exposure patterns, and aquatic animal species taxa that are regularly used for risk assessment of plant protection products. We developed a software code to automatically calibrate and validate GUTS-RED against survival measurements from 59 toxicity tests and to calculate selected model evaluation metrics. To assess whether specific survival data sets were better suited for calibration or validation, we applied a design in which all possible combinations of studies for the same species-substance combination are used for calibration and validation. We found that uncertainty of calibrated parameters was lower when the full range of effects (i.e., from high survival to high mortality) was covered by input data. Increasing the number of toxicity studies used for calibration further decreased parameter uncertainty. Including data from both acute and chronic studies as well as studies under pulsed and constant exposure in model calibrations improved model predictions on different types of validation data. Using our results, we derived a workflow, including recommendations for the sequence of modeling steps from the selection of input data to a final judgment on the suitability of GUTS-RED for the data set. Environ Toxicol Chem 2024;43:197-210. © 2023 Bayer AG and The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Treated and highly diluted, but wastewater still impacts diversity and energy fluxes of freshwater food webs

Wastewater treatment plants (WWTPs) have greatly improved water quality globally. However, treated effluents still contain a complex cocktail of pollutants whose environmental effects might go unnoticed, masked by additional stressors in the receiving waters or by spatiotemporal variability. We conducted a BACI (Before-After/Control-Impact) ecosystem manipulation experiment, where we diverted part of the effluent of a large tertiary WWTP into a small, unpolluted stream to assess the effects of a well-treated and highly diluted effluent on riverine diversity and food web dynamics. We sampled basal food resources, benthic invertebrates and fish to search for changes on the structure and energy transfer of the food web with the effluent. Although effluent toxicity was low, it reduced diversity, increased primary production and herbivory, and reduced energy fluxes associated to terrestrial inputs. Altogether, the effluent decreased total energy fluxes in stream food webs, showing that treated wastewater can lead to important ecosystem-level changes, affecting the structure and functioning of stream communities even at high dilution rates. The present study shows that current procedures to treat wastewater can still affect freshwater ecosystems and highlights the need for further efforts to treat polluted waters to conserve aquatic food webs.

Chemical Mixtures and Multiple Stressors: Same but Different?

Ecosystems are strongly influenced by multiple anthropogenic stressors, including a wide range of chemicals and their mixtures. Studies on the effects of multiple stressors have largely focussed on nonchemical stressors, whereas studies on chemical mixtures have largely ignored other stressors. However, both research areas face similar challenges and require similar tools and methods to predict the joint effects of chemicals or nonchemical stressors, and frameworks to integrate multiple chemical and nonchemical stressors are missing. We provide an overview of the research paradigms, tools, and methods commonly used in multiple stressor and chemical mixture research and discuss potential domains of cross-fertilization and joint challenges. First, we compare the general paradigms of ecotoxicology and (applied) ecology to explain the historical divide. Subsequently, we compare methods and approaches for the identification of interactions, stressor characterization, and designing experiments. We suggest that both multiple stressor and chemical mixture research are too focused on interactions and would benefit from integration regarding null model selection. Stressor characterization is typically more costly for chemical mixtures. While for chemical mixtures comprehensive classification systems at suborganismal level have been developed, recent classification systems for multiple stressors account for environmental context. Both research areas suffer from rather simplified experimental designs that focus on only a limited number of stressors, chemicals, and treatments. We discuss concepts that can guide more realistic designs capturing spatiotemporal stressor dynamics. We suggest that process-based and data-driven models are particularly promising to tackle the challenge of prediction of effects of chemical mixtures and nonchemical stressors on (meta-)communities and (meta-)food webs. We propose a framework to integrate the assessment of effects for multiple stressors and chemical mixtures. Environ Toxicol Chem 2023;42:1915-1936. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

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.

Three perspectives on the prediction of chemical effects in ecosystems

The increasing production, use and emission of synthetic chemicals into the environment represents a major driver of global change. The large number of synthetic chemicals, limited knowledge on exposure patterns and effects in organisms and their interaction with other global change drivers hamper the prediction of effects in ecosystems. However, recent advances in biomolecular and computational methods are promising to improve our capacity for prediction. We delineate three idealised perspectives for the prediction of chemical effects: the suborganismal, organismal and ecological perspective, which are currently largely separated. Each of the outlined perspectives includes essential and complementary theories and tools for prediction but captures only part of the phenomenon of chemical effects. Links between the perspectives may foster predictive modelling of chemical effects in ecosystems and extrapolation between species. A major challenge for the linkage is the lack of data sets simultaneously covering different levels of biological organisation (here referred to as biological levels) as well as varying temporal and spatial scales. Synthesising the three perspectives, some central aspects and associated types of data seem particularly necessary to improve prediction. First, suborganism- and organism-level responses to chemicals need to be recorded and tested for relationships with chemical groups and organism traits. Second, metrics that are measurable at many biological levels, such as energy, need to be scrutinised for their potential to integrate across levels. Third, experimental data on the simultaneous response over multiple biological levels and spatiotemporal scales are required. These could be collected in nested and interconnected micro- and mesocosm experiments. Lastly, prioritisation of processes involved in the prediction framework needs to find a balance between simplification and capturing the essential complexity of a system. For example, in some cases, eco-evolutionary dynamics and interactions may need stronger consideration. Prediction needs to move from a static to a real-world eco-evolutionary view.

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.

Size- and sex-related sensitivity differences of aquatic crustaceans to imidacloprid

Field collected aquatic invertebrates are often used as test organisms in the refinement of the standard Tier 1 risk assessment of various pollutants. This approach can provide insights into the effects of pollutants on the natural environment. However, researchers often pragmatically select test organisms of a specific sex and/or size, which may not represent the sensitivity of the whole population. To investigate such intraspecies sensitivity differences, we performed standard acute toxicity and toxicokinetic tests with different size classes and sex of Gammarus pulex and Asellus aquaticus. Furthermore, toxicokinetics and toxicodynamics models were used to understand the mechanism of the intraspecies sensitivity differences. We used neonates, juveniles and male and female adults in separate dedicated experiments, in which we exposed the animals to imidacloprid and its bioactive metabolite, imidacloprid-olefin. For both species, we found that neonates were the most sensitive group. For G. pulex, the sensitivity decreased linearly with size, which can be explained by the size-related uptake rate constant in the toxicokinetic process and size-related threshold value in the toxicodynamic process. For A. aquaticus, female adults were least sensitive to imidacloprid, which could be explained by a low internal biotransformation of imidacloprid to imidacloprid-olefin. Besides, imidacloprid-olefin was more toxic than imidacloprid to A. aquaticus, with differences being 8.4 times for females and 2.7 times for males. In conclusion, we established size-related sensitivity differences for G. pulex and sex-related sensitivity for A. aquaticus, and intraspecies differences can be explained by both toxicokinetic and toxicodynamic processes. Our findings suggest that to protect populations in the field, we should consider the size and sex of focal organisms and that a pragmatic selection of test organisms of equal size and/or sex can underestimate the sensitivities of populations in the field.

Moving beyond Risk Quotients: Advancing Ecological Risk Assessment to Reflect Better, More Robust and Relevant Methods

Under standard guidance for conducting Ecological Risk Assessments (ERAs), the risks of chemical exposure to diverse organisms are most often based on deterministic point estimates evaluated against safety-factor-based levels of concern (LOCs). While the science and guidance for mechanistic effect models (e.g., demographic, population, and agent-based) have long been demonstrated to provide more ecologically relevant effect endpoints upon which risk can be evaluated, their application in ERAs has been limited, particularly in the US. This special issue highlights the state of the science in effect modeling for ERAs through demonstrated application of the recently published Population modeling Guidance, Use, Interpretation, and Development for ERA (Pop-GUIDE). We introduce this issue with a perspective on why it is critical to move past the current application of deterministic endpoints and LOCs. We demonstrate how the current, widely used approaches contain extensive uncertainty that could be reduced considerably by applying models that account for species life histories and other important endogenous and exogenous factors critical to species sustainability. We emphasize that it is long past time to incorporate better, more robust, and ecologically relevant effect models into ERAs, particularly for chronic risk determination. The papers in this special issue demonstrate how mechanistic models that follow Pop-GUIDE better inform ERAs compared to the current standard practice.

Alternative animal models in predictive toxicology

Toxicity testing relies heavily on animals, especially rodents as part of the non-clinical laboratory testing of substances. However, the use of mammalians and the number of animals employed in research has become a concern for institutional ethics committees. Toxicity testing involving rodents and other mammals is laborious and costly. Alternatively, non-rodent models are used as replacement, as they have less ethical considerations and are cost-effective. Of the many alternative models that can be used as replacement models, which ones can be used in predictive toxicology? What is the correlation between these models and rodents? Are there standardized protocols governing the toxicity testing of these commonly used predictive models? This review outlines the common alternative animal models for predictive toxicology to address the importance of these models, the challenges, and their standard testing protocols.

Direct and Delayed Mortality of Ceriodaphnia dubia and Rainbow Trout Following Time-Varying Acute Exposures to Zinc

The potential for delayed mortality following short-term episodic pollution events was evaluated by exposing cladocerans (Ceriodaphnia dubia) and rainbow trout (Oncorhynchus mykiss) to zinc (Zn) in various 1- to 48-h and 1- to 96-h exposures, respectively, followed by transferring the exposed organisms to clean water for up to 47 h for C. dubia and up to 95 h for trout for additional observation. For C. dubia, 1-h exposures of up to 3790 µg Zn/L never resulted in mortality during the actual Zn exposures, but by 48 h, a 1-h exposure to 114 µg/L, a concentration similar to the present US national water quality acute criterion for the test water conditions, ultimately killed 70% of C. dubia. With C. dubia, the speed of action of Zn toxicity was faster for intermediate concentrations than for the highest concentrations tested. For rainbow trout, pronounced delayed mortalities by 96 h only occurred following ≥8-h exposures. For both species, ultimate mortalities from Zn exposures ≤8 h mostly presented as delayed mortalities, whereas for exposures ≥24 h, almost all ultimate mortalities presented during the actual exposure periods. With Zn, risks of delayed mortality following exposures to all concentrations tested were much greater for the more sensitive, small-bodied invertebrate (C. dubia) than for the less sensitive, larger-bodied fish (rainbow trout). These results, along with previous studies, show that delayed mortality is an important consideration in evaluating risks to aquatic organisms from brief, episodic exposures to some substances. Environ Toxicol Chem 2021;40:2484-2498. © 2021 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Neonicotinoid Exposures that Stimulate Predatory Stink Bug, Podisus maculiventris (Hemiptera: Pentatomidae), Reproduction Do Not Inhibit Its Behavior

Exposure to sublethal amounts of pesticide can compromise life-history traits and behavior of natural enemies thereby reducing their effectiveness as predators. However, sublethal exposures to pesticides and other stressors may also stimulate insects, a dose-response phenomenon known as hormesis. We previously reported stimulatory effects on reproduction in the beneficial insect predator Podisus maculiventris (Say) (Hemiptera: Pentatomidae) following exposure to sublethal concentrations of imidacloprid. Here we examined whether these same treatments stimulated behavior and/or predation of P. maculiventris. Stimulation of some behaviors occurred at a reproductively hormetic concentration and two additional sublethal concentrations, depending upon bioassay design and sex. We observed no substantial inhibition of behavior or predation at a reproductively hormetic concentration, demonstrating that reproductive fitness in P. maculiventris may be stimulated without compromising behaviors important in its effectiveness as a natural enemy.

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.

Disentangling Mechanisms Behind Chronic Lethality through Toxicokinetic-Toxicodynamic Modeling

Ecotoxicological profiles of the 3 insecticides imidacloprid, thiacloprid, and flupyradifurone in terms of acute and chronic effects were analyzed in Chironomus riparius. Toxicokinetic-toxicodynamic modeling revealed that chironomids would die from starvation as a result of prolonged feeding inhibition under chronic exposures. The starvation effect is an indirect cause for mortality, which, for the neonicotinoids, adds to the direct/acute mortality, although the results suggests that this additional effect is not relevant for flupyradifurone. Environ Toxicol Chem 2021;40:1706-1712. © 2021 Bayer Inc. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Chemicals with increasingly complex modes of action result in greater variation in sensitivity between earthworm species

The scale of variation in species sensitivity to toxicants has been theoretically linked to mode of action. Specifically, it has been proposed there will be greater variations for chemicals with a putative specific biological target than for toxicants with a non-specific narcotic mechanism. Here we test the hypothesis that mode of action is related to variation in sensitivity in a specifically designed experiment for species from a single ecologically important terrestrial taxa, namely earthworms. Earthworm toxicity tests were conducted with five species for four chemicals, providing a series of increasingly complex modes of action: a putative narcotic polycyclic aromatic hydrocarbon (fluoranthene), and three insecticides (chlorpyrifos, cypermethrin, imidacloprid) with known neuronal receptor targets. Across all the chemicals, the standard epigeic test species Eisenia fetida and Lumbricus rubellus, were generally among the two least sensitive, while the endogenic Aporrectodea caliginosa and Megascolecidae Amynthas gracilis were generally more sensitive (never being among the two least sensitive species). This indicates a potential for bias in the earthworm ecotoxicology literature, which is dominated by studies in epigeic Lumbricidae, but contains few endogeic or Megascolecidae data. Results confirmed the lowest range of variation in sensitivities for effects on reproduction was for fluoranthene (2.5 fold). All insecticides showed greater variation for species sensitivity (cypermethrin: 7.5 fold, chlorpyrifos: 10.3 fold, imidacloprid: 31.5 fold) consistent with the specific mechanisms of the pesticides. Difference in toxicodynamics, based on mode of action specificity and receptor complexity was reflected in the magnitude of sensitivity variation. However, measurements of tissue concentrations also indicated the potential importance of toxicokinetics in explaining species sensitivity variations for chlorpyrifos and cypermethrin.

A biological characteristic extrapolation of compound toxicity for different developmental stage species with toxicokinetic-toxicodynamic model

Intraspecific difference in toxicity brings uncertainty to ecological risk assessment (ERA) and water quality criteria (WQC) of chemicals. Here, we compared intraspecies sensitivity to toxicants for Mesocyclops leuckarti of which toxicity data was obtained from published literatures, and zebrafish Danio rerio of which toxicity data was done in this study). Due to the internal concentration of chemicals not measured, simplified toxicokinetic-toxicodynamic (TK-TD) models were used, and we investigated whether TK-TD parameters estimated by Bayesian method might represent the differences in sensitivity between life-stages of 2 species. The results demonstrated that the difference in TK-TD parameters (background mortality m₀, no effect concentration NEC, the killing rate k_s, and the dominant rate k_d) could represent the toxicity difference between life-stages of individual species. The TK-TD model could predict toxicity in individual species (Cyprinus carpio L., Enchytraeus crypticus, Folsomia candida, Hyalella Azteca) exposed to different chemical concentrations and successfully extrapolate toxicity between different life stages of Mesocyclops leuckarti and Danio rerio by scaling several TK-TD parameters. The modified TK-TD model on the extrapolation toxicity of chemicals between life stages for species could be useful for the ERA and for deriving and revising WQC for chemicals.

Sex-dependent locomotion and physiological responses shape the insecticidal susceptibility of parasitoid wasps

The adaptive fitness of insect species can be shaped by how males and females respond, both physiologically and behaviorally, to environmental challenges, such as pesticide exposure. In parasitoid wasps, most toxicological investigations focus only on female responses (e.g., survival and especially parasitism abilities), leaving the male contributions to adaptive fitness (survival, locomotion, mate search) poorly investigated. Here, we evaluated the toxicity of the spinosyn insecticide spinosad against the South American fruit fly, Anastrepha fraterculus, and we used the parasitoid wasp Diachasmimorpha longicaudata (Ashmead) to evaluate whether sex-linked locomotory and physiological responses would influence the susceptibility of these organisms to spinosad. Our results revealed that D. longicaudata males were significantly more susceptible (median lethal time (LT₅₀) = 24 h) to spinosad than D. longicaudata females (LT₅₀ = 120 h), which may reflect the differences in their locomotory and physiological (e.g., respiratory) responses to mitigate insecticide exposure. Compared to D. longicaudata females, male wasps were lighter (P < 0.001), walked for longer distances (P < 0.001) and periods (P < 0.001), and exhibited higher sensilla densities in their tarsi (P = 0.008), which may facilitate their intoxication with the insecticide. These findings indicate that male parasitoids should not be exempt from insecticide selectivity tests, as these organisms can be significantly more affected by such environmental challenges than their female conspecifics.

Towards a spatiotemporally explicit toxicokinetic-toxicodynamic model for earthworm toxicity

The aim of the environmental risk assessment of chemicals is the prevention of unacceptable adverse effects on the environment. Therefore, the risk assessment for in-soil organisms, such as earthworms, is based on two key elements: the exposure assessment and the effect assessment. In the current risk assessment scheme, these two elements are not linked. While for the exposure assessment, advanced exposure models can take the spatial and temporal scale of substances into account, the effect assessment in the lower tiers considers only a limited temporal and spatial variability. However, for soil organisms, such as earthworms, those scales play a significant role as species move through the soil in response to environmental factors. To overcome this gap, we propose a conceptual integration of pesticide exposure, ecology, and toxicological effects on earthworms using a modular modeling approach. An essential part of this modular approach is the environment module, which utilizes exposure models to provide spatially and temporally explicit information on environmental variables (e.g., temperature, moisture, organic matter content) and chemical concentrations. The behavior module uses this information and simulates the feeding and movement of different earthworm species using a trait-based approach. The resulting exposure can be processed by a toxicokinetic-toxicodynamic (TKTD) module. TKTD models are particularly suitable to make effect predictions for time-variable exposure situations as they include the processes of uptake, elimination, internal distribution, and biotransformation of chemicals and link the internal concentration to an effect at the organism level. The population module incorporates existing population models of different earthworm species. The modular approach is illustrated using a case study with an insecticide. Our results emphasize that using a modular model approach will facilitate the integration of exposure and effects and thus enhance the risk assessment of soil organisms.

Mechanistic Effect Modeling Approach for the Extrapolation of Species Sensitivity

In the higher-tier environmental risk assessment of chemicals, species sensitivity distributions (SSDs) are used to statistically describe differences in sensitivity between species and derive community level endpoints. SSDs are usually based on the results from short-term laboratory experiments performed under constant environmental conditions. However, different species may be kept at different "optimal" temperatures, which influence their apparent sensitivity and thus the derivation of endpoints. Also, the extrapolation capacity of SSDs is largely limited to the tested species and conditions. Time-variable exposures and effects at higher levels of biological organization, including biological interactions, are not considered. The quantitative effect prediction at higher tiers would ultimately require the extrapolation of toxicokinetics and toxicodynamics to untested species and the involvement of population and community modeling. In this regard, we tested a toxicokinetic-toxicodynamic modeling approach to mechanistically consider and correct endpoints for ambient temperature and demonstrate the significance for SSDs. We explored correlations in toxicokinetic-toxicodynamic model parameters which would allow for the extrapolation of sensitivities to untested species. Finally, we illustrate the applicability of the approach for higher level effect predictions using an individual-based model. Our results suggest that mechanistic effect modeling approaches can reduce the uncertainties in higher tier effect assessments related to knowledge gaps.

Toxicokinetic-toxicodynamic modeling of cadmium and lead toxicity to larvae and adult zebrafish

Toxicity of hazard materials to organism is different between larvae and adult zebrafish. However, this different effect was seldom considered in toxicological modeling. Here, we measured Cd and Pb toxicity for larvae and adult zebrafish (Danio rerio) and assessed whether metal toxicity can be better simulated by the one-compartment or two-compartment toxicokinetic (TK) and toxicodynamic (TD) models with assumption of stochastic death (SD) and individual tolerance (IT), respectively. Results showed that, for larvae, the one-compartment model generally fitted the observed accumulation and survival better than two-compartment model. In contrast, for adult, the two-compartment model simulation satisfied the observed accumulation and survival better than one-compartment model. In addition, both the SD and the IT models generally described the Cd or Pb toxicity well, although the IT model predictions were slightly better than the SD model in adult fish, the opposite phenomenon was observed in larvae. Our results suggested that variations in both TK and TD parameters might be needed to quantify the toxicity sensitivity in larvae and adult zebrafish, and accounting these variations in mechanistic toxicological effect models (e.g. TK-TD) will allow more accurate predictions of hazard materials effects to organisms.

Short-term to long-term extrapolation of lethal effects of an herbicide on the marine mysid shrimp Americamysis Bahia by use of the General Unified Threshold Model of Survival (GUTS)

Risk assessments for plant protection products and their active ingredients that are based on standard laboratory tests performed under constant exposure conditions may result in an overestimation of risks because exposure in the environment is often characterized by a few short peaks. Here, the General Unified Threshold Model of Survival (GUTS) was used to conduct a refined risk assessment for the herbicide tembotrione and its effects on the marine invertebrate Americamysis bahia, for which the standard chronic effect assessment failed. The GUTS model was first calibrated with time-to-effect and concentration-response data from 2 independent acute experiments with A. bahia. Model parameters for both toxicodynamic assumptions of stochastic death (SD) and individual tolerance (IT) were estimated with the reduced GUTS model (GUTS-RED) using the scaled internal concentration as a dose metric. Both the calibrated GUTS-RED-SD and GUTS-RED-IT models described survival dynamics well. Model validation using datasets of 2 independent chronic tests yielded robust predictions of long-term toxicity of tembotrione on A. bahia, with GUTS-RED-IT being more reliable than GUTS-RED-SD. The validated model was subsequently used to predict survival from time-variable exposure profiles, as derived from the FOrum for Co-ordination of pesticide fate models and their USe (FOCUS). Because ecotoxicological independence of peaks had not been empirically verified, the link between exposure and effects was assessed with complete exposure profiles. Effect thresholds resulting from different peak exposure concentrations and durations were determined with GUTS and directly compared with the exposure concentrations from the FOCUS surface water scenarios. The derived values were higher than the predicted FOCUS critical concentrations. Additionally, comparing the areas under the curve (AUCs) derived with GUTS for multiple peak exposure profiles to those from FOCUS revealed significant additional safety margins, demonstrating that only unrealistically high exposure concentrations would produce significant effects. In conclusion, no unacceptable effects of tembotrione on aquatic invertebrates under realistic environmental exposure conditions are expected. Integr Environ Assess Manag 2019;15:29-39. © 2018 SETAC.

Radiosensitivity and transgenerational effects in non-human species

The ALLIANCE working group on effects of ionising radiation on wildlife brings together European researchers to work on the topics of radiosensitivity and transgenerational effects in non-human biota. Differences in radiation sensitivity across species and phyla are poorly understood, but have important implications for understanding the overall effects of radiation and for radiation protection; for example, sensitive species may require special attention in monitoring and radiation protection, and differences in sensitivity between species also lead to overall effects at higher levels (community, ecosystem), since interactions between species can be altered. Hence, understanding the mechanisms of interspecies radiation sensitivity differences may help to clarify mechanisms underpinning intraspecies variation. Differences in sensitivity may only be revealed when organisms are exposed to ionising radiation over several generations. This issue of potential long-term or hereditary effects for both humans and wildlife exposed to low doses of ionising radiation is a major concern. Animal and plant studies suggest that gamma irradiation can lead to observable effects in the F1 generation that are not attributable to inheritance of a rare stable DNA mutation. Several studies have provided evidence of an increase in genomic instability detected in germ or somatic cells of F1 organisms from exposed F0 organisms. This can lead to induced radiosensitivity, and can result in phenotypic effects or lead to reproductive effects and teratogenesis. In particular, studies have been conducted to understand the possible role of epigenetic modifications, such as DNA methylation, histone modifications, or expression of non-coding RNAs in radiosensitivity, as well as in adaptation effects. As such, research using biological models in which the relative contribution of genetic and epigenetic processes can be elucidated is highly valuable.

Aquatic ecotoxicological models and their applicability in Arctic regions

Dose-response modeling is one of the most important steps of ecological risk assessment. It requires concentration-effects relationships for the species under consideration. There are very limited studies and experimental data available for the Arctic aquatic species. Lack of toxicity data hinders obtaining dose-response relationships for lethal (LC50 values), sub-lethal and carcinogenic effects. Gaps in toxicity data could be filled using a variety of in-silico ecotoxicological methods. This paper reviews the suitability of such methods for the Arctic scenario. Mechanistic approaches like toxicokinetic and toxicodynamic analysis are found to be better suited for interspecies extrapolation than statistical methods, such as Quantitative Structure-Activity Relationships/Quantitative Structure Activity-Activity Relationship, ICE, and other empirical models, such as Haber's law and Ostwald's equation. A novel approach is proposed where the effects of the toxicant exposure are quantified based on the probability of cellular damage and metabolites interactions. This approach recommends modeling cellular damage using a toxicodynamic model and physiology or metabolites interactions using a toxicokinetic model. Together, these models provide more reliable estimates of toxicity in the Arctic aquatic species, which will assist in conducting ecological risk assessment of Arctic environment.

Daphniatox - Online monitoring of aquatic pollution and toxic substances

The microcrustacean Daphnia is sensitive to many toxic substances and can be cultured easily. The Daphniatox instrument is based on computerized image analysis tracking swimming organisms in real time. The software evaluates 14 endpoints including motility, swimming velocity, orientation with respect to light and gravity as well as cell form and size. The system determines movement vectors of a large number of organisms to warrant high statistical significance and calculates mean values as well as standard deviation. Tests with K dichromate show that the toxin inhibits motility (EC₅₀ 0.75 mg/L), swimming velocity (EC₅₀ 0.70 mg/L) and even causes a significant decrease in length (16% at 4 mg/L) and changes the form of the animals, This bioassay can be used to monitor the toxicity of a large number of dissolved pollutants and toxic substances such as arsenic, dichromate and persistent organic pollutants.

Modelling survival: exposure pattern, species sensitivity and uncertainty

The General Unified Threshold model for Survival (GUTS) integrates previously published toxicokinetic-toxicodynamic models and estimates survival with explicitly defined assumptions. Importantly, GUTS accounts for time-variable exposure to the stressor. We performed three studies to test the ability of GUTS to predict survival of aquatic organisms across different pesticide exposure patterns, time scales and species. Firstly, using synthetic data, we identified experimental data requirements which allow for the estimation of all parameters of the GUTS proper model. Secondly, we assessed how well GUTS, calibrated with short-term survival data of Gammarus pulex exposed to four pesticides, can forecast effects of longer-term pulsed exposures. Thirdly, we tested the ability of GUTS to estimate 14-day median effect concentrations of malathion for a range of species and use these estimates to build species sensitivity distributions for different exposure patterns. We find that GUTS adequately predicts survival across exposure patterns that vary over time. When toxicity is assessed for time-variable concentrations species may differ in their responses depending on the exposure profile. This can result in different species sensitivity rankings and safe levels. The interplay of exposure pattern and species sensitivity deserves systematic investigation in order to better understand how organisms respond to stress, including humans.

Demographic Toxicokinetic-Toxicodynamic Modeling of Lethal Effects

The aquatic effect assessment of chemicals is largely based on standardized measures of toxicity determined in short-term laboratory tests which are designed to reduce variability. For this purpose, uniform individuals of a species are kept under environmental and chemical exposure conditions which are as constant as possible. In nature, exposure often appears to be pulsed, effects might last longer than a few days, sensitivity might vary among different sized organisms and populations are usually size or age structured and are subject to demographic processes. To overcome this discrepancy, we tested toxicokinetic-toxicodynamic models of different complexities, including body size scaling approaches, for their ability to represent lethal effects observed for Daphnia magna exposed to triphenyltin. The consequences of the different toxicokinetic and toxicodynamic assumptions for population level responses to pulsed exposure are tested by means of an individual based model and are evaluated by confronting model predictions with population data for various pulsed exposure scenarios. We provide an example where increased model complexity reduces the uncertainty in model outputs. Furthermore, our results emphasize the importance of considering population demography in toxicokinetics and toxicodynamics for understanding and predicting potential chemical impacts at higher levels of biological organization.