Evaluating dipolar source localization feasibility from intracerebral SEEG recordings

Antagonistic behavior of brain networks mediated by low-frequency oscillations: electrophysiological dynamics during internal-external attention switching

Antagonistic activity of brain networks likely plays a fundamental role in how the brain optimizes its performance by efficient allocation of computational resources. A prominent example involves externally/internally oriented attention tasks, implicating two anticorrelated, intrinsic brain networks: the default mode network (DMN) and the dorsal attention network (DAN). To elucidate electrophysiological underpinnings and causal interplay during attention switching, we recorded intracranial EEG (iEEG) from 25 epilepsy patients with electrode contacts localized in the DMN and DAN. We show antagonistic network dynamics of activation-related changes in high-frequency (> 50 Hz) and low-frequency (< 30 Hz) power. The temporal profile of information flow between the networks estimated by functional connectivity suggests that the activated network inhibits the other one, gating its activity by increasing the amplitude of the low-frequency oscillations. Insights about inter-network communication may have profound implications for various brain disorders in which these dynamics are compromised.

Predicting radiofrequency thermocoagulation surgical outcomes in refractory focal epilepsy patients using functional coupled neural mass model

Objective

The success rate of achieving seizure freedom after radiofrequency thermocoagulation surgery for patients with refractory focal epilepsy is about 20-40%. This study aims to enhance the prediction of surgical outcomes based on preoperative decisions through network model simulation, providing a reference for clinicians to validate and optimize surgical plans.

Methods

Twelve patients with epilepsy who underwent radiofrequency thermocoagulation were retrospectively reviewed in this study. A coupled model based on model subsets of the neural mass model was constructed by calculating partial directed coherence as the coupling matrix from stereoelectroencephalography (SEEG) signals. Multi-channel time-varying model parameters of excitation and inhibitions were identified by fitting the real SEEG signals with the coupled model. Further incorporating these model parameters, the coupled model virtually removed contacts destroyed in radiofrequency thermocoagulation or selected randomly. Subsequently, the coupled model after virtual surgery was simulated.

Results

The identified excitatory and inhibitory parameters showed significant difference before and after seizure onset (p < 0.05), and the trends of parameter changes aligned with the seizure process. Additionally, excitatory parameters of epileptogenic contacts were higher than that of non-epileptogenic contacts, and opposite findings were noticed for inhibitory parameters. The simulated signals of postoperative models to predict surgical outcomes yielded an area under the curve (AUC) of 83.33% and an accuracy of 91.67%.

Conclusion

The multi-channel coupled model proposed in this study with physiological characteristics showed a desirable performance for preoperatively predicting patients' prognoses.

Distributed source modeling of stereoencephalographic measurements of ictal activity

OBJECTIVES

Stereoelectroencephalography (SEEG) can define the epileptogenic zone (EZ). However, SEEG is susceptible to the sampling bias, where no SEEG recording is taken within a circumscribed EZ.

METHODS

Nine patients with medically refractory epilepsy underwent SEEG recording, and brain resection got positive outcomes. Ictal neuronal currents were estimated by distributed source modeling using the SEEG data and individual's anatomical magnetic resonance imaging. Using a retrospective leave-one-out data sub-sampling, we evaluated the sensitivity and specificity of the current estimates using MRI after surgical resection or radio-frequency ablation.

RESULTS

The sensitivity and specificity in detecting the EZ were indistinguishable from either the data from all electrodes or the sub-sampled data (rank sum test: rank sum = 23719, p = 0.13) when at least one remaining electrode contact was no more than 20 mm away.

CONCLUSIONS

The distributed neuronal current estimates of ictal SEEG data can mitigate the challenge of delineating the boundary of the EZ in cases of missing an electrode implanted within the EZ and a required second SEEG exploration.

SIGNIFICANCE

Distributed source modeling can be a tool for clinicians to infer the EZ by allowing for more flexible planning of the electrode implantation route and minimizing the number of electrodes.

Combining independent component analysis and source localization for improving spatial sampling of stereoelectroencephalography in epilepsy

Stereoelectroencephalography is a powerful intracerebral EEG recording method for the presurgical evaluation of epilepsy. It consists in implanting depth electrodes in the patient's brain to record electrical activity and map the epileptogenic zone, which should be resected to render the patient seizure-free. Stereoelectroencephalography has high spatial accuracy and signal-to-noise ratio but remains limited in the coverage of the explored brain regions. Thus, the implantation might provide a suboptimal sampling of epileptogenic regions. We investigate the potential of improving a suboptimal stereoelectroencephalography recording by performing source localization on stereoelectroencephalography signals. We propose combining independent component analysis, connectivity measures to identify components of interest, and distributed source modelling. This approach was tested on two patients with two implantations each, the first failing to characterize the epileptogenic zone and the second giving a better diagnosis. We demonstrate that ictal and interictal source localization performed on the first stereoelectroencephalography recordings matches the findings of the second stereo-EEG exploration. Our findings suggest that independent component analysis followed by source localization on the topographies of interest is a promising method for retrieving the epileptogenic zone in case of suboptimal implantation.

