Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE

User menu

  • Log out
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log out
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE
PreviousNext
Research Articles, Systems/Circuits

Activity Subspaces in Medial Prefrontal Cortex Distinguish States of the World

Silvia Maggi and Mark D. Humphries
Journal of Neuroscience 18 May 2022, 42 (20) 4131-4146; https://doi.org/10.1523/JNEUROSCI.1412-21.2022
Silvia Maggi
School of Psychology, University of Nottingham, Nottingham NG7 2RD, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mark D. Humphries
School of Psychology, University of Nottingham, Nottingham NG7 2RD, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Mark D. Humphries
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Medial prefrontal cortex (mPfC) activity represents information about the state of the world, including present behavior, such as decisions, and the immediate past, such as short-term memory. Unknown is whether information about different states of the world are represented in the same mPfC neural population and, if so, how they are kept distinct. To address this, we analyze here mPfC population activity of male rats learning rules in a Y-maze, with self-initiated choice trials to an arm end followed by a self-paced return during the intertrial interval (ITI). We find that trial and ITI population activity from the same population fall into different low-dimensional subspaces. These subspaces encode different states of the world: multiple features of the task can be decoded from both trial and ITI activity, but the decoding axes for the same feature are roughly orthogonal between the two task phases, and the decodings are predominantly of features of the present during the trial but features of the preceding trial during the ITI. These subspace distinctions are carried forward into sleep, where population activity is preferentially reactivated in post-training sleep but differently for activity from the trial and ITI subspaces. Our results suggest that the problem of interference when representing different states of the world is solved in mPfC by population activity occupying different subspaces for the world states, which can be independently decoded by downstream targets and independently addressed by upstream inputs.

SIGNIFICANCE STATEMENT Activity in the medial prefrontal cortex plays a role in representing the current and past states of the world. We show that during a maze task, the activity of a single population in medial prefrontal cortex represents at least two different states of the world. These representations were sequential and sufficiently distinct that a downstream population could separately read out either state from that activity. Moreover, the activity representing different states is differently reactivated in sleep. Different world states can thus be represented in the same medial prefrontal cortex population but in such a way that prevents potentially catastrophic interference between them.

  • decoding
  • learning
  • medial prefrontal cortex
  • population activity
  • sleep reactivation

Introduction

The medial prefrontal cortex (mPfC) plays key roles in adaptive behavior, including reshaping behavior in response to changes in a dynamic environment (Euston et al., 2012) and in response to errors in performance (Narayanan and Laubach, 2008; Laubach et al., 2015). Damage to mPfC prevents shifting behavioral strategies when the environment changes (Laskowski et al., 2016; Guise and Shapiro, 2017). Single neurons in mPfC shift the timing of spikes relative to hippocampal theta rhythms just before acquiring a new action–outcome rule (Benchenane et al., 2010). And multiple labs have reported that global shifts in mPfC population activity precede switching between behavioral strategies (Rich and Shapiro, 2009; Durstewitz et al., 2010; Karlsson et al., 2012; Powell and Redish, 2016) and the extinction of learned associations (Russo et al., 2021).

Adapting behavior depends on knowledge of both the past and the present state of the world. Deep lines of research have established that mPfC activity represents information about both. The memory of the immediate past is maintained by mPfC activity in tasks requiring explicit use of working memory (Baeg et al., 2003; Fujisawa et al., 2008; Spellman et al., 2015). The use of such memory is seen in both the impairment arising from mPfC lesions (Rich and Shapiro, 2007; Young and Shapiro, 2009; Laskowski et al., 2016) and the role of mPfC in error monitoring (Laubach et al., 2015). Representations of stimuli and features happening in the present have been reported in a variety of decision-making tasks throughout PfC (Averbeck et al., 2006; Rigotti et al., 2013; Hanks et al., 2015; Siegel et al., 2015) and specifically within rodent mPfC (Sul et al., 2010; Ito et al., 2015; Guise and Shapiro, 2017).

Little is known though about how mPfC activity represents multiple states of the world. Prior studies have shown that past and upcoming choices can both modulate activity of neurons in the same mPfC population (Baeg et al., 2003; Ito et al., 2015), but none have compared how different states of the world are represented. Thus important questions remain regarding if and how different world states are encoded in the same mPfC population and how those representations are kept distinct.

To address these questions, we reanalyze here mPfC population activity from rats learning new rules on a Y-maze (Peyrache et al., 2009). This task had distinct trial and intertrial interval phases, and we have previously shown that task features of the preceding trial can be decoded from the population activity in the intertrial interval (Maggi et al., 2018), showing that mPfC activity in this task depends on the state of the world. We can thus address our key questions here by asking whether population activity in the trials also represents the state of the same task features and, if so, how that representation is kept distinct between the trial and intertrial interval phases.

We find that the mPfC population activity occupies different subspaces between trials and intertrial intervals, providing a basis for separately representing at least these two distinct states of the world. Consistent with representing world states, task features could be decoded from activity in both the trial and intertrial interval phases but were strongly distinct: decoding was of the present features in the trial and predominately of features of the preceding trial during the intertrial interval. Decoding axes were, or close to, orthogonal between the trials and intertrial intervals, showing that the subspaces supported distinct encodings. Further, population activity of the trials and intertrial intervals preferentially reactivated in post-training sleep in different ways: preferential reactivation of trial activity uniquely occurred after learning and correlated with performance during training. Our results thus suggest that representing different world states using independently decodable axes within a mPfC population could prevent interference between them, allowing them to be separately accessed by both downstream and upstream populations.

Materials and Methods

Task description and electrophysiological data

All data in this study come from previously published data (Peyrache et al., 2009). Full details of training, spike sorting, and histology can be found in Peyrache et al. (2009). The experiments were conducted in accordance with institutional (Centre National de la Recherche Scientific Comité Opérationnel d'Éthique dans les Sciences de la Vie) and international (National Institutes of Health guidelines) standards and legal regulations (Certificate no. 7186, French Ministère de l'Agriculture et de la Pêche) regarding the use and care of animals.

Four male Long-Evans rats were implanted with tetrodes in the medial wall of prefrontal cortex, covering the prelimbic and infralimbic regions, and trained on a Y-maze task (Fig. 1a). During each session, neural activity was recorded for 20–30 min of sleep or rest epoch before the training epoch, in which rats worked at the task for 20–40 min. After that, another 20–30 min of sleep or rest epoch recording followed. During the sleep epochs, intervals of slow-wave sleep were identified off-line from the local field potential (Peyrache et al., 2009; Benchenane et al., 2010).

The Y-maze had symmetrical arms, 85 cm long, 8 cm wide, and separated by 120°, connected to a central circular platform (denoted as the choice point throughout). The two choice arms had a light at the end, one of which was lit during each trial in a pseudorandom sequence. Rats self-initiated a trial by leaving the beginning of the start arm. A trial finished when the rat reached the end of the chosen goal arm. If the chosen arm was correct according to the current rule, the rat was rewarded with drops of flavored milk. As soon as the animal reached the end of the chosen arm, an intertrial interval started and lasted until the rat completed its self-paced return to the beginning of the start arm. The central platform was raised once the rat passed it to prevent backtracking along the choice arms. The light was extinguished during the return journey; unfortunately from the data available to us it is not clear exactly when (F. Battaglia, personal communication).

Each rat was exposed to the task completely naive and had to learn each rule by trial and error. The rules were presented in the following sequence: Go to the right arm, go to the lit arm, go to the left arm, go to the dark arm. A rule was switched to the next in the sequence when the animal had achieved 10 correct trials in a row, or 11 of 12. Across the four rats, there were eight rule switches in total.

The recording sessions taken from the study of Peyrache et al. (2009) were 53 in total. Each of the four rats learned at least two rules, and they respectively contributed to 14, 14, 11, and 14 sessions. We used 49 of these sessions for our analysis, of between 7 and 51 trials each. One session was omitted for missing position data, one for consistent choice of the right arm (in a dark arm rule) preventing decoder analyses (see below), and one for missing spike data in a few trials. An additional session was excluded for having only two neurons firing in all trials. Tetrode recordings were spike sorted within each recording session. Spikes were recorded with a resolution of 0.1 ms. Simultaneous tracking of the position of the rat was recorded at 30 Hz.

Testing for separable population activity

We evaluated the difference between population activity in the trial and intertrial intervals of a session by quantifying their separability in a low-dimensional space. For consistency with our previous work, for each session we selected the N active neurons that fired at least one spike on each trial (Fig. 1e), allowing us to directly compare the decoding results obtained here (see below) with those in Maggi et al. (2018); the populations thus ranged between 4 and 22 neurons (Fig. 1i).

We used principal component analysis (PCA) to project the population vectors of a session onto a common set of dimensions. For each session, we constructed an N-length vector of neuron firing rates in each trial rt, resulting in the set of population firing rate vectors {rt(1),...,rt(T)} across the T trials of a session. We similarly constructed the set of vectors of neuron firing in each intertrial interval {rI(1),...,rI(T)}. We then constructed the data matrix X from the firing rate vectors of the population, by concatenating trials and intertrial intervals in their temporal order {rt(1),rI(1),...,rt(T),rI(T)}′ across the T trials of a session; the resulting matrix thus had dimensions of 2T rows and N (neurons) columns. Applying PCA to X, we projected the firing rate vectors on to the top d principal axes (the eigenvectors of X′X) to create the top d principal components. For each set of d components, we quantified the separation between the projected trial and intertrial interval population vectors using a linear classifier (support vector machine), and report the proportion of misclassified vectors. We repeated this for between d = 1 and d = 4 axes for each session.

Linear decoding of task features

To predict which task feature was encoded in mPfC population activity, we trained and tested a range of linear decoders (Hastie et al., 2009). Here, we report the results obtained using a logistic regression classifier, but for robustness we also tested three other decoders—linear discriminant analysis, linear support vector machines, and a nearest neighbors classifier—and found similar results. The full details of the decoding analysis can be found in Maggi et al. (2018).

Task information of each trial was binary labeled for three features, outcome (labels 0, 1), the direction of the chosen arm (labels left, right), and the arm position of the light cue (labels left, right). We used leave-one-out cross-validation to decode each feature from population activity, holding out the vector of the ith trial, rt(i), training the classifier on the N – 1 remaining trial vectors, and then using the resulting weight vector to predict the label of the feature for the held-out trial. We quantified the accuracy of the decoder as the proportion of correctly predicted labels over all T held-out trials of a session. The same approach was used for the intertrial intervals, using the set of firing rate vectors {rI(1),...,rI(T)} across the T intertrial intervals.

For decoding at different positions in the maze, we first linearized the maze, divided it into five equally sized sections, and then computed the N-length firing rate vector of the population for each position p, rpt in the trial. For each trial t = 1,…,T of a session and each section of the maze p = 1,…,5, the set of population firing rate vectors {rpt(1),...,rpt(T)} was used to train the cross-validated decoder. The same approach was used for the position-dependent vectors {rpI(1),...,rpI(T)} in the intertrial intervals of a session.

For each rat and each session, the distribution of outcomes and arm choices depended on the performance of the rats, which could differ from 50%. Therefore, we trained and cross-validated the same classifier on the same datasets but shuffling the labels of the task features across trials. In this way we obtained the accuracy of detecting the correct labels by chance. We repeated the shuffling and fitting 50 times, and we averaged the accuracy across the 50 repetitions.

Testing for independent decoding

To compare the decoding axes between the trials and intertrial intervals, we again trained the classifier separately for each of the three task features but now using all the population firing rate vectors of a session, first for the trials {rt(1),...,rt(T)} and then for the intertrial intervals {rI(1),...,rI(T)}. For a given feature f, we then computed the angle, θf, between the resulting vectors of decoding weights for the trials, wt(f), and intertrial intervals, wI(f), as follows: θ=cos−1(wt(f)·wI(f)||wt(f)||||wI(f)||). Similarly, to assess how the different features were simultaneously encoded during a task phase, we computed the angles between the decoding vectors of two features f1 and f2 within a trial, or within an intertrial interval.

We evaluated the degree of independence between the trial and intertrial interval decoding by attempting to cross-decode a task feature in one phase from the activity in the other. For a given task feature, we took the above classifier trained on all trials of a session and tested its decoding on all intertrial intervals of the same session, with performance reported as the percentage of correctly labeled intertrial intervals. We also tested this in reverse, decoding the feature in the trials from the classifier trained to decode the same feature from all the intertrial intervals. To check we were not overfitting when using a decoder trained on all T phases, we further tested cross-decoding using leave-one-out by leaving out the ith trial interval pair, training on N – 1 trials, and predicting the ith intertrial interval (and vice versa for training on intertrial intervals and testing on trials). Performance for leave-one-out cross-decoding was reported as the percentage of correctly labeled held-out trials (or intertrial intervals) over all trials of the session.

Behavioral analysis

To check whether our decoding results depended on potentially different behaviors or task demands, we divided the sessions in two different ways, by rule type and by learning type. For the rule type, we grouped sessions by whether the target rule was a direction-based rule (so putatively egocentric) or a cue-based rule (so putatively allocentric).

To group by learning type, we identified learning sessions according to the criteria of the original study (Peyrache et al., 2009) of a session with three consecutive correct trials followed by a performance of at least 80% correct. The first of the three correct trials was the learning trial. Only 10 sessions satisfied these criteria. All sessions that did not meet these criteria were labeled Other.

We quantified performance in learning sessions by fitting a piecewise linear regression model to the cumulative reward curve, using robust regression to fit lines before and after the learning trial. The slopes of the two lines gave us the rate of reward accumulation before (rbefore) and after (rafter) the learning trial (Fig. 1b). We quantified performance on all other sessions in a similar way (Fig. 1c,d). For the eight rule-change sessions, we considered the slopes of the regression lines before and after the rule-switch trial. For all remaining sessions, we looked for any performance change by fitting the piecewise linear regression model to each trial in turn (allowing a minimum of five trials before and after each tested trial). We then found the trial at which the increase in slope (rafter − rbefore) was maximized, indicating the point of steepest inflection in the cumulative reward curve.

Reactivation of task-feature representations in sleep

To quantify the reactivation of waking activity in pre-training and post-training sleep, we used the population firing rate vectors computed for the decoder {r(1),...,r(T)}. We considered here the average population vector for a feature in each session, computed across all the trials for each feature. For example, we quantified the average population firing rate vector for all the right-choice trials and separately for all the left-choice trials. This vector thus represented the region in the activity subspace (Fig. 2) occupied by that particular feature in the trial or the intertrial interval.

We then compared this feature-specific activity vector with the firing rate vector of each 1 s time bin of slow-wave sleep pre-training and post-training using Spearman's correlation coefficient. This gave us a distribution of correlations between the feature-specific vector and the population activity vectors during pre-training slow-wave sleep and a distribution of correlations between the feature-specific vector and the population activity vectors during post-training slow-wave sleep.

Spearman's coefficient was chosen specifically to compare the relative activity of the neurons in the population between training and sleep epochs, and so we call this “reactivation,” not “replay.” Replay implies that specific patterns of firing from waking, such as sequences of place cells (Skaggs and McNaughton, 1996; Lee and Wilson, 2002; O'Neill et al., 2008; Denovellis et al., 2021) or sequences of neurons in an ensemble (Euston et al., 2007; Peyrache et al., 2009, 2010), reappear during sleep or quiescence. As we use it here, reactivation is assessing how well the sleep activity aligns with the two subspaces of trial and intertrial interval activity during training, and consequently suggests whether those two are being revisited.

If a feature-specific activity vector was preferentially reactivated in post-training sleep, then we would expect the distribution of the correlation coefficients between a feature and post-training slow-wave sleep to be right shifted compared with the distribution of the correlation coefficients between the same feature and pre-training slow-wave sleep. We quantified this shift by measuring the difference in the medians (Mpost − Mpre) between the two distributions of correlation coefficients. If the difference was positive, then we had a higher correlation of the feature-specific vector with the activity in post-training slow-wave sleep than with the activity in pre-training slow-wave sleep. If negative, then the feature-specific vector was more similar to the pre-training slow-wave sleep population activity.

To control for different time scales of reactivation in sleep we repeated the same procedure, changing the time bin in the slow-wave sleep pre-training and post-training. Bin sizes from 100 ms to 10 s were chosen to range below and above the mean length of a trial (∼6.5 s).

Statistics

Quoted measurement values are mean x¯ and SEM. Differences between two paired distributions were assessed using a paired Wilcoxon signed-rank test; differences from zero were assessed using a Wilcoxon signed-rank test. In Figures 3-7 we report where p values for these tests exceeded three alpha levels (0.05, 0.01, 0.005). Differences between distributions were assessed using the Kolomogorov–Smirnov test. Throughout, we have n = 49 sessions; in some analyses we subdivide these into rule types (n = 15 direction-rule sessions and n = 34 cue-rule sessions) or learning types (n = 10 learning sessions and n = 39 Other sessions). In Figure 3 we break down the decoding results by each rat, giving n = 14, 14, 11, 14 sessions for rats 1–4, respectively.

Data availability

The spike-train and behavioral data that support the findings of this study are available at http://crcns.org/data-sets/pfc/pfc-6.

Results

