A physiologically based toxicokinetic model for the zebrafish Danio rerio

Contrasting the distribution kinetics of microplastics and nanoplastics in medaka following exposure and depuration

Direct ingestion of micro/nanoplastics (MNPs) results in significant accumulation in gastrointestinal (GI) tract of fish. The breathing process of fish makes MNPs easily retained in their gills. However, the uptake of MNPs in other fish organs remains largely unknown, let alone their kinetic processes. Herein, microplastics (MPs) and nanoplastics (NPs) in vivo imaging and precise quantification in various tissues (GI tract, gill, liver, brain, eye, and skin) of seawater (SW)- and freshwater (FW)- acclimated medaka Oryzias melastigma were achieved at an environmentally relevant concentration. Subsequently, the distribution kinetics of MNPs was investigated over a 96-h uptake and 48-h depuration period. MNPs were quickly and mostly captured in GI tract and gill of O. melastigma, and then transferred to liver and brain likely via blood circulation. Such transport was more efficient for NPs as compared to MPs, as evidenced by the consistently higher bioconcentration factors in both SW and FW conditions. The detection of MNPs in eye and skin of O. melastigma was more of an adsorption process, although the specific mechanisms of adsorption and absorption process can hardly be clearly differentiated. This study presented distribution kinetics of MNPs in O. melastigma and highlighted their possible transportation among tissues.

Imazalil and prochloraz toxicokinetics in fish probed by a physiologically based kinetic (PBK) model

Azole fungicides are highly suspected endocrine disruptors (EDs) and are frequently detected in surface water. Among them, there are prochloraz (PCZ), a commonly used  molecule for ED studies, and imazalil (IMZ), a highly suspected ED. Little is known about their toxicokinetic (TK) behavior in fish. Hence, research suggested that an improved risk assessment could be achieved by gaining insight into their TK behavior. The aim of this study is to understand and model the TK of both substances in different fish species, irrespective of the scheme of exposure. TK data from the literature were retrieved including different modes of exposure (per os and waterborne). In addition, two experiments on zebrafish exposed to either IMZ or PCZ were performed to address the lack of in vivo TK data. A physiologically based kinetic (PBK) model applied to IMZ and PCZ was developed, capable of modeling different exposure scenarios. The parameters of the PBK model were simultaneously calibrated on datasets reporting internal concentration in several organs in three fish species (original and literature datasets) by Bayesian methods (Monte Carlo Markov Chain). Model predictions were then compared to other experimental data (i.e., excluded from the calibration step) to assess the predictive performance of the model. The results strongly suggest that PCZ and IMZ are actively transported across the gills, resulting in a small fraction being effectively absorbed by the fish. The model's results also confirm that both molecules are extensively metabolized by the liver into mainly glucuronate conjugates. Overall, the model performances were satisfying, predicting internal concentrations in several key organs. On average, 90% of experimental data were predicted within a two-fold range. The PBK model allows the understanding of IMZ and PCZ kinetics profiles by accurately predicting internal concentrations in three different fish species regardless of the exposure scenario. This enables a proper understanding of the mechanism of action of EDs at the molecular initiating event (MIE) by predicting bioaccumulation in target organs, thus linking this MIE to a possible adverse outcome.

Exploring BPA alternatives - Environmental levels and toxicity review

Bisphenol A alternatives are manufactured as potentially less harmful substitutes of bisphenol A (BPA) that offer similar functionality. These alternatives are already in the market, entering the environment and thus raising ecological concerns. However, it can be expected that levels of BPA alternatives will dominate in the future, they are limited information on their environmental safety. The EU PARC project highlights BPA alternatives as priority chemicals and consolidates information on BPA alternatives, with a focus on environmental relevance and on the identification of the research gaps. The review highlighted aspects and future perspectives. In brief, an extension of environmental monitoring is crucial, extending it to cover BPA alternatives to track their levels and facilitate the timely implementation of mitigation measures. The biological activity has been studied for BPA alternatives, but in a non-systematic way and prioritized a limited number of chemicals. For several BPA alternatives, the data has already provided substantial evidence regarding their potential harm to the environment. We stress the importance of conducting more comprehensive assessments that go beyond the traditional reproductive studies and focus on overlooked relevant endpoints. Future research should also consider mixture effects, realistic environmental concentrations, and the long-term consequences on biota and ecosystems.

Predicting PFAS fate in fish: Assessing the roles of dietary, respiratory, and dermal uptake in bioaccumulation modeling

An increasing number of per- and polyfluoroalkyl substances (PFAS) exposed to the environment may pose a threat to organisms and human beings. However, there is a lack of simulations comprehensively addressing and comparing the bioaccumulation of PFAS across all three major exposure routes (oral, inhalation, and dermal), especially for dermal uptake. In this study, we proposed a physiologically based kinetic (PBK) model for PFAS, aiming to predict bioaccumulation factors (BAF) in fish by considering these diverse exposure routes. 15 PFAS were used for model validation, and 11 PFAS from Taihu Lake were used for exposure contribution modeling. Approximately 64% of estimations fell within 10-fold model bias from measurements in Taihu Lake, underscoring the potential efficacy of the developed PBK model in predicting BAFs for fish. The dermal route emerges as a contributor to short-chain PFAS exposure. For example, it ranged widely from 46% to 75% (mean) for all modeling short-chain PFAS (C6-C7) in Taihu Lake. It indicated the criticality of considering dermal exposure for PFAS in fish, highlighting a gap in field studies to unravel cutaneous intake mechanisms and contributions. For longer carbon chains of PFAS (C8-C12), dermal exposure accounted for 2%-27% for all species of aquatic organisms. The fish's lipid fraction and water content played a significant role in the contribution of PFAS intake through cutaneous exposure and inhalation. Kow had a significant positive correlation with skin intake rate (p < 0.05) and gill intake rate (p < 0.001), while having a significant negative correlation with skin intake (p < 0.05) and skin intake contribution (p < 0.001). Based on the proposed modeling approach, we have introduced a simulation spreadsheet for projecting PFAS BAFs in fish tissues, hopefully broadening the predictive operational tool for a variety of chemical species.

Physiologically based kinetic (PBK) modeling as a new approach methodology (NAM) for predicting systemic levels of gut microbial metabolites

There is a clear need to develop new approach methodologies (NAMs) that combine in vitro and in silico testing to reduce and replace animal use in chemical risk assessment. Physiologically based kinetic (PBK) models are gaining popularity as NAMs in toxico/pharmacokinetics, but their coverage of complex metabolic pathways occurring in the gut are incomplete. Chemical modification of xenobiotics by the gut microbiome plays a critical role in the host response, for example, by prolonging exposure to harmful metabolites, but there is not a comprehensive approach to quantify this impact on human health. There are examples of PBK models that have implemented gut microbial biotransformation of xenobiotics with the gut as a dedicated metabolic compartment. However, the integration of microbial metabolism and parameterization of PBK models is not standardized and has only been applied to a few chemical transformations. A challenge in this area is the measurement of microbial metabolic kinetics, for which different fermentation approaches are used. Without a standardized method to measure gut microbial metabolism ex vivo/in vitro, the kinetic constants obtained will lead to conflicting conclusions drawn from model predictions. Nevertheless, there are specific cases where PBK models accurately predict systemic concentrations of gut microbial metabolites, offering potential solutions to the challenges outlined above. This review focuses on models that integrate gut microbial bioconversions and use ex vivo/in vitro methods to quantify metabolic constants that accurately represent in vivo conditions.

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.

