Nasal respiration is necessary for ketamine-dependent high frequency network oscillations and behavioral hyperactivity in rats

Changes in oscillatory activity are widely reported after subanesthetic ketamine, however their mechanisms of generation are unclear. Here, we tested the hypothesis that nasal respiration underlies the emergence of high-frequency oscillations (130–180 Hz, HFO) and behavioral activation after ketamine in freely moving rats. We found ketamine 20 mg/kg provoked “fast” theta sniffing in rodents which correlated with increased locomotor activity and HFO power in the OB. Bursts of ketamine-dependent HFO were coupled to “fast” theta frequency sniffing. Theta coupling of HFO bursts were also found in the prefrontal cortex and ventral striatum which, although of smaller amplitude, were coherent with OB activity. Haloperidol 1 mg/kg pretreatment prevented ketamine-dependent increases in fast sniffing and instead HFO coupling to slower basal respiration. Consistent with ketamine-dependent HFO being driven by nasal respiration, unilateral naris blockade led to an ipsilateral reduction in ketamine-dependent HFO power compared to the control side. Bilateral nares blockade reduced ketamine-induced hyperactivity and HFO power and frequency. These findings suggest that nasal airflow entrains ketamine-dependent HFO in diverse brain regions, and that the OB plays an important role in the broadcast of this rhythm.

