Phase-Amplitude Coupling between Infraslow and High-Frequency Activities is a Potential Biomarker for Seizure Prediction
========================================================================================================================

* Hiroaki Hashimoto
* Hui Ming Khoo
* Takufumi Yanagisawa
* Naoki Tani
* Satoru Oshino
* Haruhiko Kishima
* Masayuki Hirata

## Abstract

**IMPORTANCE** This research describes a method to accurately predict the onset of epileptic seizures; this will help treat patients timely, prevent future seizures, and improve outcomes.

**OBJECTIVE** We aimed to assess whether the phase-amplitude coupling (PAC) between infraslow activities (ISA) and high-frequency activities (HFA) increases before seizure onset.

**DESIGN AND SETTING** This retrospective, single-center case series included patients admitted to the neurosurgery department at Osaka University Hospital in Suita, Osaka, from July 2018 to July 2019.

**PARTICIPANTS** We enrolled seven patients with drug-resistant focal epilepsy who underwent intracranial electrode placement as part of a presurgical invasive electroencephalography study.

**MAIN OUTCOMES AND MEASURES** We comparatively analyzed the ISA, HFA, and ISA-HFA PAC in the seizure onset zone (SOZ) or non-SOZ (nSOZ) in the interictal, preictal, and ictal states.

**RESULTS** We recorded 15 seizures in seven patients [1 female (14%); mean (SD) age = 26 (12) years; age range, 15-47 years]. HFA and ISA were larger in the ictal states than in the interictal and preictal states. During seizures, the HFA and ISA of the SOZ were larger and earlier than those of nSOZ. In the preictal states, the ISA-HFA PAC was larger than that of the interictal states, and it began increasing at 93 seconds before the seizure onset (95% confidence interval: −116 – −71 s). There were no differences in the values and time of ISA-HFA PAC between both zones. Our phase-based analysis revealed differences between the SOZ- and nSOZ-PAC. In SOZ, the HFA amplitudes were tuned at the trough of the ISA oscillations, and in nSOZ, the HFA amplitudes were tuned at the peak of these oscillations. The receiver-operating characteristic curve showed that the ISA-HFA PAC of the SOZ showed the highest discrimination performance in the preictal and interictal states, with an area under the curve (AUC) of 0.926. However, ISA-HFA PAC was not suitable to differentiate between SOZ and nSOZ (interictal AUC = 0.555, preictal AUC = 0.691, and ictal AUC = 0.646).

**CONCLUSION AND RELEVANCE** This study demonstrated the novel insight that ISA-HFA PAC increases before the onset of seizures, regardless of the seizure onset zone. Our findings indicate that ISA-HFA PAC is a potential biomarker for predicting the onset of seizures and may be valuable to physicians who routinely treat epileptic patients.

**Question** Is phase-amplitude coupling (PAC) between infraslow activities (ISA) and high-frequency activities (HFA) a useful biomarker for seizure prediction?

**Findings** In this case series study on 15 focal-onset seizures in seven epileptic patients who underwent intracranial electrode placement, we found that a PAC of the ISA phase and HFA amplitude achieved significantly higher values in preictal states than in the interictal states, and ISA-HFA PAC of the seizure onset zone (SOZ) began increasing at 93 seconds before seizure onset (SO), while both HFA and ISA increased after SO. The receiver-operating characteristic curve showed that the ISA-HFA PAC of the SOZ showed the highest discrimination performance in the preictal and interictal states, with an area under the curve of 0.926.

**Meaning** This study demonstrates that ISA-HFA PAC can differentiate between the preictal and interictal states of a seizure, indicating that it is a potential marker for seizure prediction.

## Introduction

Epileptic seizures are common neurological disorders. Timely detection of seizures is important for physicians to diagnose and quantitatively measure epilepsy. It was recently reported that seizures can be detected by applying deep learning methods1.

Infraslow activities (ISAs) and high-frequency activities (HFA) are key biomarkers of seizure detection2-7. They are both included in wideband electroencephalogram measurements and are measured by intracranial electroencephalogram (iEEG). ISAs, which are also referred to as direct current (DC) shifts8, are low-frequency components. Clinically, modern AC amplifiers with a time constant of 10 seconds allow filter-setting to record ISA recording (cutoff frequency of a high-pass filter is above 0.016 Hz)5. HFAs can be physiological9 or pathological10; epileptic HFAs are usually > 80 Hz4 and are pathological-type. High-frequency oscillations (HFOs), which are subgroups of HFAs7, are isolated oscillations that stand out from the background and can be further classified into ripples (80-250 Hz) and fast ripples (250-500 Hz)4,11,12.

The relationship between HFA and ISA has not yet been established. Ictal HFA and ISA are likely to occur in the same contacts2,5. Both occur earlier than conventional iEEG changes, but ictal ISA is observed more often and may occur earlier than HFA3,5,13. The ictal ISAs that precede ictal HFAs are called active DC shifts, whereas ictal ISAs that occur immediately after the onset of conventional ictal EEG patterns and HFAs are called passive DC shifts5. In interictal states, an ISA accompanied by an HFA could be a useful surrogate marker of the epileptogenic zone––this is called a “red slow”5,14.

Phase-amplitude coupling (PAC) is used to investigate the relationship between the low- and high-frequency bands on the EEG15. PAC analysis measures the degree of synchronization between the phases with low-frequency oscillation and high-frequency amplitude16. Ictal HFA amplitudes coupled with δ 17,18, θ19,20, and α19 phases, and β-HFA coupling are useful markers for seizure detection21. These pathological PACs are usually accompanied by ictal HFA, whereas physiological PACs that are involved in motor-related HFAs (high γ band) appear before the HFA increases22,23. Physiological PAC may induce a delay in physiological HFA; however, it is not known whether a preceding seizure-related PAC can induce a delay in ictal HFA or not. Moreover, a coupling between a frequency of 0.1 Hz and HFA has been previously reported17, but it is unclear whether a similar coupling between 0.016 Hz and HFA exists.

