Clinical yield of magnetoencephalography distributed source imaging in epilepsy: A comparison with equivalent current dipole method

Prestimulus functional connectivity reflects attention orientation in a prospective memory task: A magnetoencephalographic (MEG) study

Prospective Memory (PM) is the ability to encode an intention in memory and retrieve it at the right time in the future. After the intention is formed, it must be maintained in memory while simultaneously monitoring the environment until the occurrence of the stimulus associated with its retrieval. Therefore, monitoring and maintenance processes must work in conjunction to subserve PM processing (monitoring/maintenance phase). Several brain regions play a role in PM, such as the anterior prefrontal cortex, inferior parietal lobules, and precuneus. Notably, these regions belong to different brain networks and are differently involved depending on the memory and attentional requests of the PM task. In this study, we investigate the neural bases of PM from a network perspective, using functional connectivity (FC) analysis to identify the networks involved in the attentional and memory mechanisms underlying PM. To this end, we analyzed MEG data collected in two different PM conditions, enhancing either the monitoring (i.e., attention) or the maintenance (i.e., memory) loads of the PM task. To disentangle the neural correlates of these mechanisms from other processes occurring after stimulus presentation, the analysis focused on the prestimulus time window (monitoring/maintenance phase). The monitoring-load condition was characterized by increased inter-network FC of the Dorsal Attention Network (DAN) in the alpha band, a marker of increased top-down monitoring. In contrast, the maintenance-load condition was associated with increased connectivity of the Ventral Attention Network (VAN) with the FrontoParietal Control and the Default-Mode Networks (FPCN and DMN, respectively). Additionally, response times were found to correlate with prestimulus alpha connectivity of different networks in the two conditions. These differences in connectivity within and between networks support the hypothesis that different networks (DAN, or VAN and DMN) and mechanisms (top-down or bottom-up, respectively) are involved in PM processing depending on the features of the PM task.

Reliable measurement of auditory-driven gamma synchrony with a single EEG electrode: A simultaneous EEG-MEG study

OBJECTIVE

