Mechanistic Understanding of the Synergistic Potential of Azole Fungicides in the Aquatic Invertebrate Gammarus pulex

Organ-specific biotransformation in salmonids: Insight into intrinsic enzyme activity and biotransformation of three micropollutants

Aquatic ecosystems continue to be threatened by chemical pollution. To what extent organisms are able to cope with chemical exposure depends on their ability to display mechanisms of defense across different organs. Among these mechanisms, biotransformation processes represent key physiological responses that facilitate detoxification and reduce the bioaccumulation potential of chemicals. Biotransformation does not only depend on the ability of different organs to display biotransformation enzymes but also on the affinity of chemicals towards these enzymes. In the present study, we explored the ability of different organs and of two freshwater fish to support biotransformation processes through the determination of in vitro phase I and II biotransformation enzyme activity, and their role in supporting intrinsic clearance and the formation of biotransformation products. Three environmentally relevant pollutants were evaluated: the polycyclic aromatic hydrocarbon (PAH) pyrene (as recommended by the OECD 319b test guideline), the fungicide azoxystrobin, and the pharmaceutical propranolol. Comparative studies using S9 sub-cellular fractions derived from the liver, intestine, gills, and brain of brown trout (Salmo trutta) and rainbow trout (Oncorhynchus mykiss) revealed significant phase I and II enzyme activity in all organs. However, organ- and species-specific differences were found. In brown trout, significant extrahepatic biotransformation was observed for pyrene but not for azoxystrobin and propranolol. In rainbow trout, the brain appeared to biotransform azoxystrobin. In this same species, propranolol appeared to be biotransformed by the intestine and gills. Biotransformation products could be detected only from hepatic biotransformation, and their profiles and formation rates displayed species-specific patterns and occurred at different magnitudes. Altogether, our findings further contribute to the current understanding of organ-specific biotransformation capacity, beyond the expression and activity of enzymes, and its dependence on specific enzyme-chemical interactions to support mechanisms of defense against exposure.

Evidence of strobilurin fungicides and their metabolites in Dongjiang River ecosystem, southern China: Bioaccumulation and ecological risks

Despite the widespread application of strobilurin fungicides (SFs) in agriculture, little is known about their distribution and bioaccumulation in aquatic ecosystems. In this study, the concentrations of 12 SFs and two of their metabolites were determined in abiotic (water and sediment; n = 83) and biotic (plant, algae, zooplankton, and fish; n = 123) samples collected from a subtropical freshwater ecosystem, namely, Dongjiang River wetland, in southern China. Among the 12 SFs measured, azoxystrobin (AZ) was the major fungicide found in surface water (median: 2.20 ng/L) and sediment (0.064 ng/g dry wt.). Azoxystrobin acid (AZ-acid), a metabolite of AZ, was the major analyte in the plant samples and had a median concentration at 0.36 ng/g dry wt. In algae and zooplankton, (Z)-metominostrobin was the predominant fungicide and had median concentrations of 3.52 and 5.55 ng/g dry wt., respectively. In fish muscle, dimoxystrobin (DIMO) was the major SF and had a median concentration of 0.47 ng/g dry wt. The bioconcentration factor (BCF) values of AZ-acid, trifloxystrobin (TFS), and pyraclostrobin (PYR) in algae and zooplankton and AZ-acid, PYR, TFS, TFS-acid, picoxystrobin, and DIMO in fish muscle exceeded 1000 L/kg (algae, zooplankton, and fish concentrations were expressed on a dry weight basis), suggesting that these fungicides can accumulate in biota. A positive association between log BCFs of SFs in fish and logKow of SFs and a negative correlation between log BCFs and the log solubility index were observed. Additionally, the risk quotient (RQ) was calculated to evaluate the potential ecotoxicological risk of SFs to different aquatic organisms (algae, zooplankton, and fish). The PYR and DIMO concentrations at 19 sampling sites had RQ values >0.1, indicating moderate ecotoxicological risks to aquatic organisms. This study is the first to document the widespread occurrence of SFs and their metabolites in aquatic ecosystems and to elucidate the bioaccumulation potential of SFs in aquatic organisms.

Single and joint toxicity of azoxystrobin and cyproconazole to Prochilodus lineatus: Bioconcentration and biochemical responses

