# Does the Thalamo-Cortical Synchrony Play a Role in Seizure Termination?

Elisa Evangelista, Christian Bénar, Francesca Bonini, Romain Carron, Bruno Colombet, Jean Régis, Fabrice Bartolomei

► To cite this version:

Elisa Evangelista, Christian Bénar, Francesca Bonini, Romain Carron, Bruno Colombet, et al. Does the Thalamo-Cortical Synchrony Play a Role in Seizure Termination?. Frontiers in Neurology, 2015, 6, pp.192. 10.3389/fneur.2015.00192. hal-01785633

HAL Id: hal-01785633

https://amu.hal.science/hal-01785633

Submitted on 4 May 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

CrossMark

# Does the thalamo-cortical synchrony play a role in seizure termination?

Elisa Evangelista<fnref n="1" />, Christian Bénar<fnref n="2" />,³, Francesca Bonini<fnref n="1" />,<fnref n="2" />,³, Romain Carron<fnref n="2" />,³,⁴, Bruno Colombet<fnref n="2" />,³, Jean Régis<fnref n="2" />,³,⁴ and Fabrice Bartolomei<fnref n="1" />,<fnref n="2" />,³*

The mechanisms underlying seizure termination are still unclear despite their therapeutic importance. We studied thalamo-cortical connectivity and synchrony in human mesial temporal lobe seizures in order to analyze their role in seizure termination. Twenty-two seizures from 10 patients with drug-resistant mesial temporal lobe epilepsy undergoing pre-surgical evaluation were analyzed using intracerebral recordings [stereoelectroencephalography (SEEG)]. We performed a measure of SEEG signal interdependencies (non-linear correlation), to estimate the functional connectivity between thalamus and cortical regions. Then, we derived synchronization indices, namely global, thalamic, mesio-temporal, and thalamo-mesio temporal index at the onset and the end of seizures. In addition, an estimation of thalamic "outputs and inputs" connectivity was proposed. Thalamus was consistently involved in the last phase of all analyzed seizures and thalamic synchronization index was significantly more elevated at the end of seizure than at the onset. The global synchronization index at the end of seizure negatively correlated with seizure duration (p = 0.045) and in the same way the thalamic synchronization index showed an inverse tendency with seizure duration. Six seizures out of twenty-two displayed a particular thalamo-cortical spike-and-wave pattern at the end. They were associated to higher values of all synchronization indices and outputs from thalamus (p = 0.0079). SWP seizures displayed a higher and sustained increase of cortical and thalamo-cortical synchronization with a stronger participation of thalamic outputs. We suggest that thalamo-cortical oscillations might contribute to seizure termination via modulation of cortical synchronization. In the subgroup of SWP seizures, thalamus may exert a control on temporal lobe structures by inducing a stable hypersynchronization that ultimately leads to seizure termination.

**Keywords:** focal seizures, thalamus, intracerebral EEG, synchrony, seizure termination

## Introduction

Epileptic seizures are paroxysmal recurrent and self-limiting events, lasting no more than a few minutes. Significant research efforts have directed toward a better understanding of the electro-physiological and neurodynamic correlates of seizure initiation, propagation, and termination (1-6). Nevertheless, at present, the intrinsic mechanisms underlying seizure termination are still unclear despite their therapeutic importance. Identifying mechanisms that result in seizure termination would lead to the development of novel and more efficient interventional approaches aimed to

**OPEN ACCESS**

**Edited by:**
Patrick A. Forcelli,
Georgetown University, USA

**Reviewed by:**
Fernando Henrique Lopes Da Silva,
University of Amsterdam,
Netherlands

Andrei V. Medvedev,
Georgetown University, USA

\*Correspondence:
Fabrice Bartolomei,
Service de Neurophysiologie Clinique,
CHU Timone, 264 Rue st Pierre,
Marseille 13005, France
fabrice.bartolomei@ap-hm.fr

**Specialty section:**
This article was submitted to
Epilepsy, a section of the
journal Frontiers in Neurology

**Received:** 28 June 2015
**Accepted:** 18 August 2015
**Published:** 01 September 2015

**Citation:**
Evangelista E, Bénar C, Bonini F,
Carron R, Colombet B, Régis J and
Bartolomei F (2015) Does the
thalamo-cortical synchrony play a
role in seizure termination?
Front. Neurol. 6:192.
doi: 10.3389/fneur.2015.00192

<footnote n="1">¹Service de Neurophysiologie Clinique, CHU Timone, Assistance Publique des Hôpitaux de Marseille, Marseille, France,</footnote>

<footnote n="2">²UMR1106, INSERM, Marseille, France, ³Institut de Neurosciences des Systèmes Marseille, Aix Marseille Université, Marseille, France, ⁴Service de Neurochirurgie Fonctionnelle et Stéréotaxie, Assistance Publique des Hôpitaux de Marseille, Marseille, France</footnote>
induce or reproduce the same termination mechanisms that are triggered in self-limited seizures.

It is known that mechanisms of seizure termination range in scale from the regulation of transmembrane potentials at the level of single neurons to local neuron and interneuron networks and to long-range cortico-subcortical modulating networks comprising remote structures, such as the brain stem and basal ganglia. These long-range interactions may be reflected by an increased synchronization between distant cortical regions and between cortex and subcortical regulatory circuits during seizures (7).

The exact relationship between synchronization and seizure termination is not completely understood. In an extensive series of *in vivo* cat experiments including multi-site intracellular and extracellular as well as local cortical field potential recordings, it has been observed that there was a greater increase in intrahemispheric synchrony immediately before seizure termination (5, 8). Similarly, in recent studies on human partial seizures utilizing different mathematical methods, cortical correlation increased before the seizure ended (9, 10). It has been thus suggested that the hypersynchrony observed at the end of seizures may play a crucial role in seizure termination.

The role of subcortical structures in modulating seizure termination remains an unanswered question. Subcortical structures play a significant role in modulating seizure activity, in terms of cortical seizure threshold, duration, or severity. Yet, the role of synchronized thalamo-cortical oscillations has been well demonstrated in animal models of absence seizures, where thalamus was shown to provide a resonant circuitry that organizes, amplifies, and synchronizes seizure activity (4, 11).

The involvement of thalamus during partial seizures has been previously described in animal models of temporal lobe seizures (12–18). Electrophysiological recordings of spontaneous limbic seizures in chronic epileptic rats showed that circuits involving the midline thalamus and limbic structures are activated in the early stages of seizure and that the infusion of lidocaine in thalamic nuclei shortens seizure duration (16). Thus, authors suggested that thalamus may be necessary to amplify and distribute ictal activity throughout the system and may play a significant role in seizure modulation.

Regarding the role of thalamus in human partial epilepsies, morphologic and functional neuroimaging studies have found structural alterations (19–23) and metabolic or perfusion changes (24–28) of the ipsilateral thalamus predominant in the dorsomedian nucleus.

Moreover, intracerebral recordings showed that synchronized loops between thalamus and temporal structures take place during the course of seizure and may play a role in amplifying the spreading of the discharge (29).

Based on the notion that thalamus synchronizes its activity with the cortex and on the findings of an increase of synchronization at the end of seizure, one may question the role of the thalamus in seizure termination.

Thus, the purpose of this work was to study the role of the thalamus in seizure termination through its involvement in thalamo-cortical synchronization.

To this aim, we studied thalamic activity in human mesio-temporal lobe seizures (MTLEs) using intracerebral recordings

performed with stereoelectroencephalography (SEEG). We analyzed the influence of thalamic inputs and outputs in the last phase of seizures in order to investigate the potential role of thalamus in modulating seizure termination.

# Materials and Methods

## Patient Selection and SEEG Recordings

Twenty-two MTLEs from 10 patients undergoing stereotactic EEG recordings (SEEG) were studied. All patients had a pre-surgical evaluation including detailed clinical evaluation, routine MRI, and surface video-EEG (30). Patients were retrospectively selected for the present study according to the following criteria: (i) the epileptogenic zone is located in the mesial temporal lobe; (ii) at least one intracerebral electrode reached the thalamus. During all SEEG procedures, patients signed informed consent, and the study was approved by the Institutional Review board (IRB00003888) of INSERM (IORG0003254, FWA00005831). Intracranial EEG electrodes were implanted stereotaxically according to the Talairach's method (31). The location of each electrode contact was based on CT-scan/MRI data fusion according to previous studies (32). The location of the contacts in the thalamus was based on Talairach (33) and Schaltenbrand (34) atlases. Thalamic target was mainly the medial pulvinar group or the posterior part of the dorsomedian nucleus and corresponded to the most internal contact of the electrode passing through the superior temporal gyrus (electrode H, **Figure 1**) (29, 35). This type of electrode was implanted for clinical purpose, mainly for exploring the superior temporal gyrus and the insula. Signals were recorded on a 128 channel Deltamed™ system, sampled at 256 Hz with no digital filter. Two hardware filters are a high-pass filter (cut-off frequency equal to 0.16 Hz at -3 dB) and a first-order low-pass filter (cut-off frequency equal to 97 Hz at -3 dB) to avoid aliasing.

Patients' clinical data are provided in **Table 1**.

## Procedure: SEEG Signal Analysis

### Definition of Regions

We analyzed the statistical interdependencies between bipolar signals exploring different regions of interest. We selected the contacts exploring in each patient three mesial temporal regions (amygdala, hippocampus, and entorhinal cortex, respectively A, Hip, and EC), three neocortical temporal regions (NC, corresponding to the superior, middle, or inferior temporal gyrus, respectively STG, MTG, and ITG) and the thalamus (Th) (**Figures 1A,B**). Electrodes with significant artifacts were excluded.

### Periods of Interest

For each seizure, two periods of interest (both lasting 15 s) were selected (**Figure 1C**):

a. *Seizure onset (SO)* including 5 s before and 10 s after the appearance of a tonic discharge in mesial temporal structures.

b. *End of seizure (ES)* including the last 15 s of the seizure discharge.
A

<figure><img src="image_3.png" /></figure>

B

<figure><img src="image_4.png" /></figure>

C

<figure><img src="image_5.png" /></figure>

D

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

**FIGURE 1** | Intracerebral recordings using multiple contacts electrodes placed according to Talairach's stereotactic method. (A) Schematic diagram of SEEG electrodes placement on a lateral view of a three-dimensional (3D) reconstruction of the neocortical surface in a patient with MTLE. Electrode A explored the amygdala (internal leads) and the anterior part of the MTG (lateral leads); electrode B explored the anterior hippocampus (internal leads) and the middle part of MTG (lateral leads); electrode TB explored the entorhinal cortex (internal leads) and the anterior part of ITG (external leads); electrode C explored the posterior hippocampus (internal leads) and the posterior part of MTG (lateral leads); electrode H explored the thalamus (internal leads) and the posterior part of STG (lateral leads). (B) MRI axial view of standard orthogonal

electrode H exploring thalamus by its internal leads and STG by its external leads. (C) Ictal SEEG recordings with selected traces from mesio-temporal structures, i.e., amygdala (A), hippocampus (Hip), entorhinal cortex (EC), and from thalamus (Th). Seizure onset period (SO) included 5 s before and 10 s after the onset of high-frequency activity in mesial temporal structures (top panel); the end of seizure period (ES) was defined as the last 15 s of the seizure discharge (bottom panel). (D) SEEG recordings of two seizures in the ES period displaying the thalamic rhythmic pattern with a spike-wave activity (on the top) or an arched morphology (on the bottom) coincidental with cortical spike-and-wave discharges. A, amygdala; Hip, hippocampus; EC, entorhinal cortex; NC, neocortex; Th, thalamus.

Duration of each seizure was calculated in seconds starting from the beginning of SO to the end of ES.

## "Thalamic" Patterns

Firstly, intracerebral signals were visually inspected. In the ES period, 6 out of 22 seizures displayed an electrical pattern characterized by rhythmic activity with a spike-wave or an

arched morphology coincidental with a cortical spike-and-wave discharges (SWDs) (Figure 1D). This pattern was also frequently characterized by a gradual increase in the inter-spike duration (4/6 seizures). A similar thalamic pattern has been previously described by Meeren and colleagues (4) during the generalized SWDs in a genetic model of absence epilepsy, the WAG/Rij rat. We defined this spike-and-wave pattern (SWP) as pattern "A". The
<table><caption>TABLE 1 | Clinical data of patients included in the study.</caption><thead><tr><th>Patients</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th></tr></thead><tbody><tr><td>Age (years)</td><td>33</td><td>32</td><td>22</td><td>14</td><td>36</td><td>35</td><td>17</td><td>50</td><td>54</td><td>36</td></tr><tr><td>Gender</td><td>F</td><td>M</td><td>M</td><td>F</td><td>F</td><td>F</td><td>M</td><td>F</td><td>M</td><td>M</td></tr><tr><td>Age at onset (years)</td><td>3</td><td>17</td><td>12</td><td>3</td><td>16</td><td>3</td><td>5</td><td>15</td><td>35</td><td>18</td></tr><tr><td>Epilepsy duration (years)</td><td>30</td><td>15</td><td>10</td><td>11</td><td>20</td><td>32</td><td>12</td><td>35</td><td>19</td><td>18</td></tr><tr><td>Etiology</td><td>Unk</td><td>CD</td><td>As</td><td>HS</td><td>HS</td><td>CD</td><td>HS</td><td>Non-specific</td><td>Unk</td><td>Unk</td></tr><tr><td>Aura</td><td>Deja vu</td><td>Anxiety, olfactive</td><td>Deja vu</td><td>No</td><td>No</td><td>Epigastric, cold</td><td>No</td><td>Cold, fear, taste</td><td>No</td><td>Anxiety</td></tr><tr><td>Early loss of contact</td><td>Yes</td><td>Yes</td><td>Yes</td><td>Yes</td><td>No</td><td>No</td><td>Yes</td><td>Yes</td><td>Yes</td><td>Yes</td></tr><tr><td>Secondary generalization</td><td>Rare</td><td>Frequent</td><td>Rare</td><td>No</td><td>Rare</td><td>Rare</td><td>Rare</td><td>No</td><td>Frequent</td><td>Rare</td></tr><tr><td>Type of TLE</td><td>L-MTLE</td><td>L-MTLE</td><td>R-MTLE</td><td>L-MTLE</td><td>R-MTLE</td><td>R-MTLE</td><td>L-MTLE</td><td>R-MTLE</td><td>L-MTLE</td><td>L-MTLE</td></tr></tbody></table>