In our previous case report, we revealed that the PAC between ISA and HFA preceded the seizure onset (SO) and proposed that ISA-HFA PAC might be a useful marker for seizure prediction24. In this study, we hypothesized that ISA-HFA PAC might precede the SO and help us discriminate between interictal and preictal states; this may have implications for clinical use in seizure prediction. To assess this, we extracted ISA using a 0.016–1 Hz bandpass filter and HFA using an 80–250 Hz bandpass filter. We used the synchronization index (SI)16 to investigate PAC.

## Methods

### Study setting and participants

The retrospective study was performed at Osaka University Hospital in Suita, Osaka, Japan, from July 2018 to July 2019. The study design was approved by the Ethics Committee of Osaka University Hospital (approval no. 19193) and informed consent was obtained using the opt-out method on our center’s website. We included seven patients with drug-resistant focal epilepsy who underwent intracranial electrode placement as part of a presurgical invasive EEG study (Table 1).

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/12/06/2020.11.07.20226258/T1)

Table 1. 
Clinical profile of the enrolled patients.

Seizure numbers are serial numbers of seizures in each participant.

SOC, seizure onset contact; MTLE, mesial temporal lobe epilepsy; FIAS, focal-impaired awareness seizure; FBTCS, focal-to-bilateral tonic-clonic seizure; PLE, parietal lobe epilepsy; OLE, occipital lobe epilepsy; FAS, focal aware seizure

* Focal resection surgery was not performed because it was impossible to detect the seizure onset zone.

### Intracranial electrodes

To acquire iEEG data, we used a combination of subdural grids (10, 20, or 30 contacts), strips (four or six contacts), and depth electrodes (six contacts) (Unique Medical Co. Ltd., Tokyo, Japan) that were placed by conventional craniotomy (details shown in eMethod).

### Data acquisition and preprocessing

We acquired iEEG signals using a 128-channel digital EEG system (EEG 2000; Nihon Kohden Corporation, Tokyo, Japan) at a sampling rate of 1 kHz and a time constant of 10 seconds. The BESA Research 6.0 software (BESA GmbH, Grafelfing, Germany) preprocessed the raw signals using a high-cut filter at 333 Hz to prevent aliasing and a 60-Hz notch filter to eliminate the AC line artifact and exported the data as a text file. This text file containing iEEG signals was imported into MATLAB R2020a (MathWorks, Natick, MA, USA) and iEEG signals were digitally re-referenced to a common average of all electrode contacts in each patient.

We saved the iEEG data every 60 minutes; therefore, one text file contained one 60-minute signal. We applied a bandpass filter to the entire 60-minute data to prevent edge-effect artifacts.

### Seizure onset and contacts

The SO was determined by conventional visual inspection of iEEG signals25. The contacts that showed initial epileptic changes immediately after the SO were related to the seizure onset zone (SOZ) and those that showed no initial epileptic changes immediately after the SO were determined as non-SOZ (nSOZ) contacts; we randomly chose the contacts in nSOZ contacts such that we had a similar number of nSOZ and SOZ contacts (Table 1) (eFigure1).

The time before the SO was the preictal state and the time after the SO was the ictal state. We randomly extracted the time from the iEEG data without seizures such that the count of time values was ten times that of the SO; we defined this data as the interictal state.

### Infraslow activities and high-frequency activities

