Molecular Disease Stages of Oligodendrocytic and Neuronal Tau Burden in Progressive Supranuclear Palsy

Abstract

Background Progressive supranuclear palsy (PSP) is a primary tauopathy defined by the accumulation of 4R tau isoforms in neurons, oligodendrocytes and astrocytes. Despite evidence of genetic susceptibility operating through glial cell types, it remains poorly understood how cell type-specific epigenetic-transcriptional programs evolve with progression of tau pathology. Methods We conducted single-nucleus chromatin accessibility (snATACseq) and RNA sequencing (snRNAseq) on postmortem frontal cortex samples from PSP patients (n = 8) and matched controls (n = 8), yielding over 144,000 nuclei passing quality control. Tau pathology burden, including neurofibrillary tangles, coiled bodies, and tufted astrocytes, was quantified on AT8-immunostained sections from the same individuals. We integrated differential gene expression analysis, transcription factor motif enrichment, weighted gene co-expression network analysis, and pseudotime modeling anchored to cell type-specific tau pathology burden to delineate molecular pseudo-progression trajectories. Results In eight cell types, 20 subclasses, and 70 subclusters, PSP brains displayed a selective depletion of certain excitatory deep-layer neurons and oligodendrocyte subclusters, with relative preservation of inhibitory neurons and vascular cells. Genetic risk enrichment was localized to astrocytes and oligodendrocytes, whereas excitatory neurons exhibited the greatest transcriptional dysregulation. Oligodendrocyte pseudo-progression indicated a transition from homeostatic myelination programs (MBP, MOBP) through glucocorticoid-responsive stress (FKBP5, ZBTB16), to compensatory myelination (PLP1, CNP) and proteostasis stress (UCHL1, CYRAB, CLU). Neuronal pseudo- progression revealed early dysregulation of synaptic (RORB2, NRG3, NPTX1), microtubule dynamics (KIF2C, RAB27B, TUBA/B), and survival (MEG3, FTX) pathways, alongside a transient increase in neuron-glia interactions (GRIP, CNTNAP4, ERBB4), converging late on ribosomal translation and vesicular trafficking modules across all neuronal subtypes. Cross-modal integration with independent cerebrospinal fluid proteomics identified a concordant subset of glial reactivity, axonal injury, and synaptic markers jointly dysregulated in inhibitory neurons, oligodendrocytes, and excitatory deep-layer neurons. Conclusion PSP pathogenesis reflects a combination of glial genetic susceptibility and staged, cell type-specific transcriptional dysfunction. Oligodendrocytes transition from myelination-competent states to FKBP5-mediated stress states, while neurons show variably timed loss of synaptic excitability and survival programs, preceded by neuron-glia interactions and followed by convergent ribosomal-proteostatic failure. These cytopathology-anchored trajectories outline a potential pathophysiological sequence and may inform candidate selection for stage-specific therapeutic interventions in PSP.