Spatial collinearity constrains multivariate molecular-enriched network estimation.

Spatial overlap among brain receptor maps derived from PET imaging can distort analyses that model multiple receptors together. Using test-retest fMRI data, the authors show that as more receptors are included in a multivariate model, the reliability of the resulting functional connectivity networks decreases, and this degradation is driven by collinearity among the receptor maps. A univariate approach, modeling each receptor independently, produces more reliable networks and, in a study comparing LSD to placebo, better captured the known role of the 5HT-2A receptor. Spatial collinearity is a fundamental constraint on multivariate molecular-enriched network estimation, and univariate modeling is recommended as a more robust default.

Study at a glance

Abstract

Analyses of neuroimaging data increasingly leverage the distribution of neurotransmitter receptors derived from Positron Emission Tomography (PET) to bridge the gap between micro- and macro-scale brain function. However, these receptor maps are highly spatially overlapping which can give rise to interpretive and analytical challenges. Here, we systematically investigate the impact of spatial collinearity among PET maps in the context of Receptor-Enriched Analysis of functional Connectivity by Targets (REACT), a method that uses receptor maps as spatial regressors to derive subject-level molecular-enriched functional connectivity networks. Exhaustive combinatorial analysis across 19 receptor and transporter maps showed that collinearity scales rapidly with the number of receptors modelled simultaneously, and that this was relatively stable across parcellation scales, reflecting the intrinsic organisation of neurotransmitter systems. Using test-retest fMRI data from the Human Connectome Project, we demonstrate that modelling greater numbers of receptors degrades the reliability of molecular-enriched networks derived from conventional multivariate REACT models, and that collinearity among receptor maps drives this degradation. An alternative univariate approach, in which each receptor is modelled independently, yielded more reliable networks and, when applied to a within-subjects study of LSD compared to placebo, better recovered the role of the 5HT-2A receptor in LSD's neural effects. These findings identify spatial collinearity as a fundamental constraint on multivariate molecular-enriched network estimation and support univariate modelling as a more robust default for this class of analysis.