Impact Factor 2.889 | CiteScore 3.1

More on impact ›

Original Research ARTICLE

Front. Neurol., 29 September 2020 | https://doi.org/10.3389/fneur.2020.576465

Identification of Disease-Associated Variants by Targeted Gene Panel Resequencing in Parkinson's Disease

Kensuke Daida1, Manabu Funayama1,2,3*, Yuanzhe Li1, Hiroyo Yoshino2, Arisa Hayashida1, Aya Ikeda1, Kotaro Ogaki1, Kenya Nishioka1 and Nobutaka Hattori1,2,3
  • 1Department of Neurology, Juntendo University School of Medicine, Tokyo, Japan
  • 2Research Institute for Diseases of Old Age, Graduate School of Medicine, Juntendo University, Tokyo, Japan
  • 3Center for Genomic and Regenerative Medicine, Graduate School of Medicine, Juntendo University, Tokyo, Japan

Background: Recent advanced technologies, such as high-throughput sequencing, have enabled the identification of a broad spectrum of variants. Using targeted-gene-panel resequencing for Parkinson's disease (PD)-associated genes, we have occasionally found several single-nucleotide variants (SNVs), which are thought to be disease-associated, in PD patients. To confirm the significance of these potentially disease-associated variants, we performed genome association analyses, using next-generation target resequencing, to evaluate the associations between the identified SNVs and PD.

Methods: We obtained genomic DNA from 766 patients, who were clinically diagnosed with PD, and 336 healthy controls, all of Japanese origin. All data were analyzed using Ion AmpliSeq panel sequences, with 29 PD- or dementia-associated genes in a single panel. We excluded any variants that did not comply with the Hardy–Weinberg equilibrium in the control group. Variant frequencies in the PD and control groups were compared using PLINK. The identified variants were confirmed to a frequency difference of P < 0.05, after applying the Benjamini–Hochberg procedure using Fisher's exact test. The pathogenicity and prevalence of each variant were estimated based on a public gene database.

Results: We identified three rare variants that were significantly associated with PD: rs201012663/rs150500694 in SYNJ1 and rs372754391 in DJ-1, which are intronic variants, and rs7412 in ApoE, which is an exonic variant. The variants in SYNJ1 and ApoE were frequently identified in the control group, and rs201012663/rs150500694 in SYNJ1 may play a protective role against PD. The DJ-1 variant was frequently identified in the PD group, with a high odds ratio of 2.2.

Conclusion: The detected variants may represent genetic modifiers or disease-related variants in PD. Targeted-gene-panel resequencing may represent a useful method for detecting disease-causing variants and genetic association studies in PD.

Introduction

Parkinson's disease (PD) is the second-most frequent neurodegenerative disorder, associated with motor and non-motor symptoms (1). Clinical symptoms are characterized by tremor, rigidity, bradykinesia, and gait disturbances. To date, advanced genetic methods have revealed several genes associated with both familial and sporadic PD (2). Initially, genes were identified based on large pedigrees associated with Mendelian forms of PD, using positional cloning and linkage analyses, which resulted in the identification of SNCA (3), LRRK2 (4, 5), PRKN (6), and PINK1 (7). Later, next-generation sequencing (NGS) was used to identify additional causative genes, such as ATP13A2 (8), CHCHD2 (9), VPS13C (10), and PSAP (11). Furthermore, genome-wide association studies (GWASs) have identified single-nucleotide variants (SNVs) and other rare variants associated with sporadic PD (1214). The explained heritability ranged from 16 to 36%, even in the latest large GWAS, as reported in 2019 (15).

Many genome-association studies, including GWAS, have been conducted for PD; however, some unrevealed genetic background remains, referred to as “missing heritability” (13, 14, 16). Missing heritability is the difference between heritability estimated from twin studies and GWAS, as GWAS has only been able to detect some of the heritability estimated from twin studies (17). Many explanations for missing heritability have been proposed, including unrevealed variants with smaller effects, rarer variants that are poorly detected by the currently available genotyping arrays, copy number variants that cannot be detected by available arrays, and the low power to detect gene–gene interactions (16). Variants of GBA are known to be strong risk factors for sporadic PD but have not been detected by GWAS, likely due to a low minor allele frequency.

