Tap Water Contributions to Plasma Concentrations of Poly- and Perfluoroalkyl Substances (PFAS) in a Nationwide Prospective Cohort of U.S. Women.

BACKGROUND: Between 2013 and 2015, concentrations of poly- and perfluoroalkyl substances (PFAS) in public drinking water supplies serving at least six million individuals exceeded the level set forth in the health advisory established by the U.S. Environmental Protection Agency. Other than data reported for contaminated sites, no systematic or prospective data exist on the relative source contribution (RSC) of drinking water to human PFAS exposures.
OBJECTIVES: This study estimates the RSC of tap water to overall PFAS exposure among members of the general U.S.
METHODS: We measured concentrations of 15 PFAS in home tap water samples collected in 1989-1990 from 225 participants in a nationwide prospective cohort of U.S. women: the Nurses' Health Study (NHS). We used a one-compartment toxicokinetic model to estimate plasma concentrations corresponding to tap water intake of PFAS. We compared modeled results with measured plasma PFAS concentrations among a subset of 110 NHS participants.
RESULTS: Tap water perfluorooctanoic acid (PFOA) and perfluorononanoic acid (PFNA) were statistically significant predictors of plasma concentrations among individuals who consumed [Formula: see text] cups of tap water per day. Modeled median contributions of tap water to measured plasma concentrations were: PFOA 12% (95% probability interval 11%-14%), PFNA 13% (8.7%-21%), linear perfluorooctanesulfonic acid (nPFOS) 2.2% (2.0%-2.5%), branched perfluorooctanesulfonic acid (brPFOS) 3.0% (2.5%-3.2%), and perfluorohexanesulfonic acid (PFHxS) 34% (29%-39%). In five locations, comparisons of PFASs in community tap water collected in the period 2013-2016 with samples from 1989-1990 indicated increases in quantifiable PFAS and extractable organic fluorine (a proxy for unquantified PFAS).
CONCLUSIONS: Our results for 1989-1990 compare well with the default RSC of 20% used in risk assessments for legacy PFAS by many agencies. Future evaluation of drinking water exposures should incorporate emerging PFAS. https://doi.org/10.1289/EHP4093.

Introduction

Poly- and perfluoroalkyl substances (PFAS) are a broad class of fluorinated aliphatic compounds that are widely used by industry and for commercial applications and have been detected in the serum of 98% of U.S. individuals (Khalil et al. 2016). Exposure to some PFAS has been associated with developmental, metabolic, and immune disorders in humans (Grandjean and Budtz-Jørgensen 2013; Liu et al. 2018; Vaughn et al. 2013). Exposure sources for PFAS are diverse and include consumer products, food, indoor dust, and drinking water (Domingo and Nadal 2017; Haug et al. 2011; Miralles-Marco and Harrad 2015). Between 2013 and 2015, public drinking water supplies serving at least six million people exceeded the lifetime health advisory level for PFAS established by the U.S. Environmental Protection Agency (U.S. EPA) (Hu et al. 2016). However, little information is available on the relative contribution of drinking water to observed levels in human plasma among the general population and how such exposures have changed over time.

Near contaminated sites, drinking water PFAS concentrations more than two orders of magnitude higher than the U.S. EPA health advisory level guideline ( range) have been reported (Emmett et al. 2006; Gyllenhammar et al. 2015; Landsteiner et al. 2014; Worley et al. 2017). At such locations, drinking water can account for up to 75% of total PFAS exposure (Emmett et al. 2006; Hoffman et al. 2010; Seals et al. 2011; Vestergren and Cousins 2009). Away from point sources, drinking water PFAS concentrations approximately 1,000-fold lower (low range) have been more commonly measured (Ericson et al. 2008; Quiñones and Snyder 2009). Even these lower concentrations have been associated with elevated serum concentrations in U.S. women (Hurley et al. 2016). Many other sources of PFAS (consumer products, dust, food) are known to be important for overall exposure of different human populations (Sunderland et al. 2018; Tokranov et al. 2018). Additional data are thus needed to better understand the contribution of drinking water to total PFAS exposures among the general U.S. population.

Interindividual variability in consumption rates and toxicokinetics can affect the contribution of drinking water to serum PFAS concentrations. Uptake and elimination of PFAS varies depending on age, race, menstrual status, childbirth, and breastfeeding (Gribble et al. 2015; Wong et al. 2014; Zhang et al. 2013). Regulatory agencies often use a term known as the relative source contribution (RSC) during risk assessments to represent the fraction of total PFAS exposure allocated to drinking water (Table S1). Most agencies, including the U.S. EPA, have adopted a default RSC value for drinking water of 20% for all PFAS (DeWitt 2015; Minnesota Department of Health 2017).

Here we estimate the RSC of drinking water collected in 1989–1990 to measured concentrations of PFAS in plasma from women in a nationwide prospective cohort [the Nurses’ Health Study (NHS)]. We focus on plasma samples because plasma is the available matrix in the NHS. We also examine temporal changes in PFAS concentrations in drinking water by comparing selected measurements from 1989–1990 to matched sampling locations between 2013 and 2016. We use these data to better understand the changing significance of drinking water as an exposure pathway for PFAS in the U.S. general population.

Methods

Study Population

In 1976, 121,700 female registered nurses between 30 and 55 years(y) of age were enrolled in the NHS (Belanger et al. 1978). Study participants have responded to mailed questionnaires on their medical, lifestyle, and health-related history every two years since that time. As part of the 1988 questionnaire cycle, study participants provided information on their tap water consumption. Their county of residence was obtained and updated from the mailing address provided during each questionnaire cycle.

Blood and Tap Water Sample Collection

