Predicting concentrations of organic chemicals in fish by using toxicokinetic models

Application of toxicokinetic-toxicodynamic models in the aquatic ecological risk assessment of metals: A review

The issue of toxic metal pollution is a considerable environmental concern owing to its complex nature, spatial and temporal variability, and susceptibility to environmental factors. Current water quality criteria and ecological risk assessments of metals are based on single-metal toxicity data from short-term, simplified indoor exposure conditions, ignoring the complexity of actual environmental conditions. This results in increased uncertainty in predicting toxic metal toxicity and risk assessment. Using appropriate bioavailability and effect modeling of metals is critical for establishing environmental quality standards and performing risk assessments for metals. Traditional dose-effect models are based on a static statistical relationship and fall short of revealing the bioavailability and effect processes of metals and do not effectively assess ecological impacts under complex exposure conditions. This paper summarizes the toxicokinetic-toxicodynamic (TK-TD) model, which is gaining interest in environmental and ecotoxicological research. The key concepts, and theories of its construction theories, are discussed and the application of the TK-TD model in toxicity prediction and risk assessment of different metals in the aquatic environment, and trends in the development of the TK-TD model are highlighted. The findings of our review prove that the TK-TD model can effectively predict toxic metal toxicity in real time and under complex exposure conditions in the future.

Toxicokinetics of selenate in earthworm sub-tissues and potential bio-accessibility assessment of earthworm-derived selenium

Selenium (Se) pollution is mainly caused by anthropogenic activities, and the resulting biosecurity concerns have garnered significant attention in recent years. Using one-compartmental toxicokinetic (TK) modelling, this study explored the kinetic absorption, sub-tissue distribution, and elimination processes of the main Se species (selenate, Se(VI)) in the cultivated aerobic soil of the earthworm Eisenia fetida. The bio-accessibility of earthworm-derived Se was assessed using an in vitro simulated gastrointestinal digestion test to evaluate its potential trophic risk. The results demonstrated that Se accumulated in the pre-clitellum (PC) and total tissues (TT) of earthworms in a time- and dose-dependent manner. The highest Se levels in the PC, post-clitellum (PoC), and TT were 70.54, 57.93, and 64.26 mg/kg during the uptake phase, respectively. The kinetic Se contents in the earthworms PC and TT were consistent with the TK model but not with PoC. The earthworm TT exhibited a faster uptake (Kus = 0.83-1.02 mg/kg/day) and elimination rate of Se (Kee = 0.044-0.049 mg/kg/day), as well as a shorter half-life time (LT1/2 = 15.88-14.22 days) than PC at low soil Se levels (≤5 mg/kg). Conversely, the opposite trend was observed with higher Se concentrations (10 and 20 mg/kg). These results are likely attributable to the tissue specificity and concentration of the toxicant. Earthworms PC and TT exhibited a higher kinetic Se accumulation factor (BAFk) than steady-state BAF (BAFss), with values ranging from 8 to 24 and 3-13, respectively. Furthermore, the bio-accessibility of earthworm-derived Se to poultry ranged from 66.25 % to 84.35 %. As earthworms are at the bottom of the terrestrial food chain, the high bio-accessibility of earthworm-derived Se poses a potential risk to predators. This study offers data support and a theoretical foundation for understanding the biological footprint of soil Se and its toxicological impacts and ecological hazards.

Recent Progress on Physiologically Based Pharmacokinetic (PBPK) Model: A Review Based on Bibliometrics

Physiologically based pharmacokinetic/toxicokinetic (PBPK/PBTK) models are designed to elucidate the mechanism of chemical compound action in organisms based on the physiological, biochemical, anatomical, and thermodynamic properties of organisms. After nearly a century of research and practice, good results have been achieved in the fields of medicine, environmental science, and ecology. However, there is currently a lack of a more systematic review of progress in the main research directions of PBPK models, especially a more comprehensive understanding of the application in aquatic environmental research. In this review, a total of 3974 articles related to PBPK models from 1996 to 24 March 2024 were collected. Then, the main research areas of the PBPK model were categorized based on the keyword co-occurrence maps and cluster maps obtained by CiteSpace. The results showed that research related to medicine is the main application area of PBPK. Four major research directions included in the medical field were "drug assessment", "cross-species prediction", "drug-drug interactions", and "pediatrics and pregnancy drug development", in which "drug assessment" accounted for 55% of the total publication volume. In addition, bibliometric analyses indicated a rapid growth trend in the application in the field of environmental research, especially in predicting the residual levels in organisms and revealing the relationship between internal and external exposure. Despite facing the limitation of insufficient species-specific parameters, the PBPK model is still an effective tool for improving the understanding of chemical-biological effectiveness and will provide a theoretical basis for accurately assessing potential risks to ecosystems and human health. The combination with the quantitative structure-activity relationship model, Bayesian method, and machine learning technology are potential solutions to the previous research gaps.

Cross-Species Extrapolation of Biological Data to Guide the Environmental Safety Assessment of Pharmaceuticals-The State of the Art and Future Priorities

The extrapolation of biological data across species is a key aspect of biomedical research and drug development. In this context, comparative biology considerations are applied with the goal of understanding human disease and guiding the development of effective and safe medicines. However, the widespread occurrence of pharmaceuticals in the environment and the need to assess the risk posed to wildlife have prompted a renewed interest in the extrapolation of pharmacological and toxicological data across the entire tree of life. To address this challenge, a biological "read-across" approach, based on the use of mammalian data to inform toxicity predictions in wildlife species, has been proposed as an effective way to streamline the environmental safety assessment of pharmaceuticals. Yet, how effective has this approach been, and are we any closer to being able to accurately predict environmental risk based on known human risk? We discuss the main theoretical and experimental advancements achieved in the last 10 years of research in this field. We propose that a better understanding of the functional conservation of drug targets across species and of the quantitative relationship between target modulation and adverse effects should be considered as future research priorities. This pharmacodynamic focus should be complemented with the application of higher-throughput experimental and computational approaches to accelerate the prediction of internal exposure dynamics. The translation of comparative (eco)toxicology research into real-world applications, however, relies on the (limited) availability of experts with the skill set needed to navigate the complexity of the problem; hence, we also call for synergistic multistakeholder efforts to support and strengthen comparative toxicology research and education at a global level. Environ Toxicol Chem 2024;43:513-525. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Studying mixture effects on uptake and tissue distribution of PFAS in zebrafish (Danio rerio) using physiologically based kinetic (PBK) modelling

Per- and polyfluoroalkyl substances (PFAS) are ubiquitously distributed in the aquatic environment. They include persistent, mobile, bioaccumulative, and toxic chemicals and it is therefore critical to increase our understanding on their adsorption, distribution, metabolism, excretion (ADME). The current study focused on uptake of seven emerging PFAS in zebrafish (Danio rerio) and their potential maternal transfer. In addition, we aimed at increasing our understanding on mixture effects on ADME by developing a physiologically based kinetic (PBK) model capable of handling co-exposure scenarios of any number of chemicals. All studied chemicals were taken up in the fish to varying degrees, whereas only perfluorononanoate (PFNA) and perfluorooctanoate (PFOA) were quantified in all analysed tissues. Perfluorooctane sulfonamide (FOSA) was measured at concerningly high concentrations in the brain (Cmax over 15 μg/g) but also in the liver and ovaries. All studied PFAS were maternally transferred to the eggs, with FOSA and 6:2 perfluorooctane sulfonate (6,2 FTSA) showing significant (p < 0.02) signs of elimination from the embryos during the first 6 days of development, while perfluorobutane sulfonate (PFBS), PFNA, and perfluorohexane sulfonate (PFHxS) were not eliminated in embryos during this time-frame. The mixture PBK model resulted in >85 % of predictions within a 10-fold error and 60 % of predictions within a 3-fold error. At studied levels of PFAS exposure, competitive binding was not a critical factor for PFAS kinetics. Gill surface pH influenced uptake for some carboxylates but not the sulfonates. The developed PBK model provides an important tool in understanding kinetics under complex mixture scenarios and this use of New Approach Methodologies (NAMs) is critical in future risk assessment of chemicals and early warning systems.

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.

Multimodal Model to Predict Tissue-to-Blood Partition Coefficients of Chemicals in Mammals and Fish

Tissue-to-blood partition coefficients (Ptb) are key parameters for assessing toxicokinetics of xenobiotics in organisms, yet their experimental data were lacking. Experimental methods for measuring Ptb values are inefficient, underscoring the urgent need for prediction models. However, most existing models failed to fully exploit Ptb data from diverse sources, and their applicability domain (AD) was limited. The current study developed a multimodal model capable of processing and integrating textual (categorical features) and numerical information (molecular descriptors/fingerprints) to simultaneously predict Ptb values across various species, tissues, blood matrices, and measurement methods. Artificial neural network algorithms with embedding layers were used for the multimodal modeling. The corresponding unimodal models were developed for comparison. Results showed that the multimodal model outperformed unimodal models. To enhance the reliability of the model, a method considering categorical features, weighted molecular similarity density, and weighted inconsistency in molecular activities of structure-activity landscapes was used to characterize the AD. The model constrained by the AD exhibited better prediction accuracy for the validation set, with the determination coefficient, root mean-square error, and mean absolute error being 0.843, 0.276, and 0.213 log units, respectively. The multimodal model coupled with the AD characterization can serve as an efficient tool for internal exposure assessment of chemicals.

Bioconcentration of polycyclic aromatic hydrocarbons in different tissues of zebrafish (Danio rerio) investigated with PBTK model

Due to the lipophilicity, polycyclic aromatic hydrocarbons (PAHs) are easily accumulated in fish. However, the research on PAH bioaccumulation process in different fish tissues and the relevant effect mechanisms are still deficient. The bioconcentration of PAHs (phenanthrene, anthracene, fluoranthene, and pyrene) in different zebrafish tissues (skin, fish muscle, gill, digestive tract, liver, gonad, and residual) was studied. It was found that there was a difference in the PAH concentrations in different zebrafish tissues. Compared with other tissues, the PAH concentration was highest in the skin and lowest in the fish muscle. For example, the steady-state concentration of phenanthrene in the skin was nearly five times higher than that in the muscle. PAH distribution was related with the lipid contents in different zebrafish tissues; however, the correlation was not significant (p > 0.05), indicating that the lipid content was not the determining factor for the PAH distribution. The distribution was also affected by the bioconcentration kinetics of PAHs in different zebrafish tissues, and the PAH hydrophobic properties. In addition, the physiological based toxicokinetic (PBTK) model showed good performance in predicting PAH internal concentrations, and it may be used to predict the concentrations of PAHs in different fish tissues in future.

Role of Gastrointestinal Microbiota from Crucian Carp in Microbial Transformation and Estrogenicity Modification of Novel Plastic Additives

Ingestion is a major exposure route for hydrophobic organic pollutants in fish, but the microbial transformation and estrogenic modification of the novel plastic additives by the gut microbiota of fish remain obscure. Using an in vitro approach, we provide evidence that structure-related transformation of various plastic additives by the gastric and intestinal (GI) microbiota from crucian carp, with the degradation ratio of bisphenols and triphenyl phosphate faster than those of brominated compounds. The degradation kinetics for these pollutants could be limited by oxygen and cometabolic substrates (i.e., glucose). The fish GI microbiota could utilize the vast majority of carbon sources in a Biolog EcoPlate, suggesting their high metabolic potential and ability to transform various organic compounds. Unique microorganisms associated with transformation of the plastic additives including genera of Citrobacter, Klebsiella, and some unclassified genera in Enterobacteriaceae were identified by combining high-throughput genetic analyses and metagenomic analyses. Through identification of anaerobic transformation products by high-resolution mass spectrometry, alkyl-cleavage was found the common transformation mechanism, and hydrolysis was the major pathway for ester-containing pollutants. After anaerobic incubation, the estrogenic activities of triphenyl phosphate and bisphenols A, F, and AF declined, whereas that of bisphenol AP increased.

New Approach Methodologies for the Endocrine Activity Toolbox: Environmental Assessment for Fish and Amphibians

Multiple in vivo test guidelines focusing on the estrogen, androgen, thyroid, and steroidogenesis pathways have been developed and validated for mammals, amphibians, or fish. However, these tests are resource-intensive and often use a large number of laboratory animals. Developing alternatives for in vivo tests is consistent with the replacement, reduction, and refinement principles for animal welfare considerations, which are supported by increasing mandates to move toward an "animal-free" testing paradigm worldwide. New approach methodologies (NAMs) hold great promise to identify molecular, cellular, and tissue changes that can be used to predict effects reliably and more efficiently at the individual level (and potentially on populations) while reducing the number of animals used in (eco)toxicological testing for endocrine disruption. In a collaborative effort, experts from government, academia, and industry met in 2020 to discuss the current challenges of testing for endocrine activity assessment for fish and amphibians. Continuing this cross-sector initiative, our review focuses on the current state of the science regarding the use of NAMs to identify chemical-induced endocrine effects. The present study highlights the challenges of using NAMs for safety assessment and what work is needed to reduce their uncertainties and increase their acceptance in regulatory processes. We have reviewed the current NAMs available for endocrine activity assessment including in silico, in vitro, and eleutheroembryo models. New approach methodologies can be integrated as part of a weight-of-evidence approach for hazard or risk assessment using the adverse outcome pathway framework. The development and utilization of NAMs not only allows for replacement, reduction, and refinement of animal testing but can also provide robust and fit-for-purpose methods to identify chemicals acting via endocrine mechanisms. Environ Toxicol Chem 2023;42:757-777. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Uptake, distribution, and elimination of selenite in earthworm Eisenia fetida at sublethal concentrations based on toxicokinetic model

