Figures
Abstract
Objective
Variants in GBA are associated with Lewy Body (LB) pathology. We investigated whether variants in other lysosomal storage disorder (LSD) genes also contribute to disease pathogenesis.
Methods
We performed a genetic analysis of four LSD genes including GBA, HEXA, SMPD1, and MCOLN1 in 231 brain autopsies. Brain autopsies included neuropathologically defined LBD without Alzheimer Disease (AD) changes (n = 59), AD without significant LB pathology (n = 71), Alzheimer disease and lewy body variant (ADLBV) (n = 68), and control brains without LB or AD neuropathology (n = 33). Sequencing of HEXA, SMPD1, MCOLN1 and GBA followed by ‘gene wise’ genetic association analysis was performed. To determine the functional effect, a biochemical analysis of GBA in a subset of brains was also performed. GCase activity was measured in a subset of brain samples (n = 64) that included LBD brains, with or without GBA mutations, and control brains. A lipidomic analysis was also performed in brain autopsies (n = 67) which included LBD (n = 34), ADLBV (n = 3), AD (n = 4), PD (n = 9) and control brains (n = 17), comparing GBA mutation carriers to non-carriers.
Results
In a ‘gene-wise’ analysis, variants in GBA, SMPD1 and MCOLN1 were significantly associated with LB pathology (p range: 0.03–4.14 x10-5). Overall, the mean levels of GCase activity were significantly lower in GBA mutation carriers compared to non-carriers (p<0.001). A significant increase and accumulation of several species for the lipid classes, ceramides and sphingolipids, was observed in LBD brains carrying GBA mutations compared to controls (p range: p<0.05-p<0.01).
Citation: Clark LN, Chan R, Cheng R, Liu X, Park N, Parmalee N, et al. (2015) Gene-Wise Association of Variants in Four Lysosomal Storage Disorder Genes in Neuropathologically Confirmed Lewy Body Disease. PLoS ONE 10(5): e0125204. https://doi.org/10.1371/journal.pone.0125204
Academic Editor: Coro Paisan-Ruiz, Icahn School of Medicine at Mount Sinai, UNITED STATES
Received: November 13, 2014; Accepted: March 11, 2015; Published: May 1, 2015
Copyright: © 2015 Clark et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: All deidentified genotype data and related meta data underlying the findings reported in this manuscript are available at the public repository Dryad (datadryad.org). The doi is: doi:10.5061/dryad.61c8t. Genotype data from 128 Ashkenazi Jewish (AJ) healthy individuals that were used in a secondary analysis to supplement the limited number of brain autopsy controls is available upon request at the European Genome-Phenome Archive (dataset accession: EGAD00001000781).
Funding: National Institutes of Health NS050487 (LNC) NS060113 (LNC), P50AG08702 (MS), www.nih.gov; The Parkinson's Disease Foundation (www.pdf.org). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Lewy body disorders (LBD) which include Parkinson’s Disease (PD) and Dementia with Lewy bodies (DLB) are characterized by neuronal loss in the substantia nigra (SN) and the presence of neuronal cytoplasmic inclusions composed predominantly of α-synuclein termed Lewy Bodies (LBs)[1–3]. α-synuclein immunoreactivity, including LB, have been described as features seen in the neuropathology of several lysosomal storage disorders including notably Gaucher disease (GD), but also Sandhoff disease, Tay Sachs disease, and Sanfilippo syndrome[4–8]. Heterozygosity for mutations in the gene encoding glucocerebrosidase (GBA), which cause Gaucher disease (GD), has been identified as a risk factor for both PD and DLB. We and others have shown that in sporadic and familial PD, GBA mutations are associated with early-onset PD and may modify age at onset of PD[9,10] and that in brain autopsies GBA mutation status was significantly associated with the presence of cortical LB (OR = 6.48, 95% CI, 2.45–17.16, p<0.001) and a neuropathological diagnosis of DLB after adjusting for sex, age at death, and presence of APOE-4[11]. A recent study that assessed the association of specific founder mutations in each of the lysosomal storage disorder genes HEXA, SMPD1 and MCOLN1, in 938 Ashkenazi Jewish (AJ) PD patients and 282 matched AJ controls, reported SMPD1 L302P as a risk factor for PD in the AJ population[12].
To determine whether variants in other lysosomal storage disease genes, in the same pathway as GBA, are associated with LBs we conducted an independent genetic study of the lysosomal storage disorder genes GBA, HEXA, SMPD1, and MCOLN1 in 231 brain autopsies from the New York Brain Bank at Columbia University. Brain autopsies included neuropathologically defined LBD without AD changes (n = 59), AD without significant LB pathology (n = 71), ADLBV (n = 68), and control brains without LB or AD neuropathology (n = 33). The functional effect of GBA mutations was also determined by performing a biochemical analysis of GBA in a subset of brains.
Materials and Methods
Clinical material
Brain tissue samples were obtained from the New York Brain Bank at Columbia University including cases obtained through the Alzheimer’s Disease Research Center and the Center for Parkinson’s Disease and Other Movement Disorders. Brain autopsies included neuropathologically defined LBD without AD changes (n = 59), AD without significant LB pathology (n = 71), ADLBV (n = 68), and control brains without LB or AD neuropathology (n = 33). (Table 1 and S1 Table). LB and Alzheimer plaque and tangle pathology was assessed according to published guidelines as described previously[13–15] (see S1 Methods for a detailed description of neuropathological evaluation). Clinical information on dementia was available for 208 brain autopsy samples (see S1 Methods). While the primary analysis in this paper is of autopsy proven cases, an additional group of living controls (128 Ashkenazi Jewish (AJ) healthy individuals were used in a secondary analysis to supplement the limited number of brain autopsy controls (see S1 Methods for a description of the controls). Columbia University Institutional Review Board approved the protocols and consent procedures. Written informed consent was obtained from all participants in the study.
Population Stratification and Ashkenazi Jewish Ancestry
Since a founder effect for LSD gene mutations have been reported in the AJ population we also determined AJ ancestry in brain autopsy samples. Information about AJ ancestry was not available for brain autopsies. We used two methods to examine AJ ancestry and underlying population structure in brain autopsies. In the first method, Multidimensional scaling (MDS) as implemented in the program PLINK (Version 1.07) for detecting population outliers and adjusting for population stratification was used. Briefly, we used 288, 963 autosomal SNPs for brain autopsies (n = 62), augmented with 252 AJ samples with subjects from the HapMap website (http:www.hapmap.org/), which included 90 CEU, 90 Yorubans and 90 Asians. The best fitting model assumed two underlying populations with overlap of 27 white brain autopsies with the AJ cluster and the remainder of the white brain autopsies with the white CEU cluster. In the second method, principle component analysis (PCA) as implemented in the GCTA package [16] was used to examine ancestry and admixture in white brain autopsies, AJ samples together with subjects from HapMap. Projection of all the sample genotypes along the two principle components (PC2 and PC3) is shown in S1 Methods. As in the MDS analysis performed in PLINK, there is tight clustering of 27 brain autopsies with AJ sample cluster and the remainder of the white brain autopsies cluster with CEU samples.
Molecular Genetic Analysis
Frozen cerebellar tissue was used to extract DNA. Sequencing of all GBA exons was performed as described previously[9]. Sequencing of all exons of HEXA, SMPD1 and MCOLN1 was also performed. Details of PCR and sequencing primers are available upon request. APOE genotyping was performed by MALDI-TOF mass spectrometry on the Sequenom platform as described previously[11].
Analysis of functional effect of variants
The National Center for Biotechnology information (NCBI), ClinVar, the NHLBI Exome Sequencing project (ESP) exome variant server in addition to in silico prediction was used to assess the deleterious effect of variants.
Enzyme Activity Measurements
For enzyme activity measurements, a subset (n = 64) of the total autopsy sample (n = 231) for which frozen brain tissue was available was selected based on neuropathological diagnosis and GBA mutation carrier status (GBA mutation carriers (n = 16), LBD brains without GBA mutations (n = 18) and control brains (n = 30)). Brain autopsy tissue (Cerebellum, BA4 and BA9 and ScxV) samples were homogenized in water (10% wt./vol.) using a Misonix Sonic Dismembrator and centrifuge at 30,000 Xg for 20 min. Protein concentration was determined using the Lowry method. The reaction mixture for β -glucocerebrosidase determination consisted of 50ug of protein, 50ul of 20mM 4-methylumbelliferyl-β-D-glucopyranoside, 10ul of 1M Citrate-Phosphate pH 5.0 and 10ul of 2% Sodium Tauro Deoxycholate. The reaction mixture was incubated at 37°C for 2Hours and then subsequently stopped with 2 ml 0.2 M glycine buffer, pH 10.3. The Hexosamindase A enzymatic reaction mixture consisted of 10ug of protein and 100ul of 3mM 4-methylumbelliferyl-2-acetoamido-2-deoxy-b-D-glucopyranoside in Citrate-Phosphate buffer pH4.0. Samples were incubated at 37°C for 10min and 0.2M glycine buffer was also used to stop the reaction. Fluorescence was determined in fluorescence spectrophotometer (Hitachi F-2500) at an excitation wavelength of 365 nm and emission wavelength of 448 nm. Samples were compared against a 4-methylumbelliferone (4-MU) standard curve prepared in 0.2 M glycine buffer. Enzyme activities were calculated in nmoles of 4-MU hydrolyzed/mg protein/hr. LBD brains did not carry variants in any of the other LSD genes analyzed. Frozen post-mortem interval (PMI) was available for all autopsy tissue and PMI did not appear to influence GCase activity.
Lipid Profiling
For lipid profiling, a convenience subset (n = 67) of the total autopsy sample (n = 231) was selected based on neuropathological diagnosis and GBA mutation carrier status that included LBD brains from GBA mutation carriers (n = 13), LBD brains without GBA mutations (n = 33), AD brains (n = 4) and control brains (n = 17). Characteristics of the autopsy subjects are provided in S3 Table. Lipid extracts were prepared using a modified Bligh/Dyer extraction procedure, spiked with appropriate internal standards. The samples were analysed using an Agilent 1260 HPLC system coupled to an Agilent 6490 Triple Quadrupole mass spectrometer. The lipidomic profiles generated for each sample were obtained through a combination of HPLC separation and mass spectrometry in multiple reactions monitoring mode which allows for the unambiguous identification of lipids as described previously [17,18]. (S1 Methods).
Statistical Analysis
T tests and chi square tests were used to compare continuous and categorical variables respectively. To determine whether a gene, represented by multiple sequenced variants, is associated with affection status or not, we applied the sequence kernal association test (SKAT) algorithm [19]. As above, age and sex were included in one model as covariates, and permutation based p-value was computed. To determine whether multiple variants in the lysosomal disease genes are associated with LB pathology in an additive manner, after correcting for age and sex as covariates, we applied multiple logistic regression. Although this additive model is simplistic, we reasoned that this is one way to gain insight into a set of functional (i.e., nonsynonymous) variants in the common disease pathway.
For lipidomics data, Statistical analysis for the AD and LBD mutation samples was based on the one way analysis of variance followed by post hoc Fisher’s least significant difference test while the LBD wild type samples was based on Student’s T-test. In all cases, *, p < 0.05; **, p <0.01; ***, p <0.001.
Results
Demographic and Neuropathological Characteristics of Autopsy Samples
The basic demographic and neuropathologic information of autopsy samples analysed is shown in Table 1 (all autopsies, N = 231) and S1 Table (white non-Hispanic ethnicity only, N = 196). Overall, the proportion of men (71.2%) was higher in the LBD and Alzheimer disease and lewy body variant (ADLBV) group, compared with that in the AD group (36.6%). Overall, LBD patients had a significantly earlier age at onset of dementia (66.58±10.28 years vs. 70.71±8.46; p = 0.05), earlier age at death (78.26±8.66 years vs. 81.64±8.60; p = 0.04), and had more years of education (16.61±2.33 years vs. 13.84±4.34; p = 0.02) compared to patients with AD. APOE4 allele frequencies did not differ from reported population frequencies in non-AD groups.
Sequencing and Association Analysis: Variants identified and predicted impact on function
Overall, we identified 51 (22.1%) subjects with a GBA variant, 39 (16.9%) with a HEXA variant, 33 (14.3%) with an SMPD1 variant and 62 (26.8%) with an MCOLN1 variant (Table 1). Many of the variants that we identified have been reported previously as pathogenic mutations in patients with the associated lysosomal storage disorder (Table 2). LSD variants that were significantly associated in brain autopsies with a neuropathological diagnosis of LBD are shown in Table 2.
Single Gene Wise Association: Multiple variants in GBA, SMPD1 and MCOLN1 are associated with a neuropathological diagnosis of LBD
SNP-set (Sequence) Kernal Association Test (SKAT) was used to evaluate association of variants in GBA, HEXA, SMPD1 and MCOLN1 (Table 3). When evaluating all variants, strongest association was observed for GBA variants in LBD (p = 2.95 x10-5) and ADLBV (p = 3.59 x10-2) (Table 3). Risk variants in GBA, SMPD1 and MCOLN1 were also significantly associated with LBD (p range = 0.03–4.14x10-5) and ADLBV (p range = 0.02–0.01) pathology but not AD (Table 3). We also evaluated association of protective variants and observed association of variants in SMPD1 in LBD (p = 0.03) and ADLBV (p = 0.02), but not AD, and MCOLN1 variants in LBD (p = 0.02), ADLBV (p = 0.005) but not AD (Table 3).
The following secondary analyses of the same SKAT models was also performed: 1) SKAT analysis of LSD variants with MAF<5% in all samples (n = 231)(Table 4), 2) SKAT analysis of LSD variants with MAF<5% in ‘white’ subjects only (n = 196) (Table 5) and 3) SKAT analysis of LSD variants in all samples (n = 231) using a larger control group which included the brain autopsy controls (n = 33) and the AJ controls (n = 128) (S2 Table). When we restricted the analysis to variants with MAF<5% (Table 4), strongest association was observed for GBA variants in LBD (p = 1.37x10-4). Risk variants in GBA, SMPD1 and MCOLN1 remained significantly associated with LBD (p range = 0.04–1.77x10-4) and ADLBV (p range = 0.04–0.02) pathology but not AD (Table 4). When we restricted the analysis to whites only with variants with MAF<5% (Table 5) strongest association was observed for GBA variants in LBD (p = 0.0118). Risk variants in SMPD1 in ADLBV were also significant (p = 0.0274). Although not significant, there was a trend towards significance for MCOLN1 risk variants in LBD (p = 0.189) and ADLBV (p = 0.072). However, the small sample size of the ‘white’ controls (n = 20) in this stratified analysis may be a confounding factor and the results should be interpreted with caution.
SKAT analysis using the larger control group replicated the findings observed using the brain autopsy controls (n = 33) alone (S2 Table).
Additive Effect of Multiple Variants in GBA, SMPD1 and MCOLN1
In exploratory analyses logistic regression analysis was also used to determine whether multiple variants in the same disease pathway are associated with disease pathology in an additive manner after adjusting for age and sex as covariates. Strong associations (p range: 0.03–3.8x10-5) were also observed for LBD, ADLBV, and AD cases with multiple variants in GBA+SMPD1 or GBA+SMPD1+MCOLN1 (Table 3).
GCase activity is decreased in LBD GBA mutation carriers compared to LBD non-carriers
To determine whether carrier GBA mutation status was associated with reduced enzymatic activity (haploinsufficiency) we assayed GCase activity in a subset of autopsy samples (n = 64). GCase activity was measured in LBD brains from GBA mutation carriers (n = 16), LBD brains without GBA mutations (n = 18) and control brains (n = 30) (Fig 1) from the following brain regions Cerebellum, BA4 and BA9 and ScxV.
A) GCase activity was significantly reduced in LBD cases carrying GBA mutations (n = 16) compared to LBD non-GBA carriers (n = 18) and controls (n = 30). Differences in activity for HEXA were not significant in any group. B) GCase was significantly reduced in LBD cases with mutations classified as ‘severe’ type (L444P, 84insGG etc.) compared to controls, and to LBD cases with ‘mild’ mutations (N370S) or variants of unknown phenotypic effect (E326K, T369M). Differences in activity for HEXA were not significant. * p<0.05, ** p<0.01, *** p<0.001.
The enzyme activity of a second lysosomal hydrolase, α-hexosaminidase was also assayed to demonstrate specificity of decreased activity of GCase. Overall, the mean levels of GCase activity (p<0.001) and the β-glucocerebrosidase: α -hexosaminidase ratio (p<0.001) were significantly lower in GBA mutation carriers compared to non-carriers (Fig 1). We also observed significant differential enzyme activity of GCase or for the β-glucocerebrosidase:α-hexosaminidase ratio in subjects carrying GBA mutations classified phenotypically (as in Gaucher disease) as ‘severe’ type (e.g. 84insGG, L444P) (p<0.01) compared to subjects carrying GBA mutations classified phenotypically as ‘mild’ type (e.g. N370S, R496H) (p<0.05) or of unknown phenotypic effect (E326K, T369M) (p<0.001) compared to controls (Fig 1). Lastly, we examined the relation between GCase activity in subjects with a clinical diagnosis of dementia compared to cases without dementia. Overall, the mean levels of GCase activity (p = 0.0021) and the β-glucocerebrosidase: α -hexosaminidase ratio (p = 0.0014) were significantly lower in cases with dementia than in controls. The pattern of association between GBA mutation status and the GCase activity was comparable to the combined samples, suggesting that those with dementia and a neuropathological diagnosis of DLB are driving the association.
GBA mutation carriers show significant differences in lipid species and accumulation of ceramide and sphingolipids
To determine the functional effect of reduced GCase activity in brains with LBs carrying GBA mutations compared to those without GBA mutations, AD, and control brains we performed a lipidomic analysis in postmortem brain tissue (n = 67) obtained from the primary motor cortex (BA4). The a priori hypothesis was that LBD GBA mutation carriers should show significant differences in ‘specific’ lipid species (substrate and product of GBA hydrolysis) and accumulation of ceramides and sphingolipids compared to those without GBA mutations, AD, and control brains. Characteristics of autopsy subjects with lipidomic analysis is provided as supporting data (S3 Table). The cold and frozen PMIs for autopsy tissue used in our analysis is shown in S4 Table.
Several lipid classes were significantly altered in brains with LBs carrying GBA mutations compared to controls (P range: p<0.05-p<0.01) (Fig 2 and Fig 3) and this remained significant after using an false discovery rate (FDR) control to correct for multiple comparisons of lipids (q<0.05-q<0.01). Major phospholipid subclasses such as phosphatidylcholine (PC) and phosphatidylethanolamine (PE) were decreased while phosphatidylserine (PS) was increased. There were also striking changes in sphingolipid composition. A small but significant decline in the most abundant sphingolipid sphingomyelin (SM) was seen in the diseased tissue but was compensated by increased levels of select dihydrosphingomyelin (dhSM) species and total ceramide (Cer) levels. While there was a trend towards an increase in accumulation of the known GCase substrate, GluCer, the difference was not statistically significant (Fig 2 and Fig 3). However, the complex glycosphingolipid that is biosynthetically upstream of GluCer, GM3, is highly enriched in these tissues. In addition, there is also a significant accumulation of galactosylceramide (GalCer) and its biosynthetic derivative sulfatides containing hydroxy fatty acid (Sulf-OH).
A) Heat map showing statistically significant changes in major lipid subclasses in LBD GBA mutation carriers compared to LBD wildtype, AD cases and controls and B) Heat map showing statistically significant changes in lipid classes in LBD GBA mutation carriers compared to LBD wildtype, AD cases and controls. The heat map columns reflect all significant lipid changes (q<0.05) in a diseased compared to control patients. The color bar represents the log2 value of the ratio of each lipid species. Statistical analysis for the AD and LBD Mutation samples was based on the one way analysis of variance followed by post hoc Fisher’s least significant difference test while the LBD (wildtype) samples was based on Student’s T-test. A false discovery rate control was used to correct for multiple comparisons.
Comparative lipid profile of post-mortem brain tissue obtained from patients diagnosed with various neurological conditions. The individual lipid subclasses of each group of patients was expressed as relative to control group levels for 2 separate sets of experiments (i.e. AD and LBD GBA mutation carrier relative to Control S1, LBD non carrier (wildtype) relative to Control S2). Statistical analysis for the AD and LBD Mutation samples was based on the one way analysis of variance followed by post hoc Fisher’s least significant difference test while the LBD non carrier (wildtype) samples was based on Student’s T-test. A false discovery rate control was used to correct for multiple comparisons. * q<0.05, ** q<0.01, *** q<0.001. PC, phosphatidylcholine; ePC, ether phosphatidylcholine; PE, phosphatidylethanolamine; pPE, plasmalogen phosphatidylethanolamine; PS, phosphatidylserine; PI, phosphatidylinositol; PA, phosphatidic acid; PG, phosphatidylglycerol; LBPA, lysobisphosphatidic acid; Cer, ceramide; SM, sphingomyelin; dhSM, dihydrosphingomyelin; GalCer, galactosylceramide; GluCer, glucosylceramide; Sulf, sulfatide; Sulf-h, hydroxylated sulfatide; GM3, monosialodihexosylganglioside
In brains with LBs without GBA mutations, similar changes were observed in major phospholipids PC, PE and PS levels but to a lesser degree compared to those with GBA mutations, and in the sphingolipids GalCer, Sulf and Sulf-OH. Interestingly LBD brains with and without GBA mutations displayed an accumulation in lysobisphosphatidic acid LBPA (also known as bis(monoacylglycero)phosphate), a lipid that is specifically enriched in the late endosome and lysosome was observed. LBPA also showed a trend for increase in LBD brains with GBA mutations. To determine the specificity of our assay and also provide a reference point for lipid changes in the LBD tissue, we also analyzed AD brains and found no significant changes. This is in contrast to a previous study [17], although different brain regions were analyzed (i.e., prefrontal cortex and entorhinal cortex).
Discussion
In the current study we performed a genetic analysis of four lysosomal storage disorder genes including GBA, HEXA, SMPD1, MCOLN1 in 231 brain autopsies from the New York Brain Bank at Columbia University. A biochemical analysis of GBA was also performed in a subset of brains. We show that in addition to prior reported variants in GBA, variants in SMPD1 and MCOLN1 are also significantly associated with LB or ADLBV pathology. Additional gene-wise analyses for variants based on the SKAT algorithm also identified independent association of variants in GBA, SMPD1 and MCOLN1 that were significantly associated with LBD and ADLBV pathologies but not AD. Strong association and an additive effect of multiple variants in GBA+SMPD1 or GBA+SMPD1+MCOLN1 were also observed across all disease phenotypes analysed.
The importance of the lysosomal pathway in CNS function and LBD and PD is highlighted by the identification of genetic risk factors or rare variants/mutations in lysosomal genes in case-control association studies (GBA and NAGLU)[7,9–11], GWAS studies (LAMP3, SCARB2)[24,25] or linkage analysis and exome sequencing in PD families (ATP13A2, VPS35 (endolysosomal pathway)[26,27]. A large multisite study of brain autopsy samples from subjects with different forms of dementia identified GBA mutations in 7.6% (6/79) of pure DLB cases (OR, 7.6 [95% CI, 1.8–31.9]) compared to 3.6% (8/222) of ADLBV cases (OR, 4.6 [95% CI, 1.2–17.6]) [28].
In our study, multiple variants predicted to be deleterious or damaging in GBA, SMPD1 and MCOLN1 in autopsy samples were significantly associated with LB and ADLBV pathology. We identified a total of 26 variants in autopsy samples that have been previously reported as mutations in lysosomal storage disorders. These causal mutations in lysosomal storage disorders are usually observed in the homozygous or compound heterozygous state whereas in the autopsy samples that we examined these mutations were observed in the heterozygous state suggesting that haploinsufficiency of lysosomal genes may contribute to LB and ADLBV phenotype. Overall, ~15% of all LBD autopsy samples also carried a variant in more than one lysosomal storage disease gene examined suggesting that ‘multiple hits’ in the same biochemical pathway, the lysosomal pathway, might increase risk for LBD.
Our data also shows that GBA mutation status is associated with significantly reduced GCase activity and a neuropathological diagnosis of LBD suggesting that haploinsufficiency or partial enzyme activity leads to increase in α -synuclein levels and Lewy body pathology. A decrease in GCase activity has been reported previously in brain autopsies from patients with Type I Gaucher Disease and parkinsonism and more recently in brain autopsies from patients with PD that carry GBA mutations [4,29]. A decrease in GCase activity has also been reported in brain autopsies from patients with sporadic PD without GBA mutations, with the greatest reduction in the substantia nigra [29] and conflicting reports of decreased GCase activity in the frontal cortical regions [29, 30]. In LBD brain autopsy samples from the frontal cortex without GBA mutations we did not observe a decrease in GCase activity and our findings are consistent with one published study [29]. These conflicting reports of a reduction in GCase activity in different brain regions from PD or LBD autopsy samples without GBA mutations may reflect different stages of disease progression, neuronal loss or α -synuclein accumulation.
We have previously demonstrated the utility of lipidomics as a means to understand dysregulation of lipid metabolism and generate novel insights linked to AD pathogenesis[17]. Applying similar methodologies to the analysis of the motor cortex region of LBD brains with and without GBA mutations, we observed that there are significant alterations in both major phospholipid and sphingolipid subclasses compared to controls. In the case of LBD carrying GBA mutations, GluCer was not significantly accumulated as one might expect. This may be reconciled by the existence of non-lysosomal glucosylceramidase GBA2 that is significantly expressed in the brain and can compensate for the deficiency in GBA activity [31, 32]. Nevertheless, it is worth noting that in both cases of LBD with and without GBA mutation, there appears to be significant or a trend towards accumulation of other sphingolipid subclasses including Cer, GalCer, Sulf, Sulf-OH and GM3 and the unusual phospholipid LBPA. LBPA is enriched in late endosomes where it functions in biogenesis of multivesicular bodies [33] and also in lysosomes where it plays a role in stimulating the hydrolysis of membrane bound sphingolipids. The overall profile of both LBD with and without GBA mutation cases suggests a common theme of dysfunction occurring in the endolysosomal degradative pathway that ultimately lead to defects in lysosomal clearance of autophagosomes and an accumulation of α -synuclein in LBD.
Recent studies in AD suggest that disease pathogenesis may begin more than 20 years before the onset of dementia [34]. Similarly in LBD, non-motor symptoms may predate motor symptoms by decades suggesting pathophysiological changes before clinical onset [34–37]. To date, there are no effective biomarkers for LBD. Our study suggests that combined genetic and lipidomic data may prove effective in disease risk prediction, biomarker development (CSF) and targeted therapeutic strategies.
Supporting Information
S1 Table. Demographic and Neuropathological Characteristics of White Autopsy Subjects.
*APOE4 missing for 16 cases
https://doi.org/10.1371/journal.pone.0125204.s002
(DOCX)
S2 Table. Gene-wise association SKAT analysis with AJ controls.
A) Gene-wise association with SKAT analysis with AJ controls only (n = 128) and B) Gene wise association SKAT analysis with brain controls (n = 33) and AJ controls (n = 128). *Corrected for covariates. ** Indicates number of markers included in the test. 1 Risk variants are variants more frequent among cases than controls; whereas, variants are considered protective when they are more frequent in controls than cases.
https://doi.org/10.1371/journal.pone.0125204.s003
(DOCX)
S3 Table. Characteristics of Autopsy Subjects who had Lipidomic Analysis.
* APOE missing in 2 cases.
https://doi.org/10.1371/journal.pone.0125204.s004
(DOCX)
S4 Table. Cold and Frozen PMIs brain autopsies.
https://doi.org/10.1371/journal.pone.0125204.s005
(DOCX)
Author Contributions
Conceived and designed the experiments: LNC R. Chan GDP JHL. Performed the experiments: R. Chan N. Parmalee SK PAT. Analyzed the data: LNC R. Cheng R. Chan XL N. Parmalee N. Park SK PAT JHL. Contributed reagents/materials/analysis tools: PAT GMP R. Chan GDP. Wrote the paper: LNC R. Chan GDP JHL. Contributed autopsy samples: EC JPV RA KM LLH RM MS. Critical review of the manuscript: LNC R. Chan R. Cheng XL N. Parmalee N. Park SK EC PAT GMP JPV RA KM LLH SF RM MS GDP JHL.
References
- 1. McKeith IG, Perry EK, Perry RH. Report of the second dementia with Lewy body international workshop: diagnosis and treatment. Consortium on Dementia with Lewy Bodies. Neurology. 1999;53: 902–905. pmid:10496243
- 2. McKeith IG. Clinical Lewy body syndromes. Ann N Y Acad Sci. 2000;920: 1–8. pmid:11193136
- 3. McKeith IG. Spectrum of Parkinson's disease, Parkinson's dementia, and Lewy body dementia. Neurol Clin. 2000;18: 865–902. pmid:11072265
- 4. Tayebi N, Walker J, Stubblefield B, Orvisky E, LaMarca ME, Wong K, et al. Gaucher disease with parkinsonian manifestations: does glucocerebrosidase deficiency contribute to a vulnerability to parkinsonism? Mol Genet Metab. 2003;79: 104–109. pmid:12809640
- 5. Suzuki K, Iseki E, Togo T, Yamaguchi A, Katsuse O, Katsuyama K, et al. Neuronal and glial accumulation of alpha- and beta-synucleins in human lipidoses. Acta Neuropathol. 2007;114: 481–489. pmid:17653558
- 6. Choi JH, Stubblefield B, Cookson MR, Goldin E, Velayati A, Tayebi N, et al. Aggregation of α-synuclein in brain samples from subjects with glucocerebrosidase mutations. Mol Genet Metab. 2011;104: 185–188. pmid:21742527
- 7. Winder-Rhodes SE, Garcia-Reitbock P, Ban M, Evans JR, Jacques TS, Kemppinen A, et al. Genetic and pathological links between Parkinson's disease and the lysosomal disorder Sanfilippo syndrome. Mov Disord. 2012;27: 312–315. pmid:22102531
- 8. Cox TM, Cachon-Gonzalez MB. The cellular pathology of lysosomal diseases. J Pathol. 2012;226: 241–254. pmid:21990005
- 9. Clark LN, Ross BM, Wang Y, Mejia-Santana H, Harris J, Louis ED, et al. Mutations in the glucocerebrosidase gene are associated with early-onset Parkinson disease. Neurology. 2007;69: 1270–1277. pmid:17875915
- 10. Sidransky E, Nalls MA, Aasly JO, Aharon-Peretz J, Annesi G, Barbosa ER, et al. Multicenter analysis of glucocerebrosidase mutations in Parkinson's disease. N Engl J Med. 2009;361: 1651–1661. pmid:19846850
- 11. Clark LN, Kartsaklis LA, Wolf Gilbert R, Dorado B, Ross BM, Kisselev S, et al. Association of glucocerebrosidase mutations with dementia with lewy bodies. Arch Neurol.2009;66: 578–583. pmid:19433657
- 12. Gan-Or Z, Ozelius LJ, Bar-Shira A, Saunders-Pullman R, Mirelman A, Kornreich R, et al. The p.L302P mutation in the lysosomal enzyme gene SMPD1 is a risk factor for Parkinson disease. Neurology. 2013;80: 1606–1610. pmid:23535491
- 13. McKeith IG, Dickson DW, Lowe J, Emre M, O'Brien JT, Feldman H, et al. Diagnosis and management of dementia with Lewy bodies: third report of the DLB Consortium. Neurology. 2005;65: 1863–1872. pmid:16237129
- 14. Mirra SS, Heyman A, McKeel D, Sumi SM, Crain BJ, Brownlee LM, et al. The Consortium to Establish a Registry for Alzheimer's Disease (CERAD). Part II. Standardization of the neuropathologic assessment of Alzheimer's disease. Neurology. 1991;41: 479–486. pmid:2011243
- 15.
NACC Neuropathology (NP) Diagnosis Coding Guidebook. Seattle: University of Washington: National Alzheimer’s Coordinating Center. 2005–2008.
- 16. Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet. 2011;88: 76–82. pmid:21167468
- 17. Chan RB, Oliveira TG, Cortes EP, Honig LS, Duff KE, Small SA, et al. Comparative lipidomic analysis of mouse and human brain with Alzheimer disease. J Biol Chem. 2012;287: 2678–2688. pmid:22134919
- 18. Chan R, Uchil PD, Jin J, Shui G, Ott DE, Mothes W, et al. Retroviruses human immunodeficiency virus and murine leukemia virus are enriched in phosphoinositides. J Virol. 2008;82: 11228–11238. pmid:18799574
- 19. Wu MC, Lee S, Cai T, Li Y, Boehnke M, Lin X. Rare-variant association testing for sequencing data with the sequence kernel association test. Am J Hum Genet. 2011;89: 82–93. pmid:21737059
- 20. Hruska KS, LaMarca ME, Scott CR, Sidransky E. Gaucher disease: mutation and polymorphism spectrum in the glucocerebrosidase gene (GBA). Hum Mutat. 2008;29: 567–583. pmid:18338393
- 21. Irun P, Mallen M, Dominguez C, Rodriguez-Sureda V, Alvarez-Sala LA, Arslan N, et al. Identification of seven novel SMPD1 mutations causing Niemann-Pick disease types A and B. Clin Genet. 2013;84: 356–361. pmid:23252888
- 22. Simonaro CM, Desnick RJ, McGovern MM, Wasserstein MP, Schuchman EH. The demographics and distribution of type B Niemann-Pick disease: novel mutations lead to new genotype/phenotype correlations. Am J Hum Genet. 2002;71: 1413–1419. pmid:12369017
- 23. Triggs-Raine B, Richard M, Wasel N, Prence EM, Natowicz MR. Mutational analyses of Tay-Sachs disease: studies on Tay-Sachs carriers of French Canadian background living in New England. Am J Hum Genet. 1995;56: 870–879. pmid:7717398
- 24. Nalls MA, Plagnol V, Hernandez DG, Sharma M, Sheerin U-M, Saad M, 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–649. pmid:21292315
- 25. Do CB, Tung JY, Dorfman E, Kiefer AK, Drabant EM, Francke U, 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(6): e1002141. pmid:21738487
- 26. Di Fonzo A, Chien HF, Socal M, Giraudo S, Tassorelli C, Iliceto G, et al. ATP13A2 missense mutations in juvenile parkinsonism and young onset Parkinson disease. Neurology. 2007;68: 1557–1562. pmid:17485642
- 27. Vilarino-Guell C, Wider C, Ross OA, Dachsel JC, Kachergus JM, Lincoln SJ, et al. VPS35 mutations in Parkinson disease. Am J Hum Genet. 2011;89: 162–167. pmid:21763482
- 28. Tsuang D, Leverenz JB, Lopez OL, Hamilton RL, Bennett DA, Schneider JA, et al. GBA mutations increase risk for Lewy body disease with and without Alzheimer disease pathology. Neurology. 2012;79: 1944–1950. pmid:23035075
- 29. Gegg ME, Burke D, Heales SJR, Cooper JM, Hardy J, Wood NW, et al. Glucocerebrosidase deficiency in substantia nigra of parkinson disease brains. Ann Neurol. 2012;72: 455–463. pmid:23034917
- 30. Murphy KE, Gysbers AM, Abbott SK, Tayebi N, Kim WS, Sidransky E, et al. Reduced glucocerebrosidase is associated with increased α-synuclein in sporadic Parkinson's disease. Brain. 2014;137: 834–848. pmid:24477431
- 31. Boot RG, Verhoek M, Donker-Koopman W, Strijland A, van Marle J, Overkleeft HS, et al. Identification of the non-lysosomal glucosylceramidase as beta-glucosidase 2. J Biol Chem. 2007;282: 1305–1312. pmid:17105727
- 32. Yildiz Y, Matern H, Thompson B, Allegood JC, Warren RL, Ramirez DMO, et al. Mutation of beta-glucosidase 2 causes glycolipid storage disease and impaired male fertility. J Clin Invest. 2006;116: 2985–2994. pmid:17080196
- 33. Kobayashi T, Beuchat MH, Lindsay M, Frias S, Palmiter RD, Sakuraba H, et al. Late endosomal membranes rich in lysobisphosphatidic acid regulate cholesterol transport. Nat Cell Biol. 1999;1: 113–118. pmid:10559883
- 34. Bateman RJ, Xiong C, Benzinger TLS, Fagan AM, Goate A, Fox NC, et al. Clinical and biomarker changes in dominantly inherited Alzheimer's disease. N Engl J Med. 2012;367: 795–804. pmid:22784036
- 35. Braak H, Del Tredici K, Rub U, de Vos RAI, Jansen Steur ENH, Braak E. Staging of brain pathology related to sporadic Parkinson's disease. Neurobiol Aging. 2003;24: 197–211. pmid:12498954
- 36. Hawkes CH. The prodromal phase of sporadic Parkinson's disease: does it exist and if so how long is it? Mov Disord. 2008;23: 1799–1807. pmid:18759359
- 37. Koller WC. When does Parkinson's disease begin? Neurology. 1992; 42: 27–31; discussion 41–28. pmid:1350072