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

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

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

PBPK modeling to support risk assessment of pyrethroid exposure in French pregnant women

BACKGROUND

Pyrethroids are widely used pesticides and are suspected to affect children's neurodevelopment. The characterization of pyrethroid exposure during critical windows of development, such as fetal development and prenatal life, is essential to ensure a better understanding of pyrethroids potential effects within the concept of Developmental Origins of Health and Disease.

OBJECTIVE

The aim of this study was to estimate maternal exposure of French pregnant women from biomonitoring data and simulate maternal and fetal internal concentrations of 3 pyrethroids (permethrin, cypermethrin and deltamethrin) using a multi-substance pregnancy-PBPK (physiologically based pharmacokinetics) model. The estimated maternal exposures were compared to newly proposed toxicological reference values (TRV) children specific also called draft child-specific reference value to assess pyrethroid exposure risk during pregnancy i.e. during the in utero exposure period.

METHODS

A pregnancy-PBPK model was developed based on an existing adult pyrethroids model. The maternal exposure to each parent compound of pregnant women of the Elfe (French Longitudinal Study since Childhood) cohort was estimated by reverse dosimetry based on urinary biomonitoring data. To identify permethrin and cypermethrin contribution to their common urinary biomarkers of exposure, an exposure ratio based on biomarkers in hair was tested. Finally, exposure estimates were compared to current and draft child-specific reference values derived from rodent prenatal and postnatal exposure studies.

RESULTS

The main contributor to maternal pyrethroid diet intake is cis-permethrin. In blood, total internal concentrations main contributor is deltamethrin. In brain, the major contributors to internal pyrethroid exposure are deltamethrin for fetuses and cis-permethrin for mothers. Risk is identified only for permethrin when referring to the draft child-specific reference value. 2.5% of the population exceeded permethrin draft child-specific reference value.

CONCLUSIONS

A new reverse dosimetry approach using PBPK model combined with human biomonitoring data in urine and hair was proposed to estimate Elfe pregnant population exposure to a pyrethroids mixture with common metabolites.

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.

PBPK Modeling Approach to Predict the Behavior of Drugs Cleared by Metabolism in Pregnant Subjects and Fetuses

This study aimed to develop a physiologically based pharmacokinetic (PBPK) model that simulates metabolically cleared compounds' pharmacokinetics (PK) in pregnant subjects and fetuses. This model accounts for the differences in tissue sizes, blood flow rates, enzyme expression levels, plasma protein binding, and other physiological factors affecting the drugs' PK in both the pregnant woman and the fetus. The PBPKPlus™ module in GastroPlus® was used to model the PK of metoprolol, midazolam, and metronidazole for both non-pregnant and pregnant groups. For each of the three compounds, the model was first developed and validated against PK data in healthy non-pregnant volunteers and then applied to predict the PK in the pregnant groups. The model accurately described the PK in both the non-pregnant and pregnant groups and explained well the differences in the plasma concentration due to pregnancy. When available, the fetal plasma concentration, placenta, and fetal tissue concentrations were also predicted reasonably well at different stages of pregnancy. The work described the use of a PBPK approach for drug development and demonstrates the ability to predict differences in PK in pregnant subjects and fetal exposure for metabolically cleared compounds.

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.

Development of a physiologically based toxicokinetic model for lead in pregnant women: The role of bone tissue in the maternal and fetal internal exposure

Epidemiological studies have shown associations between prenatal exposure to lead (Pb) and neurodevelopmental effects in young children. Prenatal exposure is generally characterized by measuring the concentration in the umbilical cord at delivery or in the maternal blood during pregnancy. To assess internal Pb exposure during prenatal life, we developed a pregnancy physiologically based pharmacokinetic (p-PBPK) model that to simulates Pb levels in blood and target tissues in the fetus, especially during critical periods for brain development. An existing Pb PBPK model was adapted to pregnant women and fetuses. Using data from literature, both the additional maternal bone remodeling, that causes Pb release into the blood, and the Pb placental transfers were estimated by Bayesian inference. Additional maternal bone remodeling was estimated to start at 21.6 weeks. Placental transfers were estimated between 4.6 and 283 L.day-1 at delivery with high interindividual variability. Once calibrated, the p-PBPK model was used to simulate fetal exposure to Pb. Internal fetal exposure greatly varies over the pregnancy with two peaks of Pb levels in blood and brain at the end of the 1st and 3rd trimesters. Sensitivity analysis shows that the fetal blood lead levels are affected by the maternal burden of bone Pb via maternal bone remodeling and by fetal bone formation at different pregnancy stages. Coupling the p-PBPK model with an effect model such as an adverse outcome pathway could help to predict the effects on children's neurodevelopment.

Biomonitoring equivalents for perfluorooctanoic acid (PFOA) for the interpretation of biomonitoring data

BACKGROUND

Perfluorooctanoic acid (PFOA) is detected in the blood of virtually all biomonitoring study participants. Assessing health risks associated with blood PFOA levels is challenging because exposure guidance values (EGVs) are typically expressed in terms of external dose. Biomonitoring equivalents (BEs) consistent with EGVs could facilitate health-based interpretations.

OBJECTIVE

To i) derive BEs for serum/plasma PFOA corresponding to non-cancer EGVs of the U.S. Environmental Protection Agency (U.S. EPA), the Agency for Toxic Substances and Disease Registry (ATSDR) and Health Canada, and ii) compare with PFOA concentrations from national biomonitoring surveys.

METHODS

Starting from EGV points of departure, we employed pharmacokinetic data/models and uncertainty factors. Points of departure in pregnant rodents (U.S. EPA 2016, ATSDR) were converted into fetus and pup serum concentrations using an animal gestation/lactation pharmacokinetic model, and equivalent human fetus and child concentrations were converted into BEs in maternal serum using a human gestation/lactation model. The point of departure in adult rodents (Health Canada) was converted into a BE using experimental data. For epidemiology-based EGVs (U.S. EPA 2023, draft), BEs were directly based on epidemiological data or derived using a human gestation/lactation pharmacokinetic model. BEs were compared with Canadian/U.S. biomonitoring data.

RESULTS

Non-cancer BEs (ng/mL) were 684 (Health Canada, 2018) or ranged from 15 to 29 (U.S. EPA, 2016), 6-10 (ATSDR, 2021) and 0.2-0.8 (U.S. EPA, 2023, draft). Ninety-fifth percentiles of serum levels from the 2018-2019 Canadian Health Measures Survey (CHMS) and the 2017-2018 National Health and Nutrition Examination Survey (NHANES) were slightly below the BE for ATSDR, and geometric means were above the non-cancer BEs for the U.S. EPA (2023, draft).

CONCLUSION

Non-cancer BEs spanned three orders of magnitude. The lowest BEs were for EGVs based on developmental endpoints in epidemiological studies. Concentrations in Canadian/U.S. national surveys were higher than or close to BEs for the most recent non-cancer EGVs.

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.

AOP-helpFinder 2.0: Integration of an event-event searches module

To support the use of alternative methods in regulatory assessment of chemical risks, the concept of adverse outcome pathway (AOP) constitutes an important toxicological tool. AOP represents a structured representation of existing knowledge, linking molecular initiating event (MIE) initiated by a prototypical stressor that leads to a cascade of biological key event (KE) to an adverse outcome (AO). Biological information to develop such AOP is very dispersed in various data sources. To increase the chance of capturing relevant existing data to develop a new AOP, the AOP-helpFinder tool was recently implemented to assist researchers to design new AOP. Here, an updated version of AOP-helpFinder proposes novel functionalities. The main one being the implementation of an automatic screening of the abstracts from the PubMed database to identify and extract event-event associations. In addition, a new scoring system was created to classify the identified co-occurred terms (stressor-event or event-event (which represent key event relationships) to help prioritization and support the weight of evidence approach, allowing a global assessment of the strength and reliability of the AOP. Moreover, to facilitate interpretation of the results, visualization options are also proposed. The AOP-helpFinder source code are fully accessible via GitHub, and searches can be performed via a web interface at http://aop-helpfinder-v2.u-paris-sciences.fr/.

Community-facing toxicokineticmodels to estimate PFAS serum levels based on life history and drinking water exposures

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) are drinking water contaminants. Tools to assess the potential body burden associated with drinking PFAS-contaminated water may be helpful for public health assessment of exposed communities.

METHODS

