Abstract

Objective. Studies have demonstrated the utility of closed-loop neuromodulation in treating focal onset seizures. There is an utmost need of neurostimulation therapy for generalized tonic-clonic seizures. The study goals are to map the thalamocortical network dynamics during the generalized convulsive seizures and identify targets for reliable seizure detection. Methods. Local field potentials were recorded from bilateral cortex, hippocampi, and centromedian thalami in Sprague-Dawley rats. Pentylenetetrazol was used to induce multiple convulsive seizures. The performances of two automated seizure detection methods (line length and P-operators) as a function of different cortical and subcortical structures were estimated. Multiple linear correlations-Granger’s Causality was used to determine the effective connectivity. Results. Of the 29 generalized tonic-clonic seizures analyzed, line length detected 100% of seizures in all the channels while the P-operator detected only 35% of seizures. The detection latencies were shortest in the thalamus in comparison to the cortex. There was a decrease in amplitude correlation within the thalamocortical network during the seizure, and flow of information was decreased from thalamus to hippocampal-parietal nodes. Significance. The preclinical study confirms thalamus as a superior target for automated detection of generalized seizures and modulation of synchrony to increase coupling may be a strategy to abate seizures.

1. Introduction

Among the different seizure subtypes, generalized tonic-clonic seizures are most disabling due to loss of consciousness, the potential to cause physical injuries, and cardiorespiratory compromise including SUDEP (Sudden Unexpected Death in Epilepsy) [1, 2]. Generalized tonic-clonic seizures can occur in primary generalized epilepsy (now termed as Genetic Generalized Epilepsy) or in partial epilepsy where the seizures are secondarily generalized. Despite therapeutic advances in the treatment of epilepsy, the outcome of primary generalized epilepsy remains suboptimal. When antiseizure medication fails to control seizures, in partial epilepsy resective surgery can potentially cure, but the therapeutic options in generalized epilepsy are limited. Recently an FDA approved closed-loop responsive neurostimulation therapy (NeuroPace) was proven to be beneficial in interrupting seizures including convulsive seizures of partial onset [3]. We wondered if such a treatment paradigm can be extended to control primarily generalized seizures.

Experimental models of generalized epilepsy and functional imaging studies in human suggest that generalized seizures are initiated in cortical nodes, mainly frontal and parietal, before they entrain the thalamus in a bidirectional ictal network [47]. Maintenance of the integrity of the reciprocal thalamocortical reverberatory loop is necessary for sustaining ictogenesis. Therefore to design a physiology-based, rational closed-loop stimulation paradigm would require an understanding of the thalamocortical network dynamics during seizure evolution and identifying targets within the thalamocortical loop for reliable seizure detection. To address these critical knowledge gaps, we have performed a preclinical study with the focus on identifying detection metrics that can trigger feedback stimulation effectively to abate generalized convulsive seizures. Specifically, the study aims are to (a) evaluate the performance of automated seizure detectors (line length and P-operator) as a function of different corticosubcortical targets and (b) map the peri-ictal information flow dynamics involving the corticosubcortical network in a validated rodent model of generalized convulsive seizures [8, 9]. The correlation or synchronization of neural activity within seizure-generating sites is implicated in seizure genesis, and desynchronization is a potential mechanism proposed in the stimulation induced suppression of seizures [1012]. Abnormal synchrony in the corticosubcortical network has also been reported in generalized spike-wave seizures [13, 14]. Therefore, based on these studies, we speculate synchrony as a potential target for rational brain stimulation paradigm. The significance and role of synchrony are likely to depend on the extent of the coupling of interconnected neurons. Coupling measures are used to establish relationship between two EEG signals from different channels of single brain region or two brain regions. Among different methods to analyze coupling, Granger’s Causality is one popular method to analyze the coupling strength and direction or information flow of EEG signals between two electrodes in different brain regions. Therefore, using stepwise multiple linear correlations-Granger’s Causality analysis, here we explore the spatiotemporal coupling within and between the corticosubcortical networks during seizure evolution. We hypothesize that the causal sources of generalized seizures are spatially distributed within the thalamocortical network and are abnormally correlated in the time domain during seizure evolution. Furthermore, as ictogenesis of generalized seizure involves rapid recruitment of widespread corticosubcortical networks [1], we anticipate that the performance of automated seizure detectors as a function of anatomical targets will be invariant.