Fungicides are widely used across the world to protect crops and their presence in freshwater systems is increasing. However, the evaluation of their potential impacts on non-target organisms is in the minority of studies related to pesticides. In the current research, the single and joint toxicity of azoxystrobin (AZX) and cyproconazole (CYP) was investigated in juvenile fish Prochilodus lineatus. In particular, we evaluated bioconcentration and biochemical responses following a short-term exposure to environmentally relevant concentrations of the fungicides (alone and in mixture). We also determined interactions between the biological responses when the two compounds were used in mixture. Our results demonstrate that AZX and CYP pose a risk to native freshwater fish by causing deleterious effects. Both compounds, alone and in mixture, bioaccumulated in P. lineatus and triggered neurotoxicity and changes in oxidative stress biomarkers in several organs. Moreover, muscle was a target tissue for these fungicides and a synergistic interaction was observed for the mixture. Due to the lack of studies in fish assessing the effects following exposure to AZX-CYP mixtures and considering a realistic exposure situation in agriculture-impacted water bodies, these findings provide new and relevant information for future studies.

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.

Comparative insight of pesticide transformations between river and wetland systems

The widespread use of pesticides threatens the environment and ecosystems. Despite the positive effects of plant protection products, pesticides also have unexpected negative effects on nontarget organisms. The microbial biodegradation of pesticides is one of the major pathways for reducing their risks at aquatic ecosystems. The objective of this study was to compare the biodegradability of pesticides in simulated wetland and river systems. Parallel experiments were conducted with 17 pesticides based on the OECD 309 guidelines. A comprehensive analytical method, such as target screening combined with suspect and non-target screening, was performed to evaluate the biodegradation via identification of transformation products (TPs) using LC-HRMS. As evidence of biodegradation, we identified 97 TPs for 15 pesticides. Metolachlor and dimethenamid had 23 and 16 TPs, respectively, including Phase II glutathione conjugates. The analysis of 16S rRNA sequences for microbials characterized operational taxonomic units. Rheinheimera and Flavobacterium, which have the potential for glutathione S-transferase, were dominant in wetland systems. Estimation of toxicity, biodegradability, and hydrophobicity using QSAR prediction indicated lower environmental risks of detected TPs. We conclude that the wetland system is more favorable for pesticide degradation and risk mitigation mainly attributed to the abundance and variety of the microbial communities.

Visualization of azoxystrobin penetration in wheat leaves using mass microscopy imaging

Fungicides must penetrate the internal tissues of plants to kill pathogenic fungi. Mass spectrometers have been used to confirm this penetration, but conventional mass spectrometric methods cannot distinguish the fungicides in different internal tissues owing to the extraction steps. However, matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) can detect the penetration of fungicides into leaf sections through direct analysis of the sample surfaces. Therefore, the objective of this study was to establish a method for visualizing fungicide penetration in wheat leaf cross sections using MALDI-MSI. The penetration of azoxystrobin from the epidermal to the internal tissue of the leaves was observed. Moreover, azoxystrobin accumulates in the cells around the vascular bundle. This study suggests that MSI can be useful for the evaluation of fungicide penetration in plant leaves.

Speed it up: How temperature drives toxicokinetics of organic contaminants in freshwater amphipods

The acceleration of global climate change draws increasing attention towards interactive effects of temperature and organic contaminants. Many studies reported a higher sensitivity of aquatic invertebrates towards contaminant exposure with increasing or fluctuating temperatures. The hypothesis of this study was that the higher sensitivity of invertebrates is associated with the changes of toxicokinetic processes that determine internal concentrations of contaminants and consequently toxic effects. Therefore, the influence of temperature on toxicokinetic processes and the underlying mechanisms were studied in two key amphipod species (Gammarus pulex and Hyalella azteca). Bioconcentration experiments were carried out at four different temperatures with a mixture of 12 exposure relevant polar organic contaminants. Tissue and medium samples were taken in regular intervals and analysed by online solid-phase extraction liquid chromatography high-resolution tandem mass spectrometry. Subsequently, toxicokinetic rates were modelled and analysed in dependence of the exposure temperature using the Arrhenius equation. An exponential relationship between toxicokinetic rates versus temperature was observed and could be well depicted by applying the Arrhenius equation. Due to a similar Arrhenius temperature of uptake and elimination rates, the bioconcentration factors of the contaminants were generally constant across the temperature range. Furthermore, the Arrhenius temperature of the toxicokinetic rates and respiration was mostly similar. However, in some cases (citalopram, cyprodinil), the bioconcentration factor appeared to be temperature dependent, which could potentially be explained by the influence of temperature on active uptake mechanisms or biotransformation. The observed temperature effects on toxicokinetics may be particularly relevant in non-equilibrated systems, such as exposure peaks in summer as exemplified by the exposure modelling of a field measured pesticide peak where the internal concentrations increased by up to fourfold along the temperature gradient. The results provide novel insights into the mechanisms of chemical uptake, biotransformation and elimination in different climate scenarios and can improve environmental risk assessment.

Combined Application of Tacrolimus with Cyproconazole, Hymexazol and Novel {2-(3-R-1H-1,2,4-triazol-5-yl)phenyl}amines as Antifungals: In Vitro Growth Inhibition and In Silico Molecular Docking Analysis to Fungal Chitin Deacetylase

Agents with antifungal activity play a vital role as therapeutics in health care, as do fungicides in agriculture. Effectiveness, toxicological profile, and eco-friendliness are among the properties used to select suitable substances. Furthermore, a steady supply of new agents with different modes of action is required to counter the well-known potential of human and phyto-pathogenic fungi to develop resistance against established antifungals. Here, we use an in vitro growth assay to investigate the activity of the calcineurin inhibitor tacrolimus in combination with the commercial fungicides cyproconazole and hymexazol, as well as with two earlier reported novel {2-(3-R-1H-1,2,4-triazol-5-yl)phenyl}amines, against the fungi Aspergillus niger, Colletotrichum higginsianum, Fusarium oxysporum and the oomycete Phytophthora infestans, which are notoriously harmful in agriculture. When tacrolimus was added in a concentration range from 0.25 to 25 mg/L to the tested antifungals (at a fixed concentration of 25 or 50 mg/L), the inhibitory activities were distinctly enhanced. Molecular docking calculations revealed triazole derivative 5, (2-(3-adamantan-1-yl)-1H-1,2,4-triazol-5-yl)-4-chloroaniline), as a potent inhibitor of chitin deacetylases (CDA) of Aspergillus nidulans and A. niger (AnCDA and AngCDA, respectively), which was stronger than the previously reported polyoxorin D, J075-4187, and chitotriose. The results are discussed in the context of potential synergism and molecular mode of action.

How to analyse and account for interactions in mixture toxicity with toxicokinetic-toxicodynamic models

The assessment of chemical mixture toxicity is one of the major challenges in ecotoxicology. Chemicals can interact, leading to more or less effects than expected, commonly named synergism and antagonism respectively. The classic ad hoc approach for the assessment of mixture effects is based on dose-response curves at a single time point, and is limited to identifying a mixture interaction but cannot provide predictions for untested exposure durations, nor for scenarios where exposure varies in time. We here propose a new approach using toxicokinetic-toxicodynamic modelling: The General Unified Threshold model of Survival (GUTS) framework, recently extended for mixture toxicity assessment. We designed a dedicated mechanistic interaction module coupled with the GUTS mixture model to i) identify interactions, ii) test hypotheses to identify which chemical is likely responsible for the interaction, and finally iii) simulate and predict the effect of synergistic and antagonistic mixtures. We tested the modelling approach experimentally with two species (Enchytraeus crypticus and Mamestra brassicae) exposed to different potentially synergistic mixtures (composed of: prochloraz, imidacloprid, cypermethrin, azoxystrobin, chlorothalonil, and chlorpyrifos). Furthermore, we also tested the model with previously published experimental data on two other species (Bombus terrestris and Daphnia magna) exposed to pesticide mixtures (clothianidin, propiconazole, dimethoate, imidacloprid and thiacloprid) found to be synergistic or antagonistic with the classic approach. The results showed an accurate simulation of synergistic and antagonistic effects for the different tested species and mixtures. This modelling approach can identify interactions accounting for the entire time of exposure, and not only at one time point as in the classic approach, and provides predictions of the mixture effect for untested mixture exposure scenarios, including those with time-variable mixture composition.

Combined toxicities of cadmium and five agrochemicals to the larval zebrafish (Danio rerio)

Different pollutants usually co-exist in the natural environment, and the ecological and health risk assessment of agrochemicals needs to be carried out based on the combined toxicological effects of pollutants. To examine the combined toxicity to aquatic organisms, the effects of cadmium (Cd) and five pesticides (acetamiprid, carbendazim, azoxystrobin, chlorpyrifos, and bifenthrin) mixture on zebrafish (Danio rerio) larvae were assessed. The data from the 96-h toxicity test indicated that bifenthrin possessed the highest toxicity to D. rerio with the LC₅₀ value of 0.15 mg L^-1, followed by chlorpyrifos (0.36 mg L^-1) and azoxystrobin (0.63 mg L^-1). Cd (6.84 mg L^-1) and carbendazim (8.53 mg L^-1) induced the intermediate toxic responses, while acetamiprid (58.39 mg L^-1) presented the lowest toxicity to the organisms. Pesticide mixtures containing chlorpyrifos and bifenthrin or acetamiprid and carbendazim showed synergistic impacts on the zebrafish. Besides, two binary combinations of Cd-acetamiprid and Cd-chlorpyrifos also displayed a synergistic effect on D. rerio. Our results offered a better idea of the mixed ecological risk assessment of Cd and different agricultural chemicals to aquatic organisms. Our findings better interpreted how the interaction between Cd and various agrochemicals changed their toxicity to aquatic vertebrates and provided valuable insights into critical impacts on the ecological hazard of their combinations.

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

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

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/g_ww. 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.

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

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

