A physiologically based toxicokinetic (PBTK) model for moderately hydrophobic organic chemicals in the European eel (Anguilla anguilla)

Fish Physiologically Based Toxicokinetic Modeling Approach for In Vitro-In Vivo and Cross-Species Extrapolation of Endocrine-Disrupting Chemicals in Risk Assessment

High-throughput in vitro assays combined with in vitro-in vivo extrapolation (IVIVE) leverage in vitro responses to predict the corresponding in vivo exposures and thresholds of concern. The integrated approach is also expected to offer the potential for efficient tools to provide estimates of chemical toxicity to various wildlife species instead of animal testing. However, developing fish physiologically based toxicokinetic (PBTK) models for IVIVE in ecological applications is challenging, especially for plausible estimation of an internal effective dose, such as fish equivalent concentration (FEC). Here, a fish PBTK model linked with the IVIVE approach was established, with parameter optimization of chemical unbound fraction, pH-dependent ionization and hepatic clearance, and integration of temperature effect and growth dilution. The fish PBTK-IVIVE approach provides not only a more precise estimation of tissue-specific concentrations but also a reasonable approximation of FEC targeting the estrogenic potency of endocrine-disrupting chemicals. Both predictions were compared with in vivo data and were accurate for most indissociable/dissociable chemicals. Furthermore, the model can help determine cross-species variability and sensitivity among the five fish species. Using the available IVIVE-derived FEC with target pathways is helpful to develop predicted no-effect concentration for chemicals with similar mode of action and support screening-level ecological risk assessment.

Physiologically based toxicokinetic modelling of Tri(2-chloroethyl) phosphate (TCEP) in mice accounting for multiple exposure routes

Exposure routes are important for health risk assessment of chemical risks. The application of physiologically based toxicokinetic (PBTK) models to predict concentrations in vivo can determine the effects of harmful substances and tissue accumulation on the premise of saving experimental costs. In this study, Tri(2-chloroethyl) phosphate (TCEP), an organophosphate ester (OPE), was used as an example to study the PBTK model of mice exposed to different exposure doses by multiple routes. Different routes of exposure (gavage and intradermal injection) can cause differences in the concentration of chemicals in the organs. TCEP that enters the body through the mouth is mainly concentrated in the gastrointestinal tract and liver. However, the concentrations of chemicals that enter the skin into the mice are higher in skin, rest of body, and blood. In addition, TCEP was absorbed and accumulated very rapidly in mice, within half an hour after a single exposure. We have successfully established a mouse PBTK model of the TCEP accounting for multiple exposure Routes and obtained a series of kinetic parameters. The model includes blood, liver, kidney, stomach, intestine, skin, and rest of body compartments. Oral and dermal exposure route was considered for PBTK model. The PBTK model established in this study has a good predictive ability. More than 70% of the predicted values deviated from the measured values by less than 5-fold. In addition, we extrapolated the model to humans. A human PBTK model is built. We performed a health risk assessment for world populations based on human PBTK model. The risk of TCEP in dust is greater through mouth than through skin. The risk of TCEP in food of Chinese population is greater than dust.

Bio-uptake, tissue distribution and metabolism of a neonicotinoid insecticide clothianidin in zebrafish

Neonicotinoids have been often detected in aquatic environment with high concentrations; however, little is known about their risk and fate to/in fish. This study systematically investigated the bio-uptake, tissue distribution and metabolism of neonicotinoids in zebrafish, taking clothianidin (CLO) as an example. The results revealed the uptake and elimination kinetics of CLO in whole fish and different tissues was very similar, and its bioconcentration factor (<1) indicates the low bioaccumulation potential in zebrafish. The highest accumulative tissues for CLO were found to be intestine and liver. Eight biotransformation products were identified in intestine and liver, and the metabolic pathways were found to be N-demethylation and nitro-reduction. The metabolic kinetics of two products (desmethyl clothianidin and clothianidin urea) revealed the metabolism of CLO mainly occurred in liver and intestine. This suggested that the hepatobiliary system played an important role in the metabolism and elimination of CLO. This study provides a comprehensive evaluation of the toxicokinetics of CLO in zebrafish, and these results can contribute to its ecological risk assessment.

A Novel Multispecies Toxicokinetic Modeling Approach in Support of Chemical Risk Assessment

Standardized laboratory tests with a limited number of model species are a key component of chemical risk assessments. These surrogate species cannot represent the entire diversity of native species, but there are practical and ethical objections against testing chemicals in a large variety of species. In previous research, we have developed a multispecies toxicokinetic model to extrapolate chemical bioconcentration across species by combining single-species physiologically based toxicokinetic (PBTK) models. This "top-down" approach was limited, however, by the availability of fully parameterized single-species models. Here, we present a "bottom-up" multispecies PBTK model based on available data from 69 freshwater fishes found in Canada. Monte Carlo-like simulations were performed using statistical distributions of model parameters derived from these data to predict steady-state bioconcentration factors (BCFs) for a set of well-studied chemicals. The distributions of predicted BCFs for 1,4-dichlorobenzene and dichlorodiphenyltrichloroethane largely overlapped those of empirical data, although a tendency existed toward overestimation of measured values. When expressed as means, predicted BCFs for 26 of 34 chemicals (82%) deviated by less than 10-fold from measured data, indicating an accuracy similar to that of previously published single-species models. This new model potentially enables more environmentally relevant predictions of bioconcentration in support of chemical risk assessments.

Combined sediment desorption and bioconcentration model to predict levels of dioxin-like chemicals in fish

Flooding and other sediment disturbances can lead to increases in sediment resuspension. In this context, it is of central importance to understand the kinetics of release from these sediments and the uptake of pollutants, such as polychlorinated biphenyls (PCBs) and polychlorinated dioxins and furans (PCDD/Fs), into aquatic organisms. In the present study, we parameterized a sediment desorption model based on experimentally determined rapidly-desorbing fractions of dioxin-like chemicals (DLCs). We coupled this desorption model with a physiologically-based toxicokinetic model for rainbow trout. This combined model was used to predict DLC concentrations in the muscle of exposed fish. The performance of this model was evaluated using a previously published dataset on DLC uptake from sediment suspensions during simulated re-suspension events. Predictions generally differed less than 10-fold from measured values, and the model showed a good global coefficient of determination (R²) of 0.95. The root mean squared error (RMSE) for PCBs was 0.31 log units and 0.53 log units for PCDD/Fs. The results of our study demonstrate that the prediction of bioconcentration and related risk to fish resulting from sediment resuspension can be accurately predicted using coupled desorption and toxicokinetic models.

Assessment of the quality of European silver eels and tentative approach to trace the origin of contaminants - A European overview

The European eel is critically endangered. Although the quality of silver eels is essential for their reproduction, little is known about the effects of multiple contaminants on the spawning migration and the European eel management plan does not take this into account. To address this knowledge gap, we sampled 482 silver eels from 12 catchments across Europe and developed methods to assess three aspects of eel quality: muscular lipid content (N = 169 eels), infection with Anguillicola crassus (N = 482), and contamination by persistent organic pollutants (POPs, N = 169) and trace elements (TEs, N = 75). We developed a standardized eel quality risks index (EQR) using these aspects for the subsample of 75 female eels. Among 169 eels, 33% seem to have enough muscular lipids content to reach the Sargasso Sea to reproduce. Among 482 silver eels, 93% were infected by A. crassus at least once during their lifetime. All contaminants were above the limit of quantification, except the 1,2-bis(2,4,6-tribromophenoxy)ethane (BTBPE), Ag and V. The contamination by POPs was heterogeneous between catchments while TEs were relatively homogeneous, suggesting a multi-scale adaptation of management plans. The EQR revealed that eels from Warwickshire were most impacted by brominated flame-retardants and agricultural contaminants, those from Scheldt were most impacted by agricultural and construction activities, PCBs, coal burning, and land use, while Frémur eels were best characterized by lower lipid contents and high parasitic and BTBPE levels. There was a positive correlation between EQR and a human footprint index highlighting the capacity of silver eels for biomonitoring human activities and the potential impact on the suitability of the aquatic environment for eel population health. EQR therefore represents a step forward in the standardization and mapping of eel quality risks, which will help identify priorities and strategies for restocking freshwater ecosystems.

