Article Text

PDF

Original article
Organophosphate insecticide use and cancer incidence among spouses of pesticide applicators in the Agricultural Health Study
  1. Catherine C Lerro1,2,
  2. Stella Koutros1,
  3. Gabriella Andreotti1,
  4. Melissa C Friesen1,
  5. Michael C Alavanja1,
  6. Aaron Blair1,
  7. Jane A Hoppin3,
  8. Dale P Sandler4,
  9. Jay H Lubin5,
  10. Xiaomei Ma2,
  11. Yawei Zhang2,
  12. Laura E Beane Freeman1
  1. 1Occupational and Environmental Epidemiology Branch, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, Maryland, USA
  2. 2Yale University School of Public Health, New Haven, Connecticut, USA
  3. 3Department of Biological Sciences, North Carolina State University, Raleigh, North Carolina, USA
  4. 4Epidemiology Branch, National Institute of Environmental Health Sciences, Research Triangle Park, North Carolina, USA
  5. 5Biostatistics Branch, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, Maryland, USA
  1. Correspondence to Catherine C Lerro, Occupational and Environmental Epidemiology Branch, Division of Cancer Epidemiology and Genetics, National Cancer Institute, 9609 Medical Center Drive, 6E624, MSC 7991, Bethesda, MD 20892-7991, USA; lerrocc{at}mail.nih.gov

Abstract

Objectives Organophosphates (OPs) are among the most commonly used insecticides. OPs have been linked to cancer risk in some epidemiological studies, which have been largely conducted in predominantly male populations. We evaluated personal use of specific OPs and cancer incidence among female spouses of pesticide applicators in the prospective Agricultural Health Study cohort.

Methods At enrolment (1993–1997), spouses provided information about ever use of specific pesticides, including 10 OPs, demographic information, reproductive health history and other potential confounders. We used Poisson regression to estimate relative risks (RRs) and 95% CIs for all cancers diagnosed through 2010 for North Carolina and through 2011 for Iowa.

Results Among 30 003 women, 25.9% reported OP use, and 718 OP-exposed women were diagnosed with cancer during the follow-up period. Any OP use was associated with an elevated risk of breast cancer (RR=1.20, 95% CI 1.01 to 1.43). Malathion, the most commonly reported OP, was associated with increased risk of thyroid cancer (RR=2.04, 95% CI 1.14 to 3.63) and decreased risk of non-Hodgkin lymphoma (RR=0.64, 95% CI 0.41 to 0.99). Diazinon use was associated with ovarian cancer (RR=1.87, 95% CI 1.02 to 3.43).

Conclusions We observed increased risk with OP use for several hormonally-related cancers, including breast, thyroid and ovary, suggesting potential for hormonally-mediated effects. This study represents the first comprehensive analysis of OP use and cancer risk among women, and thus demonstrates a need for further evaluation.

Statistics from Altmetric.com

What this paper adds

  • Organophosphates are among the most commonly used insecticides.

  • Though organophosphates have been associated with increased cancer risk, there have been no prospective studies examining use of individual organophosphate insecticides and risk of multiple cancer sites in women.

  • We observed increased risk with organophosphate insecticide use for several hormonally-mediated cancers, including breast, thyroid and ovary.

  • Our results suggest the potential for hormonally-mediated effects of organophosphate insecticides with respect to cancer risk among women.

Background

Organophosphates (OPs) are among the most commonly sold and used active insecticide ingredients in the USA in all market sectors (ie, agriculture, home and garden, industrial, commercial and government), and currently comprise approximately 35% of insecticides used.1 Some OPs, such as malathion, are registered for outdoor residential use in the USA,2 while others, such as diazinon and chlorpyrifos, were once registered for residential use, but now only agricultural use is allowed.3 ,4 Selected OPs are used widely in the USA and abroad in public health programmes for mosquito control.5 The International Agency for Research on Cancer classifies malathion and diazinon as probably carcinogenic to humans (group 2A) and dichlorvos, parathion, and tetrachlorvinphos as possibly carcinogenic to humans (group 2B),6 with the US Environmental Protection Agency additionally classifying parathion as a possible human carcinogen.7

