BioMed Research International
Volume 2015 (2015), Article ID 678701, 6 pages
http://dx.doi.org/10.1155/2015/678701

Research Article

Evaluating LRRK2 Genetic Variants with Unclear Pathogenicity

Fathima Shaffra Refai, Shin Hui Ng, and Eng-King Tan

Department of Neurology, National Neuroscience Institute, 7 Hospital Drive, Outram Road, Singapore 169611

Received 20 September 2014; Revised 19 January 2015; Accepted 26 January 2015

Academic Editor: Hiroyuki Tomiyama

Copyright © 2015 Fathima Shaffra Refai et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Mutations in the leucine-rich repeat kinase 2 (LRRK2) have been known to be a major genetic component affecting Parkinson’s disease (PD). However, the pathogenicity of many of the LRRK2 variants is unclear because they have been detected in single patients or also in patients and controls. Here, we selected 5 exonic variants (L1165P, T1410M, M1646T, L2063X, and Y2189C) from each of the protein domain of LRRK2 and analysed their possible association with pathogenicity using in vitro functional assays. Point mutations representing each of these variants were incorporated into the LRRK2 gene, and functional aspects such as the percentage of cell survival upon application of stress and kinase activity were measured. Our results showed that all 5 variants had a significantly negative effect on the survival of cells, in both presence and absence of stress, as compared to the wild-type. In addition, there was also a slight increase in kinase activity in most of the variants in comparison to the wild-type. A negative correlation between cell survival and kinase activity was observed. These data suggest that most of the variants despite being located in different domains of LRRK2 appear to exert a potential pathogenic effect possibly through an increased kinase activity, supporting a gain of function mechanism.

1. Introduction

Over the past decade, several genes have been linked to the pathogenesis of Parkinson’s disease (PD), a neurodegenerative disorder characterized by the loss of dopaminergic neurons in the substantia nigra pars compacta of the brain [1]. Among these genes, mutations in leucine-rich repeat kinase 2 (LRRK2; PARK8; MIM number 609007) are the most common cause of Mendelian inheritance of PD, contributing to at least 5% of familial PD patients [2]. Interestingly, mutations in LRRK2 have also been found in 1-2% of sporadic PD cases, narrowing the lines distinguishing between sporadic and genetic causes of the disease [3]. Polymorphic LRRK2 variants have also been shown to modulate risk of PD [4–7].

The LRRK2 gene is encoded by a segment containing 51 exons, located in the chromosome 12q12 locus, which is highly conserved among most vertebrates [8]. It encodes a large 2527-amino-acid protein belonging to the ROCO family [8, 9] and is composed of several distinct domains: leucine-rich repeat (LRR), Roc (Ras in complex proteins) [10], COR (C-terminal of ROC), tyrosine and serine/threonine kinase (MAPKKK), and WD40 domain. Even though the precise role of LRRK2 is yet to be understood, the presence of these domains indicates its participation in cell signaling pathways [11].

To date, up to 7649 different LRRK2 sequence variants have been identified and reported worldwide, from many different ethnicities (http://www.ncbi.nlm.nih.gov/projects/SNP/). However, only 6 of these variants (N1437H, R1441C, R1441G, Y1699C, G2019S, and I2020T) are considered to be definitely disease causing, based on their absence in controls and cosegregation of the disease in families [12]. The role of many other variants and their contribution to the disease is currently unknown. One of the key questions is whether these variants located in different protein domains exert their pathogenicity (if any) via an increased kinase activity. This has important implications in maximizing therapeutic approaches in LRRK2-linked PD.

To address these gaps in knowledge, we selected 5 LRRK2 variants, one from each domain, whose pathogenicity is unclear because they are present in single patients or in patients and controls. These could be rare polymorphic variants, benign or genetic mutations. We designed an in vitro cell-based system, expressing the 5 variants, and analysed their toxicity using kinase and cell survival assays.

2. Materials and Methods

2.1. LRRK2 Expression Constructs

The pEGFP-N1-LRRK2-WT plasmid coding for the human LRRK2 cDNA (GenBank: BC117180.1) sequence was used as a template to generate point mutations, 3494T>C (L1165P), 4229C>T (T1410M), 4937T>C (M1646T), 6187_6191delCTCTA (L2063X), and 6566A>G (Y2189C) (Figure 1), using the QuickChange site-directed mutagenesis kit (Stratagene) in accordance with the manufacturer’s protocol. The authenticity of the resulting plasmids was confirmed via sequencing using a sequencer (ABI PRISM 3100 Genetic Analyzer; Applied Biosystems).

Figure 1: Schematic representation of LRRK2 protein, showing variants reported in patients with PD in relation to the domains.

2.2. Cell Culture and Transfection

Human embryonic kidney 293T (HEK 293T) cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM; Invitrogen) supplemented with 10% fetal bovine serum, 100 U/mL penicillin, and 100 mg/mL streptomycin, 1X nonessential amino acids, and 100 mg/mL sodium pyruvate (all from Gibco-BRL) at 37°C in a 5% CO₂ incubator. The plasmids were transfected using Lipofectamine 2000 (Invitrogen) with reference to the manufacturer’s protocol.

