Article Text

Original research
Relationships between asbestos exposure and pleural plaques: dose and time effects using fractional polynomials
  1. Morgane Menant1,
  2. Ilyes Benlala2,3,
  3. Isabelle Thaon4,5,
  4. Pascal Andujar6,7,8,
  5. Benoist Julia9,
  6. Patrick Brochard1,10,
  7. Christos Chouaid11,
  8. Bénédicte Clin12,13,
  9. Antoine Gislard14,15,
  10. Celine Gramond1,
  11. Christophe Paris16,17,
  12. Jean-Claude Pairon6,7,8,
  13. Fleur Delva1,10
  1. 1 Bordeaux Population Health EPICENE, Bordeaux, France
  2. 2 Faculté de Médecine, Université de Bordeaux, Talence, France
  3. 3 Service d'Imagerie Thoracique et Cardiovasculaire, CHU Bordeaux, Pessac, France
  4. 4 IINSPIIRE, Université de Lorraine, INSERM, Nancy, France
  5. 5 Centre de Consultations de Pathologie professionnelle, CHU de Nancy, Nancy, France
  6. 6 GEIC20 Team, INSERM U955, Créteil, France
  7. 7 Service de Pathologies Professionnelles et de l’Environnement, Centre Hospitalier Intercommunal Créteil, Institut Santé-Travail Paris-Est, Créteil, France
  8. 8 Institut Santé-Travail Paris-Est, Créteil, France
  9. 9 IIMTPIF, Créteil, France
  10. 10 Service Santé Travail Environnement, CHU de Bordeaux, Bordeaux, France
  11. 11 Service de Pneumologie, CHI Créteil 94000 Créteil, Créteil, France
  12. 12 Service de Santé au Travail et Pathologie Professionnelle, Caen, France
  13. 13 INSERM U1086, Caen, France
  14. 14 Faculté de Médecine, Normandie Université, UNIROUEN, UNICAEN, ABTE, Rouen, France
  15. 15 Centre de consultations de pathologie professionnelle, Rouen cedex, France
  16. 16 Service de Santé au Travail et Pathologie Professionnelle, Rennes, France
  17. 17 INSERM U1085, Rennes, France
  1. Correspondence to Fleur Delva, Bordeaux Population Health EPICENE, Bordeaux 33000, France; fleur.delva{at}chu-bordeaux.fr

Abstract

Objective The aim of this study was to confirm the relationship between several parameters of exposure to asbestos and pleural plaques (PP) using data from a large cohort of retired workers occupationally exposed to asbestos in France.

Method A large screening programme, including high-resolution CT (HRCT) examinations at inclusion and two other HRCT campaigns, was organised from 2003 to 2016 in four regions of France for voluntary, formerly asbestos-exposed workers. Exposure to asbestos has been evaluated by industrial hygienists based on the complete work history. The time since first exposure, the time since last exposure, Cumulative Exposure Index and maximum level of exposure to asbestos, were used in logistic regression using fractional polynomials to model the relationship with PP.

Results The study included 5392 subjects with at least one HRCT available. There was a significant non-linear effect of time since first exposure, time since last exposure and Cumulative Exposure Index to asbestos on the presence of PP. The risk of PP increased with increasing Cumulative Exposure Index to asbestos adjusted for time since first exposure, age and smoking status. Models also show that PP odds rise with increasing time since first exposure adjusted for cumulative index exposure, age and smoking status. PP odds decrease when time since last exposure increases.

Conclusion The study provides new data on the link between asbestos exposure and the presence of PP using fractional polynomials with non-linear relationships for time exposure parameters and asbestos exposure parameters.

  • Asbestos
  • Epidemiology
  • Occupational Health

Data availability statement

No data are available.

https://creativecommons.org/licenses/by/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

WHAT IS ALREADY KNOWN ON THIS TOPIC

  • Asbestos exposure is known to be associated with pleural plaques (PP).

WHAT THIS STUDY ADDS

  • There are non-linear relationships between time since first exposure and time since last exposure and the presence of PP.

  • Definition of populations for screening of lung cancer is largely debated at present time (level of smoke, age…) and parameters of asbestos exposure and existence of pleural plaques might also be considered in the definition of eligible subjects.

HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY

  • It appears important to carefully characterise the occurrence of PP to better define appropriate surveillance and screening programmes for these workers.

