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

Simultaneous invasive and non-invasive recordings in humans: A novel Rosetta stone for deciphering brain activity

Simultaneous noninvasive and invasive electrophysiological recordings provide a unique opportunity to achieve a comprehensive understanding of human brain activity, much like a Rosetta stone for human neuroscience. In this review we focus on the increasingly-used powerful combination of intracranial electroencephalography (iEEG) with scalp electroencephalography (EEG) or magnetoencephalography (MEG). We first provide practical insight on how to achieve these technically challenging recordings. We then provide examples from clinical research on how simultaneous recordings are advancing our understanding of epilepsy. This is followed by the illustration of how human neuroscience and methodological advances could benefit from these simultaneous recordings. We conclude with a call for open data sharing and collaboration, while ensuring neuroethical approaches and argue that only with a true collaborative approach the promises of simultaneous recordings will be fulfilled.

Effect of Inverse Solutions, Connectivity Measures, and Node Sizes on EEG Source Network: A Simultaneous EEG Study

Brain network provides an essential perspective for studying normal and pathological brain activities. Reconstructing the brain network in the source space becomes more needed, for example, as a target in non-invasive neuromodulation. Precise estimating source activities from the scalp EEG is still challenging because it is an ill-posed question and because of the volume conduction effect. There is no consensus on how to reconstruct the EEG source network. This study uses simultaneous scalp EEG and stereo-EEG to investigate the effect of inverse solutions, connectivity measures, and node sizes on the reconstruction of the source network. We evaluated the performance of different methods on both source activity and network. Numerical simulation was also carried out for comparison. The weighted phase-lag index (wPLI) method achieved significantly better performance on the reconstructed networks in source space than five other connectivity measures (directed transfer function (DTF), partial directed coherence (PDC), efficient effective connectivity (EEC), Pearson correlation coefficient (PCC), and amplitude envelope correlation (AEC)). There is no significant difference between the inverse solutions (standardized low-resolution brain electromagnetic tomography (sLORETA), weighted minimum norm estimate (wMNE), and linearly constrained minimum variance (LCMV) beamforming) on the reconstructed source networks. The source network based on signal phases can fit intracranial activities better than signal waveform properties or causality. Our study provides a basis for reconstructing source space networks from scalp EEG, especially for future neuromodulation research.

Spontaneous HFO Sequences Reveal Propagation Pathways for Precise Delineation of Epileptogenic Networks

Epilepsy, a neurological disorder affecting millions worldwide, poses great challenges in precisely delineating the epileptogenic zone - the brain region generating seizures - for effective treatment. High-frequency oscillations (HFOs) are emerging as promising biomarkers; however, the clinical utility is hindered by the difficulties in distinguishing pathological HFOs from non- epileptiform activities at single electrode and single patient resolution and understanding their dynamic role in epileptic networks. Here, we introduce an HFO-sequencing approach to analyze spontaneous HFOs traversing cortical regions in 40 drug-resistant epilepsy patients. This data- driven method automatically detected over 8.9 million HFOs, pinpointing pathological HFO- networks, and unveiled intricate millisecond-scale spatiotemporal dynamics, stability, and functional connectivity of HFOs in prolonged intracranial EEG recordings. These HFO sequences demonstrated a significant improvement in localization of epileptic tissue, with an 818.47% increase in concordance with seizure-onset zone (mean error: 2.92 mm), compared to conventional benchmarks. They also accurately predicted seizure outcomes for 90% AUC based on pre-surgical information using generalized linear models. Importantly, this mapping remained reliable even with short recordings (mean standard deviation: 3.23 mm for 30-minute segments). Furthermore, HFO sequences exhibited distinct yet highly repetitive spatiotemporal patterns, characterized by pronounced synchrony and predominant inward information flow from periphery towards areas involved in propagation, suggesting a crucial role for excitation-inhibition balance in HFO initiation and progression. Together, these findings shed light on the intricate organization of epileptic network and highlight the potential of HFO-sequencing as a translational tool for improved diagnosis, surgical targeting, and ultimately, better outcomes for vulnerable patients with drug-resistant epilepsy.

One Sentence Summary

Pathological fast brain oscillations travel like traffic along varied routes, outlining recurrently visited neural sites emerging as critical hotspots in epilepsy network.

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.

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.

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.

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.

Localizing targets for neuromodulation in drug-resistant epilepsy using intracranial EEG and computational model

