Spike Synchronization and Rate Modulation Differentially Involved in Motor Cortical Function

See allHide authors and affiliations

Science  12 Dec 1997:
Vol. 278, Issue 5345, pp. 1950-1953
DOI: 10.1126/science.278.5345.1950


It is now commonly accepted that planning and execution of movements are based on distributed processing by neuronal populations in motor cortical areas. It is less clear, though, how these populations organize dynamically to cope with the momentary computational demands. Simultaneously recorded activities of neurons in the primary motor cortex of monkeys during performance of a delayed-pointing task exhibited context-dependent, rapid changes in the patterns of coincident action potentials. Accurate spike synchronization occurred in relation to external events (stimuli, movements) and was commonly accompanied by discharge rate modulations but without precise time locking of the spikes to these external events. Spike synchronization also occurred in relation to purely internal events (stimulus expectancy), where firing rate modulations were distinctly absent. These findings indicate that internally generated synchronization of individual spike discharges may subserve the cortical organization of cognitive motor processes.

The classical view of neural coding has emphasized the importance of information carried by changes in the neurons' spike discharge rates (1). Thus, relating rate changes of individual neurons to cognitive and motor functions has yielded remarkable insights into neuronal mechanisms of sensory (2) and motor (3) processing. In parallel, however, the concept emerged that computational processes in the brain could also rely on the relative timing of spike discharges among neurons within functional groups (4)—commonly called cell assemblies (5). Indeed, it has been argued that the synaptic influence of multiple neurons converging onto others is much stronger if they fire in coincidence (6), making synchrony of firing ideally suited to raising the saliency of responses and to expressing relations among neurons with high temporal precision (4). For this temporal scheme to work, correlations in spike timing must reflect cooperative interactions among the neurons constituting an assembly. Moreover, to permit neurons to participate in different cell assemblies at different times, depending on the present stimulus context and behavioral demands, these correlations should be dynamic. Indeed, such predicted stimulus- and behavior-dependent modulations of spike correlations with differing degrees of temporal precision were observed in various sensory cortical areas, particularly visual (7), auditory (8), somatosensory (9), and frontal (10) areas. However, very little is known about the possible functional role of action potential timing in the motor cortex (11), in spite of the fact that this part of the brain directly governs overt behavior. We present evidence for the hypothesis that precise synchronization of individual action potentials among selected groups of motor cortex neurons is involved in dynamically organizing the cortical network during the planning and execution of voluntary movements. Such synchronization can occur even without modulations of the firing rates, particularly during planning and in relation to internal (“expectation”) events. Thus, our findings call for an extension of prevailing models of motor cortical function that rely exclusively on discharge rates in populations of neurons commonly assessed from averaged responses and recorded from one neuron after another (12). Preliminary results have been presented in abstract form (13).

We designed a delayed-response hand-movement task that allowed us to study neuronal cooperativity in the motor cortex under conditions in which the processing of internal, purely cognitive events could be separated from that of external events, such as presentation of visual stimuli and execution of movements (Fig.1A). During a trial, two successive visual stimuli, separated by a variable time interval [preparatory period (PP)], were presented to the animal. The first stimulus [preparatory signal (PS)] indicated the target position of the requested hand movement, and the second stimulus [response signal (RS)] instructed the monkey to make the required movement and touch the target. In accordance with earlier findings in human studies (14), the reaction time (RT) to onset of movement decreased systematically with increasing duration of the PP (Fig. 1B). We interpret this finding as evidence that the monkey updates its estimate of the probability for the RS to occur. Because this probability attains progressively higher values as the four potential RS-presentation instants (expected signal ES1, ES2, ES3, and RS) pass by, the monkey can adjust its preparatory activities accordingly. During 295 sessions, the activities of two to seven neurons were recorded simultaneously (15) in the primary motor cortex of two monkeys (16). Instances of near-coincident (5 ms) firing of two or more neurons were detected on a trial-by-trial basis by using the unitary event analysis method (Fig.2) (17). Those coincident spike constellations that occurred significantly more often than expected by chance are referred to as unitary events (UEs) (Fig. 2F, pairs of red dots) (18, 19). We interpret the neuronal composition of a UE as defining the subset of recorded neurons that were momentarily engaged in a cell assembly. The occurrence of a considerable number of UEs is revealed in Fig. 2F. Interestingly, these were not distributed randomly in time. Instead, they appeared in clusters, which were loosely (50 to 100 ms) time-locked to behaviorally relevant events, in this case to expected occurrence times of the RS (Fig. 2F, ES2 and ES3). Comparing the auto- and cross-correlograms of the spike data in epochs with significant UEs with those of the intervening epochs did not reveal any conspicuous signs of periodicities, which suggests that these precise spike synchronizations were not linked to the coherent oscillations reported in (11). A total of 359 pairs of neurons were analyzed in this way. In one-third (120 of 359) of them, clusters of precisely coincident (5 ms) UEs were detected. The precise UEs, in turn, were loosely (50 to 100 ms) time-locked to behaviorally relevant events, such as movement events and actually occurring, but also merely expected, stimuli. In a number of cases we even found coincidence precision down to the 1- to 2-ms level, with temporal clustering similar to that at a coincidence precision of 5 ms. This consistent temporal relation between UEs and behavioral events, reflected in the timing and the temporal width of the UE clusters, suggests a functional significance of precisely synchronous action potentials.

Figure 1

(A) Two macaque monkeys were trained to touch a target presented on a video display equipped with a capacitive touch screen. To start a trial, the animal had to hold down a switch. After a constant delay of 1 s (PreD), a PS forming an open circle (diameter, 3 cm) appeared on the screen, indicating the upcoming target position. After a second delay of variable duration (PP), the RS, indicated by filling the circle, instructed the animal to move its hand to the target and touch it. The PP could last 600, 900, 1200, or 1500 ms; these intervals were presented in random order and with equal probability. Thus, from one possible PP duration (and, hence, ES) to the next, the conditional probability for the RS to occur increased, reaching values of 0.25, 0.33, 0.5, and 1. RT and movement time (MT) were recorded in each trial. Criteria for obtaining the reward were (1) to keep the switch pressed during the PP, (2) to move to and touch the target after RS, and (3) to perform the movement so that neither RT nor MT exceeded 1 s. (B) Behavioral results. For each monkey, mean RTs and standard error in relation to PP duration are shown as recorded during the last training sessions (monkey 1, N = 35; monkey 2, N = 59) when they performed the task with more than 90% correct trials. By using a two-factor (PP duration × training session) analysis of variance, changes in RT were found to be highly significant for both monkeys [monkey 1, F(3,102) = 38.19,P < 0.001; monkey 2, F(3,174) = 157.81,P < 0.001].

Figure 2

UE analysis of the spiking activity of two simultaneously recorded single MI neurons. The 36 trials with longest PP duration (1500 ms) were selected in this example. Thus, the monkey could expect the RS to occur at three successive moments (ES1, ES2, and ES3) before it actually occurred at RS. Results of UE analysis over a time interval of 2100 ms, starting 300 ms before PS and lasting until 300 ms after RS, are shown. (A) Conventional raster displays of spike discharges of the two neurons (top, neuron 2; bottom, neuron 3). Each dot represents an action potential, and each horizontal line depicts the spiking activity in a single trial. (B) Spike discharge rate for each neuron (normalized to spikes per second) was computed by sliding a boxcar window of 100 ms in steps of 5 ms over the spike events [dots in (A)]. (C) In a copy of the raster displays from (A), spike coincidences of the two neurons within a precision window of 5 ms are indicated by blue dots. (D) Comparison of measured and expected coincidence rates. The measured coincidence rate (blue curve, normalized to coincidences per second) was derived in the same way as the spike rates, by sliding a boxcar window of 100 ms in steps of 5 ms over the coincident events [blue dots in (C)]. The expected coincidence rate (black curve), based on the null hypothesis of independent firing, was calculated as the product of the individual firing rates [curves in (B)]. (E) For each time window, the statistical significance for a positive difference between measured and expected coincidence rates was calculated from a Poisson distribution (with the mean set to the expected coincidence rate) as the cumulative probability Pof observing the actual number of coincidences or an even larger one by chance. The larger the number of excessive coincidences, the closerP is to 0. Similarly, the larger the number of lacking coincidences, the closer its complement, 1 – P, is to 0, while P approaches 1. To enhance visual resolution at the relevant low-probability values of P or 1 – P, we plotted a logarithmic function of the two: log10[(1 − P)/P]. For excessive coincidences, this function is dominated byP; for lacking coincidences, it is dominated by 1 – P. This procedure, derived from the “surprise” measure (20), is comparable to measuring significance on a decibel scale. It yields positive numbers for excessive coincidences (for example, 1 for P = 0.1 and 2 for P = 0.01) and negative numbers for lacking coincidences, and it changes sign at the chance level (P = 0.5). Whenever the significance value of an excess number of observed coincidences exceeded a fixed threshold, α+ (hereP = 0.05), this defined an epoch with significantly more coincidences than expected by chance. We marked these significant precise coincidences as UEs, and they are indicated by red dots in the raster displays in (F). Occasionally, the significance value dropped below the negative threshold level α(here, shortly after PS). Such cases, defining epochs with significantly fewer coincidences than expected by chance, were more rare because of the generally low expected coincidence rates in these data (21) and are not considered further here. The widths of the coincidence window and of the sliding time window as well as the significance threshold α+ could be set to a selected (fixed) value in the analysis; values of 5 ms and 100 ms, respectively, and P = 0.05 were used throughout the results reported here.

