Importance of toxicokinetics for interspecies variation in sensitivity to chemicals

Assessing ecological responses of exposure to the pyrethroid insecticide lambda-cyhalothrin in sub-tropical freshwater ecosystems

Pyrethroid insecticides are widely detected in aquatic ecosystems due to their extensive use in agriculture and horticulture, which could pose a potential risk to aquatic non-target organisms. While previous ecotoxicological studies have been conducted mainly with standard tests and local species under temperate conditions, scarce information is available on the effects of pyrethroid insecticides on communities and ecosystems under (sub-)tropical conditions. A single application of lambda-cyhalothrin at concentrations of 0, 9, 30, and 100 ng/L was evaluated in outdoor mesocosms under sub-tropical conditions. Lambda-cyhalothrin was found to dissipate rapidly in the water column, with only 11 % and 7 % of the remaining dose measured at 1 and 3 days after application, respectively. Lambda-cyhalothrin concentrations disappeared considerably faster from the water compartment compared to temperate conditions. Consistent decreases in abundance were observed for Lecane lunaris at the medium and higher treatments (NOEC = 9 ng/L) and at the highest treatment (NOEC = 30 ng/L) for Keratella tropica. On the contrary, two taxa belonging to Cladocera (i.e., Ceriodaphnia sp. and Diaphanosoma sp.) showed the most prominent increase in abundance related to the lambda-cyhalothrin treatments. At the community level, a consistent no observed effect concentrations (NOECs) of 9 ng/L could be calculated for the zooplankton community. A marginal significant overall treatment related effect was observed for the macroinvertebrate community. The results of species sensitivity distribution (SSD) analysis based on results of acute toxicity experiments conducted alongside the mesocosm experiment and obtained from the literature indicated that macroinvertebrates from temperate regions may be generally more sensitive than their counterparts in (sub-)tropical regions. Overall, these findings suggest that environmentally relevant concentrations of the pyrethroid insecticide lambda-cyhalothrin may lead to different ecological outcomes in freshwater ecosystems in the (sub-)tropics relative to temperate regions.

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.

Toxicokinetics Explain Differential Freshwater Ecotoxicity of Nanoencapsulated Imidacloprid Compared to Its Conventional Active Ingredient

Agricultural applications of nanotechnologies necessitate addressing safety concerns associated with nanopesticides, yet research has not adequately elucidated potential environmental risks between nanopesticides and their conventional counterparts. To address this gap, we investigated the risk of nanopesticides by comparing the ecotoxicity of nanoencapsulated imidacloprid (nano-IMI) with its active ingredient to nontarget freshwater organisms (embryonic Danio rerio, Daphnia magna, and Chironomus kiinensis). Nano-IMI elicited approximately 5 times higher toxicity than IMI to zebrafish embryos with and without chorion, while no significant difference was observed between the two invertebrates. Toxicokinetics further explained the differential toxicity patterns of the two IMI analogues. One-compartmental two-phase toxicokinetic modeling showed that nano-IMI exhibited significantly slower elimination and subsequently higher bioaccumulation potential than IMI in zebrafish embryos (dechorinated), while no disparity in toxicokinetics was observed between nano-IMI and IMI in D. magna and C. kiinensis. A two-compartmental toxicokinetic model successfully simulated the slow elimination of IMI from C. kiinensis and confirmed that both analogues of IMI reached toxicologically relevant targets at similar levels. Although nanopesticides exhibit comparable or elevated toxicity, future work is of utmost importance to properly understand the life cycle risks from production to end-of-life exposures, which helps establish optimal management measures before their widespread applications.

Transfer and bioaccumulation of pesticides in terrestrial arthropods and food webs: State of knowledge and perspectives for research

Arthropods represent an entry point for pesticide transfers in terrestrial food webs, and pesticide accumulation in upper chain organisms, such as predators can have cascading consequences on ecosystems. However, the mechanisms driving pesticide transfer and bioaccumulation in food webs remain poorly understood. Here we review the literature on pesticide transfers mediated by terrestrial arthropods in food webs. The transfer of pesticides and their potential for bioaccumulation and biomagnification are related to the chemical properties and toxicokinetic of the substances, the resistance and detoxification abilities of the contaminated organisms, as well as by their effects on organisms' life history traits. We further identify four critical areas in which knowledge gain would improve future predictions of pesticides impacts on terrestrial food webs. First, efforts should be made regarding the effects of co-formulants and pesticides mixtures that are currently understudied. Second, progress in the sensitivity of analytical methods would allow the detection of low concentrations of pesticides in small individual arthropods. Quantifying pesticides in arthropods preys, their predators, and arthropods or vertebrates at higher trophic level would bring crucial insights into the bioaccumulation and biomagnification potential of pesticides in real-world terrestrial food webs. Finally, quantifying the influence of the trophic structure and complexity of communities on the transfer of pesticides could address several important sources of variability in bioaccumulation and biomagnification across species and food webs. This narrative review will inspire future studies aiming to quantify pesticide transfers in terrestrial food webs to better capture their ecological consequences in natural and cultivated landscapes.

Quantitative differentiation of toxicity contributions and predicted global risk of fipronil and its transformation products to aquatic invertebrates

Biotransformation often alters chemical toxicity, yet its impacts on risk assessment are hardly quantified due to the challenges in acquiring internal exposure-based thresholds for chemicals that are readily metabolizable. Here, we integrated toxic unit and toxicokinetics to quantitatively assess toxicity contributions and potential risk of both parent compound and transformation products (TPs) to aquatic organisms, using fipronil (FIP) as a representative toxicant. In aquatic invertebrates Chironomus dilutus and Hyalella azteca, approximately 90 % of FIP was transformed to fipronil sulfone (SUL). FIP and SUL exhibited similar intrinsic toxicity to these organisms, which was contrary to conventional perception that SUL was more toxic than FIP. However, biotransformation was still important in risk assessment because the TP had 10-fold slower depuration rate than FIP. The amphipod H. azteca was found to be as sensitive to FIPs as the insect C. dilutus, which was previously considered ten times more sensitive based on external thresholds. This discrepancy has led to overlooking the toxicity of FIP to H. azteca in regional risk assessments. Lastly, we predicted the lethal risk of FIPs in global surface water. When using external thresholds for prediction, FIPs in 3.4 % of the water samples were lethally toxic to H. azteca, and the percentage of water samples at risk increased to 14 % when internal thresholds were used and SUF dominated the risk. This study presents an improved method for quantifying aquatic risk of readily metabolized toxicants. Our findings underscore the urgency of considering TPs in water quality assessments, especially for sensitive species that are at risk in the environment.

Cross-Species Extrapolation of Biological Data to Guide the Environmental Safety Assessment of Pharmaceuticals-The State of the Art and Future Priorities

The extrapolation of biological data across species is a key aspect of biomedical research and drug development. In this context, comparative biology considerations are applied with the goal of understanding human disease and guiding the development of effective and safe medicines. However, the widespread occurrence of pharmaceuticals in the environment and the need to assess the risk posed to wildlife have prompted a renewed interest in the extrapolation of pharmacological and toxicological data across the entire tree of life. To address this challenge, a biological "read-across" approach, based on the use of mammalian data to inform toxicity predictions in wildlife species, has been proposed as an effective way to streamline the environmental safety assessment of pharmaceuticals. Yet, how effective has this approach been, and are we any closer to being able to accurately predict environmental risk based on known human risk? We discuss the main theoretical and experimental advancements achieved in the last 10 years of research in this field. We propose that a better understanding of the functional conservation of drug targets across species and of the quantitative relationship between target modulation and adverse effects should be considered as future research priorities. This pharmacodynamic focus should be complemented with the application of higher-throughput experimental and computational approaches to accelerate the prediction of internal exposure dynamics. The translation of comparative (eco)toxicology research into real-world applications, however, relies on the (limited) availability of experts with the skill set needed to navigate the complexity of the problem; hence, we also call for synergistic multistakeholder efforts to support and strengthen comparative toxicology research and education at a global level. Environ Toxicol Chem 2024;43:513-525. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

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.

