Electrocorticography and stereo EEG provide distinct measures of brain connectivity: implications for network models

Supplementary motor area in speech initiation: A large-scale intracranial EEG evaluation of stereotyped word articulation

Speech production engages a distributed network of cortical and subcortical brain regions. The supplementary motor area (SMA) has long been thought to be a key hub in coordinating across these regions to initiate voluntary movements, including speech. We analyzed direct intracranial recordings from 115 patients with epilepsy as they articulated a single word in a subset of trials from a picture-naming task. We aimed to characterize the temporal dynamics of SMA relative to other cortical regions. SMA and preSMA were among the first regions to activate after cue onset, peaked in activity before articulation onset, and were the earliest regions to predict trial-to-trial response time. Neural activity at single electrodes in SMA and preSMA was closely associated with speech initiation; activity began at a highly predictable time after stimulus onset and extended until speech onset for any given trial. Our results support the idea that SMA is a key node in the speech initiation network.

Localizing the seizure onset zone and predicting the surgery outcomes in patients with drug-resistant epilepsy: A new approach based on the causal network

BACKGROUND AND OBJECTIVE

Accurate localization of the seizure onset zone (SOZ) is crucial for surgical treatment in patients with drug-resistant epilepsy (DRE). However, clinical identification of SOZ often relies on physician experience and has a certain subjectivity. Therefore, it is emergent to develop quantitative computational tools to assist clinicians in identifying SOZ.

METHODS

We conduct a retrospective study on intracranial electroencephalography (iEEG) data from 46 patients with DRE. The interactions between different brain regions are quantified by using the phase transfer entropy (PTE), based on which the causal influence index (CII) is proposed to quantify the degree of influence of nodes on the network. Subsequently, the features extracted by the CII are used to construct a random forest classification model, which the performance in identifying SOZ and the generalizability are validated in patients with successful surgeries. Then, based on the CII features of the clinically labeled SOZ, a logistic regression prediction model is constructed to predict the probability of surgical success. The statistical analysis between patients with successful and failed surgery is conducted with the Mann-Whitney U test. Finally, the consistency between the predicted SOZ and the clinically labeled SOZ is verified across different Engel classes.

RESULTS

The classification model combining the low-frequency and high-frequency features can achieve an accuracy of 82.18% (sensitivity: 85.01%, specificity: 79.69%) and an area under curve (AUC) of 0.90 in identifying SOZ. Furthermore, the model exhibits strong generalizability in identifying SOZ in patients with MRI lesional and non-lesional, as well as those implanted with electrocorticography (ECOG) and stereotactic EEG (SEEG) electrodes. Moreover, the prediction model could achieve an average accuracy of 79.8% and an AUC of 0.84. Of note, the prediction of surgical success probability is significant between patients with successful and failed surgeries (P<0.001). Correspondingly, the highest consistency between model-predicted SOZ and clinically labeled SOZ can be observed in patients with successful surgeries, but this consistency gradually decreases with increasing Engel classes.

CONCLUSIONS

These results demonstrate that the CII may be a potential biomarker for identifying the SOZ in patients with DRE, which may provide a new perspective for the treatment of epilepsy.

A Multivariate Approach to Quantifying Risk Factors Impacting Stereotactic Robotic-Guided Stereoelectroencephalography

BACKGROUND AND OBJECTIVES

Stereoelectroencephalography (SEEG) is an important method for invasive monitoring to establish surgical candidacy in approximately half of refractory epilepsy patients. Identifying factors affecting lead placement can mitigate potential surgical risks. This study applies multivariate analyses to identify perioperative factors affecting stereotactic electrode placement.

METHODS

We collected registration and accuracy data for consecutive patients undergoing SEEG implantation between May 2022 and November 2023. Stereotactic robotic guidance, using intraoperative imaging and a novel frame-based fiducial, was used for planning and SEEG implantation. Entry-point (EE), target-point (TE), and angular errors were measured, and statistical univariate and multivariate linear regression analyses were performed.

RESULTS

Twenty-seven refractory epilepsy patients (aged 15-57 years) undergoing SEEG were reviewed. Sixteen patients had unilateral implantation (10 left-sided, 6 right-sided); 11 patients underwent bilateral implantation. The mean number of electrodes per patient was 18 (SD = 3) with an average registration mean error of 0.768 mm (SD = 0.108). Overall, 486 electrodes were reviewed. Univariate analysis showed significant correlations of lead error with skull thickness (EE: P = .003; TE: P = .012); entry angle (EE: P < .001; TE: P < .001; angular error: P = .030); lead length (TE: P = .020); and order of electrode implantation (EE: P = .003; TE: P = .001). Three multiple linear regression models were used. All models featured predictors of implantation region (157 temporal, 241 frontal, 79 parietal, 9 occipital); skull thickness (mean = 5.80 mm, SD = 2.97 mm); order (range: 1-23); and entry angle in degrees (mean = 75.47, SD = 11.66). EE and TE error models additionally incorporated lead length (mean = 44.08 mm, SD = 13.90 mm) as a predictor. Implantation region and entry angle were significant predictors of error (P ≤ .05).

CONCLUSION

Our study identified 2 primary predictors of SEEG lead error, region of implantation and entry angle, with nonsignificant contributions from lead length or order of electrode placement. Future considerations for SEEG may consider varying regional approaches and angles for more optimal accuracy in lead placement.

Generalisability of epileptiform patterns across time and patients

The complexity of localising the epileptogenic zone (EZ) contributes to surgical resection failures in achieving seizure freedom. The distinct patterns of epileptiform activity during interictal and ictal phases, varying across patients, often lead to suboptimal localisation using electroencephalography (EEG) features. We posed two key questions: whether neural signals reflecting epileptogenicity generalise from interictal to ictal time windows within each patient, and whether epileptiform patterns generalise across patients. Utilising an intracranial EEG dataset from 55 patients, we extracted a large battery of simple to complex features from stereo-EEG (SEEG) and electrocorticographic (ECoG) neural signals during interictal and ictal windows. Our features (n = 34) quantified many aspects of the signals including statistical moments, complexities, frequency-domain and cross-channel network attributes. Decision tree classifiers were then trained and tested on distinct time windows and patients to evaluate the generalisability of epileptogenic patterns across time and patients, respectively. Evidence strongly supported generalisability from interictal to ictal time windows across patients, particularly in signal power and high-frequency network-based features. Consistent patterns of epileptogenicity were observed across time windows within most patients, and signal features of epileptogenic regions generalised across patients, with higher generalisability in the ictal window. Signal complexity features were particularly contributory in cross-patient generalisation across patients. These findings offer insights into generalisable features of epileptic neural activity across time and patients, with implications for future automated approaches to supplement other EZ localisation methods.

The time-evolving epileptic brain network: concepts, definitions, accomplishments, perspectives

Epilepsy is now considered a network disease that affects the brain across multiple levels of spatial and temporal scales. The paradigm shift from an epileptic focus-a discrete cortical area from which seizures originate-to a widespread epileptic network-spanning lobes and hemispheres-considerably advanced our understanding of epilepsy and continues to influence both research and clinical treatment of this multi-faceted high-impact neurological disorder. The epileptic network, however, is not static but evolves in time which requires novel approaches for an in-depth characterization. In this review, we discuss conceptual basics of network theory and critically examine state-of-the-art recording techniques and analysis tools used to assess and characterize a time-evolving human epileptic brain network. We give an account on current shortcomings and highlight potential developments towards an improved clinical management of epilepsy.

From basic sciences and engineering to epileptology: A translational approach

Collaborative efforts between basic scientists, engineers, and clinicians are enabling translational epileptology. In this article, we summarize the recent advances presented at the International Conference for Technology and Analysis of Seizures (ICTALS 2022): (1) novel developments of structural magnetic resonance imaging; (2) latest electroencephalography signal-processing applications; (3) big data for the development of clinical tools; (4) the emerging field of hyperdimensional computing; (5) the new generation of artificial intelligence (AI)-enabled neuroprostheses; and (6) the use of collaborative platforms to facilitate epilepsy research translation. We highlight the promise of AI reported in recent investigations and the need for multicenter data-sharing initiatives.

Physiological and pathological neuronal connectivity in the living human brain based on intracranial EEG signals: the current state of research

Over the past decades, studies of human brain networks have received growing attention as the assessment and modelling of connectivity in the brain is a topic of high impact with potential application in the understanding of human brain organization under both physiological as well as various pathological conditions. Under specific diagnostic settings, human neuronal signal can be obtained from intracranial EEG (iEEG) recording in epilepsy patients that allows gaining insight into the functional organisation of living human brain. There are two approaches to assess brain connectivity in the iEEG-based signal: evaluation of spontaneous neuronal oscillations during ongoing physiological and pathological brain activity, and analysis of the electrophysiological cortico-cortical neuronal responses, evoked by single pulse electrical stimulation (SPES). Both methods have their own advantages and limitations. The paper outlines available methodological approaches and provides an overview of current findings in studies of physiological and pathological human brain networks, based on intracranial EEG recordings.

Personalized dynamic network models of the human brain as a future tool for planning and optimizing epilepsy therapy

Epilepsy is a common neurological disorder, with one third of patients not responding to currently available antiepileptic drugs. The proportion of pharmacoresistant epilepsies has remained unchanged for many decades. To cure epilepsy and control seizures requires a paradigm shift in the development of new approaches to epilepsy diagnosis and treatment. Contemporary medicine has benefited from the exponential growth of computational modeling, and the application of network dynamics theory to understanding and treating human brain disorders. In epilepsy, the introduction of these approaches has led to personalized epileptic network modeling that can explore the patient's seizure genesis and predict the functional impact of resection on its individual network's propensity to seize. The application of the dynamic systems approach to neurostimulation therapy of epilepsy allows designing stimulation strategies that consider the patient's seizure dynamics and long-term fluctuations in the stability of their epileptic networks. In this article, we review, in a nontechnical fashion suitable for a broad neuroscientific audience, recent progress in personalized dynamic brain network modeling that is shaping the future approach to the diagnosis and treatment of epilepsy.

The seizure severity score: a quantitative tool for comparing seizures and their response to therapy

Objective.Epilepsy is a neurological disorder characterized by recurrent seizures which vary widely in severity, from clinically silent to prolonged convulsions. Measuring severity is crucial for guiding therapy, particularly when complete control is not possible. Seizure diaries, the current standard for guiding therapy, are insensitive to the duration of events or the propagation of seizure activity across the brain. We present a quantitative seizure severity score that incorporates electroencephalography (EEG) and clinical data and demonstrate how it can guide epilepsy therapies.Approach.We collected intracranial EEG and clinical semiology data from 54 epilepsy patients who had 256 seizures during invasive, in-hospital presurgical evaluation. We applied an absolute slope algorithm to EEG recordings to identify seizing channels. From this data, we developed a seizure severity score that combines seizure duration, spread, and semiology using non-negative matrix factorization. For validation, we assessed its correlation with independent measures of epilepsy burden: seizure types, epilepsy duration, a pharmacokinetic model of medication load, and response to epilepsy surgery. We investigated the association between the seizure severity score and preictal network features.Main results.The seizure severity score augmented clinical classification by objectively delineating seizure duration and spread from recordings in available electrodes. Lower preictal medication loads were associated with higher seizure severity scores (p= 0.018, 97.5% confidence interval = [-1.242, -0.116]) and lower pre-surgical severity was associated with better surgical outcome (p= 0.042). In 85% of patients with multiple seizure types, greater preictal change from baseline was associated with higher severity.Significance.We present a quantitative measure of seizure severity that includes EEG and clinical features, validated on gold standard in-patient recordings. We provide a framework for extending our tool's utility to ambulatory EEG devices, for linking it to seizure semiology measured by wearable sensors, and as a tool to advance data-driven epilepsy care.

Functional connectivity discriminates epileptogenic states and predicts surgical outcome in children with drug resistant epilepsy

Normal brain functioning emerges from a complex interplay among regions forming networks. In epilepsy, these networks are disrupted causing seizures. Highly connected nodes in these networks are epilepsy surgery targets. Here, we assess whether functional connectivity (FC) using intracranial electroencephalography can quantify brain regions epileptogenicity and predict surgical outcome in children with drug resistant epilepsy (DRE). We computed FC between electrodes on different states (i.e. interictal without spikes, interictal with spikes, pre-ictal, ictal, and post-ictal) and frequency bands. We then estimated the electrodes' nodal strength. We compared nodal strength between states, inside and outside resection for good- (n = 22, Engel I) and poor-outcome (n = 9, Engel II-IV) patients, respectively, and tested their utility to predict the epileptogenic zone and outcome. We observed a hierarchical epileptogenic organization among states for nodal strength: lower FC during interictal and pre-ictal states followed by higher FC during ictal and post-ictal states (p < 0.05). We further observed higher FC inside resection (p < 0.05) for good-outcome patients on different states and bands, and no differences for poor-outcome patients. Resection of nodes with high FC was predictive of outcome (positive and negative predictive values: 47-100%). Our findings suggest that FC can discriminate epileptogenic states and predict outcome in patients with DRE.

Quantitative approaches to guide epilepsy surgery from intracranial EEG

Over the past 10 years, the drive to improve outcomes from epilepsy surgery has stimulated widespread interest in methods to quantitatively guide epilepsy surgery from intracranial EEG (iEEG). Many patients fail to achieve seizure freedom, in part due to the challenges in subjective iEEG interpretation. To address this clinical need, quantitative iEEG analytics have been developed using a variety of approaches, spanning studies of seizures, interictal periods, and their transitions, and encompass a range of techniques including electrographic signal analysis, dynamical systems modeling, machine learning and graph theory. Unfortunately, many methods fail to generalize to new data and are sensitive to differences in pathology and electrode placement. Here, we critically review selected literature on computational methods of identifying the epileptogenic zone from iEEG. We highlight shared methodological challenges common to many studies in this field and propose ways that they can be addressed. One fundamental common pitfall is a lack of open-source, high-quality data, which we specifically address by sharing a centralized high-quality, well-annotated, multicentre dataset consisting of >100 patients to support larger and more rigorous studies. Ultimately, we provide a road map to help these tools reach clinical trials and hope to improve the lives of future patients.

Quantifying trial-by-trial variability during cortico-cortical evoked potential mapping of epileptogenic tissue

OBJECTIVE

Measuring cortico-cortical evoked potentials (CCEPs) is a promising tool for mapping epileptic networks, but it is not known how variability in brain state and stimulation technique might impact the use of CCEPs for epilepsy localization. We test the hypotheses that (1) CCEPs demonstrate systematic variability across trials and (2) CCEP amplitudes depend on the timing of stimulation with respect to endogenous, low-frequency oscillations.

METHODS

We studied 11 patients who underwent CCEP mapping after stereo-electroencephalography electrode implantation for surgical evaluation of drug-resistant epilepsy. Evoked potentials were measured from all electrodes after each pulse of a 30 s, 1 Hz bipolar stimulation train. We quantified monotonic trends, phase dependence, and standard deviation (SD) of N1 (15-50 ms post-stimulation) and N2 (50-300 ms post-stimulation) amplitudes across the 30 stimulation trials for each patient. We used linear regression to quantify the relationship between measures of CCEP variability and the clinical seizure-onset zone (SOZ) or interictal spike rates.

RESULTS

We found that N1 and N2 waveforms exhibited both positive and negative monotonic trends in amplitude across trials. SOZ electrodes and electrodes with high interictal spike rates had lower N1 and N2 amplitudes with higher SD across trials. Monotonic trends of N1 and N2 amplitude were more positive when stimulating from an area with higher interictal spike rate. We also found intermittent synchronization of trial-level N1 amplitude with low-frequency phase in the hippocampus, which did not localize the SOZ.

SIGNIFICANCE

These findings suggest that standard approaches for CCEP mapping, which involve computing a trial-averaged response over a .2-1 Hz stimulation train, may be masking inter-trial variability that localizes to epileptogenic tissue. We also found that CCEP N1 amplitudes synchronize with ongoing low-frequency oscillations in the hippocampus. Further targeted experiments are needed to determine whether phase-locked stimulation could have a role in localizing epileptogenic tissue.

MEG abnormalities and mechanisms of surgical failure in neocortical epilepsy

OBJECTIVE

Epilepsy surgery fails to achieve seizure freedom in 30%-40% of cases. It is not fully understood why some surgeries are unsuccessful. By comparing interictal magnetoencephalography (MEG) band power from patient data to normative maps, which describe healthy spatial and population variability, we identify patient-specific abnormalities relating to surgical failure. We propose three mechanisms contributing to poor surgical outcome: (1) not resecting the epileptogenic abnormalities (mislocalization), (2) failing to remove all epileptogenic abnormalities (partial resection), and (3) insufficiently impacting the overall cortical abnormality. Herein we develop markers of these mechanisms, validating them against patient outcomes.

METHODS

Resting-state MEG recordings were acquired for 70 healthy controls and 32 patients with refractory neocortical epilepsy. Relative band-power spatial maps were computed using source-localized recordings. Patient and region-specific band-power abnormalities were estimated as the maximum absolute z-score across five frequency bands using healthy data as a baseline. Resected regions were identified using postoperative magnetic resonance imaging (MRI). We hypothesized that our mechanistically interpretable markers would discriminate patients with and without postoperative seizure freedom.

RESULTS

Our markers discriminated surgical outcome groups (abnormalities not targeted: area under the curve [AUC] = 0.80, p = .003; partial resection of epileptogenic zone: AUC = 0.68, p = .053; and insufficient cortical abnormality impact: AUC = 0.64, p = .096). Furthermore, 95% of those patients who were not seizure-free had markers of surgical failure for at least one of the three proposed mechanisms. In contrast, of those patients without markers for any mechanism, 80% were ultimately seizure-free.

SIGNIFICANCE

The mapping of abnormalities across the brain is important for a wide range of neurological conditions. Here we have demonstrated that interictal MEG band-power mapping has merit for the localization of pathology and improving our mechanistic understanding of epilepsy. Our markers for mechanisms of surgical failure could be used in the future to construct predictive models of surgical outcome, aiding clinical teams during patient pre-surgical evaluations.

Addressing spatial bias in intracranial EEG functional connectivity analyses for epilepsy surgical planning

