Dynamic imaging of seizure activity in pediatric epilepsy patients

Estimation of seizure onset zone from ictal scalp EEG using independent component analysis in extratemporal lobe epilepsy

OBJECTIVE

The purpose of this study is to localize the seizure onset zone of patients suffering from drug-resistant epilepsy. During the last two decades, multiple studies proposed the use of Independent Component Analysis (ICA) to analyze ictal electroencephalogram (EEG) recordings. This study aims at evaluating ICA potential with quantitative measurements. In particular, we address the challenging step where the components extracted by ICA of an ictal nature must be selected.

APPROACH

We considered a cohort of 10 patients suffering from extratemporal lobe epilepsy who were rendered seizure-free after surgery. Different sets of pre-processing parameters were compared and component features were explored to help distinguish ictal components from others. Quantitative measurements were implemented to determine whether some of the components returned by ICA were located within the resection zone and thus likely to be ictal. Finally, an assistance to the component selection was proposed based on the implemented features.

MAIN RESULTS

For every seizure, at least one component returned by ICA was localized within the resection zone, with the optimal pre-processing parameters. Three features were found to distinguish components localized within the resection zone: the dispersion of their active brain sources, the ictal rhythm power and the contribution to the EEG variance. Using the implemented component selection assistance based on the features, the probability that the first proposed component yields an accurate estimation reaches 51.43% (without assistance: 24.74%). The accuracy reaches 80% when considering the best result within the first five components.

SIGNIFICANCE

This study confirms the utility of ICA for ictal EEG analysis in extratemporal lobe epilepsy, and suggests relevant features to analyze the components returned by ICA. A component selection assistance is proposed to guide clinicians in their choice for ictal components.

Ictal EEG source localization in focal epilepsy: Review and future perspectives

Electroencephalographic (EEG) source imaging localizes the generators of neural activity in the brain. During presurgical epilepsy evaluation, EEG source imaging of interictal epileptiform discharges is an established tool to estimate the irritative zone. However, the origin of interictal activity can be partly or fully discordant with the origin of seizures. Therefore, source imaging based on ictal EEG data to determine the seizure onset zone can provide precious clinical information. In this descriptive review, we address the importance of localizing the seizure onset zone based on noninvasive EEG recordings as a complementary analysis that might reduce the burden of the presurgical evaluation. We identify three major challenges (low signal-to-noise ratio of the ictal EEG data, spread of ictal activity in the brain, and validation of the developed methods) and discuss practical solutions. We provide an extensive overview of the existing clinical studies to illustrate the potential clinical utility of EEG-based localization of the seizure onset zone. Finally, we conclude with future perspectives and the needs for translating ictal EEG source imaging into clinical practice.

Electrophysiological Source Imaging: A Noninvasive Window to Brain Dynamics

Brain activity and connectivity are distributed in the three-dimensional space and evolve in time. It is important to image brain dynamics with high spatial and temporal resolution. Electroencephalography (EEG) and magnetoencephalography (MEG) are noninvasive measurements associated with complex neural activations and interactions that encode brain functions. Electrophysiological source imaging estimates the underlying brain electrical sources from EEG and MEG measurements. It offers increasingly improved spatial resolution and intrinsically high temporal resolution for imaging large-scale brain activity and connectivity on a wide range of timescales. Integration of electrophysiological source imaging and functional magnetic resonance imaging could further enhance spatiotemporal resolution and specificity to an extent that is not attainable with either technique alone. We review methodological developments in electrophysiological source imaging over the past three decades and envision its future advancement into a powerful functional neuroimaging technology for basic and clinical neuroscience applications.

Back-Projection Cortical Potential Imaging: Theory and Results

Electroencephalography (EEG) is the single brain monitoring technique that is non-invasive, portable, passive, exhibits high-temporal resolution, and gives a directmeasurement of the scalp electrical potential. Amajor disadvantage of the EEG is its low-spatial resolution, which is the result of the low-conductive skull that "smears" the currents coming from within the brain. Recording brain activity with both high temporal and spatial resolution is crucial for the localization of confined brain activations and the study of brainmechanismfunctionality, whichis then followed by diagnosis of brain-related diseases. In this paper, a new cortical potential imaging (CPI) method is presented. The new method gives an estimation of the electrical activity on the cortex surface and thus removes the "smearing effect" caused by the skull. The scalp potentials are back-projected CPI (BP-CPI) onto the cortex surface by building a well-posed problem to the Laplace equation that is solved by means of the finite elements method on a realistic head model. A unique solution to the CPI problem is obtained by introducing a cortical normal current estimation technique. The technique is based on the same mechanism used in the well-known surface Laplacian calculation, followed by a scalp-cortex back-projection routine. The BP-CPI passed four stages of validation, including validation on spherical and realistic head models, probabilistic analysis (Monte Carlo simulation), and noise sensitivity tests. In addition, the BP-CPI was compared with the minimum norm estimate CPI approach and found superior for multi-source cortical potential distributions with very good estimation results (CC >0.97) on a realistic head model in the regions of interest, for two representative cases. The BP-CPI can be easily incorporated in different monitoring tools and help researchers by maintaining an accurate estimation for the cortical potential of ongoing or event-related potentials in order to have better neurological inferences from the EEG.

Identifying epileptic source location and extent: an iterative sparse electromagnetic source imaging algorithm

In this paper we have introduced a novel electromagnetic source imaging (ESI) technique and demonstrated its validity and excellent performance in imaging the location and extent of underlying epileptic sources in patients suffering from focal epilepsy. The proposed algorithm employs ideas from sparse signal processing literature and convex optimization theories to improve source imaging results obtained from scalp-recorded electroencephalogram (EEG). EEG source imaging results generally use subjective methods to determine the extent of the underlying brain activity. The proposed technique provides significant improvement in dealing with such shortcomings and eliminates the need for thresholding. The results of our computer simulations and clinical validation study demonstrate the excellent performance of the proposed algorithm and suggest it may become a useful tool for objectively determining the location and extent of focal epileptic activity in a noninvasive fashion.

Noninvasive Electromagnetic Source Imaging and Granger Causality Analysis: An Electrophysiological Connectome (eConnectome) Approach

OBJECTIVE

Combined source-imaging techniques and directional connectivity analysis can provide useful information about the underlying brain networks in a noninvasive fashion. Source-imaging techniques have been used successfully to either determine the source of activity or to extract source time-courses for Granger causality analysis, previously. In this work, we utilize source-imaging algorithms to both find the network nodes [regions of interest (ROI)] and then extract the activation time series for further Granger causality analysis. The aim of this work is to find network nodes objectively from noninvasive electromagnetic signals, extract activation time-courses, and apply Granger analysis on the extracted series to study brain networks under realistic conditions.

METHODS

Source-imaging methods are used to identify network nodes and extract time-courses and then Granger causality analysis is applied to delineate the directional functional connectivity of underlying brain networks. Computer simulations studies where the underlying network (nodes and connectivity pattern) is known were performed; additionally, this approach has been evaluated in partial epilepsy patients to study epilepsy networks from interictal and ictal signals recorded by EEG and/or Magnetoencephalography (MEG).

RESULTS

Localization errors of network nodes are less than 5 mm and normalized connectivity errors of ∼20% in estimating underlying brain networks in simulation studies. Additionally, two focal epilepsy patients were studied and the identified nodes driving the epileptic network were concordant with clinical findings from intracranial recordings or surgical resection.

CONCLUSION

Our study indicates that combined source-imaging algorithms with Granger causality analysis can identify underlying networks precisely (both in terms of network nodes location and internodal connectivity).

SIGNIFICANCE

The combined source imaging and Granger analysis technique is an effective tool for studying normal or pathological brain conditions.

Imaging brain source extent from EEG/MEG by means of an iteratively reweighted edge sparsity minimization (IRES) strategy