2.3. Treatment with H₂O₂

HEK 293T cells transiently expressing wild-type and mutant forms of LRRK2 were incubated in culture medium containing 500 μM H₂O₂ for 24 hours. The optimal concentration for the treatment of HEK 293T cells was determined by investigating the H₂O₂ concentration that would elicit a cell death of 50% (EC₅₀) from a dose-response curve.

2.4. Analysis of Cell Survival

Upon treatment with H₂O₂, cell viability was assessed using Cell Titer 96 AQueous One Solution Reagent (MTS Assay; Promega). Cells were incubated with MTS/DMEM solution (1 : 5 ratio) at 37°C for 1.5 hours in the dark. The resulting solution was mixed well and the optical density (OD) measured at 490 nm using a spectrophotometer (Benchmark plus, Microplate spectrophotometer; Biorad). Results were expressed as a percentage of untreated wild-type control.

2.5. Assessment of Kinase Activity

In vitro LRRK2 kinase activity was measured using Transcreener ADP² FI Assay (BellBrook Labs). LRRK2 variants were purified using agarose anti-GFP beads (Vector Laboratories). The purified proteins were incubated in 25 μL of kinase buffer (15 mM HEPES, 20 mM NaCl, 1 mM EGTA, 0.02% Tween 20, 10 mM MgCl₂, and 0.1 mg/mL BGG) containing 100 μM ATP for 1 hour at 30°C in a black opaque 96-well plate. The reactants were then mixed in a ratio of 1 : 1 with ADP Detection Mixture containing 1X Stop & Detect buffer B, 8 nM ADP Alexa594 Tracer, and 94 μg/mL ADP² Antibody-IRDye QC-1, and incubated at room temperature for 30 minutes. The intensity of fluorescence emission was measured using the Synergy H1 Hybrid Microplate reader (Biotek, USA) at a wavelength of excitation of 590 nm and emission of 617 nm.

2.6. Statistical Analysis

The statistical significance was assessed using the two-tailed Student’s t-test, with a 1% level of significance to test the hypothesis. Correlation coefficients between kinase and cell survival data were used to analyse the significance of their relatedness.

3. Results

HEK 293T cells transiently expressing wild-type or mutant forms of LRRK2 were created. These include those expressing L1165P, T1410M, M1646T, L2063X, Y2189C, or wild-type LRRK2 (Figure 1).

3.1. Baseline Cell Survival

Cell lines expressing all mutants exhibited a significantly lower percentage of cell survival as compared to the wild-type LRRK2, in the absence of any form of stress () (Figure 2(a)). Mutant Y2189C from the WD40 domain induced the most amount of toxicity exhibiting a mean cell survival of only 66%. On average, cell survivability of all mutants ranged between 66 and 78%.

(a)

(b)

(c)

