Idiopathic Parkinson’s disease is characterized by a progressive loss of dopaminergic neurons, but the exact disease etiology remains largely unknown. To date, Parkinson’s disease research has mainly focused on nigral dopaminergic neurons, although recent studies suggest disease-related changes also in non-neuronal cells and in midbrain regions beyond the substantia nigra. While there is some evidence for glial involvement in Parkinson’s disease, the molecular mechanisms remain poorly understood. The aim of this study was to characterize the contribution of all cell types of the midbrain to Parkinson’s disease pathology by single-nuclei RNA sequencing and to assess the cell type-specific risk for Parkinson’s disease employing the latest genome-wide association study. We profiled >41 000 single-nuclei transcriptomes of postmortem midbrain from six idiopathic Parkinson’s disease patients and five age-/sex-matched controls. To validate our findings in a spatial context, we utilized immunolabeling of the same tissues. Moreover, we analyzed Parkinson’s disease-associated risk enrichment in genes with cell type-specific expression patterns. We discovered a neuronal cell cluster characterized by CADPS2 overexpression and low TH levels, which was exclusively present in IPD midbrains. Validation analyses in laser-microdissected neurons suggest that this cluster represents dysfunctional dopaminergic neurons. With regard to glial cells, we observed an increase in nigral microglia in Parkinson’s disease patients. Moreover, nigral idiopathic Parkinson’s disease microglia were more amoeboid, indicating an activated state. We also discovered a reduction in idiopathic Parkinson’s disease oligodendrocyte numbers with the remaining cells being characterized by a stress-induced upregulation of S100B. Parkinson’s disease risk variants were associated with glia- and neuron-specific gene expression patterns in idiopathic Parkinson’s disease cases. Furthermore, astrocytes and microglia presented idiopathic Parkinson’s disease-specific cell proliferation and dysregulation of genes related to unfolded protein response and cytokine signaling. While reactive patient astrocytes showed CD44 overexpression, idiopathic Parkinson’s disease-microglia revealed a pro-inflammatory trajectory characterized by elevated levels of IL1B, GPNMB, and HSP90AA1. Taken together, we generated the first single-nuclei RNA sequencing dataset from the idiopathic Parkinson’s disease midbrain, which highlights a disease-specific neuronal cell cluster as well as ‘pan-glial’ activation as a central mechanism in the pathology of the movement disorder. This finding warrants further research into inflammatory signaling and immunomodulatory treatments in Parkinson’s disease.

Parkinson's disease (PD) etiology is associated with genetic and environmental factors that lead to a loss of dopaminergic neurons. However, the functional interpretation of PD-associated risk variants and how other midbrain cells contribute to this neurodegenerative process are poorly understood. Here, we profiled >41,000 single-nuclei transcriptomes of postmortem midbrain tissue from 6 idiopathic PD (IPD) patients and 5 matched controls. We show that PD-risk variants are associated with glia- and neuron-specific gene expression patterns. Furthermore, Microglia and astrocytes presented IPD-specific cell proliferation and dysregulation of genes related to unfolded protein response and cytokine signalling. IPD-microglia revealed a specific pro-inflammatory trajectory. Finally, we discovered a neuronal cell cluster exclusively present in IPD midbrains characterized by CADPS2 overexpression and a high proportion of cycling cells. We conclude that elevated CADPS2 expression is specific to dysfunctional dopaminergic neurons, which have lost their dopaminergic identity and unsuccessful attempt to re-enter the cell cycle.

Mouse models are a critical tool for studying human diseases, particularly developmental disorders, as well as for advancing our general understanding of mammalian biology. However, it has long been suspected that conventional approaches for phenotyping are insufficiently sensitive to detect subtle defects throughout the developing mouse. Here we set out to establish single cell RNA sequencing (sc-RNA-seq) of the whole embryo as a scalable platform for the systematic molecular and cellular phenotyping of mouse genetic models. We applied combinatorial indexing-based sc-RNA-seq to profile 101 embryos of 26 genotypes at embryonic stage E13.5, altogether profiling gene expression in over 1.6M nuclei. The 26 genotypes include 22 mouse mutants representing a range of anticipated severities, from established multisystem disorders to deletions of individual enhancers, as well as the 4 wildtype backgrounds on which these mutants reside. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibited changes in dozens of trajectories (e.g., the pleiotropic consequences of altering the Sox9 regulatory landscape) whereas others showed phenotypes affecting specific subsets of cells. We also identify differences between widely used wildtype strains, compare phenotyping of gain vs. loss of function mutants, and characterise deletions of topological associating domain (TAD) boundaries. Intriguingly, even among these 22 mutants, some changes are shared by heretofore unrelated models, suggesting that developmental pleiotropy might be "decomposable" through further scaling of this approach. Overall, our findings show how single cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.