Invited Perspective: PFAS Bioconcentration and Biotransformation in Early Life Stage Zebrafish and Its Implications for Human Health Protection

Metabolomic changes following GenX and PFBS exposure in developing zebrafish

Short chain per- and polyfluoroalkyl substances (PFAS), including hexafluoropropylene oxide dimer acid (GenX) and perfluorobutane sulfonate (PFBS), are replacement chemicals for environmentally persistent, long-chain PFAS. Although GenX and PFBS have been detected in surface and ground water worldwide, few studies provide information on the metabolic alterations or risks associated with their exposures. In this study, larval zebrafish were used to investigate the toxicity of early-life exposure to GenX or PFBS. Zebrafish were chronically exposed from 4 h post-fertilization (hpf) to 6 days post-fertilization (dpf) to 150 µM GenX or 95.0 µM PFBS. Ultra-high-performance liquid chromatography paired with high-resolution mass spectrometry was used to quantify uptake of GenX and PFBS into zebrafish larvae and perform targeted and untargeted metabolomics. Our results indicate that PFBS was 20.4 % more readily absorbed into the zebrafish larvae compared to GenX. Additionally, PFBS exposure significantly altered 13 targeted metabolites and 21 metabolic pathways, while GenX exposure significantly altered 1 targeted metabolite and 17 metabolic pathways. Exposure to GenX, and to an even greater extent PFBS, resulted in a number of altered metabolic pathways in the amino acid metabolism, with other significant alterations in the carbohydrate, lipid, cofactors and vitamins, nucleotide, and xenobiotics metabolisms. Our results indicate that GenX and PFBS impact the zebrafish metabolome, with implications of global metabolic dysregulation, particularly in metabolic pathways relating to growth and development.

Survival, Growth, and Reproduction Responses in a Three-Generation Exposure of the Zebrafish (Danio rerio) to Perfluorooctane Sulfonate

A prior multigenerational perfluorooctane sulfonic acid (PFOS) exposure investigation in zebrafish reported adverse effects at 0.734 µg/L, among the lowest aquatic effect levels for PFOS reported to date. The present three-generation PFOS exposure quantified survival, growth, reproduction, and vitellogenin (VTG; egg yolk protein) responses in zebrafish, incorporating experimental design and procedural improvements relative to the earlier study. Exposures targeting 0.1, 0.6, 3.2, 20, and 100 µg/L in parental (P) and first filial (F1) generations lasted for 180 days post fertilization (dpf) and the second filial generation (F2) through 16 dpf. Survival decreased significantly in P and F2 generation exposures, but not in F1, at the highest PFOS treatment (100 µg/L nominal, 94-205 µg/L, measured). Significant adverse effects on body weight and length were infrequent, of low magnitude, and occurred predominantly at the highest exposure treatment. Finally, PFOS had no significant effects on P or F1 egg production and survival or whole-body VTG levels in P or F1 male fish. Overall, the predominance and magnitude of adverse PFOS effects at <1 µg/L reported in prior research were largely nonrepeatable in the present study. In contrast, the present study indicated a threshold for ecologically relevant adverse effects in zebrafish at 117 µg/L (SE 8 µg/L, n = 10) for survival and 47 µg/L (SE 11 µg/L, n = 19) for all statistically significant negative effects observed. Environ Toxicol Chem 2024;43:115-131. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Perfluorooctane Sulfonate (PFOS) Negatively Impacts Prey Capture Capabilities in Larval Zebrafish

