Brain (2002), 125, 640±655
Seizure anticipation in human neocortical partial epilepsy Vincent Navarro,1,2 Jacques Martinerie,1 Michel Le Van Quyen,1 SteÂphane Clemenceau,1,3 Claude Adam,1,2 Michel Baulac1,2 and ²Francisco Varela1 1Laboratoire
de Neurosciences Cognitives et Imagerie CeÂreÂbrale (LENA), CNRS UPR 640, 2Unite d'Epileptologie and 3Service de Neurochirurgie, HoÃpital de la PitieÂSalpeÃtrieÁre, Paris
Correspondence to:Professor Michel Baulac, Unite d'Epileptologie, Clinique Paul Castaigne, HoÃpital de la PitieÂ-SalpeÃtrieÁre, 47 boulevard de l'HoÃpital, 75651 Paris cedex 13, France E-mail:
[email protected] ²Deceased
Summary
The transition of brain activity towards an epileptic seizure is still a poorly understood phenomenon. Dynamic changes in brain activity have been detected several minutes before seizure emergence in populations of patients with mesial temporal lobe epilepsy (MTLE), using non-linear analysis of intracranial EEG recordings. Similar detection of a pre-ictal state has been obtained with standard scalp EEG recordings using a modi®ed non-linear method. Here we applied this strategy to the seizures of patients with neocortical partial epilepsy. Results obtained by non-linear similarity analysis of 41 seizures from 11 patients with refractory partial epilepsy originating from various sites of the neocortex showed that (i) a pre-ictal state was detected in 90% of the patients and in 83% of the seizures whatever their location, with a mean anticipation time of 7.5 min; (ii) similar pre-ictal dynamic changes were
May 28, 2001
detected when non-linear analysis methods were applied to either intracranial or scalp EEG recordings; (iii) the recording sites exhibiting these pre-ictal changes were distributed both within the epileptogenic focus and at remote locations; (iv) most pre-ictal dynamic changes were not correlated with linear changes in the frequency spectrum or with changes in the visually inspected EEG and the patients' behaviour. Hypotheses on the neuronal mechanisms underlying the pre-ictal period are discussed. The present results, together with those recently obtained in an MTLE population, suggest that changes in pre-ictal dynamics are a general phenomenon associated with seizure emergence in a wide population of patients with partial epilepsy, wherever the epileptogenic focus is located. The possibility of anticipating the onset of seizures has considerable therapeutic implications.
Keywords: neocortical epilepsy; scalp electroencephalogram; intracranial electroencephalogram; non-linear analysis; seizure anticipation Abbreviations: ANOVA = one-way analysis of variance; MTLE = mesial temporal lobe epilepsy
Introduction
For patients with intractable epilepsy, the unpredictability of seizure occurrence underlies an increased risk of sudden unexpected death (Cockerell et al., 1994) or morbidity (Buck et al., 1997) and is a major factor contributing to a poor quality of life (Devinsky et al., 1995). Anticipation of seizures would permit their consequences to be minimized and also allow therapeutic measures to be taken before seizures occur. Traditional approaches for detecting pre-ictal changes have encountered dif®culties. Clinical distinctions between premonitory sensations and the aura are delicate (Rajna et al., 1997), and if they are prodroma they may imply ã Guarantors of Brain 2002
very subtle and variable feelings even in the same patient. Only a few fortuitous (Baumgartner et al., 1998) or invasive (Weinand et al., 1997) studies using functional cerebral imaging have provided evidence for pre-ictal changes of regional cerebral blood ¯ow in patients with mesial temporal lobe epilepsy (MTLE). Traditional analyses of EEG signals have failed to detect speci®c changes preceding seizures, either by visual inspection of the recordings, even during intracerebral investigation, or using linear methods based on spectral analysis (Rogowski et al., 1981; Katz et al., 1991).
Table 1 Summary of clinical and surgical data of the patients Patient
Gender
Age at onset (years)
Age* (years)
Epileptic focus
Side
Brain MRI
Neuropathological ®ndings
Surgical procedure
Outcome (follow-up)
Frontal focus 1 (B.I.L.)
M
9
23
Pars orbitalis of F3
L
Normal
Astrocytic gliosis
IA (6 m)
2 (D.E.P.)
F
18
30
Anterior orbital gyrus
R
Normal
3 (M.A.Z.)
M
16
30
Middle frontal gyrus
L
Normal
Astrocytic gliosis and minor changes in cortical architecture Astrocytic gliosis
Basal±frontal external resection and F3 subpial transections Basal±frontal resection
30
45
T2
R
14
31
T1
R
External temporal scar Normal
Temporal neocortical focus 4 (Q.U.A.) M 5 (C.L.O.)
F
Temporo-occipital junction focus 6 (P.I.C.) F 2 months
Anterior frontal resection
IA (3 m)
Astrocytic gliosis
T2 cortical resection
IIB (3 y)
Astrocytic gliosis
T1 cortical resection
IVB (3.5 y)
Posterior TO parahippocampal cortical resection and amygdalohippocampectomy Partial TO cortical resection Lesionectomy
IIIA (4 y)
ID (7 m)
Lesionectomy
IA (5.5 y)
None (implication of multiple foci and functional area) Lesionectomy
IIB (1.5 y)
24
R
External TO lesion
Dysplasia
M
12
21
R
External TO scar
Astrocytic gliosis
8 (F.A.B.)
M
17
40
R
External TO lesion
9 (T.Y.R.)
M
18
23
L
Parieto-occipital focus 10 (P.I.N.) M
External TO lesion
Dysembryoplastic neuroepithelial tumour Meningioangiomatosis
13
18
L
Normal
24
26
R
Parietooccipital lesion
F
Meningioma
IB (1 y)
Surgical outcome was scored according to Engel, 1987. *At the time of the video-EEG recordings. M = male; F = female; y = years; m = months; R = right; L = left; TO = temporooccipital.
Anticipation of neocortical seizures
7 (E.L.G.)
11 (C.R.E.)
IA (3.5 y)
641
642
V. Navarro et al.
Strategies of non-linear analysis have succeeded in demonstrating pre-ictal changes in the intracerebral EEG in patients with MTLE (Lehnertz and Elger, 1998; Martinerie et al., 1998). Changes in EEG dynamics, extracted by different non-linear measures, were detected several minutes before the onset of the seizure, when patients reported no speci®c sensations and visual inspection of the recording showed no clear changes in the EEG signal (Martinerie et al., 1998). Non-linear methods were developed initially to describe the dynamics of complex physical systems with non-linear components, implying that their time course does not follow the linearity of the classical deterministic laws, but may instead exhibit non-proportional responses to speci®c inputs. Self-organizing behaviour and intermittency are other interesting properties of non-linear systems which permit transitions between states in the absence of external triggers. The epileptic process is also chronically intermittent and spontaneous, except for rare re¯ex epilepsies, and mathematical tools for the characterization of non-linear deterministic systems have been used successfully to determine changes in the state of brain activity before seizures (Schiff, 1998). Recently, we developed a new measure of non-linear analysisÐthe similarity indexÐwhich improves the robustness of the method, increases calculation speed and reduces sensitivity to EEG artefacts (Le Van Quyen et al., 1999). Applied to intracerebral recordings, this method allows accurate description of the spatiotemporal changes in the EEG before the seizure (Le Van Quyen et al., 2000). We have shown that this method can ef®ciently detect pre-ictal changes from the scalp EEG in spite of the blurring due to volume conduction and skull ®ltering (Le Van Quyen et al., 2001). These studies were mostly performed in patients with MTLE, who constitute the most homogeneous population of adults with partial epilepsy. The population of patients with neocortical partial epilepsies is much less homogeneous because of the diversity of clinical manifestations, underlying causes and electrographic characteristics (Lee et al., 2000). Neocortical partial epilepsy accounts for ~68% of the cases of partial epilepsy in adults (Semah et al., 1998). Patients with medically refractory neocortical epilepsy bene®t less from surgical treatment than those with MTLE and the postoperative prognosis is more variable (Engel, 1996; Spencer, 1998). This is due to dif®culties in the accurate localization of the epileptogenic focus, to the presence of multiple foci or to
the involvement of functional areas. The development of strategies to anticipate seizures could therefore improve the quality of life of these patients considerably. Non-linear analysis of interictal EEG recordings from patients with lesional neocortical epilepsy showed a relationship between the presumed epileptogenic area and the maximum values of a non-linear index (Widman et al., 2000). Neocortical seizures originating from the frontal lobe develop suddenly and propagate rapidly, in contrast to the slower progressive seizure emergence that occurs in MTLE, suggesting that their anticipation might be more dif®cult. Here we show, for the ®rst time, that non-linear analysis methods can detect preictal changes in EEG dynamics, from either scalp or intracerebral recordings, in patients with various types of neocortical epilepsy.
Methods Patient population
We studied 11 patients presenting medically intractable neocortical partial epilepsy (Table 1). They were selected from patients examined at the Epilepsy Unit of La PitieÂSalpeÃtrieÁre Hospital in Paris between 1993 and 2000 according to the following criteria: (i) neocortical origin of the epilepsy demonstrated by intracranial EEG recordings (n = 10) or scalp EEG recording (n = 1) in agreement with clinical and imaging data; (ii) spontaneous seizures, excluding seizures preceded by classical precipitating factors (hyperventilation and photic stimulation); (iii) a stable state of vigilance during the recording periods. EEG recordings from patients who were continuously awake (n = 22) or asleep (n = 19), on the basis of videotape recordings, were used for analysis to avoid perturbations due to changes in the clinical state. All patients required long-term continuous scalp EEG recordings and video monitoring using 21 or 27 electrodes, placed according to the extended International 10±20 System. For 10 patients, intracranial recordings were then necessary to localize accurately the structures generating seizure onset. Intracerebral electrodes with multiple contacts explored the lobe suspected of containing the epileptogenic focus and adjacent lobes. Subdural strips were added in some patients to sample the lateral or inferior cortices of the temporal or occipital lobe. Subdural grids with 36 contacts were used in two patients. The post-implantation location of the electrodes was determined by MRI. All patients gave informed consent
Fig. 1 Analysis of a scalp EEG recording from a patient with refractory partial epilepsy (seizure A of Patient 10). Visual inspection of the signal (A) did not show any modi®cation before the seizure. Signals from the seven most relevant channels are depicted during the beginning of the recording (1), during the few minutes before the seizure (2) and during the seizure (3). Spectral analysis with the fast Fourier transform did not reveal changes before the seizure. Frequency components from the signal of the Fz channel are shown (B). In contrast, non-linear analysis of the signal demonstrated changes in EEG dynamics before the seizure. The time course of the similarity index (C) for the 21 channels showed progressive deviations until the seizure. The statistical signi®cance of these deviations (D) from the reference state is represented for each channel by a colour scale (changes in the blue colour of the time window indicate signi®cance levels up to 5 SD). The anticipation time, determined for the Fz channel, was 15 min. Topographic projections of interpolated similarity indices on a two-dimensional circular head (E) are depicted for a time window during the interictal state (1), the pre-ictal state (2) and the seizure (3). The ®rst channel to reveal the seizure onset is represented by a yellow star.
Anticipation of neocortical seizures
643
644
V. Navarro et al.
to participation in this study, which was approved by the Consultative Committee for the Protection of Patients in Biomedical Research of the PitieÂ-SalpeÃtrieÁre Hospital, Paris. The video-EEG recordings were performed with a 32channel BMSI system (Nicolet-BMSI, Madison, Wis., USA). The raw data were digitized at 200 Hz with 12-bit resolution, and were passed to a channel ampli®er system with band-pass ®lter settings of 0.5±99 Hz using an external reference over linked ears. Our study was based on 41 seizures extracted from long-term EEG recordings of 60±120 min duration that included an interval of 50±90 min preceding the seizure. For all but three patients, both scalp and intracranial recordings were available for analysis. Seizure onset was de®ned as the time of the ®rst clinical symptoms (including aura sensation) or the time of the ®rst signi®cant electrographic changes (burst of spikes, sinusoidal waves or low-voltage fast activity) if they occurred earlier. Additionally, long interictal intracranial EEG periods (n = 8) of 90±450 min without seizure were analysed for ®ve patients in order to study the behaviour of the non-linear indexes far from a seizure in a stable state of vigilance.
reference period was taken from the start of the recording; it preceded the seizure by at least 50 min and had a duration of 300 s. This period contains different features of the interictal activity, possibly including artefacts and spikes, which were taken into account. The reference window started at the same time for all channels. In order to allow comparison of dynamics between this long reference window and shorter test windows (30 s), we identi®ed the most frequent occupations of the phase space ¯ow by random selection of a subset of points in the reference dynamics reconstruction. The corresponding dynamics of the two EEG windows were then compared using the normalized cross-correlation sum (Manuca and Savit, 1996; Schreiber and Schmitz, 1997) at a radius scale chosen at 30% of the cumulative neighbourhood distribution of the reference set. We used a similarity index, ranging from 0 to 1, to provide a sensitive measure of closeness of dynamics in the two EEG windows; values close to 1 indicate the greatest degree of similarity in the dynamics. The time course of this index shows how the long-term dynamics of the signals from each EEG channel change before seizure onset.
Non-linear EEG analysis
Statistical signi®cance of the pre-ictal deviations of the similarity index
Mathematical details of the non-linear method have been described previously (Le Van Quyen et al., 1999). We present here the main steps of the strategy, which was carried out independently for each EEG channel after digital ®ltering (1± 98 Hz for intracerebral recordings and 1±48 Hz for scalp recordings).
Reconstruction of EEG dynamics using time intervals
The dynamics of the EEG signal was evaluated for successive windows of 30 s duration. One novel feature of this method was to reduce the EEG signal to pure phase information. For this purpose, we used the sequence of time intervals between positive-going crossings of the signal through the electrical zero. The EEG dynamics was reconstructed by the delay method with m-dimensional embedding (m is the number of dimensions of the reconstructed dynamics), using m consecutive time intervals obtained by the zero-crossing method (m = 16) (see Fig. 1 in Le Van Quyen et al., 2001). Then, in order to retain the main components of the dynamics and to reduce the noise level, the embedding dimension was reduced to m = 4 using a singular-value decomposition (Albano et al., 1988).
Measurement of similarity between dynamics from a reference state and from test EEG windows
The dynamics corresponding to each successive EEG window were compared with that of a reference window. A long
We assessed the statistical signi®cance of pre-ictal changes by quantifying deviations from the reference state. The signi®cance S of the deviation was de®ned by the ratio (g ± m)/s (m and s correspond to the mean and standard deviation of similarity variations during the reference period), whose P value was given by the Chebyshev inequality (for any statistical distribution of g, P(|S| > k) < 1/k2, where k is the chosen statistical threshold). The temporal evolution of the statistical signi®cance S was used to detect the instant at which the similarity index reached a critical level. We characterized this instant as anticipation time when the ratio S reached a critical level of k = 5 (P = 0.04 for the bilateral test) and remained at or above this threshold until the seizure.
Traditional EEG analysis Visual inspection of the EEG recordings
All recordings were examined to de®ne the electrographic onset of the seizure accurately. We also searched for changes in EEG activity of any channel (including modi®cations of background rhythms or of interictal epileptic activity, or the presence of sustained artefacts) in comparison with activity during the reference period used for the non-linear analysis (Table 2).
Spectral analysis
Fast Fourier transforms were performed on successive EEG windows of 10 s duration to analyse the power spectrum of
Table 2 Clinical and electrographic characteristics of the ictal, interictal and pre-ictal states Patient/ Level of EEG Ictal state seizure recording Electrographic pattern of seizure onset
Initial Type Seizure Interictal EEG activities frequency of duration of seizure seizure (s)
Clinical state
1A
Scalp
13 Hz
SG
107
Sleeping 3.5
No
No
1B
Intracranial
Rhythmic alpha activity LVFA
45 Hz
SG
87
Sleeping 27
Slowing rhythm (VI, FFT)
No
As above As above LVFA
45 Hz 45 Hz 30 Hz
SP CP CP
36 33 45
Sleeping 11.5 Sleeping 5 Awake 2
As above As above No
No No No
As above As above As above HABS
30 30 30 20
CP SG CP P
24 99 23 29
Awake Sleeping Sleeping Sleeping
No
No No
2B
HABS
25 Hz
I
31
As above
Sleeping 3.5
2C 2D
HABS HABS
20 Hz 20 Hz
CP P
54 90
As above As above and brief LVFA in the focus (35±45 Hz)
Sleeping 1.5 Sleeping 13
2E
HABS
20 Hz
CP
94
2 foci of S-W (F and T), synchronized or not
Sleeping 0
3 min before the seizure, intermittent LVFA then brief 8 Hz activity (VI) 5 min before seizure, intermittent and brief LVFA on focus channels (VI); slowing rhythm on anticipation channel (VI, FFT) No 20 min before seizure, slowing rhythm (VI, FFT), more spikes (VI) No
Questioned 43 s before seizure No
1C 1D 1E 1 1 1 2
F G H A
Subdural grid
Intracranial
Interictal state
Hz Hz Hz Hz
Scalp
Beta then theta
4 Hz
CP
32
3B
Intracranial
LVFA
45 Hz
P
22
LVFA
45 Hz
P
9
Semirhythmic slow waves Rhythmic sharp waves (T2±T3 expression) Rhythmic S-W (H expression)
3 Hz
CP
55
3 Hz
CP
137
8 Hz
CP
88
3C 4A
Scalp
4B
Intracranial
4C
Frequent S-W in most frontal channels and brief LVFA Intermittent slow waves and sharp waves in some frontal contacts Frequent sharp waves close to focus Frequent spikes in hippocampus and neocortical temporal lobe As above
Awake
Anticipation Corresponding time (min) EEG changes
5.5 0 0 1.5
1
Sleeping 21
Corresponding clinical changes
No
No No
Sleeping 6.5
No
No
Awake
3
Slowing rhythm (FFT) on anticipation channel
No
Awake
16
Enhancement of spiking rhythmicity in hippocampal contacts before seizure (VI)
No
No
645
Sleeping 13
Disappearance of slow rhythms (VI, FFT) As above
Anticipation of neocortical seizures
3A
Spikes in left frontotemporal lobes Permanent S-W in the focus sometimes associated with brief LVFA As above As above Continuous poly-S-W in the focus As above As above As above 2 foci of S-W (F and T), synchronized or not
Pre-ictal state
646
Table 2 Continued
Initial Type Seizure Interictal EEG activities frequency of duration of seizure seizure (s)
Clinical state
Anticipation Corresponding time (min) EEG changes
5A
HABS (T1 expression) HABS (H expression) Semirhythmic spikes Semirhythmic spikes
25 Hz
CP
86
Awake
0
24 Hz
SP
26
2 independent foci of S-W (H and external T1) As above
Awake
1.5
10 Hz
CP
38
Few spikes close to focus
Awake
10 Hz
CP
62
Rhythmic spikes
12 Hz
CP
60
Intracranial
5B 6A
Scalp
6B
Intracranial
Interictal state
Pre-ictal state
No
4
More fast rhythms (10±35 Hz) (FFT) No
Awake
2
No
No
Awake
0
Awake
5.5
No
No
Sleeping 22
No
No
Sleeping 3
No
No
Awake
1
No
Awake
3.5
No
10
Fewer fast rhythms (25± 50 Hz)
Questioning 1 min before Phoning 3.5 min before No
7A
Scalp
Theta activity
5 Hz
CP
75
7B
Intracranial
LVFA (lingual and cuneus gyri)
45 Hz
P
107
LVFA (lingual and cuneus gyri) LVFA (02±O3 occipital sulci) Theta activity
15 Hz
SP
40
Subcontinuous interictal activity in some focus channels (15-40 Hz), then spikes for 30 min before seizure Subcontinuous interictal activity (polyspikes in 1 focus channel and rhythmic S-W in 2 other channels) 2 foci of spikes (R and L FT), not synchronized Brief discharges in focus (20 Hz) often followed by spikes in and around focus; spikes in H As above
40 Hz
SP
91
Spikes in and around focus
4 Hz
CP
113
LVFA (O3 occipital and T3 temporal gyri) As above As above As above
23 Hz
CP
75
Spikes in H and focus
Awake
30 Hz 40 Hz 30 Hz
CP CP CP
71 118 60
Awake 7 Sleeping 4.5 Sleeping 9
Theta activity No involvement of parietal channels in the seizure
5 Hz
SP SP
64 0
Spikes in H Spikes in H and focus Spikes in H and focus and brief LVFA in focus
6C
7C 7D 8A
Scalp
8B
Intracranial
8C 8D 8E 9A 9B
Scalp Subdural grid on parietal lobe
Continuous spike activity in one channel
Corresponding clinical changes
Awake Awake
0 0
No
No No No
No No No
Fewer fast rhythms (15±50 Hz) during dynamical changes that stop before seizure
No
V. Navarro et al.
Patient/ Level of EEG Ictal state seizure recording Electrographic pattern of seizure onset
Table 2 Continued
Initial Type Seizure Interictal EEG activities frequency of duration of seizure seizure (s)
Clinical state
10 A 10 B
Rhythmic spikes Poly-S-W (H activity 30 s after ®rst clinical signs) As above Rhythmic spikes (H activity 15 s after ®rst clinical signs) LVFA (precuneal activity 10 s before clinical signs)
7 Hz 3 Hz
CP SP
92 110
Spikes in R and L TO 2 foci in focus (spikes) and H (poly-S-W)
Sleeping 15 Awake 2
3 Hz 30 Hz
SP SP
99 88
As above Awake 2 foci of spikes in focus and Awake mesial or external temporal lobe
17 Hz
SP
100
Theta activity
8 Hz
CP
62
Scalp Intracranial
10 C 10 D
10 E
11 A
Scalp
Interictal state
Pre-ictal state Anticipation Corresponding time (min) EEG changes
Corresponding clinical changes
No No
No No
1 6
No No
No No
2 foci of spikes in the focus and Awake external temporal lobe, and brief LVFA (20 Hz)
13.5
Electronic game playing
Awake
12.5
Disappearance of fast rhythms (15±25 Hz) and more 5±10 Hz rhythms (FFT) in PT4; more fast rhythms (20±40 Hz) and fewer slow rhythms (1±15 Hz) in PC1 No
No
LVFA = low-voltage fast activity; HABS = high-amplitude beta spike activity; SG = secondary generalized seizure; SP = simple partial seizure; CP = complex partial seizure; P = unde®ned partial seizure occurring during sleep; I = infraclinical seizure; S-W = spike-wave; VI = visual inspection of EEG; FFT = fast Fourier transform of EEG signal; R = right; L = left; F = frontal lobe; T = temporal lobe; O = occipital lobe; H = hippocampus.
Anticipation of neocortical seizures
Patient/ Level of EEG Ictal state seizure recording Electrographic pattern of seizure onset
647
648
V. Navarro et al.
the channels that showed the most signi®cant pre-ictal changes according to the non-linear analysis.
Statistical analysis
Non-parametric tests were used to compare electroclinical data and anticipation times. The signi®cance of differences in anticipation times from scalp and intracranial EEG data or during wakefulness versus sleep was evaluated with the Mann±Whitney U-test. Kruskal±Wallis one-way analysis of variance (ANOVA) was used to compare anticipation times according to the types of clinical seizure. Spearman's coef®cient of rank correlation was determined to correlate anticipation times with the initial frequency or the duration of the epileptic discharges. Two-way ANOVA was used to study whether the location of the epileptogenic focus was correlated with the percentage of perifocal or distant EEG channels from which pre-ictal changes were detected. The criterion for statistical signi®cance was P < 0.05.
Results Representative examples of the spatiotemporal dynamics of the pre-ictal EEG changes
We ®rst present examples of the non-linear analysis of scalp (Fig. 1) and intracerebral EEG (Fig. 2) from the same patient (Patient 10), who had a left parieto-occipital focus with no underlying abnormality detected by MRI. Our analyses identi®ed a pre-ictal state preceding distinct seizures in both sets of EEG recordings. The scalp EEG recording contains a seizure that occurred during spontaneous sleep near the end of the night. The background activity showed continuous slow waves (~2.5 Hz), but no changes were apparent from visual
inspection of the recording before the seizure (Fig. 1A). The dynamics of the EEG signal from each of the 21 electrodes were analysed independently by a non-linear measure of similarity between a reference period and successive 30 s segments of the recording. The time course of the similarity indices for all channels showed high deviations from the baseline during the seizure and the post-ictal state (Fig. 1C). The EEG dynamics then returned to a state similar to the reference. Interestingly, similarity indices for several channels shifted to reduced levels for several minutes before the seizure. In order to study the spatiotemporal dynamics of these changes, the statistical signi®cance of these deviations from the baseline was calculated and dynamic changes up to 5 SD were represented for each of the channels (Fig. 1D). The statistical map demonstrated widespread pre-ictal changes in dynamics. We de®ned this period as the pre-ictal state and calculated an anticipation time from the electrode that showed the earliest of the dynamic changes that persisted until the seizure. The Fz channel allowed retrospective anticipation of seizure onset by 15 min. Intracranial EEG recordings were then performed in this patient using ®ve intracerebral electrodes to explore the left parieto-occipital junction and the temporal lobe. Intracranial EEG signals were analysed using the same non-linear method. During the recording shown in Fig. 2, the patient was watching television or playing electronic games. The time course of the similarity indices (Fig. 2C) from all 32 channels showed sharp decreases in similarity during the simple partial seizure, but decreases also occurred during two periods before the seizure (from 22 to 41 min, and from 43 to 56 min). The statistical map showed that these deviations were highly signi®cant for several channels (Fig. 2D). We measured the anticipation time from the second of the
Table 3 Summary of the spatial distribution of EEG channels with pre-ictal changes
Epileptogenic focus Distant
Number of channels
Number of channels showing pre-ictal changes
Percentage of channels showing pre-ictal changes
5.47 6 0.71 24.76 6 0.75
1.03 6 0.16 5.26 6 0.81
25.93 6 5.07 22.81 6 3.64
Data are mean 6 standard error of the mean. The percentage of channels showing pre-ictal changes was ®rst calculated for each seizure, because the total number of channels varied from 21 to 32.
Fig. 2 Analysis of an intracranial EEG recording from the same patient as in Fig. 1 (seizure E of Patient 10). Visual inspection of the signal (A) showed no changes before the seizure. Signals from the seven most relevant channels are depicted during the beginning of the recording (1), the few minutes before the seizure (2) and the seizure (3). Spectral analysis with the fast Fourier transform showed moderate changes before the seizure. A decrease of 1±15 Hz and an increase of 20±40 Hz were observed in the signal from the ®rst contact of the precuneus electrode (PC1) (B). Non-linear analysis of the signal clearly demonstrated changes in EEG dynamics preceding the seizure. The time course of the similarity index (C) for the 32 contacts showed two periods of deviation before the seizure. The statistical signi®cance of these deviations (D) from the reference state is shown for each channel using a colour scale (changes in the blue colour of time windows indicate signi®cance levels up to 5 SD). The anticipation time, determined for the PC1 channel, was 13.5 min. Channels with and without pre-ictal changes are depicted in red and blue, respectively, on a diagram of the locations of intracranial electrodes (E). The ®rst channel to show the seizure onset is represented by a yellow star. AT = anterior temporal lobe; PT = posterior temporal lobe; PC = precuneus gyrus; IC = infracalcarine; SC = supracalcarine.
Anticipation of neocortical seizures
649
650
V. Navarro et al.
Fig. 3 Percentage of channels with non-linear pre-ictal changes. Channels from the focus are shown in black and distant channels in white for each seizure for which changes in EEG dynamics were detected.
dynamic changes (43 to 56 min) because of the transient break between these periods of changed dynamics. The ®rst contact of the precuneus electrode (PC1) detected changes in dynamics 13.5 min before seizure onset. All channels with pre-ictal changes were involved in the initial spread of the seizure. They were either adjacent to the focus or in the neocortical part of the temporal lobe (Fig. 2E). No signi®cant change in dynamics occurred in the signal from the channel closest to the seizure onset (PC3). Visual inspection of the EEG revealed intermittent spikes for these channels without speci®c changes before the seizure (Fig. 2A), but fast Fourier transform analysis of the PC1 channel signal revealed that the frequency spectrum did change (Fig. 2B).
Detection of pre-ictal changes and anticipation time
Similar analyses were performed in 11 patients with medically refractory partial epilepsy whose presurgical investigations demonstrated a neocortical origin (Table 1). For eight patients, we studied several recordings that contained a seizure, from both the scalp and the intracerebral EEG (Table 2). In 10 patients, pre-ictal changes in EEG dynamics were detected before at least one seizure. Pre-ictal changes were found in 34 (83%) of the 41 seizures analysed. The mean percentage of anticipated seizures per patient was 79 6 9.5%. The mean anticipation time for these 34 seizures was 7.54 6 1.15 min.
We then studied factors that might in¯uence the detection of pre-ictal changes and the anticipation time (Table 2). First, we could detect no differences according to the location of the epileptogenic focus in the brain (Table 2). Secondly, pre-ictal changes were detected as ef®ciently from scalp EEG recordings as from intracerebral EEG recordings (eight out of 10 seizures were anticipated using scalp recordings and 26 out of 31 seizures were anticipated using intracerebral recordings) and the mean anticipation times were not statistically different (6.44 6 1.71 and 7.88 6 1.42 min for scalp and intracranial recordings, respectively). Thirdly, we examined how electroclinical data might in¯uence the variability with which pre-ictal changes were detected in the same patient. (i) The arousal state of the patient during the recording did not modify our ability to anticipate seizures (16 out of 19 seizures anticipated during sleep and 18 out of 22 seizures anticipated during wakefulness), and mean anticipation times in the two states were not statistically different (although a tendency was observed: 10 6 1.99 min during sleep and 5.33 6 1.08 min during wakefulness; P = 0.06). (ii) The type of seizure did not in¯uence seizure anticipation (eight out of 10 simple partial seizures, 18 out of 22 complex partial seizures and two out of three secondary generalized seizures were anticipated). Differences in anticipation times for the different types of seizure (2.6 6 0.73, 6.29 6 1.23, 17.25 6 2.2 and 15.25 6 11.75 min for simple, complex partial, inde®nite partial and secondarily generalized seizures, respectively) were not statistically signi®cant [Kruskal±
Anticipation of neocortical seizures
651
Fig. 4 Non-linear analysis of an intracranial EEG recording from the same patient as in Figs 1 and 2 (Patient 10) during an interictal period. Far from a seizure, the time course of the similarity index (A) for the 32 contacts is highly stable and shows only minor deviations. The statistical signi®cance of these deviations (B) from the reference state shows infrequent changes that are neither sustained nor simultaneous. A signi®cant change during more than two time windows occurred in only one of the 32 channels.
Wallis ANOVA; c2(3) = 2.94, P > 0.05]. (iii) The duration of seizures, de®ned electrographically, was not correlated with the anticipation time (Spearman's r = 0.05). However, we did observe a signi®cant correlation between the anticipation time and the initial frequency of the ictal discharge (Spearman's r = 0.36, P = 0.03), suggesting that the longer the anticipation time the faster the initial ictal frequency. We could identify several features of the minority of seizures (seven out of 41) which were not anticipated. First, we observed sudden and abrupt declines in the similarity indices at seizure onset, contrasting with the stability of the dynamics before onset (seizures G and H of Patient 1; seizure E of Patient 2). In one seizure (seizure A of Patient 5), similarity indices changed just before the seizure but not suf®ciently to reach our statistical threshold. In three cases, the dynamics of scalp and/or intracerebral EEG (seizures A and B of Patient 9; seizure C of Patient 6) changed
signi®cantly before the seizure but returned to subthreshold values several minutes before the seizure began.
Spatial distribution of the pre-ictal changes
Calculation of the non-linear similarity index independently for each channel permitted a spatiotemporal analysis of EEG dynamics. For seizures with a pre-ictal change in dynamics, the number of channels implicated in the pre-ictal state varied with a mean of 6.32 6 0.91. Pre-ictal changes were detected in fewer than three channels for seven seizures and in more than 12 channels for four seizures. We observed several spatiotemporal patterns of pre-ictal changes. We evaluated the spatial distribution of pre-ictal changes by comparing the percentage of channels in the epileptogenic focus (those channels implicated in the seizure onset) that showed preictal changes with the percentage at sites distant from the
652
V. Navarro et al.
focus (Table 3 and Fig. 3). Pre-ictal dynamical changes were detected in 25.93 6 5.07% of channels in the focus and 22.81 6 3.64% of the remote channels. Interestingly, the spatial distribution of the pre-ictal changes varied with the location of the focus (Fig. 3). For seizures originating in the temporal, parietal or occipital lobes, the higher proportion of channels showing pre-ictal changes was higher for channels in the focus (33.10 6 7.29%) than for those at a distance (15.87 6 2.85%), whereas the opposite was found for seizures of frontal origin (14.36 6 4.87 and 34.02 6 7.51%, respectively). Two-way ANOVA showed that this interaction between the site of the focus and the spatial distribution of the pre-ictal changes was statistically signi®cant [F(1,32) = 12.37, P = 0.001].
Electroclinical correlation with the pre-ictal changes
We evaluated temporal correlations between changes in EEG dynamics detected by the non-linear similarity index and modi®cations of three parameters: visual inspection of EEG recordings, the frequency spectrum of the EEG signal (a classical linear analysis) and the behaviour of the patient, as monitored on the video recordings (Table 2). No modi®cation of the EEG was observed before most seizures (n = 28). For the 13 other seizures, we detected concomitant changes of the EEG signal visually (n = 9) and/ or on examination of the frequency spectrum (n = 11). These linear changes were highly variable in terms of frequency range and dynamics (disappearance or emergence of one or several frequency bands) (Table 2). Among the awake patients (n = 22), behavioural changes sometimes occurred simultaneously with the non-linear changes in the EEG (n = 4). They involved responses to medical questions (n = 2), use of the telephone (n = 1) or playing electronic games (n = 1). Apart from the last of these behavioural changes, the correlations may have been coincidental. Changes in EEG dynamics in the game-playing patient began when he started playing, disappeared when he stopped and recommenced when he played again (Fig. 2C). Nevertheless, analysis of other video-EEG recordings from the same patient permitted us to exclude the use of the electronic game as a systematic factor precipitating the seizures.
Speci®city of the method
In order to study the speci®city of the changes in dynamics detected before the seizures, we analysed long EEG periods without a seizure. The time course of the similarity indices for all channels showed high stability of the non-linear dynamics when the vigilance state of the patient was stable (Fig. 4). Few statistically signi®cant changes occurred; they were very brief (mean duration 57.3 6 13.2 s) and largely restricted to a single channel, so that the average number of false positives
(de®ned as signi®cant changes lasting >1 min and concerning more than two channels simultaneously) was 0.31 6 0.20 per hour. On the other hand, when major changes in vigilance occurred, such as the transition from wakefulness to deep sleep, profound changes in EEG rhythms necessarily modi®ed measurements of the dynamics. If the reference state were to be taken during wakefulness, the similarity index would be expected to show deep and sustained deviations when the patient fell asleep. For this reason, we restricted our study to EEG periods before seizure in a stable state of vigilance.
Discussion Seizure anticipation in partial epilepsy
Our results demonstrate that, in addition to the interictal, ictal and post-ictal states, a pre-ictal state can be detected from the brain dynamics of patients with neocortical partial epilepsy. Analysis of EEG signals by non-linear similarity measures allowed hidden properties to be extracted from the electrical signal from the brain, showing that seizures are generally preceded by progressive changes in dynamics. This method permitted identi®cation of pre-ictal EEG changes in most seizures of patients with intractable partial epilepsy originating in the neocortex (34 out of 41 seizures) as well as the mesial temporal lobe (Le Van Quyen et al., 2000, 2001). Mean anticipation times were of the same order (7 min) for neocortical epilepsy and MTLE (Le Van Quyen et al., 2001).
Scalp and intracranial EEG analyses
Strikingly, our method possessed a similar ef®ciency when applied to recordings from either scalp or intracranial EEG. In the present study, we showed that the number of seizures anticipated and the anticipation times did not differ in analyses of scalp and later intracranial EEG data from the same patients. Similar results have been obtained recently from MTLE patients in comparisons of simultaneously recorded scalp and intracranial EEGs (Le Van Quyen et al., 2001). Possibly, the phase information contained in EEG signals, crucial to our similarity index, is transmitted faithfully to surface electrodes, whereas the amplitude of deep signals is considerably attenuated because of the interposition of the skull. The ability to identify a pre-ictal state by a standard, non-invasive route has numerous clinical applications (see below).
In¯uence of pathological parameters
Our ®nding that neither the location in the neocortex nor the aetiology of these neocortical epilepsies (developmental versus mature pathology, with or without MRI abnormalities)
Anticipation of neocortical seizures appears to in¯uence our capacity to detect pre-ictal changes indicates that the dynamic route towards the seizure is generic, even if the neuronal networks are affected differently (Spencer, 1998).
Meaning of the pre-ictal state
Two main hypotheses might explain the occurrence of dynamic changes before a seizure. The ®rst involves a `route' towards the seizure, corresponding to the progressive recruitment of neuronal activities until a critical mass has been reached, leading at this moment to the emergence of the seizure (Martinerie et al., 1998). The second hypothesis involves a `facilitating' state, corresponding to a change in brain activity that serves as a background from which a seizure could easily and suddenly appear (Le Van Quyen et al., 2001). It might consist of an impairment of the balance between neuronal excitability and its inhibitory controls (Lopes da Silva and Pijn, 1995). However, the identity of the neuronal networks involved (corticocortical and/or thalamocortical interactions) and the nature of the cellular modi®cations (synaptic or electrotonic interactions and/or changes in intrinsic cellular properties) remains to be uncovered.
Topological distribution of the pre-ictal changes in dynamics
The spatial dynamics of pre-ictal changes may permit discrimination between these two pathophysiological hypotheses. Pre-ictal changes predominant in channels in the epileptogenic focus might support the `route to seizure' hypothesis, as we observed for the patients with neocortical temporal, parietal and occipital lobe epilepsies. On the other hand, widespread pre-ictal changes remote from the focus are suggestive of the `facilitating state' hypothesis, as in patients with frontal lobe epilepsies. These results, together with those of studies of patients with temporal lobe epilepsy (Le Van Quyen et al., 2001) and studies in which measurements of cerebral blood ¯ow and metabolism showed changes remote from the epileptogenic focus, support the concept of a largescale `epileptic network' distributed in the brain beyond the epileptogenic area during the pre-ictal state (Weinand et al., 1997; Baumgartner et al., 1998) or the interictal state (Juhasz et al., 2000) in location-related epilepsy.
changes in cognitive activities may affect brain dynamics by modifying neuronal networks and thus favour seizure emergence. We have reported previously that non-linear methods can detect subtle modulations of continuous epileptic activity during the performance of various mental tasks (Le Van Quyen et al., 1997). Nevertheless, the phenomenon observed in Patient 10 differs from classical re¯ex epilepsy in that activation was not related to speci®c cognitive tasks (higher mental activities or visual stimulation), it did not lead reproducibly to a seizure and the neuropsychological EEG activation concerned idiopathic generalized epilepsies almost exclusively (Matsuoka et al., 2000). Changes in the EEG frequency spectrum preceded some seizures in sleeping patients. EEG rhythms slowed (by enhancement of the delta frequency and/or disappearance of fast frequencies), suggesting that the pre-ictal state in these cases might correspond to a transition into a deeper sleep stage. In fact, these pre-ictal changes may result from the conjunction of two phenomena: physiological changes in brain dynamics, responsible for a facilitating state, and secondary pathological changes of neuronal networks, underlying the emergence of interictal (Ferrillo et al., 2000) or ictal epileptiform abnormalities speci®cally at these moments in frontal lobe epilepsy (Crespel et al., 2000). Pathological changes in dynamics were detected by non-linear analysis in a variable number of channels, whereas physiological changes may concern all the channels. This conjunction may also explain the difference (not statistically signi®cant) between anticipation times determined in awake and sleeping patients. Even if pre-ictal changes in dynamics may sometimes be associated with modi®cations in classical clinical or electrophysiological parameters and thus could provide interesting correlations, in most seizures no changes were apparent either in the patient's behaviour or in classical EEG analysis, whereas the non-linear EEG analysis succeeded in detecting pre-ictal changes. Our comparison between pre-ictal changes in dynamics and ictal electroclinical data showed that anticipation times were signi®cantly correlated with the initial frequency of the seizures, but not with their duration or with seizure type, suggesting that seizure onset may obey the same dynamical laws as those governing the pre-ictal state, whereas other laws govern the unfolding and the end of the seizure.
Methodological improvements Clinical and electrical correlations with preictal changes in dynamics
Temporal correlations between changes in classical linear parameters (frequency spectrum) or in the patient's behaviour, and changes in the non-linear similarity index may also provide pathophysiological indications. In awake patients, behavioural changes were sometimes correlated temporally with changes in EEG dynamics. Analysis of seizure E of Patient 10 tends to support the facilitating state hypothesis:
653
Previous non-linear measures used in the analysis of brain activityÐthe neuronal complexity loss L* (Lehnertz and Elger, 1998) and the largest Lyapunov exponent (Iasemidis and Sackellares, 1996)Ðprovided quantitative measures of brain dynamics. The use of these indices to describe state changes in EEG dynamics implies a confrontation with brain dynamics that may be non-stationary for long periods, so that successive absolute indices are not easily comparable. Our novel method of non-linear analysis signi®cantly improves the strategy for detecting pre-ictal changes in several ways
654
V. Navarro et al.
(Le Van Quyen et al., 1999). First, we calculate a `relative index', which is more adapted to the non-stationarity of longterm electrophysiological recordings (Manuca and Savit, 1996), by measuring similarities between different dynamics that have been reconstructed by non-linear methods. Secondly, at the signal level we now consider only the phase data that contain the most relevant information about the dynamics (Sauer, 1994). This choice reduces both the time required for calculation and the sensitivity to spikes or artefacts, in contrast to the other non-linear measures, which are highly sensitive to variation in signal amplitude. The similarity index was developed primarily to detect preictal states. Before clinical application, our algorithms must be adapted for real-time operation and the sensitivity and the speci®city of this method must be evaluated in further studies. The factors that may in¯uence the speci®city of the method comprise the parameters of the algorithm and the statistical thresholds for time anticipation, but also the choice of the reference period. For interictal periods far from a seizure, the similarity index did not vary markedly while the state of vigilance of the patient was stable, in contrast with the deep and sustained deviations detected before most seizures. The rate of false detection (0.31 6 0.20 per hour) is close to that (1 false positive per hour) of the method for automatic seizure detection (and not anticipation) developed by Gotman (1990), which is currently used in most long-term video-EEG monitoring units. As previously pointed out by others (Katz et al., 1991), changes in the non-linear measure of dynamics, associated with changes in the patient's state of vigilance, currently limit the use of the method for automated long-term EEG monitoring. For daily analysis of brain dynamics, we plan to reset the reference period according to the state of vigilance of the patient. Alternative signal analysis strategies may provide us with more detailed information on the neuronal mechanisms that underlie the pre-ictal state. One promising strategy consists of measuring phase synchrony between the activities of different EEG channels, and preliminary results suggest that changes in synchronization in speci®c frequency ranges are associated with the pre-ictal state.
Conclusion
We have shown, in patients with medically refractory epilepsy originating in the neocortex, that non-linear changes in the electrical dynamics of the brain precede most seizures by a mean of 7.5 min. Moreover, our non-linear similarity method anticipates seizures as ef®ciently from the scalp as from intracranial EEG recordings. These ®ndings, together with recent results obtained from a population with temporal lobe epilepsy, suggest that pre-ictal changes in dynamics may be a general mechanism associated with seizure emergence and may be detectable in a wide population of patients with partial epilepsies originating from various sites of the brain. Furthermore, identi®cation of transitions towards pathological states argues that epilepsy belongs to the group of
`dynamical' diseases (Mackey and Glass, 1977). Dysregulation of physiological controls may underlie the transition (the pre-ictal state) of the brain system from a highly complex, healthy dynamics (the interictal state) towards a more periodic, pathological dynamics (the seizure). Finally, the ability to detect pre-ictal changes gives to the word `epilepsy' a sense opposite to its Greek root (attacked by surprise). The possibility of anticipating epileptic seizures offers multiple clinical and therapeutic prospects. Detecting preictal states would considerably facilitate the early ictal injection of single-photon CT tracer among patients in presurgical evaluation. More generally, ambulatory systems might be developed to give a warning of impending seizure to patients who are not candidates for surgery or for whom surgery has failed. Such systems could greatly decrease the risk of injury and improve quality of life by reducing patients' helplessness in the face of a disease that strikes with unpredictable timing. The time window opened before a seizure should also allow therapeutic intervention to counter the emergence of the seizure.
In Memoriam
On May 28, 2001, Francisco J. Varela died. He was the head of the Neurodynamics Group at the laboratory of Cognitive Neurosciences and Brain Imaging (CNRS UPR 640) at the PitieÂ-SalpeÃtrieÁre Hospital in Paris. During the later years of his career, he gave considerable impulsion to the group working in the ®eld of seizure anticipation by sharing his exceptional skills and considerable knowledge in neurobiology, neural dynamics, cognitive neuroscience and philosophy.
Acknowledgements
We wish to thank Richard Miles for critical reading of the manuscript and David Rudrauf for his help in statistical analysis. This work was supported by the Institut National de la Sante et de la Recherche MeÂdicale (INSERM) (V.N.), the Institut Electricite Sante de France (EDF) and Sano®Synthelabo.
References Albano AM, Muench J, Schwartz C, Mees AI, Rapp PE. Singularvalue decomposition and the Grassberger±Procaccia algorithm. Phys Rev A 1988; 38: 3017±27. Baumgartner C, Serles W, Leutmezer F, Pataraia E, Aull S, Czech T, et al. Pre-ictal SPECT in temporal lobe epilepsy: regional cerebral blood ¯ow is increased prior to electroencephalographyseizure onset. J Nucl Med 1998; 39: 978±82. Buck D, Baker GA, Jacoby A, Smith DF, Chadwick DW. Patients' experiences of injury as a result of epilepsy. Epilepsia 1997; 38: 439±44.
Anticipation of neocortical seizures Cockerell OC, Johnson AL, Sander JW, Hart YM, Goodridge DM, Shorvon SD. Mortality from epilepsy: results from a prospective population-based study. Lancet 1994; 344: 918±21. Crespel A, Coubes P, Baldy-Moulinier M. Sleep in¯uence on seizures and epilepsy effects on sleep in partial frontal and temporal lobe epilepsies. [Review]. Clin Neurophysiol 2000; 111 Suppl 2: S54±9. Devinsky O, Vickrey BG, Cramer J, Perrine K, Hermann B, Meador K, et al. Development of the quality of life in epilepsy inventory. Epilepsia 1995; 36: 1089±104. Engel J Jr. Outcome with respect to epileptic seizures. In: Engel J Jr, editor. Surgical treatment of the epilepsies. New York: Raven Press; 1987. p. 553±71.
655
in intracranial activities prior to human temporal lobe seizures. Eur J Neurosci 2000; 12: 2124±34. Le Van Quyen M, Martinerie J, Navarro V, Boon P, D'Have M, Adam C, et al. Anticipation of epileptic seizures from standard EEG recordings. Lancet 2001; 357: 183±8. Lopes da Silva FH, Pijn JP. Epilepsy: network models of generation. In: Arbib MA, editor. The handbook of brain theory and neural networks. Cambridge (MA): MIT Press; 1995. p. 367±9. Mackey MC, Glass L. Oscillation and chaos in physiological control system. Science 1977; 197: 287±9. Manuca R, Savit R. Stationarity and nonstationarity in time series analysis. Physica D 1996; 99: 134±61.
Engel J. Surgery for seizures. [Review]. New Engl J Med 1996, 334: 647±52.
Martinerie J, Adam C, Le Van Quyen M, Baulac M, CleÂmenceau S, Renault B, et al. Epileptic seizures can be anticipated by non-linear analysis. Nat Med 1998; 4: 1173±6.
Ferrillo F, Beelke M, Nobili L. Sleep EEG synchronization mechanisms and activation of interictal epileptic spikes. [Review]. Clin Neurophysiol 2000; 111 Suppl 2: S65±73.
Matsuoka H, Takahashi T, Sasaki M, Matsumoto K, Yoshida S, Numachi Y, et al. Neuropsychological EEG activation in patients with epilepsy. Brain 2000; 123: 318±30.
Goldberger AL. Non-linear dynamics for clinicians: chaos theory, fractals, and complexity at the bedside. Lancet 1996; 347: 1312±14.
Rajna P, Clemens B, Csibri E, Dobos E, Geregely A, Gottschal M, et al. Hungarian multicentre epidemiologic study of the warning and initial symptoms (prodrome, aura) of epileptic seizures. Seizure 1997; 6: 361±8.
Gotman J. Automatic seizure detection: improvements and evaluation. Electroencephalogr Clin Neurophysiol 1990; 76: 317± 324. Iasemidis LD, Sackellares JC. Chaos theory and epilepsy. Neuroscientist 1996; 2: 118±26. Juhasz C, Chugani DC, Muzik O, Watson C, Shah J, Shah A, et al. Electroclinical correlates of ¯umazenil and ¯uorodeoxyglucose PET abnormalities in lesional epilepsy. Neurology 2000; 55: 825± 35. Katz A, Marks DA, McCarthy G, Spencer SS. Does interictal spiking change prior to seizures? Electroencephalogr Clin Neurophysiol 1991; 79: 153±6. Lee SA, Spencer DD, Spencer SS. Intracranial EEG seizure-onset patterns in neocortical epilepsy. Epilepsia 2000; 41: 291±307. Lehnertz K, Elger CE. Can epileptic seizures be predicted? Evidence from nonlinear time series analysis of brain electrical activity. Phys Rev Lett 1998; 80: 5019±22. Le Van Quyen M, Adam C, Lachaux J-P, Martinerie J, Baulac M, Renault B, et al. Temporal patterns in human epileptic activity are modulated by perceptual discriminations. Neuroreport 1997; 8: 1703±10. Le Van Quyen M, Martinerie J, Baulac M, Varela F. Anticipating epileptic seizures in real time by a non-linear analysis of similarity between EEG recordings. Neuroreport 1999; 10: 2149±55. Le Van Quyen M, Adam C, Martinerie J, Baulac M, CleÂmenceau S, Varela F. Spatio-temporal characterizations of non-linear changes
Rogowski Z, Gath I, Bental E. On the prediction of epileptic seizures. Biol Cybern 1981; 42: 9±15. Sauer T. Reconstruction of dynamical systems from interspike intervals. Phys Rev Lett 1994; 72: 3811±14. Schiff SJ. Forecasting brain storms [letter]. Nat Med 1998; 4: 1117± 18. Schreiber T, Schmitz A. Classi®cation of time series data with nonlinear similarity measures. Phys Rev Lett 1997; 79: 1475±8. Semah F, Picot MC, Adam C, Broglin D, Arzimanoglou A, Bazin B, et al. Is the underlying cause of epilepsy a major prognostic factor for recurrence? Neurology 1998; 51: 1256±62. Spencer SS. Substrates of localization-related epilepsies: biologic implications of localizing ®ndings in humans. [Review]. Epilepsia 1998; 39: 114±23. Weinand ME, Carter LP, el-Saadany WF, Sioutos PJ, Labiner DM, Oommen KJ. Cerebral blood ¯ow and temporal lobe epileptogenicity. J Neurosurg 1997; 86: 226±32. Widman G, Lehnertz K, Urbach H, Elger CE. Spatial distribution of neuronal complexity loss in neocortical lesional epilepsies. Epilepsia 2000; 41: 811±17. Received February 22, 2001. Revised July 9, 2001. Second revision October 1, 2001. Accepted October 22, 2001