Individual performances and biochemical pathways as altered by field-realistic exposures of current-use fungicides and their mixtures in a non-target species, Gammarus fossarum

Persistent fungicides, which are widely applied to agricultural soils to protect crops, are frequently detected in freshwaters because of hydraulic transfer, possibly resulting in unintentional adverse effects on wildlife. However, the ecotoxicity of fungicides in aquatic species remains scarcely assessed at environmentally relevant concentrations, and there is scant information available on their combined effects. This study aims at investigating multi-level changes elicited by two currently used fungicides, boscalid and tebuconazole, in the amphipod Gammarus fossarum. In microcosms, gammarids were exposed for 72 h to fungicides tested individually or in binary mixture at 0.01, 0.1 and 1.0 μg/L to monitor individual performances (locomotion, respiration and amplexus formation) and biochemical parameters (involved in energy metabolism, growth, moulting and cell stress). This range of exposure concentrations was field-realistic and largely lower than local environmental quality standards for the protection of aquatic wildlife. Overall, results showed that single and combined exposures altered the mobility and respiratory activity of individuals. At the cellular level, boscalid inhibited energy-based biomarkers whereas tebuconazole led to cytotoxicity associated with reduced antioxidant defences. In binary mixtures, the biochemical responses were mainly driven by the presence of boscalid. Multi-variable analyses, integrating individual alterations and cellular impairments alike, confirmed the relevance of the multi-level approach in forecasting the toxicity of fungicides and their mixtures towards other aquatic species. This study demonstrates dose-related sublethal effects of fungicides on multiple functional traits in an ecosystem engineer under realistic exposure scenarios, highlighting the harmful signs of these toxicants. Such sublethal alterations could thus constitute reliable tools for the early diagnosis of the organisms' health and the ecological status of agriculturally impacted hydrosystems.

Grandmother's pesticide exposure revealed bi-generational effects in Daphnia magna

Man-made chemicals are a significant contributor to the ongoing deterioration of numerous ecosystems. Currently, risk assessment of these chemicals is based on observations in a single generation of animals, despite potential adverse intergenerational effects. Here, we investigate the effect of the fungicide prochloraz across three generations of Daphnia magna. We studied both the effects of continuous exposure over all generations and the effects of first-generation (F0) exposure on two subsequent generations. Effects at different levels of biological organization from genome-wide gene expression, whole organism metabolite levels, CYP enzyme activity and key phenotypic effects, such as reproduction, were monitored. Acclimation to prochloraz was found after continuous exposure. Following F0-exposure, embryonically exposed F1-offspring showed no significant effects. However, in the potentially germline exposed F2 animals, several parameters differed significantly from controls. A direct association between these F2 effects and the toxic mode of action of prochloraz was found, showing that chemicals can be harmful not only to the directly exposed generation, but also to prenatally exposed generations and in that way effects may even appear to skip a generation. This implies that current risk assessment practices are neglecting an important aspect of toxicity, such as delayed effects across generations due to a time gap between chemical exposure and emergence of effects.

Metabolomic Profiling and Toxicokinetics Modeling to Assess the Effects of the Pharmaceutical Diclofenac in the Aquatic Invertebrate Hyalella azteca

The exposure of ecologically critical invertebrate species to biologically active pharmaceuticals poses a serious risk to the aquatic ecosystem. Yet, the fate and toxic effects of pharmaceuticals on these nontarget aquatic invertebrates and the underlying mechanisms are poorly studied. Herein, we investigated the toxicokinetic (TK) processes (i.e., uptake, biotransformation, and elimination) of the pharmaceutical diclofenac and its biotransformation in the freshwater invertebrate Hyalella azteca. We further employed mass spectrometry-based metabolomics to assess the toxic effects of diclofenac on the metabolic functions of H. azteca exposed to environmentally relevant concentrations (10 and 100 μg/L). The TK results showed a quick uptake of diclofenac by H. azteca (maximum internal concentration of 1.9 μmol/kg) and rapid formation of the conjugate diclofenac taurine (maximum internal concentration of 80.6 μmol/kg), indicating over 40 times higher accumulation of diclofenac taurine than that of diclofenac in H. azteca. Depuration kinetics demonstrated that the elimination of diclofenac taurine was 64 times slower than diclofenac in H. azteca. Metabolomics results suggested that diclofenac inhibited prostaglandin synthesis and affected the carnitine shuttle pathway at environmentally relevant concentrations. These findings shed light on the significance of the TK process of diclofenac, especially the formation of diclofenac taurine, as well as the sublethal effects of diclofenac on the bulk metabolome of H. azteca. Combining the TK processes and metabolomics provides complementary insights and thus a better mechanistic understanding of the effects of diclofenac in aquatic invertebrates.

Biochemical and molecular responses of the Mediterranean mussel (Mytilus galloprovincialis) to short-term exposure to three commonly prescribed drugs

