Introduction

Cognitive control is the ability to judiciously use task-relevant information while ignoring salient concurrent information that is currently irrelevant for the task. The consensus view is that the prefrontal cortex (PFC) is crucial for cognitive control such that PFC activity is central to the underlying computations that operate through prefrontal interactions with diverse brain areas during cognitive control tasks [1–3]. The hippocampal role in the cognitive control of spatial and mnemonic information has been investigated using an active place avoidance task [4]. The task conditions rodents to avoid the location of a mild shock and uses continuous rotation of the behavioral arena to dissociate the environment into two spatial frames, one defined by the rotating arena and the other defined by the stationary room [5]. Avoiding shock at a stationary room location requires using the task-relevant spatial information and ignoring spatial information from the other, task-irrelevant frame. Dorsal hippocampus dysfunction impairs this avoidance [6–9]. During active place avoidance task variants, place cell discharge in hippocampus subfields, and head-direction cell discharge in the medial entorhinal cortical input demonstrate cognitive control of spatial representations; neural ensemble activity alternates between representing room and arena locations depending on the subject’s proximity to the frame-specific location of shock [4, 5, 10–12].

A standard “central-computation” hypothesis centralizes the computations that are necessary for cognitive control to the PFC and predicts that lesion of the PFC will impair active place avoidance [13, 14]. Recent work indicating that the hippocampus [10] and sensory thalamus [15] are also necessary for cognitive control can be reconciled with the central-computation hypothesis by assuming there are necessary interactions between modality-specific information processing in hippocampus or thalamus and the primary, control-specialized processing in the PFC [14, 16, 17]. Alternatively, the role of non-PFC brain structures in cognitive control could be explained by a “local-computation” hypothesis; the computations needed for cognitive control can be performed locally in neural networks specialized for the particular information upon which the task depends and on which cognitive control operates. The local-computation hypothesis predicts that PFC lesion can spare cognitive control. We tested these mutually exclusive hypotheses by making medial prefrontal cortex (mPFC) lesions in rats and evaluating their behavior and hippocampal physiology during active place avoidance (Fig. 1A).

Figure 1.

Results

mPFC lesion does not impair cognitive control in the active place avoidance task

Ibotenic acid, but not vehicle, targeted to mPFC caused lesions that included the cingulate cortex, prelimbic and infralimbic areas (Fig. 1B). To assess the effect of mPFC lesion on behavior, we first examined locomotion during the pretraining sessions (Fig. 1C). There is no effect of lesion on the total distance walked across the session (Fig.1D). We next examined active place avoidance learning during initial training (Days 2-3). Both groups learn to reduce their entrances into the shock zone with training (Fig. 1E). There is no behavioral difference between groups by measuring the number of entrances or time to first enter the shock zone, the clearest estimate of memory (Fig. 1F). The two groups are also indistinguishable on Day 4 during which 24-h memory is assessed by their times to first enter the shock zone and the number of entrances (Fig. 1, Day 4 Retention).

We then trained the rats to avoid the shock zone relocated 180° to assess the impact of mPFC lesion on cognitive flexibility with the additional challenge to distinguish between the current and previously learned location of shock, under the cognitive control challenge [18]. Both groups learn the conflict task variant, measured by the number of entrances and the time to enter the new shock zone (Fig. 1, Days 4-5 Conflict Training). There is a non-significant trend for greater savings in sham rats, comparing the entrances on the first and second days of conflict trials (sham = 11.12 ± 4.58 entrances, lesion = 4.3 ± 4.09 entrances; t₁₆ = 1.11; p = 0.3, d = 0.52). Consequently, we more closely examined the initial experience of the new shock zone location when the rats are first confronted with the change in the location of shock. We compared the number of entrances on the first trial of each day to estimate savings during initial training (Fig. 1F left inset) and during conflict training (Fig. 1F middle inset). During initial training, both sham (t₇ = 4.5, p = 0.003, d = 1.2) and lesion (t₉ = 4.7, p = 0.001, d = 1.5) rats improved. During conflict training, while sham rats all improved across the first two conflict trials of each day (t₇ = 3.8, p = 0.006, d = 1.4) only 6/10 lesion rats improved (t₉ = 0.87, p = 0.4, d = 0.27). We also quantified the improvement for each subject during conflict as a savings index using the first 5 minutes of the first conflict trials of each day (Fig. 1F right inset). The sham rats improved, avoiding more than the lesion rats in the first 5 minutes of the second conflict day compared to the first conflict day (t₁₆ = 2.4, one-tailed p = 0.03, d = 0.85), despite equivalent performance across eight trials on the first day. Although, we find no evidence that mPFC lesion causes a deficit in active place avoidance learning to avoid the initial location of shock, nor do we observe a deficit in memory retention up to 24-h for the initial shock location, both of which require cognitive control, the lesion is sufficient to impair conflict learning after a 24-h delay.

