Association of Dual LRRK2 G2019S and GBA Variations With Parkinson Disease Progression

Key Points

Question  What are the associations of concurrent LRRK2 G2019S and GBA variations with clinical progression of Parkinson disease (PD)?

Findings  In this cohort study combining data for 1193 participants with PD from multiple studies, individuals with dual LRRK2 G2019S and GBA variation PD had a slower rate of cognitive decline than those with GBA PD alone, and this was not different from individuals with LRRK2 G2019S PD alone, supporting the notion that there is a dominant association of the LRRK2 gene in individuals with both variations. There was also a novel statistical interaction between LRRK2 G2019S and GBA variations in cognitive decline.

Meaning  These findings suggest that there was not a convergent deleterious association of LRRK2 and GBA variations in PD progression, as would be expected based on prior cellular studies.

Importance  Despite a hypothesis that harboring a leucine-rich repeat kinase 2(LRRK2) G2019S variation and a glucocerebrosidase (GBA) variant would have a combined deleterious association with disease pathogenesis, milder clinical phenotypes have been reported in dual LRRK2 and GBA variations Parkinson disease (PD) than in GBA variation PD alone.

Objective  To evaluate the association of LRRK2 G2019S and GBA variants with longitudinal cognitive and motor decline in PD.

Design, Setting, and Participants  This longitudinal cohort study of continuous measures in LRRK2 PD, GBA PD, LRRK2/GBA PD, and wild-type idiopathic PD used pooled annual visit data ranging from 2004 to 2019 from the Mount Sinai Beth Israel, Parkinson Disease Biomarker Program, Harvard Biomarkers Study, Ashkenazi Jewish-LRRK2-Consortium, Parkinson Progression Marker Initiative, and SPOT-PD studies. Patients who were screened for GBA and LRRK2 variations and completed either a motor or cognitive assessment were included. Data were analyzed from May to July 2020.

Main Outcomes and Measures  The associations of LRRK2 G2019S and GBA genotypes on the rate of decline in Montreal Cognitive Assessment (MoCA) and Movement Disorders Society-Unified Parkinson Disease Rating Scale–Part III scores were examined using linear mixed effects models with PD duration as the time scale.

Results  Among 1193 individuals with PD (mean [SD] age, 66.6 [9.9] years; 490 [41.2%] women), 128 (10.7%) had GBA PD, 155 (13.0%) had LRRK2 PD, 21 (1.8%) had LRRK2/GBA PD, and 889 (74.5%) had idiopathic PD. Patients with GBA PD had faster decline in MoCA than those with LRRK2/GBA PD (B [SE], −0.31 [0.09] points/y; P < .001), LRRK2 PD (B [SE], −0.33 [0.09] points/y; P < .001), or idiopathic PD (B [SE], −0.23 [0.08] points/y; P = .005). There was a LRRK2 G2019S × GBA interaction in MoCA decline (B [SE], 0.22 [0.11] points/y; P = .04), but not after excluding severe GBA variations (B [SE], 0.12 [0.11] points/y; P = .28). Patients with GBA PD had significantly worse motor progression compared with those with idiopathic PD (B [SE], 0.49 [0.22] points/y; P = .03) or LRRK2 PD (B [SE], 0.77 [0.26] points/y; P = .004).

Conclusions and Relevance  These findings suggest that longitudinal cognitive decline in patients with GBA PD was more severe than in those with LRRK2/GBA PD, which more closely resembled LRRK2 PD. This further supports the notion of a dominant association of LRRK2 on GBA in individuals who carry both and raises the possibility of an LRRK2 × GBA interaction. However, the biological basis of a dominant association or interaction is not clear and is apparently contrary to basic investigations. Study of a larger cohort of individuals with severe GBA variation is warranted.