Natural and anthropogenic causes have promoted the rapid increase in environmental selenium (Se) levels, and the complex Se metabolism and dynamic in organisms make it challenging to evaluate the toxicity and ecological risks. In this study, the kinetics of selenite in earthworm Eisenia fetida were investigated based on toxicokinetic (TK) model (uptake-elimination phases: 14-14 days). The results showed the highest sub-tissue Se concentrations in pre-clitellum (PC), post-clitellum (PoC) parts, and total earthworms were 95.71, 70.40, and 79.94 mg/kg, respectively, which indicates the distinctive Se uptake capacities of E. fetida. Se kinetic rates in PCs were faster than that of the total E. fetida for both uptake (K_us = 0.30-0.80 mg/kg/day) and elimination phases (K_ee = 0.024-0.056 mg/kg/day). Longer half-life times (LT_1/2) were observed in the total earthworms (17.85-47.15 d) than PCs (12.28-29.22 d), while non-significant difference was found for the kinetic Se bioaccumulation factor (BAF_k) in PC and total earthworms (12-19), which demonstrates that Se can be efficiently bioaccumulated and eliminated in earthworm PC part. Besides, the significant increase Se concentration in PoC with rapid elimination in PC also illustrates that earthworms can alleviate the Se stress by the transformation strategy of Se from the head to tail tissues. In conclusion, the investigation of Se kinetic accumulation and elimination characteristics in this study is helpful for understanding the metabolism and detoxification processes of Se in earthworms, and also providing a theoretical basis for further Se risk assessment using TK model.

Synthesizing the impact of sea-dumped munition and related chemicals on humans and the environment

Marine environments are globally impacted by vast quantities of munition disposed following both World Wars. Dumped munitions contain conventional explosives, chemicals warfare agents as well as a variety of metals. Field monitoring studies around marine dumpsites report the presence of munition constituents in water and sediment samples. The growing interest and developments in the ocean as a new economic frontier underline the need to remediate existing dumpsites. Here, we provide a comprehensive assessment of the magnitude and potential risks associated with marine munition dumpsites. An overview of the global distribution of dumpsites identifying the most impacted areas is provided, followed by the currently available data on the detection of munition constituents in environmental samples and evidence of their toxic potential to human and environmental health. Finally, existing data gaps are identified and future research needs promoting better understanding of the impact of the dumped material on the marine environment suggested.

Maternal Transfer and Apical and Physiological Effects of Dietary Hexabromocyclododecane Exposure in Parental Fathead Minnows (Pimephales promelas)

Hexabromocyclododecane (HBCD) is a persistent organic pollutant that has been characterized as an endocrine disruptor, undergoes maternal transfer, and hinders development and growth in oviparous organisms. The present study examined the apical effects of dietary HBCD (11.5, 36.4, 106 mg/kg, wet wt) on adult fathead minnow exposed for 49 days and the subsequent accumulation and maternal transfer kinetics in adult tissue and eggs, respectively. Exposed adults displayed a significant increase in egg production in the medium treatment group, but no other significant effects were noted. Maternal transfer of dietary HBCD had a similar egg-to-muscle ratios (EMR) in the low and medium treatment groups (1.65 and 1.27 [wet wt], respectively). However, the high treatment group deviated from other treatments with an EMR of 4.2 (wet wt), potentially due to differences in total lipid content in food and/or reaching diffusion/lipid saturation limits in adult tissue, resulting in lower accumulation in the adult muscle tissue. A positive correlation was observed between egg HBCD concentration and time of exposure, which indicates that maternal transfer of HBCD is of concern in fish, and further studies should be conducted to fully elucidate the potential adverse effects that may be observed in the early life stage of oviparous organisms. Environ Toxicol Chem 2023;42:143-153. © 2022 SETAC.

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.

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.

Uptake, elimination, and toxicokinetics of selected pharmaceuticals in multiple tissues of Nile tilapia (Oreochromis niloticus) exposed to environmentally relevant concentrations

Pharmaceuticals in aquatic environment displayed adverse effects to fish. The effects are usually related to the internal levels of pharmaceuticals accumulated in specific fish tissues. In this study, we investigated the uptake, elimination, and toxicokinetics of six pharmaceuticals, e.g. naproxen (NAX), diclofenac (DCF), ibuprofen (IBU), carbamazepine (CBZ), fluoxetine (FLX), and sertraline (SER), in 11 fish tissues of Nile tilapia. The experiments were conducted in a flow-through system with an 8-day uptake/8-day elimination periods. The fish exposure groups involved the control, single FLX, and mixture of six pharmaceuticals at environmentally relevant concentration of 4 μg/L. FLX and SER showed the maximum concentrations of 145 and 201 ng/g wet weight, respectively, in fish spleen tissue, while NAX and IBU were not detected in any tissue. The mean concentrations for the pharmaceuticals in Nile tilapia tissues generally followed the order: bile> kidney, gut, stomach, liver> brain, gill, spleen> plasma, skin, muscle. The steady-state bioconcentration factors in various tissues generally range at 0.74-437.58 L/kg. The uptake and elimination toxicokinetics illustrated the rapid accumulation and depuration of pharmaceuticals in fish tissues. The results help to understand the internal bioconcentration, tissue distribution, and toxicokinetics of pharmaceuticals in multiple fish biological compartments.

One and multi-compartments toxico-kinetic modeling to understand metals' organotropism and fate in Gammarus fossarum

The use of freshwater invertebrates for biomonitoring has been increasing for several decades, but little is known about relations between external exposure concentration of metals and their biodistribution among different tissues. One and multi-compartments toxicokinetic (TK) models are powerful tools to formalize and predict how a contaminant is bioaccumulated. The aim of this study is to develop modeling approaches to improve knowledge on dynamic of accumulation and fate of Cd and Hg in gammarid's organs. Gammarids were exposed to dissolved metals (11.1 ± 1.2 µg.L^-1 of Cd or 0.27 ± 0.13 µg.L^-1 of Hg) before a depuration phase. At each sampling days, their organs (caeca, cephalon, intestine and remaining tissues) were separated by dissection before analyses. Results allowed us to determine that i) G.fossarum takes up Cd as efficiently as the mussel M.galloprovincialis, but eliminates it more rapidly, ii) organs which accumulate and depurate the most, in terms of concentrations, are caeca and intestine for both metals; iii) the one-compartment TK models is the most relevant for Hg, while the multi-compartments TK model allows a better fit to Cd data, demonstrating dynamic transfer of Cd among organs.

Antibiotics in mariculture organisms of different growth stages: Tissue-specific bioaccumulation and influencing factors

Maricultured organisms are chronically exposed to water containing antibiotics but the bioaccumulative behavior of antibiotics in exposed organisms at different growth stages has received little attention. Here, we investigated the concentrations and tissue-specific bioaccumulation characteristics of 19 antibiotics during three growth stages (youth stage, growth stage, and adult stage) of various organisms (Scophthalmus maximus, Penaeus vannamei, Penaeus japonicus, and Apostichopus japonicus) cultivated in typical marine aquaculture regions, and explored the factors that could affect the bioaccumulation of antibiotics. Tetracyclines (TCs) and fluoroquinolones (FQs) were the dominant antibiotics in all organisms, and the total concentrations of the target antibiotics in fish (S. maximus) were significantly higher than those in shrimp (P. vannamei and P. japonicus) and sea cucumber (A. japonicus) (p < 0.01). The bioaccumulation capacity of a class of statistically significant antibiotics in most samples was strongest during the youth stage and weakest during the adult stage. The antibiotics exhibited higher bioaccumulation capacity in lipid-rich tissues (fish liver and shrimp head) or respiratory organs (fish gill) than muscle. Our results also reveal significant metabolic transformation of enrofloxacin in fish. Different from previous studies, the logarithm bioaccumulation factor (log BAF) was positively correlated with log D_lipw in low-biotransformation tissues (fish gill and muscle) rather than lipid-rich tissues (fish liver). Based on the calculated hazard quotients (HQ), doxycycline in fish muscle may pose a distinct risk to human health, which deserves special attention. Overall, these results provide insight into the bioaccumulation patterns of antibiotics during different growth stages and tissues of maricultured organisms.

Toxicokinetic Models for Bioconcentration of Organic Contaminants in Two Life Stages of White Sturgeon (Acipenser transmontanus)

The white sturgeon (Acipenser transmontanus) is an endangered ancient fish species that is known to be particularly sensitive to certain environmental contaminants, partly because of the uptake and subsequent toxicity of lipophilic pollutants prone to bioconcentration as a result of their high lipid content. To better understand the bioconcentration of organic contaminants in this species, toxicokinetic (TK) models were developed for the embryo-larval and subadult life stages. The embryo-larval model was designed as a one-compartment model and validated using whole-body measurements of benzo[a]pyrene (B[a]P) metabolites from a waterborne exposure to B[a]P. A physiologically based TK (PBTK) model was used for the subadult model. The predictive power of the subadult model was validated with an experimental data set of four chemicals. Results showed that the TK models could accurately predict the bioconcentration of organic contaminants for both life stages of white sturgeon within 1 order of magnitude of measured values. These models provide a tool to better understand the impact of environmental contaminants on the health and the survival of endangered white sturgeon populations.

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.

Photocatalytic Penicillin Degradation Performance and the Mechanism of the Fragmented TiO2 Modified by CdS Quantum Dots

In this study, a novel method was adopted to construct a CdS-TiO₂ heterostructure to degrade penicillin under sunlight. A potato extract was used during the synthesis process of CdS QDs as a stabilizer and a modifier. The CdS-TiO₂ composite with a heterostructure delivers high photocatalytic degradation efficiency. In detail, 0.6 mg/mL of CdS-TiO₂ can successfully decompose penicillin after 2 h, and 5‰ CdS-TiO₂ shows the optimal degradation efficiency with the degradation rate reaching 88%. Furthermore, the underlying mechanisms of the penicillin decomposition reaction were investigated by the EPR test and trapping experiment. It was found that the high photocatalytic degradation efficiency was attributed to the heterojunction of CdS-TiO₂, which successfully suppresses the recombination of the conduction band of CdS and the valence band of TiO₂. Moreover, it was confirmed that the reaction is the O₂-consuming process, and introducing O₂ can greatly accelerate the generation of a superoxide radical during the photocatalytic degradation process, which eventually improves the degradation of penicillin and shortens the degradation time. Finally, this work provides the possible penicillin degradation pathways, which will inspire the researchers to explore and design novel photocatalysts in the field of wastewater treatment in the future.

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.

Zebrafish early life stages as alternative model to study 'designer drugs': Concordance with mammals in response to opioids

The number of new psychoactive substances (NPS) on the illicit drug market increases fast, posing a need to urgently understand their toxicity and behavioural effects. However, with currently available rodent models, NPS assessment is limited to a few substances per year. Therefore, zebrafish (Danio rerio) embryos and larvae have been suggested as an alternative model that would require less time and resources to perform an initial assessment and could help to prioritize substances for subsequent evaluation in rodents. To validate this model, more information on the concordance of zebrafish larvae and mammalian responses to specific classes of NPS is needed. Here, we studied toxicity and behavioural effects of opioids in zebrafish early life stages. Synthetic opioids are a class of NPS that are often used in pain medication but also frequently abused, having caused multiple intoxications and fatalities recently. Our data shows that fentanyl derivatives were the most toxic among the tested opioids, with toxicity in the zebrafish embryo toxicity test decreasing in the following order: butyrfentanyl>3-methylfentanyl>fentanyl>tramadol> O-desmethyltramadol>morphine. Similar to rodents, tramadol as well as fentanyl and its derivatives led to hypoactive behaviour in zebrafish larvae, with 3-methylfentanyl being the most potent. Physico-chemical properties-based predictions of chemicals' uptake into zebrafish embryos and larvae correlated well with the effects observed. Further, the biotransformation pattern of butyrfentanyl in zebrafish larvae was reminiscent of that in humans. Comparison of toxicity and behavioural responses to opioids in zebrafish and rodents supports zebrafish as a suitable alternative model for rapidly testing synthetic opioids.

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.

In Vitro Biotransformation Assays Using Liver S9 Fractions and Hepatocytes from Rainbow Trout (Oncorhynchus mykiss): Overcoming Challenges with Difficult to Test Fragrance Chemicals

