Mapping Epileptic Networks Using Simultaneous Intracranial EEG-fMRI

Abstract

Background: Potentially curative epilepsy surgery can be offered if a single, discrete epileptogenic zone (EZ) can be identified. For individuals in whom there is no clear concordance between clinical localization, scalp EEG, and imaging data, intracranial EEG (icEEG) may be needed to confirm a predefined hypothesis regarding irritative zone (IZ), seizure onset zone (SOZ), and EZ prior to surgery. However, icEEG has limited spatial sampling and may fail to reveal the full extent of epileptogenic network if predefined hypothesis is not correct. Simultaneous icEEG-fMRI has been safely acquired in humans and allows exploration of neuronal activity at the whole-brain level related to interictal epileptiform discharges (IED) captured intracranially.

Methods: We report icEEG-fMRI in eight patients with refractory focal epilepsy who had resective surgery and good postsurgical outcome. Surgical resection volume in seizure-free patients post-surgically reflects confirmed identification of the EZ. IEDs on icEEG were classified according to their topographic distribution and localization (Focal, Regional, Widespread, and Non-contiguous). We also divided IEDs by their location within the surgical resection volume [primary IZ (IZ1) IED] or outside [secondary IZ (IZ2) IED]. The distribution of fMRI blood oxygen level-dependent (BOLD) changes associated with individual IED classes were assessed over the whole brain using a general linear model. The concordance of resulting BOLD map was evaluated by comparing localization of BOLD clusters with surgical resection volume. Additionally, we compared the concordance of BOLD maps and presence of BOLD clusters in remote brain areas: precuneus, cuneus, cingulate, medial frontal, and thalamus for different IED classes.

Results: A total of 38 different topographic IED classes were identified across the 8 patients: Focal (22) and non-focal (16, Regional = 9, Widespread = 2, Non-contiguous = 5). Twenty-nine IEDs originated from IZ1 and 9 from IZ2. All IED classes were associated with BOLD changes. BOLD maps were concordant with the surgical resection volume for 27/38 (71%) IED classes, showing statistical global maximum BOLD cluster or another cluster in the surgical resection volume. The concordance of BOLD maps with surgical resection volume was greater (p < 0.05) for non-focal (87.5%, 14/16) as compared to Focal (59%, 13/22) IED classes. Additionally, BOLD clusters in remote cortical and deep brain areas were present in 84% (32/38) of BOLD maps, more commonly (15/16; 93%) for non-focal IED-related BOLD maps.

Conclusions: Simultaneous icEEG-fMRI can reveal BOLD changes at the whole-brain level for a wide range of IEDs on icEEG. BOLD clusters within surgical resection volume and remote brain areas were more commonly seen for non-focal IED classes, suggesting that a wider hemodynamic network is at play.

Introduction

Intracranial electroencephalography (icEEG) recordings are performed during presurgical evaluation to localize irritative zone (IZ), seizure onset zone (SOZ), epileptogenic zone (EZ), and eloquent cortex for patients being considered for epilepsy surgery. icEEG has better spatial resolution and sensitivity compared to scalp EEG (1), which has low sensitivity (2, 3) and can provide inaccurate localization (4) and even lateralization (5), especially in patients with frontal lobe epilepsy. icEEG, however, has limited spatial sampling, only detecting electrical activity within a 1-cm core of tissue from recording site (6), and carries surgical risk (7, 8).

Simultaneous scalp EEG and functional magnetic resonance imaging (EEG-fMRI) can map interictal epileptiform discharges (IED) and seizure-related blood-oxygen-level-dependent (BOLD) changes (9–21). In patients with focal cortical dysplasia, IED-related BOLD changes distributed over multiple lobes are associated with poor postsurgical outcome (22). The interpretation of EEG-fMRI findings is often limited by the low sensitivity of scalp EEG, low rates of IED, and an incomplete understanding of the relationship between IED and BOLD changes (12, 13, 23–25). The sensitivity of EEG-fMRI can be increased by using topographical map correlation-based comparison of EEG recorded inside and outside the scanner (26). Simultaneous icEEG and fMRI (icEEG-fMRI) has been performed following extensive safety testing and locally adapted protocol (27–30) [see (31) for review], revealing IED-related BOLD changes local and remote from the relevant intracranial electrodes (32, 33).

We used icEEG-fMRI to explore BOLD changes and their distribution at the whole brain level for different IED classes on icEEG, in patients with refractory focal epilepsy who had good postsurgical outcome and well-characterized EZ. Surgical resection volume in this group of patients with good postsurgical represents confirmed identification of the EZ [where EZ is area of the brain deemed necessary to be resected to render patient seizure free (1, 34)].

Our hypotheses were as follows: (1) widespread BOLD networks, involving the surgical resection volume and remote brain areas, can be seen for IEDs on icEEG; (2) distribution of BOLD changes in surgical resection volume and remote brain areas is different for IED classes based on their topographic localization and their relationship with surgical resection volume. We investigated the anatomical localization and level of concordance of IED-related BOLD maps with the surgical resection volume. We also evaluated the relationship between different IED classes and level of concordance of BOLD maps with the surgical resection volume, and different IED classes and presence of BOLD changes in remote healthy cortex and other brain areas.

Methods

Eight patients with refractory focal epilepsy had icEEG-fMRI during their invasive pre-surgical evaluation, who had subsequently undergone resective epilepsy surgery with a good postsurgical outcome, i.e., completely seizure free or only auras (ILAE class I or II outcome) for more than 2 years after surgery. All patients gave written informed consent. The study was approved by the joint research ethics committee of the National Hospital for Neurology and Neurosurgery, Queen Square, London (UCLH NHS Foundation Trust) and UCL Institute of Neurology, Queen Square, London, UK.

