Magnetoencephalography interictal spike clustering in relation with surgical outcome of cortical dysplasia.

OBJECTIVE: The aim of this study was to devise an objective clustering method for magnetoencephalography (MEG) interictal spike sources, and to identify the prognostic value of the new clustering method in adult epilepsy patients with cortical dysplasia (CD).
METHODS: We retrospectively analyzed 25 adult patients with histologically proven CD, who underwent MEG examination and surgical resection for intractable epilepsy. The mean postoperative follow-up period was 3.1 years. A hierarchical clustering method was adopted for MEG interictal spike source clustering. Clustered sources were then tested for their prognostic value toward surgical outcome.
RESULTS: Postoperative seizure outcome was Engel class I in 6 (24%), class II in 3 (12%), class III in 12 (48%), and class IV in 4 (16%) patients. With respect to MEG spike clustering, 12 of 25 (48%) patients showed 1 cluster, 2 (8%) showed 2 or more clusters within the same lobe, 10 (40%) showed 2 or more clusters in a different lobe, and 1 (4%) patient had only scattered spikes with no clustering. Patients who showed focal clustering achieved better surgical outcome than distributed cases (p=0.017).
CONCLUSION: This is the first study that introduces an objective method to classify the distribution of MEG interictal spike sources. By using a hierarchical clustering method, we found that the presence of focal clustered spikes predicts a better postoperative outcome in epilepsy patients with CD.

INTRODUCTION

Cortical dysplasia (CD) is the second most common pathologic entity in surgically treated epilepsy4,8,10,11,13,17). However, the postoperative seizure outcome is relatively limited when compared outcome with epilepsy associated with hippocampal sclerosis or tumor4,7,18). Recently, magnetoencephalography (MEG) has emerged as a useful tool for localization of the epileptogenic zone. MEG studies, which have attempted to determine the relationship between interictal spike source patterns and postoperative seizure outcome, have mostly dealt with children and their patient groups were mixed with several different pathologies9,15,25). Moreover, those studies used a subjective clustering method to determine clusters; none of the studies used a statistically based clustering method. One MEG study recruited children who had a single pathology of CD25). However, because they were interested in influences of complete or partial resection of areas containing MEG spike sources on surgical outcome, they did not report an association of clustering pattern with surgical outcome. There is one MEG study with adult epilepsy patients associated with CD, but that study included only four patients14). Therefore, based on the above studies alone, it is difficult to determine which features of the MEG spike source have prognostic value in CD-associated epilepsy surgery.

In this study, we developed an objective clustering method to classify MEG spike source patterns, and we determined the prognostic value of that new clustering method in adult epilepsy patients with CD. To that end, we undertook a retrospective analysis of MEG data with respect to postsurgical outcome.

MATERIALS AND METHODS

Patient population

Initially included were 58 patients with intractable epilepsy with histologically proven CD who underwent MEG examination and surgical operation between 2005 and 2010 at our hospital. Of that total, 15 patients with less than a 1 year follow-up, one patient younger than age of 15 years at surgery and 8 patients who had fewer than 5 MEG interictal spikes were excluded. Also excluded were 8 patients with medial temporal lobe epilepsy and one patient who had atrophy. As a result, 25 patients (mean age at surgery=29.3 years; SD=8.3; 11 females) were included in our retrospective study. The mean postoperative follow-up period was 3.1 years (SD=1.2 years). All patients provided informed consent. This study was approved by the Institutional Review Board of Seoul National University Hospital.

Histological classification

We classified focal CD (FCD) according to Blumcke's criteria2). FCD type I is characterized by radial (FCD type Ia) or tangential (FCD type Ib) dyslamination of the neocortex. FCD type II is characterized by cortical dyslamination and dysmorphic neurons without (type IIa) or with (type IIb) balloon cells. FCD type III was not included because, by definition, this type is associated with other pathologies.

Magnetic resonance imaging

All patients were examined using either a GE 1.5 or 3 T MRI unit (GE Horizon Echospeed) or a Siemens 1.5 T scanner (Siemens Avento system, Erlangen, Germany) both before and after surgery. The images were reconstructed offline and transferred to a UNIX workstation (HP C3750).

MEG acquisition & spike source localization

