# Respiratory dysfunction in two rodent models of chronic epilepsy and acute seizures and its link with the brainstem serotonin system

Hayet Kouchi, Michaël Ogier, Gabriel Dieuset, Anne Morales, Béatrice Georges, Jean-Louis Rouanet, Benoît Martin, Philippe Ryvlin, Sylvain Rheims, Laurent Bezin

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

## ► To cite this version:

Hayet Kouchi, Michaël Ogier, Gabriel Dieuset, Anne Morales, Béatrice Georges, et al.. Respiratory dysfunction in two rodent models of chronic epilepsy and acute seizures and its link with the brainstem serotonin system. Scientific Reports, 2022, 12 (1), pp.10248. 10.1038/s41598-022-14153-6. hal-03716372

HAL Id: hal-03716372

https://hal.science/hal-03716372v1

Submitted on 7 Jul 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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 CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

OPEN

# Respiratory dysfunction in two rodent models of chronic epilepsy and acute seizures and its link with the brainstem serotonin system

Hayet Kouchi<sup>1,2,✉</sup>, Michaël Ogier<sup>3</sup>, Gabriel Dieuset<sup>1,4</sup>, Anne Morales<sup>1,2</sup>, Béatrice Georges<sup>1,2</sup>, Jean-Louis Rouanet<sup>5</sup>, Benoît Martin<sup>1,4</sup>, Philippe Ryvlin<sup>6</sup>, Sylvain Rheims<sup>1,2,7,8</sup> & Laurent Bezin<sup>1,2,8,✉</sup>

Patients with drug-resistant epilepsy can experience respiratory alterations, notably during seizures. The mechanisms underlying long-term alterations in respiratory function remain unclear. As the brainstem 5-HT system is a prominent modulator of respiratory function, this study aimed at determining whether epilepsy is associated with alterations in both the respiratory function and brainstem serotonin (5-HT) system in rats. Epilepsy was triggered by pilocarpine-induced status epilepticus in rats. Our results showed that 30–50% of epileptic (EPI) rats exhibited a sharp decrease in oxygen consumption (SDOC), low metabolic rate of oxygen, and slow regular ventilation (EPI/SDOC + rats). These alterations were detected only in rats with chronic epilepsy, independent of behavioral seizures, were persistent over time, and not associated with death. In these rats, 5-HT fiber density in the nucleus tractus solitarius was lower than that in the control and EPI/SDOC- rats. Both EPI/SDOC + rats and DBA/2 mice that present with audiogenic-induced seizure followed by fatal respiratory arrest—a model of sudden and expected death in epilepsy—had increased transcript levels of tryptophan hydroxylase 2 and 5-HT presynaptic transporter. Thus, our data support that 5-HT alterations are associated with chronic and acute epilepsy-related respiratory dysfunction.

## Abbreviations

<table><tr><td>5-HT</td><td>Serotonin</td></tr><tr><td>AGS</td><td>Audiogenic seizures</td></tr><tr><td>EPI</td><td>Epileptic</td></tr><tr><td>CIH</td><td>Chronic intermittent hypoxia</td></tr><tr><td>EEG</td><td>Electro-encephalography</td></tr><tr><td>EPI/SDOC-</td><td>Epileptic rat without sharp decrease in oxygen consumption</td></tr><tr><td>EPI/SDOC+</td><td>Epileptic rat with sharp decrease in oxygen consumption</td></tr><tr><td>f<sub>R</sub></td><td>Respiratory frequency</td></tr><tr><td>FRA</td><td>Fatal respiratory arrest</td></tr><tr><td>HCVR</td><td>Hypercapnia ventilatory response</td></tr><tr><td>HIS</td><td>Handling-induced seizures</td></tr><tr><td>HVR</td><td>Hypoxia ventilatory response</td></tr></table>