Ketamine increases HFO power which correlated with fast sniffing. We next examined oscillatory activity within the LFP recordings from the OB (N = 8) in more detail. In control conditions, HFO (130-180 Hz) are often barely visible in spectra (Fig. 2A1). By contrast, gamma oscillations are prominent features of LFP's recorded from the OB and divided into low-  and high-gamma bands   41 . Consistent with our previous reports ketamine produced an almost immediate increase in HFO power lasting around 15 min before returning to baseline (Fig. 2A2) 15,16,35 . Time courses of the effect of intraperitoneal injection of ketamine on low-gamma and high-gamma, and HFO are shown in Fig. 2B1. Analysis of the mean 15 min before and immediately after ketamine revealed a significant increase in HFO power (p = 0.0078 paired-test). Changes in low-and high gamma did not reach significance (p = 0.098, p = 0.20, respectively, paired t-test, Fig. 2B2). Saline injection did not influence gamma or HFO power. Given the striking similarities of the time-courses of HFO, fast sniffing and locomotion we examined if there was a relationship between these measures. Pooled analysis from all rats of the first 15 min post ketamine revealed a significant correlation (Spearman) between HFO power and proportion of fast sniffing (r = 0.3; p < 0.0001) and a stronger correlation for HFO and beam breaks (r = 0.6; p < 0.0001). Analyses from individual rats revealed this correlation was significant for HFO power/fast sniffing in 7/8 rat and HFO power/beam breaks in 8/8 rats (individual R values are shown in Fig. 2C).
Bursts of ketamine-HFO were largely coherent on left and right sides, however, occasional asymmetric bursts could be observed (Fig. 2D1). Closer inspection revealed that although bursts were largely coincident they were not synchronous with troughs rarely in phase between sides (Fig. 2D2). Coherence spectra for left and right OB revealed little coherence between the sides (Fig. 2E) indicating independent HFO generators of this rhythm.
Ketamine-induced HFO is coupled to nasal respiration. Considering ketamine influences nasal respiration and that nasal respiration can entrain oscillations in the OB we investigated this relationship further (analyses based on thermocouple/LFP data; N = 8). Example simultaneous thermocouple recordings, LFPs (1-10 Hz, 70-100 Hz, and 130-180 Hz) and average spectra are shown in Fig. 3A,B. At baseline we observed significant coupling for 70-100 Hz oscillations to the thermocouple signal and local OB oscillations. The driving frequency was 5-8 Hz, which corresponded to fast sniffing (Fig. 3B1). Immediately after injection of ketamine we observed significant coupling to the HFO band at the same driving frequency (Fig. 3B2). Notably, bursts of HFO and gamma occurred at different phases of the respiratory rhythm, with HFO associated towards the peak and gamma on the ascending phase (see polar plots and Supplementary Figs. 1 and 2).
Time-courses showing significant coupling of fast rhythms to 1-10 Hz activity from an example rat are shown, along with mean modulation index (MI) and mean driving frequency for all rats in Fig. 3C-E. After injection of ketamine, respiratory and auto coupling (coupling within the same LFP signal) to high-gamma diminished and instead strong coupling to HFO was observed. As the ketamine effect wore off (around 20 min) coupling returned to pre-ketamine levels ( Fig. 3C1-3,D1-3). In the rats pretreated with haloperidol, injection of ketamine was also associated with HFO coupling but at a slower driving frequency, around 2 Hz and this coupling was stronger than ketamine alone (mean of first 15 min p < 0.05, paired t-test, Fig. 3E1-3).
Ketamine-dependent HFO bursts are coherent across ipsilateral brain regions and modulated by theta frequencies. We examined the relationship between ketamine-dependent HFO in the OB and VS and PFC, since these regions have been shown previously to display increases in HFO power after ketamine 16   www.nature.com/scientificreports/ In a separate series of rats (N = 6) implanted with electrodes in the OB, VS and PFC injection of ketamine was associated with bursts of HFO occurred at similar times in these structures (Fig. 4A1,A2). Analysis of the 130-180 Hz band-pass filtered signal revealed close to zero phase lag between OB and PFC (in one outlier rat HFO was of relatively low amplitude and we observed a greater offset; 0.24 ± 1.1 ms outlier removed, Fig. 4A3). Offsets between OB and VS were more variable (1.1 ± 0.3 ms). Analyses of the raw LFP's revealed ketamine-dependent HFO were largely coherent in the PFC and VS, with respect to OB (Fig. 4A4). Olfactory networks are closely associated with limbic areas 36 . This suggests nasal breathing may influence oscillatory activity in downstream limbic areas. After injection of ketamine significant theta-HFO coupling was observed in all structures (see Fig. 4B1 for an example and mean MI time courses in Fig. 4B2). Two-way ANOVA (group × time F 314,2370 = 4.68, p < 0.0001) revealed coupling was significant stronger in the OB vs. VS and PFC (all p < 0.001). Mean MI and driving frequency ~ 5.5 Hz for first 10 min are shown in Fig. 4B3,B4. www.nature.com/scientificreports/ Unilateral naris blockade reduces ketamine-induced HFO in cortical and subcortical areas. Since nasal respiration and ketamine-HFO are critically linked we tested the hypothesis that unilateral naris occlusion reduces ketamine-HFO power. Rats with electrodes implanted bilaterally in the OB (N = 8) and VS, PFC (N = 7) received unilateral nares blockade. Baseline respiratory rhythm was visible in OB recordings however, although it was not always as pronounced in the raw LFP in VS and PFC channels. Unilateral naris occlusion was associated with almost immediate desynchronization of the raw LFPs on the ipsilateral side ( Fig. 5A1,B1,C1) and served as a reliable marker of naris blockade. Unilateral naris occlusion was associated with diminished HFO power after ketamine on the ipsilateral side, whilst on the control side HFO remained large time × group OB (F 59,590 = 2.52, p < 0.0001), VS (F 59,708 = 2.52, p < 0.0001), PFC (F 59,708 = 2.52, p < 0.0001). Inset histograms shows the HFO modulation index for the first 15 min post ketamine injection and shows that this was also significantly reduced on the ipsilateral versus contralateral side (Fig. 5A2,B2,C2).
Bilateral naris blockade reduces ketamine-induced hyperlocomotion and produces complex effects on HFO. Respiratory rhythms, which were clearly visible in the OB LFPs at baseline were almost completely absent in both OB's after occlusion. Interestingly, we did observe a small increase in the power of a fast-oscillatory rhythm, around 100-120 Hz. In bilateral nares-occluded rats, ketamine injection was associated with an immediate increase in the power of HFO, lasting 1-2 min followed by a return to baseline levels. In controls, ketamine-induced increases in the power of HFO lasted for the standard length of time, 10-12 min. Repeated measures two-way ANOVA revealed a significant time × group interaction (F 79,1106 = 3.79, p < 0.0001). Bonferroni post hoc analyses revealed a significant reduction of HFO power in occluded rats for the first 3.5-9 min post ketamine injection (p < 0.05; Fig. 5D). Unexpectedly, dominant frequency after ketamine was markedly reduced in occluded rats (107.9 ± 2.5 Hz, mean of the first 5 min) after injection compared to 139.3 ± 1.5 Hz in controls (p < 0.001, paired t-test, Fig. 5D insert).
Considering that the power of ketamine-induced HFO and locomotion correlate positively, and that unilateral nares blockade can reduce HFO, we speculated that bilateral nares blockade would impact ketamine-induced locomotion. In rodents, bilateral nares blockade was associated with clear mouth breathing, which was present throughout the course of the recording. In nares-occluded rats, ketamine-induced behavior was qualitatively and quantitatively different compared to controls. Rather than the typical circling behavior, associated with ketamine in controls, we observed an initial ataxic response where the rat would lose balance, 'wobble' , and slowly move its head from side to side. This lasted for around 5 min and was followed by mild locomotor excitation. Compared to controls, this hyperlocomotor response was delayed in occluded rats and significantly less intense. Repeated measures two-way ANOVA revealed a significant effect of group (F 1,79 = 12.12, p = 0.0037) and time × group interaction (F 79,1106 = 11.78, p < 0.0001). Bonferroni post hoc analyses revealed a significant difference for the first 10 min post ketamine injection (p < 0.001, Fig. 5E). Notably, we also examined ketamine-induced locomotion for unilateral naris occlusion and there were no significant differences compared to non-occluded controls (F 1,79 = 0.01, P = 0.9, Supplementary Fig. 3).

Discussion
In control conditions both low-and high-gamma were present in OB LFPs, but clear theta coupling was found only for the high gamma band (around 80 Hz). A finding consistent with observations of others who have associated this coupling with exploration and sensori-motor processing 41,42 . Ketamine increased the power of the faster HFO band which was coupled to theta rhythms. HFO-theta coupling was strongest in the OB, but also present in the VS and PFC. Thus, with respect to fast oscillations, ketamine could be said to have two effects. Firstly, the natural theta-gamma coupling, considered important for neural coding 43 would be derailed. Secondly, the emergence of aberrant HFO which correlated with two key features of the ketamine model, stereotypic sniffing and behavioral hyperactivity.
Abnormal oscillatory activity have been associated with psychiatric diseases, perhaps most notably schizophrenia 44,45 . Ketamine is reported widely to affect brain rhythms in experimental animals and humans 9,10,13 . Gamma rhythms have received a lot of attention due to their role in higher order functioning and increases gamma band activity speculated to represent elevated background electrophysiological noise with translational value for models of schizophrenia 46 . Jones et al. have shown that clozapine, olanzapine, and haloperidol reduce ketamine-dependent increases in ECoG gamma power 47,48 . The HFO band also interacts with antipsychotics, chiefly second generation drugs, but in a different way, reductions in frequency rather than power occur 49 . Spontaneous HFO is extremely weak at baseline and although increased during quiet waking and REM sleep 50 clear functional links remain unknown. After NMDA receptor administration in many brain regions, increases in HFO power tend to be the dominant fast oscillatory change in many brain regions, of awake rodents, and are visible as a clear peak in the power spectra [16][17][18]20,22,51 . It is worth pointing out that ketamine-dependent HFO, although of similar frequency, occur independently of hippocampal ripples 15 . With respect to gamma, in experimental rodents, NMDA receptor antagonists increase gamma power in many brain regions, often as activity broadband rather than a well-defined peak 9,10,47,52 . However, reductions in gamma power have also been observed after ketamine in the nucleus accumbens (bipolar LFPs) 15 , after microinfusion of MK801 to the OB 53 , and also in the present study, indicating that ketamine produces regionally-specific effects on different oscillatory networks.
The mechanisms underlying presence of coherent ketamine-dependent HFO across different brain regions are unclear. This rhythm might be attributed to local HFO generators e.g. in the VS 54 dependent on OB input 35 . Indeed, mitral tufted cells, the main projection neurons from the OB, project to diverse regions and can impose oscillatory activity in their targets (e.g., gamma in the piriform cortex 55 ). Volume-conducted signals can be detected many millimeters from their site of generation 56 and have been proposed to explain the presence of gamma rhythms in VS 57 . In our study, ipsilaterally-recorded ketamine-dependent HFO were coherent but of www.nature.com/scientificreports/ smaller amplitude in the VS/PFC compared to OB. It possible part of this rhythm arises from OB volumeconducted currents. In our previous study we found muscimol inhibition of OB activity reduced but did not completely attenuated ketamine-dependent HFO in the VS 35 . Also, on in the VS, in particular we observed a phase shifts which would be unexpected for volume-conducted signals. Together, these findings suggest that the OB accounts, at least in part, for ketamine-dependent HFO in the VS/PFC. It is possibly that other brain regions that receive powerful projections from the OB, such as the piriform cortex, play a role in ketamine-dependent HFO recorded in corticolimbic areas and further studies are warranted to investigate this. Fast oscillations are ubiquitous in sensory networks and can synchronize spike discharge of neuronal ensembles 58 . Coupling of faster with slower oscillations is considered to enable communication between distant brain areas. In this study we found synchronous bursts of ketamine-dependent HFO in the OB, PFC and VS coupled to theta-driven respiration. Coupling of fast brain rhythms to theta frequencies has been linked functionally to working memory 31 , visual 59 , and olfactory processing 60,61 . Ketamine, which enables local networks to generate HFO more easily, combined with nasal respiration represents a mechanism of synchronization across brain regions, with each respiratory cycle (driven by a wave of stimuli hitting the olfactory system).
Our study shows that ketamine induced a distinct breathing pattern in rats. Unlike natural sniffing, ketamineinduced sniffing was continuous, for almost the complete duration of ketamine's action and appeared stereotyped and purposeless. This finding is consistent with the general increases in stereotypy reported after ketamine and related NMDA receptor antagonists 62 . Increases in nasal respiration are generally linked to locomotor activity and have been associated with exploration 63 , reward expectancy 64 and increased metabolic demand 65 . In line with this, we observed ketamine-dependent increases in sniffing and locomotion were correlated both with each other and the increases in HFO power. Stereotypies and hyperlocomotion are classical dopamimetic behaviors 66 and the dopamine agonist apomorphine has long been known to produce stereotypic sniffing in rodents 67 . By contrast, HFO power is not markedly affected by dopamimetics 16 , including apomorphine (unpublished) or by dopamine receptor blockade 49 , indicating NMDA receptor blockade is the key factor underlying generation of this rhythm. We found the ketamine-related behavioral changes were attenuated by haloperidol, which has important translational value and suggests theta-driven HFO after ketamine, but not HFO generation per se, is dopamine-dependent. However, haloperidol can also act on NR2B receptors 68 and we cannot exclude the possibility that actions at this site contribute to the effect we observed.
Ketamine-dependent HFO correlated with stereotypic sniffing and hyperlocomotion, both widely established measures of psychotic-like behavior in rodents 69,70 . Correlations between HFO power and hyperactivity have been reported by many groups, however a causal link has not been demonstrated. Interestingly, Hansen found that locomotion is not the major driving factor for ketamine-dependent HFO 21 , in line with our previous work 49 (and haloperidol here) showing that reduced locomotor activity does not reduce HFO power. We and others have shown that HFO power is increased in the MAM neurodevelopmental rat model of schizophrenia 18,71 and atypical antipsychotics (clozapine/risperidone) interact with this HFO to reduce its frequency by as much as 80 Hz 49 . These findings point towards relevance of ketamine-dependent HFO and preclinical models of schizophrenia. In this study, we used 20 mg/kg ketamine and it is worth pointing out that other authors have observed ketaminedependent HFO can also emerge after lower doses 18,20,51 . Although the mechanisms of HFO generation is likely to be similar we have not investigated lower doses in the context of nares blockade.
Here, we show that nasal respiration can drive the generation of ketamine-dependent HFO. This was demonstrated using simultaneous thermocouple and LFP recordings, and unilateral naris occlusion which was associated with unilateral reductions in HFO on the ipsilateral side. Ketamine-dependent faster stereotypic nasal respiration was associated with increased HFO power, but also slow basal respiration modulated HFO power in the presence of haloperidol. However, in the current study we did not examine the effect of low doses of ketamine, for example 3 mg/kg can produce weak increases in HFO without markedly altering locomotion 51 . Further studies are warranted to document the effect of other psychoactive compounds such as LSD or DOI which can increase HFO power without affecting locomotion 21 .
To date, the majority of studies examining NMDA receptor-dependent HFO and gamma rhythms have focused on their relevance to the NMDA hypofunction model of schizophrenia, which is supported by findings that both rhythms interact with antipsychotics and are changed in other models of schizophrenia. Several studies Bilateral naris occlusion was associated with a significant reduction in HFO power compared to the control (non-occluded) state (p < 0.05, 2 way-ANOVA, Bonferroni post hoc). Although the power of spontaneous HFO was equivalent in control and occluded groups, post occlusion unexpectedly we observed a small but significant increase in the power of HFO. Inset, shows the HFO frequency (initial 5 min period post ketamine) which was significantly reduced after bilateral nares occlusion. (E) Time courses of beam breaks after ketamine for bilateral naris occluded and control (non-occluded) conditions. Bilateral naris occlusion was associated with a significant reduction in beam breaks which was also notably delayed occurring several minutes after ketamine injection. Bilateral naris occlusion did not influence beam break activity prior to ketamine injection.

