Lung Cancer among Industrial Sand Workers Exposed to Crystalline Silica
Kyle Steenland and
Wayne Sanderson
From the National Institute for Occupational Safety and Health, Cincinnati, OH.
 |
ABSTRACT
|
---|
In 1997, the International Agency for Research on Cancer determined that crystalline silica was a human carcinogen but noted inconsistencies in the epidemiology. There are few exposure-response analyses. The authors examined lung cancer mortality among 4,626 industrial sand workers, estimating exposure via a job-exposure matrix based on 4,269 industrial hygiene samples collected in 19741995. The average length of employment was 9 years, and estimated average exposure was 0.05 mg/m3 (the National Institute of Occupational Safety and Health Recommended Exposure Limit). Results confirmed excess mortality from silicosis/pneumoconioses (standardized mortality ratio = 18.2, 95% confidence interval: 10.6, 29.1; 17 deaths). The lung cancer standardized mortality ratio was 1.60 (95% confidence interval: 1.31, 1.93; 109 deaths). Limited data suggested that smoking might account for 1020% of the lung cancer excess. Exposure-response analyses by quartile of cumulative exposure (15-year lag) yielded standardized rate ratios of 1.00, 0.78, 1.51, and 1.57 (p for trend = 0.07). Nested case-control analyses after exclusion of short-term workers, who had high overall morality, yielded odds ratios by quartile of cumulative exposure (15-year lag) of 1.00, 1.35, 1.63, and 2.00 (p for trend = 0.08) and odds ratios by quartile of average exposure of 1.00, 0.92, 1.44, and 2.26 (p = 0.005). These data lend support to the labeling by the International Agency for Research on Cancer of silica as a human carcinogen. There are approximately 2 million US workers exposed to silica; 100,000 are exposed to more than 0.1 mg/m3.
lung neoplasms; silicon dioxide; silicosis
Abbreviations:
IARC, International Agency for Research on Cancer; NIOSH, National Institute for Occupational Safety and Health; NDI, National Death Index; SMR, standardized mortality ratio; SRR, standardized rate ratio
 |
INTRODUCTION
|
---|
In the 1980s, there were an estimated 1.7 million US workers exposed to crystalline silica outside of the mining industry (1
), of whom approximately 100,000 are exposed to a level two or more times greater than the National Institute of Occupational Safety and Health (NIOSH) Recommended Exposure Limit of 0.05 mg/m3 (2
). Exposure occurs in foundries, stonework, sandblasting, and potteries. Besides its well-known relation to silicosis, silica exposure in recent years has been associated with lung cancer, and in 1997, the International Agency for Research on Cancer (IARC) determined that inhaled crystalline silica from occupational sources is a definite (group 1) human carcinogen (3
). However, IARC qualified its decision by noting that the evidence was not entirely consistent and that different forms of silica might have different potentials for causing lung cancer. Furthermore, there are few studies with adequate data to analyze exposure-response trends, which can provide the strongest evidence for causality, and some studies with exposure-response data were conducted among miners and may have had confounding exposures from other lung carcinogens such as radon.
We have investigated lung cancer in a cohort mortality study of workers in the industrial sand industry and have estimated past exposures, thereby permitting an exposure-response analysis. Many industrial sand plants processed the sand into silica flour, which is composed of fine particles of crystalline silica and which has been associated with high levels of silicosis in the past (4
, 5
). The industrial sand industry has not had confounding exposures to other known occupational lung carcinogens such as radon and arsenic, confounders that can affect studies of underground mines exposed to silica. Kidney disease and arthritis, other diseases recently related to silica exposure, are discussed in a separate publication (manuscript in preparation).
 |