M, male; F, female; As, astrocytoma affecting the lateral temporal cortex; CD, mesial temporal cortical dysplasia; HS, hippocampal sclerosis; MTL, mesial temporal lobe; Unk, unknown; TLE, temporal lobe epilepsy; L-MTLE, left mesial temporal lobe epilepsy; R-MTLE, right mesial temporal lobe epilepsy.

other 16 seizures with different electrical patterns from pattern A were named pattern “B”.

## Estimation of Signal Correlation

Non-linear regression analysis (estimation of $h^2$ coefficient) (36) was used to study functional connectivity between thalamus and selected cortical regions (A, Hip, EC, and NC). This method gives an estimation of the degree and the direction of the link between brain regions. Main practical applications of this method have concerned epileptic activity analysis both in animals and in humans (2, 4, 37, 38). This method is a pair-wise approach, where the amplitude of a signal Y is described as a function of the amplitude of a signal X using a non-linear regression curve and to compute the variance of Y that is explained, or predicted, by X according to this regression curve. A numerical approximation of the non-linear regression curve is obtained by describing the scatterplot of Y vs. X by segments of linear regression curves (i.e., this is a piecewise linear model).

The non-linear correlation coefficient $h^2$ is computed and takes its values in [0,1]. Low values of $h^2$ denote that signals X and Y are independent. On the opposite, high values of $h^2$ mean that signal Y may be explained by a continuous transformation (possibly non-linear) of signal X, i.e., signals X and Y are dependent.

Time delay between signals X and Y was also calculated. Signals were sampled at 256 Hz and a 2 s window sliding with steps of 1 s were chosen. This resulted in a connectivity graph, with each node corresponding to a region of interest. Between two nodes, we selected the higher $h^2$ value between X → Y and Y → X, and the sign of the corresponding time delay defined the direction of the connection.

In this work, we propose to use a network measure that estimates the number of output and input connections for each node of the graph (in graph theory terminology, these are the output and input degrees, respectively). Output and input degrees were computed for each region compared to the others combining the number of functional connections with $h^2$ values higher than 0.4 (in order to take into account only meaningful connections) and the time delay between signals. A visual representation of output and input connectivity was displayed in colored graphs showing

the variations of directionality measures during the course of seizures (Figures 2A,B). To the aim of this work, we specifically focused on the thalamic and mesio-temporal outputs and inputs; thalamic output and input measures refer to functional connections between the thalamus and other cortical regions explored; similarly, mesio-temporal output and input measures refer to functional connections between a mesio-temporal region and other regions explored, included thalamus. Computations were done with the Anywave freeware (39) available at http:meg.univ-amu.fr.

## Statistical Analysis

### Synchronization Indices

In order to obtain a global connectivity measure, an average of $h^2$ values of all pairs of signals was computed for each analysis window (2 s sliding window with steps of 1 s) on regions of interest. Thus, we computed a thalamic connectivity measure by averaging the $h^2$ values obtained from the interactions between the thalamus and cortical regions. A mesio-temporal connectivity measure was obtained by averaging the $h^2$ values between mesio-temporal structures and the remaining regions (including thalamus). Finally, a thalamo-mesio temporal connectivity measure was derived by averaging the $h^2$ values obtained from the interaction between thalamus and mesio-temporal regions (Figure 2C). All connectivity measures were estimated as a function of time during seizures as shown in Figure 2D.

We then performed a second average across time either in the seizure onset period (SO) and in the end of seizure period (ES). For each period, we computed the average across a 3 s window around the peak of connectivity, resulting in a global index, a thalamic index (Th index), a mesio-temporal index (MT index), and a thalamo-mesio temporal index (Th-MT index).

## Graph Measures

Similarly to the synchronization indices, we computed the average of thalamo and mesio-temporal output and input measures in the ES period. This resulted in graph measures in the ES period, specifically thalamic output and input measures (Th OUT, Th IN) and mesio-temporal output and input measures (MT OUT, MT IN).
<figure><img src="image_8.png" /></figure>

B

<figure><img src="image_9.png" /></figure>

D

<figure><img src="image_10.png" /></figure>

C

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

**FIGURE 2** | (A) The graph shows the direction, represented by the arrows, of functional connections between the regions of interest selected (on the left) in a 2 s window of an ictal SEEG recording (on the right). The strength of connections is indicated by the numerical values ($h^2$ values) scaling from 0 to 1 and by a colored scale (red representing the maximum degree and blue the minimum). (B) The colored graph displayed the variations of output and input connection measures for each region compared to the others during the end of seizure period (ES) in one MTLE seizure. A'2–A'3: internal contacts of the electrode exploring the left amygdala; A'10–A'11: external contacts of the same electrode exploring the anterior part of the left MTG; B'1–B'2: internal contacts of the electrode exploring the left anterior hippocampus; B'12–B'13: external contacts of the same electrode exploring the middle part of the left MTG; TB'1–TB'2: internal contacts of the

electrode exploring the left entorhinal cortex; TB'11–TB'12: external contacts of the same electrode exploring the anterior part of the left ITG; T'1–T'2: internal contacts of the electrode exploring left thalamus. (C) Connectivity measures on regions of interest: global connectivity measure (yellow arrow), thalamic connectivity measure (red arrow), mesio-temporal connectivity measure (green arrow) and thalamo-mesio temporal connectivity measure (blue arrow). Th, thalamus; MT, mesio-temporal regions; NC, neocortical regions. (D) Variations in connectivity measure profiles during the course of one patient's seizure. Arrows pointed respectively to the onset and the end of seizure. Global $h^2$, global connectivity measure; th $h^2$, thalamic connectivity measure; MT $h^2$, mesio-temporal connectivity measure; th-MT $h^2$, thalamo-mesio temporal connectivity measure; ES, end of seizure; SO, seizure onset.

