Mitochondrial complex I deficiency stratifies idiopathic Parkinson's disease.

Idiopathic Parkinson's disease (iPD) is believed to have a heterogeneous pathophysiology, but molecular disease subtypes have not been identified. Here, we show that iPD can be stratified according to the severity of neuronal respiratory complex I (CI) deficiency, and identify two emerging disease subtypes with distinct molecular and clinical profiles. The CI deficient (CI-PD) subtype accounts for approximately a fourth of all cases, and is characterized by anatomically widespread neuronal CI deficiency, a distinct cell type-specific gene expression profile, increased load of neuronal mtDNA deletions, and a predilection for non-tremor dominant motor phenotypes. In contrast, the non-CI deficient (nCI-PD) subtype exhibits no evidence of mitochondrial impairment outside the dopaminergic substantia nigra and has a predilection for a tremor dominant phenotype. These findings constitute a step towards resolving the biological heterogeneity of iPD with implications for both mechanistic understanding and treatment strategies.

DNA Methylation Age Acceleration Is Not Associated with Age of Onset in Parkinson's Disease.

BACKGROUND: Epigenetic clocks using DNA methylation (DNAm) to estimate biological age have become popular tools in the study of neurodegenerative diseases. Notably, several recent reports have shown a strikingly similar inverse relationship between accelerated biological aging, as measured by DNAm, and the age of onset of several neurodegenerative disorders, including Parkinson's disease (PD). Common to all of these studies is that they were performed without control subjects and using the exact same measure of accelerated aging: DNAm age minus chronological age. OBJECTIVE: We aimed to assess the validity of these findings in PD, using the same dataset as in the original study, blood DNAm data from the Parkinson's Progression Markers Initiative cohort, but also including control samples in the analyses. METHODS: We replicated the analyses and findings of the previous study and then reanalyzed the dataset incorporating control samples to account for underlying age-related biases. RESULTS: Our reanalysis shows that there is no correlation between age of onset and DNAm age acceleration. Conversely, there is a pattern of overestimating DNAm age in younger and underestimating DNAm age in older individuals in the dataset that entirely explains the previously reported association. CONCLUSIONS: Our findings refute the previously reported inverse relationship between DNAm age acceleration and age of onset in PD. We show that these findings are fully accounted for by an expected over/underestimation of DNAm age in younger/older individuals. Furthermore, this effect is likely to be responsible for nearly identical findings reported in other neurodegenerative diseases. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Not every estimate counts - evaluation of cell composition estimation approaches in brain bulk tissue data.

BACKGROUND: Variation in cell composition can dramatically impact analyses in bulk tissue samples. A commonly employed approach to mitigate this issue is to adjust statistical models using estimates of cell abundance derived directly from omics data. While an arsenal of estimation methods exists, the applicability of these methods to brain tissue data and whether or not cell estimates can sufficiently account for confounding cellular composition has not been adequately assessed. METHODS: We assessed the correspondence between different estimation methods based on transcriptomic (RNA sequencing, RNA-seq) and epigenomic (DNA methylation and histone acetylation) data from brain tissue samples of 49 individuals. We further evaluated the impact of different estimation approaches on the analysis of H3K27 acetylation chromatin immunoprecipitation sequencing (ChIP-seq) data from entorhinal cortex of individuals with Alzheimer's disease and controls. RESULTS: We show that even closely adjacent tissue samples from the same Brodmann area vary greatly in their cell composition. Comparison across different estimation methods indicates that while different estimation methods applied to the same data produce highly similar outcomes, there is a surprisingly low concordance between estimates based on different omics data modalities. Alarmingly, we show that cell type estimates may not always sufficiently account for confounding variation in cell composition. CONCLUSIONS: Our work indicates that cell composition estimation or direct quantification in one tissue sample should not be used as a proxy to the cellular composition of another tissue sample from the same brain region of an individual-even if the samples are directly adjacent. The highly similar outcomes observed among vastly different estimation methods, highlight the need for brain benchmark datasets and better validation approaches. Finally, unless validated through complementary experiments, the interpretation of analyses outcomes based on data confounded by cell composition should be done with great caution, and ideally avoided all together.