Cortical magnetic activity and electroencephalographic (EEG) activities were recorded simultaneously using a whole-head MEG system (VectorView™, Elekta Neuromag Oy, Helsinki, Finland) consisting of 306 channels arranged in triplets of two planar gradiometers and one magnetometer. Spontaneous MEG/EEG signals were acquired about 60 min (four 15 min epochs of data) while the patients were lying in a supine position with eyes closed in a magnetically shielded room. The sampling rate was 600 Hz with a bandpass filter of 0.1-200 Hz. A bipolar electrooculogram and an electrocardiogram were simultaneously recorded in order to monitor eye movement and cardiac artifacts. The location of the subject's head in the MEG sensor space was determined by measuring the magnetic signals from four head position indicator coils placed on the scalp. The coil locations with respect to the three anatomical landmarks, the nasion and two preauricular points, were identified by a three dimensional digitizer (FASTRAK™, Polhemus, Colchester, VT, USA). The obtained information was used to co-register the MRI, MEG, and head coordinates so that the sources of the MEG signals could be superim-posed on the subjects' MR image. After the recording, a signal-space-separation method was applied in order to reduce environmental and biological noises21,22). Each MEG interictal spike was visually inspected and localized using the single equivalent current dipole (ECD) method in Neuromag software (Elekta Neuromag, Helsinki, Finland). The ECDs were identified by a least-squares search in a spherical volume conductor model by using subsets of 20-50 channels surrounding the maximum signals. Only ECDs with a goodness-of-fit value >85% and confidence volume <3 mm3 were accepted. Accepted ECDs then were superimposed on the patients' individual MR image.

MEG spike clustering

Hierarchical clustering analysis was adopted for MEG spike source clustering using SPSS 19.0 software (SPSS, Chicago, IL, USA). Euclidean distances between ECDs in each patient were calculated as input values. The hierarchical grouping starts by defining a number of clusters equal to the total number of ECDs, N. In the first grouping step the number of clusters is reduced to N-1 by combining the two closest ECDs using a complete linkage method. In successive stages the newly formed groups are linked to other ECDs until the entire ECD population is finally linked together into one cluster. To make an appropriate number of cluster groups, a dendrogram tree, which shows the result of the hierarchical clustering analysis, was cut by using a two-step criteria. First, if the ECDs projected onto the MR image were shown as clearly different clustered locations, we divided the dendrogram where the same number of clusters are grouped (Fig. 1). Second, we calculated the standard deviation (SD) of the Euclidian distances of ECDs within the same cluster. The dendrogram tree was then cut where any pair of ECDs distances within the same cluster did not exceed the 2.5 SD. The 1 SD should not exceed 10 mm. More than 3 spikes in one group was considered a cluster.

Fig. 1

Hierarchical dendrogram (A) and equivalent current dipoles on MR images (B) for Patient #10. The individual spike sources (n=20) are represented along with the vertical axis. The ⊃-shaped lines connect the spike objects that are combined at each aggregation stage.

Surgical procedures and postsurgical seizure outcome

Extraoperative intracranial recording was done for all 25 patients to delineate the ictal-onset zone and to obtain functional mapping. The extent of invasive intracranial electrode implant was determined as follows4). If MRI showed a visible lesion, invasive intracranial electrodes were placed around the lesion. If MRI showed no visible lesion, invasive intracranial electrodes were placed to encompass the ictal-onset zone delineated by scalp video-EEG monitoring. When noninvasive results were discordant, invasive intracranial electrodes were implanted bilaterally as needed. After the invasive study, video-EEG monitoring was done until at least three typical seizures recorded. Resection included the whole extent of the visible lesion, if present, or the whole ictal-onset zone bound by sulci or functional cortex. In addition, underlying white matter was resected, if feasible. The surgical outcome was classified according to modified Engel criteria5).

Statistical analysis

All analysis was performed by using SPSS 19.0 software (SPSS). We used the χ2 test or, where appropriate, Fisher's exact test. A value of p<0.05 was considered to be statistically significant.

RESULTS

Postoperative seizure outcome was Engel class I in 6 (24%) patients, class II in 3 (12%), class III in 12 (48%), and class IV in 4 (16%) patients. Table 1 presents a summary of the patients' profiles.

Table 1

Summary of patient profiles

Rt : right, Lt : left, T : temporal, P : parietal, F : frontal, O : occipital, WNL : within normal limit, FCD : focal cortical dysplasia, Y : Clustered MEG spikes were resected, N : Clustered MEG spikes were not resected