Visual Stimuli Modulate Local Field Potentials But Drive No High-Frequency Activity in Human Auditory Cortex

Neuroimaging studies suggest cross-sensory visual influences in human auditory cortices (ACs). Whether these influences reflect active visual processing in human ACs, which drives neuronal firing and concurrent broadband high-frequency activity (BHFA; >70 Hz), or whether they merely modulate sound processing is still debatable. Here, we presented auditory, visual, and audiovisual stimuli to 16 participants (7 women, 9 men) with stereo-EEG depth electrodes implanted near ACs for presurgical monitoring. Anatomically normalized group analyses were facilitated by inverse modeling of intracranial source currents. Analyses of intracranial event-related potentials (iERPs) suggested cross-sensory responses to visual stimuli in ACs, which lagged the earliest auditory responses by several tens of milliseconds. Visual stimuli also modulated the phase of intrinsic low-frequency oscillations and triggered 15-30 Hz event-related desynchronization in ACs. However, BHFA, a putative correlate of neuronal firing, was not significantly increased in ACs after visual stimuli, not even when they coincided with auditory stimuli. Intracranial recordings demonstrate cross-sensory modulations, but no indication of active visual processing in human ACs.

Optimization of patient-specific stereo-EEG recording sensitivity

Stereo-EEG is a minimally invasive technique used to localize the origin of epileptic activity (the epileptogenic zone) in patients with drug-resistant epilepsy. However, current stereo-EEG trajectory planning methods are agnostic to the spatial recording sensitivity of implanted electrodes. In this study, we used image-based patient-specific computational models to design optimized stereo-EEG electrode configurations. Patient-specific optimized electrode configurations exhibited substantially higher recording sensitivity than clinically implanted configurations, and this may lead to a more accurate delineation of the epileptogenic zone. The optimized configurations also achieved equally good or better recording sensitivity with fewer electrodes compared with clinically implanted configurations, and this may reduce the risk for complications, including intracranial haemorrhage. This approach improves localization of the epileptogenic zone by transforming the clinical use of stereo-EEG from a discrete ad hoc sampling to an intelligent mapping of the regions of interest.

Electrical Source Imaging of Somatosensory Evoked Potentials from Intracranial EEG Signals

Stereoelectroencephalography (SEEG) records electrical brain activity with intracerebral electrodes. However, it has an inherently limited spatial coverage. Electrical source imaging (ESI) infers the position of the neural generators from the recorded electric potentials, and thus, could overcome this spatial undersampling problem. Here, we aimed to quantify the accuracy of SEEG ESI under clinical conditions. We measured the somatosensory evoked potential (SEP) in SEEG and in high-density EEG (HD-EEG) in 20 epilepsy surgery patients. To localize the source of the SEP, we employed standardized low resolution brain electromagnetic tomography (sLORETA) and equivalent current dipole (ECD) algorithms. Both sLORETA and ECD converged to similar solutions. Reflecting the large differences in the SEEG implantations, the localization error also varied in a wide range from 0.4 to 10 cm. The SEEG ESI localization error was linearly correlated with the distance from the putative neural source to the most activated contact. We show that it is possible to obtain reliable source reconstructions from SEEG under realistic clinical conditions, provided that the high signal fidelity recording contacts are sufficiently close to the source of the brain activity.

Ambulatory Local Field Potential Recordings from the Thalamus in Epilepsy: A Feasibility Study

INTRODUCTION

Stimulation of the thalamus is gaining favor in the treatment of medically refractory multifocal and generalized epilepsy. Implanted brain stimulators capable of recording ambulatory local field potentials (LFPs) have recently been introduced, but there is little information to guide their use in thalamic stimulation for epilepsy. This study sought to assess the feasibility of chronically recording ambulatory interictal LFP from the thalamus in patients with epilepsy.

METHODS

In this pilot study, ambulatory LFP was recorded from patients who underwent sensing-enabled deep brain stimulation (DBS, 2 participants) or responsive neurostimulation (RNS, 3 participants) targeting the anterior nucleus of the thalamus (ANT, 2 electrodes), centromedian nucleus (CM, 7 electrodes), or medial pulvinar (PuM, 1 electrode) for multifocal or generalized epilepsy. Time-domain and frequency-domain LFP was investigated for epileptiform discharges, spectral peaks, circadian variation, and peri-ictal patterns.

RESULTS

Thalamic interictal discharges were visible on ambulatory recordings from both DBS and RNS. At-home interictal frequency-domain data could be extracted from both devices. Spectral peaks were noted at 10-15 Hz in CM, 6-11 Hz in ANT, and 19-24 Hz in PuM but varied in prominence and were not visible in all electrodes. In CM, 10-15 Hz power exhibited circadian variation and was attenuated by eye opening.

