This study presents a systematic machine-learning approach for classifying acute pain from raw electrophysiological signals. We address binary and ternary classification tasks, leveraging Power-In-Band (PIB) and signal coherence as distinguishing features. Our method evaluates the effectiveness of traditional machine learning algorithms on a manually curated electrophysiological dataset obtained from intracranial electroencephalography (iEEG), offering valuable insights into model performance for pain detection. Furthermore, we identify critical electrode pairings associated with acute pain, providing a clearer understanding of the neural markers that differentiate pain states. This work highlights the potential of targeted feature engineering in advancing pain classification, setting the stage for future enhancements in real-time and personalized pain assessment tools. Additionally, these findings have promising applications in neuromodulation and Deep Brain Stimulation (DBS), where adaptive and closed-loop systems could leverage identified pain markers to modulate pain-related brain regions more precisely, offering improved therapeutic options for chronic pain management

BACKGROUND:In cases of medically refractory epilepsy that is multifocal or involving eloquent cortex, a combined resection and neuromodulatory approach may offer greater seizure reduction than resection alone. Here the authors describe indications, technical considerations, and outcomes of patients in whom resection and responsive neurostimulation (RNS) device implantation were performed in the same operation. OBSERVATIONS:Three patients met inclusion criteria for this study. Preoperative sEEG demonstrated multifocal or eloquent cortex-involving seizure onset zones in all patients. In 2 patients, one depth electrode in the contralateral hippocampus and a cortical strip electrode ipsilateral to the resection were implanted. One patient received bilateral hippocampal depth electrodes and resection. Technical considerations include prioritization of stereotactic accuracy, avoidance of pulse generator discharge, and distant implantation of the generator relative to the resection site. At 1 year, 2 patients achieved Engel class IIIA, and 1 patient achieved class IIB. LESSONS:The authors suggest consideration of this combined approach of simultaneous RNS device implantation and craniotomy for resection when intracranial data indicate a role for resection but it is unlikely to achieve seizure freedom alone. https://thejns.org/doi/10.3171/CASE24895.

Maturation of human brain structure has been well-studied, but developmental changes to brain physiology are not as well understood. One consistent finding is that the peak alpha rhythm frequency (PAF) increases throughout childhood. Another is that resting-state functional connectivity shifts from sensorimotor regions in children to association regions in adolescents, a reorganization along a hierarchy called the sensorimotor-to-association (S-A) axis. In mature brains, the S-A axis has been parcellated physiologically using the duration of persistent neural activity, known as the intrinsic neural timescale (INT), which increases along the hierarchy. Here, we studied the development of PAF and INT in a cohort of patients with epilepsy 3-33 yr of age undergoing intracranial electrocorticographic (ECoG) monitoring. Given the well-known developmental trajectory of PAF, and the ability to delineate hierarchy using INT, we hypothesized that changes to PAF and INT would correlate across development, but that their relationship may be influenced by hierarchy. Consequently, we predicted that age-dependent PAF increases would accompany INT decreases, and we tested whether their relationship varied between sensorimotor and association regions. We found that PAF increased significantly with age in both sensorimotor and association regions, whereas age-dependent INT decreases were only significant in association regions. Supporting this finding, we found a significant negative relationship between PAF and INT in association regions, but not sensorimotor regions. Together, our results provide further evidence that developmental divisions across the S-A axis manifest in the relationships between neurophysiological measures.NEW & NOTEWORTHY We report a novel developmental relationship between the human resting-state alpha rhythm frequency and the duration of intrinsic neural timescales. Using resting-state electrocorticography, we found that alpha frequency increased with age at either end of the sensorimotor-to-association cortical hierarchy, whereas intrinsic neural timescales only decreased with age in association regions. This negative correlation between alpha frequency and intrinsic timescale was only evident in association regions, further linking functional maturation and cortical hierarchy.