Introduction

Asbestos is a natural mineral with a fibrous texture that was mined and industrially used for many decades in the 20th century.1 Exposure to asbestos can lead to the development of respiratory diseases. Most often, it causes benign diseases such as asbestosis or pleural plaques (PP), although malignant diseases can also occur (mesothelioma and lung, larynx and ovarian cancers).1 Although its use has been banned in many countries since the 1990s, the WHO estimates that 125 million workers were still exposed to it worldwide in 2018.2

PP are circumscribed areas of fibrous thickening, typically of the parietal pleura and diaphragmatic pleura. They are the most common disease after asbestos exposure and are considered a marker of asbestos exposure.3 Although benign, consideration of PP may be important since an association has been reported between the presence of PP and pleural mesothelioma adjusted to asbestos exposure.4 In addition, studies on a French cohort of asbestos-exposed subjects (the Asbestos Diseases COhort (ARDCO) have shown an increased risk of lung cancer in the presence of PP, confirmed by high-resolution CT (HRCT) adjusted to asbestos exposure.5 6 However, this association remains controversial.7

Numerous studies based on chest X-rays (CXR) have shown a positive association between the presence of PP and time since first asbestos exposure (TSFE),8 9 intensity of exposure10–12 or cumulative exposure to asbestos.13 14 However, CXR-based imaging does not easily identify non-calcified PP, which may influence previous results from CXR-based studies. Thoracic HRCT is considered to be the most specific and sensitive tool for PP diagnosis.15 16 There are fewer HRCT-based studies investigating the same associations, with inconsistencies between some of these studies. Eisenhawer et al 17 reported a significant association between the duration of exposure to asbestos and the presence of PP, while Mastrangelo et al did not.18 Similarly, in 2012, Ameille et al reported a non-significant effect of TSFE on the presence of PP,19 while other studies reported the opposite.17 18 Paris et al reported a significant association between TSFE and cumulative exposure to asbestos and the presence of PP.20

In 2009, Paris et al reported a significant non-linear effect of TSFE and Cumulative Exposure Index (CEI) to asbestos for the presence of PP in a French cohort study of asbestos-exposed subjects based on HRCT.21 The present study reported new data, based on follow-up of the same cohort with several HRCTs and improved characterisation of PP grounded in a double reading by expert radiologists.

The aim of the study was to confirm the relationship between asbestos exposure and PP using data from a large cohort of retired workers occupationally exposed to asbestos using HRCT data: ARDCO.

Method

Study design

A screening programme was conducted from October 2003 until December 2005, and included retired workers from four regions of France (Basse-Normandie, Haute-Normandie, Aquitaine and Rhône-Alpes), previously occupationally exposed to asbestos. The main purpose of this cohort was to improve medical surveillance of workers formerly exposed to asbestos. The ARDCO21 were recruited in four regions of France through the media (unsolicited applications), mail invitations sent to beneficiaries of the asbestos workers’ allowance, to beneficiaries of post-professional follow-up not renewed for more than 2 years, to inactive workers over 55 selected according to their occupational code, to early retirees according to their declared occupational category and by healthcare professionals. Volunteers were included if their exposure to asbestos was confirmed by an industrial hygienist according to information from their complete occupational calendar. At inclusion, subjects were offered a physical examination, HRCT and pulmonary function tests. Two additional HRCT campaigns were conducted in 2010 and 2016. Thus, all subjects included in the study had at least one available HRCT stocked on CD-ROM and no missing value on asbestos exposure variables of interest and smoking status.

Data collection

A self-administered standardised questionnaire was used to collect socio-demographic characteristics (age, gender, sector of activity…), smoking status, and complete work history with title and duration for each job occupation mentioned.

HRCT scanning

For each campaign, all HRCTs were read by two radiologists with expertise in thoracic imaging, and by a third radiologist in the event of disagreement. The radiologists did not know the subject’s cumulative asbestos exposure and were asked to complete a standardised form, including the diagnosis of PP and interstitial anomalies.22 Circumscribed quadrangular pleural elevations, with sharp borders and tissue density, sometimes calcified, led to a positive diagnosis of PP. If a subject had several HRCTs available, the first one with PP was selected for the study. For a subject with no PP, the last available HRCT was selected.

Exposure assessment procedures

The complete work history of each subject was retrospectively analysed by industrial hygienists to evaluate asbestos exposure. All the information provided by the occupational questionnaire, including all successive jobs during the working life and specific questions on particular tasks, the context of the activities performed or contact with several asbestos-containing materials, was considered. The intensity of exposure was defined for each job and classified into four categories associated with a coefficient: low (0.01), low intermediate (0.1), high intermediate (1) and high (10). The maximum level of exposure for an individual corresponds to the maximum level found among all the individual’s jobs. A CEI to asbestos was then calculated for each subject over his/her working life, as the sum of exposures of each job (coefficient level × duration in years), expressed in exposure units × years since we had no metrology data. TSFE was defined as the time between the first year of exposure to asbestos and the date of the HRCT. Time since last exposure (TSLE) was defined as the time between the last year of exposure to asbestos and the date of the HRCT.

Statistical analysis

The relationship between the parameters of occupational exposure to asbestos and the presence of PP was estimated using logistic regression models adjusted for age at HRCT and smoking status (never smokers, ex-smokers—having stopped for at least 1 year and current smokers). In model I, the parameters included for occupational asbestos exposure were CEI and TSFE. In model II, the parameters were maximum level exposure and TSFE. Models III and IV were respectively the same as models I and II, to which the TSLE variable was added.

For the association between TSFE, TSLE, age and CEI, and the presence of PP, the null hypothesis of linearity against alternative regression functions was tested and the best fitting model was selected. To this end, we used fractional polynomial regression models. This technique allows to model the effect of a continuous variable (TSFE, TSLE, age and CEI) on the result (presence of PP) by a non-linear fractional polynomial function if necessary.23

In all the models, 10 years was used as the reference value for the TSFE, 0 years was used as the reference for the TSLE and 25 unit-years was the reference value for the CEI. These reference values correspond to an approximate value of the minimum observed in our study population.

Statistical analyses were conducted using R software V.3.6.1.

Results

Cohort selection

Of the 14 218 subjects included in the ARDCO cohort, 5808 subjects had at least one CT stored on CD-ROM. Subjects with incomplete data on smoking status (n=416) were excluded from the analysis. Thus, 5392 subjects were included in the study (figure 1).

Figure 1

Study flow chart. ARDCO, Asbestos Diseases Cohort; HRCT, high-resolution CT.

Descriptive analysis

The characteristics of the study subjects are presented in table 1. The mean age at the HRCT scan was 66.9 years. The majority of the subjects (95.7%) were men and most of them were ex-smokers (64.2%). Regarding asbestos exposure, the majority of subjects (43.0%) were estimated to have a high intermediate maximum exposure level. Subjects with PP had a mean CEI to asbestos twice as high as those without (104.9 vs 55.2 units-year). Moreover, 44.8% of the subjects with PP had a ‘high’ maximum exposure level, while only 22.8% of the subjects without PP had a ‘high’ maximum exposure level to asbestos. The average time since first exposure was 47.4±8.7 years.

Table 1

Study population (N=5392)

The prevalence of PP among the 5392 subjects was 25% (n=1348), reaching 39.6% among subjects with ‘high’ maximum exposure level.

Interstitial abnormalities (minor abnormalities, abnormalities incompatible with common interstitial lung disease, possible or certain common interstitial lung disease or asbestosis) were present in 833 subjects (15.4%) and, more specifically, in 540 subjects without PP (13.3%) and 293 subjects with PP (21.7%). If we consider only possible or certain common interstitial lung disease or asbestosis, they were present in 1.7% of subjects, and more precisely in 1.2% of subjects without PP and in 3.3% of subjects with PP.

Multivariate analyses

The results of the multivariate analyses show that all of the asbestos exposure parameters studied were associated with the presence of PP (table 2).

Table 2

Association between parameters of exposure to asbestos and pleural plaques (N=5392)

A non-linear effect of age, TSFE, TSLE and CEI on the presence of PP was observed and therefore modelised using fractional polynomials.

Model I shows that TSFE and CEI are significantly and independently associated with the prevalence of PP, adjusted for age and smoking status.

Specifically, Model I shows that PP odds rise with increasing CEI and reach a maximum for values of around 300 unit-years: OR300 VS 25 units-years=2.55 (95% CI 2.17 to 2.99). For values above 300 unit-years, the ORs for PP appear to remain roughly constant (figure 2).

Figure 2

Pleural plaques ORs as a function of Cumulative Exposure Index (CEI), reference=25 units-years. Model I: logistic regression model, effect of CEI to asbestos, and time since first exposure to asbestos on presence of pleural plaques, adjusted for age and smoking status, n=5392.

Model I also shows that PP odds increase with increasing TSFE up to a TSFE value of about 50 years: OR50 VS 10 years=8.22 (95% CI 3.18 to 21.26) (figure 3).

Figure 3

Pleural plaques ORs as a function of time since first exposure (TSFE), reference=10 years. Model I, logistic regression model, effect of cumulative index exposure to asbestos, and time since first exposure to asbestos on presence of pleural plaques, adjusted for age and smoking status, n=5392.

Model II shows that TSFE and the maximum levels of asbestos exposure ‘high intermediate’ and ‘high’ (reference ‘low’) are also significantly and independently associated with the prevalence of PP, adjusted for age and smoking status. The effect of TSFE on the presence of PP is similar to model I (online supplemental figure S1).

Supplemental material

Models III and IV show that PP odds decrease when TSLE increases (online supplemental figures S2 and S3) for TSLE values up to 30 years (reference 0), adjusted to CEI, TSFE, age and smoking status for model III and adjusted to TSFE, maximum exposure level, age, and smoking status for model IV.

Discussion

The main result of this study concerns the non-linear effect of asbestos exposure determinants associated with the presence of PP for both time (TSFE and TSLE) and dose (CEI or maximal exposure level) parameters. The study confirms the previous findings for this population,21 with a finer categorisation of the exposure, a longer follow-up of the subjects with an additional and double expert reading of HRCT, and a more precise characterisation of the effects. It suggests that TSLE is associated with a decrease in the likelihood of having PP.

Dose and time effects

The study showed a non-linear effect of CEI on the PP odds. The risk of PP increases with CEI and reaches a maximum before stabilising at around 300 units-years. Other HRCT-based studies have shown a non-linear effect of CEI on PP odds modelling, using classes of CEI.

Some studies, like ours, have a study population with heterogeneous exposure to asbestos, involving several sectors of activity. Mastrangelo et al showed a significant effect of CEI for dose classes above 159 fibres/mL-year compared with 43 fibres/mL-year adjusted for age and smoking status : OR43–159 fibres/mL-year = 1.99 (95% CI 1.20 to 3.28)).18 Murray et al showed that for every unit increase of asbestos exposure in fibres/mL-year, the plaque score increased by 0.074 (95% CI 0.048 to 0.099).24 The plaque score considered the circumference of the plaques in relation to the rib cage and their thickness, and was 0 in subjects without PP.

Other studies are based on cohort data with more homogeneous exposure. Rapisarda et al reported that the prevalence ratio of PP increased with increasing CEI in fluoro-edenite exposed workers.25 Eisenhawer et al 17 studied workers in electrical energy and reported a significant effect of CEI on PP odds for exposure above 25 fibres/year compared with one fibre/year, adjusted for smoking and age. However, in Eisenhawer et al’s study, this effect was not found in the analyses including temporal variables such as TSFE or exposure duration. Modelling a non-linear association using classes rather than functions such as fractional polynomials may explain the model’s difficulty in showing the effects of dose and of time.26 Lockey et al who studied workers exposed to Libby vermiculite also showed that cumulative exposure to asbestos fibres was significantly associated with pleural changes, including diffuse and localised pleural thickening for exposure above <0.15 fibre-year/cm3.27 Ameille et al, who studied transport workers with low exposure to asbestos, did not find a statistical association between CEI to asbestos and the presence of PP (p=0.26).28

Barbieri et al published autopsy analyses in 2019 that support the links demonstrated in this and other studies between cumulative asbestos exposure and the presence of PP.29

In addition, our study showed a non-linear effect of TSFE on the PP odds. PP odds increase with TSFE to reach a maximum of around 50 years. Graphically, we can thus observe a stabilisation. Interpretation of the apparent decrease in the effect of TSFE on the presence of PP for values beyond 60 years should be approached with caution given the number of subjects and thus the power of the model for these values. Other HRCT-based studies have investigated the effect of TSFE on the prevalence of PP. Eisenhawer et al showed a linear and significant effect of TSFE based on data from a cohort of electric power workers : OR=1.61 (95% CI 1.07 to 2.41).17 In a study of transport workers with low exposure to asbestos, Ameille et al did not find this association (p=0.56).28

