Long-term exposure to urban air pollution and lung cancer mortality: A 12-year cohort study in Northern China

Long-term exposure to urban air pollution and lung cancer mortality: A 12-year cohort study in Northern China

STOTEN-20446; No of Pages 7 Science of the Total Environment xxx (2016) xxx–xxx Contents lists available at ScienceDirect Science of the Total Envir...

762KB Sizes 55 Downloads 62 Views

STOTEN-20446; No of Pages 7 Science of the Total Environment xxx (2016) xxx–xxx

Contents lists available at ScienceDirect

Science of the Total Environment journal homepage: www.elsevier.com/locate/scitotenv

Long-term exposure to urban air pollution and lung cancer mortality: A 12-year cohort study in Northern China Xi Chen a,b,1, Li-wen Zhang a,1, Jia-ju Huang a, Feng-ju Song c, Luo-ping Zhang b, Zheng-min Qian d, Edwin Trevathan e, Hong-jun Mao f, Bin Han g, Michael Vaughn d, Ke-xin Chen b, Ya-min Liu h, Jie Chen i, Bao-xin Zhao j, Guo-hong Jiang k, Qing Gu k, Zhi-peng Bai g, Guang-hui Dong l,⁎, Nai-jun Tang a,⁎⁎ a

Department of Occupational and Environmental Health, School of Public Health, Tianjin Medical University, Tianjin 300070, China Division of Environmental Health Sciences, School of Public Health, University of California, Berkeley, CA 94720-7360, USA c Departments of Epidemiology and Biostatistics, Key Laboratory of Cancer Prevention and Therapy, Tianjin, National Clinical Research Center of Cancer, Tianjin Medical University Cancer Institute and Hospital, Tianjin 300060, China d Department of Epidemiology, College for Public Health and Social Justice, Saint Louis University, Saint Louis, MO 63104, USA e Department of Psychology and Neuroscience, Baylor University, Waco, TX 76798-7334, USA f College of Environmental Science and Engineering, Nankai University, No. 94 Weijin Road, Nankai District, Tianjin 300071, China g State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, No. 8 Dayangfang, Chaoyang District, Beijing, China h Institute of Medicine and Health Information, Shandong Academy of Medical Sciences, Jinan, Shandong, China i Department of Occupational and Environmental Health, School of Public Health, China Medical University, No. 77 Puhe Road, Shenbei New District, 110122 Shenyang, Liaoning, China j Department of Chronic Disease Prevention, Taiyuan Center for Disease Control and Prevention, Taiyuan, Shanxi, China k Tianjin Center for Disease Control and Prevention, No. 6 Huayue Road, Hedong District, Tianjin 300011, China l Department of Preventive Medicine, School of Public Health, Sun Yat-sen University, Guangzhou 510080, China b

H I G H L I G H T S

G R A P H I C A L

A B S T R A C T

• The relationship between exposure to air pollution and lung cancer is proposed. • Levels of PM10 and SO2 were associated with lung cancer mortality independently. • The association was different in men and women, and varied across smoking status.

a r t i c l e

i n f o

Article history: Received 29 March 2016 Received in revised form 8 July 2016 Accepted 9 July 2016 Available online xxxx Editor: D. Barcelo

a b s t r a c t Cohort evidence that links long-term exposures to air pollution and mortality comes largely from the United States and European countries. We investigated the relationship between long-term exposures to particulate matter b 10 μm in diameter (PM10), nitrogen dioxide (NO2), and sulfur dioxide (SO2) and mortality of lung cancer in Northern China. A cohort of 39,054 participants were followed during 1998–2009. Annual average concentrations for PM10, NO2, and SO2 were determined based on data collected from central monitoring stations. Lung cancer deaths (n = 140) were obtained from death certificates, and hazard ratios (HRs) were estimated using Cox proportional hazards models, adjusting for age, gender, BMI, education, marital status, smoking status,

⁎ Corresponding author. ⁎⁎ Corresponding author at: Department of Occupational and Environmental Health, School of Public Health, Tianjin Medical University, No. 22 Qixiangtai Road, Heping District, Tianjin 300070, China. E-mail addresses: [email protected] (G. Dong), [email protected] (N. Tang). 1 These authors contributed equally to this work.

http://dx.doi.org/10.1016/j.scitotenv.2016.07.064 0048-9697/© 2016 Elsevier B.V. All rights reserved.

Please cite this article as: Chen, X., et al., Long-term exposure to urban air pollution and lung cancer mortality: A 12-year cohort study in Northern China, Sci Total Environ (2016), http://dx.doi.org/10.1016/j.scitotenv.2016.07.064

2

X. Chen et al. / Science of the Total Environment xxx (2016) xxx–xxx

Keywords: Air pollution Cohort study Lung cancer Mortality Particulate matter

passive smoking, occupation, alcohol consumption, etc. Each 10 mg/m3 increase in PM10 concentrations was associated with a 3.4%–6.0% increase in lung cancer mortality in the time-varying exposure model and a 4.0%–13.6% increase in the baseline exposure model. In multi-pollutant models, the magnitude of associations was attenuated, most strongly for PM10. The association was different in men and women, also varying across age categories and different smoking status. Substantial differences exist in the risk estimates for participants based on assignment method for air pollution exposure. © 2016 Elsevier B.V. All rights reserved.

