Comparative toxicokinetics of organic micropollutants in freshwater crustaceans

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.

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.

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.

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.

Bioaccumulation, Biotransformation, and Multicompartmental Toxicokinetic Model of Antibiotics in Sea Cucumber (Apostichopus japonicus)

Extensive application of antibiotics leads to their ubiquitous occurrence in coastal aquatic environments. However, it remains largely unknown whether antibiotics can be bioaccumulated and biotransformed in major mariculture organisms such as sea cucumbers and toxicokinetic models for Echinodermata are lacking. In this study, laboratory exposure experiments on juvenile sea cucumber (Apostichopus japonicus) were performed for seven antibiotics (sulfadiazine, sulfamethoxazole, trimethoprim, enrofloxacin, ofloxacin, clarithromycin, and azithromycin). Field sea cucumber and surrounding seawater samples were also analyzed. Results show that the sea cucumbers tend to accumulate high concentrations of the antibiotics with kinetic bioconcentration factors (BCFs) up to 1719.7 L·kg^-1 for ofloxacin. The BCFs determined in the laboratory agree well with those estimated from the field measurements. Seven biotransformation products (BTPs) of the antibiotics were identified, four of which were not reported previously in aquatic organisms. The BTPs were mainly found in the digestive tract, indicating its high capacity in the biotransformation. A multicompartmental toxicokinetic model based on the principles of passive diffusion was developed, which can successfully predict time-course concentrations of the antibiotics in different compartments of the juvenile sea cucumbers. The findings may offer a scientific basis for assessing health risks and guiding healthy mariculture of sea cucumbers.

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.

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.

How Biotransformation Influences Toxicokinetics of Azole Fungicides in the Aquatic Invertebrate Gammarus pulex

Biotransformation is a key process that can greatly influence the bioaccumulation potential and toxicity of organic compounds. In this study, biotransformation of seven frequently used azole fungicides (triazoles: cyproconazole, epoxiconazole, fluconazole, propiconazole, tebuconazole and imidazoles: ketoconazole, prochloraz) was investigated in the aquatic invertebrate Gammarus pulex in a 24 h exposure experiment. Additionally, temporal trends of the whole body internal concentrations of epoxiconazole, prochloraz, and their respective biotransformation products (BTPs) were studied to gain insight into toxicokinetic processes such as uptake, elimination and biotransformation. By the use of high resolution tandem mass spectrometry in total 37 BTPs were identified. Between one (ketoconazole) and six (epoxiconazole) BTPs were identified per parent compound except for prochloraz, which showed extensive biotransformation reactions with 18 BTPs detected that were mainly formed through ring cleavage or ring loss. In general, most BTPs were formed by oxidation and conjugation reactions. Ring loss or ring cleavage was only observed for the imidazoles as expected from the general mechanism of oxidative ring openings of imidazoles, likely affecting the bioactivity of these BTPs. Overall, internal concentrations of BTPs were up to 3 orders of magnitude lower than that of the corresponding parent compound. Thus, biotransformation did not dominate toxicokinetics and only played a minor role in elimination of the respective parent compound, with the exception of prochloraz.