SEEG in 3D: Interictal Source Localization From Intracerebral Recordings

Stereo-EEG propagating source reconstruction identifies new surgical targets for epilepsy patients

Epilepsy surgery can eliminate seizures in patients with drug-resistant focal epilepsy. Surgical intervention requires proper identification of the epileptic network and often involves implanting stereo-EEG electrodes in patients where non-invasive methods are insufficient. However, only ∼60% of patients achieve seizure-freedom following surgery. Quantitative methods have been developed to help improve surgical outcomes. However, previous quantitative methods that localized interictal spike and seizure activity using stereo-EEG recordings did not account for the propagation path encoded by the temporal dynamics of stereo-EEG recordings. Reconstructing the seizure propagation path can aid in determining whether a signal originated from the seizure onset or propagation zone, which directly informs treatment decisions. We developed a novel source reconstruction algorithm, Temporally Dependent Iterative Expansion (TEDIE), that accurately reconstructs propagating and expanding neural sources over time. TEDIE iteratively optimizes the number, location and size of neural sources to minimize the differences between the reconstructed and recorded stereo-EEG signals using temporal information to refine the reconstructions. The TEDIE output comprises a movie of seizure activity projected onto patient-specific brain anatomy. We analysed data from 46 epilepsy patients implanted with stereo-EEG electrodes at Duke Hospital (12 patients) and the Hospital of the University of Pennsylvania (34 patients). We reconstructed seizure recordings and found that TEDIE's seizure onset zone reconstructions were closer to the resected brain region for Engel 1 compared to Engel 2-4 patients, retrospectively validating the clinical utility of TEDIE. We also demonstrated that TEDIE has prospective clinical value, whereby metrics that can be determined presurgically accurately predict whether a patient would achieve seizure-freedom following surgery. Furthermore, we used TEDIE to delineate new potential surgical targets in 12/23 patients who are currently Engel 2-4. We validated TEDIE by accurately reconstructing various dynamic synthetic neural sources with known locations and sizes. TEDIE generated more accurate, focal and interpretable dynamic reconstructions of seizures compared to other algorithms (sLORETA and IRES). Our findings demonstrate that TEDIE is a promising clinical tool that can greatly improve epileptogenic zone localization and epilepsy surgery outcomes.

Accuracy of SEEG Source Localization: A Pilot Study Using Corticocortical Evoked Potentials

INTRODUCTION

EEG source localization is an established technique for localizing scalp EEG in medically refractory epilepsy but has not been adequately studied with intracranial EEG (iEEG). Differences in sensor location and spatial sampling may affect the accuracy of EEG source localization with iEEG. Corticocortical evoked potentials can be used to evaluate EEG source localization algorithms for iEEG given the known source location.

METHODS

We recorded 205 sets of corticocortical evoked potentials using low-frequency single-pulse electrical stimulation in four patients with iEEG. Averaged corticocortical evoked potentials were analyzed using 11 distributed source algorithms and compared using the Wilcoxon signed-rank test ( P < 0.05). We measured the localization error from stimulated electrodes and the spatial dispersion of each solution.

RESULTS

Minimum norm, standard low-resolution electromagnetic tomography (sLORETA), LP Norm, sLORETA-weighted accurate minimum norm (SWARM), exact LORETA (eLORETA), standardized weighted LORETA (swLORETA), and standardized shrinking LORETA-FOCUSS (ssLOFO) had the least localization error (13.3-15.7 mm) and were superior to focal underdetermined system solver (FOCUSS), logistic autoregressive average (LAURA, and LORETA, 17.9-21.7, P < 0.001). The FOCUSS solution had the smallest spatial dispersion (7.4 mm), followed by minimum norm, L1 norm, LP norm, and SWARM (20.8-28.3 mm). Gray matter stimulations had less localization error than white matter (median differences 3.1-6.1 mm) across all algorithms except SWARM, LORETA, and logistic autoregressive average. A multivariate linear regression showed that distance from the source to sensors and gray/white matter stimulation had a significant effect on localization error for some algorithms but not SWARM, minimum norm, focal underdetermined system solver, logistic autoregressive average, and LORETA.

CONCLUSIONS

Our study demonstrated that minimum norm, L1 norm, LP norm, and SWARM localize iEEG corticocortical evoked potentials well with lower localization error and spatial dispersion. Larger studies are needed to confirm these findings.

Intracranial Disease Region Composite Interpretation Technology for Enhanced Source Localization in Pediatric Epilepsy Surgery

