# Effects of urethane and isoflurane on the sensory evoked response and local blood flow in the early postnatal rat somatosensory cortex

Viktoria Shumkova, Violetta Sitdikova, Ildar Rechapov, Alexey Leukhin,
Marat Minlebaev

► To cite this version:

Viktoria Shumkova, Violetta Sitdikova, Ildar Rechapov, Alexey Leukhin, Marat Minlebaev. Effects of urethane and isoflurane on the sensory evoked response and local blood flow in the early postnatal rat somatosensory cortex. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-88461-8. http-03242697

HAL Id: hal-03242697

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

Submitted on 31 May 2021

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

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

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

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

OPEN

# Effects of urethane and isoflurane on the sensory evoked response and local blood flow in the early postnatal rat somatosensory cortex

Viktoria Shumkova<fnref n="1" />, Violetta Sitdikova<fnref n="1" />, Ildar Rechapov<fnref n="1" />, Alexey Leukhin³ & Marat Minlebaev<fnref n="1" />,²,⁴✉

Functional studies in the central nervous system are often conducted using anesthesia. While the dose-dependent effects of anesthesia on neuronal activity have been extensively characterized in adults, little is known about the effects of anesthesia on cortical activity and cerebral blood flow in the immature central nervous system. Substitution of electrophysiological recordings with the less-invasive technique of optical intrinsic signal imaging (OIS) in vivo allowed simultaneous recordings of sensory-evoked functional response and local blood flow changes in the neonatal rat barrel cortex. Using OIS we characterize the effects of two widely used anesthetics—urethane and isoflurane. We found that both anesthetics suppressed the sensory-evoked optical intrinsic signal in a dose-dependent manner. Dependence of the cortical response suppression matched the exponential decay model. At experimental levels of anesthesia, urethane affected the evoked cortical response less than isoflurane, which is in agreement with the results of electrophysiological recordings demonstrated by other authors. Changes in oxygenation and local blood flow also showed negative correlation with both anesthetics. The high similarity in immature patterns of activity recorded in different regions of the developing cortex suggested similar principles of development regardless of the cortical region. Therefore the indicated results should be taken into account during functional explorations in the entire developing cortex. Our results also point to urethane as the anesthetic of choice in non-survival experimental recordings in the developing brain as it produces less prominent impairment of cortical neuronal activity in neonatal animals.

Anesthesia suppresses neuronal activity in the brain structures and impairs their communication, this is associated with changes in the local cerebral blood flow (CBF)<fnref n="1" />,². The activity of some functional networks is preserved under anesthesia, while others are strongly suppressed³. By acting through various mechanisms, general anesthetics alter cortical responsiveness to sensory inputs, alter functional connectivity, and produce a loss of consciousness⁴⁻⁶. While light anesthesia modulates cortical functioning, deep anesthesia can induce coma and the pathological burst suppression pattern⁷⁻¹¹. The mechanisms of commonly used anesthetics have been extensively studied in the adult brain. Experimental results show that the effect of anesthesia on cortical activity is drug- and dose-dependent¹²⁻¹⁴. Functional connectivity studies using different anesthesia protocols showed they modulate connectivity patterns in unique ways⁶. Compared to the adult brain, the developing brain is more prone to anesthetic complications, because of the less-effective drug-metabolism and narrow therapeutic windows for commonly used anesthetics¹⁵. Many of the anesthetic agents that showed efficiency in adult rodents (ketamine, pentobarbital, and the fentanyl-droperidol combination) are limited in use in neonatal rats because they often lead to inadequate anesthetic depth and have a high mortality rate¹⁶. Therefore, in spite of the variety of

<footnote n="1">¹Laboratory of Neurobiology, Kazan Federal University, Kazan 420000, Russia. ²OpenLab of Neurobiology, Institute of Fundamental Medicine and Biology, Kazan Federal University, Kazan 420000, Russia. ³Laboratory of Neuromorphic computing and neurosimulations, Institute of Information Technology and Intelligent Systems, Kazan Federal University, Kazan 420000, Russia. ⁴INMED, INSERM, Aix-Marseille University, Marseille 13273, France. ✉email: marat.minlebaev@inserm.fr</footnote>
anesthetic agents that can be used to produce general anesthesia or sedation during experiments, the well-known anaesthetics isoflurane and urethane are still commonly used in exploration of the developing brain in rodents.

In contrast to adults, in neonates the activity of the nervous system is characterized by a high level of discontinuity. Transient immature patterns of activity are separated by long-lasting periods of neuronal silence<sup>17-19</sup>, suggesting that the manifestations of anesthesia on immature cortical activity may be different. Electrophysiological recordings in the visual system showed strong suppression of the light evoked cortical response even at low doses of isoflurane (0.3-1%) in P10-11 rats, however, the response was moderately facilitated in P13 animals<sup>20</sup>, suggesting developmental changes in the effects of anesthesia in the immature nervous system. The suppressing effect of isoflurane on cortical activity patterns was also shown in the barrel cortex during the first two weeks of postnatal development in rats<sup>21</sup>. At a low dose of isoflurane spontaneous activity was weakly modified, while a 1.9-2.3% concentration of isoflurane strongly suppressed immature patterns of cortical activity, supporting the graded effect of anesthesia on cortical neuronal activity in neonatal rats, as in adult rats. Recordings of spontaneous cortical activity in neonatal mice showed that the administration of isoflurane or urethane reduced the active periods of neonatal activity without significant alteration of the spectral features of the cortical activity patterns and neuronal firing rate<sup>22</sup>. Therefore, while the effects of general anesthesia on the functioning of the central nervous system have been thoroughly characterized in adults, the developmental aspects of their action remain less well understood.

Cerebral blood flow is also an object of dose- and drug dependent modulation by sedatives, analgesics, and anesthetics<sup>23-26</sup>. A dose-dependent CBF reduction in the mouse brain was shown for isoflurane<sup>27</sup>. It was also demonstrated that isoflurane anesthesia evoked a higher baseline CBF during stimulation of the cerebral cortex<sup>8</sup>. Urethane shows minimal effects on the cardiovascular and respiratory systems<sup>28</sup>. While in adults anesthesia acts drug- and dose-dependently, little is known about the effects of different anesthetic protocols on the cortical activity and associated CBF changes in the developing central nervous system.

Using optical intrinsic signal (OIS) recordings to characterize the changes in cortical activity and cerebral blood flow with different types of anesthesia, we attempted to answer this question. Single light wavelength OIS is a simple and powerful less invasive technique to study brain activity<sup>29</sup>. In clinical studies, completely non-invasive implementation of OIS recordings (functional Near-Infrared Spectroscopy, fNIRS) is a routine technique for brain studies<sup>30,31</sup>. Those techniques mainly reflect the hemodynamic cortical response and metabolism, however modulation of light scattering associated with the neuronal activity also contributes to the generation of the intrinsic signal<sup>32</sup>. In studies in adults, conventional OIS recordings and fNIRS are considered to be effective substitutes for the blood oxygenation level-dependent (BOLD) signal recorded by functional magnetic resonance imaging (fMRI)<sup>33,34</sup>.

Decomposition of multispectral OIS recordings allows the extraction of both the hemodynamic and light scattering responses<sup>35-38</sup>. The contribution of different chromophores (oxy- and deoxyhemoglobins), and tissue light scattering (LS) suggested various mechanisms underlying the signal generation recorded by OIS and BOLD fMRI. This may explain the efficiency of the OIS approach in neonatal rats, compared to BOLD fMRI<sup>39-41</sup>. OIS imaging carries rich information about neuronal activity, which manifests from several processes that dictate light propagation through biological tissue (absorption and scattering). While the light absorption signal reflects the changes in blood oxygenation related to neuronal activity, LS is linked to changes in the optical properties of neural tissue as a secondary consequence of neural excitation. In the mature brain, the LS response is overlapped in time and in amplitude by the hemodynamic response. Because of the immaturity of the nervous and vascular systems, particularly the immaturity of neurovascular coupling, both neural tissue and hemodynamic responses are much slower and separated in time<sup>37,42</sup>. The early phase of OIS, during the first seconds of the cortical response, is characterized by light scattering changes<sup>40</sup>. The high correlation of the optically recorded signal with the power of the evoked neuronal oscillatory response in the immature barrel cortex strongly supports the idea that in the developing brain OIS imaging may serve as an approach for direct measurement of cortical functional response. The second phase of the OIS response, characterized by changes in oxy- and deoxyhemoglobins, indicates changes in CBF<sup>43</sup>. Therefore, in the developing nervous system OIS can be used for characterization of the changes in evoked cortical activity and cerebral blood flow.

Therefore, from the experimental point of view, it is critically important to know own effect of drugs on cortical functioning and local hemodynamics for correct interpretation of the results. In the present study, we demonstrate the effects of isoflurane and urethane on sensory evoked functional response and the associated local blood flow changes. We report that in the immature somatosensory cortex both isoflurane and urethane induce dose- and drug dependent alterations in sensory evoked functional response and the associated local blood flow changes.

## Results

**OIS dependence on anesthesia concentration in neonatal rat pups.** To characterize the amplitude-temporal parameters of the OIS evoked in the barrel cortex of neonatal rat pups, conventional OIS recording with red light illumination was used<sup>29,44</sup> (Fig. 1A). A train of rhythmic vibrissa deflections reliably evoked a drop in reflected light that spanned over tens of seconds (see Table 1, Fig. 1B,C), in agreement with the previous results<sup>39,40</sup>. To estimate the cortical area involved in the evoked cortical response we developed an approach based on kernel density estimation and detected the contours of the OIS (see "Materials and method" section for details).

In spite of the variation of the detected contours between the frames, the maximal contour density coincided with the visually detected OIS (Fig. 1B). OIS onset was slightly delayed compared to the start of stimulation (Table 1, Fig. 1C,D). Analysis of the onset times between animals of different ages does not show significant developmental changes (Table 1). The cortical area involved in OIS quickly reached a peak size dependent on
<table><caption>**Table 1.** The amplitude-temporal parameters of the OIS.</caption><thead><tr><th rowspan="2">Parameter</th><th colspan="2">Median</th><th colspan="2">25%</th><th colspan="2">75%</th><th rowspan="2">Spearmans coefficient, ρ</th><th rowspan="2">Significance, α = 0.05</th></tr><tr><th>P6</th><th>P10</th><th>P6</th><th>P10</th><th>P6</th><th>P10</th></tr></thead><tbody><tr><td>Age</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>OIS duration (s)</td><td>60.48</td><td>40.49</td><td>39.76</td><td>34.76</td><td>64.28</td><td>42.2</td><td>-0.26</td><td>Non significant</td></tr><tr><td>OIS onset (s)</td><td>1.08</td><td>0.6</td><td>0.52</td><td>0.4</td><td>1.16</td><td>1.04</td><td>-0.08</td><td>Non significant</td></tr><tr><td>OIS amplitude (%)</td><td>0.16</td><td>0.21</td><td>0.14</td><td>0.18</td><td>0.17</td><td>0.29</td><td>0.36</td><td>Significant</td></tr><tr><td>OIS rise time (s)</td><td>2.04</td><td>2.68</td><td>1.68</td><td>2.12</td><td>2.76</td><td>3.04</td><td>0.05</td><td>Non significant</td></tr><tr><td>OIS decay time (s)</td><td>10.52</td><td>11</td><td>9.32</td><td>8.56</td><td>15.44</td><td>13.4</td><td>-0.12</td><td>Non significant</td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Evoked optical intrinsic signal in the somatosensory cortex of the neonatal rat. (A) Diagram of the experimental setup for recording evoked OIS in the barrel cortex of the neonatal rat, (B) color-coded density map of the detected OIS contours overlaid on the recorded OIS, seen as a dark spot, (C) example of the temporal profile of the OIS recorded using red light, stimulation period is shown by the grey rectangle, mean OIS dynamic is shown by the red line, the shaded area corresponds to the confidence interval, (D) developmental changes in OIS onset (0 is the start of stimulation), (E) example of OIS area changes during recording is shown by the blue line, the red line is a result of convolution with a Kaiser window of 12 points size, (F) developmental changes in the evoked OIS area, amplitude, rise and decay times. The red dashed line is a linear fit. Grey dots are results from individual experiments. Box plots show median (bold black line), the bottom and top edges of the box indicate the 25th and 75th percentiles. The whiskers extend to the most extreme data points not considered outliers, and the outliers are plotted individually using the '+' symbol.</figcaption></figure>

the animal's age. In P6 rat pups the mean OIS area was less than $1 \text{ mm}^2 (0.68 \pm 0.17 \text{ mm}^2)$, by P10 evoked OIS covered larger cortical territories ($1.65 \pm 0.25 \text{ mm}^2$, $\rho = 0.59$, $p < 0.05$, n = 23 P6-10 rat pups), this is in an agreement with the developmental increase of the barrel cortex (Fig. 1E,F).

Though there is a positive correlation of OIS area and animal maturation, we have not found significant changes in the OIS amplitude during animal development (Table 1, Fig. 1F). No significant differences were found for the rise time of the OIS nor for the decay time (Table 1, Fig. 1F). We grouped recordings in spite of the different ages of the rat pups, due to the lack of the developmental changes in OIS parameters between the animals.

To examine the dependence of evoked OIS on anesthesia, experiments were done with urethane and isoflurane$^{17,18,39,45}$. We gradually increased the concentrations of the anesthetics and monitored the evoked OIS. Firstly, using Red illumination we characterized the cortical area of the OIS. The area of evoked OIS was weakly modulated by different concentrations of isoflurane or urethane ($\rho = 0.21$ and $\rho = -0.2$ for isoflurane and urethane, respectively, n = 14 P7-10 rat pups, Fig. 2A,B). Elucidation of the components underlying the OIS requires multi spectral recordings of the signal, therefore experiments were carried out using three types of illumination (Red, Green, and IR). A progressive reduction in OIS amplitude was observed with increasing concentrations of isoflurane or urethane using all three illumination wavelengths (Fig. 2C,D). The dependence of OIS amplitude on the anesthesia concentration matched best with an exponential model for both anesthetics in all types
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Concentration dependence of OIS on isoflurane and urethane. (A) Examples of OIS recorded using red light in control and under 0.25%, 0.75% isoflurane (left) and 1 g/kg, 1.5 g/kg urethane (right) anesthesia. OIS contours in control and under different concentrations of the anesthetics are colour-coded. Group data of normalized red OIS area for changes in concentration of isoflurane and urethane are shown in (B). OISs recorded using different light wavelengths (Green, Red and IR) in control and under 1.5% isoflurane (C) and 1.5 g/kg urethane (D). The stimulation period is shown by the grey rectangle. Light wavelength is colour-coded. OISs recorded in control conditions are shown by dashed line, while OISs traces under anesthesia are shown by semitransparent lines of the corresponding colour. (E) Dependencies of the normalized amplitudes of the OISs recorded using Green, Red and IR light with different concentrations of isoflurane (left) and urethane (right) are shown. Triangle marks the concentrations at which the examples of multi spectral OISs are shown on (C, D).</figcaption></figure>

