Uptake, biotransformation and elimination of selected pharmaceuticals in a freshwater invertebrate measured using liquid chromatography tandem mass spectrometry

Predicting the Accumulation of Ionizable Pharmaceuticals and Personal Care Products in Aquatic and Terrestrial Organisms

The extent to which chemicals bioaccumulate in aquatic and terrestrial organisms represents a fundamental consideration for chemicals management efforts intended to protect public health and the environment from pollution and waste. Many chemicals, including most pharmaceuticals and personal care products (PPCPs), are ionizable across environmentally relevant pH gradients, which can affect their fate in aquatic and terrestrial systems. Existing mathematical models describe the accumulation of neutral organic chemicals and weak acids and bases in both fish and plants. Further model development is hampered, however, by a lack of mechanistic insights for PPCPs that are predominantly or permanently ionized. Targeted experiments across environmentally realistic conditions are needed to address the following questions: (1) What are the partitioning and sorption behaviors of strongly ionizing chemicals among species? (2) How does membrane permeability of ions influence bioaccumulation of PPCPs? (3) To what extent are salts and associated complexes with PPCPs influencing bioaccumulation? (4) How do biotransformation and other elimination processes vary within and among species? (5) Are bioaccumulation modeling efforts currently focused on chemicals and species with key data gaps and risk profiles? Answering these questions promises to address key sources of uncertainty for bioaccumulation modeling of ionizable PPCPs and related contaminants. Environ Toxicol Chem 2024;43:502-512. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Sulfadiazine chlorination disinfection by-products in constructed wetlands: Identification of biodegradation products and inference of transformation pathways

Disinfection by-products (DBPs) formed from chlorination of antibiotics have greater toxicity than their parent compounds. Herein, this study investigated the biotransformation process of sulfadiazine Cl-DBPs in constructed wetlands (CWs). Results showed that, S atom on sulfonyl group, and N atoms on primary and secondary amine groups were the most reactive sites of sulfadiazine molecule. S1-N4 and S1-C8 of sulfadiazine are the most vulnerable bonds to cleave, followed by C14-N4 and C11-N5 bonds. In the chlorination process, sulfadiazine went through C-N bond cleavage, N-reductive alkylation, halogenation, and desulfonation to produce two aromatic Cl-DBPs. In the biodegradation process in CWs, sulfadiazine Cl-DBPs went through processes mainly including dechlorination, S-N bond cleavage, aniline-NH2 oxidation, desulfonation, phenol-OH oxidation, benzene ring cleavage, C-N bond cleavage, and β-oxidation of fatty acids under the action of a variety of oxidoreductases and hydrolases, during which a total of ten biodegradation products was identified. Moreover, sulfadiazine affected the biodegradation rather than the adsorption process in CWs. The two aromatic sulfadiazine Cl-DBPs had much higher bioaccumulation potentials than their parent sulfadiazine, but for the ten biodegradation products of sulfadiazine Cl-DBPs in CWs, 70% and almost 100% of them had lower bioaccumulation potentials than sulfadiazine and their parent sulfadiazine Cl-DBPs, respectively. The CWs were effective in reducing the environmental risk of sulfadiazine Cl-DBPs.

Transformation pathways of enrofloxacin chlorination disinfection by-products in constructed wetlands

Antibiotic residues and their chlorinated disinfection by-products (Cl-DBPs) have adverse effects on organisms in aquaculture water. Taking enrofloxacin (ENR) as target antibiotic, this study investigated the degradation and transformation of ENR Cl-DBPs in constructed wetlands (CWs). Results showed that, ENR and its Cl-DBPs affected the biodegradation of CWs at the preliminary stage, but did not affect the adsorption by plant roots, substrates, and biofilms. The piperazine group of ENR had great electronegativity, and was prone to electrophilic reactions. The carboxyl on quinolone group of ENR had strong nucleophilicity, and was prone to nucleophilic reactions. C atoms with significant negative charges on the aromatic structure of quinolone group were prone to halogenation. During the chlorination of ENR, one pathway was the reaction of quinolone group, in which nucleophilic substitution reaction by chlorine occurred at C26 atom on carboxyl group, then halogenation occurred under the action of Cl+ at C17 site on the aromatic ring; the other pathway was the reaction of piperazine group, in which N7 atom was firstly attacked by HOCl, resulting in piperazine ring cleavage, then followed by deacylation, dealkylation, and halogenation. During the biodegradation of ENR Cl-DBPs, the reactivity of piperazine structure was strong, especially at N6, N7, C13, and C14 sites, while the ring structure of quinolone group was quite stable, and only occurred decyclopropyl at N5 site. Overall, the biodegradation of ENR Cl-DBPs in CWs went through processes including piperazine ring cleavage, tertiary amine splitting, dealkylation, and aldehyde oxidation under the action of coenzymes, in which metabolites such as ketones, aldehydes, carboxylic acids, amides, primary amines, secondary amines, tertiary amines and acetaldehyde esters were produced. Most ENR Cl-DBPs had greater bioaccumulation potential and stronger toxicity than their parent compound, fortunately, CWs effectively reduced the environmental risk of ENR Cl-DBPs through the cooperation of adsorption and biodegradation.

