Information Theory enables the quantification of how much information a neuronal response carries about external stimuli and is hence a natural analytic framework for studying neural coding. The main difficulty in its practical application to spike train analysis is that estimates of neuronal information from experimental data are prone to a systematic error (called "bias"). This bias is an inevitable consequence of the limited number of stimulus-response samples that it is possible to record in a real experiment. In this paper, we first explain the origin and the implications of the bias problem in spike train analysis. We then review and evaluate some recent general-purpose methods to correct for sampling bias: the Panzeri-Treves, Quadratic Extrapolation, Best Universal Bound, Nemenman-Shafee-Bialek procedures, and a recently proposed shuffling bias reduction procedure. Finally, we make practical recommendations for the accurate computation of information from spike trains. Our main recommendation is to estimate information using the shuffling bias reduction procedure in combination with one of the other four general purpose bias reduction procedures mentioned in the preceding text. This provides information estimates with acceptable variance and which are unbiased even when the number of trials per stimulus is as small as the number of possible discrete neuronal responses.

Although the timing of single spikes is known to code for time-varying features of a sensory stimulus, it remains unclear whether time is also exploited in the neuronal coding of the spatial structure of the environment, where nontemporal stimulus features are fundamental. This report demonstrates that, in the whisker representation of rat cortex, precise spike timing of single neurons increases the information transmitted about stimulus location by 44%, compared to that transmitted only by the total number of spikes. Crucial to this code is the timing of the first spike after whisker movement. Complex, single neuron spike patterns play a smaller, synergistic role. Timing permits very few spikes to transmit high quantities of information about a behaviorally significant, spatial stimulus.

Neuronal responses to ongoing stimulation in many systems change over time, or “adapt.” Despite the ubiquity of adaptation, its effects on the stimulus information carried by neurons are often unknown. Here we examine how adaptation affects sensory coding in barrel cortex. We used spike-triggered covariance analysis of single-neuron responses to continuous, rapidly varying vibrissa motion stimuli, recorded in anesthetized rats. Changes in stimulus statistics induced spike rate adaptation over hundreds of milliseconds. Vibrissa motion encoding changed with adaptation as follows. In every neuron that showed rate adaptation, the input–output tuning function scaled with the changes in stimulus distribution, allowing the neurons to maintain the quantity of information conveyed about stimulus features. A single neuron that did not show rate adaptation also lacked input–output rescaling and did not maintain information across changes in stimulus statistics. Therefore, in barrel cortex, rate adaptation occurs on a slow timescale relative to the features driving spikes and is associated with gain rescaling matched to the stimulus distribution. Our results suggest that adaptation enhances tactile representations in primary somatosensory cortex, where they could directly influence perceptual decisions.

This study explores the nature of population coding in sensory cortex by applying information theoretic analyses to neuron pairs recorded simultaneously from rat barrel cortex. We quantified the roles of individual spikes and spike patterns in encoding whisker stimulus location. 82%-85% of the total information was contained in the timing of individual spikes: first spike time was particularly crucial. Spike patterns within neurons accounted for the remaining 15%-18%. Neuron pairs located in the same barrel column coded redundantly, whereas pairs in neighboring barrel columns coded independently. The barrel cortical population code for stimulus location appears to be the time of single neurons' first poststimulus spikes-a fast, robust coding mechanism that does not rely on "synergy" in crossneuronal spike patterns.

Rats, using their whiskers, have excellent capabilities in texture discrimination. What is the representation of texture in rat somatosensory cortex? We hypothesize that as rats "whisk" over a surface, the spatial frequency of a grooved or pebbled texture is converted to a temporal frequency of whisker vibration. Surface features such as groove depth or grain size modulate the amplitude of this vibration. Validation of the hypothesis depends on showing that vibration parameters have distinct neuronal representations in cortex. To test this, we delivered sinusoidal vibrations to the whisker shaft and analyzed cortical neuronal activity. Seven amplitudes and seven frequencies were combined to construct 49 stimuli while recording activity through a 10 x 10 microelectrode array inserted into the middle layers of barrel cortex. We find that cortical neurons do not explicitly encode vibration frequency (f) or amplitude (A) by any coding measure (average spike counts over different time windows, spike timing patterns in the peristimulus time histograms or in autocorrelograms). Instead, neurons explicitly encode the product of frequency and amplitude, which is proportional to the mean speed of the vibration. The quantity Af is an invariant because neuronal response encodes this feature independently of the values of the individual terms A and f. This was true across a wide time scale of firing rate measurements, from 5 to 500 msec. We conclude that vibration kinetics are rapidly and reliably encoded in the firing rate of cortical ensembles. Therefore, the cortical representation of vibration speed could underlie texture discrimination.