Pharmaceuticals represent a group of emerging contaminants. The short-term effect (3 and 7 days) of warfarin (1 and 10 mg L^-1), dexamethasone (0.392 and 3.92 mg L^-1) and imidazole (0.013 and 0.13 mg L^-1) exposure was evaluated on mussels (Mytilus galloprovincialis). Total antioxidant status, glutathione reductase, glutathione peroxidase (GPx) and superoxide dismutase enzyme activities, and the expression of genes involved in the xenobiotic response (ATP binding cassette subfamily B member 1 (abcb1) and several nuclear receptor family J (nr1j) isoforms), were evaluated. All nr1j isoforms are suggested to be the xenobiotic receptor orthologs of the NR1I family. All drugs increased GPx activity and altered the expression of particular nr1j isoforms. Dexamethasone exposure also decreased abcb1 expression. These findings raised some concerns regarding the release of these pharmaceuticals into the aquatic environment. Thus, further studies might be needed to perform an accurate environmental risk assessment of these 3 poorly studied drugs.

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

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

Identification of transformation products to characterize the ability of a natural wetland to degrade synthetic organic pollutants

Natural wetlands have been recognized as a natural reactor for degradation and elimination of environmental pollutants. The Upo Wetland, the largest inland wetland in Korea, is mainly surrounded by agricultural lands and it is susceptible to contamination from excess nutrient loads and synthetic organic contaminants (SOCs) (e.g., pesticides). The aim of this study was to identify major SOCs in the wetland and evaluate their degradation. We used high resolution mass spectrometry (HRMS) with a two-step analysis approach (i.e., 1^st analysis for target measurement along with suspect and non-target screening (SNTS) and 2^nd analysis for complimentary suspect screening) to identify and quantify the transformation products (TPs) of the identified parent SOCs. Quantitative analysis of 30 targets, mainly including pesticides, showed that fungicides were the major SOCs detected in the wetland, accounting for about 50% of the composition ratio of the total SOCs quantified. Orysastrobin occurred at the highest mean concentration (>700 ng/L), followed by two other fungicides, carbendazim and tricyclazole. The first analysis (SNTS) tentatively identified 39 TPs (30 by suspect, 9 by non-target screening) of 14 parent pesticides. Additionally, the second analysis (complimentary suspect screening) identified 9 more TPs. Among the 48 total TPs identified, 7 were confirmed with reference standards. The identification of the remaining TPs had a high confidence level (e.g., level 2 or 3). Regarding transport though the wetland, most TPs showed greater peak area ratios (i.e., the relative portion of chromatographic area of the TPs to the parent compound) at the outlet point of the wetland compared to the inlet point. The risk quotient, which was calculated using the concentrations of parent compounds, decreased toward the outlet, demonstrating the degradation capacity of the wetland. The estimates for biodegradability, hydrophobicity, and toxicity by an in-silico quantitative structure-activity relationship (QSAR) model indicated a lower half-life, lower logD_OW, and greater effect concentration for most TPs compared to the parent compounds. Based on these results, we conclude that natural wetlands play a role as an eco-friendly reactor for degrading SOCs to form numerous TPs that are lower risk than the parent compounds.

Retrospective screening of high-resolution mass spectrometry archived digital samples can improve environmental risk assessment of emerging contaminants: A case study on antifungal azoles

Environmental risk assessment associated with aquatic and terrestrial contamination is mostly based on predicted or measured environmental concentrations of a limited list of chemicals in a restricted number of environmental compartments. High resolution mass spectrometry (HRMS) can provide a more comprehensive picture of exposure to harmful chemicals, particularly through the retrospective analysis of digitally stored HRMS data. Using this methodology, our study characterized the contamination of various environmental compartments including 154 surface water, 46 urban effluent, 67 sediment, 15 soil, 34 groundwater, 24 biofilm, 41 gammarid and 49 fish samples at 95 sites widely distributed over the Swiss Plateau. As a proof-of-concept, we focused our investigation on antifungal azoles, a class of chemicals of emerging concern due to their endocrine disrupting effects on aquatic organisms and humans. Our results demonstrated the occurrence of antifungal azoles and some of their (bio)transformation products in all the analyzed compartments (0.1-100 ng/L or ng/g d.w.). Comparison of actual and predicted concentrations showed the partial suitability of level 1 fugacity modelling in predicting the exposure to azoles. Risk quotient calculations additionally revealed risk of exposure especially if some of the investigated rivers and streams are used for drinking water production. The case study clearly shows that the retrospective analysis of HRMS/MS data can improve the current knowledge on exposure and the related risks to chemicals of emerging concern and can be effectively employed in the future for such purposes.

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.

Characterizing biotransformation products and pathways of the flame retardant triphenyl phosphate in Daphnia magna using non-target screening

