Are Epilepsy-Related fMRI Components Dependent on the Presence of Interictal Epileptic Discharges in Scalp EEG?

Spatial independent component analysis (ICA) is increasingly being used to extract resting-state networks from fMRI data. Previous studies showed that ICA also reveals independent components (ICs) related to the seizure onset zone. However, it is currently unknown how these epileptic ICs depend on the presence of interictal epileptic discharges (IEDs) in the EEG. The goal of this study was to explore the relation between ICs obtained from fMRI epochs during the occurrence of IEDs in the EEG and those without IEDs. fMRI data sets with co-registered EEG were retrospectively selected of patients from whom the location of the epileptogenic zone was confirmed by outcome of surgery (n = 8). The fMRI data were split into two epochs: one with IEDs visible in scalp EEG and one without. Spatial ICA was applied to the fMRI data of each part separately. The maps of all resulting components were compared to the resection area and the EEG-fMRI correlation pattern by computing a spatial correlation coefficient to detect the epilepsy-related component. For all patients, except one, there was a remarkable resemblance between the epilepsy-related components selected during epochs with IEDs and those without IEDs. These findings suggest that epilepsy-related ICs are not dependent on the presence of IEDs in scalp EEG. Since these epileptic ICs showed partial overlap with resting-state networks of healthy volunteers (n = 10), our study supports the need for new ways to classify epileptic ICs.

Detection of short-lasting and ictal spike-and-wave discharges in around-the-ears EEG recordings in children with absence epilepsy.

PURPOSE: Long-term ambulatory EEG recordings can improve the monitoring of absence epilepsy in children, but signal quality and increased review workload are a concern. We evaluated the feasibility of around-the-ears EEG arrays (cEEGrids) to capture 3-Hz short-lasting and ictal spike-and-wave discharges and assessed the performance of automated detection software in cEEGrids data. We compared patterns of bilateral synchronisation between short-lasting and ictal spike-and-wave discharges. METHODS: We recruited children with suspected generalised epilepsy undergoing routine video-EEG monitoring and performed simultaneous cEEGrids recordings. We used ASSYST software to detect short-lasting 3-Hz spike-and-wave discharges (1-3 s) and ictal spike-and-wave discharges in the cEEGrids data. We assessed data quality and sensitivity of cEEGrids for spike-and-wave discharges in routine EEG. We determined the sensitivity and false detection rate for automated spike-and-wave discharge detection in cEEGrids data. We compared bihemispheric synchrony across the onset of short-lasting and ictal spike-and-wave discharges using the mean phase coherence in the 2-4 Hz frequency band. RESULTS: We included nine children with absence epilepsy (median age = 11 y, range 8-15 y, nine females) and recorded 4 h and 27 min of cEEGrids data. The recordings from seven participants were suitable for quantitative analysis, containing 82 spike-and-wave discharges. The cEEGrids captured 58 % of all spike-and-wave discharges (median individual sensitivity: 100 %, range: 47-100 %). ASSYST detected 82 % of all spike-and-wave discharges (median: 100 %, range: 41-100 %) with a false detection rate of 48/h (median: 6/h, range: 0-154/h). The mean phase coherence significantly increased during short-lasting and ictal spike-and-wave discharges in the 500-ms pre-onset to 1-s post-onset interval. CONCLUSIONS: cEEGrids are of variable quality for monitoring spike-and-wave discharges in children with absence epilepsy. ASSYST could facilitate the detection of short-lasting and ictal spike-and-wave discharges with clear periodic structures but with low specificity. A similar course of bihemispheric synchrony between short-lasting and ictal spike-and-wave discharges indicates that cortico-thalamic driving may be relevant for both types of spike-and-wave discharges.

Novel artefact removal algorithms for co-registered EEG/fMRI based on selective averaging and subtraction.

Co-registered EEG and functional MRI (EEG/fMRI) is a potential clinical tool for planning invasive EEG in patients with epilepsy. In addition, the analysis of EEG/fMRI data provides a fundamental insight into the precise physiological meaning of both fMRI and EEG data. Routine application of EEG/fMRI for localization of epileptic sources is hampered by large artefacts in the EEG, caused by switching of scanner gradients and heartbeat effects. Residuals of the ballistocardiogram (BCG) artefacts are similarly shaped as epileptic spikes, and may therefore cause false identification of spikes. In this study, new ideas and methods are presented to remove gradient artefacts and to reduce BCG artefacts of different shapes that mutually overlap in time. Gradient artefacts can be removed efficiently by subtracting an average artefact template when the EEG sampling frequency and EEG low-pass filtering are sufficient in relation to MR gradient switching (Gonçalves et al., 2007). When this is not the case, the gradient artefacts repeat themselves at time intervals that depend on the remainder between the fMRI repetition time and the closest multiple of the EEG acquisition time. These repetitions are deterministic, but difficult to predict due to the limited precision by which these timings are known. Therefore, we propose to estimate gradient artefact repetitions using a clustering algorithm, combined with selective averaging. Clustering of the gradient artefacts yields cleaner EEG for data recorded during scanning of a 3T scanner when using a sampling frequency of 2048 Hz. It even gives clean EEG when the EEG is sampled with only 256 Hz. Current BCG artefacts-reduction algorithms based on average template subtraction have the intrinsic limitation that they fail to deal properly with artefacts that overlap in time. To eliminate this constraint, the precise timings of artefact overlaps were modelled and represented in a sparse matrix. Next, the artefacts were disentangled with a least squares procedure. The relevance of this approach is illustrated by determining the BCG artefacts in a data set consisting of 29 healthy subjects recorded in a 1.5 T scanner and 15 patients with epilepsy recorded in a 3 T scanner. Analysis of the relationship between artefact amplitude, duration and heartbeat interval shows that in 22% (1.5T data) to 30% (3T data) of the cases BCG artefacts show an overlap. The BCG artefacts of the EEG/fMRI data recorded on the 1.5T scanner show a small negative correlation between HBI and BCG amplitude. In conclusion, the proposed methodology provides a substantial improvement of the quality of the EEG signal without excessive computer power or additional hardware than standard EEG-compatible equipment.

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.

A semi-automatic method to determine electrode positions and labels from gel artifacts in EEG/fMRI-studies.

The analysis of simultaneous EEG and fMRI data is generally based on the extraction of regressors of interest from the EEG, which are correlated to the fMRI data in a general linear model setting. In more advanced approaches, the spatial information of EEG is also exploited by assuming underlying dipole models. In this study, we present a semi automatic and efficient method to determine electrode positions from electrode gel artifacts, facilitating the integration of EEG and fMRI in future EEG/fMRI data models. In order to visualize all electrode artifacts simultaneously in a single view, a surface rendering of the structural MRI is made using a skin triangular mesh model as reference surface, which is expanded to a "pancake view". Then the electrodes are determined with a simple mouse click for each electrode. Using the geometry of the skin surface and its transformation to the pancake view, the 3D coordinates of the electrodes are reconstructed in the MRI coordinate frame. The electrode labels are attached to the electrode positions by fitting a template grid of the electrode cap in which the labels are known. The correspondence problem between template and sample electrodes is solved by minimizing a cost function over rotations, shifts and scalings of the template grid. The crucial step here is to use the solution of the so-called "Hungarian algorithm" as a cost function, which makes it possible to identify the electrode artifacts in arbitrary order. The template electrode grid has to be constructed only once for each cap configuration. In our implementation of this method, the whole procedure can be performed within 15 min including import of MRI, surface reconstruction and transformation, electrode identification and fitting to template. The method is robust in the sense that an electrode template created for one subject can be used without identification errors for another subject for whom the same EEG cap was used. Furthermore, the method appears to be robust against spurious or missing artifacts. We therefore consider the proposed method as a useful and reliable tool within the larger toolbox required for the analysis of co-registered EEG/fMRI data.

Correction for pulse height variability reduces physiological noise in functional MRI when studying spontaneous brain activity.

EEG correlated functional MRI (EEG-fMRI) allows the delineation of the areas corresponding to spontaneous brain activity, such as epileptiform spikes or alpha rhythm. A major problem of fMRI analysis in general is that spurious correlations may occur because fMRI signals are not only correlated with the phenomena of interest, but also with physiological processes, like cardiac and respiratory functions. The aim of this study was to reduce the number of falsely detected activated areas by taking the variation in physiological functioning into account in the general linear model (GLM). We used the photoplethysmogram (PPG), since this signal is based on a linear combination of oxy- and deoxyhemoglobin in the arterial blood, which is also the basis of fMRI. We derived a regressor from the variation in pulse height (VIPH) of PPG and added this regressor to the GLM. When this regressor was used as predictor it appeared that VIPH explained a large part of the variance of fMRI signals acquired from five epilepsy patients and thirteen healthy volunteers. As a confounder VIPH reduced the number of activated voxels by 30% for the healthy volunteers, when studying the generators of the alpha rhythm. Although for the patients the number of activated voxels either decreased or increased, the identification of the epileptogenic zone was substantially enhanced in one out of five patients, whereas for the other patients the effects were smaller. In conclusion, applying VIPH as a confounder diminishes physiological noise and allows a more reliable interpretation of fMRI results.

Burden of disease of people with epilepsy during an optimized diagnostic trajectory: costs and quality of life.

