Article thumbnail
Location of Repository

Imaging haemodynamic changes related to seizures: comparison of EEG-based general linear model, independent component analysis of fMRI and intracranial EEG

By R.C. Thornton, R. Rodionov, H. Laufs, S. Vulliemoz, A. Vaudano, D. Carmichael, S. Cannadathu, M. Guye, A. McEvoy and S. Lhatoo


Background: Simultaneous EEG-fMRI can reveal haemodynamic changes associated with epileptic activity which may contribute to understanding seizure onset and propagation. Methods: Nine of 83 patients with focal epilepsy undergoing pre-surgical evaluation had seizures during EEG-fMRI and analysed using three approaches, two based on the general linear model (GLM) and one using independent component analysis (ICA): 1. EEGs were divided into up to three phases: early ictal EEG change, clinical seizure onset and late ictal EEG change and convolved with a canonical haemodynamic response function (HRF) (canonical GLM analysis). 2. Seizures lasting three scans or longer were additionally modelled using a Fourier basis set across the entire event (Fourier GLM analysis). 3. Independent component analysis (ICA) was applied to the fMRI data to identify ictal BOLD patterns without EEG. The results were compared with intracranial EEG. Results: The canonical GLM analysis revealed significant BOLD signal changes associated with seizures on EEG in 7/9 patients, concordant with the seizure onset zone in 4/7. The Fourier GLM analysis revealed changes in BOLD signal corresponding with the results of the canonical analysis in two patients. ICA revealed components spatially concordant with the seizure onset zone in all patients (8/9 confirmed by intracranial EEG). Conclusion: Ictal EEG-fMRI visualises plausible seizure related haemodynamic changes. The GLM approach to analysing EEG-fMRI data reveals localised BOLD changes concordant with the ictal onset zone when scalp EEG reflects seizure onset. ICA provides additional information when scalp EEG does not accurately reflect seizures and may give insight into ictal haemodynamics

Year: 2010
OAI identifier:
Provided by: UCL Discovery

