Postsynaptic burst reactivation of hippocampal neurons enables associative plasticity of temporally discontiguous inputs

Abstract
Editor's evaluation
Introduction
Results
Discussion
Methods
Data availability
References
Article and author information
Metrics

Abstract

A fundamental unresolved problem in neuroscience is how the brain associates in memory events that are separated in time. Here, we propose that reactivation-induced synaptic plasticity can solve this problem. Previously, we reported that the reinforcement signal dopamine converts hippocampal spike timing-dependent depression into potentiation during continued synaptic activity (Brzosko et al., 2015). Here, we report that postsynaptic bursts in the presence of dopamine produce input-specific LTP in mouse hippocampal synapses 10 min after they were primed with coincident pre- and post-synaptic activity (post-before-pre pairing; Δt = –20 ms). This priming activity induces synaptic depression and sets an NMDA receptor-dependent silent eligibility trace which, through the cAMP-PKA cascade, is rapidly converted into protein synthesis-dependent synaptic potentiation, mediated by a signaling pathway distinct from that of conventional LTP. This synaptic learning rule was incorporated into a computational model, and we found that it adds specificity to reinforcement learning by controlling memory allocation and enabling both ‘instructive’ and ‘supervised’ reinforcement learning. We predicted that this mechanism would make reactivated neurons activate more strongly and carry more spatial information than non-reactivated cells, which was confirmed in freely moving mice performing a reward-based navigation task.

Editor's evaluation

This article contains fundamental findings that substantially advance understanding of an important research question, mostly using an appropriate and validated methodology in line with the current state-of-the-art, with good and convincing support for the claims. The message of the article will have a profound and lasting influence on neuroscience.

https://doi.org/10.7554/eLife.81071.sa0

Introduction

For an animal to successfully adjust its behavior to changing environmental demands, it needs to learn to associate a sequence of events or actions to a subsequent outcome, for example a reward, which may be experienced minutes or hours later. How this is achieved in the brain remains unknown. Even the longest timescale synaptic plasticity event reported so far, behavioral timescale synaptic plasticity (BTSP; Bittner et al., 2017; Milstein et al., 2021), cannot bridge this temporal gap. In machine learning, reinforcement learning algorithms solve this problem by storing a temporary record of the occurrence of an event, known as an eligibility trace, which can later undergo learning changes once the outcome is known (Sutton and Barto, 2018). In biology, dopamine (DA) is thought to represent such a reinforcement learning signal, converting a temporary eligibility trace into a lasting synaptic change (Ljungberg et al., 1992; Schultz et al., 1997; Frémaux and Gerstner, 2016; Gerstner et al., 2018). However, the use of a brief scalar signal to alter the relevant synaptic weights raises two fundamental problems, first, that of time scales, and second, that of specificity (or credit assignment). Whereas initial Pavlovian conditioning of reward-predicting stimuli has optimal stimulus-reward intervals of 1.5–3.0 s, other learning mechanisms occur in the second, minute, and even hour ranges (Schultz, 2007). Moreover, with a longer time delay, the relationship between the preceding behavior and the outcome is less clear, questioning what events or actions before the reward should be associated with the outcome. One possible mechanism that could link previous activity (e.g. exploration of an environment) with a specific outcome (e.g. finding a reward) is neuronal reactivation, or replay activity, which in the hippocampus is enhanced by reward (Singer and Frank, 2009; Ambrose et al., 2016).

The rodent hippocampus is important for spatial memories. Spatial representations are built during exploration of an environment, when the hippocampus shows theta activity, and is later reactivated or replayed both during sleep (Wilson and McNaughton, 1994; Nádasdy et al., 1999; Lee and Wilson, 2002) and in the awake state (Kudrimoti et al., 1999; Foster and Wilson, 2006; Csicsvari et al., 2007; Diba and Buzsáki, 2007). Reactivation occurs during sharp wave ripples (SWRs), when neurons typically fire action potentials in brief bursts (Buzsáki et al., 1992; Diba and Buzsáki, 2007). It has been reported that stimulation of dopaminergic input promotes reactivation of hippocampal cell assemblies and memory persistence (McNamara et al., 2014).

We investigated the effect of action potential bursts on synaptic plasticity in individual postsynaptic hippocampal CA1 neurons (SWR-associated ‘reactivation’) during dopaminergic modulation (‘reward signal’) after they had undergone a spike pairing protocol (prior exploration-based synaptic ‘priming’). We found that the pairing protocol set an NMDA receptor-dependent silent eligibility trace, which could be converted several minutes later by burst activity in the presence of DA into protein synthesis-dependent long-term potentiation (LTP) mediated by a signaling pathway distinct from that of conventional LTP. Using this synaptic learning rule in a computational model we show that reactivation-induced plasticity increases specificity to reinforcement learning, offering a candidate mechanism of credit assignment in neural networks. To investigate how reactivation affects the functional properties of neurons in vivo, we used chronic calcium imaging of the dorsal CA1 region of the hippocampus in freely moving mice performing a reward-location learning task. We defined neuronal reactivation as activity during immobility at the reward-location in neurons that were already previously active during the reward approach. Neurons that reactivated after finding the reward had increased calcium activity and place map peaks compared to non-reactivated neurons, suggestive of changes in synaptic weights in reactivated neurons.

Results

Reactivation during dopaminergic modulation induces LTP

To investigate the effect of burst reactivation of individual postsynaptic CA1 cells during dopaminergic modulation, we monitored the synaptic weights of afferent synapses that had previously undergone a spike timing-dependent synaptic priming protocol. For this, we used whole-cell recording of CA1 pyramidal neurons in mouse hippocampal slices (Figure 1A). To be able to distinguish between conventional LTP and reactivation-induced LTP, we used a priming protocol that induces synaptic depression (Andrade-Talavera et al., 2016). Single postsynaptic action potentials followed by presynaptic input stimulation (post-before-pre protocol, Δt = –20 ms) led to input-specific synaptic depression in the test pathway (t-LTD; 61% ± 11% vs 100%, t(9) = 3.7, p=0.005, n=10; Figure 1B and F). Application of DA alone, without resuming synaptic stimulation after this pairing protocol, did not affect the depression (53% ± 7% vs 100%, t(5) = 7.0, p=0.0009, n=6; Figure 1C and F blue trace). Strikingly, postsynaptic action potential bursts (5–6 APs) in the presence of DA, 10 min after the pairing protocol, triggered an immediate induction of synaptic potentiation (135% ± 14.9% vs 100%, t(8) = 2.4, p=0.044, n=9; Figure 1D and F red trace). Burst stimulation alone, in the absence of DA, did not prevent synaptic depression (72% ± 12% vs 100%, t(6) = 2.5, p=0.048, n=7; Figure 1E and F black trace). This result suggests that the pairing protocol sets an eligibility trace allowing activated synapses to be selectively altered minutes later by reactivation of the postsynaptic neuron in the presence of DA. To our knowledge, these are the longest-lasting synaptic eligibility traces reported in the brain.

Figure 1

Download asset Open asset

Postsynaptic burst reactivation induces LTP in the presence of dopamine (DA).

(A), Schematic of the experimental paradigm (top) and setup (bottom).↓=6 action potential bursts; Whole-cell recording in CA1 stratum pyramidale; electrical stimulation electrodes in stratum radiatum. Plasticity was induced in one pathway (Paired pathway), and a second pathway was used for stability control and for confirmation of input specificity (Control pathway). Normalized EPSP slopes were averaged and plotted as a function of time. (B), Post-before pre-pairing protocol leads to input-specific synaptic depression. Pairing protocol (Δt = –20 ms) induces t-LTD (black trace) and does not affect synaptic weights in control pathway (gray trace). (C), DA application after a post-before-pre pairing protocol (Δt = –20 ms) does not prevent t-LTD (+DA –Burst, blue trace) and does not affect synaptic weights in control pathway (gray trace). (D), Application of DA together with action potential bursts of the postsynaptic cell (indicated by black arrow) induces synaptic potentiation after a post-before-pre pairing protocol (Δt = –20 ms) (+DA +Burst, red trace) and does not affect synaptic weights in control pathway (gray trace). (E), The same protocol, without application of DA, leads to synaptic depression (–DA +Burst, black trace) and does not affect synaptic weights in control pathway (gray trace). (F), Summary of results. All traces show an EPSP before (1) and 40 min after pairing (2). Plots show averages of normalized EPSP slopes ± SEM.

Figure 1—source data 1 Normalized EPSP slopes of all recorded cells.: https://cdn.elifesciences.org/articles/81071/elife-81071-fig1-data1-v2.xlsx
Download elife-81071-fig1-data1-v2.xlsx

We next investigated the requirements for setting the eligibility trace. First, we found that synaptic potentiation was indeed not observed without prior spike pairing (93% ± 8% vs 100%, t(6) = 0.84, p=0.43, n=7; Figure 2A). Induction of hippocampal t-LTD requires metabotropic glutamate receptors (mGluRs; Andrade-Talavera et al., 2016). To investigate whether LTD or mGluR signaling is required for burst-induced potentiation, we applied the mGluR antagonist MPEP throughout the recording. This blocked t-LTD (94% ± 6% vs 100%, t(6) = 0.95, p=0.38, n=7; Figure 2B and C, blue; vs control t-LTD 69% ± 7%, n=7; t(11.5) = 2.67, p=0.02; Figure 2B and C, black) but burst-induced potentiation was intact (131% ± 10% vs 100%, t(6) = 3.07, p=0.021, n=7; Figure 2B and C, red), suggesting that setting the eligibility trace by spike pairing is distinct from the signaling mechanism that mediates t-LTD.

Figure 2

Download asset Open asset

Setting of the eligibility trace is independent of LTD but requires postsynaptic NMDARs.

(A), Application of DA and burst stimulation, but without pairing protocol, does not induce potentiation. (B), Post-before-pre pairing protocol (Δt = –20 ms) leads to synaptic depression (t-LTD) (black trace), which is blocked by MPEP (blue trace). MPEP does not block DA- and burst-induced potentiation (red trace). (C), Summary of results. (D), Application of AP5 during pairing blocks DA- and burst-induced potentiation. (E), Postsynaptic intracellular MK801 (iMK801) does not block t-LTD (blue trace) but blocks DA- and burst-induced potentiation (red trace). (F), Summary of results. All traces show an EPSP before (1) and 40 min after pairing (2). Plots show averages of normalized EPSP slopes ± SEM.

Figure 2—source data 1 Normalized EPSP slopes of all recorded cells.: https://cdn.elifesciences.org/articles/81071/elife-81071-fig2-data1-v2.xlsx
Download elife-81071-fig2-data1-v2.xlsx

Many forms of hippocampal plasticity require NMDA receptors (NMDARs; Shipton and Paulsen, 2014). We, therefore, asked if activation of NMDARs during the pairing protocol is required for the synaptic eligibility trace to be set. We found that application of the NMDAR antagonist d-AP5 during pairing (Δt = –20 ms) abolished both t-LTD and the subsequent burst-induced potentiation, resulting in no change of synaptic weights (100% ± 9% vs 100%, t(6) = 0.043, p=0.97, n=7; Figure 2D and F). We then investigated whether specifically postsynaptic NMDARs are required for the eligibility trace. Loading the postsynaptic cell with the NMDAR channel blocker MK801 through the recording pipette did not affect t-LTD (68% ± 7% vs 100%, t(5) = 4.77, p=0.005, n=6; Figure 2E and F) but completely abolished burst-induced potentiation, leaving synaptic depression instead (61% ± 10% vs 100%, t(6) = 3.85, p=0.0084, n=6; Figure 2E and F). These results show a double dissociation between t-LTD induced by mGluR and non-postsynaptic NMDAR signaling and the eligibility trace for reactivation-induced potentiation set by postsynaptic ionotropic NMDARs.