Between May 1989 and September 1990, 32,826 NHS participants between the ages of 43 and 69 y provided matched blood and tap water samples. Blood samples were collected by the participants themselves using a blood collection kit containing three Sodium Heparin Collection tubes. Tap water samples were collected concurrently with blood samples from participants’ home kitchens using high-density polyethylene (HDPE) bottles with screw caps. Blood samples and frozen water samples were couriered overnight to the Brigham and Women’s Hospital in Boston, Massachusetts. Collection instructions specified that blood should be drawn between Monday and Thursday to ensure that the kit reached the lab on a weekday for immediate processing. Across all blood samples, 97% arrived within 26 h of being drawn.

Upon return to the lab, whole blood samples were immediately separated into plasma, white blood cells, and red blood cells, and stored in the vapor phase of liquid nitrogen freezers at . Tap water samples were thawed, acidified with nitric acid to , and stored at room temperature in a warehouse in Malden, Massachusetts. Storage at room temperature is acceptable for this study because PFAS are extremely stable in solution and do not fully degrade in water under these conditions (Wang et al. 2015).

We tested for potential absorption of PFAS to sampling bottles following previous methods established for soil extractions (Guelfo and Higgins 2013; Houtz et al. 2013). Briefly, we added ammonium hydroxide with methanol to an evacuated bottle that was then vortexed and sonicated for 0.5–1 h and placed on a shaker table for 2 h. We concentrated the aqueous sample with nitrogen prior to analysis for PFAS by high-performance liquid chromatography–tandem mass spectrometry (HPLC-MS/MS). Results confirmed the sample bottles did not contain a fluoropolymer lining (Figure S1).

We selected 225 study participants living in 22 states for home tap water PFAS measurements (Figure 1, Table 1). Sample selection was intended to cover geographically diverse areas and to overlap with plasma samples from 110 individuals analyzed for PFAS as part of a separate epidemiological investigation. Individuals included in this study had demographic, biometric, and lifestyle factors that were similar to the rest of cohort (Table S2). The study protocol was approved by the institutional review boards of Brigham and Women’s Hospital and the Harvard T.H. Chan School of Public Health. Return of self-administered questionnaires was considered informed consent. We also obtained informed consent for collecting blood and water samples.

Figure 1.

Locations of 225 home tap water samples obtained in 1989–1990.

Table 1

Overview of Nurses’ Health Study (NHS) participants ( women, 1989–1990 data) included in this study.

	Tap water intakea
	0–2 cups/day	3–5 cups/day	6–9 cups/day	10 or more
n	5	26	135	59
Age, y	53.0±6.8	55.4±7.5	54.3±6.8	51.8±6.2
White	5 (100%)	25 (96%)	129 (96%)	58 (98%)
BMI, kg m^−2	23.2±4.4	23.5±3.9	25.6±4.4	26.6±5.3
Weight, lb	128.8±18.6	141.6±28.0	152.3±25.9	157.4±35.5
Parity
No birth	1 (20%)	4 (15%)	2 (1%)	3 (5%)
1–3 birth	4 (80%)	13 (50%)	86 (64%)	43 (73%)
3+ births	0 (0%)	9 (35%)	47 (35%)	13 (22%)
Breastfeeding durationc
Never	2 (40%)	12 (46%)	56 (41%)	23 (39%)
<12  months	1 (20%)	7 (27%)	50 (37%)	18 (31%)
>=12  months	2 (40%)	7 (27%)	29 (21%)	18 (31%)
Menstruation status
Premenopause	1 (20%)	4 (15%)	21 (16%)	18 (31%)
Postmenopause	4 (80%)	22 (85%)	114 (84%)	41 (69%)
Seafood, servings/day	0.3±0.3	0.3±0.3	0.3±0.3	0.4±0.3
Popcorn, servings/day	0.3±0.4	0.1±0.2	0.1±0.2	0.2±0.2
Years residing at current location
<2	0 (0%)	3 (12%)	6 (4%)	5 (8%)
2–4	1 (20%)	2 (8%)	22 (16%)	9 (15%)
4–14	2 (40%)	11 (42%)	39 (29%)	22 (37%)
>14	2 (40%)	10 (38%)	68 (50%)	23 (39%)

Daily tap water consumption calculated as the sum of tap water consumed at all locations.

Racial category is dichotomized to white and nonwhite due to the small number of women belonging to other racial categories. White is defined without regard to Hispanic ethnicity.

Breastfeeding duration based on total months spent nursing all children reported in 1986 NHS questionnaire data.

Analysis of PFAS in Tap Water

Each archived tap water sample was homogenized by shaking vigorously before subsampling for analysis of 15 PFASs by LC-MS/MS at Harvard University. The perfluorosulfonic acids (PFSAs) measured included: perfluorobutane sulfonic acid (PFBS, four carbon chain length: C-4), linear and branched perfluorohexane sulfonic acid (nPFHxS, brPFHxS: C-6), linear and branched perfluorooctanesulfonic acid (nPFOS, brPFOS: C-8), and perfluorodecanesulfonic acid (PFDS: C-10). The perfluorocarboxylic acids (PFCAs) measured included: perfluoropentanoic acid (PFPeA: C-5), perfluoroheptanoic acid (PFHpA: C-7), linear and branched PFOA (nPFOA, brPFOA: C-8), linear and branched perfluorononanoic acid (nPFNA, brPFNA: C-9), perfluorodecanoic acid (PFDA: C-10), perfluoroundecanoic acid (PFUnDA: C-11), and perfluorododecanoic acid (PFDoA: C-12). The limit of detection (LOD) was calculated as the blank plus the average concentration at which the sample signal-to-noise ratio was three (Table S3). Perfluorohexanoic acid (PFHxA), N-methyl perfluorooctanesulfonamidoacetic acid (N-MeFOSAA) and N-ethyl perfluorooctanesulfonamidoacetic acid (N-EtFOSAA) were not detected in any samples and, therefore, are not included as analytes.