Figure 2: Impact of LRRK2 wild-type and mutants on basal cell survival and survival in response to H₂O₂ insult in HEK 293T cells and their corresponding kinase activities. (a) Changes in percentage cell survival between cells expressing LRRK2 wild-type and mutants from LRR (L1165P), ROC (T1410M), COR (M1646T), kinase (L2063X), and WD40 (Y2189C) domains, shown in response to 500 μM H₂O₂ for 24 hours. (b) Autokinase activity of the wild-type and mutants in RLU measured using the Transcreener ADP² FI Assay. (c) Correlation between autokinase activity and percentage cell survival. Data are means ± standard error of the mean (SEM), three readings per data point, from at least 3 independent experiments. Two-tailed Student’s t-test; # and versus percentage change in cell survival in cell expressing LRRK2 mutants with respect to the wild-type, in the presence and absence of H₂O₂, respectively; versus kinase activity with respect to wild-type LRRK2.

3.2. Cell Survival upon H₂O₂ Application

Similarly, upon application of stress induced by H₂O₂, cell lines expressing the mutant forms of LRRK2 exhibited a significantly lower amount of cell survival than that of the wild-type () (Figure 2(a)). The highest toxicity was induced by the COR domain mutant M1646T (cell survival = 37%), while the least toxicity was induced by the WD40 domain mutant Y2189C (cell survival = 46%).

3.3. Kinase Activity

Experiments depicting the autokinase activity of LRRK2 variants showed a significantly increased activity in all variants as compared to wild-type (), except for L2063X () (Figure 2(b)). L2063X is a truncation mutation in the kinase domain, and these results suggest that amino acids 1879 to 2063 are enough to maintain baseline kinase activity.

3.4. Correlation between Kinase Activity and Cell Survival

Correlation coefficient calculations showed a strong negative correlation between autokinase activity and percentage of cell survival in the absence ( = −0.81) and presence of stress ( = −0.71) (Figure 2(c)).

4. Discussion

Extensive research has been carried out to identify the numerous LRRK2 variants in PD. Although commonly occurring mutations, for example, the G2019S, highly frequent (30–40%) among the Ashkenazi Jews and the North African Arabs [20], have been well-studied, there is limited information available on most of the other variants [1, 21]. One of the key questions is whether these variants located in different protein domains exert their pathogenicity (if any) via an increased kinase activity. Development of specific LRRK2 kinase inhibitors has been the current key approach worldwide.

Recently, Ross et al. [15] in a case-control study assessed 121 exonic variants in 15,540 individuals (8611 patients and 6929 controls) from the white, Asian, and Arab-Berber populations. They identified risk associations of variant M1646T in white, A419V and G2385R in Asian, and Y2189C in Arab-Berber populations. In our study, we tested for the risk of carrying two of these variants (M1646T and Y2189C) together with some others (L1165P, M1646T, and L2063X) using functional assays. The selection criteria were based on the fact that these variants have been detected in single patients or also in patients and controls, and there have been no functional studies previously reported on these variants; hence, their pathogenicity is unclear. We selected one variant from each representative protein domain (Figure 1).

Our results from the cell survival and kinase studies suggest that all of these variants appear to be toxic, and carrying any of them would likely increase the risk of developing PD (Figure 2). The strong negative correlation between kinase activity and cell survival is consistent with previous findings showing that increased autokinase activity of LRRK2 G2019S is harmful to the cell [22, 23]. Smith et al. [22] reported a 3-fold increase in LRRK2 G2019S autokinase activity in comparison to wild-type, contributing to a 2-fold decrease in cell viability, while, in our hands, we previously reported a 2-fold increase in LRRK2 G2019S autokinase activity as compared to wild-type, which was accompanied by a slight increase in neuronal cell toxicity [24]. In another instance, Chan et al. [25] also noted that neuronal cells transiently expressing LRRK2 G2019S are 1.1 times more toxic as compared with wild-type. It is also possible that LRRK2 mutants induce both kinase-dependent and independent forms of cell death and this should be addressed in further studies.

Nonkinase domains may have a regulatory effect on the kinase domain, either directly or indirectly via the interaction with other proteins. The ROC domain comprising of conserved motifs for GTPase activity is known to regulate kinase activity by acting as a molecular switch, alternating between GDP- and GTP-bound states [22, 26]. When bound to GTP, the switch region located outside the domain is in an active state, leading to an increase in kinase activity, while when bound to GDP, the tertiary structure of the switch region is in an inactive state, thereby leading to a decrease in kinase activity. These conformational changes in the ROC domain are conveyed to the kinase domain via the COR domain which acts as molecular hinge. LRR and WD40 domains are composed of highly conserved folds often found in signalling proteins in which they play a role in protein-protein interaction [27]. Variants in these domains may influence the LRRK2 kinase activity by mediating the interaction with its substrate, for example, by causing hyperphosphorylation of a neuroprotective protein. Jorgensen et al. [28] reported the involvement of WD40 domain in LRRK2 dimerization, known to affect its kinase activity and by that cell death.

