Mapping human preictal and ictal haemodynamic networks using simultaneous intracranial EEG-fMRI

Accurately characterising the brain networks involved in seizure activity may have important implications for our understanding of epilepsy. Intracranial EEG-fMRI can be used to capture focal epileptic events in humans with exquisite electrophysiological sensitivity and allows for identification of brain structures involved in this phenomenon over the entire brain. We investigated ictal BOLD networks using the simultaneous intracranial EEG-fMRI (icEEG-fMRI) in a 30 year-old male undergoing invasive presurgical evaluation with bilateral depth electrode implantations in amygdalae and hippocampi for refractory temporal lobe epilepsy. One spontaneous focal electrographic seizure was recorded. The aims of the data analysis were firstly to map BOLD changes related to the ictal activity identified on icEEG and secondly to compare different fMRI modelling approaches. Visual inspection of the icEEG showed an onset dominated by beta activity involving the right amygdala and hippocampus lasting 6.4 s (ictal onset phase), followed by gamma activity bilaterally lasting 14.8 s (late ictal phase). The fMRI data was analysed using SPM8 using two modelling approaches: firstly, purely based on the visually identified phases of the seizure and secondly, based on EEG spectral dynamics quantification. For the visual approach the two ictal phases were modelled as 'ON' blocks convolved with the haemodynamic response function; in addition the BOLD changes during the 30 s preceding the onset were modelled using a flexible basis set. For the quantitative fMRI modelling approach two models were evaluated: one consisting of the variations in beta and gamma bands power, thereby adding a quantitative element to the visually-derived models, and another based on principal components analysis of the entire spectrogram in attempt to reduce the bias associated with the visual appreciation of the icEEG. BOLD changes related to the visually defined ictal onset phase were revealed in the medial and lateral right temporal lobe. For the late ictal phase, the BOLD changes were remote from the SOZ and in deep brain areas (precuneus, posterior cingulate and others). The two quantitative models revealed BOLD changes involving the right hippocampus, amygdala and fusiform gyrus and in remote deep brain structures and the default mode network-related areas. In conclusion, icEEG-fMRI allowed us to reveal BOLD changes within and beyond the SOZ linked to very localised ictal fluctuations in beta and gamma activity measured in the amygdala and hippocampus. Furthermore, the BOLD changes within the SOZ structures were better captured by the quantitative models, highlighting the interest in considering seizure-related EEG fluctuations across the entire spectrum.