Physiologically based pharmacokinetic modeling of human exposure to perfluorooctanoic acid suggests historical non drinking-water exposures are important for predicting current serum concentrations

Chemical risk assessment in food animals via physiologically based pharmacokinetic modeling - Part II: Environmental pollutants on animal and human health assessments

Human activities generate a large amount of environmental pollutants, including drugs and agricultural and industrial chemicals that are released into the air, water, and soil. Environmental pollutants can enter food animals through contaminated feed and water, posing risks to human health via the food chain. Physiologically based pharmacokinetic (PBPK) modeling is used to predict the target organ dosimetry informing human health risk assessment. However, there is a lack of critical reviews concerning PBPK models for environmental pollutants in food animals in the last several years (2020-2024). This review is part of a series of reviews focusing on applications of PBPK models for drugs and environmental chemicals in food animals to inform human health and food safety assessments. Part I is focused on veterinary drugs. The present article is Part II and focuses on environmental chemicals, including pesticides, polychlorinated biphenyls (PCBs), bisphenols, and per- and polyfluoroalkyl substances (PFAS). This article discusses the existing challenges in developing PBPK models for environmental pollutants and shares our perspectives on future directions, including the combinations of in vitro to in vivo extrapolation (IVIVE), machine learning and artificial intelligence, read-across approaches, and quantitative pharmacodynamic modeling to enhance the potential applications of PBPK models in assessing human health and food safety.

Advancing understanding of human variability through toxicokinetic modeling, in vitro-in vivo extrapolation, and new approach methodologies

The merging of physiology and toxicokinetics, or pharmacokinetics, with computational modeling to characterize dosimetry has led to major advances for both the chemical and pharmaceutical research arenas. Driven by the mutual need to estimate internal exposures where in vivo data generation was simply not possible, the application of toxicokinetic modeling has grown exponentially in the past 30 years. In toxicology the need has been the derivation of quantitative estimates of toxicokinetic and toxicodynamic variability to evaluate the suitability of the tenfold uncertainty factor employed in risk assessment decision-making. Consideration of a host of physiologic, ontogenetic, genetic, and exposure factors are all required for comprehensive characterization. Fortunately, the underlying framework of physiologically based toxicokinetic models can accommodate these inputs, in addition to being amenable to capturing time-varying dynamics. Meanwhile, international interest in advancing new approach methodologies has fueled the generation of in vitro toxicity and toxicokinetic data that can be applied in in vitro-in vivo extrapolation approaches to provide human-specific risk-based information for historically data-poor chemicals. This review will provide a brief introduction to the structure and evolution of toxicokinetic and physiologically based toxicokinetic models as they advanced to incorporate variability and a wide range of complex exposure scenarios. This will be followed by a state of the science update describing current and emerging experimental and modeling strategies for population and life-stage variability, including the increasing application of in vitro-in vivo extrapolation with physiologically based toxicokinetic models in pharmaceutical and chemical safety research. The review will conclude with case study examples demonstrating novel applications of physiologically based toxicokinetic modeling and an update on its applications for regulatory decision-making. Physiologically based toxicokinetic modeling provides a sound framework for variability evaluation in chemical risk assessment.

Dissect the association between per- and polyfluoroalkyl substances (PFAS) and kidney function from the perspective of lipid molecules

Per- and polyfluoroalkyl substances (PFAS) have been linked to kidney function. Studies have shown that PFAS can cause changes in lipid metabolism and that lipids play an important role in regulating kidney function. However, few studies have explored the overall impact of PFAS mixture on kidney function. Moreover, the mechanisms by which PFAS influences kidney function remain unclear. This study was performed to investigate the overall impact of PFAS mixture on kidney function indexes, dissect the mechanism by which PFAS affect kidney function by analyzing lipid molecule profiles, and analyze the associations between different subclasses of lipids and kidney function indexes. We measured blood PFAS levels and kidney function indexes in a community population containing 278 males. Metabolomic analysis detected 332 lipid molecules. A quantile-based g-computation model was applied to assess the overall effect of PFAS mixture on kidney function index, and revealed that PFAS mixture were associated with a higher level of uric acid (UA). Linear regression analysis demonstrated a positive association between PFOA and UA, and logistic regression analysis indicated a positive association between PFOA and hyperuricemia odds. Notably, none of the PFAS were associated with the estimated glomerular filtration rate, indicating that PFAS didn't have an obvious effect on glomerular filtration. Further analysis identified 20 lipid molecules associated with both PFOA and UA. High-dimensional mediation effect analysis showed that seven lipid molecules (one glycerophospholipid, three fatty acyls, and three prenol lipids) mediated the association between PFOA and UA. Additionally, quantile-based g-computation analysis revealed positive associations between specific lipid subclasses-mainly fatty acid esters, fatty acids and conjugates, and sesquiterpenoids-and kidney function indexes. Our findings provide insights into the renal toxicity of PFAS and may also lead to more in-depth investigations using animal models and other population studies.

Development of a Physiologically Based Pharmacokinetic (PBPK) Model for F-53B in Pregnant Mice and Its Extrapolation to Humans

Chlorinated polyfluorinated ether sulfonic acid (F-53B), a commonly utilized alternative for perfluorooctane sulfonate, was detected in pregnant women and cord blood recently. However, the lack of detailed toxicokinetic information poses a significant challenge in assessing the human risk assessment for F-53B exposure. Our study aimed to develop a physiologically based pharmacokinetic (PBPK) model for pregnant mice, based on toxicokinetic experiments, and extrapolating it to humans. Pregnant mice were administered 80 μg/kg F-53B orally and intravenously on gestational day 13. F-53B concentrations in biological samples were analyzed via ultraperformance liquid chromatography-mass spectrometry. Results showed the highest F-53B accumulation in the brain, followed by the placenta, amniotic fluid, and liver in fetal mice. These toxicokinetic data were applied to F-53B PBPK model development and evaluation, and Monte Carlo simulations were used to characterize the variability and uncertainty in the human population. Most of the predictive values were within a 2-fold range of experimental data (>72%) and had a coefficient of determination (R2) greater than 0.68. The developed mouse model was then extrapolated to the human and evaluated with human biomonitoring data. Our study provides an important step toward improving the understanding of toxicokinetics of F-53B and enhancing the quantitative risk assessments in sensitive populations, particularly in pregnant women and fetuses.

Contribution of Continued Dermal Exposure of PFAS-Containing Sunscreens to Internal Exposure: Extrapolation from In Vitro and In Vivo Tests to Physiologically Based Toxicokinetic Models

Per- and polyfluoroalkyl substances (PFASs) are widely present in sunscreen products as either active ingredients or impurities. They may penetrate the human skin barrier and then pose potential health risks. Herein, we aimed to develop a physiologically based toxicokinetic (PBTK) model capable of predicting the body loading of PFASs after repeated, long-term dermal application of commercial sunscreens. Ten laboratory-prepared sunscreens, generally falling into two categories of water-in-oil (W/O) and oil-in-water (O/W) sunscreens, were subject to in vitro percutaneous penetration test to assess the impacts of four sunscreen ingredients on PFAS penetration. According to the results, two sunscreen formulas representing W/O and O/W types that mostly enhanced PFAS dermal absorption were then selected for a subsequent 30 day in vivo exposure experiment in mice. PBTK models were successfully established based on the time-dependent PFAS concentrations in mouse tissues (R2 = 0.885-0.947) and validated through another 30 day repeated exposure experiment in mice using two commercially available sunscreens containing PFASs (R2 = 0.809-0.835). The PBTK model results suggest that applying sunscreen of the same amount on a larger skin area is more conducive to PFAS permeation, thus enhancing the exposure risk. This emphasizes the need for caution in practical sunscreen application scenarios, particularly during the summer months.

PI3K/Akt/FoxO Pathway Mediates Antagonistic Toxicity in HepG2 Cells Coexposed to Deoxynivalenol and Enniatins

