Comparison of Brain Networks During Interictal Oscillations and Spikes on Magnetoencephalography and Intracerebral EEG

Comparison of beamformer and ICA for dynamic connectivity analysis: A simultaneous MEG-SEEG study

Magnetoencephalography (MEG) is a powerful tool for estimating brain connectivity with both good spatial and temporal resolution. It is particularly helpful in epilepsy to characterize non-invasively the epileptic networks. However, using MEG to map brain networks requires solving a difficult inverse problem that introduces uncertainty in the activity localization and connectivity measures. Our goal here was to compare independent component analysis (ICA) followed by dipole source localization and the linearly constrained minimum-variance beamformer (LCMV-BF) for characterizing regions with interictal epileptic activity and their dynamic connectivity. After a simulation study, we compared ICA and LCMV-BF results with intracerebral EEG (stereotaxic EEG, SEEG) recorded simultaneously in 8 epileptic patients, which provide a unique 'ground truth' to which non-invasive results can be confronted. We compared the signal time courses extracted applying ICA and LCMV-BF on MEG data to that of SEEG, both for the actual signals and the dynamic connectivity computed using cross-correlation (evolution of links in time). With our simulations, we illustrated the different effect of the temporal and spatial correlation among sources on the two methods. While ICA was more affected by the temporal correlation but robust against spatial configurations, LCMV-BF showed opposite behavior. Moreover, ICA seems more suited to retrieve the simulated networks. In case of real patient data, good MEG/SEEG correlation and good localization were obtained in 6 out of 8 patients. In 4 of them ICA had the best performance (higher correlation, lower localization distance). In terms of dynamic connectivity, the evolution in time of the cross-correlation links could be retrieved in 5 patients out of 6, however, with more variable results in terms of correlation and distance. In two patients LCMV-BF had better results than ICA. In one patient the two methods showed equally good outcomes, and in the remaining two patients ICA performed best. In conclusion, our results obtained by exploiting simultaneous MEG/SEEG recordings suggest that ICA and LCMV-BF have complementary qualities for retrieving the dynamics of interictal sources and their network interactions.

Virtual MEG sensors based on beamformer and independent component analysis can reconstruct epileptic activity as measured on simultaneous intracerebral recordings

The prevailing gold standard for presurgical determination of epileptogenic brain networks is intracerebral EEG, a potent yet invasive approach. Magnetoencephalography (MEG) is a state-of-the art non-invasive method for investigating epileptiform discharges. However, it is not clear at what level the precision offered by MEG can reach that of SEEG. Here, we present a strategy for non-invasively retrieving the constituents of the interictal network, with high spatial and temporal precision. Our method is based on MEG and a combination of spatial filtering and independent component analysis (ICA). We validated this approach in twelve patients with drug-resistant focal epilepsy, thanks to the unprecedented ground truth provided by simultaneous recordings of MEG and SEEG. A minimum variance adaptive beamformer estimated the source time series and ICA was used to further decompose these time series into network constituents (MEG-ICs), each having a time series (virtual electrode) and a topography (spatial distribution of amplitudes in the brain). We show that MEG has a considerable sensitivity of 0.80 and 0.84 and a specificity of 0.93 and 0.91 for reconstructing deep and superficial sources, respectively, when compared to the ground truth (SEEG). For each epileptic MEG-IC (n = 131), we found at least one significantly correlating SEEG contact close to zero lag after correcting for multiple comparisons. All the patients except one had at least one epileptic component that was highly correlated (Spearman rho>0.3) with that of SEEG traces. MEG-ICs correlated well with SEEG traces. The strength of correlation coefficients did not depend on the depth of the SEEG contacts or the clinical outcome of the patient. A significant proportion of the MEG-ICs (n = 83/131) were localized in proximity with their maximally correlating SEEG, within a mean distance of 20±12.18mm. Our research is the first to validate the MEG-retrieved beamformer IC sources against SEEG-derived ground truth in a simultaneous MEG-SEEG framework. Observations from the present study suggest that non-invasive MEG source components may potentially provide additional information, comparable to SEEG in a number of instances.

Assessment of Effective Network Connectivity among MEG None Contaminated Epileptic Transitory Events

Characterizing epileptogenic zones EZ (sources responsible of excessive discharges) would assist a neurologist during epilepsy diagnosis. Locating efficiently these abnormal sources among magnetoencephalography (MEG) biomarker is obtained by several inverse problem techniques. These techniques present different assumptions and particular epileptic network connectivity. Here, we proposed to evaluate performances of distributed inverse problem in defining EZ. First, we applied an advanced technique based on Singular Value Decomposition (SVD) to recover only pure transitory activities (interictal epileptiform discharges). We evaluated our technique's robustness in separation between transitory and ripples versus frequency range, transitory shapes, and signal to noise ratio on simulated data (depicting both epileptic biomarkers and respecting time series and spectral properties of realistic data). We validated our technique on MEG signal using detector precision on 5 patients. Then, we applied four methods of inverse problem to define cortical areas and neural generators of excessive discharges. We computed network connectivity of each technique. Then, we confronted obtained noninvasive networks to intracerebral EEG transitory network connectivity using nodes in common, connection strength, distance metrics between concordant nodes of MEG and IEEG, and average propagation delay. Coherent Maximum Entropy on the Mean (cMEM) proved a high matching between MEG network connectivity and IEEG based on distance between active sources, followed by Exact low-resolution brain electromagnetic tomography (eLORETA), Dynamical Statistical Parametric Mapping (dSPM), and Minimum norm estimation (MNE). Clinical performance was interesting for entire methods providing in an average of 73.5% of active sources detected in depth and seen in MEG, and vice versa, about 77.15% of active sources were detected from MEG and seen in IEEG. Investigated problem techniques succeed at least in finding one part of seizure onset zone. dSPM and eLORETA depict the highest connection strength among all techniques. Propagation delay varies in this range [18, 25]ms, knowing that eLORETA ensures the lowest propagation delay (18 ms) and the closet one to IEEG propagation delay.

Presurgical Evaluation of Epilepsy Using Resting-State MEG Functional Connectivity

Epilepsy is caused by abnormal electrical discharges (clinically identified by electrophysiological recording) in a specific part of the brain [originating in only one part of the brain, namely, the epileptogenic zone (EZ)]. Epilepsy is now defined as an archetypical hyperexcited neural network disorder. It can be investigated through the network analysis of interictal discharges, ictal discharges, and resting-state functional connectivity. Currently, there is an increasing interest in embedding resting-state connectivity analysis into the preoperative evaluation of epilepsy. Among the various neuroimaging technologies employed to achieve brain functional networks, magnetoencephalography (MEG) with the excellent temporal resolution is an ideal tool for estimating the resting-state connectivity between brain regions, which can reveal network abnormalities in epilepsy. What value does MEG resting-state functional connectivity offer for epileptic presurgical evaluation? Regarding this topic, this paper introduced the origin of MEG and the workflow of constructing source-space functional connectivity based on MEG signals. Resting-state functional connectivity abnormalities correlate with epileptogenic networks, which are defined by the brain regions involved in the production and propagation of epileptic activities. This paper reviewed the evidence of altered epileptic connectivity based on low- or high-frequency oscillations (HFOs) and the evidence of the advantage of using simultaneous MEG and intracranial electroencephalography (iEEG) recordings. More importantly, this review highlighted that MEG-based resting-state functional connectivity has the potential to predict postsurgical outcomes. In conclusion, resting-state MEG functional connectivity has made a substantial progress toward serving as a candidate biomarker included in epileptic presurgical evaluations.

Dynamic analysis on simultaneous iEEG-MEG data via hidden Markov model

BACKGROUND

Intracranial electroencephalography (iEEG) recordings are used for clinical evaluation prior to surgical resection of the focus of epileptic seizures and also provide a window into normal brain function. A major difficulty with interpreting iEEG results at the group level is inconsistent placement of electrodes between subjects making it difficult to select contacts that correspond to the same functional areas. Recent work using time delay embedded hidden Markov model (HMM) applied to magnetoencephalography (MEG) resting data revealed a distinct set of brain states with each state engaging a specific set of cortical regions. Here we use a rare group dataset with simultaneously acquired resting iEEG and MEG to test whether there is correspondence between HMM states and iEEG power changes that would allow classifying iEEG contacts into functional clusters.

METHODS

Simultaneous MEG-iEEG recordings were performed at rest on 11 patients with epilepsy whose intracranial electrodes were implanted for pre-surgical evaluation. Pre-processed MEG sensor data was projected to source space. Time delay embedded HMM was then applied to MEG time series. At the same time, iEEG time series were analyzed with time-frequency decomposition to obtain spectral power changes with time. To relate MEG and iEEG results, correlations were computed between HMM probability time courses of state activation and iEEG power time course from the mid contact pair for each electrode in equally spaced frequency bins and presented as correlation spectra for the respective states and iEEG channels. Association of iEEG electrodes with HMM states based on significant correlations was compared to that based on the distance to peaks in subject-specific state topographies.

RESULTS

Five HMM states were inferred from MEG. Two of them corresponded to the left and the right temporal activations and had a spectral signature primarily in the theta/alpha frequency band. All the electrodes had significant correlations with at least one of the states (p < 0.05 uncorrected) and for 27/50 electrodes these survived within-subject FDR correction (q < 0.05). These correlations peaked in the theta/alpha band. There was a highly significant dependence between the association of states and electrodes based on functional correlations and that based on spatial proximity (p = 5.6e-6,χ2 test for independence). Despite the potentially atypical functional anatomy and physiological abnormalities related to epilepsy, HMM model estimated from the patient group was very similar to that estimated from healthy subjects.

CONCLUSION

Epilepsy does not preclude HMM analysis of interictal data. The resulting group functional states are highly similar to those reported for healthy controls. Power changes recorded with iEEG correlate with HMM state time courses in the alpha-theta band and the presence of this correlation can be related to the spatial location of electrode contacts close to the individual peaks of the corresponding state topographies. Thus, the hypothesized relation between iEEG contacts and HMM states exists and HMM could be further explored as a method for identifying comparable iEEG channels across subjects for the purposes of group analysis.

A review of magnetoencephalography use in pediatric epilepsy: an update on best practice

Introduction: Magnetoencephalography (MEG) is a noninvasive technique that is used for presurgical evaluation of children with drug-resistant epilepsy (DRE).Areas covered: The contributions of MEG for localizing the epileptogenic zone are discussed, in particular in extra-temporal lobe epilepsy and focal cortical dysplasia, which are common in children, as well as in difficult to localize epilepsy such as operculo-insular epilepsy. Further, the authors review current evidence on MEG for mapping eloquent cortex, its performance, application in clinical practice, and potential challenges.Expert opinion: MEG could change the clinical management of children with DRE by directing placement of intracranial electrodes thereby enhancing their yield. With improved identification of a circumscribed epileptogenic zone, MEG could render more patients as suitable candidates for epilepsy surgery and increase utilization of surgery.

Large-Scale Desynchronization During Interictal Epileptic Discharges Recorded With Intracranial EEG

It is increasingly recognized that deep understanding of epileptic seizures requires both localizing and characterizing the functional network of the region where they are initiated, i. e., the epileptic focus. Previous investigations of the epileptogenic focus' functional connectivity have yielded contrasting results, reporting both pathological increases and decreases during resting periods and seizures. In this study, we shifted paradigm to investigate the time course of connectivity in relation to interictal epileptiform discharges. We recruited 35 epileptic patients undergoing intracranial EEG (iEEG) investigation as part of their presurgical evaluation. For each patient, 50 interictal epileptic discharges (IEDs) were marked and iEEG signals were epoched around those markers. Signals were narrow-band filtered and time resolved phase-locking values were computed to track the dynamics of functional connectivity during IEDs. Results show that IEDs are associated with a transient decrease in global functional connectivity, time-locked to the peak of the discharge and specific to the high range of the gamma frequency band. Disruption of the long-range connectivity between the epileptic focus and other brain areas might be an important process for the generation of epileptic activity. Transient desynchronization could be a potential biomarker of the epileptogenic focus since 1) the functional connectivity involving the focus decreases significantly more than the connectivity outside the focus and 2) patients with good surgical outcome appear to have a significantly more disconnected focus than patients with bad outcomes.

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.

Magnetoencephalography With Optically Pumped 4He Magnetometers at Ambient Temperature

