Local generation of multineuronal spike sequences in the hippocampal CA1 region
Edited by Howard Eichenbaum, Center for Memory and Brain, Boston University, Boston, MA, and accepted by the Editorial Board July 12, 2015 (received for review May 9, 2015)
Significance
Neuronal sequences in the CA1 hippocampal region during short burst events (called sharp wave ripples) reliably reflect the slower sequential order of place activation during exploration, patterned by local-circuit interneurons. We tested the hypothesis that in addition to afferent drive, local mechanisms contribute critically to the sequential patterning of pyramidal cells and the diverse interneuron types. We show that optogenetic activation of a small group of CA1 pyramidal cells in mouse hippocampus yields consistent multineuronal spike sequences, and that their sequential firing order is similar to that observed during spontaneous “ripple” and theta oscillations. Thus, multineuronal spike sequences are not exclusively inherited from upstream sources but result from local-circuit dynamics and biophysical properties of the individual neurons.
Abstract
Sequential activity of multineuronal spiking can be observed during theta and high-frequency ripple oscillations in the hippocampal CA1 region and is linked to experience, but the mechanisms underlying such sequences are unknown. We compared multineuronal spiking during theta oscillations, spontaneous ripples, and focal optically induced high-frequency oscillations (“synthetic” ripples) in freely moving mice. Firing rates and rate modulations of individual neurons, and multineuronal sequences of pyramidal cell and interneuron spiking, were correlated during theta oscillations, spontaneous ripples, and synthetic ripples. Interneuron spiking was crucial for sequence consistency. These results suggest that participation of single neurons and their sequential order in population events are not strictly determined by extrinsic inputs but also influenced by local-circuit properties, including synapses between local neurons and single-neuron biophysics.
Sign up for PNAS alerts.
Get alerts for new articles, or get an alert when an article is cited.
A hypothesized hallmark of cognition is self-organized sequential activation of neuronal assemblies (1). Self-organized neuronal sequences have been observed in several cortical structures (2–5) and neuronal models (6–7). In the hippocampus, sequential activity of place cells (8) may be induced by external landmarks perceived by the animal during spatial navigation (9) and conveyed to CA1 by the upstream CA3 region or layer 3 of the entorhinal cortex (10). Internally generated sequences have been also described in CA1 during theta oscillations in memory tasks (4, 11), raising the possibility that a given neuronal substrate is responsible for generating sequences at multiple time scales. The extensive recurrent excitatory collateral system of the CA3 region has been postulated to be critical in this process (4, 7, 12, 13).
The sequential activity of place cells is “replayed” during sharp waves (SPW) in a temporally compressed form compared with rate modulation of place cells (14–20) and may arise from the CA3 recurrent excitatory networks during immobility and slow wave sleep. The SPW-related convergent depolarization of CA1 neurons gives rise to a local, fast oscillatory event in the CA1 region (“ripple,” 140–180 Hz; refs. 8 and 21). Selective elimination of ripples during or after learning impairs memory performance (22–24), suggesting that SPW ripple-related replay assists memory consolidation (12, 13). Although the local origin of the ripple oscillations is well demonstrated (25, 26), it has been tacitly assumed that the ripple-associated, sequentially ordered firing of CA1 neurons is synaptically driven by the upstream CA3 cell assemblies (12, 15), largely because excitatory recurrent collaterals in the CA1 region are sparse (27). However, sequential activity may also emerge by local mechanisms, patterned by the different biophysical properties of CA1 pyramidal cells and their interactions with local interneurons, which discharge at different times during a ripple (28–30). A putative function of the rich variety of interneurons is temporal organization of principal cell spiking (29–32). We tested the “local-circuit” hypothesis by comparing the probability of participation and sequential firing of CA1 neurons during theta oscillations, natural spontaneous ripple events, and “synthetic” ripples induced by local optogenetic activation of pyramidal neurons.
Results
Spike Rates and Multineuronal Spike Sequences During Theta Waves and Ripple Events Are Correlated.
The firing rates of CA1 pyramidal cells (PYR) and putative interneurons (INT) during theta oscillations and ripple events were correlated (rank correlations: PYR: 0.59, P < 0.001, permutation test; INT: 0.89, P < 0.001; n = 268 PYR and 37 INT recorded from four freely moving mice; Fig. 1A). To compare the rate modulation during theta and ripples on a neuron-by-neuron basis, we computed a “gain” measure for each unit, defined as the mean spiking rate during an event (ripple or theta cycle), divided by the baseline spiking rate. Although spike gain was higher during ripples than during theta (median ratio: PYR: 2.8, P < 0.001, Wilcoxon’s signed-rank paired test; INT: 2.2, P < 0.001), the spike gain of PYR during theta and ripples was correlated (rank correlations: PYR: 0.21, P < 0.001; INT: 0.27, P = 0.06; Fig. 1B). Thus, spike rates and their modulation are correlated during ripple and theta brain states.
Fig. 1.
During ripples, the spiking of multiple simultaneously recorded neurons (median, 24; range, 9–60) occurred in recurring sequences (median self-consistency: 0.28, P < 0.001, Wilcoxon’s signed-rank test, zero-median null, n = 11 sessions from four mice; Fig. 1 C–E) during most sessions (10 of 11 sessions, 91%, with self-consistent sequences, P < 0.05, permutation test). The same neurons exhibited sequential spiking also during individual theta cycles (median self-consistency: 0.24, P < 0.001) during all sessions (11 of 11, 100%; Fig. 1 C–E). To compare the relationship between sequential ordering of neurons during individual ripple events and theta cycles (“cross-class consistency”), we used the spike sequences during ripple events as a template against which spiking during each theta cycle was compared. In 10 of 11 cases (91%, chance level, 5%, P < 0.001, exact Binomial test), the sequences during ripples and theta were similar (median cross-class consistency: 0.21, P < 0.001; Fig. 1 C–E). Similar results were observed when only INT (median, 4; range, 2–12) were used in comparing sequences (median cross-class consistency: 0.18, P = 0.03; 8 of 9 sessions, P < 0.001; Fig. 1F). However, when only PYR (median, 22; range, 6–51) were considered, the median cross-class consistency fell to chance level (P = 0.83; Fig. 1G). On a session-by-session basis, 6 of 11 (55%; P < 0.001) sessions exhibited significant rank correlations, although half were consistent (3 of 11, 27%, P = 0.02; median cross-class consistency: 0.05) and half were reversed (3 of 11, 27%, P = 0.02; median, -0.04; Fig. 1G).
These correlations suggest that the same neuronal substrates are responsible for generating sequential firing of neurons during both theta and ripple states. Because layer 3 of the entorhinal input to CA1 is not active during ripples (33), it is unlikely that the entorhinal input is responsible for generating correlated sequences across different states. Sequences may be generated by the upstream CA3 region or locally in CA1 but these results alone cannot disambiguate the contributions of the two regions.
Single Neuron Spiking Modulations Are Correlated During Spontaneous and Induced Ripples.
Upon recording from the CA1 pyramidal cell layer of freely moving mice, high frequency “ripple” oscillations (21) were readily observed (Fig. 2Ab). During spontaneous ripple events, individual neurons increased their firing rate several-fold (28): pyramidal cells (PYR) about 5.1-fold (mean over n = 268 PYR; SEM: 0.5), and interneurons 4.1-fold (SEM: 0.6; n = 37 INT recorded from four freely moving CaMKII::ChR2 mice; Fig. 2B, Inset). Units with higher baseline (nonripple) firing rates also had higher firing rates during spontaneous ripples (rank correlation coefficient: PYR: 0.55, P < 0.001, permutation test; INT: 0.86, P < 0.001; Fig. 2B). In these same animals and brain state (Fig. S1), we also induced synthetic high-frequency oscillations (26) using focal application of blue light (estimated light intensity at the center of the CA1 pyramidal cell layer, 0.01–1 mW/mm2; Fig. 2Ac). Illumination in the opsin-free (control) mouse failed to induce oscillations or spikes (Fig. S2). During synthetic ripples induced at threshold intensity (the lowest intensity that generated ripples with comparable power to spontaneous ripples; ref. 26), firing rates of individual units also increased: PYR spiking increased 4.7-fold (SEM: 0.6), and INT spiking increased 1.8-fold (SEM: 0.4; Fig. 2C). The marked increase of PYR relative to INT spiking (P = 0.01, Mann–Whitney U test) is consistent with direct activation of PYR (and indirect activation of INT) in these animals during induced ripples. For both PYR and INT, the firing rates during induced ripples were correlated with the firing rates during spontaneous ripples on a neuron-by-neuron basis (rank correlations: PYR: 0.45, INT: 0.91, P < 0.001 for both; Fig. 2D). Spike gain was correlated during induced and spontaneous ripples (rank correlations: PYR: 0.24, INT: 0.74, P < 0.001 for both; Fig. 2E). Thus, firing rates and their modulations during induced ripples are not an arbitrary outcome of opsin- or optics-dependent parameters (expression level, distance from light source). Furthermore, the tendency of a given unit to modify its spike rate during CA3-driven (SPW induced) ripples is correlated with its tendency to modify spiking during locally driven (i.e., light-induced) ripples.
Fig. 2.
Fig. S1.
Fig. S2.
Correlated rate modulations are not sufficient to generate precise spike sequences. To determine the time relations between single neuron spiking and ripple cycles, we defined a phase for each spike that occurred during a ripple. Most units (224 of 268 PYR, 84%; 34 of 37 INT, 95%) exhibited consistent locking to ripple phase (P < 0.05, Rayleigh test), and phase-locked PYR spiked consistently earlier than phase-locked INT (mean± SEM phases: PYR: 156 ± 1.7°, INT: 254 ± 6.7°; corresponding to time lag of 1.8 ± 0.2 ms; Fig. 3A; Fig. 3C, x axis). Similar results were observed during induced ripples (mean ± SEM phases: PYR: 138 ± 4.4°, INT: 288 ± 9.8°; Fig. 3B; Fig. 3C, y axis). On a neuron-by-neuron basis, the spike phases during spontaneous and induced ripples were correlated (circular–circular correlation coefficient: PYR: 0.7, P < 0.001, χ2 test; INT: 0.42, P = 0.04; n = 97 PYR and 26 INT phase-locked during both spontaneous and induced ripples; Fig. 3). Therefore, the precise temporal relations between single unit spiking and ripple cycles are preserved during spontaneous and induced ripples.
Fig. 3.
Multineuronal Spike Sequences Are Correlated During Spontaneous Oscillations and Induced Ripples.
Because both spike rates and within-ripple cycle phase preference of spiking were correlated between spontaneous and induced ripples, we next examined whether sequential order of neuronal firing within ripple events is consistent. To this end, we rank ordered the spikes of all pyramidal cells and interneurons that occurred in a given event (Fig. 4A), and determined the consistency of multineuronal spiking across events (Fig. 4B); because induced ripples were focal (one site was illuminated at a time; Fig. 2Ac), this analysis was done on a site-by-site basis. During spontaneous ripples, most multineuronal sequences were self-consistent (23 of 27 recording sites, 85%; P < 0.001, exact Binomial test), with a median rank correlation of 0.28 (P < 0.001, Wilcoxon’s signed-rank test with a zero-median null; Fig. 4C). During induced ripples, 46% (12 of 26, P < 0.001) of the sequences were self-consistent, with a median rank correlation of 0.21 (P < 0.001; 0.32 for the consistent sequences; Fig. 4D). Similar results were observed under anesthesia (separate cohort of n = 5 urethane-anesthetized mice; Fig. S3). Thus, multineuronal spike sequences during induced ripples are not random.
Fig. 4.
Fig. S3.
Sequences during spontaneous and induced ripples were not independent: 13 of 26 induced sequences (50%, P < 0.001) were consistent with their same-site spontaneous counterparts, with a median rank correlation of 0.13 (P < 0.001; 0.21 for the consistent sequences; Fig. 4D). The median rank correlation was not different for the three datasets (spontaneous-ripple self-consistency; induced-ripple self-consistency; and cross-class consistency of induced-ripple, relative to spontaneous ripple sequences; P = 0.35, Kruskal–Wallis test). Similar results were observed in all freely moving animals (Table S1). Although the probability to observe consistent sequences was highest during spontaneous ripples (compared with induced-ripples and cross-class sequences, P = 0.003 and P = 0.006, χ2 test; Fig. 4E), consistent sequences occurred at similar probabilities during induced ripples and across classes (P = 0.78, Fig. 4E). Sequences during individual theta cycles and induced ripples were also interdependent: 10 of 25 (40%, P < 0.001) sites exhibited theta vs. induced ripples cross-class consistency, with a median rank correlation of 0.17 (P < 0.001; 0.32 for the consistent sequences; Fig. S4). Thus, sequences observed during spontaneous theta or ripple oscillations and induced ripples are similar, indicating that despite being driven by very different inputs (medial septum-associated theta, CA3-associated SPW, and local optically induced ChR2 activation), spiking during these three constellations share some common generative mechanisms.
Table S1.
Preparation | Mouse number | Number of sites | Consistent sequences during spontaneous ripples | Consistent sequences during induced ripples | Cross-class consistency |
---|---|---|---|---|---|
FM | 331 | 9 | 7 | 7 | 7 |
428 | 4 | 4 | 2 | 1 | |
482 | 9 | 9 | 2 | 4 | |
552 | 5 | 3 | 3 | 3 | |
Total | 27 | 23 | 14 | 15 | |
UA | 622 | 6 | 6 | 5 | 2 |
624 | 6 | 6 | 6 | 2 | |
629 | 6 | 6 | 5 | 1 | |
763 | 5 | 2 | 3 | 1 | |
765 | 6 | 6 | 4 | 0 | |
Total | 29 | 26 | 23 | 6 |
FM, freely moving; UA, urethane-anesthetized; consistent, P < 0.05, permutation test. All results pertain to template method, PYR+INT participation.
Fig. S4.
In constructing the multineuronal spike sequences, we considered spikes of all simultaneously recorded neurons, regardless of cell type. When limiting the participation to PYR, only 4 of 27 sites exhibited consistent sequences during spontaneous ripples (15%, P = 0.04, exact Binomial test) and the median rank correlation was lower (0.04, P = 0.02, signed-rank test). Similarly, only 7 of 22 (32%, P < 0.001; median rank, 0.13, P = 0.23) sequences were consistent during induced ripples, and cross-class consistency fell to chance level (3 of 22 sites, 13%; P = 0.09; median rank, 0.06, P = 0.25; Fig. 5A). Moreover, cross-class rank correlations were 3.4 times lower than when both PYR and INT were considered (P = 0.04, Wilcoxon’s paired signed-rank test; Fig. 5A, Right). When limiting participation to INT, sequences were still consistent (Fig. 5B), and cross-class consistency was 1.75 times lower than when all neurons were considered (P = 0.55, paired signed-rank test; rank correlation: 0.54, P < 0.001), suggesting that rate difference between PYR and INT and/or intrinsic PYR-INT correlations are not accountable for the observed consistency. Thus, although interneurons comprised only a small subpopulation of the participating units (12%) and were indirectly activated during optical ripple induction, they made an important contribution to the consistency of multineuronal sequences during both spontaneous and induced ripples.
Fig. 5.
We have used a cross-validated template method for determining the consistency of multineuronal sequences. Upon using an alternative “pair-wise” method, the exact values of the rank correlations were lower than those yielded by the template method (P = 0.002, signed-rank paired test; Fig. 5C). However, the numbers of consistent sequences (spontaneous ripples: 26 of 27, 96%; induced ripples: 14 of 25, 56%; cross-class consistency: 18 of 26, 69%) were similar to those detected using the template method (P > 0.05 for all three comparisons, χ2 test), and on a site-by-site basis, the values yielded by the two methods were correlated (0.58, P = 0.006, permutation test). Thus, although rank correlations were widely distributed and rarely exceeded values of 0.5 (during theta, ripples, and induced ripples) indicating relatively loose self- and cross-class consistency of multineuronal sequences, the observations are robust to the choice of analysis method.
Discussion
We found sequential activation of CA1 neuronal populations during theta oscillations, spontaneous ripples, and optogenetically induced ripple oscillations. Both firing rates and sequences were correlated across theta and spontaneous ripples. The participating neurons overlapped during theta, spontaneous and synthetic ripples, with high rate correlations across these states. Moreover, sequential ordering of neurons was also correlated across native and optogenetically induced events. These results indicate that mechanisms confined to the CA1 region contribute to both ripple- and theta-associated neuronal sequences. We suggest that sequential spiking patterns observed in the CA1 region in the intact animal arise from the interplay between extrinsic inputs (from CA3 and the entorhinal cortex), local synaptic interactions, and intrinsic cellular conductances.
The sequential ordering of pyramidal neurons during sharp wave ripples often faithfully reflects the sequential patterning of place cells before immobility and sleep. The “replay” of waking activity has been assumed to play a critical role in memory consolidation and transfer of information from the hippocampus to the neocortex (12, 13, 15–18), because selective elimination of ripples after task learning impairs memory performance (22–24). Recent studies also emphasized the prospective and constructive role of ripple sequences (16–19). Sequential activity of CA1 pyramidal neurons during ripples and in memory tasks during theta oscillations has been implicitly assumed to be generated in the recurrent excitatory networks of the hippocampal CA3 region and inherited through synaptic activity by the CA1 neuronal population (4, 12, 16, 21, 23) although direct evidence has not been shown experimentally. Our findings demonstrate that the CA1 region also contributes to the generation of these forward and reversed sequences.
Recently, we demonstrated that local optogenetic stimulation of CA1 pyramidal neurons is sufficient to induce fast LFP oscillations with many properties similar to spontaneously occurring ripples (26). Sequential activation of neurons during such artificial events by itself is not surprising: one can assume that the distribution of ChR2 is not equal across neurons and that light scattering from the optic fiber tip to the neuronal tissue is unequal. Thus, one might expect that optogenetic induction of synthetic ripples would be highly stereotyped. However, under the hypothesis of heterogeneous ChR2 distribution and unequal light distribution, no relationship of firing patterns is expected between spontaneous and synthetic ripples. In contrast, we observed strong firing rate correlations of both pyramidal cells and putative interneurons between baseline activity, spontaneous and induced ripples. Furthermore, phase-locking of neurons to ripple waves was also correlated. Finally, the temporal ordering of pyramidal cells and interneurons during theta oscillations, spontaneous and induced ripples was weakly yet consistently correlated.
When only sequences of pyramidal neurons were considered, the correlation between spontaneous and induced oscillations was not significant. However, this was also the case when sequence patterns were correlated within the same conditions. Interneurons were embedded in the assembly sequences, similar to a previous study examining the stability of sequential activation of pyramidal cells and interneurons during slow oscillations in the neocortex (34). Single interneurons are typically controlled by a limited number of pyramidal cell assemblies (35, 36) and suppress the competing assemblies, which can explain their orderly presence in the sequential patterns. In turn, interneuron spiking may provide a scaffold – a temporal backbone – for pyramidal cell spiking. It may well be that many more pyramidal neurons took part in the sequential patterning of activity than those monitored here by the silicon probes; recordings of larger ensembles of pyramidal cells in future experiments are needed to address this hypothesis.
How can correlations between internally and externally driven events arise? One clue lies in the wide dynamic range of the excitability of pyramidal cells and the spike transfer between pyramidal cells and fast spiking interneurons (37). An expected outcome of such a wide distribution of intrinsic properties is that the participation probability and temporal patterning of CA1 neurons during spontaneous ripples can be only partially dictated by the upstream CA3 neurons, and that the intrinsic properties bias the participation probability and sequential recruitment of neurons in response to optogenetic activation. In line with this reasoning, pyramidal neurons vary extensively in their responses to in vitro current injection (38, 39). Thus, even if every neuron experienced an identical ramp of input current (from the CA3 population during spontaneous ripple events or from the light stimulus during synthetic ripples), diverse spike patterns and sequential activity would still ensue due to intrinsic neuronal/circuit heterogeneity.
The similar participation probability and sequential patterning of putative interneurons during both spontaneous and induced ripples indicate that, in addition to the hypothesized biophysical heterogeneity of neurons, local synaptic interactions are also important. Sequential firing of interneuron spikes, intermingled with pyramidal neurons during spontaneous ripples, is not unexpected given the timing relations observed among the distinct interneuron classes during ripples (28–30). Although most interneurons receive predominantly local input, others also receive synaptic input from CA3 and extrahippocampal sources (40). Our findings demonstrate that similar sequential patterning of CA1 pyramidal cells and interneurons can be induced by local mechanisms and somewhat independently of the nature of triggering inputs.
We hypothesize that the intersection of three sources of mechanisms explains the observed spike correlations between spontaneous and synthetic ripples as well as their variability: sequential input drive from CA3, wide distribution of intrinsic cellular properties, and differences in local connectivity (the last two leading to unequal responses to extrinsic activation). The preexisting sequences may themselves reflect the local circuit configurations shaped by previous experiences, by the developmental genetic program, or any combination thereof. This hypothesis does not reduce the importance of spatiotemporally structured input from CA3 and experience-induced synaptic changes, yet emphasizes that intrinsic cellular properties and local synaptic interactions in the CA1 region (41) also play an important role in shaping the sequential firing of CA1 neurons. Sequentially ordered firing patterns are not unique to the hippocampus but also occur in the neocortex during the UP state of slow oscillations (34, 42), suggesting that intrinsic neuronal and local-circuit properties may play a similar role in other systems as well.
Our findings demonstrate that circuit and intrinsic channel properties of neurons – and in particular, interactions between principal cells and interneurons – play a role in determining the probability by which particular neurons are drawn for a given task (15) and may influence their sequential activation during both theta and ripple oscillations in the intact animal. Previous work already indicated that spike rates of pyramidal cells and interneurons in CA1 (37, 43), and of head direction neurons in the thalamus and postsubiculum (44), are robustly correlated across brain states. Correlated sequences and pair-wise correlations between waking exploration and sleep are often taken as evidence for experience-induced synaptic plasticity (13–15, 45). However, preexisting correlations have also been reported between sleep epochs that precede and follow a novel experience (“preplay;” refs. 43 and 46). These observations, combined with the present findings using optical activation of hippocampal circuits, indicate that learning-induced effects should be carefully separated from preexisting conditions biased by the intrinsic properties of neurons and the local interactions between pyramidal cells and interneurons (41–43, 46, 47).
Methods
Ten mice were used in this study; nine expressed ChR2 in CA1 pyramidal cells under the CaMKII promoter, and one was opsin-free; all mice were implanted with multisite diode probes (48). Each probe had four or six shanks (intershank spacing, 200 μm) with 8 or 10 sites per shank (intersite spacing, 20 μm), and each shank was associated with an optical fiber coupled to a miniature LED or laser diode. Light stimulation (50- to 70-ms light pulses) was applied on every shank separately during spontaneous behavior in the home cage (n = 4 mice) or during urethane anesthesia (n = 6 mice; ref. 49). Offline, spikes were sorted and classified as putatitve PYR or INT using a Gaussian mixture model (49). Spike sequences were detected and compared using a cross-validated template method or an alternate pair-wise method; statistical significance was verified by resampling (permutation test). An α level of 0.05 was used throughout. All animal handling procedures were approved by the New York University Animal Care and Facilities committee. For more details, see SI Methods.
SI Methods
Animals and Surgery.
Four freely moving and six urethane anesthetized mice were used in this study. All animals expressed ChR2 under the CaMKII promoter. In two of the freely moving animals, virus (rAAV5/CaMKIIa-hChR2(h134R)-EYFP; viral titer estimated at 4 × 1012 IU/mL; University of North Carolina viral core facility, courtesy of K. Deisseroth) was injected stereotactically into CA1 (AP -1.6, ML 1.1, DV 1.1 ± 0.2 mm; 55 nl/site; Nanoject II, Drummond). The other animals were cross-bred transgenics [B6.Cg-Tg(Camk2a-cre)T29-1Stl female, Jackson Labs 005359; cross-bred with B6;129S-Gt(ROSA)26Sortm32(CAG-COP4*H134R/EYFP)Hze/J male, Jackson Labs 012569; five mice] or opsin-free (C57BL/6J, Jackson Labs; one mouse). In all animals, a craniotomy was made above the right hemisphere (AP −1.6 and ML 1.1 mm), and a 32-site/four-shank probe (Buzsaki32, NeuroNexus) equipped with four blue and/or violet diodes (48) was mounted on a movable microdrive and implanted in the neocortex under isoflurane anesthesia. After recovery from anesthesia, subsequent downward movements of the probe were made in 35- to 70-μm increments over 12- to 24-h intervals until approaching the CA1 pyramidal cell layer, recognized by the appearance of multiple high-amplitude units and spontaneous ripple events. After implantation, animals were kept one to a cage on a reversed light/dark cycle. Experiments with urethane-anesthetized animals (1.5 g/kg) utilitized 64-site, six-shank/six-LED diode-probes and followed similar procedures (49). Peak light power, measured at the tip of the shanks before implantation, was 29 ± 12 μW (mean ± SD) for 470-nm LEDs and 239 ± 61 μW for 405-nm laser diodes.
Electrophysiology and Photostimulation.
Neural activity was wide-band filtered (1–5,000 Hz, RC Electronics), amplified (20×, HST/32V-G20 headstage, Plexon; and 50×, RC Electronics), and digitized (16 bit, 20 kHz digitization, DataMax) or amplified and multiplexed on the headstage (0.3–10,000 Hz, 400×; RHA2132, Intan) and then digitized (14 bit, 20 kHz, AmpliPex). Freely moving animals were equipped with a 3-axis accelerometer (ADXL-330, Analog Devices) for monitoring head-movements, two head-mounted LEDs for online video tracking, or both. Recordings from freely moving animals were carried out in the home cage during spontaneous behavior. After each session the probe was either left in place or advanced in 35–70 μm steps and the brain was allowed to settle overnight. At each location in the brain, neuronal activity was inspected for spontaneous spiking activity, and if encountered, a baseline period of at least 15 min was recorded followed by light stimulation.
All shanks were equipped with diodes, driven by a custom-made 16-channel current source (48), controlled by a programmable DSP (25 kHz, RX6 or RX8, Tucker-Davis Technologies) via MATLAB (MathWorks). Photostimulation was performed at each shank separately at the minimal light intensity that evoked an effect detectable by visual inspection during the experiment, and at multiple intensities above and below that level. All stimuli were a-periodic and had a square pulse (50-70 ms) envelope.
Spike Sorting and Ripple Detection.
For offline analysis, spike waveforms were extracted from the wide-band recorded signals and sorted into individual units (49). Briefly, waveforms were linearly detrended, projected onto a common basis obtained by principal component analysis of the data, and sorted automatically followed by manual adjustment. Only well-isolated units (amplitude >50 μV; l-ratio <0.05; interspike interval index <0.2) were used. Subsequently each unit was tagged as excitatory/inhibitory [based on peaks/troughs in the short-time (±5 ms) pair-wise cross correlation; P < 0.001, convolution test] and/or classified as putative PYR or INT (based on a Gaussian-mixture model; P < 0.05; ref. 49). We recorded a total of 582 well isolated cells from CA1 of four freely moving and six anesthetized mice. Of these, 496 were PYR and 86 were INT.
Ripples were detected independently at each recording site (26). Briefly, the wide-band signal was band-pass filtered (80–250 Hz), and instantaneous power was computed by clipping extreme values to 5 SD, rectifying, and low-pass filtering (52.5 Hz). The mean and SD were computed from the power of the clipped signal during slow-wave sleep (defined as nontheta, nonmovement periods) in the absence of light stimulation. Subsequently, the power of the original trace was computed, and all events exceeding 5 SD from the mean were selected. Short events (duration < 15 ms) were discarded, and adjacent events (gap < 15 ms) were merged. Events were then expanded until the power fell below 2 SD and aligned by the trough (of the nonrectified signal) closest to the peak power. For each shank separately, the site with the maximal ripple amplitude was determined and defined as the center of the CA1 pyramidal cell layer.
Sequence Determination.
For the analysis of sequences, only sites in which at least 100 spontaneous ripples were observed and at least 10 optical stimuli were applied at threshold intensity were used, yielding a total of 56 distinct sites (27 in freely moving mice; 29 in urethane-anesthetized mice) recorded during 16 sessions (one to three sessions per mouse). Threshold intensity (defined as the minimal light intensity that induced synthetic ripples with amplitude not smaller than that observed during spontaneous ripples; ref. 26) was used, because at higher intensities many spikes are superimposed near the induced ripple cycle troughs and cannot be disambiguated using extracellular recording methods.
For each event class (theta cycle, spontaneous ripple, or induced ripple), all spiking units were considered, both those recorded on the same shank as the ripple and on other shanks, because interneurons on distant shanks modulate their timing during spatially confined ripple events (26). Only events with at least two spiking units were considered. For each specific event, the spiking units were rank ordered chronologically. For instance, if the considered units are 1, 2, 3, 4, and 5, and during a given event units 2 and 4 spiked, 4 first, then unit 4 received rank = 1 and unit 2 received rank = 2. If a given unit spiked more than once during an event, only the first spike was considered. Two methods were used: template and pair-wise. In the cross-validated template method, a template for each specific ripple or theta event was constructed based on the averaged rank of all units over all other events (i.e., excluding that specific event). Then the rank correlation between each specific event and its template was computed, and averaged over all events; this includes n comparisons (n being the number of valid events). In the pair-wise method, there is no template; rather, the rank correlation is computed between each pair of events (n × (n − 1)/2 comparisons).
For comparing spike sequences during spontaneous and induced events (cross-class consistency), the same two methods were used. In the template method, the template was based on the spontaneous events, whereas the events tested were the induced events. In the pair-wise method, rank correlations were computed between all possible spontaneous and synthetic ripple events (i.e., spontaneous–spontaneous and induced–induced correlations were not considered). Thus, in this application, the template method is asymmetric, whereas the pair-wise method is symmetric (i.e., the null hypothesis is the same, with a different alternative).
Statistical significance of sequence consistency was determined by a permutation test, randomly shuffling the ranks of the units that participated in each event and carrying out the entire procedure, culminating in the computation of the mean rank correlation. This was repeated 1,000 times and the observed mean rank correlation was compared with the distribution of the randomly permuted mean rank correlations (Fig. 4B).
Statistical Testing.
For a given effect size, the power of any statistical test depends on the α level. To increase the sensitivity of detecting effects, results are reported based on a significance threshold α = 0.05, and all groups included enough samples to enable rejection of the null at that level. Wilcoxon’s signed-rank test was used for paired testing (or, one-sided, for a zero-median null) and Mann–Whitney U test (two-sided) was used for comparing the medians of unpaired groups. For proportions, the exact Binomial test (single comparisons) and the χ2 test (pairwise or multiple groups) were used. Circular measures were tested using the χ2 or the Rayleigh test. All results were verified independently by resampling (permutation) tests.
Acknowledgments
We thank Karl Deisseroth for providing adeno-associated viruses, Kamran Diba for insightful comments, and Shigeyoshi Fujisawa for help maintaining transgenic animals. This work was supported by the Rothschild Foundation (to E.S.), Human Frontiers in Science Program LT-000346/2009-L (to E.S.), Machiah Foundation 20090098 (to E.S.), National Institutes of Health Grants (MH54671; MH102840; NS074015), the Mathers Foundation, The Human Frontier Science Program, and the National Science Foundation (Temporal Dynamics of Learning Center Grant SBE 0542013).
Supporting Information
Supporting Information (PDF)
Supporting Information
- Download
- 1.03 MB
References
1
DO Hebb The Organization of Behavior (Wiley & Sons, New York, 1949).
2
Y Prut, et al., Spatiotemporal structure of cortical activity: Properties and behavioral relevance. J Neurophysiol 79, 2857–2874 (1998).
3
S Fujisawa, A Amarasingham, MT Harrison, G Buzsáki, Behavior-dependent short-term assembly dynamics in the medial prefrontal cortex. Nat Neurosci 11, 823–833 (2008).
4
E Pastalkova, V Itskov, A Amarasingham, G Buzsáki, Internally generated cell assembly sequences in the rat hippocampus. Science 321, 1322–1327 (2008).
5
CD Harvey, P Coen, DW Tank, Choice-specific sequences in parietal cortex during a virtual-navigation decision task. Nature 484, 62–68 (2012).
6
M Abeles Corticonics: Neural Circuits of the Cerebral Cortex (Cambridge University Press, New York, 1991).
7
V Itskov, C Curto, E Pastalkova, G Buzsáki, Cell assembly sequences arising from spike threshold adaptation keep track of time in the hippocampus. J Neurosci 31, 2828–2834 (2011).
8
J O’Keefe, L Nadel The Hippocampus as a Cognitive Map (Oxford University Press, Oxford, 1978).
9
J O’Keefe, N Burgess, Geometric determinants of the place fields of hippocampal neurons. Nature 381, 425–428 (1996).
10
M Fyhn, S Molden, MP Witter, EI Moser, MB Moser, Spatial representation in the entorhinal cortex. Science 305, 1258–1264 (2004).
11
CJ MacDonald, KQ Lepage, UT Eden, H Eichenbaum, Hippocampal “time cells” bridge the gap in memory for discontiguous events. Neuron 71, 737–749 (2011).
12
G Buzsáki, Two-stage model of memory trace formation: A role for “noisy” brain states. Neuroscience 31, 551–570 (1989).
13
MA Wilson, BL McNaughton, Reactivation of hippocampal ensemble memories during sleep. Science 265, 676–679 (1994).
14
Z Nádasdy, H Hirase, A Czurkó, J Csicsvari, G Buzsáki, Replay and time compression of recurring spike sequences in the hippocampus. J Neurosci 19, 9497–9507 (1999).
15
AK Lee, MA Wilson, Memory of sequential experience in the hippocampus during slow wave sleep. Neuron 36, 1183–1194 (2002).
16
DJ Foster, MA Wilson, Reverse replay of behavioural sequences in hippocampal place cells during the awake state. Nature 440, 680–683 (2006).
17
K Diba, G Buzsáki, Forward and reverse hippocampal place-cell sequences during ripples. Nat Neurosci 10, 1241–1242 (2007).
18
MP Karlsson, LM Frank, Awake replay of remote experiences in the hippocampus. Nat Neurosci 12, 913–918 (2009).
19
TJ Davidson, F Kloosterman, MA Wilson, Hippocampal replay of extended experience. Neuron 63, 497–507 (2009).
20
T Pfeiffer, A Draguhn, S Reichinnek, M Both, Optimized temporally deconvolved Ca²⁺ imaging allows identification of spatiotemporal activity patterns of CA1 hippocampal ensembles. Neuroimage 94, 239–249 (2014).
21
G Buzsáki, Z Horváth, R Urioste, J Hetke, K Wise, High-frequency network oscillation in the hippocampus. Science 256, 1025–1027 (1992).
22
G Girardeau, K Benchenane, SI Wiener, G Buzsáki, MB Zugaro, Selective suppression of hippocampal ripples impairs spatial memory. Nat Neurosci 12, 1222–1223 (2009).
23
V Ego-Stengel, MA Wilson, Disruption of ripple-associated hippocampal activity during rest impairs spatial learning in the rat. Hippocampus 20, 1–10 (2010).
24
SP Jadhav, C Kemere, PW German, LM Frank, Awake hippocampal sharp-wave ripples support spatial memory. Science 336, 1454–1458 (2012).
25
D Sullivan, et al., Relationships between hippocampal sharp waves, ripples, and fast gamma oscillation: Influence of dentate and entorhinal cortical activity. J Neurosci 31, 8605–8616 (2011).
26
E Stark, et al., Pyramidal cell-interneuron interactions underlie hippocampal ripple oscillations. Neuron 83, 467–480 (2014).
27
AM Thomson, S Radpour, Excitatory connections between CA1 pyramidal cells revealed by spike triggered averaging in slices of rat hippocampus are partially NMDA receptor mediated. Eur J Neurosci 3, 587–601 (1991).
28
J Csicsvari, H Hirase, A Czurkó, A Mamiya, G Buzsáki, Oscillatory coupling of hippocampal pyramidal cells and interneurons in the behaving Rat. J Neurosci 19, 274–287 (1999).
29
T Klausberger, P Somogyi, Neuronal diversity and temporal dynamics: The unity of hippocampal circuit operations. Science 321, 53–57 (2008).
30
C Varga, P Golshani, I Soltesz, Frequency-invariant temporal ordering of interneuronal discharges during hippocampal oscillations in awake mice. Proc Natl Acad Sci USA 109, E2726–E2734 (2012).
31
TF Freund, G Buzsáki, Interneurons of the hippocampus. Hippocampus 6, 347–470 (1996).
32
G Fishell, B Rudy, Mechanisms of inhibition within the telencephalon: “where the wild things are”. Annu Rev Neurosci 34, 535–567 (2011).
33
JJ Chrobak, G Buzsáki, Selective activation of deep layer (V-VI) retrohippocampal cortical neurons during hippocampal sharp waves in the behaving rat. J Neurosci 14, 6160–6170 (1994).
34
A Luczak, P Barthó, SL Marguet, G Buzsáki, KD Harris, Sequential structure of neocortical spontaneous activity in vivo. Proc Natl Acad Sci USA 104, 347–352 (2007).
35
AP Maurer, SL Cowen, SN Burke, CA Barnes, BL McNaughton, Phase precession in hippocampal interneurons showing strong functional coupling to individual pyramidal cells. J Neurosci 26, 13485–13492 (2006).
36
C Geisler, D Robbe, M Zugaro, A Sirota, G Buzsáki, Hippocampal place cell assemblies are speed-controlled oscillators. Proc Natl Acad Sci USA 104, 8149–8154 (2007).
37
K Mizuseki, G Buzsáki, Preconfigured, skewed distribution of firing rates in the hippocampus and entorhinal cortex. Cell Reports 4, 1010–1021 (2013).
38
T Jarsky, R Mady, B Kennedy, N Spruston, Distribution of bursting neurons in the CA1 region and the subiculum of the rat hippocampus. J Comp Neurol 506, 535–547 (2008).
39
KA Dougherty, T Islam, D Johnston, Intrinsic excitability of CA1 pyramidal neurones from the rat dorsal and ventral hippocampus. J Physiol 590, 5707–5722 (2012).
40
TF Freund, M Antal, GABA-containing neurons in the septum control inhibitory interneurons in the hippocampus. Nature 336, 170–173 (1988).
41
SH Lee, et al., Parvalbumin-positive basket cells differentiate among hippocampal pyramidal cells. Neuron 82, 1129–1144 (2014).
42
A Luczak, P Barthó, KD Harris, Spontaneous events outline the realm of possible sensory responses in neocortical populations. Neuron 62, 413–425 (2009).
43
H Hirase, X Leinekugel, A Czurkó, J Csicsvari, G Buzsáki, Firing rates of hippocampal neurons are preserved during subsequent sleep episodes and modified by novel awake experience. Proc Natl Acad Sci USA 98, 9386–9390 (2001).
44
A Peyrache, MM Lacroix, PC Petersen, G Buzsáki, Internally organized mechanisms of the head direction sense. Nat Neurosci 18, 569–575 (2015).
45
M Mizunuma, et al., Unbalanced excitability underlies offline reactivation of behaviorally activated neurons. Nat Neurosci 17, 503–505 (2014).
46
G Dragoi, S Tonegawa, Preplay of future place cell sequences by hippocampal cellular assemblies. Nature 469, 397–401 (2011).
47
HS Kudrimoti, CA Barnes, BL McNaughton, Reactivation of hippocampal cell assemblies: Effects of behavioral state, experience, and EEG dynamics. J Neurosci 19, 4090–4101 (1999).
48
E Stark, T Koos, G Buzsáki, Diode probes for spatiotemporal optical control of multiple neurons in freely moving animals. J Neurophysiol 108, 349–363 (2012).
49
E Stark, et al., Inhibition-induced theta resonance in cortical circuits. Neuron 80, 1263–1276 (2013).
Information & Authors
Information
Published in
Classifications
Submission history
Published online: August 3, 2015
Published in issue: August 18, 2015
Keywords
Acknowledgments
We thank Karl Deisseroth for providing adeno-associated viruses, Kamran Diba for insightful comments, and Shigeyoshi Fujisawa for help maintaining transgenic animals. This work was supported by the Rothschild Foundation (to E.S.), Human Frontiers in Science Program LT-000346/2009-L (to E.S.), Machiah Foundation 20090098 (to E.S.), National Institutes of Health Grants (MH54671; MH102840; NS074015), the Mathers Foundation, The Human Frontier Science Program, and the National Science Foundation (Temporal Dynamics of Learning Center Grant SBE 0542013).
Notes
This article is a PNAS Direct Submission. H.E. is a guest editor invited by the Editorial Board.
Authors
Competing Interests
The authors declare no conflict of interest.
Metrics & Citations
Metrics
Citation statements
Altmetrics
Citations
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.
Cited by
Loading...
View Options
View options
PDF format
Download this article as a PDF file
DOWNLOAD PDFLogin options
Check if you have access through your login credentials or your institution to get full access on this article.
Personal login Institutional LoginRecommend to a librarian
Recommend PNAS to a LibrarianPurchase options
Purchase this article to access the full text.