Estimating extended brain sources using EEG/MEG source imaging techniques is challenging. EEG and MEG have excellent temporal resolution at millisecond scale but their spatial resolution is limited due to the volume conduction effect. We have exploited sparse signal processing techniques in this study to impose sparsity on the underlying source and its transformation in other domains (mathematical domains, like spatial gradient). Using an iterative reweighting strategy to penalize locations that are less likely to contain any source, it is shown that the proposed iteratively reweighted edge sparsity minimization (IRES) strategy can provide reasonable information regarding the location and extent of the underlying sources. This approach is unique in the sense that it estimates extended sources without the need of subjectively thresholding the solution. The performance of IRES was evaluated in a series of computer simulations. Different parameters such as source location and signal-to-noise ratio were varied and the estimated results were compared to the targets using metrics such as localization error (LE), area under curve (AUC) and overlap between the estimated and simulated sources. It is shown that IRES provides extended solutions which not only localize the source but also provide estimation for the source extent. The performance of IRES was further tested in epileptic patients undergoing intracranial EEG (iEEG) recording for pre-surgical evaluation. IRES was applied to scalp EEGs during interictal spikes, and results were compared with iEEG and surgical resection outcome in the patients. The pilot clinical study results are promising and demonstrate a good concordance between noninvasive IRES source estimation with iEEG and surgical resection outcomes in the same patients. The proposed algorithm, i.e. IRES, estimates extended source solutions from scalp electromagnetic signals which provide relatively accurate information about the location and extent of the underlying source.

Ictal EEG Source Imaging for Presurgical Evaluation of Refractory Focal Epilepsy

BACKGROUND

Electroencephalography source imaging (ESI) is a promising tool for localizing the cortical sources of both ictal and interictal epileptic activities. Many studies have shown the clinical usefulness of interictal ESI, but very few have investigated the utility of ictal ESI. The aim of this article is to examine the clinical usefulness of ictal ESI for epileptic focus localization in patients with refractory focal epilepsy, especially extratemporal lobe epilepsy.

METHODS

Both ictal and interictal ESI were performed by the use of patient-specific realistic forward models and 3 different linear distributed inverse models. Lateralization as well as concordance between ESI-estimated focuses and single-photon emission computed tomography (SPECT) focuses were assessed.

RESULTS

All the ESI focuses (both ictal and interictal) were found lateralized to the same hemisphere as ictal SPECT focuses. Lateralization results also were in agreement with the lesion sides as visualized on magnetic resonance imaging. Ictal ESI results, obtained from the best-performing inverse model, were fully concordant with the same cortical lobe as SPECT focuses, whereas the corresponding concordance rate is 87.50% in case of interictal ESI.

CONCLUSIONS

Our findings show that ictal ESI gives fully lateralized and highly concordant results with ictal SPECT and may provide a cost-effective substitute for ictal SPECT.

Thalamocortical relationship in epileptic patients with generalized spike and wave discharges--A multimodal neuroimaging study

Unlike focal or partial epilepsy, which has a confined range of influence, idiopathic generalized epilepsy (IGE) often affects the whole or a larger portion of the brain without obvious, known cause. It is important to understand the underlying network which generates epileptic activity and through which epileptic activity propagates. The aim of the present study was to investigate the thalamocortical relationship using non-invasive imaging modalities in a group of IGE patients. We specifically investigated the roles of the mediodorsal nuclei in the thalami and the medial frontal cortex in generating and spreading IGE activities. We hypothesized that the connectivity between these two structures is key in understanding the generation and propagation of epileptic activity in brains affected by IGE. Using three imaging techniques of EEG, fMRI and EEG-informed fMRI, we identified important players in generation and propagation of generalized spike-and-wave discharges (GSWDs). EEG-informed fMRI suggested multiple regions including the medial frontal area near to the anterior cingulate cortex, mediodorsal nuclei of the thalamus, caudate nucleus among others that related to the GSWDs. The subsequent seed-based fMRI analysis revealed a reciprocal cortical and bi-thalamic functional connection. Through EEG-based Granger Causality analysis using (DTF) and adaptive DTF, within the reciprocal thalamocortical circuitry, thalamus seems to serve as a stronger source in driving cortical activity from initiation to the propagation of a GSWD. Such connectivity change starts before the GSWDs and continues till the end of the slow wave discharge. Thalamus, especially the mediodorsal nuclei, may serve as potential targets for deep brain stimulation to provide more effective treatment options for patients with drug-resistant generalized epilepsy.

