Formation of PFAAs in fish through biotransformation: A PBPK approach

Role of niche and micro-habitat preferences in per- and polyfluoroalkyl substances occurrence in the gills of tropical lake fish species

The present study has investigated per- and poly-fluoroalkyl substances (PFAS) in the gill tissues of various fish species inhabiting different trophic levels within Eleyele Lake, a tropical freshwater lake in Nigeria. The mean concentrations of PFAS congeners were determined, and their trends and patterns were analyzed across different trophic species. The results revealed variations in congener abundance and species-specific patterns that was influenced by habitat and niche preferences. Multivariate associations using canonical-correlation analysis (CCA) revealed distinct trends in the relationships between gill concentrations of specific PFAS congeners and different trophic groups. The strongest congener relationships were observed in the pelagic omnivore (Oreochromic niloticus: ON) with positive associations for 4:2 FTS, 9CL-PF3ONS, PFTDA, MeFOSA and PFHxS. The differences in congener profiles for the two herbivorous fish (Sarotherodon melanotheron (SM) and Coptodon galilaeus (CG)) reflect possible divergence in microhabitat and niche preferences. Furthermore, the congener overlaps between the herbivore (CG), and benthic omnivore (Clarias gariepinus: ClG) indicate a possible niche and microhabitat overlap. Our study provides valuable insights into the congener dynamics of PFAS at Eleyele Lake. However, the dissimilarity and overlapping PFAS congener profile in fish gills reflects the interplay of species niche preference and microhabitat associations. The present study highlights the need for further research to assess ecological risks and develop effective PFAS management strategies.

Towards process-based modelling and parameterisation of bioaccumulation in humans across PFAS congeners

Per- and polyfluoroalkyl substances (PFAS) are a large class of stable toxic chemicals which have ended up in the environment and in organisms in significant concentrations. Toxicokinetic models are needed to facilitate extrapolation of bioaccumulation data across PFAS congeners and species. For the present study, we carried out an inventory of accumulation processes specific for PFAS, deviating from traditional Persistent Organic Pollutants (POPs). In addition, we reviewed toxicokinetic models on PFAS reported in literature, classifying them according to the number of compartments distinguished as a one-compartment model (1-CM), two-compartment model (2- CM) or a multi-compartment model, (multi-CM) as well as the accumulation processes included and the parameters used. As the inventory showed that simple 1-CMs were lacking, we developed a generic 1-CM of ourselves to include PFAS specific processes and validated the model for legacy perfluoroalkyl acids. Predicted summed elimination constants were accurate for long carbon chains (>C6), indicating that the model properly represented toxicokinetic processes for most congeners. Results for urinary elimination rate constants were mixed, which might be caused by the exclusion of reabsorption processes (renal reabsorption, enterohepatic circulation). The 1-CM needs to be improved further in order to better predict individual elimination pathways. Besides that, more data on PFAS-transporter specific processes are needed to extrapolate across PFAS congeners and species.

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.

A State-of-the-Science Review of Interactions of Per- and Polyfluoroalkyl Substances (PFAS) with Renal Transporters in Health and Disease: Implications for Population Variability in PFAS Toxicokinetics

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) are ubiquitous in the environment and have been shown to cause various adverse health impacts. In animals, sex- and species-specific differences in PFAS elimination half-lives have been linked to the activity of kidney transporters. However, PFAS molecular interactions with kidney transporters are still not fully understood. Moreover, the impact of kidney disease on PFAS elimination remains unclear.

OBJECTIVES

This state-of-the-science review integrated current knowledge to assess how changes in kidney function and transporter expression from health to disease could affect PFAS toxicokinetics and identified priority research gaps that should be addressed to advance knowledge.

METHODS

We searched for studies that measured PFAS uptake by kidney transporters, quantified transporter-level changes associated with kidney disease status, and developed PFAS pharmacokinetic models. We then used two databases to identify untested kidney transporters that have the potential for PFAS transport based on their endogenous substrates. Finally, we used an existing pharmacokinetic model for perfluorooctanoic acid (PFOA) in male rats to explore the influence of transporter expression levels, glomerular filtration rate (GFR), and serum albumin on serum half-lives.

