Dynamic statistical parametric mapping for analyzing ictal magnetoencephalographic spikes in patients with intractable frontal lobe epilepsy

Automatic Detection and Classification of Epileptic Seizures from EEG Data: Finding Optimal Acquisition Settings and Testing Interpretable Machine Learning Approach

Deep learning (DL) is emerging as a successful technique for automatic detection and differentiation of spontaneous seizures that may otherwise be missed or misclassified. Herein, we propose a system architecture based on top-performing DL models for binary and multigroup classifications with the non-overlapping window technique, which we tested on the TUSZ dataset. The system accurately detects seizure episodes (87.7% Sn, 91.16% Sp) and carefully distinguishes eight seizure types (95-100% Acc). An increase in EEG sampling rate from 50 to 250 Hz boosted model performance: the precision of seizure detection rose by 5%, and seizure differentiation by 7%. A low sampling rate is a reasonable solution for training reliable models with EEG data. Decreasing the number of EEG electrodes from 21 to 8 did not affect seizure detection but worsened seizure differentiation significantly: 98.24 ± 0.17 vs. 85.14 ± 3.14% recall. In detecting epileptic episodes, all electrodes provided equally informative input, but in seizure differentiation, their informative value varied. We improved model explainability with interpretable ML. Activation maximization highlighted the presence of EEG patterns specific to eight seizure types. Cortical projection of epileptic sources depicted differences between generalized and focal seizures. Interpretable ML techniques confirmed that our system recognizes biologically meaningful features as indicators of epileptic activity in EEG.

Clinical significance of ictal magnetoencephalography in patients undergoing epilepsy surgery

OBJECTIVE

The significance of ictal magnetoencephalography (MEG) is not well appreciated. We evaluated the relationships between ictal MEG, MRI, intracranial electroencephalography (ICEEG), surgery and postoperative seizure outcome.

METHODS

A total of 45 patients (46 cases) with ictal MEG who underwent epilepsy surgery was included. We examined the localization of each modality, surgical resection area and seizure freedom after surgery.

RESULTS

Twenty-one (45.7%) out of 46 cases were seizure-free at more than 6 months follow-up. Median duration of postoperative follow-up was 16.5 months. The patients in whom ictal, interictal single equivalent current dipole (SECD) and MRI lesion localization were completely included in the resection had a higher chance of being seizure-free significantly (p < 0.05). Concordance between ictal and interictal SECD localizations was significantly associated with seizure-freedom. Concordance between MRI lesion and ictal SECD, concordance between ictal ICEEG and ictal and interictal SECD, as well as concordance between ictal ICEEG and MRI lesion were significantly associated with seizure freedom.

CONCLUSIONS

Ictal MEG can contribute useful information for delineating the resection area in epilepsy surgery.

SIGNIFICANCE

Resection should include ictal, interictal SECDs and MRI lesion localization, when feasible. Concordant ictal and interictal SECDs on MEG can be a favorable predictor of seizure freedom.

A Graph Fourier Transform Based Bidirectional Long Short-Term Memory Neural Network for Electrophysiological Source Imaging

Electrophysiological source imaging (ESI) refers to the process of reconstructing underlying activated sources on the cortex given the brain signal measured by Electroencephalography (EEG) or Magnetoencephalography (MEG). Due to the ill-posed nature of ESI, solving ESI requires the design of neurophysiologically plausible regularization or priors to guarantee a unique solution. Recovering focally extended sources is more challenging, and traditionally uses a total variation regularization to promote spatial continuity of the activated sources. In this paper, we propose to use graph Fourier transform (GFT) based bidirectional long-short term memory (BiLSTM) neural network to solve the ESI problem. The GFT delineates the 3D source space into spatially high, medium and low frequency subspaces spanned by corresponding eigenvectors. The low frequency components can naturally serve as a spatially low-band pass filter to reconstruct extended areas of source activation. The BiLSTM is adopted to learn the mapping relationship between the projection of low-frequency graph space and the recorded EEG. Numerical results show the proposed GFT-BiLSTM outperforms other benchmark algorithms in synthetic data under varied signal-to-noise ratios (SNRs). Real data experiments also demonstrate its capability of localizing the epileptogenic zone of epilepsy patients with good accuracy.

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.

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.

