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.

Simulated resections and RNS placement can optimize post-operative seizure outcomes when guided by fast ripple networks.

In medication-resistant epilepsy, the objective of epilepsy surgery is to render 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 (SOZ) in the epilepsy monitoring unit (EMU) by an expert epileptologist, and incorporation of neuroimaging findings from MRI, PET, SPECT, and MEG modalities. Resecting cortical tissue generating high-frequency oscillations (HFOs) has been investigated as a more efficacious alternative to targeting the SOZ. In this study, we used support vector machines (SVMs) to compare the resection volumes of actual resections, defined using the clinical standard of care, with virtual resections of fast ripple (FR 350-600 Hz) networks. Cross-validation of the SVM that labeled patients as seizure free or not seizure free using FR metrics as factors demonstrated an accuracy of 0.78. In all the patients rendered seizure free after surgery, we found that the virtual resection, defined using FR generator sites with highest rate and greatest autonomy, was larger than the actual resection. In the patients who experienced seizures after the actual resection, a virtual resection that included the SOZ and other FR generating regions rendered half of these patients virtually seizure free. We also examined patients implanted with the responsive neurostimulator system (RNS) and virtually targeted the RNS stimulation contacts proximal to sites generating FR. We used the simulations to investigate if the likelihood of a RNS super response (>90% seizure reduction) would be increased.

ComBating inter-site differences in field strength: harmonizing preclinical traumatic brain injury MRI data.

Integrating datasets from multiple sites and scanners can increase statistical power for neuroimaging studies but can also introduce significant inter-site confounds. We evaluated the effectiveness of ComBat, an empirical Bayes approach, to combine longitudinal preclinical MRI data acquired at 4.7 or 9.4 T at two different sites in Australia. Male Sprague Dawley rats underwent MRI on Days 2, 9, 28, and 150 following moderate/severe traumatic brain injury (TBI) or sham injury as part of Project 1 of the NIH/NINDS-funded Centre Without Walls EpiBioS4Rx project. Diffusion-weighted and multiple-gradient-echo images were acquired, and outcomes included QSM, FA, and ADC. Acute injury measures including apnea and self-righting reflex were consistent between sites. Mixed-effect analysis of ipsilateral and contralateral corpus callosum (CC) summary values revealed a significant effect of site on FA and ADC values, which was removed following ComBat harmonization. Bland-Altman plots for each metric showed reduced variability across sites following ComBat harmonization, including for QSM, despite appearing to be largely unaffected by inter-site differences and no effect of site observed. Following harmonization, the combined inter-site data revealed significant differences in the imaging metrics consistent with previously reported outcomes. TBI resulted in significantly reduced FA and increased susceptibility in the ipsilateral CC, and significantly reduced FA in the contralateral CC compared with sham-injured rats. Additionally, TBI rats also exhibited a reversal in ipsilateral CC ADC values over time with significantly reduced ADC at Day 9, followed by increased ADC 150 days after injury. Our findings demonstrate the need for harmonizing multi-site preclinical MRI data and show that this can be successfully achieved using ComBat while preserving phenotypical changes due to TBI.

Patterns of subregional cerebellar atrophy across epilepsy syndromes: An ENIGMA-Epilepsy study.

