# Consciousness alteration in focal epilepsy is related to loss of signal complexity and information processing

Nada El Youssef, Aude Jegou, Julia Makhalova, Lionel Naccache, Christian Bénar, Fabrice Bartolomei

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

## To cite this version:

Nada El Youssef, Aude Jegou, Julia Makhalova, Lionel Naccache, Christian Bénar, et al. Consciousness alteration in focal epilepsy is related to loss of signal complexity and information processing. Scientific Reports, 2022, 12 (1), pp.22276. 10.1038/s41598-022-25861-4. hal-04100836

HAL Id: hal-04100836

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

Submitted on 28 Sep 2023

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_2.png" /></figure>

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

# Consciousness alteration in focal epilepsy is related to loss of signal complexity and information processing

Nada El Youssef<fnref n="1" />,⁶, Aude Jegou²,⁶, Julia Makhalova<fnref n="1" />,³, Lionel Naccache⁴, Christian Bénar²,⁶ & Fabrice Bartolomei<fnref n="1" />,²,⁵,⁶✉

Alteration of awareness is a main feature of focal epileptic seizures. In this work, we studied how the information contained in EEG signals was modified during temporal lobe seizures with altered awareness by using permutation entropy (PE) as a measure of the complexity of the signal. PE estimation was performed in thirty-six seizures of sixteen patients with temporal lobe epilepsy who underwent SEEG recordings. We tested whether altered awareness (based on the Consciousness Seizure Score) was correlated with a loss of signal complexity. We estimated global changes in PE as well as regional changes to gain insight into the mechanisms associated with awareness impairment. Our results reveal a positive correlation between the decrease of entropy and the consciousness score as well as the existence of a threshold on entropy that could discriminate seizures with no alteration of awareness from seizures with profound alteration of awareness. The loss of signal complexity was diffuse, extending bilaterally and to the associative cortices, in patients with profound alteration of awareness and limited to the temporal mesial structures in patients with no alteration of awareness. Thus PE is a promising tool to discriminate between the different subgroups of awareness alteration in TLE.

Temporal lobe epilepsy (TLE) is the most common form of focal epilepsy. Amongst the clinical manifestations of TLE seizures, alteration of consciousness is one of the most dramatic features, affecting 60–80% of patients<fnref n="1" /> and impacting quality of life. This justifies research efforts towards a better understanding of the mechanisms of loss of consciousness and towards finding more effective therapies². Consciousness is a complex concept, which depends on the interaction between bilateral cortical networks and subcortical components, forming the “consciousness system”³,⁴. Alteration of consciousness during focal seizures can be conceptualized in a two-dimensional view of consciousness² by integrating two main components: awareness, subserved by the corticocortical networks⁵, and wakefulness, associated with reticulo-thalamo-cortical projections⁶. However, previous work in focal epilepsy has mainly focused on alteration of the level of awareness. In the new ILAE classification of seizure types, “consciousness” was thus replaced by “awareness”⁷, which refers to the awareness of events occurring during a seizure⁸.

Two main theories have been proposed to account for the alteration of awareness during focal seizures: the “global workspace theory” (GWT)²,⁹,¹⁰ and “the network inhibition hypothesis”³,⁴,¹¹,¹². The GWT suggests that coordinated activity from the associative cortices (mainly the prefrontal cortex and the posterior parietal associative cortex) is needed to access consciousness from different sensory modalities¹³. Alteration of awareness in TLE could be linked to the alteration of the global workspace through hypersynchronization between cortical and subcortical structures (thalamo-cortical connections)¹⁴. The degree of thalamo-cortical synchrony has been shown to correlate with the degree of consciousness changes¹⁴⁻¹⁶ and consciousness changes have been related to hypersynchrony in the beta/gamma range¹⁰,¹⁷. The “network inhibition hypothesis”³,⁴,¹¹,¹² postulates the deactivation

<footnote n="1">¹APHM, Timone Hospital, Epileptology and Cerebral Rhythmology, Marseille, France. ²Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France. ³APHM, Timone Hospital, CEMEREM, Marseille, France. ⁴APHM, Departments of Neurology & Clinical Neurophysiology Pitié Salpêtrière Hospital, Paris, France. ⁵Service d'Epileptologie et de Rythmologie Cérébrale, Hôpital Timone, 264 Rue Saint-Pierre, 13005 Marseille, France. ⁶These authors contributed equally: Nada El Youssef, Aude Jegou, Christian Bénar and Fabrice Bartolomei. ✉email: fabrice.bartolomei@ap-hm.fr</footnote>
<table><caption>**Table 1.** Patients characteristics. HS: hippocampal sclerosis; DNET: dysembryoplastic neuroepithelial tumor; GN: ganglioglioma; PNH: periventricular nodular heterotopia; HH: hypothalamic hamartoma; PS: perinatal stroke; TBI: traumatic brain injury; FCD: focal cortical dysplasia; Bil: bilateral; Unk: unknown.</caption><thead><tr><th>Patient number</th><th>Age at time of SEEG (Y)</th><th>Gender</th><th>Epilepsy duration</th><th>Etiology</th><th>Epilepsy subtype</th><th>Side</th><th>Surgery</th><th>Engel</th><th>Number of analyzed seizure/LOC</th><th>Number of Seizure according to CSS (A/B/C)</th></tr></thead><tbody><tr><td>1</td><td>28</td><td>F</td><td>9</td><td>HS</td><td>MTLE</td><td>L</td><td>No</td><td>N/A</td><td>4/2</td><td>2/1/1</td></tr><tr><td>2</td><td>21</td><td>F</td><td>13</td><td>DNET</td><td>LTLE</td><td>R</td><td>Yes</td><td>4</td><td>3/0</td><td>3/0/0</td></tr><tr><td>3</td><td>36</td><td>F</td><td>21</td><td>HH</td><td>LTLE</td><td>Bil</td><td>No</td><td>N/A</td><td>1/0</td><td>1/0/0</td></tr><tr><td>4</td><td>30</td><td>F</td><td>9</td><td>PVHT</td><td>LTLE</td><td>R&gt;L</td><td>No</td><td>N/A</td><td>2/1</td><td>1/1/0</td></tr><tr><td>5</td><td>27</td><td>F</td><td>10</td><td>HS</td><td>MTLE</td><td>R</td><td>Yes</td><td>1a</td><td>2/0</td><td>2/0/0</td></tr><tr><td>6</td><td>26</td><td>F</td><td>3</td><td>GN</td><td>MTLE</td><td>R</td><td>Yes</td><td>1a</td><td>2/2</td><td>0/0/2</td></tr><tr><td>7</td><td>26</td><td>M</td><td>26</td><td>PS</td><td>LTLE</td><td>R&gt;L</td><td>No</td><td>N/A</td><td>2/2</td><td>0/1/1</td></tr><tr><td>8</td><td>44</td><td>M</td><td>4</td><td>TBI</td><td>LTLE</td><td>R</td><td>Yes</td><td>1a</td><td>3/3</td><td>0/3/0</td></tr><tr><td>9</td><td>45</td><td>F</td><td>18</td><td>HS</td><td>MTLE</td><td>L</td><td>Yes</td><td>1a</td><td>5/5</td><td>0/1/4</td></tr><tr><td>10</td><td>35</td><td>F</td><td>19</td><td>HS</td><td>MTLE</td><td>L</td><td>Yes</td><td>2a</td><td>1/1</td><td>0/0/1</td></tr><tr><td>11</td><td>33</td><td>M</td><td>22</td><td>TBI</td><td>LTLE</td><td>Bil</td><td>No</td><td>N/A</td><td>5/5</td><td>0/0/5</td></tr><tr><td>12</td><td>56</td><td>M</td><td>45</td><td>FCD</td><td>LTLE</td><td>R</td><td>Yes</td><td>1a</td><td>2/2</td><td>0/0/2</td></tr><tr><td>13</td><td>22</td><td>F</td><td>14</td><td>TBI</td><td>LTLE</td><td>R</td><td>Yes</td><td>3</td><td>1/1</td><td>0/0/1</td></tr><tr><td>14</td><td>44</td><td>M</td><td>15</td><td>UNK</td><td>LTLE</td><td>L</td><td>No</td><td>N/A</td><td>1/1</td><td>0/0/1</td></tr><tr><td>15</td><td>46</td><td>F</td><td>8</td><td>UNK</td><td>LTLE</td><td>L</td><td>Yes</td><td>3</td><td>1/1</td><td>0/0/1</td></tr><tr><td>16</td><td>29</td><td>F</td><td>10</td><td>UNK</td><td>LTLE</td><td>L</td><td>Yes</td><td>N/A</td><td>1/1</td><td>0/0/1</td></tr></tbody></table>

of cortical regions through inhibition of subcortical arousal systems affected by ictal discharge, resulting in neocortical inhibition and emergence of slow delta rhythms at the level of the frontoparietal neocortex.

In parallel, efforts have been made to quantify clinical and electrophysiological phenomena associated with altered consciousness in focal seizures. Three psychometric scales have been specifically proposed to measure ictal consciousness/awareness changes: the Ictal Consciousness Inventory (ICI)¹⁸, the Consciousness Seizure Scale (CSS)¹⁴ and the Responsiveness in Epilepsy Scale—versions I and II (RES-I and RES-II)¹⁹. The CSS is based on eight criteria measured by an epileptologist, with good inter and intra-observer reliability¹⁴.

Many neurophysiological methods have been used to quantify the complexity of signals across consciousness states, essentially in the domain of disorders of consciousness (DoC) such as coma or vegetative states. Entropy has been widely applied for quantification of brain activity, from the evaluation of consciousness alteration to assessment of brain networks interactions and the “ageing brain”²⁰. Permutation entropy (PE) as proposed by Bandt²¹, is a method to measure the complexity of dynamic systems. It thus allows assessment of dynamic changes in EEG recordings, with higher entropy being associated with the higher complexity of signals.

In this study based on intracerebral stereotactic-EEG recordings (SEEG), we computed delta entropy, (ΔE), corresponding to the decrease from baseline, and delta time metrics, derived from the permutation entropy, and correlated this measure to the awareness changes occurring during temporal lobe seizures. We tested whether altered awareness was correlated with the loss in signal complexity, a hypothesis previously formulated from connectivity data⁹. Finally, we estimated global as well as regional changes in PE to gain insight into the mechanisms associated with awareness impairment.

## Results

Sixteen patients with TLE were selected based on a clinical examination during the seizure that was informative enough to obtain a Consciousness Seizure Scale (CSS) score (Table 1). The CSS¹⁴ takes into account different components of consciousness leading to a composite score from 0 (no alteration of consciousness) to 9 (complete alteration of consciousness) and allows to distinguish different groups according to the extent of alteration of consciousness, with group A (score 0–1) corresponding to preserved consciousness, group B (score 2–5) to moderate alteration of consciousness and group C (score 6–9) to profound alteration of consciousness.

We excluded secondary generalized seizures and seizures that were not tested within the first 30 s after electrical onset. A total of thirty-six seizures from sixteen subjects were analyzed. According to the CSS scores, nine seizures were classified into group A, seven seizures into group B, and 20 seizures into group C, respectively.

**Comparison of permutation entropy between the three seizure groups.** Permutation entropy (PE) was computed on each channel, on a sliding window. The PE time courses of all seizures within each group were superimposed after subtraction of the mean entropy value of the baseline from the entropy value of the signal. Figure 1 illustrates changes in PE for 8 selected channels in one patient. “Delta entropy” (ΔE) is the difference between the entropy value at the baseline and the minimum entropy reached during the seizure.

As shown in Fig. 2A, three different patterns of ictal entropy dynamics were obtained for each group. A ΔE close to zero was seen in group A, while a high negative ΔE (corresponding to the lowest values of PE) was observed for Group C. Group B seizures disclosed intermediate ΔE values. There was a significant correlation
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Illustration of the Permutation Entropy estimation. (A) 3D Representation of the stereotactic-EEG electrodes in one patient. (B) Electrophysiological signal associated to 8 different selected contacts of the implanted channels (bipolar montage). (C) Permutation entropy signal of those channels showing the decrease of values after seizure onset. (D) Map Representation of the entropy signal with the different markers, as set by the scorer, corresponding to the beginning of the seizure (red), the maximal loss of consciousness (pink) and the end of the seizure (black). (E) Average of the entropy signal of all the channels during a seizure with the following markers: Beginning of the seizure (red), Loss of consciousness (pink), minimum of entropy value (cyan) (this marker is used to determine the delta entropy) and the end of seizure (black). The red and black markers are determined by the scorer based on the electrophysiological signal; the pink marker is determined by the scorer based on the CSS; the cyan marker is determined in an automated way. Arrow means "after applying".</figcaption></figure>

between the mean ΔE values and the CSS scores (Fig. 2B). The difference in the mean ΔE values was statistically significant (p < 0.001) between the groups A and C (Fig. 2C) and to a lesser extent between the groups C and B (p < 0.05). No correlation was found between the ΔE_time and the consciousness score.

**Thresholded ΔE values to discriminate between the groups A and C.** By using ROC method, we obtained two thresholds of ΔE that can distinguish the group, one at – 0.13 and the other at – 0.14 (Fig. 1S in Supplemental Material). Both these thresholds allowed perfect distinction between the group A and group C on our dataset (Fig. 2B). We thereby used the averaged threshold value of ΔE at – 0.135 for our analysis. Values in Group B were varying around the threshold (that was optimized for classification of C versus A) and showed larger distribution suggesting that this group is more heterogeneous.

**Entropy alterations in temporal and extra-temporal regions within groups A and C.** Figure 3 shows an illustrative example of a patient with two different seizures, one with no alteration of consciousness (group A) and another with severe alteration of consciousness (group C): Remarkably, these seizures disclose clearly different patterns of entropy changes, with more regions involved when alteration of consciousness is observed.

Since we observed a clear relationship between global ΔE and awareness level, we then investigated the variations of ΔE across different brain regions defined as explained above. We visualized the ΔE values in the patient's anatomy to assess which regions are more involved with regards to the alteration of awareness, by comparing Group A with Group C seizures. Figure 4A illustrates all electrodes of all patients of each respective group represented on a brain 3D mesh, with the mean ΔE values of each bipolar contact represented as spheres according to a color scale. Compared to the group A, group C seizures were characterized by a greater number of involved contacts in both hemispheres. The ΔE values for each anatomical region are indicated in Fig. 4B and C. Numerous regions of both hemispheres disclosed changes in delta entropy and thus showed high negative delta entropy values in Group C (left frontal, left insula, left lateral temporal, left mesiotemporal, left occipital, left occipito-temporal, right central, right insula, right lateral temporal, right mesiotemporal and right occipital), while few were involved in Group A.

Discussion

This is the first study to our knowledge that evaluates permutation entropy in relationship with alteration of awareness in focal seizures. We chose to study a relatively homogeneous group of patients with temporal lobe seizures. Our results reveal a positive correlation between the delta entropy (ΔE) and the consciousness score.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Visualization of the entropy changes by consciousness groups (Group A, B, C). (A) Mean entropy values averaged from all channels for all the seizures in each group. (B) Correlation between consciousness scores and mean delta entropy value with the threshold determined by ROC curve to differentiate Group C and Group A (threshold in orange at – 0.135) (C) Distribution of the delta entropy values according to each group. \* p<0.05 \*\*\* p<0.001.</figcaption></figure>

Furthermore, we were able to define a threshold of the entropy changes, measured with ΔE, allowing to discriminate between the seizures of group A (no alteration of awareness) and group C (severe alteration of awareness). Another important finding is that in patients with loss of awareness unlike in those with preserved awareness, the demonstrated changes in ΔE values showed more extended spatial distribution.

As far as various disorders of consciousness (DoC, e.g., syncope, coma, vegetative states, sleep disorders, encephalopathy and seizures) are concerned, it is the degree of impairment within the “consciousness system” that defines the extent of alteration of consciousness²³. Numerous studies have been conducted to identify the components of consciousness in patients with DoC and, most importantly, to identify physiological markers in the clinical setting and thereby to better define the patient's prognosis.

With respect to the structural imaging, findings demonstrated to be specific for patients with permanent DoC (Vegetative state, VS, and minimally conscious state, MCS) as compared to healthy controls, include diffuse abnormalities in the cortical, subcortical, and white matter volume and connectivity²⁴.

The use of entropy allows a successful evaluation of the "quantity of information" or the "complexity" of a system. Several theories predict that the complexity of information integration²⁵ or distributed processing²⁶ is elevated during conscious states. Permutation entropy assesses the regularity of the probabilistic distributions of temporal patterns in the signal²¹.

Permutation entropy has been studied by different teams in various DoC. Thul²⁷ found reduced permutation entropy on scalp EEG in patients with minimally conscious state and in the unresponsive wakefulness syndrome (UWS) compared to healthy controls. Furthermore, this reduction was most pronounced in the UWS group.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Stereotactic-EEG electrodes implantation of one subject having one seizure in Group A and one seizure in Group C. The red vertical line corresponds to seizure onset and the black vertical line, to seizure offset as determined visually by the scorer based on the SEEG signal.</figcaption></figure>

Another method developed by King<sup>28</sup> based on weighted symbolic mutual information (wSMI) allowed discrimination between VS and MCS. Sitt and colleagues<sup>29</sup> used different metrics obtained from scalp EEG to discriminate between different states of consciousness (low frequency power, EEG complexity and information exchange). They found that permutation entropy-based measures were particularly efficient in the theta frequency range to discriminate VS from the other groups. PE has been widely used as a tool to track the qualitative effect of anesthetics throughout the different stages (awake, to slightly and deeply sedated)<sup>30,31</sup> becoming a means of evaluating the depth of anesthesia monitoring<sup>32–34</sup>.

These studies confirmed that access to consciousness requires sufficient information content resulting in high signal complexity whereas altered states of consciousness are associated with low signal complexity.

Li<sup>31</sup> found that permutation entropy can be used to predict absence seizures using genetic absence epilepsy rats from Strasbourg (GAERS) model, and to demonstrate the presence of a pre-ictal state in absence seizures: a dynamic progressive decrease in PE was observed before seizure onset. Furthermore, Ferlazzo<sup>35</sup> showed higher permutation entropy in fronto-temporal areas and lower entropy values in posterior areas, during interictal and ictal periods, in patients with absence seizures in comparison to healthy controls.

The originality of our study is the use of permutation entropy in the context of alteration of consciousness occurring in patients with TLE, recorded with intracerebral electrodes (SEEG). In contrast with other disorders of consciousness, this state is transient and must be captured during the seizure periods. Taken as a whole, we observed a decrease of PE values at the seizure onset and throughout the seizure spread in the brain. This decrease was more pronounced in patients with an alteration of awareness. We estimated the decrease from baseline by a metric called “delta entropy” (ΔE) and we found a positive correlation between the ΔE value and the consciousness score: Group A seizures were characterized by a ΔE value close to zero and Group C seizures with a high negative ΔE. This means that during a seizure with no alteration of consciousness the variation in permutation entropy between the baseline and the minimal entropy value obtained is smaller in comparison to a seizure with severe alteration of consciousness. These results suggest that seizures with maximal alteration of consciousness also have maximal loss of signal complexity. We can thus observe a spectrum of decreased delta entropy as we go from group A to group C, which also applies to the extent of the regions involved as we move along these groups. Indeed, we noted an involvement of numerous brain regions in Group C seizures, while in Group A seizures, mostly the temporal lobe structures (left and right) were involved.

Thus, it appears that the loss of complexity observed in group A seizures is limited to temporal structures. The associative cortical regions keep their level of complexity, probably compatible with information exchange and conscious processing.

In group C seizures, the loss of signal complexity was more diffuse, extending bilaterally as well as to the associative cortices. In term of pathophysiology, these results suggest that the loss of signal complexity becomes incompatible with conscious processing. This is in agreement with different theories of consciousness based on information theory (IT)<sup>36</sup>, re-entrant processes (RP)<sup>37</sup> and global workspace theory (GW)<sup>13</sup>. Indeed, the relationship with the level of complexity alteration at the local level is in agreement with the IT or RP, while the extension to extra-temporal cortices is in agreement with an interruption of the information flow in the GW.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Visualization of the delta entropy values inside the brain. (A) Stereotactic-EEG electrodes of all subjects are projected into MNI template using GARDEL²² software allowing a 3D vision of the contacts. Then delta entropy value by contact for Group A and C are visualized, the size and the color being proportional to the values. (B) Visualization of delta entropy by anatomical regions for Group A and C. (C) Distribution of delta entropy values by anatomical regions with the threshold at – 0.135 (in orange).</figcaption></figure>

Noteworthy, no clear pattern was seen in group B seizures, where the delta entropy values were more heterogeneous. This can be explained by the fact that the extent of alteration of consciousness in this group is more diverse,
since the consciousness scores vary between 2 and 5 and thus these seizures can behave either like group A or group C seizures.

Furthermore, no positive correlation was found between the delta time ($\Delta_{\text{time}}$) and the consciousness score in any group. This could be compatible with the observation that the duration of a seizure does not necessarily determine the extent of alteration of consciousness: an isolated typical epileptic aura can last few minutes and usually does not impair the components of consciousness.

Finally, despite some limitations mentioned below, our results suggest that the delta entropy marker can be used as a potential predictor tool for assessing the extent of alteration of consciousness³⁸.

## Limitations

A limitation to our study is the small sample size: in fact, more data is needed to properly apply the threshold entropy marker and its usage as a predicting tool in terms of the extent of alteration of consciousness. Another limitation is the lack of sufficient electrode sampling of some regions: in the setting of temporal lobe epilepsy, considered as the main hypothesis during the pre-surgical evaluation, some cases did not have enough sampling of the extratemporal, particularly in the parietal regions. Additionally, not all patients had subcortical contacts in the thalamus. Thus, the extent of involvement of the thalamus and parietal regions remains not clearly defined by the current data.

Furthermore, the calculation of the threshold was determined on the same dataset used for the analysis; this method needs to be applied on another dataset to confirm the threshold value.

## Conclusion

Permutation entropy is a promising tool to discriminate between the different levels of consciousness alteration in seizures, and mainly between group A and group C by using a threshold value of delta entropy. Both the value of entropy changes and the spatial extent probably account for the consciousness alterations. This metric is relatively easy to use, especially when compared to connectivity methods (i.e. $h^2$ method¹⁴⁻¹⁶). It could then be implemented to estimate the therapeutic impact of certain neuromodulation techniques on improving consciousness during seizures, as observed with VNS³⁹⁻⁴² or DBS⁴³,⁴⁴. Future studies are needed to study permutation entropy and alteration of consciousness in the setting of extra-temporal epilepsies.

## Materials and methods

**Patients' selection.** Patients were retrospectively selected from our database of focal pharmaco-resistant epilepsy cases undergoing pre-surgical evaluation within the Epilepsy Unit at the Timone Hospital in Marseille using SEEG recording between October 2015 and December 2018. A systematic patient examination and interview has been conducted by nurses or physicians during seizures to assess the level of interaction as well as language impairment⁴⁵,⁴⁶. The study was approved by the institutional review board (Assistance Publique Hôpitaux de Marseille APHM 2021-179) and informed consent was obtained from each subject at the time of SEEG recording. All experiments were performed in accordance with relevant guidelines and regulations.

**SEEG recordings.** Recordings were performed using multi-contact depth electrodes (10-15 contacts, length: 2 mm, diameter: 0.8 mm, 1.5 mm apart) placed intracerebrally according to Talairach's stereotactic method⁴⁷. Signals were then recorded on a 128 or 256 channel Deltamed TM system and amplified at 256, 512 or 1024 Hz and recorded with no digital filter. Two analog filters were used at the time of recording: a high-pass filter (cut-off frequency equal to 0.16 Hz at -3 dB) and an anti-aliasing low-pass filter (cut-off frequency equal to 97 Hz at 256 Hz, or 170 Hz at 512 Hz or 340 Hz at 1024 Hz).

**Estimation of awareness impairment using the consciousness seizure scale (CSS).** We analyzed the alteration of awareness using the Consciousness Seizure Scale, CSS¹⁴. The scoring was done by two neurologists (NEY and FB) and the mean of the two scores was taken. The CSS relies on an eight-criterion scale taking into account the different components of consciousness: (1) unresponsiveness, (2) visual attention, (3) interaction, (4) consciousness of the seizure, (5) adapted behavior, (6) post-ictal amnesia and (7) interictal amnesia. Finally, the (8) "global appreciation criterion" based on the physician global assessment and scored between 0 and 2. The timing of the maximal alteration of consciousness (as determined by the scorer) was set as Tmax.

According to our previous studies¹⁴,¹⁵,⁴⁵, the seizures were classified into Group A (score 0 to 1), Group B (score 2 to 5) and Group C (score 6 to 9) with Group A corresponding to preserved awareness, Group B to intermediate alteration of awareness and Group C to profound alteration of awareness.

**Permutation entropy method.** In this study, we chose entropy to characterize the level of disorganization of the signals and to define the timing of maximal LOC as well as spatial extent of alteration of consciousness. The specific entropy method chosen was the "Permutation Entropy", PE²¹, which takes into account the complexity of the electrophysiological signals. PE is based on a measure of Shannon entropy of the relationships between neighbor values of the time series.

The method is rendered more robust than standard information-based method by calculating the entropy for n "embedding dimensions" in the range of 3-7, (i.e., different number of neighbors or samples). We used a fixed delay $\tau$ which indicates the distance between neighbors ($\tau=1$ for consecutive samples).

$$H(n) = - \sum p(\pi) \log p(\pi)$$
With $\pi$ describing the dynamics between the neighbor values (e.g., $x_t < x_{t+1}$), $p(\pi)$ probability of this order and $n$ embedding dimensions. H (n) gives a value between 0 and 1, n meaning that the signal is perfectly organized (for example in the case of a pure sine wave) and 1 completely disorganized (completely random noise). The algorithm used to apply permutation entropy on our data is the one developed by Unakafova$^{48}$ due to its algorithmic performance. The efficiency of this method stems from the use of pre-filled tables to describe the order relation between neighbor values.