Each sample was adjusted to a pH between 5 and 7 so that all target analytes were in deprotonated state by adding ammonium hydroxide solution, and spiked with of a mass-labeled PFAS mixture (Wellington Laboratories) as internal standards for quantification. PFASs were extracted using an Oasis Wax solid phase extraction (SPE) cartridge (, sorbent, particle size) following methods established in prior work (Taniyasu et al. 2005; Zhang et al. 2016). Sample detection for native PFAS was performed using an Agilent 6460 triple quadrupole LC-MS/MS instrument (Agilent Technologies) equipped with an online SPE system described elsewhere (Zhang et al. 2016). One positive control (deionized water spiked with native PFAS) and one negative control (field blank or procedural blank) were included in every batch of 10 samples. The whole method recovery was between 80% and 98% for all analytes (Table S4). We excluded PFDS from our subsequent data analysis due to a low recovery rate of 44%. Variability between duplicates was less than 20% for all samples. Our method performance was comparable to a recently published ultrasensitive method for analyzing PFAS in small sample volumes (Dasu et al. 2017).

We quantified branched isomers for PFOS, PFOA, PFNA, and PFHxS using calibration standards for the linear isomers, assuming the same instrumental response factor, following previous work (Pellizzaro et al. 2018; Ullah et al. 2011; Zhang et al. 2016). PFDA and PFBS were detected in less than 20% of the 225 archived tap water samples and were excluded from subsequent analyses (Table 2).

Table 2

PFAS Concentrations in archived tap water samples collected from the Nurses’ Health Study (NHS) participants’ residences in 1989–1990 ().

	>LOD n (%)	Median (IQR)a [ng/L]	Max [ng/L]
PFCAs
PFPeA	84 (37)	0.61 (0.13, 1.72)	62.6
PFHpA	69 (31)	0.71 (0.35, 1.15)	9.65
nPFOA	134 (60)	0.57 (0.33, 1.36)	92.8
brPFOA	64 (28)	0.39 (0.17, 0.51)	12.0
nPFNA	62 (28)	0.23 (0.13, 0.29)	17.9
brPFNA	22 (10)	0.19 (0.10, 0.34)	16.1
PFDAb	43 (19)	0.22 (0.19, 0.94)	27.6
PFUnDA	109 (48)	1.15 (0.32, 2.42)	48.4
PFDoDA	76 (34)	2.31 (0.47, 6.40)	51.0
PFSAs
PFBSb	12 (5)	0.20 (0.19, 0.25)	2.97
nPFHxS	134 (60)	0.51 (0.23, 1.43)	28.0
brPFHxS	76 (34)	0.44 (0.18, 0.99)	4.93
nPFOS	107 (48)	0.69 (0.29, 1.80)	66.5
brPFOS	77 (34)	0.94 (0.51, 2.04)	18.7
PFDSc	27 (12)	0.11 (0.08, 0.52)	35.8

Median (IQR) were calculated for samples only.

PFDA and PFBS were excluded from subsequent analyses due to the proportion of nondetects.

PFDS was excluded from further analyses due to low recovery rate, even though its detection frequency is low as well.

For the remaining PFAS, we used the Robust Regression on Order Statistics (ROS) for data containing multiple detection limits to impute values for nondetects following established methods (Lee and Helsel 2005). First, a linear model was used to regress uncensored observed concentrations against their normal quantiles (order statistics). Censored concentrations were modeled using the parameters of the linear regression and normal quantiles of the censored observations, estimated by the exceedance probability of each censoring limit and the observation's rank. ROS assumes data are left-censored, and the relationship between concentrations and normal quantiles satisfies the assumptions of a linear regression. Water PFAS concentrations were log-transformed before ROS to correct for heteroskedasticity.

Analysis of PFAS in Plasma

Plasma samples for 110 individuals collected at the same time as tap water samples (1989–1990) were analyzed for PFAS by LC-MS/MS in Denmark as part of a separate epidemiological study. For this study, plasma concentrations of 12 PFAS were determined using online SPE followed by high-pressure LC-MS/MS at the University of Southern Denmark, following the methods reported elsewhere (Haug et al. 2009). The analytical system consisted of a Thermo Scientific EQuan MAX system connected to a TSQ Quantum Ultra Triple Stage Quadrupole mass spectrometer (Thermo Fisher Scientific). The online SPE was performed on a Betasil C8 () column, and the separation was performed on a Betasil C8 () column (Thermo Fisher Scientific). The LOD for all plasma PFAS was . Splits of blinded quality control samples were inserted in the different batches and used to calculate a between-batch coefficient of variation (CV) (Table S5). Measured PFAS concentrations were adjusted for batch effects following established methods (Rosner et al. 2008). A linear model was first fit to regress PFAS concentrations on batch indicator dummy variables. PFAS concentrations were then recalibrated by subtracting the difference between the coefficient of each individual batch and the average of the coefficients of all batches. All analyses reported here were based on the recalibrated values of plasma PFAS.

Comparability of Cross-Lab Water and Plasma Measurements

Two labs were used to conduct PFAS analysis in this study (Harvard University for water samples and the University of Southern Denmark for plasma samples) using parallel analytical methods. Each lab specializes in the respective matrix analyzed to ensure the highest-quality data were produced for this study. Rigorous QA protocols were followed for both water and plasma measurements reported here, as summarized in Table S4.

Statistical Methods