Among studies with heterogeneous asbestos exposure study populations and data from HRCT, Mastrangelo et al found a non-linear effect of TSFE on the presence of PP modelled by classes: OR27–30 VS ≤26 years = 3.64 (95% CI 2.17 to 6.11).18

With regard to TSLE, the effect on the presence of PP presented in our study is adjusted for cumulative asbestos exposure or maximum exposure level, and is therefore significant irrespective of asbestos exposure dose. However, it is possible that a selective survival bias may have occurred. Indeed, compared with less-exposed individuals, those most exposed to asbestos are more likely to develop PP, but also have a higher risk of developing and dying from other diseases such as lung cancer or mesothelioma, and therefore of not being able to be followed up as long as those with less severe damage. It is therefore possible that, in this type of population, a longer TSLE is observed in subjects without PP than in those with PP, partly due to this survival bias. In such hypothesis, the conclusion about the association of TSLE and PP could not be applied to subjects highly exposed to asbestos.

Strengths and limitations

The present study modelled the non-linear effects of quantitative variables on the presence of PP using fractional polynomials, which allow for a more accurate account of their effect than classes.

Moreover, asbestos exposure measurement was assessed by industrial hygienists using the occupational calendar, which provided an expertise covering the subjects’ entire career. This collection method may involve memory bias. In the case of our study, the professional calendar data were collected before the CT examinations, we consider that the exposure history was not influenced by the diagnosis of PP.

The PP diagnosis was made by HRCT, with a double or triple reading in the event of disagreement between experts, helping to limit classification errors.

The prevalence of PP in all subjects in our study was 25%. Since exposure to asbestos of subjects in this study covered a wide range of industrial activities, the level of CEI was significantly variable. However, the prevalence of PP reached 39.6% in subjects with a ‘high’ maximum exposure level.

The cohort included subjects on a voluntary basis, so we can suspect that subjects with the highest exposure to asbestos were already being monitored by a physician or pneumologist with possible asbestos-related disease, as they were more likely to develop symptoms and so did not respond to inclusion. The results of our study may therefore not be transposable to a population highly exposed to asbestos. Finally, as HRCT was performed at a given point in time, the TSFE may have been overestimated in the event that the PP appeared before the diagnostic examination was performed. Thus, it is possible that the times after first exposure to asbestos presented in this study were slightly overestimated in a fraction of the subjects.

Conclusion

Based on a large French CT-based screening programme in asbestos-exposed subjects, the study provided new data on the link between asbestos temporal parameters and the presence of PP. Indeed, using fractional polynomials, it described non-linear relationships between TSFE and TSLE and the presence of PP. Notably, TSLE appears to show that cessation of exposure is associated with a decrease in the likelihood of developing PP. Time since cessation of exposure to asbestos could be considered in the follow-up of populations exposed to asbestos, in particular in designing the screening programme for lung cancer. Indeed, as a link is suspected between the presence of PP and excess risk of lung cancer in asbestos-exposed workers, it appears important to carefully characterise the occurrence of PP in order to better define appropriate surveillance and screening programmes for these workers.

Data availability statement

No data are available.

Ethics statements

Patient consent for publication

Ethics approval

The study was approved by the hospital ethics committee (CCPPRB Paris-Cochin n° 1946 2002), CCP Ile-de-France III, no 1946/11/02-02 (2010). All participants received information on the study and provided their written informed consent.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors MM performed literature review, all statistical analysis and drafted the first version of this manuscript. Conception or design of the work was done by CP, BC, PB, J-CP and FD. Acquisition and interpretation of data was done by MM, IB, CG, BJ, BC, AG, J-CP and FD. FD and J-CP supervised all aspects of this manuscript. MM, IB, IT, PA, BJ, PB, CC, BC, AG, CG, CP, J-CP and FD participated in the drafting, revision and correction of the final text. FD is acting as guarantor.

  • Funding This study was supported by the French Ministry of Labour and Social Relations; the French National Health Insurance (Occupational Risk Prevention Department); the French Agency for Food, Environmental and Occupational Health & Safety, Grant number: ANSES Grant 07-CRD-51 and EST 2006/1/43 and EST, and the Institut National Du Cancer SPADOC20-049.

  • Competing interests None declared.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.