We analyze here data from rats learning rules in a Y-maze that had tetrodes implanted in mPfC before the first session of training. Across sessions, animals were asked to learn one of four rules that were given in sequence (Go to the right arm, go to the lit arm, go to the left arm, go to the dark arm). Rules were switched after 10 correct choices (or 11 of 12). The animal self-initiated each trial by running along the central stem of the Y-maze and choosing one of the arms (Fig. 1a). The trial finished at the end of the arm, and a reward was delivered if the chosen arm matched the current rule being acquired. During the subsequent intertrial interval, the rat made a self-paced return to the start of the central arm to initiate the next trial. Trials were 6.5 ± 0.5 s on average; intertrial intervals were 55.6 ± 1.1 s. Throughout, population activity was recorded in the prelimbic and infralimbic cortex (Fig. 1e), for which we shall use the term mPfC here (Laubach et al., 2018, propose that these regions are equivalent to the anterior cingulate cortex in primates).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Rule learning and neural activity on the Y-maze task. a, Schematic of the Y-maze task, showing a rat at the start position. A trial is the period from the start position to the end of the chosen arm; the intertrial interval is the return from the arm end to the start position. On each trial one arm end was lit, chosen in a pseudorandom order, regardless of whether it was relevant to the currently enforced rule. Across sessions, animals were asked to learn one of the following four rules in the sequence: Go to the right arm, go to the lit arm, go to the left arm, go to the dark arm. Rules switched after 10 correct choices (or 11 of 12). b, Example reward curve from a learning session, plotting the cumulative number of correct trials. Black dashed line identifies the learning trial as the first of three consecutive correct trials followed by at least 80% correct trials. Inset, Reward rates before (light red) and after (dark red) the learning trial. Reward rates were given by the slope of linear regressions fitted to the reward curve before and after the learning trial. c, Example reward curve from an Other session, one of eight in which the rule switched; the black dashed line identifies the rule change trial. Inset, Reward rates before and after the rule change. d, Change in reward rate during all learning sessions (red) or Other sessions (blue). e, Raster plots of spiking activity in the medial prefrontal cortex during a single trial and the following intertrial interval. f, Neuron firing rates in each trial and in the following intertrial interval. Left, An example session, black dots are the data in e. Right, All sessions. g, For the same example session, the distribution of Spearman's rank coefficients between the population vectors of firing rates in the trial and intertrial interval (left) and the corresponding p values (right); p = 0.05 is indicated by the dashed line. h, As in g for all sessions. i, The number of neurons in populations analyzed here, by rat. Each symbol is a session. Boxes plot median and interquartile range (IQR); whiskers extend to 1.5 IQR.

Neural activity statistics differed between trials and intertrial intervals. Neurons that had the highest firing rates in a trial tended to have lower firing rates in the following intertrial interval (Fig. 1f). However, the vector of rates across the population was not strongly correlated between the trial and following intertrial interval (Fig. 1g,h), meaning that changes in firing rates between the two phases of the task were not systematically in one direction. This low correlation of population activity between the trial and intertrial interval is also consistent with a change in representation as we now report.

Separable subspaces of population activity between trials and intertrial intervals

We first asked whether population activity occupied different subspaces between consecutive pairs of trials and intertrial intervals as a basis for representing these as two different states of the world in the same population. To do so, we projected all population activity vectors of a session (Fig. 2a) into a low-dimensional space (Fig. 2b) and then quantified how easily we could separate them into trials and intertrial intervals. Using just one dimension for the projection was sufficient for near-perfect separation in many sessions; using two was sufficient for above-chance performance in all sessions (Fig. 2c). Population activity thus occupied a different low-dimensional subspace between the trials and the intertrial intervals.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Population activity is easily separable between trials and intertrial intervals. a, Population activity vectors for the trials (black dot) and following intertrial intervals of one session with 54 trial-ITI pairs. The heat map shows the normalized firing rate for each neuron. b, Projection of population activity vectors of that session on to two dimensions shows a complete separation of trial and intertrial interval activity. The black line is the linear separation found by the classifier. PC, Principal component. c, Summary of classification error over all sessions as a function of the number of dimensions. Each gray dot is the error for one session at that number of projecting dimensions. Dashed line gives chance performance. Box plots show medians (red line), interquartile ranges (blue box), and outliers (red pluses). d–f, Same as a–c but comparing population activity vectors for maze sections 4 and 5 in the intertrial interval. g, Schematic of maze sections. h, The average classification error in the separation of population activity vectors between each pair of maze sections within trials. Classification error is for projections of the population activity vectors in a two-dimensional space. i, As for h, within intertrial intervals. j, As for h, for the separability of the trial and intertrial interval population activity across maze sections. The white square indicates the arm-end position, where the transition from trial to intertrial interval occurs. k, Distribution of the delay between the start of the intertrial interval and the start of the outbound phase, with a median delay of 1 s (dashed line) and a mean of 1.14 s. l, Same session as a but divided by inbound and outbound phases. m, As for b, classifying vectors by inbound and outbound phases. n, As for c, for inbound and outbound vectors. o, Comparing the classification error for the population vectors divided by trial-ITI or by inbound-outbound phases. Each line is a session; black lines show medians. The p values are from paired Wilcoxon signed-rank tests.

This was not true within each phase (Fig. 2d–f): when we divided the maze into sections (Fig. 2g), the population activity at nearby maze positions was not easily separable within trials (Fig. 2h) or within intertrial intervals (Fig. 2i), with the notable exception of position 1—the return to the starting position—in the intertrial intervals. By contrast, population activity vectors at one position in the trial and another in the intertrial interval were easily separable between every pair of positions (Fig. 2j). Thus, although of course population activity changed across maze positions (Fig. 2h,i), those changes were smaller and continuous within each phase of the task but larger and discontinuous between them as they moved into a different subspace of activity.

A detailed examination of what might be driving this move into a different subspace of population activity is beyond our scope here, but we can show there are at least two plausible causes. Aligning population activity to the trial and intertrial interval implies this change is caused by reaching the arm end. But a range of other salient events may be causing this shift in population activity. To begin examining possible events, we instead divided the task into inbound and outbound phases, where the start of the outbound phase was defined from the point where the heading direction of the animal had turned toward the start of the maze for at least 400 ms, which happened on average 1.14 s later than the start of the intertrial interval (Fig. 2k). Dividing the population activity accordingly (Fig. 2l), we found the separability of inbound and outbound phase population vectors was excellent even in only a few dimensions (Fig. 2m,n).

Indeed, population separability was equally good for both the trial-ITI and inbound-outbound separations (Fig. 2o). It is thus equally plausible that the shift of subspace occupied by population activity is driven by a change in heading direction as by reaching the arm end. As these events occur close together in time (Fig. 2k), considerable further work, and likely additional experiments, would be necessary to tease apart the causal mechanisms (see Discussion). For the remainder of the article we thus continue to consider activity subspaces defined by the trial-ITI split while being mindful that exactly when the shift occurred and what causes it is unknown.

Different states of the same task features can be decoded from population activity

We then tested whether these distinct subspaces corresponded to encoding different states of the world between the trials and intertrial intervals. Using a linear decoder on the vector of population activity during each trial (see above, Materials and Methods), we decoded key features of the task during the trial, that is, the choice of arm direction of the animal in the trial, the outcome of the trial, and which arm end was lit during the trial. We trained the same decoders using the same population vectors but with features shuffled across trials (see above, Materials and Methods) to define appropriate chance levels for each decoder given the unbalanced distribution of some task features such as outcome.

We could decode all direction choice, light position, and outcome in the current trial above chance (Fig. 3a,b, left). In Figure 3a we plot the absolute accuracy of decoding; in Figure 3b we also plot the decoding accuracy relative to the shuffled data for each session, which, as it accounts for the different distributions of features (e.g., outcome) in each session, better shows the effect size of the decoding. Relative decoding accuracies well above zero could even be seen for each animal (Fig. 3d), despite the small populations (Fig. 1i, median 10 neurons), the limited numbers of trials per session (median 29 trials) available for training the decoder, and the low number (11–14) of sessions. As the outcome was not yet known during the trial, the ability to decode outcome implies anticipatory activity for outcome in mPfC neurons, as previously reported by a number of labs (Euston et al., 2012). However, we found no correlation between the performance of an animal in a session and our ability to decode the upcoming outcome (Fig. 3c), suggesting this anticipatory activity is not dependent on how frequently reward was acquired. Nonetheless, although it is unclear what this anticipatory activity reflects, the ability to decode outcome was robust across all classifiers we tested (data not shown; see above, Materials and Methods; Maggi et al., 2018). To test for the effects of past states on population activity in the trials, we also tried decoding the direction choice, outcome, and light position of the preceding trial and found that decoding was at or close to chance (Fig. 3a,b, right; e, each subject). Population activity in mPfC during the trials thus depended on features in the present state of the task and weakly or not at all on features in the past trials.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Decoding of different task states in the trials and intertrial intervals. a, Accuracy of decoding task features from population activity during trials. In black we plot the accuracy of decoding the choice of arm direction (Dir), light position (Lig), and outcome (Out) for the current trial (left) and the previous trial (right). In gray we plot the decoding accuracy of shuffled labels across trials. Differences assessed using a paired Wilcoxon signed-rank test, *p < 0.05, **p < 0.01, #p < 0.005. Symbols plot means ± SEM across 49 sessions. b, As for a, plotted as relative decoding accuracy, the difference between the decoding accuracy of the data and of the mean of the shuffled data in that session. All p values are given for a Wilcoxon signed-rank test against zero median. c, No correlation between session performance and the accuracy of decoding the upcoming outcome of a trial. Gray line is a linear regression fit to the data. d, e, Per-subject breakdown of the decoding results in b. f, As for a, but for population activity during the intertrial intervals of each session, decoding the features of the following (left) and preceding (right) trials. g, as for b, for decoding during the intertrial interval. h, As for c, for decoding the outcome preceding the intertrial interval. i, j, Per-subject breakdown of the decoding results in g.

