Rare GCH1 heterozygous variants contributing to Parkinson's disease

The Role of Cellular Defense Systems of Ferroptosis in Parkinson's Disease and Alzheimer's Disease

Parkinson's disease (PD) and Alzheimer's disease (AD) are the most common rapidly developing neurodegenerative diseases that lead to serious health and socio-economic consequences. Ferroptosis is a non-apoptotic form of cell death; there is growing evidence to support the notion that ferroptosis is involved in a variety of pathophysiological contexts, and there is increasing interest in the role of ferroptosis in PD and AD. Simultaneously, cells may have evolved four defense systems to counteract the toxic effects of ferroptosis occasioned by lipid peroxidation. This review, which focuses on the analysis of ferroptosis in the PD and AD context, outlines four cellular defense systems against ferroptosis and how each of them is involved in PD and AD.

A systematic analysis of genotype-phenotype associations with PLA2G6

BACKGROUND

PLA2G6-associated neurodegeneration (PLAN) can be categorized into infantile neuroaxonal dystrophy (INAD), atypical neuroaxonal dystrophy (aNAD), neurodegeneration with brain iron accumulation (NBIA), and early-onset parkinsonism (EOP).

OBJECTIVES

To determine the genotype-phenotype association in PLAN.

METHODS

"PLA2G6" or "PARK14" or "phospholipase A2 group VI" or "iPLA2β" were searched across MEDLINE from June 23, 1997, to March 1, 2023. A total of 391 patients were identified, and 340 patients of them were finally included in the assessment.

RESULTS

The loss of function (LOF) mutation ratios were significantly different (p < 0.001), highest in INAD, followed by NBIA, aNAD, and EOP. Four ensemble scores (i.e., BayesDel, VARITY, ClinPred, and MetaRNN) were assessed to predict the deleteriousness of missense mutations and demonstrated significant differences (p < 0.001). Binary logistic regression analyses demonstrated that LOF mutations were independently associated with brain iron accumulation (p = 0.006) and ataxia (p = 0.025).

CONCLUSIONS

LOF or more deleterious missense mutations are more likely to promote the development of serious phenotype of PLAN, and LOF mutations are independently associated with brain iron accumulation and ataxia.

Crucial neuroprotective roles of the metabolite BH4 in dopaminergic neurons

Dopa-responsive dystonia (DRD) and Parkinson's disease (PD) are movement disorders caused by the dysfunction of nigrostriatal dopaminergic neurons. Identifying druggable pathways and biomarkers for guiding therapies is crucial due to the debilitating nature of these disorders. Recent genetic studies have identified variants of GTP cyclohydrolase-1 (GCH1), the rate-limiting enzyme in tetrahydrobiopterin (BH4) synthesis, as causative for these movement disorders. Here, we show that genetic and pharmacological inhibition of BH4 synthesis in mice and human midbrain-like organoids accurately recapitulates motor, behavioral and biochemical characteristics of these human diseases, with severity of the phenotype correlating with extent of BH4 deficiency. We also show that BH4 deficiency increases sensitivities to several PD-related stressors in mice and PD human cells, resulting in worse behavioral and physiological outcomes. Conversely, genetic and pharmacological augmentation of BH4 protects mice from genetically- and chemically induced PD-related stressors. Importantly, increasing BH4 levels also protects primary cells from PD-affected individuals and human midbrain-like organoids (hMLOs) from these stressors. Mechanistically, BH4 not only serves as an essential cofactor for dopamine synthesis, but also independently regulates tyrosine hydroxylase levels, protects against ferroptosis, scavenges mitochondrial ROS, maintains neuronal excitability and promotes mitochondrial ATP production, thereby enhancing mitochondrial fitness and cellular respiration in multiple preclinical PD animal models, human dopaminergic midbrain-like organoids and primary cells from PD-affected individuals. Our findings pinpoint the BH4 pathway as a key metabolic program at the intersection of multiple protective mechanisms for the health and function of midbrain dopaminergic neurons, identifying it as a potential therapeutic target for PD.

A Putative Guanosine Triphosphate Cyclohydrolase I Named CaGCH1 Is Involved in Hyphal Branching and Fruiting Development in Cyclocybe aegerita

Guanosine triphosphate (GTP) cyclohydrolase I (GCH1) is the limiting enzyme of the tetrahydrobiopterin (BH4) synthesis pathway. The disruption of gch1 gene may cause conditional lethality due to folic acid auxotrophy in microorganisms, although the function of gch1 in basidiomycetes has not been deciphered so far. In the present study, gch1 expression in Cyclocybe aegerita (cagch1) was downregulated using the RNAi method, which resulted in growth retardation in both solid and liquid medium, with the hyphal tips exhibiting increased branching compared to that in the wild strain. The development of fruiting bodies in the mutant strains was significantly blocked, and there were short and bottle-shaped stipes. The transcriptional profile revealed that the genes of the MAPK pathway may be involved in the regulation of these effects caused by cagch1 knockdown, which provided an opportunity to study the role of gch1 in the development process of basidiomycetes.