Magnetoencephalographic Mapping of Epileptic Spike Population Using Distributed Source Analysis: Comparison With Intracranial Electroencephalographic Spikes

INTRODUCTION

This study evaluates magnetoencephalographic (MEG) spike population as compared with intracranial electroencephalographic (IEEG) spikes using a quantitative method based on distributed source analysis.

METHODS

We retrospectively studied eight patients with medically intractable epilepsy who had an MEG and subsequent IEEG monitoring. Fifty MEG spikes were analyzed in each patient using minimum norm estimate. For individual spikes, each vertex in the source space was considered activated when its source amplitude at the peak latency was higher than a threshold, which was set at 50% of the maximum amplitude over all vertices. We mapped the total count of activation at each vertex. We also analyzed 50 IEEG spikes in the same manner over the intracranial electrodes and created the activation count map. The location of the electrodes was obtained in the MEG source space by coregistering postimplantation computed tomography to MRI. We estimated the MEG- and IEEG-active regions associated with the spike populations using the vertices/electrodes with a count over 25.

RESULTS

The activation count maps of MEG spikes demonstrated the localization associated with the spike population by variable count values at each vertex. The MEG-active region overlapped with 65 to 85% of the IEEG-active region in our patient group.

CONCLUSIONS

Mapping the MEG spike population is valid for demonstrating the trend of spikes clustering in patients with epilepsy. In addition, comparison of MEG and IEEG spikes quantitatively may be informative for understanding their relationship.

Combined Infection of Human Herpes Viruses: A Risk Factor for Subsequent Colectomy in Ulcerative Colitis

BACKGROUND

Little is known about the prevalence and pathogenicity of human herpes viruses except for cytomegalovirus (CMV) in patients with inflammatory bowel disease (IBD). The aim of this study was to determine the prevalence of human herpes viruses on colonic mucosa in patients with IBD and assess the long-term outcomes in these patients.

METHODS

We examined the colonic mucosal specimens from 66 patients with ulcerative colitis (UC), 54 patients with Crohn's disease (CD), and 29 healthy patients to identify the 6 most common types of human herpes virus, using multiplex polymerase chain reaction (PCR) technique.

RESULTS

Herpes simplex virus (HSV)-1/2 and varicella-zoster virus (VZV) were not detected in any of the groups. There was a higher prevalence of Epstein-Barr virus (EBV) (21.2%) and CMV (15.1%) in patients with UC than in patients with CD (EBV 9.3%, CMV 0%) and patents in the healthy control group (EBV 0%, CMV 3.4%). The prevalence of human herpes virus (HHV)-6A/B and HHV-7 was not statistically different among the groups. Five UC patients with inflammation had coexisting CMV and EBV or HHV-6. The combined infection of CMV with EBV or HHV-6 was a significant and independent prognostic factor for subsequent colectomy in patients with UC.

CONCLUSIONS

The increased prevalence of CMV coexisting with EBV/HHV-6 infection was associated with the clinical course in patients with UC. 10.1093/ibd/izy005_video1izy005_Video_15786489376001.

Multimodal investigation of epileptic networks: The case of insular cortex epilepsy

The insula is a deep cortical structure sharing extensive synaptic connections with a variety of brain regions, including several frontal, temporal, and parietal structures. The identification of the insular connectivity network is obviously valuable for understanding a number of cognitive processes, but also for understanding epilepsy since insular seizures involve a number of remote brain regions. Ultimately, knowledge of the structure and causal relationships within the epileptic networks associated with insular cortex epilepsy can offer deeper insights into this relatively neglected type of epilepsy enabling the refining of the clinical approach in managing patients affected by it. In the present chapter, we first review the multimodal noninvasive tests performed during the presurgical evaluation of epileptic patients with drug refractory focal epilepsy, with particular emphasis on their value for the detection of insular cortex epilepsy. Second, we review the emerging multimodal investigation techniques in the field of epilepsy, that aim to (1) enhance the detection of insular cortex epilepsy and (2) unveil the architecture and causal relationships within epileptic networks. We summarize the results of these approaches with emphasis on the specific case of insular cortex epilepsy.

