Electric source imaging of interictal activity accurately localises the seizure onset zone

Visual analysis of high density EEG: As good as electrical source imaging?

•

Visual analysis of HD-EEG is an excellent tool to explore the epileptogenic focus.

•

ESI remains the gold standard for presurgical evaluation of the cortical source.

•

ESI at 50% slope/ESI at peak discordance could predict worse surgical outcome.

Objective

In this study, we sought to determine whether visual analysis of high density EEG (HD-EEG) would provide similar localizing information comparable to electrical source imaging (ESI).

Methods

HD-EEG (256 electrodes) recordings from 20 patients suffering from unifocal, drug-resistant epilepsy (13 women, mean age 29.1 ± 2.62 years, 11 with temporal lobe epilepsy) were examined. In the visual analysis condition, we identified the 5 contacts with maximal spike amplitude and determined their localization with respect to the underlying cortex. ESI was computed using the LAURA algorithm of the averaged spikes in the patient’s individual MRI. We considered the localization “correct” if all 5 contacts were concordant with the resection volume underneath or if ESI was located within the resection as determined by the postoperative MRI.

Results

Twelve patients were postoperatively seizure-free (Engel Class IA), while the remaining eight were in class IB to IV. Visual analysis and ESI showed sensitivity of 58% and 75%, specificity of 75% and 87%, and accuracy of 65% and 80%, respectively. In 70% of cases, visual analysis and ESI provided concordant results.

Conclusions

Localization of the electrodes with maximal spike amplitude provides very good estimation of the localization of the underlying source. However, ESI has a higher accuracy and adds 3D information; therefore, it should remain the tool of choice for presurgical evaluation.

Significance

The present study proposes the possibility to analyze HD-EEG visually, in tandem with ESI or alone, if ESI is not accessible.

The network integration of epileptic activity in relation to surgical outcome

OBJECTIVE

Epilepsy is a network disease with epileptic activity and cognitive impairment involving large-scale brain networks. A complex network is involved in the seizure and in the interictal epileptiform discharges (IEDs). Directed connectivity analysis, describing the information transfer between brain regions, and graph analysis are applied to high-density EEG to characterise networks.

METHODS

We analysed 19 patients with focal epilepsy who had high-density EEG containing IED and underwent surgery. We estimated cortical activity during IED using electric source analysis in 72 atlas-based cortical regions of the individual brain MRI. We applied directed connectivity analysis (information Partial Directed Coherence) and graph analysis on these sources and compared patients with good vs poor post-operative outcome at global, hemispheric and lobar level.

RESULTS

We found lower network integration reflected by global, hemispheric, lobar efficiency during the IED (p < 0.05) in patients with good post-surgical outcome, compared to patients with poor outcome. Prediction was better than using the IED field or the localisation obtained by electric source imaging.

CONCLUSIONS

Abnormal network patterns in epilepsy are related to seizure outcome after surgery.

SIGNIFICANCE

Our finding may help understand networks related to a more "isolated" epileptic activity, limiting the extent of the epileptic network in patients with subsequent good post-operative outcome.

Interictal epileptogenic zone localization in patients with focal epilepsy using electric source imaging and directed functional connectivity from low-density EEG

OBJECTIVE

Electrical source imaging (ESI) is used increasingly to estimate the epileptogenic zone (EZ) in patients with epilepsy. Directed functional connectivity (DFC) coupled to ESI helps to better characterize epileptic networks, but studies on interictal activity have relied on high-density recordings. We investigated the accuracy of ESI and DFC for localizing the EZ, based on low-density clinical electroencephalography (EEG).

METHODS

We selected patients with the following: (a) focal epilepsy, (b) interictal spikes on standard EEG, (c) either a focal structural lesion concordant with the electroclinical semiology or good postoperative outcome. In 34 patients (20 temporal lobe epilepsy [TLE], 14 extra-TLE [ETLE]), we marked interictal spikes and estimated the cortical activity during each spike in 82 cortical regions using a patient-specific head model and distributed linear inverse solution. DFC between brain regions was computed using Granger-causal modeling followed by network topologic measures. The concordance with the presumed EZ at the sublobar level was computed using the epileptogenic lesion or the resected area in postoperative seizure-free patients.

RESULTS

ESI, summed outflow, and efficiency were concordant with the presumed EZ in 76% of the patients, whereas the clustering coefficient and betweenness centrality were concordant in 70% of patients. There was no significant difference between ESI and connectivity measures. In all measures, patients with TLE had a significantly higher (P < 0.05) concordance with the presumed EZ than patients with with ETLE. The brain volume accepted for concordance was significantly larger in TLE.

SIGNIFICANCE

ESI and DFC derived from low-density EEG can reliably estimate the EZ from interictal spikes. Connectivity measures were not superior to ESI for EZ localization during interictal spikes, but the current validation of the localization of connectivity measure is promising for other applications.

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.

EEG Source Imaging: A Practical Review of the Analysis Steps

The electroencephalogram (EEG) is one of the oldest technologies to measure neuronal activity of the human brain. It has its undisputed value in clinical diagnosis, particularly (but not exclusively) in the identification of epilepsy and sleep disorders and in the evaluation of dysfunctions in sensory transmission pathways. With the advancement of digital technologies, the analysis of EEG has moved from pure visual inspection of amplitude and frequency modulations over time to a comprehensive exploration of the temporal and spatial characteristics of the recorded signals. Today, EEG is accepted as a powerful tool to capture brain function with the unique advantage of measuring neuronal processes in the time frame in which these processes occur, namely in the sub-second range. However, it is generally stated that EEG suffers from a poor spatial resolution that makes it difficult to infer to the location of the brain areas generating the neuronal activity measured on the scalp. This statement has challenged a whole community of biomedical engineers to offer solutions to localize more precisely and more reliably the generators of the EEG activity. High-density EEG systems combined with precise information of the head anatomy and sophisticated source localization algorithms now exist that convert the EEG to a true neuroimaging modality. With these tools in hand and with the fact that EEG still remains versatile, inexpensive and portable, electrical neuroimaging has become a widely used technology to study the functions of the pathological and healthy human brain. However, several steps are needed to pass from the recording of the EEG to 3-dimensional images of neuronal activity. This review explains these different steps and illustrates them in a comprehensive analysis pipeline integrated in a stand-alone freely available academic software: Cartool. The information about how the different steps are performed in Cartool is only meant as a suggestion. Other EEG source imaging software may apply similar or different approaches to the different steps.

Electrical Neuroimaging of Music Processing in Pianists With and Without True Absolute Pitch

True absolute pitch (AP), labeling of pitches with semitone precision without a reference, is classically studied using isolated tones. However, AP is acquired and has its function within complex dynamic musical contexts. Here we examined event-related brain responses and underlying cerebral sources to endings of short expressive string quartets, investigating a homogeneous population of young highly trained pianists with half of them possessing true-AP. The pieces ended regularly or contained harmonic transgressions at closure that participants appraised. Given the millisecond precision of ERP analyses, this experimental plan allowed examining whether AP alters music processing at an early perceptual, or later cognitive level, or both, and which cerebral sources underlie differences with non-AP musicians. We also investigated the impact of AP on general auditory cognition. Remarkably, harmonic transgression sensitivity did not differ between AP and non-AP participants, and differences for auditory cognition were only marginal. The key finding of this study is the involvement of a microstate peaking around 60 ms after musical closure, characterizing AP participants. Concurring sources were estimated in secondary auditory areas, comprising the planum temporale, all transgression conditions collapsed. These results suggest that AP is not a panacea to become a proficient musician, but a rare perceptual feature.

Electroencephalography, magnetoencephalography and source localization: their value in epilepsy

PURPOSE OF REVIEW

Source localization of cerebral activity using electroencephalography (EEG) or magnetoencephalography (MEG) can reveal noninvasively the generators of the abnormal signals recorded in epilepsy, such as interictal epileptic discharges (IEDs) and seizures. Here, we review recent progress showcasing the usefulness of these techniques in treating patients with drug-resistant epilepsy.

RECENT FINDINGS

The source localization of IEDs by high-density EEG and MEG has now been proved in large patient cohorts to be accurate and clinically relevant, with positive and negative predictive values rivaling those of structural MRI. Localizing seizure onsets is an emerging technique that seems to perform similarly well to the localization of interictal spikes, although there remain questions regarding the processing of signals for reliable results. The localization of somatosensory cortex using EEG/MEG is well established. The localization of language cortex is less reliable, although progress has been made regarding hemispheric lateralization. Source localization is also able to reveal how epilepsy alters the dynamics of neuronal activity in the large-scale networks that underlie cerebral function.

SUMMARY

Given the high performance of EEG/MEG source localization, these tools should find a place similar to that of established techniques like MRI in the assessment of patients for epilepsy surgery.

Subcortical electrophysiological activity is detectable with high-density EEG source imaging

Subcortical neuronal activity is highly relevant for mediating communication in large-scale brain networks. While electroencephalographic (EEG) recordings provide appropriate temporal resolution and coverage to study whole brain dynamics, the feasibility to detect subcortical signals is a matter of debate. Here, we investigate if scalp EEG can detect and correctly localize signals recorded with intracranial electrodes placed in the centromedial thalamus, and in the nucleus accumbens. Externalization of deep brain stimulation (DBS) electrodes, placed in these regions, provides the unique opportunity to record subcortical activity simultaneously with high-density (256 channel) scalp EEG. In three patients during rest with eyes closed, we found significant correlation between alpha envelopes derived from intracranial and EEG source reconstructed signals. Highest correlation was found for source signals in close proximity to the actual recording sites, given by the DBS electrode locations. Therefore, we present direct evidence that scalp EEG indeed can sense subcortical signals.

Electroencephalography (EEG) allows the measurement of electrical signals associated with brain activity, but it is unclear if EEG can accurately measure subcortical activity. Here, the authors show that source dynamics, reconstructed from scalp EEG, correlate with activity recorded from human thalamus and nucleus accumbens.

Large-Scale 3-5 Hz Oscillation Constrains the Expression of Neocortical Fast Ripples in a Mouse Model of Mesial Temporal Lobe Epilepsy

Large-scale slow oscillations allow the integration of neuronal activity across brain regions during sensory or cognitive processing. However, evidence that this form of coding also holds for pathological networks, such as for distributed networks in epileptic disorders, does not yet exist. Here, we show in a mouse model of unilateral hippocampal epilepsy that epileptic fast ripples generated in the neocortex distant from the primary focus occur during transient trains of interictal epileptic discharges. During these epileptic paroxysms, local phase-locking of neuronal firing and a phase-amplitude coupling of the epileptic discharges over a slow oscillation at 3-5 Hz are detected. Furthermore, the buildup of the slow oscillation begins in the bihippocampal network that includes the focus, which synchronizes and drives the activity across the large-scale epileptic network into the frontal cortex. This study provides the first functional description of the emergence of neocortical fast ripples in hippocampal epilepsy and shows that cross-frequency coupling might be a fundamental mechanism underlying the spreading of epileptic activity.

Diagnostic accuracy of interictal source imaging in presurgical epilepsy evaluation: A systematic review from the E-PILEPSY consortium

OBJECTIVE