Revealing a novel contributing landscape of ferroptosis-related genes in Parkinson’s disease

Transcriptomics studies have yielded great insights into disease processes by detecting differentially expressed genes (DEGs). In this study, due to the high heritability of Parkinson’s disease (PD), we performed bioinformatics analyses on nine transcriptomic datasets regarding substantia nigra from Gene Expression Omnibus database, including seven microarray datasets and two next-generation sequencing datasets. As a result, between age-matched PD patients and normal control, we identified 630 DEGs, of which 22 hub DEGs involved in PD or ferroptosis were found to be associated with each other at the transcriptional level and protein-protein interaction network, suggesting their high correlations among these hub genes. Moreover, 16 DEGs were singled out due to their comparable AUC (>0.6) in random forest classifiers, including seven PD-related genes (MAP4K4, LRP10, UCHL1, PAM, RIT2, SNCA, GCH1) and nine ferroptosis-related genes (GCH1, DDIT4, RGS4, MAPK9, CAV1, RELA, DUSP1, ATP6V1G2, ATF4 and ISCU). Furthermore, to probe the potential of those hub genes in predicting the PD progression and survival, we constructed a Cox model featured by an eight-gene signature, including four PD-related genes (SNCA, UCHL1, LRP10, and GCH1) and four ferroptosis-related genes (DDIT4, RGS4, RELA, and CAV1), and validated it successful in an independent dataset, indicating that it would be an effective tool for clinical research to predict PD progression. In conclusion, ferroptosis-related DEGs identified in this study were closely correlated with the known PD-related genes, revealing the involvement of ferroptosis in the development of PD. This study presented the potential of several ferroptosis-related genes as novel clinical biomarkers for PD.

Low-frequency and rare coding variants of NUS1 contribute to susceptibility and phenotype of Parkinson's disease

NUS1 has been recently identified as a candidate gene for Parkinson's disease (PD). Few studies have examined the association of NUS1 variants with PD susceptibility and phenotypes. In the first cohort, whole-exome sequencing was performed to identify variants in NUS1 exon-coding and exon-intron regions in 1542 cases and 1625 controls. 13 variants were totally detected, of which 10 rare variants and 3 low-frequency variants. Burden analysis showed that rare NUS1 variants significantly enriched in PD (p=0.016). We also performed a meta-analysis based on previous and our studies to correlate NUS1 mutations with PD susceptibility. Integrating our previous cohort (3210 cases and 2807 controls) and the first cohort identified the significant association of rs539668656 with PD risk (odds ratio (OR) = 2.82, p = 0.016). The genotype-phenotype association analysis showed that patients carrying rare variants, or rs539668656 were significantly associated with earlier onset age, depression, emotional impairment and severe disease condition. Our results support the role of NUS1 rare variants and rs539668656 towards PD susceptibility and phenotype.

Gene4PD: A Comprehensive Genetic Database of Parkinson's Disease