Temperature effect on perfluorooctane sulfonate toxicokinetics in rainbow trout (Oncorhynchus mykiss): Exploration via a physiologically based toxicokinetic model

Salmonids are poikilotherms, which means that their internal temperature varies with that of water. Water temperature thus controls many of their lifecycle processes and physiological functions, which could influence the mechanisms of absorption, distribution, metabolism and excretion (ADME) of many substances, including perfluorinated alkyl acids (PFAAs). However, the processes governing the fate of PFAAs are still poorly understood in fish. Here we developed a physiologically-based toxicokinetic (PBTK) model for rainbow trout (Oncorhynchus mykiss) to study changes in physiological functions and PFAA ADME at different temperatures. The model was calibrated using experimental data from dietary exposure to perfluorooctane sulfonate at 7 °C and 19 °C. Predictions of PFOS concentrations were globally satisfactory at both temperatures, when accounting for the influence of temperature on growth, ventilation rate, cardiac output, clearances, and absorption rates. Accounting for the influence of temperature on tissue-plasma partition coefficients significantly improved predicted in-organ PFOS concentrations.

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

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

Physiologically based toxicokinetic and toxicodynamic (PBTK-TD) modelling of Cd and Pb exposure in adult zebrafish Danio rerio: Accumulation and toxicity

Accurately predicting the accumulation and toxicity of metals in organisms is a challenging work in ecotoxicology. Here, we developed and validated a physiologically based toxicokinetic and toxicodynamic (PBTK-TD) model for adult zebrafish exposed to Cd and Pb. The model included the gill, liver, intestine, gonad, carcass, and brain, which were linked by blood circulation in the PBTK process and by dynamic relationships between the target organ concentrations and mortality in the TD process. Results showed that the PBTK sub-model can accurately describe and predict the uptake, distribution and disposition kinetics of Cd and Pb in zebrafish. The exchange rates and the accumulation of the metals in the organs were significantly different. For Cd, the highest exchange rate was between blood and liver, and the greatest accumulation of Cd occurred in the liver. For Pb, the greatest accumulation occurred in the gill. The TD sub-model further indicated that metal concentrations in the gill may effectively act as more suitable indicator of Cd and Pb toxic effect than whole body or other organs. The proposed PBTK-TD model is helpful to understanding the fundamental processes by which zebrafish regulate the uptake and disposition of metal and to quantitatively predicting metal toxicity.

Generic physiologically-based toxicokinetic modelling for fish: Integration of environmental factors and species variability

One of the goals of environmental risk assessment is to protect the whole ecosystem from adverse effects resulting from exposure to chemicals. Many research efforts have aimed to improve the quantification of dose-response relationships through the integration of toxicokinetics. For this purpose, physiologically-based toxicokinetic (PBTK) models have been developed to estimate internal doses from external doses in a time-dependent manner. In this study, a generic PBTK model was developed and adapted for rainbow trout (Onchorhynchus mykiss), zebrafish (Danio rerio), fathead minnow (Pimephales promelas), and three-spined stickleback (Gasterosteus aculeatus). New mechanistic approaches were proposed for including the effects of growth and temperature in the model. Physiological parameters and their inter-individual variability were estimated based on the results of extensive literature searches or specific experimental data. The PBTK model was implemented for nine environmental contaminants (with log k_ow from -0.9 to 6.8) to predict whole-body concentrations and concentrations in various fish's organs. Sensitivity analyses were performed for a lipophilic and a hydrophilic compound to identify which parameters have most impact on the model's outputs. Model predictions were compared with experimental data according to dataset-specific exposure scenarios and were accurate: 50% of predictions were within a 3-fold factor for six out of nine chemicals and 75% of predictions were within a 3-fold factor for three of the most lipophilic compounds studied. Our model can be used to assess the influence of physiological and environmental factors on the toxicokinetics of chemicals and provide guidance for assessing the effect of those critical factors in environmental risk assessment.

Biotransformation of Benzo[ a]pyrene by Three Rainbow Trout ( Onchorhynchus mykiss) Cell Lines and Extrapolation To Derive a Fish Bioconcentration Factor

