Establishing gene regulatory networks from Parkinson's disease risk loci.

The latest meta-analysis of genome-wide association studies identified 90 independent variants across 78 genomic regions associated with Parkinson's disease, yet the mechanisms by which these variants influence the development of the disease remains largely elusive. To establish the functional gene regulatory networks associated with Parkinson's disease risk variants, we utilized an approach combining spatial (chromosomal conformation capture) and functional (expression quantitative trait loci) data. We identified 518 genes subject to regulation by 76 Parkinson's variants across 49 tissues, whicih encompass 36 peripheral and 13 CNS tissues. Notably, one-third of these genes were regulated via trans-acting mechanisms (distal; risk locus-gene separated by >1 Mb, or on different chromosomes). Of particular interest is the identification of a novel trans-expression quantitative trait loci-gene connection between rs10847864 and SYNJ1 in the adult brain cortex, highlighting a convergence between familial studies and Parkinson's disease genome-wide association studies loci for SYNJ1 (PARK20) for the first time. Furthermore, we identified 16 neurodevelopment-specific expression quantitative trait loci-gene regulatory connections within the foetal cortex, consistent with hypotheses suggesting a neurodevelopmental involvement in the pathogenesis of Parkinson's disease. Through utilizing Louvain clustering we extracted nine significant and highly intraconnected clusters within the entire gene regulatory network. The nine clusters are enriched for specific biological processes and pathways, some of which have not previously been associated with Parkinson's disease. Together, our results not only contribute to an overall understanding of the mechanisms and impact of specific combinations of Parkinson's disease variants, but also highlight the potential impact gene regulatory networks may have when elucidating aetiological subtypes of Parkinson's disease.

Regulation of mitochondrial metabolism in murine skeletal muscle by the medium-chain fatty acid receptor Gpr84.

Fatty acid receptors have been recognized as important players in glycaemic control. This study is the first to describe a role for the medium-chain fatty acid (MCFA) receptor G-protein-coupled receptor (Gpr) 84 in skeletal muscle mitochondrial function and insulin secretion. We are able to show that Gpr84 is highly expressed in skeletal muscle and adipose tissue. Mice with global deletion of Gpr84 [Gpr84 knockout (KO)] exhibit a mild impairment in glucose tolerance when fed a MCFA-enriched diet. Studies in mice and pancreatic islets suggest that glucose intolerance is accompanied by a defect in insulin secretion. MCFA-fed KO mice also exhibit a significant impairment in the intrinsic respiratory capacity of their skeletal muscle mitochondria, but at the same time also exhibit a substantial increase in mitochondrial content. Changes in canonical pathways of mitochondrial biogenesis and turnover are unable to explain these mitochondrial differences. Our results show that Gpr84 plays a crucial role in regulating mitochondrial function and quality control.-Montgomery, M. K., Osborne, B., Brandon, A. E., O'Reilly, L., Fiveash, C. E., Brown, S. H. J., Wilkins, B. P., Samsudeen, A., Yu, J., Devanapalli, B., Hertzog, A., Tolun, A. A., Kavanagh, T., Cooper, A. A., Mitchell, T. W., Biden, T. J., Smith, N. J., Cooney, G. J., Turner, N. Regulation of mitochondrial metabolism in murine skeletal muscle by the medium-chain fatty acid receptor Gpr84.

Parkinson's disease-linked human PARK9/ATP13A2 maintains zinc homeostasis and promotes α-Synuclein externalization via exosomes.

α-Synuclein plays a central causative role in Parkinson's disease (PD). Increased expression of the P-type ATPase ion pump PARK9/ATP13A2 suppresses α-Synuclein toxicity in primary neurons. Our data indicate that ATP13A2 encodes a zinc pump; neurospheres from a compound heterozygous ATP13A2(-/-) patient and ATP13A2 knockdown cells are sensitive to zinc, whereas ATP13A2 over-expression in primary neurons confers zinc resistance. Reduced ATP13A2 expression significantly decreased vesicular zinc levels, indicating ATP13A2 facilitates transport of zinc into membrane-bound compartments or vesicles. Endogenous ATP13A2 localized to multi-vesicular bodies (MVBs), a late endosomal compartment located at the convergence point of the endosomal and autophagic pathways. Dysfunction in MVBs can cause a range of detrimental effects including lysosomal dysfunction and impaired delivery of endocytosed proteins/autophagy cargo to the lysosome, both of which have been observed in cells with reduced ATP13A2 function. MVBs also serve as the source of intra-luminal nanovesicles released extracellularly as exosomes that can contain a range of cargoes including α-Synuclein. Elevated ATP13A2 expression reduced intracellular α-Synuclein levels and increased α-Synuclein externalization in exosomes >3-fold whereas ATP13A2 knockdown decreased α-Synuclein externalization. An increased export of exosome-associated α-Synuclein may explain why surviving neurons of the substantia nigra pars compacta in sporadic PD patients were observed to over-express ATP13A2. We propose ATP13A2's modulation of zinc levels in MVBs can regulate the biogenesis of exosomes capable of containing α-Synuclein. Our data indicate that ATP13A2 is the first PD-associated gene involved in exosome biogenesis and indicates a potential neuroprotective role of exosomes in PD.

Identification of 27 allele-specific regulatory variants in Parkinson's disease using a massively parallel reporter assay.

Genome wide association studies (GWAS) have identified a number of genomic loci that are associated with Parkinson's disease (PD) risk. However, the majority of these variants lie in non-coding regions, and thus the mechanisms by which they influence disease development, and/or potential subtypes, remain largely elusive. To address this, we used a massively parallel reporter assay (MPRA) to screen the regulatory function of 5254 variants that have a known or putative connection to PD. We identified 138 loci with enhancer activity, of which 27 exhibited allele-specific regulatory activity in HEK293 cells. The identified regulatory variant(s) typically did not match the original tag variant within the PD associated locus, supporting the need for deeper exploration of these loci. The existence of allele specific transcriptional impacts within HEK293 cells, confirms that at least a subset of the PD associated regions mark functional gene regulatory elements. Future functional studies that confirm the putative targets of the empirically verified regulatory variants will be crucial for gaining a greater understanding of how gene regulatory network(s) modulate PD risk.

Deep sequencing of proteotoxicity modifier genes uncovers a Presenilin-2/beta-amyloid-actin genetic risk module shared among alpha-synucleinopathies.

Whether neurodegenerative diseases linked to misfolding of the same protein share genetic risk drivers or whether different protein-aggregation pathologies in neurodegeneration are mechanistically related remains uncertain. Conventional genetic analyses are underpowered to address these questions. Through careful selection of patients based on protein aggregation phenotype (rather than clinical diagnosis) we can increase statistical power to detect associated variants in a targeted set of genes that modify proteotoxicities. Genetic modifiers of alpha-synuclein (ɑS) and beta-amyloid (Aß) cytotoxicity in yeast are enriched in risk factors for Parkinson's disease (PD) and Alzheimer's disease (AD), respectively. Here, along with known AD/PD risk genes, we deeply sequenced exomes of 430 ɑS/Aß modifier genes in patients across alpha-synucleinopathies (PD, Lewy body dementia and multiple system atrophy). Beyond known PD genes GBA1 and LRRK2, rare variants AD genes (CD33, CR1 and PSEN2) and Aß toxicity modifiers involved in RhoA/actin cytoskeleton regulation (ARGHEF1, ARHGEF28, MICAL3, PASK, PKN2, PSEN2) were shared risk factors across synucleinopathies. Actin pathology occurred in iPSC synucleinopathy models and RhoA downregulation exacerbated ɑS pathology. Even in sporadic PD, the expression of these genes was altered across CNS cell types. Genome-wide CRISPR screens revealed the essentiality of PSEN2 in both human cortical and dopaminergic neurons, and PSEN2 mutation carriers exhibited diffuse brainstem and cortical synucleinopathy independent of AD pathology. PSEN2 contributes to a common-risk signal in PD GWAS and regulates ɑS expression in neurons. Our results identify convergent mechanisms across synucleinopathies, some shared with AD.

Redefining the hypotheses driving Parkinson's diseases research.

Parkinson's disease (PD) research has largely focused on the disease as a single entity centred on the development of neuronal pathology within the central nervous system. However, there is growing recognition that PD is not a single entity but instead reflects multiple diseases, in which different combinations of environmental, genetic and potential comorbid factors interact to direct individual disease trajectories. Moreover, an increasing body of recent research implicates peripheral tissues and non-neuronal cell types in the development of PD. These observations are consistent with the hypothesis that the initial causative changes for PD development need not occur in the central nervous system. Here, we discuss how the use of neuronal pathology as a shared, qualitative phenotype minimises insights into the possibility of multiple origins and aetiologies of PD. Furthermore, we discuss how considering PD as a single entity potentially impairs our understanding of the causative molecular mechanisms, approaches for patient stratification, identification of biomarkers, and the development of therapeutic approaches to PD. The clear consequence of there being distinct diseases that collectively form PD, is that there is no single biomarker or treatment for PD development or progression. We propose that diagnosis should shift away from the clinical definitions, towards biologically defined diseases that collectively form PD, to enable informative patient stratification. N-of-one type, clinical designs offer an unbiased, and agnostic approach to re-defining PD in terms of a group of many individual diseases.