We implemented a suite of one-compartment toxicokinetic models using extensively calibrated toxicokinetic parameters (half-life and volume of distribution). We implemented the models both in the R programming language for research purposes, and as a web estimator for the general public (built in typescript.js). These models simulate exposure to PFAS water concentrations for individuals with varying characteristics such as age, sex, weight, and breastfeeding history. The models account for variability and uncertainty in parameter inputs to produce Monte Carlo-based estimates of serum concentration. For children, the models additionally account for gestational exposure, lactational exposure, and potential exposure through formula feeding. For adults who have borne children, the models account for clearance through birth and breastfeeding. We ran simulations of individuals with known PFAS water and serum concentrations to evaluate the model. We then compared the predicted serum PFAS concentrations to measured data.

RESULTS

The models accurately estimate individual-level serum levels for each PFAS for most adults within ½ order of magnitude. We found that the models somewhat overestimated serum concentrations for children in the tested locations, and that these overestimates are generally within an order of magnitude.

DISCUSSION

This paper presents scientifically robust models that allow users to estimate serum PFAS concentrations based on known PFAS water concentrations and physiologic information. However, accuracy in historical water concentration inputs, exposure from non-drinking water sources, and life-history characteristics of individuals present a complex problem for individual estimation. Additional refinements to the model suite to improve the prediction of individual results may consist of including duration of exposure and additional life-history characteristics.

Determinants of maternal and neonatal PFAS concentrations: a review

Per- and polyfluoroalkyl substances (PFAS) are used for their properties such as stain and water resistance. The substances have been associated with adverse health outcomes in both pregnant mothers and infants, including pre-eclampsia and low birthweight. A growing body of research suggests that PFAS are transferred from mother to fetus through the placenta, leading to in utero exposure. A systematic review was performed using the PubMed database to search for studies evaluating determinants of PFAS concentrations in blood matrices of pregnant mothers and neonates shortly after birth. Studies were included in this review if an observational study design was utilized, exposure to at least one PFAS analyte was measured, PFAS were measured in maternal or neonatal matrices, at least one determinant of PFAS concentrations was assessed, and results such as beta estimates were provided. We identified 35 studies for inclusion in the review and evaluated the PFAS and determinant relationships among the factors collected in these studies. Parity, breastfeeding history, maternal race and country of origin, and household income had the strongest and most consistent evidence to support their roles as determinants of certain PFAS concentrations in pregnant mothers. Reported study findings on smoking status, alcohol consumption, and pre-pregnancy body mass index (BMI) suggest that these factors are not important determinants of PFAS concentrations in pregnant mothers or neonates. Further study into informative factors such as consumer product use, detailed dietary information, and consumed water sources as potential determinants of maternal or neonatal PFAS concentrations is needed. Research on determinants of maternal or neonatal PFAS concentrations is critical to estimate past PFAS exposure, build improved exposure models, and further our understanding on dose-response relationships, which can influence epidemiological studies and risk assessment evaluations. Given the potential for adverse outcomes in pregnant mothers and neonates exposed to PFAS, it is important to identify and understand determinants of maternal and neonatal PFAS concentrations to better implement public health interventions in these populations.

Systematic review of physiologically based kinetic lactation models for transfer of xenobiotic compounds to milk

Lactational elimination has been described mathematically for nearly 50 years. Over 40 published articles, containing >50 physiologically based kinetic (PBK) lactation models were included in the systematic review. These PBK models described the lactational elimination of xenobiotic compounds in humans, rats, mice, and dairy cows and goats. A total of 78 compounds have been modelled, ranging from industrial chemicals, pesticides, to pain medication, antibiotics, and caffeine. Few models included several species or compounds, and models were thus generally not translational or generic. Three dairy cow models mechanistically described the intramammary disposition of pharmaceuticals after intramammary administration, including volume changes caused by milking, while empirically describing the remaining pharmacokinetics. The remaining models were semi- or whole body PBK models, describing long-term exposure of environmental pollutants, or short-term exposure of pharmaceuticals. The absolute majority described the disposition to the mammary gland or milk with perfusion limited compartments, but permeability limited models were available as well. With long-term exposure, models often included changes in milk volume and/or consumption by the offspring, and changes in body weight of offspring. Periodic emptying of the mammary gland, as with feeding or milking, was sparsely applied. Rodent models used similar physiological parameters, while values of physiological parameters applied in human models could range widely. When milk composition was included in the models, it most often included the fat content. The review gives an extensive overview of the applied functions and modelling strategies of PBK lactation models.

