Article Text
Abstract
Objectives There is a general consensus that pesticides are involved in the aetiology of Parkinson's disease (PD), although associations between specific pesticides and the risk of developing PD have not been well studied. This study examines the risk of developing PD associated with specific organophosphate (OP) pesticides and their mechanisms of toxicity.
Methods This case–control study uses a geographic information system-based exposure assessment tool to estimate ambient exposure to 36 commonly used OPs from 1974 to 1999. All selected OPs were analysed individually and also in groups formed according to their presumed mechanisms of toxicity.
Results The study included 357 incident PD cases and 752 population controls living in the Central Valley of California. Ambient exposure to each OP evaluated separately increased the risk of developing PD. However, most participants were exposed to combinations of OPs rather than a single pesticide. Risk estimates for OPs grouped according to different presumed functionalities and toxicities were similar and did not allow us to distinguish between them. However, we observed exposure-response patterns with exposure to an increasing number of OPs.
Conclusions This study adds strong evidence that OPs are implicated in the aetiology of idiopathic PD. However, studies of OPs at low doses reflective of real-world ambient exposure are needed to determine the mechanisms of neurotoxicity.
Statistics from Altmetric.com
Introduction
Parkinson's disease (PD) is an idiopathic neurodegenerative disease associated with aging, environmental and genetic factors, and gender.1 ,2 Although many studies have found associations between pesticides and PD, much heterogeneity of findings remains to be explained.3 While a number of methodological limitations may contribute to conflicting reports in the literature, including difficulties in correctly ascertaining PD case status, inappropriate control selection and lack of statistical power, the major limitations of most studies are due to inadequate lifetime exposure assessment for pesticides. Previous studies generally assigned pesticide exposure based on self-report, which is likely affected by recall bias. In addition, many studies define occupational pesticide exposure as exposure to any type of pesticide or pesticides belonging to broad pesticide classes (ie, insecticides, fungicides, herbicides), further contributing to conflicting findings if some, but not all, pesticides in each class contribute to PD aetiology.3
Organophosphate (OP) pesticides represent the largest group of insecticides and are commonly used in farming in the USA despite being responsible for millions of poisonings and thousands of deaths worldwide.4–6 Some humans who experienced acute OP poisonings also developed signs of parkinsonism, suggesting that OPs may have an effect on the striatal dopaminergic system.7 ,8 While the main mechanism of OP toxicity is oxidative stress likely mediated by cholinesterase inhibition, there has recently been some suggestion that OPs may disrupt mitochondrial functions.5
This study aims to assess if specific OPs affect PD risk and whether OPs with presumed functional and/or neurotoxic mechanisms, such as teratogenicity or oxidative stress generation, contribute to PD risk. Since mitochondrial inhibition has been proposed as a prominent pathological pathway for PD aetiology, we assess the group of OPs presumably having this action in our human population.
Materials and methods
All of the research procedures described in this study were approved by the UCLA Internal Review Board for human subjects. Written informed consent was obtained from all participants.
Participant recruitment
Incident idiopathic PD patients were enrolled between 1 January 2001 and 1 January 2007, and population-based controls from the mostly rural agricultural tri-county area (Kern, Tulare, Fresno) in central California between 1January 2002 and early 2011. Subject recruitment methods9 ,10 and case definition criteria11 ,12 have been described in detail elsewhere.
Briefly, cases were defined as individuals who (1) had been diagnosed with PD for the first time by a physician within the past 3 years; (2) were residents of Fresno, Kern, or Tulare Counties and had lived in California for at least 5 years; (3) had been seen by UCLA movement disorder specialists and confirmed as having clinically ‘probable’ or ‘possible’ PD; (4) did not have any other diagnosed neurological condition or serious psychiatric condition; (5) were not in the last stages of a terminal illness; and (6) had agreed to participate in the study.
Of the 1167 PD patients initially identified through neurologists, large medical groups and public service announcements, 604 did not meet eligibility criteria (397 with an initial PD diagnosis more than 3 years prior to recruitment, 134 lived outside the tri-county area, 51 without idiopathic PD and 22 were too ill to participate). Of 563 eligible cases, 90 could not be examined by the study movement disorder specialists (56 declined to participate or moved away, 18 had become too ill and 16 died prior to their appointment). Of 473 patients examined, 107 did not meet eligibility criteria (94 did not meet published criteria for idiopathic PD13 ,14 when examined initially, an additional 13 were reclassified as non-PD during further follow-up).15 Additionally, six subjects withdrew before the exposure interview and three were excluded due to a PD diagnosis outside our study period, leaving a total of 357 idiopathic PD cases for this study.
Population-based controls were recruited initially from Medicare lists (2001) and, after the implementation of the Health Insurance Portability and Accountability Act (HIPAA), from residential tax assessor records for the tri-county area. Two sampling strategies, described in detail elsewhere,9 were implemented to increase enrolment success and achieve representativeness of the control population: (1) random selection of residential parcels and enrolment via mail and phone and (2) clustered random selection of five households and enrolment via in-person visits to homes.
Of the 1212 potential controls contacted through the first recruitment strategy, 457 were ineligible: 409 were <35 years of age, 44 were too ill to participate and 4 primarily resided outside the study area. Of the 755 eligible population controls, 409 declined enrolment in the study, became too ill to participate or moved out of the area after screening and prior to enrolment. A total of 346 population controls enrolled based on the first strategy. Of the 346 controls, 341 provided all information needed in this analysis. Of the 4756 individuals screened for eligibility through the second strategy, 3515 were ineligible (88% due to age criteria). From 1241 eligible population controls, 634 declined participation and 607 enrolled. We excluded 196 of the 607 enrolled controls from the analysis because 183 only completed an abbreviated questionnaire, 2 did not provide ethnicity information and 11 provided incomplete residential or occupational histories. There were 24 participants recruited at their homes that did not provide information on family history of PD and were assumed to not have a family history of PD. Additionally, more than one control was recruited from some neighbourhood clusters of five homes, which potentially may introduce biases due to correlated exposures if these controls had been neighbours between 1974 and 1999. However, after selecting one control at random and excluding all others from multiple control clusters (n=23), our results were comparable and thus we will report on all controls recruited under the second strategy. In total, 752 controls provided all information necessary for the following analyses.
For all study participants, we conducted telephone interviews to obtain demographic information, other risk factor data, as well as residential and workplace addresses.
GIS-based environmental pesticide exposure assessment
Using our Geographic Information System (GIS)-based system, we combined pesticide use reports (PUR) data, land-use maps and geocoded address information16 ,17 to produce estimates of pesticide exposure within a 500 m radius circular buffer around participants’ workplace and residential addresses as suggested in previous literature.18–20 Using a 500 m buffer also represents an intermediate distance for pesticide applications as some studies suggested a 200 m buffer while others found a 1250 m buffer to work best depending on the pesticide examined.21 ,22 Residential and workplace addresses were automatically geocoded to TigerLine files (Navteq, 2006) and then manually resolved in a multistep process similar to that described by McElroy.19 A specific address was considered precisely geocoded if it was resolved at the level of an actual address, a parcel/lot centroid, street centroid or street intersection. Less precisely geocoded addresses were resolved at the zip code, city, county, state centroids or did not have enough information to be geocoded. A technical discussion16 and a detailed description10 of our approach have been provided elsewhere.
Briefly, since 1974 the California Department of Pesticide Regulations has collected PURs for agricultural pesticide application. Each PUR record included the name of the pesticide's active ingredient, the poundage applied, the crop and acreage of the field, the application method and location and date of application. Our GIS-based computer model combines data from PURs with California Department of Water Resources (CDWR) land-use maps and historical workplace and residential addresses from study participants to estimate ambient pesticide exposure. Annual exposure estimates were calculated by adding the poundage of pesticide applied in each 500 m buffer surrounding the workplace or residential address and weighting the total poundage by the proportion of acreage treated within the buffer. For each of the pesticides examined in this study, we summed the annual pounds applied per acre to obtain 26 annual exposure values for each pesticide separately for workplace and residential addresses. These annual exposure estimates were then averaged across the 26-year study period from 1974 to 1999. We chose to include a pesticide from the OP class if five or more cases or controls were exposed to any amount of that pesticide and included a total of 36 OPs in this study (see online supplementary appendix).
For the present analysis, we considered a participant exposed to a specific OP pesticide if their 26-year average ambient exposure was equal to or greater than the medians observed in the controls at either residence or workplace addresses, respectively. Thus, ambient exposures to individual pesticides were calculated and determined separately at workplace and residential addresses. We were also able to determine if participants were exposed only at residential addresses, workplace addresses or at both residential and workplace addresses. It is important to emphasise that ambient exposures at workplace addresses are not necessarily associated with occupational exposures as a result of job tasks. Participants who work within 500 m of a field where a pesticide was applied would be considered ambiently exposed to that pesticide even if that participant's occupation was unrelated to agriculture and pesticide application. Study participants were considered unexposed to pesticides if (1) they did not work or live within 500 m of pesticide applications in the tri-county area between 1974 and 1999 and thus could not be assigned an exposure estimate or (2) were never exposed to at least the median value of the control group for any pesticide.
Statistical analysis
All selected OPs were analysed individually and also in groups according to their presumed mechanisms of toxicity, such as being a carcinogen (carcinogenic: known or possible carcinogen; non-carcinogenic: not likely or not listed as carcinogen), teratogen or endocrine disruptor, and their acute toxicity due to acetylcholinesterase inhibition (toxic: extreme or high toxicity; non-toxic: moderate or slight toxicity) based on the Pesticide Action Network Pesticide Database (see online supplementary appendix).23 We also grouped OPs according to some presumed mitochondrial disrupting potential based on some recent literature (see online supplementary appendix). However, these classifications by presumed mechanism were not mutually exclusive (ie, a pesticide may have been assigned to multiple groups).
For the analysis of the individual OPs (table 2) as well as for the groups of OPs with different mechanisms of toxicity (table 3), we created a variable for each OP or group of OPs with the same mechanism of toxicity individually, which were analysed in separate models. We took a different approach when we assessed the association between mitochondria disrupting OPs and PD risk (table 4). We created and analysed a single variable that categorised all the participants into mutually exclusive groups: (1) participants with no exposure to pesticides, (2) those exposed to non-mitochondria disrupting OPs and (3) those exposed to 1–7 and 8–14 mitochondria disrupting OPs.
Logistic regression analyses were conducted to assess associations between PD and the GIS-based estimates of individual and grouped ambient OP exposures at workplace and residential addresses. Trend tests were conducted using an ordinal count variable. We adjusted for age at diagnosis (cases) or age at interview (controls) as a dichotomous variable (≤60 or >60 years), gender, race (white vs non-white), education (<12, 12, >12 years), having a first-degree family member with PD (yes, no), smoking (current, former, never) and included an indicator variable for participants exposed to other pesticides, including organochlorines, dithiocarbamates and paraquat but not OPs.
We used SAS V.9.2 (SAS Institute Inc., Cary, North Carolina, USA) to perform unconditional logistic regression analyses.
Results
Study participants were predominantly Caucasian, over the age of 60 and a minority reported a family history of PD (table 1). Cases were slightly older than controls, more often male and had completed fewer years of formal education. They were also more likely to have never smoked cigarettes or to have stopped smoking.
Demographic characteristics of the study population
Effect estimates for each individual OP separately are shown in table 2. All OPs investigated were associated with 100% or greater increase in risk of developing PD. Generally, exposures at workplace addresses or at both residential and workplace addresses conferred a greater risk than exposures at residences only. Acephate, ethephon, phorate, naled, malathion, merphos, chlorpyrifos, disulfoton, dimethoate and monocrotophos were the only OPs that were strongly associated with increased risks of developing PD at residential addresses only, workplaces addresses only, as well as at residential and workplace addresses together.
Effect estimates (ORs and 95% CIs) for ambient exposures to organophosphates at residences and workplaces during 1974–1999
All groups of OPs classified by presumed mechanism/toxicity appeared to be associated with similarly increased odds of developing PD—that is, none of the groupings seemed more important than the others (table 3). Again, exposures at workplace addresses only or at both residential and workplace addresses were associated with a greater increase in risk of PD than exposures at residential addresses only. This pattern also held for OPs with presumed mitochondrial disrupting function—that is, we could not distinguish risk between those exposed to OPs with and without such presumed function (table 4). Moreover, our data suggested a trend where exposure to an increasing number of potential mitochondria disrupting OPs at workplaces (p-trend=0.0014), or at both residences and workplaces (p-trend < 0.0001), was associated with an increasing risk of PD.
Effect estimates for ambient exposures to subsets of organophosphates based on mechanism of toxicity at residences and workplaces during 1974–1999
Effect estimates for ambient exposures to presumed mitochondrial disrupting organophosphates at residences and workplaces during 1974–1999
Discussion
This population-based case–control study conducted in the Central Valley of California found exposure to all OPs included in this study to be associated with an increased odds of developing PD. Thus those classified as teratogens, endocrine disruptors and carcinogens exhibited similar effects as those assigned mitochondrial disrupting function, suggesting that the true mechanisms of neurotoxic action might not be known or that many OPs contributing to these groups have multiple mechanisms of toxicity. Also, we may not be able to distinguish between mechanisms due to OPs exhibiting multiple mechanisms—for example, 7 of the 12 OP endocrine disruptors were also presumed mitochondrial disruptors. Generally, ambient exposure at residential addresses resulted in weaker associations with PD than at workplace addresses or exposure in both locations, suggesting that workplace-related ambient exposures and exposure at multiple locations might result in higher personal exposure levels than exposures at residential addresses only. This may be due to the fact that most people do not work at their residence and thus are not present when the OPs are applied agriculturally during daytime work hours. The hydrolytic properties of OPs can cause the active ingredients to break down quickly outdoors and being present during or shortly after applications of OPs may contribute to the risk of developing PD.24 ,25
A strength of this study is the unique GIS-based pesticide exposure assessment method that was used to assess ambient pesticide exposure at residential and workplace addresses from drift and/or contact with dust and soils contaminated with pesticides in this heavily agricultural region.26 The specific type, amount and location of an applied pesticide active ingredient can be identified using historical PUR records, which is a vast improvement over the majority of previous studies that relied on less accurate exposure assessment, such as self-reported exposures. A recent meta-analysis of 46 studies showed a moderate association between general pesticide exposure and PD (summary RR: 1.62; 95% CI 1.40 to 1.88).3 However, it also found substantial heterogeneity of results between studies included in the analysis. The authors attributed this heterogeneity to unreliable self-reported exposure, which may have introduced non-differential misclassification that typically biases results towards the null, thereby preventing the detection of existing associations between pesticides and PD risk. Another issue with most of the studies included in the meta-analysis is that the majority of the studies defined pesticide exposure as ever/never exposure to any pesticide or to broad subgroups such as herbicides, insecticides and fungicides. This approach assumes that all pesticides or pesticides belonging to broad pesticide subgroups similarly affect the risk of developing PD, which is unlikely and may introduce bias due to exposure misclassification. Thus, estimates of the risk of developing PD after exposure to specific pesticides are probably needed to understand, at least in part, this heterogeneity.
A limitation of our exposure assessment method is that the precise resolution of the GIS-based pesticide exposure estimation relies on how accurately participants reported their workplace and residential addresses. Workplace addresses tend to be recalled with less accuracy than residential addresses leading to less precise geocoding. Indeed, workplace addresses were geocoded less precisely than residential addresses in this study. However, the geocoding quality of addresses for cases and controls was similar. We found that when assessing geocoded data in our study, 24.6% of cases and 21.2% of controls spent 50% or more of the years between 1974 and 1999 at precisely geocoded workplace addresses while 37.5% of cases and 48.4% of controls spent 50% or more of these years at precisely geocoded residential addresses. This suggests that geocoding precision is not likely to account for difference in the estimated effects.
A limitation of all pesticide research is the issue of correlated exposures. It is difficult to distinguish the effect that individual pesticides have on PD risk since people who are exposed to pesticides are often concurrently exposed to multiple pesticides. Thus, it is not possible to separate individual pesticide effects by simply adjusting for other pesticides. The ambient pesticide exposure estimates we derived contain an inherent amount of uncertainty as we were not able to measure historical pesticide serum levels in participants. Therefore, we decided to take a conservative approach and only consider participants exposed to a particular pesticide if they were assigned ambient exposure to that pesticide equal to or above the controls’ median exposure value. This increases the specificity of our exposure assessment, and participants who are less likely to be exposed will also be less likely to be misclassified as exposed.
Despite these limitations, our GIS model provides a valid and high-quality indicator of passive pesticide exposure from applications and drift in close proximity to workplaces and residences. It is unlikely that the GIS-based results are affected by selection bias because our participants did not self-report pesticide exposures and were most likely not aware of their historical exposures to specific pesticides applied within 500 m of their residences and/or workplaces. There is also no reason to suspect that cases and controls would choose to differentially participate in this study based on whether or not they lived and/or worked near agricultural plots during the 26-year exposure period we investigated.
The primary mechanism ascribed to OP toxicity is the inhibition of acetylcholinesterase (AChE) leading to an accumulation of acetylcholine at cholinergic synapses and overstimulation of muscarinic and nicotinic receptors.6 Although inhibition of AChE results in cholinergic system overstimulation and consequently cell death,27 ,28 it is unlikely to be the primary mechanism for PD development.5 OPs also may exhibit a number of non-cholinergic mechanisms of toxicity such as the inhibition of mitochondrial processes.5 Concerning the presumed mitochondrial disrupting properties of OPs, this mechanism has not been studied extensively yet and the majority of the existing studies used large experimental doses of OPs that may have caused non-specific general toxicity and are not reflective of real-world ambient exposures and potential mechanisms. Thus, the current literature is insufficient to establish whether individual OPs are mitochondrial disruptors at ambient exposure levels and our results do not support such a specific mechanism for the pesticides we evaluated.
Oxidative stress occurs when the body's antioxidant defences are not able to neutralise an excess of reactive oxygen species (ROS).29 The brain uses a large amount of oxygen that may lead to an increased production of ROS and is also particularly susceptible to oxidative damage as it is composed largely of polyunsaturated fatty acids.28 When AChE is inhibited, an excessive amount of ROS may be generated when cells are not able to maintain energy levels due to high-energy consumption and inhibition of oxidative phosphorylation.30 Oxidative stress may also originate from the mitochondria when their function is disrupted.31 ,32 Dopaminergic neurons, unlike other neurons, are autonomously active pacemakers and have a greater reliance on oxidative phosphorylation.33 However, grouping OP pesticides according to some presumed function/toxicity did not allow us to distinguish OPs from each other according to whether or not they increased PD risk differently in our human population. Besides governing aerobic respiration, mitochondria also are involved in the apoptotic neurodegenerative processes.5 Apoptosis may be caused by increased production of ROS and the translocation and inhibition of proteins involved in respiration such as cytochrome c.34 Mevinphos has been shown to disrupt oxidative phosphorylation by causing the dysfunction of complexes I through IV, which leads to cell death due to ATP depletion.35 Monocrotophos have also been shown to induce apoptosis in neurons and inhibit metabolism.34 Chlorpyrifos and chlorpyrifos-oxon exposure resulted in increased mitochondrial length, decreased number of mitochondria and decreased mitochondria movement in axon at concentrations that did not inhibit AChE.36 Taken together, the evidence that OPs cause cell death specifically mediated through the disruption of mitochondrial functioning is not adequate and our results suggest that while a large number of OPs increased PD risk the mechanisms of action remain to be explored.
In conclusion, this study adds strong evidence that OPs are implicated in the aetiology of idiopathic PD. Additionally, ambient exposure to OPs at workplaces and combined ambient exposure at residences and workplaces seem to be especially important. This is the first study to show in a human population that exposure to increasing numbers of OPs is associated with elevated risks of PD. Future studies should examine possible neurotoxic mechanisms of OPs at low doses that are reflective of real-world ambient exposure.
What this paper adds
-
Commonly used pesticide exposure estimation techniques are not sophisticated enough to assess concurrent exposure to specific pesticides, which is vital to understanding the aetiology of Parkinson's disease.
-
This study relies on a geographic information system-based approach to assign historical pesticide use to participant addresses and derive exposures to specific organophosphate pesticides.
-
All of the specific organophosphate pesticides included in the study increased the risk of Parkinson's disease, but the mechanism(s) behind this risk increase is unclear.
-
In order for epidemiological studies to assess the neurotoxicity of specific pesticides in human populations, more animal and cell model research must be conducted to examine neurotoxic mechanisms of pesticides at low doses that are reflective of real-world ambient exposure.
References
Supplementary materials
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.
Files in this Data Supplement:
- Data supplement 1 - Online appendix
Footnotes
-
Contributors AW performed all data analyses and drafted the manuscript as its first author. MC and TTL developed the geographic information system-based exposure assessment tool, supervised and conducted the geocoding and pesticide exposure assessment, and reviewed the manuscript. JMB served as the lead neurologist and performed clinical exams to confirm PD diagnosis. BR is the primary investigator of the study and was responsible for its design and conduct. She was also involved in the conceptualisation of data analysis, the presentation of the study results and the writing of the manuscript. All authors have discussed the topic, drafted and approved the final manuscript.
-
Funding This work was supported by the National Institute of Environmental Health Science (ES10544, U54ES12078, 5P30 ES07048), National Institute of Neurological Disorders and Stroke (NS 038367) and Department of Defense Prostate Cancer Research Program (051037); in addition, initial pilot funding was provided by the American Parkinson's Disease Association.
-
Competing interests None.
-
Ethics approval UCLA IRB.
-
Provenance and peer review Not commissioned; externally peer reviewed.