Per- and polyfluoroalkyl substances (PFAS) are widely used in many industrial and domestic applications, which has resulted in unintentional human exposures and bioaccumulation in blood and other organs. Perfluorooctane sulfonate (PFOS) is among the most prevalent PFAS in the environment and has been postulated to affect brain functions in exposed organisms. However, the impacts of PFOS in early neural development have not been well described. We used zebrafish larvae to assess the effects of PFOS on two fundamental complex behaviors, prey capture and learning. Zebrafish exposed to PFOS concentrations ranging from 2 to 20 µM for differing 48-h periods were viable through early larval stages. In addition, PFOS uptake was unaffected by the presence of a chorion. We employed two different experimental paradigms; first we assessed the impacts of increasing organismal PFOS bioaccumulation on prey capture and learning, and second, we probed stage-specific sensitivity to PFOS by exposing zebrafish at different developmental stages (0-2 vs. 3-5 days post fertilization). Following both assays we measured the amount of PFOS present in each larva and found that PFOS levels varied in larvae from different groups within each experimental paradigm. Significant negative correlations were observed between larval PFOS accumulation and percentage of captured prey, whereas nonsignificant negative correlations were observed between PFOS accumulation and experienced-induced prey capture learning. These findings suggest that PFOS accumulation negatively affects larval zebrafish's ability to perform complicated multisensory behaviors and highlights the potential risks of PFOS exposure to animals in the wild, with implications for human health. Environ Toxicol Chem 2024;00:1-9. © 2023 SETAC.

Differential Developmental Neurotoxicity and Tissue Uptake of the Per- and Polyfluoroalkyl Substance Alternatives, GenX and PFBS

Per- and polyfluoroalkyl substances (PFAS) are a class of synthetic chemicals with several applications. Multiple adverse health effects are reported for longer carbon chain (≤C8) PFAS. Shorter carbon chain PFAS, [e.g., hexafluoropropylene oxide dimer acid (HFPO-DA; GenX) and perfluorobutanesulfonic acid (PFBS)] were introduced as alternatives. Past studies indicate that longer-chain PFAS are neurotoxic targeting the dopamine pathway, but it is not known if shorter-chain PFAS act similarly. This study aimed to evaluate developmental neurotoxicity and tissue uptake of GenX and PFBS using the zebrafish (Danio rerio). First, acute toxicity was assessed by measuring LC50 at 120 h postfertilization (hpf). Body burden was determined after embryonic exposure (1-72 hpf) to sublethal concentrations of GenX or PFBS by LC-ESI-MS/MS. Locomotor activity using a visual motor response assay at 120 hpf and dopamine levels at 72 hpf was assessed after embryonic exposure. PFBS was more acutely toxic and bioaccumulative than GenX. GenX and PFBS caused hyperactivity at 120 hpf, but stronger behavioral alterations were observed for PFBS. An increase in whole organism dopamine occurred at 40 ppb of GenX, while a decrease was observed at 400 ppb of PFBS. Differences detected in dopamine for these two PFAS indicate differential mechanisms of developmental neurotoxicity.

PBTK modeled perfluoroalkyl acid kinetics in zebrafish eleutheroembryos suggests impacts on bioconcentrations by chorion porosity dynamics

The zebrafish eleutheroembryo (zfe) is widely used as a model to characterize the toxicity of chemicals. However, analytical methods are still missing to measure organ concentrations. Therefore, physiologically-based toxicokinetic (PBTK) modeling may overcome current limitations to help understand the relationship between toxic effects and internal exposure in various organs. A previous PBTK model has been updated to include the chorionic transport barrier and its permeabilization, hatching dynamics within a zfe population over development, and active mediated transport mechanisms. The zfe PBTK model has been calibrated using measured time-dependent internal concentrations of PFBA, PFHxS, PFOA, and PFOS in a zfe population and evaluated using external datasets from the literature. Calibration was successful with 96% of the predictions falling within a 2-fold range of the observed concentrations. The external dataset was correctly estimated with about 50% of the predictions falling within a factor of 3 of the observed data and 10% of the predictions are out of the 10-fold error. The calibrated model suggested that active mediated transport differs between PFAS with a sulfonic and carboxylic acid functional end groups. This PBTK model predicts well the fate of PFAS with various physicochemical properties in zfe. Therefore, this model may improve the use of zfe as an alternative model in toxicokinetic-toxicodynamic studies and help to refine and reduce zfe-based experiments, while giving insights into the internal kinetics of chemicals.