MATERIALS AND METHODS
|
---|
Background on the industry
In 1987, the industrial sand industry employed approximately 2,600 workers in 60 plants in the United States (6
). Industrial silica sand is obtained from a variety of sources, from a loose, unconsolidated granular state to hard, highly compacted rocks; the ore form of quartz determines how it will be mined. Hard-rock mining is usually done in an open pit quarry where holes are drilled into the rock layers, filled with explosives, and detonated to break the rock into movable pieces. Uncompacted sand is collected by using dredges, hydraulic pumps, scrapers, or clamshells. Trucks are often used to transport the sand to the plant for processing, but pipelines and belt conveyors are also used. At the plant, rocks are crushed to obtain progressively smaller sizes; if the sand is already granular, crushing is usually not necessary. The crushed and granular sand is screened and sized. The screening and sizing operation may be either a wet process or a dry one; dry screening typically has higher dust exposures. Often the sand is milled in rotating ball mills to reduce the sand to a fine powder (silica flour). The final sand products, varying in particle size and quartz content, are bagged or bulk-loaded for shipment. Silica sand has a variety of industrial uses. It is used in sandblasting, fine polishing, glassmaking, and rubber manufacture (where silica powder is used as a dry lubricant); as a filler in plastics and paints; and as a carrier in cosmetics (7
). The industrial sand industry has little potential for exposures to other known or suspected lung carcinogens. There is some potential of exposure to diesel exhaust from heavy equipment for quarry workers (11 percent of all workers); however, this exposure would be expected to be minimal because it is outside and not in a confined space.
Cohort definition
Our study covers 18 plants in 11 different states, chosen because they were part of a trade association and had adequate records. It is based on 5,086 personnel records of former and current workers collected in 19871988 at these plants. All of the plants had discarded the records of workers who terminated prior to some specific year, typically in the 1950s or 1960s. Therefore, person-time began at each of these plants no earlier than that specific year. Because personnel records were not sufficiently accurate to determine whether persons who had worked for only a few days had, in fact, been exposed to silica, all workers were required to have worked at least 1 week. A total of 249 workers (5 percent) had worked for less than 1 week. An additional 156 workers (3 percent) had inadequate data for birth date or work history dates to be included, leaving 4,681 subjects who had worked for a week or more available for analysis. Because US referent rates for certain diseases of interest (e.g., silicosis, multiple-cause analyses of other causes) were available only from 1960 onward, we began mortality follow-up in 1960. This decision led to the elimination of another 55 people (13 deaths) who could not be traced through 1960 (i.e., they had died or were lost to follow-up prior to 1960), leaving a final cohort of 4,626. Because only a small number of people were lost by beginning follow-up in 1960, results were virtually unchanged from those without this restriction on the beginning of follow-up.
Analyses for exposure-response trends were conducted on a smaller subcohort with good work history data, i.e., complete data on dates of employment and the type of job performed while employed. Exposure-response analyses were conducted by using a job-exposure matrix in which estimated exposure levels were assigned to workers on the basis of the jobs they held, the plants in which they worked, and the calendar period in which they worked, as discussed below. Information on where 13 percent of the employees worked was lacking; while employment dates were available, personnel records did not provide information on which jobs these workers held (most of these employees were at one plant that had poor records). These were eliminated for exposure-response analyses, leaving 4,027 available for those analyses. However, people who worked in an unknown job for part of their work history were included in exposure-response analyses; the unknown jobs were assigned the mean exposure level for a plantwide worker with a job title such as "general laborer" or "maintenance worker" in the same plant at the same time period. Only 215 of the 4,027 workers (5 percent) in the exposure-response analyses had an unknown job at some point during their work history, and these unknown jobs represented only 2 percent of all jobs held by these workers (workers had a mean of three jobs during their employment). A further subset of the 4,027 workers in the exposure-response analyses, those with at least 6 months employment (n = 3,361), were also analyzed for exposure-response trends.
Exposure data and job-exposure matrix
Quantitative estimates of exposure to respirable silica were estimated for each worker over time by plant and job via a job-exposure matrix, created by using existing exposure data by job and plant. A thorough description of this procedure is available elsewhere (8
). Briefly, 4,269 personal samples for respirable silica were available for 19741996, either from the Mine Safety Health Administration (available at all 18 plants) or from company data (available at seven of 18 plants) (modeling indicated no significant differences between Mine Safety Health Administration and company data). These samples covered a wide variety of jobs. These data were then restricted to 19741988 (1988 was the end of data collection for our study, so work histories ended in that year) and modeled by using general linear models. Predicted values from the model were used to create a job-exposure matrix for the period 19741988. In this matrix, we estimated exposure levels for four categories of plants (low, medium, medium-high, and high), three time periods (19741979, 19801984, and 19851988), and 10 job categories. The time periods were chosen by observing when changes in levels appeared to have taken place, and the 10 job categories were based on the main areas where subjects worked and have been used in other surveys (9
).
In addition to the silica measurements after 1974, an exposure assessment study using midget impingers was conducted in 19 silica sand plants in 1946. In this study, dust particles smaller than 5 µm in diameter were counted optically and reported in concentrations of million particles per cubic foot. A summary report of this study presented the results of the sampling as mean dust concentrations (n = 125) by job and plant (9
), within the same job categories used in the 19741988 data described above. These mean impinger-dust concentrations were converted to respirable-mass silica concentrations in µg/m3 by multiplying them by the percent of silica (quartz) found by job in the respirable dust samples from 1974 to 1988 and a conversion factor of 0.1 (10
12
). The converted, job-specific medians of the same 10 job categories considered in the 19741988 data were used to estimate job-specific exposures in 1946. The silica respirable concentrations were then extrapolated linearly between 1946 and 1974 for each job category. Exposures prior to 1947 were considered to be constant.
Follow-up and analytic methods
Follow-up was conducted via the National Death Index (NDI) for the years 19791996, the Social Security death tapes, and the Internal Revenue Service. Cause-of-death was obtained from NDI or from death certificates obtained from the states for deaths that occurred prior to 1979 (n = 333) (not covered by NDI). Social Security death tapes would be expected to cover about three quarters of the deaths occurring in working populations during the years 19601978, the follow-up years not covered by NDI in our study (13
); about 6 percent of our cohort was last known alive and was not followed through 1979 and, hence, would be susceptible to underascertainment of deaths during this period. This 6 percent was relatively young, and the number of expected deaths missed among this group would be minimalprobably less than 10, or less than 1 percent of all deaths observed.
We used standard life-table analyses for person-time data to analyze the cause-specific rates for the exposed population versus the US population (14
). Follow-up began on January 1, 1960, on the date of first exposure, or on the date when records were first available, whichever was later. Follow-up continued until the date of death for deceased subjects, until December 31, 1996 (the end of the NDI search), or until a subject was last known to be alive if that date was earlier than January 1, 1979. Subjects known to be alive after January 1, 1979 and not found in the NDI to be deceased were assumed to be alive as of December 31, 1996. Life-table analyses for cancers considered unlagged data and data lagged 15 years, a typical minimum induction period for epithelial tumors. In addition to underlying-cause analyses, we also conducted multiple-cause analyses in which cause-specific mortality events were defined as any mention on the death certificate. Such analyses are particularly important for diseases that may be prevalent at death but that are not the underlying cause (e.g., tuberculosis, silicosis) (15
).
The data (person-time and observed deaths) were divided into quartiles of cumulative exposure (for each worker,
(level of exposure x duration of exposure), over all i jobs held) for the purposes of categorical exposure-response analyses. Quartile cutpoints were chose before analysis by using the distribution of cumulative exposure for all decedents, motivated by the desire to obtain approximately equally precise estimates of effect in each quartile. Effects were measured either by indirectly standardized mortality ratios (SMRs) by using the US population as referent or by directly standardized rate ratios (SRRs) using the lowest quartile of the exposed as the referent. Tests for linear trends in SRRs or SMRs with exposure were calculated via methods outlined by Rothman (16
) or Breslow and Day (17
), respectively. In these tests for trends, the midpoints of the exposure categories were used as scores; for the last category, we used the upper category boundary plus 50 percent of the upper category boundary.
More detailed internal exposure-response analyses were conducted via a nested case-control design analyzed by conditional logistic regression. This allowed for analyses by average exposure (cumulative exposure/duration of exposure) and peak exposure in addition to cumulative exposure. Different lag periods could be considered, as well as different transformation of exposure (e.g., log exposure). In these analyses, 100 controls were chosen randomly per case from among all those who had survived to the same age as the case or later. Controls were matched to cases on race, sex, and date of birth within 5 years. An SAS procedure (PHREG) (18
) was used to perform the analysis. Quartiles of exposure in these analyses were chosen on the basis of distribution of exposure among noncases. Tests for exposure-response trend were done via use of a continuous variable for exposure in the model.
Lifetime excess risk of lung cancer was calculated, taking into account competing risk of death from other causes (19
). Age-specific male lung cancer mortality rates were taken from the Surveillance, Epidemiology, and End Results data for 19871991 (20
); we used US 1991 age-specific, all-cause mortality rates.
Data on smoking
Limited data on smoking were available for the cohort. Cross-sectional data on smoking were available from company records for 404 men in the cohort. These data were collected during the years 19781989 at four of the largest plants in the study where NIOSH gathered more extensive records pursuant to an exposure survey. We focused our analysis of these data on the 358 men aged 2564 years for comparison with age-specific US data on smoking. We estimated the impact of smoking differences on lung cancer rates for these men versus the US population, which was our comparison population in our cohort mortality analyses. Using techniques for indirect adjustment outlined by Axelson and Steenland (21
), we estimated the lung cancer rate ratio between the cohort and the national population due to smoking differences, under the assumption that the men for whom data on smoking were collected at the four plants were representative of the entire cohort. In these calculations, we assumed that the lung cancer rate ratios for heavy smokers (more than one pack a day), light smokers (one pack a day or less), and former smokers were 20, 10, and 5, respectively. Data for the United States were obtained from the 1987 National Health Survey (22
). The expected lung cancer rate ratio between the cohort and the US population due to smoking differences alone is Icohort/IUS, where the lung cancer rates for the cohort (Icohort) and for the United States (IUS) are each a weighted average of rates for never smokers (rate = I0), former smokers (rate = 5I0), light smokers (rate = 10I0), and heavy smokers (rate = 20I0), as follows:
 | (1) |
 |