Increased cancer risk has been associated with several OP insecticides in epidemiological studies, including case-control studies in the USA,8 Canada9 and Italy,10 nested case-control studies of structural pest control workers in Florida11 and farm workers in California,12 and, more recently, among licensed pesticide applicators in the prospective Agricultural Health Study (AHS) cohort. AHS investigators have linked the use of diazinon, chlorpyrifos and terbufos to lung cancer,13–15 use of diazinon, terbufos, fonofos and malathion to leukaemia,14–17 and use of terbufos to non-Hodgkin lymphoma (NHL) overall, as well as to specific NHL subtypes.18 Additionally, increases in aggressive prostate cancer have been observed among male applicators who were using terbufos, fonofos and malathion.19 Many studies have focused on occupational exposure among farmers; however, OP insecticides are also widely used by others occupationally engaged in pest control, as well as residentially in the general population.

Studies of OP use and cancer outcomes have largely been conducted in predominantly male populations. Consequently, little is known about the potential impact of personal OP use among women, specifically on the development of female cancers, despite the fact that OPs as a class are thought to have endocrine disrupting properties.20–22 Moreover, many of the cancer sites to be examined, including breast, lung, ovary, uterus and thyroid, are of major public health importance in the USA, because they are commonly diagnosed and are important contributors to cancer deaths among women.23 In this analysis, we plan to evaluate the association between self-reported personal use of OP insecticides among spouses of pesticide applicators and subsequent cancer risk.

Methods

Study population

The AHS cohort has been described elsewhere in detail.24 Briefly, 52 394 private pesticide applicators (mainly farmers) and 4916 commercial pesticide applicators were recruited and enrolled during 1993–1997 in Iowa and North Carolina when they obtained or renewed their licenses to apply restricted use pesticides. Private applicators who indicated at enrolment that they were married were asked to have their spouse complete a take-home enrolment questionnaire focusing on farm exposures and general health, and a second questionnaire focusing on reproductive health history. The 32 345 spouses of private applicators who responded to the enrolment questionnaire are the focus of this study.

Exposure assessment

Use of OP insecticides and other potential confounders were assessed at enrolment using the spouse questionnaire, available at http://aghealth.nih.gov/background/questionnaires.html. Questions about pesticide use for spouses of pesticide applicators were asked as follows: ‘During your lifetime, have you ever personally mixed or applied [pesticide]? (Include pesticides used for farm use, commercial application and personal use in your home or garden.)’. They were prompted for specific pesticides using the active ingredient name and one or more trade names. Chemicals were grouped on the questionnaire according to functional class (insecticides, fungicides, herbicides, etc). Spouses self-reported lifetime ever use of 50 pesticides, including 10 OP insecticides (chlorpyrifos, coumaphos, diazinon, dichlorvos, fonofos, malathion, parathion, phorate, terbufos, trichlorfon). If any of these OPs were reported, the spouse was considered exposed to OPs as a chemical class. If they reported no exposure to any OP they were considered unexposed. Otherwise, they were considered to be missing for exposure to the chemical class grouping.

Cancer follow-up

Incident cancer cases were ascertained via regular linkage with Iowa and North Carolina state cancer registries. Cancer site was classified according to the International Classification of Diseases for Oncology (third revision), and lymphoma subtypes were classified according to the original Surveillance, Epidemiology and End Results Lymphoma Subtype Recode. We analysed first primary cancers diagnosed from the date of enrolment interview through date of death, movement out of state, or last date of study follow-up (31 December 2011 for Iowa, 31 December 2010 for North Carolina, USA), whichever was earliest. Our analysis included cancers with malignant behaviour, as well as in situ bladder cancers, which were included in the analysis, as per the standard grouping for bladder cancer. The study protocol was approved by all relevant institutional review boards.

Statistical analysis

We excluded male spouses as they were few (n=220), and women were the focus of our evaluation. We additionally excluded women with cancer diagnoses prior to enrolment (n=907), those with missing or zero person-years of follow-up (n=110) and women with missing information for all 10 OPs (n=1105), leaving 30 003 female spouses available for analysis. We excluded persons missing information for the OP of interest for specific analyses. For analyses of ovarian and uterine cancer, women were censored at date of oophorectomy or hysterectomy, or excluded if they had an oophorectomy (n=3074) or hysterectomy (n=5208) prior to study enrolment.