The emerging mycotoxins enniatins (ENNs) and the traditional mycotoxin deoxynivalenol (DON) often co-contaminate various grain raw materials and foods. While the liver is their common target organ, the mechanism of their combined effect remains unclear. In this study, the combined cytotoxic effects of four ENNs (ENA, ENA1, ENB, and ENB1) with DON and their mechanisms were investigated using the HepG2 cell line. Additionally, a population exposure risk assessment of these mycotoxins was performed by using in vitro experiments and computer simulations. The results showed that only ENA at 1/4 IC50 and ENB1 at 1/8 IC50 coexposed with DON showed an additive effect, while ENB showed the strongest antagonism at IC50 (CI = 3.890). Co-incubation of ENNs regulated the signaling molecule levels which were disrupted by DON. Transcriptome analysis showed that ENB (IC50) up-regulated the PI3K/Akt/FoxO signaling pathway and inhibited the expression of apoptotic genes (Bax, P53, Caspase 3, etc.) via phosphorylation of FoxO, thereby reducing the cytotoxic effects caused by DON. Both types of mycotoxins posed serious health risks, and the cumulative risk of coexposure was particularly important for emerging mycotoxins.

Estimating the dynamic early life exposure to PFOA and PFOS of the HELIX children: Emerging profiles via prenatal exposure, breastfeeding, and diet

In utero and children's exposure to per- and polyfluoroalkyl substances (PFAS) is a major concern in health risk assessment as early life exposures are suspected to induce adverse health effects. Our work aims to estimate children's exposure (from birth to 12 years old) to PFOA and PFOS, using a Physiologically-Based Pharmacokinetic (PBPK) modelling approach. A model for PFAS was updated to simulate the internal PFAS exposures during the in utero life and childhood, and including individual characteristics and exposure scenarios (e.g., duration of breastfeeding, weight at birth, etc.). Our approach was applied to the HELIX cohort, involving 1,239 mother-child pairs with measured PFOA and PFOS plasma concentrations at two sampling times: maternal and child plasma concentrations (6 to 12 y.o). Our model predicted an increase in plasma concentrations during fetal development and childhood until 2 y.o when the maximum concentrations were reached. Higher plasma concentrations of PFOA than PFOS were predicted until 2 y.o, and then PFOS concentrations gradually became higher than PFOA concentrations. From 2 to 8 y.o, mean concentrations decreased from 3.1 to 1.88 µg/L or ng/mL (PFOA) and from 4.77 to 3.56 µg/L (PFOS). The concentration-time profiles vary with the age and were mostly influenced by in utero exposure (on the first 4 months after birth), breastfeeding (from 5 months to 2 (PFOA) or 5 (PFOS) y.o of the children), and food intake (after 3 (PFOA) or 6 (PFOS) y.o of the children). Similar measured biomarker levels can correspond to large differences in the simulated internal exposures, highlighting the importance to investigate the children's exposure over the early life to improve exposure classification. Our approach demonstrates the possibility to simulate individual internal exposures using PBPK models when measured biomarkers are scarce, helping risk assessors in gaining insight into internal exposure during critical windows, such as early life.

A novel method to derive a human safety limit for PFOA by gene expression profiling and modelling

Perfluorooctanoic acid (PFOA) is a persistent environmental contaminant that can accumulate in the human body due to its long half-life. This substance has been associated with liver, pancreatic, testicular and breast cancers, liver steatosis and endocrine disruption. PFOA is a member of a large group of substances also known as "forever chemicals" and the vast majority of substances of this group lack toxicological data that would enable their effective risk assessment in terms of human health hazards. This study aimed to derive a health-based guidance value for PFOA intake (ng/kg BW/day) from in vitro transcriptomics data. To this end, we developed an in silico workflow comprising five components: (i) sourcing in vitro hepatic transcriptomics concentration-response data; (ii) deriving molecular points of departure using BMDExpress3 and performing pathway analysis using gene set enrichment analysis (GSEA) to identify the most sensitive molecular pathways to PFOA exposure; (iii) estimating freely-dissolved PFOA concentrations in vitro using a mass balance model; (iv) estimating in vivo doses by reverse dosimetry using a PBK model for PFOA as part of a quantitative in vitro to in vivo extrapolation (QIVIVE) algorithm; and (v) calculating a tolerable daily intake (TDI) for PFOA. Fourteen percent of interrogated genes exhibited in vitro concentration-response relationships. GSEA pathway enrichment analysis revealed that "fatty acid metabolism" was the most sensitive pathway to PFOA exposure. In vitro free PFOA concentrations were calculated to be 2.9% of the nominal applied concentrations, and these free concentrations were input into the QIVIVE workflow. Exposure doses for a virtual population of 3,000 individuals were estimated, from which a TDI of 0.15 ng/kg BW/day for PFOA was calculated using the benchmark dose modelling software, PROAST. This TDI is comparable to previously published values of 1.16, 0.69, and 0.86 ng/kg BW/day by the European Food Safety Authority. In conclusion, this study demonstrates the combined utility of an "omics"-derived molecular point of departure and in silico QIVIVE workflow for setting health-based guidance values in anticipation of the acceptance of in vitro concentration-response molecular measurements in chemical risk assessment.

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.

Comparison of aggregated exposure to perfluorooctanoic acid (PFOA) from diet and personal care products with concentrations in blood using a PBPK model - Results from the Norwegian biomonitoring study in EuroMix

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) constitute a large group of compounds that are water, stain, and oil repellent. Numerous sources contribute to the blood levels of PFAS in the European population. The main contributor for perfluorooctanoic acid (PFOA) is food, house dust, consumer products and personal care products (PCPs).

OBJECTIVES

The purpose of the present work is to calculate the dietary and dermal external exposure to PFOA, estimate the aggregated internal exposure from diet and PCPs using a PBPK model, and compare estimates with measured concentrations.

METHODS

Detailed information on diet and PCP use from the EuroMix study is combined with concentration data of PFOA in food and PCPs in a probabilistic exposure assessment. A physiologically based pharmacokinetic model (PBPK) was further refined by incorporating a dermal exposure pathway, and changes in the kidney and faecal excretion.

RESULTS

The aggregated internal exposure using the PBPK model shows that the major contributor to the internal exposure is diet for both males and females. The estimated internal exposure of PFOA for the EuroMix population was in the same range but lower than the measured blood concentrations using the lower bound (LB) external exposure estimates, showing that the LB estimates are underestimations. For seven females the internal exposure of PFOA were higher from PCPs than from diet.

CONCLUSION

PCPs and diet contributed in the same range to the internal PFOA exposure for several women participating in EuroMix. This calls for additional studies on exposure to PFOA and possibly other PFAS from PCPs, especially for women. Overall, PBPK modelling was shown as valuable tool in understanding the sources of PFOA exposure and in guiding risk assessments and regulatory decisions.

Exploring Route-Specific Pharmacokinetics of PFAS in Mice by Coupling in Vivo Tests and Physiologically Based Toxicokinetic Models

BACKGROUND

Oral ingestion, inhalation, and skin contact are important exposure routes for humans to uptake per- and polyfluoroalkyl substances (PFAS). However, nasal and dermal exposure to PFAS remains unclear, and accurately predicting internal body burden of PFAS in humans via multiple exposure pathways is urgently required.

OBJECTIVES

We aimed to develop multiple physiologically based toxicokinetic (PBTK) models to unveil the route-specific pharmacokinetics and bioavailability of PFAS via respective oral, nasal, and dermal exposure pathways using a mouse model and sought to predict the internal concentrations in various tissues through multiple exposure routes and extrapolate it to humans.

METHODS

Mice were administered the mixed solution of perfluorohexane sulfonate, perfluorooctane sulfonate, and perfluorooctanoic acid through oral, nasal, and dermal exposure separately or jointly. The time-dependent concentrations of PFAS in plasma and tissues were determined to calibrate and validate the individual and combined PBTK models, which were applied in single- and repeated-dose scenarios.

RESULTS

The developed route-specific PBTK models successfully simulated the tissue concentrations of PFAS in mice following single or joint exposure routes as well as long-term repeated dose scenarios. The time to peak concentration of PFAS in plasma via dermal exposure was much longer (34.1-83.0 h) than that via nasal exposure (0.960 h). The bioavailability of PFAS via oral exposure was the highest (73.2%-98.0%), followed by nasal (33.9%-66.8%) and dermal exposure (4.59%-7.80%). This model was extrapolated to predict internal levels in human under real environment.