Only one case of PD, with the L1165P mutation has been reported, where the patient developed the first signs of disease at the age of 47 (Table 1) [14]. The patient displayed symptoms typical of PD, together with dementia, and, pathologically, it was classified as a Lewy body disease [13, 14]. The absence of this substitution in all of the controls screened confirmed its pathogenicity [14]. L1165 is known to be highly conserved across many species, and therefore a substitution is predicted to cause a dramatic change structurally in the Leucine-rich region [14]. SIFT analyses predicted that the L1165P substitution is not tolerated.

Table 1: LRRK2 1 mutations and associated clinical features.

Protein domain
Nucleotide change
Amino acid change
Number of patients
Number of controls
Reported clinical features
Source of information

Leucine-rich repeat (LRR)
3494T>C
L1165P
1
0
Age of onset around 47 years. Slow progressive disease with typical PD symptoms that respond well to levodopa therapy
Chen-Plotkin et al., 2008 [13]; Covy et al., 2009 [14]

Ras of complex proteins (ROC)
4229C>T
T1410M
19
11
Middle-age to late onset of disease, displaying typical symptoms of PD and a good response to levodopa therapy
Ross et al., 2011 [15]; Lesage et al., 2009 [16]; Abdalla-Carvalho et al., 2010 [17]

C-terminal of ROC (COR)
4937T>C
M1646T
4
2
Typical PD
Jasinska-Myga et al., 2010 [18]; Ross et al., 2011 [15]

MAP kinase kinase kinase (MAPKKK)
6187_6191delCTCTA
L2063X
1
2
Typical PD
Ross et al., 2011 [15]

WD40 repeat (WD40)
6566A>G
Y2189C
14
8
Middle-age to late onset of disease, showing characteristic PD symptoms with a positive response to levodopa therapy
Abdalla-Carvalho et al., 2010 [17]; Nuytemans et al., 2008 [19]; Jasinska-Myga et al., 2010 [18]; Ross et al., 2011 [15]

GenBank: BC117180.1.

There is ambiguity in the trend reported for the T1410M variant. Some have classified it to be nonpathogenic due to the fact that it is present in both cases as well as in controls (Table 1) [16], while others say that the mutation affects a highly conserved region, located in the ROC domain, and therefore may distort the tertiary structure of the protein, hence disrupting its GTP hydrolysis function [17, 26]. T1410 has also been identified as an autophosphorylation site by mass spectrometry [29, 30]. Our results together with SIFT analyses depict that T1410M tends more towards the pathogenic direction.

Not much has been mentioned about the variants M1646T and L2063X in published literature. L2063X is a truncation mutation in the kinase domain, but the kinase activity is not affected. The exact reason as to why it causes apoptosis is unclear. The Y2189C variant was first identified by Nuytemans et al. [19] in the Belgian population and he reported that in silico conservation analysis showed that the residue Tyr in itself was not evolutionarily conserved, but the aromatic nature of the amino acid was, and may therefore be important in maintaining the structure and function of the protein. They provided evidence from SIFT analyses that Y2189C is deleterious for the function of LRRK2. Other groups also reported its trend toward association due to an increase in the minor allele frequency in patients as compared to controls [18]. Our observation of the highest cellular toxicity being exerted by this mutant in the absence of H₂O₂, while the lowest toxicity in the presence of H₂O₂ depicts that WD40 domain may be less susceptible to oxidative stress.

Our study has limitations. First, we have only examined 5 variants and can only provide proof of principle findings specific to these variants only. Second, we have used autophosphorylation of LRRK2 as a surrogate marker of kinase activity. LRRK2 is likely to have physiologic substrates or protein interactors in vivo that can modulate its activity.

In conclusion, as proof of principle experiments, we investigated 5 LRRK2 variants (L1165P, T1410M, M1646T, L2063X, and Y2189C) with unclear pathogenicity located in each of the major protein domain. Our results suggest that these variants are proapoptotic and this appeared to correlate with an increased kinase activity in most variants. Further studies in animal models will be useful to further characterise their pathogenicity and this will potentially help future genetic testing programmes. Our study also suggests that some of the variants, even if they are not in the kinase domain, may act via a common toxic gain of function likely through an increased kinase activity. Specific LRRK2 inhibitors could have potential therapeutic uses in carriers of these variants.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors acknowledge the National Medical Research Council for their Parkinson’s disease Translational & Clinical Research (TCR) Programme support, Singapore Millennium Foundation, and Duke-NUS Graduate Medical School.

References

E.-K. Tan and L. M. Skipper, “Pathogenic mutations in Parkinson disease,” Human Mutation, vol. 28, no. 7, pp. 641–653, 2007.

A. di Fonzo, C. F. Rohé, J. Ferreira et al., “A frequent LRRK2 gene mutation associated with autosomal dominant Parkinson's disease,” The Lancet, vol. 365, no. 9457, pp. 412–415, 2005.

W. P. Gilks, P. M. Abou-Sleiman, S. Gandhi et al., “A common LRRK2 mutation in idiopathic Parkinson's disease,” The Lancet, vol. 365, no. 9457, pp. 415–416, 2005.

T. Peeraully and E. K. Tan, “Genetic variants in sporadic parkinson's disease: east vs west,” Parkinsonism and Related Disorders, vol. 18, no. 1, pp. S63–S65, 2012.

J. N. Foo, L. C. Tan, H. Liany et al., “Analysis of non-synonymous-coding variants of Parkinson's disease-related pathogenic and susceptibility genes in East Asian populations,” Human Molecular Genetics, vol. 23, no. 14, pp. 3891–3897, 2014.

E.-K. Tan, R. Peng, Y.-Y. Teo et al., “Multiple LRRK2 variants modulate risk of Parkinson disease: a Chinese multicenter study,” Human Mutation, vol. 31, no. 5, pp. 561–568, 2010.

U. Kumari and E. K. Tan, “LRRK2 in Parkinson's disease: genetic and clinical studies from patients,” FEBS Journal, vol. 276, no. 22, pp. 6455–6463, 2009.

B. Patra, A. J. Parsian, B. A. Racette, J. H. Zhao, J. S. Perlmutter, and A. Parsian, “LRRK2 gene G2019S mutation and SNPs [haplotypes] in subtypes of Parkinson's disease,” Parkinsonism and Related Disorders, vol. 15, no. 3, pp. 175–180, 2009.

A. Zimprich, S. Biskup, P. Leitner et al., “Mutations in LRRK2 cause autosomal-dominant parkinsonism with pleomorphic pathology,” Neuron, vol. 44, no. 4, pp. 601–607, 2004.

L. Bosgraaf and P. J. M. van Haastert, “Roc, a Ras/GTPase domain in complex proteins,” Biochimica et Biophysica Acta, vol. 1643, no. 1–3, pp. 5–10, 2003.

R. M. Sancho, B. M. H. Law, and K. Harvey, “Mutations in the LRRK2 Roc-COR tandem domain link Parkinson's disease to Wnt signalling pathways,” Human Molecular Genetics, vol. 18, no. 20, pp. 3955–3968, 2009.

A. B. Singleton, M. J. Farrer, and V. Bonifati, “The genetics of Parkinson's disease: progress and therapeutic implications,” Movement Disorders, vol. 28, no. 1, pp. 14–23, 2013.

A. S. Chen-Plotkin, W. Yuan, C. Anderson et al., “Corticobasal syndrome and primary progressive aphasia as manifestations of LRRK2 gene mutations,” Neurology, vol. 70, no. 7, pp. 521–527, 2008.

J. P. Covy, W. Yuan, E. A. Waxman, H. I. Hurtig, V. M. van Deerlin, and B. I. Giasson, “Clinical and pathological characteristics of patients with leucine-rich repeat kinase-2 mutations,” Movement Disorders, vol. 24, no. 1, pp. 32–39, 2009.

