Occupational Exposure to Pesticides and Lung Cancer Risk: A Propensity Score Analyses

Article information

Cancer Res Treat. 2022;54(1):130-139
Publication date (electronic) : 2021 March 31
doi : https://doi.org/10.4143/crt.2020.1106
1Division of Cancer Prevention and Early Detection, National Cancer Control Institute, National Cancer Center, Goyang, Korea
2Graduate School of Cancer Science and Policy and Division of Cancer Prevention, National Cancer Center, Goyang, Korea
Correspondence: Eun Young Park, Division of Cancer Prevention and Early Detection, National Cancer Control Institute, National Cancer Center, 323 Ilsan-ro, Ilsandong-gu, Goyang 10408, Korea, Tel: 82-31-920-2930, Fax: 82-31-920-2934, E-mail: goajoa@ncc.re.kr
Received 2020 October 23; Accepted 2021 March 30.

Abstract

Purpose

Occupational exposure to pesticides is thought to be associated with lung cancer, but studies have yielded conflicting results. We performed a propensity score (PS) based analyses to evaluate the relationship between occupational exposure to pesticides and lung cancer risk in the Korea National Cancer Center community-based cohort study (KNCCCS).

Materials and Methods

During the follow-up period, 123 incidental lung cancer cases were identified, of the 7,471 subjects in the final statistical analysis. Information about occupational exposure to pesticides and other factors was collected at enrollment (2003–2010). Cox proportional hazards regression analyses were conducted. Four PS-based approaches (i.e., matching, stratification, inverse probability-of-treatment weighting, and the use of the PS as a covariate) were adopted, and the results were compared. PS was obtained from the logistic regression model. Absolute standardized differences according to occupational exposure to pesticides were provided to evaluate the balance in baseline characteristics.

Results

In the Cox proportional hazards regression model, the hazard ratio (HR) for lung cancer according to occupational exposure to pesticides was 1.82 (95% confidence interval [CI], 1.11 to 2.98). With all the propensity score matching (PSM) methods, the HRs for lung cancer based on exposure to pesticides ranged from 1.65 (95% CI, 1.04 to 2.64) (continuous term with PSM) to 2.84 (95% CI, 1.81 to 4.46) (stratification by 5 strata of the PS). The results varied slightly based on the method used, but the direction and statistical significance remained the same.

Conclusion

Our results strengthen the evidence for an association between occupational exposure to pesticides and the risk of lung cancer.

Introduction

Pesticides, including herbicides, insecticides, fungicides, fumigants, and rodenticides, are important chemical agents widely used in the agricultural sector. Pesticides can, however, also be toxic to humans, and it has been estimated that approximately 250,000 people worldwide die of pesticide poisoning each year [1]. Long-term occupational exposure to pesticides has been associated with an increased risk of lung cancer. Lung cancer is the leading cause of cancer-related death worldwide, including in South Korea. According to the annual 2017 report on cause of death statistics from the Korea National Statistics Office (KNSO), lung cancer mortality increased from 28.7 per 100,000 people in 2006 to 35.1 per 100,000 people in 2016 [2].

Although smoking is a known major risk factor for lung cancer that has been identified in many studies, exposure to pesticides also contributes to the occurrence of lung cancer [3]. Increased lung cancer-related mortality among pesticide applicators has been reported [48], suggesting the possibility that exposure to pesticides may increase the risk of lung cancer among farmers [9]. In addition, previous cohort studies reported positive associations between exposure to pesticides and the occurrence of lung cancer [913].

Exposure to pesticides has caused lung tumors in rodent bioassays, but the epidemiologic data supporting an association between pesticides and lung cancer risk in humans are mixed [10,14]. Jones et al. [12] reported an increased lung cancer incidence among male pesticide applicators with high pesticide exposure in lifetime days (diazinon: rate ratio, 1.60; 95% confidence interval [CI], 1.11 to 2.31; p-trend=0.02). However, Silver et al. [13] found no association of the incidence of lung cancer with either lifetime days or intensity-weighted lifetime days of pesticide use. In addition, a follow-up study that was part of the Agricultural Health Study (AHS) provided additional evidence supporting a positive association between pesticide use and lung cancer risk [9].

Although there have been reports on the risk of lung cancer due to pesticide exposure in several studies, previous studies have used limited data with inadequate control of covariate adjustments. To reduce confounding by numerous patient characteristics, propensity score (PS) methods are increasingly being used in many studies as an alternative to conventional covariate adjustment to investigate the effects of risk factors and include stratification, matching, inverse probability weighting (IPW), and the use of the PS as a covariate in a conventional regression model [15,16]. Therefore, we examined the association between occupational exposure to pesticides and lung cancer incidence in a community-based cohort study and compared the results of various PS methods and conventional models that adjust for covariates.

Materials and Methods

1. Data source and study population

