The spatio-temporal mapping of epileptic networks: combination of EEG-fMRI and EEG source imaging

Comparison of Resting-State EEG Network Analyses With and Without Parallel MRI in Genetic Generalized Epilepsy

Genetic generalized epilepsy (GGE) is conceptualized as a brain disorder involving distributed bilateral networks. To study these networks, simultaneous EEG-fMRI measurements can be used. However, inside-MRI EEG suffers from strong MR-related artifacts; it is not established whether EEG-based metrics in EEG-fMRI resting-state measurements are suitable for the analysis of group differences at source-level. We evaluated the impact of the inside-MR measurement condition on statistical group comparisons of EEG on source-level power and functional connectivity in patients with GGE versus healthy controls. We studied the cross-modal spatial relation of statistical group differences in seed-based FC derived from EEG and parallel fMRI. We found a significant increase in power and a frequency-specific change in functional connectivity for the inside MR-scanner compared to the outside MR-scanner condition. For power, we found reduced group difference between GGE and controls both in terms of statistical significance as well as effect size. Group differences for ImCoh remained similar both in terms of statistical significance as well as effect size. We found increased seed-based FC for GGE patients from the thalamus to the precuneus cortex region in fMRI, and in the theta band of simultaneous EEG. Our findings suggest that the analysis of EEG functional connectivity based on ImCoh is suitable for MR-EEG, and that relative group difference in a comparison of patients with GGE against controls are preserved. Spatial correspondence of seed-based FC group differences between the two modalities was found for the thalamus.

Camera-based Prospective Motion Correction in Paediatric Epilepsy Patients Enables EEG-fMRI Localization Even in High-motion States

BACKGROUND

EEG-fMRI is a useful additional test to localize the epileptogenic zone (EZ) particularly in MRI negative cases. However subject motion presents a particular challenge owing to its large effects on both MRI and EEG signal. Traditionally it is assumed that prospective motion correction (PMC) of fMRI precludes EEG artifact correction.

METHODS

Children undergoing presurgical assessment at Great Ormond Street Hospital were included into the study. PMC of fMRI was done using a commercial system with a Moiré Phase Tracking marker and MR-compatible camera. For retrospective EEG correction both a standard and a motion educated EEG artefact correction (REEGMAS) were compared to each other.

RESULTS

Ten children underwent simultaneous EEG-fMRI. Overall head movement was high (mean RMS velocity < 1.5 mm/s) and showed high inter- and intra-individual variability. Comparing motion measured by the PMC camera and the (uncorrected residual) motion detected by realignment of fMRI images, there was a five-fold reduction in motion from its prospective correction. Retrospective EEG correction using both standard approaches and REEGMAS allowed the visualization and identification of physiological noise and epileptiform discharges. Seven of 10 children had significant maps, which were concordant with the clinical EZ hypothesis in 6 of these 7.

CONCLUSION

To our knowledge this is the first application of camera-based PMC for MRI in a pediatric clinical setting. Despite large amount of movement PMC in combination with retrospective EEG correction recovered data and obtained clinically meaningful results during high levels of subject motion. Practical limitations may currently limit the widespread use of this technology.

EEG-fMRI of epileptiform discharges: non-invasive investigation of the whole brain

Simultaneous EEG-fMRI is a unique and non-invasive method for investigating epileptic activity. Interictal epileptiform discharge-related EEG-fMRI provides cortical and subcortical blood oxygen level-dependent (BOLD) signal changes specific to epileptic discharges. As a result, EEG-fMRI has revealed insights into generators and networks involved in epileptic activity in different types of epilepsy, demonstrating-for instance-the implication of the thalamus in human generalized spike and wave discharges and the role of the Default Mode Network (DMN) in absences and focal epilepsy, and proposed a mechanism for the cortico-subcortical interactions in Lennox-Gastaut syndrome discharges. EEG-fMRI can find deep sources of epileptic activity not available to scalp EEG or MEG and provides critical new information to delineate the epileptic focus when considering surgical treatment or electrode implantation. In recent years, methodological advances, such as artifact removal and automatic detection of events have rendered this method easier to implement, and its clinical potential has since been established by evidence of the impact of BOLD response on clinical decision-making and of the relationship between concordance of BOLD responses with extent of resection and surgical outcome. This review presents the recent developments in EEG-fMRI methodology and EEG-fMRI studies in different types of epileptic disorders as follows: EEG-fMRI acquisition, gradient and pulse artifact removal, statistical analysis, clinical applications, pre-surgical evaluation, altered physiological state in generalized genetic epilepsy, and pediatric EEG-fMRI studies.

Optimizing EEG Source Reconstruction with Concurrent fMRI-Derived Spatial Priors

Reconstructing EEG sources involves a complex pipeline, with the inverse problem being the most challenging. Multiple inversion algorithms are being continuously developed, aiming to tackle the non-uniqueness of this problem, which has been shown to be partially circumvented by including prior information in the inverse models. Despite a few efforts, there are still current and persistent controversies regarding the inversion algorithm of choice and the optimal set of spatial priors to be included in the inversion models. The use of simultaneous EEG-fMRI data is one approach to tackle this problem. The spatial resolution of fMRI makes fMRI derived spatial priors very convenient for EEG reconstruction, however, only task activation maps and resting-state networks (RSNs) have been explored so far, overlooking the recent, but already accepted, notion that brain networks exhibit dynamic functional connectivity fluctuations. The lack of a systematic comparison between different source reconstruction algorithms, considering potentially more brain-informative priors such as fMRI, motivates the search for better reconstruction models. Using simultaneous EEG-fMRI data, here we compared four different inversion algorithms (minimum norm, MN; low resolution electromagnetic tomography, LORETA; empirical Bayes beamformer, EBB; and multiple sparse priors, MSP) under a Bayesian framework (as implemented in SPM), each with three different sets of priors consisting of: (1) those specific to the algorithm; (2) those specific to the algorithm plus fMRI task activation maps and RSNs; and (3) those specific to the algorithm plus fMRI task activation maps and RSNs and network modules of task-related dFC states estimated from the dFC fluctuations. The quality of the reconstructed EEG sources was quantified in terms of model-based metrics, namely the expectation of the posterior probability P(model|data) and variance explained of the inversion models, and the overlap/proportion of brain regions known to be involved in the visual perception tasks that the participants were submitted to, and RSN templates, with/within EEG source components. Model-based metrics suggested that model parsimony is preferred, with the combination MSP and priors specific to this algorithm exhibiting the best performance. However, optimal overlap/proportion values were found using EBB and priors specific to this algorithm and fMRI task activation maps and RSNs or MSP and considering all the priors (algorithm priors, fMRI task activation maps and RSNs and dFC state modules), respectively, indicating that fMRI spatial priors, including dFC state modules, might contain useful information to recover EEG source components reflecting neuronal activity of interest. Our main results show that providing fMRI spatial derived priors that reflect the dynamics of the brain might be useful to map neuronal activity more accurately from EEG-fMRI. Furthermore, this work paves the way towards a more informative selection of the optimal EEG source reconstruction approach, which may be critical in future studies.

When Is Simultaneous Recording Necessary? A Guide for Researchers Considering Combined EEG-fMRI

Electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) provide non-invasive measures of brain activity at varying spatial and temporal scales, offering different views on brain function for both clinical and experimental applications. Simultaneous recording of these measures attempts to maximize the respective strengths of each method, while compensating for their weaknesses. However, combined recording is not necessary to address all research questions of interest, and experiments may have greater statistical power to detect effects by maximizing the signal-to-noise ratio in separate recording sessions. While several existing papers discuss the reasons for or against combined recording, this article aims to synthesize these arguments into a flow chart of questions that researchers can consider when deciding whether to record EEG and fMRI separately or simultaneously. Given the potential advantages of simultaneous EEG-fMRI, the aim is to provide an initial overview of the most important concepts and to direct readers to relevant literature that will aid them in this decision.

Localization of Epileptic Foci Based on Simultaneous EEG-fMRI Data

Combining functional magnetic resonance imaging (fMRI) and electroencephalography (EEG) enables a non-invasive investigation of the human brain function and evaluation of the correlation of these two important modalities of brain activity. This paper explores recent reports on using advanced simultaneous EEG–fMRI methods proposed to map the regions and networks involved in focal epileptic seizure generation. One of the applications of EEG and fMRI combination as a valuable clinical approach is the pre-surgical evaluation of patients with epilepsy to map and localize the precise brain regions associated with epileptiform activity. In the process of conventional analysis using EEG–fMRI data, the interictal epileptiform discharges (IEDs) are visually extracted from the EEG data to be convolved as binary events with a predefined hemodynamic response function (HRF) to provide a model of epileptiform BOLD activity and use as a regressor for general linear model (GLM) analysis of the fMRI data. This review examines the methodologies involved in performing such studies, including techniques used for the recording of EEG inside the scanner, artifact removal, and statistical analysis of the fMRI signal. It then discusses the results reported for patients with primary generalized epilepsy and patients with different types of focal epileptic disorders. An important matter that these results have brought to light is that the brain regions affected by interictal epileptic discharges might not be limited to the ones where they have been generated. The developed methods can help reveal the regions involved in or affected by a seizure onset zone (SOZ). As confirmed by the reviewed literature, EEG–fMRI provides information that comes particularly useful when evaluating patients with refractory epilepsy for surgery.

Augmenting interictal mapping with neurovascular coupling biomarkers by structured factorization of epileptic EEG and fMRI data

EEG-correlated fMRI analysis is widely used to detect regional BOLD fluctuations that are synchronized to interictal epileptic discharges, which can provide evidence for localizing the ictal onset zone. However, the typical, asymmetrical and mass-univariate approach cannot capture the inherent, higher order structure in the EEG data, nor multivariate relations in the fMRI data, and it is nontrivial to accurately handle varying neurovascular coupling over patients and brain regions. We aim to overcome these drawbacks in a data-driven manner by means of a novel structured matrix-tensor factorization: the single-subject EEG data (represented as a third-order spectrogram tensor) and fMRI data (represented as a spatiotemporal BOLD signal matrix) are jointly decomposed into a superposition of several sources, characterized by space-time-frequency profiles. In the shared temporal mode, Toeplitz-structured factors account for a spatially specific, neurovascular 'bridge' between the EEG and fMRI temporal fluctuations, capturing the hemodynamic response's variability over brain regions. By analyzing interictal data from twelve patients, we show that the extracted source signatures provide a sensitive localization of the ictal onset zone (10/12). Moreover, complementary parts of the IOZ can be uncovered by inspecting those regions with the most deviant neurovascular coupling, as quantified by two entropy-like metrics of the hemodynamic response function waveforms (9/12). Hence, this multivariate, multimodal factorization provides two useful sets of EEG-fMRI biomarkers, which can assist the presurgical evaluation of epilepsy. We make all code required to perform the computations available at https://github.com/svaneynd/structured-cmtf.

Dynamic analysis of fMRI activation during epileptic spikes can help identify the seizure origin

OBJECTIVE

We use the dynamic electroencephalography-functional magnetic resonance imaging (EEG-fMRI) method to incorporate variability in the amplitude and field of the interictal epileptic discharges (IEDs) into the fMRI analysis. We ask whether IED variability analysis can (a) identify additional activated brain regions during the course of IEDs, not seen in standard analysis; and (b) demonstrate the origin and spread of epileptic activity. We explore whether these functional changes recapitulate the structural connections and propagation of epileptic activity during seizures.

METHODS

Seventeen patients with focal epilepsy and at least 30 IEDs of a single type during simultaneous EEG-fMRI were studied. IED variability and EEG source imaging (ESI) analysis extracted time-varying dynamic changes. General linear modeling (GLM) generated static functional maps. Dynamic maps were compared to static functional maps. The dynamic sequence from IED variability was compared to the ESI results. In a subset of patients, we investigated structural connections between active brain regions using diffusion-based fiber tractography.

RESULTS

IED variability distinguished the origin of epileptic activity from its propagation in 15 of 17 (88%) patients. This included two cases where no result was obtained from the standard GLM analysis. In both of these cases, IED variability revealed activation in line with the presumed epileptic focus. Two cases showed no result from either method. Both had very high spike rates associated with dysplasia in the postcentral gyrus. In all 15 cases with dynamic activation, the observed dynamics were concordant with ESI. Fiber tractography identified specific white matter pathways between brain regions that were active at IED onset and propagation.

SIGNIFICANCE

Dynamic techniques involving IED variability can provide additional power for EEG-fMRI analysis, compared to standard analysis, revealing additional biologically plausible information in cases with no result from the standard analysis and gives insight into the origin and spread of IEDs.

Functional connectivity of the irritative zone identified by electrical source imaging, and EEG-correlated fMRI analyses

•

The irritative zone differs on Electrical Source Imaging (ESI) and EEG-fMRI.

•

Findings differ in functional connectivity and show low temporal correlation.

•

ESI and EEG-fMRI reveal distinct aspects of the irritative zone.

•

The differences may reflect differences in temporal resolution of the techniques.

Objective

The irritative zone - the area generating epileptic spikes - can be studied non-invasively during the interictal period using Electrical Source Imaging (ESI) and simultaneous electroencephalography-functional magnetic resonance imaging (EEG-fMRI). Although the techniques yield results which may overlap spatially, differences in spatial localization of the irritative zone within the same patient are consistently observed. To investigate this discrepancy, we used Blood Oxygenation Level Dependent (BOLD) functional connectivity measures to examine the underlying relationship between ESI and EEG-fMRI findings.

Methods

Fifteen patients (age 20–54), who underwent presurgical epilepsy investigation, were scanned using a single-session resting-state EEG-fMRI protocol. Structural MRI was used to obtain the electrode localisation of a high-density 64-channel EEG cap. Electrical generators of interictal epileptiform discharges were obtained using a distributed local autoregressive average (LAURA) algorithm as implemented in Cartool EEG software. BOLD activations were obtained using both spike-related and voltage-map EEG-fMRI analysis. The global maxima of each method were used to investigate the temporal relationship of BOLD time courses and to assess the spatial similarity using the Dice similarity index between functional connectivity maps.

Results

ESI, voltage-map and spike-related EEG-fMRI methods identified peaks in 15 (100%), 13 (67%) and 8 (53%) of the 15 patients, respectively. For all methods, maxima were localised within the same lobe, but differed in sub-lobar localisation, with a median distance of 22.8 mm between the highest peak for each method. The functional connectivity analysis showed that the temporal correlation between maxima only explained 38% of the variance between the time course of the BOLD response at the maxima. The mean Dice similarity index between seed-voxel functional connectivity maps showed poor spatial agreement.

Significance

Non-invasive methods for the localisation of the irritative zone have distinct spatial and temporal sensitivity to different aspects of the local cortical network involved in the generation of interictal epileptiform discharges.

Pushing the Limits of EEG: Estimation of Large-Scale Functional Brain Networks and Their Dynamics Validated by Simultaneous fMRI

Functional magnetic resonance imaging (fMRI) is the technique of choice for detecting large-scale functional brain networks and to investigate their dynamics. Because fMRI measures brain activity indirectly, electroencephalography (EEG) has been recently considered a feasible tool for detecting such networks, particularly the resting-state networks (RSNs). However, a truly unbiased validation of such claims is still missing, which can only be accomplished by using simultaneously acquired EEG and fMRI data, due to the spontaneous nature of the activity underlying the RSNs. Additionally, EEG is still poorly explored for the purpose of mapping task-specific networks, and no studies so far have been focused on investigating networks' dynamic functional connectivity (dFC) with EEG. Here, we started by validating RSNs derived from the continuous reconstruction of EEG sources by directly comparing them with those derived from simultaneous fMRI data of 10 healthy participants, and obtaining an average overlap (quantified by the Dice coefficient) of 0.4. We also showed the ability of EEG to map the facial expressions processing network (FEPN), highlighting regions near the posterior superior temporal sulcus, where the FEPN is anchored. Then, we measured the dFC using EEG for the first time in this context, estimated dFC brain states using dictionary learning, and compared such states with those obtained from the fMRI. We found a statistically significant match between fMRI and EEG dFC states, and determined the existence of two matched dFC states which contribution over time was associated with the brain activity at the FEPN, showing that the dynamics of FEPN can be captured by both fMRI and EEG. Our results push the limits of EEG toward being used as a brain imaging tool, while supporting the growing literature on EEG correlates of (dynamic) functional connectivity measured with fMRI, and providing novel insights into the coupling mechanisms underlying the two imaging techniques.

WONOEP appraisal: Network concept from an imaging perspective

Neuroimaging techniques applied to a variety of organisms-from zebrafish, to rodents to humans-can offer valuable insights into neuronal network properties and their dysfunction in epilepsy. A wide range of imaging methods used to monitor neuronal circuits and networks during evoked seizures in animal models and advances in functional magnetic resonance imaging (fMRI) applied to patients with epilepsy were discussed during the XIV Workshop on Neurobiology of Epilepsy (XIV WONOEP) organized in 2017 by the Neurobiology Commission of the International League Against Epilepsy (ILAE). We review the growing number of technological approaches developed, as well as the current state of knowledge gained from studies applying these advanced imaging approaches to epilepsy research.