Auditory-driven gamma synchrony (GS) is linked to the function of a specific cortical circuit based on a parvalbumin+ and pyramidal neuron loop. This circuit is impaired in neuropsychiatric conditions (i.e. schizophrenia, Alzheimer's disease, stroke etc.) and its relevance in clinical practice is increasingly being recognized. Auditory stimulation at a typical gamma frequency of 40 Hz can be applied as a 'stress test' of excitation/inhibition (E/I) of the entire cerebral cortex, to drive GS and record it with magnetoencephalography (MEG) or high-density electroencephalography (EEG). However, these two techniques are costly and not widely available. Therefore, we assessed whether a single EEG electrode is sufficient to provide an accurate estimate of the auditory-driven GS level of the entire cortical surface while expecting the highest correspondence in the auditory and somatosensory cortices.

METHODS

We measured simultaneous EEG-MEG in 29 healthy subjects, utilizing 3 EEG electrodes (C4, F4, O2) and a full MEG setup. Recordings were performed during binaural exposure to auditory gamma stimulation and during silence. We compared GS measurement of each of the three EEG electrodes separately against full MEG mapping. Time-resolved phase locking value (PLVt) was computed between EEG signals and cortex reconstructed MEG signals.

RESULTS

During auditory stimulation, but not at rest, EEG captures a significant amount of GS, especially from both auditory cortices and motor-premotor regions. This was especially true for frontal (C4) and central electrodes (F4).

DISCUSSION AND CONCLUSIONS

While hd-EEG and MEG are necessary for accurate spatial mapping of GS at rest and during auditory stimulation, a single EEG channel is sufficient to detect the global level of GS. These results have great translational potential for mapping GS in standard clinical settings.

Electromagnetic Source Imaging in Presurgical Evaluation of Children with Drug-Resistant Epilepsy

For children with drug-resistant epilepsy (DRE), seizure freedom relies on the delineation and resection (or ablation/disconnection) of the epileptogenic zone (EZ) while preserving the eloquent brain areas. The development of a reliable and noninvasive localization method that provides clinically useful information for the localization of the EZ is, therefore, crucial to achieving successful surgical outcomes. Electric and magnetic source imaging (ESI and MSI) have been increasingly utilized in the presurgical evaluation of these patients showing promising findings in the delineation of epileptogenic as well as eloquent brain areas. Moreover, the combination of ESI and MSI into a single solution, namely electromagnetic source imaging (EMSI), performed on simultaneous high-density electroencephalography (HD-EEG) and magnetoencephalography (MEG) recordings has shown higher source localization accuracy than either modality alone. Despite these encouraging findings, such techniques are performed in only a few tertiary epilepsy centers, are rarely recorded simultaneously, and are underutilized in pediatric cohorts. This study illustrates the experimental setup for recording simultaneous MEG and HD-EEG data as well as the methodological framework for analyzing these data aiming to localize the irritative zone, the seizure onset zone, and eloquent brain areas in children with DRE. More specifically, the experimental setups are presented for (i) recording and localizing interictal and ictal epileptiform activity during sleep and (ii) recording visual-, motor-, auditory-, and somatosensory-evoked responses and mapping relevant eloquent brain areas (i.e., visual, motor, auditory, and somatosensory) during visuomotor task, as well as auditory and somatosensory stimulations. Detailed steps of the data analysis pipeline are further presented for performing EMSI as well as individual ESI and MSI using equivalent current dipole (ECD) and dynamic statistical parametric mapping (dSPM).

EEG/MEG source imaging of deep brain activity within the maximum entropy on the mean framework: Simulations and validation in epilepsy

Electro/Magneto-EncephaloGraphy (EEG/MEG) source imaging (EMSI) of epileptic activity from deep generators is often challenging due to the higher sensitivity of EEG/MEG to superficial regions and to the spatial configuration of subcortical structures. We previously demonstrated the ability of the coherent Maximum Entropy on the Mean (cMEM) method to accurately localize the superficial cortical generators and their spatial extent. Here, we propose a depth-weighted adaptation of cMEM to localize deep generators more accurately. These methods were evaluated using realistic MEG/high-density EEG (HD-EEG) simulations of epileptic activity and actual MEG/HD-EEG recordings from patients with focal epilepsy. We incorporated depth-weighting within the MEM framework to compensate for its preference for superficial generators. We also included a mesh of both hippocampi, as an additional deep structure in the source model. We generated 5400 realistic simulations of interictal epileptic discharges for MEG and HD-EEG involving a wide range of spatial extents and signal-to-noise ratio (SNR) levels, before investigating EMSI on clinical HD-EEG in 16 patients and MEG in 14 patients. Clinical interictal epileptic discharges were marked by visual inspection. We applied three EMSI methods: cMEM, depth-weighted cMEM and depth-weighted minimum norm estimate (MNE). The ground truth was defined as the true simulated generator or as a drawn region based on clinical information available for patients. For deep sources, depth-weighted cMEM improved the localization when compared to cMEM and depth-weighted MNE, whereas depth-weighted cMEM did not deteriorate localization accuracy for superficial regions. For patients' data, we observed improvement in localization for deep sources, especially for the patients with mesial temporal epilepsy, for which cMEM failed to reconstruct the initial generator in the hippocampus. Depth weighting was more crucial for MEG (gradiometers) than for HD-EEG. Similar findings were found when considering depth weighting for the wavelet extension of MEM. In conclusion, depth-weighted cMEM improved the localization of deep sources without or with minimal deterioration of the localization of the superficial sources. This was demonstrated using extensive simulations with MEG and HD-EEG and clinical MEG and HD-EEG for epilepsy patients.

Excitation/Inhibition balance relates to cognitive function and gene expression in temporal lobe epilepsy: a high density EEG assessment with aperiodic exponent

Patients with epilepsy are characterized by a dysregulation of excitation/inhibition balance (E/I). The assessment of E/I may inform clinicians during the diagnosis and therapy management, even though it is rarely performed. An accessible measure of the E/I of the brain represents a clinically relevant feature. Here, we exploited the exponent of the aperiodic component of the power spectrum of the electroencephalography (EEG) signal, as a non-invasive and cost-effective proxy of the E/I balance. We recorded resting-state activity with high-density EEG from 67 patients with temporal lobe epilepsy and 35 controls. We extracted the exponent of the aperiodic fit of the power spectrum from source-reconstructed EEG and tested differences between patients with epilepsy and controls. Spearman's correlation was performed between the exponent and clinical variables (age of onset, epilepsy duration and neuropsychology) and cortical expression of epilepsy-related genes derived from the Allen Human Brain Atlas. Patients with temporal lobe epilepsy showed a significantly larger exponent, corresponding to inhibition-directed E/I balance, in bilateral frontal and temporal regions. Lower E/I in the left entorhinal and bilateral dorsolateral prefrontal cortices corresponded to a lower performance of short-term verbal memory. Limited to patients with temporal lobe epilepsy, we detected a significant correlation between the exponent and the cortical expression of GABRA1, GRIN2A, GABRD, GABRG2, KCNA2 and PDYN genes. EEG aperiodic exponent maps the E/I balance non-invasively in patients with epilepsy and reveals a close relationship between altered E/I patterns, cognition and genetics.

Assessing cortical excitability with electroencephalography: A pilot study with EEG-iTBS

BACKGROUND

Cortical excitability measures neural reactivity to stimuli, usually delivered via Transcranial Magnetic Stimulation (TMS). Excitation/inhibition balance (E/I) is the ongoing equilibrium between excitatory and inhibitory activity of neural circuits. According to some studies, E/I could be estimated in-vivo and non-invasively through the modeling of electroencephalography (EEG) signals and termed 'intrinsic excitability' measures. Several measures have been proposed (phase consistency in the gamma band, sample entropy, exponent of the power spectral density 1/f curve, E/I index extracted from detrend fluctuation analysis, and alpha power). Intermittent theta burst stimulation (iTBS) of the primary motor cortex (M1) is a non-invasive neuromodulation technique allowing controlled and focal enhancement of TMS cortical excitability and E/I of the stimulated hemisphere.

OBJECTIVE

Investigating to what extent E/I estimates scale with TMS excitability and how they relate to each other.

METHODS

M1 excitability (TMS) and several E/I estimates extracted from resting state EEG recordings were assessed before and after iTBS in a cohort of healthy subjects.

RESULTS

Enhancement of TMS M1 excitability, as measured through motor-evoked potentials (MEPs), and phase consistency of the cortex in high gamma band correlated with each other. Other measures of E/I showed some expected results, but no correlation with TMS excitability measures or strong consistency with each other.

CONCLUSIONS

EEG E/I estimates offer an intriguing opportunity to map cortical excitability non-invasively, with high spatio-temporal resolution and with a stimulus independent approach. While different EEG E/I estimates may reflect the activity of diverse excitatory-inhibitory circuits, spatial phase synchrony in the gamma band is the measure that best captures excitability changes in the primary motor cortex.

Electrical Source Imaging of Somatosensory Evoked Potentials from Intracranial EEG Signals

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

The impact of ROI extraction method for MEG connectivity estimation: practical recommendations for the study of resting state data

Magnetoencephalography and electroencephalography (M/EEG) seed-based connectivity analysis requires the extraction of measures from regions of interest (ROI). M/EEG ROI-derived source activity can be treated in different ways. It is possible, for instance, to average each ROI's time series prior to calculating connectivity measures. Alternatively, one can compute connectivity maps for each element of the ROI prior to dimensionality reduction to obtain a single map. The impact of these different strategies on connectivity results is still unclear. Here, we address this question within a large MEG resting state cohort (N=113) and within simulated data. We consider 68 ROIs (Desikan-Kiliany atlas), two measures of connectivity (phase locking value-PLV, and its imaginary counterpart- ciPLV), and three frequency bands (theta 4-8 Hz, alpha 9-12 Hz, beta 15-30 Hz). We compare four extraction methods: (i) mean, or (ii) PCA of the activity within the seed or ROI before computing connectivity, map of the (iii) average, or (iv) maximum connectivity after computing connectivity for each element of the seed. Hierarchical clustering is then applied to compare connectivity outputs across multiple strategies, followed by direct contrasts across extraction methods. Finally, the results are validated by using a set of realistic simulations. We show that ROI-based connectivity maps vary remarkably across strategies in terms of connectivity magnitude and spatial distribution. Dimensionality reduction procedures conducted after computing connectivity are more similar to each-other, while PCA before approach is the most dissimilar to other approaches. Although differences across methods are consistent across frequency bands, they are influenced by the connectivity metric and ROI size. Greater differences were observed for ciPLV than PLV, and in larger ROIs. Realistic simulations confirmed that after aggregation procedures are generally more accurate but have lower specificity (higher rate of false positive connections). Though computationally demanding, after dimensionality reduction strategies should be preferred when higher sensitivity is desired. Given the remarkable differences across aggregation procedures, caution is warranted in comparing results across studies applying different methods.

An artificial intelligence-based pipeline for automated detection and localisation of epileptic sources from magnetoencephalography

Objective.Magnetoencephalography (MEG) is a powerful non-invasive diagnostic modality for presurgical epilepsy evaluation. However, the clinical utility of MEG mapping for localising epileptic foci is limited by its low efficiency, high labour requirements, and considerable interoperator variability. To address these obstacles, we proposed a novel artificial intelligence-based automated magnetic source imaging (AMSI) pipeline for automated detection and localisation of epileptic sources from MEG data.Approach.To expedite the analysis of clinical MEG data from patients with epilepsy and reduce human bias, we developed an autolabelling method, a deep-learning model based on convolutional neural networks and a hierarchical clustering method based on a perceptual hash algorithm, to enable the coregistration of MEG and magnetic resonance imaging, the detection and clustering of epileptic activity, and the localisation of epileptic sources in a highly automated manner. We tested the capability of the AMSI pipeline by assessing MEG data from 48 epilepsy patients.Main results.The AMSI pipeline was able to rapidly detect interictal epileptiform discharges with 93.31% ± 3.87% precision based on a 35-patient dataset (with sevenfold patientwise cross-validation) and robustly rendered accurate localisation of epileptic activity with a lobar concordance of 87.18% against interictal and ictal stereo-electroencephalography findings in a 13-patient dataset. We also showed that the AMSI pipeline accomplishes the necessary processes and delivers objective results within a much shorter time frame (∼12 min) than traditional manual processes (∼4 h).Significance.The AMSI pipeline promises to facilitate increased utilisation of MEG data in the clinical analysis of patients with epilepsy.

Validating MEG source imaging of resting state oscillatory patterns with an intracranial EEG atlas

BACKGROUND

Magnetoencephalography (MEG) is a widely used non-invasive tool to estimate brain activity with high temporal resolution. However, due to the ill-posed nature of the MEG source imaging (MSI) problem, the ability of MSI to identify accurately underlying brain sources along the cortical surface is still uncertain and requires validation.

METHOD

We validated the ability of MSI to estimate the background resting state activity of 45 healthy participants by comparing it to the intracranial EEG (iEEG) atlas (https://mni-open-ieegatlas.

RESEARCH

mcgill.ca/). First, we applied wavelet-based Maximum Entropy on the Mean (wMEM) as an MSI technique. Next, we converted MEG source maps into intracranial space by applying a forward model to the MEG-reconstructed source maps, and estimated virtual iEEG (ViEEG) potentials on each iEEG channel location; we finally quantitatively compared those with actual iEEG signals from the atlas for 38 regions of interest in the canonical frequency bands.

RESULTS

The MEG spectra were more accurately estimated in the lateral regions compared to the medial regions. The regions with higher amplitude in the ViEEG than in the iEEG were more accurately recovered. In the deep regions, MEG-estimated amplitudes were largely underestimated and the spectra were poorly recovered. Overall, our wMEM results were similar to those obtained with minimum norm or beamformer source localization. Moreover, the MEG largely overestimated oscillatory peaks in the alpha band, especially in the anterior and deep regions. This is possibly due to higher phase synchronization of alpha oscillations over extended regions, exceeding the spatial sensitivity of iEEG but detected by MEG. Importantly, we found that MEG-estimated spectra were more comparable to spectra from the iEEG atlas after the aperiodic components were removed.

CONCLUSION

This study identifies brain regions and frequencies for which MEG source analysis is likely to be reliable, a promising step towards resolving the uncertainty in recovering intracerebral activity from non-invasive MEG studies.

Neurophysiology, Neuropsychology, and Epilepsy, in 2022: Hills We Have Climbed and Hills Ahead. Neurophysiology in epilepsy

Since the discovery of the human electroencephalogram (EEG), neurophysiology techniques have become indispensable tools in our armamentarium to localize epileptic seizures. New signal analysis techniques and the prospects of artificial intelligence and big data will offer unprecedented opportunities to further advance the field in the near future, ultimately resulting in improved quality of life for many patients with drug-resistant epilepsy. This article summarizes selected presentations from Day 1 of the two-day symposium "Neurophysiology, Neuropsychology, Epilepsy, 2022: Hills We Have Climbed and the Hills Ahead". Day 1 was dedicated to highlighting and honoring the work of Dr. Jean Gotman, a pioneer in EEG, intracranial EEG, simultaneous EEG/ functional magnetic resonance imaging, and signal analysis of epilepsy. The program focused on two main research directions of Dr. Gotman, and was dedicated to "High-frequency oscillations, a new biomarker of epilepsy" and "Probing the epileptic focus from inside and outside". All talks were presented by colleagues and former trainees of Dr. Gotman. The extended summaries provide an overview of historical and current work in the neurophysiology of epilepsy with emphasis on novel EEG biomarkers of epilepsy and source imaging and concluded with an outlook on the future of epilepsy research, and what is needed to bring the field to the next level.

Data-driven approach for the delineation of the irritative zone in epilepsy in MEG

The reliable identification of the irritative zone (IZ) is a prerequisite for the correct clinical evaluation of medically refractory patients affected by epilepsy. Given the complexity of MEG data, visual analysis of epileptiform neurophysiological activity is highly time consuming and might leave clinically relevant information undetected. We recorded and analyzed the interictal activity from seven patients affected by epilepsy (Vectorview Neuromag), who successfully underwent epilepsy surgery (Engel > = II). We visually marked and localized characteristic epileptiform activity (VIS). We implemented a two-stage pipeline for the detection of interictal spikes and the delineation of the IZ. First, we detected candidate events from peaky ICA components, and then clustered events around spatio-temporal patterns identified by convolutional sparse coding. We used the average of clustered events to create IZ maps computed at the amplitude peak (PEAK), and at the 50% of the peak ascending slope (SLOPE). We validated our approach by computing the distance of the estimated IZ (VIS, SLOPE and PEAK) from the border of the surgically resected area (RA). We identified 25 spatiotemporal patterns mimicking the underlying interictal activity (3.6 clusters/patient). Each cluster was populated on average by 22.1 [15.0–31.0] spikes. The predicted IZ maps had an average distance from the resection margin of 8.4 ± 9.3 mm for visual analysis, 12.0 ± 16.5 mm for SLOPE and 22.7 ±. 16.4 mm for PEAK. The consideration of the source spread at the ascending slope provided an IZ closer to RA and resembled the analysis of an expert observer. We validated here the performance of a data-driven approach for the automated detection of interictal spikes and delineation of the IZ. This computational framework provides the basis for reproducible and bias-free analysis of MEG recordings in epilepsy.

Hierarchical Bayesian modeling of the relationship between task-related hemodynamic responses and cortical excitability

Investigating the relationship between task-related hemodynamic responses and cortical excitability is challenging because it requires simultaneous measurement of hemodynamic responses while applying noninvasive brain stimulation. Moreover, cortical excitability and task-related hemodynamic responses are both associated with inter-/intra-subject variability. To reliably assess such a relationship, we applied hierarchical Bayesian modeling. This study involved 16 healthy subjects who underwent simultaneous Paired Associative Stimulation (PAS10, PAS25, Sham) while monitoring brain activity using functional Near-Infrared Spectroscopy (fNIRS), targeting the primary motor cortex (M1). Cortical excitability was measured by Motor Evoked Potentials (MEPs), and the motor task-related hemodynamic responses were measured using fNIRS 3D reconstructions. We constructed three models to investigate: (1) PAS effects on the M1 excitability, (2) PAS effects on fNIRS hemodynamic responses to a finger tapping task, and (3) the correlation between PAS effects on M1 excitability and PAS effects on task-related hemodynamic responses. Significant increase in cortical excitability was found following PAS25, whereas a small reduction of the cortical excitability was shown after PAS10 and a subtle increase occurred after sham. Both HbO and HbR absolute amplitudes increased after PAS25 and decreased after PAS10. The probability of the positive correlation between modulation of cortical excitability and hemodynamic activity was 0.77 for HbO and 0.79 for HbR. We demonstrated that PAS stimulation modulates task-related cortical hemodynamic responses in addition to M1 excitability. Moreover, the positive correlation between PAS modulations of excitability and hemodynamics brought insight into understanding the fundamental properties of cortical function and cortical excitability.

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.

Presurgical accuracy of dipole clustering in MRI-negative pediatric patients with epilepsy: Validation against intracranial EEG and resection

OBJECTIVE

To assess the utility of interictal magnetic and electric source imaging (MSI and ESI) using dipole clustering in magnetic resonance imaging (MRI)-negative patients with drug resistant epilepsy (DRE).

METHODS

We localized spikes in low-density (LD-EEG) and high-density (HD-EEG) electroencephalography as well as magnetoencephalography (MEG) recordings using dipoles from 11 pediatric patients. We computed each dipole's level of clustering and used it to discriminate between clustered and scattered dipoles. For each dipole, we computed the distance from seizure onset zone (SOZ) and irritative zone (IZ) defined by intracranial EEG. Finally, we assessed whether dipoles proximity to resection was predictive of outcome.

RESULTS

LD-EEG had lower clusterness compared to HD-EEG and MEG (p < 0.05). For all modalities, clustered dipoles showed higher proximity to SOZ and IZ than scattered (p < 0.001). Resection percentage was higher in optimal vs. suboptimal outcome patients (p < 0.001); their proximity to resection was correlated to outcome (p < 0.001). No difference in resection percentage was seen for scattered dipoles between groups.

CONCLUSION

MSI and ESI dipole clustering helps to localize the SOZ and IZ and facilitate the prognostic assessment of MRI-negative patients with DRE.

SIGNIFICANCE

Assessing the MSI and ESI clustering allows recognizing epileptogenic areas whose removal is associated with optimal outcome.

Clinical Yield of Electromagnetic Source Imaging and Hemodynamic Responses in Epilepsy: Validation With Intracerebral Data

BACKGROUND AND OBJECTIVES

Accurate delineation of the seizure-onset zone (SOZ) in focal drug-resistant epilepsy often requires stereo-EEG (SEEG) recordings. Our aims were to propose a truly objective and quantitative comparison between EEG/magnetoencephalography (MEG) source imaging (EMSI), EEG/fMRI responses for similar spikes with primary irritative zone (PIZ) and SOZ defined by SEEG and to evaluate the value of EMSI and EEG/fMRI to predict postsurgical outcome.

METHODS

We identified patients with drug-resistant epilepsy who underwent EEG/MEG, EEG/fMRI, and subsequent SEEG at the Epilepsy Service from the Montreal Neurological Institute and Hospital. We quantified multimodal concordance within the SEEG channel space as spatial overlap with PIZ/SOZ and distances to the spike-onset, spike maximum amplitude and seizure core intracerebral channels by applying a new methodology consisting of converting EMSI results into SEEG electrical potentials (EMSIe-SEEG) and projecting the most significant fMRI response on the SEEG channels (fMRIp-SEEG). Spatial overlaps with PIZ/SOZ (AUCPIZ, AUCSOZ) were assessed by using the area under the receiver operating characteristic curve (AUC). Here, AUC represents the probability that a randomly picked active contact exhibited higher amplitude when located inside the spatial reference than outside.

RESULTS

Seventeen patients were included. Mean spatial overlaps with the PIZ and SOZ were 0.71 and 0.65 for EMSIe-SEEG and 0.57 and 0.62 for fMRIp-SEEG. Good EMSIe-SEEG spatial overlap with the PIZ was associated with smaller distance from the maximum EMSIe-SEEG contact to the spike maximum amplitude channel (median distance 14 mm). Conversely, good fMRIp-SEEG spatial overlap with the SOZ was associated with smaller distances from the maximum fMRIp-SEEG contact to the spike-onset and seizure core channels (median distances 10 and 5 mm, respectively). Surgical outcomes were correctly predicted by EEG/MEG in 12 of 15 (80%) patients and EEG/fMRI in 6 of 11(54%) patients.

DISCUSSION

With the use of a unique quantitative approach estimating EMSI and fMRI results in the reference SEEG channel space, EEG/MEG and EEG/fMRI accurately localized the SOZ and the PIZ. Precisely, EEG/MEG more accurately localized the PIZ, whereas EEG/fMRI was more sensitive to the SOZ. Both neuroimaging techniques provide complementary localization that can help guide SEEG implantation and select good candidates for surgery.

Novel User-Friendly Application for MRI Segmentation of Brain Resection following Epilepsy Surgery

Delineation of resected brain cavities on magnetic resonance images (MRIs) of epilepsy surgery patients is essential for neuroimaging/neurophysiology studies investigating biomarkers of the epileptogenic zone. The gold standard to delineate the resection on MRI remains manual slice-by-slice tracing by experts. Here, we proposed and validated a semiautomated MRI segmentation pipeline, generating an accurate model of the resection and its anatomical labeling, and developed a graphical user interface (GUI) for user-friendly usage. We retrieved pre- and postoperative MRIs from 35 patients who had focal epilepsy surgery, implemented a region-growing algorithm to delineate the resection on postoperative MRIs and tested its performance while varying different tuning parameters. Similarity between our output and hand-drawn gold standards was evaluated via dice similarity coefficient (DSC; range: 0-1). Additionally, the best segmentation pipeline was trained to provide an automated anatomical report of the resection (based on presurgical brain atlas). We found that the best-performing set of parameters presented DSC of 0.83 (0.72-0.85), high robustness to seed-selection variability and anatomical accuracy of 90% to the clinical postoperative MRI report. We presented a novel user-friendly open-source GUI that implements a semiautomated segmentation pipeline specifically optimized to generate resection models and their anatomical reports from epilepsy surgery patients, while minimizing user interaction.

Sensorimotor integration within the primary motor cortex by selective nerve fascicle stimulation

KEY POINTS

Cortical integration of sensory inputs is crucial for dexterous movement. Short-latency somatosensory afferent inhibition of motor cortical output is typically produced by peripheral whole-nerve stimulation. We exploited intraneural multichannel electrodes used to provide sensory feedback for prosthesis control to assess whether and how selective intraneural sensory stimulation affects sensorimotor cortical circuits in humans. The activation of the primary somatosensory cortex (S1) was explored by recording scalp somatosensory evoked potentials. Sensorimotor integration was tested by measuring the inhibitory effect of the afferent stimulation on the output of the primary motor cortex (M1) generated by transcranial magnetic stimulation. We demonstrate in humans that selective intraneural sensory stimulation elicits a measurable activation of S1 and that it inhibits the output of M1 at the same time range of whole-nerve superficial stimulation.

ABSTRACT

The integration of sensory inputs in the motor cortex is crucial for dexterous movement. We recently demonstrated that a closed-loop control based on the feedback provided through intraneural multi-channel electrodes implanted in the median and ulnar nerves of a participant with upper limb amputation improved manipulation skills and increased prosthesis embodiment. Here we assessed, in the same participant, whether and how selective intraneural sensory stimulation also elicits a measurable cortical activation and affects sensorimotor cortical circuits. After estimating the activation of the primary somatosensory cortex evoked by intraneural stimulation, sensorimotor integration was investigated by testing the inhibition of primary motor cortex (M1) output to transcranial magnetic stimulation, after both intraneural and perineural stimulation. Selective sensory intraneural stimulation evoked a low-amplitude, 16 ms-latency, parietal response in the same area of the earliest component evoked by whole-nerve stimulation, compatible with fast-conducting afferent fiber activation. For the first time, we show that the same intraneural stimulation was also capable of decreasing M1 output, at the same time range of the short-latency afferent inhibition effect of whole-nerve superficial stimulation. The inhibition generated by the stimulation of channels activating only sensory fibers was stronger than the one due to intraneural or perineural stimulation of channels activating mixed fibers. We demonstrate in a human subject that the cortical sensorimotor integration inhibiting M1 output previously described after the experimental whole-nerve stimulation is present also with a more ecological selective sensory fiber stimulation. Abstract Figure: Double-sided filament electrodes (ds-FILE), bearing 16 active sites, and perineural Cuff electrodes were implanted in the median and ulnar nerve of the arm in a hand amputee (upper left panel, single nerve represented). Selectivity of stimulation (1), evoked activity in the somatosensory cortex (2), and sensorimotor integration (3) were investigated. TMS: transcranial magnetic stimulation. This article is protected by copyright. All rights reserved.

Resting state network connectivity is attenuated by fMRI acoustic noise

INTRODUCTION

During the past decades there has been an increasing interest in tracking brain network fluctuations in health and disease by means of resting state functional magnetic resonance imaging (rs-fMRI). Rs-fMRI however does not provide the ideal environmental setting, as participants are continuously exposed to noise generated by MRI coils during acquisition of Echo Planar Imaging (EPI). We investigated the effect of EPI noise on resting state activity and connectivity using magnetoencephalography (MEG), by reproducing the acoustic characteristics of rs-fMRI environment during the recordings. As compared to fMRI, MEG has little sensitivity to brain activity generated in deep brain structures, but has the advantage to capture both the dynamic of cortical magnetic oscillations with high temporal resolution and the slow magnetic fluctuations highly correlated with BOLD signal.

METHODS

Thirty healthy subjects were enrolled in a counterbalanced design study including three conditions: a) silent resting state (Silence), b) resting state upon EPI noise (fMRI), and c) resting state upon white noise (White). White noise was employed to test the specificity of fMRI noise effect. The amplitude envelope correlation (AEC) in alpha band measured the connectivity of seven Resting State Networks (RSN) of interest (default mode network, dorsal attention network, language, left and right auditory and left and right sensory-motor). Vigilance dynamic was estimated from power spectral activity.