Nicotinamide riboside supplementation is not associated with altered methylation homeostasis in Parkinson's disease.

Replenishing nicotinamide adenine dinucleotide (NAD) via supplementation of nicotinamide riboside (NR) has been shown to confer neuroprotective effects in models of aging and neurodegenerative diseases, including Parkinson's disease (PD). Although generally considered safe, concerns have been raised that NR supplementation could impact methylation dependent reactions, including DNA methylation, because of increased production and methylation dependent breakdown of nicotinamide (NAM). We investigated the effect of NR supplementation on DNA methylation in a double blinded, placebo-controlled trial of 29 human subjects with PD, in blood cells and muscle tissue. Our results show that NR had no impact on DNA methylation homeostasis, including individuals with common pathogenic mutations in the MTHFR gene known to affect one-carbon metabolism. Pathway and methylation variance analyses indicate that there might be minor regulatory responses to NR. We conclude that short-term therapy with high-dose NR for up to 30 days has no deleterious impact on methylation homeostasis.

Altered transcriptome-proteome coupling indicates aberrant proteostasis in Parkinson's disease.

Aberrant proteostasis is thought to be implicated in Parkinson's disease (PD), but patient-derived evidence is scant. We hypothesized that impaired proteostasis is reflected as altered transcriptome-proteome correlation in the PD brain. We integrated transcriptomic and proteomic data from prefrontal cortex of PD patients and young and aged controls to assess RNA-protein correlations across samples. The aged brain showed a genome-wide decrease in mRNA-protein correlation. Genes encoding synaptic vesicle proteins showed negative correlations, likely reflecting spatial separation of mRNA and protein into soma and synapses. PD showed a broader transcriptome-proteome decoupling, consistent with a proteome-wide decline in proteostasis. Genes showing negative correlation in PD were enriched for proteasome subunits, indicating accentuated spatial separation of transcript and protein in PD neurons. In addition, PD showed positive correlations for mitochondrial respiratory chain genes, suggesting a tighter regulation in the face of mitochondrial dysfunction. Our results support the hypothesis that aberrant proteasomal function is implicated in PD pathogenesis.

No evidence of extensive non-CpG methylation in mtDNA.

While most research suggests mitochondrial DNA (mtDNA) harbors low or no methylation, a few studies claim to report evidence of high-level methylation in the mtDNA. The reasons behind these contradictory results are likely to be methodological but remain largely unexplored. Here, we critically reanalyzed a recent study by Patil et al. (2019) reporting extensive methylation in human mtDNA in a non-CpG context. Our analyses refute the original findings and show that these do not reflect the biology of the tested samples, but rather stem from a combination of methodological and technical pitfalls. The authors employ an oversimplified model that defines as methylated all reference positions with methylation proportions above an arbitrary cutoff of 9%. This substantially exacerbates the overestimation of methylated cytosines due to the selective degradation of unmethylated cytosine-rich regions. Additional limitations are the small sample sizes and lack of sample-specific controls for bisulfite conversion efficiency. In conclusion, using the same dataset employed in the original study by Patil et al., we find no evidence supporting the existence of extensive non-CpG methylation in the human mtDNA.

Ultra-deep whole genome bisulfite sequencing reveals a single methylation hotspot in human brain mitochondrial DNA.

While DNA methylation is established as a major regulator of gene expression in the nucleus, the existence of mitochondrial DNA (mtDNA) methylation remains controversial. Here, we characterized the mtDNA methylation landscape in the prefrontal cortex of neurological healthy individuals (n=26) and patients with Parkinson's disease (n=27), using a combination of whole-genome bisulphite sequencing (WGBS) and bisulphite-independent methods. Accurate mtDNA mapping from WGBS data required alignment to an mtDNA reference only, to avoid misalignment to nuclear mitochondrial pseudogenes. Once correctly aligned, WGBS data provided ultra-deep mtDNA coverage (16,723 ± 7,711) and revealed overall very low levels of cytosine methylation. The highest methylation levels (5.49 ± 0.97%) were found on CpG position m.545, located in the heavy-strand promoter 1 region. The m.545 methylation was validated using a combination of methylation-sensitive DNA digestion and quantitative PCR analysis. We detected no association between mtDNA methylation profile and Parkinson's disease. Interestingly, m.545 methylation correlated with the levels of mtDNA transcripts, suggesting a putative role in regulating mtDNA gene expression. In addition, we propose a robust framework for methylation analysis of mtDNA from WGBS data, which is less prone to false-positive findings due to misalignment of nuclear mitochondrial pseudogene sequences.