Objective.To determine the effect of epilepsy on intracranial electroencephalography (EEG) functional connectivity, and the ability of functional connectivity to localize the seizure onset zone (SOZ), controlling for spatial biases.Approach.We analyzed intracranial EEG data from patients with drug-resistant epilepsy admitted for pre-surgical planning. We calculated intracranial EEG functional networks and determined whether changes in functional connectivity lateralized the SOZ using a spatial subsampling method to control for spatial bias. We developed a 'spatial null model' to localize the SOZ electrode using only spatial sampling information, ignoring EEG data. We compared the performance of this spatial null model against models incorporating EEG functional connectivity and interictal spike rates.Main results.About 110 patients were included in the study, although the number of patients differed across analyses. Controlling for spatial sampling, the average connectivity was lower in the SOZ region relative to the same anatomic region in the contralateral hemisphere. A model using intra-hemispheric connectivity accurately lateralized the SOZ (average accuracy 75.5%). A spatial null model incorporating spatial sampling information alone achieved moderate accuracy in classifying SOZ electrodes (mean AUC = 0.70, 95% CI 0.63-0.77). A model incorporating intracranial EEG functional connectivity and spike rate data further outperformed this spatial null model (AUC 0.78,p= 0.002 compared to spatial null model). However, a model incorporating functional connectivity without spike rate data did not significantly outperform the null model (AUC 0.72,p= 0.38).Significance.Intracranial EEG functional connectivity is reduced in the SOZ region, and interictal data predict SOZ electrode localization and laterality, however a predictive model incorporating functional connectivity without interictal spike rates did not significantly outperform a spatial null model. We propose constructing a spatial null model to provide an estimate of the pre-implant hypothesis of the SOZ, and to serve as a benchmark for further machine learning algorithms in order to avoid overestimating model performance because of electrode sampling alone.

Normative intracranial EEG maps epileptogenic tissues in focal epilepsy

Planning surgery for patients with medically refractory epilepsy often requires recording seizures using intracranial EEG. Quantitative measures derived from interictal intracranial EEG yield potentially appealing biomarkers to guide these surgical procedures; however, their utility is limited by the sparsity of electrode implantation as well as the normal confounds of spatiotemporally varying neural activity and connectivity. We propose that comparing intracranial EEG recordings to a normative atlas of intracranial EEG activity and connectivity can reliably map abnormal regions, identify targets for invasive treatment and increase our understanding of human epilepsy. Merging data from the Penn Epilepsy Center and a public database from the Montreal Neurological Institute, we aggregated interictal intracranial EEG retrospectively across 166 subjects comprising >5000 channels. For each channel, we calculated the normalized spectral power and coherence in each canonical frequency band. We constructed an intracranial EEG atlas by mapping the distribution of each feature across the brain and tested the atlas against data from novel patients by generating a z-score for each channel. We demonstrate that for seizure onset zones within the mesial temporal lobe, measures of connectivity abnormality provide greater distinguishing value than univariate measures of abnormal neural activity. We also find that patients with a longer diagnosis of epilepsy have greater abnormalities in connectivity. By integrating measures of both single-channel activity and inter-regional functional connectivity, we find a better accuracy in predicting the seizure onset zones versus normal brain (area under the curve = 0.77) compared with either group of features alone. We propose that aggregating normative intracranial EEG data across epilepsy centres into a normative atlas provides a rigorous, quantitative method to map epileptic networks and guide invasive therapy. We publicly share our data, infrastructure and methods, and propose an international framework for leveraging big data in surgical planning for refractory epilepsy.

Probabilistic comparison of gray and white matter coverage between depth and surface intracranial electrodes in epilepsy

In this study, we quantified the coverage of gray and white matter during intracranial electroencephalography in a cohort of epilepsy patients with surface and depth electrodes. We included 65 patients with strip electrodes (n = 12), strip and grid electrodes (n = 24), strip, grid, and depth electrodes (n = 7), or depth electrodes only (n = 22). Patient-specific imaging was used to generate probabilistic gray and white matter maps and atlas segmentations. Gray and white matter coverage was quantified using spherical volumes centered on electrode centroids, with radii ranging from 1 to 15 mm, along with detailed finite element models of local electric fields. Gray matter coverage was highly dependent on the chosen radius of influence (RoI). Using a 2.5 mm RoI, depth electrodes covered more gray matter than surface electrodes; however, surface electrodes covered more gray matter at RoI larger than 4 mm. White matter coverage and amygdala and hippocampal coverage was greatest for depth electrodes at all RoIs. This study provides the first probabilistic analysis to quantify coverage for different intracranial recording configurations. Depth electrodes offer increased coverage of gray matter over other recording strategies if the desired signals are local, while subdural grids and strips sample more gray matter if the desired signals are diffuse.