MEG

Twelve of 25 (48%) patients showed 1 cluster, 8 (32%) showed 2 clusters, 3 (12%) showed 3 clusters, 1 (4%) showed 5 clusters, and 1 (4%) had scattered spikes with no clustering. Of the patients who showed 2 or more clusters, 2 patients' clusters were located within the same lobe and 10 patients' clusters were in different lobes. We defined the presence of 1 cluster and ≥2 clusters within the same lobe as focal clustering and ≥2 cluster in different lobes and scattered spikes as distributed clustering. Patients with focal clustering had a better surgical outcome than those with distributed clustering (p=0.017, Engel I vs. II-IV) (Table 2).

Table 2

Surgical outcome with respect to MEG cluster, CD type, resected lobe, and MR findings

MEG : magnetoencephalography, CD : cortical dysplasia

We visually checked whether the MEG clustered area was resected regardless of its extent (complete or partial) by projecting the ECDs on postoperative MR images. The area that contained a MEG cluster was removed in 5 of 6 (83%) cases in the Engel class I group, 2 of 3 (67%) in Engel II, 5 of 12 (42%) in Engel III, and 1 of 4 (25%) in the Engel IV group. Patients who had resection of a MEG cluster tended to show a better surgical outcome (p=0.097) (Table 1).

Histological subtypes and resection location

Of the 25 patients, 19 (76%) had type I and 6 (24%) had type II CD (Table 2). With respect to postoperative seizure outcome, there were no statistically significant differences were found between patients having type I CD and type II CD (p=0.074). Of the 25 patients, 14 (56%) underwent frontal lobe resection, 8 (32%) temporal lobe resection, and 3 (12%) occipital lobe resection (Table 2). There was no postoperative seizure outcome difference among patients with different lobe resection areas (p=0.278).

MRI findings

Nine (36%) of 25 patients showed a visible lesion of focal CD on their MR images (Table 2). No statistical differences were found between normal and abnormal MRI findings with respect to surgical outcome (p=0.363).

DISCUSSION

In this study, we developed an objective method of MEG interictal spike source classification, and identified the prognostic values of MEG dipole clusters in adult epilepsy patients with CD. A hierarchical clustering method was adopted to classify MEG spike sources in 25 adult patients with histologically proven CD, who underwent surgical resection for intractable epilepsy and had at least 1 year of follow-up. The results of the MEG spike source clustering method were then tested for their prognostic value toward surgical outcome. The 14 (56%) patients in the focal clustering group achieved better outcome than the 11 (44%) in the distributed clustering group. No statistical differences were found according to the CD type and MRI findings with respect to surgical outcome.

Hierarchical clustering method

Some researchers have aimed at determining the relationship between MEG interictal spike source patterns and postoperative seizure outcome in their studies9,15,25). However, the fact that each researcher used different criteria in classification of spike sources to describe their results makes it difficult to evaluate the prognostic values of MEG in epilepsy. Moreover, all of these clustering methods are somewhat subjective in that none used a statistically based clustering method. The hierarchical clustering method which we adopted in this study is a widely used data analysis tool in many fields, including machine learning, data mining, bioinformatics, and image analysis. It is a simple and objective method compared to other clustering methods used by previous researchers. To adopt this method in MEG spike source clustering, we simply entered the Euclidean distance matrix between ECDs as an input value. Subsequently, the obtained dendrogram depicts the clustering results by linking the closest ECDs together. By using this clustering method, we could clearly separate sub-groups of either focal or distributed cluster according to the dendrogram's clustering patterns. If other researchers who are interested in MEG spike source clustering use this hierarchical clustering method, we can expect an increase in the comparability of results among studies in the future.

MEG cluster with respect to surgical outcome

Studies focused on MEG clustering with respect to surgical outcome have suggested that focal clustering of MEG spike sources and resection of that area may predict a favorable postoperative outcome. However, most of these studies have dealt with children or with mixed groups with various pathologies6,9,15,16,25). Although there is one study dealing with adult epilepsy patients associated with CD, only 4 patients were included in their study14). We studied an adult group of 25 patients with histologically proven CD. This is the largest MEG epilepsy study with the single pathology of CD in an adult group with more than a 1 year follow-up period.

As a result of the MEG spike source clustering method, 14 (56%) of our patients showed focal clustering and 11 (44%) showed distributed clustering. There are statistically significant differences in postoperative seizure outcome between these two groups : the focal clustering group achieved a better outcome than the distributed clustering group. In addition to clustering type, we also compared the clustered area with the surgically resected area as seen on postoperative MR images. The area that contained the MEG cluster was removed in 5 of 6 (83%) patients in the Engel I group, 2 of 3 (67%) in Engel II, 5 of 12 (42%) in Engel III, and 1 of 4 (25%) in the Engel IV group. Four of the 5 (80%) patients in Engel III had another cluster (other than the resected cluster) in a location remote from the actual resection area compared with only 1 of 7 (14%) patients with another cluster in the Engel I and II groups. One Engel IV patient (#25) is the only patient in our group who underwent multiple area resection. Additionally, he had an awake-lesionectomy with intra-operative sensory, motor, and language mapping to avoid injury to the eloquent cortex. Multiple area resection and hindrance of complete resection caused by the eloquent area were considered as a bad prognostic factor12). Failure to achieve a seizure free outcome despite resection of the focal clustered area in Patient #25 may result from these reasons.

On the basis of our results, we propose that the presence of focal clustered spike sources and resection of that area predict a favorable surgical outcome, whereas the presence of distributed clustering predicts a poor surgical outcome, in adult epilepsy patients with CD.

MEG cluster with non-lesional cases

To draw firm conclusions about the relationship between MEG clustering and surgical outcome, explanation of the patients who had poor outcome despite the presence of focal clustering is necessary. This is especially important for the patients who showed no abnormality on their MR images. Many studies agree on the prognostic value of MRI-visible lesions1,19,23,26), but not many studies have elucidated the prognostic factors of the patients with normal MRI findings, especially in relation with associated electrophysiological characteristics16). Of our patients who showed no abnormal features on their MR images, 3 patients showed a poor surgical outcome (Engel class IV) despite the presence of focal MEG spike clustering. One patient (#25) had an awake-lesionectomy with multiple resection of frontal lobe as stated previously, while two other patients (#22, #23) had 1 cluster located in the posterior part of the resection area. We suggest that if the area containing MEG spike sources were resected, our non-lesional, single clustered patients might achieve favorable seizure outcome. Further intensive study with non-lesional patients is needed to confirm this suggestion.

Histological subtypes and MR findings

With regard to postoperative seizure outcome and CD type, no statistically significant differences were found. Even though there is a debate about surgical outcome depending on CD type, recent studies have concluded that type II CD patients have better outcome than type I patients3,4,20,24). The discrepancy between previous studies and our results may due to the limited number of patients. Tendency to show better outcome in type II CD in our results support this idea (p=0.074). Even though many studies report better outcome in lesional cases, no surgical outcome differences were observed between patients with normal and abnormal MRI findings in our patient cohort. This could be explained by the fact that the small number of patients having abnormality on MR images in our study makes it difficult to reveal statistically significant results. Only 9 of 25 patients showed lesions on preoperative MR images. Furthermore, 5 of 6 who achieved Engel class II-IV had divergent presurgical evaluation results, including PET, ictal SPECT, video-EEG monitoring, invasive EEG, MEG, and MRI. Previous studies have reported that divergent presurgical evaluation results predict poor outcome13,16). Thus, the limited number of patients and the presence of divergent presurgical evaluation results may play roles in the lack of improved outcome in the lesional cases in our study.

Limitation

Its retrospective design in a single institute is a major caveat which should be taken into account in exploiting the results of the present study. Nevertheless, in terms of the prognostic values of the MEG interictal spike cluster, this is the largest MEG study that has dealt with adult epilepsy group who had the single pathology of cortical dysplasia. In addition, this is the first study which introduces an objective method to classify the distribution of MEG interictal spike sources.

CONCLUSION

In this study, a hierarchical clustering method was adopted for MEG interictal spike source clustering in adult epilepsy patients with CD. By using this method of hierarchical clustering, we were able to delineate the value of using clustering to predict surgical outcome more clearly. The presence of focal clustered MEG interictal spike sources and subsequent resection of the clustered area predicted a favorable postoperative outcome. In contrast, the presence of a distributed cluster predicted a poor surgical outcome.