1. Neuroscience
Download icon

Hierarchical organization of context in the hippocampal episodic code

  1. Susumu Takahashi Is a corresponding author
  1. Doshisha University, Japan
  2. Japan Science and Technology Agency, Japan
  3. Kyoto Sangyo University, Japan
Research Article
Cited
9
Views
1,007
Comments
0
Cite as: eLife 2013;2:e00321 doi: 10.7554/eLife.00321

Abstract

The hippocampal system appears to be critically important in establishing episodic memory of both internal and external events within contexts as well as spatial memory, which enables flexible spatial navigation. However, the neuronal substrates that function across different memories in the hippocampal system are poorly understood. I monitored large-scale activity patterns of hippocampal neuronal ensembles in rats performing a novel, continuous task that combined one visually guided and two memory-guided types of navigations in a constant environment. I found that the activity patterns of the hippocampal ensemble represent spatiotemporal contexts (journeys) constructed by temporally ordered past, present and expected future places in tandem with visually or mnemonically guided non-spatial contexts (task-demands) to form episodes. This finding therefore suggests that the hierarchical organization of contexts based on pattern separation and completion enables the hippocampus to play a dual role in spatial navigation and recall of episodic memory.

https://doi.org/10.7554/eLife.00321.001

eLife digest

A little over 10 years ago, researchers discovered that a brain region called the hippocampus is larger in London taxi drivers than it is in the general population. This tied in with results from animal studies, which had revealed a key role for the hippocampus in spatial navigation and memory. However, other work has shown that the hippocampus is equally important for remembering personal experiences—a form of memory known as episodic memory.

Many neurons in the hippocampus display ‘place fields', which means that they fire bursts of action potentials whenever an animal is in a specific location. Place fields tend to remain stable during repeated visits to an environment: the same cells fire whenever the animal returns to a particular place. However, if the animal enters a new environment, a neuron might adopt a different place field or not show any place field at all. This phenomenon is known as remapping.

Now, Takahashi has provided further insight into the circumstances under which such remapping occurs. He recorded from large numbers of neurons in the rat hippocampus—in a subregion called CA1—as the animals moved through a maze shaped like a digital figure ‘8'. The rats had to perform three tasks within the maze: one guided by visual cues, and two that were memory-based.

In the visual task, a light informed the rats to turn either left or right to obtain a reward. In the first memory task, the rats had to alternate their choices to obtain the reward, running through the maze from right-to-left and then from left-to-right (non-delayed spatial alternation). The second memory task worked the same way, except that the rats had to wait 5 s before turning left or right (delayed spatial alternation).

Takahashi compared the responses of hippocampal CA1 neurons as rats performed the three tasks. As expected, he found that neurons tended to remap their place fields based on the animal's initial and final locations in the maze, regardless of which task the animal was performing. Surprisingly, however, neurons with specific place fields distinguished between the three tasks by firing at different rates in each.

By combining information about the locations and rates at which large assemblies of neurons fired, Takahashi found that he could accurately predict which task a rat had been performing, where it had come from, and where it had ended up, because the place field remapping was hierarchically organized. Moreover, the prediction could be made even before the rat had completed the task. Overall, these results add to our understanding of how the hippocampus performs its dual roles in spatial navigation and episodic memory.

https://doi.org/10.7554/eLife.00321.002

Introduction

