Stimulation artifact correction method for estimation of early cortico-cortical evoked potentials

Single-pulse electrical stimulation artifact removal using the novel matching pursuit-based artifact reconstruction and removal method (MPARRM)

Objective.Single-pulse electrical stimulation (SPES) has been widely used to probe effective connectivity. However, analysis of the neural response is often confounded by stimulation artifacts. We developed a novel matching pursuit-based artifact reconstruction and removal method (MPARRM) capable of removing artifacts from stimulation-artifact-affected electrophysiological signals.Approach.To validate MPARRM across a wide range of potential stimulation artifact types, we performed a bench-top experiment in which we suspended electrodes in a saline solution to generate 110 types of real-world stimulation artifacts. We then added the generated stimulation artifacts to ground truth signals (stereoelectroencephalography signals from nine human subjects recorded during a receptive speech task), applied MPARRM to the combined signal, and compared the resultant denoised signal with the ground truth signal. We further applied MPARRM to artifact-affected neural signals recorded from the hippocampus while performing SPES on the ipsilateral basolateral amygdala in nine human subjects.Main results.MPARRM could remove stimulation artifacts without introducing spectral leakage or temporal spread. It accommodated variable stimulation parameters and recovered the early response to SPES within a wide range of frequency bands. Specifically, in the early response period (5-10 ms following stimulation onset), we found that the broadband gamma power (70-170 Hz) of the denoised signal was highly correlated with the ground truth signal (R=0.98±0.02, Pearson), and the broadband gamma activity of the denoised signal faithfully revealed the responses to the auditory stimuli within the ground truth signal with94%±1.47%sensitivity and99%±1.01%specificity. We further found that MPARRM could reveal the expected temporal progression of broadband gamma activity along the anterior-posterior axis of the hippocampus in response to the ipsilateral amygdala stimulation.Significance.MPARRM could faithfully remove SPES artifacts without confounding the electrophysiological signal components, especially during the early-response period. This method can facilitate the understanding of the neural response mechanisms of SPES.

Signatures of Electrical Stimulation Driven Network Interactions in the Human Limbic System

Stimulation-evoked signals are starting to be used as biomarkers to indicate the state and health of brain networks. The human limbic network, often targeted for brain stimulation therapy, is involved in emotion and memory processing. Previous anatomic, neurophysiological, and functional studies suggest distinct subsystems within the limbic network (Rolls, 2015). Studies using intracranial electrical stimulation, however, have emphasized the similarities of the evoked waveforms across the limbic network. We test whether these subsystems have distinct stimulation-driven signatures. In eight patients (four male, four female) with drug-resistant epilepsy, we stimulated the limbic system with single-pulse electrical stimulation. Reliable corticocortical evoked potentials (CCEPs) were measured between hippocampus and the posterior cingulate cortex (PCC) and between the amygdala and the anterior cingulate cortex (ACC). However, the CCEP waveform in the PCC after hippocampal stimulation showed a unique and reliable morphology, which we term the "limbic Hippocampus-Anterior nucleus of the thalamus-Posterior cingulate, HAP-wave." This limbic HAP-wave was visually distinct and separately decoded from the CCEP waveform in ACC after amygdala stimulation. Diffusion MRI data show that the measured end points in the PCC overlap with the end points of the parolfactory cingulum bundle rather than the parahippocampal cingulum, suggesting that the limbic HAP-wave may travel through fornix, mammillary bodies, and the anterior nucleus of the thalamus (ANT). This was further confirmed by stimulating the ANT, which evoked the same limbic HAP-wave but with an earlier latency. Limbic subsystems have unique stimulation-evoked signatures that may be used in the future to help network pathology diagnosis.SIGNIFICANCE STATEMENT The limbic system is often compromised in diverse clinical conditions, such as epilepsy or Alzheimer's disease, and characterizing its typical circuit responses may provide diagnostic insight. Stimulation-evoked waveforms have been used in the motor system to diagnose circuit pathology. We translate this framework to limbic subsystems using human intracranial stereo EEG (sEEG) recordings that measure deeper brain areas. Our sEEG recordings describe a stimulation-evoked waveform characteristic to the memory and spatial subsystem of the limbic network that we term the "limbic HAP-wave." The limbic HAP-wave follows anatomic white matter pathways from hippocampus to thalamus to the posterior cingulum and shows promise as a distinct biomarker of signaling in the human brain memory and spatial limbic network.

Quantitative approaches to guide epilepsy surgery from intracranial EEG

