Increased sensitivity of fast BOLD fMRI with a subject-specific hemodynamic response function and application to epilepsy

• The HRF is more variable across subjects than brain regions.
• Subject to subject variation of HRF waveform correlates between regions.
• A subject-specific HRF model leads to more sensitive detection of task activations.
• The approach yields complementary information in epileptic patient EEG-fMRI exams.

Activation detection in functional Magnetic Resonance Imaging (fMRI) typically assumes the hemodynamic response to neuronal activity to be invariant across brain regions and subjects. Reports of substantial variability of the morphology of blood-oxygenation-level-dependent (BOLD) responses are accumulating, suggesting that the use of a single generic model of the expected response in general linear model (GLM) analyses does not provide optimal sensitivity due to model misspecification. Relaxing assumptions of the model can limit the impact of hemodynamic response function (HRF) variability, but at a cost on model parsimony. Alternatively, better specification of the model could be obtained from a priori knowledge of the HRF of a given subject, but the effectiveness of this approach has only been tested on simulation data. Using fast BOLD fMRI, we characterized the variability of hemodynamic responses to a simple event-related auditory-motor task, as well as its effect on activation detection with GLM analyses. We show the variability to be higher between subjects than between regions and variation in different regions to correlate from one subject to the other. Accounting for subject-related variability by deriving subject-specific models from responses to the task in some regions lead to more sensitive detection of responses in other regions. We applied the approach to epilepsy patients, where task-derived patient-specific models provided additional information compared to the use of a generic model for the detection of BOLD responses to epileptiform activity identified on scalp electro-encephalogram (EEG). This work highlights the importance of improving the accuracy of the model for detecting neuronal activation with fMRI, and the fact that it can be done at no cost to model parsimony through the acquisition of independent a priori information about the hemodynamic response.