First insight into the sex-dependent accumulation, tissue distribution and potential toxicities of 2-ethylhexyl diphenyl phosphate and its metabolites in adult zebrafish

The 2-ethylhexyl diphenyl phosphate (EHDPHP), a primary organophosphorus flame retardant used in various industrial products, is prone to biotransformation. However, there is a knowledge gap on the sex- and tissue-specific accumulation and potential toxicities of EHDPHP (M1) and its metabolites (M2-M16). In this study, adult zebrafish (Danio rerio) were exposed to EHDPHP (0, 5, 35 and 245 µg/L) for 21-day, which was followed by 7-day depuration. The bioconcentration factor (BCF) of EHDPHP in female zebrafish was 26.2 ± 7.7% lower than in males due to the lower uptake rate (k_u) while higher depuration rate (k_d) in the females. The regular ovulation and higher metabolic efficiency promoted elimination from female zebrafish, thus leading to much less (28-44%) accumulation of ∑(M1-M16) in female zebrafish. They exhibited the highest accumulation in the liver and intestine in both sexes, which might be regulated by tissue-specific transporters and histones evidenced by molecular docking results. Intestine microbiota analysis further revealed that female zebrafish were more susceptible to EHDPHP exposure, with more significant changes in phenotype number and KEGG pathways in female than male fish. Disease prediction results suggested that EHDPHP exposure might cause cancers, cardiovascular diseases as well as endocrine disorders in both sexes. These results provide a comprehensive understanding of the sex-dependent accumulation and toxicity of EHDPHP and its metabolites.

Uptake and depuration of three common antibiotics in benthic organisms: Sea cucumber (Holothuria tubulosa), snakelocks anemone (Anemonia sulcata) and beadlet anemone (Actinia equina)

Antibiotics are widely used drugs in human and veterinary medicine, which has attracted great attention in relation to the development of bacterial resistance, currently a problem of great concern for governments and states, as it is related to the resurgence of infectious diseases already eradicated. Understanding the bioaccumulation of antibiotics in aquatic organisms is an important key to understanding their risk assessment. The present study was designed to study the bioaccumulation of target antibiotics in relevant organisms inhabiting benthic marine environments. The uptake and elimination of ciprofloxacin (CIP), sulfamethoxazole (SMX) and trimethoprim (TMP) were investigated in sea cucumbers (Holothuria tubulosa), snakelock anemone (Anemonia sulcata) and beadlet anemone (Actinia equina) under controlled laboratory conditions. The results show that antibiotics have a particular trend over time during all periods of absorption and depuration. The tissue distribution of antibiotics in sea cucumber is strongly influenced by the structure of the compounds, while CIP is concentrated in the body wall; TMP is concentrated in the digestive tract. Two different approaches were used to estimate bioconcentration factors (BCFs) in different animal models, based on toxicokinetic data and measured steady-state concentrations. The BCF ranges were 456-2731 L/kg, 6-511 L/kg and 9-100 L/kg for TMP, CIP and SMX, respectively. The estimated BCF values obtained classify TMP as cumulative in A. equina and H. tubulosa, underlining the potential bioconcentration in these marine organisms. A correlation was observed between the BCFs of the target antibiotics and the octanol-water distribution coefficient (D_ow) (r² > 0.7). The animal-specific BCF followed the order of beadlet anemone > sea cucumber > snakelock anemone.

Enantiospecific Uptake and Depuration Kinetics of Chiral Metoprolol and Venlafaxine in Marine Medaka (Oryzias melastigma): Tissue Distribution and Metabolite Formation

The increasing use of chiral pharmaceuticals has led to their widespread presence in the environment. However, their toxicokinetics have rarely been reported. Therefore, the tissue-specific uptake and depuration kinetics of two pairs of pharmaceutical enantiomers, S-(-)-metoprolol versus R-(+)-metoprolol and S-(+)-venlafaxine versus R-(-)-venlafaxine, were studied in marine medaka (Oryzias melastigma) during a 28-day exposure and 14-day clearance period. The toxicokinetics of the studied pharmaceuticals, including uptake and depuration rate constants, depuration half-life (t_1/2), and bioconcentration factor (BCF), were reported for the first time. The whole-fish results demonstrated a higher S- than R-venlafaxine bioaccumulation potential, whereas no significant difference was observed between S- and R-metoprolol. O-desmethyl-metoprolol (ODM) and α-hydroxy-metoprolol (AHM) were the main metoprolol metabolites identified by suspect screening, and the ratios of ODM to AHM were 3.08 and 1.35 for S- and R-metoprolol, respectively. N,O-Didesmethyl-venlafaxine (NODDV) and N-desmethyl-venlafaxine (NDV) were the main venlafaxine metabolites, and the ratios of NODDV to NDV were 1.55 and 0.73 for S- and R-venlafaxine, respectively. The highest tissue-specific BCFs of the four enantiomers were all found in the eyes, meriting in-depth investigation.