Triphenyl phosphate (TPHP), one of the organophosphate flame retardants, has been widely used in manufacturing, thereby causing a gradual increase in TPHP concentrations in aquatic environments. However, the information on the biotransformation mechanism of TPHP in invertebrates is lacking. The study identified the biotransformation products of TPHP in Daphnia magna, which showed particularly high toxicity in aquatic organisms, and determined the rates of depuration. Daphnia magna, a standard species for toxicity studies, was exposed to triphenyl phosphate and transferred to the pure medium. The biotransformation products of TPHP and its depuration rates were determined by liquid chromatography-high resolution mass spectrometry. Nine biotransformation products (five in the positive mode and four in the negative mode) of triphenyl phosphate were identified in D. magna. Based on the depuration ratio, the major biotransformation mechanism is estimated to be cysteine conjugation and sulfation. Certain biotransformation products (diphenyl phosphate, hydroxylated triphenyl phosphate, and thiol triphenyl phosphate) might induce toxicity in biota. The results could be used to predict main biotransformation processes and toxic products of organophosphate flame retardants in aquatic invertebrates.

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.

In vitro biotransformation of pharmaceuticals and pesticides by trout liver S9 in the presence and absence of carbamazepine

The aim of the present study was to develop (i) a technique for identifying metabolites of organic contaminants by using an in vitro system of trout S9 and liquid chromatography-high-resolution mass spectrometry-based identification method and (ii) to apply this technique to identify the interactive potential of carbamazepine on the formation rate of other metabolites. The pharmaceuticals carbamazepine and propranolol and the pesticides azoxystrobin, diazinon, and fipronil were selected as test contaminants. As a result, a total of ten metabolites were identified for the five parent substances, six of which were confirmed using reference standards. Metabolic reactions included hydroxylation, epoxidation, S-oxidation, and dealkylation. The metabolic transformation rate ranged from 0.2 to 3.5 pmol/mg protein/min/μmol substrate. In the binary exposure experiment with increasing carbamazepine concentration, the formation rates of diazinon and fipronil metabolites (MDI2 and MFP2, respectively) increased, while formation of metabolites of propranolol and azoxystrobin (MPR1, MPR2, MPR3, and MAZ1) slowed down. Meanwhile, S9 pre-exposed to carbamazepine produced diazoxon, a toxic metabolite of diazinon, and pyrimidinol, a less toxic metabolite, more rapidly. These results suggest that carbamazepine, a perennial environmental pollutant, might modulate the toxicity of other substances such as diazinon but further in vivo studies are needed.

Environmental Stress Increases Synergistic Effects of Pesticide Mixtures on Daphnia magna

Some widely used pesticide mixtures produce more than additive effects according to conventional combined effect models. However, synergistic effects have been so far generally observed at unrealistically high pesticide concentrations. Here, we used Daphnia magna as a test organism and investigated how food limitation-a common ecological stressor-affects the mixture toxicity of a pyrethroid insecticide and azole fungicide. We also compared three models regarding the prediction of mixture effects, including concentration addition (CA), effect addition (EA), and stress addition model (SAM). We revealed that especially under low food, the strength of synergism between esfenvalerate and prochloraz increased with an increasing concentration of prochloraz independent of the null model. Under high food conditions and at concentrations of prochloraz ≥32 μg/L, we observed a marginal synergistic effect with a model deviation ratio (MDR) = 2.1 at 32 μg/L prochloraz and 2.2 at 100 μg/L prochloraz when using CA as the null model. In contrast, the combination of both pesticides and food stress caused synergistic effects shown by an MDR = 10.9 even at 1 μg/L of prochloraz that is frequently detected in the environment. The combined effects of pesticides and food stress could be predicted best with the SAM that showed the lowest mean deviation between effect observation and prediction (mean deviation SAM = 16 [SD = 28], EA = 1072 [2105], CA = 1345 [2644]). We conclude that common environmental stressors can strongly increase the synergistic effects of toxicants. This knowledge is especially relevant considering current efforts to include the additional risk of pesticide mixtures and environmental stressors into the environmental risk assessment of pesticides.

A Bayesian framework for estimating parameters of a generic toxicokinetic model for the bioaccumulation of organic chemicals by benthic invertebrates: Proof of concept with PCB153 and two freshwater species

