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

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.

Time Is Ripe for Targeting Per- and Polyfluoroalkyl Substances-Induced Hormesis: Global Aquatic Hotspots and Implications for Ecological Risk Assessment

Globally implemented ecological risk assessment (ERA) guidelines marginalize hormesis, a biphasic dose-response relationship characterized by low-dose stimulation and high-dose inhibition. The present study illuminated the promise of hormesis as a scientific dose-response model for ERA of per- and polyfluoroalkyl substances (PFAS) represented by perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS). A total of 266 hormetic dose-response relationships were recompiled from 1237 observations, covering 30 species from nine representative taxonomic groups. The standardized hormetic amplitudes followed the log-normal probability distribution, being subject to the limits of biological plasticity but independent of stress inducers. The SHapley Additive exPlanations algorithm revealed that the target endpoint was the most important variable explaining the hormetic amplitudes. Subsequently, quantitative frameworks were established to incorporate hormesis into the predicted no-effect concentration levels, with a lower induction dose and a zero-equivalent point but a broader hormetic zone for PFOS. Realistically, 10,117 observed concentrations of PFOA and PFOS were gathered worldwide, 4% of which fell within hormetic zones, highlighting the environmental relevance of hormesis. Additionally, the hormesis induction potential was identified in other legacy and emerging PFAS as well as their alternatives and mixtures. Collectively, it is time to incorporate the hormesis concept into PFAS studies to facilitate more realistic risk characterizations.

Interspecies scaling of toxicity reference values in human health versus ecological risk assessments: A critical review

Risk assessments that focus on anthropogenic chemicals in environmental media-whether considering human health or ecological effects-often rely on toxicity data from experimentally studied species to estimate safe exposures for species that lack similar data. Current default extrapolation approaches used in both human health risk assessments and ecological risk assessments (ERAs) account for differences in body weight between the test organisms and the species of interest, but the two default approaches differ in important ways. Human health risk assessments currently employ a default based on body weight raised to the three-quarters power. Ecological risk assessments for wildlife (i.e., mammals and birds) are typically based directly on body weight, as measured in the test organism and receptor species. This review describes differences in the experimental data underlying these default practices and discusses the many factors that affect interspecies variability in chemical exposures. The interplay of these different factors can lead to substantial departures from default expectations. Alternative methodologies for conducting more accurate interspecies extrapolations in ERAs for wildlife are discussed, including tissue-based toxicity reference values, physiologically based toxicokinetic and/or toxicodynamic modeling, chemical read-across, and a system of categorical defaults based on route of exposure and toxic mode of action. Integr Environ Assess Manag 2024;20:749-764. © 2023 SETAC.

Tissue-specific distributions of organic ultraviolet absorbents in free-range chickens and domestic pigeons from Guangzhou, China

The ecological risks of organic ultraviolet absorbents (UVAs) have been of increasing concern. Studies have found that these chemicals could be accumulated in terrestrial animals and pose toxicities. However, tissue distribution of UVAs in terrestrial species was far from well understood. In this study, free-range chickens and domestic pigeons were selected to investigate the occurrence and tissue distribution of UVAs. Total concentrations of eleven UVAs in muscles ranged from 778 to 2874 (mean 1413 ± 666) ng/g lipid weight, which were higher than those in aquatic species worldwide. Since low UVA concentrations in local environment were previously reported, the results implied the strong accumulation of UVAs in studied species. Brain, stomach and kidney were main target organs for studied UVAs, differentiating from the strong liver sequestration in aquatic species. The mean tissue-to-muscle ratios of 1.02-4.23 further indicated the preferential accumulation of target UVAs in these tissues. The tissue-to-blood ratios of benzophenone (BP), 2-ethylhexyl salicylate (EHS) and homosalate (HMS) in brain were 370, 1207 and 52.0, respectively, implying their preferential accumulation in brain. More research is needed to characterize the toxicokinetics and tissue distribution of UVAs in terrestrial wild species, in order to further understand their potential risks.