The NADPARK study: A randomized phase I trial of nicotinamide riboside supplementation in Parkinson's disease.

We conducted a double-blinded phase I clinical trial to establish whether nicotinamide adenine dinucleotide (NAD) replenishment therapy, via oral intake of nicotinamide riboside (NR), is safe, augments cerebral NAD levels, and impacts cerebral metabolism in Parkinson's disease (PD). Thirty newly diagnosed, treatment-naive patients received 1,000 mg NR or placebo for 30 days. NR treatment was well tolerated and led to a significant, but variable, increase in cerebral NAD levels-measured by 31phosphorous magnetic resonance spectroscopy-and related metabolites in the cerebrospinal fluid. NR recipients showing increased brain NAD levels exhibited altered cerebral metabolism, measured by 18fluoro-deoxyglucose positron emission tomography, and this was associated with mild clinical improvement. NR augmented the NAD metabolome and induced transcriptional upregulation of processes related to mitochondrial, lysosomal, and proteasomal function in blood cells and/or skeletal muscle. Furthermore, NR decreased the levels of inflammatory cytokines in serum and cerebrospinal fluid. Our findings nominate NR as a potential neuroprotective therapy for PD, warranting further investigation in larger trials.

Distinct Mitochondrial Remodeling During Mesoderm Differentiation in a Human-Based Stem Cell Model.

Given the considerable interest in using stem cells for modeling and treating disease, it is essential to understand what regulates self-renewal and differentiation. Remodeling of mitochondria and metabolism, with the shift from glycolysis to oxidative phosphorylation (OXPHOS), plays a fundamental role in maintaining pluripotency and stem cell fate. It has been suggested that the metabolic "switch" from glycolysis to OXPHOS is germ layer-specific as glycolysis remains active during early ectoderm commitment but is downregulated during the transition to mesoderm and endoderm lineages. How mitochondria adapt during these metabolic changes and whether mitochondria remodeling is tissue specific remain unclear. Here, we address the question of mitochondrial adaptation by examining the differentiation of human pluripotent stem cells to cardiac progenitors and further to differentiated mesodermal derivatives, including functional cardiomyocytes. In contrast to recent findings in neuronal differentiation, we found that mitochondrial content decreases continuously during mesoderm differentiation, despite increased mitochondrial activity and higher levels of ATP-linked respiration. Thus, our work highlights similarities in mitochondrial remodeling during the transition from pluripotent to multipotent state in ectodermal and mesodermal lineages, while at the same time demonstrating cell-lineage-specific adaptations upon further differentiation. Our results improve the understanding of how mitochondrial remodeling and the metabolism interact during mesoderm differentiation and show that it is erroneous to assume that increased OXPHOS activity during differentiation requires a simultaneous expansion of mitochondrial content.

Genome-wide histone acetylation analysis reveals altered transcriptional regulation in the Parkinson's disease brain.