1. Introduction The complex mixture of outdoor air pollutants contains a number of known carcinogens that have been found to be associated with increased risk of lung cancer in numerous studies over the past 50 years (Hamra et al., 2014). The International Agency for Research on Cancer (IARC) recently concluded that “exposure to outdoor air pollution and particulate matter (PM) in outdoor air is carcinogenic to humans (IARC Group 1) and causes lung cancer” (Hamra et al., 2014). Despite drawing on many of the latest studies conducted around the globe to reach this conclusion, few studies were conducted in China and other developing countries to demonstrate the association between outdoor air pollution and lung cancer deaths (Cao et al., 2011; Loomis et al., 2014). China did not officially monitor national-wide PM2.5 levels until Jan. of 2013. According to the Chinese Ministry of Environmental Protection, the annual average PM10 concentration in major Chinese cities increased from 87 μg/m3 in 2009 to 118 μg/m3 in 2012 (Report on the State of the Environment in China, 2009, 2012). This level of air pollution is much higher compared to levels in European and American countries. In the past three decades, China's extensive industrial development, coal-dependent energy consumption, and increasing number of vehicles have led to a dramatic rise in air pollutant emissions. The country is now experiencing two important challenges: 1) China produces a large number of major pollutants, which causes high levels of air pollution and substantially reduced visibility (Ouyang, 2013), and 2) air pollution over time has changed from the traditional coal combustion form to the mixed coal combustion/motor vehicle emission type (Kan et al., 2009). Lung cancer is now one of the major health threats experienced in China, and the burden is getting increasingly serious. In the annual report of cancer in China, lung cancer was listed as the number one cause of cancer incidence in males and the leading cause of cancer mortality for both males and females (National Office for Cancer Prevention and Control et al., 2014). In 2004–2005, the death rate from lung cancer was 30.84 per 100,000, representing a 464.8% increase over rates in 1973–1975 in China (Loomis et al., 2014). Positive exposure-response associations were reported in studies that collected quantitative pollutant exposure data (Krewski et al., 2009; Raaschou-Nielsen et al., 2013). In the previous analysis of this cohort, we evaluated the PM10 exposure–response relationship for cardiovascular mortality. We reported that long-term exposure to PM10 was independently associated with the mortality of cardiovascular disease after controlling for other risk factors (Zhang et al., 2014). In the present study, we explore the association of ambient exposures to PM10, sulfur dioxide (SO2), and nitrogen dioxide (NO2) and lung cancer mortality.

with a population of 12.3 million in 2009. Shenyang has an area of a total 13,308 km2 and had a population of 7.9 million as of 2009. In the city of Shenyang, the main industries include nonferrous metals, steel manufacturing, electric power generation machinery, and chemicaland coke-related industries. Taiyuan is the capital of Shanxi province and is the largest coal-producing area in China. The city is 800 m above sea level, and is located on the eastern side of the Loess Plateau. The population was 3.5 million in 2009. Rizhao, a merging coastal city, is located on the west coast of the Yellow China Sea and possessed a population of 2.7 million as of 2009. The number of environmental monitor stations in each city varied from one (Rizhao, the smallest city) to seven (Tianjin, the largest city). Small neighborhoods (street blocks or apartment buildings) within a 1 km radius from the monitoring stations were selected and numbered to form a sampling frame. Random samples were then drawn until the desired sample size was achieved. Each neighborhood unit contains households ranging from 500 to 700 persons. Eligible participants must have been born prior to January 1, 1975, and had to be a resident in the defined area since January 1, 1998. Among 48,114 participants, a total of 188 individuals (0.6%) were not linked to provincial administrative data maintained by the Institute for Clinical Evaluative Sciences (ICES) in Toronto. In addition, 2743 (0.6%) with missing residence location detail and 52 (0.2%) with incomplete data on covariates were excluded, yielding 39,054 participants (81.17%) in the study population. Participants completed a questionnaire during the visit, including demographic characteristics, smoking status, alcohol consumption, occupational exposures, and leisure exercise. The majority of the cohort (81.17%) did not change residences from January 1998 to the end of follow-up, whereas a move within the city was registered for the other 7353 residents. Written informed consent was obtained from all participants. Use of data on human subjects was approved by the California Health and Human Services Agency's Committee for the Protection of Human Subjects and the institutional review boards for each participating organization. 2.2. Outcome assessment

2. Methods

Family members completed the questionnaire for deceased participants. They were asked to provide information surrounding the death of the study participants including place, time, and cause of death. These data were crosschecked against death certificates obtained from the local Center for Disease Control and Prevention (CDC). More than 98% of all known causes of death were recorded in this survey. We classified deaths according to the cause of death using the International Classification of Disease-10 (ICD-10) (C33, C34 for lung cancer) coding system.

2.1. Study area and population

2.3. Estimation of exposure

As with the previous study, this cohort study was conducted from January 1998 to December 2009 in four Northern Chinese cities: Tianjin, Shenyang, Taiyuan, and Rizhao. These cities present a wide range of particle air pollution levels in northern China. The locations of the four cities and monitoring stations have been described in detail elsewhere (Zhang et al., 2014). Briefly, Tianjin comprises an area of 11,917.3 km2

Air pollutant data were acquired from the local Environmental Monitoring Centers, the government agencies in charge of the collection of data on air pollution. Measurement methods included tapered element oscillating microbalance, ultraviolet fluorescence, and chemiluminescence for PM10 (μg/m3), SO2 (μg/m3), and NO2 (μg/m3), respectively. We selected seven monitoring stations in Tianjin, five in

Please cite this article as: Chen, X., et al., Long-term exposure to urban air pollution and lung cancer mortality: A 12-year cohort study in Northern China, Sci Total Environ (2016), http://dx.doi.org/10.1016/j.scitotenv.2016.07.064

X. Chen et al. / Science of the Total Environment xxx (2016) xxx–xxx

Shenyang, two in Taiyuan, and one in Rizhao. The locations of these stations were sufficiently far from any source of emissions. The exposure levels of air pollution for the participants were estimated according to the data of the nearest monitoring stations. There was no professional ambient PM10 monitoring until 2000 in the city of Rizhao. Thus, the missing PM10 data for the first two years was converted from the concentrations of total suspended particles (TSP, the commonly used ratio is closer to 55 μg/m3 PM10 per 100 TSP (Ho and Nielsen, 2013)). Annual average concentrations were computed based on calculated 24-hour average concentrations. Previously, we estimated exposure for the study participants by taking the mean concentration of PM10 pollution over their surviving years during the cohort study (Zhang et al., 2014). This resulted in effect estimates that were biased high because pollution levels have decreased over time (as shown Fig. 1) and survivors were therefore assigned a lower exposure value, as opposed to participants who died during the study. This is not appropriate as had been observed previously in the California teacher study. For the purposes of this study, the baseline exposure model is most appropriate for evaluating lung cancer (Ostro et al., 2010; Ostro et al., 2011). Two methods were used to assign exposure to study subjects in the present analysis: 1) assigning the 1998 air pollution exposure for all and 2) time-varying exposures for both deaths and survivors. (How the time-varying exposure is defined should be described in detail.) We also added the HR estimates of lung cancer mortality and air pollution based on the average pollution for the surviving period for comparison.

2.4. Covariates Covariates we considered were: age (continuous variable), gender (male or female), marital status (single/separated/divorced/widowed or married), education (high school or more, less than high school), body mass index (BMI) (b20, 20–25, 25–30, N30 kg/m2), smoking status (with cigarette amounts per day, 0, 1–10, 11–20, N 20), alcohol consumption (yes or no), occupational exposures (self-reported exposure to dust and fumes in the workplace, yes or no) and leisure exercise (active or non-active).

3

2.5. Statistical analysis We used Cox proportional hazards models (PASW Statistics, Version 18.0, IBM) to analyze the associations between pollution concentrations and subsequent mortality in 1999–2009 (Beelen et al., 2008). We adjusted the analyses using age at baseline, gender, BMI, education, marital status, smoking status, passive smoking, alcohol consumption, occupation, diet, and leisure exercise. (not consistent with the variables listed in “Covariates”) To estimate the effects of NO2, SO2 and PM10, we used single-pollutant and bi-pollutant models. Because the average pollution of PM10 was highly correlated with SO2 (r = 0.947) from Pearson correlation analysis, the bi-pollutant models were used as PM10-andNO2 model, SO2-and-NO2 model. We performed stratified analyses to assess effect modification by the covariates. We modeled the Cox proportional hazards estimates described above separately for males and females to evaluate potential effect modification by sex. In addition, the model was stratified by age (b60 and ≥60). Furthermore, we conducted analyses that stratified cigarettes-per-day by cigarette amount to account for the impact of smoking on effect estimates and the shape of the exposure–response relationship. Because neighborhoods are usually inhabited by residents with similar characteristics (socioeconomic, health, diets) and similar environmental exposures, we also assessed confounding and clustering of the association between exposure and mortality by sensitivity analysis. We performed a shared frailty model to investigate the role of city as an additional sensitivity analysis. These models were calculated with the survival package of the Stata version 12.0 (StataCorp LP, College Station, TX) (Carey et al., 2013; Cesaroni et al., 2013). 3. Results 3.1. Description of the study subjects Table 1 shows the characteristics of the study participants from the four cities. The 39,054 members of the cohort were on average 44.29 years of age at the start of follow-up (SD = 13.95; range, 23– 89). The average BMI was 22.63. More than 25% (male 45.77%, female

Fig. 1. Average annual PM10, NO2 and SO2 concentrations (μg/m3) during the period of 1998 to 2008 in four cities in China. Seven monitoring stations in Tianjin, five in Shenyang, two in Taiyuan, and one in Rizhao were selected in the present study. Measurement methods included tapered element oscillating microbalance, ultraviolet fluorescence, and chemiluminescence for PM10, SO2, and NO2, respectively.

Please cite this article as: Chen, X., et al., Long-term exposure to urban air pollution and lung cancer mortality: A 12-year cohort study in Northern China, Sci Total Environ (2016), http://dx.doi.org/10.1016/j.scitotenv.2016.07.064

4

X. Chen et al. / Science of the Total Environment xxx (2016) xxx–xxx

Table 1 Baseline characteristics of the cohort participants from four cities. Characteristic

Tianjin

Shenyang

Taiyuan

Rizhao

Total

Total case Attrition rate Sex (%) Males Females Age (mean (SD)) BMI (mean (SD)) Income (%) b500 ≥500 Education (%) bHigh school ≥High school Smoking Amount (%) None 1–10 11–20 N20 Alcohol intake (%) Exercise (%) Occupational exposure (%) Marital status(%) Married Single/separated/divorced Total mortality (%) Lung cancer (%)

9663 20.80%

9921 17.50%

10,090 17.69%

9380 20.39%

39,054 18.83%

4661 (48.2) 5002 (51.8) 46.29 (13.00) 22.93 (3.19)

4811 (48.5) 5110 (51.5) 47.30 (14.41) 22.38 (3.36)

5051 (50.1) 5039 (49.9) 44.78 (14.00) 22.87 (2.66)

4937 (52.6) 4443 (47.4) 38.53 (12.56) 22.35 (2.52)

19,460 (49.8) 19,594 (50.2) 44.29 (13.95) 22.63 (2.97)

5325 (55.1) 4338 (44.9)

5959 (60.1) 3962 (39.9)

4894 (48.5) 5196 (51.5)

5555 (59.2) 3825 (40.8)

21,733 (55.6) 17,321 (44.4)

4642 (48.0) 4959 (51.3)

5966 (60.1) 3955 (39.9)

4926 (48.8) 5164 (51.2)

6631 (70.1) 2535 (27.0)

22,165 (56.8) 16,613 (42.5)

6565 (67.9) 1290 (13.3) 1394 (14.4) 219 (2.3) 2064 (21.4) 2584 (27.5) 712 (7.4)

6963 (70.2) 1390 (14.0) 1230 (12.4) 225 (2.3) 2027 (20.4) 4642 (48.0) 476 (4.8)

7386 (73.2) 1418 (14.1) 1111 (11.0) 142 (1.4) 1827 (18.1) 4153 (41.2) 1256 (12.4)

7483 (79.8) 1165 (12.4) 597 (6.4) 57 (0.6) 2021 (21.5) 2584 (27.5) 381 (4.1)

28,397 (72.7) 5263 (13.5) 4332 (11.1) 643 (1.6) 7939 (20.3) 14,946 (38.3) 2825 (7.2)

9001 (93.1) 632 (6.5) 177 (1.8) 18 (0.2)

8913 (89.8) 1008 (10.2) 489 (4.9) 42 (0.4)

9265 (91.8) 825 (8.2) 375 (3.7) 62 (0.6)

8757 (93.4) 623 (6.6) 312 (3.3) 18 (0.2)

35,936 (92.0) 3088 (7.9) 1353 (3.5) 140 (0.4)

The values are presented as percentages (%) or mean (SD). Some columns do not add up to 100% because of missing data. Attrition includes “moved out”, “nobody at home” and “refused to participate in the study” or “refused to complete the questionnaire”.

4.84%) were current smokers, and 20% (male 38.20%, female 4.60%) were drinkers. A total of 1353 (3.46%) participants died during the 12 years of follow-up (from 1998 to 2009, average length of followup: 11.83 years); of these, 140 (0.36%) died from lung cancer. As shown in Table 1, the four study sites were similar with respect to age, sex, smoking and drinking behavior. However, differences were evident for other characteristics (e.g. diet habits). The crude mortality rate (CMR) in the present study was 35 per 10,000 person-years. Lung cancer was responsible for 10.35% of all deaths (CMR = 4 per 10,000 person-years). 3.2. Pollution exposure We calculated annual average concentrations from the 24-hour average concentrations of PM10, SO2 and NO2 of all monitoring stations in each city (Fig. 1). The baseline average pollution levels [mean (SD, range, 25th percentile, 50th percentile, 75th percentile)] were 144.34 μg/m3 (3.63, 90.64–274.00, 125.50, 146.95, 183.88) for PM10, 40.66 μg/m3 (1.57, 18.00–78.59, 27.13, 38.18, 53.09) for NO2, and 66.90 μg/m3 (3.40, 11.00–224.40, 54.66, 70.64, 89.83) for SO2. The 1998 to 2009 average pollution of PM10 was highly correlated with SO2 (Correlation coefficient: 0.947, P = 0.000). 3.3. Association of air pollution and mortality Adjusted HRs (95% CIs) for lung cancer mortality in relation to air pollution (PM10, SO2 and NO2 concentrations) are presented in Tables 2 and 3. Two methods were used to assign exposure to study subjects: 1) assigning the 1998 air pollution exposure value for all and 2) timevarying exposures for both deaths and survivors. Additionally, average pollution for the surviving period was also used to assign exposure to subjects in Tables 2 and 3 to compare results with other two pollution-assign methods. There are overlapping CIs, the magnitudes are about the same, and there is the expected lessened precision in the time-varying model. These associations were greatly attenuated and very small after mutual adjustment.

In the fully adjusted single-pollutant models (PM10, SO2, and NO2 were modeled independently), each 10 μg/m3 increase in PM10 was associated with a 3.4%–6.0% increase in the risk of lung cancer death (HR per each 10 μg/m3 = 1.047; 95% CI: 1.034, 1.060) in time-varying exposure model, and with a 4.0%–13.6% increase in the risk of lung cancer death (HR per each 10 μg/m3 = 1.087; 95% CI: 1.040, 1.136) in 1998-exposure model. A significant positive association was also observed for SO2, with a 1.0%–2.5% increase in time-varying exposure model and a 0.8%–6.7% increase in 1998-exposure model. No association was observed between NO2 and mortality from all-cause disease and lung cancer in the fully adjusted single-pollutant model. The HRs (95% CI) of selected risk factors in the PM10 single-pollutant model were shown in Table S1 in the online data supplement. As the average pollution of PM10 was highly correlated with SO2, the concentrations of PM10 and SO2 were not included simultaneously in the bi-pollutant model. As shown in Table 3, PM10 and SO2 were statistically significant predictors of all-cause mortality and lung cancer in the bi-pollutant models.

3.4. Stratified analysis for PM10 and lung cancer mortality We evaluated the potential modification effects on the independent effect of PM10 with stratified analyses. Fig. 2 presents adjusted HRs (95% CIs) for lung cancer mortality in relation to mean PM10 concentrations stratified according to the selected characteristics at the baseline (gender, age, and daily cigarette amounts). There were some suggestions of modification effects by age (with b 60-year-old at higher risk than ≥ 60-year-old) and by gender (with men at higher risk than women) for lung cancer mortality. As reported in Fig. 2, the adjusted HRs were at higher levels in the population that smoked N 20 cigarettes/day. Among never-smokers, the adjusted HR of PM10 was higher in those exposed to second-hand smoke for lung cancer mortality than that in those not exposed to second hand smoke.