Clinical Background

Prior to implantation, all patients had undergone detailed clinical history and examination, a structural MRI as per protocol specifically designed for epilepsy (35), long-term scalp video-EEG monitoring, neuropsychological and neuropsychiatric assessments, and additional functional imaging tests including positron emission tomography (PET), magnetoencephalography (MEG), or ictal single photon emission computed tomography (ictal SPECT) as indicated (see Table 1).

In accordance with routine clinical practice at our center, implantation of intracranial electrodes was guided by a hypothesis-based consensus decision generated from the results of non-invasive investigations. The SOZ, EZ, and the extent of surgical resection (Table 2) were defined by experienced Clinical Neurophysiologists/Epileptologists (BD, TW, and MW) and members of the multidisciplinary team based on invasive (multiple grid/depth electrode contacts on icEEG) and non-invasive investigations. The implantation scheme for each patient is shown in Supplementary Table 1. Post-surgical outcome (Table 2) was assessed with the ILAE classification (36).

Intracranial EEG-fMRI Acquisition

After the clinical icEEG recordings were completed, the implanted electrodes (numbering between 56 and 128 contacts) were connected to magnetic resonance scanner-compatible cables and amplifier system (32) for icEEG-fMRI acquisition. icEEG was recorded, processed online (to reduce the scanner-related artifacts), and displayed (BrainVision Recorder, Brain Products, Germany) during the fMRI scanning.

