Investigation of Shared Genetic Risk Factors Between Parkinson's Disease and Cancers

in Movement Disorders (2023), n/a(n/a),

Abstract Background Epidemiological studies that examined the association between Parkinson's disease (PD) and cancers led to inconsistent results, but they face a number of methodological difficulties ... [more ▼]

Abstract Background Epidemiological studies that examined the association between Parkinson's disease (PD) and cancers led to inconsistent results, but they face a number of methodological difficulties. Objective We used results from genome-wide association studies (GWASs) to study the genetic correlation between PD and different cancers to identify common genetic risk factors. Methods We used individual data for participants of European ancestry from the Courage-PD (Comprehensive Unbiased Risk Factor Assessment for Genetics and Environment in Parkinson's Disease; PD, N = 16,519) and EPITHYR (differentiated thyroid cancer, N = 3527) consortia and summary statistics of GWASs from iPDGC (International Parkinson Disease Genomics Consortium; PD, N = 482,730), Melanoma Meta-Analysis Consortium (MMAC), Breast Cancer Association Consortium (breast cancer), the Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome (prostate cancer), International Lung Cancer Consortium (lung cancer), and Ovarian Cancer Association Consortium (ovarian cancer) (N comprised between 36,017 and 228,951 for cancer GWASs). We estimated the genetic correlation between PD and cancers using linkage disequilibrium score regression. We studied the association between PD and polymorphisms associated with cancers, and vice versa, using cross-phenotypes polygenic risk score (PRS) analyses. Results We confirmed a previously reported positive genetic correlation of PD with melanoma (Gcorr = 0.16 [0.04; 0.28]) and reported an additional significant positive correlation of PD with prostate cancer (Gcorr = 0.11 [0.03; 0.19]). There was a significant inverse association between the PRS for ovarian cancer and PD (odds ratio [OR] = 0.89 [0.84; 0.94]). Conversely, the PRS of PD was positively associated with breast cancer (OR = 1.08 [1.06; 1.10]) and inversely associated with ovarian cancer (OR = 0.95 [0.91; 0.99]). The association between PD and ovarian cancer was mostly driven by rs183211 located in an intron of the NSF gene (17q21.31). Conclusions We show evidence in favor of a contribution of pleiotropic genes to the association between PD and specific cancers. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA. [less ▲]

Mitochondrial DNA heteroplasmy distinguishes disease manifestation in PINK1/PRKN-linked Parkinson’s disease

in Brain (2022)

Biallelic mutations in PINK1/PRKN cause recessive Parkinson’s disease. Given the established role of PINK1/Parkin in regulating mitochondrial dynamics, we explored mitochondrial DNA (mtDNA) integrity and ... [more ▼]

Biallelic mutations in PINK1/PRKN cause recessive Parkinson’s disease. Given the established role of PINK1/Parkin in regulating mitochondrial dynamics, we explored mitochondrial DNA (mtDNA) integrity and inflammation as disease modifiers in carriers of mutations in these genes. MtDNA integrity was investigated in a large collection of biallelic (n = 84) and monoallelic (n = 170) carriers of PINK1/PRKN mutations, idiopathic Parkinson’s disease patients (n = 67) and controls (n = 90). In addition, we studied global gene expression and serum cytokine levels in a subset. Affected and unaffected PINK1/PRKN monoallelic mutation carriers can be distinguished by heteroplasmic mtDNA variant load (AUC = 0.83, CI:0.74-0.93). Biallelic PINK1/PRKN mutation carriers harbor more heteroplasmic mtDNA variants in blood (p = 0.0006, Z = 3.63) compared to monoallelic mutation carriers. This enrichment was confirmed in iPSC-derived (controls, n = 3; biallelic PRKN mutation carriers, n = 4) and postmortem (control, n = 1; biallelic PRKN mutation carrier, n = 1) midbrain neurons. Lastly, the heteroplasmic mtDNA variant load correlated with IL6 levels in PINK1/PRKN mutation carriers (r = 0.57, p = 0.0074). PINK1/PRKN mutations predispose individuals to mtDNA variant accumulation in a dose- and disease-dependent manner. [less ▲]

Forecasting of a complex microbial community using meta-omics

E-print/Working paper (2022)

Microbial communities are complex assemblages whose dynamics are shaped by abiotic and biotic factors. A major challenge concerns correctly forecasting the community behaviour in the future. In this ... [more ▼]

Microbial communities are complex assemblages whose dynamics are shaped by abiotic and biotic factors. A major challenge concerns correctly forecasting the community behaviour in the future. In this context, communities in biological wastewater treatment plants (BWWTPs) represent excellent model systems, because forecasting them is required to ultimately control and operate the plants in a sustainable manner. Here, we forecast the microbial community from the water-air interface of the anaerobic tank of a BWWTP via longitudinal meta-omics (metagenomics, metatranscriptomics and metaproteomics) data covering 14 months at weekly intervals. We extracted all the available time-dependent information, summarised it in 17 temporal signals (explaining 91.1 of the temporal variance) and linked them over time to rebuild the sequence of ecological phenomena behind the community dynamics. We forecasted the signals over the following five years and tested the predictions with 21 extra samples. We were able to correctly forecast five signals accounting for 22.5 of the time-dependent information in the system and generate mechanistic predictions on the ecological events in the community (e.g. a predation cycle involving bacteria, viruses and amoebas). Through the forecasting of the 17 signals and the environmental variables readings we reconstructed the gene abundance and expression for the following 5 years, showing a nearly perfect trend prediction (coefficient of determination >= 0.97) for the first 2 years. The study demonstrates the maturity of microbial ecology to forecast composition and gene expression of open microbial ecosystems using year-spanning interactions between community cycles and environmental parameters. [less ▲]

binny: an automated binning algorithm to recover high-quality genomes from complex metagenomic datasets

in Briefings in Bioinformatics (2022)

The reconstruction of genomes is a critical step in genome-resolved metagenomics and for multi-omic data integration from microbial communities. Here, we present binny, a binning tool that produces high ... [more ▼]

The reconstruction of genomes is a critical step in genome-resolved metagenomics and for multi-omic data integration from microbial communities. Here, we present binny, a binning tool that produces high-quality metagenome-assembled genomes (MAG) from both contiguous and highly fragmented genomes. Based on established metrics, binny outperforms or is highly competitive with commonly used and state-of-the-art binning methods and finds unique genomes that could not be detected by other methods. binny uses k-mer-composition and coverage by metagenomic reads for iterative, nonlinear dimension reduction of genomic signatures as well as subsequent automated contig clustering with cluster assessment using lineage-specific marker gene sets. When compared with seven widely used binning algorithms, binny provides substantial amounts of uniquely identified MAGs and almost always recovers the most near-complete (⁠>95% pure, >90% complete) and high-quality (⁠>90% pure, >70% complete) genomes from simulated datasets from the Critical Assessment of Metagenome Interpretation initiative, as well as substantially more high-quality draft genomes, as defined by the Minimum Information about a Metagenome-Assembled Genome standard, from a real-world benchmark comprised of metagenomes from various environments than any other tested method. [less ▲]

Genotype-phenotype correlations in SCN8A-related disorders reveal prognostic and therapeutic implications

in Brain (2022)

We report detailed functional analyses and genotype-phenotype correlations in 392 individuals carrying disease-causing variants in SCN8A, encoding the voltage-gated Na+ channel NaV1.6, with the aim of ... [more ▼]

We report detailed functional analyses and genotype-phenotype correlations in 392 individuals carrying disease-causing variants in SCN8A, encoding the voltage-gated Na+ channel NaV1.6, with the aim of describing clinical phenotypes related to functional effects. Six different clinical subgroups could be identified: 1) Benign familial infantile epilepsy (BFIE) (n = 15, normal cognition, treatable seizures), 2) intermediate epilepsy (n = 33, mild ID, partially pharmaco-responsive), 3) developmental and epileptic encephalopathy (DEE, n = 177, severe ID, majority pharmaco-resistant), 4) generalized epilepsy (n = 20, mild to moderate ID, frequently with absence seizures), 5) unclassifiable epilepsy (n = 127), and 6) neurodevelopmental disorder without epilepsy (n = 20, mild to moderate ID). Groups 1–3 presented with focal or multifocal seizures (median age of onset: four months) and focal epileptiform discharges, whereas the onset of seizures in group 4 was later (median: 42 months) with generalized epileptiform discharges. We performed functional studies expressing missense variants in ND7/23 neuroblastoma cells and primary neuronal cultures using recombinant tetrodotoxin-insensitive human NaV1.6 channels and whole-cell patch-clamping. Two variants causing DEE showed a strong gain-of-function (GOF, hyperpolarising shift of steady-state activation, strongly increased neuronal firing rate), and one variant causing BFIE or intermediate epilepsy showed a mild GOF (defective fast inactivation, less increased firing). In contrast, all three variants causing generalized epilepsy induced a loss-of-function (LOF, reduced current amplitudes, depolarising shift of steady-state activation, reduced neuronal firing). Including previous studies, functional effects were known for 170 individuals. All 136 individuals carrying a functionally tested GOF variant had either focal (97, groups 1–3), or unclassifiable epilepsy (39), whereas 34 with a LOF variant had either generalized (14), no (11) or unclassifiable (6) epilepsy; only three had DEE. Computational modeling in the GOF group revealed a significant correlation between the severity of the electrophysiological and clinical phenotypes. GOF variant carriers responded significantly better to sodium channel blockers (SCBs) than to other anti-seizure medications, and the same applied for all individuals of groups 1–3.In conclusion, our data reveal clear genotype-phenotype correlations between age at seizure onset, type of epilepsy and gain- or loss-of-function effects of SCN8A variants. Generalized epilepsy with absence seizures is the main epilepsy phenotype of LOF variant carriers and the extent of the electrophysiological dysfunction of the GOF variants is a main determinant of the severity of the clinical phenotype in focal epilepsies. Our pharmacological data indicate that SCBs present a treatment option in SCN8A-related focal epilepsy with onset in the first year of life. [less ▲]