This study was performed using data from the Korea National Cancer Center community-based cohort study (KNCCCS), a community-based cohort study conducted from 1993 to 2010 in Korea. The methods of the study have been described in detail in the published study protocol [17]. All participants (16,304 men and women) were over 20 years old, with an average age at cohort entry as follows: 58.6±12.4 years for men (n=6,302) and 57.7±13.3 for women (n=10,002). Data were collected through a face-to-face interview by well-trained interviewers using a questionnaire, and health outcomes were also investigated. The questionnaires were administered by trained interviewers and focused on demographic characteristics, past medical history, family history of cancer, history of medication use, dietary habits, smoking and alcohol consumption habits, physical activity, occupational history, history of exposure to acupuncture and blood transfusion, history of exposure to pesticides and electromagnetic fields and reproductive history for women. Trained staff, including interviewers, were provided appropriate instructions and guidelines about recruiting participants, and they asked the participants for consent prior to participation in the study. All study participants were followed until 2017 through linkage with the Korean Central Cancer Registry for cancer incidence. All participants in this cohort were linked to mortality data from Statistics Korea up to 2017.

Fig. 1 shows the flow diagram of the selection of the study population. Of the 16,304 subjects in the KNCCCS, 9,448 men and women were eligible for measurement of pesticide exposure. After exclusion of 2,095 subjects who had a history of cancer before an entry in the KNCCCS (n=244), missing data on pesticide-related variables such as frequency of pesticide use and duration (n=1,311), smoking status (n=194), and other covariates such as alcohol consumption status, education, obesity, physical activity and job (n=346), we included 7,471 subjects in the final statistical analysis. For propensity score matching (PSM) analysis, 4,112 subjects were ultimately included after excluding 3,359 unmatched subjects.

Fig. 1

Flow diagram of the selection of the study population.

2. History of exposure to pesticides

Information on exposure to pesticides was collected between 2003 and 2010 by trained interviewers. Participants were asked to report all pesticides used in the year prior to the interview and the frequency and duration of pesticide use. The history of exposure to pesticides was assessed using the following four questions: ‘In your lifetime, have you ever used pesticides?’, ‘Have you worked in agriculture in the past?’, ‘Are you currently working in agriculture?’, and ‘Have you worked in agriculture?’. Occupational exposure to pesticides was defined as having experience using pesticides and working in agriculture in the present.

3. Definition of lung cancer

The cohort has been followed up annually to identify new cancer cases from the date of the baseline survey. It has been an exclusively register-based passive follow-up, and the participants have not been recontacted. The principal outcome variable was incident lung cancer, based on data from the National Cancer Registry. The Korean Ministry of Health and Welfare started a nationwide, hospital-based cancer registry in 1980 [18]. More than 180 hospitals are currently participating, and the data cover approximately 99% of new cancer cases in Korea. Cases were identified through a personal identification number and other usual identification variables, such as names and addresses. Cancer cases were classified according to the International Classification of Diseases, 10th edition (ICD-10), and lung cancer was coded as C33–34 [18,19].

4. Statistical analyses

We examined differences of baseline characteristics such as pesticide use, current job status, pesticide-related variables such as the frequency and duration of pesticide use, smoking status, age, study areas, year of entry into the cohort, education achievement, alcohol consumption status, physical activity, and obesity.

Follow-up started at enrollment and continued until a lung cancer diagnosis or censoring. Censoring occurred at the date of death or a diagnosis of any other cancers or the end of follow-up (December 31, 2017).

We conducted crude and multiple Cox proportional hazards regression, including the following covariates that are well-known risk factors for lung cancer: age, sex, smoking status (nonsmoker, < 30 pack-years, ≥ 30 pack-years), alcohol consumption (0 g/day, < 24 g/day, ≥ 24 g/day), obesity (body mass index [BMI]; < 25 kg/m2, ≥ 25 kg/m2), physical activity (< 5 days/wk, ≥ 5 days/wk), education level (illiterate, middle school or less, high school, and college or more), and study area (Sancheong, Changwon, Chuncheon, Chungju, and Haman).

For PS-based analyses, we calculated the PS. The PS was the probability that an individual would have had a history of occupational exposure to pesticides based on personal demographic and lifestyle information and was obtained from the logistic regression model adjusted for these variables: age, sex, smoking status, alcohol consumption, obesity, physical activity, education level, and study area as follows; First, we applied a 1:1 PSM technique, which is an eighth-digit-to-first-digit greedy matching method, as elucidated previously. Unmatched cases are eliminated to ensure the best matching. To evaluate the balance in baseline characteristics between those who had ever and never used pesticides after PSM, we determined the absolute values of the standardized differences. Differences less than 0.1 were considered negligible. Second, we constructed weighted logistic regression models that showed the adjusted effect of exposure to pesticides after stratification into five strata according to the PS quintiles. Additionally, a logistic regression model using the PS as a covariate (a continuous variable and a categorical variable by quintiles), which is similar to traditional regression analyses, was constructed. Finally, two major PS-weighting methods, which are standardization methods depending on the establishment of a standard population, were adopted. One such method involved inverse probability-of-treatment weighting (IPTW). This model considers the total study subjects as the standard population and uses weights of (1/PS) for those who have used pesticides and [1/(1–PS)] for those who have not used pesticides.