Analysis of ictal magnetoencephalography using gradient magnetic-field topography (GMFT) in patients with neocortical epilepsy

OBJECTIVE

We aimed to validate the usefulness of gradient magnetic-field topography (GMFT) for analysis of ictal magnetoencephalography (MEG) in patients with neocortical epilepsy.

METHODS

We identified 13 patients presenting with an ictal event during preoperative MEG. We applied equivalent current dipole (ECD) estimation and GMFT to detect and localize the ictal MEG onset, and compared these methods with the ictal onset zone (IOZ) derived from chronic intracranial electroencephalography. The surgical resection areas and outcomes were also evaluated.

RESULTS

GMFT detected and localized the ictal MEG onset in all patients, whereas ECD estimation showed localized ECDs in only 2. The delineation of GMFT was concordant with the IOZ at the gyral-unit level in 10 of 12 patients (83.3%). The detectability and precision of delineation of ictal MEG activity by GMFT were significantly superior to those of ECD (p<0.05 and p<0.01, respectively). Complete resection of the IOZ in the concordant group provided seizure freedom in 3 patients, whereas seizures remained in 9 patients who had incomplete resections.

CONCLUSIONS

Because of its higher spatial resolution, GMFT of ictal MEG is superior to conventional ECD estimation in patients with neocortical epilepsy.

SIGNIFICANCE

Ictal MEG study is a useful tool to estimate the seizure onset in patients with neocortical epilepsy.

Localized N20 Component of Somatosensory Evoked Magnetic Fields in Frontoparietal Brain Tumor Patients Using Noise-Normalized Approaches

PURPOSE

To localize sensorimotor cortical activation in 10 patients with frontoparietal tumors using quantitative magnetoencephalography (MEG) with noise-normalized approaches.

MATERIAL AND METHODS

Somatosensory evoked magnetic fields (SEFs) were elicited in 10 patients with somatosensory tumors and in 10 control participants using electrical stimulation of the median nerve via the right and left wrists. We localized the N20m component of the SEFs using dynamic statistical parametric mapping (dSPM) and standardized low-resolution brain electromagnetic tomography (sLORETA) combined with 3D magnetic resonance imaging (MRI). The obtained coordinates were compared between groups. Finally, we statistically evaluated the N20m parameters across hemispheres using non-parametric statistical tests.

RESULTS

The N20m sources were accurately localized to Brodmann area 3b in all members of the control group and in seven of the patients; however, the sources were shifted in three patients relative to locations outside the primary somatosensory cortex (SI). Compared with the affected (tumor) hemispheres in the patient group, N20m amplitudes and the strengths of the current sources were significantly lower in the unaffected hemispheres and in both hemispheres of the control group. These results were consistent for both dSPM and sLORETA approaches.

CONCLUSION

Tumors in the sensorimotor cortex lead to cortical functional reorganization and an increase in N20m amplitude and current-source strengths. Noise-normalized approaches for MEG analysis that are integrated with MRI show accurate and reliable localization of sensorimotor function.

Source localization of the seizure onset zone from ictal EEG/MEG data

INTRODUCTION

Surgical treatment of drug-resistant epilepsy relies on the identification of the seizure onset zone (SOZ) and often requires intracranial EEG (iEEG). We have developed a new approach for non-invasive magnetic and electric source imaging of the SOZ (MSI-SOZ and ESI-SOZ) from ictal magnetoencephalography (MEG) and EEG recordings, using wavelet-based Maximum Entropy on the Mean (wMEM) method. We compared the performance of MSI-SOZ and ESI-SOZ with interictal spike source localization (MSI-spikes and ESI-spikes) and clinical localization of the SOZ (i.e., based on iEEG or lesion topography, denoted as clinical-SOZ).

METHODS

A total of 46 MEG or EEG seizures from 13 patients were analyzed. wMEM was applied around seizure onset, centered on the frequency band showing the strongest power change. Principal component analysis applied to spatiotemporal reconstructed wMEM sources (0.4-1 s around seizure onset) identified the main spatial pattern of ictal oscillations. Qualitative sublobar concordance and quantitative measures of distance and spatial overlaps were estimated to compare MSI/ESI-SOZ with MSI/ESI-Spikes and clinical-SOZ.

RESULTS

MSI/ESI-SOZ were concordant with clinical-SOZ in 81% of seizures (MSI 90%, ESI 64%). MSI-SOZ was more accurate and identified sources closer to the clinical-SOZ (P = 0.012) and to MSI-Spikes (P = 0.040) as compared with ESI-SOZ. MSI/ESI-SOZ and MSI/ESI-Spikes did not differ in terms of concordance and distance from the clinical-SOZ.

CONCLUSIONS

wMEM allows non-invasive localization of the SOZ from ictal MEG and EEG. MSI-SOZ performs better than ESI-SOZ. MSI/ESI-SOZ can provide important additional information to MSI/ESI-Spikes during presurgical evaluation. Hum Brain Mapp 37:2528-2546, 2016. © 2016 Wiley Periodicals, Inc.

Statistical parametric mapping for analyzing interictal magnetoencephalography in patients with left frontal lobe epilepsy

PURPOSE

Frontal lobe epilepsy is a common epileptic disorder and is characterized by recurring seizures that arise in the frontal lobes. The purpose of this study is to identify the epileptogenic regions and other abnormal regions in patients with left frontal lobe epilepsy (LFLE) based on the magnetoencephalogram (MEG), and to understand the effects of clinical variables on brain activities in patients with LFLE.

METHOD

Fifteen patients with LFLE (23.20 ± 8.68 years, 6 female and 9 male) and 16 healthy controls (23.13 ± 7.66 years, 6 female and 10 male) were included in resting-stage MEG examinations. Epileptogenic regions of LFLE patients were confirmed by surgery. Regional brain activations were quantified using statistical parametric mapping (SPM). The correlation between the activations of the abnormal brain regions and the clinical seizure parameters were computed for LFLE patients.

RESULTS

Brain activations of LFLE patients were significantly elevated in left superior/middle/inferior frontal gyri, postcentral gyrus, inferior temporal gyrus, insula, parahippocampal gyrus and amygdala, including the epileptogenic regions. Remarkable decreased activations were found mainly in the left parietal gyrus and precuneus. There is a positive correlation between the duration of the epilepsy (in month) and activations of the abnormal regions, while no relation was found between age of seizure onset (year), seizure frequency and the regions of the abnormal activity of the epileptic patients.

CONCLUSION

Our findings suggest that the aberrant brain activities of LFLE patients were not restricted to the epileptogenic zones. Long duration of epilepsy might induce further functional damage in patients with LFLE.

Ictal Magnetic Source Imaging in Presurgical Assessment

Ictal MEG recordings constitute rare data. The objective of this study was to evaluate ictal magnetic source localization (MSI), using two algorithms: linearly constrained minimum variance (LCMV), a beamforming technique and equivalent current dipole (ECD). Ictal MSI was studied in six patients. Three of them were undergoing post-operative re-evaluation. For all patients, results were validated by the stereoelectroencephalographic (SEEG) definition of the epileptogenic zone (EZ). EZ was quantified using the epileptogenicity index (EI) method, which accounts for both the propensity of a brain area to generate rapid discharges and the time for this area to become involved in the seizure. EI values range from 0 (no epileptogenicity) to 1 (maximal epileptogenicity). Levels of concordance between ictal MSI and EZ were determined as follows: A: ictal MSI localized the site whose value EI = 1, B: MSI localized a part of the EZ (not corresponding to the maximal value of EI = 1), C: a region could be identified on ictal MSI but not on SEEG, D: a region could be identified on SEEG but not on MSI, E: different regions were localized on MSI and SEEG. Ictal MEG pattern consisted of rhythmic activities between 10 and 20 Hz for all patients. For LCMV (first maxima), levels of concordance were A (two cases), B (two cases) and E (two cases). For ECD fitted on each time point separately (location characterized by the best goodness-of-fit value), levels of concordance were A (one case), B (one case), D (three cases) and E (one case). For ECD calculated for the whole time window, levels of concordance were A (two cases) and D (four cases). Source localization methods performed on rhythmic patterns can localize the EZ as validated by SEEG. In terms of concordance, LCMV was superior to ECD. In some cases, LCMV allows extraction of several maxima that could reflect ictal dynamics. In a medial temporal lobe epilepsy case, ictal MSI indicated an area of delayed propagation and was non-contributory to the presurgical assessment.

