Systematic rare variant analyses identify RAB32 as a susceptibility gene for familial Parkinson's disease.

Despite substantial progress, causal variants are identified only for a minority of familial Parkinson's disease (PD) cases, leaving high-risk pathogenic variants unidentified1,2. To identify such variants, we uniformly processed exome sequencing data of 2,184 index familial PD cases and 69,775 controls. Exome-wide analyses converged on RAB32 as a novel PD gene identifying c.213C > G/p.S71R as a high-risk variant presenting in ~0.7% of familial PD cases while observed in only 0.004% of controls (odds ratio of 65.5). This variant was confirmed in all cases via Sanger sequencing and segregated with PD in three families. RAB32 encodes a small GTPase known to interact with LRRK2 (refs. 3,4). Functional analyses showed that RAB32 S71R increases LRRK2 kinase activity, as indicated by increased autophosphorylation of LRRK2 S1292. Here our results implicate mutant RAB32 in a key pathological mechanism in PD-LRRK2 kinase activity5-7-and thus provide novel insights into the mechanistic connections between RAB family biology, LRRK2 and PD risk.

Novel recruitment strategy to enrich for LRRK2 mutation carriers.

The LRRK2 G2019S mutation is found at higher frequency among Parkinson disease (PD) patients of Ashkenazi Jewish (AJ) ancestry. This study was designed to test whether an internet-based approach could be an effective approach to screen and identify mutation carriers. Individuals with and without PD of AJ ancestry were recruited and consented through an internet-based study website. An algorithm was applied to a series of screening questions to identify individuals at increased risk to carry the LRRK2 G2019S mutation. About 1000 individuals completed the initial screening. Around 741 qualified for mutation testing and 650 were tested. Seventy-two individuals carried at least one LRRK2 G2019S mutation; 38 with PD (12.5%) and 34 without (10.1%). Among the AJ PD participants, each affected first-degree relative increased the likelihood the individual was LRRK2+ [OR = 4.7; 95% confidence interval = (2.4-9.0)]. The same was not observed among the unaffected AJ subjects (P = 0.11). An internet-based approach successfully screened large numbers of individuals to identify those with risk factors increasing the likelihood that they carried a LRRK2 G2019S mutation. A similar approach could be implemented in other disorders to identify individuals for clinical trials, biomarker analyses and other types of research studies.

Copy number variation in familial Parkinson disease.

Copy number variants (CNVs) are known to cause Mendelian forms of Parkinson disease (PD), most notably in SNCA and PARK2. PARK2 has a recessive mode of inheritance; however, recent evidence demonstrates that a single CNV in PARK2 (but not a single missense mutation) may increase risk for PD. We recently performed a genome-wide association study for PD that excluded individuals known to have either a LRRK2 mutation or two PARK2 mutations. Data from the Illumina370Duo arrays were re-clustered using only white individuals with high quality intensity data, and CNV calls were made using two algorithms, PennCNV and QuantiSNP. After quality assessment, the final sample included 816 cases and 856 controls. Results varied between the two CNV calling algorithms for many regions, including the PARK2 locus (genome-wide p = 0.04 for PennCNV and p = 0.13 for QuantiSNP). However, there was consistent evidence with both algorithms for two novel genes, USP32 and DOCK5 (empirical, genome-wide p-values<0.001). PARK2 CNVs tended to be larger, and all instances that were molecularly tested were validated. In contrast, the CNVs in both novel loci were smaller and failed to replicate using real-time PCR, MLPA, and gel electrophoresis. The DOCK5 variation is more akin to a VNTR than a typical CNV and the association is likely caused by artifact due to DNA source. DNA for all the cases was derived from whole blood, while the DNA for all controls was derived from lymphoblast cell lines. The USP32 locus contains many SNPs with low minor allele frequency leading to a loss of heterozygosity that may have been spuriously interpreted by the CNV calling algorithms as support for a deletion. Thus, only the CNVs within the PARK2 locus could be molecularly validated and associated with PD susceptibility.