Finally, we applied the Cox proportional hazards regression to determine the effect of occupational exposure to pesticides on lung cancer risk, using the following models: crude and adjusted models, PSM, regression adjusted with the PS, and PS weighting. Standardized differences were calculated using the macro ‘tableone’ written in R statistical software and developed by Yoshida et al. [20].

Additionally, time-lag analyses after excluding subjects with a follow-up period of equal or less than 3 years were performed.

All statistical analyses were performed with SAS software ver. 9.3 (SAS Institute Inc., Cary, NC) and R ver. 3.4.3 (R Foundation for Statistical Computing, Vienna, Austria) and were two-sided with a significance level of p < 0.05.

Results

The median follow-up period for subjects (n=7,471) eligible for the Cox proportional hazards regression analysis was 10.43 years (interquartile range [IQR], 8.5 to 12.48). A total of 123 primary lung cancer cases had been diagnosed since KNCCCS enrollment, and the median number of years of follow-up was 5.9 (IQR, 2.73 to 9.05).

The baseline characteristics of the study subjects are presented in Table 1. The mean age of the study subjects was approximately 60 years. The study population had higher proportions of, pesticides user, farm workers, subjects with BMI less than 25 kg/m2, and with a low education levels. Most of women were never-smokers (92.53%) and never-drinkers (76.05%). In contrast, most of men were former- or current smokers (79.4%) and drinkers (75.76%) (Table 1).

Baseline characteristics of the study participants

Table 2 shows the associations between occupational pesticide exposure-related variables and lung cancer risk. After stratification by sex, the association between occupational exposure to pesticides and lung cancer risk was statistically significant for men (HR, 2.40; 95% CI, 1.19 to 4.82). In addition, duration of pesticide use equal and greater to 20 years was associated with lung cancer risk (HR, 2.43; 95% CI, 1.05 to 5.65). Average frequency of pesticide use was not associated with the risk of lung cancer. In contrast, there were no statistical associations between occupational pesticide exposure-related variables and lung cancer risk for women. The findings from time-lag analyses after excluding subjects with a follow-up period of equal or less than 3 years were also not much different from primary analyses (S1 Table).

Hazard ratios for lung cancer incidence according to history of occupational pesticide exposure, and smoking status

All matched samples tended toward the characteristics of the unmatched data. Fig. 2 compares the covariate balance values for matching, stratification, and IPW by using the absolute standardized difference between the groups that did and did not use pesticides. Without the use of PS methods, the covariate balance was insufficient for almost all variables. Matching and stratification produced an excellent balance for all variables. IPTW achieved satisfactory covariate balance, except in year of cohort entry, the study area, owing to borderline cases. Due to the lack of covariate balance, the results of the analyses that did not use PS methods may be considered unreliable.

Fig. 2

Comparison of the standardized mean differences by covariates with different propensity score-based methods according to occupational exposure to pesticides. Crude, whole dataset; Matched, 1:1 matched dataset; Strata5, dataset stratified into five strata; IPTW, inverse probability-of-treatment weighting dataset. This plot shows the standardized mean differences between study subjects who had been occupationally exposed to pesticides and those who had not; a value > 0.1 indicates the imbalance of a covariate.

A comparison of the associations between occupational exposure to pesticides and lung cancer risk assessed by PS-based methods is shown in Table 3. The HRs ranged from 1.65 (1.04, 2.64) (regression adjusted for PS, as continuous term) to 2.84 (1.81, 4.46) (stratification by 5 strata of the PS).

Association between occupational exposure to pesticides and lung cancer incidence determined by PS-based methods

Discussion

To the best of our knowledge, this is the first report in which occupational exposure to pesticides has been found to be associated with lung cancer incidence using various PS methods or conventional covariate adjustment. Using data on the history of exposure to pesticides and lung cancer in the KNCCCS, we found estimated HRs (CI) ranging from 1.82 (1.11–2.98) (multiple Cox proportional hazards regression model) to 2.84 (1.81–4.46) for the stratification by five PS strata and PS-based methods yielded stronger associations than those obtained from the Cox proportional hazards regression model. There are advantages and disadvantages to PSM. PSM provides reliable results with excellent covariate balance, stratification retains data from all study subjects, and IPTW is easy to implement. In contrast, PSM leads to the exclusion of unmatched cases to ensure the best matching, stratification performs poorly when there are few outcome events, and IPTW may provide imprecise estimates of the treatment effect [15]. The meaningful findings of this study are consistent, regardless of whether PSM was used. Thus, this study provides strong evidence for the relationship between occupational exposure to pesticides and lung cancer risk.

