An evaluation of kurtosis beamforming in magnetoencephalography to localize the epileptogenic zone in drug resistant epilepsy patients

Towards the automated detection of interictal epileptiform discharges with magnetoencephalography

BACKGROUND

The analysis of clinical magnetoencephalography (MEG) in patients with epilepsy traditionally relies on visual identification of interictal epileptiform discharges (IEDs), which is time consuming and dependent on subjective criteria.

NEW METHOD

Here, we explore the ability of Independent Components Analysis (ICA) and Hidden Markov Modeling (HMM) to automatically detect and localize IEDs. We tested our pipelines on resting-state MEG recordings from 10 school-aged children with (multi)focal epilepsy.

RESULTS

In focal epilepsy patients, both pipelines successfully detected visually identified IEDs, but also revealed unidentified low-amplitude IEDs. Success was more mitigated in patients with multifocal epilepsy, as our automated pipeline missed IED activity associated with some foci-an issue that could be alleviated by post-hoc manual selection of epileptiform ICs or HMM states.

COMPARISON WITH EXISTING METHODS

We compared our results with visual IED detection by an experienced clinical magnetoencephalographer, getting heightened sensitivity and requiring minimal input from clinical practitioners.

CONCLUSIONS

IED detection based on ICA or HMM represents an efficient way to identify IED localization and timing. The development of these automatic IED detection algorithms provide a step forward in clinical MEG practice by decreasing the duration of MEG analysis and enhancing its sensitivity.

From cryogenic to on-scalp magnetoencephalography for the evaluation of paediatric epilepsy

Magnetoencephalography (MEG) is a neurophysiological technique based on the detection of brain magnetic fields. Whole-head MEG systems typically house a few hundred sensors requiring cryogenic cooling in a rigid one-size-fits-all (commonly adult-sized) helmet to keep a thermal insulation space. This leads to an increased brain-to-sensor distance in children, because of their smaller head circumference, and decreased signal-to-noise ratio. MEG allows detection and localization of interictal and ictal epileptiform discharges, and pathological high frequency oscillations, as a part of the presurgical assessment of children with refractory focal epilepsy, where electroencephalography is not contributive. MEG can also map the eloquent cortex before surgical resection. MEG also provides insights into the physiopathology of both generalized and focal epilepsy. On-scalp recordings based on cryogenic-free sensors have demonstrated their use in the field of childhood focal epilepsy and should become a reference technique for diagnosing epilepsy in the paediatric population. WHAT THIS PAPER ADDS: Magnetoencephalography (MEG) contributes to the diagnosis and understanding of paediatric epilepsy. On-scalp MEG recordings demonstrate some advantages over cryogenic MEG.

Standardized hierarchical adaptive Lp regression for noise robust focal epilepsy source reconstructions

OBJECTIVE

To investigate the ability of standardization to reduce source localization errors and measurement noise uncertainties for hierarchical Bayesian algorithms with L1- and L2-norms as priors in electroencephalography and magnetoencephalography of focal epilepsy.

METHODS

Description of the standardization methodology relying on the Hierarchical Bayesian framework, referred to as the Standardized Hierarchical Adaptive Lp-norm Regularization (SHALpR). The performance was tested using real data from two focal epilepsy patients. Simulated data that resembled the available real data was constructed for further localization and noise robustness investigation.

RESULTS

The proposed algorithms were compared to their non-standardized counterparts, Standardized low-resolution brain electromagnetic tomography, Standardized Shrinking LORETA-FOCUSS, and Dynamic statistical parametric maps. Based on the simulations, the standardized Hierarchical adaptive algorithm using L2-norm was noise robust for 10 dB signal-to-noise ratio (SNR), whereas the L1-norm prior worked robustly also with 5 dB SNR. The accuracy of the standardized L1-normed methodology to localize focal activity was under 1 cm for both patients.

CONCLUSIONS

Numerical results of the proposed methodology display improved localization and noise robustness. The proposed methodology also outperformed the compared methods when dealing with real data.

SIGNIFICANCE

The proposed standardized methodology, especially when employing the L1-norm, could serve as a valuable assessment tool in surgical decision-making.

Beamforming Seizures from the Temporal Lobe Using Magnetoencephalography

