Virtual stimulation of the interictal EEG network localizes the EZ as a measure of cortical excitability.

Introduction: For patients with drug-resistant epilepsy, successful localization and surgical treatment of the epileptogenic zone (EZ) can bring seizure freedom. However, surgical success rates vary widely because there are currently no clinically validated biomarkers of the EZ. Highly epileptogenic regions often display increased levels of cortical excitability, which can be probed using single-pulse electrical stimulation (SPES), where brief pulses of electrical current are delivered to brain tissue. It has been shown that high-amplitude responses to SPES can localize EZ regions, indicating a decreased threshold of excitability. However, performing extensive SPES in the epilepsy monitoring unit (EMU) is time-consuming. Thus, we built patient-specific in silico dynamical network models from interictal intracranial EEG (iEEG) to test whether virtual stimulation could reveal information about the underlying network to identify highly excitable brain regions similar to physical stimulation of the brain. Methods: We performed virtual stimulation in 69 patients that were evaluated at five centers and assessed for clinical outcome 1 year post surgery. We further investigated differences in observed SPES iEEG responses of 14 patients stratified by surgical outcome. Results: Clinically-labeled EZ cortical regions exhibited higher excitability from virtual stimulation than non-EZ regions with most significant differences in successful patients and little difference in failure patients. These trends were also observed in responses to extensive SPES performed in the EMU. Finally, when excitability was used to predict whether a channel is in the EZ or not, the classifier achieved an accuracy of 91%. Discussion: This study demonstrates how excitability determined via virtual stimulation can capture valuable information about the EZ from interictal intracranial EEG.

Network excitability of stimulation-induced spectral responses helps localize the seizure onset zone.

OBJECTIVE: While evoked potentials elicited by single pulse electrical stimulation (SPES) may assist seizure onset zone (SOZ) localization during intracranial EEG (iEEG) monitoring, induced high frequency activity has also shown promising utility. We aimed to predict SOZ sites using induced cortico-cortical spectral responses (CCSRs) as an index of excitability within epileptogenic networks. METHODS: SPES was conducted in 27 epilepsy patients undergoing iEEG monitoring and CCSRs were quantified by significant early (10-200 ms) increases in power from 10 to 250 Hz. Using response power as CCSR network connection strengths, graph centrality measures (metrics quantifying each site's influence within the network) were used to predict whether sites were within the SOZ. RESULTS: Across patients with successful surgical outcomes, greater CCSR centrality predicted SOZ sites and SOZ sites targeted for surgical treatment with median AUCs of 0.85 and 0.91, respectively. We found that the alignment between predicted and targeted SOZ sites predicted surgical outcome with an AUC of 0.79. CONCLUSIONS: These findings indicate that network analysis of CCSRs can be used to identify increased excitability of SOZ sites and discriminate important surgical targets within the SOZ. SIGNIFICANCE: CCSRs may supplement traditional passive iEEG monitoring in seizure localization, potentially reducing the need for recording numerous seizures.

CARLA: Adjusted common average referencing for cortico-cortical evoked potential data.

Human brain connectivity can be mapped by single pulse electrical stimulation during intracranial EEG measurements. The raw cortico-cortical evoked potentials (CCEP) are often contaminated by noise. Common average referencing (CAR) removes common noise and preserves response shapes but can introduce bias from responsive channels. We address this issue with an adjusted, adaptive CAR algorithm termed "CAR by Least Anticorrelation (CARLA)". CARLA was tested on simulated CCEP data and real CCEP data collected from four human participants. In CARLA, the channels are ordered by increasing mean cross-trial covariance, and iteratively added to the common average until anticorrelation between any single channel and all re-referenced channels reaches a minimum, as a measure of shared noise. We simulated CCEP data with true responses in 0-45 of 50 total channels. We quantified CARLA's error and found that it erroneously included 0 (median) truly responsive channels in the common average with ≤42 responsive channels, and erroneously excluded ≤2.5 (median) unresponsive channels at all responsiveness levels. On real CCEP data, signal quality was quantified with the mean R2 between all pairs of channels, which represents inter-channel dependency and is low for well-referenced data. CARLA re-referencing produced significantly lower mean R2 than standard CAR, CAR using a fixed bottom quartile of channels by covariance, and no re-referencing. CARLA minimizes bias in re-referenced CCEP data by adaptively selecting the optimal subset of non-responsive channels. It showed high specificity and sensitivity on simulated CCEP data and lowered inter-channel dependency compared to CAR on real CCEP data.

A biophysically constrained brain connectivity model based on stimulation-evoked potentials.

BACKGROUND: Single-pulse electrical stimulation (SPES) is an established technique used to map functional effective connectivity networks in treatment-refractory epilepsy patients undergoing intracranial-electroencephalography monitoring. While the connectivity path between stimulation and recording sites has been explored through the integration of structural connectivity, there are substantial gaps, such that new modeling approaches may advance our understanding of connectivity derived from SPES studies. NEW METHOD: Using intracranial electrophysiology data recorded from a single patient undergoing stereo-electroencephalography (sEEG) evaluation, we employ an automated detection method to identify early response components, C1, from pulse-evoked potentials (PEPs) induced by SPES. C1 components were utilized for a novel topology optimization method, modeling 3D electrical conductivity to infer neural pathways from stimulation sites. Additionally, PEP features were compared with tractography metrics, and model results were analyzed with respect to anatomical features. RESULTS: The proposed optimization model resolved conductivity paths with low error. Specific electrode contacts displaying high error correlated with anatomical complexities. The C1 component strongly correlated with additional PEP features and displayed stable, weak correlations with tractography measures. COMPARISON WITH EXISTING METHOD: Existing methods for estimating neural signal pathways are imaging-based and thus rely on anatomical inferences. CONCLUSIONS: These results demonstrate that informing topology optimization methods with human intracranial SPES data is a feasible method for generating 3D conductivity maps linking electrical pathways with functional neural ensembles. PEP-estimated effective connectivity is correlated with but distinguished from structural connectivity. Modeled conductivity resolves connectivity pathways in the absence of anatomical priors.

The effect of propofol on effective brain networks.

OBJECTIVE: We compared the effective networks derived from Single Pulse Electrical Stimulation (SPES) in intracranial electrocorticography (ECoG) of awake epilepsy patients and while under general propofol-anesthesia to investigate the effect of propofol on these brain networks. METHODS: We included nine patients who underwent ECoG for epilepsy surgery evaluation. We performed SPES when the patient was awake (SPES-clinical) and repeated this under propofol-anesthesia during the surgery in which the ECoG grids were removed (SPES-propofol). We detected the cortico-cortical evoked potentials (CCEPs) with an automatic detector. We constructed two effective networks derived from SPES-clinical and SPES-propofol. We compared three network measures (indegree, outdegree and betweenness centrality), the N1-peak-latency and amplitude of CCEPs between the two effective networks. RESULTS: Fewer CCEPs were observed during SPES-propofol (median: 6.0, range: 0-29) compared to SPES-clinical (median: 10.0, range: 0-36). We found a significant correlation for the indegree, outdegree and betweenness centrality between SPES-clinical and SPES-propofol (respectively rs = 0.77, rs = 0.70, rs = 0.55, p < 0.001). The median N1-peak-latency increased from 22.0 ms during SPES-clinical to 26.4 ms during SPES-propofol. CONCLUSIONS: Our findings suggest that the number of effective network connections decreases, but network measures are only marginally affected. SIGNIFICANCE: The primary network topology is preserved under propofol.

Effective connectivity relates seizure outcome to electrode placement in responsive neurostimulation.

Responsive neurostimulation is a closed-loop neuromodulation therapy for drug resistant focal epilepsy. Responsive neurostimulation electrodes are placed near ictal onset zones so as to enable detection of epileptiform activity and deliver electrical stimulation. There is no standard approach for determining the optimal placement of responsive neurostimulation electrodes. Clinicians make this determination based on presurgical tests, such as MRI, EEG, magnetoencephalography, ictal single-photon emission computed tomography and intracranial EEG. Currently functional connectivity measures are not being used in determining the placement of responsive neurostimulation electrodes. Cortico-cortical evoked potentials are a measure of effective functional connectivity. Cortico-cortical evoked potentials are generated by direct single-pulse electrical stimulation and can be used to investigate cortico-cortical connections in vivo. We hypothesized that the presence of high amplitude cortico-cortical evoked potentials, recorded during intracranial EEG monitoring, near the eventual responsive neurostimulation contact sites is predictive of better outcomes from its therapy. We retrospectively reviewed 12 patients in whom cortico-cortical evoked potentials were obtained during stereoelectroencephalography evaluation and subsequently underwent responsive neurostimulation therapy. We studied the relationship between cortico-cortical evoked potentials, the eventual responsive neurostimulation electrode locations and seizure reduction. Directional connectivity indicated by cortico-cortical evoked potentials can categorize stereoelectroencephalography electrodes as either receiver nodes/in-degree (an area of greater inward connectivity) or projection nodes/out-degree (greater outward connectivity). The follow-up period for seizure reduction ranged from 1.3-4.8 years (median 2.7) after responsive neurostimulation therapy started. Stereoelectroencephalography electrodes closest to the eventual responsive neurostimulation contact site tended to show larger in-degree cortico-cortical evoked potentials, especially for the early latency cortico-cortical evoked potentials period (10-60 ms period) in six out of 12 patients. Stereoelectroencephalography electrodes closest to the responsive neurostimulation contacts (≤5 mm) also had greater significant out-degree in the early cortico-cortical evoked potentials latency period than those further away (≥10 mm) (P < 0.05). Additionally, significant correlation was noted between in-degree cortico-cortical evoked potentials and greater seizure reduction with responsive neurostimulation therapy at its most effective period (P < 0.05). These findings suggest that functional connectivity determined by cortico-cortical evoked potentials may provide additional information that could help guide the optimal placement of responsive neurostimulation electrodes.

Delayed effective connectivity characterizes the epileptogenic zone during stereo-EEG.

OBJECTIVE: Single-pulse electrical stimulations (SPES) can elicit normal and abnormal responses that might characterize the epileptogenic zone, including spikes, high-frequency oscillations and cortico-cortical evoked potentials (CCEPs). In this study, we investigate their association with the epileptogenic zone during stereoelectroencephalography (SEEG) in 28 patients with refractory focal epilepsy. METHODS: Characteristics of CCEPs (distance-corrected or -uncorrected latency, amplitude and the connectivity index) and the occurrence of spikes and ripples were assessed. Responses within the epileptogenic zone and within the non-involved zone were compared using receiver operating characteristics curves and analysis of variance (ANOVA) either in all patients, patients with well-delineated epileptogenic zone, and patients older than 15 years old. RESULTS: We found an increase in distance-corrected CCEPs latency after stimulation within the epileptogenic zone (area under the curve = 0.71, 0.72, 0.70, ANOVA significant after false discovery rate correction). CONCLUSIONS: The increased distance-corrected CCEPs latency suggests that neuronal propagation velocity is altered within the epileptogenic network. This association might reflect effective connectivity changes at cortico-cortical or cortico-subcortico-cortical levels. Other responses were not associated with the epileptogenic zone, including the CCEPs amplitude, the connectivity index, the occurrences of induced ripples and spikes. The discrepancy with previous descriptions may be explained by different spatial brain sampling between subdural and depth electrodes. SIGNIFICANCE: Increased distance-corrected CCEPs latency, indicating delayed effective connectivity, characterizes the epileptogenic zone. This marker could be used to help tailor surgical resection limits after SEEG.

SEEG electrode shaft affects amplitude and latency of potentials evoked with single pulse electrical stimulation.

BACKGROUND: Long and thin shaft electrodes are implanted intracerebrally for stereoelectroencephalography (SEEG) in patients with pharmacoresistant focal epilepsies. Two adjacent contacts of one of such electrodes can deliver a train of single pulse electrical stimulations (SPES), and evoked potentials (EPs) are recorded on other contacts. In this study we assess if stimulating and recording on the same shaft, as opposed to different shafts, has an impact on common EP features. NEW METHOD: We leverage the large volume of SEEG data gathered in the F-TRACT database and analyze data from nearly one thousand SEEG implantations in order to verify whether stimulation and recording from the same shaft influence the EP pattern. RESULTS: We found that when the stimulated and the recording contacts were located on the same shaft, the mean and median amplitudes of an EP are greater, and its mean and median latencies are smaller than when the contacts were located on different shafts. This effect is small (Cohen's d ∼ 0.1), but robust (p-value < 10-3) across the SEEG database. COMPARISON WITH EXISTING METHOD(S): Our study is the first one to address this question. Due to the choice of commonly used EP features, our method is congruent with other studies. CONCLUSIONS: The magnitude of the reported effect does not obligate all standard analyses to correct for it, unless they aim at high precision. The source of the effect is not clear. Manufacturers of SEEG electrodes could examine it and potentially minimize the effect in their future products.

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.

A biophysically constrained brain connectivity model based on stimulation-evoked potentials.

Background: Single-pulse electrical stimulation (SPES) is an established technique used to map functional effective connectivity networks in treatment-refractory epilepsy patients undergoing intracranial-electroencephalography monitoring. While the connectivity path between stimulation and recording sites has been explored through the integration of structural connectivity, there are substantial gaps, such that new modeling approaches may advance our understanding of connectivity derived from SPES studies. New Method: Using intracranial electrophysiology data recorded from a single patient undergoing sEEG evaluation, we employ an automated detection method to identify early response components, C1, from pulse-evoked potentials (PEPs) induced by SPES. C1 components were utilized for a novel topology optimization method, modeling 3D conductivity propagation from stimulation sites. Additionally, PEP features were compared with tractography metrics, and model results were analyzed with respect to anatomical features. Results: The proposed optimization model resolved conductivity paths with low error. Specific electrode contacts displaying high error correlated with anatomical complexities. The C1 component strongly correlates with additional PEP features and displayed stable, weak correlations with tractography measures. Comparison with existing methods: Existing methods for estimating conductivity propagation are imaging-based and thus rely on anatomical inferences. Conclusions: These results demonstrate that informing topology optimization methods with human intracranial SPES data is a feasible method for generating 3D conductivity maps linking electrical pathways with functional neural ensembles. PEP-estimated effective connectivity is correlated with but distinguished from structural connectivity. Modeled conductivity resolves connectivity pathways in the absence of anatomical priors.

Localization of Epileptic Brain Responses to Single-Pulse Electrical Stimulation by Developing an Adaptive Iterative Linearly Constrained Minimum Variance Beamformer.

Delayed responses (DRs) to single pulse electrical stimulation (SPES) in patients with severe refractory epilepsy, from their intracranial recordings, can help to identify regions associated with epileptogenicity. Automatic DR localization is a large step in speeding up the identification of epileptogenic focus. Here, for the first time, an adaptive iterative linearly constrained minimum variance beamformer (AI-LCMV) is developed and employed to localize the DR sources from intracranial electroencephalogram (EEG) recorded using subdural electrodes. The prime objective here is to accurately localize the regions for the corresponding DRs using an adaptive localization method that exploits the morphology of DRs as the desired sources. The traditional closed-form linearly constrained minimum variance (CF-LCMV) solution is meant for tracking the sources with dominating power. Here, by incorporating the morphology of DRs, as a constraint, to an iterative linearly constrained minimum variance (LCMV) solution, the array of subdural electrodes is used to localize the low-power DRs, some not even visible in any of the electrode signals. The results from the cases included in this study also indicate more distinctive locations compared to those achievable by conventional beamformers. Most importantly, the proposed AI-LCMV is able to localize the DRs invisible over other electrodes.

Evoked responses to single pulse electrical stimulation reveal impaired striatal excitability in a rat model of Parkinson's disease.

BACKGROUND: Sensorimotor beta oscillations are increased in Parkinson's disease (PD) due to the alteration of dopaminergic transmission. This electrophysiological read-out is reported both in patients and in animal models such as the 6-OHDA rat model obtained with unilateral nigral injection of 6-hydroxydopamine (6-OHDA). Current treatments, based on dopaminergic replacement, transiently normalize this pathological beta activity and improve patients' quality of life. OBJECTIVES: We wanted to assess in vivo whether the abnormal beta oscillations can be correlated with impaired striatal or cortical excitability of the sensorimotor system and modulated by the pharmacological manipulation of the dopaminergic system. METHODS: In the unilateral 6-OHDA rat model and control animals, we used intra-striatal and intra-cortical single-pulse electrical stimulation (SPES) and concurrent local field potentials (LFP) recordings. In the two groups, we quantified basal cortico-striatal excitability from time-resolved spectral analyses of LFP evoked responses induced remotely by intracerebral stimulations. The temporal dependance of cortico-striatal excitability to dopaminergic transmission was further tested using electrophysiological recordings combined with levodopa injection. RESULTS: LFP evoked responses after striatal stimulation showed a transient reduction of power in a large time-frequency domain in the 6-OHDA group compared to the sham group. This result was specific to the striatum, as no significant difference was observed in cortical LFP evoked responses between the two groups. This impaired striatal excitability in the 6-OHDA group was observed in the striatum at least during the first 3 months after the initial lesion. In addition, the striatum responses to SPES during a levodopa challenge showed a transient potentiation of the decrease of responsiveness in frequencies below 40 Hz. CONCLUSION: The spectral properties of striatal responses to SPES show high sensitivity to dopaminergic transmission in the unilateral 6-OHDA rat model. We thus propose that this approach could be used in preclinical models as a time-resolved biomarker of impaired dopaminergic transmission capable of monitoring progressive neurodegeneration and/or challenges to drug intake.

Towards optimizing single pulse electrical stimulation: High current intensity, short pulse width stimulation most effectively elicits evoked potentials.

BACKGROUND: While single pulse electrical stimulation (SPES) is increasingly used to study effective connectivity, the effects of varying stimulation parameters on the resulting cortico-cortical evoked potentials (CCEPs) have not been systematically explored. OBJECTIVE: We sought to understand the interacting effects of stimulation pulse width, current intensity, and charge on CCEPs through an extensive testing of this parameter space and analysis of several response metrics. METHODS: We conducted SPES in 11 patients undergoing intracranial EEG monitoring using five combinations of current intensity (1.5, 2.0, 3.0, 5.0, and 7.5 mA) and pulse width at each of three charges (0.750, 1.125, and 1.500 µC/phase) to study how CCEP amplitude, distribution, latency, morphology, and stimulus artifact amplitude vary with each parameter. RESULTS: Stimulations with a greater charge or a greater current intensity and shorter pulse width at a given charge generally resulted in greater CCEP amplitudes and spatial distributions, shorter latencies, and increased waveform correlation. These effects interacted such that stimulations with the lowest charge and highest current intensities resulted in greater response amplitudes and spatial distributions than stimulations with the highest charge and lowest current intensities. Stimulus artifact amplitude increased with charge, but this could be mitigated by using shorter pulse widths. CONCLUSIONS: Our results indicate that individual combinations of current intensity and pulse width, in addition to charge, are important determinants of CCEP magnitude, morphology, and spatial extent. Together, these findings suggest that high current intensity, short pulse width stimulations are optimal SPES settings for eliciting strong and consistent responses while minimizing charge.

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.

Localizing seizure onset zones in surgical epilepsy with neurostimulation deep learning.

OBJECTIVE: In drug-resistant temporal lobe epilepsy, automated tools for seizure onset zone (SOZ) localization that use brief interictal recordings could supplement presurgical evaluations and improve care. Thus, the authors sought to localize SOZs by training a multichannel convolutional neural network on stereoelectroencephalography (SEEG) cortico-cortical evoked potentials. METHODS: The authors performed single-pulse electrical stimulation in 10 drug-resistant temporal lobe epilepsy patients implanted with SEEG. Using 500,000 unique poststimulation SEEG epochs, the authors trained a multichannel 1-dimensional convolutional neural network to determine whether an SOZ had been stimulated. RESULTS: SOZs were classified with mean sensitivity of 78.1% and specificity of 74.6% according to leave-one-patient-out testing. To achieve maximum accuracy, the model required a 0- to 350-msec poststimulation time period. Post hoc analysis revealed that the model accurately classified unilateral versus bilateral mesial temporal lobe seizure onset, as well as neocortical SOZs. CONCLUSIONS: This was the first demonstration, to the authors' knowledge, that a deep learning framework can be used to accurately classify SOZs with single-pulse electrical stimulation-evoked responses. These findings suggest that accurate classification of SOZs relies on a complex temporal evolution of evoked responses within 350 msec of stimulation. Validation in a larger data set could provide a practical clinical tool for the presurgical evaluation of drug-resistant epilepsy.

Cortico-cortical evoked potentials in response to varying stimulation intensity improves seizure localization.

OBJECTIVE: As single pulse electrical stimulation (SPES) is increasingly utilized to help localize the seizure onset zone (SOZ), it is important to understand how stimulation intensity can affect the ability to use cortico-cortical evoked potentials (CCEPs) to delineate epileptogenic regions. METHODS: We studied 15 drug-resistant epilepsy patients undergoing intracranial EEG monitoring and SPES with titrations of stimulation intensity. The N1 amplitude and distribution of CCEPs elicited in the SOZ and non-seizure onset zone (nSOZ) were quantified at each intensity. The separability of the SOZ and nSOZ using N1 amplitudes was compared between models using responses to titrations, responses to one maximal intensity, or both. RESULTS: At 2 mA and above, the increase in N1 amplitude with current intensity was greater for responses within the SOZ, and SOZ response distribution was maximized by 4-6 mA. Models incorporating titrations achieved better separability of SOZ and nSOZ compared to those using one maximal intensity. CONCLUSIONS: We demonstrated that differences in CCEP amplitude over a range of current intensities can improve discriminability of SOZ regions. SIGNIFICANCE: This study provides insight into the underlying excitability of the SOZ and how differences in current-dependent amplitudes of CCEPs may be used to help localize epileptogenic sites.

Time-frequency signatures evoked by single-pulse deep brain stimulation to the subcallosal cingulate.

Precision targeting of specific white matter bundles that traverse the subcallosal cingulate (SCC) has been linked to efficacy of deep brain stimulation (DBS) for treatment resistant depression (TRD). Methods to confirm optimal target engagement in this heterogenous region are now critical to establish an objective treatment protocol. As yet unexamined are the time-frequency features of the SCC evoked potential (SCC-EP), including spectral power and phase-clustering. We examined these spectral features-evoked power and phase clustering-in a sample of TRD patients (n = 8) with implanted SCC stimulators. Electroencephalogram (EEG) was recorded during wakeful rest. Location of electrical stimulation in the SCC target region was the experimental manipulation. EEG was analyzed at the surface level with an average reference for a cluster of frontal sensors and at a time window identified by prior study (50-150 ms). Morlet wavelets generated indices of evoked power and inter-trial phase clustering. Enhanced phase clustering at theta frequency (4-7 Hz) was observed in every subject and was significantly correlated with SCC-EP magnitude, but only during left SCC stimulation. Stimulation to dorsal SCC evinced stronger phase clustering than ventral SCC. There was a weak correlation between phase clustering and white matter density. An increase in evoked delta power (2-4 Hz) was also coincident with SCC-EP, but was less consistent across participants. DBS evoked time-frequency features index mm-scale changes to the location of stimulation in the SCC target region and correlate with structural characteristics implicated in treatment optimization. Results also imply a shared generative mechanism (inter-trial phase clustering) between evoked potentials evinced by electrical stimulation and evoked potentials evinced by auditory/visual stimuli and behavioral tasks. Understanding how current injection impacts downstream cortical activity is essential to building new technologies that adapt treatment parameters to individual differences in neurophysiology.

Stimulating native seizures with neural resonance: a new approach to localize the seizure onset zone.

Successful outcomes in epilepsy surgery rely on the accurate localization of the seizure onset zone. Localizing the seizure onset zone is often a costly and time-consuming process wherein a patient undergoes intracranial EEG monitoring, and a team of clinicians wait for seizures to occur. Clinicians then analyse the intracranial EEG before each seizure onset to identify the seizure onset zone and localization accuracy increases when more seizures are captured. In this study, we develop a new approach to guide clinicians to actively elicit seizures with electrical stimulation. We propose that a brain region belongs to the seizure onset zone if a periodic stimulation at a particular frequency produces large amplitude oscillations in the intracranial EEG network that propagate seizure activity. Such responses occur when there is 'resonance' in the intracranial EEG network, and the resonant frequency can be detected by observing a sharp peak in the magnitude versus frequency response curve, called a Bode plot. To test our hypothesis, we analysed single-pulse electrical stimulation response data in 32 epilepsy patients undergoing intracranial EEG monitoring. For each patient and each stimulated brain region, we constructed a Bode plot by estimating a transfer function model from the intracranial EEG 'impulse' or single-pulse electrical stimulation response. The Bode plots were then analysed for evidence of resonance. First, we showed that when Bode plot features were used as a marker of the seizure onset zone, it distinguished successful from failed surgical outcomes with an area under the curve of 0.83, an accuracy that surpassed current methods of analysis with cortico-cortical evoked potential amplitude and cortico-cortical spectral responses. Then, we retrospectively showed that three out of five native seizures accidentally triggered in four patients during routine periodic stimulation at a given frequency corresponded to a resonant peak in the Bode plot. Last, we prospectively stimulated peak resonant frequencies gleaned from the Bode plots to elicit seizures in six patients, and this resulted in an induction of three seizures and three auras in these patients. These findings suggest neural resonance as a new biomarker of the seizure onset zone that can guide clinicians in eliciting native seizures to more quickly and accurately localize the seizure onset zone.

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.

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.