Pentylenetetrazol (PTZ), a GABA-A receptor antagonist, is the most commonly used proconvulsant to induce spike and wave, myoclonus, and generalized tonic-clonic seizures [9]. The tonic-clonic seizures induced by PTZ are believed to represent generalized seizures. Unlike genetic models of generalized epilepsy where spike and wave absence seizures are predominant [15, 16], in the PTZ model after administration of the chemical, multiple generalized convulsive seizures, often in clusters, are induced. Since the goal of our study is to define stimulation metrics against generalized tonic-clonic seizures, we felt that the PTZ model suits better than genetic models for our study aims.

2. Methods

2.1. Animals

Experiments were performed on male Sprague-Dawley rats weighing 400–550 g that were obtained from Charles River Laboratories, MA, USA. Animals were housed in pairs with food and water ad libitum and kept in a 12 h light/dark cycle. All animal experimentations were approved by the Institutional Animal Care and Use Committee. All efforts were made to minimize the animal suffering and reduce the number of animals used in the experiments.

2.2. Surgery for Recording Electrodes

All procedures are performed using sterile techniques. Rats were anesthetized by inhalation of 5% isoflurane in an induction chamber. The periosteum was cleared from the cranium, and the exposed skull was cleaned with hydrogen peroxide (3%). Once the skull surface was dry, using bregma as the reference location, eight 0.25 mm craniotomies were performed with a stereotaxic drill for placement of monopolar depth electrodes that were custom-made from 175 μm tungsten microwire (California Fine Wire). Four cortical depth electrodes were placed targeting orbitofrontal and somatosensory cortex: AP: 2.0 mm, lateral: ±3.0, and depth: 1 mm; AP: −4.0 mm, lateral: ±3.0, and depth: 1.0 mm. Two hippocampal depth electrodes targeting dentate gyrus (coordinates: AP: −5.6 mm, lateral: ±4.5, and depth: 5.0 mm) and two additional depth electrodes targeting centromedial thalamic nuclei (coordinates: AP: −2.5 mm, lateral: ±1.5, and depth: 5.9 mm at a 10° lateral angle) were placed (Figure 1(b)) [17, 18]. A D-sub connector (pinsout) was then soldered to the recording electrode wires. The connector’s ground wire was exposed to the muscle tissue down the back of the neck. The entire head piece was then adhered to the skull using dental acrylic. A day following surgery, the rat was placed in a customized cage, and the head piece was connected to the video-EEG system (Natus Quantum EEG system™).

2.3. Seizure Induction and Electrophysiological Recording

Following a 7-day recovery period, seizures were induced by intraperitoneal injection of a single dose of PTZ (130 mg/kg) [19]. Seizures were induced within 20 minutes after administration of the proconvulsant. The recording was continued for 5 hours that allowed recording of multiple seizures. The video EEG was sampled at 1024 Hz.

2.4. Defining Onset of Generalized Seizures

Seizures were defined as abnormal electrographic activities lasting more than 10 seconds and associated with relatively high frequency and amplitude [19]. The local field potential (LFP) changes that were associated with generalized tonic-clonic seizures were continuous high frequency, high amplitude spike waves that were present bilaterally on all recorded channels. At times, these ictal electrographic changes transitioned to low-amplitude high-frequency activity before evolving into rhythmic high amplitude spikes during the clonic phase of the seizures. Seizure onset was defined as the earliest occurrence of high frequency, high amplitude spike waves throughout the recording channels that were distinctive from the background activity and that evolved in frequency and morphology. These changes were visually identified and marked as unequivocal electrographic onset (UEO) (Figure 1(a)). The seizure intensity was scored visually using a revised Racine scale that was validated in a previous study [20]. Only seizures that had clear clinical correlates of generalized tonic-clonic activity (stages 5 and 6) were included for analysis.

2.5. Automated Seizure Detectors
2.5.1. Line Length

The line length (LL) feature was derived as a simplification of the running fractal dimension of a signal [21]. It measures the length of the signal in a particular window and compares it to a variable threshold. The length of the signal is proportional to the amplitude and frequency of the signal, making this feature highly suitable to sense changes in amplitude and/or frequency that typically occur during seizures.

2.5.2. P-Operator

Clinical analysis of EEG is done by visual inspection. The visually appealing geometric features are therefore of primary importance in EEG signals. If we visualize a time domain EEG signal as the trajectory of a particle moving in a force field with one degree of freedom, then we will be able to trace the evolving geometry of the signal as the motion of the particle with variable acceleration along the ordinate. If is the time domain single channel EEG, then the ordinate at a time is and acceleration at is . The work done to displace the particle by an amount along the ordinate is . Time taken to accomplish this work is . So the rate at which the work was done is . In classical mechanics, is known as the power of the moving particle at the point . This is the power at which the particle is giving the specific shape to its trajectory. In other words, is the power at which the specific waveform of the signal at the point is being created. Note that this has nothing to do with the spectral power of . Since during an epileptic seizure distinct waveform changes take place in the EEG signal, of the background EEG and of the EEG during a seizure should be quite different. This way, enhances the contrast between the background EEG and the EEG during the seizure leading to convenient seizure detection by setting a suitable threshold or otherwise. The changes are no surprise because the efficacy of first- and second-order temporal difference in seizure detection has already been well documented [22, 23]. We call the power-operator or simply the P-operator. Custom written codes in Matlab (Mathworks, Natick, MA) were used for automated seizure detection.

2.6. Automated Detection of Seizures Offline and ROC Curves

A sweep of different thresholds was used to calculate the receiver operating characteristic (ROC) curves. We used the automated detection algorithms to determine true positives, true negatives, false positives, and false negatives for seizure detection. Visual markings of UEO were considered gold standard (Figure 1(a)). We calculated Sensitivity = TP/(TP + FN), Specificity = TN/(TN + FP), and false positive rate for each threshold to build an ROC curve and determine the optimal threshold for each animal and each channel (cortical, hippocampal, and thalamus).

2.7. Stepwise Multiple Linear Correlations-Granger’s Causality to Determine Effective Connectivity (EC) in the Periseizure Period

The seizure data was visually identified to determine segments preceding the seizure (labeled as PreSz, 1 min), during the seizure (Sz, the length of seizure), and immediately after the seizure (PostSz, 1 min) for analysis of effective connectivity between the channels. Using analysis of time-lagged relationships, effective connectivity is defined as the influence of one neural system exerting over another during an experimental context and therefore moves beyond describing instantaneous connections between brain regions and helps to clarify how brain areas communicate [24]. The eight channels that were used to estimate effective connectivity are placed bilaterally and spatially similar in frontal, centromedian thalami, parietal, and hippocampi. The data (sampled at 1024 Hz) was notch-filtered at 60 Hz and band-pass-filtered using finite impulse response filter between 0.01 and 100 Hz. Initially, multiple linear Pearson’s correlations were performed on spatially similar channels using a moving window Pearson Correlation (Ws = 1024 sample points, shifted by 1 sample point) with the pair of frontal, thalami, parietal, and hippocampi. This derived channel was obtained to generate a time series for each of the regions that contained the similar seizure ictal electrical activity between them and not the influence of background activity of the individual regions. Subsequently, the derived frontal, thalamus, parietal, and hippocampal channels were used to analyze effective connectivity (EC) using Granger’s Causality (GC), a form of multivariate vector autoregression model determining directed interregional coupling of a collection of time series, measured by one’s dependence over the other [25, 26]. GC was calculated across seizure and pre- and postseizure stages using the following parameters: window size of 1024 samples, 50% window overlap, AR order of 10, and number of surrogates of 100 to determine the statistically significant EC (p < 0.05, network threshold of 0.8) corresponding to using the type false discovery rate implementation [27]. The analysis resulted uniformly in 118 temporal windows in both the preseizure and the postseizure periods while the seizure window had variable windows depending on the length of the seizure. A schematic representation of the different pipeline steps is summarized in Figure 2.