We next investigated the mechanism that induces synaptic potentiation during burst stimulation. It was previously shown that activation of NMDARs after pairing is necessary to induce DA-dependent potentiation with subthreshold synaptic stimulation (Brzosko et al., 2015). However, the NMDAR antagonist d-AP5, applied after the pairing protocol but before burst stimulation, did not affect burst-induced potentiation (127% ± 9%, t(5) = 2.94, p=0.032, n=6; Figure 3A and C). In contrast, voltage-gated calcium channels (VGCCs) are required as application of nimodipine, an L-type VGCC blocker, before the burst completely abolished LTP and left synaptic depression instead (72% ± 5% vs 100%, t(5) = 5.6, p=0.0026, n=6; Figure 3B and C) indicating that, as in some other forms of synaptic potentiation (Grover and Teyler, 1990), Ca²⁺ entry through VGCCs is required for burst reactivation to induce LTP.

Figure 3

Download asset Open asset

Voltage gated calcium channels mediate DA and reactivation-induced plasticity during burst stimulation.

Normalized EPSP slopes were averaged and plotted as a function of time. (A), Application of AP5 after pairing does not prevent DA- and burst-induced potentiation. (B), Application of nimodipine during the burst prevents DA- and burst-induced potentiation leaving synaptic depression. (C), Summary of results. All traces show an EPSP before (1) and 40 min after pairing (2). Plots show averages of normalized EPSP slopes ± SEM.

Figure 3—source data 1 Normalized EPSP slopes of all recorded cells.: https://cdn.elifesciences.org/articles/81071/elife-81071-fig3-data1-v2.xlsx
Download elife-81071-fig3-data1-v2.xlsx

Signaling pathway mediating reactivation-induced LTP

These findings suggest that coincident DA signaling and postsynaptic Ca²⁺ increase enable the potentiation of previously primed synapses. Searching for a potential coincidence detector for DA and intracellular Ca²⁺, we focused on adenylyl cyclases (ACs). They are activated by Gs-coupled dopamine D1/D5 receptor stimulation (Neve et al., 2004), and subtypes AC1 and AC8 are additionally Ca²⁺- stimulated (Wayman et al., 1994; Watson et al., 2000; Ferguson and Storm, 2004). To investigate whether AC subtypes AC1 and/or AC8 are involved in the form of plasticity described here, we tested the induction protocol in AC1/AC8 double knockout (AC DKO) mice (Wong et al., 1999). When postsynaptic burst stimulation in the presence of DA was applied after a negative pairing protocol (Δt = –20 ms; Figure 4Ai) in slices from AC DKO mice, the conversion to potentiation was absent and significantly different from DA- and burst-induced potentiation in slices from wildtype mice (AC DKO, 90% ± 8%, n=6 vs WT, 132% ± 11%, n=8; t(12) = 2.8, p=0.015; Figure 4Bi and Biii), revealing a role for AC1/AC8 as coincidence detector for DA- and Ca²⁺-induced potentiation. In contrast, conventional, DA-independent t-LTP induced by a pre-before-post pairing protocol (Δt = +10 ms; Figure 4Aii) showed significant potentiation comparable to that seen in wildtype mice (AC DKO, 150% ± 19% vs 100%, t(5) = 2.6, p=0.049, n=6; WT, 163% ± 14% vs 100%, t(9) = 4.4, p=0.0015, n=10; Figure 4Bii and Biii).

Figure 4

Download asset Open asset

DA and reactivation-induced potentiation require AC1/AC8 and PKA.

Schematics show the difference between the three induction protocols. Normalized EPSP slopes were averaged and plotted as a function of time. (Ai), DA and burst stimulation after a post-before-pre pairing protocol (Δt = –20 ms) induces potentiation. (**Aii**), A pre-before-post pairing protocol induces t-LTP (Δt = +10 ms). (B), AC DKO mice do not show DA-dependent plasticity with postsynaptic bursts (Bi) but shows conventional t-LTP (**Bii**). Summary of results (**Biii**). (C), Postsynaptic application of protein kinase inhibitor-(6-22)-amide (PKI) prevents DA-dependent plasticity with postsynaptic bursts (Ci) but leaves conventional t-LTP intact (**Cii**). Summary of results (**Ciii**). All traces show an EPSP before (1) and 40 minutes after pairing (2). Plots show averages of normalized EPSP slopes ± SEM.

Figure 4—source data 1 Normalized EPSP slopes of all recorded cells.: https://cdn.elifesciences.org/articles/81071/elife-81071-fig4-data1-v2.xlsx
Download elife-81071-fig4-data1-v2.xlsx

AC activation produces an increase in cyclic adenosine monophosphate (cAMP) which activates protein kinase A (PKA; Sassone-Corsi, 2012). To test whether this signaling cascade is required for reactivation-induced potentiation, we loaded the postsynaptic cell with a PKA blocker, protein kinase inhibitor-(6-22)-amide, through the recording pipette. In this configuration burst stimulation in the presence of DA after the post-before-pre protocol failed to induce synaptic potentiation (69% ± 13% vs 100%, t(5) = 2.4, p=0.064, n=6; Figure 4Ci and Ciii). In contrast, conventional pre-before-post pairing induced significant potentiation, albeit of somewhat reduced magnitude (137% ± 14% vs 100%, t(5) = 2.63, p=0.0463, n=6; Figure 4Cii and Ciii).

The requirement of DA and PKA for burst-induced potentiation is shared with late-phase LTP (Frey et al., 1990; Frey et al., 1993), which requires protein synthesis (Frey et al., 1988). We therefore investigated whether the burst-induced rapid potentiation also requires protein synthesis by delivering the protein synthesis inhibitor anisomycin (AM) to the postsynaptic cell through the recording pipette. We found that, with anisomycin, burst stimulation in the presence of DA no longer induced conversion to potentiation but left a synaptic depression instead, which was significantly different from vehicle control (AM, 52% ± 11% vs 100%, t(5) = 4.5, p=0.0062, n=6; Figure 5Ai and Aiii, red; vs vehicle, 128% ± 17%, n=5; t(9) = 3.9, p=0.0033; Figure 5Ai and Aiii, green). In contrast, conventional t-LTP induced by pre-before-post pairing was unaffected by anisomycin (161% ± 20% vs 100%, t(5) = 3.0, p=0.029, n=6; Figure 5Aii and Aiii, black). Furthermore, the post-before-pre pairing protocol induced t-LTD under these conditions (65% ± 6% vs 100%, t(5) = 5.9, p=0.0019, n=6; Figure 5Aii and Aiii, gray). We confirmed that postsynaptically applied anisomycin did not affect synaptic responses in baseline conditions (95% ± 9.7% vs 100%, t(6) = 0.56, p=0.59, n=7; Figure 5—figure supplement 1A). Furthermore, we compared action potential properties during pairing in cells with anisomycin to cells loaded with vehicle controls. Spike amplitude (AM 112 mV ± 3 mV, Vehicle 111 mV ± 3 mV) and spike width (AM 3.3 ms ± 0.2 ms, Vehicle 3.3 ms ± 0.2 ms) showed no significant differences (amplitude t(10) = 0.09050, p=0.92; width t(10) = 0.1134, p=0.91; Figure 5—figure supplement 1B, C).

Figure 5 with 1 supplement see all

Download asset Open asset

DA and reactivation-induced plasticity requires protein synthesis.

(A), Postsynaptic anisomycin prevents DA-dependent plasticity with postsynaptic burst stimulation (Ai), but leaves conventional t-LTD (**Aii**, gray trace) and t-LTP (black trace) intact. Summary of results (**Aiii**). All traces show an EPSP before (1) and 40 min after (2) pairing. Plots show averages of normalized EPSP slopes ± SEM.

Figure 5—source data 1 Normalized EPSP slopes, spike amplitudes and spike width of all recorded cells.: https://cdn.elifesciences.org/articles/81071/elife-81071-fig5-data1-v2.xlsx
Download elife-81071-fig5-data1-v2.xlsx

Burst-dependent plasticity increases specificity in reinforcement learning models

These experimental results show that, after a priming event, burst reactivation in the presence of DA induces a rapid form of protein synthesis-dependent LTP. This mechanism would ensure that only salient neuronal activity induces long-term changes in the network. We implemented this synaptic learning rule in a computational model to explore how such plasticity would control learning in a feedforward artificial neural network resembling hippocampal CA1. The learning rule states that the change in synaptic weights $Δ w$ between input and output neurons ( $i n p \to o$ ) depends on an eligibility trace e (set during the initial priming event), the reinforcement signal (dopamine d), and reactivation (bursting activity b).

{{∆ w}_{i n p \to o} = α_{L T P} e d b}_{}

The parameter α_LTP is the learning rate. When there is no DA or bursts during the trial, the rule is updated to result in depression proportional to the eligibility trace with a learning rate α_LTD (see Methods). When using a standard reinforcement learning (RL) rule, which does not depend on burst reactivation, all previously primed synapses are potentiated after receiving the DA signal (Figure 6A). Thus, the global neuromodulatory signal in traditional RL models provides limited information. In contrast, when applying the burst-dependent learning rule (Equation 1) to the network, in which potentiation of primed synapses depends on both the reward signal and reactivation, a selected subset of inputs becomes potentiated, while inputs on non-reactivated neurons remain depressed (gray) (Figure 6B). During DA modulation, information is allocated to primed synapses by reactivation of the postsynaptic neuron, and the broader computational implications of this learning rule depend on the control of postsynaptic neuronal bursting activity. First, it is possible that neurons are recruited to a new memory trace based on their relative neuronal excitability before the training session as suggested by the memory allocation hypothesis (Yiu et al., 2014). According to this scenario, the most excitable cells would be the most likely to show action potential bursts during reactivation and, therefore, show synaptic potentiation. Alternatively, there is evidence for replay of prioritized experience (Igata et al., 2021), suggesting that cells encoding the most salient events preceding the reward would reactivate, and thereby determine which set of cells would show potentiation during reward (Figure 6B). Assuming prioritized experience reflects experience relevant to the reward, this could help credit assignment in the network. In addition, because of the exclusive requirement of postsynaptic activity for potentiation to occur, this mechanism offers another intriguing possibility, namely that other inputs active at the reward location carries additional information about the nature of the reward, e.g., food, or the reward location, such as the presence of specific landmarks, which could elicit postsynaptic bursting activity. Under this hypothesis, neuronal reactivation would serve to associate a specific outcome to the priming event. When different instructive inputs induce bursting each in distinct subsets of neurons during reward, selective increases in synaptic weights would not only enable the encoding of reward but also distinguish between different rewards (Figure 7A). Finally, we explore how the learning rule performs in a network supervised by feedback synaptic input to strengthen synapses onto specific neurons encoding reward-related features. By allowing feedback input to assign which part of the network is responding, the burst-dependent learning rule enables the network to selectively learn relevant information (Figure 7B), resulting in potentiation of those synapses (magenta and cyan in Figure 7Biii) that are active temporally separated, but less when simultaneous active (Figure 7Biv). This provides a mechanism to associate temporally separated, reward-relevant information in a neuronal network. The burst-dependent learning rule provides the network with a gating mechanism for memory allocation to an engram. Thus, burst-induced plasticity reduces the number of neurons encoding the reward location and other reward-related information, increasing the specificity of synaptic memory in a neuronal network.

Figure 6

Download asset Open asset

Dopamine-dependent burst-induced plasticity rule reduces the number of neurons in the network encoding a memory.

(**Ai-iii**), Standard reinforcement learning rule shows reward associating inputs 1–10 (blue) with potentiation of all synaptic weights (blue). (**Bi-iii**), Burst-dependent potentiation reduces the number of neurons encoding the memory, leading to potentiation of synapses exclusively onto the most excitable burst-firing neurons 3, 6 (blue).

Figure 6—source data 1 Raster plot data and synaptic weights.: https://cdn.elifesciences.org/articles/81071/elife-81071-fig6-data1-v2.csv
Download elife-81071-fig6-data1-v2.csv

Figure 7

Download asset Open asset

Dopamine-dependent burst-induced plasticity rule enables reinforcement learning (RL) models to encode a specific salient event.

(**Ai-iii**), Instructive RL rule allows two inputs that code for different information to store the memory in separate sets of neurons, thus encoding not only the reward, but also other reward-relevant features 3, 6 (magenta, cyan). (**Bi-ii**), A supervised network enables burst-eliciting feedback synaptic input to assign credit to select synapses in the network to encode a desired reward identity. (**Biii**) Time-dependent lateral inhibition at the output neurons suppress non-relevant information. When only one of the inputs is active, the animal can learn two different memories over time in neurons 3, 6 (magenta, cyan). (**Biv**) When both inputs are active at the same time they compete with each other, and synapses onto these neurons (magenta, cyan) are less potentiated.

Figure 7—source data 1 Raster plot data and synaptic weights.: https://cdn.elifesciences.org/articles/81071/elife-81071-fig7-data1-v2.csv
Download elife-81071-fig7-data1-v2.csv

Increased calcium responses and spatial information in reactivated CA1 place cells

Based on these results we predicted that, when an animal navigates toward a reward, the previously reactivated hippocampal neurons would be more strongly activated than non-reactivated neurons. To test this prediction, we monitored calcium transients in hippocampal excitatory cells with a head-mounted microscope (Ghosh et al., 2011) while mice navigated on a ‘cheeseboard’ maze (Dupret et al., 2010) with two reward locations, one new to the animal and one previously learnt (Figure 8A). We defined neuronal reactivation as activity during immobility in neurons previously active during locomotion. Cells that were active when mice moved towards the reward locations were classified as reactivated at reward if they were active again during immobility when mice consumed reward (Foster and Wilson, 2006; O’Neill et al., 2006; Csicsvari et al., 2007) and non-reactivated if no further calcium event was detected after they had reached the reward. Of the cells that were active on the maze in a given trial, 44 ± 1% were reactivated at either or both of the reward locations. There was no detectable difference between the number of cells reactivating at the previously learnt or new reward location (Figure 8B). Less frequent were the reactivations during immobility at non-rewarded locations where 15 ± 2% of cells reactivated (Figure 8B).

Figure 8 with 2 supplements see all

Download asset Open asset

Cells that reactivated at reward location have higher activity peaks than non-reactivated cells.

(A), Running path of a mouse (left) and calcium traces of example cells (right) in a representative single trial. Mice ran toward two reward locations: one previously learned (R₁) and one new (R₂). Calcium traces for two cells that reactivated at the reward location and two that did not. Arrows mark reactivations. Vertical line marks time of arrival at R₂. (B), Percentage of cells that reactivated at one or both of the reward locations or other locations. Points mark percentage in single trials. Box-and-whisker plots: median, 25 and 75th percentile, and extreme values. n=444 trials. (C), Traces centered on activity peaks for two example cells: one that reactivated at the new reward location in trial 1 (top) and one that did not reactivate in any of the trials (bottom). (D), Histogram of activity peaks normalized to maximum AUC value compared between cells that reactivated for the first time in the preceding trial and cells that did not reactivate in any of the trials. Triangles mark mean values. Permutation tests for repeated measures ANOVA, model for reactivations at the reward locations: significant effect of trial (F_{(6, 96)}=4.51, p<0.001) and reactivation (F_(1,16) = 10.15, p=0.006), model for reactivations at non-rewarded locations: significant effect of trial (F_{(6, 96)}=7.85, p=0.001) and non-significant of reactivation (F_{(1, 16)}=0.56, p=0.47). n=6785 cell trials from 1405 cells. (E), The duration from activity during locomotion to the time of first reactivation did not correlate with activity peaks in the following trial, suggesting a long-lasting eligibility trace.(F), Mean calcium peak in each trial for cells that reactivated in trial 1 or trial 4 compared to cells that did not reactivate in any trial. Ribbons extend +/- 1 SEM. Cells that reactivated in trial 1 had significantly higher normalized calcium peaks in all trials. Cells that reactivated for the first time in trial 4 had significantly higher normalized calcium peaks in trials 4, 5, 7 and 8 but not in trial 6 and trials before the reactivation. Permutation t-tests with Benjamini-Hochberg correction for multiple comparisons and between-animal random effects. n=6721 cell trials from 681 reactivated in trial 1 cells, n=769 from 96 reactivated in trial 4 cells, and n=5980 from non-reactivated 607 cells. **p<0.01.

To compare the strength of neuronal activation, we measured area-under-curve (AUC) of calcium events occurring before and after the reactivation (Figure 8C). Cells with the largest activity peaks during locomotion were the most likely to reactivate at the reward location. Following their first reactivation, they had larger activity peaks than the previously non-reactivated cells (0.54 ± 0.01 vs 0.42 ± 0.01 of the cell’s max AUC, F_(1,16) = 10.12, p=0.006, Figure 8D). The effect on activity peaks was specific to reward locations, and reactivation at other locations did not affect activity peaks in the following trials (F_{(1, 16)}=0.56, p=0.47, Figure 8D). The effect of reactivation was independent of the cell event rates: of two cell groups with matching event rates in a given trial, the one whose cells were reactivated at reward for the first time, had larger peaks in the following trial (F_{(1, 16)}=7.54, p=0.02, Figure 8—figure supplement 1A). The effect of reactivation at reward had two other similarities with the synaptic plasticity experiments: a long-lasting eligibility trace and persistence of the change. To test the first, we confirmed that the activity peaks in the trial following the reactivation did not correlate with the time from the activity during locomotion to the reactivation (Figure 8E). To test the persistence of the change, we confirmed that the reactivated cells maintained higher calcium activity peaks than non-reactivated cells throughout the later trials (Figure 8F, Figure 8—figure supplement 1B). Typically, the peaks during locomotion increased in the trial immediately before the first reactivation at reward (Figure 8—figure supplement 1B), suggesting that more excited cells are the ones that undergo the change. The calcium activity peaks we report could be affected by photobleaching of the GCaMP6f sensor. As the reactivated cells were more active than non-reactivated cells (Figure 8—figure supplement 1A), we would expect them to be more affected by photobleaching, but nevertheless, the magnitude of their calcium transients remained higher than in the non-reactivated cells.

60% of the reactivated cells (n=971 from 7 mice) and 42% of the non-reactivated cells (n=284 from 7 mice) showed activity at specific locations and were classified as place cells (see Methods). Both place cells and non-place cells showed higher calcium peaks following reactivation and the increase in place cells was not significantly different from that in non-place cells (F_(1,16) = 0.11, p=0.75, Figure 8—figure supplement 1C, D). To investigate any learning-induced changes in place maps we assessed how the location-averaged activity in place cells changed with reactivation from the first to the second half of the trials (change from trials 1–4 to trials 5–8, Figure 8—figure supplement 2A). The place map peaks increased significantly more in the place cells that reactivated in trial 4 or later than in the non-reactivated place cells (Figure 8—figure supplement 2B, place map peak ratio 1.4 ± 0.1 vs 1.1 ± 0.1, n=194 reactivated and n=306 non-reactivated place cells; permutation test for repeated measures ANOVA: F_{(1, 16)}=11.3, p=0.008). The change in spatial information also significantly differed between the two groups (Figure 8—figure supplement 2C, change by 0.011 ± 0.016 to −0.047 ± 0.014 a.u.; permutation test for repeated measures ANOVA: F_{(1, 16)}=4.4, p=0.048), but there was no significant difference in place map stability (Figure 8—figure supplement 2E, correlation of 0.45 ± 0.02 vs 0.47 ± 0.02; permutation test for repeated measures ANOVA: F_{(1, 16)}=0.29, p=0.59). The time from the activity during locomotion to the reactivation did not correlate with a change in spatial information (Figure 8—figure supplement 2D). Reactivated place cells conveyed more spatial information in late trials (trials 5–8) compared to place cells that did not reactivate (0.15 ± 0.01 vs 0.12 ± 0.01 a.u.; permutation test for repeated measures ANOVA: F_{(1, 16)}=4.4, p=0.055). This learning-associated increase in calcium response and spatial information supports a reactivation-dependent LTP-like mechanism in vivo (Cacucci et al., 2007).

Discussion

In summary, we investigated the effects of postsynaptic neuronal reactivation on hippocampal synaptic plasticity, reinforcement learning, and spatial coding. We found that postsynaptic burst reactivation of CA1 pyramidal cells in the presence of the reward signal DA rapidly potentiates synapses that have previously undergone a spike timing-dependent priming protocol, providing direct evidence for reactivation-induced synaptic plasticity. A computational model showed how this learning rule increases specificity in reinforcement learning models. Recordings from freely moving mice showed that neurons that reactivated at reward locations had enhanced CA1 place cell calcium signals and carried more spatial information than cells that did not reactivate.

The results suggest that reactivation-induced plasticity is mediated by two sequential coincidence detectors: postsynaptic NMDARs detecting coincident pre- and postsynaptic activity and AC1/AC8 as coincidence detector of DA and reactivation-induced Ca²⁺ increase. Although we used a spike pairing protocol, we cannot exclude the possibility that activation of postsynaptic NMDA receptors without postsynaptic action potentials would be sufficient to set the eligibility trace.

AC1/AC8 is synergistically activated when the two signals, Gs-coupled dopamine D1/D5 receptor activation and Ca²⁺ influx, occur at the same time (Wayman et al., 1994; Watson et al., 2000; Ferguson and Storm, 2004; Neve et al., 2004). The time course of DA signaling depends on brain area, firing mode (tonic or phasic) of dopaminergic cells, DA release and diffusion as well as time course of intracellular signaling pathways (Liu et al., 2021). Recent developments of fluorescent DA sensors (Sun et al., 2018; Sun et al., 2020; Elizarova et al., 2022) would enable monitoring the precise time course of DA in the hippocampus in future studies. Moreover, to investigate the precise timing requirements for DA-dependent reactivation-induced plasticity further, uncaging of caged DA or optogenetically-induced DA release would be suitable approaches for temporal control of the DA transient.

