Annals of Occupational Hygiene Advance Access originally published online on August 3, 2004
Annals of Occupational Hygiene 2004 48(6):565-578; doi:10.1093/annhyg/meh044
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
© 2004 British Occupational Hygiene Society Published by Oxford University Press;
Urinary 1-Hydroxypyrene and Polycyclic Aromatic Hydrocarbon Exposure Among Asphalt Paving Workers
1 Harvard School of Public Health, Boston, MA 02115, USA; 2 Occupational Safety and Health Administration, Washington, DC 20210, USA; 3 University of CaliforniaSan Francisco, San Francisco, CA 94143, USA
* Author to whom correspondence should be addressed. e-mail: mmcclean{at}hsph.harvard.edu
Received 6 January 2003; in final form 30 December 2003
| ABSTRACT |
|---|
|
|
|---|
Objectives: Using urinary 1-hydroxypyrene (1-OHP) as a measure of total absorbed dose, the primary objective of this study was to evaluate the total effect of inhalation and dermal PAH exposures while considering other factors such as age, body mass index and smoking that may also have a significant effect on urinary 1-OHP.
Methods: The study population included two groups of highway construction workers: 20 paving workers and 6 milling workers. During multiple consecutive workshifts, personal air and dermal samples were collected from each worker and analyzed for pyrene. During the same work week, urine samples were collected pre-shift, post-shift and at bedtime each day and analyzed for 1-OHP. Distributed lag models were used to evaluate the independent effect of inhalation and dermal exposures that occurred at each of several preceding exposure periods and were used to identify the relevant period of influence for each pathway.
Results: The paving workers had inhalation (mean 0.3 µg/m3) and dermal (5.7 ng/cm2) exposures to pyrene that were significantly higher than the milling workers. At pre-shift on Monday morning, following a weekend away from work, the pavers and millers had the same mean baseline urinary 1-OHP level of 0.4 µg/g creatinine. The mean urinary 1-OHP levels among pavers increased significantly from pre-shift to post-shift during each work day, while the mean urinary 1-OHP levels among millers varied little and remained near the baseline level throughout the study period. Among pavers there was a clear increase in the pre-shift data during the work week, such that the average pre-shift level on day 4 (1.4 µg/g creatinine) was 3.5 times higher than the average pre-shift results on day 1 (0.4 µg/g creatinine). The results of the distributed lag model indicated that the impact of dermal exposure was approximately eight times the impact of inhalation exposure. Furthermore, dermal exposure that occurred during the preceding 32 h had a statistically significant effect on urinary 1-OHP, while the effect of inhalation exposure was not significant.
Conclusions: We found that distributed lag models are a valuable tool for analyzing longitudinal biomarker data and our results indicate that dermal contact is the primary route of exposure to PAHs among asphalt paving workers. An exposure assessment of PAHs that does not consider dermal exposure may considerably underestimate cumulative exposure and control strategies aimed at reducing occupational exposure to asphalt-related PAHs should include an effort to reduce dermal exposure.
Keywords: asphalt dermal exposure distributed lag models hydroxypyrene paving pyrene
| INTRODUCTION |
|---|
|
|
|---|
Asphalt is a semi-solid residue that results from the non-destructive distillation of crude petroleum oil and contains a complex mixture of polycyclic aromatic hydrocarbons (PAHs), some of which are either known or suspected to be carcinogenic (Partanen et al., 1995
300 000 workers (Asphalt Institute, 1990
Exposure to PAHs can occur due to occupational sources (asphalt industry, coke plants, petrol refineries, aluminum industry) and non-occupational sources (smoking, diesel exhaust, grilled food) and can occur via inhalation, ingestion and dermal contact (Roggi et al., 1997
). Historically, the assessment of occupational exposure to PAHs has relied primarily on air monitoring (Jongeneelen, 2001
). However, there has been increasing evidence that dermal contact is the primary route of exposure to PAHs (Van Rooij et al., 1993a
,b
; Quinlan et al., 1995
; Viau and Vyskocil, 1995
; Borak et al., 2002
).
The mounting evidence regarding the importance of dermal PAH exposure has largely been due to the increased use of biological monitoring, which typically evaluates the relationship between urinary 1-hydroxypyrene (1-OHP) and environmental pyrene exposures. Pyrene is a useful indicator of environmental PAH exposure since it is abundant in most PAH mixtures and is strongly correlated with total PAHs in air (r = 0.87) and dermal (r = 0.65) samples (McClean et al., 2004
). Pyrene is absorbed through the gastrointestinal tract, lungs and/or skin and then metabolized to 1-OHP, which is ultimately excreted in urine. Urinary 1-OHP has become the most commonly used biological marker of PAH exposure, providing a convenient measure of total absorbed dose.
Urinary 1-OHP and other occupational biomarkers are often evaluated using modeling approaches that examine only the effect of the most recent environmental exposure. However, these approaches require the effect to be fully explained by the most recent exposure without considering the possible effect of preceding exposures. This can be problematic since occupational exposures are often serially correlated, such that the most recent exposure is actually capturing some information about previous exposures. Therefore, because the effects of previous exposures are not directly considered, the independent effect of the most recent exposure is overestimated while the total effect of exposure is underestimated (Pope and Schwartz, 1996
).
As an alternative method of evaluating biomarker data, this study uses distributed lag models to evaluate the effect of inhalation and dermal exposures at each of several preceding exposure periods, thereby resulting in a more accurate estimate of the total effect of each exposure pathway. Though distributed lag models have mainly been used in econometrics and social sciences, this analytical technique has recently been applied to epidemiological studies of air pollution (Pope and Schwartz, 1996
; Schwartz, 2000
; Braga et al., 2001
).
Using urinary 1-OHP as a measure of total absorbed dose, the objectives of this study were to: (i) evaluate the total effect and relative contribution of inhalation and dermal PAH exposures; (ii) assess the effect of inhalation and dermal exposure over different time periods to determine how differences in absorption and/or metabolism may affect total absorbed dose; (iii) evaluate the effect of other factors such as age, body mass index (BMI) and smoking status that may also influence urinary 1-OHP measurements. Understanding how inhalation and dermal exposure to PAHs affects total absorbed dose has significant implications with respect to designing control strategies in the paving industry, as well as to conducting future epidemiological studies of asphalt workers.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Study population
The study population included 26 highway construction workers, 20 pavers and six millers. The milling workers used grinding equipment to remove layers of asphalt from old roads in preparation for the paving crews, who later applied hot-mix asphalt while resurfacing the roads. All participants were male, worked for the same company and lived in the Greater Boston area. Written and informed consent was obtained from each study participant prior to sampling and all sampling was conducted in accordance with a standardized human subjects protocol that was approved by the Institutional Review Board at the Harvard School of Public Health.
The details of the paving operation are described elsewhere (McClean et al., 2004
). Briefly, each paving crew consisted of 68 workers performing four different tasks. The paver operator sat between the hopper and the screed while driving the paving machine. Two screed men stood on a platform attached to the back of the paving machine, each controlling one side of the screed, while 23 rakers followed closely using hand tools to fill holes and gaps. One or two rolling machines, each with its own operator, were used to smooth and compact the laid asphalt. The study population included workers from each of the four task categories (paver operators, screed men, rakers and roller operators) in three different crews.
The details of the milling operation are also described elsewhere (McClean et al., 2004
). Briefly, the milling crew consisted of six workers who used a large grinding machine, a smaller trimming machine and a bucket loader to grind away the surface of existing asphalt roads and deposit the ground asphalt into the back of a dump truck. The milling crew was included as a reference group since these highway construction workers do not work with hot-mix asphalt and have very low exposures to pyrene.
Study design
In May and June 1999, multiple exposure samples were obtained from each worker during full work shifts at job sites located within 1 h of Boston. Each day of sampling included the collection of personal air samples (particulate and vapor), dermal patch samples and urine samples. Questionnaire data regarding personal characteristics (e.g. age, height, weight) and non-occupational PAH exposure (e.g. smoking status, dietary information) were also obtained.
Personal air samples
Personal air samples were collected from each worker in accordance with NIOSH Method 5506 (NIOSH, 1998
). The air samples were obtained from each worker on three consecutive days for pavers and two consecutive days for millers. The air sampling system consisted of a Teflon filter and cassette holder to collect particulate PAHs, connected in series with an XAD-2 sorbent tube to collect the organic vapor. Using a personal sampling pump operating at 2 l/min, a 37 mm diameter filter (laminated PTFE of 2 µm pore size) was placed in a cassette and attached to each worker's lapel near the breathing zone. The sorbent tube containing XAD-2 was attached in line and downstream from each filter cassette. Flow rates were checked before, during and after sample collection using a calibrated rotameter. Opaque filter cassettes and foil-wrapped sorbent tubes were used to prevent sample degradation due to sunlight. Samples were transported in coolers and stored at 20°C.
The filter (particulate) and tube (vapor) samples were analyzed separately, each extracted with 4 ml of hexane and ultrasonicated for 60 min. Two milliliters were syringe-filtered, transferred to a clean tube and 2 ml of DMSO was added. The mixture was tumbled overnight and the layers were transferred to separate tubes. The DMSO layer was then analyzed for pyrene using high pressure liquid chromatography (HPLC).
Dermal patch samples
Personal dermal patch samples were collected from the wrists of each worker. The dermal samples were obtained from each worker on three consecutive days for pavers and two consecutive days for millers. The skin sampling method was a modification of the method described by Jongeneelen et al. (1988)
and Van Rooij et al. (1993c)
. A soft polypropylene adsorbent material (40 mm diameter) was attached to an exposure pad to create a dermal patch with an effective surface area of 8·71 cm2. Using an adhesive backing, the patches were attached to the underside of each wrist and resulted in the collection of two samples per worker per day. Following sample collection, the exposure pads were placed in foil-wrapped Petri dishes, transported in coolers and stored at 20°C.
The dermal patch samples were cut with a 33·3 mm punch, transferred to labeled culture tubes and 4 ml of DMSO was added. The tubes were then capped and sonicated for 1 h. Following sonication, 3 ml of the extract were transferred to a clean culture tube and 3 ml of hexane was added. The mixture was tumbled overnight, the layers were transferred to separate tubes and the DMSO layer was analyzed for pyrene using HPLC.
Urine samples
Three urine samples (pre-shift, post-shift and bedtime) were obtained on consecutive days and an additional pre-shift sample was obtained on the morning of the final day. This resulted in the collection of 10 urine samples from each paving worker and seven urine samples from each milling worker, with the first sample collected at pre-shift on Monday morning, with
8 h between each sample. All urine samples were collected in sterilized polypropylene specimen containers, transported in coolers and then frozen at 20°C until analysis.
Urines samples were analyzed for 1-OHP by HPLC using a modification of the methods described by Jongeneelen et al. (1987)
and Tolos et al. (1989)
. A 10 ml aliquot of urine was adjusted to pH 5.0 with 1.0 M acetic acid. After 10 ml of CH3COONa buffer and 25 µl of ß-glucuronidase (Sigma G-7770) were added, the mixture was incubated overnight at 37°C. The next day, Sep-Pak C18 cartridges (Waters, lot W8189b1) were primed with 5 ml of methanol and 10 ml of H2O and the hydrolyzed urine samples were passed through the cartridges. The retained solutes were eluted with 8 ml of methanol and the solvent was evaporated at 40°C using a constant flow of nitrogen (150 ml/min). The residue was dissolved in 2 ml of methanol, transferred to amber HPLC auto sampler vials and stored at 20°C until analysis.
The method was evaluated by pooling urine collected from non-exposed individuals (non-smoking graduate students), splitting it into three 1 l amber glass jars, spiking one of the jars with 1 ng/ml 1-OHP and another with 10 ng/ml. The third jar was left unspiked. To determine method recovery and precision, these control urines were run together with each batch of unknown samples. Based on the 21 control urines analyzed at each concentration, the average recoveries were 95 and 89%, while the relative standard deviations were 12 and 14% for the 1 and 10 ng/ml control urines, respectively.
Because HPLC may respond differently to 1-OHP in urine than to 1-OHP in methanol, a calibration curve was constructed with the pooled urine using the method of standard additions to minimize matrix interference. An intermediate standard was serially diluted in the pooled urine to prepare nominal 1-OHP concentrations of 0, 1, 5, 10, 15 and 20 ng/ml. The resulting data produced a linear standard curve with an intercept close to 0.0. The instrumental limit of detection was
0.50 ng/ml, calculated as 3 x SD of the smallest standard. All urine samples were analyzed and corrected for creatinine prior to statistical analysis.
Data analysis
All statistical analyses were conducted using SAS statistical software (SAS Institute, Cary, NC) and statistical significance is reported at the 0.05 level. The pyrene data (air and dermal samples) and 1-OHP data (urine) were analyzed using descriptive statistics, graphical displays and distributed lag models for repeated measures data. Following an algebraic transformation of the inhalation and dermal data, the MIXED procedure was used to evaluate the urinary 1-OHP data.
Air and dermal samples
Units for total air exposures are reported in µg/m3. One total air exposure estimate was calculated for each worker on each sampling day by adding the particulate and vapor measurements. Units for dermal exposure are reported in ng/cm2. One dermal exposure estimate was calculated for each worker on each sampling day by averaging the left and right wrist measurements. When only one wrist measurement was available, the result from that sample was used in place of the average.
Method limits of detection (LOD) were estimated as three times the SD of the field blanks. In cases where the mean field blanks were significantly different from zero (
= 0.05), the corresponding data were corrected by subtracting the mean field blank from each sample. Negative values that resulted from the subtraction of field blanks were included in all analyses, as were all values less than the detection limits (Gilbert, 1987
).
Six individual dermal samples collected from two workers were excluded from analysis (prior to averaging wrists) because there was sufficient evidence to suspect that the samples had been contaminated with diesel fuel. Both workers were rakers on the same paving crew who were observed to be purposefully contaminating the dermal patches with the fuel. These six patches were visibly discolored due to saturation and the resulting values were orders of magnitude higher than all other samples.
To evaluate the correlation between the total air and dermal data, traditional methods of estimating correlation coefficients (i.e. Pearson and Spearman) could not be used due to the repeated measures design of the study. Use of these traditional methods would erroneously ignore the number of subjects as the correct sample size while instead using the total number of observations as the incorrect sample size, thereby increasing the degrees of freedom (Hamlett et al., 2003
). As an alternative, all correlation coefficients were estimated using linear mixed effects models as described by Hamlett et al. (2003)
.
Urine samples
Units for urinary 1-OHP data are reported in µg/g creatinine. Since highly dilute or highly concentrated urine samples are generally not suitable for analysis, urine samples in which the creatinine concentration was either <0.3 or >3.0 g/l were excluded from the data analysis (ACGIH, 2002
). These criteria resulted in the exclusion of 10 urine samples.
ShapiroWilks tests and graphical displays indicated that the urinary 1-OHP data were not normally distributed. However, a log transformation of the data did result in an approximately normal distribution. Accordingly, all modeling was conducted using the log transformed urinary 1-OHP data.
Distributed lag models
Distributed lag models allow us to evaluate the possibility that each urinary 1-OHP measurement is not only influenced by the most recent inhalation and dermal exposures, but also by exposures that occurred during previous exposure periods. For notational and conceptual purposes, time t will refer to time intervals of 8 h each; inhalation and dermal samples were collected during each time interval, whereas urine samples were collected at the end of each time interval. Accordingly, the unconstrained lag model evaluating q + 1 exposure periods for the ith subject at time interval t is described as:
![]() |
it represents the residual error term. The total effect of exposure X is the sum of its impact during the most recent period and its impact on previous days, i.e. ß0 + ß1 + ... + ßq. However, because Xit is correlated with Xit1, Xit2,..., Xitq, the high degree of collinearity in this unconstrained model would result in unstable estimates of the values of ßj. To reduce the degree of collinearity and as an alternative to the unconstrained model, a polynomial distributed lag function can be used to constrain ßj to a polynomial function of lag j:
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
The distributed lag models required some assumptions regarding inhalation and dermal exposures that occurred away from work. Though inhalation and dermal measurements were obtained during the 8 h work shifts (corresponding to the time periods preceding the post-shift urine samples), exposure measurements were not available during the remaining 16 h of the day (corresponding to bedtime and pre-shift samples). Since the workers did not have occupational exposure during this time, inhalation and dermal exposures to pyrene were assumed to be 0. Similarly, since the first urine sample was always collected on Monday morning, the same rationale was used to assume 0 exposure during the 56 h prior to collection of the Monday morning pre-shift sample (representing the most recent exposure period and six lagged exposure periods during the weekend).
Linear mixed effects models
Since multiple urine samples were collected from each worker, linear mixed effects models were used to evaluate the distributed lag effect of inhalation and dermal exposure. This relationship was evaluated while controlling for the effect of age, BMI and smoking status (current smoker versus non-smoker). The repeated measures design and the use of linear mixed effects models allowed each worker to serve as their own control while evaluating changes in urinary 1-OHP over time. If we allow
j to represent the effect of inhalation exposure (Iij) and
j to represent the effect of dermal exposure (Dij), then the final linear mixed effects model can be described as follows:
![]() |
values represent the fixed effects for AGE (yr), BMI (kg/m2) and SMOKER (yes/no). Inhalation and exposures were initially evaluated using six lag periods. However, based on the observed results it was clear that the evaluation of six lag periods was not reasonable in all cases. For instance, dermal inhalation exposure was shown to have a positive effect on urinary1-OHP out to a maximum of three lag periods. To avoid the inclusion of irrelevant lag periods, the number of lags included in the final model for each exposure pathway was adjusted according to the biological plausibility of the observed results. Models were fitted using a heterogeneous autoregressive covariance matrix which yielded the lowest Akaike's Information Criteria diagnostic values compared with the other covariance matrices examined (e.g. independence, compound symmetry and homogeneous autoregressive). | RESULTS |
|---|
|
|
|---|
Table 1 presents the summary statistics for the pyrene exposure data. The summary statistics are stratified by exposure pathway and for paving workers are also stratified by task. In total air samples, pyrene was detected above the LOD of 0.01 µg/m3 in 97% of the paving samples (mean 0.3 µg/m3), but only 18% of the milling samples. In dermal samples, pyrene was detected above the LOD of 2.3 ng/cm2 in 61% of the paving samples (mean 5.7 ng/cm2), but only 27% of the milling samples. Exposure measurements for pavers were significantly higher than for millers for inhalation (P = 0.001), while differences were marginally significant for dermal exposure (P = 0.1). Additionally, there was a strong correlation between the inhalation and dermal exposure measurements obtained for paving workers (r = 0.66).
|
Table 2 presents the summary statistics for the pre-shift, post-shift and bedtime urinary 1-OHP data for each day of the study period. Figure 1 shows the mean urinary 1-OHP results for pavers and millers. At pre-shift on Monday morning (following a weekend of no exposure) pavers and millers were found to have the same mean 1-OHP level of 0.4 µg/g creatinine. Among paving workers there was a clear increase from pre-shift to post-shift during each day of the study period. There was also a clear increase in the pre-shift data during the work week, such that the average pre-shift results on day 4 were 3.5 times higher than the average pre-shift results on day 1. However, among the milling workers the urinary 1-OHP levels remained around the baseline value with very little change between pre-shift and post-shift and very little change during the work week.
|
|
Figure 2 shows the urinary 1-OHP results averaged by paving task. According to this graph, the screed men and rakers had the highest urinary 1-OHP results, followed by the paver operators and, finally, the roller operators. The results of a linear mixed effects model confirmed that urinary 1-OHP levels were in fact significantly different by paving task(P = 0.01).
|
Table 3 presents the results of three different linear mixed effects models evaluating the distributed lag effects of inhalation and/or dermal exposure, while controlling for age, BMI and smoking status. Model 1 evaluated the effect of inhalation exposure, Model 2 evaluated the effect of dermal exposure and Model 3 evaluated the effect of both inhalation and dermal exposure. To compare the effect of inhalation exposure and dermal exposure at the mean exposure level, the contribution of each exposure pathway was calculated by multiplying each effect estimate by the corresponding mean exposure for each lag period and exposure pathway. Ultimately, the total effect of exposure was calculated by summing the effect estimates for each exposure pathway, while the total contribution was calculated by summing the contribution terms for each exposure pathway.
|
The results of Model 1 indicate that urinary 1-OHP was increased due to inhalation exposure that occurred during the preceding 32 h, although the effect estimates were statistically significant for only the first 16 h. Similarly, the results of Model 2 indicate that urinary 1-OHP was increased due to dermal exposure that occurred during the preceding 48 h, although the effect estimates were statistically significant for only the first 32 h. Based on the results of these two separate models, the total contribution of dermal exposure (0.64) appears to be twice that of inhalation exposure (0.32).
The results of Model 3 indicate that urinary 1-OHP was increased due to inhalation exposure that occurred during the preceding 16 h, although the effect of inhalation was not statistically significant for either 8 h exposure period (Table 3). Urinary 1-OHP was increased due to dermal exposure that occurred during the preceding 56 h, although the effect estimates were statistically significant for only the first 32 h (Table 3). The largest effect of inhalation exposure was observed for the most recent exposure period (08 h) while the largest effect of dermal exposure was observed for the first lagged period (816 h). Based on the results of this model, the total contribution of dermal exposure (0.56) appears to be eight times that of inhalation exposure (0.07).
For dermal exposure, the results of Model 3 were very similar to the results of Model 2, with the only real differences being a slight reduction in the total effect of dermal exposure (from 0.41 to 0.37) and the fact that the largest effect shifted from the most recent dermal exposure (08 h) in Model 2 to the first lagged dermal exposure (816 h) in Model 3. Conversely, for inhalation exposure, the results of Model 3 were quite different from the results of Model 1. When inhalation exposure was included in the model without dermal exposure (Model 1) the results indicated that the total effect of inhalation (3.5) was significant, with a positive effect being observed during the preceding 32 h. However, when dermal exposure was added to the model (Model 3) the total effect of inhalation was reduced to 0.77 and was not significant, with a positive effect being observed during the preceding 16 h.
Figure 3 shows the contributions of inhalation exposure and dermal exposure at the mean inhalation and dermal exposure levels according to the results of Model 3. The error bars show the 95% confidence intervals for each of the estimates. Each lag period represents an 8 h period of exposure, where lag 0 represents the most recent exposure period (08 h), lag 1 represents the preceding 816 h, etc.
|
As mentioned, each of the models described above also evaluated the effect of age, BMI and smoking status. The results of Model 3 indicate that age had a positive effect on urinary 1-OHP (P = 0.08), BMI had a negative effect on urinary 1-OHP (P = 0.06) and smokers tended to have higher urinary 1-OHP levels as compared with non-smokers (P = 0.09). The effects of age, BMI and smoking were all marginally significant and were retained in the final models. It should be noted that consistent results were observed without adjusting the urinary 1-OHP concentrations for creatinine.
The effect of diet was initially included in the models and evaluated the effect of grilled food versus non-grilled food, as reported by questionnaire. The effect of diet was not significant (P = 0.5) and was excluded from the final models.
| DISCUSSION |
|---|
|
|
|---|
Although most studies evaluating PAH exposure have focused primarily on inhalation exposure, an increasing number have presented evidence suggesting that dermal exposure is actually the primary route of exposure to PAHs (Van Rooij et al., 1993a
75% of the total absorbed dose was attributable to dermal exposure. In creosote facility workers, Borak et al. (2002)The dermal data are presented in units of ng/cm2, which represents the amount of PAH per square centimeter of exposed skin (at the wrist) during each work shift. The wrist measurements are not intended to represent total body dermal exposure, since exposures of most other parts of the body are likely to be lower than of the wrist. However, the wrist samples do provide a useful tool for characterizing and comparing the intensity of daily dermal exposures.
Urinary 1-OHP by job
The paving workers had inhalation and dermal exposures that were significantly higher than the milling workers, while Fig. 1 shows that the urinary 1-OHP levels appear to be consistent with these differences in pyrene exposure.
At pre-shift on day 1, pavers and millers were found to have the same mean 1-OHP level (0.4 µg/g creatinine). This baseline level is probably due to non-occupational exposures such as smoking, diet and/or other environmental exposures. It is unlikely to be attributable to occupational exposure since these samples were always collected following a weekend away from work. Also, since we know that the pyrene exposures associated with each job are significantly different, we would not expect these 1-OHP levels to be identical for pavers and millers if they were actually job related.
While the millers showed little deviation from the baseline 1-OHP level during the work week, the pavers showed a clear increase from pre-shift to post-shift during each day of the study period. Among pavers, there was also a clear increase in the daily pre-shift data during the work week, such that the average pre-shift results on day 4 were 3.5 times higher than the average pre-shift results on day 1.
Urinary 1-OHP by task
Among paving workers, the pyrene data presented in Table 1 indicate that the task-based means ranged from 0.06 to 0.6 µg/m3 for inhalation exposure (roller operators < rakers < screed men < paving operators) while the task-based means ranged from 2.1 to 7.7 ng/cm2 for dermal exposure (roller operators < paving operators < rakers < screed men). In a separate analysis of these exposure data (McClean et al., 2004
), both the inhalation and dermal exposures were found to vary significantly by paving task.
Figure 2 shows the task-based urinary 1-OHP levels, which were also found to vary significantly by paving task (P = 0.01). The results of the linear mixed effects model indicated that the rank order of the task-based urinary 1-OHP levels was identical to the rank order of the task-based dermal exposures (roller operators < paving operators < rakers < screed men). A closer look at the paving operators is particularly interesting given that the paving operators had the highest inhalation exposures but the second lowest dermal exposures. Consistent with dermal exposure, the paving workers also had the second lowest urinary 1-OHP levels. This is evident in Fig. 2 and was confirmed by the results of the linear mixed effects model.
A work-shift increase in urinary 1-OHP levels is observed to varying extents for each paving task (Fig. 2). For screed men and rakers, the tasks with the highest dermal exposure, there is a clear increase in 1-OHP levels from pre-shift to post-shift during each day of the study period; however, the daily spike in 1-OHP level is less pronounced for the paver operators and not particularly evident for the roller operators. It is possible that the effect of non-occupational exposures is more graphically evident for the paver operators and roller operators.
Distributed lag models
The use of polynomial distributed lag models allowed the analysis of 1-OHP levels while evaluating the effect of inhalation and dermal exposure during the preceding 56 h, divided into seven 8 h time periods. Because the inhalation and dermal measurements were only obtained during the 8 h work shifts, this modeling approach required an assumption that inhalation and dermal exposure were 0 during all 8 h time periods that were spent away from work. This included the 8 h period between post-shift and bedtime, the 8 h period between bedtime and pre-shift (the next day) and all six 8 h time periods during the preceding weekend.
Figure 3 shows the mean contribution of each exposure pathway during each lag period. The results indicate that the preceding 16 h of inhalation exposure had a positive effect on urinary 1-OHP, although the effect estimates were not significant for either of the two 8 h lag periods. The fact that the largest effect estimate was observed for the most recent inhalation exposure suggests that the absorption, metabolism and excretion of inhaled pyrene was quite rapid. The results also indicate that the preceding 56 h of dermal exposure had a positive effect on urinary 1-OHP, although the effect estimates were only significant for the most recent exposure and for the first three lag periods (032 h). The fact that the largest effect estimate was observed for the first lag period (816 h), as well as the fact that the effect estimates were very similar during each of the first three exposure periods (024 h), suggests that the absorption of pyrene by dermal contact was slower than by inhalation. These results are consistent with Viau and Vyskocil (1995)
, who estimated that the excretion of urinary 1-OHP peaks between 10 and 15 h after dermal exposure to pyrene.
The results of Model 3 indicate that the impact of dermal exposure was eight times that of inhalation exposure, while controlling for the effect of age, BMI and smoking. In fact, the total effect of dermal exposure was significant while the total effect of inhalation exposure was not significant. Also, while the total effect of dermal exposure was very similar between models (Model 2 versus Model 3), the total effect of inhalation exposure was considerably different between models (Model 1 versus Model 3). Model 1 was constructed such that inhalation exposure was evaluated without controlling for the effect of dermal exposure, just as Model 2 evaluated dermal exposure without controlling for the effect of inhalation exposure. By including both inhalation and dermal exposure in the model together, Model 3 allowed us to evaluate the effect of each exposure pathway while controlling for the effect of the other.
A comparison of the models indicates that the total effect of inhalation exposure was considerably decreased (82% reduction) once dermal exposure was added to the model, while the total effect of dermal exposure was reduced by only 7% once inhalation exposure was added to the model. Since we know that inhalation exposure and dermal exposure are correlated, it appears that the effect of inhalation from Model 1 was largely due to confounding, since dermal exposure was not evaluated in the model. If dermal exposure had not been measured and included in this analysis, the total effect of inhalation would have been overestimated due to the correlation with dermal exposure. However, the results of Model 3 indicate that the dermal pathway is actually the primary route of exposure to PAHs among this group of workers.
In addition to evaluating the effect of inhalation and dermal exposure on urinary 1-OHP, each of the three models also found a marginally significant effect of age, BMI and smoking. The results of Model 3 indicate that BMI had a negative effect on urinary 1-OHP, suggesting that the excretion of urinary 1-OHP decreases as BMI increases. Since pyrene is very lipid-soluble, one possible explanation for this association is that increased BMI may result in increased storage of pyrene in fatty tissue such that less of the pyrene is metabolized to 1-OHP and excreted. This finding is consistent with Roggi et al. (1997)
, who also found that BMI significantly influences the levels of urinary 1-OHP.
| CONCLUSIONS |
|---|
|
|
|---|
This study represents one of the first applications of distributed lag models to biological monitoring data. The results of this urinary 1-OHP analysis suggest that among this group of paving workers the impact of dermal exposure was eight times the impact of inhalation exposure. Dermal exposure that occurred during the preceding 32 h had a statistically significant effect on urinary 1-OHP, while the effect of inhalation exposure was not significant. These findings are consistent with a growing body of literature suggesting that dermal contact is the primary route of occupational exposure to PAHs. Accordingly, an assessment of occupational exposure to PAHs that does not consider dermal exposure may considerably underestimate cumulative exposure. Similarly, control strategies aimed at reducing occupational exposure to asphalt-related PAHs should include an effort to reduce dermal exposure.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
We would like to thank Dr Doug Dockery for his contribution. The authors also wish to thank the National Cancer Institute for funding this project (grant no. 5R01 CA74413-03).
| REFERENCES |
|---|
|
|
|---|
ACGIH. (2002) 2002 TLVs and BEIs: threshold limit values for chemical substances and physical agents and biological exposure indices. Cincinnati, OH: American Conference of Governmental Industrial Hygienists.
Asphalt Institute. (1990) Report to OSHA and NIOSH: status of Asphalt Industry Steering Committee research program on the health effects of asphalt fumes and recommendation for a worker health standard. Lexington, KY: AsphaltInstitute.
Boffetta P, Jourenkova N, Gustavsson P. (1997) Cancer risk from occupational and environmental exposure to polycyclic aromatic hydrocarbons. Cancer Causes Control; 8: 44472.[CrossRef][Web of Science][Medline]
Boogaard PJ, van Sittert NJ. (1995) Urinary 1-hydroxypyrene as biomarker of exposure to polycyclic aromatic hydrocarbons in workers in petrochemical industries: baseline values and dermal uptake. Sci Total Environ; 163: 2039.[CrossRef][Medline]
Borak J, Sirianni G, Cohen H, Chemerynski S, Jongeneelen F. (2002) Biological versus ambient exposure monitoring of creosote facility workers. J Occup Environ Med; 44: 31019.[Web of Science][Medline]
Braga A, Saldiva P, Pereira L, Menezes J, Conceicao G, Lin C, Zanobetti A, Schwartz J, Dockery D. (2001) Health effects of air pollution exposure on children and adolescents in Sao Paulo, Brazil. Pediatr Pulmonol; 31: 10613.[CrossRef][Web of Science][Medline]
Chiazze L, Watkins DK, Amsel J. (1991) Asphalt and risk of cancer in man. Br J Ind Med; 48: 53842.[Web of Science][Medline]
Gamble JF, Nicolich MJ, Barone NJ, Vincent WJ. (1999) Exposureresponse of asphalt fumes with changes in pulmonary function and symptoms. Scand J Work Environ Health; 25: 186206.[Web of Science][Medline]
Gilbert RO. (1987) Statistical methods for environmental pollution monitoring. Van Nostrand Reinhold: 1778.
Hamlett A, Ryan L, Serrano-Trespalacios P, Wolfinger R. (2003) Mixed models for assessing correlation in the presence of replication. J Air Waste Manag Assoc; 53: 44250.[Web of Science][Medline]
Harvey A. (1990) The econometric analysis of time series. Cambridge, MA: MIT Press.
Jongeneelen FJ. (2001) Benchmark guideline for urinary 1-hydroxypyrene as biomarker of occupational exposure to polycyclic aromatic hydrocarbons. Ann Occup Hyg; 45: 313.
Jongeneelen FJ, Anzion RB, Henderson PT. (1987) Determination of hydroxylated polycyclic aromatic hydrocarbons in urine. J Chromatogr; 413: 22732.[Web of Science][Medline]
Jongeneelen FJ, Scheepers PT, Groenendijk A, Van Aerts LA, Anzion RB, Bos RP, Veenstra SJ. (1988) Airborne concentrations, skin contamination, and urinary metabolite excretion of polycyclic aromatic hydrocarbons among paving workers exposed to coal tar derived road tars. Am Ind Hyg Assoc J; 49: 6007.[Web of Science][Medline]
McClean M, Rinehart R, Ngo L, Eisen E, Herrick R. (2004) Inhalation and dermal exposure among asphalt paving workers. Ann Occup Hyg.
Moen BE, Nilsson R, Nordlinder R, Ovrebo S, Bleie K, Skorve AH, Hollund BE. (1996) Assessment of exposure to polycyclic aromatic hydrocarbons in engine rooms by measurement of urinary 1-hydroxypyrene. Occup Environ Med; 53: 6926.
NIOSH. (1998) Polycyclic aromatic hydrocarbons by HPLC: method 5506. In: NIOSH manual of analytical methods, 4th edn. Cincinnati, OH: US Department of Health and Human Services.
NIOSH. (2000) Health effects of occupational exposure to asphalt. Cincinnati, OH: US Department of Health and Human Services.
Partanen T, Boffetta P. (1994) Cancer risk in asphalt workers and roofers: review and meta-analysis of epidemiologic studies. Am J Ind Med; 26: 72140.[Web of Science][Medline]
Partanen TJ, Boffetta P, Heikkila PR, Frentzel-Beyme RR, Heederik D, Hours M, Jarvholm BG, Kauppinen TP, Kromhout H, Langard S et al. (1995) Cancer risk for European asphalt workers. Scand J Work Environ Health; 21: 2528.[Web of Science][Medline]
Pope C, Schwartz J. (1996) Time series for the analysis of pulmonary health data. Am J Respir Crit Care Med; 154: S22933.[Web of Science][Medline]
Quinlan R, Kowalczyk G, Gardiner K, Calvert IA, Hale K, Walton ST. (1995) Polycyclic aromatic hydrocarbon exposure in coal liquefaction workers: the value of urinary 1-hydroxypyrene excretion in the development of occupational hygiene control strategies. Ann Occup Hyg; 39: 32946.
Roggi C, Minoia C, Sciarra GF, Apostoli P, Maccarini L, Magnaghi S, Cenni A, Fonte A, Nidasio GF, Micoli G. (1997) Urinary 1-hydroxypyrene as a marker of exposure to pyrene: an epidemiological survey on a general population group. Sci Total Environ; 199: 24754.[CrossRef][Medline]
Schwartz J. (2000) The distributed lag between air pollution and daily deaths. Epidemiology; 11: 3206.[CrossRef][Web of Science][Medline]
Tolos WP, Lowry LK, MacKenzie BA. (1989) 1-Pyrenol in urine: a biological monitoring method to assess exposure to polycyclic aromatic hydrocarbons containing pyrene. In: Eleventh symposium on polycyclic aromatic hydrocarbons. New York: Gordon & Breach.
Van Rooij JG, Bodelier-Bade MM, Jongeneelen FJ. (1993a) Estimation of individual dermal and respiratory uptake of polycyclic aromatic hydrocarbons in 12 coke oven workers. Br J Ind Med; 50: 62332.[Web of Science][Medline]
Van Rooij JG, van Lieshout EM, Bodelier-Bade MM, Jongeneelen FJ. (1993b) Effect of the reduction of skin contamination on the internal dose of creosote workers exposed to PAH. Scand J Work Environ Health; 19: 2007.[Web of Science][Medline]
Van Rooij JG, De Roos JH, Bodelier-Bade MM, Jongeneelen FJ. (1993c) Absorption of polycyclic aromatic hydrocarbons through human skin: differences between anatomical sites and individuals. J Toxicol Environ Health; 38: 35568.[Web of Science][Medline]
Viau C, Vyskocil A. (1995) Patterns of 1-hydroxypyrene excretion in volunteers exposed to pyrene by the dermal route. Sci Total Environ; 163: 18790.[CrossRef][Medline]
This article has been cited by other articles:
![]() |
J. R. Sobus, M. D. Mcclean, R. F. Herrick, S. Waidyanatha, F. Onyemauwa, L. L. Kupper, and S. M. Rappaport Investigation of PAH Biomarkers in the Urine of Workers Exposed to Hot Asphalt Ann. Hyg., August 1, 2009; 53(6): 551 - 560. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. Sobus, M. D. McClean, R. F. Herrick, S. Waidyanatha, L. A. Nylander-French, L. L. Kupper, and S. M. Rappaport Comparing Urinary Biomarkers of Airborne and Dermal Exposure to Polycyclic Aromatic Compounds in Asphalt-Exposed Workers Ann. Hyg., August 1, 2009; 53(6): 561 - 571. [Abstract] [Full Text] [PDF] |
||||
![]() |
J R Sobus, S Waidyanatha, M D McClean, R F Herrick, T J Smith, E Garshick, F Laden, J E Hart, Y Zheng, and S M Rappaport Urinary naphthalene and phenanthrene as biomarkers of occupational exposure to polycyclic aromatic hydrocarbons Occup. Environ. Med., February 1, 2009; 66(2): 99 - 104. [Abstract] [Full Text] [PDF] |
||||
![]() |
C-H Pan, C-C Chan, Y-L Huang, and K-Y Wu Urinary 1-hydroxypyrene and malondialdehyde in male workers in Chinese restaurants Occup. Environ. Med., November 1, 2008; 65(11): 732 - 735. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Pesch, M. Kappler, K. Straif, B. Marczynski, R. Preuss, B. Rossbach, H.-P. Rihs, T. Weiss, S. Rabstein, C. Pierl, et al. Dose-Response Modeling of Occupational Exposure to Polycyclic Aromatic Hydrocarbons with Biomarkers of Exposure and Effect Cancer Epidemiol. Biomarkers Prev., September 1, 2007; 16(9): 1863 - 1873. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. KROMHOUT and I. BURSTYN Vapours and Aerosols of Bitumen: Exposure Data Obtained by the German Bitumen Forum Ann. Hyg., June 1, 2007; 51(4): 423 - 425. [Full Text] [PDF] |
||||
![]() |
M. D. McCLEAN, J. K. WIENCKE, K. T. KELSEY, A. VARKONYI, L. NGO, E. A. EISEN, and R. F. HERRICK DNA Adducts among Asphalt Paving Workers Ann. Hyg., January 1, 2007; 51(1): 27 - 34. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. OGDEN Annals of Occupational Hygiene at Volume 50: Many Achievements, a Few Mistakes, and an Interesting Future Ann. Hyg., November 1, 2006; 50(8): 751 - 764. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. RUHL and U. MUSANKE The German Bitumen Forum--Cooperation in Partnership Ann. Hyg., July 1, 2006; 50(5): 441 - 444. [Full Text] [PDF] |
||||
![]() |
J. UNWIN, J. COCKER, E. SCOBBIE, and H. CHAMBERS An Assessment of Occupational Exposure to Polycyclic Aromatic Hydrocarbons in the UK Ann. Hyg., June 1, 2006; 50(4): 395 - 403. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Marczynski, M. Raulf-Heimsoth, R. Preuss, M. Kappler, K. Schott, B. Pesch, G. Zoubek, J.-U. Hahn, T. Mensing, J. Angerer, et al. Assessment of DNA Damage in WBCs of Workers Occupationally Exposed to Fumes and Aerosols of Bitumen. Cancer Epidemiol. Biomarkers Prev., April 1, 2006; 15(4): 645 - 651. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. CAVALLO, C. L. URSINI, P. BAVAZZANO, C. CASSINELLI, A. FRATTINI, B. PERNICONI, A. DI FRANCESCO, A. CIERVO, B. RONDINONE, and S. IAVICOLI Sister Chromatid Exchange and Oxidative DNA Damage in Paving Workers exposed to PAHs Ann. Hyg., April 1, 2006; 50(3): 211 - 218. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. J. BOOGAARD and J. H. URBANUS Assessment of Asphalt Workers' Dermal Exposure Ann. Hyg., January 1, 2005; 49(1): 93 - 93. [Full Text] [PDF] |
||||
![]() |
M. D. MCCLEAN, K. T. KELSEY, and R. F. HERRICK Reply Ann. Hyg., January 1, 2005; 49(1): 93 - 94. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||













