Mapping the neuroanatomical heterogeneity of OCD using a framework integrating normative model and non-negative matrix factorization

Obsessive–compulsive disorder (OCD) is a spectrum disorder with high interindividual heterogeneity. We propose a comprehensible framework integrating normative model and non-negative matrix factorization (NMF) to quantitatively estimate the neuroanatomical heterogeneity of OCD from a dimensional perspective. T1-weighted magnetic resonance images of 98 first-episode untreated patients with OCD and matched healthy controls (HCs, n = 130) were acquired. We derived individualized differences in gray matter morphometry using normative model and parsed them into latent disease factors using NMF. Four robust disease factors were identified. Each patient expressed multiple factors and exhibited a unique factor composition. Factor compositions of patients were significantly correlated with severity of symptom, age of onset, illness duration, and exhibited sex differences, capturing sources of clinical heterogeneity. In addition, the group-level morphological differences obtained with two-sample t test could be quantitatively derived from the identified disease factors, reconciling the group-level and subject-level findings in neuroimaging studies. Finally, we uncovered two distinct subtypes with opposite morphological differences compared with HCs from factor compositions. Our findings suggest that morphological differences of individuals with OCD are the unique combination of distinct neuroanatomical patterns. The proposed framework quantitatively estimating neuroanatomical heterogeneity paves the way for precision medicine in OCD.