Genetic stratification of motor and QoL outcomes in Parkinson's disease in the EARLYSTIM study

in Parkinsonism and Related Disorders (2022)

Purpose The decision for subthalamic deep brain stimulation (STN-DBS) in Parkinson's disease (PD) relies on clinical predictors. Whether genetic variables could predict favourable or unfavourable ... [more ▼]

Purpose The decision for subthalamic deep brain stimulation (STN-DBS) in Parkinson's disease (PD) relies on clinical predictors. Whether genetic variables could predict favourable or unfavourable decisions is under investigation. Objective First, we aimed to reproduce the previous observation that SNCA rs356220 was associated with favourable STN-DBS motor response. In additional exploratory analyses, we studied if other PD risk and progression variants from the latest GWAS are associated with therapeutic outcome. Further, we evaluated the predictive value of polygenic risk scores. Methods We comprehensively genotyped patients from the EarlyStim cohort using NeuroChip, and assessed the clinico-genetic associations with longitudinal outcome parameters. Results The SNCA rs356220 variant did not predict UPDRS III outcomes. However, it was associated with quality of life improvement in secondary analyses. Several polymorphisms from previously identified GWAS hits predicted motor or quality of life outcomes in DBS patients. Polygenic risk scores did not predict any outcome parameter. Conclusions Our findings support the hypothesis that different common genetic markers are associated with favourable quality of life outcomes of STN-DBS in PD. These findings can be the basis for further validation in larger and independent cohorts. [less ▲]

Conserved patterns across ion channels correlate with variant pathogenicity and clinical phenotypes

in Brain: a Journal of Neurology (2022)

Clinically identified genetic variants in ion channels can be benign or cause disease by increasing or decreasing the protein function. Consequently, therapeutic decision-making is challenging without ... [more ▼]

Clinically identified genetic variants in ion channels can be benign or cause disease by increasing or decreasing the protein function. Consequently, therapeutic decision-making is challenging without molecular testing of each variant. Our biophysical knowledge of ion channel structures and function is just emerging, and it is currently not well understood which amino acid residues cause disease when mutated.We sought to systematically identify biological properties associated with variant pathogenicity across all major voltage and ligand-gated ion channel families. We collected and curated 3,049 pathogenic variants from hundreds of neurodevelopmental and other disorders and 12,546 population variants for 30 ion channel or channel subunits for which a high-quality protein structure was available. Using a wide range of bioinformatics approaches, we computed 163 structural features and tested them for pathogenic variant enrichment. We developed a novel 3D spatial distance scoring approach that enables comparisons of pathogenic and population variant distribution across protein structures.We discovered and independently replicated that several pore residue properties and proximity to the pore axis were most significantly enriched for pathogenic variants compared to population variants. Using our 3D scoring approach, we showed that the strongest pathogenic variant enrichment was observed for pore-lining residues and alpha-helix residues within 5Å distance from the pore axis center and not involved in gating. Within the subset of residues located at the pore, the hydrophobicity of the pore was the feature most strongly associated with variant pathogenicity. We also found an association between the identified properties and both clinical phenotypes and functional in vitro assays for voltage-gated sodium channels (SCN1A, SCN2A, SCN8A) and N-methyl-D-aspartate (NMDA) receptor (GRIN1, GRIN2A, GRIN2B) encoding genes. In an independent expert-curated dataset of 1,422 neurodevelopmental disorder pathogenic patient variants and 679 electrophysiological experiments, we show that pore axis distance is associated with seizure age of onset and cognitive performance as well as differential gain vs. loss-of-channel function.In summary, we identified biological properties associated with ion-channel malfunction and show that these are correlated with in vitro functional read-outs and clinical phenotypes in patients with neurodevelopmental disorders. Our results suggest that clinical decision support algorithms that predict variant pathogenicity and function are feasible in the future. [less ▲]