For the five major PFAS (PFOA, PFNA, nPFOS, brPFOS, PFHxS), we used a generalized additive model (GAM) to regress log-transformed plasma PFAS against log-transformed tap water PFAS with a cubic spline smoothing function (Hastie 2016). The relative importance of tap water to PFAS exposure is expected to vary according to tap water intake. We therefore stratified the GAM analysis by tap water consumption rates of greater or less than 8 cups per day. We calculated Cook’s distance statistic to identify data points with high leverage. The adjusted GAM included age, race/ethnicity, body weight, menstruation status, parity, breastfeeding history, years residing at current address, seafood consumption, and popcorn consumption. Our analysis was based on responses to the 1990 questionnaire, which were modeled as continuous variables or categorized as shown in Table 1. We calculated the duration of residence (years) at the tap water sample location based on residential addresses in 1976, 1986, 1988, and 1990. Statistical analyses were performed using R (version 3.3.2; R Development Core Team), and statistical significance was defined as . Concentrations for tap water and plasma PFAS were log-transformed (ln) before fitting a general additive model.

To illustrate the implications of high method reporting limits (MRL) in data from the U.S. EPA’s third Unregulated Contaminants Monitoring Rule (UCMR3) for human exposure, we reported effect sizes based on a hypothetical increase in tap water PFAS concentrations from the median of 1989–1990 samples to the MRL in U.S. EPA UCMR3. PFOS isomers are not reported separately in UCMR3, and here we assume 70% of PFOS is linear and 30% of PFOS is branched, based on literature values (Yu et al. 2015).

Toxicokinetic (TK) Model Simulations of Plasma PFAS Concentrations

One-compartment toxicokinetic (TK) models have been successfully used to relate external PFAS exposures from sources such as drinking water to serum PFAS concentrations (Gomis et al. 2017; Lorber and Egeghy 2011; Thompson et al. 2010). TK models can be used to assess the relative importance of different exposure sources (Lorber and Egeghy 2011; Trudel et al. 2008; Vestergren and Cousins 2009) and to establish benchmark doses for PFAS from toxicological and epidemiological studies (Goeden et al. 2019). For example, Grandjean and Budtz-Jørgensen (2013) used such modeling to suggest a health advisory level for perfluorooctanoic acid (PFOA) in drinking water of based on vaccine antibody responses in children and an uncertainty factor of 10 (Grandjean and Budtz-Jørgensen 2013).

We used the difference between measured plasma PFAS concentrations for 110 individuals and modeled contributions to plasma based on measured levels in tap water in 1989–1990 to estimate the RSC of drinking water to overall PFAS exposure. For individuals whose tap water PFAS concentrations were below detection, we assumed tap water was a negligible fraction of their overall exposure and excluded them from TK modeling. For example, consumption of of water per day for a woman at our analytical detection limit for PFOA (average ) corresponds to a plasma PFOA value of , which is 1.7% of the measured mean in the NHS plasma samples. TK modeling was conducted only for plasma PFAS with a , namely PFOA, PFNA, nPFOS, brPFOS, and PFHxS (Table S5).

For TK modeling, we assumed steady state from chronic drinking water exposure and simulated the resulting plasma concentrations as follows: where (L/day) is the tap water consumption rate reported in NHS questionnaire data; (day) is the half-life in the human body for a given PFAS; (kg) is individual body weight reported in the NHS survey data; (mL/kg) is the volume of distribution, which is the theoretical volume that reflects how PFAS distribute throughout the body; and (ng/L) is the measured drinking water PFAS concentration in this study.

We used Monte Carlo (MC) simulations to investigate how variability in input parameters, such as self-reported water consumption rates, chemical half-lives, and the volume of distribution, affected the modeled relative source contribution of tap water. We ran 300 iterations of the TK model using input parameters that were randomly sampled from specified distributions (Table S6). Means and variances for the elimination half-life in the human body and volumes of distribution were based on a literature review (search terms are provided in Supplemental Information Section S1) and modeled using a truncated () lognormal distribution (Table S7). No information is available on the volume of distribution of PFAS in humans (Thompson et al. 2010); therefore, our parameterization was based on two animal studies (Ohmori et al. 2003; Sundström et al. 2012). A precedent for the use of animal data when human TK parameters are not available has been established in prior work (Gomis et al. 2016; Zhang et al. 2013). Data on half-lives of PFAS in human plasma are available from multiple studies, and we thus used a weighted average based on the inverse of the variance reported in each study in our TK model.

NHS participants reported tap water consumption inside and outside the home in the following categories as part of the 1988 questionnaire: none (zero), 1 to 2, 3 to 5, 6 to 9, and 10 or more cups per day. We used a uniform distribution with the lower and upper limit of each category, except for the “none” (assigned zero) and “10 or more” (assigned 11) categories. We calculated the mean and 95% probability interval among the 300 iterations for the median, mean, 25th and 75th percentile of the RSC of drinking water. We quantified the contribution of different input parameters to the total variability in the RSC of drinking water as the square of the correlation coefficient normalized to the sum of the squared correlation coefficients (Wang et al. 2016).

Temporal Changes in Tap Water Exposures to PFAS

To better understand temporal shifts in drinking water exposures to PFAS, we matched the locations of each of the 225 NHS participant’s homes to the U.S. EPA’s UCMR3 database by county. For the 2013–2015 data at each location, we calculated both the maximum PFAS concentration and detection frequency. This calculation was necessary because the UCMR3 method reporting limits (range ) were much higher than the detection limits in this study (Table S3). We calculated the Spearman’s rank order correlation coefficient between tap water PFAS concentrations in 1989–1990 and those in the period 2013–2015. In our TK modeling, we replaced drinking water PFAS concentrations measured in 1989–1990 with those from the UCMR database (2013–2016) to simulate potential changes in plasma concentrations. We used two approaches to account for samples below detection in the UCMR database. First, we replaced nondetects with the MRL divided by square root of 2. Second, we replaced nondetects with a concentration of zero. We then compared the sensitivity of modeled results to the treatment of nondetects. We did not use regression-based imputation methods such as ROS because the detection frequency in UCMR database was only 4%.

Recent work indicates many precursors to PFAS are present in environmental samples that are not detectable by targeted analyses (Barzen-Hanson et al. 2017; Wang et al. 2017). We thus conducted a pilot investigation that measured total extractable organofluorine (EOF) as a proxy for the total burden of fluorinated compounds in tap water from five cities where we had participant samples. In 2016, we collected samples from the same municipal water supplies as the original participants’ homes (Table S8). At each location, the tap water was flushed for 2 to 3 min before duplicate samples were collected in prerinsed HDPE bottles. One field blank was prepared at each sampling location using HPLC-grade water. Two duplicates were taken at each of the five sampling locations. For targeted PFAS in the 2016 water samples, we followed the same analytical procedures as reported above for the 1989–1990 tap water samples, except that no pH adjustment was needed and was used in the analysis. Field blanks for the five locations sampled in the summer of 2016 were all below the LOD for the targeted 15 PFAS. The relative percentage difference identified in duplicate samples for 15 PFAS was . Average concentrations measured at each location are reported here.

EOF was measured on paired 1989–1990 and 2016 tap water samples at Örebro University in Örebro, Sweden, using extraction methods that separate inorganic and organic fluorine and combustion ion chromatography, following previously established methods (Miyake et al. 2007). Briefly, SPE and dilution procedures were optimized to adequately separate inorganic fluorine (such as fluoride used for water fluoridation) from the organic fraction. Then, all the organic fluorine was converted into hydrogen fluoride (HF) during combustion and absorbed into a hydrogen peroxide solution. The concentration of was quantified by ion chromatography. The method detection limit (MDL) for EOF () was calculated from the average concentration of three procedure blanks plus three times the SD. Samples reported for EOF have been blank corrected and additional QA/QC information is provided in Table S4. Organofluorine concentrations contributed by quantifiable PFAS are calculated as: where is the number of fluorine atoms in a single PFAS, is the molecular weight of fluorine (g/mol), is the molecular weight of the corresponding PFAS (g/mol), is the concentration of PFAS (ng/L).

To help interpret drivers of temporal differences in tap water PFAS between 1989–1990 and 2013–2016, we obtained data on the number of industries known to release PFAS between 1987 and 2015 (U.S. EPA 2018). No information on the magnitudes of PFAS releases is available from the EPA Toxic Release Inventory (TRI) database. Relevant industrial sources were thus identified following the methods outlined in previous work using the North American Industrial Classification System (NAICS) code (Table S9) (Zhang et al. 2016).

Results

Concentrations of PFAS in Matched Tap Water and Plasma Samples

Concentrations of the 15 PFAS measured in 1989–1990 tap water samples from the NHS study are shown in Table 2. The most frequently detected PFAS were nPFOA (60%), nPFHxS (60%), nPFOS (48%), PFUnDA (48%), and PFPeA (37%). We detected 11 out of the 15 PFAS targeted for analysis in 28% of the samples. Median concentrations of PFAS were in the low ng/L range. Highest concentrations measured in tap water collected in 1989–1990 were for nPFOA, for nPFOS, and for nPFHxS.

All five PFAS included in toxicokinetic modeling (PFOA, PFNA, nPFOS, brPFOS, and PFHxS) were detected in 100% of plasma samples from NHS participants (Table 3). Median concentrations measured in NHS participants were: PFOA (), PFOS (sum of linear and branched isomers, ), PFHxS (), and PFNA (). Median concentrations of PFAS in plasma for NHS participants reported here were similar to those for adult women who participated in the 1999–2000 NHANES survey (PFOA: , PFOS: , PFHxS: , PFNA: ) (Calafat et al. 2007).

Table 3

PFAS Concentrations in tap water and plasma collected in 1989–1990 from Nurses’ Health Study (NHS) participants ().

	Plasma PFAS			Tap water PFAS
	>LOD n (%)	Median (IQR) [ng/mL]	Max [ng/mL]	>LOD n (%)	Median (IQR)a [ng/L]	Max [ng/L]
PFOA	110 (100)	4.78 (3.56, 6.47)	77.72	50 (45)	0.57 (0.46, 1.65)	104.74
PFNA	110 (100)	0.61 (0.41, 0.89)	11.53	31 (28)	0.13 (0.12, 0.35)	1.98
nPFOS	110 (100)	15.86 (12.27, 22.00)	79.83	58 (53)	0.43 (0.25, 1.43)	16.56
brPFOS	110 (100)	11.99 (9.10, 16.13)	47.45	50 (45)	1.12 (0.44, 1.91)	10.38
PFHxS	110 (100)	1.89 (1.04, 2.85)	52.66	67 (61)	0.57 (0.10, 1.86)	12.19

Medians (IQR) were only calculated for samples greater than the LOD (limit of detection).

Measured tap water PFOA and PFNA concentrations were statistically significant predictors of plasma concentrations among the 66 NHS participants who consumed cups of tap water per day (Figure 2). We estimated the expected increases in plasma PFOA and PFNA associated with an increase in tap water concentrations from the 1989–1990 median to the U.S. EPA UCMR3 method reporting limit (Table 4). These associations remained statistically significant after adjustment for age, race/ethnicity, body weight, menstruation status, parity, breastfeeding history, years residing in current address, seafood consumption, and popcorn consumption. After adjustments for covariates, point estimates changed by for PFOA and for PFNA. We included and excluded one individual in the PFOA model because Cook’s distance statistic suggested this was a high leverage point. Excluding this observation did not change the statistical significance of the unadjusted association. The effect size was reduced from 250% (95% confidence interval (CI): 130%, 450%) to 80% (95% CI: 0.20%, 220%).