Relative risks (RRs) and 95% CIs were estimated using Poisson regression in SAS V.9.3 (SAS Institute, Inc, Cary, North Carolina, USA) for all cancer sites combined and specific sites where sample size allowed (n≥10 exposed cases). For the evaluation of use of any OPs as a class, no OP use was the referent category. For individual chemical analyses, persons reporting no use of the specific OP were included in the referent category. We conducted sensitivity analyses comparing those who applied an individual OP to those who never applied any OP (referent).

We adjusted all models for age (continuous), state of residence (Iowa or North Carolina), cigarette pack-years smoked as reported at enrolment (never smoker, pack-year quartiles: ≤1.5, 1.51–6.625, 6.626–18, >18, missing), race (white, other, missing), alcohol use (never, less than once per month, one to three times per month, once per week or more, missing), educational attainment (high school degree or less, some college, college graduate, missing), body mass index (BMI) (≤25, 25.1–30, >30, missing) and family history of cancer (yes, no, missing; specific cancer site where available). We also controlled for being the person who usually treats the home or lawn for pests, for ever use of specific pesticides most highly correlated with OP use (Spearman ρ>0.40; online supplementary table S1) and pesticides previously found to be associated with specific cancer outcomes in the AHS.25 In analyses of cancers of the breast, ovary and uterus, as well as all sites combined, we additionally adjusted for menopausal status at enrolment (no, before age 50 years, after age 50 years, missing), number of live births (0, 1, 2, 3, 4+, missing), and ever use of oral contraceptives at enrolment (yes, no, missing). We additionally explored inclusion of number of live births prior to the age of 30 years, and the use of hormone replacement therapy among postmenopausal women. These variables did not appreciably alter our results and thus were not included in the final models. We considered family history of cancer as a potential effect modifier, but interaction terms did not reach statistical significance in any model. We also considered adjusting for smoking using other metrics (eg, smoking duration in years, current/former/never use), but the results were similar to adjustment for pack-years smoked.

Breast cancers were examined by estrogen receptor (ER) and progesterone receptor (PR) status. For female cancer sites (breast, ovary, and uterus), we examined the statistical interaction between OP use and menopausal status at enrolment, and additionally performed stratified analyses by menopausal status at enrolment. We conducted sensitivity analyses restricting to cases diagnosed more than 5 years after enrolment, and also restricted to women who reported any pesticide application. We also performed analyses in which we did not control for home and lawn use, to ensure we were not over-adjusting for OP and correlated pesticide use. Finally, we stratified results by BMI (≤25, >25) to examine whether BMI might modify associations. All tests were two sided with α=0.05.

Results

Median follow-up time was 15.3 years. At enrolment, 25.9% of female spouses with valid information on OP use reported ever using at least one OP insecticide (online supplementary table S1). The most commonly used OP insecticides were malathion (19.5%) and diazinon (10.3%). Table 1 describes selected demographic, health and behavioural characteristics of AHS spouses. Ever users of OPs were older, from Iowa, white, more highly educated, heavier users of alcohol and more overweight than non-users. They were also more likely to have had a family history of cancer, more live births and to have gone through menopause at enrolment. Ever users of OPs were also more likely to report being the one who usually treats the home and/or lawn for pests. Lawn and home pesticide users were more likely to report herbicide use (data not shown).

Table 1

Selected characteristics of AHS spouses at enrolment with valid OP use information (n=29 325*), stratified by ever use of OP insecticides

Table 2 summarises the results for ever use of OPs and risk of cancers with n≥10 exposed cases. Use of any OP (RR=1.20, 95% CI 1.01 to 1.43) was significantly associated with breast cancer. Chlorpyrifos use (RR=1.41, 95% CI 1.00 to 1.99) and terbufos use (RR=1.52, 95% CI 0.97 to 2.36) were associated with non-significantly elevated risk of breast cancer. Chlorpyrifos was associated with a significantly increased risk of ERPR breast cancer (RR=2.26, 95% CI 1.07 to 4.75). Malathion use was associated with a significantly increased risk of thyroid cancer (RR=2.04, 95% CI 1.14 to 3.63) and decreased risk of NHL (RR=0.64, 95% CI 0.41 to 0.99). Diazinon use was associated with a significantly increased risk of ovarian cancer (RR=1.87, 95% CI 1.02 to 3.43). We observed no other associations between overall or specific OP use and cancer risk for any other site.

