Detection of interictal epileptiform discharges using signal envelope distribution modelling: application to epileptic and non-epileptic intracranial recordings

Timing matters for accurate identification of the epileptogenic zone

OBJECTIVE

Interictal biomarkers of the epileptogenic zone (EZ) and their use in machine learning models open promising avenues for improvement of epilepsy surgery evaluation. Currently, most studies restrict their analysis to short segments of intracranial EEG (iEEG).

METHODS

We used 2381 hours of iEEG data from 25 patients to systematically select 5-minute segments across various interictal conditions. Then, we tested machine learning models for EZ localization using iEEG features calculated within these individual segments or across them and evaluated the performance by the area under the precision-recall curve (PRAUC).

RESULTS

On average, models achieved a score of 0.421 (the result of the chance classifier was 0.062). However, the PRAUC varied significantly across the segments (0.323-0.493). Overall, NREM sleep achieved the highest scores, with the best results of 0.493 in N2. When using data from all segments, the model performed significantly better than single segments, except NREM sleep segments.

CONCLUSIONS

The model based on a short segment of iEEG recording can achieve similar results as a model based on prolonged recordings. The analyzed segment should, however, be carefully and systematically selected, preferably from NREM sleep.

SIGNIFICANCE

Random selection of short iEEG segments may give rise to inaccurate localization of the EZ.

An electroencephalogram microdisplay to visualize neuronal activity on the brain surface

Functional mapping during brain surgery is applied to define brain areas that control critical functions and cannot be removed. Currently, these procedures rely on verbal interactions between the neurosurgeon and electrophysiologist, which can be time-consuming. In addition, the electrode grids that are used to measure brain activity and to identify the boundaries of pathological versus functional brain regions have low resolution and limited conformity to the brain surface. Here, we present the development of an intracranial electroencephalogram (iEEG)-microdisplay that consists of freestanding arrays of 2048 GaN light-emitting diodes laminated on the back of micro-electrocorticography electrode grids. With a series of proof-of-concept experiments in rats and pigs, we demonstrate that these iEEG-microdisplays allowed us to perform real-time iEEG recordings and display cortical activities by spatially corresponding light patterns on the surface of the brain in the surgical field. Furthermore, iEEG-microdisplays allowed us to identify and display cortical landmarks and pathological activities from rat and pig models. Using a dual-color iEEG-microdisplay, we demonstrated coregistration of the functional cortical boundaries with one color and displayed the evolution of electrical potentials associated with epileptiform activity with another color. The iEEG-microdisplay holds promise to facilitate monitoring of pathological brain activity in clinical settings.

Natural language processing models reveal neural dynamics of human conversation

Through conversation, humans relay complex information through the alternation of speech production and comprehension. The neural mechanisms that underlie these complementary processes or through which information is precisely conveyed by language, however, remain poorly understood. Here, we used pretrained deep learning natural language processing models in combination with intracranial neuronal recordings to discover neural signals that reliably reflect speech production, comprehension, and their transitions during natural conversation between individuals. Our findings indicate that neural activities that encoded linguistic information were broadly distributed throughout frontotemporal areas across multiple frequency bands. We also find that these activities were specific to the words and sentences being conveyed and that they were dependent on the word's specific context and order. Finally, we demonstrate that these neural patterns partially overlapped during language production and comprehension and that listener-speaker transitions were associated with specific, time-aligned changes in neural activity. Collectively, our findings reveal a dynamical organization of neural activities that subserve language production and comprehension during natural conversation and harness the use of deep learning models in understanding the neural mechanisms underlying human language.

A platform for brain network sensing and stimulation with quantitative behavioral tracking: Application to limbic circuit epilepsy

Temporal lobe epilepsy is a common neurological disease characterized by recurrent seizures. These seizures often originate from limbic networks and people also experience chronic comorbidities related to memory, mood, and sleep (MMS). Deep brain stimulation targeting the anterior nucleus of the thalamus (ANT-DBS) is a proven therapy, but the optimal stimulation parameters remain unclear. We developed a neurotechnology platform for tracking seizures and MMS to enable data streaming between an investigational brain sensing-stimulation implant, mobile devices, and a cloud environment. Artificial Intelligence algorithms provided accurate catalogs of seizures, interictal epileptiform spikes, and wake-sleep brain states. Remotely administered memory and mood assessments were used to densely sample cognitive and behavioral response during ANT-DBS. We evaluated the efficacy of low-frequency versus high-frequency ANT-DBS. They both reduced seizures, but low-frequency ANT-DBS showed greater reductions and better sleep and memory. These results highlight the potential of synchronized brain sensing and behavioral tracking for optimizing neuromodulation therapy.

Method for cycle detection in sparse, irregularly sampled, long-term neuro-behavioral timeseries: Basis pursuit denoising with polynomial detrending of long-term, inter-ictal epileptiform activity

Numerous physiological processes are cyclical, but sampling these processes densely enough to perform frequency decomposition and subsequent analyses can be challenging. Mathematical approaches for decomposition and reconstruction of sparsely and irregularly sampled signals are well established but have been under-utilized in physiological applications. We developed a basis pursuit denoising with polynomial detrending (BPWP) model that recovers oscillations and trends from sparse and irregularly sampled timeseries. We validated this model on a unique dataset of long-term inter-ictal epileptiform discharge (IED) rates from human hippocampus recorded with a novel investigational device with continuous local field potential sensing. IED rates have well established circadian and multiday cycles related to sleep, wakefulness, and seizure clusters. Given sparse and irregular samples of IED rates from multi-month intracranial EEG recordings from ambulatory humans, we used BPWP to compute narrowband spectral power and polynomial trend coefficients and identify IED rate cycles in three subjects. In select cases, we propose that random and irregular sampling may be leveraged for frequency decomposition of physiological signals. Trial Registration: NCT03946618.

Benchmarking signal quality and spatiotemporal distribution of interictal spikes in prolonged human iEEG recordings using CorTec wireless brain interchange

Neuromodulation through implantable pulse generators (IPGs) represents an important treatment approach for neurological disorders. While the field has observed the success of state-of-the-art interventions, such as deep brain stimulation (DBS) or responsive neurostimulation (RNS), implantable systems face various technical challenges, including the restriction of recording from a limited number of brain sites, power management, and limited external access to the assessed neural data in a continuous fashion. To the best of our knowledge, for the first time in this study, we investigated the feasibility of recording human intracranial EEG (iEEG) using a benchtop version of the Brain Interchange (BIC) unit of CorTec, which is a portable, wireless, and externally powered implant with sensing and stimulation capabilities. We developed a MATLAB/SIMULINK-based rapid prototyping environment and a graphical user interface (GUI) to acquire and visualize the iEEG captured from all 32 channels of the BIC unit. We recorded prolonged iEEG (~ 24 h) from three human subjects with externalized depth leads using the BIC and commercially available clinical amplifiers simultaneously in the epilepsy monitoring unit (EMU). The iEEG signal quality of both streams was compared, and the results demonstrated a comparable power spectral density (PSD) in all the systems in the low-frequency band (< 80 Hz). However, notable differences were primarily observed above 100 Hz, where the clinical amplifiers were associated with lower noise floor (BIC-17 dB vs. clinical amplifiers <  - 25 dB). We employed an established spike detector to assess and compare the spike rates in each iEEG stream. We observed over 90% conformity between the spikes rates and their spatial distribution captured with BIC and clinical systems. Additionally, we quantified the packet loss characteristic in the iEEG signal during the wireless data transfer and conducted a series of simulations to compare the performance of different interpolation methods for recovering the missing packets in signals at different frequency bands. We noted that simple linear interpolation has the potential to recover the signal and reduce the noise floor with modest packet loss levels reaching up to 10%. Overall, our results indicate that while tethered clinical amplifiers exhibited noticeably better noise floor above 80 Hz, epileptic spikes can still be detected successfully in the iEEG recorded with the externally powered wireless BIC unit opening the road for future closed-loop neuromodulation applications with continuous access to brain activity.

Scalp and hippocampal sleep correlates of memory function in drug-resistant temporal lobe epilepsy

Seminal animal studies demonstrated the role of sleep oscillations such as cortical slow waves, thalamocortical spindles, and hippocampal ripples in memory consolidation. In humans, whether ripples are involved in sleep-related memory processes is less clear. Here, we explored the interactions between sleep oscillations (measured as traits) and general episodic memory abilities in 26 adults with drug-resistant temporal lobe epilepsy who performed scalp-intracranial electroencephalographic recordings and neuropsychological testing, including two analogous hippocampal-dependent verbal and nonverbal memory tasks. We explored the relationships between hemispheric scalp (spindles, slow waves) and hippocampal physiological and pathological oscillations (spindles, slow waves, ripples, and epileptic spikes) and material-specific memory function. To differentiate physiological from pathological ripples, we used multiple unbiased data-driven clustering approaches. At the individual level, we found material-specific cerebral lateralization effects (left-verbal memory, right-nonverbal memory) for all scalp spindles (rs > 0.51, ps < 0.01) and fast spindles (rs > 0.61, ps < 0.002). Hippocampal epileptic spikes and short pathological ripples, but not physiological oscillations, were negatively (rs > -0.59, ps < 0.01) associated with verbal learning and retention scores, with left lateralizing and antero-posterior effects. However, data-driven clustering failed to separate the ripple events into defined clusters. Correlation analyses with the resulting clusters revealed no meaningful or significant associations with the memory scores. Our results corroborate the role of scalp spindles in memory processes in patients with drug-resistant temporal lobe epilepsy. Yet, physiological and pathological ripples were not separable when using data-driven clustering, and thus our findings do not provide support for a role of sleep ripples as trait-like characteristics of general memory abilities in epilepsy.

Theta Oscillations Support Prefrontal-hippocampal Interactions in Sequential Working Memory

The prefrontal cortex and hippocampus may support sequential working memory beyond episodic memory and spatial navigation. This stereoelectroencephalography (SEEG) study investigated how the dorsolateral prefrontal cortex (DLPFC) interacts with the hippocampus in the online processing of sequential information. Twenty patients with epilepsy (eight women, age 27.6 ± 8.2 years) completed a line ordering task with SEEG recordings over the DLPFC and the hippocampus. Participants showed longer thinking times and more recall errors when asked to arrange random lines clockwise (random trials) than to maintain ordered lines (ordered trials) before recalling the orientation of a particular line. First, the ordering-related increase in thinking time and recall error was associated with a transient theta power increase in the hippocampus and a sustained theta power increase in the DLPFC (3-10 Hz). In particular, the hippocampal theta power increase correlated with the memory precision of line orientation. Second, theta phase coherences between the DLPFC and hippocampus were enhanced for ordering, especially for more precisely memorized lines. Third, the theta band DLPFC → hippocampus influence was selectively enhanced for ordering, especially for more precisely memorized lines. This study suggests that theta oscillations may support DLPFC-hippocampal interactions in the online processing of sequential information.

Mouse model of focal cortical dysplasia type II generates a wide spectrum of high-frequency activities

High-frequency oscillations (HFOs) represent an electrographic biomarker of endogenous epileptogenicity and seizure-generating tissue that proved clinically useful in presurgical planning and delineating the resection area. In the neocortex, the clinical observations on HFOs are not sufficiently supported by experimental studies stemming from a lack of realistic neocortical epilepsy models that could provide an explanation of the pathophysiological substrates of neocortical HFOs. In this study, we explored pathological epileptiform network phenomena, particularly HFOs, in a highly realistic murine model of neocortical epilepsy due to focal cortical dysplasia (FCD) type II. FCD was induced in mice by the expression of the human pathogenic mTOR gene mutation during embryonic stages of brain development. Electrographic recordings from multiple cortical regions in freely moving animals with FCD and epilepsy demonstrated that the FCD lesion generates HFOs from all frequency ranges, i.e., gamma, ripples, and fast ripples up to 800 Hz. Gamma-ripples were recorded almost exclusively in FCD animals, while fast ripples occurred in controls as well, although at a lower rate. Gamma-ripple activity is particularly valuable for localizing the FCD lesion, surpassing the utility of fast ripples that were also observed in control animals, although at significantly lower rates. Propagating HFOs occurred outside the FCD, and the contralateral cortex also generated HFOs independently of the FCD, pointing to a wider FCD network dysfunction. Optogenetic activation of neurons carrying mTOR mutation and expressing Channelrhodopsin-2 evoked fast ripple oscillations that displayed spectral and morphological profiles analogous to spontaneous oscillations. This study brings experimental evidence that FCD type II generates pathological HFOs across all frequency bands and provides information about the spatiotemporal properties of each HFO subtype in FCD. The study shows that mutated neurons represent a functionally interconnected and active component of the FCD network, as they can induce interictal epileptiform phenomena and HFOs.

The brain dynamics of visuospatial perspective-taking captured by intracranial EEG

Visuospatial perspective-taking (VPT) is the ability to imagine a scene from a position different from the one used in self-perspective judgments (SPJ). We typically use VPT to understand how others see the environment. VPT requires overcoming the self-perspective, and impairments in this process are implicated in various brain disorders, such as schizophrenia and autism. However, the underlying brain areas of VPT are not well distinguished from SPJ-related ones and from domain-general responses to both perspectives. In addition, hierarchical processing theory suggests that domain-specific processes emerge over time from domain-general ones. It mainly focuses on the sensory system, but outside of it, support for this hypothesis is lacking. Therefore, we aimed to spatiotemporally distinguish brain responses domain-specific to VPT from the specific ones to self-perspective, and domain-general responses to both perspectives. In particular, we intended to test whether VPT- and SPJ specific responses begin later than the general ones. We recorded intracranial EEG data from 30 patients with epilepsy who performed a task requiring laterality judgments during VPT and SPJ, and analyzed the spatiotemporal features of responses in the broad gamma band (50-150 Hz). We found VPT-specific processing in a more extensive brain network than SPJ-specific processing. Their dynamics were similar, but both differed from the general responses, which began earlier and lasted longer. Our results anatomically distinguish VPT-specific from SPJ-specific processing. Furthermore, we temporally differentiate between domain-specific and domain-general processes both inside and outside the sensory system, which serves as a novel example of hierarchical processing.

A review of signal processing and machine learning techniques for interictal epileptiform discharge detection

Brain interictal epileptiform discharges (IEDs), as one of the hallmarks of epileptic brain, are transient events captured by electroencephalogram (EEG). IEDs are generated by seizure networks, and they occur between seizures (interictal periods). The development of a robust method for IED detection could be highly informative for clinical treatment procedures and epileptic patient management. Since 1972, different machine learning techniques, from template matching to deep learning, have been developed to automatically detect IEDs from scalp EEG (scEEG) and intracranial EEG (iEEG). While the scEEG signals suffer from low information details and high attenuation of IEDs due to the high skull electrical impedance, the iEEG signals recorded using implanted electrodes enjoy higher details and are more suitable for identifying the IEDs. In this review paper, we group IED detection techniques into six categories: (1) template matching, (2) feature representation (mimetic, time-frequency, and nonlinear features), (3) matrix decomposition, (4) tensor factorization, (5) neural networks, and (6) estimation of the iEEG from the concurrent scEEG followed by detection and classification. The methods are compared quantitatively (e.g., in terms of accuracy, sensitivity, and specificity), and their general advantages and limitations are described. Finally, current limitations and possible future research paths related to this field are mentioned.

Timing of Allocentric and Egocentric Spatial Processing in Human Intracranial EEG

Spatial reference frames (RFs) play a key role in spatial cognition, especially in perception, spatial memory, and navigation. There are two main types of RFs: egocentric (self-centered) and allocentric (object-centered). Although many fMRI studies examined the neural correlates of egocentric and allocentric RFs, they could not sample the fast temporal dynamics of the underlying cognitive processes. Therefore, the interaction and timing between these two RFs remain unclear. Taking advantage of the high temporal resolution of intracranial EEG (iEEG), we aimed to determine the timing of egocentric and allocentric information processing and describe the brain areas involved. We recorded iEEG and analyzed broad gamma activity (50-150 Hz) in 37 epilepsy patients performing a spatial judgment task in a three-dimensional circular virtual arena. We found overlapping activation for egocentric and allocentric RFs in many brain regions, with several additional egocentric- and allocentric-selective areas. In contrast to the egocentric responses, the allocentric responses peaked later than the control ones in frontal regions with overlapping selectivity. Also, across several egocentric or allocentric selective areas, the egocentric selectivity appeared earlier than the allocentric one. We identified the maximum number of egocentric-selective channels in the medial occipito-temporal region and allocentric-selective channels around the intraparietal sulcus in the parietal cortex. Our findings favor the hypothesis that egocentric spatial coding is a more primary process, and allocentric representations may be derived from egocentric ones. They also broaden the dominant view of the dorsal and ventral streams supporting egocentric and allocentric space coding, respectively.

Impedance Rhythms in Human Limbic System

The impedance is a fundamental electrical property of brain tissue, playing a crucial role in shaping the characteristics of local field potentials, the extent of ephaptic coupling, and the volume of tissue activated by externally applied electrical brain stimulation. We tracked brain impedance, sleep-wake behavioral state, and epileptiform activity in five people with epilepsy living in their natural environment using an investigational device. The study identified impedance oscillations that span hours to weeks in the amygdala, hippocampus, and anterior nucleus thalamus. The impedance in these limbic brain regions exhibit multiscale cycles with ultradian (∼1.5-1.7 h), circadian (∼21.6-26.4 h), and infradian (∼20-33 d) periods. The ultradian and circadian period cycles are driven by sleep-wake state transitions between wakefulness, nonrapid eye movement (NREM) sleep, and rapid eye movement (REM) sleep. Limbic brain tissue impedance reaches a minimum value in NREM sleep, intermediate values in REM sleep, and rises through the day during wakefulness, reaching a maximum in the early evening before sleep onset. Infradian (∼20-33 d) impedance cycles were not associated with a distinct behavioral correlate. Brain tissue impedance is known to strongly depend on the extracellular space (ECS) volume, and the findings reported here are consistent with sleep-wake-dependent ECS volume changes recently observed in the rodent cortex related to the brain glymphatic system. We hypothesize that human limbic brain ECS changes during sleep-wake state transitions underlie the observed multiscale impedance cycles. Impedance is a simple electrophysiological biomarker that could prove useful for tracking ECS dynamics in human health, disease, and therapy.SIGNIFICANCE STATEMENT The electrical impedance in limbic brain structures (amygdala, hippocampus, anterior nucleus thalamus) is shown to exhibit oscillations over multiple timescales. We observe that impedance oscillations with ultradian and circadian periodicities are associated with transitions between wakefulness, NREM, and REM sleep states. There are also impedance oscillations spanning multiple weeks that do not have a clear behavioral correlate and whose origin remains unclear. These multiscale impedance oscillations will have an impact on extracellular ionic currents that give rise to local field potentials, ephaptic coupling, and the tissue activated by electrical brain stimulation. The approach for measuring tissue impedance using perturbational electrical currents is an established engineering technique that may be useful for tracking ECS volume.

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 Brain Electroencephalogram Microdisplay for Precision Neurosurgery

Brain surgeries are among the most delicate clinical procedures and must be performed with the most technologically robust and advanced tools. When such surgical procedures are performed in functionally critical regions of the brain, functional mapping is applied as a standard practice that involves direct coordinated interactions between the neurosurgeon and the clinical neurology electrophysiology team. However, information flow during these interactions is commonly verbal as well as time consuming which in turn increases the duration and cost of the surgery, possibly compromising the patient outcomes. Additionally, the grids that measure brain activity and identify the boundaries of pathological versus functional brain regions suffer from low resolution (3-10 mm contact to contact spacing) with limited conformity to the brain surface. Here, we introduce a brain intracranial electroencephalogram microdisplay (Brain-iEEG-microdisplay) which conforms to the brain to measure the brain activity and display changes in near real-time (40 Hz refresh rate) on the surface of the brain in the surgical field. We used scalable engineered gallium nitride (GaN) substrates with 6" diameter to fabricate, encapsulate, and release free-standing arrays of up to 2048 GaN light emitting diodes (μLEDs) in polyimide substrates. We then laminated the μLED arrays on the back of micro-electrocorticography (μECoG) platinum nanorod grids (PtNRGrids) and developed hardware and software to perform near real-time intracranial EEG analysis and activation of light patterns that correspond to specific cortical activities. Using the Brain-iEEG-microdisplay, we precisely ideFSntified and displayed important cortical landmarks and pharmacologically induced pathological activities. In the rat model, we identified and displayed individual cortical columns corresponding to individual whiskers and the near real-time evolution of epileptic discharges. In the pig animal model, we demonstrated near real-time mapping and display of cortical functional boundaries using somatosensory evoked potentials (SSEP) and display of responses to direct electrical stimulation (DES) from the surface or within the brain tissue. Using a dual-color Brain-iEEG-microdisplay, we demonstrated co-registration of the functional cortical boundaries with one color and displayed the evolution of electrical potentials associated with epileptiform activity with another color. The Brain-iEEG-microdisplay holds the promise of increasing the efficiency of diagnosis and possibly surgical treatment, thereby reducing the cost and improving patient outcomes which would mark a major advancement in neurosurgery. These advances can also be translated to broader applications in neuro-oncology and neurophysiology.

Distinct patterns of interictal intracranial EEG in focal cortical dysplasia type I and II

OBJECTIVE

Focal cortical dysplasia (FCD) is the most common malformation causing refractory focal epilepsy. Surgical removal of the entire dysplastic cortex is crucial for achieving a seizure-free outcome. Precise presurgical distinctions between FCD types by neuroimaging are difficult, mainly in patients with normal magnetic resonance imaging findings. However, the FCD type is important for planning the extent of surgical approach and counselling.

METHODS

This study included patients with focal drug-resistant epilepsy and definite histopathological FCD type I or II diagnoses who underwent intracranial electroencephalography (iEEG). We detected interictal epileptiform discharges (IEDs) and their recruitment into repetitive discharges (RDs) to compare electrophysiological patterns characterizing FCD types.

RESULTS

Patients with FCD type II had a significantly higher IED rate (p < 0.005), a shorter inter-discharge interval within RD episodes (p < 0.003), sleep influence on decreased RD periodicity (p < 0.036), and longer RD episode duration (p < 0.003) than patients with type I. A Bayesian classifier stratified FCD types with 82% accuracy.

CONCLUSION

Temporal characteristics of IEDs and RDs reflect the histological findings of FCD subtypes and can differentiate FCD types I and II.

SIGNIFICANCE

Presurgical prediction of FCD type can help to plan a more tailored surgical approach in patients with normal magnetic resonance findings.

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.

Identification and classification of pathology and artifacts for human intracranial cognitive research

Intracranial electroencephalography (iEEG) presents a unique opportunity to extend human neuroscientific understanding. However, typically iEEG is collected from patients diagnosed with focal drug-resistant epilepsy (DRE) and contains transient bursts of pathological activity. This activity disrupts performances on cognitive tasks and can distort findings from human neurophysiology studies. In addition to manual marking by a trained expert, numerous IED detectors have been developed to identify these pathological events. Even so, the versatility and usefulness of these detectors is limited by training on small datasets, incomplete performance metrics, and lack of generalizability to iEEG. Here, we employed a large annotated public iEEG dataset from two institutions to train a random forest classifier (RFC) to distinguish data segments as either 'non-cerebral artifact' (n = 73,902), 'pathological activity' (n = 67,797), or 'physiological activity' (n = 151,290). We found our model performed with an accuracy of 0.941, specificity of 0.950, sensitivity of 0.908, precision of 0.911, and F1 score of 0.910, averaged across all three event types. We extended the generalizability of our model to continuous bipolar data collected in a task-state at a different institution with a lower sampling rate and found our model performed with an accuracy of 0.789, specificity of 0.806, and sensitivity of 0.742, averaged across all three event types. Additionally, we created a custom graphical user interface to implement our classifier and enhance usability.

Closed-loop enhancement and neural decoding of cognitive control in humans

Deficits in cognitive control-that is, in the ability to withhold a default pre-potent response in favour of a more adaptive choice-are common in depression, anxiety, addiction and other mental disorders. Here we report proof-of-concept evidence that, in participants undergoing intracranial epilepsy monitoring, closed-loop direct stimulation of the internal capsule or striatum, especially the dorsal sites, enhances the participants' cognitive control during a conflict task. We also show that closed-loop stimulation upon the detection of lapses in cognitive control produced larger behavioural changes than open-loop stimulation, and that task performance for single trials can be directly decoded from the activity of a small number of electrodes via neural features that are compatible with existing closed-loop brain implants. Closed-loop enhancement of cognitive control might remediate underlying cognitive deficits and aid the treatment of severe mental disorders.

Perceived seizure risk in epilepsy â€" Chronic electronic surveys with and without concurrent EEG

Objective

Previous studies suggested that patients with epilepsy might be able to fore-cast their own seizures. We sought to assess the relationships of premonitory symptoms and perceived seizure risk with future and recent self-reported and EEG-confirmed seizures in the subjects living with epilepsy in their natural home environments.

Methods

We collected long-term e-surveys from ambulatory patients with and without concurrent EEG recordings. Information obtained from the e-surveys included medication compliance, sleep quality, mood, stress, perceived seizure risk and seizure occurrences preceding the survey. EEG seizures were identified. Univariate and multivariate generalized linear mixed-effect regression models were used to estimate odds ratios (ORs) for the assessment of the relationships. Results were compared with device seizure forecasting literature using a mathematical formula converting OR to equivalent area under the curve (AUC).

Results

Sixty-nine subjects returned 12,590 e-survey entries, with four subjects acquiring concurrent EEG recordings. Univariate analysis revealed increased stress (OR = 2.52, 95% CI = [1.52, 4.14], p < 0.001) and decreased mood (0.32, [0.13, 0.82], 0.02) were associated with increased relative odds of future self-reported seizures. On multivariate analysis, previous self-reported seizures (4.24, [2.69, 6.68], < 0.001) were most strongly associated with future self-reported seizures, and high perceived seizure risk (3.30, [1.97, 5.52], < 0.001) remained significant when prior self-reported seizures were added to the model. No significant association was found between e-survey responses and subsequent EEG seizures.

Significance

It appears that patients may tend to self-forecast seizures that occur in sequential groupings. Our results suggest that low mood and increased stress may be the result of previous seizures rather than independent premonitory symptoms. Patients in the small cohort with concurrent EEG showed no ability to self-predict EEG seizures. The conversion from OR to AUC values facilitates direct comparison of performance between survey and device studies involving survey premonition and forecasting.

Key points

Long-term e-surveys data and concurrent EEG signals were collected across three study sites to assess the ability of the patients to self-forecast their seizures.Patients may tend to self-forecast self-reported seizures that occur in sequential groupings.Factors, such as mood and stress, may not be independent premonitory symptoms but may be the consequence of recent seizures.No ability to self-forecast EEG confirmed seizures was observed in a small cohort with concurrent EEG validation.A mathematic relation between OR and AUC provides a means to compare forecasting performance between survey and device studies.

A Subpopulation of Spikes Predicts Successful Epilepsy Surgery Outcome

OBJECTIVE

Epileptic spikes are the traditional interictal electroencephalographic (EEG) biomarker for epilepsy. Given their low specificity for identifying the epileptogenic zone (EZ), they are given only moderate attention in presurgical evaluation. This study aims to demonstrate that it is possible to identify specific spike features in intracranial EEG that optimally define the EZ and predict surgical outcome.

METHODS

We analyzed spike features on stereo-EEG segments from 83 operated patients from 2 epilepsy centers (37 Engel IA) in wakefulness, non-rapid eye movement sleep, and rapid eye movement sleep. After automated spike detection, we investigated 135 spike features based on rate, morphology, propagation, and energy to determine the best feature or feature combination to discriminate the EZ in seizure-free and non-seizure-free patients by applying 4-fold cross-validation.

RESULTS

The rate of spikes with preceding gamma activity in wakefulness performed better for surgical outcome classification (4-fold area under receiver operating characteristics curve [AUC] = 0.755 ± 0.07) than the seizure onset zone, the current gold standard (AUC = 0.563 ± 0.05, p = 0.015) and the ripple rate, an emerging seizure-independent biomarker (AUC = 0.537 ± 0.07, p = 0.006). Channels with a spike-gamma rate exceeding 1.9/min had an 80% probability of being in the EZ. Combining features did not improve the results.

INTERPRETATION

Resection of brain regions with high spike-gamma rates in wakefulness is associated with a high probability of achieving seizure freedom. This rate could be applied to determine the minimal number of spiking channels requiring resection. In addition to quantitative analysis, this feature is easily accessible to visual analysis, which could aid clinicians during presurgical evaluation. ANN NEUROL 2023;93:522-535.

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.

Daily resting-state intracranial EEG connectivity for seizure risk forecasts

Forecasting seizure risk aims to detect proictal states in which seizures would be more likely to occur. Classical seizure prediction models are trained over long-term electroencephalographic (EEG) recordings to detect specific preictal changes for each seizure, independently of those induced by shifts in states of vigilance. A daily single measure-during a vigilance-controlled period-to estimate the risk of upcoming seizure(s) would be more convenient. Here, we evaluated whether intracranial EEG connectivity (phase-locking value), estimated from daily vigilance-controlled resting-state recordings, could allow distinguishing interictal (no seizure) from preictal (seizure within the next 24 h) states. We also assessed its relevance for daily forecasts of seizure risk using machine learning models. Connectivity in the theta band was found to provide the best prediction performances (area under the curve ≥ .7 in 80% of patients), with accurate daily and prospective probabilistic forecasts (mean Brier score and Brier skill score of .13 and .72, respectively). More efficient ambulatory clinical application could be considered using mobile EEG or chronic implanted devices.

Utilization of temporal autoencoder for semi-supervised intracranial EEG clustering and classification

Manual visual review, annotation and categorization of electroencephalography (EEG) is a time-consuming task that is often associated with human bias and requires trained electrophysiology experts with specific domain knowledge. This challenge is now compounded by development of measurement technologies and devices allowing large-scale heterogeneous, multi-channel recordings spanning multiple brain regions over days, weeks. Currently, supervised deep-learning techniques were shown to be an effective tool for analyzing big data sets, including EEG. However, the most significant caveat in training the supervised deep-learning models in a clinical research setting is the lack of adequate gold-standard annotations created by electrophysiology experts. Here, we propose a semi-supervised machine learning technique that utilizes deep-learning methods with a minimal amount of gold-standard labels. The method utilizes a temporal autoencoder for dimensionality reduction and a small number of the expert-provided gold-standard labels used for kernel density estimating (KDE) maps. We used data from electrophysiological intracranial EEG (iEEG) recordings acquired in two hospitals with different recording systems across 39 patients to validate the method. The method achieved iEEG classification (Pathologic vs. Normal vs. Artifacts) results with an area under the receiver operating characteristic (AUROC) scores of 0.862 ± 0.037, 0.879 ± 0.042, and area under the precision-recall curve (AUPRC) scores of 0.740 ± 0.740, 0.714 ± 0.042. This demonstrates that semi-supervised methods can provide acceptable results while requiring only 100 gold-standard data samples in each classification category. Subsequently, we deployed the technique to 12 novel patients in a pseudo-prospective framework for detecting Interictal epileptiform discharges (IEDs). We show that the proposed temporal autoencoder was able to generalize to novel patients while achieving AUROC of 0.877 ± 0.067 and AUPRC of 0.705 ± 0.154.

A consensus statement on detection of hippocampal sharp wave ripples and differentiation from other fast oscillations

Decades of rodent research have established the role of hippocampal sharp wave ripples (SPW-Rs) in consolidating and guiding experience. More recently, intracranial recordings in humans have suggested their role in episodic and semantic memory. Yet, common standards for recording, detection, and reporting do not exist. Here, we outline the methodological challenges involved in detecting ripple events and offer practical recommendations to improve separation from other high-frequency oscillations. We argue that shared experimental, detection, and reporting standards will provide a solid foundation for future translational discovery.

Computerized application for epilepsy in China: Does the era of artificial intelligence comes?

Epilepsy, one of the most common neurological diseases in China, is notorious for its spontaneous, unprovoked and recurrent seizures. The etiology of epilepsy varies among individual patients, including congenital gene mutation, traumatic injury, infections, etc. This heterogeneity partly hampered the accurate diagnosis and choice of appropriate treatments. Encouragingly, great achievements have been achieved in computational science, making it become a key player in medical fields gradually and bringing new hope for rapid and accurate diagnosis as well as targeted therapies in epilepsy. Here, we historically review the advances of computerized applications in epilepsy-especially those tremendous findings achieved in China-for different purposes including seizure prediction, localization of epileptogenic zone, post-surgical prognosis, etc. Special attentions are paid to the great progress based on artificial intelligence (AI), which is more "sensitive", "smart" and "in-depth" than human capacities. At last, we give a comprehensive discussion about the disadvantages and limitations of current computerized applications for epilepsy and propose some future directions as further stepping stones to embrace "the era of AI" in epilepsy.

Satelight: Self-attention-based model for epileptic spike detection from multi-electrode EEG

OBJECTIVE

Because of the lack of highly skilled experts, automated technologies that support electroencephalogram (EEG)-based in epilepsy diagnosis are advancing. Deep convolutional neural network-based models have been used successfully for detecting epileptic spikes, one of the biomarkers, from EEG. However, a sizeable number of supervised EEG records are required for training.

APPROACH

This study introduces the Satelight model, which uses the self-attention (SA) mechanism. The model was trained using a clinical EEG dataset labeled by  ve specialists, including 16,008 epileptic spikes and 15,478 artifacts from 50 children. The SA mechanism is expected to reduce the number of parameters and efficiently extract features from a small amount of EEG data. To validate the effectiveness, we compared various spike detection approaches with the clinical EEG data.

MAIN RESULTS

The experimental results showed that the proposed method detected epileptic spikes more effectively than other models (accuracy = 0:876 and false positive rate = 0:133).

SIGNIFICANCE

The proposed model had only one-tenth the number of parameters as the other effective model, despite having such a high detection performance. Further exploration of the hidden parameters revealed that the model automatically attended to the EEG's characteristic waveform locations of interest.

Evoked responses to note onsets and phrase boundaries in Mozart's K448

Understanding the neural correlates of perception of hierarchical structure in music presents a direct window into auditory organization. To examine the hypothesis that high-level and low-level structures—i.e. phrases and notes—elicit different neural responses, we collected intracranial electroencephalography (iEEG) data from eight subjects during exposure to Mozart’s K448 and directly compared Event-related potentials (ERPs) due to note onsets and those elicited by phrase boundaries. Cluster-level permutation tests revealed that note-onset-related ERPs and phrase-boundary-related ERPs were significantly different at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-150$$\end{document}-150, 200, and 450 ms relative to note onset and phrase markers. We also observed increased activity in frontal brain regions when processing phrase boundaries. We relate these observations to (1) a process which syntactically binds notes together hierarchically to form larger phrases; (2) positive emotions induced by successful prediction of forthcoming phrase boundaries and violations of melodic expectations at phrase boundaries.

Louis EK, Gregg NM, Lundstrom BN, Crowe CM, Attia TP, Crepeau D, Balzekas I, Marks VS, Wheeler LP, Cimbalnik J, Cook M, Janca R, Sturges BK, Leyde K, Miller KJ, Van Gompel JJ, Denison T, Worrell GA, Kremen V. Distributed brain co-processor for tracking spikes, seizures and behavior during electrical brain stimulation

Early implantable epilepsy therapy devices provided open-loop electrical stimulation without brain sensing, computing, or an interface for synchronized behavioral inputs from patients. Recent epilepsy stimulation devices provide brain sensing but have not yet developed analytics for accurately tracking and quantifying behavior and seizures. Here we describe a distributed brain co-processor providing an intuitive bi-directional interface between patient, implanted neural stimulation and sensing device, and local and distributed computing resources. Automated analysis of continuous streaming electrophysiology is synchronized with patient reports using a hand-held device and integrated with distributed cloud computing resources for quantifying seizures, interictal epileptiform spikes, and patient symptoms during therapeutic electrical brain stimulation. The classification algorithms for interictal epileptiform spikes and seizures were developed and parameterized using long-term ambulatory data from 9 humans and 8 canines with epilepsy, and then implemented prospectively in out-of-sample testing in 2 pet canines and 4 humans with drug resistant epilepsy living in their natural environments. Accurate seizure diaries are needed as the primary clinical outcome measure of epilepsy therapy and to guide brain stimulation optimization. The brain co-processor system described here enables tracking interictal epileptiform spikes, seizures, and correlation with patient behavioral reports. In the future correlation of spikes and seizures with behavior will allow more detailed investigation of the clinical impact of spikes and seizures on patients.

SleepSEEG: Automatic sleep scoring using intracranial EEG recordings only

OBJECTIVE

To perform automatic sleep scoring based only on intracranial EEG, without the need for scalp electroencephalography (EEG), electrooculography (EOG) and electromyography (EMG), in order to study sleep, epilepsy, and their interaction.

APPROACH

Data from 33 adult patients was used for development and training of the automatic scoring algorithm using both oscillatory and non-oscillatory spectral features. The first step consisted in unsupervised clustering of channels based on feature variability. For each cluster the classification was done in two steps, a multiclass tree followed by binary classification trees to distinguish the more challenging stage N1. The test data consisted in 11 patients, in whom the classification was done independently for each channel and then combined to get a single stage per epoch.

MAIN RESULTS

An overall agreement of 78% was observed in the test set between the sleep scoring of the algorithm and two human experts scoring based on scalp EEG, EOG and EMG. Balanced sensitivity and specificity were obtained for the different sleep stages. The performance was excellent for stages W, N2, and N3, and good for stage R, but with high variability across patients. The performance for the challenging stage N1 was poor, but at a similar level as for published algorithms based on scalp EEG. High confidence epochs in different stages (other than N1) can be identified with median per patient specificity >80%.

SIGNIFICANCE

The automatic algorithm can perform sleep scoring of long term recordings of patients with intracranial electrodes undergoing presurgical evaluation in the absence of scalp EEG, EOG and EMG, which are normally required to define sleep stages but are difficult to use in the context of intracerebral studies. It also constitutes a valuable tool to generate hypotheses regarding local aspects of sleep, and will be significant for sleep evaluation in clinical epileptology and neuroscience research.

Neurophysiological Evidence for Cognitive Map Formation during Sequence Learning

Humans deftly parse statistics from sequences. Some theories posit that humans learn these statistics by forming cognitive maps, or underlying representations of the latent space which links items in the sequence. Here, an item in the sequence is a node, and the probability of transitioning between two items is an edge. Sequences can then be generated from walks through the latent space, with different spaces giving rise to different sequence statistics. Individual or group differences in sequence learning can be modeled by changing the time scale over which estimates of transition probabilities are built, or in other words, by changing the amount of temporal discounting. Latent space models with temporal discounting bear a resemblance to models of navigation through Euclidean spaces. However, few explicit links have been made between predictions from Euclidean spatial navigation and neural activity during human sequence learning. Here, we use a combination of behavioral modeling and intracranial encephalography (iEEG) recordings to investigate how neural activity might support the formation of space-like cognitive maps through temporal discounting during sequence learning. Specifically, we acquire human reaction times from a sequential reaction time task, to which we fit a model that formulates the amount of temporal discounting as a single free parameter. From the parameter, we calculate each individual's estimate of the latent space. We find that neural activity reflects these estimates mostly in the temporal lobe, including areas involved in spatial navigation. Similar to spatial navigation, we find that low-dimensional representations of neural activity allow for easy separation of important features, such as modules, in the latent space. Lastly, we take advantage of the high temporal resolution of iEEG data to determine the time scale on which latent spaces are learned. We find that learning typically happens within the first 500 trials, and is modulated by the underlying latent space and the amount of temporal discounting characteristic of each participant. Ultimately, this work provides important links between behavioral models of sequence learning and neural activity during the same behavior, and contextualizes these results within a broader framework of domain general cognitive maps.

Local and distant cortical responses to single pulse intracranial stimulation in the human brain are differentially modulated by specific stimulation parameters

BACKGROUND

Electrical neuromodulation via direct electrical stimulation (DES) is an increasingly common therapy for a wide variety of neuropsychiatric diseases. Unfortunately, therapeutic efficacy is inconsistent, likely due to our limited understanding of the relationship between the massive stimulation parameter space and brain tissue responses.

OBJECTIVE

To better understand how different parameters induce varied neural responses, we systematically examined single pulse-induced cortico-cortico evoked potentials (CCEP) as a function of stimulation amplitude, duration, brain region, and whether grey or white matter was stimulated.

METHODS

We measured voltage peak amplitudes and area under the curve (AUC) of intracranially recorded stimulation responses as a function of distance from the stimulation site, pulse width, current injected, location relative to grey and white matter, and brain region stimulated (N = 52, n = 719 stimulation sites).

RESULTS

Increasing stimulation pulse width increased responses near the stimulation location. Increasing stimulation amplitude (current) increased both evoked amplitudes and AUC nonlinearly. Locally (<15 mm), stimulation at the boundary between grey and white matter induced larger responses. In contrast, for distant sites (>15 mm), white matter stimulation consistently produced larger responses than stimulation in or near grey matter. The stimulation location-response curves followed different trends for cingulate, lateral frontal, and lateral temporal cortical stimulation.

CONCLUSION

These results demonstrate that a stronger local response may require stimulation in the grey-white boundary while stimulation in the white matter could be needed for network activation. Thus, stimulation parameters tailored for a specific anatomical-functional outcome may be key to advancing neuromodulatory therapy.

Large-scale neural recordings with single neuron resolution using Neuropixels probes in human cortex

Recent advances in multi-electrode array technology have made it possible to monitor large neuronal ensembles at cellular resolution in animal models. In humans, however, current approaches restrict recordings to a few neurons per penetrating electrode or combine the signals of thousands of neurons in local field potential (LFP) recordings. Here we describe a new probe variant and set of techniques that enable simultaneous recording from over 200 well-isolated cortical single units in human participants during intraoperative neurosurgical procedures using silicon Neuropixels probes. We characterized a diversity of extracellular waveforms with eight separable single-unit classes, with differing firing rates, locations along the length of the electrode array, waveform spatial spread and modulation by LFP events such as inter-ictal discharges and burst suppression. Although some challenges remain in creating a turnkey recording system, high-density silicon arrays provide a path for studying human-specific cognitive processes and their dysfunction at unprecedented spatiotemporal resolution.

Human brain mapping with multithousand-channel PtNRGrids resolves spatiotemporal dynamics

Electrophysiological devices are critical for mapping eloquent and diseased brain regions and for therapeutic neuromodulation in clinical settings and are extensively used for research in brain-machine interfaces. However, the existing clinical and experimental devices are often limited in either spatial resolution or cortical coverage. Here, we developed scalable manufacturing processes with a dense electrical connection scheme to achieve reconfigurable thin-film, multithousand-channel neurophysiological recording grids using platinum nanorods (PtNRGrids). With PtNRGrids, we have achieved a multithousand-channel array of small (30 μm) contacts with low impedance, providing high spatial and temporal resolution over a large cortical area. We demonstrated that PtNRGrids can resolve submillimeter functional organization of the barrel cortex in anesthetized rats that captured the tissue structure. In the clinical setting, PtNRGrids resolved fine, complex temporal dynamics from the cortical surface in an awake human patient performing grasping tasks. In addition, the PtNRGrids identified the spatial spread and dynamics of epileptic discharges in a patient undergoing epilepsy surgery at 1-mm spatial resolution, including activity induced by direct electrical stimulation. Collectively, these findings demonstrated the power of the PtNRGrids to transform clinical mapping and research with brain-machine interfaces.

Thalamic deep brain stimulation modulates cycles of seizure risk in epilepsy

Chronic brain recordings suggest that seizure risk is not uniform, but rather varies systematically relative to daily (circadian) and multiday (multidien) cycles. Here, one human and seven dogs with naturally occurring epilepsy had continuous intracranial EEG (median 298 days) using novel implantable sensing and stimulation devices. Two pet dogs and the human subject received concurrent thalamic deep brain stimulation (DBS) over multiple months. All subjects had circadian and multiday cycles in the rate of interictal epileptiform spikes (IES). There was seizure phase locking to circadian and multiday IES cycles in five and seven out of eight subjects, respectively. Thalamic DBS modified circadian (all 3 subjects) and multiday (analysis limited to the human participant) IES cycles. DBS modified seizure clustering and circadian phase locking in the human subject. Multiscale cycles in brain excitability and seizure risk are features of human and canine epilepsy and are modifiable by thalamic DBS.

Spatio-temporal spike dynamics predict surgical outcome in adult focal epilepsy

OBJECTIVE

We hypothesized that spatio-temporal dynamics of interictal spikes reflect the extent and stability of epileptic sources and determine surgical outcome.

METHODS

We studied 30 consecutive patients (14 good outcome). Spikes were detected in prolonged stereo-electroencephalography recordings. We quantified the spatio-temporal dynamics of spikes using the variance of the spike rate, line length and skewness of the spike distribution, and related these features to outcome. We built a logistic regression model, and compared its performance to traditional markers.

RESULTS

Good outcome patients had more dominant and stable sources than poor outcome patients as expressed by a higher variance of spike rates, a lower variance of line length, and a lower variance of positive skewness (ps < 0.05). The outcome was correctly predicted in 80% of patients. This was better or non-inferior to predictions based on a focal lesion (p = 0.016), focal seizure-onset zone, or complete resection (ps > 0.05). In the five patients where traditional markers failed, spike distribution predicted the outcome correctly. The best results were achieved by 18-h periods or longer.

CONCLUSIONS

Analysis of spike dynamics shows that surgery outcome depends on strong, single and stable sources.

SIGNIFICANCE

Our quantitative method has the potential to be a reliable predictor of surgical outcome.

Long Time Constant May Endorses Sharp Waves and Spikes Than Sharp Transients in Scalp Electroencephalography: A Comparison of Both After-Slow Among Different Time Constant and High-Frequency Activity Analysis

Objective: To clarify whether long time constant (TC) is useful for detecting the after-slow activity of epileptiform discharges (EDs): sharp waves and spikes and for differentiating EDs from sharp transients (Sts).

Methods: We employed 68 after-slow activities preceded by 32 EDs (26 sharp waves and six spikes) and 36 Sts from 52 patients with partial and generalized epilepsy (22 men, 30 women; mean age 39.08 ± 13.13 years) defined by visual inspection. High-frequency activity (HFA) associated with the apical component of EDs and Sts was also investigated to endorse two groups. After separating nine Sts that were labeled by visual inspection but did not fulfill the amplitude criteria for after-slow of Sts, 59 activities (32 EDs and 27 Sts) were analyzed about the total area of after-slow under three TCs (long: 2 s; conventional: 0.3 s; and short: 0.1 s).

Results: Compared to Sts, HFA was found significantly more with the apical component of EDs (p < 0.05). The total area of after-slow in all 32 EDs under TC 2 s was significantly larger than those under TC 0.3 s and 0.1 s (p < 0.001). Conversely, no significant differences were observed in the same parameter of 27 Sts among the three different TCs. Regarding separated nine Sts, the total area of after-slow showed a similar tendency to that of 27 Sts under three different TCs.

Significance: These results suggest that long TC could be useful for selectively endorsing after-slow of EDs and differentiating EDs from Sts. These findings are concordant with the results of the HFA analysis. Visual inspection is also equally good as the total area of after-slow analysis.

AiED: Artificial intelligence for the detection of intracranial interictal epileptiform discharges

OBJECTIVE

Deep learning provides an appealing solution for the ongoing challenge of automatically classifying intracranial interictal epileptiform discharges (IEDs). We report results from an automated method consisting of a template-matching algorithm and convolutional neural network (CNN) for the detection of intracranial IEDs ("AiED").

METHODS

1000 intracranial electroencephalogram (EEG) epochs extracted randomly from 307 subjects with refractory epilepsy were annotated independently by two expert neurophysiologists. These annotated epochs were divided into 1062 two-second epochs with IEDs and 1428 two-second epochs without IEDs, which were transformed into spectrograms prior to training the neural network. The highest performing network was validated on an annotated external test set.

RESULTS

The final network had an F1-score of 0.95 (95% CI: 0.91-0.98) and an average Area Under the Receiver Operating Characteristic of 0.98 (95% CI: 0.96-1.00). For the external test set, it showed an overall F1-score of 0.71, correctly identifying 100% of all high-amplitude IED complexes, 96.23% of all high-amplitude isolated IEDs, and 66.15% of all IEDs of atypical morphology.

CONCLUSIONS

Template-matching combined with a CNN offers a fast, robust method for detecting intracranial IEDs.

SIGNIFICANCE

"AiED" is generalizable and achieves comparable performance to human reviewers; it may support clinical and research EEG analyses.

Intracranial EEG Biomarkers for Seizure Lateralization in Rapidly-Bisynchronous Epilepsy After Laser Corpus Callosotomy

Objective: It has been asserted that high-frequency analysis of intracranial EEG (iEEG) data may yield information useful in localizing epileptogenic foci. Methods: We tested whether proposed biomarkers could predict lateralization based on iEEG data collected prior to corpus callosotomy (CC) in three patients with bisynchronous epilepsy, whose seizures lateralized definitively post-CC. Lateralization data derived from algorithmically-computed ictal phase-locked high gamma (PLHG), high gamma amplitude (HGA), and low-frequency (filtered) line length (LFLL), as well as interictal high-frequency oscillation (HFO) and interictal epileptiform discharge (IED) rate metrics were compared against ground-truth lateralization from post-CC ictal iEEG. Results: Pre-CC unilateral IEDs were more frequent on the more-pathologic side in all subjects. HFO rate predicted lateralization in one subject, but was sensitive to detection threshold. On pre-CC data, no ictal metric showed better predictive power than any other. All post-corpus callosotomy seizures lateralized to the pathological hemisphere using PLHG, HGA, and LFLL metrics. Conclusions: While quantitative metrics of IED rate and ictal HGA, PHLG, and LFLL all accurately lateralize based on post-CC iEEG, only IED rate consistently did so based on pre-CC data. Significance: Quantitative analysis of IEDs may be useful in lateralizing seizure pathology. More work is needed to develop reliable techniques for high-frequency iEEG analysis.

Effects of hippocampal interictal discharge timing, duration, and spatial extent on list learning

Interictal epileptiform discharges (IEDs) can impair memory. The properties of IEDs most detrimental to memory, however, are undefined. We studied the impact of temporal and spatial characteristics of IEDs on list learning. Subjects completed a memory task during intracranial EEG recordings including hippocampal depth and temporal neocortical subdural electrodes. Subjects viewed a series of objects, and after a distracting task, recalled the objects from the list. The impacts of IED presence, duration, and propagation to neocortex during encoding of individual stimuli were assessed. The effects of IED total number and duration during maintenance and recall periods on delayed recall performance were also determined. The influence of IEDs during recall was further investigated by comparing the likelihood of IEDs preceding correctly recalled items vs. periods of no verbal response. Across 6 subjects, we analyzed 28 hippocampal and 139 lateral temporal contacts. Recall performance was poor, with a median of 17.2% correct responses (range 10.4-21.9%). Interictal epileptiform discharges during encoding, maintenance, and recall did not significantly impact task performance, and there was no significant difference between the likelihood of IEDs during correct recall vs. periods of no response. No significant effects of discharge duration during encoding, maintenance, or recall were observed. Interictal epileptiform discharges with spread to lateral temporal cortex during encoding did not adversely impact recall. A post hoc analysis refining model assumptions indicated a negative impact of IED count during the maintenance period, but otherwise confirmed the above results. Our findings suggest no major effect of hippocampal IEDs on list learning, but study limitations, such as baseline hippocampal dysfunction, should be considered. The impact of IEDs during the maintenance period may be a focus of future research.

Musical components important for the Mozart K448 effect in epilepsy

There is growing evidence for the efficacy of music, specifically Mozart’s Sonata for Two Pianos in D Major (K448), at reducing ictal and interictal epileptiform activity. Nonetheless, little is known about the mechanism underlying this beneficial “Mozart K448 effect” for persons with epilepsy. Here, we measured the influence that K448 had on intracranial interictal epileptiform discharges (IEDs) in sixteen subjects undergoing intracranial monitoring for refractory focal epilepsy. We found reduced IEDs during the original version of K448 after at least 30-s of exposure. Nonsignificant IED rate reductions were witnessed in all brain regions apart from the bilateral frontal cortices, where we observed increased frontal theta power during transitions from prolonged musical segments. All other presented musical stimuli were associated with nonsignificant IED alterations. These results suggest that the “Mozart K448 effect” is dependent on the duration of exposure and may preferentially modulate activity in frontal emotional networks, providing insight into the mechanism underlying this response. Our findings encourage the continued evaluation of Mozart’s K448 as a noninvasive, non-pharmacological intervention for refractory epilepsy.

Microscale dynamics of electrophysiological markers of epilepsy

OBJECTIVE

Interictal discharges (IIDs) and high frequency oscillations (HFOs) are established neurophysiologic biomarkers of epilepsy, while microseizures are less well studied. We used custom poly(3,4-ethylenedioxythiophene) polystyrene sulfonate (PEDOT:PSS) microelectrodes to better understand these markers' microscale spatial dynamics.

METHODS

Electrodes with spatial resolution down to 50 µm were used to record intraoperatively in 30 subjects. IIDs' degree of spread and spatiotemporal paths were generated by peak-tracking followed by clustering. Repeating HFO patterns were delineated by clustering similar time windows. Multi-unit activity (MUA) was analyzed in relation to IID and HFO timing.

RESULTS

We detected IIDs encompassing the entire array in 93% of subjects, while localized IIDs, observed across < 50% of channels, were seen in 53%. IIDs traveled along specific paths. HFOs appeared in small, repeated spatiotemporal patterns. Finally, we identified microseizure events that spanned 50-100 µm. HFOs covaried with MUA, but not with IIDs.

CONCLUSIONS

Overall, these data suggest that irritable cortex micro-domains may form part of an underlying pathologic architecture which could contribute to the seizure network.

SIGNIFICANCE

These results, supporting the possibility that epileptogenic cortex comprises a mosaic of irritable domains, suggests that microscale approaches might be an important perspective in devising novel seizure control therapies.

Theta oscillations synchronize human medial prefrontal cortex and amygdala during fear learning

Numerous animal studies have demonstrated that fear acquisition and expression rely on the coordinated activity of medial prefrontal cortex (mPFC) and amygdala and that theta oscillations support interregional communication within the fear network. However, it remains unclear whether these results can be generalized to fear learning in humans. We addressed this question using intracranial electroencephalography recordings in 13 patients with epilepsy during a fear conditioning paradigm. We observed increased power and inter-regional synchronization of amygdala and mPFC in theta (4 to 8 hertz) oscillations for conditioned stimulus (CS+) versus CS-. Analysis of information flow revealed that the dorsal mPFC (dmPFC) led amygdala activity in theta oscillations. Last, a computational model showed that trial-by-trial changes in amygdala theta oscillations predicted the model-based associability (i.e., learning rate). This study provides compelling evidence that theta oscillations within and between amygdala, ventral mPFC, and dmPFC constitute a general mechanism of fear learning across species.

40-Hz auditory stimulation for intracranial interictal activity: A pilot study

OBJECTIVES

To study the effects of auditory stimuli on interictal epileptiform discharge (IED) rates evident with intracranial monitoring.

MATERIALS AND METHODS

Eight subjects undergoing intracranial EEG monitoring for refractory epilepsy participated in this study. Auditory stimuli consisted of a 40-Hz tone, a 440-Hz tone modulated by a 40-Hz sinusoid, Mozart's Sonata for Two Pianos in D Major (K448), and K448 modulated by a 40-Hz sinusoid (modK448). Subjects were stratified into high- and low-IED rate groups defined by baseline IED rates. Subject-level analyses identified individual responses to auditory stimuli, discerned specific brain regions with significant reductions in IED rates, and examined the influence auditory stimuli had on whole-brain sigma power (12-16 Hz).

RESULTS

All subjects in the high baseline IED group had a significant 35.25% average reduction in IEDs during the 40-Hz tone; subject-level reductions localized to mesial and lateral temporal regions. Exposure to Mozart K448 showed significant yet less homogeneous responses. A post hoc analysis demonstrated two of the four subjects with positive IED responses had increased whole-brain power at the sigma frequency band during 40-Hz stimulation.

CONCLUSIONS

Our study is the first to evaluate the relationship between 40-Hz auditory stimulation and IED rates in refractory epilepsy. We reveal that 40-Hz auditory stimuli may be a noninvasive adjunctive intervention to reduce IED burden. Our pilot study supports the future examination of 40-Hz auditory stimuli in a larger population of subjects with high baseline IED rates.

Anterior thalamic stimulation improves working memory precision judgments

BACKGROUND

The anterior nucleus of thalamus (ANT) has been suggested as an extended hippocampal system. The circuit of ANT and hippocampus has been widely demonstrated to be associated with memory function. Both lesions to each region and disrupting inter-regional information flow can induce working memory impairment. However, the role of this circuit in working memory precision remains unknown.

OBJECTIVE

To test the role of the hippocampal-anterior thalamic pathway in working memory precision, we delivered intracranially electrical stimulation to the ANT. We hypothesize that ANT stimulation can improve working memory precision.

METHODS

Presurgical epilepsy patients with depth electrodes in ANT and hippocampus were recruited to perform a color-recall working memory task. Participants were instructed to point out the color they were supposed to recall by clicking a point on the color wheel, while the intracranial EEG data were synchronously recorded. For randomly selected half trials, a bipolar electrical stimulation was delivered to the ANT electrodes.

RESULTS

We found that compared to non-stimulation trials, working memory precision judgements were significantly improved for stimulation trials. ANT electrical stimulation significantly increased spectral power of gamma (30-100 Hz) oscillations and decreased interictal epileptiform discharges (IED) in the hippocampus. Moreover, the increased gamma power during the pre-stimulus and retrieval period predicted the improvement of working memory precision judgements.

CONCLUSION

ANT electrical stimulation can improve working memory precision judgements and modulate hippocampal gamma activity, providing direct evidence on the role of the human hippocampal-anterior thalamic axis in working memory precision.

Microscale Physiological Events on the Human Cortical Surface

Despite ongoing advances in our understanding of local single-cellular and network-level activity of neuronal populations in the human brain, extraordinarily little is known about their "intermediate" microscale local circuit dynamics. Here, we utilized ultra-high-density microelectrode arrays and a rare opportunity to perform intracranial recordings across multiple cortical areas in human participants to discover three distinct classes of cortical activity that are not locked to ongoing natural brain rhythmic activity. The first included fast waveforms similar to extracellular single-unit activity. The other two types were discrete events with slower waveform dynamics and were found preferentially in upper cortical layers. These second and third types were also observed in rodents, nonhuman primates, and semi-chronic recordings from humans via laminar and Utah array microelectrodes. The rates of all three events were selectively modulated by auditory and electrical stimuli, pharmacological manipulation, and cold saline application and had small causal co-occurrences. These results suggest that the proper combination of high-resolution microelectrodes and analytic techniques can capture neuronal dynamics that lay between somatic action potentials and aggregate population activity. Understanding intermediate microscale dynamics in relation to single-cell and network dynamics may reveal important details about activity in the full cortical circuit.

Interictal spike networks predict surgical outcome in patients with drug-resistant focal epilepsy

OBJECTIVE

To determine if properties of epileptic networks could be delineated using interictal spike propagation seen on stereo-electroencephalography (SEEG) and if these properties could predict surgical outcome in patients with drug-resistant epilepsy.

METHODS

We studied the SEEG of 45 consecutive drug-resistant epilepsy patients who underwent subsequent epilepsy surgery: 18 patients with good post-surgical outcome (Engel I) and 27 with poor outcome (Engel II-IV). Epileptic networks were derived from interictal spike propagation; these networks described the generation and propagation of interictal epileptic activity. We compared the regions in which spikes were frequent and the regions responsible for generating spikes to the area of resection and post-surgical outcome. We developed a measure termed source spike concordance, which integrates information about both spike rate and region of spike generation.

RESULTS

Inclusion in the resection of regions with high spike rate is associated with good post-surgical outcome (sensitivity = 0.82, specificity = 0.73). Inclusion in the resection of the regions responsible for generating interictal epileptic activity independently of rate is also associated with good post-surgical outcome (sensitivity = 0.88, specificity = 0.82). Finally, when integrating the spike rate and the generators, we find that the source spike concordance measure has strong predictability (sensitivity = 0.91, specificity = 0.94).

INTERPRETATIONS

Epileptic networks derived from interictal spikes can determine the generators of epileptic activity. Inclusion of the most active generators in the resection is strongly associated with good post-surgical outcome. These epileptic networks may aid clinicians in determining the area of resection during pre-surgical evaluation.

Rapid coordination of effective learning by the human hippocampus

Although the human hippocampus is necessary for long-term memory, controversial findings suggest that it may also support short-term memory in the service of guiding effective behaviors during learning. We tested the counterintuitive theory that the hippocampus contributes to long-term memory through remarkably short-term processing, as reflected in eye movements during scene encoding. While viewing scenes for the first time, short-term retrieval operative within the episode over only hundreds of milliseconds was indicated by a specific eye-movement pattern, which was effective in that it enhanced spatiotemporal memory formation. This viewing pattern was predicted by hippocampal theta oscillations recorded from depth electrodes and by shifts toward top-down influence of hippocampal theta on activity within visual perception and attention networks. The hippocampus thus supports short-term memory processing that coordinates behavior in the service of effective spatiotemporal learning.

Spatiotemporal dynamics between interictal epileptiform discharges and ripples during associative memory processing

We describe the spatiotemporal course of cortical high-gamma activity, hippocampal ripple activity and interictal epileptiform discharges during an associative memory task in 15 epilepsy patients undergoing invasive EEG. Successful encoding trials manifested significantly greater high-gamma activity in hippocampus and frontal regions. Successful cued recall trials manifested sustained high-gamma activity in hippocampus compared to failed responses. Hippocampal ripple rates were greater during successful encoding and retrieval trials. Interictal epileptiform discharges during encoding were associated with 15% decreased odds of remembering in hippocampus (95% confidence interval 6-23%). Hippocampal interictal epileptiform discharges during retrieval predicted 25% decreased odds of remembering (15-33%). Odds of remembering were reduced by 25-52% if interictal epileptiform discharges occurred during the 500-2000 ms window of encoding or by 41% during retrieval. During encoding and retrieval, hippocampal interictal epileptiform discharges were followed by a transient decrease in ripple rate. We hypothesize that interictal epileptiform discharges impair associative memory in a regionally and temporally specific manner by decreasing physiological hippocampal ripples necessary for effective encoding and recall. Because dynamic memory impairment arises from pathological interictal epileptiform discharge events competing with physiological ripples, interictal epileptiform discharges represent a promising therapeutic target for memory remediation in patients with epilepsy.