Elucidating the spatial distribution of organic contaminants and their biotransformation products in amphipod tissue by MALDI- and DESI-MS-imaging

The application of mass spectrometry imaging (MSI) is a promising tool to analyze the spatial distribution of organic contaminants in organisms and thereby improve the understanding of toxicokinetic and toxicodynamic processes. MSI is a common method in medical research but has been rarely applied in environmental science. In the present study, the suitability of MSI to assess the spatial distribution of organic contaminants and their biotransformation products (BTPs) in the aquatic invertebrate key species Gammarus pulex was studied. Gammarids were exposed to a mixture of common organic contaminants (carbamazepine, citalopram, cyprodinil, efavirenz, fluopyram and terbutryn). The distribution of the parent compounds and their BTPs in the organisms was analyzed by two MSI methods (MALDI- and DESI-HRMSI) after cryo-sectioning, and by LC-HRMS/MS after dissection into different organ compartments. The spatial distribution of contaminats in gammarid tissue could be successfully analyzed by the different analytical methods. The intestinal system was identified as the main site of biotransformation, possibly due to the presence of biotransforming enzymes. LC-HRMS/MS was more sensitive and provided higher confidence in BTP identification due to chromatographic separation and MS/MS. DESI was found to be the more sensitive MSI method for the analyzed contaminants, whereas additional biomarkers were found using MALDI. The results demonstrate the suitability of MSI for investigations on the spatial distribution of accumulated organic contaminants. However, both MSI methods required high exposure concentrations. Further improvements of ionization methods would be needed to address environmentally relevant concentrations.

Importance of Dietary Uptake for in Situ Bioaccumulation of Systemic Fungicides Using Gammarus pulex as a Model Organism

Bioaccumulation of organic contaminants from contaminated food sources might pose an underestimated risk toward shredding invertebrates. This assumption is substantiated by monitoring studies observing discrepancies of predicted tissue concentrations determined from laboratory-based experiments compared with measured concentrations of systemic pesticides in gammarids. To elucidate the role of dietary uptake in bioaccumulation, gammarids were exposed to leaf material from trees treated with a systemic fungicide mixture (azoxystrobin, cyprodinil, fluopyram, and tebuconazole), simulating leaves entering surface waters in autumn. Leaf concentrations, spatial distribution, and leaching behavior of fungicides were characterized using liquid chromatography coupled with high-resolution tandem mass spectrometry (LC-HRMS/MS) and matrix-assisted laser desorption ionization-mass spectrometric imaging. The contribution of leached fungicides and fungicides taken up from feeding was assessed by assembling caged (no access) and uncaged (access to leaves) gammarids. The fungicide dynamics in the test system were analyzed using LC-HRMS/MS and toxicokinetic modeling. In addition, a summer scenario was simulated where water was the initial source of contamination and leaves contaminated by sorption. The uptake, translocation, and biotransformation of systemic fungicides by trees were compound-dependent. Internal fungicide concentrations of gammarids with access to leaves were much higher than in caged gammarids of the autumn scenario, but the difference was minimal in the summer scenario. In food choice and dissectioning experiments gammarids did not avoid contaminated leaves and efficiently assimilated contaminants from leaves, indicating the relevance of this exposure pathway in the field. The present study demonstrates the potential impact of dietary uptake on in situ bioaccumulation for shredders in autumn, outside the main application period. The toxicokinetic parameters obtained facilitate modeling of environmental exposure scenarios. The uncovered significance of dietary uptake for detritivores warrants further consideration from scientific as well as regulatory perspectives. Environ Toxicol Chem 2023;42:1993-2006. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Simultaneously discovering the fate and biochemical effects of pharmaceuticals through untargeted metabolomics

Untargeted metabolomics is an established approach in toxicology for characterising endogenous metabolic responses to xenobiotic exposure. Detecting the xenobiotic and its biotransformation products as part of the metabolomics analysis provides an opportunity to simultaneously gain deep insights into its fate and metabolism, and to associate the internal relative dose directly with endogenous metabolic responses. This integration of untargeted exposure and response measurements into a single assay has yet to be fully demonstrated. Here we assemble a workflow to discover and analyse pharmaceutical-related measurements from routine untargeted UHPLC-MS metabolomics datasets, derived from in vivo (rat plasma and cardiac tissue, and human plasma) and in vitro (human cardiomyocytes) studies that were principally designed to investigate endogenous metabolic responses to drug exposure. Our findings clearly demonstrate how untargeted metabolomics can discover extensive biotransformation maps, temporally-changing relative systemic exposure, and direct associations of endogenous biochemical responses to the internal dose.

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.

Triazophos-induced Respiratory and Behavioral Effects and Development of Adverse Outcome Pathway (AOP) for short-term Exposed Freshwater Snail, Bellamya Bengalensis

The physiological effects of triazophos were examined using respiratory and behavioral endpoints in Bellamya bengalensis under a 96-hour acute exposure regime. Physiological manifestation of respiratory stress was measured using the rate of oxygen consumption while behavioral toxicity was measured using crawling reflexes, touch response, and mucus production. The threshold effect values for LOEC (Lowest Observed Effect Concentration), NOEC (No Observed Effect Concentration), and MATC (Maximum Acceptable Toxicant Concentration) at 96 h were 0.40, 0.60, and 0.075 mg/l, respectively. Definitive 96 h acute exposures for both respiratory and behavioral endpoints tests were determined using a control group and concentrations ranging from 0.40 to 1.60 mg/l monitored for 24, 48, 72, and 96 h. Test organisms irrespective of exposure concentration demonstrated an initial rise in oxygen consumption rate after 24 h, followed by a progressive decrease in toxicant concentration and exposure period. The in silico structural analysis presents triazophos as having an electrophilic toxic structure similar to choline esterase inhibitors, and also capable of inducing oxidative stress. The AOP highlighted neurotoxicity and oxidative stress as plausible pathways of triazophos toxicity in mollusk species.

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.

Variation in the Transcriptome Response and Detoxification Gene Diversity Drives Pesticide Tolerance in Fishes

Pesticides are critical for invasive species management but often have negative effects on nontarget native biota. Tolerance to pesticides should have an evolutionary basis, but this is poorly understood. Invasive sea lamprey (Petromyzon marinus) populations in North America have been controlled with a pesticide lethal to them at lower concentrations than native fishes. We addressed how interspecific variation in gene expression and detoxification gene diversity confer differential pesticide sensitivity in two fish species. We exposed sea lamprey and bluegill (Lepomis macrochirus), a tolerant native species, to 3-trifluoromethyl-4-nitrophenol (TFM), a pesticide commonly used in sea lamprey control. We then used whole-transcriptome sequencing of gill and liver to characterize the cellular response in both species. Comparatively, bluegill exhibited a larger number of detoxification genes expressed and a larger number of responsive transcripts overall, which likely contributes to greater tolerance to TFM. Understanding the genetic and physiological basis for pesticide tolerance is crucial for managing invasive species.

Health-related toxicity of emerging per- and polyfluoroalkyl substances: Comparison to legacy PFOS and PFOA

Per- and polyfluoroalkyl substances (PFAS) are highly persistent, manufactured chemicals used in various manufacturing processes and found in numerous commercial products. With over 9000 compounds belonging to this chemical class, there is increasing concern regarding human exposure to these compounds due to their persistent, bioaccumulative, and toxic nature. Human exposure to PFAS may occur from a variety of exposure sources, including, air, food, indoor dust, soil, water, from the transfer of PFAS from non-stick wrappers to food, use of cosmetics, and other personal care products. This critical review presents recent research on the health-related impacts of PFAS exposure, highlighting compounds other than Perfluorooctanoic acid (PFOA) and Perfluoroctane sulfonate (PFOS) that cause adverse health effects, updates the current state of knowledge on PFAS toxicity, and, where possible, elucidates cause-and-effect relationships. Recent reviews identified that exposure to PFAS was associated with adverse health impacts on female and male fertility, metabolism in pregnancy, endocrine function including pancreatic dysfunction and risk of developing Type 2 diabetes, lipid metabolism and risk of childhood adiposity, hepatic and renal function, immune function, cardiovascular health (atherosclerosis), bone health including risk for dental cavities, osteoporosis, and vitamin D deficiency, neurological function, and risk of developing breast cancer. However, while cause-and-effect relationships for many of these outcomes were not able to be clearly elucidated, it was identified that 1) the evidence derived from both animal models and humans suggested that PFAS may exert harmful impacts on both animals and humans, however extrapolating data from animal to human studies was complicated due to differences in exposure/elimination kinetics, 2) PFAS precursor kinetics and toxicity mechanism data are still limited despite ongoing exposures, and 3) studies in humans, which provide contrasting results require further investigation of the long-term-exposed population to better evaluate the biological toxicity of chronic exposure to PFAS.

Effects of glyphosate on cladocera: A synthetic review

Glyphosate [N-(phosphonomethyl) glycine] is currently the most widely used herbicide worldwide. Its application in agricultural and urban areas can lead to the dispersion and arrival to aquatic systems causing environmental deterioration with detrimental effects on the inhabiting biota. This is triggered not only by the herbicide per se but also its metabolite aminomethyl-phosphonic acid (AMPA), which can be highly toxic to many aquatic organisms. Water fleas are some of the key components in aquatic food webs, being one of the most sensitive groups to pollutants. Although being often used in standardized toxicity tests, they are comparatively less studied in relation to glyphosate exposition. Here we examine the current scientific literature regarding the acute and sublethal toxicity of glyphosate in the Cladocera taxonomic group, with special comparisons between the active ingredient (A.I) and formulations. Our results document a high variation in the lethal concentrations reported for different cladoceran species, due to the high diversity of products used in the toxicity tests. Most articles accounting for sublethal effects were performed on the standard Daphnia magna species. Reproduction, including decreased fecundity and delayed age of first reproduction, is usually one of the most severely affected individual traits. Although still scarce, studies documenting metabolic and genetic alterations might provide accurate information on the mechanisms of action of the herbicide.

Molting enhances internal concentrations of fipronil and thereby decreases survival of two estuarine resident marine crustaceans

In aquatic arthropods, molting is a pivotal physiological process for normal development, but it may also expose them to higher risks from xenobiotics, because the organism may take up additional water during that time. This study aimed to assess the effects of molting on bioconcentration and survival after 96-h exposure to insecticide fipronil with or without oxygenase (CYP450s) inhibitor piperonyl butoxide (PBO) of two estuarine resident marine crustacean species: the sand shrimp Crangon uritai and the kuruma prawn Penaeus japonicus, with 96-h LC₅₀ value of fipronil = 2.0 µg/L and 0.2 µg/L, respectively. Two graded concentrations included group high (H) (equivalent to the 96-h LC₅₀ values) and low (L) (one-tenth of the H group concentration). Molting and survival were individually checked, and internal concentrations of fipronil and its metabolites (fipronil desulfinyl, fipronil sulfide, fipronil sulfone) were measured. The results showed that, only fipronil and fipronil sulfone were detected from organism, and that internal concentrations of these insecticides in molted specimens were higher than those of unmolted ones but comparable with those of dead ones. Accordingly, mortality was more frequent in molted specimens than those that were unmolted. Furthermore, involvement of oxygenase and higher lethal body burden threshold may confer higher tolerance to fipronil in sand shrimp than in the kuruma prawn. This study is the first to demonstrate that the body-residue-based approach is useful for deciphering the causal factors underlying fipronil toxicity, but highlights the need to consider physiological factors in arthropods, which influence and lie beyond body burden, molting and drug metabolism.

Systematic Underestimation of Pesticide Burden for Invertebrates under Field Conditions: Comparing the Influence of Dietary Uptake and Aquatic Exposure Dynamics

Pesticides used in agriculture can end up in nearby streams and can have a negative impact on nontarget organisms such as aquatic invertebrates. During registration, bioaccumulation potential is often investigated using laboratory tests only. Recent studies showed that the magnitude of bioaccumulation in the field substantially differs from laboratory conditions. To investigate this discrepancy, we conducted a field bioaccumulation study in a stream known to receive pollutant loadings from agriculture. Our work incorporates measurements of stream pesticide concentrations at high temporal resolution (every 20 min), as well as sediment, leaves, and caged gammarid analyses (every 2-24 h) over several weeks. Of 49 investigated pesticides, 14 were detected in gammarids with highly variable concentrations of up to 140 ± 28 ng/gww. Toxicokinetic modeling using laboratory-derived uptake and depuration rate constants for azoxystrobin, cyprodinil, and fluopyram showed that despite the highly resolved water concentrations measured, the pesticide burden on gammarids remains underestimated by a factor of 1.9 ± 0.1 to 31 ± 3.0, with the highest underestimations occurring after rain events. Including dietary uptake from polluted detritus leaves and sediment in the model explained this underestimation only to a minor proportion. However, suspended solids analyzed during rain events had high pesticide concentrations, and uptake from them could partially explain the underestimation after rain events. Additional comparison between the measured and modeled data showed that the pesticide depuration in gammarids is slower in the field. This observation suggests that several unknown mechanisms may play a role, including lowered enzyme expression and mixture effects. Thus, it is important to conduct such retrospective risk assessments based on field investigations and adapt the registration accordingly.

Determination of In Vivo Biotransformation Kinetics Using Early-Time Biota Concentrations

Technical challenges have hampered the characterization of biotransformation kinetics-a critical link in understanding and predicting the toxicokinetics and ecotoxicology of organic compounds. A shortcut approach to characterize the in vivo biotransformation rate constant (k_M ) with incomplete pathway or metabolite details was proposed. The value of k_M can be derived as 2tln1fPC(t)) , with f_PC (t) being the molar equivalent fraction of the parent compound (PC) at an early time t in both constant exposure and decay source chemical uptake scenarios. The approximation-based k_M values agreed well with k_M values derived from rigorous fitting or toxicokinetic modeling (n = 42, root mean square error = 0.30) with accuracy exceeding those of typical toxicokinetic or partitioning models. The method is accurate when sampling time is adequately resolved (i.e., t < ln(2)/k_M ) but will likely produce biased k_M values with improper time-averaging. The approximate equation yields consistent theoretical expectations for fast and slow biotransformation reactions and is fully compatible with standard bioaccumulation and toxicity testing protocols. The simplification strategy circumvents statistical complications and numerical issues inherent in regressing or modeling the toxicokinetics of multimetabolite systems and may be adapted to similar problems at other physiological scales or ecotoxicological contexts. The method can help advance interspecies comparison of chemical metabolism and support the development of in vitro-in vivo extrapolations and in silico models needed for building next-generation ecological and health risk-assessment practices. Environ Toxicol Chem 2022;41:148-158. © 2021 SETAC.