O. A. Ross, A. I. Soto-Ortolaza, M. G. Heckman et al., “Association of LRRK2 exonic variants with susceptibility to Parkinson's disease: a case-control study,” The Lancet Neurology, vol. 10, no. 10, pp. 898–908, 2011.

S. Lesage, C. Condroyer, A. Lannuzel et al., “Molecular analyses of the LRRK2 gene in Europeanand North African autosomal dominant Parkinson's disease,” Journal of Medical Genetics, vol. 46, no. 7, pp. 458–464, 2009.

C. B. Abdalla-Carvalho, C. B. Santos-Rebouças, B. C. Guimarães et al., “Genetic analysis of LRRK2 functional domains in Brazilian patients with Parkinson's disease,” European Journal of Neurology, vol. 17, no. 12, pp. 1479–1481, 2010.

B. Jasinska-Myga, J. Kachergus, C. Vilariño-Güell et al., “Comprehensive sequencing of the LRRK2 gene in patients with familial Parkinson's disease from North Africa,” Movement Disorders, vol. 25, no. 13, pp. 2052–2058, 2010.

K. Nuytemans, R. Rademakers, J. Theuns et al., “Founder mutation p.R1441C in the leucine-rich repeat kinase 2 gene in Belgian Parkinson's disease patients,” European Journal of Human Genetics, vol. 16, no. 4, pp. 471–479, 2008.

L. J. Ozelius, G. Senthil, R. Saunders-Pullman et al., “LRRK2 G2019S as a cause of Parkinson's disease in Ashkenazi Jews,” The New England Journal of Medicine, vol. 354, no. 4, pp. 424–425, 2006.

D. G. Healy, M. Falchi, S. S. O'Sullivan et al., “Phenotype, genotype, and worldwide genetic penetrance of LRRK2-associated Parkinson's disease: a case-control study,” The Lancet Neurology, vol. 7, no. 7, pp. 583–590, 2008.

W. W. Smith, Z. Pei, H. Jiang, V. L. Dawson, T. M. Dawson, and C. A. Ross, “Kinase activity of mutant LRRK2 mediates neuronal toxicity,” Nature Neuroscience, vol. 9, no. 10, pp. 1231–1233, 2006.

E. Greggio, S. Jain, A. Kingsbury et al., “Kinase activity is required for the toxic effects of mutant LRRK2/dardarin,” Neurobiology of Disease, vol. 23, no. 2, pp. 329–341, 2006.

D. C. Angeles, B.-H. Gan, L. Onstead et al., “Mutations in LRRK2 increase phosphorylation of peroxiredoxin 3 exacerbating oxidative stress-induced neuronal death,” Human Mutation, vol. 32, no. 12, pp. 1390–1397, 2011.

S. L. Chan, L.-L. Chua, D. C. Angeles, and E.-K. Tan, “MAP1B rescues LRRK2 mutant-mediated cytotoxicity,” Molecular Brain, vol. 7, no. 1, article 29, 2014.

J. Deng, P. A. Lewis, E. Greggio, E. Sluch, A. Beilina, and M. R. Cookson, “Structure of the ROC domain from the Parkinson's disease-associated leucine-rich repeat kinase 2 reveals a dimeric GTPase,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 5, pp. 1499–1504, 2008.

B. K. Gilsbach and A. Kortholt, “Structural biology of the LRRK2 GTPase and kinase domains: implications for regulation,” Frontiers in Molecular Neuroscience, vol. 7, article 32, 2014.

N. D. Jorgensen, Y. Peng, C. C.-Y. Ho et al., “The WD40 domain is required for LRRK2 neurotoxicity,” PLoS ONE, vol. 4, no. 12, Article ID e8463, 2009.

P. P. Pungaliya, Y. Bai, K. Lipinski et al., “Identification and characterization of a leucine-rich repeat kinase 2 (LRRK2) consensus phosphorylation motif,” PLoS ONE, vol. 5, no. 10, Article ID e13672, 2010.

S. Kamikawaji, G. Ito, and T. Iwatsubo, “Identification of the autophosphorylation sites of LRRK2,” Biochemistry, vol. 48, no. 46, pp. 10963–10975, 2009.