Effect of perfluoroalkyl exposure in pregnancy and infancy on intrauterine and childhood growth and anthropometry. Sub study from COPSAC2010 birth cohort

BACKGROUND

Perfluoroalkyl substances PFOS and PFOA are persistent and bioaccumulative exogenous chemicals in the human body with a range of suspected negative health effects. It is hypothesised that exposure during prenatal and early postnatal life might have particularly detrimental effects on intrauterine and childhood growth. In a Danish longitudinal mother-child cohort we investigate effect of PFOS and PFOA in pregnancy and infancy on intrauterine and childhood growth and anthropometry.

METHODS

COPSAC₂₀₁₀ is an ongoing population based mother-child cohort of 738 pregnant women and their children followed from 24 week gestation with longitudinal deep clinical phenotyping until age 10 years. In this observational cohort sub study plasma PFOS and PFOA concentrations were semi-quantified by untargeted metabolomics in the mothers at week 24 and 1 week postpartum and in the children at ages 6 and 18 months and calibrated using a targeted pipeline. We examined associations to intrauterine and childhood growth and anthropometry, including interactions with child sex. Untargeted and targeted blood metabolomics profiles were integrated to investigate underlying mechanisms.

FINDINGS

Pregnancy plasma PFOA concentrations were associated with lower birth size -0.19 [-0.33; -0.05] BMI z-score per 1-ng/mL and increased childhood height (z-scored) at age 6: 0.18 [0.05; 0.31], but there was no association between childs' own infancy plasma PFOA concentration and height. Pregnancy plasma PFOS concentrations were also associated with lower birth BMI (-0.04 [-0.08; -0.01]), but in childhood pregnancy plasma PFOS concentration interacted with child sex on BMI and fat percentage at 6 years with negative associations in girls and positive in boys. The effect of maternal plasma PFOS concentration on lower girl BMI was borderline mediated through increasing child plasma lactosyl-ceramide levels (p-mediation=0.08). Similarly the effect of maternal plasma PFOS concentration on higher boy fat percentage was borderline mediated through increasing child plasma lactosyl-ceramide levels (p-mediation=0.07). Infancy concentrations of plasma PFOS associated with lower height in childhood, -0.06 z-score at age 6 [-0.19; -0.03].

INTERPRETATION

Higher PFOS and PFOA plasma concentrations during pregnancy had detrimental effects on fetal growth. The effects on childhood growth were not similar as PFOA increased child height, opposite of PFOS in multipollutant models suggesting a differing fetal programming effect. Sex specific growth effects were borderline mediated through an altered lactosyl-ceramide metabolism, proposing a possible mechanism of PFOS that has long-lasting health consequences in this observational study.

FUNDING

All funding received by COPSAC are listed on www.copsac.com. The Lundbeck Foundation (Grant no R16-A1694); The Novo Nordic Foundation (Grant nos NNF20OC0061029, NNF170C0025014, NNF180C0031764) The Ministry of Health (Grant no 903516); Danish Council for Strategic Research (Grant no 0603-00280B) and The Capital Region Research Foundation have provided core support to the COPSAC research center. Effort from JALS is supported by R01HL123915, R01HL141826, and R01HL155742 from NIH/NHLBI. CEW was supported by the Swedish Heart Lung Foundation (HLF 20180290, HLF 20200693). BC has received funding for this project from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 946228). The funding agencies did not have any role in design and conduct of the study; collection, management, and interpretation of the data; or preparation, review, or approval of the manuscript.

Risk Assessment of Perfluorooctane Sulfonate (PFOS) using Dynamic Age Dependent Physiologically based Pharmacokinetic Model (PBPK) across Human Lifetime

