|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1The Interdisciplinary Center for Neural Computation; 2Department of Physiology, Hadassah Medical School; 3The School of Engineering and Computer Science; and 4Department of Neurobiology, Life Science Institute, Faculty of Sciences, The Hebrew University, Jerusalem, Israel
Submitted 25 March 2008; accepted in final form 6 October 2008
|
|
ABSTRACT |
|---|
|
|
|
INTRODUCTION |
|---|
|
Most BG models use firing rate for information encoding (Albin et al. 1989
; Bar-Gad et al. 2003
; Berns and Sejnowski 1998
; DeLong 1990
; Wickens 1997
). These models rely on physiological evidence that BG neurons, as in most areas of the nervous system, change their firing rate in response to behaviorally and emotionally relevant events. However, visual examination of our data suggested that response duration is also modulated by behavioral events (Fig. 1; also see additional example in Fig. 3G of Morris et al. 2004
). This motivated us to study the response duration as an encoding strategy of BG neurons that does not necessarily rely on firing rate alone. Response duration is an attractive encoding candidate since it naturally emerges from models of the basal ganglia that emphasize the role of information gating in BG processing (Deniau and Chevalier 1985
; Hikosaka and Wurtz 1983
) and the recent suggestion of duration encoding of reward omission by SNc neurons (Bayer et al. 2007
). Unlike rate models, gating models emphasize not only the extent of the gate opening (firing rate amplitude) but also the duration of the opening (e.g., response duration). The required processing of response features, such as response duration, may result in a decrease in information. [According to the information processing inequality, processing a response feature can decrease only the amount of information that could be provided by the original feature (Cover and Thomas 1991
).] However, in gating models, the next neuronal station will not be required to estimate the duration, but to make use of the information passed/blocked through the gate. We therefore used mutual information (MI) to quantify the modulation of response duration and to compare it quantitatively with other carriers of behavioral information in the basal ganglia. The MI is a nonparametric measure of the association between two variables: in our case, the visual stimuli (with four different future reward conditions) and the neural responses. We used the MI as a measure of the association, since the MI is able to detect statistical dependence beyond the linear relationships that are detected by the correlation coefficient. When MI = 0, the two variables are statistically independent and the neuronal activity does not encode information regarding the cues. When the MI is equal to the entropy of the stimulus set, in our case about 2 bits (depending on the actual frequency of the four cues, which varied somewhat from session to session), the stimuli can be perfectly distinguished by the response on a single-trial basis. Intermediate values correspond to intermediate resolvability of the stimuli by observation of the neural responses on a single-trial basis (Cover and Thomas 1991
).
|
|
|
METHODS |
|---|
|
10 correct trials of each of the possible trial types. The behavioral paradigm, recording methods, and data preprocessing are detailed in Arkadir et al. (2004)
We limited our study to the first 500 ms of the cue epoch to avoid the confounding effects of the limb and the reward licking movements. To study the response on a trial-by-trial basis, we divided the spike trains into 50-ms bins aligned on the onset of the cue signal indicating future reward prospects. We then characterized the type of activity in each bin of each trial as either no response, decrease in firing rate, or increase in firing rate by comparing the number of spikes in that bin to baseline activity. The baseline activity was calculated as the mean spike count (firing rate) during intertrial intervals. In >43% of the responses, we found that the peak response deviated by >2SDs. However, the changes in firing rates could be gradual and, at the response onsets and offsets, the changes in firing rates were often substantially smaller. Therefore decreases/increases in firing rate were detected by changes of
1.25SD with respect to baseline activity. We repeated the same analysis with a criterion at 1.5SD, with no appreciable difference (data not shown). For each trial, we defined the duration of the negative/positive responses as the number of bins in the longest sequence during the time period between 150 and 500 ms after cue onset (allowing single or multiple gaps of no more than one bin) in which the activity was characterized as either a negative or positive response. We denoted this time period as the "firing rate response period." Figure 1 depicts this process in five different neurons from the three BG structures.
We estimated the MI between the cues and the neuronal responses using explicit bias-information loss trade-off optimization (Nelken et al. 2005
). Four response features were used: 1) response duration (as defined earlier), 2) the response spike count (the cumulative spike count during the response period), 3) the extreme response amplitude (maximal for positive response and minimal for negative response), and 4) the total spike count (the cumulative spike count during the entire tested time period; i.e., the period between 150 and 500 ms after cue onset). Although calculating the total spike count over a shorter epoch could have produced more information in some cases it might also have induced a bias between the different events and the different neural populations. We therefore preferred the conservative choice of a fixed window that enabled better comparison of the different populations. The 150- to 500-ms time window was chosen since it captures the response in the tested data.
In the SNr and GPe, but not in the SNc, neurons responded to the future direction of movement indicated by the cue as well as to the reward probability (Arkadir et al. 2004
; DeLong et al. 1983
; Wichmann and Kliem 2004
). In these structures the information on the reward probability was calculated for each of the movement directions separately and averaged according to the number of times each direction occurred. Cells were scored as response duration-informative if the information extracted from a cell's positive or negative response duration on future reward probability was significant when compared with 150 surrogate shuffled MI values. (The reward probabilities were shuffled between trials and the significance of the nonshuffled vs. shuffled results was calculated according to the Mann–Whitney test, P < 0.05.) Cells in which <20% of the trials from each type had a positive/negative response were not considered as responding by positive/negative responses. For each nucleus, we calculated the mean response duration for duration-informative neurons as a function of reward probability.
To check whether the information provided by the other response features was redundant or synergistic with the information provided by response duration and to evaluate the exclusiveness of the information to one of the different response features (duration, extreme responses, or spike count) we used a normalized synergy measure. In synergistic encoding the information provided by two response features simultaneously is greater than the sum of the information provided by each one of them separately, whereas in redundant encoding the sum of information is smaller than the information provided by both features together. The normalized synergy between the information provided by a response feature R1 on a stimulus S and the information provided by a response feature R2 on S is given by (Schneidman et al. 2003
)
![]() |
0. It ranges from 1 [total synergy, MI(S; R1) = MI(S; R2) = 0]; through 0 [MI independent variables, when MI(S; R1) + MI(S; R2) = MI(S; R1, R2)]; to –1 [total redundancy, when MI(S; R1) = MI(S; R2) = MI(S; R1, R2)].
The synergy measurement does not indicate how the information is distributed between the two variables. For example, norm syn(R1, R2) = –0.5 can occur in a number of scenarios: 1) both variables have equal amount of information on S and they share two thirds of their information [MI(S; R1) = MI(S; R2) = 0.75MI(S; R1, R2)]. 2) R1 has twice as much information on S than R2 and all the information R2 has is already provided by R1 [MI(S; R1) = MI(S; R1, R2) = 2MI(S; R2)]. In a similar way, ambiguity exists when norm syn(R1, R2) = 0 (Fig. 2). We therefore also used a different measure: the normalized contribution. The normalized contribution of R2 to the information on S carried by the pair R1 and R2 is the fraction of information that could not have been extracted without R2
![]() |
|
|
|
RESULTS |
|---|
|
Most neurons in the SNr, SNc, and GPe showed significant spike rate responses (Mann–Whitney test on spike counts of response bins vs. intertrial-interval bins, P < 0.05) to the appearance of at least one of the four possible cues (Table 1, % responding neurons). Hereafter, we refer to the significant decreases in firing rate as "negative responses" and the increases in firing rate as "positive responses."
|
All of the responding SNc neurons exhibited positive responses (increase in discharge rate) to the cue and thus all of the duration-informative neurons in the SNc showed a positive response. However, twice as many of the SNr duration-informative neurons (68%) responded with a negative rather than a positive response, whereas in the GPe the duration-informative neurons distributed evenly between positive and negative responses.
The average response duration in SNr and SNc duration-informative neurons exhibited monotonic increases/decreases with the probability of future reward (Fig. 3). Such behaviors were not observed in the average response duration of GPe duration-informative neurons. The mean duration of the SNc positive responses was longer for high reward probabilities and shorter for low reward probabilities (one-way ANOVA, P < 0.01). The mean response duration in SNr neurons presented a different behavior from the SNc responses. Whereas the SNr negative response followed the same monotonic increase as SNc neurons (duration increased for high reward probability, one-way ANOVA, P < 0.01), the positive response tended (without reaching statistical significance) to be shortened with increased reward probability. In addition, the average durations of the negative GPe responses in duration-informative neurons are significantly shorter (one-way ANOVA, P < 0.01) than the positive responses. We checked whether the negative versus positive difference seen in response duration of SNr neurons also exists in other response features discussed herein. Only the mean total spike count of negative duration-informative neurons showed a significant monotonic decrease (one-way ANOVA, P < 0.01) displaying higher values for low reward probabilities.
|
|
|
Extracting response duration from single trials is rather complex (Fig. 1), since it involves classifying each time bin in the responses of each single trial. Different ways of decoding the neuronal responses may be more informative. The standard way of quantifying responses is by counting the total spike count in a fixed, long window. Whereas total spike count encompasses time bins that may not be really important, possibly diluting the information it carries through the addition of noise, it may also take into account the responses during time bins that were wrongly classified as nonresponsive. In practice, the second effect probably dominated since we found that total spike count carried as much information as, or more information than, the response duration. This can be seen in Fig. 4 where the MI extracted from the total spike count over the fixed 350-ms window for both positive and negative responses was typically higher than the MI carried by response duration (as also reflected in Table 2). Nevertheless, since the duration did not add significantly to the information provided by the total spike count, the information provided by both the total spike count and the duration is actually provided by the total spike count alone. Therefore the amount of information given by both features exclusive to the total spike count (see "norm contribution" in Table 2: 32%, range 20–50%) is the amount of information given by the total spike count and not by the response duration, i.e., most (68%, range 50–80%) of the information provided by the total spike count was encoded by the response duration as well. Thus the difference between response duration and total spike count as a carrier of information was not substantial. In addition, we tested more detailed temporal representations of the response (such as analyzing the activity as a sequence of responsive and nonresponsive bins). However, the bias-corrected amount of information did not significantly increase with greater detail, suggesting that there is no easily decoded additional information in the temporal pattern of the 50-ms bins of the BG responses (data not shown).
|
|
DISCUSSION |
|---|
|
Evidence for the importance of response duration in the nervous system (Christensen-Dalsgaard et al. 1998
; DeBusk et al. 1997
; Kanold and Manis 2005
; Nagata et al. 2003
; Rogers and Newland 2002
; Zhang et al. 2004
) and even specifically in the responses of striatal cholinergic neurons (Aosaki et al. 1995
; Ravel et al. 2003
) and in responses of pallidal neurons (Anderson and Turner 1991
) have been reported. Recently, Steuber et al. (2007)
found that the best criterion to distinguish between different patterns of the cerebellar parallel fibers was the duration of the pause. O'Donnell and Grace (1998)
showed that a known psychotic modulator (phencyclidine) acts by reducing the frequency and duration of the spontaneously occurring depolarized plateaus observed in the membrane potential of accumbens neurons. In addition, Kepecs and Raghavachari (2007)
presented a model that can account for the generation of different duration up-state activity in striatal neurons. By simulating a bursting neuron model, Kepecs and colleagues showed that bursts of different durations can code for different stimulus features (Kepecs and Lisman 2003
; Kepecs et al. 2002
). Moreover, they found that synapses can be tuned to preferentially respond to specific burst durations and demonstrated the decodability of a neural code based on burst duration (Kepecs and Lisman 2004
). In addition, a study of dopaminergic neurons suggested that these neurons overcome the limited dynamic range available for encoding of negative values by modification of the duration of their negative responses (Bayer et al. 2007
). Finally, a recent study (Person and Perkel 2005
) revealed in zebra finch songbirds that high-frequency trains of pallidal spikes can drive activity in thalamic relay neurons by rebound excitation. The latency of this rebound is strongly affected by the duration of the pallidal spike train (their Fig. 4). Thus positive SNr responses will yield different thalamic rebound excitation as a function of the response duration and, perhaps, longer SNr pauses will cause longer and probably stronger thalamic disinhibitory effects.
Our results are in line with models of the basal ganglia in which information transmission is carried through gating mechanisms, where the neuronal activity determines the extent of gate aperture and its duration (Deniau and Chevalier 1985
; Mink 1996
). GPe is probably involved in gating the inputs to the subthalamic nucleus, striatum, and the output structures of the BG, whereas the SNr gates the reciprocal thalamocortical neuronal loops. Contrary to other models (Mink 1996
), gating models suggest that the basal ganglia output enables, rather than selects or initiates (McHaffie et al. 2005
; Mink 1996
), movements or other voluntary behaviors. We suggest that the duration of the gate aperture might take part in the competition between possible actions (giving advantage to those BG–thalamic channels with longer opened gate duration). In addition, reward expectation has been shown to influence response latency (Lauwereyns et al. 2002
); it is thus possible that response duration is part of the mechanism that modulates or affects the action's vigor through the multiple channels of the basal ganglia output to the reciprocal thalamocortical networks. The significant fraction of duration-informative SNr neurons with negative responses is compatible with this suggested role of the SNr. The SNr has GABAergic projections to the thalamus and to the superior colliculus (Hikosaka and Wurtz 1983
; Redgrave et al. 1992
); thus when SNr neurons decrease their activity, they disinhibit their targets. This may be the mechanistic implementation of the gating mechanism, where the background tonic activity and positive responses of SNr neurons close the gate, whereas negative SNr responses open the gate (Deniau and Chevalier 1985
; Hikosaka and Wurtz 1983
). Our finding that the mean negative response durations of duration-informative neurons monotonically increases (long durations for high probability and short durations for low probability; Fig. 3) is also compatible with the gating model. Negative SNr responses open the gate and we expect the gate to be opened for longer durations (and enable a movement) for the higher reward prospects. The positive responses show (nonsignificant) an opposite trend. This could also be explained by the gating model. Since closing the gate does not determine the movement, its duration is less significant. The monotonic decrease seen in the mean total spike count of negative duration-informative neurons—displaying higher values for low reward probabilities—is in line with the monotonic increase found in the mean response duration of those neurons. The total spike count consists of responsive and nonresponsive bins; thus long negative responses that have fewer bins with normal activity and more bins with low activity will have a low total spike count and vice versa. In addition, it has been shown that dopamine (D1 receptors) inhibits SNr activity (Kliem et al. 2007
). Therefore in line with SNr negative responses, the mean response duration in SNc (positive) responses is long for high reward probabilities and short for low probabilities.
A negative response (i.e., response with a decrease in firing rate) has the advantage of being detectable faster than an increase of the same magnitude (APPENDIX). However, in negative responses, when the neurons decrease their firing rate to close to zero for different time lengths, the minimal firing rate and the response spike count do not differ between different probabilities and thus provide little information on the reward probability. In this situation, the response duration remains informative (consequently, the total spike count that includes time bins in which short responses return to normal activity will also be informative). On the other hand, in positive responses, when the neurons respond with an increase in firing rate for different time lengths, the response spike count varies between reward probabilities. Hence for neurons with positive responses, duration did not have a major advantage over maximal firing rate or response spike count.
Spike train analyses usually have focused on two extreme schemes of neuronal encoding. Most studies have dealt with the highly information reduced measure of spike count or rate (Ahissar et al. 2000
; Georgopoulos et al. 1986
; Gershon et al. 1998
; Shadlen and Newsome 1998
). Other studies have focused on information in spike patterns (Arabzadeh et al. 2006
; Friedrich et al. 2004
; Ikegaya et al. 2004
; Jones et al. 2004
; Panzeri et al. 2001
; Prut et al. 1998
; Victor and Purpura 1996
). However, even when spike timing carries information, the use of simple reduced measures of the timing of the responses may be sufficient (Nelken et al. 2005
). Our observations demonstrate that coding schemes based on response duration are viable and that they use much of the information available in the responses of BG neurons. Furthermore, we show that response duration encodes reward probability not only in positive (increase in firing rate, bursts) responses, but also (in the GABAeregic BG nuclei) in negative (decrease in firing rate) responses, and that duration encoding is not a simple epiphenomenon of rate encoding. Finally, we argued that response duration may be highly relevant to the effects of BG neurons on their target structures. The significant amount of information encoded in response duration makes response duration a prominent coding element that can mechanistically control the influence of the basal ganglia on its target structures.
|
|
APPENDIX |
|---|
|
The sequential probability ratio test (SPRT; Wald 1947
) is a serial Bayesian test designed to determine which of two hypothesized distributions is used to generate data samples. In this test we obtain data samples and examine the log likelihood (the probability of a given distribution to generate these data points) ratio of the two distributions. The test sequentially inputs data samples until the ratio is above an acceptance or below a rejection threshold. At that time the test can determine (with a reliability that depends solely on the thresholds) which of the two distributions generated the data.
We assume that the spike trains are generated by a Poisson distribution and that an external/internal event can cause the neuron to switch from one firing rate (
1) to a different firing rate (
2). We seek to evaluate how many samples (spike counts/rates) are needed (using SPRT) to determine the neuron's firing rate after the switch has occurred. The likelihood ratio for the two hypotheses is
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
) denote the number of samples needed to detect a distribution switch from
= R to
= R(1 +
), nR||R(1–
) denote the number of samples needed to detect a distribution switch from
= R to
= R(1 –
), and Sth denote the ratio acceptance threshold (the threshold to detect the distribution change). We can express the required sample sizes as
![]() |
![]() |
![]() |
![]() |
= R to
= R(1 +
) and the number of samples needed to detect a distribution switch from
= R to
= R(1 –
) is always >1 and reaches 2.6 as
reaches 1. Thus detecting a decrease in firing rate will require fewer data samples—i.e., less time than detecting an increase of the same magnitude. |
|
GRANTS |
|---|
|
|
|
ACKNOWLEDGMENTS |
|---|
|
|
|
FOOTNOTES |
|---|
Address for reprint requests and other correspondence: N. Parush, Department of Physiology, The Hebrew University–Hadassah Medical School, P.O. Box 12272, Jerusalem, Israel 91120 (E-mail: naamap{at}alice.nc.huji.ac.il)
|
|
REFERENCES |
|---|
|
Albin RL, Young AB, Penney JB. The functional anatomy of basal ganglia disorders. Trends Neurosci 12: 366–375, 1989.[CrossRef][Web of Science][Medline]
Anderson ME, Turner RS. A quantitative analysis of pallidal discharge during targeted reaching movement in the monkey. Exp Brain Res 86: 623–632, 1991.[Web of Science][Medline]
Aosaki T, Kimura M, Graybiel AM. Temporal and spatial characteristics of tonically active neurons of the primate's striatum. J Neurophysiol 73: 1234–1252, 1995.
Arabzadeh E, Panzeri S, Diamond ME. Deciphering the spike train of a sensory neuron: counts and temporal patterns in the rat whisker pathway. J Neurosci 26: 9216–9226, 2006.
Arkadir D, Morris G, Vaadia E, Bergman H. Independent coding of movement direction and reward prediction by single pallidal neurons. J Neurosci 24: 10047–10056, 2004.
Averbeck BB, Latham PE, Pouget A. Neural correlations, population coding and computation. Nat Rev Neurosci 7: 358–366, 2006.[CrossRef][Web of Science][Medline]
Bar-Gad I, Morris G, Bergman H. Information processing, dimensionality reduction and reinforcement learning in the basal ganglia. Prog Neurobiol 71: 439–473, 2003.[CrossRef][Web of Science][Medline]
Barlow HB. Sensory mechanisms, the reduction of redundancy, and intelligence. In: NPL Symposium on the Mechanization of Thought Processes, No. 10. London: Her Majesty's Stationery Office, 1959, p. 535–539.
Bayer HM, Lau B, Glimcher PW. Statistics of midbrain dopamine neuron spike trains in the awake primate. J Neurophysiol 98: 1428–1439, 2007.
Berns GS, Sejnowski TJ. A computational model of how the basal ganglia produce sequences. J Cogn Neurosci 10: 108–121, 1998.[CrossRef][Web of Science][Medline]
Bialek W, Rieke F, de Ruyter van Steveninck RR, Warland D. Reading a neural code. Science 252: 1854–1857, 1991.
Bolam JP, Hanley JJ, Booth PA, Bevan MD. Synaptic organisation of the basal ganglia. J Anat 196: 527–542, 2000.[CrossRef][Web of Science][Medline]
Chechik G, Anderson MJ, Bar-Yosef O, Young ED, Tishby N, Nelken I. Reduction of information redundancy in the ascending auditory pathway. Neuron 51: 359–368, 2006.[CrossRef][Web of Science][Medline]
Christensen-Dalsgaard J, Jorgensen MB, Kanneworff M. Basic response characteristics of auditory nerve fibers in the grassfrog (Rana temporaria). Hear Res 119: 155–163, 1998.[CrossRef][Web of Science][Medline]
Cover TM, Thomas JA. Elements of Information Theory. New York: Wiley, 1991.
DeBusk BC, Debruyn EJ, Snider RK, Kabara JF, Bonds AB. Stimulus-dependent modulation of spike burst length in cat striate cortical cells. J Neurophysiol 78: 199–213, 1997.
DeLong MR. Primate models of movement disorders of basal ganglia origin. Trends Neurosci 13: 281–285, 1990.[CrossRef][Web of Science][Medline]
DeLong MR, Crutcher MD, Georgopoulos AP. Relations between movement and single cell discharge in the substantia nigra of the behaving monkey. J Neurosci 3: 1599–1606, 1983.[Web of Science][Medline]
Deniau JM, Chevalier G. Disinhibition as a basic process in the expression of striatal functions. II. The striato-nigral influence on thalamocortical cells of the ventromedial thalamic nucleus. Brain Res 334: 227–233, 1985.[CrossRef][Web of Science][Medline]
Friedrich RW, Habermann CJ, Laurent G. Multiplexing using synchrony in the zebrafish olfactory bulb. Nat Neurosci 7: 862–871, 2004.[CrossRef][Web of Science][Medline]
Georgopoulos AP, Schwartz AB, Kettner RE. Neuronal population coding of movement direction. Science 233: 1416–1419, 1986.
Gershon ED, Wiener MC, Latham PE, Richmond BJ. Coding strategies in monkey V1 and inferior temporal cortices. J Neurophysiol 79: 1135–1144, 1998.
Hikosaka O, Wurtz RH. Visual and oculomotor functions of monkey substantia nigra pars reticulata. IV. Relation of substantia nigra to superior colliculus. J Neurophysiol 49: 1285–1301, 1983.
Ikegaya Y, Aaron G, Cossart R, Aronov D, Lampl I, Ferster D, Yuste R. Synfire chains and cortical songs: temporal modules of cortical activity. Science 304: 559–564, 2004.
Jones LM, Depireux DA, Simons DJ, Keller A. Robust temporal coding in the trigeminal system. Science 304: 1986–1989, 2004.
Kanold PO, Manis PB. Encoding the timing of inhibitory inputs. J Neurophysiol 93: 2887–2897, 2005.
Kepecs A, Lisman J. Information encoding and computation with spikes and bursts. Network 14: 103–118, 2003.[Web of Science][Medline]
Kepecs A, Lisman J. How to read a burst duration code. Neurocomputing 58: 1–6, 2004.[CrossRef][Web of Science]
Kepecs A, Raghavachari S. Gating information by two-state membrane potential fluctuations. J Neurophysiol 97: 3015–3023, 2007.
Kepecs A, Wang X-J, Lisman J. Bursting neurons signal input slope. J Neurosci 22: 9053–9062, 2002.
Kliem MA, Maidment NT, Ackerson LC, Chen S, Smith Y, Wichmann T. Activation of nigral and pallidal dopamine D1-like receptors modulates basal ganglia outflow in monkeys. J Neurophysiol 98: 1489–1500, 2007.
Las L, Stern EA, Nelken I. Representation of tone in fluctuating maskers in the ascending auditory system. J Neurosci 25: 1503–1513, 2005.
Lauwereyns J, Watanabe K, Coe B, Hikosaka O. A neural correlate of response bias in monkey caudate nucleus. Nature 418: 413–417, 2002.[CrossRef][Web of Science][Medline]
McHaffie JG, Stanford TR, Stein BE, Coizet V, Redgrave P. Subcortical loops through the basal ganglia. Trends Neurosci 28: 401–407, 2005.[CrossRef][Web of Science][Medline]
Mink JW. The basal ganglia: focused selection and inhibition of competing motor programs. Prog Neurobiol 50: 381–425, 1996.[CrossRef][Web of Science][Medline]
Morris G, Arkadir D, Nevet A, Vaadia E, Bergman H. Coincident but distinct messages of midbrain dopamine and striatal tonically active neurons. Neuron 43: 133–143, 2004.[CrossRef][Web of Science][Medline]
Nagata T, Suzuki H, Zhang R, Ozaki M, Kawakami Y. Mechanical stimulation activates small fiber mediated nociceptive responses in the nucleus gigantocellularis. Exp Brain Res 149: 505–511, 2003.[Web of Science][Medline]
Nelken I, Chechik G, Mrsic-Flogel TD, King AJ, Schnupp JW. Encoding stimulus information by spike numbers and mean response time in primary auditory cortex. J Comput Neurosci 19: 199–221, 2005.[CrossRef][Web of Science][Medline]
Nevet A, Morris G, Saban G, Fainstein N, Bergman H. Discharge rate of substantia nigra pars reticulata neurons is reduced in non-parkinsonian monkeys with apomorphine-induced orofacial dyskinesia. J Neurophysiol 92: 1973–1981, 2004.
O'Donnell P, Grace AA. Phencyclidine interferes with the hippocampal gating of nucleus accumbens neuronal activity in vivo. Neuroscience 87: 823–830, 1998.[CrossRef][Web of Science][Medline]
Panzeri S, Petersen RS, Schultz SR, Lebedev M, Diamond ME. The role of spike timing in the coding of stimulus location in rat somatosensory cortex. Neuron 29: 769–777, 2001.[CrossRef][Web of Science][Medline]
Person AL, Perkel DJ. Unitary IPSPs drive precise thalamic spiking in a circuit required for learning. Neuron 46: 129–140, 2005.[CrossRef][Web of Science][Medline]
Prut Y, Vaadia E, Bergman H, Haalman I, Slovin H, Abeles M. Spatiotemporal structure of cortical activity: properties and behavioral relevance. J Neurophysiol 79: 2857–2874, 1998.
Ravel S, Legallet E, Apicella P. Responses of tonically active neurons in the monkey striatum discriminate between motivationally opposing stimuli. J Neurosci 23: 8489–8497, 2003.
Redgrave P, Marrow L, Dean P. Topographical organization of the nigrotectal projection in rat: evidence for segregated channels. Neuroscience 50: 571–595, 1992.[CrossRef][Web of Science][Medline]
Rogers SM, Newland PL. Gustatory processing in thoracic local circuits of locusts. J Neurosci 22: 8324–8333, 2002.
Schneidman E, Bialek W, Berry MJ. Synergy, redundancy, and independence in population codes. J Neurosci 23: 11539–11553, 2003.
Shadlen MN, Newsome WT. The variable discharge of cortical neurons: implications for connectivity, computation, and information coding. J Neurosci 18: 3870–3896, 1998.
Steuber V, Mittmann W, Hoebeek FE, Silver RA, De Zeeuw CI, Hausser M, De Schutter E. Cerebellar LTD and pattern recognition by Purkinje cells. Neuron 54: 121–136, 2007.[CrossRef][Medline]
Victor JD, Purpura KP. Nature and precision of temporal coding in visual cortex: a metric-space analysis. J Neurophysiol 76: 1310–1326, 1996.
Wald A. Sequental Analysis. New York: Wiley, 1947.
Wichmann T, Kliem MA. Neuronal activity in the primate substantia nigra pars reticulata during the performance of simple and memory-guided elbow movements. J Neurophysiol 91: 815–827, 2004.
Wickens J. Basal ganglia: structure and computation. Network Comput Neural Syst 8: R77–R109, 1997.[CrossRef]
Zhang R, Tomida M, Katayama Y, Kawakami Y. Response durations encode nociceptive stimulus intensity in the rat medial prefrontal cortex. Neuroscience 125: 777–785, 2004.[CrossRef][Web of Science][Medline]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |