Article Text

Original article
Pesticide use and risk of end-stage renal disease among licensed pesticide applicators in the Agricultural Health Study
Free
  1. Jill F Lebov1,
  2. Lawrence S Engel1,
  3. David Richardson1,
  4. Susan L Hogan2,
  5. Jane A Hoppin3,
  6. Dale P Sandler4
  1. 1Department of Epidemiology, University of North Carolina, Chapel Hill, North Carolina, USA
  2. 2Department of Medicine, Division of Nephrology and Hypertension, University of North Carolina, Chapel Hill, North Carolina, USA
  3. 3Department of Biological Sciences, Center for Human Health and the Environment, North Carolina State University, Raleigh, North Carolina, USA
  4. 4Epidemiology Branch/Chronic Disease Epidemiology Group, National Institute of Environmental Health Sciences, Research Triangle Park, North Carolina, USA
  1. Correspondence to Dr Dale P Sandler, Epidemiology Branch/Chronic Disease Epidemiology Group, National Institute of Environmental Health Sciences, P.O. Box 12233, MD A3-05, 111 T.W. Alexander Drive, Research Triangle Park, NC 27709, USA; sandler{at}niehs.nih.gov

Abstract

Objectives Experimental studies suggest a relationship between pesticide exposure and renal impairment, but epidemiological evidence is limited. We evaluated the association between exposure to 39 specific pesticides and end-stage renal disease (ESRD) incidence in the Agricultural Health Study, a prospective cohort study of licensed pesticide applicators in Iowa and North Carolina.

Methods Via linkage to the US Renal Data System, we identified 320 ESRD cases diagnosed between enrolment (1993–1997) and December 2011 among 55 580 male licensed pesticide applicators. Participants provided information on use of pesticides via self-administered questionnaires. Lifetime pesticide use was defined as the product of duration and frequency of use and then modified by an intensity factor to account for differences in pesticide application practices. Cox proportional hazards models, adjusted for age and state, were used to estimate associations between ESRD and: (1) ordinal categories of intensity-weighted lifetime use of 39 pesticides, (2) poisoning and high-level pesticide exposures and (3) pesticide exposure resulting in a medical visit or hospitalisation.

Results Positive exposure-response trends were observed for the herbicides alachlor, atrazine, metolachlor, paraquat, and pendimethalin, and the insecticide permethrin. More than one medical visit due to pesticide use (HR=2.13; 95% CI 1.17 to 3.89) and hospitalisation due to pesticide use (HR=3.05; 95% CI 1.67 to 5.58) were significantly associated with ESRD.

Conclusions Our findings support an association between ESRD and chronic exposure to specific pesticides, and suggest pesticide exposures resulting in medical visits may increase the risk of ESRD.

Clinical trial registration Clinicaltrials.gov NCT00352924.

Statistics from Altmetric.com

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.

What this paper adds

  • Much is known about clinical risk factors for end-stage renal disease (chronic kidney disease requiring dialysis or kidney transplant for survival), but research on environmental risk factors for kidney disease is limited.

  • In this study of male pesticide applicators, risk of end-stage renal disease increased with increasing cumulative exposure to several pesticides, including the herbicides alachlor, metolachlor, paraquat and pendimethalin, and the insecticide permethrin.

  • Risk of end-stage renal disease was significantly greater for pesticide applicators who reported multiple doctor visits due to pesticide use and hospitalisation due to pesticide use, compared with those who reported no medical visits due to pesticide use.

  • Exposure to certain pesticides may increase the risk of end-stage renal disease; however, additional studies are needed to support these findings.

Introduction

In 2011, over 600 000 US residents were receiving treatment for end-stage renal disease (ESRD), a life-threatening condition requiring dialysis or kidney transplant for survival. Although much research has been conducted on clinical precursors to ESRD, such as diabetes and hypertension, research is limited on the impact of environmental and occupational factors. Literature on the nephrotoxic effects of pesticides in humans is extremely limited, and mainly comes from case reports of both fatal and non-fatal pesticide poisoning, which have described acute nephrotoxicity with a variety of pesticide classes.1–5 Most of the evidence regarding nephrotoxicity of pesticides is restricted to experimental animal studies. Renal damage and dysfunction have been observed in experimental animal studies with exposure to specific pesticides in a dose-dependent6–9 and/or exposure duration-dependent7 ,9–11 manner. A variety of pesticide classes have been shown to cause renal damage and dysfunction in animals, including organophosphate,6 ,9 organochlorine,10 ,12 carbamate13 and pyrethroid14 insecticides and triazine15 and chlorophenoxy8 herbicides.