CONCLUSION

Chronic ambulatory recording of thalamic LFP is feasible. Common spectral peaks can be observed but vary between electrodes and across neural states. DBS and RNS devices provide a wealth of complementary data that have the potential to better inform thalamic stimulation for epilepsy.

Reconstruction and localization of auditory sources from intracerebral SEEG using independent component analysis

Stereo-electroencephalography (SEEG) is the surgical implantation of electrodes in the brain to better localize the epileptic network in pharmaco-resistant epileptic patients. This technique has exquisite spatial and temporal resolution. Still, the number and the position of the electrodes in the brain is limited and determined by the semiology and/or preliminary non-invasive examinations, leading to a large number of unexplored brain structures in each patient. Here, we propose a new approach to reconstruct the activity of non-sampled structures in SEEG, based on independent component analysis (ICA) and dipole source localization. We have tested this approach with an auditory stimulation dataset in ten patients. The activity directly recorded from the auditory cortex served as ground truth and was compared to the ICA applied on all non-auditory electrodes. Our results show that the activity from the auditory cortex can be reconstructed at the single trial level from contacts as far as ∼40 mm from the source. Importantly, this reconstructed activity is localized via dipole fitting in the proximity of the original source. In addition, we show that the size of the confidence interval of the dipole fitting is a good indicator of the reliability of the result, which depends on the geometry of the SEEG implantation. Overall, our approach allows reconstructing the activity of structures far from the electrode locations, partially overcoming the spatial sampling limitation of intracerebral recordings.

Ideal current dipoles are appropriate source representations for simulating neurons for intracranial recordings

OBJECTIVE

To determine whether dipoles are an appropriate simplified representation of neural sources for stereo-EEG (sEEG).

METHODS

We compared the distributions of voltages generated by a dipole, biophysically realistic cortical neuron models, and extended regions of cortex to determine how well a dipole represented neural sources at different spatial scales and at electrode to neuron distances relevant for sEEG. We also quantified errors introduced by the dipole approximation of neural sources in sEEG source localization using standardized low-resolution electrotomography (sLORETA).

RESULTS

For pyramidal neurons, the coefficient of correlation between voltages generated by a dipole and neuron model were > 0.9 for distances > 1 mm. For small regions of cortex (∼0.1 cm2), the error in voltages between a dipole and region was < 100 µV for all distances. However, larger regions of active cortex (>5 cm2) yielded > 50 µV errors within 1.5 cm of an electrode when compared to single dipoles. Finally, source localization errors were < 5 mm when using dipoles to represent realistic neural sources.

CONCLUSIONS

Single dipoles are an appropriate source model to represent both single neurons and small regions of active cortex, while multiple dipoles are required to represent large regions of cortex.

SIGNIFICANCE

Dipoles are computationally tractable and valid source models for sEEG.

Source localization of ictal SEEG to predict postoperative seizure outcome

OBJECTIVE

Stereo-electroencephalography (SEEG) is inherently-three-dimensional and can be modeled using source localization. This study aimed to assess the validity of ictal SEEG source localization.

METHODS

The dominant frequency at ictal onset was used for source localization in the time and frequency domains using rotating dipoles and current density maps. Validity was assessed by concordance with the epileptologist-defined seizure onset zone (conventional SOZ) and the surgical treatment volume (TV) of seizure-free versus non-seizure-free patients.

RESULTS

Source localization was performed on 68 seizures from 27 patients. Median distance to nearest contact in the conventional SOZ was 7 (IQR 6-12) mm for time-domain dipoles. Current density predicted ictal activity with up to 86 % (60-87 %) accuracy. Distance from time-domain dipoles to the TV was smaller (P = 0.045) in seizure-free (2 [0-4] mm) versus non-seizure-free (12 [2-17] mm) patients, and predicted surgical outcome with 91 % sensitivity and 63 % specificity. Removing near-field data from contacts within the TV negated outcome prediction (P = 0.51).

CONCLUSIONS

Source localization of SEEG accurately mapped ictal onset compared with conventional interpretation. Proximity of dipoles to the TV predicted seizure outcome when near-field recordings were analyzed.

SIGNIFICANCE

Ictal SEEG source localization is useful in corroborating the epileptogenic zone, assuming near-field recordings are obtained.

Inter-individual variability in dorsal stream dynamics during word production

The current standard model of language production involves a sensorimotor dorsal stream connecting areas in the temporo-parietal junction with those in the inferior frontal gyrus and lateral premotor cortex. These regions have been linked to various aspects of word production such as phonological processing or articulatory programming, primarily through neuropsychological and functional imaging group studies. Most if not all the theoretical descriptions of this model imply that the same network should be identifiable across individual speakers. We tested this hypothesis by quantifying the variability of activation observed across individuals within each dorsal stream anatomical region. This estimate was based on electrical activity recorded directly from the cerebral cortex with millisecond accuracy in awake epileptic patients clinically implanted with intracerebral depth electrodes for pre-surgical diagnosis. Each region's activity was quantified using two different metrics-intra-cerebral evoked related potentials and high gamma activity-at the level of the group, the individual and the recording contact. The two metrics show simultaneous activation of parietal and frontal regions during a picture naming task, in line with models that posit interactive processing during word retrieval. They also reveal different levels of between-patient variability across brain regions, except in core auditory and motor regions. The independence and non-uniformity of cortical activity estimated through the two metrics push the current model towards sub-second and sub-region explorations focused on individualized language speech production. Several hypotheses are considered for this within-region heterogeneity.

On the localization of oscillatory sources from (S)EEG recordings

Electrophysiological brain source localization consists in estimating the source positions and activities responsible for the (S)EEG measurements. The localization procedure is usually carried out in the time domain, however in specific situations the activities of interest can be located at well defined frequencies, e.g. in response to a rhythmic stimulation. This paper addresses the problem of sparse localization of multiple sources oscillating at the same frequency. In particular the non-unicity of the solution is emphasized, as alternative source maps involving equivalent or less number of sources can be found, challenging source localization methods based on sparsity. These limitations are illustrated under a realistic SEEG simulation framework, and the usefulness to perform localization for this modality is strengthen out.

Informed MEG/EEG source imaging reveals the locations of interictal spikes missed by SEEG

Determining the accurate locations of interictal spikes has been fundamental in the presurgical evaluation of epilepsy surgery. Stereo-electroencephalography (SEEG) is able to directly record cortical activity and localize interictal spikes. However, the main caveat of SEEG techniques is that they have limited spatial sampling (covering <5% of the whole brain), which may lead to missed spikes originating from brain regions that were not covered by SEEG. To address this problem, we propose a SEEG-informed minimum-norm estimates (SIMNE) method by combining SEEG with magnetoencephalography (MEG) or EEG. Specifically, the spike locations determined by SEEG offer as a priori information to guide MEG source reconstruction. Both computer simulations and experiments using data from five epilepsy patients were conducted to evaluate the performance of SIMNE. Our results demonstrate that SIMNE generates more accurate source estimation than a traditional minimum-norm estimates method and reveals the locations of spikes missed by SEEG, which would improve presurgical evaluation of the epileptogenic zone.

SEEG in 3D: Interictal Source Localization From Intracerebral Recordings

BACKGROUND

Stereo-electroencephalography (SEEG) uses a three-dimensional configuration of depth electrodes to localize epileptiform activity, but traditional analysis of SEEG is spatially restricted to the point locations of the electrode contacts. Interpolation of brain activity between contacts might allow for three-dimensional representation of epileptiform activity and avoid pitfalls of SEEG interpretation.

OBJECTIVE

The goal of this study was to validate SEEG-based interictal source localization and assess the ability of this technique to monitor far-field activity in non-implanted brain regions.

METHODS

Interictal epileptiform discharges were identified on SEEG in 26 patients who underwent resection, ablation, or disconnection of the suspected epileptogenic zone. Dipoles without (free) and with (scan) gray matter restriction, and current density (sLORETA and SWARM methods), were calculated using a finite element head model. Source localization results were compared to the conventional irritative zone (IZ) and the surgical treatment volumes (TV) of seizure-free vs. non-seizure-free patients.

RESULTS

The median distance from dipole solutions to the nearest contact in the conventional IZ was 7 mm (interquartile range 4-15 mm for free dipoles and 4-14 mm for scan dipoles). The IZ modeled with SWARM predicted contacts within the conventional IZ with 83% (75-100%) sensitivity and 94% (88-100%) specificity. The proportion of current within the TV was greater in seizure-free patients (P = 0.04) and predicted surgical outcome with 45% sensitivity and 93% specificity. Dipole solutions and sLORETA results did not correlate with seizure outcome. Addition of scalp EEG led to more superficial modeled sources (P = 0.03) and negated the ability to predict seizure outcome (P = 0.23). Removal of near-field data from contacts within the TV resulted in smearing of the current distribution (P = 0.007) and precluded prediction of seizure freedom (P = 0.20).

CONCLUSIONS

Source localization accurately represented interictal discharges from SEEG. The proportion of current within the TV distinguished between seizure-free and non-seizure-free patients when near-field recordings were obtained from the surgical target. The high prevalence of deep sources in this cohort likely obscured any benefit of concurrent scalp EEG. SEEG-based interictal source localization is useful in illustrating and corroborating the epileptogenic zone. Additional techniques are needed to localize far-field epileptiform activity from non-implanted brain regions.

Speech-related auditory salience detection in the posterior superior temporal region