Review: Fish bile, a highly versatile biomarker for different environmental pollutants

Ecotoxicological assessments encompass a broad spectrum of biochemical endpoints and ecological factors, allowing for comprehensive assessments concerning pollutant exposure levels and their effects on both fish populations and surrounding ecosystems. While these evaluations offer invaluable insights into the overall health and dynamics of aquatic environments, they often provide an integrated perspective, making it challenging to pinpoint the precise sources and individual-level responses to environmental contaminants. In contrast, biliary pollutant excretion assessments represent a focused approach aimed at understanding how fish at the individual level respond to environmental stressors. In this sense, the analysis of pollutant profiles in fish bile not only serves as a valuable exposure indicator, but also provides critical information concerning the uptake, metabolism, and elimination of specific contaminants. Therefore, by investigating unique and dynamic fish responses to various pollutants, biliary assessments can contribute significantly to the refinement of ecotoxicological studies. This review aims to discuss the multifaceted utility of bile as a potent biomarker for various environmental pollutants in fish in targeted monitoring strategies, such as polycyclic aromatic hydrocarbons, metals, pesticides, pharmaceuticals, estrogenic compounds, resin acids, hepatotoxins and per- and polyfluorinated substances. The main caveats of this type of assessment are also discussed, as well as future directions of fish bile studies.

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.

The posterity of Zebrafish in paradigm of in vivo molecular toxicological profiling

The aggrandised advancement in utility of advanced day-to-day materials and nanomaterials has raised serious concern on their biocompatibility with human and other biotic members. In last few decades, understanding of toxicity of these materials has been given the centre stage of research using many in vitro and in vivo models. Zebrafish (Danio rerio), a freshwater fish and a member of the minnow family has garnered much attention due to its distinct features, which make it an important and frequently used animal model in various fields of embryology and toxicological studies. Given that fertilization and development of zebrafish eggs take place externally, they serve as an excellent model organism for studying early developmental stages. Moreover, zebrafish possess a comparable genetic composition to humans and share almost 70% of their genes with mammals. This particular model organism has become increasingly popular, especially for developmental research. Moreover, it serves as a link between in vitro studies and in vivo analysis in mammals. It is an appealing choice for vertebrate research, when employing high-throughput methods, due to their small size, swift development, and relatively affordable laboratory setup. This small vertebrate has enhanced comprehension of pathobiology and drug toxicity. This review emphasizes on the recent developments in toxicity screening and assays, and the new insights gained about the toxicity of drugs through these assays. Specifically, the cardio, neural, and, hepatic toxicology studies inferred by applications of nanoparticles have been highlighted.

Multi-tissue proteogenomic analysis for mechanistic toxicology studies in non-model species

New approach methodologies (NAM), including omics and in vitro approaches, are contributing to the implementation of 3R (reduction, refinement and replacement) strategies in regulatory science and risk assessment. In this study, we present an integrative transcriptomics and proteomics analysis workflow for the validation and revision of complex fish genomes and demonstrate how proteogenomics expression matrices can be used to support multi-level omics data integration in non-model species in vivo and in vitro. Using Atlantic salmon as an example, we constructed proteogenomic databases from publicly available transcriptomic data and in-house generated RNA-Seq and LC-MS/MS data. Our analysis identified ∼80,000 peptides, providing direct evidence of translation for over 40,000 RefSeq structures. The data also highlighted 183 co-located peptide groups that supported a single transcript each, and in each case, either corrected a previous annotation, supported Ensembl annotations not present in RefSeq, or identified novel previously unannotated genes. Proteogenomics data-derived expression matrices revealed distinct profiles for the different tissue types analyzed. Focusing on proteins involved in defense against xenobiotics, we detected distinct expression patterns across different salmon tissues and observed homology in the expression of chemical defense proteins between in vivo and in vitro liver systems. Our study demonstrates the potential of proteogenomic analyses in extending our understanding of complex fish genomes and provides an advanced bioinformatic toolkit to support the further development of NAMs and their application in regulatory science and (eco)toxicological studies of non-model species.