Over the past 10 years, the drive to improve outcomes from epilepsy surgery has stimulated widespread interest in methods to quantitatively guide epilepsy surgery from intracranial EEG (iEEG). Many patients fail to achieve seizure freedom, in part due to the challenges in subjective iEEG interpretation. To address this clinical need, quantitative iEEG analytics have been developed using a variety of approaches, spanning studies of seizures, interictal periods, and their transitions, and encompass a range of techniques including electrographic signal analysis, dynamical systems modeling, machine learning and graph theory. Unfortunately, many methods fail to generalize to new data and are sensitive to differences in pathology and electrode placement. Here, we critically review selected literature on computational methods of identifying the epileptogenic zone from iEEG. We highlight shared methodological challenges common to many studies in this field and propose ways that they can be addressed. One fundamental common pitfall is a lack of open-source, high-quality data, which we specifically address by sharing a centralized high-quality, well-annotated, multicentre dataset consisting of >100 patients to support larger and more rigorous studies. Ultimately, we provide a road map to help these tools reach clinical trials and hope to improve the lives of future patients.

Quantifying trial-by-trial variability during cortico-cortical evoked potential mapping of epileptogenic tissue

OBJECTIVE

Measuring cortico-cortical evoked potentials (CCEPs) is a promising tool for mapping epileptic networks, but it is not known how variability in brain state and stimulation technique might impact the use of CCEPs for epilepsy localization. We test the hypotheses that (1) CCEPs demonstrate systematic variability across trials and (2) CCEP amplitudes depend on the timing of stimulation with respect to endogenous, low-frequency oscillations.

METHODS

We studied 11 patients who underwent CCEP mapping after stereo-electroencephalography electrode implantation for surgical evaluation of drug-resistant epilepsy. Evoked potentials were measured from all electrodes after each pulse of a 30 s, 1 Hz bipolar stimulation train. We quantified monotonic trends, phase dependence, and standard deviation (SD) of N1 (15-50 ms post-stimulation) and N2 (50-300 ms post-stimulation) amplitudes across the 30 stimulation trials for each patient. We used linear regression to quantify the relationship between measures of CCEP variability and the clinical seizure-onset zone (SOZ) or interictal spike rates.

RESULTS

We found that N1 and N2 waveforms exhibited both positive and negative monotonic trends in amplitude across trials. SOZ electrodes and electrodes with high interictal spike rates had lower N1 and N2 amplitudes with higher SD across trials. Monotonic trends of N1 and N2 amplitude were more positive when stimulating from an area with higher interictal spike rate. We also found intermittent synchronization of trial-level N1 amplitude with low-frequency phase in the hippocampus, which did not localize the SOZ.

SIGNIFICANCE

These findings suggest that standard approaches for CCEP mapping, which involve computing a trial-averaged response over a .2-1 Hz stimulation train, may be masking inter-trial variability that localizes to epileptogenic tissue. We also found that CCEP N1 amplitudes synchronize with ongoing low-frequency oscillations in the hippocampus. Further targeted experiments are needed to determine whether phase-locked stimulation could have a role in localizing epileptogenic tissue.

Developmental trajectory of transmission speed in the human brain

The structure of the human connectome develops from childhood throughout adolescence to middle age, but how these structural changes affect the speed of neuronal signaling is not well described. In 74 subjects, we measured the latency of cortico-cortical evoked responses across association and U-fibers and calculated their corresponding transmission speeds. Decreases in conduction delays until at least 30 years show that the speed of neuronal communication develops well into adulthood.

Communication dynamics in the human connectome shape the cortex-wide propagation of direct electrical stimulation

Communication between gray matter regions underpins all facets of brain function. We study inter-areal communication in the human brain using intracranial EEG recordings, acquired following 29,055 single-pulse direct electrical stimulations in a total of 550 individuals across 20 medical centers (average of 87 ± 37 electrode contacts per subject). We found that network communication models-computed on structural connectivity inferred from diffusion MRI-can explain the causal propagation of focal stimuli, measured at millisecond timescales. Building on this finding, we show that a parsimonious statistical model comprising structural, functional, and spatial factors can accurately and robustly predict cortex-wide effects of brain stimulation (R²=46% in data from held-out medical centers). Our work contributes toward the biological validation of concepts in network neuroscience and provides insight into how connectome topology shapes polysynaptic inter-areal signaling. We anticipate that our findings will have implications for research on neural communication and the design of brain stimulation paradigms.

Variability of Single Pulse Electrical Stimulation Responses Recorded with Intracranial Electroencephalography in Epileptic Patients

