Skip to main content

ORIGINAL RESEARCH article

Front. Neurol., 21 September 2021
Sec. Epilepsy
This article is part of the Research Topic Advances and Applications of the EEG-fMRI Technique on Epilepsies View all 12 articles

Mapping Epileptic Networks Using Simultaneous Intracranial EEG-fMRI

  • 1Department of Clinical and Experimental Epilepsy, University College London (UCL) Institute of Neurology, National Hospital for Neurology and Neurosurgery, London, United Kingdom
  • 2Magnetic Resonance Imaging (MRI) Unit, Epilepsy Society, Chalfont St. Peter, United Kingdom
  • 3Neurology Department, University Hospital Coventry and Warwickshire, Coventry, United Kingdom
  • 4Epilepsy Unit, Neurology Department, Hospital Clinic Barcelona, Barcelona, Spain
  • 5Imaging and Biophysics Unit, University College London (UCL) Institute of Child Health, London, United Kingdom
  • 6Clinical Neurophysiology, National Hospital for Neurology and Neurosurgery, London, United Kingdom

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 (921). 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, 2325). 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 (2730) [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).

TABLE 1
www.frontiersin.org

Table 1. Clinical characteristics.

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).

TABLE 2
www.frontiersin.org

Table 2. Invasive localization and post-surgical outcome.

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
www.frontiersin.org

Figure 1. IED-related BOLD changes for Patient #2. (A) Implantation scheme of intracranial grid, strip, and depth electrodes and the invasively defined EZ is highlighted as red-colored electrodes on the sketched diagram and red square on the 3D rendered brain. (B) Representative sample of icEEG showing different IED classes based on their topographic distribution and localization included in the general linear model. (C) SPM[F] maps (p < 0.001) for different IED types overlaid on co-registered postsurgical T1-volume. (i) Focal IED [G4, 5 (n = 72); Concordant plus BOLD map]: changes in left inferior/middle frontal gyrus (GM cluster, within 2 cm of the surgical resection volume), left posterior temporal, and cuneus/precuneus. (ii) Focal IED [G12–15 (n = 29); BOLD map with Some concordance]: changes in left medial occipital (GM cluster), left inferior/middle frontal gyrus (within 2 cm of the surgical resection volume), left parieto-temporal, and cingulate gyrus. (iii) Regional IED [G4–6 + G12, 13 + G22–24 + G28–30 (n = 80); Discordant BOLD map]: changes in left inferior frontal gyrus (GM cluster 2.6 cm from the surgical resection volume), medial superior frontal gyrus, left temporo-occipital, and left middle frontal gyrus. (iv) Regional IED [G12–15 + G21–24 + DP2–4 (n = 350); Concordant plus BOLD map]: changes in left middle/inferior frontal gyrus (GM cluster, within 2 cm of the surgical resection volume), medial superior frontal gyrus, and supplementary motor area. (v) Widespread IED [G4–8 + G12–15 + G20–24 + G28–30 + DP2–4 (n = 244); Concordant plus BOLD map]: changes in left middle frontal gyrus (GM cluster, within 2 cm of the surgical resection volume), right middle temporal gyrus, posterior supplementary motor area, and cingulate/precuneus.

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, 4346).

We performed chi-square tests (χ2) (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%).

TABLE 3
www.frontiersin.org

Table 3. BOLD changes for individual IED class.

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%), (χ2 = 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, 2225, 4749). 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, 6163), 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.

Data Availability Statement

The datasets presented in this article are not readily available because subjects consent forms explicitly mention all information will be kept secure and strictly confidential, accessible only to those involved directly in this research. No information will be passed to any third parties or outside the EU for any reason without explicit consent of the subjects. Requests to access the datasets should be directed to Prof. Louis Lemieux, bG91aXMubGVtaWV1eCYjeDAwMDQwO3VjbC5hYy51aw==.

Ethics Statement

The studies involving human participants were reviewed and approved by 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. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