Our experiments were carried out at 2 mM external calcium concentration, which is a standard calcium concentration used in most ex vivo plasticity experiments, but above the reported ionic calcium concentration in rat and human cerebrospinal fluid (Jones and Keep, 1988; Forsberg et al., 2019). Unfortunately, the extracellular calcium concentration at synaptic sites is not known, as discussed elsewhere (Lopes and Cunha, 2019). It was recently reported that burst pairing, but not pairings of single pre- and postsynaptic action potentials, induces synaptic plasticity at 1.3 mM extracellular calcium in rat hippocampal slices (Inglebert et al., 2020). It will be interesting to investigate whether, even if single-spike pairing might not induce plasticity at low calcium concentrations, it would still be sufficient for the initial priming of those synapses. If that were the case, this would set the conditions to enable reactivation-induced plasticity, which relies on bursts, and would hold also in low-calcium conditions. The signaling cascade leading to synaptic potentiation involves postsynaptic PKA and protein synthesis, which are not required for conventional early LTP (Park et al., 2014). Traditionally, LTP has been classified into early- and late-phase LTP (Frey et al., 1988; Frey et al., 1990), and it has been reported that dopaminergic signaling is required for maintenance of late-phase LTP (Frey et al., 1990; Huang and Kandel, 1995; Matthies et al., 1997). The DA-dependent form of plasticity we describe here shares properties with ‘late-phase’ LTP, including a role of postsynaptic action potentials during induction (Dudek and Fields, 2002) and a requirement of protein synthesis for expression (Frey et al., 1988; Huang and Kandel, 1995). However, it is remarkably fast, suggesting a dissociation between different signaling pathways, rather than different temporal phases of LTP. The mechanism we describe here would be compatible with some of the key concepts of the ‘revised’ synaptic tagging hypothesis (Redondo and Morris, 2011). Specifically, our findings strongly support the view that the fate of a memory is not determined at the time of encoding. This is based on the finding that plasticity induction can lead to two events: (1) expression of LTP or LTD, and (2) setting of a synapse-specific eligibility trace (‘tagging’) which allows modulation by protein synthesis-dependent signaling. A major difference is that our findings highlight a role of DA as a key modulator of the eligibility trace. The mechanistic basis for the surprising involvement of protein synthesis in this rapidly induced form of plasticity remains to be investigated.

It has been reported that anisomycin can potentiate JNKs (Iordanov et al., 1997). We cannot exclude the possibility that the drug may have affected intracellular signaling cascades that interfere with the plasticity signaling pathway described here.

Our experimental findings uncover a synaptic learning rule that could support a two-stage model of hippocampal memory formation (Buzsáki, 1989), in which eligibility traces are laid down during hippocampal theta activity with subsequent postsynaptic burst reactivation during sharp wave-ripples inducing LTP at those synapses.

We considered the activity during navigation on the maze as the animal approaches the reward resembling the STDP priming protocol. Substantial evidence supports a role of NMDAR-dependent STDP in the formation of place fields during navigation (Mehta, 2015; Moore et al., 2021). It has been postulated that both LTP and LTD are involved in place field formation. This was based on the observation that place fields shift backward with experience (Mehta et al., 1997), and a computational model predicted that without LTD place field broadening would occur. Thus LTP is required when entering the place field, and LTD when the animal exits the place field (Mehta et al., 2000). This is specific to navigation, as opposed to just walking on a linear track without task, and place field plasticity is predictive of navigational performance (Moore et al., 2021).

Exploring possible computational implications of this synaptic learning rule, we first tested it with a neural network in which the most excitable cells are reactivated. The memory allocation hypothesis suggests that learning triggers a temporary increase in neuronal excitability, enabling the linking of individual memories acquired close in time (Silva et al., 2009; Cai et al., 2016; Sehgal et al., 2018). We found that our learning rule selectively strengthens the reactivated synapses, linking together the memory allocation hypothesis with burst reactivation-induced plasticity. Interestingly, it was recently reported that DA released by locus coeruleus cells projecting to dCA1 has a key permissive role in contextual memory linking (Chowdhury et al., 2021). Moreover, the rule could also accommodate a temporally discontiguous instructive learning signal or a specific supervisory feedback signal. Thus, it is possible that DA serves as a scalar reward signal which combines with a vectorial representation of the reward identity which triggers the reactivation of a specific subset of neurons. Thus, it was suggested that the direct pathway from the entorhinal cortex could provide an instructive signal to generate accumulation of CA1 place cells at the reward location (Grienberger and Magee, 2021). An attractive possibility is that a downstream subset of neurons active during navigation serves as an instructive input onto upstream neurons during reward. This would be consistent with a dual role of DA as reinforcement signal and enhancer of reverse replay (Ambrose et al., 2016), establishing a predictive chain of potentiated synapses toward the rewarded outcome reminiscent of the successor representation (Dayan, 1993; Stachenfeld et al., 2017). It was suggested in a computational study that feedback regulation of synaptic plasticity by bursts in higher hierarchical circuits can coordinate lower-level connections (Payeur et al., 2021). Our results reveal a possible biological substrate to support such a mechanism and re-emphasize the importance of bursting activity for synaptic plasticity (Lisman, 1997; Pike et al., 1999).

The hippocampus receives dopaminergic input from two main sources, the ventral tegmental area (Scatton et al., 1980; Gasbarri et al., 1994), signaling reward, and the locus coeruleus (Smith and Greene, 2012; Kempadoo et al., 2016), thought to signal novelty (Takeuchi et al., 2016; Wagatsuma et al., 2018), but more recently implicated also in spatial reward learning (Kaufman et al., 2020). In a reward location learning task, we found that calcium signals during locomotion were higher in CA1 principal cells that reactivated at reward location. The signal increased in the trial preceding the first reactivation, following which the signal remained higher in all subsequent trials. Moreover, place cells that reactivated showed higher spatial information than non-reactivating place cells in late trials. The overall stability of calcium signals and spatial information is broadly consistent with earlier reports using one-photon imaging in freely moving mice (Ziv et al., 2013). Although it was recently reported in head-fixed mice that reactivation increases long-term stability of place cells that have fields distant from the reward after several days (Grosmark et al., 2021), we did not see a significant increase in place cell stability after a single reactivation of the cell at the reward location. Irrespectively, the difference in both calcium signal and spatial information in reactivating vs non-reactivating cells is suggestive of plasticity related to reactivation events specific to reward location. More work will be required to establish under what conditions the novel burst-induced potentiation mechanism is engaged during hippocampus-dependent learning and memory.

Methods

Key resources table

Reagent type (species) or resource	Designation	Source or reference	Identifiers	Additional information
Strain, strain background (Mus musculus)	Wild-type C57BL/6 J, Mus musculus	Harlan, Bicester, UK or Central Animal Facility, Physiological Laboratory, Cambridge University		Age range used for slice preparation: 12–19 days. Females and males were used.
Strain, strain background (Mus musculus)	Adenylate cyclase double knock-out (AC DKO), C57BL/6 J, Mus musculus	Mouse line generated by Hongbing Wang and imported from Michigan State University, MI, US		Age range used for slice preparation: 12–19 days. Females and males were used.
Strain, strain background (Mus musculus)	Thy1 – GCaMP6f, C57BL/6 J, Mus musculus	Jax Ref: Dana et al., 2014	024276	Only males were used.
Chemical compound, drug	Dopamine hydrochloride	Sigma–Aldrich	H8502	100 μM
Chemical compound, drug	D-AP5	Tocris Bioscience	0106	100 µM
Chemical compound, drug	Nimodipine	Tocris Bioscience	0600/100	10 µM
Chemical compound, drug	PKA inhibitor fragment (6-22) amide	Tocris Bioscience	1904/1	1 µM
Chemical compound, drug	Anisomycin	Tocris Bioscience	1290/10	0.5 mM
Chemical compound, drug	MK801	Tocris Bioscience	0924	1 mM
Software, algorithm	Igor Pro 6.37	WaveMetrics
Software, algorithm	Prism 8.2.0 (435)	Graphpad
Software, algorithm	Matlab R2021a	MathWorks
Software, algorithm	CaImAn software (version 1.8.5, Python)	Giovannucci et al., 2019

Mice

Experimental procedures and animal use were performed in accordance with UK Home Office regulations of the UK Animals (Scientific Procedures) Act 1986 and Amendment Regulations 2012, following ethical review by the University of Cambridge Animal Welfare and Ethical Review Body (AWERB). All animal procedures were authorized under Personal and Project licences held by the authors.

Mice were housed on a 12 hr light/dark cycle at 19–23 °C and were provided with food and water ad libitum. Experiments were carried out on wildtype C57BL/6 J mice (Harlan, Bicester, UK or Central Animal Facility, Physiological Laboratory, Cambridge University), and adenylate cyclase double knockout (AC DKO) mice, which have the genes for both AC1 and AC8 deleted globally. This mouse line was generated as described previously²⁵ and was imported from Michigan State University, MI, US. For in vivo experiments, seven adult males Thy1 – GCaMP6f transgenic mice were used (Dana et al., 2014) (Jax: 024276). Mice were housed with 2–4 cage-mates in cages with running wheels.

Electrophysiology

Slice preparation

Mice of both sexes at postnatal day (P) 12–19 were used in this study. Mice were anesthetized with isoflurane (4% isoflurane in oxygen) and decapitated. The brain was rapidly removed and immersed in ice-cold artificial cerebrospinal fluid (ACSF) containing (in mM): 126 NaCl, 3 KCl, 26.4 NaH₂CO₃, 1.25 NaH₂PO₄, 2 MgSO₄, 2 CaCl₂, and 10 glucose (pH 7.2, 270–290 mOsm/L). The ACSF solution was continuously bubbled with carbogen gas (95% O₂, 5% CO₂). Horizontal slices (350 μm thick) were sectioned with a vibrating microtome (Leica VT 1200 S, Leica Biosystems, Wetzlar, Germany). The slices were then incubated for at least 60 min in ACSF at room temperature in a submerged-style storage chamber before recording. Slices were used for 1–7 hr following sectioning.

Whole-cell patch clamp recording

For recordings, individual slices were transferred to an immersion-type recording chamber and perfused with ACSF (2 ml/min) at 24–26 °C. Neurons were visualized and selected using infrared differential interference contrast (DIC) microscopy using a 40 X water-immersion objective. The hippocampal subfields were visually identified and whole-cell patch-clamp recordings were performed on CA1 pyramidal neurons. For stimulation of Schaffer collaterals, monopolar stimulation electrodes were placed in stratum radiatum. Test and control pathway electrodes were placed at the same distance (>100 µm) from and either side of the recorded neuron. Patch pipettes (pipette resistance 4–7 MΩ) were pulled from borosilicate glass capillaries (0.68 mm inner diameter, 1.2 mm outer diameter) using a P-97 Flaming/Brown micropipette puller (Sutter Instruments Co., Novato, California, USA). Pipettes were filled with a solution containing (mM): 110 potassium gluconate, 4 NaCl, 40 HEPES, 2 ATP-Mg, 0.3 GTP (pH 7.2–7.3, 270–285 mOsm/L). The liquid junction potential was not corrected.

All experiments were performed in current-clamp mode. Cells were accepted for the experiment if their resting membrane potential was between −55 and −70 mV. The membrane potential was held at −70 mV throughout the recording by direct current application via the recording electrode. Before the start of each recording, all cells were tested for regular spiking responses to positive current steps—characteristic of pyramidal neurons.

Stimulation protocol

Excitatory postsynaptic potentials (EPSPs) were evoked alternately in two input pathways (test and control) by direct current pulses at 0.2 Hz (stimulus duration 50 μs) through metal stimulation electrodes. Control pathways were used in all experiments to ensure stability control (not always shown). The stimulation intensity was adjusted (100 μA– 500 µA) to evoke an EPSP with peak amplitude between 3 and 8 mV. After a stable EPSP baseline period of at least 10 min, STDP was induced in the test pathway by repeated pairings of single evoked EPSPs and single postsynaptic action potential elicited with the minimum somatic current pulse (1–1.8 nA, 3ms) via the recording electrode. Pairings were repeated 100 times at 0.2 Hz. Spike-timing intervals (Δt in ms) were measured between the onset of the EPSP and the onset of the action potential.