of illumination (r² for OISs recorded using Green, Red and IR lights were 0.99, 0.99, 0.99 and 0.98, 0.98, 0.98, respectively, n = 6 and n = 8 rat pups of P7-10, Fig. 2E).
<table><caption>**Table 2.** Developmental changes of OIS baseline dependence on anesthesia (Spearman's rank correlation coefficient, $\rho$).</caption><thead><tr><th></th><th>Green illumination</th><th>Red illumination</th><th>IR illumination</th></tr></thead><tbody><tr><td>Urethane</td><td>2 ± 2 × 10<sup>-5</sup></td><td>0.5 ± 1 × 10<sup>-5</sup></td><td>1.5 ± 1 × 10<sup>-5</sup></td></tr><tr><td>Isoflurane</td><td>-0.5 ± 1 × 10<sup>-6</sup></td><td>-1.3 ± 0.4 × 10<sup>-6</sup></td><td>-2 ± 11 × 10<sup>-5</sup></td></tr></tbody></table>

Since the OIS is calculated using normalization for its pre-stimulation baseline, changes in recorded light intensity (reflection for Red and Green wavelengths, and transmission for IR light) during and before stimulation may affect the calculation of OIS. In order to test the impact of the baseline changes on calculating the dependence of OIS on the anaesthesia concentration, we also estimated the amplitude of the recorded light intensity as a function of the concentration of the anesthetics. The results showed that recorded light intensity during baseline weakly depended on the concentration and type of anesthetic. We did not see developmental changes in OIS baseline dependence on anesthesia concentration (Table 2). Therefore, modulation of the OIS amplitude is exclusively associated with evoked response changes at different concentrations of urethane or isoflurane.

To characterize the developmental changes of the OIS dependencies on the anesthetics the exponential concentration constants (drop in OIS amplitude for e time) were calculated. Our results showed that concentration dependence of OIS varies between the wavelength of light used (the average exponential concentration constants for Green, Red and IR OISs were 1.1 ± 0.1%, 1.3 ± 0.2%, 1.3 ± 0.2% for isoflurane and 1.1 ± 0.1 g/kg, 1.5 ± 0.1 g/kg and 1.6 ± 0.1 g/kg for urethane). As the principal chromophores are characterized by different absorption spectra, the difference in concentration constants for Green, Red and IR light wavelengths suggests that the anesthetics affect different OIS mechanisms.

**Dose-dependent changes of the evoked cortical response.** There are several mechanisms underlying OISs, and their contribution to the generation of OIS depends on the maturation of the nervous and vascular systems. While in the mature brain the OIS is predominantly determined by the hemodynamic changes following neuronal activity, in the immature nervous system the early part of the OIS is largely defined by the tissue component$^{40}$.

To characterize the effects of urethane and isoflurane on the evoked cortical response in the immature nervous system we estimated the changes in the tissue component (light scattering, LS), local blood flow (total hemoglobin, HbT), and tissue oxygenation level (oxygenated hemoglobin, HbO) in the barrel cortex following sensory stimulation. Analysis of the multi spectral OIS recordings showed that an increase in the concentration of the anesthetics resulted in a reduction of the amplitude of the LS without changes in the cortical area involved in the cortical functional response (Fig. 3A,B). While in control conditions the mean area of the functional response was slightly bigger than $1 \text{ mm}^2 (1.37 \pm 0.12 \text{ mm}^2$, n = 14 P7-10 rats), under different concentrations of isoflurane or urethane, the area of the LS response changed non significantly by only a few percent ($\rho = -0.21$ and $\rho = -0.36, -1.90 \pm 0.766$% and $3.79 \pm 19.37\%$ from control area for 1.5% isoflurane and for 1.25 g/kg of urethane, respectively, n = 14 P7-10 rats). The exponential fit model highly matched LS amplitude dependence on anesthetic concentrations ($r^2 0.99$ and $0.99$ for urethane and isoflurane, respectively) showing non linear correlation of functional cortical response on depth of anesthesia ($e$ concentration constant $1.4 \pm 0.1$ g/kg for urethane and $1.2 \pm 0.1$% for isoflurane, n = 14 P7-10 rats, Fig. 3C,D).

Though urethane and isoflurane belong to different classes of anesthetic, we compared their efficiency on the evoked functional activity of the barrel cortex Fig. 4. Our results showed that a slight increase in isoflurane concentration in the “subexperimental” concentration range (lower than minimal alveolar concentration, MAC) strongly affected the LS amplitude (an increase of the isoflurane concentration by about $0.23 \pm 0.01\%$ reduced LS by 10%), while in the MAC range ($2.21-2.47\%$), the isoflurane modulation of LS was less powerful and required an almost threefold increase in concentration to change LS by 10% ($0.83 \pm 0.05\%/10\%$ LS, Fig. 4), demonstrating strong changes in isoflurane efficiency at different concentrations. LS dependence on the urethane concentration showed a less pronounced exponential fit. The LS changes in the experimental concentration range ($1-1.5$ g/kg) and “subexperimental” ranges were smaller ($0.2 \pm 0.01$ g/kg per 10% of LS and $0.31 \pm 0.01$ g/kg per 10% of LS for “subexperimental” and experimental ranges of urethane concentrations, respectively). We have also found that commonly accepted concentrations for isoflurane and urethane affected LS differently. While urethane at a concentration of $0.9-1.5$ g/kg reduced LS by up to three times (to $50-30\%$ from control), an experimental dose of isoflurane reduced the functional cortical response by more than seven times (LS dropped to $13-10\%$ from control) and in some experiments the tissue component was almost completely blocked.

We also analyzed the dose- and anesthetic-dependence of CBF on the different concentrations of urethane and isoflurane. The hemodynamic response following evoked cortical activity also showed negative dose-dependence on both anesthetics (Fig. 5A). Similar to the LS, increase of the anesthesia concentration suppressed local increase of the oxy, and total hemoglobins following whisker stimulation Fig. 5A. Those changes matched the exponential decay model (concentration constant $e$ for isoflurane was $1.2 \pm 0.1\%$, $1.1 \pm 0.1\%$ and for urethane $1.1 \pm 0.1$ g/kg, $1.1 \pm 0.1$ g/kg, respectively, Fig. 5B). While a direct effect of anesthesia on cerebral blood flow was shown in adults, we have seen that the hemodynamic response almost linearly correlated with the LS changes (Pearson coefficients for HbO changes under isoflurane $-0.99$ and urethane anesthesia $-0.99$; HbT changes under isoflurane $-0.99$ and urethane $-0.99$, n = 14 $p < 0.01$, Fig. 5C).
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Dose-dependence of functional response on isoflurane and urethane. LS responses in control and under 0.25%, 0.75% of isoflurane (A) and 1 g/kg, 1.5 g/kg urethane (B) anesthesia. LS responses in control and under different anesthetic concentrations are colour-coded. Group data of normalized area changes are shown by box plot graphs under the examples of LS responses. Individual dependencies of normalized LS changes on different isoflurane (C) and urethane (D) concentrations are shown in the upper row. Dots correspond to the LS amplitude at different concentrations of anesthetic agent. Exponential model fits are shown by lines of the corresponding colour. Examples of LS changes in time in control and under 1.5% of isoflurane and 1.5 g/kg urethane are shown as insets. Group data are shown in the bottom row. Exponential fits from each experiment are shown by grey lines, while mean fits are shown by blue and red lines for isoflurane and urethane, respectively.</figcaption></figure>

Discussion

The effects of anesthesia on cortical activity and cerebral blood flow in the developing nervous system need to be investigated because of their importance to functional studies in vivo and their interpretation. From this point of view, our results enriched this line of investigation and contributed novel insights into the actions of two widely known anesthetics, urethane and isoflurane. Using the less-invasive technique of intrinsic optical signal imaging in the barrel cortex of the neonatal rat model we show dose-dependent effects of urethane and isoflurane on the cortical evoked response, changes in local blood flow, and oxygenation.

Isoflurane and urethane are anesthetics widely used in neuroimaging and electrophysiological investigations in rodents in vivo, serving to eliminate physiological stress and motion effects. While several studies assessed the effect of anesthesia on immature cortical activity, the fidelity between the different anesthetics and the functional cortical response and cerebral blood flow in neonates have not previously been examined, and therefore warrants thorough investigation and discussion<sup>21,22,46</sup>.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** LS modulation by isoflurane and urethane. Functional effects of isoflurane and urethane on LS are shown. The type of anesthetic is colour-coded, the shaded areas correspond to the confidence intervals. Horizontal blue and red lines display the experimental concentration ranges for isoflurane and urethane, respectively. Note, that the experimental concentration of isoflurane suppresses LS by up to seven times, while for the experimental concentration of urethane LS decreased by only 2–3 times.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Dependence of local cerebral blood flow on anesthesia concentration and LS changes. (A) Oxygenated hemoglobin (HbO, red colour) and total hemoglobin (HbT, blue colour) changes in control and under 1.25% isoflurane (left) and 1.5 g/kg urethane (right). The grey box corresponds to the stimulation period. Control is shown by dashed line of corresponding colour, while semitransparent lines correspond to the anesthesia condition. (B) Group data of normalized concentration dependence of HbO, and HbT on isoflurane (left) and urethane (right). Individual examples are shown as insets. Triangles marks the concentrations at which the examples of HbO and HbT are shown on (A). (C) HbO, HbT dependencies on the LS modified by isoflurane (left pair) or urethane (right pair).</figcaption></figure>

The early postnatal period of somatosensory cortex development is characterized by the expression of unique patterns of immature oscillatory activity. These can be spontaneous, or evoked by stimulation of the sensory receptive fields and observed in the corresponding cortical representations of the body in the somatosensory cortex<sup>17,18,20,47</sup>. These patterns are also seen in the developing visual cortex<sup>20,48</sup>, suggesting that cortical development undergoes similar principles regardless of the cortical region. Based on the high correlation of the OIS
with electrophysiological recordings in the developing barrel cortex⁴⁰, we suggest that the observed changes in OIS may reflect the common principles of anesthesia action both on local blood flow and the cortical response, regardless of the cortical region.

We show that both of the anesthetic agents produced a consistent dose-dependent alteration of the cortical evoked response to vibrissa stimulation. This was evidenced by dose-dependent decreases of light scattering of the barrel cortex when receiving multiple sensory stimulation. The anesthetic dose-dependent reduction in evoked and spontaneous activity is shown in previous electrophysiological findings in the adult nervous system⁸,¹⁴,⁴⁹,⁵⁰.

The cellular and molecular pharmacology of anesthetics has been reviewed extensively, yet the mechanisms underlying the actions of anesthesia are not completely understood. Dose-dependent depression of glutamatergic transmission, and potentiation of GABAa transmission has been shown for both isoflurane and urethane²⁸,⁵¹⁻⁵³. Though isoflurane targeting of both excitatory and inhibitory systems is debatable. Isoflurane (as well as other volatile anesthetics such as halothane and sevoflurane) was shown to potentiate GABAa and glycine receptors⁵⁶,⁵⁷, while AMPA receptors were less affected⁵⁸. However, compared to such intravenous general anesthetics as propofol and etomidate, whose primary target is GABAa receptors, the amplifying effect of isoflurane on GABAa function was less prominent⁵⁹,⁶⁰. This supports the idea that anesthesia effects produced by isoflurane involves extra mechanisms. Isoflurane was also reported to potently inhibit neuronal nACh receptors⁶¹,⁶² and activated two-pore potassium channels, such as TREK-1 (K2P2.1) and TASK (K2P3.1)⁶³ that contribute to its anesthetic action. In contrast to isoflurane, urethane showed only a limited effect on nACh receptors²⁸. But urethane was shown to depress the intrinsic excitability of layer 4 neurons and increased the shunting potassium background leak conductance⁶⁴. Urethane-dependent current appeared to be through the opening of Ba2+-sensitive potassium channels, whose properties are consistent with two-pore K+ channels such as TWIK (K2P1.1)⁶⁵. Isoflurane-sensitive TREK-1 and TASK potassium channels did not appear to contribute to the suppressive action of urethane on neuronal activity. Therefore isoflurane and urethane are distinct in their spectrum of action that may explain their differences in anesthesia. In spite of the immaturity of the nervous system and the progressing development of the cortical circuitry, these mechanisms of anesthesia may also play role in the dose dependent reduction of evoked cortical network activity. The exponential model fits observed for both anesthetics display the process of tissue saturation by the anesthesia, which is not agent-specific. A progressive decrease in the apparent saturation level of the LS reflects a progressive suppression of cortical neurons responding simultaneously to sensory stimuli until the neuronal involvement is minimal.

The anesthetic action on thalamic cells (shown during deep anesthesia⁶⁶) which have been hyperpolarized because of the anesthetic action on GABAa or activation of potassium channels⁵³,⁶⁷, may also provide the conditions for the suppression of thalamic input into the barrel cortex. The emergence of the arousal system that occurs around this age also has a direct effect via reduced acetylcholine level²⁰,⁶⁸. Compared to urethane, isoflurane much more efficiently suppressed evoked responses even at low concentrations. While the light scattering component of the OIS indicates neuronal activity, the difference between isoflurane and urethane is likely to result from the unique modulation of cortical response by different anesthetics. Electrophysiological studies showed that immature patterns of evoked activity were weakly modified by urethane in the concentration range of 0.5–1.5 g/kg³⁹. But strong suppression of the evoked cortical response by isoflurane was also demonstrated in the developing barrel and visual cortices²⁰,²¹. Interestingly isoflurane suppresses the oscillatory part of the evoked response while the sensory evoked potential was left almost intact. These unique mechanisms of isoflurane action on the sensory evoked potential and the oscillatory part of the response require further investigation. However, inhibition of the late phase of the evoked response may be the result of the combination of suppression of cortical activity and recurrent inhibition from the reticular nucleus to the thalamic principal cells⁵³. This had an added effect on the generation of thalamus driven immature oscillatory activity during this developmental period³⁹. The less prominent effect of urethane suggests better preservation of the cortical functioning that makes urethane the anesthetic of choice in studies in the developing nervous system. Comparative fMRI studies with different anesthetics also showed that functional connectivity pattern under urethane anesthesia was close to the non-anesthetized animal (only 6% of all connections were significantly suppressed)⁶. Therefore urethane anesthesia exhibits the fewest differences in functional connectivity in the cortex supporting our conclusion of preferential use of urethane.

The previously shown discrepancy of the dose-dependence effects on the evoked cortical activity in mice and rat pups may be due to the maturity of the central nervous system. The advanced maturation of mice compared to the rats explains the difference in the action of anesthesia in these studies⁶⁹. The developmental decline in the efficiency of anesthesia is also supported by the age-dependent decrease in the effect of anaesthesia on the spontaneous and evoked cortical activity in the barrel and visual cortices²⁰,²¹.

In contrast to urethane, which has restricted effect on the cardiovascular system, isoflurane is also known as a vasodilator⁷⁰,⁷¹, it is likely that the cerebrovascular reactivity to isoflurane is dose-dependent, even in neonates. However, independently of the type of anesthetic agent used, a linear dependence of the HbO and HbT changes on the LS decrease was observed. Therefore, there were not consistent difference between the effects of isoflurane or urethane on the intensity of the neurally coupled hemodynamic response. There may be developmental changes in the isoflurane provoked vasodilation that are not yet matured at this age in the rat pup. This concept is consistent with studies demonstrating developmental changes in vasomotor tone of the pulmonary artery in the neonatal rat⁷². The suppression of the hemodynamic changes under anesthesia shown are also in agreement with other authors. Early optical imaging experiments revealed that in the awake animal, hemodynamic signals are substantially larger in response to whisker stimulation⁷³ or limb stimulation⁷⁴,⁷⁵ than in anesthetized animals.

Our study revealed that the sensory evoked cortical functional response and local CBF were uniquely modulated by isoflurane and urethane. Our results highlight the necessity of careful monitoring of the anesthesia protocols for reliable comparison between experimental results. We also want to point out that urethane allows better control over the depth of anesthesia, as opposed to isoflurane, and affects the cortical functions in the developing brain less.
## Methods

**Surgery.** All animal-use protocols followed the guidelines of Kazan Federal University on the use of laboratory animals (ethical approval by the Institutional Animal Care and Use Committee of Kazan State Medical University N9-2013). Wistar rats of both sexes from postnatal days [P] 6–10 were used (P0 was the day of birth). Surgery was performed under isoflurane anesthesia (5% for induction and 1–2% during surgery). The rat skull was cleaned of skin and periosteum and covered with an acrylic (Meliodent RR) except for a 5 × 5 mm² window above the barrel cortex. A metal ring with an inner diameter of 5 mm was attached to the cement on the head of the animal to make a chamber. Subsequently, the metal ring was fixed into the stereotaxic apparatus. To increase the transparency of the bone in the bottom of the chamber the opened skull part was polished and the chamber was filled with saline and covered with a coverslip. The rats were warmed, surrounded by a cotton nest, heated via a thermal pad (35–37°C) and left for an hour to recover from anesthesia.

**Optical intrinsic signal recordings and analysis.** OIS was recorded in head-fixed rat pups using a video acquisition system. A CCD camera (QICAM Fast 1394) was positioned above the barrel cortex using stereotaxic coordinates⁷⁶. The camera was focused 400–800 µm below the skull, which approximately corresponds to the depth of the granular layer of the barrel cortex as indicated by electrophysiological activity⁴⁷. The cortex was illuminated by light emitting diodes (LED) positioned around the microscope objective, managed and synchronized with the camera by an Arduino Uno microcontroller. Video frames were recorded at a resolution of 130 × 174 (1 px = 35 µm) and a frame rate of 6 Hz. Single wavelength illumination with a 625 nm LED was used to characterize developmental changes in OIS, while three types of diodes (528 nm, 625 nm and 850 nm herein referred to as Green, Red and IR, respectively) were used to record OIS for analysis of the functional response and cerebral blood flow changes. The actual spectra of the diodes were measured using a Thorlabs CCS175 spectrometer. The effects of anesthesia on the evoked cortical response and CBF were estimated using alternate illumination of the barrel cortex with the three types of diodes (6 Hz/LED). While the Green and Red diodes were positioned above the skull around the head to produce spatially uniform illumination, the IR diode was placed below the head of the animal facing towards the camera. Different wavelengths were used for post-hoc estimation of the volumetric and oximetric changes in hemodynamic response and tissue light scattering changes as previously described³⁷,⁴⁰. OIS was evoked by vibrissa deflections. A train of brief pulses (10 ms pulse duration, 10 s train duration, 90 s inter-train interval, repeated 12–20 times) was applied at a frequency of 2 Hz via a piezo deflector (Noliac). A piece of soft foam on the tip of the piezo bender touched most of the vibrissa on the animal’s contralateral snout, providing simultaneous deflection of multiple whiskers. Data pre-processing included per-frame spatial filtering with a 2D Gaussian filter (sigma = 2 px) and subsequent illumination artifact correction. The illumination correction consisted of, first, calculating the reference time profile averaged over the reference region, followed by its subtraction from each pixel intensity⁷⁷. The pre-processed frames were then averaged across all trials (n = 10–20). The experimental protocol was one-conditioned, so the OIS map was calculated using the first-frame subtraction approach as described previously²⁹,⁴². Preliminary detection of the OIS signal was carried out by an operator based on visual detection of a group of pixels with conjointly changed brightness during the stimulation. The OIS time profile was calculated for the manually chosen OIS center. To facilitate visualization of OIS in reflection mode, the signals were inverted. In the experiments with OIS decomposition for hemodynamic and functional responses, the data were recorded using alternate lighting with three different light wavelengths. To calculate the averaged extinction coefficients of oxy- and deoxyhemoglobin (HbO and HbR), spectrophotometric data on rat blood absorption was used⁷⁸. ΔHbO, ΔHbR, and the tissue component were calculated using the modified Beer-Lambert law as described previously³⁸,⁴⁰.

**Cortical response contour detection.** To calculate OIS contours the fragmentation analysis approach was used. Firstly, a threshold level of 95–98% was estimated based on the intensity distribution of all pixels of the entire image, followed by detection of the pixels with intensity exceeding the threshold. Second, kernel density estimation (KDE) was used to estimate the probability density function of these pixels. A group of united pixels with density exceeding half of the probability density function were considered as OIS and the area was calculated. Two additional conditions were used to increase the accuracy of detection: (1) the manually chosen center of OIS should be inside the contour of the automatically detected OIS, (2) the minimal fullness value of the detected contours was also set at 50–90%. Fulfilment of these conditions served as confirmation of the automatically detected OIS. The code for OIS contour detection is provided at the following address (http://github.com/research-team/memristive-spinal-cord/blob/master/analysis/signal_recognition.py). While the amplitude of the recorded OIS was detectable at up to 3 g/kg of urethane or 4% isoflurane, significant detection of OIS contours was possible only at lower concentrations, 1.5 g/kg for urethane or 1.25% for isoflurane, therefore the concentration range of anesthetics for spatial analysis of OIS and its components is limited to 1.5 g/kg of urethane and 1.25% for isoflurane in the text.

**Anesthesia.** Two anesthetics (urethane and isoflurane) were used in this study. These anesthetics have different mechanisms of action and are different in delivery method, both are widely used in neuroscience. During experiments the urethane (Sigma-Aldrich, USA) was delivered via intraperitoneal injection of 15% solution. The volume of injection was calculated based on the weight of rat pups, measured before the experiment (median weight 23.8 g, Supplementary Fig. 1, Supplementary Table 1). The volume of urethane solution was calculated for every animal, because of the individual and developmental variances (to increase urethane concentration for 0.5 g/kg, the median single injection volume was 0.08 ml (25% 0.07, 75% 0.083, Supplementary Fig. 1). Weight comparison between the groups of animals used for characterization of the effects of urethane and isoflurane anesthesia did not show a significant difference between the groups.
Inhalation anesthesia by isoflurane (Baxter, USA) was performed using a rodent inhalant anesthesia apparatus (E-Z Anesthesia, USA). The flow rate of isoflurane was set as 0.6-0.8 l/min. The OIS acquisitions were started at a 10 min delay from the change in the anesthesia concentration. Previous studies on cortical activity in the developing mouse cortex using urethane anesthesia, showed a prominent decrease of cortical oscillatory activity observed 5-15 min after urethane injection, followed by partial recovery towards the baseline²². A rapid decrease of the period of activity following urethane injection was observed in different brain areas regardless of the concentration of anesthesia. The temporal sequence of the urethane quick onset and long-lasting effect of urethane was suggested to be linked to the pharmacokinetics of urethane. Isoflurane also has a short duration of action. The inductive concentration of isoflurane (5%) resulted in a significant decline of heart and respiratory rates within 10 min⁷⁹.

**Data analysis.** Statistical analysis was performed using the MATLAB Statistics toolbox. To avoid bias that may be produced by high variation in the data, comparisons between the experiments and group statistics were done using normalization of each experiment to its control value. The control values were calculated using the OIS recorded prior to the first application of the anesthetic agent. To detect the amplitude, the OIS signal was firstly convolved with a 50 point Kaiser window with shape factor 2, followed by the detection of the maximal peak value. The type of the model fit of the data was estimated based on the coefficient of determination (r²). The significance of developmental changes was estimated using Spearman's rank correlation coefficient (ρ) with the significance level set at *p* = 0.05. The linearity of dependence was tested by estimation of the Pearson correlation coefficient with the significance level set at *p* = 0.05. Group data are expressed by box plots with the central mark indicating the median, and the bottom and top edges of the box indicating the 25th and 75th percentiles, respectively. The whiskers extend to the most extreme data points not considered outliers, and the outliers are plotted individually using the '+' symbol. The confidence interval (CI) was shown as a shaded area and calculated by using 2.5 Jackknife standard deviation that corresponds to a significance threshold of *p* = 0.05.

**Animal research: reporting of in vivo experiments (ARRIVE).** We state that the study was carried out in compliance with the ARRIVE guidelines. Although there is little information in the literature regarding specific signs of stress and pain in rodent pups, during our experiments we visually controlled the presence of continuous mass movements. Presence of vocalization was also used as a sign of stress of the neonatal rat pup. If these signs were present, the experimenter rechecked the positioning of the animal to eliminate the source of discomfort. If the animal continued active movements and vocalization, the experiment was stopped. However, mimicking the natural conditions for the neonatal rat maintained the animal in a stressless state in most of the experiments.

**Data availability**

Original and processed data, signal processing and analysis routines are available on request from the corresponding author.

Received: 9 October 2020; Accepted: 31 March 2021

Published online: 05 May 2021

**References**

1. Hudetz, A. G. General anesthesia and human brain connectivity. *Brain Connectivity* 2, 291–302. https://doi.org/10.1089/brain.2012.0107(MaryAnnLiebertInc (2012).

2. Gao, Y.-R.R. et al. Time to wake up: Studying neurovascular coupling and brain-wide circuit function in the un-anesthetized animal. *Neuroimage* 153, 382–398. https://doi.org/10.1101/077024 (2016).

3. Nallasamy, N. & Tsao, D. Y. Functional connectivity in the brain: Effects of anesthesia. *Neuroscientist* 17, 94–106. https://doi.org/10.1177/1073858410374126 (2011).

4. Grandjean, J., Schroeter, A., Batata, I. & Rudin, M. Optimization of anesthesia protocol for resting-state fMRI in mice based on differential effects of anesthetics on functional connectivity patterns. *Neuroimage* 102, 838–847. https://doi.org/10.1016/j.neuroimage.2014.08.043 (2014).

5. Jonckers, E. et al. Different anesthesia regimes modulate the functional connectivity outcome in mice. *Magn. Reson. Med.* 72, 1103–1112. https://doi.org/10.1002/mrm.24990 (2014).

6. Paasonen, J., Stenroos, P., Salo, R. A., Kiviniemi, V. & Gröhn, O. Functional connectivity under six anesthesia protocols and the awake condition in rat brain. *Neuroimage* 172, 9–20. https://doi.org/10.1016/j.neuroimage.2018.01.014 (2018).

7. Hamilton, C., Ma, Y. & Zhang, N. Global reduction of information exchange during anesthetic-induced unconsciousness. *Brain Struct. Funct.* 222, 3205–3216. https://doi.org/10.1007/s00429-017-1396-0 (2017).

8. Masamoto, K., Fukuda, M., Vazquez, A. & Kim, S. G. Dose-dependent effect of isoflurane on neurovascular coupling in rat cerebral cortex. *Eur. J. Neurosci.* 30, 242–250. https://doi.org/10.1111/j.1460-9568.2009.06812.x (2009).

9. Devonshire, I. M., Grandy, T. H., Dommett, E. J. & Greenfield, S. A. Effects of urethane anaesthesia on sensory processing in the rat barrel cortex revealed by combined optical imaging and electrophysiology. *Eur. J. Neurosci.* 32, 786–797. https://doi.org/10.1111/j.1460-9568.2010.07322.x (2010).

10. Ferron, J.-F., Kroeger, D., Chever, O. & Amzica, F. Cortical inhibition during burst suppression induced with isoflurane anesthesia. *J. Neurosci.* 29, 9850–9860. https://doi.org/10.1523/JNEUROSCI.5176-08.2009 (2009).

11. Kroeger, D. & Amzica, F. Hypersensitivity of the anesthesia-induced comatose brain. *J. Neurosci.* 27, 10597–10607. https://doi.org/10.1523/JNEUROSCI.3440-07.2007 (2007).

12. Hutchison, R. M., Hutchison, M., Manning, K. Y., Menon, R. S. & Everling, S. Isoflurane induces dose-dependent alterations in the cortical connectivity profiles and dynamic properties of the brain's functional architecture. *Hum. Brain Mapp.* 35, 5754–5775. https://doi.org/10.1002/hbm.22583 (2014).

13. Field, K. J., White, W. J. & Lang, C. M. Anaesthetic effects of chloral hydrate, pentobarbitone and urethane in adult male rats. *Lab. Animals* 27, 258–269. https://doi.org/10.1258/002367793780745471 (1993).
14. Erchova, I. A., Lebedev, M. A. & Diamond, M. E. Somatosensory cortical neuronal population activity across states of anaesthesia. *Eur. J. Neurosci.* **15**, 744–752. https://doi.org/10.1046/j.0953-816x.2002.01898.x (2002).

15. Clowry, G. J. & Flecknell, P. A. The successful use of fentanyl/fluanisone ('Hypnorm') as an anaesthetic for intracranial surgery in neonatal rats. *Lab. Animals* **34**, 260–264. https://doi.org/10.1258/002367700780384771 (2000).

16. Danneman, P. J. & Mandrell, T. D. Evaluation of five agents/methods for anesthesia of neonatal rats. *Lab. Animal Sci.* **47**, 386–395 (1997).

17. Khazipov, R. et al. Early motor activity drives spindle bursts in the developing somatosensory cortex. *Nature* **432**, 758–761. https://doi.org/10.1038/nature03132 (2004).

18. Yang, J.-W.J.W., Hanganu-Opatz, I. L., Sun, J. J. & Luhmann, H. J. Three patterns of oscillatory activity differentially synchronize developing neocortical networks in vivo. *J. Neurosci.* **29**, 9011–9025. https://doi.org/10.1523/JNEUROSCI.5646-08.2009 (2009).

19. Minlebaev, M., Ben-Ari, Y. & Khazipov, R. Network mechanisms of spindle-burst oscillations in the neonatal rat barrel cortex in vivo. *J. Neurophysiol.* **97**, 692–700. https://doi.org/10.1152/jn.00759.2006 (2007).

20. Colonnese, M. T. et al. A conserved switch in sensory processing prepares developing neocortex for vision. *Neuron* **67**, 480–498. https://doi.org/10.1016/j.neuron.2010.07.015 (2010).

21. Sitdikova, G. F. et al. Isoflurane suppresses early cortical activity. *Ann. Clin. Transl. Neurol.* **1**, 15–26. https://doi.org/10.1002/acn3.16 (2014).

22. Chini, M. et al. Neural correlates of anesthesia in newborn mice and humans. *Front. Neural Circ.* **13**, 2243–2252. https://doi.org/10.3389/fncir.2019.00038 (2019).

23. Van Aken, H. & Van Hemelrijck, J. Influence of anesthesia on cerebral blood flow and cerebral metabolism: An overview. *Agressol. revue internationale de physio-biologie et de pharmacologie appliquées aux effets de l'agression* **32**, 303–306 (1991).

24. Alkire, M. T., Hudetz, A. G. & Tononi, G. Consciousness and anesthesia NIH public access. *Sci.* **322**, 876–880. https://doi.org/10.1126/science.1149213.Consciousness (2008).

25. Kaisti, K. K. et al. Effects of surgical levels of propofol and sevoflurane anesthesia on cerebral blood flow in healthy subjects studied with positron emission tomography. *Anesthesiology* **96**, 1358–1370. https://doi.org/10.1097/00000542-200206000-00015 (2002).

26. Liu, X., Li, R., Yang, Z., Hudetz, A. G. & Li, S. J. Differential effect of isoflurane, medetomidine, and urethane on BOLD responses to acute levo-tetrahydropalmatine in the rat. *Magnet. Resonance Med.* **68**, 552–559. https://doi.org/10.1002/mrm.23243 (2012).

27. Munting, L. P. et al. Influence of different isoflurane anesthesia protocols on murine cerebral hemodynamics measured with pseudo-continuous arterial spin labeling. *NMR Biomed.* https://doi.org/10.1002/nbm.4105 (2019).

28. Hara, K. & Harris, R. A. The anesthetic mechanism of urethane: The effects on neurotransmitter-gated ion channels. *Anesth. Analg.* **94**, 313–318. https://doi.org/10.1213/00000539-200202000-00015 (2002).

29. Grinvald, A., Lieke, E., Frostig, R. D., Gilbert, C. D. & Wiesel, T. N. Functional architecture of cortex revealed by optical imaging of intrinsic signals. *Nature* **324**, 361–364. https://doi.org/10.1038/324361a0 (1986).

30. Wyatt, J. S., Delpy, D. T., Cope, M., Wray, S. & Reynolds, E. O. R. Quantification of cerebral oxygenation and haemodynamics in sick newborn infants by near infrared spectrophotometry. *Lancet* **328**, 1063–1066. https://doi.org/10.1016/S0140-6736(86)90467-8 (1986).

31. Obrig, H. NIRS in clinical neurology—A ‘promising’ tool?. *Neuroimage* **85**, 535–546. https://doi.org/10.1016/j.neuroimage.2013.03.045 (2014).

32. Cope, M. The application of near infrared spectroscopy to non invasive monitoring of cerebral oxygenation in the newborn infant. PhD thesis 342 (1991).

33. Strangman, G., Culver, J. P., Thompson, J. H. & Boas, D. A. A quantitative comparison of simultaneous BOLD fMRI and NIRS recordings during functional brain activation. *Neuroimage* **17**, 719–731. https://doi.org/10.1006/nimg.2002.1227 (2002).

34. Huppert, T. J., Hoge, R. D., Diamond, S. G., Franceschini, M. A. & Boas, D. A. A temporal comparison of BOLD, ASL, and NIRS hemodynamic responses to motor stimuli in adult humans. *Neuroimage* **29**, 368–382. https://doi.org/10.1016/j.neuroimage.2005.08.065 (2006).

35. Kocsis, L., Herman, P. & Eke, A. The modified Beer–Lambert law revisited. *Phys. Med. Biol.* https://doi.org/10.1088/0031-9155/51/5/N02 (2006).

36. Ma, Y. et al. Wide-field optical mapping of neural activity and brain haemodynamics: Considerations and novel approaches. *Philos Trans. R. Soc. Lond. Ser. B Biol. Sci.* **371**, 20150360. https://doi.org/10.1098/rstb.2015.0360 (2016).

37. Kozberg, M. G., Chen, B. R., DeLeo, S. E., Bouchard, M. B. & Hillman, E. M. C. Resolving the transition from negative to positive blood oxygen level-dependent responses in the developing brain. *Proc. Natl. Acad. Sci. USA* **110**, 4380–4385. https://doi.org/10.1073/pnas.1212785110 (2013).

38. Kozberg, M. & Hillman, E. Neurovascular coupling and energy metabolism in the developing brain. *Prog. Brain Res.* **225**, 213–242. https://doi.org/10.1016/bs.pbr.2016.02.002 (2016).

39. Minlebaev, M., Colonnese, M., Tsintsadze, T., Sirota, A. & Khazipov, R. Early gamma oscillations synchronize developing thalamus and cortex. *Science* **334**, 226–229. https://doi.org/10.1126/science.1210574 (2011).

40. Sintsov, M., Suchkov, D. D., Khazipov, R. & Minlebaev, M. Developmental changes in sensory-evoked optical intrinsic signals in the rat barrel cortex. *Front. Cell. Neurosci.* **11**, 392. https://doi.org/10.3389/fncel.2017.00392 (2017).

41. Colonnese, M. T., Phillips, M. A., Constantine-Paton, M., Kaila, K. & Jasanoff, A. Development of hemodynamic responses and functional connectivity in rat somatosensory cortex. *Nat. Neurosci.* **11**, 72–79. https://doi.org/10.1038/nn2017 (2008).

42. Sintsov, M. Y., Suchkov, D. & Minlebaev, M. Detection of intrinsic optical signals in the somatosensory cortex of neonatal rats by principal components analysis. *Neurosci. Behav. Physiol.* **48**, 551–556. https://doi.org/10.1007/s11055-018-0598-0 (2018).

43. Grubb, B., Colacino, J. M. & Schmidt-Nielsen, K. Cerebral blood flow in birds: Effect of hypoxia. *Am. J. Physiol. Hear. Circ. Physiol.* **3**, 230–234. https://doi.org/10.1152/ajpheart.1978.234.3.h230 (1978).

44. Chen-Bee, C. H., Agoncillo, T., Lay, C. C. & Frostig, R. D. Intrinsic signal optical imaging of brain function using short stimulus delivery intervals. *J. Neurosci. Methods* **187**, 171–182. https://doi.org/10.1016/j.jneumeth.2010.01.009 (2010).

45. Orliaguet, G. et al. Minimum alveolar concentration of volatile anesthetics in rats during postnatal maturation. *Anesthesiology* **95**, 734–739. https://doi.org/10.1097/00000542-200109000-00028 (2001).

46. Marcar, V. L., Schwarz, U., Martin, E. & Loenneker, T. How depth of anesthesia influences the blood oxygenation level-dependent signal from the visual cortex of children. *Am. J. Neuroradiol.* **27**, 799–805 (2006).

47. Mitrukhina, O., Suchkov, D., Khazipov, R. & Minlebaev, M. Imprecise whisker map in the neonatal rat barrel cortex. *Cereb. Cortex* **25**, 3458–3467. https://doi.org/10.1093/cercor/bhu169 (2015).

48. Hanganu-Opatz, I. L., Ben-Ari, Y. & Khazipov, R. Retinal waves trigger spindle bursts in the neonatal rat visual cortex. *J. Neurosci.* **26**, 6728–6736. https://doi.org/10.1523/JNEUROSCI.0752-06.2006 (2006).

49. Angel, A. & Gratton, D. A. The effect of anaesthetic agents on cerebral cortical responses in the rat. *Br. J. Pharmacol.* **76**, 541–549. https://doi.org/10.1111/j.1476-5381.1982.tb09252.x (1982).

50. Armstrong-James, M. & George, M. J. Influence of anesthesia on spontaneous activity and receptive field size of single units in rat Sm1 neocortex. *Exp. Neurol.* **99**, 369–387. https://doi.org/10.1016/0014-4886(88)90155-0 (1988).

51. Dildy-Mayfield, J. E., Eger, E. I. & Harris, R. A. Anesthetics produce subunit-selective actions on glutamate receptors. *J. Pharmacol. Exp. Ther.* **276**, 1058–1065 (1996).

52. Kimbro, J. R., Kelly, P. J., Drummond, J. C., Cole, D. J. & Patel, P. M. Isoflurane and pentobarbital reduce AMPA toxicity in vivo in the rat cerebral cortex. *Anesthesiology* **92**, 806–812. https://doi.org/10.1097/00000542-20003000-00024 (2000).
53. Franks, N. P. General anaesthesia: From molecular targets to neuronal pathways of sleep and arousal. *Nat. Rev. Neurosci.* **9**, 370–386. https://doi.org/10.1038/nrn2372 (2008).

54. Accorsi-Mendonça, D., Leão, R. M., Aguiar, J. F., Varanda, W. A. & Machado, B. H. Urethane inhibits the GABAergic neurotransmission in the nucleus of the solitary tract of rat brain stem slices. *Am. J. Physiol. Regul. Integr. Comp. Physiol.* **292**, 396–402. https://doi.org/10.1152/ajpregu.00776.2005 (2007).

55. Campagna, J. A., Miller, K. W. & Forman, S. A. Mechanisms of actions of inhaled anesthetics. *New Engl. J. Med.* **348**, 2110–2124. https://doi.org/10.1056/nejmra021261 (2003).

56. Nishikawa, K., Jenkins, A., Paraskevakis, I. & Harrison, N. L. Volatile anesthetic actions on the GABAA receptors: Contrasting effects of a1(S270) and b2(N265) point mutations. *Neuropharmacology* **42**, 337–345. https://doi.org/10.1016/S0028-3908(01)00189-7 (2002).

57. Yamakura, T. & Harris, R. A. Effects of gaseous anesthetics nitrous oxide and xenon on ligand-gated ion channels: Comparison with isoflurane and ethanol. *Anesthesiology* **93**, 1095–1101. https://doi.org/10.1097/00000542-200010000-00034 (2000).

58. Yamakura, T., Bertaccini, E., Trudell, J. R. & Harris, R. A. Anesthetics and ion channels: Molecular models and sites of action. *Annu. Rev. Pharmacol. Toxicol.* **41**, 23–51. https://doi.org/10.1146/annurev.pharmtox.41.1.23 (2001).

59. Proctor, W. R., Mynlieff, M. & Dunwiddie, T. V. Facilitatory action of etomidate and pentobarbital on recurrent inhibition in rat hippocampal pyramidal neurons. *J. Neurosci.* **6**, 3161–3168. https://doi.org/10.1523/jneurosci.06-11-03161.1986 (1986).

60. Hales, T. G. & Lambert, J. J. The actions of propofol on inhibitory amino acid receptors of bovine adrenomedullary chromaffin cells and rodent central neurones. *Br. J. Pharmacol.* **104**, 619–628. https://doi.org/10.1111/j.1476-5381.1991.tb12479.x (1991).

61. Cardoso, R. A., Yamakura, T., Brozowski, S. J., Chavez-Noriega, L. E. & Harris, R. A. Human neuronal nicotinic acetylcholine receptors expressed in Xenopus oocytes predict efficacy of halogenated compounds that disobey the Meyer-Overton rule. *Anesthesiol.* **91**, 1370–1377. https://doi.org/10.1097/00000542-199911000-00029 (1999).

62. Flood, P. & Role, L. W. Neuronal nicotinic acetylcholine receptor modulation by general anesthetics. *Toxicol. Lett.* **100–101**, 149–153. https://doi.org/10.1016/S0378-4274(98)00179-9(ToxicolLett) (1998).

63. Patel, A. J. et al. Inhalational anesthetics activate two-pore-domain K+ channels. *Nat. Neurosci.* **2**, 422–426. https://doi.org/10.1038/8084 (1999).

64. Sceniak, M. P. & MacIver, M. B. Cellular actions of urethane on rat visual cortical neurons in vitro. *J. Neurophysiol.* **95**, 3865–3874. https://doi.org/10.1152/jn.01196.2005 (2006).

65. Patel, A. J. & Honoré, E. Properties and modulation of mammalian 2P domain K+ channels. *Trends Neurosci.* **24**, 339–346. https://doi.org/10.1016/S0166-2236(00)01810-5 (2001).

66. Aguilar, J. R. & Castro-Alamancos, M. A. Spatiotemporal gating of sensory inputs in thalamus during quiescent and activated states. *J. Neurosci.* **25**, 10990–11002. https://doi.org/10.1523/JNEUROSCI.3229-05.2005 (2005).

67. Franks, N. P. & Lieb, W. R. Volatile general anaesthetics activate a novel neuronal K+ current. *Nature* **333**, 662–664. https://doi.org/10.1038/333662a0 (1988).

68. Castro-Alamancos, M. A. Dynamics of sensory thalamocortical synaptic networks during information processing states. *Prog. Neurobiol.* **74**, 213–247. https://doi.org/10.1016/j.pneurobio.2004.09.002 (2004).

69. Clancy, B., Darlington, R. B. & Finlay, B. L. Translating developmental time across mammalian species. *Neuroscience* **105**, 7–17. https://doi.org/10.1016/S0306-4522(01)00171-3 (2001).

70. Flynn, N. M., Buljubasic, N., Bosnjak, Z. J. & Kampine, J. P. Isoflurane produces endothelium-independent relaxation in canine middle cerebral arteries. *Anesthesiology* **76**, 461–467. https://doi.org/10.1097/00000542-199203000-00021 (1992).

71. Farber, N. E. et al. Region-specific and agent-specific dilation of intracerebral microvessels by volatile anesthetics in rat brain slices. *Anesthesiology* **87**, 1191–1198. https://doi.org/10.1097/00000542-199711000-00024 (1997).

72. Chicoine, L. G. et al. Maturational changes in the regulation of pulmonary vascular tone by nitric oxide in neonatal rats. *Am. J. Physiol. Lung Cell. Mol. Physiol.* **293**, L1261–L1270. https://doi.org/10.1152/ajplung.00235.2006 (2007).

73. Martin, C., Jones, M., Martindale, J. & Mayhew, J. Haemodynamic and neural responses to hypercapnia in the awake rat. *Eur. J. Neurosci.* **24**, 2601–2610. https://doi.org/10.1111/j.1460-9568.2006.05135.x (2006).

74. Lahti, K. M., Ferris, C. F., Li, F., Sotak, C. H. & King, J. A. Comparison of evoked cortical activity in conscious and propofol anesthetized rats using functional MRI. *Magn. Reson. Med.* **41**, 412–416. https://doi.org/10.1002/(SICI)1522-2594(199902)41:2%3c412::AID-MRM28%3e3.0.CO;2-3 (1999).

75. Peeters, R. R., Tindemans, I., De Schutter, E. & Van der Linden, A. Comparing BOLD fMRI signal changes in the awake and anesthetized rat during electrical forepaw stimulation. *Magn. Reson. Imaging* **19**, 821–826. https://doi.org/10.1016/S0730-725X(01)00391-5 (2001).

76. Khazipov, R. et al. Atlas of the postnatal rat brain in stereotaxic coordinates. *Front. Neuroanat.* https://doi.org/10.3389/fnana.2015.00161 (2015).

77. Sintsov, M., Suchkov, D., Khazipov, R. & Minlebaev, M. Improved recordings of the optical intrinsic signals in the neonatal rat barrel cortex. *BioNanoScience* **7**, 333–337. https://doi.org/10.1007/s12668-016-0359-x (2017).

78. Zijlstra, W. G. & Buursma, A. Spectrophotometry of hemoglobin: Absorption spectra of bovine oxyhemoglobin, deoxyhemoglobin, carboxyhemoglobin, and methemoglobin. *Comp. Biochem. Physiol. B Biochem. Mol. Biol.* **118**, 743–749. https://doi.org/10.1016/S0305-0491(97)00230-7 (1997).

79. Huss, M. K., Chum, H. H., Chang, A. G., Jampachairsi, K. & Pacharinsak, C. The physiologic effects of isoflurane, sevoflurane, and hypothermia used for anesthesia in neonatal rats (rattus norvegicus). *J. Am. Assoc. Lab. Animal Sci.* **55**, 83–88 (2016).

## Acknowledgements

This work was supported by RSF Grant 16-15-10174 (MM, experimentation and analysis), the subsidy allocated to Kazan Federal University for the state assignment in the sphere of scientific activities 0671-2020-0059 (development of the analytical software) and Fondation de France (consumables). We thank Alexander Yamoldin and Arina Kul'mamet'eva for their participation in the experimental part of this work. We thank Rustem Khazipov and Igor Medyna for their advice and comments on the manuscript.

## Author contributions

V.Sh., V.Si.: Acquisition of data and analysis. I.R.: Development and design of OIS acquisition setup. A.L.: Development of signal contour detection analysis. M.M.: Study concept, supervision and writing.

## 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-021-88461-8.

**Correspondence** and requests for materials should be addressed to M.M.

**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) 2021