Large-scale dynamic modeling of task-fMRI signals via subspace system identification.

OBJECTIVE: We analyze task-based fMRI time series to produce large-scale dynamical models that are capable of approximating the observed signal with good accuracy. APPROACH: We extend subspace system identification methods for deterministic and stochastic state-space models with external inputs. The dynamic behavior of the generated models is characterized using control-theoretic analysis tools. To validate their effectiveness, we perform a probabilistic inversion of the identified input-output relationships via joint state-input maximum likelihood estimation. Our experimental setup explores a large dataset generated using state-of-the-art acquisition and pre-processing methods from the Human Connectome Project. MAIN RESULTS: We analyze both anatomically parcellated and spatially dense time series, and propose an efficient algorithm to address the high-dimensional optimization problem resulting from the latter. Our results enable the quantification of input-output transfer functions between each task condition and each region of the cortex, as exemplified by a motor task. Further, the identified models produce impulse response functions between task conditions and cortical regions that are compatible with typical hemodynamic response functions. We then extend subspace methods to account for multi-subject experimental configurations, identifying models that capture common dynamical characteristics across subjects. Finally, we show that system inversion via maximum-likelihood allows the time-of-occurrence of the task stimuli to be estimated from the observed outputs. SIGNIFICANCE: The ability to produce dynamical input-output models might have an impact in the expanding field of neurofeedback. In particular, the models we produce allow the partial quantification of the effect of external task-related inputs on the metabolic response of the brain, conditioned on its current state. Such a notion provides a basis for leveraging control-theoretic approaches to neuromodulation and self-regulation in therapeutic applications.

Spectral mapping of brain functional connectivity from diffusion imaging.

Understanding the relationship between the dynamics of neural processes and the anatomical substrate of the brain is a central question in neuroscience. On the one hand, modern neuroimaging technologies, such as diffusion tensor imaging, can be used to construct structural graphs representing the architecture of white matter streamlines linking cortical and subcortical structures. On the other hand, temporal patterns of neural activity can be used to construct functional graphs representing temporal correlations between brain regions. Although some studies provide evidence that whole-brain functional connectivity is shaped by the underlying anatomy, the observed relationship between function and structure is weak, and the rules by which anatomy constrains brain dynamics remain elusive. In this article, we introduce a methodology to map the functional connectivity of a subject at rest from his or her structural graph. Using our methodology, we are able to systematically account for the role of structural walks in the formation of functional correlations. Furthermore, in our empirical evaluations, we observe that the eigenmodes of the mapped functional connectivity are associated with activity patterns associated with different cognitive systems.