Pathway-Anchored Multimodal Clustering for Parkinson's Disease

Pathway-Anchored Multimodal Clustering for Parkinson's Disease

Ashwin Vinod, Aditya Sai Ellendula, Shubham Bhardwaj, Aparna Dev, Aaron Dominic, Chandrajit Bajaj

Department of Computer Science, The University of Texas at Austin

Department of Computer Science, Amrita Vishwa Vidyapeetham, Coimbatore, India

bioRxiv preprint   |   Paper


Annotated pathway-anchored Parkinson's disease circuit map

Figure 1. Annotated pathway-anchored circuit map showing six PD-relevant pathway families, their major anatomical regions, and their clinical interpretation. The figure emphasizes that multimodal imaging features are interpreted through predefined disease-relevant circuits rather than through an unstructured connectome.

Pathway-anchored multimodal framework for circuit-level Parkinson's disease imaging analysis

Figure 2. Pathway-anchored multimodal framework for circuit-level Parkinson's disease imaging analysis. (A) Six PD-relevant pathway families are mapped onto color-coded brain regions and linked through known cross-system interactions. Solid lines indicate anatomical localization, double-headed arrows indicate pathway interactions, and dashed arrows indicate the modifying influence of microvascular burden. Insets highlight deep nuclei, white-matter hyperintensity burden, and cerebellar structures. (B) T1 MRI, diffusion MRI-derived microstructural measures, and DaT-SPECT provide ROI-level structural, microstructural, and dopaminergic features. These features are assembled into a subject-by-feature matrix, mapped to predefined pathway bins, feature-wise z-scored, and summarized as Multimodal Pathway Integrity Scores (MPIS). (C) SRVCC co-clustering identifies imaging-derived patient groups and feature modules. Pathway-level MPIS profiles are then associated with motor, cognitive, and behavioral measures, supporting circuit-level interpretation of patient heterogeneity.


TL;DR

This work organizes multimodal Parkinson's imaging by neurobiological pathway rather than treating the brain as an undifferentiated feature bank. Using structural MRI, diffusion MRI, and DaT-SPECT from the PPMI cohort, the paper combines SRVCC co-clustering with a Multimodal Pathway Integrity Score (MPIS) to recover circuit-level signatures that track motor and cognitive burden while remaining interpretable at the pathway level.


Why this matters

Parkinson's disease is heterogeneous, but many multimodal analyses still produce diffuse latent factors that are hard to interpret biologically. This paper instead asks whether patient structure becomes clearer when imaging features are grouped into known circuits such as nigrostriatal, frontostriatal, limbic, sensory-visuospatial, cerebello-thalamo-cortical, and microvascular pathways.

The resulting representation is intended to support statements like "lower nigrostriatal integrity aligns with greater motor burden" rather than opaque cluster labels with no mechanistic interpretation.


Method

The dataset includes 294 participants from PPMI: 185 with Parkinson's disease, 72 healthy controls, and 37 SWEDD participants. The primary multimodal view combines structural MRI volumetry, free-water-corrected diffusion MRI measures, and DaT-SPECT specific binding ratios.

For each predefined pathway, ROI-level features are standardized and aggregated into a Multimodal Pathway Integrity Score (MPIS). Higher MPIS reflects a more intact pathway signature through higher FA, higher volume and SBR, and lower mean diffusivity.

The subject-by-feature matrix is then analyzed with Scalable Robust Variational Compositional Co-clustering (SRVCC), which jointly learns patient clusters and feature clusters while preserving a pathway-aware interpretation.

Multimodal segmentation and quantification pipeline

Figure 3. Structural MRI defines the reference space, diffusion MRI and DaT-SPECT are aligned into it, and regional pathway features are extracted for downstream co-clustering.

Cerebello-thalamo-cortical, limbic, microvascular, and sensory pathway schematics

Figure 4. Pathway schematics for (A) cerebello-thalamo-cortical, (B) limbic/mesolimbic, (C) microvascular-burden, and (D) sensory/visual-auditory/visuospatial-attention pathways.

Frontostriatal and nigrostriatal pathway schematics

Figure 5. Pathway schematics of (A) frontostriatal (cognitive/executive/attention) and (B) nigrostriatal (motor basal ganglia-thalamo-cortical).


Key findings

  • Model selection favored a compact SRVCC configuration of (Kr,Kc)=(5,5)(K_r, K_c) = (5, 5), balancing reconstruction quality and cluster structure.
  • Cluster stability was strong across both random seeds and bootstrap resampling: median ARI/NMI of 0.82/0.88 across seeds and 0.76/0.84 across bootstrap runs.
  • Nigrostriatal MPIS showed the clearest motor association, with lower integrity linked to higher MDS-UPDRS III burden (ρ0.201\rho \approx -0.201, q6.8×104q \approx 6.8 \times 10^{-4}).
  • Frontostriatal MPIS also tracked motor burden (ρ0.191\rho \approx -0.191, q0.0014q \approx 0.0014), supporting the role of executive-attention circuitry in Parkinsonian heterogeneity.
  • Sensory/visuospatial MPIS had the strongest cognition association, with higher integrity linked to better MoCA performance (ρ0.163\rho \approx 0.163, q0.0071q \approx 0.0071).
  • Limbic/mesolimbic MPIS showed a smaller but still positive association with cognition (ρ0.119\rho \approx 0.119, q0.045q \approx 0.045).

Nigrostriatal MPIS versus motor burden

Figure 6. Nigrostriatal pathway results: lower MPIS aligns with worse motor burden and strong separation across imaging-defined clusters.

Sensory-visuospatial MPIS versus cognition

Figure 7. Sensory/visuospatial pathway results: higher MPIS aligns with better global cognition and strong cluster separation.


Robustness and Interpretation

The paper explicitly frames these effects as coherent but modest pathway-level associations rather than a turnkey diagnostic system. That matters: the value here is not a single headline classifier, but an interpretable multimodal summary that links dopaminergic, microstructural, and morphometric changes back to specific circuits.

The reported pathway patterns were also stable to common MPIS design choices. Across normalization, modality reweighting, and MD sign variants, the primary MPIS remained highly concordant with alternatives, with median Pearson correlation around 0.96 and directionally consistent clinical associations.


Citation

@article{vinod2025pathway,
  title   = {Pathway Anchored Multimodal Clustering Reveals Circuit Level Signatures in Parkinsons Disease},
  author  = {Vinod, Ashwin and Ellendula, Aditya Sai and Bhardwaj, Shubham and Dev, Aparna and Dominic, Aaron and Bajaj, Chandrajit},
  journal = {bioRxiv},
  year    = {2025},
  doi     = {10.64898/2025.12.15.694278}
}