Quantitative analysis of high-frequency oscillations (80-500 Hz) recorded in human epileptic hippocampus and entorhinal cortex

A cortical field theory - dynamics and symmetries

We characterise cortical dynamics using partial differential equations (PDEs), analysing various connectivity patterns within the cortical sheet. This exploration yields diverse dynamics, encompassing wave equations and limit cycle activity. We presume balanced equations between excitatory and inhibitory neuronal units, reflecting the ubiquitous oscillatory patterns observed in electrophysiological measurements. Our derived dynamics comprise lowest-order wave equations (i.e., the Klein-Gordon model), limit cycle waves, higher-order PDE formulations, and transitions between limit cycles and near-zero states. Furthermore, we delve into the symmetries of the models using the Lagrangian formalism, distinguishing between continuous and discontinuous symmetries. These symmetries allow for mathematical expediency in the analysis of the model and could also be useful in studying the effect of symmetrical input from distributed cortical regions. Overall, our ability to derive multiple constraints on the fields - and predictions of the model - stems largely from the underlying assumption that the brain operates at a critical state. This assumption, in turn, drives the dynamics towards oscillatory or semi-conservative behaviour. Within this critical state, we can leverage results from the physics literature, which serve as analogues for neural fields, and implicit construct validity. Comparisons between our model predictions and electrophysiological findings from the literature - such as spectral power distribution across frequencies, wave propagation speed, epileptic seizure generation, and pattern formation over the cortical surface - demonstrate a close match. This study underscores the importance of utilizing symmetry preserving PDE formulations for further mechanistic insights into cortical activity.

Preventive effects of transcranial photobiomodulation on epileptogenesis in a kainic acid-induced rat epilepsy model

OBJECTIVE

Temporal lobe epilepsy affects nearly 50 million people worldwide and is a major burden to families and society. A significant portion of patients are living in developing countries with limited access to therapeutic resources. This highlights the urgent need to develop more readily available, noninvasive treatments for seizure control. This research explored the effectiveness of transcranial photobiomodulation (tPBM), a non-invasive method utilizing photon-tissue interactions, for preventing epileptogenesis and controlling seizures.

METHODS

In a kainic acid (KA)-induced rat model of epilepsy, two different wavelengths of tPBM, 808 nm and 940 nm, were applied separately in two groups of animals (KA+808 and KA+940). The ability of tPBM for seizure control was evaluated by comparing the occurrence rate of interictal epileptiform discharges (IED) and behavioral seizures among three groups: KA, KA+808, KA+940. Prevention of epileptogenesis was assessed by comparing the occurrence rate of high frequency oscillations (HFOs), especially fast ripple (FR) rate, among the three groups. Nissl staining and immunostaining for the apoptosis marker caspase-3 were used as indications of neuroprotection.

RESULTS

The KA+808 group and the KA+940 group showed significantly lower FR and IED rates compared to the KA group. Weekly FR rates started to drop during the first week of tPBM treatment. The KA+808 and KA+940 groups also displayed milder seizure behaviors and less neuronal loss in hippocampal areas compared to KA rats without tPBM treatment. Similarly, lower caspase-3 levels in the KA+808 and KA+940 compared with the KA group suggested effectiveness of tPBM in reducing cell death.

SIGNIFICANCE

tPBM of 808 nm/940 nm showed effectiveness in suppressing epileptogenesis and ictogenesis in the KA-induced rat epilepsy model. This effectiveness of tPBM can be linked to the neuroprotection benefits of photon-tissue interactions. Further studies are warranted to elucidate the fundamental mechanism of tPBM protection, determine optimal treatment parameters and validate its effectiveness in other epilepsy models.

Simulated resections and responsive neurostimulator placement can optimize postoperative seizure outcomes when guided by fast ripple networks

In medication-resistant epilepsy, the goal of epilepsy surgery is to make a patient seizure free with a resection/ablation that is as small as possible to minimize morbidity. The standard of care in planning the margins of epilepsy surgery involves electroclinical delineation of the seizure-onset zone and incorporation of neuroimaging findings from MRI, PET, single-photon emission CT and magnetoencephalography modalities. Resecting cortical tissue generating high-frequency oscillations has been investigated as a more efficacious alternative to targeting the seizure-onset zone. In this study, we used a support vector machine (SVM), with four distinct fast ripple (FR: 350-600 Hz on oscillations, 200-600 Hz on spikes) metrics as factors. These metrics included the FR resection ratio, a spatial FR network measure and two temporal FR network measures. The SVM was trained by the value of these four factors with respect to the actual resection boundaries and actual seizure-free labels of 18 patients with medically refractory focal epilepsy. Leave-one-out cross-validation of the trained SVM in this training set had an accuracy of 0.78. We next used a simulated iterative virtual resection targeting the FR sites that were of highest rate and showed most temporal autonomy. The trained SVM utilized the four virtual FR metrics to predict virtual seizure freedom. In all but one of the nine patients who were seizure free after surgery, we found that the virtual resections sufficient for virtual seizure freedom were larger in volume (P < 0.05). In nine patients who were not seizure free, a larger virtual resection made five virtually seizure free. We also examined 10 medically refractory focal epilepsy patients implanted with the responsive neurostimulator system and virtually targeted the responsive neurostimulator system stimulation contacts proximal to sites generating FR at highest rates to determine if the simulated value of the stimulated seizure-onset zone and stimulated FR metrics would trend towards those patients with a better seizure outcome. Our results suggest the following: (i) FR measures can accurately predict whether a resection, defined by the standard of care, will result in seizure freedom; (ii) utilizing FR alone for planning an efficacious surgery can be associated with larger resections; (iii) when FR metrics predict the standard-of-care resection will fail, amending the boundaries of the planned resection with certain FR-generating sites may improve outcome and (iv) more work is required to determine whether targeting responsive neurostimulator system stimulation contact proximal to FR generating sites will improve seizure outcome.

Overlap of spike and ripple propagation onset predicts surgical outcome in epilepsy

OBJECTIVE

Interictal biomarkers are critical for identifying the epileptogenic focus. However, spikes and ripples lack specificity while fast ripples lack sensitivity. These biomarkers propagate from more epileptogenic onset to areas of spread. The pathophysiological mechanism of these propagations is elusive. Here, we examine zones where spikes and high frequency oscillations co-occur (SHFO), the spatiotemporal propagations of spikes, ripples, and fast ripples, and evaluate the spike-ripple onset overlap (SRO) as an epilepsy biomarker.

METHODS

We retrospectively analyzed intracranial EEG data from 41 patients with drug-resistant epilepsy. We mapped propagations of spikes, ripples, and fast ripples, and identified their onset and spread zones, as well as SHFO and SRO. We then estimated the SRO prognostic value in predicting surgical outcome and compared it to onset and spread zones of spike, ripple, and fast ripple propagations, and SHFO.

RESULTS

We detected spikes and ripples in all patients and fast ripples in 12 patients (29%). We observed spike and ripple propagations in 40 (98%) patients. Spike and ripple onsets overlapped in 35 (85%) patients. In good outcome patients, SRO showed higher specificity and precision (p < 0.05) in predicting resection compared to onset and zones of spikes, ripples, and SHFO. Only SRO resection predicted outcome (p = 0.01) with positive and negative predictive values of 82% and 57%, respectively.

INTERPRETATION

SRO is a specific and precise biomarker of the epileptogenic zone whose removal predicts outcome. SRO is present in most patients with drug-resistant epilepsy. Such a biomarker may reduce prolonged intracranial monitoring and improve outcome.

Connectivity of high-frequency bursts as SOZ localization biomarker

For patients with refractory epilepsy, the seizure onset zone (SOZ) plays an essential role in determining the specific regions of the brain that will be surgically resected. High-frequency oscillations (HFOs) and connectivity-based approaches have been identified among the potential biomarkers to localize the SOZ. However, there is no consensus on how connectivity between HFO events should be estimated, nor on its subject-specific short-term reliability. Therefore, we propose the channel-level connectivity dispersion (CLCD) as a metric to quantify the variability in synchronization between individual electrodes and to identify clusters of electrodes with abnormal synchronization, which we hypothesize to be associated with the SOZ. In addition, we developed a specialized filtering method that reduces oscillatory components caused by filtering broadband artifacts, such as sharp transients, spikes, or direct current shifts. Our connectivity estimates are therefore robust to the presence of these waveforms. To calculate our metric, we start by creating binary signals indicating the presence of high-frequency bursts in each channel, from which we calculate the pairwise connectivity between channels. Then, the CLCD is calculated by combining the connectivity matrices and measuring the variability in each electrode's combined connectivity values. We test our method using two independent open-access datasets comprising intracranial electroencephalography signals from 89 to 15 patients with refractory epilepsy, respectively. Recordings in these datasets were sampled at approximately 1000 Hz, and our proposed CLCDs were estimated in the ripple band (80-200 Hz). Across all patients in the first dataset, the average ROC-AUC was 0.73, and the average Cohen's d was 1.05, while in the second dataset, the average ROC-AUC was 0.78 and Cohen's d was 1.07. On average, SOZ channels had lower CLCD values than non-SOZ channels. Furthermore, based on the second dataset, which includes surgical outcomes (Engel I-IV), our analysis suggested that higher CLCD interquartile (as a measure of CLCD distribution spread) is associated with favorable outcomes (Engel I). This suggests that CLCD could significantly assist in identifying SOZ clusters and, therefore, provide an additional tool in surgical planning for epilepsy patients.