Biochemical, physiological (haematological, oxygen-consumption rate) and behavioural effects of mercury exposures on the freshwater snail, Bellamya bengalensis

The widespread occurrence of Mercury (Hg) and its derivatives in the aquatic environment and risks to the health of local populations has necessitated investigations into its toxic effects on sessile species. The toxicity of Mercury was observed sequentially from 96 h acute exposure regime (behavioural endpoints) to chronic durations (haematological and biochemical toxicity endpoints) in Bellamya bengalensis. Time-dependent lethal endpoints for acute toxicity (LC50) of mercury i.e., 24,48,72 and 96 h were estimated as 0.94, 0.88, 0.69 and 0.40 mg/l respectively. Threshold effect values i.e., LOEC (Lowest Observed Effect Concentration), NOEC (No Observed Effect Concentration) and MATC (Maximum Acceptable Toxicant Concentration) at 96 h were found to be 0.10, 0.05, 0.039 mg/l respectively. The study of oxygen consumption rate and behavioural changes during acute toxicity and haematological and biochemical responses during chronic toxicity to sublethal concentrations (10% and 20% of 96 h LC50) of mercury to the snail were also conducted. The organisms showed initial elevation at 24 h but later gradual decrease in oxygen consumption rate with the increase of concentration of mercury and time of exposure. For behavioural studies, variable test concentrations from 0.00 to 1.00 mg/l were used for 24, 48, 72 and 96 h. The crawling activity and clumping tendency decreased with the progress of time at all treatment periods and stopped ultimately at 96 h of exposure from 0.7 mg/l onwards whereas touch reflex was not observed at 96 h exposure at all treatments except at 0.09 mg/l. In haemocyte count, no significant variation was observed among control values between various exposure periods (p > 0.05) though variations were observed in sub-lethal concentrations versus control at all treatment duration (7, 14, 21, 28d, p < 0.05). In biochemical response study, the protein content in hepatopancreas of the snails treated at sublethal concentrations of mercury (10% and 20% of 96 h LC50) reduced significantly versus control after 21d of exposure (p < 0.05). In gonads, the protein content of the treated snails significantly reduced at all treatment concentrations versus control at all exposure times (p < 0.05). Based on the safe levels indicated above, the concentration of 0.01 to 0.04 ppm of mercury can be considered safe for Bellamya bengalensis and any less-hardy aquatic species. These responses elicited by our molluscan model will not only help in biomonitoring of environmental mercury contamination in water bodies but will also provide support to ecological health and risk assessment.

Differences in marine megafauna in vitro sensitivity highlights the need for species-specific chemical risk assessments

Sea turtles, dolphins and dugongs can be exposed to large mixtures of contaminants due to the proximity of foraging locations to anthropogenic inputs. Differences in accumulation and effect result in differences of chemical risk to these species. However, little is known about the effect of contaminants in marine wildlife. Cell-based, or in vitro, exposure experiments offer an ethical alternative to investigate the effect of contaminants in wildlife. Data from in vitro studies can then be placed in an environmental context, by using screening risk assessments, comparing effect data with accumulation data from the literature, to identify risk to populations of marine wildlife. Cytotoxicity of Cr⁶⁺, Cd²⁺, Hg²⁺, 4,4'-DDE, and PFNA were investigated in primary skin fibroblasts of green turtles, loggerhead turtles, hawksbill turtles, dugongs, Burrunan dolphins, and common bottlenose dolphins. The general order of toxicity for all species was Hg²⁺> Cr⁶⁺ > Cd²⁺> 4,4'-DDE > PFNA, and significant differences in cytotoxicity were found between species for Cr⁶⁺, Cd²⁺ and PFNA. For Cd²⁺, in particular, cells from turtle species were less sensitive than mammalian species, and dugong cells were by far the most sensitive. The results from the cytotoxicity assay were then used in combination with published data on tissue contaminant concentrations to calculate risk quotients for identifying populations of each species most at risk from these chemicals. Cr, Cd and Hg were identified as posing risk in all six species. Dugongs were particularly at risk from Cd accumulation and dolphin species were particularly at risk from Hg accumulation. These results demonstrate the importance of using species-specific effect and accumulation data for developing chemical risk assessments and can be used to inform managers of priority contaminants, species, or populations. Development of additional in vitro endpoints, and improving links between in vitro and in vivo effects, would further improve this approach to understanding chemical risk in marine megafauna.

An Extensive Metabolomics Workflow to Discover Cardiotoxin-Induced Molecular Perturbations in Microtissues

Discovering modes of action and predictive biomarkers of drug-induced structural cardiotoxicity offers the potential to improve cardiac safety assessment of lead compounds and enhance preclinical to clinical translation during drug development. Cardiac microtissues are a promising, physiologically relevant, in vitro model, each composed of ca. 500 cells. While untargeted metabolomics is capable of generating hypotheses on toxicological modes of action and discovering metabolic biomarkers, applying this technology to low-biomass microtissues in suspension is experimentally challenging. Thus, we first evaluated a filtration-based approach for harvesting microtissues and assessed the sensitivity and reproducibility of nanoelectrospray direct infusion mass spectrometry (nESI-DIMS) measurements of intracellular extracts, revealing samples consisting of 28 pooled microtissues, harvested by filtration, are suitable for profiling the intracellular metabolome and lipidome. Subsequently, an extensive workflow combining nESI-DIMS untargeted metabolomics and lipidomics of intracellular extracts with ultra-high performance liquid chromatography-mass spectrometry (UHPLC-MS/MS) analysis of spent culture medium, to profile the metabolic footprint and quantify drug exposure concentrations, was implemented. Using the synthetic drug and model cardiotoxin sunitinib, time-resolved metabolic and lipid perturbations in cardiac microtissues were investigated, providing valuable data for generating hypotheses on toxicological modes of action and identifying putative biomarkers such as disruption of purine metabolism and perturbation of polyunsaturated fatty acid levels.

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.

High Tolerance and Delayed Responses of Daphnia magna to Neonicotinoid Insecticide Imidacloprid: Toxicokinetic and Toxicodynamic Modeling

Species sensitivity to neonicotinoids has been shown to be highly variable among aquatic invertebrates. Toxicokinetic and toxicodynamic (TKTD) models were constructed to mechanistically elucidate the susceptibility of Daphnia magna to imidacloprid. D. magna was highly tolerant to single short-term exposure to imidacloprid (96-h LC50 of 8.47 μg/mL), but delayed and carry-over toxicity occurred under repeated pulse exposures. Kinetic distribution of imidacloprid between exoskeleton and soft tissues of D. magna was evaluated using a newly developed method. Approximately 84% imidacloprid was distributed to soft tissues but was rapidly depurated from the tissue (t_1/2 of 1.2 h), resulting in low bioaccumulation and high tolerance. TKTD modeling also successfully simulated the survival of D. magna after pulsed exposures. The calculated recovery time was 45 d, indicating significant delayed and carry-over toxicity of the insecticide. While complete elimination of imidacloprid only took about 5 h (TK), slow damage recovery (45 d) caused slow organism recovery (TD). Consequently, although D. magna was tolerant to imidacloprid due to fast depuration from soft tissue, long damage recovery time significantly enhanced the toxicity under repeated pulse exposures. Our study highlights the necessity of integrating delayed and carry-over toxicity quantification in assessing the risk of neonicotinoids to aquatic invertebrates.

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