The impact of long-term pesticide exposure on human kidney function remains largely unknown. Generally, studies conducted in El Salvador, Nicaragua and Sri Lanka indicate an elevated prevalence of chronic kidney disease (CKD) among agricultural workers16–19 compared to those who have never worked in agriculture, particularly among male agricultural workers;20 pesticide exposure is postulated to be a contributor to kidney disease in these regions, but existing evidence has not confirmed this hypothesis.16–18 These studies lack specificity with regard to pesticide type and have not been able to adequately assess the long-term effects of chronic low-level or acute high-level pesticide exposure on ESRD risk. To our knowledge, the only study to assess pesticide exposure and ESRD found self-reported work in a place with frequent or daily exposure to insect or plant spray to be associated with increased ESRD risk.21

The Agricultural Health Study (AHS) is the largest prospective study of pesticide applicators in the USA. Linking the AHS to the US Renal Data System (USRDS) provides a unique opportunity to evaluate the relationship between pesticide use and ESRD risk. Using this linkage, we evaluated associations between chronic and acute pesticide exposure and ESRD risk.

Methods

Population and case definition

Details of the AHS design have been described previously.22 Briefly, the AHS recruited private pesticide applicators (mainly farmers) (N=52 394) in North Carolina and Iowa, and commercial pesticide applicators (N=4916) in Iowa who applied for or renewed a restricted-use pesticide license between 1993 and 1997. At the licensing facility, each pesticide applicator was asked to complete a brief enrolment questionnaire. Participating applicators were also given a packet of additional questionnaires to complete at home and mail back (take-home questionnaire). Approximately 82% of eligible private applicators and 47% of eligible commercial applicators enrolled in the study. Farm workers were not included in the AHS. At enrolment, applicators provided information on lifetime pesticide use and pesticide use practices, demographic characteristics, lifestyle factors, farm information, and medical history in a self-administered questionnaire.22 Of enrolled applicators, 44% also completed the take-home questionnaire with additional questions about medical history and pesticide use.23 Questionnaires are available on the AHS web site: http://aghealth.nih.gov/collaboration/questionnaires.html.

We identified ESRD cases diagnosed between study enrolment and end of follow-up (31 December 2011) through linkage with the USRDS. The USRDS collects data on all ESRD cases in the USA through Medical Evidence Form CMS-2728, which is required for all new patients with ESRD, regardless of Medicare eligibility. The USRDS derives the first ESRD service date (FSD) by taking the earliest of: (1) the date of the start of dialysis for chronic renal failure, as reported on the Medical Evidence form; (2) the date of a kidney transplant or (3) the date of the first Medicare dialysis claim.24 The FSD was used to estimate age at ESRD diagnosis. Date of death was obtained from state mortality files and the National Death Index. Because the distribution of ESRD risk factors differs by gender, and because few applicators were female, we excluded female applicators (N=1562; 2.7%) from this analysis. We also excluded applicators under age 18 (N=127; <1.0%) and ESRD cases diagnosed prior to enrolment (N=42; 11.5% of cases). This left us with 55 580 participants for analyses of enrolment questionnaire variables, and 24 565 participants for analyses of take-home questionnaire variables.

Exposure assessment

Information provided on the enrolment and take-home questionnaires was used to estimate lifetime pesticide exposure. Though information on pesticide use is available in later phases of the AHS study, we limited these analyses to enrolment information because recent exposure was anticipated to have a minor impact on ESRD risk, given the typical decades-long progression of this disease from chronic stage 1 to ESRD, and recent exposure information was not available for all participants. Participants provided information on years of use (duration) and average days per year of use (frequency) for 22 pesticides on the enrolment questionnaire. Duration and frequency of use data were obtained on the take-home questionnaire for 28 additional pesticides. For each pesticide, an intensity-weighted lifetime exposure metric was generated by multiplying lifetime-days of use (product of duration and frequency of use) by an intensity score that accounts for differences in exposure resulting from variation in pesticide application methods, repair of pesticide application equipment, and use of personal protective equipment.25 The intensity score algorithm was developed using AHS-specific pesticide exposure monitoring data, in conjunction with expert judgment from published studies on pesticide exposure, including information from the Pesticide Handlers Exposure Database.25 We used intensity-weighted lifetime-days as our primary exposure metric. Owing to the relatively small number of cases, we used the distribution of use among cases to create cut points for intensity-weighted lifetime use of specific pesticides. For pesticides used by ≥15% of cases, we categorised non-zero intensity-weighted lifetime-days into tertiles with non-users as the referent group. For less frequently used pesticides, non-zero intensity-weighted lifetime-days of use were split at the median (<median vs ≥ median), with non-users as the referent group. Analyses were restricted to pesticides for which there were at least five cases in each exposure stratum. This left 39 pesticides available for analyses.

To assess overall pesticide use, we evaluated risk related to duration, frequency and lifetime-days of use of any pesticide. Duration and frequency of use were categorised into three levels (lowest category of use (referent), >lowest category of use to the median value, and >median). Cumulative lifetime-days of use was categorised into quartiles.

Participant report of medical visits due to pesticide use (enrolment questionnaire), unusually high personal exposure to any pesticide (take-home questionnaire), and doctor-diagnosed pesticide poisoning (take-home questionnaire) were also evaluated in relation to ESRD risk.

Statistical analysis

We used Cox proportional hazards models to calculate HRs for risk of ESRD, using age as the timescale and adjusting for state as a covariate in all models. Person-time was accrued from the date of study enrolment until the earliest of ESRD diagnosis, death, or the end of study follow-up (31 December 2011). The proportional hazards assumption was evaluated for each model by entering a product term (exposure of interest×time on study) into each model. A product term coefficient which differed significantly from zero (χ2 p value <0.10) indicated a potential violation of the proportional hazards assumption.

Private and commercial applicators were analysed together because there were too few ESRD cases among the latter to analyse them separately. Race and education level were identified as additional potential confounders through directed acyclic graph analyses (DAGs) informed by a review of prior literature. Because adjustment for these factors did not substantially change HR estimates, and power was reduced due to incomplete ascertainment of education and race data, we did not adjust for these factors in the final analyses. Although diabetes and body mass index (BMI) were associated with ESRD risk in this and other studies,26 ,27 it was unknown whether these conditions affected pesticide use, and it is possible that use of specific pesticides may increase the risk of diabetes and BMI,28–31 which may be on the causal pathway to ESRD. HR estimates obtained from models adjusted for diabetes and BMI were not meaningfully different from crude estimates; therefore, these two conditions were not included in final adjusted models. Hypertension was not adjusted for because it is largely asymptomatic and therefore is unlikely to have affected pesticide use practices. Additionally, in our study, there was no association between cumulative pesticide use and hypertension (data not shown).

To assess potential confounding by other pesticides, we examined pairwise correlations between pesticides that were strongly (HR in any strata ≥1.5 or ≤0.65) or significantly associated with ESRD in single pesticide-adjusted exposure-response models. For pesticides with a Spearman correlation coefficient ≥0.3, we constructed models with both pesticides, using the intensity-weighted variables included in the main analyses. For pesticides that were correlated with more than one pesticide, we first evaluated each pesticide pair and then added correlated pesticides one at a time into subsequent models. We assessed model fit using Akaike information criteria (AIC) and selected that with the lowest AIC as the final model.

To assess linear exposure-response trends in intensity-weighted lifetime use, we used within-category medians as the score for each level of use for each chemical. Exposure-response trends were also evaluated for duration, frequency and cumulative lifetime-days of use of any pesticide, and number of doctor visits related to pesticide use.

ESRD is the final stage of CKD, which is often debilitating in later stages of the disease. Cases may have already experienced the effects of CKD prior to study enrolment, which could have influenced their pesticide use. If those with earlier stages of renal disease have reduced exposure due to modified application practices, effect estimates for specific pesticide use would be biased towards the null. This bias is commonly referred to as the healthy worker survivor effect.32 To evaluate the potential for this effect to influence our findings, we repeated analyses, excluding person-time for all participants for the first 5 years after enrolment under the assumption that ESRD cases diagnosed within 5 years following study enrolment likely had poor renal health at enrolment.