Usefulness of multiple frequency band source localizations in ictal MEG

OBJECTIVE

We evaluated the diagnostic value of multiple frequency band MEG source localization within a wide time window during the preictal period.

METHODS

Data for 13 epilepsy patients who showed an ictal event during MEG were analyzed. Several seconds of preictal data were localized in the theta, alpha, beta, and gamma bands by using wavelet transformation and the sLORETA algorithm. The same analysis was performed with narrow time and frequency band. Localization concordances to the surgically resected area were compared.

RESULTS

Source localization in the gamma band for a 10s window before ictal onset showed best concordance to the resection cavity. Eight of 13 patients showed sub-lobar concordance in the 10s gamma band localization, whereas 3 showed concordance in the narrow time and frequency analysis. Four of 7 patients with focal cortical dysplasia (FCD) achieved seizure-free outcome, and all 4 showed sub-lobar concordance.

CONCLUSIONS

A 10s time window gamma source localization method can be used to delineate the epileptogenic zone.

SIGNIFICANCE

The use of a long period during preictal gamma source localization has the potential to become a localizing biomarker of the epileptogenic zone in candidates for surgical intervention, especially in MRI-suspected FCD.

Sparse MEG source imaging for reconstructing dynamic sources of interictal spikes in partial epilepsy

PURPOSE

The present study aimed to test the feasibility of a novel neuroimaging technique, that is, variation-based sparse cortical current density (VB-SCCD) imaging algorithm, in noninvasively estimating location and extent of epileptic sources from interictal magnetoencephalography (MEG) data.

METHODS

A total of 108 interictal spikes from 3 partial epilepsy patients were selected to perform VB-SCCD source analysis. Cortical sources were identified at spike peaks, rising phases, and entire spikes, respectively, from all interictal spikes in each patient, to estimate source locations and extents, and validated using presurgical evaluation data. Other source analysis methods, that is, minimum norm estimate and sparse source imaging were also performed for comparison.

RESULTS

Cortical sources reconstructed by VB-SCCD that are consistent with clinical presurgical evaluation outcomes have detection rates of 65.8% at spike peaks, 85.1% during rising phases, and 92.6% in entire spikes. Stable spatiotemporal patterns of reconstructed cortical sources were also obtained using VB-SCCD, which provide more insights about the formation and propagation of interictal epileptic activity.

CONCLUSIONS

Our present results suggest that the VB-SCCD technique has the capability in estimating location and extent of epileptic sources of interictal spikes and is promising to become a valuable noninvasive tool in assisting presurgical planning for partial epilepsy patients.

Clinical value of magnetoencephalographic spike propagation represented by spatiotemporal source analysis: correlation with surgical outcome

OBJECTIVE

To investigate the correlation between spike propagation represented by spatiotemporal source analysis of magnetoencephalographic (MEG) spikes and surgical outcome in patients with temporal lobe epilepsy.

METHODS

Thirty-seven patients were divided into mesial (n=27) and non-mesial (n=10) groups based on the presurgical evaluation. In each patient, ten ipsilateral spikes were averaged, and spatiotemporal source maps of the averaged spike were obtained by using minimum norm estimate. Regions of interest (ROIs) were created including temporoparietal, inferior frontal, mesial temporal, anterior and posterior part of the lateral temporal cortex. We extracted activation values from the source maps and the threshold was set at half of the maximum activation at the peak latency. The leading and propagated areas of the spike were defined as those ROIs with activation reaching the threshold at the earliest and at the peak latencies, respectively. Surgical outcome was assessed based on Engel's classification. Binary variables were created from leading areas (restricted to the anterior and mesial temporal ROIs or not) and from propagation areas (involving the temporoparietal ROI or not), and for surgical outcome (Class I or not). Fisher's exact test was used for significance testing.

