Abstract
Biallelic variants in NPC1, a lysosomal gene coding for a transmembrane protein involved in cholesterol trafficking, may cause Niemann-Pick disease type C (NPC). A few cases of NPC1 mutation carriers have been reported with a Parkinson’s disease (PD) presentation. In addition, pathological studies demonstrated phosphorylated alpha-synuclein and Lewy pathology in brains of NPC patients. Therefore, we aimed to examine whether NPC1 genetic variants may be associated with PD. Full sequencing of NPC1 was performed in 2,657 PD patients and 3,647 controls from three cohorts, using targeted sequencing with molecular inversion probes. A total of 9 common variants and 126 rare variants were identified across the three cohorts. To examine association with PD, regression models adjusted for age, sex and origin were performed for common variants, and optimal sequence Kernel association test (SKAT-O) was performed for rare variants. After correction for multiple comparisons, common and rare NPC1 variants were not associated with PD. Our results do not support a link between heterozygous variants in NPC1 and PD.
1. Introduction
Parkinson’s Disease (PD) is a common neurodegenerative movement disorder, caused by a combination of genetic and environmental factors, with the normal process of aging (Ball et al., 2019; Blauwendraat et al.). The heritability of PD, based on common genetic variance, was calculated in the range of 26-36% (Nalls et al., 2019), yet rare genetic variants also contribute to PD. Numerous genes with common or rare variants associated with PD have a role in the autophagy-lysosome pathway (ALP), including GBA, LRRK2, SNCA, SMPD1, TMEM175, SCARB2 and others (Senkevich and Gan-Or, 2019).
Variants in several genes that are associated with lysosomal storage disorders have been implicated in PD. Homozygous or compound heterozygous (biallelic) mutations in the glucocerebrosidase (GBA) gene result in Gaucher disease, and heterozygous carriers of GBA variants are at increased risk for PD (Sidransky et al., 2009). Similarly, biallelic mutations in the lysosomal enzyme gene SMPD1 may cause Niemann-Pick disease types A and B, and heterozygous variants are associated with risk for PD (Alcalay et al., 2019; Gan-Or et al., 2013). Genetic variants in other genes involved in lysosomal storage disorders, such as ASAH1, GALC, GLA and others, may contribute to PD risk as well (Alcalay et al., 2018; Nalls et al., 2019; Robak et al., 2017).
NPC1 is a lysosomal gene that codes a transmembrane protein involved in cholesterol trafficking (Cruz et al., 2000). Biallelic mutations in NPC1 may cause Niemann-Pick disease type C, a lysosomal storage disorder with multi-organ involvement (Chiba et al., 2014). Several reported PD cases carrying a mutation in NPC1 have suggested that NPC1 mutations could be a risk factor for developing PD (Josephs et al., 2004; Kluenemann et al., 2013; Schneider et al., 2019). The presence of Synuclein in Lewy bodies, the pathological hallmark of PD, have also been documented in post-mortem brain studies of several Niemann-Pick type C patients (Chiba et al., 2014; Saito et al., 2004), further supporting a potential association between NPC1 mutations and PD.
To study the association between heterozygous variants in the NPC1 gene and risk for PD, we fully sequenced NPC1 in a total of 2,657 PD patients and 3,647 controls, and examined the association of both common and rare variants with PD.
2. Methods
2.1. Population
The study population included 2,657 patients with PD and 3,647 controls from three cohorts, collected at McGill University (Quebec, Canada and Montpellier, France), Columbia University (New York, NY), and Sheba Medical Center (Israel). The McGill cohort includes 1,026 PD patients and 2,588 controls, the Columbia cohort includes 1,026 PD patients and 525 controls, and the Sheba cohort includes 605 PD patients and 534 controls (detailed in Table 1). In all cohorts, to address the age and sex differences, we performed adjustment for age and sex in the statistical analysis. Participants in the McGill cohort were recruited in Québec, Canada (partially through the Quebec Parkinson Network, QPN (Gan-Or et al., 2020)) and in France, and ethnicity were confirmed using principal component analysis. The Columbia cohort has been described previously (Alcalay et al., 2016), and includes patients and controls of mixed ethnicity (European, Ashkenazi [AJ] descent and a minority of Hispanics and blacks). We adjusted for ethnicity when analyzing this cohort. The Sheba cohort was recruited in Israel and includes patients and controls of full AJ origin (all four grandparents are AJ). All three cohorts were sequenced in the same lab (McGill University), following the same protocol. All PD patients were diagnosed by movement disorder specialists according to the UK brain bank criteria (Hughes et al., 1992) or the MDS clinical diagnostic criteria (Postuma et al., 2015)..
2.2. Standard Protocol Approvals, Registrations, and Patient Consents
The institutional review board (McGill University Health Center Research Ethics Board - MUHC REB) approved the study protocols (reference number IRB00010120). All local IRBs approved the protocols and informed consent was obtained from all individual participants before entering the study.
2.3. Targeted next-generation sequencing by molecular inversion probes
The entire coding sequence of the NPC1 gene, including exon-intron boundaries (±50bps) and the 5’ and 3’ untranslated regions (UTRs), was targeted using molecular inversion probes (MIPs) as described earlier (O’Roak et al., 2012). All MIPs used to sequence NPC1 in the present study are included in Supplementary Table 1. Targeted DNA capture and amplification was done as previously described (Ross et al., 2016), and the full protocol is available upon request. The library was sequenced using Illumina HiSeq 2500/4000 platform at the McGill University and Genome Quebec Innovation Centre. Reads were mapped to the human reference genome (hg19) with Burrows-Wheeler Aligner (Li and Durbin, 2009). Genome Analysis Toolkit (GATK, v3.8) was used for post-alignment quality control and variant calling (McKenna et al., 2010), and ANNOVAR was used for annotation (Khan et al., 2019). Data on the frequency of each NPC1 variant were also extracted from the public database Genome Aggregation Database (GnomAD) (Lek et al., 2016).
2.4. Quality control
To only include high-quality variants and samples, we performed quality control (QC) by filtering out variants and samples with reduced quality, using the PLINK software v1.9 (Purcell et al., 2007). Samples with average genotyping rate of less than 87% for minimal coverage 15x and 85% for minimal coverage 50x were excluded. SNPs with genotyping rate lower than 90% were excluded. Common variants that deviated from the Hardy-Weinberg equilibrium set at p=0.001 threshold in the control populations were filtered out. Threshold for missingness difference between cases and controls was set at Bonferroni corrected p=0.05. To be included in the analysis, the minimum quality score (GQ) of genetic variants was set to 30. Rare variants (minor allele frequency, MAF<0.01) had to have a minimal coverage of >50x to be included, and common variants had to have a minimal coverage of >15x to be included in the analysis. The composition of the three cohorts following these QC steps is described in Table 1.
2.5. Statistical Analysis
The association between common NPC1 variants and PD was examined by using logistic regression models using PLINK v1.9, with the status (patient or control) as a dependent variable, age and sex as covariates in all cohorts, and AJ ancestry as an additional covariate in the NY cohort. To analyze rare variants (MAF < 0.01), an optimized sequence Kernel association test (SKAT-O, R package) (Lee et al., 2012) was performed. SKAT-O was performed for the whole gene, coding and non-coding variants and by comparing synonymous, nonsynonymous, stop frameshift and splicing variants. SKAT-O was also done separately on rare variants with Combined Annotation Dependent Depletion (CADD) score of ≥ 12.37 (Amendola et al., 2015) representing the top 2% of potentially deleterious variants.
3. Results
The average coverage of the NPC1 gene was 880X, with 100% of nucleotides covered at >15X, and 95.7% at >50X in the McGill cohort, 832X, with 100% at >15X and 95,7% at >50X in the Columbia cohort and 1108X, with 100% at >15X and 95.7% at 50X in the Sheba cohort. Supplementary Table 2 details all common variants found in the three cohorts, and their association with risk of PD. None of the common variants was associated with PD in all 3 cohorts after Bonferroni correction (Supplementary Table 2).
In the McGill cohort, 74 rare variants were included in the analysis, including 26 nonsynonymous variants, 4 splice and 3 loss-of-function variants. In the Columbia cohort, 52 rare variants were included in the analysis, including 16 nonsynonymous variants and 1 splice variant. In the Sheba cohort, 27 rare variants were included in the analysis, including 7 nonsynonymous, 1 splice and 4 loss-of-function variants. In the McGill cohort we found 3 variants that are known as pathogenic for NPC, and 4 variants with conflicting data but likely pathogenic. All these variants except one intronic splice variant were found only in controls. In the Sheba cohort we found 1 likely pathogenic variant for NPC in a control subject. All nonsynonymous, loss-of-function and rare variants affecting splice sites from the three cohorts are detailed in Supplementary Table 3. SKAT-O revealed a nominal association between all rare variants and PD (p=0.04) in the McGill cohort, which was mainly driven by synonymous intronic variants. This association was not statistically significant after correction for multiple comparisons. When including only nonsynonymous, splice and loss-of-function variants, or variants with high CADD score, no association was found in any of the cohorts (Supplementary Table 4).
4. Discussion
In the current study, we did not find any associations between common and rare variants in NPC1 and risk of PD, suggesting that NPC1 does not have a major role in PD in our cohorts. Other lysosomal genes involved in lysosomal storage disorders and in PD, such as SMPD1, ASAH1, GALC and potentially GLA, are all directly involved in the GBA glycosphingolipid metabolism pathway within the lysosome (Senkevich and Gan-Or, 2019). The lack of association with PD of NPC1 and other genes outside of this pathway (Robak et al., 2017) may suggest that PD may be specifically associated with the GBA glycosphingolipid metabolism pathway. NPC1 is involved in cholesterol trafficking and metabolism; thus far genes from this pathway have not been implicated in PD (Senkevich and Gan-Or, 2019). In our study, all but one of the NPC1 pathogenic mutations known to cause NPC or potentially pathogenic variants were identified only in controls. Nevertheless, we cannot rule out that very rare disease-causing mutations in NPC1, may be associated with PD. Much larger studies will be required to study the role of such mutations in PD.
A recent review (Schneider et al., 2019) summarized cases of patients with neurodegenerative diseases including PD who carried NPC1 variants. We found in our cohorts a total of 66 carriers of p.Asn222Ser (rs55680026) and p.Ser1004Leu (rs150334966), that have been previously suggested to be associated with PD. Both variants are not very rare and were either equally common or more common in our controls (Table 2). This likely rules out a role for these variants in PD and suggests that their presence in previously reported PD patients may represent chance findings.
Our study has several limitations. In our cohorts, differences between PD patients and controls in sex and age are significant. To address this limitation, we adjusted the regression model with age and sex as covariates. In addition, the cohorts studied here consist of relatively homogeneous populations such as Ashkenazi Jews and French Canadians. Thus, additional studies in other populations are required.
To conclude, our study does not support an important role for common and rare NPC1 variants in PD. Larger studies will be required to determine whether very rare, specific variants that are pathogenic for Niemann-Pick type C may have a role in PD.
Data Availability
All data in the manuscript is available upon request
6. Conflict of interests
Ziv Gan-Or has received consulting fees from Lysosomal Therapeutics Inc., Idorsia, Prevail Therapeutics, Denali, Ono Therapeutics, Deerfield and Inception Sciences (now Ventus). Alberto Espay has received grant support from the NIH and the Michael J Fox Foundation; personal compensation as a consultant/scientific advisory board member for Abbvie, Neuroderm, Neurocrine, Amneal, Adamas, Acadia, Acorda, InTrance, Sunovion, Lundbeck, and USWorldMeds; publishing royalties from Lippincott Williams & Wilkins, Cambridge University Press, and Springer; and honoraria from USWorldMeds, Acadia, and Sunovion. Dr. Alcalay received consultation fees for Sanofi, Roche, Janssen and Restorbio. None of these companies were involved in any parts of preparing, drafting and publishing this study.
5. Acknowledgements
We thank the participants for contributing to the study. This work was financially supported by grants from the Michael J. Fox Foundation, the Canadian Consortium on Neurodegeneration in Aging (CCNA), the Canada First Research Excellence Fund (CFREF), awarded to McGill University for the Healthy Brains for Healthy Lives initiative (HBHL), and Parkinson Canada. The Columbia University cohort is supported by the Parkinson’s Foundation, the National Institutes of Health (K02NS080915, and UL1 TR000040) and the Brookdale Foundation. KS is supported by a post-doctoral fellowship from the Canada First Research Excellence Fund (CFREF), awarded to McGill University for the Healthy Brains for Healthy Lives initiative (HBHL). GAR holds a Canada Research Chair in Genetics of the Nervous System and the Wilder Penfield Chair in Neurosciences. EAF is supported by a Foundation Grant from the Canadian Institutes of Health Research (FDN grant – 154301). ZGO is supported by the Fonds de recherche du Québec - Santé (FRQS) Chercheurs-boursiers award, in collaboration with Parkinson Quebec, and by the Young Investigator Award by Parkinson Canada. OM holds the Canada Research Chair in non-motor symptoms of Parkinson’s disease and the Tourmaline Chair in Parkinson’s disease. The access to part of the participants for this research has been made possible thanks to the Quebec Parkinson’s Network (http://rpq-qpn.ca/en/). We thank Armaghan Alam, Daniel Rochefort, Helene Catoire, Clotilde Degroot and Vessela Zaharieva for their assistance.