Electroencephalography (EEG) based source localization (ESL) is a useful method to localize the epileptogenic zone in epilepsy surgery. However, previous techniques only perform 3-dimensional (3D) reconstruction, and do not conduct delineation on the cortex surface as a resection guidance, and there is very little data on intracranial EEG and pediatric cases. This study proposes an Intracranial Disease-region Composite-interpretation (IDC) EEG-based source localization (ESL) scheme that uses 3D extended reality (XR) edge computing to enhance visualization and comprehensive interpretation of intracranial EEG-based source localization (iESL) for patients with pediatric epilepsy. The proposed IDC-ESL method was effective in predicting the surgical outcome in patients with focal epilepsy, which can be effectively used for epilepsy surgery. Seizure freedom was clearly associated with complete resection of combined EEG features of interictal spike, high-frequency oscillation (HFO), and seizure onset zone (SOZ), and it had the highest significance in localizing the epileptogenic zone. However, for patients with Lennox-Gastaut syndrome (LGS), IDC-ESL was not performed effectively because of a deeply seated lesion and multifocal abnormalities. It could only roughly estimate the affected area, mainly because of insular involvement. Cautious interpretation based on intraoperative electrocorticography (ECoG) is required for accurate insular resection, particularly for LGS cases.

A spatial perturbation framework to validate implantation of the epileptogenic zone

Stereo-electroencephalography (SEEG) is the gold standard to delineate surgical targets in focal drug-resistant epilepsy. SEEG uses electrodes placed directly into the brain to identify the seizure-onset zone (SOZ). However, its major constraint is limited brain coverage, potentially leading to misidentification of the 'true' SOZ. Here, we propose a framework to assess adequate SEEG sampling by coupling epileptic biomarkers with their spatial distribution and measuring the system's response to a perturbation of this coupling. We demonstrate that the system's response is strongest in well-sampled patients when virtually removing the measured SOZ. We then introduce the spatial perturbation map, a tool that enables qualitative assessment of the implantation coverage. Probability modelling reveals a higher likelihood of well-implanted SOZs in seizure-free patients or non-seizure free patients with incomplete SOZ resections, compared to non-seizure-free patients with complete resections. This highlights the framework's value in sparing patients from unsuccessful surgeries resulting from poor SEEG coverage.

Distributed source modeling of stereoencephalographic measurements of ictal activity

OBJECTIVES

Stereoelectroencephalography (SEEG) can define the epileptogenic zone (EZ). However, SEEG is susceptible to the sampling bias, where no SEEG recording is taken within a circumscribed EZ.

METHODS

Nine patients with medically refractory epilepsy underwent SEEG recording, and brain resection got positive outcomes. Ictal neuronal currents were estimated by distributed source modeling using the SEEG data and individual's anatomical magnetic resonance imaging. Using a retrospective leave-one-out data sub-sampling, we evaluated the sensitivity and specificity of the current estimates using MRI after surgical resection or radio-frequency ablation.

RESULTS

The sensitivity and specificity in detecting the EZ were indistinguishable from either the data from all electrodes or the sub-sampled data (rank sum test: rank sum = 23719, p = 0.13) when at least one remaining electrode contact was no more than 20 mm away.

CONCLUSIONS

The distributed neuronal current estimates of ictal SEEG data can mitigate the challenge of delineating the boundary of the EZ in cases of missing an electrode implanted within the EZ and a required second SEEG exploration.

SIGNIFICANCE

Distributed source modeling can be a tool for clinicians to infer the EZ by allowing for more flexible planning of the electrode implantation route and minimizing the number of electrodes.

Semiautomated electric source imaging determines epileptogenicity of encephaloceles in temporal lobe epilepsy

OBJECTIVE

We aimed to assess the ability of semiautomated electric source imaging (ESI) from long-term video-electroencephalographic (EEG) monitoring (LTM) to determine the epileptogenicity of temporopolar encephaloceles (TEs) in patients with temporal lobe epilepsy.

METHODS

We conducted a retrospective study involving 32 temporal lobe epilepsy patients with TEs as potentially epileptogenic lesions in structural magnetic resonance imaging scans. Findings were validated through invasive intracerebral stereo-EEG in six of 32 patients and postsurgical outcome after tailored resection of the TE in 17 of 32 patients. LTM (mean duration = 6 days) was performed using the 10/20 system with additional T1/T2 for all patients and sphenoidal electrodes in 23 of 32 patients. Semiautomated detection and clustering of interictal epileptiform discharges (IEDs) were carried out to create IED types. ESI was performed on the averages of the two most frequent IED types per patient, utilizing individual head models, and two independent inverse methods (sLORETA [standardized low-resolution brain electromagnetic tomography], MUSIC [multiple signal classification]). ESI maxima concordance and propagation in spatial relation to TEs were quantified for sources with good signal quality (signal-to-noise ratio > 2, explained signal > 60%).