Genomewide linkage study of modifiers of LRRK2-related Parkinson's disease.

Mutations in the leucine-rich repeat kinase 2 gene, located at 12q12, are the most common known genetic causes of Parkinson's disease. Studies of leucine-rich repeat kinase 2 mutation carriers have shown incomplete and age-dependent penetrance, and previous studies have suggested that inherited susceptibility factors may modify the penetrance of leucine-rich repeat kinase 2 mutations. Genomewide linkage to age of onset of leucine-rich repeat kinase 2-related Parkinson's disease was evaluated in a sample of 113 leucine-rich repeat kinase 2 mutation carriers from 64 families using single-nucleotide polymorphism data from the Illumina HumanCNV370 genotyping array. Association between onset age and single-nucleotide polymorphisms under suggestive linkage peaks was also evaluated. The top logarithmic odds score for onset age (logarithmic odds score = 2.43) was in the chromosome 1q32.1 region. Moderate linkage to onset was also identified at 16q12.1 (logarithmic odds score = 1.58). Examination of single-nucleotide polymorphism association to Parkinson's disease onset under the linkage peaks revealed no statistically significant single-nucleotide polymorphism associations. The 2 novel genomic regions identified may harbor modifiers of leucine-rich repeat kinase 2-related Parkinson's disease onset age or penetrance, and further study of these regions may provide important insight into leucine-rich repeat kinase 2-related Parkinson's disease.

Maternal inheritance and mitochondrial DNA variants in familial Parkinson's disease.

BACKGROUND: Mitochondrial function is impaired in Parkinson's disease (PD) and may contribute to the pathogenesis of PD, but the causes of mitochondrial impairment in PD are unknown. Mitochondrial dysfunction is recapitulated in cell lines expressing mitochondrial DNA (mtDNA) from PD patients, implicating mtDNA variants or mutations, though the role of mtDNA variants or mutations in PD risk remains unclear. We investigated the potential contribution of mtDNA variants or mutations to the risk of PD. METHODS: We examined the possibility of a maternal inheritance bias as well as the association between mitochondrial haplogroups and maternal inheritance and disease risk in a case-control study of 168 multiplex PD families in which the proband and one parent were diagnosed with PD. 2-tailed Fisher Exact Tests and McNemar's tests were used to compare allele frequencies, and a t-test to compare ages of onset. RESULTS: The frequency of affected mothers of the proband with PD (83/167, 49.4%) was not significantly different from the frequency of affected females of the proband generation (115/259, 44.4%) (Odds Ratio 1.22; 95%CI 0.83-1.81). After correcting for multiple tests, there were no significant differences in the frequencies of mitochondrial haplogroups or of the 10398G complex I gene polymorphism in PD patients compared to controls, and no significant associations with age of onset of PD. Mitochondrial haplogroup and 10398G polymorphism frequencies were similar in probands having an affected father as compared to probands having an affected mother. CONCLUSIONS: These data fail to demonstrate a bias towards maternal inheritance in familial PD. Consistent with this, we find no association of common haplogroup-defining mtDNA variants or for the 10398G variant with the risk of PD. However, these data do not exclude a role for mtDNA variants in other populations, and it remains possible that other inherited mitochondrial DNA variants, or somatic mDNA mutations, contribute to the risk of familial PD.

Genomewide association study for onset age in Parkinson disease.