Interictal high resolution (HR-) electric source imaging (ESI) and magnetic source imaging (MSI) are non-invasive tools to aid epileptogenic zone localization in epilepsy surgery candidates. We carried out a systematic review on the diagnostic accuracy and quality of evidence of these modalities.

METHODS

Embase, Pubmed and the Cochrane database were searched on 13 February 2017. Diagnostic accuracy studies taking post-surgical seizure outcome as reference standard were selected. Quality appraisal was based on the QUADAS-2 framework.

RESULTS

Eleven studies were included: eight MSI (n = 267), three HR-ESI (n = 127) studies. None was free from bias. This mostly involved: selection of operated patients only, interference of source imaging with surgical decision, and exclusion of indeterminate results. Summary sensitivity and specificity estimates were 82% (95% CI: 75-88%) and 53% (95% CI: 37-68%) for overall source imaging, with no statistical difference between MSI and HR-ESI. Specificity is higher when partially concordant results were included as non-concordant (p < 0.05). Inclusion of indeterminate test results as non-concordant lowered sensitivity (p < 0.05).

CONCLUSIONS

Source imaging has a relatively high sensitivity but low specificity for identification of the epileptogenic zone.

SIGNIFICANCE

We need higher quality studies allowing unbiased test evaluation to determine the added value and diagnostic accuracy of source imaging in the presurgical workup of refractory focal epilepsy.

Assessing the localization accuracy and clinical utility of electric and magnetic source imaging in children with epilepsy

OBJECTIVE

To evaluate the accuracy and clinical utility of conventional 21-channel EEG (conv-EEG), 72-channel high-density EEG (HD-EEG) and 306-channel MEG in localizing interictal epileptiform discharges (IEDs).

METHODS

Twenty-four children who underwent epilepsy surgery were studied. IEDs on conv-EEG, HD-EEG, MEG and intracranial EEG (iEEG) were localized using equivalent current dipoles and dynamical statistical parametric mapping (dSPM). We compared the localization error (E_Loc) with respect to the ground-truth Irritative Zone (IZ), defined by iEEG sources, between non-invasive modalities and the distance from resection (D_res) between good- (Engel 1) and poor-outcomes. For each patient, we estimated the resection percentage of IED sources and tested whether it predicted outcome.

RESULTS

MEG presented lower E_Loc than HD-EEG and conv-EEG. For all modalities, D_res was shorter in good-outcome than poor-outcome patients, but only the resection percentage of the ground-truth IZ and MEG-IZ predicted surgical outcome.

CONCLUSIONS

MEG localizes the IZ more accurately than conv-EEG and HD-EEG. MSI may help the presurgical evaluation in terms of patient's outcome prediction. The promising clinical value of ESI for both conv-EEG and HD-EEG prompts the use of higher-density EEG-systems to possibly achieve MEG performance.

SIGNIFICANCE

Localizing the IZ non-invasively with MSI/ESI facilitates presurgical evaluation and surgical prognosis assessment.

Interictal spike connectivity in human epileptic neocortex

OBJECTIVE

Interictal spikes are a biomarker of epilepsy, yet their precise roles are poorly understood. Using long-term neocortical recordings from epileptic patients, we investigated the spatial-temporal propagation patterns of interictal spiking.

METHODS

Interictal spikes were detected in 10 epileptic patients. Short time direct directed transfer function was used to map the spatial-temporal patterns of interictal spike onset and propagation across different cortical topographies.

RESULTS

Each patient had unique interictal spike propagation pattern that was highly consistent across times, regardless of the frequency band. High spiking brain regions were often not spike onset regions. We observed frequent spike propagations to shorter distances and that the central sulcus forms a strong barrier to spike propagation. Spike onset and seizure onset seemed to be distinct networks in most cases.

CONCLUSIONS

Patients in epilepsy have distinct and unique network of causal propagation pattern which are very consistent revealing the underlying epileptic network. Although spike are epileptic biomarkers, spike origin and seizure onset seems to be distinct in most cases.

SIGNIFICANCE

Understanding patterns of interictal spike propagation could lead to the identification patient-specific epileptic networks amenable to surgical or other treatments.

Estimating EEG Source Dipole Orientation Based on Singular-value Decomposition for Connectivity Analysis

In the last decade, the use of high-density electrode arrays for EEG recordings combined with the improvements of source reconstruction algorithms has allowed the investigation of brain networks dynamics at a sub-second scale. One powerful tool for investigating large-scale functional brain networks with EEG is time-varying effective connectivity applied to source signals obtained from electric source imaging. Due to computational and interpretation limitations, the brain is usually parcelled into a limited number of regions of interests (ROIs) before computing EEG connectivity. One specific need and still open problem is how to represent the time- and frequency-content carried by hundreds of dipoles with diverging orientation in each ROI with one unique representative time-series. The main aim of this paper is to provide a method to compute a signal that explains most of the variability of the data contained in each ROI before computing, for instance, time-varying connectivity. As the representative time-series for a ROI, we propose to use the first singular vector computed by a singular-value decomposition of all dipoles belonging to the same ROI. We applied this method to two real datasets (visual evoked potentials and epileptic spikes) and evaluated the time-course and the frequency content of the obtained signals. For each ROI, both the time-course and the frequency content of the proposed method reflected the expected time-course and the scalp-EEG frequency content, representing most of the variability of the sources (~ 80%) and improving connectivity results in comparison to other procedures used so far. We also confirm these results in a simulated dataset with a known ground truth.

Compressibility of High-Density EEG Signals in Stroke Patients

Stroke is a critical event that causes the disruption of neural connections. There is increasing evidence that the brain tries to reorganize itself and to replace the damaged circuits, by establishing compensatory pathways. Intra- and extra-cellular currents are involved in the communication between neurons and the macroscopic effects of such currents can be detected at the scalp through electroencephalographic (EEG) sensors. EEG can be used to study the lesions in the brain indirectly, by studying their effects on the brain electrical activity. The primary goal of the present work was to investigate possible asymmetries in the activity of the two hemispheres, in the case one of them is affected by a lesion due to stroke. In particular, the compressibility of High-Density-EEG (HD-EEG) recorded at the two hemispheres was investigated since the presence of the lesion is expected to impact on the regularity of EEG signals. The secondary objective was to evaluate if standard low density EEG is able to provide such information. Eighteen patients with unilateral stroke were recruited and underwent HD-EEG recording. Each EEG signal was compressively sensed, using Block Sparse Bayesian Learning, at increasing compression rate. The two hemispheres showed significant differences in the compressibility of EEG. Signals acquired at the electrode locations of the affected hemisphere showed a better reconstruction quality, quantified by the Structural SIMilarity index (SSIM), than the EEG signals recorded at the healthy hemisphere (p < 0.05), for each compression rate value. The presence of the lesion seems to induce an increased regularity in the electrical activity of the brain, thus an increased compressibility.

Directed functional connections underlying spontaneous brain activity

Neuroimaging studies have shown that spontaneous brain activity is characterized as changing networks of coherent activity across multiple brain areas. However, the directionality of functional interactions between the most active regions in our brain at rest remains poorly understood. Here, we examined, at the whole-brain scale, the main drivers and directionality of interactions that underlie spontaneous human brain activity by applying directed functional connectivity analysis to electroencephalography (EEG) source signals. We found that the main drivers of electrophysiological activity were the posterior cingulate cortex (PCC), the medial temporal lobes (MTL), and the anterior cingulate cortex (ACC). Among those regions, the PCC was the strongest driver and had both the highest integration and segregation importance, followed by the MTL regions. The driving role of the PCC and MTL resulted in an effective directed interaction directed from posterior toward anterior brain regions. Our results strongly suggest that the PCC and MTL structures are the main drivers of electrophysiological spontaneous activity throughout the brain and suggest that EEG-based directed functional connectivity analysis is a promising tool to better understand the dynamics of spontaneous brain activity in healthy subjects and in various brain disorders.

Automatic ictal onset source localization in presurgical epilepsy evaluation

OBJECTIVE

To test the diagnostic accuracy of a new automatic algorithm for ictal onset source localization (IOSL) during routine presurgical epilepsy evaluation following STARD (Standards for Reporting of Diagnostic Accuracy) criteria.

METHODS

We included 28 consecutive patients with refractory focal epilepsy (25 patients with temporal lobe epilepsy (TLE) and 3 with extratemporal epilepsy) who underwent resective epilepsy surgery. Ictal EEG patterns were analyzed with a novel automatic IOSL algorithm. IOSL source localizations on a sublobar level were validated by comparison with actual resection sites and seizure free outcome 2 years after surgery.

RESULTS

Sensitivity of IOSL was 92.3% (TLE: 92.3%); specificity 60% (TLE: 50%); positive predictive value 66.7% (TLE: 66.7%); and negative predictive value 90% (TLE: 85.7%). The likelihood ratio was more than ten times higher for concordant IOSL results as compared to discordant results (p = 0.013).

CONCLUSIONS

We demonstrated the clinical feasibility of our IOSL approach yielding reasonable high performance measures on a sublobar level.

SIGNIFICANCE

Our IOSL method may contribute to a correct localization of the seizure onset zone in temporal lobe epilepsy and can readily be used in standard epilepsy monitoring settings. Further studies are needed for validation in extratemporal epilepsy.

Electrical Neuroimaging of Music Processing Reveals Mid-Latency Changes with Level of Musical Expertise

This original research focused on the effect of musical training intensity on cerebral and behavioral processing of complex music using high-density event-related potential (ERP) approaches. Recently we have been able to show progressive changes with training in gray and white matter, and higher order brain functioning using (f)MRI [(functional) Magnetic Resonance Imaging], as well as changes in musical and general cognitive functioning. The current study investigated the same population of non-musicians, amateur pianists and expert pianists using spatio-temporal ERP analysis, by means of microstate analysis, and ERP source imaging. The stimuli consisted of complex musical compositions containing three levels of transgression of musical syntax at closure that participants appraised. ERP waveforms, microstates and underlying brain sources revealed gradual differences according to musical expertise in a 300–500 ms window after the onset of the terminal chords of the pieces. Within this time-window, processing seemed to concern context-based memory updating, indicated by a P3b-like component or microstate for which underlying sources were localized in the right middle temporal gyrus, anterior cingulate and right parahippocampal areas. Given that the 3 expertise groups were carefully matched for demographic factors, these results provide evidence of the progressive impact of training on brain and behavior.

Dysfunctional gaze processing in bipolar disorder

Gaze conveys emotional information, and humans present sensitivity to its direction from the earliest days of life. Bipolar disorder is a disease characterized by fluctuating states of emotional and cognitive dysregulation. To explore the role of attentional control on face processing in bipolar patients (BP) we used gaze direction as an emotion modulation parameter in a two-back Working Memory (WM) task while high-density EEG data were acquired. Since gaze direction influences emotional attributions to faces with neutral expressions as well, we presented neutral faces with direct and averted gaze. Nineteen euthymic BP and a sample of age- and gender-matched controls were examined.

In BP we observed diminished P200 and augmented P300 evoked responses, differentially modulated by non-repeated or repeated faces, as well as by gaze direction. BP showed a reduced P200 amplitude, significantly stronger for faces with direct gaze than averted gaze. Source localization of P200 indicated decreased activity in sensory-motor regions and frontal areas suggestive of abnormal affective processing of neutral faces.

The present study provides neurophysiological evidence for abnormal gaze processing in BP and suggests dysfunctional processing of direct eye contact as a prominent characteristic of bipolar disorder.

•

This ERP study identified abnormalities in gaze processing in bipolar patients.

•

We observed functional anomalies in the P200 and P300 evoked responses.

•

BP showed a strong suppression of the P200 for faces with direct gaze.

•

Source localization indicated decreased activity in sensory-motor regions.

Automated long-term EEG analysis to localize the epileptogenic zone

Objective

We investigated the performance of automatic spike detection and subsequent electroencephalogram (EEG) source imaging to localize the epileptogenic zone (EZ) from long-term EEG recorded during video-EEG monitoring.

Methods

In 32 patients, spikes were automatically detected in the EEG and clustered according to their morphology. The two spike clusters with most single events in each patient were averaged and localized in the brain at the half-rising time and peak of the spike using EEG source imaging. On the basis of the distance from the sources to the resection and the known patient outcome after surgery, the performance of the automated EEG analysis to localize the EZ was quantified.

Results

In 28 out of the 32 patients, the automatically detected spike clusters corresponded with the reported interictal findings. The median distance to the resection in patients with Engel class I outcome was 6.5 and 15 mm for spike cluster 1 and 27 and 26 mm for cluster 2, at the peak and the half-rising time of the spike, respectively. Spike occurrence (cluster 1 vs. cluster 2) and spike timing (peak vs. half-rising) significantly influenced the distance to the resection (p < 0.05). For patients with Engel class II, III, and IV outcomes, the median distance increased to 36 and 36 mm for cluster 1. Localizing spike cluster 1 at the peak resulted in a sensitivity of 70% and specificity of 100%, positive prediction value (PPV) of 100%, and negative predictive value (NPV) of 53%. Including the results of spike cluster 2 led to an increased sensitivity of 79% NPV of 55% and diagnostic OR of 11.4, while the specificity dropped to 75% and the PPV to 90%.

Significance

We showed that automated analysis of long-term EEG recordings results in a high sensitivity and specificity to localize the epileptogenic focus.

Seizure Onset Zone Localization from Ictal High-Density EEG in Refractory Focal Epilepsy

Epilepsy surgery is the most efficient treatment option for patients with refractory epilepsy. Before surgery, it is of utmost importance to accurately delineate the seizure onset zone (SOZ). Non-invasive EEG is the most used neuroimaging technique to diagnose epilepsy, but it is hard to localize the SOZ from EEG due to its low spatial resolution and because epilepsy is a network disease, with several brain regions becoming active during a seizure. In this work, we propose and validate an approach based on EEG source imaging (ESI) combined with functional connectivity analysis to overcome these problems. We considered both simulations and real data of patients. Ictal epochs of 204-channel EEG and subsets down to 32 channels were analyzed. ESI was done using realistic head models and LORETA was used as inverse technique. The connectivity pattern between the reconstructed sources was calculated, and the source with the highest number of outgoing connections was selected as SOZ. We compared this algorithm with a more straightforward approach, i.e. selecting the source with the highest power after ESI as the SOZ. We found that functional connectivity analysis estimated the SOZ consistently closer to the simulated EZ/RZ than localization based on maximal power. Performance, however, decreased when 128 electrodes or less were used, especially in the realistic data. The results show the added value of functional connectivity analysis for SOZ localization, when the EEG is obtained with a high-density setup. Next to this, the method can potentially be used as objective tool in clinical settings.

Directed Functional Brain Connectivity Based on EEG Source Imaging: Methodology and Application to Temporal Lobe Epilepsy

OBJECTIVE

The importance of functional brain connectivity to study physiological and pathological brain activity has been widely recognized. Here, we aimed to 1) review a methodological pipeline to investigate directed functional connectivity between brain regions using source signals derived from high-density EEG; 2) elaborate on some methodological challenges; and 3) apply this pipeline to temporal lobe epilepsy (TLE) patients and healthy controls to investigate directed functional connectivity differences in the theta and beta frequency bands during EEG epochs without visible pathological activity.

METHODS

The methodological pipeline includes: EEG acquisition and preprocessing, electrical-source imaging (ESI) using individual head models and distributed inverse solutions, parcellation of the gray matter in regions of interest, fixation of the dipole orientation for each region, computation of the spectral power in the source space, and directed functional connectivity estimation using Granger-causal modeling. We specifically analyzed how the signal-to-noise ratio (SNR) changes using different approaches for the dipole orientation fixation. We applied this pipeline to 20 left TLE patients, 20 right TLE patients, and 20 healthy controls.

RESULTS

Projecting each dipole to the predominant dipole orientation leads to a threefold SNR increase as compared to the norm of the dipoles. By comparing connectivity in TLE versus controls, we found significant frequency-specific outflow differences in physiologically plausible regions.

CONCLUSION

The results suggest that directed functional connectivity derived from ESI can help better understand frequency-specific resting-state network alterations underlying focal epilepsy.

SIGNIFICANCE

EEG-based directed functional connectivity could contribute to the search of new biomarkers of this disorder.

EEG Derived Brain Activity Reflects Treatment Response from Vagus Nerve Stimulation in Patients with Epilepsy

The mechanism of action of vagus nerve stimulation (VNS) is yet to be elucidated. To that end, the effects of VNS on the brain of epileptic patients were studied. Both when VNS was switched "On" and "Off", the brain activity of responders (R, seizure frequency reduction of over 50%) was compared to the brain activity of nonresponders (NR, seizure frequency reduction of less than 50%). Using EEG recordings, a significant increase in P300 amplitude for R and a significant decrease in P300 amplitude for NR were found. We found biomarkers for checking the efficacy of VNS with accuracy up to 94%. The results show that P300 features recorded in nonmidline electrodes are better P300 biomarkers for VNS efficacy than P300 features recorded in midline electrodes. Using source localization and connectivity analyses, the activity of the limbic system, insula and orbitofrontal cortex was found to be dependent on VNS switched "On" versus "Off" or patient group (R versus NR). The results suggest an important role for these areas in the mechanism of action of VNS, although a larger patient study should be done to confirm the findings.

Magnetoencephalographic signatures of insular epileptic spikes based on functional connectivity

Failure to recognize insular cortex seizures has recently been identified as a cause of epilepsy surgeries targeting the temporal, parietal, or frontal lobe. Such failures are partly due to the fact that current noninvasive localization techniques fare poorly in recognizing insular epileptic foci. Our group recently demonstrated that magnetoencephalography (MEG) is sensitive to epileptiform spikes generated by the insula. In this study, we assessed the potential of distributed source imaging and functional connectivity analyses to distinguish insular networks underlying the generation of spikes. Nineteen patients with operculo-insular epilepsy were investigated. Each patient underwent MEG as well as T1-weighted magnetic resonance imaging (MRI) as part of their standard presurgical evaluation. Cortical sources of MEG spikes were reconstructed with the maximum entropy on the mean algorithm, and their time courses served to analyze source functional connectivity. The results indicate that the anterior and posterior subregions of the insula have specific patterns of functional connectivity mainly involving frontal and parietal regions, respectively. In addition, while their connectivity patterns are qualitatively similar during rest and during spikes, couplings within these networks are much stronger during spikes. These results show that MEG can establish functional connectivity-based signatures that could help in the diagnosis of different subtypes of insular cortex epilepsy. Hum Brain Mapp 37:3250-3261, 2016. © 2016 Wiley Periodicals, Inc.

Epileptic networks are strongly connected with and without the effects of interictal discharges

OBJECTIVE

Epilepsy is increasingly considered as the dysfunction of a pathologic neuronal network (epileptic network) rather than a single focal source. We aimed to assess the interactions between the regions that comprise the epileptic network and to investigate their dependence on the occurrence of interictal epileptiform discharges (IEDs).

METHODS

We analyzed resting state simultaneous electroencephalography-functional magnetic resonance imaging (EEG-fMRI) recordings in 10 patients with drug-resistant focal epilepsy with multifocal IED-related blood oxygen level-dependent (BOLD) responses and a maximum t-value in the IED field. We computed functional connectivity (FC) maps of the epileptic network using two types of seed: (1) a 10-mm diameter sphere centered in the global maximum of IED-related BOLD map, and (2) the independent component with highest correlation to the IED-related BOLD map, named epileptic component. For both approaches, we compared FC maps before and after regressing out the effect of IEDs in terms of maximum and mean t-values and percentage of map overlap.

RESULTS

Maximum and mean FC maps t-values were significantly lower after regressing out IEDs at the group level (p < 0.01). Overlap extent was 85% ± 12% and 87% ± 12% when the seed was the 10-mm diameter sphere and the epileptic component, respectively.

SIGNIFICANCE

Regions involved in a specific epileptic network show coherent BOLD fluctuations independent of scalp EEG IEDs. FC topography and strength is largely preserved by removing the IED effect. This could represent a signature of a sustained pathologic network with contribution from epileptic activity invisible to the scalp EEG.

The Role of Skull Modeling in EEG Source Imaging for Patients with Refractory Temporal Lobe Epilepsy

We investigated the influence of different skull modeling approaches on EEG source imaging (ESI), using data of six patients with refractory temporal lobe epilepsy who later underwent successful epilepsy surgery. Four realistic head models with different skull compartments, based on finite difference methods, were constructed for each patient: (i) Three models had skulls with compact and spongy bone compartments as well as air-filled cavities, segmented from either computed tomography (CT), magnetic resonance imaging (MRI) or a CT-template and (ii) one model included a MRI-based skull with a single compact bone compartment. In all patients we performed ESI of single and averaged spikes marked in the clinical 27-channel EEG by the epileptologist. To analyze at which time point the dipole estimations were closer to the resected zone, ESI was performed at two time instants: the half-rising phase and peak of the spike. The estimated sources for each model were validated against the resected area, as indicated by the postoperative MRI. Our results showed that single spike analysis was highly influenced by the signal-to-noise ratio (SNR), yielding estimations with smaller distances to the resected volume at the peak of the spike. Although averaging reduced the SNR effects, it did not always result in dipole estimations lying closer to the resection. The proposed skull modeling approaches did not lead to significant differences in the localization of the irritative zone from clinical EEG data with low spatial sampling density. Furthermore, we showed that a simple skull model (MRI-based) resulted in similar accuracy in dipole estimation compared to more complex head models (based on CT- or CT-template). Therefore, all the considered head models can be used in the presurgical evaluation of patients with temporal lobe epilepsy to localize the irritative zone from low-density clinical EEG recordings.

Brain imaging in the assessment for epilepsy surgery

Brain imaging has a crucial role in the presurgical assessment of patients with epilepsy. Structural imaging reveals most cerebral lesions underlying focal epilepsy. Advances in MRI acquisitions including diffusion-weighted imaging, post-acquisition image processing techniques, and quantification of imaging data are increasing the accuracy of lesion detection. Functional MRI can be used to identify areas of the cortex that are essential for language, motor function, and memory, and tractography can reveal white matter tracts that are vital for these functions, thus reducing the risk of epilepsy surgery causing new morbidities. PET, SPECT, simultaneous EEG and functional MRI, and electrical and magnetic source imaging can be used to infer the localisation of epileptic foci and assist in the design of intracranial EEG recording strategies. Progress in semi-automated methods to register imaging data into a common space is enabling the creation of multimodal three-dimensional patient-specific datasets. These techniques show promise for the demonstration of the complex relations between normal and abnormal structural and functional data and could be used to direct precise intracranial navigation and surgery for individual patients.

