DNA methylation biomarker for cumulative lead exposure is associated with Parkinson’s disease
Clinical Epigenetics volume 13, Article number: 59 (2021)
Lead, a known neurotoxicant, has previously received attention in Parkinson’s disease (PD) research, but epidemiologic studies have been limited in sample size and findings are equivocal. We generated two methylation-based biomarkers for cumulative tibia and patella bone-measured lead exposure in 1528 PD patients and 1169 controls. PD status was associated with increased levels of the DNAm biomarker for tibia-lead levels. We estimated a meta-OR for PD of 1.89 per unit DNAm tibia-lead increase (95% CI 1.59, 2.24; p = 8.1E−13). The current study supports the notion that chronic and long-term lead exposure tracked via DNAm may contribute to PD pathogenesis.
Parkinson’s disease (PD) is a progressive neurodegenerative disorder characterized by gradual loss of dopaminergic neurons in the substantia nigra. PD etiology is complex and multi-factorial, likely involving a combination of genetic, aging, and environmental components. Lead, a known neurotoxicant, has previously received attention in PD, but findings in epidemiologic studies have been equivocal.
Lead passes through the blood–brain barrier and enters neurons and glia through calcium channels, affecting neurotransmitter release, energy metabolism, and inducing reactive oxygen species (1). Long-term, cumulative lead levels, which depend on both exposure levels and an individual’s uptake, storage, and clearance, may affect the aging brain the most. Interestingly, in epidemiologic studies, results for PD association seem to depend on the type of exposure assessment employed. Specifically, studies using bone-lead measurements found positive associations (2, 3), while those relying on self-report of occupational exposures or emissions data, with one exception (4), did not (5,6,7,8). Suggesting that errors in exposure measurement and issues of timing may obscure the role lead plays in neurodegeneration.
Clearly, bone-lead-measurement is the gold-standard for long-term lead exposure, given that the half-life of lead in the tibia is approximately 20 years (9). However, cost and feasibility have often precluded this option, especially in larger studies. Recent advances in high-throughput epigenetic technology are paving the way for biomarker discovery in many fields, including exposure science. Two methylation-based biomarkers for cumulative lead exposure measured in the tibia and patella have recently been developed (10). Using these biomarkers, for the first time we evaluate the influence of methylation-based cumulative lead exposure on PD risk in two, independent population-based studies with over 2600 participants.
For our analysis, we used genome-wide DNA methylation (DNAm) data from two publicly available PD studies (Gene Expression Omnibus (GEO) accession numbers GSE145361 (SGPD), GSE72774 and GSE72776 (PEG)). DNAm was measured with the Human Methylation 450 k BeadChip in whole blood samples.
The System Genomics of Parkinson’s Disease (SGPD) is a consortium of three studies from across Australia and New Zealand (11). The data available on GEO data consists of 1889 samples (959 PD patients and 930 controls) of European ancestry. Prevalent PD patients were recruited (PD duration 2–40 + years), and controls consisted primarily of community-based age-matched volunteers from the same communities, as well as some patient’s spouses and siblings. Cohort details, DNA extraction methods, quality control procedures, and normalization methods (quantile-normalized and normalization adjusted for batch, slide, cohort, sentrix row/column, sex, and age) have been described (11).
The Parkinson’s Environment and Genes (PEG) study is population-based study from three agricultural counties of Central California (12). GEO data is available for 807 samples (569 PD patients and 238 controls) of European and Hispanic ancestry. Patients early in disease (mean PD duration = 2.9 years (SD = 2.3)) were diagnosed in-person by UCLA Movement Disorder Specialists (J.B.). Population-based controls from the same communities were randomly sampled from Medicare lists and via residential tax assessor's records. Cohort details, DNA extraction methods, quality control procedures, and normalization methods have been previously described (13).
We generated two epigenetic biomarkers for cumulative lead exposure (tibia and patella), developed in the Normative Aging Study (NAS). The epigenetic biosensors of patella and tibia are linear combinations of 59 and 138 CpGs, respectively, identified with site-by-site analysis and combined via machine-learning algorithms trained on K x-ray fluorescence (KXRF) in-vivo measures of bone-lead (10). To determine lead biomarker levels in SGPD and PEG, we extracted the published regression coefficients from NAS (10), and applied them to the corresponding DNAm beta matrices. At time of publication, the specificity of the DNAm biosensors has not been validated beyond the initial study.
In order to quantify the relationship between the DNAm lead and PD, we used logistic regression to estimate odds ratios (ORs) and 95% CIs for PD. The OR is a measure of association that compares the odds of PD based on exposure, with an OR greater than 1 indicating that the odds of PD increase with exposure. Regression diagnostics showed that model assumptions were met. Based on data availability from GEO, we controlled for age (estimated with the Horvath DNAmAge in SGPD (14)), sex, ancestry (PEG only), blood cell composition, smoking history (PEG only), and mean methylation by sample to account for global methylation. We imputed blood cell proportions using the Houseman method (15). We assessed between-study heterogeneity for SGPD and PEG results with Cochran’s Q, and calculated a meta OR using a fixed-effects model with weights based on precision.
Characteristics of SGPD and PEG participants can be found in Table 1. In both studies, PD patients were a few years older than their respective controls and more often men. In PEG, there were less former/current smokers among patients. Additional file 1: Fig. 1 shows box plots of the data, including all data points and the mean values for both DNAm biomarkers by PD status and stratified by study.
PD status was strongly associated with the DNAm biomarker for tibia-lead levels (Table 2) in both cohorts (Cochran’s Q indicated no significant heterogeneity study results, p = 0.25). Specifically, with minimal adjustment (age, sex, ancestry), we estimated in SGPD, an OR for PD of 2.06 per unit DNAm tibia-lead increase (95% CI 1.66, 2.56; p = 5.4E−11) and in PEG of 1.60 (95% CI 1.20, 2.15; p = 1.0E−03), with a meta-OR of 1.89 (95% CI 1.59, 2.24; p = 8.1E−13). Further adjustment for blood cell composition, mean methylation by sample, and smoking status (PEG only), resulted in a meta OR of 1.52 (95% CI 1.25, 1.86; p = 2.8E−5). The DNAm biomarker for patella-lead levels was inversely associated with PD in SGPD (fully adjusted model: OR = 0.70 (95% CI 0.52, 0.93)), but not in PEG (OR = 1.18 (95% CI 0.80, 1.75) (heterogeneity between study results p = 0.006). Results for both DNAm markers stratified by sex can be found in Additional file 3: Table 1. Additionally, results were similar in models both with or without adjustment for smoking (Additional file 3: Table 2).
In two large population-based PD studies, we used DNAm to estimate cumulative, bone-lead levels by applying an externally developed epigenetic (DNAm) biosensor (10), and found that PD was strongly associated with the DNAm pattern estimated for tibia-lead levels. Positive associations were seen independently in both studies and persisted after multiple adjustments. Our findings for DNAm tibia-lead are supported by positive associations in two prior epidemiologic studies that relied on actual tibia bone-lead measurements but had more limited sample sizes (121 and 330 PD patients (2, 3)).
We found no association between DNAm patella-lead and PD in PEG, and even an inverse-association in the SGPD study. Interestingly, one of the previous epidemiologic studies also only observed tibia- but not patella-lead measurements to be associated with PD (2). This is perhaps due to the half-life of lead, which is approximately 8 years in the patella (representing trabecular bone) versus decades in the tibia (representing cortical bone) (2, 16). Thus, tibia-lead measurements represent cumulative exposure from a more distant past than patella. This shorter half-life may also explain why an inverse-association was observed in SGPD, where the patients had a longer PD duration at blood-draw (2–42 years), while PEG patients were enrolled early in disease (mean 2.9 years). Should a primary source of lead exposure among adults be occupationally related, exposure would be less likely to have occurred among patients in the years around PD diagnosis, as patients stop working several years earlier on average than controls due to disease (17). Industries with occupational lead exposure include construction jobs with lead-based paint activities, including removal or demolition, metal and chemical product manufacturing, battery manufacturing and recycling, and electronics manufacturing (18). Other sources of chronic lead exposure, such as drinking water, leaded gasoline, and household paint, however, may be less influenced by PD onset. Overall, multiple sources of evidence from previous and our current study support a positive association for long-term, cumulative lead exposure that extended over many decades, as measured in the tibia, with PD, rather than shorter-term and more recent exposures reflected in patella-lead levels.
The inordinate advantage of our study is the availability of DNAm-lead biomarkers, which also allows for future applications to all PD studies, or more generally studies of age-related disease, with methylation data and larger-scale systematic research into the long-term neurotoxicity of lead. The technology and methodology to track many health-related factors through methylation is rapidly developing. However, currently the DNAm-lead biomarkers have two limitations in our opinion, first, that they were developed in an all-male cohort, and second, that the specificity has not yet been validated beyond the initial study (10). Future validation of the DNAm-lead biomarkers would be incredibly beneficial to the research community, given the widespread adverse health effects from chronic lead exposure. Our study provides a compelling first glance into the utility of DNAm-lead exposure measures in PD and neurodegenerative research in general, and solid evidence for the involvement of this DNAm-lead exposure pattern in PD pathogenesis.
Availability of data and materials
The datasets analyzed during the current study are available in the GEO repository, Accession Numbers GSE145361 (SGPD), GSE72774 and GSE72776 (PEG).
Reuben A. Childhood Lead Exposure and Adult Neurodegenerative Disease. J Alzheimer’s Disease. 2018;64(1):17–42.
Weisskopf MG, Weuve J, Nie H, Saint-Hilaire MH, Sudarsky L, Simon DK, et al. Association of cumulative lead exposure with Parkinson’s disease. Environ Health Perspect. 2010;118(11):1609–13.
Coon S, Stark A, Peterson E, Gloi A, Kortsha G, Pounds J, et al. Whole-body lifetime occupational lead exposure and risk of Parkinson’s disease. Environ Health Perspect. 2006;114(12):1872–6.
Gorell JM, Johnson CC, Rybicki BA, Peterson EL, Kortsha GX, Brown GG, et al. Occupational exposures to metals as risk factors for Parkinson’s disease. Neurology. 1997;48(3):650–8.
Seidler A, Hellenbrand W, Robra BP, Vieregge P, Nischan P, Joerg J, et al. Possible environmental, occupational, and other etiologic factors for Parkinson’s disease: a case-control study in Germany. Neurology. 1996;46(5):1275–84.
Firestone JA, Lundin JI, Powers KM, Smith-Weller T, Franklin GM, Swanson PD, et al. Occupational factors and risk of Parkinson’s disease: a population-based case-control study. Am J Ind Med. 2010;53(3):217–23.
Palacios N, Fitzgerald K, Roberts AL, Hart JE, Weisskopf MG, Schwarzschild MA, et al. A prospective analysis of airborne metal exposures and risk of Parkinson disease in the nurses’ health study cohort. Environ Health Perspect. 2014;122(9):933–8.
Willis AW, Evanoff BA, Lian M, Galarza A, Wegrzyn A, Schootman M, et al. Metal emissions and urban incident Parkinson disease: a community health study of medicare beneficiaries by using geographic information systems. Am J Epidemiol. 2010;172(12):1357–63.
Hu H, Shih R, Rothenberg S, Schwartz BS. The epidemiology of lead toxicity in adults: Measuring dose and consideration of other methodologic issues. Environ Health Perspect. 2007;115(3):455–62.
Colicino E, Just A, Kioumourtzoglou MA, Vokonas P, Cardenas A, Sparrow D, et al. Blood DNA methylation biomarkers of cumulative lead exposure in adults. J Expo Sci Environ Epidemiol. 2019;31(1):108–16.
Vallerga CL, Zhang F, Fowdar J, McRae AF, Qi T, Nabais MF, et al. Analysis of DNA methylation associates the cystine–glutamate antiporter SLC7A11 with risk of Parkinson’s disease. Nat Commun. 2020;
Ritz BR, Paul KC, Bronstein JM. Of pesticides and men: a California story of genes and environment in Parkinson’s disease. Curr Environ Heal reports. 2016;3(1):40–52.
Chuang Y-H, Paul KC, Bronstein JM, Bordelon Y, Horvath S, Ritz B. Parkinson’s disease is associated with DNA methylation levels in human blood and saliva. Genome Med. 2017;9(1):76.
Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013;14(10):96.
Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformat. 2012;
Kim R, Landrigan C, Mossmann P, Sparrow D, Hu H. Age and secular trends in bone lead levels in middle-aged and elderly men: three-year longitudinal follow-up in the normative aging study. Am J Epidemiol. 1997;146(7):586–91.
Koerts J, König M, Tucha L, Tucha O. Working capacity of patients with Parkinson’s disease—a systematic review. Parkinsonism Relat Disord. 2016;27:9–24.
Koh DH, Locke SJ, Chen YC, Purdue MP, Friesen MC. Lead exposure in US worksites: A literature review and development of an occupational lead exposure database from the published literature. Am J Ind Med. 2015;58(6):605–16.
This work was supported by the National Institute of Environmental Health Science (Grant Numbers R21ES024356, 2R01ES010544, U54ES012078).
Ethics approval and consent to participate
The PEG study was approved by the UCLA Institutional Review Board (IRB# 11-001530) and informed consent was obtained from all individuals. Our research conformed to the Declaration of Helsinki.
Consent for publication
Not applicable, this manuscript does not contain any individual person’s data in any form, e.g. individual details, images, or videos.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Box plot figures of the DNAm lead-biomarkers by PD and stratified by study.
Table showing DNAm estimated cumulative lead exposure and PD risk, stratified by sex.
Table showing DNAm estimated cumulative lead exposure and PD risk in PEG, with and without controlling for smoking.
About this article
Cite this article
Paul, K.C., Horvath, S., Del Rosario, I. et al. DNA methylation biomarker for cumulative lead exposure is associated with Parkinson’s disease. Clin Epigenet 13, 59 (2021). https://doi.org/10.1186/s13148-021-01051-3
- Parkinson’s Disease
- DNA methylation