We have developed a targeted-gene-panel resequencing protocol to screen 29 PD-associated genes, simultaneously. Panel resequencing has both advantages and disadvantages because it can identify multiple types of variants, including pathogenic variants, risk-associated variants, and rare variants of uncertain significance. Therefore, determining which variants are disease-associated can be difficult. A previous report describing Mendelian genes showed that rare functional variants occurred more frequently in sporadic PD cases than in control cases, indicating that Mendelian genes may be associated not only with familial PD but also with sporadic PD, which may be assessable using panel resequencing (18). In our analyses, through targeted-gene-panel resequencing, rare variants were identified in ~40% of PD patients with a family history or early-onset PD (data not shown), and pathogenic variants were found in an even smaller percentage of patients. We also identified several putative disease-associated variants in PD patients. We hypothesized that these variants may play a role in PD onset and could account for some degree of missing heritability. Thus, we aimed to implement target-panel resequencing, to identify associations between SNVs and familial or early-onset PD. Our method contributes to expanding the understanding of missing heritability among familial and early-onset PD patients.

Materials and Methods

Participants

The present study was approved by the ethics committee of Juntendo University, Tokyo, Japan, and all participants provided written informed consent to participate in the genetic research. We collected DNA samples from the Juntendo PD DNA bank, which included 766 patients with PD, who were clinically diagnosed using standard criteria (1), and 336 healthy control subjects. Among these, 407 PD patients had a family history of PD (average age at onset: 54.6 ± 15.77 years, range 6–88), and the remaining 359 PD patients were without family history (average age at onset: 42.0 ± 11.22 years, range 9–83). We also collected data regarding the Hoehn and Yahr stages for each PD patient. The healthy controls were defined as individuals without any individual or family history of neurodegenerative disorders. An overview of the clinical characteristics of the included PD patients and healthy controls is shown in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Demographic data of the analyzed subjects.

Processing Data Output From the Ion Torrent System

The sequencing analysis of the Ion AmpliSeq panel (Thermo Fisher Scientific, Waltham, MA, USA) was performed using the Ion Chef System (Thermo Fisher Scientific) and the Ion S5 Sequencer (Thermo Fisher Scientific), according to the manufacturer's instructions. Our Ion AmpliSeq panel (Thermo Fisher Scientific, IAD103177_182) included 29 PD- and dementia-related genes (Table 2), and its coverage was 98.34% (829 amplicons, missed: 1,646 bp) (manuscript in preparation). The output data were obtained as a variant call format (VCF) file from the Ion torrent system. VCF files were processed using vcftools (19).

TABLE 2
www.frontiersin.org

Table 2. PD- and dementia-related genes analyzed by resequencing.

Statistical Analysis to Compare the Frequencies of Non-rare Variants

We confirmed all samples with a mean depth > 100 and excluded those amplicons with read depths smaller than 10. The analyzed variants were confirmed to exist among the target sequences and to have read depth of coverages >45. We also calculated the coverage percentage. During the variant-screening stage, we excluded all variants that did not comply with Hardy–Weinberg equilibrium (HWE; P < 0.05) within the control group (Figure 1). We analyzed only the control group during the variant-screening stage because performing HWE analysis while including PD patients would introduce bias. During the analysis stage, the variant frequencies observed for the PD and healthy non-PD groups were compared using PLINK 1.9 (20). To verify this comparison, variants with a frequency difference of P < 0.05, based on the performance of the Benjamini–Hochberg procedure and Fisher's exact test, were analyzed using the genotyping data available in 4.7KJPN, from the Japanese Multi Omics Reference Panel (jMorp) (21), and a genome aggregation database (gnomAD) (22). The scheme used for the analysis is presented in Figure 1. To confirm the presence of significant variants identified during the association study, we conducted Sanger sequencing on three cases with the variant and three cases without the variant, during the panel resequencing experiment.

FIGURE 1
www.frontiersin.org

Figure 1. Scheme of analysis. PD, Parkinson's disease; SYNJ1, Synaptojanin 1; APOE, Apolipoprotein E.

Results

The percentage of coverage was calculated, showing that 99.7% of the total dataset was read at a depth of least 1×, 98.9% at 20×, 97.9% at 100×, and 84% at 500×. During the variant screening stage, we identified 796 variants in our healthy controls, of which 749 were retained after screening for HWE compliance (P < 0.05) and were included in the association analysis performed using PLINK. We conducted Sanger sequencing on nine significant variants with p-values below 0.05 after performing the Benjamini–Hochberg procedure for Fisher's exact test, and five of them (chr1:65830299 T>G, chr1:65830300 T>G, chr3:184033555, chr2:233620927-233620929, and chr1:205743943) were not validated and excluded from the analysis. All of the false-positive variants were positioned around the tandem repeat of mononucleotides that was considered to cause false positives. Table 3 shows the top 15 variants that had the lowest p-values based on Fisher's exact test.

