Independent component analysis reveals dynamic ictal BOLD responses in EEG-fMRI data from focal epilepsy patients

INTRODUCTION Seizures occur rarely during EEG-fMRI acquisitions of epilepsy patients, but can potentially offer a better estimation of the epileptogenic zone than interictal activity. Independent component analysis (ICA) is a data-driven method that imposes minimal constraints on the hemodynamic response function (HRF). In particular, the investigation of HRFs with clear peaks, but varying latency, may be used to differentiate the ictal focus from propagated activity. METHODS ICA was applied on ictal EEG-fMRI data from 15 patients. Components related to seizures were identified by fitting an HRF to the component time courses at the time of the ictal EEG events. HRFs with a clear peak were used to derive maps of significant BOLD responses and their associated peak delay. The results were then compared with those obtained from a general linear model (GLM) method. Concordance with the presumed epileptogenic focus was also assessed. RESULTS The ICA maps were significantly correlated with the GLM maps for each patient (Spearman's test, p<0.05). The ictal BOLD responses identified by ICA always included the presumed epileptogenic zone, but were also more widespread, accounting for 20.3% of the brain volume on average. The method provided a classification of the components as a function of peak delay. BOLD response clusters associated with early HRF peaks were concordant with the suspected epileptogenic focus, while subsequent HRF peaks may correspond to ictal propagation. CONCLUSION ICA applied to EEG-fMRI can detect areas of significant BOLD response to ictal events without having to predefine an HRF. By estimating the HRF peak time in each identified region, the method could also potentially provide a dynamic analysis of ictal BOLD responses, distinguishing onset from propagated activity.