mPFC lesion alters functional relationships amongst related brain areas

Cytochrome oxidase (CO), a sensitive metabolic marker for neuronal function [19], was used to evaluate whether lesion effects were restricted to the mPFC. In a subset of rats (lesion n = 8, sham n = 8), CO activity was evaluated in 14 functionally-related brain regions (Fig. 2A). CO activity in the central amygdala, but not elsewhere, is altered by the lesion, with activity increased in lesion rats (Table S1).

Figure 2.

Pair-wise interregional correlations of CO activity were evaluated to assess whether extra-mPFC functional relationships are altered after mPFC lesion. After FDR correction (0.01) 15 of the 91 correlations are significant in the sham sample but only 7 in the lesion sample (Fig. 2B, S2b; test of proportions: z = 2.26, p = 0.03). The correlations within the ventral hippocampus are preserved after mPFC lesion, but some correlations are lost within the dorsal hippocampal formation (Fig. 2B). Correlations between the nucleus reuniens and both the dorsal and ventral hippocampus are also lost after mPFC lesion. Correlations between the basolateral amygdala and nucleus reuniens and between the basomedial amygdala and ventral hippocampus are lost after mPFC lesion, and a correlation between dorsal CA1 and the central nucleus of the amygdala appears after mPFC lesion. These data indicate that the mPFC lesion causes functional changes beyond mPFC.

mPFC lesion does not alter cognitive control of hippocampal neural representations

Cognitive control requires the subject to process and use a class of information purposefully, at the expense of other competing information, which we previously demonstrated in the population dynamics of dorsal hippocampus spatial discharge as rats and mice navigate on a rotating arena [4, 5, 10, 20]. Thus, we examined spatial frame ensemble preference (SFEP), the representational dynamics of hippocampus CA1 discharge on the rotating arena during pretraining, as well as retention of the conditioned place avoidance of the initial and conflicting shock zone locations (Fig. 3A). Note that the physical conditions were identical in these three trials because the shock was off.

Figure 3.

We first examined basic discharge properties of individual CA1 principal cells, which do not differ between the sham and mPFC lesion groups (Fig 3B). We then examined overdispersion of place cell firing, which is known to be reduced by prefrontal inactivation [21]. In contrast, overdispersion is expected to increase with cognitive control and other extra-positional processes that increase discharge non-stationarity [22]. We find that during pretraining, overdispersion of hippocampal discharge is reduced by a factor of two in lesion compared to sham rats, consistent with an effective lesion (Fig. 3C left). Might this also indicate reduced cognitive control in the lesion rats? Remarkably, during conditioned place avoidance the overdispersion in lesion rats increases to the level observed in sham rats, consistent with the increased demand for cognitive control (Fig. 3C middle and right). These findings confirm the effectiveness of the prefrontal lesion and also provide electrophysiological evidence consistent with intact cognitive control in the lesion rats.

We next examined the spatial frame ensemble preference (SFEP), the electrophysiological signature of cognitive control, during conditioned place avoidance. CA1 ensemble discharge alternates between preferentially representing room locations and arena locations (Fig. 4A; runs tests for all recordings were significant; z’s ≥ 12.3, p’s ≤ 10^-33). During pretraining, SFEP is equally likely to signal the current room location or alternatively, the current arena location, in both the sham and lesion rats (Fig. 4A). Importantly, during avoidance of the initial shock location, SFEP favors the arena frame because the rat avoids the room frame location of shock (Fig. 4B) and SFEP is room-preferring near the shock zone and arena-preferring far from the shock zone (Fig. 4C,D), as previously demonstrated [4, 5]. Relocating the shock zone 180° for conflict training changes SFEP, but only in the lesion group such that it is equally probable to be room- or arena-preferring, consistent with needing to decide between avoiding the current and previous shock location throughout the environment (Fig. 4B) [23]. The ability to manipulate SFEP of hippocampal representational dynamics by the presence and location of shock directly demonstrates cognitive control at the level of the hippocampus. Despite an influence of mPFC lesion on cognitive control and the non-stationarity of hippocampal representations of space in hippocampal discharge, the lesion does not impair cognitive control in the place avoidance task.

Figure 4.

Discussion

Although permanent lesion of the mPFC alters baseline metabolic coupling of mPFC- and hippocampus-related brain areas (Fig.2, S2), as well as alters hippocampal representational dynamics during expression of cognitive flexibility in an active place avoidance task (Fig.3), we find no behavioral or electrophysiological evidence that the lesion impairs cognitive control in the basic active place avoidance task (Fig.1, 3, 4), despite the lesion being sufficient to cause a cognitive flexibility deficit when the shock zone was relocated for the conflict training trials (Fig. 1F insets). We also observed reduced bursting in the discharge of the underlying thalamus in the lesion group compared to sham (Table S2). Thus, CO imaging and electrophysiological evidence identify changes in the brain beyond the directly damaged mPFC area. In particular, the dorsal hippocampus loses the inhibitory input from mPFC [24, 25] and loses the metabolic correlation with the nucleus reuniens, which is thought to be a relay between the mPFC and the dorsal hippocampus [26, 27].

We have previously demonstrated cognitive control in the active place avoidance task variant we used (Fig. 1) because the rats must ignore local rotating place cues to avoid the stationary shock zone. Even when the arena does not rotate, rats distinctly learn to avoid the location of shock according to distal visual room cues and local olfactory arena cues, such that the distinct place memories can be independently manipulated using probe trials [28, 29]. When the arena rotates as in the present studies, neural manipulations that impair the place avoidance no longer impair when the irrelevant arena cues are hidden by shallow water [8, 9, 30, 31]. Furthermore, persistent hippocampal neural circuit changes caused by active place avoidance training are not detected when shallow water hides the irrelevant arena cues to reduce the cognitive control demand [4, 32, 33]. While these findings unequivocally demonstrate the salience of relevant stationary room cues to use for avoiding shock and irrelevant arena cues to ignore during active place avoidance, the most compelling evidence of cognitive control comes from recording hippocampal ensemble discharge. Hippocampal ensemble discharge purposefully represents current position using stationary room information when the subject is close to the stationary shock zone and alternatively represents rotating room information when the mouse is far from the stationary shock zone [Fig. 4; 4]. Despite this evidence from task design, behavioral observations, and direct electrophysiological demonstration of cognitive control, one might still argue that it is logically possible that the active place avoidance task does not require cognitive control and this is why the mPFC lesion was not impairing. While we reject such circular reasoning that only tasks that require mPFC can be cognitive control tasks, we nonetheless acknowledge that this may be the opinion of others.

We assert the evidence is compelling, and together these findings require rejecting the central-computation hypothesis that the mPFC is essential for the neural computations that are necessary for all cognitive control tasks. The present findings do not rule out a role for mPFC and hippocampal interactions in other cognitive processes [2, 34–37], especially those that rely on the ventral hippocampus and perhaps the direct excitatory mPFC-ventral hippocampus and direct mPFC-dorsal hippocampus inhibitory connections [24, 25]. Nonetheless, the present data favor the local-computation hypothesis because bilateral and unilateral inactivation as well as numerous other manipulations of the dorsal hippocampus impair learning, consolidation, and retrieval of the conditioned active place avoidance [6–8, 38], and optogenetic silencing of dentate gyrus impairs conflict learning. Indeed, we directly confirmed that mPFC lesion does not reduce the goal-directed representational switching of dorsal hippocampus (Fig.4) that is a sine qua non for cognitive control of spatial information [4, 5, 10, 11, 20, 23], despite detecting other behavioral and electrophysiological effects of the lesion. It is especially important that mPFC lesion impaired conflict avoidance of the relocated shock zone even after eight trials during which both sham and lesion animals learned and performed at a behavioral asymptote (Fig. 1F). We suggest this identifies a key role for the mPFC in cognitive flexibility, specifically the ability to judiciously select between different memories of the location of shock, and that this is demonstrably distinct from judiciously selecting between room- and arena-based classes of spatial information.

It is possible that acute manipulations of the mPFC, such as optogenetic, chemogenetic or pharmacological inactivation would lead to different results. In which case, post-lesion reorganization of neural circuits is sufficient to compensate for the loss of the mPFC functions that support cognitive control as required by the active place avoidance task. However, this would not change the conclusion to reject the central-computation hypothesis. It is also possible that the role of the mPFC is best evaluated in tasks that rely on egocentric spatial information [39], rather than the allocentric spatial information upon which the present active place avoidance task variant relies, but this possibility would also require rejecting, or at least severely limiting, the central-computation hypothesis. Finally, we observed impaired performance characterized by cognitive flexibility in the conflict task variant, demonstrating not only that the lesion was effective, but a dissociation of the neural substrates needed for different components of cognitive control, as has been previously reported [40, 41]. The mPFC has a crucial role in judiciously selecting between distinct memories of the room-frame shock location, which arguably is a form of memory-guided cognitive control that also depends on hippocampus [5, 23, 42], but not a form of cognitive control that has been demonstrated in the representational discharge of hippocampal cells [5]. Taken together, the findings are consistent with the local-computation hypothesis that diverse (but probably not all) brain circuits can perform the computations needed for cognitive control of the information and/or behavior for which that circuit is specialized to process.

Materials and methods

All methods complied with the Public Health and Service Policy on Humane Care and Use of Laboratory Animals and were approved by NYU’s University Animal Welfare Committee under protocol 12-1383, which follow National Institutes of Health guidelines.

Subjects

Thirty adult, male, Long-Evans rats were purchased from Charles River to arrive at New York University at approximately 40 days old. The rats were given at least one week to acclimate to the facility and were housed two per cage on a 12:12 light:dark cycle with free access to food and water. Fifteen rats were assigned to the lesion group and 14 to the sham group. One animal was used to confirm the coordinates, by injecting fluorogold (Fluorochrome, Denver, CO, 2.5% w/v in distilled water) at the same rate and for the same duration as for the lesions. Three of the sham and three of the mPFC lesions rats were used for the electrophysiology recordings. The behavior of these rats was not included in the behavioral assessment due to the requirement for extended training sessions to record neuronal activity.

Behavioral analysis followed by cytochrome oxidase activity

One rat in the lesion group was excluded because the lesion was inadequate and one had to be removed from the study due to a skin irritation, resulting in 10 rats in the lesion group. All of these rats underwent behavioral assessment (lesion n=10), eight of which were stained for cytochrome oxidase. The tissue from one lesion and one sham rat could not be used for CO because the tissues were lost to thawing and tissue from one rat from each group was not processed. For the control group, one rat was excluded from performing behavioral experiments because his lesioned cage mate had a skin irritation that prevented the lesion rat from behavioral testing; this control brain was still processed for cytochrome oxidase activity. One rat was excluded from the behavioral experiments due to equipment malfunction during training, and one rat had tissue damage during the sham surgery that appeared as a lesion, so this rat was excluded from the study. Thus, the sham group had eight rats for behavioral assessment and eight for cytochrome oxidase assessment. Our final group numbers for the behavioral assessment were therefore sham group n = 8, lesion group n = 10, and our final group numbers for cytochrome oxidase activity were sham group n = 8 and lesion group n = 8.

Lesion Surgery

Between the ages of 48-56 days, the rat received prefrontal cortex or sham lesion under sodium pentobarbital (50 mg/kg, i.p.) anesthesia. Bilateral lesion targeted two sites in each hemisphere at the following stereotaxic coordinates: from Bregma, Site 1: A.P +2.5, ML ±0.6, DV - 5.0, relative to the skull surface, Site 2: AP +3.5, ML ± 0.6, DV - 5.2, relative to skull surface. Ibotenic acid (0.06 M in PBS) was used to create an excitotoxic lesion. It was injected with a micro-infusion pump at a flow rate 0.2 μL/min through a stainless-steel cannula (0.25 mm outer diameter). Rats in the sham group underwent the same surgery procedure, but only PBS was injected. The rats were given one week to recover.

Active place avoidance training

After recovery from surgery, the rats were handled for 5 min per day for 5 days. On the first day of behavior training, animals had a pretraining session that consisted of two trials during which the rats were allowed to explore the stationary arena to habituate to the environment. The following two days, the rats underwent eight training trials per day on a rotating arena, learning to avoid entering a 60° sector shock zone in which the rats received a mild foot shock (500ms, 60Hz, 0.4mA). This shock zone was defined by the stationary cues in the room. On the fourth day of the behavioral training, the rats had a single trial with the shock on to test retention of the training. The rats were then trained on days 4 and 5 in a conflict training session in which the rats received eight trials per day with the shock zone located 180° from the initial training position. All trials were 10 min with 10 min intertrial intervals. The total distance walked during the session measured locomotion, the time to first enter the shock zone location measured avoidance memory, and the number of entrances into the shock zone was used to measure place learning. Savings in place learning on the conflict trials when the shock zone was relocated 180° was assessed for each rat by a savings index, which compares the number of entrances during the first 5 minutes of the first and second day of conflict training, Conflict Trials 1 and 9 respectively. The difference is normalized by the number of entrances on the first conflict trial and reported as a percentage:

Electrophysiology procedures

The following procedures have been previously described in detail[12]. Briefly, rats were anesthetized (pentobarbital 50 mg/kg) and received mPFC or sham lesions as described above. A Neuropixel electrode array was implanted in the brain to traverse AP −3.8 mm, ML 2.5 mm for recording dorsal hippocampus and the tip targeted DV 9.3 mm, which allowed recording from the underlying thalamus. The Neuropixel was stabilized and protected by mounting it in a custom 3-D printed appliance. The appliance and five bone screws were cemented to the skull (Unifast, GC America Inc., IL). The rats were allowed at least 1 week to recover before further manipulation.

The rats received active place avoidance training as described above. In addition, to permit neural ensemble recordings, after the last pretraining, training, and conflict training sessions, shock was turned off and the session was extended for 10 - 20 min during which CA1 neural ensemble discharge was recorded. Signals were filtered between 300 Hz and 10 kHz and sampled at 30 kHz for single unit recording. Because of these extended sessions, the behavior from these animals was not included in our behavioral assessment.

Electrophysiology data analysis

Single units were automatically sorted using Kilosort 2[43]. Units were only studied if the estimated contamination rate was < 20% with spikes from other neurons. This was estimated from refractory period violations relative to expected. Units with non-characteristic or noisy waveforms were also excluded. We computed a “burst ratio” to characterize burstiness of discharge as: number of spikes with ISI ≤ 30 ms divided by number of spikes with 100 ms ≥ ISI ≤ 130 ms.

Hippocampus single units were classified as complex-spike or theta cells as in prior work[10, 44, 45]. Complex-spike cells appear to be pyramidal cells, having long-duration waveforms (> 250 µs), low discharge rate (< 5 AP/s) and a tendency to fire in bursts. Theta cells are likely local interneurons[46], and had short-duration waveforms (< 250 µs), high discharge rate (> 2 AP/s), and did not tend to fire in bursts. Only hippocampal units classified as pyramidal cells were studied. Thalamic units fired at high rates and were less likely to discharge in bursts (Table S2).

