SleepSEEG: Automatic sleep scoring using intracranial EEG recordings only

Network signatures define consciousness state during focal seizures

OBJECTIVE

Epilepsy is a common neurological disorder affecting 1% of the global population. Loss of consciousness in focal impaired awareness seizures (FIASs) and focal-to-bilateral tonic-clonic seizures (FBTCSs) can be devastating, but the mechanisms are not well understood. Although ictal activity and interictal connectivity changes have been noted, the network states of focal aware seizures (FASs), FIASs, and FBTCSs have not been thoroughly evaluated with network measures ictally.

METHODS

We obtained electrographic data from 74 patients with stereoelectroencephalography (SEEG). Sliding window band power, functional connectivity, and segregation were computed on preictal, ictal, and postictal data. Five-minute epochs of wake, rapid eye movement sleep, and deep sleep were also extracted. Connectivity of subcortical arousal structures was analyzed in a cohort of patients with both SEEG and functional magnetic resonance imaging (fMRI). Given that custom neuromodulation of seizures is predicated on detection of seizure type, a convolutional neural network was used to classify seizure types.

RESULTS

We found that in the frontoparietal association cortex, an area associated with consciousness, both consciousness-impairing seizures (FIASs and FBTCSs) and deep sleep had increases in slow wave delta (1-4 Hz) band power. However, when network measures were employed, we found that only FIASs and deep sleep exhibited an increase in delta segregation and a decrease in gamma segregation. Furthermore, we found that only patients with FIASs had reduced subcortical-to-neocortical functional connectivity with fMRI versus controls. Finally, our deep learning network demonstrated an area under the curve of .75 for detecting consciousness-impairing seizures.

SIGNIFICANCE

This study provides novel insights into ictal network measures in FASs, FIASs, and FBTCSs. Importantly, although both FIASs and FBTCSs result in loss of consciousness, our results suggest that ictal network changes in FIASs uniquely resemble those that occur during deep sleep. Our results may inform novel neuromodulation strategies for preservation of consciousness in epilepsy.

TSFAN: tensorized spatial-frequency attention network with domain adaptation for cross-session EEG-based biometric recognition

Objective.Electroencephalogram (EEG) signals are promising biometrics owning to their invisibility, adapting to the application scenarios with high-security requirements. However, It is challenging to explore EEG identity features without the interference of device and state differences of the subject across sessions. Existing methods treat training sessions as a single domain, affected by the different data distribution among sessions. Although most multi-source unsupervised domain adaptation (MUDA) methods bridge the domain gap between multiple source and target domains individually, relationships among the domain-invariant features of each distribution alignment are neglected.Approach.In this paper, we propose a MUDA method, Tensorized Spatial-Frequency Attention Network (TSFAN), to assist the performance of the target domain for EEG-based biometric recognition. Specifically, significant relationships of domain-invariant features are modeled via a tensorized attention mechanism. It jointly incorporates appropriate common spatial-frequency representations of pairwise source and target but also cross-source domains, without the effect of distribution discrepancy among source domains. Additionally, considering the curse of dimensionality, our TSFAN is approximately represented in Tucker format. Benefiting the low-rank Tucker Network, the TSFAN can scale linearly in the number of domains, providing us the great flexibility to extend TSFAN to the case associated with an arbitrary number of sessions.Main results.Extensive experiments on the representative benchmarks demonstrate the effectiveness of TSFAN in EEG-based biometric recognition, outperforming state-of-the-art approaches, as verified by cross-session validation.Significance.The proposed TSFAN aims to investigate the presence of consistent EEG identity features across sessions. It is achieved by utilizing a novel tensorized attention mechanism that collaborates intra-source transferable information with inter-source interactions, while remaining unaffected by domain shifts in multiple source domains. Furthermore, the electrode selection shows that EEG-based identity features across sessions are distributed across brain regions, and 20 electrodes based on 10-20 standard system are able to extract stable identity information.

Clinical Utility of Sleep Recordings During Presurgical Epilepsy Evaluation With Stereo-Electroencephalography: A Systematic Review

SUMMARY

Although the role of sleep in modulating epileptic activity is well established, many epileptologists overlook the significance of considering sleep during presurgical epilepsy evaluations in cases of drug-resistant epilepsy. Here, we conducted a comprehensive literature review from January 2000 to May 2023 using the PubMed electronic database and compiled evidence to highlight the need to revise the current clinical approach. All articles were assessed for eligibility by two independent reviewers. Our aim was to shed light on the clinical value of incorporating sleep monitoring into presurgical evaluations with stereo-electroencephalography. We present the latest developments on the important bidirectional interactions between sleep and various forms of epileptic activity observed in stereo-electroencephalography recordings. Specifically, epileptic activity is modulated by different sleep stages, peaking in non-rapid eye movement sleep, while being suppressed in rapid eye movement sleep. However, this modulation can vary across different brain regions, underlining the need to account for sleep to accurately pinpoint the epileptogenic zone during presurgical assessments. Finally, we offer practical solutions, such as automated sleep scoring algorithms using stereo-electroencephalography data alone, to seamlessly integrate sleep monitoring into routine clinical practice. It is hoped that this review will provide clinicians with a readily accessible roadmap to the latest evidence concerning the clinical utility of sleep monitoring in the context of stereo-electroencephalography and aid the development of therapeutic and diagnostic strategies to improve patient surgical outcomes.

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.

Utility of intracranial EEG networks depends on re-referencing and connectivity choice