Non-negligible vector effect of micro(nano)plastics on tris(1,3-dichloro-2-propyl) phosphate in zebrafish quantified by toxicokinetic model

Micro(nano)plastics (MNPs) inevitably interact with coexisting contaminants and can act as vectors to affect their fate in organisms. However, the quantitative contribution of MNPs in the in vivo bioaccumulation and distribution of their coexisting contaminants remains unclear. Here, by selecting tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) as the typical coexisting contaminant, we quantified the contribution of MNPs to bioaccumulation and distribution of TDCIPP with toxicokinetic models. Results indicated that MNPs differentially facilitated TDCIPP bioaccumulation and distribution, and NPs slowed down TDCIPP depuration more significantly than MPs. Model analysis further revealed increasing contributions of MNPs to whole-fish TDCIPP bioaccumulation over time, with NPs (33-42%) contributing more than MPs (12-32%) at 48 h exposure. NPs contributed more than MPs to TDCIPP distribution in the liver (13-19% for MPs; 36-52% for NPs) and carcass (24-45% for MPs; 57-71% for NPs). The size-dependent vector effect might be attributed to the fact that MNPs promote contaminant transfer by damaging biofilm structure and increasing tissue membrane permeability, with NPs exerting stronger effects. This work demonstrated the effectiveness of using modeling tools to understand the relative importance of MNPs as contaminant vectors in the TK process and highlighted the higher contaminant transfer potential of NPs under combined exposure scenarios.

Tissue-specific accumulation, depuration and histopathological effects of 3,6-dichlorocarbazole and 2,7-dibromocarbazole in adult zebrafish (Danio rerio)

Although polyhalogenated carbazoles have been detected with increasing frequency in aquatic ecosystems, their bioaccumulation in fish and corresponding pathological effects related to bioaccumulation are still unclear. Here, we investigated the tissue-specific accumulation, depuration, and histopathological effects of two typical PHCZs, 3,6-dichlorocarbazole (36-CCZ) and 2,7-dibromocarbazole (27-BCZ), in adult zebrafish at three levels (0, 0.15 μg/L (5 × environmentally relevant level), and 50 μg/L (1/10 LC50). The lowest concentrations of 36-CCZ (1.2 μg/g ww) and 27-BCZ (1.4 μg/g ww) were observed in muscle, and the greatest concentrations of 36-CCZ (3.6 μg/g ww) and 27-BCZ (4 μg/g ww) were detected in intestine among the tested tissues. BCFww of 36-CCZ and 27-BCZ in zebrafish ranged from 172.9 (muscle) to 606.6 (intestine) and 285.2 (muscle) to 987.5 (intestine), respectively, indicating that both 36-CCZ and 27-BCZ have high potential of bioaccumulation in aquatic system. The 0.15 μg/L level of 36-CCZ or 27-BCZ caused lipid accumulation in liver, while 50 μg/L of 36-CCZ or 27-BCZ induced liver lesions such as fibrous septa, cytolysis, and nuclear dissolution. Brain damage such as multinucleated cells and nuclear solidification were only observed at 50 μg/L of 27-BCZ. This study provided valuable information in assessing the health and ecological risks of 36-CCZ and 27-BCZ.

Diel Pattern of Microplastic Residues in Zebrafish

As one of the emerging pollutants, microplastics (MPs) can be taken up by aquatic organisms through ingestion. However, little is known about the uptake pattern in organisms over time and the associated mechanisms of retention patterns. The present study aims to elucidate these patterns in fish, their relationship with light/dark conditions, and examine the uptake kinetic process of small-sized plastic pollutants, especially during the long-neglected dark period. Zebrafish were sampled every 2 h during the light and dark periods after exposure to an environmentally relevant concentration (100 items/L) of MPs. The results demonstrated that MP residues in zebrafish decreased during the dark period rather than increased over time. The MP retention rhythm and the swimming behavior of exposed zebrafish displayed a statistically significant light/dark variation. Moreover, a very strong and statistically significant positive correlation was found between the swimming speed of zebrafish and the number of MP residues in the gastrointestinal tracts of zebrafish. These results clearly demonstrate that fibrous MP residues in the fish have a discernible diel pattern. This work improves the understanding of the dynamic residual process of MPs in organisms and calls for further in-depth circadian toxicokinetic studies to better suit particle pollutants.

Tissue-specific accumulation, bioaccumulation, and depuration of liquid crystal monomers (LCMs) in adult zebrafish (Danio rerio)

This study investigated the properties of bioaccumulation, tissue-specific accumulation, and depuration of liquid crystal monomers (LCMs) in adult zebrafish (Danio rerio) exposed to a mixture of 39 LCMs for 19 days followed by depuration for 12 days. Tissue-specific accumulation of LCMs was examined at the time point of day 19, and we observed that the distribution of LCMs varied among different tissues with the following order of Σ₃₉LCM concentrations, the intestine > brain > gill > liver > muscle. We also observed that the bioaccumulation potential of LCMs varied among LCM groups with different functional groups, and LCMs with the cyan group were prone to accumulate in zebrafish. Among the 39 target LCMs, nine exhibited bioconcentration factors (BCFs) >1000, indicating their great bioaccumulation potential in aquatic environments. The experimental BCFs of 22 LCMs in the present study were lower than the theoretical values predicted by the Estimation Programs Interface (EPI) Suite software developed by U.S. Environmental Protection Agency (USEPA), suggesting that their bioaccumulation potential might be overestimated by theoretical estimation techniques. Another interesting finding was the significant positive correlation relationship in both sexes of zebrafish (p < 0.01, r² = 0.66 for male; p < 0.01, r² = 0.41 for female) between logBCF_ww and logK_ow values of LCMs. Overall, this study provides fundamental information regarding the bioaccumulation potentials of LCMs, which could be helpful for further investigating the health risks of LCMs in aquatic environments.

A PBPK model to evaluate zebrafish eleutheroembryos' actual exposure: bisphenol A and analogs' (AF, F, and S) case studies

The zebrafish eleutheroembryo model is increasingly used to assess the toxicity and developmental adverse effects of xenobiotics. However, the actual exposure is seldom measured (poorly accessible), while a predictive model could estimate these concentrations. The predictions with a new eleutheroembryo physiologically based pharmacokinetic (PBPK) model have been evaluated using datasets obtained from literature data for several bisphenols. The model simulated the toxicokinetics of bisphenols A (BPA), AF, F, and S through the eleutheroembryo tissues while considering the body and organ growth. We further improved the predictions by adding dynamic flows through the embryo and/or its chorion, impact of experimental temperature, metabolic clearance, and saturation of the absorption by Bayesian calibration. The model structure was determined using the BPA dataset and generalized to the other bisphenols. This model revealed the central role of the chorion in the compound uptake in the first 48 h post-fertilization. The predictions for the BPA substitutes estimated by our PBPK model were compared to available toxicokinetics data for zebrafish embryos, and 63% and 88% of them were within a twofold and fivefold error intervals of the corresponding experimental values, respectively. This model provides a tool to design new eleutheroembryo assays and evaluate the actual exposure.

Enhanced prediction of internal concentrations of phenolic endocrine disrupting chemicals and their metabolites in fish by a physiologically based toxicokinetic incorporating metabolism (PBTK-MT) model

Bisphenol A (BPA), 4-nonylphenol (4-NP), and triclosan (TCS) are phenolic endocrine disrupting chemicals (EDCs), which are widely detected in aquatic environments and further bioaccumulated and metabolized in fish. Physiologically based toxicokinetic (PBTK) models have been used to describe the absorption, distribution, metabolism, and excretion (ADME) of parent compounds in fish, whereas the metabolites are less explored. In this study, a PBTK incorporating metabolism (PBTK-MT) model for BPA, 4-NP, and TCS was established to enhance the performance of the traditional PBTK model. The PBTK-MT model comprised 16 compartments, showing great accuracy in predicting the internal concentrations of three compounds and their glucuronidated and sulfated conjugates in fish. The impact of typical hepatic metabolism on the PBTK-MT model was successfully resolved by optimizing the mechanism for deriving the partition coefficients between the blood and liver. The PBTK-MT model exhibited a potential data gap-filling capacity for unknown parameters through a backward extrapolation approach of parameters. Model sensitivity analysis suggested that only five parameters were sensitive in at least two PBTK-MT models, while most parameters were insensitive. The PBTK-MT model will contribute to a well understanding of the environmental behavior and risks of pollutants in aquatic biota.

Toxicokinetic Modeling of Per- and Polyfluoroalkyl Substance Concentrations within Developing Zebrafish (Danio rerio) Populations

Per- and polyfluoroalkyl substances (PFAS) are pervasive environmental contaminants, and their relative stability and high bioaccumulation potential create a challenging risk assessment problem. Zebrafish (Danio rerio) data, in principle, can be synthesized within a quantitative adverse outcome pathway (qAOP) framework to link molecular activity with individual or population level hazards. However, even as qAOP models are still in their infancy, there is a need to link internal dose and toxicity endpoints in a more rigorous way to further not only qAOP models but adverse outcome pathway frameworks in general. We address this problem by suggesting refinements to the current state of toxicokinetic modeling for the early development zebrafish exposed to PFAS up to 120 h post-fertilization. Our approach describes two key physiological transformation phenomena of the developing zebrafish: dynamic volume of an individual and dynamic hatching of a population. We then explore two different modeling strategies to describe the mass transfer, with one strategy relying on classical kinetic rates and the other incorporating mechanisms of membrane transport and adsorption/binding potential. Moving forward, we discuss the challenges of extending this model in both timeframe and chemical class, in conjunction with providing a conceptual framework for its integration with ongoing qAOP modeling efforts.

Environmental health and risk assessment metrics with special mention to biotransfer, bioaccumulation and biomagnification of environmental pollutants

The environment pollutants, which are landed up in environment because of human activities like urbanization, mining and industrializations, affects human health, plants and animals. The living organisms present in environment are constantly affected by the toxic pollutants through direct contact or bioaccumulation of chemicals from the environment. The toxic and hazardous pollutants are easily transferred to different environmental matrices like land, air and water bodies such as surface and ground waters. This comprehensive review deeply discusses the routes and causes of different environmental pollutants along with their toxicity, impact, occurrences and fate in the environment. Environment health and risk assessment tools that are used to evaluate the harmfulness, exposure of living organisms to pollutants and the amount of pollutant accumulated are explained with help of bio-kinetic models. Biotransfer, toxicity factor, biomagnification and bioaccumulation of different pollutants in the air, water and marine ecosystems are critically addressed. Thus, the presented survey would be collection of correlations those addresses the factors involved in assessing the environmental health and risk impacts of distinct environmental pollutants.

Physiologically Based Toxicokinetic Modeling of Bisphenols in Zebrafish (Danio rerio) Accounting for Variations in Metabolic Rates, Brain Distribution, and Liver Accumulation

Bisphenol A (BPA) is an industrial chemical, which has raised human health and environmental concerns due to its endocrine-disrupting properties. BPA analogues are less well-studied despite their wide use in consumer products. These analogues have been detected in water and aquatic organisms around the world, with some analogues showing toxic effects in various species including fish. Here, we present novel organ-specific time-course distribution data of bisphenol Z (BPZ) in female zebrafish (Danio rerio), including concentrations in the ovaries, liver, and brain, a rarely sampled organ with high toxicological relevance. Furthermore, fish-specific in vitro biotransformation rates were determined for 11 selected bisphenols. A physiologically based toxicokinetic (PBTK) model was adapted for four of these bisphenols, which was able to predict levels in the gonads, liver, and brain as well as the whole body within a 2-5-fold error with respect to experimental data, covering several important target organs of toxicity. In particular, predicted liver concentrations improved compared to currently available PBTK models. Predicted data indicate that studied bisphenols mainly distribute to the carcass and gonads and less to the brain. Our model provides a tool to increase our understanding on the distribution and kinetics of a group of emerging pollutants.

The toxicokinetics of bisphenol A and its metabolites in fish elucidated by a PBTK model

Bisphenol A (BPA) is a chemical of major concern due to its endocrine disrupting function, high production volume, and persistence in the aquatic environment. Consequently, organisms such as fish are subject to chronic exposure to BPA. However, physiologically-based toxicokinetic (PBTK) models, which are valuable tools to improve the understanding of a chemical's fate in an organism, have never been specifically adapted to model BPA toxicokinetics (TK) in fish. In our work, an existing PBTK developed for four different fish species was modified to model BPA ADME processes (absorption, distribution, metabolization and excretion). The metabolization of BPA into BPA-monoglucuronide (BPA gluc) and BPA-monosulfate (BPA sulf) and their TK in various organs was taking into account in the model. Experiments were performed to generate BPA TK data in a model species commonly used in ecotoxicology, the stickleback. The model structure had to include two sites of metabolization to simulate BPA TK accurately in stickleback organs. Thus, the fish liver may not be the only site of the metabolization of BPA: plasma or gills could also play a role in BPA metabolization. The PBTK model predictive performance evaluated on literature data in zebrafish and rainbow trout concurs with this conclusion. Finally, a calibration mixing data from the three species was compared to the calibration on stickleback data only.

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).

A Generalized Physiologically Based Kinetic Model for Fish for Environmental Risk Assessment of Pharmaceuticals

An increasing number of pharmaceuticals found in the environment potentially impose adverse effects on organisms such as fish. Physiologically based kinetic (PBK) models are essential risk assessment tools, allowing a mechanistic approach to understanding chemical effects within organisms. However, fish PBK models have been restricted to a few species, limiting the overall applicability given the countless species. Moreover, many pharmaceuticals are ionizable, and fish PBK models accounting for ionization are rare. Here, we developed a generalized PBK model, estimating required parameters as functions of fish and chemical properties. We assessed the model performance for five pharmaceuticals (covering neutral and ionic structures). With biotransformation half-lives (HLs) from EPI Suite, 73 and 41% of the time-course estimations were within a 10-fold and a 3-fold difference from measurements, respectively. The performance improved using experimental biotransformation HLs (87 and 59%, respectively). Estimations for ionizable substances were more accurate than any of the existing species-specific PBK models. The present study is the first to develop a generalized fish PBK model focusing on mechanism-based parameterization and explicitly accounting for ionization. Our generalized model facilitates its application across chemicals and species, improving efficiency for environmental risk assessment and supporting an animal-free toxicity testing paradigm.