TABLE 3
www.frontiersin.org

Table 3. Top 15 variants detected in genome association studies.

Four variants were significantly associated with PD: rs201012663 and rs150500694 in SYNJ1, rs372754391 in DJ-1, and rs7412 in ApoE (Table 4). The two SYNJ1 variants, rs201012663 and rs150500694, were considered to represent a single variant because they are located four bases apart and demonstrated the same frequency in our subjects and public gene databases, which suggests that these variants are strongly linked (Tables 3, 4). The SYNJ1 variants are both located in an intron, with an odds ratio of 0.37. The DJ-1 variant (rs372754391) was also intronic and was more frequently identified in the PD cohort than in controls, with an odds ratio of 2.2. However, its frequency in the public database was quite large compared with the frequency in our data. The ApoE variant was exonic and was more frequently observed in the control group than in the PD group, with an odds ratio of 0.39. The ApoE variant was one of the single-nucleotide polymorphisms (SNPs) that determine the ApoE genotype. The E2 ApoE genotype was more frequently observed in the control group, whereas the E4 genotype was more frequently observed in the PD group (Table 5). No significant differences in age, age at onset, or Hoehn and Yahr scores were observed between patients with and without detected variants (Table 6).

TABLE 4
www.frontiersin.org

Table 4. Details of the detected variants with significant value.

TABLE 5
www.frontiersin.org

Table 5. Allele frequencies in patients, according to APOE genotype.

TABLE 6
www.frontiersin.org

Table 6. Clinical characteristics of patients, according to the presence of the identified variants.

We do not have data for four variants (rs16856139, rs11931532, rs11931074, and rs1994090) that were previously identified in a GWAS performed in Japanese PD patients because these variants were absent from our target panel (13). LRRK2 G2385R (rs34778348), which is a risk factor for PD in East Asian individuals, was the 21st most significant variant identified among our cohort (23). Except for rs34778348, none of the currently known risk variants for PD were detected.

Discussion

We performed a genetic case–control analysis, using NGS data from our Ion AmpliSeq panel. We identified three variants in three different genes: the combination of rs201012663 and rs150500694 in SYNJ1, rs372754391 in DJ-1, and rs7412 in ApoE. None of these three variants were reported as PD-related variants when we searched a GWAS catalog on June 8, 2020 (24). Our identified variants might account for missing heritability in PD. Targeted resequencing could perform deeper reads of selected genes associated with phenotypes than the microarrays that are normally used in GWAS. Thus, targeted resequencing-based association studies may be able to identify risk variants that have not been previously identified by GWAS (17).

The three identified variants have never previously been reported as variants associated with PD. In our study, variants in SYNJ1 (rs201012663 and rs150500694) showed a higher frequency in the control group than in the PD group. SYNJ1 is known to be a causative gene for early-onset Parkinsonism, with atypical characteristics, such as seizures, dystonia, and dementia, with an autosomal-recessive inheritance pattern (25, 26). This gene encodes the protein Synaptojanin 1, a polyphosphoinositide phosphatase that is concentrated at synapses (27, 28). Synaptojanin 1 is associated with synaptic vesicle endocytosis. The variants identified in SYNJ1 (rs201012663/rs150500694) in this study have not previously been reported to be pathogenic variants. Synaptojanin 1 is also known to play a role in the pathogenesis of Alzheimer's disease (AD), associated with a PI (4, 5)P2 imbalance. The haploinsufficiency of SYNJ1 protects cells from the neurotoxic actions of Aβ42 (29). The variants rs201012663/rs150500694 might play a similarly protective role against alpha synuclein-mediated neurotoxicity.

