The Hemodynamic Response Function Varies Across Anatomical Location and Pathology in the Epileptic Brain

The hemodynamic response function (HRF) links neuronal activity to functional magnetic resonance imaging (fMRI) signals. While most fMRI studies use a "canonical" HRF, increasing evidence from studies of healthy subjects suggests that the HRF depends on anatomical location and disease states. Here, we investigate how HRF variability relates to anatomical location and pathology in the epileptic brain, using a large simultaneous electroencephalogram and fMRI dataset. Applying HRF deconvolution and temporal decomposition, we built the first whole-brain HRF library specific to epilepsy, identifying four distinct shape groups. We mapped HRF features across parcellations of two atlases using novel Bayesian hierarchical models. In non-epileptogenic regions, HRF shape and spatial distributions align with findings from healthy subjects. Within pathological regions, they vary significantly according to pathology. Our results indicate that HRF variability is associated with pathology, in addition to its dependence on anatomical location, motivating region- and pathology-based HRF modulation in epilepsy studies.