RESULTS

In total and mesial group, restricted anterior/mesial temporal leading areas were correlated with Class I (p<0.05). Temporoparietal propagation was correlated with Class II-IV (p<0.05). For the non-mesial group, no significant relation was found.

CONCLUSIONS

Spike propagation patterns represented by spatiotemporal source analysis of MEG spikes may provide useful information for prognostic implication in presurgical evaluation of epilepsy.

Electroencephalographic inverse localization of brain activity in acute traumatic brain injury as a guide to surgery, monitoring and treatment

OBJECTIVE

To inverse-localize epileptiform cortical electrical activity recorded from severe traumatic brain injury (TBI) patients using electroencephalography (EEG).

METHODS

Three acute TBI cases were imaged using computed tomography (CT) and multimodal magnetic resonance imaging (MRI). Semi-automatic segmentation was performed to partition the complete TBI head into 25 distinct tissue types, including 6 tissue types accounting for pathology. Segmentations were employed to generate a finite element method model of the head, and EEG activity generators were modeled as dipolar currents distributed over the cortical surface.

RESULTS

We demonstrate anatomically faithful localization of EEG generators responsible for epileptiform discharges in severe TBI. By accounting for injury-related tissue conductivity changes, our work offers the most realistic implementation currently available for the inverse estimation of cortical activity in TBI.

CONCLUSION

Whereas standard localization techniques are available for electrical activity mapping in uninjured brains, they are rarely applied to acute TBI. Modern models of TBI-induced pathology can inform the localization of epileptogenic foci, improve surgical efficacy, contribute to the improvement of critical care monitoring and provide guidance for patient-tailored treatment. With approaches such as this, neurosurgeons and neurologists can study brain activity in acute TBI and obtain insights regarding injury effects upon brain metabolism and clinical outcome.

Minimum norm estimates in MEG can delineate the onset of interictal epileptic discharges: A comparison with ECoG findings

The analysis of epileptic discharges in magnetoencephalography with minimum norm estimates (MNE) is expected to provide more precise localization of epileptic discharges compared with electroencephalographic estimations. However, the clinical feasibility of MNE remains unclear. In this study, we aimed to elucidate the onset and propagation patterns of interictal spikes using MNE. Seven patients with intractable epilepsy whose epileptogenicity was assumed to exist in the convexity of the cerebral cortex were studied. For MNE and electrocorticography (ECoG), we characterized the propagation patterns of interictal epileptic discharges according to the area in which they originated and where they extended; we then examined whether the propagation patterns observed in MNE were identified by ECoG. We also examined the relationship between the positions of spikes estimated by the equivalent current dipole (ECD) method and MNE. Among the seven patients, nine propagation patterns of epileptic discharges were observed by MNE, all of which were also identified by ECoG. In seven patterns, the epileptic activity propagated around the initial portion. However, in two patterns, the center of activities moved according to propagation with maintained activity of the initial portion. The locations of spikes identified by the ECD method were within the areas estimated by MNE when the epileptic activity propagated. However, the ECD method failed to detect onset activities identified by MNE in three of nine patterns. Thus, MNE is more useful as a means of presurgical evaluation for epilepsy than the ECD method because it can delineate the onset of epileptic activities as shown in ECoG.

Language lateralization represented by spatiotemporal mapping of magnetoencephalography

BACKGROUND AND PURPOSE

Determination of hemispheric language dominance is critical for planning epilepsy surgery. We assess the usefulness of spatiotemporal source analysis of magnetoencephalography for determining language laterality.

MATERIALS AND METHODS