DISCUSSION

Based on these data, we predict the following: PFAS were absorbed quickly via nasal exposure, whereas a distinct hysteresis effect was observed for dermal exposure. Almost all the PFAS to which mice were exposed via gastrointestinal route were absorbed into plasma, which exhibited the highest bioavailability. Exhalation clearance greatly depressed the bioavailability of PFAS via nasal exposure, whereas the lowest bioavailability in dermal exposure was because of the interception of PFAS within the skin layers. https://doi.org/10.1289/EHP11969.

Development and application of an interactive generic physiologically based pharmacokinetic (igPBPK) model for adult beef cattle and lactating dairy cows to estimate tissue distribution and edible tissue and milk withdrawal intervals for per- and polyfluoroalkyl substances (PFAS)

Humans can be exposed to per- and polyfluoroalkyl substances (PFAS) through dietary intake from milk and edible tissues from food animals. This study developed a physiologically based pharmacokinetic (PBPK) model to predict tissue and milk residues and estimate withdrawal intervals (WDIs) for multiple PFAS including PFOA, PFOS and PFHxS in beef cattle and lactating dairy cows. Results showed that model predictions were mostly within a two-fold factor of experimental data for plasma, tissues, and milk with an estimated coefficient of determination (R2) of >0.95. The predicted muscle WDIs for beef cattle were <1 day for PFOA, 449 days for PFOS, and 69 days for PFHxS, while the predicted milk WDIs in dairy cows were <1 day for PFOA, 1345 days for PFOS, and zero day for PFHxS following a high environmental exposure scenario (e.g., 49.3, 193, and 161 ng/kg/day for PFOA, PFOS, and PFHxS, respectively, for beef cattle for 2 years). The model was converted to a web-based interactive generic PBPK (igPBPK) platform to provide a user-friendly dashboard for predictions of tissue and milk WDIs for PFAS in cattle. This model serves as a foundation for extrapolation to other PFAS compounds to improve safety assessment of cattle-derived food products.

Can oral toxicity data for PFAS inform on toxicity via inhalation?

Per- and poly-fluoroalkyl substances (PFAS) are ubiquitous in the environment and are detected in wildlife and humans. With respect to human exposure, studies have shown that ingestion is the primary route of exposure; however, in certain settings, exposure via inhalation could also be a significant source of exposure. While many studies examined toxicity of PFAS via ingestion, limited information is available for PFAS toxicity via the inhalation route, translating into a lack of exposure guidelines. Consequently, this article examined whether route-to-route extrapolation to derive guidelines for inhalation exposure is appropriate for PFAS. Perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS) were used as exemplary PFAS given the abundance of toxicity data for these two compounds. Our evaluation determined that available toxicity and toxicokinetic data support route-to-route extrapolation for PFAS in order to derive inhalation-based standards. Results from this analysis suggest that an air concentration of 7.0 × 10-5  mg/m3 (or 0.07 μg/m3 ) would be an appropriate RfC for PFOA and PFOS assuming the 2016 EPA RfD of 0.00002 mg/kg-day, whereas use of the interim RfDs proposed in 2022 of 1.5 × 10-9 and 7.9 × 10-9  mg/kg would yield much lower RfCs of 5.25 × 10-9 and 2.77 × 10-8  mg/m3 (or 5.25 × 10-6 and 2.77 × 10-5 μg/m3 ) for PFOA and PFOS, respectively.

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.

Folate concentrations and serum perfluoroalkyl and polyfluoroalkyl substance concentrations in adolescents and adults in the USA (National Health and Nutrition Examination Study 2003-16): an observational study

BACKGROUND

Perfluoroalkyl and polyfluoroalkyl substances (PFAS) are a family of highly fluorinated aliphatic compounds, which are widely used in commercial applications, including food packaging, textiles, and non-stick cookware. Folate might counteract the effects of environmental chemical exposures. We aimed to explore the relationship between blood folate biomarker concentrations and PFAS concentrations.

METHODS

This observational study pooled cross-sectional data from the National Health and Nutrition Examination Survey (NHANES) 2003 to 2016 cycles. NHANES is a population-based national survey that measures the health and nutritional status of the US general population every 2 years by means of questionnaires, physical examination, and biospecimen collection. Folate concentrations in red blood cells and in serum, and perfluorooctanoic acid (PFOA), perfluorooctanesulfonic acid (PFOS), perfluorononanoic acid (PFNA), and perfluorohexane sulfonic acid (PFHxS) concentrations in serum were examined. We used multivariable regression models to assess the percentage change in serum PFAS concentrations in relation to changes in folate biomarker concentrations. We additionally used models with restricted cubic splines to investigate the shape of these associations.

FINDINGS

This study included 2802 adolescents and 9159 adults who had complete data on PFAS concentrations, folate biomarkers, and covariates, were not pregnant, and had never had a cancer diagnosis at the time of the survey. The mean age was 15·4 years (SD 2·3) for adolescents and 45·5 years (17·5) for adults. The proportion of male participants was slightly higher in adolescents (1508 [54%] of 2802 participants) than in adults (3940 [49%] of 9159 participants). We found negative associations between red blood cell folate concentrations and serum concentrations of PFOS (percentage change for a 2·7 fold-increase in folate level -24·36%, 95% CI -33·21 to -14·34) and PFNA (-13·00%, -21·87 to -3·12) in adolescents, and PFOA (-12·45%, -17·28 to -7·35), PFOS (-25·30%, -29·67 to -20·65), PFNA (-21·65%, -26·19 to -16·82), and PFHxS (-11·70%, -17·32 to 5·70) in adults. Associations for serum folate concentrations and PFAS were in line with those found for red blood cell folate levels, although the magnitude of the effects was lower. Restricted cubic spline models suggested linearity of the observed associations, particularly for associations in adults.

INTERPRETATION

In this large-scale, nationally representative study, we found consistent inverse associations for most examined serum PFAS compounds in relation to folate concentrations measured in either red blood cells or serum among both adolescents and adults. These findings are supported by mechanistic in-vitro studies that show the potential of PFAS to compete with folate for several transporters implicated in PFAS toxicokinetics. If confirmed in experimental settings, these findings could have important implications for interventions to reduce the accumulated PFAS body burden and mitigate the related adverse health effects.

FUNDING

United States National Institute of Environmental Health Sciences.

Estimating historical exposure to perfluoroalkyl acids in Security, Fountain, and Widefield Colorado: use of water-infrastructure blending and toxicokinetic models

Drinking water can be a major source of poly- and perfluoroalkyl substance (PFAS) exposure for humans. The lack of historic data on PFAS drinking-water concentrations and consumption patterns are a limiting factor for developing estimates of past exposure. Here, in contribution to a community-scale PFAS health effects study near fire training facilities that contaminated a local aquifer with PFASs, we present a novel water-infrastructure, mass-balance mixing model coupled to a non-steady state, single-compartment toxicokinetic model that used Monte Carlo simulations to estimate the start of PFAS exposure in drinking water for individuals within three PFAS-impacted communities in El Paso County, Colorado. Our modeling focused on perfluorohexane sulfonic acid (PFHxS) because median serum PFHxS concentrations in a sample of local residents (n = 213) were twelve times the median observed in the U.S. National Health and Nutrition Examination Survey (2015-2016). Modeling results for study participants were grouped according to their community of residence, revealing a median start of exposure for the town of Fountain of 1998 (25-75% interquartile range [IQR], 1992 to 2010), 2006 (IQR 1995 to 2012) for Security, and 2009 (IQR 1996-2012) for Widefield. Based on the towns' locations relative to an identified hydraulically upgradient PFAS source, the modeled exposure sequencing does not completely align with this conceptual flow model, implying the presence of an additional PFAS source for the groundwater between Widefield and Fountain.

Mechanistic Middle-Out Physiologically Based Toxicokinetic Modeling of Transporter-Dependent Disposition of Perfluorooctanoic Acid in Humans

Perfluorooctanoic acid (PFOA) is an environmental toxicant exhibiting a years-long biological half-life (t_1/2) in humans and is linked with adverse health effects. However, limited understanding of its toxicokinetics (TK) has obstructed the necessary risk assessment. Here, we constructed the first middle-out physiologically based toxicokinetic (PBTK) model to mechanistically explain the persistence of PFOA in humans. In vitro transporter kinetics were thoroughly characterized and scaled up to in vivo clearances using quantitative proteomics-based in vitro-to-in vivo extrapolation. These data and physicochemical parameters of PFOA were used to parameterize our model. We uncovered a novel uptake transporter for PFOA, highly likely to be monocarboxylate transporter 1 which is ubiquitously expressed in body tissues and may mediate broad tissue penetration. Our model was able to recapitulate clinical data from a phase I dose-escalation trial and divergent half-lives from clinical trial and biomonitoring studies. Simulations and sensitivity analyses confirmed the importance of renal transporters in driving extensive PFOA reabsorption, reducing its clearance and augmenting its t_1/2. Crucially, the inclusion of a hypothetical, saturable renal basolateral efflux transporter provided the first unified explanation for the divergent t_1/2 of PFOA reported in clinical (116 days) versus biomonitoring studies (1.3-3.9 years). Efforts are underway to build PBTK models for other perfluoroalkyl substances using similar workflows to assess their TK profiles and facilitate risk assessments.

Association of Early Pregnancy Perfluoroalkyl and Polyfluoroalkyl Substance Exposure With Birth Outcomes

Importance

Prenatal perfluoroalkyl and polyfluoroalkyl substances (PFAS) have been linked to adverse birth outcomes. Previous research showed that higher folate concentrations are associated with lower blood PFAS concentrations in adolescents and adults. Further studies are needed to explore whether prenatal folate status mitigates PFAS-related adverse birth outcomes.

Objective

To examine whether prenatal folate status modifies the negative associations between pregnancy PFAS concentrations, birth weight, and gestational age previously observed in a US cohort.

Design, Setting, and Participants

In a prospective design, a prebirth cohort of mothers or pregnant women was recruited between April 1999 and November 2002, in Project Viva, a study conducted in eastern Massachusetts. Statistical analyses were performed from May 24 and October 25, 2022.

Exposure

Plasma concentrations of 6 PFAS compounds were measured in early pregnancy (median gestational week, 9.6). Folate status was assessed through a food frequency questionnaire and measured in plasma samples collected in early pregnancy.

Main Outcomes and Measures

Birth weight and gestational age, abstracted from delivery records; birth weight z score, standardized by gestational age and infant sex; low birth weight, defined as birth weight less than 2500 g; and preterm birth, defined as birth at less than 37 completed gestational weeks.

Results

The cohort included a total of 1400 mother-singleton pairs. The mean (SD) age of the mothers was 32.21 (4.89) years. Most of the mothers were White (73.2%) and had a college degree or higher (69.1%). Early pregnancy plasma perfluorooctanoic acid concentration was associated with lower birth weight and birth weight z score only among mothers whose dietary folate intake (birth weight: β, -89.13 g; 95% CI, -166.84 to -11.42 g; birth weight z score: -0.13; 95% CI, -0.26 to -0.003) or plasma folate concentration (birth weight: -87.03 g; 95% CI, -180.11 to 6.05 g; birth weight z score: -0.14; 95% CI, -0.30 to 0.02) were below the 25th percentile (dietary: 660 μg/d, plasma: 14 ng/mL). No associations were found among mothers in the higher folate level groups, although the tests for heterogeneity did not reject the null. Associations between plasma perfluorooctane sulfonic acid and perfluorononanoate (PFNA) concentrations and lower birth weight, and between PFNA and earlier gestational age were noted only among mothers whose prenatal dietary folate intake or plasma folate concentration was in the lowest quartile range. No associations were found among mothers in higher folate status quartile groups.

Conclusions and Relevance

In this large, US prebirth cohort, early pregnancy exposure to select PFAS compounds was associated with adverse birth outcomes only among mothers below the 25th percentile of prenatal dietary or plasma folate levels.

Human health risk assessment of 6:2 Cl-PFESA through quantitative in vitro to in vivo extrapolation by integrating cell-based assays, an epigenetic key event, and physiologically based pharmacokinetic modeling

Human health risk assessment of chemicals is essential but often relies on time-consuming and animal and labor-extensive procedures. Here, we develop a population-based, quantitative in vitro to in vivo extrapolation (QIVIVE) approach which depended on cellular effects monitored by in vitro assays, considered chemical internal concentration determined by LC-MS/MS, extrapolated into in vivo target tissue concentration through physiologically based pharmacokinetic (PBPK) modelling, and assessed populational health risk using in silico modelling. By applying this QIVIVE approach to 6:2 chlorinated polyfluorinated ether sulfonate (6:2 Cl-PFESA), as a representative of the emerging pollutants, we find that 6:2 Cl-PFESA disturbed lipid homeostasis in HepG2 cells through enhancement of lipid accumulation and fatty acid β-oxidation, during which miR-93-5p served as a key event towards toxicity and thus, could serve as an efficient toxicity marker for risk assessment; further, the disruption potency of lipid homeostasis of 6:2 Cl-PFESA for the most of studied populations in China might be of moderate concern. Together, our approach improved the reliability of QIVIVE during human health risk assessment, which can readily be used for other chemicals.

Red Blood Cell Folate Modifies the Association between Serum Per- and Polyfluoroalkyl Substances and Antibody Concentrations in U.S. Adolescents

Per- and polyfluoroalkyl substances (PFAS) exposure has been associated with reduced antibody levels. Higher red blood cell (RBC) folate was previously associated with lower serum PFAS concentrations in adolescents. This study included 819 adolescents aged 12-19 years who had detectable rubella and measles antibody levels in serum from the U.S. National Health and Nutrition Examination Survey 2003-2004 and 2009-2010 cycles. We found inverse associations between serum PFOS and PFHxS and rubella antibodies, between PFOA and mumps antibodies, and between PFAS mixtures and rubella and mumps antibodies, only among adolescents with RBC folate concentrations <66th percentile (lower folate group) while not among adolescents with higher RBC folate levels (upper folate group). Specifically, per quartile increase in serum concentrations of the total PFAS mixture was associated with a 9.84% (95% CI: -15.57%, -3.74%) decrease in rubella antibody and an 8.79% (95% CI: -14.39%, -2.82%) decrease in the mumps antibody concentrations only in the lower folate group, while null associations were found for the upper folate group. If confirmed in mechanistic studies or prospective epidemiologic studies, these findings may have important implications for using folate as a mitigation measure against immune-related PFAS effects.

An Overview of Physiologically-Based Pharmacokinetic Models for Forensic Science

A physiologically-based pharmacokinetic (PBPK) model represents the structural components of the body with physiologically relevant compartments connected via blood flow rates described by mathematical equations to determine drug disposition. PBPK models are used in the pharmaceutical sector for drug development, precision medicine, and the chemical industry to predict safe levels of exposure during the registration of chemical substances. However, one area of application where PBPK models have been scarcely used is forensic science. In this review, we give an overview of PBPK models successfully developed for several illicit drugs and environmental chemicals that could be applied for forensic interpretation, highlighting the gaps, uncertainties, and limitations.

The association of perfluoroalkyl substances (PFAS) exposure and kidney function in Korean adolescents using data from Korean National Environmental Health Survey (KoNEHS) cycle 4 (2018–2020): a cross-sectional study

Background

Perfluoroalkyl substances (PFAS) are chemicals widely used in various products in everyday life. Due to its unique strong binding force, the half-life of PFAS is very long, so bioaccumulation and toxicity to the human body are long-standing concerns. In particular, effects on kidney function have recently emerged and there are no studies on the effect of PFAS on kidney function through epidemiological investigations in Korea. From 2018 to 2020, the Korean National Environmental Health Survey (KoNEHS) cycle 4, conducted an epidemiological investigation on the blood concentration of PFAS for the first time in Korea. Based on this data, the relationship between PFAS blood concentration and kidney function was analyzed for adolescents.