Cohort studies of brain stimulations performed with stereo-electroencephalographic (SEEG) electrodes in epileptic patients allow to derive large scale functional connectivity. It is known, however, that brain responses to electrical or magnetic stimulation techniques are not always reproducible. Here, we study variability of responses to single pulse SEEG electrical stimulation. We introduce a second-order probability analysis, i.e. we extend estimation of connection probabilities, defined as the proportion of responses trespassing a statistical threshold (determined in terms of Z-score with respect to spontaneous neuronal activity before stimulation) over all responses and derived from a number of individual measurements, to an analysis of pairs of measurements.Data from 445 patients were processed. We found that variability between two equivalent measurements is substantial in particular conditions. For long ( > ~ 90 mm) distances between stimulating and recording sites, and threshold value Z = 3, correlation between measurements drops almost to zero. In general, it remains below 0.5 when the threshold is smaller than Z = 4 or the stimulating current intensity is 1 mA. It grows with an increase of either of these factors. Variability is independent of interictal spiking rates in the stimulating and recording sites.We conclude that responses to SEEG stimulation in the human brain are variable, i.e. in a subject at rest, two stimulation trains performed at the same electrode contacts and with the same protocol can give discrepant results. Our findings highlight an advantage of probabilistic interpretation of such results even in the context of a single individual.

Advances in human intracranial electroencephalography research, guidelines and good practices

Since the second-half of the twentieth century, intracranial electroencephalography (iEEG), including both electrocorticography (ECoG) and stereo-electroencephalography (sEEG), has provided an intimate view into the human brain. At the interface between fundamental research and the clinic, iEEG provides both high temporal resolution and high spatial specificity but comes with constraints, such as the individual's tailored sparsity of electrode sampling. Over the years, researchers in neuroscience developed their practices to make the most of the iEEG approach. Here we offer a critical review of iEEG research practices in a didactic framework for newcomers, as well addressing issues encountered by proficient researchers. The scope is threefold: (i) review common practices in iEEG research, (ii) suggest potential guidelines for working with iEEG data and answer frequently asked questions based on the most widespread practices, and (iii) based on current neurophysiological knowledge and methodologies, pave the way to good practice standards in iEEG research. The organization of this paper follows the steps of iEEG data processing. The first section contextualizes iEEG data collection. The second section focuses on localization of intracranial electrodes. The third section highlights the main pre-processing steps. The fourth section presents iEEG signal analysis methods. The fifth section discusses statistical approaches. The sixth section draws some unique perspectives on iEEG research. Finally, to ensure a consistent nomenclature throughout the manuscript and to align with other guidelines, e.g., Brain Imaging Data Structure (BIDS) and the OHBM Committee on Best Practices in Data Analysis and Sharing (COBIDAS), we provide a glossary to disambiguate terms related to iEEG research.

Simultaneous stereo-EEG and high-density scalp EEG recordings to study the effects of intracerebral stimulation parameters

BACKGROUND

Cortico-cortical evoked potentials (CCEPs) recorded by stereo-electroencephalography (SEEG) are a valuable tool to investigate brain reactivity and effective connectivity. However, invasive recordings are spatially sparse since they depend on clinical needs. This sparsity hampers systematic comparisons across-subjects, the detection of the whole-brain effects of intracortical stimulation, as well as their relationships to the EEG responses evoked by non-invasive stimuli.

OBJECTIVE

To demonstrate that CCEPs recorded by high-density electroencephalography (hd-EEG) provide additional information with respect SEEG alone and to provide an open, curated dataset to allow for further exploration of their potential.

METHODS

The dataset encompasses SEEG and hd-EEG recordings simultaneously acquired during Single Pulse Electrical Stimulation (SPES) in drug-resistant epileptic patients (N = 36) in whom stimulations were delivered with different physical, geometrical, and topological parameters. Differences in CCEPs were assessed by amplitude, latency, and spectral measures.

RESULTS

While invasively and non-invasively recorded CCEPs were generally correlated, differences in pulse duration, angle and stimulated cortical area were better captured by hd-EEG. Further, intracranial stimulation evoked site-specific hd-EEG responses that reproduced the spectral features of EEG responses to transcranial magnetic stimulation (TMS). Notably, SPES, albeit unperceived by subjects, elicited scalp responses that were up to one order of magnitude larger than the responses typically evoked by sensory stimulation in awake humans.

CONCLUSIONS

CCEPs can be simultaneously recorded with SEEG and hd-EEG and the latter provides a reliable descriptor of the effects of SPES as well as a common reference to compare the whole-brain effects of intracortical stimulation to those of non-invasive transcranial or sensory stimulations in humans.

Tinnitus Perception in Light of a Parietal Operculo-Insular Involvement: A Review

In tinnitus literature, researchers have increasingly been advocating for a clearer distinction between tinnitus perception and tinnitus-related distress. In non-bothersome tinnitus, the perception itself can be more specifically investigated: this has provided a body of evidence, based on resting-state and activation fMRI protocols, highlighting the involvement of regions outside the conventional auditory areas, such as the right parietal operculum. Here, we aim to conduct a review of available investigations of the human parietal operculo-insular subregions conducted at the microscopic, mesoscopic, and macroscopic scales arguing in favor of an auditory-somatosensory cross-talk. Both the previous literature and new results on functional connectivity derived from cortico-cortical evoked potentials show that these subregions present a dense tissue of interconnections and a strong connectivity with auditory and somatosensory areas in the healthy brain. Disrupted integration processes between these modalities may thus result in erroneous perceptions, such as tinnitus. More precisely, we highlight the role of a subregion of the right parietal operculum, known as OP3 according to the Jülich atlas, in the integration of auditory and somatosensory representation of the orofacial muscles in the healthy population. We further discuss how a dysfunction of these muscles could induce hyperactivity in the OP3. The evidence of direct electrical stimulation of this area eliciting auditory hallucinations further suggests its involvement in tinnitus perception. Finally, a small number of neuroimaging studies of therapeutic interventions for tinnitus provide additional evidence of right parietal operculum involvement.

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