Physiology-informed toxicokinetic model for the zebrafish embryo test developed for bisphenols

Zebrafish embryos (ZFE) is a widely used model organism, employed in various research fields including toxicology to assess e.g., developmental toxicity and endocrine disruption. Variation in effects between chemicals are difficult to compare using nominal dose as toxicokinetic properties may vary. Toxicokinetic (TK) modeling is a means to estimate internal exposure concentration or dose at target and to enable extrapolation between experimental conditions and species, thereby improving hazard assessment of potential pollutants. In this study we advance currently existing TK models for ZFE with physiological ZFE parameters and novel experimental bisphenol data, a class of chemicals with suspected endocrine activity. We developed a five-compartment model consisting of water, plastic, chorion, yolk sack and embryo in which surface area and volume changes as well as the processes of biotransformation and blood circulation influence mass fluxes. For model training and validation, we measured internal concentrations in ZFE exposed individually to BPA, bisphenol AF (BPAF) and Z (BPZ). Bayesian inference was applied for parameter calibration based on the training data set of BPZ. The calibrated TK model predicted internal ZFE concentrations of the majority of external test data within a 5-fold error and half of the data within a 2-fold error for bisphenols A, AF, F, and tetrabromo bisphenol A (TBBPA). We used the developed model to rank the hazard of seven bisphenols based on predicted internal concentrations and measured in vitro estrogenicity. This ranking indicated a higher hazard for BPAF, BPZ, bisphenol B and C (BPB, BPC) than for BPA.

TKPlate 1.0: An Open-access platform for toxicokinetic and toxicodynamic modelling of chemicals to implement new approach methodologies in chemical risk assessment

This publication is linked to the following EFSA Supporting Publications articles: http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2023.EN-8441/full, http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2023.EN-8440/full, http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2023.EN-8437/full.

A p21-GFP zebrafish model of senescence for rapid testing of senolytics in vivo

Senescence drives the onset and severity of multiple ageing-associated diseases and frailty. As a result, there has been an increased interest in mechanistic studies and in the search for compounds targeting senescent cells, known as senolytics. Mammalian models are commonly used to test senolytics and generate functional and toxicity data at the level of organs and systems, yet this is expensive and time consuming. Zebrafish share high homology in genes associated with human ageing and disease. They can be genetically modified relatively easily. In larvae, most organs develop within 5 days of fertilisation and are transparent, which allows tracking of fluorescent cells in vivo in real time, testing drug off-target toxicity and assessment of cellular and phenotypic changes. Here, we have generated a transgenic zebrafish line that expresses green fluorescent protein (GFP) under the promoter of a key senescence marker, p21. We show an increase in p21:GFP+ cells in larvae following exposure to ionising radiation and with natural ageing. p21:GFP+ cells display other markers of senescence, including senescence-associated β-galactosidase and IL6. The observed increase in senescent cells following irradiation is associated with a reduction in the thickness of muscle fibres and mobility, two important ageing phenotypes. We also show that quercetin and dasatinib, two senolytics currently in clinical trials, reduce the number of p21:GFP+ cells, in a rapid 5-day assay. This model provides an important tool to study senescence in a living organism, allowing the rapid selection of senolytics before moving to more expensive and time-consuming mammalian systems.

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.

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 meta-analysis of ecotoxicological models used for plant protection product risk assessment before their placing on the market