Methods

We investigated 5 types of PFAS and their total blood concentration in 811 middle and high school students, living in Korea and included in KoNEHS cycle 4, and tried to find changes in kidney function in relation to PFAS concentration. After dividing the concentration of each of the 5 PFAS and the total concentration into quartiles, multivariable linear regression was performed to assess the correlation with kidney function. The bedside Schwartz equation was used as an indicator of kidney function.

Results

As a result of multivariable linear regression, when observing a change in kidney function according to the increase in the concentration of each of the 5 PFAS and their total, a significant decrease in kidney function was confirmed in some or all quartiles.

Conclusions

In this cross-sectional study of Korean adolescents based on KoNEHS data, a negative correlation between serum PFAS concentration and kidney function was found. A well-designed longitudinal study and continuous follow-up are necessary.

Bayesian Estimation of Human Population Toxicokinetics of PFOA, PFOS, PFHxS, and PFNA from Studies of Contaminated Drinking Water

BACKGROUND

Setting health-protective standards for poly- and perfluoroalkyl substances (PFAS) exposure requires estimates of their population toxicokinetics, but existing studies have reported widely varying PFAS half-lives (T_½) and volumes of distribution (V_d).

OBJECTIVES

We combined data from multiple studies to develop harmonized estimates of T_½ and V_d, along with their interindividual variability, for four PFAS commonly found in drinking water: perfluorooctanoic acid (PFOA), perfluorooctane sulfonate (PFOS), perfluorononanoic acid (PFNA), and perfluorohexane sulfonate (PFHxS).

METHODS

We identified published data on PFAS concentrations in human serum with corresponding drinking water measurements, separated into training and testing data sets. We fit training data sets to a one-compartment model incorporating interindividual variability, time-dependent drinking water concentrations, and background exposures. Use of a hierarchical Bayesian approach allowed us to incorporate informative priors at the population level, as well as at the study level. We compared posterior predictions to testing data sets to evaluate model performance.

RESULTS

Posterior median (95% CI) estimates of T_½ (in years) for the population geometric mean were 3.14 (2.69, 3.73) for PFOA, 3.36 (2.52, 4.42) for PFOS, 2.35 (1.65, 3.16) for PFNA, and 8.30 (5.38, 13.5) for PFHxS, all of which were within the range of previously published values. The extensive individual-level data for PFOA allowed accurate estimation of population variability, with a population geometric standard deviation of 1.57 (95% CI: 1.42, 1.73); data from other PFAS were also consistent with this degree of population variability. V_d estimates ranged from 0.19 to 0.43L/kg across the four PFAS, which tended to be slightly higher than previously published estimates.

DISCUSSION

These results have direct application in both risk assessment (quantitative interspecies extrapolation and uncertainty factors for interindividual variability) and risk communication (interpretation of monitoring data). In addition, this study provides a rigorous methodology for further refinement with additional data, as well as application to other PFAS. https://doi.org/10.1289/EHP10103.

Biomonitoring of per- and polyfluoroalkyl substances in minority angler communities in central New York State

Onondaga Lake in central New York State was listed as a Superfund site in 1994 due to industrial disposal of pollutants. A biomonitoring program was conducted to assess exposure to over 70 legacy contaminants and contaminants of emerging concern in populations disproportionately at risk for exposure residing near Onondaga Lake and to educate these communities on how to reduce exposures. The populations of focus were refugees from Burma and Bhutan and low-income, primarily African American, anglers (urban anglers). These communities consume locally caught fish for economic as well as cultural reasons and therefore may be at higher risk of exposure. This study focuses on assessment of exposure to per- and polyfluoroalkyl substances (PFAS) and associations with local fish consumption. Using respondent driven sampling, 311 refugees and 89 urban anglers were enrolled in the study. Following informed consent, study participants provided blood and urine specimens and completed a questionnaire. Percentiles of locally caught fish meals in the past 12 months by race/ethnicity groups showed that the Burmese participants of Karen ethnicity were the highest consumers, with a median of 135 meals compared to 103 meals for the other Burmese participants, 70 meals for the urban anglers, and 44 meals for the Bhutanese participants. Compared to the National Health and Nutrition Examination Survey (NHANES) 2015-16 sample of the general U.S. population, the Karen participants had markedly elevated perfluorooctane sulfonic acid (PFOS) and perfluorodecanoic acid (PFDA) levels with median serum concentrations 9.5 times greater (41.6 ng/mL vs. 4.4 ng/mL) and 26.9 times greater (2.69 ng/mL vs. 0.10 ng/mL), respectively; the other Burmese participants had moderately elevated levels of PFOS and PFDA with median serum concentrations 3.0 times greater (13.3 ng/mL vs. 4.4 ng/mL) and 7.3 greater times greater (0.73 ng/mL vs. 0.10 ng/mL), respectively; and, PFAS levels were not elevated in the Bhutanese or urban angler cohorts. Male gender was consistently the strongest predictor of PFAS exposure among all study cohorts. A positive association between local fish consumption was indicated only for PFOS among urban anglers. An association between local fish consumption and PFAS was not statistically significant among the refugee cohorts, perhaps due to the lack of 'lower-end' exposure or exposure variability. Community events were held by the program staff to present the biomonitoring results and distribute community outreach materials with visual aids specific for the study populations to promote safe fish eating.

Physiologically based pharmacokinetic (PBPK) modeling of perfluorohexane sulfonate (PFHxS) in humans

Per- and polyfluoroalkyl substances (PFAS) are persistent, man-made compounds prevalent in the environment and consistently identified in human biomonitoring samples. In particular, perfluorooctane sulfonic acid (PFOS), perfluorooctanoic acid (PFOA), and perfluorohexane sulfonic acid (PFHxS) have been identified at U.S. Air Force installations. The study of human toxicokinetics and physiologically based pharmacokinetic (PBPK) modeling of PFHxS has been less robust and has been limited in scope and application as compared to PFOS and PFOA. The primary goal of the current effort was to develop a PBPK model describing PFHxS disposition in humans that can be applied to retrospective, current, and future human health risk assessment of PFHxS. An existing model developed for PFOS and PFOA was modified and key parameter values for exposure and toxicokinetics were calibrated for PFHxS prediction based on human biomonitoring data, particularly general population serum levels from the U.S. Centers for Disease Prevention and Control (CDC) National Health and Nutrition Examination Survey (NHANES). Agreement between the model and the calibration and evaluation data was excellent and recapitulated observed trends across sex, age, and calendar years. Confidence in the model is greatest for application to adults in the 2000-2018 time frame and for shorter-term future projections.

Advancing per- and polyfluoroalkyl substances (PFAS) research: an overview of ATSDR and NCEH activities and recommendations

The National Center for Environmental Health (NCEH), part of the Centers for Disease Control and Prevention (CDC), and the Agency for Toxic Substances and Disease Registry (ATSDR) support and conduct research advancing national, state, and local public health response to per- and polyfluoroalkyl substances (PFAS). PFAS are a group of manufactured chemicals used in industry and consumer products that persist in the environment. Given the growing evidence linking PFAS with adverse health effects in humans, NCEH and ATSDR developed a public health research framework to capture the broad range of PFAS research activities being conducted and supported by the agency to determine future research priorities and identify opportunities for interagency collaboration. The framework was conceptualized via a multidisciplinary visioning process designed to identify compelling questions and research activities that span five scientific domains: toxicology, exposure, human health, public health action, and cross-cutting priorities. This paper presents a framework, compelling questions and research activities to help NCEH and ATSDR advance scientific discovery in partnership with federal, state, and local stakeholders as part of a comprehensive public health response to PFAS contamination.

Why is elevation of serum cholesterol associated with exposure to perfluoroalkyl substances (PFAS) in humans? A workshop report on potential mechanisms