The representational dynamics of CA1 discharge alternating between preferentially signaling the current location in the stationary room or the current location on the rotating arena was investigated to directly evaluate cognitive control during place avoidance behavior. Spatial-frame-specific momentary positional information (I_pos) is used to estimate the location-specific information from the activity of a cell during a brief interval (Δt = 133 ms).

is the probability of observing i activity at location x, and p_i is the overall probability of observing i activity. I_pos(room) estimates the information about the current location x in the room, whereas I_pos(arena) separately estimates the information about current location x on the rotating arena[47]. I_pos(t) can be positive or negative, but the absolute value is large whenever the cell’s activity at the current location is distinct or ‘‘surprising’’ compared to the location-independent probability of observing the same activity. The value |I_pos(t)| is abbreviated I_pos A unit’s spatial-frame preference is estimated at each Δt moment as the difference: ΔI_pos = I_pos(room) - I_pos(arena), where positive values indicate a momentary preference for signaling the room location, and negative values indicate a preference for signaling the arena location. ΔI_pos timeseries were evaluated for significance compared to chance fluctuations around the mean using the z value from a runs test. The spatial frame ensemble preference (SFEP) at each Δt was computed as the average ΔI_pos over all cells. This ensemble ΔI_pos timeseries (assessed for significance by runs test) was averaged over all times and locations to estimate an overall SFEP as the proportion of time the ensemble activity was room (or arena) preferring[4, 5]. To evaluate if SFEP was purposeful, the arena was divided into approximately two halves based on where individual rats avoided during a retention test with shock off. The “near” half was the largest sector that includes at least 50% of the recording time and surrounds the least visited 20° sector of the room frame (within the shock zone). The “far” half is the remaining sector. The probability of observing room-preferring ΔI_pos was computed at each location to make a representational preference map and the average probability in the near and far halves were compared by statistical evaluation.

Overdispersion of hippocampal discharge was computed as previously described[22]. Briefly, The spike and position timeseries is divided into 5-s episodes and the expected firing (exp) is computed assuming an inhomogeneous Poisson process. The standardized rate is computed as the normalized difference between the observed (obs) and expected firing: . Episodes with exp = 0 are undefined, and only episodes in which exp is greater than the cell’s mean firing were analyzed because this criterion selects episodes during which the rat passes through the central region of the firing field. The variance of the distribution of z values quantifies overdispersion, which can be compared between conditions using a F test.

Tissue processing and histochemistry

Verification of lesion coordinates

Immediately following surgery to inject fluorogold, the rat was transcardially perfused with 0.9% saline and 10% formalin, and the brain extracted and postfixed in 10% formalin for 24h at 4°C. The brain was then cryoprotected in 30% sucrose (w/v in 1x phosphate buffered saline) and stored at 4°C until cut on a cryostat (40 μm). The sections were mounted onto gelatin coated slides and scanned using an Olympus VS120 microscope (fluorescence, 10x). The slides were then placed in e-pure water (2 times, 1 min each) before being dehydrated (50, 70, 80, 95, 100, 100% ethanol, 1 min each) and defatted in 50:50 chloroform ethanol for 20 min. The slides were then rehydrated (100, 100, 90, 80, 70, 50% ethanol, 1 min each) before being Nissl stained with Cresyl violet (1 min), rinsed in e-pure water, and dehydrated (50% ethanol, 1 min, 70% ethanol until the white matter was white, 95% ethanol with acetic acid until desired intensity, 100, 100% ethanol, 2 min each). The tissues were then cleared in xylenes (3 times, 5 min each) before being coverslipped. The slides were then scanned again (light microscope, 10x). Images were manipulated using Adobe Photoshop CS6 to perform auto-contrast on each image and to remove background from the images. The fluorescent images were then all thresholded in a single operation to remove the background and the result was superimposed on the corresponding Nissl-stained sections using Adobe Illustrator.

Cytochrome oxidase activity and Nissl staining

On the day following the completion of behavior training, the rats were anesthetized with isoflurane, immediately decapitated, and the brains were extracted. The brains were rapidly frozen in isopentane on dry ice and stored at −80 °C. Sets of brains consisting of two to four animals per group were cut simultaneously on a cryostat (40 μm), and sorted into the three series, one of which was Nissl stained and one used for cytochrome oxidase histochemistry. The series used for Nissl staining was stored at room temperature and the other two series stored at −80°C until processed for cytochrome oxidase histochemistry. To control for variability across batches of histochemical staining, 20, 40, 60 and 80 μm sections of fresh rat brain tissue homogenate (prepared as in [48]) were included. Cytochrome oxidase staining was performed according to [49]. Stained slides were scanned with an Olympus VS 120 light microscope (2 x) and the optical densities measured from captured images. Images were converted to 8-bit gray scale using ImageJ [50] and optical densities read from the standard slides and 14 brain regions (Figure 1). These brain regions included the dysgranular and granular retrosplenial cortices (RSD and RSG, respectively), the nucleus reuniens (RE), the central nucleus of the amygdala (CEA), basomedial and basolateral amgydala (BMA, and BLA, respectively), the dorsal hippocampal CA1, CA2, CA3 and dentate gyrus areas (dCA1, dCA2, dCA3, dDG, respectively) the ventral hippocampal CA1, CA3 and dentate gyrus areas (vCA1,vCA3, vDG, respectively) and the dorsal subiculum (DS). The optical densities were measured using ImageJ (NIH) and cytochrome oxidase activity was normalized as in [51]. Optical densities were measured while blind to the group identity. Three to six optical density readings were taken for each brain region, from both hemispheres and averaged for each individual subject.

To confirm the lesion site, one series of sections was stained with Cresyl violet. The slides were placed in the e-pure water (2 times, 1 min each) and dehydrated in a series of ethanol baths (50, 70, 80, 90, 100, 100%, 1 min each) prior to clearing the fats in a 1:1 mixture of ethanol:chloroform for 20 min. The slides were then rehydrated (100, 100, 90, 80, 70, 50% ethanol, 1 min each), rinsed in e-pure water, and placed in the Cresyl violet for 10 min. The slides were again dehydrated and cleared in xylenes (3 times, 5 min each) before being coverslipped. Images were captured with an Olympus VS120 light microscope at 10x. Only rats with verified bilateral lesions were included in the data analysis.

Statistical analysis

Cytochrome oxidase activity

Group averages of the optical densities were calculated for each brain regions and expressed as mean ± SEM relative cytochrome oxidase activity/ μm of tissue. Interregional metabolic covariation was examined by calculating Pearson correlations between each brain region. The statistical significance of interregional cytochrome oxidase activity correlations was determined by a 0.01 False Discovery Rate. To determine if the interregional correlations were significantly different between groups, we transformed the r value of the correlation to Fisher’s z-scores. Significant correlations (p < 0.05) were used to generate graph theoretical networks using the Brain Connectivity Toolbox in MATLAB (https://sites.google.com/site/bctnet/).

Behavior and electrophysiology

Place avoidance task performance and SFEP electrophysiology measures were compared using multivariate analysis of variance. Repeated measure analyses were performed using the Greenhouse-Geisser correction. One-factor group comparisons were performed by t test. Runs tests were used to evaluate the likelihood of observing the ensemble ΔI_pos timeseries by chance. Statistical significance was set at 0.05 for all comparisons. Test statistics, degrees of freedom and effect sizes are provided throughout, in addition to exact p values.

Additional information

Author Contributions

EP, KCO, and AAF designed the research, EP, KCO, GG, DT, KN, ASA, and NR performed the research, EP, NK, SS-C, and KCO analyzed the data and AAF, KCO, and EP wrote the paper.

Inclusion and Ethics

This research was conducted with respect for the animal subjects, the individual researchers and their unique perspectives, and our collective effort to learn and generate individually useful knowledge.

Funding Disclosure

Supported by NIH grants R01NS105472, R01MH115304, and R01MH132204 to AAF.