Sensitivity of the resting-state haemodynamic response function estimation to autonomic nervous system fluctuations

The haemodynamic response function (HRF) is a key component of the blood oxygen level-dependent (BOLD) signal, providing the mapping between neural activity and the signal measured with functional magnetic resonance imaging (fMRI). Most of the time the HRF is associated with task-based fMRI protocols, in which its onset is explicitly included in the design matrix. On the other hand, the HRF also mediates the relationship between spontaneous neural activity and the BOLD signal in resting-state protocols, in which no explicit stimulus is taken into account. It has been shown that resting-state brain dynamics can be characterized by looking at sparse BOLD 'events', which can be retrieved by point process analysis. These events can be then used to retrieve the HRF at rest. Crucially, cardiac activity can also induce changes in the BOLD signal, thus affecting both the number of these events and the estimation of the haemodynamic response. In this study, we compare the resting-state haemodynamic response retrieved by means of a point process analysis, taking the cardiac fluctuations into account. We find that the resting-state HRF estimation is significantly modulated in the brainstem and surrounding cortical areas. From the analysis of two high-quality datasets with different temporal and spatial resolution, and through the investigation of intersubject correlation, we suggest that spontaneous point process response durations are associated with the mean interbeat interval and low-frequency power of heart rate variability in the brainstem.