OBJECTIVE: The intricate neuroanatomical structure of the cerebellum is of longstanding interest in epilepsy, but has been poorly characterized within the current corticocentric models of this disease. We quantified cross-sectional regional cerebellar lobule volumes using structural magnetic resonance imaging in 1602 adults with epilepsy and 1022 healthy controls across 22 sites from the global ENIGMA-Epilepsy working group. METHODS: A state-of-the-art deep learning-based approach was employed that parcellates the cerebellum into 28 neuroanatomical subregions. Linear mixed models compared total and regional cerebellar volume in (1) all epilepsies, (2) temporal lobe epilepsy with hippocampal sclerosis (TLE-HS), (3) nonlesional temporal lobe epilepsy, (4) genetic generalized epilepsy, and (5) extratemporal focal epilepsy (ETLE). Relationships were examined for cerebellar volume versus age at seizure onset, duration of epilepsy, phenytoin treatment, and cerebral cortical thickness. RESULTS: Across all epilepsies, reduced total cerebellar volume was observed (d = .42). Maximum volume loss was observed in the corpus medullare (dmax = .49) and posterior lobe gray matter regions, including bilateral lobules VIIB (dmax = .47), crus I/II (dmax = .39), VIIIA (dmax = .45), and VIIIB (dmax = .40). Earlier age at seizure onset ( η ρ max 2 = .05) and longer epilepsy duration ( η ρ max 2 = .06) correlated with reduced volume in these regions. Findings were most pronounced in TLE-HS and ETLE, with distinct neuroanatomical profiles observed in the posterior lobe. Phenytoin treatment was associated with reduced posterior lobe volume. Cerebellum volume correlated with cerebral cortical thinning more strongly in the epilepsy cohort than in controls. SIGNIFICANCE: We provide robust evidence of deep cerebellar and posterior lobe subregional gray matter volume loss in patients with chronic epilepsy. Volume loss was maximal for posterior subregions implicated in nonmotor functions, relative to motor regions of both the anterior and posterior lobe. Associations between cerebral and cerebellar changes, and variability of neuroanatomical profiles across epilepsy syndromes argue for more precise incorporation of cerebellar subregional damage into neurobiological models of epilepsy.

Epilepsy phenotype and its reproducibility after lateral fluid percussion-induced traumatic brain injury in rats: Multicenter EpiBioS4Rx study project 1.

OBJECTIVE: This study was undertaken to assess reproducibility of the epilepsy outcome and phenotype in a lateral fluid percussion model of posttraumatic epilepsy (PTE) across three study sites. METHODS: A total of 525 adult male Sprague Dawley rats were randomized to lateral fluid percussion-induced brain injury (FPI) or sham operation. Of these, 264 were assigned to magnetic resonance imaging (MRI cohort, 43 sham, 221 traumatic brain injury [TBI]) and 261 to electrophysiological follow-up (EEG cohort, 41 sham, 220 TBI). A major effort was made to harmonize the rats, materials, equipment, procedures, and monitoring systems. On the 7th post-TBI month, rats were video-EEG monitored for epilepsy diagnosis. RESULTS: A total of 245 rats were video-EEG phenotyped for epilepsy on the 7th postinjury month (121 in MRI cohort, 124 in EEG cohort). In the whole cohort (n = 245), the prevalence of PTE in rats with TBI was 22%, being 27% in the MRI and 18% in the EEG cohort (p > .05). Prevalence of PTE did not differ between the three study sites (p > .05). The average seizure frequency was .317 ± .725 seizures/day at University of Eastern Finland (UEF; Finland), .085 ± .067 at Monash University (Monash; Australia), and .299 ± .266 at University of California, Los Angeles (UCLA; USA; p < .01 as compared to Monash). The average seizure duration did not differ between UEF (104 ± 48 s), Monash (90 ± 33 s), and UCLA (105 ± 473 s; p > .05). Of the 219 seizures, 53% occurred as part of a seizure cluster (≥3 seizures/24 h; p >.05 between the study sites). Of the 209 seizures, 56% occurred during lights-on period and 44% during lights-off period (p > .05 between the study sites). SIGNIFICANCE: The PTE phenotype induced by lateral FPI is reproducible in a multicenter design. Our study supports the feasibility of performing preclinical multicenter trials in PTE to increase statistical power and experimental rigor to produce clinically translatable data to combat epileptogenesis after TBI.

Successful harmonization in EpiBioS4Rx biomarker study on post-traumatic epilepsy paves the way towards powered preclinical multicenter studies.

OBJECTIVE: Project 1 of the Preclinical Multicenter Epilepsy Bioinformatics Study for Antiepileptogenic Therapy (EpiBioS4Rx) consortium aims to identify preclinical biomarkers for antiepileptogenic therapies following traumatic brain injury (TBI). The international participating centers in Finland, Australia, and the United States have made a concerted effort to ensure protocol harmonization. Here, we evaluate the success of harmonization process by assessing the timing, coverage, and performance between the study sites. METHOD: We collected data on animal housing conditions, lateral fluid-percussion injury model production, postoperative care, mortality, post-TBI physiological monitoring, timing of blood sampling and quality, MR imaging timing and protocols, and duration of video-electroencephalography (EEG) follow-up using common data elements. Learning effect in harmonization was assessed by comparing procedural accuracy between the early and late stages of the project. RESULTS: The animal housing conditions were comparable between the study sites but the postoperative care procedures varied. Impact pressure, duration of apnea, righting reflex, and acute mortality differed between the study sites (p < 0.001). The severity of TBI on D2 post TBI assessed using the composite neuroscore test was similar between the sites, but recovery of acute somato-motor deficits varied (p < 0.001). A total of 99% of rats included in the final cohort in UEF, 100% in Monash, and 79% in UCLA had blood samples taken at all time points. The timing of sampling differed on day (D)2 (p < 0.05) but not D9 (p > 0.05). Plasma quality was poor in 4% of the samples in UEF, 1% in Monash and 14% in UCLA. More than 97% of the final cohort were MR imaged at all timepoints in all study sites. The timing of imaging did not differ on D2 and D9 (p > 0.05), but varied at D30, 5 months, and ex vivo timepoints (p < 0.001). The percentage of rats that completed the monthly high-density video-EEG follow-up and the duration of video-EEG recording on the 7th post-injury month used for seizure detection for diagnosis of post-traumatic epilepsy differed between the sites (p < 0.001), yet the prevalence of PTE (UEF 21%, Monash 22%, UCLA 23%) was comparable between the sites (p > 0.05). A decrease in acute mortality and increase in plasma quality across time reflected a learning effect in the TBI production and blood sampling protocols. SIGNIFICANCE: Our study is the first demonstration of the feasibility of protocol harmonization for performing powered preclinical multi-center trials for biomarker and therapy discovery of post-traumatic epilepsy.

Pathological neurons generate ripples at the UP-DOWN transition disrupting information transfer.