Studies of intracranial EEG networks have been used to reveal seizure generators in patients with drug-resistant epilepsy. Intracranial EEG is implanted to capture the epileptic network, the collection of brain tissue that forms a substrate for seizures to start and spread. Interictal intracranial EEG measures brain activity at baseline, and networks computed during this state can reveal aberrant brain tissue without requiring seizure recordings. Intracranial EEG network analyses require choosing a reference and applying statistical measures of functional connectivity. Approaches to these technical choices vary widely across studies, and the impact of these technical choices on downstream analyses is poorly understood. Our objective was to examine the effects of different re-referencing and connectivity approaches on connectivity results and on the ability to lateralize the seizure onset zone in patients with drug-resistant epilepsy. We applied 48 pre-processing pipelines to a cohort of 125 patients with drug-resistant epilepsy recorded with interictal intracranial EEG across two epilepsy centres to generate intracranial EEG functional connectivity networks. Twenty-four functional connectivity measures across time and frequency domains were applied in combination with common average re-referencing or bipolar re-referencing. We applied an unsupervised clustering algorithm to identify groups of pre-processing pipelines. We subjected each pre-processing approach to three quality tests: (i) the introduction of spurious correlations; (ii) robustness to incomplete spatial sampling; and (iii) the ability to lateralize the clinician-defined seizure onset zone. Three groups of similar pre-processing pipelines emerged: common average re-referencing pipelines, bipolar re-referencing pipelines and relative entropy-based connectivity pipelines. Relative entropy and common average re-referencing networks were more robust to incomplete electrode sampling than bipolar re-referencing and other connectivity methods (Friedman test, Dunn-Šidák test P < 0.0001). Bipolar re-referencing reduced spurious correlations at non-adjacent channels better than common average re-referencing (Δ mean from machine ref = -0.36 versus -0.22) and worse in adjacent channels (Δ mean from machine ref = -0.14 versus -0.40). Relative entropy-based network measures lateralized the seizure onset hemisphere better than other measures in patients with temporal lobe epilepsy (Benjamini-Hochberg-corrected P < 0.05, Cohen's d: 0.60-0.76). Finally, we present an interface where users can rapidly evaluate intracranial EEG pre-processing choices to select the optimal pre-processing methods tailored to specific research questions. The choice of pre-processing methods affects downstream network analyses. Choosing a single method among highly correlated approaches can reduce redundancy in processing. Relative entropy outperforms other connectivity methods in multiple quality tests. We present a method and interface for researchers to optimize their pre-processing methods for deriving intracranial EEG brain networks.

Automated sleep staging on reduced channels in children with epilepsy

Objectives

This study aimed to validate a sleep staging algorithm using in-hospital video-electroencephalogram (EEG) in children without epilepsy, with well-controlled epilepsy (WCE), and with drug-resistant epilepsy (DRE).

Methods

Overnight video-EEG, along with electrooculogram (EOG) and chin electromyogram (EMG), was recorded in children between 4 and 18 years of age. Classical sleep staging was performed manually as a ground truth. An end-to-end hierarchical recurrent neural network for sequence-to-sequence automatic sleep staging (SeqSleepNet) was used to perform automated sleep staging using three channels: C4-A1, EOG, and chin EMG.

Results

In 176 children sleep stages were manually scored: 47 children without epilepsy, 74 with WCE, and 55 with DRE. The 5-class sleep staging accuracy of the automatic sleep staging algorithm was 84.7% for the children without epilepsy, 83.5% for those with WCE, and 80.8% for those with DRE (Kappa of 0.79, 0.77, and 0.73 respectively). Performance per sleep stage was assessed with an F1 score of 0.91 for wake, 0.50 for N1, 0.83 for N2, 0.84 for N3, and 0.86 for rapid eye movement (REM) sleep.

Conclusion

We concluded that the tested algorithm has a high accuracy in children without epilepsy and with WCE. Performance in children with DRE was acceptable, but significantly lower, which could be explained by a tendency of more time spent in N1, and by abundant interictal epileptiform discharges and intellectual disability leading to less recognizable sleep stages. REM sleep time, however, significantly affected in children with DRE, can be detected reliably by the algorithm.Clinical trial registration: ClinicalTrials.gov, identifier NCT04584385.

Interictal intracranial EEG asymmetry lateralizes temporal lobe epilepsy

Patients with drug-resistant temporal lobe epilepsy often undergo intracranial EEG recording to capture multiple seizures in order to lateralize the seizure onset zone. This process is associated with morbidity and often ends in postoperative seizure recurrence. Abundant interictal (between-seizure) data is captured during this process, but these data currently play a small role in surgical planning. Our objective was to predict the laterality of the seizure onset zone using interictal (between-seizure) intracranial EEG data in patients with temporal lobe epilepsy. We performed a retrospective cohort study (single-center study for model development; two-center study for model validation). We studied patients with temporal lobe epilepsy undergoing intracranial EEG at the University of Pennsylvania (internal cohort) and the Medical University of South Carolina (external cohort) between 2015 and 2022. We developed a logistic regression model to predict seizure onset zone laterality using interictal EEG. We compared the concordance between the model-predicted seizure onset zone laterality and the side of surgery between patients with good and poor surgical outcomes. 47 patients (30 women; ages 20-69; 20 left-sided, 10 right-sided, and 17 bilateral seizure onsets) were analyzed for model development and internal validation. 19 patients (10 women; ages 23-73; 5 left-sided, 10 right-sided, 4 bilateral) were analyzed for external validation. The internal cohort cross-validated area under the curve for a model trained using spike rates was 0.83 for a model predicting left-sided seizure onset and 0.68 for a model predicting right-sided seizure onset. Balanced accuracies in the external cohort were 79.3% and 78.9% for the left- and right-sided predictions, respectively. The predicted concordance between the laterality of the seizure onset zone and the side of surgery was higher in patients with good surgical outcome. In conclusion, interictal EEG signatures are distinct across seizure onset zone lateralities. Left-sided seizure onsets are easier to distinguish than right-sided onsets. A model trained on spike rates accurately identifies patients with left-sided seizure onset zones and predicts surgical outcome.

Passive and active markers of cortical excitability in epilepsy

Electroencephalography (EEG) has been the primary diagnostic tool in clinical epilepsy for nearly a century. Its review is performed using qualitative clinical methods that have changed little over time. However, the intersection of higher resolution digital EEG and analytical tools developed in the past decade invites a re-exploration of relevant methodology. In addition to the established spatial and temporal markers of spikes and high-frequency oscillations, novel markers involving advanced postprocessing and active probing of the interictal EEG are gaining ground. This review provides an overview of the EEG-based passive and active markers of cortical excitability in epilepsy and of the techniques developed to facilitate their identification. Several different emerging tools are discussed in the context of specific EEG applications and the barriers we must overcome to translate these tools into clinical practice.

Sleep, oscillations, and epilepsy

Sleep and wake are defined through physiological and behavioral criteria and can be typically separated into non-rapid eye movement (NREM) sleep stages N1, N2, and N3, rapid eye movement (REM) sleep, and wake. Sleep and wake states are not homogenous in time. Their properties vary during the night and day cycle. Given that brain activity changes as a function of NREM, REM, and wake during the night and day cycle, are seizures more likely to occur during NREM, REM, or wake at a specific time? More generally, what is the relationship between sleep-wake cycles and epilepsy? We will review specific examples from clinical data and results from experimental models, focusing on the diversity and heterogeneity of these relationships. We will use a top-down approach, starting with the general architecture of sleep, followed by oscillatory activities, and ending with ionic correlates selected for illustrative purposes, with respect to seizures and interictal spikes. The picture that emerges is that of complexity; sleep disruption and pathological epileptic activities emerge from reorganized circuits. That different circuit alterations can occur across patients and models may explain why sleep alterations and the timing of seizures during the sleep-wake cycle are patient-specific.

Differential cortical network engagement during states of un/consciousness in humans

What happens in the human brain when we are unconscious? Despite substantial work, we are still unsure which brain regions are involved and how they are impacted when consciousness is disrupted. Using intracranial recordings and direct electrical stimulation, we mapped global, network, and regional involvement during wake vs. arousable unconsciousness (sleep) vs. non-arousable unconsciousness (propofol-induced general anesthesia). Information integration and complex processing we`re reduced, while variability increased in any type of unconscious state. These changes were more pronounced during anesthesia than sleep and involved different cortical engagement. During sleep, changes were mostly uniformly distributed across the brain, whereas during anesthesia, the prefrontal cortex was the most disrupted, suggesting that the lack of arousability during anesthesia results not from just altered overall physiology but from a disconnection between the prefrontal and other brain areas. These findings provide direct evidence for different neural dynamics during loss of consciousness compared with loss of arousability.

The Differing Effects of Sleep on Ictal and Interictal Network Dynamics in Drug-Resistant Epilepsy

OBJECTIVE

Sleep has important influences on focal interictal epileptiform discharges (IEDs), and the rates and spatial extent of IEDs are increased in non-rapid eye movement (NREM) sleep. In contrast, the influence of sleep on seizures is less clear, and its effects on seizure topography are poorly documented. We evaluated the influences of NREM sleep on ictal spatiotemporal dynamics and contrasted these with interictal network dynamics.

METHODS

We included patients with drug-resistant focal epilepsy who underwent continuous intracranial electroencephalography (iEEG) with depth electrodes. Patients were selected if they had 1 to 3 seizures from each vigilance state, wakefulness and NREM sleep, within a 48-hour window, and under the same antiseizure medication. A 10-minute epoch of the interictal iEEG was selected per state, and IEDs were detected automatically. A total of 25 patients (13 women; aged 32.5 ± 7.1 years) were included.

RESULTS

The seizure onset pattern, duration, spatiotemporal propagation, and latency of ictal high-frequency activity did not differ significantly between wakefulness and NREM sleep (all p > 0.05). In contrast, IED rates and spatial distribution were increased in NREM compared with wakefulness (p < 0.001, Cliff's d = 0.48 and 0.49). The spatial overlap between vigilance states was higher for seizures (57.1 ± 40.1%) than IEDs (41.7 ± 46.2%; p = 0.001, Cliff's d = 0.51).

INTERPRETATION

In contrast to its effects on IEDs, NREM sleep does not affect ictal spatiotemporal dynamics. This suggests that once the brain surpasses the seizure threshold, it will follow the underlying epileptic network irrespective of the vigilance state. These findings offer valuable insights into neural network dynamics in epilepsy and have important clinical implications for localizing seizure foci. ANN NEUROL 2023.

Unsupervised Multitaper Spectral Method for Identifying REM Sleep in Intracranial EEG Recordings Lacking EOG/EMG Data

A large number of human intracranial EEG (iEEG) recordings have been collected for clinical purposes, in institutions all over the world, but the vast majority of these are unaccompanied by EOG and EMG recordings which are required to separate Wake episodes from REM sleep using accepted methods. In order to make full use of this extremely valuable data, an accurate method of classifying sleep from iEEG recordings alone is required. Existing methods of sleep scoring using only iEEG recordings accurately classify all stages of sleep, with the exception that wake (W) and rapid-eye movement (REM) sleep are not well distinguished. A novel multitaper (Wake vs. REM) alpha-rhythm classifier is developed by generalizing K-means clustering for use with multitaper spectral eigencoefficients. The performance of this unsupervised method is assessed on eight subjects exhibiting normal sleep architecture in a hold-out analysis and is compared against a classical power detector. The proposed multitaper classifier correctly identifies 36±6 min of REM in one night of recorded sleep, while incorrectly labeling less than 10% of all labeled 30 s epochs for all but one subject (human rater reliability is estimated to be near 80%), and outperforms the equivalent statistical-power classical test. Hold-out analysis indicates that when using one night's worth of data, an accurate generalization of the method on new data is likely. For the purpose of studying sleep, the introduced multitaper alpha-rhythm classifier further paves the way to making available a large quantity of otherwise unusable IEEG data.

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.

Regional variability in intracerebral properties of NREM to REM sleep transitions in humans

Transitions between wake and sleep states show a progressive pattern underpinned by local sleep regulation. In contrast, little evidence is available on non-rapid eye movement (NREM) to rapid eye movement (REM) sleep boundaries, considered as mainly reflecting subcortical regulation. Using polysomnography (PSG) combined with stereoelectroencephalography (SEEG) in humans undergoing epilepsy presurgical evaluation, we explored the dynamics of NREM-to-REM transitions. PSG was used to visually score transitions and identify REM sleep features. SEEG-based local transitions were determined automatically with a machine learning algorithm using features validated for automatic intra-cranial sleep scoring (10.5281/zenodo.7410501). We analyzed 2988 channel-transitions from 29 patients. The average transition time from all intracerebral channels to the first visually marked REM sleep epoch was 8 s ± 1 min 58 s, with a great heterogeneity between brain areas. Transitions were observed first in the lateral occipital cortex, preceding scalp transition by 1 min 57 s ± 2 min 14 s (d = -0.83), and close to the first sawtooth wave marker. Regions with late transitions were the inferior frontal and orbital gyri (1 min 1 s ± 2 min 1 s, d = 0.43, and 1 min 1 s ± 2 min 5 s, d = 0.43, after scalp transition). Intracranial transitions were earlier than scalp transitions as the night advanced (last sleep cycle, d = -0.81). We show a reproducible gradual pattern of REM sleep initiation, suggesting the involvement of cortical mechanisms of regulation. This provides clues for understanding oneiric experiences occurring at the NREM/REM boundary.

Spectro-spatial features in distributed human intracranial activity proactively encode peripheral metabolic activity

Mounting evidence demonstrates that the central nervous system (CNS) orchestrates glucose homeostasis by sensing glucose and modulating peripheral metabolism. Glucose responsive neuronal populations have been identified in the hypothalamus and several corticolimbic regions. However, how these CNS gluco-regulatory regions modulate peripheral glucose levels is not well understood. To better understand this process, we simultaneously measured interstitial glucose concentrations and local field potentials in 3 human subjects from cortical and subcortical regions, including the hypothalamus in one subject. Correlations between high frequency activity (HFA, 70-170 Hz) and peripheral glucose levels are found across multiple brain regions, notably in the hypothalamus, with correlation magnitude modulated by sleep-wake cycles, circadian coupling, and hypothalamic connectivity. Correlations are further present between non-circadian (ultradian) HFA and glucose levels which are higher during awake periods. Spectro-spatial features of neural activity enable decoding of peripheral glucose levels both in the present and up to hours in the future. Our findings demonstrate proactive encoding of homeostatic glucose dynamics by the CNS.

Spatial and temporal correlations in human cortex are inherently linked and predicted by functional hierarchy, vigilance state as well as antiepileptic drug load

The ability of neural circuits to integrate information over time and across different cortical areas is believed an essential ingredient for information processing in the brain. Temporal and spatial correlations in cortex dynamics have independently been shown to capture these integration properties in task-dependent ways. A fundamental question remains if temporal and spatial integration properties are linked and what internal and external factors shape these correlations. Previous research on spatio-temporal correlations has been limited in duration and coverage, thus providing only an incomplete picture of their interdependence and variability. Here, we use long-term invasive EEG data to comprehensively map temporal and spatial correlations according to cortical topography, vigilance state and drug dependence over extended periods of time. We show that temporal and spatial correlations in cortical networks are intimately linked, decline under antiepileptic drug action, and break down during slow-wave sleep. Further, we report temporal correlations in human electrophysiology signals to increase with the functional hierarchy in cortex. Systematic investigation of a neural network model suggests that these dynamical features may arise when dynamics are poised near a critical point. Our results provide mechanistic and functional links between specific measurable changes in the network dynamics relevant for characterizing the brain's changing information processing capabilities.

Integration of white matter architecture to stereo-EEG better describes epileptic spike propagation

OBJECTIVE

Stereo-electroencephalography (SEEG)-derived epilepsy networks are used to better understand a patient's epilepsy; however, a unimodal approach provides an incomplete picture. We combine tractography and SEEG to determine the relationship between spike propagation and the white matter architecture and to improve our understanding of spike propagation mechanisms.

METHODS

Probablistic tractography from diffusion imaging (dMRI) of matched subjects from the Human Connectome Project (HCP) was combined with patient-specific SEEG-derived spike propagation networks. Two regions-of-interest (ROIs) with a significant spike propagation relationship constituted a Propagation Pair.

RESULTS

In 56 of 59 patients, Propagation Pairs were more often tract-connected as compared to all ROI pairs (p < 0.01; d = -1.91). The degree of spike propagation between tract-connected ROIs was greater (39 ± 21%) compared to tract-unconnected ROIs (31 ± 18%; p < 0.0001). Within the same network, ROIs receiving propagation earlier were more often tract-connected to the source (59.7%) as compared to late receivers (25.4%; p < 0.0001).

CONCLUSIONS

Brain regions involved in spike propagation are more likely to be connected by white matter tracts. Between nodes, presence of tracts suggests a direct course of propagation, whereas the absence of tracts suggests an indirect course of propagation.

SIGNIFICANCE

We demonstrate a logical and consistent relationship between spike propagation and the white matter architecture.