BACKGROUND: Age at onset in Parkinson disease (PD) is a highly heritable quantitative trait for which a significant genetic influence is supported by multiple segregation analyses. Because genes associated with onset age may represent invaluable therapeutic targets to delay the disease, we sought to identify such genetic modifiers using a genomewide association study in familial PD. There have been previous genomewide association studies (GWAS) to identify genes influencing PD susceptibility, but this is the first to identify genes contributing to the variation in onset age. METHODS: Initial analyses were performed using genotypes generated with the Illumina HumanCNV370Duo array in a sample of 857 unrelated, familial PD cases. Subsequently, a meta-analysis of imputed SNPs was performed combining the familial PD data with that from a previous GWAS of 440 idiopathic PD cases. The SNPs from the meta-analysis with the lowest p-values and consistency in the direction of effect for onset age were then genotyped in a replication sample of 747 idiopathic PD cases from the Parkinson Institute Biobank of Milan, Italy. RESULTS: Meta-analysis across the three studies detected consistent association (p < 1 x 10(-5)) with five SNPs, none of which reached genomewide significance. On chromosome 11, the SNP with the lowest p-value (rs10767971; p = 5.4 x 10(-7)) lies between the genes QSER1 and PRRG4. Near the PARK3 linkage region on chromosome 2p13, association was observed with a SNP (rs7577851; p = 8.7 x 10(-6)) which lies in an intron of the AAK1 gene. This gene is closely related to GAK, identified as a possible PD susceptibility gene in the GWAS of the familial PD cases. CONCLUSION: Taken together, these results suggest an influence of genes involved in endocytosis and lysosomal sorting in PD pathogenesis.

Alpha-synuclein and familial Parkinson's disease.

Whole gene duplications and triplications of alpha-synuclein (SNCA) can cause Parkinson's disease (PD), and variation in the promoter region (Rep1) and 3' region of SNCA has been reported to increase disease susceptibility. Within our cohort, one affected individual from each of 92 multiplex PD families showing the greatest evidence of linkage to the region around SNCA was screened for dosage alterations and sequence changes; no dosage or non-synonymous sequence changes were found. In addition, 737 individuals (from 450 multiplex PD families) that met strict diagnostic criteria for PD and did not harbor a known causative mutation, as well as 359 neurologically normal controls, were genotyped for the Rep1 polymorphism and four SNPs in the 3' region of SNCA. The four SNPs were in high LD (r(2) > 0.95) and were analyzed as a haplotype. The effects of the Rep1 genotype and the 3' haplotype were evaluated using regression models employing only one individual per family. Cases had a 3% higher frequency of the Rep1 263 bp allele compared with controls (OR = 1.54; empirical P-value = 0.02). There was an inverse linear relationship between the number of 263 bp alleles and age of onset (empirical P-value = 0.0004). The 3' haplotype was also associated with disease (OR = 1.29; empirical P-value = 0.01), but not age of onset (P = 0.40). These data suggest that dosage and sequence changes are a rare cause of PD, but variation in the promoter and 3' region of SNCA convey an increased risk for PD.

Genomewide association study for susceptibility genes contributing to familial Parkinson disease.

Five genes have been identified that contribute to Mendelian forms of Parkinson disease (PD); however, mutations have been found in fewer than 5% of patients, suggesting that additional genes contribute to disease risk. Unlike previous studies that focused primarily on sporadic PD, we have performed the first genomewide association study (GWAS) in familial PD. Genotyping was performed with the Illumina HumanCNV370Duo array in 857 familial PD cases and 867 controls. A logistic model was employed to test for association under additive and recessive modes of inheritance after adjusting for gender and age. No result met genomewide significance based on a conservative Bonferroni correction. The strongest association result was with SNPs in the GAK/DGKQ region on chromosome 4 (additive model: p = 3.4 x 10(-6); OR = 1.69). Consistent evidence of association was also observed to the chromosomal regions containing SNCA (additive model: p = 5.5 x 10(-5); OR = 1.35) and MAPT (recessive model: p = 2.0 x 10(-5); OR = 0.56). Both of these genes have been implicated previously in PD susceptibility; however, neither was identified in previous GWAS studies of PD. Meta-analysis was performed using data from a previous case-control GWAS, and yielded improved p values for several regions, including GAK/DGKQ (additive model: p = 2.5 x 10(-7)) and the MAPT region (recessive model: p = 9.8 x 10(-6); additive model: p = 4.8 x 10(-5)). These data suggest the identification of new susceptibility alleles for PD in the GAK/DGKQ region, and also provide further support for the role of SNCA and MAPT in PD susceptibility.

Clinical correlates of depressive symptoms in familial Parkinson's disease.

Depression is one of the most common nonmotor complications of Parkinson's disease (PD) and has a major impact on quality of life. Although several clinical factors have been associated with depression in PD, the relationship between depression and stage of illness as well as between depression and degree of disability remains controversial. We have collected clinical data on 1,378 PD cases from 632 families, using the Unified Parkinson's Disease Rating Scale (UPDRS) Parts II (activities of daily living) & III (motor), the Mini-Mental State Exam, the Geriatric Depression Scale (GDS), and the Blessed Functional Activity Scale (Blessed). Analyses were performed using the 840 individuals with verified PD and without evidence of cognitive decline. Logistic regression was used to identify study variables that individually and collectively best predicted the presence of depressive symptoms (GDS >or= 10). After correcting for multiple tests, depressive symptoms were significantly associated with Hoehn and Yahr stage and other clinical measures but not with any genetic variant (parkin, LRRK2, APOE). The Blessed score, education, presence of a first degree relative with signs of depression, and UPDRS Part II were found to best predict depressive symptomatology (R(2) = 0.33; P = 4 x 10(-48)). Contrary to several reports, the results from this large study indicate that stage of illness, motor impairment, and functional disability are strongly correlated with depressive symptoms.

R1514Q substitution in Lrrk2 is not a pathogenic Parkinson's disease mutation.

Mutations in LRRK2 were first reported as causing Parkinson's disease (PD) in late 2004. Since then, approximately a dozen LRRK2 substitutions have been identified that are believed to be pathogenic mutations. The substitution of adenine for guanine at nucleotide 4541 (4541G>A) in LRRK2 was recently reported. This substitution resulted in the replacement of an arginine at position 1514 with a glutamine (R1514Q). Although this substitution was not found in a large cohort of controls, its pathogenicity could not be verified. We have now genotyped the R1514Q substitution in a sample of 954 PD patients from 429 multiplex PD families. This substitution was identified in 1.8% of the PD patients; however, the majority of the PD sibships segregating this substitution were discordant for this putative mutation. In addition, the R1514Q substitution was detected in 1.4% of neurologically evaluated, control individuals. These data suggest that the R1514Q variant is not a pathogenic LRRK2 mutation. We believe it is imperative that the causative nature of any newly identified genetic variant be determined before it is included in any panel for diagnostic testing.

Mutations in LRRK2 other than G2019S are rare in a north American-based sample of familial Parkinson's disease.

A total of 956 individuals with Parkinson's disease (PD) from 430 multiplex PD pedigrees were screened for 12 previously reported, pathogenic LRRK2 mutations: R793M, L1114L, I1371V, R1441C, R1441G, R1441H, Y1699C, M1869T, I2012T, I2020T, G2385R, and IVS31 +3G > A. Previous screening identified the LRRK2 G2019S mutation in 5% of our families. Only 1 of the 12 newly screened mutations, R1441C, was detected in a single family in our patient cohort. These results indicate that, although the G2019S mutation remains the most common mutation identified in familial PD patients, other mutations in LRRK2 are infrequent.

Mutations in DJ-1 are rare in familial Parkinson disease.

Mutations in DJ-1 (PARK7) are one cause of early-onset autosomal-recessive parkinsonism. We screened for DJ-1 mutations in 93 affected individuals from the 64 multiplex Parkinson disease (PD) families in our sample that had the highest family-specific multipoint LOD scores at the DJ-1 locus. In addition to sequencing all coding exons for alterations, we used multiplex ligation-dependent probe amplification (MLPA) to examine the genomic copy number of DJ-1 exons. A known polymorphism (R98Q) was found in five PD subjects, once as a homozygote and in the other four cases as heterozygotes. No additional missense mutations and no exon deletions or duplications were detected. Our results, in combination with those of previous studies, suggest that alterations in DJ-1 are not a common cause of familial PD.

Presence of an APOE4 allele results in significantly earlier onset of Parkinson's disease and a higher risk with dementia.

The epsilon4 allele of the apolipoprotein E gene (APOE4) has been consistently associated with a greater risk of Alzheimer's disease (AD) as well as an earlier onset of AD. It is possible that APOE4 may also play a role in the etiology of other neurodegenerative disorders, such as Parkinson's disease (PD). APOE genotype, age of onset, disease duration, smoking history, and dementia status were collected for families with PD, yielding 324 Caucasian families with complete information. Logistic regression employing one individual per family and including age of onset and disease duration as covariates demonstrated a significantly increased risk of dementia for those individuals having inherited at least one epsilon4 allele (OR=3.37; P=0.002). Survival analyses also demonstrated a significantly earlier age of onset for those subjects with at least one epsilon4 allele (59.7 years) as compared with those homozygous for the more common epsilon3 allele (62.4 years; P=0.009). Thus, consistent with previous studies, we find evidence that the presence of an epsilon4 allele results in significantly earlier onset of PD and a greater likelihood of dementia. It appears the similarities between PD and AD may be due to an overlap in the diseases' genetic etiology.

Genetic screening for a single common LRRK2 mutation in familial Parkinson's disease.

Mutations in the leucine-rich repeat kinase 2 (LRRK2) gene cause some forms of autosomal dominant Parkinson's disease. We measured the frequency of a novel mutation (Gly2019 ser) in familial Parkinson's disease by screening genomic DNA of patients and controls. Of 767 affected individuals from 358 multiplex families, 35 (5%) individuals were either heterozygous (34) or homozygous (one) for the mutation, and had typical clinical findings of idiopathic Parkinson's disease. Thus, our results suggest that a single LRRK2 mutation causes Parkinson's disease in 5% of individuals with familial disease. Screening for this mutation should be a component of genetic testing for Parkinson's disease.

Evaluation of the role of Nurr1 in a large sample of familial Parkinson's disease.

Parkinson's disease (PD) is a common neurodegenerative disorder in humans with wide variability in the age of disease onset. Although the disease has been thought previously to occur sporadically in most patients, there is increasing evidence of a genetic contribution to the disorder. Recently, a polymorphic variant within intron 6 of the Nurr1 gene was reported to be associated with sporadic and familial PD. In an effort to identify susceptibility genes for PD, we have collected 783 PD patients from 372 families and 397 healthy controls from 217 families. PD patients and healthy controls were genotyped for the intron 6 insertion polymorphism by BseRI restriction endonuclease digestion. No significant difference in either homozygosity or heterozygosity for the 7048G7049 (IVS6 1361 +16insG) polymorphism was detected in the PD patient cohort as compared with the panel of healthy controls. Moreover, direct sequencing of exon 1 of the Nurr1 gene in PD patients failed to detect either of the two recently reported Nurr1 mutations identified in a small subset of a PD patient cohort. Taken together, these data suggest that genetic alteration at the Nurr1 locus is not a significant risk factor for the development of Parkinson's disease in our large sample of familial PD patients.

Genome-wide linkage analysis and evidence of gene-by-gene interactions in a sample of 362 multiplex Parkinson disease families.

Parkinson disease (PD) is the second most common neurodegenerative disorder. We studied 754 affected individuals, comprising 425 sibling pairs, to identify PD susceptibility genes. Screening of the parkin gene was performed in a subset of the sample having earlier age of PD onset or a positive LOD score with a marker in the parkin gene. All subjects were evaluated using a rigorous neurological assessment. Two diagnostic models were considered for genome-wide, non-parametric linkage analyses. Model I included only those individuals with a more stringent diagnosis of verified PD (216 sibling pairs) and resulted in a maximum LOD score of 3.4 on chromosome 2. Model II included all affected individuals (425 sibling pairs) and yielded a LOD score of 3.1 on the X chromosome. Our large sample was then employed to test for gene-by-gene (epistatic) interactions. A genome screen using the 23 families with PD patients having a mutation in only one allele of the parkin gene detected evidence of linkage to chromosome 10 (LOD=2.3). The 85 families with a very strong family history of PD were employed in a genome screen and, in addition to strong evidence of linkage to chromosome 2 (LOD=4.9), also produced a LOD of 2.4 on chromosome 14. A genome screen performed in the 277 families without a strong family history of PD detected linkage to chromosomes 10 (LOD=2.4) and X (LOD=3.2). These findings demonstrate consistent evidence of linkage to chromosomes 2 and X and also support the hypothesis that gene-by-gene interactions are important in PD susceptibility.