Table 2

RRs and 95% CIs* for ever use of OP insecticides, compared to never use, for various cancers among AHS spouses

We stratified analyses for cancers of the breast, ovary and uterus based on self-reported menopausal status at enrolment (table 3), with 15 144 women classified as premenopausal and 12 216 as postmenopausal. Among postmenopausal women, we observed significantly elevated risk of breast cancer associated with use of any OP (RR=1.27, 95% CI 1.00 to 1.62), and non-significantly elevated breast cancer risk associated with chlorpyrifos (RR=1.53, 95% CI 0.96 to 2.44) and terbufos (RR=1.73, 95% CI 0.93 to 3.21). Among women who used diazinon, we observed significantly elevated risk of ovarian cancer among premenopausal women (RR=3.26, 95% CI 1.31 to 8.13), but not postmenopausal women (RR=1.18, 95% CI 0.46 to 3.03, Pinteraction=0.06). We observed significant interactions with menopausal status and malathion for ovarian cancer risk (Pinteraction=0.04), and with menopausal status and diazinon for uterine cancer risk (Pinteraction=0.03). The stratum-specific risk estimates were not statistically significant, but the RRs were elevated among premenopausal women.

Table 3

RRs and 95% CIs* for ever-use of OP insecticides stratified by menopausal status at enrolment, compared to never use among AHS spouses for selected cancer sites

When analyses were restricted to cancer cases (n≥10) diagnosed at least 5 years after study enrolment (n=29 244), the results were mostly unchanged, with a few exceptions (table 4). The association between diazinon and ovarian cancer was no longer significant (RR=1.88, 95% CI 0.93 to 3.78) but remained elevated (n=10 exposed cases). We noted a statistically significant association between any OP use and multiple myeloma (RR=3.00, 95% CI 1.08 to 8.34). Additionally, diazinon (RR=1.24, 95% CI 0.99 to 1.56), coumaphos (RR=1.64, 95% CI 0.98 to 2.74) and parathion (RR=1.72, 95% CI 0.99 to 2.99) were all associated with non-significantly elevated risk of breast cancer.

Table 4

RRs and 95% CIs* for ever use of OP insecticides, compared to never use, for various cancers diagnosed five or more years after study enrolment among AHS spouses

When we restricted our study population to spouses who reported any pesticide application (n=16 685), the results remained unchanged. Results were similar in sensitivity analyses in evaluation of individual OPs and using those who had never used any OP as the referent group. We also evaluated the impact of controlling for home and lawn pesticide use; the results were similar with and without control. Stratification by BMI revealed that the significant results in the models were more pronounced overall among normal weight women (BMI≤25), with the exception of the association with malathion use and NHL, which was stronger among women with BMI >25 (results not shown).

Discussion

This is the first study, to our knowledge, to prospectively evaluate use of OPs and cancer at multiple sites among women. It also provides the first epidemiological evaluation of many female cancers, such as of the ovary and uterus, with this important chemical pesticide class. We observed increased risk of several hormonally-related cancers, including thyroid, ovary and breast.

We observed a strong association between malathion use and thyroid cancer. A previous study of male AHS private applicators found that malathion was associated with an increased prevalence of hypothyroidism;26 however, a similar study among female AHS spouses found no association with malathion and hypothyroidism, hyperthyroidism or other thyroid disease.27 Although hypothyroidism and hyperthyroidism have been hypothesised to be associated with thyroid cancer risk, the evidence has been somewhat inconsistent.28 ,29 Low thyroid stimulating hormone levels are thought to be associated with future thyroid cancer risk;30 however, laboratory studies in rats found that malathion was associated with increased thyroid stimulating hormone secretion.31 In agricultural areas, nitrate exposure via diet and drinking water has been associated with thyroid cancer and hypothyroidism;32 we were not able to control for nitrate intake in our analyses.

Increased risk of ovarian cancer was associated with diazinon use, with a significantly increased risk among women who were premenopausal, but not postmenopausal, at enrolment. We also noted a significant interaction between menopausal status and malathion use for risk of ovarian cancer, with elevated risk among premenopausal women. However, this may be a chance finding, as there was no overall association between malathion and ovarian cancer. An excess of ovarian cancer has been reported among female pesticide applicators in the AHS,33 but the small number of female applicators precluded evaluation by specific pesticides. Diazinon has been shown to alter DNA methylation patterns in the promoter regions of several genes associated with cancer, and has been correlated with decreased DNA excision repair in vitro.34 ,35 Diazinon use has been associated with shortened relative telomere length in male AHS pesticide applicators.36 Diazinon has also been shown to exhibit estrogenic properties, and to have a genotoxic effect on human mucosal cells.37 ,38

Use of any OP, terbufos and chlorpyrifos was associated with non-significantly increased breast cancer risk for each, in our analyses. We also noted significantly increased risk associated with chlorpyrifos use and ERPR breast cancer. An AHS study with follow-up for breast cancers through 2000, with 309 cases, saw no significantly increased risk for personal use of any OP.39 Our analysis of the updated cohort included 1059 accrued female breast cancer cases. A small registry-based case–control study of Hispanic farm workers in California examined use of pesticides and risk of breast cancer (n=128 cases), finding no association with diazinon and a suggestion of an association with malathion but no exposure–response.40 Many laboratory studies have noted associations between OPs and breast cancer in vitro and in vivo. OPs, particularly malathion and parathion, have been shown to induce malignant transformation of breast cells,41 ,42 alter estrogen activity and ER transactivity,43 ,44 and upregulate genes associated with carcinogenesis, sometimes in combination with estrogen.45

We observed a statistically significant inverse association with NHL and malathion use. A recent AHS analysis of male applicators found null associations with malathion use and NHL risk overall.18 Additionally, some case–control studies have shown a positive association for malathion use and NHL.9 ,46 ,47 Adjustment for other variables shown to be associated with NHL risk in farming populations, including whether or not our study participants grew up on a farm, their self-reported history of physician diagnosed allergies and contact with farm animals, did not alter the relationships in our study. Given these conflicting results and a lack of a plausible biological mechanism, it is unclear whether our observed inverse association between malathion and NHL is real or a chance finding. In sensitivity analyses restricting our population to those diagnosed more than 5 years after study enrolment, we noted a positive association with use of any OP and multiple myeloma. No association was observed with any individual OP and multiple myeloma in a recent study of AHS applicators,18 though excesses of multiple myeloma have been noted among pesticide applicators and in farming populations.33 ,48 The findings in our study were based on only a few exposed cases; therefore, further evaluations are needed to confirm these results.

The mechanism of pesticidal action of OPs involves inhibition of acetylcholinesterase activity.49 Excess acetylcholine as a result of OP exposure may act on cervical sympathetic neuronal nicotinic receptors, and activation of these neurones can promote thyroid hormone secretion via release of norepinephrine from the interfollicular adrenergic nerve endings.50–53 However, the potential mechanism of carcinogenesis may be unrelated to the mechanism of pesticidal action. Hypothesised carcinogenic mechanisms of OPs include increased cellular proliferation,42 oxidative stress54–56 and immunotoxicity.57 Given our findings with several hormonally-related cancers, it is also of note that OPs are thought to have endocrine disrupting properties. OPs may influence sex steroid hormone homeostasis, causing alterations in the levels of circulating and bioavailable sex hormones,58–61 and potentially impacting cellular proliferation and risk for hormone-related cancers.62 We noted that the associations for hormonally-related cancers were strongest among women with BMI ≤25. While there is some evidence that exposure to endocrine disrupting chemicals may impact body size, the relationship is not clear, and there are issues surrounding timing of exposure and reverse causation that make interpretation of these studies difficult.63 There is little information about a relationship between OP insecticide use and body size. A previous analysis in AHS examined the potential modifying effect of pesticides on the BMI-cancer association. There was a significant positive association between BMI and breast cancer in postmenopausal women who did not use OPs, but no association with BMI among those who did; however, the test for interaction was not significant.64