In contrast, from population activity during the intertrial interval we could decode the direction choice, outcome, and light position of the immediately preceding trial well above chance (Fig. 3f,g, right), which could even be seen for each rat despite the relatively small number (11–14) of sessions each performed (Fig. 3i). One caveat is that although the light was extinguished during the intertrial interval, precisely when is not clear from the data we have available (see above, Materials and Methods); consequently, it is possible that the decoding of the light could represent in part its ongoing state. Decoding of the past outcome also did not depend on performance in the session (Fig. 3h). Decoding the same features of the immediately following trial was at chance (Fig. 3 f,g, left, j). Thus, trial and intertrial activity both represented distinct states of the world. Moreover, the evidence suggests trial activity represented the present, and intertrial activity predominantly represented the past.

Independent decoding axes between the trials and intertrial intervals

Having found evidence that the activity of a single mPfC population occupies different subspaces encoding distinct states of the world, we could now ask whether and how the representations are kept distinct to downstream targets.

To compare the population coding between the trial and intertrial interval, we determined the decoding axis of trial activity for each of the present features and the decoding axis of intertrial interval activity for those same features in the preceding trial (see above, Materials and Methods). These decoding axes were close to orthogonal for all three features: the angles cluster at or close to π/2 (or, equivalently, their dot-product clusters at or around zero; Fig. 4a). And although the decoding axes for direction choice and light position departed from purely orthogonal, the median departure was small, being 0.067π for direction and 0.045π for light position. These differences between trial and intertrial decoding axes were also consistently and substantially larger than the differences within the same phase (Fig. 4b). Thus, the state of the world in the trial and intertrial interval can be independently decoded from the same mPfC population.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Independent decoding of task features between trials and intertrial intervals. a, Angles between the decoding weight vectors for trials and intertrial intervals for each task feature. For reference, we also compute the angle between trial and intertrial interval decoding vectors obtained by training on shuffled label data (gray). Box plots show median (line), interquartile range (box), and 95% interval (tails) over all sessions. The p values are from Wilcoxon rank sum tests for the difference from π/2. b, Variation of angles between decoding axes within the same task phase. We plot here the average standard deviation (STD) of the angles between decoding axes for the same feature in either the trials (left) or intertrial intervals (right) of a session. Within-session decoding axes for each feature were taken from the classifiers fitted to each training subset of trials or intertrial intervals during the cross-validation. c, As for a, but comparing the decoding weight vectors between features within trials. d, As for c, within intertrial intervals. e, Cross-decoding performance for each task feature of the current trial, using the decoding vectors in a. Left, Performance when the decoder was trained on activity during trials and tested on activity in the intertrial intervals. Black dashed line shows the chance level obtained training the classifier on shuffled labels for the trials and testing on intertrial intervals given the same shuffled labels. Right, Performance when the decoder was trained on activity from the intertrial intervals, and tested on activity in the trials; *p < 0.05, **p < 0.01, #p < 0.005. f, As for e, using leave-one-out cross-decoding.

By contrast, within the trial and the intertrial interval, pairs of decoding axes for different features were not close to orthogonal, except the direction and outcome axes (Fig. 4c,d). This neatly demonstrates that the near orthogonality of the decoding axes between trials and intertrial intervals is not then a trivial consequence of the decoding axes being random vectors drawn from the same distribution because the decoding axes of the same dimension within each phase are not orthogonal. Notably, the distributions of angles between the decoding axes for a given pair of features were preserved between the trials and the intertrial intervals, with outcome direction around π/2, light direction centered below π/2, and light outcome centered above π/2. Thus, although each decoding axis rotated close to orthogonal between the trial and intertrial interval, the relationships between the feature decoding axes were preserved.

To quantify how distinct these independent axes made the decoding of the trial and intertrial states, we cross-decoded one from the other: for each feature type, we trained the classifier on all trials of a session and tested its ability to decode the same feature from the following intertrial intervals. We found that cross-decoding was at chance level for both outcome and light position, and significant but weak for direction (Fig. 4e), consistent with the angles between their decoding axes in the trials and intertrial intervals (Fig. 4a). This result was robust whether we trained on trials and tested on intertrial intervals or vice versa. Cross-decoding was also weak or at chance if we used leave-one-out testing instead (Fig. 4f) by leaving out the ith trial and its following intertrial interval, training on N – 1 trials, and predicting the ith intertrial interval. Thus the near independent decoding axes (Fig. 4a) indeed imply that downstream targets could independently read out either the trial or the intertrial state of the task from mPfC population activity.

Decoding and cross-decoding are robust across types of session

We explored the extent to which this decoding depended on what occurred during each session. We first split the sessions by whether the target rule was direction based (15 sessions) or cue based (34 sessions). For trials, the present direction choice and outcome could still be significantly decoded for both types of rule, despite the considerable drop in power from the reduced number of sessions (Fig. 5a). For intertrial intervals, the preceding direction choice, outcome, and light position could still be decoded well above chance for both types of rule (Fig. 5b).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Breakdown of decoding and cross-decoding by session type. a, Breakdown of the trial decoding results in Figure 3b by the rule type of each session (15 direction rule sessions, 34 cue rule sessions); *p < 0.05, **p < 0.01, #p < 0.005. b, As for a, breakdown of the intertrial interval decoding results by the rule type of each session. c, Breakdown of the trial decoding results in Figure 3b by whether a rule was learned in a session (10 identified learning sessions, 39 other sessions). d, As for c, breakdown of the intertrial interval decoding results by learning and other sessions. e, f, Breakdown of the cross-decoding results in Figure 4d by rule type (left) and performance type (right) for each cross-decoding direction.

To determine whether learning itself affected any dependence on the task state, we then separated the sessions into two behavioral groups, putative learning sessions (n = 10), identified by a step change in task performance (see above, Materials and Methods), and the remaining sessions, called here Other (n = 39). We found decoding of task features was similar when comparing learning sessions and all Other sessions for both trials (Fig. 5c) and intertrial intervals (Fig. 5d). The sole exception, decoding the current light position during trials of Other sessions but not learning sessions, could be due either to a real effect or to the low power for decoding from 10 learning sessions.

For completeness, we also examined the breakdown of the cross-decoding results in Figure 4d by types of session. Figure 5, e and f, shows that cross-decoding of most features between trial and intertrial activity remained at chance, with again significant but weak cross-decoding of direction.

Evolution of decoding within trials and intertrial intervals

It is likely that the decoding of task features from mPfC activity is partly dependent on maze position (Ito et al., 2015; Spellman et al., 2015). To further examine the evolution of decoding over the trial and intertrial interval, we again divided the maze into five equally sized sections (Fig. 2g) and constructed population firing rate vectors for each position. Although the trials averaged only ∼6.5 s in duration, and so each position was occupied for ∼1 s, we still obtained clear evidence for decoding the direction choice, outcome, and light position of the current trial across multiple contiguous locations (Fig. 6a, left). The contrast between the strong decoding of the features of the current trial and the weak decoding of the features of the previous trial was even clearer across maze positions (Fig. 6b, right).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Population vector decoding across the maze. a, Decoder performance at each maze location. We plot here decoding at consecutive locations across the trial and following ITI, with the decoding target for each location being the task feature in the current trial (so decoding the preceding trial during the intertrial interval). Blue bars give significance values of p < 0.05, thin blue line; p < 0.01, medium thickness blue line; p < 0.005, thick blue line; for Wilcoxon signed-rank test for difference from zero. Symbols plot means ± SEM across n = 49 sessions. b, As per a, for decoding the features of the following trial from the current intertrial interval (left) and decoding the features of the previous trial from the current trial (right). c, Cross-decoding performance at each maze location is at chance.

This evolution means that there is contiguous decoding from the trial to the intertrial interval for all three features (Fig. 6a). Despite this contiguity, the cross-decoding between the same position in the two phases was at chance (Fig. 6c). In particular, cross-decoding at the arm end (position 5) was at chance, despite the rat continuously occupying this position during the transition from the trial to the intertrial interval. This suggests that the distinct decoding of the trial and intertrial states of the same feature appeared immediately at the arm end, or close to it (Fig. 2).

Figure 7 shows that these position-dependent decoding and cross-decoding results for trials are broadly robust to breaking them down by the type of rule or by learning behavior. Breakdowns of position decoding by session type in the intertrial intervals are given in Maggi et al. (2018), their Figure 5. In particular, we note here that the decoding of the state of the light during the intertrial interval only significantly occurs at position 5 when taken over all sessions (Fig. 6a, right), and as these data do not specify precisely when the light was extinguished during the interval, it is unknown whether that reflects the ongoing state of the light or the past state.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Robustness of location-dependent decoding in trials across session types. a, b, Breakdown of decoding performance in Figure 6a, according to the type of rule in each session, direction (15 sessions) or cued (34 sessions). c, d, Breakdown of decoding performance in Figure 6a, according to if the session contained putative rule learning (10 sessions, red symbols and lines) or not (39 sessions, blue symbols and lines). e, f, Breakdown of cross-decoding performance at each maze location (Figure 6c) by rule type (e) and performance type (f). Similar results were obtained for cross-decoding intertrial intervals from trials (data not shown).

Population representations of features reactivate in sleep

That the population activity occupies linearly separable subspaces between the trial and intertrial intervals (or the inbound and outbound phases) strongly suggests that the mPfC populations can be driven to either one or the other by upstream inputs. In turn, this implies that the representations of these two world states were independently addressable. To explore this question further, we turned to activity of the same populations during sleep.

Prior reports showed that patterns of mPfC population activity during training are preferentially repeated in post-training slow-wave sleep (Euston et al., 2007; Peyrache et al., 2009; Singh et al., 2019), consistent with a role in memory consolidation. However, these analyses looked only at specific templates or the reappearance of correlations between neurons, so it is unknown what task states these repeated patterns represented. Thus, we took advantage of the fact that our mPfC populations were also recorded during both pre-training and post-training sleep to ask whether their activity during sleep was specifically driven to either or both of the activity subspaces occupied by the population during the trials and intertrial intervals.

We first tested whether population activity representing features in the trials reactivated during slow-wave sleep. For each feature of the task happening in the present (e.g., choosing the left arm), we created the mean vector of population activity specific to that feature during trials in a session. This average population vector thus represented the region of the activity subspace (Fig. 2) occupied during trials with that feature. To seek reactivation of this region of the subspace in slow-wave sleep, we computed population firing rate vectors in pre-training and post-training slow-wave sleep in time bins of 1 s duration and correlated each sleep vector with the feature-specific trial vector (Fig. 8a). We thus obtained a distribution of correlations between the trial vector and all pre-training sleep vectors and a similar distribution between the trial vector and all post-training sleep vectors. Greater correlation with post-training sleep activity would then be evidence of preferential reactivation of feature-specific activity in post-training sleep.

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8.

Reactivation of population subspaces in post-training sleep. a, Example population activity vectors. Top, From one learning session, we plot the average firing rate vector for correct trials (Trial 1). For comparison, we also plot examples of firing rate vectors from pre-training and post-training slow-wave sleep (1 s bins). Neurons are ranked in order of their firing rates in the trial vector. Bottom, as for the top for an example session not classified as learning. b, Example distributions of Spearman's rank correlations between trial and sleep population activity. Top, For the same learning session as a, we plot the distributions of correlations between each vector of feature-specific trial activity and the population activity vectors in pre-training and post-training slow-wave sleep. Bottom, As for the top for the example nonlearning session in a. R, Right arm; 1 indicates rewarded trial. c, Summary of reactivations of trial activity across all sessions. For each task feature, we plot the difference between the medians of the pre-training and post-training correlation distributions. A difference greater than zero indicates greater correlation between activity in trials and post-training sleep. Each symbol is a session. Empty symbols are sessions with significantly different correlation distributions at p < 0.05 (Kolmogorov–Smirnov test). Gray filled symbols are not significantly different. One black circle for learning and one for nonlearning sessions identify the two example sessions in a and b. d, As for c, but plotting the median differences between distributions for paired task features within the same sleep epoch. For example, in the far-left column, we plot the difference between the correlations with pre-training sleep activity for right-choice (PreR) and left-choice (PreL) specific trial vectors. e, f, As for c and d, for intertrial interval activity. g, Relationship between population size and the strength of preferential reactivation of trial activity. Symbols show sessions; red are learning, blue are Other sessions. Solid line gives the best-fit linear regression; text gives Spearman's rank correlation. h, As for g, for intertrial interval activity.

We examined reactivation separately between learning and Other sessions, seeking consistency with previous reports that reactivation of waking population activity in mPfC most clearly occurs immediately after rule acquisition (Peyrache et al., 2009; Singh et al., 2019). Figure 8b (top) shows an example of a learning session with preferential reactivation. For all trial features, the distribution of correlations between the trial and post-training sleep population activity is right shifted from the distribution for pre-training sleep. For example, the population activity vector for choosing the right arm is more correlated with activity vectors in post-training (PostR) than pre-training (PreR) sleep.

Such post-training reactivation was not inevitable. In Figure 8b (bottom), we plot another example in which the trial-activity vector equally correlates with population activity in pre-training and post-training sleep. Although specific pairs of features (such as the left and right light positions) differed in their overall correlation between sleep and trial activity, no feature shows preferential reactivation in post-training sleep.

These examples were recapitulated across the data (Fig. 8c). In learning sessions, feature-specific activity vectors were consistently more correlated with activity in post-training than pre-training sleep. By contrast, the Other sessions showed less consistent preferential reactivation of any feature-specific activity vector in post-training sleep. As a control for statistical artifacts in our reactivation analysis, we looked for differences in reactivation between paired features (e.g., left vs right arm choice) within the same sleep epoch and found these all center on zero (Fig. 8d). Thus, population representations of present task features in the trials were preferentially reactivated in post-training sleep, and this most consistently occurred after a learning session.

We repeated the same analyses using feature-specific population vectors from the intertrial interval activity and also found evidence of preferential reactivation in some sessions (Fig. 8e,f). However, in contrast to trial activity, there was no consistent preferential reactivation of intertrial interval activity after a learning session.

Neither the preferential reactivation of trial nor intertrial activity was explained by significantly higher correlations between waking and sleep activity vectors from smaller populations (Fig. 8g,h).

As our measure of reactivation is asking whether and when the activity of the mPfC population revisits the trial and/or intertrial activity subspaces, it could do so on a range of time scales. These patterns of preferential reactivation were consistent across a range of bin sizes used to construct the activity vectors during sleep (Fig. 9). Notably, across these time scales, trial activity showed two independent properties from intertrial interval activity—consistent preferential reactivation after learning sessions, and preferential reactivation in those sessions was stronger at smaller bin sizes. These results are consistent with trial and intertrial activity subspaces being independently addressable; we thus sought further evidence of their independence.

Figure 9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 9.

Robust preferential reactivation across time scales of sleep activity. a, Preferential reactivation of trial activity across time scales. At each time bin used to construct population activity vectors in sleep, we plot the distribution over sessions of the median differences between pre-training and post-training correlation distributions for learning (top), and other (bottom) sessions. Distributions are plotted as the mean (thick lines) ± 2 SEM (thin lines); at the 1 s bin, these summarize the distributions shown in full in Figure 8c. Left, middle, and right columns plot two distributions, one per pair of features; lighter colors indicate left (L) or error trials (0), whereas darker colors indicate right (R) or correct trials (1). Time bins range from 100 ms to 10 s, tested every 150 ms. Dotted lines at the top of each panel indicate bins with reactivation significantly above zero (Wilcoxon signed-rank test, p < 0.05 small dot; p < 0.005 midsize dot; p < 0.001 large dot; N = 10 learning, N = 39 Other sessions). b, As for a, for intertrial activity.

Independent properties of trial and intertrial activity reactivation in sleep

We asked whether the amount of reactivation of population activity differed between trial and intertrial activity. The reactivation of trial population activity was strongly correlated between pre-training and post-training sleep (Fig. 10a), but the reactivation of intertrial interval activity was less correlated (Fig. 10b), and this was consistent across time scales used to construct the sleep activity vectors (Fig. 10c). Thus the overall reactivation of trial and intertrial interval activity was consistently different, again suggestive that the two subspaces of activity were independently addressable.

Figure 10.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 10.

Trial activity reactivation is more correlated between pre-training and post-training sleep. a, Correlation between the reactivation of trial activity in pre-training and post-training sleep for each task feature. Reactivation is as defined in Figure 8a, the median rank correlation between the waking and sleep activity vectors. Note the preferential reactivation is shown by symbols falling above the dashed identity line. Sleep activity vectors constructed using 1 s bins. CC, Correlation coefficient. Red indicates learning, blue indicates Other. b, As for a, for intertrial activity. c, Difference between trial and intertrial activity reactivation correlations across bin sizes used to construct sleep activity vectors.

Given the above evidence that reactivation of trial and intertrial interval activity could be independently controlled, we further asked whether they differed in how preferential reactivation correlated with behavior. Following the differences in reactivation after learning sessions (Fig. 9), we looked at the degree of learning in a session, which we quantified by the size of the change in reward rate in that session (see above, Materials and Methods). We found preferential reactivation of trial activity correlated with the change in reward rate (Fig. 11a), but preferential reactivation of intertrial activity did not (Fig. 11b). Again, this difference between trial and intertrial activity reactivation was consistent across a wide range of time scales used to construct the sleep activity vectors (Fig. 11c,d).

Figure 11.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 11.

Preferential reactivation of trial activity correlates with performance. a, Correlation between performance and preferential reactivation of trial activity. Sleep activity vectors constructed using 1 s bins. Values are Spearman's rank correlations. We plot here reactivation of vectors corresponding to left (direction and light) or correct; correlations for other feature-specific vectors are similar in magnitude, 0.37 (choose right), 0.35 (cue on right), 0.2 (error trials). b, As for a, for intertrial interval activity. Correlations for other feature-specific vectors are similar in magnitude, −0.004 (choose right), 0.02 (cue on right), −0.08 (error trials). c, Correlations between performance and preferential reactivation of trial activity across bin sizes used to construct sleep activity vectors for the features in a. d, As for c, for intertrial interval activity.

Discussion

Activity in the prefrontal cortex is known to represent different states of the world, including the immediate past or present state in a range of tasks (Baeg et al., 2003; Averbeck et al., 2006; Fujisawa et al., 2008; Sul et al., 2010; Rigotti et al., 2013; Hanks et al., 2015; Ito et al., 2015; Siegel et al., 2015; Spellman et al., 2015; Guise and Shapiro, 2017). How the representations of the different states relate to each other, and whether they coexist in the same population of neurons, has been unclear. Consequently, it is unknown how downstream readouts of prefrontal cortex activity can distinguish activity representing different states of the world.

Here, we have shown one potential solution in the medial prefrontal cortex of rats learning rules in a Y-maze: different states are encoded in the same population in such a way that linear decoders can read out different states of multiple features of the task. That encoding had two notable features. First, population activity is linearly separable between the trial and intertrial interval in as little as one dimension and so exists in different subspaces during these two phases of the task. Second, the decoding was roughly orthogonal between the trial and intertrial activity. These two features allow a simple solution to the interference problem.

The interference problem

Any neural population whose activity contains information about multiple states of the world faces the problem of interference (Libby and Buschman, 2021), of how downstream populations can distinguish the activity that depends on each state, so that the sequence and causality of world events is clear. The inverse problem is how inputs to the population can selectively recall only the activity that depends on a particular state.

As we have shown here, because trial and intertrial activity occupies different subspaces of the population activity, a downstream target using a linear decoder can distinguish the two (Semedo et al., 2019). This suggests a simple solution to the interference problem of having two downstream populations, one whose input weights from the mPfC population match the decoding axis for the trial state and another whose input weights from the mPfC population match the decoding axis for the intertrial state. Then the first downstream population only responds to activity representing the state of the trial and the other only to activity representing the state of the intertrial.

Key here is that the decoding axes are orthogonal, or close to it, although the population activity in mPfC is not. In Figure 12 we show this by plotting the angles between the mean activity vectors representing each feature in trials and intertrial intervals: we see that the activity representing each feature is more closely aligned between trials and intertrial intervals than are the corresponding decoding axes. Despite this alignment, because the activity sits in different linearly separable subspaces between trials and intertrial intervals, the different states of the task in the trial and in the intertrial interval are easily distinguishable by a linear decoder.

Figure 12.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 12.

Population activity is not orthogonal, only decoding axes. For each feature (e.g., choose left), we plot the distribution of angles between the corresponding mean activity vector in trials and intertrial intervals. We also replot from Figure 4a the distribution of angles between the trial and intertrial interval decoding axes for that feature. L, Left; R, right. The number 1 indicates reward, 0 indicates no reward.

We also found evidence here of a solution to the inverse problem, as the existence of different subspaces of activity between trials and intertrial intervals means that upstream inputs could in principle separately drive the population activity to either subspace. The sequential decoding also strongly suggests that the same mPfC population can be driven into different representations by upstream inputs.

To explore this further we looked at activity of the same mPfC populations during sleep to ask whether trial and intertrial intervals representations of the task features are reactivated differently. Trial and intertrial activity were both preferentially reactivated in post-training slow-wave sleep, yet we found evidence that preferential reactivation of trial activity differed in the following four ways: the time scales at which it occurred most strongly, it occurred after learning sessions, the strength of reactivation was more consistent between pre-training and post-training sleep, and it correlated with the performance of the rats in the sessions. Together, these differences between the reactivation of trial and intertrial interval activity are consistent with upstream inputs to the mPfC population being able to separately address the representations of these states.

The consistency of preferential reactivation across broad time scales suggests that it is the changes to the relative excitability of neurons within the mPfC population that are carried forward into sleep (Singh et al., 2019). Thus, this consistency across broad time scales implies that whenever the neurons of the population are active, they are active together with approximately the same ordering of firing rates.

Mixed population coding in mPfC

Our finding that small mPfC populations can sustain mixed encoding of two or more of direction choice, light position, and outcome of the current trial is consistent with prior reports of mixed or multiplexed coding by single neurons in the prefrontal cortex (Jung et al., 1998; Horst and Laubach, 2012; Rigotti et al., 2013; Fusi et al., 2016; Aoi et al., 2020). These encodings were also position dependent. Decoding of direction choice reliably occurred from the choice point of the maze onward, but it is unclear whether this represents a causal role in the choice itself or an ongoing representation of a choice being made.

Indeed, we are not claiming that the specific task features we decoded are necessarily explicitly represented in mPfC population activity. Rather, throughout we have interpreted the decoding of these features as evidence that mPfC population activity is at least representing the state of the world, similar to reinforcement learning views of PfC representations (Wang et al., 2018), because these features are a part of that state; and, hence, any change in one of those features, such as arm choice, would thus be a different state of the world.

Previous studies have reported that past choices modulated mPfC population activity during trials (Baeg et al., 2003; Sul et al., 2010). In contrast to the robust decoding of the present, we found weak evidence that mPfC activity during a trial depended on the light position of the previous trial and weak evidence that it depended on the direction choice of the previous trial only during direction-based rules. Moreover, these features of the past could only be decoded at one or two locations on the maze. Thus, during trials, population activity in the prefrontal cortex had robust, sustained dependence on multiple features of the present but at best weakly and transiently depended on one feature of the past.

Indeed, we have evidence here that the trial and intertrial activity represent not just different task states but, respectively, the present and past state of the task. Trial activity decoded present but not past features; intertrial activity decoded features of the preceding trial. The latter is consistent with well-established roles for the prefrontal cortex in short-term memory (Funahashi et al., 1989; Machens et al., 2010; Constantinidis et al., 2018; Lundqvist et al., 2018). However, the limitations of the Y-maze task data mean we cannot rule out that the intertrial activity also represented some features of the present during that interval, which is a question to be pursued further. Nonetheless, we have strong evidence that mPfC activity represents distinct states in the trials and in the intertrial intervals.

What could drive changes in mPfC population activity

The evolution of activity within trials and intertrial intervals was continuous, with adjacent maze sections containing more similar population activity, yet the transition from the trial to the intertrial interval was discontinuous, with population activity moving to a different subspace, linearly separable from the trials' subspace. What might be driving this shift from the trial to the intertrial interval subspace of activity and hence its decodability?

The division into trials and intertrial intervals or the inbound and outbound phases in Figure 2 both distinguish the two legs of the journey in the maze. During the return trip to the starting position, the change in context and direction of movement would likely change the signals available to the mPfC. It does not automatically follow though that changes in context and movement cause the observed changes in population activity in mPfC; those changes to sensory and movement information could have changed mPfC population activity so that it did not encode anything about the immediately preceding trial, in the same way, for example, that we showed the intertrial activity encodes nothing about the immediately upcoming trial, even when the decision of the trial could be known in advance. Thus, our finding that we could still decode the state of the immediately preceding task features from intertrial activity despite the changes in context and movement information is nontrivial. Indeed, it implies that those changes could be the drivers of the observed changes in population activity.

This suggests multiple lines of fruitful further work here. One open question is which inputs to the mPfC drive the move from one activity subspace to another. Given the switch in context and heading direction, a likely candidate is the direct input from region CA1 of the hippocampus (Jones and Wilson, 2005; Benchenane et al., 2010, 2011). Another open question is precisely when the change in activity subspace happens. We showed here preliminary results that the larger, discontinuous change in population activity could be plausibly either on reaching the arm end or on initiating the outbound trip back to the starting position. Another is the precise function of the representations of the trial and intertrial interval; one possibility is they respectively reflect reward prediction and reward processing. One way to tackle this question would be to examine how much the clean independence between the decoding of task states depends on the behavioral task. For example, tasks where the future choice of arm depends on recent history, such as double-ended T-mazes (Jones and Wilson, 2005), multiarm sequence mazes (Poucet et al., 1991), or delayed nonmatch to place (Spellman et al., 2015), blur the separation of the present and the past. Comparing population-level decoding of the states in such tasks would give useful insights into when they are or are not independently coded within mPfC.

