Skip to main content

Mitochondrial dysfunction is a key pathological driver of early stage Parkinson’s



The molecular drivers of early sporadic Parkinson’s disease (PD) remain unclear, and the presence of widespread end stage pathology in late disease masks the distinction between primary or causal disease-specific events and late secondary consequences in stressed or dying cells. However, early and mid-stage Parkinson’s brains (Braak stages 3 and 4) exhibit alpha-synuclein inclusions and neuronal loss along a regional gradient of severity, from unaffected-mild-moderate-severe. Here, we exploited this spatial pathological gradient to investigate the molecular drivers of sporadic PD.


We combined high precision tissue sampling with unbiased large-scale profiling of protein expression across 9 brain regions in Braak stage 3 and 4 PD brains, and controls, and verified these results using targeted proteomic and functional analyses.


We demonstrate that the spatio-temporal pathology gradient in early-mid PD brains is mirrored by a biochemical gradient of a changing proteome. Importantly, we identify two key events that occur early in the disease, prior to the occurrence of alpha-synuclein inclusions and neuronal loss: (i) a metabolic switch in the utilisation of energy substrates and energy production in the brain, and (ii) perturbation of the mitochondrial redox state. These changes may contribute to the regional vulnerability of developing alpha-synuclein pathology. Later in the disease, mitochondrial function is affected more severely, whilst mitochondrial metabolism, fatty acid oxidation, and mitochondrial respiration are affected across all brain regions.


Our study provides an in-depth regional profile of the proteome at different stages of PD, and highlights that mitochondrial dysfunction is detectable prior to neuronal loss, and alpha-synuclein fibril deposition, suggesting that mitochondrial dysfunction is one of the key drivers of early disease.


Neurodegenerative diseases are progressive, devastating and incurable, and are becoming increasingly prevalent in our ageing population [1]. Parkinson’s disease (PD) is the most common neurodegenerative movement disorder, affecting 1–2% of people aged over 65, with an estimated > 6.1 million people affected worldwide [2]. PD presents with clinical symptoms at an advanced stage of disease, when intraneuronal inclusions known as Lewy bodies (LB), composed primarily of insoluble aggregates of the protein alpha-synuclein, mitochondria, and lipid membranes, and significant neuronal loss are evident [3, 4]. Based on a cross-sectional study of post-mortem PD brain, Braak proposed that nigral and extranigral pathology do not occur simultaneously, but rather, LB pathology develops sequentially across the brain in anatomically connected regions in a caudal-rostral manner [5]. Whilst this progression of pathology may be controversial, several studies have confirmed that PD brains demonstrate alpha-synuclein positive lesions and neuronal loss in medullary, pontine and midbrain nuclei, with further pathology in the nucleus basalis of Meynert (90–98%), limbic cortex (50–60%), cingulate cortex (32–46%), frontal cortex (29–31%), and amygdala (25%), illustrating that the susceptible regions in PD and their interconnectivity remains consistent across studies [6, 7].

Genetic studies of familial and sporadic PD strongly suggest that altered expression, concentration, or mutant forms of alpha-synuclein are central to the pathogenesis of PD [8,9,10,11]. Furthermore, experimental and cellular studies suggest that alpha-synuclein may spread from neuron to neuron in a prion-like manner, inducing seeding and toxicity in neighbouring cells [12, 13]. Despite this, the mechanisms that cause toxicity, and lead to the pathological spread of disease throughout the brain, are still largely unknown. As no cell or animal model fully recapitulates the spatio-temporal spread of protein aggregates and neuronal loss of PD, such mechanisms can only be studied in PD patient brain tissue. Although, post-mortem brain is a powerful source of understanding the human disease, it is unable to distinguish disease-specific signals from the non-disease related ‘noise’ from stressed or dying cells in regions with advanced stage pathology. In regions with end-stage disease, there are large amounts of cell death, protein aggregation, chronic inflammation, oxidative stress and neurotoxicity and therefore many analyses of end stage tissue predominantly capture these pathways. Distinguishing whether these processes are causal, specific to Parkinson’s pathogenesis, and occur from the start of the condition, or whether they result from secondary consequences of dead or dying neurons is important in the search for the early events driving disease. Unbiased large-scale profiling of protein expression across thousands of proteins in specific brain regions, with sequential levels of pathology, therefore, has the potential to distinguish these signals, and uncover the mechanisms underlying brain region vulnerability and spread of pathology. However, proteomic studies investigating PD often focus on minimal brain regions (predominantly substantia nigra or frontal cortex) and/or focus on late stage disease where extensive neuronal loss can be observed [14,15,16,17,18,19,20]. Dixet et al [21]. reviewed the need for additional brain regions to be studied simultaneously to understand how protein expression changes with disease progression.

In this study, we address this gap in knowledge through the comprehensive investigation of the earliest changes in cerebral protein expression in PD. To achieve this, we selected cases with LB pathology corresponding to Braak stage 3 or 4, in which pathology is restricted to the brainstem, striatum and antero-medial temporal cortex, and not found in frontal and parietal cortices. We sampled across the 9 brain regions affected in PD, within any individual case, to capture the disease signature at mild, moderate, and severe disease stages, according to Braak staging. This approach enabled us to identify disease-specific convergent mechanisms underlying cell specific vulnerability. We utilised early-mid stage PD tissue to identify pathways common across disease regions, including those with minimal pathology, thus suggesting their involvement early in disease pathogenesis. We were able to identify proteins that exhibit altered expression prior to the appearance of pathology (unaffected regions), highlighting pathways that are dysregulated at the earliest stages of PD, and therefore likely driving the disorder. Together, the use of anatomically driven proteomic analyses, and the validation of targets through targeted proteomic and functional investigation, has identified critical disease mechanisms likely underpinning the pathogenesis of PD.

Materials and methods


Post-mortem brain tissue was obtained from the Neurological Tissue Bank, IDIBAPS-HC-Biobanc, Barcelona; Human Brain Tissue Bank, Budapest; UK Parkinson’s Disease Society Tissue Bank, Imperial College London; the London Neurodegenerative Diseases Brain Bank, Institute of Psychiatry, King’s College, London; Netherlands Brain Bank, Amsterdam; and the Newcastle Brain Tissue Resource. Informed consent was given in all cases. Cases and controls were matched as close as possible for age and sex and all had a post-mortem delay of less than 20 h. Details of cases used are listed in Table 1. Control subjects did not have a diagnosis of neurological disease in life and at post-mortem displayed age-related changes only. All cases were assessed for alpha-synuclein pathology and rated according to Braak staging [5]. PD cases were categorised as either Braak stage 3 or 4 (early) or Braak stage 6 (late). Ethical approval for the study was obtained from the Local Research Ethics Committee of the National Hospital for Neurology and Neurosurgery.

Table 1 Case demographics table

Dissection of regions of interest

Samples were taken where available (outlined in Table 1) from the Braak staging areas; substantia nigra, caudate nucleus, putamen, cingulate cortex, parahippocampal cortex, temporal neocortex, frontal cortex, parietal cortex, and as a control region, the cerebellum. To ensure samples contained only the region of interest and were not contaminated with white matter or surrounding nuclei, tissue was sampled from the appropriate frozen tissue blocks using a punch technique [22]. Accordingly, approximately 300–500 μm thick sections were cut in a cryostat at − 7 °C and for the removal of specific cortical regions and nuclei with the help of special hollow needles, ranging from 3 to 5 mm in diameter. Between 20 and 60 μg tissue was collected per region and samples were subsequently stored at − 80 °C until required [23].

Tissue preparation

Samples were prepared in two fractions a ‘soluble fraction’ and an ‘insoluble fraction’. Tissue was homogenised on ice in 500 μl of 50 mM Ambic buffer, 2% ASB-14 using the TissueRuptor mechanical homogeniser (Qiagen). After sonication, protein concentration was assessed using the bicinchoninic acid protein assay kit (Thermo Scientific). For initial LC–MS/MS analysis, a pool was created of equal amounts of PD cases or controls for each region to make a total pool of 1 mg of protein. At the same time, 300 μg of protein for each individual sample was collected for individual analyses. Pooled and individual samples then underwent acetone precipitation and the pellet fractions were treated with 70% formic acid as detailed previously [24].

In-solution digestion

Pooled freeze dried brain tissue samples were subjected to in-solution digestion as described previously [25, 26]. Briefly samples were resuspended and reduced using dithioerythritol then subsequently carboamindomethylated using iodoacetamide and digested using 1 μg of sequence grade Lys-C and trypsin (Promega).

High pH fractionation

Pooled peptide samples were fractionated at high pH prior to proteomic analysis to provide a greater depth of coverage of the proteome. 400 μl of 0.2% ammonia was added to 100 μl of the peptide solution, vortexed and spun at 13,000 g for 5 min. 50 mg Isolute C18 columns (Biotage) were primed with 1 ml of 50% acetonitrile containing 0.1% ammonia, followed by two 1 ml washes with 0.1% ammonia solution. The peptide mixture was next added to the column and allowed to flow-through under gravity. The initial breakthrough was collected, re-applied to the column and the subsequent breakthrough collected. 500 μl of 0.1% ammonia was used to wash the column and the peptides were subsequently eluted into 7 fractions with 500 μl of increasing concentrations of acetonitrile containing 0.1% ammonia; 5%, 8%, 10%, 15%, 25%, 100% and finally 100% methanol. All fractions were dried by centrifugal evaporation and then reconstituted in 3% acetonitrile, 0.1% trifluoroacetic acid. Prior to mass-spec analysis, peptides were spiked with 100fmol MassPREP™ Enolase (Waters UK).

Label-free quantitative MS

The first proteomic analysis was performed using MSe label free proteomics on a nanoAcquity nano-LC and QTOF Premier (Waters Corporation, Manchester, UK) as described previously [27] on pooled samples in 9 brain regions (substantia nigra, putamen, caudate nucleus, parahippocampal cortex, temporal cortex, cingulate cortex, frontal cortex, parietal cortex and cerebellum) in the first instance. Peptides were trapped, desalted, and separated as described previously [27].

The second proteomic analysis was performed on four Braak stage 3/4 cases and four controls as individual samples from the substantia nigra, putamen, parahippocampus and frontal cortex as previously described [28] using a nanoAquity coupled to a Synapt-G2-Si mass spectrometer with high definition ion mobility capability (Waters, UK) with online 2D fractionation. This enables deeper phenotyping than the first run using a high definition more sensitive mass spectrometer. Additionally, two Braak stage 6 samples were run from the putamen and frontal cortex. Samples were prepared and run using the same methods as above other than high pH fractionation procedures.

Data analysis of brain samples analysed by LC–MS/MS

First proteomic analysis: ProteinLynx GlobalServer version 2.4 (Waters Corporation) was used to process all data acquired from the first run as described previously [27]. Protein identifications were obtained by searching UniProt human reference proteome canonical database (June 2016): 1 missed cleavage, 4% false discovery rate and fixed modifications of carboamidomethylation of cysteines and dynamic modifications of oxidation of methionine.

Second proteomic analysis: Data acquired using the Synapt G2 MS were analysed using Progenesis LC–MS (Nonlinear Dynamics Limited, Newcastle, UK) raw data was processed as described previously but with the Uniprot database (June 2019) [24]. Briefly peptides were searched as described above except with 1% false discovery rate. Protein data for identifications with a confidence score > 20 and more than one unique peptide were exported for further analysis.


As a label-free approach was taken, all genes or proteins that met the threshold set (> twofold change in one region, and also changing in four other regions) were subject to gene ontology analysis using Webgestalt [29, 30] and DAVID (version 6.8, [31]). GOview was used to compare terms that were over-represented between regions or disease groups ( Ingenuity Pathway Analysis software (Qiagen) was used to perform in depth canonical pathway analysis and determine biological functions altered in the datasets.

Multivariate analysis on the second data analysis on the Synapt G2 MS was performed using SIMCA v15 (Umetrics, Sweden). The mitochondrial sub proteome analysis was done using Mitominer to annotate all proteins in the dataset as to mitochondria association by GO [32]. Mitochondrial proteins were grouped according to main mitochondrial functions using GO Panther [33].

Analysis of the cell type-specific protein expression was performed in a number of available datasets: mouse brain proteomics generated by Sharma et al. [34], human brain transcriptomics and proteomics [35, 36] and mouse PD transcriptomics [37] to assess for enrichment of the significantly associated proteins/genes in particular brain cell types. Enrichment analysis was performed using Bioconductor R package GeneOverlap [38] that returns the Fisher’s exact test p-value and odds ratio as previously performed in Wingo et al, [39]. Multiple comparison test was performed using Benjamini & Hochberg method and only overlaps with FDR < 5% were considered significant.

Targeted proteomic MRM LC–MS/MS assay

Potential candidates identified in the label free proteomics were selected for validation on the basis of  a change in expression over five of the nine brain regions with at least one region having a > twofold change. Secondly, proteins that had a > fivefold change in expression in any one region were added. Selected proteins were not included in the validation list eg. keratins, structural proteins, ribonucleoproteins. The assay was developed in the same way as shown in Heywood et al. [40]. A full list of peptides is given in Additional file 4: Table S1. Intact heavy labelled peptides were used as internal standards. Additionally, twenty nanograms of a generic yeast enolase whole protein standard (Sigma Aldrich, UK) was used as another level of internal standard. This was added to brain samples before they were processed from the individual protein extractions described above. Digested peptides were separated using a 50 mm C18+ UPLC column attached to a C18+ Vanguard and analysed using a Xevo TQ-S (Waters, Manchester, UK). UPLC and MS tune conditions were performed as previously described [40]. Samples were run in duplicate. Chromatograms were analyzed using Waters TargetLynx software. Peptides were standardized to an internal standard. Absolute levels were obtained from standard curves of custom synthesised peptides spiked into pooled brain digest.

Mitochondrial enzyme activity assays

Punch dissections from putamen, frontal cortex and parietal cortex were also from 4–6 of the controls and 4–6 of the Braak stage 3/4 PD cases. These cases were age/gender (male (M) and female (F)) and post-mortem delay (PMD) matched (See Table 1). Each brain region was homogenised and mitochondrial respiratory chain complex I, II–III, IV, II and citrate synthase activities were assayed spectrophotometrically at 30 °C as described in the study by Hargreaves et al. [41].


8 µm mounted FFPE tissue sections were incubated at 60 °C overnight. Sections were deparaffinised in Xylene and rehydrated in decreasing grades of alcohol. Slides were incubated in a methanol/hydrogen peroxide (0.3%) solution for 10 min to block endogenous peroxidase activity. Room temperature epitope retrieval and heat-induced epitope retrieval were utilised. Sections were first incubated in 98% formic acid for 10 min and then transferred to a boiling 0.1 M citrate buffer (pH 6.0) solution and pressure cooked at maximum pressure for 10 min. To block non-specific binding, sections were incubated in 10% non-fat milk in TBST for 30 min at room temperature. Mouse monoclonal alpha-synuclein [4D6] (abcam ab1903,1:10,000) primary antibody was diluted in TBST and sections were incubated in 200 µl of antibody solution for 1 h at room temperature. Sections were washed 3 times for 5 min in TBST and then were incubated for 30 min in 200 µl of biotinylated goat anti-mouse IgG secondary antibody (Vector Laboratories BA 9200, 1:200). Slides were washed as in the previous step and then incubated in pre-conjugated Strept(avidin)–Biotin Complex (ABC; DAKO) for signal amplification. Slides were then washed as in previous and were submerged in 5% 3,3ʹ-Diaminobenzidine (DAB) chromogen in TBST for 5 min then counterstained in Mayer’s haematoxylin (BDH). Finally, slides were dehydrated in increasing grades of alcohol (70, 90 and 100% IMS), cleared in xylene and mounted with DPX.


Detection of proteins in human post-mortem brain

Brains from 11 early-mid stage PD (Braak 3/4) cases were collected and specific regions were dissected according to tissue availability. Average age of death, average disease duration, and average post-mortem delay were calculated for all cases and matched to controls. The average post-mortem delay for cases was 7.92 h (range 3.5–16.5) and 9.21 h (range 3–23) for controls. Mean age of death for cases was 79 yrs (range 70–88) and 74 yrs (range 61–88) for controls. Average disease duration was 11.75 yrs (range 3–30). Mean pH in cases was 6.35 (range 6–6.67) and 6.30 in controls (range 5.51–6.63). Label free proteomic analysis was performed on 9 brain regions for PD and controls (pathological and clinical data outlined in Table 1). Brain regions consisted of substantia nigra, caudate nucleus, putamen, parahippocampal cortex, cingulate cortex, temporal cortex, frontal cortex, parietal cortex, and cerebellum. The degree of pathology was based on the Braak staging system and defined as ‘severe’ (substantia nigra), ‘moderate’ (caudate nucleus and putamen), ‘mild’ (parahippocampal cortex, cingulate cortex and temporal cortex), ‘unaffected’ (frontal cortex, parietal cortex and cerebellum). This is highlighted with representative pathology as shown in Fig. 1.

Fig. 1
figure 1

Case selection according to Braak staging system. a Diagram indicating the amount of Lewy body pathology in the regions used in the study according to Braak staging criteria. Severity of pathology at early-mid stage Braak 3/4 and late Braak stage 6 is highlighted by colours shown in reference to the key. Table indicates how many cases used for each region and disease group. b Panel of representative pathology images from a Braak stage 3/4 brain, highlighting the level of alpha-synuclein pathology in each region at this stage of disease at 20 × and 40 × magnification. Regions are labelled as follows: SN Substantia nigra; Cau Caudate; Put Putamen; Temp Temporal cortex; Para Parahippocampal gyrus; Cing Cingulate cortex; Fron Frontal cortex; Pari Parietal cortex. Black arrows show Lewy body presence and red arrows Lewy neurites. All scale bars represent 10 μm

In the first proteomic analysis, samples were pooled by region (workflow shown in Fig. 2a). For each pool, the optimised fractionation and solubilisation methods enabled the detection of a total of 1147 proteins: consisting of 1004 supernatant proteins (soluble) and 531 pellet proteins (insoluble) (Fig. 2b). The level of insoluble alpha-synuclein, determined by mass spectrometry, followed the same spatial gradient as the LB pathology using standard immunohistochemical methods in the Braak staging system; with the substantia nigra showing the highest expression, and the levels decreasing in moderate, mild, and unaffected regions, with the lowest expression found in the cerebellum (Fig. 2c). This confirms the use of the staging system to assess pathology severity.

Fig. 2
figure 2

Workflow of methodology. a Flowchart highlighting methods used within the experiment. Tissue was micro-dissected per region, tissue homogenised and proteins were extracted and digested into peptides. Samples were proteome profile using mass spectrometry. Candidates were selected before confirmatory targeted Multiple Reaction Monitoring LCMS/MS analysis was completed. Functional assays were performed. b Venn diagram indicating how many proteins were detected uniquely in the supernatant and pellet fractions and how many overlapped. c Heat map of alpha-synuclein expression as fold change in Braak stage 3/4 compared to control indicated by colour. Red indicated upregulation and green downregulation compared to controls with deeper colour indicating a higher level of change. White indicates no change in expression compared to controls. d Overall change in the brain proteome for each region and condition expressed as a percentage of total proteome either upregulated (red) or downregulated (green) by more than 1.5 fold. S represents the supernatant fraction and P represents the pellet fraction. Graph created in GraphPad Prism v8

We calculated the % proteome change for each region in the soluble and insoluble fractions and found that the severely affected region (substantia nigra) was associated with the greatest proteome change, whilst moderate and mild regions had fewer changes in protein expression (Fig. 2d). Of note, the insoluble fraction had consistently more changes in protein levels than the soluble fraction across all regions.

Bioinformatic analysis highlights mitochondrial dysfunction across the PD brain

The fold change for all proteins in the PD brain (compared to controls) was analysed using IPA (see methods). The top canonical pathway, that was most significantly represented throughout the data across all regions, was mitochondrial dysfunction (see Fig. 3a). When observing the protein expression as fold change compared to control for the mitochondrial pathway, many of the proteins were upregulated (depicted in red), or down regulated (depicted in green; Additional file 1: Figure S1b).

Fig. 3
figure 3

Pathway analysis of proteomics. Heatmap from IPA software showing the most significantly altered canonical pathways found in a the first mass spectrometry dataset and b the second mass spectrometry dataset with Mitochondrial dysfunction the top pathway significantly altered across all regions in Braak stage 3/4 compared to controls and in both datasets. The p-value is indicated by purple colour. The deeper the colour the lower the p-value

Proteins that were changing in expression at > 1.5 fold (Table 2) were further assessed using Webgestalt [29, 30] and David v6.8 [31] to generate a cumulative list of gene ontology (GO) terms for each region (any duplicate terms across the software programs were excluded) in the sections Biological Process (BP), Molecular Function (MF) and Cellular Components (CC). In the substantia nigra (severe), there were a larger number of GO terms (276 BPs/135 MF/144 CCs) compared to mildly affected or unaffected regions (parietal region had 132 BPs/41MF/80 CCs), demonstrating that the degree of altered protein expression reflects the pathological severity of the region in PD brain (Additional file 1: Figure S1a). The full list of BPs/MF’s and CCs found across all regions are shown in Additional file 5: Tables S2. Additional file 6: Tables S3, Additional file 7: Tables S4 respectively.

Table 2 Differently expressed proteins across brain regions in early PD (Braak stage 3/4) compared to controls

Next, we were interested in common pathways that are affected across multiple regions, and therefore ranked pathways according to their detection in different regions. 17 GO terms for BPs were affected in all regions; 9 GO terms for MF in all regions; 16 GO terms for CC in all regions. These common overlapping pathways in multiple regions (5 or more) included: translation, protein folding and assembly, metabolic and energy dependent processes (for BPs); activity of antioxidants, metabolism, and small molecule binding (for MFs). In the cellular compartment analysis, the top three compartments in all regions included extracellular space, cytoplasm, and mitochondria.

Validation of proteome hits

To validate the hits generated from the list of differentially expressed candidate proteins (Table 3), we developed and performed targeted MRM LC–MS/MS analysis on 86 proteins across all early-mid PD cases (Braak stage 3/4), and controls, and brain regions as shown in Table 1. From this analysis we were able to quantify the expression fold changes of 73 proteins. Of these, 57 proteins were statistically significantly elevated in at least one region, 17 of which have mitochondrial localisation. Notably, the largest number of mitochondrial proteins (13) showed statistically significant changes in the parahippocampus, a region that exhibits only mild pathology. These are shown in Fig. 4, and include NDUA2 (4.12 fold increase, p = 0.033), M2OM (4.73 fold increase, p = 0.01), ACON (1.8 fold increase, p = 0.015), VDAC2 (2.64 fold increase, p = 0.004), ODO2 (2.04 fold increase, p = 0.006), IDH3A (twofold increase, p = 0.011), NU5M (1.53 fold increase, p = 0.009), GRP75 (4.18 fold increase, p = 0.03), SCOT1 (3.6 fold increase, p = 0.01), ADT1 (0.38 fold increase, p = 0.04), ATPK (2.62 fold increase, p = 0.035), ATP5E (4.36 fold increase, p = 0.0067), EHCM (2.15 fold increase, p = 0.007). Other proteins changing in the parahippocampus, not related to the mitochondria, included UCLH1, TDP43, cold inducible RNA binding protein, ATP13A2, AATC. Changes in a mildly affected region, at the boundary of pathology vs no pathology, reflect alterations in pathways that occur early in the disease process, here indicated to be mitochondrial metabolism, bioenergetics, ATP synthesis, and the mitochondrial stress response.

Table 3 Candidate list of proteins to validate
Fig. 4
figure 4

Validation of increased mitochondrial proteins in PD cases using multiple reaction monitoring LC–MS/MS. Graphs highlighting differences in fold change as a ratio against internal standard an for multiple mitochondrial proteins in the parahippocampus, o metallothionein-2 across the putamen, parahippocampus, cingulate cortex and frontal cortex, and p mitochondrial 2-oxoglutarate malate carrier across the putamen, parahippocampus and cingulate cortex. Regions are colour-coded for severity of pathology present at that stage of disease as shown in key. T-tests were done for each pairing and statistically significant results at p < 0.05 are highlighted on the graphs. Graphs and statistics completed with GraphPad Prism v8

Certain proteins demonstrated changes across several brain regions with varying pathology. As shown in Fig. 4, metallothionein-2 showed an increase in regions with moderate pathology (putamen, p = 0.03), mild pathology (cingulate cortex, p = 0.022; parahippocampus, p = 011) and no pathology (frontal cortex, p = 0.008). Mitochondrial 2-oxoglutarate malate carrier exhibited a fivefold increase in a region with moderate pathology (putamen, p = 0.03), a 4.7 fold increase in a region with mild pathology (parahippocampus, p = 0.01) and a 1.8 fold increase in a second region with mild pathology (cingulate cortex, p = 0.01), Fig. 4p. Microsomal glutathione S transferase demonstrated changes in a moderate pathology region, putamen (2.93, p = 0.04), as well in the mildly affected regions, cingulate cortex (0.99, p = 0.006), and temporal cortex (1.54, p = 0.0025). Such proteins suggest pathways affecting multiple regions susceptible to PD, independent of disease stage, may include mitochondrial metabolism, redox balance, and metal ion homeostasis.

Proteins that showed statistically significant changes in the frontal and parietal cortex (no pathology at this early-mid stage, but predicted to be vulnerable in PD), may highlight the very earliest changes in cells seen before the emergence of a synucleinopathy. These include Sirtuin-2, metallothionein 2, and glycogen phosphorylase.

Proteins that were highlighted to change (HADHB, ACAT1, VDAC2, IDH3A, GPR75, SIRT2, G6PD, ATP13A2, NDUFA2, SLC25A11, ACO2, DLST, OXCT1, SLC25A4, NDUFS5, ATP5MF, ATP5F1E, ECHS1, MT2A, TARDBP, GOT1, MGST1, PYGL) were mined in the available cell specific datasets, and those detected are shown in Additional file 8: Table S5. Protein enrichment for the selected targets (23) were analysed using the GeneOverlap Bioconductor package in cell types and a significant association was found with Sharma et al (2015) oligodendrocyte dataset (FDR = 0.03) with overlap occurring for 3 target proteins: SIRT2, ACAT1 and PYGL. Other proteins overlapped with additional other cell types e.g MT2A overlapped with astrocyte cell type data (Additional file 8: Table S5) but these were non-significant. These results indicate that further investigation to determine how these key proteins act within different cell types is needed.

Detailed characterisation of the mitochondrial proteome in 4 stages of disease

We sought to understand in greater depth, the nature of the mitochondrial dysfunction in early-mid post-mortem brain. Therefore, we repeated the HPLC–ESI–MS/MS analysis on the higher resolution Synapt-G2 Si machine, using a sub-selection of the original samples, performed individually across 4 selected regions: substantia nigra, putamen, parahippocampus and frontal cortex (representing severe, moderate, mild and unaffected regions). We detected a combined (supernatant and pellet) total of 2343 proteins. A principal components analysis (Fig. 5) shows that the variation between samples is primarily accounted for by brain region. Furthermore, the differences in the proteome between control and early PD are modest, suggesting that early disease varies subtly from control brain.

Fig. 5
figure 5

Principle component plot of Synapt G2-Si data. Variation is accounted primarily by brain region. Sub region variation can be observed between late PD and control and early PD but no overall proteomic variation is observed between control and early PD indicating the changes are subtle. Each brain region is mapped on hemi-dissected hemisphere and coronal slice, in its respective annotated colour

In the first analysis (pooled samples, 8 PD brain regions), we illustrated that 202 proteins with altered expression (> 1.5 fold) in early PD vs control were associated with mitochondria (human MitoCarta 2.0 [42]). In the second analysis (individual samples, 4 brain regions), 140 mitochondrial proteins exhibited a change in expression (> 1.5 fold) in early PD vs controls (79 of these overlapped with the proteins identified in the first analysis). The second data set (2–4 individuals, 4 brain regions) was analysed using the IPA software. The three most significantly changing pathways across the 4 regions were ‘mitochondrial dysfunction’, ‘oxidative phosphorylation’, and ‘sirtuin signalling’, thereby replicating our original finding (Fig. 3b). The altered expression of mitochondrial proteins in each region are shown in Table 4 and as a heatmap in Additional file 1: Figure S1c.

Table 4 Mitochondrial proteins that show altered expression

The mitochondrial proteome was further analysed to assess the specific mitochondrial processes affected. Mitochondrial proteins were sub-grouped into key mitochondrial functions and the percentage of significantly affected proteins in each function is illustrated in Fig. 6c and g. Regions with mild and severe pathology (late frontal cortex, and putamen of early/late PD), demonstrated changes in the ‘pyruvate dehydrogenase complex’, ‘TCA cycle’, ‘Fatty acid catabolism’, and ‘stress response’. Interestingly, regions with no pathology (early frontal) still exhibited alterations in the TCA cycle. Fatty acid oxidation is not the primary source of energy in the brain, however, neuronal cells utilise fatty acids as fuel when in a diseased state [43]. These changes may indicate altered metabolic flux of energy substrates, such as acetyl-coA. Fatty acid oxidation is associated with increased expression of ACAT1, which is increased in the early frontal cortex as well as the late frontal cortex. This indicates that there may be a switch in the metabolic use towards fatty acid oxidation. Notably, we also observe decreased expression of mitochondrial trifunctional enzyme subunit beta (HADHB), in early and late PD putamen, compared to controls (Fig. 6d).

Fig. 6
figure 6figure 6

Mitochondrial analysis of label free proteomic data. Oxidative phosphorylation pathway diagram with per region overlay of proteomic expression for the putamen and frontal cortex in early PD compared to controls, a and late PD compared to controls, e. Red represents upregulation and green represents downregulation. Intensity of colour indicates level of altered expression. b, f Complex proteins totalled and ratioed to citrate synthase across brain regions and disease stages. Two-way ANOVA with Sidek’s multiple comparison post-hoc testing was performed for each complex/citrate synthase and statistically significant results are demonstrated on graphs with * representing p < 0.03, ** representing p < 0.002, *** representing p < 0.0002 and **** representing p < 0.0001. Bars represent SEM. c, g Table indicating the % of significantly altered proteins per pathway for the sub-functions of the mitochondria as determined through Mitominer and Panther. Regions are colour-coded for severity of pathology present at that stage of disease as shown in key. d HADHB and ACAT1 expression across multiple brain regions, Graphs and statistics were completed using GraphPad Prism v8

Mitochondrial complex protein abundance values were extracted, totalled, and ratioed to protein abundance of citrate synthase to compare between conditions (Fig. 6b and f). Notably, an upregulation of complex I, II and V occurs in early putamen. Changes in complex II and/or IV appear to only occur in late putamen and frontal cortex. To determine if changes could be due to variation in mitochondrial number, the levels of citrate synthase were examined. Citrate synthase levels did not change between disease groups (p = 0.1647) but did alter between regions (p < 0.0001, Additional file 3: Figure S3).

Overall, this data demonstrates that multiple mitochondrial pathways may be simultaneously affected in PD, and these include the TCA cycle and fatty acid oxidation, as well as oxidative phosphorylation and ATP synthesis.

Mitochondrial pathology worsens with disease progression

In addition to assessing Braak stage 3/4 brains, we assessed late Braak stage 6 PD brain (compared to controls) for two selected regions—the frontal cortex and the putamen, shown in the PCA plot (Fig. 5). This enabled us to study the same region across two disease stages: that is, the frontal cortex unaffected (Braak stage 3, termed early frontal) and mildly affected (Braak stage 6, termed late frontal). We studied the overlap between the early brain changes (early vs control) and the late brain changes (late vs control) (Fig. 7). Together, this analysis demonstrated that the progression to late stage brain, within a region, is associated with common pathways to early PD brain, and the acquisition of more, and new pathways in late PD.

Fig. 7
figure 7

Early stage PD overlap against Late stage PD. Pie charts show the number of GO terms (Biological processes, molecular functions or cellular components) that are overlapping between early PD (Braak stage 3/4) and late stage PD (Braak stage 6) in either the putamen, a region moderately affected at Braak stage 3/4 but severely affected at Braak stage 6, and the frontal cortex, a region not affected at Braak stage 3/4 but mildly affected at Braak stage 6. Colour coded for severity of pathology present at that stage of disease as shown in key

Exclusive to early frontal cortex where no pathology is observed, are ‘mitochondrion organisation’, ‘protein localisation to mitochondrion’ and ‘reactive oxygen species metabolism’ (Additional file 10: Table S7). Overlapping between early frontal cortex and late frontal cortex includes ‘response to reactive oxygen species’, ‘mitochondrial electron transport’, ‘mitochondrial complex 1’, and ‘mitochondrial apoptosis’. Exclusive to late frontal cortex includes ‘mitochondrial fission’, ‘release of cytochrome c from mitochondria’, ‘acetyl-coA and pyruvate metabolism’, ‘TCA cycle’, and ‘mitochondrial electron transport’.

The overlapping BPs in early putamen (Braak stage 3) vs late putamen (Braak stage 6) included ‘response to reactive oxygen species’, ‘fatty acid beta-oxidation’ and ‘insulin receptor signalling pathway’ (Additional file 11: Table S8). Exclusive to the late stage putamen includes ‘mitochondrial genome maintenance’, ‘mitochondrial fission’, ‘acetyl-CoA metabolism’, ‘pyruvate metabolism’, ‘TCA cycle’, and ‘mitochondrial electron transport’. This suggests the existence of mitochondrial pathway disruption early in putamen pathology, with progression to prominent and persistent mitochondrial dysfunction in late putamen pathology.

Together, this data suggests that mitochondrial pathways are disrupted in unaffected regions in early PD, and that these processes progress to affect further mitochondrial and metabolic pathways over the course of the disease.

Validation of mitochondrial activity in PD brain

Due to the high prevalence of pathways involving mitochondrial metabolism and ATP synthesis, we reasoned that mitochondrial dysfunction may be detectable in early PD brain. Respiratory chain activity can be studied in post-mortem brain, and has been previously implicated in PD pathology, primarily relating to complex I dysfunction [44, 45]. We took four brain regions reflecting moderate pathology, mild pathology, and no pathology (putamen, parahippocampal cortex, frontal cortex, and parietal cortex), and investigated their respiratory chain activity, normalising to citrate synthase (Fig. 8). In the moderately affected region in this study (putamen), we found a significant decrease in complex I activity (p = 0.026), complex II plus III activity (p = 0.0043), and in complex IV activity (p = 0.0281). Interestingly, this reveals complete impairment of the respiratory chain. Notably, this decrease in complex activity was associated with an increase in complex protein levels (Fig. 6b and f), perhaps indicative of a compensatory mechanism or that the increased complex protein is damaged and non-functional. In the mildly affected cortical region (parahippocampal cortex), we found a significant reduction in complex I function only (p = 0.0079), with preservation of the other complexes. However, in unaffected regions (frontal and parietal cortex) we found no impairment of respiratory chain activity, whilst isolated complex II activity was unaffected in all cases. This data demonstrates that in early PD brain, the appearance of respiratory chain dysfunction and metabolism occurs in association with mild aggregation pathology and is associated with a progressive and severe respiratory chain dysfunction in mildly and moderately affected regions.

Fig. 8
figure 8

Mitochondrial assays show altered mitochondrial function in early affected regions. Functional mitochondrial assays across the putamen, parahippocampus, frontal cortex and parietal cortex for a Complex I, b Complex II, c Complex II/III and d Complex IV. All are normalised to citrate synthase levels and the ArcSin (ASIN) of data used for statistical analysis. Regions are colour-coded for severity of pathology present at that stage of disease as shown in key. Mann–Whitney U-tests were performed with * representing p < 0.03, ** representing p < 0.002, *** representing p < 0.0002 and **** representing p < 0.0001. Error bars represent SEM. Graphs and statistics were completed using GraphPad Prism v8


This study identifies critical proteins and pathways altered in early sporadic PD. Importantly, we focus on Braak stage 3–4 brain, in which LB pathology and neuronal loss are found on a gradient from severe-moderate-mild-unaffected [5]. Sampling across this gradient allows the capture of the early molecular changes as well as the late molecular changes resulting from both PD pathology and the bystander effects of neuronal loss. Here, we performed an unbiased capture of data across the PD brain, harnessing the spatio-temporal gradient of pathology, to investigate 2 key disease signatures: (i) common proteins and pathways across disease regions independent of pathological stage, thus reflecting shared causal mechanisms and/or shared molecular signatures of cell vulnerability, (ii) proteins and pathways with altered expression in regions with mild or unaffected pathology, likely reflecting the early drivers of the disease.

Proteomics provides large scale data of protein expression in complex biological samples, enabling the inference of biologically meaningful functional information. However, tissue heterogeneity and technical variability has led to a lack of consistency and reproducibility [46]. Here, punch dissection of carefully delineated anatomical regions enabled precision of sampling to the level of cortical areas and brain nuclei, whilst still capturing multiple cell types in a region [22] and retaining sufficient throughput. Our proteomic data confirmed that the expression of insoluble alpha-synuclein in these brains existed in the same biochemical gradient as that characterised immunohistochemically by LB pathology. We demonstrate that the number of altered biological processes or molecular function terms matches the severity of pathology. As expected, the highest number of altered pathways occurs in severely affected regions, with the smallest number found in mild or unaffected regions. The spatio-temporal gradient of pathology in early PD brains is thus matched by a biochemical gradient of a changing proteome. This observation further supports the approach of studying the small number of altered pathways from pathologically unaffected or mildly affected regions to understand the driving events in the disease process.

Identifying common pathways across early PD brain

Ingenuity pathway analysis revealed mitochondrial dysfunction as the ‘top pathway’ in early PD, with altered expression across all regions. The comparison of common pathways across the early PD brain, filtering for proteins with a fold change > 1.5 and common to > 5 brain regions, resulted in the clear emergence of mitochondrial proteins (33% of total) and mitochondrial pathways such as metabolism, oxidative phosphorylation, and oxidative stress. Next, we compared two specific brain regions at different stages of disease: (i) the frontal cortex of Braak stage 3 (unaffected) versus Braak stage 6 (mildly affected), in which the convergent processes in the same region with different stages of pathology included mitochondrial dysfunction. (ii) the putamen from Braak stage 3 (moderate pathology) versus Braak stage 6 (severe pathology) also demonstrated mitochondrial metabolism and oxidative stress as convergent pathways. These approaches, that is, within early PD brain across pathology stages, as well as across PD brain of different stages, highlight the importance of mitochondrial function early in the disease, prior to pathology, and in late disease.

Specific mitochondrial located proteins that demonstrated clearly altered expression in a single or multiple affected regions of the early stage brain in both proteomic studies and the validation study predict critical processes that are affected in early PD:

  1. (i)

    Oxidative phosphorylation and ATP synthesis, due to specific changes in components of the respiratory chain, notably NADH-ubiquinone oxidoreductase chain 5 (Core subunit of the mitochondrial membrane respiratory chain NADH dehydrogenase (Complex I)), and NADH dehydrogenase (ubiquinone) 1 alpha subcomplex subunit 2 (Accessory subunit of the mitochondrial membrane respiratory chain NADH dehydrogenase (Complex I)), ATP synthase F(0) complex subunit B1, and ATP synthase subunit epsilon (both key components of mitochondrial membrane ATP synthase (F1F0 ATP synthase or Complex V) that produces ATP from ADP).

  2. (ii)

    Mitochondrial metabolism, due to specific changes in Isocitrate dehydrogenase (NAD) subunit alpha (Catalytic subunit of the enzyme which catalyses the decarboxylation of isocitrate (ICT) into alphaketoglutarate); Succinyl-CoA:3-ketoacid coenzyme A transferase 1 (Key enzyme for ketone body catabolism which transfers the CoA moiety from succinate to acetoacetate), Dihydrolipoyllysine-residue succinyltransferase component of 2-oxoglutarate dehydrogenase complex (The 2-oxoglutarate dehydrogenase complex catalyses the overall conversion of 2-oxoglutarate to succinyl-CoA and CO2), Aconitate hydratase (Catalyses the isomerisation of citrate to isocitrate via cis-aconitate), Mitochondrial 2-oxoglutarate malate carrier (Catalyses the transport of 2-oxoglutarate across the inner mitochondrial membrane in an electroneutral exchange for malate or other dicarboxylic acids). Other affected mitochondrial proteins include the Voltage dependent anion-selective channel protein 2, which mediates transport across the mitochondrial outer membrane, and Stress-70 protein, a mitochondrial chaperone protein, and ADP/ATP translocase which mediates the transport of ATP out of the mitochondria. Interestingly, Microsomal glutathione S-transferase 3, changed its expression across 3 PD regions, and is the enzyme that catalyses oxidation of hydroxy-fatty acids, and participates in lipid metabolism.

  3. (iii)

    Glutathione metabolism and redox regulation, glucose-6-phosphate dehydrogenase (G6PD) is reduced in early frontal cortex. G6PD is involved in the pentose-phosphate pathway and reduces glutathione [47]. This is in line with previous studies that have commented on the pentose-phosphate pathway involvement in PD and identified reduction in G6PD in the putamen of PD post-mortem brain [48]. Reduced Glutathione combats and protects against oxidative stress [49, 50]. The enzyme that catalyses this reaction, glutathione reductase is also reduced in the substantia nigra and putamen. Glutathione-S-transferase omega-1 is also reduced in early putamen and early frontal cortex. This enzyme is responsible along with glutathione for protecting the mitochondrial membranes from oxidative damage [51]. Alteration of this pathway early in disease would explain why an increase in oxidative stress is observed. As a consequence of increased oxidative stress, protein mixed disulfides can form and this can cause a number of adverse effects for other proteins and the way they function [52]. Additionally oxidative stress can reduce transcription of TCA cycle enzymes [53]. Hence the disruption of these pathways in early PD may lead to increased mitochondrial dysfunction and progression of disease.

Notably, respiratory chain enzyme activity can be assessed in frozen post-mortem tissue, with this approach leading to the discovery that impaired complex I function is associated with PD pathology [44, 45]. Here, we discovered that regions of moderate pathology, including the putamen, are associated with severe and comprehensive functional abnormality, affecting all complexes of the respiratory chain. In areas of mild pathology, such as the parahippocampus, in which significant changes of mitochondrial proteins, including complex I were validated, only complex I is significantly impaired. Conversely, in pathologically unaffected regions, which were associated with minimal changes in mitochondrial located proteins, no functional deficit was observed. Therefore, we concluded that in early stage PD brain, the mitochondrial respiratory chain function follows a gradient of activity from no impairment, to severely impaired, along the gradient of α-synuclein pathology. Interestingly, the primary induction of complex I dysfunction in mice results in reduced dopamine synthesis and motor deficits [54], further highlighting the likely importance of mitochondrial bioenergetic impairment early in the disease. It should be noted that no loss of isolated complex II activity was demonstrated in this study. In areas with reduced complex activity, we observed an increase in expression levels of complexes. It is therefore possible that the presence of oxidative stress leads to post translational modifications that damage the protein complexes, altering their function. Further work is required to ascertain why this exclusively nuclear encoded component of the respiratory chain is preserved, i.e. the susceptibility of nuclear vs mitochondrial DNA and their gene products, to oxidative stress [55].

Nigral alterations in mitochondrial respiration and oxidative damage have been previously identified in post-mortem brain [56, 57]. Analysis by Werner et al [58]. of the substantia nigra, of 5 PD cases, detected 321 differentially expressed proteins, including those involved in iron metabolism and glutathione related redox metabolism. Another study of the substantia nigra found 119 differentially expressed proteins in PD patients, including proteins involved in mitochondrial dysfunction, oxidative stress, protein degradation, and neuroinflammation [15]. Notably, oxidative damage is not limited to the substantia nigra in post-mortem PD brain, but is evident in the frontal cortex, with proteins such as alpha-synuclein, UCLH1, SOD1, SOD2, and DJ-1 undergoing oxidative modifications [59, 60]. Furthermore, early stage PD has been associated with oxidatively modified proteins implicated in glycolysis and energy metabolism [61]. Therefore, the current literature suggests that mitochondrial dysfunction and oxidative stress are critical features of PD pathology.

Our data confirms that alterations in mitochondrial pathways and function are key features of the PD brain. Importantly, we find that significant changes in mitochondrial protein expression occur early in the disease process in pathologically unaffected regions (frontal and parietal cortex), indicating the likely role of mitochondrial dysfunction in PD pathogenesis. We next examined further changes in unaffected regions that may occur prior to the development of alpha-synuclein pathology. Our study validated metallothionein-2 and sirtuin-2 as 2 proteins that change significantly in regions unaffected pathologically, but predicted to be affected in PD. Metallothionein-2 has a high content of cysteine residues that bind various heavy metals, altering the bioavailability of metals associated with PD, such as iron, copper, and zinc [62, 63]. Interestingly, both cell and animal models have demonstrated that metallothionein-2 protects against the nitrative and peroxynitrite stress induced by alpha-synuclein [64, 65], whilst its treatment prevents paraquat-induced loss of dopaminergic neurons [66]. We show an overlap of MT2A in astrocyte proteomic cell type data. It has been previously shown to be expressed in astrocytes and further to be transferred to neurons during the process of axon regeneration [67, 68]. Here, metallothionein-2 exhibited increased expression along the gradient of pathology in the early stage PD brain and remained elevated in unaffected frontal cortex. This highlights a protein involved in heavy metals and oxidative stress as playing a key role in the early stages of PD and indicates the potential for astrocytic involvement in early disease processes.

NAD-dependent protein deacetylase sirtuin-2 is considered a master regulator of ageing processes, modulating multiple and diverse biological processes through the deacetylation of internal lysines on many proteins. Sirtuin-2 is implicated in a number of pathways critical in PD pathogenesis, including mitochondrial metabolism, mitophagy, autophagy, and oxidative stress [69,70,71,72,73,74]. Furthermore, sirtuin-2 has been shown to deacetylate alpha-synuclein, driving its aggregation and toxicity [75]. Thus, the inhibition or deletion of sirtuin-2 is shown to be neuroprotective in PD models [76,77,78]. However, reduced sirtuin-2 has also been shown to have deleterious effects [73, 79, 80], demonstrating that whilst its loss may reduce alpha-synuclein aggregation, it may also dysregulate other key pathways, leading to reduced ATP levels and increased oxidative stress.

Given its role in key PD pathways, it is of major interest to see that sirtuin-2, in the proteomic study, was reduced across six affected PD regions, and validated in 3 regions. Importantly, sirtuin-2 is reduced in the frontal cortex (unaffected region) and temporal cortex (mildly affected region) of early stage PD brain. Interestingly, Sirtuin-2 is predominantly expressed in oligodendroglia, and is critical to their differentiation and myelination [81]. Here, we find SIRT2, ACAT1 and PYGL to be significantly enriched in oligodendrocyte cell-type specific proteomic data. Oligodendrocyte-specific proteins, such as the myelinating proteins MBP, MOG, MOBP, PLP, OMG, MPM2, and CNP, are also downregulated across disease stages, including in unaffected regions by > 1.5 fold compared to controls. This alteration in oligodendrocyte-specific proteins in ‘yet to be affected’ regions in PD brains, supports recent transcriptomic studies implicating oligodendrocyte dysfunction as an early and important event in PD [82,83,84]. The reduction in sirtuin-2 in early PD may therefore (i) be linked to oligodendrocyte dysfunction and represent an early mechanism by which mitochondrial metabolism and oxidative stress occurs, or (ii) may in some way be linked to a protective mechanism against, and prior to the emergence of, alpha-synuclein aggregation.


In summary, we provide a detailed and comprehensive assessment of the human PD brain proteome across disease stages to capture the early and common signatures of cells vulnerable to PD pathology. Here, we demonstrate the presence of mitochondrial dysfunction throughout the PD brain, irrespective of the level of pathology (Fig. 9). Critically, we show that this dysfunction occurs prior to, or concomitantly with, the earliest stages of alpha-synuclein aggregation, and certainly prior to neuronal loss, implicating mitochondrial dysfunction as an early driver of the disease, and not just an end-stage phenomenon. Furthermore, we highlight two proteins, sirtuin-2 and metallothionein 2, that significantly change in vulnerable PD regions prior to the appearance of alpha-synuclein pathology, and thus require further investigation as to their pathologic or compensatory role in PD. Together, our use of comparative brain proteomics in carefully selected regions from early stage disease reveals the pathways driving the pathogenesis of sporadic PD. Importantly, our identification of these pathways, rather than those implicated later in the disease, may highlight novel targets for therapeutic intervention that could slow the underlying progression of PD.

Fig. 9
figure 9

Progression of mitochondrial protein expression changes through the early PD brain. Schematic showing how mitochondrial proteins change in expression through the areas of the brain that are affected by alpha-synuclein pathology in early stage PD. Green represents downregulated proteins and red represents upregulated proteins when compared to controls. Arrows represent severity of the pathology across regions of the brain. Middle brain image represents how the pathology has spread through the brain at Braak stage 3/4 (adapted from [85]) Mitochondrial expression diagrams were made using IPA software. Figure made with

Availability of data and materials

The datasets generated and/or analysed during the current study are available from the corresponding author on reasonable request.



Parkinson’s disease


Lewy bodies


Liquid chromatography with tandem mass spectrometry


Gene ontology


Multiple reaction monitoring


Ultra performance liquid chromatography


Post-mortem delay


Formalin-fixed paraffin-embedded


Tris buffered saline with tween




Industrial methylated spirit


Ingenuity pathway analysis


Biological processes


Molecular function


Cellular component


High performance liquid chromatography coupled with electrospray ionisation with tandem mass spectrometry


Tricarboxylic acid


Principle components analysis


Adenosine triphosphate


Adenosine diphosphate


Nicotinamide adenine dinucleotide


Glucose-6-phosphate dehydrogenase


  1. Hou Y, Dan X, Babbar M et al (2019) Ageing as a risk factor for neurodegenerative disease. Nat Rev Neurol 15:565–581.

    Article  PubMed  Google Scholar 

  2. GBD (2015) Neurological disorders collaborator group (2017) Global, regional, and national burden of neurological disorders during 1990–2015: a systematic analysis for the Global Burden of Disease Study 2015. Lancet Neurol, 16:877–897,

  3. McCann H, Stevens CH, Cartwright H, Halliday GM (2014) α-Synucleinopathy phenotypes. Parkinsonism Relat Disord 20(Suppl 1):S62–S67.

    Article  PubMed  Google Scholar 

  4. Shahmoradian SH, Lewis AJ, Genoud C et al (2019) Lewy pathology in Parkinson’s disease consists of crowded organelles and lipid membranes. Nat Neurosci 22:1099–1109.

    Article  CAS  PubMed  Google Scholar 

  5. Braak H, Tredici KD, Rüb U et al (2003) Staging of brain pathology related to sporadic Parkinson’s disease. Neurobiol Aging 24:197–211.

    Article  PubMed  Google Scholar 

  6. Ulusoy A, Rusconi R, Pérez-Revuelta BI et al (2013) Caudo-rostral brain spreading of α-synuclein through vagal connections. EMBO Mol Med 5:1119–1127.

    Article  CAS  PubMed  Google Scholar 

  7. Jellinger KA (2003) Alpha-synuclein pathology in Parkinson’s and Alzheimer’s disease brain: incidence and topographic distribution–a pilot study. Acta Neuropathol 106:191–201.

    Article  PubMed  Google Scholar 

  8. Sahay S, Ghosh D, Dwivedi S et al (2015) Familial parkinson disease-associated mutations alter the site-specific microenvironment and dynamics of α-synuclein. J Biol Chem 290:7804–7822.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Wood SJ, Wypych J, Steavenson S et al (1999) alpha-synuclein fibrillogenesis is nucleation-dependent. Implications for the pathogenesis of Parkinson’s disease. J Biol Chem 274:19509–19512.

    Article  CAS  PubMed  Google Scholar 

  10. Conway KA, Harper JD, Lansbury PT (1998) Accelerated in vitro fibril formation by a mutant alpha-synuclein linked to early-onset Parkinson disease. Nat Med 4:1318–1320.

    Article  CAS  PubMed  Google Scholar 

  11. Kim C, Lee S-J (2008) Controlling the mass action of alpha-synuclein in Parkinson’s disease. J Neurochem 107:303–316.

    Article  CAS  PubMed  Google Scholar 

  12. Olanow CW, Brundin P (2013) Parkinson’s disease and alpha synuclein: is Parkinson’s disease a prion-like disorder? Mov Disord 28:31–40.

    Article  CAS  PubMed  Google Scholar 

  13. Li J-Y, Englund E, Holton JL et al (2008) Lewy bodies in grafted neurons in subjects with Parkinson’s disease suggest host-to-graft disease propagation. Nat Med 14:501–503.

    Article  CAS  PubMed  Google Scholar 

  14. Basso M, Giraudo S, Corpillo D et al (2004) Proteome analysis of human substantia nigra in Parkinson’s disease. Proteomics 4:3943–3952.

    Article  CAS  PubMed  Google Scholar 

  15. Kitsou E, Pan S, Zhang J et al (2008) Identification of proteins in human substantia nigra. Proteomics Clin Appl 2:776–782.

    Article  CAS  PubMed  Google Scholar 

  16. Licker V, Turck N, Kövari E et al (2014) Proteomic analysis of human substantia nigra identifies novel candidates involved in Parkinson’s disease pathogenesis. Proteomics 14:784–794.

    Article  CAS  PubMed  Google Scholar 

  17. Licker V, Côte M, Lobrinus JA et al (2012) Proteomic profiling of the substantia nigra demonstrates CNDP2 overexpression in Parkinson’s disease. J Proteomics 75:4656–4667.

    Article  CAS  PubMed  Google Scholar 

  18. Dutta D, Ali N, Banerjee E et al (2018) Low levels of prohibitin in substantia nigra makes dopaminergic neurons vulnerable in parkinson’s disease. Mol Neurobiol 55:804–821.

    Article  CAS  PubMed  Google Scholar 

  19. Dumitriu A, Golji J, Labadorf AT et al (2016) Integrative analyses of proteomics and RNA transcriptomics implicate mitochondrial processes, protein folding pathways and GWAS loci in Parkinson disease. BMC Med Genomics 9:5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Ping L, Duong DM, Yin L et al (2018) Global quantitative analysis of the human brain proteome in Alzheimer’s and Parkinson’s Disease. Sci Data 5:180036.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Dixit A, Mehta R, Singh AK (2019) Proteomics in human parkinson’s disease: present scenario and future directions. Cell Mol Neurobiol 39:901–915.

    Article  PubMed  Google Scholar 

  22. Palkovits M (1973) Isolated removal of hypothalamic or other brain nuclei of the rat. Brain Res 59:449–450.

    Article  CAS  PubMed  Google Scholar 

  23. Palkovits M (1986) Microdissection of individual brain nuclei and areas. General Neurochemical Techniques. Humana Press, New Jersey, pp 1–18

    Google Scholar 

  24. Toomey CE, Heywood W, Benson BC et al (2020) Investigation of pathology, expression and proteomic profiles in human TREM2 variant postmortem brains with and without Alzheimer’s disease. Brain Pathol 30:794–810.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Bennett K, Callard R, Heywood W et al (2010) New role for LEKTI in skin barrier formation: label-free quantitative proteomic identification of caspase 14 as a novel target for the protease inhibitor LEKTI. J Proteome Res 9:4289–4294.

    Article  CAS  PubMed  Google Scholar 

  26. Bostanci N, Heywood W, Mills K et al (2010) Application of label-free absolute quantitative proteomics in human gingival crevicular fluid by LC/MS E (gingival exudatome). J Proteome Res 9:2191–2199.

    Article  CAS  PubMed  Google Scholar 

  27. Coats CJ, Heywood WE, Virasami A et al (2018) Proteomic analysis of the myocardium in hypertrophic obstructive cardiomyopathy. Circ Genom Precis Med 11:e001974.

    Article  CAS  PubMed  Google Scholar 

  28. Bliss E, Heywood WE, Benatti M et al (2016) An optimised method for the proteomic profiling of full thickness human skin. Biol Proced Online 18:15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Zhang B, Kirov S, Snoddy J (2005) WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res 33:W741–W748.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Wang J, Vasaikar S, Shi Z et al (2017) WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucl Acid Res 45:W130–W137.

    Article  CAS  Google Scholar 

  31. Huang DW, Sherman BT, Lempicki RA (2009) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4:44–57.

    Article  CAS  Google Scholar 

  32. Smith AC, Robinson AJ (2016) MitoMiner v3.1, an update on the mitochondrial proteomics database. Nucl Acid Res 44:D1258–D1261.

    Article  CAS  Google Scholar 

  33. Thomas PD, Campbell MJ, Kejariwal A et al (2003) PANTHER: a library of protein families and subfamilies indexed by function. Genome Res 13:2129–2141.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Sharma K, Schmitt S, Bergner CG et al (2015) Cell type- and brain region-resolved mouse brain proteome. Nat Neurosci 18:1819–1831.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Wang D, Liu S, Warrell J et al (2018) Comprehensive functional genomic resource and integrative model for the human brain. Science.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Zhu B, Park J-M, Coffey S et al (2022) Single-cell transcriptomic and proteomic analysis of Parkinson’s disease Brains. BioRxiv.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Zhong J, Tang G, Zhu J et al (2021) Single-cell brain atlas of Parkinson’s disease mouse model. J Genet Genomics 48:277–288.

    Article  CAS  PubMed  Google Scholar 

  38. GitHub - shenlab-sinai/GeneOverlap: R package for testing and visualizing gene list overlaps. Accessed 29 Jul 2022

  39. Wingo AP, Dammer EB, Breen MS et al (2019) Large-scale proteomic analysis of human brain identifies proteins associated with cognitive trajectory in advanced age. Nat Commun 10:1619.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Heywood WE, Galimberti D, Bliss E et al (2015) Identification of novel CSF biomarkers for neurodegeneration and their validation by a high-throughput multiplexed targeted proteomic assay. Mol Neurodeg 10:64.

    Article  CAS  Google Scholar 

  41. Hargreaves IP, Heales SJR, Briddon A et al (2002) Blood mononuclear cell coenzyme Q10 concentration and mitochondrial respiratory chain succinate cytochrome-c reductase activity in phenylketonuric patients. J Inherit Metab Dis 25:673–679.

    Article  CAS  PubMed  Google Scholar 

  42. Calvo SE, Clauser KR, Mootha VK (2016) MitoCarta2.0: an updated inventory of mammalian mitochondrial proteins. Nucl Acid Res 44:D1251–D1257.

    Article  CAS  Google Scholar 

  43. Schönfeld P, Reiser G (2013) Why does brain metabolism not favor burning of fatty acids to provide energy? Reflections on disadvantages of the use of free fatty acids as fuel for brain. J Cereb Blood Flow Metab 33:1493–1499.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Schapira AH, Cooper JM, Dexter D et al (1989) Mitochondrial complex I deficiency in Parkinson’s disease. Lancet 1:1269.

    Article  CAS  PubMed  Google Scholar 

  45. Parker WD, Parks JK, Swerdlow RH (2008) Complex I deficiency in Parkinson’s disease frontal cortex. Brain Res 1189:215–218.

    Article  CAS  PubMed  Google Scholar 

  46. Piehowski PD, Petyuk VA, Orton DJ et al (2013) Sources of technical variability in quantitative LC-MS proteomics: human brain tissue sample analysis. J Proteome Res 12:2128–2137.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Ghergurovich JM, García-Cañaveras JC, Wang J et al (2020) A small molecule G6PD inhibitor reveals immune dependence on pentose phosphate pathway. Nat Chem Biol 16:731–739.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Dunn L, Allen GF, Mamais A et al (2014) Dysregulation of glucose metabolism is an early event in sporadic Parkinson’s disease. Neurobiol Aging 35:1111–1115.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Owen JB, Butterfield DA (2010) Measurement of oxidized/reduced glutathione ratio. Methods Mol Biol 648:269–277.

    Article  CAS  PubMed  Google Scholar 

  50. Fico A, Paglialunga F, Cigliano L et al (2004) Glucose-6-phosphate dehydrogenase plays a crucial role in protection from redox-stress-induced apoptosis. Cell Death Differ 11:823–831.

    Article  CAS  PubMed  Google Scholar 

  51. Marí M, Morales A, Colell A et al (2009) Mitochondrial glutathione, a key survival antioxidant. Antioxid Redox Signal 11:2685–2700.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Cumming RC, Andon NL, Haynes PA et al (2004) Protein disulfide bond formation in the cytoplasm during oxidative stress. J Biol Chem 279:21749–21758.

    Article  CAS  PubMed  Google Scholar 

  53. Janero DR, Hreniuk D (1996) Suppression of TCA cycle activity in the cardiac muscle cell by hydroperoxide-induced oxidant stress. Am J Physiol 270:C1735–C1742.

    Article  CAS  PubMed  Google Scholar 

  54. González-Rodríguez P, Zampese E, Stout KA et al (2021) Disruption of mitochondrial complex I induces progressive parkinsonism. Nature 599:650–656.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Lemarie A, Grimm S (2011) Mitochondrial respiratory chain complexes: apoptosis sensors mutated in cancer? Oncogene 30:3985–4003.

    Article  CAS  PubMed  Google Scholar 

  56. Faucheux BA, Martin M-E, Beaumont C et al (2003) Neuromelanin associated redox-active iron is increased in the substantia nigra of patients with Parkinson’s disease. J Neurochem 86:1142–1148.

    Article  CAS  PubMed  Google Scholar 

  57. Jenner P, Olanow CW (1998) Understanding cell death in Parkinson’s disease. Ann Neurol 44:S72-84.

    Article  CAS  PubMed  Google Scholar 

  58. Werner CJ, Heyny-von Haussen R, Mall G, Wolf S (2008) Proteome analysis of human substantia nigra in Parkinson’s disease. Proteome Sci 6:8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Garcia-Esparcia P, Llorens F, Carmona M, Ferrer I (2014) Complex deregulation and expression of cytokines and mediators of the immune response in Parkinson’s disease brain is region dependent. Brain Pathol 24:584–598.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Garcia-Esparcia P, Koneti A, Rodríguez-Oroz MC et al (2018) Mitochondrial activity in the frontal cortex area 8 and angular gyrus in Parkinson’s disease and Parkinson’s disease with dementia. Brain Pathol 28:43–57.

    Article  CAS  PubMed  Google Scholar 

  61. Ferrer I, Martinez A, Blanco R et al (2011) Neuropathology of sporadic Parkinson disease before the appearance of parkinsonism: preclinical Parkinson disease. J Neural Transm 118:821–839.

    Article  PubMed  Google Scholar 

  62. Wong DL, Merrifield-MacRae ME, Stillman MJ (2017) Lead(ii) binding in metallothioneins. Met Ions Life Sci.

    Article  PubMed  Google Scholar 

  63. Baird SK, Kurz T, Brunk UT (2006) Metallothionein protects against oxidative stress-induced lysosomal destabilization. Biochem J 394:275–283.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Ebadi M, Sharma S (2006) Metallothioneins 1 and 2 attenuate peroxynitrite-induced oxidative stress in Parkinson disease. Exp Biol Med 231:1576–1583.

    Article  CAS  Google Scholar 

  65. Sharma SK, Ebadi M (2003) Metallothionein attenuates 3-morpholinosydnonimine (SIN-1)-induced oxidative stress in dopaminergic neurons. Antioxid Redox Signal 5:251–264.

    Article  CAS  PubMed  Google Scholar 

  66. Mohamad Najib NH, Yahaya MF, Das S, Teoh SL (2021) The effects of metallothionein in paraquat-induced Parkinson disease model of zebrafish. Int J Neurosci.

    Article  PubMed  Google Scholar 

  67. Chung RS, Penkowa M, Dittmann J et al (2008) Redefining the role of metallothionein within the injured brain: extracellular metallothioneins play an important role in the astrocyte-neuron response to injury. J Biol Chem 283:15349–15358.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Waller R, Murphy M, Garwood CJ et al (2018) Metallothionein-I/II expression associates with the astrocyte DNA damage response and not Alzheimer-type pathology in the aging brain. Glia 66:2316–2323.

    Article  PubMed  Google Scholar 

  69. Sauve AA, Youn DY (2012) Sirtuins: NAD(+)-dependent deacetylase mechanism and regulation. Curr Opin Chem Biol 16:535–543.

    Article  CAS  PubMed  Google Scholar 

  70. Yang T, Sauve AA (2006) NAD metabolism and sirtuins: metabolic regulation of protein deacetylation in stress and toxicity. AAPS J 8:E632–E643.

    Article  PubMed  PubMed Central  Google Scholar 

  71. North BJ, Verdin E (2004) Sirtuins: Sir2-related NAD-dependent protein deacetylases. Genome Biol 5:224.

    Article  PubMed  PubMed Central  Google Scholar 

  72. Yeong KY, Berdigaliyev N, Chang Y (2020) Sirtuins and their implications in neurodegenerative diseases from a drug discovery perspective. ACS Chem Neurosci 11:4073–4091.

    Article  CAS  PubMed  Google Scholar 

  73. Liu G, Park S-H, Imbesi M et al (2017) Loss of NAD-Dependent Protein Deacetylase Sirtuin-2 Alters Mitochondrial Protein Acetylation And Dysregulates Mitophagy. Antioxid Redox Signal 26:849–863.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Chamberlain KA, Huang N, Xie Y et al (2021) Oligodendrocytes enhance axonal energy metabolism by deacetylation of mitochondrial proteins through transcellular delivery of SIRT2. Neuron 109:3456-3472.e8.

    Article  CAS  PubMed  Google Scholar 

  75. de Oliveira RM, Vicente Miranda H, Francelle L et al (2017) The mechanism of sirtuin 2-mediated exacerbation of alpha-synuclein toxicity in models of Parkinson disease. PLoS Biol 15:e2000374.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Outeiro TF, Kontopoulos E, Altmann SM et al (2007) Sirtuin 2 inhibitors rescue alpha-synuclein-mediated toxicity in models of Parkinson’s disease. Science 317:516–519.

    Article  CAS  PubMed  Google Scholar 

  77. Chen X, Wales P, Quinti L et al (2015) The sirtuin-2 inhibitor AK7 is neuroprotective in models of Parkinson’s disease but not amyotrophic lateral sclerosis and cerebral ischemia. PLoS ONE 10:e0116919.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Wang X, Wang M, Yang L et al (2014) Inhibition of sirtuin 2 exerts neuroprotection in aging rats with increased neonatal iron intake. Neural Regen Res 9:1917–1922.

    Article  PubMed  PubMed Central  Google Scholar 

  79. Zhang J, Wang C, Nie H et al (2016) SIRT2 plays a significant role in maintaining the survival and energy metabolism of PIEC endothelial cells. Int J Physiol Pathophysiol Pharmacol 8:120–127

    PubMed  PubMed Central  Google Scholar 

  80. Singh CK, Chhabra G, Ndiaye MA et al (2018) The role of sirtuins in antioxidant and redox signaling. Antioxid Redox Sig 28:643–661.

    Article  CAS  Google Scholar 

  81. Li W, Zhang B, Tang J et al (2007) Sirtuin 2, a mammalian homolog of yeast silent information regulator-2 longevity regulator, is an oligodendroglial protein that decelerates cell differentiation through deacetylating alpha-tubulin. J Neurosci 27:2606–2616.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Agarwal D, Sandor C, Volpato V et al (2020) A single-cell atlas of the human substantia nigra reveals cell-specific pathways associated with neurological disorders. Nat Commun 11:4183.

    Article  PubMed  PubMed Central  Google Scholar 

  83. Feleke R, Reynolds RH, Smith AM et al (2021) Cross-platform transcriptional profiling identifies common and distinct molecular pathologies in Lewy body diseases. Acta Neuropathol 142:449–474.

    Article  PubMed  PubMed Central  Google Scholar 

  84. Bryois J, Skene NG, Hansen TF et al (2020) Genetic identification of cell types underlying brain complex traits yields insights into the etiology of Parkinson’s disease. Nat Genet 52:482–493.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Olanow CW, Stern MB, Sethi K (2009) The scientific and clinical basis for the treatment of parkinson disease (2009). Neurology 72:S1-136.

    Article  PubMed  Google Scholar 

Download references


The study, CT and SP were funded by Parkinson’s UK. SG was supported by Wellcome (100172/Z/12/2), and is an MRC Senior Clinical Fellow (MR/T008199/1). TR is partly supported by the National Institute for Health Research (NIHR) Queen Square Biomedical Research Unit in Dementia based at University College London Hospitals (UCLH), University College London (UCL). WH and KM conduct research at Great Ormond Street Hospital NHS Foundation Trust and UCL Great Ormond Street Institute of Child Health which is made possible by the NIHR Great Ormond Street Hospital Biomedical Research Centre. The work was undertaken at Queen Square Brain Bank for Neurological disorders. We are grateful for the help and guidance from Prof Janice Holton. Tissue was collected from: the Parkinson’s UK Brain Bank at Imperial College London (Registered charity in England and Wales (258197) and in Scotland (SC037554); Neurodegenerative Diseases Brain Bank, Institute of Psychiatry, Psychology & Neuroscience, King's College London; Newcastle Brain Tissue Resource, Institute of Neuroscience and Newcastle University Institute for Ageing; Neurological Tissue bank, University of Barcelona; Human Brain Tissue Bank, Budapest, Semmelweis University; and Netherlands Brain Bank, Amsterdam. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health.


The research leading to these results received funding from Parkinson’s UK under Grant Agreement G-1004.

Author information

Authors and Affiliations



All authors contributed to the study conception and design. Material preparation were performed by CET, SP, CT, JA, EG, MP, SG and TR. Data collection and analysis were performed by CET, WH, JRE, JL, SP, MAS, IH, MO and KD. The first draft of the manuscript was written by CET, WH, SP, KM and SG and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Christina E. Toomey or Sonia Gandhi.

Ethics declarations

Ethical approval and consent to participate

Informed consent was given for all cases. Ethical approval was obtained for the use of post-mortem human brain tissue from the Local Research Ethics Committee of the National Hospital for Neurology and Neurosurgery.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

GO terms altered in early PD compared to controls and mitochondrial gene expression patterns (a) Bar charts indicating number of Gene Ontology (GO) terms that are represented in the dataset for each region as determined by Webgestalt and DAVID databases. Top chart shows Biological processes, middle the molecular functions and lower the cellular components. The pie chart adjacent to each chart shows the number of GO terms that overlap or are uniquely represented in the dataset for a region of severe pathology (substantia nigra), mild pathology (parahippocampus) and a region unaffected at Braak stage 3/4 (frontal cortex) as determined by GOview. Colours indicate the level each region is affected at Braak stage 3/4 as determined in Figure 1b. Protein heatmaps from (b) first mass spectrometry run and (c) second mass spectrometry run from IPA showing the level of expression change per protein in the Mitochondrial Dysfunction pathway across each brain region in Braak stage 3/4 compared to controls. Red indicated upregulation and green indicated downregulation compared to controls. Intensity of colour shows level of expression change with deeper colour indicating higher up- or down- regulation.

Additional file 2: Figure S2.

Candidates for validation workflow A pie chart showing the proportion of mitochondrial proteins out of the total proteins that were detected. The box outlines the method used to select candidate proteins for validation. The final pie chart shows the proportion of mitochondrial proteins within the list of proteins that were selected for validation.

Additional file 3: Figure S3.

Citrate synthase protein levels across the brain. Graph highlighting how much citrate synthase protein was detected per disease group and brain region. Two-way ANOVA with Sidak’s multiple comparisons determined that citrate synthase levels across disease groups were non-significant (p = 0.1647) whilst those between regions were significant (p < 0.0001). Figure and analysis made with GraphPad Prism 9.02.

Additional file 4:Table S1.

Peptide list for candidates. Peptide chosen for each protein that was validated with multiple reaction monitoring.

Additional file 5:Table S2.

Biological processes represented across brain regions for early PD cases compared to controls. Shaded box reflects  GO term represented in that brain region.

Additional file 6:Table S3.

Molecular functions represented across brain regions for early PD cases compared to controls. Shaded box reflects GO term represented in that brain region.

Additional file 7:Table S4.

Cellular components represented across brain regions for early PD cases compared to controls Shaded box reflects GO term represented in that brain region.

Additional file 8: Table S5.

Comparison of candidate proteins to cell type transcriptomic and proteomic data. Each protein that was compared to downloaded datasets from the articles listed to determine which cell type they associated with. * indicates < 5% FDR significant association in enrichment analysis.

Additional File 9: Table S6.

Proteins with > 1.5 fold expression change within each region for early PD compared to controls. Protein ID is shown for each protein in the list. Red shading indicates upregulation and green shading indicates downregulation. Each region is highlighted to show the severity of pathology present in that region at Braak stage 3/4 with deepest severely affected > moderately affected > mildly affected > not affected.

Additional file 10: Table S7.

List of GO terms represented between early and late PD cases in the frontal cortex. Biological processes, molecular functions and cellular components are split into different columns. Top 30 GO terms that overlap between both sets of brains compared to controls are given as well as those only found in early PD and only found in late PD.

Additional file 11: Table S8.

List of GO terms represented between early and late PD cases in putamen. Biological processes, molecular functions and cellular components are split into different columns. Top 30 GO terms that overlap between both sets of brains compared to controls are given as well as those only found in early PD and only found in late PD.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Toomey, C.E., Heywood, W.E., Evans, J.R. et al. Mitochondrial dysfunction is a key pathological driver of early stage Parkinson’s. acta neuropathol commun 10, 134 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: