Document Detail

Slow modulations of high-frequency activity (40-140 Hz) discriminate preictal changes in human focal epilepsy.
Jump to Full Text
MedLine Citation:
PMID:  24686330     Owner:  NLM     Status:  In-Data-Review    
Abstract/OtherAbstract:
Recent evidence suggests that some seizures are preceded by preictal changes that start from minutes to hours before an ictal event. Nevertheless an adequate statistical evaluation in a large database of continuous multiday recordings is still missing. Here, we investigated the existence of preictal changes in long-term intracranial recordings from 53 patients with intractable partial epilepsy (in total 531 days and 558 clinical seizures). We describe a measure of brain excitability based on the slow modulation of high-frequency gamma activities (40-140 Hz) in ensembles of intracranial contacts. In prospective tests, we found that this index identified preictal changes at levels above chance in 13.2% of the patients (7/53), suggesting that results may be significant for the whole group (p < 0.05). These results provide a demonstration that preictal states can be detected prospectively from EEG data. They advance understanding of the network dynamics leading to seizure and may help develop novel seizure prediction algorithms.
Authors:
C Alvarado-Rojas; M Valderrama; A Fouad-Ahmed; H Feldwisch-Drentrup; M Ihle; C A Teixeira; F Sales; A Schulze-Bonhage; C Adam; A Dourado; S Charpier; V Navarro; M Le Van Quyen
Publication Detail:
Type:  Journal Article     Date:  2014-04-01
Journal Detail:
Title:  Scientific reports     Volume:  4     ISSN:  2045-2322     ISO Abbreviation:  Sci Rep     Publication Date:  2014  
Date Detail:
Created Date:  2014-04-01     Completed Date:  -     Revised Date:  -    
Medline Journal Info:
Nlm Unique ID:  101563288     Medline TA:  Sci Rep     Country:  England    
Other Details:
Languages:  eng     Pagination:  4545     Citation Subset:  IM    
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Descriptor/Qualifier:

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Sci Rep
Journal ID (iso-abbrev): Sci Rep
ISSN: 2045-2322
Publisher: Nature Publishing Group
Article Information
Download PDF
Copyright © 2014, Macmillan Publishers Limited. All rights reserved
open-access:
Received Day: 10 Month: 04 Year: 2013
Accepted Day: 12 Month: 03 Year: 2014
Electronic publication date: Day: 01 Month: 04 Year: 2014
collection publication date: Year: 2014
Volume: 4E-location ID: 4545
PubMed Id: 24686330
ID: 3971396
Publisher Item Identifier: srep04545
DOI: 10.1038/srep04545

Slow modulations of high-frequency activity (40–140 Hz) discriminate preictal changes in human focal epilepsy
C. Alvarado-Rojas1
M. Valderrama2
A. Fouad-Ahmed1
H. Feldwisch-Drentrup3456
M. Ihle5
C. A. Teixeira7
F. Sales8
A. Schulze-Bonhage35
C. Adam19
A. Dourado7
S. Charpier1
V. Navarro19
M. Le Van Quyena1
1Centre de Recherche de l'Institut du Cerveau et de la Moelle épinière (CRICM) INSERM UMRS 975 - CNRS UMR 7225-UPMC, Hôpital de la Pitié-Salpêtrière, Paris, FRANCE
2Department of Biomedical Engineering, Universidad de Los Andes, Bogotá, COLOMBIA
3Bernstein Center Freiburg, University of Freiburg, Freiburg, GERMANY
4Freiburg Center for Data Analysis and Modeling, University of Freiburg, GERMANY
5Epilepsy Center, University Hospital of Freiburg, GERMANY
6Department of Neurobiology and Biophysics, Faculty of Biology, University of Freiburg, Freiburg, GERMANY
7CISUC-Centro de Informatica e Sistemas da Universidade de Coimbra, Faculty of Sciences and Technology, University of Coimbra, Coimbra, PORTUGAL
8Hospitais da Universidade de Coimbra, PORTUGAL
9AP-HP, Epilepsy Unit, Groupe Hospitalier Pitié-Salpêtrière, FRANCE
aquyen@T-Online.de

In spite of available drug and surgical treatment options, a third of individuals with partial epilepsy have intractable seizures. The unpredictability of seizure occurrence limits their daily life, and underlies an enhanced risk of sudden unexpected death or morbidity1. Multiple quantitative analyses of the electroencephalogram (EEG) have aimed to define independent features and reproducible patterns that herald seizure onset2, 3, 4, 5, 6. Several studies have suggested that epileptic seizures do not occur randomly, but rather emerge from slow preictal changes in brain excitability that evolve over long timescales and predispose the brain to seizure6. In particular, a recent study - using intracranial electrodes connected to a telemetry unit implanted in ambulatory patients - showed that a prospective seizure prediction is possible with the knowledge of the patient's pre-seizure EEG patterns7. In most of the investigated patients, the warning device can successfully identify seizure occurrences in advance (from minutes to hours) and at a level better than would be expected on the basis of chance alone. However, despite encouraging results, this study has been limited by the small number of patients and complex machine learning algorithms making difficult to explore possible underlying physiological mechanisms.