RESULTS
|
---|
Table 1 gives descriptive data on the cohort. The mean length of follow-up was 24 years, and the mean year of first employment was 1967. The mean length of employment was 9 years, with a wide range: 20 percent were employed for less than 6 months, while 31 percent were employed for more than 10 years. Twenty-four percent of the cohort had died, and the cause of death was available for 95 percent of the deceased. If missing causes of death were distributed equally among different causes, then one would expect that cause-specific mortality rates in our cohort would be underestimated by approximately 5 percent. There were 101,177 person-years at risk for the entire cohort (n = 4,021).
Table 2 gives exposure data for the 10 job categories used in the analysis for the period 19741996. It is apparent from the table that high exposures were not uncommon in these plants. The adjusted geometric means from the model for plants with high, medium-high, medium, and low exposure in the period 19741996 were 41.6, 24.8, 15.4, and 11.1 µg/m3, respectively. Exposures have dropped considerably over time. The geometric mean exposure in 19741979 (based on 1,278 samples) was 51 µg/m3, dropping to 11.6 µg/m3 for the period 19851988 (based on 680 samples). The median exposure in 19461947 was 78 µg/m3 (based on 125 sample means).
Table 2 also shows the distribution of job categories for workers in their last job for the 4,021 workers in the exposure-response analysis, who had an mean of 3.8 jobs. The most common job category was "other," which included "laborers" (20 percent of cohort), general maintenance (7 percent), and mechanics (6 percent). The second most common category was "administrative," which included managers and clerical staff (8 percent), testers (2 percent), and guards (2 percent). Only 0.6 percent of the cohort had undefined or missing job information for their last job.
Table 3 gives data for selected causes of death for the cohort versus the US population, for both underlying and multiple causes. As expected, there are large excesses for silicosis, unspecified pneumoconiosis (probably silicosis), and tuberculosis (often increased among populations exposed to silica). Lung cancer shows a 60 percent excess. There is a general increase of 2030 percent in mortality for this cohort from all causes combined, cancer, and heart disease. It is likely that lifestyle characteristics such as smoking account for some of this increase, although exposure to silica is likely to play a role for some specific causes, as discussed below.
View this table:
[in this window]
[in a new window]
|
TABLE 3. Standardized mortality ratios and 95% confidence intervals for selected causes of death, underlying and multiple-cause analyses in 18 industrial sand plants in the United States
|
|
Table 4 presents exposure-response analyses for selected causes for the subcohort with good work histories (eliminating 13 percent of the original cohort who worked in unknown jobs only). SMRs for lung cancer show elevations in both the first and fourth quartiles of cumulative exposure. Internal SRR analyses for lung cancer using the lowest quartile as referent show a positive trend in SRRs (p = 0.07) for cumulative exposure lagged 15 years. Silicosis and other unspecified pneumoconioses (probably silicosis in this cohort) as underlying causes show a positive trend in SRRs with cumulative exposure (p < 0.00001). Multiple-cause analyses of respiratory tuberculosis, emphysema, and other nonmalignant respiratory disease (all nonmalignant disease excluding emphysema, bronchitis, asthma, and pneumonia; 20 percent of deaths in this category had silicosis mentioned on the death certificates) also show positive trends in SRRs with increasing exposure. These categories, especially tuberculosis and nonspecific chronic obstructive pulmonary disease, would be expected to be related to increasing silica exposure.
View this table:
[in this window]
[in a new window]
|
TABLE 4. Exposure-response analyses: standardized mortality ratios (observed deaths) and standardized rate ratios (observed deaths) by quartile of cumulative exposure to respirable silica for selected causes (n = 4,027)* in 18 industrial sand plants in the United States
|
|
Analysis of the lung cancer mortality data by duration of employment did not show any consistent trend, but did show elevated lung cancer for those with less than 6 months exposure and for those with 1020 years and 20 or more years exposure. The lung cancer SMRs for less than 6 months; 6 months to 1 year; and 12, 25, 510, 1020, and 20 or more years were 2.38 (95 percent CI: 1.17, 4.22), 1.67 (95 percent CI: 0.54, 3.90), 2.27 (95 percent CI: 0.91, 4.68), 1.11 (95 percent CI: 0.44, 2.29), 0.84 (95 percent CI: 0.27, 1.96), 1.67 (95 percent CI: 0.99, 2.64), and 1.54 (95 percent CI: 1.05, 2.17), respectively (based on 11, 5, 7, 7, 5, 18, and 32 deaths, respectively).
Analysis of mortality for all causes combined showed an similar increased mortality for short-term workers with less than 6 months duration of employment. The SMRs for all causes for under 6 months, 6 months to 1 year, and 12, 25, 510, 1020, and 20 or more years employment were 1.47, 1.03, 1.10, 1.17, 1.00, 1.35, and 1.30, respectively, based on 98, 48, 50, 103, 80, 179, and 302 deaths, respectively.
Short-term workers are known to have high excess mortality, often for reasons unrelated to occupational exposures (23
, 24
). We therefore reanalyzed the data restricted to those with 6 or more months employment (n = 3,361). SMR results for these analyses are shown in table 5. While high SMRs for the lowest quartile are now somewhat reduced, they are still elevated. The U-shaped curve of SMRs from table 4 is still apparent, although less marked.
View this table:
[in this window]
[in a new window]
|
TABLE 5. Exposure-response analyses: lung cancer standardized mortality ratios (observed deaths) by quartile of cumulative exposure to respirable silica in a cohort with 6 or more months employment (n = 3,361)* at 18 industrial sand plants in the United States
|
|
Results of nested case-control analyses of lung cancer for those with more than 6 months employment are shown in table 6. The results for cumulative exposure generally parallel those in table 5. Average exposure (defined as cumulative exposure divided by duration of exposure) showed somewhat more pronounced trends than did cumulative exposure. There was a twofold excess risk for those in the highest quartile with either exposure. Log transformation of exposure variables did not improve the fit of the model. Neither duration of exposure nor lagged duration of exposure showed a significant relation to lung cancer risk. Inclusion of a dichotomous variable for the state in which the plant was located, dichotomized based on state lung cancer rates, had little effect on results.
View this table:
[in this window]
[in a new window]
|
TABLE 6. Odds ratios and 95% confidence intervals for a nested case-control study of lung cancer among workers with more than 6 months employment in 18 industrial sand plants in the United States
|
|
Review of the death certificates of lung cancer cases found mention of silicosis on only two of 109 death certificates, less than the approximately four instances of joint occurrence expected (43 of 1,073 deaths had silicosis mentioned on the death certificates). This analysis is limited because mortality data underestimate silicosis and because of the small numbers involved.
Assuming a mean exposure level of 0.05 mg/m3, exposure at the NIOSH Recommended Exposure Limit from ages 2065 years, and results for the best model in the nested case-control study (average exposure), the lifetime excess lung cancer risk through age 75 for a male would be 1.3 percent (95 percent CI: 0.4, 2.4). With the model for cumulative exposure lagged 15 years, the estimated lifetime excess lung cancer risk was similar (1.4 percent, 95 percent CI: 0.3, 3.8), although confidence intervals were wider. For comparison, the background lung cancer lifetime risk by age 75 years for a US male not exposed to silica is 5.3 percent (smokers and nonsmokers combined). The Occupational Safety and Health Administration generally seeks to regulate levels so that lifetime excess risk does not excess 1 in 1,000.
Data on smoking collected on 346 men aged 2564 years employed at four plants during the years 19781989 are shown in table 7. The most common year of data collection (the mode) was 1987. Table 7 also presents the data for the cohort collected in 1987 and the comparable data for US males in 1987. Use of calendar time- and age-specific data to compare smoking habits between the cohort and the US population is important for a time when smoking habits were changing, as they were in the United States during the 1980s. It is clear that the cohort, especially the younger workers, smoked somewhat more than did the comparable US population. However, differences were not dramatic. We used the 1987 data for the cohort and the United States and applied an indirect adjustment to estimate the effects of smoking differences between the cohort and the United States on expected lung cancer rate ratios. Under the assumption that the smoking differences seen in table 7 are representative of the entire cohort, one would expect a lung cancer rate ratio for the cohort versus the United States of 1.18 based on the data for men aged 2544 years and a lung cancer rate ratio of 1.09 based on the data for men aged 4564 years.
View this table:
[in this window]
[in a new window]
|
TABLE 7. Data on smoking collected in a cohort of 346 men aged 2564 years at four industrial sand plants in the United States during 19781989
|
|
We also considered the effect of smoking differences within our cohort to see whether they might account for our positive internal exposure-response trends. Categorical analyses of workers aged 2564 years for smoking status (never, current, former) by four cumulative exposure categories (those used in table 6), controlling for age category, showed little association between smoking status and cumulative exposure (p = 0.25). A regression of cumulative exposure (continuous variable) on smoking status controlling for age found that neither current (p = 0.37) nor former smokers (p = 0.28) had significantly higher cumulative dose than did never smokers. We also divided our sample of workers with data on smoking into two categories of cumulative exposure (<0.59 and
0.59 mg/m3-years) and applied an Axelson-type adjustment. We found that the high-exposure group would be expected to have a lung cancer rate ratio of 1.04 compared with the low-exposure group for those aged 45 years and above because of smoking differences. The corresponding rate ratio for those under age 45 years was 1.10.
 |
