Multineuronal spike sequences repeat with millisecond precision

Introduction

The brain uses a limited number of neurons to process virtually unlimited patterns of information from external environments. Therefore, individual neurons are unlikely to independently process specific information, and it is more plausible that they cooperatively form subgroups that work as functional units. This idea, called the “cell assembly” hypothesis (Hebb, 1949; Harris, 2005; Buzsaki, 2010), leads to two important predictions about neuronal circuit operation. First, a given neuron can participate in two or more cell assemblies. Second, the synapse weight between two given neurons is modifiable over time. More specifically, the weight is strengthened when the neurons work cooperatively, otherwise it is weakened. This bidirectional synaptic plasticity has been experimentally proven by studies showing that cortical synapses are capable of exhibiting long-term potentiation (LTP) and long-term depression (LTD) of synaptic transmission. These two features suggest that neuronal networks self-organize through reorganization of neuronal connectivity due to ongoing external stimuli and are thereby functionally compartmentalized to form cell assemblies.

Synchronous activity among multiple neurons is regarded as one of the simplest aspects of cell assembly dynamics, not only because it triggers the induction of synaptic plasticity but also because it is realized through synchronization between groups of presynaptic neurons (Takahashi et al., 2010). At a more microscopic level, synchronized activity is often composed of sequential activation of multiple neurons. Indeed, such multineuronal spike sequences (MSSs) are known to exist during sensory-evoked and spontaneous network activity at frequencies greater than chance (Abeles and Gerstein, 1988; Prut et al., 1998; Abeles and Gat, 2001; Ikegaya et al., 2004; Luczak et al., 2007). MSSs are conceptually consistent with the so-called “synfire chains” hypothesis, in which neurons form several layers through which synchronized spikes can stably propagate with temporal precision, therefore generating a chain of spikes across neurons (Abeles, 1991). Although MSSs have been reported in different brain preparations, including the neocortex of monkeys and rodents in vivo and in vitro, some reports still cast doubt on the existence of MSSs (Oram et al., 1999; Baker and Lemon, 2000; Mokeichev et al., 2007).