The widespread use of Perfluorooctane sulfonate (PFOS) in everyday life, its long half-life, and the lipophilicity that makes it easily accumulate in the body, raises the question of its safe exposure among different population groups. There are currently enough epidemiological studies showing evidence of PFOS exposure and its associated adverse effects on humans. Moreover, it is already known that physiological changes along with age e.g. organ volume, renal blood flow, cardiac output and albumin concentrations affect chemicals body burden. Human biomonitoring cohort studies have reported PFOS concentrations in blood and autopsy tissue data with PFOS present in sensitive organs across all human lifespan. However, to interpret such biomonitoring data in the context of chemical risk assessment, it is necessary to have a mechanistic framework that explains show the physiological changes across age affects the concentration of chemical inside different tissues of the human body. PBPK model is widely and successfully used in the field of risk assessment. The objective of this manuscript is to develop a dynamic age-dependent PBPK model as an extension of the previously published adult PFOS model and utilize this model to predict and compare the PFOS tissue distribution and plasma concentration across different age groups. Different cohort study data were used for exposure dose reconstruction and evaluation of time-dependent concentration in sensitive organs. Predicted plasma concentration followed trends observed in biomonitoring data and model predictions showed the increased disposition of PFOS in the geriatric population. PFOS model is sensitive to parameters governing renal resorption and elimination across all ages, which is related to PFOS half-life in humans. This model provides an effective framework for improving the quantitative risk assessment of PFOS throughout the human lifetime, particularly in susceptible age groups. The dynamic age-dependent PBPK model provides a step forward for developing such kind of dynamic model for other perfluoroalkyl substances.

Accumulation of PFAS by livestock - determination of transfer factors from water to serum for cattle and sheep in Australia

In this study accumulation has been determined of several per- and polyfluoroalkyl substances (PFAS) from ingested water to steady state serum concentration for adult beef cattle and sheep raised on a hobby farm impacted by PFAS contamination. PFAS concentrations in stock water were stable for more than a year, they were non-measurable in grass but present at very low levels in soil which equated to just 1% of the intake from water. Prior to quantifying PFAS in cattle serum there had been no breeding for 18 months. Although there were high concentrations of several PFAS in the water, only perfluorooctane sulphonate (PFOS) and perfluorohexane sulphonate (PFHxS) were in cattle serum in appreciable amounts; perfluoroheptane sulphonate (PFHpS), perfluorononanoic acid (PFNA) and perfluorodecanoic acid (PFDA) were much lower. Transfer factors (TFs) for cattle were calculated by dividing steady state serum concentration (ng PFAS/mL) by water concentration (µg PFAS/L). Average and upper estimate TF values for cattle were calculated; the former were 140 (total PFOS, i.e. tPFOS), 130 (PFHpS), 65 (PFHxS), 170 (PFNA), and 120 (PFDA). Previous investigation campaigns at the farm provided relative steady state serum PFAS concentrations for sheep and cattle that allowed adjustment of the cattle TFs. The resulting average estimate TFs for ewes (non-pregnant and not lactating) were 20 (tPFOS) and 30 (PFHxS), other PFAS were not measurable in sheep serum. Discussion on using these TFs in human health risk assessments is provided. With certain assumptions/caveats the TFs allow estimations of PFAS steady state serum concentrations for use in preliminary human health risk assessments (HHRAs) when only PFAS in stock water is known.

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.

Multistate models of developmental toxicity: Application to valproic acid-induced malformations in the zebrafish embryo

For the determination of acute toxicity of chemicals in zebrafish (Danio rerio) embryos, the OECD test guideline 236, relative to the Fish Embryo Toxicity Test (FET), stipulates a dose-response analysis of four lethal core endpoints and a quantitative characterization of abnormalities including their time-dependency. Routinely, the data are analyzed at the different observation times separately. However, observations at a given time strongly depend on the previous effects and should be analyzed jointly with them. To solve this problem, we developed multistate models for occurrence of developmental malformations and live events in zebrafish embryos exposed to eight concentrations of valproic acid (VPA) the first five days of life. Observations were recorded daily per embryo. We statistically infer on model structure and parameters using a numerical Bayesian framework. Hatching probability rate changed with time and we compared five forms of its time-dependence; a constant rate, a piecewise constant rate with a fixed hatching time at 48 h post fertilization, a piecewise constant rate with a variable hatching time, as well as a Hill and Gaussian form. A piecewise constant function of time adequately described the hatching data. The other transition rates were conditioned on the embryo body concentration of VPA, obtained using a physiologically-based pharmacokinetic model. VPA impacted mostly the malformation probability rate in hatched and non-hatched embryos. Malformation reversion probability rates were lowered by VPA. Direct mortality was low at the concentrations tested, but increased linearly with internal concentration. The model makes full use of data and gives a finer grain analysis of the teratogenic effects of VPA in zebrafish than the OECD-prescribed approach. We discuss the use of the model for obtaining toxicological reference values suitable for inter-species extrapolation. A general result is that complex multistate models can be efficiently evaluated numerically.

