EEG Reactivity Predicts Individual Efficacy of Vagal Nerve Stimulation in Intractable Epileptics

Background: Chronic vagal nerve stimulation (VNS) is a well-established non-pharmacological treatment option for drug-resistant epilepsy. This study sought to develop a statistical model for prediction of VNS efficacy. We hypothesized that reactivity of the electroencephalogram (EEG) to external stimuli measured during routine preoperative evaluation differs between VNS responders and non-responders.

Materials and Methods: Power spectral analyses were computed retrospectively on pre-operative EEG recordings from 60 epileptic patients with VNS. Thirty five responders and 25 non-responders were compared on the relative power values in four standard frequency bands and eight conditions of clinical assessment—eyes opening/closing, photic stimulation, and hyperventilation. Using logistic regression, groups of electrodes within anatomical areas identified as maximally discriminative by n leave-one-out iterations were used to classify patients. The reliability of the predictive model was verified with an independent data-set from 22 additional patients.

Results: Power spectral analyses revealed significant differences in EEG reactivity between responders and non-responders; specifically, the dynamics of alpha and gamma activity strongly reflected VNS efficacy. Using individual EEG reactivity to develop and validate a predictive model, we discriminated between responders and non-responders with 86% accuracy, 83% sensitivity, and 90% specificity.

Conclusion: We present a new statistical model with which EEG reactivity to external stimuli during routine presurgical evaluation can be seen as a promising avenue for the identification of patients with favorable VNS outcome. This novel method for the prediction of VNS efficacy might represent a breakthrough in the management of drug-resistant epilepsy, with wide-reaching medical and economic implications.

Introduction

Resective surgery is currently the best therapeutic option for treatment of patients with drug-resistant epilepsy, but a substantial number of intractable patients remains who are ineligible for such treatment or for whom resective surgery fails to abolish seizures. Chronic vagal nerve stimulation (VNS) has become a well-established alternative, offering a palliative method of treatment for drug-resistant epilepsy; it rarely results in complete seizure freedom (~5% of treated patients), but provides substantial (≥50%) seizure reduction in 50–60% of individuals. Unfortunately, however, seizure frequency remains unchanged after VNS therapy in ~25% of patients (1, 2).

Identifying individuals who will benefit from VNS therapy prior to the implantation of the VNS device would improve patient selection, minimize unnecessary surgical procedures, and reduce associated financial expenses dramatically; and yet there exists no method with which to predict individual efficacy pre-intervention (2). Achieving a pre-operative classification of individual patients as VNS responders or non-responders (i.e., patients with ≥50% or <50% seizure reduction, respectively) would represent a major breakthrough in the treatment of drug-resistant epilepsy.

It is presumed that VNS increases seizure threshold by activating neuronal networks in the thalamus and other limbic structures (3, 4), but the precise mechanism of VNS action is not yet understood fully. Both synchronization and desynchronization of the electroencephalogram (EEG) has been proposed as a possible mechanism behind the antiepileptic effect of VNS (5, 6), and recent neurophysiological studies focusing on EEG parameters lend support to this: Fraschini et al. report a significant correlation between VNS-induced global desynchronization in gamma bands and positive clinical outcome in temporal lobe epilepsy patients (7). Similarly, Bodin et al. revealed a lower level of global EEG synchronization in delta and alpha frequency bands during the ON phase of VNS in responders (8). Theoretically, differential alterations in brain rhythms from VNS therapy between responders and non-responders might reflect inter-individual variability in the (non-specific) susceptibility of EEG to be synchronized or desynchronized by external stimulation. It follows that differences in this susceptibility might underlie individual VNS efficacy.

We tested the hypothesis that EEG reactivity to standard external stimuli used during routine pre-operative EEG assessment differs between VNS responders and non-responders, with the aim of developing a reliable statistical model for prediction of VNS efficacy.

Materials and Methods

Study Design

We performed retrospective analyses of EEG data collected from all adult patients implanted with a VNS device for drug-resistant epilepsy in the Brno Epilepsy Center between 2005 and 2015. Data from patients implanted between 2005 and 2012 were used for investigation of EEG reactivity to external stimuli and subsequently for development of the statistical model (Cohort 1). Additional data from patients implanted with VNS between 2013 and 2015 were used as independent data-set for validation of the statistical model (Cohort 2). All the data were acquired during routine outpatient pre-operative assessment, 20 min recording at morning with two standard eyes opening and closing activation procedures (i.e., 10 s period with eyes open), photic stimulation (PS), and hyperventilation (HV). Each EEG recording was filtered into individual frequency bands and segmented into specific intervals representing the eyes opening and closing, PS and HV periods. The relative powers of EEG spectrum in distinct time intervals for a particular frequency band and brain area were then calculated.

