Pyrethroid pesticide exposures may be associated with the onset of depression in later life via disruption of dopaminergic, serotonergic, and neurological functioning. We sought to investigate the association between living near agricultural pyrethroid pesticide applications and depression measures in central California, using two waves (PEG 1&2, total N = 1,654) of a case control study of Parkinson's disease (PD). At enrollment, participants self-reported history of use of depression medications and dates of MD-diagnosed depression and anxiety. Participants also completed a Geriatric Depression Scale-Short Form upon enrollment. We used the California Pesticide Use Registry to assign estimated ambient pyrethroid pesticide exposures at participant's home addresses over the 5 years before the index date (date of outcome, or an age-matched year for participants without the outcome). We used logistic and linear regression to evaluate associations between living near any pyrethroid applications over the 5-year index period and measures of depression and anxiety. We also evaluated modification by study wave and PD status. We observed associations of pyrethroids with depression, depression medications, and anxiety (adjusted odds ratio [aOR] depression = 1.54, 95% confidence interval [CI] 1.14, 2.07; aOR depression medications = 1.68, 95% CI 1.25, 2.25; aOR anxiety = 1.60, 95% CI 1.17, 2.18). However, we observed no associations with mild/moderate depressive symptoms according to the GDS score at enrollment (aOR = 1.04, 95% CI 0.77, 1.42). We did not observe a consistent modification of the pyrethroid-depression associations by study wave and PD status. Ambient pyrethroid pesticide exposures may be associated with measures of depression in later life.

Background

Pyrethroid pesticide use is increasing worldwide, although the full extent of associated health effects is unknown. An epigenome-wide association study (EWAS) with exploratory pathway analysis may help identify potential pyrethroid-related health effects.

Methods

We performed an exploratory EWAS of chronic ambient pyrethroid exposure using control participants' blood in the Parkinson's Environment and Genes Study in the Central Valley of California (N = 237). We estimated associations of living and working near agricultural pyrethroid pesticide applications in the past 5 years (binary) with site-specific differential methylation, and used a false discovery rate (FDR) cut off of 0.05 for significance. We controlled for age, sex, education, cell count, and an ancestral marker for Hispanic ethnicity. We normalized methylation values for Type I/II probe bias using Beta-Mixture Quantile (BMIQ) normalization, filtered out cross-reactive probes, and evaluated for remaining bias with Surrogate Variable Analysis (SVA). We also evaluated the effects of controlling for cell count and normalizing for Type I/II probe bias by comparing changes in effect estimates and p-values for the top hits across BMIQ and GenomeStudio normalization methods, and controlling for cell count. To facilitate broader interpretation, we annotated genes to the CpG sites and performed gene set overrepresentation analysis, using genes annotated to CpG sites that were associated with pyrethroids at a raw p < 0.05, and controlling for background representation of CpG sites on the chip. We did this for both a biological process context (Gene Ontology terms) using missMethyl, and a disease set context using WebGestalt. For these gene set overrepresentation analyses we also used an FDR cut off of 0.05 for significance of gene sets.

Results

After controlling for cell count and applying BMIQ normalization, 4 CpG sites were differentially methylated in relation to pyrethroid exposures. When using GenomeStudio's Illumina normalization, 415 CpG sites were differentially methylated, including all four identified with the BMIQ method. In the gene set overrepresentation analyses, we identified 6 GO terms using BMIQ normalization, and 76 using Illumina normalization, including the 6 identified by BMIQ. For disease sets, we identified signals for Alzheimer's disease, leukemia and several other cancers, diabetes, birth defects, and other diseases, for both normalization methods. We identified minimal changes in effect estimates after controlling for cell count, and controlling for cell count generally weakened p-values. BMIQ normalization, however, resulted in different beta coefficients and weakened p-values.

Conclusions

Chronic ambient pyrethroid exposure is associated with differential methylation at CpG sites that annotate to a wide variety of disease states and biological mechanisms that align with prior research. However, this EWAS also implicates several novel diseases for future investigation, and highlights the relative importance of different background normalization methods in identifying associations.

Pesticides are widely used in the agricultural Central Valley region of California. Historically, this has included organophosphates (OPs), organochlorines (OCs), and pyrethroids (PYRs). This study aimed to identify perturbations of the serum metabolome in response to each class of pesticide and mutual associations between groups of metabolites and multiple pesticides. We conducted high-resolution metabolomic profiling of serum samples from 176 older adults living in the California Central Valley using liquid chromatography with high-resolution mass spectrometry. We estimated chronic pesticide exposure (from 1974 to year of blood draw) to OPs, OCs, and PYRs from ambient sources at homes and workplaces with a geographic information system (GIS)-based model. Based on partial least-squares regression and pathway enrichment analysis, we identified metabolites and metabolic pathways associated with one or multiple pesticide classes, including mitochondrial energy metabolism, fatty acid and lipid metabolism, and amino acid metabolism. Utilizing an integrative network approach, we found that the fatty acid β-oxidation pathway is a common pathway shared across all three pesticide classes. The disruptions of the serum metabolome suggested that chronic pesticide exposure might result in oxidative stress, inflammatory reactions, and mitochondrial dysfunction, all of which have been previously implicated in a wide variety of diseases. Overall, our findings provided a comprehensive view of the molecular mechanisms of chronic pesticide toxicity, and, for the first time, our approach informs exposome research by moving from macrolevel population exposures to microlevel biologic responses.

Although Parkinson's Disease (PD) is typically described in terms of motor symptoms, depression is a common feature. We explored whether depression influences blood-based genome-wide DNA methylation (DNAm) in 692 subjects from a population-based PD case-control study, using both a history of clinically diagnosed depression and current depressive symptoms measured by the geriatric depression scale (GDS). While PD patients in general had more immune activation and more accelerated epigenetic immune system aging than controls, the patients experiencing current depressive symptoms (GDS≥5) showed even higher levels of both markers than patients without current depressive symptoms (GDS<5). For PD patients with a history of clinical depression compared to those without, we found no differences in immune cell composition. However, a history of clinical depression among patients was associated with differentially methylated CpGs. Epigenome-wide association analysis (EWAS) revealed 35 CpGs associated at an FDR≤0.05 (569 CpGs at FDR≤0.10, 1718 CpGs at FDR≤0.15). Gene set enrichment analysis implicated immune system pathways, including immunoregulatory interactions between lymphoid and non-lymphoid cells (p-adj = 0.003) and cytokine-cytokine receptor interaction (p-adj = 0.004). Based on functional genomics, 25 (71%) of the FDR≤0.05 CpGs were associated with genetic variation at 45 different methylation quantitative trait loci (meQTL). Twenty-six of the meQTLs were also expression QTLs (eQTLs) associated with the abundance of 53 transcripts in blood and 22 transcripts in brain (substantia nigra, putamen basal ganglia, or frontal cortex). Notably, cg15199181 was strongly related to rs823114 (SNP-CpG p-value = 3.27E-310), a SNP identified in a PD meta-GWAS and related to differential expression of PM20D1, RAB29, SLC41A1, and NUCKS1. The entire set of genes detected through functional genomics was most strongly overrepresented for interferon-gamma-mediated signaling pathway (enrichment ratio = 18.8, FDR = 4.4e-03) and T cell receptor signaling pathway (enrichment ratio = 13.2, FDR = 4.4e-03). Overall, the current study provides evidence of immune system involvement in depression among Parkinson's patients.

Background

Stochastic epigenetic mutations (SEM) reflect a deviation from normal site-specific methylation patterns. Epigenetic mutation load (EML) captures the accumulation of SEMs across an individual's genome and may reflect dysfunction of the epigenetic maintenance system in response to epigenetic challenges.

Objective

We investigate whether EML is associated with PD risk and time to events (i.e., death and motor symptom decline).

Methods

We employed logistic regression and Cox proportional hazards regression to assess the association between EML and several outcomes. Our analyses are based on 568 PD patients and 238 controls from the Parkinson's disease, Environment and Genes (PEG) study, for whom blood-based methylation data was available.

Results

We found an association for PD onset and EML in all genes (OR = 1.90; 95%CI 1.52-2.37) and PD-related genes (OR = 1.87; 95%CI 1.50-2.32). EML was also associated with time to a minimum score of 35 points on the motor UPDRS exam (OR = 1.28; 95%CI 1.06-1.56) and time to death (OR = 1.29, 95%CI 1.11-1.49). An analysis of PD related genes only revealed five intragenic hotspots of high SEM density associated with PD risk.

Conclusion

Our findings suggest an enrichment of methylation dysregulation in PD patients in general and specifically in five PD related genes. EML may also be associated with time to death and motor symptom progression in PD patients.