BACKGROUND: Parkinson's disease (PD) is a complex, age-related neurodegenerative disorder of largely unknown etiology. PD is strongly associated with mitochondrial respiratory dysfunction, which can lead to epigenetic dysregulation and specifically altered histone acetylation. Nevertheless, and despite the emerging role of epigenetics in age-related brain disorders, the question of whether aberrant histone acetylation is involved in PD remains unresolved. METHODS: We studied fresh-frozen brain tissue from two independent cohorts of individuals with idiopathic PD (n = 28) and neurologically healthy controls (n = 21). We performed comprehensive immunoblotting to identify histone sites with altered acetylation levels in PD, followed by chromatin immunoprecipitation sequencing (ChIP-seq). RNA sequencing data from the same individuals was used to assess the impact of altered histone acetylation on gene expression. RESULTS: Immunoblotting analyses revealed increased acetylation at several histone sites in PD, with the most prominent change observed for H3K27, a marker of active promoters and enhancers. ChIP-seq analysis further indicated that H3K27 hyperacetylation in the PD brain is a genome-wide phenomenon with a strong predilection for genes implicated in the disease, including SNCA, PARK7, PRKN and MAPT. Integration of the ChIP-seq with transcriptomic data from the same individuals revealed that the correlation between promoter H3K27 acetylation and gene expression is attenuated in PD patients, suggesting that H3K27 acetylation may be decoupled from transcription in the PD brain. Strikingly, this decoupling was most pronounced among nuclear-encoded mitochondrial genes, corroborating the notion that impaired crosstalk between the nucleus and mitochondria is involved in the pathogenesis of PD. Our findings independently replicated in the two cohorts. CONCLUSIONS: Our findings strongly suggest that aberrant histone acetylation and altered transcriptional regulation are involved in the pathophysiology of PD. We demonstrate that PD-associated genes are particularly prone to epigenetic dysregulation and identify novel epigenetic signatures associated with the disease.

Differential transcript usage in the Parkinson's disease brain.

Studies of differential gene expression have identified several molecular signatures and pathways associated with Parkinson's disease (PD). The role of isoform switches and differential transcript usage (DTU) remains, however, unexplored. Here, we report the first genome-wide study of DTU in PD. We performed RNA sequencing following ribosomal RNA depletion in prefrontal cortex samples of 49 individuals from two independent case-control cohorts. DTU was assessed using two transcript-count based approaches, implemented in the DRIMSeq and DEXSeq tools. Multiple PD-associated DTU events were detected in each cohort, of which 23 DTU events in 19 genes replicated across both patient cohorts. For several of these, including THEM5, SLC16A1 and BCHE, DTU was predicted to have substantial functional consequences, such as altered subcellular localization or switching to non-protein coding isoforms. Furthermore, genes with PD-associated DTU were enriched in functional pathways previously linked to PD, including reactive oxygen species generation and protein homeostasis. Importantly, the vast majority of genes exhibiting DTU were not differentially expressed at the gene-level and were therefore not identified by conventional differential gene expression analysis. Our findings provide the first insight into the DTU landscape of PD and identify novel disease-associated genes. Moreover, we show that DTU may have important functional consequences in the PD brain, since it is predicted to alter the functional composition of the proteome. Based on these results, we propose that DTU analysis is an essential complement to differential gene expression studies in order to provide a more accurate and complete picture of disease-associated transcriptomic alterations.

Using urine to diagnose large-scale mtDNA deletions in adult patients.

OBJECTIVE: The aim of this study was to evaluate if urinary sediment cells offered a robust alternative to muscle biopsy for the diagnosis of single mtDNA deletions. METHODS: Eleven adult patients with progressive external ophthalmoplegia and a known single mtDNA deletion were investigated. Urinary sediment cells were used to isolate DNA, which was then subjected to long-range polymerase chain reaction. Where available, the patient`s muscle DNA was studied in parallel. Breakpoint and thus deletion size were identified using both Sanger sequencing and next generation sequencing. The level of heteroplasmy was determined using quantitative polymerase chain reaction. RESULTS: We identified the deletion in urine in 9 of 11 cases giving a sensitivity of 80%. Breakpoints and deletion size were readily detectable in DNA extracted from urine. Mean heteroplasmy level in urine was 38% ± 26 (range 8 - 84%), and 57% ± 28 (range 12 - 94%) in muscle. While the heteroplasmy level in urinary sediment cells differed from that in muscle, we did find a statistically significant correlation between these two levels (R = 0.714, P = 0.031(Pearson correlation)). INTERPRETATION: Our findings suggest that urine can be used to screen patients suspected clinically of having a single mtDNA deletion. Based on our data, the use of urine could considerably reduce the need for muscle biopsy in this patient group.

Common gene expression signatures in Parkinson's disease are driven by changes in cell composition.

The etiology of Parkinson's disease is largely unknown. Genome-wide transcriptomic studies in bulk brain tissue have identified several molecular signatures associated with the disease. While these studies have the potential to shed light into the pathogenesis of Parkinson's disease, they are also limited by two major confounders: RNA post-mortem degradation and heterogeneous cell type composition of bulk tissue samples. We performed RNA sequencing following ribosomal RNA depletion in the prefrontal cortex of 49 individuals from two independent case-control cohorts. Using cell type specific markers, we estimated the cell type composition for each sample and included this in our analysis models to compensate for the variation in cell type proportions. Ribosomal RNA depletion followed by capture by random primers resulted in substantially more even transcript coverage, compared to poly(A) capture, in post-mortem tissue. Moreover, we show that cell type composition is a major confounder of differential gene expression analysis in the Parkinson's disease brain. Accounting for cell type proportions attenuated numerous transcriptomic signatures that have been previously associated with Parkinson's disease, including vesicle trafficking, synaptic transmission, immune and mitochondrial function. Conversely, pathways related to endoplasmic reticulum, lipid oxidation and unfolded protein response were strengthened and surface as the top differential gene expression signatures in the Parkinson's disease prefrontal cortex. Our results indicate that differential gene expression signatures in Parkinson's disease bulk brain tissue are significantly confounded by underlying differences in cell type composition. Modeling cell type heterogeneity is crucial in order to unveil transcriptomic signatures that represent regulatory changes in the Parkinson's disease brain and are, therefore, more likely to be associated with underlying disease mechanisms.

Rare genetic variation in mitochondrial pathways influences the risk for Parkinson's disease.

BACKGROUND: Mitochondrial dysfunction plays a key role in PD, but the underlying molecular mechanisms remain unresolved. We hypothesized that the disruption of mitochondrial function in PD is primed by rare, protein-altering variation in nuclear genes controlling mitochondrial structure and function. OBJECTIVE: The objective of this study was to assess whether genetic variation in genes associated with mitochondrial function influences the risk of idiopathic PD. METHODS: We employed whole-exome sequencing data from 2 independent cohorts of clinically validated idiopathic PD and controls, the Norwegian ParkWest cohort (n = 411) and the North American Parkinson's Progression Markers Initiative (n = 640). We applied burden-based and variance-based collapsing methods to assess the enrichment of rare, nonsynonymous, and damaging genetic variants on genes, exome-wide, and on a comprehensive set of mitochondrial pathways, defined as groups of genes controlling specific mitochondrial functions. RESULTS: Using the sequence kernel association test, we detected a significant polygenic enrichment of rare, nonsynonymous variants in the gene-set encoding the pathway of mitochondrial DNA maintenance. Notably, this was the strongest association in both cohorts and survived multiple testing correction (ParkWest P = 6.3 × 10-3 , Parkinson's Progression Markers Initiative P = 6.9 × 10-5 , metaanalysis P = 3.2 × 10-6 ). CONCLUSIONS: Our results show that the enrichment of rare inherited variation in the pathway controlling mitochondrial DNA replication and repair influences the risk of PD. We propose that this polygenic enrichment contributes to the impairment of mitochondrial DNA homeostasis, thought to be a key mechanism in the pathogenesis of PD, and explains part of the disorder's "missing heritability." © 2018 The Authors. Movement Disorders published by Wiley Periodicals, Inc. on behalf of International Parkinson and Movement Disorder Society.

Ultradeep mapping of neuronal mitochondrial deletions in Parkinson's disease.

Mitochondrial DNA (mtDNA) deletions accumulate with age in postmitotic cells and are associated with aging and neurodegenerative disorders such as Parkinson's disease. Although the exact mechanisms by which deletions form remain elusive, the dominant theory is that they arise spontaneously at microhomologous sites and undergo clonal expansion. We characterize mtDNA deletions at unprecedented resolution in individual substantia nigra neurons from individuals with Parkinson's disease, using ultradeep sequencing. We show that the number of deleted mtDNA species per neuron is substantially higher than previously reported. Moreover, each deleted mtDNA species shows significant differences in sequence composition compared with the remaining mtDNA population, which is highly consistent with independent segregation and clonal expansion. Deletion breakpoints occur consistently in regions of sequence homology, which may be direct or interrupted stretches of tandem repeats. While our results support a crucial role for misannealing in deletion generation, we find no overrepresentation of the 3'-repeat sequence, an observation that is difficult to reconcile with the current view of replication errors as the source of mtDNA deletions.

Defective mitochondrial DNA homeostasis in the substantia nigra in Parkinson disease.

Increased somatic mitochondrial DNA (mtDNA) mutagenesis causes premature aging in mice, and mtDNA damage accumulates in the human brain with aging and neurodegenerative disorders such as Parkinson disease (PD). Here, we study the complete spectrum of mtDNA changes, including deletions, copy-number variation and point mutations, in single neurons from the dopaminergic substantia nigra and other brain areas of individuals with Parkinson disease and neurologically healthy controls. We show that in dopaminergic substantia nigra neurons of healthy individuals, mtDNA copy number increases with age, maintaining the pool of wild-type mtDNA population in spite of accumulating deletions. This upregulation fails to occur in individuals with Parkinson disease, however, resulting in depletion of the wild-type mtDNA population. By contrast, neuronal mtDNA point mutational load is not increased in Parkinson disease. Our findings suggest that dysregulation of mtDNA homeostasis is a key process in the pathogenesis of neuronal loss in Parkinson disease.

Dynamical properties of gene regulatory networks involved in long-term potentiation.

The long-lasting enhancement of synaptic effectiveness known as long-term potentiation (LTP) is considered to be the cellular basis of long-term memory. LTP elicits changes at the cellular and molecular level, including temporally specific alterations in gene networks. LTP can be seen as a biological process in which a transient signal sets a new homeostatic state that is "remembered" by cellular regulatory systems. Previously, we have shown that early growth response (Egr) transcription factors are of fundamental importance to gene networks recruited early after LTP induction. From a systems perspective, we hypothesized that these networks will show less stable architecture, while networks recruited later will exhibit increased stability, being more directly related to LTP consolidation. Using random Boolean network (RBN) simulations we found that the network derived at 24 h was markedly more stable than those derived at 20 min or 5 h post-LTP. This temporal effect on the vulnerability of the networks is mirrored by what is known about the vulnerability of LTP and memory itself. Differential gene co-expression analysis further highlighted the importance of the Egr family and found a rapid enrichment in connectivity at 20 min, followed by a systematic decrease, providing a potential explanation for the down-regulation of gene expression at 24 h documented in our preceding studies. We also found that the architecture exhibited by a control and the 24 h LTP co-expression networks fit well to a scale-free distribution, known to be robust against perturbations. By contrast the 20 min and 5 h networks showed more truncated distributions. These results suggest that a new homeostatic state is achieved 24 h post-LTP. Together, these data present an integrated view of the genomic response following LTP induction by which the stability of the networks regulated at different times parallel the properties observed at the synapse.

Quantifying the evolutionary divergence of protein structures: the role of function change and function conservation.

The molecular clock hypothesis, stating that protein sequences diverge in evolution by accumulating amino acid substitutions at an almost constant rate, played a major role in the development of molecular evolution and boosted quantitative theories of evolutionary change. These studies were extended to protein structures by the seminal paper by Chothia and Lesk, which established the approximate proportionality between structure and sequence divergence. Here we analyse how function influences the relationship between sequence and structure divergence, studying four large superfamilies of evolutionarily related proteins: globins, aldolases, P-loop and NADP-binding. We introduce the contact divergence, which is more consistent with sequence divergence than previously used structure divergence measures. Our main findings are: (1) Small structure and sequence divergences are proportional, consistent with the molecular clock. Approximate validity of the clock is also supported by the analysis of the clustering coefficient of structure similarity networks. (2) Functional constraints strongly limit the structure divergence of proteins performing the same function and may allow to identify incomplete or wrong functional annotations. (3) The rate of structure versus sequence divergence is larger for proteins performing different functions than for proteins performing the same function. We conjecture that this acceleration is due to positive selection for new functions. Accelerations in structure divergence are also suggested by the analysis of the clustering coefficient. (4) For low sequence identity, structural diversity explodes. We conjecture that this explosion is related to functional diversification. (5) Large indels are almost always associated with function changes.