Identifying robust and reliable volatile organic compounds in human sebum for biomarker discovery.

Sebum from sebaceous glands is a rich source of volatile organic compounds (VOCs) that can readily be sampled non-invasively from the surface of skin. The VOC profiles of sebum can then be used to obtain information regarding different medical conditions including diabetes and Parkinson's Disease. However, the effects of sampling approaches and environmental factors on sebum VOC profiles are not established and the confident attribution of VOCs to disease states needs to be free of extraneous influences such as sampling materials and preparatory conditions. Here, we investigated a more standardised skin swab sampling approach for profiling sebum VOCs from healthy human subjects using thermal desorption-gas chromatography-mass spectrometry (TD-GC-MS). Using a standard GC-MS method for the chemical analysis of sebum swabs, a surprisingly high number of VOCs originate from 'blank' medical swab material alone (up to 74 VOCs) and from the ambient environment (up to 29 VOCs) based on control experiments. We found that heat-treatment of medical swabs prior to GC-MS reduced the number of VOCs detected from 'blank' swabs and improved the reproducibility of VOC profiling, however significant VOC absorption can still occur from environmental exposure to ambient air. VOCs identified in 'blank' swabs consisted predominantly of hydrocarbons, esters, and silicon-based compounds and depended strongly on the material used (cotton and polyester-rayon). Environmental VOCs found to absorb to swabs from the ambient air during sampling included 1-butylheptyl-benzene and hexadecanoic acid methyl ester as well as exogenous VOCs such as isopropyl palmitate and isopropyl myristate. In contrast, sebum VOCs consisted primarily of esters, alcohols, ketones, and aldehydes. 23 and 18 VOCs were identified in sebum collected using polyester-rayon and cotton-based medical swabs, respectively, with 14 VOCs common to both swabs. The effect of subject bathing prior to sebum sampling had minimal impact on the VOC profiles. However, individual differences owing to external factors such as skin type, diet, and exercise will likely influence sebum production. This study highlights the importance of using rigorous controls in sebum sampling, and recommendations are provided for future research involving sebum VOC analysis. For example, the use of sebum sample replicates across multiple days, and the use of control swabs during sample collection is required to confirm the origin and reliability of sebum VOCs. It is anticipated that these recommendations in conjunction with a library of well-established VOCs from medical swabs will further strengthen biomarker identification resulting from sebum VOC analysis.

Pathogenic effects of novel mutations in the P-type ATPase ATP13A2 (PARK9) causing Kufor-Rakeb syndrome, a form of early-onset parkinsonism.

Kufor-Rakeb syndrome (KRS) is a rare form of autosomal recessive juvenile or early-onset, levodopa responsive parkinsonism and has been associated with mutations in ATP13A2(also known as PARK9), a lysosomal type 5 P-type ATPase. Recently, we identified novel compound heterozygous mutations, c.3176T>G (p.L1059R) and c.3253delC (p.L1085WfsX1088) in ATP13A2 of two siblings affected with KRS. When overexpressed, wild-type ATP13A2 localized to Lysotracker-positive and LAMP2-positive lysosomes while both truncating and missense mutated ATP13A2 were retained in the endoplasmic reticulum (ER). Both mutant proteins were degraded by the proteasomal but not the lysosomal pathways. In addition, ATP13A2 mRNA with c.3253delC was degraded by nonsense-mediated mRNA decay (NMD), which was protected by cycloheximide treatment. To validate our findings in a biologically relevant setting, we used patient-derived human olfactory neurosphere cultures and fibroblasts and demonstrated persistent ER stress by detecting upregulation of unfolded protein response-related genes in the patient-derived cells. We also confirmed NMD degraded ATP13A2 c.3253delC mRNA in the cells. These findings indicate that these novel ATP13A2 mutations are indeed pathogenic and support the notion that mislocalization of the mutant ATP13A2, resultant ER stress, alterations in the proteasomal pathways and premature degradation of mutant ATP13A2 mRNA contribute to the aetiology of KRS.

Machine Learning Identifies Six Genetic Variants and Alterations in the Heart Atrial Appendage as Key Contributors to PD Risk Predictivity.