Processing auditory human speech requires both detection (early and transient) and analysis (sustained). We analyzed high gamma (70-110 Hz) activity of intracranial electroencephalography waveforms acquired during an auditory task that paired forward speech, reverse speech, and signal correlated noise. We identified widespread superior temporal sites with sustained activity responding only to forward and reverse speech regardless of paired order. More localized superior temporal auditory onset sites responded to all stimulus types when presented first in a pair and responded in recurrent fashion to the second paired stimulus in select conditions even in the absence of interstimulus silence; a novel finding. Auditory onset activity to a second paired sound recurred according to relative salience, with evidence of partial suppression during linguistic processing. We propose that temporal lobe auditory onset sites facilitate a salience detector function with hysteresis of 200 ms and are influenced by cortico-cortical feedback loops involving linguistic processing and articulation.

Data-driven method to infer the seizure propagation patterns in an epileptic brain from intracranial electroencephalography

Surgical interventions in epileptic patients aimed at the removal of the epileptogenic zone have success rates at only 60-70%. This failure can be partly attributed to the insufficient spatial sampling by the implanted intracranial electrodes during the clinical evaluation, leading to an incomplete picture of spatio-temporal seizure organization in the regions that are not directly observed. Utilizing the partial observations of the seizure spreading through the brain network, complemented by the assumption that the epileptic seizures spread along the structural connections, we infer if and when are the unobserved regions recruited in the seizure. To this end we introduce a data-driven model of seizure recruitment and propagation across a weighted network, which we invert using the Bayesian inference framework. Using a leave-one-out cross-validation scheme on a cohort of 45 patients we demonstrate that the method can improve the predictions of the states of the unobserved regions compared to an empirical estimate that does not use the structural information, yet it is on the same level as the estimate that takes the structure into account. Furthermore, a comparison with the performed surgical resection and the surgery outcome indicates a link between the inferred excitable regions and the actual epileptogenic zone. The results emphasize the importance of the structural connectome in the large-scale spatio-temporal organization of epileptic seizures and introduce a novel way to integrate the patient-specific connectome and intracranial seizure recordings in a whole-brain computational model of seizure spread.

Distributed source modeling of intracranial stereoelectro-encephalographic measurements

Intracranial stereoelectroencephalography (sEEG) provides unsurpassed sensitivity and specificity for human neurophysiology. However, functional mapping of brain functions has been limited because the implantations have sparse coverage and differ greatly across individuals. Here, we developed a distributed, anatomically realistic sEEG source-modeling approach for within- and between-subject analyses. In addition to intracranial event-related potentials (iERP), we estimated the sources of high broadband gamma activity (HBBG), a putative correlate of local neural firing. Our novel approach accounted for a significant portion of the variance of the sEEG measurements in leave-one-out cross-validation. After logarithmic transformations, the sensitivity and signal-to-noise ratio were linearly inversely related to the minimal distance between the brain location and electrode contacts (slope≈-3.6). The signa-to-noise ratio and sensitivity in the thalamus and brain stem were comparable to those locations at the vicinity of electrode contact implantation. The HGGB source estimates were remarkably consistent with analyses of intracranial-contact data. In conclusion, distributed sEEG source modeling provides a powerful neuroimaging tool, which facilitates anatomically-normalized functional mapping of human brain using both iERP and HBBG data.

A P300-Based BCI System Using Stereoelectroencephalography and Its Application in a Brain Mechanistic Study

Stereoelectroencephalography (SEEG) signals can be obtained by implanting deep intracranial electrodes. SEEG depth electrodes can record brain activity from the shallow cortical layer and deep brain structures, which is not achievable through other recording techniques. Moreover, SEEG has the advantage of a high signal-to-noise ratio (SNR). Therefore, it provides a potential way to establish a highly efficient brain-computer interface (BCI) and aid in understanding human brain activity. In this study, we implemented a P300-based BCI using SEEG signals. A single-character oddball paradigm was applied to elicit P300. To predict target characters, we fed the feature vectors extracted from the signals collected by five SEEG contacts into a Bayesian linear discriminant analysis (BLDA) classifier. Thirteen epileptic patients implanted with SEEG electrodes participated in the experiment and achieved an average online spelling accuracy of 93.85%. Moreover, through single-contact decoding analysis and simulated online analysis, we found that the SEEG-based BCI system achieved a high performance even when using a single signal channel. Furthermore, contacts with high decoding accuracies were mainly distributed in the visual ventral pathway, especially the fusiform gyrus (FG) and lingual gyrus (LG), which played an important role in building P300-based SEEG BCIs. These results might provide new insights into P300 mechanistic studies and the corresponding BCIs.

Optimization of the Cortical Traveling Wave Analysis framework for feasibility in Stereo-Electroencephalography

The study of brain waves propagation is of interest to understand the neural involvement in both physiological and pathological events, such as interictal epileptic spikes (IES). The possibility to track the trajectory of IESs could be useful to better characterize the role of the involved structures in the epileptic network, adding valuable information to the epileptic focus localization. Methods for the cortical traveling wave analysis (CTWA) have been proposed to trace the preferred propagation path of sleep slow waves, using scalp high-density EEG and reconstructing the trajectories both in the sensors and in the sources space. In this work, we propose a feasibility study of the application of these concepts to Stereo-EEG (SEEG) data for the analysis of IES. Through simulations, we selected the best performing Electrical Source Imaging inverse solution for our purpose and illustrate the CTWA procedure. We further show an exemplary application on real data and discuss advantages and pitfalls of the application of CTWA in SEEG.