## Statistical Tests

In order to analyze how thalamus, mesio-temporal regions, and neocortex synchronize their activity during the course of seizures, global, thalamo, mesio-temporal, and thalamo-mesio temporal synchronization indices (global index, Th index, MT index, and Th-MT index) were compared between the SO and ES period by a Wilcoxon test.

A non-parametric correlation test (Spearman test) was performed to find possible correlations between synchronization indices in ES period (global index, Th index, MT index, and Th-MT index) and seizure duration.

In order to compare seizure duration, synchronization, and direction indices between the two groups (pattern A and B), a non-parametric Mann-Whitney test has been applied.

Electrophysiological data are provided in Table 2.

## Results

### Dynamic Changes According to Synchronization Indices

Significant differences of all synchronization indices were observed between the onset of seizure (SO) and the end of seizure (ES). Indeed, global, thalamo, mesio-temporal, and thalamo-mesio temporal synchronization indices (global index, Th index, MT index, and Th-MT index) displayed a marked increase in the ES period compared to the SO period ($p < 0.0001$ for each comparison, **Figure 3A**).
<table><caption>TABLE 2 | Electrophysiological data.</caption><thead><tr><th>Patients</th><th>Seizures</th><th>Pattern</th><th>Seizure duration(s)</th><th>ES global index</th><th>ES Th index</th><th>ES MT index</th><th>ES Th-MT index</th><th>SO global index</th><th>SO Th index</th><th>SO MT index</th><th>SO Th-MT index</th><th>Th OUT</th><th>Th IN</th><th>MT OUT</th><th>MT IN</th></tr></thead><tbody><tr><td>1</td><td>1</td><td>A</td><td>55</td><td>0.49</td><td>0.51</td><td>0.52</td><td>0.82</td><td>0.19</td><td>0.22</td><td>0.24</td><td>0.51</td><td>4.25</td><td>0.38</td><td>3.13</td><td>1.5</td></tr><tr><td>2</td><td>2</td><td>A</td><td>90</td><td>0.48</td><td>0.5</td><td>0.53</td><td>0.81</td><td>0.22</td><td>0.25</td><td>0.18</td><td>0.20</td><td>2.17</td><td>1.00</td><td>3.67</td><td>1.00</td></tr><tr><td>3</td><td>3</td><td>A</td><td>150</td><td>0.46</td><td>0.52</td><td>0.42</td><td>0.47</td><td>0.17</td><td>0.12</td><td>0.17</td><td>0.14</td><td>3.30</td><td>2.20</td><td>1.00</td><td>3.50</td></tr><tr><td>4</td><td>4</td><td>B</td><td>540</td><td>0.26</td><td>0.21</td><td>0.24</td><td>0.2</td><td>0.15</td><td>0.10</td><td>0.22</td><td>0.14</td><td>1.18</td><td>0.73</td><td>1.55</td><td>0.18</td></tr><tr><td>5</td><td>5</td><td>B</td><td>59</td><td>0.38</td><td>0.41</td><td>0.34</td><td>0.32</td><td>0.26</td><td>0.22</td><td>0.30</td><td>0.21</td><td>0.50</td><td>0.40</td><td>1.60</td><td>0.60</td></tr><tr><td>6</td><td>6</td><td>B</td><td>63</td><td>0.33</td><td>0.36</td><td>0.33</td><td>0.42</td><td>0.19</td><td>0.16</td><td>0.24</td><td>0.13</td><td>0.56</td><td>1.22</td><td>0.56</td><td>0.44</td></tr><tr><td>7</td><td>7</td><td>B</td><td>88</td><td>0.44</td><td>0.45</td><td>0.44</td><td>0.36</td><td>0.22</td><td>0.17</td><td>0.26</td><td>0.18</td><td>1.90</td><td>1.80</td><td>3.00</td><td>1.70</td></tr><tr><td>8</td><td>8</td><td>A</td><td>60</td><td>0.59</td><td>0.55</td><td>0.59</td><td>0.6</td><td>0.19</td><td>0.12</td><td>0.15</td><td>0.10</td><td>2.91</td><td>1.64</td><td>2.55</td><td>1.82</td></tr><tr><td>9</td><td>9</td><td>A</td><td>53</td><td>0.49</td><td>0.53</td><td>0.5</td><td>0.53</td><td>0.22</td><td>0.14</td><td>0.19</td><td>0.13</td><td>0.67</td><td>3.17</td><td>2.67</td><td>0.08</td></tr><tr><td>10</td><td>10</td><td>A</td><td>70</td><td>0.49</td><td>0.58</td><td>0.42</td><td>0.47</td><td>0.14</td><td>0.12</td><td>0.10</td><td>0.13</td><td>2.00</td><td>1.86</td><td>1.29</td><td>2.00</td></tr><tr><td>11</td><td>11</td><td>B</td><td>126</td><td>0.41</td><td>0.43</td><td>0.4</td><td>0.52</td><td>0.13</td><td>0.14</td><td>0.11</td><td>0.12</td><td>1.75</td><td>1.75</td><td>2.63</td><td>1.25</td></tr><tr><td>12</td><td>12</td><td>B</td><td>54</td><td>0.36</td><td>0.33</td><td>0.32</td><td>0.29</td><td>0.15</td><td>0.10</td><td>0.12</td><td>0.09</td><td>0.8</td><td>0.70</td><td>2.10</td><td>0.80</td></tr><tr><td>13</td><td>13</td><td>B</td><td>65</td><td>0.35</td><td>0.37</td><td>0.43</td><td>0.45</td><td>0.17</td><td>0.11</td><td>0.17</td><td>0.12</td><td>0.22</td><td>1.78</td><td>3.33</td><td>0.22</td></tr><tr><td>14</td><td>14</td><td>B</td><td>200</td><td>0.33</td><td>0.33</td><td>0.38</td><td>0.57</td><td>0.13</td><td>0.12</td><td>0.13</td><td>0.11</td><td>1.50</td><td>1.70</td><td>2.40</td><td>1.00</td></tr><tr><td>15</td><td>15</td><td>B</td><td>290</td><td>0.29</td><td>0.29</td><td>0.32</td><td>0.4</td><td>0.17</td><td>0.15</td><td>0.16</td><td>0.18</td><td>0.5</td><td>0.50</td><td>0.17</td><td>1.33</td></tr><tr><td>16</td><td>16</td><td>B</td><td>455</td><td>0.27</td><td>0.27</td><td>0.34</td><td>0.25</td><td>0.11</td><td>0.12</td><td>0.13</td><td>0.29</td><td>2.20</td><td>0.20</td><td>2.20</td><td>0.60</td></tr><tr><td>17</td><td>17</td><td>B</td><td>70</td><td>0.27</td><td>0.34</td><td>0.39</td><td>0.42</td><td>0.14</td><td>0.18</td><td>0.16</td><td>0.47</td><td>0.25</td><td>1.25</td><td>1.63</td><td>0.38</td></tr><tr><td>18</td><td>18</td><td>B</td><td>113</td><td>0.31</td><td>0.33</td><td>0.34</td><td>0.55</td><td>0.11</td><td>0.14</td><td>0.10</td><td>0.14</td><td>1.11</td><td>0.78</td><td>1.22</td><td>0.78</td></tr><tr><td>19</td><td>19</td><td>B</td><td>124</td><td>0.31</td><td>0.42</td><td>0.35</td><td>0.54</td><td>0.16</td><td>0.17</td><td>0.16</td><td>0.19</td><td>1.89</td><td>0.89</td><td>1.22</td><td>1.78</td></tr><tr><td>20</td><td>20</td><td>B</td><td>64</td><td>0.46</td><td>0.38</td><td>0.42</td><td>0.46</td><td>0.37</td><td>0.28</td><td>0.45</td><td>0.37</td><td>0.60</td><td>1.20</td><td>1.30</td><td>0.80</td></tr><tr><td>21</td><td>21</td><td>B</td><td>50</td><td>0.37</td><td>0.43</td><td>0.39</td><td>0.38</td><td>0.32</td><td>0.26</td><td>0.42</td><td>0.35</td><td>0.70</td><td>2.50</td><td>2.30</td><td>0.20</td></tr><tr><td>22</td><td>22</td><td>B</td><td>55</td><td>0.28</td><td>0.38</td><td>0.31</td><td>0.37</td><td>0.22</td><td>0.17</td><td>0.20</td><td>0.12</td><td>1.00</td><td>1.60</td><td>1.70</td><td>0.30</td></tr></tbody></table>

Global index, global synchronization index; Th index, thalamic synchronization index; MT index, mesio-temporal synchronization index; Th-MT index, thalamo-mesio temporal synchronization index; ES, end of seizure period; SO, seizure onset period; Th OUT, thalamic output measure; Th IN, thalamic input index; MT OUT, mesio-temporal output measure; MT IN, mesio-temporal input index.

In order to estimate the influence of synchronization changes on seizure duration, we tested for a correlation between synchronization indices at ES period (global index, Th index, MT index, and Th-MT index) and duration of seizures.

The global index was negatively correlated with seizure duration ($p = 0.045$), namely higher synchronization values corresponded to lower seizure duration. In the same way, we found a trend for an inverse tendency between Th index and seizure duration ($p = 0.052, R^2 = 0.40$) (Figure 3B). In contrast, no significant correlations or tendency were found between other synchronization indices (Th-MT index and MT index) and seizure duration.

## Pattern Differences

We then studied how connectivity measures, namely global, thalamic, mesio-temporal, and thalamo-mesio temporal, vary as a function of time during the last phase of seizure. Comparing group A (seizures with thalamic spike-and-wave pattern) with group B (seizures with other thalamic patterns), thalamic connectivity disclosed different courses in the ES period. Indeed, pattern A seizures showed an overall stable, higher, and longer increase of thalamo-cortical synchronization clearly before the end of seizure. On the contrary, pattern B seizures showed a heterogeneous, unstable, and shorter increase of thalamic synchronization, often localized at the beginning of ES period and followed by a slow synchrony decrease (Figure 3C).

We compared the synchronization indices and the direction indices derived from ES period between pattern A and pattern B seizures. All synchronization indices were significantly higher

in pattern A than pattern B group (global index $p = 0.0004$, Th index $p = 0.0004$, Th-MT index $p = 0.0063$, and MT index $p = 0.0015$) (Figure 4A).

Concerning the direction measures, only thalamic output measure (Th OUT) showed significant higher values ($p = 0.0079$) in pattern A compared with pattern B group (Figure 4B). On the contrary, no significant differences were found comparing thalamic input, mesio-temporal output, and mesio-temporal input measures (Th IN, MT OUT, and MT IN). Therefore, in the last phase of seizure, synchronization levels and thalamic outputs were higher in pattern A than in pattern B (Figure 4C).

Taken as a whole, no statistical significant difference was found in seizure duration between A and B groups but pattern A group included only short seizures, while pattern B comprised a wide range of seizure duration (Figure 4D). Moreover, in two patients showing seizures with both patterns, pattern A seizures were shorter than pattern B seizures (150 s vs. 540 s, 70 s vs. 126 s).

## Discussion

In the present study, we analyzed the ictal functional connectivity between mesio-temporal regions and thalamus by non-linear correlation analysis in 22 drug-resistant MTLEs. The aim of this work was to explore the dynamic involvement of thalamus in partial seizures and its potential contribution in seizure termination.

In all the 22 analyzed mesio-temporal seizures, thalamic activity was affected by the epileptic discharge during the course of seizure. Results of non-linear correlation analysis showed a
A

<figure><img src="image_13.png" /></figure>

B

<figure><img src="image_14.png" /></figure>

C

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

**FIGURE 3** | (A) Increase of the thalamic synchronization index (Th index) in the ES period compared to the SO period in all seizures analyzed (*p* < 0.0001 for each comparison) by Wilcoxon test. (B) Negative correlation between global index and seizure duration (*p* = 0.045) and inverse tendency between Th index and seizure duration (*p* = 0.052, R² 0.40) in the ES period by Spearman test.

(C) Different patterns of the increase of thalamic connectivity between groups A and B in the ES period; pattern A seizures (N<sub>A</sub> = 6) showed an overall stable, higher, and longer increase of thalamo-cortical synchronization clearly before the end of seizure compared to pattern B seizures (N<sub>B</sub> = 16). The mean and the average deviation of values are shown. Dashed arrow indicates the end of seizure.

significant increase of synchronization between thalamus, mesio-temporal region and neocortex (*p* < 0.0001) in the last phase of seizure, compared to the beginning.

These results are in agreement with previous animal and human studies that reported an increase of synchronization at the end of seizure. Indeed, in an extensive series of *in vivo* cat experiments, including multi-site intracellular and extracellular as well as local cortical field potential recordings (5, 8, 40), it has been observed a progressive increase of cortical synchrony during the course of seizure, thus promoting the recruitment of neurons. However, when all the affected neuronal pool was involved into highly synchronous paroxysmal activity, the seizure stopped. As a possible mechanism, the authors proposed that the increased depolarization during the full-blown seizure was associated with the activation of the hyperpolarizing sodium- and calcium-activated potassium currents, overcoming the depolarizing influence of ionic currents. Moreover, an additional contribution may be represented by a dramatic decrease in the extracellular calcium

concentration resulting in the decrease of effectiveness of synaptic transmission.

Recent intracranial studies on human partial seizures, making use of different signal interdependencies analysis, described a progressive increase of cortical correlation, maximal at the end of seizure (9, 10, 29, 41).

Although the exact relationship between increase of synchronization and seizure termination is not yet completely understood, the highly synchronous paroxysmal activity in large neuronal networks observed at the end of seizure may be considered as an active seizure termination mechanism.

The role of thalamo-cortical interactions in modulating epi-leptic activity has been well demonstrated in absence epilepsy. Indeed, highly synchronized thalamo-cortical oscillations repre-sent the neuronal substrate for the generation and maintenance of SWDs and thalamus probably provides a resonant circuitry to amplify and sustain the rhythmic discharges (4).

Although the thalamic involvement in partial seizures has been described in animal models (13–16, 42), its role in modulating
<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

**FIGURE 4** | Comparison between group A and B. (A) Mean and standard deviation of synchronization indices derived from ES period for the two groups; all synchronization indices were significantly higher in pattern A than pattern B seizures (Mann-Whitney test, global index *p* = 0.0004, Th index *p* = 0.0004, Th-MT index *p* = 0.0063, and MT index *p* = 0.0015). (B) Mean and standard deviation of thalamic output graph measures (Th OUT) in group A and B

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /></figure>

(*p* = 0.0079). (C) Pattern A seizures showed higher thalamic synchronization indices (Th index) and thalamic output measures (Th OUT) than pattern B seizures in the ES period. (D) Mean of seizure duration in group A and group B. No statistical difference between seizure A and seizure B duration, but pattern A group had only short seizures, while pattern B included a wide range of seizure duration.

the epileptic activity, namely influencing generation, propagation, and termination of seizures, has been less characterized in humans.

Indeed, only few intracerebral studies explored thalamic activity in human partial seizures (29, 35, 43, 44). We previously investigated the thalamic activity during temporal lobe seizures using a non-linear correlation analysis (29). We demonstrated that thalamo-cortical loop takes place during the course of mesiotemporal seizures after a variable latency and tends to be more strongly involved between thalamus and neocortical regions, while seizures started from mesiotemporal regions.

The potential role of thalamus in seizure termination is however not well known. Taken as a whole, subcortical networks, including the thalamo-cortical loop, may modulate different dynamic aspects of epileptic seizure (7, 45).

In this study, we found a constant thalamic recruitment and a marked increase of thalamo-cortical synchronization in the last phase of each analyzed seizure, thus confirming the concept that thalamo-cortical network is highly involved in modulating the epileptic discharge in human mesiotemporal lobe epilepsy.

Moreover, we found an inverse tendency between seizure duration and respectively global and thalamo-cortical synchronization levels at the end of seizure. Therefore, an increase of final global and thalamo-cortical synchronization corresponded

to a decrease in seizure duration. We speculated that thalamo-cortical oscillations might exert a seizure terminating role via modulation of cortical synchronization. In fact, thalamo-cortical oscillations promote a progressive increase of hypersynchronous activity in large neuronal networks during the course of seizure, thus likely leading to a massive and simultaneous neuronal hyperpolarization and finally to seizure termination [see Ref. (5)].

Recently, Langlois and colleagues examined the involvement of the thalamic parafascicular nucleus (PF) in the occurrence of hippocampal paroxysmal discharges (HPDs), in *in vivo* intracellular and extracellular as well as local field potential recordings in a chronic mouse model of MTLE. They showed that PF neurons fire synchronously just before the end of the HPDs (18). Moreover, authors demonstrated that the high-frequency PF stimulation (130 Hz) interrupted the ongoing HPS and the pharmacological inhibition or activation of the PF neurons respectively suppresses or aggravates HPDs, suggesting that the PF nucleus plays a role in the active termination of seizure.

The medial pulvinar group (PuM) is the main explored part of the thalamus in our study. Recording this structure does not require electrodes in addition to those exploring regions potentially involved in the epileptogenic zone. Recent intracerebral-evoked responses' recordings obtained after PuM and cortical electrical stimulation in seven epileptic patients undergoing pre-surgical evaluation, confirmed the developed PuM connectivity
with mesial temporal structures (44), as previously reported by anatomical studies in monkeys (46, 47).

Moreover, in human thalamus, PuM is one of the most concentrated nucleus in calbindin matrix cells (48, 49). These thalamo-cortical matrix cells receive inputs from subcortical pathways and project more diffusely upon the cortex. In turn, cortico-thalamic fibers project to new thalamic nuclei and via matrix cells to other cortical areas, supporting the spread of activity (50, 51). Such an organization of connectivity between PuM and cortex would allow the modulation of cortico-cortical networks (52).

By visual inspection of intracerebral ictal activity at the end of seizures, a definite subgroup of seizures displaying a spike-and-wave pattern has been identified (SWP-pattern A seizures). A similar pattern has been previously described in regard to generalized SWDs in animal model of absence epilepsy and reflects the highly synchronized oscillations of thalamo-cortical networks (4). In the last phase of spike-and-wave seizures, we found much higher overall synchronization values, especially for the global and thalamo-cortical synchronization (p = 0.0004) compared to seizures with different patterns (pattern B seizures). Moreover, this group displayed a higher and sustained increase of thalamo-cortical synchronization clearly before seizure termination with a stronger participation of thalamic outputs. In addition, we found that this pattern was clearly related to a significant higher value of thalamic output measure (Th OUT), indicating that the driving input in the loop is thalamus.

Although there was not a significant difference in duration between the two groups, the spike-and-wave seizures were always short in duration, while the other pattern seizures varied in a

wider range of durations. Furthermore, in patients showing both seizure patterns, spike-and-wave pattern seizures were consistently shorter.

Altogether, these results suggest that MTLEs disclosing this spike-and-wave pattern have stronger thalamo-cortical synchronized activity before the end of seizure, driven by thalamic outputs. Moreover, we speculate that thalamus may in this case exert a control on cortical structures by inducing a stable hypersynchronization that ultimately leads to seizure termination.

Finally, a better understanding of the role of subcortical circuits in partial seizure termination may improve insights into brain autoregulatory mechanisms. Moreover, it may offer new venues for developing novel treatments and for advancing protocols of thalamic deep brain stimulation (DBS) in patients with epilepsy. Recent development of DBS in epilepsy mainly focused on the anterior nucleus (53), showing an efficacy for reducing seizures in a significant percentage of patients. Our data suggest that the PuM could be an interesting target in this perspective.

## Author Contributions

EE and CB contributed to study design, data collection, data analysis, and writing manuscript. FBo and RC contributed to data collection and writing manuscript. RC and JR performed stereotactic explorations of patients. BC has contributed to software development and signal analysis. FBa has contributed to study design, data collection and analysis, and writing manuscript. All authors had reviewed and agreed the final manuscript.

## References