Please cite this article as: Chen, X., et al., Long-term exposure to urban air pollution and lung cancer mortality: A 12-year cohort study in Northern China, Sci Total Environ (2016), http://dx.doi.org/10.1016/j.scitotenv.2016.07.064

X. Chen et al. / Science of the Total Environment xxx (2016) xxx–xxx

3.5. Sensitivity analysis Confounding and clustering by city in the association between exposure and mortality was investigated in the sensitivity analysis, and the effect estimates in the shared frailty model were comparable to those obtained from the standard Cox model. 4. Discussion Using data from the four cities cohort in northern China, we found an association between increased lung cancer mortality and long-term exposure to urban air pollution. Each 10 μg/m3 increase in PM10 was associated with a significant 3.4%–6.0% increase in the risk of lung cancer death in time-varying exposure model. In addition, the average exposure of the cohort for 12 years (144.34 μg/m3 for PM10, 40.66 μg/m3 for NO2 and 66.90 μg/m3 for SO2) was higher than that in other study populations from Europe or North America. The mean concentrations of PM2.5 ranged from 9 to 21 μg/m3 in America during 1999–2009, and from 13.6 to 48.1 μg/m3 in Europe from 2008 to 2011 (Lepeule et al., 2012; Raaschou-Nielsen et al., 2013). In England in 2002, the mean concentrations of NO2 and SO2 were 22.5 and 3.9 μg/m3, respectively (Carey et al., 2013). The mean concentration of NO2 from 1988 to 2002 was 12.27 ppb (25.19 μg/m3) in North America (Jerrett et al., 2013). The atmosphere contains many compounds with known carcinogenic potential, such as polycyclic aromatic hydrocarbons and heavy metals. In such a complex mixture, it is difficult to ascribe the responsibility of causing lung cancer to a single component. Until a clear pathogenetic mechanism is defined, it would be prudent to consider the pollutants mentioned above as possible substitute markers of sources of air pollution, such as traffic or industrial pollution. PM can induce oxidative stress mediated by a particle-induced inflammation causing macrophages to release reactive oxygen species (ROS), transition metals on the particle surface capable of generating ROS through the Fenton reaction, or quinones in the particles that produce ROS through redox cycling. The variable concentrations of transition metals and organic compounds detected in the collected PM samples might be the active agents leading to a cumulative DNA damage, critical for carcinogenesis. NO2 has been shown to be a good indicator of traffic-related pollution, and SO2 can be used to estimate diesel and industrial emissions. It is important to note that large cohort studies on lung cancer and air pollution rely on measures of specific pollutants (Cui et al., 2015; Fajersztajn et al., 2013; Tseng et al., 2012). Previous cohort studies found associations that were statistically significant or of borderline significance (Beeson et al., 1998; Cao et al., 2011; Cesaroni et al., 2013; Cesaroni et al., 2012; Hales et al., 2012; Heinrich et al., 2013; Katanoda et al., 2011; Pope et al., 2002; Turner et al., 2011). Our study provides additional evidence to the current body of knowledge regarding the PM effects on lung cancer mortality. The present study, with 140 lung cancer cases, estimated an HR of 1.087 (95% CI: 1.040, 1.136) per 10 μg/m3 of PM10 in the single-pollutant model and 1.090 (95% CI: 1.044, 1.138) per 10 μg/m3 of PM10 in the

5

bi-pollutant model with NO2 using 1998 exposure for all. This finding is slightly lower than the estimate of an extended follow-up of the Harvard Six Cities study in a US cohort (351 cases) of 1.37 HR (95% CI: 1.07, 1.75) per 10 μg/m3, and the original and extended study of ACS CPS-II cohort (with a 15–27% increase in lung cancer mortality per 10 μg/ m3 PM2.5) (Lepeule et al., 2012; Turner et al., 2011). The result was also higher than the meta-analysis of 17 cohorts in nine European countries from the European Study of Cohorts for Air Pollution Effects (ESCAPE) (HR = 1.22 (95% CI: 1.03, 1.45) for PM10, HR = 1.18 (95% CI: 0.96, 1.46) for PM2.5, and HR = 1.09 (95% CI: 0.88, 1.33) for PMcoarse), from studies in Japan (HR = 1.24, 95% CI: 1.12, 1.37 with a 10-unit increase in PM2.5) and China (Cao et al., 2011; Katanoda et al., 2011; Raaschou-Nielsen et al., 2013). The HRs for lung cancer were lower in the present cohort, indicating that the dose-response curves were not linear, especially at higher concentrations than existing European Union and American air quality values for PM. The associations of lung cancer mortality with NO2 and SO2 were also estimated in this cohort. In the single-pollutant models, the effects of PM10 and SO2 exposure on lung cancer mortality were both statistically significant. In the bi-pollutant models, the effects of PM10 and SO2 were also significant on lung cancer mortality with the combined effect of NO2, with HRs of 1.051 (95% CI: 1.032, 1.069) for PM10 and 1.013 (95% CI: 1.006, 1.021) for SO2. Results from recent studies conducted in Rome and the Netherlands observed that the spatial distribution of air pollution is stable across the 10-year study periods (Cesaroni et al., 2012; Rosenlund et al., 2008). Since the latency period for lung cancer can be several decades and long-term exposure to air pollution is important in the development of lung cancer, high exposure levels of heating-related SO2 in the early years of the study period are especially important (Gulliver et al., 2011; Nicolich and Gamble, 2001). Our results reveal associations between air pollution and lung cancer in participants who did not move during the follow-up. Although the increasing number of vehicles has led to a rise in emissions of air pollutants, the effects of PM10, NO2 and SO2 in the present study suggest that China's substantial coal-dependent energy consumption still played an important role in health problems, since heating-related SO2, rather than traffic-related NO2, was consistently the stronger risk indicator. Identifying population subgroups susceptible to air pollution is an additional focus of this study. Some research has suggested that socioeconomic status, sex, and health characteristics, which are usually treated as confounders, could modify the association. In the present analysis, we identified the characteristics of participants who are most susceptible to the effects of air pollution by stratifying analyses (Cesaroni et al., 2013; Puett et al., 2009; Puett et al., 2008). Because active smoking and tobacco smoke are well-known risk factors for lung cancer, we also evaluated effect modification by daily amount of cigarette smoking and passive smoking. Among never-smokers, the adjusted relative risk (RR) for lung cancer mortality due to PM10 was higher in individuals exposed to second-hand smoke than in those who did not. Also, cigarette quantity caused risk ratios to fluctuate indicating that the quantity of cigarettes is

Table 2 Adjusteda HR (95%CI) for lung cancer mortality in relation to each 10 μg/m3 increase in air pollution according to selected risk factors, multiplicative scale, follow-up 1998–2009, China (single-pollutant models). Methods

Cause of death

No. of death

TVE

All-cause Lung cancer All-cause Lung cancer All-cause Lung cancer

1353 140 1353 140 1353 140

1998 E APSP

PM10 (10 μg/m3)

SO2(10 μg/m3)

NO2(10 μg/m3)

HR

95% CI

HR

95% CI

HR

95% CI

1.010 1.047 1.013 1.087 1.233 1.653

1.008, 1.012 1.034, 1.060 0.999, 1.027 1.040, 1.136 1.208, 1.260 1.527, 1.791

0.997 1.018 0.971 1.037 1.028 1.224

0.995, 0.999 1.010, 1.025 0.961, 0.981 1.008, 1.067 1.009, 1.047 1.164, 1.286

0.977 0.969 0.925 0.898 0.889 0. 837

0.972, 0.982 0.953, 0.984 0.899, 0.952 0.821, 0.982 0.857, 0.922 0.746, 0.938

TVE, time-varying exposures for both deaths and controls; 1998E, assigning the 1998 air pollution exposure value for all; APSP, average pollution for the surviving period was used to assign exposure to subjects. CI, confidence interval. HR, hazard ratio. PM10, particulate matter b10 μm in diameter; NO2, nitrogen dioxide; SO2, sulfur dioxide. a Adjusted for age at baseline, gender, BMI, education, marital status, smoking status, passive smoking, household income, alcohol consumption, occupation, diet, and leisure exercise.

Please cite this article as: Chen, X., et al., Long-term exposure to urban air pollution and lung cancer mortality: A 12-year cohort study in Northern China, Sci Total Environ (2016), http://dx.doi.org/10.1016/j.scitotenv.2016.07.064

6

X. Chen et al. / Science of the Total Environment xxx (2016) xxx–xxx

Table 3 Adjusteda HR (95%CI) for lung cancer mortality in relation to each 10 μg/m3 increase in air pollution according to selected risk factors, multiplicative scale, follow-up 1998–2009, China (bipollutant models). Methods

TVE 1998 E APSP

Cause of death

All-cause Lung cancer All-cause Lung cancer All-cause Lung cancer

No. of death

1353 140 1353 140 1353 140

Bi-pollutant model 1

Bi-pollutant model 2

PM10 (10 μg/m3)

NO2(10 μg/m3)

SO2(10 μg/m3)

NO2(10 μg/m3)

HR

95% CI

HR

95% CI

HR

95% CI

HR

95% CI

1.007 1.051 1.017 1.090 1.234 1.813

1.005, 1.009 1.032, 1.069 1.003, 1.031 1.044, 1.138 1.205, 1.264 1.657, 1.984

0.980 1.007 0.922 0.888 1.003 1.304

0.976, 0.985 0.985, 1.028 0.896, 0.948 0.812, 0.972 0.963, 1.044 1.148, 1.481

0.997 1.013 0.975 1.033 1.021 1.210

0.994, 0.999 1.006, 1.021 0.965, 0.985 1.005, 1.061 1.003, 1.038 1.145, 1.278

0.976 0.977 0.933 0.902 0.893 0.947

0.971, 0.981 0.961, 0.992 0.907, 0.960 0.823, 0.989 0.861, 0.926 0.842, 1.067

TVE, time-varying exposures for both deaths and controls; 1998E, assigning the 1998 air pollution exposure value for all; APSP, average pollution for the surviving period was used to assign exposure to subjects. CI, confidence interval. HR, hazard ratio. PM10, particulate matter b10 μm in diameter; NO2, nitrogen dioxide; SO2, sulfur dioxide. Bi-pollutant model 1 and 2: because the average pollution of PM10 was highly correlated with SO2, the bi-pollutant models were used as PM10-and-NO2 model (bi-pollutant model 1), SO2-and-NO2 model (bi-pollutant model 2), in turn. a Adjusted for age at baseline, gender, BMI, education, marital status, smoking status, passive smoking, household income, alcohol consumption, occupation, diet, and leisure exercise.