Scientific Reports
| (2020) 10:18981 | https://doi.org/10.1038/s41598-020-75641-1 www.nature.com/scientificreports/ have shown gamma power increases in humans after ketamine, but frequencies above the gamma band are notoriously hard to record without invasive techniques. However, there is some evidence from MEG that ketaminedependent HFO does occur in the human brain. Intriguingly, clinical reports have shown that the psychosislike actions of ketamine correlate with later antidepressant effects 7 and can even predict clinical responses in treatment-resistant depression 72 . The short-lasting changes in ketamine-dependent oscillations (around 15 min) may represent this initial stage of ketamine's actions which lead to longer-term plastic changes. HFO in the brain tends to denote spiking activity 73 and thus HFO synchronized at theta frequencies by respiration would be expected to reflect acute widespread firing of projection neurons. Indeed, ketamine-dependent HFO, at least in the OB, associates with the firing of projection neurons which appears essential for relay to limbic regions 35 . Electrical high frequency stimulation is an emerging method used for the treatment of depression with effective frequencies commonly 130 Hz 74 strikingly similar to ketamine-dependent HFO. Against the notion that HFO reflect direct/indirect therapeutic effects, Zanos has shown that the antidepressant actions of the ketamine metabolite (2R,6R)-HNK are independent of NMDA receptor blockade and devoid of ketamine-related side effects 75 . These data would argue that NMDA receptor-dependent actions of ketamine, including HFO, are not required for antidepressant effect, but rather may be associated with adverse drug effects. Thus although the functional relevance of ketamine-dependent oscillatory changes, whether therapeutic or adverse, remain to the explained the presence of fast oscillatory activity in the brain would be predicted to have physiological consequences. For example, in vitro electrical high frequency stimulation (HFS > 100 Hz) has long been known to produce long-term potentiation (LTP) 76 . Potentially more relevant to the ketamine model is that theta-burst high frequency stimulation also produces LTP 77 possibly in a more robust way 78 .
It is also worth noting that bulbectomy is a widely used rodent model of depression which produces downstream changes in the same limbic circuits as affected in patients 36 . Ketamine-HFO are dependent on OB activity 35 and this issue gains important since intranasal delivery is the preferred route of administration for ketamine and its analogues, in depression. Thus, OB effects appear important which may work in parallel with other networks, for example Yang showed recently that ketamine inhibited burst activity in the habenula an "antireward" area 79 .
Although an OB-habenula pathway is believed to exist 80 it is unknown if these ketamine-dependent changes are related.
In summary, as noted by Moberly et al., voluntary or involuntary changes in nasal breathing modulate brain activity 26 . Olfactory sampling involves multiple brain regions and coherent oscillatory activity across areas may serve as a mechanism for temporal coordination of cortical and limbic networks 36 . We conclude that nasal airflow is necessary for the emergence of ketamine-dependent HFO in multiple brain regions and behavioral hyperactivity. Zarate has proposed that mechanistic similarities may exist between ketamine-induced depersonalization and antidepressant response 72 . We speculate that ketamine's HFO represent converging electrophysiological activity which could account for initial psychotic-like effects and later antidepressant effects and further studies are warranted to address this issue.