Behavioral alterations and gills damage in Mytilus galloprovincialis exposed to an environmental concentration of delorazepam

Psychoactive compounds, and benzodiazepines (BZPs) in particular, represent an important class of emerging pollutants due to their large (ab)use and high resistance to degradation. Nowadays it is known that sewage treatment does not completely eliminate these substances and, therefore, BZPs and their metabolites reach concern levels in most aquatic environments all over Europe, ranging from µg/L to ng/L. In this study, we investigated the effects of delorazepam on Mytilus galloprovincialis, a model organism in toxicity testing and a key species in coastal marine ecosystems. Given its psychoactive activity, the study primarily addressed discovering the effects on behavior, by conventional valve opening and closure tests. Possible cytotoxic activity was also investigated by analyzing valve abductor muscles, gills histology, and correlated oxygen consumption. Results demonstrate negative effects on mussel behavior, interference with metabolism, and alteration of gill morphology and protein content. In conclusion, delorazepam confirms its toxicity to aquatic environments, highlighting the possibility that BZDs can ultimately affect the structure of the food web and the functions of the coastal ecosystems.

Intestinal microbiota perturbations in the gastropod Trochus niloticus concurrently exposed to ocean acidification and environmentally relevant concentrations of sulfamethoxazole

Ocean acidification (OA) and antibiotic pollution pose severe threats to the fitness of keystone species in marine ecosystems. However, the combined effects of OA and antibiotic pollution on the intestinal microbiota of marine organisms are still not well known. In this study, we exposed the herbivorous gastropod Trochus niloticus, a keystone species to maintains the stability of coral reef ecosystems, to acidic seawater (pH 7.6) and/or sulfamethoxazole (SMX, 100 ng/L, 1000 ng/L) for 28 days and determined their impacts on (1) the accumulation of SMX in the intestine of T. niloticus; (2) the characteristics of the intestinal microbiota in T. niloticus; (3) the relative abundances of sulfonamide resistance genes (i.e., sul1 and sul2) and intI1 in the intestinal microbiota of T. niloticus. Our results show that OA exposure leads to dramatic microbiota dysbiosis in the intestine of T. niloticus, including changes in bacterial community diversity and structure, decreased abundances of dominant species, existences of characteristic taxa, and altered functional predictions. In addition, SMX exposure at environmentally relevant concentrations had little effect on the intestinal microbiota of T. niloticus, whether in isolation or in combination with OA. However, after exposure to the higher SMX concentration (1000 ng/L), the accumulation of SMX in the intestine of T. niloticus could induce an increase in the copies of sul2 in the intestinal microbiota. These results suggest that the intestinal health of T. niloticus might be affected by OA and SMX, which might lead to fitness loss of the keystone species in coral reef ecosystems.

The influence of bromide and iodide ions on the sulfamethoxazole (SMX) halogenation during chlorination

Disinfection by-products (DBPs) were produced during the chlorination process, posing a threat to drinking water safety and human health. In the presence of bromide and iodide ions, brominated and iodinated DBPs will be generated, which might be more toxic than the parent compound. However, there are few studies on brominated and iodinated DBPs of antibiotics. Therefore, in this study, the fates of sulfamethoxazole (SMX) during chlorination in different systems (Blank; SMX + NaClO; SMX+ NaClO+ Br^-; SMX+ NaClO+I^-; SMX+ NaClO+ Br^- + I^-) were investigated. In different systems, all the reaction followed a pseudo-first-order kinetics, while the reaction rates of NaClO with SMX were different, the reaction rates were in order of SMX + NaClO + Br^- + I^- > SMX + NaClO + Br^- > SMX + NaClO + I^- > SMX + NaClO. When Br^- and I^- existed simultaneously, the reaction rate was the fastest. Iodide played an important role in oxidation and promoted the chlorination of SMX. SMX mainly underwent S-C cleavage, S-N hydrolysis, desulfonation, and substitution reactions. Nine disinfection by-products, including three reported for the first time, were identified using a non-targeted approach, and degradation pathways were proposed. Furthermore, EPI Suite software was applied to predict the environmental accumulation potential and environmental persistence of the degradation products. The results indicated that SMX and degradation products had little environmental accumulative potential and environmental persistence.

Prioritization based on risk assessment to study the bioconcentration and biotransformation of pharmaceuticals in glass eels (Anguilla anguilla) from the Adour estuary (Basque Country, France)

The presence of contaminants of emerging concern in the aquatic environment directly impacts water-living organisms and can alter their living functions. These compounds are often metabolized and excreted, but they can also be accumulated and spread through the food chain. The metabolized contaminants can also lead to the formation of new compounds with unknown toxicity and bioaccumulation potential. In this work, we have studied the occurrence, bioconcentration, and biotransformation of CECs in glass eels (Anguilla anguilla) using UHPLC-HRMS. To select the target CECs, we first carried out an environmental risk assessment of the WWTP effluent that releases directly into the Adour estuary (Bayonne, Basque Country, France). The risk quotients of every detected contaminant were calculated and three ecotoxicologically relevant contaminants were chosen to perform the exposure experiment: propranolol, diazepam, and irbesartan. An experiment of 14 days consisting of 7 days of exposure and 7 days of depuration was carried out to measure the bioconcentration of the chosen compounds. The quantitative results of the concentrations in glass eel showed that diazepam and irbesartan reached BCF ≈10 on day 7, but both compounds were eliminated after 7 days of depuration. On the other hand, propranolol's concentration remains constant all along with the experiment, and its presence can be detected even in the non-exposed control group, which might suggest environmental contamination. Two additional suspect screening strategies were used to identify metabolization products of the target compounds and other xenobiotics already present in wild glass eels. Only one metabolite was identified, nordiazepam, a well-known diazepam metabolite, probably due to the low metabolic rate of glass eels at this stage. The xenobiotic screening confirmed the presence of more xenobiotics in wild glass eels, prominent among them, the pharmaceuticals exemestane, primidone, iloprost, and norethandrolone.

Bioaccumulation and bioamplification of pharmaceuticals and endocrine disruptors in aquatic insects

Environmental fate of emerging contaminants such as pharmaceuticals and endocrine disrupting compounds at the aquatic terrestrial boundary are largely unexplored. Aquatic insects connect aquatic and terrestrial food webs as their life cycle includes aquatic and terrestrial life stages, thus they represent an important inter-habitat linkage not only for energy and nutrient flow, but also for contaminant transfer to terrestrial environments. We measured the concentrations of pharmaceuticals and endocrine disrupting compounds in the larval and adult tissues (last larval stages and teneral adults) of five Odonata species sampled in a wastewater-impacted river, in order to examine their bioaccumulation and bioamplification at different taxonomic levels. Twenty different compounds were bioaccumulated in insect tissues, with majority having higher concentrations (up to 90% higher) in aquatic larvae compared to terrestrial adults (reaching 88 ng/g for 1H-benzotriazole). However, increased concentration in adults was observed for seven compounds in at least one suborder (41% of the accumulated), confirming contaminants bioamplification across the metamorphosis. Both, bioaccumulation and bioamplification differed at various taxa levels; the order (Odonata), suborder (Anisoptera and Zygoptera) and species level. Highest variability was observed between Anisoptera and Zygoptera, due to the underlying differences in their ecology. Generally, Zygoptera had higher concentrations of contaminants in both larvae and adults. Additionally, we aimed at predicting effects of contaminant properties on bioaccumulation and bioamplification patterns using the commonly used physicochemical and pharmacokinetic descriptors on both order and suborder levels, however, neither of the two processes could be consistently predicted with simple linear models. Our study highlights the importance of taxonomy in studies aiming at advancing the understanding of contaminant exchange between aquatic and terrestrial food webs, as higher taxonomic categories include ecologically diverse groups, whose contribution to "the dark side of subsidies" could substantially differ.

Simultaneous multi-residue analytical method for anesthetics and sedatives in seafood samples by LC-ESI/MSMS

A simultaneous multi-residue analytical method for 27 regulated and unregulated anesthetics and sedatives in seafood using liquid chromatography/tandem mass spectrometry with electrospray ionization (LC-ESI/MSMS) was developed and tested on flatfish, eels, and shrimp. To optimize the efficiency of the method, extraction and clean-up procedures with various solvents and sorbents were tested. The most efficient pretreatment methods were extraction using acetonitrile (ACN) only (for flatfish and eel) and 0.1 % ammonium acetate in ACN (for shrimp). Validation was performed under the guidelines of CODEX Alimentarius (CAC/GL-71) and the Korean National Institute of Food and Drug Safety Evaluation (NIFDS). The limit of detection (LOD), limit of quantification (LOQ), accuracy, and precision for all compounds ranged from 0.0002 to 0.002 mg/kg, 0.0005-0.005 mg/kg, 64.7-112.5 %, and 1.0-8.6 %, respectively. The coefficient of determination (linearity, R²) was over 0.98. Therefore, the method meets the requirements of both the domestic and international guidelines.

Transcriptome responses in copepods Calanus finmarchicus, Calanus glacialis and Calanus hyperboreus exposed to phenanthrene and benzo[a]pyrene

Arctic and sub-arctic pelagic organisms can be exposed to effluents and spills from offshore petroleum-related activities and thus it is important to understand how they respond to crude oil related contaminants such as polycyclic aromatic hydrocarbons (PAHs). The copepod species Calanus finmarchicus, Calanus glacialis and Calanus hyperboreus represent key links in the arctic marine food web. We performed a transcriptome analysis of the three species exposed to phenanthrene (Phe) and benzo[a]pyrene (BaP) representing low and high molecular weight PAHs, respectively. Differential expression of several genes involved in many cellular pathways was observed after 72 h exposure to Phe (0.1 μM) and BaP (0.1 μM). In C. finmarchicus and C. glacialis, the exposure resulted in up-regulation of genes encoding enzymes in xenobiotic biotransformation, particularly the phase II cytosolic sulfonation system that include 3'-phosphoadenosine 5'-phosphosulfate synthase (PAPSS) and sulfotransferases (SULTs). The sulfonation pathway genes were more strongly induced by BaP than Phe in C. finmarchicus and C. glacialis but were not affected in C. hyperboreus. However, a larger number of genes and pathways were modulated in C. hyperboreus by the PAHs including genes encoding xenobiotic biotransformation and lipid metabolism enzymes, suggesting stronger responses in this species. The results suggest that the cytosolic sulfonation is a major phase II conjugation pathway for PAHs in C. finmarchicus and C. glacialis. Some of the biotransformation systems affected are known to be involved in metabolism of endogenous compounds such as ecdysteroids, which may suggest potential interference with physiological and developmental processes of the copepod species.

Understanding the bioaccumulation of pharmaceutical active compounds by clams Ruditapes decussatus exposed to a UWWTP discharge

Twenty-four pharmaceutical active compounds (PhACs) were evaluated in the soft tissues of clams Ruditappes decussatus exposed along a 1.5-km dispersal gradient of the treated effluent from an urban wastewater treatment plant discharging in Ria Formosa, and compared with those in the marine waters and discharged effluents. The clams were exposed for 1 month, in June-July 2016, 2017 and 2018. PhACs were quantified by high performance liquid chromatography coupled to tandem mass spectrometry after the quick, easy, cheap, effective, rugged and safe (QuEChERS) method (clams) or solid-phase extraction (water samples). The most representative PhACs in the effluents and receiving waters (regardless of the tidal dilution effect) were diclofenac, carbamazepine and caffeine (on average ≤ 2 μg/L) and only caffeine exhibited significant inter-annual differences, with higher values in 2017. In turn, the most bioaccumulated PhACs in clams were caffeine (0.54-27 ng/g wet weight, significantly higher in 2016) and acetaminophen (0.37-3.7 ng/g wet weight, significant lower in 2016). A multivariate principal component analysis showed (i) PhAC bioaccumulation primarily depended on biotic factors (clams length and weight), (ii) PhAC physicochemical properties Log K_ow, pKa and water solubility interplaying with water abiotic variables were more relevant for explaining data variability in water than the physical dilution/tidal mixing, (iii) this process, reflected by the salinity gradient, had a tertiary role in data variation, responsible for spatial discrimination of marine waters. This study provides a better understanding of PhACs bioaccumulation by clams Ruditapes decussatus in real environmental conditions, under the influence of urban treated effluent dispersal in Ria Formosa coastal lagoon, a major producer of bivalves, ultimately disentangling key factors of PhAC bioaccumulation.

Water temperature affects the biotransformation and accumulation of a psychoactive pharmaceutical and its metabolite in aquatic organisms

Pharmaceutically active compounds (PhACs) have been shown to accumulate in aquatic and riparian food-webs. Yet, our understanding of how temperature, a key environmental factor in nature, affects uptake, biotransformation, and the subsequent accumulation of PhACs in aquatic organisms is limited. In this study, we tested to what extent bioconcentration of an anxiolytic drugs (temazepam and oxazepam) is affected by two temperature regimes (10 and 20 °C) and how the temperature affects the temazepam biotransformation and subsequent accumulation of its metabolite (oxazepam) in aquatic organisms. We used European perch (Perca fluviatilis) and dragonfly larvae (Sympetrum sp.), which represent predator and prey species of high ecological relevance in food chains of boreal and temperate aquatic ecosystems. Experimental organisms were exposed to target pharmaceuticals at a range of concentrations (0.2-6 µg L^-1) to study concentration dependent differences in bioconcentration and biotransformation. We found that the bioconcentration of temazepam in perch was significantly reduced at higher temperatures. Also, temperature had a strong effect on temazepam biotransformation in the fish, with the production and subsequent accumulation of its metabolite (oxazepam) being two-fold higher at 20 °C compared to 10 °C. In contrast, we found no temperature dependency for temazepam bioconcentration in dragonfly larvae and no detectable biotransformation of the parent compound that would result in measurable concentrations of oxazepam in this organism. Our results highlight that while organisms may share the same aquatic ecosystem, their exposure to PhACs may change differently across temperature gradients in the environment.

The added value of Bayesian inference for estimating biotransformation rates of organic contaminants in aquatic invertebrates

Toxicokinetic (TK) models refer to the process of contaminant bioaccumulation as a balance between rate of uptake from different sources (e.g., water or food), and rate of elimination via different processes such as excretion, growth and/or biotransformation. Biotransformation can considerably modify the fate of chemicals in an organism, especially their bioavailability, residence time, and toxicity. Invertebrate models generally neglect this process as they assume a low metabolic activity. However, some species such as Gammarus sp. amphipods are able to metabolize a vast range of organic compounds. Some recent TK models include biotransformation, but they prove limited for estimating related parameters by giving negative values and/or large uncertainties for biotransformation rate(s). Here we propose a generic TK model accounting for biotransformation using a Bayesian framework for simultaneously estimating the parameters. We illustrated the added value of our method by fitting this generic TK model to 22 published datasets of several benthic invertebrate species exposed to different chemicals. All parameters are estimated simultaneously for all datasets and showed narrow estimates. Furthermore, the median model predictions and their 95% credibility intervals showed that the model confidently fitted the data. In most cases the uncertainties around biotransformation rate(s) were reduced in comparison to the original studies. From a methodology standpoint, this paper reflects that Bayesian inference has real added value for simultaneously estimating all TK parameters for parent chemicals and their metabolite(s) based on all available data, while accounting for different types of data and the correlation between parameters. Bayesian inference was able to overcome the limits of previous methods, since no parameters were fixed and no irrelevant negative values were obtained. Moreover, the 95% credibility intervals around model predictions, which are core uncertainties for Environmental Risk Assessment, were easily acquired.

Developing in vitro models to assess fish gill excretion of emerging contaminants

Advances in analytical methods have enabled the detection of emerging contaminants at ever lower concentrations in freshwaters. However, such measurements need to be linked to effect-based assays to identify risks. The bioconcentration factor (BCF) forms part of a chemical's environmental risk assessment (ERA), and current regulatory testing guidelines to calculate fish BCFs use hundreds of fish per chemical. Due to ethical concerns a reduction in the numbers of animals used is desired, and there is a need to identify in vitro or in silico alternatives which meet regulatory acceptance. This study describes the successful demonstration of a FIsh Gill Cell culture System (FIGCS) to assess an often overlooked parameter in pharmacokinetics: the excretion of drugs across the gill. The FIGCS tolerates the application of natural waters on its apical surface, mimicking the situation of the live fish, and thus in combination with advanced analytical methods, offers an opportunity to take lab-based testing used for ERA, such as compound uptake, biotransformation or excretion directly into field for validation with natural waters. Here we used the basic drug propranolol and the acidic ibuprofen as a demonstration of the FIGCS utility in three separate experiments. Excretion across the apical membrane showed saturation kinetics, suggesting the involvement of carrier-mediated processes. Both propranolol and ibuprofen were excreted across the epithelium from the media (internal blood equivalent) to the water, with ibuprofen excretion being considerably slower than propranolol excretion. Further studies indicate that ibuprofen may be complexing with fetal bovine serum (FBS) reducing bioavailability; in contrast propranolol efflux rate was unaffected, indicating that drugs behave differently in the presence of FBS and other plasma proteins. A key issue in future ERA is to better understand the effects of mixtures of different pollutant classes found in environmental samples, and this model offers an ethical path to do this.

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.

Bioconcentration and behavioral effects of four benzodiazepines and their environmentally relevant mixture in wild fish

We studied the adverse effects of four benzodiazepines frequently measured in European surface waters. We evaluated bioaccumulation potential of oxazepam, bromazepam, temazepam, and clobazam in freshwater fish species - perch (Perca fluviatilis) and we conducted a series of behavioral trials to assess their potential to alter boldness, activity, and social behavior. All selected endpoints were studied individually for each target benzodiazepine and as a mixture of all tested compounds to assess possible combinatory effects. We used a three-dimensional automated tracking system to quantify the fish behavior. The four compounds bioconcentrated differently in fish muscle (temazepam > clobazam > oxazepam > bromazepam) at high exposure (9.1, 6.9, 5.7, 8.1 µg L^-1, respectively) and low exposure (0.5, 0.5, 0.3, 0.4 µg L^-1, respectively) concentrations. A significant amount of oxazepam was also measured in fish exposed to temazepam, most likely because of the metabolic transformation of temazepam within the fish. Bromazepam, temazepam, and clobazam significantly affected fish behavior at high concentration, while no statistically significant changes were registered for oxazepam. The studied benzodiazepines affected behavior in combination, because the mixture treatment significantly changed several important behavioral traits even at low concentration, while no single compound exposure had such an effect at that dose. Based on our results, we conclude that effects of pharmaceuticals on aquatic environments could be underestimated if risk assessments only rely on the evaluation of single compounds. More studies focused on the combinatory effects of environmentally relevant mixtures of pharmaceuticals are necessary to fill the gaps in this knowledge.

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.