Effect of rare earth element lanthanum on lipid deposition and Wnt10b signaling in the liver of male zebrafish

This paper investigates the effect of lanthanum (La) on lipid deposition and Wnt10b signaling in the liver of male zebrafish with exposure of 0, 10, 20, and 30 μmol/L La. It suggests that La can be accumulated in liver, and its treatments decrease the activities and gene expression of enzymes related to fatty acid synthesis. The levels of total cholesterol (TC), triglyceride (TG), and nonesterified fatty acids (NEFA) as well as the size of lipid droplets are decreased by La treatments. Moreover, La treatments affect the composition of fatty acids and the content of nutrient elements. Meanwhile, they also induce the gene expression of wnt10b, β-catenin, pparα, and pparγ, but inhibit gsk-3β gene expression in liver. Further study on the result of wnt10b gene interference shows that Wnt10b/β-catenin signaling plays a crucial role in the regulatory process of hepatic lipid deposition. Taken together, our observations suggest that La accumulation affects lipid deposition in the liver of male zebrafish, and Wnt10b signaling pathway may be involved in this process.

Tissue-Specific Accumulation, Biotransformation, and Physiologically Based Toxicokinetic Modeling of Benzotriazole Ultraviolet Stabilizers in Zebrafish (Danio rerio)

Benzotriazole ultraviolet stabilizers (BUVSs) are high-production-volume chemicals with ubiquitous occurrence in the aquatic environment. However, little is known about their bioconcentration and biotransformation, and physiologically based toxicokinetic (PBTK) models for BUVSs are lacking. This study selected six BUVSs for which experiments were performed with zebrafish (Danio rerio) exposed to two different levels (0.5 and 10 μg·L^-1). Higher kinetic bioconcentration factors (BCFs) were observed at the lower exposure level with environmental relevance, with BCF of 3.33 × 10³ L·kg^-1 for 2-(2-hydroxy-3,5-di-tert-butylphenyl)-5-chlorobenzotriazole (UV-327). This phenomenon was interpreted by a nonlinear adsorption mechanism, where binding with specific protein sites contributes to bioconcentration. Muscle exhibited the lowest accumulation, in which depuration half-life of UV-327 was 19.5 d. In kidney, muscle, ovary, gill, and skin, logBCF increased with increase in log K_OW of the BUVSs until log K_OW was ca. 6.5, above which logBCF decreased. However, the trend was not observed in the liver and intestine. Six biotransformation products were identified and mainly accumulated in the liver and intestine. Considering the nonlinear adsorption mechanism in the PBTK model, the prediction accuracy of the model was improved, highlighting the binding of xenobiotics with specific protein sites in assessing the bioconcentration of chemicals for their 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.

Modeling acetylcholine esterase inhibition resulting from exposure to a mixture of atrazine and chlorpyrifos using a physiologically-based kinetic model in fish

Aquatic organisms are exposed to mixtures of chemicals that may interact. Mixtures of atrazine (ATR) and chlorpyrifos (CPF) may elicit synergic effects on the permanent inhibition of acetylcholinesterase (AChE) in certain aquatic organisms, causing severe damage. Mechanistic mathematical models of toxicokinetics and toxicodynamics (TD) may be used to better characterize and understand the interactions of these two chemicals. In this study, a previously published generic physiologically-based toxicokinetic (PBTK) model for fish was adapted to ATR and CPF. A sub-model of the kinetics of one of the main metabolites of CPF, chlorpyrifos-oxon (CPF-oxon), was included, as well as a TD model. Inhibition of two esterases, AChE and carboxylesterase, by ATR, CPF and CPF-oxon, was modeled using TD modeling of quantities of total and inactive esterases. Specific attention was given to the parameterization and calibration of the model to accurately predict the concentration and effects observed in the fish using Bayesian inference and published data from fathead minnow (Pimephales promelas), zebrafish (Danio rerio) and common carp (Cyprinus carpio L.). A PBTK-TD for mixtures was used to predict dose-response relationships for comparison with available adult fish data. Synergistic effects of a joint exposure to ATR and CPF could not be demonstrated in adult fish.

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.

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.

Development of a generic zebrafish embryo PBPK model and application to the developmental toxicity assessment of valproic acid analogs

In order to better explain, predict, or extrapolate to humans the developmental toxicity effects of chemicals to zebrafish (Danio rerio) embryos, we developed a physiologically-based pharmacokinetic (PBPK) model designed to predict organ concentrations of neutral or ionizable chemicals, up to 120 h post-fertilization. Chemicals' distribution is modeled in the cells, lysosomes, and mitochondria of ten organs of the embryo. The model's partition coefficients are calculated with sub-models using physicochemical properties of the chemicals of interest. The model accounts for organ growth and changes in metabolic clearance with time. We compared ab initio model predictions to data obtained on culture medium and embryo concentrations of valproic acid (VPA) and nine analogs during continuous dosing under the OECD test guideline 236. We further improved the predictions by estimating metabolic clearance and partition coefficients from the data by Bayesian calibration. We also assessed the performance of the model at reproducing data published by Brox et al. (2016) on VPA and 16 other chemicals. We finally compared dose-response relationships calculated for mortality and malformations on the basis of predicted whole embryo concentrations versus those based on nominal water concentrations. The use of target organ concentrations substantially shifted the magnitude of dose-response parameters and the relative toxicity ranking of chemicals studied.

Nanoplastics Cause Neurobehavioral Impairments, Reproductive and Oxidative Damages, and Biomarker Responses in Zebrafish: Throwing up Alarms of Wide Spread Health Risk of Exposure

Plastic pollution is a growing global emergency and it could serve as a geological indicator of the Anthropocene era. Microplastics are potentially more hazardous than macroplastics, as the former can permeate biological membranes. The toxicity of microplastic exposure on humans and aquatic organisms has been documented, but the toxicity and behavioral changes of nanoplastics (NPs) in mammals are scarce. In spite of their small size, nanoplastics have an enormous surface area, which bears the potential to bind even bigger amounts of toxic compounds in comparison to microplastics. Here, we used polystyrene nanoplastics (PS-NPs) (diameter size at ~70 nm) to investigate the neurobehavioral alterations, tissue distribution, accumulation, and specific health risk of nanoplastics in adult zebrafish. The results demonstrated that PS-NPs accumulated in gonads, intestine, liver, and brain with a tissue distribution pattern that was greatly dependent on the size and shape of the NPs particle. Importantly, an analysis of multiple behavior endpoints and different biochemical biomarkers evidenced that PS-NPs exposure induced disturbance of lipid and energy metabolism as well as oxidative stress and tissue accumulation. Pronounced behavior alterations in their locomotion activity, aggressiveness, shoal formation, and predator avoidance behavior were exhibited by the high concentration of the PS-NPs group, along with the dysregulated circadian rhythm locomotion activity after its chronic exposure. Moreover, several important neurotransmitter biomarkers for neurotoxicity investigation were significantly altered after one week of PS-NPs exposure and these significant changes may indicate the potential toxicity from PS-NPs exposure. In addition, after ~1-month incubation, the fluorescence spectroscopy results revealed the accumulation and distribution of PS-NPs across zebrafish tissues, especially in gonads, which would possibly further affect fish reproductive function. Overall, our results provided new evidence for the adverse consequences of PS-NPs-induced behavioral dysregulation and changes at the molecular level that eventually reduce the survival fitness of zebrafish in the ecosystem.