OBJECTIVE: To confirm and investigate why pathological high-frequency oscillations (pHFOs), including ripples (80-200 Hz) and fast ripples (200-600 Hz), are generated during the UP-DOWN transition of the slow wave and if information transmission mediated by ripple temporal coupling is disrupted in the seizure-onset zone (SOZ). METHODS: We isolated 217 total units from 175.95 intracranial electroencephalography (iEEG) contact-hours of synchronized macro- and microelectrode recordings from 6 patients. Sleep slow oscillation (.1-2 Hz) epochs were identified in the iEEG recording. iEEG HFOs that occurred superimposed on the slow wave were transformed to phasors and adjusted by the phase of maximum firing in nearby units (i.e., maximum UP). We tested whether, in the SOZ, HFOs and associated action potentials (APs) occur more often at the UP-DOWN transition. We also examined ripple temporal correlations using cross-correlograms. RESULTS: At the group level in the SOZ, HFO and HFO-associated AP probability was highest during the UP-DOWN transition of slow wave excitability (p < < .001). In the non-SOZ, HFO and HFO-associated AP was highest during the DOWN-UP transition (p < < .001). At the unit level in the SOZ, 15.6% and 20% of units exhibited more robust firing during ripples (Cohen's d = .11-.83) and fast ripples (d = .36-.90) at the UP-DOWN transition (p < .05 f.d.r. corrected), respectively. By comparison, also in the SOZ, 6.6% (d = .14-.30) and 8.5% (d = .33-.41) of units had significantly less firing during ripples and fast ripples at the UP-DOWN transition, respectively. Additional data shows that ripple and fast ripple temporal correlations, involving global slow waves, between the hippocampus, entorhinal cortex, and parahippocampal gyrus were reduced by >50% in the SOZ compared to the non-SOZ (N = 3). SIGNIFICANCE: The UP-DOWN transition of slow wave excitability facilitates the activation of pathological neurons to generate pHFOs. Ripple temporal correlations across brain regions may be important in memory consolidation and are disrupted in the SOZ, perhaps by pHFO generation.

Patterns of subregional cerebellar atrophy across epilepsy syndromes: An ENIGMA-Epilepsy study.

Objective: The intricate neuroanatomical structure of the cerebellum is of longstanding interest in epilepsy, but has been poorly characterized within the current cortico-centric models of this disease. We quantified cross-sectional regional cerebellar lobule volumes using structural MRI in 1,602 adults with epilepsy and 1,022 healthy controls across twenty-two sites from the global ENIGMA-Epilepsy working group. Methods: A state-of-the-art deep learning-based approach was employed that parcellates the cerebellum into 28 neuroanatomical subregions. Linear mixed models compared total and regional cerebellar volume in i) all epilepsies; ii) temporal lobe epilepsy with hippocampal sclerosis (TLE-HS); iii) non-lesional temporal lobe epilepsy (TLE-NL); iv) genetic generalised epilepsy; and (v) extra-temporal focal epilepsy (ETLE). Relationships were examined for cerebellar volume versus age at seizure onset, duration of epilepsy, phenytoin treatment, and cerebral cortical thickness. Results: Across all epilepsies, reduced total cerebellar volume was observed (d=0.42). Maximum volume loss was observed in the corpus medullare (dmax=0.49) and posterior lobe grey matter regions, including bilateral lobules VIIB (dmax= 0.47), Crus I/II (dmax= 0.39), VIIIA (dmax=0.45) and VIIIB (dmax=0.40). Earlier age at seizure onset (ηρ2max=0.05) and longer epilepsy duration (ηρ2max=0.06) correlated with reduced volume in these regions. Findings were most pronounced in TLE-HS and ETLE with distinct neuroanatomical profiles observed in the posterior lobe. Phenytoin treatment was associated with reduced posterior lobe volume. Cerebellum volume correlated with cerebral cortical thinning more strongly in the epilepsy cohort than in controls. Significance: We provide robust evidence of deep cerebellar and posterior lobe subregional grey matter volume loss in patients with chronic epilepsy. Volume loss was maximal for posterior subregions implicated in non-motor functions, relative to motor regions of both the anterior and posterior lobe. Associations between cerebral and cerebellar changes, and variability of neuroanatomical profiles across epilepsy syndromes argue for more precise incorporation of cerebellum subregions into neurobiological models of epilepsy.

High-rate leading spikes in propagating spike sequences predict seizure outcome in surgical patients with temporal lobe epilepsy.

Inter-ictal spikes aid in the diagnosis of epilepsy and in planning surgery of medication-resistant epilepsy. However, the localizing information from spikes can be unreliable because spikes can propagate, and the burden of spikes, often assessed as a rate, does not always correlate with the seizure onset zone or seizure outcome. Recent work indicates identifying where spikes regularly emerge and spread could localize the seizure network. Thus, the current study sought to better understand where and how rates of single and coupled spikes, and especially brain regions with high-rate and leading spike of a propagating sequence, informs the extent of the seizure network. In 37 patients with medication-resistant temporal lobe seizures, who had surgery to treat their seizure disorder, an algorithm detected spikes in the pre-surgical depth inter-ictal EEG. A separate algorithm detected spike propagation sequences and identified the location of leading and downstream spikes in each sequence. We analysed the rate and power of single spikes on each electrode and coupled spikes between pairs of electrodes, and the proportion of sites with high-rate, leading spikes in relation to the seizure onset zone of patients seizure free (n = 19) and those with continuing seizures (n = 18). We found increased rates of single spikes in mesial temporal seizure onset zone (ANOVA, P < 0.001, η2 = 0.138), and increased rates of coupled spikes within, but not between, mesial-, lateral- and extra-temporal seizure onset zone of patients with continuing seizures (P < 0.001; η2 = 0.195, 0.113 and 0.102, respectively). In these same patients, there was a higher proportion of brain regions with high-rate leaders, and each sequence contained a greater number of spikes that propagated with a higher efficiency over a longer distance outside the seizure onset zone than patients seizure free (Wilcoxon, P = 0.0172). The proportion of high-rate leaders in and outside the seizure onset zone could predict seizure outcome with area under curve = 0.699, but not rates of single or coupled spikes (0.514 and 0.566). Rates of coupled spikes to a greater extent than single spikes localize the seizure onset zone and provide evidence for inter-ictal functional segregation, which could be an adaptation to avert seizures. Spike rates, however, have little value in predicting seizure outcome. High-rate spike sites leading propagation could represent sources of spikes that are important components of an efficient seizure network beyond the clinical seizure onset zone, and like the seizure onset zone these, too, need to be removed, disconnected or stimulated to increase the likelihood for seizure control.

Fast ripples reflect increased excitability that primes epileptiform spikes.

The neuronal circuit disturbances that drive inter-ictal and ictal epileptiform discharges remain elusive. Using a combination of extra-operative macro-electrode and micro-electrode inter-ictal recordings in six pre-surgical patients during non-rapid eye movement sleep, we found that, exclusively in the seizure onset zone, fast ripples (200-600 Hz), but not ripples (80-200 Hz), frequently occur <300 ms before an inter-ictal intra-cranial EEG spike with a probability exceeding chance (bootstrapping, P < 1e-5). Such fast ripple events are associated with higher spectral power (P < 1e-10) and correlated with more vigorous neuronal firing than solitary fast ripple (generalized linear mixed-effects model, P < 1e-9). During the intra-cranial EEG spike that follows a fast ripple, action potential firing is lower than during an intra-cranial EEG spike alone (generalized linear mixed-effects model, P < 0.05), reflecting an inhibitory restraint of intra-cranial EEG spike initiation. In contrast, ripples do not appear to prime epileptiform spikes. We next investigated the clinical significance of pre-spike fast ripple in a separate cohort of 23 patients implanted with stereo EEG electrodes, who underwent resections. In non-rapid eye movement sleep recordings, sites containing a high proportion of fast ripple preceding intra-cranial EEG spikes correlate with brain areas where seizures begin more than solitary fast ripple (P < 1e-5). Despite this correlation, removal of these sites does not guarantee seizure freedom. These results are consistent with the hypothesis that fast ripple preceding EEG spikes reflect an increase in local excitability that primes EEG spike discharges preferentially in the seizure onset zone and that epileptogenic brain regions are necessary, but not sufficient, for initiating inter-ictal epileptiform discharges.

Pathological neurons generate ripples at the UP-DOWN transition disrupting information transfer.

Objective: To confirm and investigate why pathological HFOs (pHFOs), including Ripples [80-200 Hz] and fast ripples [200-600 Hz], are generated during the UP-DOWN transition of the slow wave and if pHFOs interfere with information transmission. Methods: We isolated 217 total units from 175.95 iEEG contact-hours of synchronized macro- and microelectrode recordings from 6 patients. Sleep slow oscillation (0.1-2 Hz) epochs were identified in the iEEG recording. iEEG HFOs that occurred superimposed on the slow wave were transformed to phasors and adjusted by the phase of maximum firing in nearby units (i.e., maximum UP). We tested whether, in the seizure onset zone (SOZ), HFOs and associated action potentials (AP) occur more often at the UP-DOWN transition. We also examined ripple temporal correlations using cross correlograms. Results: At the group level in the SOZ, HFO and HFO-associated AP probability was highest during the UP-DOWN transition of slow wave excitability (p<<0.001). In the non-SOZ, HFO and HFO-associated AP was highest during the DOWN-UP transition (p<<0.001). At the unit level in the SOZ, 15.6% and 20% of units exhibited more robust firing during ripples (Cohen's d=0.11-0.83) and fast ripples (d=0.36-0.90) at the UP-DOWN transition (p<0.05 f.d.r corrected), respectively. By comparison, also in the SOZ, 6.6% (d=0.14-0.30) and 8.5% (d=0.33-0.41) of units had significantly less firing during ripples and fast ripples at the UP-DOWN transition, respectively. Additional data shows ripple temporal correlations, involving global slow waves, between the hippocampus, entorhinal cortex, and parahippocampal gyrus were reduced by ~50-80% in the SOZ compared to the non-SOZ (N=3). Significance: The UP-DOWN transition of slow wave excitability facilitates the activation of pathological neurons to generate pHFOs. The pathological neurons and pHFOs disrupt ripple temporal correlations across brain regions that transfer information and may be important in memory consolidation.

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.

Image data harmonization tools for the analysis of post-traumatic epilepsy development in preclinical multisite MRI studies.

Preclinical MRI studies have been utilized for the discovery of biomarkers that predict post-traumatic epilepsy (PTE). However, these single site studies often lack statistical power due to limited and homogeneous datasets. Therefore, multisite studies, such as the Epilepsy Bioinformatics Study for Antiepileptogenic Therapy (EpiBioS4Rx), are developed to create large, heterogeneous datasets that can lead to more statistically significant results. EpiBioS4Rx collects preclinical data internationally across sites, including the United States, Finland, and Australia. However, in doing so, there are robust normalization and harmonization processes that are required to obtain statistically significant and generalizable results. This work describes the tools and procedures used to harmonize multisite, multimodal preclinical imaging data acquired by EpiBioS4Rx. There were four main harmonization processes that were utilized, including file format harmonization, naming convention harmonization, image coordinate system harmonization, and diffusion tensor imaging (DTI) metrics harmonization. By using Python tools and bash scripts, the file formats, file names, and image coordinate systems are harmonized across all the sites. To harmonize DTI metrics, values are estimated for each voxel in an image to generate a histogram representing the whole image. Then, the Quantitative Imaging Toolkit (QIT) modules are utilized to scale the mode to a value of one and depict the subsequent harmonized histogram. The standardization of file formats, naming conventions, coordinate systems, and DTI metrics are qualitatively assessed. The histograms of the DTI metrics were generated for all the individual rodents per site. For inter-site analysis, an average of the individual scans was calculated to create a histogram that represents each site. In order to ensure the analysis can be run at the level of individual animals, the sham and TBI cohort were analyzed separately, which depicted the same harmonization factor. The results demonstrate that these processes qualitatively standardize the file formats, naming conventions, coordinate systems, and DTI metrics of the data. This assists in the ability to share data across the study, as well as disseminate tools that can help other researchers to strengthen the statistical power of their studies and analyze data more cohesively.

Outcome of stereo-electroencephalography with single-unit recording in drug-refractory epilepsy.

OBJECTIVE: The aim of this study was to evaluate the utility and safety of "hybrid" stereo-electroencephalography (SEEG) in guiding epilepsy surgery and in providing information at single-neuron levels (i.e., single-unit recording) to further the understanding of the mechanisms of epilepsy and the neurocognitive processes unique to humans. METHODS: The authors evaluated 218 consecutive patients undergoing SEEG procedures from 1993 through 2018 at a single academic medical center to assess the utility and safety of this technique in both guiding epilepsy surgery and providing single-unit recordings. The hybrid electrodes used in this study contained macrocontacts and microwires to simultaneously record intracranial EEG and single-unit activity (hybrid SEEG). The outcomes of SEEG-guided surgical interventions were examined, as well as the yield and scientific utility of single-unit recordings in 213 patients who participated in the research involving single-unit recordings. RESULTS: All patients underwent SEEG implantation by a single surgeon and subsequent video-EEG monitoring (mean of 10.2 electrodes per patient and 12.0 monitored days). Epilepsy networks were localized in 191 (87.6%) patients. Two clinically significant procedural complications (one hemorrhage and one infection) were noted. Of 130 patients who underwent subsequent focal epilepsy surgery with a minimum 12-month follow-up, 102 (78.5%) underwent resective surgery and 28 (21.5%) underwent closed-loop responsive neurostimulation (RNS) with or without resection. Seizure freedom was achieved in 65 (63.7%) patients in the resective group. In the RNS group, 21 (75.0%) patients achieved 50% or greater seizure reduction. When the initial period of 1993 through 2013 before responsive neurostimulator implantation in 2014 was compared with the subsequent period of 2014 through 2018, the proportion of SEEG patients undergoing focal epilepsy surgery grew from 57.9% to 79.7% due to the advent of RNS, despite a decline in focal resective surgery from 55.3% to 35.6%. A total of 18,680 microwires were implanted in 213 patients, resulting in numerous significant scientific findings. Recent recordings from 35 patients showed a yield of 1813 neurons, with a mean yield of 51.8 neurons per patient. CONCLUSIONS: Hybrid SEEG enables safe and effective localization of epileptogenic zones to guide epilepsy surgery and provides unique scientific opportunities to investigate neurons from various brain regions in conscious patients. This technique will be increasingly utilized due to the advent of RNS and may prove a useful approach to probe neuronal networks in other brain disorders.

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.

Fast ripples reflect increased excitability that primes epileptiform spikes.

The neuronal circuit disturbances that drive interictal and ictal epileptiform discharges remains elusive. Using a combination of extraoperative macro- and micro-electrode interictal recordings in six presurgical patients during non-rapid eye movement (REM) sleep we found that, exclusively in the seizure onset zone, fast ripples (FR; 200-600Hz), but not ripples (80-200 Hz), frequently occur <300 msec before an interictal intracranial EEG (iEEG) spike with a probability exceeding chance (bootstrapping, p<1e-5). Such FR events are associated with higher spectral power (p<1e-10) and correlated with more vigorous neuronal firing than solitary FR (generalized linear mixed-effects model, GLMM, p<1e-3) irrespective of FR power. During the iEEG spike that follows a FR, action potential firing is lower than during a iEEG spike alone (GLMM, p<1e-10), reflecting an inhibitory restraint of iEEG spike initiation. In contrast, ripples do not appear to prime epileptiform spikes. We next investigated the clinical significance of pre-spike FR in a separate cohort of 23 patients implanted with stereo EEG electrodes who underwent resections. In non-REM sleep recordings, sites containing a high proportion of FR preceding iEEG spikes correlate with brain areas where seizures begin more than solitary FR (p<1e-5). Despite this correlation, removal of these sites does not guarantee seizure freedom. These results are consistent with the hypothesis that FR preceding EEG spikes reflect an increase in local excitability that primes EEG spike discharges preferentially in the seizure onset zone and that epileptogenic brain regions are necessary, but not sufficient, for initiating interictal epileptiform discharges.

Early preclinical plasma protein biomarkers of brain trauma are influenced by early seizures and levetiracetam.

OBJECTIVE: We used the lateral fluid percussion injury (LFPI) model of moderate-to-severe traumatic brain injury (TBI) to identify early plasma biomarkers predicting injury, early post-traumatic seizures or neuromotor functional recovery (neuroscores), considering the effect of levetiracetam, which is commonly given after severe TBI. METHODS: Adult male Sprague-Dawley rats underwent left parietal LFPI, received levetiracetam (200 mg/kg bolus, 200 mg/kg/day subcutaneously for 7 days [7d]) or vehicle post-LFPI, and were continuously video-EEG recorded (n = 14/group). Sham (craniotomy only, n = 6), and naïve controls (n = 10) were also used. Neuroscores and plasma collection were done at 2d or 7d post-LFPI or equivalent timepoints in sham/naïve. Plasma protein biomarker levels were determined by reverse phase protein microarray and classified according to injury severity (LFPI vs. sham/control), levetiracetam treatment, early seizures, and 2d-to-7d neuroscore recovery, using machine learning. RESULTS: Low 2d plasma levels of Thr231 -phosphorylated tau protein (pTAU-Thr231 ) and S100B combined (ROC AUC = 0.7790) predicted prior craniotomy surgery (diagnostic biomarker). Levetiracetam-treated LFPI rats were differentiated from vehicle treated by the 2d-HMGB1, 2d-pTAU-Thr231 , and 2d-UCHL1 plasma levels combined (ROC AUC = 0.9394) (pharmacodynamic biomarker). Levetiracetam prevented the seizure effects on two biomarkers that predicted early seizures only among vehicle-treated LFPI rats: pTAU-Thr231 (ROC AUC = 1) and UCHL1 (ROC AUC = 0.8333) (prognostic biomarker of early seizures among vehicle-treated LFPI rats). Levetiracetam-resistant early seizures were predicted by high 2d-IFNγ plasma levels (ROC AUC = 0.8750) (response biomarker). 2d-to-7d neuroscore recovery was best predicted by higher 2d-S100B, lower 2d-HMGB1, and 2d-to-7d increase in HMGB1 or decrease in TNF (P < 0.05) (prognostic biomarkers). SIGNIFICANCE: Antiseizure medications and early seizures need to be considered in the interpretation of early post-traumatic biomarkers.

Stimulation better targets fast-ripple generating networks in super responders to the responsive neurostimulator system.

How responsive neurostimulation (RNS) decreases seizure frequency is unclear. Stimulation may alter epileptic networks during inter-ictal epochs. Definitions of the epileptic network vary but fast ripples (FRs) may be an important substrate. We, therefore, examined whether stimulation of FR-generating networks differed in RNS super responders and intermediate responders. In 10 patients, with subsequent RNS placement, we detected FRs from stereo-electroencephalography (SEEG) contacts during pre-surgical evaluation. The normalized coordinates of the SEEG contacts were compared with those of the eight RNS contacts, and RNS-stimulated SEEG contacts were defined as those within 1.5 cm3 of the RNS contacts. We compared the post-RNS placement seizure outcome to (1) the ratio of stimulated SEEG contacts in the seizure-onset zone (SOZ stimulation ratio [SR]); (2) the ratio of FR events on stimulated contacts (FR SR); and (3) the global efficiency of the FR temporal correlational network on stimulated contacts (FR SGe). We found that the SOZ SR (p = .18) and FR SR (p = .06) did not differ in the RNS super responders and intermediate responders, but the FR SGe did (p = .02). In super responders, highly active desynchronous sites of the FR network were stimulated. RNS that better targets FR networks, as compared to the SOZ, may reduce epileptogenicity more.

Graph theoretical measures of fast ripple networks improve the accuracy of post-operative seizure outcome prediction.

Fast ripples (FR) are a biomarker of epileptogenic brain, but when larger portions of FR generating regions are resected seizure freedom is not always achieved. To evaluate and improve the diagnostic accuracy of FR resection for predicting seizure freedom we compared the FR resection ratio (RR) with FR network graph theoretical measures. In 23 patients FR were semi-automatically detected and quantified in stereo EEG recordings during sleep. MRI normalization and co-registration localized contacts and relation to resection margins. The number of FR, and graph theoretical measures, which were spatial (i.e., FR rate-distance radius) or temporal correlational (i.e., FR mutual information), were compared with the resection margins and with seizure outcome We found that the FR RR did not correlate with seizure-outcome (p > 0.05). In contrast, the FR rate-distance radius resected difference and the FR MI mean characteristic path length RR did correlate with seizure-outcome (p < 0.05). Retesting of positive FR RR patients using either FR rate-distance radius resected difference or the FR MI mean characteristic path length RR reduced seizure-free misclassifications from 44 to 22% and 17%, respectively. These results indicate that graph theoretical measures of FR networks can improve the diagnostic accuracy of the resection of FR events for predicting seizure freedom.