We focus the present work on the hippocampus. One of the reasons for this selection is that the hippocampus spontaneously emits sharp waves-ripples (SW-Rs), a transient form of high-frequency field oscillations. SW-Rs primarily occur during slow-wave sleep and quiet awake states (Buzsaki et al., 1983). During SW-Rs, neurons that were previously involved in behavioral exploration, called place cells, increase their firing rates sequentially in the same or opposite order of which those neurons were activated during the behavioral exploration (Lee and Wilson, 2002; Harris et al., 2003; Foster and Wilson, 2006; O'Neill et al., 2006; Diba and Buzsaki, 2007; Pastalkova et al., 2008; Davidson et al., 2009). Therefore, the hippocampus may be a good model to study MSSs. In vitro slices prepared from the hippocampus are also reported to emit spontaneous SW-Rs (Norimoto et al., 2012; Sun et al., 2012), and the slices may produce MSSs. If it were the case, it would be easy to analyze the properties of MSSs because the in vitro experimental system is more accessible and manipulatable than the in vivo system.

Functional multineuron calcium imaging (fMCI) is an optical technique that monitors spikes of neurons in situ through spike-evoked somatic calcium transients. We have optimized the method for hippocampal slice cultures (Takahashi et al., 2007, 2011). In the CA3 region of organotypically cultured slices, the probability of a synaptic connection between randomly selected adjacent pyramidal cells is approximately 25%. This ratio is higher than that reported in acute slice preparations (2–8%) (Miles and Wong, 1983; Smith et al., 1995). In acute hippocampal slices, however, 75–90% of the axons of CA3 pyramidal neurons are cut, even in slices as thick as 500 μm (Gomez-Di Cesare et al., 1997). In contrast, cultured networks self-restore their complexity to a realistic extent. Indeed, levels of spontaneous excitatory or inhibitory postsynaptic currents are similar between ex vivo and in vivo hippocampal neurons (Takahashi et al., 2010). Moreover, we have demonstrated, by using an optical synapse mapping technique, that in such restored CA3 networks, pyramidal cells are nonrandomly connected to generate diverse repertoires of synchronized activity, like the in vivo conditions (Takahashi et al., 2010).

Materials and Methods

Animal

Experiments were performed with the approval of the animal experiment ethics committee at the University of Tokyo (approval number: 19-43, P21-6) according to the University of Tokyo guidelines for the care and use of laboratory animals.

Slice Culture Preparations

Entorhino-hippocampal organotypic slices were prepared from 7-day-old Wistar/ST rats (SLC, Shizuoka, Japan) as previously described (Koyama et al., 2007). Briefly, rat pups were anesthetized by hypothermia and decapitated. The brains were removed and placed in aerated, ice-cold Gey's balanced salt solution supplemented with 25 mM glucose. Horizontal entorhino-hippocampal slices were made at a thickness of 300 μm using a vibratome (DTK-1500, Dosaka, Kyoto, Japan). The slices were placed on Omnipore membrane filters (JHWP02500, Millipore, Bedford, MA, USA) and incubated in 5% CO₂ at 37°C. The culture medium, which was composed of 50% minimal essential medium (Invitrogen, Gaithersburg, MD, USA), 25% Hanks' balanced salt solution, 25% horse serum (Cell Culture Laboratory, Cleveland, OH, USA), and antibiotics, was changed every 3.5 days. Experiments were performed on days 7–11 in vitro.

Although slice cultures are known to form abnormal connections that very rarely exist under normal conditions, such as CA1-to-CA1, CA1-to-CA3, and CA3-to-dentate gyrus connections (Gahwiler et al., 1997; De Simoni et al., 2003), there are few of these aberrant connections in our slice culture preparations. Investigation using reverse optical trawling, a synapse mapping technique (Sasaki et al., 2009), demonstrated that these abnormal connections are less than 0.5% of the total connections and that an overwhelming number of connections project to their normal targets. This result is most likely because the entorhinal cortex was not dissected out in our preparations. Lesions of the entorhinal cortex are known to result in abnormal sprouting and reorganization of hippocampal networks in vivo and ex vivo (Laurberg and Zimmer, 1981; West and Dewey, 1984).

Functional Multineuron Calcium Imaging

Slices were incubated with 2 ml of dye solution at 37°C for 1 h (Takahashi et al., 2011). The dye solution was aCSF containing 0.0005% Oregon Green 488 BAPTA-1AM (OGB-1AM), 0.01% Pluronic F-127, and 0.005% Cremophor EL. After a 1-h recovery, a slice was transferred to a recording chamber. Images were acquired at 500 frames/s with a Nipkow-disk confocal unit (CSUX-1, Yokogawa Electric, Tokyo, Japan), a back-illuminated electron-multiplying charge-coupled device (EM-CCD) camera (iXon DV860, Andor, Belfast, Northern Ireland, UK), a water-immersion objective lens (16×, 0.80 NA, Nikon, Tokyo, Japan), and Solis software (Andor). Fluorophores were excited at 488 nm with an argon laser (10–15 mW, 532-BS-AO4, Omnichrome, Chino, CA, USA) and visualized with a 507-nm long-pass emission filter. In each cell body, the fluorescence change ΔF/F was calculated as (F_t–F₀)/F₀, where F_t is the fluorescence intensity at frame time t, and F₀ is baseline (Figure 1A). Spike timing was defined as the onset of individual calcium transients with an automatic machine-learning algorithm that can accurately detect spike timing within one frame jitter (Sasaki et al., 2008). In some experiments, picrotoxin was bath-applied to prevent fast inhibitory synaptic transmission. Picrotoxin (purchased from Sigma-Aldrich, St. Louis, MO) was dissolved in aCSF at the final concentration of 50 μM and perfused to slices. Imaging was started 30 min after the perfusion onset.

FIGURE 1

Figure 1. High-speed fMCI of spontaneous CA3 network activity. (A) Top: representative confocal image of the CA3 stratum pyramidale in an OGB1-loaded hippocampal slice culture. Bottom: Time changes in the OGB1 fluorescence intensity in the somata of 7 randomly selected neurons indicated by the colored circles in the top photograph. (B) Representative rastergram of all 110 CA3 neurons (top) and time histogram of the percentage of active neurons (bin = 10 ms). A large synchronization is time-expanded in the inset and is displayed together with 150–300 Hz band-passed LFPs simultaneously recorded from the CA1 pyramidal stratum. Synchronization was accompanied by SW-R-like oscillations.

Local Field Potential Recordings and Ripple Detection

In a single experiment, CA1 local field potentials (LFP) were recorded during fMCI monitoring of the calcium activity of CA3 neurons. Glass pipettes were filled with 2 M NaCl and placed in CA1 stratum pyramidale. To extract the SW-R activity, the recorded data were band-pass filtered at 150–300 Hz. SW-R-like events were automatically detected based on their oscillatory powers and durations; the root mean square (3-ms window) of the band-passed signal was used to detect SW-Rs of 10 ms duration with a power threshold of 5 standard deviations (SDs).

MSS Detection

We used a template-matching algorithm to search for MSSs (Ikegaya et al., 2004). We first selected cells that showed more than one calcium transient. After determining the reference calcium levels of a reference cell (cell₁), we designated a vector consisting of a set of cells and relative timings of their calcium events as follows: (cell₂, …, cell_N, t₂, t₃, …, t_N), where t_i denotes the delay of the event in cell_i after the reference event. t_i was limited to less than 500 ms. This vector was used as a template and was slid along the successive events of cell₁ throughout the recording session. If more than two elements were identical between any template pairs, we regarded the matched elements as an MSS. Each mismatched spike configuration was used as another template in a subsequent scan. Thus, every event was considered part of a template MSS, and each template occurred at least once. Unless otherwise specified, one frame jitter (2 ms) was allowed so that the total number of MSSs satisfied statistical demands (Kendall et al., 1994).

Surrogate Data

To determine whether MSSs or their structures can arise from a stochastic process, we created surrogate raster plots using a Monte Carlo resampling method (Ikegaya et al., 2004). Each calcium event was exchanged between a pair of cells, maintaining their relative timings (Figure 3A). This procedure was repeated for all calcium transients in each raster plot. This randomization preserves the event frequencies of individual cells and the population modulation of event timings, such as network synchronization. In each shuffled surrogate, we searched MSSs using the same detection algorithm. Twenty surrogates were generated for each dataset, and the averages across the 20 surrogates were defined as the chance level.

Asynchrony Index

The temporal sparseness of network activity during the observation period was captured by a normalized Shannon index, termed asynchrony index (Usami et al., 2008; Mizunuma et al., 2009; Ujita et al., 2011). The Shannon index quantifies the dispersion of components in a histogram and is generally defined as −∑_i (k_i/K) log₂ (k_i/K), where K is the total number of components, and k_i is the number of components in the i-th bin. This definition of diversity is conceptually equivalent to Shannon's entropy. Because this index is sensitive to K and the bin size, SI has often been normalized with the maximal value and other standard values to compare groups. Here we normalized Shannon index with the maximal (SImax) and minimal values (SImin) that can be taken in the distribution of the same number of spikes in the raster plot. SImax and SImin were obtained through data shuffle with maintaining K and bin; SImax is given when components are as evenly redistributed over the time axis as possible, whereas SImin is given when components are as temporally biased as possible. Then normalized Shannon index is defined as (SI–SImin)/(SImax–SImin). Thus, it takes a value from 0 to 1, with higher values being more dispersive.

Data Representation

We reported all averaged values as the mean ± SDs.

Results

High-Speed Imaging of Spiking CA3 Networks ex vivo

Hippocampal slice cultures were incubated in OGB1AM, and OGB1-loaded neurons were imaged from the CA3 stratum pyramidale at 500 Hz using a spinning-disk confocal microscope and a high-speed EM-CCD camera (Takahashi et al., 2011). The microscopic field covered an area of approximately 350 × 200 μm (Figure 1A) and contained an average of 91.7 ± 26.7 neurons (mean ± SD of 9 videos, ranging from 60 to 137 neurons). Each video was 130 s in length, and a total of 9 videos were recorded from 9 slices (n = 9 rat pups borned from 9 mothers).

In all 9 videos, spontaneous activity was evident; among a total of 825 neurons in 9 videos, 757 neurons (91.8%) exhibited at least one spike during the observation period. The mean firing rates of active neurons were 0.25 ± 0.29 Hz (n = 9 videos), ranging from 0.008 to 2.33 Hz. Therefore, spontaneous activity was sparse as a whole. Nevertheless, neurons occasionally exhibited synchronization at the population level (Figure 1B); synchronization that recruited more than 5% or 10% of the total neurons occurred at frequencies of 2.46 ± 3.0 or 0.26 ± 0.52 per min, respectively (n = 9 videos, bin = 2 ms = one frame).

We succeeded in simultaneous LFP recording in one single video; note that in general, LFP recording is technically difficult in slice cultures, most likely because neurite reorganization during cultivation slightly alters the fine layer structure of the hippocampus, collapsing the net dipole moment generated by synaptic activity or because cultured networks may be heterogenous from preparation to preparation (e.g., see Figure 4). In the single dataset, we found that synchronous activity was always accompanied by SW-R-like high-frequency oscillations (Figure 1B, inset).

MSSs with Millisecond Precision

We detected MSSs. The maximal length of an MSS, i.e., the time interval between the first spike and the last spike, was set to be 500 ms. Three examples of MSSs found in a raster plot are shown in Figure 2. Not surprisingly, the number of MSSs depended on the spike time jitters allowed to detect MSSs (Figure 3A left). MSSs increased in number when jitters were increased from 0 ms (0 frame) to 32 ms (16 frames); note that 0-ms jitter (= 0 frames) had a time window of 2 ms due to 500-Hz imaging.

FIGURE 2

Figure 2. Examples of MSSs. Three representative MSSs are merged in a single rastergram and magnified in time in the bottom panels.

FIGURE 3

Figure 3. Monte Carlo evaluation for the significance of MSSs. (A) Left: the total numbers of MSSs in the rastergram shown in (A) are plotted as a function of spike-time jitters allowed for MSS detection. Right: MSSs were identified in the same way in spike-shuffled rastergrams. Data are the mean ± SDs of 20 surrogates. (B) Relationship between the CVs of the number of MSSs detected in surrogates and the number of surrogates generated. (C) Ratios of the numbers of MSSs in the original (real) rastergram to the mean numbers of MSSs in 20 randomized surrogates. Data are the same as (B).

The total number of MSSs varied among videos, probably because the network states fluctuate over time and are not identical among preparations (Sasaki et al., 2007). The total number of MSSs in each video (2-ms jitter allowed) was plotted against the total number of neurons recorded in the corresponding videos (Figure 4A), the mean activity frequency of those neurons (Figure 4B), and the asynchrony index of the corresponding videos (Figure 4C); note that the asynchrony index is a normalized Shannon index for spike dispersion along the time axis (Usami et al., 2008; Mizunuma et al., 2009; Ujita et al., 2011). Thus, we investigated the statistical significance of MSSs by comparing the number of MSSs to its stochastic level. To estimate the stochastic level, we created surrogate raster plots by randomly exchanging spikes across neurons. Specifically, a single spike of a randomly selected neuron was swapped with a spike of another randomly selected neuron without changing their absolute spike timings, and this swapping procedure was repeated until all spikes in the original dataset were exchanged (Figure 3A right, inset). This shuffling method preserved both the firing rates of individual neurons and the level and frequency of network synchronization. Therefore, comparing an original dataset with these surrogates makes it possible to examine whether MSSs are actively generated by complex network dynamics or are merely a mathematically natural consequence of the firing rates of individual neurons.

FIGURE 4

Figure 4. Across-video variations of the number of MSSs. The total number of MSSs varied among datasets and thus is plotted against the total number of neurons recorded in the corresponding videos (A), the mean frequency of activity of those neurons (B), and the synchrony index of the corresponding videos (C). Each dot indicates a single video. n = 9 videos.

The coefficient of variation (CV) for the number of MSSs detected in surrogates depended on the number of surrogates generated. The CV rapidly decreased as a function of the number of surrogates; it dropped to about 0.05 by 20 surrogates, and thereafter, it was kept almost constant (Figure 3B, 2-ms jitter allowed). To reduce computation burden with preserving statistical stability, therefore, we searched MSSs in 20 surrogates generated from each raster plot. MSSs in the surrogates also increased in number with spike jitter timing (Figure 3A right). We thus plotted the ratio of the MSS number in the real dataset to that in the surrogates as a function of jitters. The ratio peaked at a jitter of 0 ms and decreased gradually with increasing jitters (Figure 3C). We repeated this procedure for all other datasets and found similar results in the pooled data (n = 9 videos; Figure 5), that is, the mean real-to-surrogate ratio is a simple reduction function of spike jitters. Based on these results, we reached two fundamental conclusions: (1) MSSs emerge more frequently than expected by chance and therefore cannot be explained by a stochastic process, and (2) relative spike times within MSSs are maintained at the millisecond level.

FIGURE 5

Figure 5. Millisecond precision of MSSs. Summary of 9 videos. Color-edged circles are 5% significant data points between real datasets and surrogates (paired t-test). The data are the mean ± SDs.

Internal Structures of MSSs

In the following analyses, we used MSSs detected at a spike jitter of 2 ms. Although the real-to-surrogate ratios of MSSs was maximal at a jitter of 0 ms, the absolute number of MSSs detected at 0 ms was substantially smaller than at 2 ms. To enhance statistical accuracy, we selected the 2-ms jitter in this study. Note that the MSS properties described below were fundamentally consistent and exhibited similar tendencies at jitters of 0, 2, 4, and 8 ms (data not shown).

At a jitter of 2 ms, we identified a total of 27,804 MSSs in 9 videos. Among the 825 total neurons, 624 neurons (75.6%) participated in at least one MSS. On average, single MSSs contained 3.4 ± 1.1 neurons, persisted for 107.9 ± 158.6 ms (i.e., the MSS length, defined as the time interval between the first and last spike), and were repeated 2.1 ± 1.1 times during the observation periods. We compared these values to those found in surrogate datasets. The real-to-surrogate ratio of the number of neurons involved in single MSSs was consistently higher than 1 and peaked at 5–8 neurons (Figure 6A). As to the MSS length, the real-to-surrogate ratios of the number of MSSs were significantly higher than 1, and roughly in the range of 40–280 ms, peaking at 60 ms (Figure 6B). The real-to-surrogate ratios of the number of MSSs were independent of MSS repeat numbers and were consistently higher than 1 (Figure 6C).

FIGURE 6

Figure 6. Characterization of MSSs. (A–C) The number of cells that participated in individual MSSs (A), the duration that individual MSS events persisted (i.e., the time intervals between the first and last spikes within MSSs) (B), and the number of repetitions of MSSs (C) are shown as ratios of values in real datasets to the mean value of 20 randomized surrogates. Color-edged circles are 5% significant data points between real datasets and surrogates (paired t-test). The data are the mean ± SDs of 9 videos.

Two examples of MSSs (Figure 7A) are shown in a cell map (Figure 7B). We examined whether the neurons involved in each MSS were spatially clustered. For each MSS, we computed the center of gravity for the locations of all neurons that participated in the MSS and calculated the mean distance from these neurons to the center of gravity. If the MSS is spatially clustered, the mean distance will be smaller than chance. The chance values were estimated by 100 surrogates of pseudo-MSSs generated by the same number of neurons randomly selected from the original map. The mean distance from the center of gravity did not differ between real and surrogate MSSs [t_{(27, 804)} = 1.08, P = 0.28, paired t-test; Figure 7C]. Thus, MSS-participating neurons are spatially dispersed at the stochastic level.

FIGURE 7

Figure 7. Spatial dispersion of MSSs. (A,B) Two representative MSSs in a rastergram (A) and in a cell map (B). (C) The mean Euclidean distance from individual cells that participated in a MSS to the center of gravity of all cells that participated in the MSS was compared surrogates in which the same number of cells in the original MSS were selected randomly from the same cell map. This comparison was repeated for 27,804 MSSs in 9 videos, and data are shown as the mean ± SDs. (D) The number of MSSs was still overrepresented in surrogates generated by activity shuffling in which spatial configurations were preserved. **P < 0.01, paired t-test. The data are the mean ± SDs of 9 videos.

Here we returned to the statistical issues about MSS existence. Our spike shuffling method did not consider the spatial organization of cells that participated in MSSs. We thus created surrogates using another shuffling method; for each MSS, spikes were preserved for only one MSS appearance, while the other spikes constituting the MSS were randomly exchanged within the MSS. This procedure did not collapse the spatial organization of cells that exhibited MSSs. Using this shuffling method, we again confirmed that MSSs were overrepresented relative to chance (Figure 7D, n = 9 videos). However, because these Monte Carlo-shuffling evaluations did not completely abolish statistically innate problems, such as false-negative errors and false-positive errors, we also adopted a completely different approach to validate MSSs without spike shuffling; MSSs are embossed in pairwise correlograms among spike triplets occurring across three distinct cells (Luczak et al., 2007). For each cell trio, one cell was designated the trigger for calculation of the joint distribution of spike times of the other two (Figure 8A). A clear dense peak was observed at t₂–t₁ = −231 ms and t₃–t₁ = −12 ms in this distribution (n = 180 triplet spikes), suggesting that a particular sequence occurred preferentially (Figure 8B).

FIGURE 8

Figure 8. Pairwise correlogram of precisely repeating triplets. (A) For every trio of neurons, a spike triplet is described by two inter-spike intervals (t₂–t₁ and t₃–t₁). (B) Count matrix for a representative triplet of neurons, indicating the probability of different inter-spike interval combinations.

We sought to examine the timing of MSS appearance relative to the entire network activity. Figure 9A demonstrates the peri-MSS time histogram of the mean firing rates of all neurons in the video. In the peri-MSS time histogram, we aligned the first spikes in individual MSSs at time 0. Data were pooled from 27,804 MSSs in 9 videos. The histogram revealed that the neuronal network transiently increased the global firing rate during MSSs. This transient synchronization persisted for approximately 100 ms, which corresponded to the durations of SW-Rs (Buzsaki et al., 1983).

FIGURE 9

Figure 9. Mechanisms of MSSs. (A) Peri-MSS time histogram of the mean firing rates of all neurons monitored in 9 videos. The firing rates were aligned at the time relative to the onset of 27,804 individual MSSs. (B) The number of MSSs relative to chance were decreased by bath application of 50 μM picrotoxin. **P < 0.01 versus control, Student's t-test. The data are the mean ± SDs of 9 videos.

During SW-Rs, excitatory neurons and inhibitory neurons were both activated at particular phases in a cell type-specific manner (Klausberger et al., 2003). To examine the involvement of GABAergic transmission, we perfused slices with 50 μM picrotoxin, a GABA_A receptor channel inhibitor. Bath application of picrotoxin reduced MSSs, as compared to control solution (Figure 9B). Thus, MSSs do not represent simple spike chains via excitatory synapses, but rather, they are more likely to emerge actively from network-coordinated excitatory and inhibitory balance.

Biased MSS Participation of Individual Neurons

The same neurons were often recruited in different MSSs. On average, individual cells participated in 81.8 ± 282.2 MSSs; however, the frequency of MSS participation varied among neurons. The representative cell map shown in Figure 10A represents the real-to-surrogate ratios of the number of MSS participations with a pseudo-colored scale, indicating that some neurons frequently participated in MSSs, whereas others did not. In Figure 10B, we plotted the Lorenz curve of these frequency ratios (n = 825 neurons in 9 videos). The Gini coefficient was 0.44, indicating that neurons are not homogeneous in terms of MSS participation.

FIGURE 10

Figure 10. Relationship between the firing rates of neurons and their MSS participations. (A) Top: a cell map representing the number of individual cells participating in MSSs during the observation period. Data are shown as the ratio of the real rastergrams to their randomized controls. Bottom: The mean Euclidean distance to the center of gravity from individual cells with the highest x% of the real-to-surrogate ratio. The ranges of x values are indicated in the abscissa. For example, in the case of 5–10%, we selected cells with the real-to-surrogate ratios that ranged from the top 5 to the top 10%, calculated the center of gravity of those selected cells, and measured the mean distance from those cells to their center of gravity. (B) A Lorenz curve representing the proportion of real-to-surrogate ratios assumed by the proportion of cells with the lowest ratios. A Lorenz curve is used to show the degree of “inequality” of a distribution defined by two variables (herein, the real-to-surrogate ratio of individual cells versus the number of cells with a given real-to-surrogate ratio). The inequality can be represented by a Gini coefficient, a number between 0 and 1, where perfect equality has a Gini coefficient of zero, and absolute inequality yields a Gini coefficient of 1. In this case, the Gini coefficient was 0.44. (C) The ratios of the numbers of MSS participations of individual cells in real datasets to the mean values in 20 randomized surrogates are plotted against their mean firing rates. Color-edged circles are 5% significant data points between real datasets and surrogates (paired-t test). The data are the mean ± SDs of 9 videos.

To investigate the spatial organization of MSS-participating neurons, we selected neurons that scored in the top 5% of the frequency ratios and computed the mean distance from these cells to their center of gravity (Figure 10A bottom). We also calculated the mean distances for the top 5–10, 10–15, 15–20, and 20–25% of the neurons and found that they did not differ among these groups [P = 0.43, F_{(4, 40)} = 0.97, One-Way ANOVA, Figure 10A bottom]. Thus, the spatial distribution of neurons is unlikely to depend on the frequency of MSS participation.

We plotted the frequency of MSS participation against the firing rates of the corresponding neurons (Figure 10C). The real-to-surrogate ratios of the MSS participation frequency decreased with the firing rates. These results indicate that frequently firing neurons do not necessarily participate frequently in MSSs, whereas rarely firing neurons seem to contribute more significantly to MSSs.

We then conducted a similar analysis for neuron pairs. Joint participation was defined as when two given neurons simultaneously participated in the same MSS. We counted the number of joint participations for all possible pairs of neurons in each video. Figure 11A shows a representative matrix of the real-to-surrogate ratios of the frequency of joint participation of 96 neurons in a single video. Like the behaviors by single neurons, some neuron pairs co-participated frequently in MSSs, whereas other pairs did not. Moreover, the real-to-surrogate ratios of the joint participation frequency decreased with the joint firing rates of the neuron pairs (Figure 11B); note that the joint firing rates were defined as (f_i × f_j)^1/2, where f_i and f_j are the firing rates of cell_i and cell_j.

FIGURE 11

Figure 11. Relationship between the joint firing rates of neuron pairs and their MSS joint participations. (A) A pseudo-colored matrix representing the numbers of joint participations by two given cells. Data are shown as the ratio of the real rastergrams to their randomized controls. (B) The ratios of the number of MSS co-participations of cell pairs in real datasets to the mean value of 20 randomized surrogates were plotted against their mean joint firing rates. Color-edged circles are 5% significant data points between real datasets and surrogates (paired-t test). The data are the mean ± SDs of 9 videos.

MSS Chains

The significant joint participation in MSSs, which was described above, is of importance in understanding MSS dynamics. Given that single neurons were recruited to different MSSs, the joint participation suggests the existence of “core” neuron groups that were shared with different MSSs. Indeed, we often encountered MSS series in which parts of MSSs were replayed by parts of other MSSs (Figure 12A).

FIGURE 12

Figure 12. Flexibility and complexity of MSSs. (A) An example of an MSS chain. Spikes that consisted of MSSs were partially shared (i.e., relayed) by other MSSs. Arrows indicate the shared spikes, referred to here as core patterns. (B) Two examples of the time evolution of the chains that originated from single MSSs. (C) A rastergram of the onset of individual MSSs, the order of which were arranged along the dendrogram clustered by Ward's method.

We thus defined a MSS chain as a sequence of MSSs that shared at least two spikes. Under this definition, a transition from one MSS to another MSS through the shared spikes was referred to herein as ‘relay’. We also defined a core pattern as shared spikes at each relay step of a MSS chain. Figure 12A shows an example MSS chain in which three MSSs were relayed (from blue to red and from red to green) through two core patterns (indicated by arrows). Of 27,804 MSSs in 9 videos, 22,234 MSSs (80.0%) contained at least one core pattern. Single MSSs contained 1.44 ± 0.85 core patterns, and single core patterns were shared by 4.0 ± 4.2 MSSs (n = 3264 core patterns). The mean firing rate of neurons that were involved core patterns was 0.38 ± 0.26 Hz (n = 401), which was significantly higher than that of non-core neurons [0.08 ± 0.10 Hz, n = 427, P = 4.6 × 10⁻⁸, t₍₈₂₆₎ = 5.57]. Figure 12B shows the whole dynamics of two representative MSS chains; each circle indicates a single MSS, and each line indicates a relay between two MSSs. Single videos included 26.1 ± 11.6 independent MSS chains (n = 9 videos). Single MSS chains consisted of 86.5 ± 284.3 MSSs, contained 53.6 ± 158.1 relay steps, and were 7.9 ± 4.4 s in length (n = 234 chains).

Relay steps did not always reflect simple relays between MSSs. The steps often exhibited divergent relays, in which core patterns in one MSS were subsequently used in two or more MSSs, or convergent relays, in which core patterns in two or more MSSs were simultaneously used in a subsequent MSS (see Figure 12B). Single chains contained 17.8 ± 56.6 divergent relays and 19.8 ± 61.7 convergent relays (n = 234 chains), suggesting that MSS chains constitute MSS subgroups. Therefore, based on core patterns between MSSs, we analyzed MSSs using Ward's method. Figure 12C depicts a dendrogram of MSSs in a representative video, indicating that MSSs were clustered into subgroups. Similar MSS cliqueness was observed in the other videos. As a whole, therefore, MSSs are not mere repetitions of precise firing patterns within specific sets of neurons, but they are parts of larger complex and flexible network dynamics.

Discussion

In this work, we used fMCI with 2 ms temporal resolution and monitored CA3 network activity in cultured hippocampal slices. We searched MSSs using a template-matching method and analyzed them by comparing them to surrogates. We generated the surrogates using a spike-exchanging shuffling procedure, a randomization method that is believed to most reliably minimize false-positive errors because it does not collapse the firing rates of individual neurons or population modulation. We demonstrated that the temporal precision of MSSs in the hippocampus is high at the millisecond level. MSSs consisted of heterogeneous neurons or neuronal subsets. Moreover, we found core patterns across multiple MSSs, and the core patterns served as hubs through which MSSs are replaced with other MSSs and are sometimes split or joined together with other MSSs.

MSS-like fixed firing patterns have been reported in the neocortex of monkeys and rats (Abeles and Gerstein, 1988; Aertsen et al., 1991; Abeles et al., 1993; Prut et al., 1998; Mao et al., 2001; Cossart et al., 2003; Shmiel et al., 2006; Luczak et al., 2007). They have been thought to be experimental evidence supporting the synfire chain hypothesis, a theoretical framework for efficient spike propagation; however, the statistical significance of MSSs has been questioned by several studies (Oram et al., 1999, 2001; Baker and Lemon, 2000; Mokeichev et al., 2007). In this work, we simultaneously monitored spikes of approximately 100 neurons and confirmed that MSSs occurred in spontaneous activity more than chance. Importantly, the existence of MSSs was most significant at a spike jitter of zero frames. Thus, MSSs were repeated with spike precision of less than 2 ms. However, it is statistically difficult to find the true null hypothesis of temporal structures of spikes (Luczak et al., 2007). Therefore, we also tried to demonstrate the existence of MSSs by examining the properties of MSS-participating neurons. If MSSs are a stochastic product, neurons or neuron pairs with higher firing rates would be expected to participate more frequently in MSSs. Conversely, we found that, in the ratio scale, neurons with lower firing rates contributed more to MSSs. This result implies that MSSs are generated by network dynamics independent of the firing rates of individual neurons. Moreover, the number of neurons involved in single MSSs and the MSS length peaked at 6 neurons and 60 ms. These data also suggest that MSSs reflect organized network dynamics.

The majority of MSSs were generated by specific sets of neurons. Similar heterogeneous participations have been suggested by multi-electrode array recordings from cultures of dissociated mouse neocortical neurons (Sun et al., 2010). In this work, we further identified the core patterns that mediated MSS relays. The core patterns seemed to serve as hubs that generated a wide variety of MSSs that were dynamically associated with each other. The clustering of MSSs in the dendrogram suggests that single cells participate in different MSSs and that single MSSs are also involved in larger-scale MSS chains. This “MSS-family” concept is consistent with cell assembly dynamics and also with phase sequences of synfire chains, both of which occur together with network synchronization. However, the complexity of MSSs must be interpreted with caution. For instance, the mean number of neurons involved in single MSSs was 3.4. The mean number of total neurons monitored simultaneously was 91.7, whereas approximately 5000 CA3 neurons exist in a hippocampal slice culture (Kimura et al., 2011). This result means that we missed the vast majority of MSSs. Therefore, we underestimate both the true MSS size and the true complexity of MSS dynamics.

Although MSSs are theoretically accepted to be important for stable spike propagation through neuronal microcircuits that consist of weak and stochastic synapses (Abeles, 1991; Diesmann et al., 1999; Reyes, 2003), the physiological roles of MSSs in brain function remain unclear. Our work did not address this fundamental question, but it is intriguing to find that MSSs occurred during SW-Rs. We believe that MSSs underlie memory replay during SW-Rs. Manipulations of MSSs, i.e., artificial controls to increase or decrease MSSs, will help our understanding of the behavioral function of MSSs.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by Grants-in-Aid for Science Research on Innovative Areas, “Mesoscopic Neurocircuitry” (No. 22115003), from the Ministry of Education, Culture, Sports, Science and Technology of Japan and by the Funding Program for Next Generation World-Leading Researchers (LS023).

References

Abeles, M. (1991). Corticonics: Neural Circuits of the Cerebral Cortex. Cambridge, MA: Cambridge University Press. doi: 10.1017/CBO9780511574566

CrossRef Full Text

Abeles, M., Bergman, H., Margalit, E., and Vaadia, E. (1993). Spatiotemporal firing patterns in the frontal cortex of behaving monkeys. J. Neurophysiol. 70, 1629–1638.

Pubmed Abstract | Pubmed Full Text

Abeles, M., and Gat, I. (2001). Detecting precise firing sequences in experimental data. J. Neurosci. Methods 107, 141–154. doi: 10.1016/S0165-0270(01)00364-8

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Abeles, M., and Gerstein, G. L. (1988). Detecting spatiotemporal firing patterns among simultaneously recorded single neurons. J. Neurophysiol. 60, 909–924.

Pubmed Abstract | Pubmed Full Text

Aertsen, A., Vaadia, E., Abeles, M., Ahissar, E., Bergman, H., Karmon, B., et al. (1991). Neural interactions in the frontal cortex of a behaving monkey: signs of dependence on stimulus context and behavioral state. J. Hirnforsch. 32, 735–743.

Pubmed Abstract | Pubmed Full Text

Baker, S. N., and Lemon, R. N. (2000). Precise spatiotemporal repeating patterns in monkey primary and supplementary motor areas occur at chance levels. J. Neurophysiol. 84, 1770–1780.

Pubmed Abstract | Pubmed Full Text

Buzsaki, G. (2010). Neural syntax: cell assemblies, synapsembles, and readers. Neuron 68, 362–385. doi: 10.1016/j.neuron.2010.09.023

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Buzsaki, G., Leung, L. W., and Vanderwolf, C. H. (1983). Cellular bases of hippocampal EEG in the behaving rat. Brain Res. 287, 139–171.

Pubmed Abstract | Pubmed Full Text

Cossart, R., Aronov, D., and Yuste, R. (2003). Attractor dynamics of network UP states in the neocortex. Nature 423, 283–288. doi: 10.1038/nature01614

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Davidson, T. J., Kloosterman, F., and Wilson, M. A. (2009). Hippocampal replay of extended experience. Neuron 63, 497–507. doi: 10.1016/j.neuron.2009.07.027

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

De Simoni, A., Griesinger, C. B., and Edwards, F. A. (2003). Development of rat CA1 neurones in acute versus organotypic slices: role of experience in synaptic morphology and activity. J. Physiol. 550, 135–147. doi: 10.1113/jphysiol.2003.039099

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Diba, K., and Buzsaki, G. (2007). Forward and reverse hippocampal place-cell sequences during ripples. Nat. Neurosci. 10, 1241–1242. doi: 10.1038/nn1961

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Diesmann, M., Gewaltig, M. O., and Aertsen, A. (1999). Stable propagation of synchronous spiking in cortical neural networks. Nature 402, 529–533. doi: 10.1038/990101

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Foster, D. J., and Wilson, M. A. (2006). Reverse replay of behavioural sequences in hippocampal place cells during the awake state. Nature 440, 680–683. doi: 10.1038/nature04587

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Gahwiler, B. H., Capogna, M., Debanne, D., Mckinney, R. A., and Thompson, S. M. (1997). Organotypic slice cultures: a technique has come of age. Trends Neurosci. 20, 471–477. doi: 10.1016/S0166-2236(97)01122-3

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Gomez-Di Cesare, C. M., Smith, K. L., Rice, F. L., and Swann, J. W. (1997). Axonal remodeling during postnatal maturation of CA3 hippocampal pyramidal neurons. J. Comp. Neurol. 384, 165–180. doi: 10.1002/(SICI)1096-9861(19970728)384:2<165::AID-CNE1>3.0.CO;2-#

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Harris, K. D. (2005). Neural signatures of cell assembly organization. Nat. Rev. Neurosci. 6, 399–407. doi: 10.1038/nrn1669

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Harris, K. D., Csicsvari, J., Hirase, H., Dragoi, G., and Buzsaki, G. (2003). Organization of cell assemblies in the hippocampus. Nature 424, 552–556. doi: 10.1038/nature01834

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hebb, D. O. (1949). The Organization of Behavior: A Neuropsychological Theory. New York, NY: Wiley.

Ikegaya, Y., Aaron, G., Cossart, R., Aronov, D., Lampl, I., Ferster, D., et al. (2004). Synfire chains and cortical songs: temporal modules of cortical activity. Science 304, 559–564. doi: 10.1126/science.1093173

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kendall, M. G., Stuart, A., Ord, J. K., Arnold, S. F., and O'Hagan, A. (1994). Kendall's Advanced Theory of Statistics. New York, NY: Halsted Press.

Kimura, R., Kang, S., Takahashi, N., Usami, A., Matsuki, N., Fukai, T., et al. (2011). Hippocampal polysynaptic computation. J. Neurosci. 31, 13168–13179. doi: 10.1523/JNEUROSCI.1920-11.2011

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Klausberger, T., Magill, P. J., Marton, L. F., Roberts, J. D., Cobden, P. M., Buzsaki, G., et al. (2003). Brain-state- and cell-type-specific firing of hippocampal interneurons in vivo. Nature 421, 844–848. doi: 10.1038/nature01374

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Koyama, R., Muramatsu, R., Sasaki, T., Kimura, R., Ueyama, C., Tamura, M., et al. (2007). A low-cost method for brain slice cultures. J. Pharmacol. Sci. 104, 191–194. doi: 10.1254/jphs.SC0070119

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Laurberg, S., and Zimmer, J. (1981). Lesion-induced sprouting of hippocampal mossy fiber collaterals to the fascia dentata in developing and adult rats. J. Comp. Neurol. 200, 433–459. doi: 10.1002/cne.902000310

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lee, A. K., and Wilson, M. A. (2002). Memory of sequential experience in the hippocampus during slow wave sleep. Neuron 36, 1183–1194. doi: 10.1016/S0896-6273(02)01096-6

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Luczak, A., Bartho, P., Marguet, S. L., Buzsaki, G., and Harris, K. D. (2007). Sequential structure of neocortical spontaneous activity in vivo. Proc. Natl. Acad. Sci. U.S.A. 104, 347–352. doi: 10.1073/pnas.0605643104

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Mao, B. Q., Hamzei-Sichani, F., Aronov, D., Froemke, R. C., and Yuste, R. (2001). Dynamics of spontaneous activity in neocortical slices. Neuron 32, 883–898. doi: 10.1016/S0896-6273(01)00518-9

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Miles, R., and Wong, R. K. (1983). Single neurones can initiate synchronized population discharge in the hippocampus. Nature 306, 371–373. doi: 10.1038/306371a0

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Mizunuma, M., Takahashi, N., Usami, A., Matsuki, N., and Ikegaya, Y. (2009). High-temperature, but not high-pressure, conditions alter neuronal activity. J. Pharmacol. Sci. 110, 117–121. doi: 10.1254/jphs.09031SC

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Mokeichev, A., Okun, M., Barak, O., Katz, Y., Ben-Shahar, O., and Lampl, I. (2007). Stochastic emergence of repeating cortical motifs in spontaneous membrane potential fluctuations in vivo. Neuron 53, 413–425. doi: 10.1016/j.neuron.2007.01.017

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Norimoto, H., Mizunuma, M., Ishikawa, D., Matsuki, N., and Ikegaya, Y. (2012). Muscarinic receptor activation disrupts hippocampal sharp wave-ripples. Brain Res. 1461, 1–9. doi: 10.1016/j.brainres.2012.04.037

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

O'Neill, J., Senior, T., and Csicsvari, J. (2006). Place-selective firing of CA1 pyramidal cells during sharp wave/ripple network patterns in exploratory behavior. Neuron 49, 143–155. doi: 10.1016/j.neuron.2005.10.037

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Oram, M. W., Hatsopoulos, N. G., Richmond, B. J., and Donoghue, J. P. (2001). Excess synchrony in motor cortical neurons provides redundant direction information with that from coarse temporal measures. J. Neurophysiol. 86, 1700–1716.

Pubmed Abstract | Pubmed Full Text

Oram, M. W., Wiener, M. C., Lestienne, R., and Richmond, B. J. (1999). Stochastic nature of precisely timed spike patterns in visual system neuronal responses. J. Neurophysiol. 81, 3021–3033.

Pubmed Abstract | Pubmed Full Text

Pastalkova, E., Itskov, V., Amarasingham, A., and Buzsaki, G. (2008). Internally generated cell assembly sequences in the rat hippocampus. Science 321, 1322–1327. doi: 10.1126/science.1159775

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Prut, Y., Vaadia, E., Bergman, H., Haalman, I., Slovin, H., and Abeles, M. (1998). Spatiotemporal structure of cortical activity: properties and behavioral relevance. J. Neurophysiol. 79, 2857–2874.

Pubmed Abstract | Pubmed Full Text

Reyes, A. D. (2003). Synchrony-dependent propagation of firing rate in iteratively constructed networks in vitro. Nat. Neurosci. 6, 593–599. doi: 10.1038/nn1056

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sasaki, T., Matsuki, N., and Ikegaya, Y. (2007). Metastability of active CA3 networks. J. Neurosci. 27, 517–528. doi: 10.1523/JNEUROSCI.4514-06.2007

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sasaki, T., Minamisawa, G., Takahashi, N., Matsuki, N., and Ikegaya, Y. (2009). Reverse optical trawling for synaptic connections in situ. J. Neurophysiol. 102, 636–643. doi: 10.1152/jn.00012.2009

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sasaki, T., Takahashi, N., Matsuki, N., and Ikegaya, Y. (2008). Fast and accurate detection of action potentials from somatic calcium fluctuations. J. Neurophysiol. 100, 1668–1676. doi: 10.1152/jn.00084.2008

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Shmiel, T., Drori, R., Shmiel, O., Ben-Shaul, Y., Nadasdy, Z., Shemesh, M., et al. (2006). Temporally precise cortical firing patterns are associated with distinct action segments. J. Neurophysiol. 96, 2645–2652. doi: 10.1152/jn.00798.2005

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Smith, K. L., Szarowski, D. H., Turner, J. N., and Swann, J. W. (1995). Diverse neuronal populations mediate local circuit excitation in area CA3 of developing hippocampus. J. Neurophysiol. 74, 650–672.

Pubmed Abstract | Pubmed Full Text

Sun, J. J., Kilb, W., and Luhmann, H. J. (2010). Self-organization of repetitive spike patterns in developing neuronal networks in vitro. Eur. J. Neurosci. 32, 1289–1299. doi: 10.1111/j.1460-9568.2010.07383.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sun, Y., Norimoto, H., Pu, X. P., Matsuki, N., and Ikegaya, Y. (2012). Cannabinoid receptor activation disrupts the internal structure of hippocampal sharp wave-ripple complexes. J. Pharmacol. Sci. 118, 288–294. doi: 10.1254/jphs.11199FP

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Takahashi, N., Oba, S., Yukinawa, N., Ujita, S., Mizunuma, M., Matsuki, N., et al. (2011). High-speed multineuron calcium imaging using Nipkow-type confocal microscopy. Curr. Protoc. Neurosci. 2, 14. doi: 10.1002/0471142301.ns0214s57

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Takahashi, N., Sasaki, T., Matsumoto, W., Matsuki, N., and Ikegaya, Y. (2010). Circuit topology for synchronizing neurons in spontaneously active networks. Proc. Natl. Acad. Sci. U.S.A. 107, 10244–10249. doi: 10.1073/pnas.0914594107

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Takahashi, N., Sasaki, T., Usami, A., Matsuki, N., and Ikegaya, Y. (2007). Watching neuronal circuit dynamics through functional multineuron calcium imaging (fMCI). Neurosci. Res. 58, 219–225. doi: 10.1016/j.neures.2007.03.001

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ujita, S., Mizunuma, M., Matsuki, N., and Ikegaya, Y. (2011). Asynchronously enhanced spiking activity of ischemic neuronal networks. Biol. Pharmaceut. Bull. 34, 764–767. doi: 10.1248/bpb.34.764

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Usami, A., Matsuki, N., and Ikegaya, Y. (2008). Spontaneous plasticity of multineuronal activity patterns in activated hippocampal networks. Neural Plast. 2008:108969. doi: 10.1155/2008/108969

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

West, J. R., and Dewey, S. L. (1984). Mossy fiber sprouting in the fascia dentata after unilateral entorhinal lesions: quantitative analysis using computer-assisted image processing. Neuroscience 13, 377–384. doi: 10.1016/0306-4522(84)90237-9

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text