Strengths of our study include the longitudinal design with regular linkage to population registries for cancer and mortality outcomes, and little or no loss-to-follow-up, as well as information about the use of specific pesticides. Many epidemiological studies examine OPs as a class because of a small number of exposed cases, or because exposure to individual active ingredients is not evaluated. Owing to the prospective design of the AHS, there is no risk of differential reporting of pesticide use based on cancer outcome; any non-differential recall bias would bias the results toward the null. Blair et al65 assessed reliability of the AHS questionnaire among pesticide applicators; the level of agreement for ever pesticide use is quite high, ranging from 70% to greater than 90%. Though this work was carried out among applicators, we believe the spouses’ responses are similarly reliable. Spouses were prompted in the survey to mark all pesticides ever applied in their lifetime. Pesticide active ingredient and common trade names were listed in the survey. Because many of these women grew up on farms (60%), it is likely they are familiar with regularly used pesticides. In focusing on spouses, we were able to examine cancer outcomes that are unique to (eg, ovary, uterus) or most common among women (eg, breast, thyroid).

A limitation of our analysis was sample size; only about one quarter of our sample reported OP use at enrolment. Owing to a small number of cases, we were unable to evaluate very rare cancer sites and may have limited power to evaluate cancer sites with low incidence in AHS. Although information was collected on known risk factors for female cancers (eg, menopausal status at enrolment, oral contraceptive use, parity), certain important details were either not provided (type of oral contraceptive or hormone replacement therapy), or available for only a portion of the cohort (time-varying menopausal status), and thus could not be assessed as potential confounders. Many spouses in our cohort applied more than one pesticide in their lifetime. We controlled for use of pesticides that were highly correlated with the OP of interest, as well as pesticides that had been associated with specific cancer sites in previous analyses, to minimise these possible sources of confounding by use of multiple pesticides. Because we examined the use of several OP insecticides and cancer outcomes, it is possible that the findings could be due to chance. We were only able to examine self-reported lifetime personal ever use, and had no information about duration or time period of use. The assumption that all exposures are equivalent may be incorrect, as patterns of OP use and chemical formulation may have changed over time. The inability to differentiate between high and low use may mask potential associations. We also only evaluated personal use of pesticides in this analysis and not to exposure from other sources. Based on how pesticide information was ascertained, we were not able to distinguish between occupational OP use on the farm versus residential indoor and outdoor uses. Many women reported being the person who applies pesticides to the home and lawn, but did not report personal use of specific pesticides or pesticides overall. We controlled for being the person applying home and lawn pesticides in order to capture this use. We were concerned about potential for over-adjustment for OP use, however, insecticides applied in the home at this time were primarily pyrethroids, and insecticides used on the lawn were primarily OPs.66 However, lawn use in our study reflected primarily herbicide and not insecticide use. Thus we feel confident we were not over-adjusting for OP use.

In the first study to prospectively examine use of OP insecticides and risk of multiple cancer sites among women, we observed associations with several cancer sites including thyroid, ovary and breast. Previous studies examining OP insecticide use and cancer have focused primarily on men, making ours a unique evaluation. The increased risks that we observed for hormonally-related cancers are consistent with the hypothesis that OPs might act as endocrine disruptors, although additional studies exploring this and other possible mechanisms are needed. Future studies should continue to consider use of individual OPs to fully understand their impact on cancer risk. Because of the ubiquitous use of OP insecticides in agricultural and in residential settings, future research should attempt to confirm these findings by assessing exposure–response trends, non-occupational environmental sources of OP exposure and hormonal changes in women exposed to OPs.

Acknowledgments

Data in this analysis are based on the AHS data releases P1REL201209.00, P2REL201209 and P3REL201209.

References

View Abstract
  • Supplementary Data

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

Footnotes

  • Funding This work was supported by the Intramural Research Program of the National Institutes of Health, the National Cancer Institute (Z01-CP010119) and the National Institute of Environmental Health Sciences (Z01-ES049030).

  • Competing interests None declared.

  • Ethics approval National Cancer Institute Special Studies Institutional Review Board, Westat Institutional Review Board and the University of Iowa Institutional Review Board-01.

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

Request permissions

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