Analysis of common and rare VPS13C variants in late-onset Parkinson disease

Association of Rare Variants in ARSA with Parkinson's Disease

BACKGROUND

Several lysosomal genes are associated with Parkinson's disease (PD), yet the association between PD and ARSA remains unclear.

OBJECTIVES

To study rare ARSA variants in PD.

METHODS

To study rare ARSA variants (minor allele frequency < 0.01) in PD, we performed burden analyses in six independent cohorts with 5801 PD patients and 20,475 controls, followed by a meta-analysis.

RESULTS

We found evidence for associations between functional ARSA variants and PD in four cohorts (P ≤ 0.05 in each) and in the meta-analysis (P = 0.042). We also found an association between loss-of-function variants and PD in the United Kingdom Biobank cohort (P = 0.005) and in the meta-analysis (P = 0.049). These results should be interpreted with caution as no association survived multiple comparisons correction. Additionally, we describe two families with potential co-segregation of ARSA p.E382K and PD.

CONCLUSIONS

Rare functional and loss-of-function ARSA variants may be associated with PD. Further replications in large case-control/familial cohorts are required. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Association of rare variants in ARSA with Parkinson's disease

Background

Several lysosomal genes are associated with Parkinson's disease (PD), yet the association between PD and ARSA , which encodes for the enzyme arylsulfatase A, remains controversial.

Objectives

To evaluate the association between rare ARSA variants and PD.

Methods

To study possible association of rare variants (minor allele frequency<0.01) in ARSA with PD, we performed burden analyses in six independent cohorts with a total of 5,801 PD patients and 20,475 controls, using optimized sequence Kernel association test (SKAT-O), followed by a meta-analysis.

Results

We found evidence for an association between functional ARSA variants and PD in four independent cohorts (P≤0.05 in each) and in the meta-analysis (P=0.042). We also found an association between loss-of-function variants and PD in the UKBB cohort (P=0.005) and in the meta-analysis (P=0.049). However, despite replicating in four independent cohorts, these results should be interpreted with caution as no association survived correction for multiple comparisons. Additionally, we describe two families with potential co-segregation of the ARSA variant p.E384K and PD.

Conclusions

Rare functional and loss-of-function ARSA variants may be associated with PD. Further replication in large case-control cohorts and in familial studies is required to confirm these associations.

ATG2 and VPS13 proteins: molecular highways transporting lipids to drive membrane expansion and organelle communication

Communication between organelles is an essential process that helps maintain cellular homeostasis and organelle contact sites have recently emerged as crucial mediators of this communication. The emergence of a class of molecular bridges that span the inter-organelle gaps has now been shown to direct the flow of lipid traffic from one lipid bilayer to another. One of the key components of these molecular bridges is the presence of an N-terminal Chorein/VPS13 domain. This is an evolutionarily conserved domain present in multiple proteins within the endocytic and autophagy trafficking pathways. Herein, we discuss the current state-of-the-art of this class of proteins, focusing on the role of these lipid transporters in the autophagy and endocytic pathways. We discuss the recent biochemical and structural advances that have highlighted the essential role Chorein-N domain containing ATG2 proteins play in driving the formation of the autophagosome and how lipids are transported from the endoplasmic reticulum to the growing phagophore. We also consider the VPS13 proteins, their role in organelle contacts and the endocytic pathway and highlight how disease-causing mutations disrupt these contact sites. Finally, we open the door to discuss other Chorein_N domain containing proteins, for instance, UHRF1BP1/1L, their role in disease and look towards prokaryote examples of Chorein_N-like domains. Taken together, recent advances have highlighted an exciting opportunity to delve deeper into inter-organelle communication and understand how lipids are transported between membrane bilayers and how this process is disrupted in multiple diseases.

Clinical Manifestations and Molecular Backgrounds of Parkinson's Disease Regarding Genes Identified From Familial and Population Studies

Over the past 20 years, numerous robust analyses have identified over 20 genes related to familial Parkinson's disease (PD), thereby uncovering its molecular underpinnings and giving rise to more sophisticated approaches to investigate its pathogenesis. α-Synuclein is a major component of Lewy bodies (LBs) and behaves in a prion-like manner. The discovery of α-Synuclein enables an in-depth understanding of the pathology behind the generation of LBs and dopaminergic neuronal loss. Understanding the pathophysiological roles of genes identified from PD families is uncovering the molecular mechanisms, such as defects in dopamine biosynthesis and metabolism, excessive oxidative stress, dysfunction of mitochondrial maintenance, and abnormalities in the autophagy–lysosome pathway, involved in PD pathogenesis. This review summarizes the current knowledge on familial PD genes detected by both single-gene analyses obeying the Mendelian inheritance and meta-analyses of genome-wide association studies (GWAS) from genome libraries of PD. Studying the functional role of these genes might potentially elucidate the pathological mechanisms underlying familial PD and sporadic PD and stimulate future investigations to decipher the common pathways between the diseases.

Rare PSAP Variants and Possible Interaction with GBA in REM Sleep Behavior Disorder

BACKGROUND

PSAP encodes saposin C, the co-activator of glucocerebrosidase, encoded by GBA. GBA mutations are associated with idiopathic/isolated REM sleep behavior disorder (iRBD), a prodromal stage of synucleinopathy.

OBJECTIVE

To examine the role of PSAP mutations in iRBD.

METHODS

We fully sequenced PSAP and performed Optimized Sequence Kernel Association Test in 1,113 iRBD patients and 2,324 controls. We identified loss-of-function (LoF) mutations, which are very rare in PSAP, in three iRBD patients and none in controls (uncorrected p = 0.018).

RESULTS

Two variants were stop mutations, p.Gln260Ter and p.Glu166Ter, and one was an in-frame deletion, p.332_333del. All three mutations have a deleterious effect on saposin C, based on in silico analysis. In addition, the two carriers of p.Glu166Ter and p.332_333del mutations also carried a GBA variant, p.Arg349Ter and p.Glu326Lys, respectively. The co-occurrence of these extremely rare PSAP LoF mutations in two (0.2%) GBA variant carriers in the iRBD cohort, is unlikely to occur by chance (estimated co-occurrence in the general population based on gnomAD data is 0.00035%). Although none of the three iRBD patients with PSAP LoF mutations have phenoconverted to an overt synucleinopathy at their last follow-up, all manifested initial signs suggestive of motor dysfunction, two were diagnosed with mild cognitive impairment and all showed prodromal clinical markers other than RBD. Their probability of prodromal PD, according to the Movement Disorder Society research criteria, was 98% or more.

CONCLUSION

These results suggest a possible role of PSAP variants in iRBD and potential genetic interaction with GBA, which requires additional studies.

Multi-modality machine learning predicting Parkinson's disease

Personalized medicine promises individualized disease prediction and treatment. The convergence of machine learning (ML) and available multimodal data is key moving forward. We build upon previous work to deliver multimodal predictions of Parkinson's disease (PD) risk and systematically develop a model using GenoML, an automated ML package, to make improved multi-omic predictions of PD, validated in an external cohort. We investigated top features, constructed hypothesis-free disease-relevant networks, and investigated drug-gene interactions. We performed automated ML on multimodal data from the Parkinson's progression marker initiative (PPMI). After selecting the best performing algorithm, all PPMI data was used to tune the selected model. The model was validated in the Parkinson's Disease Biomarker Program (PDBP) dataset. Our initial model showed an area under the curve (AUC) of 89.72% for the diagnosis of PD. The tuned model was then tested for validation on external data (PDBP, AUC 85.03%). Optimizing thresholds for classification increased the diagnosis prediction accuracy and other metrics. Finally, networks were built to identify gene communities specific to PD. Combining data modalities outperforms the single biomarker paradigm. UPSIT and PRS contributed most to the predictive power of the model, but the accuracy of these are supplemented by many smaller effect transcripts and risk SNPs. Our model is best suited to identifying large groups of individuals to monitor within a health registry or biobank to prioritize for further testing. This approach allows complex predictive models to be reproducible and accessible to the community, with the package, code, and results publicly available.

Monogenic Parkinson’s Disease: Genotype, Phenotype, Pathophysiology, and Genetic Testing

Parkinson’s disease may be caused by a single pathogenic variant (monogenic) in 5–10% of cases, but investigation of these disorders provides valuable pathophysiological insights. In this review, we discuss each genetic form with a focus on genotype, phenotype, pathophysiology, and the geographic and ethnic distribution. Well-established Parkinson’s disease genes include autosomal dominant forms (SNCA, LRRK2, and VPS35) and autosomal recessive forms (PRKN, PINK1 and DJ1). Furthermore, mutations in the GBA gene are a key risk factor for Parkinson’s disease, and there have been major developments for X-linked dystonia parkinsonism. Moreover, atypical or complex parkinsonism may be due to mutations in genes such as ATP13A2, DCTN1, DNAJC6, FBXO7, PLA2G6, and SYNJ1. Furthermore, numerous genes have recently been implicated in Parkinson’s disease, such as CHCHD2, LRP10, TMEM230, UQCRC1, and VPS13C. Additionally, we discuss the role of heterozygous mutations in autosomal recessive genes, the effect of having mutations in two Parkinson’s disease genes, the outcome of deep brain stimulation, and the role of genetic testing. We highlight that monogenic Parkinson’s disease is influenced by ethnicity and geographical differences, reinforcing the need for global efforts to pool large numbers of patients and identify novel candidate genes.