Pyrethroid insecticides are known to be highly toxic to most aquatic nontarget organisms, but little is known about the mechanisms causing some species to be highly sensitive while others are hardly affected by the pyrethroids. The aim of the present study was to measure the sensitivity (EC₅₀-values) of 10 aquatic invertebrates toward a 24 h pulse of the pyrethroid cypermethrin and subsequently test if the difference in sensitivity could be explained by measured morphological and physiological traits and modeled toxicokinetic (TK) and toxicodynamic (TD) parameters. Large differences were observed for the measured uptake and elimination kinetics, with bioconcentration factors (BCFs) ranging from 53 to 2337 at the end of the exposure. Similarly, large differences were observed for the TDs, and EC₅₀-values after 168 h varied 120-fold. Modeling the whole organism cypermethrin concentrations indicated compartmentation into a sorbed fraction and two internal fractions: a bioavailable and non-bioavailable internal fraction. Strong correlations between surface/volume area and the TK parameters (sorption and uptake rate constants and the resulting BCF) were found, but none of the TK parameters correlated with sensitivity. The only parameter consistently correlating with sensitivity across all species was the killing rate constant of the GUTS-RED-SD model (the reduced general unified threshold models of survival assuming stochastic death), indicating that sensitivity toward cypermethrin is more related to the TD parameters than to TK parameters.

Single and combined effects of insecticides on multi-level biomarkers in the non-target amphipod Gammarus fossarum exposed to environmentally realistic levels

Aquatic media are ultimate recipients of various contaminants including pesticides pervasively applied in agrosystems. Characterizing the ecotoxicity of pesticides and their mixtures to aquatic wildlife at field-realistic levels is thus crucial for environmental risk assessment. This study aims at assessing the effects of two current-use insecticides, imidacloprid and chlorpyrifos, on Gammarus fossarum using multi-level biomarkers. In microcosms, gammarids were exposed for 72 h to insecticides tested individually or in mixture at 0.01, 0.1 and 1 μg/L of each chemical. Multi-metric responses were assessed at the individual level (behavioural traits: locomotion, respiration and amplexus formation) and the cellular level (enzymes involved in growth, moulting, digestion and cell stress). The results showed insecticide-elicited behavioural and biochemical responses from the lowest concentration of 0.01 μg/L. Overall, single exposures stimulated behavioural traits and inhibited enzymatic activities, highlighting subtle impacts at different organizational levels but these were not dose related. For binary mixtures, antagonistic effects (i.e. less-than-additive) on biomarkers were mainly observed when compared with single exposures. Multi-variable analyses indicated the complementarity of behavioural and biochemical biomarkers in identifying sublethal biological alterations and dose-dependent multiple action sites of insecticides. Besides, the mortality observed only for the mixture at 1 μg/L demonstrated a high lethal potential of insecticides in a simple binary combination. To conclude, this study demonstrates disturbances in individual performances and cellular impairments occurring at environmentally realistic exposure levels in a non-target wild species. Since the sublethal effects, such as those identified with this multi-biomarker approach, could lead to long-term alterations in population dynamics of agricultural areas, they constitute promising early endpoints for risk assessment of insecticides.

Metabolomic profiles associated with exposure to per- and polyfluoroalkyl substances (PFASs) in aquatic environments

Per- and polyfluoroalkyl substances (PFASs) are frequently detected in aquatic environments. Longer chained perfluoroalkyl acids (PFAAs), in particular, have been found to bioaccumulate in a broad range of aquatic biota. PFAAs have a physiochemical similarity to naturally occurring fatty acids and could potentially disrupt metabolic processes, however, there has been limited study in this area, especially in aquatic species. In this study, the associations between PFAAs and metabolite profiles were investigated in crustaceans. Eastern School Prawn (Metapenaeus macleayi) were obtained from three different locations (n = 15 per location) with similar environmental conditions but different levels of PFAA contamination. The concentrations of PFAAs, fatty acids and amino acids were analysed and differences in PFAA and metabolite profiles were evaluated. Different PFAA profiles were mirrored by significant differences in the composition of both fatty acid and amino acid profiles, indicating a potential association between PFAA concentration and the composition of metabolites in prawns. These results highlight a need for further research investigating the impacts of PFAA exposure, with the current study providing a foundation for further investigation of the relationship between PFAA bioaccumulation and organism metabolism.

Diurnal variation around an optimum and near-critically high temperature does not alter the performance of an ectothermic aquatic grazer

The growing threat of global climate change has led to a profusion of studies examining the effects of warming on biota. Despite the potential importance of natural variability such as diurnal temperature fluctuations, most experimental studies on warming are conducted under stable temperatures. Here, we investigated whether the responses of an aquatic invertebrate grazer (Lymnaea stagnalis) to an increased average temperature differ when the thermal regime is either constant or fluctuates diurnally. Using thermal response curves for several life-history and immune defense traits, we first identified the optimum and near-critically high temperatures that Lymnaea potentially experience during summer heat waves. We then exposed individuals that originated from three different populations to these two temperatures under constant or fluctuating thermal conditions. After 7 days, we assessed growth, reproduction, and two immune parameters (phenoloxidase-like activity and antibacterial activity of hemolymph) from each individual. Exposure to the near-critically high temperature led to increased growth rates and decreased antibacterial activity of hemolymph compared to the optimum temperature, whilst temperature fluctuations had no effect on these traits. The results indicate that the temperature level per se, rather than the variability in temperature was the main driver altering trait responses in our study species. Forecasting responses in temperature-related responses remains challenging, due to system-specific properties that can include intraspecific variation. However, our study indicates that experiments examining the effects of warming using constant temperatures can give similar predictions as studies with fluctuating thermal dynamics, and may thus be useful indicators of responses in nature.

Occurrence of neonicotinoids and fipronil in estuaries and their potential risks to aquatic invertebrates

This study aimed to evaluate and qualify field-based potential risks of seven neonicotinoid and phenylpyrazole (fipronil) insecticides on aquatic invertebrates, including estuary-resident marine crustaceans. One hundred and ninety-three estuarine water samples, with salinity ranging from 0.5 to 32.7, were collected from four estuarine sites in the Seto Inland Sea of Japan, in 2015-2018 and the insecticide levels were measured. Five neonicotinoid and fipronil insecticides were successfully identified, and their occurrence varied temporally. Marine crustaceans were simultaneously harvested every month from one of the estuarine water sampling sites in 2015-2017. Three predominant crustacean species, kuruma prawn (Penaeus japonicus), sand shrimp (Crangon uritai), and mysid (Neomysis awatschensis), were captured and their seasonal presence was species independent. A 96-h laboratory toxicity study with the insecticides using kuruma prawn, sand shrimp, and a surrogate mysid species (Americamysis bahia) indicated that fipronil exerted the highest toxicity to the three crustaceans. Using both toxicity data and insecticide occurrence in estuarine water (salinity ≥10, n = 169), the potential risks on the three marine crustaceans were quantified by calculating the proportion of mixture toxicity effects (P_mix). The P_mix of seven neonicotinoids on the crustaceans was less than 0.8%, which is likely to be too low to indicate adverse effects caused by the insecticides. However, short temporal detection of fipronil (exclusively in June and July) significantly affected the P_mix, which presented the maximal P_mix values of 21%, 3.4%, and 72% for kuruma prawn, sand shrimp, and mysid, respectively, indicating a significant effect on the organisms. As for estuarine water (salinity <10), some water samples contained imidacloprid and fipronil exceeding the freshwater benchmarks for aquatic invertebrates. The present study provides novel insights into the seasonally varying risks of insecticides to estuarine crustaceans and highlights the importance of considering whether ecological risk periods coincide with crustacean presence.