In accordance with our icEEG-fMRI protocol (29) echo planar images (EPI: TR/TE/flip angle = 3,000 ms/40 ms/90°, 64 × 64 acquisition matrix, 38 × 2.5 mm slices with a 0.5-mm gap) were acquired using a 1.5-T Siemens Avanto scanner (Erlangen, Germany) with a standard transmit/receive head coil and low specific absorption rate sequences ( ≤ 0.1 W/kg, head average) to reduce the risk of health hazards. One (for patients #2 and 4) or two (for patients #1, 3, 5, 6, 7, and 8) 10-min resting-state EPI time series (depending on patient comfort inside the scanner and time constraints) and T1-weighted structural scans were acquired.

Intracranial EEG Pre-processing and IED Classification

icEEG recorded during fMRI was corrected offline for scanner-related artifact (37) and reviewed by expert users (UC and MC) to identify and classify all IED using BrainVision Analyzer2 (Brain Products GmbH, Germany) and compared with clinical long-term icEEG recording and reports.

The identified IEDs were classified for the purpose of fMRI modeling according to the topographic distribution and localization (Table 2). For this topographic scheme, IED were classified according to the number of electrodes involved, their spatial location, field extent, and propagation (38) (Figure 1, Table 2) as either Focal: if they involved 2–4 contiguous electrode contacts and had similar field; Regional: if they involved 5–10 contiguous electrode contacts that may span up to two gyri; Widespread: if they involved more than 10 contiguous electrode contacts; or Non-contiguous: if they had a focal or regional field but also propagated to non-contiguous electrode contacts. The Regional, Widespread, and Non-contiguous classes taken together formed the non-focal IED class.

Figure 1

Furthermore, to assess our second hypothesis, we divided IEDs for their relationship with the surgical resection volume (i.e., confirmed EZ): IED classes overlapping the surgical resection volume were identified as IZ1 (primary irritative zone) and IED classes outside the surgical resection volume were identified as IZ2 (secondary irritative zone) (39) (see Table 2).

For patients #3 and #7 one of the two icEEG-fMRI sessions had to be excluded: for patient #3, icEEG had scanning-related artifacts and patient #7 had a subclinical seizure during one of the sessions (40).

fMRI Processing and Modeling

The fMRI data were analyzed using Statistical Parametric Mapping 8 (www.fil.ion.ucl.ac.uk) after discarding the first two volumes to account for the T1-saturation effect. Functional imaging data were corrected for slice acquisition time, realigned to the mean, and spatially smoothed using an isotropic Gaussian kernel of 8-mm FWHM (41).

A general linear model (GLM) was built to map IED-related hemodynamic changes. For patients who underwent two EPI series, these were included in a single GLM as separate sessions. Each IED was represented either as a zero-duration stick function (individual IED) or blocks (runs of IED). Each IED class was modeled as a separate effect and the corresponding time series of stick functions or blocks convolved with the canonical hemodynamic response function and its temporal and dispersion derivatives. In line with previous analyses, 24 inter-scan realignment parameters [6 realignment parameters from image pre-processing and a Volterra-expansion of these (42)] were included in the GLM as confounds to account for motion-related effects similar to our previous work (40).

Assessment of IED-Related BOLD Changes

For each IED class, the presence of significant BOLD clusters was assessed over whole brain using SPM[F]-maps at a statistical threshold of p < 0.001 (uncorrected for family-wise error) and a cluster threshold of five contiguous voxels, as in previous studies from our group and others (15, 22, 26, 32, 40, 43). The resulting SPMs were co-registered with pre- and post-surgical T1-weighted MRI scans using rigid-body registration in SPM. The localization of BOLD clusters, for each IED class, was visually assessed in relation to the surgical resection volume. Clusters of activity were also assessed in remote areas including cuneus, precuneus, cingulate gyrus, medial frontal lobe, and other brain areas such as basal ganglia, thalamus. The fitted BOLD time course for each cluster was plotted and classified as increases, decreases, or biphasic (consisting of both increases and decreases) according to the sign of the peak change relative to baseline.

The concordance of BOLD maps with surgical resection volume (i.e., confirmed EZ) was assessed for each IED class on icEEG using a concordance classification scheme in line with our previous work (15, 22, 26, 40, 43) as either:

• Entirely concordant: All BOLD clusters overlapping with/or located within 2 cm of the surgical resection volume in the same lobe.

• Concordant plus: The global statistical maximum BOLD cluster (GM cluster) overlapping with/or located within 2 cm of the surgical resection volume in the same lobe and other clusters were remote (i.e., >2 cm away in different lobe or opposite hemisphere) from the surgical resection volume.

• Some concordance: The GM cluster was remote from the surgical resection volume and at least one of the other clusters overlapped with or located within 2 cm of the surgical resection volume in the same lobe.

• Discordant: all clusters were remote, i.e., more than 2 cm from the surgical resection volume in the same lobe or were in a different lobe or opposite hemisphere from the surgical resection volume.

BOLD clusters confined to the ventricular system, vascular tree, edges, and base of brain and cerebellum were not considered further in this analysis (20, 43–46).

We performed chi-square tests (χ²) (SPSS Statistics) to assess the association between (1) topographic IED classes and level of concordance of BOLD maps/presence of BOLD clusters in remote brain areas/presence of balloon cells in patients with FCD; (2) IZ1/IZ2 IED classes and level of concordance of BOLD maps/presence of BOLD clusters in remote brain areas/presence of balloon cells in patients with FCD; (3) presence of balloon cells in FCD patients and level of concordance of BOLD maps.

Results

The clinical details for eight patients fulfilling selection criteria are summarized in Table 1. There were six males; the median age at the time of icEEG-fMRI was 32 years and the median age at seizure onset was 7.5 years. Six patients had frontal lobe epilepsy, one had temporal lobe epilepsy, and one had parietal lobe epilepsy. The median follow-up time with ILAE class I/II postsurgical outcome was 6 years.

Classification of IED

All patients had a mixture of different topographic IED classes (Table 2). The number of IED classes in any given patient ranged between 3 and 7 (median 4.5). There was a total of 38 IED classes across the group. Out of these 38 IED classes, 22 were Focal, and 16 were non-focal: Regional = 9, Widespread = 2, and Non-contiguous = 5 according to topographic classification scheme (see Table 2). In terms of irritative zones, 29/38 IEDs originated from IZ1 (Focal = 18, Regional = 8, Widespread = 2, and Non-contiguous = 1) and 9/38 from IZ2 (Focal = 4, Regional = 1, Widespread = 0, and Non-contiguous = 4).

Distribution of IED-Related BOLD Changes for IED Classes

All IED classes were associated with significant BOLD clusters (Table 3) that were both co-located with recording electrodes but also in regions remote from them (Figure 1). BOLD clusters were seen within the surgical resection volume (Concordant) in 71% (27/38) of IED-related BOLD maps. The cluster of concordance corresponded to the statistical global maxima in 8 maps (Entirely Concordant = 1, Concordant plus = 7) and to the second or other significant cluster in 19 maps (Some concordance). At least two maps were concordant in every patient, with a mean of 70% of the maps per patient being concordant (range 50–100%).

All maps except one contained more than one BOLD cluster (see Table 3). Across the group, BOLD clusters were distributed in the ipsi/contralateral hemisphere remote cortical or other brain areas including precuneus, medial superior frontal gyrus, cingulate, basal ganglia, and thalamus in 32/38 (84%) of BOLD maps (see Table 3).

Relationship With IED Topographic Classification

BOLD maps for non-focal IEDs (Regional, Non-contiguous, and Widespread) were more commonly concordant with surgical resection volume (Entirely Concordant = 1, Concordant plus = 5, Some Concordance = 8; 14/16, 87.5%) than for Focal IEDs (Concordant plus = 2, Some Concordance = 11; 13/22, 59%), (χ² = 7.08, p < 0.05). Presence of BOLD clusters in remote cortical and/or other brain areas, i.e., precuneus, medial frontal, cingulate, and thalamus, was more frequent for non-focal IED maps (15/16, 93%: Regional = 8, Non-contiguous = 5, and Widespread = 2) as compared to Focal IED maps (17/22, 77%) but the difference did not reach statistical significance.

Relationship With Irritative Zones

BOLD maps were concordant with the surgical resection volume for 20/29 (68%) IED classes from IZ1 (Entirely Concordant = 1; Concordant plus: 6; Some concordance: 13), and 7/9 (77%) IZ2 IED (Concordant plus: 1; Some concordance: 6). The map's level of concordance or presence of BOLD changes in remote cortical and/or other brain areas did not differ significantly between IZ1 and IZ2 IED classes.

Structural Abnormalities

Seven patients had structural abnormalities seen on MRI (Table 1). In six patients, these were in the EZ [focal cortical dysplasia (FCD) = 5, hippocampal sclerosis (HS) = 1], and one patient had an incidental finding of HS unrelated to the EZ.

In the subgroup of patients with FCD, four patients had FCD type IIB with balloon cells and one patient had FCD type IIA with no balloon cells (see Table 2). All patients with FCD showed at least one map with a BOLD cluster overlying the lesion. We did not find a statistically significant association between presence/absence of balloon cells and different IED classes and level of concordance of BOLD maps in this small subgroup of patients.

The patient with hippocampal sclerosis (#7) did not show any BOLD cluster directly overlying HS; however, BOLD clusters were seen within 2 cm of the structural abnormality.

Discussion

Scalp EEG-fMRI studies have shown that IED-related BOLD changes in EZ can predict good postsurgical outcome (9, 13, 14, 22–25, 47–49). One of the limiting factors for these studies has been low sensitivity of scalp EEG to capture the whole spectrum of epileptiform activity that can be revealed by invasive recordings (50, 51). Therefore, what is considered the baseline (“non-epileptic” state) in scalp EEG-fMRI studies must in fact contain a significant amount of epileptic discharges. Simultaneous icEEG-fMRI allows us to overcome this problem by exploring whole-brain changes for epileptiform discharges recorded directly from the cortex using icEEG. This study revealed significant BOLD signal changes for a wide range of IEDs using simultaneous icEEG-fMRI. Furthermore, we found that:

• Significant BOLD clusters for IEDs on icEEG were localized both within the surgical resection volume and remote cortical and other brain areas;

• More than 70% of IED classes showed BOLD maps concordant with the surgical resection volume, where BOLD clusters were seen within the surgical resection volume;

• IED with wider topographic distribution and localization: non-focal IED classes on icEEG were associated with the presence of BOLD clusters within the surgical resection volume.

Previous studies using icEEG-fMRI have shown BOLD changes related to IEDs and seizures (32, 33, 40, 52, 53). Comparison of visual and automated IED classification on icEEG (53) presented a more objective interpretation of icEEG, but there was no statistically significant difference in concordance of the BOLD maps for two IED classification techniques. The relationship of BOLD clusters in surgical resection volume and in remote brain areas for different IED classes has not been explored in previous studies. For this study, we think that visual classification of IEDs based on their topographic distribution and localization and IZs, which reflects clinical insight of the expert user, facilitates clinical interpretation of resulting BOLD maps. We compared distribution of BOLD clusters in surgical resection volume (i.e., level of concordance) and remote brain areas for different IED classes using icEEG-fMRI in the largest group of patients to date who had undergone epilepsy surgery and had a good postsurgical outcome with a long follow-up time (median: 6 years). The surgical resected volume can be rendered confirmed EZ considering long postsurgical seizure freedom.

Methodological Considerations

The feasibility and safety of simultaneous icEEG-fMRI has been established (27, 28, 30, 31). Signal degradation can be observed within up to 1 cm (often less at 1.5 T as in this study) of the electrode contacts and is orientation dependent (27), therefore limiting interpretation of the BOLD maps in the electrode contact's immediate vicinity. However, BOLD effects are generally more widespread (27, 54). BOLD maps revealed significant clusters for different IED classes on icEEG, which were concordant with the surgical resection volume (confirmed EZ) and other non-invasive and invasive investigations, and were also seen in distant areas known to be related to resting-state networks associated with interictal discharges (55, 56). Therefore, it is unlikely that these changes are false positive. Also, icEEG has high sensitivity to show IEDs from smaller generators as compared to scalp EEG (50, 51) and, thus, may be associated with relatively weaker BOLD changes from smaller brain regions; this is in line with previous icEEG-fMRI studies (29, 32, 33, 40).

IEDs were represented as single or series of events in separate regressors for each different class to evaluate specific BOLD pattern in a GLM framework (12, 22, 24, 32) using a standard hemodynamic response function and its derivatives as a hemodynamic kernel (41) to account for a degree of variability in hemodynamic peak delay and duration (12, 24).

In patients with a single seizure onset zone, there may be more than one IED class reflecting different topographic localization and distribution and IZs and not all of these require removal for good surgical outcome. Our interpretation of the IED classes took into account spatial localization, field distribution, propagation, and their relationship with the EZ (1, 38, 39).

In the concordance classification scheme, the first two levels of concordance—Entirely concordant and Concordant plus—are defined based on the location of GM cluster overlapping with or within 2 cm of the surgical resection volume ± presence of other BOLD clusters. For BOLD maps with Some concordance, a cluster other than GM cluster was overlapping with or within 2 cm of the surgical resection volume. In this retrospective study, confirmed EZ was known, and this other cluster in the surgical resection volume was identified. However, during prospective pre-surgical evaluation of patients with a presumed EZ, this cluster can be identified by a consensus agreement, for example, if it is concordant with the structural lesion such as FCD and/or other non-invasive/invasive localization techniques. Our choice of 2 cm as a distance threshold (within a single lobe) to ascertain concordant BOLD clusters reflects the uncertainties associated with implantation and co-registration-related brain shift and the anticipated spatial dislocation of two classes of signals due to neurovascular coupling (57, 58). We evaluated the level of concordance of IED-related BOLD maps irrespective of sign of BOLD change, as both BOLD increases and decreases can be found in the EZ (20, 22, 59, 60).

Neurophysiological and Neurobiological Significance

BOLD clusters were seen in multiple areas for all IEDs on icEEG, and these areas included surgical resection volume (i.e., confirmed EZ) and adjacent/remote apparently healthy cortex. This suggests the possibility of common underlying brain areas or networks recruited as propagation nodes or even generators for different IED classes (50, 51, 61–63), or these widespread BOLD changes may be secondary to extensive underlying pathology (64). We suggest that BOLD changes in cortex and other brain areas remote from the surgical resection volume (i.e., confirmed EZ) may represent propagated epileptic activity in agreement with scalp EEG-fMRI (16, 22, 24, 65, 66) and electric source imaging studies (67). Also, this propagated epileptic activity in remote cortical or other brain areas such as precuneus, medial frontal, cuneus, and thalamus may represent an interaction with resting-state networks in line with previous scalp EEG-fMRI studies (55, 56), which can have implications on level of cognition and consciousness (55, 68) at some level and grants further research. We noted that changes in these areas that are part of default mode network were deactivations, but activations were also seen for some IEDs (see Table 3). Though a complete picture of underlying neuronal activity for IEDs may not be seen on icEEG (48) due to its limited spatial sampling, it is difficult to further elucidate whether these activations represent propagation of epileptic activity and deactivations represent involvement of default mode network. Future investigations correlating IEDs on icEEG with topographic maps of IEDs on scalp EEG and their associated BOLD changes will be required to understand the full pathologic nature of such networks.

The presence of BOLD clusters in surgical resection volume (confirmed EZ), as reflected by level of concordance of BOLD maps, was associated with topographic and field distribution of IED on icEEG. Non-focal IEDs on icEEG with wider topographic and field distribution (Regional, Widespread, and Non-contiguous) showed BOLD clusters in surgical resection volume more commonly, compared to focal IEDs on icEEG. This finding is similar to a recent scalp EEG-fMRI study (66) in which widespread epileptic discharges were more likely to show BOLD activation in seizure onset areas. The significance of this finding raises interesting questions about the BOLD effect, for example: is there a spatial scale of neural activity below which the strength of the BOLD change reflects only the local intensity of that activity, in contrast to its spatial extent? Hemodynamic changes may be limited to the activation of a minimum neuronal volume and its synchronization on EEG (69); this may explain the more common presence of BOLD changes in surgical resection volume and remote areas for IEDs with more widespread field extent. It is possible that signal dropout in the local vicinity of icEEG contact (28, 29) can limit to show BOLD change for IED with a very focal field extent. Future imaging sequence development with less signal dropout around implanted electrodes may be able to localize BOLD changes for very focal IEDs on icEEG. In addition, duration of underlying field potentials for epileptic discharges reflected by the sharp wave width can also affect amplitude of the BOLD signal (70), and event parameterization (amplitude, frequency content and duration) may be a useful way forward to further investigate BOLD changes for IED on icEEG.

Clinical Significance

We found that icEEG-fMRI has greater sensitivity: all patients showed IED-related BOLD changes, whereas previously published scalp EEG-fMRI studies have shown IED-related BOLD changes in 30–78% of patients (12, 24, 26). We suggest that this partly reflects the high sensitivity, specificity, and spatial resolution of icEEG (when placed judiciously) compared to scalp EEG (1) and the possibility that this results in more accurate definition of the BOLD baseline. Also, there could possibly be selection bias; patients with a clear focus and IEDs on scalp EEG are more likely to proceed for invasive icEEG.

The strength of our data is that the surgical resection volume represents confirmed EZ as reflected by long seizure freedom after surgery (1). This level of confidence is lacking in previous studies. We found BOLD clusters located in the surgical resection volume in 70% of the maps for different IED classes on icEEG. As icEEG-fMRI can reveal BOLD network across the whole brain and does not suffer from limited spatial sampling of icEEG, it is possible that BOLD clusters remote from surgical resection volume may represent other generator or propagator areas of epileptic activity that are not covered by icEEG. Small sample size and heterogenous underlying pathology could be considered limitations of this study, restraining generalized application of these findings to all patients undergoing epilepsy surgery. It will be interesting to compare in the future, in larger sample size, if there is any difference of BOLD patterns for IEDs on icEEG between seizure-free patients and patients who did not achieve seizure freedom after epilepsy surgery, and if it can inform epilepsy surgery approach.

In conclusion, icEEG-fMRI studies constitute a significant step toward the better understanding of hemodynamic changes related to epileptic activity. It can provide localization of BOLD network at whole-brain level with high sensitivity for different classes of interictal discharges on icEEG originating from focal areas. In addition, BOLD clusters in surgical resection volume (confirmed EZ) were seen more commonly for non-focal epileptiform discharges on icEEG.

Funding

This work was partly funded through grants and bursaries from the Medical Research Council (MRC Grant No. G0301067), Action Medical Research, Swiss National Science Foundation (SNF grant 320030-141165 and 33CM30-140332, SPUM Epilepsy), University of Modena, Reggio Emilia, and UCL Institute of Neurology. This work was undertaken at UCLH/UCL, which received a proportion of funding from the Department of Health's NIHR Biomedical Research Centers funding scheme.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

• 1.

  LudersHONajmINairDWiddess-WalshPBingmanW. The epileptogenic zone: general principles. Epileptic Disord. (2006) 8(Suppl. 2):S1–9.

  • Google Scholar

• 2.

  BinnieCDStefanH. Modern electroencephalography: its role in epilepsy management. Clin Neurophysiol. (1999) 110:1671–97. 10.1016/S1388-2457(99)00125-X

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3.

  SmithSJ. EEG in the diagnosis, classification, and management of patients with epilepsy. J Neurol Neurosurg Psychiatry. (2005) 76(Suppl. 2): ii2-ii7. 10.1136/jnnp.2005.069245.0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4.

  RemiJVollmarCdeMAHeinlinJPeraudANoachtarS. Congruence and discrepancy of interictal and ictal EEG with MRI lesions in focal epilepsies. Neurology. (2011) 77:1383–90. 10.1212/WNL.0b013e31823152c3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5.

  CatarinoCBVollmarCNoachtarS. Paradoxical lateralization of non-invasive electroencephalographic ictal patterns in extra-temporal epilepsies. Epilepsy Res. (2011) 99:147–55. 10.1016/j.eplepsyres.2011.11.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6.

  LachauxJPRudraufDKahaneP. Intracranial EEG and human brain mapping. J Physiol Paris. (2003) 97:613–628. 10.1016/j.jphysparis.2004.01.018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7.

  WellmerJvon derGFKlarmannUWeberCElgerCEUrbachHet al. Risks and benefits of invasive epilepsy surgery workup with implanted subdural and depth electrodes. Epilepsia. (2012) 53:1322–32. 10.1111/j.1528-1167.2012.03545.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8.

  CardinaleFRizziMVignatiECossuMCastanaLd'OrioPet al. Stereoelectroencephalography: retrospective analysis of 742 procedures in a single centre. Brain. (2019) 142:2688–704. 10.1093/brain/awz196

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9.

  Caballero-GaudesCVan de VilleDGrouillerFThorntonRLemieuxLSeeckMet al. Mapping interictal epileptic discharges using mutual information between concurrent EEG and fMRI. Neuroimage. (2013) 68:248–62. 10.1016/j.neuroimage.2012.12.011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10.

  DonaireACapdevilaACarrenoMSetoainXRumiaJAparicioJet al. Identifying the cortical substrates of interictal epileptiform activity in patients with extratemporal epilepsy: an EEG-fMRI sequential analysis and FDG-PET study. Epilepsia. (2013) 54:678–90. 10.1111/epi.12091

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11.

  PittauFDubeauFGotmanJ. Contribution of EEG/fMRI to the definition of the epileptic focus. Neurology. (2012) 78:1479–87. 10.1212/WNL.0b013e3182553bf7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12.

  Salek-HaddadiADiehlBHamandiKMerschhemkeMListonAFristonKet al. Hemodynamic correlates of epileptiform discharges: an EEG-fMRI study of 63 patients with focal epilepsy. Brain Res. (2006) 1088:148–66. 10.1016/j.brainres.2006.02.098

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13.

  ZijlmansMHuiskampGHersevoortMSeppenwooldeJHvan HuffelenACLeijtenFS. EEG-fMRI in the preoperative work-up for epilepsy surgery. Brain. (2007) 130:2343–53. 10.1093/brain/awm141

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14.

  CoanACChaudharyUJFredericGCamposBMPeraniSDeCAet al. EEG-fMRI in the presurgical evaluation of temporal lobe epilepsy. J Neurol Neurosurg Psychiatry. (2016) 87:642–9. 10.1136/jnnp-2015-310401

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15.

  MarkoulaSChaudharyUJPeraniSDeCAYadeeTDuncanJSet al. The impact of mapping interictal discharges using EEG-fMRI on the epilepsy presurgical clinical decision making process: a prospective study. Seizure. (2018) 61:30–7. 10.1016/j.seizure.2018.07.016

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16.

  PittauFFerriLFahoumFDubeauFGotmanJ. Contributions of EEG-fMRI to assessing the epileptogenicity of focal cortical dysplasia. Front Comput Neurosci. (2017) 11:8. 10.3389/fncom.2017.00008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17.

  DonaireABargalloNFalconCMaestroICarrenoMSetoainJet al. Identifying the structures involved in seizure generation using sequential analysis of ictal-fMRI data. Neuroimage. (2009) 47:173–83. 10.1016/j.neuroimage.2009.03.059

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18.

  ThorntonRCRodionovRLaufsHVulliemozSVaudanoACarmichaelDet al. Imaging haemodynamic changes related to seizures: comparison of EEG-based general linear model, independent component analysis of fMRI and intracranial EEG. Neuroimage. (2010) 53:196–205. 10.1016/j.neuroimage.2010.05.064

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19.

  TyvaertLLeVanPDubeauFGotmanJ. Noninvasive dynamic imaging of seizures in epileptic patients. Hum Brain Mapp. (2009) 30:3993–4011. 10.1002/hbm.20824

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20.

  ChaudharyUJCarmichaelDWRodionovRThorntonRCBartlettPVulliemozSet al. Mapping pre-ictal and ictal haemodynamic networks using video-electroencephalography and functional imaging. Brain. (2012) 10.1093/brain/aws302

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21.

  Sierra-MarcosAMaestroIFalconCDonaireASetoainJAparicioJet al. Ictal EEG-fMRI in localization of epileptogenic area in patients with refractory neocortical focal epilepsy. Epilepsia. (2013) 54:1688–98. 10.1111/epi.12329

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22.

  ThorntonRVulliemozSRodionovRCarmichaelDWChaudharyUJDiehlBet al. Epileptic networks in focal cortical dysplasia revealed using electroencephalography-functional magnetic resonance imaging. Ann Neurol. (2011) 70:822–37. 10.1002/ana.22535

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23.

  Al-AsmiABenarCGGrossDWKhaniYAAndermannFPikeBet al. fMRI activation in continuous and spike-triggered EEG-fMRI studies of epileptic spikes. Epilepsia. (2003) 44:1328–39. 10.1046/j.1528-1157.2003.01003.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24.

  ThorntonRLaufsHRodionovRCannadathuSCarmichaelDWVulliemozSet al. EEG correlated functional MRI and postoperative outcome in focal epilepsy. J Neurol Neurosurg Psychiatry. (2010) 81:922–7. 10.1136/jnnp.2009.196253

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25.

  TyvaertLHawcoCKobayashiELeVanPDubeauFGotmanJ. Different structures involved during ictal and interictal epileptic activity in malformations of cortical development: an EEG-fMRI study. Brain. (2008) 131:2042–60. 10.1093/brain/awn145

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26.

  GrouillerFThorntonRCGroeningKSpinelliLDuncanJSSchallerKet al. With or without spikes: localization of focal epileptic activity by simultaneous electroencephalography and functional magnetic resonance imaging. Brain. (2011) 134:2867–86. 10.1093/brain/awr156

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27.

  CarmichaelDWThorntonJSRodionovRThorntonRMcEvoyAAllenPJet al. Safety of localizing epilepsy monitoring intracranial electroencephalograph electrodes using MRI: radiofrequency-induced heating. J Magn Reson Imaging. (2008) 28:1233–44. 10.1002/jmri.21583

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28.

  CarmichaelDWThorntonJSRodionovRThorntonRMcEvoyAWOrdidgeRJet al. Feasibility of simultaneous intracranial EEG-fMRI in humans: a safety study. Neuroimage. (2010) 49:379–90. 10.1016/j.neuroimage.2009.07.062

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29.

  CarmichaelDWVulliemozSRodionovRThorntonJSMcEvoyAWLemieuxL. Simultaneous intracranial EEG-fMRI in humans: protocol considerations and data quality. Neuroimage. (2012) 63:301–9. 10.1016/j.neuroimage.2012.05.056

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30.

  BoucousisSMBeersCACunninghamCJGaxiola-ValdezIPittmanDJGoodyearBGet al. Feasibility of an intracranial EEG-fMRI protocol at 3T: risk assessment and image quality. Neuroimage. (2012) 63:1237–48. 10.1016/j.neuroimage.2012.08.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31.

  HawsawiHBCarmichaelDW. Lemieux L. Safety of simultaneous scalp or intracranial EEG during MRI: a review. Front Phy. (2017) 5:42. 10.3389/fphy.2017.00042

  • CrossRef
  • Google Scholar

• 32.

  VulliemozSCarmichaelDWRosenkranzKDiehlBRodionovRWalkerMCet al. Simultaneous intracranial EEG and fMRI of interictal epileptic discharges in humans. Neuroimage. (2011) 54:182–90. 10.1016/j.neuroimage.2010.08.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33.

  CunninghamCBGoodyearBGBadawyRZaamoutFPittmanDJBeersCAet al. Intracranial EEG-fMRI analysis of focal epileptiform discharges in humans. Epilepsia. (2012) 53:1636–48. 10.1111/j.1528-1167.2012.03601.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34.

  RosenowFLudersH. Presurgical evaluation of epilepsy. Brain. (2001) 124:1683–700. 10.1093/brain/124.9.1683

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35.

  DuncanJS. Imaging in the surgical treatment of epilepsy. Nat Rev Neurol. (2010) 6:537–50. 10.1038/nrneurol.2010.131

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36.

  WieserHGBlumeWTFishDGoldensohnEHufnagelAKingDet al. ILAE commission report. Proposal for a new classification of outcome with respect to epileptic seizures following epilepsy surgery. Epilepsia. (2001) 42:282–6. 10.1046/j.1528-1157.2001.4220282.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37.

  AllenPJJosephsOTurnerR. A method for removing imaging artifact from continuous EEG recorded during functional MRI. Neuroimage. (2000) 12:230–9. 10.1006/nimg.2000.0599

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38.

  LudersHNoachtarSBensonJK. Atlas and Classification of Electroencephalography. Philadelphia, PA: Saunders (2000).

  • Google Scholar

• 39.

  BettusGRanjevaJPWendlingFBenarCGConfort-GounySRegisJet al. Interictal functional connectivity of human epileptic networks assessed by intracerebral EEG and BOLD signal fluctuations. PLoS ONE. (2011) 6:e20071. 10.1371/journal.pone.0020071

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40.

  ChaudharyUJCentenoMThorntonRCRodionovRVulliemozSMcEvoyAWet al. Mapping human preictal and ictal haemodynamic networks using simultaneous intracranial EEG-fMRI. Neuroimage Clin. (2016) 11:486–93. 10.1016/j.nicl.2016.03.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41.

  FristonKJHolmesAPWorsleyKJPolineJBFirthCDFrackowiakRSJ. Statistical parametric maps in functional imaging: a general linear approach. Hum Brain Mapp. (1995) 2:189–210. 10.1002/hbm.460020402

  • CrossRef
  • Google Scholar

• 42.

  FristonKJWilliamsSHowardRFrackowiakRSTurnerR. Movement-related effects in fMRI time-series. Magn Reson Med. (1996) 35:346–55. 10.1002/mrm.1910350312

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43.

  ChaudharyUJRodionovRCarmichaelDWThorntonRCDuncanJSLemieuxL. Improving the sensitivity of EEG-fMRI studies of epileptic activity by modelling eye blinks, swallowing and other video-EEG detected physiological confounds. Neuroimage. (2012) 61:1383–93. 10.1016/j.neuroimage.2012.03.028

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44.

  BirnRMCoxRWBandettiniPA. Experimental designs and processing strategies for fMRI studies involving overt verbal responses. Neuroimage. (2004) 23:1046–58. 10.1016/j.neuroimage.2004.07.039

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45.

  JansenMWhiteTPMullingerKJLiddleEBGowlandPAFrancisSTet al. Motion-related artefacts in EEG predict neuronally plausible patterns of activation in fMRI data. Neuroimage. (2012) 59:261–70. 10.1016/j.neuroimage.2011.06.094

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46.

  JohnstoneTOres WalshKSGreischarLLAlexanderALFoxASDavidsonRJet al. Motion correction and the use of motion covariates in multiple-subject fMRI analysis. Hum Brain Mapp. (2006) 27:779–88. 10.1002/hbm.20219

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47.

  AnDFahoumFHallJOlivierAGotmanJDubeauF. Electroencephalography/functional magnetic resonance imaging responses help predict surgical outcome in focal epilepsy. Epilepsia. (2013) 54:2184–94. 10.1111/epi.12434

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48.

  van HoudtPJDe MunckJCLeijtenFSHuiskampGJColonAJBoonPAet al. EEG-fMRI correlation patterns in the presurgical evaluation of focal epilepsy: a comparison with electrocorticographic data and surgical outcome measures. Neuroimage. (2013) 75:238–48. 10.1016/j.neuroimage.2013.02.033

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49.

  KowalczykMAOmidvarniaAAbbottDFTailbyCVaughanDNJacksonGD. Clinical benefit of presurgical EEG-fMRI in difficult-to-localize focal epilepsy: a single-institution retrospective review. Epilepsia. (2020) 61:49–60. 10.1111/epi.16399

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50.

  TaoJXRayAHawes-EbersoleSEbersoleJS. Intracranial EEG substrates of scalp EEG interictal spikes. Epilepsia. (2005) 46:669–76. 10.1111/j.1528-1167.2005.11404.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51.

  TaoJXBaldwinMHawes-EbersoleSEbersoleJS. Cortical substrates of scalp EEG epileptiform discharges. J Clin Neurophysiol. (2007) 24:96–100. 10.1097/WNP.0b013e31803ecdaf

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52.

  AghakhaniYBeersCAPittmanDJGaxiola-ValdezIGoodyearBGFedericoP. Co-localization between the BOLD response and epileptiform discharges recorded by simultaneous intracranial EEG-fMRI at 3 T. Neuroimage Clin. (2015) 7:755–63. 10.1016/j.nicl.2015.03.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53.

  SharmaNKPedreiraCChaudharyUJCentenoMCarmichaelDWYadeeTet al. BOLD mapping of human epileptic spikes recorded during simultaneous intracranial EEG-fMRI: the impact of automated spike classification. Neuroimage. (2019) 184:981–92. 10.1016/j.neuroimage.2018.09.065

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54.

  DeichmannRJosephsOHuttonCCorfieldDRTurnerR. Compensation of susceptibility-induced BOLD sensitivity losses in echo-planar fMRI imaging. Neuroimage. (2002) 15:120–35. 10.1006/nimg.2001.0985

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55.

  FahoumFZelmannRTyvaertLDubeauFGotmanJ. Epileptic discharges affect the default mode network–FMRI and intracerebral EEG evidence. PLoS ONE. (2013) 8:e68038. 10.1371/journal.pone.0068038

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56.

  LaufsHHamandiKSalek-HaddadiAKleinschmidtAKDuncanJSLemieuxL. Temporal lobe interictal epileptic discharges affect cerebral activity in “default mode” brain regions. Hum Brain Mapp. (2007) 28:1023–32. 10.1002/hbm.20323

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57.

  DisbrowEASlutskyDARobertsTPKrubitzerLA. Functional MRI at 1.5 tesla: a comparison of the blood oxygenation level-dependent signal and electrophysiology.Proc Natl Acad Sci USA. (2000) 97:9718–23. 10.1073/pnas.170205497

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58.

  NimskyCGanslandtOCernySHastreiterPGreinerGFahlbuschR. Quantification of, visualization of, and compensation for brain shift using intraoperative magnetic resonance imaging. Neurosurgery. (2000) 47:1070–9. 10.1097/00006123-200011000-00008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59.

  JacobsJKobayashiEBoorRMuhleHStephanWHawcoCet al. Hemodynamic responses to interictal epileptiform discharges in children with symptomatic epilepsy. Epilepsia. (2007) 48:2068–78. 10.1111/j.1528-1167.2007.01192.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60.

  PittauFFahoumFZelmannRDubeauFGotmanJ. Negative BOLD response to interictal epileptic discharges in focal epilepsy. Brain Topogr. (2013) 26:627–40. 10.1007/s10548-013-0302-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61.

  RayATaoJXHawes-EbersoleSMEbersoleJS. Localizing value of scalp EEG spikes: a simultaneous scalp and intracranial study. Clin Neurophysiol. (2007) 118:69–79. 10.1016/j.clinph.2006.09.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62.

  FlanaganDBadawyRAJacksonGD. EEG-fMRI in focal epilepsy: Local activation and regional networks. Clin Neurophysiol. (2013) 125:21–31. 10.1016/j.clinph.2013.06.182

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63.

  LaufsHRichardsonMPSalek-HaddadiAVollmarCDuncanJSGaleKet al. Converging PET and fMRI evidence for a common area involved in human focal epilepsies. Neurology. (2011) 77:904–10. 10.1212/WNL.0b013e31822c90f2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64.

  MoranNFLemieuxLKitchenNDFishDRShorvonSD. Extrahippocampal temporal lobe atrophy in temporal lobe epilepsy and mesial temporal sclerosis. Brain. (2001) 124:167–75. 10.1093/brain/124.1.167

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65.

  FahoumFLopesRPittauFDubeauFGotmanJ. Widespread epileptic networks in focal epilepsies: EEG-fMRI study. Epilepsia. (2012) 53:1618–27. 10.1111/j.1528-1167.2012.03533.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66.

  YamazoeTvonENKhooHMHuangYHZazubovitsNDubeauFet al. Widespread interictal epileptic discharge more likely than focal discharges to unveil the seizure onset zone in EEG-fMRI. Clin Neurophysiol. (2019) 130:429–38. 10.1016/j.clinph.2018.12.014

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 67.

  VulliemozSLemieuxLDaunizeauJMichelCMDuncanJS. The combination of EEG source imaging and EEG-correlated functional MRI to map epileptic networks. Epilepsia. (2010) 51:491–505. 10.1111/j.1528-1167.2009.02342.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68.

  ChaudharyUJCentenoMCarmichaelDWVollmarCRodionovRBonelliSet al. Imaging the interaction: Epileptic discharges, working memory, and behavior. Hum Brain Mapp. (2012) 34:2910–7. 10.1002/hbm.22115

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69.

  AbreuRLealALopes daSFFigueiredoP. EEG synchronization measures predict epilepsy-related BOLD-fMRI fluctuations better than commonly used univariate metrics. Clin Neurophysiol. (2018) 129:618–35. 10.1016/j.clinph.2017.12.038

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70.

  MurtaTHuLTierneyTMChaudharyUJWalkerMCCarmichaelDWet al. A study of the electro-haemodynamic coupling using simultaneously acquired intracranial EEG and fMRI data in humans. Neuroimage. (2016) 142:371–80. 10.1016/j.neuroimage.2016.08.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar