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-Sidá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.

Alcohol for seizure induction in the epilepsy monitoring unit.

RATIONALE: Seizure induction techniques are used in the epilepsy monitoring unit (EMU) to increase diagnostic yield and reduce length of stay. There are insufficient data on the efficacy of alcohol as an induction technique. METHODS: We performed a retrospective cohort study using six years of EMU data at our institution. We compared cases who received alcohol for seizure induction to matched controls who did not. The groups were matched on the following variables: age, reason for admission, length of stay, number of antiseizure medications (ASM) at admission, whether ASMs were tapered during admission, and presence of interictal epileptiform discharges. We used both propensity score and exact matching strategies. We compared the likelihood of epileptic seizures and nonepileptic events in cases versus controls using Kaplan-Meier time-to-event analysis, as well as odds ratios for these outcomes occurring at any time during the admission. RESULTS: We analyzed 256 cases who received alcohol (median dose 2.5 standard drinks) and 256 propensity score-matched controls. Cases who received alcohol were no more likely than controls to have an epileptic seizure (X2(1) = 0.01, p = 0.93) or nonepileptic event (X2(1) = 2.1, p = 0.14) in the first 48 h after alcohol administration. For the admission overall, cases were no more likely to have an epileptic seizure (OR 0.89, 95 % CI 0.61-1.28, p = 0.58), nonepileptic event (OR 0.97, CI 0.62-1.53, p = 1.00), nor require rescue benzodiazepine (OR 0.63, CI 0.35-1.12, p = 0.15). Stratified analyses revealed no increased risk of epileptic seizure in any subgroups. Sensitivity analysis using exact matching showed that results were robust to matching strategy. CONCLUSIONS: Alcohol was not an effective induction technique in the EMU. This finding has implications for counseling patients with epilepsy about the risks of drinking alcohol in moderation in their daily lives.

Resting-state background features demonstrate multidien cycles in long-term EEG device recordings.

BACKGROUND: Longitudinal EEG recorded by implanted devices is critical for understanding and managing epilepsy. Recent research reports patient-specific, multi-day cycles in device-detected epileptiform events that coincide with increased likelihood of clinical seizures. Understanding these cycles could elucidate mechanisms generating seizures and advance drug and neurostimulation therapies. OBJECTIVE/HYPOTHESIS: We hypothesize that seizure-correlated cycles are present in background neural activity, independent of interictal epileptiform spikes, and that neurostimulation may temporarily interrupt these cycles. METHODS: We analyzed regularly-recorded seizure-free data epochs from 20 patients implanted with a responsive neurostimulation (RNS) device for at least 1.5 years, to explore the relationship between cycles in device-detected interictal epileptiform activity (dIEA), clinician-validated interictal spikes, background EEG features, and neurostimulation. RESULTS: Background EEG features tracked the cycle phase of dIEA in all patients (AUC: 0.63 [0.56-0.67]) with a greater effect size compared to clinically annotated spike rate alone (AUC: 0.55 [0.53-0.61], p < 0.01). After accounting for circadian variation and spike rate, we observed significant population trends in elevated theta and beta band power and theta and alpha connectivity features at the cycle peaks (sign test, p < 0.05). In the period directly after stimulation we observe a decreased association between cycle phase and EEG features compared to background recordings (AUC: 0.58 [0.55-0.64]). CONCLUSIONS: Our findings suggest that seizure-correlated dIEA cycles are not solely due to epileptiform discharges but are associated with background measures of brain state; and that neurostimulation may temporarily interrupt these cycles. These results may help elucidate mechanisms underlying seizure generation, provide new biomarkers for seizure risk, and facilitate monitoring, treating, and managing epilepsy with implantable devices.

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.

Quantifying interictal intracranial EEG to predict focal epilepsy.

Introduction: Intracranial EEG (IEEG) is used for 2 main purposes, to determine: (1) if epileptic networks are amenable to focal treatment and (2) where to intervene. Currently these questions are answered qualitatively and sometimes differently across centers. There is a need for objective, standardized methods to guide surgical decision making and to enable large scale data analysis across centers and prospective clinical trials. Methods: We analyzed interictal data from 101 patients with drug resistant epilepsy who underwent presurgical evaluation with IEEG. We chose interictal data because of its potential to reduce the morbidity and cost associated with ictal recording. 65 patients had unifocal seizure onset on IEEG, and 36 were non-focal or multi-focal. We quantified the spatial dispersion of implanted electrodes and interictal IEEG abnormalities for each patient. We compared these measures against the "5 Sense Score (5SS)," a pre-implant estimate of the likelihood of focal seizure onset, and assessed their ability to predict the clinicians' choice of therapeutic intervention and the patient outcome. Results: The spatial dispersion of IEEG electrodes predicted network focality with precision similar to the 5SS (AUC = 0.67), indicating that electrode placement accurately reflected pre-implant information. A cross-validated model combining the 5SS and the spatial dispersion of interictal IEEG abnormalities significantly improved this prediction (AUC = 0.79; p<0.05). The combined model predicted ultimate treatment strategy (surgery vs. device) with an AUC of 0.81 and post-surgical outcome at 2 years with an AUC of 0.70. The 5SS, interictal IEEG, and electrode placement were not correlated and provided complementary information. Conclusions: Quantitative, interictal IEEG significantly improved upon pre-implant estimates of network focality and predicted treatment with precision approaching that of clinical experts. We present this study as an important step in building standardized, quantitative tools to guide epilepsy surgery.

Resting-state background features demonstrate multidien cycles in long-term EEG device recordings.

Background: Longitudinal EEG recorded by implanted devices is critical for understanding and managing epilepsy. Recent research reports patient-specific, multi-day cycles in device-detected epileptiform events that coincide with increased likelihood of clinical seizures. Understanding these cycles could elucidate mechanisms generating seizures and advance drug and neurostimulation therapies. Objective/Hypothesis: We hypothesize that seizure-correlated cycles are present in background neural activity, independent of interictal epileptiform spikes, and that neurostimulation may disrupt these cycles. Methods: We analyzed regularly-recorded seizure-free data epochs from 20 patients implanted with a responsive neurostimulation (RNS) device for at least 1.5 years, to explore the relationship between cycles in device-detected interictal epileptiform activity (dIEA), clinician-validated interictal spikes, background EEG features, and neurostimulation. Results: Background EEG features tracked the cycle phase of dIEA in all patients (AUC: 0.63 [0.56 - 0.67]) with a greater effect size compared to clinically annotated spike rate alone (AUC: 0.55 [0.53-0.61], p < 0.01). After accounting for circadian variation and spike rate, we observed significant population trends in elevated theta and beta band power and theta and alpha connectivity features at the cycle peaks (sign test, p < 0.05). In the period directly after stimulation we observe a decreased association between cycle phase and EEG features compared to background recordings (AUC: 0.58 [0.55-0.64]). Conclusions: Our findings suggest that seizure-correlated dIEA cycles are not solely due to epileptiform discharges but are associated with background measures of brain state; and that neurostimulation may disrupt these cycles. These results may help elucidate mechanisms underlying seizure generation, provide new biomarkers for seizure risk, and facilitate monitoring, treating, and managing epilepsy with implantable devices.

A comprehensive wireless neurological and cardiopulmonary monitoring platform for pediatrics.

Neurodevelopment in the first 10 years of life is a critical time window during which milestones that define an individual's functional potential are achieved. Comprehensive multimodal neurodevelopmental monitoring is particularly crucial for socioeconomically disadvantaged, marginalized, historically underserved and underrepresented communities as well as medically underserved areas. Solutions designed for use outside the traditional clinical environment represent an opportunity for addressing such health inequalities. In this work, we present an experimental platform, ANNE EEG, which adds 16-channel cerebral activity monitoring to the existing, USA FDA-cleared ANNE wireless monitoring platform which provides continuous electrocardiography, respiratory rate, pulse oximetry, motion, and temperature measurements. The system features low-cost consumables, real-time control and streaming with widely available mobile devices, and fully wearable operation to allow a child to remain in their naturalistic environment. This multi-center pilot study successfully collected ANNE EEG recordings from 91 neonatal and pediatric patients at academic quaternary pediatric care centers and in LMIC settings. We demonstrate the practicality and feasibility to conduct electroencephalography studies with high levels of accuracy, validated via both quantitative and qualitative metrics, compared against gold standard systems. An overwhelming majority of parents surveyed during studies indicated not only an overall preference for the wireless system, but also that its use would improve their children's physical and emotional health. Our findings demonstrate the potential for the ANNE system to perform multimodal monitoring to screen for a variety of neurologic diseases that have the potential to negatively impact neurodevelopment.

Long-term epilepsy outcome dynamics revealed by natural language processing of clinic notes.

OBJECTIVE: Electronic medical records allow for retrospective clinical research with large patient cohorts. However, epilepsy outcomes are often contained in free text notes that are difficult to mine. We recently developed and validated novel natural language processing (NLP) algorithms to automatically extract key epilepsy outcome measures from clinic notes. In this study, we assessed the feasibility of extracting these measures to study the natural history of epilepsy at our center. METHODS: We applied our previously validated NLP algorithms to extract seizure freedom, seizure frequency, and date of most recent seizure from outpatient visits at our epilepsy center from 2010 to 2022. We examined the dynamics of seizure outcomes over time using Markov model-based probability and Kaplan-Meier analyses. RESULTS: Performance of our algorithms on classifying seizure freedom was comparable to that of human reviewers (algorithm F1 = .88 vs. human annotator κ = .86). We extracted seizure outcome data from 55 630 clinic notes from 9510 unique patients written by 53 unique authors. Of these, 30% were classified as seizure-free since the last visit, 48% of non-seizure-free visits contained a quantifiable seizure frequency, and 47% of all visits contained the date of most recent seizure occurrence. Among patients with at least five visits, the probabilities of seizure freedom at the next visit ranged from 12% to 80% in patients having seizures or seizure-free at the prior three visits, respectively. Only 25% of patients who were seizure-free for 6 months remained seizure-free after 10 years. SIGNIFICANCE: Our findings demonstrate that epilepsy outcome measures can be extracted accurately from unstructured clinical note text using NLP. At our tertiary center, the disease course often followed a remitting and relapsing pattern. This method represents a powerful new tool for clinical research with many potential uses and extensions to other clinical questions.

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.

A pharmacokinetic model of antiseizure medication load to guide care in the epilepsy monitoring unit.

OBJECTIVE: Evaluating patients with drug-resistant epilepsy often requires inducing seizures by tapering antiseizure medications (ASMs) in the epilepsy monitoring unit (EMU). The relationship between ASM taper strategy, seizure timing, and severity remains unclear. In this study, we developed and validated a pharmacokinetic model of total ASM load and tested its association with seizure occurrence and severity in the EMU. METHODS: We studied 80 patients who underwent intracranial electroencephalographic recording for epilepsy surgery planning. We developed a first order pharmacokinetic model of the ASMs administered in the EMU to generate a continuous metric of overall ASM load. We then related modeled ASM load to seizure likelihood and severity. We determined the association between the rate of ASM load reduction, the length of hospital stay, and the probability of having a severe seizure. Finally, we used modeled ASM load to predict oncoming seizures. RESULTS: Seizures occurred in the bottom 50th percentile of sampled ASM loads across the cohort (p < .0001, Wilcoxon signed-rank test), and seizures requiring rescue therapy occurred at lower ASM loads than seizures that did not require rescue therapy (logistic regression mixed effects model, odds ratio = .27, p = .01). Greater ASM decrease early in the EMU was not associated with an increased likelihood of having a severe seizure, nor with a shorter length of stay. SIGNIFICANCE: A pharmacokinetic model can accurately estimate ASM levels for patients in the EMU. Lower modeled ASM levels are associated with increased seizure likelihood and seizure severity. We show that ASM load, rather than ASM taper speed, is associated with severe seizures. ASM modeling has the potential to help optimize taper strategy to minimize severe seizures while maximizing diagnostic yield.

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.

Quantitative artifact reduction and pharmacologic paralysis improve detection of EEG epileptiform activity in critically ill patients.

OBJECTIVE: Epileptiform activity is common in critically ill patients, but movement-related artifacts-including electromyography (EMG) and myoclonus-can obscure EEG, limiting detection of epileptiform activity. We sought to determine the ability of pharmacologic paralysis and quantitative artifact reduction (AR) to improve epileptiform discharge detection. METHODS: Retrospective analysis of patients who underwent continuous EEG monitoring with pharmacologic paralysis. Four reviewers read each patient's EEG pre- and post- both paralysis and AR, and indicated the presence of epileptiform discharges. We compared the interrater reliability (IRR) of identifying discharges at baseline, post-AR, and post-paralysis, and compared the performance of AR and paralysis according to artifact type. RESULTS: IRR of identifying epileptiform discharges at baseline was slight (N = 30; κ = 0.10) with a trend toward increase post-AR (κ = 0.26, p = 0.053) and a significant increase post-paralysis (κ = 0.51, p = 0.001). AR was as effective as paralysis at improving IRR of identifying discharges in those with high EMG artifact (N = 15; post-AR κ = 0.63, p = 0.009; post-paralysis κ = 0.62, p = 0.006) but not with primarily myoclonus artifact (N = 15). CONCLUSIONS: Paralysis improves detection of epileptiform activity in critically ill patients when movement-related artifact obscures EEG features. AR improves detection as much as paralysis when EMG artifact is high, but is ineffective when the primary source of artifact is myoclonus. SIGNIFICANCE: In the appropriate setting, both AR and paralysis facilitate identification of epileptiform activity in critically ill patients.

Spike patterns surrounding sleep and seizures localize the seizure-onset zone in focal epilepsy.

OBJECTIVE: Interictal spikes help localize seizure generators as part of surgical planning for drug-resistant epilepsy. However, there are often multiple spike populations whose frequencies change over time, influenced by brain state. Understanding state changes in spike rates will improve our ability to use spikes for surgical planning. Our goal was to determine the effect of sleep and seizures on interictal spikes, and to use sleep and seizure-related changes in spikes to localize the seizure-onset zone (SOZ). METHODS: We performed a retrospective analysis of intracranial electroencephalography (EEG) data from patients with focal epilepsy. We automatically detected interictal spikes and we classified different time periods as awake or asleep based on the ratio of alpha to delta power, with a secondary analysis using the recently published SleepSEEG algorithm. We analyzed spike rates surrounding sleep and seizures. We developed a model to localize the SOZ using state-dependent spike rates. RESULTS: We analyzed data from 101 patients (54 women, age range 16-69). The normalized alpha-delta power ratio accurately classified wake from sleep periods (area under the curve = .90). Spikes were more frequent in sleep than wakefulness and in the post-ictal compared to the pre-ictal state. Patients with temporal lobe epilepsy had a greater wake-to-sleep and pre- to post-ictal spike rate increase compared to patients with extra-temporal epilepsy. A machine-learning classifier incorporating state-dependent spike rates accurately identified the SOZ (area under the curve = .83). Spike rates tended to be higher and better localize the seizure-onset zone in non-rapid eye movement (NREM) sleep than in wake or REM sleep. SIGNIFICANCE: The change in spike rates surrounding sleep and seizures differs between temporal and extra-temporal lobe epilepsy. Spikes are more frequent and better localize the SOZ in sleep, particularly in NREM sleep. Quantitative analysis of spikes may provide useful ancillary data to localize the SOZ and improve surgical planning.

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.

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.

Extracting seizure frequency from epilepsy clinic notes: a machine reading approach to natural language processing.

OBJECTIVE: Seizure frequency and seizure freedom are among the most important outcome measures for patients with epilepsy. In this study, we aimed to automatically extract this clinical information from unstructured text in clinical notes. If successful, this could improve clinical decision-making in epilepsy patients and allow for rapid, large-scale retrospective research. MATERIALS AND METHODS: We developed a finetuning pipeline for pretrained neural models to classify patients as being seizure-free and to extract text containing their seizure frequency and date of last seizure from clinical notes. We annotated 1000 notes for use as training and testing data and determined how well 3 pretrained neural models, BERT, RoBERTa, and Bio_ClinicalBERT, could identify and extract the desired information after finetuning. RESULTS: The finetuned models (BERTFT, Bio_ClinicalBERTFT, and RoBERTaFT) achieved near-human performance when classifying patients as seizure free, with BERTFT and Bio_ClinicalBERTFT achieving accuracy scores over 80%. All 3 models also achieved human performance when extracting seizure frequency and date of last seizure, with overall F1 scores over 0.80. The best combination of models was Bio_ClinicalBERTFT for classification, and RoBERTaFT for text extraction. Most of the gains in performance due to finetuning required roughly 70 annotated notes. DISCUSSION AND CONCLUSION: Our novel machine reading approach to extracting important clinical outcomes performed at or near human performance on several tasks. This approach opens new possibilities to support clinical practice and conduct large-scale retrospective clinical research. Future studies can use our finetuning pipeline with minimal training annotations to answer new clinical questions.

Implanting intracranial electrodes does not affect spikes or network connectivity in nearby or connected brain regions.

To determine the effect of implanting electrodes on electrographic features of nearby and connected brain regions in patients with drug-resistant epilepsy, we analyzed intracranial EEG recordings from 10 patients with drug-resistant epilepsy who underwent implant revision (placement of additional electrodes) during their hospitalization. We performed automated spike detection and measured EEG functional networks. We analyzed the original electrodes that remained in place throughout the full EEG recording, and we measured the change in spike rates and network connectivity in these original electrodes in response to implanting new electrodes. There was no change in overall spike rate pre- to post-implant revision (t(9) = 0.1, p = 0.95). The peri-revision change in the distribution of spike rate and connectivity across electrodes was no greater than chance (Monte Carlo method, spikes: p = 0.40, connectivity: p = 0.42). Electrodes closer to or more functionally connected to the revision site had no greater change in spike rate or connectivity than more distant or less connected electrodes. Changes in electrographic features surrounding electrode implantation are no greater than baseline fluctuations occurring throughout the intracranial recording. These findings argue against an implant effect on spikes or network connectivity in nearby or connected brain regions.

Using Generalized Polyspike Train to Predict Drug-Resistant Idiopathic Generalized Epilepsy.

INTRODUCTION: The authors tested the hypothesis that the EEG feature generalized polyspike train (GPT) is associated with drug-resistant idiopathic generalized epilepsy (IGE). METHODS: The authors conducted a single-center case-control study of patients with IGE who had outpatient EEGs performed between 2016 and 2020. The authors classified patients as drug-resistant or drug-responsive based on clinical review and in a masked manner reviewed EEG data for the presence and timing of GPT (a burst of generalized rhythmic spikes lasting less than 1 second) and other EEG features. A relationship between GPT and drug resistance was tested before and after controlling for EEG duration. The EEG duration needed to observe GPT was also calculated. RESULTS: One hundred three patients were included (70% drug-responsive and 30% drug-resistant patients). Generalized polyspike train was more prevalent in drug-resistant IGE (odds ratio, 3.8; 95% confidence interval, 1.3-11.4; P = 0.02). This finding persisted when controlling for EEG duration both with stratification and with survival analysis. A median of 6.5 hours (interquartile range, 0.5-12.7 hours) of EEG recording was required to capture the first occurrence of GPT. CONCLUSIONS: The findings support the hypothesis that GPT is associated with drug-resistant IGE. Prolonged EEG recording is required to identify this feature. Thus, >24-hour EEG recording early in the evaluation of patients with IGE may facilitate prognostication.

A Full-Stack Application for Detecting Seizures and Reducing Data During Continuous Electroencephalogram Monitoring.

Continuous electroencephalogram monitoring is associated with lower mortality in critically ill patients; however, it is underused due to the resource-intensive nature of manually interpreting prolonged streams of continuous electroencephalogram data. Here, we present a novel real-time, machine learning-based alerting and monitoring system for epilepsy and seizures that dramatically reduces the amount of manual electroencephalogram review. METHODS: We developed a custom data reduction algorithm using a random forest and deployed it within an online cloud-based platform, which streams data and communicates interactively with caregivers via a web interface to display algorithm results. We developed real-time, machine learning-based alerting and monitoring system for epilepsy and seizures on continuous electroencephalogram recordings from 77 patients undergoing routine scalp ICU electroencephalogram monitoring and tested it on an additional 20 patients. RESULTS: We achieved a mean seizure sensitivity of 84% in cross-validation and 85% in testing, as well as a mean specificity of 83% in cross-validation and 86% in testing, corresponding to a high level of data reduction. This study validates a platform for machine learning-assisted continuous electroencephalogram analysis and represents a meaningful step toward improving utility and decreasing cost of continuous electroencephalogram monitoring. We also make our high-quality annotated dataset of 97 ICU continuous electroencephalogram recordings public for others to validate and improve upon our methods.