Variants in the leucine-rich repeat kinase 2 (LRRK2; OMIM 609007) and glucocerebrosidase (GBA; OMIM 606463) genes are frequent and worldwide genetic contributors to Parkinson disease (PD) susceptibility.^1,2 LRRK2 variations occur in approximately 2% of PD overall, and up to 40% of PD in certain ethnic groups, particularly individuals of Ashkenazi Jewish or Arabian Berber descent.¹ GBA variation and risk variants in GBA may be present in as many of 10% of PD worldwide and in approximately 15% or greater in PD among individuals of Ashkenazi Jewish descent.^2-4 Wide ranges in conferred disease risk exist based on allelic heterogeneity. The penetrance of LRRK2 variations (specifically G2019S), or the chance that a carrier of this variation will manifest PD by late adulthood, has been estimated at 26% to 43% and is lower in individuals with GBA variations, ranging from 9% to 19%.^5-8 Therefore, additional risk factors are clearly at play.^9-11 One potential critical interaction is the association of GBA and LRRK2 variations in carriers of both variations (hereafter, LRRKR2/GBA).

Several lines of basic experimental data demonstrate convergent deleterious associations of LRRK2 and GBA variations.^12,13 These are consistent with limited clinical reports suggesting earlier age at onset of PD and increased penetrance in carriers of LRRK2/GBA.^14,15 However, cross-sectional analyses of patients with established PD with both LRRK2 and GBA variations suggest that carriers of LRRK2/GBA may instead have a less severe clinical course than those harboring GBA variations only,^15,16 including less severe motor symptoms, and less dementia, and less severe nonmotor features in LRRK2/GBA PD compared with GBA PD.^15,16 Because there are founder variations for LRRK2 G2019S and GBA variations in individuals of Ashkenazi Jewish descent, the likelihood of identifying carriers of LRRK2/GBA is increased in this ethnic group. We expanded on these studies by evaluating longitudinal PD progression, with particular focus on continuous cognitive and motor scores in individuals with LRRK2/GBA PD compared with those with LRRK2 PD, GBA PD, and idiopathic PD (ie, individuals without GBA or LRRK2 variations), and assessing the LRRK2 × GBA gene interaction on cognitive and motor decline. In addition, we extended the sample to include cohorts with individuals not of Ashkenazi Jewish descent as well.

Ethical approval for the Mount Sinai Beth Israel (MSBI) Biomarker study was received from the Mount Sinai institutional review board, and all other sites and studies had local institutional review board approval. All participants provided written informed consent. This study is reported following the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guideline for cohort studies.

To maximize the number of dual LRRK2 G2019S and GBA variant carriers with PD, we combined data from subgroups of prospective longitudinal cohorts of participants with or without Ashkenazi Jewish ancestry and included LRRK2/GBA PD and compared these with LRRK2 PD, GBA PD, and idiopathic PD.

Data from MSBI Biomarker study, National Institute of Neurological Disorders and Stroke (NINDS) Parkinson Disease Biomarker Program (PDBP),¹⁷ the Harvard Biomarkers Study (HBS),^18-20 the Michael J. Fox Foundation (MJFF), LRRK2 Ashkenazi Jewish Cohort Consortium (AJLCC, including Columbia, Tel Aviv Sourasky Medical Center, and MSBI),²¹ the SPOT study,²² and Parkinson Progression Marker Initiative (PPMI)²³ were included. Study visits were performed between February 25, 2004, and December 31, 2019, depending on the study. The goal was to include the maximum number of individuals with LRRK2/GBA PD. Participants were screened for LRRK2 and GBA variations and were included if they completed either a motor rating (Movement Disorder Society–sponsored revision of the Unified Parkinson Disease Rating Scale–Part III [MDS-UPDRS-III],²⁴ or the original unrevised version [UPDRS-III]²⁵) or a cognitive task (Montreal Cognitive Assessment [MoCA]²⁶ or Mini-Mental State Examination [MMSE]²⁷), with demographic information available. MSBI and AJLCC participants were genotyped for LRRK2 G2019S and the 11 most common GBA variations among individuals of Ashkenazi Jewish descent (ie, N370S, 84GG, IVS2 + 1, V394L, D409G, L444P, A456P, RecNcil, R496H, E326K, and T369M). The PDBP and PPMI studies genotyped participants via NeuroX-derived genotyping or the Immunochip Array, including GBA T369M, E326K, N370S, F298L, C62Y, K13R, F255Y, E150K, E427K, D419N, or A309V, and the HBS study identified variations through targeted next-generation sequencing as previously described.^28,29 PPMI genetic data were queried May 2020.

Our comparison groups focused on the subset of LRRK2 PD previously described for progression (AJLCC, including follow-up at MSBI), and the subset of GBA PD not previously reported from the PDBP and MSBI studies (GBA PD and idiopathic PD). To prevent any double representation of participant enrollment, only unique data from a single cohort were included if a participant was enrolled in multiple cohorts. A total of 534 individuals with LRRK2 or idiopathic PD from MSBI and the AJ LCC cohort were also included in prior description of longitudinal progression,³⁰ and 3 participants with LRRK2/GBA PD in the PPMI study were previously reported.¹⁶ GBA variants were categorized³¹ into 4 variation severity categories: severe (including 84GG, L444P, L444R, R120W, RecNcil, V394L, and biallelic), mild (including N370S and R496H), variant (including E326K and T369M), and wild type.

Criteria for PD and clinical, demographic, motor, and nonmotor data were collected per each study’s protocol, as previously described^17,18,21-23 (Table 1). Longitudinal visits occurred a mean (SD) of 12 (3) months apart. Data were harmonized according to published conversions, with discordant scales converted to a single scale for analysis: the UPDRS-III³² subsection score was converted to the MDS-UPDRS III^24,25 for the MSBI and AJLCC cohorts, and continuous MMSE²⁷ scores were converted to MoCA²⁶ scores for the HBS cohort as previously described.³³

The primary goals were the evaluation of the differences in the longitudinal rate of change in cognition (MoCA) and motor function (MDS-UPDRS III) among idiopathic PD, LRRK2 PD, GBA PD (all variant types included in the main analysis), and LRRK2/GBA PD. Linear mixed effects models with robust variance estimates using PD duration (from age of onset) as the time scale were used to examine the associations of LRRK2 and GBA genotypes on the rate of decline in MoCA and MDS-UPDRS-III scores. The primary models included LRRK2 G2019S and GBA genotypes and their interaction term as the primary independent variables. Cognitive and motor progressions were modeled adjusting for sex, age at baseline, PD duration at baseline, and cohort as fixed effects. Models included a participant-specific random effect to account for the correlation in repeated measurements within the same participant. In addition, motor progression models adjusted for levodopa equivalent dose and for presence of deep brain stimulation implants at time of visit. All analyses were performed using Stata statistical software version 16 (StataCorp). Because there is heterogeneity among GBA variants,^31,34 sensitivity analyses excluding severe and biallelic GBA variation carriers were also performed. Individual waves with missing data were excluded. All tests were 2-sided, and P < .05 was considered significant. Data were analyzed from May to July 2020.

Among 1193 participants with PD (mean [SD] age, 66.6 [9.9] years; 490 [41.2%] women) included, 128 (10.7%) had GBA PD, 155 (13.0%) had LRRK2 PD, 21 (1.8%) had LRRK2/GBA PD, and 889 (74.5%) had idiopathic PD. The mean (SD) duration of study participation was 2.85 (1.43) years.

Site information and baseline clinical and genotype characteristics are presented in Table 1. Participants with LRRK2 PD were older at baseline assessment (mean [SD] age, 68.4 [9.2] years) compared with individuals with idiopathic PD (mean [SD] age, 66.5 [10.0] years; P = .03) or GBA PD (mean [SD] age, 64.8 [9.7] years; P < .001) but not LRRK2/GBA PD (mean [SD] age, 65.7 [9.0] years; P = .18). Similarly, participants with LRRK2 PD had longer mean (SD) duration of motor symptoms at baseline assessment (8.3 [6.5] years) compared with participants with idiopathic PD (5.5 [4.9] years; P < .001) or GBA PD (6.9 [5.9] years; P = .046) but not LRRK2/GBA PD (6.0 [6.7] years; P = .13). Participants with GBA PD were younger at motor symptoms onset (mean [SD] age, 57.5 [10.2] years) than those with LRRK2 PD (mean [SD] age, 60.1 [9.6] years), LRRK2/GBA PD (mean [SD] age, 59.6 [9.9]), or idiopathic PD (mean [SD] age, 61.1 [10.0] years). Idiopathic PD was less common in women (342 women of 889 participants [38.5%]) compared with GBA PD (58 women of 128 participants [45.3%]), LRRK2 PD (78 women of 155 participants [50.3%]), and LRRK2/GBA PD (78 women of 155 participants [57.1%]). There were no baseline differences in MDS-UPDRS III or MoCA scores among groups.

Participants with GBA PD were less likely to have mild variations and more likely to have high-risk variations than participants with LRRK2/GBA PD (biallelic: 8 participants [6.3%] vs 1 participant [4.8%]; severe: 11 participants [8.6%] vs 2 participants [9.5%]; mild: 49 participants [38.3%] vs 16 participants [76.2%]; variant: 60 participants [46.9%] vs 2 participants [9.5%]; overall P = .01) (Table 2); however, when mild or high-risk variant and severe or biallelic variations were combined, the distribution in GBA PD did not differ from LRRK2/GBA PD (109 participants [85.2%] vs 18 participants [85.7%]; P = .62). Information regarding deep brain stimulation was collected in all but 1 cohort, and included 6 participants with idiopathic PD (0.7%), 6 participants with GBA PD (4.7%), 5 participants with LRRK2 PD (3.2%), and 2 participants with LRRK2/GBA PD (9.5%).

The estimated (SE) rates of decline in total MoCA score were −0.29 (0.05) points/y for idiopathic PD, −0.52 (0.09) points/y for GBA PD, −0.19 (0.06) points/y for LRRK2 PD, and −0.21 (0.06) points/y for LRRK2/GBA PD (Figure 1). There was a slower, although not statistically significant, rate of cognitive decline for carriers with a LRRK2 G2019S variation compared with individuals with idiopathic PD for LRRK2 PD alone (B [SE], 0.10 [0.06] points/y; P = .08) and LRRK2/GBA PD (B [SE], 0.08 [0.05] points/y; P = .12). As anticipated, there was a worse rate of cognitive decline in participants with GBA PD compared with participants with idiopathic PD (B [SE], −0.23 [0.08] points/y; P = .005) and LRRK2 PD (B [SE], −0.33 [0.09] points/y; P < .001), as well as those with LRRK2/GBA PD (B [SE], −0.31 [0.09] points/y; P < .001). However, there was no difference in rate of cognitive change between participants with LRRK2 PD and those with LRRK2/GBA PD (B [SE], 0.01 [0.07] points/y; P = .85).

There was a significant GBA × LRRK2 G2019S interaction in decline in cognitive scores (B [SE], 0.22 [0.11] points/y; P = .04). This interaction is described in detail in the eTable and eFigure in the Supplement. Specifically, among participants without the LRRK2 G2019S variation, participants with GBA experienced faster decline compared with those without (ie, GBA PD vs idiopathic PD), but among participants with LRRK2 G2019S, those with a GBA variation had no difference compared with participants without a GBA variation (ie, LRRK2/GBA PD vs LRRK2 PD). Similarly, among participants without a GBA variation, carrying a LRRK2 G2019S variation was associated with a slower cognitive decline (LRRK2 PD vs idiopathic PD) as well as among individuals with LRRK2/GBA variations (LRRK2/GBA PD vs GBA PD). In sensitivity analysis excluding participants with severe GBA and dual variations the direction was the same, but the interaction was no longer significant (B [SE], 0.12 [0.11] points/y; P = .28).