Ecstatic Epileptic Seizures: A Glimpse into the Multiple Roles of the Insula

Ecstatic epileptic seizures are a rare but compelling epileptic entity. During the first seconds of these seizures, ecstatic auras provoke feelings of well-being, intense serenity, bliss, and "enhanced self-awareness." They are associated with the impression of time dilation, and can be described as a mystic experience by some patients. The functional neuroanatomy of ecstatic seizures is still debated. During recent years several patients presenting with ecstatic auras have been reported by others and us (in total n = 52); a few of them in the setting of presurgical evaluation including electrical brain stimulation. According to the recently recognized functions of the insula, and the results of nuclear brain imaging and electrical stimulation, the ecstatic symptoms in these patients seem to localize to a functional network centered around the anterior insular cortex, where we thus propose to locate this rare ictal phenomenon. Here we summarize the role of the multiple sensory, autonomic, affective, and cognitive functions of the insular cortex, which are integrated into the creation of self-awareness, and we suggest how this system may become dysfunctional on several levels during ecstatic aura.

Electrical source imaging of interictal spikes using multiple sparse volumetric priors for presurgical epileptogenic focus localization

Electrical source imaging of interictal spikes observed in EEG recordings of patients with refractory epilepsy provides useful information to localize the epileptogenic focus during the presurgical evaluation. However, the selection of the time points or time epochs of the spikes in order to estimate the origin of the activity remains a challenge. In this study, we consider a Bayesian EEG source imaging technique for distributed sources, i.e. the multiple volumetric sparse priors (MSVP) approach. The approach allows to estimate the time courses of the intensity of the sources corresponding with a specific time epoch of the spike. Based on presurgical averaged interictal spikes in six patients who were successfully treated with surgery, we estimated the time courses of the source intensities for three different time epochs: (i) an epoch starting 50 ms before the spike peak and ending at 50% of the spike peak during the rising phase of the spike, (ii) an epoch starting 50 ms before the spike peak and ending at the spike peak and (iii) an epoch containing the full spike time period starting 50 ms before the spike peak and ending 230 ms after the spike peak. To identify the primary source of the spike activity, the source with the maximum energy from 50 ms before the spike peak till 50% of the spike peak was subsequently selected for each of the time windows. For comparison, the activity at the spike peaks and at 50% of the peaks was localized using the LORETA inversion technique and an ECD approach. Both patient-specific spherical forward models and patient-specific 5-layered finite difference models were considered to evaluate the influence of the forward model. Based on the resected zones in each of the patients, extracted from post-operative MR images, we compared the distances to the resection border of the estimated activity. Using the spherical models, the distances to the resection border for the MSVP approach and each of the different time epochs were in the same range as the LORETA and ECD techniques. We found distances smaller than 23 mm, with robust results for all the patients. For the finite difference models, we found that the distances to the resection border for the MSVP inversions of the full spike time epochs were generally smaller compared to the MSVP inversions of the time epochs before the spike peak. The results also suggest that the inversions using the finite difference models resulted in slightly smaller distances to the resection border compared to the spherical models. The results we obtained are promising because the MSVP approach allows to study the network of the estimated source-intensities and allows to characterize the spatial extent of the underlying sources.

•

A Bayesian ESI technique is evaluated to localize interictal spike activity.

•

Averaged spikes in six patients were used that were seizure free after surgery.

•

We compared the technique with the LORETA an ECD technique.

•

We evaluated both spherical and 5-layered finite difference forward models.

•

Our approach is potentially useful to delineate the irritative zone.

All-in-one interictal presurgical imaging in patients with epilepsy: single-session EEG/PET/(f)MRI

PURPOSE

In patients with pharmacoresistant focal epilepsy, resection of the epileptic focus can lead to freedom from seizures or significant improvement in well-selected candidates. Localization of the epileptic focus with multimodal concordance is crucial for a good postoperative outcome. Beyond the detection of epileptogenic lesions on structural MRI and focal hypometabolism on FDG PET, EEG-based Electric Source Imaging (ESI) and simultaneous EEG and functional MRI (EEG-fMRI) are increasingly applied for mapping epileptic activity. We here report presurgical multimodal interictal imaging using a hybrid PET/MR scanner for single-session FDG PET, MRI, EEG-fMRI and ESI.

METHODS

This quadrimodal imaging procedure was performed in a single session in 12 patients using a high-density (256 electrodes) MR-compatible EEG system and a hybrid PET/MR scanner. EEG was used to exclude subclinical seizures during uptake of the PET tracer, to compute ESI on interictal epileptiform discharges and to guide fMRI analysis for mapping haemodynamic changes correlated with interictal epileptiform activity.

RESULTS

The whole multimodal recording was performed in less than 2 hours with good patient comfort and data quality. Clinically contributory examinations with at least two modalities were obtained in nine patients and with all modalities in five patients.

CONCLUSION

This single-session quadrimodal imaging procedure provided reliable and contributory interictal clinical data. This procedure avoids multiple scanning sessions and is associated with less radiation exposure than PET-CT. Moreover, it guarantees the same medication level and medical condition for all modalities. The procedure improves workflow and could reduce the duration and cost of presurgical epilepsy evaluations.

Dynamic directed interictal connectivity in left and right temporal lobe epilepsy

OBJECTIVE

There is increasing evidence that epileptic activity involves widespread brain networks rather than single sources and that these networks contribute to interictal brain dysfunction. We investigated the fast-varying behavior of epileptic networks during interictal spikes in right and left temporal lobe epilepsy (RTLE and LTLE) at a whole-brain scale using directed connectivity.

METHODS

In 16 patients, 8 with LTLE and 8 with RTLE, we estimated the electrical source activity in 82 cortical regions of interest (ROIs) using high-density electroencephalography (EEG), individual head models, and a distributed linear inverse solution. A multivariate, time-varying, and frequency-resolved Granger-causal modeling (weighted Partial Directed Coherence) was applied to the source signal of all ROIs. A nonparametric statistical test assessed differences between spike and baseline epochs. Connectivity results between RTLE and LTLE were compared between RTLE and LTLE and with neuropsychological impairments.

RESULTS

Ipsilateral anterior temporal structures were identified as key drivers for both groups, concordant with the epileptogenic zone estimated invasively. We observed an increase in outflow from the key driver already before the spike. There were also important temporal and extratemporal ipsilateral drivers in both conditions, and contralateral only in RTLE. A different network pattern between LTLE and RTLE was found: in RTLE there was a much more prominent ipsilateral to contralateral pattern than in LTLE. Half of the RTLE patients but none of the LTLE patients had neuropsychological deficits consistent with contralateral temporal lobe dysfunction, suggesting a relationship between connectivity changes and cognitive deficits.

SIGNIFICANCE

The different patterns of time-varying connectivity in LTLE and RTLE suggest that they are not symmetrical entities, in line with our neuropsychological results. The highest outflow region was concordant with invasive validation of the epileptogenic zone. This enhanced characterization of dynamic connectivity patterns could better explain cognitive deficits and help the management of epilepsy surgery candidates.

Sleep affects cortical source modularity in temporal lobe epilepsy: A high-density EEG study

OBJECTIVE

Interictal epileptiform discharges (IEDs) constitute a perturbation of ongoing cerebral rhythms, usually more frequent during sleep. The aim of the study was to determine whether sleep influences the spread of IEDs over the scalp and whether their distribution depends on vigilance-related modifications in cortical interactions.

METHODS

Wake and sleep 256-channel electroencephalography (EEG) data were recorded in 12 subjects with right temporal lobe epilepsy (TLE) differentiated by whether they had mesial or neocortical TLE. Spikes were selected during wake and sleep. The averaged waking signal was subtracted from the sleep signal and projected on a bidimensional scalp map; sleep and wake spike distributions were compared by using a t-test. The superimposed signal of sleep and wake traces was obtained; the rising phase of the spike, the peak, and the deflections following the spike were identified, and their cortical generator was calculated using low-resolution brain electromagnetic tomography (LORETA) for each group.

RESULTS

A mean of 21 IEDs in wake and 39 in sleep per subject were selected. As compared to wake, a larger IED scalp projection was detected during sleep in both mesial and neocortical TLE (p<0.05). A series of EEG deflections followed the spike, the cortical sources of which displayed alternating activations of different cortical areas in wake, substituted by isolated, stationary activations in sleep in mesial TLE and a silencing in neocortical TLE.

CONCLUSION

During sleep, the IED scalp region increases, while cortical interaction decreases.

SIGNIFICANCE

The interaction of cortical modules in sleep and wake in TLE may influence the appearance of IEDs on scalp EEG; in addition, IEDs could be proxies for cerebral oscillation perturbation.

Mapping epileptic activity: sources or networks for the clinicians?

Epileptic seizures of focal origin are classically considered to arise from a focal epileptogenic zone and then spread to other brain regions. This is a key concept for semiological electro-clinical correlations, localization of relevant structural lesions, and selection of patients for epilepsy surgery. Recent development in neuro-imaging and electro-physiology and combinations, thereof, have been validated as contributory tools for focus localization. In parallel, these techniques have revealed that widespread networks of brain regions, rather than a single epileptogenic region, are implicated in focal epileptic activity. Sophisticated multimodal imaging and analysis strategies of brain connectivity patterns have been developed to characterize the spatio-temporal relationships within these networks by combining the strength of both techniques to optimize spatial and temporal resolution with whole-brain coverage and directional connectivity. In this paper, we review the potential clinical contribution of these functional mapping techniques as well as invasive electrophysiology in human beings and animal models for characterizing network connectivity.

An interictal EEG spectral metric for temporal lobe epilepsy lateralization

OBJECTIVE

Visually-obvious abnormalities in the resting baseline EEG--slowing, spiking and high-frequency oscillations (HFOs)--are cardinal, though incompletely understood, features of the seizure onset zone in focal epilepsy. We hypothesized that evidence of cortical network dysfunction in temporal lobe epilepsy (TLE) would persist in the absence of visually-classifiable abnormalities in the baseline EEG recorded within the conventional passband, and that metrics of such dysfunction could serve as a lateralizing diagnostic in TLE.

METHODS

Epochs of resting EEG without significant abnormalities in light sleep over several days were compared between a group of 10 patients with proven TLE and 10 subjects without epilepsy. A novel laterality metric computed from the line length of normalized power spectra from the temporal channels was compared between the two groups.

RESULTS

Significant group differences in spectral line length laterality metric were found between the TLE and control group. At the individual level, seven of 10 TLE patients had highly significant laterality metrics, all concordant with the known laterality of their disease.

SIGNIFICANCE

Detailed spectral analysis offers novel insight into TLE network behavior, independent of the orthodox abnormalities of EEG slowing, spikes or HFOs. The results may be deployed in a practical diagnostic manner, offer insight into the EEG manifestations of disordered cellular network architecture in TLE, and maybe understood through simple analogy with the theory of linear time-invariant physical systems.