Brain-wide ongoing activity is responsible for significant cross-trial BOLD variability

A notorious issue of task-based functional magnetic resonance imaging (fMRI) is its large cross-trial variability. To quantitatively characterize this variability, the blood oxygenation level-dependent (BOLD) signal can be modeled as a linear summation of a stimulation-relevant and an ongoing (i.e. stimulation-irrelevant) component. However, systematic investigation on the spatiotemporal features of the ongoing BOLD component and how these features affect the BOLD response is still lacking. Here we measured fMRI responses to light onsets and light offsets in awake rats. The neuronal response was simultaneously recorded with calcium-based fiber photometry. We established that between-region BOLD signals were highly correlated brain-wide at zero time lag, including regions that did not respond to visual stimulation, suggesting that the ongoing activity co-fluctuates across the brain. Removing this ongoing activity reduced cross-trial variability of the BOLD response by ~30% and increased its coherence with the Ca2+ signal. Additionally, the negative ongoing BOLD activity sometimes dominated over the stimulation-driven response and contributed to the post-stimulation BOLD undershoot. These results suggest that brain-wide ongoing activity is responsible for significant cross-trial BOLD variability, and this component can be reliably quantified and removed to improve the reliability of fMRI response. Importantly, this method can be generalized to virtually all fMRI experiments without changing stimulation paradigms.