The estimated (SE) rates of worsening in MDS-UPDRS III scores were 1.15 (0.16) points/y for idiopathic PD, 1.64 (0.26) points/y for GBA PD, 0.87 ± 0.20) points/y for LRRK2 PD, and 1.00 (0.35) points/y for LRRK2/GBA PD. There was no significant difference in rate of motor change between participants with LRRK2 PD and those with idiopathic PD or LRRK2/GBA PD, nor between participants with LRRK2/GBA PD and those with idiopathic PD. Participants with GBA PD had worse motor progression compared with those with idiopathic PD (estimate [SE], 0.49 [0.22] points/y; P = .03) and those with LRRK2 PD (estimate [SE] 0.77 [0.26] points/y; P = .004), and a not statistically significant difference from those with LRRK2/GBA PD (estimate [SE], 0.64 [0.36] points/y; P = .07). There was no significant interaction association between LRRK2 G2019S and GBA variations for motor progression in the main analysis (estimate [SE], −0.36 [0.40] points/y; P = .38) (Table 3).

In this cohort study systematically evaluating longitudinal cognitive scores, LRRK2 G2019S and GBA variations did not have a combined deleterious association with disease progression in participants with both mutations. Cognitive decline in participants with LRRK2/GBA PD was less severe than in those with GBA PD and more closely resembled that observed in participants with LRRK2 PD. Our findings are overall consistent with those of 2 major studies evaluating LRRK2/GBA^15,16 and are in alignment by finding that nonmotor anomalies were less prominent in individuals with LRRK2/GBA variations than in those with GBA variations. In a 2019 study by Yahalom et al in an Ashkenazi Jewish cohort,¹⁵ individuals with LRRK2/GBA PD had significantly less cognitive impairment, more rapid eye movement sleep behavior disorder, and less dementia and psychosis than those with GBA PD. In a larger and separate Israeli study evaluating individuals with LRRK2/GBA PD,¹⁶ there were no cross-sectional differences in MoCA performance, but there was less depression in participants with LRRK2/GBA PD than among those with severe variation GBA PD and more preserved olfaction in participants with LRRK2/GBA PD, as well as greater motor function than in those with GBA PD. Omer et al¹⁶ suggested a dominant association of the LRRK2 variation in participants with variations in LRRK2 and GBA. Our work extends these prior studies^15,16 through longitudinal evaluation of continuous measures of cognition and motor function and with inclusion of participants with and without Ashkenazi ancestry. We specifically demonstrated that the longitudinal cognitive decline in individuals with LRRK2/GBA PD, as measured by MoCA, was more similar to that observed in individuals with LRRK2 PD, consistent with the cross-sectional data from the study by Yahalom et al.¹⁵

Our study also expands prior work by identifying a novel statistical interaction between LRRK2 G2019S variations and GBA variations in cognitive decline. However, we cannot be certain that the statistical interaction translates to a biological one. Genetic heterogeneity, including GBA allelic differences, as well as differences in GBA and LRRK2 penetrance, clinical expression, and distribution of types of variations in the GBA and LRRK2/GBA groups may all play roles in explaining the statistical interaction. Sensitivity analyses excluding individuals with severe GBA variations no longer demonstrated significance in the interaction association for cognitive decline. However, it should be noted that the power to detect an interaction was subsequently reduced in this restricted sensitivity sample. Furthermore, although the magnitude of the association was reduced, the direction of the association was maintained, suggesting that the overall association of GBA is still diminished in those with the LRRK2/GBA variation. In a 2020 cross-sectional evaluation of the largest sample of individuals with LRRK2/GBA dual variation,¹⁶ the dominant association of LRRK2 G2019S persisted with analysis restricted to individuals with mild GBA variations.

As anticipated from prior studies of single variation data,^35,36 we found that harboring a GBA variation was associated with faster motor deterioration. While the slope of motor decline was worse in participants with GBA PD than in those with LRRK2/GBAPD, the difference was not statistically significant. The difference between participants with LRRK2 PD and those with idiopathic PD also did not reach statistical significance, but the magnitude and direction of the difference was similar to our prior work.³⁰

Our overall finding, that cognition was not worse in participants with LRRK2/GBA PD than in those with only the LRRK2 variation, challenges cellular and clinical data suggesting convergent deleterious associations of dual LRRK2 and GBA variations. Studies in GBA1 variant mouse astrocytes¹³ and in human LRRK2 variation–induced pluripotent stem cell¹² derived dopamine neurons support a convergent deleterious association of GBA and LRRK2 variations on some phenotypes; that is, inhibiting the overactivity of LRRK2 kinase improved function in the lysosomal glucocerebrosidase GBA pathway and ameliorated GBA variation–mediated cell loss. However, these findings are in contrast to the elevated peripheral glucocerebrosidase activity observed from dried blood spots from individuals with LRRK2 G2019S PD.²²

Furthermore, not all studies that evaluated PD penetrance and age at onset among individuals of Ashkenazi Jewish descent with PD suggested a combined deleterious association in carriers of both LRRK2 G2019S and GBA variations or risk variants. Some have demonstrated a higher odds ratio of developing PD, as well as an earlier age of PD onset, than in dual carriers vs those of only a LRRK2 G2019S or GBA variant.^14,15 If these penetrance and age at onset findings are replicated, one explanation that may reconcile these with the finding of less severe nonmotor progression in carriers of dual LRRK2/GBA variations is that LRRK2 and GBA may differentially be associated with the development of PD (penetrance and onset) and the PD disease course (clinical course). We speculate that once PD is established, the joint association might occur in 1 of 2 ways. On the one hand, individuals with dual LRRK2/GBA variations might be more likely to express LRRK2-associated pathologic and clinical expression when LRRK2 G2019S variations and associated gain-of-function renders otherwise deleterious GBA variations inert with respect to GBA progression phenotypes (Figure 2). Alternatively, because most GBA variations have significantly lower penetrance than LRRK2 G2019S, GBA variations in these individuals with LRRK2/GBA variations may be nonpenetrant carriers, and the clinical course may be driven by the LRRK2 variation, with the GBA variation not demonstrating great association, particularly for mild or high-risk variations. However, the biological underpinnings of these possibilities are not clear, and further interrogation of biological mechanisms, including biomarker analysis in a larger sample of individuals with PD harboring dual LRRK2 and severe GBA variations, and an analysis of the consequence of LRRK2 and/or GBA variations on such a biomarker is needed.

A strength of our study is that we have pooled and analyzed standardized, prospective data with up to 7 years of follow-up. Unlike previous studies that have published data on individuals with LRRK2/GBA, our analysis allowed longitudinal assessment to estimate the rate of motor and cognitive progression in LRRK2/GBA PD. Furthermore, to increase the sample size and enrich the genetic cohorts, our study included data collected from multiple cohorts and increased the proportion and variation of GBA variation types represented. However, because Ashkenazi Jewish ethnicity was not queried in 2 sites, we could not stratify the entire study. In sensitivity analysis limited to individuals with known Ashkenazi Jewish ethnicities, the slope of cognitive decline in LRRK2/GBA PD remained less than GBA PD.

Our study has some weaknesses, such as the noted genetic heterogeneity across sites and variation groups, including the few individuals with dual LRRK2/GBA variation and severe GBA variations, and the small sample size of individuals with LRRK2/GBA PD, which limits generalizability. As we did not have equivalent GBA sequencing data for all sites, some individuals with GBA variations were likely miscoded as having idiopathic PD. However, this would bias toward the null rather than toward our finding in cognition. Not all data were available or used from all cohorts. For example, HBS data were limited to those associated with the MSBI study, and to avoid any potential overlap, we limited PPMI data to LRRK2/GBA PD and did not include the entire PPMI cohort data. In sensitivity analyses excluding the 7 individuals with dual LRRK2/GBA variation from PPMI, all findings, including group differences in the rate of change in MoCA score and MDS-UPDRS III rating, and an interaction association in MoCA were maintained. Additionally, there was a need to harmonize measures across sites, and furthermore, MoCA scores may be lower in non-English versions.³⁷

The findings of this cohort study support prior, primarily cross-sectional data, suggesting that individuals with dual LRRK2/GBA variations do not have a worse PD clinical course than individuals with GBA PD or LRRK2 PD alone. We also found a statistical interaction between LRRK2 G2019S and GBA that may point to a possible biological association of LRRK2 G2019S canceling the deleterious associations of GBA. However, because the association was no longer significant in sensitivity analyses excluding individuals with severe or dual variations, we cannot exclude that the interaction may be a statistical reflection of genetic heterogeneity in this group, and that the differential associations of GBA variations and the interaction’s biological relevance are not certain. Because it was challenging to identify individuals that harbor both GBA and LRRK2 G2019S variations and we combined diverse groups, further investigation with systematic larger groups, especially including individuals with severe GBA variations, is warranted. This will facilitate better understanding of putative interaction mechanisms and potential therapeutic implications for cognitive progression.

Accepted for Publication: February 20, 2021.

Published: April 21, 2021. doi:10.1001/jamanetworkopen.2021.5845

Open Access: This is an open access article distributed under the terms of the CC-BY License. © 2021 Ortega RA et al. JAMA Network Open.

Corresponding Authors: Roberto A. Ortega, MS (roberto.ortega@mountsinai.org); Rachel Saunders-Pullman, MD, MPH (rachel.saunders-pullman@mountsinai.org), Department of Neurology, Mount Sinai Beth Israel and Icahn School of Medicine at Mount Sinai, 10 Union Square E, Ste 5J, New York, NY 10003.

Author Contributions: Drs Wang and Saunders-Pullman had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.

Concept and design: Giladi, Bressman, Saunders-Pullman.

Acquisition, analysis, or interpretation of data: Ortega, Wang, Raymond, Bryant, Scherzer, Thaler, Alcalay, West, Mirelman, Kuras, Marder, Ozelius, Bressman, Saunders-Pullman.

Drafting of the manuscript: Ortega, Wang, West, Bressman, Saunders-Pullman.

Critical revision of the manuscript for important intellectual content: All authors.

Statistical analysis: Ortega, Wang, Saunders-Pullman.

Obtained funding: Scherzer, Alcalay, Marder, Giladi, Ozelius, Saunders-Pullman.

Administrative, technical, or material support: Ortega, Raymond, Scherzer, Thaler, Alcalay, Mirelman, Kuras, Giladi, Bressman, Saunders-Pullman.

Supervision: Scherzer, Bressman, Saunders-Pullman.

Conflict of Interest Disclosures: Dr Wang reported grants from the National Institutes of Health (NIH) outside the submitted work. Dr Thaler reported receiving personal fees from Abbvie outside the submitted work. Dr Alcalay reported receiving personal fees from Sanofi and Janssen outside the submitted work. Dr West reported receiving personal fees from the Michael J. Fox foundation, grants from the Michael J. Fox foundation, Biogen, and Escape Bio, and having a patent for EP2844660B1 outside the submitted work; in addition. Dr Mirelman reported receiving personal fees from Neuroderm Consulting outside the submitted work. Dr Marder reported receiving grants from the Michael J Fox Foundation, the NIH, and Parkinson’s Foundation outside the submitted work. Dr Giladi reported receiving personal fees from Sionara, NeuroDerm, Pharma2B, Denali, Neuron23, Sanofi-Genzyme, Biogen, Abbvie, Lysosomal Therapeutics, Movement Disorder Society and grants from the Michael J Fox Foundation, Parkinson’s Foundation, European Union, Israel Science Foundation, Teva National Network of Excellence program, Biogen, Ionis, Sieratzki Family Foundation, and Aufzien Academic Center in Tel-Aviv University outside the submitted work. Dr Bressman reported grants from Michael J Fox Foundation outside the submitted work. No other disclosures were reported.

Funding/Support: This study was supported by National Institutes of Health (NIH) National Institute of Neurological Disorders and Stroke (NINDS) (grant Nos. U01 NS107016 and U01 NS094148), and the Bigglesworth Family Foundation, Bachman-Strauss Chair.

Role of the Funder/Sponsor: The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Additional Contributions: Data and biospecimens used in preparation of this manuscript were obtained from the Parkinson’s Disease Biomarkers Program (PDBP) Consortium, supported by the National Institute of Neurological Disorders and Stroke at the National Institutes of Health. Investigators in the PDBP Consortium include Roger Albin, Roy Alcalay, Alberto Ascherio, Thomas Beach, Sarah Berman, Bradley Boeve, F. DuBois Bowman, Shu Chen, Alice Chen-Plotkin, William Dauer, Ted Dawson, Paula Desplats, Richard Dewey, Ray Dorsey, Jori Fleisher, Kirk Frey, Douglas Galasko, James Galvin, Dwight German, Lawrence Honig, Xuemei Huang, David Irwin, Kejal Kantarci, Anumantha Kanthasamy, Daniel Kaufer, James Leverenz, Carol Lippa, Irene Litvan, Oscar Lopez, Jian Ma, Lara Mangravite, Karen Marder, Laurie Orzelius, Vladislav Petyuk, Judith Potashkin, Liana Rosenthal, Rachel Saunders-Pullman, Clemens Scherzer, Michael Schwarzschild, Tanya Simuni, Andrew Singleton, David Standaert, Debby Tsuang, David Vaillancourt, David Walt, Andrew West, Cyrus Zabetian, Jing Zhang, and Wenquan Zou. The PDBP Investigators, with exception of Drs Alcalay, Marder, Ozelius, Saunders-Pullman, Scherzer and West have not participated in reviewing the data analysis or content of the manuscript. The Harvard Biomarkers Study (HBS) is a collaborative initiative of Brigham and Women’s Hospital and Massachusetts General Hospital, codirected by Clemens Scherzer and Bradley T. Hyman. HBS investigators are listed at https://www.bwhparkinsoncenter.org. The HBS Investigators have not participated in reviewing the current manuscript. The HBS Study Investigators include codirectors: Clemens R. Scherzer (Brigham and Women’s Hospital) and Bradley T. Hyman (Massachusetts General Hospital); investigators and study coordinators: Yuliya Kuras, Karbi Choudhury, Michael T. Hayes, Aleksandar Videnovic, Nutan Sharma, Vikram Khurana, Claudio Melo De Gusmao, and Reisa Sperling (Brigham and Women’s Hospital); and John H. Growdon, Michael A. Schwarzschild, Albert Y. Hung, Alice W. Flaherty, Deborah Blacker, Anne-Marie Wills, Steven E. Arnold, Ann L. Hunt, Nicte I. Mejia, Anand Viswanathan, Stephen N. Gomperts, Mark W. Albers, Maria Allora-Palli, David Hsu, Alexandra Kimball, Scott McGinnis, John Becker, Randy Buckner, Thomas Byrne, Maura Copeland, Bradford Dickerson, Matthew Frosch, Theresa Gomez-Isla, Steven Greenberg, Julius Hedden, Elizabeth Hedley-Whyte, Keith Johnson, Raymond Kelleher, Aaron Koenig, Maria Marquis-Sayagues, Gad Marshall, Sergi Martinez-Ramirez, Donald McLaren, Olivia Okereke, Elena Ratti, Christopher William, Koene Van Dij, Shuko Takeda, Anat Stemmer-Rachaminov, Jessica Kloppenburg, Catherine Munro, Rachel Schmid, Sarah Wigman, Sara Wlodarcsyk (Massachusetts General Hospital); data coordination: Thomas Yi (Brigham and Women’s Hospital); and biobank management staff: Idil Tuncali (Brigham and Women’s Hospital). Data used in preparation of this study were obtained from the Michael J. Fox Foundation–sponsored LRRK2 Cohort Consortium. The LRRK2 Cohort Consortium is coordinated and funded by The Michael J. Fox Foundation for Parkinson’s Research. Data used in the preparation of this article were obtained from the Parkinson’s Progression Markers Initiative (PPMI) database PPMI, a public-private partnership, is funded by the Michael J. Fox Foundation for Parkinson’s Research and funding partners, and the full names of all of the PPMI funding partners can be found at found at www.ppmi-info.org/fundingpartners.