Epigenetic changes may be biomarkers of health. Epigenetic age acceleration (EAA), the discrepancy between epigenetic age measured via epigenetic clocks and chronological age, is associated with morbidity and mortality. However, the intersection of epigenetic clocks with microRNAs (miRNAs) and corresponding miRNA-based health implications have not been evaluated. We analyzed DNA methylation and miRNA profiles from blood sampled among 332 individuals enrolled across 2 U.S.-based firefighter occupational studies (2015-2018 and 2018-2020). We considered 7 measures of EAA in leukocytes (PhenoAge, GrimAge, Horvath, skin-blood, and Hannum epigenetic clocks, and extrinsic and intrinsic epigenetic age acceleration). We identified miRNAs associated with EAA using individual linear regression models, adjusted for sex, race/ethnicity, chronological age, and cell type estimates, and investigated downstream effects of associated miRNAs with miRNA enrichment analyses and genomic annotations. On average, participants were 38 years old, 88% male, and 75% non-Hispanic white. We identified 183 of 798 miRNAs associated with EAA (FDR q < 0.05); 126 with PhenoAge, 59 with GrimAge, 1 with Horvath, and 1 with the skin-blood clock. Among miRNAs associated with Horvath and GrimAge, there were 61 significantly enriched disease annotations including age-related metabolic and cardiovascular conditions and several cancers. Enriched pathways included those related to proteins and protein modification. We identified miRNAs associated with EAA of multiple epigenetic clocks. PhenoAge had more associations with individual miRNAs, but GrimAge and Horvath had greater implications for miRNA-associated pathways. Understanding the relationship between these epigenetic markers could contribute to our understanding of the molecular underpinnings of aging and aging-related diseases.

Parkinson's disease (PD) is a complex neurodegenerative disease with etiology rooted in genetic vulnerability and environmental factors. Here we combine quantitative epidemiologic study of pesticide exposures and PD with toxicity screening in dopaminergic neurons derived from PD patient induced pluripotent stem cells (iPSCs) to identify Parkinson's-relevant pesticides. Agricultural records enable investigation of 288 specific pesticides and PD risk in a comprehensive, pesticide-wide association study. We associate long-term exposure to 53 pesticides with PD and identify co-exposure profiles. We then employ a live-cell imaging screening paradigm exposing dopaminergic neurons to 39 PD-associated pesticides. We find that 10 pesticides are directly toxic to these neurons. Further, we analyze pesticides typically used in combinations in cotton farming, demonstrating that co-exposures result in greater toxicity than any single pesticide. We find trifluralin is a driver of toxicity to dopaminergic neurons and leads to mitochondrial dysfunction. Our paradigm may prove useful to mechanistically dissect pesticide exposures implicated in PD risk and guide agricultural policy.

Pesticides have been associated with Parkinson's disease (PD), and protective gloves and workplace hygiene can reduce pesticide exposure. We assessed whether use of gloves and workplace hygiene modified associations between pesticides and PD. The Farming and Movement Evaluation (FAME) study is a nested case-control study within the Agricultural Health Study. Use of protective gloves, other PPE, and hygiene practices were determined by questionnaire (69 cases and 237 controls were included). We considered interactions of gloves and hygiene with ever-use of pesticides for all pesticides with ≥5 exposed and unexposed cases and controls in each glove-use stratum (paraquat, permethrin, rotenone, and trifluralin). 61% of respondents consistently used protective gloves and 87% consistently used ≥2 hygiene practices. Protective glove use modified the associations of paraquat and permethrin with PD: neither pesticide was associated with PD among protective glove users, while both pesticides were associated with PD among non-users (paraquat OR 3.9 [95% CI 1.3, 11.7], interaction p=0.15; permethrin OR 4.3 [95% CI 1.2, 15.6] interaction p=0.05). Rotenone was associated with PD regardless of glove use. Trifluralin was associated with PD among participants who used <2 hygiene practices (OR 5.5 [95% CI 1.1, 27.1]) but was not associated with PD among participants who used 2 or more practices (interaction p=0.02). Although sample size was limited in the FAME study, protective glove use and hygiene practices appeared to be important modifiers of the association between pesticides and PD and may reduce risk of PD associated with certain pesticides.

Objective

The aim of the study is to examine associations between years of firefighting service and eight chronological age-adjusted measures of blood leukocyte epigenetic age acceleration: Horvath, Hannum, SkinBloodClock, Intrinsic, Extrinsic, PhenoAge, GrimAge, and DNAm telomere length.

Methods

The study used a repeated measures analysis of data from 379 incumbent firefighters from eight career departments and 100 recruit firefighters from two of the departments, across the United States.

Results

Incumbent firefighters had on average greater epigenetic age acceleration compared with recruit firefighters, potentially due to the cumulative effect of occupational exposures. However, among incumbent firefighters, additional years of service were associated with epigenetic age deceleration, particularly for GrimAge, a strong predictor of mortality.

Conclusions

Long-term studies with more specific occupational exposure classification are needed to better understand the relationship between years of service and aging biomarkers.