Dynamic coordination among three large-scale functional brain networks – the default mode network (DMN),
the central executive network (CEN), and the salience network (SN) – has been found to be aberrant in many
neuropsychiatric disorders, including obsessive-compulsive disorder (OCD). However, this line of research
typically relies on a group-based definition of the networks of interest: a standardized anatomical or functional
atlas or parcellation or a group-based independent component analysis (g-ICA). These group approaches
do not allow for a full accounting of individual variation in the structure of and relationships between these
large-scale networks. Discounting inter-individual heterogeneity, especially spatial variation (which has
been shown to be greater in several clinical populations, though not previously in OCD), may lead to failure
to detect significant effects not because they are not present, but because we did not target the right nodes
for every individual. We aim to address this problem in a systematic way, using existing data from patients
with OCD and in matched healthy controls.
Two approaches to remedy this problem have been proposed. An individualized probabilistic ICA
approach defines brain networks individually for each subject and then enters the individualized measures
into a second-level random effects analysis. A hierarchical probabilistic group ICA approach provides
model-based estimation of brain functional networks at both the population and subject level
simultaneously. Properties of the data (i.e., relations between subject-level and population-level variance
in variables of interest) determine whether hierarchical or non-hierarchical modeling will produce superior
results. We will analyze previously collected resting-state fMRI data from six studies (total 253 individuals
diagnosed with OCD, 148 of them were not on any medication, and 271 healthy controls without any
psychopathology, HC). We will derive subject-specific network maps and time courses using both individual
and hierarchical methods and use them to examine how subjects’ diagnosis, continuous clinical measures,
and non-clinical characteristics relate to (1) individual variability in topographical and functional organization
of DMN, SN, and CEN and (2) individual variability in functional coupling among these networks.
Prior research has demonstrated that rs-fMRI metrics of these brain-wide networks vary across OCD
dimensions and can serve as predictors of treatment response. However, calculation of these metrics does not
generally account for cross-subject topological heterogeneity, which can be misinterpreted as variations in
coupling. Our research will produce new and more precise markers of cross-patient heterogeneity, and de-
bias promising existing biomarkers of treatment response.