Several epidemiological studies have found associations between exposure to pesticides and lung cancer [810]. In a cumulative meta-analysis of several cohort studies, the pooled relative risk estimate for lung cancer due to exposure to pesticides was 1.18 (standardized mortality ratio, 1.03 to 1.35; p=0.014) [9]. In a small, nested case-control study of structural pesticide applicators in Florida, Pesatori et al. [8] observed suggestive positive associations between exposure to pesticides and lung cancer (odds ratio, 2.4; 95% CI, 1.0 to 5.9). The AHS reported positive associations between select pesticides and the risk of lung cancer. Similar results were replicated in later studies with the AHS cohort [10]. The AGRIculture and CANcer (AGRICAN) cohort study of French farmers also suggested associations between the risk of small-cell lung cancer and exposure to pesticides (HR, 2.38; 95% CI, 1.07 to 5.28) [11]. These significant trends in the risk of lung cancer were also observed in our adjusted results. However, a number of studies have shown nonsignificant associations [2123] or negative relationships [2427] between occupational exposure to pesticides and lung cancer. Although several studies controlled for some important risk factors for lung cancer, such as indoor/outdoor air pollutants [28,29], lifestyle and psychosocial factors [30], and genetic predisposition [31], when assessing the association between pesticide exposure and lung cancer, such factors have not been routinely taken into account. In addition, impurities or promoting agents in pesticide formulations, such as dioxin and dioxin-like contaminants in phenoxy herbicides [5,32], might have contributed to the significant association found between exposure to some pesticides and lung cancer [5,8,32].

We found an association between the use of pesticides and lung cancer in the present analysis, although we did not obtain detailed exposure information for each type of pesticide. In vivo studies have shown that pesticides induce biochemical changes, leading to the development of tumor cells in the lung [33]. Since aryl hydrocarbon receptor activation mediates the tumorigenicity of dioxin-like compounds, microRNA expression induced by pesticides regulates lung cancer promotion at a wide range of workplace sites [34]. There is also experimental mechanistic evidence that pesticides such as chlorpyrifos and carbofuran can induce oxidative stress and oxidative DNA damage [35,36] and that pesticides may be genotoxic [37].

This study has several limitations. The number of lung cancer patients who had been exposed to pesticides was small, which negatively affected the precision of the results and our ability to evaluate risk stratified by the histologic type of lung cancer. Second, we did not have detailed exposure information, including years of use, applications per year, and applications in a lifetime, for each pesticide. Third, we could not fully control the effect of the study area as a confounding variable, although it was included multiple models. We attempted to consider stratification analysis by study area. However, the sample size according to study area was insufficient to obtain statistical power for stratification analysis. Forth, the selection biases that are present in most observational studies may impair causal inference in this study. To improve the ability to detect a causal relationship between exposure to pesticides and lung cancer incidence, we performed PS-based analyses in a community-based cohort. A major advantage of using PS-based methods in observational studies is that selection biases can be minimized by balancing nonrandomized individuals’ data to reach the level of causality determined in randomized controlled trials. Finally, there is the possibility of information bias due to the use of a self-reported questionnaire. However, the relationship between exposure to pesticides and lung cancer incidence was consistent across several PS methods, including stratification, matching, IPW, and use of the PS as a covariate in a logistic regression model, and the relationship revealed was stronger than that demonstrated in previous studies.

We found that occupational exposure to pesticides is associated with lung cancer risk. Inadequate control of the effects of covariates may have masked these effects in earlier studies. Although this study did not measure specific pesticides, we did observe a strong effect of occupational exposure to pesticides on the development of lung cancer. Further evaluations using PS methods of the associations between specific exposure to pesticides and a detailed history of exposure to pesticides and lung cancer risk are needed.

Supplementary Information