Toxicokinetic (TK) models are relevant and widely used to predict chemical concentrations in biological organisms. The importance of dietary uptake for aquatic invertebrates has been increasingly assessed in recent years. However, the model parameters are estimated on limited specific laboratory data sets that are bounded by several uncertainties. The aim of this study was to implement a Bayesian framework for simultaneously estimating the parameters of a generic TK model for benthic invertebrate species from all data collected. We illustrate our approach on the bioaccumulation of PCB153 by two species with different life traits and therefore exposure routes: Chironomus riparius larvae exposed to spiked sediment for 7 days and Gammarus fossarum exposed to spiked sediment and/or leaves for 7 days and then transferred to a clean media for 7 more days. The TK models assuming first-order kinetics were fitted to the data using Bayesian inference. The median model predictions and their 95% credibility intervals showed that the model fit the data well. From a methodological point of view, this paper illustrates that simultaneously estimating all model parameters from all available data by Bayesian inference, while considering the correlation between parameters and different types of data, is a real added value for TK modeling. Moreover, we demonstrated the ability of a generic TK model considering uptake and elimination routes as modules to add according to the availability of the data measured. From an ecotoxicological point of view, we show differences in PCB153 bioaccumulation between chironomids and gammarids, explained by the different life traits of these two organisms.

Imidacloprid induces adverse effects on fish early life stages that are more severe in Japanese medaka (Oryzias latipes) than in zebrafish (Danio rerio)

Neonicotinoids are widely used insecticides that have frequently been found in freshwater with concentrations ranging from ng to μg/L. It is known that these compounds impact non-target invertebrates, such as bees and gammaridae, in terms of toxicity and behavior, but impacts and species differences on vertebrates such as fish are little explored. The aim of this study was to investigate and compare the effects of one widely used neonicotinoid, imidacloprid, on development and behavior of two fish model species: Zebrafish (Danio rerio) and Japanese medaka (Oryzias latipes). Fish were exposed for 5 (zebrafish) and 14 (medaka) days from 0.2 to 2000 μg/L imidacloprid by aqueous exposure. Survival, development, behavior and histological features were monitored and organism-internal concentrations and biotransformation products measured. Imidacloprid caused sublethal effects in both species but the effects were much stronger in medaka with deformities, lesions and reduced growth being the most prominent impacts. Due to the overall longer time of development, time-integrated exposure of medaka was about 2-fold higher compared to zebrafish, potentially accounting for parts of the sensitivity differences. Our results underline the importance of taking species sensitivity differences into account especially when considering that medaka responded at imidacloprid concentrations that have been measured in the environment.

The synergistic potential of azole fungicides does not directly correlate to the inhibition of cytochrome P450 activity in aquatic invertebrates

The ability of azole fungicides to inhibit cytochrome P450 dependent metabolism is proposed to be the main mechanism for their synergizing effect on pyrethroid insecticide toxicity in aquatic invertebrates. This study investigates the correlation between inhibition strength and synergistic potential of azole fungicides in the crustacean Daphnia magna and the insect larvae Chironomus riparius. Inhibition strength was measured in vivo toward the cytochrome P450 catalysed conversion of 7-ethoxycoumarin to 7-hydroxycoumarin (ECOD). Synergistic potentials were determined as the ratio between predicted and observed toxicity of mixtures based on the model of concentration addition (CA) and independent action (IA). Azoles (n = 9-11) enhanced the toxicity of α-cypermethrin in D. magna (Synergy ratios CA: 0.8 - 16; IA: 1.1 - 22) and inhibited cytochrome P450 activity by different degrees (IC50: 0.0023 - 36 μM for D. magna and 0.08 - 24 μM for C. riparius). Inhibition strengths were strongly correlated in the two organisms (r: 0.937 p: 0.019 for triazoles and r: 0.903 p: 0.097 for imidazoles). Lipophilicity governed the inhibition strength of triazoles in both species (r > 0.9, p < 0.05). No correlation was observed between inhibition strengths and synergistic potentials. Several reasons for the apparent lack of correlation were discussed.

Metabolic Activity of Cytochrome P450s Towards Four Pyrethroids in Midgut Tissue From Locusta migratoria (Orthoptera: Acrididae)

Cytochrome P450 monooxygenases (P450s) play important roles in metabolizing various insecticides and often contribute to the development of insecticide resistance in insects and other arthropod species. The objective of this study was to compare the metabolism of four commonly used pyrethroids including deltamethrin, fluvalinate, fenvalerate, and permethrin in the midgut tissue of Locusta migratoria Linnaeus (Orthoptera: Acrididae) by using synergism bioassay and ultra-performance liquid chromatography (UPLC)-mass spectrometer (MS) analyses. Our study showed that piperonyl butoxide (PBO, P450 enzyme inhibitor) can significantly synergize the toxicity of deltamethrin, fluvalinate, and fenvalerate with synergism ratios ranging from 1.30 to 1.70 folds. Preincubations of the midgut tissue with PBO followed by incubations with each of the four pyrethroids resulted in significantly higher amounts of unmetabolized deltamethrin and fluvalinate than those in the control (preincubation without PBO) as well as preincubations with other two detoxification enzyme inhibitors. These results indicate that P450s play important roles in metabolizing deltamethrin and fluvalinate in the midgut tissue. Our further study using deltamethrin as a representative pyrethroid and UPLC-MS techniques confirmed that the reduced amount of deltamethrin in the control (preincubation without PBO) was due to the metabolism of deltamethrin to yield hydroxydeltamethrin which is a major metabolite produced by P450-mediated aromatic hydroxylation of deltamethrim. These results provide new insights into differential metabolic activity of P450s towards different pyrethroids in the midgut tissue of L. migratoria.

Bioaccumulation, Biotransformation, and Synergistic Effects of Binary Fungicide Mixtures in Hyalella azteca and Gammarus pulex: How Different/Similar are the Two Species?

Aquatic organisms are consistently exposed to a mixture of micropollutants that can bioaccumulate, undergo biotransformation, and may exert mixture effects. However, little is known on the underlying mechanisms and species-specificity. Herein we investigated bioaccumulation, biotransformation and synergistic effects of azole (i.e., prochloraz) and strobilurin (i.e., azoxystrobin) fungicides in the two aquatic invertebrate species, Hyalella azteca and Gammarus pulex. Bioaccumulation of azoxystrobin was similar, whereas bioaccumulation of prochloraz was slightly different in the two species but was still significantly below the REACH criteria for bioaccumulative substances. Similar biotransformation patterns were observed in both species, and only a few unique biotransformation reactions were detected in H. azteca such as malonyl-glucose and taurine conjugation. Toxicokinetic modeling additionally indicated that biotransformation is a more important elimination pathway in H. azteca. In mixtures, no-observed-adverse-effect levels of prochloraz decreased the LC₅₀s of azoxystrobin in both species which correlated well with increased internal azoxystrobin concentrations. This synergistic effect is partly due to the inhibition of cytochrome P450 monooxygenases by prochloraz which subsequently triggered the reduced biotransformation of azoxystrobin (lower by five folds in H. azteca). The largely similar responses in both species suggest that the easier-to-cultivate H. azteca is a promising representative of invertebrates for toxicity testing.

Comparisons of microsomal cytochrome P450 content and enzymatic activity towards selected model substrates and insecticides in different tissues from the migratory locust (Locusta migratoria)

Cytochrome P450 monooxygenases (P450s) are important enzymes for biotransformations of various endogenous and xenobiotic substances in various organisms. In this study, we examined microsomal P450 protein content and enzymatic activity in four major detoxification tissues dissected from fifth-instar nymphs of the migratory locust (Locusta migratoria). The highest microsomal P450 protein content was found in the gastric caeca (a part of the midgut), followed by the midgut, Malpighian tubules and fat bodies. Microsomal P450s showed the highest aromatic hydroxylation, O-dealkylation and O-dearylation activities towards six of the seven model substrates examined in the fat bodies. Although the gastric caeca showed the highest P450 protein content, the enzymatic activities towards six of the seven model substrates were the lowest in this tissue. Further, the midgut, gastric caeca and fat bodies showed significant metabolic activities towards two pyrethroid insecticides (deltamethrin and fluvalinate), but no significant activities towards the other four insecticides (malathion, chlorpyrifos, carbaryl and methoprene). These results support our conclusions: 1) total P450 protein content alone cannot be reliably used to predict its enzymatic activity, and 2) insect P450 enzymatic activity is both tissue and substrate dependent.

Lethal and sub-lethal effects of cyproconazole on freshwater organisms: a case study with Chironomus riparius and Dugesia tigrina

The fungicide cyproconazole (CPZ) inhibits the biosynthesis of ergosterol, an essential sterol component in fungal cell membrane and can also affect non-target organisms by its inhibitory effects on P450 monooxygenases. The predicted environmental concentration of CPZ is up to 49.05 μg/L and 145.89 μg/kg in surface waters and sediments, respectively, and information about CPZ toxicity towards non-target aquatic organisms is still limited. This study aimed to address the lack of ecotoxicological data for CPZ, and thus, an evaluation of the lethal and sub-lethal effects of CPZ was performed using two freshwater invertebrates (the midge Chironomus riparius and the planarian Dugesia tigrina). The estimated CPZ 48 h LC₅₀ (95% CI) was 17.46 mg/L for C. riparius and 47.38 mg/L for D. tigrina. The emergence time (EmT₅₀) of C. riparius was delayed by CPZ exposure from 0.76 mg/L. On the other hand, planarians showed higher tolerance to CPZ exposure. Sub-lethal effects of CPZ on planarians included reductions in locomotion (1.8 mg/L), delayed photoreceptors regeneration (from 0.45 mg/L), and feeding inhibition (5.6 mg/L). Our results confirm the moderate toxicity of CPZ towards aquatic invertebrates but sub-lethal effects observed also suggest potential chronic effects of CPZ with consequences for population dynamics.