BACKGROUND

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Imaging the effective networks associated with cortical function through intracranial high-frequency stimulation

Direct electrical stimulation (DES) is considered to be the gold standard for mapping cortical function. A careful mapping of the eloquent cortex is key to successful resective or ablative surgeries, with a minimal postoperative deficit, for treatment of drug‐resistant epilepsy. There is accumulating evidence suggesting that not only local, but also remote activations play an equally important role in evoking clinical effects. By introducing a new intracranial stimulation paradigm and signal analysis methodology allowing to disambiguate EEG responses from stimulation artifacts we highlight the spatial extent of the networks associated with clinical effects. Our study includes 26 patients that underwent stereoelectroencephalographic investigations for drug‐resistant epilepsy, having 337 depth electrodes with 4,351 contacts sampling most brain structures. The routine high‐frequency electrical stimulation protocol for eloquent cortex mapping was altered in a subtle way, by alternating the polarity of the biphasic pulses in a train, causing the splitting the spectral lines of the artifactual components, exposing the underlying tissue response. By performing a frequency‐domain analysis of the EEG responses during DES we were able to capture remote activations and highlight the effect's network. By using standard intersubject averaging and a fine granularity HCP‐MMP parcellation, we were able to create local and distant connectivity maps for 614 stimulations evoking specific clinical effects. The clinical value of such maps is not only for a better understanding of the extent of the effects' networks guiding the invasive exploration, but also for understanding the spatial patterns of seizure propagation given the timeline of the seizure semiology.

A 32-Channel Time-Multiplexed Artifact-Aware Neural Recording System

This paper presents a low-power, low-noise microsystem for the recording of neural local field potentials or intracranial electroencephalographic signals. It features 32 time-multiplexed channels at the electrode interface and offers the possibility to spatially delta encode data to take advantage of the large correlation of signals captured from nearby channels. The circuit also implements a mixed-signal voltage-triggered auto-ranging algorithm which allows to attenuate large interferers in digital domain while preserving neural information. This effectively increases the system dynamic range and avoids the onset of saturation. A prototype, fabricated in a standard 180 nm CMOS process, has been experimentally verified in-vitro with cellular cultures of primary cortical neurons from mice. The system shows an integrated input-referred noise in the 0.5-200 Hz band of 1.4 μV_rms for a spot noise of about 85 nV /√{Hz}. The system draws 1.5 μW per channel from 1.2 V supply and obtains 71 dB + 26 dB dynamic range when the artifact-aware auto-ranging mechanism is enabled, without penalising other critical specifications such as crosstalk between channels or common-mode and power supply rejection ratios.

Energy-Efficient Integrated Circuit Solutions Toward Miniaturized Closed-Loop Neural Interface Systems

Miniaturized implantable devices play a crucial role in neural interfaces by monitoring and modulating neural activities on the peripheral and central nervous systems. Research efforts toward a compact wireless closed-loop system stimulating the nerve automatically according to the user's condition have been maintained. These systems have several advantages over open-loop stimulation systems such as reduction in both power consumption and side effects of continuous stimulation. Furthermore, a compact and wireless device consuming low energy alleviates foreign body reactions and risk of frequent surgical operations. Unfortunately, however, the miniaturized closed-loop neural interface system induces several hardware design challenges such as neural activity recording with severe stimulation artifact, real-time stimulation artifact removal, and energy-efficient wireless power delivery. Here, we will review recent approaches toward the miniaturized closed-loop neural interface system with integrated circuit (IC) techniques.

Single-pulse electrical stimulation methodology in freely moving rat

BACKGROUND

Cortico-cortical evoked potentials (CCEP) are becoming popular to infer brain connectivity and cortical excitability in implanted refractory epilepsy patients. Our goal was to transfer this methodology to the freely moving rodent.

NEW METHOD

CCEP were recorded on freely moving Sprague-Dawley rats, from cortical and subcortical areas using depth electrodes. Electrical stimulation was applied using 1 ms biphasic current pulse, cathodic first, at a frequency of 0.5 Hz, with intensities ranging from 0.2 to 0.8 mA. Data were then processed in a similar fashion to human clinical studies, which included epoch selection, artefact correction and smart averaging.

RESULTS

For a large range of tested intensities, we recorded CCEPs with very good signal to noise ratio and reproducibility between animals, without any behavioral modification. The CCEP were composed of different components according to recorded and stimulated sites, similarly to human recordings.

COMPARISON WITH EXISTING METHODS

We minimally adapted a clinically-motivated methodology to a freely moving rodent model to achieve high translational relevance of future preclinical studies.

CONCLUSIONS

Our results indicate that the CCEP methodology can be applied to freely moving rodents and transferred to preclinical research. This will be of interest to address various neuroscientific questions, in physiological and pathological conditions.

Optimal Multichannel Artifact Prediction and Removal for Neural Stimulation and Brain Machine Interfaces

Neural implants that deliver multi-site electrical stimulation to the nervous systems are no longer the last resort but routine treatment options for various neurological disorders. Multi-site electrical stimulation is also widely used to study nervous system function and neural circuit transformations. These technologies increasingly demand dynamic electrical stimulation and closed-loop feedback control for real-time assessment of neural function, which is technically challenging since stimulus-evoked artifacts overwhelm the small neural signals of interest. We report a novel and versatile artifact removal method that can be applied in a variety of settings, from single- to multi-site stimulation and recording and for current waveforms of arbitrary shape and size. The method capitalizes on linear electrical coupling between stimulating currents and recording artifacts, which allows us to estimate a multi-channel linear Wiener filter to predict and subsequently remove artifacts via subtraction. We confirm and verify the linearity assumption and demonstrate feasibility in a variety of recording modalities, including in vitro sciatic nerve stimulation, bilateral cochlear implant stimulation, and multi-channel stimulation and recording between the auditory midbrain and cortex. We demonstrate a vast enhancement in the recording quality with a typical artifact reduction of 25-40 dB. The method is efficient and can be scaled to arbitrary number of stimulus and recording sites, making it ideal for applications in large-scale arrays, closed-loop implants, and high-resolution multi-channel brain-machine interfaces.

Signal recovery from stimulation artifacts in intracranial recordings with dictionary learning

OBJECTIVE

Electrical stimulation of the human brain is commonly used for eliciting and inhibiting neural activity for clinical diagnostics, modifying abnormal neural circuit function for therapeutics, and interrogating cortical connectivity. However, recording electrical signals with concurrent stimulation results in dominant electrical artifacts that mask the neural signals of interest. Here we develop a method to reproducibly and robustly recover neural activity during concurrent stimulation. We concentrate on signal recovery across an array of electrodes without channel-wise fine-tuning of the algorithm. Our goal includes signal recovery with trains of stimulation pulses, since repeated, high-frequency pulses are often required to induce desired effects in both therapeutic and research domains. We have made all of our code and data publicly available.

APPROACH

We developed an algorithm that automatically detects templates of artifacts across many channels of recording, creating a dictionary of learned templates using unsupervised clustering. The artifact template that best matches each individual artifact pulse is subtracted to recover the underlying activity. To assess the success of our method, we focus on whether it extracts physiologically interpretable signals from real recordings.

MAIN RESULTS

We demonstrate our signal recovery approach on invasive electrophysiologic recordings from human subjects during stimulation. We show the recovery of meaningful neural signatures in both electrocorticographic (ECoG) arrays and deep brain stimulation (DBS) recordings. In addition, we compared cortical responses induced by the stimulation of primary somatosensory (S1) by natural peripheral touch, as well as motor cortex activity with and without concurrent S1 stimulation.

SIGNIFICANCE

Our work will enable future advances in neural engineering with simultaneous stimulation and recording.

Quantifying volume conducted potential using stimulation artefact in cortico-cortical evoked potentials

BACKGROUND

Cortico-cortical evoked potentials (CCEP) are a technique using low frequency stimulation to infer regions of cortical connectivity in patients undergoing Stereo-electroencephalographic (SEEG) monitoring for refractory epilepsy. Little attention has been given to volume conducted components of CCEP responses, and how they may inflate CCEP connectivity.

NEW METHOD

Using data from 37 SEEG-CCEPs patients, a novel method was developed to quantify stimulation artefact by measuring the peak-to-peak voltage difference in the first 10 ms after CCEP stimulation. Early responses to CCEP stimulation were also quantified by calculating the root mean square of the 10-100 ms period after each stimulation pulse. Both the early CCEP responses and amplitude of stimulation artefact were regressed by physical distance, stimulation waveform, stimulation intensity and tissue type to identify conduction related properties.

RESULTS

Both stimulation artefact and early responses were correlated strongly with the inverse square of the distance from the stimulating electrode. Once corrected for the inverse square distance from the electrode, stimulation artefact and CCEP responses showed a linear relationship, indicating a volume conducted component.

COMPARISON WITH EXISTING METHODS

This is the first study to use stimulation artefact to quantify volume conducted potentials, and is the first to quantify volume conducted potentials in SEEG. A single prior study utilizing electrocorticography has shown that parts of early CCEP responses are due to volume conduction.

CONCLUSIONS

The linear relationship between stimulation artefact amplitude and CCEP early responses, once corrected for distance, suggests that stimulation artefact can be used as a measure to quantify the volume conducted components.

A Method for Suppressing Electrical Stimulation Artifacts from Electromyography

When surface electromyography (EMG) signal is used in a real-time functional electrical stimulation (FES) system for feedback control, the artifact from electrical stimulation is a key challenge for EMG signal processing. To address this challenge, this study proposes a novel method to suppress stimulation artifacts in the EMG-driven closed-loop FES system. The proposed method is inspired by an experimental study that compares artifacts generated by electrical stimulations with different current intensities. It is found that (1) spikes of stimulation artifacts are susceptible to the current intensity and (2) tailing components are similar under different current intensities. Based on these observations, the proposed method combines the blanking and template subtracting strategies for suppressing stimulation artifact. The length of blanking window for suppressing the stimulation spike is adaptively determined by a spike detection algorithm and the first-order derivative analysis of signal. An autoregressive model is used to estimate the tailing part of stimulation artifact, which is an adaptive template for subtracting the artifact. The proposed method is evaluated on both semi-synthetic and experimental datasets. Verified on the semi-synthetic dataset, the proposed method achieves better performance than the classic blanking method. Validated on the experimental dataset, the proposed method substantially decreases the power of stimulation artifact in the EMG. These results indicate that the proposed method can effectively suppress the stimulation artifact while retains the useful EMG signal for an EMG-driven FES system.

Localization of Epileptogenic Zone Based on Cortico-Cortical Evoked Potential (CCEP): A Feature Extraction and Graph Theory Approach

Objective

Epilepsy is a chronic brain disease, which is prone to relapse and affects individuals of all ages worldwide, particularly the very young and elderly. Up to one-third of these patients are medically intractable and require resection surgery. However, the outcomes of epilepsy surgery rely upon the clear identification of epileptogenic zone (EZ). The combination of cortico-cortical evoked potential (CCEP) and electrocorticography (ECoG) provides an opportunity to observe the connectivity of human brain network and more comprehensive information that may help the clinicians localize the epileptogenic focus more precisely. However, there is no standard analysis method in the clinical application of CCEPs, especially for the quantitative analysis of abnormal connectivity of epileptic networks. The aim of this paper was to present an approach on the batch processing of CCEPs and provide information relating to the localization of EZ for clinical study.

Methods

Eight medically intractable epilepsy patients were included in this study. Each patient was implanted with subdural grid electrodes and electrical stimulations were applied directly to their cortex to induce CCEPs. After signal preprocessing, we constructed three effective brain networks at different spatial scales for each patient, regarding the amplitudes of CCEPs as the connection weights. Graph theory was then applied to analyze the brain network topology of epileptic patients, and the topological metrics of EZ and non-EZ (NEZ) were compared.

Results

The effective connectivity network reconstructed from CCEPs was asymmetric, both the number and the amplitudes of effective CCEPs decreased with increasing distance between stimulating and recording sites. Besides, the distribution of CCEP responses was associated with the locations of EZ which tended to have higher degree centrality (DC) and nodal shortest path length (NLP) than NEZ.

Conclusion

Our results indicated that the brain networks of epileptics were asymmetric and mainly composed of short-distance connections. The DC and NLP were highly consistent to the distribution of the EZ, and these topological parameters have great potential to be readily applied to the clinical localization of the EZ.

Diffusion geometry approach to efficiently remove electrical stimulation artifacts in intracranial electroencephalography

OBJECTIVE

Cortical oscillations, electrophysiological activity patterns, associated with cognitive functions and impaired in many psychiatric disorders can be observed in intracranial electroencephalography (iEEG). Direct cortical stimulation (DCS) may directly target these oscillations and may serve as therapeutic approaches to restore functional impairments. However, the presence of electrical stimulation artifacts in neurophysiological data limits the analysis of the effects of stimulation. Currently available methods suffer in performance in the presence of nonstationarity inherent in biological data.

APPROACH

Our algorithm, shape adaptive nonlocal artifact removal (SANAR) is based on unsupervised manifold learning. By estimating the Euclidean median of k-nearest neighbors of each artifact in a nonlocal fashion, we obtain a faithful representation of the artifact which is then subtracted. This approach overcomes the challenges presented by nonstationarity.

MAIN RESULTS

SANAR is effective in removing stimulation artifacts in the time domain while preserving the spectral content of the endogenous neurophysiological signal. We demonstrate the performance in a simulated dataset as well as in human iEEG data. Using two quantitative measures, that capture how much of information from endogenous activity is retained, we demonstrate that SANAR's performance exceeds that of one of the widely used approaches, independent component analysis, in the time domain as well as the frequency domain.

