- Open Access
Characterization of a pathogenic variant in GBA for Parkinson’s disease with mild cognitive impairment patients
Molecular Brain volume 13, Article number: 102 (2020)
Parkinson’s disease (PD) is the second most common neurodegenerative disease, and mild cognitive impairment (MCI) is a well-established risk factor for the development of dementia in PD. A growing body of evidence suggests that low expression of glucocerebrosidase (GBA) promotes the transmission of α-synuclein (α-Syn) interpolymers and the progression of PD. However, how GBA mutations affect the pathogenesis of PD via abnormal aggregation of α-Syn is unclear, and no clinically valid PD-MCI genetic markers have been identified. Here, we first located a GBA eQTL, rs12411216, by analysing DHS, eQTL SNP, and transcription factor binding site data using the UCSC database. Subsequently, we found that rs12411216 was significantly associated with PD-MCI (P < 0.05) in 306 PD patients by genotyping. In exploring the relationship between rs12411216 and GBA expression, the SNP was found to be associated with GBA expression in 50 PD patients through qPCR verification. In a further CRISPR/Cas9-mediated genome editing module, the SNP was identified to cause a decrease in GBA expression, weaken enzymatic activity and enhance the abnormal aggregation of α-Syn in SH-SY5Y cells. Additionally, using an electrophoretic mobility shift assay, we confirmed that the binding efficiency of transcription factor E2F4 was affected by the rs12411216 SNP. In conclusion, our results showed that rs12411216 regulated GBA expression, supporting its potential role as a PD-MCI genetic biomarker and highlighting novel mechanisms underlying Parkinson’s disease.
Parkinson’s disease (PD) is the second most common degenerative disease after Alzheimer’s disease. It involves the degeneration of dopaminergic neurons in the substantia nigra pars compacta and diffuse Lewy body deposition. The clinical manifestations of PD include resting tremor, bradykinesia, myotonia, abnormal posture and gait disorder. Additionally, this disease is accompanied by non-motor symptoms such as depression, constipation and sleep disorders, which impose a heavy burden on patients, families and society . The motor symptoms of PD included motor retardation, tremor, and myotonia. Both posture and gait abnormalities are well documented, but non-motor symptoms are often overlooked.
Cognitive impairments, including Parkinson’s disease-mild cognitive impairment (MCI) and Parkinson’s disease dementia (PDD), were the most common and disabling non-motor symptoms . MCI is an early warning signal of late dementia in PD. Research has shown that nearly 30% of newly diagnosed PD patients have MCI, more than 40% of PD patients with normal cognitive function will develop MCI within 6 years, and nearly 80% of patients may develop PDD at the late stage of PD . Moreover, the early diagnosis and prevention of PD-MCI can prevent the development of PDD. However, effective clinical genetic markers or targets for early diagnosis are currently not available.
The GBA-encoded protein glucocerebrosidase can decompose glucocerebroside into glucose and ceramide. Recent studies showed that low GBA expression promoted the prion-like spread of α-Syn interpolymer complexes, facilitated the pathogenesis of PD and increased cognitive damage [3, 4]. GBA gene expression was high in the brain tissue of healthy individuals but was decreased in that of PD patients [3, 5]. Genetic and clinical studies confirmed that GBA mutations were more prevalent in PD patients than in other populations [6, 7]. Similar to the common N370S and L444P mutations, GBA heterozygous mutations that were associated with the onset of PD were defined as deleterious [8, 9]. In most PD populations, the GBA mutation frequency was higher than that of other PD-related genes. A study consisting of 5691 PD patients and 4988 controls showed that the odds ratio of GBA mutations (L444P & N370S) in PD patients was 5.43 (95% confidence interval: 3.89–7.57), indicating that GBA mutations were important factors in PD development .
Regardless of ethnicity, age, or sex, PD patients had a higher frequency of GBA mutations than the control individuals. In a study containing 99 patients with idiopathic PD from Israel and 1543 Jewish controls from Germany, a significantly higher proportion of Parkinson’s disease patients (31.3%) were found to carry one or more GBA mutations than control groups (6.2%, P < 0.001) by screening six common GBA mutations . Similar findings could be observed from studies of different ethnic groups in other regions . Furthermore, a large US cohort study revealed that the GBA mutation frequency in familial PD was 4.1%, while in the control population, only 1.1% of individuals harboured such mutations . Taken together, for both idiopathic and familial Parkinson’s disease, the frequency of GBA mutations was higher in case groups than in control populations . Accordingly, a recent study showed that knocking out the GBA gene in human neuroblastoma cells (SH-SY5Y cells) led to neuronal glycosphingolipid accumulation .
Current studies have shown that abnormal aggregation of α-Syn is a significant pathological feature of PD, and the expression of SNCA, a gene encoding α-Syn, is directly related to the pathogenesis of PD [12, 14]. Additionally, the decrease in the proportion of α-Syn tetramers and the increase in the proportion of α-Syn monomers (which can form abnormal aggregates) in cells could easily lead to the onset of PD [13, 15]. Moreover, the loss of GBA function increased α-Syn levels, promoted its prion-like spread and aggregation and aggravated cognitive impairment in PD [16, 17]. However, the relationship between GBA mutations and abnormal α-Syn aggregation and their roles in the pathogenesis of PD remains unclear.
In this study, by retrieving data from the UCSC database, we hypothesized that the rs12411216 SNP might regulate the functions of the GBA gene. The rs12411216 SNP was found to overlap with the core motif of transcription factor E2F4 and served as the candidate for downstream analysis. For our downstream validation, a total of 306 PD outpatients were recruited from two hospitals for cognitive impairment assessment and genotyping. GBA expression, enzymatic activity, and the aggregation of α-Syn were measured in CRISPR/Cas9-edited SH-SY5Y cells. Mechanistically, electrophoretic mobility shift assays (EMSAs) were applied to preliminarily explore declines in transcription factor binding.
Materials and methods
Definition of GBA regulatory regions
Using the NCBI database, we obtained the coordinates of the GBA gene in the human genome sequence (hg19 version). The regulatory regions of the GBA gene were delineated as sequences ±50 kb away from the transcriptional start site.
Identification of DHSs potentially regulating GBA
Total DHS data (track: wgEncodeRegDnaseClusteredV3) were obtained from the UCSC website. All these data were converted to hg19 version using the UCSC liftover tool for consistency. Single base-pair DHS data were removed as they were not recognized in subsequent processing procedures. The intersection of the processed DHS data was retained, and regulatory regions of the GBA gene were decided using BEDOPS with the “bedops -element-of 1” parameter. Of these identified DHSs, those that reacted to the GBA promoter were selected according to Hi-C data, which were obtained from a previous study . A total of 143 DHSs located in the regulatory regions of GBA were identified and termed GBA-related DHSs that potentially regulated GBA expression.
Biological analysis of rs12411216 regulatory function
eQTL data (filename: gtexEqtlCluster) were downloaded from the UCSC website (version hg19. 8), and eQTL SNPs lying within GBA-related DHSs were selected as primary candidates. To examine their regulatory function in GBA, after searching against the transcription factor binding site database (filename: wgEncodeRegTfbsClustered), the rs12411216 SNP was found to overlap with the core motif of transcription factor E2F4 and served as the candidate for downstream analysis.
A total of 306 PD patients from Guangdong Provincial People's Hospital (259 patients) and Tianjin Huanhu Hospital (47 patients) underwent an rs12411216 SNP typing test and a “PD-NC&PD-MCI typing test”. Global cognitive and functional status was assessed using the Mini-Mental State Examination (MMSE) and the Montreal Cognitive Assessment (MOCA). All the patients were assessed with the Animal Fluency Test and Picture Arrangement Test (Wechsler Adult Intelligence Scale, WAIS-R) for executive functions, the Digit Span and Digit Symbol tests for attention and working memory, the Logical Memory Test and Immediate Memory Test (Wechsler Memory Scale WMS) for memory, the Vocabulary and Similarities Test (WAIS-R) for language, and the Block Design and Object Assembly tests for visuospatial function.
Patients with PD were divided into those with normal cognition (PD-NC) and those with mild cognitive impairment (PD-MCI) according to the Level II criteria of the 2012 Movement Disorder Society (MDS) Task Force Guidelines. A total of 50 PD patients from Guangdong Provincial People's Hospital underwent genotyping for the rs12411216 SNP and the GBA expression test.
CRISPR/Cas9-mediated knockout of the DHS and the replacement of rs12411216
Based on the DNase hypersensitive site (DHS) sequence of PD (chr12: 131621602–131,623,539), we designed single guide RNAs (sgRNAs) to knock out this regulatory region in the SH-SY5Y cell line. pCMV-Cas9, a vector that has a selectable neomycin marker, was obtained from Addgene (41815). The following sgRNAs were designed via the Zhang Laboratory website (https://zlab.bio/guide-design-resources) for the DHS-knockout experiment: sgRNA1, 5′-GGAGTGCAATGGCGCGATCTCGG-3′ and sgRNA2, 5′- CTCCGTCTCAAAACAAAAAAGGG-3′. The following sgRNA was designed for the GBA-SNP replacement experiment: SNP-sgRNA: GTGCGGCGAGACCCTGGGGCAGG, along with donor DNA (Supplementary Material). These sgRNAs were then cloned into U6-sgRNA plasmid vectors. A total of 1 μg of DNA (0.25 μg sgRNA1 vector or SNP-sgRNA vector + 0.25 μg sgRNA2 vector or donor DNA vector + 0.5 μg pCMV-Cas9-vector) was diluted in 50 μl Dulbecco’s modified Eagle’s medium (DMEM). Next, 3 μl GenJet In Vitro DNA Transfection Reagent (SL100489; SignaGen Laboratories) was diluted in 50 μl Dulbecco’s modified Eagle’s medium (DMEM), and two solutions were mixed and incubated for 20 min. After that, sgRNA1, sgRNA2 and pCMV-Cas9 vectors (0.25 μg sgRNA1 vector + 0.25 μg sgRNA2 vector + 0.5 μg pCMV-Cas9-vector) were cotransfected into the SH-SY5Y cell line using GenJet In Vitro DNA Transfection Reagent. To select cells expressing the neomycin resistance protein, the antibiotic G418 (800 μg/ml) was added to the medium and DMEM 48 h after the transfection, with G418 thereafter replaced every 48 h. After the addition of G418 for 8 days, single colonies were selected using the limiting dilution and verified by both PCR and sequencing methods .
GCase enzyme activity assay
GCase enzyme activity was assessed using the following methods. First, enriched lysosomes in 50 μl of activity assay buffer were added to a 96-well plate. This composition of this buffer was 0.25% Triton X-100 (TX; Sigma-Aldrich), 0.25% taurocholic acid (Sigma-Aldrich) and 1 mM ethylenediaminetetraacetic acid (EDTA) in citrate-phosphate buffer (pH 5.4). Second, 50 μl of 1% bovine serum albumin with 1 mM of the GBA substrate 4-methylumbelliferyl β-glucoside (M3633; Sigma-Aldrich) and/or 10 mM of the GBA inhibitor conduritol B epoxide (Sigma-Aldrich) were added, followed by incubation at 37 °C for 40 min. Last, 50 μl of 1 M glycine (pH 12.5) was added to terminate the reaction, and the fluorescence was measured using a plate reader (excitation = 355 nm, emission = 460 nm, 0.1 s; PerkinElmer) .
Electrophoretic mobility shift assay (EMSA)
First, nuclear extract was isolated from SH-SY5Y cells using Biyuntian kit P0028. Second, 3 μg of nucleoprotein was mixed with a biotin-labelled probe (Supplementary Material) containing rs12411216 (50 fmol) and incubated for 20 min at 25 °C. Third, the mixture (some of which was attached to the biotin-labelled probe) was separated using 6% polyacrylamide gel in Tris-glycine-EDTA running buffer at 4 °C for 2 h. Last, molecules attached on the gel were then transferred and assessed using an enhanced chemiluminescence detection system (Sage Creation, Beijing, China) .
We first placed SH-SY5Y cells on coverslips with poly-D-lysine coating. Next, 4% paraformaldehyde was used to fix cells, followed by blocking with 5% normal donkey serum (Jackson ImmunoResearch, Bar Harbor, ME, USA) and 0.1% TX (Sigma-Aldrich) for 1 h at room temperature. To assess the expression levels of p-α-Syn, a series of incubations with the primary antibody anti-p-α-Syn (1:1000, Abcam) were conducted at 4 °C overnight. The samples were washed off with 0.1% TX in phosphate-buffered saline, followed by 1 h of incubation with a mixture of fluorescein isothiocyanate (FITC)-conjugated secondary antibody (Jackson ImmunoResearch Laboratories, Inc.) and cyanine 3 (CY3)-conjugated secondary antibody (Jackson ImmunoResearch Laboratories, Inc.) at room temperature. Fluorescent images were acquired using a confocal microscope (Zeiss Confocal LSM 710) after coverslips were mounted. All images were processed using Zeiss Zen software. The range of signal intensity of the threshold in the selected areas was measured using ImageJ software.
Quantitative PCR (qPCR)
GBA and SNCA primer pairs are listed in the Supplementary Material. The β-actin gene was used as a reference gene for the qPCR analysis. Reagents for qPCR were obtained from Takara Biotechnology (DRR096A; Dalian, China). Expression data of GBA and SNCA genes in both wild-type SH-SY5Y (eight replicates) and knockout SH-SY5Y cells (six replicates) were collected for downstream analysis. Relative expression was calculated using the following formula: relative expression = 2−ΔΔCt; the relative expression was normalized based on the expression level of knockout SH-SY5Y cells. P-values for the differences between the wild-type and knockout cell groups were calculated using the t-test.
SH-SY5Y cells were lysed in 300 μl of radioimmunoprecipitation assay buffer (Biotech Well; Shanghai, China) in a 96-well plate containing 1 mM phenylmethylsulfonyl fluoride (PMSF) per well. Total protein concentrations were determined by bicinchoninic acid assays (Biotech Well). Samples were normalized with lysis buffer to ensure that each sample had the same protein concentration. Next, samples were mixed with an equal volume of 2× Laemmli sample buffer and solubilized by boiling for 10 min at 99 °C. Proteins were separated by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), and the proteins were then detected using mouse monoclonal antibodies against GBA, α-Syn and phospho-Ser129 α-Syn (p-α-Syn, a marker of pathologic α-Syn) at the dilution recommended by the manufacturer (EarthOx, San Francisco, CA, USA). β-actin was detected using a mouse monoclonal antibody against β-actin (EarthOx, San Francisco, CA, USA).
The chi-squared test was used to verify the association between rs12411216 and PD-MCI in clinical samples. The data are depicted as means ± standard deviations and were analysed using an unpaired Student’s t-test. A P value smaller than 0.05 was defined as statistically significant. All graphs and statistical calculations were performed using GraphPad Prism (Version 6.0). Other supplementary materials were available online.
The rs12411216 SNP was significantly associated with PD-MCI
In exploring the presence of a regulatory SNP of the GBA gene, the eQTL SNP rs12411216 was identified as a candidate in the distal regulatory element of the GBA gene using the UCSC database (Fig. 1a). Allelic and genotypic comparisons between 122 PD-MCI patients and 184 PD-NC patients who had Parkinson’s disease but no cognitive impairment showed that 5.56% (n = 17) of PD patients (n = 306) had PD-MCI and a homozygous (A/A) rs12411216 genotype, while 16.67% (n = 51) had PD-MCI and a heterozygous (C/A) rs12411216 genotype. The allele frequencies among the PD-MCI patients were 34.84% for A and 65.16% for C and were consistent with the 1000 Genomes Project datasets. The rs12411216 genotype was significantly associated with PD-MCI (chi-squared test, P = 0.0478).
The rs12411216 SNP was significantly associated with GBA expression in PD patients
After genotyping 50 PD patients from Guangdong Provincial People's Hospital, we found that, consistent with the 1000 Genomes Project data, GBA allele frequencies were 72.00% for C/C, 22.00% for C/A, and 6.00% for A/A. Additionally, the GBA mRNA level was significantly decreased in patients with homozygous mutations (AA genotype) compared to patients with the wild-type genotype (CC) (Fig. 1b).
The rs12411216 SNP affected the expression of GBA by decreasing the binding efficiency of transcription factor E2F4
To verify whether the rs12411216 SNP affected GBA expression in vivo, the rs12411216 SNP was replaced, and the upstream DHS sequence of GBA (chr1: 155164145–155,164,895) was knocked out in SH-SY5Y cells using the CRISPR/Cas9-mediated genome editing method (Fig. 2a). Our results showed that the DHS sequence significantly affected the expression of GBA (Fig. 2b–c), and rs12411216 was a functional SNP of the GBA gene (Fig. 2e–g). To determine whether GBA mutation affected the expression of SNCA, an important gene implicated in PD, we investigated the mRNA expression of SNCA and the protein expression of α-Syn and p-α-Syn in the knockout and replacement cells (Fig. 2h). We found no significant difference in the mRNA or protein expression of the SNCA gene, but the level of p-α-Syn was significantly increased both in the knockout and in the replacement cells. Together, these results indicated that rs12411216 was a regulatory SNP of the GBA gene and that the DHS fragment might serve as a regulatory element affecting the expression of the GBA gene.
To uncover how the regulatory SNP affects the expression of the GBA gene, we further analysed the presence of a core binding motif of transcription factor E2F4 within this DHS fragment. Based on the EMSA, we observed that when the wild-type SNP biotin-labelled probe (rs12411216: C) was applied, the probe bound to the nuclear protein, and a shifted band could be observed (Fig. 3, Line 2). After adding E2F4 antibody, a super-shift band could also be observed (Fig. 3, Line 4). However, when using the mutant SNP biotin-labelled probe (rs12411216: A), the binding ability of the probe to the nuclear protein was much weaker (Fig. 3, Line 6) than that of the wild-type biotin-labelled probe, and the binding ability to E2F4 was also reduced (Fig. 3, Line 8). In conclusion, our results showed that the rs12411216 mutation decreased the binding capacity of transcription factor E2F4.
rs12411216 affects GCase enzyme activity and α-Syn aggregation
As shown in Fig. 2d, GCase enzyme activity was significantly decreased in both the DHS knockout and rs12411216 replacement cells compared to that in the wild-type cells. By using anti-p-Syn antibody in a cellular immunofluorescence analysis, we observed that intracellular p-Syn levels were increased by 2.13- and 1.52-fold after the knockout of the DHS and the replacement of the SNP, respectively (Fig. 4). The pathological aggregation of α-Syn was more severe in cells with rs12411216 replacement. These results indicated that replacing rs12411216 could cause pathological aggregation of α-Syn and increase the risk of PD.
This study chose the GBA gene, because mutations in the GBA gene are recognized as a strong risk factor for the development of PD and dementia with Lewy bodies (DLB). The reason for choosing the rs1241126 SNP was that the rs12411216 SNP is located near the core binding domain of transcription factor E2F4 in the regulatory element of the GBA gene (the corresponding DHS).
Our screening process and specific criteria for the DHS and SNP have been added to the discussion section, and the details are listed as follows. In this study, DNase-seq data in the UCSC database were used to obtain chromatin opening (DHSs) sites of Parkinson’s disease-related genes. These DHSs may regulate the expression of Parkinson’s disease-related genes (such as SNCA and GBA). Though the DHS data revealed several potential regulatory regions of GBA, they cannot directly suggest the phenotypic changes. Therefore, we filtered the DHS data combining Hi-C data and eQTL SNP data, which provide chromatin reaction data and the connection between SNPs and gene expression.
First, after mining the Hi-C data of Parkinson’s disease from the ENCODE database, we obtained the gene locations of chromatin that directly acted on the promoter of Parkinson’s disease. Based on the binding sites screened in the DHS and Hi-C data, we compared eQTL data and found that the rs12411216 SNP was located near the core binding domain of transcription factor E2F4, which was the regulatory element of the GBA gene (the corresponding DHS).
Additionally, according to a large number of brain tissue sample expression datasets, mRNAs harbouring this SNP were found to be highly expressed, suggesting that this SNP might be a brain-tissue-specific regulatory marker. Thus, we hypothesized that this SNP mutation might affect the binding of transcription factors to both enhancers and promoters by changing the binding efficiency, affecting the expression of the GBA gene as a result.
Next, we recruited 306 Parkinson’s patients (including PD-NC & PD-MCI based on MMSE and MOCA scores) and determined their cognitive impairment status. After genotyping, we found that the rs12411216 SNP was marginally significantly associated with PD-MCI (P = 0.0478). Subsequently, by using both in vivo and in vitro assays, we demonstrated that the rs12411216 SNP could decrease GBA expression, enhance abnormal aggregation of α-Syn, and affect the binding efficiency of transcription factor E2F4.
Our results of 306 PD patients showed that the rs12411216 SNP was significantly associated with PD-MCI. However, the P-value (P = 0.0478) for such an association is near the borderline, which was likely due to the small sample size in this study. In another independent study in which 790 PD patients were genotyped for the association of GBA variants with PD, PD patients with GBA mutations had a higher frequency of cognitive dysfunction than non-carriers . Considering these findings, we concluded that the rs12411216 SNP is a functional variant that can serve as a genetic marker for PD-MCI.
Interestingly, our quantitative expression analysis revealed that PD patients (n = 3) with homozygous mutations (AA genotype) had lower GBA expression than patients with the other two genotypes (i.e., CA or CC). Of these three patients, one underwent cognitive impairment assessment (mainly using the MMSE and MoCA) and was identified to have PD-MCI. In addition, more studies demonstrated that GBA mutations were significantly associated with PD-MCI [3, 22], indicating that PD patients with the AA genotype for the rs12411216 SNP were more likely to develop PD-MCI.
Moreover, a recent GWAS study found a significant association between the rs12411216 SNP and occipital lobe volume (PEuropean ancestry-only = 3.9 × 10− 8) . The occipital lobe is part of the cerebral cortex and is related to cognitive ability. Using neuroimaging tests, some studies have reported that PD-MCI patients may develop local brain atrophy (including occipital atrophy) at an early stage [24, 25]. These studies suggested that occipital lobe volume might be a factor in the pathogenesis of PD-MCI.
Studies focusing on GBA in PD have pointed out that several processes might be involved in GBA-related Parkinson’s disease, including lysosomal dysfunction , α-Syn-related mechanisms , sphingolipid dysregulation , autophagy defects and protein transport defects . However, identifying specific pathways and mechanisms requires further investigation. Therefore, the goal of our future studies is to investigate the potential mechanism of GBA in PD. An increasing number of reports confirmed that most disease-related mutations were located in non-coding regulatory regions . Genome editing in human neural progenitors suggested that the distal SNP rs1191151 (downstream 700 kb) regulated FOXG1 expression, supporting its potential role as a schizophrenia risk gene . Our results showed that the DHS of rs12411216 regulates GBA expression despite being located 40 kb upstream of the GBA gene.
Genome editing and EMSA results showed that the binding activity of E2F4 was weakened, as proposed in our hypothesis. In SH-SY5Y cells, the rs12411216 mutation may not affect the expression of SNCA, but it may cause pathological aggregation of α-syn. The rs12411216 SNP is a causal SNP that affects GBA expression, and GBA activity was altered in the edited SH-SY5Y cells. GBA deficiency caused by decreased GBA enzymatic activity primarily leads to pathological accumulation of glycosphingolipids in lysosomes. Therefore, we hypothesized that the rs12411216 mutation would lead to a decrease in glucocerebrosidase (GBA protein, GCase) and enzyme activity, after which abnormally expressed GBA would affect the transport and degradation of a-Syn, thereby accelerating the formation of Lewy bodies and promoting the pathogenesis of PD (Fig. 5).
However, the current results were only based on knockout and replacement experiments on cell lines using CRISPR/Cas9-mediated genome editing methods. Evidence for the molecular genetics of PD-MCI is not currently available. A better explanation for the relationship between the rs12411216 mutation and GBA expression might be available if the brain tissue of an animal model for GCase enzyme activity could be determined. In our future research, we intend to construct an iPSC model (PD patient source) or a Parkinson’s disease mouse model involving an A53T and GBA rs12411216 double mutation to explore the relationship between rs12411216 and PD-MCI.
In summary, we showed that the distal regulatory SNP rs12411216 could affect the expression and enzymatic activity of the GBA gene and enhance α-Syn aggregation, indicating its strong relevance to the pathological progression of Parkinson’s disease. Overall, this study illustrated the relationship between GBA mutation and α-Syn aggregation and demonstrated that the rs12411216 SNP was a causative variant that could serve as a de novo biomarker for PD-MCI prognosis.
Availability of data and materials
Supplementary data are available online at Molecular Brain. All authors have reviewed the contents of the manuscript, approved its contents and validated the accuracy of the data.
Ascherio A, Schwarzschild MA. The epidemiology of Parkinson’s disease: risk factors and prevention. Lancet Neurol. 2016;15(12):1257–72.
Moore SF, Barker RA. Predictors of Parkinson's disease dementia: towards targeted therapies for a heterogeneous disease. Parkinsonism Relat Disord. 2014;20(Suppl 1):S104–7.
Mata IF, Leverenz JB, Weintraub D, Trojanowski JQ, Chen-Plotkin A, Van Deerlin VM, Ritz B, Rausch R, Factor SA, Wood-Siverio C, et al. GBA variants are associated with a distinct pattern of cognitive deficits in Parkinson's disease. Mov Disord. 2016;31(1):95–102.
Sardi SP, Clarke J, Viel C, Chan M, Tamsett TJ, Treleaven CM, Bu J, Sweet L, Passini MA, Dodge JC, et al. Augmenting CNS glucocerebrosidase activity as a therapeutic strategy for parkinsonism and other Gaucher-related synucleinopathies. Proc Natl Acad Sci U S A. 2013;110(9):3537–42.
Fagerberg L, Hallstrom BM, Oksvold P, Kampf C, Djureinovic D, Odeberg J, Habuka M, Tahmasebpoor S, Danielsson A, Edlund K, et al. Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics. Mol Cell Proteomics. 2014;13(2):397–406.
Sidransky E, Nalls MA, Aasly JO, Aharon-Peretz J, Annesi G, Barbosa ER, Bar-Shira A, Berg D, Bras J, Brice A, et al. Multicenter analysis of glucocerebrosidase mutations in Parkinson's disease. N Engl J Med. 2009;361(17):1651–61.
Hruska KS, LaMarca ME, Scott CR, Sidransky E. Gaucher disease: mutation and polymorphism spectrum in the glucocerebrosidase gene (GBA). Hum Mutat. 2008;29(5):567–83.
Migdalska-Richards A, Wegrzynowicz M, Rusconi R, Deangeli G, Di Monte DA, Spillantini MG, Schapira AHV. The L444P Gba1 mutation enhances alpha-synuclein induced loss of nigral dopaminergic neurons in mice. Brain. 2017;140(10):2706–21.
Rodriguez-Traver E, Diaz-Guerra E, Rodriguez C, Fernandez P, Arenas F, Arauzo-Bravo M, Orera M, Kulisevsky J, Moratalla R, Vicario C. A collection of integration-free iPSCs derived from Parkinson's disease patients carrying mutations in the GBA1 gene. Stem Cell Res. 2019;38:101482.
Aharon-Peretz J, Rosenbaum H, Gershoni-Baruch R. Mutations in the glucocerebrosidase gene and Parkinson's disease in Ashkenazi Jews. N Engl J Med. 2004;351(19):1972–7.
Nichols WC, Pankratz N, Marek DK, Pauciulo MW, Elsaesser VE, Halter CA, Rudolph A, Wojcieszek J, Pfeiffer RF, Foroud T, et al. Mutations in GBA are associated with familial Parkinson disease susceptibility and age at onset. Neurology. 2009;72(4):310–6.
Cullen V, Sardi SP, Ng J, Xu YH, Sun Y, Tomlinson JJ, Kolodziej P, Kahn I, Saftig P, Woulfe J, et al. Acid beta-glucosidase mutants linked to Gaucher disease, Parkinson disease, and Lewy body dementia alter alpha-synuclein processing. Ann Neurol. 2011;69(6):940–53.
Kim S, Yun SP, Lee S, Umanah GE, Bandaru VVR, Yin X, Rhee P, Karuppagounder SS, Kwon SH, Lee H, et al. GBA1 deficiency negatively affects physiological alpha-synuclein tetramers and related multimers. Proc Natl Acad Sci U S A. 2018;115(4):798–803.
Spinelli KJ, Taylor JK, Osterberg VR, Churchill MJ, Pollock E, Moore C, Meshul CK, Unni VK. Presynaptic alpha-synuclein aggregation in a mouse model of Parkinson's disease. J Neurosci. 2014;34(6):2037–50.
Ulusoy A, Musgrove RE, Rusconi R, Klinkenberg M, Helwig M, Schneider A, Di Monte DA. Neuron-to-neuron alpha-synuclein propagation in vivo is independent of neuronal injury. Acta Neuropathol Commun. 2015;3:13.
Barkhuizen M, Anderson DG, Grobler AF. Advances in GBA-associated Parkinson's disease--pathology, presentation and therapies. Neurochem Int. 2016;93:6–25.
Hsieh PC, Wang CC, Tsai CL, Yeh YM, Lee YS, Wu YR. POLG R964C and GBA L444P mutations in familial Parkinson's disease: case report and literature review. Brain Behav. 2019;9(5):e01281.
Gittelman RM, Hun E, Ay F, Madeoy J, Pennacchio L, Noble WS, Hawkins RD, Akey JM. Comprehensive identification and analysis of human accelerated regulatory DNA. Genome Res. 2015;25(9):1245–55.
Lu Y, Wang X, Yu H, Li J, Jiang Z, Chen B, Lu Y, Wang W, Han C, Ouyang Y, et al. Evolution and comprehensive analysis of DNaseI hypersensitive sites in regulatory regions of primate brain-related genes. Front Genet. 2019;10:152.
Yang B, Wu A, Hu Y, Tao C, Wang JM, Lu Y, Xing R. Mucin 17 inhibits the progression of human gastric cancer by limiting inflammatory responses through a MYH9-p53-RhoA regulatory feedback loop. J Exp Clin Cancer Res. 2019;38(1):283.
Neumann J, Bras J, Deas E, O'Sullivan SS, Parkkinen L, Lachmann RH, Li A, Holton J, Guerreiro R, Paudel R, et al. Glucocerebrosidase mutations in clinical and pathologically proven Parkinson's disease. Brain. 2009;132(Pt 7):1783–94.
D'Souza T, Rajkumar AP. Systematic review of genetic variants associated with cognitive impairment and depressive symptoms in Parkinson’s disease. Acta Neuropsychiatr. 2020;32(1):10–22.
van der Lee SJ, Knol MJ, Chauhan G, Satizabal CL, Smith AV, Hofer E, Bis JC, Hibar DP, Hilal S, van den Akker EB, et al. A genome-wide association study identifies genetic loci associated with specific lobar brain volumes. Commun Biol. 2020;32(1):10-22.
Gao Y, Nie K, Huang B, Mei M, Guo M, Xie S, Huang Z, Wang L, Zhao J, Zhang Y, et al. Changes of brain structure in Parkinson's disease patients with mild cognitive impairment analyzed via VBM technology. Neurosci Lett. 2017;658:121–32.
Zhang Z, Zhang C, Yao J, Chen X, Gao F, Jiang S, Chen W, Zhou J, Wang G. Protein-based amide proton transfer-weighted MR imaging of amnestic mild cognitive impairment. Neuroimage Clin. 2019;25:102153.
Li H, Ham A, Ma TC, Kuo SH, Kanter E, Kim D, Ko HS, Quan Y, Sardi SP, Li A, et al. Mitochondrial dysfunction and mitophagy defect triggered by heterozygous GBA mutations. Autophagy. 2019;15(1):113–30.
Fernandes HJ, Ryan BJ, Wade-Martins R. Commentary: Parkinson disease-linked GBA mutation effects reversed by molecular chaperones in human cell and fly models. Front Neurosci. 2016;10:578.
Gatto EM, Da Prat G, Etcheverry JL, Drelichman G, Cesarini M. Parkinsonisms and Glucocerebrosidase Deficiency: A Comprehensive Review for Molecular and Cellular Mechanism of Glucocerebrosidase Deficiency. Brain Sci. 2019;9(2):30.
Ikuno M, Yamakado H, Akiyama H, Parajuli LK, Taguchi K, Hara J, Uemura N, Hatanaka Y, Higaki K, Ohno K, et al. GBA haploinsufficiency accelerates alpha-synuclein pathology with altered lipid metabolism in a prodromal model of Parkinson’s disease. Hum Mol Genet. 2019;28(11):1894–904.
Lindholm Carlstrom E, Halvardson J, Etemadikhah M, Wetterberg L, Gustavson KH, Feuk L. Linkage and exome analysis implicate multiple genes in non-syndromic intellectual disability in a large Swedish family. BMC Med Genet. 2019;12(1):156.
Won H, de la Torre-Ubieta L, Stein JL, Parikshak NN, Huang J, Opland CK, Gandal MJ, Sutton GJ, Hormozdiari F, Lu D, et al. Chromosome conformation elucidates regulatory relationships in developing human brain. Nature. 2016;538(7626):523–7.
Mingding Li from Zhejiang University is acknowledged for his constructive suggestions for the writing of this manuscript.
This work was supported by the 2019 Major Projects from the Education Department of Guangdong Province (No. 2019KZDXM009), the Guangdong Provincial Science and Technology Plan Project (No. 2019A030317013), the Guangdong Provincial Natural Science Foundation – Project (No. 2019A1515010703), the High-level Hospital Construction Project (No. DFJH201907), the Key Program of Natural Science Foundation of Guangdong Province, China (No. 2017B030311015), the Guangzhou Municipal People’s Livelihood Science and Technology Project (No. 201803010085), the Tianjin Science and Technology Project (No. 16ZXMJSY00010), and the Undergraduate Innovation Fund of South China University of Technology.
Ethics approval and consent to participate
The study was approved by the Ethics Committee of Guangdong Provincial People's Hospital. All subjects or authorized representatives signed informed consent forms prior to the study procedures. The study was conducted in compliance with the Declaration of Helsinki and the International Conference on Harmonization Good Clinical Practice Guideline (World Medical Association, 1997). Written informed consent for analysis of blood sample data was obtained either from patients themselves or from their next of kin for all patients.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Jiang, Z., Huang, Y., Zhang, P. et al. Characterization of a pathogenic variant in GBA for Parkinson’s disease with mild cognitive impairment patients. Mol Brain 13, 102 (2020). https://doi.org/10.1186/s13041-020-00637-x
- Parkinson’s disease-mild cognitive impairment