A Scoping Assessment of Implemented Toxicokinetic Models of Per- and Polyfluoro-Alkyl Substances, with a Focus on One-Compartment Models

Toxicokinetic (TK) models have been used for decades to estimate concentrations of per-and polyfluoroalkyl substances (PFAS) in serum. However, model complexity has varied across studies depending on the application and the state of the science. This scoping effort seeks to systematically map the current landscape of PFAS TK models by categorizing different trends and similarities across model type, PFAS, and use scenario. A literature review using Web of Science and SWIFT-Review was used to identify TK models used for PFAS. The assessment covered publications from 2005-2020. PFOA, the PFAS for which most models were designed, was included in 69 of the 92 papers, followed by PFOS with 60, PFHxS with 22, and PFNA with 15. Only 4 of the 92 papers did not include analysis of PFOA, PFOS, PFNA, or PFHxS. Within the corpus, 50 papers contained a one-compartment model, 17 two-compartment models were found, and 33 used physiologically based pharmacokinetic (PBTK) models. The scoping assessment suggests that scientific interest has centered around two chemicals-PFOA and PFOS-and most analyses use one-compartment models in human exposure scenarios.

A Machine Learning Model to Estimate Toxicokinetic Half-Lives of Per- and Polyfluoro-Alkyl Substances (PFAS) in Multiple Species

Per- and polyfluoroalkyl substances (PFAS) are a diverse group of man-made chemicals that are commonly found in body tissues. The toxicokinetics of most PFAS are currently uncharacterized, but long half-lives (t_½) have been observed in some cases. Knowledge of chemical-specific t_½ is necessary for exposure reconstruction and extrapolation from toxicological studies. We used an ensemble machine learning method, random forest, to model the existing in vivo measured t_½ across four species (human, monkey, rat, mouse) and eleven PFAS. Mechanistically motivated descriptors were examined, including two types of surrogates for renal transporters: (1) physiological descriptors, including kidney geometry, for renal transporter expression and (2) structural similarity of defluorinated PFAS to endogenous chemicals for transporter affinity. We developed a classification model for t_½ (Bin 1: <12 h; Bin 2: <1 week; Bin 3: <2 months; Bin 4: >2 months). The model had an accuracy of 86.1% in contrast to 32.2% for a y-randomized null model. A total of 3890 compounds were within domain of the model, and t_½ was predicted using the bin medians: 4.9 h, 2.2 days, 33 days, and 3.3 years. For human t_½, 56% of PFAS were classified in Bin 4, 7% were classified in Bin 3, and 37% were classified in Bin 2. This model synthesizes the limited available data to allow tentative extrapolation and prioritization.

Sublethal Exposure of Per- and Polyfluoroalkyl Substances of Varying Chain Length and Polar Functionality Results in Distinct Metabolic Responses in Daphnia magna

Per- and polyfluoroalkyl substances (PFAS) are a class of persistent organic pollutants used in industrial applications because of their physicochemical properties, which results in their ubiquitous presence across environmental matrices. To date, legacy PFAS have been well studied; however, the concentration of alternative PFAS may exceed the concentration of legacy pollutants, and more information is needed regarding the sublethal toxicity at the molecular level of aquatic model organisms, such as Daphnia magna. Perfluorobutanoic acid (PFBA), perfluorohexanoic acid (PFHxA), perfluorohexanesulfonic acid (PFHxS), and perfluorononanoic acid (PFNA) are four widely detected PFAS alternatives of varying chain length and polar functionality that are quantified in aquatic environments. The present study examines the metabolic perturbations of PFAS with varying chemistries to D. magna using targeted mass spectrometry-based metabolomics. Daphnia were acutely exposed to sublethal concentrations of PFBA, PFHxA, PFHxS, and PFNA before the polar metabolite profile was extracted from single organisms. Multivariate analysis demonstrated significant separation between the sublethal concentrations of PFHxA, PFHxS, and PFNA relative to the controls; in sum, longer chain lengths demonstrated greater overall perturbations to the extracted metabolic profiles. Univariate statistics revealed significant perturbations in the concentrations of several amino acids, nucleotides/nucleosides, and neurotransmitters with exposure to PFAS. These metabolic perturbations are consistent with disruptions in energy metabolism (pantothenate and coenzyme A metabolism, histidine metabolism) and protein synthesis (aminoacyl-transfer RNA biosynthesis and amino acid metabolism), which were identified through biochemical pathway analysis. These results provide evidence that although PFAS chemistry (chain length and polar functional group) invokes unique metabolic responses, there is also an underlying toxic mode of action that is common with select PFAS exposure. Overall, the present study highlights the capabilities of environmental metabolomics to elucidate the molecular-level perturbations of pollutants within the same chemical class to model aquatic organisms, which can be used to prioritize risk assessment of substituted PFAS alternatives. Environ Toxicol Chem 2023;42:242-256. © 2022 SETAC.