Permanent fish cell lines constitute a promising complement or substitute for fish in the environmental risk assessment of chemicals. We demonstrate the potential of a set of cell lines originating from rainbow trout ( Oncorhynchus mykiss) to aid in the prediction of chemical bioaccumulation in fish, using benzo[ a]pyrene (BaP) as a model chemical. We selected three cell lines from different tissues to more fully account for whole-body biotransformation in vivo: the RTL-W1 cell line, representing the liver as major site of biotransformation, and the RTgill-W1 (gill) and RTgutGC (intestine) cell lines, as important environment-organism interfaces, which likely influence chemical uptake. All three cell lines were found to effectively biotransform BaP. However, rates of in vitro clearance differed, with the RTL-W1 cell line being most efficient, followed by RTgutGC. Co-exposures with α-naphthoflavone as potent inhibitor of biotransformation, assessment of CYP1A catalytic activity, and the progression of cellular toxicity upon prolonged BaP exposure revealed that BaP is handled differently in the RTgill-W1 compared to the other two cell lines. Application of the cell-line-derived in vitro clearance rates into a physiology-based toxicokinetic model predicted a BaP bioconcentration factor (BCF) of 909-1057 compared to 920 reported for rainbow trout in vivo.

Toxicokinetic models and related tools in environmental risk assessment of chemicals

Environmental risk assessment of chemicals for the protection of ecosystems integrity is a key regulatory and scientific research field which is undergoing constant development in modelling approaches and harmonisation with human risk assessment. This review focuses on state-of-the-art toxicokinetic tools and models that have been applied to terrestrial and aquatic species relevant to environmental risk assessment of chemicals. Both empirical and mechanistic toxicokinetic models are discussed using the results of extensive literature searches together with tools and software for their calibration and an overview of applications in environmental risk assessment. These include simple tools such as one-compartment models, multi-compartment models to physiologically-based toxicokinetic (PBTK) models, mostly available for aquatic species such as fish species and a number of chemical classes including plant protection products, metals, persistent organic pollutants, nanoparticles. Data gaps and further research needs are highlighted.

Evidence for High Concentrations and Maternal Transfer of Substituted Diphenylamines in European Eels Analyzed by Two-Dimensional Gas Chromatography-Time-of-Flight Mass Spectrometry and Gas Chromatography-Fourier Transform Ion Cyclotron Resonance Mass Spectrometry

Chemical pollution is hypothesized to be one of the factors driving the strong decline of the critically endangered European eel population. Specifically, the impact of contaminants on the quality of spawning eels and subsequent embryo survival and development has been discussed as crucial investigation point. However, so far, only very limited information on potential negative effects of contaminants on the reproduction of eels is available. Through the combination of nontargeted ultrahigh-resolution mass spectrometry and multidimensional gas chromatography, combined with more-conventional targeted analytical approaches and multimedia mass-balance modeling, compounds of particular relevance, and their maternal transfer in artificially matured European eels from the German river Ems have been identified. Substituted diphenylamines were, unexpectedly, found to be the primary organic contaminants in the eel samples, with concentrations in the μg g^-1 wet weight range. Furthermore, it could be shown that these contaminants, as well as polychlorinated biphenyls (PCBs), organochlorine pesticides, and polyaromatic hydrocarbons (PAHs), are not merely stored in lipid rich tissue of eels but maternally transferred into gonads and eggs. The results of this study provide unique information on both the fate and behavior of substituted diphenylamines in the environment as well as their relevance as contaminants in European eels.

Fipronil and two of its transformation products in water and European eel from the river Elbe

Fipronil is an insecticide which, based on its mode of action, is intended to be predominantly toxic towards insects. Fipronil bioaccumulates and some of its transformation products were reported to be similar or even more stable in the environment and to show an enhanced toxicity against non-target organisms compared to the parent compound. The current study investigated the occurrence of Fipronil and two of its transformation products, Fipronil-desulfinyl and Fipronil-sulfone, in water as well as muscle and liver samples of eels from the river Elbe (Germany). In water samples total concentrations of FIP, FIP-d and FIP-s ranged between 0.5-1.6ngL(-1) with FIP being the main component in all water samples followed by FIP-s and FIP-d. In contrast, FIP-s was the main component in muscle and liver tissues of eels with concentrations of 4.05±3.73ngg(-1) ww and 19.91±9.96ngg(-1) ww, respectively. Using a physiologically based toxicokinetic (PBTK) model for moderately hydrophobic organic chemicals, the different distributions of FIP, FIP-d and FIP-s in water and related tissue samples could be attributed to metabolic processes of eels. The measured concentrations in water of all analytes and their fractional distribution did not reflect the assumed seasonal application of FIP and it seems that the water was constantly contaminated with FIP, FIP-d and FIP-s.