The performance of the spatiotemporal Kalman filter and LORETA in seizure onset localization

The assumption of spatial-smoothness is often used to solve the bioelectric inverse problem during electroencephalographic (EEG) source imaging, e.g., in low resolution electromagnetic tomography (LORETA). Since the EEG data show a temporal structure, the combination of the temporal-smoothness and the spatial-smoothness constraints may improve the solution of the EEG inverse problem. This study investigates the performance of the spatiotemporal Kalman filter (STKF) method, which is based on spatial and temporal smoothness, in the localization of a focal seizure's onset and compares its results to those of LORETA. The main finding of the study was that the STKF with an autoregressive model of order two significantly outperformed LORETA in the accuracy and consistency of the localization, provided that the source space consists of a whole-brain volumetric grid. In the future, these promising results will be confirmed using data from more patients and performing statistical analyses on the results. Furthermore, the effects of the temporal smoothness constraint will be studied using different types of focal seizures.

Effect of EEG electrode number on epileptic source localization in pediatric patients

OBJECTIVE

To investigate the relationship between EEG source localization and the number of scalp EEG recording channels.

METHODS

128 EEG channel recordings of 5 pediatric patients with medically intractable partial epilepsy were used to perform source localization of interictal spikes. The results were compared with surgical resection and intracranial recordings. Various electrode configurations were tested and a series of computer simulations based on a realistic head boundary element model were also performed in order to further validate the clinical findings.

RESULTS

The improvement seen in source localization substantially decreases as the number of electrodes increases. This finding was evaluated using the surgical resection, intracranial recordings and computer simulation. It was also shown in the simulation that increasing the electrode numbers could remedy the localization error of deep sources. A plateauing effect was seen in deep and superficial sources with further increasing the electrode number.

CONCLUSION

The source localization is improved when electrode numbers increase, but the absolute improvement in accuracy decreases with increasing electrode number.

SIGNIFICANCE

Increasing the electrode number helps decrease localization error and thus can more ably assist the physician to better plan for surgical procedures.

Assessing dynamic spectral causality by lagged adaptive directed transfer function and instantaneous effect factor

It is of significance to assess the dynamic spectral causality among physiological signals. Several practical estimators adapted from spectral Granger causality have been exploited to track dynamic causality based on the framework of time-varying multivariate autoregressive (tvMVAR) models. The nonzero covariance of the model's residuals has been used to describe the instantaneous effect phenomenon in some causality estimators. However, for the situations with Gaussian residuals in some autoregressive models, it is challenging to distinguish the directed instantaneous causality if the sufficient prior information about the "causal ordering" is missing. Here, we propose a new algorithm to assess the time-varying causal ordering of tvMVAR model under the assumption that the signals follow the same acyclic causal ordering for all time lags and to estimate the instantaneous effect factor (IEF) value in order to track the dynamic directed instantaneous connectivity. The time-lagged adaptive directed transfer function (ADTF) is also estimated to assess the lagged causality after removing the instantaneous effect. In this study, we first investigated the performance of the causal-ordering estimation algorithm and the accuracy of IEF value. Then, we presented the results of IEF and time-lagged ADTF method by comparing with the conventional ADTF method through simulations of various propagation models. Statistical analysis results suggest that the new algorithm could accurately estimate the causal ordering and give a good estimation of the IEF values in the Gaussian residual conditions. Meanwhile, the time-lagged ADTF approach is also more accurate in estimating the time-lagged dynamic interactions in a complex nervous system after extracting the instantaneous effect. In addition to the simulation studies, we applied the proposed method to estimate the dynamic spectral causality on real visual evoked potential (VEP) data in a human subject. Its usefulness in time-variant spectral causality assessment was demonstrated through the mutual causality investigation of brain activity during the VEP experiments.

Noninvasive imaging of the high frequency brain activity in focal epilepsy patients

High-frequency (HF) activity represents a potential biomarker of the epileptogenic zone in epilepsy patients, the removal of which is considered to be crucial for seizure-free surgical outcome. We proposed a high frequency source imaging (HFSI) approach to noninvasively image the brain sources of the scalp-recorded HF EEG activity. Both computer simulation and clinical patient data analysis were performed to investigate the feasibility of using the HFSI approach to image the sources of HF activity from noninvasive scalp EEG recordings. The HF activity was identified from high-density scalp recordings after high-pass filtering the EEG data and the EEG segments with HF activity were concatenated together to form repetitive HF activity. Independent component analysis was utilized to extract the components corresponding to the HF activity. Noninvasive EEG source imaging using realistic geometric boundary element head modeling was then applied to image the sources of the pathological HF brain activity. Five medically intractable focal epilepsy patients were studied and the estimated sources were found to be concordant with the surgical resection or intracranial recordings of the patients. The present study demonstrates, for the first time, that source imaging from the scalp HF activity could help to localize the seizure onset zone and provide a novel noninvasive way of studying the epileptic brain in humans. This study also indicates the potential application of studying HF activity in the presurgical planning of medically intractable epilepsy patients.

Lateralization and localization of epilepsy related hemodynamic foci using presurgical fMRI

OBJECTIVE

The aim was to develop a method for the purpose of localizing epilepsy related hemodynamic foci for patients suffering intractable focal epilepsy using task-free fMRI alone.

METHODS

We studied three groups of subjects: patients with intractable focal epilepsy, healthy volunteers performing motor tasks, and healthy volunteers in resting state. We performed spatial independent component analysis (ICA) on the fMRI alone data and developed a set of IC selection criteria to identify epilepsy related ICs. The method was then tested in the two healthy groups.

RESULTS

In seven out of the nine surgery patients, identified ICs were concordant with surgical resection. Our results were also consistent with presurgical evaluation of the remaining one patient without surgery and may explain why she was not suitable for resection treatment. In the motor task study of ten healthy subjects, our method revealed components with concordant spatial and temporal features as expected from the unilateral motor tasks. In the resting state study of seven healthy subjects, the method successfully rejected all components in four out of seven subjects as non-epilepsy related components.

CONCLUSION

These results suggest the lateralization and localization value of fMRI alone in presurgical evaluation for patients with intractable unilateral focal epilepsy.

SIGNIFICANCE

The proposed method is noninvasive in nature and easy to implement. It has the potential to be incorporated in current presurgical workup for treating intractable focal epilepsy patients.

Source localization of rhythmic ictal EEG activity: a study of diagnostic accuracy following STARD criteria

PURPOSE

Although precise identification of the seizure-onset zone is an essential element of presurgical evaluation, source localization of ictal electroencephalography (EEG) signals has received little attention. The aim of our study was to estimate the accuracy of source localization of rhythmic ictal EEG activity using a distributed source model.

METHODS

Source localization of rhythmic ictal scalp EEG activity was performed in 42 consecutive cases fulfilling inclusion criteria. The study was designed according to recommendations for studies on diagnostic accuracy (STARD). The initial ictal EEG signals were selected using a standardized method, based on frequency analysis and voltage distribution of the ictal activity. A distributed source model-local autoregressive average (LAURA)-was used for the source localization. Sensitivity, specificity, and measurement of agreement (kappa) were determined based on the reference standard-the consensus conclusion of the multidisciplinary epilepsy surgery team. Predictive values were calculated from the surgical outcome of the operated patients. To estimate the clinical value of the ictal source analysis, we compared the likelihood ratios of concordant and discordant results. Source localization was performed blinded to the clinical data, and before the surgical decision.

KEY FINDINGS

Reference standard was available for 33 patients. The ictal source localization had a sensitivity of 70% and a specificity of 76%. The mean measurement of agreement (kappa) was 0.61, corresponding to substantial agreement (95% confidence interval (CI) 0.38-0.84). Twenty patients underwent resective surgery. The positive predictive value (PPV) for seizure freedom was 92% and the negative predictive value (NPV) was 43%. The likelihood ratio was nine times higher for the concordant results, as compared with the discordant ones.

SIGNIFICANCE

Source localization of rhythmic ictal activity using a distributed source model (LAURA) for the ictal EEG signals selected with a standardized method is feasible in clinical practice and has a good diagnostic accuracy. Our findings encourage clinical neurophysiologists assessing ictal EEGs to include this method in their armamentarium.