A P300-based Brain Computer Interface Using Stereo-electroencephalography Signals

Stereo-electroencephalography (SEEG) signals can be obtained by implanting deep intracranial electrodes, which are currently used for epileptic diagnosis. In this study, we implemented a P300-based Brain Computer Interface (BCI) using SEEG signals. 40 buttons corresponding to 40 numbers displayed in a graphical user interface (GUI) were intensified in a random order. To select a number, the user could focus on the corresponding button when it was flashing. Five epileptic patients implanted with SEEG electrodes attended the experiment and achieved an average online accuracy of 97.33%. Moreover, through single contact decoding and simulated online analysis, we found that these subjects achieved an average accuracy of 82.00% using a single channel of signal. These results indicated that our SEEG-based BCI had a high performance, which was mainly because of the high quality of SEEG signals.

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

BACKGROUND

Cortico-Cortical Evoked Potentials (CCEPs) are a novel low frequency stimulation method used for brain mapping during intracranial epilepsy investigations. Only a handful of metrics have been applied to CCEP data to infer connectivity, and no comparison as to which is best has been performed.

NEW METHOD

We implement a novel method which involved superimposing synthetic cortical responses onto stereoelectroencephalographic (SEEG) data, and use this to compare several metric's ability to detect the simulated patterns. In this we compare two commonly employed metrics currently used in CCEP analysis against eight time series similarity metrics (TSSMs), which have been widely used in machine learning and pattern matching applications.

RESULTS

Root Mean Square (RMS), a metric commonly employed in CCEP analysis, was sensitive to a wide variety of response patterns, but insensitive to simulated epileptiform patterns. Autoregressive (AR) coefficients calculated by Burg's method were also sensitive to a wide range of patterns, but were extremely sensitive to epileptiform patterns. Other metrics which employed elastic warping techniques were less sensitive to the simulated response patterns.

COMPARISON WITH EXISTING METHODS

Our study is the first to compare CCEP connectivity metrics against one-another. Our results found that RMS, which has been used in many CCEP studies previously, was the most sensitive metric across a wide range of patterns.

CONCLUSIONS

Our novel method showed that RMS is a robust and sensitive measure, validating much of the findings of the SEEG-CCEP literature to date. Autoregressive coefficients may also be a useful metric to investigate epileptic networks.

Intracranial EEG recording in autoimmune epilepsy dataset

The appearance of seizures in autoimmune epilepsy on intracranial recordings has not been previously demonstrated. The following data shows a multifocal epilepsy in a patient with seronegative autoimmune epilepsy (reported here; “Electroclinical Insights into Autoimmune Epilepsy”, Gillinder, 2019). Independent seizures were seen to arise from 5 separate foci. These all began with slow repetitive spiking in a highly restricted area. Only after many minutes would this activity spread to other regions. Despite arising from different locations, all foci affected the posterior insula resulting in clinical symptoms.

Automatic detection of the epileptogenic zone: An application of the fingerprint of epilepsy

BACKGROUND

The successful delineation of the epileptogenic zone in epilepsy monitoring is crucial for achieving seizure freedom after epilepsy surgery.

NEW METHOD

We aim to improve epileptogenic zone localization by utilizing a computer-assisted tool for the automated grading of the seizure activity recorded in various locations for 20 patients undergoing stereo electroencephalography. Their epileptic seizures were processed to extract two potential biomarkers. The concentration of these biomarkers from within each patient's implantation were then graded to identify their epileptogenic zone and were compared to the clinical assessment.

RESULTS

Our technique was capable of ranking the clinically defined epileptogenic zone with high accuracy, above 95%, with a true to false positive ratio of 1:1.52, and was effective with both temporal and extra-temporal onset epilepsies.

COMPARISON WITH EXISTING METHOD

We compared our method to two other groups performing localization using similar biomarkers. Our classification metrics, sensitivity and precision together were comparable to both groups and our overall accuracy from a larger population was also higher then both.

CONCLUSIONS

Our method is highly accurate, automated and non-parametric providing clinicians another tool that can be used to help identify the epileptogenic zone in patients undergoing the stereo electroencephalography procedure for epilepsy monitoring.

Utility of stereo-electroencephalography recording guided by magnetoencephalography in the surgical treatment of epilepsy patients with negative magnetic resonance imaging results