RESULTS

fMRI and White acoustic noise consistently reduced connectivity of cortical networks. The effects were widespread, but noise and network specificities were also present. For fMRI noise, decreased connectivity was found in the right auditory and sensory-motor networks. Progressive increase of slow theta-delta activity related to drowsiness was found in all conditions, but was significantly higher for fMRI . Theta-delta significantly and positively correlated with variations of cortical connectivity.

DISCUSSION

rs-fMRI connectivity is biased by unavoidable environmental factors during scanning, which warrant more careful control and improved experimental designs. MEG is free from acoustic noise and allows a sensitive estimation of resting state connectivity in cortical areas. Although underutilized, MEG could overcome issues related to noise during fMRI, in particular when investigation of motor and auditory networks is needed.

Imaging the extent and location of spatiotemporally distributed epileptiform sources from MEG measurements

Non-invasive MEG/EEG source imaging provides valuable information about the epileptogenic brain areas which can be used to aid presurgical planning in focal epilepsy patients suffering from drug-resistant seizures. However, the source extent estimation for electrophysiological source imaging remains to be a challenge and is usually largely dependent on subjective choice. Our recently developed algorithm, fast spatiotemporal iteratively reweighted edge sparsity minimization (FAST-IRES) strategy, has been shown to objectively estimate extended sources from EEG recording, while it has not been applied to MEG recordings. In this work, through extensive numerical experiments and real data analysis in a group of focal drug-resistant epilepsy patients' interictal spikes, we demonstrated the ability of FAST-IRES algorithm to image the location and extent of underlying epilepsy sources from MEG measurements. Our results indicate the merits of FAST-IRES in imaging the location and extent of epilepsy sources for pre-surgical evaluation from MEG measurements.

Cortical Thickness of Brain Areas Beyond Stroke Lesions and Sensory-Motor Recovery: A Systematic Review

Background: The clinical outcome of patients suffering from stroke is dependent on multiple factors. The features of the lesion itself play an important role but clinical recovery is remarkably influenced by the plasticity mechanisms triggered by the stroke and occurring at a distance from the lesion. The latter translate into functional and structural changes of which cortical thickness might be easy to quantify one of the main players. However, studies on the changes of cortical thickness in brain areas beyond stroke lesion and their relationship to sensory-motor recovery are sparse. Objectives: To evaluate the effects of cerebral stroke on cortical thickness (CT) beyond the stroke lesion and its association with sensory-motor recovery. Materials and Methods: Five electronic databases (PubMed, Embase, Web of Science, Scopus and the Cochrane Library) were searched. Methodological quality of the included studies was assessed with the Newcastle-Ottawa Scale for non-randomized controlled trials and the Risk of Bias Cochrane tool for randomized controlled trials. Results: The search strategy retrieved 821 records, 12 studies were included and risk of bias assessed. In most of the included studies, cortical thinning was seen at the ipsilesional motor area (M1). Cortical thinning can occur beyond the stroke lesion, typically in regions anatomically connected because of anterograde degeneration. Nonetheless, studies also reported cortical thickening of regions of the unaffected hemisphere, likely related to compensatory plasticity. Some studies revealed a significant correlation between changes in cortical thickness of M1 or somatosensory (S1) cortical areas and motor function recovery. Discussion and Conclusions: Following a stroke, changes in cortical thickness occur both in regions directly connected to the stroke lesion and in contralateral hemisphere areas as well as in the cerebellum. The underlying mechanisms leading to these changes in cortical thickness are still to be fully understood and further research in the field is needed. Systematic Review Registration: https://www.crd.york.ac.uk/prospero/display_record.php?ID=CRD42020200539; PROSPERO 2020, identifier: CRD42020200539.