BACKGROUND: Diagnosing epilepsy can be lengthy and stressful, potentially leading to increased use of healthcare resources and a reduction in quality of life. AIM: This study aims to determine cost and quality of life before and after an optimized diagnostic procedure for people suspected of having epilepsy from a societal perspective with a follow-up of 12 months. In addition, this study aims to differentiate between people diagnosed with epilepsy during the follow-up of the study and the people who are diagnosed as not having epilepsy or for whom diagnosis is still uncertain. METHODS: A questionnaire regarding the use of healthcare resources was used accompanied by the EQ-5D-3 L. Multiple imputations by chained equations with predictive mean matching was used to account for missing data. To investigate the uncertainty of the results, non-parametric bootstrapped (1000 times) was used. RESULTS: In total, 116 people were included in the study. Total average costs per patient made in the previous 3 months had decreased from €4594 before the optimized diagnostic trajectory to €2609 in the 12 months after the optimized diagnostic trajectory. Healthcare costs were the largest expense group (52-66%) and had decreased significantly from baseline measurement to 12 months after baseline (€2395 vs €1581). Productivity costs had decreased from €1367 to €442 per 3 months. Total annual costs were similar between people diagnosed with epilepsy during the follow-up of the study and the people who are diagnosed as not having epilepsy or for whom diagnosis is still uncertain. Quality of Life had significantly increased over the course of 12 months from 0.80 to 0.84 (Dutch tariff). DISCUSSION: This study indicates that an optimized diagnostic trajectory has positively influenced the use of healthcare resources and the quality of life in people with epilepsy. As chronic care patients make diverse costs, future research should identify the long-term costs after an optimized diagnostic trajectory for patients with epilepsy, possibly identifying patients who are at high risk of becoming high-cost users in the future for early intervention.

Spatiotemporal mapping of interictal epileptiform discharges in human absence epilepsy: A MEG study.

PURPOSE: Although absence epilepsy is considered to be a prototypic type of generalized epilepsy, it is still under debate whether generalized 3 Hz spike-and-wave discharges (SWDs) might have a cortical focal origin. Here it is investigated whether focal interictal epileptiform discharges (IEDs), which typically occur in the electro- (EEG) and magnetoencephalogram (MEG) in case of focal epilepsy, are present in the MEG of children with absence epilepsy. Next, the location of the sources of the IEDs is established, and it is investigated whether the location is concordant to the earlier established focal cortical regions involved in the generalized SWDs of these children. METHODS: Whole head MEG recordings of seven children with absence epilepsy were reviewed with respect to the presence of IEDs (spikes and sharp waves). These IEDs were grouped into distinct clusters, in which each contribution to a cluster yields a comparable magnetic field distribution. Source localization was then performed onto the average signal of each cluster using an equivalent current dipole model and a realistic head model of the cortical surface. RESULTS: IEDs were detected in 6 out of 7 patients. Source reconstruction indicated most often frontal, central or parietal origins of the IED in either the left and or right hemisphere. Spatiotemporal assessment of the IEDs indicated a stable location of the averages of these discharges, indicating a single underlying cortical source. DISCUSSION: The outcome of this pilot study shows that MEG is well suited for the detection of IEDs and suggests that their estimated sources coincide rather well with the cortical regions involved during the spikes of the SWDs. It is discussed whether the presence of IEDs, classically seen as a marker of focal epilepsies, indicate that absence epilepsy should be considered as a focal type of epilepsy, in which changes in the network are evolving rapidly.

Comparison of analytical strategies for EEG-correlated fMRI data in patients with epilepsy.

The simultaneous recording of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) can be used to localize interictal epileptiform discharges (IEDs). Previous studies have reported varying degrees of concordance of EEG-fMRI with electroclinical findings. The aim of the present study is to evaluate to what extent this variability is determined by the analytical strategy or by the properties of the EEG data. For that purpose, 42 IED sets obtained in 29 patients with epilepsy were reanalyzed using a finite impulse response approach, which estimates the hemodynamic response function (HRF) from the data and allows non-causal effects. Cardiac effects were treated as additional confounders in the model. This approach was compared to the classical approach assuming a fixed HRF for each voxel in the brain. The performance of each method was assessed by comparing the fMRI results to the EEG focus. The flexible model revealed more significantly activated voxels, which resulted in more activated brain regions concordant with the EEG focus (26 vs. 16). Correction for cardiac effects improved the results in 7 out of the 42 data sets. Furthermore, design theory for event-related experiments was applied in order to determine the influence of the number of IEDs and their temporal distribution on the success of an experiment. It appeared that this success is highly dependent upon the number of IEDs present during the recording and less on their temporal spacing. We conclude that the outcome of EEG-fMRI can be improved by using an optimized analytical strategy, but also depends on the number of IEDs occurring during the recording.