Early-to-mid idiopathic Parkinson’s disease shows a more cytotoxic but declined CD8-regulatory peripheral immune profile

E-print/Working paper (2022)

Parkinson’s disease (PD) is the second most common neurodegenerative disease. Brain neuroinflammation plays a role in PD pathogenesis. However, the involvement of the peripheral immune system has not been ... [more ▼]

Parkinson’s disease (PD) is the second most common neurodegenerative disease. Brain neuroinflammation plays a role in PD pathogenesis. However, the involvement of the peripheral immune system has not been systematically investigated. Here we analyzed >700 combinatorial immunological features in fresh blood of 28 early-to-mid-stage PD patients and 24 matched controls. We found an enhanced cytotoxic immune profile in idiopathic PD patients (iPD), with a higher frequency of terminally-differentiated effector CD8 T (TEMRA), late-differentiated CD8+ natural killer T cells and neutrophils. This immune profile was intensified by elevated serum granzyme A, reduced percentages of CD8+FOXP3+ regulatory T cells and group 2 innate lymphoid cells with immunosuppressive or tolerance-inducing functions. The frequency of CD8 TEMRA was negatively correlated with disease duration, suggesting a contribution to PD pathogenesis. Our work provides a comprehensive map on disturbed peripheral adaptive and innate immune cells in early-to-mid iPD, proposing easily-accessible candidates for early diagnosis and treatments. [less ▲]

The Interaction between HLA-DRB1 and Smoking in Parkinson's Disease Revisited

in Movement Disorders (2022)

Abstract Background Two studies that examined the interaction between HLA-DRB1 and smoking in Parkinson's disease (PD) yielded findings in opposite directions. Objective To perform a large-scale ... [more ▼]

Abstract Background Two studies that examined the interaction between HLA-DRB1 and smoking in Parkinson's disease (PD) yielded findings in opposite directions. Objective To perform a large-scale independent replication of the HLA-DRB1 × smoking interaction. Methods We genotyped 182 single nucleotide polymorphism (SNPs) associated with smoking initiation in 12 424 cases and 9480 controls to perform a Mendelian randomization (MR) analysis in strata defined by HLA-DRB1. Results At the amino acid level, a valine at position 11 (V11) in HLA-DRB1 displayed the strongest association with PD. MR showed an inverse association between genetically predicted smoking initiation and PD only in absence of V11 (odds ratio, 0.74, 95 confidence interval, 0.59–0.93, PInteraction = 0.028). In silico predictions of the influence of V11 and smoking-induced modifications of α-synuclein on binding affinity showed findings consistent with this interaction pattern. Conclusions Despite being one of the most robust findings in PD research, the mechanisms underlying the inverse association between smoking and PD remain unknown. Our findings may help better understand this association. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society [less ▲]

Genome-wide meta-analysis of over 29,000 people with epilepsy reveals 26 loci and subtype-specific genetic architecture

E-print/Working paper (2022)

Epilepsy is a highly heritable disorder affecting over 50 million people worldwide, of which about one-third are resistant to current treatments. Here, we report a trans-ethnic GWAS including 29,944 cases ... [more ▼]

Epilepsy is a highly heritable disorder affecting over 50 million people worldwide, of which about one-third are resistant to current treatments. Here, we report a trans-ethnic GWAS including 29,944 cases, stratified into three broad- and seven sub-types of epilepsy, and 52,538 controls. We identify 26 genome-wide significant loci, 19 of which are specific to genetic generalized epilepsy (GGE). We implicate 29 likely causal genes underlying these 26 loci. SNP-based heritability analyses show that common variants substantially close the missing heritability gap for GGE. Subtype analysis revealed markedly different genetic architectures between focal and generalized epilepsies. Gene-set analysis of GGE signals implicate synaptic processes in both excitatory and inhibitory neurons in the brain. Prioritized candidate genes overlap with monogenic epilepsy genes and with targets of current anti-seizure medications. Finally, we leverage our results to identify alternate drugs with predicted efficacy if repurposed for epilepsy treatment. [less ▲]