Exploring inter-species sensitivity to a model hydrocarbon, 2-Methylnaphtalene, using a process-based model

We compared inter-species sensitivity to a model narcotic compound, 2-Methylnaphthalene, to test if taxonomical relatedness, feeding guilds, and trophic level govern species sensitivities on species distributed in different regions. We fitted a toxicokinetic-toxicodynamic model to survival patterns over time for 26 species using new and raw data from the literature. Species sensitivity distributions provided little insight into understanding patterns in inter-species sensitivity. The range of no-effect concentrations (NEC) obtained for 26 species showed little variation (mean 0.0081 mM; SD 0.009). Results suggest that the NEC alone does not explain the complexity of the species tolerances. The dominant rate constant and the derived time to observe an effect (t₀), a function of concentration, might provide the means for depicting patterns in sensitivity and better ecotoxicological testing. When comparing the t₀ functions, we observed that Arctic species have shorter time frames to start showing effects. Mollusks and second trophic level species took longer to build up a lethal body burden than the rest. Coupling our results with fate and transport models would allow forecasting narcotic compounds toxicity in time and thus improve risk assessment.

Neonicotinoid insecticides are potential substrates of the multixenobiotic resistance (MXR) mechanism in the non-target invertebrate, Dreissena sp

Mussels are among the most frequently used invertebrate animals in aquatic toxicology to detect toxic exposure in the environment. The presence and activity of a cellular defence system, the multixenobiotic resistance (MXR) mechanism, was also established in these organisms. In isolated gill tissues of dreissenid mussels (D. bugensis) the MXR activity was assayed after treatment by commercially available insecticides (formulated products) which contain neonicotinoids as their active ingredients: Actara (thiamethoxam), Apacs (clothianidin), Calypso (thiacloprid) and Kohinor (imidacloprid), respectively. While applying the accumulation assay method, 0.5 μM rhodamine B was used as model substrate and 20 μM verapamil as model inhibitor of the MXR mechanism. In acute (in vitro) experiments when isolated gills were co-incubated in graded concentrations of insecticides and rhodamine B simultaneously, Calypso and Kohinor treatment resulted increasing rhodamine accumulation. Chemical analysis of gills in vitro incubated in insecticides demonstrated higher tissue concentrations of thiamethoxam, clothianidin and thiacloprid in the presence of verapamil suggesting that the active ingredients of Actara, Apacs and Calypso are potential substrates of the MXR mediated cellular efflux. In contrast, verapamil did significantly alter the accumulated imidacloprid concentrations in gills, suggesting that the active component of Kohinor is not transported by the MXR mechanism. Chronic (in vivo) exposures of the intact animals in lower, 1, 10 mg/L concentration of neonicotinoid products, resulted in a decreased level of both rhodamine accumulation and verapamil inhibition by the 12th-14th days of treatment. These results suggest an enhancement of MXR activity (chemostimulation), building up gradually in the animals exposed to Actara, Apacs and Kohinor, respectively. Neonicotinoid-type insecticides are generally considered as selective neurotoxins for insects, targeting the nicotinic type acetylcholine receptors (nAChRs) in their central nervous system. Our present results provide the first evidences that neonicotinoid insecticides are also able to alter the transmembrane transport mechanisms related to the MXR system.

Effects of florfenicol and oxytetracycline on the tropical cladoceran Ceriodaphnia silvestrii: A mixture toxicity approach to predict the potential risks of antimicrobials for zooplankton

Antimicrobials are commonly used in aquaculture to treat infectious diseases in fish. The overuse of these chemicals, however, has made them a contamination source for the aquatic environments. In this study, single and combined effects of florfenicol (FLO) and oxytetracycline (OTC), two antimicrobials widely used in the fish farming, were evaluated in acute and chronic toxicity tests using the tropical cladoceran Ceriodaphnia silvestrii as a model species. Also, a preliminary risk characterization of FLO and OTC for zooplankton was carried out, taking into account different exposure scenarios. The results obtained revealed that FLO and OTC have adverse effects on the mobility, reproduction and population growth rate of C. silvestrii in single exposures. In addition, mixture effects on the C. silvestrii were more severe than predicted effects based on the Concentration Addition model, showing a synergistic deviation for the mobility and a dose-level dependent deviation for the reproduction (synergism at higher levels than EC₆₀). In relation to the risk characterization, risk quotients (RQs) exceeded 1 for chronic toxicity data obtained in both OTC and mixture exposures, indicating that the concentrations of these chemicals in Brazilian freshwater bodies could potentially present risks for the reproduction of zooplankton species in tropical regions. The RQs obtained for the mixtures were higher than those obtained for each chemical separately. Therefore, it is highly recommended that RQs are derived from single and mixture exposure data in order to obtain a more accurate risk characterization.

Internal Concentrations in Gammarids Reveal Increased Risk of Organic Micropollutants in Wastewater-Impacted Streams

Internal concentrations link external exposure to the potential effect, as they reflect what the organisms actually take up and experience physiologically. In this study, we investigated whether frequently detected risk-driving substances in water were found in the exposed organisms and if they are classified the same based on the whole body internal concentrations. Field gammarids were collected upstream and downstream of ten wastewater treatment plants in mixed land use catchments. The sampling was conducted in autumn and winter, during low flow conditions when diffuse agricultural input was reduced. The field study was complemented with laboratory and flume experiments to determine the bioaccumulation potentials of selected substances. For 32 substances, apparent bioaccumulation factors in gammarids were determined for the first time. With a sensitive multiresidue method based on online-solid phase extraction followed by liquid chromatography coupled to high resolution mass spectrometry, we detected 63 (semi-) polar organic substances in the field gammarids, showing higher concentrations downstream than upstream. Interestingly, neonicotinoids, which are particularly toxic toward invertebrates, were frequently detected and were further determined as major contributors to the toxic pressure based on the toxic unit approach integrating internal concentration and toxic potency. The total toxic pressure based on internal concentrations was substantially higher compared to when external concentrations were used. Thus, internal concentrations may add more value to the current environmental risk assessment that is typically based solely on external exposure.

Comparison of the behavioural effects of pharmaceuticals and pesticides on Diamesa zernyi larvae (Chironomidae)

