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

Simultaneous EEG–fMRI acquisitions in patients with epilepsy often reveal distributed patterns of Blood Oxygen Level Dependant (BOLD) change correlated with epileptiform discharges. We investigated if electrical source imaging (ESI) performed on the interictal epileptiform discharges (IED) acquired during fMRI acquisition could be used to study the dynamics of the networks identified by the BOLD effect, thereby avoiding the limitations of combining results from separate recordings. Nine selected patients (13 IED types identified) with focal epilepsy underwent EEG–fMRI. Statistical analysis was performed using SPM5 to create BOLD maps. ESI was performed on the IED recorded during fMRI acquisition using a realistic head model (SMAC) and a distributed linear inverse solution (LAURA). ESI could not be performed in one case. In 10/12 remaining studies, ESI at IED onset (ESIo) was anatomically close to one BOLD cluster. Interestingly, ESIo was closest to the positive BOLD cluster with maximal statistical significance in only 4/12 cases and closest to negative BOLD responses in 4/12 cases. Very small BOLD clusters could also have clinical relevance in some cases. ESI at later time frame (ESIp) showed propagation to remote sources co-localised with other BOLD clusters in half of cases. In concordant cases, the distance between maxima of ESI and the closest EEG–fMRI cluster was less than 33 mm, in agreement with previous studies. We conclude that simultaneous ESI and EEG–fMRI analysis may be able to distinguish areas of BOLD response related to initiation of IED from propagation areas. This combination provides new opportunities for investigating epileptic networks

Simultaneous intracranial EEG and fMRI of interictal epileptic discharges in humans

Simultaneous scalp EEG–fMRI measurements allow the study of epileptic networks and more generally, of the coupling between neuronal activity and haemodynamic changes in the brain. Intracranial EEG (icEEG) has greater sensitivity and spatial specificity than scalp EEG but limited spatial sampling. We performed simultaneous icEEG and functional MRI recordings in epileptic patients to study the haemodynamic correlates of intracranial interictal epileptic discharges (IED). Two patients undergoing icEEG with subdural and depth electrodes as part of the presurgical assessment of their pharmaco-resistant epilepsy participated in the study. They were scanned on a 1.5 T MR scanner following a strict safety protocol. Simultaneous recordings of fMRI and icEEG were obtained at rest. IED were subsequently visually identified on icEEG and their fMRI correlates were mapped using a general linear model (GLM). On scalp EEG–fMRI recordings performed prior to the implantation, no IED were detected. icEEG–fMRI was well tolerated and no adverse health effect was observed. intra-MR icEEG was comparable to that obtained outside the scanner. In both cases, significant haemodynamic changes were revealed in relation to IED, both close to the most active electrode contacts and at distant sites. In one case, results showed an epileptic network including regions that could not be sampled by icEEG, in agreement with findings from magneto-encephalography, offering some explanation for the persistence of seizures after surgery. Hence, icEEG–fMRI allows the study of whole-brain human epileptic networks with unprecedented sensitivity and specificity. This could help improve our understanding of epileptic networks with possible implications for epilepsy surgery

Independent component analysis of interictal fMRI in focal epilepsy: comparison with general linear model-based EEG-correlated fMRI

The general linear model (GLM) has been used to analyze simultaneous EEG–fMRI to reveal BOLD changes linked to interictal epileptic discharges (IED) identified on scalp EEG. This approach is ineffective when IED are not evident in the EEG. Data-driven fMRI analysis techniques that do not require an EEG derived model may offer a solution in these circumstances. We compared the findings of independent components analysis (ICA) and EEG-based GLM analyses of fMRI data from eight patients with focal epilepsy. Spatial ICA was used to extract independent components (IC) which were automatically classified as either BOLD-related, motion artefacts, EPI-susceptibility artefacts, large blood vessels, noise at high spatial or temporal frequency. The classifier reduced the number of candidate IC by 78%, with an average of 16 BOLD-related IC. Concordance between the ICA and GLM-derived results was assessed based on spatio-temporal criteria. In each patient, one of the IC satisfied the criteria to correspond to IED-based GLM result. The remaining IC were consistent with BOLD patterns of spontaneous brain activity and may include epileptic activity that was not evident on the scalp EEG. In conclusion, ICA of fMRI is capable of revealing areas of epileptic activity in patients with focal epilepsy and may be useful for the analysis of EEG–fMRI data in which abnormalities are not apparent on scalp EEG

EEG correlated functional MRI and postoperative outcome in focal epilepsy

Background: The main challenge in assessing patients with epilepsy for resective surgery is localising seizure onset. Frequently, identification of the irritative and seizure onset zones requires invasive EEG. EEG correlated functional MRI (EEG-fMRI) is a novel imaging technique which may provide localising information with regard to these regions. In patients with focal epilepsy, interictal epileptiform discharge (IED) correlated blood oxygen dependent level (BOLD) signal changes were observed in approximately 50% of patients in whom IEDs are recorded. In 70%, these are concordant with expected seizure onset defined by non-invasive electroclinical information. Assessment of clinical validity requires post-surgical outcome studies which have, to date, been limited to case reports of correlation with intracranial EEG. The value of EEG-fMRI was assessed in patients with focal epilepsy who subsequently underwent epilepsy surgery, and IED correlated fMRI signal changes were related to the resection area and clinical outcome. Methods: Simultaneous EEG-fMRI was recorded in 76 patients undergoing presurgical evaluation and the locations of IED correlated preoperative BOLD signal change were compared with the resected area and postoperative outcome. Results: 21 patients had activations with epileptic activity on EEG-fMRI and 10 underwent surgical resection. Seven of 10 patients were seizure free following surgery and the area of maximal BOLD signal change was concordant with resection in six of seven patients. In the remaining three patients, with reduced seizure frequency post-surgically, areas of significant IED correlated BOLD signal change lay outside the resection. 42 of 55 patients who had no IED related activation underwent resection. Conclusion: These results show the potential value of EEG-fMRI in presurgical evaluation

EEG–fMRI mapping of asymmetrical delta activity in a patient with refractory epilepsy is concordant with the epileptogenic region determined by intracranial EEG

We studied a patient with refractory focal epilepsy using continuous EEG-correlated fMRI. Seizures were characterized by head turning to the left and clonic jerking of the left arm, suggesting a right frontal epileptogenic region. Interictal EEG showed occasional runs of independent nonlateralized slow activity in the delta band with right frontocentral dominance and had no lateralizing value. Ictal scalp EEG had no lateralizing value. Ictal scalp EEG suggested right-sided central slow activity preceding some seizures. Structural 3-T MRI showed no abnormality. There was no clear epileptiform abnormality during simultaneous EEG–fMRI. We therefore modeled asymmetrical EEG delta activity at 1–3 Hz near frontocentral electrode positions. Significant blood oxygen level-dependent (BOLD) signal changes in the right superior frontal gyrus correlated with right frontal oscillations at 1–3 Hz but not at 4–7 Hz and with neither of the two frequency bands when derived from contralateral or posterior electrode positions, which served as controls. Motor fMRI activations with a finger-tapping paradigm were asymmetrical: they were more anterior for the left hand compared with the right and were near the aforementioned EEG-correlated signal changes. A right frontocentral perirolandic seizure onset was identified with a subdural grid recording, and electric stimulation of the adjacent contact produced motor responses in the left arm and after discharges. The fMRI localization of the left hand motor and the detected BOLD activation associated with modeled slow activity suggest a role for localization of the epileptogenic region with EEG–fMRI even in the absence of clear interictal discharges

Combined EEG-fMRI and tractography to visualise propagation of epileptic activity

In a patient with refractory temporal lobe epilepsy, EEG-fMRI showed activation in association with left anterior temporal interictal discharges, in the left temporal, parietal and occipital lobes. Dynamic causal modelling suggested propagation of neural activity from the temporal focus to the area of occipital activation. Tractography showed connections from the site of temporal lobe activation to the site of occipital activation. This demonstrates the principle of combining EEG-fMRI and tractography to delineate the pathways of propagation of epileptic activity

Coupling between gamma-band power and cerebral blood volume during recurrent acute neocortical seizures

Characterization of neural and hemodynamic biomarkers of epileptic activity that can be measured using non-invasive techniques is fundamental to the accurate identification of the epileptogenic zone (EZ) in the clinical setting. Recently, oscillations at gamma-band frequencies and above (>30 Hz) have been suggested to provide valuable localizing information of the EZ and track cortical activation associated with epileptogenic processes. Although a tight coupling between gamma-band activity and hemodynamic-based signals has been consistently demonstrated in non-pathological conditions, very little is known about whether such a relationship is maintained in epilepsy and the laminar etiology of these signals. Confirmation of this relationship may elucidate the underpinnings of perfusion-based signals in epilepsy and the potential value of localizing the EZ using hemodynamic correlates of pathological rhythms. Here, we use concurrent multi-depth electrophysiology and 2-dimensional optical imaging spectroscopy to examine the coupling between multi-band neural activity and cerebral blood volume (CBV) during recurrent acute focal neocortical seizures in the urethane-anesthetized rat. We show a powerful correlation between gamma-band power (25-90 Hz) and CBV across cortical laminae, in particular layer 5, and a close association between gamma measures and multi-unit activity (MUA). Our findings provide insights into the laminar electrophysiological basis of perfusion-based imaging signals in the epileptic state and may have implications for further research using non-invasive multi-modal techniques to localize epileptogenic tissue

Modelling large motion events in fMRI studies of patients with epilepsy

EEG-correlated fMRI can provide localisation information on the generators of epileptiform discharges in patients with focal epilepsy. To increase the technique's clinical potential, it is important to consider ways of optimising the yield of each experiment while minimizing the risk of false-positive activation. Head motion can lead to severe image degradation and result in false-positive activation and is usually worse in patients than in healthy subjects. We performed general linear model fMRI data analysis on simultaneous EEG–fMRI data acquired in 34 cases with focal epilepsy. Signal changes associated with large inter-scan motion events (head jerks) were modelled using modified design matrices that include ‘scan nulling’ regressors. We evaluated the efficacy of this approach by mapping the proportion of the brain for which F-tests across the additional regressors were significant. In 95% of cases, there was a significant effect of motion in 50% of the brain or greater; for the scan nulling effect, the proportion was 36%; this effect was predominantly in the neocortex. We conclude that careful consideration of the motion-related effects in fMRI studies of patients with epilepsy is essential and that the proposed approach can be effective

Noncanonical spike-related BOLD responses in focal epilepsy

Till now, most studies of the Blood Oxygen Level-Dependent (BOLD) response to interictal epileptic discharges (IED) have assumed that its time course matches closely to that of brief physiological stimuli, commonly called the canonical event-related haemodynamic response function (canonical HRF). Analyses based on that assumption have produced significant response patterns that are generally concordant with prior electroclinical data. In this work, we used a more flexible model of the event-related response, a Fourier basis set, to investigate the presence of other responses in relation to individual IED in 30 experiments in patients with focal epilepsy. We found significant responses that had a noncanonical time course in 37% of cases, compared with 40% for the conventional, canonical HRF-based approach. In two cases, the Fourier analysis suggested activations where the conventional model did not. The noncanonical activations were almost always remote from the presumed generator of epileptiform activity. In the majority of cases with noncanonical responses, the noncanonical responses in single-voxel clusters were suggestive of artifacts. We did not find evidence for IED-related noncanonical HRFs arising from areas of pathology, suggesting that the BOLD response to IED is primarily canonical. Noncanonical responses may represent a number of phenomena, including artefacts and propagated epileptiform activity

Implementation and evaluation of simultaneous video-electroencephalography and functional magnetic resonance imaging

The objective of this study was to demonstrate that the addition of simultaneous and synchronised video to electroencephalography (EEG)-correlated functional magnetic resonance imaging (fMRI) could increase recorded information without data quality reduction. We investigated the effect of placing EEG, video equipment and their required power supplies inside the scanner room, on EEG, video and MRI data quality, and evaluated video-EEG-fMRI by modelling a hand motor task. Gradient-echo, echo-planner images (EPI) were acquired on a 3-T MRI scanner at variable camera positions in a test object [with and without radiofrequency (RF) excitation], and human subjects. EEG was recorded using a commercial MR-compatible 64-channel cap and amplifiers. Video recording was performed using a two-camera custom-made system with EEG synchronization. An in-house script was used to calculate signal to fluctuation noise ratio (SFNR) from EPI in test object with variable camera positions and in human subjects with and without concurrent video recording. Five subjects were investigated with video-EEG-fMRI while performing hand motor task. The fMRI time series data was analysed using statistical parametric mapping, by building block design general linear models which were paradigm prescribed and video based. Introduction of the cameras did not alter the SFNR significantly, nor did it show any signs of spike noise during RF off conditions. Video and EEG quality also did not show any significant artefact. The Statistical Parametric Mapping{T} maps from video based design revealed additional blood oxygen level-dependent responses in the expected locations for non-compliant subjects compared to the paradigm prescribed design. We conclude that video-EEG-fMRI set up can be implemented without affecting the data quality significantly and may provide valuable information on behaviour to enhance the analysis of fMRI data