With the development and characterization of biomarkers that may reflect neural network state as well as a patient’s clinical deficits, there is growing interest in more complex stimulation designs. While current implantable neuromodulation systems offer pathways to expand the design and application of adaptive stimulation paradigms, technological drawbacks of these systems limit adaptive neuromodulation exploration. In this paper, we discuss the implementation of a phase-triggered stimulation paradigm using a research platform composed of an investigational system known as the CorTec Brain Interchange (CorTec GmbH, Freiburg, Germany), and an open-source software tool known as OMNI-BIC. We then evaluate the stimulation paradigm’s performance in both benchtop and in vivo human demonstrations. Our findings indicate that the Brain Interchange and OMNI-BIC platform is capable of reliable administration of phase-triggered stimulation and has the potential to help expand investigation within the adaptive neuromodulation design space.

Background Electroencephalography (EEG) and electrocorticography (ECoG) recordings have been used to decode finger movements by analyzing brain activity. Traditional methods focused on single bandpass power changes for movement decoding, utilizing machine learning models requiring manual feature extraction. New Method This study introduces a 3D convolutional neural network (3D-CNN) model to decode finger movements using ECoG data. The model employs adaptive, explainable AI (xAI) techniques to interpret the physiological relevance of brain signals. ECoG signals from epilepsy patients during awake craniotomy were processed to extract power spectral density across multiple frequency bands. These data formed a 3D matrix used to train the 3D-CNN to predict finger trajectories. Results The 3D-CNN model showed significant accuracy in predicting finger movements, with root-mean-square error (RMSE) values of 0.26-0.38 for single finger movements and 0.20-0.24 for combined movements. Explainable AI techniques, Grad-CAM and SHAP, identified the high gamma (HG) band as crucial for movement prediction, showing specific cortical regions involved in different finger movements. These findings highlighted the physiological significance of the HG band in motor control. Comparison with Existing Methods The 3D-CNN model outperformed traditional machine learning approaches by effectively capturing spatial and temporal patterns in ECoG data. The use of xAI techniques provided clearer insights into the model's decision-making process, unlike the "black box" nature of standard deep learning models. Conclusions The proposed 3D-CNN model, combined with xAI methods, enhances the decoding accuracy of finger movements from ECoG data. This approach offers a more efficient and interpretable solution for brain-computer interface (BCI) applications, emphasizing the HG band's role in motor control.

[This corrects the article DOI: 10.3389/fnins.2021.769872.].

Response inhibition refers to the ability to suppress inappropriate actions that interfere with goal-driven behavior. The inferior frontal gyrus (IFG) is known to be associated with inhibition of a motor response by assuming executive control over motor cortex outputs. This study aimed to evaluate the pediatric development of response inhibition through subdural electrocorticography (ECoG) recording. Subdural ECoG recorded neural activities simultaneously during a Go/No-Go task, which was optimized for children. Different frequency power [theta: 4-8 Hz; beta: 12-40 Hz; high-gamma (HG): 70-200 Hz] was estimated within the IFG and motor cortex. Age-related analysis was computed by each bandpass power ratio between Go and No-Go conditions, and phase-amplitude coupling (PAC) over IFG by using the modulating index metric in two conditions. For all the eight pediatric patients, HG power was more activated in No-Go trials than in Go trials, in either right- or left-side IFG when available. In the IFG region, the power over theta and HG in No-Go conditions was higher than those in Go conditions, with significance over the right side (p < 0.05). The age-related lateralization from both sides to the right side was observed from the ratio of HG power and PAC value between the No-Go and Go trials. In the pediatric population, the role of motor inhibition was observed in both IFG, with age-related lateralization to the right side, which was proved in the previous functional magnetic resonance imaging studies. In this study, the evidence correlation of age and response inhibition was observed directly by the evidence of cortical recordings.

Electrical stimulation is a powerful tool for targeted neurorehabilitation, and recent work in adaptive stimulation where stimulation can be adjusted in real-time has shown promise in improving stimulation outcomes and reducing stimulation-induced side effects. Mapping the relationship between electrical stimulation input and neural activity response can help reveal their interactions and can give us tools to iterate and improve on our stimulation protocols. Here, we introduce methods for identifying low-dimensional subspaces of human intracortical responses to electrical stimulation in invasive electroencephalography. In epilepsy patients (n=4) undergoing clinical monitoring, we applied a stimulation protocol of varying amplitude and frequency in 5-second intervals to capture a range of responses to different stimulation conditions. We characterized these responses using time-frequency spectral power, applied baseline subtraction and outlier removal procedures, and performed principal component analysis across frequencies. We identified that intracortical responses to different stimulation conditions can be represented in a 3-dimensional subspace, accounting for more than 95% of the variance. Using support vector machine classification, we demonstrated separability of intracortical responses in different stimulation conditions across subjects, where this separability was contingent on applying baseline subtraction and outlier removal. Our results represent a first step towards building a predictive model of neural response from stimulation input, an important prerequisite for adaptive closed-loop stimulation for targeted neurorehabilitation.