Several studies have indicated the presence of contaminants in Alpine aquatic ecosystems. Even if measured concentrations are far below those that cause acute effects, continuous exposure to sub-lethal concentrations may have detrimental effects on the aquatic species present in these remote environments. This may lead to a cascade of indirect effects at higher levels of the ecological hierarchy (i.e., the community). To improve the determination of ecologically relevant risk endpoints, behavioural alterations in organisms due to pollutants are increasingly studied in ecotoxicology. In fact, behaviour links physiological function with ecological processes, and can be very sensitive to environmental stimuli and chemical exposure. This is the first study on behavioural alteration in a wild population of an Alpine species. In the present study, a video tracking system was standardized and subsequently used to identify contaminant-induced behavioural alterations in Diamesa zernyi larvae (Diptera, Chironomidae). Diamesa zernyi larvae, collected in an Italian Alpine stream (Rio Presena, Trentino Region), were acclimated for 24 h and successively exposed to several aquatic contaminants (pesticides: chlorpyrifos, metolachlor, boscalid, captan; pharmaceuticals: ibuprofen, furosemide, trimethoprim) at concentrations corresponding to their Lowest Observed Effect Concentration (LOEC). After 24, 48, 72, and 96 h of exposure, changes in the distance moved, the average speed, and the frequency of body bends were taken to reflect contaminant- and time-dependent effects on larval behaviour. In general, metolachlor, captan, and trimethoprim tended to reduce all the endpoints under consideration, whereas chlorpyrifos, boscalid, ibuprofen, and furosemide seemed to increase the distances moved by the larvae. This could be related to the different mechanisms of action of the investigated chemicals. Independently of the contaminant, after 72 h a general slowing down of all the behavioural activities occurred. Finally, we propose a behavioural stress indicator to compare the overall behavioural effects induced by the various contaminants.

Freshwater shrimps as sensitive test species for the risk assessment of pesticides in the tropics

The aquatic risk assessment of pesticides in tropical areas has often been disputed to rely on toxicity data generated from tests performed with temperate species. Given the differences in ecosystem structure between temperate and tropical ecosystems, test species other than those used in temperate regions have been proposed as surrogates for tropical aquatic effect assessments. Freshwater shrimps, for example are important components of tropical freshwater ecosystems, both in terms of their role in ecosystem functioning and their economic value. In the present study, available toxicity data of (tropical and sub-tropical) freshwater shrimps for insecticides and fungicides were compiled and compared with those available for Daphnia magna and other aquatic invertebrates. Freshwater shrimps appeared to be especially sensitive to GABA-gated chloride channel antagonist and sodium channel modulator insecticides. However, shrimp taxa showed a moderate and low sensitivity to acetylcholinesterase inhibiting insecticides and fungicides respectively. Implications for the use of freshwater shrimps in tropical pesticide effect assessments and research needs are discussed.

Using extirpation to evaluate ionic tolerance of freshwater fish

Field data of fish occurrences and specific conductivity were used to estimate the tolerance of freshwater fish to elevated ion concentrations and to compare the differences between species- and genus-level analyses for individual effects. We derived extirpation concentrations at the 95th percentile (XC95) of a weighted cumulative frequency distribution for fish species inhabiting streams of the central and southern Appalachians by customizing methods used previously with macroinvertebrate genera. Weighting factors were calculated based on the number of sites in basins where each species occurred, reducing overweighting observations of species restricted to fewer basins. Comparing the species- and genus-level fish XC95 values, XC95s for fish genera were near the XC95s for the most salt-tolerant species in the genus. Therefore, a genus-level effect threshold is not reliably predictive of species-level extirpation, unless the genus is monospecific in the assessed assemblage. Of the 101 fish species XC95 values, 5% were <509 and 10% were <565 µS/cm. The lowest XC95 for a species was 322 µS/cm, which is >300 µS/cm, the exposure estimated to extirpate 5% of macroinvertebrate genera in the central Appalachians. Above 509 µS/cm, 41 of the 101 species are expected to decline in occurrence. Environ Toxicol Chem 2018;37:871-883. Published 2017 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

Charting a path towards non-destructive biomarkers in threatened wildlife: A systematic quantitative literature review

Threatened species are susceptible to irreversible population decline caused by adverse sub-lethal effects of chemical contaminant exposure. It is therefore vital to develop the necessary tools to predict and detect these effects as early as possible. Biomarkers of contaminant exposure and effect are widely applied to this end, and a significant amount of research has focused on development and validation of sensitive and diagnostic biomarkers. However, progress in the use biomarkers that can be measured using non-destructive techniques has been relatively slow and there are still many difficulties to overcome in the development of sound methods. This paper systematically quantifies and reviews studies that have aimed to develop or validate non-destructive biomarkers in wildlife, and provides an analysis of the successes of these methods based on the invasiveness of the methods, the potential for universal application, cost, and the potential for new biomarker discovery. These data are then used to infer what methods and approaches appear the most effective for successful development of non-destructive biomarkers of contaminant exposure in wildlife. This review highlights that research on non-destructive biomarkers in wildlife is severely lacking, and suggests further exploration of in vitro methods in future studies.

Defence mechanisms: the role of physiology in current and future environmental protection paradigms

Ecological risk assessments principally rely on simplified metrics of organismal sensitivity that do not consider mechanism or biological traits. As such, they are unable to adequately extrapolate from standard laboratory tests to real-world settings, and largely fail to account for the diversity of organisms and environmental variables that occur in natural environments. However, an understanding of how stressors influence organism health can compensate for these limitations. Mechanistic knowledge can be used to account for species differences in basal biological function and variability in environmental factors, including spatial and temporal changes in the chemical, physical and biological milieu. Consequently, physiological understanding of biological function, and how this is altered by stressor exposure, can facilitate proactive, predictive risk assessment. In this perspective article, existing frameworks that utilize physiological knowledge (e.g. biotic ligand models, adverse outcomes pathways and mechanistic effect models), are outlined, and specific examples of how mechanistic understanding has been used to predict risk are highlighted. Future research approaches and data needs for extending the incorporation of physiological information into ecological risk assessments are discussed. Although the review focuses on chemical toxicants in aquatic systems, physical and biological stressors and terrestrial environments are also briefly considered.

The relative sensitivity of freshwater species to antimony(III): Implications for water quality guidelines and ecological risk assessments

Antimony (Sb) is a pollutant in many jurisdictions, yet its threat to aquatic biota is unclear. Water quality guidelines (WQGs) for Sb are not well established and large uncertainty factors are commonly applied in derivation. We constructed freshwater species sensitivity distributions (SSDs) for Sb(III) using available acute toxicity data sourced from temperate and tropical regional studies. A tiered ecological risk assessment (ERA) approach using risk quotients (RQs) was applied for characterisation of risks presented by Sb(III) concentrations measured in the freshwater environment. Multiple parametric models were fitted for each SSD, with the optimal model used to derive the 5% hazardous concentration (HC5), defined as protective of 95% of species, and the corresponding predicted no effect concentration (PNEC). The HC5 values for whole and temperate SSDs were estimated at 781 and 976 μg L^-1 Sb(III), respectively, while the PNECs for both datasets were 156 and 195 μg L^-1 Sb(III), respectively. Due to limited tropical data, a temperate-to-tropic extrapolation factor of 10 was used to estimate an interim PNEC for tropical regions of 20 μg L^-1 Sb(III). Based on published freshwater Sb(III) concentration values across a range of locations, potential ecological risks posed by Sb(III) in some freshwater systems studied would be classified as medium to high risk, but the majority of locations sampled would fall into the low ecological risk category. Our results facilitate the understanding of toxic effects of Sb(III) to freshwater species but also demonstrate that data for Sb ERA are extremely limited.

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.

Does Waterborne Exposure Explain Effects Caused by Neonicotinoid-Contaminated Plant Material in Aquatic Systems?

Neonicotinoids are increasingly applied on trees as protection measure against insect pests. Consequently, neonicotinoids are inevitably transferred into aquatic environments either via spray drift or surface runoff or (due to neonicotinoids' systemic nature) via senescent leaves. There particularly leaf-shredding invertebrates may be exposed to neonicotinoids through both the water phase and the consumption of contaminated leaves. In 7 day bioassays (n = 30), we examined ecotoxicological differences between these two exposure scenarios for an amphipod and an insect nymph with their feeding rate as the response variable. Organisms either experienced waterborne neonicotinoid (i.e., imidacloprid, thiacloprid, and acetamiprid) exposure only or a combined exposure (waterborne and dietary) through both the consumption of contaminated leaves and neonicotinoids leaching from leaves into water. The amphipod (7 day EC₅₀s from 0.3 to 8.4 μg/L) was more sensitive than the insect nymph (7 day EC₅₀s from 7.0 to 19.4 μg/L). Moreover, for both species, concentration-response models derived from water concentrations indicated higher effects under the combined exposure. Together with the observed inability of shredders to avoid neonicotinoid-contaminated leaves, our results emphasize the relevance of dietary exposure (e.g., via leaves) for systemic insecticides. Thus, it would be prudent to consider dietary exposure during the registration of systemic insecticides to safeguard ecosystem integrity.

Toxic Mixtures in Time-The Sequence Makes the Poison

"The dose makes the poison". This principle assumes that once a chemical is cleared out of the organism (toxicokinetic recovery), it no longer has any effect. However, it overlooks the other process of re-establishing homeostasis, toxicodynamic recovery, which can be fast or slow depending on the chemical. Therefore, when organisms are exposed to two toxicants in sequence, the toxicity can differ if their order is reversed. We test this hypothesis with the freshwater crustacean Gammarus pulex and four toxicants that act on different targets (diazinon, propiconazole, 4,6-dinitro-o-cresol, 4-nitrobenzyl chloride). We found clearly different toxicity when the exposure order of two toxicants was reversed, while maintaining the same dose. Slow toxicodynamic recovery caused carry-over toxicity in subsequent exposures, thereby resulting in a sequence effect-but only when toxicodynamic recovery was slow relative to the interval between exposures. This suggests that carry-over toxicity is a useful proxy for organism fitness and that risk assessment methods should be revised as they currently could underestimate risk. We provide the first evidence that carry-over toxicity occurs among chemicals acting on different targets and when exposure is several days apart. It is therefore not only the dose that makes the poison but also the exposure sequence.

Toward refined environmental scenarios for ecological risk assessment of down-the-drain chemicals in freshwater environments

Current regulatory practice for chemical risk assessment suffers from the lack of realism in conventional frameworks. Despite significant advances in exposure and ecological effect modeling, the implementation of novel approaches as high-tier options for prospective regulatory risk assessment remains limited, particularly among general chemicals such as down-the-drain ingredients. While reviewing the current state of the art in environmental exposure and ecological effect modeling, we propose a scenario-based framework that enables a better integration of exposure and effect assessments in a tiered approach. Global- to catchment-scale spatially explicit exposure models can be used to identify areas of higher exposure and to generate ecologically relevant exposure information for input into effect models. Numerous examples of mechanistic ecological effect models demonstrate that it is technically feasible to extrapolate from individual-level effects to effects at higher levels of biological organization and from laboratory to environmental conditions. However, the data required to parameterize effect models that can embrace the complexity of ecosystems are large and require a targeted approach. Experimental efforts should, therefore, focus on vulnerable species and/or traits and ecological conditions of relevance. We outline key research needs to address the challenges that currently hinder the practical application of advanced model-based approaches to risk assessment of down-the-drain chemicals. Integr Environ Assess Manag 2017;13:233-248. © 2016 SETAC.

Embryonic exposure to model naphthenic acids delays growth and hatching in the pond snail Lymnaea stagnalis

Naphthenic acids (NAs), a class of structurally diverse carboxylic acids with often complex ring structures and large aliphatic tail groups, are important by-products of many petrochemical processes including the oil sands mining activity of Northern Alberta. While it is evident that NAs have both acute and chronic harmful effects on many organisms, many aspects of their toxicity remain to be clarified. Particularly, while substantive data sets have been collected on NA toxicity in aquatic prokaryote and vertebrate model systems, to date, nothing is known about the toxic effects of these compounds on the embryonic development of aquatic invertebrate taxa, including freshwater mollusks. This study examines under laboratory conditions the toxicity of NAs extracted from oil sands process water (OSPW) and the low-molecular weight model NAs cyclohexylsuccinic acid (CHSA), cyclohexanebutyric acid (CHBA), and 4-tert-butylcyclohexane carboxylic acid (4-TBCA) on embryonic development of the snail Lymnaea stagnalis, a common freshwater gastropod with a broad Palearctic distribution. Evidence is provided for concentration-dependent teratogenic effects of both OSPW-derived and model NAs with remarkably similar nominal threshold concentrations between 15 and 20 mg/L and 28d EC₅₀ of 31 mg/L. In addition, the data provide evidence for substantial toxicokinetic differences between CHSA, CHBA and 4-TBCA. Together, our study introduces Lymnaea stagnalis embryonic development as an effective model to assay NA-toxicity and identifies molecular architecture as a potentially important toxicokinetic parameter in the toxicity of low-molecular weight NA in embryonic development of aquatic gastropods.

Sensitivity of the early-life stages of freshwater mollusks to neonicotinoid and butenolide insecticides

Neonicotinoid insecticides can be transported from agricultural fields, where they are used as foliar sprays or seed treatments, to surface waters by surface or sub-surface runoff. Few studies have investigated the toxicity of neonicotinoid or the related butenolide insecticides to freshwater mollusk species. The current study examined the effect of neonicotinoid and butenolide exposures to the early-life stages of the ramshorn snail, Planorbella pilsbryi, and the wavy-rayed lampmussel, Lampsilis fasciola. Juvenile P. pilsbryi were exposed to imidacloprid, clothianidin, or thiamethoxam for 7 or 28 d and mortality, growth, and biomass production were measured. The viability of larval (glochidia) L. fasciola was monitored during a 48 h exposure to six neonicotinoids (imidacloprid, thiamethoxam, clothianidin, acetamiprid, thiacloprid, or dinotefuran), or a butenolide (flupyradifurone). The 7-d LC50s of P. pilsbryi for imidacloprid, clothianidin, and thiamethoxam were ≥4000 μg/L and the 28-d LC50s were ≥182 μg/L. Growth and biomass production were considerably more sensitive endpoints than mortality with EC50s ranging from 33.2 to 122.0 μg/L. The 48-h LC50s for the viability of glochidia were ≥456 μg/L for all seven insecticides tested. Our data indicate that neonicotinoid and butenolide insecticides pose less of a hazard with respect to mortality of the two species of mollusk compared to the potential hazard to other non-target aquatic insects.

Acute Toxicity Prediction to Threatened and Endangered Species Using Interspecies Correlation Estimation (ICE) Models

Evaluating contaminant sensitivity of threatened and endangered (listed) species and protectiveness of chemical regulations often depends on toxicity data for commonly tested surrogate species. The U.S. EPA's Internet application Web-ICE is a suite of Interspecies Correlation Estimation (ICE) models that can extrapolate species sensitivity to listed taxa using least-squares regressions of the sensitivity of a surrogate species and a predicted taxon (species, genus, or family). Web-ICE was expanded with new models that can predict toxicity to over 250 listed species. A case study was used to assess protectiveness of genus and family model estimates derived from either geometric mean or minimum taxa toxicity values for listed species. Models developed from the most sensitive value for each chemical were generally protective of the most sensitive species within predicted taxa, including listed species, and were more protective than geometric means models. ICE model estimates were compared to HC5 values derived from Species Sensitivity Distributions for the case study chemicals to assess protectiveness of the two approaches. ICE models provide robust toxicity predictions and can generate protective toxicity estimates for assessing contaminant risk to listed species.