SIGNIFICANCE

This approach allows for the analysis of iEEG data, single channel or multiple channels, during DCS, a crucial step in advancing our understanding of the effects of periodic stimulation and developing new therapies.

A multi-channel approach for cortical stimulation artefact suppression in depth EEG signals using time-frequency and spatial filtering

OBJECTIVE

The stereo electroencephalogram (SEEG) recordings are the sate of the art tool used in pre-surgical evaluation of drug-unresponsive epileptic patients. Coupled with SEEG, electrical cortical stimulation (CS) offer a complementary tool to investigate the lesioned/healthy brain regions and to identify the epileptic zones with precision. However, the propagation of this stimulation inside the brain masks the cerebral activity recorded by nearby multi-contact SEEG electrodes. The objective of this paper is to propose a novel filtering approach for suppressing the CS artifact in SEEG signals using time, frequency as well as spatial information.

METHODS

The method combines spatial filtering with tunable-Q wavelet transform (TQWT). SEEG signals are spatially filtered to isolate the CS artifacts within a few number of sources/components. The artifacted components are then decomposed into oscillatory background and sharp varying transient signals using tunable-Q wavelet transform (TQWT). The CS artifact is assumed to lie in the transient part of the signal. Using prior known time-frequency information of the CS artifacts, we selectively mask the wavelet coefficients of the transient signal and extract out any remaining significant electrophysiological activity.

RESULTS

We have applied our proposed method of CS artifact suppression on simulated and real SEEG signals with convincing performance. The experimental results indicate the effectiveness of the proposed approach.

CONCLUSION

The proposed method suppresses CS artifacts without affecting the background SEEG signal.

SIGNIFICANCE

The proposed method can be applied for suppressing both low and high frequency CS artifacts and outperforms current methods from the literature.

Monitoring Corticocortical Evoked Potentials During Intracranial Vascular Surgery

BACKGROUND

Monitoring of corticocortical evoked potentials (CCEPs) during brain tumor surgery of patients under anesthesia was recently reported to be effective in assisting in preservation of speech function. The aim of this study was to investigate whether CCEPs can be reproducibly measured between the frontal and temporal lobes during standard intracranial vascular surgery under general anesthesia; whether dynamic changes in CCEPs caused by reduced focal cerebral blood flow can be measured; and whether CCEPs can be used to monitor speech function, particularly associated with the left side of the brain.

METHODS

We monitored CCEPs during 58 vascular surgeries (42 clipping procedures; 15 bypasses, 1 of which overlapped with clipping; and 2 hematoma removals from the left frontal and temporal lobe) at Kashiwaba Neurosurgical Hospital from October 2016 to January 2018.

RESULTS

CCEPs could be reproducibly and routinely monitored in bilateral vascular surgeries. None of the patients experienced any postoperative symptoms or showed any ischemic lesions on postoperative magnetic resonance imaging; however, 5 patients temporarily demonstrated reduced CCEPs intraoperatively that were caused by transient obstructions of blood flow. Motor evoked potentials and somatosensory evoked potentials were simultaneously monitored intraoperatively and did not show any changes.

CONCLUSIONS

The results of our pilot study show that CCEPs can be routinely monitored during bilateral intracranial vascular surgery and that they are sensitive to ischemia. CCEPs on the left side could serve as unique intraoperative monitoring of speech function under anesthesia.

A quantitative method for evaluating cortical responses to electrical stimulation

BACKGROUND

Electrical stimulation of the cortex using subdurally implanted electrodes can causally reveal structural connectivity by eliciting cortico-cortical evoked potentials (CCEPs). While many studies have demonstrated the potential value of CCEPs, the methods to evaluate them were often relatively subjective, did not consider potential artifacts, and did not lend themselves to systematic scientific investigations.

NEW METHOD

We developed an automated and quantitative method called SIGNI (Stimulation-Induced Gamma-based Network Identification) to evaluate cortical population-level responses to electrical stimulation that minimizes the impact of electrical artifacts. We applied SIGNI to electrocorticographic (ECoG) data from eight human subjects who were implanted with a total of 978 subdural electrodes. Across the eight subjects, we delivered 92 trains of approximately 200 discrete electrical stimuli each (amplitude 4-15 mA) to a total of 64 electrode pairs.

RESULTS

We verified SIGNI's efficacy by demonstrating a relationship between the magnitude of evoked cortical activity and stimulation amplitude, as well as between the latency of evoked cortical activity and the distance from the stimulated locations.

CONCLUSIONS

SIGNI reveals the timing and amplitude of cortical responses to electrical stimulation as well as the structural connectivity supporting these responses. With these properties, it enables exploration of new and important questions about the neurophysiology of cortical communication and may also be useful for pre-surgical planning.

Probabilistic functional tractography of the human cortex revisited