Comparative toxicokinetics and toxicity of PFOA and its replacement GenX in the early stages of zebrafish

PER: and poly-fluoroalkyl substances (PFAS) are receiving attention due to their persistence, and potential adverse effects on environmental and human health. Efforts to reduce long-chained PFAS (≥C8) compounds were implemented in 2006 as a part of "PFOA Stewardship Program Initiative" (PFOA-perfluorooctanoic acid). Short-chained PFAS (<C8) were introduced as replacements, which were believed to have lower potential for environmental persistence and bioaccumulation. Little is known about the uptake and elimination, and potential toxic effects of these replacement compounds. Hence, it is important to compare toxicokinetics and toxicity of long-chain PFAS to their replacement compounds. To this end, zebrafish (ZF), Danio rerio, embryos were exposed to PFOA and its short-chain replacement perfluoro (2-methyl-3-oxahexanoic) acid (GenX) with the aim to assess uptake and elimination kinetics, hatching success, morphology, startle response, and survival. At 24 hpf, LC50 was 82 μM for PFOA and 170 μM for GenX. At 54 hpf, GenX but not PFOA showed an increase in hatching success. At 120 hpf, no statistically significant differences were seen in white light startle response below the LC50. PFAS internal concentrations were measured at 72 and 120 hpf during exposure phase, and at 168 hpf during depuration phase. GenX and PFOA internal concentrations in 120 hpf larvae exposed to highest concentration (20 μM) were 35.02 and 44.51 μM, respectively. Concentrations were eliminated almost completely at 168 hpf for GenX up to 95%, while for PFOA up to 50%. As steady-state was not reached, we estimated kinetic bioconcentration factors (BCFkin). BCFkin for GenX was lower than PFOA at equimolar concentrations. However, bioconcentration factors were higher at the lower exposure concentrations for both chemicals, suggesting a concentration-dependent uptake of PFASs. The predicted internal effect concentrations, accounting for the differences in bioconcentration factors, were 229 μM for GenX and 226 μM for PFOA, suggesting similar toxic potencies.

Kinetics and toxicity of an environmentally relevant mixture of halogenated organic compounds in zebrafish embryo

Persistent and semi-persistent halogenated compounds cause health problems for the animals occupying the upper level of the food web in the Baltic Sea. Atlantic salmon (Salmo salar), being a top piscivore in the Baltic Sea, has been observed to carry a large body burden of halogenated toxins. Here, a mixture of nine halogenated compounds belonging to different groups was created, based on the observed composition of halogenated toxins in salmon serum. The toxicokinetic properties of the compounds were studied in zebrafish (Danio rerio) embryos to achieve the same proportions between the internal doses of the compounds in the zebrafish as in the salmon. Toxicity was evaluated for the compounds dosed individually as well as in a mixture. Perfluorooctanesulfonic acid (PFOS) was the dominant compound in the salmon and was observed to be the driving force for effects on swimbladder inflation caused by the mixture with a 50% effect concentration of 4.8 µM nominal dose, or 1300 µMD based on the area under the internal concentration-time curve (AUC). The driving compound for other severe effects caused by the mixture, including lethality, spinal deformity, and edemas, was the hydroxylated polybrominated diphenyl ether 6-OH-BDE47, which was observed to have a 50% lethality concentration of 93 nM, corresponding to 94 µMD based on internal dose (AUC). The individual compounds were observed to act additively on most of the documented outcomes when dosed as a mixture.

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

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

