Localizing hidden Interictal Epileptiform Discharges with simultaneous intracerebral and scalp high-density EEG recordings.

BACKGROUND: Scalp EEG is one of the main tools in the clinical evaluation of epilepsy. In some cases intracranial Interictal Epileptiform Discharges (IEDs) are not visible from the scalp. Recent studies have shown the feasibility of revealing them in the EEG if their timings are extracted from simultaneous intracranial recordings, but their potential for the localization of the epileptogenic zone is not yet well defined. NEW METHOD: We recorded simultaneous high-density EEG (HD-EEG) and stereo-electroencephalography (SEEG) during interictal periods in 8 patients affected by drug-resistant focal epilepsy. We identified IEDs in the SEEG and systematically analyzed the time-locked signals on the EEG by means of evoked potentials, topographical analysis and Electrical Source Imaging (ESI). The dataset has been standardized and is being publicly shared. RESULTS: Our results showed that IEDs that were not clearly visible at single-trials could be uncovered by averaging, in line with previous reports. They also showed that their topographical voltage distributions matched the position of the SEEG electrode where IEDs had been identified, and that ESI techniques can reconstruct it with an accuracy of ∼2 cm. Finally, the present dataset provides a reference to test the accuracy of different methods and parameters. COMPARISON WITH EXISTING METHODS: Our study is the first to systematically compare ESI methods on simultaneously recorded IEDs, and to share a public resource with in-vivo data for their evaluation. CONCLUSIONS: Simultaneous HD-EEG and SEEG recordings can unveil hidden IEDs whose origins can be reconstructed using topographical and ESI analyses, but results depend on the selected methods and parameters.

Connectome spectrum electromagnetic tomography: A method to reconstruct electrical brain source networks at high-spatial resolution.

Connectome spectrum electromagnetic tomography (CSET) combines diffusion MRI-derived structural connectivity data with well-established graph signal processing tools to solve the M/EEG inverse problem. Using simulated EEG signals from fMRI responses, and two EEG datasets on visual-evoked potentials, we provide evidence supporting that (i) CSET captures realistic neurophysiological patterns with better accuracy than state-of-the-art methods, (ii) CSET can reconstruct brain responses more accurately and with more robustness to intrinsic noise in the EEG signal. These results demonstrate that CSET offers high spatio-temporal accuracy, enabling neuroscientists to extend their research beyond the current limitations of low sampling frequency in functional MRI and the poor spatial resolution of M/EEG.

Interictal EEG source connectivity to localize the epileptogenic zone in patients with drug-resistant epilepsy: A machine learning approach.

OBJECTIVE: To deconstruct the epileptogenic networks of patients with drug-resistant epilepsy (DRE) using source functional connectivity (FC) analysis; unveil the FC biomarkers of the epileptogenic zone (EZ); and develop machine learning (ML) models to estimate the EZ using brief interictal electroencephalography (EEG) data. METHODS: We analyzed scalp EEG from 50 patients with DRE who had surgery. We reconstructed the activity (electrical source imaging [ESI]) of virtual sensors (VSs) across the whole cortex and computed FC separately for epileptiform and non-epileptiform EEG epochs (with or without spikes). In patients with good outcome (Engel 1a), four cortical regions were defined: EZ (resection) and three non-epileptogenic zones (NEZs) in the same and opposite hemispheres. Region-specific FC features in six frequency bands and three spatial ranges (long, short, inner) were compared between regions (Wilcoxon sign-rank). We developed ML classifiers to identify the VSs in the EZ using VS-specific FC features. Cross-validation was performed using good outcome data. Performance was compared with poor outcomes and interictal spike localization. RESULTS: FC differed between EZ and NEZs (p < .05) during non-epileptiform and epileptiform epochs, showing higher FC in the EZ than its homotopic contralateral NEZ. During epileptiform epochs, the NEZ in the epileptogenic hemisphere showed higher FC than its contralateral NEZ. In good outcome patients, the ML classifiers reached 75% accuracy to the resection (91% sensitivity; 74% specificity; distance from EZ: 38 mm) using epileptiform epochs (gamma and beta frequency bands) and 62% accuracy using broadband non-epileptiform epochs, both outperforming spike localization (accuracy = 47%; p < .05; distance from EZ: 57 mm). Lower performance was seen in poor outcomes. SIGNIFICANCE: We present an FC approach to extract EZ biomarkers from brief EEG data. Increased FC in various frequencies characterized the EZ during epileptiform and non-epileptiform epochs. FC-based ML models identified the resection better in good than poor outcome patients, demonstrating their potential for presurgical use in pediatric DRE.

High-density electroencephalographic functional networks in genetic generalized epilepsy: Preserved whole-brain topology hides local reorganization.

OBJECTIVE: Genetic generalized epilepsy (GGE) accounts for approximately 20% of adult epilepsy cases and is considered a disorder of large brain networks, involving both hemispheres. Most studies have not shown any difference in functional whole-brain network topology when compared to healthy controls. Our objective was to examine whether this preserved global network topology could hide local reorganizations that balance out at the global network level. METHODS: We recorded high-density electroencephalograms from 20 patients and 20 controls, and reconstructed the activity of 118 regions. We computed functional connectivity in windows free of interictal epileptiform discharges in broad, delta, theta, alpha, and beta frequency bands, characterized the network topology, and used the Hub Disruption Index (HDI) to quantify the topological reorganization. We examined the generalizability of our results by reproducing a 25-electrode clinical system. RESULTS: Our study did not reveal any significant change in whole-brain network topology among GGE patients. However, the HDI was significantly different between patients and controls in all frequency bands except alpha (p < .01, false discovery rate [FDR] corrected, d < -1), and accompanied by an increase in connectivity in the prefrontal regions and default mode network. This reorganization suggests that regions that are important in transferring the information in controls were less so in patients. Inversely, the crucial regions in patients are less so in controls. These findings were also found in delta and theta frequency bands when using 25 electrodes (p < .001, FDR corrected, d < -1). SIGNIFICANCE: In GGE patients, the overall network topology is similar to that of healthy controls but presents a balanced local topological reorganization. This reorganization causes the prefrontal areas and default mode network to be more integrated and segregated, which may explain executive impairment associated with GGE. Additionally, the reorganization distinguishes patients from controls even when using 25 electrodes, suggesting its potential use as a diagnostic tool.

High-density electric source imaging patterns and outcomes following temporal lobectomy in patients with hippocampal sclerosis.

OBJECTIVE: The objective of this study was to ascertain specific patterns of electrical source imaging (ESI) that are associated with a good surgical outcome (no seizure recurrence) using 256-channel high-density (HD) electroencephalography (EEG) in patients with temporal lobe epilepsy with hippocampal sclerosis (TLE-HS) who underwent temporal lobectomy. METHODS: Adult patients (age ≥ 18 years) were prospectively recruited from September 2016 to May 2020 at the authors' center. All patients underwent phase I presurgical evaluation and were subsequently advised to proceed with surgery based on consensus from a multidisciplinary epilepsy conference, without knowing HD-ESI results. All recruited patients were followed for at least 12 months after surgery. The outcome of interest was a status of no seizure recurrence, which was assessed at the end of the study. The association between ESI patterns and outcome was assessed using the chi-square or Fisher exact test. Associated p values as well as odds ratios and 95% CIs were reported. The diagnostic performance of the significant pattern associated with the outcome was also evaluated. RESULTS: Fifty-eight patients with known predictors for either good or worse surgical outcomes were recruited. The mean postoperative follow-up period was 33.34 (SD 13.88) months. Forty-seven patients had sufficient interictal epileptiform discharges for HD-ESI analysis. Thirteen of these 47 patients experienced seizure recurrence. The most common source localizations were at Brodmann area (BA) 20 (inferior temporal area) and BA 21 (middle temporal area). A specific ESI pattern of BA 21 without extratemporal sources was significantly associated with no seizure recurrence (p = 0.047). This pattern had a high positive predictive value of 100% and false-positive rate of 0% associated with no seizure recurrence following the surgery. CONCLUSIONS: A specific ESI pattern that was highly associated with no seizure recurrence following surgery was demonstrated by a 256-channel HD-EEG. If this pattern can be reproducibly proven in further studies, some TLE-HS patients may be able to proceed with surgery without further investigations.

Sleep Spindle Generation Before and After Epilepsy Surgery: A Source Imaging Study in Children with Drug-Resistant Epilepsy.

INTRODUCTION: Literature lacks studies investigating the cortical generation of sleep spindles in drug-resistant epilepsy (DRE) and how they evolve after resection of the epileptogenic zone (EZ). Here, we examined sleep EEGs of children with focal DRE who became seizure-free after focal epilepsy surgery, and aimed to investigate the changes in the spindle generation before and after the surgery using low-density scalp EEG and electrical source imaging (ESI). METHODS: We analyzed N2-sleep EEGs from 19 children with DRE before and after surgery. We identified slow (8-12 Hz) and fast spindles (13-16 Hz), computed their spectral features and cortical generators through ESI and computed their distance from the EZ and irritative zone (IZ). We performed two-way ANOVA testing the effect of spindle type (slow vs. fast) and surgical phase (pre-surgery vs. post-surgery) on each feature. RESULTS: Power, frequency and cortical activation of slow spindles increased after surgery (p < 0.005), while this was not seen for fast spindles. Before surgery, the cortical generators of slow spindles were closer to the EZ (57.3 vs. 66.2 mm, p = 0.007) and IZ (41.3 vs. 55.5 mm, p = 0.02) than fast spindle generators. CONCLUSIONS: Our data indicate alterations in the EEG slow spindles after resective epilepsy surgery. Fast spindle generation on the contrary did not change after surgery. Although the study is limited by its retrospective nature, lack of healthy controls, and reduced cortical spatial sampling, our findings suggest a spatial relationship between the slow spindles and the epileptogenic generators.

Effects of inverse methods and spike phases on interictal high-density EEG source reconstruction.

OBJECTIVE: To determine the effect of inverse methods and timepoints of interictal epileptic discharges (IEDs) used for high-density electric source imaging (hd-ESI) in pharmacoresistant focal epilepsies. METHODS: We retrospectively evaluated the hd-ESI and [18F]fluorodeoxyglucose positron emission tomography (18FDG-PET) of 21 operated patients with pharmacoresistant focal epilepsy (Engel I). Volumetric hd-ESI was performed with three different inverse methods such as the inverse solution linearly constrained minimum variance (LCMV, a beamformer method), standardized low resolution electromagnetic tomography (sLORETA) and weighted minimum-norm estimation (wMNE) and at different IED phases. Hd-ESI accuracy was determined by volumetric overlap and distance between hd-ESI source maximum, as well as 18FDG-PET hypometabolic region relative to the resection zone (RZ). RESULTS: In our cohort, the shortest distances and greatest volumetric overlaps to the RZ were found in the half-rise and peak-phase for all inverse methods. The distance to the RZ was not different between the centroid of the clinical hypothesis-based cluster and the source maximum in peak-phase. However, the distance of the hypothesis-based cluster was significantly shorter compared to the cluster selected by the smallest p-value. CONCLUSIONS: Hd-ESI provides the greatest accuracy in determining the RZ at the IED half-rise and peak-phase for all applied inverse methods, whereby sLORETA and LCMV were equally accurate. SIGNIFICANCE: Our results offer guidance in selecting inverse methods and IED phases for hd-ESI, compare the performance of hd-ESI and 18FDG-PET and encourage future studies in investigating the relationship between interictal ESI and 18FDG-PET hypometabolism.

Dynamic networks differentiate the language ability of children with cochlear implants.

Background: Cochlear implantation (CI) in prelingually deafened children has been shown to be an effective intervention for developing language and reading skill. However, there is a substantial proportion of the children receiving CI who struggle with language and reading. The current study-one of the first to implement electrical source imaging in CI population was designed to identify the neural underpinnings in two groups of CI children with good and poor language and reading skill. Methods: Data using high density electroencephalography (EEG) under a resting state condition was obtained from 75 children, 50 with CIs having good (HL) or poor language skills (LL) and 25 normal hearing (NH) children. We identified coherent sources using dynamic imaging of coherent sources (DICS) and their effective connectivity computing time-frequency causality estimation based on temporal partial directed coherence (TPDC) in the two CI groups compared to a cohort of age and gender matched NH children. Findings: Sources with higher coherence amplitude were observed in three frequency bands (alpha, beta and gamma) for the CI groups when compared to normal hearing children. The two groups of CI children with good (HL) and poor (LL) language ability exhibited not only different cortical and subcortical source profiles but also distinct effective connectivity between them. Additionally, a support vector machine (SVM) algorithm using these sources and their connectivity patterns for each CI group across the three frequency bands was able to predict the language and reading scores with high accuracy. Interpretation: Increased coherence in the CI groups suggest overall that the oscillatory activity in some brain areas become more strongly coupled compared to the NH group. Moreover, the different sources and their connectivity patterns and their association to language and reading skill in both groups, suggest a compensatory adaptation that either facilitated or impeded language and reading development. The neural differences in the two groups of CI children may reflect potential biomarkers for predicting outcome success in CI children.