BACKGROUND

Surgical treatment of drug-resistant temporal lobe epilepsy (TLE) depends on proper identification of the seizure onset zone (SOZ) and differentiation of mesial, temporolimbic seizure onsets from temporal neocortical seizure onsets. Noninvasive source imaging using electroencephalography (EEG) and magnetoencephalography (MEG) can provide accurate information on interictal spike localization; however, EEG and MEG have low sensitivity for epileptiform activity restricted to deep temporolimbic structures. Moreover, in mesial temporal lobe epilepsy (MTLE), interictal spikes frequently arise in neocortical foci distant from the SOZ, rendering interictal spike localization potentially misleading for presurgical planning.

METHODS

In this study, we used two different beamformer techniques applied to the MEG signal of ictal events acquired during EEG-MEG recordings in six patients with TLE (three neocortical, three MTLE) in whom the ictal source localization results could be compared to ground truth SOZ localizations determined from intracranial EEG and/or clinical, neuroimaging, and postsurgical outcome evidence.

RESULTS

Beamformer analysis proved to be highly accurate in all cases and was able to identify focal SOZs in mesial, temporolimbic structures. In three patients, interictal spikes were absent, too complex for dipole modeling, or localized to anterolateral temporal neocortex distant to a mesial temporal SOZ, and thus unhelpful in presurgical investigation.

CONCLUSIONS

MEG beamformer source reconstruction is suitable for analysis of ictal events in TLE and can complement or supersede the traditional analysis of interictal spikes. The method outlined is applicable to any type of epileptiform event, expanding the information value of MEG and broadening its utility for presurgical recording in epilepsy.

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.

Pier K, Mullinger KJ, Marsh ED, Brookes MJ, Gaetz W. Mapping Interictal activity in epilepsy using a hidden Markov model: A magnetoencephalography study

Epilepsy is a highly heterogeneous neurological disorder with variable etiology, manifestation, and response to treatment. It is imperative that new models of epileptiform brain activity account for this variability, to identify individual needs and allow clinicians to curate personalized care. Here, we use a hidden Markov model (HMM) to create a unique statistical model of interictal brain activity for 10 pediatric patients. We use magnetoencephalography (MEG) data acquired as part of standard clinical care for patients at the Children's Hospital of Philadelphia. These data are routinely analyzed using excess kurtosis mapping (EKM); however, as cases become more complex (extreme multifocal and/or polymorphic activity), they become harder to interpret with EKM. We assessed the performance of the HMM against EKM for three patient groups, with increasingly complicated presentation. The difference in localization of epileptogenic foci for the two methods was 7 ± 2 mm (mean ± SD over all 10 patients); and 94% ± 13% of EKM temporal markers were matched by an HMM state visit. The HMM localizes epileptogenic areas (in agreement with EKM) and provides additional information about the relationship between those areas. A key advantage over current methods is that the HMM is a data-driven model, so the output is tuned to each individual. Finally, the model output is intuitive, allowing a user (clinician) to review the result and manually select the HMM epileptiform state, offering multiple advantages over previous methods and allowing for broader implementation of MEG epileptiform analysis in surgical decision-making for patients with intractable epilepsy.

Validating EEG, MEG and Combined MEG and EEG Beamforming for an Estimation of the Epileptogenic Zone in Focal Cortical Dysplasia

MEG and EEG source analysis is frequently used for the presurgical evaluation of pharmacoresistant epilepsy patients. The source localization of the epileptogenic zone depends, among other aspects, on the selected inverse and forward approaches and their respective parameter choices. In this validation study, we compare the standard dipole scanning method with two beamformer approaches for the inverse problem, and we investigate the influence of the covariance estimation method and the strength of regularization on the localization performance for EEG, MEG, and combined EEG and MEG. For forward modelling, we investigate the difference between calibrated six-compartment and standard three-compartment head modelling. In a retrospective study, two patients with focal epilepsy due to focal cortical dysplasia type IIb and seizure freedom following lesionectomy or radiofrequency-guided thermocoagulation (RFTC) used the distance of the localization of interictal epileptic spikes to the resection cavity resp. RFTC lesion as reference for good localization. We found that beamformer localization can be sensitive to the choice of the regularization parameter, which has to be individually optimized. Estimation of the covariance matrix with averaged spike data yielded more robust results across the modalities. MEG was the dominant modality and provided a good localization in one case, while it was EEG for the other. When combining the modalities, the good results of the dominant modality were mostly not spoiled by the weaker modality. For appropriate regularization parameter choices, the beamformer localized better than the standard dipole scan. Compared to the importance of an appropriate regularization, the sensitivity of the localization to the head modelling was smaller, due to similar skull conductivity modelling and the fixed source space without orientation constraint.

Sensitivity of a 29-Channel MEG Source Montage

In this paper, we study the performance of a source montage corresponding to 29 brain regions reconstructed from whole-head magnetoencephalographic (MEG) recordings, with the aim of facilitating the review of MEG data containing epileptiform discharges. Test data were obtained by superposing simulated signals from 100-nAm dipolar sources to a resting state MEG recording from a healthy subject. Simulated sources were placed systematically to different cortical locations for defining the optimal regularization for the source montage reconstruction and for assessing the detectability of the source activity from the 29-channel MEG source montage. The signal-to-noise ratio (SNR), computed for each source from the sensor-level and source-montage signals, was used as the evaluation parameter. Without regularization, the SNR from the simulated sources was larger in the sensor-level signals than in the source montage reconstructions. Setting the regularization to 2% increased the source montage SNR to the same level as the sensor-level SNR, improving the detectability of the simulated events from the source montage reconstruction. Sources producing a SNR of at least 15 dB were visually detectable from the source-montage signals. Such sources are located closer than about 75 mm from the MEG sensors, in practice covering all areas in the grey matter. The 29-channel source montage creates more focal signals compared to the sensor space and can significantly shorten the detection time of epileptiform MEG discharges for focus localization.

Testing covariance models for MEG source reconstruction of hippocampal activity

Beamforming is one of the most commonly used source reconstruction methods for magneto- and electroencephalography (M/EEG). One underlying assumption, however, is that distant sources are uncorrelated and here we tested whether this is an appropriate model for the human hippocampal data. We revised the Empirical Bayesian Beamfomer (EBB) to accommodate specific a-priori correlated source models. We showed in simulation that we could use model evidence (as approximated by Free Energy) to distinguish between different correlated and uncorrelated source scenarios. Using group MEG data in which the participants performed a hippocampal-dependent task, we explored the possibility that the hippocampus or the cortex or both were correlated in their activity across hemispheres. We found that incorporating a correlated hippocampal source model significantly improved model evidence. Our findings help to explain why, up until now, the majority of MEG-reported hippocampal activity (typically making use of beamformers) has been estimated as unilateral.

MEG detection of high frequency oscillations and intracranial-EEG validation in pediatric epilepsy surgery

OBJECTIVE

To assess the feasibility of automatically detecting high frequency oscillations (HFOs) in magnetoencephalography (MEG) recordings in a group of ten paediatric epilepsy surgery patients who had undergone intracranial electroencephalography (iEEG).

METHODS

A beamforming source-analysis method was used to construct virtual sensors and an automatic algorithm was applied to detect HFOs (80-250 Hz). We evaluated the concordance of MEG findings with the sources of iEEG HFOs, the clinically defined seizure onset zone (SOZ), the location of resected brain structures, and with post-operative outcome.

RESULTS

In 8/9 patients there was good concordance between the sources of MEG HFOs and iEEG HFOs and the SOZ. Significantly more HFOs were detected in iEEG relative to MEG t(71) = 2.85, p < .05. There was good concordance between sources of MEG HFOs and the resected area in patients with good and poor outcome, however HFOs were also detected outside of the resected area in patients with poor outcome.

CONCLUSION

Our findings demonstrate the feasibility of automatically detecting HFOs non-invasively in MEG recordings in paediatric patients, and confirm compatibility of results with invasive recordings.

SIGNIFICANCE

This approach provides support for the non-invasive detection of HFOs to aid surgical planning and potentially reduce the need for invasive monitoring, which is pertinent to paediatric patients.

Delineation of epileptogenic zones with high frequency magnetic source imaging based on kurtosis and skewness

BACKGROUND