To evaluate whether patterns of association were consistent across states, we entered a product term for state into pesticide use models for those pesticides for which there were at least five cases in each stratum of use in both states.

We used the AHS data set releases P1REL201209, P3REL201209.00, and AHSREL201304.00. All statistical analyses were done using SAS V.9.3 (Cary, North Carolina, USA).

Results

Of the 55 580 participants eligible for analysis, 320 (308 private and 12 commercial) were diagnosed with ESRD over an average 15.7-year follow-up period (incidence rate: 36.6 ESRD cases per 100 000 person-years). Among the subset of 24 565 participants who returned the take-home questionnaire, there were 136 cases (incidence rate: 35.1 ESRD cases per 100 000 person-years). ESRD incidence was significantly higher in North Carolina compared with Iowa, regardless of age, which follows the pattern of ESRD incidence in the general population.24 In age-adjusted and state-adjusted models, education level of high school or lower and obesity at enrolment, were associated with increased risk of ESRD (table 1). Self-reported doctor diagnosis of diabetes, high blood pressure and kidney disease were significantly associated with increased risk of ESRD. There was a suggestive but non-significant association between pack-years of cigarettes smoked and ESRD, but applicator type, number of years living on a farm, and alcohol consumption at enrolment were not associated with ESRD risk.

Table 1

Association between ESRD and demographic and medical conditions among private and commercial applicators, adjusted for age and state, Agricultural Health Study (1993–1997)

More than one doctor visit due to pesticide use and hospitalisation due to pesticide use were both significantly associated with ESRD with a significant trend observed for increasing number of pesticide-related doctor visits (p for trend=0.038) (table 2). ESRD risk was not associated with either self-reported unusually high personal pesticide exposure or pesticide poisoning, though only five ESRD cases reported a pesticide poisoning diagnosis. No exposure-response relationships were observed for duration and frequency of general pesticide use (data not shown), or for cumulative lifetime-days of general pesticide use (table 2).

Table 2

Association between ESRD risk and acute and cumulative pesticide exposure, adjusted for age and state, among male pesticide applicators, Agricultural Health Study (1993–1997)

In intensity-weighted cumulative use analyses, positive associations were observed primarily among herbicides (table 3). ESRD risk was associated with the highest tertile of intensity-weighted use of five herbicides: atrazine, metolachlor, alachlor, paraquat and pendimethalin, compared to no use. We observed a significant (p for trend <0.05) exposure-response trend with increasing use levels for all these herbicides. Although exposure-response trends were not seen for the herbicides petroleum oil or imazethapyr, ever use of these chemicals was significantly associated with risk (HR=1.63; 95% CI 1.11 to 2.41 and HR=1.46; 95% CI 1.08 to 1.99, respectively; data not shown). The proportional hazards assumption held for all exposures of interest.

Table 3

Intensity-weighted lifetime-days of use of specific pesticides and ESRD risk, adjusted for age and state, among male pesticide applicators in the Agricultural Health Study (1993–1997)

Among non-herbicide pesticides, associations for the insecticides coumaphos, parathion and phorate (organophosphates), aldicarb (carbamate), chlordane (organochlorine), and for the fungicide metalaxyl were elevated (ie, >1.6), but did not reach statistical significance (table 3). A significant positive exposure-response trend was present for permethrin (for crops) (p for trend=0.0311).

In analyses of correlated pesticides, we found 13 pesticides to be correlated (ie, had a Spearman correlation coefficient ≥0.30) with at least one other pesticide. Adjustment for correlated pesticides resulted in reduced overall sample size due to missing data for each chemical. Adjusted estimates were generally similar in magnitude and direction, but were attenuated and less precise. Patterns of exposure-response also did not change. After adjustment for correlated pesticides, the association between ESRD risk and the top tertile of intensity-weighted use remained significant only for pendimethalin. The association with the top tertile of alachlor remained elevated but no longer significantly associated with ESRD risk after adjustment for correlated pesticides. We did not observe a correlation coefficient ≥0.3 for the following pesticides: 2,4,5-Trichlorophenoxyacetic acid (2,4,5-T) 2,4,5T, chlorimuron ethyl, paraquat, petroleum oil, coumaphos, fonofos, parathion and permethrin (crops) (data not shown).