Objective: It is challenging for neurosurgeons to perform surgeries on patients without detectable structural lesions. Therefore, this retrospective study aimed to explore the outcome of stereo-electroencephalography (SEEG) in suspicious areas guided by magnetoencephalography (MEG)-magnetic resonance imaging (MRI) reconstruction in MRI-negative epilepsy patients. Methods: This study included 47 patients with negative-MRI epilepsy. Seizure outcome at 24 months was assessed using a modified Engel's classification. Accordingly, class I and II were considered favorable outcomes, whereas classes III and IV were unfavorable. Furthermore, patients were classified into a consistent group if the results of MEG and SEEG indicated the same area of the brain. The relationship between surgical outcome and the concordance of MEG and SEEG was analyzed. Results: A complete seizure-free condition was achieved in 22 (47%) patients. Sex, handedness, age and duration of illness were not significantly associated with seizure-free outcome (p = .187 [Pearson chi-squared test]). The number of patients with favorable outcome (Engle I and II) was as high as 68% at the time of follow-up. Furthermore, more seizure-free patients were found in the SEEG and MEG consistent group. Conclusions: SEEG is a valuable tool in the pre-evaluation for resective epilepsy surgery, particularly in negative-MRI epilepsy patients; MEG greatly facilitates localization for SEEG electrode implantation. However, none of these tools are absolutely sensitive and reliable; therefore, collecting as much information as possible is necessary to achieve satisfactory results in epilepsy surgery.

Characterizing the role of the structural connectome in seizure dynamics

How does the human brain's structural scaffold give rise to its intricate functional dynamics? This is a central question in translational neuroscience that is particularly relevant to epilepsy, a disorder affecting over 50 million subjects worldwide. Treatment for medication-resistant focal epilepsy is often structural-through surgery or laser ablation-but structural targets, particularly in patients without clear lesions, are largely based on functional mapping via intracranial EEG. Unfortunately, the relationship between structural and functional connectivity in the seizing brain is poorly understood. In this study, we quantify structure-function coupling, specifically between white matter connections and intracranial EEG, across pre-ictal and ictal periods in 45 seizures from nine patients with unilateral drug-resistant focal epilepsy. We use high angular resolution diffusion imaging (HARDI) tractography to construct structural connectivity networks and correlate these networks with time-varying broadband and frequency-specific functional networks derived from coregistered intracranial EEG. Across all frequency bands, we find significant increases in structure-function coupling from pre-ictal to ictal periods. We demonstrate that short-range structural connections are primarily responsible for this increase in coupling. Finally, we find that spatiotemporal patterns of structure-function coupling are highly stereotyped for each patient. These results suggest that seizures harness the underlying structural connectome as they propagate. Mapping the relationship between structural and functional connectivity in epilepsy may inform new therapies to halt seizure spread, and pave the way for targeted patient-specific interventions.

Electrophysiological Brain Connectivity: Theory and Implementation

We review the theory and algorithms of electrophysiological brain connectivity analysis. This tutorial is aimed at providing an introduction to brain functional connectivity from electrophysiological signals, including electroencephalography (EEG), magnetoencephalography (MEG), electrocorticography (ECoG), stereoelectroencephalography (SEEG). Various connectivity estimators are discussed, and algorithms introduced. Important issues for estimating and mapping brain functional connectivity with electrophysiology are discussed.

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.

Electromagnetic source imaging using simultaneous scalp EEG and intracranial EEG: An emerging tool for interacting with pathological brain networks

OBJECTIVE

The goal of this study is to investigate the performance, merits and limitations of source imaging using intracranial EEG (iEEG) recordings and to compare its accuracy to the results of EEG source imaging. Accuracy in this study, is measured both by determining the location and inter-nodal connectivity of underlying brain networks.

METHODS

Systematic computer simulation studies are conducted to evaluate iEEG-based source imaging vs. EEG-based source imaging, and source imaging using both EEG and iEEG. To test the source imaging models, networks of inter-connected nodes (in terms of activity) are simulated. The location of the network nodes is randomly selected within a realistic geometry head model and a connectivity link is created among these nodes based on a multi-variate auto-regressive (MVAR) model. Then the forward problem is solved to calculate the potentials at the electrodes and noise (white and correlated) is added to these simulated potentials to simulate realistic measurements. Subsequently, the inverse problem is solved and an algorithm based on principle component analysis is performed on the estimated source activities to determine the location of the simulated network nodes. The activity of these nodes (over time), is then extracted, and used to estimate the connectivity links among the mentioned nodes using Granger causality analysis.

RESULTS

Source imaging based on iEEG recordings may or may not improve the accuracy in localization, depending on the number and location of active nodes relative to iEEG electrodes and to other nodes within the network. However, our simulation results suggest that combining EEG and iEEG modalities (simultaneous scalp and intracranial recordings) can improve the imaging accuracy significantly.

CONCLUSIONS

While iEEG source imaging is useful in estimating the exact location of sources near the iEEG electrodes, combining EEG and iEEG recordings can achieve a more accurate imaging due to the high spatial coverage of the scalp electrodes and the added near field information provided by the iEEG electrodes.