Based on their individual responses to VNS, patients were categorized as Responders or Non-responders. In Cohort 1, Responders and Non-responders were first compared on the relative power values, and then a statistical model for prediction of VNS efficacy was developed. Subsequently, the validity of the statistical model was tested in Cohort 2. The study was conducted in St. Anne's University Hospital and approved by the local ethics committee. All patients gave their informed written consent for the use of their pre-operative data.

Patients' Description

All patients were implanted with a VNS system (Cyberonics, Houston) according to a standard implantation procedure (9). Before implantation, all patients underwent a comprehensive assessment protocol for epilepsy surgery candidates, including a detailed history and neurological examination, magnetic resonance imaging (MRI), interictal PET, neuropsychological testing, and scalp video-EEG monitoring. In some patients, ictal and interictal SPECT (SISCOM) and invasive video-EEG monitoring have been completed if necessary. Based on the results of all the investigations, patients indicated for VNS and included in this study were ruled out as suitable candidates for resective epilepsy surgery. Our analyses were applied to data acquired from patients who fulfilled the following criteria: The duration of VNS treatment was at least 2 years; the efficacy of VNS treatment was determined in regular visits every 3 or 6 months; and artifact-free pre-operative interictal EEG was available for eye opening and closing, PS and HV periods.

Demographic information and data regarding the type and number of antiepileptic drugs (AEDs) at the time of implantation were obtained by review of patients' charts. The efficacy of VNS was categorized using a classification system reported by McHugh (10). The cut-off value for seizure-reduction between responders and non-responders was 50%. Patients were defined as Responders or Non-responder only if they were categorized as such for the entire follow-up period.

EEG Analysis

First, we compared relative EEG powers between Responders and Non-responders in Cohort 1. Interictal scalp EEG was recorded on a 64-channel Alien Deymed system with international 10–20 electrode placement and a sampling frequency of 128 Hz. Standard antialiasing filters were applied before digitalization. Occasional artifacts were rejected manually and further processing was performed with artifact-free EEG periods. The resulting EEG signals were filtered into four frequency bands: theta (4–7.5 Hz), alpha (8–12 Hz), beta (14–30 Hz), and gamma (31–45 Hz). A Hilbert transform was then used to estimate the envelopes of pre-defined pass-band frequency oscillations as a function of time (Figure 1). The EEG records were segmented into the following conditions (time intervals):

1 –Rest#1 (2 min)

2 –Eyes opening/closing (10 s)

3 –Rest#2 (immediately after eye closure; 10 s)

4 –Photic stimulation (PS; 2.5 min)

5 –Hyperventilation (HV; 4 min)

6 –Eyes opening/closing (10 s)

7 –Rest#3 (immediately after eye closure; 10 s)

8 –Rest#4 (2 min)

FIGURE 1

Figure 1. Pre-processing of EEG signal. The EEG was segmented into eight time intervals, and then filtered into four frequency bands: theta, alpha, beta, and gamma. Subsequent analyses focused on the oscillatory power changes within these intervals, which were analyzed separately for each frequency band.

Further analysis was focused on oscillatory power changes in these conditions. Absolute mean power of the EEG spectrum was computed as a mean value of the passband power envelope inside each interval separately, for each scalp electrode. Subsequently, relative mean power (RPW) was calculated as a percentage decrease or increase of mean power relative to baseline, i.e., event-related desynchronization or synchronization, respectively. As a baseline we selected Rest#1 (11). We then evaluated differences between Responders and Non-responders by comparing relative power in seven conditions of clinical assessment: Open/Close#1, Rest#2, PS, HV, Open/Close#2, Rest#3, and Rest#4.

Statistical Analysis

Demographic data were compared between Responders and Non-responders using Fisher's exact test or the Mann-Whitney test. Statistical comparisons of RPW between Responders and Non-responders were also performed with Mann-Whitney tests. Using false discovery rate (FDR) (12), p-values for all electrodes were corrected for multiple comparisons in each time interval separately. Differences were considered significant when p ≤ 0.05.

Prediction of Response to VNS

Developing the Statistical Model

The statistical model for prediction of VNS efficacy was developed from data obtained in Cohort 1 in the following steps: Firstly, electrodes were grouped into seven anatomical regions as follows: (1) left frontal—Fp1, F3, Fz; (2) right frontal—Fp2, F4, Fz; (3) left anterotemporal—F7, T3; (4) right anterotemporal—F8, T4; (5) central—C3, Cz, C4; (6) left posterior quadrant—P3, Pz, T5, O1; (7) right posterior quadrant—P4, Pz, T6, O2. Specifically, mean values were calculated for the respective groups of electrodes. This resulted in 196 electrode group variables (7 conditions × 4 frequency bands × 7 anatomical regions). In the second step, the maximally discriminative electrode group variables were then selected using stepwise logistic regression performed in leave-one-out (LOO) manner; specifically, with n patients, logistic regression was performed over n iterations, each with one subject omitted. This approach allowed us to avoid overestimating the classification results, which occurs when classification is performed on electrode groups selected using all subjects simultaneously. Thirdly, electrode group variables identified most frequently as maximally discriminative after the n LOO iterations were used for classification using three classifiers; namely, logistic regression (LR), linear support vector machines (SVM), and linear discriminant analysis (LDA). In this step, we used a LOO cross-validation to split the data into training and testing sets: one patient was chosen randomly as a testing subject and the remaining patients were employed for training the classifier. The testing subject was then classified as belonging to the Responder or to the Non-responder class, and the resulting class label was compared with the true classification label. This procedure was repeated using each of the subjects as the testing subject sequentially, and the overall classification performance measures of accuracy, sensitivity, and specificity were calculated (see the black schematics in Figure 2). Fourth, using one-sample binomial test we performed a comparison of the achieved classification accuracies with a classification by chance. We also attempted to predict VNS efficacy based on 532 possible single electrode variables (7 conditions × 4 frequency bands × 19 electrodes; see Supplementary Material 1).

FIGURE 2

Figure 2. The statistical model for prediction of VNS efficacy based on EEG analysis. Schematics representing the development of the statistical model on data from Cohort 1 outlining the selection of the most discriminative electrodes groups as well as classifier training (CT) and testing (T) performed in LOO manner (the black schematics); and two approaches taken for the validation procedure, whereby our statistical model was applied to an independent data-set from Cohort 2 (the gray schematics). Equivalent steps taken in each of the two validation approaches are visualized using solid lines, including the reduction of Cohort 2 by selecting the maximally discriminative electrodes identified in Cohort 1. The two approaches differ in the steps represented by dashed lines—i.e., the classification (C) of the reduced dataset using: (A) a single classifier trained on all n subjects from Cohort 1, and (B) voting of n classifiers trained on n-1 subsets of Cohort 1.

Validation of the Statistical Model

The validity of the aforementioned statistical model was verified on the independent data-set obtained from Cohort 2. This validation cohort consisted of patients suffering from drug-resistant epilepsy who were implanted with VNS in our Center between 2013 and 2015. We obtained EEG data recorded using an identical EEG system and pre-surgical assessment protocol, and processed mathematically in the exact same way described above. Patients were classified as Responders and Non-responders with our statistical model using two approaches (see the gray schematics in Figure 2). Both began with data reduction, whereby we selected only data corresponding to the maximally discriminative electrode group variables identified based on Cohort 1. In the first approach, this data reduction was followed subsequently by a classification of m subjects from Cohort 2 as Responders or Non-responders based on a single classifier, which was trained on all n subjects from the Cohort 1. In contrast, in the second approach we classified the reduced data of Cohort 2 using majority voting of n classifiers trained on n-1 subsets of Cohort 1. The resulting indices for accuracy, sensitivity, and specificity were then compared to those from the classification of Cohort 1.

Results

Participants

Cohort 1

The VNS device was implanted in 110 patients in our center between 2005 and 2012. We excluded 50 patients (45%) from further analyses−17 (15%) because of poor-quality pre-surgical EEG, 13 (12%) due to attrition, and 20 (18%) who switched between VNS outcomes during the follow-up period (see Supplementary Material 2 for details). Demographics for the 60 patients included in the analyses are summarized in Table 1. According to the response to VNS, patients were subdivided into 35 (58%) Responders and 25 (42%) Non-responders. When examining demographic data, we observed significant differences between Responders and Non-responders in patients' age at epilepsy onset, duration of epilepsy before VNS implantation, and treatment by valproic acid (Table 1).

TABLE 1

Table 1. Demographic and treatment data for Cohort 1 and Cohort 2.

Cohort 2

In total, the VNS device was implanted in 25 patients between 2013 and 2015. Three patients were excluded from our analyses due to poor-quality pre-surgical EEG recordings. Data from the remaining 22 patients were used for independent validation: 12 (54.5%) Responders and 10 (45.5%) Non-responders. The demographic data are summarized in Table 1. The Cohort 2 did not differ from the Cohort 1 in the demographic characteristics, apart from the anticipated difference in duration of VNS (p < 0.001) and subtle differences in antiepileptic treatment (specifically, eslicarbazepine, lacosamide, and zonisamide were more frequent in Cohort 2; p = 0.029, p = 0.003, and p = 0.012, respectively).

RPW Differences Between Responders and Non-responders

In Cohort 1, computation of RPWs and their dynamics in the recording sites revealed significant changes for both Responders and Non-responders in all investigated frequency bands, all defined conditions, and the majority of electrodes (Figure 3). These changes were largely equivalent for both Responders and Non-responders within specific frequency bands and conditions. In the alpha and gamma frequency bands, however, there were striking dissimilarities between the two patient groups in conditions, with the most prominent differences during photic stimulation (PS) and hyperventilation (HV). Comparing Responders and Non-responders with respect to the RPW in all pre-defined frequency bands and conditions, we revealed significant differences in alpha and gamma frequency bands (Figure 4).

FIGURE 3

Figure 3. Relative mean powers (RPWs) for all frequency bands. Each head represents RPW in a given condition and frequency band. The RPW percentage scale is displayed on the right-hand side. Where a statistically significant difference exists in a given condition relative to baseline (Rest#1), the electrode is marked by a white dot. If no significant difference exists, the electrode is marked by a black dot.

FIGURE 4

Figure 4. Relative mean power (RPW) differences between Responders and Non-responders. Statistically significant or non-significant differences in a given electrode are marked by a white or black dot, respectively. The significant differences between Responders and Non-responders were identified in alpha and gamma frequency ranges only (heads in red boxes).

When analyzing the alpha band, significant differences between Responders and Non-responders were found in the following conditions: 4-PS, 5-HV, and 8–Rest#4. The differences were present over different brain regions in each condition: In 4-PS, there were differences over central and anterior areas (absence of desynchronization in Responders); in the 5-HV interval we observed differences over central and posterior regions (significantly higher synchronization in Responders); and in 8-R, differences were localized to right anterior and left posterior areas (higher power in Responders but persisting desynchronization in Non-responders).

When focusing on the gamma frequency band, we observed significant differences (gamma synchronization in Responders and desynchronization in Non-responders) in RPW within 4-PS, 5-HV and 7-Rest#3: In the first two conditions, differences were distributed across almost the whole scalp. The differences in RPWs within 7-Rest#3 were observed within right anterior and left posterior areas.

Prediction of VNS Response—Statistical Model

Based on the EEG data of Cohort 1, the statistical model was developed. Eight groups of electrodes were selected as the most discriminative in this statistical model (visualized in Figure 5). The best classification results based on these eight most discriminative groups of electrodes were obtained using the LR classifier, achieving 86.7% accuracy, with 88.6% sensitivity and 84% specificity (Table 2). This classification accuracy was significantly higher than those achieved by chance (p < 0.001). The SVM classifier achieved lower accuracy (75%) but it was still significantly higher than that achieved by chance (p = 0.004). The lowest classification performance measures were obtained in classification using LDA (accuracy 65%). Detailed visualizations of classification accuracy are provided in Figure 6. The final statistical model is described in more detail in Supplementary Material 3.

FIGURE 5

Figure 5. Electrode groups and individual electrodes selected for statistical modeling. Electrode groups (A) and individual electrodes (B) selected for the development of the statistical model, in each condition and frequency band.

TABLE 2

Table 2. Classification performance indices for electrode group variables.

FIGURE 6

Figure 6. Classification accuracy based on a number of selected groups of electrodes. The bottom row of numbers shows a minimum number of LOO iterations (# iterations) in which a given number of electrode groups (# groups) were selected as most discriminative (e.g., eight electrode groups were selected as maximally discriminative in 35 or more LOO iterations). The best LR classification accuracy and corresponding results for SMV and LDA (shown in Table 2) are depicted by circles and a dashed line.

Validation of the Statistical Model

The data of patients in Cohort 2 were used for independent validation. When comparing the results of our classification model applied to Cohort 2 against real-life outcome, the best results were again obtained for the model using the LR classifier achieving an accuracy of 86.4%, sensitivity of 83.3%, and specificity of 90%. The complete results achieved using classifier voting are summarized in Table 2. Classification performance based on a single classifier was lower (data not shown).

Discussion

Over 80,000 epilepsy patients have been treated worldwide using vagal nerve stimulation (VNS), a standard of modern non-pharmacological treatment. Importantly, the number of epileptic patients who are eligible for VNS therapy is approximately three million subjects worldwide. In this light, a method for the reliable prediction of individual efficacy of VNS therapy is desperately needed.

Recently, several authors have attempted to identify predictors of VNS outcome (1, 2, 4, 13). Despite these efforts, however, clinical predictors of individual responsiveness to VNS therapy remain elusive. Similar estimates of efficacy are reported for diverse neurostimulation techniques in the treatment of drug-resistant epilepsy (i.e., VNS, Deep Brain Stimulation of Anterior Thalamic Nuclei, Brain Responsive RNS System, transcutaneous VNS, or external Trigeminal Nerve Stimulation), which might suggest a more consequential impact of external stimuli per se on epileptic activity. For this reason, we retrospectively evaluated routine EEG data acquired before implantation in a large cohort of VNS patients. Using standard computations of power spectral analyses of interictal EEG, we reveal significant differences between responders and non-responders in two pre-defined frequency bands (alpha and gamma) and four conditions of standard clinical assessment. Based on RPWs and their dynamics, we have developed and validated a statistical model for prediction of VNS efficacy that discriminated between responders and non-responders with almost 90% accuracy.

Our primary finding is that VNS responders and non-responders differ significantly in EEG power dynamics within alpha and gamma frequency bands prior to therapy. Whilst both patient groups demonstrated equivalent alpha desynchronization during eyes opening, they differed in alpha reactivity to photic stimulation and hyperventilation; specifically, responders showed no decrease in alpha power during the former but an enormous increase during the latter. This reactivity pattern stands in contrast to that observed in healthy individuals, in whom photic stimulation typically leads to alpha attenuation and standardized hyperventilation has been shown to decrease alpha power (14). Interestingly, we also observed significant increases in gamma power during both photic stimulation and hyperventilation in responders relative to non-responders. Hyperventilation-induced physiological changes are thought to be a consequence of increased neuronal excitability resulting from the hypocapnia-induced alkalosis (15, 16). Significantly enhanced alpha and gamma activities during hyperventilation in responders might reflect distinct properties of responders' brains vis-a-vis neuronal excitability and synaptic transmission. Nevertheless, hyperventilation is a complicated practice which, in epileptic subjects, results in unpredictable responses. The reasons for this remain unclear (e.g., metabolic hypersynchronization, altered neurotransmitters, etc.). It also remains unclear why the alpha (and less expressed gamma) desynchronization persist in our non-responders at the end of EEG recordings (Rest#4). It seems that the protocol itself, particularly photic stimulation and hyperventilation, induces some change from the resting-state baseline for non-responders. Differential dynamics of power changes after the stimulation—faster in responders and slower in non-responders—might represent another characteristic in which these groups of patients differ substantially.

Interestingly, valproic acid (VPA) was used more frequently in pharmacological treatment prior to implantation in our VNS responders compared to the non-responders. We might therefore ask whether differences in EEG reactivity observed in our study is related to some pharmacological imprint. Although this speculative explanation cannot be fully excluded, it seems unlikely that pharmacological impact on EEG is the main factor driving our findings. Still more unlikely is substantial VPA impact on the prediction of VNS efficacy, bearing in mind that only one third of responders were treated with VPA and the accuracy of individual VNS efficacy prediction was almost 90%.

Since Hans Berger's initial observation in 1933, the best known example of EEG reactivity is alpha attenuation (alpha blocking and desynchronization). This is observed typically when subjects open their eyes (the Berger effect), but alpha also disappears when subjects become drowsy and it can be blocked by numerous kinds of external stimuli (visual or auditory) or mental operations (e.g., imagery, visualization, mental arithmetics). This implies that alpha reactivity as a modality-independent, general phenomenon, reflecting the functional modes of thalamo-cortical and cortico-cortical loops that facilitate/inhibit the transmision of information in the brain (17, 18). Particularly noteworthy is the high inter-individual variability in alpha reactivity; it has been shown to differ between extraverts and introverts, for example, and there is less pronounced alpha desynchronization in people with high intelligent quotient during several cognitive tasks (19, 20). Alpha reactivity has also been reported to decrease with aging (21–24) and in patients with mild cognitive impairment and Alzheimer disease (25), and a lack of alpha reactivity has been used to predict the long-term deterioration of higher functions in subjects with cognitive decline (26). From this viewpoint, our discovery of differential alpha reactivity in VNS responders vs. non-responders further emphasizes the multifold functions of the diffuse alpha system (27).

In our study, significant differences in EEG power dynamics within the gamma frequency band between VNS responders and non-responders very likely reflects distinct reactivity of true brain gamma. When interpreting thoroughly our results in this frequency band, however we shall keep in mind the study of Whitham et al. (28); these authors showed recently that even the normal resting EEG might reveal significant contamination with electromyography (EMG) activity in this frequency band. Further, the level of EMG contamination increased dramatically when subjects perform various experimental tasks and are sitting during EEG recordings (29). More recently, however, Boytsova et al. showed that EMG contamination does not necessarily hide high-frequency EEG and does not preclude qualitative detections of electroencephalographic correlates of mental activities in beta and low gamma frequency ranges (30). Unfortunately, despite the availability of several techniques for the reduction of muscular artifacts from EEG traces, at a present none are able to guarantee that analyzed data are completely free of high-frequency artifacts (29). In our study, EEG was recorded under standard conditions: patients laid comfortably and were instructed to relax (including their facial muscles), no task was performed, and recordings from patients containing artifacts identified with careful visual inspection (e.g., muscular activity) were excluded (17 out of 110 patients). As such, we strongly believe an increased gamma power during both photic stimulation and hyperventilation in responders relative to non-responders is not resulting from distinct muscle artifacts contaminations.

Both brain rhythms—alpha and gamma—are considered to represent a kind of universal code consistent with their putative role in brain signaling (27). Both are generated in a widely distributed system, with a major role of thalamocortical circuits in their origin (31–33). The differential impact of external stimuli on alpha and gamma between VNS responders and non-responders might be mediated by differences in neuronal interconnectivity and different levels of neurotransmitters within underlying cerebral matrices (24). Such differences might influence the effect of external stimuli delivered to thalamocortical circuits and other brain networks via the vagal nerve (34, 35). Consistent with this notion, many consider the mechanism of VNS action to be modulation of synaptic activity in the thalamus and thalamocortical projections, increased plasticity in GABA receptors, and modulation of GABAergic activity that is related directly to gamma oscillations (6, 36). Indeed, one of the most plausible factors involved in gamma variation between our populations of responders vs. non-responders is represented by the role of GABA neurotransmission and divergent dynamics of inhibitory interneuron networks within the central nervous system (37, 38). The close relationship between different types of external stimulation and brain reactivity we have observed might be indicative of a common mechanism underlying the various forms of neurostimulation in epilepsy treatment.

Finally, as seen in the developed and validated statistical model, for accurate prediction of VNS efficacy distinct recording electrode groups, different conditions (especially hyperventilation, but also eyes opening/closing and resting periods), and all frequency bands were selected as the most discriminative ones. Using our approach we achieved high accuracy, sensitivity and specificity in both well-defined and predominant VNS patients. We also evaluated a statistical model based on single electrodes, which had slightly superior classification performance in Cohort 1 and predominant VNS patients but failed in the independent validation set of Cohort 2. This indicates that groups of electrodes are better suited for VNS efficacy prediction—by integrating data from larger regions, such groupings appear to produce more robust results.

To conclude, we have revealed that EEG reactivity to external stimuli used during routine pre-operative EEG investigation differs between VNS responders and non-responders. Moreover, we have developed and validated a statistical tool that can predict with extremely high accuracy whether or not individual drug-resistant epileptic patients will benefit from VNS treatment (Patent Number EP3437692-A1). This electrophysiological marker could prove invaluable when providing patients with expected postoperative prognosis. Further research is required before this can be achieved, however; our findings come from a retrospective and monocenter study, with a limited number of patients comprising the independent validation. Our results must therefore be replicated in prospective, multicentre, and well-designed clinical study in order to obtain a clear-cut statistical power.

Ethics Statement

Ethic committee in St. Anne's University hospital approved this study. All patients gave their informed consent.

Author Contributions

MB, RR, PJ, JChr, and DS: identification of research topic, preparation of study design, interpretation of results. ID and MP: preparation of manuscript, design of study, interpretation of results. JChl and PJ: mathematical analysis. EK: statistical analysis.

Funding

The results of this research have been acquired within CEITEC 2020 (LQ1601) project with financial contribution made by the Ministry of Education, Youth and Sports of the Czech Republic (MEYS CR) within special support paid from the National Programme for Sustainability II funds, MEYS CR project LO1212, and the Ministry of Health of the Czech Republic, grant nr. NV19-04-00343.

Conflict of Interest Statement

The patent application was applied (patent applicant Masaryk University; inventors are the authors of the manuscript: MB, ID, EK, JChl, RR, JChr, MP).

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to thank Matthew C. Walker for his valuable help and advice with manuscript.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2019.00392/full#supplementary-material

References

1. Englot DJ, Chang EF, Auguste KI. Vagus nerve stimulation for epilepsy: a meta-analysis of efficacy and predictors of response. J Neurosurg. (2011) 115:1248–55. doi: 10.3171/2011.7.JNS11977

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Englot DJ, Rolston JD, Wright CW, Hassnain KH, Chang EF. Rates and predictors of seizure freedom with vagus nerve stimulation for intractable epilepsy. Neurosurgery. (2016) 79:345–53. doi: 10.1227/NEU.0000000000001165

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Theodore WH, Fisher RS. Brain stimulation for epilepsy. Lancet Neurol. (2004) 3:111–8. doi: 10.1016/S1474-4422(03)00664-1

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Alexander GM, McNamara JO. Vagus nerve stimulation elevates seizure threshold in the kindling model. Epilepsia. (2012) 53:2043–52. doi: 10.1111/j.1528-1167.2012.03646.x

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Jaseja H. EEG-desynchronization as the major mechanism of anti-epileptic action of vagal nerve stimulation in patients with intractable seizures: clinical neurophysiological evidence. Med Hypotheses. (2010) 74:855–6. doi: 10.1016/j.mehy.2009.11.031

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Marrosu F, Santoni F, Puligheddu M, Barberini L, Maleci A, Ennas F, et al. Increase in 20-50 Hz (gamma frequencies) power spectrum and synchronization after chronic vagal nerve stimulation. Clin Neurophysiol. (2005) 116:2026–36. doi: 10.1016/j.clinph.2005.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Fraschini M, Puligheddu M, Demuru M, Polizzi L, Maleci A, Tamburini G, et al. VNS induced desynchronization in gamma bands correlates with positive clinical outcome in temporal lobe pharmacoresistant epilepsy. Neurosci Lett. (2013) 536:14–8. doi: 10.1016/j.neulet.2012.12.044

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Bodin C, Aubert S, Daquin G, Carron R, Scavarda D, McGonigal A, et al. Responders to vagus nerve stimulation (VNS) in refractory epilepsy have reduced interictal cortical synchronicity on scalp EEG. Epilepsy Res. (2015) 113:98–103. doi: 10.1016/j.eplepsyres.2015.03.018

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Landre E. Vagus nerve stimulation and refractory partial epilepsies. Rev Neurol. (2004) 160:S280–7.

PubMed Abstract | Google Scholar

10. McHugh JC, Singh HW, Phillips J, Murphy K, Doherty CP, Delanty N. Outcome measurement after vagal nerve stimulation therapy: proposal of a new classification. Epilepsia. (2007) 48:375–8. doi: 10.1111/j.1528-1167.2006.00931.x

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Pfurtscheller G, Aranibar A. Event-related cortical desynchronization detected by power measurements of scalp EEG. Electroencephalogr Clin Neurophysiol. (1977) 42:817–26.

PubMed Abstract | Google Scholar

12. Benjamini Y, Hochberg Y. Controlling the false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc B. (1995) 57:289–300.

Google Scholar

13. Liu H, Yang Z, Huang L, Qu W, Hao H, Li L. Heart-rate variability indices as predictors of the response to vagus nerve stimulation in patients with drug-resistant epilepsy. Epilepsia. (2017) 58:1015–22. doi: 10.1111/epi.13738

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Van der Worp HB, Kraaier V, Wieneke GH, Van Huffelen AC. Quantitative EEG during progressive hypocarbia and hypoxia. Hyperventilation-induced EEG changes reconsidered. Electroencephalogr Clin Neurophysiol. (1991) 79:335–41.

PubMed Abstract | Google Scholar

15. Wyke BD. Brain function and blood sugar: observations based on a case of islet cell adenoma of the pancreas. Electroencephalogr Clin Neurophysiol. (1952) 4:339–50.

PubMed Abstract | Google Scholar

16. Neurbauer AC, Sange G, Pfurtscheller G. Psychometric intelligence and event-related desynchronization during performance of a letter matching task. In: Pfurtscheller G, Lopes da Silva FH, editors. Handbook of Electroencephalography nad Clinical Neurophysiology. Amsterdam: Willy (1999). p. 219–31.

Google Scholar

17. Steriade M, Llinas RR. The functional states of the thalamus and the associated neuronal interplay. Physiol Rev. (1988) 68:649–742.

PubMed Abstract | Google Scholar

18. Pfurtscheller G, Lopes da Silva FH. Event-related EEG/MEG synchronization and desynchronization: basic principles. Clin Neurophysiol. (1999) 110:1842–57.

PubMed Abstract | Google Scholar

19. Venturini R, De Pascalis V, Imperiali MG, Martini PS. EEG alpha reactivity and extraversion-introversion. Pers Indiv Differ. (1981) 2:215–20.

Google Scholar

20. Neubauer A, Freudenthaler HH, Pfurtscheller G. Intelligence and spatiotemporal patterns of event-related desynchronization (ERD). Intelligence. (1995) 20:249–66.

Google Scholar

21. Duffy FH, Albert MS, McAnulty G. Brain electrical activity in patients with presenile and senile dementia of the Alzheimer type. Ann Neurol. (1984) 16:439–48.

PubMed Abstract | Google Scholar

22. Kononen M, Partanen JV. Blocking of EEG alpha activity during visual performance in healthy adults. A quantitative study. Electroencephalogr Clin Neurophysiol. (1993) 87:164–6.

PubMed Abstract | Google Scholar

23. Marciani MG, Gotman J. Effect of drug-withdrawal on location of seizure onset. Epilepsia. (1986) 27:423–31.

Google Scholar

24. Gaal ZA, Boha R, Stam CJ, Molnar M. Age-dependent features of EEG-reactivity–spectral, complexity, and network characteristics. Neurosci Lett. (2010) 479:79–84. doi: 10.1016/j.neulet.2010.05.037

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Babiloni C, Lizio R, Vecchio F, Frisoni GB, Pievani M, Geroldi C, et al. Reactivity of cortical alpha rhythms to eye opening in mild cognitive impairment and Alzheimer's disease: an EEG study. J Alzheimers Dis. (2010) 22:1047–64. doi: 10.3233/JAD-2010-100798

PubMed Abstract | CrossRef Full Text | Google Scholar

26. van der Hiele K, Bollen EL, Vein AA, Reijntjes RH, Westendorp RG, van Buchem MA, et al. EEG markers of future cognitive performance in the elderly. J Clin Neurophysiol. (2008) 25:83–9. doi: 10.1097/WNP.0b013e31816a5b25

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Basar E, Schurmann M, Basar-Eroglu C, Karakas S. Alpha oscillations in brain functioning: an integrative theory. Int J Psychophysiol. (1997) 26:5–29.

PubMed Abstract | Google Scholar

28. Whithain EM, Pope KJ, Fitzgibbon SP, Lewis T, Clark CR, Loveless S, et al. Scalp electrical recording during paralysis: Quantitative evidence that EEG frequencies above 20 Hz are contaminated by EMG. Clin Neurophysiol. (2007) 118:1877–88. doi: 10.1016/j.clinph.2007.04.027

CrossRef Full Text | Google Scholar

29. Muthukumaraswamy SD. High-frequency brain activity and muscle artifacts in MEG/EEG: a review and recommendations. Front Hum Neurosci. (2013) 7:138. doi: 10.3389/fnhum.2013.00138

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Boytsova JA, Danko SG, Medvedev SV. When EMG contamination does not necessarily hide high-frequency EEG: scalp electrical recordings before and after Dysport injections. Exp Brain Res. (2016) 234:3091–106. doi: 10.1007/s00221-016-4708-3

CrossRef Full Text | Google Scholar

31. Andersen P, Andersson SA, Lomo T. Thalamo-cortical relations during spontaneous barbiturate spindles. Electroencephalogr Clin Neurophysiol. (1968) 24:90.

PubMed Abstract | Google Scholar

32. Lopes da Silva FH, van Lierop THMT, Schrijer CF, Storm van Leeuwen W. Organization of thalamic and cortical alpha rhythms: Spectra and coherences. Electroencephalogr Clin Neurophysiol. (1973) 35:627–39.

Google Scholar

33. Timofeev I, Steriade M. Fast (mainly 30-100 Hz) oscillations in the cat cerebellothalamic pathway and their synchronization with cortical potentials. J Physiol. (1997) 504:153–68.

PubMed Abstract | Google Scholar

34. Bartolomei F, Bonini F, Vidal E, Trebuchon A, Lagarde S, Lambert I, et al. How does vagal nerve stimulation (VNS) change EEG brain functional connectivity? Epilepsy Res. (2016) 126:141–6. doi: 10.1016/j.eplepsyres.2016.06.008

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Wostyn S, Staljanssens W, De Taeye L, Strobbe G, Gadeyne S, Van Roost D, et al. EEG derived brain activity reflects treatment response from vagus nerve stimulation in patients with epilepsy. Int J Neural Syst. (2017) 27:1650048. doi: 10.1142/S0129065716500489

PubMed Abstract | CrossRef Full Text | Google Scholar

36. De Herdt V, De Waele J, Raedt R, Wyckhuys T, El Tahry R, Vonck K, et al. Modulation of seizure threshold by vagus nerve stimulation in an animal model for motor seizures. Acta Neurol Scand. (2010) 121:271–6. doi: 10.1111/j.1600-0404.2009.01223.x

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Whittington MA, Traub RD, Kopell N, Ermentrout B, Buhl EH. Inhibition-based rhythms: experimental and mathematical observations on network dynamics. Int J Psychophysiol. (2000) 38:315–36. doi: 10.1016/S0167-8760(00)00173-2

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Traub RD, Whittington MA, Colling SB, Buzsáki G, Jefferys JG. Analysis of gamma rhythms in the rat hippocampus in vitro and in vivo. J Physiol. (1996) 493(Pt 2):471–84.

PubMed Abstract | Google Scholar