Serum concentrations of cholesterol are positively correlated with exposure to perfluorooctanoic acid (PFOA) and perfluorooctane sulfonic acid (PFOS) in humans. The associated change in cholesterol is small across a broad range of exposure to PFOA and PFOS. Animal studies generally have not indicated a mechanism that would account for the association in humans. The extent to which the relationship is causal is an open question. Nonetheless, the association is of particular importance because increased serum cholesterol has been considered as an endpoint to derive a point of departure in at least one recent risk assessment. To gain insight into potential mechanisms for the association, both causal and non-causal, an expert workshop was held Oct 31 and Nov 1, 2019 to discuss relevant data and propose new studies. In this report, we summarize the relevant background data, the discussion among the attendees, and their recommendations for further research.

Derivation of a Human In Vivo Benchmark Dose for Perfluorooctanoic Acid From ToxCast In Vitro Concentration-Response Data Using a Computational Workflow for Probabilistic Quantitative In Vitro to In Vivo Extrapolation

A computational workflow which integrates physiologically based kinetic (PBK) modeling, global sensitivity analysis (GSA), approximate Bayesian computation (ABC), and Markov Chain Monte Carlo (MCMC) simulation was developed to facilitate quantitative in vitro to in vivo extrapolation (QIVIVE). The workflow accounts for parameter and model uncertainty within a computationally efficient framework. The workflow was tested using a human PBK model for perfluorooctanoic acid (PFOA) and high throughput screening (HTS) in vitro concentration–response data, determined in a human liver cell line, from the ToxCast/Tox21 database. In vivo benchmark doses (BMDs) for PFOA intake (ng/kg BW/day) and drinking water exposure concentrations (µg/L) were calculated from the in vivo dose responses and compared to intake values derived by the European Food Safety Authority (EFSA). The intake benchmark dose lower confidence limit (BMDL₅) of 0.82 was similar to 0.86 ng/kg BW/day for altered serum cholesterol levels derived by EFSA, whereas the intake BMDL₅ of 6.88 was six-fold higher than the value of 1.14 ng/kg BW/day for altered antibody titer also derived by the EFSA. Application of a chemical-specific adjustment factor (CSAF) of 1.4, allowing for inter-individual variability in kinetics, based on biological half-life, gave an intake BMDL₅ of 0.59 for serum cholesterol and 4.91 (ng/kg BW/day), for decreased antibody titer, which were 0.69 and 4.31 the EFSA-derived values, respectively. The corresponding BMDL₅ for drinking water concentrations, for estrogen receptor binding activation associated with breast cancer, pregnane X receptor binding associated with altered serum cholesterol levels, thyroid hormone receptor α binding leading to thyroid disease, and decreased antibody titer (pro-inflammation from cytokines) were 0.883, 0.139, 0.086, and 0.295 ng/ml, respectively, with application of no uncertainty factors. These concentrations are 5.7-, 36-, 58.5-, and 16.9-fold lower than the median measured drinking water level for the general US population which is approximately, 5 ng/ml.

PFAS soil and groundwater contamination via industrial airborne emission and land deposition in SW Vermont and Eastern New York State, USA

In order to understand the extent to which airborne PFAS emission can impact soil and groundwater, we conducted a sampling campaign in areas of conserved forest lands near Bennington, VT/Hoosick Falls, NY. This has been home to sources of PFAS air-emissions from Teflon-coating operations for over 50 years. Since 2015, the Vermont and New York Departments of Environmental Conservation have documented ∼1200 residential wells and two municipal water systems across a 200 km2 area contaminated with perfluorooctanoic acid (PFOA). Given the large areal extent of the plume, and the fact that much of the contaminated area lies up-gradient and across rivers from manufactures, we seek to determine if groundwater contamination could have resulted primarily from air-emission, land deposition, and subsequent leaching to infiltrating groundwater. Sampling of soils and groundwater in the Green Mountain National Forest (GMNF) downwind of factories shows that both soil and groundwater PFOA contamination extend uninterrupted from inhabited areas into conserved forest lands. Groundwater springs and seeps in the GMNF located 8 km downwind, but >300 meters vertically above factories, contain up to 100 ppt PFOA. Our results indicate that air-emitted PFAS can contaminate groundwater and soil in areas outside of those normally considered down-gradient of a source with respect to regional groundwater flow.

Per- and polyfluoroalkyl substances and kidney function: Follow-up results from the Diabetes Prevention Program trial

Per- and polyfluoroalkyl substances (PFAS) are ubiquitously detected in populations worldwide and may hinder kidney function. The objective of the study was to determine longitudinal associations of plasma PFAS concentrations with estimated glomerular filtration rate (eGFR) and evaluate whether a lifestyle intervention modify the associations. We studied 875 participants initially randomized to the lifestyle or placebo arms in the Diabetes Prevention Program (DPP, 1996-2002) trial and Outcomes Study (DPPOS, 2002-2014). We ran generalized linear mixed models accounting a priori covariates to evaluate the associations between baseline PFAS concentrations and repeated measures of eGFR, separately, for six PFAS (PFOS, PFOA, PFHxS, EtFOSAA, MeFOSAA, PFNA); then used quantile-based g-computation to evaluate the effects of the six PFAS chemicals as a mixture. The cohort was 64.9% female; 73.4% 40-64 years-old; 29.4% with hypertension; 50.5% randomized to lifestyle intervention and 49.5% to placebo and had similar plasma PFAS concentrations as the general U.S. population in 1999-2000. Most participants had normal kidney function (eGFR > 90 mL/min/1.73 m²) over the approximately 14 years of follow-up. We found that plasma PFAS concentrations during DPP were inversely associated with eGFR during DPPOS follow-up. Each quartile increase in baseline plasma concentration of the 6 PFAS as a mixture was associated with 2.26 mL/min/1.73 m² lower eGFR (95% CI: -4.12, -0.39) at DPPOS Year 5, approximately 9 years since DPP randomization and PFAS measurements. The lifestyle intervention did not modify associations, but inverse associations were stronger among participants with hypertension at baseline. Among prediabetic adults, we found inverse associations between baseline plasma PFAS concentrations and measures of eGFR throughout 14 years of follow-up. The lifestyle intervention of diet, exercise and behavioral changes did not modify the associations, but persons with hypertension may have heightened susceptibility.

Development of a Gestational and Lactational Physiologically Based Pharmacokinetic (PBPK) Model for Perfluorooctane Sulfonate (PFOS) in Rats and Humans and Its Implications in the Derivation of Health-Based Toxicity Values

BACKGROUND

There is a great concern on potential adverse effects of exposure to perfluorooctane sulfonate (PFOS) in sensitive subpopulations, such as pregnant women, fetuses, and neonates, due to its reported transplacental and lactational transfer and reproductive and developmental toxicities in animals and humans.

OBJECTIVES

This study aimed to develop a gestational and lactational physiologically based pharmacokinetic (PBPK) model in rats and humans for PFOS to aid risk assessment in sensitive human subpopulations.

METHODS

Based upon existing PBPK models for PFOS, the present model addressed a data gap of including a physiologically based description of basolateral and apical membrane transporter-mediated renal reabsorption and excretion in kidneys during gestation and lactation. The model was calibrated with published rat toxicokinetic and human biomonitoring data and was independently evaluated with separate data. Monte Carlo simulation was used to address the interindividual variability.

RESULTS

Model simulations were generally within 2-fold of observed PFOS concentrations in maternal/fetal/neonatal plasma and liver in rats and humans. Estimated fifth percentile human equivalent doses (HEDs) based on selected critical toxicity studies in rats following U.S. Environmental Protection Agency (EPA) guidelines ranged from 0.08 to 0.91 μ g / kg  per day . These values are lower than the HEDs estimated in U.S. EPA guidance (0.51 - 1.6 μ g / kg  per day ) using an empirical toxicokinetic model in adults.

CONCLUSIONS

The results support the importance of renal reabsorption/excretion during pregnancy and lactation in PFOS dosimetry and suggest that the derivation of health-based toxicity values based on developmental toxicity studies should consider gestational/lactational dosimetry estimated from a life stage-appropriate PBPK model. This study provides a quantitative tool to aid risk reevaluation of PFOS, especially in sensitive human subpopulations, and it provides a basis for extrapolating to other per- and polyfluoroalkyl substances (PFAS). All model codes and detailed tutorials are provided in the Supplemental Materials to allow readers to reproduce our results and to use this model. https://doi.org/10.1289/EHP7671.