Before their placing on the market, the safety of plant protection products (PPP) towards both human and animal health, and the environment has to be assessed using experimental and modelling approaches. Models are crucial tools for PPP risk assessment and some even help to avoid animal testing. This review investigated the use of modelling approaches in the ecotoxicology section of PPP active substance assessment reports prepared by the authorities and opened to consultation from 2011 to 2021 in the European Union. Seven categories of models (Structure-Activity, ToxicoKinetic, ToxicoKinetic-ToxicoDynamic, Species Sensitivity Distribution, population, community, and mixture) were searched for into the reports of 317 active substances. At least one model category was found for 44 % of the investigated active substances. The most detected models were Species Sensitivity Distribution, Structure-Activity and ToxicoKinetic for 27, 21 and 15 % of the active substances, respectively. The use of modelling was of particular importance for conventional active substances such as sulfonylurea or carbamates contrary to microorganisms and plant derived substances. This review also highlighted a strong imbalance in model usage among the biological groups considered in the European Regulation (EC) No 1107/2009. For example, models were more often used for aquatic than for terrestrial organisms (e.g., birds, mammals). Finally, a gap between the set of models used in reports and those existing in the literature was observed highlighting the need for the implementation of more sophisticated models into PPP regulation.

Organ-specific accumulation of cadmium and zinc in Gammarus fossarum exposed to environmentally relevant metal concentrations

One of the best approaches for improving the assessment of metal toxicity in aquatic organisms is to study their organotropism (i.e., the distribution of metals among organs) through a dynamical approach (i.e., via kinetic experiments of metal bioaccumulation), to identify the tissues/organs that play a key role in metal regulation (e.g., storage or excretion). This study aims at comparing the organ-specific metal accumulation of a non-essential (Cd) and an essential metal (Zn), at their environmentally relevant exposure concentrations, in the gammarid Gammarus fossarum. Gammarids were exposed for 7 days to ¹⁰⁹Cd- or ⁶⁵Zn-radiolabeled water at a concentration of 52.1 and 416 ng.L^-1 (stable equivalent), respectively, and then placed in clean water for 21 days. At different time intervals, the target organs (i.e., caeca, cephalons, intestines, gills, and remaining tissues) were collected and ¹⁰⁹Cd or ⁶⁵Zn contents were quantified by gamma-spectrometry. A one-compartment toxicokinetic (TK) model was fitted by Bayesian inference to each organ/metal dataset in order to establish TK parameters. Our results indicate: i) a contrasting distribution pattern of concentrations at the end of the accumulation phase (7^th day): gills > caeca ≈ intestines > cephalons > remaining tissues for Cd and intestines > caeca > gills > cephalons > remaining tissues for Zn; ii) a slower elimination of Cd than of Zn by all organs, especially in the gills in which the Cd concentration remained constant during the 21-day depuration phase, whereas Zn concentrations decreased sharply in all organs after 24 h in the depuration phase; iii) a major role of intestines in the uptake of waterborne Cd and Zn at environmentally relevant concentrations.

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.

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.

Liver-Based Probabilistic Risk Assessment of Exposure to Organophosphate Esters via Dust Ingestion Using a Physiologically Based Toxicokinetic (PBTK) Model

Organophosphate esters (OPEs) are widely used and harmful to organisms and human health. Dust ingestion is an important exposure route for OPEs to humans. In this study, by integrating ToxCast high-throughput in vitro assays with in vitro to in vivo extrapolation (IVIVE) via physiologically based Toxicokinetic (PBTK) modeling, we assessed the hepatocyte-based health risk for humans around the world due to exposure to two typical OPEs (TPHP and TDCPP) through the dust ingestion exposure route. Results showed that the health guidance value of TPHP and TCDPP obtained in this study was lower than the value obtained through animal experiments. In addition, probabilistic risk assessment results indicate that populations worldwide are at low risk of exposure to TPHP and TDCPP through dust ingestion due to low estimated daily intakes (EDIs) which are much lower than the reference dose (RfDs) published by the US EPA, except in some regional cases. Most margin of exposure (MOE) ranges of TDCPP for children are less than 100, which indicates a moderately high risk. Researchers should be concerned about exposure to TDCPP in this area. The method proposed in this study is expected to be applied to the health risk assessment of other chemicals.

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.

The role of Dynamic Energy Budgets in conservation physiology

The contribution of knowledge, concepts and perspectives from physiological ecology to conservation decision-making has become critical for understanding and acting upon threats to the persistence of sensitive species. Here we review applications of dynamic energy budget (DEB) theory to conservation issues and discuss how this theory for metabolic organization of all life on earth (from bacteria to whales) is well equipped to support current and future investigations in conservation research. DEB theory was first invented in 1979 in an applied institution for environmental quality assessment and mitigation. The theory has since undergone extensive development and applications. An increasing number of studies using DEB modelling have provided valuable insights and predictions in areas that pertain to conservation such as species distribution, evolutionary biology, toxicological impacts and ecosystem management. We discuss why DEB theory, through its mechanistic nature, its universality and the wide range of outcomes it can provide represents a valuable tool to tackle some of the current and future challenges linked to maintaining biodiversity, ensuring species survival, ecotoxicology, setting water and soil quality standards and restoring ecosystem structure and functioning in a changing environment under the pressure of anthropogenic driven changes.

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 population approach for the estimation of methylmercury ToxicoKinetics in red mullets

It is known that, as the vast majority of the anthropogenically emitted mercury can be found in aquatic ecosystems, where several methylating bacteria are present, fish consumption represents the most critical intake source of the most toxic form of mercury, the methylmercury (MeHg). The aim of this work is to predict MeHg levels in the fish muscles which, being the edible portion, are part of the human diet. A physiologically based toxicokinetics model was developed to evaluate the kinetics of MeHg in red mullets. Fishes were described by means of a multi-compartment model including stomach, gut, blood, muscles and an additional compartment virtually encompassing all the remaining organs. Absorption, distribution and excretion were modelled considering different MeHg routes of administration and excretion: intake by ingestion of contaminated food, intake and elimination through inhalation-exhalation and excretion through feces. The model has been firstly validated on Terapon jarbua fish (using the weighted least squares method for parameter estimation) to be subsequently readapted to predict methylmercury concentrations in the muscle of red mullets (using an approximate Bayesian computation approach). This simple multicompartmental model could be considered part, a link in the chain, of a wider more complex project aiming at tracking the fate of MeHg from polluted seawater to the human end consumer. The present study could be useful to surveillance organizations in order to carry out a more comprehensive and informed risk assessment analysis and to take appropriate preventive measures by evaluating possible new MeHg concentration thresholds to minimize public health hazards.

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.

Effect of Water pH on the Uptake of Acidic (Ibuprofen) and Basic (Propranolol) Drugs in a Fish Gill Cell Culture Model

Water pH is predicted to affect the uptake of ionizable pharmaceuticals in fish. The current study used an in vitro primary fish gill cell culture system to assess the effect of pH values in the range of 4.5-8.75 on the uptake rates of the base propranolol (pKa 9.42) and the acid ibuprofen (pKa 4.59). The rate-limiting step in the uptake was the diffusive supply flux of the unionized form from the water to the apical membrane, with subsequent rapid transfer across the epithelium. Computed uptake rate based on the unionized fraction best described the uptake of propranolol and ibuprofen over the range of pH values 5-8 and 6-8.75, respectively. For ibuprofen, the computed uptake rate overestimated the uptake below pH 6 where the unionized fraction increased from 4% at pH 6 to 55% at pH 4.5. As the unionized fraction increased, the uptake rate plateaued suggesting a saturation of the transport process. For both drugs, large variations in the uptake occur with only small fluctuations in pH values. This occurs between pH values 6 and 8, which is the pH range acceptable in regulatory test guidelines and seen in most of our freshwaters.

Muscle Glycogen Phosphorylase and Its Functional Partners in Health and Disease

Glycogen phosphorylase (PG) is a key enzyme taking part in the first step of glycogenolysis. Muscle glycogen phosphorylase (PYGM) differs from other PG isoforms in expression pattern and biochemical properties. The main role of PYGM is providing sufficient energy for muscle contraction. However, it is expressed in tissues other than muscle, such as the brain, lymphoid tissues, and blood. PYGM is important not only in glycogen metabolism, but also in such diverse processes as the insulin and glucagon signaling pathway, insulin resistance, necroptosis, immune response, and phototransduction. PYGM is implicated in several pathological states, such as muscle glycogen phosphorylase deficiency (McArdle disease), schizophrenia, and cancer. Here we attempt to analyze the available data regarding the protein partners of PYGM to shed light on its possible interactions and functions. We also underline the potential for zebrafish to become a convenient and applicable model to study PYGM functions, especially because of its unique features that can complement data obtained from other approaches.