Parkinson’s disease (PD) is a complex neurodegenerative disorder with a strong genetic component. A growing number of variants and genes have been reported to be associated with PD; however, there is no database that integrate different type of genetic data, and support analyzing of PD-associated genes (PAGs). By systematic review and curation of multiple lines of public studies, we integrate multiple layers of genetic data (rare variants and copy-number variants identified from patients with PD, associated variants identified from genome-wide association studies, differentially expressed genes, and differential DNA methylation genes) and age at onset in PD. We integrated five layers of genetic data (8302 terms) with different levels of evidences from more than 3,000 studies and prioritized 124 PAGs with strong or suggestive evidences. These PAGs were identified to be significantly interacted with each other and formed an interconnected functional network enriched in several functional pathways involved in PD, suggesting these genes may contribute to the pathogenesis of PD. Furthermore, we identified 10 genes were associated with a juvenile-onset (age ≤ 30 years), 11 genes were associated with an early-onset (age of 30–50 years), whereas another 10 genes were associated with a late-onset (age > 50 years). Notably, the AAOs of patients with loss of function variants in five genes were significantly lower than that of patients with deleterious missense variants, while patients with VPS13C (P = 0.01) was opposite. Finally, we developed an online database named Gene4PD (http://genemed.tech/gene4pd) which integrated published genetic data in PD, the PAGs, and 63 popular genomic data sources, as well as an online pipeline for prioritize risk variants in PD. In conclusion, Gene4PD provides researchers and clinicians comprehensive genetic knowledge and analytic platform for PD, and would also improve the understanding of pathogenesis in PD.

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.

Characterizing the Expression Patterns of Parkinson's Disease Associated Genes

Background

The expression pattern represents a quantitative phenotype that provides an in-depth view of the molecular mechanism in Parkinson’s disease (PD); however, the expression patterns of PD-associated genes (PAGs) and their relation to age at onset (AAO) remain unclear.

Methods

The known PD-causing genes and PD-risk genes, which were collected from latest published authoritative meta-analysis, were integrated as PAGs. The expression data from Genotype-Tissue Expression database, Allen Brian Map database, and BrainSpan database, were extracted to characterize the tissue specificity, inhibitory-excitatory neuron expression profile, and spatio-temporal expression pattern of PAGs, respectively. The AAO information of PD-causing gene was download from Gene4PD and MDSgene database.

Results

We prioritized 107 PAGs and found that the PAGs were more likely to be expressed in brain-related tissues than non-brain tissues and that more PAGs had higher expression levels in excitatory neurons than inhibitory neurons. In addition, we identified two spatio-temporal expression modules of PAGs in human brain: the first module showed a higher expression level in the adult period than in the prenatal period, and the second module showed the opposite features. It showed that more PAGs belong to the first module that the second module. Furthermore, we found that the median AAO of patients with mutations in PD-causing genes of the first module was lower than that of the second module.

Conclusion

In conclusion, this study provided comprehensive landscape of expression patterns, AAO features and their relationship for the first time, improving the understanding of pathogenesis, and precision medicine in PD.

GCH1 mutations in hereditary spastic paraplegia

GCH1 mutations have been associated with dopa-responsive dystonia (DRD), Parkinson's disease (PD) and tetrahydrobiopterin (BH₄ )-deficient hyperphenylalaninemia B. Recently, GCH1 mutations have been reported in five patients with hereditary spastic paraplegia (HSP). Here, we analyzed a total of 400 HSP patients (291 families) from different centers across Canada by whole exome sequencing (WES). Three patients with heterozygous GCH1 variants were identified: monozygotic twins with a p.(Ser77_Leu82del) variant, and a patient with a p.(Val205Glu) variant. The former variant is predicted to be likely pathogenic and the latter is pathogenic. The three patients presented with childhood-onset lower limb spasticity, hyperreflexia and abnormal plantar responses. One of the patients had diurnal fluctuations, and none had parkinsonism or dystonia. Phenotypic differences between the monozygotic twins were observed, who responded well to levodopa treatment. Pathway enrichment analysis suggested that GCH1 shares processes and pathways with other HSP-associated genes, and structural analysis of the variants indicated a disruptive effect. In conclusion, GCH1 mutations may cause HSP; therefore, we suggest a levodopa trial in HSP patients and including GCH1 in the screening panels of HSP genes. Clinical differences between monozygotic twins suggest that environmental factors, epigenetics, and stochasticity could play a role in the clinical presentation.

Contribution of coding/non-coding variants in NUS1 to late-onset sporadic Parkinson's disease

INTRODUCTION

A recent study reported that rare variants in NUS1 were associated with Parkinson's disease (PD). We aimed to assess the relative contribution of rare and common coding/non-coding variants of NUS1 to late-onset PD patients (LOPD).

METHODS

Whole genome sequencing data were analyzed for target NUS1 regions, derived from a cohort of 1962 cases and 1279 controls. The genetic association analyses were performed using logistic regression analysis and Sequence Kernel association test. Expression quantitative trait loci (eQTL) analysis was conducted to further explore the association of variants with NUS1 expression based on the data from GTEx database.

RESULTS

We identified 18 rare coding variants. p.Y131C was first identified in LOPD. However, no significant burden of rare NUS1 coding variants in LOPD was found. The rare variant sets of two regulatory elements (GH06J117605 and GH06J117674) were significantly enriched in LOPD even after Bonferroni correction (adjusted P = 0.013; adjusted P = 0.010). Considering the joint effect of rare and common variants, all variant sets within GH06J117605 and GH06J117674 showed association with LOPD but were no longer significant after Bonferroni correction. None of the common variants within coding/non-coding regions were significant after Bonferroni correction. The eQTL results suggested these variants in GH06J117605 and GH06J117674 could potentially have eQTL effects on the brain tissues.

CONCLUSIONS

These findings provide novel insight into the role of NUS1 regulatory regions in the development of LOPD and indicate that the variants in regulatory elements of NUS1 may be associated with LOPD by influencing the gene expression level.

The role of genetics in Parkinson's disease: a large cohort study in Chinese mainland population

This study aimed to determine the mutational spectrum of familial Parkinson's disease and sporadic early-onset Parkinson's disease (sEOPD) in a mainland Chinese population and the clinical features of mutation carriers. We performed multiplex ligation-dependent probe amplification assays and whole-exome sequencing for 1676 unrelated patients with Parkinson's disease in a mainland Chinese population, including 192 probands from families with autosomal-recessive Parkinson's disease, 242 probands from families with autosomal-dominant Parkinson's disease, and 1242 sEOPD patients (age at onset ≤ 50). According to standards and guidelines from the American College of Medical Genetics and Genomics, pathogenic/likely pathogenic variants in 23 known Parkinson's disease-associated genes occurred more frequently in the autosomal-recessive Parkinson's disease cohort (65 of 192, 33.85%) than in the autosomal-dominant Parkinson's disease cohort (10 of 242, 4.13%) and the sEOPD cohort (57 of 1242, 4.59%), which leads to an overall molecular diagnostic yield of 7.88% (132 of 1676). We found that PRKN was the most frequently mutated gene (n = 83, 4.95%) and present the first evidence of an SNCA duplication and LRRK2 p.N1437D variant in mainland China. In addition, several novel pathogenic/likely pathogenic variants including LRRK2 (p.V1447M and p.Y1645S), ATP13A2 (p.R735X and p.A819D), FBXO7 (p.G67E), LRP10 (c.322dupC/p.G109Rfs*51) and TMEM230 (c.429delT/p.P144Qfs*2) were identified in our cohort. Furthermore, the age at onset of the 132 probands with genetic diagnoses (median, 31.5 years) was about 14.5 years earlier than that of patients without molecular diagnoses (i.e. non-carriers, median 46.0 years). Specifically, the age at onset of Parkinson's disease patients with pathogenic/likely pathogenic variants in ATP13A2, PLA2G6, PRKN, or PINK1 was significantly lower than that of non-carriers, while the age at onset of carriers with other gene pathogenic/likely pathogenic variants was similar to that of non-carriers. The clinical spectrum of Parkinson's disease-associated gene carriers in this mainland Chinese population was similar to that of other populations. We also detected 61 probands with GBA possibly pathogenic variants (3.64%) and 59 probands with GBA p.L444P (3.52%). These results shed insight into the genetic spectrum and clinical manifestations of Parkinson's disease in mainland China and expand the existing repertoire of pathogenic or likely pathogenic variants involved in known Parkinson's disease-associated genes. Our data highlight the importance of genetic testing in Parkinson's disease patients with age at onset < 40 years, especially in those from families with a recessive inheritance pattern, who may benefit from early diagnosis and treatment.

GCH1 variants contribute to the risk and earlier age-at-onset of Parkinson's disease: a two-cohort case-control study

Background

Common and rare variants of guanosine triphosphate cyclohydrolase 1 (GCH1) gene may play important roles in Parkinson’s disease (PD). However, there is a lack of comprehensive analysis of GCH1 genotypes, especially in non-coding regions. The aim of this study was to explore the genetic characteristics of GCH1, including rare and common variants in coding and non-coding regions, in a large population of PD patients in Chinese mainland, as well as the phenotypic characteristics of GCH1 variant carriers.

Methods

In the first cohort of this case-control study, we performed whole-exome sequencing in 1555 patients with early-onset or familial PD and 2234 healthy controls; then in the second cohort, whole-genome sequencing was performed in sporadic late-onset PD samples (1962 patients), as well as 1279 controls. Variants at target GCH1 regions were extracted, and then genetic and detailed phenotypic data were analyzed using regression models and the sequence kernel association test. We also performed a meta-analysis to correlate deleterious GCH1 variants with age at onset (AAO) in PD patients.

Results

For coding variants, we identified a significant burden of GCH1 deleterious variants in early-onset or familial PD cases compared to controls (1.2% vs 0.1%, P < 0.0001). In the analysis of possible regulatory variants in GCH1 non-coding regions, rs12323905 (P = 0.001, odds ratio = 1.19, 95%CI 1.07–1.32) was significantly associated with PD, and variant sets in untranslated regions and intron regions, GCH1 brain-specific expression quantitative trait loci, and two possible promoter/enhancer (GH14J054857 and GH14J054880) were suggestively associated with PD. Genotype-phenotype correlation analysis revealed that the carriers of GCH1 deleterious variants manifested younger AAO (P < 0.0001), and had milder motor symptoms, milder fatigue symptoms and more autonomic nervous dysfunctions. Meta-analysis of six studies demonstrated 6.4-year earlier onset in GCH1 deleterious variant carriers (P = 0.0009).

Conclusions

The results highlight the importance of deleterious variants and non-coding variants of GCH1 in PD in Chinese mainland and suggest that GCH1 mutation can influence the PD phenotype, which may help design experimental studies to elucidate the mechanisms of GCH1 in the pathogenesis of PD.

Genetic analysis of N6-methyladenosine modification genes in Parkinson's disease

Parkinson's disease (PD) is a common neurodegenerative disease with a relatively unclear etiology. Previous studies have shown that N6-methyladenosine (m6A) is a vital RNA modification enriched in brain tissue, and that the genes involved in m6A modification are implicated in various neurologic diseases. Here, we conducted a comprehensive genetic analysis using targeted sequencing with molecular inversion probes (MIPs) to identify m6A-modification genes (including METTL3, METTL14, WTAP, FTO, ALKBH5, YTHDF1, YTHDF2, YTHDF3, HNRNPC, and ELAVL1) in a total of 1647 sporadic PD patients and 1372 controls of Han Chinese origin. PD patients were divided into early-onset PD (EOPD) and late-onset PD (LOPD) based on whether the onset of motor symptoms occurred before or after 50 years of age. Rare variants were subjected to gene-based burden tests and common variants were subjected to single-variant association analyses. As a result, we identified 214 rare variants in all 10 m6A-modification genes and 16 common variants in 7 genes. Gene-wise association analyses of rare variants in each m6A-modification gene did not achieved a p value of less than 0.05 in either total cohorts or 2 age groups. In fact, p values greater than 0.05 were found when conducting single-variant association analyses on common variants of these genes between PD and control patients. Our comprehensive analyses of m6A-modification genes suggest that there is no significant association between these 10 m6A-modification genes and the risk of sporadic PD.

A Survey of Gene Prioritization Tools for Mendelian and Complex Human Diseases

Modern high-throughput experiments provide us with numerous potential associations between genes and diseases. Experimental validation of all the discovered associations, let alone all the possible interactions between them, is time-consuming and expensive. To facilitate the discovery of causative genes, various approaches for prioritization of genes according to their relevance for a given disease have been developed. In this article, we explain the gene prioritization problem and provide an overview of computational tools for gene prioritization. Among about a hundred of published gene prioritization tools, we select and briefly describe 14 most up-to-date and user-friendly. Also, we discuss the advantages and disadvantages of existing tools, challenges of their validation, and the directions for future research.

Genetic analysis of DNA methylation and hydroxymethylation genes in Parkinson's disease

DNA methylation is an important regulatory mechanism of Parkinson's disease (PD). To investigate the relationship between DNA methylation and hydroxymethylation genes and PD, we performed gene-targeted sequencing using molecular inversion probes in a Chinese PD population. We sequenced 12 genes related to DNA methylation and hydroxymethylation in 1657 patients and 1394 control subjects. We conducted genewise association analyses of rare variants detected in the present study and identified the TET1 gene as important in PD (p = 0.0037738, 0.013, 0.019521 (b.collapse test, variable threshold test, and skat-o test, respectively; sex + age as covariates). However, no positive results were observed when conducting association analyses on common variants in these genes. We performed a comprehensive analysis of associations between variants of DNA methylation and hydroxymethylation genes and PD, resulting in determination that TET1 might play a role in PD.

Common and rare GCH1 variants are associated with Parkinson's disease

GCH1 encodes the enzyme guanosine triphospahte (GTP) cyclohydrolase 1, essential for dopamine synthesis in nigrostriatal cells, and rare mutations in GCH1 may lead to Dopa-responsive dystonia (DRD). While GCH1 is implicated in genomewide association studies in Parkinson's disease (PD), only a few studies examined the role of rare GCH1 variants in PD, with conflicting results. In the present study, GCH1 and its 5' and 3' untranslated regions were sequenced in 1113 patients with PD and 1111 controls. To examine the association of rare GCH1 variants with PD, burden analysis was performed. Three rare GCH1 variants, which were previously reported to be pathogenic in DRD, were found in five patients with PD and not in controls (sequence Kernel association test, p = 0.024). A common haplotype, tagged by rs841, was associated with a reduced risk for PD (OR = 0.71, 95% CI = 0.61-0.83, p = 1.24 × 10^-4), and with increased GCH1 expression in brain regions relevant for PD (www.gtexportal.org). Our results support a role for rare, DRD-related variants, and common GCH1 variants in the pathogenesis of PD.