Figure 2.

Associations between tap water and plasma PFAS concentrations among NHS participants in 1989–1990, estimated with a Generalized Additive Model (GAM) with a cubic spline smoothing function. Model estimates were not adjusted for covariates. Shaded area showed 95% confidence interval (CI). Participants who consume cups of tap water per day () are shown in light gray and participants who consume cups of tap water per day () are shown in blue. The PFOA model was run twice after removing one individual with high leverage. For a description of the models, see Table 4.

Table 4

Modeled difference in plasma PFAS concentrations (%) with a hypothetical increase in tap water PFAS concentrations.

	MRL	Exposure changeb	Unadjusted				Adjustede
			≥8  cups/day		<8  cups/day		≥8  cups/day		<8  cups/day
	ng/L	ng/L	% Change (95% CI)	p-value	% Change (95% CI)	p-value	% Change (95% CI)	p-value	% Change (95% CI)	p-value
PFOA	20.0	19.4	253 (127, 448)	<0.001	73.9 (−3.9, 215)	0.3	227 (23.3, 766)	<0.001	23.5 (−52.3, 220)	0.4
PFOAc	20.0	19.4	79.8 (0.2, 223)	0.04	70.2 (−4.2, 202)	0.3	55.2 (−46, 346)	0.2	24.4 (−51.5, 219)	0.3
PFNA	20.0	19.9	2720 (444, 14522)	<0.01	618 (91.4, 2592)	0.03	1540 (98.1, 13483)	0.06	275.2 (−30.2, 1917)	0.2
nPFOSd	28.0	27.6	13.3 (−29.3, 81.6)	0.8	105 (−11.9, 375)	0.4	−0.4 (−62.6, 165)	1.0	75 (−55.4, 587)	0.4
brPFOS	12.0	10.9	21.1 (−14.0, 70.6)	0.3	8.7 (−34.1, 79.2)	0.6	12.8 (−57, 194)	0.6	−22.9 (−70.9, 104)	0.3
PFHxS	30.0	29.4	157 (−41.3, 1025)	0.5	54.6 (−17.5, 190)	0.1	134 (−69, 1660)	0.5	28 (−71.2, 470)	0.4

Exposure change estimated based on a hypothetical increase in tap water PFAS concentrations from the median of 1989–1990 samples to the method reporting limits (MRL) in U.S. EPA UCMR3. Concentrations for tap water and plasma PFAS were log-transformed before fitting a general additive model.

The difference between the MRL and median tap water PFAS concentrations in 1989–1990.

PFOA model with one high leverage data point removed. This individual had a tap water PFOA concentration of and a plasma PFOA concentration of . She maintained the same residential location between 1986 and 1990 but moved at least once in the preceding decade (1976 to 1986).

PFOS isomers are not reported separately in UCMR3. Here we assume based on literature values that 70% of PFOS is linear and 30% of PFOS is branched (Yu et al, 2015).

Adjusted for age, race/ethnicity, body weight, menstruation status, parity, breastfeeding history, years residing at current address, seafood consumption, and popcorn consumption.

TK Modeling Results

One-compartment TK modeling suggested that, in 1989–1990, the median relative contribution of tap water consumption to plasma concentrations was 12% [interquartile range (IQR): 7.7%– 20%] for PFOA, 13% (IQR: 6.4%–21%) for PFNA, 2.2% (IQR: 0.9%–.8%) for nPFOS, 3.0% (IQR: 1.2%–6.5%) for brPFOS, and 34% (IQR: 15%–61%) for PFHxS (Figure 3, Table S10). There was no statistically significant difference between the RSC of tap water to plasma PFAS concentrations across different groups of residential years at the same address (Table S11).

Figure 3.

Estimated relative source contribution (RSC) of tap water to overall PFAS exposure using a one-compartment toxicokinetic model. The dashed line represents the default RSC (20%) used in risk assessment to derive drinking water advisory levels. The upper panel shows the RSC among 110 Nurses’ Health Study (NHS) participants in 1989–1990. The lower-panel shows the RSC stratified by number of years at the same residential location. Box and whisker plots show the fifth, 25th, 50th, 75th and 95th percentiles among each group. Individual estimates are denoted by gray dots with small random variation added to their horizontal position for better separation. All data used to generate this figure are provided in Tables S8 and S9.

Median RSCs based on MC simulations were similar () to deterministic (single central point) estimates. However, the probabilistically estimated mean RSCs for PFNA and PFHxS increased by more than 5% (Table S12). The modeled RSCs for PFOA and PFOS showed less variability than did PFNA and PFHxS (Figure S2). For PFNA, uncertainty in the chemical half-life accounted for 60% of the variability in the modeled RSC. For PFHxS, the volume of distribution and drinking water intake were the two most important determinants of overall variability, contributing to 47% and 43%, respectively (Figure S3).

Temporal Changes in PFAS Concentrations in Tap Water

