Study-level wavelet cluster analysis and data-driven signal models in pharmacological MRI

In pharmacological MRI (phMRI) studies tracking signal changes following the acute administration of a compound, the spatiotemporal pattern of response is often unknown a priori. Moreover, when analysed within a general linear model (GLM) framework, the experimental paradigm of a single injection point under-informs the construction of an appropriate signal model, and information from pharmacokinetics or ancillary in vivo studies may be unavailable or insufficient to accurately describe the dynamic signal changes observed following injection of the drug. Here, we extend the application of a data-driven clustering algorithm, wavelet cluster analysis (WCA), to phMRI data from one or more groups of subjects in a study. A WCA decomposition of spatially concatenated time series’ provides a compact overview of spatiotemporal response patterns across cohorts, highlighting typical temporal signatures, brain regions implicated in the response and inter-subject variability. Further, we demonstrate the use of regressors based on selected temporal components as suitable signal models in GLM-based analyses, resulting in a close fit to dynamic phMRI signal changes. This approach is illustrated with simulated data and two representative in vivo phMRI studies in the rat (nicotine and apomorphine challenges).