A multi-taxonomic framework for assessing relative petrochemical vulnerability of marine biodiversity in the Gulf of Mexico

A fundamental understanding of the impact of petrochemicals and other stressors on marine biodiversity is critical for effective management, restoration, recovery, and mitigation initiatives. As species-specific information on levels of petrochemical exposure and toxicological response are lacking for the majority of marine species, a trait-based assessment to rank species vulnerabilities to petrochemical activities in the Gulf of Mexico can provide a more comprehensive and effective means to prioritize species, habitats, and ecosystems for improved management, restoration and recovery. To initiate and standardize this process, we developed a trait-based framework, applicable to a wide range of vertebrate and invertebrate species, that can be used to rank relative population vulnerabilities of species to petrochemical activities in the Gulf of Mexico. Through expert consultation, 18 traits related to likelihood of exposure, individual sensitivity, and population resilience were identified and defined. The resulting multi-taxonomic petrochemical vulnerability framework can be adapted and applied to a wide variety of species groups and geographic regions. Additional recommendations and guidance on the application of the framework to rank species vulnerabilities under specific petrochemical exposure scenarios, management needs or data limitations are also discussed.

The Zebrafish Xenograft Models for Investigating Cancer and Cancer Therapeutics

Simple Summary

The identification and development of new anti-cancer drugs requires extensive testing in animal models to establish safety and efficacy of drug candidates. The transplantation of human tumor tissue into mouse (tumor xenografts) is commonly used to study cancer progression and to test potential drugs for their anti-cancer activity. Mouse models do not afford the ability to test a large number of drug candidates quickly as it takes several weeks to conduct these experiments. In contrast, tumor xenograft studies in zebrafish provide an efficient platform for rapid testing of safety and efficacy in less than two weeks.

Abstract

In order to develop new cancer therapeutics, rapid, reliable, and relevant biological models are required to screen and validate drug candidates for both efficacy and safety. In recent years, the zebrafish (Danio rerio) has emerged as an excellent model organism suited for these goals. Larval fish or immunocompromised adult fish are used to engraft human cancer cells and serve as a platform for screening potential drug candidates. With zebrafish sharing ~80% of disease-related orthologous genes with humans, they provide a low cost, high-throughput alternative to mouse xenografts that is relevant to human biology. In this review, we provide background on the methods and utility of zebrafish xenograft models in cancer research.

EFSA's OpenFoodTox: An open source toxicological database on chemicals in food and feed and its future developments

Since its creation in 2002, the European Food Safety Authority (EFSA) has produced risk assessments for over 5000 substances in >2000 Scientific Opinions, Statements and Conclusions through the work of its Scientific Panels, Units and Scientific Committee. OpenFoodTox is an open source toxicological database, available both for download and data visualisation which provides data for all substances evaluated by EFSA including substance characterisation, links to EFSA's outputs, applicable legislations regulations, and a summary of hazard identification and hazard characterisation data for human health, animal health and ecological assessments. The database has been structured using OECD harmonised templates for reporting chemical test summaries (OHTs) to facilitate data sharing with stakeholders with an interest in chemical risk assessment, such as sister agencies, international scientific advisory bodies, and others. This manuscript provides a description of OpenFoodTox including data model, content and tools to download and search the database. Examples of applications of OpenFoodTox in chemical risk assessment are discussed including new quantitative structure-activity relationship (QSAR) models, integration into tools (OECD QSAR Toolbox and AMBIT-2.0), assessment of environmental footprints and testing of threshold of toxicological concern (TTC) values for food related compounds. Finally, future developments for OpenFoodTox 2.0 include the integration of new properties, such as physico-chemical properties, exposure data, toxicokinetic information; and the future integration within in silico modelling platforms such as QSAR models and physiologically-based kinetic models. Such structured in vivo, in vitro and in silico hazard data provide different lines of evidence which can be assembled, weighed and integrated using harmonised Weight of Evidence approaches to support the use of New Approach Methodologies (NAMs) in chemical risk assessment and the reduction of animal testing.