RESULTS

The literature search identified nine human and eight rat kidney transporters that were previously investigated for their ability to transport PFAS, as well as seven human and three rat transporters that were confirmed to transport specific PFAS. We proposed a candidate list of seven untested kidney transporters with the potential for PFAS transport. Model results indicated PFOA toxicokinetics were more influenced by changes in GFR than in transporter expression.

DISCUSSION

Studies on additional transporters, particularly efflux transporters, and on more PFAS, especially current-use PFAS, are needed to better cover the role of transporters across the PFAS class. Remaining research gaps in transporter expression changes in specific kidney disease states could limit the effectiveness of risk assessment and prevent identification of vulnerable populations. https://doi.org/10.1289/EHP11885.

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.

Field-Based Distribution and Bioaccumulation Factors for Cyclic and Aliphatic Per- and Polyfluoroalkyl Substances (PFASs) in an Urban Sedentary Waterbird Population

The field-based distribution and bioaccumulation factor (BAF) for per- and polyfluoroalkyl substances (PFASs) were determined in residential Black Swans (Cygnus atratus) from an urban lake (Melbourne, Australia). The concentrations of 46 aliphatic and cyclic PFASs were determined by HPLC-MS/MS in serum and excrement from swans, and water, sediment, aquatic macrophytes, soil, and grass samples in and around the lake. Elevated concentrations of ∑₄₆PFASs were detected in serum (120 ng mL^-1) and excrement (110 ng g^-1 dw) were strongly related indicating a potential noninvasive sampling methodology. Environmental concentrations of PFASs were consistent with a highly impacted ecosystem and notably high concentrations of perfluoro-4-ethylcyclohexanesulfonate (PFECHS, 67584-42-3; C₈HF₁₅SO₃) were detected in water (27 ng L^-1) and swan serum (16 ng mL^-1). In the absence of credible putative alternative sources of PFECHS input to the lake, we propose that the use of high-performance motorsport vehicles is a likely source of contamination to this ecosystem. The BAF of perfluorocarboxylic acids increased with each additional CF₂ moiety from PFOA (15.7 L kg^-1 ww) to PFDoDA (3615 L kg^-1 ww). The BAF of PFECHS was estimated as 593 L kg^-1 ww, which is lower compared with that of PFOS (1097 L kg^-1 ww).

The impact of precursors on aquatic exposure assessment for PFAS: Insights from bioaccumulation modeling

Risk assessment for per- and polyfluoroalkyl substances (PFAS) is complicated by the fact that PFAS include several thousand compounds. Although new analytical methods have increased the number that can be identified in environmental samples, a significant fraction of them remain uncharacterized. Perfluorooctane sulfonate (PFOS) is the PFAS compound of primary interest when evaluating risks to humans and wildlife owing to the consumption of aquatic organisms. The exposure assessment for PFOS is complicated by the presence of PFOS precursors and their transformation, which can occur both in the environment and within organisms. Thus, the PFOS to which wildlife or people are exposed may consist of PFOS that was discharged directly into the environment and/or other PFOS precursors that were transformed into PFOS. This means that exposure assessment and the development of remedial strategies may depend on the relative concentrations and properties not only of PFOS but also of other PFAS that are transformed into PFOS. A bioaccumulation model was developed to explore these issues. The model embeds toxicokinetic and bioenergetic components within a larger food web calculation that accounts for uptake from both food and water, as well as predator-prey interactions. Multiple chemicals are modeled, including parent-daughter reactions. A series of illustrative simulations explores how chemical properties can influence exposure assessment and remedial decision making. Integr Environ Assess Manag 2021;17:705-715. © 2021 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Bioaccumulation, Biodistribution, Toxicology and Biomonitoring of Organofluorine Compounds in Aquatic Organisms