RESULTS

ESI maxima correctly colocalized with a TE in 20 of 32 patients (62.5%) either at the onset or half-rising flank of at least one IED type per patient. ESI maxima showed propagation from the temporal pole to other temporal or extratemporal regions in 14 of 32 patients (44%), confirming propagation originating in the area of the TE. The findings from both inverse methods validated each other in 14 of 20 patients (70%), and sphenoidal electrodes exhibited the highest signal amplitudes in 17 of 23 patients (74%). The concordance of ESI with the TE predicted a seizure-free postsurgical outcome (Engel I vs. >I) with a diagnostic odds ratio of 2.1.

SIGNIFICANCE

Semiautomated ESI from LTM often successfully identifies the epileptogenicity of TEs and the IED onset zone within the area of the TEs. Additionally, it shows potential predictive power for postsurgical outcomes in these patients.

Optimization of patient-specific stereo-EEG recording sensitivity

Stereo-EEG is a minimally invasive technique used to localize the origin of epileptic activity (the epileptogenic zone) in patients with drug-resistant epilepsy. However, current stereo-EEG trajectory planning methods are agnostic to the spatial recording sensitivity of implanted electrodes. In this study, we used image-based patient-specific computational models to design optimized stereo-EEG electrode configurations. Patient-specific optimized electrode configurations exhibited substantially higher recording sensitivity than clinically implanted configurations, and this may lead to a more accurate delineation of the epileptogenic zone. The optimized configurations also achieved equally good or better recording sensitivity with fewer electrodes compared with clinically implanted configurations, and this may reduce the risk for complications, including intracranial haemorrhage. This approach improves localization of the epileptogenic zone by transforming the clinical use of stereo-EEG from a discrete ad hoc sampling to an intelligent mapping of the regions of interest.

Electrical Source Imaging of Somatosensory Evoked Potentials from Intracranial EEG Signals

Stereoelectroencephalography (SEEG) records electrical brain activity with intracerebral electrodes. However, it has an inherently limited spatial coverage. Electrical source imaging (ESI) infers the position of the neural generators from the recorded electric potentials, and thus, could overcome this spatial undersampling problem. Here, we aimed to quantify the accuracy of SEEG ESI under clinical conditions. We measured the somatosensory evoked potential (SEP) in SEEG and in high-density EEG (HD-EEG) in 20 epilepsy surgery patients. To localize the source of the SEP, we employed standardized low resolution brain electromagnetic tomography (sLORETA) and equivalent current dipole (ECD) algorithms. Both sLORETA and ECD converged to similar solutions. Reflecting the large differences in the SEEG implantations, the localization error also varied in a wide range from 0.4 to 10 cm. The SEEG ESI localization error was linearly correlated with the distance from the putative neural source to the most activated contact. We show that it is possible to obtain reliable source reconstructions from SEEG under realistic clinical conditions, provided that the high signal fidelity recording contacts are sufficiently close to the source of the brain activity.

Reconstruction and localization of auditory sources from intracerebral SEEG using independent component analysis

Stereo-electroencephalography (SEEG) is the surgical implantation of electrodes in the brain to better localize the epileptic network in pharmaco-resistant epileptic patients. This technique has exquisite spatial and temporal resolution. Still, the number and the position of the electrodes in the brain is limited and determined by the semiology and/or preliminary non-invasive examinations, leading to a large number of unexplored brain structures in each patient. Here, we propose a new approach to reconstruct the activity of non-sampled structures in SEEG, based on independent component analysis (ICA) and dipole source localization. We have tested this approach with an auditory stimulation dataset in ten patients. The activity directly recorded from the auditory cortex served as ground truth and was compared to the ICA applied on all non-auditory electrodes. Our results show that the activity from the auditory cortex can be reconstructed at the single trial level from contacts as far as ∼40 mm from the source. Importantly, this reconstructed activity is localized via dipole fitting in the proximity of the original source. In addition, we show that the size of the confidence interval of the dipole fitting is a good indicator of the reliability of the result, which depends on the geometry of the SEEG implantation. Overall, our approach allows reconstructing the activity of structures far from the electrode locations, partially overcoming the spatial sampling limitation of intracerebral recordings.

Automated seizure onset zone locator from resting-state functional MRI in drug-resistant epilepsy

OBJECTIVE