Multicompartmental Toxicokinetic Modeling of Discrete Dietary and Continuous Waterborne Uptake of Two Polycyclic Aromatic Hydrocarbons by Zebrafish Danio rerio

In the present study, we developed a multicompartmental toxicokinetic model for two polycyclic aromatic hydrocarbons (phenanthrene and anthracene) in their deuterated form (PAHs-d₁₀) in zebrafish considering continuous waterborne uptake and discrete dietary uptake. We quantified the bioconcentration, bioaccumulation, and depuration of these two PAHs-d₁₀ in zebrafish, and then estimated the kinetic parameters by fitting the model into the experimental data. The experimental and fitting results both showed that there was a peak concentration in each compartment of zebrafish after every dietary uptake, while the peak value depended on the ingestion amount of the PAH-d₁₀ and varied among different compartments. The PAH-d₁₀ amount in the blood reached 20-27% of the total amount bioaccumulated in zebrafish at steady-state, followed by skin (20-26%), and fillet (16-22%). The rank of PAH-d₁₀ steady-state concentrations in each compartment showed inconsistency with its lipid contents, indicating that the distribution of the PAHs-d₁₀ in zebrafish was not merely affected by the lipid content in each compartment, but also affected by their kinetics and biotransformation. This study suggests that discrete dietary uptake caused by intermittent food ingestion significantly affects the bioaccumulation of PAHs in fish. Further studies are needed to investigate such effect on other toxicants that are more resistant to biotransformation than PAHs in fish.

Different factors determined the toxicokinetics of organic chemicals and nanomaterials exposure to zebrafish (Danio Rerio)

Little is known about how the chemical properties (molecular structure, such as the hydrophobic and hydrophilic end group for organic chemical, and particle size for nanomaterials (NMs)) quantitatively affect the toxicokinetics (TK) in organisms especially in short-term, single-species studies. A novel method based on a first-order one compartment TK model which described the monophasic uptake pattern and two-compartment TK model which adequately described the biphasic metabolism pattern was used to determine the bioconcentration and TK rate constants of organic compounds (n = 17) and nanomaterials (NMs, n = 7) in zebrafish. For both one and two compartment model, the uptake (k_in) and elimination (k_out) rate constants were fitted using a one- and two-compartment first-order kinetic model, and bioconcentration factors (BCF) and 95% depuration times (t₉₅) for all tested chemicals were calculated, respectively. The results showed that there was significant difference in TK parameters k_in, k_out, and BCF between organic chemicals and nano metal oxides. For organic compounds, significant correlations were found between the k_in and BCF and the octanol-water partition coefficient (K_ow) and molecular mass. For nano metal oxides, there was a significant negative correlation between the k_in or BCF and particle size, but a positive correlation between k_in and Zeta potential of nanoparticles and also a significant positive correlation between k_out and particle size or specific surface area. Those findings indicated that NMs particle size does matter in biological influx and efflux processes. Our results suggest that the TK process for organic compound and NMs are correlated by different chemical properties and highlight that the K_ow, the absorption k_in, metabolism k₁₂ and k₂₁, elimination rate k_out, and all the parameters that enable the prediction and partitioning of chemicals need to be precisely determined in order to allow an effective TK modeling. It would therefore appear that the TK process of untested chemicals by a fish may be extrapolated from known chemical properties.

Elucidating the fate of perfluorooctanoate sulfonate using a rainbow trout (Oncorhynchus mykiss) physiologically-based toxicokinetic model

Per- and poly-fluorinated substances (PFAS) are widely found in freshwater ecosystems because of their resistance to degradation. Among them, several long-chain perfluoroalkyl acids bioaccumulate in aquatic vertebrates, but our understanding of the mechanisms of absorption, distribution and elimination is still limited in fish. For this purpose, we developed a 10-compartment physiologically-based toxicokinetic (PBTK) model to elucidate perfluorooctane sulfonate (PFOS) kinetics in adult rainbow trout. This PBTK model included various physiological characteristics: blood perfusion to each organ, plasmatic fraction, PFOS free fraction, and growth of individuals. The parameters were optimized using Bayesian inferences. First, only PFOS absorption by diet was considered in the model as well as its elimination by urine, bile and feces. Then two mechanistic hypotheses, assumed to govern PFOS toxicokinetics in fish, namely the enterohepatic cycle and the absorption and elimination though gills, were tested. Improvement of the model's fit to the data was studied in each organ by comparing predictions with observed data using relative error. The experimental data set was obtained from an exposure experiment, where adult rainbow trout were fed with a PFOS-spiked diet for 42 days, followed by a 35-day depuration period. In all cases, PFOS concentrations were accurately predicted in organs and feces by the model. The results of this PBTK model demonstrated that feces represented the major elimination route for PFOS while urine was a minor route. Also, PFOS branchial uptake can be substantial despite low concentrations of the compound in water, and elimination through gills should not be neglected. Finally, the enterohepatic cycle is likely to play a minor role in PFOS toxicokinetics. Overall, this PBTK model accurately described PFOS distribution in rainbow trout and provides information on the relative contribution of absorption and elimination pathways.

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.

Physiologically based toxicokinetics (PBTK) models for pharmaceuticals and personal care products in wild common carp (Cyprinus carpio)

Pharmaceutical and personal care products (PPCPs) are regarded as "pseudo-persistent" pollutants due to their being continuously loaded into the aquatic environment. Physiologically based toxicokinetics (PBTK) models that can quantitatively describe absorption, distribution, metabolism and excretion processes of chemicals in biota are of importance to predict internal exposure (e.g. doses at specific target tissues/organs) from external exposure concentrations. In this study, PBTK models with up to six compartments including brain, liver, kidney, gills, richly perfused tissues and poorly perfused tissues, were developed for predicting internal distribution of 10 PPCPs in wild common carp (Cyprinus carpio). The PBTK predicted concentrations were close to the measured ones, with deviations less than 1 log unit for most of PPCPs. Sensitivity analysis showed that various partition coefficients of the chemicals exerted significant influence on model outputs.

Investigating the interaction between melamine and cyanuric acid using a Physiologically-Based Toxicokinetic model in rainbow trout

Following outbreaks of feed and food adulterations with a melamine and cyanuric acid mixture in 2007 and melamine in 2008 respectively, the kinetics and toxicodynamics of the mixture have been investigated particularly in sensitive species such as the rainbow trout. Tissue concentrations and intensity of the adverse effect, melamine-cyanurate crystal formation in kidney, were reported in similar experimental conditions. Here, a recent PBTK model for rainbow trout has been applied to model the kinetics of both single compounds based on residue levels in tissues. Both PBTK models for the single compounds were combined and a model of crystal formation for the mixture melamine-cyanuric acid was also added to predict the intensity of crystal formation under the assumptions that crystals formed either in urine or in kidney tissue. Modelling the kinetics of melamine and cyanuric acid provided a better understanding and prediction of intensity of crystal formation in case of sequential exposures with varying intensity or co-exposure. This study demonstrates, for the first time, how fish PBTK models can play a key role in the understanding and prediction of toxicokinetics and toxicodynamics of mixtures. This study also illustrates how adverse effects may potentially occur even when the compounds are not administered together as a mixture.

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.

Statistical relationship between metabolic decomposition and chemical uptake predicts bioconcentration factor data for diverse chemical exposures

BACKGROUND

A challenge of in vitro to in vivo extrapolation (IVIVE) is to predict the physical state of organisms exposed to chemicals in the environment from in vitro exposure assay data. Although toxicokinetic modeling approaches promise to bridge in vitro screening data with in vivo effects, they are often encumbered by a need for redesign or re-parameterization when applied to different tissues or chemicals.

RESULTS

We demonstrate a parameterization of reverse toxicokinetic (rTK) models developed for the adult zebrafish (Danio rerio) based upon particle swarm optimizations (PSO) of the chemical uptake and degradation rates that predict bioconcentration factors (BCF) for a broad range of chemicals. PSO reveals a relationship between chemical uptake and decomposition parameter values that predicts chemical-specific BCF values with moderate statistical agreement to a limited yet diverse chemical dataset, and all without a need to retrain the model to new data.

CONCLUSIONS

The presented model requires only the octanol-water partitioning ratio to predict BCFs to a fidelity consistent with existing QSAR models. This success begs re-evaluation of the modeling assumptions; specifically, it suggests that chemical uptake into arterial blood may be limited by transport across gill membranes (diffusion) rather than by counter-current flow between gill lamellae (convection). Therefore, more detailed molecular modeling of aquatic respiration may further improve predictive accuracy of the rTK approach.

Importance of Toxicokinetics to Assess the Utility of Zebrafish Larvae as Model for Psychoactive Drug Screening Using Meta-Chlorophenylpiperazine (mCPP) as Example

The number of new psychoactive substances (NPS) increases rapidly, harming society and fuelling the need for alternative testing strategies. These should allow the ever-increasing number of drugs to be tested more effectively for their toxicity and psychoactive effects. One proposed strategy is to complement rodent models with zebrafish (Danio rerio) larvae. Yet, our understanding of the toxicokinetics in this model, owing to the waterborne drug exposure and the distinct physiology of the fish, is incomplete. We here explore the toxicokinetics and behavioral effects of an NPS, meta-chlorophenylpiperazine (mCPP), in zebrafish larvae. Uptake kinetics of mCPP, supported by toxicokinetic modeling, strongly suggested the existence of active transport processes. Internal distribution showed a dominant accumulation in the eye, implying that in zebrafish, like in mammals, melanin could serve as a binding site for basic drugs. We confirmed this by demonstrating significantly lower drug accumulation in two types of hypo-pigmented fish. Comparison of the elimination kinetics between mCPP and previously characterized cocaine demonstrated that drug affinities to melanin in zebrafish vary depending on the structure of the test compound. As expected from mCPP-elicited responses in rodents and humans, zebrafish larvae displayed hypoactive behavior. However, significant differences were seen between zebrafish and rodents with regard to the concentration-dependency of the behavioral response and the comparability of tissue levels, corroborating the need to consider the organism-internal distribution of the chemical to allow appropriate dose modeling while evaluating effects and concordance between zebrafish and mammals. Our results highlight commonalities and differences of mammalian versus the fish model in need of further exploration.

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.

Early life exposure to ethinylestradiol enhances subsequent responses to environmental estrogens measured in a novel transgenic zebrafish

Estrogen plays fundamental roles in a range of developmental processes and exposure to estrogen mimicking chemicals has been associated with various adverse health effects in both wildlife and human populations. Estrogenic chemicals are found commonly as mixtures in the environment and can have additive effects, however risk analysis is typically conducted for single-chemicals with little, or no, consideration given for an animal’s exposure history. Here we developed a transgenic zebrafish with a photoconvertable fluorophore (Kaede, green to red on UV light exposure) in a skin pigment-free mutant element (ERE)-Kaede-Casper model and applied it to quantify tissue-specific fluorescence biosensor responses for combinations of estrogen exposures during early life using fluorescence microscopy and image analysis. We identify windows of tissue-specific sensitivity to ethinylestradiol (EE2) for exposure during early-life (0–5 dpf) and illustrate that exposure to estrogen (EE2) during 0–48 hpf enhances responsiveness (sensitivity) to different environmental estrogens (EE2, genistein and bisphenol A) for subsequent exposures during development. Our findings illustrate the importance of an organism’s stage of development and estrogen exposure history for assessments on, and possible health risks associated with, estrogen exposure.

Evaluating parameter availability for physiologically based pharmacokinetic (PBPK) modeling of perfluorooctanoic acid (PFOA) in zebrafish

Physiologically based pharmacokinetic (PBPK) models are considered useful tools to describe the absorption, distribution, metabolism and excretion of xenobiotics. For accurate predictions, PBPK models require species-specific and compound-specific parameters. Zebrafish are considered an appropriate vertebrate model for investigating the toxicity of a wide variety of compounds. However, no specific mechanistic model exists for the pharmacokinetics of perfluoroalkyl acids (PFAAs) in zebrafish, despite growing concern about this class of ubiquitous environmental contaminants. The purpose of this study was to evaluate the current state of knowledge for the parameters that would be needed to construct such a model for zebrafish. We chose perfluorooctanoic acid (PFOA) as a model PFAA with greater data availability. We have updated a previous PBPK model for rainbow trout to simulate PFOA fate in zebrafish following waterborne exposure. For the first time, the model considers hepatobiliary circulation. In order to evaluate the availability of parameters to implement this model, we performed an extensive literature review to find zebrafish-specific parameters. As in previous approaches, we broadened our search to include mammalian and other fish studies when zebrafish-specific data were lacking. Based on the method used to measure or estimate parameters, or based on their species-specific origin, we scored and ranked the quality of available parameters. These scores were then used in Monte Carlo and partial rank correlation analyses to identify the most critical data gaps. The liver, where fatty acid binding proteins (FABPs) and plasma proteins are considered, represented the best model-data agreement. Lack of agreement in other tissues suggest better parameters are needed. The results of our study highlight the lack of zebrafish-specific parameters. Based on sensitivity and uncertainty analysis, parameters associated with PFAA-protein interactions and passive diffusion need further refinement to enable development of predictive models for these emerging chemicals in zebrafish.

In vitro to in vivo extrapolation for high throughput prioritization and decision making

In vitro chemical safety testing methods offer the potential for efficient and economical tools to provide relevant assessments of human health risk. To realize this potential, methods are needed to relate in vitro effects to in vivo responses, i.e., in vitro to in vivo extrapolation (IVIVE). Currently available IVIVE approaches need to be refined before they can be utilized for regulatory decision-making. To explore the capabilities and limitations of IVIVE within this context, the U.S. Environmental Protection Agency Office of Research and Development and the National Toxicology Program Interagency Center for the Evaluation of Alternative Toxicological Methods co-organized a workshop and webinar series. Here, we integrate content from the webinars and workshop to discuss activities and resources that would promote inclusion of IVIVE in regulatory decision-making. We discuss properties of models that successfully generate predictions of in vivo doses from effective in vitro concentration, including the experimental systems that provide input parameters for these models, areas of success, and areas for improvement to reduce model uncertainty. Finally, we provide case studies on the uses of IVIVE in safety assessments, which highlight the respective differences, information requirements, and outcomes across various approaches when applied for decision-making.

Bioaccumulation mechanism of organophosphate esters in adult zebrafish (Danio rerio)

Although organophosphate esters (OPEs) have been detected with growing frequency in water ecosystems, the underlying accumulation mechanisms of these compounds in fish are still unknown. Here, we investigated the tissue-specific accumulation and depuration of seven OPEs in adult zebrafish at three levels (0, 1/150 LC₅₀ (environmentally relevant level), and 1/30 LC₅₀ per OPE congener) in laboratory after 19 days exposure and 3 days depuration. The bioaccumulation of OPEs varied among tissues. Muscle contained the lowest level of OPEs and liver had the highest level of two (TPP and TCEP) of the seven OPEs at steady state. The high levels and slow depuration rates of TDCIPP, TPHP, and TCP observed in roe indicated that the accumulated OPEs were potentially stored in roe and transferred to the next generation. After examination of the major metabolites (organophosphate diesters) in selected tissues, a physiologically based toxicokinetic (PBTK) model used in fish was adopted to explore the key factors affecting the bioaccumulation of OPEs in zebrafish. Biotransformation of OPEs with polychlorinated alkyl moieties (i.e. TDCIPP) and aryl moieties (i.e. TPHP and TCP) has more significant impacts on the accumulation than those of OPEs with alkyl or short chain chlorinated alkyl moieties. Furthermore, the partition process between tissues and blood was also investigated, and was demonstrated to be the dominant process for OPEs accumulation in zebrafish. This study provides critical information on the bioaccumulation, tissue distribution, and metabolization of OPEs in relation with OPE structures in fish, as well as the underlying bioaccumulation mechanisms/pathways of OPEs in aquatic life.

Aquatic ecotoxicological models and their applicability in Arctic regions

Dose-response modeling is one of the most important steps of ecological risk assessment. It requires concentration-effects relationships for the species under consideration. There are very limited studies and experimental data available for the Arctic aquatic species. Lack of toxicity data hinders obtaining dose-response relationships for lethal (LC50 values), sub-lethal and carcinogenic effects. Gaps in toxicity data could be filled using a variety of in-silico ecotoxicological methods. This paper reviews the suitability of such methods for the Arctic scenario. Mechanistic approaches like toxicokinetic and toxicodynamic analysis are found to be better suited for interspecies extrapolation than statistical methods, such as Quantitative Structure-Activity Relationships/Quantitative Structure Activity-Activity Relationship, ICE, and other empirical models, such as Haber's law and Ostwald's equation. A novel approach is proposed where the effects of the toxicant exposure are quantified based on the probability of cellular damage and metabolites interactions. This approach recommends modeling cellular damage using a toxicodynamic model and physiology or metabolites interactions using a toxicokinetic model. Together, these models provide more reliable estimates of toxicity in the Arctic aquatic species, which will assist in conducting ecological risk assessment of Arctic environment.

A toxicokinetic model for fish including multiphase sorption features

In the scientific field of physiologically based toxicokinetic modeling the complexity of the model used depends on the complexity of the problem to be handled, leading to a broad range of existing models from simple 1-box models to complex multicompartment models. Most of these models work with lumped parameters, for example, an uptake efficiency parameter that can only be obtained with a fit of experimental data. The authors' goal was a model that is completely based on well-defined physiological and physicochemical parameters. Lumped parameters fitted on training data sets would limit the model's applicability. This would enable a new view on process understanding for uptake, distribution, and elimination procedures. Eventual goals are a better localization of chemicals within the organism itself, and to set the stage for future extensions toward ionic compounds and active transport across membranes. The model evaluation reported in the present study has shown that uptake, clearance, and bioaccumulation data for nonpolar chemicals are well predicted. Environ Toxicol Chem 2017;36:1538-1546. © 2016 SETAC.

Comparison of the In Vivo Biotransformation of Two Emerging Estrogenic Contaminants, BP2 and BPS, in Zebrafish Embryos and Adults

Zebrafish embryo assays are increasingly used in the toxicological assessment of endocrine disruptors. Among other advantages, these models are 3R-compliant and are fit for screening purposes. Biotransformation processes are well-recognized as a critical factor influencing toxic response, but major gaps of knowledge exist regarding the characterization of functional metabolic capacities expressed in zebrafish. Comparative metabolic studies between embryos and adults are even scarcer. Using ³H-labeled chemicals, we examined the fate of two estrogenic emerging contaminants, benzophenone-2 (BP2) and bisphenol S (BPS), in 4-day embryos and adult zebrafish. BPS and BP2 were exclusively metabolized through phase II pathways, with no major qualitative difference between larvae and adults except the occurrence of a BP2-di-glucuronide in adults. Quantitatively, the biotransformation of both molecules was more extensive in adults. For BPS, glucuronidation was the predominant pathway in adults and larvae. For BP2, glucuronidation was the major pathway in larvae, but sulfation predominated in adults, with ca. 40% conversion of parent BP2 and an extensive release of several conjugates into water. Further larvae/adults quantitative differences were demonstrated for both molecules, with higher residue concentrations measured in larvae. The study contributes novel data regarding the metabolism of BPS and BP2 in a fish model and shows that phase II conjugation pathways are already functional in 4-dpf-old zebrafish. Comparative analysis of BP2 and BPS metabolic profiles in zebrafish larvae and adults further supports the use of zebrafish embryo as a relevant model in which toxicity and estrogenic activity can be assessed, while taking into account the absorption and fate of tested substances.

Physiological fidelity or model parsimony? The relative performance of reverse-toxicokinetic modeling approaches

Background

Physiologically-based toxicokinetic (PBTK) models are often developed to facilitate in vitro to in vivo extrapolation (IVIVE) using a top-down, compartmental approach, favoring architectural simplicity over physiological fidelity despite the lack of general guidelines relating model design to dynamical predictions. Here we explore the impact of design choice (high vs. low fidelity) on chemical distribution throughout an animal’s organ system.

Results

We contrast transient dynamics and steady states of three previously proposed PBTK models of varying complexity in response to chemical exposure. The steady states for each model were determined analytically to predict exposure conditions from tissue measurements. Steady state whole-body concentrations differ between models, despite identical environmental conditions, which originates from varying levels of physiological fidelity captured by the models. These differences affect the relative predictive accuracy of the inverted models used in exposure reconstruction to link effects-based exposure data with whole-organism response thresholds obtained from in vitro assay measurements.

Conclusions

Our results demonstrate how disregarding physiological fideltiy in favor of simpler models affects the internal dynamics and steady state estimates for chemical accumulation within tissues, which, in turn, poses significant challenges for the exposure reconstruction efforts that underlie many IVIVE methods. Developing standardized systems-level models for ecological organisms would not only ensure predictive consistency among future modeling studies, but also ensure pragmatic extrapolation of in vivo effects from in vitro data or modeling exposure-response relationships.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0407-3) contains supplementary material, which is available to authorized users.

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.

Tissue-Specific Accumulation, Depuration, and Transformation of Triphenyl Phosphate (TPHP) in Adult Zebrafish (Danio rerio)

Understanding bioaccumulation and metabolism is critical for evaluating the fate and potential toxicity of compounds in vivo. We recently investigated, for the first time, the bioconcentration and tissue distribution of triphenyl phosphate (TPHP) and its main metabolites in selected tissues of adult zebrafish. To further confirm the metabolites, deuterated TPHP (d₁₅-TPHP) was used in the exposure experiments at an environmentally relevant level (20 μg/L) and at 1/10 LC₅₀ (100 μg/L). After 11-14 days of exposure to 100 μg/L of d₁₅-TPHP, the accumulation and excretion of d₁₅-TPHP reached equilibrium, at which point the intestine contained the highest d₁₅-TPHP (μg/g wet weight, ww) concentration (3.12 ± 0.43), followed by the gills (2.76 ± 0.12) > brain (2.58 ± 0.19) > liver (2.30 ± 0.34) ≫ muscle (0.53 ± 0.04). The major metabolite of d₁₅-TPHP, d₁₀-diphenyl phosphate (d₁₀-DPHP), was detected at significantly higher contents in the liver and intestine, at levels up to 3.0-3.5 times those of d₁₅-TPHP. The metabolic pathways of TPHP were elucidated, including hydrolysis, hydroxylation, and glucuronic acid conjugation after hydroxylation. Finally, a physiologically based toxicokinetic (PBTK) model was used to explore the key factors influencing the bioaccumulation of d₁₅-TPHP in zebrafish. These results provide important information for the understanding of the metabolism, disposition, and toxicology of TPHP in aquatic organisms.

Bioavailability and bioconcentration potential of perfluoroalkyl-phosphinic and -phosphonic acids in zebrafish (Danio rerio): Comparison to perfluorocarboxylates and perfluorosulfonates

Currently, information regarding bioavailability and bioconcentration potential of perfluoroalkyl phosphinic acids (PFPiAs) in aquatic organisms does not exist. The main objective of the present study was to assess uptake and elimination kinetics of PFPiAs in zebrafish (Danio rerio) following aqueous exposure. The results showed that PFPiA exposure can result in very high steady-state bioconentration factors (BCFss), compared to perfluorocarboxylates and perfluorosulfonates.C6/C10 PFPiA exhibited the highest BCFss, ranging between 10(7) and 10(10), orders of magnitude higher than those for long-chain perfluorocarboxylates. Strong positive relationships were observed between BCFss versus the membrane-water distribution coefficient (Dmw) and the protein-water partition coefficient (Kpw) of the studied perfluoroalkyl substances. However, BCFss exhibited a substantial drop for the very hydrophobic PFPiAs (C8/C10 and C6/C12 PFPiAs). The reduced BCFss of these long-chain PFPiAs (perfluoroalkyl chain length=18; Dmw=10(9)) is likely the result of reduced bioavailability due to interaction with solute molecules/organic matter present in the water phase and/or reduced gill membrane permeability. While PFPiAs can be metabolized to perfluoroalkyl phosphonic acids, the metabolic transformation rate seems insufficient to counteract the high degree of uptake across gill membranes. These findings help to better understand exposure pathways and bioaccumulation behavior of these important perfluorinated acids in aquatic systems.

Perspectives for integrating human and environmental exposure assessments

Integrated Risk Assessment (IRA) has been defined by the EU FP7 HEROIC Coordination action as "the mutual exploitation of Environmental Risk Assessment for Human Health Risk Assessment and vice versa in order to coherently and more efficiently characterize an overall risk to humans and the environment for better informing the risk analysis process" (Wilks et al., 2015). Since exposure assessment and hazard characterization are the pillars of risk assessment, integrating Environmental Exposure assessment (EEA) and Human Exposure assessment (HEA) is a major component of an IRA framework. EEA and HEA typically pursue different targets, protection goals and timeframe. However, human and wildlife species also share the same environment and they similarly inhale air and ingest water and food through often similar overlapping pathways of exposure. Fate models used in EEA and HEA to predict the chemicals distribution among physical and biological media are essentially based on common properties of chemicals, and internal concentration estimations are largely based on inter-species (i.e. biota-to-human) extrapolations. Also, both EEA and HEA are challenged by increasing scientific complexity and resources constraints. Altogether, these points create the need for a better exploitation of all currently existing data, experimental approaches and modeling tools and it is assumed that a more integrated approach of both EEA and HEA may be part of the solution. Based on the outcome of an Expert Workshop on Extrapolations in Integrated Exposure Assessment organized by the HEROIC project in January 2014, this paper identifies perspectives and recommendations to better harmonize and extrapolate exposure assessment data, models and methods between Human Health and Environmental Risk Assessments to support the further development and promotion of the concept of IRA. Ultimately, these recommendations may feed into guidance showing when and how to apply IRA in the regulatory decision-making process for chemicals.

A fugacity-based toxicokinetic model for narcotic organic chemicals in fish

A novel dynamic fugacity-based model is described, developed, and tested that simulates the uptake of narcotic organic chemicals in fish from water as occurs in aquatic bioconcentration and toxicity tests. The physiologically based toxicokinetic model treats the time course of chemical distribution in 4 compartments (tissue groups) in the fish, including the liver, in which biotransformation may occur. In addition to calculating bioconcentration and toxicokinetics, 5 possible toxic endpoints are defined corresponding to chemical concentration, fugacity, or activity reaching a critical value that causes 50% mortality. The mathematical description of multicompartment uptake is simplified by expressing the equations in the fugacity format. The model is parameterized and tested against reported empirical data for the bioconcentration of pentachloroethane in rainbow trout and for uptake and mortality from aquatic exposures to naphthalene and 1,2,4-trichlorobenzene in fathead minnows. Model performance is evaluated, and it is concluded that with suitable parameterization it has potential for application for assessment of both bioconcentration and toxicity expressed as median lethal concentrations, critical body residues, and chemical activity as a function of time to death.

In situ experiments to assess effects of constraints linked to caging on ecotoxicity biomarkers of the three-spined stickleback (Gasterosteus aculeatus L.)

The aim of this study was to evaluate the effects of caging constraints on multiple fish biomarkers used during ecotoxicological studies (biometric data, immune and antioxidant systems, and energetic status). Two of these constraints were linked to caging: starvation and fish density in cages, and one in relation to the post-caging handling: a short transport. Three in situ experiments were conducted with three-spined sticklebacks (Gasterosteus aculeatus). The first experiment compared the effects of three densities (low, medium, and high). The second experiment compared effects of starvation in fish fed every two days with fish that were not fed. Finally comparisons between sticklebacks which have suffered a short car transport after caging and sticklebacks killed without preliminary transport were made. The lack of food had no effect on fish energetic reserves but negatively affected their condition index and their immune system. Transport and high density induced oxidative stress, defined as an overproduction of reactive oxygen species and a stimulation of the antioxidant system. These two constraints also harmed the leucocyte viability. In order not to have any impact on ecotoxicity biomarkers during in situ experiments, it is preferable to decrease fish density in cages, prevent transport before dissections, and feed fish when the caging lasts more than two weeks.