Occup. Hyg., pp. 1–13 Ann. Occup. Hyg., Vol. 56, Ann. No. 9, pp. 1025–1037, 2012 ©The TheAuthor Author2012. 2012.Published Publishedby by Oxford Oxford University University Press Press on onbehalf behalfofofthe theBritish British Occupational Occupational Hygiene Hygiene Society Society doi:10.1093/annhyg/mes023 doi:10.1093/annhyg/mes023 Retrospective Exposure Assessment of Perfluorooctanoic Acid Serum Concentrations at a Fluoropolymer Manufacturing Plant SUSAN R. WOSKIE1*, REBECCA GORE1 and KYLE STEENLAND2 1 Department of Work Environment, University of Massachusetts Lowell, One University Avenue, Lowell, MA 01854, USA; 2Department of Environmental Health, Emory University, Atlanta, GA 30322, USA Received 21 October 2011; in final form 27 February 2012 ; published online 26 April 2012 Perfluorooctanoic acid (PFOA) is a suspect human carcinogen, causes neonatal loss, liver enlargement, and a variety of tumors in rodents, and has been associated with increased cholesterol levels in humans. Mortality analyses of worker cohorts have not been conclusive or consistent. As part of a series of epidemiologic studies of workers in a West Virginia plant that manufactures fluoropolymers, estimates of serum PFOA for the worker cohort were developed for the period of 1950–2004. An existing database of 2125 worker biomarker measurements of serum PFOA was used to model retrospective exposures. Historical PFOA serum levels for eight job category/job group combinations were modeled using linear mixed models to account for repeated measures, along with exposure determinants such as cumulative years worked in potentially exposed jobs, the amount of C8 used or emitted by the plant over time, as well as a four-knot restricted cubic spline function to reflect the influence of process changes over calendar time on exposure. The modeled biomarker levels matched well with measured levels, including those collected independently as part of a community study of PFOA levels (Spearman correlations of 0.8 for internal data comparisons and 0.6 for external data comparisons). These annualized PFOA serum estimates will be used in a series of morbidity and mortality studies of this worker cohort. Keywords: exposure assessment; exposure determinants; exposure science; historical exposure assessment; retrospective exposure assessment exposure metric; Teflon INTRODUCTION Ammonium perfluorooctanoate (APFO), the ammonium salt of perfluorooctanoic acid (PFOA), is used as a surfactant in the polymerization of tetrafluoroethylene to make certain fluoropolymers such as polytetrafluoroethylene (PTFE). This perfluorooctanoate is also known as C8. PTFE products marketed as Teflon and Gore-tex among others have a variety of uses including a coatings for non-stick pans, waterrepellent clothing, dome materials, wire and cable, inert tubing and laboratory supplies, and semiconductor applications among many others. *Author to whom correspondence should be addressed. Tel: þ978-934-3295; fax: þ978-452-5711; e-mail: [email protected] PFOA is a suspect human carcinogen and causes neonatal loss in mice (US EPA, 2005). PFOA has been shown to cause tumors of the pancreas, liver, and testes in rodents. It also causes liver enlargement in rodents and non-human primates (Steenland et al., 2010). There are two worker cohorts exposed to PFOA that have been studied for mortality and morbidity outcomes: one at a 3M plant in Minnesota and one at a Dupont plant in West Virginia. Studies of the Dupont cohort have reported associations with higher cholesterol (Sakr et al., 2007a,b). However, mortality analyses have not been conclusive or consistent across the two cohorts (Steenland et al., 2010). Elevated standardized mortality ratios for the Dupont cohort for diabetes and a kidney cancer were found when the referent group was workers at other Dupont 1025 1 of 13 2 of 13 1026 S. Woskie, R. Gore and K. Steenland plants (Leonard et al., 2008). A positive exposure– response trend for heart disease was reported when cumulative exposure cutpoints were chosen from all exposed workers at Dupont and a 10-year lag was used (Sakr et al., 2009). Although not significant, suggestive positive trends in internal exposure–response analyses for diabetes, stroke, prostate cancer, and pancreatic cancer were reported in the 3M cohort by Lundin et al. (2009). Non-differential exposure misclassification can impact the validity and reduce the precision of a cohort study. In studies where exposures are dichotomous, it can bias estimates toward the null. In studies using exposure–response analysis with exposure measured on a continuous scale, the impact of mis-measured exposure is more complex but often results in biased exposure–response trends, making interpretation of exposure–disease associations more difficult (Armstrong, 1990; Dosemeci et al., 1990; Steenland et al., 2000). In theory, improvements in exposure assessment that eliminate misclassification can have dramatic effects on relative risk estimates; however, it is difficult to evaluate the impact in realworld studies (Blair et al., 2007). Previous PFOA worker cohort exposure assessments The historical exposure reconstruction described here attempts to improve the exposure assessments used to date for the occupational studies of PFOA workers. The PFOA exposure assessment for the mortality study of the 3M plant was confined to categorizing workers into not exposed, probably exposed, and definitely exposed (Lundin et al., 2009). A more extensive exposure assessment was conducted for the Dupont plant (Kreckmann et al., 2009). Using a dataset from a cross-sectional survey of serum levels for 1000 workers in 2004, 60 job groups were assigned to one of three exposure categories, primarily based on their median 2004 job serum level. Then cohort members were given a cumulative serum PFOA level calculated as the sum of the number of years (1950–2004) in each of the three job categories times the job category mean exposure level for all of the jobs in the category in the 2004 survey (Kreckmann et al., 2009). The exposure assessment reported here differs in that it does not use the 2004 serum measurements to assign jobs to similar exposure groups (which we define via job category and job groups within a job category, see below); rather, work process and potential contact with PFOA were used to assign jobs to similar exposure groups (job category/job group combinations). In addition, this job exposure matrix increases the number of similar exposure groups from three to eight, and, perhaps most important, unlike the previous assessment, it uses observed serum levels by job from 1979 to 2004 to model predicted serum levels over time. METHODS Worker blood samples The first measurements of worker exposures to PFOA were in 1972. From 1972 to 1981, the blood analysis for PFOA was done on whole blood by Dupont using the Wickbold torch method. This method measured total organic fluorine using an oxyhydrogen torch/bomb decomposition method followed by a spectrophotometric method of detection. During that time, all organic fluorine was assumed to be related to PFOA. The limit of detection (LOD) for this method was not available in the Dupont documentation; however, for this method, the lowest concentration reported in our modeling dataset was 0.015 ppm and no non-detectable values were reported. In 1981, the blood analysis was converted to a gas chromatography with electron capture detector (GC-ECD) method. Whole-blood samples underwent lyophilization, and then the dried residue was combined with methanolic HCl followed by extraction of the methyl esters into hexane. An internal standard was used to correct for analytical variations. The GC-ECD LOD was reported as 0.01 ppm by Dupont. However, out of 1248 analyzed by this method, in our modeling dataset, only one was reported as nondetectable. Based on replicate analysis of 114 samples by two different methods, correction factors (average 1.28) were developed by Dupont and applied to the collected blood data to convert samples analyzed by Wickbold torch to GC-ECD PFOA analysis equivalents. Beginning in 2003, worker samples were analyzed by an high pressure liquid chromatography with tandem mass spectrometry (LC-MS/MS) method, which uses serum, not whole blood. The LOD for the LC/MS/MS method was reported as 0.005 ppm and none of the samples in the modeling dataset were non-detectable. PFOA binds to serum proteins (Han et al., 2003; Jones et al., 2003; Wu et al., 2009). A study comparing serum, plasma, and whole-blood samples of 12 humans found that the average ratio of PFOA in serum to ethylenediaminetetraacetic acid (EDTA) treated whole blood was 2.1 (1.7–2.8) (Ehresman et al., 2007). Another study with five human samples found a mean PFOA ratio of 1.37 plasma/whole blood (range 1.19–1.82) (Karrman et al., 2006). When both datasets are combined, the Retrospective workplace PFOA serum concentrations 3 of 13 1027 average ratio is 1.8 (1.19–2.8). All whole-blood data collected before 2003 were adjusted to the serum equivalent by multiplying by 1.8. Results are reported in parts per million (1 ppm 1 lg ml1 PFOA). For the model development, the biomarker dataset of serum PFOA concentrations was confined to those samples where the subject had been in the job group for at least 1 year before the sample was collected (2125 samples from 1308 workers). Similar exposure groups (job category/job group combinations) Through a series of interviews and group discussions with current and retired Dupont staff familiar with plant processes, work history entries (division, department, and job title) were categorized into similarly exposed job groups (30) and larger job categories (5). Within the five larger job categories, the job groups were evaluated by stem and leaf plots and multiple comparison least squares means (LS means) tests, from a mixed model of the log serum concentration that included calendar year, job group, and total years in potentially C8-exposed jobs prior to a specific serum sample. Based on sample size and the LS means tests, only three of the five larger job categories retained subdivisions of job groups (for a total of eight job category/job group combinations). Fine powder/granular PTFE job category (direct exposure to PFOA). Fluropolymers have been produced since 1951 at the Dupont chemical plant in Washington, West Virginia. Within this area, the chemical operator job group runs the reactors and were in charge of measuring the powdered C8 (APFO) for introduction into the reactors. Fine powder/ dispersion used C8 from the start and granular began use in 1965. Both used powdered C8 up until 1989 when introduction of a premixed liquid C8 was begun. The finishing operator job group works in the powder drying and packaging area and transfer liquid polymer dispersions into shipping totes. Finishing operators also do dryer cleanout and help with maintenance of the dryers. In the analysis, this job category included a dichotomous variable for working in the chemical operator job group, resulting in two job category/job group combinations. Fluorinated ethylene propylene and perfluoroalkoxy fluoropolymer job category (direct exposure to PFOA). The first co-polymer [fluorinated ethylene propylene (FEP)] began pilot operations in 1953 but only used C8 in powder form from 1975 to 1991, and then they converted to the premixed liquid form of C8. The second co-polymer [perfluoroalkoxy (PFA)] was co-manufactured with FEP and used powdered C8 from 1973 to 1982. Then production was transferred to a non-aqueous, non-C8 process. In 1998, PFA manufacture was again switched back to an aqueous C8 process. Beginning in 1995, the FEP/PFA chemical operators and finishing operators rotated jobs, eliminating any finer distinction in potential exposures based on job tasks. Non-PFOA (C8) use in Teflon polymer/co-polymer production job category (intermittent direct or plant background PFOA exposure). This job category included a number of jobs that were associated with and housed in buildings adjacent to the Teflon polymer/co-polymer production departments. These included a combined job group of Teflon polymer/ co-polymer laboratorians, engineers, upper level supervisiors, and clerks as well as production jobs in Teflon polymer/co-polymer production operations that never used C8, such as the ethylenetetrafluoroethylene fluoropolymer (Tefzel) and fluorotelomer (telomer) co-polymer operations. The tetrafluorethylene (TFE) monomer operation, which never used C8, began in 1950, was treated as a separate job group due to their work location adjacent to the PTFE operations. In the analysis, this job category included a dichotomous variable for working in the monomer operator job group, resulting in two job category/job group combinations. Maintenance job category (intermittent direct or plant background PFOA exposure). Until 1983, most maintenance workers were assigned to a centralized ‘mechanical’ or ‘utility’ division, even though much of their work may have been in the Teflon polymer/co-polymer area. From 1983 on, workers could be assigned either to a centralized division job group or to a Teflon polymer/co-polymer division job group. In the analysis, this job category included a dichotomous variable for having been assigned to Teflon/co-polymer maintenance job group, resulting in two job category/job group combinations. Non-Teflon/co-polymer production division with no PFOA use job category (plant background PFOA exposure). This group included the other production divisions in the plant that did not use C8 (acrylics, Butacite, Delrin, engineering polymers, compounding, nylon filaments, specialty compounding, as well as power services, utility pool, and non-Teflon polymer/co-polymer-associated administrative jobs in engineering, business services, and other plant services). Jobs considered potentially exposed to C8 included the job categories of fine powder/granular, FEP/PFA, non-PFOA (C8) use jobs in Teflon and co-polymer production, and the Teflon/co-polymer maintenance job group. 4 of 13 1028 S. Woskie, R. Gore and K. Steenland Exposure modeling using biological monitoring data The goal of the modeling reported here was to develop the best fitting predictive models using serum data collected on a subset of individual workers over time, in order to predict serum levels for each individual worker in the cohort over time. Therefore, we developed five separate regression models, one for each of the job categories representing similar exposure conditions, to best capture the timerelated process changes in each job category. Since serum concentrations were highly skewed and more lognormal than normal based on evaluation of probability plots and histograms, the natural log of the serum concentrations were used in the modeling. The natural log of the serum concentrations were examined using a variety of methods of modeling time trends; either mixed models with a continuous measure for time or broken stick regression (two or three piece linear) modeling of time, lowess smoothing, or restricted cubic splines. Continuous measures were unsatisfactory because serum samples were clustered in time due to a sampling campaign approach. Piece-wise linear models were unsatisfactory since process change in manufacturing is more likely to result in a smoother and more gradual continuous change in serum levels, since the half-life of PFOA in serum is estimated to be 3.5 years (Olsen et al., 2007). Smoothing was judged to be a more realistic approach to modeling time; however, lowess models do not produce model predictors that can be used with other datasets to estimate an outcome based on independently available X values. Since the goal of our modeling was to estimate serum levels for individuals for each year of their work history, lowess models did not work. Instead, we used a restricted cubic spline (RCS) model, which produces model predictions by forcing the first and second derivatives of the function to agree at the knots. By using the RCS model, the function is also constrained to be linear at the tails, avoiding some of the problems with using unconstrained cubic spline models. These RCS models were developed using the RCSPLINE SAS macro developed by Frank Harrell who also reported that the location of the knots in an RCS model is not as important as the number of knots (Harrell et al., 1988; Harrell, 2001). For the job categories with direct C8 exposure, process changes over calendar time were modeled with a four-knot RCS function that utilized a linear term (T1) and two additional terms defined by the knots and cubic powers (T2 and T3). The number of knots (three to five) was determined by the model fit (Akaike information criterion using restricted maximum likelihood) as well as by comparing the predicted to the measured values across 5-year time windows. Although we tried a series of knot locations for each number of knots, in the end, the quantiles for knot location suggested by Harrell (2001) for each number of knots provided the best fit. For four knots, this was the 5, 35, 65, and 95 percentile of the distribution of number of years after 1970 that the serum samples were collected. For fine powder/granular, the knots were at 11, 13, 17, and 34 years after 1970. For FEP/ PFA, the knots were at 9, 13, 25, and 34 years after 1970. For each model, covariates that most improved model fit and therefore predictive power were chosen, at the cost of some cross-model consistency. These covariates included the cumulative number of prior years in potentially C8-exposed job, job group, the annual amount of PFOA (product C8) used at the plant (units of 1000 lbs year1), or the estimated annual PFOA emissions from the plant (Paustenbach et al., 2007; Shin et al., 2011a) (Fig. 1). All models were developed using the SAS mixed model procedure where the natural log of the serum PFOA concentration was the outcome, and a repeated measures covariance structure was used to account for the presence of multiple samples on some of the subjects within a job category. Several covariance structures were examined including compound symmetry, autoregressive, and unstructured. Due to the structure of the data, some covariance structures would not allow model convergence and there was little difference in the results of those structures that were fit. Therefore, a compound symmetry structure was used. Final model diagnostics from SAS mixed influence and residual analyses were examined and where highly influential values were identified by the press statistic and restricted likelihood distance, they were removed so that models would have maximum generalizability. For the fine powder and granular job group, four samples were removed and one person who moved between chemical operator and finisher jobs was removed (10 samples). For PFOA non-production jobs, two samples were removed and for maintenance jobs one sample was removed. To evaluate the impact of the modeling effort on variance components, mixed models with only the random effect were run to determine the baseline total, between- and within-worker variance estimates. Then the variance components for the models that included the fixed effects were determined and the percent change in the between- and within-worker variance calculated. Retrospective workplace PFOA serum concentrations 5 of 13 1029 Fig. 1. The annual amount of PFOA (product C8) used at the plant (lbs/year) and the estimated annual PFOA emissions from the plant (Paustenbach, et al. 2007 and Shin, et al. 2011a). Annual exposure estimates To estimate each worker’s work life exposure, for each work history record, a job category/job group assignment was made, cumulative prior years in potentially PFOA-exposed jobs were calculated for each year, and then the retrospective serum estimates were made for each job held in each year. For years when multiple jobs were held, an annual weighted average serum level was calculated as the timeweighted combination of the serum level in each job times the time in that job for the year. For people missing job information within a year, we used a background serum level of 0.03 ppm, based on a median value for community residents living near the plant in 2005, who were drinking water contaminated with PFOA (Steenland et al., 2009a). If an entire year is missing and the person had an annual exposure estimate greater than 0.03 in the previous year, then the data gap year is assigned 0.82 times the previous estimate; otherwise, the year is assigned 0.03. The decay to 82% of the previous year is based on the half-life estimate of 3.5 years (Olsen et al., 2007). When a worker left employment at the West Virginia plant, their serum level in the last year worked was allowed to decay in subsequent years using a half-life of 3.5 years. If they did not work for a full year during their last year, they were given the time-weighted average serum level using their working time serum level plus 0.82 times their work- ing time serum level for their non-working time during that year. RESULTS Prediction models The overall dataset had 2125 samples, covering 1308 workers from the years of 1979–2004. Twenty-six percent of the workers had two or more samples, with the range being 1–15 samples per person (Table 1). The serum concentration models were developed individually for each job category (Table 2). Where job groups within a job category showed significant differences in serum levels in the model, the job group was included as a predictor. Examples of job groups modeled within job categories included finish versus chemical operators within the fine powder/granular PTFE job category, monomer operators versus all other workers within the non-C8 use Teflon/co-polymer production job category, and Teflon/co-polymer versus general maintenance workers within the maintenance job category. With one exception (FEP/PFA), the models predicting serum levels all used either the amount of C8 used in the plant in the year the sample was collected or the amount of C8 emitted by the plant in the year the sample was collected as a predictor (r 5 0.57). Both these were significant positive variables with a consistent beta resulting in a 3% increase in serum PFOA level per increase in 1000 pounds of C8 emitted per year 0.007–4.14 0.06–6.81 0.50 (0.89) 0.16 (0.24) 0–30 0–22 3.6 (5.2) 1.9 (4.3) 2004 1982 1–3 1–4 7 200 463 246 504 Maintenance Non-Teflon/co-polymer production division jobs with no PFOA use 19 0.008–14.58 0.13–14.04 1.69 (2.53) 0.44 (0.84) 1–35 1–35 12.7 (8.4) 11.6 (8.0) 2000 1986 1–9 1–8 626 48 96 480 208 FEP/PFA Non-PFOA (C8) use jobs in Teflon and co-polymer production 20 0.007–59.40 0.09–59.40 0.58 (2.05) 2.88 (5.47) 0–36 1–36 8.5 (8.6) 11.8 (9.0) 1995 1985 1–15 1–15 1308 170 2125 541 Total Fine powder and granular PTFE 26 57 Median (mean) of serum PFOA samples (ppm) Min-max of cumulative years in potentially exposed jobs Mean (SD) of cumulative years in potentially exposed jobs Year by which 50% samples collected Number of samples per worker (min–max) Percent workers with 2 samples Number of serum samples Number of workers sampled S. Woskie, R. Gore and K. Steenland Job categories Table 1. Description of PFOA serum dataset used in retrospective modeling. All subjects were in their job category/job group for at least 1 year prior to serum sample. Range of serum PFOA samples (ppm) 6 of 13 1030 for models using the amount of C8 emitted as a covariate. For the model with the amount of C8 used in the plant, the positive beta resulted in an increase in serum PFOA level of 1% per 1000 pounds C8 used per year. The resulting models were used to predict serum PFOA for each worker in each year for the cohort. The earliest samples used to develop the models were in 1978 and 1979, although PFOA began use in 1951 at the West Virginia plant. For the job categories with direct exposure to C8 that used the RCS models to account for process changes over time, we set calendar year to 1979 for the pre-1979 years. This prevented overestimation of early serum levels where we do not know how process changes may have influenced serum concentrations. For these job categories, changes in serum levels predicted prior to 1979 are largely influenced by the annual amount of PFOA (product C8) used at the plant (units of 1000 lbs year1) as well as the individuals’ tenure in jobs with potential C8 exposure. Since the FEP/PFA operation started in 1953, but did not begin use of C8 until 1975, the level in 1974 was due to plant background exposures, not direct use. Therefore, the average for the intermittent direct/plant background-exposed jobs in the Teflon/co-polymer area (Teflon engineers, laboratory, supervisor, office workers, monomer operators, and Teflon/co-polymer maintenance) in 1974 was used to extrapolate back to a zero concentration in 1953 when FEP began pilot operations. For job categories with intermittent or plant background PFOA exposures, the estimated annual PFOA emissions from the plant were a significant predictor in the models. However, the models did not predict zero exposures in 1950 when C8 began use at the plant. Examination of the annual amount of C8 used/emitted from the plant showed a relatively linear increase until about 1964 when fluctuations in use/emission began (Fig. 1). Therefore, retrospective estimation was modified so that model predictions were used retrospectively through 1964, and then a linear extrapolation back to zero in 1950 was used for these job groups. Since the retrospective serum exposure estimates include an individual-level predictor (cumulative prior years in potentially exposed jobs), determining the pattern of serum levels over time required calculating the median serum exposure level for all workers in each job category/job group in a given year. This was done in 5-year increments for each job category/job group for those directly exposed to C8 (Fig. 2) and those with intermittent or plant background PFOA exposures (Fig. 3). For the fine powder/granular job category, the chemical operators run the reactors and were in charge of measuring Retrospective workplace PFOA serum concentrations 7 of 13 1031 Table 2. Mixed models using repeated measures for natural log of serum PFOA concentration (ppm). Separate mixed models job categories ln PFOA serum concentration (ppm) as outcome Beta Standard error P value Fine powder/granular PTFE job category (direct exposure to PFOA)a Intercept T1 time (number of years after 1970 when serum sample taken) T2 time (spline function) T3 Time (spline function) Cumulative prior years in potentially exposed jobs Annual amount of PFOA (C8) used (units 5 1000 lbs year�1) Chemical operator job group (0/1: yes 5 1) FEP/PFA job category (direct exposure to PFOA)b Intercept T1 time (number of years after 1970 when serum sample taken) T2 time (spline function) T3 time (spline function) Cumulative years in potentially exposed jobs 1.61 0.32 ,0.0001 �0.13 0.03 ,0.0001 0.46 0.63 0.46 �0.47 1.01 0.64 0.02 0.01 ,0.01 0.01 ,0.01 1.03 0.11 0.001 ,0.0001 2.81 0.39 ,0.0001 �0.23 0.03 ,0.0001 1.61 0.19 ,0.0001 �2.54 0.28 ,0.0001 0.01 ,0.0001 0.03 Non-PFOA (C8) use in Teflon and co-polymer production job category (intermittent direct or plant background PFOA exposure)c 0.08 ,0.0001 0.03 ,0.01 ,0.0001 Annual amount of PFOA (C8) emitted from plant (units 5 1000 lbs year�1) 0.03 ,0.01 ,0.0001 TFE monomer operator job group (0/1: yes 5 1) 0.94 0.12 ,0.0001 �1.30 0.11 ,0.0001 0.03 0.01 0.01 Annual amount of PFOA (C8) emitted from plant (units 5 1000 lbs year�1) 0.03 0.01 0.0001 Teflon/co-polymer maintenance job group (0/1: yes 5 1) 0.47 0.14 0.002 �2.24 0.05 ,0.0001 0.07 0.01 ,0.0001 0.03 ,0.01 ,0.0001 Intercept Cumulative prior years in potentially exposed jobs Maintenance job category (intermittent direct or plant background PFOA exposure)c Intercept Cumulative prior years in potentially exposed jobs Non-Teflon/co-polymer production division job category with no PFOA use (plant background PFOA exposure)d Intercept Cumulative prior years in potentially exposed jobs Annual amount of PFOA (C8) emitted from plant (units 5 1000 lbs year�1) �1.74 a Model for fine powder/granular PTFE job category: Yij 5b0 þ b1 Cum yearsij þ b2 C8 Usedij þ b3 Chem operatorij þ b4 T1 ðtÞij þ b5 T2 � ðtÞij þ b6 T3 � ðtÞij þ wi þ eij ; where Yij 5 the natural log of the serum concentration of the jth sample for the ith worker. The b’s are fixed effects, wi is the random effect for the worker and eij is the error term. T1 (t)ij 5 t, where t 5 number of years post-1970. ðt�11Þ 3 ðt�17Þ 3 ðt�34Þ 3 ðt�13Þ 3 ðt�17Þ 3 ðt�34Þ 3 T2 � t 5 kd � þ � 1:35 kd � þ þ 0:358 kd � þ T3 � t 5 kd � þ � 1:24 kd � þ þ 0:24 kd � þ 2 weighting for numerical stability5kd � 5ð34 � 11Þ3 ; ðt � cÞþ 5t � c if t � c . 0; elseðt � cÞþ 50; where c 5 knot in units of years post-1970. b Model for FEP/PFA job category: Yij 5b0 þ b1 Cum yearsij þ b4 T1 ðtÞij þ b5 T2 � ðtÞij þ b6 T3 � ðtÞij þ wi þ eij ; where Yij 5 the natural log of the serum concentration of the jth sample for the ith worker. The b’s are fixed effects, wi is the random effect for the worker and eij is the error term. T1 (t)ij 5 t, where t 5 number of years post-1970. ðt�9Þ 3 ðt�25Þ 3 ðt�34Þ 3 ðt�13Þ 3 ðt�25Þ 3 ðt�34Þ 3 T2 t 5 kd þ � 2:78 kd þ þ 1:78 kd þ T3 t 5 kd þ � 2:33 kd þ þ 1:33 kd þ 2 weighting for numerical stability5kd5ð34 � 9Þ3 ; ðt � cÞþ 5t � c if t � c . 0; elseðt � cÞþ 50; where c 5 knot in units of years post-1970. c Models for non-PFOA (C8) use in Teflon and co-polymer production job category or maintenance job category: Yij 5b0 þ b1 Cum yearsij þ b2 C8 Emittedij þ b3 jobgroupij þ wi þ eij ; where Yij 5 the natural log of the serum concentration of the jth sample for the ith worker. The b’s are fixed effects, wi is the random effect for the worker and eij is the error term. d Model for non-Teflon/co-polymer production division job category with no PFOA use: Yij 5b0 þ b1 Cum yearsij þ b2 C8 Emittedij þ wi þ eij ; where Yij 5 the natural log of the serum concentration of the jth sample for the ith worker. The b’s are fixed effects, wi is the random effect for the worker and eij is the error term. 8 of 13 1032 S. Woskie, R. Gore and K. Steenland Fig. 2 Median serum levels estimated from models for workers in job groups with direct exposure to C8 during that year. Fig. 3 Median serum levels estimated from models for workers in a job groups with intermittent direct or plant background PFOA exposure to C8 during that year. Models used to estimate levels until 1965 then from 1964 to 1950 linear decline from 1965 median level used. the powdered C8 (APFO) for introduction into the reactors had higher serum levels than the finish operators who work in the powder drying and packaging area and transfer liquid polymer dispersions into shipping totes. No significant differences were found between chemical and finish operators in the FEP/ PFA job category in part due to a staffing change in 1995 so that workers rotated between these jobs. We were unable to differentiate those workers in FEP operations and those in PFA operations; all were considered co-polymer operations. Within the nonPFOA (C8) use jobs in Teflon and co-polymer production job category, the TFE monomer operators had significantly higher serum levels than the other jobs including engineers, laboratorians, supervisors, and office workers. Within the maintenance job category, the Teflon/co-polymer maintenance workers had higher serum levels than general maintenance jobs. Evaluation of the change in the variance components attributable to the modeling effort are shown in Table 3. For all job categories, the betweenworker variance exceeded the within-worker variance both before and after inclusion of the fixed effects, including covariates that reflect calendar time changes in the plant operations. However, the inclusion of fixed effects explained 27–31% of the variance of the job categories with potential PFOA exposure and resulted in reductions in the betweenworker variance of 27–38% while also reducing the within-worker variance component. Predicted versus observed serum levels Comparison of the model predictions to sample data job category/job group overall and by decade found the paired predicted and sample values were not significantly different overall (average difference 0.08 ppm, P , 0.001) for about 75% of the 30 job/decade analyses. In the job/decade analyses where there was a significant difference, predicted geometric mean PFOA serum concentrations were lower than Retrospective workplace PFOA serum concentrations 9 of 13 1033 sample geometric mean (GM) concentrations in all but one case (Teflon/co-polymer maintenance in 1980–1989). Overall, the Spearman correlation coefficient for the predicted and measured PFOA concentrations was 0.8. As part of the settlement of a community class action lawsuit, blood samples were collected and measured for serum PFOA in Ohio and West Virginia residents who lived in water districts surrounding the Dupont chemical plant (Shin et al., 2011b). Of those residents, 2034 were Dupont workers in our cohort. We compared the serum PFOA concentration predicted by our models for 2004 with that measured during the community study (2005–2006). The overall Spearman correlation for this comparison was 0.6 with a median predicted of 0.11 ppm (mean 0.26 ppm) and a median measured of 0.11 ppm (mean 0.32 ppm). This comparison included 1025 subjects who were no longer working at the plant in 2004 (median predicted and measured 5 0.08 ppm). For the 1009 subjects working at the plant in 2004, the median predicted was 0.13 ppm (mean 0.36 ppm) and the median measured was 0.16 ppm (mean 0.47 ppm). Annual exposures The study cohort of over 6000 workers had over 56,000 work history records. The model-based retrospective serum estimates were made for each job held by a worker in each year. For years when multiple jobs were held, an annual weighted average serum level was calculated as the time-weighted combination of the serum level in each job times the time in that job for the year. Fig. 4 shows stem and leaf plots of weighted annual average serum PFOA levels (parts per million) for the full cohort of workers by those in job category/job groups with potential C8 exposure during each year and those in job categories without potential C8 exposure during each year. DISCUSSION Interpretation of the retrospective estimates The trends in predicted serum PFOA levels over time shown in Figs 2 and 3 reflect a complex plant history including changes in production levels, engineering controls, and personal protective equipment use instituted to decrease C8 exposures. For the fine powder/granular job category (finish operator and chemical operator job groups), the increases up until 1980 (Fig. 2) reflect a gradual increase in C8 use (Fig. 1). However, between 1980 and 2000, despite the increase in use, serum levels declined due to the incorporation of a number of exposure controls including a change to liquid injection pumps for addition of C8 into reactors, replacement of weighing powdered C8 with premixed liquid C8, addition of a dryer scrubber, use of personal protective equipment (Lin et al., 2005) when weighing C8, entering the dryers, cleaning tanks, changing dispersion filters or loading liquid dispersions into totes loading operations, and when changing filters in the dispersion operation (Fig. 2). For the FEP/PFA job category, the amount of C8 used in the plant was not a useful predictor in part because between 1982 and 1997, the PFA operation became a non-aqueous/non-C8 operation and in part because the proportion of C8 used in the plant by this operation changed as production increased with the addition of a new reactor and the new bead facility. The lower serum levels in the 1980s appear to be a result of the elimination of C8 in the PFA operation as well as the change to premixed C8 and the sealing of the dryer room for the FEP operation. Higher levels Table 3. Variance components for job category serum PFOA models. Model with only worker as random effect Full models from Table 2a Total variance r2T % Between-worker variance r2B % Within-worker variance r2W % Between-worker variance r2B % Within-worker variance r2W % Fixed effect variance r2F Fine powder/granular 1.55 81 19 58 11 31 FEP/PFA 1.05 78 22 58 13 30 Non-PFOA use in Teflon and co-polymer production Maintenance 1.14 66 34 41 33 27 0.92 58 42 41 33 27 Non-Teflon/co-polymer production 0.82 73 27 68 13 18 Job category a The total variance from the model containing only worker as a random effect was used to calculate the contribution of fixed effects to partitioning of the variance components. 10 of 13 1034 S. Woskie, R. Gore and K. Steenland Fig. 4 Stem and leaf plots of the model estimated weighted annual average serum PFOA levels (ppm) for the full cohort of workers by those in job category/job groups with potential C8 exposure during the year (A) and those in job categories without potential C8 exposure during the year (B). Box plots show the minimum and maximum on the whiskers, the 25 th and 75 th percentile in the box and the bar in the box is the median. in the 1990s can be attributed to the startup of the bead facility, the addition of a new reactor, and the re-addition of C8 into the PFA operation. The plateau of levels in the late 1990s and 2000s may reflect the addition of dryer scrubbers that offset the increasing production levels. For workers in job groups with intermittent or plant background PFOA exposures (Fig. 3), serum levels reached a peak in 2000 that corresponded to the period of highest use of C8 in the plant (Fig. 1), the subsequent decrease reflects lower use and also the addition of deep bed scrubbers in the dryer operations. The Teflon polymer/co-polymer maintenance job group within the maintenance job category had a peak level in 1980 but that is the one case where the exposure models predicted higher levels than the serum sample data. For the nonPFOA (C8) Teflon polymer/co-polymer production job category, the job group of TFE monomer operator had higher exposures than other jobs with intermittent exposures to C8 during visits to the production areas or use in QC/QA procedures including engineers, laboratorians, higher level supervisors, and clerks. This is likely due to the location of the monomer operation adjacent to the fine powder/ granular operations, with the monomer control room directly above the fine powder dryers. Little is known about the toxicokinetics of PFOA in humans. In this study, we found a consistent 2–3% increase in serum PFOA level per cumulative prior years of work in a potentially PFOA-exposed job for all models. This suggests that PFOA may result in a significant body burden in chronically exposed workers. The only other report of body build up of PFOA is in the study of community members exposed through contaminated water supplies. An average of 1% increase in serum PFOA per year of residence in an exposed water district has been reported for the area of the plant (Seals et al., 2011). It has been suggested that the longer elimination half-life for PFOA in humans than in other species may be due to differences in biliary excretion and gut resorption (Olsen et al., 2007). These mechanisms may also explain the potential for accumulation of PFOA in humans. Comparison of the serum levels measured here with other datasets is difficult. The occupational data from studies of 3M employees have been reported cross-sectionally for the years of data collection, which ranged from 1993 to 2000 (Olsen et al., 2000; Olsen and Zobel, 2007). The data reported here cover the time period of 1979–2004. Samples collected for this study from 2000 to 2004 represent all job categories and have a median serum level of 0.24 ppm (mean 0.77 ppm) with a range of 0.007–16.20 ppm. Samples from a 2000 survey of three 3M plants where PFOA was manufactured and used in fluoroelastomer production had a median of 1.10 ppm (mean 2.21 ppm with a range of 0.01– 92.03 ppm; Olsen and Zobel, 2007). Lower overall levels in the present study may be a reflection of the high percentage of workers either unexposed (27%) or intermittently exposed (42%) in the 2000–2004 samples (n 5 731). While at the 3M plants, about 75% of the males at two of the plants were reported to be in production jobs (Olsen et al., 2003). Workers unexposed to C8 in this study were those in non-Teflon and non-co-polymer production Retrospective workplace PFOA serum concentrations 11 of 13 1035 divisions and business staff. Their overall exposures were a median of 0.16 ppm (mean of 0.24) with a range of 0.007–4.14 ppm. In the 2000–2004 period levels, this job category had a GM serum level of 0.14 ppm. These values are higher than those reported for the local community, who had a median serum level in 2005–2006 of 0.04 ppm for current residents of any water district in the area of the plant (Steenland et al., 2009a) and higher than the National Health and Nutrition Examination Survey results for the general US population from the 1999–2000 or 2003–2004 surveys which ranged from 0.004 to 0.005 ppm for subjects aged 20–39 and 40–59 years (Kato et al., 2011). Retrospective exposure assessment using biomarkers Although much has been said regarding the future of molecular epidemiology and the utility of biomarkers in accounting for all routes of exposure, few studies have been able to use biomarkers to look at health effects caused by chronic exposures. In part, this is due to the paucity of biomarker data for target cohorts over the extended periods of time necessary to account for disease latency. External exposure data are more commonly available. Where there is matched external and biomarker data, toxicokinetic models can be used to retrospectively estimate internal dose (biomarker level) from intake of external exposures (Steenland et al., 2001). This sub-cohort of workers represents a unique situation with over 2100 samples, covering over 1300 workers over a 25-year period with 26% of the workers having two or more samples. In addition, detailed work histories were available for this sub-cohort and the full cohort of over 6000 people. This enabled the development of retrospective models for serum levels based on time and job which could be applied to the full cohort. A unique feature of this approach is the use of RCS models to account for non-linear and non-step changes in serum levels over time that are reflective of the variation in fluctuating usage/ emissions of C8 by the plant and the progressive introduction of engineering and work practice controls over time. The X values (T1–3) from the RCS function were then used in linear mixed models to develop terms to use with calendar time for the estimation of historical exposures for the cohort. RCS models in conjunction with linear mixed effect models were also used to evaluate atrazine exposures among pesticide applicators (Hines et al., 2006). However, those models were not used for developing retrospective exposure estimates. The modeling described here found more between-worker than within-worker variance. This has also been noted by the few other occupational studies that have examined this topic for more persistent chemicals like mercury and lead (Symanski and Greeson, 2002). Several authors have described how adding fixed effects to models alter variance components for air samples, but this has not been widely done for biomarkers (Peretz et al., 2002; McClean et al., 2004). In this study, the fixed effects including covariates that reflect calendar time changes in the plant operations explained 27–31% of the variance of the jobs with potential PFOA exposure. One measure that has been used to describe the potential bias introduced when modeling exposure response relationships is the estimated within-worker to be tween-worker variance ratio k5r2W r2B . When k is smaller, the potential bias is less. A study of 127 datasets found that biomarker data had a lower median variance ratio (1.04) than air sampling data (2.40), suggesting that biomarker can provide a less biased surrogate for exposures than typical air sampling (Lin et al., 2005). In our own data, lambda ranged from 0.20 to 0.80 across the five regression models. Use of biomonitoring data collected for surveillance purposes to model retrospective exposures has inherent limitations. In this study, the frequency of monitoring varied over time and exposure category and although attempts at adjustment have been made, the methods of analysis changed over time as well. Nevertheless, we believe the approach described here represents an improvement in the specificity of exposure characterization for this cohort. We believe that the use of smoothing better accounts for the gradual changes in serum levels over time due to process changes; however, it cannot predict levels beyond the range of the data. Thus, for the direct exposure jobs where we lack data, we have presumed that process changes (our interpretation the spline function) did not have a significant impact on exposures before 1979 (earliest serum data). Nevertheless, there remained some aspect of time in the prediction models since the amount of C8 used in the plant did remain an important predictor of exposure before 1979 for the fine powder/granular PTFE jobs. In addition, the models described here do not include important personal characteristics that may influence serum levels. For example, studies have reported that age, race, and gender are associated with PFOA serum levels (Calafat et al., 2007; Steenland et al., 2009a). Although we do not include age in our models, to some extent, age is implicitly included since workers tended to stay at the plant and all serum prediction models included a covariate for cumulative years in potentially exposed jobs. 12 of 13 1036 S. Woskie, R. Gore and K. Steenland Another influence on worker serum levels may have come from personal exposures via water in communities surrounding the plant (Emmett et al., 2006; Steenland et al., 2009a). The exposure estimates reported here do not explicitly account for residential exposures over time, although it is believed that relative to workplace exposures these are relatively small. For example, current workers were reported to have a median serum PFOA level of 0.147 versus 0.074 ppm for former workers and 0.028 ppm for current/former residents in the study of the nearby community member PFOA levels (Steenland et al., 2009b). To evaluate the health effects of chronic exposure to PFOA in this cohort, it would be ideal to have personal biomarker data for the full historical period for each cohort member. However, since that is not possible, a method to retrospectively estimate exposures with minimal misclassification is needed. We believe the method described here has strengths in accommodating gradual changes over time in serum PFOA levels by tying exposure levels to process changes and plant usage patterns and maximizing the number of job groups for which historical estimates are made. Our comparisons of predicted with internal and external measurement data were favorable (Spearman correlations of 0.8 and 0.6, respectively), suggesting the models perform well in predicting worker exposures. The annualized PFOA serum estimates developed here will be used in a series of morbidity and mortality studies of this worker cohort. FUNDING C8 Class Action Settlement Agreement (Circuit Court of Wood County, West Virginia) between DuPont and Plaintiffs. Acknowledgements—Many thanks to the Teflon employees and retirees who assisted in gathering historical information to help us understand historical process changes and how to assign work record department/division/job combinations to job groups similarly exposed to C8. Thanks to Dr James Deddens for advice and assistance with the RCS approach and SAS Macro. This research was funded by the C8 Class Action Settlement Agreement (Circuit Court of Wood County, West Virginia) between DuPont and Plaintiffs, which resulted from releases into drinking water of the chemical PFOA (C8). Funds are administered by an agency that reports to the court. The work and conclusions are independent of either party to the lawsuit. REFERENCES Armstrong BG. (1990) The effects of measurement errors on relative risk regressions. Am J Epidemiol; 132: 1176–84. Blair A, Stewart P, Lubin JH et al. (2007) Methodological issues regarding confounding and exposure misclassification in epidemiological studies of occupational exposures. Am J Ind Med; 50: 199–207. Calafat AM, Kuklenyik Z, Reidy JA et al. (2007) Serum concentrations of 11 polyfluoroalkyl compounds in the U.S. population: data from the national health and nutrition examination survey (nhanes). Environ Sci Technol; 41: 2237–42. Dosemeci M, Wacholder S, Lubin JH. (1990) Does nondifferential misclassification of exposure always bias the true effect toward the null? Am J Epidemiol; 105: 746–8. Ehresman DJ, Froehlich JW, Olsen GW et al. (2007) Comparison of human whole blood, plasma, and serum matrices for the determination of perfluorooctanesulfonate (pfos), perfluorooctanoate (pfoa), and other fluorochemicals. Environ Res; 103: 176–84. Emmett EA, Zhang H, Shofer FS et al. (2006) Community exposure to perfluorooctanoate: relationships between serum levels and certain health parameters. J Occup Environ Med; 48: 771–9. Han X, Snow TA, Kemper RA et al. (2003) Binding of perfluorooctanoic acid to rat and human plasma proteins. Chem Res Toxicol; 16: 775–81. Harrell FE. (2001) Regression modeling strategies. New York, NY: Springer-Verlag. Harrell FE Jr, Lee KL, Pollock BG. (1988) Regression models in clinical studies: determining relationships between predictors and response. J Natl Cancer Inst; 80: 1198–202. Hines CJ, Deddens JA, Lu C et al. (2006) Mixed-effect models for evaluating multiple measures of atrazine exposure among custom applicators. J Occup Environ Hyg; 3: 274– 83. Jones PD, Hu W, De Coen W et al. (2003) Binding of perfluorinated fatty acids to serum proteins. Environ Toxicol Chem; 22: 2639–49. Karrman A, van Bavel B, Jarnberg U et al. (2006) Perfluorinated chemicals in relation to other persistent organic pollutants in human blood. Chemosphere; 64: 1582–91. Kato K, Basden BJ, Needham LL et al. (2011) Improved selectivity for the analysis of maternal serum and cord serum for polyfluoroalkyl chemicals. J Chromatogr A; 1218: 2133–7. Kreckmann KH, Sakr CJ, Leonard RC et al. (2009) Estimation and validation of biomarker-based exposures for historical ammonium perfluorooctanoate. J Occup Environ Hyg; 6: 511–6. Leonard RC, Kreckmann KH, Sakr CJ et al. (2008) Retrospective cohort mortality study of workers in a polymer production plant including a reference population of regional workers. Ann Epidemiol; 18: 15–22. Lin YS, Kupper LL, Rappaport SM. (2005) Air samples versus biomarkers for epidemiology. Occup Environ Med; 62: 750–60. Lundin JI, Alexander BH, Olsen GW et al. (2009) Ammonium perfluorooctanoate production and occupational mortality. Epidemiology; 20: 921–8. McClean MD, Rinehart RD, Ngo L et al. (2004) Inhalation and dermal exposure among asphalt paving workers. Ann Occup Hyg; 48: 663–71. Olsen GW, Burris JM, Burlew MM et al. (2000) Plasma cholecystokinin and hepatic enzymes, cholesterol and lipoproteins in ammonium perfluorooctanoate production workers. Drug Chem Toxicol; 23: 603–20. Olsen GW, Burris JM, Burlew MM et al. (2003) Epidemiologic assessment of worker serum perfluorooctanesulfonate (pfos) and perfluorooctanoate (pfoa) concentrations and medical surveillance examinations. J Occup Environ Med; 45: 260–70. Retrospective workplace PFOA serum concentrations 13 of 13 1037 Olsen GW, Burris JM, Ehresman DJ et al. (2007) Half-life of serum elimination of perfluorooctanesulfonate, perfluorohexanesulfonate, and perfluorooctanoate in retired fluorochemical production workers. Environ Health Perspect; 115: 1298–305. Olsen GW, Zobel LR. (2007) Assessment of lipid, hepatic, and thyroid parameters with serum perfluorooctanoate (pfoa) concentrations in fluorochemical production workers. Int Arch Occup Environ Health; 81: 231–46. Paustenbach DJ, Panko JM, Scott PK et al. (2007) A methodology for estimating human exposure to perfluorooctanoic acid (pfoa): a retrospective exposure assessment of a community (1951-2003). J Toxicol Environ Health A; 70: 28–57. Peretz C, Goren A, Smid T et al. (2002) Application of mixedeffects models for exposure assessment. Ann Occup Hyg; 46: 69–77. Sakr CJ, Kreckmann KH, Green JW et al. (2007a) Cross-sectional study of lipids and liver enzymes related to a serum biomarker of exposure (ammonium perfluorooctanoate or apfo) as part of a general health survey in a cohort of occupationally exposed workers. J Occup Environ Med; 49: 1086–96. Sakr CJ, Leonard RC, Kreckmann KH et al. (2007b) Longitudinal study of serum lipids and liver enzymes in workers with occupational exposure to ammonium perfluorooctanoate. J Occup Environ Med; 49: 872–9. Sakr CJ, Symons JM, Kreckmann KH et al. (2009) Ischaemic heart disease mortality study among workers with occupational exposure to ammonium perfluorooctanoate. Occup Environ Med; 66: 699–703. Seals R, Bartell SM, Steenland K. (2011) Accumulation and clearance of perfluorooctanoic acid (pfoa) in current and former residents of an exposed community. Environ Health Perspect; 119: 119–24. Shin H-M, Vieira V, Ryan PB et al. (2011a) Environmental fate and transport modeling for perfluorooctanoic acid emitted from the Washington works facility in West Virginia. Environ Sci Technol; 45: 1435–42. Shin H-M, Vieira VM, Ryan PB et al. (2011b) Retrospective exposure estimation and predicted versus observed serum perfluorooctanoic acid concentrations for participants in the c8 health project. Environ Health Perspect; 119: 1760–1765. Steenland K, Deddens J, Piacitelli L. (2001) Risk assessment for 2,3,7,8-tetrachlorodibenzo-p-dioxin (tcdd) based on an epidemiologic study. Am J Epidemiol; 154: 451–8. Steenland K, Deddens JA, Zhao S. (2000) Biases in estimating the effect of cumulative exposure in log-linear models when estimated exposure levels are assigned. Scand J Work Environ Health; 26: 37–43. Steenland K, Fletcher T, Savitz DA. (2010) Epidemiologic evidence on the health effects of perfluorooctanoic acid (pfoa). Environ Health Perspect; 118: 1100–8. Steenland K, Jin C, MacNeil J et al. (2009a) Predictors of pfoa levels in a community surrounding a chemical plant. Environ Health Perspect; 117: 1083–8. Steenland K, Tinker S, Frisbee S et al. (2009b) Association of perfluorooctanoic acid and perfluorooctane sulfonate with serum lipids among adults living near a chemical plant. Am J Epidemiol; 170: 1268–78. Symanski E, Greeson NM. (2002) Assessment of variability in biomonitoring data using a large database of biological measures of exposure. AIHA J (Fairfax, Va); 63: 390–401. US EPA. (2005) Draft risk assessment of the potential human health effects associated with exposure to perfluorooctanoic acid and its salts. Washington, DC: U.S. Environmental Protection Agency, Office of Pollution Prevention and Toxics, Risk Assessment Division. Available at http:// www.epa.gov/oppt/pfoa/pubs/pfoarisk.pdf. Accessed 20 October 2011. Wu LL, Gao HW, Gao NY et al. (2009) Interaction of perfluorooctanoic acid with human serum albumin. BMC Struct Biol; 9: 31.
© Copyright 2026 Paperzz