Biomonitoring of pesticides, pharmaceuticals and illicit drugs in a freshwater invertebrate to estimate toxic or effect pressure

Multiple classes of environmental contaminants have been found in aquatic environments, globally. Understanding internalised concentrations in the organism could further improve the risk assessment process. The present study is concerned with the determination of several contaminant classes (107 compounds) in Gammarus pulex collected from 15 sites covering 5 river catchments across Suffolk, UK. Quantitative method performance was acceptable for 67 compounds including pharmaceuticals, pesticides, illicit drugs and drugs of abuse. A total of 56 compounds were detectable and ranged from Collapse

Key Words

• Environmental risk assessment
• Exposome
• Pesticides
• Pharmaceuticals

Collapse

MESH Headings

• Animals
• Environmental Monitoring/methods
• Fresh Water
• Humans
• Illicit Drugs/chemistry
• Illicit Drugs/toxicity
• Invertebrates/chemistry
• Pesticides/chemistry
• Pesticides/toxicity
• Toxicity Tests
• Water Pollutants, Chemical/chemistry

Collapse

Grants

• MR/S019669/1 Medical Research Council

Collapse

Affiliation(s)

Thomas H Miller Department of Analytical, Environmental & Forensic Sciences, School of Population Health & Environmental Sciences, Faculty of Life Sciences and Medicine, King's College London, 150 Stamford Street, London SE1 9NH, UK.
Keng Tiong Ng Department of Analytical, Environmental & Forensic Sciences, School of Population Health & Environmental Sciences, Faculty of Life Sciences and Medicine, King's College London, 150 Stamford Street, London SE1 9NH, UK
Samuel T Bury St Olaves Grammer School, Goddington Lane, Orpington, BR6 9SH, UK
Sophie E Bury Department of Pyschology, Royal Holloway, University of London, Egham, Surrey TW20 0EX, UK
Nicolas R Bury School of Science, Technology and Engineering, University of Suffolk, James Hehir Building, University Avenue, Ipswich, Suffolk IP3 0FS, UK; Suffolk Sustainability Institute, University of Suffolk, Waterfront Building, Neptune Quay, Ipswich IP4 1QJ, UK
Leon P Barron Department of Analytical, Environmental & Forensic Sciences, School of Population Health & Environmental Sciences, Faculty of Life Sciences and Medicine, King's College London, 150 Stamford Street, London SE1 9NH, UK

Collapse

Enantiomer stability and combined toxicity of duloxetine and econazole on Daphnia magna using real concentrations determined by capillary electrophoresis

Enantiomer stability was investigated in this work for the first time for duloxetine and econazole in individual solutions and their mixtures under the standardized ecotoxicity test experimental conditions for Daphnia magna and abiotic conditions. Real (and not nominal) enantiomer concentrations were employed for calculations since their determination was achieved by Capillary Electrophoresis. Relevant differences were found in stability profiles for both drugs in any case. Toxicity was evaluated for the first time in this work for mixtures of duloxetine and econazole on Daphnia magna. Dose-effect parameters were calculated at different exposure times (24, 48, and 72 h) showing a significant inhibition of daphnids mobility when increasing the incubation time. Combination index values enabled to obtain the type and level of interaction of drugs with the organism. A strong synergism was observed at 48 h exposure time and any effect level, which demonstrated the high toxicity of the drug mixture compared with the individual drug solutions. These results were corroborated when evaluating the oxidative stress using fluorescence images.

Transformation of norfloxacin during the chlorination of marine culture water in the presence of iodide ions