If synchronization of neuronal activity is relevant in task-specific computations, not only the occurrence of synchronization but also its strength and composition should be correlated with behavioral aspects such as stimulus probability and RT, motivation and performance, and expectancy. Two examples demonstrate that such relationships were indeed observed (Fig.3). In the first example (Fig. 3A), we consider UEs occurring in relation to signal expectancy. Two neurons are shown that synchronized their activities each time about 150 to 50 ms before the RS was expected (ES1 to ES3, RS) (top panel). It is clear from the increasing numbers of UEs in successive clusters that the rate of coincident firing grew as the duration of the PP increased and, hence, the probability for the RS to occur increased. This phenomenon, which was observed in both monkeys, suggests a positive correlation between the growing stimulus expectancy and the occurrence of UEs. A similar correlation was observed when trials for a given PP duration were subdivided into two groups, one characterized by short (top) and the other by long (bottom) RTs. In line with the prediction from the behavioral RT results (Fig. 1B), many more coincidences appeared in the short than in the long RT group. The second example (Fig. 3B) shows a relation between synchronized spike activity and the motivation of the animal. During the trials summarized in the top panel, the monkey performed the task correctly and efficiently; during the trials shown in the bottom panel (collected at the end of the experimental session), the animal had stopped responding to the RS but still initiated the trials, so that stimuli were presented as usual. Observe that epochs of excessive synchronized spike events occurred systematically in relation to the PS during the correctly performed trials (top), whereas such synchronization is completely lacking in the nonresponding trials (bottom). Note that in the latter condition the firing rates had also dropped, but this reduction was taken into account by the UE-analysis technique (17). Finally, it is also shown that the composition of the set of the jointly firing neurons may systematically alternate between distinct constellations during successive stages of the trial (Fig. 3B). Before PS, neurons 4 and 6 synchronized their spike events; immediately after PS (from 0 to 250 ms), neurons 2 and 4 synchronized; then (from 250 to 350 ms), neurons 4 and 6 synchronized again; and, finally, all three neurons (2, 4, and 6) coordinated their spiking activity. Such alternating neuron composition suggests a sequential engagement of three different cell assemblies (or cycling through a single, composite cell assembly). In 20 of 43 sessions (46%), the analysis of three or more neurons in parallel revealed such rapid changes in the constellation of synchronously firing cells in systematic relation to behaviorally relevant events.

Figure 3

Synchronized neuronal activity in relation to behavioral aspects of the task. (A) UEs (indicated by red dots in the raster displays) in relation to both expectancy and RT during selected trials in which the longest PP was presented. Trials were divided in three approximately equally sized groups: trials with shortest RTs (top; mean = 270 ms, N = 25), intermediate RTs (not shown), and longest RTs (bottom; mean = 340 ms, N = 19). (B) UEs in relation to performance of the animal. Trials were divided into two groups: trials during which the animal performed the task correctly (top; N = 107), and trials that were correctly initiated by the animal and, hence, the visual stimuli (RS and PS) occurred in the usual way, but the animal did not respond to the RS (bottom; N = 94). Further details are given in Fig. 2.

