- Open Access
- Open Peer Review
Mortality and morbidity study of petrochemical employees in a polluted site
Environmental Healthvolume 11, Article number: 34 (2012)
The area of Gela was included among the 57 Italian polluted sites of national interest for environmental remediation because of its widespread contamination from a petrochemical complex. The present study investigates mortality and morbidity of the cohort of Gela petrochemical workers with the aim of disentangling occupational from residential risk.
Mortality was assessed for 5,627 men hired from 1960, year of the plant start-up, to 1993; it was followed up for vital status in the period 1960–2002. Morbidity was analysed for 5,431 workers neither dead nor lost to follow-up from 1960 to 2001 and was based on Hospital Discharge Records in the period 2001–2006. The work experience was classified in terms of job categories such as blue collars, white collars, and both – workers who shifted from blue to white collar (95%) or vice versa. An ad hoc mobility model was applied to define qualitative categories of residence in Gela, as residents and commuters. Standardized Mortality Ratios (SMRs) and Mortality Rate Ratios (MRRs) were computed, the latter by using a Poisson regression model. Morbidity was analyzed in terms of Hospital Discharge Odds Ratios (HDORs) through a logistic regression model. While performing the internal comparisons, white collars was the reference category for the job analysis, and commuters was the reference category for the residential analysis.
In the light of epidemiological evidence about health risk from petrochemical industries in both occupational and environmental settings, and/or on the basis of information about occupational and residential contamination and health risk in the area of Gela, noteworthy results are shown for lung cancer [MRR: 2.11 (CI 90%; 0.96-4.63) in blue collars; 1.71 (1.09-2.69) in residents], respiratory diseases [HDOR: 2.0 (1.0-3.0) in blue collars; 1.4 (0.96-2.06) in residents] and genitourinary diseases [HDOR: 1.34 (1.06-1.68) in blue collars; 1.23 (1.04-1.45) in residents].
The results support a role of the exposures in the occupational and residential settings, the latter due to the local ascertained contamination, in affecting the workers’ health. These results underline the urgent need of water, soil, air and food-chain monitoring programs, to discover active sources of exposure and consequently define public health interventions.
Gela is a town of approximately 80,000 inhabitants located by the sea in the South West coast of Sicily Region, in Southern Italy. Crude oil was found in the underground in the late ‘50s; since the early ‘60s a vast area of Gela municipality, nearby Gela town, has hosted a large oil refinery, together with a thermoelectric power plant and petrochemical plants for production of organic (ethylene, acrylonitrile) and inorganic (sulphuric acid, ammonia, chlorine, urea) chemicals. From 1998, a part of Gela municipal area — the entire petrochemical complex and an extended sea portion — was included among the 57 Italian Polluted Sites (IPSs) of national interest for environmental remediation. Data collected in this site after year 2000 documented groundwater, soil and air contamination. Maximum concentrations of arsenic, vinyl chloride, mercury and benzene were some orders of magnitude higher than the threshold values [1, 2]. It was documented that soil and shallow water in Gela IPS were severely contaminated by metals and by organohalogenated compounds . On the basis of the Italian National Registry of Emissions and Sources data , the emissions of NOx, SOx, and benzene in Gela area represented respectively 11%, 30%, and 31% of the total emissions occurring in Sicily. Furthermore, two monitoring studies of pine needles and road dust samples showed that Gela town was heavily affected by metals and metalloids [4, 5]. Finally, a recent biomonitoring study in a group of Gela residents showed high levels of total arsenic in urine . Urinary arsenic levels higher than reference values, although indicative of recent exposures , are of interest because inorganic arsenic or its metabolites may cause adverse health effects including lung cancer [8, 9].
In 1970, the case of Gela was described as an example of “industrialization without development” , because the expected evolution of local cultural and societal structures, figured on the basis of the capitalistic industrialization in Northern Europe and North America, had not taken place. The situation in Gela recalled the one that can be seen in developing countries, where an industrialization process externally promoted, not deriving from the evolution of the local society, could not induce local territorial and population development , but produced environmental contamination. In this context, the implementation of environmental monitoring and reclamation activities can be hampered. In 1995 a Decree called for an air quality monitoring program ; still, at the end of 2010 such program was not implemented yet . As far as the health status of residents is concerned, since the mid ‘90s ecological mortality and morbidity studies showing probable health risks from environmental and occupational exposures were conducted [14–16].
After year 2000, the public prosecutor in Gela started an inquiry on potential industrial pollution-related health risks concerning both the workers and the general population. Workers employment rosters were seized, and a mortality study of Gela petrochemical workers was carried out [17, 18]. The cohort study was successively endorsed and funded by the WHO European Centre for Environment and Health as part of the Project “Technical assistance to the Sicilian Region, Special Office for risk areas to develop plans for environmental remediation”, in three Sicilian areas at high environmental risk, including Gela. The present paper aims at studying mortality and morbidity of Gela petrochemical workers to disentangle occupational from residential health risks.
Workers employed in the petrochemical plant in Gela from January 1, 1960, namely from the start of plant operation, to January 31, 2002, were 7,147 (6,961 men, 186 women). The cohort considered in the present study comprised 5,627 workers for the mortality analysis and 5,431 for the morbidity analysis.
In the mortality study, all men hired from 1960 to 1993 and born in Sicily were included. The vital status was ascertained for the period 1960–2002 through a request to the municipalities’ registry office; the workers hired after 1993 were excluded to allow a minimum 10 years latency time for delayed exposures effects . The causes of death were coded by an expert nosologist according to International Classification of Diseases (ICD) ruling in the year of death. The vital status was unknown for 4.1% of the total cohort and the cause of death was missing for 27 deaths (4.6%) in the restricted cohort.
In the morbidity study, 5,431 men hired from 1960 to 2000 not deceased or lost to follow-up, and born in Sicily were included. The Hospital Discharge Records (HDRs) available for the Sicily Region in the period 2001–2006 were analysed; the outcome was hospitalization. For each subject the first hospitalization for any given cause of a priori interest was selected; the information on main diagnosis was used in the analysis. The Sicilian HDRs Registry covers all the residents’ hospital admissions in or out of their own Region. For each HDR, the fiscal code (a univocal identification code for Italian citizens) was reconstituted. HDRs containing no valid personal data were excluded. Day-hospital admissions, rehabilitation, long stays and duplicate records were also excluded. HDRs were linked to cohort subjects using fiscal code as a key.
Twenty employment rosters of eight different companies were used to collect workers’ personal and employment data. Data from most of the workers were available in different employment rosters; in each roster, besides personal data, the date of start of employment and the job category in terms of blue collar or white collar were registered. The majority of workers were employed in different companies inside the complex; the duration of employment could not be determined as the only reliable information was the hiring date, and the termination date of employment was not recorded. No data were available about occupational exposures at individual or group level.
Workers had different origins: 85% of them came from Sicily and approximately 15% from other Regions . The latter were hired chiefly in the first 10 years of plant operation , mainly to train unskilled local workers , and usually lived in Gela for short periods; then they moved to other places in Italy where some petrochemical plants newly constructed needed skilled labour to start operation . Workers born outside Sicily were excluded from the present study to ensure that Sicilian population was an appropriate comparison in the analysis with an external reference, and to minimize exposure heterogeneity and unknown sources of potential confounding.
Outcomes of a priori interest
Some main diseases groups and single specific diseases, for which the epidemiological evidence had shown an association with occupation in petrochemicals or residence in the neighbourhoods, were analysed [21, 22]. The following neoplastic (a) and non-neoplastic (b) diseases were selected: a) malignant neoplasms of stomach, colon-rectum, liver, trachea, bronchus and lung, pleura, skin, bladder, kidney, central nervous system (CNS), Hodgkin disease, non-Hodgkin lymphoma, multiple myeloma, leukemias; b) chronic obstructive pulmonary disease (COPD), asthma, nephritis, nephrotic syndrome and nephrosis.
Job and residence categories
Workers were classified into one of the following three job categories: only white collars if they were registered only as white collar in all the employment rosters, only blue collars if they were registered only as blue collar in all the employment rosters, both white and blue collars if they changed from blue collar to white collar or vice versa in different employment rosters; from now on white collars, blue collars, both.
The reconstruction of the residential history was unfeasible. The companies’ rosters provided information on residence at hiring for about 90% of the workers, but this information was not used as it was specific for a given point in time. Moreover, the hiring might have happened before the official change of residence (and usually firms do not need to update internal documents). The companies’ rosters registered also workers’ birthplace abstracted from the identity document or driving license; the birthplace was available for almost all the workers (98.5%). This information had already been used in a previous study to classify workers as residents in Gela or commuters, but the criterion to define commuters was not robust . In the present study the workers were classified as residents in Gela or commuters using the procedure described as follows.
Firstly, birthplace data were used on the basis of two assumptions: a) the workers born in Gela were likely to have been residents in Gela at least during the employment period; b) the workers born in Sicilian municipalities other than Gela might have been or not been residents in Gela during the employment period. The assumption b) was based on anecdotal information from several documents [10, 23] successively confirmed by a qualitative sociological study [11, 20].
Secondly, a gravity-type mobility model was adopted to identify commuting workers. The mobility model has been described elsewhere ; it was implemented to estimate the probability that a "central locality", where a specific plant is established, attracts commuting flows from other places. The general model was applied to the Gela case, thus defining the probability of commuting to Gela from other Sicilian municipalities .
Finally, the mobility model results were applied to define two qualitative categories of likelihood of residence in Gela:
Commuters - workers born in Sicilian municipalities, excluding Gela with a probability of commuting defined by the model as >=0.5 (N. 3,234);
Residents - workers born in Gela (N. 1,684) and in Sicilian municipalities with a probability of commuting defined by the model as <0.5 (N. 709). The latter were fixed as “moved to Gela at hiring”.
To compute Standardized Mortality Ratios (SMRs) the expected number of deaths was based on the Sicilian death rates specific for gender, 5-year age and calendar time periods. Mortality data are registered by the Italian National Bureau of Statistics, and are electronically available starting from 1970; therefore the mortality rates for the period 1970–75 were applied to the period 1960–1970. The unknown causes of death contributed to the SMRs calculation; that is, a distribution of unknown causes equal to the known ones was assumed and the number of unknown deaths was added to each group on a proportional basis .
For the internal analysis, Mortality Rate Ratios (MRRs) were calculated applying a Poisson regression model using as predictive variables job and residence categories, age, calendar period, and time since first employment. Hospital Discharge Odds Ratios (HDORs) were calculated applying a Logistic regression model, using as predictive variables job and residence categories, age, and time since first employment. White collars was the reference category for job and Commuters was the one for residence. A multiple imputation strategy was adopted for missing data on job category (7% in the mortality study; 5% in the morbidity study) using 20 imputed datasets in both the regression models [25, 26] under the reasonable assumption of their missing at random (MAR). A sensitivity analysis made by attributing to different job categories the subjects with missing job information, did not show substantial differences in the main results.
Only causes with at least 10 observed cases were considered for the internal analysis.
The Confidence Intervals (CIs) were estimated by the maximum likelihood method. The 90% level was chosen to limit the use of CIs as significance test thus avoiding the selection of relevant results on the sole basis of “statistical significance” .
The statistical analysis was performed using STATA 11.0 software .
Table 1 shows the main descriptive data of the cohort.
In the main diseases groups, the SMRs analysis shows that the observed mortality is lower than expected for most of the causes of death (Table 2). A slight increase is observed for primary liver cancer and kidney malignant neoplasms; the increase for asthma is negligible.
Table 3 shows the mortality results by job and residence categories.
For the majority of the analysed causes, MMRs are higher among blue collars than among workers belonged to the both job category. Among blue collars an increase is detected for all causes, all malignant neoplasms and digestive system malignancies. For colon and lung cancer, disease of respiratory system and chronic obstructive pulmonary diseases (COPD) an increase is suggested. An increased risk with imprecise estimates for the category of all lymphatic and haematopoietic tissue neoplasms is observed in occupational and residential mortality analysis.
MRRs for residents are increased for all malignant neoplasms of the respiratory system and lung cancer, and decreased for genitourinary neoplasms.
Table 4 presents the results for Hospital Discharge Records by job category and residence.
An increased risk among blue collars for all neoplasms, diseases of circulatory system, COPD, diseases of genitourinary system, nephritis and nephrosis is shown. An increase is suggested for respiratory malignancies.
Among workers in both job category an increase in bladder cancer is observed. In the same category an increase for nephritis is suggested.
HDORs among residents are augmented for CNS neoplasms and genitourinary diseases; a decreased risk is shown for skin and bladder cancer. For kidney cancer and COPD an increase is suggested.
The procedures adopted for cohort enumeration, follow-up, vital status ascertainment and causes of death coding complied with recommended standards in occupational epidemiology . Based on a previous validity study , Sicilian Region’s HDRs were retained of good quality and suitable for this study.
The main limits of the present investigation are the lack of information on individual residential history, occupational and environmental exposures. Sometimes defining occupational exposure in quantitative terms in a retrospective cohort study is not possible . This could be a major problem in analysing the risk for workers in petrochemical complexes, as the one in Gela, where chemical mixtures are common, and their intensity and duration vary over-time . Moreover, when the population is not characterized with respect to specific exposures , the risk estimates for some groups could suffer from dilution effect; consequently, identifying risk increases could be more difficult. Finally, in the mortality study of Gela petrochemical workers, the MRR for chronic and long latency diseases might be underestimated because the median age at the end of follow-up is about 59 and 78% of the workers are younger than 65.
The observed deficit in SMR estimates can be mainly attributed to the Healthy Worker Effect (HWE). This effect is common in occupational cohorts and typically characterized by lower mortality from all combined causes of death, cardiovascular, respiratory, and neoplastic diseases . The HWE is a phenomenon detected in different countries and in cohorts experiencing different exposures [32–36]. It has been suggested that the HWE bias can be minimized by the selection of an appropriate comparison such as a working population [32, 36] or by analytical methods such as lag time analyses or analyses restricted to newly employed or long term workers . In the present study data limitations did not allow the application of such approaches.
Work and residence categorizations suffer from intrinsic data limitations. As in most occupational cohort studies, data on confounding from lifestyle habits such as smoking, were not available. As for smoking, its prevalence is usually higher in blue collars than in white collars, but this has not been documented in the petrochemical sector where smoking is usually strictly prohibited in several working areas (in Gela plant since 1963). Though, sensitivity analyses have shown that plausible smoking differences will rarely result in a lung cancer “relative risk” greater than 1.5 and will be even lower for diseases that are less strongly associated with smoking . This was also confirmed by several analyses with empirical data in occupational studies [39, 40].
For the both job category, changes from blue collar to white collar (95%) or vice versa might not always correspond to a change in the work environment because it could only be a change in the professional level.
The job classification adopted in the present study has consequences for the internal analysis which deserve specific comments. On one hand, an overestimation of MRR in blue collars and an underestimation in the both job category result from the classification of person-years; on the other hand, an underestimation of risk estimates in blue collars cannot be excluded if they had a major risk due to the work environment, because changes from blue collar to white collar might not always correspond to a change in the work environment.
The categorization of residence is intrinsically affected by misclassification, since the likelihood that commuters have been resident in Gela is greater than the likelihood that residents have been resident elsewhere; this is because the majority of residents (70%) were defined only on the basis of place of birth, while all commuters were identified using both information on place of birth and applying a mobility model that might include further uncertainties in the classification. As a consequence, if differences in risk are really due to a residential risk in Gela, when comparing residents vs. commuters underestimation of MRRs and HDORs are expected. For some causes, confounding could be due to residential risky exposures (e.g. urban air pollution) other than to industrial pollution in Gela, but there are no data to infer this point. As far as confounding from lifestyle habits, it is a minor problem in residential analysis because there are no reasons to assume differences between residents and commuters. Internal residential analyses are adjusted for job categories, thus taking into account, at least in part, possible confounding factors associated with the occupational context. Moreover, confounding for some causes might still occur because of differences in lifestyles between blue collars and white collars.
Causes selected for analysis in the present study differ in terms of prognosis and duration. Some conditions are highly fatal, others have a more favourable prognosis and may cause repeated hospital admissions. For these reasons, mortality and morbidity data are complementary, and provide a deeper and more realistic description of the population health profile.
The restriction criteria adopted for the morbidity analysis implied the exclusion of subjects deceased before 2003, possibly resulting in a survival cohort selected because of good health status. This could lead to an underestimation of the occupational risk, when performing the morbidity analysis, because of major risk observed in the mortality study in the blue collars category. The cohort was restricted to workers born in Sicily; the average age was 55 at the start of the morbidity follow-up. This indicated that most of them were likely to be resident active workers in Sicily in the period at study. The number of the subjects possibly lost to the morbidity follow up because of migration outside Sicily, might have been negligible and there are no reasons to believe that the hypothetical losses were differential for the job and residence categories.
The periods of hiring and follow-up considered in the mortality study coincide with the ones of previous analyses [17, 18]. Nevertheless, the present investigation has novelties in several elements: restriction only to Sicilian workers, application of a robust mobility model to define residential categories; consideration of all causes of a priori interest; computation through a multiple imputation approach of missing data on job categories. Furthermore, in the present study results from the mortality study were considered together with the ones of the morbidity analysis, to disentangle occupational from residential risk. Results of previous mortality studies, whose internal analyses were limited to all causes, all cancer and lung cancer, are consistent with the present ones.
In the light of the observed increased risks, some considerations and comments on specific diseases follow.
The available epidemiological evidence of health effects of residence in the neighbourhood of petrochemical plants is coherently pointing to an increased risk for respiratory diseases, including lung cancer . An augmented risk for lung cancer among petrochemical workers was observed in some occupational cohort studies . The main risk factor for lung cancer is tobacco smoke, but also environmental and occupational factors can contribute to its etiopathogenesis . In the present study, the increase in lung cancer mortality among Gela residents and among blue collars point to a causal link with environmental and occupational exposures respectively. The increased occurrence of lung cancer among blue collars is also observed in hospital discharge analysis.
As for lung cancer, the main risk factor for COPD is tobacco smoke, but outdoor pollution and occupational exposure may contribute to the risk of developing the disease [42, 43]. COPD mortality is increased in Gela residents and in blue collars; the latter also show higher HDORs.
Lymphohematopoietic tissue neoplasms showed a 10-15% increased risk for employment in chemical industry . In 1989, the International Agency for Research on Cancer classified occupational exposures in petroleum refining as probably carcinogenic to humans (Group 2A), on the basis of limited evidence for skin cancer and leukaemia . Subsequent petrochemical workers cohort studies confirmed an increased risk for lymphoematopoietic tissue neoplasms . Evidence for residential risk for leukemia is suggested by some studies, while it is scarce and conflicting for other lymphohematopoietic cancers . The present study suggests a moderate increase in mortality, but not in hospital discharges, from lymphohematopoietic neoplasms in the study population.
Because of documented air and water contamination from mercury and other heavy metals in the area of Gela , and given the presence of a chlor-alkali unit in the petrochemical plant in the period 1970–1994 , renal diseases are also of particular relevance in the present investigation. In fact, heavy metals, namely mercury, cadmium and lead, have an ascertained toxic effect on the kidney . Exposure to mercury might have contributed to the increase in nephritis, nephrotic syndrome and nephrosis among blue collars, shown by hospital discharges.
According to Wrensch et al. review of primary brain tumours , the question of whether employment in petrochemical, petroleum and oil production industry entails an increased risk remains unresolved. Nevertheless, a relationship between certain parts of the industry and increased risk is documented by the fact that most studies have shown an increased mortality of 20% to 80% . In recent years residential exposure to petrochemicals was investigated in studies with conflicting results [48, 49]. The increased risk of CNS cancer in Gela residents shown by hospital discharge analysis might thus be associated with environmental exposures.
For bladder cancer, tobacco smoking is the major established risk factor. Occupational exposure to aromatic amines (for example 2-napthylamine and benzidine) and polycyclic aromatic hydrocarbons are found to increase the risk . Other potential risks are represented by working in dyestuff manufacture, painting, rubber, leather, and aluminium industry. Carcinogens in drinking water (chlorination by-products and arsenic) may increase risk of bladder cancer . It has been estimated that exposure to occupational carcinogens causes up to 25% of all bladder tumours . Some evidence of risk is arising in the petrochemical industry . The increased hospital discharge odds ratio described by the present study among subjects who had been both blue and white collars, but not among blue collars alone, is not easily interpretable.
In the light of epidemiological evidence about health risk from petrochemical industries in both occupational and environmental settings, and/or on the basis of information about occupational and residential contamination and health risk in the area of Gela, noteworthy results are shown for lung cancer, respiratory and genitourinary diseases. They support a role of exposures in the occupational and residential settings, the latter due to the local ascertained contamination, in affecting the workers’ health.
These results underline the urgent need of water, soil, air and food-chain monitoring programs, in order to look for active sources of exposure and to define adequate preventive strategies. In particular, as recommended since 1995 , an ad hoc air pollution program should be implemented.
central nervous system
chronic obstructive pulmonary disease
Hospital Discharge Odds Ratios
Hospital Discharge Records
Healthy Worker Effect
International Classification of Diseases
Italian Polluted Sites
Mortality Rate Ratios
Standardized Mortality Ratios.
Musmeci L, Bianchi F, Carere M, Cori L: Environment and health in Gela (Sicily): state of the knowledge and prospects for future studies. Epidemiol Prev. 2009, 33 (suppl 1): 1-159.
Directive 2006/118/EC of the European Parliament and of the Council of 12 December 2006 on the Protection of groundwater against pollution and deterioration. Official Journal of the European Union. L 372/19
Italian National Inventory of Emissions and Sources. http://www.eper.sinanet.apat.it/site/it-IT/. (in Italian)
Bosco ML, Varrica D, Dongorrà G: Case study: inorganic pollutants associated with particulate matter from an area near a petrochemical plant. Environ Res. 2005, 99: 18-30. 10.1016/j.envres.2004.09.011.
Manno E, Varrica D, Dongorrà G: Metal distribution in road dust samples collected in an urban area close to a petrochemical plant at Gela, Sicily. Atmos Environ. 2006, 40: 5929-41. 10.1016/j.atmosenv.2006.05.020.
Bianchi F, Bustaffa E, Cori L, Imbriani M, Minichilli F, Migliore S, Minoia C, Ronchi A, Turci R: Human biomonitoring in the area around the petrochemical site of Gela, Sicily, Italy [abstract]. EUROEPI2010. Epidemiology and Public Health in an Evolving Europe. Epidemiol Prev. 2010, 82-suppl 1
Hughes MF: Biomarkers of exposure: a case study with inorganic arsenic. Environ Health Perspect. 2006, 114: 1790-6.
Smith AH, Marshall G, Yuan Y, Ferreccio C, Liaw J, von Ehrenstein O, Steinmaus C, Bates MN, Selvin S: Increased mortality from lung cancer and bronchiectasis in young adults after exposure to arsenic in utero and in early childhood. Environ Health Perspect. 2006, 114: 1293-6. 10.1289/ehp.8832.
Straif K, Benbrahim-Tallaa L, Baan R, Grosse Y, Secretan B, El Ghissassi F, Bouvard V, Guha N, Freeman C, Galichet L, Cogliano V: WHO International Agency for Research on Cancer Monograph Working Group. A review of human carcinogens--part C: metals, arsenic, dusts, and fibres. The Lancet. 2009, 10: 453-454. 10.1016/S1470-2045(09)70134-2.
Hytten E, Marchioni M: Industrialization without development. 1970, Franco Angeli, Gela a southern history. Milan, (in Italian)
Saitta P: Spaces and Societies at Risk. Ecology, Oil and Change in Gela, Sicily. 2009, Think Thanks, Naples, (in Italian
Decree of the President of the Republic 17.01.1995: Remediation plan for the rehabilitation of the area in the Province of Caltanissetta, South-eastern Sicily. G.U. No 02/05/1995 100, ordinary supplement No 51. (in Italian)
Settimo G, Mudu P, Viviano G: Air pollution: general problems in the area of Gela. Epidemiol Prev. 2009, 33 (suppl 1): 37-42.
Martuzzi M, Mitis F, Biggeri A, Terracini B, Bertollini R: Environment and health status of the population in areas with high risk of environmental crisis in Italy. Epidemiol Prev. 2002, 26 (suppl 6): 1-53.
Fano V, Cernigliaro A, Scondotto S, Addario SP, Caruso S, Mira A, Forastiere F, Perucci CA: Mortality (1995–2000) and hospital admissions (2001–2003) in the industrial area of Gela. Epidemiol Prev. 2006, 30: 27-32.
Cernigliaro A, Pollina Addario S, Cesaroni G, Fano V, Fantaci G, Tavormina E, Marras A, Dardanoni G, Forastiere F, Perucci CA, Scondotto S: Health status in high risk areas in Sicily. Update of mortality (1995–2002) and morbidity (2001–2006). 2008, Sicilian Epidemiological Observatory Bullettin, Palermo, in Italian
Pasetto R, Biggeri A, Comba P, Pirastu R: Mortality in the cohort of workers of the petrochemical plant in Gela (Sicily) 1960–2002. Epidemiol Prev. 2007, 31: 39-45.
Pasetto R, Comba P, Pirastu R: Lung cancer mortality in a cohort of workers in a petrochemical plant: Occupational or residential risk?. Int J Occup Environ Health. 2008, 14: 124-8.
Checkoway H, Pearce N, Kriebel D: Research methods in occupational epidemiology. 2004, Oxford University Press, Oxford, 2
Pasetto R, Saitta P, Bracci C: Development of an epidemiological study in a polluted site: the unexpected contribution of sociological and medico-legal approaches. Epidemiol Prev. 2008, 32: 325-8.
Pirastu R, Ancona C, Iavarone I, Mitis F, Zona A, Comba P: SENTIERI Project –Mortality study of residents in Italian polluted sites: evaluation of the epidemiological evidence. Epidemiol Prev. 2010, 5-6 (Suppl 3): 1-96.
Pirastu R, Pasetto R: Review of epidemiological evidence on health effects of residence near petrochemical plants. Human health in areas with industrial contamination. Edited by: Mudu P, Terracini B, Martuzzi M. 2012, WHO Regional Office for Europe, Copenhagen, in press
Vasta C: Gela…than the petrochemical plant came. 1998, Tipografia Lussografica, Caltanissetta, (In Italian)
Signorino G, Pasetto R, Gatto E, Mucciardi M, La Rocca M, Mudu P: Gravity models to classify commuting vs. resident workers. An application to the analysis of residential risk in a contaminated area. Int J Health Geogr. 2011, 10: 11-10.1186/1476-072X-10-11.
Sterne JA, White IR, Carlin JB, Spratt M, Royston P, Kenward MG, Wood AM, Carpenter JR: Multiple imputation for missing data in epidemiological and clinical research: potential and pitfalls. BMJ. 2009, 338: b2393-10.1136/bmj.b2393.
Stata Multiple Imputation reference manual. 2009, StataCorp LP College Station, Texas
Sterne JA, Davey Smith G: Sifting the evidence-what's wrong with significance tests?. BMJ. 2001, 322: 226-31. 10.1136/bmj.322.7280.226.
StataCorp: Stata Statistical Software: Release 11. College Station. 2009, StataCorp LP, Texas
Nicolosi A, Perriera S, Di Gaetano V, Di Gaetano P, Allegra G, Montagna M, Barone AP, Buffa S, Drago S, Bruno G, Chibbaro RM, Colucci AT, Arcà M: Sicilian Regional Program of evaluation of hospital discharge records validity. 2007, Epidemiological Observatory Bullettin, Palermo, in Italian
Kriebel D, Wegman DH, Moure-Eraso R, Punnett L: Limitations of meta-analysis: cancer in the petroleum industry. Am J Ind Med. 1990, 17: 269-71. 10.1002/ajim.4700170212.
Pearce N, Checkoway H, Kriebel D: Bias in occupational studies. Occup Environ Med. 2007, 64: 562-8. 10.1136/oem.2006.026690.
Valenty M, Tirmarche M, Mitton N, Laurier D, Gelas JM, Quesne B: Causes of death of COGEMA active workers, 1980–1995. Rev Epidemiol Sante Publique. 2003, 51: 461-8.
Su WL, Chen YH, Liou SH, Wu CP: Meta-analysis of standard mortality ratio in cotton textile workers. Eur J Epidemiol. 2004, 19: 989-97. 10.1007/s10654-004-0917-3.
Blair A, Sandler DP, Tarone R, Lubin J, Thomas K, Hoppin JA, Samanic C, Coble J, Kamel F, Knott C, Dosemeci M, Zahm SH, Lynch CF, Rothman N, Alavanja MC: Mortality among participants in the agricultural health study. Ann Epidemiol. 2005, 15: 279-85. 10.1016/j.annepidem.2004.08.008.
Mikoczy Z, Tinnerberg H, Björk J, Albin M: Cancer incidence and mortality in Swedish exposed to ethylene oxide: updated cohort study findings 1972–2006. Int J Environ Res Public Health. 2011, 8: 2009-19. 10.3390/ijerph8062009.
Thygesen LC, Hvidtfeldt UA, Mikkelsen S, Brønnum-Hansen H: Quantification of the healthy worker effect: a nationwide cohort study among electricians in Denmark. BMC Public Health. 2011, 11: 571-10.1186/1471-2458-11-571.
Stellman SD, Resnicow K: Tobacco smoking, cancer and social class. Social inequalities and cancer. Edited by: Kogevinas M, Pearce N, Susser M, Boffetta P. 1997, International Agency for Research on Cancer, Lyon, 229-50. IARC Scientific Publications No 138
Axelson O: Confounding from smoking in occupational epidemiology. Br J Ind Med. 1989, 46: 505-507.
Kriebel D, Zeka A, Eisen E, Wegman D: Quantitative evaluation of the effects of uncontrolled confounding by alcohol and tobacco in occupational cancer studies. Int J Epidemiol. 2004, 33: 1040-5. 10.1093/ije/dyh151.
Blair A, Stewart P, Lubin JH, Forastiere F: Methodological issues regarding confounding and exposure misclassification in epidemiological studies of occupational exposures. Am J Ind Med. 2007, 50: 199-207. 10.1002/ajim.20281.
Alberg AJ, Ford JG, Samet JM, American College of Chest Physicians: Epidemiology of lung cancer: ACCP evidence-based clinical practice guidelines (2nd edition). Chest. 2007, 132: 29S-55S. 10.1378/chest.07-1347.
Pauwels RA, Buist AS, Calverley PM, Jenkins CR, Hurd SS, GOLD Scientific Committee: Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease. NHLBI/WHO global initiative for chronic obstructive lung disease (GOLD) workshop summary. Am J Respir Crit Care Med. 2001, 163: 1256-76.
Anto JM, Vermeire P, Vestbo J, Sunyer J: Epidemiology of chronic obstructive pulmonary disease. Eur Respir J. 2001, 17: 982-94. 10.1183/09031936.01.17509820.
Greenberg RS, Mandel JS, Pastides H, Britton NL, Rudenko L, Starr TB: A meta-analysis of cohort studies describing mortality and cancer incidence among chemical workers in the United States and Western Europe. Epidemiology. 2001, 12: 727-40. 10.1097/00001648-200111000-00023.
International Agency for Research in Cancer: Occupational exposures in petroleum refining: crude oil and major petroleum fuels. 1989, IARC Monographs, Lion
Murunga E, Zawade ET: Environmental and occupational causes of toxic injury to the kidney and urinary tract. Environmental and Occupational Medicine. Edited by: Rom WN, Markowitz SB. 2007, Lippincott Williams & Wilkins, New York, 800-12. 4
Wrensch M, Minn Y, Chew T, Bondy M, Berger MS: Epidemiology of primary brain tumors: current concepts and review of the literature. Neuro-Oncology. 2002, 4: 278-99.
Yu CL, Wang SF, Pan PC, Wu MT, Ho CK, Smith TJ, Li Y, Pothier LJ, Christiani DC: No Association Between Residential Exposure to Petrochemicals and Brain Tumor Risk. Cancer Epidemiol Biomarkers Prev. 2005, 14: 3007-9. 10.1158/1055-9965.EPI-05-0594.
Liu CC, Chen CC, Wu TN, Yang CY: Association of Brain Cancer with Residential Exposure to Petrochemical Air Pollution in Taiwan Brain Cancer and Petrochemical Air Pollution. J Toxicol Environ Health. 2008, Part A 71: 310-314.
Janković S, Radosavljević V: Risk factors for bladder cancer. Tumori. 2007, 93: 4-12.
Negri E, La Vecchia C: Epidemiology and prevention of bladder cancer. Eur J Cancer Prev. 2001, 10: 7-14. 10.1097/00008469-200102000-00002.
WHO: World Health Organization Classification of Tumours. Pathology and Genetics of Tumours of the Urinary System and Male Genital Organs. Edited by: Eble JN, Sauter G, Epstein JI, Sesterhenn IA. 2004, IARC Press, Lyon
Baena AV, Allam MF, Díaz-Molina C, Del Castillo AS, Abdel-Rahman AG, Navajas RF: Urinary bladder cancer and the petroleum industry: a quantitative review. Eur J Cancer Prev. 2006, 15: 493-7. 10.1097/01.cej.0000198895.94560.38.
The authors wish to thank Michela Baccini for her useful suggestions in specifying the multiple imputation strategy.
The study was funded by WHO, European Centre for Environment and Health, project “Technical assistance to the Sicilian Region, special office for risk areas, to develop plans for environmental remediation”; National Institute for Occupational Safety and Prevention and Italian Ministry of Health, the “National Epidemiological Surveillance System for the Prevention of Occupational Cancer Cases”.
Roberta Pirastu and Pietro Comba served as technical consultants for the Gela public prosecutor in the first phase of the study.
RP contributed to all phases of the study: design, data collection, data analysis, manuscript preparation; AZ contributed to design, data collection and management of the morbidity study, and in manuscript preparation; RPirastu contributed to design and data collection of the mortality study, and in manuscript preparation; AC contributed to design and data collection of the morbidity study; GD contributed in discussing study design; SPA contributed in the morbidity data collection and management; SS contributed in discussing study design; PC contributed in study design and manuscript preparation. All authors read and approved the final manuscript.