Human Molecular Genetics Advance Access originally published online on November 3, 2005
Human Molecular Genetics 2005 14(23):3549-3556; doi:10.1093/hmg/ddi376
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Comprehensive evaluation of common genetic variation within LRRK2 reveals evidence for association with sporadic Parkinson's disease
1Population Genetics and 2Information and Mathematical Sciences, Genome Institute of Singapore, Singapore 138672, 3Department of Neurology and 4Division of Research, SingHealth, Singapore General Hospital, Singapore 169611 and 5National Neuroscience Institute, Singapore, 308433
* To whom correspondence should be addressed at: Population Genetics, Genome, Genome Institute of Singapore, #02-01 60 Biopolis Street, Singapore 138672, Singapore. Tel: +65 64788088; Fax: +65 64789059; Email: liuj3{at}gis.a-star.edu.sg
Received August 15, 2005; Accepted September 30, 2005
GenBank accession no. AY792511
| ABSTRACT |
|---|
|
|
|---|
Parkinson's disease (PD) is a complex neurodegenerative disorder whose aetiologies are largely unknown. To date, mutations in six genes have been found causal for some rare familial forms of the disease and common variation within at least three of these is associated with the more common sporadic forms of PD. LRRK2 is the most recently identified familial PD gene, although its role in sporadic disease is unknown. In this study, we have performed the first comprehensive evaluation of common genetic variation within LRRK2 and investigated its contribution to risk of sporadic PD. We first characterized the linkage disequilibrium within LRRK2 using a panel of densely spaced SNPs across the gene. We then identified a subset of tagging-SNPs (tSNP) that capture the majority of common variation within LRRK2. Both single tSNP and tSNP haplotype analyses, using a large epidemiologically matched sporadic casecontrol series comprising 932 individuals, yielded significant evidence for disease association. We identified a haplotype that dramatically increases disease risk when present in two copies (OR=5.5, 95%CI=2.114.0, P=0.0001). Thus, we provide the first evidence that common genetic variation within LRRK2 contributes to the risk of sporadic PD in the Chinese population.
| INTRODUCTION |
|---|
|
|
|---|
Parkinson's disease (PD) is a complex neurodegenerative disorder currently thought to affect almost one million in the USA and as many as 1.7 million in China (1
The first familial PD gene identified was alpha-synuclein (SNCA). Mis-sense mutations are responsible for a small number of autosomal dominant PD cases (5
7
). A central role for SNCA in PD stems from the finding that fibrillar forms of the protein are a major structural component of LBs in both familial and sporadic cases (8
). Genomic multiplications leading to increased expression of SNCA found in other autosomal dominant PD families [PARK4 (MIM 605543
[OMIM]
)] (9
,10
) suggest that protein levels may be important in pathogenesis. Polymorphic variability within the SNCA promoter is also associated with increased risk for sporadic PD (11
).
Mutations in the parkin gene (PRKN) are responsible for about 50% of autosomal recessive (AR) PD cases [PARK2 (MIM 600166
[OMIM]
)] (12
). The genetic aberrations in this gene are numerous ranging from single base-pair substitutions and deletions to multiple exonic deletions and duplications. Although most parkin-linked disease is inherited in AR fashion, mutations have been described in sporadic cases and haploinsufficiency may also predispose to common PD (13
). A common promoter variant that decreases PRKN transcription is also associated with sporadic PD (14
). A function of the parkin protein is to act as an E3 ubiquitin protein ligase within the ubiquitin proteasomal pathway (UPP) (15
); one of the principle mechanisms by which damaged/misfolded proteins are cleared in the cell. The protein product of another gene initially implicated in familial PD, UCH-L1 (16
), also plays a role within the UPP. Impairment of UPP function is observed in sporadic cases (17
) and a common polymorphism that alters UCH-L1 enzymatic activity is inversely associated with risk of sporadic PD (18
).
The leucine-rich repeat kinase 2 gene (LRRK2) is the most recently identified causative gene for autosomal dominant PD (19
). Although its biological function is unknown at present, peptide sequence analogy suggests the LRRK2 protein is a member of the ROCO family (20
) that contains five conserved domains within the C-terminal half of the protein. Most LRRK2 mutations thus far identified are within exons encoding these putative functional domains and they account for up to 6% of autosomal dominant cases (21
). The most common of these, G2019S, may also account up to 2% of sporadic cases (22
). As variation in genes initially implicated in familial PD is becoming increasingly relevant to PD in the wider community, we sought to investigate whether common variation within LRRK2 may also predispose risk to sporadic PD. Here, we report our findings in the first study to address this question, using a large (n=932) epidemiologically matched sporadic casecontrol series.
| RESULTS |
|---|
|
|
|---|
LRRK2 linkage disequilibrium and tagging-SNPs
Linkage disequilibrium (LD) is defined as the non-random association of alleles at adjacent loci. Loci and thus their allelic variants are expected to segregate independently in a randomly mating population. However, if a particular allele at one locus is found together on the same chromosome with a specific allele at a second locus, more often than expected by chance, then the loci are in linkage disequilibrium (23
We characterized the LD and haplotype structure of LRRK2 (including 10 kb of sequence upstream from coding start, based on Genbank accession number AY792511) in 92 control individuals. Twenty-five SNPs (Table 1) with minor allele frequencies (MAF) ranging from 0.03 to 0.49 were genotyped, yielding an average density of 1 SNP/6 kb (Fig. 1B). Seventy-eight unique multilocus haplotypes were constructed and used to calculate pair-wise D' and r2 values (Fig. 1C and D).
|
|
Complete LD exists throughout much of the gene overall, suggestive of minimal recombination in this region. In addition, r2 values suggest that information for some SNPs is redundant. We therefore, sought to identify tSNPs that would efficiently capture common variation within LRRK2. As there is no unique definition of what constitutes a haplotype block and the number of tSNPs can be further reduced if inter-block disequilibrium is also considered (25
Association analyses
The tSNPset was genotyped in all cases (n=466) and remaining control individuals (n=374). We initially carried out single marker tests of association for each LRRK2 tSNP, by calculating odds ratios under recessive, additive and dominant models (Table 2). There was a marginally significant association between rs10506151 and disease under a dominant model (OR=1.3, 95%CI=1.01.7, P=0.03). However, carrying two copies greatly increased disease risk under both recessive (OR=3.7, 95%CI=1.87.8, P=0.0001) and additive models (OR=4.0, 95%CI=1.98.3, P=0.0001). This association remains highly significant after a conservative Bonferroni correction for multiple testing (P=0.002).
|
We next performed tSNP haplotype association analysis, as haplotype analysis may provide more power to detect association than single marker analyses alone (27
|
|
In addition to the identification of a risk haplotype, our results also raise the possibility that protective variants exist within LRRK2. Haplotype 3 (OR=0.8, 95%CI=0.60.9, P=0.04) and to a greater extent Haplotype 11 (OR=0.2, 95%CI=0.10.7, P=0.01) are both under-represented in the case group (Table 3).
In silico analysis of rs10506151
SNP rs10506151 is positioned within the intron following exon 16, 2203 bp downstream of exon 16 and 692 bp upstream of exon 17. In silico analysis shows that the wild-type allele (C) is conserved in vertebrates through to the rat. In addition, the SNP is within a putative 6-mer hnRNP K binding site (ACCCAA; SNP position shown in bold).
| DISCUSSION |
|---|
|
|
|---|
This study is the first to comprehensively evaluate common genetic variation within LRRK2 and assess its potential role in sporadic PD. By using a powerful tSNP approach, we have identified a common haplotype that is highly over-represented within cases (P=0.005) and, when present in two copies, significantly increases disease risk (OR=5.5, 95%CI=2.114.0, P=0.0001). Interestingly, this haplotype carries a SNP (rs10506151) that is within a putative hnRNP K binding site and whose wild-type allele is conserved in vertebrates. hnRNP K is a component of the heterogeneous ribonucleoprotein (hnRNP) particle that is implicated in a variety of cellular processes including transcription and translation (28
We acknowledge the possibility that our findings may represent spurious associationsa common pitfall of many association studies. However, we have taken care to minimize potential confounders, in particular by using a large casecontrol series well matched in terms of age, gender and ethnicity. This strategy is also believed sufficient to attenuate the possible, and debatable effect of population stratification (27
).
Paisan-Ruiz et al. (29
) recently reported negative association of four LRRK2 SNPs (two coding non-synonymous and two coding synonymous) in a smaller and unmatched Caucasian casecontrol series. However, patients in this study were selected on the basis of early onset and/or family history of PD and the level of LRRK2 variation assessed was extremely limited. It is also possible that sporadic disease risk associated with common variation in LRRK2 is influenced by other genetic and non-genetic factors, some of which may be population specific. It will therefore now be of value to investigate our findings in other populations.
In conclusion, our results strongly suggest that, in addition to the established role in familial disease, LRRK2 genetic variation may indeed make a significant contribution to common forms of PD in the general community. Isolation of the responsible factors will provide further clues to understanding the pathogenesis of this complex disorder and further study is warranted.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Samples
From the movement disorder clinics of two major tertiary institutions in Singapore (Singapore General Hospital and National Neuroscience Institute) 466 PD patients were included in this study. PD diagnosis was made in accordance with the UK PD Society Brain Bank Clinical Diagnostic Criteria of PD, by two movement disorders specialists (EK Tan and L Tan). The mean age of disease onset (AOO) was 60±12 years (range 2481), mean current age was 65±11 years (range 3091) and 56% were male. All patients were considered sporadic, as no family history of PD was evident. All were of Chinese ethnicity. None carry the common G2019S mutation (30
SNP selection and genotyping
Aiming for an average marker density of one SNP/5 kb, 30 SNPs with a MAF >0.02 in the Chinese population (according to HapMap Project population data, URL http://www.hapmap. org/) were selected for genotyping in 92 control samples. SNPs were genotyped using MALDI-TOF mass spectrometry using the Sequenom MassARRAYTM system (San Diego, CA, USA) as previously described (19
). Briefly, multiplex genotyping assays were designed using the Sequenom DESIGNER software. Initial PCR (5 ng of genomic DNA) and primer extension reactions were carried out according to the Sequenom MassEXTEND protocol. Primers and allele codes are shown in Table 1. After purification, 15 nl of primer extension product was analysed with a MassARRAY Sequenom-Bruker Spectrometer (Bruker Biosciences). SNPs with overall call rates <0.85 (percentage of samples giving a reliable genotype call) and/or those that significantly deviated from HardyWeinberg equilibrium (at
=0.05) were excluded from further analyses (n=5).
LD characterization, tSNP selection and performance evaluation
Pair-wise D' and r2 values were calculated using multilocus haplotypes inferred from phase-unknown genotype data using PHASE version 2.1 (31
). tSNPs were identified using the TSSA algorithm. This is a block-independent approach that finds a minimum set of tSNPs which uniquely distinguish a certain percentage of all possible haplotypes (25
).
was set at 0.95 in order that the tSNPset returned was able to uniquely distinguish 95% of all possible haplotypes.
Performance of the selected tSNPset to capture haplotypic variation within LRRK2 was assessed by calculating Rh2 between tSNP haplotypes and common (frequency >0.05) full haplotypes (comprising 25 SNPs) (32
). This value can range from 0 to 1, with 1 indicating that tSNP haplotypes can perfectly predict all common 25 SNP haplotypes. To evaluate tSNP performance in capturing unobserved genotypic variation within LRRK2 we performed the leave-one-out analysis (33
). In brief, each of the 25 genotyped SNPs (k) was dropped in turn and tSNPs were selected from the remaining SNPs (k1). Haplotype r2 (possible value range 01) was then calculated between each tSNP set and remaining SNPs. This evaluation can provide an unbiased and accurate estimate of tSNP performance in predicting genotypes that are not directly assessed (34
).
Statistical analysis
Odds ratios, 95%CIs and P-values were calculated using web-based simple interactive statistical analysis (SISA) tools (URL http://home.clara.net/sisa/index.htm). KruskalWallis tests were carried out using Microsoft® Excel 2002.
In silico analysis
Sequence conservation around rs10506151 was analysed using conservation tracks available at the University of California at Santa Cruz Human Genome Browser Gateway (May 2004 assembly) (URL: http://genome.ucsc.edu/cgi-bin/hg). Intronic regulatory sequence analysis was performed using the web-based Alternative Splicing Workbench (URL: http://www. ebi.ac.uk/asd-srv/wb.cgi).
| ACKNOWLEDGEMENTS |
|---|
The authors gratefully acknowledge all participants in this study. They also would like to thank Jason Ong and Seok Shin Ng for genotyping and Marie Wong for construction of Figure 1. This study was supported by funding from the Agency for Science and Technology and Research of Singapore.
Conflict of Interest statement: The authors declare no conflicts of interest.
| REFERENCES |
|---|
|
|
|---|
- Rocca, W.A. (2005) Prevalence of Parkinson's disease in China. Lancet Neurol., 4, 328329.[CrossRef][ISI][Medline]
-
Gelb, D.J., Oliver, E. and Gilman, S. (1999) Diagnostic criteria for Parkinson disease. Arch. Neurol., 56, 3339.
[Abstract/Free Full Text] - Forno, L.S. (1996) Neuropathology of Parkinson's disease. J. Neuropathol. Exp. Neurol., 55, 259272.[ISI][Medline]
- Vila, M. and Przedborski, S. (2004) Genetic clues to the pathogenesis of Parkinson's disease. Nat. Med., 10 (suppl.), S58S62.
-
Polymeropoulos, M.H., Lavedan, C., Leroy, E., Ide, S.E., Dehejia, A., Dutra, A., Pike, B., Root, H., Rubenstein, J., Boyer, R. et al. (1997) Mutation in the alpha-synuclein gene identified in families with Parkinson's disease. Science, 276, 20452047.
[Abstract/Free Full Text] - Kruger, R., Kuhn, W., Muller, T., Woitalla, D., Graeber, M., Kosel, S., Przuntek, H., Epplen, J.T., Schols, L. and Riess, O. (1998) Ala30Pro mutation in the gene encoding alpha-synuclein in Parkinson's disease. Nat. Genet., 18, 106108.[CrossRef][ISI][Medline]
- Zarranz, J.J., Alegre, J., Gomez-Esteban, J.C., Lezcano, E., Ros, R., Ampuero, I., Vidal, L., Hoenicka, J., Rodriguez, O., Atares, B. et al. (2004) The new mutation, E46K, of alpha-synuclein causes Parkinson and Lewy body dementia. Ann. Neurol., 55, 164173.[CrossRef][ISI][Medline]
-
Spillantini, M.G., Crowther, R.A., Jakes, R., Hasegawa, M. and Goedert, M. (1998) alpha-Synuclein in filamentous inclusions of Lewy bodies from Parkinson's disease and dementia with lewy bodies. Proc. Natl. Acad. Sci. USA, 95, 64696473.
[Abstract/Free Full Text] -
Singleton, A.B., Farrer, M., Johnson, J., Singleton, A., Hague, S., Kachergus, J., Hulihan, M., Peuralinna, T., Dutra, A., Nussbaum, R. et al. (2003) alpha-Synuclein locus triplication causes Parkinson's disease. Science, 302, 841.
[Free Full Text] - Chartier-Harlin, M.C., Kachergus, J., Roumier, C., Mouroux, V., Douay, X., Lincoln, S., Levecque, C., Larvor, L., Andrieux, J., Hulihan, M. et al. (2004) alpha-Synuclein locus duplication as a cause of familial Parkinson's disease. Lancet, 364, 11671169.[CrossRef][ISI][Medline]
- Pals, P., Lincoln, S., Manning, J., Heckman, M., Skipper, L., Hulihan, M., Van Den Broeck, M., De Pooter, T., Cras, P., Crook, J. et al. (2004) alpha-Synuclein promoter confers susceptibility to Parkinson's disease. Ann. Neurol., 56, 591595.[CrossRef][ISI][Medline]
-
Lucking, C.B., Durr, A., Bonifati, V., Vaughan, J., De Michele, G., Gasser, T., Harhangi, B.S., Meco, G., Denefle, P., Wood, N.W. et al. (2000) Association between early-onset Parkinson's disease and mutations in the parkin gene. French Parkinson's Disease Genetics Study Group. N. Engl. J. Med., 342, 15601567.
[Abstract/Free Full Text] -
Mata, I.F., Lockhart, P.J. and Farrer, M.J. (2004) Parkin genetics: one model for Parkinson's disease. Hum. Mol. Genet., 13 (Spec no. 1), R127R133.
[Abstract/Free Full Text] -
West, A.B., Maraganore, D., Crook, J., Lesnick, T., Lockhart, P.J., Wilkes, K.M., Kapatos, G., Hardy, J.A. and Farrer, M.J. (2002) Functional association of the parkin gene promoter with idiopathic Parkinson's disease. Hum. Mol. Genet., 11, 27872792.
[Abstract/Free Full Text] - Shimura, H., Hattori, N., Kubo, S., Mizuno, Y., Asakawa, S., Minoshima, S., Shimizu, N., Iwai, K., Chiba, T., Tanaka, K. et al. (2000) Familial Parkinson disease gene product, parkin, is a ubiquitin-protein ligase. Nat. Genet., 25, 302305.[CrossRef][ISI][Medline]
- Leroy, E., Boyer, R., Auburger, G., Leube, B., Ulm, G., Mezey, E., Harta, G., Brownstein, M.J., Jonnalagada, S., Chernova, T. et al. (1998) The ubiquitin pathway in Parkinson's disease. Nature, 395, 451452.[CrossRef][Medline]
- Mcnaught, K.S., Belizaire, R., Isacson, O., Jenner, P. and Olanow, C.W. (2003) Altered proteasomal function in sporadic Parkinson's disease. Exp. Neurol., 179, 3846.[CrossRef][ISI][Medline]
- Liu, Y., Fallon, L., Lashuel, H.A., Liu, Z. and Lansbury, P.T., Jr. (2002) The UCH-L1 gene encodes two opposing enzymatic activities that affect alpha-synuclein degradation and Parkinson's disease susceptibility. Cell, 111, 209218.[CrossRef][ISI][Medline]
- Zimprich, A., Biskup, S., Leitner, P., Lichtner, P., Farrer, M., Lincoln, S., Kachergus, J., Hulihan, M., Uitti, R.J., Calne, D.B. et al. (2004) Mutations in LRRK2 cause autosomal-dominant parkinsonism with pleomorphic pathology. Neuron, 44, 601607.[CrossRef][ISI][Medline]
- Bosgraaf, L. and Van Haastert, P.J. (2003) Roc, a Ras/GTPase domain in complex proteins. Biochim. Biophys. Acta, 1643, 510.[Medline]
- Pardo, L.M. and Van Duijn, C.M. (2005) In search of genes involved in neurodegenerative disorders. Mutat. Res., in press.
- Gilks, W.P., Abou-Sleiman, P.M., Gandhi, S., Jain, S., Singleton, A., Lees, A.J., Shaw, K., Bhatia, K.P., Bonifati, V., Quinn, N.P. et al. (2005) A common LRRK2 mutation in idiopathic Parkinson's disease. Lancet, 365, 415416.[ISI][Medline]
- Cardon, L.R. and Abecasis, G.R. (2003) Using haplotype blocks to map human complex trait loci. Trends Genet., 19, 135140.[CrossRef][ISI][Medline]
- Johnson, G.C., Esposito, L., Barratt, B.J., Smith, A.N., Heward, J., Di Genova, G., Ueda, H., Cordell, H.J., Eaves, I.A., Dudbridge, F. et al. (2001) Haplotype tagging for the identification of common disease genes. Nat. Genet., 29, 233237.[CrossRef][ISI][Medline]
- Liu, T., Sung, W.-K., Li, Y., Liu, J., Mittal, A. and Mao, P. (2005) Effective algorithms for Tag SNP selection. J. Bioinf. Comp. Biol., in press.
- Ahmadi, K.R., Weale, M.E., Xue, Z.Y., Soranzo, N., Yarnall, D.P., Briley, J.D., Maruyama, Y., Kobayashi, M., Wood, N.W., Spurr, N.K. et al. (2005) A single-nucleotide polymorphism tagging set for human drug metabolism and transport. Nat. Genet., 37, 8489.[ISI][Medline]
- Newton-Cheh, C. and Hirschhorn, J.N. (2005) Genetic association studies of complex traits: design and analysis issues. Mutat. Res., 573, 5469.[ISI][Medline]
- Bomsztyk, K., Denisenko, O. and Ostrowski, J. (2004) hnRNP K: one protein multiple processes. Bioessays, 26, 629638.[CrossRef][ISI][Medline]
- Paisan-Ruiz, C., Lang, A.E., Kawarai, T., Sato, C., Salehi-Rad, S., Fisman, G.K., Al-Khairallah, T., St George-Hyslop, P., Singleton, A. and Rogaeva, E. (2005) LRRK2 gene in Parkinson disease. Neurology, 44, 595600.
- Tan, E.K., Shen, H., Tan, L.C., Farrer, M., Yew, K., Chua, E., Jamora, R.D., Puvan, K., Puong, K.Y., Zhao, Y. et al. (2005) The G2019S LRRK2 mutation is uncommon in an Asian cohort of Parkinson's disease patients. Neurosci. Lett., 384, 327329.[CrossRef][ISI][Medline]
- Stephens, M., Smith, N.J. and Donnelly, P. (2001) A new statistical method for haplotype reconstruction from population data. Am. J. Hum. Genet., 68, 978989.[CrossRef][ISI][Medline]
- Stram, D.O., Haiman, C.A., Hirschhorn, J.N., Altshuler, D., Kolonel, L.N., Henderson, B.E. and Pike, M.C. (2003) Choosing haplotype-tagging SNPS based on unphased genotype data using a preliminary sample of unrelated subjects with an example from the Multiethnic Cohort Study. Hum. Hered., 55, 2736.[CrossRef][ISI][Medline]
- Weale, M.E., Depondt, C., Macdonald, S.J., Smith, A., Lai, P.S., Shorvon, S.D., Wood, N.W. and Goldstein, D.B. (2003) Selection and evaluation of tagging SNPs in the neuronal-sodium-channel gene SCN1A: implications for linkage-disequilibrium gene mapping. Am. J. Hum. Genet., 73, 551565.[CrossRef][ISI][Medline]
-
Iles, M.M. (2005) Obtaining unbiased estimates of tagging SNP performance. Ann. Hum. Genet., in press.
This article has been cited by other articles:
![]() |
L. Wang, C. Xie, E. Greggio, L. Parisiadou, H. Shim, L. Sun, J. Chandran, X. Lin, C. Lai, W.-J. Yang, et al. The Chaperone Activity of Heat Shock Protein 90 Is Critical for Maintaining the Stability of Leucine-Rich Repeat Kinase 2 J. Neurosci., March 26, 2008; 28(13): 3384 - 3391. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Marin The Parkinson Disease Gene LRRK2: Evolutionary and Structural Insights Mol. Biol. Evol., December 1, 2006; 23(12): 2423 - 2433. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. H. V. Schapira The Importance of LRRK2 Mutations in Parkinson Disease. Arch Neurol, September 1, 2006; 63(9): 1225 - 1228. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



