The Surgical Method of Craniectomy Differentially Affects Acute Seizures, Brain Deformation and Behavior in a TBI Animal Model.

Traumatic brain injury (TBI) is the leading cause of morbidity and mortality worldwide. Multiple injury models have been developed to study this neurological disorder. One such model is the lateral fluid-percussion injury (LFPI) rodent model. The LFPI model can be generated with different surgical procedures that could affect the injury and be reflected in neurobehavioral dysfunction and acute EEG changes. A craniectomy was performed either with a trephine hand drill or with a trephine electric drill that was centered over the left hemisphere of adult, male Sprague Dawley rats. Sham craniectomy groups were assessed by hand-drilled (ShamHMRI) and electric-drilled (ShamEMRI) to evaluate by MRI. Then, TBI was induced in separate groups (TBIH) and (TBIE) using a fluid-percussion device. Sham-injured rats (ShamH/ShamE) underwent the same surgical procedures as the TBI rats. During the same surgery session, rats were implanted with screw and microwire electrodes positioned in the neocortex and hippocampus and the EEG activity was recorded 24 hours for the first 7 days after TBI for assessing the acute EEG seizure and Gamma Event Coupling (GEC). The electric drilling craniectomy induced greater tissue damage and sensorimotor deficits compared to the hand drill. Analysis of the EEG revealed acute seizures in at least one animal from each group after the procedure. Both TBI and Sham rats from the electric drill groups had a significant greater total number of seizures than the animals that were craniectomized manually (p<0.05). Similarly, EEG functional connectivity was lower in ShamE compared to ShamH rats. These results suggest that electrical versus hand drilling craniectomies produce cortical injury in addition to the LFPI which increases the likelihood for acute post-traumatic seizures. Differences in the surgical approach could be one reason for the variability in the injury that makes it difficult to replicate results between preclinical TBI studies.

Discovering Neurophysiological Characteristics of Pathological High-Frequency Oscillations in Epilepsy with an Explainable Deep Generative Model.

Interictal high-frequency oscillation (HFO) is a promising biomarker of the epileptogenic zone (EZ). However, objective definitions to distinguish between pathological and physiological HFOs have remained elusive, impeding HFOs' clinical applications. We employed self-supervised deep generative variational autoencoders to learn such discriminative HFO features directly from their morphologies in a data-driven manner. We studied a large retrospective cohort of 185 patients who underwent intracranial monitoring and analyzed 686,410 candidate HFO events collected from 18,265 brain contacts across diverse brain regions. The model automatically clustered HFOs into distinct morphological groups in the latent space. One cluster consisted of putative morphologically defined pathological HFOs (mpHFOs): HFOs in that cluster were observed to be associated with spikes and exhibited high signal intensity both in the HFO band (>80 Hz) at detection and in the sub-HFO band (10-80 Hz) surrounding the detection and were primarily localized in the seizure onset zone (SOZ). Moreover, resection of brain regions based on a higher prevalence of interictal mpHFOs better predicted postoperative seizure outcomes than current clinical standards based on SOZ removal. Our self-supervised, explainable, deep generative model distills pathological HFOs and thus potentially helps delineate the EZ purely from interictal intracranial EEG data.

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.

Simulated resections and RNS placement can optimize post-operative 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 (SOZ) 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 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 (RR), 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 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 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 (RNS) and virtually targeted the RNS stimulation contacts proximal to sites generating FR at highest rates to determine if the simulated value of the stimulated SOZ and stimulated FR metrics would trend toward those patients with a better seizure outcome. Our results suggest: 1) FR measures can accurately predict whether a resection, defined by the standard of care, will result in seizure freedom; 2) utilizing FR alone for planning an efficacious surgery can be associated with larger resections; 3) 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 4) more work is required to determine if targeting RNS stimulation contact proximal to FR generating sites will improve seizure outcome.

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.