Comparing connectivity metrics in cortico-cortical evoked potentials using synthetic cortical response patterns

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.

Macroscale traveling waves evoked by single-pulse stimulation of the human brain

Understanding the spatiotemporal dynamics of neural signal propagation is fundamental to unraveling the complexities of brain function. Emerging evidence suggests that cortico-cortical evoked potentials (CCEPs) resulting from single-pulse electrical stimulation may be used to characterize the patterns of information flow between and within brain networks. At present, the basic spatiotemporal dynamics of CCEP propagation cortically and subcortically are incompletely understood. We hypothesized that single-pulse electrical stimulation evokes neural traveling waves detectable in the three-dimensional space sampled by intracranial stereoelectroencephalography. Across a cohort of 21 adult patients with intractable epilepsy, we delivered 17,631 stimulation pulses and recorded CCEP responses in 1,019 electrode contacts. The distance between each pair of electrode contacts was approximated using three different metrics (Euclidean distance, path length, and geodesic distance), representing direct, tractographic, and transcortical propagation, respectively. For each robust CCEP, we extracted amplitude-, spectral-, and phase-based features to identify traveling waves emanating from the site of stimulation. Many evoked responses to stimulation appear to propagate as traveling waves (∼14-28%), despite sparse sampling throughout the brain. These stimulation-evoked traveling waves exhibited biologically plausible propagation velocities (range 0.1-9.6 m/s). Our results reveal that direct electrical stimulation elicits neural activity with variable spatiotemporal dynamics, including the initiation of neural traveling waves.

Significance Statement

Using single-pulse stimulation, we identify a subset of intracranial evoked potentials that propagate as neural traveling waves. Our results were robust across a range of distinct but complementary analysis methods. The identification of stimulation-evoked traveling waves may help to better characterize the pathways traversed by spontaneous, pathological, or task-evoked traveling waves and distinguish biologically plausible propagation from volume-conducted signals.

Data processing techniques impact quantification of cortico-cortical evoked potentials

BACKGROUND

Cortico-cortical evoked potentials (CCEPs) are a common tool for probing effective connectivity in intracranial human electrophysiology. As with all human electrophysiology data, CCEP data are highly susceptible to noise. To address noise, filters and re-referencing are often applied to CCEP data, but different processing strategies are used from study to study.

NEW METHOD

We systematically compare how common average re-referencing and filtering CCEP data impacts quantification.

RESULTS

We show that common average re-referencing and filters, particularly filters that cut out more frequencies, can significantly impact the quantification of CCEP magnitude and morphology. We identify that high cutoff high pass filters (> 0.5 Hz), low cutoff low pass filters (< 200 Hz), and common average re-referencing impact quantification across subjects. However, we also demonstrate that the presence of noise may impact CCEP quantification, and preprocessing is necessary to mitigate this. We show that filtering is more effective than re-referencing or averaging across trials for reducing most common types of noise.

COMPARISON WITH EXISTING METHODS

These results suggest that existing CCEP processing methods must be applied with care to maximize noise reduction and minimize changes to the data. We do not test every available processing strategy; rather we demonstrate that processing can influence the results of CCEP studies. We emphasize the importance of reporting all processing methods, particularly re-referencing methods.

CONCLUSIONS

We propose a general framework for choosing an appropriate processing pipeline for CCEP data, taking into consideration the noise levels of a specific dataset. We suggest that minimal gentle filtering is preferable.

Localizing seizure onset zone by a cortico-cortical evoked potentials-based machine learning approach in focal epilepsy

OBJECTIVE

We aimed to develop a new approach for identifying the localization of the seizure onset zone (SOZ) based on corticocortical evoked potentials (CCEPs) and to compare the connectivity patterns in patients with different clinical phenotypes.

METHODS

Fifty patients who underwent stereoelectroencephalography and CCEP procedures were included. Logistic regression was used in the model, and six CCEP metrics were input as features: root mean square of the first peak (N1RMS) and second peak (N2RMS), peak latency, onset latency, width duration, and area.

RESULTS

The area under the curve (AUC) for localizing the SOZ ranged from 0.88 to 0.93. The N1RMS values in the hippocampus sclerosis (HS) group were greater than that of the focal cortical dysplasia (FCD) IIa group (p < 0.001), independent of the distance between the recorded and stimulated sites. The sensitivity of localization was higher in the seizure-free group than in the non-seizure-free group (p = 0.036).

CONCLUSIONS

This new method can be used to predict the SOZ localization in various focal epilepsy phenotypes.

SIGNIFICANCE

This study proposed a machine-learning approach for localizing the SOZ. Moreover, we examined how clinical phenotypes impact large-scale abnormality of the epileptogenic networks.

An orexigenic subnetwork within the human hippocampus

Only recently have more specific circuit-probing techniques become available to inform previous reports implicating the rodent hippocampus in orexigenic appetitive processing1-4. This function has been reported to be mediated at least in part by lateral hypothalamic inputs, including those involving orexigenic lateral hypothalamic neuropeptides, such as melanin-concentrating hormone5,6. This circuit, however, remains elusive in humans. Here we combine tractography, intracranial electrophysiology, cortico-subcortical evoked potentials, and brain-clearing 3D histology to identify an orexigenic circuit involving the lateral hypothalamus and converging in a hippocampal subregion. We found that low-frequency power is modulated by sweet-fat food cues, and this modulation was specific to the dorsolateral hippocampus. Structural and functional analyses of this circuit in a human cohort exhibiting dysregulated eating behaviour revealed connectivity that was inversely related to body mass index. Collectively, this multimodal approach describes an orexigenic subnetwork within the human hippocampus implicated in obesity and related eating disorders.

The Interictal Suppression Hypothesis in focal epilepsy: network-level supporting evidence

Why are people with focal epilepsy not continuously having seizures? Previous neuronal signalling work has implicated gamma-aminobutyric acid balance as integral to seizure generation and termination, but is a high-level distributed brain network involved in suppressing seizures? Recent intracranial electrographic evidence has suggested that seizure-onset zones have increased inward connectivity that could be associated with interictal suppression of seizure activity. Accordingly, we hypothesize that seizure-onset zones are actively suppressed by the rest of the brain network during interictal states. Full testing of this hypothesis would require collaboration across multiple domains of neuroscience. We focused on partially testing this hypothesis at the electrographic network level within 81 individuals with drug-resistant focal epilepsy undergoing presurgical evaluation. We used intracranial electrographic resting-state and neurostimulation recordings to evaluate the network connectivity of seizure onset, early propagation and non-involved zones. We then used diffusion imaging to acquire estimates of white-matter connectivity to evaluate structure-function coupling effects on connectivity findings. Finally, we generated a resting-state classification model to assist clinicians in detecting seizure-onset and propagation zones without the need for multiple ictal recordings. Our findings indicate that seizure onset and early propagation zones demonstrate markedly increased inwards connectivity and decreased outwards connectivity using both resting-state (one-way ANOVA, P-value = 3.13 × 10-13) and neurostimulation analyses to evaluate evoked responses (one-way ANOVA, P-value = 2.5 × 10-3). When controlling for the distance between regions, the difference between inwards and outwards connectivity remained stable up to 80 mm between brain connections (two-way repeated measures ANOVA, group effect P-value of 2.6 × 10-12). Structure-function coupling analyses revealed that seizure-onset zones exhibit abnormally enhanced coupling (hypercoupling) of surrounding regions compared to presumably healthy tissue (two-way repeated measures ANOVA, interaction effect P-value of 9.76 × 10-21). Using these observations, our support vector classification models achieved a maximum held-out testing set accuracy of 92.0 ± 2.2% to classify early propagation and seizure-onset zones. These results suggest that seizure-onset zones are actively segregated and suppressed by a widespread brain network. Furthermore, this electrographically observed functional suppression is disproportionate to any observed structural connectivity alterations of the seizure-onset zones. These findings have implications for the identification of seizure-onset zones using only brief electrographic recordings to reduce patient morbidity and augment the presurgical evaluation of drug-resistant epilepsy. Further testing of the interictal suppression hypothesis can provide insight into potential new resective, ablative and neuromodulation approaches to improve surgical success rates in those suffering from drug-resistant focal epilepsy.

Human anterior thalamic stimulation evoked cortical potentials align with intrinsic functional connectivity

Characterizing human thalamocortical network is fundamental for understanding a vast array of human behaviors since the thalamus plays a central role in cortico-subcortical communication. Over the past few decades, advances in functional magnetic resonance imaging have allowed for spatial mapping of intrinsic resting-state functional connectivity (RSFC) between both cortical regions and in cortico-subcortical networks. Despite these advances, identifying the electrophysiological basis of human thalamocortical network architecture remains challenging. By leveraging stereoelectroencephalography electrodes temporarily implanted into distributed cortical regions and the anterior nucleus of the thalamus (ANT) of 10 patients with refractory focal epilepsy, we tested whether ANT stimulation evoked cortical potentials align with RSFC from the stimulation site, derived from a normative functional connectome (n=1000). Our study identifies spatial convergence of ANT stimulation evoked cortical potentials and normative RSFC. Other than connections to the Papez circuit, the ANT was found to be closely connected to several distinct higher-order association cortices, including the precuneus, angular gyrus, dorsal lateral prefrontal cortex, and anterior insula. Remarkably, we found that the spatial distribution and magnitude of cortical-evoked responses to single-pulse electrical stimulation of the ANT aligned with the spatial pattern and strength of normative RSFC of the stimulation site. The present study provides electrophysiological evidence that stimulation evoked electrical activity flows along intrinsic brain networks connected on a thalamocortical level.

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.

Canonical Response Parameterization: Quantifying the structure of responses to single-pulse intracranial electrical brain stimulation

Single-pulse electrical stimulation in the nervous system, often called cortico-cortical evoked potential (CCEP) measurement, is an important technique to understand how brain regions interact with one another. Voltages are measured from implanted electrodes in one brain area while stimulating another with brief current impulses separated by several seconds. Historically, researchers have tried to understand the significance of evoked voltage polyphasic deflections by visual inspection, but no general-purpose tool has emerged to understand their shapes or describe them mathematically. We describe and illustrate a new technique to parameterize brain stimulation data, where voltage response traces are projected into one another using a semi-normalized dot product. The length of timepoints from stimulation included in the dot product is varied to obtain a temporal profile of structural significance, and the peak of the profile uniquely identifies the duration of the response. Using linear kernel PCA, a canonical response shape is obtained over this duration, and then single-trial traces are parameterized as a projection of this canonical shape with a residual term. Such parameterization allows for dissimilar trace shapes from different brain areas to be directly compared by quantifying cross-projection magnitudes, response duration, canonical shape projection amplitudes, signal-to-noise ratios, explained variance, and statistical significance. Artifactual trials are automatically identified by outliers in sub-distributions of cross-projection magnitude, and rejected. This technique, which we call "Canonical Response Parameterization" (CRP) dramatically simplifies the study of CCEP shapes, and may also be applied in a wide range of other settings involving event-triggered data.

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.

Network of autoscopic hallucinations elicited by intracerebral stimulations of periventricular nodular heterotopia: An SEEG study

Periventricular nodular heterotopias (PVNH) are areas of neurons abnormally located in the white matter that might be involved in physiological cortical functions. Autoscopic hallucinations are changes in self-consciousness determined by a mismatch in integration of multiple sensory inputs. Our goal is to highlight the brain network involved in generation of autoscopic hallucination elicited by electrical stimulation of a PVNH in a drug resistant epilepsy patient. Our patient was explored using stereo-electroencephalography with electrodes covering the right posterior temporal PVNH and the adjacent cortex. Direct electrical high frequency stimulation of the PVNH elicited autoscopic hallucinations mainly involving the face and upper trunk. We then used multiple modalities to determine brain connectivity: single pulse electrical stimulation of the PVNH and stimulation-evoked potentials were used to highlight resting state effective connectivity. High-frequency stimulation using alternating polarity pulses enabled us to identify the network involved, time-locked to the clinical effect and to map symptom-related effective connectivity. Functional connectivity using a non-linear regression method was used to determine dependencies between different cortical regions following the stimulation. Finally, structural connectivity was highlighted using deterministic fiber tracking. Multi-modal connectivity analysis identified a network involving the PVNH, occipital and temporal neocortex, fusiform gyrus and parietal cortex.

Effects of stimulation intensity on intracranial cortico-cortical evoked potentials: A titration study

OBJECTIVE

The aim of the present study was to investigate the optimal stimulation parameters for eliciting cortico-cortical evoked potentials (CCEPs) for mapping functional and epileptogenic networks.

METHODS

We studied 13 patients with refractory epilepsy undergoing intracranial EEG monitoring. We systematically titrated the intensity of single-pulse electrical stimulation at multiple sites to assess the effect of increasing current on salient features of CCEPs such as N1 potential magnitude, signal to noise ratio, waveform similarity, and spatial distribution of responses. Responses at each incremental stimulation setting were compared to each other and to a final set of responses at the maximum intensity used in each patient (3.5-10 mA, median 6 mA).

RESULTS

We found that with a biphasic 0.15 ms/phase pulse at least 2-4 mA is needed to differentiate between non-responsive and responsive sites, and that stimulation currents of 6-7 mA are needed to maximize amplitude and spatial distribution of N1 responses and stabilize waveform morphology.

CONCLUSIONS

We determined a minimum stimulation threshold necessary for eliciting CCEPs, as well as a point at which the current-dependent relationship of several response metrics all saturate.

SIGNIFICANCE

This titration study provides practical, immediate guidance on optimal stimulation parameters to study specific features of CCEPs, which have been increasingly used to map both functional and epileptic brain networks in humans.

Effects of depth electrode montage and single-pulse electrical stimulation sites on neuronal responses and effective connectivity

OBJECTIVE

To determine the optimal depth electrode montages for the assessment of effective connectivity based on single-pulse electrical stimulation (SPES). To determine the effect of SPES locations on the extent of resulting neuronal propagations.

METHODS

We studied 14 epilepsy patients who underwent invasive monitoring with depth electrodes and measurement of cortico-cortical evoked potentials (CCEPs) and cortico-cortical spectral responses (CCSRs). We determined the effects of electrode montage and stimulus sites on the CCEP/CCSR amplitudes.

RESULTS

Bipolar and Laplacian montages effectively reduced the degree of SPES-related signal deflections at extra-cortical levels, including outside the brain, while maintaining those at the cortical level. SPES of structures more proximal to the deep white matter, compared to the cortical surface, elicited greater CCEPs and CCSRs.

CONCLUSIONS

On depth electrode recording, bipolar and Laplacian montages are suitable for measurement of near-field CCEPs and CCSRs. SPES of the white matter axons may induce neuronal propagations to extensive regions of the cerebral cortex.

SIGNIFICANCE

This study helps to establish the practical guidelines on the diagnostic use of CCEPs/CCSRs.