Neuromodulation has emerged as a promising technique for the treatment of epilepsy. The target for neuromodulation is critical for the effectiveness of seizure control. About 30% of patients with drug-resistant epilepsy (DRE) fail to achieve seizure freedom after surgical intervention. It is difficult to find effective brain targets for neuromodulation in these patients because brain regions are damaged during surgery. In this study, we propose a novel approach for localizing neuromodulatory targets, which uses intracranial EEG and multi-unit computational models to simulate the dynamic behavior of epileptic networks through external stimulation. First, we validate our method on a multivariate autoregressive model and compare nine different methods of constructing brain networks. Our results show that the directed transfer function with surrogate analysis achieves the best performance. Intracranial EEGs of 11 DRE patients are further analyzed. These patients all underwent surgery. In three seizure-free patients, the localized targets are concordant with the resected regions. For the eight patients without seizure-free outcome, the localized targets in three of them are outside the resected regions. Finally, we provide candidate targets for neuromodulation in these patients without seizure-free outcome based on virtual resected epileptic network. We demonstrate the ability of our approach to locate optimal targets for neuromodulation. We hope that our approach can provide a new tool for localizing patient-specific targets for neuromodulation therapy in DRE.

EEG Source Imaging using GANs with Deep Image Prior

Brain source localization from electroencephalogram (EEG) signals is an challenging problem for noninvasively localizing the brain activity. Conventional methods use handcrafted regularization terms based on neural-physiological assumptions by exploiting the spatial-temporal structure on the source signals. In recent years, deep learning frameworks have demonstrated superior performance for solving the inverse problems in the natural and medical imaging field. This study proposes a novel unsupervised learning training-free framework based on Generative Adversarial Networks and deep image prior (GANs-DIP) as a generative model simulating spatially structured source signal. The proposed framework can faithfully recover extended source patches activation patterns of the brain in an unsupervised manner. Numerical experiments on a realistic brain model are performed under different levels of signal-to-noise ratio (SNR). The proposed model shows satisfactory performance in recovering the underlying source activation.

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.

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.

Probabilistic Structure Learning for EEG/MEG Source Imaging With Hierarchical Graph Priors

Brain source imaging is an important method for noninvasively characterizing brain activity using Electroencephalogram (EEG) or Magnetoencephalography (MEG) recordings. Traditional EEG/MEG Source Imaging (ESI) methods usually assume the source activities at different time points are unrelated, and do not utilize the temporal structure in the source activation, making the ESI analysis sensitive to noise. Some methods may encourage very similar activation patterns across the entire time course and may be incapable of accounting the variation along the time course. To effectively deal with noise while maintaining flexibility and continuity among brain activation patterns, we propose a novel probabilistic ESI model based on a hierarchical graph prior. Under our method, a spanning tree constraint ensures that activity patterns have spatiotemporal continuity. An efficient algorithm based on an alternating convex search is presented to solve the resulting problem of the proposed model with guaranteed convergence. Comprehensive numerical studies using synthetic data on a realistic brain model are conducted under different levels of signal-to-noise ratio (SNR) from both sensor and source spaces. We also examine the EEG/MEG datasets in two real applications, in which our ESI reconstructions are neurologically plausible. All the results demonstrate significant improvements of the proposed method over benchmark methods in terms of source localization performance, especially at high noise levels.

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.

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.

Electromagnetic Brain Source Imaging by Means of a Robust Minimum Variance Beamformer

OBJECTIVE

Adaptive beamformer methods that have been extensively used for functional brain imaging using EEG/MEG (magnetoencephalography) signals are sensitive to model mismatches. We propose a robust minimum variance beamformer (RMVB) technique, which explicitly incorporates the uncertainty of the lead field matrix into the estimation of spatial-filter weights that are subsequently used to perform the imaging.

METHODS

The uncertainty of the lead field is modeled by ellipsoids in the RMVB method; these hyperellipsoids (ellipsoids in higher dimensions) define regions of uncertainty for a given nominal lead field vector. These ellipsoids are estimated empirically by sampling lead field vectors surrounding each point of the source space, or more generally by building several forward models for the source space. Once these uncertainty regions (ellipsoids) are estimated, they are used to perform the source-imaging task. Computer simulations are conducted to evaluate the performance of the proposed RMVB technique.

RESULTS

Our results show that robust beamformers can outperform conventional beamformers in terms of localization error, recovering source dynamics, and estimation of the underlying source extents when uncertainty in the lead field matrix is properly determined and modeled.

CONCLUSION

The RMVB can be substituted for conventional beamformers, especially in applications where source imaging is performed off-line, and computational speed and complexity are not of major concern.

SIGNIFICANCE

A high-quality source imaging can be utilized in various applications, such as determining the epileptogenic zone in medically intractable epilepsy patients or estimating the time course of activity, which is a required step for computing the functional connectivity of brain networks.