SIGNIFICANCE

The present results suggest the feasibility of localizing brain electrical sources from iEEG recordings and improving EEG source localization using simultaneous EEG and iEEG recordings to cover the whole brain. The hybrid EEG and iEEG source imaging can assist the clinicians when unequivocal decisions about determining the epileptogenic zone cannot be reached using a single modality.

SEEG dipole source localization based on an empirical Bayesian approach taking into account forward model uncertainties

Electromagnetic brain source localization consists in the inversion of a forward model based on a limited number of potential measurements. A wide range of methods has been developed to regularize this severely ill-posed problem and to reduce the solution space, imposing spatial smoothness, anatomical constraint or sparsity of the activated source map. This last criteria, based on physiological assumptions stating that in some particular events (e.g., epileptic spikes, evoked potential) few focal area of the brain are simultaneously actives, has gained more and more interest. Bayesian approaches have the ability to provide sparse solutions under adequate parametrization, and bring a convenient framework for the introduction of priors in the form of probabilistic density functions. However the quality of the forward model is rarely questioned while this parameter has undoubtedly a great influence on the solution. Its construction suffers from numerous approximation and uncertainties, even when using realistic numerical models. In addition, it often encodes a coarse sampling of the continuous solution space due to the computational burden its inversion implies. In this work we propose an empirical Bayesian approach to take into account the uncertainties of the forward model by allowing constrained variations around a prior physical model, in the particular context of SEEG measurements. We demonstrate on simulations that the method enhance the accuracy of the source time-course estimation as well as the sparsity of the resulting source map. Results on real signals prove the applicability of the method in real contexts.

Assessing human brain impedance using simultaneous surface and intracerebral recordings

Most of the literature on the brain impedance proposes a frequency-independent resistive model. Recently, this conclusion was tackled by a series of papers (Bédard et al., 2006; Bédard and Destexhe, 2009; Gomes et al., 2016), based on microscopic sale modeling and measurements. Our paper aims to investigate the impedance issue using simultaneous in vivo depth and surface signals recorded during intracerebral electrical stimulation of epileptic patients, involving a priori different tissues with different impedances. Our results confirm the conclusions from Logothethis et al. (2007): there is no evidence of frequency dependence of the brain tissue impedance (more precisely, there is no difference, in terms of frequency filtering, between the brain and the skull bone), at least at a macroscopic scale. In order to conciliate findings from both microscopic and macroscopic scales, we recall different neural/synaptic current generators' models from the literature and we propose an original computational model, based on fractional dynamics.

Combining bayesian source imaging with equivalent dipole approach to solve the intracranial EEG source localization problem

The brain source localization problem has been extensively studied in the past years, yielding a large panel of methodologies, each bringing their own strengths and weaknesses. Combining several of these approaches might help in enhancing their respective performance. Our study is carried out in the particular context of intracranial recordings, with the objective to explain the measurements based on a reduced number of dipolar activities. We take benefit of the sparse nature of the Bayesian approaches to separate the noise from the source space, and to distinguish between several source contributions on the electrodes. This first step provides accurate estimates of the dipole projections, which can be used as an entry to an equivalent current dipole fitting procedure. We demonstrate on simulations that the localization results are significantly enhanced by this post-processing step when up to five dipoles are activated simultaneously.

Data integration: Combined imaging and electrophysiology data in the cloud

There has been an increasing effort to correlate electrophysiology data with imaging in patients with refractory epilepsy over recent years. IEEG.org provides a free-access, rapidly growing archive of imaging data combined with electrophysiology data and patient metadata. It currently contains over 1200 human and animal datasets, with multiple data modalities associated with each dataset (neuroimaging, EEG, EKG, de-identified clinical and experimental data, etc.). The platform is developed around the concept that scientific data sharing requires a flexible platform that allows sharing of data from multiple file formats. IEEG.org provides high- and low-level access to the data in addition to providing an environment in which domain experts can find, visualize, and analyze data in an intuitive manner. Here, we present a summary of the current infrastructure of the platform, available datasets and goals for the near future.

Influence of the stereo-EEG sensors setup and of the averaging on the dipole localization problem

While scalp EEG/MEG source imaging have been extensively studied in the last two decades, the case of source localization from invasive measurements has resulted in few works to date. Yet there is a lot to gain from stereo-electroencephalographic (SEEG) recordings, providing high signal to noise ratio measurements of the explored brain structures. The SEEG setup consists in multi-contact electrodes inserted in the brain volume, each containing a dozen of collinear measuring contacts. This particular setup raises the question of the conditioning of the inverse problem. In recent works, we have evaluated the feasibility to localize a single dominant equivalent dipole facing different sensors and noise configurations. We deepen here the analysis by evaluating the influence of the chosen subset of sensors and of the number of averaged time samples on the accuracy of the localization. We conduct experiments on simulated data as well as on real epileptic spikes, illustrating the trade off to be made between these two factors.