DISCUSSION
|
---|
This cohort is valuable for studying silica exposure and lung cancer because 1) it was unlikely to have had significant exposure to other occupational lung carcinogens, which can occur with silica-exposed workers (e.g., miners, foundry workers), 2) it is reasonably large (109 lung cancer deaths), and 3) it has industrial hygiene sampling data over time, which enabled us to estimate historical levels of silica exposure by job, plant, and calendar time. These exposure estimates enabled us to conduct exposure-response analyses.
Such analyses are less subject to confounding by smoking than are simple comparisons of lung cancer rates in exposed workers versus the general population because workers with a high exposure may be compared with those with a low exposure who are more likely to share similar smoking and other lifestyle habits. Exposure-response trends, or the lack of them, also provide valuable evidence with which to draw conclusions about whether the agent in question truly causes the disease in question.
Weaknesses in our study include 1) the availability of only limited data on smoking, restricted to about 10 percent of the cohort, gathered in the 1980s, and 2) the lack of historic industrial hygiene data for the period 19471974, which required use of linear extrapolation between these time points.
We had no gold standard to confirm the validity of our estimates of historical exposure. However, partial confirmation was provided by the significant positive exposure-response trends observed for mortality from diseases known to be caused solely or partly by silica, i.e., silicosis, unspecified pneumoconiosis, tuberculosis, and other nonmalignant respiratory disease.
We found a 60 percent excess of lung cancer overall in the exposed cohort compared with the US population (SMR = 1.60, 95 percent CI: 1.31, 1.93). The cohort experienced high mortality for all causes combined (SMR = 1.23, 95 percent CI: 1.16, 1.31), in part due to elevations in causes related to silica exposure, but probably also partly due to nonoccupational reasons. We found increased smoking in our sample of workers from four plants in the 1980s. However, the amount of increased smoking, compared with the smoking habits of the US population, would have been expected to cause an excess lung cancer mortality on the order of 1020 percent. Lung cancer rate ratios on the order of 1.101.20 due to smoking differences between workers and the general population are what one would expect based on both empirical data (25
, 26
) and theoretical data (21
). Such increased rate ratios due to confounding by smoking can only partially explain the observed 60 percent excess lung cancer rate found in the workers exposed to silica. Similarly, the relatively small smoking differences we found between high- and low-exposure groups, assuming that they are representative of the entire cohort, would be unlikely to explain much of our positive exposure-response trends.
Increased lung cancer risk due to increased smoking could also result from effect modification if the lung cancer rate due to smoking and silica together was markedly extremely high compared with that for nonsmokers exposed to silica (21
). There are only limited data to date on this point, with inference limited by the number of nonsmoker lung cancers as usual; however, what data do exist do not suggest effect modification by smoking (27
).
We also noted that there was an increased all-cause mortality for short-term workers with less than 6 months of employment, which may have reflected nonoccupational lifestyle factors. Elimination of these persons (27 percent of the cohort) reduced, but did not eliminate, the overall lung cancer excess (SMR for lung cancer = 1.49, 95 percent CI: 1.17, 1.87).
Positive exposure-response trends were found for lung cancer with cumulative silica exposure lagged 15 years and with average silica exposure. Most of the literature on silica has focused on cumulative exposure, but there has been one report in the literature of workers exposed to silica in which average exposure, but not cumulative exposure, predicted lung cancer risk (28
). In our data, the correlation between average and cumulative silica exposures was 0.45. In contrast to cumulative and average exposures, we found no significant positive trend for lung cancer with duration of exposure, which is a component of cumulative exposure but ignores the level or intensity of exposure.
Exposure-response trends were not always monotonic; workers with the lowest exposures included short-term workers who had relatively high lung cancer rates and high overall mortality rates. Trends for cumulative exposure were stronger and more consistent when those with less than 6 months employment were eliminated.
In conclusion, our findings tend to support the 1997 judgment by IARC that crystalline silica is a lung carcinogen.
 |