Accurate localization of a seizure onset zone (SOZ) from independent components (IC) of resting-state functional magnetic resonance imaging (rs-fMRI) improves surgical outcomes in children with drug-resistant epilepsy (DRE). Automated IC sorting has limited success in identifying SOZ localizing ICs in adult normal rs-fMRI or uncategorized epilepsy. Children face unique challenges due to the developing brain and its associated surgical risks. This study proposes a novel SOZ localization algorithm (EPIK) for children with DRE.

METHODS

EPIK is developed in a phased approach, where fMRI noise-related biomarkers are used through high-fidelity image processing techniques to eliminate noise ICs. Then, the SOZ markers are used through a maximum likelihood-based classifier to determine SOZ localizing ICs. The performance of EPIK was evaluated on a unique pediatric DRE dataset (n = 52). A total of 24 children underwent surgical resection or ablation of an rs-fMRI identified SOZ, concurrently evaluated with an EEG and anatomical MRI. Two state-of-art techniques were used for comparison: (a) least squares support-vector machine and (b) convolutional neural networks. The performance was benchmarked against expert IC sorting and Engel outcomes for surgical SOZ resection or ablation. The analysis was stratified across age and sex.

RESULTS

EPIK outperformed state-of-art techniques for SOZ localizing IC identification with a mean accuracy of 84.7% (4% higher), a precision of 74.1% (22% higher), a specificity of 81.9% (3.2% higher), and a sensitivity of 88.6% (16.5% higher). EPIK showed consistent performance across age and sex with the best performance in those < 5 years of age. It helped achieve a ~5-fold reduction in the number of ICs to be potentially analyzed during pre-surgical screening.

SIGNIFICANCE

Automated SOZ localization from rs-fMRI, validated against surgical outcomes, indicates the potential for clinical feasibility. It eliminates the need for expert sorting, outperforms prior automated methods, and is consistent across age and sex.

Ideal current dipoles are appropriate source representations for simulating neurons for intracranial recordings

OBJECTIVE

To determine whether dipoles are an appropriate simplified representation of neural sources for stereo-EEG (sEEG).

METHODS

We compared the distributions of voltages generated by a dipole, biophysically realistic cortical neuron models, and extended regions of cortex to determine how well a dipole represented neural sources at different spatial scales and at electrode to neuron distances relevant for sEEG. We also quantified errors introduced by the dipole approximation of neural sources in sEEG source localization using standardized low-resolution electrotomography (sLORETA).

RESULTS

For pyramidal neurons, the coefficient of correlation between voltages generated by a dipole and neuron model were > 0.9 for distances > 1 mm. For small regions of cortex (∼0.1 cm2), the error in voltages between a dipole and region was < 100 µV for all distances. However, larger regions of active cortex (>5 cm2) yielded > 50 µV errors within 1.5 cm of an electrode when compared to single dipoles. Finally, source localization errors were < 5 mm when using dipoles to represent realistic neural sources.

CONCLUSIONS

Single dipoles are an appropriate source model to represent both single neurons and small regions of active cortex, while multiple dipoles are required to represent large regions of cortex.

SIGNIFICANCE

Dipoles are computationally tractable and valid source models for sEEG.

Source localization of ictal SEEG to predict postoperative seizure outcome

OBJECTIVE

Stereo-electroencephalography (SEEG) is inherently-three-dimensional and can be modeled using source localization. This study aimed to assess the validity of ictal SEEG source localization.

METHODS

The dominant frequency at ictal onset was used for source localization in the time and frequency domains using rotating dipoles and current density maps. Validity was assessed by concordance with the epileptologist-defined seizure onset zone (conventional SOZ) and the surgical treatment volume (TV) of seizure-free versus non-seizure-free patients.

RESULTS

Source localization was performed on 68 seizures from 27 patients. Median distance to nearest contact in the conventional SOZ was 7 (IQR 6-12) mm for time-domain dipoles. Current density predicted ictal activity with up to 86 % (60-87 %) accuracy. Distance from time-domain dipoles to the TV was smaller (P = 0.045) in seizure-free (2 [0-4] mm) versus non-seizure-free (12 [2-17] mm) patients, and predicted surgical outcome with 91 % sensitivity and 63 % specificity. Removing near-field data from contacts within the TV negated outcome prediction (P = 0.51).

CONCLUSIONS

Source localization of SEEG accurately mapped ictal onset compared with conventional interpretation. Proximity of dipoles to the TV predicted seizure outcome when near-field recordings were analyzed.

SIGNIFICANCE

Ictal SEEG source localization is useful in corroborating the epileptogenic zone, assuming near-field recordings are obtained.