Footnotes

  • This work was supported by the Medical Research Council Grants MR/J008648/1, MR/P005659/1, and MR/S025944/1. The original data collection was supported by the European Union Framework (FP6) ICEA grant. We thank Adrien Peyrache for the data, discussions, and comments on early drafts of this manuscript; Hazem Toutounji and Martin O'Neill for comments on drafts; and the Humphries' lab, past and present (Abhinav Singh, Javier Caballero, Mat Evans, Francois Cinotti, Tomas Fiers), for discussion.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Mark D. Humphries at mark.humphries{at}nottingham.ac.uk

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

References

  1. ↵
    1. Aoi MC,
    2. Mante V,
    3. Pillow JW
    (2020) Prefrontal cortex exhibits multidimensional dynamic encoding during decision-making. Nat Neurosci 23:1410–1420. doi:10.1038/s41593-020-0696-5 pmid:33020653
    OpenUrlCrossRefPubMed
  2. ↵
    1. Averbeck BB,
    2. Sohn JW,
    3. Lee D
    (2006) Activity in prefrontal cortex during dynamic selection of action sequences. Nat Neurosci 9:276–282. doi:10.1038/nn1634 pmid:16429134
    OpenUrlCrossRefPubMed
  3. ↵
    1. Baeg EH,
    2. Kim YB,
    3. Huh K,
    4. Mook-Jung I,
    5. Kim HT,
    6. Jung MW
    (2003) Dynamics of population code for working memory in the prefrontal cortex. Neuron 40:177–188. doi:10.1016/s0896-6273(03)00597-x pmid:14527442
    OpenUrlCrossRefPubMed
  4. ↵
    1. Benchenane K,
    2. Peyrache A,
    3. Khamassi M,
    4. Tierney PL,
    5. Gioanni Y,
    6. Battaglia FP,
    7. Wiener SI
    (2010) Coherent theta oscillations and reorganization of spike timing in the hippocampal- prefrontal network upon learning. Neuron 66:921–936. doi:10.1016/j.neuron.2010.05.013 pmid:20620877
    OpenUrlCrossRefPubMed
  5. ↵
    1. Benchenane K,
    2. Tiesinga PH,
    3. Battaglia FP
    (2011) Oscillations in the prefrontal cortex: a gateway to memory and attention. Curr Opin Neurobiol 21:475–485. doi:10.1016/j.conb.2011.01.004 pmid:21429736
    OpenUrlCrossRefPubMed
  6. ↵
    1. Constantinidis C,
    2. Funahashi S,
    3. Lee D,
    4. Murray JD,
    5. Qi XL,
    6. Wang M,
    7. Arnsten AFT
    (2018) Persistent spiking activity underlies working memory. J Neurosci 38:7020–7028. doi:10.1523/JNEUROSCI.2486-17.2018 pmid:30089641
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Denovellis EL,
    2. Gillespie AK,
    3. Coulter ME,
    4. Sosa M,
    5. Chung JE,
    6. Eden UT,
    7. Frank LM
    (2021) Hippocampal replay of experience at real-world speeds. Elife 10:e64505. doi:10.7554/eLife.64505
    OpenUrlCrossRef
  8. ↵
    1. Durstewitz D,
    2. Vittoz NM,
    3. Floresco SB,
    4. Seamans JK
    (2010) Abrupt transitions between prefrontal neural ensemble states accompany behavioral transitions during rule learning. Neuron 66:438–448. doi:10.1016/j.neuron.2010.03.029 pmid:20471356
    OpenUrlCrossRefPubMed
  9. ↵
    1. Euston DR,
    2. Tatsuno M,
    3. McNaughton BL
    (2007) Fast-forward playback of recent memory sequences in prefrontal cortex during sleep. Science 318:1147–1150. doi:10.1126/science.1148979 pmid:18006749
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Euston DR,
    2. Gruber AJ,
    3. McNaughton BL
    (2012) The role of medial prefrontal cortex in memory and decision making. Neuron 76:1057–1070. doi:10.1016/j.neuron.2012.12.002 pmid:23259943
    OpenUrlCrossRefPubMed
  11. ↵
    1. Fujisawa S,
    2. Amarasingham A,
    3. Harrison MT,
    4. Buzsáki G
    (2008) Behavior-dependent short-term assembly dynamics in the medial prefrontal cortex. Nat Neurosci 11:823–833. doi:10.1038/nn.2134 pmid:18516033
    OpenUrlCrossRefPubMed
  12. ↵
    1. Funahashi S,
    2. Bruce CJ,
    3. Goldman-Rakic PS
    (1989) Mnemonic coding of visual space in the monkey's dorsolateral prefrontal cortex. J Neurophysiol 61:331–349. doi:10.1152/jn.1989.61.2.331 pmid:2918358
    OpenUrlCrossRefPubMed
  13. ↵
    1. Fusi S,
    2. Miller EK,
    3. Rigotti M
    (2016) Why neurons mix: high dimensionality for higher cognition. Curr Opin Neurobiol 37:66–74. doi:10.1016/j.conb.2016.01.010 pmid:26851755
    OpenUrlCrossRefPubMed
  14. ↵
    1. Guise KG,
    2. Shapiro ML
    (2017) Medial prefrontal cortex reduces memory interference by modifying hippocampal encoding. Neuron 94:183–192.e8. doi:10.1016/j.neuron.2017.03.011 pmid:28343868
    OpenUrlCrossRefPubMed
  15. ↵
    1. Hanks TD,
    2. Kopec CD,
    3. Brunton BW,
    4. Duan CA,
    5. Erlich JC,
    6. Brody CD
    (2015) Distinct relationships of parietal and prefrontal cortices to evidence accumulation. Nature 520:220–223. doi:10.1038/nature14066 pmid:25600270
    OpenUrlCrossRefPubMed
  16. ↵
    1. Hastie T,
    2. Tibshirani R,
    3. Friedman J
    (2009) The elements of statistical learning. Berlin: Springer.
  17. ↵
    1. Horst NK,
    2. Laubach M
    (2012) Working with memory: evidence for a role for the medial prefrontal cortex in performance monitoring during spatial delayed alternation. J Neurophysiol 108:3276–3288. doi:10.1152/jn.01192.2011
    OpenUrlCrossRefPubMed
  18. ↵
    1. Ito HT,
    2. Zhang SJ,
    3. Witter MP,
    4. Moser EI,
    5. Moser MB
    (2015) A prefrontal–thalamo–hippocampal circuit for goal-directed spatial navigation. Nature 522:50–55. doi:10.1038/nature14396 pmid:26017312
    OpenUrlCrossRefPubMed
  19. ↵
    1. Jones MW,
    2. Wilson MA
    (2005) Theta rhythms coordinate hippocampal-prefrontal interactions in a spatial memory task. PLoS Biol 3:e402. doi:10.1371/journal.pbio.0030402 pmid:16279838
    OpenUrlCrossRefPubMed
  20. ↵
    1. Jung MW,
    2. Qin Y,
    3. McNaughton BL,
    4. Barnes CA
    (1998) Firing characteristics of deep layer neurons in prefrontal cortex in rats performing spatial working memory tasks. Cereb Cortex 8:437–450. doi:10.1093/cercor/8.5.437 pmid:9722087
    OpenUrlCrossRefPubMed
  21. ↵
    1. Karlsson MP,
    2. Tervo DG,
    3. Karpova AY
    (2012) Network resets in medial prefrontal cortex mark the onset of behavioral uncertainty. Science 338:135–139. doi:10.1126/science.1226518 pmid:23042898
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Laskowski CS,
    2. Williams RJ,
    3. Martens KM,
    4. Gruber AJ,
    5. Fisher KG,
    6. Euston DR
    (2016) The role of the medial prefrontal cortex in updating reward value and avoiding perseveration. Behav Brain Res 306:52–63. doi:10.1016/j.bbr.2016.03.007 pmid:26965571
    OpenUrlCrossRefPubMed
  23. ↵
    1. Laubach M,
    2. Caetano MS,
    3. Narayanan NS
    (2015) Mistakes were made: neural mechanisms for the adaptive control of action initiation by the medial prefrontal cortex. J Physiol Paris 109:104–117. doi:10.1016/j.jphysparis.2014.12.001 pmid:25636373
    OpenUrlCrossRefPubMed
  24. ↵
    1. Laubach M,
    2. Amarante LM,
    3. Swanson K,
    4. White SR
    (2018) What, if anything, is rodent prefrontal cortex? eNeuro 5:ENEURO.0315-18.2018. doi:10.1523/ENEURO.0315-18.2018
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Lee AK,
    2. Wilson MA
    (2002) Memory of sequential experience in the hippocampus during slow wave sleep. Neuron 36:1183–1194. doi:10.1016/s0896-6273(02)01096-6 pmid:12495631
    OpenUrlCrossRefPubMed
  26. ↵
    1. Libby A,
    2. Buschman TJ
    (2021) Rotational dynamics reduce interference between sensory and memory representations. Nat Neurosci 24:715–726. doi:10.1038/s41593-021-00821-9
    OpenUrlCrossRef
  27. ↵
    1. Lundqvist M,
    2. Herman P,
    3. Miller EK
    (2018) Working memory: delay activity, yes! persistent activity? maybe not. J Neurosci 38:7013–7019. doi:10.1523/JNEUROSCI.2485-17.2018 pmid:30089640
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Machens CK,
    2. Romo R,
    3. Brody CD
    (2010) Functional, but not anatomical, separation of “what” and “when” in prefrontal cortex. J Neurosci 30:350–360. doi:10.1523/JNEUROSCI.3276-09.2010
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Maggi S,
    2. Peyrache A,
    3. Humphries MD
    (2018) An ensemble code in medial prefrontal cortex links prior events to outcomes during learning. Nat Commun 9:2204. doi:10.1038/s41467-018-04638-2
    OpenUrlCrossRefPubMed
  30. ↵
    1. Narayanan NS,
    2. Laubach M
    (2008) Neuronal correlates of post-error slowing in the rat dorsomedial prefrontal cortex. J Neurophysiol 100:520–525. doi:10.1152/jn.00035.2008 pmid:18480374
    OpenUrlCrossRefPubMed
  31. ↵
    1. O'Neill J,
    2. Senior TJ,
    3. Allen K,
    4. Huxter JR,
    5. Csicsvari J
    (2008) Reactivation of experience-dependent cell assembly patterns in the hippocampus. Nat Neurosci 11:209–215. doi:10.1038/nn2037
    OpenUrlCrossRefPubMed
  32. ↵
    1. Peyrache A,
    2. Khamassi M,
    3. Benchenane K,
    4. Wiener SI,
    5. Battaglia FP
    (2009) Replay of rule-learning related neural patterns in the prefrontal cortex during sleep. Nat Neurosci 12:919–926. doi:10.1038/nn.2337 pmid:19483687
    OpenUrlCrossRefPubMed
  33. ↵
    1. Peyrache A,
    2. Benchenane K,
    3. Khamassi M,
    4. Wiener S,
    5. Battaglia F
    (2010) Principal component analysis of ensemble recordings reveals cell assemblies at high temporal resolution. J Comput Neurosci 29:309–325. doi:10.1007/s10827-009-0154-6 pmid:19529888
    OpenUrlCrossRefPubMed
  34. ↵
    1. Poucet B,
    2. Lucchessi H,
    3. Thinus-Blanc C
    (1991) What information is used by rats to update choices in the radial-arm maze? Behav Processes 25:15–26. doi:10.1016/0376-6357(91)90042-X pmid:24923799
    OpenUrlCrossRefPubMed
  35. ↵
    1. Powell NJ,
    2. Redish AD
    (2016) Representational changes of latent strategies in rat medial prefrontal cortex precede changes in behaviour. Nat Commun 7:12830. doi:10.1038/ncomms12830 pmid:27653278
    OpenUrlCrossRefPubMed
  36. ↵
    1. Rich EL,
    2. Shapiro M
    (2009) Rat prefrontal cortical neurons selectively code strategy switches. J Neurosci 29:7208–7219. doi:10.1523/JNEUROSCI.6068-08.2009 pmid:19494143
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Rich EL,
    2. Shapiro ML
    (2007) Prelimbic/infralimbic inactivation impairs memory for multiple task switches, but not flexible selection of familiar tasks. J Neurosci 27:4747–4755. doi:10.1523/JNEUROSCI.0369-07.2007 pmid:17460087
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Rigotti M,
    2. Barak O,
    3. Warden MR,
    4. Wang XJ,
    5. Daw ND,
    6. Miller EK,
    7. Fusi S
    (2013) The importance of mixed selectivity in complex cognitive tasks. Nature 497:585–590. doi:10.1038/nature12160 pmid:23685452
    OpenUrlCrossRefPubMed
  39. ↵
    1. Russo E,
    2. Ma T,
    3. Spanagel R,
    4. Durstewitz D,
    5. Toutounji H,
    6. Köhr G
    (2021) Coordinated prefrontal state transition leads extinction of reward-seeking behaviors. J Neurosci 41:2406–2419. doi:10.1523/JNEUROSCI.2588-20.2021 pmid:33531416
    OpenUrlAbstract/FREE Full Text
  40. ↵
    1. Semedo JD,
    2. Zandvakili A,
    3. Machens CK,
    4. Yu BM,
    5. Kohn A
    (2019) Cortical areas interact through a communication subspace. Neuron 102:249–259.e4. doi:10.1016/j.neuron.2019.01.026 pmid:30770252
    OpenUrlCrossRefPubMed
  41. ↵
    1. Siegel M,
    2. Buschman TJ,
    3. Miller EK
    (2015) Cortical information flow during flexible sensorimotor decisions. Science 348:1352–1355. doi:10.1126/science.aab0551 pmid:26089513
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Singh A,
    2. Peyrache A,
    3. Humphries MD
    (2019) Medial prefrontal cortex population activity is plastic irrespective of learning. J Neurosci 39:3470–3483. doi:10.1523/JNEUROSCI.1370-17.2019 pmid:30814311
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Skaggs WE,
    2. McNaughton BL
    (1996) Replay of neuronal firing sequences in rat hippocampus during sleep following spatial experience. Science 271:1870–1873. doi:10.1126/science.271.5257.1870 pmid:8596957
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Spellman T,
    2. Rigotti M,
    3. Ahmari SE,
    4. Fusi S,
    5. Gogos JA,
    6. Gordon JA
    (2015) Hippocampal-prefrontal input supports spatial encoding in working memory. Nature 522:309–314. doi:10.1038/nature14445 pmid:26053122
    OpenUrlCrossRefPubMed
  45. ↵
    1. Sul JH,
    2. Kim H,
    3. Huh N,
    4. Lee D,
    5. Jung MW
    (2010) Distinct roles of rodent orbitofrontal and medial prefrontal cortex in decision making. Neuron 66:449–460. doi:10.1016/j.neuron.2010.03.033 pmid:20471357
    OpenUrlCrossRefPubMed
  46. ↵
    1. Wang JX,
    2. Kurth-Nelson Z,
    3. Kumaran D,
    4. Tirumala D,
    5. Soyer H,
    6. Leibo JZ,
    7. Hassabis D,
    8. Botvinick M
    (2018) Prefrontal cortex as a meta-reinforcement learning system. Nat Neurosci 21:860–868. doi:10.1038/s41593-018-0147-8 pmid:29760527
    OpenUrlCrossRefPubMed
  47. ↵
    1. Young JJ,
    2. Shapiro ML
    (2009) Double dissociation and hierarchical organization of strategy switches and reversals in the rat pfc. Behav Neurosci 123:1028–1035. doi:10.1037/a0016822 pmid:19824768
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 42 (20)
Journal of Neuroscience
Vol. 42, Issue 20
18 May 2022
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Activity Subspaces in Medial Prefrontal Cortex Distinguish States of the World
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Activity Subspaces in Medial Prefrontal Cortex Distinguish States of the World
Silvia Maggi, Mark D. Humphries
Journal of Neuroscience 18 May 2022, 42 (20) 4131-4146; DOI: 10.1523/JNEUROSCI.1412-21.2022

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Activity Subspaces in Medial Prefrontal Cortex Distinguish States of the World
Silvia Maggi, Mark D. Humphries
Journal of Neuroscience 18 May 2022, 42 (20) 4131-4146; DOI: 10.1523/JNEUROSCI.1412-21.2022
Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • decoding
  • learning
  • medial prefrontal cortex
  • population activity
  • sleep reactivation

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Regional Excitatory–Inhibitory Balance Relates to Self-Reference Effect on Recollection via the Precuneus/Posterior Cingulate Cortex–Medial Prefrontal Cortex Connectivity
  • Attention Alters Population Spatial Frequency Tuning
  • Sequence Diversity and Encoded Enzymatic Differences of Monocistronic L1 ORF2 mRNA Variants in the Aged Normal and Alzheimer's Disease Brain
Show more Research Articles

Systems/Circuits

  • Hippocampal Sharp-Wave Ripples Decrease during Physical Actions Including Consummatory Behavior in Immobile Rodents
  • Developmental Olfactory Dysfunction and Abnormal Odor Memory in Immune-Challenged Disc1+/− Mice
  • Functional Roles of Gastrin-Releasing Peptide-Producing Neurons in the Suprachiasmatic Nucleus: Insights into Photic Entrainment and Circadian Regulation
Show more Systems/Circuits
  • Home
  • Alerts
  • Follow SFN on BlueSky
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Notice
  • Contact
  • Accessibility
(JNeurosci logo)
(SfN logo)

Copyright © 2025 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.