To extract ISA and HFA, we used a 0.016–1 Hz and 80–250 Hz bandpass filter, respectively. A bandpass device with a two-way, least-square, finite-impulse response filter (pop_eegfiltnew.m from the EEGLAB toolbox, [https://sccn.ucsd.edu/eeglab/index.php](https://sccn.ucsd.edu/eeglab/index.php)) was applied to the iEEG signals. The oscillations of the 0.016–1 Hz signals were representative of ISA (eFigure2B and eFigure3B). We calculated the amplitude of HFA, which was an envelope of 80–250 Hz signals, in combination with the Hilbert transformation.

We combined a permutation test with a family-wise error (FWE)-corrected threshold to detect significant HFA changes. The oscillations under −1 mV or over +1 mV2,13 were defined as significant ISA changes (eMethod, eFigure2A, and eFigure4A, 4B).

### PAC analysis

An SI16 was used to measure the strength of coupling between the HFA amplitude and the ISA phase. Hilbert transformation was performed on the bandpass-filtered signals to obtain complex-valued analytic signals [Z(t)]. The amplitude [A(t)] and phase [φ(t)] were calculated from the complex-valued signals using Equation 1. ![Formula][1]</img>  The ISA phase was calculated using the angle of the Hilbert transformation in the 0.016–1 Hz bandpass-filtered signal. The HFA amplitude was calculated using the squared magnitude of the Hilbert transformation in the 80–250 Hz bandpass filtered signal. Then, the phase of this amplitude was computed using the Hilbert transformation. The SI was calculated using Equation 2. ![Formula][2]</img>  where “n” is the number of data points. Since SI is a complex number, we used the magnitude of SI (SIm) in our calculations. SIm varies between 0 and 1, with 0 indicating completely desynchronized phases and 1 indicating perfectly synchronized phases.

The preferred phase of synchronization (SIp) was calculated by the arctan (image [SI] / real [SI]). SIp varies between −180° and +180°.

### Phase-conditioned analysis

We calculated the mean vector and performed the Rayleigh test to evaluate the nonuniformity of SIp using the CircStat toolbox26. Next, to identify the ISA phase to which the HFA amplitude was coupled, we calculated the average oscillations of ISA and the normalized amplitude of HFA within each ISA phase bin of 30°: −180° – −150°, −150° – −120°, …, and 150° – 180°.

### Correlation analysis

Using all the implanted contacts from both SOZ and nSOZ, we calculated Pearson correlation coefficients between SIm and the normalized amplitude of HFA.

### Classification

We implemented a classification to distinguish between the three states using HFA, ISA, and PAC of SOZ or nSOZ. We also differentiated between SOZ and nSOZ was using the HFA, ISA, and PAC values of the interictal, preictal, and ictal states. The receiver-operating characteristic (ROC) curve and its area under the curve (AUC) were used to compare the performance of different classifiers.

### Statistical analysis

We compared SOZ with nSOZ within HFA, ISA, or PAC, and conducted paired comparisons of HFA with ISA, ISA with PAC, and HFA with PAC within SOZ or nSOZ using the Wilcoxon rank-sum test. The AUC was also compared using the Wilcoxon signed-rank test. The results were corrected using the Bonferroni correction for multiple comparisons.

## Results

Overall, we recorded 15 seizures, all of which were focal. We observed focal to bilateral tonic-clonic seizure (FBTCS) in 12 episodes (80%). There were 43 SOZ contacts (Table 1).

Representative seizures of S1 in P1 indicated that ISA and HFA changes occurred after the SO, whereas before the SO, there were no clear changes in them. However, SIm started to increase even a few minutes before the SO (eFigures 2 and 3).

### Profiles of HFA, ISA, and PAC

Within a 30-second time-window in the three states, we averaged and compared the HFA, ISA, and ISA-HFA PAC (SIm) of SOZ and nSOZ (Figure1A, 1B, and 1C, eMethod). In ictal-HFA and ictal-ISA, SOZ showed significantly higher values than nSOZ (corrected p < 0.001); however, there were no significant differences in ictal-PAC between SOZ and nSOZ. Both ictal-HFA and the ISA of SOZ achieved their maximum value, which was significantly higher than the interictal- and preictal-HFA and the ISA of SOZ (corrected p < 0.001). On the other hand, the ictal-PAC of SOZ was significantly lower than the preictal-PAC of SOZ (corrected p < 0.001). The preictal-PAC of SOZ achieved its maximum value and was significantly higher than the interictal-PAC of SOZ (corrected p < 0.001). Similarly, in the preictal stage, SOZ-ISA and nSOZ-PAC were significantly higher than those in the interictal stage; therefore, we believed that the PAC of SOZ and nSOZ and the ISA of SOZ might be potential biomarkers for seizure prediction.

Contrary to our expectation, the interictal- and preictal-HFA of SOZ were significantly lower than that of nSOZ (corrected p < 0.001 and = 0.005 for each).

### The time point when significant changes occurred

We compared the time points when HFA, ISA, and PAC of SOZ or nSOZ showed significant changes related to seizures (criteria for significant changes in eMethod and eFigure 4) (Figure 1D). In HFA and ISA, the SOZ occurred significantly earlier than nSOZ (corrected p < 0.001), whereas, in PAC, there were no significant differences between SOZ and nSOZ (corrected p = 0.31). Similarly, there were no significant differences between SOZ-HFA and SOZ-ISA (corrected p = 0.077); however, SOZ-PAC changed at a significantly earlier time (−93.82 s, 95% confidence interval: −116.40 – −71.23 s) than SOZ-HFA and SOZ-ISA (corrected p < 0.001). The average time of HFA and ISA was immediately after the occurrence of SO, whereas the average time of PAC was before the SO. We believe that these results also support the feasibility of ISA-HFA PAC for seizure prediction.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/06/2020.11.07.20226258/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/12/06/2020.11.07.20226258/F1)

Figure 1. 
Profiles of HFA, ISA, and PAC.

Paired results calculated from contacts placed in the SOZ and nSOZ are shown. (A) HFA amplitude of ictal-SOZ achieved its maximum value and was significantly higher than that of interictal-SOZ, preictal-SOZ, and ictal-nSOZ. (B) ISA oscillations of ictal-SOZ achieved their maximum values and were significantly higher than those of interictal-SOZ, preictal-SOZ, and ictal-nSOZ. (C) For ISA-HFA PAC, the SIm of preictal-SOZ achieved its maximum value and was significantly higher than that of interictal-SOZ and ictal-SOZ. Error bars in A, B, and C indicate the standard deviation (SD). (D) Time points at which we compared significant seizure-related changes, as shown in the box and whisker plots. For HFA and ISA, SOZ was significantly before nSOZ, and there were no significant differences between the times of SOZ and nSOZ-PAC. The average time of PAC was earlier than that of SO (0 seconds). (E) We observed the rate at which significant changes occurred. For HFA and PAC, almost all contacts of SOZ and nSOZ showed significant changes. For ISA, 89% of the SOZ contacts and 67% of the contacts of nSOZ showed significant changes.

In panels A, B, C, and D, the average values and SD are indicated as figures.

*p < 0.05, **p < 0.01, Wilcoxon rank-sum test, Bonferroni corrected.

### Rate at which significant changes were observed

Almost all the SOZ and nSOZ values showed significant changes in HFA and PAC. In ISA, 89% of SOZ and 67% of nSOZ showed significant changes. In the SOZ, there were no significant differences between HFA, ISA, and PAC (Figure 1E).

### Phase-based analysis

Figure 1 shows that in the ictal states, HFA- and ISA-SOZ had significantly higher values than nSOZ, and they changed earlier than nSOZ. However, PAC showed no significant differences between SOZ and nSOZ; therefore, we used phase-based analyses to investigate the differences between SOZ- and nSOZ-PAC.

Using a 30-second duration in the interictal or preictal states, we calculated the mean vectors of SIp in SOZ or nSOZ (Figure 2A); there were no similarities between SOZ and nSOZ in these two states. We observed significant nonuniformity only in the nSOZ-interictal state (p = 0.001, Rayleigh test).

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/06/2020.11.07.20226258/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/12/06/2020.11.07.20226258/F2)