1. Bertram EH, Zhang DX, Mangan P, Fountain N, Rempe D. Functional anatomy of limbic epilepsy: a proposal for central synchronization of a diffusely hyperexcitable network. *Epilepsy Res* (1998) **32**:194-205. doi:10.1016/S0920-1211(98)00051-5

2. Bartolomei F, Wendling F, Bellanger JJ, Regis J, Chauvel P. Neural networks involving the medial temporal structures in temporal lobe epilepsy. *Clin Neurophysiol* (2001) **112**:1746-60. doi:10.1016/S1388-2457(01)00591-0

3. Wendling F, Bartolomei F. Modeling EEG signals and interpreting measures of relationship during temporal-lobe seizures: an approach to the study of epileptogenic networks. *Epileptic Disord* (2001) *Special Issue*:67-78.

4. Meeren HK, Pijn JP, Van Luijtelaar EL, Coenen AM, Lopes Da Silva FH. Cortical focus drives widespread corticothalamic networks during spontaneous absence seizures in rats. *J Neurosci* (2002) **22**:1480-95.

5. Timofeev I, Steriade M. Neocortical seizures: initiation, development and cessation. *Neuroscience* (2004) **123**:299-336. doi:10.1016/j.neuroscience.2003.08.051

6. Huguenard JR, McCormick DA. Thalamic synchrony and dynamic regulation of global forebrain oscillations. *Trends Neurosci* (2007) **30**:350-6. doi:10.1016/j.tins.2007.05.007

7. Lado FA, Moshe SL. How do seizures stop? *Epilepsia* (2008) **49**:1651-64. doi:10.1111/j.1528-1167.2008.01669.x

8. Topolnik L, Steriade M, Timofeev I. Partial cortical deafferentation promotes development of paroxysmal activity. *Cereb Cortex* (2003) **13**:883-93. doi:10.1093/cercor/13.8.883

9. Ponten SC, Bartolomei F, Stam CJ. Small-world networks and epilepsy: graph theoretical analysis of intracerebrally recorded mesial temporal lobe seizures. *Clin Neurophysiol* (2007) **118**:918-27. doi:10.1016/j.clinphys.2006.12.002

10. Schindler K, Leung H, Elger CE, Lehnertz K. Assessing seizure dynamics by analysing the correlation structure of multichannel intracranial EEG. *Brain* (2007) **130**:65-77. doi:10.1093/brain/awl304

11. McCormick DA, Contreras D. On the cellular and network bases of epileptic seizures. *Annu Rev Physiol* (2001) **63**:815-46. doi:10.1146/annurev.physiol.63.1.815

12. Turski WA, Cavalheiro EA, Bortolotto ZA, Mello LM, Schwarz M, Turski L. Seizures produced by pilocarpine in mice: a behavioral, electroencephalographic and morphological analysis. *Brain Res* (1984) **321**:237-53. doi:10.1016/0006-8993(84)90177-X

13. Clifford DB, Olney JW, Maniotis A, Collins RC, Zorumski CF. The functional anatomy and pathology of lithium-pilocarpine and high-dose pilocarpine seizures. *Neuroscience* (1987) **23**:953-68. doi:10.1016/0306-4522(87)90171-0

14. Cassidy RM, Gale K. Mediodorsal thalamus plays a critical role in the development of limbic motor seizures. *J Neurosci* (1998) **18**:9002-9.

15. Bertram EH, Scott C. The pathological substrate of limbic epilepsy: neuronal loss in the medial dorsal thalamic nucleus as the consistent change. *Epilepsia* (2000) **41**(Suppl 6):S3-8. doi:10.1111/j.1528-1157.2000.tb01548.x

16. Bertram EH, Mangan PS, Zhang D, Scott CA, Williamson JM. The midline thalamus: alterations and a potential role in limbic epilepsy. *Epilepsia* (2001) **42**:967-78. doi:10.1046/j.1528-1157.2001.042008967.x

17. Bertram EH, Zhang D, Williamson JM. Multiple roles of midline dorsal thalamic nuclei in induction and spread of limbic seizures. *Epilepsia* (2008) **49**:256-68. doi:10.1111/j.1528-1167.2007.01408.x

18. Langlois M, Polack PO, Bernard H, David O, Charpier S, Depaulis A, et al. Involvement of the thalamic parafascicular nucleus in mesial temporal lobe epilepsy. *J Neurosci* (2010) **30**:16523-35. doi:10.1523/JNEUROSCI.1109-10.2010

19. DeCarli C, Hatta J, Fazilat S, Gaillard WD, Theodore WH. Extratemporal atrophy in patients with complex partial seizures of left temporal origin. *Ann Neurol* (1998) **43**:41-5. doi:10.1002/ana.410430110

20. Dreifuss S, Vingerhoets FJ, Lazeyras F, Andino SG, Spinelli L, Delavelle J, et al. Volumetric measurements of subcortical nuclei in patients with temporal lobe epilepsy. *Neurology* (2001) **57**:1636-41. doi:10.1212/WNL.57.9.1636
21. Bernasconi A, Bernasconi N, Natsume J, Antel SB, Andermann F, Arnold DL. Magnetic resonance spectroscopy and imaging of the thalamus in idiopathic generalized epilepsy. *Brain* (2003) **126**:2447-54. doi:10.1093/brain/awg249

22. Natsume J, Bernasconi N, Andermann F, Bernasconi A. MRI volumetry of the thalamus in temporal, extratemporal, and idiopathic generalized epilepsy. *Neurology* (2003) **60**:1296-300. doi:10.1212/01.WNL.0000058764.34968.C2

23. Bonilha L, Rorden C, Castellano G, Pereira F, Rio PA, Cendes F, et al. Voxel-based morphometry reveals gray matter network atrophy in refractory medial temporal lobe epilepsy. *Arch Neurol* (2004) **61**:1379-84. doi:10.1001/archneurol.61.9.1379

24. Khan N, Leenders KL, Hajek M, Maguire P, Missimer J, Wieser HG. Thalamic glucose metabolism in temporal lobe epilepsy measured with 18F-FDG positron emission tomography (PET). *Epilepsy Res* (1997) **28**:233-43. doi:10.1016/S0920-1211(97)00049-1

25. Juhasz C, Nagy F, Watson C, Da Silva EA, Muzik O, Chugani DC, et al. Glucose and [11C]flumazenil positron emission tomography abnormalities of thalamic nuclei in temporal lobe epilepsy. *Neurology* (1999) **53**:2037-45. doi:10.1212/WNL.53.9.2037

26. Newberg AB, Alavi A, Berlin J, Mozley PD, O'Connor M, Sperling M. Ipsilateral and contralateral thalamic hypometabolism as a predictor of outcome after temporal lobectomy for seizures. *J Nucl Med* (2000) **41**:1964-8.