There have been two frameworks suggested for the principal role of the hippocampus in memory: spatial memory and episodic memory (Scoville and Milner, 1957; O'Keefe and Nadel, 1978; Vargha-Khadem et al., 1997; Eichenbaum et al., 1999). During spatial behaviors, hippocampal pyramidal neurons discharge specifically at a certain location in the environment (the ‘place field'), suggesting that hippocampal neurons function in spatial memory (O'Keefe and Dostrovsky, 1971; McNaughton and Morris, 1987; Muller et al., 1987). In contrast, episodic memory provides a record of past experiences in relation to a broad range of materials (Wood et al., 1999) and is structured by spatial, temporal, and non-spatial contexts (Tulving, 1983). Recent studies indicate that on the basis of the pattern separation and completion proposed by Marr (Marr, 1969), the place field of a given neuron changes (remapping) relative to three different contexts: spatial context characterized by environmental features (Bostock et al., 1991; Kentros et al., 1998), spatiotemporal context (journey) affected by the experienced origin and expected destination (Frank et al., 2000; Wood et al., 2000), and internally and externally guided non-spatial context (task-demand) (Gothard et al., 1996; Anderson and Jeffery, 2003; Smith and Mizumori, 2006). Moreover, during a temporal gap between events, the place cells encode successive moments or fire at a specific part of an episode (Pastalkova et al., 2008; MacDonald et al., 2011). These lines of evidence suggest that the hippocampal place code is not only closely related to the ‘where' of episodic memory but also to the ‘when and how'. To test whether the changes and dynamics of the hippocampal place code are the neuronal basis of spatial and episodic-like memory, it is necessary to conduct a multifaceted experiment in which many hippocampal neurons are simultaneously monitored during repeated exposures to either spatiotemporal or non-spatial contexts in a constant spatial environment.

Results

A continuous task with one visually guided and two memory-guided decisions

Rats were trained to navigate their way through a figure ‘8' maze in a continuous task that incorporated both visual discrimination and two types of memory-guided responses (see ‘Materials and methods', Figure 1). The task was divided into three subtasks: visual discrimination (VD), non-delayed spatial alternation (SA) and delayed spatial alternation (DA). In the VD subtask, a cue light on the left or right side was illuminated to indicate which direction to turn at a decision point (DP, Figure 1A) to receive a reward (medial forebrain bundle [MFB] electrical stimulation). In both the SA and DA subtasks, cue lights on both sides were illuminated and the rats had to alternate between left and right at the decision point (Figure 1B). The difference between the SA and DA subtasks was that in the DA subtask, a barrier appeared at the middle of the central stem for 5 s to interpose a delay period. During the delay period, the rats paused steadily in the forward direction in front of the barrier. Some lines of evidence suggest that the SA-type task is often affected by hippocampal lesions (O'Keefe and Nadel, 1978). By contrast, a recent study implied that the SA-type task is independent on the hippocampus (Ainge et al., 2007b). Unlike the visually cued VD and procedural memory-guided SA subtasks, the episodic-like memory-guided DA type of task is critically dependent on the hippocampus (Ainge et al., 2007b; Ferbinteanu et al., 2011). The rats continuously performed these subtasks at least twice in a task over the course of approximately 1 hr (Figure 1D). This task allowed for examinations of the influence of internal and external events on hippocampal neuronal activity and hippocampus-dependence both between and within different journeys and task-demands.

Task configuration.

(A) The rat has to decide which direction to turn at the decision point (DP). If correct, the rat gets reward signals (MFB stimulations) at one of the reward zones (RZ). In the DA subtask, a barrier appears for 5 s in the central stem (blue line) so that the rat has to wait in the delay zone (DZ). An example of a running trajectory of a rat (G111125) is superimposed on the maze. (B) Possible journeys for three subtasks. For the VD subtask, one of the light cues is illuminated to indicate a direction. For the memory-guided (SA and DA) subtasks, both of the light cues are illuminated to alternate the direction based on memory. (C) Examples of running trajectories of the rat (G111125). The red dots indicate the location of the onset of the estimated turn for each lap. (D) The subtask sequence in a task. Each subtask consisted of 20 laps, except the second VD subtask in each session, which only consisted of 10 laps.

https://doi.org/10.7554/eLife.00321.003

Coding of place, journey and task-demand in the hippocampal place cells

Ensemble activity was recorded from a total of 1119 pyramidal cells in the CA1 of the dorsal hippocampus using arrays of 10 extracellular dodecatrodes (Takahashi and Sakurai, 2005) in five rats running in the maze. All rats were familiarized with the task before recording, and performed the task with a high level of accuracy (overall: 99 ± 2%, VD: 99 ± 1%, SA: 97 ± 3%, DA: 99 ± 1% correctly performed [mean ± SEM]; Table 1). Laps containing errors of cue interpretation were excluded from the following analyses.

Table 1

Variability between rats: behavioral and electrophysiological measurements

https://doi.org/10.7554/eLife.00321.004
Rat No.B120303F120112G111125K120224L120221Total
No. of erroneous laps01051
Errors (%)00.703.50.7
Running speed (cm/sec ± SD) in the central stem33.13 ± 6.1745.61 ± 11.7322.09 ± 15.0841.77 ± 11.7752.03 ± 16.95
No. of pyramidal cells238744121932021119
No. of place cells (spatial information > 1.0 bits/spike)1272119959116522
Spatial information (bits/spike, mean ± SEM)2.39 ± 0.211.26 ± 0.362.47 ± 0.211.21 ± 0.182.50 ± 0.24
Unit isolation quality (isolation distance, mean ± SEM)22.42 ± 1.7027.69 ± 2.8145.41 ± 7.8623.26 ± 1.5933.60 ± 2.11
No. of interneurons19524151073

Several studies suggest that the hippocampal place code is influenced by running speed and head direction (McNaughton et al., 1983; Wiener et al., 1989; Muller et al., 1994; Czurko et al., 1999). In addition, even small differences in location on a maze can strongly influence firing rates (Muller et al., 1994). Although in the present work, these factors were considered to be minimized in the central stem of the maze because the running speed, head direction, and lateral position appeared highly consistent, an analysis of the place cells was employed between different journeys within a given subtask or between different task-demands within a journey in order to quantitatively account for any remaining effects of these potentially confounding factors. In this analysis, for each place cell in which the spatial information in at least one trial type was above 1.0 bits/spike, differences in firing rates associated with those differences in the central stem were analyzed using a one-way ANCOVA. Running speed, head direction and lateral position were included as covariates in the ANCOVA model. Firing rates of over 268 of the 522 place cells in the central stem differed significantly between different journeys within a given subtask or between different task-demands within a journey (p<0.05; Figure 2). To quantitatively examine the remapping of place fields without influences of the running speed, head direction, and lateral position, only place cells which fired at significantly different rates in the central stem on different journeys within a given subtask or different task-demands within a journey, even when all of these factors were taken into account, were examined in the following analyses.

Proportion of differentially firing place cells.

The proportion of place cells showing significantly different firing rates between different journeys within a given subtask and between different task-demands within a journey when running speed, head direction and lateral position were taken into account. Numbers indicate percentages.

https://doi.org/10.7554/eLife.00321.005

The place codes could be changed in terms of firing rates (rate remapping) or both firing rates and location (global/complete remapping; Bostock et al., 1991; Leutgeb et al., 2005). For each place cell, I independently compared the firing location and rate of the place fields during task performance. The spatial similarity was measured by calculating the spatial correlation between a pair of maps. The rate similarity was measured as one minus the normalized difference between the firing rates. Thus, the spatial similarity measurement was independent of rate similarity. In an identical trial type (i.e., combination of journey and task-demand), the distributions of the measurements between first and second exposures (i.e., first vs second exposure in the right-to-left or left-to-right journey in the VD, SA or DA subtasks) were highly similar (median Spearman correlations [rs] > 0.67; median rate similarity indices [rr] > 0.80; Figures 3A, 4B,C), suggesting that, as in numerous previous studies on the hippocampal place code (O'Keefe and Dostrovsky, 1971; McNaughton and Morris, 1987; Muller et al., 1987), the pyramidal cells are place-specific during identical visually or mnemonically guided behaviors in a constant environment. Within a given task-demand, the distributions of the spatial correlation and the change in firing rates in the place fields in the central stem between different journeys (i.e., a right-to-left vs left-to-right journey in the VD, SA, or DA subtasks) were significantly lower than those seen under control conditions (repeated exposures to identical trial types) (median rs < 0.45, Z < −22.04, p<0.001; median rr < 0.68, Z < −8.96, p<0.001, Wilcoxon rank sum test; Figures 3B and 4B,C), suggesting that place-specific activity is strongly dependent on the journey (Frank et al., 2000; Wood et al., 2000; Ferbinteanu and Shapiro, 2003; Bower et al., 2005; Ainge et al., 2007a). In contrast, in an identical journey, the distributions of the spatial correlations among the VD, SA and DA subtasks (i.e., VD vs SA subtask, VD vs DA subtask, SA vs DA subtask in the right-to-left or left-to-right journey) were similar or higher than under control conditions (median rs > 0.75, Z > 1.07, n.s. or p<0.001, Wilcoxon rank sum test; Figures 3B and 4B); however, the distributions of the change in firing rates in the place fields were clearly lower (median rr < 0.80, Z < −2.79, p<0.01, Wilcoxon rank sum test; Figures 3B and 4C), indicating that only the rate of a place-specific activity is modulated with a given task-demand (Smith and Mizumori, 2006). Taken together, the results suggest that irrespective of hippocampus-dependence (Morris and Frey, 1997; Ainge et al., 2007b; Ferbinteanu et al., 2011), the hippocampal pyramidal cells in the maze were not only activated by a specific location, but also by the experienced past and expected future places accompanied by the visually guided or memory-guided task-demands (see Figure 3—figure supplement 1).

Figure 3 with 1 supplement see all
Journey-dependent coding of a place-specific activity modulated with task-demand.

Representative color-coded rate maps for eight possible trial types. (A) The rate maps with the first cell in (B) show stable place-specific activity in repeated sessions of an identical journey within a given task-demand. The rate was coded on a color scale from blue (silent) to red (maximum rate). The pixels that were not sampled are white. The symbols above the map indicate the journey, task-demand and illuminated light-cue. For each cell, the rate scale corresponds to the peak firing rate in Hz in that condition (indicated to the right of the rate maps). (B) The rate maps for six CA1 pyramidal cells. The maps are averages for repeated sessions. Each row shows the data from one pyramidal cell. Note that different journeys in the same place resulted in place fields similar in locations, but different task-demands caused different firing intensities in those place fields.

https://doi.org/10.7554/eLife.00321.006
A quantitative assessment of the spatial and rate similarities between a pair of trial types.

(A) The similarities were measured between a pair of trial types (VD1, VD2, VD3, VD4, SA1, SA2, DA1, and DA2) or between repeated exposures to an identical trial type. Because the central stem of the maze was the common running route among the different journeys within a given subtask, the region of interest (ROI) was set at the region highlighted in the green shaded box. Box plots of the spatial (B) and rate (C) similarities between repeated exposures to an identical journey and task-demand (red), between different journeys within a given task-demand (green), and between different task-demands within a journey (blue) (median, first and third quartiles, minimum, and maximum indicated). Repeated exposures are marked with an apostrophe.

https://doi.org/10.7554/eLife.00321.008

Episode-specific neuronal trajectory of the hippocampal ensemble activity

The activity patterns relative to the place, journey and task-demand during tasks suggested that the dynamics of the ensemble activity pattern rather than just the firing pattern of single cells should be considered. The dynamics of neuronal ensemble activity as a trajectory through a state space were therefore analyzed (neuronal trajectory; see ‘Materials and methods'; Figure 5) (Briggman et al., 2005; Mazor and Laurent, 2005; Harvey et al., 2012). At each specific location, the activity state of the ensemble containing n simultaneously recorded neurons was defined as a point in an n-dimensional space, with each dimension representing the activity of a single neuron. In the trajectory of a journey within a given subtask, the neuronal ensemble activity tended to move in similar orbits that started and ended at similar positions. Moreover, the trajectories of different journeys exhibited unique orbits (Figure 6A,B). To quantify the trajectory specificity, I used a classifier based on the distance from an individual lap trajectory to the mean right-to-left/left-to-right trajectories at single locations. In each subtask, the neuronal trajectory, consisting of all of the simultaneously recorded place cells, was sufficient to separate the left-to-right journey from the right-to-left journey with high accuracy, and vice versa (Figure 6C; p<0.001, binomial test, chance = 0.5). Surprisingly, even in visually guided situations in which memory-guided decisions could not be utilized (i.e., left-to-left/right-to-right journey in the VD subtask), the accuracy was sufficiently better than chance levels before the rat behaviorally chose its next direction (prospective coding; Figure 6D; p<0.001, chance = 0.5; turn onset: Figure 1C). Furthermore, this trajectory could be used to predict all possible trial types (i.e., eight possibilities; VD1, VD2, VD3, VD4, SA1, SA2, DA1 and DA2) on single laps at better than chance levels even during periods of running through the central stem of the maze in which the head direction, speed and place factors were almost constant for each rat (Figure 7, p<0.001, binomial test, chance = 0.125; Table 1). Thus, the activity in the hippocampus could be considered to be divergent, episode-specific trajectories through a state space of neuronal ensemble activity (see Figure 7—figure supplement 1).

A schematic of a classification scheme based on the distances to the mean trajectories.

At each location, the Euclidian distance from the trajectory on a single lap (black) to the mean trial type trajectories (red and blue) was calculated (dA, dB). If dA was less than dB, the lap at the location bin was classified as trial type A, and vice versa for trial type B.

https://doi.org/10.7554/eLife.00321.009
Neuronal trajectory of different journeys and different future directions.

(A) An example time course of mean, journey-specific trajectories for right-to-left (red) and left-to-right (blue) laps in the SA subtask, plotted on the first to third principal component space (Rat #: B120303, n = 91 cells). The points marked 1, 2 and 3 correspond to the mean locations where the rat entered the lower arm, and turned and entered the reward zone, respectively. (B) Example individual trajectories for laps in the subtask in (A). (C) The classification accuracy of determining the journeys (shown in right) at different locations in the VD, SA, and DA subtasks (thin black line, n = 5 rats; solid black line, mean; blue dashed line, p=0.001, binomial test; black dashed line, chance level; red vertical line, mean turn onset). The classifier was based on a distance-dependent classification scheme. (D) As (C), except that the subtask condition under which the future direction is not predictable until arriving at the decision point. A high level of accuracy can be observed even before the onset of the turn (arrow heads).

https://doi.org/10.7554/eLife.00321.010
Figure 7 with 1 supplement see all
Episode-specific neuronal trajectories.

The classification accuracy for determining a trial type from all of the possible trial types (episodes) (thin line, n = 5 rats; thick black line, mean; blue dashed line, p=0.001, binomial test; black dashed line, chance level). The classifier was based on a distance-dependent classification scheme. The average accuracy for each rat was significantly higher than chance at the respective specific locations.

https://doi.org/10.7554/eLife.00321.011

Relative ability for distinguishing journeys, task-demands, and trial types in the hippocampal episodic code

The above results suggest that hippocampal ensemble activity contains sufficient information concerning journey, task-demand and trial type. To compare the relative information content among them for neuronal ensemble activity of each rat, I quantitatively assessed the performances of the binary classification of neuronal trajectories in the central stem of the maze between a pair of journeys irrespective of subtask, between a pair of task-demands irrespective of journey, and between a pair of trial types. For all categories, the average output of the binary classifiers showed a high level of accuracy at the ensemble level (Figure 8, p <0.001, binomial test, chance = 0.5) and no significant differences among them (p >0.05, Kruskal–Wallis test with Turkey honestly significant difference [HSD] post-hoc test). The results suggest that the abilities of the hippocampal ensemble to identify journeys, task-demands, and trial types are sufficiently high and relatively similar to each other.

Comparison of accuracy of distinguishing journeys, task-demands and trial types from neuronal ensemble activity.

Average binary classification accuracies of distinguishing pairs of journeys (i.e., VD1, SA1, and DA1 vs VD3, SA2, and DA2), pairs of task-demands (i.e., VD vs SA, VD vs DA, and SA vs DA), and pairs of trial types (i.e., all possible pairings of VD1, VD2, VD3, VD4, SA1, SA2, DA1, and DA2). For each rat, the accuracy was calculated based on the most frequently classified choices in the central stem. The error bars indicate SD.

https://doi.org/10.7554/eLife.00321.013

Hierarchical representation of journey and task-demand in hippocampal ensemble activity

The results obtained showing the similarities of the place fields may imply the pattern completion/generalization of the journey represented in the hippocampal ensemble activity in the spatial domain. If this assumption is correct, knowing the neuronal trajectories for right-to-left and left-to-right journeys in one subtask should allow classification into right-to-left and left-to-right journeys in the other subtasks. To test this, the similarities between the ensemble activity patterns in the central stem of the maze were quantitatively analyzed using a distance-dependent classification scheme. As expected, based on the mean right-to-left/left-to-right trajectories in any one of three subtasks, the classifier outputs could predict the corresponding individual journeys within the other remaining two subtasks at a high level of accuracy at the ensemble level (see ‘Materials and methods'; Figure 9, p<0.001, binomial test, chance = 0.5, n = 128 cells). These results suggest that hippocampal ensemble activity generalizes a journey over the given task-demands.

Generalization accuracy of determining a journey over task-demands from neuronal ensemble activity.

(A) Schematic of a classification scheme for estimating generalization accuracy. At each location, the distance from the trajectory in two target subtasks on a single trial to the mean right-to-left (R-L) and left-to-right (L-R) journey trajectories in a reference subtask was measured (d-to-RL and d-to-LR). If the actual journey was R-L, and d-to-RL was less than d-to-LR, the lap at that location was classified as a correct generalization, and vice versa. The figure shows an example in which a reference subtask is VD and a neuronal trajectory in the R-L journey of an SA subtask was examined. (B) The generalization accuracy for determining a R-L/L-R journey from the neuronal trajectory in each lap in two of the three subtasks based on the distance to the mean neuronal trajectories in a remaining reference subtask as a function of the number of cells (red: reference: VD, target: SA and DA; green: reference: SA, target: VD and DA; blue: reference DA, target: VD and SA). Cells recorded from five rats were combined. The accuracy was calculated based on the most frequently classified journeys in the central stem. The error bars show the SD for 1000 random choices of the cells used; the dashed line shows the chance level, and the bar next to the dashed line shows the range of accuracies using random shuffling (control).

https://doi.org/10.7554/eLife.00321.014

Discussion

To summarize, although the same prominent environmental features were encountered in the same locations, the spatial representation of the hippocampal CA1 at various points in the maze were dissimilar in both firing location and rate among different journeys, irrespective of either visually or mnemonically guided demands, in a global remapping manner. The result supports the view that hippocampal pyramidal cells show journey-dependent coding during internally or externally guided goal-directed behaviors (Ferbinteanu et al., 2011). However, when the rats experienced changes in prominent, non-spatial features with internal and external events (i.e., subtask differences), the hippocampus only modulated the journey-dependent activity by primarily changing the intensity of the firing rates. The latter finding indicates that the extent to which the place fields remap is demand-specific (Smith and Mizumori, 2006) in a rate remapping manner. Neither the journey-dependent nor demand-specific coding was chiefly influenced by the running speed, head direction, or lateral position. In addition to the evidence from the place field measurements, even when the given task-demands were different, the classifier could sufficiently predict the journey from the ensemble activity patterns, suggesting that the journey representation is generalized and that non-spatial demand-specific representation is hierarchically ranked at a lower level. The neuronal trajectory of the ensemble activity differed enough to classify any trial type; prior to its behavioral choice, significant accuracy could be achieved even in a visually-guided situation in which memory-guided decision could not be made. The neuronal ensemble code for distinguishing trial types in the hippocampal ensemble activity was relatively equal to that for journeys and task-demands. Thus, the results imply that based not only on externally induced but also internally generated events, the hippocampus processes the combination of journey and task-demand as a single episode.

Previous studies have reported that hippocampal neurons might convey information related to where an animal is located (place-specific coding) (O'Keefe and Dostrovsky, 1971; McNaughton and Morris, 1987; Muller et al., 1987; Bostock et al., 1991; Kentros et al., 1998), where it had just been (retrospective coding) (Frank et al., 2000; Wood et al., 2000; Ferbinteanu and Shapiro, 2003; Ferbinteanu et al., 2011), where it is about to go next (prospective coding) (Frank et al., 2000; Wood et al., 2000; Ferbinteanu and Shapiro, 2003; Ainge et al., 2007a; Ferbinteanu et al., 2011) and what demand it currently faces (demand-specific coding) (Gothard et al., 1996; Anderson and Jeffery, 2003; Smith and Mizumori, 2006). However, the relationship among them in the hippocampal neuronal ensemble codes was unclear in these studies because each piece of information was independently examined. I report that all of this information can be simultaneously preserved in the hippocampal output to form episodes. The information for journey and task-demand is not horizontally but hierarchically organized; expected future events in the hippocampus can trigger episode separation. From the viewpoint of context, journey and task-demand can be considered to be spatiotemporal and non-spatial contexts, respectively. Thus, the hierarchical organization of spatiotemporal and non-spatial contexts in the hippocampal episodic code, which is based on the context version of pattern separation and completion through the remapping of place fields, not only enables flexible spatial navigation but also contributes to the recall of episodic memory.

Materials and methods

Animals

Five male Wister rats (300–400 g) were housed individually in cages (20 × 25 × 23 cm). All of them were implanted with a custom-made microdrive with 10 dodecatrodes and 2 reference electrodes. The rats were kept at 80% of free-feeding body weight. All rats were maintained on a 12-hr light/12-hr dark schedule. Testing occurred in the light phase.

Surgery, electrode preparation and recording

Under isoflurane anesthesia, two stimulation electrodes were inserted into the MFB in the right lateral hypothalamus (AP 2.5, ML 1.0, DV 9.5); a custom-made microdrive with 10 independently movable 12-wire bundled dodecatrodes and 2 tetrodes were fixed to the skull above the left hippocampus (AP 3.8, ML 3.0, DV 0.5). Two tetrodes were used to record a reference signal. Each stimulation electrode consisted of two insulated stainless steel wires (coated diameter, 0.2 mm) together with ∼0.5 mm of insulation removed from one tip. The 10 dodecatrodes were lowered into the hippocampal CA1 pyramidal layer within 2 weeks of surgery. In all cases, the unit-recording dodecatrodes and reference tetrodes were constructed from 8 μm and 12.5 μm HML-coated tungsten (99.95%) wires, respectively.

After the rats had made a full recovery (∼2 weeks after surgery), they were allowed to explore the task maze and identify the reward zones (Figure 1A, RZ). During the physiological recording, the data from all channels of the dodecatrodes were unity-gain buffered, filtered (600 Hz–6 kHz), amplified (gain = 5000) and continuously sampled at 25 kHz on a 128-channel custom-made recording system with Lynx-8 amplifiers (Neuralynx, Bozeman, MT) and AD converters (16-bit resolution, National Instruments, Austin, TX). One channel of each dodecatrode was branched away, filtered (0.1 Hz–6 kHz) and amplified (gain = 500) to detect local field potentials (LFP). The recordings were included in the data analysis if sharp wave ripple events in the LFP in the corresponding dodecatrode were identified during the immobility period. Two sets of small red and green light-emitting diodes (10-cm separation), mounted above the headstage amplifier, were recorded by an overhead digital video camera and sampled at 60 Hz to track the rat's head position and direction. The neuronal data were recorded on a single recording day for each rat.

MFB stimulation

After the surgery, the optimal stimulation parameters were obtained when the rat performed nose poke responses. The optimal MFB stimulation consisted of a train of 2 ms wide, 100–200 μA, biphasic current pulses, delivered at 100 Hz for 200 ms. In the correct laps, the optimal MFB simulation was obtained 5–10 times in the reward zone. The neuronal recordings during the stimulation period were removed prior to analysis.

Task and behavioral training

The task demanded a visually or memory-guided decision in VD, SA and DA at a decision point in the figure 8-shaped maze (overall: 100 cm by 140 cm, 20 cm height; path: 20 cm width; matte black acryl resin, Figure 1A). To prevent the rats from accessing unnecessary distal room cues, the task was performed under dim light. Initially, two barriers were put into the maze to construct an O-shaped track. After electrode implantation, the rats were trained to run along the left/right O-shaped track unidirectionally for reward signals (MFB stimulation) at a reward zone (Figure 1A, RZ). Once the rats were running smoothly along both left and right O-shaped tracks (∼3 to 7 days), all of them were trained to run in the VD subtask for a reward in the figure 8-shaped maze. In the VD subtask, one of the two visual cues (LED-lights at the right and left corners; Figure 1A) was illuminated randomly; the visual cue at the decision point of the maze (Figure 1A, DP) indicated which direction to turn to receive reward signals. The rats did not access the visual cue until reaching the decision point (Figure 1A) because the height of the walls of the maze (20 cm) was much higher than the level of the rat's eyes (∼5 cm) and the rats did not stand up during the experiment. The rats were trained for ∼3 to 7 days until they achieved a threshold of 80% decisions correct for >20 laps at a constant running speed (>20 cm/s). Next, they were trained to run in the SA subtask for >20 laps. In the SA subtask, both LEDs were illuminated so that the rats could not rely on the visual cue. Instead, the rat had to choose a direction opposite to the previous one. Those alternating behaviors spontaneously occurred and the rats easily achieved the threshold of 80% decisions correct for >20 laps (∼1 to 3 days). Finally, the rats were trained to run in the DA subtask for 20 laps. The DA subtask was almost identical to the SA subtask, except that a 5-s delay period was incorporated. During the delay period, a barrier appeared for 5 s, 20 cm ahead of the entrance to the central stem (Figure 1A, blue line). Because the rats initially tended to move back from the barrier during the first few laps of the DA subtask, two other barriers were set at both sides of the return rail to prevent such behavior until it ceased. The rats achieved a threshold of 90% decisions correct for >20 laps (1 day).

After the rats correctly performed all three subtasks, they continuously performed all the subtasks within about 1 hr in the following sequence—VD (20 laps), SA (20 laps), VD (10 laps), DA (20 laps), VD (20 laps), SA (20 laps), VD (10 laps), and DA (20 laps)—and were rewarded each time they arrived at the reward zone at the correct return rail. The unit-recordings were made within 2–6 training days after each rat conducted this final sequence of subtasks. The number of possible journeys in the VD subtask was twice that of the other subtasks. Thus, to increase the sample size for the VD subtask for analysis, and to clearly distinguish the behaviorally similar SA and DA subtasks, 10 laps of the VD subtask were interpolated between the SA and DA subtasks. The rats achieved a 95% correct decision rating from the beginning of the continuous task.

Spike sorting

Spike sorting was performed offline using an automatic spike sorting software program with KlustaKwik (Harris et al., 2000) and FastICA (Hyvarinen, 1999) called ‘ICSort', which is a custom-made program using Matlab and the C++ programming language. ICSort completed two steps: clustering and separation. Initially, in the clustering step, ICSort automatically sorted the first, second and third principal components' feature vectors of all the extracted spikes in the high-dimensional clustering space using KlustaKwik. Then, the extracted spike waveforms in each sorted cluster were concatenated. Next, in the separation step, ICSort used FastICA to find spatial filters to separate the spike waveforms into spatially fixed and distinct, maximally independent components (ICs). To determine which ICs were indicative of neuronal activity, spikes of >50 μV were extracted from each IC and checked to determine whether there was a clear refractory period in the inter-spike interval (ISI) histogram. Finally, the ICSort aggregated these distributed ICs into single units based on each ICA basis vector, which represents the relative distance between the source of the IC and the tip of the electrode. This approach has been previously described in detail (Takahashi et al., 2003a, 2003b; Takahashi and Sakurai, 2007, 2009a, 2009b).

Place maps

The spatial firing rate distributions (‘place fields') for each cell were constructed in the standard manner by summing the total number of spikes that occurred in a given location bin (5 cm by 5 cm), dividing by the amount of time that the rat spent in that location, and smoothing with a Gaussian centered on each bin (Leutgeb et al., 2005). The firing rate at a point, x, was calculated as:

f(x)=i=1nw(sixh)/0Tw(y(t)xh)dt,

where n is the number of spikes, si the position of the i-th spike, y(t) the position of the rat at time t, and (0, T) the period of the recording. The kernel function, w, was a Gaussian of h = 5 cm width. Positions more than 5 cm away from the tracked path were regarded as unvisited.

Cell identification and recording periods

After spike sorting, putative pyramidal cells were distinguished from putative fast spiking interneurons by the spike width (0.4 ms) and the average firing rate (5 Hz). The low rate cells (<0.1 Hz) were excluded. Only cells identified as pyramidal were used in the analyses outlined below.

It is well known that during immobility periods, including delay periods, hippocampal activity is not place-specific (McNaughton et al., 1983; Pastalkova et al., 2008). Thus, spikes during immobility periods in the maze were excluded.

Spatial information

The spatial information was calculated as the number of bits per spike according to the formula (Skaggs et al., 1993):

Spatial Information=iPi(Ri/R)log2(Ri/R),

where i indexes over the position bins in the trial type, Pi is the probability that the rat was in bin i, Ri the mean firing rate in bin i, and R the overall mean firing rate over the trial type. More spatially specific firing leads to a larger value for this measurement.

Statistical comparison of firing rates using ANCOVA

To analyze the activity of cells in the central stem of the maze, the central stem was divided into spatial bins of equal length (approximately 0.35 cm). The following parameters were calculated for each traversal through each of the bins: (1) firing rate was calculated as the spike count in the bin divided by the time spent in the bin, and convolved with the Gaussian (SD = 5 cm); (2) running speed was measured as the time to traverse the constant length bin; (3) head direction was calculated as the mean of the animal's head directions (in degrees) within the bin; (4) lateral position was represented by the mean of the x coordinates of the rat's head within the bin. For each cell, a one-way ANCOVA was performed between different journeys within a given subtask or between different task-demands within a journey as independent (fixed) variables, firing rate as the dependent measure, and with running speed, head direction and lateral position as covariates (R version 2.15.2; The R Project for Statistical Computing, Vienna, Austria).

Spatial similarity

The spatial similarity of the place fields between two trial types or repeated exposures was calculated using the spatial correlation between the maps (Leutgeb et al., 2005). The region of interest (ROI) is defined along a common running route (i.e., the central stem; Figure 4A) between the trial types or sessions.

Rate similarity

Changes in the firing rates between two trial types were expressed by calculating a modified difference/sum score (Leutgeb et al., 2005). The score for a pair of trial types or repeated exposures was obtained by calculating the unsigned difference between the rates in the two trial types or repeated exposures and dividing the difference by the sum of the two rates. To maintain consistency across spatial similarity, the rate similarity was obtained by subtracting the score from one. The possible scores ranged from an asymptotic value of 0–1. The ROI is defined in the same way as spatial similarity.

Neuronal trajectory analysis

The dynamics of the neuronal ensemble activity were analyzed as a trajectory through a state space (neuronal trajectory; Figure 5) (Briggman et al., 2005; Mazor and Laurent, 2005; Harvey et al., 2012). Initially, the path of the rat was linearlized for each journey by projecting the actual trajectory followed by the rat on that lap onto a user-defined idealized path using nearest neighbor Delaunay triangulation. Spatial bins had a resolution of approximately 0.35 cm. At each location on the linearized path, the spike count that convolved with the Gaussian (SD = 5 cm) of the ensemble containing n simultaneously monitored cells was defined as a point in an n-dimensional space, with each dimension representing the activity of a single neuron. Since the neuronal trajectory is in n-dimensional space, to visualize the neuronal trajectory in three-dimensional space, principal component analysis (PCA) was used as a dimension-reduction technique. The data were organized into m locations by an n-cell data set.

The distances between the trajectories were measured as Euclidian distances between the corresponding locations. The classification was performed for individual laps based on the distances to the mean lap trajectories. If the distance to the mean target trajectory was less than the distance to the other mean trajectories, the location for that lap was classified as the targeted choice. Classification accuracy was calculated at each location by averaging across laps for the individual trial types. In all of the distance classifications, the mean trajectories were calculated excluding the test trajectory (i.e., leave-one-out cross validation). In Figure 8, the classification accuracy was measured for each rat based on the most frequently classified choices at the central stem of the maze. The target cells were all simultaneously monitored place cells in each rat. In Figure 9, the classification accuracy for generalization was calculated based on the most frequently classified choices at the central stem of the maze; the combination of target cells for the classification were randomly chosen 1000 times. For control conditions, shuffling was performed by randomizing the sequence of laps 1000 times.

Behavior analysis

The rats could make a decision as soon as the visual cue was detected. However, the exact point where the rats mentally detect the visual cue cannot be estimated because measuring the awareness of the visual cue for rats is difficult to perform. Instead, the timing of the behavioral output was estimated based on the location of turn onset. The turn onset was defined as the last location when the rat was in the central stem or the upper junction of the maze (i.e., DP in Figure 1A) and when the angular velocity was below a threshold value (0.12 rev/s) (Harvey et al., 2012) (Figure 1C). These locations of turn onset were considered the signature of the rats' behavioral decision in the analyses.

Statistical analysis

The spatial and rate similarities, behavior analysis and neuronal trajectory analysis were performed using custom-made programs based on Matlab functions (Matlab 7.14; MathWorks, Natic, MA). All comparisons for the spatial and rate similarities were completed using the Wilcoxon rank sum test. The statistical test for classifier accuracy was performed using the binomial test. The statistical ANCOVA comparison was performed using Matlab functions in conjunction with R (R version 2.15.2; The R Project for Statistical Computing).

Histology

After the task was performed, the rats were sacrificed by pentobarbital sodium overdose and perfused with formalin. The brains were cut coronally at 30 μm and stained with cresyl violet. Each section through the relevant part of the hippocampus was retained for analysis. All dodecatrodes were identified and the tip of each electrode was found by comparing it with the adjacent sections. The recordings from a given dodecatrode were included in the data analysis if its deepest position was under the CA1 pyramidal cell layer.

References

  1. 1
    Hippocampal CA1 place cells encode intended destination on a maze with multiple choice points
    1. JA Ainge
    2. M Tamosiunaite
    3. F Woergoetter
    4. PA Dudchenko
    (2007a)
    J Neurosci 27:9769–9779.
  2. 2
    Exploring the role of context-dependent hippocampal activity in spatial alternation behavior
    1. JA Ainge
    2. MA van der Meer
    3. RF Langston
    4. ER Wood
    (2007b)
    Hippocampus 17:988–1002.
  3. 3
    Heterogeneous modulation of place cell firing by changes in context
    1. MI Anderson
    2. KJ Jeffery
    (2003)
    J Neurosci 23:8827–8835.
  4. 4
    Experience-dependent modifications of hippocampal place cell firing
    1. E Bostock
    2. RU Muller
    3. JL Kubie
    (1991)
    Hippocampus 1:193–205.
  5. 5
    Sequential-context-dependent hippocampal activity is not necessary to learn sequences with repeated elements
    1. MR Bower
    2. DR Euston
    3. BL McNaughton
    (2005)
    J Neurosci 25:1313–1323.
  6. 6
    Optical imaging of neuronal populations during decision-making
    1. KL Briggman
    2. HD Abarbanel
    3. WB Kristan Jr
    (2005)
    Science 307:896–901.
  7. 7
    Sustained activation of hippocampal pyramidal cells by ‘space clamping' in a running wheel
    1. A Czurko
    2. H Hirase
    3. J Csicsvari
    4. G Buzsaki
    (1999)
    Eur J Neurosci 11:344–352.
  8. 8
    The hippocampus, memory, and place cells: is it spatial memory or a memory space?
    1. H Eichenbaum
    2. P Dudchenko
    3. E Wood
    4. M Shapiro
    5. H Tanila
    (1999)
    Neuron 23:209–226.
  9. 9
    Prospective and retrospective memory coding in the hippocampus
    1. J Ferbinteanu
    2. ML Shapiro
    (2003)
    Neuron 40:1227–1239.
  10. 10
    Memory modulates journey-dependent coding in the rat hippocampus
    1. J Ferbinteanu
    2. P Shirvalkar
    3. ML Shapiro
    (2011)
    J Neurosci 31:9135–9146.
  11. 11
    Trajectory encoding in the hippocampus and entorhinal cortex
    1. LM Frank
    2. EN Brown
    3. M Wilson
    (2000)
    Neuron 27:169–178.
  12. 12
    Dynamics of mismatch correction in the hippocampal ensemble code for space: interaction between path integration and environmental cues
    1. KM Gothard
    2. WE Skaggs
    3. BL McNaughton
    (1996)
    J Neurosci 16:8027–8040.
  13. 13
    Accuracy of tetrode spike separation as determined by simultaneous intracellular and extracellular measurements
    1. K Harris
    2. D Henze
    3. J Csicsvari
    4. H Hirase
    5. G Buzsaki
    (2000)
    J Neurophysiol 84:401–414.
  14. 14
    Choice-specific sequences in parietal cortex during a virtual-navigation decision task
    1. CD Harvey
    2. P Coen
    3. DW Tank
    (2012)
    Nature 484:62–68.
  15. 15
    Fast and robust fixed-point algorithms for independent component analysis
    1. A Hyvarinen
    (1999)
    IEEE Trans Neural Netw 10:626–634.
  16. 16
    Abolition of long-term stability of new hippocampal place cell maps by NMDA receptor blockade
    1. C Kentros
    2. E Hargreaves
    3. RD Hawkins
    4. ER Kandel
    5. M Shapiro
    6. RV Muller
    (1998)
    Science 280:2121–2126.
  17. 17
    Independent codes for spatial and episodic memory in hippocampal neuronal ensembles
    1. S Leutgeb
    2. JK Leutgeb
    3. CA Barnes
    4. EI Moser
    5. BL McNaughton
    6. MB Moser
    (2005)
    Science 309:619–623.
  18. 18
    Hippocampal “time cells” bridge the gap in memory for discontiguous events
    1. CJ MacDonald
    2. KQ Lepage
    3. UT Eden
    4. H Eichenbaum
    (2011)
    Neuron 71:737–749.
  19. 19
    A theory of cerebellar cortex
    1. D Marr
    (1969)
    J Physiol 202:437–470.
  20. 20
    Transient dynamics versus fixed points in odor representations by locust antennal lobe projection neurons
    1. O Mazor
    2. G Laurent
    (2005)
    Neuron 48:661–673.
  21. 21
    The Contributions of position, direction, and velocity to single unit-activity in the hippocampus of freely-moving rats
    1. BL McNaughton
    2. CA Barnes
    3. J O'Keefe
    (1983)
    Exp Brain Res 52:41–49.
  22. 22
    Hippocampal synaptic enhancement and information storage within a distributed memory system
    1. BL McNaughton
    2. RGM Morris
    (1987)
    Trends Neurosci 10:408–415.
  23. 23
    Hippocampal synaptic plasticity: role in spatial learning or the automatic recording of attended experience?
    1. RG Morris
    2. U Frey
    (1997)
    Philos Trans R Soc Lond B Biol Sci 352:1489–1503.
  24. 24
    On the directional firing properties of hippocampal place cells
    1. RU Muller
    2. E Bostock
    3. JS Taube
    4. JL Kubie
    (1994)
    J Neurosci 14:7235–7251.
  25. 25
    Spatial firing patterns of hippocampal complex-spike cells in a fixed environment
    1. RU Muller
    2. JL Kubie
    3. JB Ranck Jr
    (1987)
    J Neurosci 7:1935–1950.
  26. 26
    The hippocampus as a spatial map. Preliminary evidence from unit activity in the freely-moving rat
    1. J O'Keefe
    2. J Dostrovsky
    (1971)
    Brain Res 34:171–175.
  27. 27
    The Hippocampus as a Cognitive Map
    1. L O'Keefe
    2. L Nadel
    (1978)
    Oxford: Clarendon.
  28. 28
    Internally generated cell assembly sequences in the rat hippocampus
    1. E Pastalkova
    2. V Itskov
    3. A Amarasingham
    4. G Buzsaki
    (2008)
    Science 321:1322–1327.
  29. 29
    Loss of recent memory after bilateral hippocampal lesions
    1. WB Scoville
    2. B Milner
    (1957)
    J Neurol Neurosurg Psychiatry 20:11–21.
  30. 30
    An information-theoretic approach to deciphering the hippocampal code
    1. WE Skaggs
    2. BL McNaughton
    3. KM Gothard
    4. EJ Markus
    (1993)
    Adv Neural Inf Process Syst 5:1031–1037.
  31. 31
    Learning-related development of context-specific neuronal responses to places and events: the hippocampal role in context processing
    1. DM Smith
    2. SJ Mizumori
    (2006)
    J Neurosci 26:3154–3163.
  32. 32
    Automatic sorting for multi-neuronal activity recorded with tetrodes in the presence of overlapping spikes
    1. S Takahashi
    2. Y Anzai
    3. Y Sakurai
    (2003a)
    J Neurophysiol 89:2245–2258.
  33. 33
    A new approach to spike sorting for multi-neuronal activities recorded with a tetrode—how ICA can be practical
    1. S Takahashi
    2. Y Anzai
    3. Y Sakurai
    (2003b)
    Neurosci Res 46:265–272.
  34. 34
    Real-time and automatic sorting of multi-neuronal activity for sub-millisecond interactions in vivo
    1. S Takahashi
    2. Y Sakurai
    (2005)
    Neuroscience 134:301–315.
  35. 35
    Coding of spatial information by soma and dendrite of pyramidal cells in the hippocampal CA1 of behaving rats
    1. S Takahashi
    2. Y Sakurai
    (2007)
    Eur J Neurosci 26:2033–2045.
  36. 36
    Information in small neuronal ensemble activity in the hippocampal CA1 during delayed non-matching to sample performance in rats
    1. S Takahashi
    2. Y Sakurai
    (2009a)
    BMC Neurosci 10:115.
  37. 37
    Sub-millisecond firing synchrony of closely neighboring pyramidal neurons in hippocampal CA1 of rats during delayed non-matching to sample task
    1. S Takahashi
    2. Y Sakurai
    (2009b)
    Front Neural Circuits 3:9.
  38. 38
    Elements of episodic memory
    1. E Tulving
    (1983)
    New York: Oxford University Press.
  39. 39
    Differential effects of early hippocampal pathology on episodic and semantic memory
    1. F Vargha-Khadem
    2. DG Gadian
    3. KE Watkins
    4. A Connelly
    5. W Van Paesschen
    6. M Mishkin
    (1997)
    Science 277:376–380.
  40. 40
    Spatial and behavioral correlates of hippocampal neuronal activity
    1. SI Wiener
    2. CA Paul
    3. H Eichenbaum
    (1989)
    J Neurosci 9:2737–2763.
  41. 41
    The global record of memory in hippocampal neuronal activity
    1. ER Wood
    2. PA Dudchenko
    3. H Eichenbaum
    (1999)
    Nature 397:613–616.
  42. 42
    Hippocampal neurons encode information about different types of memory episodes occurring in the same location
    1. ER Wood
    2. PA Dudchenko
    3. RJ Robitsek
    4. H Eichenbaum
    (2000)
    Neuron 27:623–633.

Decision letter

  1. Howard Eichenbaum
    Reviewing Editor; Boston University, United States

eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see review process). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.

Thank you for choosing to send your work entitled “Hierarchical organization of context in the hippocampal episodic code” for consideration at eLife. Your article has been evaluated by a Senior editor and 3 reviewers, one of whom is a member of our Board of Reviewing Editors.

The following individuals responsible for the peer review of your submission want to reveal their identity: Howard Eichenbaum (Reviewing editor), Jozsef Csicsvari (peer reviewer), and Emma Wood (peer reviewer).

The Reviewing editor and the other reviewers discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments based on the reviewers' reports.

* The reviewers agree that the paper could be useful and important, but requires additional analyses and improved presentation of the findings. There is a major concern about the need to control for potential differences between location, speed, and head direction in making claims about differential firing on the stem of the maze prior to left and right turn trials (see, for example, Wood et al., 2000). This will require substantial further analyses and likely restriction of the analyses to the central portion of the maze stem.

* There are also concerns about the language and presentation of data. Part of the problem is the use of too much jargon.

* In addition, the description of generalization is difficult to follow. In particular, Figure 7 is hard to interpret, so perhaps the data can be presented in a more intuitive way. The premise is understandable – that knowing the R-L and L-R neuronal trajectories in one task allows classification into R-L and L-R for the other tasks, but the figure does not aid that understanding. Also, it would be useful that provide some analysis and discussion of the relative ability to classify into trajectories (L-R vs R-L) vs into tasks (VD, SA, DA) vs into specific trial types (i.e., task and trajectory information combined).

[Editors' note: the following comments were also sent to the author before acceptance.]

The author attempted to address the concerns of all the reviewers by performing ANCOVAs to eliminate cells for which differences in direction, speed, or lateral position could account for differences in firing patterns associated with distinct trajectories. However, nearly all the examples are taken from areas of the maze stem at the beginning or end of the stem where the animal's behavior, location, and direction of movement clearly differ. This is evident in the nearly all the examples shown in Figure 3 (exceptions are RatF120112-cell 53 and RatG111125-cell 326). We fear that somehow the ANCOVAs are not detecting this confound and this could undermine the conclusions. We recommend that you select a middle region of the maze stem that excludes the top and bottom areas where the trajectory of the rat separates, and use only the data from the middle region for all analyses.

https://doi.org/10.7554/eLife.00321.015

Author response

* The reviewers agree that the paper could be useful and important, but requires additional analyses and improved presentation of the findings.

The entire manuscript has been thoroughly rewritten and reconstructed. According to the reviewers' suggestions, two additional analyses were performed. To improve the presentation of the findings, some comments and figures for explanation were added.

There is a major concern about the need to control for potential differences between location, speed, and head direction in making claims about differential firing on the stem of the maze prior to left and right turn trials (see, for example, Wood et al., 2000). This will require substantial further analyses and likely restriction of the analyses to the central portion of the maze stem.

In the revised manuscript, I performed an additional analysis using ANCOVA as Wood et al., 2000 had conducted. The results are shown in Figure 4 of the revised manuscript and suggest that the differential firing in the central stem of the maze was significant even when running speed, head direction, and lateral position are taken into account. This is consistent with the findings of Wood et al., 2000. I added a paragraph describing this to the Results section and a few sentences to the Discussion section. To describe my methods in detail, I added a paragraph to the Materials and methods section.

* There are also concerns about the language and presentation of data. Part of the problem is the use of too much jargon.

In line with the reviewers' suggestions, I have thoroughly rewritten and reconstructed the previous manuscript. In particular, I have replaced the terms “spatiotemporal context”, “non-spatial context”, and “combination of spatiotemporal and non-spatial contexts” with “journey”, “task-demand”, and “trial type”, respectively, throughout the Results, Discussion, and Materials and methods sections, and the figure legends in the revised version of the manuscript. In addition, I have replaced “spatial measure” with “spatial correlation”.

* In addition, the description of generalization is difficult to follow. In particular, Figure 7 is hard to interpret, so perhaps the data can be presented in a more intuitive way. The premise is understandable – that knowing the R-L and L-R neuronal trajectories in one task allows classification into R-L and L-R for the other tasks, but the figure does not aid that understanding.

To clearly demonstrate the presentation of the data for generalization accuracy, I have added a subfigure to Figure 9 (Figure 7 in the previous version of the manuscript) and some comments to the legend of Figure 9. In addition, to clearly describe the generalization analysis, I have added a few sentences to the Results section.

Also, it would be useful that provide some analysis and discussion of the relative ability to classify into trajectories (L-R vs R-L) versus into tasks (VD, SA, DA) versus into specific trial types (i.e., task and trajectory information combined).

In line with the reviewers' suggestions, I have performed an additional analysis to assess the relative ability to classify journeys, versus task-demands, versus specific trial types using a binary classification scheme. Figure 8 of the revised manuscript shows these results. The results suggest that the abilities of the hippocampal ensemble activity to identify journeys, task-demands, and trial types are sufficiently high and relatively similar to each other. I added details pertaining to my methods to the Materials and methods section. Moreover, I have added some text addressing these results to the Results and Discussion sections.

[Editors' note: the following comments were also sent to the author before acceptance.]

The author attempted to address the concerns of all the reviewers by performing ANCOVAs to eliminate cells for which differences in direction, speed, or lateral position could account for differences in firing patterns associated with distinct trajectories. However, nearly all the examples are taken from areas of the maze stem at the beginning or end of the stem where the animal's behavior, location, and direction of movement clearly differ. This is evident in the nearly all the examples shown in Figure 3 (exceptions are RatF120112-cell 53 and RatG111125-cell 326). We fear that somehow the ANCOVAs are not detecting this confound and this could undermine the conclusions.

The reviewers pointed out that the detection error of the ANCOVAs could undermine the main conclusions because many place fields in Figure 3 in the previous manuscript can be seen at the beginning or end of the stem. I have addressed this concern in three ways.

1) Results of the reanalysis: I completely agree with the reviewers' concern about Figure 3. However, although I repeatedly checked place fields in Figure 3 in the previous manuscript, the results suggested that they fired at significantly different rates in the central stem even when head direction, running speed, and lateral position in the central stem were taken into account. For example, although RAT# G111125 Cell 109 in Figure 3b in the previous manuscript showed a large part of the place fields at the beginning or end of the stem (dotted red box in Author response image 1), all analyses in the previous manuscript excluded such regions. Instead, only other remaining parts of the place fields in the central stem (i.e., locations where the rat's head direction, running speed, and lateral position do not clearly differ) were analyzed in the previous manuscript (solid red box in Author response image 1; arrows: analyzed place fields) and showed significant rate differences. Therefore, I concluded that there is no serious concern about the detection error of the ANCOVAs.

Author response image 1
Red solid and dotted boxes indicate the analyzed and unanalyzed regions of the central stem, respectively.

Arrows indicate the location of the analyzed place fields. The rate was coded on a color scale from blue (silent) to red (maximum rate). The pixels that were not sampled are white. The symbols above the map indicate the journey, task-demand and illuminated light-cue. For each cell, the rate scale corresponds to the peak firing rate in Hz in that condition (indicated to the right of the rate maps).

2) The bias of the place fields location shown in Figure 3: As the journey-dependent coding was changed in a global remapping manner, if the place fields of a place cell in one of two journeys were located at the beginning or end of the central stem, the fields in another journey, if any, tended to be located at the another side (e.g., Author response Image 2, RAT# G111125 Cell 109 in Figure 3b in the previous manuscript). The length of the central stem was 60cm and the place field width was approximately 20–30cm so that the beginning and end of the central stem tended to be exclusive to each other. In addition to the tendency, to demonstrate the place fields in the central stem in both left-to-right and right-to-left journeys, and the main finding (i.e., the journey-dependent coding was changed in a global remapping manner), the place fields in Figure 3 in the previous manuscript were selected. This could have caused the reviewers' concern.

Author response image 2
The journey-dependent global remapping may make the place fields be exclusively located at the beginning or end of the central stem.

The rate was coded on a color scale from blue (silent) to red (maximum rate). The pixels that were not sampled are white. The symbols above the map indicate the journey, task-demand, and illuminated light-cue. For each cell, the rate scale corresponds to the peak firing rate in Hz in that condition (indicated to the right of the rate maps).

We recommend that you select a middle region of the maze stem that excludes the top and bottom areas where the trajectory of the rat separates, and use only the data from the middle region for all analyses.

3) Analyses of the place fields in the middle region of the central stem: In line with the reviewers' suggestion, I selected cells which showed the place fields restricted in a middle region of the central stem and fired at significantly different rates between different journeys or between different subtasks even when running speed, head direction, and lateral position were taken into account. The results shown in Author response images 4–9 did not undermine the main conclusions, with the exceptions that only the distributions of the change in firing rates in the place fields between VD3 and SA2 were similar to control condition (n.s., P=0.44, Wilcoxon rank sum test, Author response image 5), and that the relative ability of distinguishing trial types was slightly higher than that of distinguishing task-demands (P=0.03, Kruskal-Wallis test with Tukey honestly significant difference (HSD) post-hoc test, Author response image 8). However, as the number of selected cells for each ensemble (K: 3 cells, F: 4 cells, L: 9 cells, B: 16 cells, G: 27 cells; Only 5% (59/1119) of the recorded cells) was too few to statistically analyze at the ensemble level, these minor exceptions might be due to the sampling bias.

Taken together, as I concluded that the reason that nearly all the examples in Figure 3b are taken from areas of the maze stem at the beginning or end of the stem was not due to the detection error of the ANCOVAs, but due to the selection bias influenced by the global remapping of place fields, I would instead fear that the restricted selection of cells in the middle of the central stem causes the sampling bias and could undermine the main conclusions. Therefore, the data in this response letter were not included in the revised manuscript.

Also, in all examples, indicate that middle region on illustrations show only examples of cells that differ in firing pattern in that middle section.

As I totally agree with the reviewers' concern about Figure 3b, I removed ambiguous cells and replaced five cells in Figure 3b in the previous manuscript with cells that fired at significantly different rates in the middle of the central stem in the revised manuscript (see Author response image 3; boxes indicate the middle of the central stem).

Author response image 3
Journey-dependent coding of a place-specific activity modulated with task-demand.

The rate maps for five CA1 pyramidal cells. The maps are averages for repeated sessions. Each row shows the data from one pyramidal cell. Note that different journeys in the same place resulted in place fields similar in locations, but different task-demands caused different firing intensities in those place fields. Boxes indicate the middle region of the central stem. The rate was coded on a color scale from blue (silent) to red (maximum rate). The pixels that were not sampled are white. The symbols above the map indicate the journey, task-demand, and illuminated light-cue. For each cell, the rate scale corresponds to the peak firing rate in Hz in that condition (indicated to the right of the rate maps).

Author response image 4
Proportion of place cells that fired at significantly different rates in the middle region of the central stem.
Author response image 5
A quantitative assessment of the spatial and rate similarities between a pair of trial types for place cells that fired at significantly different rates in the middle region of the central stem.

(a) The region of interest (ROI) was set at the region highlighted in the green shaded box. (c) In contrast to figure 4, only the distributions of the change in firing rates in the place fields between VD3 and SA2 were similar to control condition (n.s., P=0.44, Wilcoxon rank sum test).

Author response image 6
Neuronal trajectory of different journeys and different future directions for place cells that fired at significantly different rates in the middle region of the central stem.
Author response image 7
Episode-specific neuronal trajectories for place cells that fired at significantly different rates in the middle region of the central stem.
Author response image 8
Comparison of accuracy of distinguishing journeys, task-demands, and trial types from neuronal ensemble activity of place cells that fired at significantly different rates in the middle region of the central stem.

In contrast to figure 8, the average classification accuracy of distinguishing trial types was different from that of distinguishing task-demands (P=0.03, Kruskal-Wallis test with Tukey honestly significant difference (HSD) post-hoc test).

Author response image 9
Generalization accuracy of determining a journey over task-demands from neuronal ensemble activity of place cells that fired at significantly different rates in the middle region of the central stem.
https://doi.org/10.7554/eLife.00321.016

Article and author information

Author details

  1. Susumu Takahashi

    1. Laboratory of Neural Circuitry, Graduate School of Brain Science, Doshisha University, Kizugawa, Japan
    2. Precursory Research for Embryonic Science and Technology (PRESTO), Japan Science and Technology Agency, Kawaguchi, Japan
    3. Faculty of Computer Science and Engineering, Kyoto Sangyo University, Kyoto, Japan
    Contribution
    ST, Designed the tasks, set up the behaviour and electrophysiology apparatus, performed the surgery, behaviour and electrophysiology experiments, data collection and analyses, and wrote the paper.
    For correspondence
    stakahas@mail.doshisha.ac.jp
    Competing interests
    The author declares that no competing interests exist

Funding

JST PRESTO

  • Susumu Takahashi

JSPS KAKENHI (23650158, 24300148)

  • Susumu Takahashi

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

I would like to thank Y. Sakurai for the initial design of the behavioral task and for discussion and comments on the manuscript, all the members of the Sakurai and Fujiyama laboratories for technical assistance, and M. Kawato and all of the advisors and members of the JST PRESTO program ‘Decoding and controlling brain information' for their fruitful comments and encouragement.

Ethics

Animal experimentation: All procedures were approved by the Doshisha University (approved number: 1229) and Kyoto Sangyo University (approved numbers: 2011-04 and 2011-05) Institutional Animal Care and Use Committees.

Reviewing Editor

  1. Howard Eichenbaum, Reviewing Editor, Boston University, United States

Publication history

  1. Received: October 12, 2012
  2. Accepted: December 26, 2012
  3. Version of Record published: February 5, 2013 (version 1)

Copyright

© 2013, Takahashi

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

  • 1,007
    Page views
  • 151
    Downloads
  • 9
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.

Comments

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Download citations (links to download the citations from this article in formats compatible with various reference manager tools)

Open citations (links to open the citations from this article in various online reference manager services)

Further reading

    1. Neuroscience
    Ting-hao Huang et al.
    Tools and Resources Updated