a modifying influence of the relationship between exposure to PM10 and lung cancer mortality (Beelen et al., 2008; Gallus et al., 2008; Pope et al., 2011). Our study has several limitations. First, individual-level exposures were estimated from participants' residence and nearest monitor data, which assumed that PM10 levels were relatively uniform in the surrounding areas. This limitation warrants the need for a more detailed exposure assessment (the land-use regression model) and the consideration of sources of exposure assessment (e.g. influence of indoor sources of air pollution). Second, owing to the nature of this retrospective cohort study, recall bias may be introduced for certain covariates, such as diet habits and physical activities. Additionally, data for histological subtypes of lung cancer was not obtained, which is a potential limitation of the study because the associations of air pollution and lung cancer were stronger for adenocarcinomas of the lung in some studies. For this reason the ICD-10 was used for classification instead of International Classification of Diseases for Oncology, 3rd edition [ICDO3].

5. Conclusion The association between lung cancer mortality and PM10 was pronounced in this extended follow-up study. Long-term exposure to SO2 was also associated with increased risk of lung cancer mortality in single-pollution models. The combined effect of NO2 and PM10 resulted in a significant increase in mortality risk for lung cancer. Significant differences exist in the risk estimates for participants based on assignment method for air pollution exposure. The present study results add to the limited literature of the long-term effects of high levels of air pollution on lung cancer in developing countries. Supplementary data to this article can be found online at http://dx. doi.org/10.1016/j.scitotenv.2016.07.064.

Acknowledgements This work was supported by the Special Environmental Research Fund for Public Welfare (No. 200709048) from Ministry of Environmental Protection of the People's Republic of China and the National Basic Research Program (973) of China (No. 2011CB503801).

References Beelen, R., Hoek, G., van den Brandt, P.A., Goldbohm, R.A., Fischer, P., Schouten, L.J., et al., 2008. Long-term effects of traffic-related air pollution on mortality in a Dutch cohort (NLCS-AIR study). Environ. Health Perspect. 116, 196–202. Beeson, W.L., Abbey, D.E., Knutsen, S.F., 1998. Long-term concentrations of ambient air pollutants and incident lung cancer in California adults: results from the AHSMOG study. Adventist Health Study on Smog. Environ. Health Perspect. 106, 813–822. Cao, J., Yang, C., Li, J., Chen, R., Chen, B., Gu, D., et al., 2011. Association between long-term exposure to outdoor air pollution and mortality in China: a cohort study. J. Hazard. Mater. 186, 1594–1600. Carey, I.M., Atkinson, R.W., Kent, A.J., van Staa, T., Cook, D.G., Anderson, H.R., 2013. Mortality associations with long-term exposure to outdoor air pollution in a national English cohort. Am. J. Respir. Crit. Care Med. 187, 1226–1233. Cesaroni, G., Porta, D., Badaloni, C., Stafoggia, M., Eeftens, M., Meliefste, K., et al., 2012. Nitrogen dioxide levels estimated from land use regression models several years apart and association with mortality in a large cohort study. Environ. Health 11, 48. Cesaroni, G., Badaloni, C., Gariazzo, C., Stafoggia, M., Sozzi, R., Davoli, M., et al., 2013. Longterm exposure to urban air pollution and mortality in a cohort of more than a million adults in Rome. Environ. Health Perspect. 121, 324–331. Cui, P., Huang, Y., Han, J., Song, F., Chen, K., 2015. Ambient particulate matter and lung cancer incidence and mortality: a meta-analysis of prospective studies. Eur. J. Pub. Health 25, 324–329. Fajersztajn, L., Veras, M., Barrozo, L.V., Saldiva, P., 2013. Air pollution: a potentially modifiable risk factor for lung cancer. Nat. Rev. Cancer 13, 674–678. Gallus, S., Negri, E., Boffetta, P., McLaughlin, J.K., Bosetti, C., La Vecchia, C., 2008. European studies on long-term exposure to ambient particulate matter and lung cancer. Eur. J. Cancer Prev. 17, 191–194. Gulliver, J., Morris, C., Lee, K., Vienneau, D., Briggs, D., Hansell, A., 2011. Land use regression modeling to estimate historic (1962–1991) concentrations of black smoke and sulfur dioxide for Great Britain. Environ. Sci. Technol. 45, 3526–3532.

Fig. 2. Adjusted hazard ratio (95% CIs) for lung cancer mortality per 10 μg/m3 elevation in PM10 concentrations, stratified by population characteristics of gender, age, quantity of cigarettes, and passive smoking in never smokers. The passive smoking population (Yes + No) was equal to the nonsmoking population.

Please cite this article as: Chen, X., et al., Long-term exposure to urban air pollution and lung cancer mortality: A 12-year cohort study in Northern China, Sci Total Environ (2016), http://dx.doi.org/10.1016/j.scitotenv.2016.07.064

X. Chen et al. / Science of the Total Environment xxx (2016) xxx–xxx Hales, S., Blakely, T., Woodward, A., 2012. Air pollution and mortality in New Zealand: cohort study. J. Epidemiol. Community Health 66, 468–473. Hamra, G.B., Guha, N., Cohen, A., Laden, F., Raaschou-Nielsen, O., Samet, J.M., et al., 2014. Outdoor particulate matter exposure and lung cancer: a systematic review and metaanalysis. Environ. Health Perspect. 122, 906–911. Heinrich, J., Thiering, E., Rzehak, P., Kramer, U., Hochadel, M., Rauchfuss, K.M., et al., 2013. Long-term exposure to NO2 and PM10 and all-cause and cause-specific mortality in a prospective cohort of women. Occup. Environ. Med. 70, 179–186. Ho, M.S., Nielsen, C.P., 2013. Clearing the Air: The Health and Economic Damages of Air Pollution in China (ISBN-13: 978-0262083584). Jerrett, M., Burnett, R.T., Beckerman, B.S., Turner, M.C., Krewski, D., Thurston, G., et al., 2013. Spatial analysis of air pollution and mortality in California. Am. J. Respir. Crit. Care Med. 188, 593–599. Kan, H., Chen, B., Hong, C., 2009. Health impact of outdoor air pollution in China: current knowledge and future research needs. Environ. Health Perspect. 117, A187. Katanoda, K., Sobue, T., Satoh, H., Tajima, K., Suzuki, T., Nakatsuka, H., et al., 2011. An association between long-term exposure to ambient air pollution and mortality from lung cancer and respiratory diseases in Japan. J. Epidemiol. 21, 132–143. Krewski, D., Jerrett, M., Burnett, R.T., Ma, R., Hughes, E., Shi, Y., et al., 2009. Extended follow-up and spatial analysis of the American Cancer Society study linking particulate air pollution and mortality. Res. Rep. Health Eff. Inst. 5–114 (discussion 115-36). Lepeule, J., Laden, F., Dockery, D., Schwartz, J., 2012. Chronic exposure to fine particles and mortality: an extended follow-up of the Harvard Six Cities study from 1974 to 2009. Environ. Health Perspect. 120, 965–970. Loomis, D., Huang, W., Chen, G., 2014. The International Agency for Research on Cancer (IARC) evaluation of the carcinogenicity of outdoor air pollution: focus on China. Chin. J. Cancer 33, 189–196. Ministry of Environmental Protection of the People's Republic of China, 2009l. Report on the State of the Environment in China. Ministry of Environmental Protection of the People's Republic of China, 2012l. Report on the State of the Environment in China. National Office for Cancer Prevention and Control, National Center for Cancer Registry, Disease Prevention and Control Bureau, National Health and Family Planning Commission of the People's Republic of China, 2014. Chinese Cancer Registry Annual Report (2013). Military Medical Science Press, Beijing. Nicolich, M.J., Gamble, J.F., 2001. Urban air pollution and lung cancer in Stockholm. Epidemiology 12, 590–592.

7

Ostro, B., Lipsett, M., Reynolds, P., Goldberg, D., Hertz, A., Garcia, C., et al., 2010. Long-term exposure to constituents of fine particulate air pollution and mortality: results from the California Teachers Study. Environ. Health Perspect. 118, 363–369. Ostro, B., Tobias, A., Querol, X., Alastuey, A., Amato, F., Pey, J., et al., 2011. The effects of particulate matter sources on daily mortality: a case-crossover study of Barcelona, Spain. Environ. Health Perspect. 119, 1781–1787. Ouyang, Y., 2013. China wakes up to the crisis of air pollution. Lancet Respir. Med. 1, 12. Pope III, C.A., Burnett, R.T., Thun, M.J., Calle, E.E., Krewski, D., Ito, K., et al., 2002. Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution. JAMA 287, 1132–1141. Pope III, C.A., Burnett, R.T., Turner, M.C., Cohen, A., Krewski, D., Jerrett, M., et al., 2011. Lung cancer and cardiovascular disease mortality associated with ambient air pollution and cigarette smoke: shape of the exposure-response relationships. Environ. Health Perspect. 119, 1616–1621. Puett, R.C., Schwartz, J., Hart, J.E., Yanosky, J.D., Speizer, F.E., Suh, H., et al., 2008. Chronic particulate exposure, mortality, and coronary heart disease in the nurses' health study. Am. J. Epidemiol. 168, 1161–1168. Puett, R.C., Hart, J.E., Yanosky, J.D., Paciorek, C., Schwartz, J., Suh, H., et al., 2009. Chronic fine and coarse particulate exposure, mortality, and coronary heart disease in the Nurses' Health Study. Environ. Health Perspect. 117, 1697–1701. Raaschou-Nielsen, O., Andersen, Z.J., Beelen, R., Samoli, E., Stafoggia, M., Weinmayr, G., et al., 2013. Air pollution and lung cancer incidence in 17 European cohorts: prospective analyses from the European Study of Cohorts for Air Pollution Effects (ESCAPE). Lancet Oncol. 14, 813–822. Rosenlund, M., Forastiere, F., Stafoggia, M., Porta, D., Perucci, M., Ranzi, A., et al., 2008. Comparison of regression models with land-use and emissions data to predict the spatial distribution of traffic-related air pollution in Rome. J. Expo. Sci. Environ. Epidemiol. 18, 192–199. Tseng, C.Y., Huang, Y.C., Su, S.Y., Huang, J.Y., Lai, C.H., Lung, C.C., et al., 2012. Cell type specificity of female lung cancer associated with sulfur dioxide from air pollutants in Taiwan: an ecological study. BMC Public Health 12, 4. Turner, M.C., Krewski, D., Pope 3rd, C.A., Chen, Y., Gapstur, S.M., Thun, M.J., 2011. Longterm ambient fine particulate matter air pollution and lung cancer in a large cohort of never-smokers. Am. J. Respir. Crit. Care Med. 184, 1374–1381. Zhang, L.W., Chen, X., Xue, X.D., Sun, M., Han, B., Li, C.P., et al., 2014. Long-term exposure to high particulate matter pollution and cardiovascular mortality: a 12-year cohort study in four cities in northern China. Environ. Int. 62, 41–47.

Please cite this article as: Chen, X., et al., Long-term exposure to urban air pollution and lung cancer mortality: A 12-year cohort study in Northern China, Sci Total Environ (2016), http://dx.doi.org/10.1016/j.scitotenv.2016.07.064