Comparison of sub-lethal metabolic perturbations of select legacy and novel perfluorinated alkyl substances (PFAS) in Daphnia magna

Per- and polyfluoroalkyl substances (PFAS) are a class of pollutants of concern due to their ubiquitous presence, persistence, and toxicity in aquatic environments. Legacy PFAS pollutants such as perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) have been more widely studied in aquatic environments. However, replacement PFAS, such as ammonium perfluoro (2-methyl-3-oxahexanoate; GenX) are increasingly being detected with little known information surrounding their toxicity. Here, Daphnia magna, a model organism for freshwater ecotoxicology was used to compare the acute sub-lethal toxicity of PFOS, PFOA, GenX, and PFAS mixtures. Using liquid chromatography with tandem mass spectrometry (LC-MS/MS), the targeted polar metabolic profile extracted from single Daphnia was quantified to investigate perturbations in the exposure groups versus the unexposed organisms. Multivariate statistical analyses demonstrated significant non-monotonic separation in PFOA, GenX, and PFAS mixture exposures. Sub-lethal exposure to concentrations of PFOS did not lead to significant separation in multivariate analyses. Univariate statistics and pathway analyses were used to elucidate the mode of action of PFAS exposure. Exposure to all individual PFAS led to significant perturbations in many amino acids including cysteine, histidine, tryptophan, glycine, and serine. These perturbations are consistent with biochemical pathway disruptions in the pantothenate and Coenzyme A (CoA) biosynthesis, thiamine metabolism, histidine metabolism, and aminoacyl-tRNA biosynthesis pathways. Overall, the collected metabolomic data is consistent with disruptions in energy metabolism and protein synthesis as the primary mode of action of sub-lethal PFAS exposure. Secondary modes of action among individual pollutant exposures demonstrated that the structural properties (carboxylic acid vs. sulfonic acid group) may play a role in the metabolic perturbations observed. Sub-lethal exposure to PFAS mixtures highlighted a mixed response when compared to the individual pollutants (PFOS, PFOA, and GenX). Overall, this study emphasizes the niche capability of environmental metabolomics to differentiate secondary modes of action from metabolic perturbations in both single pollutant and pollutant mixtures within the same chemical class.

Global Exposure to Per- and Polyfluoroalkyl Substances and Associated Burden of Low Birthweight

Low birthweight (LBW) is a worldwide public health concern, while the global burden of LBW attributable to endocrine-disrupting chemicals, such as per- and polyfluoroalkyl substances (PFAS), has not yet been evaluated. Here, we established a large dataset for the biomonitoring of seven representative congeners of PFAS by examining data from 2325 publications. Global exposure to perfluorooctanesulfonic acid (PFOS) was the highest, followed by perfluorohexanesulfonic acid (PFHxS) and perfluorooctanoic acid (PFOA). Spatiotemporal exposure to PFAS varied considerably, with daily intake estimated in the range of 0.01-1.7 ng/kg/day. Moreover, decreasing trends in PFOS, PFHxS, and PFOA exposure were noted in most regions of the world over the past two decades, but such trends were not observed for other PFAS with long carbon chains, especially in East Asia. Furthermore, we estimated that human exposure to PFOA contributed to approximately 461,635 (95% confidence interval: 57,418 to 854,645) cases per year of LBW during the past two decades, predominantly from Asian regions. Although our estimation may be constrained by uncertainties from the dose-response curve and data availability, this study has unveiled that PFAS might be a contributor to global LBW prevalence during 2000-2019, supporting continuous actions to mitigate PFAS contamination.