The antibacterial agent norfloxacin (NOR) and sodium hypochlorite (NaClO), which are both widely used in marine culture, react with each other to form the halogenated disinfection byproducts (X-DBPs). The effects of the water characteristics and iodide concentration on the reaction kinetics were investigated. The results showed that the reaction rate of NOR with NaClO increases from 0.0586 min^-1 to 0.1075 min^-1 when the iodide concentration was changed from 0 μg^-1 to 50 μg^-1. This demonstrated the enhancement of NOR oxidation in the presence of iodide ions. Four novel iodinated DBPs (I-DBPs) were identified in the marine culture water. Iodine substitutions occurred at the C3 and C8 positions of NOR. The formation mechanisms of X-DBPs in the marine culture water were proposed based on the intermediate and final products. NOR may undergo a ring-opening reaction, a de-carbonyl reaction and substitution to form intermediates and finally generate the X-DBPs. Furthermore, the predicted logK_OW and logBCF values of the I-DBPs were higher than that of the Br-DBPs and Cl-DBPs. The AOX concentration in the synthetic water samples decreased in the following order: seawater (8.49 mg L^-1) > marine culture water (4.05 mg L^-1) > fresh water (1.89 mg L^-1). The amount of AOX also increased with the increase in iodide concentration. These results indicated that the I-DBPs were more toxic than their brominated and chlorinated analogues.

Prediction of bioconcentration factors in fish and invertebrates using machine learning

The application of machine learning has recently gained interest from ecotoxicological fields for its ability to model and predict chemical and/or biological processes, such as the prediction of bioconcentration. However, comparison of different models and the prediction of bioconcentration in invertebrates has not been previously evaluated. A comparison of 24 linear and machine learning models is presented herein for the prediction of bioconcentration in fish and important factors that influenced accumulation identified. R² and root mean square error (RMSE) for the test data (n = 110 cases) ranged from 0.23-0.73 and 0.34-1.20, respectively. Model performance was critically assessed with neural networks and tree-based learners showing the best performance. An optimised 4-layer multi-layer perceptron (14 descriptors) was selected for further testing. The model was applied for cross-species prediction of bioconcentration in a freshwater invertebrate, Gammarus pulex. The model for G. pulex showed good performance with R² of 0.99 and 0.93 for the verification and test data, respectively. Important molecular descriptors determined to influence bioconcentration were molecular mass (MW), octanol-water distribution coefficient (logD), topological polar surface area (TPSA) and number of nitrogen atoms (nN) among others. Modelling of hazard criteria such as PBT, showed potential to replace the need for animal testing. However, the use of machine learning models in the regulatory context has been minimal to date and is critically discussed herein. The movement away from experimental estimations of accumulation to in silico modelling would enable rapid prioritisation of contaminants that may pose a risk to environmental health and the food chain.

Spatio-temporal bioaccumulation and trophic transfer of ionizable pharmaceuticals in a semi-arid urban river influenced by snowmelt

Bioaccumulation of pharmaceuticals in aquatic organisms is increasingly reported in the peer-reviewed literature. However, seasonal instream dynamics including occurrence and bioaccumulation across trophic positions are rarely studied, particularly in semiarid streams with flows influenced by seasonal snowmelt and municipal effluent discharges. Thus, we selected East Canyon Creek in Park City, Utah, USA to examine spatio-temporal bioaccumulation of select ionizable pharmaceuticals across trophic positions using trophic magnification factors calculated at incremental distances (0.15, 1.4, 13 miles) downstream from a municipal effluent discharge during spring (May), Summer (August), and fall (October). Nine target analytes were detected in all species during all sampling events. Trophic dilution was consistently observed for amitriptyline, caffeine, diphenhydramine, diltiazem, fluoxetine, and sertraline, regardless of seasonal instream flows or distance from effluent discharge. Calculated TMFs ranged from 0.01-0.71 with negative slopes observed for all regressions of chemical residue in tissue and trophic position. We further presents the first empirical investigation of normalizing pharmaceutical concentrations to lipid, phospholipid or protein fractions using pair matched fish samples. Empirical results identify that normalization of ionizable pharmaceutical residues in aquatic tissues to neutral lipids, polar lipids, or the total protein fraction is inappropriate, though bioaccumulation studies examining influences of internal partitioning (e.g., plasma proteins) are needed.

A review of the pharmaceutical exposome in aquatic fauna

Pharmaceuticals have been considered 'contaminants of emerging concern' for more than 20 years. In that time, many laboratory studies have sought to identify hazard and assess risk in the aquatic environment, whilst field studies have searched for targeted candidates and occurrence trends using advanced analytical techniques. However, a lack of a systematic approach to the detection and quantification of pharmaceuticals has provided a fragmented literature of serendipitous approaches. Evaluation of the extent of the risk for the plethora of human and veterinary pharmaceuticals available requires the reliable measurement of trace levels of contaminants across different environmental compartments (water, sediment, biota - of which biota has been largely neglected). The focus on pharmaceutical concentrations in surface waters and other exposure media have therefore limited both the characterisation of the exposome in aquatic wildlife and the understanding of cause and effect relationships. Here, we compile the current analytical approaches and available occurrence and accumulation data in biota to review the current state of research in the field. Our analysis provides evidence in support of the 'Matthew Effect' and raises critical questions about the use of targeted analyte lists for biomonitoring. We provide six recommendations to stimulate and improve future research avenues.