Intracranial EEGs evidenced visual object processing in the human medial temporal lobe subregions

The perirhinal cortex (PRC) and parahippocampal cortex (PHC) are core regions along the visual dual-stream. The specific functional roles of the PRC and PHC and their interactions with the downstream hippocampus cortex (HPC) are crucial for understanding visual memory. Our research used human intracranial EEGs to study the neural mechanism of the PRC, PHC, and HPC in visual object encoding. Single-regional function analyses found evidence that the PRC, PHC, and HPC are activated ∼100 ms within the broad-gamma band and that the PRC was more strongly activated than either the PHC or the HPC after an object stimulus. Inter-regional analyses showed strong bidirectional interactions of the PRC with both the PHC and HPC in the low-frequency band, whereas the interactions between the PHC and HPC were not significant. These findings demonstrated the core role of the PRC in encoding visual object information and supported the hypothesis of PRC-HPC-ventral object pathway. The recruitment of the PHC and its interaction with the PRC in visual object encoding also provide new insights beyond the traditional dorsal-stream hypothesis.

High frequency oscillations in human memory and cognition: a neurophysiological substrate of engrams?

Despite advances in understanding the cellular and molecular processes underlying memory and cognition, and recent successful modulation of cognitive performance in brain disorders, the neurophysiological mechanisms remain underexplored. High frequency oscillations beyond the classic electroencephalogram spectrum have emerged as a potential neural correlate of fundamental cognitive processes. High frequency oscillations are detected in the human mesial temporal lobe and neocortical intracranial recordings spanning gamma/epsilon (60-150 Hz), ripple (80-250 Hz) and higher frequency ranges. Separate from other non-oscillatory activities, these brief electrophysiological oscillations of distinct duration, frequency and amplitude are thought to be generated by coordinated spiking of neuronal ensembles within volumes as small as a single cortical column. Although the exact origins, mechanisms and physiological roles in health and disease remain elusive, they have been associated with human memory consolidation and cognitive processing. Recent studies suggest their involvement in encoding and recall of episodic memory with a possible role in the formation and reactivation of memory traces. High frequency oscillations are detected during encoding, throughout maintenance, and right before recall of remembered items, meeting a basic definition for an engram activity. The temporal coordination of high frequency oscillations reactivated across cortical and subcortical neural networks is ideally suited for integrating multimodal memory representations, which can be replayed and consolidated during states of wakefulness and sleep. High frequency oscillations have been shown to reflect coordinated bursts of neuronal assembly firing and offer a promising substrate for tracking and modulation of the hypothetical electrophysiological engram.

Net synaptic drive of fast-spiking interneurons is inverted towards inhibition in human FCD I epilepsy

Focal cortical dysplasia type I (FCD I) is the most common cause of pharmaco-resistant epilepsy with the poorest prognosis. To understand the epileptogenic mechanisms of FCD I, we obtained tissue resected from patients with FCD I epilepsy, and from tumor patients as control. Using whole-cell patch clamp in acute human brain slices, we investigated the cellular properties of fast-spiking interneurons (FSINs) and pyramidal neurons (PNs) within the ictal onset zone. In FCD I epilepsy, FSINs exhibited lower firing rates from slower repolarization and action potential broadening, while PNs had increased firing. Importantly, excitatory synaptic drive of FSINs increased progressively with the scale of cortical activation as a general property across species, but this relationship was inverted towards net inhibition in FCD I epilepsy. Further comparison with intracranial electroencephalography (iEEG) from the same patients revealed that the spatial extent of pathological high-frequency oscillations (pHFO) was associated with synaptic events at FSINs.

The hemodynamic response to co-occurring interictal epileptiform discharges and high-frequency oscillations localizes the seizure-onset zone

OBJECTIVE

To use intracranial electroencephalography (EEG) to characterize functional magnetic resonance imaging (fMRI) activation maps associated with high-frequency oscillations (HFOs) (80-250 Hz) and examine their proximity to HFO- and seizure-generating tissue.

METHODS

Forty-five patients implanted with intracranial depth electrodes underwent a simultaneous EEG-fMRI study at 3 T. HFOs were detected algorithmically from cleaned EEG and visually confirmed by an experienced electroencephalographer. HFOs that co-occurred with interictal epileptiform discharges (IEDs) were subsequently identified. fMRI activation maps associated with HFOs were generated that occurred either independently of IEDs or within ±200 ms of an IED. For all significant analyses, the Maximum, Second Maximum, and Closest activation clusters were identified, and distances were measured to both the electrodes where the HFOs were observed and the electrodes involved in seizure onset.

RESULTS

We identified 108 distinct groups of HFOs from 45 patients. We found that HFOs with IEDs produced fMRI clusters that were closer to the local field potentials of the corresponding HFOs observed within the EEG than HFOs without IEDs. In addition to the fMRI clusters being closer to the location of the EEG correlate, HFOs with IEDs generated Maximum clusters with greater z-scores and larger volumes than HFOs without IEDs. We also observed that HFOs with IEDs resulted in more discrete activation maps.

SIGNIFICANCE

Intracranial EEG-fMRI can be used to probe the hemodynamic response to HFOs. The hemodynamic response associated with HFOs that co-occur with IEDs better identifies known epileptic tissue than HFOs that occur independently.

Optimizing automated detection of high frequency oscillations using visual markings does not improve SOZ localization

OBJECTIVE

High frequency oscillations (HFOs) are a biomarker of the seizure onset zone (SOZ) and can be visually or automatically detected. In theory, one can optimize an automated algorithm's parameters to maximize SOZ localization accuracy; however, there is no consensus on whether or how this should be done. Therefore, we optimized an automated detector using visually identified HFOs and evaluated the impact on SOZ localization accuracy.

METHODS

We detected HFOs in intracranial EEG from 20 patients with refractory epilepsy from two centers using (1) unoptimized automated detection, (2) visual identification, and (3) automated detection optimized to match visually detected HFOs.

RESULTS

SOZ localization accuracy based on HFO rate was not significantly different between the three methods. Across patients, visually optimized detector settings varied, and no single set of settings produced universally accurate SOZ localization. Exploratory analysis suggests that, for many patients, detection settings exist that would improve SOZ localization.

CONCLUSIONS

SOZ localization accuracy was similar for all three methods, was not improved by visually optimizing detector settings, and may benefit from patient-specific parameter optimization.

SIGNIFICANCE

Visual HFO marking is laborious, and optimizing automated detection using visual markings does not improve localization accuracy. New patient-specific detector optimization methods are needed.

Intan Technologies integrated circuits can produce analog-to-digital conversion artifacts that affect neural signal acquisition

Objective.Intan Technologies' integrated circuits (ICs) are valuable tools for neurophysiological data acquisition, providing signal amplification, filtering, and digitization from many channels (up to 64 channels/chip) at high sampling rates (up to 30 kSPS) within a compact package (⩽9× 7 mm). However, we found that the analog-to-digital converters (ADCs) in the Intan RHD2000 series ICs can produce artifacts in recorded signals. Here, we examine the effects of these ADC artifacts on neural signal quality and describe a method to detect them in recorded data.Approach.We identified two types of ADC artifacts produced by Intan ICs: 1) jumps, resulting from missing output codes, and 2) flatlines, resulting from overrepresented output codes. We identified ADC artifacts in neural recordings acquired with Intan RHD2000 ICs and tested the repeated performance of 17 ICsin vitro. With the on-chip digital-signal-processing disabled, we detected the ADC artifacts in each test recording by examining the distribution of unfiltered ADC output codes.Main Results.We found larger ADC artifacts in recordings using the Intan RHX data acquisition software versions 3.0-3.2, which did not run the necessary ADC calibration command when the inputs to the Intan recording controller were rescanned. This has been corrected in the Intan RHX software version 3.3. We found that the ADC calibration routine significantly reduced, but did not fully eliminate, the occurrence and size of ADC artifacts as compared with recordings acquired when the calibration routine was not run (p< 0.0001). When the ADC calibration routine was run, we found that the artifacts produced by each ADC were consistent over time, enabling us to sort ICs by performance.Significance.Our findings call attention to the importance of evaluating signal quality when acquiring electrophysiological data using Intan Technologies ICs and offer a method for detecting ADC artifacts in recorded data.

Spike ripples localize the epileptogenic zone best: an international intracranial study

We evaluated whether spike ripples, the combination of epileptiform spikes and ripples, provide a reliable and improved biomarker for the epileptogenic zone compared with other leading interictal biomarkers in a multicentre, international study. We first validated an automated spike ripple detector on intracranial EEG recordings. We then applied this detector to subjects from four centres who subsequently underwent surgical resection with known 1-year outcomes. We evaluated the spike ripple rate in subjects cured after resection [International League Against Epilepsy Class 1 outcome (ILAE 1)] and those with persistent seizures (ILAE 2-6) across sites and recording types. We also evaluated available interictal biomarkers: spike, spike-gamma, wideband high frequency oscillation (HFO, 80-500 Hz), ripple (80-250 Hz) and fast ripple (250-500 Hz) rates using previously validated automated detectors. The proportion of resected events was computed and compared across subject outcomes and biomarkers. Overall, 109 subjects were included. Most spike ripples were removed in subjects with ILAE 1 outcome (P < 0.001), and this was qualitatively observed across all sites and for depth and subdural electrodes (P < 0.001 and P < 0.001, respectively). Among ILAE 1 subjects, the mean spike ripple rate was higher in the resected volume (0.66/min) than in the non-removed tissue (0.08/min, P < 0.001). A higher proportion of spike ripples were removed in subjects with ILAE 1 outcomes compared with ILAE 2-6 outcomes (P = 0.06). Among ILAE 1 subjects, the proportion of spike ripples removed was higher than the proportion of spikes (P < 0.001), spike-gamma (P < 0.001), wideband HFOs (P < 0.001), ripples (P = 0.009) and fast ripples (P = 0.009) removed. At the individual level, more subjects with ILAE 1 outcomes had the majority of spike ripples removed (79%, 38/48) than spikes (69%, P = 0.12), spike-gamma (69%, P = 0.12), wideband HFOs (63%, P = 0.03), ripples (45%, P = 0.01) or fast ripples (36%, P < 0.001) removed. Thus, in this large, multicentre cohort, when surgical resection was successful, the majority of spike ripples were removed. Furthermore, automatically detected spike ripples localize the epileptogenic tissue better than spikes, spike-gamma, wideband HFOs, ripples and fast ripples.

Different Trajectories of Functional Connectivity Captured with Gamma-Event Coupling and Broadband Measures of EEG in the Rat Fluid Percussion Injury Model

Functional connectivity (FC) after TBI is affected by an altered excitatory-inhibitory balance due to neuronal dysfunction, and the mechanistic changes observed could be reflected differently by contrasting methods. Local gamma event coupling FC (GEC-FC) is believed to represent multiunit fluctuations due to inhibitory dysfunction, and we hypothesized that FC derived from widespread, broadband amplitude signal (BBA-FC) would be different, reflecting broader mechanisms of functional disconnection. We tested this during sleep and active periods defined by high delta and theta EEG activity, respectively, at 1,7 and 28d after rat fluid-percussion-injury (FPI) or sham injury (n=6/group) using 10 indwelling, bilateral cortical and hippocampal electrodes. We also measured seizure and high-frequency oscillatory activity (HFOs) as markers of electrophysiological burden. BBA-FC analysis showed early hyperconnectivity constrained to ipsilateral sensory-cortex-to-CA1-hippocampus that transformed to mainly ipsilateral FC deficits by 28d compared to shams. These changes were conserved over active epochs, except at 28d when there were no differences to shams. In comparison, GEC-FC analysis showed large regions of hyperconnectivity early after injury within similar ipsilateral and intrahemispheric networks. GEC-FC weakened with time, but hyperconnectivity persisted at 28d compared to sham. Edge- and global connectivity measures revealed injury-related differences across time in GEC-FC as compared to BBA-FC, demonstrating greater sensitivity to FC changes post-injury. There was no significant association between sleep fragmentation, HFOs, or seizures with FC changes. The within-animal, spatial-temporal differences in BBA-FC and GEC-FC after injury may represent different mechanisms driving FC changes as a result of primary disconnection and interneuron loss.

Significance statement

The present study adds to the understanding of functional connectivity changes in preclinical models of traumatic brain injury. In previously reported literature, there is heterogeneity in the directionality of connectivity changes after injury, resulting from factors such as severity of injury, frequency band studied, and methodology used to calculate FC. This study aims to further clarify differential mechanisms that result in altered network topography after injury, by using Broadband Amplitude-Derived FC and Gamma Event Coupling-Derived FC in EEG. We found post-injury changes that differ in complexity and directionality between measures at and across timepoints. In conjunction with known results and future studies identifying different neural drivers underlying these changes, measures derived from this study could provide useful means from which to minimally-invasively study temporally-evolving pathology after TBI.

PyHFO: lightweight deep learning-powered end-to-end high-frequency oscillations analysis application

Objective. This study aims to develop and validate an end-to-end software platform, PyHFO, that streamlines the application of deep learning (DL) methodologies in detecting neurophysiological biomarkers for epileptogenic zones from EEG recordings.Approach. We introduced PyHFO, which enables time-efficient high-frequency oscillation (HFO) detection algorithms like short-term energy and Montreal Neurological Institute and Hospital detectors. It incorporates DL models for artifact and HFO with spike classification, designed to operate efficiently on standard computer hardware.Main results. The validation of PyHFO was conducted on three separate datasets: the first comprised solely of grid/strip electrodes, the second a combination of grid/strip and depth electrodes, and the third derived from rodent studies, which sampled the neocortex and hippocampus using depth electrodes. PyHFO demonstrated an ability to handle datasets efficiently, with optimization techniques enabling it to achieve speeds up to 50 times faster than traditional HFO detection applications. Users have the flexibility to employ our pre-trained DL model or use their EEG data for custom model training.Significance. PyHFO successfully bridges the computational challenge faced in applying DL techniques to EEG data analysis in epilepsy studies, presenting a feasible solution for both clinical and research settings. By offering a user-friendly and computationally efficient platform, PyHFO paves the way for broader adoption of advanced EEG data analysis tools in clinical practice and fosters potential for large-scale research collaborations.

Automatic Detection of Scalp High-Frequency Oscillations Based on Deep Learning

Scalp high-frequency oscillations (sHFOs) are a promising non-invasive biomarker of epilepsy. However, the visual marking of sHFOs is a time-consuming and subjective process, existing automatic detectors based on single-dimensional analysis have difficulty with accurately eliminating artifacts and thus do not provide sufficient reliability to meet clinical needs. Therefore, we propose a high-performance sHFOs detector based on a deep learning algorithm. An initial detection module was designed to extract candidate high-frequency oscillations. Then, one-dimensional (1D) and two-dimensional (2D) deep learning models were designed, respectively. Finally, the weighted voting method is used to combine the outputs of the two model. In experiments, the precision, recall, specificity and F1-score were 83.44%, 83.60%, 96.61% and 83.42%, respectively, on average and the kappa coefficient was 80.02%. In addition, the proposed detector showed a stable performance on multi-centre datasets. Our sHFOs detector demonstrated high robustness and generalisation ability, which indicates its potential applicability as a clinical assistance tool. The proposed sHFOs detector achieves an accurate and robust method via deep learning algorithm.

Effects of midazolam on high-frequency oscillations in amygdala and hippocampus of epilepsy patients

High-frequency oscillations (HFOs) are associated with normal brain function, but are also increasingly recognized as potential biomarkers of epileptogenic tissue. Considering the important role of interneuron activity in physiological HFO generation, we studied their modulation by midazolam (MDZ), an agonist of γ-aminobutyric acid type A (GABAA)-benzodiazepine receptors. Here, we analyzed 80 intracranial electrode contacts in amygdala and hippocampus of 13 patients with drug-refractory focal epilepsy who had received MDZ for seizure termination during presurgical monitoring. Ripples (80-250 Hz) and fast ripples (FRs; 250-400 Hz) were compared before and after seizures with MDZ application, and according to their origin either within or outside the individual seizure onset zone (SOZ). We found that MDZ distinctly suppressed all HFOs (ripples and FRs), whereas the reduction of ripples was significantly less pronounced inside the SOZ compared to non-SOZ contacts. The rate of FRs inside the SOZ was less affected, especially in hippocampal contacts. In a few cases, even a marked increase of FRs following MDZ administration was seen. Our results demonstrate, for the first time, a significant HFO modulation in amygdala and hippocampus by MDZ, thus giving insights into the malfunction of GABA-mediated inhibition within epileptogenic areas and its role in HFO generation.

Photothrombosis induced cortical stroke produces electrographic epileptic biomarkers in mice

Objective

Interictal epileptiform spikes, high-frequency ripple oscillations, and their co-occurrence (spike ripples) in human scalp or intracranial voltage recordings are well-established epileptic biomarkers. While clinically significant, the neural mechanisms generating these electrographic biomarkers remain unclear. To reduce this knowledge gap, we introduce a novel photothrombotic stroke model in mice that reproduces focal interictal electrographic biomarkers observed in human epilepsy.

Methods

We induced a stroke in the motor cortex of C57BL/6 mice unilaterally (N=7) using a photothrombotic procedure previously established in rats. We then implanted intracranial electrodes (2 ipsilateral and 2 contralateral) and obtained intermittent local field potential (LFP) recordings over several weeks in awake, behaving mice. We evaluated the LFP for focal slowing and epileptic biomarkers - spikes, ripples, and spike ripples - using both automated and semi-automated procedures.

Results

Delta power (1-4 Hz) was higher in the stroke hemisphere than the non-stroke hemisphere in all mice ( p <0.001). Automated detection procedures indicated that compared to the non-stroke hemisphere, the stroke hemisphere had an increased spike ripple ( p =0.006) and spike rates ( p =0.039), but no change in ripple rate ( p =0.98). Expert validation confirmed the observation of elevated spike ripple rates ( p =0.008) and a trend of elevated spike rate ( p =0.055) in the stroke hemisphere. Interestingly, the validated ripple rate in the stroke hemisphere was higher than the non-stroke hemisphere ( p =0.031), highlighting the difficulty of automatically detecting ripples. Finally, using optimal performance thresholds, automatically detected spike ripples classified the stroke hemisphere with the best accuracy (sensitivity 0.94, specificity 0.94).

Significance

Cortical photothrombosis-induced stroke in commonly used C57BL/6 mice produces electrographic biomarkers as observed in human epilepsy. This model represents a new translational cortical epilepsy model with a defined irritative zone, which can be broadly applied in transgenic mice for cell type specific analysis of the cellular and circuit mechanisms of pathologic interictal activity.

Key Points

Cortical photothrombosis in mice produces stroke with characteristic intermittent focal delta slowing.Cortical photothrombosis stroke in mice produces the epileptic biomarkers spikes, ripples, and spike ripples.All biomarkers share morphological features with the corresponding human correlate.Spike ripples better lateralize to the lesional cortex than spikes or ripples.This cortical model can be applied in transgenic mice for mechanistic studies.

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.

High frequency oscillation network dynamics predict outcome in non-palliative epilepsy surgery

High frequency oscillations are a promising biomarker of outcome in intractable epilepsy. Prior high frequency oscillation work focused on counting high frequency oscillations on individual channels, and it is still unclear how to translate those results into clinical care. We show that high frequency oscillations arise as network discharges that have valuable properties as predictive biomarkers. Here, we develop a tool to predict patient outcome before surgical resection is performed, based on only prospective information. In addition to determining high frequency oscillation rate on every channel, we performed a correlational analysis to evaluate the functional connectivity of high frequency oscillations in 28 patients with intracranial electrodes. We found that high frequency oscillations were often not solitary events on a single channel, but part of a local network discharge. Eigenvector and outcloseness centrality were used to rank channel importance within the connectivity network, then used to compare patient outcome by comparison with the seizure onset zone or a proportion within the proposed resected channels (critical resection percentage). Combining the knowledge of each patient's seizure onset zone resection plan along with our computed high frequency oscillation network centralities and high frequency oscillation rate, we develop a Naïve Bayes model that predicts outcome (positive predictive value: 100%) better than predicting based upon fully resecting the seizure onset zone (positive predictive value: 71%). Surgical margins had a large effect on outcomes: non-palliative patients in whom most of the seizure onset zone was resected ('definitive surgery', ≥ 80% resected) had predictable outcomes, whereas palliative surgeries (<80% resected) were not predictable. These results suggest that the addition of network properties of high frequency oscillations is more accurate in predicting patient outcome than seizure onset zone alone in patients with most of the seizure onset zone removed and offer great promise for informing clinical decisions in surgery for refractory epilepsy.

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.

Sharp detection of oscillation packets in rich time-frequency representations of neural signals

Brain oscillations most often occur in bursts, called oscillation packets, which span a finite extent in time and frequency. Recent studies have shown that these packets portray a much more dynamic picture of synchronization and transient communication between sites than previously thought. To understand their nature and statistical properties, techniques are needed to objectively detect oscillation packets and to quantify their temporal and frequency extent, as well as their magnitude. There are various methods to detect bursts of oscillations. The simplest ones divide the signal into band limited sub-components, quantifying the strength of the resulting components. These methods cannot by themselves cope with broadband transients that look like genuine oscillations when restricted to a narrow band. The most successful detection methods rely on time-frequency representations, which can readily show broadband transients and harmonics. However, the performance of such methods is conditioned by the ability of the representation to localize packets simultaneously in time and frequency, and by the capabilities of packet detection techniques, whose current state of the art is limited to extraction of bounding boxes. Here, we focus on the second problem, introducing two detection methods that use concepts derived from clustering and topographic prominence. These methods are able to delineate the packets' precise contour in the time-frequency plane. We validate the new approaches using both synthetic and real data recorded in humans and animals and rely on a super-resolution time-frequency representation, namely the superlets, as input to the detection algorithms. In addition, we define robust tests for benchmarking and compare the new methods to previous techniques. Results indicate that the two methods we introduce shine in low signal-to-noise ratio conditions, where they only miss a fraction of packets undetected by previous methods. Finally, algorithms that delineate precisely the border of spectral features and their subcomponents offer far more valuable information than simple rectangular bounding boxes (time and frequency span) and can provide a solid foundation to investigate neural oscillations' dynamics.

Improving the accuracy of epileptogenic zone localization in stereo EEG with machine learning algorithms

The precise identification of the epileptogenic zone (EZ) is paramount in the presurgical evaluation of epilepsy patients to ensure successful surgical outcomes. The analysis of Stereo EEG, an instrumental tool for EZ localization, poses considerable challenges even for experienced epileptologists. Consequently, the development of machine learning (ML)-based computational tools for enhanced EZ localization is imperative. In this investigation, we developed ML models utilizing Stereo EEG from 15 patients, who remained seizure-free (Engel 1 a-d) following EZ resection, over an average follow-up period of 44.4 months. Utilizing Delphos and MNI detectors, spikes and High Frequency Oscillations (HFOs) were identified from Stereo EEG in Resected Zone (RZ) and non-Resected Zone (non-RZ). Linear and non-linear features were estimated from each modality using MATLAB. A total of 27,744 spikes, 7,790 ripples, and 7,632 fast ripples, along with their combinations, were employed to train the ML models. The Gradient Boosting classifier demonstrated the highest prediction accuracy of 98.5% for EZ localization in Mesial Temporal Lobe Epilepsy (MTLE) when trained with features derived from the spike-ripple combination. In the case of Neocortical Epilepsy (NE), the Extra Trees classifier achieved an accuracy of 87.6% when utilizing features from fast ripples. The Random Forest, Extra Trees, and Gradient Boosting algorithms were the most effective for predicting the RZ. Linear features outperformed non-linear features in predicting epileptogenic zones within the epileptic brain. Our study establishes the capability of ML methodologies in localizing epileptogenic zones with high accuracy. Future studies that focus on increasing the training sample size and incorporating more advanced machine learning (ML) algorithms have the potential to significantly improve the accuracy of these models in pinpointing epileptogenic networks. Additionally, implementing this ML approach across multiple research centers would contribute to the broader validation and generalizability of this technique.

Interictal invasive very high-frequency oscillations in resting awake state and sleep

Interictal very high-frequency oscillations (VHFOs, 500-2000 Hz) in a resting awake state seem to be, according to a precedent study of our team, a more specific predictor of a good outcome of the epilepsy surgery compared to traditional interictal high-frequency oscillations (HFOs, 80-500 Hz). In this study, we retested this hypothesis on a larger cohort of patients. In addition, we also collected patients' sleep data and hypothesized that the occurrence of VHFOs in sleep will be greater than in resting state. We recorded interictal invasive electroencephalographic (iEEG) oscillations in 104 patients with drug-resistant epilepsy in a resting state and in 35 patients during sleep. 21 patients in the rest study and 11 patients in the sleep study met the inclusion criteria (interictal HFOs and VHFOs present in iEEG recordings, a surgical intervention and a postoperative follow-up of at least 1 year) for further evaluation of iEEG data. In the rest study, patients with good postoperative outcomes had significantly higher ratio of resected contacts with VHFOs compared to HFOs. In sleep, VHFOs were more abundant than in rest and the percentage of resected contacts in patients with good and poor outcomes did not considerably differ in any type of oscillations. In conclusion, (1) our results confirm, in a larger patient cohort, our previous work about VHFOs being a specific predictor of the area which needs to be resected; and (2) that more frequent sleep VHFOs do not further improve the results.

Exploring Ripple Waves in the Human Brain

Ripples are brief (<150 ms) high-frequency oscillatory neural activities in the brain with a range of 140 to 200 Hz in rodents and 80 to 140 Hz in humans. Ripples are regarded as playing an essential role in several aspects of memory function, mainly in the hippocampus. This type of ripple generally occurs with sharp waves and is called a sharp-wave ripple (SPW-R). Extensive research of SPW-Rs in the rodent brain while actively awake has also linked the function of these SPW-Rs to navigation and decision making. Although many studies with rodents unveiled SPW-R function, research in humans on this subject is still sparse. Therefore, unveiling SPW-R function in the human hippocampus is warranted. A certain type of ripples may also be a biomarker of epilepsy. This type of ripple is called a pathological ripple (p-ripple). p-ripples have a wider range of frequency (80-500 Hz) than SPW-Rs, and the range of frequency is especially higher in brain regions that are intrinsically linked to epilepsy onset. Brain regions producing ripples are too small for scalp electrode recording, and intracranial recording is typically needed to detect ripples. In addition, SPW-Rs in the human hippocampus have been recorded from patients with epilepsy who may have p-ripples. Differentiating SPW-Rs and p-ripples is often not easy. We need to develop more sophisticated methods to record SPW-Rs to differentiate them from p-ripples. This paper reviews the general features and roles of ripple waves.

A scoping review of the relationship between alcohol, memory consolidation and ripple activity: An overview of common methodologies to analyse ripples

Alcohol abuse is not only responsible for 5.3% of the total deaths in the world but also has a substantial impact on neurological and memory disabilities throughout the population. One extensively studied brain area involved in cognitive functions is the hippocampus. Evidence in several rodent models has shown that ethanol produces cognitive impairment in hippocampal-dependent tasks and that the damage is varied according to the stage of development at which the rodent was exposed to ethanol and the dose. To the authors' knowledge, there is a biomarker for cognitive processes in the hippocampus that remains relatively understudied in association with memory impairment by alcohol administration. This biomarker is called sharp wave-ripples (SWRs) which are synchronous neuronal population events that are well known to be involved in memory consolidation. Methodologies for facilitated or automatic identification of ripples and their analysis have been reported for a wider bandwidth than SWRs. This review is focused on communicating the state of the art about the relationship between alcohol, memory consolidation and ripple activity, as well as the use of the common methodologies to identify SWRs automatically.

Developmental atlas of phase-amplitude coupling between physiologic high-frequency oscillations and slow waves

We investigated the developmental changes in high-frequency oscillation (HFO) and Modulation Index (MI) - the coupling measure between HFO and slow-wave phase. We generated normative brain atlases, using subdural EEG signals from 8251 nonepileptic electrode sites in 114 patients (ages 1.0-41.5 years) who achieved seizure control following resective epilepsy surgery. We observed a higher MI in the occipital lobe across all ages, and occipital MI increased notably during early childhood. The cortical areas exhibiting MI co-growth were connected via the vertical occipital fasciculi and posterior callosal fibers. While occipital HFO rate showed no significant age-association, the temporal, frontal, and parietal lobes exhibited an age-inversed HFO rate. Assessment of 1006 seizure onset sites revealed that z-score normalized MI and HFO rate were higher at seizure onset versus nonepileptic electrode sites. We have publicly shared our intracranial EEG data to enable investigators to validate MI and HFO-centric presurgical evaluations to identify the epileptogenic zone.

Optimizing detection and deep learning-based classification of pathological high-frequency oscillations in epilepsy

OBJECTIVE

This study aimed to explore sensitive detection methods for pathological high-frequency oscillations (HFOs) to improve seizure outcomes in epilepsy surgery.

METHODS

We analyzed interictal HFOs (80-500 Hz) in 15 children with medication-resistant focal epilepsy who underwent chronic intracranial electroencephalogram via subdural grids. The HFOs were assessed using the short-term energy (STE) and Montreal Neurological Institute (MNI) detectors and examined for spike association and time-frequency plot characteristics. A deep learning (DL)-based classification was applied to purify pathological HFOs. Postoperative seizure outcomes were correlated with HFO-resection ratios to determine the optimal HFO detection method.

RESULTS

The MNI detector identified a higher percentage of pathological HFOs than the STE detector, but some pathological HFOs were detected only by the STE detector. HFOs detected by both detectors had the highest spike association rate. The Union detector, which detects HFOs identified by either the MNI or STE detector, outperformed other detectors in predicting postoperative seizure outcomes using HFO-resection ratios before and after DL-based purification.

CONCLUSIONS

HFOs detected by standard automated detectors displayed different signal and morphological characteristics. DL-based classification effectively purified pathological HFOs.

SIGNIFICANCE

Enhancing the detection and classification methods of HFOs will improve their utility in predicting postoperative seizure outcomes.

Challenges of neural interfaces for stroke motor rehabilitation

More than 85% of stroke survivors suffer from different degrees of disability for the rest of their lives. They will require support that can vary from occasional to full time assistance. These conditions are also associated to an enormous economic impact for their families and health care systems. Current rehabilitation treatments have limited efficacy and their long-term effect is controversial. Here we review different challenges related to the design and development of neural interfaces for rehabilitative purposes. We analyze current bibliographic evidence of the effect of neuro-feedback in functional motor rehabilitation of stroke patients. We highlight the potential of these systems to reconnect brain and muscles. We also describe all aspects that should be taken into account to restore motor control. Our aim with this work is to help researchers designing interfaces that demonstrate and validate neuromodulation strategies to enforce a contingent and functional neural linkage between the central and the peripheral nervous system. We thus give clues to design systems that can improve or/and re-activate neuroplastic mechanisms and open a new recovery window for stroke patients.

Characterization of High-Gamma Activity in Electrocorticographic Signals

Introduction

Electrocorticographic (ECoG) high-gamma activity (HGA) is a widely recognized and robust neural correlate of cognition and behavior. However, fundamental signal properties of HGA, such as the high-gamma frequency band or temporal dynamics of HGA, have never been systematically characterized. As a result, HGA estimators are often poorly adjusted, such that they miss valuable physiological information.

Methods

To address these issues, we conducted a thorough qualitative and quantitative characterization of HGA in ECoG signals. Our study is based on ECoG signals recorded from 18 epilepsy patients while performing motor control, listening, and visual perception tasks. In this study, we first categorize HGA into HGA types based on the cognitive/behavioral task. For each HGA type, we then systematically quantify three fundamental signal properties of HGA: the high-gamma frequency band, the HGA bandwidth, and the temporal dynamics of HGA.

Results

The high-gamma frequency band strongly varies across subjects and across cognitive/behavioral tasks. In addition, HGA time courses have lowpass character, with transients limited to 10 Hz. The task-related rise time and duration of these HGA time courses depend on the individual subject and cognitive/behavioral task. Task-related HGA amplitudes are comparable across the investigated tasks.

Discussion

This study is of high practical relevance because it provides a systematic basis for optimizing experiment design, ECoG acquisition and processing, and HGA estimation. Our results reveal previously unknown characteristics of HGA, the physiological principles of which need to be investigated in further studies.

Etomidate-Induced Myoclonus in Sprague-Dawley Rats Involves Neocortical Glutamate Accumulation and N -Methyl- d -Aspartate Receptor Activity

BACKGROUND

Etomidate-induced myoclonus, a seizure-like movement, is of interest to anesthetists. However, its origin in the brain and its underlying mechanism remain unclear.

METHODS

Adult male Sprague-Dawley rats were anesthetized with etomidate, propofol, or lidocaine plus etomidate. We assessed the incidence of myoclonus, behavioral scores, and levels of glutamate and γ-aminobutyric acid (GABA) in the neocortex and hippocampus. To determine the origin and how N -methyl- d -aspartate receptors (NMDARs) modulate etomidate-induced neuroexcitability, the local field potential and muscular tension were monitored. Calcium imaging in vitro and immunoblotting in vivo were conducted to investigate the mechanisms underlying myoclonus.

RESULTS

The incidence of etomidate (1.5 mg/kg in vivo)-induced myoclonus was higher than that of propofol (90% vs 10%, P = .0010) and lidocaine plus etomidate (90% vs 20%, P = .0050). Etomidate at doses of 3.75 and 6 mg/kg decreased the mean behavioral score at 1 (mean difference [MD]: 1.80, 95% confidence interval [CI], 0.58-3.02; P = .0058 for both), 2 (MD: 1.60, 95% CI, 0.43-2.77; P = .0084 and MD: 1.70, 95% CI, 0.54-2.86; P = .0060), 3 (MD: 1.60, 95% CI, 0.35-2.85; P = .0127 and MD: 1.70, 95% CI, 0.46-2.94; P = .0091) minutes after administration compared to etomidate at a dose of 1.5 mg/kg. In addition, 0.5 and 1 µM etomidate in vitro increased neocortical intracellular calcium signaling; this signaling decreased when the concentration increased to 5 and 10 μM. Etomidate increased the glutamate level compared to propofol (mean rank difference: 18.20; P = .003), and lidocaine plus etomidate (mean rank difference: 21.70; P = .0002). Etomidate in vivo activated neocortical ripple waves and was positively correlated with muscular tension amplitude (Spearman's r = 0.785, P < .0001). Etomidate at 1.5 mg/kg decreased the K-Cl cotransporter isoform 2 (KCC2) level compared with propofol (MD: -1.15, 95% CI, -1.47 to -0.83; P < .0001) and lidocaine plus etomidate (MD: -0.64, 95% CI, -0.96 to -0.32; P = .0002), DL-2-amino-5-phosphopentanoic acid (AP5) suppressed these effects, while NMDA enhanced them.

CONCLUSIONS

Etomidate-induced myoclonus or neuroexcitability is concentration dependent. Etomidate-induced myoclonus originates in the neocortex. The underlying mechanism involves neocortical glutamate accumulation and NMDAR modulation and myoclonus correlates with NMDAR-induced downregulation of KCC2 protein expression.

Flexible, high-resolution cortical arrays with large coverage capture microscale high-frequency oscillations in patients with epilepsy

OBJECTIVE

Effective surgical treatment of drug-resistant epilepsy depends on accurate localization of the epileptogenic zone (EZ). High-frequency oscillations (HFOs) are potential biomarkers of the EZ. Previous research has shown that HFOs often occur within submillimeter areas of brain tissue and that the coarse spatial sampling of clinical intracranial electrode arrays may limit the accurate capture of HFO activity. In this study, we sought to characterize microscale HFO activity captured on thin, flexible microelectrocorticographic (μECoG) arrays, which provide high spatial resolution over large cortical surface areas.

METHODS

We used novel liquid crystal polymer thin-film μECoG arrays (.76-1.72-mm intercontact spacing) to capture HFOs in eight intraoperative recordings from seven patients with epilepsy. We identified ripple (80-250 Hz) and fast ripple (250-600 Hz) HFOs using a common energy thresholding detection algorithm along with two stages of artifact rejection. We visualized microscale subregions of HFO activity using spatial maps of HFO rate, signal-to-noise ratio, and mean peak frequency. We quantified the spatial extent of HFO events by measuring covariance between detected HFOs and surrounding activity. We also compared HFO detection rates on microcontacts to simulated macrocontacts by spatially averaging data.

RESULTS

We found visually delineable subregions of elevated HFO activity within each μECoG recording. Forty-seven percent of HFOs occurred on single 200-μm-diameter recording contacts, with minimal high-frequency activity on surrounding contacts. Other HFO events occurred across multiple contacts simultaneously, with covarying activity most often limited to a .95-mm radius. Through spatial averaging, we estimated that macrocontacts with 2-3-mm diameter would only capture 44% of the HFOs detected in our μECoG recordings.

SIGNIFICANCE

These results demonstrate that thin-film microcontact surface arrays with both highresolution and large coverage accurately capture microscale HFO activity and may improve the utility of HFOs to localize the EZ for treatment of drug-resistant epilepsy.

The role of superficial and deep layers in the generation of high frequency oscillations and interictal epileptiform discharges in the human cortex

Describing intracortical laminar organization of interictal epileptiform discharges (IED) and high frequency oscillations (HFOs), also known as ripples. Defining the frequency limits of slow and fast ripples. We recorded potential gradients with laminar multielectrode arrays (LME) for current source density (CSD) and multi-unit activity (MUA) analysis of interictal epileptiform discharges IEDs and HFOs in the neocortex and mesial temporal lobe of focal epilepsy patients. IEDs were observed in 20/29, while ripples only in 9/29 patients. Ripples were all detected within the seizure onset zone (SOZ). Compared to hippocampal HFOs, neocortical ripples proved to be longer, lower in frequency and amplitude, and presented non-uniform cycles. A subset of ripples (≈ 50%) co-occurred with IEDs, while IEDs were shown to contain variable high-frequency activity, even below HFO detection threshold. The limit between slow and fast ripples was defined at 150 Hz, while IEDs' high frequency components form clusters separated at 185 Hz. CSD analysis of IEDs and ripples revealed an alternating sink-source pair in the supragranular cortical layers, although fast ripple CSD appeared lower and engaged a wider cortical domain than slow ripples MUA analysis suggested a possible role of infragranularly located neural populations in ripple and IED generation. Laminar distribution of peak frequencies derived from HFOs and IEDs, respectively, showed that supragranular layers were dominated by slower (< 150 Hz) components. Our findings suggest that cortical slow ripples are generated primarily in upper layers while fast ripples and associated MUA in deeper layers. The dissociation of macro- and microdomains suggests that microelectrode recordings may be more selective for SOZ-linked ripples. We found a complex interplay between neural activity in the neocortical laminae during ripple and IED formation. We observed a potential leading role of cortical neurons in deeper layers, suggesting a refined utilization of LMEs in SOZ localization.

Optimizing Detection and Deep Learning-based Classification of Pathological High-Frequency Oscillations in Epilepsy

Objective

This study aimed to explore sensitive detection methods and deep learning (DL)-based classification for pathological high-frequency oscillations (HFOs).

Methods

We analyzed interictal HFOs (80-500 Hz) in 15 children with medication-resistant focal epilepsy who underwent resection after chronic intracranial electroencephalogram via subdural grids. The HFOs were assessed using the short-term energy (STE) and Montreal Neurological Institute (MNI) detectors and examined for pathological features based on spike association and time-frequency plot characteristics. A DL-based classification was applied to purify pathological HFOs. Postoperative seizure outcomes were correlated with HFO-resection ratios to determine the optimal HFO detection method.

Results

The MNI detector identified a higher percentage of pathological HFOs than the STE detector, but some pathological HFOs were detected only by the STE detector. HFOs detected by both detectors exhibited the most pathological features. The Union detector, which detects HFOs identified by either the MNI or STE detector, outperformed other detectors in predicting postoperative seizure outcomes using HFO-resection ratios before and after DL-based purification.

Conclusions

HFOs detected by standard automated detectors displayed different signal and morphological characteristics. DL-based classification effectively purified pathological HFOs.

Significance

Enhancing the detection and classification methods of HFOs will improve their utility in predicting postoperative seizure outcomes.

HIGHLIGHTS

HFOs detected by the MNI detector showed different traits and higher pathological bias than those detected by the STE detectorHFOs detected by both MNI and STE detectors (the Intersection HFOs) were deemed the most pathologicalA deep learning-based classification was able to distill pathological HFOs, regard-less of the initial HFO detection methods.

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.

Detection of pathological high-frequency oscillations in refractory epilepsy patients undergoing simultaneous stereo-electroencephalography and magnetoencephalography

BACKGROUND

Stereo-electroencephalography (SEEG) and magnetoencephalography (MEG) have generally been used independently as part of the pre-surgical evaluation of drug-resistant epilepsy (DRE) patients. However, the possibility of simultaneously employing these recording techniques to determine whether MEG has the potential of offering the same information as SEEG less invasively, or whether it could offer a greater spatial indication of the epileptogenic zone (EZ) to aid surgical planning, has not been previously evaluated.

METHODS

Data from 24 paediatric and adult DRE patients, undergoing simultaneous SEEG and MEG as part of their pre-surgical evaluation, was analysed employing manual and automated high-frequency oscillations (HFOs) detection, and spectral and source localisation analyses.

RESULTS

Twelve patients (50%) were included in the analysis (4 males; mean age=25.08 years) and showed interictal SEEG and MEG HFOs. HFOs detection was concordant between the two recording modalities, but SEEG displayed higher ability of differentiating between deep and superficial epileptogenic sources. Automated HFO detector in MEG recordings was validated against the manual MEG detection method. Spectral analysis revealed that SEEG and MEG detect distinct epileptic events. The EZ was well correlated with the simultaneously recorded data in 50% patients, while 25% patients displayed poor correlation or discordance.

CONCLUSION

MEG recordings can detect HFOs, and simultaneous use of SEEG and MEG HFO identification facilitates EZ localisation during the presurgical planning stage for DRE patients. Further studies are necessary to validate these findings and support the translation of automated HFO detectors into routine clinical practice.

Oscillations in the dentate gyrus as a tool for the performance of the hippocampal functions: Healthy and epileptic brain

The dentate gyrus (DG) is part of the hippocampal formation and is essential for important cognitive processes such as navigation and memory. The oscillatory activity of the DG network is believed to play a critical role in cognition. DG circuits generate theta, beta, and gamma rhythms, which participate in the specific information processing performed by DG neurons. In the temporal lobe epilepsy (TLE), cognitive abilities are impaired, which may be due to drastic alterations in the DG structure and network activity during epileptogenesis. The theta rhythm and theta coherence are especially vulnerable in dentate circuits; disturbances in DG theta oscillations and their coherence may be responsible for general cognitive impairments observed during epileptogenesis. Some researchers suggested that the vulnerability of DG mossy cells is a key factor in the genesis of TLE, but others did not support this hypothesis. The aim of the review is not only to present the current state of the art in this field of research but to help pave the way for future investigations by highlighting the gaps in our knowledge to completely appreciate the role of DG rhythms in brain functions. Disturbances in oscillatory activity of the DG during TLE development may be a diagnostic marker in the treatment of this disease.

Evolution of interictal activity in models of mesial temporal lobe epilepsy

Interictal activity and seizures are the hallmarks of focal epileptic disorders (which include mesial temporal lobe epilepsy, MTLE) in humans and in animal models. Interictal activity, which is recorded with cortical and intracerebral EEG recordings, comprises spikes, sharp waves and high-frequency oscillations, and has been used in clinical practice to identify the epileptic zone. However, its relation with seizures remains debated. Moreover, it is unclear whether specific EEG changes in interictal activity occur during the time preceding the appearance of spontaneous seizures. This period, which is termed "latent", has been studied in rodent models of MTLE in which spontaneous seizures start to occur following an initial insult (most often a status epilepticus induced by convulsive drugs such as kainic acid or pilocarpine) and may mirror epileptogenesis, i.e., the process leading the brain to develop an enduring predisposition to seizure generation. Here, we will address this topic by reviewing experimental studies performed in MTLE models. Specifically, we will review data highlighting the dynamic changes in interictal spiking activity and high-frequency oscillations occurring during the latent period, and how optogenetic stimulation of specific cell populations can modulate them in the pilocarpine model. These findings indicate that interictal activity: (i) is heterogeneous in its EEG patterns and thus, presumably, in its underlying neuronal mechanisms; and (ii) can pinpoint to the epileptogenic processes occurring in focal epileptic disorders in animal models and, perhaps, in epileptic patients.

Interictal epileptiform discharges affect memory in an Alzheimer's Disease mouse model

Interictal epileptiform discharges (IEDs) are transient abnormal electrophysiological events commonly observed in epilepsy patients but are also present in other neurological disease, such as Alzheimer's Disease (AD). Understanding the role IEDs have on the hippocampal circuit is important for our understanding of the cognitive deficits seen in epilepsy and AD. We characterize and compare the IEDs of human epilepsy patients from microwire hippocampal recording with those of AD transgenic mice with implanted multi-layer hippocampal silicon probes. Both the local field potential features and firing patterns of pyramidal cells and interneurons were similar in mouse and human. We found that as IEDs emerged from the CA3-1 circuits, they recruited pyramidal cells and silenced interneurons, followed by post-IED suppression. IEDs suppressed the incidence and altered the properties of physiological sharp-wave ripples (SPW-Rs), altered their physiological properties, and interfered with the replay of place field sequences in a maze. In addition, IEDs in AD mice inversely correlated with daily memory performance. Together, our work implicates that IEDs may present a common and epilepsy-independent phenomenon in neurodegenerative diseases that perturbs hippocampal-cortical communication and interferes with memory.

Significant Statement

Prevalence of neurodegenerative diseases and the number of people with dementia is increasing steadily. Therefore, novel treatment strategies for learning and memory disorders are urgently necessary. IEDs, apart from being a surrogate for epileptic brain regions, have also been linked to cognitive decline. Here we report that IEDs in human epilepsy patients and AD mouse models have similar local field potential characteristics and associated firing patterns of pyramidal cells and interneurons. Mice with more IEDs displayed fewer hippocampal SPW-Rs, poorer replay of spatial trajectories, and decreased memory performance. IED suppression is an unexplored target to treat cognitive dysfunction in neurodegenerative diseases.

Diverse nature of interictal oscillations: EEG-based biomarkers in epilepsy

Interictal electroencephalogram (EEG) patterns, including high-frequency oscillations (HFOs), interictal spikes (ISs), and slow wave activities (SWAs), are defined as specific oscillations between seizure events. These interictal oscillations reflect specific dynamic changes in network excitability and play various roles in epilepsy. In this review, we briefly describe the electrographic characteristics of HFOs, ISs, and SWAs in the interictal state, and discuss the underlying cellular and network mechanisms. We also summarize representative evidence from experimental and clinical epilepsy to address their critical roles in ictogenesis and epileptogenesis, indicating their potential as electrophysiological biomarkers of epilepsy. Importantly, we put forwards some perspectives for further research in the field.

Epileptic seizure focus detection from interictal electroencephalogram: a survey

Electroencephalogram (EEG) is one of most effective clinical diagnosis modalities for the localization of epileptic focus. Most current AI solutions use this modality to analyze the EEG signals in an automated manner to identify the epileptic seizure focus. To develop AI system for identifying the epileptic focus, there are many recently-published AI solutions based on biomarkers or statistic features that utilize interictal EEGs. In this review, we survey these solutions and find that they can be divided into three main categories: (i) those that use of biomarkers in EEG signals, including high-frequency oscillation, phase-amplitude coupling, and interictal epileptiform discharges, (ii) others that utilize feature-extraction methods, and (iii) solutions based upon neural networks (an end-to-end approach). We provide a detailed description of seizure focus with clinical diagnosis methods, a summary of the public datasets that seek to reduce the research gap in epilepsy, recent novel performance evaluation criteria used to evaluate the AI systems, and guidelines on when and how to use them. This review also suggests a number of future research challenges that must be overcome in order to design more efficient computer-aided solutions to epilepsy focus detection.

Hippocampal area CA2 controls seizure dynamics, interictal EEG abnormalities and social comorbidity in mouse models of temporal lobe epilepsy

Temporal lobe epilepsy (TLE) is characterized by spontaneous recurrent seizures, abnormal activity between seizures, and impaired behavior. CA2 pyramidal neurons (PNs) are potentially important because inhibiting them with a chemogenetic approach reduces seizure frequency in a mouse model of TLE. However, whether seizures could be stopped by timing inhibition just as a seizure begins is unclear. Furthermore, whether inhibition would reduce the cortical and motor manifestations of seizures are not clear. Finally, whether interictal EEG abnormalities and TLE comorbidities would be improved are unknown. Therefore, real-time optogenetic silencing of CA2 PNs during seizures, interictal activity and behavior were studied in 2 mouse models of TLE. CA2 silencing significantly reduced seizure duration and time spent in convulsive behavior. Interictal spikes and high frequency oscillations were significantly reduced, and social behavior was improved. Therefore, brief focal silencing of CA2 PNs reduces seizures, their propagation, and convulsive manifestations, improves interictal EEG, and ameliorates social comorbidities.

HIGHLIGHTS

Real-time CA2 silencing at the onset of seizures reduces seizure durationWhen CA2 silencing reduces seizure activity in hippocampus it also reduces cortical seizure activity and convulsive manifestations of seizuresInterictal spikes and high frequency oscillations are reduced by real-time CA2 silencingReal-time CA2 silencing of high frequency oscillations (>250Hz) rescues social memory deficits of chronic epileptic mice.

High-frequency oscillations (250-500 Hz) in animal models of Alzheimer's disease and two animal models of epilepsy

OBJECTIVE

To test the hypothesis that high-frequency oscillations (HFOs) between 250 and 500 Hz occur in mouse models of Alzheimer's disease (AD) and thus are not unique to epilepsy.

METHODS

Experiments were conducted in three mouse models of AD: Tg2576 mice that simulate a form of familial AD, presenilin 2 knock-out (PS2KO) mice, and the Ts65Dn model of Down's syndrome. We recorded HFOs using wideband (0.1-500 Hz, 2 kHz) intra-hippocampal and cortical surface electroencephalography (EEG) at 1 month until 24 months of age during wakefulness, slow wave sleep (SWS), and rapid eye movement (REM) sleep. In addition, interictal spikes (IISs) and seizures were analyzed for the possible presence of HFOs. Comparisons were made to the intra-hippocampal kainic acid and pilocarpine models of epilepsy.

RESULTS

We describe for the first time that hippocampal and cortical HFOs are a new EEG abnormality in AD mouse models. HFOs occurred in all transgenic mice but no controls. They were also detectable as early as 1 month of age and prior to amyloid beta plaque neuropathology. HFOs were most frequent during SWS (vs REM sleep or wakefulness). Notably, HFOs in the AD and epilepsy models were indistinguishable in both spectral frequency and duration. HFOs also occurred during IISs and seizures in the AD models, although with altered spectral properties compared to isolated HFOs.

SIGNIFICANCE

Our data demonstrate that HFOs, an epilepsy biomarker with high translational value, are not unique to epilepsy and thus not disease specific. Our findings also strengthen the idea of hyperexcitability in AD and its significant overlap with epilepsy. HFOs in AD mouse models may serve as an EEG biomarker, which is detectable from the scalp and thus amenable to noninvasive detection in people at risk for AD.

Intracranial electrophysiological recordings on a swine model of mesial temporal lobe epilepsy

Objective

To test the feasibility and reliability of intracranial electrophysiological recordings in an acute status epilepticus model on laboratory swine.

Method

Intrahippocampal injection of kainic acid (KA) was performed on 17 male Bama pigs (Sus scrofa domestica) weighing between 25 and 35 kg. Two stereoelectroencephalography (SEEG) electrodes with a total of 16 channels were implanted bilaterally along the sensorimotor cortex to the hippocampus. Brain electrical activity was recorded 2 h daily for 9-28 days. Three KA dosages were tested to evaluate the quantities capable of evoking status epilepticus. Local field potentials (LFPs) were recorded and compared before and after the KA injection. We quantified the epileptic patterns, including the interictal spikes, seizures, and high-frequency oscillations (HFOs), up to 4 weeks after the KA injection. Test-retest reliability using intraclass correlation coefficients (ICCs) were performed on interictal HFO rates to evaluate the recording stability of this model.

Results

The KA dosage test suggested that a 10 μl (1.0 μg/μl) intrahippocampal injection could successfully evoke status epilepticus lasting from 4 to 12 h. At this dosage, eight pigs (50% of total) had prolonged epileptic events (tonic-chronic seizures + interictal spikes n = 5, interictal spikes alone n = 3) in the later 4 weeks of the video-SEEG recording period. Four pigs (25% of total) had no epileptic activities, and another four (25%) had lost the cap or did not complete the experiments. Animals that showed epileptiform events were grouped as E + (n = 8) and the four animals showing no signs of epileptic events were grouped as E- (n = 4). A total of 46 electrophysiological seizures were captured in the 4-week post-KA period from 4 E + animals, with the earliest onset on day 9. The seizure durations ranged from 12 to 45 s. A significant increase of hippocampal HFOs rate (num/min) was observed in the E+ group during the post-KA period (weeks 1, 2,4, p < 0.05) compared to the baseline. But the E-showed no change or a decrease (in week 2, p = 0.43) compared to their baseline rate. The between-group comparison showed much higher HFO rates in E + vs. E - (F = 35, p < 0.01). The high ICC value [ICC (1, k) = 0.81, p < 0.05] quantified from the HFO rate suggested that this model had a stable measurement of HFOs during the four-week post-KA periods.

Significance

This study measured intracranial electrophysiological activity in a swine model of KA-induced mesial temporal lobe epilepsy (mTLE). Using the clinical SEEG electrode, we distinguished abnormal EEG patterns in the swine brain. The high test-retest reliability of HFO rates in the post-KA period suggests the utility of this model for studying mechanisms of epileptogenesis. The use of swine may provide satisfactory translational value for clinical epilepsy research.

Processing of cell assemblies in the lateral entorhinal cortex

There is evidence that olfactory cortex responds to its afferent input with the generation of cell assemblies: collections of principal neurons that fire together over a time scale of tens of ms. If such assemblies form an odor representation, then a fundamental question is how each assembly then induces neuronal activity in downstream structures. We have addressed this question in a detailed model of superficial layers of lateral entorhinal cortex, a recipient of input from olfactory cortex and olfactory bulb. Our results predict that the response of the fan cell subpopulation can be approximated by a relatively simple Boolean process, somewhat along the lines of the McCulloch/Pitts scheme; this is the case because of the sparsity of recurrent excitation amongst fan cells. However, because of recurrent excitatory connections between layer 2 and layer 3 pyramidal cells, synaptic and probably also gap junctional, the response of pyramidal cell subnetworks cannot be so approximated. Because of the highly structured anatomy of entorhinal output projections, our model suggests that downstream targets of entorhinal cortex (dentate gyrus, hippocampal CA3, CA1, piriform cortex, olfactory bulb) receive differentially processed information.

Characterizing physiological high-frequency oscillations using deep learning

Objective.Intracranially-recorded interictal high-frequency oscillations (HFOs) have been proposed as a promising spatial biomarker of the epileptogenic zone. However, HFOs can also be recorded in the healthy brain regions, which complicates the interpretation of HFOs. The present study aimed to characterize salient features of physiological HFOs using deep learning (DL).Approach.We studied children with neocortical epilepsy who underwent intracranial strip/grid evaluation. Time-series EEG data were transformed into DL training inputs. The eloquent cortex (EC) was defined by functional cortical mapping and used as a DL label. Morphological characteristics of HFOs obtained from EC (ecHFOs) were distilled and interpreted through a novel weakly supervised DL model.Main results.A total of 63 379 interictal intracranially-recorded HFOs from 18 children were analyzed. The ecHFOs had lower amplitude throughout the 80-500 Hz frequency band around the HFO onset and also had a lower signal amplitude in the low frequency band throughout a one-second time window than non-ecHFOs, resembling a bell-shaped template in the time-frequency map. A minority of ecHFOs were HFOs with spikes (22.9%). Such morphological characteristics were confirmed to influence DL model prediction via perturbation analyses. Using the resection ratio (removed HFOs/detected HFOs) of non-ecHFOs, the prediction of postoperative seizure outcomes improved compared to using uncorrected HFOs (area under the ROC curve of 0.82, increased from 0.76).Significance.We characterized salient features of physiological HFOs using a DL algorithm. Our results suggested that this DL-based HFO classification, once trained, might help separate physiological from pathological HFOs, and efficiently guide surgical resection using HFOs.

Source-sink connectivity: a novel interictal EEG marker for seizure localization

Over 15 million epilepsy patients worldwide have drug-resistant epilepsy. Successful surgery is a standard of care treatment but can only be achieved through complete resection or disconnection of the epileptogenic zone, the brain region(s) where seizures originate. Surgical success rates vary between 20% and 80%, because no clinically validated biological markers of the epileptogenic zone exist. Localizing the epileptogenic zone is a costly and time-consuming process, which often requires days to weeks of intracranial EEG (iEEG) monitoring. Clinicians visually inspect iEEG data to identify abnormal activity on individual channels occurring immediately before seizures or spikes that occur interictally (i.e. between seizures). In the end, the clinical standard mainly relies on a small proportion of the iEEG data captured to assist in epileptogenic zone localization (minutes of seizure data versus days of recordings), missing opportunities to leverage these largely ignored interictal data to better diagnose and treat patients. IEEG offers a unique opportunity to observe epileptic cortical network dynamics but waiting for seizures increases patient risks associated with invasive monitoring. In this study, we aimed to leverage interictal iEEG data by developing a new network-based interictal iEEG marker of the epileptogenic zone. We hypothesized that when a patient is not clinically seizing, it is because the epileptogenic zone is inhibited by other regions. We developed an algorithm that identifies two groups of nodes from the interictal iEEG network: those that are continuously inhibiting a set of neighbouring nodes ('sources') and the inhibited nodes themselves ('sinks'). Specifically, patient-specific dynamical network models were estimated from minutes of iEEG and their connectivity properties revealed top sources and sinks in the network, with each node being quantified by source-sink metrics. We validated the algorithm in a retrospective analysis of 65 patients. The source-sink metrics identified epileptogenic regions with 73% accuracy and clinicians agreed with the algorithm in 93% of seizure-free patients. The algorithm was further validated by using the metrics of the annotated epileptogenic zone to predict surgical outcomes. The source-sink metrics predicted outcomes with an accuracy of 79% compared to an accuracy of 43% for clinicians' predictions (surgical success rate of this dataset). In failed outcomes, we identified brain regions with high metrics that were untreated. When compared with high frequency oscillations, the most commonly proposed interictal iEEG feature for epileptogenic zone localization, source-sink metrics outperformed in predictive power (by a factor of 1.2), suggesting they may be an interictal iEEG fingerprint of the epileptogenic zone.

Reduction of Hippocampal High-Frequency Activity in Wag/Rij Rats with a Genetic Predisposition to Absence Epilepsy

In temporal lobe epilepsy, high frequency oscillations serve as electroencephalographic (EEG) markers of epileptic hippocampal tissue. In contrast, absence epilepsy and other idiopathic epilepsies are known to result from thalamo-cortical abnormalities, with the hippocampus involvement considered to be only indirect. We aimed to uncover the role of the hippocampus in absence epilepsy using a genetic rat model of absence epilepsy (WAG/Rij rats), in which spike-wave discharges (SWDs) appear spontaneously in cortical EEG. We performed simultaneous recordings of local field potential from the hippocampal dentate gyrus using pairs of depth electrodes and epidural cortical EEG in freely moving rats. Hippocampal ripples (100-200 Hz) and high frequency oscillations (HFO, 50-70 Hz) were detected using GUI RIPPLELAB in MatLab (Navarrete et al., 2016). Based on the dynamics of hippocampal ripples, SWDs were divided into three clusters, which might represent different seizure types in reference to the involvement of hippocampal processes. This might underlie impairment of hippocampus-related cognitive processes in some patients with absence epilepsy. A significant reduction to nearly zero-ripple-density was found 4-8 s prior to SWD onset and during 4 s immediately after SWD onset. It follows that hippocampal ripples were not just passively blocked by the onset of SWDs, but they were affected by spike-wave seizure initiation mechanisms. Hippocampal HFO were reduced during the preictal, ictal and postictal periods in comparison to the baseline. Therefore, hippocampal HFO seemed to be blocked with spike-wave seizures. All together, this might underlie impairment of hippocampus-related cognitive processes in some patients with absence epilepsy. Further investigation of processes underlying SWD-related reduction of hippocampal ripples and HFO oscillations may help to predict epileptic attacks and explain cognitive comorbidities in patients with absence epilepsy.

Finger movement and coactivation predicted from intracranial brain activity using extended block-term tensor regression

Objective.We introduce extended Block-Term Tensor Regression (eBTTR), a novel regression method designed to account for the multilinear nature of human intracranial finger movement recordings.Approach.The proposed method relies on recursive Tucker decomposition combined with automatic component extraction.Main results.eBTTR outperforms state-of-the-art regression approaches, including multilinear and deep learning ones, in accurately predicting finger trajectories as well as unintentional finger coactivations.Significance.eBTTR rivals state-of-the-art approaches while being less computationally expensive which is an advantage when intracranial electrodes are implanted acutely, as part of the patient's presurgical workup, limiting time for decoder development and testing.

Multi-Scale Deep Learning of Clinically Acquired Multi-Modal MRI Improves the Localization of Seizure Onset Zone in Children With Drug-Resistant Epilepsy

The present study investigates the effectiveness of a deep learning neural network for non-invasively localizing the seizure onset zone (SOZ) using multi-modal MRI data that are clinically acquired from children with drug-resistant epilepsy. A cortical parcellation was applied to localize the SOZ in cortical nodes of the epileptogenic hemisphere. At each node, the laminar surface analysis was followed to sample 1) the relative intensity of gray matter and white matter in multi-modal MRI and 2) the neighboring white matter connectivity using diffusion tractography edge strengths. A cross-validation was employed to train and test all layers of a multi-scale residual neural network (msResNet) that can classify SOZ node in an end-to-end fashion. A prediction probability of a given node belonging to the SOZ class was proposed as a non-invasive MRI marker of seizure onset likelihood. In an independent validation cohort, the proposed MRI marker provided a very large effect size of Cohen's d = 1.21 between SOZ and non-SOZ, and classified SOZ with a balanced accuracy of 0.75 in lesional and 0.67 in non-lesional MRI groups. The subsequent multi-variate logistic regression found the incorporation of the proposed MRI marker into interictal intracranial EEG (iEEG) markers further improves the differentiation between the epileptogenic focus (defined as SOZ resected during surgery) and non-epileptogenic sites (i.e., non-SOZ sites preserved during surgery) up to 15 % in non-lesional MRI group, suggesting that the proposed MRI marker could improve the localization of epileptogenic foci for successful pediatric epilepsy surgery.

Tracking Multisite Seizure Propagation Using Ictal High-Gamma Activity

PURPOSE

Spatial patterns of long-range seizure propagation in epileptic networks have not been well characterized. Here, we use ictal high-gamma activity (HGA) as a proxy of intense neuronal population firing to map the spatial evolution of seizure recruitment.

METHODS

Ictal HGA (80-150 Hz) was analyzed in 13 patients with 72 seizures recorded by stereotactic depth electrodes, using previously validated methods. Distinct spatial clusters of channels with the ictal high-gamma signature were identified, and seizure hubs were defined as stereotypically recruited nonoverlapping clusters. Clusters correlated with asynchronous seizure terminations to provide supportive evidence for independent seizure activity at these sites. The spatial overlap between seizure hubs and interictal ripples was compared.

RESULTS

Ictal HGA was detected in 71% of seizures and 10% of implanted contacts, enabling tracking of contiguous and noncontiguous seizure recruitment. Multiple seizure hubs were identified in 54% of cases, including 43% of patients thought preoperatively to have unifocal epilepsy. Noncontiguous recruitment was associated with asynchronous seizure termination (odds ratio = 19.7; p = 0.029). Interictal ripples demonstrated greater spatial overlap with ictal HGA in cases with single seizure hubs compared with those with multiple hubs (100% vs. 66% per patient; p = 0.03).

CONCLUSIONS

Ictal HGA may serve as a useful adjunctive biomarker to distinguish contiguous seizure spread from propagation to remote seizure sites. High-gamma sites were found to cluster in stereotyped seizure hubs rather than being broadly distributed. Multiple hubs were common even in cases that were considered unifocal.