Early life exposure to per- and polyfluoroalkyl substances (PFAS) and latent health outcomes: A review including the placenta as a target tissue and possible driver of peri- and postnatal effects

Per- and polyfluoroalkyl substances (PFAS) are ubiquitous drinking water contaminants of concern due to mounting evidence implicating adverse health outcomes associated with exposure, including reduced kidney function, metabolic syndrome, thyroid disruption, and adverse pregnancy outcomes. PFAS have been produced in the U.S. since the 1940s and now encompass a growing chemical family comprised of diverse chemical moieties, yet the toxicological effects have been studied for relatively few compounds. Critically, exposures to some PFAS in utero are associated with adverse outcomes for both mother and offspring, such as hypertensive disorders of pregnancy (HDP), including preeclampsia, and low birth weight. Given the relationship between HDP, placental dysfunction, adverse health outcomes, and increased risk for chronic diseases in adulthood, the role of both developmental and lifelong exposure to PFAS likely contributes to disease risk in complex ways. Here, evidence for the role of some PFAS in disrupted thyroid function, kidney disease, and metabolic syndrome is synthesized with an emphasis on the placenta as a critical yet understudied target of PFAS and programming agent of adult disease. Future research efforts must continue to fill the knowledge gap between placental susceptibility to environmental exposures like PFAS, subsequent perinatal health risks for both mother and child, and latent health effects in adult offspring.

Derivation of a chronic reference dose for perfluorohexane sulfonate (PFHxS) for reproductive toxicity in mice

Perfluorohexane sulfonate (PFHxS) is a six-carbon perfluoroalkyl sulfonic acid that was used as an industrial surfactant, but is now found as an environmental contaminant worldwide. In addition to its use as an industrial surfactant, it is a legacy contaminant from the use of aqueous film-forming foams. Despite its widespread occurrence in the environment and evidence of biological activity associated with PFHxS and similar perfluoroalkyl sulfonic acids in rodents, there is no oral toxicity value currently available from the IRIS Database. To derive an oral reference dose (RfD) for PFHxS, available toxicity studies were reviewed using a weight-of-evidence approach. A 42-day mouse reproductive study was chosen as the critical study for the derivation of the oral RfD. Benchmark dose modeling was utilized to derive a point of departure (POD) for a reduction in litter size. A 95% lower confidence limit on the benchmark dose (BMDL) of 13,900 ng/mL (serum PFHxS) was modeled for a reduction in litter size. An oral RfD for PFHxS of 4.0 ng/kg/d was calculated by conversion of the BMDL to a human equivalent oral dose using a human half-life adjusted dosimetric conversion factor and the application of a total uncertainty factor of 300. Additional research is needed to better characterize the toxicity associated with oral exposure to PFHxS and refine the development of toxicity values.

Bayesian evaluation of a physiologically based pharmacokinetic (PBPK) model for perfluorooctane sulfonate (PFOS) to characterize the interspecies uncertainty between mice, rats, monkeys, and humans: Development and performance verification

A challenge in the risk assessment of perfluorooctane sulfonate (PFOS) is the large interspecies differences in its toxicokinetics that results in substantial uncertainty in the dosimetry and toxicity extrapolation from animals to humans. To address this challenge, the objective of this study was to develop an open-source physiologically based pharmacokinetic (PBPK) model accounting for species-specific toxicokinetic parameters of PFOS. Considering available knowledge about the toxicokinetic properties of PFOS, a PBPK model for PFOS in mice, rats, monkeys, and humans after intravenous and oral administrations was created. Available species-specific toxicokinetic data were used for model calibration and optimization, and independent datasets were used for model evaluation. Bayesian statistical analysis using Markov chain Monte Carlo (MCMC) simulation was performed to optimize the model and to characterize the uncertainty and interspecies variability of chemical-specific parameters. The model predictions well correlated with the majority of datasets for all four species, and the model was validated with independent data in rats, monkeys, and humans. The model was applied to predict human equivalent doses (HEDs) based on reported points of departure in selected critical toxicity studies in rats and monkeys following U.S. EPA's guidelines. The lower bounds of the model-derived HEDs were overall lower than the HEDs estimated by U.S. EPA (e.g., 0.2 vs. 1.3 μg/kg/day based on the rat plasma data). This integrated and comparative analysis provides an important step towards improving interspecies extrapolation and quantitative risk assessment of PFOS, and this open-source model provides a foundation for developing models for other perfluoroalkyl substances.

Prediction of maternal and foetal exposures to perfluoroalkyl compounds in a Spanish birth cohort using toxicokinetic modelling

Prenatal exposures to perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) have been associated with child health outcomes, but many of these associations remain poorly characterized. The aim of this work was to provide new indicators of foetal exposure for the Spanish INMA birth cohort. First, a pregnancy and lactation physiologically based pharmacokinetic (PBPK) model was calibrated in a population framework to provide quantitative estimates for the PFOA and PFOS placental transfers in humans. The estimated distributions indicated that PFOA crosses the placental barrier at a rate three times higher than PFOS and shows a higher variability between mothers. The PBPK model was then used to back-calculate the time-varying daily intakes of the INMA mothers corrected for their individual history from a spot maternal concentration. We showed the importance of accounting for the mothers' history as different dietary intakes can result in similar measured concentrations at one time point. Finally, the foetal exposure was simulated in target organs over pregnancy using the PBPK model and the estimated maternal intakes. We showed that the pattern of PFOA and PFOS exposures varies greatly among the foetuses. About a third has levels of either one compound always higher than the levels of the other compound. The other two thirds showed different ranking of PFOA and PFOS in terms of concentrations in the target organs. Our simulated foetal exposures bring additional information to the measured maternal spot concentrations and can help to better characterize the prenatal exposure in target organs during windows of susceptibility.

Challenges Associated With Applying Physiologically Based Pharmacokinetic Modeling for Public Health Decision-Making

The development and application of physiologically based pharmacokinetic (PBPK) models in chemical toxicology have grown steadily since their emergence in the 1980s. However, critical evaluation of PBPK models to support public health decision-making across federal agencies has thus far occurred for only a few environmental chemicals. In order to encourage decision-makers to embrace the critical role of PBPK modeling in risk assessment, several important challenges require immediate attention from the modeling community. The objective of this contemporary review is to highlight 3 of these challenges, including: (1) difficulties in recruiting peer reviewers with appropriate modeling expertise and experience; (2) lack of confidence in PBPK models for which no tissue/plasma concentration data exist for model evaluation; and (3) lack of transferability across modeling platforms. Several recommendations for addressing these 3 issues are provided to initiate dialog among members of the PBPK modeling community, as these issues must be overcome for the field of PBPK modeling to advance and for PBPK models to be more routinely applied in support of public health decision-making.

Risk to human health related to the presence of perfluorooctane sulfonic acid and perfluorooctanoic acid in food

The European Commission asked EFSA for a scientific evaluation on the risks to human health related to the presence of perfluorooctane sulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) in food. Regarding PFOS and PFOA occurrence, the final data set available for dietary exposure assessment contained a total of 20,019 analytical results (PFOS n = 10,191 and PFOA n = 9,828). There were large differences between upper and lower bound exposure due to analytical methods with insufficient sensitivity. The CONTAM Panel considered the lower bound estimates to be closer to true exposure levels. Important contributors to the lower bound mean chronic exposure were 'Fish and other seafood', 'Meat and meat products' and 'Eggs and egg products', for PFOS, and 'Milk and dairy products', 'Drinking water' and 'Fish and other seafood' for PFOA. PFOS and PFOA are readily absorbed in the gastrointestinal tract, excreted in urine and faeces, and do not undergo metabolism. Estimated human half-lives for PFOS and PFOA are about 5 years and 2-4 years, respectively. The derivation of a health-based guidance value was based on human epidemiological studies. For PFOS, the increase in serum total cholesterol in adults, and the decrease in antibody response at vaccination in children were identified as the critical effects. For PFOA, the increase in serum total cholesterol was the critical effect. Also reduced birth weight (for both compounds) and increased prevalence of high serum levels of the liver enzyme alanine aminotransferase (ALT) (for PFOA) were considered. After benchmark modelling of serum levels of PFOS and PFOA, and estimating the corresponding daily intakes, the CONTAM Panel established a tolerable weekly intake (TWI) of 13 ng/kg body weight (bw) per week for PFOS and 6 ng/kg bw per week for PFOA. For both compounds, exposure of a considerable proportion of the population exceeds the proposed TWIs.