Methods
Surgery: 25 male Wistar rats (250-350 g) were used in this study. Group 1 (thermocouple study): twisted stainless steel electrodes (125 µm, Science Products, Germany) were implanted in the OB (AP + 7.5, ML ± 0.5, DV 3-3.5 mm) along with two precision fine bare wire temperature sensors (80 µm diameter, 5TC-TT-KI-40-1M, Omega Engineering Inc., Czech Republic) into the right and left nasal cavity for monitoring nasal airflow (N = 8). These rats were used to examine the effect of 20 mg/kg ketamine on nasal respiration and oscillatory activity in the OB. A subgroup (N = 5 rats) were used to examine the effect of 1.0 mg/kg haloperidol pretreatment. Group 2 (multistructure study): stainless steel electrodes were implanted in the OB, PFC (AP + 3.2, ML 0.5, DV 3.0), and VS (AP 1.6, ML + 1.0, DV 7.0) for simultaneous recording of LFPs used to examine HFO in multi-structures simultaneously (N = 6). Two rats had electrodes bilaterally implanted in the OB, PFC, VS and were also used in naris blockade experiments (group 3 and 4). Four rats had additional electrodes targeted to the hippocampus and/or amygdala and piriform cortex for a supplementary experiment (ketamine-dependent HFO was inconsistent in these areas and the data not shown). Group 3 (unilateral naris blockade in PFC and VS): five rats were implanted with twisted stainless steel electrodes bilaterally in the PFC and VS to determine the effect of unilateral naris blockade on these regions. Additionally, 2 rats from group 2 (bilaterally implanted in OB, PFC and VS) received naris blockade giving a total N = 7. Group 4 (unilateral and bilateral naris blockade in OB): stainless steel electrodes were implanted bilaterally in the OB of 6 rats for unilateral and bilateral naris occlusion experiments. Two rats from group 2 (bilateral implantation in OB, PFC, and VS) also received naris blockade giving a total N = 8. A screw posterior to the bregma was used as a reference/ground in all cases.
One week after surgery, rats were placed in an arena (44 × 50 × 42 cm). LFPs and thermocouple recording were recorded through a JFET preamplifier, amplified 1000×, filtered 0.1-1000 Hz (A-M Systems, USA), digitized at 5 kHz (Micro1401, CED, Cambridge, UK). Horizontal locomotor activity was assessed by photocell beam breaks (Columbus Instruments, USA). Thermocouple experiments were performed according to the Latinsquare design, whereby each rat was injected twice in a pseudorandomized order with either ketamine 20 mg/ kg (Sigma, Poland) or saline. A subgroup of rats (N = 5) were preinjected with 1 mg/kg haloperidol followed 15 min later by 20 mg/kg ketamine. Effect of 20 mg/kg ketamine on LFP oscillations was examined in rats with electrodes implanted in the OB, VS, PFC (N = 6). For naris blockade experiments occlusion was achieved using a silicon occluder 81 either unilaterally for rats implanted with electrodes bilaterally in the VS and PFC (N = 7). Rats with electrodes bilaterally in the OB received both unilateral and bilateral nares occlusion (OB; N = 8). For naris blockade experiment, rats were baselined for 20 min and then briefly anesthetized using isoflurane to allow insertion of the silicon occluder(s). Rats were recorded for a further 60 min to ensure sufficient time for isoflurane washout, and then injected with ketamine 20 mg/kg. Electrode locations were determined on 40 µm Analyses. Mean power spectra of the LFP were computed on successive 30 s data blocks using a fast Fourier transform (4096 points) to calculate integrated, dominant power and frequency of low-gamma (30-60 Hz), high-gamma (70-100 Hz) and HFO (130-180 Hz). Thermocouple signals (1-10 Hz) were used to determine the dominant sniffing frequency and the proportion of fast (4-10 Hz) sniffing behavior. Waveform correlations of 130-180 Hz bandpass filtered LFPs ipsilateral (OB, PFC, VS) and contralateral OB were computed and maximum correlation and offset were used to calculate similarity and synchronicity. Coherence was computed using the coherence function from Python script library to estimate the magnitude squared coherence of discrete-time signals. In this method power spectral density is estimated with Welch's method. Coherence was computed shortly after ketamine injection (1 min). To establish phase difference between HFO bursts we searched for maximum cross-correlation score between two signals in range of ± 5 ms (50 time samples), which covers one cycle of the HFO wave.
Phase-amplitude coupling analysis. Cross frequency phase-amplitude coupling was examined between the phase of oscillations in the thermocouple signal or low-frequency (1-10 Hz) LFP, and amplitude of high-frequency LFP oscillations (https ://githu b.com/Gabri elaJu rkiew icz/ePAC). It enables analysis of coupling with the Modulation Index (MI) method 82,83 which relies on measuring the distance between the obtained distribution of high-frequency amplitude across low-frequency phases from the uniform one. The preprocessing of the LFPs and thermocouple data consisted of high-pass filtering (cutoff 0.1 Hz), lowpass filtering (cutoff frequency 250 Hz; both Butterworth 2nd order), downsampling to 625 Hz, and dividing data into 20 s fragments. The investigated low-frequency band ranged from 2 to 10 Hz (1 Hz step and 2 Hz filtration bandwidth) and high-frequency band ranged from 30 to 200 Hz (5 Hz step). Remaining parameters for the ePAC toolbox were set to default values (nbCycles = 3, w = 5, nbBins = 18, Nboot = 200, peMI = 95, pPhaseCom = 95, Athresh = 0.1).
Statistical analysis. Data are expressed as mean ± SEM or median ± interquartile range according to data normality. For multiple group analyses data were analyzed using repeated-measures ANOVA followed by the Bonferroni post hoc test or the Kruskal-Wallis test followed by Dunn's multiple comparison test. Spearman's rank correlation coefficient was used to examine the relationship between sniffing, locomotion and oscillatory activity. p < 0.05 were considered statistically significant.