- Research article
- Open Open Peer Review
A systematic review and integrative approach to decode the common molecular link between levodopa response and Parkinson’s disease
BMC Medical Genomicsvolume 10, Article number: 56 (2017)
PD is a progressive neurodegenerative disorder commonly treated by levodopa. The findings from genetic studies on adverse effects (ADRs) and levodopa efficacy are mostly inconclusive. Here, we aim to identify predictive genetic biomarkers for levodopa response (LR) and determine common molecular link with disease susceptibility. A systematic review for LR was conducted for ADR, and drug efficacy, independently. All included articles were assessed for methodological quality on 14 parameters. GWAS of PD were also reviewed. Protein-protein interaction (PPI) analysis using STRING and functional enrichment using WebGestalt was performed to explore the common link between LR and PD.
From 37 candidate studies on levodopa toxicity, 18 genes were found associated, of which, CAn STR 13, 14 (DRD2) was most significantly associated with dyskinesia, followed by rs1801133 (MTHFR) with hyper-homocysteinemia, and rs474559 (HOMER1) with hallucination. Similarly, 8 studies on efficacy resulted in 4 genes in which rs28363170, rs3836790 (SLC6A3) and rs4680 (COMT), were significant. To establish the molecular connection between LR with PD, we identified 35 genes significantly associated with PD. With 19 proteins associated with LR and 35 with PD, two independent PPI networks were constructed. Among the 67 nodes (263 edges) in LR, and 62 nodes (190 edges) in PD pathophysiology, UBC, SNCA, FYN, SRC, CAMK2A, and SLC6A3 were identified as common potential candidates.
Our study revealed the genetically significant polymorphism concerning the ADRs and levodopa efficacy. The six common genes may be used as predictive markers for therapy optimization and as putative drug target candidates.
Parkinson’s disease (PD) is a second most common progressive neurodegenerative disorder followed by Alzheimer’s disease . It affects 1.5% of the global population over the age of 65 years . Characterised by motor symptoms, like gait dysfunctioning, bradykinesia, rigidity, and resting tremors, PD has been believed to be caused due to loss of dopamine at the dopaminergic neurons in the substantia nigra pars compacta . Along with the dopaminergic disruption, other non-motor dysfunctioning like depression, sleep disorder, dementia are also observed in PD patients which can be a plausible consequence of both dopaminergic and non-dopaminergic systems. Pathological confirmation is obtained by the presence of Lewy bodies- fibrillar aggregates, mostly consisting of protein alpha synuclein, in the affected neurons of the brain .
Levodopa (or L-Dopa), ever since its discovery, has been used as a potent anti-Parkinson’s medication and functions as symptoms alleviating therapy, by maintaining the dopamine concentration at the synapse and reduce the motor fluctuations observed in PD patients . Almost 15–20% of the patients do not respond to the therapy or show adverse profiles primarily, levodopa-induced dyskinesia  after 5 years of therapy. Managing ADR is thus one of the most challenging aspects of PD. Carriers of specific genetic polymorphisms of drug metabolising enzymes, drug transporters, drug receptors and proteins involved in drug pathway of anti-Parkinson’s drugs may predispose to adverse reactions or altered efficacy.
Several susceptibility loci have been studied already with the familial cases of PD, like SNCA (PARK1), LRRK2, PRKN (PARK2), PINK1 (PARK6), DJ-1 (PARK7) . However very less has been elucidated about the genetic background of the sporadic cases of PD. Neurodegenerative diseases including PD are multifactorial in nature. Mechanisms like mitochondrial dysfunction, Lewy body formation, oxidative stress, altered protein handling, and inflammatory change are considered to lead to cell dysfunction and death by apoptosis or autophagy. Ageing is one of the most studied risk factor for PD, and the biochemical changes that are a consequence of aging amplify these abnormalities in PD patients’ brain . Candidate studies have pin- pointed genes like NAT2, MAOB, GST, mitochondrial tRNA, S18Y variant of UCHL1, SNCA, MAPT H1 haplotype and LRRK2 . GWA studies have identified more risk loci: BST1, GAK, HLA-DR, ACMSD, STK39, MCCC1/LAMP3, SYT11, PARK16, FGF20, and GPNMB, but with lower significance to establish a valid association for clinical management . Also, since the mechanism of development and progression of PD have not been elucidated fully, current treatment options are only targeted at providing symptomatic respite. Understanding of these multiple aspects of PD may potentially reward this field of study for clinical intervention.
The aim of the present article is to summarize all the studies carried out on polymorphism-association of administration of levodopa on sporadic PD patients and its treatment outcome as ADR and the altered efficacy of the drug. We, also describe the interplay of the molecular pathways involved in the mechanism of levodopa induced ADRs, LR and the disease pathology. This is an attempt to identify the molecular targets as genes and if the polymorphisms in such genes predispose certain patient population susceptible to causing ADRs and altered efficacy. For this purpose, we perform a systematic review through several online databases, select the relevant articles on the basis of pre-defined inclusion and exclusion criteria based on the focus of our study, separately, for LR, and PD disease susceptibility. These articles are further assessed for their methodological quality and finally the data was extracted for the list of genes (and its variants) associated with the drug response and disease risk. This effort has been further elaborated using computational approaches like network modelling to rule out the systematic biases from high-throughput multiple datasets and identify if there is any molecular mechanism involved in LR and PD susceptibility that intersect each other. Such proteins can be plausible targets to minimize toxicity, elaborate the therapeutic efficacy and capture disease risk.
All the methodologies performed in the study were drawn following the Human Genome Epidemiology Network for the systematic review of genetic association studies [11,12,13,14] and the PRISMA guidelines .
Data source and search strategy
A systematic search in Medline  and Web of science  was performed using standard MeSH terms “Parkinsons’s disease”, “variant”, “Polymorphism”, “SNP”, “single nucleotide polymorphism”, “pharmacogenomics”, “response” with AND/OR Boolean operators to identify all the human studies on genetics of Parkinson’s disease and/or on drug response by anti-Parkinson drugs. Also, a check for the studies that were not identified by the previous search, of pharmacogenetic relevance from the PharmGKB database , were added using search term “Parkinson’s disease”. The searches were limited to human studies.
Study selection criteria
The study selection was carried out independently in two stages by two different authors (DG and MKM) from relevant articles published up to March 9, 2016. All the articles that are reviews, commentary, erratums, editorials, technical reports, news, evaluation studies were initially screened (n = 173). Articles that were in duplicate (n = 260) and published in languages other than English (n = 36) were also excluded. At first, the articles were screened by titles based on relevance, as obtained by the search. Secondly, the abstracts of all primarily screened articles were retrieved and assessed according to the inclusion and exclusion criteria provided in Additional file 1: s4. Further, the articles were distinctly segregated into ADR of levodopa and efficacy of the drug. Only full text articles were included in the final study corpus. In case of disagreement regarding the screening of the articles, an independent reviewer (PT) was consulted to resolve the discrepancies. The cross-references of the finally selected articles were also searched for additional relevant articles. Further, the MeSH term for the adverse effect with ‘levodopa’ or ‘L-Dopa’ were again retrieved to double-check for any missing articles of purpose. All the baseline univariate significant allelic/genotypic associations with ADR/s and with L-dopa efficacy are reported in the Additional file 2: Table S1(8a) and Additional file 3: Table S2(8b), respectively. Similar search was performed for the drug response related articles as well.
Data extraction and quality assessment
Data were extracted by DG and MKM and checked by PT and RK. In case of sequential or multiple publications from the same group of authors, only the recent article has been included or studies which report exclusive findings. Data extracted from each eligible publication is provided in Table 1 (complete in Additional file 2: Table S1(8a)) for ADR articles and Table 2 (complete in Additional file 3: Table S2(8b) for drug efficacy studies. Ethnicity was classified as African, Asian or Caucasian [19,20,21,22,23]. If the ethnicity was not reported, the source population based on the country in which the study was conducted was considered e.g. Chicago. The genetic associations were stratified by ethnicity/population to explore the inter-ethnic variations. All the populations of the subjects have been categorised into its respective super-populations based on the 1000genome project (Phase II).
On a systematic analysis of the GWAS, conducted so far on PD susceptibility risk, included twenty GWAS studies (consisting of 45,465 cases and 173,222 controls). Sixty-one loci have been identified as significantly associated with the disease risk (p ≤ 0.01 × 10−8) (Table 3). A detailed summary of these significant genetic variants obtained from GWAS in the field of PD has been represented in Additional file 4: Table S1. The genes BST1, CCDC62/HIP1R, DGKQ/GAK, GBA, ITGA8, LRRK2, MAPT, MCCC1/LAMP3, PARK16, SNCA, STK39, and SYT11/RAB25 are disease risk loci following the collaborative meta-analyses. SNCA (p ≤ 4.16 × 10−73) and MAPT (p ≤ 2.37 × 10−48) have been studied in 9 and 7 GWA studies establishing the functional relevance in the disease physiology, making them the most prominent loci. Followed by LRRK2 and GAK in 4 studies, GBA/ SYT11 and MCCC1/ LAMP3 in 3 studies respectively.
Two reviewers (DG and MKM) independently assessed the methodological quality of all the selected articles using a predefined set of criteria. All included studies were assessed for the quality of data presented by using modified criteria suggested by Wells K. et al. . The quality assessment was scored on 14 parameters (Additional file 1: s5 and s6), with a positive score awarded for each detail present in study, the lack of detail was described as either NA (not applicable) or NR (not reported). NA was assessed with an equal positive score and was given only when the study was deemed independent of the parameter; NR was equated to no scoring and was independently awarded by the authors if the methodology was found insufficient or unreported. The detailed list of the 14 parameters used for the quality assessment has been discussed in Additional file 1: s7. Conflicting scores were reached to a consensus upon discussing (RK and PT). If the score was obtained as 11 or higher, the study was ranked as high quality.
Protein-protein interaction network
To decipher the connecting molecular link between the roles of the genes studied for LR and disease risk would help us rule out the bias if any to identify the interacting proteins of drug response thus elucidate the genetic landscape of the disease. Two independent protein-protein interaction (PPI) networks were constructed using STRING application with minimum required interaction score of 0.7 (high confidence), active interaction sources were set for only known interactions (databases and experiments) and maximum number of interactors to show in the 1st shell- no more than 50, and no 2nd shell interaction . Pathway enrichment analysis was conducted by WebGestalt . Using pathway commons enrichment analysis, default GO slim classification, 0.001 significance level, and minimum no. of genes in a category was set at 5. The enrichment analysis was run adjusting the false discovery rate (FDR) using the Benjamini-Hochberg (BH) procedure to obtain the results independently for the two set of genes.
Search and study selection
The workflow of the search and study selection has been represented in Fig. 1. A total of 1041 articles were obtained, of which 469 studies were excluded that includes duplicates, reviews, and articles in other languages (n = 469). The remaining 572 articles were screened by title and abstract following the inclusion and exclusion criteria resulting in further exclusion of 498 articles: 189 studies discussing other diseases or comorbid conditions, 16 studies on familial PD were removed as current study focuses on sporadic form of the disease, 76 studies were based on animal or in vitro models for PD, 40 studies were not on genetic association and 166 others did not discuss any drug response, finally 11 papers on drugs other than levodopa were also excluded to narrow down the scope of current study to the most widely prescribed medication. A total of 74 eligible publications were further divided into studies that are on adverse effect of levodopa or the efficacy of the drug in the patient cohort. In case of unavailable full text articles, the authors were contacted (n = 22). Responses received (n = 16) were included in the study, rest excluded (n = 6). From the cross-references of the included studies, six additional articles fulfilled the inclusion criteria [19, 27,28,29,30,31]. Thus, finally 38 eligible publications on levodopa induced ADR and 8 on drug response had sufficient data available for extraction to carry forward the systematic review.
The methodological and demographic characteristics of the ADR studies and the drug efficacy related studies of levodopa are summarised in Table 1 and Table 2 respectively [complete in Additional file 2: Table S1(8a) and Additional file 3: Table S2(8b)]. A total of 4127 subjects with ADRs were enrolled in the 37 studies. Most of the studies included were cohort studies except four which were case-control studies [29, 32,33,34]. Apart from original research articles, four letters [19, 21, 35, 36], two brief reports [37, 38], and two short communications [33, 39] were included as they had sufficient data pertaining to our criteria of inclusion and exclusion, one randomized control trial  was included. All the patients recruited in the independent studies were primarily administered with levodopa alone (n = 21), dopamine agonist (n = 4), DDC inhibitor/ carbidopa (n = 4), COMT inhibitor/ entacapone/tolcapone (n = 5) or MAO-B inhibitor (n = 1). The dose of levodopa administered ranged between 200.00 to 805.14 mg/day. The range of follow up period for the recruited subjects in the studies was between 6 to 10.3 months. The PD subjects recruited were diagnosed by UK Brain Bank Criteria (UK BBC)  in 26 articles, by Gelb’s criteria  in 1 article, by CAPIT  (Core assessment programme for intra-cerebral transplantations) in 1 article or by an experienced neurologist in 2 articles. The different motor functioning assessment scales used are provided in Additional file 2: Table S1(8a). The variability in HY assessment scale of 3.5 or more have been used.
Out of thirty-seven studies, thirteen focussed on levodopa induced dyskinesia exclusively, three on other motor fluctuation exclusively, one on wearing on/off, three on hyper-homocysteinemia, five studies were on hallucination, one study each on COMT inhibitor induced toxicity and elevated liver transaminase levels, and twelve studies discussed multiple ADR in the same cohort of recruited patients. Motor fluctuations were observed to be the most common adverse effect of anti-Parkinson’s medications. Dyskinesia (subjects of tardive dyskinesia, peak dose dyskinesia, diaphasic dyskinesia were grouped together) being the most prevalent among the subjects with ADRs, was present in 45.72% of patients. This group had an early age at onset of motor symptoms, longer disease duration and 560.96 ± 321.97 mg/day levodopa daily mean doses. In addition to levodopa, around one-third of the total patients were administered with COMT or MAO-B inhibitor or DDC inhibitor like carbidopa, entacapone or tolcapone. In addition to dyskinesia, adverse effects like other motor complications (motor impulsivity, wearing on-off, chorea, dystonia) were observed in 35.74% of total ADR subjects, and hyper-homocysteinemia 2.62%. Hallucinations occur as a consequence of psychosis, hence the two have been synonymously used and their subjects have been summed up together constituting 12.38%. One paper each discussed about COMT inhibitor induced toxicity (0.26%) and elevated transaminase level (3.27%).
Europeans (EUR) constituted as the major population of the studied subjects comprising 47.47%, followed by East Asian (EAS) at 26.52%, South Asian (SAS) at 13.32% closely followed by American (AMR) 12.66%. Most of the studies recruited subjects of the same population except Cheshire P. et al. (2013)  included patients from UK and Australia, Foltynie T. et al. (2014)  included all UK Caucasians, one Afro-Carribean, two Asian-Indian and one half Caucasian and half Asian-Indian, Greenbaum L. et al. (2013)  included Jewish and Italian subjects and Ziegler DA. et al. (2014)  included 122 American and 1 Asian.
A total of 645 subjects have been include in the 8 studies of levodopa efficacy. Of these, three studies are brief communications [48,49,50]; three were cohort studies [51,52,53] and one each of RCT , and letter . Three studies recruited patients administered with levodopa alone [49, 53, 54], four studies prescribed levodopa with DDC inhibitors like Benserazide [50,51,52], Pyridoxine , and one on Dopamine agonist (Pramipexole)  and COMT inhibitor (carbidopa) . The mean levodopa dose administered is 356.8 mg/day. All the patients were diagnosed by UK BBC , expect Moreau C. et al.  by Gibb’s criteria , Tan EK et al.  and Xie T et al.  by neurologists and Devos D. et al.  do not mention the diagnosing criteria. Two hundred twenty seven subjects responded to the therapy assessed by the scoring criteria. The study populations included Asians (54.57%) and Europeans (45.42%).
The cumulative quality assessment score obtained by individual ADR studies are represented in Additional file 2: Table S1(8a) and that of drug efficacy in Additional file 3: Table S2(8b). In ADR studies, the mean methodological assessment score was calculated to be 10.56 (SD 2.15), range 7 to 13. On the modified scale, twenty of thirty seven articles were deemed as good quality with a cut off score of ≥11, thirteen articles scoring ≥9–10 were categorised under moderate quality and finally any scores below 9 were judged as poor quality which included four articles. In L-dopa efficacy studies, the mean methodological score was 11.13 (SD 1.86), ranging 7 to 13. Six articles qualified to be good, and one each as moderate and poor quality.
Genetic factors in levodopa induced adverse effects
From total publications, 40 variants within 18 genes (HOMER1, ADORA2A, ANKK1, MTHFR, DRD2, SLC6A3, COMT, UGT1A, ACE2, BDNF, ABCC8, RYR1, DRD3, GRIN2A, SLC6A4, HTR2A, CYP2D6, CCK) were found to have significant association (p ≤ 0.05) with any type of levodopa induced ADR in PD (Additional file 2: Table S1(8a)).
Europeans (EUR) studies included 4258 subjects and reported 30 variants associated with ADRs. Among EUR, 1347 subjects were Italian from 8 studies [32, 35, 36, 46, 57,58,59,60]. With 13, 14 CAn STR repeats in DRD2 gene being the most significant variant associated with dyskinesia, followed by rs1801133 (MTHFR) associated with hyper-homocysteinemia, rs474559 (HOMER1) with hallucinations. Other polymorphisms including rs886292 (ABCC8), rs11880894 (RYR1), rs1800497 (DRD2), rs6265 (BDNF), rs11646587, rs7192557 and rs8057394 (GRIN2A) were found to be associated with dyskinesia in patients administered with levodopa medication. 5-HTTLPR and rs6313 (5-HTR2A) were reported in higher frequency in patients with psychosis and Acuña G. et al. (2002)  suggested the association of 10 UGT1A SNPs with elevated liver transaminase level. Among East Asians (EAS), two Chinese studies [62, 63] constituting 746 subjects concluded, rs4680 (COMT) and I/D polymorphism (ACE2) associated in patients with wearing on-off and psychosis respectively. A Japanese study  showed significant association with -45C/T (CCK) in patients with hallucination and rs1801131, rs1801133 (MTHFR) with elevated plasma homocysteine levels in a Taiwanese study . In American (AMR) population, Brazilians constituted four studies [28, 64,65,66] determining rs4704559 (HOMER1), rs2298383 (ADORA2A), rs1800497 (ANKK1) associated with dyskinesia, rs3761422 (ADORA2A) with motor fluctuation and rs28363170 (DAT1) with hallucination. Among South Asian (SAS) population, Israelis were most abundantly studied representing association of rs393795 (DAT1), rs886292 (ABCC8), rs11880894 (RYR1) and rs1800497 (DRD2) with dyskinesia [39, 46, 67, 68].
Genetic factors in other LR
On elaborate systematic extraction of published literature on LR, in terms of efficacy, of the drug, only eight studies deemed our defined inclusion and exclusion criteria. The enzymes directly involved in the metabolism and activity of levodopa is evidently been mostly studied with the altered LR. rs4680 (COMT) , rs6280 (DRD3) , rs921451, rs3837091 (DDC) , rs28363170, rs3836790 (SLC6A3)  were the significant variants with reduced LR (Additional file 3: Table S2(8b)). However, no conclusive results could be drawn from this systematic analysis due to large variability and low significance.
Genetics of PD susceptibility
Employing GWAS dataset to stratify disease susceptibility loci, we follow an unbiased approach to identify such loci in sporadic PD cases. Nalls et al. (2014)  recently conducted a large scale meta-analysis to identify the associated loci with disease risk. Keeping this study as the base of the systematic review of all the GWAS on PD risk, and adding the recent studies to it. Twenty studies were included for the systematic review with 45,465 cases and 173,222 controls, mostly from including Caucasian population followed by Jewish, Chinese and Japanese. Sixty one loci in genes like BST1, CCDC62/HIP1R, TMEM175/DGKQ/GAK, GBA, ITGA8, LRRK2, MAPT, MCCC1/LAMP3, PARK16, SNCA, STK39, and SYT11/RAB25 were associated with the disease susceptibility as shown in Table 3. Additional file 4: Table S1 tabulates all the loci found to be associated with disease susceptibility and the significant single nucleotide polymorphism (SNPs) in bold (p value ≤1.0 × 10−8).
Protein-protein interaction network
We performed PPI analysis using genes obtained from the systematic review of LR and disease risk in order to understand the functional association among the genes in the respective gene modules. With 19 proteins associated with LR and 35 with PD, two independent PPI networks were constructed respectively using STRING database to identify critical candidate genes/proteins (Fig. 2). In LR, the 67 nodes, represent the genes, and 263 edges weight the likelihood of nodes in common biological functions.) In PD susceptibility, 62 nodes linked with 190 edges are depicted. Functional enrichment analysis using WebGestalt identified three common pathways (Alpha synuclein signalling, ADP-ribosylation factor 6 (Arf6) downstream pathway, and Insulin-like growth factor 1 (IGF1) pathway) in the top 20 pathways for LR and PD (Fig. 3). The common six proteins, UBC, SNCA, FYN, SRC, SLC6A3, CAMK2A, between LR and PD were thus obtained. The functional significance of the identified nodes in PPIs are clearly substantiated by the biological processes obtained from functional enrichment analysis.
PD is a progressive brain disease which causes significant movement disability . The treatment is aimed at symptomatic management rather than complete cure. However, the challenge is large clinical variability in drug response and adverse effects on prolonged therapy. Discerning the genetic factors responsible for this variability to the drug toxicity and efficacy can provide better clinical management. This study identifies such genetic variants in genes involved in L-dopa metabolism and in the disease etiology by a systematic review approach. Further from the limited number of genes obtained from the systematic review, we extended our effort to integrated computational approaches like network modelling and functional enrichment to identify the other interacting proteins and thereby distinguish the common proteins and molecular pathways that participate in LR and the disease. We additionally show the limitations of the published literature and give insights that may be useful to future studies.
We have implemented a modified scale of Wells K et al. (2009)  criteria with five additional parameters, to assess the quality of articles included in the systematic review. To the best of our knowledge our study has incorporated the most comprehensive methodological quality assessment scoring for screening articles of systematic review. Candidate gene studies have been screened for the systematic review of LR.
A total of 18 genes from the 37 ADR studies and 4 genes from the 8 efficacy studies were retrieved after the systematic review. Most of the genes are related to dopaminergic pathway and their role have been depicted in Fig. 4. Most of the studies included the genes related to dopaminergic pathway. For instance, among ADR studies, CAn STR 13, 14 (DRD2) was found to be most significantly associated with dyskinesia, rs1801133 (MTHFR) with hyper-homocysteinemia, and rs474559 (HOMER1) with hallucination. Carriers of 13, 14 alleles are found to have lower risk of developing dyskinesia but role of this repeat is still unknown. Patients with the TT677 (rs1801133, 677C > T) genotype exhibit 50% reduced activity of MTHFR enzyme, consequently elevating the plasma homocysteine levels . rs474559 G allele (HOMER1) have lower prevalence of dyskinesia as it might disrupt the glutamatergic transmission . In efficacy related studies, rs28363170, rs3836790 (SLC6A3) and rs4680 (COMT), were important. Individuals with rs3836790 6/6 or rs28363170 10/10 (SLC6A3) genotypes have higher transporter expression leading to lower dopamine levels at the synapse . The haplotype structure formed by four SNPs (rs6269: A > G, rs4633: C > T, rs4818: C > G, rs4680:A > G) characterises the COMT enzyme activity to low (ACCG), medium (ATCA) and high (GCGG) [72, 73]. Accordingly, the levodopa metabolism is affected, altering the synaptic dopamine concentration. Also we observed that SLC6A3, COMT, and DRD3 genes were common between the ADR and efficacy studies resulting in 19 exclusive genes from LR studies.
In addition, to identify the disease genes involved in drug response and vice versa, genes found implicated in PD susceptibility from GWAS were also retrieved. This led to obtaining 61 significantly associated SNPs (p value ≤1.0 × 10−8) pertaining to 35 genes (Additional file 4: Table S1). Then an integrated network analysis resulted in six common molecular targets (SNCA, FYN, SRC, UBC, CAMK2A, and SLC6A3) from the overlap between 67 nodes (and 263 edges) in LR and 62 nodes (and 190 edges) in PD pathophysiology, respectively. Among the six common molecular targets, SNCA has been widely established to be a major player in PD susceptibility as it a major component of Lewy bodies and mutant SNCA has a greater tendency to acquire misfolding [70, 71]. Aggregation of SNCA has been shown to be neurotoxic for the cell through the formation of intermediate aggregates called protofibrils . In a recent report the distinct role of alpha-synuclein forming fibrils as the major toxic, resulting in progressive motor impairment and cell death leading to neurotoxic phenotypes in PD is demonstrated . FYN, a tyrosine kinase family protein, found inside nerve cells and helps in communicating signals or chemical instructions between different cellular components. This protein has been observed to get modified on levodopa administration, causing dyskinesia . Wang et al. (2016) validated neuro-inflammation inhibition by SRC (SRC proto-oncogene, non-receptor tyrosine kinase) signalling pathway to be a potential drug and disease candidate which supports our finding SRC as the common molecular bridge to both drug response and disease pathology . UBC belongs to the Ubiquitin family C which carries out the ubiquitin mediated proteolysis and aggregate ubiquitin monomers in the diseased brain. The ubiquitin proteins cause aberrations in the ubiquitin proteasome system (UPS) leading to PD pathogenesis . The neuronal protein, CAMK2A alters with intracellular calcium ion concentration change that is abnormally activated following dopamine depletion thus modulating the neuronal function in striatum . Zhang et al. (2014) also established an interaction between CaMK2A and Dopamine D2 receptors in striatal neurons, sensitive to long-term levodopa administration to PD rats . Finally, the SLC6A3/DAT1 variants have a significant effect on striatal activation and performance in PD as suggested by Habak et al. These results have furnished evidences on the role of these candidates in both levodopa metabolism impairment and disease risk. Further, these plausible biomarkers might bridge the path between levodopa metabolism and disease pathology resulting in reduced ADRs, optimum efficacy and, accurate diagnosis.
Functional enrichment analysis revealed prominently , alpha-synuclein pathway to be the most significant candidate with the set of disease and response related genes respectively followed by other growth factor signalling pathways like Afr6 downstream pathway, IGF1 pathway, and so on. Arf6, ADP-ribosylation factor, signalling plays a role in the Ras- mediated cell signalling . It is also responsible in the intracellular trafficking of DRD2 by GRK and PKC proteins . The potential role of IGF1, Insulin-like growth factor 1, signalling has been studied with neurodegeneration in human, participating in functions like brain neuron survival, synaptic transmission as well as plasticity . Bernhard FP. et al. (2016)  also established that IGF1 might serve as a PD prediction marker, observing elevated levels of IGF1 in PD patients. Additional file 5: s2 and Additional file 6: s3 tabulate the enriched functions obtained by levodopa response genes and PD related genes. We highlighted the potential usefulness of these biological functions in PD treatment which can be affirmed by in vitro and/or in vivo model systems.
Although significant findings have been observed in our study, several limitations exist. The papers included in the systematic review presented high heterogeneity in terms of diagnosis, response criteria, drugs administered with different doses and genotyping techniques. As suggested by Schumacher-Schuh et al. (2014)  the phenotypic heterogeneity in terms of adverse effects lacks clinical instrument to adequately measure the ADR, whereas in terms of efficacy, several response rating scales have been incorporated. In GWAS studies, the assayed SNPs are usually to mark a genome region that influences the studied phenotype. However, we have picked up the annotated genes corresponding to the significantly associated SNPs from the respective studies, to identify the proteins that play a role in the biological processes which ultimately influences the phenotype. Motor fluctuation, a common ADR of levodopa, lacks clear clinical classification and hence assessment. A regular record of patient motor state could be preferred. Genetic heterogeneity is another source of variability between studies because different markers in the same genes were employed for these associations; moreover, patients with different genetic backgrounds may not be strictly comparable. One major limitation of network biology is the quality and the coverage of the interactions. The rate of discovery of false positives and false negatives are high which shows the need to rank the reported interactions for further validation.
In summary, the present study provides a framework for better understanding of the molecular interplay between L-Dopa metabolism with PD pathophysiology and also a means to evaluate putative biomarkers to bridge the gap in treatment outcome and disease risk. We propose the above six genes could be useful in predicting both the LR and disease risk, simultaneously. This however warrants further experimental validations to develop into a targeted therapy. Translating these evidences into future validation would present pre-diagnostic marker development which can be applicable in clinical manifestation. A definitive role of these molecular targets in the disease progression can also lead to substantive advancement in PD treatment.
- ABCC8 :
ATP binding cassette subfamily C member 8
- ACE2 :
Angiotensin I converting enzyme 2
- ACMSD :
Aminocarboxymuconate semialdehyde decarboxylase
- ADORA2A :
Adenosine A2a receptor
Adverse drug reactions
- ANKK1 :
Ankyrin repeat and kinase domain containing 1
- Arf6 :
ADP-ribosylation factor 6
- BDNF :
Brain derived neurotrophic factor
- BST1 :
Bone marrow stromal cell antigen 1
- CAMK2A :
Calcium/calmodulin-dependent protein kinase II alpha
Core assessment programme for intra-cerebral transplantations
- CCDC62/HIP1R :
Coiled-coil domain containing 62
- CCK :
- COMT :
- CYP2D6 :
Cytochrome P450 family 2 subfamily D member 6
- DDC :
- DJ-1 :
- DRD2 :
Dopamine Receptor D2
- DRD3 :
Dopamine receptor D3
False discovery rate
- FGF20 :
Fibroblast growth factor 20
- FYN :
FYN proto-oncogene, Src family tyrosine kinase
- GAK :
Cyclin G associated kinase
- GBA :
- GPNMB :
- GRIN2A :
Glutamate ionotropic receptor NMDA type subunit 2A
G protein-coupled receptor kinase
- GST :
Genome-wide association study
- HLA-DR :
Major histocompatibility complex, class II, DR
- HOMER1 :
Homer scaffolding protein 1
- HTR2A :
5-hydroxytryptamine receptor 2A
Hoehn and Yahr scale
- IGF1 :
Insulin-like Growth Factor 1
- ITGA8 :
Integrin subunit alpha 8
- LRRK2 :
Leucine rich repeat kinase 2
- MAOB :
Monoamine oxidase B
- MAPT :
Microtubule associated protein tau
- MCCC1/LAMP3 :
Methylcrotonoyl-CoA carboxylase 1
Medical subject headings
- MTHFR :
- NAT2 :
- PARK16 :
Parkinson disease 16 (susceptibility)
- PARK2 :
Parkinson disease (autosomal recessive, juvenile) 2, parkin
- PARK6 :
PTEN induced putative kinase 1
- PARK7 :
Parkinsonism associated deglycase
- PINK1 :
PTEN induced putative kinase 1
Protein kinase C
- PRKN :
Parkin RBR E3 ubiquitin protein ligase
- RYR1 :
Ryanodine receptor 1
- SLC6A3 :
Solute carrier family 6 member 3
- SLC6A4 :
Solute carrier family 6 member 4
- SNCA :
Single nucleotide polymorphisms
- SRC :
SRC proto-oncogene, non-receptor tyrosine kinase
- STK39 :
Serine/threonine kinase 39
Short tandem repeat
- SYT11 :
- TMEM175/DGKQ/GAK :
Transmembrane protein 175
- UBC :
- UCHL1 :
Ubiquitin C-terminal hydrolase L1
- UGT1A :
UDP glucuronosyltransferase family 1 member A complex locus
- UK BBC:
United Kingdom’s brain bank criteria
Gibrat CS-PM, Bousquet M, Lévesque D, Rouillard C, Cicchetti F. Differences between subacute and chronic MPTP mice models: investigation of dopaminergic neuronal degeneration and alpha-synuclein inclusions. J Neurochem. 2009;109:1469–82.
Blesa JLJ, Obeso JA. Parkinson's disease: cell vulnerability and disease progression. Front Neuroanat. 2015;15:125.
J J. Parkinson's disease: clinical features and diagnosis. J Neurol Neurosurg Psychiatry. 2008;79:368–76.
Schneider SAOJ. Clinical and pathological features of Parkinson's disease. Curr Top Behav Neurosci. 2015;22:205–20.
Group TPS. Levodopa and the progression of Parkinson's disease. N Engl J Med. 2004;351:2498–508.
Gilgun-Sherki Y, Djaldetti R, Melamed E, Offen D. Polymorphism in candidate genes: implications for the risk and treatment of idiopathic Parkinson's disease. Pharmacogenomics J. 2004;4:291–306.
M D, M B, M K. Pharmacogenetics of Parkinson’s disease – through mechanisms of drug actions. Curr Genomics. 2013;14:568–77.
Schapira AH, P J. Etiology and pathogenesis of Parkinson’s disease. Mov Disord. 2011;26:1049–55.
Tan EK. The role of common genetic risk variants in Parkinson disease. Clin Genet. 2007;72:387–93.
AB S, Farrer MJ, V B. The genetics of Parkinson's disease: progress and therapeutic implications. Mov Disord. 2013;28:14–23.
FK K, McQueen MB, Khoury MJ, Tanzi RE, Bertram L, JP I. Evaluation of the potential excess of statistically significant findings in published genetic association studies: application to Alzheimer's disease. Am J Epidemiol. 2008;168:855–65.
Ioannidis JP, Trikalinos TA, MJ K. Implications of small effect sizes of individual genetic variants on the design and interpretation of genetic association studies of complex diseases. Am J Epidemiol. 2006;164:609–14.
Little J, Higgins JP, Ioannidis JP, Moher D, Gagnon F, von Elm E, Khoury MJ, Cohen B, Davey-Smith G, Grimshaw J, et al. STrengthening the REporting of genetic association studies (STREGA)--an extension of the STROBE statement. Eur J Clin Investig. 2009;39:247–66.
Sagoo GS, Little J, JP H. Systematic reviews of genetic association studies. Human genome epidemiology network. PLoS Med. 2009;6:e28.
Hutton B, Salanti G, Caldwell DM, Chaimani A, Schmid CH, Cameron C, Ioannidis JP, Straus S, Thorlund K, Jansen JP, et al. The PRISMA extension statement for reporting of systematic reviews incorporating network meta-analyses of health care interventions: checklist and explanations. Ann Intern Med. 2015;162:777–84.
MEDLINE Database. https://www.ncbi.nlm.nih.gov/pubmed/. Accessed 9 Mar 2016.
Web of Science. http://apps.webofknowledge.com/WOS_GeneralSearch_input.do?%20product=WOS&search_mode=GeneralSearch&SID=X1cFbN52ilAsnVB8QjI&preferencesSaved=. Accessed 9 Mar 2016.
Hewett M OD, Rubin DL, Easton KL, Stuart JM, Altman RB, Klein TE: PharmGKB: the Pharmacogenetics Knowledge Base. Nucleic Acids Res 2002 Jan 2002, 30:163-165.
Kiferle L, Ceravolo R, Petrozzi L, Rossi C, Frosini D, Rocchi A, Siciliano G, Bonuccelli U, Murri L. Visual hallucinations in Parkinson's disease are not influenced by polymorphisms of serotonin 5-HT2A receptor and transporter genes. Neurosci Lett. 2007;422:228–31.
Ivanova SA, Loonen AJ, Pechlivanoglou P, Freidin MB, Al Hadithy AF, Rudikov EV, Zhukova IA, Govorin NV, Sorokina VA, Fedorenko OY, et al. NMDA receptor genotypes associated with the vulnerability to develop dyskinesia. Transl Psychiatry. 2012;2:e67.
Religa D, Czyzewski K, Styczynska M, Peplonska B, Lokk J, Chodakowska-Zebrowska M, Stepien K, Winblad B, Barcikowska M. Hyperhomocysteinemia and methylenetetrahydrofolate reductase polymorphism in patients with Parkinson's disease. Neurosci Lett. 2006;404:56–60.
Kaiser R, Hofer A, Grapengiesser A, Gasser T, Kupsch A, Roots I, Brockmoller J. L -dopa-induced adverse effects in PD and dopamine transporter gene polymorphism. Neurology. 2003;60:1750–5.
Strong JA, Dalvi A, Revilla FJ, Sahay A, Samaha FJ, Welge JA, Gong J, Gartner M, Yue X, L Y. Genotype and smoking history affect risk of levodopa-induced dyskinesias in Parkinson's disease. Mov Disord. 2006;21:654–9.
Wells K. Study quality assessment in systematic reviews of research on intervention effects. Res Soc Work Pract. 2009;19:52–62.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, Simonovic M, Roth A, Santos A, Tsafou KP, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43:D447–52.
Wang J, Duncan D, Shi Z, Zhang B. WEB-based GEne SeT AnaLysis toolkit (WebGestalt): update 2013. Nucleic Acids Res. 2013;41:W77–83.
Kim JS, Kim J-Y, Kim J-M, Kim JW, Chung SJ, Sung R, Kim RN, Mi JK, Kim H-T, Choi K-G, Shin D-I, et al. No correlation between COMT genotype and entacapone benefits in Parkinson’s disease. Neurol Asia. 2011;16:211–6.
Schumacher-Schuh AF, Francisconi C, Altmann V, Monte TL, Callegari-Jacques SM, Rieder CR, Hutz MH. Polymorphisms in the dopamine transporter gene are associated with visual hallucinations and levodopa equivalent dose in Brazilians with Parkinson's disease. Int J Neuropsychopharmacol. 2013:1–8.
Stefanovic M, Topic E, Ivanisevic AM, Relja M, Korsic M. Genotyping of CYP2D6 in Parkinson's disease. Clin Chem Lab Med. 2000;38:929–34.
Fujii C, Harada S, Ohkoshi N, Hayashi A, Yoshizawa K, Ishizuka C, Nakamura T. Association between polymorphism of the cholecystokinin gene and idiopathic Parkinson's disease. Clin Genet. 1999;56:394–9.
Yuan RY, Sheu JJ, Yu JM, Hu CJ, Tseng IJ, Ho CS, Yeh CY, Hung YL, Chiang TR. Methylenetetrahydrofolate reductase polymorphisms and plasma homocysteine in levodopa-treated and non-treated Parkinson's disease patients. J Neurol Sci. 2009;287:64–8.
Gorgone G, Curro M, Ferlazzo N, Parisi G, Parnetti L, Belcastro V, Tambasco N, Rossi A, Pisani F, Calabresi P, et al. Coenzyme Q10, hyperhomocysteinemia and MTHFR C677T polymorphism in levodopa-treated Parkinson's disease patients. NeuroMolecular Med. 2012;14:84–90.
Wu H, Dong F, Wang Y, Xiao Q, Yang Q, Zhao J, Quinn TJ, Chen SD, Liu J. Catechol-O-methyltransferase Val158Met polymorphism: modulation of wearing-off susceptibility in a Chinese cohort of Parkinson's disease. Parkinsonism Relat Disord. 2014;20:1094–6.
Goetz CG, Burke PF, Leurgans S, Berry-Kravis E, Blasucci LM, Raman R, Zhou L. Genetic variation analysis in parkinson disease patients with and without hallucinations: case-control study. Arch Neurol. 2001;58:209–13.
Pascale E, Purcaro C, Passarelli E, Guglielmi R, Vestri AR, Passarelli F, Meco G. Genetic polymorphism of Angiotensin-converting enzyme is not associated with the development of Parkinson's disease and of L-dopa-induced adverse effects. J Neurol Sci. 2009;276:18–21.
De Bonis ML, Tessitore A, Pellecchia MT, Longo K, Salvatore A, Russo A, Ingrosso D, Zappia V, Barone P, Galletti P, Tedeschi G. Impaired transmethylation potential in Parkinson's disease patients treated with L-Dopa. Neurosci Lett. 2010;468:287–91.
de Lau LM, Verbaan D, Marinus J, Heutink P, van Hilten JJ. Catechol-O-methyltransferase Val158Met and the risk of dyskinesias in Parkinson's disease. Mov Disord. 2012;27:132–5.
Paus S, Gadow F, Knapp M, Klein C, Klockgether T, Wullner U. Motor complications in patients form the German competence network on Parkinson's disease and the DRD3 Ser9Gly polymorphism. Mov Disord. 2009;24:1080–4.
Yahalom G, Kaplan N, Vituri A, Cohen OS, Inzelberg R, Kozlova E, Korczyn AD, Rosset S, Friedman E, Hassin-Baer S. Dyskinesias in patients with Parkinson's disease: effect of the leucine-rich repeat kinase 2 (LRRK2) G2019S mutation. Parkinsonism Relat Disord. 2012;18:1039–41.
Corvol JC, Bonnet C, Charbonnier-Beaupel F, Bonnet AM, Fievet MH, Bellanger A, Roze E, Meliksetyan G, Ben Djebara M, Hartmann A, et al. The COMT Val158Met polymorphism affects the response to entacapone in Parkinson's disease: a randomized crossover clinical trial. Ann Neurol. 2011;69:111–8.
Daniel SE, Lees AJ. Parkinson's disease society brain Bank, London: overview and research. J Neural Transm Suppl. 1993;39:165–72.
Gelb DJ, Oliver E, Gilman S. Diagnostic criteria for Parkinson disease. Arch Neurol. 1999;56:33–9.
Quinn NBR, Craufurd D, Goldman S, Hodges J, Kieburtz K, Lindvall O, MacMillan J, Roos R. Core assessment program for Intracerebral transplantation in Huntington's disease (CAPIT-HD). Mov Disord. 1996;11:143–50.
Cheshire P, Bertram K, Ling H, O'Sullivan SS, Halliday G, McLean C, Bras J, Foltynie T, Storey E, Williams DR. Influence of single nucleotide polymorphisms in COMT, MAO-A and BDNF genes on dyskinesias and levodopa use in Parkinson's disease. Neurodegener Dis. 2014;13:24–8.
Foltynie T, Cheeran B, Williams-Gray CH, Edwards MJ, Schneider SA, Weinberger D, Rothwell JC, Barker RA, Bhatia KP. BDNF val66met influences time to onset of levodopa induced dyskinesia in Parkinson's disease. J Neurol Neurosurg Psychiatry. 2009;80:141–4.
Greenbaum L, Goldwurm S, Zozulinsky P, Lifschytz T, Cohen OS, Yahalom G, Cilia R, Tesei S, Asselta R, Inzelberg R, et al. Do tardive dyskinesia and L-dopa induced dyskinesia share common genetic risk factors? An exploratory study. J Mol Neurosci. 2013;51:380–8.
Ziegler DA, Ashourian P, Wonderlick JS, Sarokhan AK, Prelec D, Scherzer CR, Corkin S. Motor impulsivity in Parkinson disease: associations with COMT and DRD2 polymorphisms. Scand J Psychol. 2014;55:278–86.
Tan EK, Cheah SY, Fook-Chong S, Yew K, Chandran VR, Lum SY, Yi Z. Functional COMT variant predicts response to high dose pyridoxine in Parkinson's disease. Am J Med Genet B Neuropsychiatr Genet. 2005;137b:1–4.
Xie T, Ho SL, Li LS, Ma OC. G/A1947 polymorphism in catechol-O-methyltransferase (COMT) gene in Parkinson's disease. Mov Disord. 1997;12:426–7.
Contin M, Martinelli P, Mochi M, Riva R, Albani F, Baruzzi A. Genetic polymorphism of catechol-O-methyltransferase and levodopa pharmacokinetic-pharmacodynamic pattern in patients with Parkinson's disease. Mov Disord. 2005;20:734–9.
Liu YZ, Tang BS, Yan XX, Liu J, Ouyang DS, Nie LN, Fan L, Li Z, Ji W, Hu DL, et al. Association of the DRD2 and DRD3 polymorphisms with response to pramipexole in Parkinson's disease patients. Eur J Clin Pharmacol. 2009;65:679–83.
Devos D, Lejeune S, Cormier-Dequaire F, Tahiri K, Charbonnier-Beaupel F, Rouaix N, Duhamel A, Sablonniere B, Bonnet AM, Bonnet C, et al. Dopa-decarboxylase gene polymorphisms affect the motor response to L-dopa in Parkinson's disease. Parkinsonism Relat Disord. 2014;20:170–5.
Bialecka M, Drozdzik M, Klodowska-Duda G, Honczarenko K, Gawronska-Szklarz B, Opala G, Stankiewicz J. The effect of monoamine oxidase B (MAOB) and catechol-O-methyltransferase (COMT) polymorphisms on levodopa therapy in patients with sporadic Parkinson's disease. Acta Neurol Scand. 2004;110:260–6.
Moreau C, Meguig S, Corvol JC, Labreuche J, Vasseur F, Duhamel A, Delval A, Bardyn T, Devedjian JC, Rouaix N, et al. Polymorphism of the dopamine transporter type 1 gene modifies the treatment response in Parkinson's disease. Brain. 2015;138:1271–83.
Lee MS, Lyoo CH, Ulmanen I, Syvanen AC, Rinne JO. Genotypes of catechol-O-methyltransferase and response to levodopa treatment in patients with Parkinson's disease. Neurosci Lett. 2001;298:131–4.
Ward CDGW. Research diagnostic criteria for Parkinson's disease. Adv Neurol. 1990;53:245–9.
Oliveri RL, Annesi G, Zappia M, Civitelli D, Montesanti R, Branca D, Nicoletti G, Spadafora P, Pasqua AA, Cittadella R, et al. Dopamine D2 receptor gene polymorphism and the risk of levodopa-induced dyskinesias in PD. Neurology. 1999;53:1425–30.
Ferrari M, Martignoni E, Blandini F, Riboldazzi G, Bono G, Marino F, Cosentino M. Association of UDP-glucuronosyltransferase 1A9 polymorphisms with adverse reactions to catechol-O-methyltransferase inhibitors in Parkinson's disease patients. Eur J Clin Pharmacol. 2012;68:1493–9.
Zappia M, Annesi G, Nicoletti G, Arabia G, Annesi F, Messina D, Pugliese P, Spadafora P, Tarantino P, Carrideo S, et al. Sex differences in clinical and genetic determinants of levodopa peak-dose dyskinesias in Parkinson disease: an exploratory study. Arch Neurol. 2005;62:601–5.
De Luca V, Annesi G, De Marco EV, de Bartolomeis A, Nicoletti G, Pugliese P, Muscettola G, Barone P, Quattrone A. HOMER1 promoter analysis in Parkinson's disease: association study with psychotic symptoms. Neuropsychobiology. 2009;59:239–45.
Acuna G, Foernzler D, Leong D, Rabbia M, Smit R, Dorflinger E, Gasser R, Hoh J, Ott J, Borroni E, et al. Pharmacogenetic analysis of adverse drug effect reveals genetic variant for susceptibility to liver toxicity. Pharmacogenomics J. 2002;2:327–34.
Lin JJ, Yueh KC, Lin SZ, Harn HJ, Liu JT. Genetic polymorphism of the angiotensin converting enzyme and L-dopa-induced adverse effects in Parkinson's disease. J Neurol Sci. 2007;252:130–4.
Wang J, Liu ZL, Chen B. Dopamine D5 receptor gene polymorphism and the risk of levodopa-induced motor fluctuations in patients with Parkinson's disease. Neurosci Lett. 2001;308:21–4.
Schumacher-Schuh AF, Altmann V, Rieck M, Tovo-Rodrigues L, Monte TL, Callegari-Jacques SM, Medeiros MS, Rieder CR, MH H: Association of common genetic variants of HOMER1 gene with levodopa adverse effects in Parkinson's disease patients. Pharmacogenomics J 2014, 14:289-294.
Rieck M, Schumacher-Schuh AF, Callegari-Jacques SM, Altmann V, Schneider Medeiros M, Rieder CR, MH H. Is there a role for ADORA2A polymorphisms in levodopa-induced dyskinesia in Parkinson's disease patients? Pharmacogenomics. 2015;16:573–82.
Rieck M, Schumacher-Schuh AF, Altmann V, Francisconi CL, Fagundes PT, Monte TL, Callegari-Jacques SM, Rieder CR, MH H. DRD2 haplotype is associated with dyskinesia induced by levodopa therapy in Parkinson's disease patients. Pharmacogenomics J. 2012;13:1701–10.
Molchadski I, Korczyn AD, Cohen OS, Katzav A, Nitzan Z, Chapman J, S H-B. The role of apolipoprotein E polymorphisms in levodopa-induced dyskinesia. Acta Neurol Scand. 2011;123:117–21.
Kaplan N, Vituri A, Korczyn AD, Cohen OS, Inzelberg R, Yahalom G, Kozlova E, Milgrom R, Laitman Y, Friedman E, et al. Sequence variants in SLC6A3, DRD2, and BDNF genes and time to levodopa-induced dyskinesias in Parkinson's disease. J Mol Neurosci. 2014;53:183–8.
Nalls MA, Pankratz N, Lill CM, Do CB, Hernandez DG, Saad M, DeStefano AL, Kara E, Bras J, Sharma M, et al. Large-scale meta-analysis of genome-wide association data identifies six new risk loci for Parkinson's disease. Nat Genet. 2014;46:989–93.
Lewis S. Neurodegenerative disorders: Parkinson's disease reveals hidden depths. Nat Rev Neurosci. 2013;14
Dreher J-C KP, Kolachana B, Weinberger DR, Berman KF. Variation in dopamine genes influences responsivity of the human reward system. Proc Natl Acad Sci. 2009;106:617–22.
Diatchenko LSG, Nackley AG, Bhalang K, Sigurdsson A, Belfer I, Goldman D, Xu K, Shabalina SA, Shagin D, Max MB, Makarov SS, Maixner W. Genetic basis for individual variations in pain perception and the development of a chronic pain condition. Hum Mol Genet. 2005;14:135–43.
Nackley AGSS, Tchivileva IE, Satterfield K, Korchynskyi O, Makarov SS, Maixner W, Diatchenko L. Human catechol-O-methyltransferase haplotypes modulate protein expression by altering mRNA secondary structure. Science. 2006;314:1930–3.
Habak C, Noreau A, Nagano-Saito A, Mejia-Constain B, Degroot C, Strafella AP, Chouinard S, Lafontaine AL, Rouleau GA, Monchi O. Dopamine transporter SLC6A3 genotype affects cortico-striatal activity of set-shifts in Parkinson's disease. Brain. 2014;137:3025–35.
Sanz-Blasco S, Avale E, Campana S, Damianich A, Saborido MD, Gomez G, Taravini I, Gershanik O, Ferrario JE. Exploring Fyn as a novel molecule in levodopa induced dyskinesias. Mov Disord. 2015;30
Wang YD, Bao XQ, Xu S, Yu WW, Cao SN, Hu JP, Li Y, Wang XL, Zhang D, Yu SS. A Novel Parkinson's Disease Drug Candidate with Potent Anti-neuroinflammatory Effects through the Src Signaling Pathway. J Med Chem. 2016;59:9062–79.
Lim KLTJ. Role of the ubiquitin proteasome system in Parkinson's disease. BMC Biochem. 2007;22
Brown AM, Deutch AY, Colbran RJ. Dopamine depletion alters phosphorylation of striatal proteins in a model of Parkinsonism. Eur J Neurosci. 2005;22(1):247–56.
Zhang SF, Xie CL, Wang Q, Liu ZG. Interactions of CaMKII with dopamine D2 receptors: roles in levodopa-induced dyskinesia in 6-hydroxydopamine lesioned Parkinson's rats. Sci Rep. 2015;4(1)
Donaldson JG. Multiple roles for Arf6: sorting, structuring, and signaling at the plasma membrane. J Biol Chem. 2003;278:41573–6.
Cho DI ZM, Min C, Kwon KJ, Shin CY, Choi HK, Kim KM. ARF6 and GASP-1 are post-endocytic sorting proteins selectively involved in the intracellular trafficking of dopamine D2 receptors mediated by GRK and PKC in transfected cells. Br J Pharmacol. 2013;168:1355–74.
Jan homolak ij, Maša filipovi. The role of IGF-1 in neurodegenerative diseases. Gyrus 2015, 3:162–167.
Felix P, Bernhard SH, Binder G, Weber K, Apel A, Roeben B, Deuschle C, Maechtel M, Heger T, Nussbaum S, Gasser T, Maetzler W, Berg D. Insulin-like growth factor 1 (IGF-1) in Parkinson's disease: potential as trait-, progression- and prediction marker and confounding factors. PLoS One. 2016;11:e0150552.
Schumacher-Schuh AF, Rieder CR, Hutz MH. Parkinson's disease pharmacogenomics: new findings and perspectives. Pharmacogenomics. 2014;15:1253–71.
Edwards TL, Scott WK, Almonte C, Burt A, Powell EH, Beecham GW, Wang L, Zuchner S, Konidari I, Wang G, et al. Genome-wide association study confirms SNPs in SNCA and the MAPT region as common risk factors for Parkinson disease. Ann Hum Genet. 2010;74:97–109.
Satake W, Nakabayashi Y, Mizuta I, Hirota Y, Ito C, Kubo M, Kawaguchi T, Tsunoda T, Watanabe M, Takeda A, et al. Genome-wide association study identifies common variants at four loci as genetic risk factors for Parkinson's disease. Nat Genet. 2009;41:1303–7.
Simon-Sanchez J, Schulte C, Bras JM, Sharma M, Gibbs JR, Berg D, Paisan-Ruiz C, Lichtner P, Scholz SW, Hernandez DG, et al. Genome-wide association study reveals genetic risk underlying Parkinson's disease. Nat Genet. 2009;41:1308–12.
Do CB, Tung JY, Dorfman E, Kiefer AK, Drabant EM, Francke U, Mountain JL, Goldman SM, Tanner CM, Langston JW, et al. Web-based genome-wide association study identifies two novel loci and a substantial genetic component for Parkinson's disease. PLoS Genet. 2011;7:e1002141.
Hill-Burns EM, Wissemann WT, Hamza TH, Factor SA, Zabetian CP, Payami H. Identification of a novel Parkinson's disease locus via stratified genome-wide association study. BMC Genomics. 2014;15:118.
Nalls MA, Plagnol V, Hernandez DG, Sharma M, Sheerin UM, Saad M, Simon-Sanchez J, Schulte C, Lesage S, Sveinbjornsdottir S, et al. Imputation of sequence variants for identification of genetic risks for Parkinson's disease: a meta-analysis of genome-wide association studies. Lancet. 2011;377:641–9.
Pankratz N, Beecham GW, DeStefano AL, Dawson TM, Doheny KF, Factor SA, Hamza TH, Hung AY, Hyman BT, Ivinson AJ, et al. Meta-analysis of Parkinson's disease: identification of a novel locus, RIT2. Ann Neurol. 2012;71:370–84.
Hamza TH, Zabetian CP, Tenesa A, Laederach A, Montimurro J, Yearout D, Kay DM, Doheny KF, Paschall J, Pugh E, et al. Common genetic variation in the HLA region is associated with late-onset sporadic Parkinson's disease. Nat Genet. 2010;42:781–5.
Spencer CC, Plagnol V, Strange A, Gardner M, Paisan-Ruiz C, Band G, Barker RA, Bellenguez C, Bhatia K, Blackburn H, et al. Dissection of the genetics of Parkinson's disease identifies an additional association 5′ of SNCA and multiple associated haplotypes at 17q21. Hum Mol Genet. 2011;20:345–53.
Vacic V, Ozelius LJ, Clark LN, Bar-Shira A, Gana-Weisz M, Gurevich T, Gusev A, Kedmi M, Kenny EE, Liu X, et al. Genome-wide mapping of IBD segments in an Ashkenazi PD cohort identifies associated haplotypes. Hum Mol Genet. 2014;23:4693–702.
We would like to acknowledge the Director, CSIR-IGIB, Dr. Sanjay Kumar, for his vision and support. The authors would like to thank, Nikhil Sahai, B.Tech student, Jaypee University of Information Technology, Himachal Pradesh for partially compiling the GWAS gene set and tabulate the same as a part of his dissertation.
Financial support from Council of Scientific and Industrial Research (CSIR) and Indian Council of Medical Research (ICMR) is duly acknowledged. Financial support from CSIR- funded projects - GENCODE-A (BSC0123) and GOMED (MLP1601) are duly acknowledged to carry out this study. DG acknowledges CSIR projects (BSC0123 and MLP1601), MKM acknowledges MLP0901, PT acknowledges CSIR (Research Associateship) and CR acknowledges UGC, Govt. of India for providing fellowship.
Availability of data and materials
The data supporting the results of this research paper are included within this article and its additional supplementary files.
Ethics approval and consent to participate
No patient samples were collected and analysed during this study. All GWAS data were provided as summary statistics by the consortia acknowledged in this study having been collected in accordance with ethical regulations in the partner countries and as defined in original research publications by such consortia.
Consent for publication
The authors declare no potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Including Supplementary Material such as Supplementary text S4, S7, Supplementary Tables S5–S6. (DOCX 34 kb)
Characteristics of included studies for assessment of association between genetic variants and ADRs in PD (DOCX 123 kb)
Characteristics of included studies for assessment of association between genetic variants and LR in PD (DOCX 52 kb)
GWAS studies on Parkinson’s disease genetics. (XLS 92 kb)
Functional enrichment analysis (by WebGestalt) result for each enriched function obtained using levodopa response associated genes. (XLS 48 kb)
Functional enrichment analysis (by WebGestalt) result for each enriched function obtained using Parkinson’s disease associated genes. (XLS 52 kb)