Image data and computational grids for computing brain shift and solving the electrocorticography forward problem.

This article describes the dataset applied in the research reported in NeuroImage article "Patient-specific solution of the electrocorticography forward problem in deforming brain" [1] that is available for download from the Zenodo data repository (https://zenodo.org/record/7687631) [2]. Preoperative structural and diffusion-weighted magnetic resonance (MR) and postoperative computed tomography (CT) images of a 12-year-old female epilepsy patient under evaluation for surgical intervention were obtained retrospectively from Boston Children's Hospital. We used these images to conduct the analysis at The University of Western Australia's Intelligent Systems for Medicine Laboratory using SlicerCBM [3], our open-source software extension for the 3D Slicer medical imaging platform. As part of the analysis, we processed the images to extract the patient-specific brain geometry; created computational grids, including a tetrahedral grid for the meshless solution of the biomechanical model and a regular hexahedral grid for the finite element solution of the electrocorticography forward problem; predicted the postoperative MRI and DTI that correspond to the brain configuration deformed by the placement of subdural electrodes using biomechanics-based image warping; and solved the patient-specific electrocorticography forward problem to compute the electric potential distribution within the patient's head using the original preoperative and predicted postoperative image data. The well-established and open-source file formats used in this dataset, including Nearly Raw Raster Data (NRRD) files for images, STL files for surface geometry, and Visualization Toolkit (VTK) files for computational grids, allow other research groups to easily reuse the data presented herein to solve the electrocorticography forward problem accounting for the brain shift caused by implantation of subdural grid electrodes.

From theory to practical fundamentals of electroencephalographic source imaging in localizing the epileptogenic zone.

With continued advancement in computational technologies, the analysis of electroencephalography (EEG) has shifted from pure visual analysis to a noninvasive computational technique called EEG source imaging (ESI), which involves mathematical modeling of dipolar and distributed sources of a given scalp EEG pattern. ESI is a noninvasive phase I test for presurgical localization of the seizure onset zone in focal epilepsy. It is a relatively inexpensive modality, as it leverages scalp EEG and magnetic resonance imaging (MRI) data already collected typically during presurgical evaluation. With an adequate number of electrodes and combined with patient-specific MRI-based head models, ESI has proven to be a valuable and accurate clinical diagnostic tool for localizing the epileptogenic zone. Despite its advantages, however, ESI is routinely used at only a minority of epilepsy centers. This paper reviews the current evidence and practical fundamentals for using ESI of interictal and ictal epileptic activity during the presurgical evaluation of drug-resistant patients. We identify common errors in processing and interpreting ESI studies, describe the differences in approach needed for localizing interictal and ictal EEG discharges through practical examples, and describe best practices for optimizing the diagnostic information available from these studies.

Virtual implantation using conventional scalp EEG delineates seizure onset and predicts surgical outcome in children with epilepsy.

OBJECTIVE: Delineation of the seizure onset zone (SOZ) is required in children with drug resistant epilepsy (DRE) undergoing neurosurgery. Intracranial EEG (icEEG) serves as gold standard but has limitations. Here, we examine the utility of virtual implantation with electrical source imaging (ESI) on ictal scalp EEG for mapping the SOZ and predict surgical outcome. METHODS: We retrospectively analyzed EEG data from 35 children with DRE who underwent surgery and dichotomized into seizure-free (SF) and non-seizure-free (NSF). We estimated virtual sensors (VSs) at brain locations that matched icEEG implantation and compared ictal patterns at VSs vs icEEG. We calculated the agreement between VSs SOZ and clinically defined SOZ and built receiver operating characteristic (ROC) curves to test whether it predicted outcome. RESULTS: Twenty-one patients were SF after surgery. Moderate agreement between virtual and icEEG patterns was observed (kappa = 0.45, p < 0.001). Virtual SOZ agreement with clinically defined SOZ was higher in SF vs NSF patients (66.6% vs 41.6%, p = 0.01). Anatomical concordance of virtual SOZ with clinically defined SOZ predicted outcome (AUC = 0.73; 95% CI: 0.57-0.89; sensitivity = 66.7%; specificity = 78.6%; accuracy = 71.4%). CONCLUSIONS: Virtual implantation on ictal scalp EEG can approximate the SOZ and predict outcome. SIGNIFICANCE: SOZ mapping with VSs may contribute to tailoring icEEG implantation and predict outcome.

Source EEG reveals that Rolandic epilepsy is a regional epileptic encephalopathy.

Rolandic epilepsy is the most common form of epileptic encephalopathy, characterized by sleep-potentiated inferior Rolandic epileptiform spikes, seizures, and cognitive deficits in school-age children that spontaneously resolve by adolescence. We recently identified a paucity of sleep spindles, physiological thalamocortical rhythms associated with sleep-dependent learning, in the Rolandic cortex during the active phase of this disease. Because spindles are generated in the thalamus and amplified through regional thalamocortical circuits, we hypothesized that: 1) deficits in spindle rate would involve but extend beyond the inferior Rolandic cortex in active epilepsy and 2) regional spindle deficits would better predict cognitive function than inferior Rolandic spindle deficits alone. To test these hypotheses, we obtained high-resolution MRI, high-density EEG recordings, and focused neuropsychological assessments in children with Rolandic epilepsy during active (n = 8, age 9-14.7 years, 3F) and resolved (seizure free for > 1 year, n = 10, age 10.3-16.7 years, 1F) stages of disease and age-matched controls (n = 8, age 8.9-14.5 years, 5F). Using a validated spindle detector applied to estimates of electrical source activity in 31 cortical regions, including the inferior Rolandic cortex, during stages 2 and 3 of non-rapid eye movement sleep, we compared spindle rates in each cortical region across groups. Among detected spindles, we compared spindle features (power, duration, coherence, bilateral synchrony) between groups. We then used regression models to examine the relationship between spindle rate and cognitive function (fine motor dexterity, phonological processing, attention, and intelligence, and a global measure of all functions). We found that spindle rate was reduced in the inferior Rolandic cortices in active but not resolved disease (active P = 0.007; resolved P = 0.2) compared to controls. Spindles in this region were less synchronous between hemispheres in the active group (P = 0.005; resolved P = 0.1) compared to controls; but there were no differences in spindle power, duration, or coherence between groups. Compared to controls, spindle rate in the active group was also reduced in the prefrontal, insular, superior temporal, and posterior parietal regions (i.e., "regional spindle rate", P < 0.039 for all). Independent of group, regional spindle rate positively correlated with fine motor dexterity (P < 1e-3), attention (P = 0.02), intelligence (P = 0.04), and global cognitive performance (P < 1e-4). Compared to the inferior Rolandic spindle rate alone, models including regional spindle rate trended to improve prediction of global cognitive performance (P = 0.052), and markedly improved prediction of fine motor dexterity (P = 0.006). These results identify a spindle disruption in Rolandic epilepsy that extends beyond the epileptic cortex and a potential mechanistic explanation for the broad cognitive deficits that can be observed in this epileptic encephalopathy.