27. Benedek K, Juhasz C, Muzik O, Chugani DC, Chugani HT. Metabolic changes of subcortical structures in intractable focal epilepsy. *Epilepsia* (2004) **45**:1100-5. doi:10.1111/j.0013-9580.2004.43303.x

28. Blumenfeld H, McNally K, Vanderhill S, Paige A, Chung R, Davis K, et al. Positive and negative network correlations in temporal lobe epilepsy. *Cereb Cortex* (2004) **14**:892-902. doi:10.1093/cercor/bhh048

29. Guye M, Regis J, Tamura M, Wendling F, Mcgonigal A, Chauvel P, et al. The role of corticothalamic coupling in human temporal lobe epilepsy. *Brain* (2006) **129**:1917-28. doi:10.1093/brain/awl151

30. Maillard L, Vignal JP, Gavaret M, Guye M, Biraben A, Mcgonigal A, et al. Semiologic and electrophysiologic correlations in temporal lobe seizure subtypes. *Epilepsia* (2004) **45**:1590-9. doi:10.1111/j.0013-9580.2004.09704.x

31. Talairach J, Bancaud J, Bonis A, Szikla G, Trottier S, Vignal JP, et al. Surgical therapy for frontal epilepsies. *Adv Neurol* (1992) **57**:707-32.

32. Bartolomei F, Wendling F, Regis J, Gavaret M, Guye M, Chauvel P. Pre-ictal synchronicity in limbic networks of mesial temporal lobe epilepsy. *Epilepsy Res* (2004) **61**:89-104. doi:10.1016/j.eplepsyres.2004.06.006

33. Talairach J, Tournoux P. *Co-Planar Stereotaxic Atlas of the Human Brain*. New York, NY: Thieme (1988).

34. Schaltenbrand G, Wahren W. *Atlas for Stereotaxy of the Human Brain. With an Accompanying Guide*. 2nd ed. Stuttgart: Thieme (1977).

35. Arthuis M, Valton L, Regis J, Chauvel P, Wendling F, Naccache L, et al. Impaired consciousness during temporal lobe seizures is related to increased long-distance cortical-subcortical synchronization. *Brain* (2009) **132**:2091-101. doi:10.1093/brain/awp086

36. Pijn J, Lopes Da Silva F. Propagation of electrical activity: nonlinear associations and time delays between EEG signals. In: Zschocke S, Speckmann E-J, editors. *Basic Mechanisms of the EEG*. Boston, MA: Birkauser (1993). p. 41-61.

37. Wendling F, Bartolomei F, Bellanger JJ, Chauvel P. Interpretation of interdependencies in epileptic signals using a macroscopic physiological model of the EEG. *Clin Neurophysiol* (2001) **112**:1201-18. doi:10.1016/S1388-2457(01)00547-8

38. Chavez M, Le Van Quyen M, Navarro V, Baulac M, Martinerie J. Spatio-temporal dynamics prior to neocortical seizures: amplitude versus phase couplings. *IEEE Trans Biomed Eng* (2003) **50**:571-83. doi:10.1109/TBME.2003.810696

39. Colombet B, Woodman M, Badier JM, Benar CG. AnyWave: a cross-platform and modular software for visualizing and processing electrophysiological signals. *J Neurosci Methods* (2015) **242**:118-26. doi:10.1016/j.jneumeth.2015.01.017

40. Steriade M, Amzica F. Dynamic coupling among neocortical neurons during evoked and spontaneous spike-wave seizure activity. *J Neurophysiol* (1994) **72**:2051-69.

41. Schiff SJ, Sauer T, Kumar R, Weinstein SL. Neuronal spatiotemporal pattern discrimination: the dynamical evolution of seizures. *Neuroimage* (2005) **28**:1043-55. doi:10.1016/j.neuroimage.2005.06.059

42. Turski L, Cavalheiro EA, Schwarz M, Turski WA, De Moraes Mello LE, Bortolotto ZA, et al. Susceptibility to seizures produced by pilocarpine in rats after microinjection of isoniazid or gamma-vinyl-GABA into the substantia nigra. *Brain Res* (1986) **370**:294-309. doi:10.1016/0006-8993(86)90484-1

43. Rosenberg DS, Mauguiere F, Demarquay G, Ryvlin P, Isnard J, Fischer C, et al. Involvement of medial pulvinar thalamic nucleus in human temporal lobe seizures. *Epilepsia* (2006) **47**:98-107. doi:10.1111/j.1528-1167.2006.00375.x

44. Rosenberg DS, Mauguiere F, Catenoix H, Faillenot I, Magnin M. Reciprocal thalamocortical connectivity of the medial pulvinar: a depth stimulation and evoked potential study in human brain. *Cereb Cortex* (2009) **19**:1462-73. doi:10.1093/cercor/bhn185

45. Deransart C, Depaulis A. The control of seizures by the basal ganglia? A review of experimental data. *Epileptic Disord* (2002) **4**(Suppl 3):S61-72.

46. Insausti R, Amaral DG, Cowan WM. The entorhinal cortex of the monkey: III. Subcortical afferents. *J Comp Neurol* (1987) **264**:396-408. doi:10.1002/cne.902640307

47. Grieve KL, Acuna C, Cudeiro J. The primate pulvinar nuclei: vision and action. *Trends Neurosci* (2000) **23**:35-9. doi:10.1016/S0166-2236(99)01482-4

48. Jones EG. Viewpoint: the core and matrix of thalamic organization. *Neuroscience* (1998) **85**:331-45. doi:10.1016/S0306-4522(97)00581-2

49. Munkle MC, Waldvogel HJ, Faull RL. The distribution of calbindin, calretinin and parvalbumin immunoreactivity in the human thalamus. *J Chem Neuroanat* (2000) **19**:155-73. doi:10.1016/S0891-0618(00)00060-0

50. Jones EG. *The Thalamus*. 2nd ed. Cambridge: Cambridge University Press (2007).

51. Jones EG. Synchrony in the interconnected circuitry of the thalamus and cerebral cortex. *Ann N Y Acad Sci* (2009) **1157**:10-23. doi:10.1111/j.1749-6632.2009.04534.x

52. Shipp S. The functional logic of cortico-pulvinar connections. *Philos Trans R Soc Lond B Biol Sci* (2003) **358**:1605-24. doi:10.1098/rstb.2002.1213

53. Fisher R, Salanova V, Witt T, Worth R, Henry T, Gross R, et al. Electrical stimulation of the anterior nucleus of thalamus for treatment of refractory epilepsy. *Epilepsia* (2010) **51**:899-908. doi:10.1111/j.1528-1167.2010.02536.x

**Conflict of Interest Statement:** The 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.

Copyright © 2015 Evangelista, Bénar, Bonini, Carron, Colombet, Régis and Bartolomei. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.