The U.S. EPA’s UCMR3 database includes data on six legacy PFAS measured between 2013 and 2015 in U.S. drinking water supplies (Hu et al. 2016). Four of the UCMR3 compounds (PFOA, PFNA, PFOS, and PFHxS) overlap with PFAS detected in this work. Comparing the 1989–1990 NHS tap water sampling locations to UCMR3 data revealed that 46% (representing 53 out of 144 counties) overlapped with sampling locations that had detectable concentrations of at least one PFAS in the UCMR3 database (U.S. EPA 2015). PFOS concentrations measured in 1989–1990 NHS samples were significantly and positively correlated with both the maximum concentration (Spearman correlation , ), and the detection frequency (, ) in the 2013–2015 UCMR3 samples. PFNA concentrations measured in 1989–1990 NHS samples were significantly and positively correlated with the maximum concentration in the 2013–2015 UCMR3 samples (, ). We did not find significant correlations for the maximum concentration for PFOA (, ), and PFHxS (, ), or the detection frequency for PFOA (, ), PFNA (, ), and PFHxS (, ).

Figure 4 shows modeled median plasma concentrations are sensitive to the treatment of nondetects in the UCMR3 database due to its relatively high MRL of in comparison with our work. We replaced nondetects for PFOA in the UCMR3 database by the MRL divided by the square root of 2. This replacement results in a median modeled plasma concentration in 2013–2015 that is more than 10-fold higher () than modeled contributions of drinking water in 1989–1990 (). Alternately, replacing nondetects by zero results in median modeled plasma concentrations in 2013–2015 that are also zero.

Figure 4.

Estimated plasma PFAS concentrations for 225 individuals using a one-compartment toxicokinetic model and paired recent (2013–2015) measurements of PFASs in tap water with those from 1989–1990. Data under “2013–2015(a)” show nondetects replaced by the method ; Data under “2013–2015(b)” show nondetects replaced by zero. Box and whisker plots indicate fifth, 25th, 50th, 75th and 95th percentiles among each group. Individual estimates are shown as gray dots with small random variation added to their horizontal position for better separation.

We measured a 5- to 320-fold increase in total EOF among the tap water samples collected from the same city in 1989–1990 and 2016 (Table 5). The fraction of unquantified organic fluorine (unknown PFAS) was above 60% (range: 60%–94%) in all the tap water samples collected in 2016, in comparison with 8%– 89% in 1989–1990 samples. It is possible that some degradation of polyfluoroalkyl precursors to PFAS occurred over time in the archived tap water samples from 1989–1990, but this degradation would not affect the total EOF values reported.

Table 5

Extractable organic fluorine levels from tap water samples matched by city at five locations in Massachusetts (MA1 – MA5) in 1989–1990 and 2016

Extractable organic fluorine (ng/L)	MA1		MA2		MA3		MA4		MA5
	1989–1990	2016	1989–1990	2016	1989–1990	2016	1989–1990	2016	1989–1990	2016
PFOA	0.2	6.2	0.5	1.7	0.9	4.8	0.6	0.9	1.3	0.9
PFOS	0.4	1.6	0.4	0.8	1.2	4.2	0.5	0.3	0.6	0.3
Other PFCAs	0.1	7.4	0.8	4.2	1.3	9.6	0.6	1.7	0.0	5.1
Other PFSAs	0.3	4.3	0.3	1.7	1.5	5.6	0.2	0.7	0.4	0.1
PFOS precursors	0.0	0.0	0.3	0.0	0.0	0.0	0.0	0.0	0.0	0.0
Unknown	6.7	135.6	19.8	105.2	2.9	39.4	0.2	58.5	5.4	9.6

1989–1990 tap water samples were collected from five participants’ home addresses, one sample at each location.

2016 tap water samples were collected in the cities from the same municipal water supplies as original participant’s homes. Two samples were collected at each location.

Discussion

In U.S. women without known occupational exposure to PFAS, we found significant positive associations between tap water and plasma PFOA and PFNA concentrations in samples collected in 1989–1990. These relationships were strongest among individuals who reported consuming cups of tap water per day. A default RSC of 20% is used by U.S. federal and state agencies in risk assessments for deriving drinking water advisory levels for legacy PFAS (Post et al. 2017; U.S. EPA 2016). Our results suggest that this RSC was reasonable for legacy PFAS in 1989–1990. However, recent increases in legacy PFAS concentrations have been reported in many of the same municipal water supplies in the U.S. EPA UCMR database. Our analysis shows large increases in unknown EOF (PFAS not quantified in the targeted analysis) in five paired tap water samples from 1989–1990 and 2016.

Tap water PFAS concentrations measured in 1989–1990 in this study were mostly in the low ng/L range and below current drinking water guidelines. Only 14 out of 225 samples exceeded the state of Michigan’s guideline for PFOS () and six samples exceeded the state of New Jersey’s advisory level for PFOA () and PFNA (). Concentrations reported in this study for 1989–1990 are much lower than those more recently measured near contaminated sites, and thus these archived samples provide unique insight into past drinking water exposures among the U.S. general population (Emmett et al. 2006; Hoffman et al. 2010; Hölzer et al. 2008).

Concentrations and detection of PFOS in drinking water samples collected in 2013–2015 as part of the U.S. EPA UCMR3 database were positively correlated with tap water concentrations measured in 1989–1990. No statistically significant temporal correlations were observed for the other three PFAS examined (PFHxS, PFNA, PFOA). This observation may reflect the large quantities of PFOS released to the environment between 1951 and 2002 in comparison with the other compounds. For example, global production of perfluorooctanesulfonyl fluoride (POSF), the parent chemical to PFOS, was 5 to 37 times greater than that of PFCAs over that time period (Earnshaw et al. 2014; Paul et al. 2009; Wang et al. 2014).

