Identification of novel risk loci, causal insights, and heritable risk for Parkinson's disease: a meta-analysis of genome-wide association studies

Background Genome-wide association studies (GWAS) in Parkinson's disease have increased the scope of biological knowledge about the disease over the past decade. We aimed to use the largest aggregate of GWAS data to identify novel risk loci and gain further insight into the causes of Parkinson's disease. Methods We did a meta-analysis of 17 datasets from Parkinson's disease GWAS available from European ancestry samples to nominate novel loci for disease risk. These datasets incorporated all available data. We then used these data to estimate heritable risk and develop predictive models of this heritability. We also used large gene expression and methylation resources to examine possible functional consequences as well as tissue, cell type, and biological pathway enrichments for the identified risk factors. Additionally, we examined shared genetic risk between Parkinson's disease and other phenotypes of interest via genetic correlations followed by Mendelian randomisation. Findings Between Oct 1, 2017, and Aug 9, 2018, we analysed 7·8 million single nucleotide polymorphisms in 37 688 cases, 18 618 UK Biobank proxy-cases (ie, individuals who do not have Parkinson's disease but have a first degree relative that does), and 1·4 million controls. We identified 90 independent genome-wide significant risk signals across 78 genomic regions, including 38 novel independent risk signals in 37 loci. These 90 variants explained 16–36% of the heritable risk of Parkinson's disease depending on prevalence. Integrating methylation and expression data within a Mendelian randomisation framework identified putatively associated genes at 70 risk signals underlying GWAS loci for follow-up functional studies. Tissue-specific expression enrichment analyses suggested Parkinson's disease loci were heavily brain-enriched, with specific neuronal cell types being implicated from single cell data. We found significant genetic correlations with brain volumes (false discovery rate-adjusted p=0·0035 for intracranial volume, p=0·024 for putamen volume), smoking status (p=0·024), and educational attainment (p=0·038). Mendelian randomisation between cognitive performance and Parkinson's disease risk showed a robust association (p=8·00 × 10−7). Interpretation These data provide the most comprehensive survey of genetic risk within Parkinson's disease to date, to the best of our knowledge, by revealing many additional Parkinson's disease risk loci, providing a biological context for these risk factors, and showing that a considerable genetic component of this disease remains unidentified. These associations derived from European ancestry datasets will need to be followed-up with more diverse data. Funding The National Institute on Aging at the National Institutes of Health (USA), The Michael J Fox Foundation, and The Parkinson's Foundation (see appendix for full list of funding sources)

Identification of candidate Parkinson disease genes by integrating genome-wide association study, expression, and epigenetic data sets

Importance Substantial genome-wide association study (GWAS) work in Parkinson disease (PD) has led to the discovery of an increasing number of loci shown reliably to be associated with increased risk of disease. Improved understanding of the underlying genes and mechanisms at these loci will be key to understanding the pathogenesis of PD. Objective To investigate what genes and genomic processes underlie the risk of sporadic PD. Design and Setting This genetic association study used the bioinformatic tools Coloc and transcriptome-wide association study (TWAS) to integrate PD case-control GWAS data published in 2017 with expression data (from Braineac, the Genotype-Tissue Expression [GTEx], and CommonMind) and methylation data (derived from UK Parkinson brain samples) to uncover putative gene expression and splicing mechanisms associated with PD GWAS signals. Candidate genes were further characterized using cell-type specificity, weighted gene coexpression networks, and weighted protein-protein interaction networks. Main Outcomes and Measures It was hypothesized a priori that some genes underlying PD loci would alter PD risk through changes to expression, splicing, or methylation. Candidate genes are presented whose change in expression, splicing, or methylation are associated with risk of PD as well as the functional pathways and cell types in which these genes have an important role. Results Gene-level analysis of expression revealed 5 genes (WDR6 [OMIM 606031], CD38 [OMIM 107270], GPNMB [OMIM 604368], RAB29 [OMIM 603949], and TMEM163 [OMIM 618978]) that replicated using both Coloc and TWAS analyses in both the GTEx and Braineac expression data sets. A further 6 genes (ZRANB3 [OMIM 615655], PCGF3 [OMIM 617543], NEK1 [OMIM 604588], NUPL2 [NCBI 11097], GALC [OMIM 606890], and CTSB [OMIM 116810]) showed evidence of disease-associated splicing effects. Cell-type specificity analysis revealed that gene expression was overall more prevalent in glial cell types compared with neurons. The weighted gene coexpression performed on the GTEx data set showed that NUPL2 is a key gene in 3 modules implicated in catabolic processes associated with protein ubiquitination and in the ubiquitin-dependent protein catabolic process in the nucleus accumbens, caudate, and putamen. TMEM163 and ZRANB3 were both important in modules in the frontal cortex and caudate, respectively, indicating regulation of signaling and cell communication. Protein interactor analysis and simulations using random networks demonstrated that the candidate genes interact significantly more with known mendelian PD and parkinsonism proteins than would be expected by chance. Conclusions and Relevance Together, these results suggest that several candidate genes and pathways are associated with the findings observed in PD GWAS studies

Peripheral inflammatory immune response differs among sporadic and familial Parkinson’s disease

Peripheral inflammatory immune responses are thought to play a major role in the pathogenesis of Parkinson’s disease (PD). The neutrophil-to-lymphocyte ratio (NLR), a biomarker of systemic inflammation, has been reported to be higher in patients with PD than in healthy controls (HCs). The present study was aimed at determining if the peripheral inflammatory immune response could be influenced by the genetic background of patients with PD. We included a discovery cohort with 222 patients with PD (132 sporadic PD, 44 LRRK2-associated PD (with p.G2019S and p.R1441G variants), and 46 GBA-associated PD), as well as 299 HCs. Demographic and clinical data were recorded. Leukocytes and their subpopulations, and the NLR were measured in peripheral blood. Multivariate lineal regression and post-hoc tests were applied to determine the differences among the groups. Subsequently, a replication study using the Parkinson’s Progression Markers Initiative cohort was performed which included 401 patients with PD (281 sPD patients, 66 LRRK2-PD patients, 54 GBA-PD patients) and a group of 174 HCs. Patients with sporadic PD and GBA-associated PD showed a significantly lower lymphocyte count, a non-significantly higher neutrophil count and a significantly higher NLR than HCs. The peripheral inflammatory immune response of patients with LRRK2-associated PD did not differ from HCs. Our study supports the involvement of a peripheral inflammatory immune response in the pathophysiology of sPD and GBA-associated PD. However, this inflammatory response was not found in LRRK2-associated PD, probably reflecting different pathogenic inflammatory mechanisms

Non-motor symptom burden in patients with Parkinson’s disease with impulse control disorders and compulsive behaviours: results from the COPPADIS cohort

The study was aimed at analysing the frequency of impulse control disorders (ICDs) and compulsive behaviours (CBs) in patients with Parkinson’s disease (PD) and in control subjects (CS) as well as the relationship between ICDs/CBs and motor, nonmotor features and dopaminergic treatment in PD patients. Data came from COPPADIS-2015, an observational, descriptive, nationwide (Spain) study. We used the validated Questionnaire for Impulsive-Compulsive Disorders in Parkinson's Disease-Rating Scale (QUIP-RS) for ICD/CB screening. The association between demographic data and ICDs/CBs was analyzed in both groups. In PD, this relationship was evaluated using clinical features and treatment-related data. As result, 613 PD patients (mean age 62.47 ± 9.09 years, 59.87% men) and 179 CS (mean age 60.84 ± 8.33 years, 47.48% men) were included. ICDs and CBs were more frequent in PD (ICDs 12.7% vs. 1.6%, p < 0.001; CBs 7.18% vs. 1.67%, p = 0.01). PD patients had more frequent previous ICDs history, premorbid impulsive personality and antidepressant treatment (p < 0.05) compared with CS. In PD, patients with ICDs/CBs presented younger age at disease onset, more frequent history of previous ICDs and premorbid personality (p < 0.05), as well as higher comorbidity with nonmotor symptoms, including depression and poor quality of life. Treatment with dopamine agonists increased the risk of ICDs/CBs, being dose dependent (p < 0.05). As conclusions, ICDs and CBs were more frequent in patients with PD than in CS. More nonmotor symptoms were present in patients with PD who had ICDs/CBs compared with those without. Dopamine agonists have a prominent effect on ICDs/CBs, which could be influenced by dose

Epidemiology of ataxia and hereditary spastic paraplegia in Spain: a cross-sectional study

Introducción: Las ataxias (AT) y paraparesias espásticas hereditarias (PEH) son síndromes neurodegenerativos raros. Nos proponemos conocer la prevalencia de las AT y PEH en Espana˜ en 2019. Pacientes y métodos: Estudio transversal, multicéntrico, descriptivo y retrospectivo de los pacientes con AT y PEH, desde marzo de 2018 a diciembre de 2019 en toda Espana. ˜ Resultados: Se obtuvo información de 1933 pacientes procedentes de 11 Comunidades Autónomas, de 47 neurólogos o genetistas. Edad media: 53,64 anos ˜ ± 20,51 desviación estándar (DE); 938 varones (48,5%), 995 mujeres (51,5%). En 920 pacientes (47,6%) no se conoce el defecto genético. Por patologías, 1.371 pacientes (70,9%) diagnosticados de AT, 562 diagnosticados de PEH (29,1%). La prevalencia estimada de AT es 5,48/100.000 habitantes, y la de PEH es 2,24 casos/100.000 habitantes. La AT dominante más frecuente es la SCA3. La AT recesiva más frecuente es la ataxia de Friedreich (FRDA). La PEH dominante más frecuente es la SPG4, y la PEH recesiva más frecuente es la SPG7. Conclusiones: La prevalencia estimada de AT y PEH en nuestra serie es de 7,73 casos/100.000 habitantes. Estas frecuencias son similares a las del resto del mundo. En el 47,6% no se ha conseguido un diagnóstico genético. A pesar de las limitaciones, este estudio puede contribuir a estimar los recursos, visibilizar estas enfermedades, detectar las mutaciones más frecuentes para hacer los screenings por comunidades, y favorecer los ensayos clínicos.Introduction: Ataxia and hereditary spastic paraplegia are rare neurodegenerative syndromes. We aimed to determine the prevalence of these disorders in Spain in 2019. Patients and methods: We conducted a cross-sectional, multicentre, retrospective, descriptive study of patients with ataxia and hereditary spastic paraplegia in Spain between March 2018 and December 2019. Results: We gathered data from a total of 1933 patients from 11 autonomous communities, provided by 47 neurologists or geneticists. Mean (SD) age in our sample was 53.64 (20.51) years; 938 patients were men (48.5%) and 995 were women (51.5%). The genetic defect was unidentified in 920 patients (47.6%). A total of 1371 patients (70.9%) had ataxia and 562 (29.1%) had hereditary spastic paraplegia. Prevalence rates for ataxia and hereditary spastic paraplegia were estimated at 5.48 and 2.24 cases per 100 000 population, respectively. The most frequent type of dominant ataxia in our sample was SCA3, and the most frequent recessive ataxia was Friedreich ataxia. The most frequent type of dominant hereditary spastic paraplegia in our sample was SPG4, and the most frequent recessive type was SPG7. Conclusions: In our sample, the estimated prevalence of ataxia and hereditary spastic paraplegia was 7.73 cases per 100 000 population. This rate is similar to those reported for other countries. Genetic diagnosis was not available in 47.6% of cases. Despite these limitations, our study provides useful data for estimating the necessary healthcare resources for these patients, raising awareness of these diseases, determining the most frequent causal mutations for local screening programmes, and promoting the development of clinical trials

The Genetic Architecture of Parkinson Disease in Spain: Characterizing Population‐Specific Risk, Differential Haplotype Structures, and Providing Etiologic Insight

The Iberian Peninsula stands out as having variable levels of population admixture and isolation, making Spain an interesting setting for studying the genetic architecture of neurodegenerative diseases. To perform the largest PD genome-wide association study restricted to a single country. We performed a GWAS for both risk of PD and age at onset in 7,849 Spanish individuals. Further analyses included population-specific risk haplotype assessments, polygenic risk scoring through machine learning, Mendelian randomization of expression, and methylation data to gain insight into disease-associated loci, heritability estimates, genetic correlations, and burden analyses. We identified a novel population-specific genome-wide association study signal at PARK2 associated with age at onset, which was likely dependent on the c.155delA mutation. We replicated four genome-wide independent signals associated with PD risk, including SNCA, LRRK2, KANSL1/MAPT, and HLA-DQB1. A significant trend for smaller risk haplotypes at known loci was found compared to similar studies of non-Spanish origin. Seventeen PD-related genes showed functional consequence by two-sample Mendelian randomization in expression and methylation data sets. Long runs of homozygosity at 28 known genes/loci were found to be enriched in cases versus controls. Our data demonstrate the utility of the Spanish risk haplotype substructure for future fine-mapping efforts, showing how leveraging unique and diverse population histories can benefit genetic studies of complex diseases. The present study points to PARK2 as a major hallmark of PD etiology in Spain. © 2019 International Parkinson and Movement Disorder Society

The Genetic Architecture of Parkinson Disease in Spain: Characterizing Population-Specific Risk, Differential Haplotype Structures, and Providing Etiologic Insight.

The Iberian Peninsula stands out as having variable levels of population admixture and isolation, making Spain an interesting setting for studying the genetic architecture of neurodegenerative diseases. To perform the largest PD genome-wide association study restricted to a single country. We performed a GWAS for both risk of PD and age at onset in 7,849 Spanish individuals. Further analyses included population-specific risk haplotype assessments, polygenic risk scoring through machine learning, Mendelian randomization of expression, and methylation data to gain insight into disease-associated loci, heritability estimates, genetic correlations, and burden analyses. We identified a novel population-specific genome-wide association study signal at PARK2 associated with age at onset, which was likely dependent on the c.155delA mutation. We replicated four genome-wide independent signals associated with PD risk, including SNCA, LRRK2, KANSL1/MAPT, and HLA-DQB1. A significant trend for smaller risk haplotypes at known loci was found compared to similar studies of non-Spanish origin. Seventeen PD-related genes showed functional consequence by two-sample Mendelian randomization in expression and methylation data sets. Long runs of homozygosity at 28 known genes/loci were found to be enriched in cases versus controls. Our data demonstrate the utility of the Spanish risk haplotype substructure for future fine-mapping efforts, showing how leveraging unique and diverse population histories can benefit genetic studies of complex diseases. The present study points to PARK2 as a major hallmark of PD etiology in Spain. © 2019 International Parkinson and Movement Disorder Society

Identification of novel risk loci, causal insights, and heritable risk for Parkinson's disease: a meta-analysis of genome-wide association studies

BackgroundGenome-wide association studies (GWAS) in Parkinson's disease have increased the scope of biological knowledge about the disease over the past decade. We aimed to use the largest aggregate of GWAS data to identify novel risk loci and gain further insight into the causes of Parkinson's disease.MethodsWe did a meta-analysis of 17 datasets from Parkinson's disease GWAS available from European ancestry samples to nominate novel loci for disease risk. These datasets incorporated all available data. We then used these data to estimate heritable risk and develop predictive models of this heritability. We also used large gene expression and methylation resources to examine possible functional consequences as well as tissue, cell type, and biological pathway enrichments for the identified risk factors. Additionally, we examined shared genetic risk between Parkinson's disease and other phenotypes of interest via genetic correlations followed by Mendelian randomisation.FindingsBetween Oct 1, 2017, and Aug 9, 2018, we analysed 7·8 million single nucleotide polymorphisms in 37 688 cases, 18 618 UK Biobank proxy-cases (ie, individuals who do not have Parkinson's disease but have a first degree relative that does), and 1·4 million controls. We identified 90 independent genome-wide significant risk signals across 78 genomic regions, including 38 novel independent risk signals in 37 loci. These 90 variants explained 16-36% of the heritable risk of Parkinson's disease depending on prevalence. Integrating methylation and expression data within a Mendelian randomisation framework identified putatively associated genes at 70 risk signals underlying GWAS loci for follow-up functional studies. Tissue-specific expression enrichment analyses suggested Parkinson's disease loci were heavily brain-enriched, with specific neuronal cell types being implicated from single cell data. We found significant genetic correlations with brain volumes (false discovery rate-adjusted p=0·0035 for intracranial volume, p=0·024 for putamen volume), smoking status (p=0·024), and educational attainment (p=0·038). Mendelian randomisation between cognitive performance and Parkinson's disease risk showed a robust association (p=8·00 × 10-7).InterpretationThese data provide the most comprehensive survey of genetic risk within Parkinson's disease to date, to the best of our knowledge, by revealing many additional Parkinson's disease risk loci, providing a biological context for these risk factors, and showing that a considerable genetic component of this disease remains unidentified. These associations derived from European ancestry datasets will need to be followed-up with more diverse data.FundingThe National Institute on Aging at the National Institutes of Health (USA), The Michael J Fox Foundation, and The Parkinson's Foundation (see appendix for full list of funding sources)