Figure 2. 
Phase-based analyses of PAC

(A) We calculated average values of the SIp of SOZ or nSOZ in the interictal and preictal states. There were no similarities between the SOZ and nSOZ. The angle, length, and standard deviation of the mean vector are indicated as θ, R, and SD and p-values calculated by the Rayleigh test are shown. (B) We sequentially plotted the angles of the mean vector from −5 minutes to 1 minute. We observed continuous and periodic shapes in SOZ and discontinuous shapes in nSOZ. (C) The phase-tuning HFA-normalized amplitude was at its peak at the trough of the ISA phase and the bottom during the peak of the ISA phase in the preictal states; this trend was reversed in nSOZ. In the ictal states, these clear patterns became obscure. The error bars indicate 95% confidence intervals.

We sequentially plotted the angle (°) of mean vectors calculated from SIp in SOZ or nSOZ from −5 minutes to +1 minute around the SO (Figure 2B). In the SOZ, we documented continuous and periodic angle changes from −180° to +180°, while angle changes in nSOZ showed discontinuous and collapsed periodicity.

Figure 2C depicts the phases of ISA coupled with the amplitude of HFA, where “Preictal” indicates an average of 5 minutes (−5 to 0 minutes) and “Ictal” indicates an average of 1 minute (0 to 1 minute). In preictal-SOZ states, the HFA amplitude peaked at the trough of the ISA oscillation and reached its bottom-most value at the peak of the ISA oscillation, whereas we found an opposite trend in the nSOZ state––this was a contrasting result. These clear patterns were disrupted in the ictal state.

### Correlation between HFA and PAC

HFA increased in the ictal states, whereas PAC increased and displayed clear patterns in the preictal states rather than in the ictal states. After observing these contrasting results, we investigated the correlation between HFA and PAC.

There was no significant correlation between the preictal-HFA and preictal-PAC, and ictal-HFA and ictal-PAC; however, the correlation between the ictal-HFA and preictal-PAC was significantly positive (corrected p < 0.001) (Figure 3A). We calculated correlation coefficients (r) and their corrected p values in combination with sequential HFA and sequential PAC from −5 minutes to +2 minutes around the SO and displayed them as a matrix (Figure 3B). We observed an obvious positive correlation between ictal-HFA (after 0 minutes) and preictal-PAC (before 0 minutes). This demonstrates that contacts that indicated increased PAC before SO showed more of an increase in HFA after SO.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/06/2020.11.07.20226258/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2020/12/06/2020.11.07.20226258/F3)

Figure 3. 
Correlation between HFA and PAC.

(A) There was no significant correlation between HFA and PAC in the preictal or ictal state. However, ictal-HFA and preictal-PAC showed a significantly positive correlation (corrected p < 0.001). (B) Correlation coefficients and their corresponding significant p values (corrected p < 0.01) are shown for all combinations of sequential HFA and sequential PAC, from −5 minutes to +2 minutes around the SO (0 minutes). We observed a significant positive correlation in the range after the SO of HFA and before SO of PAC.

### Classification

To evaluate whether our method was accurately predicting seizures, we compared the preictal states with interictal states; to evaluate seizure detection, we compared the ictal states with interictal, and the ictal states with the preictal states (Figure 4A). After classifying the preictal and interictal states, we found that the AUC of SOZ-PAC was at its maximum and significantly higher than those of SOZ-ISA, SOZ-HFA, and nSOZ-PAC (corrected p < 0.001). SOZ-ISA showed the best performance in terms of seizure detection, better than SOZ-HFA, SOZ-PAC, and nSOZ-ISA (corrected p < 0.001, except for nSOZ-ISA in ictal vs. preictal; corrected p = 0.46).

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/06/2020.11.07.20226258/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2020/12/06/2020.11.07.20226258/F4)

Figure 4. 
Average receiver operating characteristic curves of the results following stage classification.

(A); Preictal vs. interictal, ictal vs. interictal, and ictal vs. preictal stages were classified using HFA, ISA, and PAC of SOZ or nSOZ. SOZ-PAC showed the highest performance between the preictal and interictal states, while SOZ-ISA had higher potential to separate the ictal and interictal, and ictal and preictal states. (B) On classifying SOZ and nSOZ using HFA, ISA, and PAC, we found that HFA had the maximum AUC in the interictal and preictal states. ISA performed best in the ictal states.

Finally, we classified SOZ and nSOZ using HFA, ISA, and PAC (Figure 4B). In the interictal states, PAC and ISA showed a chance level; however, HFA performed better than the others (corrected p < 0.001). Similarly, HFA showed the best performance in the preictal states (corrected p < 0.001). However, in the ictal states, ISA showed the best performance (corrected p < 0.001).

## Discussion

HFA and ISA are often observed during seizures13,27-29 and they are known to be useful biomarkers to detect SOZ5,30,31. In this study, SOZ achieved more ictal-HFA and ictal-ISA than nSOZ. The onset of ictal-ISA is typically earlier than that of the HFA13,31; however, we did not observe any differences between the onsets of SOZ-HFA and SOZ-ISA in this study. We defined the onset of ISA at the time point at > +1 mV or less −1 mV and the onset of ISA as the initial time of ISA reported in previous studies13. Therefore, our ISA onset time might be later than the values reported in previous studies. The percentage of SOZ-ISA in this study (88.89%) was concordant with that of a previous report (87%)5.

Our proposed PAC method brought new insight into the relationship between ISA and HFA. In this study, ISA-HFA PAC reached its maximum in the preictal state, began changing before SO, and showed no difference between SOZ and nSOZ. In the phase-based analyses, we observed differences between preictal-SOZ and preictal-nSOZ. The HFA amplitude in SOZ was tuned according to the trough of the ISA, which was concordant with the results of previous reports19; however, in nSOZ, the HFA amplitude was tuned at the peak of the ISA. These contrasting trends might be useful for differentiating between SOZ and nSOZ in preictal states.

During the seizure-dependent HFA increase, the PAC using a lower frequency (except for ISA) simultaneously increases19,20,24. However, in this study, PAC with ISA increased before the increase in HFA, which was concordant with the profile of physiological PAC22,23. Moreover, the PAC before SO was positively correlated with the later HFA; therefore, ISA-HFA PAC might play an essential role in inducing an HFA burst during or just before seizures.

Various features that are calculated using EEG have been proposed for seizure prediction32, and machine learning methods have been used to detect preictal states33,34. We observed a statistical significance between the interictal and preictal stages in SOZ-PAC, nSOZ-PAC, and SOZ-ISA (Figure 1); the AUC of preictal vs. interictal were big in this order (Figure 4). The AUC of SOZ-PAC was 0.926; therefore, SOZ-PAC could be a potential, novel biomarker for seizure prediction. Using ISA-HFA PAC, it is possible to achieve accurate and high-performing markers without complex algorithms and long training times. Accurate seizure prediction carries advantages for responsive neurostimulation with implantable devices35,36 and ultimately allows for early treatment and prevention of seizures.

Several features and algorithms have been used to detect seizures32, including the PAC of the β band21. Previously, ISA and HFA were simultaneously observed during seizures5,13,31, and our study added to the results by demonstrating that SOZ-ISA had a higher performance in differentiating between ictal states than SOZ-HFA. ISA and HFA are known to be useful for detecting SOZ5,13,31; we showed that in the ictal state, ISA performed better than HFA for SOZ differentiation.

Contrary to one expectation in our study, SOZ-HFA values were significantly lower than nSOZ-HFA in the interictal states. Reportedly, interictal-HFO is less correlated with the ictal onset zone31; therefore, interictal-HFA may be suppressed more in SOZ than in nSOZ. This surprising result further contributed to the interictal discrimination between SOZ and nSOZ.

This study has some limitations. Firstly, we obtained the result using mostly iEEG, and it is unclear whether the same results can be obtained using scalp EEG, which is the most commonly used method. Secondly, we evaluated only focal-onset seizures; therefore, it remains unclear whether generalized-onset seizures, such as absence or myoclonic seizures, would indicate the same results. Thirdly, we performed off-line analyses; however, in clinical situations, on-line analysis is fundamental for confirming the feasibility of seizure prediction using ISA-HFA PAC. We believe that further large-scale investigation is needed to strengthen our results.

## Data Availability

All data that were generated by or analyzed in this study are available from the corresponding authors upon reasonable request and after additional ethics approvals regarding the data provision to individual institutions.

## Author Contributions

H.H. conceived the study, collected the data, created the MATLAB program, analyzed the data, created all figures, and was primarily responsible for writing the manuscript. H.M.K., N.T., S.O., H.K., and M.H. performed the epileptic surgery. All authors clinically cared for and evaluated the patient. H.M.K., and T.Y. advised H.H. on scientific matters. H.K. and M.H. supervised this study. All authors have reviewed the manuscript.

## Conflict of Interest Disclosures

No author has any conflict of interest to disclose.

## Funding/Support

This study was supported by Grants-in-Aid for Early-Career Scientists (KAKENHI; grant no. 18K18366), which is funded by the Japan Society for the Promotion of Science (JSPS; Tokyo, Japan).

## Supplementary Online Content

### eMethods

#### Intracranial electrodes

The diameter of each contact was 3 or 5 mm, and the inter-contact distances were 5, 7, or 10 mm for the grid and strip electrodes. The diameter was 1.5 mm and the inter-contact distance was 5 mm for the depth electrodes.

#### Topographies of HFA, ISA, and PAC

The 0.016–1 Hz bandpass-filtered signals for HFA, 80–250 Hz bandpass-filtered signals for ISA, and SIm signals for PAC were analyzed from 5 minutes before to 5 minutes after the SO. The power of the HFA was the square of the HFA amplitude (envelope) and the series was normalized by HFA power in the initial 1 minute (from -5 to -4 minutes before the SO). The HFA normalized power is shown as a topography (eFigure 3A). We defined the active seizure-related ISA as oscillations of the 0.016–1 Hz bandpass-filtered signals less than -1 mV or more than 1 mV; the active ISAs were displayed in grayscale topography (eFigure 3B). Statistically significant SIm values, acquired by a combination of a bootstrapped technique and familywise error (FWE)-corrected threshold (see below), were indicated in the topography (eFigure 4C).

#### Bootstrapped technique and FWE-corrected threshold

For the statistical assessment of SIm, we randomly shifted the phase-time series of the HFA amplitude and calculated the bootstrapped SIm (SImb) using the lower frequency phase. We repeated this procedure 1000 times to create the distribution of SImb1. The maximum values of SImb were stored at each iteration in its distribution. The values at 95% of the distribution of the maximum FWE-corrected threshold have been applied to the observed SIm to obtain the solution to multiple comparisons2. SIm values over the FWE-corrected threshold were statistically significant.

#### Profiles of HFA, ISA, and PAC (Figure 1A, 1B, and 1C)

For within-HFA comparisons, we used the amplitude of HFA signals, which were an envelope of 80–250 Hz bandpass-filtered signals. For within-ISA comparisons, the absolute values of oscillation of 0.016–1 Hz bandpass-filtered signals were used because the ISAs were slow-negative or slow-positive activities. For within-PAC comparisons, we used the SIm between the ISA phase and HFA amplitude.

We averaged values with a 30-second time window. The time interval of 30 seconds before the SO corresponded to preictal states, and that of 30 seconds after the SO corresponded to ictal states. The time of onset of the interictal states was randomly selected from iEEG data with no seizures.

#### Criteria for significant HFA, ISA, and PAC changes

To investigate the time when significant seizure-related changes occurred, we evaluated seizure-related signals from 5 minutes before to 2 minutes after the SO.

To detect significant changes in HFA, we used a permutation test3 to compare the initial 10-second data and the next 10-second data of HFA-normalized amplitude, which were acquired following normalization using the initial 10-second data. Each permutation test produced a set of differences between the initial 10-second data and the next sequential 10-second data. The maximum value of the differences from each permutation test was stored, and the values at 95% of the distribution of these maximum values were taken as the FWE-corrected threshold. The values above the FWE-corrected threshold were statistically significant2 and the time points when the HFA amplitude first crossed over the FWE-corrected threshold were defined as significant changes in HFA (eFigure 4A). However, if the HFA amplitude temporally crossed over the FWE-corrected threshold due to a single-spiking activity, we excluded it from the analysis.

For ISA, the time when iEEG signals first crossed under -1 mV or over +1 mV4,5 were defined as significant changes in ISA (eFigure 4B).

For PAC, we used significant SIm values, which were obtained in combination with the bootstrapped technique and FWE-corrected threshold. If the significant SIm formed a cluster, i.e., significant SIm was seen several times in 30-second time duration, we considered the first time point as the significant time of PAC changes (eFigure 4C).

![eFigure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/06/2020.11.07.20226258/F5.medium.gif)

[eFigure 1.](http://medrxiv.org/content/early/2020/12/06/2020.11.07.20226258/F5)

eFigure 1. 
Contacts related to SOZ and nSOZ (P5-S3).

(A) Illustration of implanted electrodes and the brain are displayed. The numbers correspond to the contacts’ number. (B) The high-intensity lesion in the right occipital lobe is shown on fluid-attenuated inversion recovery MRI (red wedge arrows) (B). This lesion is pointed by the arrow in the panel A. (C) SOZ contacts including A31, A32, A36, and A37 that were placed on the lesion showed initial infraslow activities immediately after the SO (0 s), followed by low-voltage fast waves. At that time, nSOZ contacts including A40, A41, A54, and A55 showed no epileptic changes. The number of SOZ and nSOZ contacts are same.

![eFigure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/06/2020.11.07.20226258/F6.medium.gif)

[eFigure 2.](http://medrxiv.org/content/early/2020/12/06/2020.11.07.20226258/F6)

eFigure 2. 
A representative seizure of S1 in P1.

Representative multimodal signals related to seizure (S1 in P1 shown in Table 1) were calculated from one contact placed in the seizure onset zone (SOZ). From top to bottom column, raw iEEG signals in black (A), 0.016–1 Hz band-pass filtered signal as ISA in green (B), 80–250 Hz band-pass filtered signals as HFA in red (C), and SIm as ISA-HFA PAC in blue (D) are displayed. The 0 min corresponded to the SO. In iEEG signals, ± 1 mV were indicated as horizontal dotted lines, which were threshold for significant ISA related to seizures. After the SO, obvious activities over ± 1 mV were observed in iEEG signals, which were seizure-related ISA (A), and the same shape of signals was confirmed in ISA (B). Therefore, we knew that the 0.016–1 Hz band-pass filter enabled us to extract ISA activities. At the same time, an increase in HFA was observed, which was seizure-related HFA. However, before the SO, there were no obvious activities of either ISA or HFA. However, SIm started to increase from 15 min before the SO and reached a peak at the SO.

![eFigure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/06/2020.11.07.20226258/F7.medium.gif)

[eFigure 3.](http://medrxiv.org/content/early/2020/12/06/2020.11.07.20226258/F7)

eFigure 3. 
Topographies of ISA, HFA, and PAC of S1 in P1.

The 0 minute-time point corresponded to the SO. All contacts of P1 were used; in eFigure 2 we only used the A2 contact that was placed in the SOZ. (A) The HFA-normalized power began increasing from the A2 contacts after SO and spreading to other contacts. (B) Active seizure-related ISAs, which were defined as more than 1 mV or less than -1 mV, are indicated on grayscale. There were obvious activities of both HFA and ISA after SO; however, no clear activities of either HFA or ISA were observed before SO. However, (C) we observed significant SIm with ISA-HFA, which was obtained using the FWE-corrected threshold that was applied to raw SIm, even before SO.

![eFigure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/06/2020.11.07.20226258/F8.medium.gif)

[eFigure 4.](http://medrxiv.org/content/early/2020/12/06/2020.11.07.20226258/F8)

eFigure 4. 
Period of significant change in S1 of P1.

The seizure-related data from 5 min before to 2 min after the SO (0 min) of S1 in P1 are shown. (A) The HFA-normalized amplitude (red line) and FWE-corrected threshold (dashed blue line) are indicated. At 39.40 seconds, the HFA-normalized amplitude crossed over the FWE-corrected threshold for the first time, making this timepoint significant. (B) The iEEG signals (green lines) and range from -1 mV to +1 mV (gray mesh) are indicated. The time point of 7.80 seconds was the significant time of ISA change when the iEEG signal first crossed under -1 mV. (C) Statistically significant SIm appeared at -171.31 seconds, which was the significant time of PAC changes.

## Footnotes

*   Hui Ming Khoo khoo{at}nsurg.med.osaka-u.ac.jp, Takufumi Yanagisawa tyanagisawa{at}nsurg.med.osaka-u.ac.jp, Naoki Tani n-tani{at}nsurg.med.osaka-u.ac.jp, Satoru Oshino s-oshino{at}nsurg.med.osaka-u.ac.jp, Haruhiko Kishima hkishima{at}nsurg.med.osaka-u.ac.jp, Masayuki Hirata mhirata{at}ndr.med.osaka-u.ac.jp

*   New eFigure was added and other eFigures were revised. New; eFigure 1 Revised; eFigure 2-4

*   Received November 7, 2020.
*   Revision received December 6, 2020.
*   Accepted December 6, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  1.Emami A, Kunii N, Mtasuo T, Shinozaki T, Kawai K, Takahashi HJNC. Seizure detection by convolutional neural network-based analysis of scalp electroencephalography plot images. 2019:101684.
    
    

2.  2.Modur PN, Vitaz TW, Zhang S. Seizure localization using broadband EEG: comparison of conventional frequency activity, high frequency oscillations and infraslow activity. Journal of clinical neurophysiology: official publication of the American Electroencephalographic Society. 2012;29(4):309.
    
    

3.  3.Imamura H, Matsumoto R, Inouchi M, et al. Ictal wideband ECoG: direct comparison between ictal slow shifts and high frequency oscillations. Clinical Neurophysiology. 2011;122(8):1500–1504.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.clinph.2010.12.060&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21353632&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 

4.  4.Zijlmans M, Jiruska P, Zelmann R, Leijten FS, Jefferys JG, Gotman J. High - frequency oscillations as a new biomarker in epilepsy. Annals of neurology. 2012;71(2):169–178.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.22548&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22367988&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 

5.  5.Ikeda A, Takeyama H, Bernard C, et al. Active direct current (DC) shifts and “Red slow”: two new concepts for seizure mechanisms and identification of the epileptogenic zone. Neurosci Res. 2020;156:95–101.
    
    

6.  6.Jirsch J, Urrestarazu E, LeVan P, Olivier A, Dubeau F, Gotman J. High-frequency oscillations during human focal seizures. Brain. 2006;129(6):1593–1608.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awl085&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16632553&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000238224900024&link_type=ISI) 

7.  7.Ayoubian L, Lacoma H, Gotman J. Automatic seizure detection in SEEG using high frequency activities in wavelet domain. Medical engineering & physics. 2013;35(3):319–328.
    
    

8.  8.Ikeda A, Terada K, Mikuni N, et al. Subdural recording of ictal DC shifts in neocortical seizures in humans. Epilepsia. 1996;37(7):662–674.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1528-1157.1996.tb00631.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8681899&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996UY99900011&link_type=ISI) 

9.  9.Hashimoto H, Hasegawa Y, Araki T, et al. Non-invasive detection of language-related prefrontal high gamma band activity with beamforming MEG. Sci Rep. 2017;7(1):14262.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-017-14452-3&link_type=DOI) 

10. 10.Matsumoto A, Brinkmann BH, Matthew Stead S, et al. Pathological and physiological high-frequency oscillations in focal human epilepsy. Journal of neurophysiology. 2013;110(8):1958–1964.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/jn.00341.2013&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23926038&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000325815000020&link_type=ISI) 

11. 11.Modur PN, Zhang S, Vitaz TW. Ictal high-frequency oscillations in neocortical epilepsy: implications for seizure localization and surgical resection. Epilepsia. 2011;52(10):1792–1801.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1528-1167.2011.03165.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21762451&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000296067000013&link_type=ISI) 

12. 12.Jacobs J, LeVan P, Châtillon C-É, Olivier A, Dubeau F, Gotman J. High frequency oscillations in intracranial EEGs mark epileptogenicity rather than lesion type. Brain. 2009;132(4):1022–1037.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awn351&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19297507&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000265269600025&link_type=ISI) 

13. 13.Kanazawa K, Matsumoto R, Imamura H, et al. Intracranially recorded ictal direct current shifts may precede high frequency oscillations in human epilepsy. Clinical Neurophysiology. 2015;126(1):47–59.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.clinph.2014.05.028&link_type=DOI) 

14. 14.Inoue T, Inouchi M, Matsuhashi M, et al. Interictal Slow and High-Frequency Oscillations: Is it an Epileptic Slow or Red Slow? Journal of Clinical Neurophysiology. 2019;36(2):166–170.
    
    

15. 15.Canolty RT, Edwards E, Dalal SS, et al. High gamma power is phase-locked to theta oscillations in human neocortex. Science. 2006;313(5793):1626–1628.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMTMvNTc5My8xNjI2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDYvMjAyMC4xMS4wNy4yMDIyNjI1OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

16. 16.Cohen MX. Assessing transient cross-frequency coupling in EEG data. J Neurosci Methods. 2008;168(2):494–499.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jneumeth.2007.10.012&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18061683&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000253824400023&link_type=ISI) 

17. 17.Nariai H, Matsuzaki N, Juhász C, et al. Ictal high-frequency oscillations at 80?200 Hz coupled with delta phase in epileptic spasms. Epilepsia. 2011;52(10):e130–e134.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1528-1167.2011.03263.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21972918&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000296067000003&link_type=ISI) 

18. 18.Iimura Y, Jones K, Takada L, et al. Strong coupling between slow oscillations and wide fast ripples in children with epileptic spasms: Investigation of modulation index and occurrence rate. Epilepsia. 2018;59(3):544–554.
    
    

19. 19.Ibrahim GM, Wong SM, Anderson RA, et al. Dynamic modulation of epileptic high frequency oscillations by the phase of slower cortical rhythms. Experimental neurology. 2014;251:30–38.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.expneurol.2013.10.019&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24211781&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 

20. 20.Hashimoto H, Khoo HM, Yanagisawa T, et al. Visualization of Phase-Amplitude Coupling Using Rhythmic High-Frequency Activity. BioRxiv. 2020.
    
    

21. 21.Edakawa K, Yanagisawa T, Kishima H, et al. Detection of Epileptic Seizures Using Phase-Amplitude Coupling in Intracranial Electroencephalography. Sci Rep-Uk. 2016;6.
    
    

22. 22.Yanagisawa T, Yamashita O, Hirata M, et al. Regulation of Motor Representation by Phase-Amplitude Coupling in the Sensorimotor Cortex. Journal of Neuroscience. 2012;32(44):15467–15475.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjExOiIzMi80NC8xNTQ2NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzEyLzA2LzIwMjAuMTEuMDcuMjAyMjYyNTguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

23. 23.Hashimoto H, Takahashi K, Kameda S, et al. Motor and Sensory Cortical Processing of Neural Oscillatory Activities revealed by Human Swallowing using Intracranial Electrodes. bioRxiv. 2020:2020.2007.2021.213868.
    
    

24. 24.Hashimoto H, Khoo HM, Yanagisawa T, et al. Coupling between infraslow activities and high-frequency oscillations precedes seizure onset. Epilepsia Open. 2020;5(3):501–506.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/epi4.12425&link_type=DOI) 

25. 25.Perucca P, Dubeau F, Gotman J. Intracranial electroencephalographic seizure-onset patterns: effect of underlying pathology. Brain. 2014;137(Pt 1):183–196.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awt299&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24176980&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000330434000018&link_type=ISI) 

26. 26.Berens P. CircStat: a MATLAB toolbox for circular statistics. J Stat Softw. 2009;31(10):1–21.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v031.i10&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25402854&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 

27. 27.Ikeda A, Taki W, Kunieda T, et al. Focal ictal direct current shifts in humanepilepsy as studied by subdural and scalp recording. Brain. 1999;122(5):827–838.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/122.5.827&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10355669&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000080528000005&link_type=ISI) 

28. 28.Akiyama T, Chan DW, Go CY, et al. Topographic movie of intracranial ictal high-frequency oscillations with seizure semiology: epileptic network in Jacksonian seizures. Epilepsia. 2011;52(1):75–83.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1528-1167.2010.02776.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21070217&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 

29. 29.Ochi A, Otsubo H, Donner EJ, et al. Dynamic changes of ictal high-frequency oscillations in neocortical epilepsy: using multiple band frequency analysis. Epilepsia. 2007;48(2):286–296.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1528-1167.2007.00923.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17295622&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 

30. 30.Wang S, Wang IZ, Bulacio JC, et al. Ripple classification helps to localize the seizure-onset zone in neocortical epilepsy. Epilepsia. 2013;54(2):370–376.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1528-1167.2012.03721.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23106394&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000314750200021&link_type=ISI) 

31. 31.Wu S, Kunhi Veedu HP, Lhatoo SD, Koubeissi MZ, Miller JP, Lüders HO. Role of ictal baseline shifts and ictal high-frequency oscillations in stereo -electroencephalography analysis of mesial temporal lobe seizures. Epilepsia. 2014;55(5):690–698.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24725106&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 

32. 32.Ramgopal S, Thome-Souza S, Jackson M, et al. Seizure detection, seizure prediction, and closed-loop warning systems in epilepsy. Epilepsy & behavior. 2014;37:291–307.
    
    

33. 33.D’Alessandro M, Vachtsevanos G, Esteller R, et al. A multi-feature and multi-channel univariate selection process for seizure prediction. Clinical neurophysiology. 2005;116(3):506–516.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.clinph.2004.11.014&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15721064&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000227895400003&link_type=ISI) 

34. 34.Petrosian A, Prokhorov D, Homan R, Dasheiff R, Wunsch II D. Recurrent neural network based prediction of epileptic seizures in intra-and extracranial EEG. Neurocomputing. 2000;30(1-4):201–218.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0925-2312(99)00126-5&link_type=DOI) 

35. 35.Sun FT, Morrell MJ, Wharen RE. Responsive cortical stimulation for the treatment of epilepsy. Neurotherapeutics. 2008;5(1):68–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.nurt.2007.10.069&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18164485&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000252532300008&link_type=ISI) 

36. 36.Cook MJ, O’Brien TJ, Berkovic SF, 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. The Lancet Neurology. 2013;12(6):563–571.
    
    

## eReferences

1.  1.Cohen MX. Assessing transient cross-frequency coupling in EEG data. J Neurosci Methods. 2008;168(2):494–499.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jneumeth.2007.10.012&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18061683&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000253824400023&link_type=ISI) 

2.  2.Cohen MX. Analyzing neural time series data: theory and practice. MIT press; 2014.
    
    

3.  3.Maris E, Oostenveld R. Nonparametric statistical testing of EEG- and MEG-data. J Neurosci Methods. 2007;164(1):177–190.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jneumeth.2007.03.024&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17517438&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248170300019&link_type=ISI) 

4.  4.Modur PN, Vitaz TW, Zhang S. Seizure localization using broadband EEG: comparison of conventional frequency activity, high frequency oscillations and infraslow activity. Journal of clinical neurophysiology: official publication of the American Electroencephalographic Society. 2012;29(4):309.
    
    

5.  5.Ikeda A, Taki W, Kunieda T, et al. Focal ictal direct current shifts in humanepilepsy as studied by subdural and scalp recording. Brain. 1999;122(5):827–838.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/122.5.827&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10355669&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F06%2F2020.11.07.20226258.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000080528000005&link_type=ISI)

 [1]: /embed/graphic-2.gif
 [2]: /embed/graphic-3.gif