Accumulation of neurotoxic organochlorines and trace elements in brain of female European eel (Anguilla anguilla)

Xenobiotics such as organochlorine compounds (OCs) and metals have been suggested to play a significant role in the collapse of European eel stocks in the last decades. Several of these pollutants could affect functioning of the nervous system. Still, no information is so far available on levels of potentially neurotoxic pollutants in eel brain. In present study, carried out on female eels caught in Belgian rivers and canals, we analyzed brain levels of potentially-neurotoxic trace elements (Ag, Al, As, Cd, Co, Cr, Cu, Fe, Hg, MeHg, Mn, Ni, Pb, Sn, Sb, Zn) and OCs (Polychlorinated biphenyls, PCBs; Hexachlorocyclohexanes, HCHs; Dichlorodiphenyltrichloroethane and its metabolites, DDTs). Data were compared to levels in liver and muscle tissues. Eel brain contained very high amounts of OCs, superior to those found in the two other tissues. Interestingly, the relative abundance of PCB congeners markedly differed between tissues. In brain, a predominance of low chlorinated PCBs was noted, whereas highly chlorinated congeners prevailed in muscle and liver. HCHs were particularly abundant in brain, which contains the highest amounts of β-HCH and ϒ-HCH. p,p'-DDTs concentration was similar between brain and muscle (i.e., about twice that of liver). A higher proportion of p,p'-DDT was noticed in brain. Except for Cr and inorganic Hg, all potentially neurotoxic metals accumulated in brain to levels equal to or lower than hepatic levels. Altogether, results indicate that eel brain is an important target for organic and, to a lesser extent, for inorganic neurotoxic pollutants.

Dioxin in the Elbe river basin: policy and science under the water framework directive 2000-2015 and toward 2021

A critical review of the last 25 years of dioxin policy in the Elbe river catchment is presented along seven main theses of the River Basin Community (RBC)-Elbe background document "Pollutants" for the Management Plan 2016-2021. In this period, polychlorinated dibenzodioxins/-furans (PCDD/Fs) and dioxin-like polychlorinated biphenyls (dl-PCBs) will play a major role: (i) as new priority substances for which environmental quality standards (EQSs) need to be derived (Directive 2013/39/EC); (ii) in the search for innovative solutions in sediment remediation (i.e., respecting the influence of mechanical processes; Flood Risk Directive 2007/60/EC); and (iii) as indicators at the land-sea interface (Marine Strategy Framework Directive 2008/56/EC). In the Elbe river catchment, aspects of policy and science are closely connected, which became particularly obvious in a classic example of dioxin hot spot contamination, the case of the Spittelwasser creek. Here, the "source-first principle" of the first cycle of the European Water Framework Directive (WFD) had to be confirmed in a controversy on the dioxin hot spots with Saxony-Anhalt's Agency for Contaminated Sites (LAF). At the Spittelwasser site, the move from "inside the creek" to "along the river banks" goes parallel to a general paradigm shift in retrospective risk assessment frameworks and remediation techniques for organic chemicals (Ortega-Calvo et al. 2015). With respect to dioxin, large-scale stabilization applying activated carbon additions is particularly promising. Another important aspect is the assessment of the ecotoxicology of dioxins and dl- PCBs in context of sediment mobility and flood risk assessment, which has been studied in the project framework FloodSearch. Currently, the quality goals of the WFD to reach a "good chemical status" are not met in many catchment areas because substances such as mercury do and others probably will (PCDD/Fs and dl-PCB) exceed biota-EQS values catchment area-wide. So far, relating biota-EQS values to sediment-EQSs is not possible. To overcome these limitations, the DioRAMA project was initiated, which has led to improved approaches for the assessment of dioxin-contaminated sediment using in vitro bioassays and to a robust dataset on the interrelation between dioxins and dioxin-like compounds in sediments and biota.