Neuromagnetic high frequency brain signals (HFBS, > 80 Hz) are a new biomarker for localization of epileptogenic zones (EZs) for pediatric epilepsy.

METHODS

Twenty three children with drug-resistant epilepsy and age/sex matched healthy controls were studied with magnetoencephalography (MEG). Epileptic HFBS in 80-250 Hz and 250-600 Hz were quantitatively determined by comparing with normative controls in terms of kurtosis and skewness. Magnetic sources of epileptic HFBS were localized and then compared to clinical EZs determined by invasive recordings and surgical outcomes.

RESULTS

Kurtosis and skewness of HFBS were significantly elevated in epilepsy patients compared to healthy controls (p < 0,001 and p < 0.0001, respectively). Sources of elevated MEG signals in comparison to normative data were co-localized to EZs for 22 (22/23, 96 %) patients.

CONCLUSIONS

The results indicate, for the first time, that epileptic HFBS can be noninvasively quantified by measuring kurtosis and skewness in MEG data. Magnetic source imaging based on kurtosis and skewness can accurately localize EZs.

SIGNIFICANCE

Source imaging of kurtosis and skewness of MEG HFBS provides a novel way for preoperative localization of EZs for epilepsy surgery.

Interictal spike localization for epilepsy surgery using magnetoencephalography beamforming

OBJECTIVE

Magnetoencephalography (MEG) kurtosis beamforming is an automated localization method for focal epilepsy. Visual examination of virtual sensors, which are source activities reconstructed by beamforming, can improve performance but can be time-consuming for neurophysiologists. We propose a framework to automate the method and evaluate its effectiveness against surgical resections and outcomes.

METHODS

We retrospectively analyzed MEG recordings of 13 epilepsy surgery patients who had one-year minimum post-operative follow-up. Kurtosis beamforming was applied and manual inspection was confined to morphological clusters. The region with the Maximum Interictal Spike Frequency (MISF) was validated against prospectively modelled sLORETA solutions and surgical resections linked to outcome.

RESULTS

Our approach localized spikes in 12 out of 13 patients. In eight patients with Engel I surgical outcomes, beamforming MISF regions were concordant with surgical resection at overlap level for five patients and at lobar level for three patients. The MISF regions localized to spike onset and propagation modelled by sLORETA in two and six patients, respectively.

CONCLUSIONS

Automated beamforming using MEG can predict postoperative seizure freedom at the lobar level but tends to localize propagated MEG spikes.

SIGNIFICANCE

MEG beamforming may contribute to non-invasive procedures to predict surgical outcome for patients with drug-refractory focal epilepsy.

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.

Comparison of beamformer implementations for MEG source localization

Beamformers are applied for estimating spatiotemporal characteristics of neuronal sources underlying measured MEG/EEG signals. Several MEG analysis toolboxes include an implementation of a linearly constrained minimum-variance (LCMV) beamformer. However, differences in implementations and in their results complicate the selection and application of beamformers and may hinder their wider adoption in research and clinical use. Additionally, combinations of different MEG sensor types (such as magnetometers and planar gradiometers) and application of preprocessing methods for interference suppression, such as signal space separation (SSS), can affect the results in different ways for different implementations. So far, a systematic evaluation of the different implementations has not been performed. Here, we compared the localization performance of the LCMV beamformer pipelines in four widely used open-source toolboxes (MNE-Python, FieldTrip, DAiSS (SPM12), and Brainstorm) using datasets both with and without SSS interference suppression.

We analyzed MEG data that were i) simulated, ii) recorded from a static and moving phantom, and iii) recorded from a healthy volunteer receiving auditory, visual, and somatosensory stimulation. We also investigated the effects of SSS and the combination of the magnetometer and gradiometer signals. We quantified how localization error and point-spread volume vary with the signal-to-noise ratio (SNR) in all four toolboxes.

When applied carefully to MEG data with a typical SNR (3–15 ​dB), all four toolboxes localized the sources reliably; however, they differed in their sensitivity to preprocessing parameters. As expected, localizations were highly unreliable at very low SNR, but we found high localization error also at very high SNRs for the first three toolboxes while Brainstorm showed greater robustness but with lower spatial resolution. We also found that the SNR improvement offered by SSS led to more accurate localization.