Quantitative determination of concordance in localizing epileptic focus by component-based EEG-fMRI

BACKGROUND AND OBJECTIVE

Accurate seizure onset zone (SOZ) localization is an essential step in pre-surgical assessment of patients with refractory focal epilepsy. Complex pathophysiology of epileptic cerebral structures, seizure types and frequencies have not been considered as influential features for accurate identification of SOZ using EEG-fMRI. There is a crucial need to quantitatively measure concordance between presumed SOZ and IED-related BOLD response in different brain regions to improve SOZ delineation.

METHODS

A novel component-based EEG-fMRI approach is proposed to measure physical distance between BOLD clusters and selected component dipole location using patient-specific high resolution anatomical images. The method is applied on 18 patients with refractory focal epilepsy to localize epileptic focus and determine concordance quantitatively and compare between maximum BOLD cluster with identified component dipole. To measure concordance, distance from a voxel with maximal z-score of maximum BOLD to center of extracted component dipole is measured.

RESULTS

BOLD clusters to spikes distances for concordant (<25 mm), partially concordant (25-50 mm), and discordant (>50 mm) groups were significantly different (p < 0.0001). The results showed full concordance in 17 IED types (17.85 ± 4.69 mm), partial concordance in 4 (36.47 ± 8.84 mm), and nodiscordance, which is a significant rise compared to the existing literature. The proposed method is premised on the cross-correlation between the spike template outside the scanner and the highly-ranked extracted components. It successfully surpasses the limitations of conventional EEG-fMRI studies which are largely dependent on inside-scanner spikes. More significantly, the proposed method improves localization accuracy to 97% which marks a dramatic rise compared to conventional works.

CONCLUSIONS

This study demonstrated that BOLD changes were related to epileptic spikes in different brain regions in patients with refractory focal epilepsy. In a systematic quantitative approach, concordance levels based on the distance between center of maximum BOLD cluster and dipole were determined by component-based EEG-fMRI method. Therefore, component-based EEG-fMRI can be considered as a reliable predictor of SOZ in patients with focal epilepsy and included as part of clinical evaluation for patients with medically resistant epilepsy.

The impact of mapping interictal discharges using EEG-fMRI on the epilepsy presurgical clinical decision making process: A prospective study

PURPOSE

We set out to establish the clinical utility of EEG-correlated fMRI as part of the presurgical evaluation, by measuring prospectively its effects on the clinical decision.

METHODS

Patients with refractory extra-temporal focal epilepsy, referred for presurgical evaluation were recruited in a period of 18 months. The EEG-fMRI based localization was presented during a multi-disciplinary meeting after the team had defined the presumed RESULTS: Sixteen patients (six women), with a median age of 28 years, were recruited. Interpretable EEG-fMRI results were available in 13: interictal epileptic discharges (IEDs) were recorded in eleven patients and seizures were recorded in two patients. In three patients, no epileptic activity was captured during EEG-fMRI acquisition and in two of those an IED topographic map correlation was performed (between EEG recorded inside the scanner and long-term video EEG monitoring). EEG-fMRI results presentation had no impact on the initial clinical decision in three patients (23%) of the thirteen and resulted in a modification of the initial surgical plan in ten patients (77%) of the thirteen finally presented in MDT; in eight patients the impact was on the planned placement of invasive electrodes and in two patients the EEG-fMRI led to additional non-invasive tests before proceeding further with surgery.

CONCLUSION

The study is a prospective observational cohort study specifically designed to assess the impact of EEG-fMRI on the clinical decision making process, suggesting a significant influence of EEG-fMRI on epilepsy surgery planning.

EEG synchronization measures predict epilepsy-related BOLD-fMRI fluctuations better than commonly used univariate metrics

OBJECTIVE

We hypothesize that the hypersynchronization associated with epileptic activity is best described by EEG synchronization measures, and propose to use these as predictors of epilepsy-related BOLD fluctuations.

METHODS

We computed the phase synchronization index (PSI) and global field synchronization (GFS), within two frequency bands, a broadband (1-45 Hz) and a narrower band focused on the presence of epileptic activity (3-10 Hz). The associated epileptic networks were compared with those obtained using conventional unitary regressors and two power-weighted metrics (total power and root mean square frequency), on nine simultaneous EEG-fMRI datasets from four epilepsy patients, exhibiting inter-ictal epileptiform discharges (IEDs).

RESULTS

The average PSI within 3-10 Hz achieved the best performance across several measures reflecting reliability in all datasets. The results were cross-validated through electrical source imaging of the IEDs. The applicability of PSI when no IEDs are recorded on the EEG was evaluated on three additional patients, yielding partially plausible networks in all cases.

CONCLUSIONS

Epileptic networks can be mapped based on the EEG PSI metric within an IED-specific frequency band, performing better than commonly used EEG metrics.

SIGNIFICANCE

This is the first study to investigate EEG synchronization measures as potential predictors of epilepsy-related BOLD fluctuations.

Combined electroencephalography-functional magnetic resonance imaging and electrical source imaging improves localization of pediatric focal epilepsy

OBJECTIVE

Surgical treatment in epilepsy is effective if the epileptogenic zone (EZ) can be correctly localized and characterized. Here we use simultaneous electroencephalography-functional magnetic resonance imaging (EEG-fMRI) data to derive EEG-fMRI and electrical source imaging (ESI) maps. Their yield and their individual and combined ability to (1) localize the EZ and (2) predict seizure outcome were then evaluated.

METHODS

Fifty-three children with drug-resistant epilepsy underwent EEG-fMRI. Interictal discharges were mapped using both EEG-fMRI hemodynamic responses and ESI. A single localization was derived from each individual test (EEG-fMRI global maxima [GM]/ESI maximum) and from the combination of both maps (EEG-fMRI/ESI spatial intersection). To determine the localization accuracy and its predictive performance, the individual and combined test localizations were compared to the presumed EZ and to the postsurgical outcome.

RESULTS

Fifty-two of 53 patients had significant maps: 47 of 53 for EEG-fMRI, 44 of 53 for ESI, and 34 of 53 for both. The EZ was well characterized in 29 patients; 26 had an EEG-fMRI GM localization that was correct in 11, 22 patients had ESI localization that was correct in 17, and 12 patients had combined EEG-fMRI and ESI that was correct in 11. Seizure outcome following resection was correctly predicted by EEG-fMRI GM in 8 of 20 patients, and by the ESI maximum in 13 of 16. The combined EEG-fMRI/ESI region entirely predicted outcome in 9 of 9 patients, including 3 with no lesion visible on MRI.

INTERPRETATION

EEG-fMRI combined with ESI provides a simple unbiased localization that may predict surgery better than each individual test, including in MRI-negative patients. Ann Neurol 2017;82:278-287.

Contributions of EEG-fMRI to Assessing the Epileptogenicity of Focal Cortical Dysplasia

Purpose: To examine the ability of the BOLD response to EEG spikes to assess the epileptogenicity of the lesion in patients with focal cortical dysplasia (FCD).

Method: Patients with focal epilepsy and FCD who underwent 3T EEG-fMRI from 2006 to 2010 were included. Diagnosis of FCD was based on neuroradiology (MRI+), or histopathology in MRI-negative cases (MRI−). Patients underwent 120 min EEG-fMRI recording session. Spikes similar to those recorded outside the scanner were marked in the filtered EEG. The lesion (in MRI+) or the removed cortex (in MRI−) was marked on the anatomical T1 sequence, blindly to the BOLD response, after reviewing the FLAIR images. For each BOLD response we assessed the concordance with the spike field and with the lesion in MRI+ or the removed cortex in MRI−. BOLD responses were considered “concordant” if the maximal t-value was inside the marking. Follow-up after resection was used as gold-standard.

Results: Twenty patients were included (13 MRI+, 7 MRI−), but in seven the EEG was not active or there were artifacts during acquisition. In all 13 studied patients, at least one BOLD response was concordant with the spike field; in 9/13 (69%) at least one BOLD response was concordant with the lesion: in 6/7 (86%) MRI+ and in 3/6 (50%) MRI− patients.

Conclusions: Our study shows a high level of concordance between FCD and BOLD response. This data could provide useful information especially for MRI negative patients. Moreover, it shows in almost all FCD patients, a metabolic involvement of remote cortical or subcortical structures, corroborating the concept of epileptic network.

EEG-fMRI in the presurgical evaluation of temporal lobe epilepsy

OBJECTIVE

Drug-resistant temporal lobe epilepsy (TLE) often requires thorough investigation to define the epileptogenic zone for surgical treatment. We used simultaneous interictal scalp EEG-fMRI to evaluate its value for predicting long-term postsurgical outcome.

METHODS

30 patients undergoing presurgical evaluation and proceeding to temporal lobe (TL) resection were studied. Interictal epileptiform discharges (IEDs) were identified on intra-MRI EEG and used to build a model of haemodynamic changes. In addition, topographic electroencephalographic correlation maps were calculated between the average IED during video-EEG and intra-MRI EEG, and used as a condition. This allowed the analysis of all data irrespective of the presence of IED on intra-MRI EEG. Mean follow-up after surgery was 46 months. International League Against Epilepsy (ILAE) outcomes 1 and 2 were considered good, and 3-6 poor, surgical outcome. Haemodynamic maps were classified according to the presence (Concordant) or absence (Discordant) of Blood Oxygen Level-Dependent (BOLD) change in the TL overlapping with the surgical resection.

RESULTS

The proportion of patients with good surgical outcome was significantly higher (13/16; 81%) in the Concordant than in the Discordant group (3/14; 21%) (χ(2) test, Yates correction, p=0.003) and multivariate analysis showed that Concordant BOLD maps were independently related to good surgical outcome (p=0.007). Sensitivity and specificity of EEG-fMRI results to identify patients with good surgical outcome were 81% and 79%, respectively, and positive and negative predictive values were 81% and 79%, respectively.

INTERPRETATION

The presence of significant BOLD changes in the area of resection on interictal EEG-fMRI in patients with TLE retrospectively confirmed the epileptogenic zone. Surgical resection including regions of haemodynamic changes in the TL may lead to better postoperative outcome.

Towards motion insensitive EEG-fMRI: Correcting motion-induced voltages and gradient artefact instability in EEG using an fMRI prospective motion correction (PMC) system

The simultaneous acquisition of electroencephalography and functional magnetic resonance imaging (EEG-fMRI) is a multimodal technique extensively applied for mapping the human brain. However, the quality of EEG data obtained within the MRI environment is strongly affected by subject motion due to the induction of voltages in addition to artefacts caused by the scanning gradients and the heartbeat. This has limited its application in populations such as paediatric patients or to study epileptic seizure onset. Recent work has used a Moiré-phase grating and a MR-compatible camera to prospectively update image acquisition and improve fMRI quality (prospective motion correction: PMC). In this study, we use this technology to retrospectively reduce the spurious voltages induced by motion in the EEG data acquired inside the MRI scanner, with and without fMRI acquisitions. This was achieved by modelling induced voltages from the tracking system motion parameters; position and angles, their first derivative (velocities) and the velocity squared. This model was used to remove the voltages related to the detected motion via a linear regression. Since EEG quality during fMRI relies on a temporally stable gradient artefact (GA) template (calculated from averaging EEG epochs matched to scan volume or slice acquisition), this was evaluated in sessions both with and without motion contamination, and with and without PMC. We demonstrate that our approach is capable of significantly reducing motion-related artefact with a magnitude of up to 10mm of translation, 6° of rotation and velocities of 50mm/s, while preserving physiological information. We also demonstrate that the EEG-GA variance is not increased by the gradient direction changes associated with PMC. Provided a scan slice-based GA template is used (rather than a scan volume GA template) we demonstrate that EEG variance during motion can be supressed towards levels found when subjects are still. In summary, we show that PMC can be used to dramatically improve EEG quality during large amplitude movements, while benefiting from previously reported improvements in fMRI quality, and does not affect EEG data quality in the absence of large amplitude movements.

Presurgical brain mapping in epilepsy using simultaneous EEG and functional MRI at ultra-high field: feasibility and first results

OBJECTIVES

The aim of this study was to demonstrate that eloquent cortex and epileptic-related hemodynamic changes can be safely and reliably detected using simultaneous electroencephalography (EEG)-functional magnetic resonance imaging (fMRI) recordings at ultra-high field (UHF) for clinical evaluation of patients with epilepsy.

MATERIALS AND METHODS

Simultaneous EEG-fMRI was acquired at 7 T using an optimized setup in nine patients with lesional epilepsy. According to the localization of the lesion, mapping of eloquent cortex (language and motor) was also performed in two patients.

RESULTS

Despite strong artifacts, efficient correction of intra-MRI EEG could be achieved with optimized artifact removal algorithms, allowing robust identification of interictal epileptiform discharges. Noise-sensitive topography-related analyses and electrical source localization were also performed successfully. Localization of epilepsy-related hemodynamic changes compatible with the lesion were detected in three patients and concordant with findings obtained at 3 T. Local loss of signal in specific regions, essentially due to B 1 inhomogeneities were found to depend on the geometric arrangement of EEG leads over the cap.

CONCLUSION

These results demonstrate that presurgical mapping of epileptic networks and eloquent cortex is both safe and feasible at UHF, with the benefits of greater spatial resolution and higher blood-oxygenation-level-dependent sensitivity compared with the more traditional field strength of 3 T.

Incorporating priors for EEG source imaging and connectivity analysis

Electroencephalography source imaging (ESI) is a useful technique to localize the generators from a given scalp electric measurement and to investigate the temporal dynamics of the large-scale neural circuits. By introducing reasonable priors from other modalities, ESI reveals the most probable sources and communication structures at every moment in time. Here, we review the available priors from such techniques as magnetic resonance imaging (MRI), functional MRI (fMRI), and positron emission tomography (PET). The modality's specific contribution is analyzed from the perspective of source reconstruction. For spatial priors, EEG-correlated fMRI, temporally coherent networks (TCNs) and resting-state fMRI are systematically introduced in the ESI. Moreover, the fiber tracking (diffusion tensor imaging, DTI) and neuro-stimulation techniques (transcranial magnetic stimulation, TMS) are also introduced as the potential priors, which can help to draw inferences about the neuroelectric connectivity in the source space. We conclude that combining EEG source imaging with other complementary modalities is a promising approach toward the study of brain networks in cognitive and clinical neurosciences.

Distributions of Irritative Zones Are Related to Individual Alterations of Resting-State Networks in Focal Epilepsy

Alterations in the connectivity patterns of the fMRI-based resting-state networks (RSNs) have been reported in several types of epilepsies. Evidence pointed out these alterations might be associated with the genesis and propagation of interictal epileptiform discharges (IEDs). IEDs also evoke blood-oxygen-level dependent (BOLD) responses, which have been used to delineate irritative zones during preoperative work-up. Therefore, one may expect a relationship between the topology of the IED-evoked BOLD response network and the altered spatial patterns of the RSNs. In this study, we used EEG recordings and fMRI data obtained simultaneously from a chronic model of focal epilepsy in Wistar rats to verify our hypothesis. We found that IED-evoked BOLD response networks comprise both cortical and subcortical structures with a rat-dependent topology. In all rats, IEDs evoke both activation and deactivation types of BOLD responses. Using a Granger causality method, we found that in many cases areas with BOLD deactivation have directed influences on areas with activation (p<0.05). We were able to predict topological properties (i.e., focal/diffused, unilateral/bilateral) of the IED-evoked BOLD response network by performing hierarchical clustering analysis on major spatial features of the RSNs. All these results suggest that IEDs and disruptions in the RSNs found previously in humans may be different manifestations of the same transient events, probably reflecting altered consciousness. In our opinion, the shutdown of specific nodes of the default mode network may cause uncontrollable excitability in other functionally connected brain areas. We conclude that IED-evoked BOLD responses (i.e., activation and deactivation) and alterations of RSNs are intrinsically related, and speculate that an understanding of their interplay is necessary to discriminate focal epileptogenesis and network propagation phenomena across different brain modules via hub-based connectivity.

Intrinsic coupling modes in source-reconstructed electroencephalography

Intrinsic coupling of neuronal assemblies constitutes a key feature of ongoing brain activity, yielding the rich spatiotemporal patterns observed in neuroimaging data and putatively supporting cognitive processes. Intrinsic coupling has been investigated in electrophysiological recordings using two types of functional connectivity measures: amplitude and phase coupling. These two coupling modes differ in their likely causes and functions, and have been proposed to provide complementary insights into intrinsic neuronal interactions. Here, we investigate the relationship between amplitude and phase coupling in source-reconstructed electroencephalography (EEG). Volume conduction is a key obstacle for connectivity analysis in EEG-we therefore also test the envelope correlation of orthogonalized signals and the phase lag index. Functional connectivity between six seed source regions (bilateral visual, sensorimotor, and auditory cortices) and all other cortical voxels was computed. For all four measures, coupling between homologous sensory areas in both hemispheres was significantly higher than with other voxels at the same physical distance. The frequency of significant coupling differed between sensory areas: 10 Hz for visual, 30 Hz for auditory, and 40 Hz for sensorimotor cortices. By contrasting envelope correlations and phase locking values, we observed two distinct clusters of voxels showing a different relationship between amplitude and phase coupling. Large clusters contiguous to the seed regions showed an identity (1:1) relationship between amplitude and phase coupling, whereas a cluster located around the contralateral homologous regions showed higher phase than amplitude coupling. These results show a relationship between intrinsic coupling modes that is distinct from the effect of volume conduction.

Advantages and disadvantages of a fast fMRI sequence in the context of EEG-fMRI investigation of epilepsy patients: A realistic simulation study

EEG-fMRI is an established technique to allow mapping BOLD changes in response to interictal discharges recorded in the EEG of epilepsy patients. Traditional fMRI experiments rely on an echo planar imaging (EPI) sequence with a time resolution given by its time-to-repetition (TR) of ~2 s. Recently, multiple fast fMRI sequences have been developed to get around the limited temporal resolution of the EPI sequence, and achieved a TR in the 100 ms range or lower. One such sequence is called Magnetic Resonance EncephaloGraphy (MREG). Its high temporal resolution should offer increased detection sensitivity and statistical power in the context of epilepsy studies and in fMRI experiments in general. The aim of this work was to investigate the advantages and disadvantages offered by MREG. This was done by superimposing artificial event-related BOLD responses on EPI and MREG background signals, from 5 epileptic patients, that were free of epileptic discharges (spikes) on simultaneously recorded EEG. These functional datasets simulated different spiking rates and hemodynamic response amplitudes, and were analyzed with the commonly used General Linear Model (GLM) with the canonical hemodynamic response function (HRF) as a fixed model of the response. Robustness to violation of the assumptions of the GLM was additionally assessed with similar simulations using variable spike-to-spike response amplitudes and 8 non-canonical HRFs. Consistent with previous work, MREG yields higher maximum statistical t-values than EPI, but our simulations showed these statistics to be inflated, as the false positive rate at a standard threshold was high. At thresholds set to appropriately control specificity, EPI showed better true positive rate and larger cluster size than MREG. However, the lack of an appropriate calibration of the amplitude of the responses across the sequences precludes definitive judgment on their relative sensitivity. In addition, we show that a mismatch between the assumed and actual HRF impairs more MREG detection performance, but that EPI is more affected by non-modeled spike-to-spike variations of response amplitude. Filtering-out physiological noise, which is not aliased at the fast sampling rate of MREG, and the modeling of temporal autocorrelation are advantageous in increasing the detection power of MREG. This simulation study 1) warrants care when interpreting statistical t-values from fast fMRI sequences, 2) proposes thresholds for valid inferences and processing methods for maximal sensitivities, and 3) demonstrates the relative robustness/susceptibility of MREG and EPI to violation of the GLM's assumptions.

A Comparison of Independent Component Analysis (ICA) of fMRI and Electrical Source Imaging (ESI) in Focal Epilepsy Reveals Misclassification Using a Classifier

Interictal epileptiform discharges (IEDs) can produce haemodynamic responses that can be detected by electroencephalography-functional magnetic resonance imaging (EEG-fMRI) using different analysis methods such as the general linear model (GLM) of IEDs or independent component analysis (ICA). The IEDs can also be mapped by electrical source imaging (ESI) which has been demonstrated to be useful in presurgical evaluation in a high proportion of cases with focal IEDs. ICA advantageously does not require IEDs or a model of haemodynamic responses but its use in EEG-fMRI of epilepsy has been limited by its ability to separate and select epileptic components. Here, we evaluated the performance of a classifier that aims to filter all non-BOLD responses and we compared the spatial and temporal features of the selected independent components (ICs). The components selected by the classifier were compared to those components selected by a strong spatial correlation with ESI maps of IED sources. Both sets of ICs were subsequently compared to a temporal model derived from the convolution of the IEDs (derived from the simultaneously acquired EEG) with a standard haemodynamic response. Selected ICs were compared to the patients' clinical information in 13 patients with focal epilepsy. We found that the misclassified ICs clearly related to IED in 16/25 cases. We also found that the classifier failed predominantly due to the increased spectral range of fMRIs temporal responses to IEDs. In conclusion, we show that ICA can be an efficient approach to separate responses related to epilepsy but that contemporary classifiers need to be retrained for epilepsy data. Our findings indicate that, for ICA to contribute to the analysis of data without IEDs to improve its sensitivity, classification strategies based on data features other than IC time course frequency is required.

Patient-specific detection of cerebral blood flow alterations as assessed by arterial spin labeling in drug-resistant epileptic patients

Electrophysiological and hemodynamic data can be integrated to accurately and precisely identify the generators of abnormal electrical activity in drug-resistant focal epilepsy. Arterial Spin Labeling (ASL), a magnetic resonance imaging (MRI) technique for quantitative noninvasive measurement of cerebral blood flow (CBF), can provide a direct measure of variations in cerebral perfusion associated with the epileptic focus. In this study, we aimed to confirm the ASL diagnostic value in the identification of the epileptogenic zone, as compared to electrical source imaging (ESI) results, and to apply a template-based approach to depict statistically significant CBF alterations. Standard video-electroencephalography (EEG), high-density EEG, and ASL were performed to identify clinical seizure semiology and noninvasively localize the epileptic focus in 12 drug-resistant focal epilepsy patients. The same ASL protocol was applied to a control group of 17 healthy volunteers from which a normal perfusion template was constructed using a mixed-effect approach. CBF maps of each patient were then statistically compared to the reference template to identify perfusion alterations. Significant hypo- and hyperperfused areas were identified in all cases, showing good agreement between ASL and ESI results. Interictal hypoperfusion was observed at the site of the seizure in 10/12 patients and early postictal hyperperfusion in 2/12. The epileptic focus was correctly identified within the surgical resection margins in the 5 patients who underwent lobectomy, all of which had good postsurgical outcomes. The combined use of ESI and ASL can aid in the noninvasive evaluation of drug-resistant epileptic patients.

Neuronal Networks during Burst Suppression as Revealed by Source Analysis

Introduction

Burst-suppression (BS) is an electroencephalography (EEG) pattern consisting of alternant periods of slow waves of high amplitude (burst) and periods of so called flat EEG (suppression). It is generally associated with coma of various etiologies (hypoxia, drug-related intoxication, hypothermia, and childhood encephalopathies, but also anesthesia). Animal studies suggest that both the cortex and the thalamus are involved in the generation of BS. However, very little is known about mechanisms of BS in humans. The aim of this study was to identify the neuronal network underlying both burst and suppression phases using source reconstruction and analysis of functional and effective connectivity in EEG.

Material/Methods

Dynamic imaging of coherent sources (DICS) was applied to EEG segments of 13 neonates and infants with burst and suppression EEG pattern. The brain area with the strongest power in the analyzed frequency (1–4 Hz) range was defined as the reference region. DICS was used to compute the coherence between this reference region and the entire brain. The renormalized partial directed coherence (RPDC) was used to describe the informational flow between the identified sources.

Results/Conclusion

Delta activity during the burst phases was associated with coherent sources in the thalamus and brainstem as well as bilateral sources in cortical regions mainly frontal and parietal, whereas suppression phases were associated with coherent sources only in cortical regions. Results of the RPDC analyses showed an upwards informational flow from the brainstem towards the thalamus and from the thalamus to cortical regions, which was absent during the suppression phases. These findings may support the theory that a “cortical deafferentiation” between the cortex and sub-cortical structures exists especially in suppression phases compared to burst phases in burst suppression EEGs. Such a deafferentiation may play a role in the poor neurological outcome of children with these encephalopathies.

Multimodal effective connectivity analysis reveals seizure focus and propagation in musicogenic epilepsy

Dynamic causal modeling (DCM) is a method to non-invasively assess effective connectivity between brain regions. 'Musicogenic epilepsy' is a rare reflex epilepsy syndrome in which seizures can be elicited by musical stimuli and thus represents a unique possibility to investigate complex human brain networks and test connectivity analysis tools. We investigated effective connectivity in a case of musicogenic epilepsy using DCM for fMRI, high-density (hd-) EEG and MEG and validated results with intracranial EEG recordings. A patient with musicogenic seizures was examined using hd-EEG/fMRI and simultaneous '256-channel hd-EEG'/'whole head MEG' to characterize the epileptogenic focus and propagation effects using source analysis techniques and DCM. Results were validated with invasive EEG recordings. We recorded one seizure with hd-EEG/fMRI and four auras with hd-EEG/MEG. During the seizures, increases of activity could be observed in the right mesial temporal region as well as bilateral mesial frontal regions. Effective connectivity analysis of fMRI and hd-EEG/MEG indicated that right mesial temporal neuronal activity drives changes in the frontal areas consistently in all three modalities, which was confirmed by the results of invasive EEG recordings. Seizures thus seem to originate in the right mesial temporal lobe and propagate to mesial frontal regions. Using DCM for fMRI, hd-EEG and MEG we were able to correctly localize focus and propagation of epileptic activity and thereby characterize the underlying epileptic network in a patient with musicogenic epilepsy. The concordance between all three functional modalities validated by invasive monitoring is noteworthy, both for epileptic activity spread as well as for effective connectivity analysis in general.