Alternate stimulation of EPSPs was resumed immediately after the pairing protocol and monitored for at least 40 min. For the burst stimulation protocol, stimulation of EPSPs was not resumed for an additional 10 min, and at the end of that period, six bursts, each of five action potentials at 50 Hz, were elicited with an inter-burst interval of 0.1 Hz by somatic current pulses (1.8 nA, 10 ms) via the recording electrode. In a subset of experiments, only five bursts were applied, which led to potentiation of a similar magnitude. Immediately after the bursts, stimulation of EPSPs was resumed and monitored for at least 30 min, however, a small subset of recordings were stopped at 28 min. Presynaptic stimulation frequency to evoke EPSPs remained constant throughout the experiment. The unpaired pathway served to verify input-specificity and as a stability control. The burst stimulation protocol is summarized in Figure 1a (top).

Drugs

Drugs were bath-applied to the whole slice through the perfusion system by dilution of concentrated stock solutions (prepared in water or DMSO) in ACSF, or by adding the drugs to the patch pipette solution when it was applied intracellularly to the postsynaptic cell only. If the drug was not water-soluble, vehicle control experiments were carried out. For each set of recordings, interleaved control and drug conditions were carried out and were pseudorandomly chosen. The following drugs were used in this study: 100 μM dopamine hydrochloride (Sigma–Aldrich, Dorset, United Kingdom), 100 μM D-AP5 (Tocris Bioscience, Bristol, United Kingdom), 10 µM nimodipine (Tocris Bioscience), 1 µM PKA inhibitor fragment (6-22) amide (Tocris Bioscience), 0.5 mM anisomycin (stock solution in EtOH; Tocris Bioscience), and 1 mM MK801 (Tocris Bioscience).

Data acquisition and data analysis of slice recordings

Data were collected using an Axon Multiclamp 700B amplifier (Molecular Devices, Sunnyvale, California, USA) and filtered at 2 kHz. Data were acquired and digitized at 5 kHz using an Instrutech ITC-18 A/D interface board (Instrutech, Port Washington, New York, USA) and custom-made acquisition procedures in Igor Pro (WaveMetrics, Lake Oswego, Oregon, USA).

All experiments were carried out in current clamp (‘bridge’) mode, and only cells with an initial series resistance between 9 and 16 MΩ were included. Series resistance was compensated for by adjusting the bridge balance, and data were discarded if series resistance changed by more than 30%. Offline analyses were done using custom-made procedures in Igor Pro. EPSP slopes were measured on the rising phase of the EPSP as a linear fit between the time points corresponding to 25–30% and 70–75% of the peak amplitude. For statistical analysis, the mean EPSP slope per minute of the recording was calculated from 12 consecutive sweeps and normalized to the baseline (each data point in source data files is the mean of 12 averaged EPSPs). Normalized EPSP slopes from the last 5 min of the baseline (immediately before pairing) and from the last 5 min of the recording were averaged. The magnitude of plasticity, as an indicator of change in synaptic weights, was defined as the average EPSP slope after pairing expressed as a percentage of the average EPSP slope during baseline.

Statistical analysis of slice recordings

Statistical comparisons were performed using one-sample two-tailed, paired two-tailed, or unpaired two-tailed Student’s t-test, with a significance level of α=0.05. Data are presented as mean ± SEM. Significance levels are indicated by *p<0.05, **p<0.01, ***p<0.001. All datasets passed the test for normality using the Shapiro-Wilk test (α=0.05).

Computational modeling

We simulated a set of n=10 output neurons, which each received input from 10 input neurons. When an instructive input was added, output neurons additionally received input from two out of 10 instructive neurons. Each output neuron projected uniquely onto one readout neuron which again projected back to the output neuron in a one-to-one mapping. Each output neuron was modeled as an Integrate-and-Fire neuron and 100 trials were simulated, where the voltage v is described by

τ_{v} \frac{d v}{d t} = - v + w_{i n p \to o} I_{i n p} + w_{i n s t \to o} I_{i n s t} + w_{r e a d \to o} I_{r e a d} + I_{i n t r i n s i c}

where $τ_{v} = 10 m s$ is the membrane time constant, $I_{i n p}$ are the spike trains of the input neurons, $I_{i n s t}$ are the spike trains from the instructive neurons, $I_{r e a d}$ are the spike trains from the read-out neurons, $w_{i n p \to o}$ are the weights from the input to the output neurons, $w_{i n s t \to o} = 1 / N$ are the weights from the instructive neurons to the output neurons, and $w_{r e a d \to o} = 1$ are the weights from the readout to the output neurons. In addition, each neuron was receiving an intrinsic current $I_{i n t r i n s i c} = α_{i n t r i n s i c} η$ , where $η$ is simply white noise drawn from a uniform distribution between 0 and 1 and $α_{i n t r i n s i c}$ is the magnitude of the current ( $α_{i n t r i n s i c}$ was applied to all neurons in Figure 6B, while in the other configurations no intrinsic current was applied). When the voltage crosses the firing threshold = 0.4, the neuron is reset to the resting potential v=0. Each read-out neuron was also modeled as an Integrate-and-Fire neuron

τ_{v} \frac{d v_{r e a d}}{d t} = - v_{r e a d} + w_{o \to r e a d} I_{o} + I_{}

where $I_{o}$ are the spike-trains of the output neurons, $w_{o \to r e a d}$ = $1 / N$ are the weights from the output to the readout neurons, and $I_{}$ are the spike trains from the supervised neurons. Similarly, if the voltage crossed the firing threshold = 0.4, the neurons emit a spike and the voltage is reset to 0. $w_{i n p \to o}$ were plastic under the following rule. Every time the input neurons are firing (at $t^{p r e}$ ), they are leaving a trace $x^{p r e},$ $τ_{S T D P} \frac{d x^{p r e}}{d t} = - x^{p r e} + δ (t - t^{p r e}),$ where $τ_{S T D P} = 10 m s$ is the trace time constant. Similarly, every time the output neuron spikes (at $t^{p o s t}$ ), it is leaving a trace $x^{p o s t},$ $τ_{S T D P} \frac{d x^{p o s t}}{d t} = - x^{p o s t} + δ (t - t^{p o s t})$ .

The eligibility trace e is described as

τ_{e} \frac{d e}{d t} = - e + x^{p r e} x^{p o s t}

where $τ_{e} = 10$ min is the eligibility time constant. This value was based on the experimental protocol, where burst reactivation was applied 10 min after priming. The $w_{i n p \to o}$ are potentiated if there is an eligibility trace together with dopamine and a postsynaptic burst:

w_{i n p \to o} = α_{L T P} d e b

where $α_{L T P} = 0.0002$ is the learning rate, d=0 if there is no dopamine, and d=1 if there is dopamine.

The postsynaptic burst b is detected as follows. We first computed a trace $x^{b u r s t}$ as

τ_{b u r s t} \frac{d x^{b u r s t}}{d t} = - x^{b u r s t} + δ (t - t^{p o s t}),

where $τ_{b u r s t} = 5 m s$ is the burst time constant. We set b=1 if $x^{b u r s t} > b u r s t_{t h r e s h o l d}$ , to 0 otherwise, with ${b u r s t}_{t h r e s h o l d}$ = 1.1. If there was no dopamine nor bursts during the whole trial, then the updated rule resulted in a depression

w_{i n p \to o} = α_{L T D} e

where $α_{L T D} = α_{L T P}$ /500 is the learning rate for the depression. Weights are bound to stay positive. The weights $w_{i n p \to o}$ are initialized to 1 /N. We simulated the network for 2 s in Figure 6 and 2.5 s in Figure 7. The input neurons were firing Poisson statistics at 20 Hz for the first 500 ms of the trial (pairing phase). Dopamine was present during the last 500 ms of the trials in Figure 6 and the last 1 s in Figure 7. The network was simulated for 100 trials (error bars are standard deviations).

In Figure 6A, a standard reinforcement learning rule was used: $w_{i n p \to o} = w_{i n p \to o} + α_{R L} d e$ , where $α_{R L} = α_{L T P}$ /50. In Figure 6B, a burst dependent rule was used. Output neurons 3 and 6 had an increased excitability. We modeled that by decreasing the firing threshold to 0.25. All neurons also received an intrinsic current with $α_{i n t r i n s i c} = 0.06$ during the last 500 ms of each trial. In Figure 7A, the instructive neuron 3 fired Poisson statistics at 500 Hz from time 1.5–2 s and the instructive neuron 6 fired Poisson statistics at 500 Hz from time 2–2.5 s. In Figure 7B ii/iii, the read-out neuron 3 received additional Poisson inputs at 60 Hz from time 1.5–2 s and the read-out neuron 6 from time 2–2.5 s. For Figure 7B iv, however, the read-out neurons 3 and 6 received their additional Poisson inputs at the same time, from 1.5–2.5 s. In Figure 7B, we added lateral inhibition, where each read-out neuron received an inhibitory filtered version (with a time constant of 50 ms) of the spike trains of other read-out neurons, with weights of –1 /N (without self-connection).

Behavior

Surgery

Mice underwent two surgeries: the first one to implant a GRIN lens directly above the cells of interest, and the other to fix an aluminum baseplate above the GRIN lens for later attachment of the miniature microscope. The procedures followed the protocol described before by Resendez et al., 2016. Surgeries were carried out following minimal standard for aseptic surgery. Analgesic (Meloxicam, 2 mg^.kg-1 intraperitoneal) was administered 30 min prior to surgery. Mice were anesthetized with isoflurane (5% induction, 1–2% maintenance, Abbott Ltd, Maidenhead, UK) mixed with oxygen as carrier gas (flow rate 1.0–2.0 L^.min^-1) and placed in a stereotaxic frame (David Kopf Instruments, Tujunga, CA, USA). The skull was exposed by making skin incision and Bregma and Lambda were aligned horizontally. A 1.5–2 mm-wide craniotomy was drilled above the implantation site. The brain tissue above the implantation site was aspirated. Buffered ACSF was applied throughout the aspiration to prevent desiccation of the tissue. A GRIN lens (1 mm diameter, 4.3 mm length, 0.4 pitch, 0.50 numerical aperture, Grintech) was stereotaxically lowered at coordinates –1.75 AP, 1.75 ML, 1.35–1.40 DV (in mm from Bregma) and fixed to the skull surface with ultraviolet-light curable glue (Loctite 4305) and further fixed with dental adhesive (Metabond, Sun Medical) and dental acrylic cement (Simplex Rapid, Kemdent). A metal head bar was attached to the cranium using dental acrylic cement for head-fixing the animal during the microscope mounting.

If the GCaMP6f expression was visible in the implanted mouse, 4 weeks later the animals were anesthetized for the purpose of attaching a baseplate for the microscope above the top of the GRIN lens. The baseplate was cemented into place and the miniscope was unlocked and detached from the baseplate.

Behavioral learning task

The mice performed a rewarded spatial navigation task on a round-shaped maze (cheeseboard maze; Dupret et al., 2010). The 120 cm diameter cheeseboard had a total of 177 evenly spaced wells. The rewarded wells were baited with ∼100 μL of condensed milk mixed 1:1 with water.

