Abstract
Most studies on asbestos-related diseases are based on chest radiographs, and dose–response relationships are still controversial. The aim of this study was to describe the most relevant parameters of asbestos exposure linked to pleural plaques and asbestosis diagnosed by high-resolution computed tomography (HRCT).
A large screening programme including systematic HRCT examinations was organised from 2003 to 2005 in France for formerly asbestos-exposed workers. The time since first exposure (TSFE), level, duration and cumulative exposure to asbestos were used in adjusted unconditional logistic regression to model the relationships of the two diseases.
Analysis of a population of 5,545 subjects demonstrated that TSFE (p<0.0001) and cumulative exposure (p = 0.02) (or level, depending on the models used), were independently associated with the frequency of pleural plaques. Only cumulative exposure (p<0.0001) or level of exposure (p = 0.02) were significantly associated with asbestosis. All trend tests were significant for these parameters.
To date and to our knowledge, this study is the largest programme based on HRCT for the screening of asbestos-exposed subjects. Both time–response and dose–response relationships were demonstrated for pleural plaques, while only dose–response relationships were demonstrated for asbestosis. These parameters must be included in the definition of high-risk populations for HRCT-based screening programmes.
Asbestos exposure can lead to the development of benign and malignant respiratory diseases. The commonest asbestos-related diseases are benign diseases and many studies have examined the relationships between asbestos exposure and these diseases 1. Overall, the prevalence of both pleural plaques and asbestosis is associated with time since first exposure (TSFE) to asbestos, intensity level, duration or cumulative exposure to asbestos depending on the studies 1–4, but these criteria remain poorly defined. Certain differences have been suggested between these two diseases, as TSFE seems to be the best predictor for pleural plaques whereas cumulative exposure was reported to be the major determinant for asbestosis 1, 2. Moreover, chest radiograph (CXR) features of asbestosis have repeatedly been shown to be associated with heavy asbestos exposure, with a cut-off close to 25 fibres·mL−1·yr−1 1.
However, the majority of published studies are based on CXR data, resulting in difficulties for estimation of dose–response relationships for asbestos-related diseases due to imprecise assessment of asbestos exposure but also imprecise radiographic diagnosis of asbestosis and pleural plaques 3. Thoracic high-resolution computed tomography (HRCT) has been clearly demonstrated to be more sensitive and specific than CXR for the diagnosis of asbestos-related diseases 5. To our knowledge, only four studies have used HRCT for the assessment of dose–response relationships for pleural plaques and asbestosis 6–9. Due to various study-dependent conditions, these studies did not provide evidence of clear relationships between asbestosis and pleural plaques and asbestos exposure. In 2001, Bégin and Christman 10 highlighted the need for an objective and independent measurement of disease activity, such as HRCT, but emphasised that computed tomography (CT) scanning was not yet a “gold standard”, despite its higher sensitivity compared with CXR, due to the absence of clearly established exposure–response relationships.
The aim of the present study was to describe the relationships between asbestos exposure and pleural plaques and asbestosis in a large cohort of formerly asbestos-exposed workers, and to assess asbestos exposure parameters linked to the presence of HRCT features of these two diseases by means of multivariate analysis.
METHODS
Study design
A large-scale screening programme for asbestos-related diseases, namely the Asbestos Post EXposure Survey (APEXS), was organised at the request of the French Ministry of Labour between October 2003 and December 2005 in four regions of France (Aquitaine, Rhône-Alpes, Lower Normandy and Upper Normandy). The study was initiated following a national consensus statement published in 1999 that recommended the use of chest CT scan for the surveillance of workers with previous moderate-to-high occupational exposure to asbestos 11. It is noteworthy that, in France, workers are entitled to financial compensation and early retirement in the case of asbestos-related disease, including pleural plaques. The target group for this screening programme comprised unemployed or retired asbestos-exposed workers covered by French National Health Insurance. Various recruitment procedures (letters, newspapers, radio, etc.) were used in each region.
All subjects with supposed occupational asbestos exposure were invited to undergo a free screening programme. The screening programme included physical examination, pulmonary function tests performed by the respiratory physician of the subject’s choice and spiral HRCT performed by a programme-approved radiologist. Medical data and a CD-ROM copy of the HRCT were sent directly to the investigators by the patient’s practitioner or radiologist. All data were then rated and recorded in a coordination centre for each region.
Exposure assessment procedures
Volunteers for the screening programme had to fill out a standardised questionnaire describing all jobs held throughout the subject’s working life as well as specific asbestos-exposing tasks. The questionnaires were analysed, blinded to radiological data, by industrial hygienists in Aquitaine and Normandy and by trained national health insurance agents in Rhône-Alpes. Asbestos exposure was assessed using an a priori simplified job-exposure matrix according to occupation and industrial activities. The level of exposure was defined for each subject’s entire career and classified into four classes: high (defined as continuous exposure for at least 1 yr or discontinuous exposure for at least 10 yrs), low (passive exposure), moderate (all other occupational exposure) and nil (no exposure). A cumulative exposure index (CEI) was then calculated by multiplying this level (3, 2, 1 and 0, respectively) by the total exposure duration (in years).
Radiological procedures
Recommended HRCT acquisition parameters were defined in accordance with the French Thoracic Imaging Society guidelines. The main recommendations were slices with a thickness of 1.5–5 mm, a pitch of 1.5–2.0, 120 kV, 60–150 mA maximum, window width of 1600 HU and window level of -600 HU for parenchymal slices, window width of 400 HU and window level of 50 HU for mediastinal slices, with differences depending on mono- or multi-detector equipment. Acquisition in the prone position was only required in the case of gravitational abnormalities. All CT scans were performed without contrast material.
Radiologists who participated in the programme (nearly 300) had to sign an agreement statement. They received guidelines to perform HRCT for the diagnosis of asbestos-related benign diseases. Pleural plaques were defined as focal areas of pleural thickening with soft tissue attenuation, which could contain calcification, typically located on the posterior wall of the lower half of the pleural spaces, often in the paravertebral regions, on the anterior wall between the third and the fifth costal interspaces, and/or on the diaphragmatic pleura. Their parietal nature was defined by their sharply defined edge, well demarcated from adjacent lung, usually with a typical tableland shape, easily recognised regardless of their thickness, or sometimes with less typical tapered margins. In the case of atypical pleural features, only bilateral pleural changes, located in typical areas, were taken into account. They had to be differentiated from the transversus thoracis muscle in anterior parasternal areas, and from subcostalis muscles along the posterior chest wall.
Diffuse pleural thickening was not defined by dimensional criteria but by the presence of parenchymal bands and/or rounded atelectasis, showing the visceral nature of the pleural changes.
Asbestosis was defined by the following pulmonary interstitial abnormalities initially affecting subpleural posterobasal pulmonary areas: early subpleural dot-like opacities, intralobular reticulation, subpleural curvilinear line, interlobular septal thickening and ground-glass opacity (the most common findings); and honeycombing (in the case of advanced asbestosis). These features had to affect both lungs and had to remain present on prone views to establish the diagnosis of asbestosis.
Radiologists also received specific training in the interpretation of HRCT from experienced radiologists and occupational physicians.
HRCT reports and films were both sent to coordination centres. For the purposes of the present study, the presence or absence of radiological abnormalities was rated by only using the HRCT reports according to the above definitions by physicians blinded to clinical data.
The project was approved by the Cochin Hospital ethics committee in Paris. All patients received information on the study and gave their written informed consent to the radiologist for the increased radiation dose delivered by HRCT.
Statistical analysis
The variables used to characterise asbestos exposure were duration of exposure, TSFE, level of exposure and CEI. All variables apart from level of exposure were ranked in classes following the quartile distribution of values. The level of exposure was ranked in four classes (nil, low, moderate and high).
Statistical associations between determinants of asbestos exposure and the presence of pleural plaques or asbestosis on HRCT were studied using stepwise unconditional logistic regression. As CEI was a combination of level and duration of exposure, two distinct models were used with either CEI or the other variables (i.e. level and duration of exposure).
In view of the low prevalence of diffuse pleural thickening, this disease was not specifically analysed, but these patients were included in all analyses concerning pleural plaques and asbestosis.
Crude odds ratios (ORs) and adjusted ORs for confounding factors (body mass index (BMI) and smoking status) were obtained separately for pleural plaques and asbestosis. Trend tests for logistic regressions were obtained using categorical variables as quantitative variables and by replacing ordinal variables by median values.
A cross-validation of a random sample of occupational questionnaires (n = 400) and HRCT reports (n = 200) was organised in the four regions, yielding agreements between hygienists or physicians rated as good to excellent by κ statistics.
All statistical analyses were performed with SAS Software version 9.1 (SAS Institute, Inc., Cary, NC, USA).
RESULTS
Cohort selection
A total of 20,157 subjects volunteered to participate in the APEXS and 16,885 subjects returned the completed occupational questionnaire; 13,859 (82.1%) of these subjects were eligible for the screening programme in terms of their National Health Insurance cover and more than half of them (7,275; 52.5%) underwent thoracic HRCT. Subjects whose HRCT reports were not sent to the coordination centre (n = 734) or who presented incomplete data (n = 709), as well as females (n = 255) and nonexposed subjects (n = 32), were excluded from analysis.
Descriptive analysis
Global characteristics of the APEXS population (n = 5,545) are given in table 1⇓. The mean±sd age was 63.5±5.7 yrs (range 39–85 yrs). Among the subjects in whom smoking status was known, 55% (2,066) were former smokers and only 9% (327) were current smokers. Only 26.2% (956) of the population had a normal BMI. The main jobs and industrial activities were metal working, construction and repairs.
Age, smoking status, body mass index and main occupational status# in the Asbestos Post EXposure Survey
Table 2⇓ describes the exposure characteristics. Overall, the mean±sd TSFE was 41.9±8.1 yrs with an average duration of exposure of 27.4±11.1 yrs. 80% (4,443) of the study population was estimated to have a high level of exposure throughout the entire career.
Asbestos exposure parameters in the Asbestos Post EXposure Survey
Based on HRCT data, the overall prevalence was 15.9% (882 out of 5,545) for pleural plaques, 6.8% (375) for asbestosis and 0.7% (37) for diffuse pleural thickening. A significant relationship was observed for both pleural plaques and asbestosis with age and a borderline significant relationship was observed for smoking status (table 3⇓). Higher prevalences of pleural plaques and asbestosis were observed in insulation workers, handlers or shipyard repair workers. Strong relationships were noted for pleural plaques and, to a lesser extent, for asbestosis with the four parameters of asbestos exposure (table 4⇓).
Prevalence of pleural plaques and asbestosis by age, smoking status, body mass index and main occupational status# in the Asbestos Post EXposure Survey¶
Prevalence of pleural plaques and asbestosis by asbestos exposure parameters in the Asbestos Post EXposure Survey#
Multivariate analyses: pleural plaques
On multivariate analyses, all exposure parameters except for duration of exposure showed positive relationships with the prevalence of pleural plaques (table 5⇓). Significant and independent associations were observed between pleural plaques and TSFE and CEI (model I) and between pleural plaques and TSFE and level of exposure (model II). Trend tests were also highly significant for TSFE (p<0.0001), CEI (p<0.0003) and level of asbestos exposure (p<0.0001). Taking into account the adjustment variables did not modify any of these relationships and had only a minimal effect on ORs.
Crude and adjusted odds ratios(ORs) for pleural plaques on high-resolution computed tomography in multivariate analyses using unconditional logistic regressions#
Asbestosis
Crude analyses of asbestosis (table 6⇓) showed a positive association with CEI (model I) and TSFE and level of exposure (model II). Only CEI (model I) and level of exposure (model II) remained significant after adjustment for smoking status. However, dose–response curves appeared to be weaker than for pleural plaques.
Crude and adjusted adjusted odds ratios(ORs) for asbestosis on high-resolution computed tomography in multivariate analyses using unconditional logistic regressions#
DISCUSSION
This large-scale screening programme using HRCT demonstrated strong relationships between asbestos exposure and pleural plaques and, to a lesser extent, between asbestos exposure and asbestosis. Exposure determinants associated with the presence of pleural plaques were both time parameters (TSFE) and dose parameters (CEI or level of exposure depending on the model used). A different pattern was observed for asbestosis, as only dose parameters were significantly associated with the prevalence of this disease. It is noteworthy that no relationship was observed between duration of exposure and these two asbestos-related diseases.
Numerous authors have studied the relationships between asbestos exposure and pleural plaques and the majority have reported the role of TSFE (as latency) in the prevalence of this disease 12–16. In a recent statement, the American Thoracic Society underlined that TSFE was the parameter most frequently correlated with this disease 1. For example, Koskinen et al. 12, in a large survey of Finnish construction, shipyard and asbestos industry workers, reported a stronger relationship for pleural plaques with TSFE than with duration of exposure. However, multivariate analyses were not performed in this study. In a radiological follow-up of 386 workers exposed to amosite, Ehrlich et al. 14 also found that TSFE was the best predictor of the prevalence of pleural plaques in multivariate analysis, with a significant OR of 1.67 (95% confidence interval 1.15–2.42) by decade of the variable. Järvholm 17 consequently included TSFE to model the relationship between asbestos exposure and prevalence of pleural plaques. Conversely, Copes et al. 18, in a study in Quebec chrysotile workers, found that exposure peaks might be the best exposure determinant of pleural plaques. In a multivariate analysis, Finkelstein and Vingilis 19 suggested that cumulative exposure was the best predictor of pleural plaques. It is noteworthy that all these results are based on CXR. As mentioned in the introductory section, only a few studies using HRCT for the diagnosis of pleural plaques have been published. Algranti et al. 7 found a relationship between TSFE and pleural plaques in 828 former asbestos cement workers, but no multivariate analysis was available in their study. Schaeffner et al. 8, in a series of 102 cases of lung cancer, found a significant relationship between CEI and the simultaneous presence of pleural plaques and asbestosis. Van Cleemput et al. 9, in a subset of 51 subjects, showed no statistical relationship between cumulative exposure to asbestos or TSFE and the surface area of pleural plaques, as ascertained by HRCT. The present study demonstrates, for the first time, that TSFE and CEI (or level of exposure) are both significantly and independently associated with the prevalence of pleural plaques. These findings partially agree with some of the previous studies discussed and discrepancies with published data may be explained by the statistical power and the use of multivariate analysis in the present study, and by the sensitivity of HRCT compared with CXR for the diagnosis of pleural plaques.
Several studies in which the diagnosis of asbestosis was based on CXR found significant associations between duration of exposure or TSFE and this disease 12, 20, 21, but the strongest associations were generally reported with cumulative exposure 1, 18, 19, 22, by studies including two with multivariate analyses 14, 23. The study conducted by Algranti et al. 7 using HRCT reported a stronger association between TSFE and asbestosis than between CEI and asbestosis. In a previous and independent study including 706 subjects screened by HRCT 6, it was found that CEI was the best predictor of early asbestosis. The current results are consistent with the hypothesis that cumulative exposure is more important than TSFE for the development of asbestosis 4, as no significant relationship was observed for TSFE in adjusted models. These results are in accordance with biological data on retention of asbestos 24. However, although significant correlations were observed between CEI (or level of exposure) and asbestosis, these correlations appeared to be weak. One explanation for this feature may be the low prevalence of asbestosis observed in this cohort (6.8%), which is in agreement with recent data on the time-course of asbestosis epidemiology 25.
This study presents a number of limitations. First, the study population was selected for screening purposes and only volunteers were included. The variety of occupations and industrial activities described in this cohort were similar to those reported in other studies conducted in the general population in France 26. However, despite the large number of subjects included in this study, the study population may not be representative of the general population.
Secondly, a possibly more important weakness is the procedure used for reading and rating CT scans. In order to validate our results, a random subsample of 1,600 CT scans was independently reviewed by two of the panel of 12 physicians, with a third reading in the case of disagreement. Unilateral pleural plaques were considered relevant only in the case of typical pleural features and, for the diagnosis of pulmonary asbestosis, a bilateral extent of the interstitial abnormalities was needed. This is in accordance with publications dealing with asbestos-related pleural abnormalities and with the weakness of the data about unilateral asbestosis 27, 28. After review, the same models were confirmed for pleural plaques, but were not significant for asbestosis due to an insufficient number of cases (data not shown). An independent re-interpretation of all CT scans with the same procedure and definitions has since been undertaken and is at present ongoing.
Thirdly, as all asbestos exposure assessments were based exclusively on assessments by industrial hygienists, the accuracy of the parameters used may be questionable, particularly the CEI. Dose–response curves for the prevalence of pleural plaques and asbestosis assessed by HRCT were recently published for a distinct population with available atmospheric measurements 29. This study corroborated our present findings, albeit the previous study also found a slight but significant association between TSFE and asbestosis prevalence.
A fourth limitation of the study is the presence of numerous missing data, especially for confounding factors such as smoking status or BMI. However, comparisons of subjects with or without these data failed to demonstrate any marked differences between the groups in terms of the distribution of asbestos exposure parameters (data not shown). Moreover, the role of smoking status on parenchymal abnormalities or BMI on pleural plaques appears to predominantly concern results based on CXR analyses 30, 31 rather than CT 32.
Finally, as emphasised by Stayner et al. 23, one of the major limitations of this type of approach to the modelling of exposure is the uncertainty of the models. We used log-linear models to describe the relationship between asbestos-related diseases and exposure. This model is widely used for this purpose in the literature 19 and allows the results to be expressed simply in terms of ORs and numerical classes. However, Stayner et al. 23 reported that the best predictor model for asbestosis appears to be the power model and Järvholm 17 also used a similar model to estimate the prevalence of pleural plaques. In accordance with these recommendations, we re-analysed our data using this alternative model and did not observe any differences in terms of significant variables in all models used for either pleural plaques or asbestosis.
Conclusion
To date, this study is the largest programme based on HRCT for the screening of asbestos-exposed subjects. We report, for the first time, strong and independent correlations between time parameters of asbestos exposure (as expressed by TSFE) and pleural plaques, and between dose parameters (namely cumulative exposure or level of exposure) and both pleural plaques and asbestosis. These findings indicate that these time and dose parameters must be included in the definition of high-risk populations if further screening programmes are considered.
Support statement
This work was supported by grants from the French National Health Insurance (Dept of Occupational Risk Prevention) and the French Ministry of Labour and Social Relations (both Paris, France).
Statement of interest
None declared.
Acknowledgments
The members of the APEXS are: J. Ameille (Dept of Occupational Diseases, AP-HP, Raymond Poincaré Hospital, Garches, France), B. Aubert and J.L. Rehel (IRSN, Fontenay-aux-Roses, France), J. Baron, P. Brochard and P. Reungoat (LSTE, Dept of Occupational Diseases, University Hospital, Bordeaux, France), J. Benichou (Dept of Biostatistics, University Hospital, Rouen, France), A. Caillet, N. Lestang, M. Letourneux and M. Pinet (Dept of Occupational Diseases, University Hospital, Caen, France), P. Catilina (Dept of Occupational Diseases, University Hospital, Clermont-Ferrand, France), S. Chamming’s and F. Conso (IIMTPIF, Paris, France), G. Christ de Blasi, A. Gislard, M. Maurel, A. Porte and L. Thorel (Dept of Occupational Diseases, University Hospital, Rouen), E. Guichard and E. Schorle (ERSM Rhone-Alpes, Regional Security Insurance, Lyon, France), B. Millet and F.X. Thomas (Dept of Occupational Risk Prevention, Regional Security Insurance, Lyon), L. Mouchot, C. Paris and A. Stoufflet (INSERM U954, Vandoeuvre-les-Nancy, France), J.C. Pairon (Dept of Occupational Diseases, Centre Hospitalier Intercommunal, Creteil, France) and A. Sobaszek (Dept of Occupational Diseases, University Hospital, Lille, France). The authors would like to thank all members of the programme for their contribution to this survey.
The authors also thank A. Saul (Paris, France) for his valuable advice in editing this manuscript.
- Received June 20, 2008.
- Accepted November 17, 2008.
- © ERS Journals Ltd