To normalize and compare the entropy time course across patients and seizures, we selected seizures with a minimum baseline of 30 s before seizure onset. Before applying permutation entropy on a given seizure signal, we performed careful selection of channels by keeping those in the grey matter and removing the white matter channels. A bipolar montage was used (i.e., difference of consecutive channels). Channels with artefact or high level of noise were removed. Permutation entropy was computed on each channel, on a sliding window of 10 s with an overlap of 5 s, an order of 3 and a delay of 1 sample. We studied entropy at two different scales: on a global level by averaging all SEEG contacts, and on a regional level by using anatomical regions (see Supplementary material for more details). For these two different scales, each entropy signal was normalized by subtracting the mean of the baseline channel by channel. To measure changes in entropy, we used two different variables named “delta entropy” (ΔE), and “delta time” (Δtime) (Fig. 1). “Delta entropy” is the difference between the entropy value at the baseline and the minimum entropy reached during the seizure, “Delta time” is the time between the end of baseline (corresponding to the beginning of the seizure) and the time at which the minimum entropy occurs. Figure 1 shows these two variables, the signal markers that were set by the scorer corresponding to the beginning and end of the seizure (T0 and Tend) and the time of maximal alteration of consciousness (for groups B and C) as defined by the CSS.

For the global level analysis, we calculated the mean entropy by using all channels of the montage, allowing to test for a global difference across the groups and to determine a threshold separating the seizures of Group A and Group C. Using the classification, A versus C, we computed the false positive rate, the true positive rate and the F1 score of the different threshold values (-0.3: -0.01:0) of delta entropy. We thus obtained a ROC curve leading to an optimal threshold (supplemental fig. S1). All the seizures with a delta entropy value inferior to the threshold were labelled as Group C (Fig. 3).

CT-scan with electrodes/preoperative MRI data co-registration was performed using the open-source software GARDEL$^{22}$ and projected into the MNI template for 3D representation shown in Fig. 4. For the regional level analysis, each bipolar SEEG contact was assigned to a brain region by using an automated parcellation according to the VEP atlas$^{49,50}$ with accuracy visually checked by a neurologist. For the purpose of this study, the following anatomical regions (left and right hemispheric) were defined by regrouping several VEP atlas regions: mesio-temporal, lateral temporal, central, frontal, parietal, occipital, insular, occipito-temporal, parieto-occipital, deep gray matter, thalamus (see Supplementary material for details). Only the regions sampled in at least two seizures were kept for the analysis. For each anatomical region and each subject, the delta entropy values (ΔE) from all contributing channels were averaged, resulting in one value per region. By assessing the ΔE per anatomical region for all the seizures of each respective group (A, B and C), we estimated which brain regions are more involved during the change in awareness. However, due to the limitations inherent to the SEEG method, the number of contacts per region differed between the patients and the number of contributing subjects differed between the regions, resulting in a lack of statistical power of such analysis (Fig. 1, Supplemental material).