Targeted sequencing of Parkinson's disease loci genes highlights SYT11, FGF20 and other associations

Genome-wide association studies (GWAS) have identified numerous loci associated with Parkinson's disease. The specific genes and variants that drive the associations within the vast majority of these loci are unknown. We aimed to perform a comprehensive analysis of selected genes to determine the potential role of rare and common genetic variants within these loci. We fully sequenced 32 genes from 25 loci previously associated with Parkinson's disease in 2657 patients and 3647 controls from three cohorts. Capture was done using molecular inversion probes targeting the exons, exon-intron boundaries and untranslated regions (UTRs) of the genes of interest, followed by sequencing. Quality control was performed to include only high-quality variants. We examined the role of rare variants (minor allele frequency < 0.01) using optimized sequence Kernel association tests. The association of common variants was estimated using regression models adjusted for age, sex and ethnicity as required in each cohort, followed by a meta-analysis. After Bonferroni correction, we identified a burden of rare variants in SYT11, FGF20 and GCH1 associated with Parkinson's disease. Nominal associations were identified in 21 additional genes. Previous reports suggested that the SYT11 GWAS association is driven by variants in the nearby GBA gene. However, the association of SYT11 was mainly driven by a rare 3' UTR variant (rs945006601) and was independent of GBA variants (P = 5.23 × 10-5 after exclusion of all GBA variant carriers). The association of FGF20 was driven by a rare 5' UTR variant (rs1034608171) located in the promoter region. The previously reported association of GCH1 with Parkinson's disease is driven by rare non-synonymous variants, some of which are known to cause dopamine-responsive dystonia. We also identified two LRRK2 variants, p.Arg793Met and p.Gln1353Lys, in 10 and eight controls, respectively, but not in patients. We identified common variants associated with Parkinson's disease in MAPT, TMEM175, BST1, SNCA and GPNMB, which are all in strong linkage disequilibrium with known GWAS hits in their respective loci. A common coding PM20D1 variant, p.Ile149Val, was nominally associated with reduced risk of Parkinson's disease (odds ratio 0.73, 95% confidence interval 0.60-0.89, P = 1.161 × 10-3). This variant is not in linkage disequilibrium with the top GWAS hits within this locus and may represent a novel association. These results further demonstrate the importance of fine mapping of GWAS loci, and suggest that SYT11, FGF20, and potentially PM20D1, BST1 and GPNMB should be considered for future studies as possible Parkinson's disease-related genes.

LRRK2 p.M1646T is associated with glucocerebrosidase activity and with Parkinson's disease

The LRRK2 p.G2019S Parkinson's disease (PD) variant is associated with elevated glucocerebrosidase (GCase) activity in peripheral blood. We aimed to evaluate the association of other LRRK2 variants with PD and its association with GCase activity. LRRK2 and GBA were fully sequenced in 1123 PD patients and 576 controls from the Columbia and PPMI cohorts, in which GCase activity was measured in dried blood spots by liquid chromatography-tandem mass spectrometry. LRRK2 p.M1646T was associated with increased GCase activity in the Columbia University cohort (β = 1.58, p = 0.0003), and increased but not significantly in the PPMI cohort (β = 0.29, p = 0.58). p.M1646T was associated with PD (odds ratio = 1.18, 95% confidence interval = 1.09-1.28, p = 7.33E-05) in 56,306 PD patients and proxy-cases, and 1.4 million controls. Our results suggest that the p.M1646T variant is associated with risk of PD with a small effect and with increased GCase activity in peripheral blood.

Contribution of rare homozygous and compound heterozygous VPS13C missense mutations to dementia with Lewy bodies and Parkinson's disease

Dementia with Lewy bodies (DLB) and Parkinson's disease (PD) are clinically, pathologically and etiologically disorders embedded in the Lewy body disease (LBD) continuum, characterized by neuronal α-synuclein pathology. Rare homozygous and compound heterozygous premature termination codon (PTC) mutations in the Vacuolar Protein Sorting 13 homolog C gene (VPS13C) are associated with early-onset recessive PD. We observed in two siblings with early-onset age (< 45) and autopsy confirmed DLB, compound heterozygous missense mutations in VPS13C, p.Trp395Cys and p.Ala444Pro, inherited from their healthy parents in a recessive manner. In lymphoblast cells of the index patient, the missense mutations reduced VPS13C expression by 90% (p = 0.0002). Subsequent, we performed targeted resequencing of VPS13C in 844 LBD patients and 664 control persons. Using the optimized sequence kernel association test, we obtained a significant association (p = 0.0233) of rare VPS13C genetic variants (minor allele frequency ≤ 1%) with LBD. Among the LBD patients, we identified one patient with homozygous missense mutations and three with compound heterozygous missense mutations in trans position, indicative for recessive inheritance. In four patients with compound heterozygous mutations, we were unable to determine trans position. The frequency of LBD patient carriers of proven recessive compound heterozygous missense mutations is 0.59% (5/844). In autopsy brain tissue of two unrelated LBD patients, the recessive compound heterozygous missense mutations reduced VPS13C expression. Overexpressing of wild type or mutant VPS13C in HeLa or SH-SY5Y cells, demonstrated that the mutations p.Trp395Cys or p.Ala444Pro, abolish the endosomal/lysosomal localization of VPS13C. Overall, our data indicate that rare missense mutations in VPS13C are associated with LBD and recessive compound heterozygous missense mutations might have variable effects on the expression and functioning of VPS13C. We conclude that comparable to the recessive inherited PTC mutations in VPS13C, combinations of rare recessive compound heterozygous missense mutations reduce VPS13C expression and contribute to increased risk of LBD.

Analysis of Heterozygous PRKN Variants and Copy-Number Variations in Parkinson's Disease

BACKGROUND

Biallelic PRKN mutation carriers with Parkinson's disease (PD) typically have an earlier disease onset, slow disease progression, and, often, different neuropathology compared to sporadic PD patients. However, the role of heterozygous PRKN variants in the risk of PD is controversial.

OBJECTIVES

Our aim was to examine the association between heterozygous PRKN variants, including single-nucleotide variants and copy-number variations (CNVs), and PD.

METHODS

We fully sequenced PRKN in 2809 PD patients and 3629 healthy controls, including 1965 late-onset (63.97 ± 7.79 years, 63% men) and 553 early-onset PD patients (43.33 ± 6.59 years, 68% men). PRKN was sequenced using targeted next-generation sequencing with molecular inversion probes. CNVs were identified using a combination of multiplex ligation-dependent probe amplification and ExomeDepth. To examine whether rare heterozygous single-nucleotide variants and CNVs in PRKN are associated with PD risk and onset, we used optimized sequence kernel association tests and regression models.

RESULTS

We did not find any associations between all types of PRKN variants and risk of PD. Pathogenic and likely-pathogenic heterozygous single-nucleotide variants and CNVs were less common among PD patients (1.52%) than among controls (1.8%, false discovery rate-corrected P = 0.55). No associations with age at onset and in stratified analyses were found.

CONCLUSIONS

Heterozygous single-nucleotide variants and CNVs in PRKN are not associated with PD. Molecular inversion probes allow for rapid and cost-effective detection of all types of PRKN variants, which may be useful for pretrial screening and for clinical and basic science studies targeting specifically PRKN patients. © 2020 International Parkinson and Movement Disorder Society.

Mutation screening and burden analysis of VPS13C in Chinese patients with early-onset Parkinson's disease

Homozygous and compound heterozygous mutations in the vacuolar protein sorting 13C (VPS13C) gene can cause autosomal recessive parkinsonism via mitochondrial pathway. The present study aimed to screen the mutations of VPS13C in a cohort of Chinese patients with early-onset Parkinson's disease (EOPD) and further explore its pathogenicity via burden analysis. A total of 669 patients with EOPD were sequenced with whole-exome sequencing and analyzed homozygous or compound heterozygous mutations in VPS13C. Moreover, rare variants with minor allele frequency <0.1% were included in the burden analysis. In total, 7 (1.05%) patients with EOPD carried compound heterozygous mutations in VPS13C, including 3 patients with novel compound heterozygous missense mutations and 4 patients with at least 1 nonsense or splicing-site mutations. Furthermore, burden analysis indicated that patients with EOPD had an enrichment of rare variants in VPS13C. In conclusion, our findings of compound missense mutations expanded the mutation spectrum of VPS13C in EOPD. Burden analysis further elucidated the importance of VPS13C in the pathogenesis of PD.