Comparison of hemodynamic response functions estimated from resting-state functional MRI and invasive electrophysiological recordings in rats

Resting-state functional MRI (rs-fMRI) is a popular technology that has enriched our understanding of brain and spinal cord functioning, including how different regions communicate (connectivity). But fMRI is an indirect measure of neural activity capturing blood hemodynamics. The hemodynamic response function (HRF) interfaces between the unmeasured neural activity and measured fMRI time series. The HRF is variable across brain regions and individuals, and is modulated by non-neural factors. Ignoring this HRF variability causes errors in FC estimates. Hence, it is crucial to reliably estimate the HRF from rs-fMRI data. Robust techniques have emerged to estimate the HRF from fMRI time series. Although such techniques have been validated non-invasively using simulated and empirical fMRI data, thorough invasive validation using simultaneous electrophysiological recordings, the gold standard, has been elusive. This report addresses this gap in the literature by comparing HRFs derived from invasive intracranial electroencephalogram recordings with HRFs estimated from simultaneously acquired fMRI data in six epileptic rats. We found that the HRF shape parameters (HRF amplitude, latency and width) were not significantly different (p>0.05) between ground truth and estimated HRFs. In the single pathological region, the HRF width was marginally significantly different (p=0.03). Our study provides preliminary invasive validation for the efficacy of the HRF estimation technique in reliably estimating the HRF non-invasively from rs-fMRI data directly. This has a notable impact on rs-fMRI connectivity studies, and we recommend that HRF deconvolution be performed to minimize HRF variability and improve connectivity estimates.