Cumulative asbestos exposure and mortality from asbestos related diseases in a pooled analysis of 21 asbestos cement cohorts in Italy
Environmental Health volume 18, Article number: 71 (2019)
Despite the available information on cancer risk, asbestos is used in large areas in the world, mostly in the production of asbestos cement. Moreover, questions are raised regarding the shape of the dose response relation, the relation with time since exposure and the association with neoplasms in various organs. We conducted a study on the relationship between cumulative asbestos exposure and mortality from asbestos related diseases in a large Italian pool of 21 cohorts of asbestos-cement workers with protracted exposure to both chrysotile and amphibole asbestos.
The cohort included 13,076 workers, 81.9% men and 18.1% women, working in 21 Italian asbestos-cement factories, with over 40 years of observation. Exposure was estimated by plant and period, and weighted for the type of asbestos used. Data were analysed with consideration of cause of death, cumulative exposure and time since first exposure (TSFE), and by gender. SMRs were computed using reference rates by region, gender and calendar time. Poisson regression models including cubic splines were used to analyse the effect of cumulative exposure to asbestos and TSFE on mortality for asbestos-related diseases. 95% Confidence Intervals (CI) were computed according to the Poisson distribution.
Mortality was significantly increased for ‘All Causes’ and ‘All Malignant Neoplasm (MN)’, in both genders. Considering asbestos related diseases (ARDs), statistically significant excesses were observed for MN of peritoneum (SMR: men 14.19; women 15.14), pleura (SMR: 22.35 and 48.10), lung (SMR: 1.67 and 1.67), ovary (in the highest exposure class SMR 2.45), and asbestosis (SMR: 507 and 1023). Mortality for ARDs, in particular pleural and peritoneal malignancies, lung cancer, ovarian cancer and asbestosis increased monotonically with cumulative exposure. Pleural MN mortality increased progressively in the first 40 years of TSFE, then reached a plateau, while peritoneal MN showed a continuous increase. The trend of lung cancer SMRs also showed a flattening after 40 years of TSFE. Attributable proportions for pleural, peritoneal, and lung MN were respectively 96, 93 and 40%.
Mortality for ARDs was associated with cumulative exposure to asbestos. Risk of death from pleural MN did not increase indefinitely with TSFE but eventually reached a plateau, consistently with reports from other recent studies.
Asbestos fibres in their different mineralogical forms of chrysotile and amphiboles are a well known carcinogen acting on the respiratory tract and other organs. In 2009, the International Agency for Research on Cancer (IARC) updated the evaluation of asbestos fibres and confirmed that all types of asbestos cause malignant mesothelioma (MM), and cancer of lung, larynx and ovary in humans (Class 1), while the evidence was lower for pharynx, stomach and colorectal cancer (Class 2A) . The evidence of association shows a complex relation between the amount of exposure (dose and duration) and its time pattern (latency and time since cessation)  that still deserves clarification from large studies investigating dose, latency and other time related factors.
The asbestos cement industry was the largest asbestos consumer, using 85% of asbestos produced or imported in European countries . It employed a large number of workers: asbestos cement workers in Italy were estimated in 9000 in 1979 and 5000 in 1987, active in a large number of plants, most of which were of small or medium size [4, 5].
Several mortality studies on asbestos cement workers have been conducted in the past decades in Canada , United States , Israel , Italy [9,10,11,12,13,14,15,16,17,18], and in other European countries [19,20,21,22,23,24,25,26,27,28]. Only few studies analysed the mortality by gender [12, 14, 18, 19, 21, 25, 27, 28], generally including small numbers of female workers and observed deaths. Moreover, some relevant questions, such as the trend of the risk of mesothelial MN with cumulative dose and latency, are still debated [14, 29,30,31] and are of interest for both preventive and compensation purposes. The results of cohort studies on asbestos cement workers are also of interest for the countries where the product is still made. In all countries, the basic production process is the “Hatschek process”, first patented in late nineteenth century. The information on present asbestos cement industry are scanty and refer of poor working conditions even in recent years [32,33,34,35].
The present multicentre investigation is part of a larger project on epidemiological surveillance of asbestos workers . It aimed at investigating mortality by cause, gender, and time-related variables, in particular cumulative dose and latency, with particular attention to the neoplasms with sufficient or limited evidence of association with asbestos .
Material and methods
This study includes 21 asbestos cement factories, part of the Italian pooled cohorts of asbestos workers . The inclusion criteria for individual cohorts were availability of a data set with updated follow-up and a period of observation longer than 40 years. The main characteristics of the cohorts are reported in Additional file 1: Table S1.
The initial study base included 14,779 workers: 12,136 men (82.1%) and 2643 women (17.9%). During the preliminary analyses 11.5% of records (1703) were excluded, distributed over most of the cohorts, for the following reasons: first employment after 1992 (approval of the asbestos ban in Italy), conflicting dates, incomplete working periods, or unlikely age at hiring or end of employment (Hiring < 13 or > 70; End> 70). The two cohorts of the Eternit-Naples and Fibronit-Broni factories were limited to the workers hired after January 1st, 1950 because the quality of follow-up (FU) data was too limited before that date. Workers employed in more than one plant were identified, and their separate work histories were merged, except for cohort specific analyses. Figure 1 shows the flow chart for the construction of the pooled cohort. The final data set includes 12,578 workers, 81.7% men and 18.3% women.
Vital status was ascertained through the Registrar Office of the municipality of last residence. The FU date was December 31st, 2010 for the three cohorts of Tuscany and the cohort of Lazio, and December 31st, 2012 or later for the other 17 cohorts. For subjects lost to FU, the last contact date was used. Causes of death until 1985 were obtained from the municipalities, afterwards from Local Health Authorities (LHA). Causes were coded according to the 8th, 9th, or 10th Revisions of the International Classification of Diseases (ICD), according to the date of death. The ICD 10th revision is used in Italy from 2003 . Additional file 1: Table S2 presents the causes of death considered and the corresponding codes. We included ‘a priori’ in the analyses the causes of death associated with asbestos exposure, in particular asbestosis, malignant neoplasms (MN) of pleura and peritoneum, larynx, lung, and ovary, and the other neoplasms with limited evidence of association according to IARC classification . Causes of death reported in association with the Healthy Worker Effect (HWE) were considered in the analyses, in particular cardiovascular diseases, digestive diseases, respiratory diseases, psychiatric diseases, neurological diseases, digestive diseases and genitourinary diseases, as well as the total mortality . Cardiovascular diseases and respiratory diseases were also evaluated as related to smoking.
The production cycle of asbestos cement industries in Italy was similar in the different companies and was based on the “Hatschek process” . The dry mixture used for most productions contained about 13–15% of asbestos , mostly chrysotile from Italy, Canada and former Soviet Union. Amphibole fibres, mostly crocidolite from Australia and South Africa and a small percentage of amosite from South Africa, were also used depending on the product type. In the production of high-pressure pipes the percentage of asbestos was higher, reaching 20% (usually 1/3 of crocidolite, 1/3 of long chrysotile from Canada, 1/3 Balangero chrysotile). Asbestos exposure was measured in some companies at different times and the results have been reported in scientific papers [10,11,12,13,14, 16, 18], in courts during litigations and in reports from industrial hygiene laboratories of Italian public health agencies or universities. The reported levels of exposure were very high before 1980 (0.2–45 fibres/cc), intermediate between 1980 and 1989 (0.2–11) and lower after 1989 (< 0.1–0.3).
The main products were plain or corrugated sheets, pipes, slabs, tanks, chimneys and other products. In many factories, production also included the manual manufacture of small pieces. The asbestos cement manufacturing process can be reassumed in three stages: mixing, moulding and finishing. Asbestos exposure was high in mixing and finishing while in moulding area, exposure was lower due to the materials being wet. Diffusion of dust from higher to lower exposure areas was common, especially in the plants without physical barriers to separate the areas. These two main modes of exposure have been defined as “direct” for the high exposure areas and “indirect” for the lower ones. Indirect exposure regarded in particular the workers of the moulding workshops, while direct exposure affected the workers employed in mixing and finishing operations.
Two expert industrial hygienists (AA and SS) collected and evaluated the information regarding the use of asbestos, the work process, the plant layout and the measurements of airborne asbestos fibres, for each plant and year of activity. Considered data sources included both published and unpublished reports, in particular company reports, surveys of exposure, judicial examinations, and reports collected from workers . For each factory and calendar period, the original investigators had provided available information, including estimates of the proportion of exposed workers, the percentage of typical working time in tasks with asbestos exposure and the range of minimum and maximum level of asbestos concentration separately for direct and indirect exposures. Company-specific data were checked against those from other factories included in our pooled study and literature data to identify possible inconsistencies (i.e. differences that could not be explained by plant-specific features) and to fill data gaps (e.g. by deriving time-trends for asbestos concentration).
For each plant and periods, the experts estimated the proportion of workers exposed, the percentage of typical working time in tasks with asbestos exposure and the range of minimum and maximum concentration of asbestos airborne fibres (f/ml), separately for direct and indirect exposure.
Tasks and jobs of individual workers were not known, therefore plant and period-specific data were used to compute for each plant and year an Average Exposure Index (AEI) to be applied to all members of a given cohort. First, the ranges of concentration for direct and indirect exposures were summarized in single values by computing the geometric mean between minimum and maximum levels, adjusting for the average proportion of time in tasks with asbestos exposure. Such geometric means were taken as estimates of the time-weighted average level of exposure for direct and indirect exposures, specific for factory and period. The AEI value was calculated for each plant and year as the average of direct and indirect exposures, weighted by the respective proportional size of the workforce, according to the following formula:
where E = Exposure geometric mean, w = proportional size of the workforce and d = direct, i = indirect exposure, and p = plant, y = year.
From the AEI a Cumulative Average Exposure Index (CEI) was computed for the occupational history of each worker summing the contribution of all periods of activity:
Additionally, a fibre-type-weighted-AEI was computed based on the proportion of chrysotile (CHpy), amosite (Apy), and crocidolite (CRpy) used yearly in each plant. The weights were the MM potency factors for chrysotile, amosite and crocidolite (respective 1:14:71) as estimated by Hodgson and Darnton . The fibre-type-weighted-AEI provided the average chrysotile equivalent asbestos concentration in fibres per ml. It was computed as:
A fibre-type-weighted-Cumulative Exposure Index (fibre-type-weighted-CEI) was computed for each worker summing the fibre-type-weighted AEI for all periods of activity.
The dimension of AEI was a concentration, expressed in fibres/ml and the dimension of the fibre-type-weighted-AEI was the equivalent concentration of chrysotile asbestos fibres, in fibres/ml. CEI and fibre-type-weighted-CEI had the dimension of concentration times years (f/ml*year), the latter being the equivalent concentration of chrysotile asbestos fibres times years. Additional file 1: Table S3 data presents the exposure indexes by plant and period of activity.
Standardized Mortality Rates (SMR) were stratified by gender, and ‘a priori’ defined classes of calendar time, cumulative exposure, duration and Time Since First Exposure (TSFE) (also known as ‘latency’). For the purpose of present analyses, CEI was categorized in tertiles, based on the distribution of the total cohort. Class limits were: 5.0 and 48.5 f/ml*years for CEI (median 22.0); 54.0 and 620.0 (median 248.0) f/ml*years for fibre-type-weighted-CEI. Duration of exposure was calculated by summing up all the work periods since the date of first employment. TSFE was calculated from the date of first employment until the most recent date of observation. The same categories were used for Poisson regression models.
The numbers of expected deaths (used for SMRs) were based on regional mortality rates provided by the National Institute of Statistics – ISTAT (Rome, Italy) by cause, gender, and year, available from 1970 . Correspondingly, all analyses (including Poisson regression) were restricted to person years (PY) and events occurring after January 1st 1970; subjects lost to follow-up or dead before that date were excluded. For each cohort, expected deaths were calculated with the application of the corresponding regional rates. Multivariable analyses were carried out using Poisson regression . Chi square test for linear trend of SMRs was calculated when appropriate . 95% Confidence Intervals (CI) were computed according to the Poisson distribution of observed deaths . Statistical significance was set at 5%.
Models including restricted cubic splines (4 degrees) on fibre-type-weighted-CEI, TSFE and period were fitted to take into account the possible non linear relation between the predictive variables and mortality . In these models, fibre-type-weighted-CEI and TSFE were analysed as continuous variables.
Data were prepared using MS Access and SAS 9.2. Analyses were carried out using OCMAP plus, STATA 11, SAS 9.2 and R 3.2.5 with Survival and RMS packages .
Table 1 presents the distribution of cohort members (10,275 male and 2303 female workers) and person-years (PYs) (309,675 in men and 79,241 in women), by vital status, year, age at first exposure (defined as the beginning of employment period) and duration of employment. Forty four percent of men and 54.1% of women were alive at the end of FU; 54.4 and 44.9% died, and 1.3 and 0.9% emigrated or were lost at FU. The cause of death was known for 95.9% of male and 97.6% of female decedents. Year of first exposure was before 1970 for 70.2% of workers; 60.1% were younger than 30 at the beginning of activity; duration of employment was shorter than 20 years in 76%. Additional file 1: Table S4 presents the distribution of PY by age class and gender.
Table 2 presents mortality figures (observed and expected deaths, SMRs and 95% CI), by gender. Overall mortality and mortality for all MN were significantly increased in both sexes. Mortality was significantly higher than expected in both sexes for the MN of the respiratory tract category, and for MN of the lung, pleura, and peritoneum, and also for MN of unspecified site among men. Mortality from ovarian MN was higher than expected even if the excess was not statistically significant. Mortality from MN of larynx, pharynx, stomach, colon and rectum was not significantly higher than expected. Mortality from neoplasms of digestive tract excluding peritoneal MM was very close to expected in both genders (Men: SMR = 1.00; 95% CI: 0.92–1.09; Women: SMR = 1.18; 95% CI: 0.96–1.43). Mortality from asbestosis was extremely high (348 deaths in men - SMR: 507- and 51 in women -SMR: 1023), and determined an excess for the general categories of pneumoconiosis and respiratory diseases. Thirty three deaths were attributed to ‘other pneumoconioses’ (32 in men and 1 in women), including 21 silicoses and 12 ‘unspecified pneumoconioses’; the excess was statistically significant in men. In men, mortality was significantly lower than expected for neurological and cardiovascular diseases. In women for no causes of death mortality was significantly lower than expected. Mortality from unspecified causes (0.1% of total deaths in men and 0.2% in women) was increased in both sexes.
Mortality by tertile of cumulative exposure is presented in Table 3. The table also includes the p-value of the Chi-square test for linear trend. A statistically significant increasing trend was observed in both sexes for the MN of the pleura and of the peritoneum and for asbestosis, as well as for total mortality and total malignancies. The trend for lung cancer mortality was statistically significant in men, while it was irregular in women. Women showed a statistically significant trend also for ovarian cancer, and the SMR was increased in the upper tertile of cumulative exposure. A negative trend was observed in men for cardiovascular and for accidental deaths. The pattern of mortality from ‘other pneumoconioses’ by tertile of cumulative exposure suggested that at least a proportion of those cases are asbestosis cases: we observed an increasing trend, similar to that observed for the ‘asbestosis’ category but with lower SMRs: I tertile: 1 obs, 1.18 exp., SMR: 0.85 (95% CI: 0.02–4.72); II tertile: 6 obs, 2.46 exp., SMR: 2,44 (95% CI: 0.90–5.31); III tertile: 26 obs, 4.39 exp., SMR: 5.92 (95% CI: 3.87–8.68) (data not tabulated).
Mortality by duration of employment is presented in Additional file 1: Table S5. Results were similar to the trends observed for cumulative exposure and are not discussed in details.
Table 4 shows mortality in relation to TSFE (“latency”). The analysis was limited to the causes selected ‘a priori’ because associated with asbestos exposure, and to some categories related to the HWE. In both genders, ‘All causes’ and ‘All MN’ mortality was lower or equal to expected in the first 19 years of TSFE, and increased afterwards. No case of pleural MN was observed in the first 10 years of TSFE. The trend for pleural MN showed in both sexes increasing SMRs up to the 30–39 years class of TSFE, and a plateau or even a reduction at longer TSFE periods. No deaths were observed for MN of the peritoneum in the first two decades, while SMRs were increased from 20 to 29 years of TSFE in men and 40–49 years in women (p < 0.01 in both sexes) with no evidence of plateau. Lung cancer showed a trend similar to pleural MN: mortality was lower than expected in the first 10 years of TSFE, then increased to the highest SMR in the 30–39 year class, followed by a decrease. Laryngeal cancer mortality was below the expected values in the first 4 periods, and increased after 40–49 years of TSFE. Only 2 female cases of laryngeal cancer were observed. Deaths from ovarian cancer were observed only after 20 years of TSFE and the SMR was increased after 50 years of TSFE. There were no deaths from asbestosis in the first 10 years of TSFE in men and in the first 30 years in women. Subsequently, there was a continuous increase over the entire observation period, with SMRs higher in women. Mortality from cardiovascular diseases was lower than expected in all age categories for men but not for women.
Table 5 presents the results of Poisson regression analyses for pleural and peritoneal malignancies and lung MN. Several models were tested; the presented results correspond to the best fitting models, including different sets of variables according to the different diseases and sexes. All best fitting models included cumulative exposure (fibre-type-weighted-CEI). For pleural MN, the RR increased with cumulative exposure, while with TSFE it showed an increase in the first four decades, followed by a plateau, in both genders. For peritoneal MN, both genders showed a monotonic trend of increasing RR with increasing TSFE and cumulative exposure, albeit the small number of cases showed some irregularities for women. For lung neoplasm, in men RR showed an increase with cumulative exposure and a curvilinear trend for TSFE; for women, the model showed an increase of RR with cumulative exposure only from the lowest to the intermediate categories, with no further increase, and TSFE did not contribute to the model fit.
Figure 2 presents the cubic spline describing the relation of mortality rate from pleural (a), peritoneal (b) and lung MN (c) with cumulative exposure (fibre-type-weighted-CEI) and TSFE. For pleural neoplasm, mortality continuously increased with cumulative exposure, while it reached a plateau after 30 years of TSFE. Similar trends were observed for lung cancer, with a plateau becoming apparent after 40 years since first exposure. On the opposite, an increasing trend of mortality for peritoneal cancer was observed along the all range of fibre-type-weighted-CEI and TSFE.
Analyses for cumulative exposure were repeated using both the unweighted CEI and the fibre-type-weighted-CEI computed using the weights estimated by Garabrant et al. . None of these analyses (not in details) showed relevant differences with the analyses presented in Tables 3 and 5 and in Fig. 2 and all confirmed the dose response trends with increasing exposure.
Analyses by period of first exposure are reported in Additional file 1: Table S6. Among workers with first employment in 1980–89, a statistically significant excess was observed for pleural MN and for asbestosis (1 case) among men (employed in 1980, aged 36, until 1994), and one case of MN of the pleura was observed in women.
Our study is a pooled analysis of 21 cohorts of asbestos cement workers, which accounted for a wide fraction of workers employed in the asbestos cement production in Italy . At our knowledge, this is the largest study in the world on asbestos cement workers. Other specific features of this study include a very long follow-up, more than 40 years, and a significant number of women, providing a substantial contribution to gender-specific analyses.
The follow-up results are satisfactory, with only 1.3% of subjects with unknown status (lost or emigrated). Causes of death are known for over 95% of decedents in both sexes. SMR analyses are based on regional mortality rates, in order to increase the comparability between cohorts and reference rates. The decision to restrict the analyses to 1970 onwards depended on the availability of reference mortality rates  and not on the quality of cohort data, which is overall high. For consistency of the design and analyses, we decided to apply the restriction also to the internal analyses with Poisson regression, even if it was not strictly necessary.
Asbestos exposure could not be assessed on an individual basis, because of the lack of individual data on jobs and work activities for members of almost all cohorts. Therefore, we based our exposure assessment on the average index (AEI) representing the plant and period workforce average exposure, obtained from the estimation of exposure for workers with direct and indirect asbestos exposure. We calculated the individual cumulative exposure index (CEI) of cohort members by applying the AEI specific for plant and period to their duration and timing of employment. To take into account differences in the use of amphibole and chrysotile asbestos by plant and period, indices were weighted for the proportion of the different fibre types and their estimated carcinogenic potency factor for pleural MM . More recent estimates of fibre type potency were used in sensitivity analyses , although the data base used for that review was criticized for incompleteness , and no relevant difference was observed.
The study was based on mortality data, as cancer registries in Italy do not cover the population and period of interest for the study. There was no specific code for MM of peritoneum and pleura in the 8th and 9th ICD, and ICD 10 th is in use from 2003 in Italy. The use of mortality data could have caused a misclassification of MM with other cancers, in particular metastasis or lung cancer, in both directions. Therefore, literature evidence was searched to evaluate this possible bias. The sensitivity of death certificates for the identification of MM was explored in a meta-analysis by Kopylev et al. , who observed an underestimation of MM incidence from mortality data. This observation was supported by other studies not included in that revision: 74,5% of pleural MM cases could be identified from mortality records in Italy  and 87% in Southern England . Similar results were observed by Conti et al., who compared mortality and incidence for peritoneal MM in Italy . Some studies in our cohort [12, 14, 16] had performed a record linkage with the Italian Mesothelioma Registry data: results were satisfactory and showed that SMRs did not over-estimate SIRs of MM.
Mortality for “All causes” and “All MN” showed in both sexes a statistically significant increase, in general and in stratified analyses, in particular according to cumulative exposure and TSFE. The pattern was observed also for the workers who started employment in the more recent periods. This overall result, which is not affected by questions regarding classification of causes of deaths, shows macroscopically the consequences of the exposure in the asbestos-cement production industries. In total, the cohorts in study showed an excess of 1255.6 deaths, corresponding to a 19% increase of overall mortality. This large excess should be evaluated also with consideration of the evidence of a relevant HWE, as shown by the low mortality for cardiovascular and neurological diseases and by the low mortality observed in the first ten years of TSFE.
A statistically significant increase in mortality was observed in both sexes for pleural, peritoneal, and lung MN, consistent with the recent IARC evaluation : overall the number of deaths were 394, 133 and 858, and the corresponding attributable proportions  for these conditions were respectively 96, 93 and 40%, similar for men and women. The point estimate of SMRs in women are higher than in men for peritoneal and pleural MN. These results do not reflect a higher risk or an increased sensitivity to asbestos, but rather denote the lower female reference mortality rates in general population .
Analyses by cumulative exposure showed a statistically significant trend for pleural and peritoneal MN and for lung cancer, providing additional support to the evidence of dose response relation for MM  and for lung cancer [53, 54]. Poisson regression analyses showed similar RRs in men and women for pleural MN. For peritoneal MN, men had higher RRs than women at the same levels of cumulative exposure, however the number of female cases was limited and confidence intervals were wide and overlapping between the two group.
In our study, SMRs for both pleural and lung MN increased up to 30–39 years of TSFE and showed a plateau thereafter, with some fluctuations in women. On the contrary, the SMR for peritoneal MN increased monotonically in both genders. The same trends were observed also in the Poisson regression analyses using splines. As regards MM, our results contradict the traditional models of asbestos carcinogenesis, that predicted an unlimited increase of MM rates over time from the beginning of asbestos exposure, according to a power function of TSFE [55, 56] and did not differentiate between the anatomical location of MM [55, 56]. On the contrary, our results supported the evidence of a flattening in the trend of MM rates with TSFE and showed a difference between pleural and peritoneal MN, in agreement with the more recent evidence showing a reduction of the increasing incidence trend after long latency for pleural but not for peritoneal MM [14, 29,30,31, 57, 58]. The only exception in the recent literature was a follow-up study of former Polish asbestos workers participating in a health surveillance program, among whom a continued increase in the risk of pleural MM by TSFE was reported . The interpretation of such study is difficult as results may have been biased by the different participation of subjects with symptoms or otherwise interested in health assessment. A biological interpretation of the evidence is that asbestos fibres are slowly cleared from the lungs. For amphibole asbestos, a constant annual elimination rate in the order of up to 10% was estimated in different studies [60,61,62,63,64,65,66]. Toxicological  and histological studies  showed a different clearance of the fibres, especially in subjects with chrysotile only exposure. The same explanation holds also for lung cancer, showing results consistent with the reduction of risk following a reduction of lung burden. The difference between pleura and peritoneum is also consistent with reports from other studies [14, 29, 30, 68], and can be explained by the different distribution of asbestos fibers in the body, that interest the pleura more directly and the peritoneum only after internal transportation .
Ovarian cancer overall showed a modest increase in mortality, with SMR similar to the RR measured in meta-analyses [70, 71] and in a preliminary analysis of the present cohorts, including other industrial sectors as well as the asbestos cement . However, we observed a significant increase in women in the III tertile of cumulative exposure, in those who started employment before 1950, after more than 30 years of employment, or TSFE longer than 50 years.
We observed only a slight increase of laryngeal cancer mortality in both sexes. In men, the SMRs were higher in the higher classes of TSFE (40–49 and 50+ years) but we did not see a relation with cumulative exposure. In women, the two observed cases were in the mid and upper tertile of cumulative exposure and in the upper categories of TSFE. In summary, our cohort results provide little support to the association with laryngeal cancer, contrary to the two recent reviews by the Institute of Medicine  and by Peng et al. . Our study was limited by the use of mortality data and the absence of information on smoking and drinking. Mortality analyses are not very sensitive for diseases with long survival, such as laryngeal cancer in recent periods . A possible (negative) confounding by smoking and drinking habits cannot be ruled out: the mortality analyses showed a reduced mortality for cardiovascular deaths, suggesting that smoking habits were not worse than in the general population, although the evidence in this respect is limited.
Overall this large cohort did not show an excess of MN of digestive organs other than the peritoneal MN and provide little support to the causality of asbestos exposure for these neoplasms, contrary to recent results from meta-analyses or other cohorts [72, 75]. However, it should be noticed that for colon and rectum MN an excess was shown after TSFE longer than 50 years. As discussed for larynx cancer, these tumours have long survival and mortality could not be a sensitive indicator of risk.
For asbestosis, we observed in both sexes a clear association with cumulative exposure, a result that corroborates on the quality of our exposure assessment. A specific work on asbestosis in this cohort study is in progress, however present analyses already provide some relevant results. The increase in mortality from asbestosis was large, in particular among women. We believe that these data do not express a real gender difference but reflect the very low female reference mortality rates from these diseases in Italy . Five cases of death attributed to asbestosis were observed in men in the lowest tertile of cumulative exposure, corresponding to the cumulative exposure ‘up to 54.0’ f/ml*years for fibre-type-weighted-CEI. Mortality from asbestosis by period of first employment declined but remained higher than expected until the most recent periods, with one case observed among workers who started activity after 1980, possibly suggesting high exposures also in the last periods of industrial use of asbestos in Italy. Mortality from asbestosis was assessed on the basis of the underlying cause of death and no best-evidence assignment of the cause of death has been attempted. Previous reports from some of the cohorts included in the present study [14, 16], showed that additional cases had been reported when concomitant causes of death were considered, therefore an underestimation of SMRs is possible in this study. With consideration to the possible misclassification of outcome, we underline that 32 deaths in men and 1 in women had been recorded as “other pneumoconiosis” in the death certificate, possibly attributable to incorrect reporting of asbestosis, as SMRs increased with increasing cumulative exposure.
The possible occurrence of HWE was analysed with consideration of total mortality and of the causes of death reported in the literature as most often associated with it , namely cardiovascular, respiratory, neurological, psychiatric, digestive and genitourinary diseases. Overall mortality was higher than expected, because of the large increase of ARD mortality. The pattern of specific causes of death showed a marked HWE effect; in particular results for cardiovascular diseases mortality were corresponding to a HWE from selection at hire and possibly an additionally healthy worker survivor effect : cardiovascular disease mortality was not associated with cumulative exposure and remained lower than expected in the analyses by TSFE.
We do not have any individual information on smoking habits for the study, as often happens in occupational cohort investigations. However, the risk for cardiovascular diseases, strongly associated to smoking, is decreased in the cohort even after long TSFE, suggesting that the proportion of smokers was at least not larger than in the general population. Our findings on lung cancer are, therefore, unlikely due to uncontrolled confounding from smoking. We were, however, unable to assess the combined effect of tobacco smoking and asbestos [54, 76].
Our study contributes to the evidence of increased mortality risk due to asbestos exposure for MN of pleura, peritoneum, lung and ovary, as well as for asbestosis, all increasing with cumulative exposure. We have not detected an excess risk of larynx, pharynx, and stomach MN, except for larynx cancer after very long latencies. In summary, ARD mortality was higher than expected in both genders. Mortality from cardiovascular diseases was significantly lower than expected in men, suggesting an HWE but also that the prevalence of smokers among asbestos-cement workers was not larger than in the general population. Mortality for smoking-related neoplasms different from lung and larynx cancer was not increased.
Pleural and peritoneal MN showed a clear exposure-response relation and had different trends with TSFE. The TSFE analysis showed a monotonic increase of mortality in relation to TSFE for peritoneal MN, while pleura and lung MN mortality reached a plateau about 40 years after first exposure. The increase of MM rates by TSFE predicted by traditional models might be attenuated, due to the effect of the clearance of fibres.
These results underline the increased risk of ARDs in the asbestos cement production and are therefore of special relevance also for the countries that have not banned the use of asbestos.
Availability of data and materials
The datasets used and analysed during the current study are available from the corresponding author on reasonable request.
Average Exposure Index
Asbestos Related Disease
Cumulative Average Exposure Index
Healthy Worker Effect
International Agency for Research on Cancer
International Classification of Diseases
Local Health Authority
Standardized Mortality Rate
Time Since First Exposure
IARC International Agency for Research on Cancer. Arsenic, metals, fibres, and dusts. IARC Working Group on the Evaluation of Carcinogenic Risks to Humans. IARC Monogr Eval Carcinog Risks Hum. 2012;100(Pt C):11–465.
Health Effects Institute. Asbestos in public and commercial buildings: a literature review and synthesis of current knowledge. Cambridge: Health Effects Institute-Asbestos Research; 1991.
Albin M, Magnani C, Krstev S, et al. Asbestos and cancer: an overview of current trends in Europe. Environ Health Perspect. 1999;107(Suppl 2):289–98.
Mirabelli D, Kauppinen T. Occupational exposure to carcinogens in Italy: an update of CAREX database. Int J Occup Environ Health. 2005;11:53–63.
Mirabelli D, Cavone D, Luberto F, et al. Il comparto della produzione di Cemento Amianto. In: Marinaccio A, Binazzi A, Di Marzio D, et al., editors. Registro Nazionale dei Mesoteliomi Terzo Rapporto. Roma: Istituto Superiore per la Prevenzione e la Sicurezza sul Lavoro; 2010. p. 105–23.
Finkelstein MM. Absence of radiographic asbestosis and the risk of lung cancer among asbestos-cement workers: extended follow-up of a cohort. Am J Ind Med. 2010;53:1065–9.
Hughes JM, Weill H, Hammad YY. Mortality of workers employees in two asbestos cement manufacturing plant. Brit J Ind Med. 1987;44:161–74.
Tulchinsky TH, Ginsberg GM, Iscovich J, et al. Cancer in ex-asbestos cement workers in Israel, 1953-1992. Am J Ind Med. 1999;35:1–8.
Sarto F, Zambon P, Mastrangelo G, et al. Studio epidemiologico prospettico storico sulla mortalità per tumori di una coorte di soggetti esposti a cemento-asbesto. Epidemiol Prev. 1982;17–18:58–9.
Pettinari A, Mengucci R, Belli S, Comba P. Mortality of workers employed at an asbestos cement manufacturing plant in Senigallia. Med Lav. 1994;85:223–30.
Coviello V, Carbonara M, Bisceglia L, et al. Mortality in a cohort of asbestos cement workers in Bari. Epidemiol Prev. 2002;26:65–70.
Luberto F, Amendola P, Belli S, et al. Mortality study of asbestos cement workers in Emilia-Romagna. Epidemiol Prev. 2004;28:239–46.
Raffaelli I, Festa G, Costantini AS, Leva G, Gorini G. Mortality in a cohort of asbestos cement workers in Carrara, Italy. Med Lav. 2007;98:156–63.
Magnani C, Ferrante D, Barone-Adesi F, et al. Cancer risk after cessation of asbestos exposure. A cohort study of Italian asbestos cement workers. Occup Environ Med. 2008;65:164–70.
Ferrante D, Bertolotti M, Todesco A, et al. Mortality among asbestos cement workers: the cohort of the S.a.C.a. plant in Cavagnolo (Italy). Biomed Stat Clin Epidemiol. 2008;2:171–7.
Menegozzo S, Comba P, Ferrante D, et al. Mortality study in an asbestos cement factory in Naples, Italy. Ann Ist Super Sanita. 2011;47:296–304.
Fazzo L, Nicita C, Cernigliaro A, et al. Mortalità per cause asbesto-correlate e incidenza del mesotelioma fra i lavoratori del cemento-amianto di San Filippo del Mela (Messina). Epidemiol Prev. 2010;34:87–92.
Oddone E, Ferrante D, Tunesi S, Magnani C. Mortality in asbestos cement workers in Pavia, Italy: a cohort study. Am J Ind Med. 2017;60:852–66.
Thomas HF, Benjamin IT, Elwood PC, Sweetnam PM. Further follow-up study of workers from an asbestos cement factory. Brit J Ind Med. 1982;39:273–6.
Alies-Patin AM, Valleron AJ. Mortality of workers in a French asbestos cement factory, 1940-82. Brit J Ind Med. 1985;42:219–25.
Ohlson C-G, Hogstedt C. Lung cancer among asbestos cement workers. A Swedish cohort study and a review. Brit J Ind Med. 1985;42:397–402.
Gardner MJ, Winter PD, Pannett B, Powell CA. Follow up study of workers manufacturing chrysotile asbestos cement products. Brit J Ind Med. 1986;43:726–32.
Raffn E, Lynge E, Juel K, Korsgaard B. Incidence of cancer and mortality among employees in the asbestos cement industry in Denmark. Brit J Ind Med. 1989;46:90–6.
Albin M, Jakobsson K, Attewell R, Johansson L, Welinder H. Mortality and cancer morbidity in cohorts of asbestos cement workers and referents. Brit J Ind Med. 1990;47:602–10.
Szeszenia Dabrowska N, Wilczynska U, Szymczak W. Mortality of workers in two asbestos-cement plants in Poland. Int J Occup Med Environ Health. 2000;13:121–30.
Neuberger M, Vutuc C. Three decades of pleural cancer and mesothelioma registration in Austria where asbestos cement was invented. Int Arch Occup Environ Health. 2003;76:161–6.
Smailyte G, Kurtinaitis J, Andersen A. Cancer mortality and morbidity among Lithuanian asbestos-cement producing workers. Scand J Work Environ Health. 2004;30:64–70.
Sichletidis L, Chloros D, Spyratos D, et al. Mortality from occupational exposure to relatively pure chrysotile: a 39-years study. Respiration. 2009;78:63–8.
Barone-Adesi F, Ferrante D, Bertolotti M, et al. Long-term mortality from pleural and peritoneal cancer after exposure to asbestos: possible role of asbestos clearance. Int J Cancer. 2008;123:912–6.
Reid A, de Klerk NH, Magnani C, et al. Mesothelioma risk after 40 years since first exposure to asbestos: a pooled analysis. Thorax. 2014;69:843–50.
Pira E, Romano C, Violante FS, et al. Updated mortality study of a cohort of asbestos textile workers. Cancer Med. 2016;5:2623–8.
Berman DM. Asbestos and health in the third world: the case of Brazil. Int J Health Serv. 1986;16:253–63.
Giannasi F, Thébaud-Mony A. Occupational Exposures to Asbestos in Brazil. Int J Occup Environ Health. 1997;3:150–7.
Ansari FA, Bihari V, Rastogi SK, et al. Environmental health survey in asbestos cement sheets manufacturing industry. Indian J Occup Environ Med. 2007;11:15–20.
Burki T. Health experts concerned over India's asbestos industry. Lancet. 2010;375:626–7.
Ferrante D, Chellini E, Merler E, et al. Italian pool of asbestos workers cohorts: mortality trends of asbestos related neoplasms after long latency. Occup Environ Med. 2017;74:897–8.
Pirastu R, Ranucci A, Consonni D, et al. Reference rates for cohort studies in Italy: an essential tool in occupational and residential cohort studies. Med Lav. 2016;107:473–7.
Checkoway H, Pearce NE, Crawford-Brown DJ. Research methods in occupational epidemiology. New York: Oxford University Press; 1989. p. 78–9. and 86–91
Patroni M, Trimarchi R, Andreoletti F, Foà V, Chiappino G. Occupational risk in asbestos cement industry. Med Lav. 1987;78:351–9.
Stewart PA, Stewart WF. Occupational case-control studies: II. Recommendations for exposure assessment. Am J Ind Med. 1994;26:313–26.
Hodgson JT, Darnton A. Mesothelioma risk from chrysotile. Occup Environ Med. 2010;67:432.
Garabrant DH, Pastula ST. A comparison of asbestos fiber potency and elongate mineral particle (EMP) potency for mesothelioma in humans. Toxicol Appl Pharmacol. 2018;361:127–36.
Finkelstein MM. A comparison of asbestos fiber potency and elongate mineral particle (EMP) potency for mesothelioma in humans. Toxicol Appl Pharmacol. 2019;371:1–2.
Breslow NE, Day NE. Statistical methods in cancer research. Volume II - The design and analysis of cohort studies. Lyon: IARC Sci Publ 82; 1987. p. 1–406.
Harrell F. Regression modeling strategies: with applications to linear models, logistic and ordinal regression, and survival analysis. Berlin: Springer; 2015.
Team, R. C. R language definition. Vienna: R foundation for statistical computing; 2000.
Kopylev L, Sullivan PA, Vinikoor LC, Bateson TF. Monte Carlo analysis of impact of underascertainment of mesothelioma cases on underestimation of risk. Open Epidemiol J. 2011;4:45–53.
Bruno C, Comba P, Maiozzi P, et al. Accuracy of death certification of pleural mesothelioma in Italy. Eur J Epidemiol. 1996;12:421–3.
Okello C, Treasure T, Nicholson AG, Peto J, Møller H. Certified causes of death in patients with mesothelioma in south East England. BMC Cancer. 2009;9:28.
Conti S, Minelli G, Ascoli V, et al. Peritoneal mesothelioma in Italy: trends and geography of mortality and incidence. Am J Ind Med. 2015;58:1050–8.
Rothman K. Modern Epidemiology. 1st ed. Boston & Toronto: Little, Brown and co.; 1989. p. 38–40.
Novello S, Pinto C, Torri V, et al. The third Italian consensus conference for malignant pleural mesothelioma: state of the art and recommendations. Crit Rev Oncol Hematol. 2016;104:9–20.
Clin B, Morlais F, Launoy G, et al. Cancer incidence within a cohort occupationally exposed to asbestos: a study of dose--response relationships. Occup Environ Med. 2011;68:832–6.
Olsson AC, Vermeulen R, Schüz J, et al. Exposure-response analyses of Asbestos and lung Cancer subtypes in a pooled analysis of case-control studies. Epidemiology. 2017;28:288–99.
Newhouse ML, Berry G. Predictions of mortality from mesothelial tumours in asbestos factory workers. Br J Ind Med. 1976;33:147–51.
Peto J, Seidman H, Selikoff IJ. Mesothelioma mortality in asbestos workers: implications for models of carcinogenesis and risk assessment. Brit J Cancer. 1982;45:124–35.
Berry G. Models for mesothelioma incidence following exposure to fibers in terms of timing and duration of exposure and the biopersistence of the fibers. Inhal Toxicol. 1999;11:111–30.
Berry G, Reid A, Aboagye-Sarfo P, et al. Malignant mesotheliomas in former miners and millers of crocidolite at Wittenoom (Western Australia) after more than 50 years follow-up. Br J Cancer. 2012;106:1016–20.
Świątkowska B, Szeszenia-Dąbrowska N. Mesothelioma continues to increase even 40 years after exposure - evidence from long-term epidemiological observation. Lung Cancer. 2017;108:121–5.
Du Toit RSJ. An estimate of the rate at which crocidolite asbestos fibres are cleared from the lung. Ann Occup Hyg. 1991;35:433–8.
Churg A, Wright JL. Persistence of natural mineral fibers in human lungs: an overview. Environ Health Perspect. 1994;102(suppl 5):229–33.
Tossavainen A, Karjalainen A, Karhunen PJ. Retention of asbestos fibers in the human body. Environ Health Perspect. 1994;102(suppl 5):253–5.
De Klerk NH, Musk AW, Williams V, et al. Comparison of measures of exposure to asbestos in former crocidolite workers from Wittenoom Gorge, W. Australia. Am J Ind Med. 1996;30:579–87.
Berry G, de Klerk NH, Reid A, et al. Malignant pleural and peritoneal mesotheliomas in former miners and millers of crocidolite at Wittenoom, Western Australia. Occup Environ Med. 2004;61:e14.
Berry G, Pooley F, Gibbs A, Harris JM, McDonald JC. Lung fiber burden in the Nottingham gas mask cohort. Inhal Toxicol. 2009;21:168–72.
Merler E, Somigliana A, Girardi P, Barbieri PG. Residual fibre lung burden among patients with pleural mesothelioma who have been occupationally exposed to asbestos. Occup Environ Med. 2017;74:218–27.
Churg A, Wright JL. Fibre content of lung in amphibole- and chrysotile-induced mesothelioma: implications for environmental exposure. IARC Sci Publ. 1989;90:314–8.
Hodgson JT, Darnton A. The quantitative risks of mesothelioma and lung cancer in relation to asbestos exposure. Ann Occup Hyg. 2000;44:565–601.
Miserocchi G, Sancini G, Mantegazza F, Chiappino G. Translocation pathways for inhaled asbestos fibers. Environ Health. 2008;7:4. https://doi.org/10.1186/1476-069X-7-4.
Camargo MC, Stayner LT, Straif K, et al. Occupational exposure to asbestos and ovarian cancer: a meta-analysis. Environmen Health Perspect. 2011;119:1211–7.
Reid A, de Klerk N, Musk AW. Does exposure to asbestos cause ovarian cancer? A systematic literature review and meta-analysis. Cancer Epidemiol Biomark Prev. 2011;20:1287–95.
Institute of Medicine (US). Committee on Asbestos: selected Health Effects. Washington: National Academies Press (US); 2006. https://www.ncbi.nlm.nih.gov/books/NBK20326/
Peng WJ, Mi J, Jiang YH. Asbestos exposure and laryngeal cancer mortality. Laryngoscope. 2016;126:1169–74.
Gatta G, Botta L, Sánchez MJ, Anderson LA, Pierannunzio D, Licitra L, EUROCARE Working Group. Prognoses and improvement for head and neck cancers diagnosed in Europe in early 2000s: The EUROCARE-5 population-based study. Eur J Cancer. 2015;51:2130–43.
Boulanger M, Morlais F, Bouvier V, et al. Digestive cancers and occupational asbestos exposure: incidence study in a cohort of asbestos plant workers. Occup Environ Med. 2015;72:792–7.
Wraith D, Mengersen K. A Bayesian approach to assess interaction between known risk factors: the risk of lung cancer from exposure to asbestos and smoking. Stat Methods Med Res. 2008;17:171–89.
The authors wish to thank Prof. Loredana Musmeci from ISS for their contribution to the research.
The working group includes: Lisa Bauleo18, Antonio Baldassarre19, Carol Brentisci7, Barbara Cortini11, Stefania Curti20, Manuela Gangemi7, Giuseppe Gorini11, Patrizia Legittimo20, Francesco Marinelli20, Pasqualina Marinilli21, Stefano Mattioli20, Marina Musti22, Chiara Panato14, Venere Leda Mara Pavone21, Alessandra Ranucci2,3, Elisa Romeo18, Corrado Scarnato21, Cinzia Storchi1, Antonella Stura7, Simona Verdi11
1Epidemiology Unit, Azienda Unità Sanitaria Locale - IRCCS di Reggio Emilia, Reggio Emilia, Italy
2Unit of Medical Statistics and Cancer Epidemiology, Department of Translational Medicine, University of Eastern Piedmont, Novara, Italy
3CPO-Piedmont, Novara, Italy
7Unit of Cancer Epidemiology, CPO Piedmont and University of Turin, Turin, Italy
9National Cancer Institute IRCCS Fondazione Pascale, Naples, Italy
11Occupational & Environmental Epidemiology Unit - Institute for Cancer Research, Prevention and Clinical Network (ISPRO), Florence, Italy
14Mesothelioma Register of the Veneto Region, Regional Epidemiologic System, Local Health Unit 6, Padua, Italy
18Department of Epidemiology, Lazio Regional Health Service, Rome, Italy
19Occupational Medicine Unit, Careggi University Hospital, Florence, Italy
20Department of Medical and Surgical Sciences, University of Bologna, Bologna, Italy
21Department of Public Health, Prevention and Security Area Work Environments, Local Health Authority, Bologna, Italy
22Interdisciplinary Department of Medicine, Occupational Medicine “B. Ramazzini”, University of Bari, Bari, Italy
The cohort study was conducted thanks to grants from the ‘Asbestos Project’ of the Italian National Institute of Health (ISS), [Current research 2012: asbestos project. Operative Unit 2 “Amedeo Avogadro” University of Eastern Piedmont, Novara, Research Line 1,2.] and from the INAIL [Piano Ricerca 2016–2018, “Programma Speciale Amianto”, Ricerca BRIC id 55 and Ricerca BRIC id 59].
Ethics approval and consent to participate
The main study was approved by the University of Eastern Piedmont’s Ethical Review Board (authorization CE 112/13, July 12th, 2013), while the principal investigators (PI) of each cohort submitted the individual studies to the corresponding local boards.
The authors declare that they have no competing interests.
The following authors reported that they served as expert witness in court trials on asbestos related diseases: AA, AB, FB-A, LB, CM, LM, SM, EM, DM, MM, EO, VP, SS.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Asbestos cement pooled Italian study. Summary of cohort data, including all subjects, including those with follow-up terminated before 1970. Table S2. Asbestos cement pooled Italian study. Cause of death codes according to ICD in use at the date of death. Table S3. Asbestos cement pooled Italian study. Average asbestos exposure index, crude and adjusted for carcinogenic potency of asbestos types used, by cohort. Table S4. Asbestos cement pooled Italian study. Person-years at risk by gender and age class. Table S5. Asbestos cement pooled Italian study. Mortality by gender and duration of employment for selected causes of death. Table S6. Asbestos cement pooled Italian study. Mortality by gender and period of first exposure for selected causes of death. (PDF 646 kb)
About this article
Cite this article
Luberto, F., Ferrante, D., Silvestri, S. et al. Cumulative asbestos exposure and mortality from asbestos related diseases in a pooled analysis of 21 asbestos cement cohorts in Italy. Environ Health 18, 71 (2019). https://doi.org/10.1186/s12940-019-0510-6