For the first three days, the mice foraged for rewards baited in randomly selected wells. A different, random set of wells was baited in each trial. Next, we performed a spatial learning task. The mice had to learn two locations with baited wells. The baited wells had fixed locations that were at least 40 cm apart, chosen pseudo randomly for each mouse. Mice started the trial in one of the three locations on the maze: south, east or west. The maze was rotated and wiped with a disinfectant (Dettol) in-between the trials to discourage the use of intra-maze cues. Landmarks of black and white cues were installed on the walls surrounding the maze. The trials were terminated once the mice ate both rewards or after 300 s. Each learning day consisted of 8 trials with 2–4-minute-long breaks between the trials. To minimize the effects of the novel environment and task structure, we analyzed the neural activity on the first learning day after one of the previous reward locations was moved. After the 5-day-long learning, the memory retention was tested on the next day in a 4 to 5-minute-long unbaited trial. Following the learning and testing of the memory for the first set of locations, we translocated one of the reward locations. The new location was a pseudo-randomly chosen to be at least 40 cm away from the current and previous reward locations. The learning of the new set of locations was performed over two days and tested in an unbaited trial as described above.

The trials were recorded with an overhead webcam video camera. The video was recorded at 24 Hz frame rate. The mice body location was tracked with DeepLabCut software (Mathis et al., 2018), a custom-written software was written to map the mouse coordinates to relative location on the maze. The extracted tracks were smoothed with a Gaussian kernel. Periods of running were identified when velocity of the mouse smoothed with a moving average 0.5 s window exceeded 4 cm/s. Immobility was defined as periods of not running that exceeded a duration of 0.5 s.

Calcium imaging

CaImAn software (version 1.8.5, Python) was used to motion-correct any movements between the calcium imaging frames, identify the cells and extract their fluorescence signal from the video recordings (Giovannucci et al., 2019). The method for the cell and signal detection is based on constrained non-negative matrix factorization (CNMF-E; Pnevmatikakis et al., 2016)⁠. CaImAn extracted background-subtracted calcium fluorescence values and the deconvolved the signal. The deconvolved signal can be interpreted as a scaled probability of a neuron being active. The calcium imaging videos recorded in the same-day trials were concatenated and motion-corrected to a common template frame. Signal extraction and further processing were performed on the resulting long video, allowing to detect the cells and signal present across the trials. To improve the computational performance, the video frames were cropped to a rectangle containing the regions of interest, and the video width and height were downsampled by a factor of 2.

The identified putative cells were automatically filtered using CaImAn. The results were visually inspected and the filtering parameters adjusted to exclude non-cell like shapes and traces from the filtered components. The criteria used for the filtering included a threshold for signal to noise ratio of the trace, the minimum and maximum size of the component‘s spatial footprint, threshold for consistency of the spatial footprint at different times of the component‘s activation, and a threshold for component‘s resemblance to a neuronal soma as evaluated by a convolutional neural network provided with CaImAn software.

The deconvolved trace was time binned, averaging the values in 200 ms bins. A calcium event was detected whenever the cell‘s deconvolved signal crossed 20% of its day-maximum value. A cell was classified as active during locomotion if it had at least one calcium event. A cell was classified as reactivated if it had at least one calcium event during immobility period and it was active during preceding locomotion. If during the immobility period mice were located within 6 cm of the reward, the reactivation was classified as reactivation at reward, otherwise as reactivation at non-rewarded location. Activity peaks were quantified by their area-under-curve (AUC) which was calculated as a convolution of the preprocessed calcium signal with a 2-s-long flat kernel. The preprocessing of calcium signal subtracted the cell’s median value and truncated the values below 0, so that only the above-median calcium signal is integrated in the AUC calculation. We excluded any samples from cells whose maximum value AUC in a given trial did not exceed 0.

Place cell detection and analysis

To assess how spatial locations modulated activity of a cell, we considered periods of running as described in the 'Behavioral learning task' section and calculated place maps — mean neural activity per spatial bin. The total activity inside 6×6 cm bins was summed from the smoothed deconvolved signal. The mean neural activity in the spatial bin was then calculated as a ratio of the total activity to the total occupancy in the bin after both maps were smoothed across the space using a 2D Gaussian kernel with σ=12 cm. The place map was filtered to include spatial bins with total occupancy that exceeded 1 s (5 time bins, thresholded on unsmoothed total occupancy).

Spatial information of a cell‘s activity was calculated using the place map values. Spatial information (Markus et al., 1994) was defined as:

S p a t i a l I n f o r m a t i o n = \sum_{i = 1}^{N} p_{i} \frac{λ_{i}}{\bar{λ}} l o g_{2} (\frac{λ_{i}}{\bar{λ}})

where $\bar{λ}$ represents the mean value of the neural signal, $p_{i}$ represents probability of the occupancy of the i-th bin, and $λ_{i}$ represents its mean neural activity. Dividing by $\bar{λ}$ ensures the metric is independent of the cell‘s average activity. The units of spatial information calculated on calcium fluorescence can be reported as bits per action potential (Climer and Dombeck, 2021). However, because the actual action potentials were not measured, we report them as arbitrary units.

Spatial information was compared to the value expected by chance. The chance level was calculated by circularly shifting the activity with regards to the actual location. For each cell, the activity was circularly shifted within the trial by a time offset chosen randomly (minimum offset of 10 s). If the cell‘s spatial information exceeded 95% values calculated on 1000 random shifts of its activity, it was defined as a place cell.

A limited number of neuronal responses sampled per spatial bin can lead to an upward bias in the estimated spatial information (Treves and Panzeri, 1995). To correct this bias, we subtracted its estimated value from the estimated spatial information. The bias was estimated as the mean spatial information from the time-shifting procedure used for place cell detection. This procedure does not require binning of the neuronal responses from the calcium imaging as required by analytical estimation (Panzeri et al., 2007), and has been used previously to estimate mutual information bias (Akrami et al., 2018).

Statistical testing of in vivo results

To compare the activity peaks between the reactivated and non-reactivated cells throughout the day, we used permutation tests for repeated measures ANOVA (Figure 8E, Figure 8—figure supplement 1A,B, Figure 8—figure supplement 2). The ANOVA modeled fixed effects of trial ordinal and reactivation and random effects of within animal-session factors. Significance level was set to α=0.05. To compare the effect of reactivation in each trial on cells grouped by the trial of their first activation, we used permutation t-tests (Figure 8F, Figure 8—figure supplement 1C). Multiple comparisons were corrected with Benjamini-Hochberg method with the type I error rate set to 0.05. The permutations were restricted to within animal-session permutations. Both ANOVA and t-test statistics were computed based on and 10,000 permutations using ‘permuco‘ R package. Statistical analysis was performed in R version 3.6.3.

Data availability

Data availability Experimental data and code are available at: Code for computational model and code for in vivo analysis (including a link to in vivo data) are available at: https://github.com/przemyslawj/dCA1-reactivations copy archived at swh:1:rev:22a4e82293f6c36c6fef8c0f06c3f6c68c4045ad. Data of plasticity experiments and of simulation data from computational model are available at: https://data.mendeley.com/datasets/dx7cdgpcz3/1.

The following data sets were generated

1. Jarzebowski P
(2022) Github
ID github.com/przemyslawj/dCA1-reactivations. dCA1-reactivations.

https://github.com/przemyslawj/dCA1-reactivations
1. Fuchsberger T
(2022) Mendeley Data
Postsynaptic burst reactivation of hippocampal neurons enables associative plasticity of temporally discontiguous inputs.

https://doi.org/10.17632/dx7cdgpcz3.1

References

1. Akrami A
2. Kopec CD
3. Diamond ME
4. Brody CD
(2018) Posterior parietal cortex represents sensory history and mediates its effects on behaviour
Nature 554:368–372.

https://doi.org/10.1038/nature25510
- PubMed
- Google Scholar
(2016) Reverse replay of hippocampal place cells is uniquely modulated by changing reward
Neuron 91:1124–1136.

https://doi.org/10.1016/j.neuron.2016.07.047
- PubMed
- Google Scholar
(2016) Presynaptic spike timing-dependent long-term depression in the mouse hippocampus
Cerebral Cortex 26:3637–3654.

https://doi.org/10.1093/cercor/bhw172
- PubMed
- Google Scholar
(2017) Behavioral time scale synaptic plasticity underlies CA1 place fields
Science 357:1033–1036.

https://doi.org/10.1126/science.aan3846
- PubMed
- Google Scholar
(2015) Retroactive modulation of spike timing-dependent plasticity by dopamine
eLife 4:e09685.

https://doi.org/10.7554/eLife.09685
- PubMed
- Google Scholar
1. Buzsáki G
(1989) Two-stage model of memory trace formation: a role for “noisy” brain states
Neuroscience 31:551–570.

https://doi.org/10.1016/0306-4522(89)90423-5
- PubMed
- Google Scholar
1. Buzsáki G
2. Horváth Z
3. Urioste R
4. Hetke J
5. Wise K
(1992) High-frequency network oscillation in the hippocampus
Science 256:1025–1027.

https://doi.org/10.1126/science.1589772
- PubMed
- Google Scholar
1. Cacucci F
2. Wills TJ
3. Lever C
4. Giese KP
5. O’Keefe J
(2007) Experience-dependent increase in CA1 place cell spatial information, but not spatial reproducibility, is dependent on the autophosphorylation of the alpha-isoform of the calcium/calmodulin-dependent protein kinase II
The Journal of Neuroscience 27:7854–7859.

https://doi.org/10.1523/JNEUROSCI.1704-07.2007
- PubMed
- Google Scholar
1. Cai DJ
2. Aharoni D
3. Shuman T
4. Shobe J
5. Biane J
6. Song W
7. Wei B
8. Veshkini M
9. La-Vu M
10. Lou J
11. Flores SE
12. Kim I
13. Sano Y
14. Zhou M
15. Baumgaertel K
16. Lavi A
17. Kamata M
18. Tuszynski M
19. Mayford M
20. Golshani P
21. Silva AJ
(2016) A shared neural ensemble links distinct contextual memories encoded close in time
Nature 534:115–118.

https://doi.org/10.1038/nature17955
- PubMed
- Google Scholar
(2021) A locus coeruleus- dorsal CA1 dopaminergic circuit modulates memory linking
Neuron 110:3374–3388.

https://doi.org/10.2139/ssrn.3985190
- Google Scholar
1. Climer JR
2. Dombeck DA
(2021) Information theoretic approaches to deciphering the neural code with functional fluorescence imaging
ENeuro 8:ENEURO.0266-21.2021.

https://doi.org/10.1523/ENEURO.0266-21.2021
- PubMed
- Google Scholar
(2007) Place-selective firing contributes to the reverse-order reactivation of CA1 pyramidal cells during sharp waves in open-field exploration
The European Journal of Neuroscience 26:704–716.

https://doi.org/10.1111/j.1460-9568.2007.05684.x
- PubMed
- Google Scholar
1. Dana H
2. Chen TW
3. Hu A
4. Shields BC
5. Guo C
6. Looger LL
7. Kim DS
8. Svoboda K
(2014) Thy1-GCaMP6 transgenic mice for neuronal population imaging in vivo
PLOS ONE 9:e108697.

https://doi.org/10.1371/journal.pone.0108697
- PubMed
- Google Scholar
1. Dayan P
(1993) Improving generalization for temporal difference learning: the successor representation
Neural Computation 5:613–624.

https://doi.org/10.1162/neco.1993.5.4.613
- Google Scholar
1. Diba K
2. Buzsáki G
(2007) Forward and reverse hippocampal place-cell sequences during ripples
Nature Neuroscience 10:1241–1242.