Suggested articles


  1. (1999). 2) A single block may not accurately represent dynamic processes such as seizures (Niedermeyer et al
  2. (1998). a Method for Its Subtraction.
  3. (2000). A Method for Removing Imaging Artifact from Continuous EEG Recorded during Functional MRI.
  4. (2001). A Proposed Diagnostic Scheme for People with Epileptic Seizures and with Epilepsy:
  5. (1999). Abnormal EEG Patterns: Epileptic and Paroxysmal.
  6. Analysis (ICA) Spatial ICA (Formisano et al., 2004; McKeown et al., 2003) as implemented in Voyager QX software (Brain Innovation,
  7. (2007). Analysis ICA is limited by the large number of components generated, resulting in low specificity, which was addressed by the use of an automatic component classification scheme (De Martino et al.,
  8. (2001). Anticipation of epileptic seizures from standard EEG recordings.
  9. (2009). Causal hierarchy within the thalamo-cortical network in spike and wave discharges.
  10. (1995). Characterizing Dynamic Brain Responses with fMRI: A Multivariate Approach.
  11. (2007). Classification of fMRI independent components using IC-fingerprints and support vector machine classifiers.
  12. (2004). Cortex-based independent component analysis of fMRI time series.
  13. (2001). Current trends in electroencephalography.
  14. (2008). Different structures involved during ictal and interictal epileptic activity in malformations of cortical development: an EEG-fMRI study.
  15. (2006). EEG-fMRI of epileptic spikes: Concordance with EEG source localization and intracranial EEG.
  16. (2006). EEG-fMRI of idiopathic and secondarily generalized epilepsies.
  17. (2000). EEG-triggered functional MRI in patients with pharmacoresistant epilepsy.
  18. (2003). Electroencephalographic signatures of attentional and cognitive default modes in spontaneous brain activity fluctuations at rest.
  19. (2007). Electrophysiological signatures of resting state networks in the human brain.
  20. (1996). Extraction of spatio-temporal signatures from depth EEG seizure signals based on objective matching in warped vectorial observations.
  21. (2003). fMRI "deactivation" of the posterior cingulate during generalized spike and wave.
  22. (2007). Focal increases in perfusion and decreases in hemoglobin oxygenation precede seizure onset in spontaneous human epilepsy.
  23. (2005). Functional MRI of the pre-ictal state.
  24. (1999). Given that ictal EEG does not always reflect seizure onset
  25. (2006). Hemodynamic correlates of epileptiform discharges: an EEG-fMRI study of 63 patients with focal epilepsy. Brain Res.
  26. (2009). Identifying the structures involved in seizure generation using sequential analysis of ictal-fMRI data.
  27. (2006). In most EEG-fMRI studies of ictal activity, seizures have been represented as a single block across their duration convolved with an HRF
  28. (2001). Inaugural Article: A default mode of brain function.
  29. (2008). Independent component analysis as a model-free approach for the detection of BOLD changes related to epileptic spikes: A simulation study.
  30. (2003). Independent component analysis of functional MRI: what is signal and what is noise?
  31. (2007). Independent component analysis of interictal fMRI in focal epilepsy: Comparison with general linear modelbased EEG-correlated fMRI.
  32. (2010). Independent component analysis reveals dynamic ictal BOLD responses in EEG-fMRI data from focal epilepsy patients.
  33. (1995). Localization of subclinical ictal activity by functional magnetic resonance imaging: correlation with invasive monitoring.
  34. (2006). Modelling cardiac signal as a confound in EEG-fMRI and its application in focal epilepsy studies.
  35. (1999). Modern electroencephalography: its role in epilepsy management.
  36. (1996). Motion-related effects were included in the GLM as 24 regressors (6 scan realignment parameters and a Volterra expansion of these (Friston
  37. (1991). Network modulation of a slow intrinsic oscillation of cat thalamocortical neurons implicated in sleep delta waves: cortically induced synchronization and brainstem cholinergic suppression.
  38. (2006). Neurovascular coupling and oximetry during epileptic events.
  39. (1998). Non-invasive epileptic focus localization using EEG-triggered functional MRI and electromagnetic tomography.
  40. (2008). Noncanonical spike-related BOLD responses in focal epilepsy.
  41. (2009). Noninvasive dynamic imaging of seizures in epileptic patients. Hum.Brain Mapp.
  42. (2010). Our findings of ICs spatially concordant with the seizure onset zone, complement a recent study in which spatially concordant ICs were identified in patients in whom a GLM analysis had previously revealed seizure-related changes (Levan et al.,
  43. Presurgical evaluation for partial epilepsy: Relative contributions of chronic depth-electrode recordings versus FDG-PET and scalp-sphenoidal ictal EEG.
  44. (2004). Probabilistic independent component analysis for functional magnetic resonance imaging.
  45. (2000). Quantification of, visualization of, and compensation for brain shift using intraoperative magnetic resonance imaging.
  46. (2002). Simultaneous EEG-Correlated Ictal fMRI.
  47. (2009). Spatiotemporal Dynamics of Perfusion and Oximetry during Ictal Discharges in the Rat Neocortex.
  48. (2005). Surgical outcome and prognostic factors of cryptogenic neocortical epilepsy.
  49. (2006). Temporal lobe interictal epileptic discharges affect cerebral activity in "default mode" brain regions. Hum.Brain Mapp.
  50. (2001). The localizing value of ictal EEG in focal epilepsy.
  51. (2004). The mean numbers of IC classified as BOLD and EPI per 20 min EEG-fMRI run were 7 (range 3-14) and 24 (range 14-36) respectively. Over both runs a mean of 6 ICs were characteristic of normal physiological activity
  52. (2009). The spatio-temporal mapping of epileptic networks: combination of EEG-fMRI and EEG source imaging.
  53. (1999). the technique constrains the analysis to grey matter and uses a fixed point ICA algorithm to separate signals generated by different sources (FastICA, Hyvärinen,
  54. (1990). The technique is unrestricted by the limited sampling and invasive nature of intracranial
  55. (2007). was designed to be inclusive with respect to BOLD components to reduce the probability of misclassification of a BOLD-related IC.
  56. (2007). We used ICA combined with a component classifier [De Martino
  57. (1995). while the canonical HRF appears to be a suitable model in normal physiological conditions (Friston
  58. (2006). Widespread and intense BOLD changes during brief focal electrographic seizures.

To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.