UC, MC, and DC made substantial contributions to the conception, design, data collection, data analysis, interpretation of results, and writing the article. BD made substantial contributions to the recruitment of patients, data interpretation, and revising the manuscript critically for important intellectual content. MW made substantial contributions to the data interpretation and revising the manuscript critically for important intellectual content. JD made substantial contributions to the recruitment of patients, interpretation of results, and revising the manuscript critically for important intellectual content. LL made substantial contributions to the conception, design, data analysis, and revising the manuscript critically for important intellectual content. All authors contributed to the article and approved the submitted version.

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.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

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.

Acknowledgments

We are thankful to Dr. Tim Wehner (Neurophysiologist and Epileptologist), Catherine Scott (Electrophysiologist), and the Department of Clinical Neurophysiology, National Hospital for Neurology and Neurosurgery, Queen Square, London, UK for their expert help for obtaining intracranial EEG. We thank Prof. Maria Thom at the Department of Neuropathology, Institute of Neurology, Queen Square for her expert help with the interpretation of neuropathological diagnoses. We are grateful to the radiographers Lisa Strycharczuk, Bruce Metheringham, and Alison Duncan, and the MR physicists Mark White and Laura Mancini of the Lysholm Department of Neuroradiology and Neurophysics at the National Hospital for Neurology and Neurosurgery (UCLH NHS Foundation Trust), for their expert scanning assistance.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2021.693504/full#supplementary-material

References

1. Luders HO, Najm I, Nair D, Widdess-Walsh P, Bingman W. The epileptogenic zone: general principles. Epileptic Disord. (2006) 8(Suppl. 2):S1–9.

Google Scholar