<footnote n="1">¹CRNL, Lyon Neuroscience Research Center, Inserm U1028, CNRS UMR5292, Lyon 1 University, 95 bd Pinel, 69500 Bron, France. ²Epilepsy Institute IDEE, Bron, France. ³IRBA, Institut de Recherche Biomédicale des Armées, Brétigny-sur-Orge, France. ⁴Univ Rennes, CHU Rennes, Inserm, LTSI (Laboratoire de Traitement du Signal et de L'Image), UMR-1099, Rennes, France. ⁵Laboratoire d'Ecologie des Hydrosystèmes Naturels et Anthropisés, U5023, Lyon 1 University, Lyon, France. ⁶Department of Clinical Neurosciences, Centre Hospitalier Universitaire Vaudois, Université de Lausanne, Lausanne, Switzerland. ⁷Department of Functional Neurology and Epileptology, Hospices Civils de Lyon, Lyon, France. ⁸These authors jointly supervised this work: Sylvain Rheims and Laurent Bezin. ✉email: hayet.kouchi@univ-lyon1.fr; laurent.bezin@univ-lyon1.fr</footnote>
<table><tr><td>NFRA</td><td>Non-fatal respiratory arrest</td></tr><tr><td>GTCS</td><td>Generalized tonic-clonic seizures</td></tr><tr><td>LTSI</td><td>Laboratoire de Traitement du Signal et de l'Image</td></tr><tr><td>NTS</td><td>Nucleus tractus solitarius</td></tr><tr><td>OSA</td><td>Obstructive sleep apnea syndrome</td></tr><tr><td>Pilo-SE</td><td>Pilocarpine-induced status epilepticus</td></tr><tr><td>SE</td><td>status epilepticus</td></tr><tr><td>SERT</td><td>Serotonin transporter</td></tr><tr><td>SSRI</td><td>Selective serotonin reuptake inhibitor</td></tr><tr><td>SUDEP</td><td>Sudden unexpected death in epilepsy</td></tr><tr><td>T<sub>E</sub></td><td>Expiration time</td></tr><tr><td>T<sub>I</sub></td><td>Inspiration time</td></tr><tr><td>TLE</td><td>Temporal lobe epilepsy</td></tr><tr><td>Tph2</td><td>Tryptophan hydroxylase 2</td></tr><tr><td>T<sub>tot</sub></td><td>Total respiratory cycle duration</td></tr><tr><td>V̇<sub>E</sub></td><td>Minute volume</td></tr><tr><td>VLM</td><td>Ventrolateral medulla</td></tr><tr><td>V̇O<sub>2</sub></td><td>Oxygen consumption</td></tr><tr><td>V<sub>T</sub></td><td>Tidal volume</td></tr></table>

The alteration of respiratory function is a chronic condition in epileptic patients that may represent a significant health concern, particularly for those with drug-resistant epilepsy¹,². For example, 33–86% of seizures expressed by patients with drug-resistant epilepsy are subsequently associated with an oxygen desaturation of < 90%³,⁴. In addition, the prevalence of obstructive sleep apnea (OSA) syndrome is twofold higher in epileptic patients compared to the general population⁵,⁶. Finally, it may be that preexisting respiratory dysfunction, particularly during the interictal period, represents a risk factor for sudden unexpected death in epilepsy (SUDEP), which results primarily from postictal central apnea³. The first objective of our study was to characterize in a preclinical model of drug-resistant temporal lobe epilepsy (TLE) respiratory disorders that might exist during the interictal period. Epilepsy following pilocarpine-induced status epilepticus (Pilo-SE) in rats is one of the most studied models of TLE. Like patients with TLE, rats that became epileptic after Pilo-SE are notably characterized by chronic, spontaneous, and drug-resistant seizures⁷⁻⁹. Respiratory function of such rats has been studied previously after the onset of epilepsy and independent of seizure activity. However, this required placing the animals under anesthesia and mechanical ventilation, a condition that can affect respiratory outcomes and limit interpretations¹⁰. We therefore first explored in unanesthetized and unrestrained rats subjected to Pilo-SE whether respiratory alterations occurred during epileptogenesis or after the onset of epilepsy

It is also crucial to fill in the knowledge gap regarding epilepsy-related respiratory alterations and the underlying pathological mechanisms. In this regard, the serotonin (5-HT) system may play a key role, since it has been involved in both physiological and pathophysiological mechanisms of breathing and epilepsy¹¹,¹². In addition, post-mortem examination of the brainstem of patients who died from SUDEP showed changes in the ventrolateral medulla (VLM)¹³, a region containing different populations of respiratory neurons, including those of the pre-Bötzinger complex, a kernel generator of the inspiratory rhythm¹⁴. Specifically, SUDEP patients presented with an alteration of the brainstem 5-HT system with a notable decrease in both tryptophan hydroxylase 2 (TPH2; rate-limiting enzyme in 5-HT synthesis) and the 5-HT presynaptic transporter (SERT; a protein that uptakes 5-HT from the synaptic cleft back to the presynaptic neurons) within the medullary raphe nuclei and/or the VLM¹³. In preclinical studies, DBA/1 and DBA/2 mouse strains are frequently used as SUDEP models, because they exhibit fatal respiratory arrest (FRA) after induction of an audiogenic seizure (AGS). In these mice, brainstem 5-HT abnormalities have been reported¹⁵⁻¹⁷, and FRAs were prevented when mice were pre-treated with fluoxetine before AGS, a selective 5-HT reuptake inhibitor (SSRI)¹⁸,¹⁹. These mice prone to AGS-induced FRA present however some limitations to perform a complete characterization of respiratory dysfunction in chronic epilepsy and its potential associated 5-HT abnormalities, because they are models of induced seizures and do not develop spontaneous recurrent seizures. Thus, the second objective of our study was to investigate brainstem 5-HT system alterations related to epilepsy in Pilo-SE rats by measuring the number of 5-HT-positive cell bodies, the density of 5-HT-positive terminals, and TPH2 and SERT mRNA levels. Given the preventive effect of SSRIs on FRA in DBA/2 mice, we also investigated whether TPH2 and SERT mRNA levels could distinguish DBA/2 mice prone to FRA from those that are not.

## Results

**Development of behavioral spontaneous recurrent seizures (SRSs) in rats after Pilo-SE.** Rats that did not develop Pilo-SE or died during SE were excluded from experiments and their numbers are reported in Supplementary Information S1. SRSs were defined as detailed in the Methods section. Briefly, out of the 82 rats subjected to Pilo-SE in experiments 1–3, 59 developed at least two consecutive SRSs by the end of the 2nd week post-SE, and these were then defined as EPI rats.

According to the evaluation of two experimenters, the number of recorded seizures (stages 3–5 of severity according to Racine's scale) per 12 h of daylight period was 1.1 ± 0.1 in EPI rats at 3–8 weeks after SE. In experiment 1, none of the Pilo-SE rats developed two consecutive handling-induced seizures (HIS) at 1 week post-SE, in agreement with previous studies.
<figure><img src="image_4.png" /><figcaption>B. EPI rat at one week post-SE</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>C. EPI rat with SDOC at three weeks post-SE</figcaption></figure>

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

**Figure 1.** Detection of SDOC events during a 24-h monitoring of oxygen consumption in a subgroup of epileptic rats. Oxygen consumption (VO₂) (light period from 06:00 am to 6:00 pm) was measured using respiratory thermochemistry. Grey bar above each x-axis indicates the dark period. A single profile was observed in control rats (A), while two distinct profiles were distinguished in epileptic rats (B,C). The first profile was similar to that observed in control rats (B), the second one was characterized by several episodes of sharp decrease in oxygen consumption (SDOC) (C), occasionally reaching zero. These SDOC events were only observed in epileptic rats, thereby dichotomizing two groups of epileptic rats (EPI/SDOC+ and EPI/SDOC-).

**A subset of EPI rats presented with abnormal pattern of oxygen consumption.** In experiment 1, the metabolic rate of oxygen (VO₂) in 20 EPI rats and 8 controls was recorded for 24 h each week, 1–8 weeks after SE. A thorough screening of VO₂ was then performed for each rat.

*Presence of sharp decreases in oxygen consumption (SDOCs).* Control rats (n=8) presented with slight up and down fluctuations of VO₂ during both light and dark periods, as illustrated in one representative rat (Fig. 1A). For better characterization of daily VO₂ fluctuations, the average VO₂ for each rat during the light and dark periods was calculated. Subsequently, the lowest value below the average VO₂ was determined for each rat; this value ranged between 38 and 70% of their average VO₂, independently of the circadian rhythm. This pattern remained unchanged during recordings over the next 7 weeks. In rats subjected to Pilo-SE, the pattern of daily VO₂ was similar to that of controls during the first 2 weeks post-SE, as illustrated in one representative Pilo-SE
rat (Fig. 1B). However, by the end of the 2nd week post-SE, two groups of EPI rats were identified according to $\dot{V}O_2$ fluctuations. A subset of EPI rats presented with a large amplitude of $\dot{V}O_2$ decrease ("Sharp Decrease(s) in Oxygen Consumption"; SDOCs), defined as a decrease greater than -62%, a percentage that corresponds to the maximal decrease of $\dot{V}O_2$ in control rats. Such EPI rats were thus defined as EPI/SDOC+ rats and their lowest $\dot{V}O_2$ values ranged between 0 and 32% of the average $\dot{V}O_2$, as illustrated at 3 weeks post-SE in one representative rat (Fig. 1C). In contrast, the remaining EPI rats (EPI/SDOC- rats) were similar to controls with regard to the lowest $\dot{V}O_2$ values, which ranged between 41 and 67% of the average $\dot{V}O_2$. We cannot exclude the fact that some SDOCs could be associated with the presence of focal seizures (with subclinical manifestation). However, video monitoring confirmed that SDOCs were not associated with behavioral seizures, as illustrated by an SDOC event that occurred while the rat was in a sleep-like position (Fig. 2A). In contrast, behavioral seizures were associated with a long-lasting increase in $\dot{V}O_2$, as illustrated during a stage-5 seizure (Fig. 2B).

The highest $\dot{V}O_2$ in control rats ranged from 164 to 229% of the average $\dot{V}O_2$ (median = 183%). The highest fluctuations within the daily $\dot{V}O_2$ pattern in EPI/SDOC+ and EPI/SDOC- rats were not statistically different from controls, and ranged from 144 to 227% of the average $\dot{V}O_2$ (median = 180%). However, greater increases in $\dot{V}O_2$ were observed in EPI/SDOC+ rats than in EPI/SDOC- rats (+93 ± 5% vs +75 ± 4%; P = 0.01).

As a good indicator of the stability of the $\dot{V}O_2$ pattern, the intra-individual variation (coefficient of variation, CV) of the daily $\dot{V}O_2$ pattern (CV-$\dot{V}O_2$) was calculated independently of light and dark periods. The mean CV-$\dot{V}O_2$ of the control group was estimated to be 20.8 ± 0.5% and 22 ± 0.7% for the light and the dark periods, respectively. The mean CV-$\dot{V}O_2$ of EPI/SDOC- rats was almost similar to that of controls with 19.1 ± 0.8% and 22.8 ± 0.7% for the light and dark periods, respectively. However, the mean CV-$\dot{V}O_2$ was significantly higher in EPI/SDOC+ rats, than in both the control (P = 0.002 and P = 0.004 for light and dark periods, respectively) and EPI/SDOC- groups (P < 0.0001 and P = 0.017 for light and dark periods, respectively), reaching 27.4 ± 2.9% and 26.5 ± 1.3% for light and dark periods, respectively.

*Low metabolic rate of oxygen in EPI/SDOC+ rats.* Since $\dot{V}O_2$ was measured weekly for all rats during 8 weeks, it was possible to retrospectively assign each EPI rat into either the EPI/SDOC+ or the EPI/SDOC- group from the first week after SE.

We investigated the progression of daily $\dot{V}O_2$ in each group. We first verified that SDOC events in EPI/SDOC+ rats did not affect the average value of daily $\dot{V}O_2$, both during the light (P = 0.86) and the dark (P = 0.93) periods. In order to compare the daily $\dot{V}O_2$ during the 8 weeks of measurement between the three groups of rats, we used the data without SDOCs in EPI/SDOC+ rats. We showed, in control and EPI/SDOC- rat groups, that daily $\dot{V}O_2$ decreased progressively by 21.3 ± 2.4% (P = 0.05) and 22.0 ± 2.2% (P = 0.02), respectively, between weeks 1–2 and 8 post-SE. In contrast, the decrease observed in EPI/SDOC+ rats was much larger during the same period (-69.5 ± 1.1%; P = 0.002), and reached -53 ± 6.0% (P = 0.02) by the fifth week post-SE (Fig. 3A). We also found that the decrease of daily $\dot{V}O_2$ coincided with the onset of SDOC events (Fig. 3B).

In each group, we verified that daily $\dot{V}O_2$ did not change significantly between 5 and 8 weeks post-SE; therefore, the average daily $\dot{V}O_2$ has been calculated during that period for each group. We found that daily $\dot{V}O_2$ was greatly reduced in EPI/SDOC+ rats at 5–8 weeks post-SE compared to both EPI/SDOC- rats and controls (-39% ± 7% and -34% ± 7%, respectively). This decrease was observed during both the light period (EPI/SDOC+ vs. EPI/SDOC- rats: P < 0.001; EPI/SDOC+ vs. control rats: P = 0.002) and the dark period (EPI/SDOC+ vs. EPI/SDOC- rats: P < 0.001; EPI/SDOC+ vs. control rats: P = 0.002) (See supplementary Fig. S1A online). At weeks 1–2 post-SE, no difference was observed between groups during both the light and dark periods (See supplementary Fig. S1A online).

In the three groups of rats (controls, EPI/SDOC+ and EPI/SDOC-), we evaluated the effect of the following three factors known to influence $\dot{V}O_2^{20}$: time ($t_0$ corresponds to the time of induction of SE in EPI rats), body weight, and circadian rhythm. A full description of these effects is provided in Supplementary Information S1.

*Characterization of SDOC.* The number of SDOC events during the 24-h measurement ranged between 12 and 414 events (median = 44), with no difference between the dark and light periods. The median duration of SDOCs was 5 min (range: 20 s to 13 min). The longest SDOC event is illustrated in Fig. 2A; it lasted 13 min and occurred, according to the video, during the daylight period while the animal was in a sleep-like position without any evidence of behavioral seizure. Full description of the $\dot{V}O_2$ pattern during seizures is provided in Supplementary Information S1.

Body position during SDOCs was different from the freezing posture frequently observed before rats have a behavioral seizure (stages 3 to 5 on Racine's scale). Notably, rats always exhibited a quiet state during SDOC events. The proportion of EPI/SDOC+ rats increased between weeks 3 and 5 post-SE and then stabilized (from weeks 5 to 8 post-SE). By the end of the 5th week post-SE, 30% (6/20) of EPI rats exhibited SDOCs (Fig. 3C).

**Ventilatory and cardiac functions in epileptic rats.** Experiment 2 aimed to further investigate the respiratory function by evaluating the ventilatory and cardiac functions of EPI rats during the chronic phase of epilepsy. It should be noted that we could not record respiratory variables using whole-body plethysmography before the 10th week post-SE because rats were not able to remain quiet in the apparatus chamber, likely owing to enhanced reactivity to stressful situations.

*EPI rats presented with a low steady respiratory pattern and enhanced ventilatory response to hypoxia and hypercapnia.* Under normoxia, EPI rats (n = 16) presented a slow respiratory pattern with a significant decrease in minute ventilation ($\dot{V}_E$), compared to controls (n = 5) (-12 ± 3%, P = 0.03) (See supplementary Fig. S2A online). The reduction of $\dot{V}_E$ in EPI rats resulted from a decrease in respiratory frequency ($f_R$) (-12 ± 2% vs.
<figure><img src="image_7.png" /><figcaption>B. Rat undergoing stage-5 seizure</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 2.** SDOC can occur independently of behavioral spontaneous seizures. (A) Example of an SDOC event recorded in an EPI rat at 5 weeks post-SE, when it was adopting a sleep posture. This SDOC event occurred during the light period and was identified with the decline in VO₂ values. No behavioral seizure was observed before, during, and after SDOC occurrence. In contrast, (B) a sharp increase in VO₂ was observed after a 43-s stage-5 seizure according to Racine's scale in an EPI/SDOC+ rat with VO₂ levels over the pre-ictal level for approximately 30 min.</figcaption></figure>

controls, P = 0.001), while tidal volume (V~T~) was similar between EPI and control rats (P = 0.92). The decrease in f~R~ in EPI rats was related to a simultaneous increase in inspiratory time (T~I~) and expiratory time (T~E~), which reached +20 ± 3% (P = 0.01) and +10 ± 2% (P = 0.001) compared to the respective values measured in con-
A. Time course of the daily $\dot{V}O_2$ in experimental groups

<table><thead><tr><th rowspan="2">Weeks post-SE</th><th colspan="3">Daily $\dot{V}O_2$ (mL x min<sup>-1</sup> x kg<sup>-1</sup>)</th></tr><tr><th>Control rats (n=8)</th><th>EPI/SDOC- rats (n=14)</th><th>EPI/SDOC+ rats (n=6)</th></tr></thead><tbody><tr><td>8</td><td>19.2 ± 0.6 *</td><td>18.0 ± 0.5 *</td><td>8.1 ± 0.3 **</td></tr><tr><td>7</td><td>18.1 ± 0.4 *</td><td>19.8 ± 0.5</td><td>10.1 ± 1.6 *</td></tr><tr><td>6</td><td>20.0 ± 0.4</td><td>19.9 ± 0.8</td><td>11.0 ± 2.5 <i>p</i>=0.09</td></tr><tr><td>5</td><td>19.8 ± 0.1</td><td>22.3 ± 1.2</td><td>12.5 ± 1.6 *</td></tr><tr><td>4</td><td>20.8 ± 0.5</td><td>20.1 ± 0.6</td><td></td></tr><tr><td>3</td><td>21.2 ± 0.9</td><td>21.0 ± 0.8</td><td></td></tr><tr><td>1-2</td><td>24.4 ± 1.1</td><td>23.1 ± 1.7</td><td>26.6 ± 1.0</td></tr></tbody></table>

B. Time course of the daily $\dot{V}O_2$ in EPI/SDOC+ rats

<table><thead><tr><th rowspan="2">Weeks post-SE</th><th colspan="2">Rats identified as EPI/SDOC+</th><th colspan="2">Rats not yet identified as EPI/SDOC+</th></tr><tr><th>n</th><th>Daily $\dot{V}O_2$ (mL x min<sup>-1</sup> x kg<sup>-1</sup>)</th><th>n</th><th>Daily $\dot{V}O_2$ (mL x min<sup>-1</sup> x kg<sup>-1</sup>)</th></tr></thead><tbody><tr><td>4</td><td>3/6</td><td>13.5 ± 3.5</td><td>3/6</td><td>20.1 ± 1.8</td></tr><tr><td>3</td><td>1/6</td><td>17.2</td><td>5/6</td><td>21.1 ± 1.9</td></tr></tbody></table>

C. Time course of SDOC appearance by EPI rats

<figure><img src="image_9.png" /><figcaption>**Figure 3.** Time course of daily $\dot{V}O_2$ pattern. Respiratory thermochemistry was recorded once a week for each rat, from the 1st to the 8th week post-SE. (A) Daily $\dot{V}O_2$ was calculated for each group (control rats (n=8), EPI/SDOC- rats (n=14), EPI/SDOC+ rats (n=6)) from 1 to 8 weeks post-SE. (B) Daily $\dot{V}O_2$ of animals was calculated at 3-4 weeks post-SE in animals identified as EPI/SDOC+ rats or at 1-2 weeks post-SE in animals pending identification as EPI/SDOC+ rats. Paired comparison of daily $\dot{V}O_2$ over weeks showed statistical significance with only 1-2 weeks post-SE. Results are expressed as mean ± SEM of $\dot{V}O_2$, *P* < 0.05; **P** = 0.01 in comparison to daily $\dot{V}O_2$ at 1-2 weeks post-SE. (C) The first SDOC was detected in one rat (1/20) during the 3rd post-SE week. The maximum number of EPI rats exhibiting SDOCs (EPI/SDOC+, 30%) was reached by the end of the 5th week post-SE. The proportion of EPI/SDOC+ rats was expressed in percent of the total number of rats developing epilepsy (yellow curve) or in percent of the total number of rats presenting with SDOCs by the end of the 8th week post-SE (red curve).</figcaption></figure>
trols. Consequently, total respiratory cycle duration ($T_{tot}$) was significantly higher in EPI rats than in controls (+14 ± 2%; P = 0.002) (See supplementary Fig. S2A online).

Regarding the stability of the breathing pattern, the coefficient of variation of breathing, apnea frequency, and duration were similar between EPI and control rats (P = 0.16, P = 0.61, and P = 0.09, respectively) (see supplementary Fig. S2B online).

Compared to control rats (n = 13), EPI rats (n = 18) showed greater hypoxia ventilatory response (HVR) with higher $\dot{V}_E$ (+23 ± 6% vs controls, P = 0.007), higher $f_R$ (+16 ± 2% vs controls, P < 0.001) and similar $V_T$ (P = 0.06). Similarly, the hypercapnia ventilatory response (HCVR) was higher in EPI rats (n = 18) than in controls (n = 12) with greater $\dot{V}_E$ (+14 ± 4% vs controls, P = 0.01) and greater $f_R$ (+10 ± 3% vs controls, P = 0.02), while $V_T$ remained similar to that of controls (P = 0.09).

**Abnormal $\dot{V}O_2$ pattern is associated with pronounced decrease in ventilatory rhythm at rest.** Since oxygen is strictly provided during the ventilatory process, we aimed to compare ventilation between EPI rats with or without abnormal $\dot{V}O_2$ patterns. In experiment 2, we were able to discriminate eight EPI/SDOC+ rats, which presented with lower daily $\dot{V}O_2$ compared to that of both EPI/SDOC- rats (n = 15, -51% ± 9%; P < 0.001) and control rats (n = 16, -41% ± 11%; P = 0.003).

Under normoxia, ventilation of EPI/SDOC- rats was significantly altered compared to that of control rats (Table 1A), but to a lesser extent than that of EPI/SDOC+ rats. Indeed, compared to control rats, EPI/SDOC- rats presented a slight but not significant 6 ± 3% of reduction in $\dot{V}_E$, associated with an 8 ± 2% of reduction (P = 0.015) in $f_R$, a 13 ± 1% of increase (P < 0.001) in $T_I$ and a 9 ± 2% of increase (P = 0.025) in $T_{tot}$ (Table 1A). Furthermore, compared to EPI/SDOC- rats, EPI/SDOC+ rats presented an additional 13 ± 3% of reduction (P = 0.013) in $\dot{V}_E$, that resulted from an 8 ± 2% of reduction (P = 0.036) in $f_R$, a 13 ± 3% of increase (P < 0.001) in $T_I$ and an 8 ± 3% of increase (P = 0.013) in $T_{tot}$ (Table 1A). Tidal volume ($V_T$), the coefficient of variation of the breathing frequency, apnea frequency, and duration were similar between all groups (Table 1A). Under hypoxic and hypercapnic conditions, $\dot{V}_E$, $V_T$, and $f_R$ were similar between EPI/SDOC+ (n = 6) and EPI/SDOC- rats (n = 12) (Table 1B).

**No alteration in cardiac frequency and body temperature in EPI rats.** The respiratory system, cardiac function, and body temperature are strongly related to each other. In this experiment, we showed that cardiac frequency and body temperature were similar between EPI/SDOC+ rats (n = 2), EPI/SDOC- rats (n = 3), and control rats (n = 3) during both the light and dark periods (Table 1C,D).

**Brainstem 5-HT system of EPI rats.** In experiment 3, 20 rats developed epilepsy. Among these rats, we identified eight EPI/SDOC+ and eight EPI/SDOC- rats. Eight control rats were used in this study. The brainstem was processed for either histological (n = 3/8 rats in each group) or gene expression (n = 5/8 rats in each group) analyses.

**Alterations in the density of brainstem 5-HT immunopositive processes in EPI rats.** We determined whether metabolic and ventilatory alterations developed by EPI/SDOC+ rats were associated with dysfunction of the brainstem 5-HT system. Therefore, we focused our investigation on caudal raphe nuclei, which are known to innervate neuronal groups closely implicated in the regulation of respiratory rhythm generation, such as the NTS and the VLM<sup>14,21</sup>. 5-HT-immunolabeling was detected in both processes and cell bodies. Within the caudal raphe nuclei (pallidus, magnus, and obscurus), although we observed a downward trend in epileptic rats, there was no statistical difference in the number of 5-HT-positive neuronal cell bodies between the three groups of rats (control rats: 22 ± 6, EPI/SDOC-: 18 ± 4, EPI/SDOC+: 13 ± 5, P = 0.339) (See supplementary Fig. S3 online). In addition, there was no difference in the surface area of 5-HT positive cell bodies (control rats: 383 ± 21 pixels<sup>2</sup>, EPI/SDOC- rats: 445 ± 37 pixels<sup>2</sup>, EPI/SDOC+ rats: 397 ± 54 pixels<sup>2</sup>, P = 0.383), or in the index of 5-HT neuronal concentration (control rats: 109 ± 19 A.U., EPI/SDOC- rats: 125 ± 23 A.U., EPI/SDOC+ rats: 119 ± 18 A.U., P = 0.864) between the three rat groups (See supplementary Fig. S3 online). The density of 5-HT immunopositive processes within both the VLM (Fig. 4A) and NTS (Fig. 4B) were significantly higher in EPI/SDOC- rats than in EPI/SDOC+ rats and control rats (Fig. 4A4 and 4B4). Furthermore, the density of 5-HT-immunolabeled processes within the NTS was significantly lower in EPI/SDOC+ rats compared to controls (Fig. 4B4).

**Alterations in brainstem transcript levels of TPH2 and SERT in EPI/SDOC+ rats.** TPH2 transcript levels were significantly increased in EPI/SDOC+ compared to both EPI/SDOC- rats and control rats. There was no statistical difference in the transcript levels of TPH2 between EPI/SDOC- and control rats (Fig. 5A). SERT transcript levels were significantly increased in EPI/SDOC+ compared to EPI/SDOC- rats, but not when compared to control rats. In addition, SERT transcript levels did not differ statistically between control and EPI/SDOC- rats (Fig. 5B).

**Brainstem 5-HT system in DBA/2 mice prone to AGS.** Previous studies demonstrated that almost 75% of DBA/2 mice developing AGS died from fatal respiratory arrest (FRA) following tonic hindlimb extension<sup>22</sup>. Among the remaining mice, 15% developed non-fatal respiratory arrest (NFRA), indicating that they spontaneously recovered from post-AGS respiratory arrest. The brainstem 5-HT system in DBA/2 mice exhibiting FRA in response to AGS was altered in comparison with C57Black/6 mice that did not exhibit AGS<sup>16</sup>. Thus, in experiment 4, we examined whether the difference in TPH2 and SERT gene expression that we observed between EPI/SDOC+ rats and EPI/SDOC- rats was consistent in FRA and NFRA DBA/2 mice after AGS. Out of 11 male DBA/2 mice that underwent AGS with respiratory arrest, 64% presented with FRA (n = 7) and 36% with NFRA
<table><caption>**Table 1.** Monitoring of cardiorespiratory functions in control and EPI rats between 12 and 15 weeks post-SE. Using plethysmography, ventilatory variables were measured under normoxic conditions (21% O₂) at 12 weeks post-SE (A) and in response to hypoxia (9% O₂) or hypercapnia (5% CO₂) exposure between 13–14 weeks post-SE (B). Cardiac frequency (C) and body-temperature (D) were concomitantly evaluated at 15 weeks post-SE. Fc cardiac frequency, fR respiratory frequency, T<sub>B</sub> body-temperature, T<sub>E</sub> expiration time, T<sub>I</sub> inspiration time, T<sub>tot</sub> total respiratory cycle duration, V̇<sub>E</sub> minute volume, V<sub>T</sub> tidal volume. Statistics: NT not statistically significant.</caption><thead><tr><th></th><th>Control (n=5)</th><th>EPI/SDOC- (n=8)</th><th>EPI/SDOC+ (n=8)</th><th>Ctrl vs SDOC- (P value)</th><th>Ctrl vs SDOC+ (P value)</th><th>SDOC+ vs SDOC- (P value)</th><th></th></tr></thead><tbody><tr><td colspan="8">(A) Basal respiratory variables under normoxia at 12 weeks post-SE</td></tr><tr><td>V̇<sub>E</sub> (mL min<sup>-1</sup> × 100 g<sup>-1</sup>)</td><td>32 ± 1</td><td>30 ± 1</td><td>26 ± 1</td><td>0.274</td><td>0.002</td><td>0.013</td><td></td></tr><tr><td>fR (breath × min<sup>-1</sup>)</td><td>101 ± 3</td><td>93 ± 2</td><td>87 ± 2</td><td>0.015</td><td>&lt;0.001</td><td>0.036</td><td></td></tr><tr><td>VT (mL × 100 g<sup>-1</sup>)</td><td>0.32 ± 0.01</td><td>0.33 ± 0.01</td><td>0.31 ± 0.01</td><td>NT</td><td>NT</td><td>NT</td><td></td></tr><tr><td>TI (ms)</td><td>207 ± 7</td><td>233 ± 3</td><td>263 ± 8</td><td>0.015</td><td>&lt;0.001</td><td>&lt;0.001</td><td></td></tr><tr><td>TE (ms)</td><td>385 ± 12</td><td>413 ± 9</td><td>438 ± 10</td><td>0.088</td><td>0.003</td><td>0.083</td><td></td></tr><tr><td>Ttot (ms)</td><td>593 ± 14</td><td>646 ± 11</td><td>700 ± 17</td><td>0.025</td><td>&lt;0.001</td><td>0.013</td><td></td></tr><tr><td>Apnea frequency × min<sup>-1</sup></td><td>0.32 ± 0.17</td><td>0.27 ± 0.13</td><td>0.10 ± 0.07</td><td>NT</td><td>NT</td><td>NT</td><td></td></tr><tr><td>Apnea duration (ms)</td><td>1339 ± 67</td><td>1435 ± 65</td><td>1477 ± 75</td><td>NT</td><td>NT</td><td>NT</td><td></td></tr><tr><td>Coefficient of variation</td><td>0.49 ± 0.01</td><td>0.53 ± 0.01</td><td>0.52 ± 0.02</td><td>NT</td><td>NT</td><td>NT</td><td></td></tr><tr><td>% of baseline</td><td>Control (n = 13)</td><td>EPI/SDOC- (n = 12)</td><td>EPI/SDOC+ (n = 6)</td><td>Ctrl vs SDOC- (P value)</td><td>Ctrl vs SDOC+ (P value)</td><td>SDOC+ vs SDOC- (P value)</td><td></td></tr><tr><td colspan="8">(B) Respiratory response to hypoxia and hypercapnia at 13-14 weeks post-SE</td></tr><tr><td rowspan="3">Hypoxia 9% of O2</td><td>V̇<sub>E</sub> (mL min<sup>-1</sup> × 100g<sup>-1</sup>)</td><td>225 ± 9</td><td>273 ± 18</td><td>278 ± 19</td><td>0.019</td><td>0.046</td><td>0.847</td><td></td></tr><tr><td>fR (breath × min<sup>-1</sup>)</td><td>146 ± 5</td><td>171 ± 4</td><td>165 ± 3</td><td>&lt;0.001</td><td>0.024</td><td>0.463</td><td></td></tr><tr><td>VT (mL × 100 g<sup>-1</sup>)</td><td>153 ± 5</td><td>169 ± 8</td><td>170 ± 7</td><td>NT</td><td>NT</td><td>NT</td><td></td></tr><tr><td rowspan="3">Hypercapnia 5% of CO2</td><td>V̇<sub>E</sub> (mL min<sup>-1</sup> × 100g<sup>-1</sup>)</td><td>187 ± 6</td><td>211 ± 8</td><td>219 ± 14</td><td>0.028</td><td>0.023</td><td>0.562</td><td></td></tr><tr><td>fR (breath × min<sup>-1</sup>)</td><td>128 ± 4</td><td>141 ± 4</td><td>141 ± 8</td><td>NT</td><td>NT</td><td>NT</td><td></td></tr><tr><td>VT (mL × 100 g<sup>-1</sup>)</td><td>141 ± 3</td><td>148 ± 5</td><td>153 ± 6</td><td>NT</td><td>NT</td><td>NT</td><td></td></tr><tr><td colspan="2"></td><td>Control (n = 3)</td><td>EPI/SDOC- (n = 3)</td><td>EPI/SDOC+ (n = 2)</td><td>Ctrl vs SDOC- (P value)</td><td>Ctrl vs SDOC+ (P value)</td><td>SDOC+ vs SDOC- (P value)</td><td></td></tr><tr><td colspan="8">(C) Cardiac frequency at 15 weeks post-SE</td></tr><tr><td>Day</td><td>FC (beats × min<sup>-1</sup>)</td><td>315 ± 17</td><td>338 ± 11</td><td>333 ± 8</td><td>NT</td><td>NT</td><td>NT</td><td></td></tr><tr><td>Night</td><td>FC (beats × min<sup>-1</sup>)</td><td>332 ± 19</td><td>355 ± 9</td><td>351 ± 12</td><td>NT</td><td>NT</td><td>NT</td><td></td></tr><tr><td colspan="2"></td><td>Control (n = 3)</td><td>EPI/SDOC- (n = 3)</td><td>EPI/SDOC+ (n = 2)</td><td>Ctrl vs SDOC- (P value)</td><td>Ctrl vs SDOC+ (P value)</td><td>SDOC+ vs SDOC- (P value)</td><td></td></tr><tr><td colspan="8">D- Body-temperature at 15 weeks post-SE</td></tr><tr><td>Day</td><td>TB (°C)</td><td>37.5 ± 0.1</td><td>37.6 ± 0.1</td><td>37.7 ± 0.2</td><td>NT</td><td>NT</td><td>NT</td><td></td></tr><tr><td>Night</td><td>TB (°C)</td><td>38 ± 0.2</td><td>37.9 ± 0.1</td><td>38.1 ± 0.1</td><td>NT</td><td>NT</td><td>NT</td><td></td></tr></tbody></table>

(n = 4). We found that the levels of TPH2 and SERT transcripts were significantly elevated in FRA DBA/2 mice compared to those in NFRA DBA/2 mice (Fig. 5A,B). Strikingly, the levels of TPH2 and SERT transcripts in FRA DBA/2 mice followed a trend similar to that of EPI/SDOC+ rats when both were compared to their respective controls (NFRA DBA/2 mice and EPI/SDOC- rats, respectively).

Discussion

In this study, we report chronic alterations in the respiratory system in epileptic (EPI) rats associated with 5-HT dysfunction within the brainstem.

**Exploration of respiratory function in epilepsy models.** To date, studies in animal models of seizures, aimed at understanding respiratory dysfunction in the context of epilepsy, were mostly performed during the period surrounding seizures, i.e., during ictal and post-ictal periods²²⁻³². Respiratory alterations were detected in these studies during short recording periods (2–30 min). In view of the respiratory disorders that can affect patients suffering from temporal lobe epilepsy¹,²,⁵,⁶, our work, which was performed on a widely used preclinical model of this type of epilepsy (status epilepticus induced by intraperitoneal administration of pilocarpine, Pilo-SE rats), had two major objectives.

The first was to study the respiratory function in rats with a long history of epilepsy and at distance of brain remodeling period that follows the severe brain damage induced by Pilo-SE³³. After a 10-week period of active
<figure><img src="image_10.png" /><figcaption>B. Level expression of mRNA encoding SERT</figcaption></figure>

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

**Figure 5.** TPH2 and SERT mRNA transcript levels were measured in the brainstem of EPI rats and DBA/2 mice. Changes in brainstem levels of TPH2 (A) and SERT (B) transcripts were measured in EPI rats and DBA/2 mice at 8 weeks post-SE and 2–3 min after audiogenic seizure (AGS), respectively. In contrast to EPI rats, DBA/2 mice had no history of seizure before AGS. Some mice developed fatal respiratory arrest (DBA/2 with FRA) or non-fatal respiratory arrest (DBA/2 with NFRA) after AGS. TPH2 and SERT transcript levels in EPI/SDOC+rats and in DBA/2 mice with FRA are expressed as percentage (mean ± SEM) of EPI/SDOC– rats and DBA/2 mice with NFRA, respectively (n = 5 in each group). †P < 0.05, ††: P < 0.01, †††: P < 0.001, in comparison to EPI/SDOC– rats and DBA/2 mice with FRA. *P < 0.05; **P < 0.01, in comparison to control rats.

epilepsy (12 weeks post-SE), we observed a decrease in breathing rhythm in epileptic (EPI) rats at rest. This result is different from that reported in a recent study showing no alteration in resting ventilation in rats 2–4 weeks after Pilo-SE³⁴, probably due to the fact that we did not use the same strain of rats (Wistar vs. Sprague-Dawley in our study), or the same route of pilocarpine administration (intrahippocampally vs. intraperitoneally in our study). Nonetheless, we agree with the authors that the mechanisms leading to alterations in breathing may depend on the stage of epilepsy; in their study, the assessment was performed at an earlier stage compared to that in our study (chronic epilepsy stage).

The second objective was to perform respiratory function recording over a long duration (24 h) in order to extend the investigation beyond the period surrounding seizures. Using indirect respiratory thermochemistry, we performed a 24-h monitoring to determine the amount of oxygen consumed by the whole body. Subsequently, we provided evidence of the presence of abnormal VO₂ patterns in 30–50% of epileptic rats, which started with the installation of epilepsy (2–3 weeks post-SE). These EPI rats (EPI/SDOC + group) presented chronic VO₂ abnormalities characterized by low daily patterns, with intermittent and repetitive SDOC events.

Until now, both preclinical and clinical studies have investigated changes in the metabolic rate of O₂ in epilepsy mostly at the brain level, reporting a decline in the activity of mitochondrial respiratory complexes in rats subjected to Pilo-SE and in patients with TLE³⁵,³⁶. In our study, we measured whole-body O₂ intake, which is mostly related to gas exchange during the ventilatory process³⁷. However, whole-body O₂ intake can also be affected directly or indirectly by several factors, including age, circadian rhythm, body weight, body temperature, cardiac function, and activity²⁰,³⁷⁻⁴². None of these factors were found to be associated with abnormal pattern of
V̇O₂ in our study. We then focused on ventilatory variables that may explain abnormal V̇O₂ patterns in a subset of EPI rats.

**Respiratory system at rest in a subset of EPI rats.** The identification of heterogeneous phenotypes within a single disease is a common feature in human epilepsy, but has rarely been described in preclinical studies. One study dichotomized two groups of rats subjected to kainic acid-induced SE according to the extent of brain damage and brain metabolic alterations⁴³. In our study, the EPI/SDOC+ group presented with a decrease in daily V̇O₂ and a severe decrease in V̇E compared to the EPI/SDOC- group.

During a 20-min monitoring period at rest, the decrease in both V̇E and V̇O₂ has already been described as an inherent dysfunction in a genetic model of audiogenic seizure in the WAR rat strain³¹. Because of technical limitations (24-h V̇O₂ monitoring), we could not jointly assess V̇E and V̇O₂ in our study, which prevented us from interrogating the potential direct relationship between V̇E and V̇O₂. Nevertheless, we estimated the V̇E/V̇O₂ ratio for each rat, as an indicator of ventilatory efficiency; this ratio was significantly higher in EPI/SDOC+ rats compared to EPI/SDOC- (+52%, P = 0.007) and healthy control rats (+35%, P = 0.04), and not statistically different between the latter rat groups.

High V̇E/V̇O₂ has been described as an index of poor gas exchange during physical exercise in individuals with heart failure disease⁴⁴. Further exploration using blood oxygen measurements are required to determine the direct relationship between V̇E and V̇O₂ in epileptic rats following Pilo-SE. Another limitation of this study was the short duration of the ventilatory recording compared with the 24 h of V̇O₂ assessment. A significant fluctuation of V̇E could have possibly been observed, similar to those observed for daily V̇O₂, if the recording of V̇E had been longer. Further studies with devices that can concomitantly measure both V̇E and V̇O₂ over a long period of time in free-moving rats will help to establish a direct link between V̇O₂ and V̇E profiles and more specifically during SDOC events.

Indeed, these events could represent a disruption in the O₂ uptake balance. To our knowledge, events similar to SDOCs have never been reported, and could only be interpreted in light of our knowledge regarding dyspneic patterns. One possible explanation for SDOCs could be the presence of ventilatory pauses; however, they have to be extremely numerous and with a long duration to induce such a drop in V̇O₂. Another explanation could be a combination of several factors; for example, the lack of gasping, especially in hypoxic conditions, could lead to poor tissue oxygenation⁴⁵. In EPI/SDOC+ rats, V̇E recording was too short to determine whether the absence of gasping could be part of their ventilatory dysfunction.

**Ventilatory regulation in EPI rats.** We observed an increased HVR in EPI rats, in contrast to some studies that reported no alteration in HVR in Pilo-SE rats¹⁰,³⁴. In addition, EPI/SDOC+ and EPI/SDOC- rats exhibited sustained HCVR, whereas previous studies have reported a diminished response in Pilo-SE-induced epilepsy, rapid amygdala kindling model, and WAR rat model³¹,³⁴,⁴⁶. It is very difficult to compare our data with those in the literature owing to differences in the models used, timing of respiratory recordings, and severity of exposure to hypoxia and hypercapnia. Given that low HVR and HCVR are generally indicative of alterations in peripheral or central respiratory components, their maintenance in high levels in our study may indicate that if alterations occurred in early phases, they were compensated for, or repaired, in the late stages of epilepsy at which we performed our investigation.

**Brainstem 5-HT alterations.** 5-HT is mostly described as an excitatory drive towards the respiratory network, and subsequently the motor output²¹. Dysregulation of brainstem neuronal groups involved in the control of breathing, particularly those of the 5-HT system, has been found in patients with drug-resistant epilepsy and has been associated with an increased risk of SUDEP¹³. We found that brainstem 5-HT innervation was impaired in EPI rats compared to control rats, with notable differences between EPI/SDOC+ and EPI/SDOC- rats. Indeed, EPI/SDOC- and EPI/SDOC+ rats presented an opposite 5-HT innervation pattern within the NTS. Since the impairment of respiratory function seems to be moderate in EPI/SDOC- rats, in contrast to EPI/SDOC+ rats, it is conceivable that the increase in 5-HT innervation in both the NTS and VLM may be part of the protective mechanisms of the respiratory function. However, EPI/SDOC+ rats presented an increase in brainstem levels of both TPH2 and SERT transcripts. In EPI/SDOC+ rats, such a mechanism may support an increase in 5-HT concentration in neuron terminals by increasing 5-HT synthesis and re-uptake of intact 5-HT molecules, and thus may represent an adaptive mechanism to partially compensate for the decreased brainstem 5-HT innervation. This conclusion is supported by results reported in conditional knock-out Lmx1b mice and WAR rats, which present respiratory alterations at rest and a concomitant decrease in 5-HT drive³¹,⁴⁷, similar to that seen in EPI/SDOC+ rats. We also believe that according to the difference in 5-HT innervation between VLM and NTS observed in the same EPI/SDOC+ rats, the decrease in 5-HT could be structure-specific and requires further investigation in other primary regions of respiratory regulation.

We also compared the brainstem 5-HT system of DBA/2 mice with FRA to relevant control mice—DBA/2 mice with NFRA, in contrast to previous studies that compared DBA/2 mice with FRA to C57BL/6 mice that did not exhibit AGS¹⁶. We found an increased level of transcripts for TPH2 and SERT in DBA/2 mice with FRA. Surprisingly, 5-HT profile of EPI/SDOC+ rats was similar to that of DBA/2 mice with FRA. Accordingly, alterations in the 5-HT system are likely involved in respiratory dysfunction in both acute and chronic seizure models and more likely associated to more severe respiratory phenotype in both strains.

In conclusion, the Pilo-SE rat model, which shows both chronic respiratory alterations and spontaneous seizures, could be a relevant model to further explore the neurobiological mechanisms underlying chronic alterations in respiratory regulation in epilepsy. Identifying the mechanisms that protect EPI/SDOC+ rats from mortality may also be an interesting approach for new therapeutic strategies to prevent SUDEP.
**Methods**

Protocols and methods in this study were performed in compliance with relevant institutional and ARRIVE guidelines and regulations. All experiments have been approved by the ethical committee of the Claude Bernard Lyon 1 University (protocol # BH-2008-11) and done in accordance with the European Communities Council for animal care (directive 2010-63), according to the French law (decree 2013-118).

**Animals.** Five-week-old male Sprague-Dawley rats (Envigo, The Netherlands) were delivered to the animal facility of the Lyon 1 University (Villeurbanne, France) and housed in groups of five in transparent polycarbonate cages. Male DBA/2 mice of breeding colony of the LTSI laboratory were housed in groups of five. Rats and mice were maintained under a controlled 12-h/12-h light–dark cycle in a temperature-controlled room (24–26 °C and 21–23 °C for rats and mice, respectively), with free access to food and water.

**Experimental design.** A full description of methods has been described in detail elsewhere<sup>17,33,48–50</sup>. Rats subjected to Pilo-SE were observed daily for spontaneous seizures, until development of epilepsy (EPI rats). The respiratory thermochemistry setup was used to measure oxygen consumption ($\dot{V}$O₂) during 24 consecutive hours (daily $\dot{V}$O₂) for each rat. Animals in this study were subjected to one of the following experiments (See supplementary Fig. S4 online):

*Experiment 1.* To establish the time course of respiratory alterations in Pilo-SE rats, 4 series with 10 rats each were used in this experiment. In each series, Pilo-SE was induced in eight rats, and two rats were used as controls. For each series, the daily $\dot{V}$O₂ for each rat was measured once a week from the 1st to the 8th week post-SE.

*Experiment 2.* To further investigate cardiorespiratory functions in EPI rats during the chronic phase of epilepsy, Pilo-SE was triggered in 30 rats, while 16 rats were used as controls. During the chronic period of epilepsy (12 weeks post-SE), resting ventilation (normoxia, 21% O₂) was measured using whole-body plethysmography. Between 13–14 weeks post-SE, a subset of EPI and control rats were used to evaluate ventilatory responses to hypoxia (HVR, FiO₂=9%) and hypercapnia (HCVR, FiCO₂=5%), while the others underwent implantation of a telemetric device to evaluate cardiac frequency and body temperature at 15 weeks post-SE. As $\dot{V}$O₂ and ventilation are intermingled processes, to visualize the global state of the respiratory system in EPI rats, we performed an in-depth analysis of our data by characterizing ventilatory patterns according to $\dot{V}$O₂ patterns. Notably, daily $\dot{V}$O₂ in EPI and control rats was measured shortly before plethysmography at 7–9 weeks post-SE.

*Experiment 3.* To determine whether the integrity of the brainstem 5-HT system is altered in EPI rats with respiratory alterations, Pilo-SE was induced in 20 rats, while eight rats were used as controls. Daily $\dot{V}$O₂ was evaluated between the 6th and 8th weeks post-SE. Subsequently, animals were transcardially perfused with paraformaldehyde or NaCl and their brains were collected and processed for histology (5-HT immunodetection; n=3 in each group), or for gene expression analysis (SERT and TPH2 mRNA; n=5 in each group), respectively.

*Experiment 4.* To determine whether the integrity of the brainstem 5-HT system is compromised in DBA/2 mice with FRA, 11 DBA/2 mice underwent AGS followed by respiratory arrest. Brains from both FRA and NFRA DBA/2 mice were removed immediately after AGS. As for NFRA mice, mice with FRA were rescued and were injected with a lethal dose of pentobarbital and then perfused with NaCl. Brainstem SERT and TPH2 gene expression was measured and compared between mice with FRA and NFRA.

**In vivo procedures.** *Pilocarpine-induced SE.* Epilepsy was triggered following Pilo-SE in 180–200 g male rats at the age of 7 weeks, as previously described<sup>33,49</sup>. Scopolamine methylnitrate (1 mg×kg⁻¹, s.c.; S-2250, Sigma-Aldrich) was administered 30 min prior to pilocarpine hydrochloride (350 mg×kg⁻¹, i.p.; P-6503, Sigma-Aldrich). After an initial period of immobility, the onset of SE was characterized by repetitive tonic–clonic activity of the trunk and the limbs, occurring after repeated rearing with forelimb clonus and falling. Because SE event is a life-threatening condition in rats, there is a need to stop it after 2 h by a single injection of diazepam (Valium®, 10 mg×kg⁻¹, i.p.; Roche). Rats were then hydrated with 2 mL of saline solution (0.9% NaCl; s.c.) and transferred to individual cages. After a one-week recovery period, all rats (control and treated rats) were housed in groups of 5 to avoid any stress that may result from social isolation.

**Animal care after Pilo-SE.** Control and treated rats were weighted daily during the first two weeks following Pilo-SE, and then every week until termination of the experiment. Daily abdominal massage was performed twice a day during the first week to activate intestinal motility, which is disrupted following Pilo-SE.

**Detection of spontaneous recurrent seizures after Pilo-SE.** To maintain as low as possible the level of anxiety in rats subjected to Pilo-SE, electro-encephalographic recordings (EEG) were not performed to detect the onset of epilepsy in our experiments. Indeed, conventional EEG set-up requires that rats be housed individually in a restrictive environment, which is known to enhance stress and anxiety and worsens epilepsy phenotype in Pilo-SE models<sup>51</sup>. In addition, as mentioned below in the “whole-body plethysmography” section, it is difficult to perform ventilatory recordings in animals with high anxiety level.

It is now clearly established that rats of the Sprague-Dawley strain from the Janvier Laboratories develop spontaneous recurrent seizures between 10 and 36 days after induction of Pilo-SE at 42 days of age<sup>52,53</sup>. Because we were using Sprague-Dawley rats from Envigo Laboratories, we had to monitor the time of onset of spontaneous
recurrent seizures (SRSs) in our experiments. To determine whether rats were in the chronic phase of epilepsy, they were manipulated daily between 09:00 and 11:00 a.m. to define the onset of handling-induced seizures (HIS), which was induced by restraining rats at the level of the chest with gentle pressure for 10 s⁹,⁵². HIS were used in our study to help experimenters to define the time to start detecting the first behavioural spontaneous recurrent seizures (SRSs). Thus, upon detection of 2 HIS, rats were observed by experimenters for 5 consecutive hours (between 01:00 and 06:00 p.m.) over the following days to detect the presence or absence of behavioural SRSs. Following Pilo-SE, rats are considered epileptic once they developed at least two SRSs during the periods of the daily observations. The severity of these seizures was scored from 3 to 5 on the Racine's scale⁵⁴. Seizure severity of less than 3 was difficult to detect visually in this study, which was conducted without daily video monitoring. Stage 3 corresponds to forelimb clonus, stage 4 to bilateral forelimb clonus with rearing, and stage 5 to rearing and falling or to tonic-clonic seizures.

During the 24 h-VO₂ recordings, rats were video-monitored for 12 h during the light period, making it possible to evaluate the number of seizures and their severity per session.

**Evaluation of respiratory and cardiac functions.** *Respiratory thermochemistry.* To detect respiratory alterations in epileptic rats, we used respiratory thermochemistry as an alternative approach to plethysmography, allowing respiratory function to be recorded for a long period of time on fed, unrestrained, and unanesthetized animal. Respiratory thermochemistry mostly reflects gas exchanges and allows more specifically the measurement of the oxygen consumption (VO₂) with no bias effect from animal's body movements on the accuracy of respiratory acquisition.

Oxygen consumption (VO₂) was measured with sample rate of 3 acquisitions per minute, under standard conditions (variation in temperature = 0 °C; Ambient pressure = 760 mm Hg; dry air) using a thermochemistry open system, as previously described⁵⁵. VO₂ was continuously recorded during 24 consecutive hours on the same thermochemistry system; therefore, only one rat could be recorded each day in our experimental setting. Briefly, rats were placed individually in a thermostatic recording chamber (12 L) with steady airflow (4 L × min⁻¹), under a controlled 12-h/12-h light–dark cycle (light from 06:00 a.m. to 06:00 p.m.), and were allowed free access to food and water. Ambient temperature inside the thermostatic chamber was monitored with copper-constantan thermocouples and was maintained at 24–26 °C using a thermostatic container. Airflow rate was measured using a Platon volumeter and converted to standard STPD (standard temperature and pressure, dry) values. Fractional O₂ concentration was measured using a Servomex 1100 paramagnetic gas analyzer (Taylor Instrument Analytics Ltd, Sussex, UK) and VO₂ was calculated. The O₂ analyzer was calibrated daily with pure nitrogen gas (0% O₂) and room air (20.93% O₂). Rats were weighted before and after recording then VO₂ values were normalized to the average body weight (calculated from the weight measured before and after each recording session). Rats were video-monitored during the daylight period.

*Whole-body plethysmography.* Plethysmographic chamber allows collecting ventilatory variables from acquired breathing pattern signal. Ventilatory function was evaluated between 09:00 a.m. and 04:00 p.m. in unrestrained, unanesthetized rats under BTPS (body temperature and pressure, saturated) conditions using whole-body plethysmography (EMKA; technology). Calibration of the system was performed before each recording session by injecting 1 mL air pulse inside the recording chamber. Briefly, each rat was individually placed in the recording chamber (600 mL) that is connected to a constant bias flow supply providing continuous inflow of fresh air (2.5 L × min⁻¹). Minute volume (V̇E; mL × min⁻¹) was calculated as the product of tidal volume (V̇T; mL) and respiratory frequency (ḟR; breath × min⁻¹). Total respiratory cycle duration (Ṫtot) was calculated as the sum of the inspiratory time (Ṫi; msec) and expiratory time (ṪE; msec). We also calculated for each animal the coefficient of variation of the respiratory frequency, which is an index of the variability of the respiratory pattern⁵⁶ and the frequency of spontaneous apnea, which is defined as a cessation of breathing lasting more than two respiratory cycles (>2 × Ṫtot). Apnea number and duration were calculated. Apnea frequency was normalized to recording duration and expressed per min of breathing. Recordings were started after a 40-min habituation period. Then, resting ventilation was monitored during 20 min under normoxic condition (21% O₂). After recordings, rats were brought back to their home cage.

Ventilatory responses to hypoxia (FiO₂=9%) and to hypercapnia (FiCO₂=5%) were evaluated during two distinct recording sessions, which were performed a week apart from each other. After a 40-min habituation period, plethysmographic recordings under basal condition (normoxia) were started 10 min before gas variations. Ventilatory responses to a 10-min exposure to hypoxia or hypercapnia were subsequently evaluated. All plethysmography recordings were made between 09:00 a.m. and 04:00 p.m.; because, beyond 04:00 p.m., EPI rats were more stressed in the plethysmography chamber, requiring longer time to reach a quiet state in this novel and restricted environment. At the end of each recording session, rats were brought back to their home cage. Repeated measurements are randomly performed for each rat.

*Abdominal implantation of telemetry transmitters to evaluate cardiac function and body-temperature.* Rats were anesthetized with a mix of ketamine and xylazine (75 mg × kg⁻¹ and 25 mg × kg⁻¹, respectively; i.p.). A mid-line abdominal incision was made and a sterile telemetry transmitter (Data Sciences International, France) was placed inside the abdominal cavity. A pair of biocompatible wire electrodes originating from the transmitter were tunnelled subcutaneously and implanted to the left intercostal muscles and between the right shoulder blades. The transmitter was attached to the abdominal wall during wound closure with absorbable sutures. After surgery, animals received an intramuscular injection of penicillin (15,000 units) once a day for 3 consecutive days and were then returned to their home cage. Rats were monitored daily to control general health state. Two weeks after surgery, telemetry transmitters were magnetically activated and freely moving animals were placed
in individual cages placed under platform-receiver units. Variation in cardiac frequency and body-temperature
were recorded during 24 consecutive hours.

**Induction of audiogenic seizures in DBA/2 mice.** AGS were induced as previously described<sup>17,22</sup>, with slight modifications. Briefly, at age of 23 days, mice were individually placed in a plastic cylinder (40 cm diameter, 20 cm height) and exposed to a broadband acoustic stimulus (intensity of 110 dB SPL with a peak frequency at 12 kHz) generated by a mechanical sonicator (Deltasonic model Delta 920). The stimulus was given for a maximum duration of 60 s or until the mouse exhibited AGS, characterized by an initial wild running and ending by a tonic hindlimb extension. All mice exhibited a tonic hindlimb extension leading to a respiratory arrest identified by pinnae relaxation and a complete cessation of chest movements. Respiratory arrest was followed by spontaneous recovery of the respiratory function within 5 s in mice identified as NFRA, or was followed by death for the other mice (FRA). Mouse brain was systematically removed within 2-3 min after acoustic stimulus was triggered.

**Ex vivo procedures.** DBA/2 mice and rats were injected with a lethal dose of pentobarbital (250 mg × kg⁻¹; i.p.) before being sacrificed. For gene expression analyses, animals were transcardially perfused (30 mL min⁻¹) for 2 min with ice-cold NaCl and their brain was rapidly removed. Then, the brainstem region was isolated, frozen in liquid nitrogen and stored at – 80 °C. For immunohistochemistry procedure, animals were transcardially perfused (30 mL min⁻¹) for 2 min with ice-cold NaCl and then for 10 min with ice-cold 4% paraformaldehyde in 0.2 M phosphate buffer (30 mL min⁻¹). After cryoprotection in a 25% sucrose solution, brains were frozen at – 40 °C in isopentane and stored at – 80 °C.

**5-HT immunohistochemistry.** Coronal sections (40 µm thick) were cryostat-cut throughout the medulla oblongata and were processed free-floating for 5-HT immunostaining. Using the rat brain atlas, sections from each experimental group were selected at interaural coordinates – 2.60 mm, and were all processed together for 5-HT immunofluorescent labeling. Selected sections were rinsed in 0.1 M Phosphate Buffered Saline (PBS) and incubated in PBS containing 0.3% Triton X-100 and 3% normal donkey serum. Then, sections were incubated overnight at 4 °C with a rabbit polyclonal antibody against 5-HT (Sigma; S5545; final dilution = 1:3000). After few washes, sections were incubated overnight at 4 °C with an AlexaFluor-488 conjugated donkey anti-rabbit antibody (Molecular Probes; A-21206; final dilution = 1:1000). At the end of the immunostaining protocol, sections were mounted onto SuperFrost®Plus slides and cover slipped using Vectashield mounting medium (H-1500, Vector).

*Semi-quantification of brainstem 5-HT immunostaining.* Images (1152 × 1024 pixels) were captured with a HLX PL APO×63/0.70 objective and with the same conditions of photomultiplier gain, offset and pinhole aperture using a TCS SP2 confocal microscopy system (Leica), allowing the comparison of fluorescence intensity between sections. Image analysis was conducted on 8 z-compressed images stacks (each image=4 µm depth), which include the whole region of interest. Using the image J software (NIH), we first measured the average grey value within each 5-HT-positive neuronal cell body of caudal raphe nuclei (pallidus, magnus and obscurus), as an index of 5-HT neuronal concentration (A.U.). Using the same software, we also measured the density of 5-HT-immunopositive processes within the surface area of the NTS and the VLM, expressed as 5-HT-positive pixel/pixels².

*Quantification of transcript level variation by RT-qPCR.* Total RNA was extracted from brainstem samples using Tri-reagent LS (Molecular Research Center, Inc.). Contaminant genomic DNA was subsequently removed from the samples by treatment with Turbo DNA-free™ kit (Ambion). The messenger RNAs (mRNAs) contained in 1 µg of purified RNA extracts were reverse-transcribed using the reverse transcriptase RNase H minus (Promega) and oligo-d(T)15. To normalize the RT step, a synthetic external and non-homologous poly(A) standard RNA (SmRNA; Morales and Bezin, patent WO2004.092414) was added to the RT reaction mix (150,000 copies in each experimental sample), as previously described<sup>33,48–50</sup>. PCR amplification of selected cDNAs, i.e. those complementary to the mRNA encoding rate-limiting enzyme of brainstem 5-HT synthesis (tryptophan hydroxylase: TPH2) and the 5-HT transporter (SERT), were performed using the Rotor-Gene Q system (Qiagen) and the QuantiTect SYBR Green PCR Kit (Qiagen). Sequences of the different primer pairs used for PCR amplification of mouse and rat TPH2 and SERT cDNAs are listed in Supplementary Table S1 online.

**Statistical analysis.** The data for each experimental group were pooled and expressed as the mean ± SEM in the text and figures, unless otherwise specified. All statistics were performed using Sigma plot software 3.5 (Systat Software Inc., Point Richmond, California, USA). A Kolmogorov–Smirnov test was used to test if the data follows a normal distribution, then paired or unpaired t-test was performed for the two-group comparison. For multiple group comparisons, one-way ANOVA followed by post-hoc Fisher's LSD test was used. Statistical significance was predefined as P < 0.05.

**Data availability**

The data that support the findings of this study are available from the corresponding author upon reasonable request.
Received: 5 November 2021; Accepted: 2 June 2022

Published online: 17 June 2022

## References

1. Lacuey, N. et al. The incidence and significance of periictal apnea in epileptic seizures. *Epilepsia* **59**, 573–582. https://doi.org/10.1111/epi.14006 (2018).

2. Vilella, L. et al. Incidence, recurrence, and risk factors for peri-ictal Central Apnea and sudden unexpected death in epilepsy. *Front Neurol.* **10**, 166. https://doi.org/10.3389/fneur.2019.00166 (2019).

3. Ryvlin, P. et al. Incidence and mechanisms of cardiorespiratory arrests in epilepsy monitoring units (MORTEMUS): A retrospective study. *Lancet. Neurol.* **12**, 966–977. https://doi.org/10.1016/S1474-4422(13)70214-X (2013).

4. Bourdillon, P., Rheims, S., Isnard, J. & Guenot, M. Letter to the Editor. Temporal lobe epilepsy: open or stereotactic surgery? *J. Neurosurg.* **131**, 989. https://doi.org/10.3171/2018.11.JNS183107 (2019).

5. Sivathamboo, S. et al. Sleep-disordered breathing in epilepsy: epidemiology, mechanisms, and treatment. *Sleep* **41**. https://doi.org/10.1093/sleep/zsy015 (2018).

6. Somboon, T., Grigg-Damberger, M. M. & Foldvary-Schaefer, N. Epilepsy and sleep-related breathing disturbances. *Chest* **156**, 172–181. https://doi.org/10.1016/j.chest.2019.01.016 (2019).

7. Levesque, M., Avoli, M. & Bernard, C. Animal models of temporal lobe epilepsy following systemic chemoconvulsant administration. *J. Neurosci. Methods* **260**, 45–52. https://doi.org/10.1016/j.jneumeth.2015.03.009 (2016).

8. Curia, G., Longo, D., Biagini, G., Jones, R. S. & Avoli, M. The pilocarpine model of temporal lobe epilepsy. *J. Neurosci. Methods* **172**, 143–157. https://doi.org/10.1016/j.jneumeth.2008.04.019 (2008).

9. Matovu, D. & Cavalheiro, E. A. Differences in evolution of epileptic seizures and topographical distribution of tissue damage in selected limbic structures between male and female rats submitted to the pilocarpine model. *Front. Neurol.* **13**, 802587. https://doi.org/10.3389/fneur.2022.802587 (2022).

10. Campos, R. R., Tolentino-Silva, F. R. & Mello, L. E. Respiratory pattern in a rat model of epilepsy. *Epilepsia* **44**, 712–717. https://doi.org/10.1046/j.1528-1157.2003.38402.x (2003).

11. Petrucci, A. N., Joyal, K. G., Purnell, B. S. & Buchanan, G. F. Serotonin and sudden unexpected death in epilepsy. *Exp. Neurol.* **325**, 113145. https://doi.org/10.1016/j.expneurol.2019.113145 (2020).

12. Zhan, Q. et al. Impaired serotonergic brainstem function during and after seizures. *J. Neurosci.* **36**, 2711–2722. https://doi.org/10.1523/JNEUROSCI.4331-15.2016 (2016).

13. Patodia, S. et al. The ventrolateral medulla and medullary raphe in sudden unexpected death in epilepsy. *Brain J. Neurol.* **141**, 1719–1733. https://doi.org/10.1093/brain/awy078 (2018).

14. Garcia, A. J., 3rd, Zanella, S., Koch, H., Doi, A. & Ramirez, J. M. Chapter 3--networks within networks: the neuronal control of breathing. *Prog. Brain Res.* **188**, 31–50. https://doi.org/10.1016/B978-0-444-53825-3.00008-5 (2011).

15. Faingold, C. L., Randall, M., Mhaskar, Y. & Uteshev, V. V. Differences in serotonin receptor expression in the brainstem may explain the differential ability of a serotonin agonist to block seizure-induced sudden death in DBA/2 vs. DBA/1 mice. *Brain Res.* **1418**, 104–110. https://doi.org/10.1016/j.brainres.2011.08.043 (2011).

16. Uteshev, V. V., Tupal, S., Mhaskar, Y. & Faingold, C. L. Abnormal serotonin receptor expression in DBA/2 mice associated with susceptibility to sudden death due to respiratory arrest. *Epilepsy Res.* **88**, 183–188. https://doi.org/10.1016/j.eplepsyres.2009.11.004 (2010).

17. Martin, B., Dieuset, G., Pawluski, J. L., Costet, N. & Biraben, A. Audiogenic seizure as a model of sudden death in epilepsy: A comparative study between four inbred mouse strains from early life to adulthood. *Epilepsia* https://doi.org/10.1111/epi.16432 (2020).

18. Faingold, C. L., Tupal, S. & Randall, M. Prevention of seizure-induced sudden death in a chronic SUDEP model by semichronic administration of a selective serotonin reuptake inhibitor. *Epilepsy Behav. E&B* **22**, 186–190. https://doi.org/10.1016/j.yebeh.2011.06.015 (2011).

19. Tupal, S. & Faingold, C. L. Evidence supporting a role of serotonin in modulation of sudden death induced by seizures in DBA/2 mice. *Epilepsia* **47**, 21–26. https://doi.org/10.1111/j.1528-1167.2006.00365.x (2006).

20. Feher, J. J. *Quantitative human physiology: an introduction* (Elsevier/AP, Academic Press is an imprint of Elsevier, 2017).

21. Hodges, M. R. & Richerson, G. B. Medullary serotonin neurons and their roles in central respiratory chemoreception. *Respir. Physiol. Neurobiol.* **173**, 256–263. https://doi.org/10.1016/j.resp.2010.03.006 (2010).

22. Venit, E. L., Shepard, B. D. & Seyfried, T. N. Oxygenation prevents sudden death in seizure-prone mice. *Epilepsia* **45**, 993–996. https://doi.org/10.1111/j.0013-9580.2004.02304.x (2004).

23. Aiba, I. & Noebels, J. L. Spreading depolarization in the brainstem mediates sudden cardiorespiratory arrest in mouse SUDEP models. *Sci. Transl. Med.* **7**, 282ra246. https://doi.org/10.1126/scitranslmed.aaa4050 (2015).

24. Jefferys, J. G. R., Arafat, M. A., Irazoqui, P. P. & Lovick, T. A. Brainstem activity, apnea, and death during seizures induced by intrahippocampal kainic acid in anaesthetized rats. *Epilepsia* **60**, 2346–2358. https://doi.org/10.1111/epi.16374 (2019).

25. Faingold, C. L., Randall, M. & Tupal, S. DBA/1 mice exhibit chronic susceptibility to audiogenic seizures followed by sudden death associated with respiratory arrest. *Epilepsy Behav.* **17**, 436–440. https://doi.org/10.1016/j.yebeh.2010.02.007 (2010).

26. Johnston, S. C., Siedenberg, R., Min, J. K., Jerome, E. H. & Laxer, K. D. Central apnea and acute cardiac ischemia in a sheep model of epileptic sudden death. *Ann. Neurol.* **42**, 588–594. https://doi.org/10.1002/ana.410420409 (1997).

27. Kim, Y. et al. Severe peri-ictal respiratory dysfunction is common in Dravet syndrome. *J. Clin. Investig.* **128**, 1141–1153. https://doi.org/10.1172/JCI94999 (2018).

28. Lertwittayanon, W., Devinsky, O. & Carlen, P. L. Cardiorespiratory depression from brainstem seizure activity in freely moving rats. *Neurobiol. Dis.* **134**, 104628. https://doi.org/10.1016/j.nbd.2019.104628 (2019).

29. Nakase, K. et al. Laryngospasm, central and obstructive apnea during seizures: Defining pathophysiology for sudden death in a rat model. *Epilepsy Res.* **128**, 126–139. https://doi.org/10.1016/j.eplepsyres.2016.08.004 (2016).

30. St-John, W. M., Rudkin, A. H., Homes, G. L. & Leiter, J. C. Changes in respiratory-modulated neural activities, consistent with obstructive and central apnea, during fictive seizures in an in situ anaesthetized rat preparation. *Epilepsy Res.* **70**, 218–228. https://doi.org/10.1016/j.eplepsyres.2006.05.005 (2006).

31. Totola, L. T., Takakura, A. C., Oliveira, J. A., Garcia-Cairasco, N. & Moreira, T. S. Impaired central respiratory chemoreflex in an experimental genetic model of epilepsy. *J. Physiol.* **595**, 983–999. https://doi.org/10.1113/JP272822 (2017).

32. Villiere, S. M. et al. Seizure-associated central apnea in a rat model: Evidence for resetting the respiratory rhythm and activation of the diving reflex. *Neurobiol. Dis.* **101**, 8–15. https://doi.org/10.1016/j.nbd.2017.01.008 (2017).

33. Sanchez, P. E. et al. Optimal neuroprotection by erythropoietin requires elevated expression of its receptor in neurons. *Proc. Natl. Acad. Sci. U.S.A.* **106**, 9848–9853. https://doi.org/10.1073/pnas.0901840106 (2009).

34. Maia, O. A. C. et al. Pilocarpine-induced status epilepticus reduces chemosensory control of breathing. *Brain Res. Bull.* **161**, 98–105. https://doi.org/10.1016/j.brainresbull.2020.05.002 (2020).

35. Kudin, A. P. et al. Seizure-dependent modulation of mitochondrial oxidative phosphorylation in rat hippocampus. *Eur. J. Neurosci.* **15**, 1105–1114. https://doi.org/10.1046/j.1460-9568.2002.01947.x (2002).
36. Kunz, W. S. et al. Mitochondrial complex I deficiency in the epileptic focus of patients with temporal lobe epilepsy. *Ann. Neurol.* **48**, 766–773 (2000).

37. Gastinger, S., Sorel, A., Nicolas, G., Gratas-Delamarche, A. & Prioux, J. A comparison between ventilation and heart rate as indicator of oxygen uptake during different intensities of exercise. *J. Sports Sci. Med.* **9**, 110–118 (2010).

38. Gautier, H. Interactions among metabolic rate, hypoxia, and control of breathing. *J. Appl. Physiol.* **81**, 521–527. https://doi.org/10.1152/jappl.1996.81.2.521 (1996).

39. Refinetti, R. Metabolic heat production, heat loss and the circadian rhythm of body temperature in the rat. *Exp. Physiol.* **88**, 423–429. https://doi.org/10.1113/eph8802521 (2003).

40. Refinetti, R. The circadian rhythm of body temperature. *Front. Biosci. (Landmark Ed)* **15**, 564–594. https://doi.org/10.2741/3634 (2010).

41. Seifert, E. L. & Mortola, J. P. Circadian pattern of ventilation during acute and chronic hypercapnia in conscious adult rats. *Am. J. Physiol. Regul. Integr. Comp. Physiol.* **282**, R244–251. https://doi.org/10.1152/ajpregu.00290.2001 (2002).

42. Stephenson, R., Liao, K. S., Hamrahi, H. & Horner, R. L. Circadian rhythms and sleep have additive effects on respiration in the rat. *J. Physiol.* **536**, 225–235. https://doi.org/10.1111/j.1469-7793.2001.00225.x (2001).

43. Pearce, P. S. et al. Metabolic injury in a variable rat model of post-status epilepticus. *Epilepsia* **57**, 1978–1986. https://doi.org/10.1111/epi.13588 (2016).

44. Luks, A., Glenny, R. W. & Robertson, H. T. *Introduction to cardiopulmonary exercise testing*. (Springer, 2013).

45. Poets, C. F., Meny, R. G., Chobanian, M. R. & Bonofiglo, R. E. Gasping and other cardiorespiratory patterns during sudden infant deaths. *Pediatr. Res.* **45**, 350–354. https://doi.org/10.1203/00006450-199903000-00010 (1999).

46. Totola, L. T. et al. Amygdala rapid kindling impairs breathing in response to chemoreflex activation. *Brain Res.* **1718**, 159–168. https://doi.org/10.1016/j.brainres.2019.05.015 (2019).

47. Hodges, M. R. & Richerson, G. B. Interaction between defects in ventilatory and thermoregulatory control in mice lacking 5-HT neurons. *Respir. Physiol. Neurobiol.* **164**, 350–357. https://doi.org/10.1016/j.resp.2008.08.003 (2008).

48. Morales, A. et al. Unexpected expression of orexin-B in basal conditions and increased levels in the adult rat hippocampus during pilocarpine-induced epileptogenesis. *Brain Res.* **1109**, 164–175. https://doi.org/10.1016/j.brainres.2006.06.075 (2006).

49. Nadam, J. et al. Neuroprotective effects of erythropoietin in the rat hippocampus after pilocarpine-induced status epilepticus. *Neurobiol. Dis.* **25**, 412–426. https://doi.org/10.1016/j.nbd.2006.10.009 (2007).

50. Navarro, F. P. et al. Brain heparanase expression is up-regulated during postnatal development and hypoxia-induced neovascularization in adult rats. *J. Neurochem.* **105**, 34–45. https://doi.org/10.1111/j.1471-4159.2007.05116.x (2008).

51. Manouze, H. et al. Effects of single cage housing on stress, cognitive, and seizure parameters in the rat and mouse pilocarpine models of epilepsy. *eNeuro* **6**, https://doi.org/10.1523/ENEURO.0179-18.2019 (2019).

52. Dube, C., Boyet, S., Marescaux, C. & Nehlig, A. Relationship between neuronal loss and interictal glucose metabolism during the chronic phase of the lithium-pilocarpine model of epilepsy in the immature and adult rat. *Exp. Neurol.* **167**, 227–241. https://doi.org/10.1006/exnr.2000.7561 (2001).

53. Leroy, C. et al. Temporal changes in mRNA expression of the brain nutrient transporters in the lithium-pilocarpine model of epilepsy in the immature and adult rat. *Neurobiol. Dis.* **43**, 588–597. https://doi.org/10.1016/j.nbd.2011.05.007 (2011).

54. Racine, R. J. Modification of seizure activity by electrical stimulation II. Motor seizure. *Electroencephalogr. Clin. Neurophysiol.* **32**, 281–294. https://doi.org/10.1016/0013-4694(72)90177-0 (1972).

55. Belouze, M. et al. Leanness of Lou/C rats does not require higher thermogenic capacity of brown adipose tissue. *Physiol. Behav.* **104**, 893–899. https://doi.org/10.1016/j.physbeh.2011.05.029 (2011).

56. Ogier, M. et al. Brain-derived neurotrophic factor expression and respiratory function improve after ampakine treatment in a mouse model of Rett syndrome. *J. Neurosci.* **27**, 10912–10917. https://doi.org/10.1523/JNEUROSCI.1869-07.2007 (2007).

## Acknowledgements

We thank Frédérique Cohen-Adad (scientific documentation; TIGER team, CRNL) and Denis Ressnikoff (confocal microscopy; Centre d'Imagerie Quantitative Lyon-Est of the Claude Bernard Lyon 1 University) for their assistance. We thank Freya Cracknell for her contribution in editing the manuscript.

## Author contributions

L.B. and S.R. contributed equally to this work. L.B., S.R., and P.R. conceived and designed the experiments. A.M., J.R., and H.K. performed respiratory thermochemistry experiments. H.K. performed plethysmography and telemetry experiments. B.M., G.D., L.B., and H.K. performed AGS in DBA/2 mice. H.K. and B.G. performed histochemistry and RT-qPCR experiments. L.B., S.R., M.O., and H.K. analyzed the data. L.B., S.R., and H.K. wrote the manuscript. All authors reviewed the manuscript.

## 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-14153-6.

**Correspondence** and requests for materials should be addressed to H.K. or L.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_12.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