https://doi.org/10.1038/nn1961
- PubMed
- Google Scholar
1. Dudek SM
2. Fields RD
(2002) Somatic action potentials are sufficient for late-phase LTP-related cell signaling
PNAS 99:3962–3967.

https://doi.org/10.1073/pnas.062510599
- PubMed
- Google Scholar
(2010) The reorganization and reactivation of hippocampal maps predict spatial memory performance
Nature Neuroscience 13:995–1002.

https://doi.org/10.1038/nn.2599
- PubMed
- Google Scholar
1. Elizarova S
2. Chouaib AA
3. Shaib A
4. Hill B
5. Mann F
6. Brose N
7. Kruss S
8. Daniel JA
(2022) A fluorescent nanosensor paint detects dopamine release at axonal varicosities with high spatiotemporal resolution
PNAS 119:e2202842119.

https://doi.org/10.1073/pnas.2202842119
- PubMed
- Google Scholar
1. Ferguson GD
2. Storm DR
(2004) Why calcium-stimulated adenylyl cyclases?
Physiology 19:271–276.

https://doi.org/10.1152/physiol.00010.2004
- PubMed
- Google Scholar
1. Forsberg M
2. Seth H
3. Björefeldt A
4. Lyckenvik T
5. Andersson M
6. Wasling P
7. Zetterberg H
8. Hanse E
(2019) Ionized calcium in human cerebrospinal fluid and its influence on intrinsic and synaptic excitability of hippocampal pyramidal neurons in the rat
Journal of Neurochemistry 149:452–470.

https://doi.org/10.1111/jnc.14693
- PubMed
- Google Scholar
1. Foster DJ
2. Wilson MA
(2006) Reverse replay of behavioural sequences in hippocampal place cells during the awake state
Nature 440:680–683.

https://doi.org/10.1038/nature04587
- PubMed
- Google Scholar
1. Frémaux N
2. Gerstner W
(2016) Neuromodulated spike-timing-dependent plasticity, and theory of three-factor learning rules
Frontiers in Neural Circuits 9:85.

https://doi.org/10.3389/fncir.2015.00085
- PubMed
- Google Scholar
1. Frey U
2. Krug M
3. Reymann KG
4. Matthies H
(1988) Anisomycin, an inhibitor of protein synthesis, blocks late phases of LTP phenomena in the hippocampal CA1 region in vitro
Brain Research 452:57–65.

https://doi.org/10.1016/0006-8993(88)90008-x
- PubMed
- Google Scholar
(1990) Dopaminergic antagonists prevent long-term maintenance of posttetanic LTP in the CA1 region of rat hippocampal slices
Brain Research 522:69–75.

https://doi.org/10.1016/0006-8993(90)91578-5
- PubMed
- Google Scholar
(1993) Effects of cAMP simulate a late stage of LTP in hippocampal CA1 neurons
Science 260:1661–1664.

https://doi.org/10.1126/science.8389057
- PubMed
- Google Scholar
(1994) Mesolimbic dopaminergic neurons innervating the hippocampal formation in the rat: a combined retrograde tracing and immunohistochemical study
Brain Research 668:71–79.

https://doi.org/10.1016/0006-8993(94)90512-6
- PubMed
- Google Scholar
1. Gerstner W
2. Lehmann M
3. Liakoni V
4. Corneil D
5. Brea J
(2018) Eligibility traces and plasticity on behavioral time scales: experimental support of neoHebbian three-factor learning rules
Frontiers in Neural Circuits 12:53.

https://doi.org/10.3389/fncir.2018.00053
- PubMed
- Google Scholar
1. Ghosh KK
2. Burns LD
3. Cocker ED
4. Nimmerjahn A
5. Ziv Y
6. Gamal AE
7. Schnitzer MJ
(2011) Miniaturized integration of a fluorescence microscope
Nature Methods 8:871–878.

https://doi.org/10.1038/nmeth.1694
- PubMed
- Google Scholar
1. Giovannucci A
2. Friedrich J
3. Gunn P
4. Kalfon J
5. Brown BL
6. Koay SA
7. Taxidis J
8. Najafi F
9. Gauthier JL
10. Zhou P
11. Khakh BS
12. Tank DW
13. Chklovskii DB
14. Pnevmatikakis EA
(2019) CaImAn an open source tool for scalable calcium imaging data analysis
eLife 8:e38173.

https://doi.org/10.7554/eLife.38173
- PubMed
- Google Scholar
Preprint
1. Grienberger C
2. Magee JC
(2021) Entorhinal Cortex Directs Learning-Related Changes in CA1 Representations
bioRxiv.

https://doi.org/10.1101/2021.12.10.472158
- Google Scholar
(2021) Reactivation predicts the consolidation of unbiased long-term cognitive maps
Nature Neuroscience 24:1574–1585.

https://doi.org/10.1038/s41593-021-00920-7
- PubMed
- Google Scholar
1. Grover LM
2. Teyler TJ
(1990) Two components of long-term potentiation induced by different patterns of afferent activation
Nature 347:477–479.

https://doi.org/10.1038/347477a0
- PubMed
- Google Scholar
1. Huang YY
2. Kandel ER
(1995) D1/D5 receptor agonists induce a protein synthesis-dependent late potentiation in the CA1 region of the hippocampus
PNAS 92:2446–2450.

https://doi.org/10.1073/pnas.92.7.2446
- PubMed
- Google Scholar
(2021) Prioritized experience replays on a hippocampal predictive map for learning
PNAS 118:e2011266118.

https://doi.org/10.1073/pnas.2011266118
- PubMed
- Google Scholar
(2020) Synaptic plasticity rules with physiological calcium levels
PNAS 117:33639–33648.

https://doi.org/10.1073/pnas.2013663117
- PubMed
- Google Scholar
1. Iordanov MS
2. Pribnow D
3. Magun JL
4. Dinh TH
5. Pearson JA
6. Chen SL
7. Magun BE
(1997) Ribotoxic stress response: activation of the stress-activated protein kinase JNK1 by inhibitors of the peptidyl transferase reaction and by sequence-specific RNA damage to the alpha-sarcin/ricin loop in the 28S rRNA
Molecular and Cellular Biology 17:3373–3381.

https://doi.org/10.1128/MCB.17.6.3373
- PubMed
- Google Scholar
1. Jones HC
2. Keep RF
(1988) Brain fluid calcium concentration and response to acute hypercalcaemia during development in the rat
The Journal of Physiology 402:579–593.

https://doi.org/10.1113/jphysiol.1988.sp017223
- PubMed
- Google Scholar
(2020) A role for the locus coeruleus in hippocampal CA1 place cell reorganization during spatial reward learning
Neuron 105:1018–1026.

https://doi.org/10.1016/j.neuron.2019.12.029
- PubMed
- Google Scholar
(2016) Dopamine release from the locus coeruleus to the dorsal hippocampus promotes spatial learning and memory
PNAS 113:14835–14840.

https://doi.org/10.1073/pnas.1616515114
- PubMed
- Google Scholar
(1999) Reactivation of hippocampal cell assemblies: effects of behavioral state, experience, and EEG dynamics
The Journal of Neuroscience 19:4090–4101.

https://doi.org/10.1523/JNEUROSCI.19-10-04090.1999
- PubMed
- Google Scholar
1. Lee AK
2. Wilson MA
(2002) Memory of sequential experience in the hippocampus during slow wave sleep
Neuron 36:1183–1194.

https://doi.org/10.1016/s0896-6273(02)01096-6
- PubMed
- Google Scholar
1. Lisman JE
(1997) Bursts as a unit of neural information: making unreliable synapses reliable
Trends in Neurosciences 20:38–43.

https://doi.org/10.1016/S0166-2236(96)10070-9
- PubMed
- Google Scholar
1. Liu C
2. Goel P
3. Kaeser PS
(2021) Spatial and temporal scales of dopamine transmission
Nature Reviews. Neuroscience 22:345–358.

https://doi.org/10.1038/s41583-021-00455-7
- PubMed
- Google Scholar
(1992) Responses of monkey dopamine neurons during learning of behavioral reactions
Journal of Neurophysiology 67:145–163.

https://doi.org/10.1152/jn.1992.67.1.145
- PubMed
- Google Scholar
1. Lopes JP
2. Cunha RA
(2019) What is the extracellular calcium concentration within brain synapses?
Journal of Neurochemistry 149:435–437.

https://doi.org/10.1111/jnc.14696
- PubMed
- Google Scholar
(1994) Spatial information content and reliability of hippocampal CA1 neurons: effects of visual input
Hippocampus 4:410–421.

https://doi.org/10.1002/hipo.450040404
- PubMed
- Google Scholar
1. Mathis A
2. Mamidanna P
3. Cury KM
4. Abe T
5. Murthy VN
6. Mathis MW
7. Bethge M
(2018) DeepLabCut: markerless pose estimation of user-defined body parts with deep learning
Nature Neuroscience 21:1281–1289.

https://doi.org/10.1038/s41593-018-0209-y
- PubMed
- Google Scholar
1. Matthies H
2. Becker A
3. Schröeder H
4. Kraus J
5. Höllt V
6. Krug M
(1997) Dopamine D1-deficient mutant mice do not express the late phase of hippocampal long-term potentiation
Neuroreport 8:3533–3535.

https://doi.org/10.1097/00001756-199711100-00023
- PubMed
- Google Scholar
(2014) Dopaminergic neurons promote hippocampal reactivation and spatial memory persistence
Nature Neuroscience 17:1658–1660.

https://doi.org/10.1038/nn.3843
- PubMed
- Google Scholar
(1997) Experience-dependent, asymmetric expansion of hippocampal place fields
PNAS 94:8918–8921.

https://doi.org/10.1073/pnas.94.16.8918
- PubMed
- Google Scholar
(2000) Experience-dependent asymmetric shape of hippocampal receptive fields
Neuron 25:707–715.

https://doi.org/10.1016/s0896-6273(00)81072-7
- PubMed
- Google Scholar
1. Mehta MR
(2015) From synaptic plasticity to spatial maps and sequence learning
Hippocampus 25:756–762.

https://doi.org/10.1002/hipo.22472
- PubMed
- Google Scholar
1. Milstein AD
2. Li Y
3. Bittner KC
4. Grienberger C
5. Soltesz I
6. Magee JC
7. Romani S
(2021) Bidirectional synaptic plasticity rapidly modifies hippocampal representations
eLife 10:e73046.

https://doi.org/10.7554/eLife.73046
- PubMed
- Google Scholar
1. Moore JJ
2. Cushman JD
3. Acharya L
4. Popeney B
5. Mehta MR
(2021) Linking hippocampal multiplexed tuning, Hebbian plasticity and navigation
Nature 599:442–448.

https://doi.org/10.1038/s41586-021-03989-z
- PubMed
- Google Scholar
(1999) Replay and time compression of recurring spike sequences in the hippocampus
The Journal of Neuroscience 19:9497–9507.

https://doi.org/10.1523/JNEUROSCI.19-21-09497.1999
- PubMed
- Google Scholar
(2004) Dopamine receptor signaling
Journal of Receptor and Signal Transduction Research 24:165–205.

https://doi.org/10.1081/rrs-200029981
- PubMed
- Google Scholar
(2006) Place-selective firing of CA1 pyramidal cells during sharp wave/ripple network patterns in exploratory behavior
Neuron 49:143–155.

https://doi.org/10.1016/j.neuron.2005.10.037
- PubMed
- Google Scholar
(2007) Correcting for the sampling bias problem in spike train information measures
Journal of Neurophysiology 98:1064–1072.