Supplementary materials are available at Cancer Research and Treatment website (https://www.e-crt.org).

Notes

Ethical Statement

The study protocols were approved by the Institutional Review Board of the National Cancer Center of Korea (IRB number: NCC2016-0300; NCC2017-0217). Signed written informed consent was obtained from all participants.

Author Contributions

Conceived and designed the analysis: Park EY.

Collected the data: Park EY, Oh JK, Lim MK.

Contributed data or analysis tools: Park EY, Kim J, Park E.

Performed the analysis: Kim J, Park E.

Wrote the paper: Kim B.

Conflicts of Interest

Conflict of interest relevant to this article was not reported.

Acknowledgements

We appreciate our collaborators at the public health centers of Haman-gun, Changwon-si, Chungju-si, and Chuncheon-si and the Health Center and County Hospital of Sancheong-gun for their support with regard to the implementation of the field surveys and data collection.

This study has been financially supported by the National Cancer Center of Korea (grant number: NCC-1710240-2).

References

1. van der Plaat DA, de Jong K, de Vries M, van Diemen CC, Nedeljkovic I, Amin N, et al. Occupational exposure to pesticides is associated with differential DNA methylation. Occup Environ Med 2018;75:427–35.
2. Statistics Korea. Annual report on the cause of death statistics Daejeon: Statistics Korea; 2016.
3. World Health Organization. WHO handbook on indoor radon: a public health perspective Geneva: World Health Organization; 2009.
4. Barthel E. Increased risk of lung cancer in pesticide-exposed male agricultural workers. J Toxicol Environ Health 1981;8:1027–40.
5. Becher H, Flesch-Janys D, Kauppinen T, Kogevinas M, Steindorf K, Manz A, et al. Cancer mortality in German male workers exposed to phenoxy herbicides and dioxins. Cancer Causes Control 1996;7:312–21.
6. Blair A, Dosemeci M, Heineman EF. Cancer and other causes of death among male and female farmers from twenty-three states. Am J Ind Med 1993;23:729–42.
7. MacMahon B, Monson RR, Wang HH, Zheng TZ. A second follow-up of mortality in a cohort of pesticide applicators. J Occup Med 1988;30:429–32.
8. Pesatori AC, Sontag JM, Lubin JH, Consonni D, Blair A. Cohort mortality and nested case-control study of lung cancer among structural pest control workers in Florida (United States). Cancer Causes Control 1994;5:310–8.
9. Bonner MR, Freeman LE, Hoppin JA, Koutros S, Sandler DP, Lynch CF, et al. Occupational exposure to pesticides and the incidence of lung cancer in the agricultural health study. Environ Health Perspect 2017;125:544–51.
10. Alavanja MC, Dosemeci M, Samanic C, Lubin J, Lynch CF, Knott C, et al. Pesticides and lung cancer risk in the agricultural health study cohort. Am J Epidemiol 2004;160:876–85.
11. Boulanger M, Tual S, Lemarchand C, Guizard AV, Delafosse P, Marcotullio E, et al. Lung cancer risk and occupational exposures in crop farming: results from the AGRIculture and CANcer (AGRICAN) cohort. Occup Environ Med 2018;75:776–85.
12. Jones RR, Barone-Adesi F, Koutros S, Lerro CC, Blair A, Lubin J, et al. Incidence of solid tumours among pesticide applicators exposed to the organophosphate insecticide diazinon in the Agricultural Health Study: an updated analysis. Occup Environ Med 2015;72:496–503.
13. Silver SR, Bertke SJ, Hines CJ, Alavanja MC, Hoppin JA, Lubin JH, et al. Cancer incidence and metolachlor use in the Agricultural Health Study: an update. Int J Cancer 2015;137:2630–43.
14. US Environmental Protection Agency. Office of pesticide programs: list of chemicals evaluated for carcinogenic potential Washington, DC: Environmental Protection Agency; 2002.
15. Elze MC, Gregson J, Baber U, Williamson E, Sartori S, Mehran R, et al. Comparison of propensity score methods and covariate adjustment: evaluation in 4 cardiovascular studies. J Am Coll Cardiol 2017;69:345–57.
16. Heinze G, Juni P. An overview of the objectives of and the approaches to propensity score analyses. Eur Heart J 2011;32:1704–8.
17. Oh JK, Lim MK, Yun EH, Choi MH, Hong ST, Chang SH, et al. Cohort profile: community-based prospective cohort from the National Cancer Center, Korea. Int J Epidemiol 2017;46:e14.
18. Shin HR, Won YJ, Jung KW, Kong HJ, Yim SH, Lee JK, et al. Nationwide cancer incidence in Korea, 1999–2001: first result using the national cancer incidence database. Cancer Res Treat 2005;37:325–31.
19. Riaz SP, Horton M, Kang J, Mak V, Luchtenborg M, Moller H. Lung cancer incidence and survival in England: an analysis by socioeconomic deprivation and urbanization. J Thorac Oncol 2011;6:2005–10.
20. Yoshida K, Bartel A, Chipman JJ, Bohn J, McGowan LD, Barrett M, et al. Create ‘Table 1’ to describe baseline characteristics with or without propensity score weights [Internet] Vienna: The Comprehensive R Archive Network; 2020. [cited 2018 Apr 4]. Available from: https://cran.r-project.org/web/packages/tableone/tableone.pdf .
21. Christensen CH, Platz EA, Andreotti G, Blair A, Hoppin JA, Koutros S, et al. Coumaphos exposure and incident cancer among male participants in the Agricultural Health Study (AHS). Environ Health Perspect 2010;118:92–6.
22. Freeman LE, Rusiecki JA, Hoppin JA, Lubin JH, Koutros S, Andreotti G, et al. Atrazine and cancer incidence among pesticide applicators in the agricultural health study (1994–2007). Environ Health Perspect 2011;119:1253–9.
23. Kang D, Park SK, Beane-Freeman L, Lynch CF, Knott CE, Sandler DP, et al. Cancer incidence among pesticide applicators exposed to trifluralin in the Agricultural Health Study. Environ Res 2008;107:271–6.
24. Greenburg DL, Rusiecki J, Koutros S, Dosemeci M, Patel R, Hines CJ, et al. Cancer incidence among pesticide applicators exposed to captan in the Agricultural Health Study. Cancer Causes Control 2008;19:1401–7.
25. Lynch SM, Mahajan R, Beane Freeman LE, Hoppin JA, Alavanja MC. Cancer incidence among pesticide applicators exposed to butylate in the Agricultural Health Study (AHS). Environ Res 2009;109:860–8.
26. Mahajan R, Blair A, Lynch CF, Schroeder P, Hoppin JA, Sandler DP, et al. Fonofos exposure and cancer incidence in the agricultural health study. Environ Health Perspect 2006;114:1838–42.
27. Rusiecki JA, Patel R, Koutros S, Beane-Freeman L, Landgren O, Bonner MR, et al. Cancer incidence among pesticide applicators exposed to permethrin in the Agricultural Health Study. Environ Health Perspect 2009;117:581–6.
28. Enomoto M, Tierney WJ, Nozaki K. Risk of human health by particulate matter as a source of air pollution: comparison with tobacco smoking. J Toxicol Sci 2008;33:251–67.
29. Samet JM, Eradze GR. Radon and lung cancer risk: taking stock at the millenium. Environ Health Perspect 2000;108(Suppl 4):635–41.
30. Yang P. Lung cancer in never smokers. Semin Respir Crit Care Med 2011;32:10–21.
31. Brennan P, Hainaut P, Boffetta P. Genetics of lung-cancer susceptibility. Lancet Oncol 2011;12:399–408.
32. Kogevinas M, Becher H, Benn T, Bertazzi PA, Boffetta P, Bueno-de-Mesquita HB, et al. Cancer mortality in workers exposed to phenoxy herbicides, chlorophenols, and dioxins. An expanded and updated international cohort study. Am J Epidemiol 1997;145:1061–75.
33. Walker NJ, Crockett PW, Nyska A, Brix AE, Jokinen MP, Sells DM, et al. Dose-additive carcinogenicity of a defined mixture of “dioxin-like compounds”. Environ Health Perspect 2005;113:43–8.
34. Boffetta P, Mundt KA, Adami HO, Cole P, Mandel JS. TCDD and cancer: a critical review of epidemiologic studies. Crit Rev Toxicol 2011;41:622–36.
35. Ojha A, Srivastava N. In vitro studies on organophosphate pesticides induced oxidative DNA damage in rat lymphocytes. Mutat Res Genet Toxicol Environ Mutagen 2014;761:10–7.
36. Zafiropoulos A, Tsarouhas K, Tsitsimpikou C, Fragkiadaki P, Germanakis I, Tsardi M, et al. Cardiotoxicity in rabbits after a low-level exposure to diazinon, propoxur, and chlorpyrifos. Hum Exp Toxicol 2014;33:1241–52.
37. Mladinic M, Zeljezic D, Shaposhnikov SA, Collins AR. The use of FISH-comet to detect c-Myc and TP 53 damage in extended-term lymphocyte cultures treated with terbuthylazine and carbofuran. Toxicol Lett 2012;211:62–9.

Article information Continued

Fig. 1

Flow diagram of the selection of the study population.

Fig. 2

Comparison of the standardized mean differences by covariates with different propensity score-based methods according to occupational exposure to pesticides. Crude, whole dataset; Matched, 1:1 matched dataset; Strata5, dataset stratified into five strata; IPTW, inverse probability-of-treatment weighting dataset. This plot shows the standardized mean differences between study subjects who had been occupationally exposed to pesticides and those who had not; a value > 0.1 indicates the imbalance of a covariate.

Table 1

Baseline characteristics of the study participants

Characteristic Overall (n=7,471) Men (n=2,748) Women (n=4,723)
Cases of lung cancer 123 (1.7) 89 (3.2) 34 (0.7)
Incidence rate per 100,000 person years 162.0 333.7 69.1
Pesticide use
 No 2,553 (34.2) 627 (22.8) 1,926 (40.8)
 Yes 4,918 (65.8) 2,121 (77.2) 2,797 (59.2)
Current job status
 Non-agriculture 2,912 (39.0) 853 (31.0) 2,059 (43.6)
 Agriculture 4,559 (61.0) 1,895 (68.0) 2,664 (56.4)
Occupational exposure to pesticides
 No 3,508 (46.9) 984 (35.8) 2,524 (53.4)
 Yes 3,963 (53.1) 1,764 (64.2) 2,199 (46.6)
Average frequency of pesticide use per year 6.2±10.7 8.2±12.3 5.1±9.5
 < 5 4,465 (59.8) 1,384 (50.4) 3,081 (65.2)
 5–9 1,039 (13.9) 401 (14.6) 638 (13.5)
 10–19 1,342 (18.0) 611 (22.2) 731 (15.5)
 ≥ 20 625 (8.4) 352 (12.8) 273 (5.8)
Duration of pesticide use (yr) 17.8±17.9 20.9±17.5 16.0±17.9
 < 5 2,932 (39.2) 785 (28.6) 2,147 (45.4)
 5–9 335 (4.5) 142 (5.2) 193 (4.1)
 10–19 714 (9.6) 290 (10.5) 424 (9.0)
 ≥ 20 3,490 (46.7) 1,531 (55.7) 1,959 (41.5)
Smoking status (pack-years) 9.9±18.8 25.2±23.3 1.0±5.1
 Never smoker 4,936 (66.1) 566 (20.6) 4,370 (92.5)
 < 30 1,412 (18.9) 1,100 (40.0) 312 (6.6)
 ≥ 30 1,123 (15.0) 1,082 (39.4) 41 (0.9)
Age (yr) 60.0±11.1 60.3±10.9 59.9±11.3
 20–29 31 (0.4) 12 (0.4) 19 (0.4)
 30–39 355 (4.7) 105 (3.8) 250 (5.3)
 40–49 1,025 (13.7) 356 (12.0) 669 (14.2)
 50–59 1,753 (23.5) 704 (25.6) 1,049 (22.2)
 60–69 2,858 (38.3) 1,021 (37.2) 1,837 (38.9)
 ≥ 70 1,449 (19.4) 550 (20.0) 899 (19.0)
Study area
 Sancheong 4,007 (53.6) 1,511 (55.0) 2,496 (52.9)
 Changwon 1,527 (20.4) 582 (21.2) 945 (20.0)
 Chuncheon 477 (6.4) 109 (4.0) 368 (7.8)
 Chungju 964 (12.9) 348 (12.6) 616 (13.0)
 Haman 496 (6.7) 198 (7.2) 298 (6.3)
Year of cohort entry
 2003 1,167 (15.6) 377 (13.7) 790 (16.7)
 2004 1,673 (22.4) 636 (23.1) 1,037 (22.0)
 2005 971 (13.0) 362 (13.2) 609 (12.9)
 2006 1,201 (16.1) 477 (17.4) 724 (15.3)
 2007 837 (11.2) 300 (10.9) 537 (11.4)
 2008 685 (9.2) 283 (10.3) 402 (8.5)
 2009 470 (6.3) 161 (5.9) 309 (6.5)
 2010 467 (6.2) 152 (5.5) 315 (6.7)
Education levels
 None 1,765 (23.6) 265 (9.6) 1,500 (31.8)
 Elementary or middle school 4,289 (57.4) 1,695 (61.7) 2,594 (54.9)
 High school 1,034 (13.9) 551 (20.1) 483 (10.2)
 College or more 383 (5.1) 237 (8.6) 146 (3.1)
Alcohol consumption status
 Never drinkers 4,258 (57.0) 666 (24.2) 3,592 (76.1)
 Former drinkers 461 (6.2) 322 (11.7) 139 (2.9)
 Current drinkers 2,752 (36.8) 1,760 (64.1) 992 (21.0)
Frequency of physical activity per week
 < 5 days 3,155 (42.2) 1,392 (50.7) 1,763 (37.3)
 ≥ 5 days 4,316 (57.8) 1,356 (49.3) 2,960 (62.7)
Obesity (BMI) 23.8±3.2 23.4±2.9 24.1±3.0
 No (BMI < 25 kg/m2) 4,968 (66.5) 1,967 (71.6) 3,001 (63.5)
 Yes (BMI ≥ 25 kg/m2) 2,503 (33.5) 781 (28.4) 1,722 (36.5)

Values are presented as number (%) or mean±SD. BMI, body mass index; SD, standard deviation.

Table 2

Hazard ratios for lung cancer incidence according to history of occupational pesticide exposure, and smoking status

Overall (n=7,471) Men (n=2,748) Women (n=4,723)



Cases (%) HR (95% CI) p-value Cases (%) HR (95% CI) p-value Cases (%) HR (95% CI) p-value
Model 1
 Occupational pesticide use

  No 31 (0.9) 1.000 (reference) 15 (1.5) 1.000 (reference) 16 (0.6) 1.000 (reference)

  Yes 92 (2.3) 3.125 (2.056–4.748) < 0.001 74 (4.2) 3.233 (1.825–5.727) < 0.001 18 (0.8) 1.542 (0.775–3.068) 0.218

 Average frequency of pesticide use per year

  Continuous 1.014 (1.003–1.026) 0.017 1.007 (0.992–1.022) 0.396 1.010 (0.982–1.039) 0.473

  < 5 61 (1.4) 1.000 (reference) 39 (2.8) 1.000 (reference) 22 (0.7) 1.000 (reference)

  5–9 13 (1.3) 0.966 (0.530–1.761) 0.910 11 (2.7) 1.024 (0.524–2.003) 0.945 2 (0.3) 0.467 (0.109–1.992) 0.304

  10–19 36 (2.7) 2.016 (1.332–3.052) 0.001 29 (4.8) 1.688 (1.041–2.737) 0.034 7 (1.0) 1.383 (0.588–3.254) 0.457

  ≥ 20 13 (2.1) 1.598 (0.876–2.913) 0.126 10 (2.8) 1.027 (0.511–2.061) 0.941 3 (1.1) 1.602 (0.478–5.373) 0.445

 Duration of pesticide use (yr)

  Continuous 1.030 (1.020–1.040) < 0.001 1.025 (1.013–1.037) < 0.001 1.028 (1.010–1.046) 0.003

  < 5 23 (0.8) 1.000 (reference) 11 (1.4) 1.000 (reference) 12 (0.6) 1.000 (reference)

  5–9 7 (2.1) 3.202 (1.361–7.535) 0.008 5 (3.5) 2.942 (1.007–8.596) 0.049 2 (1.0) 2.242 (0.496–10.127) 0.294

  10–19 10 (1.4) 2.117 (0.997–4.496) 0.051 9 (3.1) 2.642 (1.075–6.492) 0.034 1 (0.2) 0.497 (0.064–3.847) 0.503

  ≥ 20 83 (2.4) 3.738 (2.314–6.036) < 0.001 64 (4.2) 3.572 (1.839–6.938) < 0.001 19 (1.0) 2.172 (1.030–4.580) 0.042

 Smoking status (pack-years)

  Continuous 1.030 (1.025–1.034) < 0.001 1.023 (1.017–1.029) < 0.001 1.035 (1.001–1.071) 0.043

  Never smoker 27 (0.6) 1.000 (reference) 2 (0.4) 1.000 (reference) 25 (0.6) 1.000 (reference)

  < 30 32 (2.3) 4.452 (2.667–7.430) < 0.001 24 (2.2) 5.986 (1.414–25.333) 0.015 8 (2.6) 5.144 (2.318–11.419) < 0.001

  ≥ 30 64 (5.7) 11.924 (7.594–18.723) < 0.001 63 (5.8) 17.272 (4.226–70.587) < 0.001 1 (2.4) 4.582 (0.621–33.824) 0.130

Model 2

 Occupational pesticide use

  No 31 (0.9) 1.000 (reference) 15 (1.5) 1.000 (reference) 16 (0.6) 1.000 (reference)

  Yes 92 (2.3) 1.819 (1.109–2.984) 0.018 74 (4.2) 2.395 (1.191–4.817) 0.014 18 (0.8) 1.268 (0.570–2.820) 0.560

 Average frequency of pesticide use per year

  Continuous 0.998 (0.982–1.014) 0.814 0.995 (0.977–1.014) 0.626 1.009 (0.977–1.042) 0.593

  < 5 61 (1.4) 1.000 (reference) 39 (2.8) 1.000 (reference) 22 (0.7) 1.000 (reference)

  5–9 13 (1.3) 0.708 (0.379–1.321) 0.278 11 (2.7) 0.847 (0.420–1.708) 0.643 2 (0.3) 0.424 (0.095–1.884) 0.260

  10–19 36 (2.7) 1.241 (0.777–1.984) 0.366 29 (4.8) 1.279 (0.736–2.222) 0.383 7 (1.0) 1.155 (0.455–2.931) 0.762

  ≥ 20 13 (2.1) 0.799 (0.414–1.544) 0.505 10 (2.8) 0.694 (0.322–1.496) 0.351 3 (1.1) 1.329 (0.371–4.762) 0.662

 Duration of pesticide use (yr)

  Continuous 1.009 (0.998–1.020) 0.124 1.004 (0.990–1.017) 0.596 1.020 (1.001–1.040) 0.038

  < 5 23 (0.8) 1.000 (reference) 11 (1.4) 1.000 (reference) 12 (0.6) 1.000 (reference)

  5–9 7 (2.1) 2.690 (1.103–6.563) 0.030 5 (3.5) 3.894 (1.215–12.483) 0.022 2 (1.0) 2.075 (0.440–9.776) 0.356

  10–19 10 (1.4) 1.609 (0.721–3.590) 0.246 9 (3.1) 2.747 (0.984–7.671) 0.054 1 (0.2) 0.467 (0.059–3.719) 0.472

  ≥ 20 83 (2.4) 1.935 (1.093–3.425) 0.023 64 (4.2) 2.434 (1.048–5.654) 0.039 19 (1.0) 1.822 (0.789–4.208) 0.160

 Smoking status (pack-years)

  Continuous 1.020 (1.013–1.026) < 0.001 1.020 (1.013–1.026) < 0.001 1.025 (0.987–1.064) 0.197

  Never smoker 27 (0.6) 1.000 (reference) 2 (0.4) 1.000 (reference) 25 (0.6) 1.000 (reference)

  < 30 32 (2.3) 3.440 (1.833–6.454) < 0.001 24 (2.2) 6.533 (1.537–27.762) 0.011 8 (2.6) 3.782 (1.612–8.876) 0.002

  ≥ 30 64 (5.7) 7.447 (3.858–14.375) < 0.001 63 (5.8) 14.532 (3.541–59.638) < 0.001 1 (2.4) 3.429 (0.439–26.810) 0.240

Model 1: unadjusted. Model 2: Adjusted by age (continuous), sex, smoking status (pack-years, continuous), alcohol consumption, physical activity, obesity, education level, year of cohort entry and study area. CI, confidence interval; HR, hazard ratio.

Table 3

Association between occupational exposure to pesticides and lung cancer incidence determined by PS-based methods

PS methods Cases/Controls HR (95% CI) p-value
Matching for PS 36/4,112 1.749 (0.886–3.453) 0.107
Stratification into 5 strata by PS 116/6,813 2.838 (1.807–4.455) < 0.001
Regression adjusted for PS
 As a continuous term 123/7,348 1.652 (1.035–2.637) 0.035
 As quintiles 123/7,348 1.851 (1.126–3.041) 0.015
Weighted model
 IPTW 123/7,348 2.359 (1.737–3.204) < 0.001

CI, confidence interval; HR, hazard ratio; IPTW, inverse probability-of-treatment weighting; PS, propensity score.