Low density electrical source imaging of the ictal onset zone in the surgical evaluation of children with epilepsy.

PURPOSE: To investigate the utility of Low Density (LD) Electrical Source Imaging (ESI) to model the ictal onset zone (IOZ) for the surgical work up of children with medically refractory epilepsy. METHODS: This was a retrospective review of 12 patients from a district and regional pediatric epilepsy center, who underwent focal resections between 2014 and 2019. ESI was generated using the Curry 8 software, incorporating T1 Magnetic Resonance Imaging (MRI) scans and scalp electroencephalogram (EEG) recordings. Concordance of the ictal LD-ESI localizations to the epileptogenic zone was assessed by comparing the location of the ictal LD-ESI to the focal resection margins on neuroimaging and noting the post-operative outcomes at one year. Localizations determined by ictal LD-ESI were also compared to interictal LD-ESI, positron emission tomography (FDG-PET) and interictal magnetoencephalography (MEG). RESULTS: Ictal ESI correctly localized the ictal onset zone in 4/6 patients, with all four being seizure free at one year. Similarly, interictal ESI localized the irritative zone in 7/9 patients with focal resections, with 6/7 being seizure free at one year. Additionally, we observed ictal ESI to be concordant to interictal ESI in 5/6 patients. Ictal ESI and interictal ESI were concordant to interictal MEG in 3/6 patients. Ictal ESI was concordant with FDG-PET in 6/7 cases. CONCLUSION: IOZ source localization through LD-ESI is a promising complementary method of assessing the epileptogenic focus in children. These findings may support the inclusion of ictal LD-ESI within the pre-surgical evaluation of children to supplement current diagnostic tools.

Automated electrical source imaging with scalp EEG to define the insular irritative zone: Comparison with simultaneous intracranial EEG.

OBJECTIVE: To evaluate the accuracy of automatedinterictallow-density electrical source imaging (LD-ESI) to define the insular irritative zone (IZ) by comparing the simultaneous interictal ESI localization with the SEEG interictal activity. METHODS: Long-term simultaneous scalp electroencephalography (EEG) and stereo-EEG (SEEG) with at least one depth electrode exploring the operculo-insular region(s) were analyzed. Automated interictal ESI was performed on the scalp EEG using standardized low-resolution brain electromagnetic tomography (sLORETA) and individual head models. A two-step analysis was performed: i) sublobar concordance betweencluster-based ESI localization and SEEG-based IZ; ii) time-locked ESI-/SEEG analysis. Diagnostic accuracy values were calculated using SEEG as reference standard. Subgroup analysis wascarried out, based onthe involvement of insular contacts in the seizure onset and patterns of insular interictal activity. RESULTS: Thirty patients were included in the study. ESI showed an overall accuracy of 53% (C.I. 29-76%). Sensitivity and specificity were calculated as 53% (C.I. 29-76%), 55% (C.I. 23-83%) respectively. Higher accuracy was found in patients with frequent and dominant interictal insular spikes. CONCLUSIONS: LD-ESI defines with good accuracy the insular implication in the IZ, which is not possible with classical interictalscalpEEG interpretation. SIGNIFICANCE: Automated LD-ESI may be a valuable additional tool to characterize the epileptogenic zone in epilepsies with suspected insular involvement.

Investigation of dynamic functional connectivity of the source reconstructed epileptiform discharges in focal epilepsy: A graph theory approach.

OBJECTIVE: The aim of the present study is to investigate with noninvasive methods the modulation of dynamic functional connectivity during interictal epileptiform discharge (IED). METHOD: We reconstructed the cortical source of the EEG recorded IED of 17 patients with focal epilepsy. We then computed dynamic connectivity using the time resolved phase locking value (PLV). We derived graph theory indices (i.e. degree, strength, local efficiency, clustering coefficient and global efficiency). Finally, we selected the atlas node with the maximum activation as the IED cortical source investigating the graph indices dynamics in theta, alpha, beta and gamma frequency bands. RESULTS: We observed IED-locked modulations of the graph indexes depending on the frequency bands. We detected a modulation of the strength, clustering coefficient, local and global efficiency both in theta and in alpha bands, which also displayed modulations of the degree index. In the beta band only the global efficiency was modulated by the IED, while no effects were detected in the gamma band. Finally, we found a correlation between alpha and theta local efficiency, as well as alpha global efficiency, and the epilepsy duration. SIGNIFICANCE: Our findings suggest that the neural synchronization is not limited to the IED cortical source, but implies a phase synchronization across multiple brain areas. We hypothesize that the aberrant electrical activity originating from the IED locus is spread amongst the other network nodes throughout the low frequency bands (i.e. theta and alpha). Moreover, IED-dependent increase in the global efficiency indicates that the IED interfere with the whole network functioning. We finally discussed possible application of this methodology for future investigation.

ConvDip: A Convolutional Neural Network for Better EEG Source Imaging.

The electroencephalography (EEG) is a well-established non-invasive method in neuroscientific research and clinical diagnostics. It provides a high temporal but low spatial resolution of brain activity. To gain insight about the spatial dynamics of the EEG, one has to solve the inverse problem, i.e., finding the neural sources that give rise to the recorded EEG activity. The inverse problem is ill-posed, which means that more than one configuration of neural sources can evoke one and the same distribution of EEG activity on the scalp. Artificial neural networks have been previously used successfully to find either one or two dipole sources. These approaches, however, have never solved the inverse problem in a distributed dipole model with more than two dipole sources. We present ConvDip, a novel convolutional neural network (CNN) architecture, that solves the EEG inverse problem in a distributed dipole model based on simulated EEG data. We show that (1) ConvDip learned to produce inverse solutions from a single time point of EEG data and (2) outperforms state-of-the-art methods on all focused performance measures. (3) It is more flexible when dealing with varying number of sources, produces less ghost sources and misses less real sources than the comparison methods. It produces plausible inverse solutions for real EEG recordings from human participants. (4) The trained network needs <40 ms for a single prediction. Our results qualify ConvDip as an efficient and easy-to-apply novel method for source localization in EEG data, with high relevance for clinical applications, e.g., in epileptology and real-time applications.

Fast oscillations >40 Hz localize the epileptogenic zone: An electrical source imaging study using high-density electroencephalography.

OBJECTIVE: Fast Oscillations (FO) >40 Hz are a promising biomarker of the epileptogenic zone (EZ). Evidence using scalp electroencephalography (EEG) remains scarce. We assessed if electrical source imaging of FO using 256-channel high-density EEG (HD-EEG) is useful for EZ identification. METHODS: We analyzed HD-EEG recordings of 10 focal drug-resistant epilepsy patients with seizure-free postsurgical outcome. We marked FO candidate events at the time of epileptic spikes and verified them by screening for an isolated peak in the time-frequency plot. We performed electrical source imaging of spikes and FO within the Maximum Entropy of the Mean framework. Source localization maps were validated against the surgical cavity. RESULTS: We identified FO in five out of 10 patients who had a superficial or intermediate deep generator. The maximum of the FO maps was localized inside the cavity in all patients (100%). Analysis with a reduced electrode coverage using the 10-10 and 10-20 system showed a decreased localization accuracy of 60% and 40% respectively. CONCLUSIONS: FO recorded with HD-EEG localize the EZ. HD-EEG is better suited to detect and localize FO than conventional EEG approaches. SIGNIFICANCE: This study acts as proof-of-concept that FO localization using 256-channel HD-EEG is a viable marker of the EZ.

Functional connectivity of the irritative zone identified by electrical source imaging, and EEG-correlated fMRI analyses.

OBJECTIVE: The irritative zone - the area generating epileptic spikes - can be studied non-invasively during the interictal period using Electrical Source Imaging (ESI) and simultaneous electroencephalography-functional magnetic resonance imaging (EEG-fMRI). Although the techniques yield results which may overlap spatially, differences in spatial localization of the irritative zone within the same patient are consistently observed. To investigate this discrepancy, we used Blood Oxygenation Level Dependent (BOLD) functional connectivity measures to examine the underlying relationship between ESI and EEG-fMRI findings. METHODS: Fifteen patients (age 20-54), who underwent presurgical epilepsy investigation, were scanned using a single-session resting-state EEG-fMRI protocol. Structural MRI was used to obtain the electrode localisation of a high-density 64-channel EEG cap. Electrical generators of interictal epileptiform discharges were obtained using a distributed local autoregressive average (LAURA) algorithm as implemented in Cartool EEG software. BOLD activations were obtained using both spike-related and voltage-map EEG-fMRI analysis. The global maxima of each method were used to investigate the temporal relationship of BOLD time courses and to assess the spatial similarity using the Dice similarity index between functional connectivity maps. RESULTS: ESI, voltage-map and spike-related EEG-fMRI methods identified peaks in 15 (100%), 13 (67%) and 8 (53%) of the 15 patients, respectively. For all methods, maxima were localised within the same lobe, but differed in sub-lobar localisation, with a median distance of 22.8 mm between the highest peak for each method. The functional connectivity analysis showed that the temporal correlation between maxima only explained 38% of the variance between the time course of the BOLD response at the maxima. The mean Dice similarity index between seed-voxel functional connectivity maps showed poor spatial agreement. SIGNIFICANCE: Non-invasive methods for the localisation of the irritative zone have distinct spatial and temporal sensitivity to different aspects of the local cortical network involved in the generation of interictal epileptiform discharges.