Maternal-Fetal Pharmacology of Drugs: A Review of Current Status of the Application of Physiologically Based Pharmacokinetic Models

Pregnancy and the postpartum period are associated with several physiological changes that can alter the pharmacokinetics (PK) and pharmacodynamics (PD) of drugs. For certain drugs, dosing changes may be required during pregnancy and postpartum to achieve drug exposures comparable to what is observed in non-pregnant subjects. There is very limited data on fetal exposure of drugs during pregnancy, and neonatal exposure through transfer of drugs via human milk during breastfeeding. Very few systematic clinical pharmacology studies have been conducted in pregnant and postpartum women due to ethical issues, concern for the fetus safety as well as potential legal ramifications. Over the past several years, there has been an increase in the application of modeling and simulation approaches such as population PK (PopPK) and physiologically based PK (PBPK) modeling to provide guidance on drug dosing in those special patient populations. Population PK models rely on measured PK data, whereas physiologically based PK models incorporate physiological, preclinical, and clinical data into the model to predict drug exposure during pregnancy. These modeling strategies offer a promising approach to identify the drugs with PK changes during pregnancy to guide dose optimization in pregnancy, when there is lack of clinical data. PBPK modeling is also utilized to predict the fetal exposure of drugs and drug transfer via human milk following maternal exposure. This review focuses on the current status of the application of PBPK modeling to predict maternal and fetal exposure of drugs and thereby guide drug therapy during pregnancy.

Assessing the impacts on fetal dosimetry of the modelling of the placental transfers of xenobiotics in a pregnancy physiologically based pharmacokinetic model

The developmental origin of health and diseases theory supports the critical role of the fetal exposure to children's health. We developed a physiologically based pharmacokinetic model for human pregnancy (pPBPK) to simulate the maternal and fetal dosimetry throughout pregnancy. Four models of the placental exchanges of chemicals were assessed on ten chemicals for which maternal and fetal data were available. These models were calibrated using non-animal methods: in vitro (InV) or ex vivo (ExV) data, a semi-empirical relationship (SE), or the limitation by the placental perfusion (PL). They did not impact the maternal pharmacokinetics but provided different profiles in the fetus. The PL and InV models performed well even if the PL model overpredicted the fetal exposure for some substances. The SE and ExV models showed the lowest global performance and the SE model a tendency to underprediction. The comparison of the profiles showed that the PL model predicted an increase in the fetal exposure with the pregnancy age, whereas the ExV model predicted a decrease. For the SE and InV models, a small decrease was predicted during the second trimester. All models but the ExV one, presented the highest fetal exposure at the end of the third trimester. Global sensitivity analyses highlighted the predominant influence of the placental transfers on the fetal exposure, as well as the metabolic clearance and the fraction unbound. Finally, the four transfer models could be considered depending on the framework of the use of the pPBPK model and the availability of data or resources to inform their parametrization.

Well-tempered MCMC simulations for population pharmacokinetic models

A full Bayesian statistical treatment of complex pharmacokinetic or pharmacodynamic models, in particular in a population context, gives access to powerful inference, including on model structure. Markov Chain Monte Carlo (MCMC) samplers are typically used to estimate the joint posterior parameter distribution of interest. Among MCMC samplers, the simulated tempering algorithm (TMCMC) has a number of advantages: it can sample from sharp multi-modal posteriors; it provides insight into identifiability issues useful for model simplification; it can be used to compute accurate Bayes factors for model choice; the simulated Markov chains mix quickly and have assured convergence in certain conditions. The main challenge when implementing this approach is to find an adequate scale of auxiliary inverse temperatures (perks) and associated scaling constants. We solved that problem by adaptive stochastic optimization and describe our implementation of TMCMC sampling in the GNU MCSim software. Once a grid of perks is obtained, it is easy to perform posterior-tempered MCMC sampling or likelihood-tempered MCMC (thermodynamic integration, which bridges the joint prior and the posterior parameter distributions, with assured convergence of a single sampling chain). We compare TMCMC to other samplers and demonstrate its efficient sampling of multi-modal posteriors and calculation of Bayes factors in two stylized case-studies and two realistic population pharmacokinetic inference problems, one of them involving a large PBPK model.

