Abstract:Large-scale dynamics of the brain are routinely modelled using systems of nonlinear dynamical equations that describe the evolution of population-level activity, with distinct neural populations often coupled according to an empirically measured structural connection matrix. This modelling approach has been used to generate insights into the neural underpinnings of spontaneous brain dynamics, as recorded with techniques such as resting state functional MRI (fMRI). In fMRI, researchers have many degrees of freedom in the way that they can process the data and recent evidence indicates that the choice of pre-processing steps can have a major effect on empirical estimates of functional connectivity. However, the potential influence of such variations on modelling results are seldom considered. Here we show, using three popular whole-brain dynamical models, that different choices during fMRI preprocessing can dramatically affect model fits and interpretations of findings. Critically, we show that the ability of these models to accurately capture patterns in fMRI dynamics is mostly driven by the degree to which they fit global signals rather than interesting sources of coordinated neural dynamics. We show that widespread deflections can arise from simple global synchronisation. We introduce a simple two-parameter model that captures these fluctuations and which performs just as well as more complex, multi-parameter biophysical models. From our combined analyses of data and simulations, we describe benchmarks to evaluate model fit and validity. Although most models are not resilient to denoising, we show that relaxing the approximation of homogeneous neural populations by more explicitly modelling inter-regional effective connectivity can improve model accuracy at the expense of increased model complexity. Our results suggest that many complex biophysical models may be fitting relatively trivial properties of the data, and underscore a need for tighter integration between data quality assurance and model development.

Stochastic brain dynamics exhibits differential regional distribution and maturation-related changes

Characterization of regional differences in resting-state fMRI with a data-driven network model of brain dynamics

Quantitative Evaluation of Activation State in Functional Brain Imaging

Integrating group-wise functional brain activities via point processes

Dynamic nonreversibility view of intrinsic brain organization and brain dynamic analysis of repetitive transcranial magnitude stimulation

Criticality in Large-Scale Brain fMRI Dynamics Unveiled by a Novel Point Process Analysis

On the intersection between data quality and dynamical modelling of large-scale fMRI signals

Brain dynamics supported by a hierarchy of complex correlation patterns defining a robust functional architecture

Exploring Physiological Differences in Brain Areas Using Statistical Complexity Analysis of BOLD Signals

Stochastic neurodynamical computation of brain functions

State-transition dynamics of resting-state functional magnetic resonance imaging data: Model comparison and test-to-retest analysis

Identification of brain states, transitions, and communities using functional MRI

Metastability, fractal scaling, and synergistic information processing: What phase relationships reveal about intrinsic brain activity

Evolutionarily conserved fMRI network dynamics in the mouse, macaque, and human brain

Parameter Estimation in Brain Dynamics Models from Resting-State fMRI Data using Physics-Informed Neural Networks

On the Stability of BOLD fMRI Correlations

Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging

Estimating effective connectivity in linear brain network models

The physics of functional magnetic resonance imaging (fMRI)