ACKNOWLEDGMENTS
|
---|
The authors thank Dr. Harlan Amandus for his advice throughout this project and Drs. Harvey Checkoway, Eva Hnizdo, and Faye Rice for their comments on earlier drafts.
 |
NOTES
|
---|
Reprint requests to Dr. Kyle Steenland, National Institute for Occupational Safety and Health, R13, 4676 Columbia Parkway, Cincinnati, OH 45226 (e-mail:kns1{at}cdc.gov).
 |
REFERENCES
|
---|
-
National Institute for Occupational Safety and Health. Work-related diseases surveillance report. Cincinnati, OH: National Institute for Occupational Safety and Health, 1991. (NIOSH publication no. 91113).
-
Linch K, Miller W, Althouse R, et al. Surveillance of respirable crystalline silica dust using OSHA compliance data (19791995). Am J Ind Med 1998;34:54758.[ISI][Medline]
-
International Agency for Research on Cancer. Silica, some silicates, coal dust, and para-aramid fibrils. Monograph 68. Evaluation of carcinogenic risks to humans. Lyon, France: International Agency for Research on Cancer, 1997.
-
Banks D, Morring K, Boehlecke B. Silicosis in the 1980s. J Am Ind Hyg Assoc 1981:42:779.
-
Banks D, Morring K, Boehlecke B, et al. Silicosis in silica flour workers, Am Rev Respir Dis 1981;124:44550.[ISI][Medline]
-
Amandus H. Protocol for a study of the health status of industrial sand workers. Morgantown, WV: Division of Respiratory Disease Studies, National Institute for Occupational Safety and Health, 1987.
-
Davis G. Silica-section V agents causing interstitial disease. In: Harber P, Schenker M, Balmes J, eds. Occupational and environmental respiratory disease. St. Louis, MO: J. Mosby, 1996.
-
Sanderson W, Steenland K, Deddens J. Historical respirable quartz exposures of industrial sand workers, 19461996. Am J Ind Med 200;38:38998.[ISI][Medline]
-
Hatch T, Holden F, Haines G. Investigations of dust exposure and controls in industrial sand production, Report to the National Industrial Sand Association. Washington, DC: National Industrial Sand Association, 1947.
-
Rice C, Harris R, Lumsden J, et al. Reconstruction of silica exposure in the North Carolina dusty trades. J Am Ind Hyg Assoc 1984;45:68996.
-
Ayer H, Dement J, Busch K, et al. A monumental studyreconstruction of a 1920 granite shed. J Am Ind Hyg Assoc 1973:34:20611.
-
Sheehy J, McJilton C. Development of a model to aid in reconstruction of historical silica dust exposures in the taconite industry. J Am Ind Hyg Assoc 1987;48:91418.
-
Schnorr T, Steenland K. Identifying deaths before 1979 using the SSA Death Master File. Epidemiology 1997;8:3213.[ISI][Medline]
-
Steenland K, Beaumont J, Spaeth S, et al. New developments in the life table analysis system of the National Institute for Occupational Safety and Health. J Occup Med 1990;32:10918.[ISI][Medline]
-
Steenland K, Nowlin S, Ryan B, et al. Use of multiple-cause mortality data in epidemiologic analyses: US rate and proportion files developed by the National Institute of Occupational Safety and Health and the National Cancer Institute. Am J Epidemiol 1992;136:85562.[Abstract]
-
Rothman KJ. Modern epidemiology. Boston, MA: Little, Brown & Co., 1986.
-
Breslow N, Day N, eds. Statistical methods in cancer research. Vol. 2. The design and analysis of cohort studies. Lyon, France: International Agency for Research on Cancer, 1987. (IARC scientific publication no. 82).
-
Statistical Analysis System. SAS user's guide: statistics. (Version 6.07). Cary NC: SAS Institute, Inc, 1991.
-
Gail M. Measuring the benefits of reduced exposure to environmental carcinogens. J Chronic Dis 1975;28:13547.[ISI][Medline]
-
Surveillance, Epidemiology, and End Results Program. SEER cancer statistics review (19731991). Washington, DC: Surveillance, Epidemiology, and End Results Program, 1994. (NIH publication no. 942789).
-
Axelson O, Steenland K. Indirect methods of assessing tobacco use in occupational studies. Am J Ind Med 1988;13:10518.[ISI][Medline]
-
National Health Survey. Smoking and tobacco use US, 1987. (Series 10, no. 169). Washington, DC: National Center for Health Statistics, 1989.
-
Steenland K, Deddens J, Salvan A, et al. Negative bias in exposure-response trends in occupational studies: modeling the healthy worker effect. Am J Epidemiol 1996;143:20210.[Abstract]
-
Kolstad HA, Olsen J. Why do short-term workers have high mortality? Am J Epidemiol 1999;149:34752.[Abstract]
-
Blair A, Hoar S, Walrath J. Comparison of crude and smoking-adjusted standardized mortality ratios. J Occup Med 1985;27:8814.[ISI][Medline]
-
Siemiatycki J, Wacholder S, Dewar R, et al. Degrees of confounding bias related to smoking, ethnic group, and socioeconomic status in estimates of the association between occupation and cancer. J Occup Med 1988;30:61725.[ISI][Medline]
-
Steenland K, Stayner L. Silica, asbestos, man-made mineral fibers, and cancer. Cancer Causes Control 1997;8:491503.[ISI][Medline]
-
Cherry N, Burgess G, Turner S, et al. Crystalline silica and risk of lung cancer in the potteries. Occup Environ Med 1998;55:77985.[Abstract]
Received for publication March 8, 2000.
Accepted for publication June 28, 2000.