The identified variant in DJ-1 might be interesting, due to the high odds ratio of 2.2. However, this variant may be specific to ethnicity because the frequency of this variant among our healthy controls was lower than that observed in public databases. This variant was not recorded in jMorp, one of the largest genomic databases in Japan, suggesting its rarity in the Japanese population. DJ-1 was initially identified as an oncogene and was later found to cause familial PD (30). DJ-1 has also been associated with other disorders, including stroke, familial amyloidotic polyneuropathy, and type 2 diabetes (3033). DJ-1 has several functions, including transcriptional regulation, antioxidative stress reactions, chaperone, protease, and mitochondrial regulation (30). DJ-1 is expressed in almost all cells, including neurons and glial cells. DJ-1 protein contains three cysteine residues, C46, C56, and C106. C106 is likely to be influenced by oxidative stress and oxidized into SOH, SO2H, and SO3H (3436). DJ-1 containing a C106 residue that has been oxidized to SO3H is thought to represent an inactive form (37). In the brains of PD patients, excessively oxidized forms of DJ-1 have been observed (38). The identified mutation might facilitate oxidation, inactivating DJ-1.

APOE genotypes have previously been associated with an increased risk of AD (39, 40). rs7412 is one of two SNVs that have been defined in common allelic APOE variants. APOE4 is known to represent a strong risk factor for AD. The variant (rs7412) identified in our study is included in APOE1 or APOE2, which are known to decrease the risk of AD. rs7412 was significantly rare in the PD group in our study. In our study, APOE2 was significantly rare in the PD group, whereas APOE4 was significantly frequent in the PD group. Larger research studies have concluded that APOE epsilon had no association with PD onset (41). Differences between our study and past studies may be due to the smaller sample size included in our study and differences in the ethnicities of the participants.

In our study, SNVs detected in previous GWAS were not identified in our cohort because most of the reported risk-associated SNVs have been identified in non-coding regions, which were not included in our targeted panel (42). Targeted resequencing can cover more SNVs within the targeted exons than DNA microarrays, which are commonly used in GWAS. Our method might enable the detection of SNVs in exons or near exons that are not included in the SNP chips used for GWAS. Our target panel was designed to include all exons and the 25 bp up- and downstream of the exon–intron boundaries. Therefore, our method allowed the discovery of PD-related variants that were not detected by GWAS. The inclusion of patients with a family history or early-onset PD in our cohort might facilitate the detection of susceptibility-associated variants, with deep genetic backgrounds. For example, mutations in GBA are more frequently identified in familial PD patients than in sporadic PD patients (43). However, our panel resequencing approach also has several disadvantages. This approach cannot be used to identify novel genes associated with PD and does not cover the majority of introns and transcriptional regulatory regions. The variants detected in this study may also be associated with sporadic PD, similar to GWASs that identified causative genes associated with sporadic PD that were previously reported to be causative genes for familial PD (SNCA, MAPT, and LRRK2) (44).

Our study includes the following limitations: (i) the sample size is too small to satisfy genome-wide significance, (ii) the lack of a second cohort to confirm our results, (iii) the possibility of sampling bias in the control group because the allele frequencies of variants in the public database were different from those identified in our healthy control group, (iv) the absence of any functional analysis to support our results, and (v) the lack of copy number variant evaluations.

We developed a new approach for surveying susceptibility-associated variants by using targeted resequencing, which may represent an effective method for revealing hidden disease-associated variants. Further studies that include additional patients remain necessary to confirm the suitability of this approach for the identification of disease-associated variants.

Data Availability Statement

The DNA sequence data of 1,102 participants used in this study are based on the informed consent of genetic testing from all participants according to the formal procedure approved by the Juntendo University School of Medicine Ethics Committee. However, some participants have refused to publish their DNA sequence data in public databases. Therefore, if the reader wishes to use the raw data used in this paper, please request directly to the corresponding author Manabu Funayama, funayama@juntendo.ac.jp.

Ethics Statement

The studies involving human participants were reviewed and approved by Juntendo University School of Medicine Institutional Review Board (No. 2019227). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