Thirty-five patients with epilepsy were studied. The patients performed a semantic word-processing task during MEG recording. Epochs containing language-related neuromagnetic activity were averaged after preprocessing. The averaged data between 250 and 550 ms after stimulus were analyzed by using dynamic statistical parametric mapping. ROIs were obtained in the opercular and triangular parts of the inferior frontal gyrus, superior temporal gyrus, and supramarginal gyrus in both hemispheres. We calculated laterality indices according to 1) dSPM-amplitude method, based on the amplitude of activation in the ROIs, and 2) dSPM-counting method, based on the number of unit dipoles with activation over a threshold in the ROIs. The threshold was determined as half of the maximum value in all ROIs for each patient. A LI ≥0.10 or ≤-0.10 was considered left- or right-hemisphere dominance, respectively; a LI between -0.10 and 0.10 was considered bilateral. All patients underwent an intracarotid amobarbital procedure as part of presurgical evaluation.

RESULTS

The dSPM-counting method demonstrated laterality consistent with the IAP in 32 of 35 patients (91.4%), the remaining 3 (8.6%) demonstrated bilateral language representation, whereas the dSPM-amplitude method showed 18 (51.4%) concordant and 17 (48.6%) bilateral. No laterality opposite to the IAP was found.

CONCLUSIONS

Spatiotemporal mapping of language lateralization with the dSPM-counting method may reduce the necessity for an IAP in as many as 90% of patients.

Ictal MEG onset source localization compared to intracranial EEG and outcome: improved epilepsy presurgical evaluation in pediatrics

PURPOSE

Magnetoencephalography (MEG) has been shown a useful diagnostic tool for presurgical evaluation of pediatric medically intractable partial epilepsy as MEG source localization has been shown to improve the likelihood of seizure onset zone (SOZ) sampling during subsequent evaluation with intracranial EEG (ICEEG). We investigated whether ictal MEG onset source localization further improves results of interictal MEG in defining the SOZ.

METHODS

We identified 20 pediatric patients with one habitual seizure during MEG recordings between October 2007 and April 2011. MEG was recorded with sampling rates of 600Hz and 4000Hz for 10 and 2min respectively. Continuous head localization (CHL) was applied. Source localization analyses were applied using multiple algorithms, both at the beginning of ictal onset and for interictal MEG discharges. Ictal MEG onsets were identified by visual inspection and power spectrum using short-time Fourier transform (STFT). Source localizations were compared with ICEEG, surgical procedure and outcome.

KEY FINDINGS

Eight patients met all inclusion criteria. Five of the 8 patients (63%) had concordant ictal MEG onset source localization and interictal MEG discharge source localizations in the same lobe, but the source of ictal MEG onset was closer to the SOZ defined by ICEEG.

SIGNIFICANCE

Although the capture of seizures during MEG recording is challenging, the source localization for ictal MEG onset proved to be a useful tool for presurgical evaluation in our pediatric population with medically intractable epilepsy.

Different cortical involvement pattern of generalized and localized spasms: a magnetoencephalography study

We report successful magnetoencephalography (MEG) recording in a child who had generalized epileptic spasms (ESs) as well as ESs involving the legs only during the recording. MEG source localization results demonstrated that (1) the interictal epileptiform discharges and both types of ESs had the same origin, that is, the right parietal region, and (2) the two types of ESs had different cortical spread patterns, that is, epileptic involvement localized to the right parietal region in spasms of the legs and rapid diffuse involvement in generalized spasms. In this case, MEG provided new insight into the mechanisms underlying the two types of ESs: both types were generated from the same focus, and in generalized ESs, abnormal excitation spread to cortical areas diffusely.

Concepts of connectivity and human epileptic activity

This review attempts to place the concept of connectivity from increasingly sophisticated neuroimaging data analysis methodologies within the field of epilepsy research. We introduce the more principled connectivity terminology developed recently in neuroimaging and review some of the key concepts related to the characterization of propagation of epileptic activity using what may be called traditional correlation-based studies based on EEG. We then show how essentially similar methodologies, and more recently models addressing causality, have been used to characterize whole-brain and regional networks using functional MRI data. Following a discussion of our current understanding of the neuronal system aspects of the onset and propagation of epileptic discharges and seizures, we discuss the most advanced and ambitious framework to attempt to fully characterize epileptic networks based on neuroimaging data.

Source localization of ictal epileptic activity investigated by high resolution EEG and validated by SEEG