Greater insights into the neuronal mechanisms of population behavior during the preictal period could not only improve seizure prediction but also help understand pathophysiological process of seizure initiation in the human brain8, 9. Recent wide-bandwidth records of local field potentials made with intracranial electrodes during presurgical evaluation have identified several new classes of electrographic activity at high-frequencies. These low-amplitude activities occur at frequencies >40 Hz, are confined to millimeter-scale tissue volumes, and thus are less visible and often confused with eye movement or muscle artifacts in standard scalp electroencephalography (EEG)10, 11. High-frequency activities, including gamma oscillations at 60–150 Hz and fast ripples at 250–500 Hz, have been associated with the initiation of epileptiform potentials and seizures in human temporal lobe and neocortical epilepsies12, 13, 14, 15. Furthermore, while high-frequency oscillations are present intermittently throughout interictal period16, 17, 18, the occurrence of high-frequency activities increases significantly from several seconds to minutes before seizure onset19, 20, 21, 22, 23. Thus, high-frequency activities may help identify periods of increased predisposition to clinical seizures4. Nevertheless, recent investigations suggested that the relationship between high-frequency activities and seizure activity several minutes before the onset is variable, and no clear trend was observed24, especially in light of the many confounding variables such as fluctuating patient state4. Epileptic and physiological high-frequency oscillations are difficult to distinguish25, 26. Both fluctuate strongly in amplitude over a time scale of several hours17 and the energy of high frequency activity changes with behavior with a maximal expression during seizure-free slow-wave sleep27, 28, 29.

The coupling between slow cortical potentials and high-frequency activities may provide novel insights into cortical network excitability30. Slow oscillations seem to be able to trigger and group local high-frequency oscillations. These cross-frequency couplings have been termed phase-amplitude coupling or “nested” oscillations and they have been suggested to represent a signature of cortical activation31, 32, 33, 34, 35. In this context, human intracranial EEG studies have identified a spatially distributed modulation of cortical high frequency oscillations in the gamma band (40–120 Hz) by theta oscillations (4–8 Hz)36, 37 and slow waves (0.5–3 Hz)29. Cross-frequency coupling of high-frequency sub-bands in the gamma range to low-frequency electrical stimulation recorded with intracranial electrodes in patients with temporal lobe epilepsy has provided important clues on seizure susceptibility38, 39. These observations led us to make a systematic study of preictal changes in cross-frequency coupling in intracranial EEG records from patients with medication-resistant partial epilepsy.

We assessed the sensitivity of cross-frequency coupling changes – defined as the ratio of correctly predicted seizures to the number of seizures investigated – and also their specificity from long-term interictal data3. Analyses were made from a large database of 53 patients with continuous multi-day intracranial EEG, collected in three European epilepsy centers (Freiburg, Germany; Paris, France; Coimbra, Portugal)40. Our analysis focused on interactions between the phase of low frequency rhythms (slow waves and theta) and the amplitude of different sub-bands of gamma rhythms. This analysis of spatial fluctuations in the coupling phase of ensembles of intracranial contacts revealed preictal changes, occurring at levels greater than chance in a small but significant number of patients. To validate our results, we compared them with predictions based on the power in individual frequency bands (delta, theta and gamma). We found that spectral powers did not lead to similar performances, suggesting that the coupling of different frequency bands is specific to the reported preictal changes.


Results

Figure 1a shows a representative subdural electrocorticographic signal recorded 2 minutes before a seizure that occurred during slow-wave sleep (Id: 13, Table 1). The raw and filtered signals in the low gamma (LG) range show the recurrence of bursts of high-frequency activity occurring preferentially at the depth-positive phase of cortical slow waves. The strongest modulation for signals from this electrode was detected between the slow wave phase and LG amplitude (Fig. 1b), supporting previous findings that slow oscillations modulate gamma power in cortical recordings28, 29, 30, 31. This cross-frequency coupling was quantified by constructing a histogram of LG amplitude over the phase of slow activity, showing a unimodal distribution41 (Fig. 1c). The single peak confirmed that gamma oscillations occurred preferentially at a mean coupling phase (ϕc) close to 0 rad with respect to the intracranial slow waves. We tracked temporal changes in the coupling phase over ensembles of intracranial contacts by extracting ϕc for each contact in consecutive non-overlapped 1-minute windows. Figure 1d illustrates a map of the temporal evolution of the number of contacts with a specific phase bin during an interictal state recorded 5 hours before a seizure. It shows that most contacts had a similar phase of coupling around 0.5 rad, which remained stable over several hours. However, tens of minutes before the seizure, large and consistent preictal deviations towards 0 rad were observed at multiple locations. Conventional power spectral analyses were less sensitive to these preictal changes (Fig. 1d). Visual inspection of the signals confirmed that high-frequency activities were strongly modulated, during the preictal state, by slow waves in many channels (Fig. 1e). Projected on post implantation MRI scans, contacts with these specific preictal phase fluctuations were distributed over spatially broad regions both within and outside the focal region (Fig. 1f).

We investigated long-term coupling between fast and slow oscillations (see Methods for a detailed description of the strategy). First, the proportion of intracranial contacts at a specific mean phase interval was estimated over multi-day recordings. Second, an alarm was raised when this proportion exceeded a critical threshold. The threshold was adjusted and parameters including preictal durations and frequency bands were optimized on a first training part of the data obtained from each patient. For the selected optimal parameters, the last part of data was used to test prospectively discrimination between preictal from non-preictal periods. Performance results are reported for this test period, with parameters optimized from the training period. We tested the statistical significance of any changes, by comparing performance with a random predictor, which does not exploit any information contained in the EEG data42. Figure 2a–b shows a long-term analysis of data from the same patient as in Figure 1. Using a preictal duration of 60 minutes, optimal for the training period, prospective performance was above chance level. Our analysis gave values of sensitivity (SS) of 73% and false prediction rate (FPR) of 0.29/h, suggesting that most seizures were preceded by a specific preictal state and the false prediction rate was low. In this patient, seizures occurred predominantly during the night (10 pm–7 am), suggesting that sleep modulated seizure susceptibility. Nevertheless, preictal changes were present before daytime seizures (see sz13 in Fig. 2b) and such changes were not always evident before night seizures (see for example the testing phase between 18 and 24 hours in Fig. 2a). Moreover, the preictal phase distribution shifted to another phase configuration (from 0.5 to 0 rad), absent from seizure-free night periods, suggesting that this phenomenon did not reflect a simple sleep stage (Fig. 2c). To further confirm that our results did not mainly reflect power fluctuations related to different states of vigilance, we compared them with predictions based on the relative power in individual frequency bands (delta: 0.1–4 Hz, theta: 4–8 Hz, gamma: 30–140 Hz). Fig. 2a (Bottom) shows the alarms triggered by applying a threshold to the spectral powers averaged over all channels and following the same strategy as that using cross-frequency coupling. For this patient, power-based predictions were less sensitive and specific to preictal changes. In particular, it can be observed that delta activity was continuously dominant in the EEG and only slightly modulated by sleep (Fig. 2a).

Table 1 shows the parameters selected for the whole group of 53 patients, as well as the performance (SS and FPR) of the algorithm. Performance varied strongly between patients, presumably due to diverse etiologies, exogenous triggers, and the inherent physiological variability. Even so, compared to a random predictor, significant preictal changes were identified prospectively at above chance level in 13.2% of the patients (7/53). Based on a binomial approximation (see Methods for details), this result is significant for the whole group of patients at the 5% level (Pbinom{7,53,0.05} = 1.6% which is less than 13.2%). For the 7 patients with statistical significant results, that average sensitivity was 68% (range: 36–100%) and the FPR was 0.33/h (range: 0.08–0.72/h). There was a preference for slow wave/LG coupling in 5/7 patients where preictal phase identification was successful (Table 1). All significant cases corresponded to epilepsies with temporal or frontal foci. A third of cases (33%) were significant for patients with frontal lobe epilepsy, and only 10% were significant for cases of temporal lobe epilepsy. Significant performances during the test period, were comparable for patients with temporal lobe (SS = 66% and FPR = 0.33/h) and frontal lobe epilepsy (SS = 62% and FPR = 0.27/h). For all data sets from the different patients, we varied the duration of preictal window lengths on the training period and found an optimal duration of 60 minutes in 6/7 patients with significant performances. Finally, we found no substantial differences in prediction of seizures occurring during the night and during the day for the significant patients: 44% of seizures occurring during the day and 56% of night-time seizures were predicted). Similar performances were observed for nocturnal or day-time seizures (10 pm-7 am; SS = 73% and FPR = 0.25/h during the night; SS = 64% and FPR = 0.19/h during the day). To confirm the specificity of our results, we compared them with predictions based on the power in individual frequency bands. Only 4 patients had significant results for delta (SS = 64% and FPR = 0.15/h), 2 for theta (SS = 50% and FPR = 0.07/h) and 1 for gamma band (SS = 62% and FPR = 0.3/h). These results suggest that the performances of our detector are relatively independent from the different states of vigilance. In addition, the coupling between these different rhythms carries more information about preictal changes than the information of individual frequency bands taken separately.

Finally, we made a detailed analysis of the coupling phase distribution during preictal/interictal states in the 7 patients with statistical significant results. Two distinct types of phenomena emerged. In most cases (5/7), the preictal and interictal phase distributions were very similar except that wider brain regions were involved during preictal periods. Figure 3a shows a representative case (Id: 20 in Table 1). In this patient with mesial temporal epilepsy, seizures occurred preferentially at night, and analysis of slow wave/LG phase coupling revealed a low FPR (0.1/h) and a high sensitivity (56%). The preictal phase distribution was centered around −π, close to the one detected in day or night-time interictal states but with a larger involvement of intracranial contacts (Fig. 3a, right). In contrast, the preictal phase distribution shifted towards a different preferred phase in a minority of cases (2/7). Figure 3b shows a representative case (Id: 46 in Table 1). In this patient with frontal lobe epilepsy, the delta/HG phase coupling was selected to give the best performances with a remarkable high specificity (SS = 36% and FPR = 0.08/h). As for the patient of Fig. 12, pre- and interictal phase dynamics had distinct modulation properties and preictal shifts towards -0.9 rad were observed at multiple locations. In all significant cases, preictal changes occurred in regions close to but outside the seizure onset zone, indicating a widespread and diffuse spatial distribution (Fig. 3c, d).


Discussion

Our study evaluated the feasibility of long-term seizure forecasting using intracranial EEG in a multicenter group of epilepsy patients with continuous long-term recordings (total analyzed duration of 531 days with total 558 clinical seizures). Bias introduced by selecting epochs for analysis was removed by using continuous multi-day EEG records from each patient3. Furthermore, as suggested3, 43, we used a first part of available data for training and adjusting parameters for each patient and a second part to test prediction quasi-prospectively on unselected data. Finally, since our data was obtained with different EEG acquisition systems in three different hospitals, the evaluation provided a realistic sample of the heterogeneity of clinical circumstances and focal epileptic syndromes44. We found that a measure of brain excitability based on the coupling between low-frequency phase and high-frequency amplitude was able to identify preictal states for a significant number of patients (13.2%). Therefore, our observations demonstrated that a quasi-prospective analysis can distinguish between preictal and non-preictal states on long-term intracranial EEG recordings.

The existence of a preictal state in partial epilepsy is debated6. Despite several encouraging results, only very few studies have shown that a (quasi-) prospective prediction algorithm can perform above chance level3, 4. For example, a forecasting algorithm using multiple spectral power bands as features and support-vector machine classification was tested on intracranial recordings of 18 patients (433 hours; 80 seizures) and reported a high sensitivity (97.5%) and low false prediction rate (0.27 per hour)43. In particular, the power changes in gamma bands (including the low and high gamma: 30–128 Hz) have been shown to be very relevant for prediction. Similarly, another study reported sensitivities of 0–100% and false prediction rates of 0–1.67/h in 6 patients based on the energy and entropy of high frequency activity (50–450 Hz)45. Nevertheless, the EEG recordings in these two studies were not continuous, with gaps between interictal and preictal segments, and data containing artifacts were sometimes removed. These procedures weaken long-term evaluation of the prediction performance during changing physiological or epileptic states. Therefore, although these studies may present higher sensitivity and/or specificity than ours, these results are not directly comparable. Recently, using a chronically implanted device running an algorithm which was trained on the basis of each patient's EEG signals, Cook et al.7 demonstrated that prospective seizure prediction is possible on long-term recordings of several months. The algorithms used in this work relied purely on classification of power measures in different frequency bands (8–128 Hz) and derived from subdurally recorded neural activity. During the prospective evaluation, the warning device worked to a level better than would be expected by chance in 8 over 11 patients. Using the same implanted device, but in dogs with naturally occurring epilepsy, a similar seizure forecasting algorithm performed significantly better than a random predictor in all investigated cases46. In our study, we used data from a large group of patients obtained during clinical monitoring for epilepsy surgery. Based on this database, we were only able to find a small group of patients having seizures that can be forecasted better than chance. Our observation seems to contrast with Cook's optimistic observations using long-term outpatient “real life” recordings and reporting high prediction performances in most of the investigated patients. One reason of this discrepancy could be the known limitations related to use of inpatient intracranial EEG data to evaluate prediction algorithms, including relatively short duration of EEG recordings (in particular for training purposes) and the effects of drug changes or surgery. Also, differences in statistical thresholds could explain different performances. Nevertheless, both studies support the feasibility of long-term seizure forecasting.

Statistically significant results were obtained in a small group of patients, which does not facilitate a full identification of pathophysiological mechanisms. Nevertheless, we resolved two types of preictal phenomena. In the first type (see data from patients of Fig. 2 or Fig. 3b), the preictal phase distribution shifted toward a distinct preferred phase. This situation may reflect the existence of a specific pathological state distinct from physiological states and independent of changes in the state of vigilance. Alterations in the modulation of high-frequency power may provide an indirect access to pathological spiking activities47 and reflect abnormal interactions between coupled networks48. In the second type of preictal phenomena (such as that shown in Fig. 3a), their phase distributions were very similar than those during interictal states but involved more distant brain regions. This may reflect a widespread increase in brain excitability, thereby a global state that is more susceptible to seizures. Clinical evidence suggests that certain normal states, particularly related to sleep and arousal, can favor seizure occurrence in epilepsy patients. For example, most nocturnal partial seizures occur during slow wave sleep while few ictal events occur during rapid eye movement (REM) sleep49, 50. These preictal states may not reflect a deterministic, transitional state that inevitably leads to seizure, but a permissive “pro-ictal” state conforming to a normal brain state with a higher probability of seizure generation51. Thus, the high-level of false-positives detected by our algorithm may reflect underlying physiological mechanisms rather than short-comings of our detection that tend to label distinct events as similar. Nevertheless, as shown by our comparison with the power fluctuations in several frequency bands, there is not one to one correspondence between preictal changes and specific states of vigilance; this complex relationship is an area of ongoing research52.

We have shown that cross-frequency coupling can provide new insights into the transition from interictal to ictal states, but clinical applicability will depend on several other factors. False seizure prediction rates greater than about 0.15 FP/h are generally felt to be unacceptable for clinical application53, 54. Our most favorable analyses provided an average sensitivity of 68% and a specificity of 0.33 FP/h, insufficient for clinical application. In only one of seven patients with significant performance, we observed both a remarkably high specificity and sensitivity (Fig. 3a). Nevertheless, we are encouraged by these results and believe that the discrimination framework we present here may be improved in the future towards a clinically acceptable performance if combined with other seizure prediction techniques55.


Methods
Database

Long-term intracranial EEG recordings from 53 epilepsy patients (26 males; age range, 3–63 years; mean age: 32 years) suffering from medically intractable partial epilepsy were analyzed. The data were recorded in three different epilepsy units: Epilepsy Center, University Medical Centre of Freiburg, Germany; Unité d'Épilepsie of the Pitié-Salpêtrière Hospital, Paris, France; Hospitais da Universidade de Coimbra, Portugal. Data acquisition and database were performed in accordance with the approved guidelines of the local ethic committees (Ethik-Kommission der Albert-Ludwigs-Universität Freiburg; Comité consultatif sur le traitement de l'information en matière de recherche dans le domaine de la santé, Pitié-Salpêtrière University Hospital; and Ethics Committee of the Coimbra University Hospital, respectively). Patients gave a written consent for a research use of these data. Funded by the European Union, this database was established to compile high-quality, long-term continuous EEG data, enriched with clinical metadata40 (http://www.epilepsiae.eu/). A standardized EEG annotation protocol was developed to ensure the comparability and reliability of seizure onset times at all project sites56. A total of 531 days of EEG records (~10 days/patient) included 558 seizures (~11 seizures/patient). Records include at least 5 days of continuous EEG from each patient with at least 5 clinically manifest epileptic seizures separated by more than 1.5 hours. EEG data were recorded using different digital video EEG systems with different sampling rates in the three different centres; Nicolet, 256 Hz; Micromed, 400 Hz; Compumedics, 512 Hz and Neurofile NT, 1024 Hz. Subclinical intracranial electrographic events were not analyzed. Table 1 shows patient characteristics, recordings details and locations of the seizure onset zone.

Measuring the mean coupling phases

Cross-frequency coupling was defined by a statistical measure of interaction between different frequency bands. Specifically the phase of a slow oscillation was compared with the amplitude of a fast oscillation. The highest amplitude occurred at the so-called mean coupling phase. In accordance with previous human studies29, 36, 37, we investigated the coupling between delta (0.5–3 Hz) or theta (3–8 Hz) for low frequency phases, and low gamma (LG: 40–70 Hz) or high gamma (HG: 70–140 Hz or 70–120 Hz when sampling frequency was 256 Hz) for high frequency amplitude envelopes. We extracted coupling by implementing an algorithm previously proposed41, where raw signals were first filtered in the bands of interest, and then phase and amplitude were extracted using the Hilbert transform. Raw signals were filtered for frequency bands of interest using an 8-order Butterworth forward-backward IIR filter to avoid phase delays and improve frequency selectivity. Instantaneous phase ϕ(t) and instantaneous amplitude envelope A(t) were then extracted for delta (ϕδ) and theta (ϕθ), and for LG (ALG) and HG (AHG), respectively. The Hilbert Transform produces an analytical representation of the signal, with real part x(t) corresponding to the filtered signal and complex part y(t), its phase-delayed version by π/2:

[Formula ID: m1]
In this complex representation, the angle and magnitude of the transform correspond respectively to the instantaneous phase and envelope of the filtered signal (see Fig. 1a for an example). The modulation between four possible low phase – high amplitude pairs (ϕx(t) Ax(t)), was computed as follows. First, the range of the phase signal ϕx(t), with values within the interval (−π,π), was divided into bins ϕiwith i = 1,…,40. Then, the time indexes kiwhen the phase fell in the interval ϕi ≤ ϕx(ki) < ϕi+1, were determined. For each bin, amplitudes at time ki of the high frequency time series were averaged (<Ax(ki)>). The distribution of mean amplitudes versus phase defined the phase of the slower activity, where high frequencies tend to appear (Fig. 1b). We quantified this tendency by approximating the distribution to a Von Mises function (circular Gaussian), and extracting the mean phase representing the preferred phase of coupling ϕc. A statistical threshold (based on a phase randomization procedure for instance) was not used41. We thus maintained computationally efficiency in treating data from data epochs of 1 minute. Codes are freely available in the EPILAB platform55 at: http://www.epilepsiae.eu/project_outputs/epilab_software/.

Identifying fluctuations of mean coupling phase patterns

We followed the time course of mean coupling phases of multiple contacts, during consecutive non-overlapping time windows of 1 minute, by determining the proportion of contacts with a mean coupling phase at a specific phase interval [ϕc1,ϕc2]. Figure 1c shows an analysis for a representative patient during interictal and preictal periods. Long-term trends of coupling phase were determined by applying a first-order Kalman filter to the proportion of electrodes within [ϕc1,ϕc2]. Thus short-time fluctuations were smoothed, as previously proposed43. A variable threshold was then set to produce an alarm when the number of intracranial contacts exceeded a defined value (Fig. 2a, b). After an alarm, a new alarm could be generated only after a refractory period corresponding to the defined preictal time.

Prospective evaluation of preictal changes

Application of highly optimized methods to small and preselected data may result in overestimating the results that are poorly reproducible on unselected, longer continuous data3. While algorithm parameters may be adjusted to data from an individual patient, this optimization must be performed only on one part of the data (the training set) excluded from the data set used subsequently to assess algorithm performance (the testing set). For training, we used continuous data including the first 4 seizures and at least 10 hours of recording. This training data permitted parameter selection to optimize sensitivity and specificity. We determined optimal values for the coupling band (among 4 pairs: slow waves or theta versus LG and HG), phase interval [ϕc1,ϕc2] (among 11 intervals) and preictal duration (10, 30, 60 min). Statistical performance was calculated, as usual, in terms of sensitivity (SS, ratio between correct detections and total seizures), and false prediction rate (FPR, number of false predictions per hour). Optimization was performed such that the performance was the closest to the optimal value (SS of 100%, FPR of 0), with distance measured by . Finally, the test set containing unknown data and seizures was used to derive a prospective statistical assessment of performance of the optimized algorithm. FPR was defined as the number of false alarms during the test set divided by the duration in which false alarms could be triggered (with refractory periods excluded), which is obtained by subtracting the time under false warning from the total interictal period3. Mathematically, the FPR is expressed as:

[Formula ID: m2]

Comparison with a random predictor

Once the sensitivity and specificity of a prospective prediction algorithm has been assessed on the test dataset, it is necessary to ask whether it is indeed superior to the chance level. We evaluated statistical significance by comparing results to chance established by a random predictor, which produced alarms in the absence of information from the EEG signal42. With random alarm generation, the probability to predict at least k out of K seizures was given by a binomial distribution:

[Formula ID: m3]
Where P indicated the probability (Poisson process) of that a single alarm was triggered during the preictal time. A unique predictor was used for the final evaluation of the testing phase of the prospective approach. Therefore the critical level of sensibility for the random predictor σrand was stated as follows:
[Formula ID: m4]
Where α represented the significance level of the predictor (α = 0.01 in this case). For the same value of FPR, the sensitivity of our method must be above σrand to be superior to the random predictor, and thus reach statistical significance.

Statistical evaluation on the whole group of patient

In order to test whether the observed results can be considered significant for a whole group of patients, the number of patients with statistical significant results can be calculated for the null hypothesis of no true predictive performance. For a significance level of α, the probability to observe, for at least n of N patients, sensitivities larger than σrand follows a binomial distribution:

[Formula ID: m5]

Comparison to prediction based on the power of delta, theta or gamma bands

To evaluate whether our observations based on cross-frequency coupling did not simply reflect power fluctuations related to different states of vigilance, a similar strategy was developed to generate alarms based on the relative power of different bands: delta (0.1–4 Hz), theta (4–8 Hz) and gamma (30–140 Hz). After filtering the raw signals, the power spectral density was computed by Burg method over consecutive windows of 5 seconds. Then, the power of each band relative to the power in the whole frequency range was extracted. The averaged relative power over all channels was used to determine an optimal threshold above which an alarm of preictal period was triggered. As in our cross-frequency coupling strategy, the predictor based on the individual bands was first optimized on the training data (in sample definition of preictal times and thresholds). Finally, the algorithm performance (sensitivity and FPR) and the comparison with a random prediction were assessed on independent data (out-of-sample).


Author Contributions

M.L.V.Q. designed the experiment; C.A.R., V.N. and M.L.V.Q. performed research; C.A.R., M.V., A.F. and H.F.D. analyzed data; C.A.R. and M.L.V.Q. wrote the manuscript; M.I., C.A.T., F.S., A.S.B., C.A., A.D. and S.C. contributed data/analytic tools. All authors reviewed the manuscript.


This work was supported by the European Union-FP7 Project EPILEPSIAE (Evolving Platform for Improving Living Expectation of Subjects Suffering from IctAl Events, Grant No 211713) and by the program Investissements d'avenir ANR-10-IAIHU-06, France. C.A.R thanks the Administrative Department for Science, Technology and Innovation (COLCIENCIAS), Colombia, for financial support. V.N. was supported by Contrat d'Interface INSERM. We also thank Richard Miles for critical reading of the manuscript.


References
Cockerell O. C.,et al.Mortality from epilepsy: results from a prospective population-based study. Lancet344, 918–921 (Year: 1994).7934347
Jachan M.,, Feldwisch-Drentrup H.,, Schelter B., & Timmer J.. The history of seizure prediction. In: Schelter B., Timmer J., & Schulze-Bonhage A. Seizure prediction in epilepsy: from basic mechanisms to clinical applications (Wiley, New York, Year: 2008).
Mormann F.,, Andrzejak R. G.,, Elger C. E., & Lehnertz K.,Seizure prediction: the long and winding road. Brain130, 314–333 (Year: 2007).17008335
Stacey W.,, Le Van Quyen M.,, Mormann F., & Schulze-Bonhage A.,What is the present-day EEG evidence for a preictal state?Epilepsy Res.97, 243–251 (Year: 2011).21885253
Zaveri H. P.,, Frei M. G.,, Arthurs S., & Osorio I.,Seizure prediction: The Fourth International Workshop. Epilepsy Behav.19, 1–3 (Year: 2010).20674508
Richardson M., & Jefferys J.,Introduction – Epilepsy Research UK Workshop 2010 on “Preictal Phenomena”. Epilepsy Res.97, 229–230 (Year: 2011).22094447
Cook M. J.,et al.Prediction of seizure likelihood with a long-term, implanted seizure advisory system in patients with drug-resistant epilepsy: a first-in-man study. Lancet Neurol.12, 563–71 (Year: 2013).23642342
Truccolo W.,et al.Single-neuron dynamics in human focal epilepsy. Nat. Neurosci.14, 635–41 (Year: 2011).21441925
Stead M.,et al.Microseizures and the spatiotemporal scales of human partial epilepsy. Brain133, 2789–97 (Year: 2010).20685804
Le Van Quyen M.,Editorial “Special issue on High-frequency oscillations in cognition and epilepsy”. Prog. Neurobiol.98, 239–318 (Year: 2012).22749711
Worrell G.,High-Frequency Oscillations Recorded on Scalp EEG. Epilepsy Curr.12(2), 57–58 (Year: 2012).22473542
Fisher R. S.,, Webber W. R.,, Leeser R. P.,, Arroyo S., & Uematsu S.,High frequency EEG activity at the start of seizures. J. Clin. Neurophysiol.9, 441–448 (Year: 1992).1517412
Alarcon G.,, Binnie C. D.,, Elwes R. D., & Polkey C. E.,Power spectrum and intracranial EEG patterns at seizure onset in partial epilepsy. Electroencephalogr. Clin. Neurophysiol.94, 326–337 (Year: 1995).7774519
Traub R. D.,et al.A possible role for Gap junctions in generation of very fast EEG oscillations preceding the onset of, and perhaps initiating, seizures. Epilepsia42, 153–170 (Year: 2001).11240585
Jirsch J. D.,et al.High-frequency oscillations during human focal seizures. Brain129, 1593–1608 (Year: 2006).16632553
Bragin A.,, Engel J.,, Wilson C. L.,, Fried I., & Buzsaki G.,High-frequency oscillations in human brain. Hippocampus9, 137–142 (Year: 1999).10226774
Crépon B.,et al.Mapping interictal oscillations greater than 200 Hz recorded with intracranial macroelectrodes in human epilepsy. Brain133, 33–45 (Year: 2010).19920064
Blanco J. A.,et al.Data mining neocortical high-frequency oscillations in epilepsy and controls. Brain134, 2948–959 (Year: 2011).21903727
Worrell G. A.,et al.High-frequency oscillations and seizure generation in neocortical epilepsy. Brain127, 1496–1506 (Year: 2004).15155522
Khosravani H.,et al.Spatial localization and time-dependant changes of electrographic high frequency oscillations in human temporal lobe epilepsy. Epilepsia50, 605–616 (Year: 2009).18717704
Jiruska P.,et al.High-frequency network activity, global increase in neuronal activity and synchrony expansion precede epileptic seizures in vitro. J. Neurosci.30, 5690–5701 (Year: 2010).20410121
Zijlmans M.,et al.Ictal and interictal high frequency oscillations in patients with focal epilepsy. Clin Neurophysiol.122, 664–671 (Year: 2011).21030302
Grenier F.,, Timofeev I., & Steriade M.,Neocortical Very Fast Oscillations (Ripples, 80–200 Hz) During Seizures: Intracellular Correlates. J. Neurophysiol.89, 841–852 (Year: 2003).12574462
Jacobs J.,et al.High frequency oscillations (80–500 Hz) in the preictal period in patients with focal seizures. Epilepsia50, 1780–1792 (Year: 2009).19400871
Engel J.,, Bragin A.,, Staba R., & Mody I.,High-frequency oscillations: what is normal and what is not?Epilepsia50, 598–604 (Year: 2009).19055491
Menendez de la Prida L., & Trevelyan A.,Cellular mechanisms of high frequency oscillations in epilepsy: on the diverse sources of pathological actitivies. Epilepsy Res.97, 308–317 (Year: 2011).21482073
Bagshaw A. P.,, Jacobs J.,, LeVan P.,, Dubeau F., & Gotman J.,Effect of sleep stage on interictal high frequency oscillations recorded from depth macroelectrodes in patients with focal epilepsy. Epilepsia50, 617–628 (Year: 2009).18801037
Le Van Quyen M.,et al.Large-scale microelectrode recordings of high frequency gamma oscillations in human cortex during sleep. J Neurosci.30, 7770–7782 (Year: 2010).20534826
Valderrama M.,et al.Cortical mapping of gamma oscillations during human slow wave sleep. PLoS One7, e33477 (Year: 2012).22496749
Steriade M.,Neuronal Substrates of Sleep and Epilepsy (Cambridge University Press, Cambridge, Year: 2003).
Steriade M.,, Amzica F., & Contreras D.,Synchronization of fast (30–40 Hz) spontaneous cortical rhythms during brain activation. J. Neurosci16, 392–417 (Year: 1996).8613806
Lakatos P.,et al.An oscillatory hierarchy controlling neuronal excitability and stimulus processing in the auditory cortex. J. Neurophysiol.94, 1904–1911 (Year: 2005).15901760
Canolty R. T., & Knight R. T.,The functional role of cross-frequency coupling. Trends Cogn. Sci.14, 506–514 (Year: 2010).20932795
Fell J., & Axmacher N.,The role of phase synchronization in memory processes. Nat. Rev. Neurosci.12, 105–118 (Year: 2011).21248789
Belluscio M. A.,, Mizuseki K.,, Schmidt R.,, Kempter R., & Buzsaki G.,Cross-frequency phase-phase coupling between theta and gamma oscillations in the hippocampus. J. Neurosci.32, 422–435 (Year: 2012).
Canolty R. T.,et al.High gamma power is phase-locked to theta oscillations in human neocortex. Science313, 1626–1628 (Year: 2006).16973878
van der Meij R.,, Kahana M., & Maris E.,Phase-amplitude coupling in human electrocorticography is spatially distributed and phase diverse. J. Neurosci.32, 111–123 (Year: 2012).22219274
Kalitzin S.,, Velis D.,, Suffczynski P.,, Parra J., & Lopes da Silva F.,Electrical brain-stimulation paradigm for estimating the seizure onset site and the time to ictal transition in temporal lobe epilepsy. Clin. Neurophysiol.116, 718–28 (Year: 2005).15721087
Freestone D. R.,et al.Probing for cortical excitability. Conf. Proc. IEEE Eng. Med. Biol. 20112011, 1644–1647 (Year: 2011).
Klatt J.,et al.The EPILEPSIAE database: An extensive electroencephalography database of epilepsy patients. Epilepsia53, 1669–1676 (Year: 2012).22738131
Tort M.,, Komorowski R.,, Eichenbaum H., & Kopell N.,Measuring phase-amplitude coupling between neuronal oscillations of different frequencies. J. Neurophysiol104, 1195–1210 (Year: 2010).20463205
Schelter B.,et al.Testing statistical significance of multivariate time series methods for epileptic seizure prediction. Chaos16, 013108 (Year: 2006).16599739
Park Y.,, Luo L.,, Parhi K. K., & Netoff T.,Seizure prediction with spectral power of EEG using cost-sensitive support vector machines. Epilepsia52, 1761–70 (Year: 2011).21692794
Lehnertz K., & Litt B.,The First International Collaborative Workshop on Seizure Prediction: summary and data description. Clin. Neurophysiol.116, 493–505 (Year: 2005).15721063
Gadhoumi K.,, Lina J. M., & Gotman J.,Discriminating preictal and interictal states in patients with temporal lobe epilepsy using wavelet analysis of intracerebral EEG. Clin. Neurophysiol.123, 1906–1916 (Year: 2012).22480601
Howbert J. J.,et al.Forecasting Seizures in Dogs with Naturally Occurring Epilepsy. PLoS ONE9, e81920 (Year: 2014).24416133
Buzsáki G.,, Anastassiou C. A., & Koch C.,The origin of extracellular fields and currents - EEG, ECoG, LFP and spikes. Nat. Rev. Neurosci.13, 407–420 (Year: 2012).22595786
Froriep U. P.,et al.Altered theta coupling between medial entorhinal cortex and dentate gyrus in temporal lobe epilepsy. Epilepsia11, 1937–47 (Year: 2012).22984867
Bazil C. W., & Walczak T. S.,Effects of sleep and sleep stage on epileptic and nonepileptic seizures. Epilepsia38, 56–62 (Year: 1997).9024184
Crespel A.,, Baldy-Moulinier M., & Coubes P.,The relationship between sleep and epilepsy in frontal and temporal lobe epilepsies: practical and physiopathologic considerations. Epilepsia39, 150–157 (Year: 1998).9577994
Wong S.,, Gardner A. B.,, Krieger A. M., & Litt B.,A stochastic framework for evaluating seizure prediction algorithms using hidden Markov models. J. Neurophysiol.97, 2525–2532 (Year: 2007).17021032
Sedigh-Sarvestani M.,et al.Rapid eye movement sleep and hippocampal theta oscillations precede seizure onset in the tetanus toxin model of temporal lobe epilepsy. J. Neurosci.34, 1105–14 (Year: 2014).24453303
Aschenbrenner-Scheibe R.,et al.How well can epileptic seizures be predicted? An evaluation of a nonlinear method. Brain126, 2616–2626 (Year: 2003).14506067
Schulze-Bonhage A.,et al.Views of patients with epilepsy on seizure prediction devices. Epilepsy Behav.18, 388–396 (Year: 2010).20624689
Teixeira C.,et al.EPILAB: A software package for studies on the prediction of epileptic seizures. J. Neurosci. Methods200, 257–271 (Year: 2011).21763347
Ihle M.,et al.EPILEPSIAE – A European epilepsy database. Comput. Methods Programs Biomed.106, 127–138 (Year: 2011).20863589

Article Categories:
  • Article


Previous Document:  Soft-solid deformation mechanics at the tip of an embedded needle.
Next Document:  Clinical Application of the "Pedicle Exposure Technique" for Atlantoaxial Instability Patients With ...