In patients with pharmaco-resistant focal epilepsies investigated with intracranial electroencephalography (iEEG), direct electrical stimulations of a cortical region induce cortico-cortical evoked potentials (CCEP) in distant cerebral cortex, which properties can be used to infer large scale brain connectivity. In 2013, we proposed a new probabilistic functional tractography methodology to study human brain connectivity. We have now been revisiting this method in the F-TRACT project (f-tract.eu) by developing a large multicenter CCEP database of several thousand stimulation runs performed in several hundred patients, and associated processing tools to create a probabilistic atlas of human cortico-cortical connections. Here, we wish to present a snapshot of the methods and data of F-TRACT using a pool of 213 epilepsy patients, all studied by stereo-encephalography with intracerebral depth electrodes. The CCEPs were processed using an automated pipeline with the following consecutive steps: detection of each stimulation run from stimulation artifacts in raw intracranial EEG (iEEG) files, bad channels detection with a machine learning approach, model-based stimulation artifact correction, robust averaging over stimulation pulses. Effective connectivity between the stimulated and recording areas is then inferred from the properties of the first CCEP component, i.e. onset and peak latency, amplitude, duration and integral of the significant part. Finally, group statistics of CCEP features are implemented for each brain parcel explored by iEEG electrodes. The localization (coordinates, white/gray matter relative positioning) of electrode contacts were obtained from imaging data (anatomical MRI or CT scans before and after electrodes implantation). The iEEG contacts were repositioned in different brain parcellations from the segmentation of patients' anatomical MRI or from templates in the MNI coordinate system. The F-TRACT database using the first pool of 213 patients provided connectivity probability values for 95% of possible intrahemispheric and 56% of interhemispheric connections and CCEP features for 78% of intrahemisheric and 14% of interhemispheric connections. In this report, we show some examples of anatomo-functional connectivity matrices, and associated directional maps. We also indicate how CCEP features, especially latencies, are related to spatial distances, and allow estimating the velocity distribution of neuronal signals at a large scale. Finally, we describe the impact on the estimated connectivity of the stimulation charge and of the contact localization according to the white or gray matter. The most relevant maps for the scientific community are available for download on f-tract. eu (David et al., 2017) and will be regularly updated during the following months with the addition of more data in the F-TRACT database. This will provide an unprecedented knowledge on the dynamical properties of large fiber tracts in human.

Recovery of early neural spikes from stimulation electrodes using a DC-coupled low gain high resolution data acquisition system

BACKGROUND

Neural responses to electrical stimulation provide valuable information to probe and study the network function. Especially, recording neural responses from the stimulated site provides improved neural interfacing method. However, it is difficult to measure short-delayed responses at the stimulated electrode due to the saturation of the amplifier after stimulation which is called "stimulus artifact". Despite the advances in handling stimulation artifacts, it is still very challenging to deal with the artifacts if one tries to stimulate and record from the same electrode.

NEW METHOD

In this paper, we developed a system consisting of 24 bit ADC and low gain DC-amplifier which allows us to record the entire responses including saturation-free stimulus artifact and neural responses with excellent resolution.

RESULTS

Our approach showed saturation-free recording after stimulation, which makes it possible to recover neural spike as early as in 2 ms at the stimulating electrode with digital elimination methods.

COMPARISON WITH EXISTING METHODS

With our system we could record neural signals after stimulation that was difficult with high gain and high pass filtered recording system due to amplifier saturation.

CONCLUSIONS

Our new system can enhance interface performance with its higher robustness and with simple system configuration.

Discrimination of a medial functional module within the temporal lobe using an effective connectivity model: A CCEP study

The temporal lobe is classically divided in two functional systems: the ventral visual pathway and the medial temporal memory system. However, their functional separation has been challenged by studies suggesting that the medial temporal lobe could be best understood as an extension of the hierarchically organized ventral visual pathway. Our purpose was to investigate (i) whether cerebral regions within the temporal lobe could be grouped into distinct functional assemblies, and (ii) which regions were central within these functional assemblies. We studied low intensity and low frequency electrical stimulations (0.5 mA, 1 Hz, 4 ms) performed during sixteen pre-surgical intracerebral EEG investigations in patients with medically intractable temporal or temporo-occipital lobe epilepsies. Eleven regions of interest were delineated per anatomical landmarks such as gyri and sulci. Effective connectivity based on electrophysiological feature (amplitude) of cortico-cortical evoked potentials (CCEPs) was evaluated and subjected to graph metrics. The amplitudes discriminated one medial module where the hippocampus could act as a signal amplifier. Mean amplitudes of CCEPs in regions of the temporal lobe showed a generalized Pareto distribution of probability suggesting neural synchronies to be self-organized critically. Our description of effective interactions within the temporal lobe provides a regional electrophysiological model of effective connectivity which is discussed in the context of the current hypothesis of pattern completion.