Do background levels of the pesticide pirimiphosmethyl in plant-based aquafeeds affect food safety of farmed Atlantic salmon?

The substitution of fish oil and fishmeal with plant-based ingredients in commercial aquafeeds for Atlantic salmon, may introduce novel contaminants that have not previously been associated with farmed fish. The organophosphate pesticide pirimiphos-methyl (PM) is one of the novel contaminants that is most prevalent in commercial salmon feed. In this study, the feed-to-fillet transfer of dietary PM and its main metabolites was investigated in Atlantic salmon fillet. Based on the experimental determined PM and metabolite uptake, metabolisation, and elimination kinetics, a physiologically based toxicokinetic (PBTK) compartmental model was developed. Fish fed PM had a relatively low (~4%) PM retention and two main metabolites (2-DAMP and Desethyl-PM) were identified in liver, muscle, kidney and bile. The absence of more metabolised forms of 2-DAMP and Desethyl-PM in Atlantic salmon indicates different metabolism in cold-water fish compared to previous studies on ruminants. The model was used to simulate the long term (>1.5 years) feed-to-fillet transfer of PM + metabolite in Atlantic salmon under realistic farming conditions including seasonal fluctuations in feed intake, growth, and fat deposition in muscle tissue. The model predictions show that with the constant presence of the highest observed PM concentration in commercial salmon feed, fillet PM+ metabolite levels were approximately 5 nmol kg^-1, with highest levels for the metabolite 2-DAMP. No EU maximum residue levels (MRL) for PM and its main metabolites exist in seafood to date, but the predicted levels were lower than the MRL for PM in swine of 32.7 nmol kg^-1.

Zebrafish as an Animal Model for Ocular Toxicity Testing: A Review of Ocular Anatomy and Functional Assays

Xenobiotics make their way into organisms from diverse sources including diet, medication, and pollution. Our understanding of ocular toxicities from xenobiotics in humans, livestock, and wildlife is growing thanks to laboratory animal models. Anatomy and physiology are conserved among vertebrate eyes, and studies with common mammalian preclinical species (rodent, dog) can predict human ocular toxicity. However, since the eye is susceptible to toxicities that may not involve a histological correlate, and these species rely heavily on smell and hearing to navigate their world, discovering visual deficits can be challenging with traditional animal models. Alternative models capable of identifying functional impacts on vision and requiring minimal amounts of chemical are valuable assets to toxicology. Human and zebrafish eyes are anatomically and functionally similar, and it has been reported that several common human ocular toxicants cause comparable toxicity in zebrafish. Vision develops rapidly in zebrafish; the tiny larvae rely on visual cues as early as 4 days, and behavioral responses to those cues can be monitored in high-throughput fashion. This article describes the comparative anatomy of the zebrafish eye, the notable differences from the mammalian eye, and presents practical applications of this underutilized model for assessment of ocular toxicity.

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.

An open source physiologically based kinetic model for the chicken (Gallus gallus domesticus): Calibration and validation for the prediction residues in tissues and eggs

Xenobiotics from anthropogenic and natural origin enter animal feed and human food as regulated compounds, environmental contaminants or as part of components of the diet. After dietary exposure, a chemical is absorbed and distributed systematically to a range of organs and tissues, metabolised, and excreted. Physiologically based kinetic (PBK) models have been developed to estimate internal concentrations from external doses. In this study, a generic multi-compartment PBK model was developed for chicken. The PBK model was implemented for seven compounds (with log K_ow range -1.37-6.2) to quantitatively link external dose and internal dose for risk assessment of chemicals. Global sensitivity analysis was performed for a hydrophilic and a lipophilic compound to identify the most sensitive parameters in the PBK model. Model predictions were compared to measured data according to dataset-specific exposure scenarios. Globally, 71% of the model predictions were within a 3-fold change of the measured data for chicken and only 7% of the PBK predictions were outside a 10-fold change. While most model input parameters still rely on in vivo experiments, in vitro data were also used as model input to predict internal concentration of the coccidiostat monensin. Future developments of generic PBK models in chicken and other species of relevance to animal health risk assessment are discussed.

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.

A physiologically based pharmacokinetic model of doxycycline for predicting tissue residues and withdrawal intervals in grass carp (Ctenopharyngodon idella)

The extensive use of doxycycline in aquaculture results in drug residue violations that negatively impact human food safety. This study aimed to develop a physiologically based pharmacokinetic (PBPK) model for doxycycline to predict drug residues and withdrawal times (WTs) in grass carp (Ctenopharyngodon idella) after daily oral administration for 3 days. Physiological parameters including cardiac output and organ weights were measured experimentally. Chemical-specific parameters were obtained from the literature or estimated by fitting to the observed data. The model properly captured the observed kinetic profiles of doxycycline in tissues (i.e., liver, kidney, muscle + skin and gill). The predicted WT in muscle + skin by Monte Carlo analysis based on sensitive parameters identified at 24 h after drug administration was 41 d, which was similar to 43 d calculated using the tolerance limit method. Sensitivity analysis identified two additional sensitive parameters at 6 weeks: intestinal transit rate constant and urinary elimination rate constant. The predicted WT in muscle + skin based on sensitive parameters identified at 6 weeks was 54 d. This model provides a useful tool to estimate tissue residues and withdrawal times for doxycycline in grass carp and also serves a foundation for extrapolation to other fish species and other tetracyclines.

Modelling BPA effects on three-spined stickleback population dynamics in mesocosms to improve the understanding of population effects

Bisphenol A (BPA), a well-known endocrine-disrupting chemical, is ubiquitously present in the aquatic environment. Its impacts at the population level on three-spined sticklebacks (Gasterosteus aculeatus) have been studied in artificial streams with low-dose BPA concentrations. The causes explaining the observed effects remained unclear. Here, we used an individual-based model coupled to a Dynamic Energy Budget model to (i) assess the potential of modelling to predict impacts at the population level using individual level laboratory ecotoxicological endpoints and (ii) provide insight on the mechanisms of BPA toxicity in these mesocosms. To do that, both direct and indirect effects of BPA on three-spined sticklebacks were incorporated in the model. Indeed, direct BPA effects on fish have been identified based on literature data whereas indirect effects on sticklebacks have been taken into account using sampling data of their prey from the exposed artificial streams. Results of the modelling showed that direct BPA effects on fish (impacts on gonad formation, growth, male reproductive behavior, eggs and larvae survival) mainly explained the three-spined stickleback population structure in the mesocosms, but indirect effects were not negligible. Hence, this study showed the potential of modelling in risk assessment to predict the impacts on fish population viability from behavioral and physiological effects measured on organisms.

Use of Zebrafish in Drug Discovery Toxicology

Unpredicted human safety events in clinical trials for new drugs are costly in terms of human health and money. The drug discovery industry attempts to minimize those events with diligent preclinical safety testing. Current standard practices are good at preventing toxic compounds from being tested in the clinic; however, false negative preclinical toxicity results are still a reality. Continual improvement must be pursued in the preclinical realm. Higher-quality therapies can be brought forward with more information about potential toxicities and associated mechanisms. The zebrafish model is a bridge between in vitro assays and mammalian in vivo studies. This model is powerful in its breadth of application and tractability for research. In the past two decades, our understanding of disease biology and drug toxicity has grown significantly owing to thousands of studies on this tiny vertebrate. This Review summarizes challenges and strengths of the model, discusses the 3Rs value that it can deliver, highlights translatable and untranslatable biology, and brings together reports from recent studies with zebrafish focusing on new drug discovery toxicology.