This review is a survey of recent advances in studies concerning the impact of poly- and perfluorinated organic compounds in aquatic organisms. After a brief introduction on poly- and perfluorinated compounds (PFCs) features, an overview of recent monitoring studies is reported illustrating ranges of recorded concentrations in water, sediments, and species. Besides presenting general concepts defining bioaccumulative potential and its indicators, the biodistribution of PFCs is described taking in consideration different tissues/organs of the investigated species as well as differences between studies in the wild or under controlled laboratory conditions. The potential use of species as bioindicators for biomonitoring studies are discussed and data are summarized in a table reporting the number of monitored PFCs and their total concentration as a function of investigated species. Moreover, biomolecular effects on taxonomically different species are illustrated. In the final paragraph, main findings have been summarized and possible solutions to environmental threats posed by PFCs in the aquatic environment are discussed.

Evaluation of Published Bioconcentration Factor (BCF) and Bioaccumulation Factor (BAF) Data for Per- and Polyfluoroalkyl Substances Across Aquatic Species

Per- and polyfluoroalkyl substances (PFAS) are a group of man-made chemicals of concern across the globe, and some of the PFAS chemicals are known to be bioaccumulative in aquatic species. A literature search for bioconcentration factors (BCFs) and bioaccumulation factors (BAFs) for PFAS has been done, and data for 22 taxonomic classes were assembled. The assembled data were evaluated for quality, and for gaps and limitations in bioaccumulation information for the PFAS universe of chemicals. In general, carbonyl and sulfonyl PFAS classes are relatively data rich, whereas phosphate, fluorotelomer, and ether PFAS classes are data limited for fish and nonexistent for most other taxonomic classes. Taxonomic classes with the most measurements were, in descending order, Teleostei (fish), Bivalvia, and Malacostraca. For fish, median whole-body log BAFs (L/kg wet wt) for perfluorooctane sulfonic acid and perfluorooctanoic acid were 3.55 (standard deviation [SD] = 0.83, n = 84) and 2.16 (SD = 0.85, n = 48) using all measurements, respectively. In comparison with freshwater species, data are limited for marine species, and further research is needed to determine whether the BAFs for freshwater and marine species should be the same or different. The BAFs for some PFAS appear to be consistent with the BCFs developed with laboratory experiments, in which values decline with increasing concentrations in water. Environ Toxicol Chem 2021;40:1530-1543. Published 2021. This article is a U.S. Government work and is in the public domain in the USA.

Elevated bioaccumulation of PFAAs in Oryzias melastigma following the increase of salinity is associated with the up-regulated expression of PFAA-binding proteins

Perfluoroalkyl acids (PFAAs) are widely detected in the environment, especially in estuarine and coastal areas where fluctuation of salinity occurs. Salinity alteration affected the distribution of PFAAs and even the bioaccumulation in organisms. However, the inner mechanism is still unclear. In this study, the marine medaka (Oryzias melastigma), a euryhaline fish model, was exposed to four PFAAs congeners under three different salinities (0, 15 and 35 psu). Results showed that the bioaccumulation of PFAAs increased in fish as the water salinity increased. PFAAs with longer lengths of carbon‑fluorine bond showed higher bioaccumulation in the fish. Salinity did not alter the levels of PFAAs in water media, however, the uptake rate of PFAAs from gills did increase with the salinity. Further analysis of the mechanism showed that PFAA bound to branchial proteins as confirmed by fluorescence spectroscopy. Higher expressions of proteins binding to PFAAs including organic anion transporter 1 (OAT1) and fatty acid-binding protein (FABP) facilitated the uptake of PFAAs through gills in fish culturing under higher salinity. In all, our study showed that elevation of salinity can induce the expression of proteins binding to PFAAs in gills, thus facilitate the uptake of water PFAAs. Salinity fluctuation should be taken into consideration when assessing the chemical risk in the estuarine and coastal areas.

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.

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.