Temporal and Spatial Dynamics of Different Interictal Epileptic Discharges: A Time-Frequency EEG Approach in Pediatric Focal Refractory Epilepsy.

Objective: Characterization of the spatial and temporal dynamics of interictal epileptic discharges (IED) using time-frequency analysis (TFA) and electrical-source localization (ESL). Methods: TFA was performed on IED (spikes, spike waves, and polyspike waves) recorded by high-density-EEG (HD-EEG) in 19 refractory focal epileptic children. Temporal modulations related to IEDs were analyzed in a time window around the IED peaks [-1,000 to 1,000 ms]. Spatial modulations were analyzed by ESL in the time-frequency and time domains. Results: IED were associated with complex power spectral modulations. We observed increases in power spectrum (IPS) patterns specific to IED type. For spikes, the TFA pattern consisted of an IPS (-100 to +100 ms, 4-50 Hz). For spike waves, the IPS was followed by a second IPS (+100 to +400 ms, 4-10 Hz), corresponding to the slow wave. IPS patterns were preceded (-400 to -100 ms, 4-40 Hz), and followed (+100 to +400 ms) by a decrease in the power spectrum (DPS) (n = 8). For 14 out of 19 patients, at least one ESL method was concordant with the epileptogenic area. For the remaining five patients, all of them had temporal epilepsies. ESL in the time-frequency domain (DPS/IPS) provided concordant (n = 6) or complementary (n = 4) information to the ESL in the time domain concerning the epileptogenic zone. ESL in time-frequency domain (DPS/IPS) was the only method to provide concordant information concerning the epileptogenic zone in three patients. Significance: TFA demonstrates complex time-frequency modulations of the neuronal networks around IED, suggesting that the pathological mechanisms are initiated well before onset of the classical hyper-synchronization of the IED. Combining time and time-frequency analysis of the ESL provides complementary information to define the epileptogenic zone in refractory focal epilepsy.

Pushing the Limits of EEG: Estimation of Large-Scale Functional Brain Networks and Their Dynamics Validated by Simultaneous fMRI.

Functional magnetic resonance imaging (fMRI) is the technique of choice for detecting large-scale functional brain networks and to investigate their dynamics. Because fMRI measures brain activity indirectly, electroencephalography (EEG) has been recently considered a feasible tool for detecting such networks, particularly the resting-state networks (RSNs). However, a truly unbiased validation of such claims is still missing, which can only be accomplished by using simultaneously acquired EEG and fMRI data, due to the spontaneous nature of the activity underlying the RSNs. Additionally, EEG is still poorly explored for the purpose of mapping task-specific networks, and no studies so far have been focused on investigating networks' dynamic functional connectivity (dFC) with EEG. Here, we started by validating RSNs derived from the continuous reconstruction of EEG sources by directly comparing them with those derived from simultaneous fMRI data of 10 healthy participants, and obtaining an average overlap (quantified by the Dice coefficient) of 0.4. We also showed the ability of EEG to map the facial expressions processing network (FEPN), highlighting regions near the posterior superior temporal sulcus, where the FEPN is anchored. Then, we measured the dFC using EEG for the first time in this context, estimated dFC brain states using dictionary learning, and compared such states with those obtained from the fMRI. We found a statistically significant match between fMRI and EEG dFC states, and determined the existence of two matched dFC states which contribution over time was associated with the brain activity at the FEPN, showing that the dynamics of FEPN can be captured by both fMRI and EEG. Our results push the limits of EEG toward being used as a brain imaging tool, while supporting the growing literature on EEG correlates of (dynamic) functional connectivity measured with fMRI, and providing novel insights into the coupling mechanisms underlying the two imaging techniques.