Exposure and hazard of bisphenol A, S and F: a multi-biomarker approach in three-spined stickleback.

Due to the estrogenic behavior of bisphenol (BP) A, industries have developed many substitutes, such as BPS and BPF. However, due to their structural similarities, adverse effects on reproduction are currently observed in various organisms, including fish. Even if new results have shown impacts of these bisphenols on many other physiological functions, their mode of action remains unclear. In this context, we proposed to better understand the impact of BPA, BPS, and BPF on immune responses (leucocyte sub-populations, cell death, respiratory burst, lysosomal presence, and phagocytic activity) and on biomarkers of metabolic detoxification (ethoxyresorufin-O-deethylase, EROD, and glutathione S-transferase, GST) and oxidative stress (glutathione peroxidase, GPx, and lipid peroxidation with thiobarbituric acid reactive substance method, TBARS) in an adult sentinel fish species, the three-spined stickleback. In order to enhance our understanding of how biomarkers change over time, it is essential to determine the internal concentration responsible for the observed responses. Therefore, it is necessary to explore the toxicokinetics of bisphenols. Thus, sticklebacks were exposed either to 100 µg/L of BPA, BPF or BPS for 21 days, or for seven days to 10 and 100 µg/L of BPA or BPS followed by seven days of depuration. Although BPS has very different TK, due to its lower bioaccumulation compared to BPA and BPF, BPS affect oxidative stress and phagocytic activity in the same way. For those reasons, the replacement of BPA by any substitute should be made carefully in terms of risk assessment on aquatic ecosystems.

PBTK-TD model of the phagocytosis activity in three-spined stickleback exposed to BPA.

Due to the high production volume and persistence in the environment of bisphenol A (BPA) and its substitutes, realistic exposure scenarii were proposed in some species to better understand the relationship between external and internal concentrations. For example, a recent PBTK model has been developed and adapted to BPA ADME (Absorption, Distribution, Metabolization, and Excretion) processes in three-spined stickleback. These substances have an impact on organism physiology including reproductive and immune functions. In this context, physiologically-based toxicokinetic models coupled with toxicodynamics (PBTK-TD) have proven to be valuable tools to fill the knowledge gap between external exposure and effect dynamics. The aim of the current work was to explain the impact of BPA on the immune response by determining its temporality. In addition, the relationship between BPA dose and these responses was investigated using a PBTK-TD model. Two experiments were performed on stickleback to characterize their biomarker responses, (i) a short exposure (14 days) at 0, 10 and 100 µg/L, including a depuration phase (7 days), and (ii) a long exposure (21 days) at 100 µg/L to measure the immunomarker dynamic over a long period. The fish spleens were sampled to analyze immune responses of stickleback at various times of exposure and depuration: leucocyte distribution, phagocytic capacity and efficiency, lysosomal presence and leucocyte respiratory burst index. At the same date, blood, muscle, and liver were sampled to quantify BPA and their metabolites (BPA monoglucuronide and BPA monosulfate). All these data enabled the development of the indirect pharmacodynamic models (PBTK-TD) by implementing the responses of biomarkers in the existing BPA PBTK of stickleback. The results shown a high induction of phagocytosis activity by BPA in the two exposure conditions. Furthermore, the immunomarkers exhibit very different temporal dynamics. This study demonstrates the need of a thorough characterization of biomarker response for a further use in Environmental Biomonitoring.

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.

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.