In a sensitivity analysis, evaluating the potential for a ‘healthy worker survivor effect’, we excluded 53 cases that were diagnosed with ESRD within 5 years after enrolment and 277 900 person-years. The greatest per cent reductions of case numbers were observed in the ‘None’ use category for all pesticides. In general, associations for intensity-weighted lifetime-use were in the same direction and of very similar magnitude compared with estimates in the main analyses. Of note, age-adjusted and state-adjusted estimates for the highest quantile of intensity-weighted coumaphos use became significant (HR=1.81, 95% CI 1.03 to 3.17, and we observed a significant positive exposure-response trend for chlordane (p for trend=0.0177) (see online supplementary table S1).

Results did not change substantially when we restricted analyses to private applicators (data not shown). Thirteen pesticides had five exposed cases in each exposure stratum in both states, and were therefore included in analyses of interaction with state. P-values for interaction with state were consistently >0.10 (data not shown), suggesting no differences by state.

Discussion

To the best of our knowledge, this is the first study to evaluate the association between ESRD risk and cumulative lifetime use of specific pesticides. Among pesticide applicators in the AHS, we found significant positive associations between intensity-weighted use of several specific pesticides and ESRD; excluding cases that arose within 5 years after enrolment strengthened some of these associations, though estimates were less precise due to the reduction in sample size. This is also the first epidemiological study of ESRD risk associated with non-fatal pesticide poisoning, acute high-level exposure, and pesticide exposure requiring medical attention. Participants who reported doctor visits and hospitalisation due to pesticide use had a significantly higher risk of ESRD diagnosis compared to those who did not, but we did not observe increased risk with applicator report of doctor-diagnosed pesticide poisoning or unusually high personal pesticide exposure.

Prior published epidemiological research on pesticide exposure and kidney disease is minimal. Results from several cross-sectional studies evaluating the relationship between agricultural work and CKD suggest a potential association between agricultural work, particularly field work, and CKD prevalence.17 ,20 Studies that have evaluated overall pesticide exposure have found positive associations with CKD.17 ,18 ,33 ,34 The only study to assess the relationship between ESRD risk and agricultural exposures observed an increased risk of ESRD among a large population of insured patients in the San Francisco Bay area who reported that they worked in a place with ‘frequent or daily exposure to insect or plant spray’ (unadjusted HR=1.78; 95% CI 1.36 to 2.34).21 By contrast, results from our analyses of general overall pesticide use did not show an association with ESRD; however, few participants in this licensed applicator cohort reported no pesticide use, and evaluation of overall pesticide use may obscure associations because only some pesticides appear to be associated with ESRD.

Epidemiologic studies of renal effects of specific pesticides are rare. Hernandez et al35 found no difference in serum creatinine levels among greenhouse workers with higher versus lower levels of apparent cholinesterase inhibition (used as a marker for organophosphate pesticide exposure). Serum levels of several organochlorine insecticides among patients with CKD were inversely associated with kidney function, potentially indicating a renal filtration deficiency resulting in an accumulation of organochlorine pesticides in the body.36 In our study, ESRD risk was elevated for three cholinesterase-inhibiting insecticides (the carbamate aldicarb and the organophosphates coumaphos and parathion), with a moderate positive trend observed for coumaphos. No associations were seen with organochlorine use, except for chlordane, which was significantly positively associated with ESRD risk in analyses excluding cases diagnosed within 5 years after enrolment. Glyphosate was recently partially banned in Sri Lanka due to its hypothesised association with kidney disease, though the ban has since been lifted. Studies that informed this partial ban suggested that glyphosate exposure leads to renal failure only when combined with high-level exposure to heavy metals.37 We found no evidence of an association between ESRD risk and glyphosate exposure.

Experimental evidence supports our findings of positive associations with exposure to the herbicides atrazine, alachlor, paraquat and pendimethalin, the insecticide permethrin, and the fungicide metalaxyl, with evidence of dose-response as well as renal damage and dysfunction at low dose levels. Glomerular lesions and renal tubular necrosis due to oxidative stress-induced cell damage have been observed in animal models with exposure to metalaxyl and paraquat,38 ,39 and kidney damage and dysfunction have been observed in rats exposed to atrazine15 and fish exposed to alachlor.40 There have been no reports of renal effects of pendimethalin among mammals; however, at least one formulation of pendimethalin contains monochlorobenzene as an inert ingredient, which has been shown to cause kidney damage in rats.41 Dose-related renal damage was observed among animals exposed to permethrin.42 ,43 Sub-acute tubulointerstitial and glomerular damage, such as that observed in animal studies with prolonged low dosing of pesticides,7 ,11 can initiate a feed-forward loop of kidney injury and progressive loss of renal function.44 In humans, pesticide poisoning can lead to acute kidney injury,45 which has been associated with increased risk of subsequent CKD and ESRD.46 ,47

When we adjusted for correlated pesticides, we found that many estimates were attenuated, and alachlor was positively associated but no longer statistically significant. We lacked statistical power to formally test interaction among specific pesticides. In addition, we saw significant associations for chemicals from a number of chemical groups. Given that there is no required renal toxicity testing for pesticides, we have no information regarding which chemicals we would expect to be associated with CKD. A large portion of commercial pesticide products are ‘other ingredients’ including solvents; therefore, it is possible that one or more of these ‘other ingredients’ are driving the risk observed here for unrelated chemicals. Unfortunately, we do not have access to what these ‘other ingredients’ are because they are regarded as confidential business information.

A significant positive exposure-response trend was observed for pesticide-related doctor visits, and participants who reported being hospitalised due to pesticide use had three times the risk of ESRD compared to those who did not. Although information about the route and type of pesticide involved in these exposures was not available, these findings support the hypothesis that frequent and/or severe pesticide exposures may increase the risk of ESRD. We did not see an association between pesticide poisoning or self-reported unusually high pesticide exposure and ESRD; however, power to detect an association was limited because information for those exposures was available only for participants who returned the take-home questionnaire, and pesticide poisoning was rarely diagnosed in the cohort. Pesticide poisoning is frequently under-diagnosed;48 thus, doctor visits or hospitalisation related to pesticide use may represent a more sensitive indicator for acute high-level exposure than pesticide poisoning diagnosis.

The study improves on prior research in several ways. First, we were able to evaluate associations of ESRD risk with a wide range of specific chemicals that vary in toxicity and extent of use. The prospective design of the study mitigates concerns about differential misclassification of exposure. Whereas prior studies were limited by small sample size and exposure to few chemicals, the large size of the AHS cohort allowed us to assess exposure-response trends for many individual pesticides. Also, use of a validated exposure-intensity metric25 allowed for a better estimate of each participant's likely pesticide exposure as opposed to non-specific pesticide use. Until now, the relationship between short-term high-level pesticide exposures and kidney disease has been evaluated only with respect to the immediate effects of pesticide poisonings. Here, we were able to evaluate measures of non-poisoning short-term high-level pesticide exposures, thereby providing an important contribution to the scientific literature regarding occupational risk factors for kidney disease. Additionally, the fact that almost all ESRD cases in the USA are captured in the USRDS reduces concerns about loss to follow-up or outcome misclassification.

Because exposure data were collected prior to disease onset and ESRD diagnosis data were obtained from a third party linkage rather than participant report, any exposure misclassification due to self-report is likely to be non-differential with respect to the outcome, which would bias estimates towards the null. Additionally, evidence suggests that report of pesticide exposure by AHS participants is reasonably reliable49 and plausible.50 The accuracy and reliability of reporting acute pesticide exposures has not been investigated.

Analyses of lifetime use of pesticides that were assessed only on the take-home questionnaire could be subject to selection bias if applicators who returned the take-home questionnaire were significantly different from those who did not by exposure or outcome, or by factors associated with the exposure or the outcome. State of enrolment, age and pesticide use characteristics were similar for those who returned the take-home questionnaire compared to those who did not, and the prevalence of ESRD was essentially the same for take-home questionnaire respondents as it was for non-respondents (0.59% vs 0.57%). Still, differences in unmeasured factors remain a possibility, and we had limited power to evaluate associations between ESRD and lifetime use of the pesticides for which duration and frequency information was collected only on the take-home questionnaire.

Lifetime use estimates could also be biased if participants with prevalent pre-end-stage kidney disease modified their pesticide application practices in the period prior to enrolment. This ‘healthy worker survivor effect’ is a common problem in occupational health studies, including among pesticide applicators,51 frequently biasing estimates toward the null.32 Results of our sensitivity analysis suggest a minimal impact of this potential effect on HR estimates. The median time from enrolment to ESRD diagnosis was 9.7 years. It is possible that unmeasured post-enrolment pesticide use may differ by case status; if recent exposures are stronger contributors to ESRD risk than pre-enrolment exposures, then our results would still be biased towards the null. However, progression of renal disease from chronic stage 1 to ESRD can take several decades; if pesticide use does contribute to kidney disease incidence, it is probable that this pathway would have been initiated prior to enrolment.

Because the AHS is a cohort of US farmers in North Carolina and Iowa, the risks for this cohort may be different from the risks experienced by farm workers in the USA or elsewhere. However, the AHS, with its detailed exposure characterisation, limited loss to follow-up over 20 years, and the ability to link to a population-based kidney disease registry, provides important human data for evaluation of the role of pesticides and kidney disease.

Last, in these analyses, we estimated associations between a large number of exposures and ESRD. We did not employ statistical methods to adjust for multiple comparisons, such as the Bonferroni correction. Previous authors have cautioned against employing Bonferroni-type corrections in an epidemiological analysis of associations between multiple environmental or occupational exposures and disease,52–54 and such methods have largely fallen out of use in such settings because Bonferroni adjustments are concerned with testing of a general null hypothesis which is rarely of interest. Alternatives to Bonferroni methods exist for inference in settings where multiple exposure effects are estimated, including Bayesian methods;55 however, we have not employed such methods here, in preference for simply describing what statistical quantities have been estimated.

Conclusions

Our study provides evidence for an association between ESRD risk and chronic exposure to specific chemicals among pesticide applicators in Iowa and North Carolina. Results from this study also suggest that pesticide exposures resulting in medical visits increase the risk of incident ESRD, raising concerns that multiple high-level pesticide exposures may contribute to irreversible kidney damage and resultant disease. Efforts to better characterise the pathway between pesticide exposure and kidney disease should include assessments of earlier disease stages, rate of progression from CKD to ESRD, and other potential routes of pesticide exposure, such as spray drift and carry-home exposures. Caution should be taken in interpreting results of such studies when diagnosis dates or disease severity information is not available, because the healthy worker survivor effect may bias estimates towards the null. Additional epidemiological studies are needed to confirm the findings of our study, given the limited research on the role of pesticide exposure in the development of renal disease, and research on the direct renal toxicity of specific chemicals must be expanded to facilitate interpretation of epidemiological results.

Acknowledgments

The authors thank the participants of the AHS for their contribution to this research. They also thank Aaron Blair, Freya Kamel, Honglei Chen, Cynthia J Hines and Laura Beane Freeman for reviewing this manuscript, and Stuart Long at Westat for his assistance with data management.

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.

Footnotes

  • Correction notice This paper has been amended since it was published Online First. The data for imazethapyr in Table 3 was misaligned and the cutpoints for cumulative use (left most column) were not there and these have now been added.

  • Contributors JFL conducted the statistical analyses, interpreted the results and drafted the manuscript. As principal investigators of the Agricultural Health Study, DPS and JAH generated the underlying cohort study data, and conceived and designed the present study. DPS secured the funding for this project. LSE, DR, SLH, DPS and JAH provided general oversight for the project and contributed to the interpretation of results. Westat provided data management services. All authors provided critical insights and editing on drafts and read and approved the final manuscript.

  • Funding This research was supported in part by the Intramural Research Program of the NIH, National Institute of Environmental Health Sciences (ZO1 ES 049030) and by an NIH, National Institute of Diabetes and Digestive and Kidney Diseases Ruth L Kirschstein National Research Service Award (NRSA) Institutional Research Training Grant (T32 DK 007750; PI RJ Falk).

  • Competing interests None declared.

  • Ethics approval Institutional Review Boards at UNC Chapel Hill, the NIH and the participating field sites.

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

Linked Articles