- 1Department of Neurology, University Hospital Ulm, Ulm, Germany
- 2Department of Neurology, LMU University Hospital, Ludwig-Maximilians-Universität (LMU), Munich, Germany
- 3German Center for Neurodegenerative Diseases (DZNE), Munich, Germany
- 4Munich Cluster for Systems Neurology (SyNergy), Munich, Germany
- 5German Center for Neurodegenerative Diseases (DZNE), Ulm, Germany
- 6Clinical Neuroanatomy, Department of Neurology, University of Ulm, Ulm, Germany
Background and objective: A neuropathological cerebral staging concept for progressive supranuclear palsy (PSP) has been proposed that tau inclusions in PSP may progress in a sequential regional pattern. The objective was to develop a hypothesis-guided region/tract of interest-based (ROI/TOI) approach to use diffusion tensor imaging (DTI) targeted to analyze in vivo the regions that are prone to be involved at each neuropathological stage of PSP.
Methods: Two data cohorts were analyzed: cohort A of 78 PSP patients [55 Richardson’s syndrome (PSP-RS) and 23 PSP with predominant parkinsonism (PSP-P)] and 63 controls, recorded at 3.0T at multiple sites, and a single-site cohort B constituted by 1.5T data of 66 PSP patients (46 PSP-RS and 20 PSP-P) and 44 controls. In cohort A, 21 PSP patients (13 PSP-RS and 8 PSP-P) and 17 controls obtained a follow-up scan after 17 months. Whole brain-based spatial statistics (WBSS) was used to identify the alterations in PSP patients vs. controls. The combined ROI- and TOI-based approach targeted structures that are prone to be involved during the course of PSP.
Results: WBSS demonstrated alterations predominantly in brainstem/midbrain, basal ganglia, and frontal lobe, more pronounced in the longitudinal data. Statistical analyses of the ROIs/TOIs showed a sequential pattern of structures that were assigned to previously defined neuropathological steps.
Conclusion: The combined ROI- and TOI-based DTI approach was able to map the disease stages of PSP in vivo cross-sectionally and longitudinally, lending support to DTI as a technical marker for imaging disease progression according to PSP stages. This approach might be useful as a tool for stratification of PSP patients MRI with respect to its proposed neuropathological progression in future longitudinal and autopsy-controlled studies.
Introduction
Progressive supranuclear palsy (PSP), a 4R-tauopathy characterized by subcortical tau inclusions in neurons, astrocytes, and oligodendroglia, is associated with various clinical phenotypes (Kovacs et al., 2020; Höglinger et al., 2017a; Shi et al., 2021; Stamelou et al., 2021). The “classical” phenotype, Richardson’s syndrome, presents with supranuclear gaze palsy, postural instability, pseudobulbar palsy, and frontal cognitive dysfunction (Steele et al., 1964; Litvan et al., 1996). Further PSP phenotypes have been identified which differ in terms of clinical symptoms, the intensity of their manifestation and the course of their onset (Respondek et al., 2017). The phenotypes Richardson’s syndrome (PSP-RS) and predominant Parkinsonism (PSP-P) occur most frequently (Höglinger et al., 2017a; Respondek et al., 2017). To capture the multifaceted phenotypical presentations, four functional domains (ocular motor dysfunction, postural instability, akinesia, cognitive dysfunction) are defined as clinical predictors in the current MDS diagnostic criteria (Höglinger et al., 2017a). Within each domain, three clinical features contribute different levels of diagnostic certainty (probable, possible, suggestive). Neuropathologically, sequences of PSP-related tau pathology are stratified for accumulation of neuronal, astroglial, and oligodendroglial abnormal tau (Kovacs et al., 2020) into six steps for postmortem diagnosis: step 1 is characterized by intraneuronal tau inclusions in the globus pallidus, subthalamic nucleus, and substantia nigra; step 2 by accumulation of intraneuronal tau in the midbrain tegmentum, medulla oblongata, and pontine base, and astroglial tau pathology in the striatum. Step 3 is marked by abnormal tau inclusions in the striatum, dentate nucleus, and amygdala; step 4 by increased intraneuronal tau in the frontal lobe, whereas astroglial tau accumulates in the amygdala, parietal, and temporal lobes. Step 5 is marked by intraneuronal tau pathology in the parietal and temporal lobes, and step 6 by the increase of intraneuronal tau in the occipital cortex, whereas astroglial tau accumulates in the brainstem, globus pallidus, and cerebellar dentate nucleus. This distribution pattern concept was validated in a pathology study by Briggs and colleagues (Briggs et al., 2021).
Magnetic resonance imaging (MRI) signs of PSP have been assigned certain levels of utility in diagnostic workups and in progression monitoring (Whitwell et al., 2017) and include prominent midbrain and superior cerebellar peduncle (SCP) atrophy. Volumetric/morphometric MRI studies show frontal atrophy in PSP (Whitwell et al., 2017; Quattrone et al., 2018; Shoeibi et al., 2019; Quattrone et al., 2024) and tissue loss in fronto-/mesiotemporal cortices and prefrontal regions in addition to the central midbrain and basal ganglia (Agosta et al., 2010). Diffusion tensor imaging (DTI) is a robust MRI tool for in vivo investigations of white matter (WM) neuronal tracts via fractional anisotropy (FA) mapping (Basser and Pierpaoli, 1996), combined with reconstruction of fiber pathways (Mori and van Zijl, 2002). In applications to PSP, micro-structural WM abnormalities have been shown by DTI in the superior cerebellar peduncle, corpus callosum, internal capsule, pons, and basal ganglia structures (Whitwell et al., 2017; Agosta et al., 2010; Seppi and Poewe, 2010; Knake et al., 2010; Bârlescu et al., 2021). To monitor disease progression in PSP, longitudinal studies investigated MRI abnormalities of involved structures, mostly the midbrain (Quattrone et al., 2018; Höglinger et al., 2017b); Whitwell et al., 2019; Quattrone et al., 2020; Kannenberg et al., 2021) and frontal lobe. Progression of WM alterations associated with the progression of the clinical symptoms have been shown both for PSP-RS and PSP-P (Agosta et al., 2018; Caso et al., 2018).
The question arises whether longitudinal MRI changes correlate with the neuropathological progression scheme proposed by Kovacs and colleagues (Kovacs et al., 2020). The aim of the current study was the transformation of neuropathologically-defined patterns into a neuroimaging concept for in vivo application to individual patients. First, axonal loss and myelin degradation were mapped by unbiased voxelwise FA analysis and, second, tract-wise analysis of FA values, i.e., reconstruction of tracts-of-interest (TOI), was used for tract-specific analysis of WM integrity along disease-specific tracts (Müller et al., 2007b). In a second analysis step, a hypothesis-based analysis of DTI data, guided by the neuropathological pattern (Kovacs et al., 2020), was used to map microstructural correlates of abnormal protein propagation in the brain.
Materials and methods
Subjects and patient characteristics
The analysis included two cohorts of which cohort B constituted the validation sample for cohort A’s cross-sectional results. The two cohorts of PSP patients (only including the most common subtypes PSP-RS and PSP-P) fulfilled MDS diagnostic criteria (Höglinger et al., 2017a; Table 1). All patients had a certainty level of probable PSP according to the MDS criteria (Höglinger et al., 2017a) with a combination of clinical features with high specificity. Cohort A included 78 PSP patients (mean age 70.0 ± 7.4 years, male/female = 38/40); of those, 21 PSP patients underwent a follow–up scan after 17.5 months on average (range 11.3–34.8 months), and 63 controls (mean age 68.6 ± 7.7 years, male/female = 30/33, with 17 controls undergoing a follow–up scan after 17.0 months on average (range 10.0–37.4 months). Data were recorded at 3.0T from the DESCRIBE and DANCER studies of the German Center for Neurodegenerative Diseases at twelve sites (Supplementary information I); all centers used tomographs with 3.0T from the identical vendor and had implemented the identical acquisition protocol for all subjects. Cohort B included a single-site data-set of 66 patients (70.5 ± 9.1 years, male/female = 38/28) and 44 healthy controls (68.5 ± 5.3 years, male/female = 25/19), recorded at 1.5T at Ulm University. PSP patients and controls were age- and sex-matched, there were no statistical differences at baseline as well as at follow-up (see Table 1).
Patients were characterized by disease severity [PSP Rating Scale (PSPRS); Golbe and Ohman-Strickland, 2007], clinical stage (Golbe et al., 2020), and phenotype according to MDS Diagnostic Criteria of PSP and the MAX rules (Höglinger et al., 2017a; Grimm et al., 2019). Only PSP patients who underwent MRI scans without relevant artifacts and without imaging abnormalities compromising the accurate assessment of the scans (e.g., extended vascular lesions) were considered for the study.
This retrospective study was conducted in compliance with the declaration of Helsinki and its later amendments. All subjects provided written informed consent according to institutional guidelines approved by the DZNE Ethics Committee (DESCRIBE, Deep Phenotyping of PSP) for cohort A (reference 311/14), and Ulm University Ethics Committee for cohort B (reference 279/19).
DTI acquisition
DTI scanning was performed by use of the following protocols: first, multicentric data (cohort A) were acquired at 3.0T scanners with 70 gradient directions (GD), including ten b = 0, 72 slices, 2.0 × 2.0 × 2.0 mm3 (matrix 120 × 120 × 72), TE = 88 ms, TR = 12,100 ms, 30 GD with b = 700 s/mm2; 30 GD with b = 1,000 s/mm2. This protocol was identical for all scanners used in the multicentric study due to the prospective planning of the data acquisition at all sites. Second, cohort B, recorded at a 1.5 Tesla scanner (Magnetom Symphony, Siemens Medical, Erlangen, Germany), contained 52 GD, including four b = 0, 64 slices, 2.0 × 2.0 × 2.8 mm3 (matrix 128 × 128 × 72), TE = 95 ms, TR = 8,000 ms, b = 1,000 s/mm2.
Methodological concept
A combination of data-driven and hypothesis-guided approaches was applied: first, an unbiased voxelwise approach (whole brain-based spatial statistics, WBSS) was used to map axonal loss and myelin degradation and, second, hypothesis-guided region-of-interest (ROI) and tract-of-interest (TOI) analysis was used for analysis of WM integrity along disease-specific tracts; in a second analysis step, a hypothesis-based analysis of DTI data, guided by the neuropathological pattern (Kovacs et al., 2020), was used to map microstructural correlates of abnormal protein propagation in the brain.
Data analysis: pre- and post-processing
Pre- and postprocessing was performed using the analysis software Tensor Imaging and Fiber Tracking (TIFT) (Müller et al., 2007a). All data were assessed for completeness, and—according to an established quality control protocol (Müller et al., 2014)
—data-sets with corrupted GD or motion artifacts were excluded from further analysis prior to correction of eddy current-induced geometric distortions. Data were transferred onto a 1 mm iso-grid for all further analyses (Müller et al., 2013). Non-linear spatial normalization to the Montreal Neurological Institute (MNI) stereotaxic standard space (Brett et al., 2002) was performed using study-specific templates and preserving directional information (Müller et al., 2007a; Andersson et al., 2003). This procedure was performed for baseline data (for subjects with only a baseline scan) and for baseline and follow-up data (for subjects with follow-up scans) using a prior intra-subject alignment (Menke et al., 2014). Fractional anisotropy (FA) maps were calculated from MNI-normalized DTI data, and a Gaussian smoothing filter of 8 mm full-width at half-maximum was applied to normalized individual FA maps (Unrath et al., 2010), providing a good balance between sensitivity and specificity. Finally, resulting FA maps were corrected for age (Behler et al., 2021). For both cohorts, a WBSS analysis was performed at the group level (see the following section “Data analysis: whole brain-based spatial statistics”), comparing males vs. females for PSP patients and for controls, respectively. No sex-dependency of FA values could be detected, thus, no correction for sex was performed in the subsequent analyses.
Briefly, PSP-RS and PSP-P display a similar sequential disease course in involved brain structures, although PSP-RS is more severe (Kovacs et al., 2020). For statistical reasons of sample size, no distinction was made in the following longitudinal analysis between PSP-RS and PSP-P, i.e., both subtypes were analyzed as one group. In order to acknowledge that progression rates are different in PSP-RS and PSP-P (Street et al., 2021), the cross-sectional comparison of PSP-RS vs. controls and PSP-P vs. controls was performed. In order to acknowledge that progression rates are different in PSP-RS and PSP-P (Street et al., 2021), the cross-sectional comparison of PSP-RS vs. controls and PSP-P vs. controls was performed.
Data analysis: whole brain-based spatial statistics
Statistical comparisons by Student’s t-test were performed voxel-wise for FA values to detect changes between subject groups (whole brain-based spatial statistics, WBSS). Voxels with FA values below 0.2 were not considered for statistical comparison (Kunimatsu et al., 2004). Statistical results were corrected for multiple comparisons using the false-discovery-rate (FDR) algorithm at a significance level of p < 0.05 (Genovese et al., 2002). Further reduction of the alpha error was performed by a spatial correlation algorithm that eliminated isolated voxels or small isolated groups of voxels in the size-range of the smoothing kernel, leading to a threshold cluster size of 512 voxels (Müller et al., 2013); that way, only relevant structures that show statistically significant alterations survive the correction for multiple comparisons (Müller et al., 2012).
Data analysis: tract of interest analysis
An averaged DTI data-set was calculated from the controls’ data-sets by arithmetic averaging of the MNI transformed data and was then used to identify tracts and pathways by fiber tracking (FT) for defined brain structures with a seed-to-target approach. For FT, seed and target regions were defined; all potential tracts originating in the seed region and ending in the target region (Eigenvector scalar product threshold, 0.9) of a given pathway define the corresponding TOI. As FT technique, a deterministic streamline tracking approach was used (Müller et al., 2007b). FA values within a given ROI were arithmetically averaged for each subject. The technique of tract-wise fractional anisotropy statistics (TFAS) was applied to quantify tractography results: by using the TOI, for each subject, FA values underlying the tracts were selected for arithmetic averaging. Region of interest (ROI) analysis was performed by arithmetically averaging FA values within a given ROI for each subject, considering only voxels with a FA value higher than 0.2 (Kunimatsu et al., 2004).
Student’s t-test (in case of normal distribution) or Mann-Whitney-U-test (in case of not normal distribution) were used to compare ROI- and TOI-based FA values statistically at the group level. Significance was defined at p < 0.05, FDR-corrected for multiple comparisons. Supplementary information II summarizes localization and size of the ROIs/TOIs studied.
ROI/TOI definitions
Multiple brain structures known from imaging and neuropathological studies to be affected in PSP were analyzed. The structures were divided into three phases of disease progression. The early phase (pattern 1) corresponds to steps 1-2 of the sequential neuropathological involvement described by Kovacs and colleagues (Kovacs et al., 2020); the intermediate (pattern 2) and late (pattern 3) phases correspond to steps 3-4 and steps 5-6, respectively (Table 2). Figure 1 illustrates the ROIs/TOIs which were assigned to the three PSP patterns.

Table 2. ROIs and TOIs selected for PSP phases 1-3 based on previous findings in imaging and neuropathological studies.

Figure 1. Regions-of-interest (ROIs) and tracts-of-interest (TOIs) selected for PSP patterns 1-3 based on previous findings in imaging and neuropathological studies. Sagittal projectional ROI (left panel) and TOI (right panel) representations for the three patterns corresponding to sequential involvement of brain structures during the course of PSP.
Pattern categorization and classification of PSP patients
In regions where a (highly) significant difference in FA values at the group level was identified, there is nevertheless—as an intrinsic property of the voxelwise comparison at the group level
—an overlap of the FA values of PSP and controls (see Supplementary information III for an example of overlap values). In order to take this property into account in the transition to an analysis at the individual level, a threshold value of μ−0.47σ was defined (Note: with a normal distribution, this value corresponds to a probability of 68% for controls to be above the threshold value). That way, thresholds were defined so that the specificity of each pattern reached 77.8% for cohort A and 77.3% for cohort B, i.e., 49 of 63 controls for cohort A and 34 of 44 controls for cohort B were defined as “different from PSP.” This approach was chosen in analogy to FA analyses in other neurodegenerative diseases, e.g., amyotrophic lateral sclerosis (Kassubek et al., 2014). For each ROI/TOI structure, FA values were normalized to mean and standard deviation of FA values of the controls (z-transformation) in order to obtain values comparable between different anatomical structures. Then, z-values for the ROI/TOI structures of each pattern were arithmetically averaged to obtain a pattern-specific < FA > value for each subject (pattern 1, 2, or 3). A threshold was defined so that specificity of each pattern was 77.8% for cohort A and 77.3% for cohort B, i.e., 49 of 63 controls for cohort A and 34 of 44 controls for cohort B were defined as “different from PSP.” This procedure was chosen in analogy to FA analyses in other neurodegenerative diseases, e.g., amyotrophic lateral sclerosis (Kassubek et al., 2014).
A categorization was then performed that classified each patient into one of three categorization steps (CS) according to the neuropathological staging concept (Kovacs et al., 2020), where CSα corresponds to steps 1-2, CSβ to steps 3-4, and pattern CSχ to steps 5-6 (Kovacs et al., 2020). For PSP patients, CSα is reached when z(< FA >) in pattern 1 is below the threshold of pattern 1 and z(< FA >) in pattern 2 is above the threshold of pattern 2; CSβ is reached when z(< FA >) in pattern 1 is below the threshold of pattern 1, z (< FA >) in pattern 2 is below the threshold of pattern 2, and z(< FA >) in CSχ is above the threshold of pattern 3; CSχ is reached when z(< FA >) in all patterns is below the threshold of the respective pattern. Figure 2 provides an analysis scheme.

Figure 2. Scheme for data analysis and categorization of PSP patients. Fractional anisotropy (FA) maps were calculated from DTI data, i.e., 36 parameters were calculated for each subject from combined region-of-interest (ROI) and tract-of-interest (TOI) analyses (see Supplementary information II). Parameters were sorted into 3 patterns and for each pattern (after “z-normalization”) a threshold-based decision of affectation (1 or 0) was performed. The combinations of pattern-specific decisions [the individual imaging-based categorization step (CS)] lead to a categorization into CSα, CSβ, and CSχ.
Data availability
The dataset used and analyzed during the current study will be made available by the corresponding author upon reasonable request to qualified researchers.
Results
Unbiased whole brain-based statistics
Statistical comparison of 78 PSP patients vs. 63 controls (cohort A) at baseline showed widespread alterations in multiple WM areas (Figure 3a). The comparison between the 66 patients vs. 44 controls from cohort B at baseline confirmed these results (Figure 3a). A separation into PSP-P and PSP-RS subtypes (Table 1) and separate comparison to controls demonstrated similar patterns in cohort B for both subgroups, whereas, in cohort A, PSP-RS displayed more pronounced alterations than PSP-P, especially in the putamen. In both cohorts, the PSP-RS variant was significantly more affected than PSP-P (Figure 3c and Supplementary information IV).

Figure 3. Whole brain-based spatial statistics (WBSS) of fractional anisotropy (FA) maps. (a) Exemplary slicewise representation of WBSS of FA maps of 78 PSP patients vs. 63 controls at baseline (cohort A) and 66 PSP patients vs. 44 controls at baseline (cohort B). A decrease of FA in PSP patients vs. controls is displayed in hot colors, an increase in cold colors. (b) Differences at the group level (illustrated by variation of thresholded t-values) in projectional sagittal views for cohorts A and B. (c) Subdivision of WBSS into the subtypes PSP with Richardson’s syndrome (PSP-RS) and PSP with predominant parkinsonism (PSP-P). MNI, Montreal Neurological Institute coordinate frame.
Variation of the t-value threshold showed a pattern of significant clusters that reminded of the hypothesized spreading concept of PSP: alterations with the highest t-values were observed in the midbrain and premotor WM, moderate alterations were observed in structures assigned to pattern 2, and the additional alterations at more lenient thresholds were observed in the parietal, temporal, and occipital WM (pattern 3) (Figure 3b and Supplementary information IV). This pattern was the same for cohorts A and B. A difference for PSP-P vs. controls existed in the putamen (Figure 3c)—this was bilaterally significant for cohort B (FA increase), whereas for cohort A, no significant putaminal alterations were found. However, closer inspection showed that, for cohort A, putaminal changes were present but disappeared in the analysis owing to the small cluster size. A direct voxelwise comparison of the two groups (PSP-RS and PSP-P patients) showed no significant result clusters for statistical reasons – the comparison of baseline and follow-up of patients vs. controls at baseline has the advantage of an improved signal to noise ratio compared to direct comparison of patient groups, for an analytical proof please refer to Müller and Kassubek (2024). In summary, cohort B is a validation sample for A’s cross-sectional results.
The longitudinal data comprised the subgroup of 21 patients from cohort A who had obtained follow-up scans. Baseline and follow-up data of these patients were compared with data from 63 controls at baseline by WBSS. Baseline comparisons showed a pattern of FA alterations similar to the comparison of the whole data samples of cohorts A and B at baseline, while follow-up scans showed extended alteration clusters in the longitudinal comparison compared to baseline comparison (Supplementary information IV). At a significance level of p < 0.0001, baseline data showed significant midbrain clusters (together 13,613 voxels), while longitudinal data showed extensive significant clusters in midbrain and premotor and motor areas (64,397 voxels) (Figure 4), indicating an expansion of the altered structures from baseline to follow-up.

Figure 4. Whole brain-based spatial statistics (WBSS). (a) Slicewise representation of WBSS (p < 0.05, corrected for multiple comparisons) of FA maps of 21 PSP patients vs. 63 controls at baseline and at follow-up. (b) Results at a thresholded amplitude of difference (t-values) in projectional sagittal views for baseline and at follow-up. MNI, Montreal Neurological Institute coordinate frame.
Hypothesis-guided ROI/TOI analyses
Corresponding to anatomical localization and connectivities, selected WM tracts classified to be associated with one of the aforementioned phases (refer to Table 2 and references therein) were analyzed.
From cohort A, 17 controls received a follow-up scan. Their data at baseline and follow-up were used to calculate the variability of FA parameters; the mean differences between baseline and follow-up of all ROIs/TOIs were Δ < FA > = 0.006 ± 0.007. This value constitutes the lower limit of FA-based detectability of pathological changes. As FA maps were age-corrected, this analysis was used to determine the influence of environmental noise and other measurement-specific factors: the average differences between baseline and follow-up of all ROIs/TOIs was Δ < FA > = 0.006 ± 0.007. Supplementary information V shows differences of < FA > for ROIs/TOIs between patients and controls at baseline (for cohorts A and B), subdivided into pattern 1-related structures, pattern 2-related structures, and pattern 3-related structures, respectively. Most structures (depending on the cohort) and especially the grand average over patterns 1, 2, and 3 displayed a significant reduction of < FA > in patients compared to controls; the exception was the putamen where an increase of < FA > in patients was observed.
Z-values for ROI/TOI structures were arithmetically averaged for each pattern to obtain a pattern-specific < FA > value for each subject (pattern 1, 2, or 3, respectively). Table 3I summarizes the group-averaged results for the three patterns for cohorts A and B. At group level, all patterns for cohorts A and B showed significant differences; based thereon, the patients were categorized into steps CSα, CSβ, CSχ (Table 3II). Only eight cohort A patients and one cohort B patient eluded categorization. Because both cohorts showed similar results, we confirmed that cohort B could be regarded as a validation sample for the cross-sectional results from cohort A.

Table 3. Pattern analysis results for cohorts A and B and pattern categorization into steps CSα, CSβ, CSχ.
Longitudinal analysis
Of the 21 patients in cohort A with at least one follow-up scan, seven had a substantial clinical progression (i.e., increase of the PSPRS stage). The corresponding z(< FA >) values are summarized in Table 3III and Figure 5. A significant decrease in z(< FA >) from baseline to follow-up compared to controls at baseline was observed for all patterns, and a tendency (albeit not significant) of decreased z(< FA >) was observed in patients from baseline to follow-up. Individual examples of longitudinal alterations of z(< FA >) for the three patterns are demonstrated in Supplementary information VI. As progression rates are different in PSP-RS and PSP-P (Street et al., 2021), the cross-sectional comparison of 13 PSP-RS vs. controls and eight PSP-P vs. controls (for whom a follow-up scan was available) was performed at baseline and at follow-up (although with low statistical impact due to low subject numbers), in order to analyze the extent to which corresponding brain regions were affected. The analysis PSP-RS vs. controls showed more extended significant FA decrease compared to the PSP-P vs. controls at baseline; at follow-up, both groups showed an increase of affection clusters. Due to statistical reasons (low subject numbers), FA increases in the bilateral putamen for PSP-P and in the right putamen for PSP-RS were below the detection threshold (Supplementary information VII).

Figure 5. Pattern analysis results [z (< FA >)] for cohort A for baseline and follow-up (N = 21 PSP-patients) vs. controls (N = 63). The transparent red interval indicates an error estimation obtained from baseline and follow-up controls’ scans. **p < 0.0001.
Imaging-based step categorization
Table 3II summarizes the categorization of individual subjects into PSP steps by in vivo ROI/TOI analysis for cohorts A and B, separated for cross-sectional and longitudinal analysis. The specificity of 77.8% (cohort A) and 77.3% (cohort B) was associated with a high sensitivity of 89.7% (cohort A) and 98.5% (cohort B), respectively. Of the total number of patients (cohorts A and B; N = 144), 93.8% could be assigned to imaging-based categorization steps CSα, CSβ, and CSχ. The categorization for the 21 patients with baseline and follow-up scans showed an increase to higher categorization steps from baseline to follow-up. Supplementary information VIII shows the percentage distribution. Effect sizes and sample sizes for the patterns 1, 2, and 3 were summarized in Supplementary information IX.
Categorization according to the PSP rating scale
According to the PSPRS, patients were categorized into PSPRS stages 1, 2, 3 + 4; one patient from cohort A and six patients from cohort B had PSPRS stage 4, and no patient had PSPRS stage 5. WBSS of PSPRS subgroups displayed more marked alterations corresponding to more marked FA decreases from PSPRS 1 to PSPRS 3 + 4 for both cohorts (Supplementary information X). Cohort A showed the highest z(< FA >) differences from controls [z (controls = 0.50; Table 3I] for the PSPRS 3 + 4 group in pattern 1, decreasing to pattern 3 and also decreasing to the PSPRS 1 group (thereby indicating highest FA alterations in the PSPRS 3 + 4 group for all patterns) and a decrease from pattern 1 to pattern 3 for all PSPRS groups. Cohort B had similar results; however, the FA decrease in cohort B was not as clearly pronounced with respect to the PSPRS classification as in cohort A.
Discussion
We investigated whether DTI can track progression of PSP and if this in vivo imaging-based categorization follows the proposed neuropathological PSP progression concept (Kovacs et al., 2020), in order to identify potential neuroimaging correlates of PSP-associated brain pathology at the group level. Structural connectivity analysis by DTI was considered to be an excellent candidate for investigating microstructural alterations during the course of neurodegenerative diseases (Kassubek and Müller, 2020; Kassubek et al., 2018). Corresponding to the progression of PSP pathology according to Kovacs et al. (2020), a sequential pattern of regional changes was demonstrated in two independent cohorts. A combination of analysis approaches to capture these changes has been used: first, cross-sectional unbiased WBSS was applied—the highest differences at the group level (in terms of t-values) resulted in areas involved at disease onset, whereas lower differences at the group level corresponded to later disease stages. The logic that an association with disease stages exists depending on the alteration differences is explained by simulations by Gorges and colleagues (Gorges et al., 2018): sequential involvement is a possible model for the disease alteration pattern revealed by WBSS analysis. Second, hypothesis-guided ROI/TOI analyses showed groups of ROIs/TOIs that could be summarized into imaging-based categorization steps CSα, CSβ, CSχ which displayed sequential involvement of disease stage-related patterns at the group level for the PSP patients. Longitudinally, both the WBSS and the ROI/TOI analyses demonstrated an increase of the brain alteration pattern after a 17-month observation period. A sophisticated model (which already has been applied in the analysis of atrophy patterns in PSP) is the “Subtype and Stage Inference” (SuStaIn) model (Saito et al., 2022); this model was able to identify distinct temporal progression patterns of brain atrophy for corticobasal syndrome and PSP-RS with high reproducibility and high accuracy based on cross-sectional structural brain MRI data. Given that PSP-RS and PSP-P displayed a similar sequential disease course for involved brain structures (Kovacs et al., 2020), the two subtypes were analyzed as one group. Nevertheless, the cross-sectional comparison of PSP-RS vs. controls and PSP-P vs. controls showed fewer alterations in the PSP-P group, albeit in the same brain regions.
Longitudinal scans of controls were used to estimate the variability of FA at the group level; this estimation can serve as a “threshold” for alterations detected in PSP patients, which were significantly higher than in controls so that FA alterations in patients at the group level can be considered to be pathological. However, due to the diverse sources of variability in in vivo DTI scans, we reduced the six neuropathological steps to three imaging-based categorization steps (CSα, CSβ, CSχ) so that a higher degree of stability in the resulting patterns could be obtained. Furthermore, the large number of ROIs/TOIs helped to overcome DTI variability and the inhomogeneity in clinical presentation of PSP. One strength of this study was that cross-check could be performed because two independently obtained patient and control cohorts were evaluated, and the results from both cohorts displayed a high degree of agreement and consistency.
In this study, longitudinal outcomes were used as a validation to identify disease-related structures and their involvement in PSP over the course of the disease; accurate modeling of disease progression would require a significant higher amount of cross-sectional and longitudinal DTI (and structural MRI) data. Approaches to longitudinal atrophy of brain structures in PSP (Scotton et al., 2023) and tractography-associated changes in PSP (Costa et al., 2025) address disease progression and can be modeled as temporal progression patterns by current approaches such as “Subtype and Stage Inference” (SuStaIn) (Saito et al., 2022).
The 4R-tauopathy PSP is characterized by subcortical tau inclusions in neurons, astrocytes, and oligodendroglia leading to the concept of PSP-associated sequences of tau progression (Kovacs et al., 2020). WM damage measured using diffusion tensor imaging (DTI) is also a consistent feature of PSP and is thought to result from a tau-driven neurodegenerative process partially independent of gray matter degeneration (Costa et al., 2025; Armstrong, 2013). DTI and its resulting metrics address a physical parameter (diffusion) and, thus, indirectly measure a grade of integrity especially in WM tracts. Alterations in DTI metrics could result from fiber loss, synaptic loss, or loss of fiber integrity; while FA is sensitive to microstructural changes, it does not indicate a specific type of lesion. On the other hand, AD tends to be strongly affected by axonal injury whereas RD is sensitive to white matter damage due to demyelination and less to changes in the axonal density or size (Song et al., 2002; Winklewski et al., 2018; Wang et al., 2015). That way, the advantage of application of DTI as a non-invasive in vivo technique is shadowed by the fact that DTI metrics do only indirectly measure tau pathology and, thus, could be affected by further diffusion-influencing factors. Diffusion properties represented by the voxelwise diffusion ellipsoid could be parameterized by a series of parameters which are essentially divided into voxelwise parameters (e.g., FA, mean diffusivity, radial and axial diffusivity) and tractwise parameters (resulting from the feasibility of tract reconstruction; Schilling et al., 2021). In the first case, the voxelwise parameters result from a parameterization by different combination of the Eigenvalues of the ellipsoid; therefore, in order to avoid redundancy, it is crucial to choose one of these parameters when mapping disease-related diffusion alterations in vivo. In the second case, fiber tractography bundle segmentation depends on scanner effects, vendor effects, acquisition resolution, diffusion sampling scheme, and diffusion sensitization. FA is particularly suitable because of its intrinsic normalization to values between 0 and 1, which then indicate the strength of the directionality of the diffusion.
This study aimed to investigate regional DTI changes longitudinally corresponding to neuropathological stages of PSP-associated sequences of tau progression (Kovacs et al., 2020). Regions with neuropathologically proven tau inclusions or degenerative changes served as a basis for ROIs/TOIs. Some of these involved gray matter, others WM regions and tracts. Selection was based on the involvement reported in previous neuropathology or imaging studies and/or the relationship to related structures that displayed involvement (mostly WM regions and tracts) (Table 2 and Figures 3, 4). Prior research showed that brain connectivity in PSP is linked to tau deposition (Franzmeier et al., 2022) in support of DTI findings of WM damage; moreover, the possible contribution of these anatomical regions to PSP symptomatology was a rationale for their inclusion (Kovacs et al., 2020; Dickson et al., 2007). In that context, the findings in the putamen need to be commented on, given that FA decrease is caused by axonal loss or myelin degradation. Still, FA elevation in the putamen remains an unresolved challenge, as gray matter has a more isotropic diffusion profile, making DTI metrics such as FA less meaningful due to its complex cytoarchitecture. DTI metrics are influenced by iron accumulation which occurs physiologically with aging so that FA levels in the putamen and caudate have been shown to be higher in healthy older subjects (Pfefferbaum et al., 2010); however, the age correction approach (Behler et al., 2021) used in our study accounts for both negative and positive age-associated FA changes (which in parallel could occur in different brain regions) due to voxelwise correction. Structural neuroimaging of PSP already showed degeneration throughout a network of gray matter structures, including midbrain, thalamus, subthalamic nucleus, globus pallidus, striatum and frontal cortex, and especially for PSP-RS and PSP-P in the putamen (Whitwell et al., 2020; Spiegel et al., 2025); furthermore, among further tracts, the corticostriatal tract showed progression in the course of disease (Costa et al., 2025). However, as tau pathology is not the only underlying pathological mechanisms (synaptic dysfunction and loss may precede tau accumulation have to be considered), the DTI-based results of the current study have to be regarded as an indicator for putaminal involvement and require additional confirmation by future studies. The FA increase pattern detected in the putamen most probably resulted from pathological processes that modify tissue integrity, e.g., neuronal remodeling and loss, secondary astrocytosis, loss of specific fiber tracts, alterations in fiber composition, or changes in cell permeability (Rosas et al., 2006; Douaud et al., 2009; Douaud et al., 2011). Therefore, a specifically pronounced loss of fiber integrity/directionality might result in an apparent FA increase.
Our findings must be considered within the context of several limitations. The challenges of in vivo DTI in neurodegeneration are manifold: (1) different sources of variability may contribute to the diffusion data of subjects, e.g., scanner variability, environmental noise, individual subject factors, which cannot be singled out but may bias the data (Müller et al., 2013). We are aware that there was an imbalance of participants (PSP and controls) of cohort A across the scanners (Supplementary information I) which may impact the group comparisons. However, given that the identical acquisition protocol had been implemented on all scanners (all of the same vendor) and we have additionally tested for the inter subject-homogeneity of whole-brain FA values across sites/scanners (Supplementary information Figure I-2), it seems safe to assume that the inter-site-variability of the scans is lower than the inter-individual variability of DTI metrics of participants. (2) As in other DTI data-sets, our study was retrospective; thus, scan repetitions (in the event of variable measurement results) were not possible. A major limitation of our in vivo neuroimaging approach was the lack of autopsy-confirmed data. Thus, the ROI- and TOI-based analysis can only provide a plausible surrogate in vivo “staging” pattern for the presumed PSP pathology in the PSP cohorts investigated here. A further limitation was that an overlap in FA values existed between controls and patients. In future studies, this overlap could be compensated by a large number of controls and subsequent exclusion of control data that lay outside the interquartile range. Due to the intrinsic property of DTI-based FA analysis (Supplementary information III), in order to account for the (slight) overlap of FA values of PSP patients and controls in pattern 3, the decision was made to categorize to CSα when z (< FA >) in pattern 1 was below the threshold of pattern 1 and z (< FA >) in pattern 2 was above the threshold of pattern 2, regardless of z (< FA >) in pattern 3. This constellation was only present in a few PSP patients with z (< FA >) in pattern 3 only slightly above the threshold of pattern 3. Given the character of this ex post facto approach, the sequential results in cohort A could not be validated in cohort B because no longitudinal data were available for cohort B. Furthermore, cohort B patients showed a shorter disease duration compared to cohort A. (3) In the current data analysis, we have chosen FA as the appropriate parameter to demonstrate the utility of diffusion-weighted imaging/DTI for the detection of disease-related changes and progression in the brains of PSP patients in terms of imaging-associated neuropathological progression in vivo. Further studies might be performed in the future which might substitute FA by other diffusion metrics (e.g., mean diffusivity, radial diffusivity, and axial diffusivity) or tractwise parameters (e.g., fiber density) and could also include volumetric alterations to finally enable a comprehensive AI-based neuroimaging categorization for PSP (Volkmann et al., 2025).
Conclusion
In conclusion, the novel ROI/TOI-based technique employed here enabled us to analyze specific PSP-associated brain structures categorized into distinct disease steps, thereby rendering it possible to image the proposed neuropathological PSP stages. This procedure constitutes an in vivo confirmation of the neuropathological propagation concept of disease stages in PSP. In future applications, it might broaden the spectrum of potential non-invasive surrogate markers as a neuroimaging-based read-out for PSP studies at the group level within a clinical context.
Data availability statement
The data analyzed in this study is subject to the following licenses/restrictions: The dataset used and analyzed during the current study will be made available by the corresponding author upon reasonable request to qualified researchers. Requests to access these datasets should be directed to amFuLmthc3N1YmVrQHVuaS11bG0uZGU=.
Ethics statement
The studies involving humans were approved by DZNE Ethics Committee, Ulm University Ethics Committee. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
LB: Data curation, Formal Analysis, Writing – original draft. GH: Data curation, Writing – review & editing. HV: Formal Analysis, Writing – review & editing. AL: Conceptualization, Writing – review & editing. KD: Conceptualization, Formal Analysis, Writing – review & editing. HB: Conceptualization, Writing – review & editing. H-PM: Conceptualization, Formal Analysis, Methodology, Visualization, Writing – original draft. JK: Conceptualization, Formal Analysis, Methodology, Supervision, Writing – original draft. Members of the DESCRIBE-PSP Study Group (Appendix).
Funding
The author(s) declare that no financial support was received for the research and/or publication of this article.
Acknowledgments
Data were provided by the study groups of the DESCRIBE studies of the Clinical Research of the German Center for Neurogenerative Diseases (DZNE).
Conflict of interest
Günter Höglinger has ongoing research collaborations with Roche, UCB, Abbvie; serves as a consultant for Abbvie, Alzprotect, Amylyx, Aprinoia, Asceneuron, Bayer, Bial, Biogen, Biohaven, Epidarex, Ferrer, Kyowa Kirin, Lundbeck, Novartis, Retrotope, Roche, Sanofi, Servier, Takeda, Teva, UCB; received honoraria for scientific presentations from Abbvie, Bayer, Bial, Biogen, Bristol Myers Squibb, Esteve, Kyowa Kirin, Pfizer, Roche, Teva, UCB, Zambon. holds a patent “Treatment of Synucleinopathies” United States Patent No. US 10,918,628 B2/European Patent Patent No.: EP 17 787 904.6-1109/3 525 788; received publication royalties from Academic Press, Kohlhammer, and Thieme. Günter Höglinger was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy within the framework of the Munich Cluster for Systems Neurology (EXC 2145 SyNergy – ID 390857198); European Joint Programme on Rare Diseases (Improve-PSP); Petermax-Müller Foundation (Etiology and Therapy of Synucleinopathies and Tauopathies). Jan Kassubek has received honoraria or consultation fees from AbbVie, Bial, Biogen, Desitin, Esteve, Licher MT, Medtronic, NeuroDerm, Novartis, STADA, UCB Pharma, and Zambon; in addition, he is Specialty Chief Editor for Frontiers in Neurology (section Applied Neuroimaging) and Associate Editor (Neurology) for Therapeutic Advances in Chronic Disease.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2025.1569302/full#supplementary-material
References
Agosta, F., Caso, F., Ječmenica-Lukić, M., Petrović, I., Valsasina, P., Meani, A., et al. (2018). Tracking brain damage in progressive supranuclear palsy: A longitudinal MRI study. J. Neurol. Neurosurg. Psychiatry 89, 696–701. doi: 10.1136/jnnp-2017-317443
Agosta, F., Kostić, V., Galantucci, S., Mesaros, S., Svetel, M., Pagani, E., et al. (2010). The in vivo distribution of brain tissue loss in Richardson’s syndrome and PSP-parkinsonism: A VBM-DARTEL study. Eur. J. Neurosci. 32, 640–647. doi: 10.1111/j.1460-9568.2010.07304.x
Albrecht, F., Bisenius, S., Neumann, J., Whitwell, J., and Schroeter, M. (2019). Atrophy in midbrain & cerebral/cerebellar pedunculi is characteristic for progressive supranuclear palsy - A double-validation whole-brain meta-analysis. Neuroimage Clin. 22:101722. doi: 10.1016/j.nicl.2019.101722
Andersson, J., Skare, S., and Ashburner, J. (2003). How to correct susceptibility distortions in spin-echo echo-planar images: Application to diffusion tensor imaging. Neuroimage 20, 870–888. doi: 10.1016/S1053-8119(03)00336-7
Archer, D., Mitchell, T., Burciu, R., Yang, J., Nigro, S., Quattrone, A., et al. (2020). Magnetic resonance imaging and neurofilament light in the differentiation of parkinsonism. Mov. Disord. 35, 1388–1395. doi: 10.1002/mds.28060
Armstrong, M. (2018). Progressive supranuclear palsy: An update. Curr. Neurol. Neurosci. Rep. 18:12. doi: 10.1007/s11910-018-0819-5
Armstrong, R. (2013). White matter pathology in progressive supranuclear palsy (PSP): A quantitative study of 8 cases. Clin. Neuropathol. 32, 399–405. doi: 10.5414/NP300608
Bârlescu, L., Müller, H., Uttner, I., Ludolph, A., Pinkhardt, E., Huppertz, H., et al. (2021). Segmental alterations of the corpus callosum in progressive supranuclear palsy: A multiparametric magnetic resonance imaging study. Front. Aging Neurosci. 13:720634. doi: 10.3389/fnagi.2021.720634
Basser, P., and Pierpaoli, C. (1996). Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. J. Magn. Reson. B 111, 209–219. doi: 10.1006/jmrb.1996.0086
Behler, A., Kassubek, J., and Müller, H. (2021). Age-related alterations in DTI metrics in the human brain-consequences for age correction. Front. Aging Neurosci. 13:682109. doi: 10.3389/fnagi.2021.682109
Brenneis, C., Seppi, K., Schocke, M., Benke, T., Wenning, G., and Poewe, W. (2004). Voxel based morphometry reveals a distinct pattern of frontal atrophy in progressive supranuclear palsy. J Neurol Neurosurg. Psychiatry 75, 246–249. doi: 10.1136/jnnp.2003.015297
Brett, M., Johnsrude, I., and Owen, A. (2002). The problem of functional localization in the human brain. Nat. Rev. Neurosci. 3, 243–249. doi: 10.1038/nrn756
Briggs, M., Allinson, K., Malpetti, M., Spillantini, M., Rowe, J., and Kaalund, S. (2021). Validation of the new pathology staging system for progressive supranuclear palsy. Acta Neuropathol. 141, 787–789. doi: 10.1007/s00401-021-02298-z
Burciu, R., Ofori, E., Shukla, P., Planetta, P., Snyder, A., Li, H., et al. (2015). Distinct patterns of brain activity in progressive supranuclear palsy and Parkinson’s disease. Mov. Disord. 30, 1248–1258. doi: 10.1002/mds.26294
Canu, E., Agosta, F., Baglio, F., Galantucci, S., Nemni, R., and Filippi, M. (2011). Diffusion tensor magnetic resonance imaging tractography in progressive supranuclear palsy. Mov. Disord. 26, 1752–1755. doi: 10.1002/mds.23739
Caso, F., Agosta, F., Ječmenica-Lukić, M., Petrović, I., Meani, A., Kostic, V., et al. (2018). Progression of white matter damage in progressive supranuclear palsy with predominant parkinsonism. Parkinsonism Relat. Disord. 49, 95–99. doi: 10.1016/j.parkreldis.2018.01.001
Chougar, L., Pyatigorskaya, N., and Lehéricy, S. (2021). Update on neuroimaging for categorization of Parkinson’s disease and atypical parkinsonism. Curr. Opin. Neurol. 34, 514–524. doi: 10.1097/WCO.0000000000000957
Clark, H., Tosakulwong, N., Weigand, S., Ali, F., Botha, H., Pham, N., et al. (2021). Gray and white matter correlates of dysphagia in progressive supranuclear palsy. Mov. Disord. 36, 2669–2675. doi: 10.1002/mds.28731
Costa, F., Gatto, R., Pham, N., Ali, F., Clark, H., Stierwalt, J., et al. (2025). Longitudinal assessment of white matter alterations in progressive supranuclear palsy variants using diffusion tractography. Parkinsonism Relat. Disord. 132:107272. doi: 10.1016/j.parkreldis.2025.107272
Dickson, D., Rademakers, R., and Hutton, M. (2007). Progressive supranuclear palsy: Pathology and genetics. Brain Pathol. 17, 74–82. doi: 10.1111/j.1750-3639.2007.00054.x
Douaud, G., Behrens, T., Poupon, C., Cointepas, Y., Jbabdi, S., Gaura, V., et al. (2009). In vivo evidence for the selective subcortical degeneration in Huntington’s disease. Neuroimage 46, 958–966. doi: 10.1016/j.neuroimage.2009.03.044
Douaud, G., Jbabdi, S., Behrens, T., Menke, R., Gass, A., Monsch, A., et al. (2011). DTI measures in crossing-fibre areas: Increased diffusion anisotropy reveals early white matter alteration in MCI and mild Alzheimer’s disease. Neuroimage 55, 880–890. doi: 10.1016/j.neuroimage.2010.12.008
Franzmeier, N., Brendel, M., Beyer, L., Slemann, L., Kovacs, G., Arzberger, T., et al. (2022). Tau deposition patterns are associated with functional connectivity in primary tauopathies. Nat. Commun. 13:1362. doi: 10.1038/s41467-022-28896-3
Genovese, C., Lazar, N., and Nichols, T. (2002). Thresholding of statistical maps in functional neuroimaging using the false discovery rate. Neuroimage 15, 870–878. doi: 10.1006/nimg.2001.1037
Golbe, L., and Ohman-Strickland, P. A. (2007). A clinical rating scale for progressive supranuclear palsy. Brain 130, 1552–1565. doi: 10.1093/brain/awm032
Golbe, L., Ohman-Strickland, P., Beisser, E., and Elghoul, F. T. A. (2020). Convenient prognostic tool and staging system for progressive supranuclear palsy. Mov. Disord. Clin. Pract. 7, 664–671. doi: 10.1002/mdc3.13010
Gorges, M., Del Tredici, K., Dreyhaupt, J., Braak, H., Ludolph, A., Müller, H., et al. (2018). Corticoefferent pathology distribution in amyotrophic lateral sclerosis: In vivo evidence from a meta-analysis of diffusion tensor imaging data. Sci. Rep. 8:15389. doi: 10.1038/s41598-018-33830-z
Gorges, M., Maier, M., Rosskopf, J., Vintonyak, O., Pinkhardt, E., Ludolph, A., et al. (2017). Regional microstructural damage and patterns of eye movement impairment: A DTI and video-oculography study in neurodegenerative parkinsonian syndromes. J. Neurol. 264, 1919–1928. doi: 10.1007/s00415-017-8579-8
Grimm, M., Respondek, G., Stamelou, M., Arzberger, T., Ferguson, L., Gelpi, E., et al. (2019). How to apply the movement disorder society criteria for diagnosis of progressive supranuclear palsy. Mov. Disord. 34, 1228–1232. doi: 10.1002/mds.27666
Hofer, S., and Frahm, J. (2006). Topography of the human corpus callosum revisited–comprehensive fiber tractography using diffusion tensor magnetic resonance imaging. Neuroimage 32, 989–994. doi: 10.1016/j.neuroimage.2006.05.044
Höglinger, G., Respondek, G., Stamelou, M., Kurz, C., Josephs, K., Lang, A., et al. (2017a). Clinical diagnosis of progressive supranuclear palsy: The movement disorder society criteria. Mov. Disord. 32, 853–864. doi: 10.1002/mds.26987
Höglinger, G., Schöpe, J., Stamelou, M., Kassubek, J., Del Ser, T., Boxer, A., et al. (2017b). Longitudinal magnetic resonance imaging in progressive supranuclear palsy: A new combined score for clinical trials. Mov. Disord. 32, 842–852. doi: 10.1002/mds.26973
Josephs, K., Xia, R., Mandrekar, J., Gunter, J., Senjem, M., Jack, C., et al. (2013). Modeling trajectories of regional volume loss in progressive supranuclear palsy. Mov. Disord. 28, 1117–1124. doi: 10.1002/mds.25437
Kannenberg, S., Caspers, J., Dinkelbach, L., Moldovan, A., Ferrea, S., Südmeyer, M., et al. (2021). Investigating the 1-year decline in midbrain-to-pons ratio in the differential diagnosis of PSP and IPD. J. Neurol. 268, 1526–1532. doi: 10.1007/s00415-020-10327-2
Kassubek, J., and Müller, H. (2020). Advanced neuroimaging approaches in amyotrophic lateral sclerosis: Refining the clinical diagnosis. Expert. Rev. Neurother. 20, 237–249. doi: 10.1080/14737175.2020.1715798
Kassubek, J., Müller, H., Del Tredici, K., Brettschneider, J., Pinkhardt, E., Lulé, D., et al. (2014). Diffusion tensor imaging analysis of sequential spreading of disease in amyotrophic lateral sclerosis confirms patterns of TDP-43 pathology. Brain 137, 1733–1740. doi: 10.1093/brain/awu090
Kassubek, J., Müller, H., Del Tredici, K., Lulé, D., Gorges, M., Braak, H., et al. (2018). Imaging the pathoanatomy of amyotrophic lateral sclerosis in vivo: Targeting a propagation-based biological marker. J. Neurol. Neurosurg. Psychiatry 89, 374–381. doi: 10.1136/jnnp-2017-316365
Knake, S., Belke, M., Menzler, K., Pilatus, U., Eggert, K., Oertel, W., et al. (2010). In vivo demonstration of microstructural brain pathology in progressive supranuclear palsy: A DTI study using TBSS. Mov. Disord. 25, 1232–1238. doi: 10.1002/mds.23054
Kovacs, G., Lukic, M., Irwin, D., Arzberger, T., Respondek, G., Lee, E., et al. (2020). Distribution patterns of tau pathology in progressive supranuclear palsy. Acta Neuropathol. 140, 99–119. doi: 10.1007/s00401-020-02158-2
Kunimatsu, A., Aoki, S., Masutani, Y., Abe, O., Hayashi, N., Mori, H., et al. (2004). The optimal trackability threshold of fractional anisotropy for diffusion tensor tractography of the corticospinal tract. Magn. Reson. Med. Sci. 3, 11–17. doi: 10.2463/mrms.3.11
Lee, J., and Lee, M. (2019). Brain iron accumulation in atypical Parkinsonian syndromes: In vivo MRI evidences for distinctive patterns. Front. Neurol. 10:74. doi: 10.3389/fneur.2019.00074
Lenka, A., Pasha, S., Mangalore, S., George, L., Jhunjhunwala, K., Bagepally, B., et al. (2017). Role of corpus callosum volumetry in differentiating the subtypes of progressive supranuclear palsy and early Parkinson’s disease. Mov. Disord. Clin. Pract. 4, 552–558. doi: 10.1002/mdc3.12473
Litvan, I., Agid, Y., Calne, D., Campbell, G., Dubois, B., Duvoisin, R., et al. (1996). Clinical research criteria for the diagnosis of progressive supranuclear palsy (Steele-Richardson-Olszewski syndrome): Report of the NINDS-SPSP international workshop. Neurology 47, 1–9. doi: 10.1212/wnl.47.1.1
Menke, R., Körner, S., Filippini, N., Douaud, G., Knight, S., Talbot, K., et al. (2014). Widespread grey matter pathology dominates the longitudinal cerebral MRI and clinical landscape of amyotrophic lateral sclerosis. Brain 137, 2546–2555. doi: 10.1093/brain/awu162
Mori, S., and van Zijl, P. (2002). Fiber tracking: Principles and strategies - a technical review. NMR Biomed. 15, 468–480. doi: 10.1002/nbm.781
Müller, H., Unrath, A., Sperfeld, A., Ludolph, A., Riecker, A., and Kassubek, J. (2007b). Diffusion tensor imaging and tractwise fractional anisotropy statistics: Quantitative analysis in white matter pathology. Biomed. Eng. Online 6:42. doi: 10.1186/1475-925X-6-42
Müller, H., and Kassubek, J. (2024). Toward diffusion tensor imaging as a biomarker in neurodegenerative diseases: Technical considerations to optimize recordings and data processing. Front. Hum. Neurosci. 18:1378896. doi: 10.3389/fnhum.2024.1378896
Müller, H., Grön, G., Sprengelmeyer, R., Kassubek, J., Ludolph, A., Hobbs, N., et al. (2013). Evaluating multicenter DTI data in Huntington’s disease on site specific effects: An ex post facto approach. Neuroimage Clin. 2, 161–167. doi: 10.1016/j.nicl.2012.12.005
Müller, H., Kassubek, J., Grön, G., Sprengelmeyer, R., Ludolph, A., Klöppel, S., et al. (2014). Impact of the control for corrupted diffusion tensor imaging data in comparisons at the group level: An application in Huntington disease. Biomed. Eng. Online 13:128. doi: 10.1186/1475-925X-13-128
Müller, H., Unrath, A., Huppertz, H., Ludolph, A., and Kassubek, J. (2012). Neuroanatomical patterns of cerebral white matter involvement in different motor neuron diseases as studied by diffusion tensor imaging analysis. Amyotroph. Lateral Scler. 13, 254–264. doi: 10.3109/17482968.2011.653571
Müller, H., Unrath, A., Ludolph, A., and Kassubek, J. (2007a). Preservation of diffusion tensor properties during spatial normalization by use of tensor imaging and fibre tracking on a normal brain database. Phys. Med. Biol. 52, N99–N109. doi: 10.1088/0031-9155/52/6/N01
Nguyen, T., Cheng, J., Chen, Y., Lin, Y., Tsai, C., Lu, C., et al. (2021). Fixel-based analysis of white matter degeneration in patients with progressive supranuclear palsy or multiple system atrophy, as compared to Parkinson’s disease. Front. Aging Neurosci. 13:625874. doi: 10.3389/fnagi.2021.625874
Nicastro, N., Rodriguez, P., Malpetti, M., Bevan-Jones, W., Simon Jones, P., Passamonti, L., et al. (2020). 18F-AV1451 PET imaging and multimodal MRI changes in progressive supranuclear palsy. J. Neurol. 267, 341–349. doi: 10.1007/s00415-019-09566-9
Padovani, A., Borroni, B., Brambati, S., Agosti, C., Broli, M., Alonso, R., et al. (2006). Diffusion tensor imaging and voxel based morphometry study in early progressive supranuclear palsy. J. Neurol. Neurosurg. Psychiatry 77, 457–463. doi: 10.1136/jnnp.2005.075713
Paviour, D., Price, S., Jahanshahi, M., Lees, A., and Fox, N. (2006). Longitudinal MRI in progressive supranuclear palsy and multiple system atrophy: Rates and regions of atrophy. Brain 129, 1040–1049. doi: 10.1093/brain/awl021
Pfefferbaum, A., Adalsteinsson, E., Rohlfing, T., and Sullivan, E. (2010). Diffusion tensor imaging of deep gray matter brain structures: Effects of age and iron concentration. Neurobiol. Aging 31, 482–493. doi: 10.1016/j.neurobiolaging.2008.04.013
Potrusil, T., Krismer, F., Beliveau, V., Seppi, K., Müller, C., Troger, F., et al. (2020). Diagnostic potential of automated tractography in progressive supranuclear palsy variants. Parkinsonism Relat. Disord. 72, 65–71. doi: 10.1016/j.parkreldis.2020.02.007
Quattrone, A., Bianco, M., Antonini, A., Vaillancourt, D., Seppi, K., Ceravolo, R., et al. (2022). Development and validation of automated magnetic resonance parkinsonism index 2.0 to distinguish progressive supranuclear palsy-parkinsonism from Parkinson’s disease. Mov. Disord. 37, 1272–1281. doi: 10.1002/mds.28992
Quattrone, A., Franzmeier, N., Huppertz, H., Klietz, M., Roemer, S., Boxer, A., et al. (2024). Magnetic resonance imaging measures to track atrophy progression in progressive supranuclear palsy in clinical trials. Mov. Disord. 39, 1329–1342. doi: 10.1002/mds.29866
Quattrone, A., Morelli, M., Nigro, S., Quattrone, A., Vescio, B., Arabia, G., et al. (2018). A new MR imaging index for differentiation of progressive supranuclear palsy-parkinsonism from Parkinson’s disease. Parkinsonism Relat. Disord. 54, 3–8. doi: 10.1016/j.parkreldis.2018.07.016
Quattrone, A., Morelli, M., Quattrone, A., Vescio, B., Nigro, S., Arabia, G., et al. (2020). Magnetic Resonance Parkinsonism Index for evaluating disease progression rate in progressive supranuclear palsy: A longitudinal 2-year study. Parkinsonism Relat. Disord. 72, 1–6. doi: 10.1016/j.parkreldis.2020.01.019
Rau, A., Jost, W., Demerath, T., Kellner, E., Reisert, M., and Urbach, H. (2022). Diffusion microstructure imaging in progressive supranuclear palsy: Reduced axonal volumes in the superior cerebellar peduncles, dentato-rubro-thalamic tracts, ventromedial thalami, and frontomesial white matter. Cereb. Cortex 32, 5628–5636. doi: 10.1093/cercor/bhac041
Respondek, G., Kurz, C., Arzberger, T., Compta, Y., Englund, E., Ferguson, L., et al. (2017). Which ante mortem clinical features predict progressive supranuclear palsy pathology? Mov. Disord. 32, 995–1005. doi: 10.1002/mds.27034
Robinson, J., Yan, N., Caswell, C., Xie, S., Suh, E., Van Deerlin, V., et al. (2020). Primary tau pathology, not copathology, correlates with clinical symptoms in PSP and CBD. J. Neuropathol. Exp. Neurol. 79, 296–304. doi: 10.1093/jnen/nlz141
Rosas, H., Tuch, D., Hevelone, N., Zaleta, A., Vangel, M., Hersch, S., et al. (2006). Diffusion tensor imaging in presymptomatic and early Huntington’s disease: Selective white matter pathology and its relationship to clinical measures. Mov. Disord. 21, 1317–1325. doi: 10.1002/mds.20979
Rosskopf, J., Müller, H., Huppertz, H., Ludolph, A., Pinkhardt, E., and Kassubek, J. (2014). Frontal corpus callosum alterations in progressive supranuclear palsy but not in Parkinson’s disease. Neurodegener. Dis. 14, 184–193. doi: 10.1159/000367693
Saito, Y., Kamagata, K., Wijeratne, P., Andica, C., Uchida, W., Takabayashi, K., et al. (2022). Temporal progression patterns of brain atrophy in corticobasal syndrome and progressive supranuclear palsy revealed by subtype and stage inference (SuStaIn). Front. Neurol. 13:814768. doi: 10.3389/fneur.2022.814768
Sakae, N., Josephs, K., Litvan, I., Murray, M., Duara, R., Uitti, R., et al. (2019). Neuropathologic basis of frontotemporal dementia in progressive supranuclear palsy. Mov. Disord. 34, 1655–1662. doi: 10.1002/mds.27816
Sakurai, K., Tokumaru, A., Shimoji, K., Murayama, S., Kanemaru, K., Morimoto, S., et al. (2017). Beyond the midbrain atrophy: Wide spectrum of structural MRI finding in cases of pathologically proven progressive supranuclear palsy. Neuroradiology 59, 431–443. doi: 10.1007/s00234-017-1812-4
Sarica, A., Quattrone, A., Mechelli, A., Vaccaro, M., Morelli, M., and Quattrone, A. (2021). Corticospinal tract abnormalities and ventricular dilatation: A transdiagnostic comparative tractography study. Neuroimage Clin. 32:102862. doi: 10.1016/j.nicl.2021.102862
Schilling, K., Tax, C., Rheault, F., Hansen, C., Yang, Q., Yeh, F., et al. (2021). Fiber tractography bundle segmentation depends on scanner effects, vendor effects, acquisition resolution, diffusion sampling scheme, diffusion sensitization, and bundle segmentation workflow. Neuroimage 242:118451. doi: 10.1016/j.neuroimage.2021.118451
Scotton, W., Shand, C., Todd, E., Bocchetta, M., Cash, D., VandeVrede, L., et al. (2023). Uncovering spatiotemporal patterns of atrophy in progressive supranuclear palsy using unsupervised machine learning. Brain Commun. 5:fcad048. doi: 10.1093/braincomms/fcad048
Seppi, K., and Poewe, W. (2010). Brain magnetic resonance imaging techniques in the diagnosis of parkinsonian syndromes. Neuroimaging Clin. N. Am. 20, 29–55. doi: 10.1016/j.nic.2009.08.016
Shi, Y., Zhang, W., Yang, Y., Murzin, A., Falcon, B., Kotecha, A., et al. (2021). Structure-based classification of tauopathies. Nature 598, 359–363. doi: 10.1038/s41586-021-03911-7
Shoeibi, A., Olfati, N., and Litvan, I. (2019). Frontrunner in translation: Progressive supranuclear palsy. Front. Neurol. 10:1125. doi: 10.3389/fneur.2019.01125
Song, S., Sun, S., Ramsbottom, M., Chang, C., Russell, J., and Cross, A. (2002). Dysmyelination revealed through MRI as increased radial (but unchanged axial) diffusion of water. Neuroimage 17, 1429–1436. doi: 10.1006/nimg.2002.1267
Spiegel, C., Marotta, C., Bertram, K., Vivash, L., and Harding, I. (2025). Brainstem and cerebellar radiological findings in progressive supranuclear palsy. Brain Commun. 7:fcaf051. doi: 10.1093/braincomms/fcaf051
Spotorno, N., Hall, S., Irwin, D., Rumetshofer, T., Acosta-Cabronero, J., Deik, A., et al. (2019). Diffusion tensor MRI to distinguish progressive supranuclear palsy from α-synucleinopathies. Radiology 293, 646–653. doi: 10.1148/radiol.2019190406
Stamelou, M., Knake, S., Oertel, W., and Höglinger, G. (2011). Magnetic resonance imaging in progressive supranuclear palsy. J. Neurol. 258, 549–558. doi: 10.1007/s00415-010-5865-0
Stamelou, M., Respondek, G., Giagkou, N., Whitwell, J., Kovacs, G., and Höglinger, G. (2021). Evolving concepts in progressive supranuclear palsy and other 4-repeat tauopathies. Nat. Rev. Neurol. 17, 601–620. doi: 10.1038/s41582-021-00541-5
Steele, J., Richardson, J., and Olszewski, J. (1964). Progressive Supranuclear Palsy. A Heterogeneous degeneration involving the brain stem, basal ganglia and cerebellum with vertical gaze and pseudobulbar palsy, nuchal dystonia and dementia. Arch. Neurol. 10, 333–359. doi: 10.1001/archneur.1964.00460160003001
Street, D., Malpetti, M., Rittman, T., Ghosh, B., Murley, A., Coyle-Gilchrist, I., et al. (2021). Clinical progression of progressive supranuclear palsy: Impact of trials bias and phenotype variants. Brain Commun. 3:fcab206. doi: 10.1093/braincomms/fcab206
Unrath, A., Müller, H., Riecker, A., Ludolph, A., Sperfeld, A., and Kassubek, J. (2010). Whole brain-based analysis of regional white matter tract alterations in rare motor neuron diseases by diffusion tensor imaging. Hum. Brain Mapp. 31, 1727–1740. doi: 10.1002/hbm.20971
Volkmann, H., Höglinger, G., Grön, G., Bârlescu, L., and Müller, H. (2025). MRI classification of progressive supranuclear palsy, Parkinson disease and controls using deep learning and machine learning algorithms for the identification of regions and tracts of interest as potential biomarkers. Comput. Biol. Med. 185:109518. doi: 10.1016/j.compbiomed.2024.109518
Wang, Y., Sun, P., Wang, Q., Trinkaus, K., Schmidt, R., Naismith, R., et al. (2015). Differentiation and quantification of inflammation, demyelination and axon injury or loss in multiple sclerosis. Brain 138, 1223–1238. doi: 10.1093/brain/awv046
Whitwell, J., Avula, R., Master, A., Vemuri, P., Senjem, M., Jones, D., et al. (2011). Disrupted thalamocortical connectivity in PSP: A resting-state fMRI, DTI, and VBM study. Parkinsonism Relat. Disord. 17, 599–605. doi: 10.1016/j.parkreldis.2011.05.013
Whitwell, J., Höglinger, G., Antonini, A., Bordelon, Y., Boxer, A., Colosimo, C., et al. (2017). Radiological biomarkers for diagnosis in PSP: Where are we and where do we need to be? Mov. Disord. 32, 955–971. doi: 10.1002/mds.27038
Whitwell, J., Tosakulwong, N., Botha, H., Ali, F., Clark, H., Duffy, J., et al. (2020). Brain volume and flortaucipir analysis of progressive supranuclear palsy clinical variants. Neuroimage Clin. 25:102152. doi: 10.1016/j.nicl.2019.102152
Whitwell, J., Tosakulwong, N., Clark, H., Ali, F., Botha, H., Weigand, S., et al. (2021). Diffusion tensor imaging analysis in three progressive supranuclear palsy variants. J. Neurol. 268, 3409–3420. doi: 10.1007/s00415-020-10360-1
Whitwell, J., Tosakulwong, N., Schwarz, C., Botha, H., Senjem, M., Spychalla, A., et al. (2019). MRI outperforms [18F]AV-1451 PET as a longitudinal biomarker in progressive supranuclear palsy. Mov. Disord. 34, 105–113. doi: 10.1002/mds.27546
Winklewski, P., Sabisz, A., Naumczyk, P., Jodzio, K., Szurowska, E., and Szarmach, A. (2018). Understanding the physiopathology behind axial and radial diffusivity changes-what do we know? Front. Neurol. 9:92. doi: 10.3389/fneur.2018.00092
Worker, A., Blain, C., Jarosz, J., Chaudhuri, K., Barker, G., Williams, S., et al. (2014). Diffusion tensor imaging of Parkinson’s disease, multiple system atrophy and progressive supranuclear palsy: A tract-based spatial statistics study. PLoS One 9:e112638. doi: 10.1371/journal.pone.0112638
Zanigni, S., Evangelisti, S., Testa, C., Manners, D., Calandra-Buonaura, G., Guarino, M., et al. (2017). White matter and cortical changes in atypical parkinsonisms: A multimodal quantitative MR study. Parkinsonism Relat. Disord. 39, 44–51. doi: 10.1016/j.parkreldis.2017.03.001
Appendix
Members of the DESCRIBE-PSP Study Group
Moritz Brandt (Department of Neurology, University Hospital Carl Gustav Carus, Technische Universität Dresden, Dresden, Germany), Katharina Buerger (Institute for Stroke and Dementia Research, University Hospital, LMU Munich, Munich, Germany), Emrah Düzel (Institute of Cognitive Neurology and Dementia Research, Otto-von-Guericke University, Magdeburg, Germany), Björn Falkenburger (Department of Neurology, University Hospital Carl Gustav Carus, Technische Universität Dresden, Dresden, Germany), Agnes Flöel (Department of Neurology, University Medicine Greifswald, Greifswald, Germany), Wenzel Glanz (Clinic for Neurology, Medical Faculty, University Hospital Magdeburg, Magdeburg, Germany), Günter U. Höglinger (Department of Neurology, University Hospital, LMU Munich, Munich, Germany), Daniel Janowitz (Institute for Stroke and Dementia Research, University Hospital, LMU Munich, Munich, Germany), Sabrina Katzdobler (Department of Neurology, University Hospital of Munich, Ludwig-Maximilians-Universität (LMU) Munich, Munich, Germany), Ingo Kilimann (Department of Psychosomatic Medicine, Rostock University Medical Center, Rostock, Germany), Okka Kimmich (Department of Neurology, University of Bonn, Bonn, Germany), Johannes Levin (Department of Neurology, University Hospital of Munich, Ludwig-Maximilians-Universität (LMU) Munich, Munich, Germany), Oliver Peters (Charité – Universitätsmedizin Berlin, corporate member of Freie Universität Berlin and Humboldt-Universität zu Berlin-Institute of Psychiatry and Psychotherapy), Josef Priller (Department of Psychiatry and Psychotherapy, Charité, Berlin, Germany), Johannes Prudlo (Department of Neurology, University Medical Centre, Rostock, Germany), Luisa-Sophie Schneider (Charité – Universitätsmedizin Berlin, corporate member of Freie Universität Berlin and Humboldt-Universität zu Berlin-Institute of Psychiatry and Psychotherapy), Annika Spottke (Department of Neurology, University of Bonn, Bonn, Germany), Eike Jakob Spruth (Department of Psychiatry and Psychotherapy, Charité, Berlin, Germany), Matthis Synofzik (Department of Neurodegenerative Diseases, Hertie Institute for Clinical Brain Research and Center of Neurology, University of Tübingen, Tübingen, Germany), Stefan Teipel (Department of Psychosomatic Medicine, Rostock University Medical Center, Rostock, Germany), Carlo Wilke (Department of Neurodegenerative Diseases, Hertie Institute for Clinical Brain Research and Center of Neurology, University of Tübingen, Tübingen, Germany).
Keywords: progressive supranuclear palsy, sequential pattern, diffusion tensor imaging (DTI), neuropathology, tau protein
Citation: Bârlescu LA, Höglinger GU, Volkmann H, Ludolph AC, Del Tredici K, Braak H, Müller H-P and Kassubek J (2025) Diffusion tensor imaging of sequential neuropathological patterns in progressive supranuclear palsy. Front. Aging Neurosci. 17:1569302. doi: 10.3389/fnagi.2025.1569302
Received: 31 January 2025; Accepted: 07 May 2025;
Published: 06 June 2025.
Edited by:
Rodolfo Gabriel Gatto, Mayo Clinic, United StatesReviewed by:
Lydia Chougar, Université Paris-Sorbonne, FranceHongli Fan, Siemens Healthineers, United States
Francesca Orlandi, San Raffaele Hospital (IRCCS), Italy
Copyright © 2025 Bârlescu, Höglinger, Volkmann, Ludolph, Del Tredici, Braak, Müller and Kassubek. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jan Kassubek, amFuLmthc3N1YmVrQHVuaS11bG0uZGU=
†These authors have contributed equally to this work
‡ORCID: Günter U. Höglinger, orcid.org/0000-0001-7587-6187; Heiko Volkmann, orcid.org/0009-0004-5885-5498; Kelly Del Tredici, orcid.org/0000-0002-9130-7838; Heiko Braak, orcid.org/0000-0002-6776-1045; Hans-Peter Müller, orcid.org/0000-0001-5706-9381; Jan Kassubek, orcid.org/0000-0002-7106-9270