2.8. Statistical Analysis of the Effective Connectivity Results

Once the GC was calculated for the three stages and across the channels, we tried to establish if the directional connectivity information was dependent on the channels involved or whether it was dependent on the stage of the seizure. Classification of the connectivity across stages and derivations was performed using a data-driven approach combining principle component analysis (PCA) using Varimax rotation and Kaiser Normalization [28, 29]. Finally, the difference in the EC between the three stages was also tested using repeated measures ANOVA with post hoc analysis. A corrected value less than 0.05 was considered significant.

3. Results

A total of 29 generalized tonic-clonic seizures (stage 5 and 6 seizure intensity) were analyzed out of 42 seizures recorded from two rats. Thirteen seizures were excluded due to seizure intensity below 5. The duration of the seizure ranged between 14 and 90 seconds. Our first strategy was to compare automated seizure detection against visual detection and plot ROC curves.

3.1. Automated Detection of Generalized Convulsive Seizures

The line length algorithm detected seizures in all the channels for all the 29 seizures (Figure 1(a)). The mean latency for detection was as follows: frontal cortex, seconds; CM thalamus, seconds; and hippocampus, seconds. The performance of the P-operator was suboptimal as it detected only 35% of seizures (10 each in the thalamus, cortex, and hippocampus) with the latencies varying between them (cortex: seconds; thalamus: seconds; and hippocampus: seconds) (Figure 3). ANOVA was performed to test the difference in seizure detection latency for the two automated detectors (Figure 4). The detection latencies varied as a function of anatomical targets with thalamus superior to frontal cortex for both line length and P-operator.

3.2. Stepwise Multiple Linear Correlations-Granger’s Causality to Determine Effective Connectivity (EC) in the Periseizure Period

There was a decrease in the linear amplitude correlation among the spatially identical channels (frontal, thalamic, parietal, and hippocampal derivations) during the seizure period. As the baseline (defined in this study as PreSz state) transitioned to seizure, there was a reduction in the coupling strength between all the corticosubcortical channels (Figure 5(a)). In the postseizure state, the coupling between the channels increased but this was not similar to the PreSz state ( for all the derivations) (Figure 5(b)). A data-driven approach classifying the changes in coupling and different seizure states (preSz-Sz-PostSz) showed that the net Granger’s Causality could be discriminated into three groups based on the seizures states (Table 1). Repeated measure ANOVA with post hoc analysis was done to evaluate the effective connectivity that had the greatest difference in mean GC between the three stages. The decrease in the flow of information during the seizure was directed from thalamus to hippocampus and from thalamus to parietal nodes. There was no significant difference in the information flow within frontothalamic reciprocal connectivity.

4. Discussion

In this study involving mapping of ictal network dynamics in a preclinical model of generalized tonic-clonic seizure, we highlight two pertinent findings that are of translational importance while designing a closed-loop neuromodulation. First, the performance of automated seizure detectors for detecting generalized seizures was variable with subcortical structures (thalami and hippocampi) superior to the cortical structures (frontoparietal). Second, during the generalized tonic-clonic seizure there was a decrease in coupling within the thalamocortical network, and the decrease in information flow was maximum from thalamus to the hippocampus and parietal network.

4.1. Subcortical Structures Are Superior Anatomic Targets for Automated Seizure Detection of Generalized Seizures