2. Binnie CD, Stefan H. Modern electroencephalography: its role in epilepsy management. Clin Neurophysiol. (1999) 110:1671–97. doi: 10.1016/S1388-2457(99)00125-X

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Remi J, Vollmar C, de MA, Heinlin J, Peraud A, Noachtar S. Congruence and discrepancy of interictal and ictal EEG with MRI lesions in focal epilepsies. Neurology. (2011) 77:1383–90. doi: 10.1212/WNL.0b013e31823152c3

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Catarino CB, Vollmar C, Noachtar S. Paradoxical lateralization of non-invasive electroencephalographic ictal patterns in extra-temporal epilepsies. Epilepsy Res. (2011) 99:147–55. doi: 10.1016/j.eplepsyres.2011.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Lachaux JP, Rudrauf D, Kahane P. Intracranial EEG and human brain mapping. J Physiol Paris. (2003) 97:613–628. doi: 10.1016/j.jphysparis.2004.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Wellmer J, von der GF, Klarmann U, Weber C, Elger CE, Urbach H, et al. Risks and benefits of invasive epilepsy surgery workup with implanted subdural and depth electrodes. Epilepsia. (2012) 53:1322–32. doi: 10.1111/j.1528-1167.2012.03545.x

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Cardinale F, Rizzi M, Vignati E, Cossu M, Castana L, d'Orio P, et al. Stereoelectroencephalography: retrospective analysis of 742 procedures in a single centre. Brain. (2019) 142:2688–704. doi: 10.1093/brain/awz196

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Caballero-Gaudes C, Van de Ville D, Grouiller F, Thornton R, Lemieux L, Seeck M, et al. Mapping interictal epileptic discharges using mutual information between concurrent EEG and fMRI. Neuroimage. (2013) 68:248–62. doi: 10.1016/j.neuroimage.2012.12.011

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Donaire A, Capdevila A, Carreno M, Setoain X, Rumia J, Aparicio J, et 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. doi: 10.1111/epi.12091

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Pittau F, Dubeau F, Gotman J. Contribution of EEG/fMRI to the definition of the epileptic focus. Neurology. (2012) 78:1479–87. doi: 10.1212/WNL.0b013e3182553bf7

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Salek-Haddadi A, Diehl B, Hamandi K, Merschhemke M, Liston A, Friston K, et al. Hemodynamic correlates of epileptiform discharges: an EEG-fMRI study of 63 patients with focal epilepsy. Brain Res. (2006) 1088:148–66. doi: 10.1016/j.brainres.2006.02.098

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Zijlmans M, Huiskamp G, Hersevoort M, Seppenwoolde JH, van Huffelen AC, Leijten FS. EEG-fMRI in the preoperative work-up for epilepsy surgery. Brain. (2007) 130:2343–53. doi: 10.1093/brain/awm141

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Coan AC, Chaudhary UJ, Frederic G, Campos BM, Perani S, De CA, et al. EEG-fMRI in the presurgical evaluation of temporal lobe epilepsy. J Neurol Neurosurg Psychiatry. (2016) 87:642–9. doi: 10.1136/jnnp-2015-310401

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Markoula S, Chaudhary UJ, Perani S, De CA, Yadee T, Duncan JS, et 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. doi: 10.1016/j.seizure.2018.07.016

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Pittau F, Ferri L, Fahoum F, Dubeau F, Gotman J. Contributions of EEG-fMRI to assessing the epileptogenicity of focal cortical dysplasia. Front Comput Neurosci. (2017) 11:8. doi: 10.3389/fncom.2017.00008

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Donaire A, Bargallo N, Falcon C, Maestro I, Carreno M, Setoain J, et al. Identifying the structures involved in seizure generation using sequential analysis of ictal-fMRI data. Neuroimage. (2009) 47:173–83. doi: 10.1016/j.neuroimage.2009.03.059

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Thornton RC, Rodionov R, Laufs H, Vulliemoz S, Vaudano A, Carmichael D, et 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. doi: 10.1016/j.neuroimage.2010.05.064

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Tyvaert L, LeVan P, Dubeau F, Gotman J. Noninvasive dynamic imaging of seizures in epileptic patients. Hum Brain Mapp. (2009) 30:3993–4011. doi: 10.1002/hbm.20824

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Chaudhary UJ, Carmichael DW, Rodionov R, Thornton RC, Bartlett P, Vulliemoz S, et al. Mapping pre-ictal and ictal haemodynamic networks using video-electroencephalography and functional imaging. Brain. (2012) doi: 10.1093/brain/aws302

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Sierra-Marcos A, Maestro I, Falcon C, Donaire A, Setoain J, Aparicio J, et al. Ictal EEG-fMRI in localization of epileptogenic area in patients with refractory neocortical focal epilepsy. Epilepsia. (2013) 54:1688–98. doi: 10.1111/epi.12329

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Thornton R, Vulliemoz S, Rodionov R, Carmichael DW, Chaudhary UJ, Diehl B, et al. Epileptic networks in focal cortical dysplasia revealed using electroencephalography-functional magnetic resonance imaging. Ann Neurol. (2011) 70:822–37. doi: 10.1002/ana.22535

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Al-Asmi A, Benar CG, Gross DW, Khani YA, Andermann F, Pike B, et al. fMRI activation in continuous and spike-triggered EEG-fMRI studies of epileptic spikes. Epilepsia. (2003) 44:1328–39. doi: 10.1046/j.1528-1157.2003.01003.x

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Thornton R, Laufs H, Rodionov R, Cannadathu S, Carmichael DW, Vulliemoz S, et al. EEG correlated functional MRI and postoperative outcome in focal epilepsy. J Neurol Neurosurg Psychiatry. (2010) 81:922–7. doi: 10.1136/jnnp.2009.196253

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Tyvaert L, Hawco C, Kobayashi E, LeVan P, Dubeau F, Gotman J. Different structures involved during ictal and interictal epileptic activity in malformations of cortical development: an EEG-fMRI study. Brain. (2008) 131:2042–60. doi: 10.1093/brain/awn145

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Grouiller F, Thornton RC, Groening K, Spinelli L, Duncan JS, Schaller K, et al. With or without spikes: localization of focal epileptic activity by simultaneous electroencephalography and functional magnetic resonance imaging. Brain. (2011) 134:2867–86. doi: 10.1093/brain/awr156

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Carmichael DW, Thornton JS, Rodionov R, Thornton R, McEvoy A, Allen PJ, et al. Safety of localizing epilepsy monitoring intracranial electroencephalograph electrodes using MRI: radiofrequency-induced heating. J Magn Reson Imaging. (2008) 28:1233–44. doi: 10.1002/jmri.21583

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Carmichael DW, Thornton JS, Rodionov R, Thornton R, McEvoy AW, Ordidge RJ, et al. Feasibility of simultaneous intracranial EEG-fMRI in humans: a safety study. Neuroimage. (2010) 49:379–90. doi: 10.1016/j.neuroimage.2009.07.062

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Carmichael DW, Vulliemoz S, Rodionov R, Thornton JS, McEvoy AW, Lemieux L. Simultaneous intracranial EEG-fMRI in humans: protocol considerations and data quality. Neuroimage. (2012) 63:301–9. doi: 10.1016/j.neuroimage.2012.05.056

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Boucousis SM, Beers CA, Cunningham CJ, Gaxiola-Valdez I, Pittman DJ, Goodyear BG, et al. Feasibility of an intracranial EEG-fMRI protocol at 3T: risk assessment and image quality. Neuroimage. (2012) 63:1237–48. doi: 10.1016/j.neuroimage.2012.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Hawsawi HB, Carmichael DW. Lemieux L. Safety of simultaneous scalp or intracranial EEG during MRI: a review. Front Phy. (2017) 5:42. doi: 10.3389/fphy.2017.00042

CrossRef Full Text | Google Scholar

32. Vulliemoz S, Carmichael DW, Rosenkranz K, Diehl B, Rodionov R, Walker MC, et al. Simultaneous intracranial EEG and fMRI of interictal epileptic discharges in humans. Neuroimage. (2011) 54:182–90. doi: 10.1016/j.neuroimage.2010.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Cunningham CB, Goodyear BG, Badawy R, Zaamout F, Pittman DJ, Beers CA, et al. Intracranial EEG-fMRI analysis of focal epileptiform discharges in humans. Epilepsia. (2012) 53:1636–48. doi: 10.1111/j.1528-1167.2012.03601.x

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Rosenow F, Luders H. Presurgical evaluation of epilepsy. Brain. (2001) 124:1683–700. doi: 10.1093/brain/124.9.1683

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Wieser HG, Blume WT, Fish D, Goldensohn E, Hufnagel A, King D, et al. ILAE commission report. Proposal for a new classification of outcome with respect to epileptic seizures following epilepsy surgery. Epilepsia. (2001) 42:282–6. doi: 10.1046/j.1528-1157.2001.4220282.x

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Allen PJ, Josephs O, Turner R. A method for removing imaging artifact from continuous EEG recorded during functional MRI. Neuroimage. (2000) 12:230–9. doi: 10.1006/nimg.2000.0599

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Luders H, Noachtar S, Benson JK. Atlas and Classification of Electroencephalography. Philadelphia, PA: Saunders (2000).

39. Bettus G, Ranjeva JP, Wendling F, Benar CG, Confort-Gouny S, Regis J, et al. Interictal functional connectivity of human epileptic networks assessed by intracerebral EEG and BOLD signal fluctuations. PLoS ONE. (2011) 6:e20071. doi: 10.1371/journal.pone.0020071

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Chaudhary UJ, Centeno M, Thornton RC, Rodionov R, Vulliemoz S, McEvoy AW, et al. Mapping human preictal and ictal haemodynamic networks using simultaneous intracranial EEG-fMRI. Neuroimage Clin. (2016) 11:486–93. doi: 10.1016/j.nicl.2016.03.010

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Friston KJ, Holmes AP, Worsley KJ, Poline JB, Firth CD, Frackowiak RSJ. Statistical parametric maps in functional imaging: a general linear approach. Hum Brain Mapp. (1995) 2:189–210. doi: 10.1002/hbm.460020402

CrossRef Full Text | Google Scholar

42. Friston KJ, Williams S, Howard R, Frackowiak RS, Turner R. Movement-related effects in fMRI time-series. Magn Reson Med. (1996) 35:346–55. doi: 10.1002/mrm.1910350312

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Chaudhary UJ, Rodionov R, Carmichael DW, Thornton RC, Duncan JS, Lemieux L. 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. doi: 10.1016/j.neuroimage.2012.03.028

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Birn RM, Cox RW, Bandettini PA. Experimental designs and processing strategies for fMRI studies involving overt verbal responses. Neuroimage. (2004) 23:1046–58. doi: 10.1016/j.neuroimage.2004.07.039

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Jansen M, White TP, Mullinger KJ, Liddle EB, Gowland PA, Francis ST, et al. Motion-related artefacts in EEG predict neuronally plausible patterns of activation in fMRI data. Neuroimage. (2012) 59:261–70. doi: 10.1016/j.neuroimage.2011.06.094

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Johnstone T, Ores Walsh KS, Greischar LL, Alexander AL, Fox AS, Davidson RJ, et al. Motion correction and the use of motion covariates in multiple-subject fMRI analysis. Hum Brain Mapp. (2006) 27:779–88. doi: 10.1002/hbm.20219

PubMed Abstract | CrossRef Full Text | Google Scholar

47. An D, Fahoum F, Hall J, Olivier A, Gotman J, Dubeau F. Electroencephalography/functional magnetic resonance imaging responses help predict surgical outcome in focal epilepsy. Epilepsia. (2013) 54:2184–94. doi: 10.1111/epi.12434

PubMed Abstract | CrossRef Full Text | Google Scholar

48. van Houdt PJ, De Munck JC, Leijten FS, Huiskamp GJ, Colon AJ, Boon PA, et 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. doi: 10.1016/j.neuroimage.2013.02.033

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Kowalczyk MA, Omidvarnia A, Abbott DF, Tailby C, Vaughan DN, Jackson GD. Clinical benefit of presurgical EEG-fMRI in difficult-to-localize focal epilepsy: a single-institution retrospective review. Epilepsia. (2020) 61:49–60. doi: 10.1111/epi.16399

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Tao JX, Ray A, Hawes-Ebersole S, Ebersole JS. Intracranial EEG substrates of scalp EEG interictal spikes. Epilepsia. (2005) 46:669–76. doi: 10.1111/j.1528-1167.2005.11404.x

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Tao JX, Baldwin M, Hawes-Ebersole S, Ebersole JS. Cortical substrates of scalp EEG epileptiform discharges. J Clin Neurophysiol. (2007) 24:96–100. doi: 10.1097/WNP.0b013e31803ecdaf

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Aghakhani Y, Beers CA, Pittman DJ, Gaxiola-Valdez I, Goodyear BG, Federico P. Co-localization between the BOLD response and epileptiform discharges recorded by simultaneous intracranial EEG-fMRI at 3 T. Neuroimage Clin. (2015) 7:755–63. doi: 10.1016/j.nicl.2015.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Sharma NK, Pedreira C, Chaudhary UJ, Centeno M, Carmichael DW, Yadee T, et al. BOLD mapping of human epileptic spikes recorded during simultaneous intracranial EEG-fMRI: the impact of automated spike classification. Neuroimage. (2019) 184:981–92. doi: 10.1016/j.neuroimage.2018.09.065

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Deichmann R, Josephs O, Hutton C, Corfield DR, Turner R. Compensation of susceptibility-induced BOLD sensitivity losses in echo-planar fMRI imaging. Neuroimage. (2002) 15:120–35. doi: 10.1006/nimg.2001.0985

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Fahoum F, Zelmann R, Tyvaert L, Dubeau F, Gotman J. Epileptic discharges affect the default mode network–FMRI and intracerebral EEG evidence. PLoS ONE. (2013) 8:e68038. doi: 10.1371/journal.pone.0068038

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Laufs H, Hamandi K, Salek-Haddadi A, Kleinschmidt AK, Duncan JS, Lemieux L. Temporal lobe interictal epileptic discharges affect cerebral activity in “default mode” brain regions. Hum Brain Mapp. (2007) 28:1023–32. doi: 10.1002/hbm.20323

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Disbrow EA, Slutsky DA, Roberts TP, Krubitzer LA. 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. doi: 10.1073/pnas.170205497

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Nimsky C, Ganslandt O, Cerny S, Hastreiter P, Greiner G, Fahlbusch R. Quantification of, visualization of, and compensation for brain shift using intraoperative magnetic resonance imaging. Neurosurgery. (2000) 47:1070–9. doi: 10.1097/00006123-200011000-00008

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Jacobs J, Kobayashi E, Boor R, Muhle H, Stephan W, Hawco C, et al. Hemodynamic responses to interictal epileptiform discharges in children with symptomatic epilepsy. Epilepsia. (2007) 48:2068–78. doi: 10.1111/j.1528-1167.2007.01192.x

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Pittau F, Fahoum F, Zelmann R, Dubeau F, Gotman J. Negative BOLD response to interictal epileptic discharges in focal epilepsy. Brain Topogr. (2013) 26:627–40. doi: 10.1007/s10548-013-0302-1

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Ray A, Tao JX, Hawes-Ebersole SM, Ebersole JS. Localizing value of scalp EEG spikes: a simultaneous scalp and intracranial study. Clin Neurophysiol. (2007) 118:69–79. doi: 10.1016/j.clinph.2006.09.010

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Flanagan D, Badawy RA, Jackson GD. EEG-fMRI in focal epilepsy: Local activation and regional networks. Clin Neurophysiol. (2013) 125:21–31. doi: 10.1016/j.clinph.2013.06.182

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Laufs H, Richardson MP, Salek-Haddadi A, Vollmar C, Duncan JS, Gale K, et al. Converging PET and fMRI evidence for a common area involved in human focal epilepsies. Neurology. (2011) 77:904–10. doi: 10.1212/WNL.0b013e31822c90f2

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Moran NF, Lemieux L, Kitchen ND, Fish DR, Shorvon SD. Extrahippocampal temporal lobe atrophy in temporal lobe epilepsy and mesial temporal sclerosis. Brain. (2001) 124:167–75. doi: 10.1093/brain/124.1.167

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Fahoum F, Lopes R, Pittau F, Dubeau F, Gotman J. Widespread epileptic networks in focal epilepsies: EEG-fMRI study. Epilepsia. (2012) 53:1618–27. doi: 10.1111/j.1528-1167.2012.03533.x

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Yamazoe T, von EN, Khoo HM, Huang YH, Zazubovits N, Dubeau F, et 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. doi: 10.1016/j.clinph.2018.12.014

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Vulliemoz S, Lemieux L, Daunizeau J, Michel CM, Duncan JS. The combination of EEG source imaging and EEG-correlated functional MRI to map epileptic networks. Epilepsia. (2010) 51:491–505. doi: 10.1111/j.1528-1167.2009.02342.x

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Chaudhary UJ, Centeno M, Carmichael DW, Vollmar C, Rodionov R, Bonelli S, et al. Imaging the interaction: Epileptic discharges, working memory, and behavior. Hum Brain Mapp. (2012) 34:2910–7. doi: 10.1002/hbm.22115

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Abreu R, Leal A, Lopes da SF, Figueiredo P. EEG synchronization measures predict epilepsy-related BOLD-fMRI fluctuations better than commonly used univariate metrics. Clin Neurophysiol. (2018) 129:618–35. doi: 10.1016/j.clinph.2017.12.038

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Murta T, Hu L, Tierney TM, Chaudhary UJ, Walker MC, Carmichael DW, et al. A study of the electro-haemodynamic coupling using simultaneously acquired intracranial EEG and fMRI data in humans. Neuroimage. (2016) 142:371–80. doi: 10.1016/j.neuroimage.2016.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: intracranial EEG, EEG-fMRI, IED/spikes, BOLD, post-surgical outcome

Citation: Chaudhary UJ, Centeno M, Carmichael DW, Diehl B, Walker MC, Duncan JS and Lemieux L (2021) Mapping Epileptic Networks Using Simultaneous Intracranial EEG-fMRI. Front. Neurol. 12:693504. doi: 10.3389/fneur.2021.693504

Received: 11 April 2021; Accepted: 20 July 2021;
Published: 21 September 2021.

Edited by:

Clarissa Lin Yasuda, State University of Campinas, Brazil

Reviewed by:

Umit Aydin, King's College London, United Kingdom
Marco Ciavarro, Istituto Neurologico Mediterraneo Neuromed (IRCCS), Italy
William B. Lo, Birmingham Children's Hospital, United Kingdom

Copyright © 2021 Chaudhary, Centeno, Carmichael, Diehl, Walker, Duncan and Lemieux. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Louis Lemieux, bG91aXMubGVtaWV1eCYjeDAwMDQwO3VjbC5hYy51aw==

Disclaimer: 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.