High resolution electroencephalography (HR-EEG) combined with source localization methods has mainly been used to study interictal spikes and there have been few studies comparing source localization of scalp ictal patterns with depth EEG. To address this issue, 10 patients with four different scalp ictal patterns (ictal spikes, rhythmic activity, paroxysmal fast activity, obscured) were investigated by both HR-EEG and stereoelectroencephalography (SEEG). Sixty-four scalp-EEG sensors and a sampling rate of 1kHz were used to record scalp ictal patterns. Five different source models (moving dipole, rotating dipole, MUSIC, LORETA, and sLORETA) were used in order to perform source localization. Seven to 10 intracerebral electrodes were implanted during SEEG investigations. For each source model, the concordance between ictal source localization and epileptogenic zone defined by SEEG was assessed. Results were considered to agree if they localized in the same sublobar area as defined by a trained epileptologist. Across the study population, the best concordance between source localization methods and SEEG (9/10) was obtained with equivalent current dipole modeling. MUSIC and LORETA had a concordance of 7/10 whereas sLORETA had a concordance of only 5/10. Four of our patients classified into different groups (ictal spikes, paroxysmal fast activity, obscured) had complete concordance between source localization methods and SEEG. A high signal to noise ratio, a short time window of analysis (<1s) and bandpass filtering around the frequency of rhythmic activity allowed improvement of the source localization results. A high level of agreement between source localization methods and SEEG can be obtained for ictal spike patterns and for scalp-EEG paroxysmal fact activities whereas scalp rhythmic discharges can be accurately localized but originated from seizure propagation network.

Propagation of epileptic spikes reconstructed from spatiotemporal magnetoencephalographic and electroencephalographic source analysis

The purpose of this study is to assess the accuracy of spatiotemporal source analysis of magnetoencephalography (MEG) and scalp electroencephalography (EEG) for representing the propagation of frontotemporal spikes in patients with partial epilepsy. This study focuses on frontotemporal spikes, which are typically characterized by a preceding anterior temporal peak followed by an ipsilateral inferior frontal peak. Ten patients with frontotemporal spikes on MEG/EEG were studied. We analyzed the propagation of temporal to frontal epileptic spikes on both MEG and EEG independently by using a cortically constrained minimum norm estimate (MNE). Spatiotemporal source distribution of each spike was obtained on the cortical surface derived from the patient's MRI. All patients underwent an extraoperative intracranial EEG (IEEG) recording covering temporal and frontal lobes after presurgical evaluation. We extracted source waveforms of MEG and EEG from the source distribution of interictal spikes at the sites corresponding to the location of intracranial electrodes. The time differences of the ipsilateral temporal and frontal peaks as obtained by MEG, EEG and IEEG were statistically compared in each patient. In all patients, MEG and IEEG showed similar time differences between temporal and frontal peaks. The time differences of EEG spikes were significantly smaller than those of IEEG in nine of ten patients. Spatiotemporal analysis of MEG spikes models the time course of frontotemporal spikes as observed on IEEG more adequately than EEG in our patients. Spatiotemporal source analysis may be useful for planning epilepsy surgery, by predicting the pattern of IEEG spikes.

New directions in clinical imaging of cortical dysplasias

Neuroimaging is essential in the work-up of patients with intractable epilepsy. In pediatric patients with medically refractory epilepsy, cortical dysplasias account for a large percentage of the epileptogenic substrate. Unfortunately, these are also the most subtle lesions to identify. For this reason, there has been ongoing interest in utilizing new advanced magnetic resonance imaging (MRI) techniques to improve the ability to identify, diagnose, characterize, and delineate cortical dysplasias. Technologic gains such as multichannel coils (32 phased array and beyond) and higher field strengths (3T, 7T, and greater) coupled with newer imaging sequences such as arterial spin labeling (ASL), susceptibility weighted imaging (SWI) and diffusion tensor/spectrum imaging (DTI/DSI) are likely to increase yield. Improved MRI techniques coupled with a multimodality approach including magnetoencephalography (MEG), positron emission tomography (PET), and other techniques will increase sensitivity and specificity for identifying cortical dysplasias.