Modeled prenatal exposure to per- and polyfluoroalkyl substances in association with child autism spectrum disorder: A case-control study

BACKGROUND/OBJECTIVE

Per- and polyfluoroalkyl substances (PFAS) display neurobehavioral toxicity in laboratory animal studies. We examined associations of modeled prenatal maternal exposure to PFAS with child diagnosis of autism spectrum disorder (ASD).

METHODS

Participants were 453 mother-child pairs from CHARGE (CHildhood Autism Risk from Genetics and Environment), a population-based case-control study. Children underwent psychometric testing and were clinically confirmed for ASD (n = 239) or typical development (TD, n = 214). At the end of the clinic visit, maternal blood specimens were collected. We quantified nine PFAS in maternal serum samples collected when their child was 2-5 years old. As surrogate in utero exposure, we used a model built from external prospective data in pregnancy and 24 months post-partum and then reconstructed maternal PFAS serum concentrations during pregnancy in this case-control sample. We used logistic regression to evaluate associations of modeled prenatal maternal PFAS concentrations with child ASD.

RESULTS

Modeled prenatal maternal perfluorohexane sulfonate (PFHxS) and perfluorooctane sulfonate (PFOS) were borderline associated with increased odds of child diagnosis of ASD (per nanogram per milliliter increase: odds ratio [OR] = 1.46; 95% confidence interval [CI]: 0.98, 2.18 for PFHxS, OR = 1.03; 95% CI: 0.99, 1.08 for PFOS). When compared to the lowest quartile (reference category), the highest quartile of modeled prenatal maternal PFHxS was associated with increased odds of child diagnosis of ASD (OR = 1.95; 95% CI: 1.02, 3.72).

CONCLUSIONS

In analyses where modeled prenatal maternal PFAS serum concentrations served as in utero exposure, we observed that prenatal PFHxS and PFOS exposure, but not other PFAS, were borderline associated with increased odds of child diagnosis of ASD. Further studies in which PFAS concentrations are prospectively measured in mothers and children at a range of developmental stages are needed to confirm these findings.

Integrative Strategy of Testing Systems for Identification of Endocrine Disruptors Inducing Metabolic Disorders-An Introduction to the OBERON Project

Exposure to chemical substances that can produce endocrine disrupting effects represents one of the most critical public health threats nowadays. In line with the regulatory framework implemented within the European Union (EU) to reduce the levels of endocrine disruptors (EDs) for consumers, new and effective methods for ED testing are needed. The OBERON project will build an integrated testing strategy (ITS) to detect ED-related metabolic disorders by developing, improving and validating a battery of test systems. It will be based on the concept of an integrated approach for testing and assessment (IATA). OBERON will combine (1) experimental methods (in vitro, e.g., using 2D and 3D human-derived cells and tissues, and in vivo, i.e., using zebrafish at different stages), (2) high throughput omics technologies, (3) epidemiology and human biomonitoring studies and (4) advanced computational models (in silico and systems biology) on functional endpoints related to metabolism. Such interdisciplinary framework will help in deciphering EDs based on a mechanistic understanding of toxicity by providing and making available more effective alternative test methods relevant for human health that are in line with regulatory needs. Data generated in OBERON will also allow the development of novel adverse outcome pathways (AOPs). The assays will be pre-validated in order to select the test systems that will show acceptable performance in terms of relevance for the second step of the validation process, i.e., the inter-laboratory validation as ring tests. Therefore, the aim of the OBERON project is to support the organization for economic co-operation and development (OECD) conceptual framework for testing and assessment of single and/or mixture of EDs by developing specific assays not covered by the current tests, and to propose an IATA for ED-related metabolic disorder detection, which will be submitted to the Joint Research Center (JRC) and OECD community.