https://doi.org/10.1152/jn.00559.2007
- PubMed
- Google Scholar
1. Park P
2. Volianskis A
3. Sanderson TM
4. Bortolotto ZA
5. Jane DE
6. Zhuo M
7. Kaang BK
8. Collingridge GL
(2014) NMDA receptor-dependent long-term potentiation comprises a family of temporally overlapping forms of synaptic plasticity that are induced by different patterns of stimulation
Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 369:20130131.

https://doi.org/10.1098/rstb.2013.0131
- PubMed
- Google Scholar
1. Payeur A
2. Guerguiev J
3. Zenke F
4. Richards BA
5. Naud R
(2021) Burst-dependent synaptic plasticity can coordinate learning in hierarchical circuits
Nature Neuroscience 24:1010–1019.

https://doi.org/10.1038/s41593-021-00857-x
- PubMed
- Google Scholar
(1999) Postsynaptic bursting is essential for “ hebbian” induction of associative long-term potentiation at excitatory synapses in rat hippocampus
The Journal of Physiology 518 (Pt 2):571–576.

https://doi.org/10.1111/j.1469-7793.1999.0571p.x
- PubMed
- Google Scholar
1. Pnevmatikakis EA
2. Soudry D
3. Gao Y
4. Machado TA
5. Merel J
6. Pfau D
7. Reardon T
8. Mu Y
9. Lacefield C
10. Yang W
11. Ahrens M
12. Bruno R
13. Jessell TM
14. Peterka DS
15. Yuste R
16. Paninski L
(2016) Simultaneous denoising, deconvolution, and demixing of calcium imaging data
Neuron 89:285–299.

https://doi.org/10.1016/j.neuron.2015.11.037
- PubMed
- Google Scholar
1. Redondo RL
2. Morris RG
(2011) Making memories last: the synaptic tagging and capture hypothesis
Nat Rev Neurosci 12:17–30.

https://doi.org/10.1038/nrn2963
- Google Scholar
1. Resendez SL
2. Jennings JH
3. Ung RL
4. Namboodiri VMK
5. Zhou ZC
6. Otis JM
7. Nomura H
8. McHenry JA
9. Kosyk O
10. Stuber GD
(2016) Visualization of cortical, subcortical and deep brain neural circuit dynamics during naturalistic mammalian behavior with head-mounted microscopes and chronically implanted lenses
Nature Protocols 11:566–597.

https://doi.org/10.1038/nprot.2016.021
- PubMed
- Google Scholar
1. Sassone-Corsi P
(2012) The cyclic AMP pathway
Cold Spring Harbor Perspectives in Biology 4:a011148.

https://doi.org/10.1101/cshperspect.a011148
- PubMed
- Google Scholar
(1980) Origin of dopaminergic innervation of the rat hippocampal formation
Neuroscience Letters 18:125–131.

https://doi.org/10.1016/0304-3940(80)90314-6
- PubMed
- Google Scholar
(1997) A neural substrate of prediction and reward
Science 275:1593–1599.

https://doi.org/10.1126/science.275.5306.1593
- PubMed
- Google Scholar
1. Schultz W
(2007) Behavioral dopamine signals
Trends in Neurosciences 30:203–210.

https://doi.org/10.1016/j.tins.2007.03.007
- PubMed
- Google Scholar
1. Sehgal M
2. Zhou M
3. Lavi A
4. Huang S
5. Zhou Y
6. Silva AJ
(2018) Memory allocation mechanisms underlie memory linking across time
Neurobiology of Learning and Memory 153:21–25.

https://doi.org/10.1016/j.nlm.2018.02.021
- PubMed
- Google Scholar
1. Shipton OA
2. Paulsen O
(2014) Glun2A and GluN2B subunit-containing NMDA receptors in hippocampal plasticity
Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 369:20130163.

https://doi.org/10.1098/rstb.2013.0163
- PubMed
- Google Scholar
1. Silva AJ
2. Zhou Y
3. Rogerson T
4. Shobe J
5. Balaji J
(2009) Molecular and cellular approaches to memory allocation in neural circuits
Science 326:391–395.

https://doi.org/10.1126/science.1174519
- PubMed
- Google Scholar
1. Singer AC
2. Frank LM
(2009) Rewarded outcomes enhance reactivation of experience in the hippocampus
Neuron 64:910–921.

https://doi.org/10.1016/j.neuron.2009.11.016
- PubMed
- Google Scholar
1. Smith CC
2. Greene RW
(2012) CNS dopamine transmission mediated by noradrenergic innervation
The Journal of Neuroscience 32:6072–6080.

https://doi.org/10.1523/JNEUROSCI.6486-11.2012
- PubMed
- Google Scholar
(2017) The hippocampus as a predictive map
Nature Neuroscience 20:1643–1653.

https://doi.org/10.1038/nn.4650
- PubMed
- Google Scholar
1. Sun F
2. Zeng J
3. Jing M
4. Zhou J
5. Feng J
6. Owen SF
7. Luo Y
8. Li F
9. Wang H
10. Yamaguchi T
11. Yong Z
12. Gao Y
13. Peng W
14. Wang L
15. Zhang S
16. Du J
17. Lin D
18. Xu M
19. Kreitzer AC
20. Cui G
21. Li Y
(2018) A genetically encoded fluorescent sensor enables rapid and specific detection of dopamine in flies, fish, and mice
Cell 174:481–496.

https://doi.org/10.1016/j.cell.2018.06.042
- PubMed
- Google Scholar
1. Sun F
2. Zhou J
3. Dai B
4. Qian T
5. Zeng J
6. Li X
7. Zhuo Y
8. Zhang Y
9. Wang Y
10. Qian C
11. Tan K
12. Feng J
13. Dong H
14. Lin D
15. Cui G
16. Li Y
(2020) Next-Generation GRAB sensors for monitoring dopaminergic activity in vivo
Nature Methods 17:1156–1166.

https://doi.org/10.1038/s41592-020-00981-9
- PubMed
- Google Scholar
Book
1. Sutton RS
2. Barto AG
(2018)
Reinforcement Learning: An Introduction

Cambridge, MA: The MIT Press.
- Google Scholar
(2016) Locus coeruleus and dopaminergic consolidation of everyday memory
Nature 537:357–362.

https://doi.org/10.1038/nature19325
- PubMed
- Google Scholar
1. Treves A
2. Panzeri S
(1995) The upward bias in measures of information derived from limited data samples
Neural Computation 7:399–407.

https://doi.org/10.1162/neco.1995.7.2.399
- Google Scholar
1. Wagatsuma A
2. Okuyama T
3. Sun C
4. Smith LM
5. Abe K
6. Tonegawa S
(2018) Locus coeruleus input to hippocampal CA3 drives single-trial learning of a novel context
PNAS 115:E310–E316.

https://doi.org/10.1073/pnas.1714082115
- PubMed
- Google Scholar
1. Watson EL
2. Jacobson KL
3. Singh JC
4. Idzerda R
5. Ott SM
6. DiJulio DH
7. Wong ST
8. Storm DR
(2000) The type 8 adenylyl cyclase is critical for Ca²⁺ stimulation of camp accumulation in mouse parotid acini
The Journal of Biological Chemistry 275:14691–14699.

https://doi.org/10.1074/jbc.275.19.14691
- PubMed
- Google Scholar
1. Wayman GA
2. Impey S
3. Wu Z
4. Kindsvogel W
5. Prichard L
6. Storm DR
(1994)
Synergistic activation of the type I adenylyl cyclase by Ca²⁺ and Gs-coupled receptors in vivo

The Journal of Biological Chemistry 269:25400–25405.
- PubMed
- Google Scholar
1. Wilson MA
2. McNaughton BL
(1994) Reactivation of hippocampal ensemble memories during sleep
Science 265:676–679.

https://doi.org/10.1126/science.8036517
- PubMed
- Google Scholar
1. Wong ST
2. Athos J
3. Figueroa XA
4. Pineda VV
5. Schaefer ML
6. Chavkin CC
7. Muglia LJ
8. Storm DR
(1999) Calcium-stimulated adenylyl cyclase activity is critical for hippocampus-dependent long-term memory and late phase LTP
Neuron 23:787–798.

https://doi.org/10.1016/s0896-6273(01)80036-2
- PubMed
- Google Scholar
1. Yiu AP
2. Mercaldo V
3. Yan C
4. Richards B
5. Rashid AJ
6. Hsiang H-LL
7. Pressey J
8. Mahadevan V
9. Tran MM
10. Kushner SA
11. Woodin MA
12. Frankland PW
13. Josselyn SA
(2014) Neurons are recruited to a memory trace based on relative neuronal excitability immediately before training
Neuron 83:722–735.

https://doi.org/10.1016/j.neuron.2014.07.017
- PubMed
- Google Scholar
1. Ziv Y
2. Burns LD
3. Cocker ED
4. Hamel EO
5. Ghosh KK
6. Kitch LJ
7. El Gamal A
8. Schnitzer MJ
(2013) Long-term dynamics of CA1 hippocampal place codes
Nature Neuroscience 16:264–266.

https://doi.org/10.1038/nn.3329
- PubMed
- Google Scholar

Article and author information

Author details

Tanja Fuchsberger

Department of Physiology, Development and Neuroscience, Physiological Laboratory, University of Cambridge, Cambridge, United Kingdom

Contribution
Conceptualization, Investigation, Writing – original draft, Writing – review and editing, Designed the experiments, Conducted the experiments and analyzed the data, Wrote the manuscript

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-4751-8806
Claudia Clopath

Department of Bioengineering, Imperial College London, London, United Kingdom

Contribution
Software, Investigation, Writing – original draft, Writing – review and editing, Developed the computational model, Wrote the manuscript

Contributed equally with
Przemyslaw Jarzebowski

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-4507-8648
Przemyslaw Jarzebowski

Department of Physiology, Development and Neuroscience, Physiological Laboratory, University of Cambridge, Cambridge, United Kingdom

Contribution
Software, Formal analysis, Investigation, Writing – original draft, Writing – review and editing, Conducted the experiments and analyzed the data, Wrote the manuscript

Contributed equally with
Claudia Clopath

Competing interests
No competing interests declared
Zuzanna Brzosko

Department of Physiology, Development and Neuroscience, Physiological Laboratory, University of Cambridge, Cambridge, United Kingdom

Contribution
Investigation, Writing – review and editing, Conducted the experiments and analyzed the data

Competing interests
No competing interests declared
Hongbing Wang

Department of Physiology, Michigan State University, East Lansing, United States

Contribution
Resources, Provided transgenic AC DKO mice

Competing interests
No competing interests declared
Ole Paulsen

Department of Physiology, Development and Neuroscience, Physiological Laboratory, University of Cambridge, Cambridge, United Kingdom

Contribution
Conceptualization, Supervision, Funding acquisition, Writing – original draft, Writing – review and editing, Designed the experiments, Wrote the manuscript

For correspondence
op210@cam.ac.uk

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-2258-5455

Funding

Biotechnology and Biological Sciences Research Council (BB/N019008/1)

Tanja Fuchsberger

Biotechnology and Biological Sciences Research Council (BB/P019560/1)

Tanja Fuchsberger

Biotechnology and Biological Sciences Research Council (Studentship)

Przemyslaw Jarzebowski

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

Acknowledgements

This research was supported by the Biotechnology and Biological Sciences Research Council, U.K. We are grateful for discussions of this project with other members of the Neuronal Oscillations Group.

Ethics

Copyright

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.