**Statistical analysis.** We computed the linear correlation between the consciousness score and delta entropy using Pearson method. Then a linear mixed-effect model (LMM) using the “lme4” package (version 1.1-26$^{51}$, of the statistical software R version 4.0.3 (R Core Team, 2020) was used to study the effect of the group (A, b or C) on the ΔE value, with subject as random effect.

Restricted Maximum Likelihood (REML) was used to fit the LMM and Satterthwaite's method was used to estimate the degrees of freedom. Post-hoc pairwise comparisons were performed with the Least Square Means (LSM) method. Importantly, the use of a mixed model allows to take into account variability across patients (Table 2, Supplemental material).

Concerning the different anatomical regions, due to the above-mentioned limitations, it was impossible to perform a statistically powerful analysis of the ΔE per region to see the effect of the group, neither a statistical comparison between the regions and the seizure groups.

**Data availability**

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Received: 28 August 2022; Accepted: 6 December 2022

Published online: 24 December 2022

**References**

1. Maillard, L. et al. Semiologic and electrophysiologic correlations in temporal lobe seizure subtypes. *Epilepsia* **45**(12), 11 (2004).

2. Lambert, I. & Bartolomei, F. Why do seizures impair consciousness and how can we reverse this?. *Curr. Opin. Neurol.* **33**(2), 173-178. https://doi.org/10.1097/WCO.0000000000000794 (2020).

3. Yu, L. & Blumenfeld, H. Theories of impaired consciousness in epilepsy. *Ann. N. Y. Acad. Sci.* **1157**, 48-60. https://doi.org/10.1111/j.1749-6632.2009.04472 (2009).

4. Blumenfeld, H. Impaired consciousness in epilepsy. *Lancet Neurol.* **11**(9), 814-826. https://doi.org/10.1016/S1474-4422(12)70188-6 (2012).
5. Bensaid, S., Modolo, J., Merlet, I., Wendling, F. & Benquet, P. COALIA: A computational model of human EEG for consciousness research. *bioRxiv* **2019**, 575043 (2019).

6. Parvizi, J. & Damasio, A. Consciousness and the brainstem. *Cognition* **79**(1), 135–160. https://doi.org/10.1016/S0010-0277(00)00127-X (2001).

7. Fisher, R. S. et al. Operational classification of seizure types by the International League Against Epilepsy: Position Paper of the ILAE Commission for Classification and Terminology. *Epilepsia* **58**(4), 522–530. https://doi.org/10.1111/epi.13670 (2017).

8. Zhao, C. W. et al. Reliability of patient self-report of cognition, awareness, and consciousness during seizures. *Ann. Clin. Transl. Neurol.* **9**(1), 16–29. https://doi.org/10.1002/acn3.51485 (2022).

9. Bartolomei, F. & Naccache, L. The global workspace (GW) theory of consciousness and epilepsy. *Behav. Neurol.* **24**(1), 67–74. https://doi.org/10.1155/2011/127864 (2011).

10. Bartolomei, F. Coherent neural activity and brain synchronization during seizure-induced loss of consciousness. *Arch. Ital. Biol.* **150**(2–3), 164–171. https://doi.org/10.4449/aib.v150i2.1252 (2012).

11. Danielson, N. B., Guo, J. N. & Blumenfeld, H. The default mode network and altered consciousness in epilepsy. *Behav. Neurol.* **24**(1), 55–65. https://doi.org/10.3233/BEN-2011-0310 (2011).

12. Blumenfeld, H. Arousal and consciousness in focal seizures. *Epilepsy Curr.* **21**, 353–359 (2021).

13. Dehaene, S. & Naccache, L. Towards a cognitive neuroscience of consciousness: Basic evidence and a workspace framework. *Cognition* **79**(1–2), 1–37. https://doi.org/10.1016/s0010-0277(00)00123-2 (2001).

14. Arthuis, M. et al. Impaired consciousness during temporal lobe seizures is related to increased long-distance cortical-subcortical synchronization. *Brain J. Neurol.* **132**, 2091–2101 (2009).

15. Lambert, I., Arthuis, M., McGonigal, A., Wendling, F. & Bartolomei, F. Alteration of global workspace during loss of consciousness: A study of parietal seizures: Alteration of global workspace during loss of consciousness. *Epilepsia* **53**(12), 2104–2110. https://doi.org/10.1111/j.1528-1167.2012.03690.x (2012).

16. Bonini, F., Lambert, I., Wendling, F. & McGonigal, A. Altered synchrony and loss of consciousness during frontal lobe seizures. *Clin. Neurophysiol.* **127**(2), 1170–1175. https://doi.org/10.1016/j.clinph.2015.02.050 (2016).

17. Dheer, P., Pati, S., Chowdhury, K. K. & Majumdar, K. K. Enhanced gamma band mutual information is associated with impaired consciousness during temporal lobe seizures. *Heliyon* **6**(12), e05769. https://doi.org/10.1016/j.heliyon.2020.e05769 (2020).

18. Cavanna, A. E. et al. Measuring the level and content of consciousness during epileptic seizures: The Ictal Consciousness Inventory. *Epilepsy Behav. EB* **13**(1), 184–188. https://doi.org/10.1016/j.yebeh.2008.01.009 (2008).

19. Nani, A. & Cavanna, A. E. The quantitative measurement of consciousness during epileptic seizures. *Epilepsy Behav EB.* **30**, 2–5. https://doi.org/10.1016/j.yebeh.2013.09.007 (2014).

20. Keshmiri, S. Entropy and the brain: An overview. *Entropy (Basel)* **22**(9), 917. https://doi.org/10.3390/e22090917 (2020).

21. Bandt, C. & Pompe, B. Permutation entropy: A natural complexity measure for time series. *Phys. Rev. Lett.* **88**(17), 174102. https://doi.org/10.1103/PhysRevLett.88.174102 (2002).

22. Medina-Villalon, S. et al. EpiTools, A software suite for presurgical brain mapping in epilepsy: Intracerebral EEG. *J. Neurosci.* Methods **303**, 7–15. https://doi.org/10.1016/j.jneumeth.2018.03.018 (2018).

23. Boly, M., Seth, A. K., Wilke, M. & Ingmundson, P. Consciousness in humans and non-human animals: Recent advances and future directions. *Front. Psychol.* **4**, 55. https://doi.org/10.3389/fpsyg.2013.00625 (2013).

24. Snider, S. B. & Edlow, B. L. MRI in disorders of consciousness. *Curr. Opin. Neurol.* **33**(6), 676–683. https://doi.org/10.1097/WCO.0000000000000873 (2020).

25. Tononi, G. Consciousness as integrated information: A provisional manifesto. *Biol. Bull.* **215**(3), 216–242. https://doi.org/10.2307/25470707 (2008).

26. Dehaene, S. & Changeux, J.-P. Experimental and theoretical approaches to conscious processing. *Neuron* **70**(2), 200–227. https://doi.org/10.1016/j.neuron.2011.03.018 (2011).

27. Thul, A. et al. EEG entropy measures indicate decrease of cortical information processing in disorders of consciousness. *Clin. Neurophysiol. Off. J. Int. Fed Clin. Neurophysiol.* **127**(2), 1419–1427. https://doi.org/10.1016/j.clinph.2015.07.039 (2016).

28. King, J. R., Sitt, J. D., Faugeras, F. & Rohaut, B. Information sharing in the brain indexes consciousness in noncommunicative patients. *Curr. Biol.* **23**(19), 1914–1919. https://doi.org/10.1016/j.cub.2013.07.075 (2013).

29. Sitt, J. D., King, J. R., El-Karoui, I. & Rohaut, B. Large scale screening of neural signatures of consciousness in patients in a vegetative or minimally conscious state. *Brain* **137**, 2258–2270. https://doi.org/10.1093/brain/awu141 (2014).

30. Olofsen, E., Sleigh, J. W. & Dahan, A. Permutation entropy of the electroencephalogram: A measure of anaesthetic drug effect. *Br. J. Anaesth.* **101**(6), 810–821. https://doi.org/10.1093/bja/aen290 (2008).

31. Li, X., Cui, S. & Voss, L. J. Using permutation entropy to measure the electroencephalographic effects of sevoflurane. *Anesthesiology* **109**(3), 448–456. https://doi.org/10.1097/ALN.0b013e318182a91b (2008).

32. Kim, P.-J., Kim, H.-G., Noh, G.-J., Koo, Y.-S. & Shin, T.-J. Usefulness of permutation entropy as an anesthetic depth indicator in children. *J. Pharmacokinet. Pharmacodyn.* **42**(2), 123–134. https://doi.org/10.1007/s10928-015-9405-5 (2015).

33. Kreuzer, M., Kochs, E. F., Schneider, G. & Jordan, D. Non-stationarity of EEG during wakefulness and anaesthesia: Advantages of EEG permutation entropy monitoring. *J. Clin. Monit. Comput.* **28**(6), 573–580. https://doi.org/10.1007/s10877-014-9553-y (2014).

34. Lange, N. et al. Permutation entropy in intraoperative aCoG of brain tumour patients in awake tumour surgery—a robust parameter to separate consciousness from unconsciousness. *Sci. Rep.* **9**(1), 16482. https://doi.org/10.1038/s41598-019-52949-1 (2019).

35. Ferlazzo, E. et al. Permutation entropy of scalp EEG: A tool to investigate epilepsies: Suggestions from absence epilepsies. *Clin. Neurophysiol. Off. J. Int. Fed Clin. Neurophysiol.* **125**(1), 13–20. https://doi.org/10.1016/j.clinph.2013.06.023 (2014).

36. Tononi, G. Consciousness as integrated information: A provisional manifesto. *Biol. Bull.* **215**(3), 216–242. https://doi.org/10.2307/25470707 (2020).

37. Lamme, V. A. & Roelfsema, P. R. The distinct modes of vision offered by feedforward and recurrent processing. *Trends Neurosci.* **23**(11), 571–579. https://doi.org/10.1016/s0166-2236(00)01657-x (2000).

38. Lüders, H. et al. Proposal: Different types of alteration and loss of consciousness in epilepsy. *Epilepsia* **55**(8), 1140–1144. https://doi.org/10.1111/epi.12595 (2014).

39. Bodin, C. et al. Responders to vagus nerve stimulation (VNS) in refractory epilepsy have reduced interictal cortical synchronicity on scalp EEG. *Epilepsy Res.* **113**, 98–103. https://doi.org/10.1016/j.elepsresys.2015.03.018 (2015).

40. Bartolomei, F. et al. How does vagal nerve stimulation (VNS) change EEG brain functional connectivity?. *Epilepsy Res.* **126**, 141–146. https://doi.org/10.1016/j.elepsresys.2016.06.008 (2016).

41. Ravan, M., Sabesan, S. & D'Cruz, O. On quantitative biomarkers of VNS therapy using EEG and ECG signals. *IEEE Trans. Biomed. Eng.* **64**(2), 419–428. https://doi.org/10.1109/TBME.2016.2554559 (2017).

42. Sangare, A. et al. The effectiveness of vagus nerve stimulation in drug-resistant epilepsy correlates with vagus nerve stimulation-induced electroencephalography desynchronization. *Brain Connect.* **10**(10), 566–577. https://doi.org/10.1089/brain.2020.0798 (2020).

43. Yu, T. et al. High-frequency stimulation of anterior nucleus of thalamus desynchronizes epileptic network in humans. *Brain* **141**(9), 2631–2643 (2018).

44. Deutschová, B. et al. Thalamic oscillatory activity may predict response to deep brain stimulation of the anterior nuclei of the thalamus. *Epilepsia* **62**(5), e70–e75. https://doi.org/10.1111/epi.16883 (2021).
45. Bonini, F. et al. Altered synchrony and loss of consciousness during frontal lobe seizures. *Clin. Neurophysiol.* **127**(2), 1170–1175.
https://doi.org/10.1016/j.clinph.2015.02.050 (2016).

46. Trebuchon, A. et al. The different patterns of seizure-induced aphasia in temporal lobe epilepsies. *Epilepsy Behav.* **EB 78**, 256–264.
https://doi.org/10.1016/j.yebeh.2017.08.022 (2018).

47. Lagarde, S. et al. Interictal stereotactic-EEG functional connectivity in refractory focal epilepsies. *Brain* **141**(10), 2966–2980.
https://doi.org/10.1093/brain/awy214 (2018).

48. Unakafova, V. A. & Keller, K. Efficiently measuring complexity on the basis of real-world data. *Entropy* **15**(10), 4392–4415. https://doi.org/10.3390/e15104392 (2013).

49. Wang, H. E. et al. VEP atlas: An anatomic and functional human brain atlas dedicated to epilepsy patients. *J. Neurosci. Methods* **348**, 108983. https://doi.org/10.1016/j.jneumeth.2020.108983 (2021).

50. Makhalova, J. et al. Virtual epileptic patient brain modeling: Relationships with seizure onset and surgical outcome. *Epilepsia* **63**, 1942–1955. https://doi.org/10.1111/epi.17310 (2022).

51. Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. *J. Stat. Softw.* **67**, 1–48. https://doi.org/10.18637/jss.v067.i01 (2015).

## Acknowledgements

We thank Dr Romain Carron and Pr Didier Scavarda for stereotactic and surgical procedures.

## Author contributions

(1) Conception and design of the study were done by F.B. and L.N. (2) Acquisition and analysis of data; N.E., J.S., A.J., C.B. (3) Writing the manuscript and elaboration of figures: N.E., F.B., L.N., J.S., C.B., A.J.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-25861-4.

**Correspondence** and requests for materials should be addressed to F.B.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022