It is conceivable that synchronization of spiking activity and modulation of discharge rates represent, at least partly, independent computational strategies used by the brain. We tested this hypothesis for all 120 pairs of neurons in which UEs occurred in temporal relation to behaviorally relevant events. We found that spike synchronization could be observed both with and without simultaneous phasic changes in neuronal firing rates (Fig. 4A). Because both PS and RS were predictable, synchronization of spiking activity and changes in discharge rate were both analyzed in relation to the expectancy as well as to the actual occurrence of the PS and RS. The results are summarized (Fig. 4, B and C). Several observations can be made. First, precise spike synchronization was found in relation to both expected and actually occurring stimuli, albeit more often with the latter (70 versus 110). However, the relations between spike synchronization and rate changes were very different in the two cases. If the stimulus was expected but did not occur, spike synchronization was never (PS, 0 of 28) or rarely (RS, 6 of 42) accompanied by a change in discharge frequency of any of the participating neurons. By contrast, when the stimulus did occur, excess synchronized spiking was often (PS, 19 of 45) or almost always (RS, 61 of 65) accompanied by a modulation of discharge rate. We wish to emphasize, however, that our procedure to identify UEs takes changes of firing rate into consideration (Fig. 2) (17). Hence, the observed UEs are not simply a consequence of rate changes. Instead, these findings indicate that the association between spike synchronization and rate modulation is tightly connected to the actual occurrence of the stimulus. Thus, both internal and external cognitive events are associated with precisely coordinated spike activity, but only external events are also accompanied by discharge rate variation.

Figure 4

Relation between synchronization of spike events and modulation of discharge rate. (A) Discharge patterns of three neurons recorded simultaneously during 96 trials. UEs are indicated by red dots in the raster diagrams. Data from three types of trials (PP durations, 900, 1200, and 1500 ms) were pooled. The visualization was interrupted 800 ms after PS and continued for the last 600 ms of the trial (ES1, first expected RS; MO, movement onset; ME, movement end) (further details are presented in Fig. 2). (B) Distribution of activity modulation indices during epochs with UEs computed for that neuron in the pair that exhibited the strongest rate modulation. Phasic changes in neuronal discharge rate were detected by quantitative analysis of the peristimulus-time histograms in fine-temporal resolution (bin width, 5 ms). The activity modulation index, a quantitative measure adopted from the visual contrast measure, was calculated for each neuron as follows. During epochs with UEs, the minimum (Min) and maximum (Max) values of the discharge rate in each 5-ms bin and the mean (Mean) discharge rate over the entire UE epoch were determined. The activity modulation index was then calculated as (Max–Min)/Mean. Values exceeding 0.6 (for Poisson spike rates of 10 to 50 spikes per second, this corresponds to rate fluctuations beyond 1 to 2 standard deviations from the mean) were considered to reflect significant modulations of discharge rate (yellow bar); for values below 0.6 (green bar) such changes were considered absent. Four distinct epochs during the trial were considered: when either PS or RS was expected and when either PS or RS actually occurred. (C) Numbers of pairs of neurons with significant UEs in relation to the visual stimuli (PS and RS) either when these were expected (left) or when they actually occurred (right). Numbers of neuron pairs in which at least one of the neurons changed its discharge rate during UE epochs, as judged from the activity modulation index, are shown in yellow; numbers of neuron pairs in which none of the neurons exhibited such a simultaneous firing rate modulation are shown in green.

Taken together, our findings demonstrate the existence of precise (5 ms) synchronization of individual spike discharges among selected groups of neurons in the motor cortex. This synchronization is associated with distinct phases in the planning and execution of voluntary movements, indicating that it plays a functional role. Moreover, our findings suggest that, under behavioral conditions as studied here, the brain may use different strategies in different contextual situations. To process a purely cognitive (that is, an internal and behaviorally relevant) event, neurons preferentially synchronize their spike occurrences without changing their firing rates at the same time. By contrast, when processing an external, behaviorally relevant event, neurons tend to synchronize their spike occurrences and modulate their firing rates at the same time. Thus, precise synchronization of spike events and modulation of discharge rate may serve different and complementary functions; they act in conjunction at some times but not at others, depending on the behavioral context.

  • * To whom correspondence should be addressed. E-mail: ariehle{at}


View Abstract

Navigate This Article