Risk to human health related to the presence of perfluorooctane sulfonic acid and perfluorooctanoic acid in food

The European Commission asked EFSA for a scientific evaluation on the risks to human health related to the presence of perfluorooctane sulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) in food. Regarding PFOS and PFOA occurrence, the final data set available for dietary exposure assessment contained a total of 20,019 analytical results (PFOS n = 10,191 and PFOA n = 9,828). There were large differences between upper and lower bound exposure due to analytical methods with insufficient sensitivity. The CONTAM Panel considered the lower bound estimates to be closer to true exposure levels. Important contributors to the lower bound mean chronic exposure were 'Fish and other seafood', 'Meat and meat products' and 'Eggs and egg products', for PFOS, and 'Milk and dairy products', 'Drinking water' and 'Fish and other seafood' for PFOA. PFOS and PFOA are readily absorbed in the gastrointestinal tract, excreted in urine and faeces, and do not undergo metabolism. Estimated human half-lives for PFOS and PFOA are about 5 years and 2-4 years, respectively. The derivation of a health-based guidance value was based on human epidemiological studies. For PFOS, the increase in serum total cholesterol in adults, and the decrease in antibody response at vaccination in children were identified as the critical effects. For PFOA, the increase in serum total cholesterol was the critical effect. Also reduced birth weight (for both compounds) and increased prevalence of high serum levels of the liver enzyme alanine aminotransferase (ALT) (for PFOA) were considered. After benchmark modelling of serum levels of PFOS and PFOA, and estimating the corresponding daily intakes, the CONTAM Panel established a tolerable weekly intake (TWI) of 13 ng/kg body weight (bw) per week for PFOS and 6 ng/kg bw per week for PFOA. For both compounds, exposure of a considerable proportion of the population exceeds the proposed TWIs.

Perfluorinated Chemicals as Emerging Environmental Threats to Kidney Health: A Scoping Review

BACKGROUND AND OBJECTIVES

Per- and polyfluoroalkyl substances (PFASs) are a large group of manufactured nonbiodegradable compounds. Despite increasing awareness as global pollutants, the impact of PFAS exposure on human health is not well understood, and there are growing concerns for adverse effects on kidney function. Therefore, we conducted a scoping review to summarize and identify gaps in the understanding between PFAS exposure and kidney health.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

We systematically searched PubMed, EMBASE, EBSCO Global Health, World Health Organization Global Index, and Web of Science for studies published from 1990 to 2018. We included studies on the epidemiology, pharmacokinetics, or toxicology of PFAS exposure and kidney-related health, including clinical, histologic, molecular, and metabolic outcomes related to kidney disease, or outcomes related to the pharmacokinetic role of the kidneys.

RESULTS

We identified 74 studies, including 21 epidemiologic, 13 pharmacokinetic, and 40 toxicological studies. Three population-based epidemiologic studies demonstrated associations between PFAS exposure and lower kidney function. Along with toxicology studies (n=10) showing tubular histologic and cellular changes from PFAS exposure, pharmacokinetic studies (n=5) demonstrated the kidneys were major routes of elimination, with active proximal tubule transport. In several studies (n=17), PFAS exposure altered several pathways linked to kidney disease, including oxidative stress pathways, peroxisome proliferators-activated receptor pathways, NF-E2-related factor 2 pathways, partial epithelial mesenchymal transition, and enhanced endothelial permeability through actin filament modeling.

CONCLUSIONS

A growing body of evidence portends PFASs are emerging environmental threats to kidney health; yet several important gaps in our understanding still exist.

Perfluorinated Chemicals as Emerging Environmental Threats to Kidney Health: A Scoping Review

BACKGROUND AND OBJECTIVES

Per- and polyfluoroalkyl substances (PFASs) are a large group of manufactured nonbiodegradable compounds. Despite increasing awareness as global pollutants, the impact of PFAS exposure on human health is not well understood, and there are growing concerns for adverse effects on kidney function. Therefore, we conducted a scoping review to summarize and identify gaps in the understanding between PFAS exposure and kidney health.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

We systematically searched PubMed, EMBASE, EBSCO Global Health, World Health Organization Global Index, and Web of Science for studies published from 1990 to 2018. We included studies on the epidemiology, pharmacokinetics, or toxicology of PFAS exposure and kidney-related health, including clinical, histologic, molecular, and metabolic outcomes related to kidney disease, or outcomes related to the pharmacokinetic role of the kidneys.

RESULTS

We identified 74 studies, including 21 epidemiologic, 13 pharmacokinetic, and 40 toxicological studies. Three population-based epidemiologic studies demonstrated associations between PFAS exposure and lower kidney function. Along with toxicology studies (n=10) showing tubular histologic and cellular changes from PFAS exposure, pharmacokinetic studies (n=5) demonstrated the kidneys were major routes of elimination, with active proximal tubule transport. In several studies (n=17), PFAS exposure altered several pathways linked to kidney disease, including oxidative stress pathways, peroxisome proliferators-activated receptor pathways, NF-E2-related factor 2 pathways, partial epithelial mesenchymal transition, and enhanced endothelial permeability through actin filament modeling.

CONCLUSIONS

A growing body of evidence portends PFASs are emerging environmental threats to kidney health; yet several important gaps in our understanding still exist.

Can profiles of poly- and Perfluoroalkyl substances (PFASs) in human serum provide information on major exposure sources?

BACKGROUND

Humans are exposed to poly- and perfluoroalkyl substances (PFASs) from diverse sources and this has been associated with negative health impacts. Advances in analytical methods have enabled routine detection of more than 15 PFASs in human sera, allowing better profiling of PFAS exposures. The composition of PFASs in human sera reflects the complexity of exposure sources but source identification can be confounded by differences in toxicokinetics affecting uptake, distribution, and elimination. Common PFASs, such as perfluorooctanoic acid (PFOA), perfluorooctane sulfonic acid (PFOS) and their precursors are ubiquitous in multiple exposure sources. However, their composition varies among sources, which may impact associated adverse health effects.

METHODS

We use available PFAS concentrations from several demographic groups in a North Atlantic seafood consuming population (Faroe Islands) to explore whether chemical fingerprints in human sera provide insights into predominant exposure sources. We compare serum PFAS profiles from Faroese individuals to other North American populations to investigate commonalities in potential exposure sources. We compare individuals with similar demographic and physiological characteristics and samples from the same years to reduce confounding by toxicokinetic differences and changing environmental releases.

RESULTS

Using principal components analysis (PCA) confirmed by hierarchical clustering, we assess variability in serum PFAS concentrations across three Faroese groups. The first principal component (PC)/cluster consists of C9-C12 perfluoroalkyl carboxylates (PFCAs) and is consistent with measured PFAS profiles in consumed seafood. The second PC/cluster includes perfluorohexanesulfonic acid (PFHxS) and the PFOS precursor N-ethyl perfluorooctane sulfonamidoacetate (N-EtFOSAA), which are directly used or metabolized from fluorochemicals in consumer products such as carpet and food packaging. We find that the same compounds are associated with the same exposure sources in two North American populations, suggesting generalizability of results from the Faroese population.

CONCLUSIONS

We conclude that PFAS homologue profiles in serum provide valuable information on major exposure sources. It is essential to compare samples collected at similar time periods and to correct for demographic groups that are highly affected by differences in physiological processes (e.g., pregnancy). Information on PFAS homologue profiles is crucial for attributing adverse health effects to the proper mixtures or individual PFASs.