Parkinson's disease (PD) is a complex neurodegenerative disease with a range of causes and clinical presentations. Over 76 genetic loci (comprising 90 SNPs) have been associated with PD by the most recent GWAS meta-analysis. Most of these PD-associated variants are located in non-coding regions of the genome and it is difficult to understand what they are doing and how they contribute to the aetiology of PD. We hypothesised that PD-associated genetic variants modulate disease risk through tissue-specific expression quantitative trait loci (eQTL) effects. We developed and validated a machine learning approach that integrated tissue-specific eQTL data on known PD-associated genetic variants with PD case and control genotypes from the Wellcome Trust Case Control Consortium. In so doing, our analysis ranked the tissue-specific transcription effects for PD-associated genetic variants and estimated their relative contributions to PD risk. We identified roles for SNPs that are connected with INPP5P, CNTN1, GBA and SNCA in PD. Ranking the variants and tissue-specific eQTL effects contributing most to the machine learning model suggested a key role in the risk of developing PD for two variants (rs7617877 and rs6808178) and eQTL associated transcriptional changes of EAF1-AS1 within the heart atrial appendage. Similarly, effects associated with eQTLs located within the Brain Cerebellum were also recognized to confer major PD risk. These findings were replicated in two additional, independent cohorts (the UK Biobank, and NeuroX) and thus warrant further mechanistic investigations to determine if these transcriptional changes could act as early contributors to PD risk and disease development.

Defining early changes in Alzheimer's disease from RNA sequencing of brain regions differentially affected by pathology.

Tau pathology in Alzheimer's disease (AD) spreads in a predictable pattern that corresponds with disease symptoms and severity. At post-mortem there are cortical regions that range from mildly to severely affected by tau pathology and neuronal loss. A comparison of the molecular signatures of these differentially affected areas within cases and between cases and controls may allow the temporal modelling of disease progression. Here we used RNA sequencing to explore differential gene expression in the mildly affected primary visual cortex and moderately affected precuneus of ten age-, gender- and RNA quality-matched post-mortem brains from AD patients and healthy controls. The two regions in AD cases had similar transcriptomic signatures but there were broader abnormalities in the precuneus consistent with the greater tau load. Both regions were characterised by upregulation of immune-related genes such as those encoding triggering receptor expressed on myeloid cells 2 and membrane spanning 4-domains A6A and milder changes in insulin/IGF1 signalling. The precuneus in AD was also characterised by changes in vesicle secretion and downregulation of the interneuronal subtype marker, somatostatin. The 'early' AD transcriptome is characterised by perturbations in synaptic vesicle secretion on a background of neuroimmune dysfunction. In particular, the synaptic deficits that characterise AD may begin with the somatostatin division of inhibitory neurotransmission.

Misfolded proteins traffic from the endoplasmic reticulum (ER) due to ER export signals.

Most misfolded secretory proteins remain in the endoplasmic reticulum (ER) and are degraded by ER-associated degradation (ERAD). However, some misfolded proteins exit the ER and traffic to the Golgi before degradation. Using model misfolded substrates, with or without defined ER exit signals, we found misfolded proteins can depart the ER by continuing to exhibit the functional export signals present in the corresponding correctly folded proteins. Anterograde transport of misfolded proteins utilizes the same machinery responsible for exporting correctly folded proteins. Passive ER retention, in which misfolded proteins fail to exit the ER due to the absence of exit signals or the inability to functionally present them, likely contributes to the retention of nonnative proteins in the ER. Intriguingly, compromising ERAD resulted in increased anterograde trafficking of a misfolded protein with an ER exit signal, suggesting that ERAD and ER exit machinery can compete for binding of misfolded proteins. Disabling ERAD did not result in transport of an ERAD substrate lacking an export signal. This is an important distinction for those seeking possible therapeutic approaches involving inactivating ERAD in anticipation of exporting a partially active protein.

Not all hERG pore domain mutations have a severe phenotype: G584S has an inactivation gating defect with mild phenotype compared to G572S, which has a dominant negative trafficking defect and a severe phenotype.

INTRODUCTION: Mutations in the pore domain of the human ether-a-go-go-related gene (hERG) potassium channel are associated with higher risk of sudden death. However, in many kindreds clinical presentation is variable, making it hard to predict risk. We hypothesized that in vitro phenotyping of the intrinsic severity of individual mutations can assist with risk stratification. METHODS AND RESULTS: We analyzed 2 hERG pore domain mutations, G572S and G584S. Similar to 90% of hERG missense mutations, G572S-hERG subunits did not traffic to the plasma membrane but could coassemble with WT subunits and resulted in a dominant negative suppression of hERG current density. The G584S-hERG subunits traffic normally but have abnormal inactivation gating. Computer models of human ventricular myocyte action potentials (AP), incorporating Markov models of the hERG mutants, indicate that G572S-hERG channels would cause more severe AP prolongation than that seen with G584S-hERG channels. CONCLUSIONS: hERG-G572S and -G584S are 2 pore domain mutations that involve the same change in sidechain but have very different in vitro phenotypes; G572S causes a dominant negative trafficking defect, whereas G584S is the first hERG missense mutation where the cause of disease can be exclusively attributed to enhanced inactivation. The G572S mutation is intrinsically more severe than the G584S mutation, consistent with the overall clinical presentation in the 2 small kindreds studied here. Further investigation, involving a larger number of cohorts, to test the hypothesis that in vitro phenotyping of the intrinsic severity of a given mutation will assist with risk stratification is therefore warranted.

An HRD/DER-independent ER quality control mechanism involves Rsp5p-dependent ubiquitination and ER-Golgi transport.

We have identified a new pathway of ER-associated degradation in Saccharomyces cerevisiae that functions separately from the HRD/DER pathway comprised of Hrd1p, Hrd3p, Der1p, and Ubc7p. This pathway, termed Hrd1p independent-proteolysis (HIP), is capable of recognizing and degrading both lumenal (CPY* and PrA*), and integral membrane proteins (Sec61-2p) that misfold in the ER. CPY* overexpression likely saturates the HRD/DER pathway and activates the HIP pathway, so the slowed degradation kinetics of CPY* in a hrd1 Delta strain is restored to a wild-type rate when CPY* is overexpressed. Substrates of HIP require vesicular trafficking between the ER and Golgi apparatus before degradation by the ubiquitin-proteasome system. Ubiquitination of HIP substrates does not involve the HRD/DER pathway ubiquitin ligase Hrd1p, but instead uses another ubiquitin ligase, Rsp5p. HIP is regulated by the unfolded protein response as Ire1p is necessary for the degradation of CPY* when overexpressed, but not when CPY* is expressed at normal levels. Both the HIP and HRD/DER pathways contribute to the degradation of CPY*, and only by eliminating both is CPY* degradation completely blocked.

Publisher Correction: Enhancers active in dopamine neurons are a primary link between genetic variation and neuropsychiatric disease.

In the version of this article initially published, the legends for Supplementary Figs. 4-8 and 10-14 contained errors. The Supplementary Figure legends have been corrected in the HTML and PDF versions of the article.

Enhancers active in dopamine neurons are a primary link between genetic variation and neuropsychiatric disease.

Enhancers function as DNA logic gates and may control specialized functions of billions of neurons. Here we show a tailored program of noncoding genome elements active in situ in physiologically distinct dopamine neurons of the human brain. We found 71,022 transcribed noncoding elements, many of which were consistent with active enhancers and with regulatory mechanisms in zebrafish and mouse brains. Genetic variants associated with schizophrenia, addiction, and Parkinson's disease were enriched in these elements. Expression quantitative trait locus analysis revealed that Parkinson's disease-associated variants on chromosome 17q21 cis-regulate the expression of an enhancer RNA in dopamine neurons. This study shows that enhancers in dopamine neurons link genetic variation to neuropsychiatric traits.

ERADicate ER stress or die trying.

Stress within the endoplasmic reticulum (ER) induces a sophisticated network of pathways termed the unfolded protein response (UPR), which is mediated through the ER transmembrane sensors PERK, ATF6, and IRE1. The UPR coordinates the temporary downregulation of protein translation, the upregulation of ER chaperones and folding machinery, and the enhanced expression of components necessary for ER-associated degradation (ERAD) essential for decreasing ER stress by clearing terminally misfolded proteins from the ER. Repetitive but futile folding attempts not only prolong ER stress but can also result in reactive oxygen species (ROS) generation, both of which may result in cell death. Additional mechanisms for decreasing stress and the protein load in the ER have been recently revealed. They include a newly identified function of IRE1 in degradation of select secretory protein mRNAs, a "preemptive" quality control responsible for averting translocation of select secretory proteins into the ER, upregulation of forward trafficking to allow misfolded proteins with intact exit signals to exit the ER, and upregulation of autophagy. The saturation or failure of some or all of these mechanisms can result in cell death and disease, including diabetes and a number of late-onset neurologic diseases.