Accurate anatomical localization of intracranial electrodes is important for identifying the seizure foci in patients with epilepsy and for interpreting effects from cognitive studies employing intracranial electroencephalography. Localization is typically performed by coregistering postimplant computed tomography (CT) with preoperative magnetic resonance imaging (MRI). Electrodes are then detected in the CT, and the corresponding brain region is identified using the MRI. Many existing software packages for electrode localization chain together separate preexisting programs or rely on command line instructions to perform the various localization steps, making them difficult to install and operate for a typical user. Further, many packages provide solutions for some, but not all, of the steps needed for confident localization. We have developed software, Locate electrodes Graphical User Interface (LeGUI), that consists of a single interface to perform all steps needed to localize both surface and depth/penetrating intracranial electrodes, including coregistration of the CT to MRI, normalization of the MRI to the Montreal Neurological Institute template, automated electrode detection for multiple types of electrodes, electrode spacing correction and projection to the brain surface, electrode labeling, and anatomical targeting. The software is written in MATLAB, core image processing is performed using the Statistical Parametric Mapping toolbox, and standalone executable binaries are available for Windows, Mac, and Linux platforms. LeGUI was tested and validated on 51 datasets from two universities. The total user and computational time required to process a single dataset was approximately 1 h. Automatic electrode detection correctly identified 4362 of 4695 surface and depth electrodes with only 71 false positives. Anatomical targeting was verified by comparing electrode locations from LeGUI to locations that were assigned by an experienced neuroanatomist. LeGUI showed a 94% match with the 482 neuroanatomist-assigned locations. LeGUI combines all the features needed for fast and accurate anatomical localization of intracranial electrodes into a single interface, making it a valuable tool for intracranial electrophysiology research.

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is a clinically effective tool for treating medically refractory Parkinson's disease (PD), but its neural mechanisms remain debated. Previous work has demonstrated that STN DBS results in evoked potentials (EPs) in the primary motor cortex (M1), suggesting that modulation of cortical physiology may be involved in its therapeutic effects. Due to technical challenges presented by high-amplitude DBS artifacts, these EPs are often measured in response to low-frequency stimulation, which is generally ineffective at PD symptom management. This study aims to characterize STN-to-cortex EPs seen during clinically relevant high-frequency STN DBS for PD. Intraoperatively, we applied STN DBS to 6 PD patients while recording electrocorticography (ECoG) from an electrode strip over the ipsilateral central sulcus. Using recently published techniques, we removed large stimulation artifacts to enable quantification of STN-to-cortex EPs. Two cortical EPs were observed - one synchronized with DBS onset and persisting during ongoing stimulation, and one immediately following DBS offset, here termed the "start" and the "end" EPs respectively. The start EP is, to our knowledge, the first long-latency cortical EP reported during ongoing high-frequency DBS. The start and end EPs differ in magnitude (p < 0.05) and latency (p < 0.001), and the end, but not the start, EP magnitude has a significant relationship (p < 0.001, adjusted for random effects of subject) to ongoing high gamma (80-150 Hz) power during the EP. These contrasts may suggest mechanistic or circuit differences in EP production during the two time periods. This represents a potential framework for relating DBS clinical efficacy to the effects of a variety of stimulation parameters on EPs.

OBJECTIVE:Laser interstitial thermal therapy (LITT) is a minimally invasive surgery for mesial temporal lobe epilepsy (mTLE), but the effects of individual patient anatomy and location of ablation volumes affect seizure outcomes. The purpose of this study is to see if features of individual patient structural connectomes predict surgical outcomes after LITT for mTLE. METHODS:This is a retrospective analysis of seizure outcomes of LITT for mTLE in 24 patients. We use preoperative diffusion tensor imaging (DTI) to simulate changes in structural connectivity after laser ablation. A two-step machine-learning algorithm is applied to predict seizure outcomes from the change in connectomic features after surgery. RESULTS:Although node-based network features such as clustering coefficient and betweenness centrality have some predictive value, changes in connection strength between mesial temporal regions predict seizure outcomes significantly better. Changes in connection strength between the entorhinal cortex (EC), and the insula, hippocampus, and amygdala, as well as between the temporal pole and hippocampus, predict Engel Class I outcomes with an accuracy of 88%. Analysis of the ablation location, as well as simulated, alternative ablations, reveals that a more medial, anterior, and inferior ablation volume is associated with a greater effect on these connections, and potentially on seizure outcomes. SIGNIFICANCE:Our results indicate (1) that seizure outcomes can be retrospectively predicted with excellent accuracy using changes in structural connectivity, and (2) that favorable connectomic changes are associated with an ablation volume involving relatively mesial, anterior, and inferior locations. These results may provide a framework whereby individual pre-operative structural connectomes can be used to optimize ablation volumes and improve outcomes in LITT for mTLE.

Studying the neural correlates of sleep can lead to revelations in our understanding of sleep and its interplay with different neurological disorders. Sleep research relies on manual annotation of sleep stages based on rules developed for healthy adults. Automating sleep stage annotation can expedite sleep research and enable us to better understand atypical sleep patterns. Our goal was to create a fully unsupervised approach to label sleep and wake states in human electro-corticography (ECoG) data from epilepsy patients. Here, we demonstrate that with continuous data from a single ECoG electrode, hidden semi-Markov models (HSMM) perform best in classifying sleep/wake states without excessive transitions, with a mean accuracy (n=4) of 85.2% compared to using K-means clustering (72.2%) and hidden Markov models (81.5%). Our results confirm that HSMMs produce meaningful labels for ECoG data and establish the groundwork to apply this model to cluster sleep stages and potentially other behavioral states.

ABSTRACTBackgroundConverging literatures suggest that deep brain stimulation (DBS) in Parkinson's disease affects multiple circuit mechanisms. One proposed mechanism is the normalization of primary motor cortex (M1) pathophysiology via effects on the hyperdirect pathway.ObjectivesWe hypothesized that DBS would reduce the current intensity necessary to modulate motor‐evoked potentials from focally applied direct cortical stimulation (DCS).MethodsIntraoperative subthalamic DBS, DCS, and preoperative diffusion tensor imaging data were acquired in 8 patients with Parkinson's disease.ResultsIn 7 of 8 patients, DBS significantly reduced the M1 DCS current intensity required to elicit motor‐evoked potentials. This neuromodulation was specific to select DBS bipolar configurations. In addition, the volume of activated tissue models of these configurations were significantly associated with overlap of the hyperdirect pathway.ConclusionsDBS reduces the current necessary to elicit a motor‐evoked potential using DCS. This supports a circuit mechanism of DBS effectiveness, potentially involving the hyperdirect pathway that speculatively may underlie reductions in hypokinetic abnormalities in Parkinson's disease. © 2020 International Parkinson and Movement Disorder Society

OBJECTIVE:Electrocorticography is an indispensable tool in identifying the epileptogenic zone in the presurgical evaluation of many epilepsy patients. Traditional electrocorticographic features (spikes, ictal onset changes, and recently high-frequency oscillations [HFOs]) rely on the presence of transient features that occur within or near epileptogenic cortex. Here the authors report on a novel corticography feature of epileptogenic cortex-covariation of high-gamma and beta frequency band power profiles. Band-limited power was measured from each recording site based on native physiological signal differences without relying on clinical ictal or interictal epileptogenic features. In this preliminary analysis, frequency windowed power correlation appears to be a specific marker of the epileptogenic zone. The authors' overall aim was to validate this observation with the location of the eventual resection and outcome. METHODS:The authors conducted a retrospective analysis of 13 adult patients who had undergone electrocorticography for surgical planning at their center. They quantified the correlation of high-gamma (70-200 Hz) and beta (12-18 Hz) band frequency power per electrode site during a cognitive task. They used a sliding window method to correlate the power of smoothed, Hilbert-transformed high-gamma and beta bands. They then compared positive and negative correlations between power in the high-gamma and beta bands in the setting of a hand versus a tongue motor task as well as within the resting state. Significant positive correlations were compared to surgically resected areas and outcomes based on reviewed records. RESULTS:Positive high-gamma and beta correlations appeared to predict the area of eventual resection and, preliminarily, surgical outcome independent of spike detection. In general, patients with the best outcomes had well-localized positive correlations (high-gamma and beta activities) to areas of eventual resection, while those with poorer outcomes displayed more diffuse patterns. CONCLUSIONS:Data in this study suggest that positive high-gamma and beta correlations independent of any behavioral metric may have clinical applicability in surgical decision-making. Further studies are needed to evaluate the clinical potential of this methodology. Additional work is also needed to relate these results to other methods, such as HFO detection or connectivity with other cortical areas.

Over the past decade, numerous neuroimaging studies based on hemodynamic markers of brain activity have examined the feeling of body ownership using perceptual body-illusions in humans. However, the direct electrophysiological correlates of body ownership at the cortical level remain unexplored. To address this, we studied the rubber hand illusion in 5 patients (3 males and 2 females) implanted with intracranial electrodes measuring cortical surface potentials. Increased high-γ (70-200 Hz) activity, an index of neuronal firing rate, in premotor and intraparietal cortices reflected the feeling of ownership. In both areas, high-γ increases were intimately coupled with the subjective illusion onset and sustained both during and in-between touches. However, intraparietal activity was modulated by tactile stimulation to a higher degree than the premotor cortex through effective connectivity with the hand-somatosensory cortex, which suggests different functional roles. These findings constitute the first intracranial electrophysiological characterization of the rubber hand illusion and extend our understanding of the dynamic mechanisms of body ownership.

Direct cortical stimulation (DCS) of primary somatosensory cortex (S1) could help restore sensation and provide task-relevant feedback in a neuroprosthesis. However, the psychophysics of S1 DCS is poorly studied, including any comparison to cutaneous haptic stimulation. We compare the response times to DCS of human hand somatosensory cortex through electrocorticographic grids with response times to haptic stimuli delivered to the hand in four subjects. We found that subjects respond significantly slower to S1 DCS than to natural, haptic stimuli for a range of DCS train durations. Median response times for haptic stimulation varied from 198 ms to 313 ms, while median responses to reliably perceived DCS ranged from 254 ms for one subject, all the way to 528 ms for another. We discern no significant impact of learning or habituation through the analysis of blocked trials, and find no significant impact of cortical stimulation train duration on response times. Our results provide a realistic set of expectations for latencies with somatosensory DCS feedback for future neuroprosthetic applications and motivate the study of neural mechanisms underlying human perception of somatosensation via DCS.

Brain–computer interfaces (BCIs) benefit greatly from performance feedback, but current systems lack automatic, task-independent feedback. Cortical responses elicited from user error have the potential to serve as state-based feedback to BCI decoders. To gain a better understanding of local error potentials, we investigate responsive cortical power underlying error-related potentials (ErrPs) from the human cortex during a one-dimensional center-out BCI task, tracking the topography of high-gamma (70–100 Hz) band power (HBP) specific to BCI error. We measured electrocorticography (ECoG) in three human subjects during dynamic, continuous control over BCI cursor velocity. Subjects used motor imagery and rest to move the cursor toward and subsequently dwell within a target region. We then identified and labeled epochs where the BCI decoder incorrectly moved the cursor in the direction opposite of the subject’s expectations (i.e., BCI error). We found increased HBP in various cortical areas 100–500 ms following BCI error with respect to epochs of correct, intended control. Significant responses were noted in primary somatosensory, motor, premotor, and parietal areas and generally regardless of whether the subject was using motor imagery or rest to move the cursor toward the target. Parts of somatosensory, temporal, and parietal areas exclusively had increased HBP when subjects were using motor imagery. In contrast, only part of the parietal cortex near the angular gyrus exclusively had an increase in HBP during rest. This investigation is, to our knowledge, the first to explore cortical fields changes in the context of continuous control in ECoG BCI. We present topographical changes in HBP characteristic specific to the generation of error. By focusing on continuous control, instead of on discrete control for simple selection, we investigate a more naturalistic setting and provide high ecological validity for characterizing error potentials. Such potentials could be considered as design elements for co-adaptive BCIs in the future as task-independent feedback to the decoder, allowing for more robust and individualized BCIs.

OBJECTIVE:The activation of the sensorimotor cortex as measured by electrocorticographic (ECoG) signals has been correlated with contralateral hand movements in humans, as precisely as the level of individual digits. However, the relationship between individual and multiple synergistic finger movements and the neural signal as detected by ECoG has not been fully explored. The authors used intraoperative high-resolution micro-ECoG (µECoG) on the sensorimotor cortex to link neural signals to finger movements across several context-specific motor tasks. METHODS:Three neurosurgical patients with cortical lesions over eloquent regions participated. During awake craniotomy, a sensorimotor cortex area of hand movement was localized by high-frequency responses measured by an 8 × 8 µECoG grid of 3-mm interelectrode spacing. Patients performed a flexion movement of the thumb or index finger, or a pinch movement of both, based on a visual cue. High-gamma (HG; 70-230 Hz) filtered µECoG was used to identify dominant electrodes associated with thumb and index movement. Hand movements were recorded by a dataglove simultaneously with µECoG recording. RESULTS:In all 3 patients, the electrodes controlling thumb and index finger movements were identifiable approximately 3-6-mm apart by the HG-filtered µECoG signal. For HG power of cortical activation measured with µECoG, the thumb and index signals in the pinch movement were similar to those observed during thumb-only and index-only movement, respectively (all p > 0.05). Index finger movements, measured by the dataglove joint angles, were similar in both the index-only and pinch movements (p > 0.05). However, despite similar activation across the conditions, markedly decreased thumb movement was observed in pinch relative to independent thumb-only movement (all p < 0.05). CONCLUSIONS:HG-filtered µECoG signals effectively identify dominant regions associated with thumb and index finger movement. For pinch, the µECoG signal comprises a combination of the signals from individual thumb and index movements. However, while the relationship between the index finger joint angle and HG-filtered signal remains consistent between conditions, there is not a fixed relationship for thumb movement. Although the HG-filtered µECoG signal is similar in both thumb-only and pinch conditions, the actual thumb movement is markedly smaller in the pinch condition than in the thumb-only condition. This implies a nonlinear relationship between the cortical signal and the motor output for some, but importantly not all, movement types. This analysis provides insight into the tuning of the motor cortex toward specific types of motor behaviors.

Prior studies using functional magnetic resonance imaging, electroencephalography, and magnetoencephalography have observed both structured patterns in resting-state functional connectivity and spontaneous longitudinal variation in connectivity patterns independent of a task. In this first study using electrocorticography (ECoG), we characterized spontaneous, intersession variation in resting-state functional connectivity not linked to a task. We evaluated pairwise connectivity between electrodes using three measures (phase locking value [PLV], amplitude correlation, and coherence) for six canonical frequency bands, capturing different characteristics of time-evolving signals. We grouped electrodes into 10 functional regions and used intraclass correlation (ICC) to estimate pairwise longitudinal stability. We found that stronger PLV (PLV ≥0.4) in theta through gamma bands and strong correlation in all bands (R2’s ≥0.6) are linked to substantial stability (ICC ≥0.6), but that stability does not imply strong phase locking or amplitude correlation. There was no notable link between strong coherence and high ICC. All within-region PLVs are markedly stable across frequencies. In addition, we highlight interaction patterns across several regions: parahippocampal/entorhinal cortex is characterized by stable, weak functional connectivity except self-connections. Dorsolateral prefrontal cortex connectivity is weak and unstable, except self-connections. Inferior parietal lobule has little stability despite narrow connectivity bounds. We confirm prior studies linking functional connectivity strength and intersession variability, extending into higher frequencies than other modalities, with greater spatial specificity than scalp electrophysiology. We suggest further studies quantitatively compare ECoG to other modalities and/or use these findings as a baseline to capture functional connectivity and dynamics linked to perturbations with a task or disease state.