KD designed and performed the experiments, analyzed the data, wrote the first manuscript, and revised the manuscript. MF designed the study, wrote the first manuscript, and revised the manuscript. YL, HY, AH, AI, KO, and KN performed the experiments, analyzed the data, and revised the manuscript. NH directed the research project and revised the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by JSPS KAKENHI (Grant Numbers 19K08003 to MF, 18K07536 to HY, 20K16504 to AI, 19K17047 to KO, 20K07893 to KN, and 18H04043 to NH), by the Brain/MINDS Beyond program from AMED (JP20dm0307024 to KO), and by the Advanced Genome Research and Bioinformatics Study to Facilitate Medical Innovation (GRIFIN) from AMED (20km0405206h0005 to NH).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was performed, in part, at the Intractable Disease Research Center Juntendo University Graduate School of Medicine. We thank Lisa Giles, PhD, from Edanz Group (https://en-author-services.edanzgroup.com/ac) for editing a draft of this manuscript.

References

1. Gibb WR, Lees AJ. The relevance of the Lewy body to the pathogenesis of idiopathic Parkinson's disease. J Neurol Neurosurg Psychiatry. (1988) 51:745–52. doi: 10.1136/jnnp.51.6.745

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Deng H, Wang P, Jankovic J. The genetics of Parkinson disease. Ageing Res Rev. (2017) 42:72–85. doi: 10.1016/j.arr.2017.12.007

CrossRef Full Text | Google Scholar

3. Polymeropoulos MH, Lavedan C, Leroy E, Ide SE, Dehejia A, Dutra A, et al. Mutation in the alpha-synuclein gene identified in families with Parkinson's disease. Science. (1997) 276:2045–7. doi: 10.1126/science.276.5321.2045

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Paisan-Ruiz C, Jain S, Evans EW, Gilks WP, Simon J, Van Der Brug M, et al. Cloning of the gene containing mutations that cause PARK8-linked Parkinson's disease. Neuron. (2004) 44:595–600. doi: 10.1016/j.neuron.2004.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Zimprich A, Biskup S, Leitner P, Lichtner P, Farrer M, Lincoln S, et al. Mutations in LRRK2 cause autosomal-dominant parkinsonism with pleomorphic pathology. Neuron. (2004) 44:601–7. doi: 10.1016/j.neuron.2004.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Kitada T, Asakawa S, Hattori N, Matsumine H, Yamamura Y, Minoshima S, et al. Mutations in the parkin gene cause autosomal recessive juvenile parkinsonism. Nature. (1998) 392:605–8. doi: 10.1038/33416

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Valente EM, Abou-Sleiman PM, Caputo V, Muqit MM, Harvey K, Gispert S, et al. Hereditary early-onset Parkinson's disease caused by mutations in PINK1. Science. (2004) 304:1158–60. doi: 10.1126/science.1096284

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Ramirez A, Heimbach A, Grundemann J, Stiller B, Hampshire D, Cid LP, et al. Hereditary parkinsonism with dementia is caused by mutations in ATP13A2, encoding a lysosomal type 5 P-type ATPase. Nat Genet. (2006) 38:1184–91. doi: 10.1038/ng1884

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Funayama M, Ohe K, Amo T, Furuya N, Yamaguchi J, Saiki S, et al. CHCHD2 mutations in autosomal dominant late-onset Parkinson's disease: a genome-wide linkage and sequencing study. Lancet Neurol. (2015) 14:274–82. doi: 10.1016/S1474-4422(14)70266-2

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Lesage S, Drouet V, Majounie E, Deramecourt V, Jacoupy M, Nicolas A, et al. Loss of VPS13C function in autosomal-recessive parkinsonism causes mitochondrial dysfunction and increases PINK1/Parkin-dependent mitophagy. Am J Hum Genet. (2016) 98:500–13. doi: 10.1016/j.ajhg.2016.01.014

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Oji Y, Hatano T, Ueno SI, Funayama M, Ishikawa KI, Okuzumi A, et al. Variants in saposin D domain of prosaposin gene linked to Parkinson's disease. Brain. (2020) 143:1190–205. doi: 10.1093/brain/awaa064

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Maraganore DM, De Andrade M, Lesnick TG, Strain KJ, Farrer MJ, Rocca WA, et al. High-resolution whole-genome association study of Parkinson disease. Am J Hum Genet. (2005) 77:685–93. doi: 10.1086/496902

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Satake W, Nakabayashi Y, Mizuta I, Hirota Y, Ito C, Kubo M, 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. doi: 10.1038/ng.485

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Chang D, Nalls MA, Hallgrimsdottir IB, Hunkapiller J, Van Der Brug M, Cai F, et al. A meta-analysis of genome-wide association studies identifies 17 new Parkinson's disease risk loci. Nat Genet. (2017) 49:1511–6. doi: 10.1038/ng.3955

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Nalls MA, Blauwendraat C, Vallerga CL, Heilbron K, Bandres-Ciga S, Chang D, et al. Identification of novel risk loci, causal insights, and heritable risk for Parkinson's disease: a meta-analysis of genome-wide association studies. Lancet Neurol. (2019) 18:1091–102. doi: 10.1016/S1474-4422(19)30320-5

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, et al. Finding the missing heritability of complex diseases. Nature. (2009) 461:747–53. doi: 10.1038/nature08494

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Maher B. Personal genomes: the case of the missing heritability. Nature. (2008) 456:18–21. doi: 10.1038/456018a

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Spataro N, Calafell F, Cervera-Carles L, Casals F, Pagonabarraga J, Pascual-Sedano B, et al. Mendelian genes for Parkinson's disease contribute to the sporadic forms of the disease. Hum Mol Genet. (2015) 24:2023–34. doi: 10.1093/hmg/ddu616

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, Depristo MA, et al. The variant call format and VCFtools. Bioinformatics. (2011) 27:2156–8. doi: 10.1093/bioinformatics/btr330

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. (2007) 81:559–75. doi: 10.1086/519795

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Tadaka S, Saigusa D, Motoike IN, Inoue J, Aoki Y, Shirota M, et al. jMorp: Japanese multi omics reference panel. Nucl Acids Res. (2018) 46:D551–7. doi: 10.1093/nar/gkx978

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. (2020) 581:434–43. doi: 10.1038/s41586-020-2308-7

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Funayama M, Li Y, Tomiyama H, Yoshino H, Imamichi Y, Yamamoto M, et al. Leucine-rich repeat kinase 2 G2385R variant is a risk factor for Parkinson disease in Asian population. Neuroreport. (2007) 18:273–5. doi: 10.1097/WNR.0b013e32801254b6

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Buniello A, Macarthur JAL, Cerezo M, Harris LW, Hayhurst J, Malangone C, et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucl Acids Res. (2019) 47:D1005–D12. doi: 10.1093/nar/gky1120

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Krebs CE, Karkheiran S, Powell JC, Cao M, Makarov V, Darvish H, et al. The Sac1 domain of SYNJ1 identified mutated in a family with early-onset progressive Parkinsonism with generalized seizures. Hum Mutat. (2013) 34:1200–7. doi: 10.1002/humu.22372

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Quadri M, Fang M, Picillo M, Olgiati S, Breedveld GJ, Graafland J, et al. Mutation in the SYNJ1 gene associated with autosomal recessive, early-onset Parkinsonism. Hum Mutat. (2013) 34:1208–15. doi: 10.1002/humu.22373

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Mcpherson PS, Garcia EP, Slepnev VI, David C, Zhang X, Grabs D, et al. A presynaptic inositol-5-phosphatase. Nature. (1996) 379:353–7. doi: 10.1038/379353a0

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Mani M, Lee SY, Lucast L, Cremona O, Di Paolo G, De Camilli P, et al. The dual phosphatase activity of synaptojanin1 is required for both efficient synaptic vesicle endocytosis and reavailability at nerve terminals. Neuron. (2007) 56:1004–18. doi: 10.1016/j.neuron.2007.10.032

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Berman DE, Dall'armi C, Voronov SV, Mcintire LB, Zhang H, Moore AZ, et al. Oligomeric amyloid-beta peptide disrupts phosphatidylinositol-4,5-bisphosphate metabolism. Nat Neurosci. (2008) 11:547–54. doi: 10.1038/nn.2100

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Ariga H, Takahashi-Niki K, Kato I, Maita H, Niki T, Iguchi-Ariga SM. Neuroprotective function of DJ-1 in Parkinson's disease. Oxid Med Cell Longev. (2013) 2013:683920. doi: 10.1155/2013/683920

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Aleyasin H, Rousseaux MW, Phillips M, Kim RH, Bland RJ, Callaghan S, et al. The Parkinson's disease gene DJ-1 is also a key regulator of stroke-induced damage. Proc Natl Acad Sci USA. (2007) 104:18748–53. doi: 10.1073/pnas.0709379104

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Koide-Yoshida S, Niki T, Ueda M, Himeno S, Taira T, Iguchi-Ariga SM, et al. DJ-1 degrades transthyretin and an inactive form of DJ-1 is secreted in familial amyloidotic polyneuropathy. Int J Mol Med. (2007) 19:885–93. doi: 10.3892/ijmm.19.6.885

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Jain D, Jain R, Eberhard D, Eglinger J, Bugliani M, Piemonti L, et al. Age- and diet-dependent requirement of DJ-1 for glucose homeostasis in mice with implications for human type 2 diabetes. J Mol Cell Biol. (2012) 4:221–30. doi: 10.1093/jmcb/mjs025

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Mitsumoto A, Nakagawa Y, Takeuchi A, Okawa K, Iwamatsu A, Takanezawa Y. Oxidized forms of peroxiredoxins and DJ-1 on two-dimensional gels increased in response to sublethal levels of paraquat. Free Radic Res. (2001) 35:301–10. doi: 10.1080/10715760100300831

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Canet-Avilés RM, Wilson MA, Miller DW, Ahmad R, Mclendon C, Bandyopadhyay S, et al. The Parkinson's disease protein DJ-1 is neuroprotective due to cysteine-sulfinic acid-driven mitochondrial localization. Proc Natl Acad Sci USA. (2004) 101:9103–8. doi: 10.1073/pnas.0402959101

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Kinumi T, Kimata J, Taira T, Ariga H, Niki E. Cysteine-106 of DJ-1 is the most sensitive cysteine residue to hydrogen peroxide-mediated oxidation in vivo in human umbilical vein endothelial cells. Biochem Biophys Res Commun. (2004) 317:722–8. doi: 10.1016/j.bbrc.2004.03.110

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Zhou W, Zhu M, Wilson MA, Petsko GA, Fink AL. The oxidation state of DJ-1 regulates its chaperone activity toward alpha-synuclein. J Mol Biol. (2006) 356:1036–48. doi: 10.1016/j.jmb.2005.12.030

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Bandopadhyay R, Kingsbury AE, Cookson MR, Reid AR, Evans IM, Hope AD, et al. The expression of DJ-1 (PARK7) in normal human CNS and idiopathic Parkinson's disease. Brain. (2004) 127:420–30. doi: 10.1093/brain/awh054

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Corder EH, Saunders AM, Strittmatter WJ, Schmechel DE, Gaskell PC, Small GW, et al. Gene dose of apolipoprotein E type 4 allele and the risk of Alzheimer's disease in late onset families. Science. (1993) 261:921–3. doi: 10.1126/science.8346443

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Verghese PB, Castellano JM, Holtzman DM. Apolipoprotein E in Alzheimer's disease and other neurological disorders. Lancet Neurol. (2011) 10:241–52. doi: 10.1016/S1474-4422(10)70325-2

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Federoff M, Jimenez-Rolando B, Nalls MA, Singleton AB. A large study reveals no association between APOE and Parkinson's disease. Neurobiol Dis. (2012) 46:389–92. doi: 10.1016/j.nbd.2012.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Nalls MA, Pankratz N, Lill CM, Do CB, Hernandez DG, Saad 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. doi: 10.1038/ng.3043

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Lesage S, Anheim M, Condroyer C, Pollak P, Durif F, Dupuits C, et al. Large-scale screening of the Gaucher's disease-related glucocerebrosidase gene in Europeans with Parkinson's disease. Hum Mol Genet. (2011) 20:202–10. doi: 10.1093/hmg/ddq454

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Pankratz N, Wilk JB, Latourelle JC, Destefano AL, Halter C, Pugh EW, et al. Genomewide association study for susceptibility genes contributing to familial Parkinson disease. Hum Genet. (2009) 124:593–605. doi: 10.1007/s00439-008-0582-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: next generation sequencing, panel resequencing, genetic association study, Parkinsion's disease, missing heritability

Citation: Daida K, Funayama M, Li Y, Yoshino H, Hayashida A, Ikeda A, Ogaki K, Nishioka K and Hattori N (2020) Identification of Disease-Associated Variants by Targeted Gene Panel Resequencing in Parkinson's Disease. Front. Neurol. 11:576465. doi: 10.3389/fneur.2020.576465

Received: 26 June 2020; Accepted: 20 August 2020;
Published: 29 September 2020.

Edited by:

Chin-Hsien Lin, National Taiwan University Hospital, Taiwan

Reviewed by:

Beisha Tang, Central South University, China
Bruno Lopes Santos-Lobato, Federal University of Pará, Brazil

Copyright © 2020 Daida, Funayama, Li, Yoshino, Hayashida, Ikeda, Ogaki, Nishioka and Hattori. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Manabu Funayama, funayama@juntendo.ac.jp