How cerebral cortex protects itself from interictal spikes: The alpha/beta inhibition mechanism

Interactions between interictal epileptiform discharges (IEDs) and distant cortical regions subserve potential effects on cognition of patients with focal epilepsy. We hypothesize that "healthy" brain areas at a distance from the epileptic focus may respond to the interference of IEDs by generating inhibitory alpha and beta oscillations. We predict that more prominent alpha-beta oscillations can be found in patients with less impaired neurocognitive profile. We performed a source imaging magnetoencephalography study, including 41 focal epilepsy patients: 21 with frontal lobe epilepsy (FLE) and 20 with mesial temporal lobe epilepsy. We investigated the effect of anterior (i.e., frontal and temporal) IEDs on the oscillatory pattern over posterior head regions. We compared cortical oscillations (5-80 Hz) temporally linked to 3,749 IEDs (1,945 frontal and 1,803 temporal) versus an equal number of IED-free segments. We correlated results from IED triggered oscillations to global neurocognitive performance. Only frontal IEDs triggered alpha-beta oscillations over posterior head regions. IEDs with higher amplitude triggered alpha-beta oscillations of higher magnitude. The intensity of posterior head region alpha-beta oscillations significantly correlated with a better neuropsychological profile. Our study demonstrated that cerebral cortex protects itself from IEDs with generation of inhibitory alpha-beta oscillations at distant cortical regions. The association of more prominent oscillations with a better cognitive status suggests that this mechanism might play a role in determining the cognitive resilience in patients with FLE.

Clinical Validation of the Champagne Algorithm for Evoked Response Source Localization in Magnetoencephalography

Magnetoencephalography (MEG) is a robust method for non-invasive functional brain mapping of sensory cortices due to its exceptional spatial and temporal resolution. The clinical standard for MEG source localization of functional landmarks from sensory evoked responses is the equivalent current dipole (ECD) localization algorithm, known to be sensitive to initialization, noise, and manual choice of the number of dipoles. Recently many automated and robust algorithms have been developed, including the Champagne algorithm, an empirical Bayesian algorithm, with powerful abilities for MEG source reconstruction and time course estimation (Wipf et al. 2010; Owen et al. 2012). Here, we evaluate automated Champagne performance in a clinical population of tumor patients where there was minimal failure in localizing sensory evoked responses using the clinical standard, ECD localization algorithm. MEG data of auditory evoked potentials and somatosensory evoked potentials from 21 brain tumor patients were analyzed using Champagne, and these results were compared with equivalent current dipole (ECD) fit. Across both somatosensory and auditory evoked field localization, we found there was a strong agreement between Champagne and ECD localizations in all cases. Given resolution of 8mm voxel size, peak source localizations from Champagne were below 10mm of ECD peak source localization. The Champagne algorithm provides a robust and automated alternative to manual ECD fits for clinical localization of sensory evoked potentials and can contribute to improved clinical MEG data processing workflows.

Clinical Validation of the Champagne Algorithm for Epilepsy Spike Localization

Magnetoencephalography (MEG) is increasingly used for presurgical planning in people with medically refractory focal epilepsy. Localization of interictal epileptiform activity, a surrogate for the seizure onset zone whose removal may prevent seizures, is challenging and depends on the use of multiple complementary techniques. Accurate and reliable localization of epileptiform activity from spontaneous MEG data has been an elusive goal. One approach toward this goal is to use a novel Bayesian inference algorithm-the Champagne algorithm with noise learning-which has shown tremendous success in source reconstruction, especially for focal brain sources. In this study, we localized sources of manually identified MEG spikes using the Champagne algorithm in a cohort of 16 patients with medically refractory epilepsy collected in two consecutive series. To evaluate the reliability of this approach, we compared the performance to equivalent current dipole (ECD) modeling, a conventional source localization technique that is commonly used in clinical practice. Results suggest that Champagne may be a robust, automated, alternative to manual parametric dipole fitting methods for localization of interictal MEG spikes, in addition to its previously described clinical and research applications.

Noninvasive Mapping of Ripple Onset Predicts Outcome in Epilepsy Surgery

OBJECTIVE

Intracranial electroencephalographic (icEEG) studies show that interictal ripples propagate across the brain of children with medically refractory epilepsy (MRE), and the onset of this propagation (ripple onset zone [ROZ]) estimates the epileptogenic zone. It is still unknown whether we can map this propagation noninvasively. The goal of this study is to map ripples (ripple zone [RZ]) and their propagation onset (ROZ) using high-density EEG (HD-EEG) and magnetoencephalography (MEG), and to estimate their prognostic value in pediatric epilepsy surgery.

METHODS

We retrospectively analyzed simultaneous HD-EEG and MEG data from 28 children with MRE who underwent icEEG and epilepsy surgery. Using electric and magnetic source imaging, we estimated virtual sensors (VSs) at brain locations that matched the icEEG implantation. We detected ripples on VSs, defined the virtual RZ and virtual ROZ, and estimated their distance from icEEG. We assessed the predictive value of resecting virtual RZ and virtual ROZ for postsurgical outcome. Interictal spike localization on HD-EEG and MEG was also performed and compared with ripples.

RESULTS

We mapped ripple propagation in all patients with HD-EEG and in 27 (96%) patients with MEG. The distance from icEEG did not differ between HD-EEG and MEG when mapping the RZ (26-27mm, p = 0.6) or ROZ (22-24mm, p = 0.4). Resecting the virtual ROZ, but not virtual RZ or the sources of spikes, was associated with good outcome for HD-EEG (p = 0.016) and MEG (p = 0.047).

INTERPRETATION

HD-EEG and MEG can map interictal ripples and their propagation onset (virtual ROZ). Noninvasively mapping the ripple onset may augment epilepsy surgery planning and improve surgical outcome of children with MRE. ANN NEUROL 2021;89:911-925.

Source imaging of seizure onset predicts surgical outcome in pediatric epilepsy

Objective:

To assess whether ictal electric source imaging (ESI) on low-density scalp EEG can approximate the seizure onset zone (SOZ) location and predict surgical outcome in children with refractory epilepsy undergoing surgery.

Methods:

We examined 35 children with refractory epilepsy. We dichotomized surgical outcome into seizure- and non-seizure-free. We identified ictal onsets recorded with scalp and intracranial EEG and localized them using equivalent current dipoles and standardized low-resolution magnetic tomography (sLORETA). We estimated the localization accuracy of scalp EEG as distance of scalp dipoles from intracranial dipoles. We also calculated the distances of scalp dipoles from resection, as well as their resection percentage and compared between seizure-free and non-seizure-free patients. We built receiver operating characteristic curves to test whether resection percentage predicted outcome.

Results:

Resection distance was lower in seizure-free patients for both dipoles (p = 0.006) and sLORETA (p = 0.04). Resection percentage predicted outcome with a sensitivity of 57.1% (95% CI, 34–78.2%), a specificity of 85.7% (95% CI, 57.2–98.2%) and an accuracy of 68.6% (95% CI, 50.7–83.5%) (p = 0.01).

Conclusion:

Ictal ESI performed on low-density scalp EEG can delineate the SOZ and predict outcome.

Significance:

Such an application may increase the number of children who are referred for epilepsy surgery and improve their outcome.

Ant Colony System Optimization for Spatiotemporal Modelling of Combined EEG and MEG Data

Electroencephalography/Magnetoencephalography (EEG/MEG) source localization involves the estimation of neural activity inside the brain volume that underlies the EEG/MEG measures observed at the sensor array. In this paper, we consider a Bayesian finite spatial mixture model for source reconstruction and implement Ant Colony System (ACS) optimization coupled with Iterated Conditional Modes (ICM) for computing estimates of the neural source activity. Our approach is evaluated using simulation studies and a real data application in which we implement a nonparametric bootstrap for interval estimation. We demonstrate improved performance of the ACS-ICM algorithm as compared to existing methodology for the same spatiotemporal model.

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.

Redefining the role of Magnetoencephalography in refractory epilepsy

Magnetoencephalography (MEG) possesses a number of features, including excellent spatiotemporal resolution, that lend itself to the functional imaging of epileptic activity. However its current use is restricted to specific scenarios, namely in the diagnosis refractory focal epilepsies where electroencephalography (EEG) has been inconclusive. This review highlights the recent progress of MEG within epilepsy, including advances in the technique itself such as simultaneous EEG/MEG and intracranial EEG/MEG recording and room temperature MEG recording using optically pumped magnetometers, as well as improved post processing of the data during interictal and ictal activity for accurate source localisation of the epileptogenic focus. These advances should broaden the scope of MEG as an important part of epilepsy diagnostics in the future.

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

•

The irritative zone differs on Electrical Source Imaging (ESI) and EEG-fMRI.

•

Findings differ in functional connectivity and show low temporal correlation.

•

ESI and EEG-fMRI reveal distinct aspects of the irritative zone.

•

The differences may reflect differences in temporal resolution of the techniques.

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.

Enhanced Fast-VESTAL for Magnetoencephalography Source Imaging: From Theory to Clinical Application in Epilepsy

A novel magnetoencephalography source imaging approach called Fast Vector-based Spatio-Temporal Analysis (Fast-VESTAL) has been successfully applied in creating source images from evoked and resting-state data from both healthy subjects and individuals with neurological and/or psychiatric disorders, but its reconstructed source images may show false-positive activations, especially under low signal-to-noise ratio conditions. Here, to effectively reduce false-positive artifacts, we introduced an enhanced Fast-VESTAL (eFast-VESTAL) approach that adopts generalized second-order cone programming. We compared the spatiotemporal characteristics of the eFast-VESTAL approach to those of the popular distributed source approaches (e.g., the minimum L2-norm/ mixed-norm methods) using computer simulations and auditory experiments. More importantly, we applied eFast-VESTAL to the presurgical evaluation of epilepsy. Our results demonstrated that eFast-VESTAL exhibited a lower dipole localization error and/or a higher correlation coefficient (CC) between the estimated source time series and ground truth under various conditions of source waveforms. Experimentally, eFast-VESTAL displayed more focal activation maps and a higher CC between the raw and predicted sensor data in response to auditory stimulation. Notably, eFast-VESTAL was the most accurate method for noninvasively detecting the epileptic zones determined using more invasive stereo-electroencephalography in the comparison.

Pearl and pitfalls in brain functional analysis by event-related potentials: a narrative review by the Italian Psychophysiology and Cognitive Neuroscience Society on methodological limits and clinical reliability-part II

This review focuses on new and/or less standardized event-related potentials methods, in order to improve their knowledge for future clinical applications. The olfactory event-related potentials (OERPs) assess the olfactory functions in time domain, with potential utility in anosmia and degenerative diseases. The transcranial magnetic stimulation-electroencephalography (TMS-EEG) could support the investigation of the intracerebral connections with very high temporal discrimination. Its application in the diagnosis of disorders of consciousness has achieved recent confirmation. Magnetoencephalography (MEG) and event-related fields (ERF) could improve spatial accuracy of scalp signals, with potential large application in pre-surgical study of epileptic patients. Although these techniques have methodological limits, such as high inter- and intraindividual variability and high costs, their diffusion among researchers and clinicians is hopeful, pending their standardization.

Magnetoencephalography resting state connectivity patterns as indicatives of surgical outcome in epilepsy patients

OBJECTIVE

Focal epilepsy is a disorder affecting several brain networks; however, epilepsy surgery usually targets a restricted region, the so-called epileptic focus. There is a growing interest in embedding resting state (RS) connectivity analysis into pre-surgical workup.

APPROACH

In this retrospective study, we analyzed Magnetoencephalography (MEG) long-range RS functional connectivity patterns in patients with drug-resistant focal epilepsy. MEG recorded prior to surgery from seven seizure-free (Engel Ia) and five non seizure-free (Engel III or IV) patients were analyzed (minimum 2-years post-surgical follow-up). MEG segments without any detectable epileptic activity were source localized using wavelet-based Maximum Entropy on the Mean method. Amplitude envelope correlation in the theta (4-8 Hz), alpha (8-13 Hz), and beta (13-26 Hz) bands were used for assessing connectivity.

MAIN RESULTS

For seizure-free patients, we found an isolated epileptic network characterized by weaker connections between the brain region where interictal epileptic discharges (IED) are generated and the rest of the cortex, when compared to connectivity between the corresponding contralateral homologous region and the rest of the cortex. Contrarily, non seizure-free patients exhibited a widespread RS epileptic network characterized by stronger connectivity between the IED generator and the rest of the cortex, in comparison to the contralateral region and the cortex. Differences between the two seizure outcome groups concerned mainly distant long-range connections and were found in the alpha-band.

SIGNIFICANCE

Importantly, these connectivity patterns suggest specific mechanisms describing the underlying organization of the epileptic network and were detectable at the individual patient level, supporting the prospect use of MEG connectivity patterns in epilepsy to predict post-surgical seizure outcome.

Effects of Independent Component Analysis on Magnetoencephalography Source Localization in Pre-surgical Frontal Lobe Epilepsy Patients

Objective: Magnetoencephalography source imaging (MSI) of interictal epileptiform discharges (IED) is a useful presurgical tool in the evaluation of drug-resistant frontal lobe epilepsy (FLE) patients. Yet, failures in MSI can arise related to artifacts and to interference of background activity. Independent component analysis (ICA) is a popular denoising procedure but its clinical application remains challenging, as the selection of multiple independent components (IC) is controversial, operator dependent, and time consuming. We evaluated whether selecting only one IC of interest based on its similarity with the average IED field improves MSI in FLE.

Methods: MSI was performed with the equivalent current dipole (ECD) technique and two distributed magnetic source imaging (dMSI) approaches: minimum norm estimate (MNE) and coherent Maximum Entropy on the Mean (cMEM). MSI accuracy was evaluated under three conditions: (1) ICA of continuous data (Cont_ICA), (2) ICA at the time of IED (IED_ICA), and (3) without ICA (No_ICA). Localization performance was quantitatively measured as actual distance of the source maximum in relation to the focus (Dmin), and spatial dispersion (SD) for dMSI.

Results: After ICA, ECD Dmin did not change significantly (p > 0.200). For both dMSI techniques, ICA application worsened the source localization accuracy. We observed a worsening of both MNE Dmin (p < 0.05, consistently) and MNE SD (p < 0.001, consistently) for both ICA approaches. A similar behaviour was observed for cMEM, for which, however, Cont_ICA seemed less detrimental.

Conclusion: We demonstrated that a simplified ICA approach selecting one IC of interest in combination with distributed magnetic source imaging can be detrimental. More complex approaches may provide better results but would be rather difficult to apply in real-world clinical setting. In a broader perspective, caution should be taken in applying ICA for source localization of interictal activity. To ensure optimal and useful results, effort should focus on acquiring good quality data, minimizing artifacts, and determining optimal candidacy for MEG, rather than counting on data cleaning techniques.

Accuracy and spatial properties of distributed magnetic source imaging techniques in the investigation of focal epilepsy patients

Source localization of interictal epileptiform discharges (IEDs) is clinically useful in the presurgical workup of epilepsy patients. We aimed to compare the performance of four different distributed magnetic source imaging (dMSI) approaches: Minimum norm estimate (MNE), dynamic statistical parametric mapping (dSPM), standardized low-resolution electromagnetic tomography (sLORETA), and coherent maximum entropy on the mean (cMEM). We also evaluated whether a simple average of maps obtained from multiple inverse solutions (Ave) can improve localization accuracy. We analyzed dMSI of 206 IEDs derived from magnetoencephalography recordings in 28 focal epilepsy patients who had a well-defined focus determined through intracranial EEG (iEEG), epileptogenic MRI lesions or surgical resection. dMSI accuracy and spatial properties were quantitatively estimated as: (a) distance from the epilepsy focus, (b) reproducibility, (c) spatial dispersion (SD), (d) map extension, and (e) effect of thresholding on map properties. Clinical performance was excellent for all methods (median distance from the focus MNE = 2.4 mm; sLORETA = 3.5 mm; cMEM = 3.5 mm; dSPM = 6.8 mm, Ave = 0 mm). Ave showed the lowest distance between the map maximum and epilepsy focus (Dmin lower than cMEM, MNE, and dSPM, p = .021, p = .008, p < .001, respectively). cMEM showed the best spatial features, with lowest SD outside the focus (SD lower than all other methods, p < .001 consistently) and high contrast between the generator and surrounding regions. The average map Ave provided the best localization accuracy, whereas cMEM exhibited the lowest amount of spurious distant activity. dMSI techniques have the potential to significantly improve identification of iEEG targets and to guide surgical planning, especially when multiple methods are combined.

Distributed source analysis of magnetoencephalography using a volume head model combined with statistical methods improves focus diagnosis in epilepsy surgery

Deep-seated epileptic focus estimation using magnetoencephalography is challenging because of its low signal-to-noise ratio and the ambiguity of current sources estimated by interictal epileptiform discharge (IED). We developed a distributed source (DS) analysis method using a volume head model as the source space of the forward model and standardized low-resolution brain electromagnetic tomography combined with statistical methods (permutation tests between IEDs and baselines and false discovery rate between voxels to reduce variation). We aimed to evaluate the efficacy of the combined DS (cDS) analysis in surgical cases. In total, 19 surgical cases with adult and pediatric focal epilepsy were evaluated. Both cDS and equivalent current dipole (ECD) analyses were performed in all cases. The concordance rates of the two methods with surgically identified epileptic foci were calculated and compared with surgical outcomes. Concordance rates from the cDS analysis were significantly higher than those from the ECD analysis (68.4% vs. 26.3%), especially in cases with deep-seated lesions, such as in the interhemispheric, fronto-temporal base, and mesial temporal structures (81.8% vs. 9.1%). Furthermore, the concordance rate correlated well with surgical outcomes. In conclusion, cDS analysis has better diagnostic performance in focal epilepsy, especially with deep-seated epileptic focus, and potentially leads to good surgical outcomes.

Scalp ripples as prognostic biomarkers of epileptogenicity in pediatric surgery

Objective

To assess the ability of high‐density Electroencephalography (HD‐EEG) and magnetoencephalography (MEG) to localize interictal ripples, distinguish between ripples co‐occurring with spikes (ripples‐on‐spike) and independent from spikes (ripples‐alone), and evaluate their localizing value as biomarkers of epileptogenicity in children with medically refractory epilepsy.

Methods

We retrospectively studied 20 children who underwent epilepsy surgery. We identified ripples on HD‐EEG and MEG data, localized their generators, and compared them with intracranial EEG (icEEG) ripples. When ripples and spikes co‐occurred, we performed source imaging distinctly on the data above 80 Hz (to localize ripples) and below 70 Hz (to localize spikes). We assessed whether missed resection of ripple sources predicted poor outcome, separately for ripples‐on‐spikes and ripples‐alone. Similarly, predictive value of spikes was calculated.

Results

We observed scalp ripples in 16 patients (10 good outcome). Ripple sources were highly concordant to the icEEG ripples (HD‐EEG concordance: 79%; MEG: 83%). When ripples and spikes co‐occurred, their sources were spatially distinct in 83‐84% of the cases. Removing the sources of ripples‐on‐spikes predicted good outcome with 90% accuracy for HD‐EEG (P = 0.008) and 86% for MEG (P = 0.044). Conversely, removing ripples‐alone did not predict outcome. Resection of spike sources (generated at the same time as ripples) predicted good outcome for HD‐EEG (P = 0.036; accuracy = 87%), while did not reach significance for MEG (P = 0.1; accuracy = 80%).

Interpretation

HD‐EEG and MEG localize interictal ripples with high precision in children with refractory epilepsy. Scalp ripples‐on‐spikes are prognostic, noninvasive biomarkers of epileptogenicity, since removing their cortical generators predicts good outcome. Conversely, scalp ripples‐alone are most likely generated by non‐epileptogenic areas.

Cortical gamma-synchrony measured with magnetoencephalography is a marker of clinical status and predicts clinical outcome in stroke survivors

•

The outcome of stroke survivors is difficult to anticipate.

•

Gamma synchrony is a reliable measure of brain function and reserve.

•

Gamma synchrony is measured with MEG in stroke survivors undergoing rehab.

•

Auditory-entrained gamma synchrony correlates with clinical status and outcome.

Background

The outcome of stroke survivors is difficult to anticipate. While the extent of the anatomical brain lesion is only poorly correlated with the prognosis, functional measures of cortical synchrony, brain networks and cortical plasticity seem to be good predictors of clinical recovery. In this field, gamma (>30 Hz) cortical synchrony is an ideal marker of brain function, as it plays a crucial role for the integration of information, it is an indirect marker of Glutamate/GABA balance and it directly estimates the reserve of parvalbulin-positive neurons, key players in synaptic plasticity. In this study we measured gamma synchronization driven by external auditory stimulation with magnetoencephalography and tested whether it was predictive of the clinical outcome in stroke survivors undergoing intensive rehabilitation in a tertiary rehabilitation center.

Material and methods

Eleven stroke survivors undergoing intensive rehabilitation were prospectively recruited. Gamma synchrony was measured non-invasively within one month from stroke onset with magnetoencephalography, both at rest and during entrainment with external 40 Hz amplitude modulated binaural sounds. Lesion location and volume were quantitatively assessed through a high-resolution anatomical MRI. Barthel index (BI) and Functional Independence Measure (FIM) scales were measured at the beginning and at the end of the admission to the rehabilitation unit.

Results

The spatial distribution of cortical gamma synchrony was altered, and the physiological right hemispheric dominance observed in healthy controls was attenuated or lost. Entrained gamma synchronization (but not resting state gamma synchrony) showed a very high correlation with the clinical status at both admission and discharge (both BI and FIM). Neither clinical status nor gamma synchrony showed a correlation with lesion volume.

Conclusions

Cortical gamma synchrony related to auditory entrainment can be reliably measured in stroke patients. Gamma synchrony is strongly associated with the clinical outcome of stroke survivors undergoing rehabilitation.

Magnetoencephalography imaging of high frequency oscillations strengthens presurgical localization and outcome prediction

In patients with medically refractory epilepsy, resective surgery is the mainstay of therapy to achieve seizure freedom. However, ∼20-50% of cases have intractable seizures post-surgery due to the imprecise determination of epileptogenic zone. Recent intracranial studies suggest that high frequency oscillations between 80 and 200 Hz could serve as one of the consistent epileptogenicity biomarkers for localization of the epileptogenic zone. However, these high frequency oscillations are not adopted in the clinical setting because of difficult non-invasive detection. Here, we investigated non-invasive detection and localization of high frequency oscillations and its clinical utility in accurate pre-surgical assessment and post-surgical outcome prediction. We prospectively recruited 52 patients with medically refractory epilepsy who underwent standard pre-surgical workup including magnetoencephalography (MEG) followed by resective surgery after determination of the epileptogenic zone. The post-surgical outcome was assessed after 22.14 ± 10.05 months. Interictal epileptic spikes were expertly identified, and interictal epileptic oscillations across the neural activity frequency spectrum from 8 to 200 Hz were localized using adaptive spatial filtering methods. Localization results were compared with epileptogenic zone and resected cortex for congruence assessment and validated against the clinical outcome. The concordance rate of high frequency oscillations sources (80-200 Hz) with the presumed epileptogenic zone and the resected cortex were 75.0% and 78.8%, respectively, which is superior to that of other frequency bands and standard dipole fitting methods. High frequency oscillation sources corresponding with the resected cortex, had the best sensitivity of 78.0%, positive predictive value of 100% and an accuracy of 78.84% to predict the patient's surgical outcome, among all other frequency bands. If high frequency oscillation sources were spatially congruent with resected cortex, patients had an odds ratio of 5.67 and 82.4% probability of achieving a favourable surgical outcome. If high frequency oscillations sources were discordant with the epileptogenic zone or resection area, patient has an odds ratio of 0.18 and only 14.3% probability of achieving good outcome, and mostly tended to have an unfavourable outcome (χ2 = 5.22; P = 0.02; φ = -0.317). In receiver operating characteristic curve analyses, only sources of high-frequency oscillations demonstrated the best sensitivity and specificity profile in determining the patient's surgical outcome with area under the curve of 0.76, whereas other frequency bands indicate a poor predictive performance. Our study is the first non-invasive study to detect high frequency oscillations, address the efficacy of high frequency oscillations over the different neural oscillatory frequencies, localize them and clinically validate them with the post-surgical outcome in patients with medically refractory epilepsy. The evidence presented in the current study supports the fact that HFOs might significantly improve the presurgical assessment, and post-surgical outcome prediction, where it could widely be used in a clinical setting as a non-invasive biomarker.

Transcranial direct current stimulation over the sensory-motor regions inhibits gamma synchrony

Transcranial direct current stimulation (tDCS) is a noninvasive brain stimulation technique able to induce plasticity phenomena. Although tDCS application has been spreading over a variety of neuroscience domains, the mechanisms by which the stimulation acts are largely unknown. We investigated tDCS effects on cortical gamma synchrony, which is a crucial player in cortical function. We performed a randomized, sham-controlled, double-blind study on healthy subjects, combining tDCS and magnetoencephalography. By driving brain activity via 40 Hz auditory stimulation during magnetoencephalography, we experimentally tuned cortical gamma synchrony and measured it before and after bilateral tDCS of the primary sensory-motor hand regions (anode left, cathode right). We demonstrated that the stimulation induces a remarkable decrease of gamma synchrony (13 out of 15 subjects), as measured by gamma phase at 40 Hz. tDCS has strong remote effects, as the cortical region mostly affected was located far away from the stimulation site and covered a large area of the right centro-temporal cortex. No significant differences between stimulations were found for baseline gamma synchrony, as well as early transient auditory responses. This suggests a specific tDCS effect on externally driven gamma synchronization. This study sheds new light on the effect of tDCS on cortical function showing that the net effect of the stimulation on cortical gamma synchronization is an inhibition.

Characterizing hippocampal dynamics with MEG: A systematic review and evidence-based guidelines

The hippocampus, a hub of activity for a variety of important cognitive processes, is a target of increasing interest for researchers and clinicians. Magnetoencephalography (MEG) is an attractive technique for imaging spectro-temporal aspects of function, for example, neural oscillations and network timing, especially in shallow cortical structures. However, the decrease in MEG signal-to-noise ratio as a function of source depth implies that the utility of MEG for investigations of deeper brain structures, including the hippocampus, is less clear. To determine whether MEG can be used to detect and localize activity from the hippocampus, we executed a systematic review of the existing literature and found successful detection of oscillatory neural activity originating in the hippocampus with MEG. Prerequisites are the use of established experimental paradigms, adequate coregistration, forward modeling, analysis methods, optimization of signal-to-noise ratios, and protocol trial designs that maximize contrast for hippocampal activity while minimizing those from other brain regions. While localizing activity to specific sub-structures within the hippocampus has not been achieved, we provide recommendations for improving the reliability of such endeavors.

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.

Eslicarbazepine Acetate Modulates EEG Activity and Connectivity in Focal Epilepsy

Introduction: Eslicarbazepine acetate (ESL) is an antiepileptic drug approved as monotherapy or add-on for the treatment of epilepsy with seizures of focal onset. ESL owns a good profile in terms of efficacy and tolerability, but its effects on EEG activity and connectivity are unknown. The purpose of this study was to investigate EEG activity and connectivity changes after ESL treatment in persons with focal epilepsy (PFE).

Material and Methods: We performed a multicentre, longitudinal, retrospective, quantitative EEG study on a population of 22 PFE, and a group of 40 controls. We investigated the ESL-related changes of EEG power spectral activity and global connectivity [phase locking value (PLV), amplitude envelope correlation (AEC) and amplitude envelope correlation of orthogonalized signals (Ortho-AEC)] for standard frequency bands (delta to gamma). Seizure frequency was evaluated to assess ESL efficacy in our cohort.

Results: ESL significantly enhanced both global power spectral density and connectivity for all frequency bands, similarly for all connectivity measures. When compared to the control group, Post-ESL power was significantly higher in theta and gamma band. Pre-ESL connectivity values were significantly lower than control for all frequency bands. Post-ESL connectivity increased and the gap between the two groups was no longer significant. ESL induced a 52.7 ± 41.1% reduction of seizure frequency, with 55% of clinical responders (reduction of seizures ≥50%).

Discussion: ESL therapy induces significant enhancement of brain activity and connectivity. Post-ESL connectivity profile of epilepsy patients was similar to the one of healthy controls.

Theta Activity in the Left Dorsal Premotor Cortex During Action Re-Evaluation and Motor Reprogramming

The ability to rapidly adjust our actions to changes in the environment is a key function of human motor control. Previous work implicated the dorsal premotor cortex (dPMC) in the up-dating of action plans based on environmental cues. Here we used electroencephalography (EEG) to identify neural signatures of up-dating cue-action relationships in the dPMC and connected frontoparietal areas. Ten healthy subjects performed a pre-cued alternate choice task. Simple geometric shapes cued button presses with the right or left index finger. The shapes of the pre-cue and go-cue differed in two third of trials. In these incongruent trials, the go-cue prompted a re-evaluation of the pre-cued action plan, slowing response time relative to trials with identical cues. This re-evaluation selectively increased theta band activity without modifying activity in alpha and beta band. Source-based analysis revealed a widespread theta increase in dorsal and mesial frontoparietal areas, including dPMC, supplementary motor area (SMA), primary motor and posterior parietal cortices (PPC). Theta activity scaled positively with response slowing and increased more strongly when the pre-cue was invalid and required subjects to select the alternate response. Together, the results indicate that theta activity in dPMC and connected frontoparietal areas is involved in the re-adjustment of cue-induced action tendencies.

Reproducibility of EEG-MEG fusion source analysis of interictal spikes: Relevance in presurgical evaluation of epilepsy

Fusion of electroencephalography (EEG) and magnetoencephalography (MEG) data using maximum entropy on the mean method (MEM-fusion) takes advantage of the complementarities between EEG and MEG to improve localization accuracy. Simulation studies demonstrated MEM-fusion to be robust especially in noisy conditions such as single spike source localizations (SSSL). Our objective was to assess the reliability of SSSL using MEM-fusion on clinical data. We proposed to cluster SSSL results to find the most reliable and consistent source map from the reconstructed sources, the so-called consensus map. Thirty-four types of interictal epileptic discharges (IEDs) were analyzed from 26 patients with well-defined epileptogenic focus. SSSLs were performed on EEG, MEG, and fusion data and consensus maps were estimated using hierarchical clustering. Qualitative (spike-to-spike reproducibility rate, SSR) and quantitative (localization error and spatial dispersion) assessments were performed using the epileptogenic focus as clinical reference. Fusion SSSL provided significantly better results than EEG or MEG alone. Fusion found at least one cluster concordant with the clinical reference in all cases. This concordant cluster was always the one involving the highest number of spikes. Fusion yielded highest reproducibility (SSR EEG = 55%, MEG = 71%, fusion = 90%) and lowest localization error. Also, using only few channels from either modality (21EEG + 272MEG or 54EEG + 25MEG) was sufficient to reach accurate fusion. MEM-fusion with consensus map approach provides an objective way of finding the most reliable and concordant generators of IEDs. We, therefore, suggest the pertinence of SSSL using MEM-fusion as a valuable clinical tool for presurgical evaluation of epilepsy.