The number of potential PFAS point sources across the United States has fluctuated between 1990 and 2015 (Figure S4). Source abundance peaked in both 1990 and the early 2000s but declined in the later parts of both decades. Over the same time period, there has been a phaseout in production and use of legacy PFAS such as PFOS and PFOA in North America and Europe and rapid growth in the number and diversity of replacement compounds (Stockholm Convention 2009; U.S. EPA 2006; Wang et al. 2014; Wang et al. 2017). It is currently unclear how these changes have affected exposures to PFAS from drinking water.

Data collected in this work provide suggestive evidence of an increase in the estimated plasma PFAS concentrations of legacy PFAS (PFOS, PFOA, PFNA, and PFHxS) through tap water consumption between 1989–1990 and 2013–2015 (Figure 4). This finding is sensitive to the treatment of nondetects in the UCMR3 database. We also report an increase in the burden of unquantified organofluorine compounds based on pilot data of five pairs of tap water samples (Table 5). These pilot data provide preliminary evidence of a potential increase in accumulation of unknown organofluorine compounds in drinking water that are not included in routine environmental monitoring and targeted analysis. Such patterns are consistent with increasing numbers of PFAS produced and released to the environment in recent years (Ritscher et al. 2018; Strynar et al. 2015; Wang et al. 2013). This increase contrasts with serum trends in PFOS and PFOA in many human populations in the United States and Europe that have been widely observed to decline since the early 2000s, reflecting the success of phaseouts in chemical production and stewardship programs (Dassuncao et al. 2018).

Tap water concentrations were statistically significant predictors of plasma PFOA and PFNA concentrations for NHS participants, but not of PFOS and PFHxS. The high median RSC of tap water for PFHxS (34.1%) among individuals with detectable tap water concentrations likely reflects its long half-life (7.3 y) in human plasma (Olsen et al. 2007).

TK modeling suggests median concentrations of the five major PFAS measured in the NHS drinking water samples from 1989–1990 would result in exposures equivalent to between 0.01 and in plasma, assuming consumption of of water per day for a woman. This finding is relevant for the U.S. general population because plasma PFAS concentrations among NHS participants in 1989–1990 and adult women in the 1999–2000 NHANES survey are almost identical (Calafat et al. 2007). The median RSC of drinking water for PFOA exposures in 1989–1990 was 12%, with an IQR between 7.7% to 20%. This finding is similar to values reported in prior work (Post et al. 2012; Post et al. 2017), although to the best of our knowledge our study represents the earliest time period for general population exposures.

Our results are consistent with existing studies on the relative importance of different exposure pathways (diet, drinking water, dust, consumer product, inhalation). Dietary ingestion was thought to be the major exposure pathway for PFOS among the general adult populations in North America and Europe between 2011 and 2016. Estimated dietary PFOS exposures ranged from 66% to 99%, whereas drinking water contributed 0.10% to 22% of total exposure (Egeghy and Lorber 2011; Gebbink et al. 2015; Shan et al. 2016). Drinking water has been estimated to account for between 0.70% and 37% of overall exposure to PFOA (Tian et al. 2016; Vestergren and Cousins 2009). Our study is the first to estimate the relative importance of drinking water for PFHxS and PFNA.

The RSC of drinking water to total PFAS exposure is important for designing guidelines that are health protective. If the RSC of drinking water is overestimated, total exposures of the general population may exceed the reference dose, even when drinking water PFAS concentrations fall below the guideline. For PFAS other than PFOS and PFOA, the RSCs of tap water have not been resolved. Our analyses suggest tap water can be an equally important contributor for PFNA (median RSC is 13%), and an even more important contributor for PFHxS (median RSC is 34%). Pilot data on EOF suggest that unknown PFAS in drinking water may have increased and that static RSC values may be problematic.

This study has several strengths and limitations. It is nested in a well-established cohort study that includes comprehensive information on diet, lifestyle, and medical history collected using validated questionnaires. A detailed residential history is available for participants and the cohort employs stringent laboratory quality control procedures, allowing the development of a novel combination of both multivariate statistical and toxicokinetic models. Data presented in this study are limited in terms of sample size, which restricts their statistical power to inform the potential temporal changes of PFAS in tap water over the past 27 years. We measured PFAS in matched tap water and plasma samples at a single point in time. However, PFAS have long half-lives in human body, and plasma PFAS concentrations are not prone to short-term, within-person variability. Finally, information needed to quantify the magnitudes of nonwater sources of PFAS exposure, such as consumer product use, dietary consumption preferences for PFAS-containing items, and indoor exposure, is not available in the NHS. Such data are needed to evaluate the relative importance of other exposure pathways for PFAS.

Conclusions

The default RSC value (20%) for tap water currently used in PFAS risk assessments of drinking water are consistent with findings from this study, based on archived tap water and plasma samples from 1989–1990. We estimated that tap water contributed between 2.2% and 34% of plasma concentrations for the five PFAS examined in 1989–1990 among a subsample of participants in the NHS. Temporal and geographic shifts in PFAS production and use may have altered the RSC from drinking water since that time.

The ongoing presence of PFOS in drinking water samples collected between 2013–2016 reinforces the concern that accumulation of legacy PFAS in the environment can elevate human exposure long after the cessation of global production. The rapid decline in human exposure to legacy PFAS since 2000 has been predominantly driven by the phaseout of PFOS and its precursors (Dassuncao et al. 2018). If exposures from consumer product use decrease due to regulatory actions and voluntary phaseouts, environmental sources may drive future exposures (Gomis et al. 2017; Vestergren and Cousins 2009). We conclude that a more holistic approach that targets the entire class of fluorinated compounds would enable more health-protective drinking water guidelines (Blum et al. 2015). Our pilot data showing temporal changes in EOF suggest that consideration of drinking water exposures to emerging PFAS is warranted for the general U.S. population (Sunderland et al. 2018).

Click here for additional data file.

Click here for additional data file.