Abstract
Postinjury epilepsy (PIE) is a devastating sequela of various brain insults. While recent studies offer novel insights into the mechanisms underlying epileptogenesis and discover potential preventive treatments, the lack of PIE biomarkers hinders the clinical implementation of such treatments. Here we explored the biomarker potential of different electrographic features in five models of PIE. Electrocorticographic or intrahippocampal recordings of epileptogenesis (from the insult to the first spontaneous seizure) from two laboratories were analyzed in three mouse and two rat PIE models. Time, frequency, and fractal and nonlinear properties of the signals were examined, in addition to the daily rate of epileptiform spikes, the relative power of five frequency bands (theta, alpha, beta, low gamma, and high gamma) and the dynamics of these features over time. During the latent pre-seizure period, epileptiform spikes were more frequent in epileptic compared with nonepileptic rodents; however, this feature showed limited predictive power due to high inter- and intra-animal variability. While nondynamic rhythmic representation failed to predict epilepsy, the dynamics of the theta band were found to predict PIE with a sensitivity and specificity of >90%. Moreover, theta dynamics were found to be inversely correlated with the latency period (and thus predict the onset of seizures) and with the power change of the high-gamma rhythm. In addition, changes in theta band power during epileptogenesis were associated with altered locomotor activity and distorted circadian rhythm. These results suggest that changes in theta band during the epileptogenic period may serve as a diagnostic biomarker for epileptogenesis, able to predict the future onset of spontaneous seizures.
SIGNIFICANCE STATEMENT Postinjury epilepsy is an unpreventable and devastating disorder that develops following brain injuries, such as traumatic brain injury and stroke, and is often associated with neuropsychiatric comorbidities. As PIE affects as many as 20% of brain-injured patients, reliable biomarkers are imperative before any preclinical therapeutics can find clinical translation. We demonstrate the capacity to predict the epileptic outcome in five different models of PIE, highlighting theta rhythm dynamics as a promising biomarker for epilepsy. Our findings prompt the exploration of theta dynamics (using repeated electroencephalographic recordings) as an epilepsy biomarker in brain injury patients.
Introduction
Postinjury epilepsy (PIE) is a devastating consequence of common brain insults, including traumatic brain injury (TBI) and stroke. It is often resistant to antiepileptic drugs and is associated with further neuropsychiatric comorbidities. Recent advancements in the understanding of PIE pathogenesis (i.e., epileptogenesis) have given rise to the hope that preventive strategies may soon be ready for clinical translation, with several studies having already demonstrated the ability to prevent experimental PIE by targeting injury-induced vascular pathology (more specifically dysfunction of the blood–brain barrier) and subsequent neuroinflammatory signaling (Zeng et al., 2009; van Vliet et al., 2012; Bar-Klein et al., 2014; Weissberg et al., 2015). However, as <20% of injured patients will develop epilepsy within a 2 year follow-up period (Pitkänen and Immonen, 2014), the large number of TBI patients required to demonstrate the statistical efficacy of these preventive treatments in the absence of biomarker-based inclusion criteria has hindered the clinical phase transition (with trials estimated to cost >$40,000,000; White et al., 2010). Moreover, all-inclusive trials are also probable to result in a higher number of patients experiencing side effects (without having needed the treatment) than patients who would benefit from PIE prevention (Friedman et al., 2014). Not only would a reliable biomarker allow adequately powered and reasonably priced trials focusing on a relevant patient population (Pitkänen et al., 2016), it would also allow sensitive identification of patients prone to future seizures as candidates for preventive treatments.
Early attempts at identifying relevant epilepsy biomarkers trace back to the middle of the previous century (Roseman and Woodhall, 1946; Cant, 1976), focusing mainly on electrographic features of neuronal hyperexcitability. While additional approaches, including genetic or anatomical screening, and blood oxygenation monitoring were explored for PIE prediction (Pitkänen et al., 2016), electrographic signals not only seem to have the greatest clinical accessibility but also the most promising potential to reflect aberrant network functionality and connectivity. Recent studies have suggested the diagnostic potential of electrographic features such as epileptiform spikes (White et al., 2010) or high-frequency oscillations (HFOs; Engel and da Silva, 2012), yet these approaches have so far failed to find clinical implementation, presumably due to limited predictive specificity (Pitkänen et al., 2016).
Having previously examined the capacity of 22 different electrographic features (reflecting time, frequency, and fractal and nonlinear signal properties) to represent ictal changes for robust detection of epileptic seizures, we demonstrated the pathological relevance of brain rhythm fluctuations (Bar-Klein et al., 2014; Levy et al., 2015; Weissberg et al., 2015; Rizzi et al., 2016). Further encouraged by findings correlating the theta rhythm in particular with epilepsy-related cognitive decline (Chauvière et al., 2009), here we set out to question the predictive potential of the previously examined 22 electrographic features and additional features of rhythm dynamics during the preseizure epileptogenic period in three models of PIE. In these three vascular injury models (Levy et al., 2015; Weissberg et al., 2011, 2015; Fig. 1), our analysis revealed that changes in theta (3–8 Hz) activity over time show great potential for epilepsy prediction. These results were further successfully validated in a rat stroke model (Lippmann et al., 2016) and in a rat electrically induced status epilepticus (EISE) model (Noè et al., 2013).
Materials and Methods
Animals.
Experiments of the mice model and the rat stroke model were conducted following approval from the ethics committee of Ben-Gurion University of the Negev, Beer-Sheva, Israel. Some of the data derive from previous studies (Levy et al., 2015; Weissberg et al., 2015; Lippmann et al., 2016). EISE experiments were performed following approval from the ethics committee of Mario Negri Institute, Milan, Italy, as published (Noè et al., 2013), and the recordings were supplied by Noè and Vezzani.
Telemetric electrocorticography and intracerebroventricular pump implantation in mice.
The procedures were performed as described previously (Levy et al., 2015; Weissberg et al., 2011, 2015). Briefly, 2- to 3-month-old mice (FVB/N mice, Harlan Biotech; and C57BL/6 mice, The Jackson Laboratory) were anesthetized with isoflurane (0.8–2%). In a stereotaxic frame, two 0.7-mm-diameter holes were drilled through the skull (all coordinates herein are relative to bregma, 3 mm posterior and 2 mm lateral) for screws serving as electrodes. Screws were connected to a wireless transmitter (Data Science International) placed in a pocket formed in the subcutaneous region of the back of the animal. One 0.7-mm-diameter hole was drilled (0.5 mm posterior, 1 mm lateral) and served for solution perfusion by an osmotic pump (ALZET) placed in the same subcutaneous pocket. The pumps were removed 7 d following the procedure. Pumps for FVB/N mice were filled with 0.4 mm bovine serum albumin (BSA; Sigma-Aldrich), 100 ng/ml transforming growth factor (TGF)-β1 (PeproTech), or 0.4 mm BSA plus 300 μm SJN2511 (a selective blocker of the TGF-β type I receptor/ALK5; Tocris Bioscience), herein, blocker. Pumps for C57BL/6 mice were filled with 5 μg/ml interleukin-6 (IL-6; PeproTech) or 0.4 mm BSA. All solutes were dissolved in artificial CSF (ACSF) prepared as described previously (Seiffert et al., 2004). Sham control animals were implanted with pumps containing ACSF or 0.4 mm fluorescein isothiocyanate-conjugated 70 kDa dextran (Sigma-Aldrich) in ACSF. After recovery, animals were placed in a behavior room maintained on a 12 h circadian cycle for 14 d with food and water accessible ad libitum. A receiver placed under the cage of each animal collected (sampling rate, 1 KH) and transmitted the signals to a personal computer.
Locomotor activity measurement.
Locomotor activity in freely moving mice was measured as described previously (Tarasiuk et al., 2014; Jukic et al., 2015) using the implanted wireless transmitter (Data Science International) as mentioned above. Briefly, the derivative of the signal strength received from the implanted transmitter is correlated with the locomotion of the animal. To assess circadian rhythmicity, the derivative was divided into epochs of 2 min with 1 min overlap, and the mean was calculated for each epoch. Then, a moving average with a constant of 500 was applied (Fig. 1A,B).
Telemetric deep recordings and photothrombosis-induced cortical ischemia.
Procedures were performed as described previously (Lippmann et al., 2016). Briefly, for electrode implantation, rats (males; weight, 200–300 g; Harlan Biotech) were anesthetized by ketamine-xylazine (1.6 and 0.6 mg/kg body weight, i.p., respectively) and maintained by isoflurane (0.8%) via a mask in a stereotaxic frame. After drilling three holes in the skull, a Teflon-coated platinum iridium wire 33 μm in diameter (Science Products) connected to the DSI transmitter wire was located 2.5 mm lateral, 3.5 mm caudal, and 2.2 mm ventral into the CA1 region of the right hippocampus (Swanson, 2004). The second electrode was fixed to a screw above the left cerebellum located 4.5 mm lateral and 5 mm dorsal to serve as an epidural reference. After a 4 d recovery period, a 48 h baseline recording was performed followed by the induction of stroke using Rose Bengal (Watson et al., 1985; Lapilover et al., 2012). Rose Bengal was administered (20 mg/kg, i.v.) to the tail vein, and a halogen light beam (3.5 mm diameter) was directed on the exposed skull (1 mm posterior and 1 mm lateral to bregma) for 15 min. The skin was sutured and treated with xylocain gel. Buprenorphine (0.05 mg/kg, s.c.) was administered. Recording continued for 9 d after photothrombosis (PT). A single animal was excluded due to multiple acute seizures occurring <24 h after the procedure.
Telemetric deep recordings and EISE.
For experimental details, see Noè et al. (2013). Briefly, male adult rats (weight, ∼300 g; Charles River Laboratories) were implanted under anesthesia with two bipolar depth electrodes inserted into the temporal poles of both hippocampi (posterior, 4.7 mm; lateral, ±5; ventral, 5). For ground and reference electrodes, screws were placed in the nasal sinus and cerebellum, respectively. Following 7–10 d of recovery, 24 h of baseline recording was performed in freely moving rats. Then, rats were stimulated unilaterally in the CA3 region of the ventral hippocampus for 90 min to induce status epilepticus (SE) as described previously (50 Hz; 400 μA peak-to-peak; 1 ms biphasic square waves in 10 s trains delivered every 11 s; De Simoni et al., 2000; Noe et al., 2008). Following the induction, rats were recorded for 8–21 d at 400 Hz. For the present study, recordings from 13 rats were analyzed, and 4 were excluded due to seizures detected before the fifth day after the SE (acute seizures).
Signal processing for seizure detection.
For the unbiased detection of seizures and determination of epilepsy onset electrographic signals were analyzed using an in-house MatLab software as described previously (Bar-Klein et al., 2014) followed by a blind manual review. ECOG analysis was based on features extraction and artificial neuronal network (ANN). To train the detector, a 1.5 h ECOG segment composed of background activity, noise, artifacts, and 24 seizures lasting 10–60 s each recorded from three models of epilepsy: Pilocarpine, genetic (synapsin triple knockout), and albumin (Becker et al., 2008; Ketzef et al., 2011; Weissberg et al., 2011). Data were bandpass filtered (2–100 Hz) and buffered into 2 s epochs with 1 s overlap. References ANN vectors of both seizures and baseline activities were assigned with values of 1 and −1, respectively. Twenty-two features previously used in brain signal analysis were extracted from each epoch. Then, a forward sequential selection algorithm (Whitney, 1971) was implemented to find the optimal features combination that will enable both minimal computational complexity and maximal quality of detection. The chosen subset includes energy, 12–20 Hz relative power (β), 20–40 Hz relative power (low gamma), SD, curve length. These features were used to train an ANN (single hidden layer, n = 50, back-propagation ANN; Rumelhat et al., 1986) with Powell–Beale conjugate gradient training (Powell, 1977). Following receiver operating characteristic (ROC) analysis, a threshold of 0.85 was set. A seizure was defined when at least five suprathreshold consecutive epochs were detected.
Signal processing for epileptogenesis detection.
We examined the 22 features explored for seizure detection, among them the mean relative power in five frequency bands: theta (3–8 Hz); alpha (8–12 Hz); beta (12–20 Hz); low gamma (20–40 Hz); and high gamma (30–100 Hz; 30–70 Hz in the EISE model) for every 60 s of ECOG segment. We also tested the dynamic changes in these power bands during the epileptogenic period using linear regression, with the calculated slopes representing the extent of the change. Recordings were analyzed from the following: mice: intracerebroventricular albumin, FVB/N (N = 10); C57BL/6 (N = 4); intracerebroventricular TGF-β (N = 3); albumin with TGF-β blocker (N = 7; Weissberg et al., 2015); IL-6 (N = 5); and ACSF/dextran (N = 9; Levy et al., 2015; Weissberg et al., 2015); rats: PT rats (N = 13; Lippmann et al., 2016); and EISE rats (N = 9; Noè et al., 2013). One control animal was excluded as an outlier (due to noise resulting in a theta slope magnitude larger than the mean ± 2.7 SDs of its group).
Signal processing for spike detection.
ECOG recordings were bandpass filtered between 1 and 45 Hz. Events higher than 2 times and lower than 10 times of the signal baseline were detected for each window of 1024 ms. Detected events contained both epileptic spikes and other types of physiological/nonphysiological noise. To detect epileptic spikes, a wavelet transform (WT)-based method was implemented (Adeli et al., 2003). Automatic detection and clustering was followed by a blind human review. Detection was performed on days 2–4 after the surgery. Four animals from the Alb plus blocker group were excluded from all the presented analyses due to a highly variable and noisy signal that did not allow a reliable spike detection.
Statistics.
ROC analyses and k-means clustering were performed using MatLab. Other statistical tests were performed using GraphPad Prism software. Relative power bands of ECOG and duration of epileptogenesis were compared by Mann–Whitney test. Survival analyses were performed using the Kaplan–Meier method with curves comparison using the Mantel–Cox test. In correlation tests Pearson's r is presented. Fisher's exact test was implemented to compare the early-onset and late-onset seizures of the EISE model. All statistical tests are two tailed.
Results
Intracerebroventricular administration of albumin, TGF-β, or IL-6 as mechanism-based models of postinjury epilepsy
To identify a disease-relevant biomarker, we set out to explore electrographic features in mechanism-driven models of PIE. In addition to discovering the epileptogenic role of albumin, which is shown to activate TGF-β signaling and a neuroinflammatory response (Cacheaux et al., 2009) with predominant IL-6 upregulation (Fig. 1A; Levy et al., 2015; Weissberg et al., 2011, 2015) and mimicking different stages of the pathogenic process following injury-induced blood–brain barrier (BBB) dysfunction, we show that intracerebroventricular infusion of albumin, TGF-β, or IL-6 results in seizures in 71%, 100%, and 60% of the animals, respectively. Notably, while a few animals presented isolated seizures within 2 d after treatment, in most of the animals, spontaneous seizures developed after a latent period of >4 d. No seizures were detected in control (sham-implanted) animals (N = 9), and in six of seven mice cotreated with albumin and a TGF-β blocker (Fig. 1B–D). Interestingly, when considering all animals in the albumin-treated group, mice of the FVB/N strain showed significantly shorter epileptogenic periods compared with the C57BL/6 strain [4.94 ± 0.53 d (N = 12) vs 12.08 ± 1.4 d (N = 3), respectively; p = 0.004]. Seizures detected before 4 d following the insult were regarded as “acute seizures”. An animal was considered epileptic if recurrent seizures were detected 4 d after perfusion onset or later. Seizures were documented over 30 d after perfusion onset. During most of the seizures, animals showed staring behavior and only rarely presented with generalized tonic–clonic manifestation (Weissberg et al., 2011, 2015). Bilateral cortical recordings performed in a few animals confirmed that seizures first occurred in the right (implanted) hemisphere and propagated to the contralateral cortex (data not shown).
Mechanism-based models of PIE. A, The mechanisms mediating epileptogenesis in the studied PIE models include the following: serum albumin (Alb) extravasation (in regions with dysfunctional BBB) → activation of astrocytic TGF-β receptors → alk5-mediated smad2/3 phosphorylation → transcriptional inflammatory changes promoting further TGF-β1 and IL-6 secretion → reorganization of the neural network → spontaneous seizures. B, A representative spontaneous seizure, as recorded in a mouse treated intraventricularly with albumin. C, Plot of the accumulating number of seizures per day in the six treatment groups demonstrates a minimum duration of 4 d as the epileptogenic period. D, Distribution of seizure onsets in the different models. dex, Dextran; B6, C57BL/6.
Daily rate of epileptiform spikes shows limited potential as a PIE predictive biomarker
The latent period greatly varies between models of PIE (Pitkänen and Immonen, 2014), and it is not known when network modification of epileptogenesis starts. Evidence shows molecular and electrophysiological changes as early as 72 h after the injury (El-Hassar et al., 2007; Yaari et al., 2007; Levy et al., 2015; Weissberg et al., 2015). Thus, taking a clinically relevant approach that strives to detect a pathogenesis in its earliest stage, we focused our analysis on the electrographic changes during days 2–4 after treatment onset.
As the epileptogenic period is often characterized by epileptiform spikes (White et al., 2010), we set out to test the predictive potential of a feature representing the daily rate of such spikes. ECOG spikes were detected automatically using WT, and the mean number of spikes per day (during days 2–4; i.e., the preseizure period) was calculated for each animal (Fig. 2D). Generally, epileptic animals had significantly higher spike rates (123.3 ± 45.46 vs 274.3 ± 46.04 spikes/d; p = 0.026); however, rate values fluctuated greatly and lacked consistency throughout the epileptogenic period. ROC analyses of the predictive accuracy of spike rates at day 2, 3, and 4 revealed the following values: day 2: area under the curve (AUC) = 0.55, p = 0.71, and 95% confidence interval (CI) = 0.03–0.79; day 3: AUC = 0.74, p = 0.013, and CI = 0.58–0.91; day 4: AUC = 0.71, p = 0.034, and CI = 0.53–0.89 (Fig. 2E).
Spike rate and nondynamic representation of brain rhythms do not predict epilepsy. A–C, The daily rate of epileptiform spikes (see inset) for each animal calculated for days 2, 3, and 4 after treatment initiation in each treatment group. D, Mean number of spikes per day for each animal. E, ROC analysis reveals AUCs up to 0.74 for days 3 and 4 and the mean of days 2–4, with only a 0.55 AUC for day 2 (not significantly different from randomized classification; AUC of 0.5 dotted black line). F, Spectral analysis during the fourth day shows no significant differences in the tested frequency bands when comparing epileptic to nonepileptic animals. Solid lines represent the mean value, and background color depicts the SD. Dex, Dextran; ALB, albumin; Lγ, low gamma; Hγ, high gamma.
Features of rhythm-specific activity show meager potential as PIE biomarkers
In addition to reports of theta band changes in SE models of epilepsy (Arabadzisz et al., 2005; Chauvière et al., 2009; Marcelin et al., 2009), we next calculated the relative power of ECOG activity in the delta, theta, alpha, beta, and low- and high-gamma rhythms (during the preseizure period of days 2–4). Overall, we found high variability within the groups with no significant differences between epileptic (N = 17) and nonepileptic (N = 19) animals. An exemplary depiction of rhythm-specific activity at day 4 (light time) demonstrates the overlap in values between the epileptic and nonepileptic animals (delta, p = 0.68; theta, p = 0.82; alpha, p = 0.71; beta, p = 0.35; low-gamma, p = 0.27; high-gamma, p = 0.42; Fig. 2F).
Theta dynamics as a biomarker for epilepsy
Despite the inability of frequency-band features to reliably reflect pathology, we hypothesized that changes in rhythmic activity over time may offer additional insight into disease development (expecting consistent activity in nonepileptic animals and rhythmic shifts associated with the epileptic transformation of the neuronal network). Among the examined bands, theta dynamics were found to offer the greatest promise, with small variability in nonepileptic animals and prominent changes in epileptic mice (Fig. 3A). Linear regression analyses further revealed a “safety range” of nonepileptic theta slopes (ranging around zero), with epileptic animals clearly falling outside this range (Fig. 3C). A ROC analysis was next used to evaluate the predictive potential of two measures of theta dynamics: the absolute slope value of each animal (sensitivity I) and its distance from the mean of control group (sensitivity II), revealing AUC = 0.9102 with CI = 0.81–1.01 and AUC = 0.9969 with CI = 0.99–1.01, respectively (p < 0.0001 for both; Fig. 3B).
Theta dynamics as a biomarker for epilepsy. A, Theta activity during the epileptogenic period from representative epileptic and nonepileptic animals (N = 3 in both groups). The smoothed raw traces and the corresponding linear fits are represented by continuous vs dashed lines, respectively. B, ROC analysis was used to assess the predictive potential of two measures of theta dynamics: absolute value of the slope (sensitivity I: AUC = 0.9102, p < 0.0001) and the distance of each slope from the mean slope of control group (sensitivity II: AUC = 0.9969, p < 0.0001). The dotted black line indicates an AUC of 0.5, a reference to randomized classification. C, Slopes of theta activity, superimposed on the margins of the safety range (dashed lines) calculated based on the OOP of sensitivity II. D, Kaplan–Meier analysis shows a significant difference between the OOP-classified groups (p < 0.0001, based on sensitivity I or II). W/I, Within safety range; O/O, out of safety range. E, Theta slopes are inversely correlated with the duration of the epileptogenic period (N = 17, r = 0.66, p = 0.0036). F, Theta and high-gamma slopes are inversely correlated in both epileptic (N = 17, r = 0.924, p < 0.0001) and nonepileptic (N = 19, r = 0.51, p = 0.021) animals. dex, Dextran. G, ROC analyses of a single sample per day (1-h-long segment): AUC = 0.66 (p = 0.1 CI = 0.46–0.86) and AUC = 0.75 (p = 0.018, CI = 0.58–0.92), respectively, to sensitivity I and II. H, ROC analyses of two samples per day (1-h-long segment, every 12 h) AUC = 0.78 (p = 0.004, CI = 0.63–0.94) and AUC = 0.873 (p = 0.0001, CI = 0.76–0.99).
The optimal operational point (OOP) was then calculated for each ROC curve (as the curve point with maximal distance from randomized classification; Fig. 3B, dashed line). Comparison of Kaplan–Meier survival curves of animals divided based on the identified OOPs thresholds further confirmed significant differences between the epileptic and nonepileptic mice (Fig. 3D; p < 0.0001).
As the false-negative results identified animals that had relatively late onsets, we also examined the correlation between the theta slopes and the duration of the epileptogenic period. Indeed, higher slope values were found to correlate with earlier onset of epilepsy (spontaneous seizures, r = 0.66, p = 0.0036; Fig. 3E). An inverse correlation was observed between the theta and high-gamma slopes but not lower-frequency bands (delta, alpha, or beta). The coefficient correlation was greater among epileptic (N = 17, r = 0.924, p < 0.0001; Fig. 3F) compared with nonepileptic animals, regardless of treatment (N = 19, r = 0.51, p = 0.021; Fig. 3F). These results further highlight the relevance of theta activity to the epileptogenic transformation of the neural networks. Interestingly, 23.5% (N = 4) of epileptic animals showed a positive rather than a negative theta slope during days of early epileptogenesis. In these animals, this increase was followed by a decrease during the next few days of recordings, such that a few days before the first seizure, theta relative power decreased (mean ± SEM slope, −0.0276 ± 0.0079; data not shown).
Since continuous EEG recordings during the latent period in human patients is difficult to implement, we next repeated the analysis using a “clinically relevant” approach by sampling 1-h-long ECOG segments (every 24 h starting from the first hour of the second day for 4 consecutive days). Analysis revealed an AUC of 0.66 (p = 0.1, CI = 0.46–0.86) and an AUC of 0.75 (p = 0.0118, CI = 0.58–0.92), respectively, for sensitivity I and II (Fig. 3G). When using 1-h-long segments every 12 h (seven segments), the AUCs were 0.78 (p = 0.004, CI = 0.63–0.94) and 0.873 (p = 0.0001, CI = 0.76–0.99), respectively, to sensitivity I and II (Fig. 3H), supporting the translational potential of repeated theta power measure.
Circadian rhythm is disturbed during epileptogenesis
Changes in theta rhythm power are known to be associated with sleep/awake cycle and behavioral state in both rodents (Winson, 1974, 1978; O'Keefe and Recce, 1993) and human (Demenet and Kleitman, 1957; Kahana et al., 1999). We thus evaluated locomotor activity and compared the diurnal rhythmicity of theta changes with locomotor activity during epileptogenesis (N = 8) and in nonepileptic animals (N = 8). Nonepileptic animals (Fig. 4A) demonstrated a significantly greater correlation between theta relative power and locomotor activity signal compared with animals undergoing epileptogenesis (Fig. 4B; median values, −0.41 vs −0.044, respectively; p = 0.0003; data not shown). This finding is consistent with the notion that slow rhythms are abundant during sleep with respect to wakefulness in preserved circadian rhythm (Rechtschaffen and Kales, 1968). Fast Fourier transform (FFT) analysis of locomotor activity revealed the dominant frequency of 1 cycle/d in nonepileptic mice, with no significant reduction in animals undergoing epileptogenesis (p = 0.19; Fig. 4C). A similar analysis for frequency changes in theta rhythm revealed a significant reduction in the dominant frequency in animals undergoing epileptogenesis (p = 0.0378; Fig. 4D), suggesting a greater sensitivity of the theta band to disruption of the circadian cycle. Notably, the mean day and night differences in theta power (days 2–4, pick to pick, as seen in Fig. 4A, top trace) among nonepileptic (N = 8) animals was smaller compared with the differences in theta power among animals undergoing epileptogenesis (N = 8, start to end, as seen in Fig. 4B, top trace; p = 0.0378; Fig. 4E).
Circadian rhythm is disturbed during epileptogenesis. A, B, Representative examples of smoothed theta relative power and locomotor activity signals in nonepileptic mice (A) or in mice undergoing epileptogenesis (B). C, D, Comparison of FFT analysis applied over the locomotor activity (C) and theta relative power (D) signals plotted in A and B between nonepileptic mice (N = 8) and mice undergoing epileptogenesis (N = 8). Comparison of the mean pick-to-pick theta relative power difference among nonepileptic animals (N = 8) and the difference between the start and end of epileptogenesis among animals undergoing epileptogenesis. E, Presented are the median values plus the range. *p = 0.0378.
Theta dynamics as an epileptogenesis biomarker in additional PIE models
To test whether our findings in mice models of vascular PIE could be generalized to other PIE models, we next analyzed recordings from the PT model of poststroke epilepsy (Lippmann et al., 2016) and the EISE model (Noè et al., 2013).
As the PT rat model is associated with early-onset spontaneous seizures (Lippmann et al., 2016), theta dynamics were quantified for the first 24 h after PT, and only in rats without seizures during that period (Fig. 5A,B). No seizures were detected in the baseline period recorded before PT. Consistent with our findings in mice, nonepileptic animals showed little daily variability, whereas theta slopes of epileptic rats fluctuated significantly (Fig. 5C) and, similarly, were inversely correlated with high-gamma dynamics (r = 0.833, p = 0.0004; data not shown). The distance of the calculated slopes from the nonepileptic mean was used for ROC analysis and resulted in an AUC of 1 (p = 0.003, CI = 1–1, Fig. 5D). Kaplan–Meier survival analysis revealed significantly different curves of animals within and out of the safety range derived from the OOP of the ROC analysis (p = 0.04; data not shown).
Theta dynamics as a biomarker for epilepsy in a stroke model of PIE. A, An example of a hippocampal seizure in a model of poststroke epilepsy. B, A heat map of daily seizure rates shows that 5 of the 14 stroke animals remained nonepileptic throughout the recorded period. Among the epileptic animals, the seizure-free period ranged between 1 and 7 d. C, D, Theta slopes (during the first 24 h from treatment initiation) of nonepileptic animals fall within a distinctive range (C), the distance from which can accurately classify disease outcome, as evidenced in the ROC analysis (D), revealing an AUC of 1 (p = 0.003). The dotted black line indicates an AUC of 0.5, a reference to randomized classification.
In the EISE model (rats monitored with bilateral intrahippocampal electrodes at the Milan laboratory), theta slopes were calculated for the ipsilateral electrode for days 3–5 after SE induction. As this robust model induces seizures in all animals, theta slopes were used to classify animals into early versus late epilepsy onset (7–8 vs 11–18 d, respectively, classified based on the overall onset median at day 11). No seizures were detected in the baseline period recorded before the electrical stimulation. Unsupervised k-means clustering divided the slopes into a group of six animals with zero-ranging slopes, consisting predominantly of late-onset animals (five of six; p = 0.046, Fisher exact test; or AUC = 0.85, p = 0.086, CI = 0.56–1.14; ROC analysis; data not shown); and a group of higher absolute slopes (two negative slopes and one positive slope; Fig. 6A), representing early-onset rats. Linear correlation between seizure onset and absolute theta slope yielded a trend like that obtained in the mice models (Fig. 6E; r = 0.653 and p = 0.054, data not shown) and revealed a negative theta–high-gamma correlation (p < 0.0001, data not shown). Interestingly, when comparing slopes between the two hippocampi, most rats showed similar theta fluctuations and seizure propagation from one hippocampus to the other (Fig. 6C–E, outlined square). In contrast, in a rat with stable theta values (outlined circle) seizures did not necessarily propagate to the contralateral hippocampus (Fig. 6E), supporting the hypothesis that theta changes are indicative of network modifications associated with a lower seizure threshold.
Theta dynamics as a biomarker for epilepsy in an EISE model of PIE. A, B, The k-means classification (dashed lines) of theta slopes (as recorded from the electrode ipsilateral to stimulation (A) resulted in a division corresponding to the latency of the epileptogenic period (B), with higher absolute slopes associated with earlier epilepsy onsets, whereas slow theta changes corresponded to slower disease development (p = 0.035). C, Comparison of ipsilateral and contralateral theta slopes in early-onset animals highlights one animal (outlined circle) with remarkably different values between the two hemispheres. D, E, Interestingly, this animal (E) also had evident differences in intrahemispheric seizure propagation compared with an animal with relatively consistent slopes [D (outlined square in C)].
Discussion
We tested electrophysiological features in five rodent models of epilepsy (mice and rats) and found that theta activity dynamics show great promise as a diagnostic biomarker for epileptogenesis, a prognostic biomarker for PIE, and a pharmacodynamics biomarker to evaluate the efficacy of antiepileptic drugs. We further show that the degree of theta shift may also serve as a prognostic indicator for the duration of the epileptogenic period, with steeper slopes reflecting the earlier occurrence of seizures.
The potential of electrographic features to predict epilepsy has been explored in both patients and animals (Roseman and Woodhall, 1946; Cant, 1976; Gotman, 1991; White et al., 2010). Nonseizure epileptiform spikes were previously suggested to show predictive promise to seizures or epilepsy (Gotman, 1991; White et al., 2010). Our results confirm that spikes are more frequent during epileptogenesis; however, the high daily variability in spike rate (among both epileptic and nonepileptic animals) leads to a false-positive rate that limits the predictive specificity of this feature (Fig. 2). The detection of spikes in nonepileptic animals is consistent with previous studies (Bragin et al., 1995; White et al., 2010) and has been previously explained by noise or physiologic sources with a spike-like shape such as physiological dentate spikes (Bragin et al., 1995). Limited reliability of spike detectors may also underlie the inability of spikes to serve as a biomarker for epileptogenesis. In light of the limited frequency range of our recording systems, we were unable to assess the previously suggested diagnostic potential of HFOs (80–600 Hz; Engel and da Silva, 2012). However, while such oscillations were shown to benefit focus localization once epilepsy is confirmed (Engel and da Silva, 2012), their potential for the prediction of epilepsy was suggested to be limited (Pitkänen et al., 2016).
Interestingly, several studies have attracted our attention to the pathological relevance of theta activity, showing its decrease during epileptogenesis in the pilocarpine and kainic acid models (Arabadzisz et al., 2005; Chauvière et al., 2009; Marcelin et al., 2009) and implicating this decrease in cognitive failure (Chauvière et al., 2009). Importantly, our results point to a limited biomarker potential of theta activity when assessed at a single point in time, as a single measurement is susceptible to differences arising from electrode placement and variability between models. In contrast, we demonstrated predictive generalization when taking several measurements from each animal, quantifying the dynamic shift in theta activity rather than independent values. In most cases, a decrease in theta oscillations was found during the epileptogenic period, with greater decreases correlated with earlier onset of epilepsy (Figs. 3E, 6A,B).
The potential of theta dynamics as a general biomarker for epileptogenesis was shown by demonstrating a similar pattern in mice and rats as well as in intrahippocampal recordings following PT-induced stroke or EISE (Fig. 6A,B), with greater theta slopes associated with significantly shorter latent periods. Moreover, in the EISE model, the ipsilateral hemisphere that underwent epileptogenesis showed a high theta slope, whereas in the contralateral hemisphere the theta slope was near zero (Fig. 6C–E). This finding strengthens the linkage between theta slope and epileptogenesis. Notably, the theta/PIE onset association was evident in four of the five models, excluding stroke-induced PIE. Such a lack of correlation may be due to the relatively short latent period in this model or because this model was the only one using ketamine anesthesia in addition to isoflurane. These reasons may underlie the necessity of comparing the theta slope of the average slope of nonepileptic animals in the PT model as opposed to the other models where the absolute value of the slope is informative as is. The discovered association between greater theta shifts and earlier PIE onsets is complemented by the previously suggested correlation between the loss of theta activity and memory deficits in experimental epileptic animals (Chauvière et al., 2009). In addition, a cognitive deficit itself was suggested as a predictive biomarker of PIE (Pascente et al., 2016). Interestingly, they both showed it in the pilocarpine-induced SE model of epilepsy. Together, these findings support the ability of theta rhythm to reflect clinically relevant network modifications that may not only underlie seizures but also contribute to cognitive decline, a hallmark of PIE and BBB dysfunction (Helmstaedter et al., 2003; Montagne et al., 2015). Future studies in models that do not result in 100% epilepsy are required to test whether changes in theta rhythm that are associated with cognitive decline can be distinguished from theta changes associated with epileptogenesis.
The onset of seizures in the PT model was relatively early (24–144 h after the stroke; Fig. 5B). In a previous study of the PT model, Kelly et al. (2001) reported seizures onset of 2 months after stroke. While early continues recordings were not performed by Kelly et al. (2001), and no long-term recordings were made in the data available for the present study, we cannot rule out that the detected seizures represent an acute response to the ischemic injury. However, the theta dynamics observed during the short postinjury period suggest that these represent a true reduction in seizure threshold within the local cortical network, regardless of how fast they develop.
The lack of theta changes in animals treated with intracerebroventricular albumin together with TGF-β blockers (Fig. 3C) further suggests that theta dynamics may also serve as a pharmacodynamics biomarker, allowing the follow-up of preventive treatment efficacy. The observation that the sampling of 1 h recordings was sensitive in predicting epilepsy further indicates the translational potential of our findings.
Interestingly, four of the epileptic animals (23.5%) in the minipump models showed an increase, rather than a decrease, in theta during days 2–4. As mentioned, the theta increase was followed by a decrease that preceded the first seizure. Previous reports showed that an increase in theta is a marker of sleep propensity under or following sleep deprivation (Aeschbach et al., 1997; Strijkstra et al., 2003; Vyazovskiy and Tobler, 2005). Our results are showing disturbed circadian rhythm during epileptogenesis (Fig. 3) that may explain the initial increase in theta.
Circadian distortions were previously described during the chronic period of epilepsy in both patients and rodents (Bastlund et al., 2005; Pung and Schmitz, 2006; Choi et al., 2016; Quigg et al., 2016). However, to our knowledge, this is the first report of disturbed circadian rhythms in animals undergoing epileptogenesis compared with nonepileptic animals. These results, while they should be confirmed in additional models and in human patients, call for the attention of clinicians and scientists to additional disturbances in brain activity among patients at risk for PIE.
While we show that theta changes predict PIE, we also found an increase in high-gamma activity (30–100 Hz) and a robust inverse correlation between the two bands. Under physiological conditions, gamma and theta oscillations are expected to behave in a phase- and amplitude-coupled manner (Buzsáki et al., 1983), and stimulation at theta frequencies was shown to induce gamma oscillations (Butler et al., 2016). Thus, our findings of strong negative theta–gamma correlation (r = 0.92; Fig. 3F) further highlight a transition from physiological to pathological network activity. An inverse correlation between theta and gamma relative power slopes was also observed in the nonepileptic group, but to a much lesser extent (r = 0.51; Fig. 3F). This may derive from the effects of surgical trauma or anesthesia rather than epileptogenesis. Importantly, increased gamma activity was previously associated with seizures (Köhling et al., 2000; Medvedev et al., 2000; Medvedev, 2002; Willoughby et al., 2003; Kobayashi et al., 2004, 2010; Herrmann and Demiralp, 2005; Traub et al., 2005; Wu et al., 2008; Inoue et al., 2008; Yamazaki et al., 2009; Andrade-Valenca et al., 2012; Maheshwari et al., 2016) via the glutamatergic (Medvedev et al., 2000; Medvedev, 2002) NMDA pathway (Maheshwari et al., 2016). Thus, our evidence of increased gamma activity, together with the previously reported glutamatergic synaptogenesis in the albumin/TGF-β PIE models (Weissberg et al., 2015) further implicate a process distorting physiological brain rhythms and connectivity, potentially through the modification of synaptic connections and reorganization of the neural network.
Our finding that C57BL/6 mice showed a slower theta decline and delayed seizure onsets compared with FVB/N mice (Figs. 1C,D, 3E) is consistent with our previous report on low IL-6 response to TGF-β in C57BL/6 mice (Levy et al., 2015) and their higher threshold for seizures and seizure-induced cell death (Borges et al., 2003; Yang et al., 2005; McLin and Steward, 2006; Bankstahl et al., 2012; Schauwecker, 2012). This relative strain-dependent resilience to an epileptogenic insult calls for future studies into the mechanisms underlying such a resilience.
In summary, we propose that theta rhythm alterations represent network modifications associated with a lower threshold for seizures and that this simple feature has a potential to become a reliable, noninvasive biomarker for epileptogenesis following brain insult. To actualize this potential, prospective clinical trials of repetitive EEG recordings in patients following brain insults are needed. We assert that these recordings have a high potential to reveal relative shifts in brain rhythms, serving as a biomarker for antiepileptogenic clinical studies.
Footnotes
This study was supported by the European Union Seventh Framework Programme (FP7/2007-2013; Grant agreement 602102, EPITARGET; to A.F.), the Israel Science Foundation (to A.F.), the Binational Israel-USA Foundation (to A.F.), The German Research Foundation (to A.F.), the CURE Epilepsy Foundation, the Nova Scotia Health Research Foundation, and the Canada Foundation for Innovation.
The authors declare no competing financial interests.
- Correspondence should be addressed to Alon Friedman, Department of Medical Neuroscience, Dalhousie University, Halifax, NS B3H 4R2, Canada. alon.friedman{at}dal.ca
