{'fraction_non_alphanumeric': 0.033218184330523697, 'fraction_numerical': 0.015918198148404034, 'mean_word_length': 4.95360315893386, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Prefrontal and hippocampal theta rhythm show anxiolytic‐like changes during periaqueductal‐elicited “panic” in rats

Abstract Anxiety and panic are both elicited by threat and co‐occur clinically. But, at the neural level, anxiety appears to inhibit the generation of panic; and vice versa. Anxiety and panic are thought to engage more anterior (a) and mid‐posterior (m) parts of the periaqueductal gray (PAG), respectively. Anxiety also engages the hippocampus and medial prefrontal cortex. Here, we tested if mPAG but not aPAG stimulation would suppress prefrontal and hippocampal theta rhythm as do anxiolytic drugs. Twelve male rats with implanted electrodes were stimulated alternately (30 s interval) in the left PAG or right reticular formation (reticularis pontis oralis [RPO]—as a positive control) with recording in the left prelimbic cortex and left and right hippocampus. PAG stimulation was set to produce freezing and RPO to produce 7–8 Hz theta rhythm before tests lasting 10 min on each of 5 days. mPAG stimulation decreased, and aPAG increased, theta power at all sites during elicited freezing. mPAG, but not aPAG, stimulation decreased prefrontal theta frequency. Stimulation did not substantially change circuit dynamics (pairwise phase consistency and partial directed coherence). Together with previous reports, our data suggest that panic‐ and anxiety‐control systems are mutually inhibitory, and neural separation of anxiety and panic extends down to the aPAG and mPAG, respectively. Our findings are consistent with recent proposals that fear and anxiety are controlled by parallel neural hierarchies extending from PAG to the prefrontal cortex.

suppress prefrontal and hippocampal theta rhythm as do anxiolytic drugs. Twelve male rats with implanted electrodes were stimulated alternately (30 s interval) in the left PAG or right reticular formation (reticularis pontis oralis [RPO]-as a positive control) with recording in the left prelimbic cortex and left and right hippocampus. PAG stimulation was set to produce freezing and RPO to produce 7-8 Hz theta rhythm before tests lasting 10 min on each of 5 days. mPAG stimulation decreased, and aPAG increased, theta power at all sites during elicited freezing. mPAG, but not aPAG, stimulation decreased prefrontal theta frequency. Stimulation did not substantially change circuit dynamics (pairwise phase consistency and partial directed coherence). Together with previous reports, our data suggest that panic-and anxiety-control systems are mutually inhibitory, and neural separation of anxiety and panic extends down to the aPAG and mPAG, respectively. Our findings are consistent with recent proposals that fear and anxiety are controlled by parallel neural hierarchies extending from PAG to the prefrontal cortex. There is a reason to see the archicortex of the hippocampus as a key mediator between these higher and lower levels. Papez (1937) saw the hippocampus is a key node for processing emotion (see also Vann & Nelson, 2015;McNaughton & Vann, 2022). Then, the hippocampus was seen as more dedicated to memory. But there is a evidence that hippocampal damage alters emotion (Davidson et al., 2005(Davidson et al., , 2009Tracy et al., 2001) and is a key control center for stress and hormones (Lathe, 2001(Lathe, , 2004. The ventral hippocampus has been most clearly linked to anxiety (Bannerman et al., 2004); but anxiolytics impair learning in the Morris water maze (McNaughton & Morris, 1987, 1992 and reliably reduce hippocampal theta rhythmicity (McNaughton et al., 2007) on which spatial learning depends (McNaughton et al., 2006). Here, we investigate the interactions of basic mechanisms of fear and anxiety with hippocampal and prefrontal theta rhythmicity.
In rodents, the most basic expression of high levels of fear (e.g., generated by a proximal predator) involves panic in the form of true freezing (i.e., not sensitive to anxiolytics; Blanchard & Blanchard, 1990a) if escape is not available; and in the form of running if escape is available (Blanchard & Blanchard, 1990b). Controlled freezing and undirected running are elicited simply by stimulation of the periaqueductal gray (PAG; De Molina & Hunsperger, 1962), which contains fundamental control centers for aversive and defensive responses (Silva & McNaughton, 2019). If the states induced by activation of the PAG are reflective of panic, then understanding them should illuminate and broaden our understanding of general fear-and specific panic-generating systems in the brain 1 and of how to treat their disorders in the clinic.
Increasing dorsal PAG stimulation intensity evokes a hierarchy of behavioral and autonomic responses like that elicited by increasing predatory threat (Evans et al., 2018). Moderate stimulation elicits freezing while stronger intensities elicit undirected escape responses consisting of running and/or jumping, similar to animals close to predators (Sudre et al., 1993). Pharmacological studies with panicolytic drugs show that dorsal PAG stimulation-elicited escape behavior is a reliable proxy for fear and panic in animals (Jenck et al., 1995). However, there has been less attention to dorsal PAG-induced freezing.
In contrast, the most basic expression of high levels of anticipatory anxiety (e.g., generated by the smell or memory of a predator) is a form of anxiolytic-sensitive immobility-a defensive quiescence that can involve distinctive postures such as "stretch attend" (Blanchard & Blanchard, 1990a). Particularly in operant conditioning to aversive stimuli, anticipatory immobility is often called "freezing." But unlike panic-related freezing, this anticipatory "contextual" form of freezing (whether innate or conditioned) is responsive to anxiolytic drugs (Conti et al., 1990) and is a model for human anxiety-related disorders in general (Tovote et al., 2016) and generalized anxiety in particular (Luyten et al., 2011).
Fear and anxiety appear to be mutually antagonistic. Defensive quiescence requires inhibition of any co-activated tendency to escape; while escape requires inhibition of both freezing and defensive quiescence. Conditioned anxiety that generates defensive quiescence inhibits escape-like responses induced by PAG activation (Magierek et al., 2003). Conversely, in humans with panic disorder, anxiety reduction can release panic attacks (Cohen et al., 1985;Klein, 1993;Mellman & Uhde, 1989). Furthermore, in rodents and humans, anxiety enhances pain sensitivity while fear reduces it (Rhudy & Meagher, 2000); and anxiety and fear are affected in opposite directions by the 5-HT system (Deakin & Graeff, 1991).
Here, we investigate forebrain local field potential (LFP) activity during immobility elicited by PAG stimulation. We contrast PAG stimulation with reticular stimulation that elicits hippocampal (HPC) theta rhythm, which is a reliable model for testing anxiolytic drug efficacy (McNaughton et al., 2007). Although reticular stimulation does not elicit an emotional response, HPC theta appears to be a key mechanism controlling anxiety ; so we expect that increases and decreases in HPC theta will reflect respective changes in anxiety systems in areas such as prefrontal cortex (LeDoux, 2015) that can become entrained to HPC theta during activities such as risk assessment (Adhikari et al., 2010;Gordon, 2011;Padilla-Coreano et al., 2019;Young & McNaughton, 2009).
We also test if the anterior PAG can generate distinct forebrain responses from more posterior PAG, as previously anticipated (Silva & McNaughton, 2019). Although stimulation of both parts of the PAG evokes immobility, anatomical, biochemical, and functional evidence suggests that these two PAG regions are not similar. Could the immobility elicited by different levels of the PAG be electrophysiologically different (reflecting defensive quiescence versus true freezing)? Our hypothesis of panic-anxiety mutual antagonism predicts that electrically stimulating the middle PAG freezing would reduce forebrain theta while anterior PAG freezing would increase forebrain theta.

| Animal housing and handling
The subjects were 15 male Sprague-Dawley rats obtained from the University of Otago Hercus-Taieri Resource Unit and weighed between 150 and 200 g on arrival. Two were removed from the study due to health complications, and a third because it had stimulating electrodes located outside of the regions of interest. Data below are reported for the remaining 12 rats, divided, after histology, into two groups with stimulating electrodes in the anterior (N = 5) or midposterior (N = 7) PAG, respectively. On arrival, rats were caged in pairs, and given ad libitum access to food and water in a temperaturecontrolled room (20-22 C), on a 12 h light-dark cycle with lights on at 6 a.m. All data were collected during the light period of the cycle.
Rats had 10 days to acclimatize to the laboratory with daily handling, and were scheduled for surgery once body weight was >300 g. All experiments were conducted in accordance with the New Zealand Animal Welfare Legislation, approved by the University of Otago Animal Ethics Committee (approval numbers 29/12 and 10/16).

| Local field potential recording and stimulation electrodes
The animals used in the current experiments were intended for use in a separate drug experiment not reported here and therefore were 1 Throughout this article, we distinguish anxiety and fear as detailed in Gray and McNaughton (2000; for updated neurology see figure 12 in Silva & McNaughton, 2019) and so include the specific reactions of "panic" (and their control by the PAG) within the more general reactions of "fear" (and their control by a hierarchy of structures running from the PAG at the lowest level to prefrontal cortex at the highest). This differs at key points from alternative uses of these words, for example, Perusini and Fanselow (2015). implanted with one of two different types of recording electrode arrays targeting the left prelimbic cortex (PrL) and the left and right hippocampus. All animals had bipolar stimulating electrodes implanted in the region of the right reticularis pontis oralis (RPO) and left PAG and a guide cannula for intracerebral drug injections aimed at the midline nucleus incertus. All testing reported here was conducted during an initial saline injection session before any drug injections had been tested. Recordings were referenced to a common screw electrode behind lambda and, during post-processing, subtraction of an adjacent pair within an array (based on histology) delivered a local "bipolar" signal for analysis. The ground was a length of bare silver wire (0.005 00 diameter, AM Systems), and was wrapped around the outside of the implant prior to suturing. Both RPO and PAG stimulating electrodes were 0.005 00 bipolar twisted stainless-steel wires (AM-Systems) with a 0.5 mm tip separation.
One type of recording array was multipolar (up to 11 electrodes per array; each electrode made from 0.001 00 diameter nichrome wires, California Fine Wire Company) connected to a 32-pin mini array socket (Ironwood). For the HPC arrays, eight recording tips were spaced 200 μm apart to cover the dorso-ventral axis. For PrL, 11 tips were spaced 500 μm apart to cover most of the medial wall of the prefrontal cortex. For further details see Banstola et al. (2021).
The other recording type was bipolar or tripolar and built from twisted, PFA-insulated stainless steel wires with a 0.005 00 thickness (AM-Systems) with gold pins that were inserted into a McIntyre connector (McIntyre & Molino, 1972). For the HPC arrays, two twisted wires with a tip separation of 0.6 mm were aimed at the hippocampal fissure and stratum oriens of CA1. For PrL, three twisted wires with a tip separation of 1 mm were aimed to span the dorso-ventral extent of the PrL.

| Surgical methods
Thirty minutes prior to the beginning of surgery all animals received a prophylactic dose of antibiotic Amphoprim (30 mg/kg). The rats were given ketamine and medetomidine (75 mg/kg and 0.5 mg/kg, s.c.) as anesthetics for surgery and given additional injections of ketamine if necessary, as determined by toe pinch reflex. As soon as voluntary movement ceased, atropine (0.05 mg/kg, s.c.) was administered to aid breathing. For inhalation anesthesia, an isoflurane vaporizer (E-Z Anesthesia) was used. An oxygen-isoflurane mixture was delivered at a rate of 1.5 L/min. The percentage of isoflurane in the mixture was regulated between 2% and 3% depending on the animal's toe-pinch reflex and breathing rate. Buprenorphine (0.1 mg/kg, s.c.) was administered for additional pain control.
The animal's scalp was then shaved, sterilized, and infiltrated with a combined solution of Marcaine (0.4 ml/kg) and Lidocaine (0.2 ml/kg) subcutaneously. Tricin (Jurox, Australia) was applied to the animal's eyes, which were also covered with a wet cotton pad for protection.
The animal's head was secured to the stereotaxic frame using nontraumatic ear bars. During surgery, all animals rested on a heated pad in order to maintain body temperature.
After placement of the rat in the stereotaxic frame, a single incision was made in the scalp, the skin was retracted and the skull exposed, where holes for the electrodes were drilled. The coordinates for implantation are shown in Table 1. These animals also had a cannula implanted at the nucleus incertus experiments reported here.
The dorsal-ventral coordinates used for implanting the guide cannula (Plastic One) were calculated so the tip would be 1 mm dorsal to the nucleus incertus.
Seven stainless steel screws were fixed to the periphery of the skull in order to secure the electrodes and cannula. The reference electrode consisted of an eighth screw posterior and lateral to lambda.
A bare silver wire tightly wrapped around outside the screws acted as earth. Dental cement held the connector in position.
For animals that received injectable anesthesia, Antisedan (atipamezole, 2.5 mg/kg s.c.) was administered to reverse the effects of the medetomidine to facilitate recovery. For animals that received gaseous anesthesia, the isoflurane-oxygen ratio was dropped to 0% and animals received pure oxygen for 10 min. In both cases, carprofen (5 mg/kg, s.c.) was given to relieve post-operative pain. During surgery, 2 ml of 0.9% physiological saline (up to 10 ml in total) was periodically administered s.c. to maintain hydration.
The animals were then individually housed in a clean cage with free access to food and water. The animal's well-being and recovery were monitored for 7 days according to the University of Otago Animal Welfare guidelines using a standard monitoring sheet. After approximately 10 days, rats were grouped back into their original pairs.

| Behavioral apparatus
All animals were tested in a modified operant chamber (width: 24.5 cm, height: 18.5 cm, and depth: 30 cm), with its levers, lights, and feeding mechanisms removed. All walls were transparent acrylic except one, which was aluminum; the floor was composed of cylindrical metal bars with a tray underneath. This operant box was placed inside a second chamber, which apart from a small hole on the ceiling through which the recording cable ran, was fully sealed from external light and noise. Illumination inside the chamber was provided by a red light, and a fan provided air recirculation and background noise. A small video camera inside the larger chamber provided a live video feed for experimenter observation.
Prior to any data collection, the rats were individually placed inside the observation box and allowed to explore for 50 min for five consecutive days. Before each exposure, the walls and floor of the observation box were cleaned with a 10% ethanol solution and allowed to dry.

| Electrical stimulation response thresholds
Stimulating currents for the PAG and RPO electrodes were delivered by a custom-built, constant current isolated stimulator controlled by custom software programmed in Visual Basic 6. Stimulation was set at 100 Hz with monophasic pulses with a width of 0.1 ms. Trains for the RPO stimulation were set to 1 s duration, and PAG trains were set to 8 s maximum duration.
In order to determine the ideal currents necessary to drive the desired responses from the PAG and RPO, rats were tested for 2 days before data acquisition. They were individually placed in the observation chamber and given 30 min to habituate with the tethers before testing.
First, the current for appropriate RPO responses was established.
When animals were immobile, an initial current of 10 μA was delivered in a 1-s train. High-intensity currents can provoke undesired motor reactions in free-moving animals, so during the delivery of the current, the animal's behavior and real-time LFPs were observed. We defined the ideal response as one where repeated stimulation elicited clear HPC theta during the stimulation train with no evoked motor response. If no induced theta or motor responses were observed, the stimulating current was increased in 5 μA steps until the highest current that caused no motor responses was reached.
Following the determination of RPO-stimulation values, the animals were tested for PAG-evoked responses. The goal for the PAG stimulation was to establish the minimum current intensity necessary to induce the emotional response of freezing. Stimulation trains were 8 s long, and the first current was set at 10 μA, with increases in 5 μA steps until the desired freezing behavior was evoked.
After repeated PAG stimulations at the freezing threshold, RPO stimulations were conducted again in order to verify that the desired hippocampal responses were still being elicited. In all cases, the currents established for the RPO-stimulation before PAG activations were still effective in inducing hippocampal theta rhythmicity when tested while not moving.

| Data acquisition and processing
Local field potentials were acquired by a Micro1401 (CED) at a sampling rate of 512 Hz. The signal was passed through Grass Model 15 amplifiers, and hardware gain was set at Â5000 for all acquired data. The signals acquired from the recording electrodes were referenced to the recording screw positioned at the posterior part of the rat's skull behind lambda in an electrically silent zone-but note that, for analysis, local bipolar re-referencing was used as described below.
For data processing, the files were down-sampled to 256 Hz using the in-built cubic spline interpolation in Spike2 (CED), bandpassed at 2-30 Hz and exported to Matlab (Mathworks).
Referencing from a common, distant electrode has been shown to bias the interpretation of signals and to cause spurious correlations (Lalla et al., 2017;Shirhatti et al., 2016). In order to minimize this, subtractions of the signal in the PrL and HPC were performed between local pairs of electrodes. Based on the histological reconstruction, the electrode centered in each region of interest (left HPC, right HPC, and PrL) was re-referenced to its closest neighbor. All data analysis was then performed on the three resulting localized signals from the left HPC, right HPC, and PrL, respectively.

| Recording sessions
Rats were tested for 10 min at the beginning of a single recording session. They were consecutively stimulated at the RPO and PAGfreezing thresholds established before. The first 1-s RPO train was delivered, followed 30 s later by an 8-s PAG train (Figure 3b), and at the end of the recording period, the rat had received 10 stimulations of each type ( Figure 3a).

| Data processing
Given that PAG stimulation periods are eight times longer than RPO stimulation, 1-s segments of EEG were extracted from each stimulation event in order to adequately compare the effects of the two stimulation types. One-second periods of data before the onset of each event were extracted from each epoch, and treated as a control "prestimulation period," and the difference between this and the 1-s periods after the onset "during stimulation" was used for assessing stimulation effects ( Figure 3b).
After segmentation of the data, each 1-s epoch was individually assessed for power, instantaneous frequency, coherence, pairwise phase consistency (PPC), and partial directed coherence (PDC). All 10 trials of each type within the session were then averaged for each measure within the "hippocampal theta" range (5-10 Hz).

| Spectral analysis
The time-frequency analyses reported here were done with two 1-s windows with 500 ms overlap between them that contains a trailing or leading 250 ms after or before the start of stimulation, respectively.
All spectral analyses were carried out using the multi-taper method, with three tapers and a numerical bandwidth product of 5 from the Chronux package (http://chronux.org/; Bokil et al., 2010;Mitra, 2007) and with a cut-off frequency set to 30 Hz. Note that the primary (central) taper of the three constitutes a Hanning window (i.e., cosine wave) that captures maximum power in the central 500 ms and little power in the leading/trailing 250 ms on either side. The secondary and tertiary tapers extend the width of this capture to some extent.
Averaging across the three provides important smoothing of the spectral density estimation and coherence estimation. Here, we defined "theta" as 5-10 Hz and z-scored power spectral density for this band was calculated relative to the rest of the frequency spectrum.

| Pairwise phase consistency analysis
Supplementary to coherence analysis, PPC was also applied. PPC calculates the distribution of the relative angular distances between two Hilbert-transformed signals and estimates the phase coupling from the two based on the outcome. Unlike spectral coherence or phaselock value estimates, PPC has been demonstrated to be free of the biases that come from a small number of observations (Vinck et al., 2010)-that is, a tendency for phase-lock values to not only be overestimated as sample size decreases but for the overestimates to converge at lower N (their figs. 2 and 4). In numerical terms, the PPC is equivalent to the square of the phase-locking value or coherence (and so is directly related to proportion of variance accounted for).

| Partial directed coherence analysis
To assess the direction of theta modulation between HPC and PrL, we employed the PDC method developed by Baccala and Sameshima (2001) and as used previously in our laboratory with hippocampal recordings Young et al., 2017)-also see Young and Eggermont (2009) for discussion.

| Experimental design and statistical analysis
Each measure of interest had its variance distribution analyzed for normality. The correlations between the means and the variances of the samples were calculated, and the data were power transformed, when necessary, before repeated measures ANOVAs were performed. For the ANOVAs, orthogonal polynomial contrasts were extracted for dimensional factors using SPSS (IBM) with the structure of some organized to test specific a priori hypotheses (see below). The interactions between these contrasts were assessed, and the linear and quadratic components were derived, where the significant interactions between the different levels of treatments can be captured by significant polynomial trends.
The use of these methods allows for 1 degree-of-freedom contrasts, which prevent the effects of sphericity; 1 degree-of-freedom contrasts are also statistical measures that are recommended in the Publication Manual of the American Psychological Association (Kirk, 2013

| Histological methods
At the end of all testing, the rats were euthanized with pentobarbi-

| Recording electrode placements
All rats had recording electrodes located in the left PrL and in the HPC, as well as stimulating electrodes located in the PAG and RPO.
Histological locations for recording electrodes are detailed separately for rats with aPAG stimulating electrodes and mPAG stimulating electrodes (gray dots vs. black , respectively in Figure 2). 3.2.2 | Effects on local theta oscillations and nondirected connectivity Figure 4 shows both raw spectral measures for pre-stimulation and during stimulation and the differences between real and pre-stimulation, that is, the effect of the stimulation. There were no major differences between mPAG and RPO during the pre-stimulation periods (i.e., immediately before stimulation) and so the text below focuses on the difference scores.
Importantly, neither the power nor frequency reductions in PrL were accompanied by any coherence reduction.

Theta pairwise phase consistency
The left and right HPC were generally more phase coherent with each other than with the PrL (region [quadratic], F(1,5) = 7.61, p = .02). No differences in theta PPC were seen between the two hippocampi and the PrL, that is, the left HPC shows similar theta phase consistency with the PrL to that of the right HPC (region [linear], F(1,5) = 0.02, p = .87). There were no significant differences between the prestimulation and stimulated periods with either mPAG or RPO stimulation.

| Effects on directed connectivity (partial directed coherence)
To assess the dynamics of causal influence between the HPCs and PrL through theta oscillations, PDC analysis was conducted. Figure

| Comparison of the effects of anterior periaqueductal gray with middle periaqueductal gray stimulation
In this section, we assess whether PAG stimulation of the five excluded rats with more anterior (aPAG) stimulating electrodes produced the same results as for the mPAG rats. All the results are reported as a direct comparison of aPAG with mPAG. Note that the    The first of each pair of panels shows the raw data and the second shows the difference between during stimulation and pre-stimulation, that is, the effect of stimulation. The observed effects and differences are detailed in the text. Here and in the subsequent figures, statistically significant effects are illustrated with dashed lines and curves. Where a main effect applies to sets of data with no higher order interaction, this is illustrated with the same size curve for each included case, the curves are linked with straight lines and only a single p value is indicated (as in panel d). Separate post hoc effects (panels a, b) have unconnected curves each with their own p indicator. p values are indicated as *p < .05 or **p < .01. PAG, periaqueductal gray; PPC, pairwise phase consistency; RPO, reticularis pontis oralis raw mPAG data are the same as in Section 3.2 and are repeated for ease of graphical comparison and understanding of the difference scores.

| Behavioral and neural effects of stimulation: overview
The behaviors elicited by stimulation in the sites were partially distinct.
As noted in Section 3.2.1, mPAG animals froze during stimulation, with occasional head scanning and faster breathing in some animals. aPAG stimulation also produced a form of freezing, but with a more consistent and intense occurrence of higher rates of breathing. aPAG stimulation also produced strong sniffing and occasional facial contractions during stimulation, which were completely absent in the mPAG rats.
With aPAG, the average minimum electric current that caused freezing was 428 μA, with a variation of ±6.9 μA (SEM) in different days of testing. This is considerably higher than with mPAG (68 μA, ±6.5 μA; t [10] = 3.88, p = .003). As noted above, higher levels of mPAG group stimulation generated escape behaviors, characterized by forward running in circles while the stimulation train was on. aPAG rats were more resistant to change of behaviors with higher currents; most consistently failing to show escape with stronger stimulation and usually remain frozen. Occasionally, higher currents in the aPAG group would cause backward locomotion, with an exaggerated arched back and lowered head. In some other animals, higher currents evoked strong, contralateral head-turning.
These two behaviors are unique to aPAG animals, as no animals in the mPAG group displayed backing, hunching, or head-turning during stimulation. Comparison of e and f shows the difference between the spectrograms with the two types of stimulation.
In the mPAG rat, both left and right HPC had a dominating power band toward the bottom end of the 5-10 Hz range before stimulation.
The low frequency is consistent with the rat being restricted in an operant chamber, and so confined to slower and smaller movements.
During mPAG stimulation this band is still visible, but power is greatly reduced. Once stimulation is turned off, 5-10 Hz activity returns, with a rebound effect, in the HPC, of showing stronger power than during pre-stimulation. In contrast to HPC, for the left PrL, there is a clear power band below 5 Hz-although some modest oscillatory activity is also present at higher frequencies. During mPAG stimulation, as with HPC, PrL power is greatly reduced across all frequencies. Immediately after stimulation is turned off, theta power transiently recovers but, especially in the following seconds, is generally below the prestimulation level.
In the aPAG rat, pre-stimulation activity is somewhat different, with the rat shown having a much wider band of left HPC activity that is strong above 5 Hz and additional power increases at 15 Hz. aPAG stimulation produces a clear increase in 15 Hz HPC power with a modest reduction in the 5-10 Hz range in both hippocampi and no obvious rebound. With PrL recording, aPAG stimulation produces a modest increase in power in the region of 10 Hz, with a rebound below 5 Hz.
3.3.2 | Effects on local theta oscillations and nondirected connectivity Figure 6 shows both raw theta measures for pre-stimulation and during stimulation and the difference between the two, that is, the effect of the stimulation. There were differences between aPAG and mPAG rats The first of each pair of panels shows the raw data and the second shows the dur-pre difference, that is, the effect of stimulation (stim. effect). The observed effects and differences are detailed in the text. Statistical coding as in Figure 4. PAG, periaqueductal gray; PPC, pairwise phase consistency during the pre-stimulation periods (i.e., immediately before stimulation). For simplicity, and comparison with Section 3.2, the text below focuses on the difference scores and then notes the prestimulation differences.

| DISCUSSION
Stimulation of the mPAG that induced freezing without escape reduced the power of both HPC and PrL theta by $30%. The reductions ceased as soon as stimulation was terminated, with a rebound in HPC theta power to a level above the pre-stimulation period. mPAG stimulation reduced PrL theta frequency to below 6 Hz (close to the HPC minimum for free moving animals) but did not affect HPC theta frequency.
In contrast to mPAG, RPO stimulation tended to produce a slight increase in HPC (but not PrL) power; and produced a small but significant reduction in frequency (from $7.7 to $7.3 Hz) in both HPC and PrL. While small, this drop in frequency contrasts with the usual increase (Green & Arduini, 1954;McNaughton & Sedgwick, 1978;Vertes, 1981). In free-moving but stationary rats, 100 Hz 0.1 ms stimulation produces HPC theta of 7-8 Hz with 50 μA pulses (McNaughton & Sedgwick, 1978, fig. 1); our pre-test adjustments of RPO stimulation also used $50 μA, and generated 7-8 Hz HPC theta against a background of no theta. However, during our comparison of mPAG and RPO, HPC theta was already at 7.5 Hz before the onset of stimulation. We did not explicitly monitor movement but, even without movement, our repeated activation of the mPAG (see also discussion of aPAG) at the threshold of freezing could have conditioned an emotional state that resulted in high basal theta activity in the hippocampus (Mikulovic et al., 2018;Sainsbury, Heynen, & Montoya, 1987 connections. This is consistent with the finding that mPAG stimulation reduces PrL but not HPC theta frequency. Our aPAG stimulation was similar to mPAG stimulation in generating freezing without escape. As with mPAG, aPAG stimulation did not produce any clear changes in coherence, PPC, or phase-directed coherence. However, the aPAG group results differed from mPAG in other respects as detailed below and in ways that are consistent with an antagonistic relationship between panic (and perhaps fear in general) with anxiety (Deakin & Graeff, 1991).
First, it is likely that aPAG stimulation produced a different form of freezing to mPAG, with more high-rate breathing, strong sniffing, and occasional facial contractions during stimulation. Critically, high levels of aPAG stimulation continued to elicit freezing without generating escape.
Second, aPAG had different pre-stimulation recordings to mPAG.
aPAG rats had higher power at higher frequencies (10 Hz and above) compared to mPAG rats ( Figure 3). Also, in the 5-10 Hz range analyzed, the aPAG group had higher baseline HPC power (but somewhat less PrL power) and a tendency to a lower HPC frequency. Since aPAG and mPAG rats were only separated after histology, it seems likely (as argued above) that the presence of substantial theta in the absence of stimulation was the result of conditioning (or some after effect of PAG stimulation) with the state generated by aPAG stimulation being different from mPAG and so resulting in different after-effects. Note that RPO stimulation does not generate such post-stimulus activation; and, with standard RPO testing, there is a tendency for the rats to become somnolent with long periods of testing and so produce no background theta.
Third, aPAG stimulation (like RPO stimulation) caused increases in power while mPAG caused decreases. The increases in power with aPAG stimulation were greatest for PrL, less for R-HPC, and negligible for L-HPC. The PrL change reflected a shift to equality with HPC from its before-stimulation lower power.
Fourth, aPAG stimulation (like RPO stimulation) did not produce the clear reduction in prefrontal frequency seen with mPAG. However, as with power, there appears to be a trend to equalization of PrL with HPC compared to before stimulation.
Theta therefore appears to control a process that is important for sustaining clinical anxiety but, at least when elicited by RPO stimulation, does not represent an anxious state as such.
aPAG stimulation in the current experiments had much the same electrophysiological effects as RPO stimulation but, in addition, produced a form of freezing coupled with a range of other reactions consistent with an elicited state of anxiety. The data also suggest that this elicited anxiety became conditioned to the apparatus to at least some extent: producing theta that would not otherwise have been observed (e.g., as in the RPO-alone stimulation we used in preliminary assessment). mPAG stimulation, by contrast, seemed to produce a somewhat different form of freezing and, at higher intensities, undirected escape reactions-consistent with the conventional interpretation that it elicits panic/fear (Deakin & Graeff, 1991;Evans et al., 2018). Contextual fear conditioning is normally seen as generating anxiety (i.e., anxiolytic-sensitive reactions) rather than fear (Le on et al., 2017; Luyten et al., 2011) and so it is not surprising that, like aPAG, mPAG stimulation appears to condition a state with high background theta.
However, the mPAG unstimulated theta has less variation between the recoding sites (with PrL having lower power and higher frequency in the mPAG case than aPAG). These differences could be because the unconditioned states are qualitatively different or simply because they produce different amounts of contextual conditioning. Critically, whatever the specific causal chains involved, the data are consistent with anxiety (or anxiety processing states) being linked to theta and suppression of such anxiety-related processes by panic/fear in general and mPAG stimulation in particular, including reports from humans (see Carrive & Morgan, 2012). Our data are also consistent with theories of the parallel hierarchical neural organization of fear and anxiety systems extending from the prefrontal cortex down to the PAG (Silva & McNaughton, 2019).