In this paper, we present the first proof of concept confirming the possibility to record magnetoencephalographic (MEG) signals with optically pumped magnetometers (OPMs) based on the parametric resonance of ⁴He atoms. The main advantage of this kind of OPM is the possibility to provide a tri-axis vector measurement of the magnetic field at room-temperature (the ⁴He vapor is neither cooled nor heated). The sensor achieves a sensitivity of 210 fT/ √ Hz in the bandwidth [2-300 Hz]. MEG simulation studies with a brain phantom were cross-validated with real MEG measurements on a healthy subject. For both studies, MEG signal was recorded consecutively with OPMs and superconducting quantum interference devices (SQUIDs) used as reference sensors. For healthy subject MEG recordings, three MEG proofs of concept were carried out: auditory evoked fields, visual evoked fields, and spontaneous activity. M100 peaks have been detected on evoked responses recorded by both OPMs and SQUIDs with no significant difference in latency. Concerning spontaneous activity, an attenuation of the signal power between 8-12 Hz (alpha band) related to eyes opening has been observed with OPM similarly to SQUID. All these results confirm that the room temperature vector ⁴He OPMs can record MEG signals and provide reliable information on brain activity.

Integration of stationary wavelet transform on a dynamic partial reconfiguration for recognition of pre-ictal gamma oscillations

To define the neural networks responsible of an epileptic seizure, it is useful to perform advanced signal processing techniques. In this context, electrophysiological signals present three types of waves: oscillations, spikes, and a mixture of both. Recent studies show that spikes and oscillations should be separated properly in order to define the accurate neural connectivity during the pre-ictal, seizure and inter-ictal states. Retrieving oscillatory activity is a sensitive task due to the frequency overlap between oscillations and transient activities. Advanced filtering techniques have been proposed to ensure a good separation between oscillations and spikes. It would be interesting to apply them in real time for instantaneous monitoring, seizure warning or neurofeedback systems. This requires improving execution time. This constraint can be overcome using embedded systems that combine hardware and software in an optimized architecture.

We propose here to implement a stationary wavelet transform (SWT) as an adaptive filtering technique retaining only pre-ictal gamma oscillations, as validated in previous work, on a partial dynamic configuration. Then, the same architecture is used with further modifications to integrate spatio temporal mapping for an early recognition of seizure build-up.

Data that contains transient, pre-ictal gamma oscillations and a seizure was simulated. the method on real intracerebral signals was also tested. The SWT was integrated on an embedded architecture. This architecture permits a spatio temporal mapping to detect the accurate time and localization of seizure build-up, while reducing computation time by a factor of around 40. Embedded systems are a promising venue for real-time applications in clinical systems for epilepsy.

The interictal mesial temporal lobe epilepsy network

OBJECTIVE

Identification of patient-specific epileptogenic networks is critical to designing successful treatment strategies. Multiple noninvasive methods have been used to characterize epileptogenic networks. However, these methods lack the spatiotemporal resolution to allow precise localization of epileptiform activity. We used intracranial recordings, at much higher spatiotemporal resolution, across a cohort of patients with mesial temporal lobe epilepsy (MTLE) to delineate features common to their epileptogenic networks. We used interictal rather than seizure data because interictal spikes occur more frequently, providing us greater power for analyzing variances in the network.

METHODS

Intracranial recordings from 10 medically refractory MTLE patients were analyzed. In each patient, hour-long recordings were selected for having frequent interictal discharges and no ictal events. For all possible pairs of electrodes, conditional probability of the occurrence of interictal spikes within a 150-millisecond bin was computed. These probabilities were used to construct a weighted graph between all electrodes, and the node degree was estimated. To assess the relationship of the highly connected regions in this network to the clinically identified seizure network, logistic regression was used to model the regions that were surgically resected using weighted node degree and number of spikes in each channel as factors. Lastly, the conditional spike probability was normalized and averaged across patients to visualize the MTLE network at group level.

RESULTS

We generated the first graph of connectivity across a cohort of MTLE patients using interictal activity. The most consistent connections were hippocampus to amygdala, anterior fusiform cortex to hippocampus, and parahippocampal gyrus projections to amygdala. Additionally, the weighted node degree and number of spikes modeled the brain regions identified as seizure networks by clinicians.

SIGNIFICANCE

Apart from identifying interictal measures that can model patient-specific epileptogenic networks, we also produce a group map of network connectivity from a cohort of MTLE patients.