Mapping epileptic activity: sources or networks for the clinicians?

Epileptic seizures of focal origin are classically considered to arise from a focal epileptogenic zone and then spread to other brain regions. This is a key concept for semiological electro-clinical correlations, localization of relevant structural lesions, and selection of patients for epilepsy surgery. Recent development in neuro-imaging and electro-physiology and combinations, thereof, have been validated as contributory tools for focus localization. In parallel, these techniques have revealed that widespread networks of brain regions, rather than a single epileptogenic region, are implicated in focal epileptic activity. Sophisticated multimodal imaging and analysis strategies of brain connectivity patterns have been developed to characterize the spatio-temporal relationships within these networks by combining the strength of both techniques to optimize spatial and temporal resolution with whole-brain coverage and directional connectivity. In this paper, we review the potential clinical contribution of these functional mapping techniques as well as invasive electrophysiology in human beings and animal models for characterizing network connectivity.

Survey of encoding and decoding of visual stimulus via FMRI: an image analysis perspective

A variety of exciting scientific achievements have been made in the last few decades in brain encoding and decoding via functional magnetic resonance imaging (fMRI). This trend continues to rise in recent years, as evidenced by the increasing number of published papers in this topic and several published survey papers addressing different aspects of research issues. Essentially, these survey articles were mainly from cognitive neuroscience and neuroimaging perspectives, although computational challenges were briefly discussed. To complement existing survey articles, this paper focuses on the survey of the variety of image analysis methodologies, such as neuroimage registration, fMRI signal analysis, ROI (regions of interest) selection, machine learning algorithms, reproducibility analysis, structural and functional connectivity, and natural image analysis, which were employed in previous brain encoding/decoding research works. This paper also provides discussions of potential limitations of those image analysis methodologies and possible future improvements. It is hoped that extensive discussions of image analysis issues could contribute to the advancements of the increasingly important brain encoding/decoding field.

Pulse artifact detection in simultaneous EEG-fMRI recording based on EEG map topography

One of the major artifact corrupting electroencephalogram (EEG) acquired during functional magnetic resonance imaging (fMRI) is the pulse artifact (PA). It is mainly due to the motion of the head and attached electrodes and wires in the magnetic field occurring after each heartbeat. In this study we propose a novel method to improve PA detection by considering the strong gradient and inversed polarity between left and right EEG electrodes. We acquired high-density EEG-fMRI (256 electrodes) with simultaneous electrocardiogram (ECG) at 3 T. PA was estimated as the voltage difference between right and left signals from the electrodes showing the strongest artifact (facial and temporal). Peaks were detected on this estimated signal and compared to the peaks in the ECG recording. We analyzed data from eleven healthy subjects, two epileptic patients and four healthy subjects with an insulating layer between electrodes and scalp. The accuracy of the two methods was assessed with three criteria: (i) standard deviation, (ii) kurtosis and (iii) confinement into the physiological range of the inter-peak intervals. We also checked whether the new method has an influence on the identification of epileptic spikes. Results show that estimated PA improved artifact detection in 15/17 cases, when compared to the ECG method. Moreover, epileptic spike identification was not altered by the correction. The proposed method improves the detection of pulse-related artifacts, particularly crucial when the ECG is of poor quality or cannot be recorded. It will contribute to enhance the quality of the EEG increasing the reliability of EEG-informed fMRI analysis.

Sensitivity and Specificity of Interictal EEG-fMRI for Detecting the Ictal Onset Zone at Different Statistical Thresholds

There is currently a lack of knowledge about electroencephalography (EEG)-functional magnetic resonance imaging (fMRI) specificity. Our aim was to define sensitivity and specificity of blood oxygen level dependent (BOLD) responses to interictal epileptic spikes during EEG-fMRI for detecting the ictal onset zone (IOZ). We studied 21 refractory focal epilepsy patients who had a well-defined IOZ after a full presurgical evaluation and interictal spikes during EEG-fMRI. Areas of spike-related BOLD changes overlapping the IOZ in patients were considered as true positives; if no overlap was found, they were treated as false-negatives. Matched healthy case-controls had undergone similar EEG-fMRI in order to determine true-negative and false-positive fractions. The spike-related regressor of the patient was used in the design matrix of the healthy case-control. Suprathreshold BOLD changes in the brain of controls were considered as false positives, absence of these changes as true negatives. Sensitivity and specificity were calculated for different statistical thresholds at the voxel level combined with different cluster size thresholds and represented in receiver operating characteristic (ROC)-curves. Additionally, we calculated the ROC-curves based on the cluster containing the maximal significant activation. We achieved a combination of 100% specificity and 62% sensitivity, using a Z-threshold in the interval 3.4–3.5 and cluster size threshold of 350 voxels. We could obtain higher sensitivity at the expense of specificity. Similar performance was found when using the cluster containing the maximal significant activation. Our data provide a guideline for different EEG-fMRI settings with their respective sensitivity and specificity for detecting the IOZ. The unique cluster containing the maximal significant BOLD activation was a sensitive and specific marker of the IOZ.

EEG-microstate dependent emergence of perceptual awareness

We investigated whether the differences in perceptual awareness for stimuli at the threshold of awareness can arise from different global brain states before stimulus onset indexed by the EEG microstate. We used a metacontrast backward masking paradigm in which subjects had to discriminate between two weak stimuli and obtained measures of accuracy and awareness while their EEG was recorded from 256 channels. Comparing targets that were correctly identified with and without awareness allowed us to contrast differences in awareness while keeping performance constant for identical physical stimuli. Two distinct pre-stimulus scalp potential fields (microstate maps) dissociated correct identification with and without awareness, and their estimated intracranial generators were stronger in primary visual cortex before correct identification without awareness. This difference in activity cannot be explained by differences in alpha power or phase which were less reliably linked with differential pre-stimulus activation of primary visual cortex. Our results shed a new light on the function of pre-stimulus activity in early visual cortex in visual awareness and emphasize the importance of trial-by-trials analysis of the spatial configuration of the scalp potential field identified with multichannel EEG.

EEG/fMRI contributions to our understanding of genetic generalized epilepsies

The first reports of combined EEG and fMRI used for evaluation of epileptic spikes date back to the mid-90s. At that time, the technique was called EEG-triggered fMRI--the "triggered" corresponded to an epilepsy specialist reviewing live EEG while the patient was located in the scanner; after the spike was identified, a scan was initiated to collect the data. Since then major progress has been made in combined EEG/fMRI data collection and analyses. These advances allow studying the electrophysiology of genetic generalized epilepsies (GGEs) in vivo in greater detail than ever. In addition to continuous data collection, we now have better methods for removing physiologic and fMRI-related artifacts, more advanced understanding of the hemodynamic response functions, and better computational methods to address the questions regarding the origins of the epileptiform discharge generators in patients with GGEs. These advances have allowed us to examine numerous cohorts of children and adults with GGEs while not only looking for spike and wave generators but also examining specific types of GGEs (e.g., juvenile myoclonic epilepsy or childhood absence epilepsy), drug-naïve patients, effects of medication resistance, or effects of epileptiform abnormalities and/or seizures on brain connectivity. While the discussion is ongoing, the prevailing thought is that the GGEs as a group are a network disorder with participation from multiple nodes including the thalami and cortex with the clinical presentation depending on which node of the participating network is affected by the disease process. This review discusses the contributions of EEG/fMRI to our understanding of GGEs.

The role of functional neuroimaging in pre-surgical epilepsy evaluation

The prevalence of epilepsy is about 1% and one-third of cases do not respond to medical treatment. In an eligible subset of patients with drug-resistant epilepsy, surgical resection of the epileptogenic zone is the only treatment that can possibly cure the disease. Non-invasive techniques provide information for the localization of the epileptic focus in the majority of cases, whereas in others invasive procedures are required. In the last years, non-invasive neuroimaging techniques, such as simultaneous recording of functional magnetic resonance imaging and electroencephalogram (EEG-fMRI), positron emission tomography (PET), single photon emission computed tomography (SPECT), electric and magnetic source imaging (MSI, ESI), spectroscopy (MRS), have proved their usefulness in defining the epileptic focus. The combination of these functional techniques can yield complementary information and their concordance is crucial for guiding clinical decision, namely the planning of invasive EEG recordings or respective surgery. The aim of this review is to present these non-invasive neuroimaging techniques, their potential combination, and their role in the pre-surgical evaluation of patients with pharmaco-resistant epilepsy.

The role of functional neuroimaging in pre-surgical epilepsy evaluation

The prevalence of epilepsy is about 1% and one-third of cases do not respond to medical treatment. In an eligible subset of patients with drug-resistant epilepsy, surgical resection of the epileptogenic zone is the only treatment that can possibly cure the disease. Non-invasive techniques provide information for the localization of the epileptic focus in the majority of cases, whereas in others invasive procedures are required. In the last years, non-invasive neuroimaging techniques, such as simultaneous recording of functional magnetic resonance imaging and electroencephalogram (EEG-fMRI), positron emission tomography (PET), single photon emission computed tomography (SPECT), electric and magnetic source imaging (MSI, ESI), spectroscopy (MRS), have proved their usefulness in defining the epileptic focus. The combination of these functional techniques can yield complementary information and their concordance is crucial for guiding clinical decision, namely the planning of invasive EEG recordings or respective surgery. The aim of this review is to present these non-invasive neuroimaging techniques, their potential combination, and their role in the pre-surgical evaluation of patients with pharmaco-resistant epilepsy.

Spatial correlation of hemodynamic changes related to interictal epileptic discharges with electric and magnetic source imaging

INTRODUCTION

Blood oxygenation level-dependent (BOLD) signal changes at the time of interictal epileptic discharges (IEDs) identify their associated vascular/hemodynamic responses. BOLD activations and deactivations can be found within the epileptogenic zone but also at a distance. Source imaging identifies electric (ESI) and magnetic (MSI) sources of IEDs, with the advantage of a higher temporal resolution. Therefore, the objective of our study was to evaluate the spatial concordance between ESI/MSI and BOLD responses for similar IEDs.

METHODS

Twenty-one patients with similar IEDs in simultaneous electroencephalogram/functional magnetic resonance imaging (EEG/fMRI) and in simultaneous EEG/magnetoencephalogram (MEG) recordings were studied. IEDs in EEG/fMRI acquisition were analyzed in an event-related paradigm within a general linear model (GLM). ESI/MSI of averaged IEDs was performed using the Maximum Entropy on the Mean. We assessed the spatial concordance between ESI/MSI and clusters of BOLD activations/deactivations with surface-based metrics.

RESULTS

ESI/MSI were concordant with one BOLD cluster for 20/21 patients (concordance with activation: 14/21 patients, deactivation: 6/21 patients, no concordance: 1/21 patients; concordance with MSI only: 3/21, ESI only: 2/21). These BOLD clusters exhibited in 19/20 cases the most significant voxel. BOLD clusters that were spatially concordant with ESI/MSI were concordant with IEDs from invasive recordings in 8/11 patients (activations: 5/8, deactivations: 3/8).

CONCLUSION

As the results of BOLD, ESI and MSI are often concordant, they reinforce our confidence in all of them. ESI and MSI confirm the most significant BOLD cluster within BOLD maps, emphasizing the importance of these clusters for the definition of the epileptic focus.

Causality within the Epileptic Network: An EEG-fMRI Study Validated by Intracranial EEG

Accurate localization of the Seizure Onset Zone (SOZ) is crucial in patients with drug-resistance focal epilepsy. EEG with fMRI recording (EEG-fMRI) has been proposed as a complementary non-invasive tool, which can give useful additional information in the pre-surgical work-up. However, fMRI maps related to interictal epileptiform activities (IED) often show multiple regions of signal change, or "networks," rather than highly focal ones. Effective connectivity approaches like Dynamic Causal Modeling (DCM) applied to fMRI data potentially offers a framework to address which brain regions drives the generation of seizures and IED within an epileptic network. Here, we present a first attempt to validate DCM on EEG-fMRI data in one patient affected by frontal lobe epilepsy. Pre-surgical EEG-fMRI demonstrated two distinct clusters of blood oxygenation level dependent (BOLD) signal increases linked to IED, one located in the left frontal pole and the other in the ipsilateral dorso-lateral frontal cortex. DCM of the IED-related BOLD signal favored a model corresponding to the left dorso-lateral frontal cortex as driver of changes in the fronto-polar region. The validity of DCM was supported by: (a) the results of two different non-invasive analysis obtained on the same dataset: EEG source imaging (ESI), and "psycho-physiological interaction" analysis; (b) the failure of a first surgical intervention limited to the fronto-polar region; (c) the results of the intracranial EEG monitoring performed after the first surgical intervention confirming a SOZ located over the dorso-lateral frontal cortex. These results add evidence that EEG-fMRI together with advanced methods of BOLD signal analysis is a promising tool that can give relevant information within the epilepsy surgery diagnostic work-up.

Electroencephalography/functional magnetic resonance imaging responses help predict surgical outcome in focal epilepsy

PURPOSE

Simultaneous electroencephalography/functional magnetic resonance imaging (EEG/fMRI) recording can noninvasively map in the whole brain the hemodynamic response following an interictal epileptic discharge. EEG/fMRI is gaining interest as a presurgical evaluation tool. This study aims to determine how hemodynamic responses related to epileptic activity can help predict surgical outcome in patients considered for epilepsy surgery.

METHODS

Thirty-five consecutive patients with focal epilepsy who had significant hemodynamic responses and eventually surgical resection, were studied. The statistical map of hemodynamic responses were generated and co-registered to postoperative anatomic imaging. Patients were classified into four groups defined by the relative relationship between the location of the maximum hemodynamic response and the resection: group 1, fully concordant; group 2, partially concordant; group 3, partially discordant; and group 4, fully discordant. These findings were correlated with surgical outcome with at least 12-month follow-up.

KEY FINDINGS

Ten patients in group 1 had the maximum t value (t-max) inside the resection; nine in group 2 had the t-max outside but close to the resection and the cluster with t-max overlapped the resection; five in group 3 had the t-max remote from resection, but with another less significant cluster in the resection; and 11 in group 4 had no response in the resection. The degree of concordance correlated largely with surgical outcome: a good surgical outcome (Engel's class I) was found in 7 of 10 patients of group 1, 4 of 9 of group 2, 3 of 5 of group 3, and only 1 of 11 of group 4. These results indicate that the partially concordant and partially discordant groups are best considered as inconclusive. In contrast, in the fully concordant and fully discordant groups, the sensitivity, specificity, positive predictive value, and negative predictive value were high, 87.5%, 76.9%, 70%, and 90.9%, respectively.

SIGNIFICANCE

This study demonstrates that hemodynamic responses related to epileptic activity can help delineate the epileptogenic region. Full concordance between maximum response and surgical resection is indicative of seizure freedom, whereas a resection leaving the maximum response intact is likely to lead to a poor outcome. EEG/fMRI is noninvasive but is limited to patients in whom interictal epileptic discharges can be recorded during the 60-90 min scan.

EEG-fMRI correlation patterns in the presurgical evaluation of focal epilepsy: a comparison with electrocorticographic data and surgical outcome measures

EEG-correlated functional MRI (EEG-fMRI) visualizes brain regions associated with interictal epileptiform discharges (IEDs). This technique images the epileptiform network, including multifocal, superficial and deeply situated cortical areas. To understand the role of EEG-fMRI in presurgical evaluation, its results should be validated relative to a gold standard. For that purpose, EEG-fMRI data were acquired for a heterogeneous group of surgical candidates (n=16) who were later implanted with subdural grids and strips (ECoG). The EEG-fMRI correlation patterns were systematically compared with brain areas involved in IEDs ECoG, using a semi-automatic analysis method, as well as to the seizure onset zone, resected area, and degree of seizure freedom. In each patient at least one of the EEG-fMRI areas was concordant with an interictally active ECoG area, always including the early onset area of IEDs in the ECoG data. This confirms that EEG-fMRI reflects a pattern of onset and propagation of epileptic activity. At group level, 76% of the BOLD regions that were covered with subdural grids, were concordant with interictally active ECoG electrodes. Due to limited spatial sampling, 51% of the BOLD regions were not covered with electrodes and could, therefore, not be validated. From an ECoG perspective it appeared that 29% of the interictally active ECoG regions were missed by EEG-fMRI and that 68% of the brain regions were correctly identified as inactive with EEG-fMRI. Furthermore, EEG-fMRI areas included the complete seizure onset zone in 83% and resected area in 93% of the data sets. No clear distinction was found between patients with a good or poor surgical outcome: in both patient groups, EEG-fMRI correlation patterns were found that were either focal or widespread. In conclusion, by comparison of EEG-fMRI with interictal invasive EEG over a relatively large patient population we were able to show that the EEG-fMRI correlation patterns are spatially accurate at the level of neurosurgical units (i.e. anatomical brain regions) and reflect the underlying network of IEDs. Therefore, we expect that EEG-fMRI can play an important role for the determination of the implantation strategy.

Mapping preictal and ictal haemodynamic networks using video-electroencephalography and functional imaging

Ictal patterns on scalp-electroencephalography are often visible only after propagation, therefore rendering localization of the seizure onset zone challenging. We hypothesized that mapping haemodynamic changes before and during seizures using simultaneous video-electroencephalography and functional imaging will improve the localization of the seizure onset zone. Fifty-five patients with ≥2 refractory focal seizures/day, and who had undergone long-term video-electroencephalography monitoring were included in the study. 'Preictal' (30 s immediately preceding the electrographic seizure onset) and ictal phases, 'ictal-onset'; 'ictalestablished' and 'late ictal', were defined based on the evolution of the electrographic pattern and clinical semiology. The functional imaging data were analysed using statistical parametric mapping to map ictal phase-related haemodynamic changes consistent across seizures. The resulting haemodynamic maps were overlaid on co-registered anatomical scans, and the spatial concordance with the presumed and invasively defined seizure onset zone was determined. Twenty patients had typical seizures during functional imaging. Seizures were identified on video-electroencephalography in 15 of 20, on electroencephalography alone in two and on video alone in three patients. All patients showed significant ictal-related haemodynamic changes. In the six cases that underwent invasive evaluation, the ictal-onset phase-related maps had a degree of concordance with the presumed seizure onset zone for all patients. The most statistically significant haemodynamic cluster within the presumed seizure onset zone was between 1.1 and 3.5 cm from the invasively defined seizure onset zone, which was resected in two of three patients undergoing surgery (Class I post-surgical outcome) and was not resected in one patient (Class III post-surgical outcome). In the remaining 14 cases, the ictal-onset phase-related maps had a degree of concordance with the presumed seizure onset zone in six of eight patients with structural-lesions and five of six non-lesional patients. The most statistically significant haemodynamic cluster was localizable at sub-lobar level within the presumed seizure onset zone in six patients. The degree of concordance of haemodynamic maps was significantly better (P < 0.05) for the ictal-onset phase [entirely concordant/concordant plus (13/20; 65%) + some concordance (4/20; 20%) = 17/20; 85%] than ictal-established [entirely concordant/concordant plus (5/13; 38%) + some concordance (4/13; 31%) = 9/13; 69%] and late ictal [concordant plus (1/9; 11%) + some concordance (4/9; 44%) = 5/9; 55%] phases. Ictal propagation-related haemodynamic changes were also seen in symptomatogenic areas (9/20; 45%) and the default mode network (13/20; 65%). A common pattern of preictal changes was seen in 15 patients, starting between 98 and 14 s before electrographic seizure onset, and the maps had a degree of concordance with the presumed seizure onset zone in 10 patients. In conclusion, preictal and ictal haemodynamic changes in refractory focal seizures can non-invasively localize seizure onset at sub-lobar/gyral level when ictal scalp-electroencephalography is not helpful.

Mapping interictal epileptic discharges using mutual information between concurrent EEG and fMRI

OBJECTIVE

The mapping of haemodynamic changes related to interictal epileptic discharges (IED) in simultaneous electroencephalography (EEG) and functional MRI (fMRI) studies is usually carried out by means of EEG-correlated fMRI analyses where the EEG information specifies the model to test on the fMRI signal. The sensitivity and specificity critically depend on the accuracy of EEG detection and the validity of the haemodynamic model. In this study we investigated whether an information theoretic analysis based on the mutual information (MI) between the presence of epileptic activity on EEG and the fMRI data can provide further insights into the haemodynamic changes related to interictal epileptic activity. The important features of MI are that: 1) both recording modalities are treated symmetrically; 2) no requirement for a-priori models for the haemodynamic response function, or assumption of a linear relationship between the spiking activity and BOLD responses, and 3) no parametric model for the type of noise or its probability distribution is necessary for the computation of MI.

METHODS

Fourteen patients with pharmaco-resistant focal epilepsy underwent EEG-fMRI and intracranial EEG and/or surgical resection with positive postoperative outcome (seizure freedom or considerable reduction in seizure frequency) was available in 7/14 patients. We used nonparametric statistical assessment of the MI maps based on a four-dimensional wavelet packet resampling method. The results of MI were compared to the statistical parametric maps obtained with two conventional General Linear Model (GLM) analyses based on the informed basis set (canonical HRF and its temporal and dispersion derivatives) and the Finite Impulse Response (FIR) models.

RESULTS

The MI results were concordant with the electro-clinically or surgically defined epileptogenic area in 8/14 patients and showed the same degree of concordance as the results obtained with the GLM-based methods in 12 patients (7 concordant and 5 discordant). In one patient, the information theoretic analysis improved the delineation of the irritative zone compared with the GLM-based methods.

DISCUSSION

Our findings suggest that an information theoretic analysis can provide clinically relevant information about the BOLD signal changes associated with the generation and propagation of interictal epileptic discharges. The concordance between the MI, GLM and FIR maps support the validity of the assumptions adopted in GLM-based analyses of interictal epileptic activity with EEG-fMRI in such a manner that they do not significantly constrain the localization of the epileptogenic zone.

Epileptic Encephalopathies with Status Epilepticus during Sleep: New Techniques for Understanding Pathophysiology and Therapeutic Options

Encephalopathy with status epilepticus during sleep (ESES) is an epileptic encephalopathy, as defined by the International League Against Epilepsy (ILAE) Task Force on Classification and Terminology, that is, a condition in which the epileptic processes themselves are believed to contribute to the disturbance in cerebral function. Clinical manifestations of ESES are heterogeneous: apart from different seizure types, they consist in combinations of cognitive, motor, and behavioural disturbances associated with a peculiar electroencephalographic pattern of paroxysmal activity significantly activated during slow sleep, which culminates in a picture of continuous spikes and waves during sleep (CSWS). The pathophysiological mechanisms underlying this condition are still incompletely understood. Establishing a clear-cut correlation between EEG abnormalities and clinical data, though interesting, is very complex. Computer-assisted EEG analyses especially if combined with functional magnetic resonance imaging (EEG-fMRI) and metabolic neuroimaging have recently emerged as useful approaches to better understand the pathophysiological processes underlying ESES. Treatment of ESES is not just limited to seizures control but it should be focused on controlling neuropsychological outcome through an improvement of the continuous epileptiform activity. General agreement on treatment guidelines is still lacking. Implementation of new techniques might allow a better understanding of the pathophysiology of ESES and could enhance therapeutics options.

Neuronal networks in west syndrome as revealed by source analysis and renormalized partial directed coherence

West syndrome is a severe epileptic encephalopathy of infancy with a poor developmental outcome. This syndrome is associated with the pathognomonic EEG feature of hypsarrhythmia. The aim of the study was to describe neuronal networks underlying hypsarrhythmia using the source analysis method (dynamic imaging of coherent sources or DICS) which represents an inverse solution algorithm in the frequency domain. In order to investigate the interaction within the detected network, a renormalized partial directed coherence (RPDC) method was also applied as a measure of the directionality of information flow between the source signals. Both DICS and RPDC were performed for EEG delta activity (1-4 Hz) in eight patients with West syndrome and in eight patients with partial epilepsies (control group). The brain area with the strongest power in the given frequency range was defined as the reference region. The coherence between this reference region and the entire brain was computed using DICS. After that, the RPDC was applied to the source signals estimated by DICS. The results of electrical source imaging were compared to results of a previous EEG-fMRI study which had been carried out using the same cohort of patients. As revealed by DICS, delta activity in hypsarrhythmia was associated with coherent sources in the occipital cortex (main source) as well as the parietal cortex, putamen, caudate nucleus and brainstem. In patients with partial epilepsies, delta activity could be attributed to sources in the occipital, parietal and sensory-motor cortex. In West syndrome, RPDC showed the strongest and most significant direction of ascending information flow from the brainstem towards the putamen and cerebral cortex. The neuronal network underlying hypsarrhythmia in this study resembles the network which was described in previous EEG-fMRI and PET studies with involvement of the brainstem, putamen and cortical regions in the generation of hypsarrhythmia. The RPDC suggests that brainstem could have a key role in the pathogenesis of West syndrome. This study supports the theory that hypsarrhythmia results from ascending brainstem pathways that project widely to basal ganglia and cerebral cortex.