In vitro metabolic stability assays using rainbow trout (Oncorhynchus mykiss) isolated hepatocytes (RT-HEP) or hepatic S9 fractions (RT-S9) were introduced to provide biotransformation rate data for the assessment of chemical bioaccumulation in fish. The present study explored the suitability of the RT-HEP and RT-S9 assays for difficult test chemicals, and the in vitro-based predictions were compared to in silico-based predictions and in vivo-measured bioconcentration factors (BCFs). The results show that volatile or reactive chemicals can be tested with minor modifications of the in vitro protocols. For hydrophobic chemicals, a passive dosing technique was developed. Finally, a design-of-experiment approach was used to identify optimal in vitro assay conditions. The modified assay protocols were applied to 10 fragrances with diverse physicochemical properties. The in vitro intrinsic clearance rates were higher in the S9 than in the hepatocyte assay, but the in vitro-in vivo (IVIV) predictions were comparable between the 2 assays. The IVIV predictions classified the test chemicals as nonbioaccumulative (BCF < 2000), which was in agreement with the in vivo data but in contrast to the in silico-based predictions. The findings from the present study provide strong evidence that the RT-HEP and RT-S9 assays can provide reliable estimates of in vivo biotransformation rates for test chemicals with difficult physicochemical properties. Environ Toxicol Chem 2020;39:2396-2408. © 2020 SETAC.

In vitro-in vivo and cross-life stage extrapolation of uptake and biotransformation of benzo[a]pyrene in the fathead minnow (Pimephales promelas)

Understanding internal dose metrics is integral to adequately assess effects environmental contaminants might have on aquatic wildlife, including fish. In silico toxicokinetic (TK) models are a leading approach for quantifying internal exposure metrics for fishes; however, they often do not adequately consider chemicals that are actively biotransformed and have not been validated against early-life stages (ELS) that are often considered the most sensitive to the exposure to contaminants. To address these uncertainties, TK models were parameterized for the rapidly biotransformed chemical benzo[a]pyrene (B[a]P) in embryo-larval and adult life stages of fathead minnows. Biotransformation of B[a]P was determined through measurements of in vitro clearance. Using in vitro-in vivo extrapolation, in vitro clearance was integrated into a multi-compartment TK model for adult fish and a one-compartment model for ELS. Model predictions were validated using measurements of B[a]P metabolites from in vivo flow-through exposures to graded concentrations of water-borne B[a]P. Significantly greater amounts of B[a]P metabolites were observed with exposure to greater concentrations of parent compound in both life stages. However, when assessing biotransformation capacity, no differences in phase I or phase II biotransformation were observed with greater exposures to B[a]P. Results of modelling suggested that biotransformation of B[a]P can be successfully implemented into in silico models to accurately predict life stage-specific abundances of B[a]P metabolites in either whole-body larvae or the bile of adult fish. Models developed increase the scope of applications in which TK models can be used to support environmental risk assessments.

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.

A Machine Learning Approach in Analyzing Bioaccumulation of Heavy Metals in Turbot Tissues

Metals are considered to be one of the most hazardous substances due to their potential for accumulation, magnification, persistence, and wide distribution in water, sediments, and aquatic organisms. Demersal fish species, such as turbot (Psetta maxima maeotica), are accepted by the scientific communities as suitable bioindicators of heavy metal pollution in the aquatic environment. The present study uses a machine learning approach, which is based on multiple linear and non-linear models, in order to effectively estimate the concentrations of heavy metals in both turbot muscle and liver tissues. For multiple linear regression (MLR) models, the stepwise method was used, while non-linear models were developed by applying random forest (RF) algorithm. The models were based on data that were provided from scientific literature, attributed to 11 heavy metals (As, Ca, Cd, Cu, Fe, K, Mg, Mn, Na, Ni, Zn) from both muscle and liver tissues of turbot exemplars. Significant MLR models were recorded for Ca, Fe, Mg, and Na in muscle tissue and K, Cu, Zn, and Na in turbot liver tissue. The non-linear tree-based RF prediction models (over 70% prediction accuracy) were identified for As, Cd, Cu, K, Mg, and Zn in muscle tissue and As, Ca, Cd, Mg, and Fe in turbot liver tissue. Both machine learning MLR and non-linear tree-based RF prediction models were identified to be suitable for predicting the heavy metal concentration from both turbot muscle and liver tissues. The models can be used for improving the knowledge and economic efficiency of linked heavy metals food safety and environment pollution studies.

Determination of the Bioaccumulative Potential Risk of Emerging Contaminants in Fish Muscle as an Environmental Quality Indicator in Coastal Lagoons of the Central Mexican Pacific

This research proposes an environmental quality indicator to detect, predict and scientifically evaluate the environmental impact generated by chemical substances within the pollutant group of nonsteroidal anti-inflammatory drugs (NSAIDs) that are categorized as emerging contaminants (ECs) with endocrine disruptive action. The present study was carried out in two coastal lagoons affected by wastewater produced by urban and rural settlements in the states of Colima and Jalisco. Four pharmaceutical compounds were analyzed: diclofenac, ibuprofen, ketorolac and naproxen. The muscle tissues of 14 fish species were analyzed; all had measurable concentrations of the four contaminants. The presence of the ECs was confirmed in all the samples collected. The bioaccumulative potential risk (BPR) of the ECs is estimated by calculating the environmental risk factors and the potential risks to human health, evaluating the concentration and assessing the risk involved in the incorporation of the pollutants into the environment. The BPR indicates the potential impact of NSAIDs on the ecology of these coastal lagoons, and predicts whether a contaminant is likely to act and persist in the environment and bioaccumulate in organisms. Additionally, it estimates the possibility of contamination and risks to human health caused by residues of the chemical contaminants.

Mechanisms of pH-Dependent Uptake of Ionizable Organic Chemicals by Fish from Oil Sands Process-Affected Water (OSPW)

Uptake and effects of ionizable organic chemicals (IOCs) that are weak acids in aqueous solution by fish can differ as a function of pH. While the pH-dependent behavior of select IOCs is well-understood, complex mixtures of IOCs, e.g., from oil sands process-affected water (OSPW), have not yet been studied systematically. Here, we established an in vitro screening method using the rainbow trout gill cell line, RTgill-W1, to investigate pH-dependent cytotoxicity and permeation of IOCs across cultured epithelia using ultra-high-performance liquid chromatography with high-resolution mass spectrometry (UPLC-HRMS). The assay was benchmarked using model chemicals and technical mixtures, and then used to characterize fractions and reconstituted extracts of field-collected OSPW. Significant pH-dependent cytotoxicity of individual IOCs, acidic fractions, and reconstituted extracts of OSPW was observed. In vitro data were in good agreement with data from a 96 h in vivo exposure experiment with juvenile rainbow trout. Permeation of some IOCs from OSPW was mediated by active transport, as revealed by studies in which inhibitors of these active transport mechanisms were applied. We conclude that the RTgill-W1 in vitro assay is useful for the screening of pH-dependent uptake of IOCs in fish, and has applications for in vitro-in vivo extrapolation, and prioritization of chemicals in nontarget screenings.

Tissue concentrations, trophic transfer and human risks of antibiotics in freshwater food web in Lake Taihu, China

The objective of this study was to evaluate the tissue distributions of antibiotics in the fish, the bioaccumulation and trophic transfer in freshwater food web in Taihu Lake, a large shallow freshwater lake. Twenty four out of 41 antibiotics were detected in the biotas of the food web; and antibiotic concentrations followed the orders: fish plasma ~ fish muscle < fish liver ~ fish bile and fish < invertebrates ~ plankton. Antibiotic concentrations in the liver of piscivores were higher than those in omnivores and planktivores. Most bioaccumulation factors (BAFs) of sulfonamides (SAs), macrolides (MLs), ionophores (IPs) and lincomycin (LIN) were less than 2000 L/kg, indicating low bioaccumulation ability of these compounds in fish. Fluoroquinolones (FQs) were frequently detected in fish liver, invertebrates and plankton with much of BAFs great than 5000 L/kg, indicating that FQs have the potential of bioaccumulation in fish. Relationship analysis between BAFs and physicochemical properties of antibiotics showed that the bioaccumulation of antibiotics in the biota was related with their adsorption ability. Generally, the antibiotics in the food web of Lake Taihu including plankton, invertebrates and fish showed trophic dilution. The normalized estimated daily intake (EDI) values are less than the acceptable daily intake (ADI) values, and then hazard quotients were much less than 1. This result suggests the consumption of fish, crab and shrimp in Lake Taihu would probably not pose direct detrimental effects on humans.

Toxicokinetics and bioaccumulation characteristics of bisphenol analogues in common carp (Cyprinus carpio)

Toxicokinetics and bioconcentration of eight common bisphenol analogues, including bisphenol A (BPA), -B, -C, -E, -S, -Z, -AF, and -AP in common carp (Cyprinus carpio) were investigated. Both free (BP_free) and total forms (BP_total) of the bisphenols were measured in various fish tissues. The conjugated forms of bisphenols were calculated based on BP_free and BP_total. The calculated bioconcentration factors (BCFs) based on the total bisphenols (BP_total) in the carp whole body were in the range of 0.3-320, agreeing with previous field results from Taihu Lake, China. The elimination rate constant (k_e) positively correlated with the fraction of conjugated form (f_conjugated), which displayed negative correlation with their log K_ow (r = -0.861, p < 0.05), indicating that conjugation facilitated their elimination and those with higher hydrophobicity were more difficult to be eliminated. Except BPA, the concentrations of all bisphenols in the carp tissues were in the order of kidney > liver ≫ muscle. The uptake rate constants (k_u) in kidney (r = 0.836, p < 0.05) and in liver (r = 0.863, p < 0.05) displayed significantly positive correlations with BCFs, and k_u in kidney was greater than in liver except BPA. These results indicated that kidney and liver played important roles in accumulating bisphenols in carp, and kidney made more contribution than liver for most bisphenols. Biliary excretion predominated for elimination of most bisphenols while BPA and BPS were mainly through urinary excretion.

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.

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.

Intestinal Fish Cell Barrier Model to Assess Transfer of Organic Chemicals in Vitro: An Experimental and Computational Study

We studied the role of the fish intestine as a barrier for organic chemicals using the epithelial barrier model built on the rainbow trout (Oncorhynchus mykiss) intestinal cell line, RTgutGC and the newly developed exposure chamber, TransFEr, specifically designed to work with hydrophobic and volatile chemicals. Testing 11 chemicals with a range of physicochemical properties (logK_OW: 2.2 to 6.3, logHLC: 6.1 to 2.3) and combining the data with a mechanistic kinetic model enabled the determination of dominant processes underlying the transfer experiments and the derivation of robust transfer rates. Against the current assumption in chemical uptake modeling, chemical transfer did not strictly depend on the logK_OW but resulted from chemical-specific intracellular accumulation and biotransformation combined with paracellular and active transport. Modeling also identified that conducting elaborate measurements of the plastic parts, including the polystyrene insert and the PET filter, is unnecessary and that stirring in the TransFEr chamber reduced the stagnant water layers compared to theoretical predictions. Aside from providing insights into chemical uptake via the intestinal epithelium, this system can easily be transferred to other cell-based barrier systems, such as the fish gill or mammalian intestinal models and may improve in vitro-in vivo extrapolation and prediction of chemical bioaccumulation into organisms.

Bioaccumulation, metabolism, and risk assessment of phenolic endocrine disrupting chemicals in specific tissues of wild fish

Phenolic endocrine disrupting chemicals (EDCs) may pose a great hazard to wildlife and humans, owing to their ubiquitous presence in the environment and potential bioaccumulation ability. We investigated the bioaccumulation, metabolism, and human health risks of six phenolic EDCs, including bisphenol A (BPA), 4-tert-octylphenol (4-t-OP), 4-nonylphenol (4-NP), estrone (E1), 17β-estradiol (E2), and 17α-ethinylestradiol (EE2), in wild fish from the Pearl River system, South China. Except EE2, the other five EDCs were detected in at least one of the four fish tissues (bile, liver, plasma, and muscle). The concentrations of BPA and 4-NP were greater than those of 4-t-OP, E1, and E2 in all tissues. The median values of log bioaccumulation factors for EDCs at the range of 3.86-4.52 in bile, 2.06-3.16 in liver, 2.69-3.87 in plasma, and 1.34-2.30 in muscle, indicating a higher bioaccumulation potential in fish bile than in other tissues. Greater levels of glucuronide/sulfate conjugated EDCs were found in fish bile and liver than in the plasma and muscle, suggesting that the liver and bile played an important role in the metabolism and excretion of phenolic EDCs in fish. The calculated hazard quotient values were below 1 for each compound, implying low risk to human health by intake of edible fish.

Improving a fish intestinal barrier model by combining two rainbow trout cell lines: epithelial RTgutGC and fibroblastic RTgutF

An in vitro model of the fish intestine is of interest for research and application in diverse fields such as fish physiology, aquaculture and chemical risk assessment. The recently developed epithelial barrier model of the fish intestine relies on the RTgutGC cell line from rainbow trout (Oncorhynchus mykiss), cultured in inserts on permeable membranes. Our aim was to extend the current system by introducing intestinal fibroblasts as supportive layer in order to reconstruct the epithelial-mesenchymal interface as found in vivo. We therefore initiated and characterized the first fibroblast cell line from the intestine of rainbow trout, which has been termed RTgutF. Co-culture studies of RTgutGC and RTgutF were performed on commercially available electric cell substrate for impedance sensing (ECIS) and on newly developed ultrathin, highly porous alumina membranes to imitate the cellular interaction with the basement membrane. Cellular events were examined with non-invasive impedance spectroscopy to distinguish between barrier tightness and cell density in the ECIS system and to determine transepithelial electrical resistance for cells cultured on the alumina membranes. We highlight the relevance of the piscine intestinal fibroblasts for an advanced intestinal barrier model, particularly on ultrathin alumina membranes. These membranes enable rapid crosstalk of cells cultured on opposite sides, which led to increased barrier tightening in the fish cell line-based epithelial-mesenchymal model.

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.

Influence of in Vitro Assay Setup on the Apparent Cytotoxic Potency of Benzalkonium Chlorides

The nominal concentration is generally used to express concentration–effect relationships in in vitro toxicity assays. However, the nominal concentration does not necessarily represent the exposure concentration responsible for the observed effect. Surfactants accumulate at interphases and likely sorb to in vitro system components such as serum protein and well plate plastic. The extent of sorption and the consequences of this sorption on in vitro readouts is largely unknown for these chemicals. The aim of this study was to demonstrate the effect of sorption to in vitro components on the observed cytotoxic potency of benzalkonium chlorides (BAC) varying in alkyl chain length (6–18 carbon atoms, C_6–18) in a basal cytotoxicity assay with the rainbow trout gill cell line (RTgill-W1). Cells were exposed for 48 h in 96-well plates to increasing concentration of BACs in exposure medium containing 0, 60 μM bovine serum albumin (BSA) or 10% fetal bovine serum (FBS). Before and after exposure, BAC concentrations in exposure medium were analytically determined. Based on freely dissolved concentrations at the end of the exposure, median effect concentrations (EC₅₀) decreased with increasing alkyl chain length up to 14 carbons. For BAC with alkyl chains of 12 or more carbons, EC₅₀’s based on measured concentrations after exposure in supplement-free medium were up to 25-times lower than EC₅₀’s calculated using nominal concentrations. When BSA or FBS was added to the medium, a decrease in cytotoxic potency of up to 22 times was observed for BAC with alkyl chains of eight or more carbons. The results of this study emphasize the importance of expressing the in vitro readouts as a function of a dose metric that is least influenced by assay setup to compare assay sensitivities and chemical potencies.

Bioaccumulation and risks of 24 personal care products in plasma of wild fish from the Yangtze River, China

We used a hybrid precipitation method to simultaneously extract and analyze 24 personal care products (PCPs), including 16 biocides, 4 synthetic musks, and 4 benzotriazoles, in the plasma of fish. The method's performance was validated for plasma samples with and without β-glucuronidase/aryl-sulfatase hydrolysis. The recoveries were in the range of 70-120% for most of the PCPs, except N,N-diethyl-3-methylbenzamide (DEET), clotrimazole (CTZ), miconazole and itraconazole at spiking concentration of 20 and 5 ng/mL. The quantification limits ranged between 0.89 and 17.9 ng/mL (hydrolyzed plasma) and 0.85-18.5 ng/mL (non-hydrolyzed plasma), except CTZ at 77.5 ng/mL and 76.3 ng/mL. Totally, 13 PCPs were detected in plasma samples of fish collected from the Yangtze River, with a maximum concentration of 58.4 ng/mL (galaxolide). Compounds with the phenol hydroxyl groups of parabens or triclosan in hydrolyzed plasma showed higher concentrations than those in unhydrolyzed plasma with the ratio of conjugation (glucuronides + sulfates) forms up to 86%. The median values for the logarithm of bioaccumulation factors were between 1.39 and 4.15, which were 2-3 orders of magnitude higher than the theoretical logarithm of bioconcentration factors. Using the fish plasma model, the effect ratios (effect concentration/measured plasma concentration ratios) of tonalide, galaxolide, benzotriazole, triclosan, and DEET reached 0.35, 4.15, 3.78, 7.52, and 9.24, respectively. These are recognized as priority chemicals for further risk assessment.

A Petri Net Approach to Physiologically Based Toxicokinetic Modeling

Physiologically based toxicokinetic (PBTK) modeling enables researchers to predict internal tissue concentrations for various species exposed to exogenous compounds through different routes at varying concentrations without having to run in vivo experiments for each scenario. Parameters for the models may be gathered from in vivo or in vitro measurements, cross-species or cross-chemical extrapolations, literature reviews, or other models. The PBTK models, described using ordinary differential equations (ODEs), are then simulated using these parameters for a given compound/exposure/species scenario. Although they are potentially useful for regulatory toxicology, the complexity of ODE programming and simulation remains a barrier for many would-be researchers. Petri nets, a graphical modeling framework, offers a more intuitive approach to PBTK modeling. To demonstrate their utility and ease of use, we present a model of waterborne fluoranthene exposure to rainbow trout (Oncorhynchus mykiss) written and simulated in Snoopy, a graphical Petri net development and simulation software package. We converted an existing ODE PBTK model and evaluated the Petri net model against the ODE model results. The simulated tissue concentrations of the Petri net model closely mirrored the simulated concentrations of the ODE model. To convert the ODE model to a Petri net model, we introduced a new parameter, blood volume (V _BLOOD ). Sensitivity analysis found V _BLOOD to be very robust when varied over an order of magnitude. The resulting Petri net PBTK model has a number of advantages over ODE models, while maintaining equivalent predictive functionality. Environ Toxicol Chem 2019;00:1-10. © 2019 SETAC.

Evaluation of toxicokinetics of nonylphenol in the adult female Sprague-Dawley rats using a physiologically based toxicokinetic model

The physiologically based toxicokinetic (PBTK) model was firstly developed in female rats to quantitatively evaluate toxicokinetics of nonylphenol (NP). Changes in NP serum concentrations over time of single oral NP administration experiments in Sprague-Dawley rats and literature data were collected to establish and calibrate the PBTK model in the SimBiology framework. The calibrated model predicted the serum and tissue NP concentrations of repeat oral NP administration for model evaluation. NP concentrations in serum and tissues (liver, kidneys, adipose, brain, uterus and ovaries) were quantified using ultra-high performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). The model output of the time course data (values are within the standard deviation defined for each data point) indicated proximity of predictions to reality. The coefficients of determination r² were all greater than 0.9, and the root mean squared error (RMSE) were within 0.177-2.027, which means the model predicted and observed serum NP concentrations were in excellent agreement. The results indicated that the model could contribute to a simplification of the future exposure risk assessments of NP in a more realistic scenario and provided a better understanding of the disposition process of NP in human.

Pregnant Inuit Women's Exposure to Metals and Association with Fetal Growth Outcomes: ACCEPT 2010⁻2015

Environmental contaminants such as heavy metals are transported to the Arctic regions via atmospheric and ocean currents and enter the Arctic food web. Exposure is an important risk factor for health and can lead to increased risk of a variety of diseases. This study investigated the association between pregnant women's levels of heavy and essential metals and the birth outcomes of the newborn child. This cross-sectional study is part of the ACCEPT birth cohort (Adaption to Climate Change, Environmental Pollution, and dietary Transition) and included 509 pregnant Inuit women ≥18 years of age. Data were collected in five Greenlandic regions during 2010⁻2015. Population characteristics and birth outcomes were obtained from medical records and midwives, respectively, and blood samples were analyzed for 13 metals. Statistical analysis included one-way ANOVA, Spearman's rho, and multiple linear and logistic regression analyses. The proportion of current smokers was 35.8%. The levels of cadmium, chromium, and nickel were higher compared to reported normal ranges. Significant regional differences were observed for several metals, smoking, and parity. Cadmium and copper were significantly inversely related to birth outcomes. Heavy metals in maternal blood can adversely influence fetal development and growth in a dose⁻response relationship. Diet and lifestyle factors are important sources of toxic heavy metals and deviant levels of essential metals. The high frequency of smokers in early pregnancy is of concern, and prenatal exposure to heavy metals and other environmental contaminants in the Greenlandic Inuit needs further research.

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.

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.

Substituted diphenylamine antioxidants and benzotriazole UV stabilizers in blood plasma of fish, turtles, birds and dolphins from North America

Substituted diphenylamine antioxidants (SDPAs) and benzotriazole UV stabilizers (BZT-UVs) are additives used in industrial and commercial applications to prevent degradation by oxidation and are contaminants of emerging environmental concern. Little is known about the fate of these contaminants in wildlife, particularly in reptiles, birds and marine mammals. Nine SDPAs and six BZT-UVs were measured in blood plasma of seven fish species, snapping turtles (Chelydra serpentina), double-crested cormorants (Phalacrocorax auritus), and bottlenose dolphins (Tursiops truncatus) from various locations in North America. Plasma SDPAs were more frequently (90-100%) detected and with higher concentrations (median: 25-270 pg g^-1, wet weight (ww)) in organisms from urban areas than rural locations (median: <method limits of quantification -136 pg g^-1). The concentrations of most SDPAs generally followed the order of fish ≥ snapping turtles > double-crested cormorants > bottlenose dolphins. Of the three quantifiable BZT-UVs, 2-(2H-benzotriazol-2-yl)-4,6-di-tert-pentylphenol (UV328) showed higher detection frequency in most species of fish, bird and turtle (range of 0-67%), indicating the widespread distribution of UV328 in the aquatic environment of lower Great Lakes region.

Measurement of kinetic parameters for biotransformation of polycyclic aromatic hydrocarbons by trout liver S9 fractions: Implications for bioaccumulation assessment

In vitro substrate depletion methods developed by the pharmaceutical industry are being used with increasing frequency to support chemical bioaccumulation assessments for fish. However, the application of these methods to high log K _ow chemicals poses special challenges. Biotransformation of three polycyclic aromatic hydrocarbons (PAHs) was measured using trout liver S9 fractions. Measured activity declined with incubation time and was reduced by acetone (used as a spiking solvent) at concentrations greater than 0.5%. Addition of alamethicin, a pore-forming peptide used to support UDP-glucuronosyltransferase activity, also reduced activity in a concentration-dependent manner. The substrate concentration dependence of activity was evaluated to estimate K _M and V _max values for each compound. Derived kinetic constants suggested that all three PAHs are transformed by the same reaction pathway and indicated an inverse correlation between K _M and chemical log K _ow. Binding effects on activity were evaluated by measuring unbound chemical concentrations across a range of S9 protein levels. Reaction rates were proportional to the unbound concentration except when these concentrations approached saturating levels, providing a direct demonstration of the free chemical hypothesis. These findings suggest that previous in vitro work with high log K _ow compounds was conducted at inappropriately high substrate concentrations resulting in underestimation of true in vivo activity. Preliminary calculations also indicate that PAH metabolism in fish may approach saturation during standardized in vivo testing efforts, potentially resulting in concentration-dependent accumulation and/or steady-state levels of accumulation greater than those which occur in a natural setting.

Uptake, distribution and elimination of chemicals in fish - Which physiological parameters are the most relevant for toxicokinetics?

Bioconcentration and toxicity studies are regularly conducted for the risk assessment of chemicals. If such tests yield different results for different fish species, this can either be due to differences in toxicokinetics or to differences in toxicodynamics. Here we investigate which physiological parameters could cause major differences in the toxicokinetics in fish. To this end it is important to distinguish physiological parameters that affect the sorption capacity of the fish from those that affect kinetic processes. Variability in the lipid content of a fish is the most influential parameter for the sorption capacity of fish and therefore most relevant for the total concentration in fish under steady-state conditions when metabolism is not relevant. In terms of kinetics, ventilation rate, uptake efficiency from food and metabolism are the most influential factors. While ventilation rate can roughly be estimated from allometric scaling equations, little general information is available on the uptake efficiency from food. The metabolism rate constant appears to be the single most influential toxicokinetic factor. This information cannot be estimated but must be determined experimentally, preferably from in vitro experiments.

Multi-compartmental toxicokinetic modeling of fipronil in tilapia: Accumulation, biotransformation and elimination

Bioaccumulation and biotransformation are critical processes modifying toxicity of easily metabolizable chemicals to aquatic organisms. In this study, tissue-specific accumulation, biotransformation and elimination of a current-use pesticide fipronil in tilapia (Oreochromis niloticus) were quantified by combining in vivo measurements and a newly developed multi-compartmental toxicokinetic model. Waterborne fipronil was taken up via gills and metabolized rapidly and solely to fipronil sulfone. Significant decrease of fipronil residues in liver and intestine during exposure period strongly suggested the induction of metabolism in these two organs. Significant transport of fipronil and fipronil sulfone in the liver-bile-intestine system implied that hepatobiliary excretion and enterohepatic re-absorption played important roles in fipronil metabolism and system circulation of the parent compound and the metabolite. The multi-compartmental model quantitatively described the highly dynamic inter-compartmental transport and rapid branchial clearance of fipronil in fish. Modeling results also suggested that uptake and biotransformation were the stronger driving forces for the inter-compartmental transport of fipronil in fish than the inherent partitioning capacity. Overall, our findings highlight the importance of biotransformation on internal disposition of fipronil in fish, which helps to improve aquatic toxicity assessment of this pesticide.