Performance of an automated detector of seizure onset may be influenced by variability in electrographic signatures of seizures, states of vigilance, site (like scalp versus intracranial EEG), and duration of recording to estimate performance [3032]. The onset of generalized tonic-clonic seizure involves widespread network incorporating corticosubcortical nodes, and the LFP is stereotyped by high frequency, high amplitude spike waves. Therefore, the results of the ROC curve demonstrating variability in performance between corticosubcortical structures were unanticipated. Our study demonstrated that, even for generalized seizure, the detection latency as a function of anatomical targets is not uniform and that thalamus or hippocampus is a better target for automated seizure detection. Although the differences in seizure detection latency were statistically nonsignificant, from the clinical standpoint a 4-second difference in detection latency can be significant if seizure suppression by intervention is the primary goal of the automated detector. The origin of generalized seizure is debated between cortical focus, thalamus, and thalamocortical network [4, 3335]. A recent study demonstrated bidirectional control of generalized epilepsy network by switching thalamocortical phasic firing to the tonic state using optogenetic modulation of the thalamus [36]. Our finding of thalamus being superior in automated seizure detection adds to the growing evidence that thalamus may be an attractive target for closed-loop neuromodulation in generalized epilepsy.

4.2. Disrupted Connectivity within the Thalamocortical Network during Generalized Tonic-Clonic Seizure

Synchronization (i.e., events occurring at the same time) in epilepsy is conceptually complex, and both decreases and increases in synchrony have been reported with ictogenesis [10]. There are multiple statistical tools to quantify correlation strength and causality among multivariate time series. In this study, we have adopted amplitude correlation (a measure of synchrony) and Granger’s Causality to estimate the coupling strength and the direction of information flow. Previous studies with magnetoencephalography (MEG) in patients with generalized seizures have reported fluctuation in synchrony with seizure progression [13, 37]. In generalized spike-wave absence seizure, there was long-range desynchrony at onset followed by local and long-range synchronization as seizure progressed. For generalized motor seizure, the increase in global synchrony (as measured by phase synchrony) was lesser compared to absence seizure. In the present study, there was a decrease in the coupling of field potentials within the thalamocortical reverberatory loop with seizures and the decrease in connectivity was directed from thalamus to hippocampus and parietal network. Mapping the ictal network dynamics might influence the temporospatial selection of stimulation parameters. Although debated, one proposed mechanism by which high frequency (>150 Hz) stimulation suppresses seizures is through desynchrony [11, 38]. If decrease in coupling within the thalamocortical network is the predominant change during ictogenesis of generalized convulsive seizure, one might speculate if increasing coupling by altering frequency (low or high frequency) or phase resetting stimulation can effectively abate seizure [3941]. Indeed low- and high-frequency stimulation have been shown to decrease seizure activity in several animal models of epilepsy [4244].

5. Study Limitations

The goal of this study was to identify potential targets for closed-loop neuromodulation of generalized tonic-clonic seizures. Accordingly, we have tested our hypotheses in a chemical model of generalized convulsive seizures and not genetic models of generalized epilepsy where spike-wave absence seizures are frequently present. Finally, we have analyzed 29 seizures recorded from two rats. Since the focus of this study was to map the dynamics of one seizure subtype (i.e., generalized tonic-clonic seizures), a higher number of subjects are less likely to impact the result.

6. Conclusions

In an acute chemical model of generalized tonic-clonic seizure, automated seizure detectors performed better for subcortical structures like thalamus or hippocampus than in cortex. The mean detection latency in thalamus by line length outperformed cortex by 4 seconds. Multiple linear correlations-Granger’s Causality revealed a decrease in coupling within the thalamocortical network during generalized tonic-clonic seizures and the decrease in information flow was significant from thalamus to hippocampal-parietal nodes. Overall, the results of this preclinical study indicate that thalamus is a superior target for automated detection of generalized seizures and modulation of synchrony to increase coupling may be a strategy for rational designing of electrical stimulation.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Puneet Dheer and Ganne Chaitanya contributed equally to this work. Kaushik Majumdar and Sandipan Pati shared senior authorship.

Acknowledgments

The authors would like to thank the UAB EEG core for access to EEG recording using the Natus Xltek machine. Sandipan Pati and Diana Pizarro would like to acknowledge the support of this study by NSF grant (NSF RII-2FECOIA 1632891). Puneet Dheer and Kaushik Majumdar would like to acknowledge the support by a DBT grant, Government of India (BT/PR7666/MED/30/936/2013).