Reliability of reported age at onset for Parkinson's disease.

An individual's age at onset of Parkinson disease (PD) can be collected through a variety of sources, including medical records, family report, and clinical observation. The most common source of PD age at onset information in the research setting is family-report, which is then typically used to classify a subject as juvenile, young, or late age at onset. The reliability of the family-reported age at onset of PD has not been rigorously examined. The present study used data from individuals diagnosed with PD to evaluate the reliability of age at onset information by comparing data obtained from three sources: 1) the subject's medical records, 2) a Family History Questionnaire, and 3) a Subject History Questionnaire. Among the 149 subjects with data for all three age at onset sources, the estimated reliability was R = 0.94. Similar reliability was observed when the sample was stratified based on gender, age at examination, disease duration, first symptom of PD, and years of education. The three measures of age at onset of PD show excellent agreement, strengthening confidence in the reliability of the reported age of clinical onset for PD.

Significant linkage of Parkinson disease to chromosome 2q36-37.

Parkinson disease (PD) is the second most common neurodegenerative disorder, surpassed in frequency only by Alzheimer disease. Elsewhere we have reported linkage to chromosome 2q in a sample of sibling pairs with PD. We have now expanded our sample to include 150 families meeting our strictest diagnostic definition of verified PD. To further delineate the chromosome 2q linkage, we have performed analyses using only those pedigrees with the strongest family history of PD. Linkage analyses in this subset of 65 pedigrees generated a LOD score of 5.1, which was obtained using an autosomal dominant model of disease transmission. This result strongly suggests that variation in a gene on chromosome 2q36-37 contributes to PD susceptibility.

Genome screen to identify susceptibility genes for Parkinson disease in a sample without parkin mutations.

Parkinson disease (PD) is a common neurodegenerative disorder characterized by bradykinesia, resting tremor, muscular rigidity, and postural instability, as well as by a clinically significant response to treatment with levodopa. Mutations in the alpha-synuclein gene have been found to result in autosomal dominant PD, and mutations in the parkin gene produce autosomal recessive juvenile-onset PD. We have studied 203 sibling pairs with PD who were evaluated by a rigorous neurological assessment based on (a) inclusion criteria consisting of clinical features highly associated with autopsy-confirmed PD and (b) exclusion criteria highly associated with other, non-PD pathological diagnoses. Families with positive LOD scores for a marker in an intron of the parkin gene were prioritized for parkin-gene testing, and mutations in the parkin gene were identified in 22 families. To reduce genetic heterogeneity, these families were not included in subsequent genome-screen analysis. Thus, a total of 160 multiplex families without evidence of a parkin mutation were used in multipoint nonparametric linkage analysis to identify PD-susceptibility genes. Two models of PD affection status were considered: model I included only those individuals with a more stringent diagnosis of verified PD (96 sibling pairs from 90 families), whereas model II included all examined individuals as affected, regardless of their final diagnostic classification (170 sibling pairs from 160 families). Under model I, the highest LOD scores were observed on chromosome X (LOD score 2.1) and on chromosome 2 (LOD score 1.9). Analyses performed with all available sibling pairs (model II) found even greater evidence of linkage to chromosome X (LOD score 2.7) and to chromosome 2 (LOD score 2.5). Evidence of linkage was also found to chromosomes 4, 5, and 13 (LOD scores >1.5). Our findings are consistent with those of other linkage studies that have reported linkage to chromosomes 5 and X.