Abstract
Brain states fluctuate between exploratory and consummatory phases of behavior. These state changes affect both internal computation and the organism’s responses to sensory inputs. Understanding neuronal mechanisms supporting exploratory and consummatory states and their switching requires experimental control of behavioral shifts and collecting sufficient amounts of brain data. To achieve this goal, we developed the ThermoMaze, which exploits the animal’s natural warmth-seeking homeostatic behavior. By decreasing the floor temperature and selectively heating unmarked areas, mice avoid the aversive state by exploring the maze and finding the warm spot. In its design, the ThermoMaze is analogous to the widely used water maze but without the inconvenience of a wet environment and, therefore, allows the collection of physiological data in many trials. We combined the ThermoMaze with electrophysiology recording, and report that spiking activity of hippocampal CA1 neurons during sharp-wave ripple events encode the position of the animal. Thus, place-specific firing is not confined to locomotion and associated theta oscillations but persist during waking immobility and sleep at the same location. The ThermoMaze will allow for detailed studies of brain correlates of immobility, preparatory-consummatory transitions and open new options for studying behavior-mediated temperature homeostasis.
Introduction
All behaviors can be considered as parts of a sequence of action-rest transition1. Brain states in vertebrates fall into dichotomous categories, and correspond roughly to what early behavioral research referred to as “preparative” (or “exploratory”) and “consummatory” (or “terminal”) classes2. In mammals, these two fundamental brain states can be readily identified by basic electrophysiological monitoring of various brain structures3. They are also referred to as voluntary and non-voluntary or conscious and non-conscious brain states3. Switching between these states is correlated with high and low release of subcortical neuromodulators4–9. Consummatory behaviors include feeding and drinking, resting and its extreme form, non-rapid eye movement (NREM) sleep. Preparatory and consummatory behaviors in the hippocampus are associated with theta oscillations and sharp wave ripples (SPW-Rs), respectively10.
Deciphering the physiological underpinnings of these categories and revealing the significance of brain state transitions for cognition requires sufficient sampling of the relevant brain states. This is usually achieved by extended repeated recordings or, when possible, recording large numbers of neurons simultaneously. Prolongation of explorative behavior can be readily achieved by placing the animal in novel environments, by food or water deprivation or introducing delays in choice behavior tasks11,12. Recently, the honeycomb maze paradigm was introduced to extend the observation periods of explorative deliberation13.
In contrast, the experimental control of consummatory classes of behavior is more difficult. Sleep provides an opportunity for long recordings. Comparison of sleep before and after learning is a standard paradigm to examine experience-induced brain plasticity14,15. Consummatory brain states associated with eating, drinking and sex change rapidly with satiety and requires prolonged periods of deprivation16–19. Controlling periods of awake immobility is most difficult20–22, mainly because forced immobilization of the animal is stressful23 and is accompanied by altered physiological states24.
Here we introduce the ThermoMaze, a behavioral paradigm that allows for the collection of large amounts of physiological data while the animal rests at distinct experimenter-controlled locations. In standard laboratory environments (20-24 °C)25, both housing and data collection take place below the thermoneutral zone of mice (26-34 °C)26–28. The ThermoMaze exploits the animal’s behavioral thermoregulation mechanisms29,30 and promotes thermotaxis (i.e., movement in response to environmental temperature)31. Searching for a warmer environment, social crowding and nest building are natural behavioral components of heat homeostasis31–33. The ThermoMaze allows the experimenter to guide small rodents to multiple positions in a two-dimensional environment. Decreasing the maze floor temperature induces heat seeking behavior and after finding a warm spot, the animal stays immobile at that spot for extended periods of time, allowing for recording large amounts of neurophysiological data in immobility-related brain states. We report on both behavioral control and hippocampal electrophysiological correlates of heat seeking activity to illustrate the versatile utility of the ThermoMaze.
Results
Design and Construction of the ThermoMaze
The ThermoMaze is designed to guide small rodents to warm spatial locations in a two-dimensional cold environment, consisting of a box (width, length, height: 20, 20, 40 cm, respectively) made from an acrylic plexiglass sheet (Fig. 1A, top). The floor of the maze is constructed from 25 Peltier elements (40 x 40 x 3.6 mm) that are attached to aluminum water cooling block heatsinks (40 x 40 x 12 mm, n = 25) with heat-conductive epoxy and are insulated from each other by wood epoxy (Fig. 1A, dashed inset). Each Peltier element is controlled by an electrically operated switch (relay) that opens and closes high-current circuits by receiving transistor-transistor-logic (TTL) signals from outside sources (Fig. 1B). Peltier elements can be heated individually up to 30 °C to provide a warm spot for the animal when other regions of the floor are under cooling (Fig. 1B, active heating of one Peltier element is shown). The ambient temperature of the maze is controlled by water circulated from the water tank through the water-cooling blocks. We set the floor temperature to either ∼25°C (room temperature) or to ∼10 °C (cooling, Fig. 1C and Suppl. Fig. 1), but a range of ambient temperatures (5-30 °C) could be employed. The water temperature is monitored by a K-type thermocouple placed inside the water tank (Fig. 1A bottom). The floor temperature of the ThermoMaze is monitored using a thermal camera (FLIR C5) providing continuous registration of real-time temperature changes (Fig. 1A).
Prior to the experiments, the thermal camera, which continuously measures the surface temperature of the floor of the ThermoMaze is calibrated by thermocouples placed directly on Peltier elements (Fig. 2). The accuracy of the FLIR 5C infrared camera is ±3 °C. With proper calibration and attention to emissivity (an object’s ability to emit rather than reflect infrared energy) the margin of error can be less than 1 °C34.
Mice seek out hidden warm spots in the ThermoMaze
To illustrate the novel advantages of the ThermoMaze on behavior and brain activity, we tested 11 mice (n = 3 male and 8 female mice) with silicon probe recordings from the hippocampus (Suppl. Table 1). One wall was marked by a prominent visual cue (black tape and blinking light-emitting diode; LED) to provide a distinct spatial cue in the box (Fig. 1A)35. On each experimental day, the mouse was placed in the ThermoMaze and allowed to explore it for 10 min at room temperature (“Pre-cooling” sub-session; Figure 3A). Next, the ThermoMaze temperature was decreased to around 14 °C for 80 min and four Peltier elements (“warm spots”; typically, in the corners) were sequentially and repeatedly turned on and heated up to 30 °C. One Peltier element was turned on for 5 minutes in a sequential order (1-2-3-4) and the sequence was repeated four times (“Cooling” sub-session; Figure 3B). The Cooling sub-session was divided into 5-minute “warm spot epochs” for analysis. The daily experimental session ended with a “Post-cooling” sub-session (free exploration at room temperature for 10 min). In addition, all mice were recorded in the home cage both before and after the experimental session (Fig. 3A). During Pre- and Post-cooling sub-sessions, the animal explored the maze relatively evenly with a moderate movement speed (Figure 3B-D), although thigmotaxis was the dominant pattern, with corners as highly preferred sites of both movement and immobility (Suppl. Fig. 2). The animals readily found the location of the warm spot after a few training sessions (median = 3). Changing the warm spot locations during Cooling induced exploration until the mouse found another warm spot and stayed on it for prolonged periods (Figure 3B and C, n = 17 session in 7 mice). Duration spent on the warm spot roughly followed a bimodal distribution with a median = 2.85 minutes (Suppl. Fig. 2A). Compared to Pre and Post sub-sessions, during Cooling, mice spent a smaller proportion of time in movement (Pre: 40 ± 19%, Post: 34 ± 16%, Cooling: 23 ± 12%, mean ± SD, defined as speed > 2.5 cm/s, n = 20 sessions from 7 mice; Figure 3D) and more time in immobility (Pre: 59 ± 19%, Post: 66 ± 16%, Cooling: 76.74 ± 12.41%, mean ± SD, defined as speed ≤ 2.5 cm/s; n = 20 sessions from 7 mice; Figure 3D). The mice spent most of the time in the corners of the ThermoMaze where heat was provided (Suppl. Fig. 2B), compared to Pre- and Post-cooling (Figure 3C). Once the heating of the Peltier element was turned off, the animal quickly left its location (median duration = 12.99 s, n = 20 sessions from 7 mice; Figure 3E) and searched for a new source of warmth. Mice increased their speed from 0 cm/s to 2.5 cm/s within 12.28 s after a warm spot was turned off (median, n = 20 sessions from 7 mice; Figure 3F) and found the new warm spot within 23.45 s (median, n = 20 sessions from 7 mice; Fig. 3G). In two additional male mice, we examined brain temperature changes during the Cooling sub-session by implanting a thermistor in the hippocampus (Suppl. Fig. 3A). In support of previous findings, we found brain state-dependent fluctuation of brain temperature (Suppl. Fig. 3B)36–38. However, cooling the environment per se did not correlate with brain temperature changes (Suppl. Fig. 3C-E), confirmation that brain temperature is strongly regulated and is largely independent of the ambient temperature38. The ThermoMaze provides an affordance for mice to select their environmental temperature through the activation of behavioral thermoregulation39.
One of the objectives in developing the ThermoMaze was to induce immobility at several locations repeatedly and for extended time periods. To confirm that this objective was achieved, we ran control sessions with the same duration as the Cooling sub-session but at room temperature (80 minutes; Suppl. Fig. 4). Under room temperature condition (3 sessions in 3 mice), mice first explored the ThermoMaze and settled in one of the corners for an extended period of time. Although mice spent a similar total amount of time immobile under both conditions, the spatial distribution of immobility durations was more uniform in the Cooling sub-session (Suppl. Fig. 4) because the ThermoMaze paradigm forced the animals to leave their chosen spot and move to the experimenter-designated locations, i.e., the new warm spots away from the corner (Suppl. Fig. 5).
Firing rate maps of hippocampal neurons in the ThermoMaze
Compared to spatial learning and memory paradigms such as the Morris water maze40, the ThermoMaze has a non-aqueous environment and thus allows for an easy setup of electrophysiological recording. We recorded neurons from the CA1 hippocampal region by multi-shank silicon probes and separated them into putative pyramidal cells and interneurons (Methods - Unit isolation and classification section). We separated behavioral states (movement or immobility) based on movement speed (speed ≧ 2.5 cm/s = movement and speed < 2.5 cm/s = immobility).
To construct spike count maps for comparing sub-sessions, the ThermoMaze was divided into 25 x 25 bins and the number of spikes emitted by a neuron in each bin was counted and normalized by the time the mouse spent in each spatial bin. The impact of cooling during movement (theta state) was compared by calculating the correlation coefficients between Pre and Post, Pre and Cooling, and Cooling and Post spike count maps (Suppl Fig. 6A). The correlation coefficients decreased significantly across all sub-sessions, with the largest change observed between Pre-cooling and Post-cooling spike count maps in the experimental mice (Suppl Fig. 6B). Thus, the Cooling sub-session in the ThermoMaze induced a moderate decorrelation of pyramidal cells’ rate maps. Such observation constrained our ability to decode spatial information from the spiking activity during SWP-Rs in the Cooling sub-session using firing rate maps constructed during Pre- and Post-cooling sub-sessions41, because the Bayesian approaches have an underlying assumption that the spatial representation (tuning functions, or rate maps) is temporally stable.
In principle, comparison of place maps during the first and last 10 min of a 100 min session at room temperature should serve as controls. However, at room temperature mice “designate” one of the corners as home base after a few minutes of exploration and stay in that corner for the rest of the session (Suppl. Fig. 4A). Thus, exploration of the maze at the end of the session was not available.
Place-selective neuronal firing during SPW-Rs at experimenter-designated locations
As expected, SPW-Rs occurred predominantly in the corners (Fig. 4A), where the mice spent most of their time resting (Fig. 3C). Compared to room-temperature control sessions where animals spent most of their time in one corner, the spatial distribution of SPW-Rs in the Cooling sub-session was more uniform (Suppl. Fig. 4A-D), indicating that the ThermoMaze paradigm successfully biased where SPW-Rs were generated. The duration and amplitude of SPW-Rs were comparable in the ThermoMaze and the homecage (Fig. 4 B, C), whereas the mean peak frequency of SPW-Rs were significantly lower (Fig. 4D). This decrease can be explained by the lower brain temperate during sleep, a state in which the animals spent most of their time in the home cage36.
To quantify spatial tuning features of neuronal firing during SPW-Rs in the ThermoMaze during the Cooling subsession, we defined a metric referred to as “spatial tuning score” (STS). We first binned the floor of the ThermoMaze into four quadrants (2x2). For each neuron, we calculated its average firing rate within SPW-Rs in each quadrant. STS was then defined by the firing rate in the quadrant with the highest within-SPW-R firing rate divided by the sum of the within-SPW-R firing rates in all four quadrants (yielding a value between 0 and 1; Fig. 5A). To test the significance of STS, we compared the STS values with their shuffled versions by randomly assigning one of the four quadrants to each SPW-R. The distribution of the STS in actual SPW-Rs was significantly higher compared to shuffled controls (Fig. 5B). Additionally, pyramidal cells exhibited higher STSs compared to interneurons (medians: pyramidal cells = 0.3432; interneurons = 0.2934; one-sided Wilcoxon rank sum test, p < 0.001). In summary, both excitatory and inhibitory neuronal populations exhibit place-selective firing during SPW-Rs, while the excitatory neurons demonstrate a stronger place-specific firing.
To quantify how well CA1 neurons encode spatial information during SPW-Rs at the population level, we carried out a Bayesian decoding analysis to read out the current position of the animal from spiking activity41. We constructed firing rate map templates using spikes within SWP-Rs in the training dataset and determined animal positions that maximized the likelihood of observing the spike train during SWP-Rs in the testing dataset (see Method). Spiking activity during SPW-Rs reliably identified the quadrant that the animal was in above chance level (Figure 5C, D) irrespective whether we incorporated the spatial distribution priors into the decoder in an example session (Figure 5C) or used a uniform prior (Figure 5D).
To relate spatial content of spikes during SPW-Rs and locomotion, we examined whether the same or different groups of neurons contributed to the place-specific firing during SPW-R and locomotion by calculating the firing rate ratios within preferred quadrant versus all quadrants. These ratios during SPW-Rs and movement were positively correlated (Fig. 5E; n = 1150 pyramidal cells in 20 sessions from 7 mice), suggesting that place cells42 during movement preserved their spatial properties during SPW-Rs (see also Suppl. Fig. 7 for further analysis and findings on interneurons).
Finally, we tested whether the preservation of spatial features of neuronal spiking also holds at the population level by constructing population vectors separately during movement and SPW-Rs. We then computed the pairwise correlation coefficients between these two conditions. As was the case for individual pyramidal cells, population vectors for the same quadrant during movement were similar to those during SPW-Rs (Figure 5F). Overall, these findings support and extend the observation that spiking activity during SPW-Rs continue to be influenced by the animal’s current position43.
To test specifically whether perceptual sensing of environmental features is critical in position-specific firing of neurons during SPW-Rs, we prolonged the duration of warm spots. After the Pre-cooling sub-session, the ThermoMaze temperature was decreased to 16 °C for 80 min and two Peltier elements were heated in an alternate fashion to 30 °C for 20 min (Figure 6A). As expected, mice spent most of the time immobile on the warm spots (Figure 6A,B). Similar to the 5-minute protocol (Fig. 4A), SPW-Rs occurred predominantly on the warm spots (Fig. 6C). The increased duration of stay on the warm spot facilitated the occurrence of sleep, as quantified by our brain state scoring algorithm (Fig. 6D, SPW-Rs). REM sleep was not detected since REM state typically emerges after 20-30 min of NREM episodes44. Mice spent a higher fraction of their time in sleep during the 20 min, compared to the 5 min sub-session (p = 0.003, n = 19 sessions in 7 mice and n = 7 sessions in 4 mice, Suppl. Table 1). The average inter-NREM interval was 1000 seconds (Fig. 6F, n = 7 sessions in 4 mice). Comparing the spike content of SPW-Rs during awake immobility and NREM sleep, we found that Pearson’s correlation coefficients between population vectors constructed during waking movement and waking SPW-Rs were higher than between movement and NREM SPW-Rs (Fig. 6G). These findings further support the view that sensory inputs during waking SPW-Rs can affect spiking content of SPW-Rs.
Discussion
To investigate the importance of brain state transitions in a controlled manner, we developed the ThermoMaze, a behavioral paradigm that allows for the collection of large amounts of physiological data while the animal rests at distinct experimenter-controlled locations. Since the paradigm exploits natural behavior, no training or handling is necessary. We demonstrate that mice regularly explore a cold environment until a warm spot is identified. They spend most of the time on a warm spot and even fall asleep, thus exhibiting a high degree of comfort. We exploited the long immobility epochs following exploration and showed how neurons active during hippocampal sharp wave ripples (SPW-R) replay waking experience. The ThermoMaze will allow for detailed studies of brain correlates of preparatory-consummatory transitions and open new options for studying temperature homeostasis.
Warmth-seeking homeostatic behavior
There is a renewed interest in exploiting natural learning patterns, as opposed to training animals for performing complex arbitrary signal-action associations45–52. In poikilotherm animals (species whose internal temperature varies with environmental temperature), energy homeostasis is one of the most fundamental homeostatic processes. Heat homeostasis involves multiple levels of coordination from cellular to systems, from peripheral to central53,54. To maintain core body temperature, thermogenic tissues rapidly increase glucose utilization by brown adipose tissue and shivering by skeletal muscle55,56. The hypothalamic preoptic area (POA) is regarded as the most important thermoregulatory “center” in the brain57,58. Connecting this area of research to learning, the POA is bidirectionally connected with the limbic system and multiple cortical areas which assist both online maintenance of body temperature and preparing the body for future expected changes (“allostasis”)23,59,60. These allostatic mechanisms induce exploratory behavior, searching for a warmer environment61,62. A location that provides a warm shelter needs to be remembered and generalized for future strategies. Our paradigm offers means to investigate exploratory-consummatory transitions, wake-sleep continuity in the same physical location and, in the reverse direction, the physiological processes that evaluate discomfort levels, motivate behavioral transition from rest to exploration and the circuit mechanisms that give rise to overt behaviors.
Mice, and rodents in general, are acrophobic and agoraphobic and tend to avoid open areas. Instead, they tend to move close to the wall and spend most of their non-exploration time in corners63. Thus, while we were able to train mice to seek out and stay on warm spots in the center of the maze after extensive training, their evolutionary “counter-preparedness”47 to stay in predator-prone open areas competed with the reward of warming. While these trained mice did stay transiently on the central warm spot, they spent more time returning to the corners. Our mice were on a normal day-light schedule thus their training during the day coincided with their sleep cycle. This explains why after 5-10 min spent on the safe and temperature-comfortable corner warm spots they regularly fell asleep. Yet, we noticed that mice did not simply transition from walking to immobility but, instead, even after finding the warm spot they regularly and repeatedly explored the rest of the maze before returning to the newly identified home base. By changing the temperature difference between the environment and the warm spot, it will be possible to generate psychophysical curves to quantify the competition between homeostatic and exploratory drives in future experiments. These measures, in turn, could be used to study the impact of perturbing peripheral and central energy-regulating mechanisms.
For several applications, it is not needed to tile the entire floor of the maze with Peltier elements. For example, a radial-arm maze with cooled floors or placed in a cold box can be equipped with heating Peltier elements at the ends of maze arms and center, allowing the experimenter to induce ambulation in the 1-dimensional arms, followed by extended immobility and sleep at designated areas. In a way, the ThermoMaze is analogous to the water maze40, also an avoidance task, but many more trials can be achieved in a single session and without the inconvenience of a wet environment.
SPW-R spiking content biased by current position of animal depending on brain state
We demonstrate the utility of the ThermoMaze for addressing long-standing questions in hippocampal physiology. Preparatory and consummatory behaviors in the hippocampus are associated with theta oscillations and SPW-Rs, respectively10. SPW-Rs also occur during NREM sleep but studying the differences between waking and sleep SPW-Rs has been hampered by the paucity of SPW-Rs in typical learning paradigms21,22,64–67. Neural activity during SPW-Rs has been shown to replay activity patterns observed during previous spatial navigation experiences21,43,64 and can even be predictive of activity during future experiences68–70. However, the extent to which SWP-R spiking context is biased by the current position of the animal is less known, as systematic control of position during rest/sleep has posed difficulty. The ThermoMaze enables the experimenter to control the animal’s position during SWP-R states. In agreement with previous studies43,65,67, we found that neurons whose place fields overlapped with the quadrant of the maze had a higher participation probability in SPW-Rs occurring at that location compared to other neurons. This observation supports the notion that waking replay events can be biased by perceiving features of the surrounding environment43. However, when the mouse fell asleep at the same location this relationship was weakened but did not disappear. Another potential explanation for the decreased correlation between sleep SPW-R and waking exploration is deterioration of replay as the function of time15. Alternatively, the persisting significant correlation between sleep SPW-Rs and previous exploration may also indicate that factors other than the perception of the animal’s vicinity is responsible for sleep replay70–72. Continuity of waking experience replay in waking and sleep SPW-Rs have been hypothesized previously but not yet tested73. Using the ThermoMaze, this and other related questions can now be addressed quantitatively.
Acknowledgements
We thank Daniel Levenstein for useful comments on the manuscript. We thank Yiyao Zhang, Anna Maslarova and Leeor Alon for their help with different aspects related to the experiments. Supported by MH122391, and U19NS107616.
Experimental methods
Animals and surgery
All experiments were approved by the Institutional Animal Care and Use Committee at New York University Langone Medical Center. Animals were handled daily and accommodated to the experimenter and the ThermoMaze before the surgery and electrophysiological recordings. Mice (adult female n = 8, 22 g and male n = 5, 26 g) were kept in a vivarium on a 12-hour light/dark cycle and were housed two per cage before surgery and individually after it. Atropine (0.05 mg/kg, s.c.) was administered after isoflurane anesthesia induction to reduce saliva production. The body temperature was monitored and kept constant at 36–37 °C with a DC temperature controller (TCAT-LV; Physitemp, Clifton, NJ). Stages of anesthesia were maintained by confirming the lack of a nociceptive reflex. The skin of the head was shaved, and the surface of the skull was cleaned by hydrogen peroxide (2%). A custom 3D-printed baseplate74 (Form2 printer, FormLabs, Sommerville, MA) was attached to the skull using C&B Metabond dental cement (Parkell, Edgewood, NY). The location of the craniotomy was marked and a stainless-steel ground screw was placed above the cerebellum. Silicon probe (Suppl. Table 1) attached to a metal microdrive75 was implanted into the dorsal CA1 of the hippocampus (2 mm posterior from Bregma and 1.5 mm lateral to midline) and a copper mesh protective cap was built around the probe. Animals received ketoprofen (5.2 mg/kg, s.c.) at the end of the surgery and on the following two days. Each animal recovered at least 5 days prior to experiments. The electrophysiology data was digitized at 20000 samples/s using an RHD2000 recording system (Intan technologies, Los Angeles, CA). The number of recorded sessions from each animal is summarized in Supplementary Table 1.
Construction of ThermoMaze
The ThermoMaze is a box (width, length, height: 20, 20, 40 cm, respectively), made from acrylic plexiglass sheet (8505K743, McMaster, Elmhurst, IL). The floor of the maze was constructed from 25 Peltier elements (40, 40, 3.6 mm, Model: TEC1-12706, voltage: 12V, Umax (V): 15V, Imax (A): 5.8A, ΔTmax(Qc=0): up to 65 °C). Each Peltier element was glued inside a custom 3D-printed frame (file can be downloaded from https://github.com/misiVoroslakos/3D_printed_designs/tree/main/ThermoMaze) using dental cement (Unifast LC, GC America, Alsip, IL) and wood epoxy (Quick-Cure, product number: BSI201, Bob Smith Industries, Atascadero, CA). Once Peltier elements were secured in the 3D-printed frame, an aluminum water cooling block heatsink (40, 40, 12 mm; a19112500ux0198, Amazon.com) was attached to each Peltier element using heat-conductive epoxy (8349TFM, MG Chemicals, Ontario, Canada). A variable voltage source (E36102A Power Supply, Keysight Technologies, Santa Rosa, CA) was attached to four Peltier elements using a relay system (4-Channel Relay Module, product number: 101-70-101, SainSmart, Lenexa, KS). The relays were controlled by an Arduino Mega (Arduino Mega 2560 Rev3) running a custom written code. Five aluminum water cooling block heatsinks were connected together using silicon tubes (5/16” ID x 7/16” OD, product number: 5233K59, McMaster, Elmhurst, IL). One of the five heatsinks was connected to a mini submersible electric brushless water pump (240L/H, 3.6W, Ledgle, ASIN: B085NQ5VVJ) using silicon tubes and another one was routed to the water tank. We used 5 water pumps to circulate water through the 25 cooling blocks. The water pumps were placed inside a water tank (40, 40, 60 cm acrylic box) and were powered using a DC power supply (E3620A, Keysight Technologies, Santa Rosa, CA). The temperature of the water tank was monitored by a K-type thermocouple (5SC-TT-K-40-72, Omega, Norwalk, CT) attached to a handheld thermometer (HH800, Omega, Norwalk, CT) and recorded by a K-type thermocouple (5SC-TT-K-40-72, Omega, Norwalk, CT) attached to an AD595 interface chip (1528-1407-ND, Digi-Key, Thief River Falls, MN) connected to an analog input of the RHD2000 USB Eval system (Intan Technologies, Los Angeles, CA). To monitor the floor temperature of the ThermoMaze, a thermal camera (C5, Flir, Thousand Oaks, CA) was used.
Behavior
The ThermoMaze setup provides a customized temperature landscape, which the animal can freely explore and choose where to settle. Without any training or shaping, a mouse will search and find the unmarked warm spot and stay on it for extended periods due to thermotaxis (movement towards locations with preferred temperature around 26–29°C; Figure 3)29,30. When the heating Peltier element is turned off, the animal quickly leaves the spot and explores the maze again until it finds another warm spot.
On each experimental day, the mouse is taken from the animal facility during their light cycle. The animal is first recorded in its homecage for 1-2 hours (pre-home). It is then transferred into the ThermoMaze under room temperature to freely explore for 10 minutes (Pre-cooling). During the Pre-cooling sub-session, the water circulation system is circulating room temperature water and the Peltier elements are not activated. After the Pre-cooling sub-session, 4 kg of ice and two ice packs (25201, Igloo) are added into the water tank while the animal remains in the ThermoMaze. Within 1 minute, the temperature of the water in the tank stabilizes at 10-13 ℃. We then turn on the pump to cool down the ThermoMaze setup (it takes ∼120 seconds to cool down the floor to 10-13 ℃). At the same time, the Arduino-controlled Peltier element heating system is turned on to heat one of the four 4 x 4 cm2 for 5 minutes, followed by another Peltier device in a fixed sequence (Fig. 2). Such sequence is repeated four times (total of 80 min) during a Cooling sub-session. After the subsession, the animal explores again at room temperature for 10 minutes (Post-cooling sub-session). To increase the temperature back to ∼20 ℃, the ice packs are removed, and 6.5 L of 55 ℃ water is added into the tank. The temperature in the ThermoMaze returns to room temperature within 2 minutes. After the Post sub-session, recording of electrophysiological activity continues in the homecage for an additional 1–2 hours (post-homecage; Fig. 3A).
To quantify the behavior of the animal within the ThermoMaze, video is recorded using a Basler camera (a2A2590-60ucBAS Basler ACE2) using the mp4 format with a framerate of 25 Hz. TTL pulses are sent from the camera to the Intan recording system to synchronize the video and the electrophysiological recordings. The animal’s location is detected within a 25x25 cm region of interest (ROI), using a custom trained DeepLabCut neural network76. Detections with a likelihood below 0.5 are discarded. The occasionally missing trajectory detections are filled using MATLAB function “fillmissing” with method “pchip” which is a shape-preserving piecewise cubic spline interpolation and are then smoothed using a 7th-order one-dimensional median filter “medfilt1”. The detection quality is visually examined by superimposing the detected animal location in each frame on the video.
Brain temperature measurement
To examine the effects of changing environmental temperature on brain temperature homeostasis, we implanted one male and one female wild type mice (C57Bl6, 28 g) with a thermistor (Semitec, 223Fu3122-07U015) in the hippocampus (2 mm posterior from bregma and 1.5 mm lateral to midline)36. After 5 days of postsurgical recovery, the animal was placed inside the ThermoMaze and brain temperature and behavior were monitored (n = 5 sessions, each session consisted of pre-homecage, Pre, Cooling, Post and post-homecage epochs).
Quantification and statistical analysis
SPW-R detection and properties
SPW-Rs were detected as described previously from manually selected channels located in the center of the CA1 pyramidal layer (https://github.com/buzsakilab/buzcode/blob/master/detectors/detectEvents/bz_FindRipples. m). Broadband LFP was bandpass-filtered between 130 and 200 Hz using a third-order Chebyshev filter, and the normalized squared signal was calculated. SPW-R peaks were detected by thresholding the normalized squared signal at 5×SDs above the mean, and the surrounding SPW-R begin, and end times were identified as crossings of 2×SDs around this peak. SPW-R duration limits were set to be between 20 and 200 ms. An exclusion criterion was provided by manually designating a ‘noise’ channel (no detectable SPW-Rs in the LFP), and events detected on this channel were interpreted as false positives (e.g., EMG artifacts). The ripple detection quality was visually examined by superimposing the detected timestamps on the raw LFP traces in NeuroScope2 software suite77.
Sleep state scoring
Brain state scoring was performed as described in the study by Watson et al.,44. In short, spectrograms were constructed with a 1-s sliding 10-s window fast Fourier transform of 1,250 Hz data at log-spaced frequencies between 1 Hz and 100 Hz. Three types of signals were used to score states: broadband LFP, narrowband high frequency LFP and electromyogram (EMG) calculated from the LFP. For broadband LFP signal, principal component analysis was applied to the Z-transformed (1–100 Hz) spectrogram. The first principal component in all cases was based on power in the low (32 Hz) frequencies. Dominance was taken to be the ratio of the power at 5–10 Hz and 2–16 Hz from the spectrogram. All states were inspected and curated manually, and corrections were made when discrepancies between automated scoring and user assessment occurred.
Unit isolation and classification
A concatenated signal file was prepared by merging all recordings from a single animal from a single day. Putative single units were first sorted using Kilosort78 and then manually curated using Phy (https://phy-contrib.readthedocs.io/). After extracting timestamps of each putative single unit activity, the spatial tuning properties, identification of 2D place cells and place fields, and participation in SPW-Rs events were analyzed using customized MATLAB (Mathworks, Natick, MA) scripts.
In the processing pipeline, cells were classified into three putative cell types: narrow interneurons, wide interneurons, and pyramidal cells. Interneurons were selected by 2 separate criteria; narrow interneurons were assigned if the waveform trough-to-peak latency was less than 0.425 ms. Wide interneuron was assigned if the waveform trough-to-peak latency was more than 0.425 ms and the rise time of the autocorrelation histogram was more than 6 ms. The remaining cells were assigned as pyramidal cells77. We have isolated 1438 putative single units from 7 animals in 20 sessions (n = 1150 putative pyramidal cells, n = 288 putative interneurons) during the ThermoMaze behavior. We also collected 228 putative pyramidal cells from 2 animals in 3 control sessions (Suppl. Fig. 4) and 434 putative single units from 4 mice in 7 sessions using the 20-minute warmth paradigm (Fig. 6).
Pyramidal cells firing rate maps and SPW-R rate maps
To visualize and compare the spatial tuning properties of neurons across sub-sessions (Pre, Cooling and Post) during movement (speed ≥ 2.5 cm/s), we first binned the ThermoMaze ROI into 25 by 25 bins (each with size 1 x 1 cm) and counted the number of spikes of a neuron that occurred in each bin when the animal was actively moving (“movement spike-count map”). Next, we summed the total duration of time (in seconds) that the animal spent moving in each spatial bin to construct the “movement occupancy map”. The sub-session rate map of a cell during movement was computed by dividing the spike-count map by the occupancy map bin-wise. Similarly, we computed the SPW-R rate map within a subsession by dividing the number of ripples that occurred in each bin by the total duration of immobility (speed < 2.5 cm/s) that the animal spent in each bin. Both firing rate maps and SPW-R rate maps were spatially smoothed using a 2-bin smoothing window (https://github.com/buzsakilab/buzcode/blob/6418ba3b4307c673988bcf6ca44b15927fef5a7d/ externalPackages/FMAToolbox/Analyses/bz_Map.m).
Spatial tuning of spikes during SPW-Rs
To quantify spatial tuning of neurons during SPW-Rs (Figure 5), we defined a metric called “within-ripple spatial tuning score” which is a value between 0 and 1. The higher score indicates stronger spatial tuning of a neuron during SPW-Rs. We first binned the ThermoMaze ROI into four quadrants (2x2) and determined the firing rate of the neuron in each quadrant within SPW-Rs (i.e., total number of spikes of the cell divided by the total duration of SPW-R in that quadrant). For each SPW-R, a 300 ms time window surrounding the ripple’s power peak time was taken and the temporal overlaps between SPW-Rs were removed. Next, the within-SPW-R firing rate ratio in a given quadrant (e.g., in quadrant A), is defined to be the firing rate of the neuron during SPW-Rs in quadrant A divided by the sum of the within-SPW-R firing rate in all four quadrants. Finally, the within-ripple spatial tuning score (Figure 5) of a neuron is defined to be the maximum within-SPW-R firing rate ratio of the cell among all quadrants. To test the hypothesis that such spatial tuning exists beyond chance level, we generated shuffled within-SPW-R firing rate maps by randomly assigning one of the four quadrants to each SPW-R. Specifically, we randomly permuted the location of the SPW-Rs so that the number of SPW-Rs per quadrant was kept fixed for the shuffled condition.
Bayesian decoding of the animal position
Bayesian decoding of the animal’s position was based on the method provided by Zhang et. al., 1998)41. In short, we utilized the spatial firing rate maps constructed to find the location that maximally explains the observation of spiking within a certain time window. Because SPW-Rs occurred mainly in the corners where the warm spots were, we simplified the analysis and binned the ThermoMaze into 2x2 quadrants, which yielded four maze areas. We constructed the firing rate map templates fi(x) of each neuron during SPW-Rs (300 ms time window surrounding the peak of each SPW-Rs) within the Cooling sub-session. The decoded position was then determined to be the quadrant that maximizes the posterior likelihood given the observed spike counts:
where x was the quadrant index, n was the spike counts vector observed surrounding the frame time, τ was the time window size and equals 300 ms, C(τ,n) was a normalization factor and was taken to be 1, P(x) was the prior probability distribution of animal location and was taken to be 1 in the case of Figure 5D, i was the index of each cell, fi(x) was the average firing rate of cell i at position x, and N was the total number of pyramidal cells recorded in the session. For the purpose of cross-validation, we divided the SPW-Rs in each session into 100 folds. For each fold (testing dataset), the firing rate map templates were constructed using SPW-Rs from the other 99 folds (training dataset), and the decoding accuracy for the omitted fold was computed as the proportion of SPW-Rs whose corresponding quadrant was correctly decoded over the total number of SPW-Rs in the fold. For each session, we report the average decoding accuracy of test datasets.
Comparison of spatial tuning during SPW-Rs and movement
To quantify the similarity between spatial tuning of neurons during SPW-R and movement (theta oscillation), we calculated the firing rate ratios during movement in a similar way as we calculated the within-SPW-R firing rate ratios (see section “Spatial tuning during SPW-Rs” above). The ThermoMaze ROI was again binned into quadrants and firing rate maps (2x2) of each neuron during movement were calculated. The firing rate ratio of a neuron in each quadrant during movement was defined as the quadrant with the actual firing rate in that quadrant divided by its mean firing rate in all quadrants. Next, the Pearson correlation between the firing rate ratios during SPW-Rs and movement in each quadrant for each cell within the Cooling sub sessions were calculated.
We also studied the correlation between pyramidal cells’ spatial tuning during SPW-Rs and movement at a population level (Fig. 6G). In each session, we first constructed population vectors in each quadrant by concatenating the firing rate ratio of each cell in a quadrant into a vector during SPW-R or movement. We then computed the pairwise correlation coefficients between correlation matrix among the four population vectors between each condition and took the average across sessions.
References
- 1.Cognition from the Body-Brain Partnership: Exaptation of MemoryAnnu. Rev. Neurosci 46:191–210https://doi.org/10.1146/annurev-neuro-101222-110632
- 2.Levels of Processing: A Framework for Memory ResearchJ. Verbal Learning Verbal Behav 11:671–684https://doi.org/10.4324/9781315440446
- 3.Concomitant Cervical and Lumbar Stenosis: Strategies for Treatment and OutcomesElectroencephalogr. Clin. Neurophysiol 26:407–418https://doi.org/10.1053/j.semss.2007.06.005
- 4.Nucleus basalis and thalamic control of neocortical activity in the freely moving ratJ. Neurosci 8:4007–4026https://doi.org/10.1523/jneurosci.08-11-04007.1988
- 5.The effects of tonic locus ceruleus output on sensory-evoked responses of ventral posterior medial thalamic and barrel field cortical neurons in the awake ratJ. Neurosci 24:10773–10785https://doi.org/10.1523/JNEUROSCI.1573-04.2004
- 6.Neocortical activation: Modulation by multiple pathways acting on central cholinergic and serotonergic systemsExp. Brain Res 116:160–174https://doi.org/10.1007/PL00005736
- 7.Cortical state and attentionNat. Rev. Neurosci 12:509–523https://doi.org/10.1038/nrn3084
- 8.Cellular bases of neocortical activation: Modulation of neural oscillations by the nucleus basalis and endogenous acetylcholineJ. Neurosci 12:4701–4711https://doi.org/10.1523/jneurosci.12-12-04701.1992
- 9.Neuromodulation of Brain State and BehaviorAnnu. Rev. Neurosci 43:391–415https://doi.org/10.1146/annurev-neuro-100219-105424
- 10.Cellular bases of hippocampal EEG in the behaving ratBrain Res 287:139–171https://doi.org/10.1016/0165-0173(83)90037-1
- 11.Probing perceptual decisions in rodentsNat. Neurosci 16:824–831https://doi.org/10.1038/nn.3410
- 12.Exploratory behavior as a function of environmental novelty and complexity in male and female ratsPsychol. Rep 97:631–638https://doi.org/10.2466/pr0.97.2.631-638
- 13.Hippocampal place cells have goal-oriented vector fields during navigationNature 607:741–746https://doi.org/10.1038/s41586-022-04913-9
- 14.Three brain states in the hippocampus and cortexHippocampus 29:184–238https://doi.org/10.1002/hipo.22956
- 15.Reactivation of Hippocampal Ensemble Memories During SleepScience (80-. ) 265:676–679
- 16.Thirst regulates motivated behavior through modulation of brainwide neural population dynamicsScience (80-. ) 364https://doi.org/10.1126/science.aav3932
- 17.Food and Water Restriction Protocols: Physiological and Behavioral ConsiderationsContemp. Top. Lab. Anim. Sci 39:9–17
- 18.Health effects of water restriction to motivate lever-pressing in ratsLab. Anim. Sci 44:135–140
- 19.Defense of Water Balance in Rats: Behavioral and Physiological Responses To DepletionJ. Comp. Physiol. Psychol 64:59–67https://doi.org/10.1037/h0024800
- 20.Awake hippocampal reactivations project onto orthogonal neuronal assembliesScience 353:1280–1283https://doi.org/10.1126/science.aaf3319
- 21.Reverse replay of behavioural sequences in hippocampal place cells during the awake stateNature 440:680–683https://doi.org/10.1038/nature04587
- 22.A hippocampal network for spatial coding during immobility and sleepNature 531:185–190https://doi.org/10.1038/nature17144
- 23.Stress and anxiety: Structural plasticity and epigenetic regulation as a consequence of stressNeuropharmacology 62:3–12https://doi.org/10.1016/j.neuropharm.2011.07.014
- 24.Spatial selectivity of rat hippocampal neurons: Dependence on preparedness for movementScience (80-. ) 244https://doi.org/10.1126/science.2740902
- 25.Garber, J.C., Barbee, R.W., Bielitzki, J.T., Clayton, L.A., Donovan, J.C., Hendriksen, C.F.M., Kohn, D.F., Lipman, N.S., Locke, P.A., Melcher, J., et al. (1996). Guide for the care and use of laboratory animals Eighth. (The National Academic Press).Guide for the care and use of laboratory animals Eighth. (The National Academic Press)
- 26.Mouse Thermoregulation: Introducing the Concept of the Thermoneutral PointCell Rep 31:139–148https://doi.org/10.1053/j.gastro.2016.08.014.CagY
- 27.Translating animal model research: Does it matter that our rodents are cold?Physiology 29:413–420https://doi.org/10.1152/physiol.00029.2014
- 28.Gordon, C.J. (1993). Temperature Regulation in Laboratory Rodents (Cambridge University Press) 10.1017/cbo9780511565595.Temperature Regulation in Laboratory Rodents (Cambridge University Press) https://doi.org/10.1017/cbo9780511565595
- 29.Some like it hot: Mouse temperature preferences in laboratory housingAppl. Anim. Behav. Sci 116:279–285https://doi.org/10.1016/j.applanim.2008.10.002
- 30.Neuronal networks controlling thermoregulatory effectorsProg. Brain Res 115:49–62https://doi.org/10.1016/s0079-6123(08)62029-4
- 31.Working with what you’ve got: Changes in thermal preference and behavior in mice with or without nesting materialJ. Therm. Biol 36:193–199https://doi.org/10.1016/j.jtherbio.2011.02.004
- 32.Behavioral thermoregulatory responses of single- and group-housed micePhysiol. Behav 65:255–262https://doi.org/10.1016/S0031-9384(98)00148-6
- 33.New apparatus for studying behavioral thermoregulation in ratsPhysiol. Behav 64:419–424https://doi.org/10.1016/S0031-9384(98)00094-8
- 34.FLIR (2016). Infrared Camera Accuracy and Uncertainty in Plain Language. 1–4.Infrared Camera Accuracy and Uncertainty in Plain Language :1–4
- 35.The effects of changes in the environment on the spatial firing of hippocampal complex-spike cellsJ. Neurosci 7:1951–1968https://doi.org/10.1523/jneurosci.07-07-01951.1987
- 36.Brain temperature affects quantitative features of hippocampal sharp wave ripplesJ. Neurophysiol 127:1417–1425https://doi.org/10.1152/jn.00047.2022
- 37.Association between brain temperature and dentate field potentials in exploring and swimming ratsScience 259https://doi.org/10.1126/science.8446900
- 38.Brain temperature homeostasis: Physiological fluctuations and pathological shiftsFront. Biosci 15:73–92https://doi.org/10.2741/3608
- 39.Relationship between autonomic and behavioral thermoregulation in the golden hamsterPhysiol. Behav 34:687–690https://doi.org/10.1152/ajpregu.1986.251.2.r320
- 40.Developments of a water-maze procedure for studying spatial learning in the ratJ. Neurosci. Methods 11:47–60https://doi.org/10.1016/0165-0270(84)90007-4
- 41.Interpreting neuronal population activity by reconstruction: Unified framework with application to hippocampal place cellsJ. Neurophysiol 79:1017–1044https://doi.org/10.1152/jn.1998.79.2.1017
- 42.O’Keefe, J., and Nadel, L. (1978). The Hippocampus as a Cognitive Map (Oxford University Press) 10.5840/philstudies19802725.The Hippocampus as a Cognitive Map (Oxford University Press) https://doi.org/10.5840/philstudies19802725
- 43.Place-selective firing of CA1 pyramidal cells during sharp wave/ripple network patterns in exploratory behaviorNeuron 49:143–155https://doi.org/10.1016/j.neuron.2005.10.037
- 44.Network Homeostasis and State Dynamics of Neocortical SleepNeuron 90:839–852https://doi.org/10.1016/j.neuron.2016.03.036
- 45.The misbehavior of organismsAm. Psychol 16:681–684https://doi.org/10.1037/h0040090
- 46.The “Where is it?” reflex: Autoshaping the orienting responseJ. Exp. Anal. Behav 37:461–484
- 47.On the generality of the laws of learningPsychol. Rev 70:406–418https://doi.org/10.1037/h0069842
- 48.A study of misbehavior: Token reinforcement in the ratJ. Exp. Anal. Behav 29:115–134
- 49.Action-driven remapping of hippocampal neuronal populations in jumping ratsProc. Natl. Acad. Sci. U. S. A 119https://doi.org/10.1073/pnas.2122141119
- 50.Neuroscience Needs Behavior: Correcting a Reductionist BiasNeuron 93:480–490https://doi.org/10.1016/j.neuron.2016.12.041
- 51.Is coding a relevant metaphor for the brain?Behav. Brain Sci https://doi.org/10.1017/S0140525X19000049
- 52.A sensorimotor alternative to coding is possibleBehav. Brain Sci 42https://doi.org/10.1017/S0140525X19001468
- 53.Central Mechanisms for ThermoregulationAnnu. Rev. Physiol 81:285–308https://doi.org/10.1146/annurev-physiol-020518-114546
- 54.Recent advances in thermoregulationAdv. Physiol. Educ 39:139–148https://doi.org/10.1152/advan.00126.2014
- 55.Cold exposure potentiates the effect of insulin on in vivo glucose uptakeAm. J. Physiol. - Endocrinol. Metab 253https://doi.org/10.1152/ajpendo.1987.253.2
- 56.The sympathetic nervous system as a homeostatic mechanismI. Absolute need for sympathetic nervous function in body temperature maintenance of cold-exposed rats. J. Pharmacol. Exp. Ther 157
- 57.A Neuronal Hub Binding Sleep Initiation and Body Cooling in Response to a Warm External StimulusCurr. Biol 28:2263–2273https://doi.org/10.1016/j.cub.2018.05.054
- 58.The TRPM2 channel is a hypothalamic heat sensor that limits fever and can drive hypothermiaScience (80-. ) 353:1393–1398https://doi.org/10.1126/science.aaf7537
- 59.Allostasis: A new paradigm to explain arousal pathology. In Handbook of Life Stress, Cognition and Health, S. Fisher and J. Reason, edsJohn Wiley & Sons, Ltd
- 60.Principles of allostasis: Optimal design, predictive regulation, pathophysiology, and rational therapeutics. Allostasis, HomeostasisCosts Physiol. Adapt https://doi.org/10.1017/CBO9781316257081.004
- 61.Allostasis: A Brain-Centered, Predictive Mode of Physiological RegulationTrends Neurosci 42:740–752https://doi.org/10.1016/j.tins.2019.07.010
- 62.Regulation of Body Temperature by the Nervous SystemNeuron 98:31–48https://doi.org/10.1016/j.neuron.2018.02.022
- 63.Animal models of anxiety disorders in rats and mice: Some conceptual issuesDialogues Clin. Neurosci 13:495–506https://doi.org/10.31887/dcns.2011.13.4/tsteimer
- 64.Forward and reverse hippocampal place-cell sequences during ripplesNat. Neurosci 10:1241–1242https://doi.org/10.1038/nn1961
- 65.The reorganization and reactivation of hippocampal maps predict spatial memory performanceNat. Neurosci 13:995–1002https://doi.org/10.1038/nn.2599
- 66.Trajectory events across hippocampal place cells require previous experienceNat. Neurosci 18:1772–1779https://doi.org/10.1038/nn.4151
- 67.Hippocampal place-cell sequences depict future paths to remembered goalsNature 497:74–79https://doi.org/10.1038/nature12112
- 68.Hippocampal Replay of Extended ExperienceNeuron 63:497–507https://doi.org/10.1016/j.neuron.2009.07.027
- 69.Hippocampal Replay Is Not a Simple Function of ExperienceNeuron 65:695–705https://doi.org/10.1016/j.neuron.2010.01.034
- 70.Awake replay of remote experiences in the hippocampusNat. Neurosci 12:913–918https://doi.org/10.1038/nn.2344
- 71.Preplay of future place cell sequences by hippocampal cellular assembliesNature 469:397–401https://doi.org/10.1038/nature09633
- 72.Diversity in neural firing dynamics supports both rigid and learned hippocampal sequencesScience (80-. ) 351:1440–1443https://doi.org/10.1126/science.aad1935
- 73.Level of arousal during the small irregular activity state in the rat hippocampal EEGJ. Neurophysiol 91:2649–2657https://doi.org/10.1152/jn.01082.2003
- 74.3D-printed Recoverable Microdrive and Base Plate System for Rodent ElectrophysiologyBio-Protocol 11
- 75.Metal microdrive and head cap system for silicon probe recovery in freely moving rodentElife 10:1–21https://doi.org/10.7554/eLife.65859
- 76.DeepLabCut: markerless pose estimation of user-defined body parts with deep learningNat. Neurosci 21:1281–1289https://doi.org/10.1038/s41593-018-0209-y
- 77.CellExplorer: A framework for visualizing and characterizing single neuronsNeuron 109:3594–3608https://doi.org/10.1016/j.neuron.2021.09.002
- 78.Kilosort: realtime spike-sorting for extracellular electrophysiology with hundreds of channels. bioRxiv
Article and author information
Author information
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
Copyright
© 2023, Vöröslakos et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 537
- downloads
- 7
- citation
- 1
Views, downloads and citations are aggregated across all versions of this paper published by eLife.