•

Different beamformer implementations are reported to sometimes yield differing source estimates for the same MEG data.

•

We compared beamformers in four major open-source MEG analysis toolboxes.

•

All toolboxes provide consistent and accurate results with 3–15-dB input SNR.

•

However, localization errors are high at very high input SNR for the tested scalar beamformers.

•

We discuss the critical differences between the implementations.

Kurtosis and skewness of high-frequency brain signals are altered in paediatric epilepsy

Intracranial studies provide solid evidence that high-frequency brain signals are a new biomarker for epilepsy. Unfortunately, epileptic (pathological) high-frequency signals can be intermingled with physiological high-frequency signals making these signals difficult to differentiate. Recent success in non-invasive detection of high-frequency brain signals opens a new avenue for distinguishing pathological from physiological high-frequency signals. The objective of the present study is to characterize pathological and physiological high-frequency signals at source levels by using kurtosis and skewness analyses. Twenty-three children with medically intractable epilepsy and age-/gender-matched healthy controls were studied using magnetoencephalography. Magnetoencephalographic data in three frequency bands, which included 2–80 Hz (the conventional low-frequency signals), 80–250 Hz (ripples) and 250–600 Hz (fast ripples), were analysed. The kurtosis and skewness of virtual electrode signals in eight brain regions, which included left/right frontal, temporal, parietal and occipital cortices, were calculated and analysed. Differences between epilepsy and controls were quantitatively compared for each cerebral lobe in each frequency band in terms of kurtosis and skewness measurements. Virtual electrode signals from clinical epileptogenic zones and brain areas outside of the epileptogenic zones were also compared with kurtosis and skewness analyses. Compared to controls, patients with epilepsy showed significant elevation in kurtosis and skewness of virtual electrode signals. The spatial and frequency patterns of the kurtosis and skewness of virtual electrode signals among the eight cerebral lobes in three frequency bands were also significantly different from that of the controls (2–80 Hz, P < 0.001; 80–250 Hz, P < 0.00001; 250–600 Hz, P < 0.0001). Compared to signals from non-epileptogenic zones, virtual electrode signals from epileptogenic zones showed significantly altered kurtosis and skewness (P < 0.001). Compared to normative data from the control group, aberrant virtual electrode signals were, for each patient, more pronounced in the epileptogenic lobes than in other lobes(kurtosis analysis of virtual electrode signals in 250–600 Hz; odds ratio = 27.9; P < 0.0001). The kurtosis values of virtual electrode signals in 80–250 and 250–600 Hz showed the highest sensitivity (88.23%) and specificity (89.09%) for revealing epileptogenic lobe, respectively. The combination of virtual electrode and kurtosis/skewness measurements provides a new quantitative approach to distinguishing pathological from physiological high-frequency signals for paediatric epilepsy. Non-invasive identification of pathological high-frequency signals may provide novel important information to guide clinical invasive recordings and direct surgical treatment of epilepsy.

Interictal magnetoencephalography in parietal lobe epilepsy - Comparison of equivalent current dipole and beamformer (SAMepi) analysis

•

MEG may aid in the localization of the epileptogenic zone in the parietal lobe.

•

SAMepi – a novel kurtosis beamformer – results in localizations similar to those of the ECD analysis.

•

A unifocal result in both the ECD and the SAMepi analysis is associated with a good clinical outcome.

Objective

To evaluate a novel analysis method (SAMepi) in the localization of interictal epileptiform magnetoencephalographic (MEG) activity in parietal lobe epilepsy (PLE) patients in comparison with equivalent current dipole (ECD) analysis.

Methods

We analyzed the preoperative interictal MEG of 17 operated PLE patients utilizing visual analysis and: (1) ECD with a spherical conductor model; (2) ECD with a boundary element method (BEM) conductor model; and (3) SAMepi – a kurtosis beamformer method. Localization results were compared between the three methods, to the location of the resection and to the clinical outcome.

Results

Fourteen patients had an epileptiform finding in the visual analysis; SAMepi detected spikes in 11 of them. A unifocal finding in both the ECD and in the SAMepi analysis was associated with a better chance of seizure-freedom (p = 0.02). There was no significant difference in the distances from the unifocal MEG localizations to the nearest border of the resection between the different analysis methods.

Conclusions

Localizations of unifocal interictal spikes detected by SAMepi did not significantly differ from the conventional ECD localizations.

Significance

SAMepi – a novel semiautomatic analysis method – is useful in localizing interictal epileptiform MEG activity in the presurgical evaluation of parietal lobe epilepsy patients.

Hyperactive frontolimbic and frontocentral resting-state gamma connectivity in major depressive disorder

BACKGROUND

Major depressive disorder (MDD) is a system-level disorder affecting multiple functionally integrated cerebral networks. Nevertheless, their temporospatial organization and potential disturbance remain mostly unknown. The present report tested the hypothesis that deficient temporospatial network organization separates MDD and healthy controls (HC), and is linked to symptom severity of the disorder.

METHODS

Eyes-closed resting-state magnetoencephalographic (MEG) recordings were obtained from twenty-two MDD and twenty-two HC subjects. Beamforming source localization and functional connectivity analysis were applied to identify frequency-specific network interactions. Then, a novel virtual cortical resection approach was used to pinpoint putatively critical network controllers, accounting for aberrant cerebral connectivity patterns in MDD.

RESULTS

We found significantly elevated frontolimbic and frontocentral connectivity mediated by gamma (30-48 Hz) activity in MDD versus HC, and the right amygdala was the key differential network controller accounting for aberrant cerebral connectivity patterns in MDD. Furthermore, this frontolimbic and frontocentral gamma-band hyper-connectivity was positively correlated with depression severity.

LIMITATIONS

The overall sample size was small, and we found significant effects in the deep limbic regions with resting-state MEG, the reliability of which was difficult to corroborate further.

CONCLUSIONS

Overall, these findings support a notion that the right amygdala critically controls the exaggerated gamma-band frontolimbic and frontocentral connectivity in MDD during the resting-state condition, which potentially constitutes pre-established aberrant pathways during task processing and contributes to MDD pathology.

Synthetic aperture magnetometry and excess kurtosis mapping of Magnetoencephalography (MEG) is predictive of epilepsy surgical outcome in a large pediatric cohort

OBJECTIVE

Resective surgery is the most effective treatment option for patients with refractory epilepsy; however identification of patients who will benefit from epilepsy surgery remains challenging. Synthetic aperture magnetometry and excess kurtosis mapping (SAM(g2)) of magnetoencephalography (MEG) is a non-invasive tool that warrants further examination in the pediatric epilepsy population. Here, we examined the utility of MEG with SAM(g2) to determine if MEG epileptiform foci correlates with surgical outcome and to develop a predictive model incorporating MEG information to best assess likelihood of seizure improvement/freedom from resective surgery.

METHODS

564 subjects who had MEG at the Children's Hospital of Philadelphia between 2010-2015 were screened. Clinical epilepsy history and prior electrographic records were extracted and reviewed and correlated with MEG findings. MEG assessments were made by both a neurologist and neuroradiologist. Predictive models were developed to assess the utility of MEG in determining Engel class at one year and five years after resective epilepsy surgery.

RESULTS

The number of MEG spike foci was highly associated with Engel class outcome at both one year and five years; however, using MEG data in isolation was not significantly predictive of 5 year surgical outcome. When combined with clinical factors; scalp EEG (single ictal onset zone), MRI (lesional or not), age and sex in a logistic regression model MEG foci was significant for Engel class outcome at both 1 year (p = 0.03) and 5 years (0.02). The percent correctly classified for Engel class at one year was 78.43% and the positive predictive value was 71.43.

SIGNIFICANCE

MEG using SAM(g2) analysis in an important non-invasive tool in the identification of those patients who will benefit most from surgery. Integrating MEG data analysis into pre-surgical evaluation can help to predict epilepsy outcome after resective surgery in the pediatric population if utilized with skilled interpretation.

Automated Source Estimation of Scalp EEG Epileptic Activity Using eLORETA Kurtosis Analysis

OBJECTIVES

eLORETA (exact low-resolution brain electromagnetic tomography) is a technique created by Pascual-Marqui et al. [Int J Psychophysiol. 1994 Oct; 18(1): 49-65] for the 3-dimensional representation of current source density in the brain by electroencephalography (EEG) data. Kurtosis analysis allows for the identification of spiky activity in the brain. In this study, we focused on the evaluation of the reliability of eLORETA kurtosis analysis. For this purpose, the results of eLORETA kurtosis source localization of paroxysmal activity in EEG were compared with those of eLORETA current source density (CSD) analysis of EEG data in 3 epilepsy patients with partial seizures.

METHODS

EEG was measured using a digital EEG system with 19 channels. We set the bandpass filter at traditional frequency band settings (1-4, 4-8, 8-15, 15-30, and 30-60 Hz) and 5-10 and 20-70 Hz and performed eLORETA kurtosis to compare the source localization of paroxysmal activity with that of visual interpretation of EEG data and CSD analysis of eLORETA in focal epilepsy patients.

RESULTS

The eLORETA kurtosis analysis of EEG data preprocessed by bandpass filtering from 20 to 70 Hz and traditional frequency band settings did not show any discrete paroxysmal source activity compatible with the results of CSD analysis of eLORETA. In all 3 cases, eLORETA kurtosis analysis filtered at 5-10 Hz showed paroxysmal activities in the theta band, which were all consistent with the visual inspection results and the CSD analysis results.

DISCUSSION

Our findings suggested that eLORETA kurtosis analysis of EEG data might be useful for the identification of spiky paroxysmal activity sources in epilepsy patients. Since EEG is widely used in the clinical practice of epilepsy, eLORETA kurtosis analysis is a promising method that can be applied to epileptic activity mapping.

An enriched granger causal model allowing variable static anatomical constraints

The anatomical connectivity constrains but does not fully determine functional connectivity, especially when one explores into the dynamics over the course of a trial. Therefore, an enriched granger causal model (GCM) integrated with anatomical prior information is proposed in this study, to describe the dynamic effective connectivity to distinguish the depression and explore the pathogenesis of depression. In the proposed frame, the anatomical information was converted via an optimized transformation model, which was then integrated into the normal GCM by variational bayesian model. Magnetoencephalography (MEG) signals and diffusion tensor imaging (DTI) of 24 depressive patients and 24 matched controls were utilized for performance comparison. Together with the sliding windowed MEG signals under sad facial stimuli, the enriched GCM was applied to calculate the regional-pair dynamic effective connectivity, which were repeatedly sifted via feature selection and fed into different classifiers. From the aspects of model errors and recognition accuracy rates, results supported the superiority of the enriched GCM with anatomical priors over the normal GCM. For the effective connectivity with anatomical priors, the best subject discrimination accuracy of SVM was 85.42% (the sensitivity was 87.50% and the specificity was 83.33%). Furthermore, discriminative feature analysis suggested that the enriched GCM that detect the variable anatomical constraint on function could better detect more stringent and less dynamic brain function in depression. The proposed approach is valuable in dynamic functional dysfunction exploration in depression and could be useful for depression recognition.

Lateralization Value of Low Frequency Band Beamformer Magnetoencephalography Source Imaging in Temporal Lobe Epilepsy

Objective: In presurgical evaluation of temporal lobe epilepsy (TLE), selection of the resection side is challenging when bilateral temporal epileptiform discharges or structural abnormalities are present. We aim to evaluate the lateralization value of beamformer analysis of magnetoencephalography (MEG) in TLE. Methods: MEG data from 14 TLE patients were analyzed through beamformer analysis. We measured the hemispherical power distribution of beamformer sources and calculated the lateralization index (LI). We calculated the LI at multiple frequencies to explore the frequency dependency and at the delta frequency to define laterality. LI values ranging from -1 to -0.05 indicated right hemispheric dominance. LI values ranging from 0.05 to 1 indicated left hemispheric dominance. LI values ranging from -0.05 to 0.05 defined bilaterality. We measured the power of beamformer sources with a 9-s duration to explore time dependency. Results: The beamformer analysis showed that 10/14 patients had power dominance ipsilateral to resection. The delta frequency band had a higher lateralization value than other frequency bands. A time-dependent power fluctuation was found in the delta frequency band. Conclusions: MEG beamformer analysis, especially in the delta band, might efficiently provide additional information regarding lateralization in TLE.